• Journal of the Korean Geosynthetics Society
• /
• v.16 no.2
• /
• pp.79-87
• /
• 2017

In this paper, for 43 sites neighboring to western area of Gangwondo where disaster of debris flow occurred from 2006 to 2013, magnitude of debris flow was estimated from results of site investigation and correlation analysis between influencing factors to its magnitude was performed. Magnitude of channelized debris flow was found greater by 6.5 times of that of hill slope debris flow and approximately 5% of total volume was occurred at initiation part of channelized debris flow. As results of analyzing yield rate of debris flow, for channelized debris flow, yield rate values of $19m^3/m$ and $8m^3/m$ were obtained for total volume being over $10,000m^3/m$ as the large scale of debris flow and less than $10,000m^3/m$ respectively, and value of $5m^3/m$ was estimated for hill slope debris flow. As results of correlation analysis of influencing factors to magnitude of debris flow, runoff distance and erosion width were very highly correlated to its magnitude whereas average slope of basin and erosion depth showed relatively low correlation. In particular, value of erosion depth was in the range of 0.5-2.6 m, being similar range to the value proposed by Ikeya (1981). Triggering rainfall to debris flow such as continuous rainfall and maximum intensity of hour rainfall were analyzed to have low correlation with magnitude of debris flow.

• Journal of Soil and Groundwater Environment
• /
• v.26 no.6
• /
• pp.106-117
• /
• 2021

The purpose of this study is prediction of landslide occurrence reflecting the subsurface flow characteristics within the soil layer in the future due to climate change in a large scale watershed. To do this, we considered the infinite slope stability theory to evaluate the landslide occurrence with predicted soil moisture content by SWAT model based on monitored data (rainfall-soil moisture-discharge). The correlation between the SWAT model and the monitoring data was performed using the coefficient of determination (R²) and the model's efficiency index (Nash and Sutcliffe model efficiency; NSE) and, an accuracy analysis of landslide prediction was performed using auROC (area under Receiver Operating Curve) analysis. In results comparing with the calculated discharge-soil moisture content by SWAT model vs. actual observation data, R² was 0.9 and NSE was 0.91 in discharge and, R² was 0.7 and NSE was 0.79 in soil moisture, respectively. As a result of performing infinite slope stability analysis in the area where landslides occurred in the past based on simulated data (SWAT analysis result of 0.7~0.8), AuROC showed 0.98, indicating that the suggested prediction method was resonable. Based on this, as a result of predicting the characteristics of landslide occurrence by 2050 using climate change scenario (RCP 8.5) data, it was calculated that four landslides could occur with a soil moisture content of more than 75% and rainfall over 250 mm/day during simulation. Although this study needs to be evaluated in various regions because of a case study, it was possible to determine the possibility of prediction through modeling of subsurface flow mechanism, one of the most important attributes in landslide occurrence.

• Journal of Korea Water Resources Association
• /
• v.43 no.9
• /
• pp.801-811
• /
• 2010

A three-dimensional numerical model called NEWTANK is employed to investigate solitary wave run-up with an internal wave-maker on a steep slope. The numerical model solves the spatially averaged Navier-Stokes equations for two-phase flows. The LES (large-eddy-simulation) approach is adopted to model the turbulence effect by using the Smagorinsky SGS (sub-grid scale) closure model. A two-step projection method is adopted in numerical solutions, aided by the Bi-CGSTAB (Bi-Conjugate Gradient Stabilized) method to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate VOF (volume-of-fluid) method is used to track the distorted and broken free surface. A solitary wave is first internally generated and propagated over a constant water depth in the three-dimensional domain. Numerically predicted results are compared with analytical solutions and numerical errors are analyzed in detail. The model is then applied to study solitary wave run-up on a steep slope and the obtained results are compared with available laboratory measurements.

• Journal of Korean Society of Forest Science
• /
• v.83 no.2
• /
• pp.175-190
• /
• 1994

