Integrated modeling as a decision‐aiding tool for groundwater management in a Mediterranean agricultural watershed

A decision‐aiding methodology for agricultural groundwater management is presented; it is based on the combination of a watershed model, a groundwater flow model, and an optimization model. This methodology was applied to an agricultural watershed in northeastern Greece. The watershed model used was the Soil and Water Assessment Tool (SWAT), which provided recharge rates for the aquifers. These recharge rates were imported in the well‐known MODFLOW groundwater flow model. Both models were calibrated and verified using field data. Then, the nonlinear optimization problem was solved by a piecewise linearization process, in which the Simplex algorithm was applied sequentially. Apart from several pumping and climate change sensitivity scenarios, a land use change scenario and a climate change scenario, combining the three models, were tested, showing the ability of this methodology to be used in the decision‐making process. Copyright © 2012 John Wiley & Sons, Ltd.     

Seasonal variation in cascade‐driven hyporheic exchange,northern Honduras

A characterization of hyporheic exchange for dry and wet season baseflow, as well as partially dewatered discharge, was done in Prieta Creek, a first‐order cascade in northern Honduras. The cascade had discharges from 1 to 15 l s^?1, had average slopes of 12%, pool spacing of 3 m, and shallow substrate of sand and gravel. Tracer tests were conducted in a 15‐m sub‐reach, a length considered to be adequate for the experiment based on the DaI test, a ratio of exchange and transport processes. In the three tests, between 9 and 18% of tracer was not recovered, possibly due to entrainment in flowpaths passing beneath the downstream monitoring location. Tracer data were analysed by the one‐dimensional transport with inflow and storage (OTIS) transient storage model (TSM) to derive standard exchange parameters, and by the solute transport in rivers (STIR) model to examine hyporheic residence time distributions (RTDs). The best fit of the observed tracer breakthrough curves was obtained by using the STIR model with a combination of two exponential RTDs to represent hyporheic retention. With increasing discharge, the OTIS model predicted increasing storage exchange fluxes and exchange coefficients and decreasing storage zone areas and transient storage times, which are trends supported by riparian and streambed piezometric head data. Riparian water levels rose during the transition from the dry to wet season, which could constrict the hyporheic storage zone. Thirteen of the 19 streambed piezometers recorded seasonal changes in hydraulic gradients and flux direction, with fewer yet stronger upwelling zones during higher discharges. The MODFLOW model missed the observed seasonal changes, possibly due to subtle changes in the seasonal change in water surface profiles. We conclude that partially dewatered dry season exchange, compared to wet season exchange, was initiated and terminated with smaller pressure gradients and, in different streambed locations, was smaller in volume, had longer residence times, and may connect with deeper and longer flow paths. Copyright © 2010 John Wiley & Sons, Ltd.     

MODFLOW在北京某工程抗浮分析的应用

简要介绍地下结构浮力研究方面的一些成果。应用MODFLOW软件进行数值模拟,结合具体工程进行分析论证,提出抗浮优化设防水位,降低抗浮设防投入。     

Groundwater discharge mechanism in semi‐arid regions and the role of evapotranspiration

The present study examined groundwater recharge/discharge mechanisms in the regional Central Sudan Rift Basins (CSRB). Aquifers in CSRB constitute poorly sorted silisiclastics of sand, clay and gravels deposited in closed hydrologic systems of the Cretaceous–Pleistocene fluviolacustrine environments. CSRB are bounded to the north by the highlands of the Central African Shear Zone (CAZS) that represents the surface and groundwater divides. Sporadic recharge in the peripheries of the basins along the CASZ occurs subsequent to decadal and centennial storm events. Inflow from the Nile into the aquifers represents an additional source of recharge. Thus, groundwater resources cannot be labelled fossil nor can they be readily recharged. Closed hydrologic troughs located adjacent to the influent Nile system mark areas of main groundwater discharge characterized by lower hydraulic heads. This study has examined mechanisms that derive the discharge of the groundwater in these closed basins and concluded that only evapotranspirative discharge can provide a plausible explanation. Groundwater abstraction is mainly through deep‐rooted trees and effective evaporation. The increase of TDS along the flow indicates local recharge at the peripheries of basins and shows the influence of evaporation and rock/water interaction. The decline in groundwater level along a flow path was calculated using Darcy's law to estimate average recharge and evapotranspirative discharge, which are equal under natural equilibrium and make the only fluxes in CSRB. Steady‐state 2D flow modelling has demonstrated that an average recharge of 4–8 mm yr^?1 and evapotranspirative discharge of 1–22 mm yr^?1 will maintain natural equilibrium in CSRB. Sporadic storms provide recharge in the highlands to preserve the current hydraulic gradient and maintain aquifer dynamics. Simulated recharge from the Nile totals about 17·5 mm yr^?1 and is therefore a significant contributor to the water balance. Copyright © 2007 John Wiley & Sons, Ltd.     

地下水数值模拟过程中的误差分析

