Assessing distributed groundwater recharge rate using integrated surface water-groundwater modelling: application to Mihocheon watershed, South Korea

A method of estimating groundwater recharge, based on water-balance components using the SWAT-MODFLOW model (an integrated surface water-groundwater model), is described. A multi-reservoir storage routing module is suggested instead of a single storage routing module in SWAT; this represents a more realistic delay in the travel of water through the vadose zone. By using this module, the parameter related to the delay time can be optimized by checking the correlation between simulated recharge and observed groundwater levels. The final step of this procedure is to compare simulated groundwater levels as well as the simulated watershed stream flow with the observed groundwater levels and watershed stream flow. This method is applied to the Mihocheon watershed in South Korea to estimate spatio-temporal groundwater recharge distribution. The computed annual recharge rate is compared with the independently estimated recharge rate using BFLOW. The hydrologic modelling results show that the annual average recharge rate should be estimated by a long-term continuous simulation with a distributed hydrologic modelling technique.     

Determination of areas vulnerable to pollution in a karstic river basin in Turkey via a decision support system based on DRASTIC,SWAT and isotopes analysis

Determination of areas vulnerable to water pollution in river basins helps to generate appropriate water management protection plans. This study aims to define areas vulnerable to pollutants in a data-scarce karstic river basin in Turkey by using a holistic approach integrating the Soil and Water Assessment Tool (SWAT), the DRASTIC framework, and selected isotopes within a decision support system based on a geographic information system to delineate vulnerable areas. DRASTIC was used to show groundwater vulnerability to pollutants. The concentrations of isotopes ¹⁸O, ²H, and ³H in groundwater were used to define the vulnerable areas of the karst region. SWAT was utilized to show watershed vulnerability to pollutants in shallow aquifers. The recharge rate parameter in DRASTIC was obtained from SWAT. This methodological approach was applied to the Yuvacık Dam Basin in Kocaeli, part of the Marmara River Basin, as it is a good example of a karstic watershed. According to this study, each approach provides different vulnerabilities when applied separately. The final map obtained from the integrated approach shows that drinking water supplies in the northeast and northwest parts of the basin are highly vulnerable to pollution. All the karst spring catchments and areas near the basin outlet are highly vulnerable. Moreover, across all water samples taken across the basin, those exhibiting the highest concentrations in nitrate were all found in the areas mapped as highly vulnerable. The methodology was validated by analyzing nitrate concentration in 22 groundwater and surface-water samples.

     

东北半干旱地区流域分布式水文模拟——以洮儿河流域为例

以东北半干旱地区典型流域-洮儿河流域为研究对象,应用SWAT模型对流域水文过程进行了模拟研究;选择流域上游子流域和中下游子流域分别进行参数敏感性分析,识别出影响模拟结果的敏感参数,研究发现部分参数敏感性存在空间变异性,分析主要原因在于气候和下垫面的空间异质性导致了流域上下游产流模式存在差异。采用1988-1997年水文气象数据进行模型率定和验证,结果表明：干流水文站月流量过程率定期Nash-Sutcliffe 效率系数平均值为0.78,验证期为0.72,相关系数都达到0.86以上,水量误差大多在20%以内,对日过程的模拟也有较高的精度;枯水年模拟结果较差,主要是因为流域降水站数量不够,难以反映降水的时空分布。对于水文、气象等资料相对缺乏的东北半干旱地区,SWAT模型的模拟结果总体令人满意,可以应用于与流域径流相关的各种模拟分析,研究成果对进一步加强洮儿河流域水资源综合管理提供了依据和手段。     

Numerical modelling of the saline interface in coastal karstic aquifers within a conceptual model uncertainty framework

