• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2000.11a
• /
• pp.74-77
• /
• 2000

Infiltration and water flow in the upper soil layer of a deep water table aquifer are modeled for multistorm runoff events. The infiltration process is developed using the sharp wetting front model of Green and Ampt, and the following redistribution process is modeled using the gravity drained rectangular approximation. The Brooks-Corey model [Brooks and Corey, 1966] is adopted to relate the effective soil saturation, the tension head, and the unsaturated hydraulic conductivity Firstly, the infiltration and redistribution model is developed for a single stom runoff event. Then a couple of events combined for multistorm runoff events. In the later case, infiltration rate of the second rainfall is strongly influenced by the length of the rainfall hiatus and soil moisture profile.

• Journal of Korean Society of Environmental Engineers
• /
• v.27 no.9
• /
• pp.995-1005
• /
• 2005

The study with laboratory sandbox model has been carried out to address potential use of reclaimed water, as a countermeasure artificially recharging the coastal aquifer, to effectively prevent from seawater intrusion due to overexploitation. It also investigated plausibility for either preserving or recovering the freshwater interface facing with seawater intrusion. To do this, we assessed hydraulic properties in artificial aquifer seawater/freshwater interface) depending upon the variation of extraction, storage and injection of reclaimed water. The variation of interface between freshwater and seawater were visualized by Surfer 8(Golden Software, USA) according to given experimental conditions. The interface between seawater and freshwater has been sensitively influenced by the change of extraction rate, where seawater zone migrated much faster into freshwater zone even though extraction rate became decreased. However, decreasing recharge rate could slow down moving of saline water zone toward freshwater zone. When the recharge was solely introduced into the sand box model, saline water intrusion was retarded than those of recharge and extraction working together. And also, the level of salinity of saline water was diluted by artificial recharge. It finally revealed that the artificial recharge would hydraulically avoid seawater intrusion while the freshwater sources could be conservatively utilized.

• The Journal of Engineering Geology
• /
• v.9 no.2
• /
• pp.147-159
• /
• 1999

Geostatistical methods were used for the groundwater flow analysis in a heterogeneous anisotropic aquifer. This study area is located at Sonbul-myeon in Hampyong-gun of Cheonnam Province which is a hydrogeological project area of KORES(Korea Resources Cooperation). Linear regression analysis shows that the topographic elevation and groundwater level of this area have very high correlation. Groundwater-level contour maps produced by ordinary kriging and cokringing have large differences in mountain areas, but small differences in hill and plain areas near the West Sea. Comparing two maps on the basis of an elevation contour map, a groundwater-level contour map using cokriging is more accurate. Analyzing the groundwater flow on two groundwater-level contour maps, the groundwater of study area flows from the high mountain areas to the plain areas near the West Sea. To verify the enffectiveness of geostatistical methods for the groundwater flow analysis in a heterogeneous anisotropic aquifer, the flow directions of groundwater were measured at two groundwater boreholes by a groundwater flowmeter system(model 200 $GeoFlo^{R}$). The measured flow directions of groundwater almost accord with those estimated on two groundwater-level contour maps produced by geostatistical methods.

• Journal of Korea Water Resources Association
• /
• v.42 no.11
• /
• pp.1005-1015
• /
• 2009

Saltwater intrusion in coastal aquifer was investigated using a laboratory model. Salt-wedge profiles were reproduced in a porous media tank 140 cm long, 70 cm high, and 10 cm wide. The experiments were performed with various conditions of porous media hydraulic conductivity, salinity, and ground surface slope to assess relationships on salt wedge location and inclination. Salt-wedge profiles induced by saltwater intrusion were observed in porous media equilibrium state, and compared with previously derived formulas of the Glover (1959), Henry (1959) and Strack (1976). It was found that salt-wedge shape and formations were affected by the water level ratio ($H_F/H_S$) due to high hydraulic conductivity, saltwater salinity and ground surface slope. High $H_F/H_S$ of porous media having high hydraulic conductivity shifted the saltwater interface toward the saltwater reservoir. Increasing surface slope of the porous media caused the salt-wedge profile inclination to decrease. Saltwater salinity also contributed to the location of saltwater interface, yet the impact was not more significant than hydraulic conductivity.

• Journal of Korea Water Resources Association
• /
• v.53 no.11
• /
• pp.985-997
• /
• 2020