In order to get fundamental information for prediction of landslide hazard, both forest and site factors affecting slope stability were investigated in many areas of active landslides. Twelve descriptors were identified and quantified to develop the prediction model by multivariate statistical analysis. The main results obtained could be summarized as follows : The main factors influencing a large scale of landslide were shown in order of precipitation, age group of forest trees, altitude, soil texture, slope gradient, position of slope, vegetation, stream order, vertical slope, bed rock, soil depth and aspect. According to partial correlation coefficient, it was shown in order of age group of forest trees, precipitation, soil texture, bed rock, slope gradient, position of slope, altitude, vertical slope, stream order, vegetation, soil depth and aspect. The main factors influencing a landslide occurrence were shown in order of age group of forest trees, altitude, soil texture, slope gradient, precipitation, vertical slope, stream order, bed rock and soil depth. Two prediction models were developed by magnitude and frequency of landslide. Particularly, a prediction method by magnitude of landslide was changed the score for the convenience of use. If the total store of the various factors mark over 9.1636, it is evaluated as a very dangerous area. The mean score of landslide and non-landslide group was 0.1977 and -0.1977, and variance was 0.1100 and 0.1250, respectively. The boundary value between the two groups related to slope stability was -0.02, and its predicted rate of discrimination was 73%. In the score range of the degree of landslide hazard based on the boundary value of discrimination, class A was 0.3132 over, class B was 0.3132 to -0.1050, class C was -0.1050 to -0.4196, class D was -0.4195 below. The rank of landslide hazard could be divided into classes A, B, C and D by the boundary value. In the number of slope, class A was 68, class B was 115, class C was 65, and class D was 52. The rate of landslide occurrence in class A and class B was shown at the hige prediction of 83%. Therefore, dangerous areas selected by the prediction method of landslide could be mapped for land-use planning and criterion of disaster district. And also, it could be applied to an administration index for disaster prevention.

• Journal of The Korean Astronomical Society
• /
• v.46 no.3
• /
• pp.125-131
• /
• 2013

We present the relation between the genus in cosmology and the Betti numbers for excursion sets of three- and two-dimensional smooth Gaussian random fields, and numerically investigate the Betti numbers as a function of threshold level. Betti numbers are topological invariants of figures that can be used to distinguish topological spaces. In the case of the excursion sets of a three-dimensional field there are three possibly non-zero Betti numbers; ${\beta}_0$ is the number of connected regions, ${\beta}_1$ is the number of circular holes (i.e., complement of solid tori), and ${\beta}_2$ is the number of three-dimensional voids (i.e., complement of three-dimensional excursion regions). Their sum with alternating signs is the genus of the surface of excursion regions. It is found that each Betti number has a dominant contribution to the genus in a specific threshold range. ${\beta}_0$ dominates the high-threshold part of the genus curve measuring the abundance of high density regions (clusters). ${\beta}_1$ dominates the genus near the median thresholds which measures the topology of negatively curved iso-density surfaces, and ${\beta}_2$ corresponds to the low-threshold part measuring the void abundance. We average the Betti number curves (the Betti numbers as a function of the threshold level) over many realizations of Gaussian fields and find that both the amplitude and shape of the Betti number curves depend on the slope of the power spectrum n in such a way that their shape becomes broader and their amplitude drops less steeply than the genus as n decreases. This behaviour contrasts with the fact that the shape of the genus curve is fixed for all Gaussian fields regardless of the power spectrum. Even though the Gaussian Betti number curves should be calculated for each given power spectrum, we propose to use the Betti numbers for better specification of the topology of large scale structures in the universe.

• Journal of Korea Water Resources Association
• /
• v.42 no.9
• /
• pp.715-724
• /
• 2009

In case of inclined crest groyne, flow constriction with water is different. Therefore, it is proper to apply to narrow channel or there is a risk of overflow with water level rising caused by installation of groynes. This study were conducted experiments on inclined crest groyne. Main factors of inclined crest groyne are angles of crests slope and changes of water level. Velocity profile around groyne was measured by LSPIV (Large Scale Particle Image Velocimetry). Flow fields around groyne were analyzed focused on main channel and recirculation area. From the results, Thalweg change has little an effect on rate length and height of recirculation rises on increasing rate length. Length of recirculation area was about 12$\sim$16 times of rate length of groyne. Maximum velocity in main channel area was measured about 1.45$\sim$2.1 times of approach velocity and has little an effect on crest angle of groyne. Back water velocity recirculation area was decreased on approach velocity. This result presents to make stable flow to bank protection.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.1
• /
• pp.81-96
• /
• 2007

Nonlinear shoaling characteristics over surf zone are numerically investigated based on spatially averaged NavierStokes equation. We also test the validity of gradient model for turbulent stresses due to wave breaking using the data acquainted during SUPERTANK LABORATORY DATA COLLECTION PROJECT(Krauss et al., 1992). It turns out that the characteristics length scale of breaking induced current is not negligible, which firmly stands against ever popular gradient model, ${\kappa}-{\varepsilon}$ model, but favors Large Eddy Simulation with finer grid. Based on these observations, we model the residual stress of spatially averaged NavierStokes equation after Lagrangian Dynamic Smagorinsky(Meneveau et al., 1996). We numerically integrate newly proposed wave equations using SPH with Gaussian kernel function. Severely deformed water surface profile, free falling water particle, queuing splash after landing of water particle on the free surface and wave finger due to structured vortex on rear side of wave crest(Narayanaswamy and Dalrymple, 2002) are successfully duplicated in the numerical simulation of wave propagation over uniform slope beach, which so far have been regarded very difficult features to mimic in the computational fluid mechanics.