Numerical modelling is increasingly used as a tool for improving management strategies in aquifers and to support the design of comprehensive projects considering natural and anthropogenic processes. Overall, numerical simulation in karstic aquifers poses a major scientific challenge due to the non-Darcian groundwater flow dynamics. In specific cases, the equivalent porous medium approach has shown acceptable results, particularly in poorly karstified aquifers with regional/subregional scales such as this case. The Yucatan coastal karstic aquifer (Mexico) has been defined as a complex regional heterogeneous system, partially confined, thus allowing the discussion of multiple conceptual models. In this research, a two-dimensional numerical model of flow and transport was implemented using SEAWAT for the NW Yucatan aquifer. Four likely conceptual models were audited, calibrated and verified using hydrogeological field data, to select the best one, considering their fit and complexity. The numerical model accuracy was evaluated using the root-mean-square error, Nash Sutcliffe efficiency and the Pearson coefficient. The Akaike information criterion and Bayesian information criterion were included for evaluating the complexity of the numerical models. In addition, the signal of tide propagation into the aquifer was assessed as a proxy to improve the numerical calibration process. Results show that the most complex numerical model has a better calibration than the simpler models, but the model accuracy is worse when compared to less complex numerical models in the verification exercise. This research offers enhancement in the knowledge of numerical modelling in heterogeneous coastal aquifers within a conceptual-model uncertainty setting.

     

Evaluation of the SWAT model for water balance study of a mountainous snowfed river basin of Nepal

In this study, a semi-distributed hydrologic model Soil and Water Assessment Tool (SWAT) has been employed for the Karnali River basin, Nepal to test its applicability for hydrological simulation. Further, model was evaluated to carry out the water balance study of the basin and to determine the snowmelt contribution in the river flow. Snowmelt Runoff Model (SRM) was also used to compare the snowmelt runoff simulated from the SWAT model. The statistical results show that performance of the SWAT model in the Karnali River basin is quite good (p-factor = 0.88 and 0.88, for daily calibration and validation, respectively; r-factor = 0.76 and 0.71, for daily calibration and validation, respectively). Baseflow alpha factor (ALPHA_BF) was found most sensitive parameter for the flow simulation. The study revealed that the average annual runoff volume available at the basin outlet is about 47.16 billion cubic metre out of which about 12% of runoff volume is contributed by the snowmelt runoff. About 25% of annual precipitation seems to be lost as evapotranspiration. The results revealed that both the models, SWAT and SRM, can be efficiently applied in the mountainous river basins of Nepal for planning and management of water resources.     

考虑径流成分的雅鲁藏布江月径流模拟

通过构建雅鲁藏布江流域分布式水文模型SWAT,并以1989~1998年为率定期,1999~2005年为验证期对3个水文站(奴各沙、羊村、奴下)的月径流进行了模拟。从模型参数组里选取径流成分模拟不同的2组参数,以Nash-Sutcliffe效率系数(NS)、相关系数(R~2)及相对误差(PBIAS)为评价指标,来分析模型的模拟结果好坏。结果表明:率定期内,在径流成分1下奴各沙、羊村站和奴下站的R~2值高达0.9以上,NS系数分别为0.9、0.92和0.72,PBIAS的值为7.1%、18.5%及34.2%,但径流成分与实际情况出入较大。而径流成分2下,模拟得到的评价指标虽然不是很高,但径流成分比较合理。验证期总体结论一致。因此,建议在模型率定验证中要充分考虑径流成分的作用,从而提高模型模拟的精度。     

SWAT模型在东江流域的应用研究

以SWAT模型为研究工具,对该模型在东江流域的3个子流域的适用性进行研究,研究表明,SWAT模型在东江流域3个子流域的适用性较好,校准期（19701975年）日径流模拟和月径流模拟的相对误差Re均在10%以内,日径流模拟的决定系数R2均在70%以上,月径流模拟的决定系数R2均在80%以上,Nash-Suttcliffe效率系数基本达到70%;验证期（19761985年及19962005年）,月径流模拟的决定系数R2和Nash-Suttcliffe效率系数均在70%以上,相对误差Re基本在±20%以内,可以满足该流域的水资源评价与规划的要求。     

SWAT manual calibration and parameters sensitivity analysis in a semi-arid watershed in North-western Morocco

