Groundwater Head Data Research Articles

Evolution of spatiotemporal ground deformation over 30 years in Xi’an, China, with multi-sensor SAR interferometry

Since the 1960s, Xi'an, China, has suffered from severe ground subsidence due to groundwater extraction and continuous urban construction. During the urbanization process, ground subsidence has threatened and will continue to threaten the stability of the city's infrastructure. Additionally, parts of Xi'an have recently experienced a period of rebound. In this context, we must monitor the spatial and long-term evolution of deformation in Xi'an city. First, archived SAR datasets from 1992 to the present, including ERS1/2, Envisat/ASAR, TerraSAR-X(TSX), ALSO/PALSAR-2(ALOS-2), and Sentinel-1(S1), provided independent historical deformation data in Xi’an city; different SAR interferometry methods can obtain the deformation rates and deformation time-series of the region. Second, to fuse the deformation time-series acquired by different sensors, the Quantile-Quantile Adjustment (QQA) algorithm was used to generate the long-term deformation time-series from October 2003 to May 2021, which has a higher precision than when using the Singular Value Decomposition (SVD) method. Quantitatively, the standard deviation (STD) of the QQA method was only 22 mm at maximum, considerably smaller than the 54 mm obtained with the SVD method. Third, ERS1/2 was linked by the segmental interpolation method to obtain the 30-year deformation time-series of Xi'an, which showed that ground deformation in Xi'an from 1992 to 2021 has experienced a “deceleration-acceleration-deceleration-rebound” process, with a maximum cumulative subsidence of over 2 m. Finally, land use type data and hydraulic head change data were also collected to analyze the deformation information. At the same time, the elastic and inelastic water storage coefficients of the aquifer were calculated by combining the long-term deformation time-series and groundwater head data; the deformation mechanism of the aquifer was analyzed accordingly.

Different Ground Subsidence Contributions Revealed by Integrated Discussion of Sentinel-1 Datasets, Well Discharge, Stratigraphical and Geomorphological Data: The Case of the Gioia Tauro Coastal Plain (Southern Italy)

Groundwater is the main water supply for agricultural and industrial needs in many coastal plains worldwide. Groundwater depletion often triggers land subsidence, which threatens manmade infrastructure and activities and aggravates other geohazards. We applied a multi-temporal interferometric synthetic aperture radar technique to Sentinel-1 datasets to detect ground motion in the Gioia Tauro plain (Calabria, Southern Italy) from 2018 to 2021. The InSAR data were analysed through the integrated use of groundwater head, stratigraphical and geomorphological data, and land use information to distinguish the potential subsidence divers. The results show that subsiding areas, with a mean rate of about 10 mm/yr, are in the middle of the plain, and their location is influenced by the spatial distribution of compressible sediments included in the shallow aquifer. Furthermore, the subsidence arrangement is spatially accordant with the main groundwater depression area, which can be ascribed to the ongoing and increasing water pumping for predominantly agricultural usage. We also observed that subsidence (up to 10 mm/yr) affects the western dock of the Gioia Tauro harbour, in front of which, in very shallow water, are two submarine canyon heads already affected by slides in the past.

A multi-compartment hydrologic model to estimate groundwater recharge in an alluvial-karst system

Estimation of groundwater recharge is essential for sustainable groundwater management planning in arid and semiarid areas. Accurate recharge estimation in a distributed scale is a great challenge, especially when multiple recharge sources exist and available data are limited and uncertain. There are various sophisticated direct and indirect methods such as water balance method to estimate this parameter. Conventional spatially and temporally lumped models cannot take into account spatial variation of hydraulic properties of an aquifer, land use, and hydrometeorological conditions in the water balance estimation. In this study, a semi-distributed hydrologic model is developed and applied in a complex alluvial-karst system in Firouzabad catchment (Iran) to estimate groundwater recharge from precipitation, subsurface inflow from adjacent karst aquifer, return flows, and riverbed infiltration. The proposed multi-compartment model considers groundwater balance in alluvial and karstic aquifers and surface soil layer. Furthermore, it contains a multi-cell aquifer model which is used to understand the macro mechanism of the alluvial aquifer. A parameter calibration procedure is implemented wherein groundwater head data are used in parameter estimation. The performance of the proposed model is assessed by evaluating the coefficient of determination (R2) and the coefficient of efficiency. The results of the groundwater balance model are analyzed to quantify the seasonal long-term (16 years, from 1992 to 2008) groundwater recharge components. The proposed model is straightforward, efficient, and easy to apply, and its results are reasonable and acceptable.

