• Journal of Soil and Groundwater Environment
• /
• v.22 no.6
• /
• pp.66-73
• /
• 2017

Empirical erosion models like Universal Soil Loss Equation (USLE) models have been widely used to make spatially distributed soil erosion vulnerability maps. Even if the models detect vulnerable sites relatively well utilizing big data related to climate, geography, geology, land use, etc within study domains, they do not adequately describe the physical process of soil erosion on the ground surface caused by rainfall or overland flow. In other words, such models are still powerful tools to distinguish the erosion-prone areas at large scale, but physics-based models are necessary to better analyze soil erosion and deposition as well as the eroded particle transport. In this study a physics-based soil erosion modeling system was developed to produce both runoff and sediment yield time series at watershed scale and reflect them in the erosion and deposition maps. The developed modeling system consists of 3 sub-systems: rainfall pre-processor, geography pre-processor, and main modeling processor. For modeling system validation, we applied the system for various erosion cases, in particular, rainfall-runoff-sediment yield simulation and estimation of probable maximum sediment (PMS) correlated with probable maximum rainfall (PMP). The system provided acceptable performances of both applications.

• Economic and Environmental Geology
• /
• v.51 no.2
• /
• pp.121-130
• /
• 2018

The KURT (KAERI Underground Research Tunnel) is a research tunnel which is located in KAERI (Korea Atomic Energy Research Institute) site. At KURT, researches on engineering and natural barrier system, which are the most important components for geological disposal system for high level radioactive waste, have been conducted. In this study, we synthesized the site characteristics obtained by various types of site investigation to introduce the geological model for KURT site, and induced the 3-D hydrogeological model for KURT site from the geological model. From the geological investigation at the surface and boreholes, four geological elements such as subsurface weathered zone, upper fractured rock, lower fractured rock and fracture zones were determined for the geological model. In addition, the geometries of these geological elements were also analyzed for the geological model to be three-dimensional. The results from 3-D geological model were used to construct the hydro-geological model for KURT site, which is one of the input data for groundwater flow modeling and safety assessment.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.275-284
• /
• 2007

For precise interpretation of magnetotelluric (MT) data distorted by irregular surface terrain, topography effects are investigated by computing apparent resistivities, phases, tippers and induction vectors for a three-dimensional (3D) hill-and-valley model. To compute MT responses for the 3D surface topography model, we use a 3D MT modeling algorithm based on an edge finite-element method which is free from vector parasites. Distortions on the apparent resistivity and phase are mainly caused by distorted currents that flow along surface topography. The distribution of tipper amplitudes over both hill and valley are the same, while the tipper points toward the center of hill and the base of the valley. The real part of induction vector also points in the same direction as that of tipper, while the imaginary part in the opposite direction.

• Nuclear Engineering and Technology
• /
• v.20 no.2
• /
• pp.132-136
• /
• 1988

A strategic question in groundwater contaminant transport modeling is whether we need to treat waste packages or drums as individual, discrete sources or as approximately lumped sources. In this paper we present analyses of array sources in porous media. We analyze a planar array of sources in porous media with groundwater flow. We compare the concentration field predicted by a detailed model of individual point sources to concentration fields predicted by an infinite plane source and a single point source, all of the same equivalent strength. From this study we identified three regions: (1) a region close to the sources where the effects of adjacent sources are significant and individual source models should be used, (2) a region extending from a few meters to hundreds to thousands of meters downstream, where an equivalent source of infinite extent gives accurate results, and (3) a far-field region, where in an equivalent source of finite extent gives accurate results.

• Journal of Korea Soil Environment Society
• /
• v.1 no.1
• /
• pp.89-102
• /
• 1996

An integrated model is presented to describe underground flow and mass transport, using a multicomponent multiphase approach. The comprehensive governing equation is derived considering mass and force balances of chemical species over four phases(water, oil, air, and soil) in a schematic elementary volume. Compact and systemati notations of relevant variables and equations are introduced to facilitate the inclusion of complex migration and transformation processes, and variable spatial dimensions. The resulting nonlinear system is solved by a multidimensional finite element code. The developed code with dynamic array allocation, is sufficiently flexible to work across a wide spectrum of computers, including an IBM ES 9000/900 vector facility, SP2 cluster machine, Unix workstations and PCs, for one-, two and three-dimensional problems. To reduce the computation time and storage requirements, the system equations are decoupled and solved using a banded global matrix solver, with the vector and parallel processing on the IBM 9000. To avoide the numerical oscillations of the nonlinear problems in the case of convective dominant transport, the techniques of upstream weighting, mass lumping, and elementary-wise parameter evaluation are applied. The instability and convergence criteria of the nonlinear problems are studied for the one-dimensional analogue of FEM and FDM. Modeling capacity is presented in the simulation of three dimensional composite multiphase TCE migration. Comprehesive simulation feature of the code is presented in a companion paper of this issue for the specific groundwater or flow and contamination problems.

• Tunnel and Underground Space
• /
• v.11 no.2
• /
• pp.156-166
• /
• 2001