Being a laborious approach, manual calibration of hydrologic model in a semi-arid context requires in-depth knowledge of the watershed and as much as possible field input data to obtain reliable simulations. In this study, manual calibration and relative sensitivity analysis approaches of the SWAT model (Soil and Water Assessment Tool) were applied for water balance in a 1993 km² watershed (on the R’dom river) located in North-western Morocco. The watershed is located in a semi-arid area dominated by agro-forestry activities. The objectives of this study were (i) to perform a local sensitivity analysis of the SWAT model taking into consideration the watershed characteristics and (ii) to implement a detailed methodology of manual calibration and validation of the model in a semi-arid context. Sensitivity analysis has been carried out on 12 different SWAT input parameters, and has revealed that 4 input parameters only were the most influential ones on flow components of the R’dom watershed. Model manual calibration was conducted along 2006 and 2007 by comparing measured and predicted monthly and daily discharges and taking Nash-Sutcliffe coefficient (NSE), determination coefficient (R ²), and percent bias (PBIAS) as goodness-of-fit indicators. Validation has been performed by the same approach through 2008 and 2009 period. All final NSE values were above 0.5, R ² values exceeded 0.7, and PBIAS lower than 25% demonstrating satisfactory model performances over the study watershed conditions. The SWAT model set-up with measured input data, manually calibrated and validated, reflects well the real hydrologic processes occurring in the R’dom watershed and can be used to assess current and future conditions and to evaluate alternative management practices.     

SWAT模型在海拉尔河流域径流模拟中的应用研究

以海拉尔河上游流域作为研究区域,基于Arc GIS构建SWAT分布式水文模型对流域水文过程进行模拟,通过对流域的基础数据整合,模型采用1999~2003年实测径流数据进行参数率定,将2004~2010年实测径流数据作为模型的验证期,对模型在海拉尔河上游的适用性进行研究。通过对月和年径流模拟值和实测值的比较,率定期和验证期的Nash系数Ens和相关系数R2分别在0.861~0.873和0.877~0.899之间。基于这两个评价标准可知:SWAT模型在海拉尔河上游流域有良好的适用性,可以为该流域的水资源管理提供依据。     

Parameters optimization based on the combination of localization and auto-calibration of SWAT model in a small watershed in Chinese Loess Plateau

This study simulated the watershed flow and sediment responses based on calibration of the SWAT model in the semi-arid Chinese Loess Plateau (LP) where soil erosion intensively occurs. After the model’s initiation and manual modification, a 7-year inconsecutively observed flow and sediment data from 1984 to 1990 was used to analyze the model’s application in the selected watershed called AJW in the Chinese LP region. The model procedure included sensitivity analysis, parameter calibration and model validation. The best parameter set was finally determined based on the combination of parameter localization and auto-calibration. Then the model was assessed for its accuracy based on the NSE estimation, resulting in 0.77 and 0.67 for calibration and 0.46 and 0.32 for validation on simulations for flow and sediment, respectively, which is a moderately satisfactory accuracy among the applications of the SWAT model. Annual watershed assessment on flow and sediment with the calibrated SWAT model resulted in a multiyear averaged annual runoff coefficient of about 2.7% and an erosion modulus of 797 t/(km²·a⁻¹) in the AJW, indicating a beneficial consequence from the implementation of the historical soil and water conservations.     

A semi-distributed groundwater recharge model for estimating water-table and water-balance variables

A semi-distributed groundwater recharge model is presented, which estimates water-table fluctuation and water-balance variables. The model is expressed by the water-balance concept linking atmospheric and hydrogeological parameters to different water uses (industrial, agricultural, domestic, etc). It was calibrated and validated using 5 years of data collected in the Dogo Plain in Japan. A 3-year dataset, from 2000 to 2002, was used in the calibration, while a 2-year dataset, from 2003 to 2004, was used for the validation. Calibration of the model was achieved by the shuffled complex evolution automatic optimization of model parameters to match simulated results with measured water-table depth. Square roots of relative error (R²) are 0.88 and 0.90 for calibration and validation processes, respectively. Monthly evolution of water storage change was then estimated and the water-table drawdown in different pumping scenarios was simulated. Finally, the groundwater-pumping amount planned by the government for future sustainable groundwater utilization was evaluated. The government-planned groundwater-pumping amount is feasible in most regions while the midstream region should be paid more attention. This study offers a scientific basis to control and prevent depletion of groundwater resources.     