Groundwater Quality Impacts from a Full‐Scale Integrated Constructed Wetland

AbstractThe concept of integrated constructed wetlands (ICW) promotes in‐situ soils to construct and line wetland cells. The integrity of soil material, however, may provide a potential pathway for contaminants to flow into the underlying groundwater. This study assessed the extent of groundwater quality deterioration due to the establishment of a full‐scale ICW system treating domestic wastewater in Ireland. The ICW is located at Glaslough in Co. Monaghan, Ireland. It consists of two sedimentation ponds and a sequence of five shallow vegetated wetland cells. The ICW cells were lined with 500‐mm thick local subsoil material, which comprised a mixture of alluvium, organic soils, tills, and gravel. Groundwater samples and head data were collected from eight piezometers, which were installed around the ICW cells. The groundwater and wetland water samples were analysed for water quality parameters such as bulk organic matter, nutrients, and pathogens. Overall, the quality of groundwater underlying the ICW system recorded some contamination with bulk organic matter and some inorganic nutrients. Significantly higher contaminant concentrations were recorded in monitoring wells upgradient and near to the distal wetland cells than downgradient ones, which were near to the proximal cells. For the downgradient piezometers, concentrations seldomly exceeded the natural background levels. Detailed analyses through the application of chemometrics models indicated that the source of contamination was largely of geogenic origin. Findings suggest that ICW systems pose a minimal risk to the groundwater quality; the greatest risk was associated with the distal wetland cells.

A Hierarchical Bayesian Model Averaging Framework for Groundwater Prediction under Uncertainty

Groundwater prediction models are subjected to various sources of uncertainty. This study introduces a hierarchical Bayesian model averaging (HBMA) method to segregate and prioritize sources of uncertainty in a hierarchical structure and conduct BMA for concentration prediction. A BMA tree of models is developed to understand the impact of individual sources of uncertainty and uncertainty propagation to model predictions. HBMA evaluates the relative importance of different modeling propositions at each level in the BMA tree of model weights. The HBMA method is applied to chloride concentration prediction for the "1,500-foot" sand of the Baton Rouge area, Louisiana from 2005 to 2029. The groundwater head data from 1990 to 2004 is used for model calibration. Four sources of uncertainty are considered and resulted in 180 flow and transport models for concentration prediction. The results show that prediction variances of concentration from uncertain model elements are much higher than the prediction variance from uncertain model parameters. The HBMA method is able to quantify the contributions of individual sources of uncertainty to the total uncertainty.

Surface–groundwater interactions in hard rocks in Sardon Catchment of western Spain: An integrated modeling approach

Summary The structural and hydrological complexity of hard rock systems (HRSs) affects dynamics of surface–groundwater interactions. These complexities are not well described or understood by hydrogeologists because simplified analyses typically are used to study HRSs. A transient, integrated hydrologic model (IHM) GSFLOW (Groundwater and Surface water FLOW) was calibrated and post-audited using 18 years of daily groundwater head and stream discharge data to evaluate the surface–groundwater interactions in semi-arid, ∼80 km 2 granitic Sardon hilly catchment in Spain characterized by shallow water table conditions, relatively low storage, dense drainage networks and frequent, high intensity rainfall. The following hydrological observations for the Sardon Catchment, and more generally for HRSs were made: (i) significant bi-directional vertical flows occur between surface water and groundwater throughout the HRSs; (ii) relatively large groundwater recharge represents 16% of precipitation ( P , 562 mm . y −1 ) and large groundwater exfiltration (∼11% of P ) results in short groundwater flow paths due to a dense network of streams, low permeability and hilly topographic relief; deep, long groundwater flow paths constitute a smaller component of the water budget (∼1% of P ); quite high groundwater evapotranspiration (∼5% of P and ∼7% of total evapotranspiration); low permeability and shallow soils are the main reasons for relatively large components of Hortonian flow and interflow (15% and 11% of P , respectively); (iii) the majority of drainage from the catchment leaves as surface water; (iv) declining 18 years trend (4.44 mm . y −1 ) of groundwater storage; and (v) large spatio-temporal variability of water fluxes. This IHM study of HRSs provides greater understanding of these relatively unknown hydrologic systems that are widespread throughout the world and are important for water resources in many regions.