SWAT (Soil and Water Assessment Tool)-MODFLOW (Modular Groundwater Flow) is a coupled model that linking semi-distributed watershed hydrology with fully-distributed groundwater behavior. In this study, the groundwater simulation results of SWAT and SWAT-MODFLOW were compared for Bokhacheon watershed in Namhan river basin. The models were calibrated and validated with 9 years (2009~2017) daily streamflow (Q) data of Heungcheon (HC) water level gauge station and the daily groundwater level observation data of Yulheon (YH). For SWAT, the groundwater parameters of GW_DELAY, GWQMN, and ALPHA_BF affecting baseflow and recession phase were treated. The SWAT results showed the coefficient of determination (R²) of 0.7 and Nash-Sutcliffe model efficiencies (NES_Q, NSE_inQ) for Q and 1/Q with 0.73 and -0.1 respectively. For SWAT-MODFLOW, the spatio-temporal aquifer hydraulic conductivity (K, m/day), specific storage (S_s, 1/m), and specific yield (S_y) were applied. The SWAT-MODFLOW showed R², NSE_Q, and NSE_inQ of 0.69, 0.74, and 0.51 respectively. The SWAT-MODFLOW considerably enhanced the low flow simulation with the help of aquifer physical information. The total streamflow of SWAT and SWAT-MODFLOW were 718.6 mm and 854.9 mm occupying baseflow of 342.9 mm and 423.5 mm respectively.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.227-227
• /
• 2015

The Wairarapa Valley occupies a predominantly rural area in the lower North Island of New Zealand. It supports a mix of intensive farming (dairy), dry stock farming (sheep and beef cattle) and horticulture (including wine grapes). The valley floor is traversed by the Ruamahanga River, the largest river in the Wellington region with a total catchment area of 3,430 km2. Environmental, cultural and recreational values associated with this Ruamahanga River are very high. The alluvial gravel and sand aquifers of the Wairarapa Valley, support productive groundwater aquifers at depths of up to 100 metres below ground while the Ruamahanga River and its tributaries present a further source of water for users. Water is allocated to users via resource consents by Greater Wellington Regional Council (GWRC). With intensifying land use, demand from the surface and groundwater resources of the Wairarapa Valley has increased substantially in recent times and careful management is needed to ensure values are maintained. This paper describes the approach being taken to manage water resources in the Wairarapa Valley and redefine appropriate limits of sustainable water use. There are three key parts: Quantifying the groundwater resource. A FEFLOW numerical groundwater flow model was developed by GWRC. This modelling phase provided a much improved understanding of aquifer recharge and abstraction processes. It also began to reveal the extent of hydraulic connection between aquifer and river systems and the importance of moving towards an integrated (conjunctive) approach to allocating water. Development of a conjunctive management framework. The FEFLOW model was used to quantify the stream flow depletion impacts of a range of groundwater abstraction scenarios. From this, three abstraction categories (A, B and C) that describe diminishing degrees of hydraulic connection between ground and surface water resources were mapped in 3 dimensions across the Valley. Interim allocation limits have been defined for each of 17 discrete management units within the valley based on both local scale aquifer recharge and stream flow depletion criteria but also cumulative impacts at the valley-wide scale. These allocation limits are to be further refined into agreed final limits through a community-led decision making process. Community involvement in the limit setting process. Historically in New Zealand, limits for sustainable resource use have been established primarily on the basis of 'hard science' and the decision making process has been driven by regional councils. Community involvement in limit setting processes has been through consultation rather than active participation. Recent legislation in the form of a National Policy Statement on Freshwater Management (2011) is reforming this approach. In particular, collaborative consensus-based decision making with active engagement from stakeholders is now expected. With this in mind, a committee of Wairarapa local people with a wide range of backgrounds was established in 2014. The role of this committee is to make final recommendations about resource use limits (including allocation of water) that reflect the aspirations of the communities they represent. To assist the committee in taking a holistic view it is intended that the existing numerical groundwater flow models will be coupled with with surface flow, contaminant transport, biological and economic models. This will provide the basis for assessing the likely outcomes of a range of future land use and resource limit scenarios.

• Journal of Environmental Science International
• /
• v.28 no.9
• /
• pp.737-749
• /
• 2019