Development of a coupled surface-groundwater-pipe network model for the sustainable management of karstic groundwater

This paper considers the hydrogeological simulation of groundwater movement in karstic regions using a hydrological modelling system (SHETRAN) which has been adapted for modelling flow in karstic aquifers. Flow and transport through karstic aquifers remains poorly understood, yet quantitative hydrogeological models are essential for developing and implementing groundwater protection policies. The new model has been developed and used within the STALAGMITE (Sustainable Management of Groundwater in Karstic Environments) project, funded by the European Commission. The SHETRAN model is physically based insofar as most of the parameters have some physical meaning. The SHETRAN model represents all of the key processes in the hydrological cycle, including subsurface flow in the saturated and unsaturated zones, surface flow over the ground surface and in channels, rainfall interception by vegetation canopies, evapotranspiration, snow-pack development and snowmelt. The modifications made to SHETRAN to simulate karstic aquifers are (1) the coupling of a pipe network model to a variably saturated, three-dimensional groundwater component (the VSS-NET component), to simulate flow under pressure in saturated conduits; (2) the coupling of surface water features (e.g. sinking streams or "ponors", and spring discharges) to the conduit system; (3) the addition of a preferential "bypass" flow mechanism to represent vertical infiltration through a high-conductivity epikarst zone. Lastly, a forward particle tracking routine has been developed to trace the path of hypothetical particles with matrix and pipe flow to springs or other discharge points. This component allows the definition of groundwater protection zones around a source for areas of the catchment (watershed) which are vulnerable to pollution from non-point sources (agriculture and forestry).     

Investigating the hydrogeological properties of springs in a karstic aquifer in Dorfak region (Guilan Province,Iran)

Karstic aquifers are considered as the main sources of groundwater in the northeast of Rudbar, Iran. The present study was conducted to evaluate the hydrogeological properties of karstic springs in this region. For this purpose, saturation indices (SI values) were calculated using the geochemical PHREEQC model for a number of minerals in the groundwater in the karstic aquifer. Moreover, AqQA-RockWare software packages were used to prepare hydrogeochemical plots for the aquifer, using which the sources of the ions in the water were identified. The origin of bicarbonate, calcium, and magnesium ions in water was determined using chloro-alkaline indices. Moreover, through plotting a Piper diagram for spring water samples, it was discovered that water type of all springs is the Ca-HCO₃ type, confirming the karstic characteristic of springs in the area. A Durov diagram also suggests that the water composition of the springs is of the bicarbonate type with the dominant Ca cation, suggesting the calcareous effects of the region on the quality of groundwater and exhibiting a single source for the springs. The calculated saturation indices show that most of the water samples are undersaturated with respect to calcite, dolomite, and CO₂. The stable isotopes (δ¹⁸O and δ²H) and deuterium excess values were used to get information about transport pathways in groundwater, atmospheric moisture, and the degree of interaction between these reservoirs. The degree of karstification of the recharge area of the karst aquifer was determined to be 5.5 from an analysis of the hydrograph Sefidab Spring.     

NARX neural network approach for the monthly prediction of groundwater levels in Sylhet Sadar,Bangladesh

Groundwater is important for managing the water supply in agricultural countries like Bangladesh. Therefore, the ability to predict the changes of groundwater level is necessary for jointly planning the uses of groundwater resources. In this study, a new nonlinear autoregressive with exogenous inputs(NARX) network has been applied to simulate monthly groundwater levels in a well of Sylhet Sadar at a local scale. The Levenberg-Marquardt(LM) and Bayesian Regularization(BR) algorithms were used to train the NARX network, and the results were compared to determine the best architecture for predicting monthly groundwater levels over time. The comparison between LM and BR showed that NARX-BR has advantages over predicting monthly levels based on the Mean Squared Error(MSE), coefficient of determination(R~2), and Nash-Sutcliffe coefficient of efficiency(NSE). The results show that BR is the most accurate method for predicting groundwater levels with an error of ± 0.35 m. This method is applied to the management of irrigation water source, which provides important information for the prediction of local groundwater fluctuation at local level during a short period.     