Conceptual modelling of groundwater–surface water interactions in the Lake Sibayi Catchment, Eastern South Africa

Lake Sibayi, a topographically closed fresh water lake in northern KwaZulu-Natal, South Africa, and coastal groundwater aquifers surrounding the lake, are important water resources and are used extensively for domestic water supplies. Both the lake and groundwater support important and ecologically sensitive wetland system in the area. Surface and subsurface geological information, groundwater head, hydrochemical and environmental isotope data were analysed to develop a conceptual model of aquifer–lake interaction, upon which a three dimensional numerical model will be based. Local geologic, groundwater head distribution, lake level, hydrochemical and environmental isotope data confirm direct hydraulic link between groundwater and the lake. In the western section of the catchment, groundwater flows to the lake where groundwater head is above lake stage, whereas along the eastern section, the presence of mixing between lake and groundwater isotopic compositions indicates that the lake recharges the aquifer. Stable isotope signals further revealed the movement of lake water through and below the coastal dune cordon before eventually discharging into the Indian Ocean. Quantification of the 15year monthly water balance for the lake shows strong seasonal variations of the water balance components and average residence time of water in the lake is about 7.5years. Recent increase in rate of water abstraction from the lake combined with decreasing precipitation and rapidly increasing pine plantations appears to be responsible for lake level reduction which would have dramatic negative effects on the neighbouring ecosystem and a potential seawater invasion of the coastal aquifer.

A New Zonation Algorithm with Parameter Estimation Using Hydraulic Head and Subsidence Observations

Parameter estimation codes such as UCODE_2005 are becoming well-known tools in groundwater modeling investigations. These programs estimate important parameter values such as transmissivity (T) and aquifer storage values (Sa ) from known observations of hydraulic head, flow, or other physical quantities. One drawback inherent in these codes is that the parameter zones must be specified by the user. However, such knowledge is often unknown even if a detailed hydrogeological description is available. To overcome this deficiency, we present a discrete adjoint algorithm for identifying suitable zonations from hydraulic head and subsidence measurements, which are highly sensitive to both elastic (Sske) and inelastic (Sskv) skeletal specific storage coefficients. With the advent of interferometric synthetic aperture radar (InSAR), distributed spatial and temporal subsidence measurements can be obtained. A synthetic conceptual model containing seven transmissivity zones, one aquifer storage zone and three interbed zones for elastic and inelastic storage coefficients were developed to simulate drawdown and subsidence in an aquifer interbedded with clay that exhibits delayed drainage. Simulated delayed land subsidence and groundwater head data are assumed to be the observed measurements, to which the discrete adjoint algorithm is called to create approximate spatial zonations of T, Sske , and Sskv . UCODE-2005 is then used to obtain the final optimal parameter values. Calibration results indicate that the estimated zonations calculated from the discrete adjoint algorithm closely approximate the true parameter zonations. This automation algorithm reduces the bias established by the initial distribution of zones and provides a robust parameter zonation distribution.

Evaluation of groundwater-surface water interaction through groundwater modelling: simulation of the effects of removal of a dam along a river at a contaminated site in Northern France

A numerical groundwater flow model has been developed for an industrial site bounded by a river in in Basse Normandie, Northern France. The scope of the work was the optimisation of the existing groundwater pump and treat system and the prediction of possible effects on groundwater circulation after the future removal of a dam located along the river. The model has been implemented with the finite difference code MODFLOW 2005 and represents an area with an extension of approximately 800 x 500 m. It has been calibrated using static conditions groundwater head data (wells deactivated) and verified with 1) abstracting conditions (wells abstracting) head data, 2) simulating pumping tests with transient simulations and 3) comparing measured average river baseflow with modelled river drainage. The model indicates that the hydraulic barrier in the present abstraction scenario has some problematic areas and needs some improvements, as confirmed by the hydrochemical data of the river water. A first predictive scenario has been developed to optimise the barrier, indicating that a flow rate of 0.5 m3/h each at three new barrier wells, in addition to the present abstraction scenario, should ensure the hydraulic containment of the site. A second predictive scenario simulates the optimised groundwater abstractions without the presence of the dam along the neighbouring river. In these conditions, the river will increase the drainage effect on the aquifer, requiring a further increase in the rate of abstraction from the existing and new wells to ensure the hydraulic containment. With this paper we would like to present an example of what we think is a correct professional approach, with the design of the simplest model as possible depending on the hydrogeological conceptual model complexity, the abundance of data and the model objectives, and where multiple confirmations of the correctness of groundwater model results have been searched for.