In this study, the sensitivity analysis of hydraulic conductivity and separation distance (distance between injection well and pumping well) was analyzed by establishing a conceptual model considering the hydrogeologic characteristics of facility agricultural complex in Korea. In the conceptual model, natural characteristics (topography and geology, precipitation, hydraulic conductivity, etc.) and artificial characteristics (separation distance from injection well to pumping well, injection rate and pumping rate, etc.) is entered, and sensitivity analysis was performed 12 scenarios using a combination of hydraulic conductivity ($10^{-1}cm/sec$, $10^{-2}cm/sec$, $10^{-3}cm/sec$, $10^{-4}cm/sec$) and separation distance (10 m, 50 m, 100 m). Groundwater drawdown at the monitoring well was increased as the hydraulic conductivity decreased and the separation distance increased. From the regression analysis of groundwater drawdown as a hydraulic conductivity at the same separation distance, it was found that the groundwater level fluctuation of artificial recharge aquifer was dominantly influenced by hydraulic conductivity. In the condition that the hydraulic conductivity of artificial recharge aquifer was $10^{-2}cm/sec$ or more, the radius of influence of groundwater level was within 20 m, but In the condition that the hydraulic conductivity is $10^{-3}cm/sec$ or less, it is confirmed that the radius of influence of groundwater increases sharply as the separation distance increases.

• Journal of Soil and Groundwater Environment
• /
• v.13 no.4
• /
• pp.8-21
• /
• 2008

A series of three-dimensional numerical simulations using a hydrodynamic dispersion numerical model is performed to analyze quantitatively impacts of layered heterogeneity of geologic media and groundwater pumping schemes on groundwater flow and salt transport in coastal aquifer systems. A two-layer heterogeneous coastal aquifer system composed of a lower sand layer (aquifer) and an upper clay layer (aquitard) and a corresponding single-layer homogeneous coastal aquifer system composed of an equivalent lumped material are simulated to evaluate impacts of layered heterogeneity on seawater intrusion. In addition, a continuous groundwater pumping scheme and two different periodical groundwater pumping schemes, which withdraw the same amount of groundwater during the total simulation time, are applied to the above two coastal aquifer systems to evaluate impacts of groundwater pumping schemes on seawater intrusion. The results of the numerical simulations show that the periodical groundwater pumping schemes have more significant adverse influences on groundwater flow and salt transport not only in the lower sand layer but also in the upper clay layer, and groundwater salinization becomes more intensified spatially and temporally as the pumping intensity is higher under the periodical groundwater pumping schemes. These imply that the continuous groundwater pumping scheme may be more suitable to minimize groundwater salinization due to seawater intrusion. The results of the numerical simulations also show that groundwater salinization in the upper clay layer occurs significantly different from that in the lower sand layer under the periodical groundwater pumping schemes. Such differences in groundwater salinization between the two adjacent layers may result from layered heterogeneity of the layered coastal aquifer system.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2023.05a
• /
• pp.287-287
• /
• 2023

Surface water-groundwater interaction (SWGI) is an important hydrological process that influences both the quantity and quality of water resources. However, regional scale SWGI model calibration and uncertainty analysis have been a challenge because integrated models inherently carry a vast number of parameters, modeling assumptions, and inputs, potentially leaving little time and budget to explore questions related to model performance and forecasting. In this study, we have proposed the application of iterative ensemble smoother (IES) for uncertainty analysis and calibration of the widely used integrated surface-subsurface model, SWAT-MODFLOW. SWAT-MODFLOW integrates Soil and Water Assessment Tool (SWAT) and a three-dimensional finite difference model (MODFLOW). The model was calibrated using a parameter estimation tool (PEST). The major advantage of the employed IES is that the number of model runs required for the calibration of an ensemble is independent of the number of adjustable parameters. The pilot point approach was followed to calibrate the aquifer parameters, namely hydraulic conductivity, specific storage, and specific yield. The parameter estimation process for the SWAT model focused primarily on surface-related parameters. The uncertainties both in the streamflow and groundwater level were assessed. The work presented provides valuable insights for future endeavors in coupled surface-subsurface modeling, data collection, model development, and informed decision-making.

• Journal of Soil and Groundwater Environment
• /
• v.20 no.3
• /
• pp.51-64
• /
• 2015

In the present study, we develop two history matching techniques based on Markov chain Monte Carlo method where radial basis function and Gaussian distribution generated by unconditional geostatistical simulation are employed as the random walk transition kernels. The Bayesian inverse methods for aquifer characterization as the developed models can be effectively applied to the condition even when the targeted information such as hydraulic conductivity is absent and there are transient hydraulic head records due to imposed stress at observation wells. The model which uses unconditional simulation as random walk transition kernel has advantage in that spatial statistics can be directly associated with the predictions. The model using radial basis function network shares the same advantages as the model with unconditional simulation, yet the radial basis function network based the model does not require external geostatistical techniques. Also, by employing radial basis function as transition kernel, multi-scale nested structures can be rigorously addressed. In the validations of the developed models, the overall predictabilities of both models are sound by showing high correlation coefficient between the reference and the predicted. In terms of the model performance, the model with radial basis function network has higher error reduction rate and computational efficiency than with unconditional geostatistical simulation.