分析了从模型建立到输出过程中的误差来源,得出概念模型误差来源极大依靠资料的详实程度,三维的数学模型仿真程度最高,数值模型是将数学模型的精确解转化为近似解的过程,存在着不可避免的截断误差,模型构建过程中由于剖分尺度的不同而产生误差的程度不同,模型输入过程中定解条件和各项参数的选取较大影响了模型的输出精度。     

Model estimates of runoff in the closed,semiarid Estancia basin,central New Mexico,USA

Abstract

The 5000 km² topographically closed Estancia basin in central New Mexico has been the focus of several palaeoclimatic studies based on changes in the level of late Pleistocene Lake Estancia. A large, unknown volume of surface runoff and groundwater from adjacent mountains contributed to the hydrological balance during highstands and lowstands. The US Department of Agriculture hydrological model SWAT (Soil and Water Assessment Tool) and the US Geological Survey groundwater flow model MODFLOW, with the LAK2 package, were used in this study to estimate runoff and water balance under present climate. A Geographic Information Systems (GIS) interface was used for SWAT, digitized data were applied for soils and vegetation, and limited streamflow data were used to obtain an approximate calibration for the model. Simulated streamflow is generally within 30% of observed values, and simulated runoff for the entire basin is about 8% of the annual inflow volume needed to support lowstands of the former Lake Estancia. Results from the combined models suggest application to other palaeoclimate investigations in semiarid lake basins.     

Spatio-temporal patterns of the interaction between groundwater and surface water in plains

The study of the dynamics of anthropic disturbances that have an effect on the hydrological systems in plains requires integral simulation tools for their diagnosis. The objective of this article is, first, to analyse and reproduce the spatio-temporal interactions between groundwater (GW) and surface water, net recharge, GW level, surface run-off, and evapotranspiration in the upper creek basin of Del Azul, which is located in the centre of the province of Buenos Aires, Argentina, and second, to obtain insights to apply the methodology to other similar situations. For this purpose, a model coupling the semidistributed hydrological model (Soil and Water Assessment Tool [SWAT]) and the hydrogeological model (MODFLOW) has been used. A simulation was carried out for a period of 13 years (2003–2015) on a daily scale. The application of the SWAT–MODFLOW coupling gave good results based on the adjustment between the calculated flows and levels, reaching a Nash–Sutcliffe of 0.6 and R²0.6 at the Seminario hydrometric station located at the watershed outlet point. According to the annual average balance, out of the total rainfall, evapotranspiration accounts for 85%, recharge accounts for 10.2%, and surface run-off accounts for 4.8%. Annual and monthly trends of the stream–aquifer interaction were determined, obtaining on average an annual GW discharge of 34 mm and an annual average recharge of the stream to the aquifer of 1.4 mm. Monthly GW discharges are higher in winter–spring (July to December with an average of 3.3 mm) and lower in summer–autumn (January to June with an average of 2.8 mm). The monthly average recharge of the stream towards the aquifer varies from 0.02 to 0.36 mm and is higher in March, May, and August, when water excess is produced in the basin. Through the analysis of coupled modelling, it is possible to analyse and reproduce the spatio-temporal transitions of flow existing between the stream, the hyporheic zone, and the aquifer.     

MODFLOW不同算法及参数设定对计算精度的影响

为了提高MODFLOW中模型计算的精确度,减少不收敛情况的发生,通过建立数值模型算例,分析了不同解算器对计算精度的影响差异、主要参数值的设定对计算精度产生的影响以及收敛指标与迭代次数的定量关系。结果表明,PCG、SAMG和GMG计算精度较高,WHS和SOR计算精度较低。PCG算法中,低收敛指标使得计算精度有所提高,但稳定性却明显降低;计算精度对PCG算法中的阻尼因子的敏感性比对SIP算法中的加速因子敏感性低,二者的减小均导致计算精度降低;SIP算法中,计算过程中的迭代次数、运行时间与水头变化的收敛指标值呈明显的反比关系。当收敛指标设定小于一定值时,迭代次数和收敛指标的对数值呈现明显的线性关系。     

Modeling surface and ground water mixing in the hyporheic zone using MODFLOW and MT3D

We used a three-dimensional MODFLOW model, paired with MT3D, to simulate hyporheic zones around debris dams and meanders along a semi-arid stream. MT3D simulates both advective transport and sink/source mixing of solutes, in contrast to particle tracking (e.g. MODPATH), which only considers advection. We delineated the hydrochemically active hyporheic zone based on a new definition, specifically as near-stream subsurface zones receiving a minimum of 10% surface water within a 10-day travel time. Modeling results indicate that movement of surface water into the hyporheic zone is predominantly an advective process. We show that debris dams are a key driver of surface water into the subsurface along the experimental reach, causing the largest flux rates of water across the streambed and creating hyporheic zones with up to twice the cross-sectional area of other hyporheic zones. Hyporheic exchange was also found in highly sinuous segments of the experimental reach, but flux rates are lower and the cross-sectional areas of these zones are generally smaller. Our modeling approach simulated surface and ground water mixing in the hyporheic zone, and thus provides numerical approximations that are more comparable to field-based observations of surface–groundwater exchange than standard particle-tracking simulations.