Groundwater-surface water interaction in Lake Nasser, Southern Egypt

A cross-sectional model, based on the two dimensional groundwater flow equation of Edelman, was applied at seven transects distributed over four geological cross sections to estimate groundwater heads and recharge from/or groundwater discharge to Lake Nasser. The lake with a length of 500 km and an average width of 12 km was created over the period 1964–1970, the time for constructing the Aswan High Dam (AHD). The model, constrained by regional-scale groundwater flow and groundwater head data in the vicinity of the lake, was successfully calibrated to timeseries of piezometeric heads collected at the cross sections in the period 1965–2004. Inverse modeling yielded high values for the horizontal hydraulic conductivity in the range of 6.0 to 31.1 m day−1 and storage coefficient between 0.01 and 0.40. The results showed the existence of a strong vertical anisotropy of the aquifer. The calibrated horizontal permeability is systematically higher than the vertical permeability (≈1000:1). The calibrated model was used to explore the recharge from/or groundwater discharge to Lake Nasser at the seven transects for a 40-year period, i.e. from 1965 to 2004. The analysis for the last 20-year period, 1985–2004, revealed that recharge from Lake Nasser reduced by 37% compared to the estimates for the first 20-year period, 1965–1984. In the period 1965–2004, seepage of Lake Nasser to the surrounding was estimated at 1.15 × 109 m3 year−1. This led to a significant rise of the groundwater table. Variance-based sensitivity and uncertainty analysis on the Edelman results were conducted applying quasi-Monte Carlo sequences (Latin Hypercube sampling). The maximum standard deviation of the total uncertainty on the groundwater table was 0.88 m at Toshka (west of the lake). The distance from the lake, followed by the storage coefficient and hydraulic conductivity, were identified as the most sensitive parameters. Copyright © 2012 John Wiley & Sons, Ltd.

Assessment of hydrological model predictive ability given multiple conceptual geological models

In this study six hydrological models that only differ with respect to their conceptual geological models are established for a 465 km2 area. The performances of the six models are evaluated in differential split‐sample tests against a unique data set with well documented groundwater head and discharge data for different periods with different groundwater abstractions. The calibration results of the six models are comparable, with no model being superior to the others. Though, the six models make very different predictions of changes in groundwater head and discharges as a response to changes in groundwater abstraction. This confirms the utmost importance of the conceptual geological model for making predictions of variables and conditions beyond the calibration situation. In most cases the observed changes in hydraulic head and discharge are within the range of the changes predicted by the six models implying that a multiple modeling approach can be useful in obtaining more robust assessments of likely prediction errors. We conclude that the use of multiple models appear to be a good alternative to traditional differential split‐sample schemes. A model averaging analysis shows that model weights estimated from model performance in the calibration or validation situation in many cases are not optimal for making other predictions. Hence, the critical assumption that is always made in model averaging, namely that the model weights derived from the calibration situation are also optimal for model predictions, cannot be assumed to be generally valid.

Large-scale groundwater modeling using global datasets: a test case for the Rhine-Meuse basin

Abstract. The current generation of large-scale hydrological models does not include a groundwater flow component. Large-scale groundwater models, involving aquifers and basins of multiple countries, are still rare mainly due to a lack of hydro-geological data which are usually only available in developed countries. In this study, we propose a novel approach to construct large-scale groundwater models by using global datasets that are readily available. As the test-bed, we use the combined Rhine-Meuse basin that contains groundwater head data used to verify the model output. We start by building a distributed land surface model (30 arc-second resolution) to estimate groundwater recharge and river discharge. Subsequently, a MODFLOW transient groundwater model is built and forced by the recharge and surface water levels calculated by the land surface model. Results are promising despite the fact that we still use an offline procedure to couple the land surface and MODFLOW groundwater models (i.e. the simulations of both models are separately performed). The simulated river discharges compare well to the observations. Moreover, based on our sensitivity analysis, in which we run several groundwater model scenarios with various hydro-geological parameter settings, we observe that the model can reasonably well reproduce the observed groundwater head time series. However, we note that there are still some limitations in the current approach, specifically because the offline-coupling technique simplifies the dynamic feedbacks between surface water levels and groundwater heads, and between soil moisture states and groundwater heads. Also the current sensitivity analysis ignores the uncertainty of the land surface model output. Despite these limitations, we argue that the results of the current model show a promise for large-scale groundwater modeling practices, including for data-poor environments and at the global scale.