• Journal of Korean Society of Forest Science
• /
• v.107 no.4
• /
• pp.393-401
• /
• 2018

The aim of this study was to analyze the characteristics of ground subsidence and land-creeping in four mining areas in Korea. Land creeping were occurred by large-scale sinkholes and tensile cracks in the mining areas. This ares showed the shape as large-scale landslides by the collapsed stone debris pushed in the direction of gravity. The rock type in the No. 4 survey area was sandstone and shale, whereas that of other three areas was limestone. The limestone was generally vulnerable to land creeping. The ground subsidence and land creeping in four mining areas were located near the ridge of the mountain. The land creeping by ground subsidence in the mining areas showed complex shapes due to the cavity of underground tunnel and the colluvial debris by gravity in the mountain area, whereas the land creeping in Korea are in the shape of landslide. The average slope of the mining area was ranged from $30.4^{\circ}$ to $33.7^{\circ}$. The values were higher ($1.5^{\circ}{\sim}4.8^{\circ}$) than the average slopes for landslide areas.

• Economic and Environmental Geology
• /
• v.55 no.5
• /
• pp.497-510
• /
• 2022

Monthly monitoring, 3D scan survey, and electrical resistivity survey were conducted from January 2018 to August 2022 to identify the cause of rockfall occurring in Daljeon-ri Columnar Joint (Natural Monument No. 415), Pohang. A total of 3,231 rocks fell from the columnar joint over the past 5 years, and 1,521 (47%) of the falling rocks were below 20 cm in length, 978 (30.3%) of 20-30 cm, and 732 (22.7%) of rocks over 30 cm. While the number of rockfalls by year has decreased since 2018, the frequency of rockfalls bigger than 30 cm tends to increase. Large-scale rockfalls occurred mainly during the thawing season (March-April) and the rainy season (June-July), and the analysis of the relationship between cumulative rainfall and rockfall occurrence showed that cumulative rainfall for 3 to 4 days is also closely related to the occurrence of rockfall. Smectite and illite, which are expansible clay minerals, were observed in XRD analysis of the slope material (filling minerals) in the columnar joint, and the presence of a fault fracture zone was confirmed in the electrical resistivity survey. In addition, the confirmed fault fracture zone and the maximum erosion point analyzed through 3D precision measurement coincided with the main rockfall occurrence point observed by the BTC-6PXD camera. Therefore, the main cause of rockfall at Daljeon-ri columnar joint in Pohang is a combination of internal factors (development of fault fracture zones and joints, weathering of rocks, presence of expansive clay minerals) and external factors (precipitation, rapid thawing phenomenon), resulting in large-scale rockfall. Meanwhile, it was also confirmed that the Pohang-Gyeongju earthquake, which was continuously raised, was not the main cause.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.41 no.6
• /
• pp.727-736
• /
• 2021

Road tunnel lengths are increasing. Some 1,300 tunnels with 1,102 km in length had been increased till 2019 from 2010. There are 64 tunnels over 3,000 m in length, with their total length adding up to 276.7 km. Safety facilities in the event of a tunnel fire are critical so as to prevent large-scale casualties. Standards for installing disaster prevention facilities are being proposed based on the guidelines of the Ministry of Land, Infrastructure and Transport, but they may be limited to deep underground tunnels. This study was undertaken to provide guidelines for the spacing of evacuation connection passages and the widths of evacuation connection doors. Evacuation with various spacing and widths was simulated in regards to evacuation time, which is the measure of safety, using the evacuation analysis simulation software EXODUS Ver.6.3 and the fire/smoke analysis software SMARTFIRE Ver.4.1. Evacuation connection gates with widths of 0.9 m and 1.2 m, and spacings of 150 m to 250 m, were set to every 20 m. In addition, longitudinal slopes of 6 % and 0 % were considered. It was determined to be safe when the evacuation completion time was shorter than the delay diffusion time. According to the simulation results, all occupants could complete evacuation before smoke spread regardless of the width of the evacuation connection door when the longitudinal slope was 6 % and the interval of evacuation connection passage was 150 m. When the evacuation connection passage spacing was 200 m and the evacuation connection gate width was 1.2 m, all occupants could evacuate when the longitudinal slope was 0 %. Due to difference in evacuation speed according to the longitudinal slope, the evacuation time with a 6 % slope was 114 seconds shorter (with the 190 m connection passage) than with a 0 % slope. A shorter spacing of evacuation connection passages may reduce the evacuation time, but this is difficult to implement in practice because of economic and structural limitations. If the width of the evacuation junction is 1.2 m, occupants could evacuate faster than with a 0.9 m width. When the width of a connection door is 1.2 m with appropriate connection passage spacing, it might provide a means to increase economic efficiency and resolve structural limitations while securing evacuation safety.