济宁地区北部岩溶地下水资源计算与潜力评价

济宁地区北部赋存有较丰富的岩溶地下水,目前已纳入当地开发利用规划。为防止过量开采引发不良环境地质问题,对该岩溶系统的资源潜力进行了科学评价。在查明济宁地区北部水文地质条件的基础上,利用水量均衡法与数值模拟法对该区岩溶地下水资源量进行了计算,通过分析地下水资源潜力模数,该区岩溶地下水资源潜力较大。     

Determination of water balance equation components in irrigated agricultural watersheds using SWAT and MODFLOW models : A case study of Samalqan plain in Iran

Increasing water demands,especially in arid and semi-arid regions,continuously exacerbate groundwater as the only reliable water resources in these regions.Samalqan watershed,Iran,is a groundwater-based irrigation watershed,so that increased aquifer extraction,has caused serious groundwater depletion.So that the catchment consists of surface water,the management of these resources is essential in order to increase the groundwater recharge.Due to the existence of rivers,the low thickness of the alluvial sediments,groundwater level fluctuations and high uncertainty in the calculation of hydrodynamic coefficients in the watershed,the SWAT and MODFLOW models were used to assess the impact of irrigation return flow on groundwater recharge and the hydrological components of the basin.For this purpose,the irrigation operation tool in the SWAT model was utilized to determine the fixed amounts and time of irrigation for each HRU(Hydrological Response Unit)on the specified day.Since the study area has pressing challenges related to water deficit and sparsely gauged,therefore,this investigation looks actual for regional scale analysis.Model evaluation criteria,RMSE and NRMSE for the simulated groundwater level were 1.8 m and 1.1%respectively.Also,the simulation of surface water flow at the basin outlet,provided satisfactory prediction(R²=0.92,NSE=0.85).Results showed that,the irrigation has affected the surface and groundwater interactions in the watershed,where agriculture heavily depends on irrigation.Annually 11.64 Mm3 water entered to the aquifer by surface recharge(precipitation,irrigation),transmission loss from river and recharge wells 5.8 Mm3 and ground water boundary flow(annually 20.5 Mm³).Water output in the watershed included ground water extraction and groundwater return flow(annually 46.4 Mm³)and ground water boundary flow(annually 0.68 Mm³).Overally,the groundwater storage has decreased by 9.14 Mm3 annually in Samalqan aquifer.This method can be applied to simulate the effects of surface water fluxes to groundwater recharge and river-aquifer interaction for areas with stressed aquifers where interaction between surface and groundwater cannot be easily assessed.     

随机森林模型预测岩溶区酸性煤矿井水锰污染

酸性煤矿井水严重威胁地下水的水质。如何更有效对受影响区域的地下水源进行动态监测是当前的一个关键问题。采用随机森林中的回归模型,利用自变量(采空区水位、岩溶水位、pH值、泉水流量、电导率)和因变量(污染离子浓度)的相关性,建立回归模型;使用测试数据进行误差分析,结果证明模型准度较高,所得预测值具有参考价值;得出各自变量对因变量影响的重要程度,分析结果与实际情况相符合。试验表明,随机森林回归模型在酸性煤矿井水污染预测方面具有适用性,可作为辅助手段监测水质污染情况,对今后工作有一定的指导意义和经济价值。     

Modeling the risk of nitrate leaching and nitrate runoff loss from intensive farmland in the Baiyangdian Basin of the North China Plain