Evaluation of Climate Input Biases and Water Balance Issues Using a Coupled Surface–Subsurface Model

A physically based, coupled and distributed hydrologic model has been set up for the Ringkøbing Fjord catchment, Denmark. This transient model, built with the MIKE SHE/MIKE 11 code, comprises all major components of the terrestrial water cycle, including a three‐dimensional finite‐difference groundwater model. The dynamic coupling of the hydrologic processes secures physically sound feedback and makes the model an ideal tool for evaluating the overall water balance and quantifying potential water balance issues. Historically, failure to obtain water balance closure has been a persistent and much debated problem in Denmark, presumably arising mainly from uncertainties in precipitation, actual evapotranspiration, and groundwater flow to the sea. In this study, the water balance issues were addressed within the modeling framework through analysis of different rain gauge catch corrections and potential evapotranspiration input. The analysis focused on the effect of different rain gauge catch corrections on the model performance, the optimized parameter sets, and state variables not included in the model calibration. The results suggest that water balance problems can be reduced by using a dynamic rain gauge catch correction based on daily wind speed and temperature fields. The model optimization and performance evaluation revealed, however, that several parameter sets gave similar performance compared with the observed groundwater head and river discharge data but still resulted in significant differences with respect to internal water fluxes such as evapotranspiration, groundwater recharge, and stream flow components. New observational data are needed to constrain the model further and thus reduce the water balance uncertainties convincingly.

Importance of considering intraborehole flow in solute transport modeling under highly dynamic flow conditions

Correct interpretation of tracer test data is critical for understanding transport processes in the subsurface. This task can be greatly complicated by the presence of intraborehole flows in a highly dynamic flow environment. At a new tracer test site (Hanford IFRC) a dynamic flow field created by changes in the stage of the adjacent Columbia River, coupled with a heterogeneous hydraulic conductivity distribution, leads to considerable variations in vertical hydraulic gradients. These variations, in turn, create intraborehole flows in fully-screened (6.5 m) observation wells with frequently alternating upward and downward movement. This phenomenon, in conjunction with a highly permeable aquifer formation and small horizontal hydraulic gradients, makes modeling analysis and model calibration a formidable challenge. Groundwater head data alone were insufficient to define the flow model boundary conditions, and the movement of the tracer was highly sensitive to the dynamics of the flow field. This study shows that model calibration can be significantly improved by explicitly considering (a) dynamic flow model boundary conditions and (b) intraborehole flow. The findings from this study underscore the difficulties in interpreting tracer tests and understanding solute transport under highly dynamic flow conditions.

Bayesian model averaging assessment on groundwater management under model structure uncertainty

This study introduces Bayesian model averaging (BMA) to deal with model structure uncertainty in groundwater management decisions. A robust optimized policy should take into account model parameter uncertainty as well as uncertainty in imprecise model structure. Due to a limited amount of groundwater head data and hydraulic conductivity data, multiple simulation models are developed based on different head boundary condition values and semivariogram models of hydraulic conductivity. Instead of selecting the best simulation model, a variance-window-based BMA method is introduced to the management model to utilize all simulation models to predict chloride concentration. Given different semivariogram models, the spatially correlated hydraulic conductivity distributions are estimated by the generalized parameterization (GP) method that combines the Voronoi zones and the ordinary kriging (OK) estimates. The model weights of BMA are estimated by the Bayesian information criterion (BIC) and the variance window in the maximum likelihood estimation. The simulation models are then weighted to predict chloride concentrations within the constraints of the management model. The methodology is implemented to manage saltwater intrusion in the “1,500-foot” sand aquifer in the Baton Rouge area, Louisiana. The management model aims to obtain optimal joint operations of the hydraulic barrier system and the saltwater extraction system to mitigate saltwater intrusion. A genetic algorithm (GA) is used to obtain the optimal injection and extraction policies. Using the BMA predictions, higher injection rates and pumping rates are needed to cover more constraint violations, which do not occur if a single best model is used.