A stochastic continuum(SC) modeling technique was developed to simulate the groundwater flow pathway in fractured rocks. This model was developed to overcome the disadvantageous points of discrete fracture network(DFN) modes which has the limitation of fracture numbers. Besides, SC model is able to perform probabilistic analysis and to simulate the conductive groundwater pathway as discrete fracture network model. The SC model was formulated based on the discrete fracture network(DFN) model. The spatial distribution of permeability in the stochastic continuum model was defined by the probability distribution and variogram functions defined from the permeabilities of subdivided smaller blocks of the DFN model. The analysis of groundwater travel time was performed to show the consistency between DFN and SC models by the numerical experiment. It was found that the stochastic continuum modes was an appropriate way to provide the probability density distribution of groundwater velocity which is required for the probabilistic safety assessment of a radioactive waste disposal facility.

• Journal of Korea Water Resources Association
• /
• v.33 no.1
• /
• pp.61-72
• /
• 2000

The purpose of this study is to improve the method of evaluating groundwater recharge by using grid-based soil moisture routing technique. A model which predicts temporal variation and spatial distribution of soil moisture on a daily time step was developed. The model uses ASCII-formatted map data supported by the irregular gridded map of the GRASS(Geographic Resources Analysis Support System)-GIS and can generate daily and monthly spatial distribution map of surface runoff, soil moisture content, evapotranspiration within the watershed. The model was applied to Ipyunggyo watershed($75.6\;\textrm{km}^2$) located in the upstream of Bocheongchun watershed. Seven maps; DEM(Digital Elevation Mode]), stream, flow path, soil, land use, Thiessen network and free groundwater level, were used for input data. Predicted streamflows resulting from two years (l995, 1996) daily data were compared with the observed values at the watershed outlet. The results of temporal variations and spatial distributions of soil moisture are presented by using GRASS GIS. As a final result, the monthly predicted groundwater recharge was presented.sented.

• Journal of Korea Water Resources Association
• /
• v.31 no.2
• /
• pp.177-187
• /
• 1998

The spatially integrated model (SIM) which can evaluate temporal variation of groundwater quality is proposed in the stream-aquifer setting entered by nonpoint source contaminants. And the developed SIM included unsaturated soil zone and was tested against the spatially distributed model (SDM) of the coupled advection-dispersion and Richards equations for the various hydrologic and aquifer simulating conditions. The result of the comparison showed that the average concentration responses of saturated aquifer and groundwater outflow between the SIM and the SDM was in good agreement, except for the case of the large dispersivity ratio and thick aquifer system. And it is shown that for the cases of the large dispersivity ratio and thick aquifer system the performance of the nonlinear SIM is better than that of the linear SIM for evaluating the average concentration of groundwater outflow response.

• Journal of Soil and Groundwater Environment
• /
• v.13 no.1
• /
• pp.77-91
• /
• 2008

In this study, we developed a technique of applying DRASTIC, which is the most widely used tool for estimation of groundwater vulnerability to the aqueous phase contaminant infiltrated from the surface, and a groundwater flow model jointly to assess groundwater contamination potential. The developed technique is then applied to Buyeo-eup area in Buyeo-gun, Chungcheongnam-do, Korea. The input thematic data of a depth to water required in DRASTIC model is known to be the most sensitive to the output while only a few observations at a few time schedules are generally available. To overcome this practical shortcoming, both steady-state and transient groundwater level distributions are simulated using a finite difference numerical model, MODFLOW. In the application for the assessment of groundwater vulnerability, it is found that the vulnerability results from the numerical simulation of a groundwater level is much more practical compared to cokriging methods. Those advantages are, first, the results from the simulation enable a practitioner to see the temporally comprehensive vulnerabilities. The second merit of the technique is that the method considers wide variety of engaging data such as field-observed hydrogeologic parameters as well as geographic relief. The depth to water generated through geostatistical methods in the conventional method is unable to incorporate temporally variable data, that is, the seasonal variation of a recharge rate. As a result, we found that the vulnerability out of both the geostatistical method and the steady-state groundwater flow simulation are in similar patterns. By applying the transient simulation results to DRASTIC model, we also found that the vulnerability shows sharp seasonal variation due to the change of groundwater recharge. The change of the vulnerability is found to be most peculiar during summer with the highest recharge rate and winter with the lowest. Our research indicates that numerical modeling can be a useful tool for temporal as well as spatial interpolation of the depth to water when the number of the observed data is inadequate for the vulnerability assessments through the conventional techniques.

• Journal of the Korean Society of Groundwater Environment
• /
• v.6 no.2
• /
• pp.101-106
• /
• 1999

Open dump of refuse causes groundwater and soil contamination by leachate and air pollution by LFG(Landfill Gas). In this paper, in order to perform a study about reduction of high leachate and LFG collection & control, using a 3-D, 2-phase, transient FDM model, the analysis of simultaneous flow of leachate and LFG has been carried out. In present numerical analysis it is assumed that 58 percents of LFG will evaporate to the ambient air and the recharge rate of a landfill be 12 percent of the average precipitation per year. All other data were excerpted at the point of 1995 when three refuse layers had been buried. From numerical analysis we concluded that maximum head value is approximately 26 mH2O<-에이치투오 (2.52 atm) in the center of the system and that installing venting trench plays an important role in landfill stabilization. Evan with the assumption of three layers constructed and low recharge rate applied, it is found that cumulative leachate and LFG productions will be 15.1 million 세제곱미터, 5.58 billion 세제곱미터, respectively after 40 years.