The water movement and soil nitrogen cycle of the Baiyangdian Basin were simulated, and the risk of nitrate leaching and nitrate runoff loss from intensive farmland was assessed by using the distributed hydrological soil and water assessment tool (SWAT) in this study. The model assessment showed that SWAT was able to simulate water and nitrate movement in the region with satisfactory results. The modeling analysis indicated that fertilizer application was the overriding source of soil nitrogen and might result in a large amount of nitrate accumulation in soils; this nitrate might be lost by leaching or runoff driven by water movement. In 2009, nitrate nitrogen leaching represented 19.5 % of the total amount of nitrogen fertilizer application, while nitrate nitrogen runoff represented 1.7 % only. Thus, it showed that the nitrate leaching was the main approach of soil nitrogen movement in farmland because of strong percolation. It also showed a significant variation of nitrate leaching from different soil depths, with the largest amount leached from surface soil layers and the smallest amount leached from lower soil layers. Therefore, it could be further revealed that the nitrate concentration was very low at soil layers lower than the root zone of crops (1.2 m). Validated by groundwater observations, groundwater pollution by nitrate derived from fertilizers was not serious because of the deep groundwater level in the study plain. However, the risk of groundwater pollution would increase significantly if precipitation increased.     

Numerical simulation and environmental impact prediction of karst groundwater in Sangu Spring Basin,China

The changes of development and utilization of karst groundwater in Sangu Spring Basin have made the original groundwater resource evaluation unable to meet the needs of future economic development. Based on analysis of existing data, combined with the characteristics of supplement, runoff and draining of regional karst groundwater, the Visual Modelflow software was used to build a numerical simulation model of Sangu spring Basin. The amount of karst groundwater resource and groundwater environment of the Basin were evaluated under different exploitation schemes, and the changes of karst groundwater environment in the future ten years were also predicted. The fitting error which is less than 0.5 m between the calculated value and measured value of the water level in the fitted borehole accounts for 93%. For the lithologically and structurally complex Sangu Spring Basin, the fitting effect of numerical simulation model was ideal. On the basis of the current mining amount of 111.80 million m³/a, the total redistributed exploited amount in the spring region was 61.79 million m³/a. Under the condition that the quantity of recoverable resources reached 173.59 million m³/a and under different precipitation schemes, all constraint conditions were satisfied, such as regional water level drawdown, maximum allowable water level drawdown in every simulated water source area and the flow rate of Guobi Spring. The results will provide a scientific basis for the rational development and utilization of karst groundwater in Sangu Spring Basin.     

Hydrological characteristics of karstic groundwater in the northeast Viet Nam as studied by isotopic techniques

The stable water isotopic composition (δ²H and δ¹⁸Ο), tritium (³H) activity, dissolved organic carbon, alkalinity, as well as the composition of carbon 13 (δ¹³C) in dissolved inorganic carbon (DIC) of 36 water samples taken from 16 resurgences in the northeast provinces of Viet Nam in the dry (Mar 2008) and rainy (June 2008) seasons were analyzed to elucidate hydrological characteristics of the karstic aquifers in the area. The stable water isotopic composition of the water samples collected clearly demonstrated that the karstic groundwater in the region was recharged from the local meteorological water. The tritium activity in the samples was found to be in between 3 and 4 TU, falling in the range of the ³H activity in the local precipitation and thus meaning that the traveling of recharge water to the resurgences was very short. Concentrated and diffuse allogenic recharges seem to be important sources of karstic groundwater in the study region. Water in the karstic aquifers could be classified into three types as: (a) water from karstic areas with dense vegetation cover that causes DIC be depleted in carbon 13 (¹³δ<?12‰ vs. Pee Dee Belemnite standard of Vienna, VPDB); (b) water from karstic areas with poor vegetation cover that originates DIC with carbon 13 composition ranging from ?11 to ?12‰; and (c) surface water from lakes, springs and rivers that has DIC with enriched carbon 13 (δ¹³C >?10‰). This implies that there are several sources of carbon dioxide contributing to the DIC in water of the karstic aquifers in the study region. Among other potential sources, the atmospheric CO₂, CO₂ from carbonate rock dissolution, biomineralization of soil organic matters and plant roots respiration seem to be important sources of the DIC in the waters of this region. The results show high vulnerability towards anthropogenic contaminants of karstic groundwater in the study region.