Evaluation of aquifer heterogeneity effects on river flow loss using a transition probability framework

River‐aquifer exchange is considered within a transition probability framework along the Rio Grande in Albuquerque, New Mexico, to provide a stochastic estimate of aquifer heterogeneity and river loss. Six plausible hydrofacies configurations were determined using categorized drill core and wetland survey data processed through the TPROGS geostatistical package. A base case homogeneous model was also constructed for comparison. River loss was simulated for low, moderate, and high Rio Grande stages and several different riverside drain stage configurations. Heterogeneity effects were quantified by determining the mean and variance of the K field for each realization compared to the root‐mean‐square (RMS) error of the observed groundwater head data. Simulation results showed that the heterogeneous models produced smaller estimates of loss than the homogeneous approximation. Differences between heterogeneous and homogeneous model results indicate that the use of a homogeneous K in a regional‐scale model may result in an overestimation of loss but comparable RMS error. We find that the simulated river loss is dependent on the aquifer structure and is most sensitive to the volumetric proportion of fines within the river channel.

Use of heat to estimate streambed fluxes during extreme hydrologic events

Using heat as a tracer, quantitative estimates of streambed fluxes and the critical stage for flow reversal were calculated for high‐flow events that occurred on the Bogue Phalia (a tributary of the Mississippi River) following the 2005 Hurricanes Katrina and Rita. In June 2005, piezometers were installed in the Bogue Phalia upstream from the stream gage near Leland, Mississippi, to monitor temperature. Even with the hurricanes, precipitation in the Bogue Phalia Basin for the months of June to October 2005 was below normal, and consequently, streamflow was below the long‐term average. Temperature profiles from the piezometers indicate that the Bogue Phalia was a gaining stream during most of this time, but relatively static streambed temperatures suggested long‐term data was warranted for heat‐based estimates of flux. However, the hurricanes caused a pair of sharp rises in stream stage over short periods of time, increasing the potential for rapid heat‐based modeling and for identification of the critical stage for flow reversal into the streambed. Heat‐based modeling fits of simulated‐to‐measured sediment temperatures show that once a critical stage was surpassed, flow direction reversed into the streambed. Results of this study demonstrate the ability to constrain estimates of streambed water flux and the critical stage of flow reversal, with little available groundwater head data, by using heat as a tracer during extreme stage events.

갑천 유역의 지하수 유동 평가

지하수 유동 모델에 반영되는 유역의 수리지질학적, 수문학적 특징들은 지하수 유동 특성에 주요한 영향을 미칠 수 있는 현장 특성들이다. 상기 특징들을 최대한 반영하기 위하여 438개의 관정자료를 토대로 갑천 유역(유역 면적<TEX>$648.3km^2$</TEX>)의 대수층 구조 특성을 분석하고, 24개의 하천 자료를 이용하여 지하수 유동 모델을 구축하였다. 그리고 검 보정된 준분포형 유출모형(SWAT)의 소유역별 지하수 함양량 결과를 3차원 지하수 유동 모델(MODFLOW)과 연계하여 갑천 유역의 광역 지하수 유동 특성을 평가하였다. 모의된 지하수위와 86개 지하수 관측정의 지하수위 비교에서 결정계수는 0.99, 유역 전체의 물교환량의 상대오차율이 약 0.57%로 갑천 유역의 지하수 유동 특성을 잘 반영하였다. 갑천 유역의 지하수는 지형 및 대수층 특성과 하천의 영향에 의하여 전반적으로 유역 남쪽에서 북쪽으로 유동하는데 산지 지역에 존재하는 하천들은 손실과 이득 하천의 형태를 반복하는 반면, 갑천 중 하류 부분과 갑천-유등천 합류 부근은 이득하천의 형태를 보이고 있다. 그리고 소유역별 지하수 물교환량 분석한 결과, 상기 지역을 포함하는 소유역의 하천 공급원이 해당 유역의 지하수 함양량보다 지하수 유동 체계에 따라 인근 소유역에서 공급되는 지하수 유입량 비율이 훨씬 높은 것으로 분석되었다. 반면 유등천 중 상류 지역의 소유역을 제외한 산지 지역의 경우는 하천 공급원이 지하수 유동 체계에 의한 유입량보다 해당 유역 내에서 공급되는 지하수 함양량 비율이 높은 것으로 분석되었다. 따라서 지표수와 지하수의 상호작용 및 하천 유출에 지하수 유동 특성이 주요한 영향을 미치는 것으로 평가되었다. Groundwater flow in a basin is greatly affected by many hydrogeological and hydrological characteristics of the basin. A groundwater flow model for the Kap-cheon basin (<TEX>$area=648.3km^2$</TEX>) in the Geum river basin was established using MODFLOW by fully considering major features obtained from observed data of 438 wells and 24 streams. Furthermore, spatial groundwater recharge distribution was estimated employing accurately calibrated watershed model developed using SWAT, a physically semi-distributed hydrological model. Model calibration using observed groundwater head data at 86 observation wells yielded the deterministic coefficient of 0.99 and the water budget discrepancy of 0.57%, indicating that the model well represented the regional groundwater flow in the Kap-cheon basin. Model simulation results showed that groundwater flow in the basin was strongly influenced by such factors as topological features, aquifer characteristics and streams. The streams in mountainous areas were found to alternate gaining and losing steams, while the streams in the vicinity of the mid-stream and down-stream, especially near the junction of Kap-cheon and Yudeong-cheon, areas were mostly appeared as gaining streams. Analysis of water budget showed that streams in mountainous areas except for the mid-stream and up-stream of Yudeong-cheon were mostly fed by groundwater recharge while the streams in the mid and down-stream areas were supplied from groundwater inflows from adjacent sub-basins. Hence, it was concluded that the interactions between surface water-groundwater in the Kap-cheon basin would be strongly inter-connected with not only streams but also groundwater flow system itself.

Impact of the shift from groundwater to surface water irrigation on aquifer dynamics and hydrochemistry in a semi-arid region in the south of Portugal

In semi-arid regions, the shift from locally extracted groundwater to regionally supplied surface water for irrigation can have a large impact on the aquifer dynamics and hydrochemistry, as is shown for an area of intensive citrus culture in the south of Portugal. Studying the changes that occur is important, particularly in the light of European and national policies that seek to preserve the quality of groundwater and the involved ecosystems. In September 2005, the application of Nitrates Directive 91/676/CEE to the study area led to its designation as a vulnerable zone and the obligation of implementing measures that reduce the risk of nitrate leaching to groundwater. In addition, the area is bordered to the south by the internationally recognized ecosystem of the Ria Formosa lagoon, which requires a careful management of water quantity and quality. This article shows that the substitution of the irrigation source indirectly caused the reduction of the nitrogen load on soil and groundwater. The study is based on a large number of groundwater head and quality data gathered over the past decade. It is observed that irrigation with surface water triggered freshening of the upper aquifer, inverting a long-existing trend of increasing groundwater salinities caused by irrigation with local groundwater. The occurrence of both phenomena in the same aquifer is rather unique and is simulated with the hydrogeochemical software PHREEQC. In terms of aquifer dynamics, a sharp rise of the water table is observed, together with an almost complete attenuation of its interseasonal oscillations. A further rise, which is to be expected in years with average or higher rainfall, could cause the water table to enter the root zone of the citrus trees in some areas, with potentially negative consequences. Therefore, an integrated mixed-source irrigation system using both surface and groundwater is highly recommended.

Time series forecasting by combining the radial basis function network and the self‐organizing map

AbstractBased on a combination of a radial basis function network (RBFN) and a self‐organizing map (SOM), a time‐series forecasting model is proposed. Traditionally, the positioning of the radial basis centres is a crucial problem for the RBFN. In the proposed model, an SOM is used to construct the two‐dimensional feature map from which the number of clusters (i.e. the number of hidden units in the RBFN) can be figured out directly by eye, and then the radial basis centres can be determined easily. The proposed model is examined using simulated time series data. The results demonstrate that the proposed RBFN is more competent in modelling and forecasting time series than an autoregressive integrated moving average (ARIMA) model. Finally, the proposed model is applied to actual groundwater head data. It is found that the proposed model can forecast more precisely than the ARIMA model. For time series forecasting, the proposed model is recommended as an alternative to the existing method, because it has a simple structure and can produce reasonable forecasts. Copyright © 2005 John Wiley &amp; Sons, Ltd.