Search

• Stream-function-based (SFB) method is developed for 2D particle tracking along exact streamlines. • The branch cut effect of multivariable stream function around wells is avoided. • Unlike existing methods, the SFB method does not lead to accumulative errors in the particle position. Particle tracking is an essential process in analyzing groundwater flow and solute transport in aquifers. A drawback of existing particle tracking methods, which are based on approximate estimations of the velocity in a time step, is the accumulative errors of numerous time steps would lead to deviation of particles away from the exact streamlines, especially when pumping and injecting wells exist. In this study, to improve the accuracy of particle tracking for the two-dimensional (2D) steady-state flow around wells, we propose a stream-function-based (SFB) particle tracking method in which a particle is traced along the exact streamline with a given value of the stream function, and the branch cut effect caused by the multivalued stream function around a well is eliminated. Three synthetic cases are used to demonstrate the advantages of the SFB method over the existing Pollock method and the fourth-order Runge-Kutta method. In the cases of a single pumping well and an injecting-pumping well pair, the Pollock and Runge-Kutta methods yield pathways that may be far away from the exact streamline, but the SFB method performs much better. In the case with a pumping well and two injecting wells, the SFB method successfully delineated a narrow capture zone of the pumping well, which requires a high accuracy particle tracking process. Moreover, the accuracy of particle positions obtained by the SFB method is insensitive to the step-travel distance. The study would be useful for accurate investigation of flow fields in aquifers disturbed by wells. [ABSTRACT FROM AUTHOR]

The Ordos Basin is one of the most intensively studied groundwater basins in China. The Ordos Plateau, located in the north part of the Ordos Basin, is ideal to study the pattern of regional groundwater circulation induced by water-table undulations due to the wavy topography and the relatively simple aquifer systems with macroscopically homogeneous sandstone. In catchments located near the first-order divide, the water table is found to be a subdued replica of the topography, and the nonclosed water-table contours in topographic highs of a catchment are indicative of regional groundwater outflow to other catchments. In topographic lows, groundwater-fed lakes/rivers, topography-driven flowing wells, water-loving and/or salt-tolerant vegetation, and soap holes are all indicative of discharge areas. In discharge areas, although groundwater inflow from recharge areas is relatively stable, seasonal variations in groundwater recharge and evapotranspiration lead to significant seasonal fluctuations in the water table, which can be used to estimate groundwater inflow and evapotranspiration rates based on water balance at different stages of water-table change. In the lowest reaches of a complex basin, superposition of local flow systems on regional flow systems has been identified based on groundwater samples collected from wells with different depths and geophysical measurements of apparent resistivity, both of which can be used for characterizing groundwater flow systems. This study enhances understanding of the pattern of regional groundwater circulation in the Ordos Plateau, and also tests the  effectiveness of methods for groundwater flow-system characterization. [ABSTRACT FROM AUTHOR]

The Badain Jaran Desert (BJD) in China is a desert with impressive sand dunes and a groundwater situation that has attracted numerous researchers. This paper gives an overview of the mysteries of groundwater in the BJD that are exhibited as five key problems identified in previous studies. These problems relate to the origin of the groundwater, the hydrological connection between the BJD and the Heihe River Basin (HRB), the infiltration recharge, the lake-groundwater interactions, and the features of stable isotope analyses. The existing controversial analyses and hypotheses have caused debate and have hindered effective water resources management in the region. In recent years, these problems have been partly addressed by additional surveys. It has been revealed that the Quaternary sandy sediments and Neogene-Cretaceous sandstones form a thick aquifer system in the BJD. Groundwater flow at the regional scale is dominated by a significant difference in water levels between the surrounding mountains and lowlands at the western and northern edges. Discharge of groundwater from the BJD to the downstream HRB occurs according to the regional flow. Seasonal fluctuations of the water level in lakes are less than 0.5 m due to the quasi-steady groundwater discharge. The magnitude of infiltration recharge is still highly uncertain because significant limitations existed in previous studies. The evaporation effect may be the key to interpreting the anomalous negative deuterium-excess in the BJD groundwater. Further investigations are expected to reveal the hydrogeological conditions in more detail. [ABSTRACT FROM AUTHOR]

Abstract: The occurrence of flowing wells in basins has been found to be closely related to the discharge area with an upward hydraulic gradient. Unfortunately, previous studies on upward gradient induced wellbore flow with equaling total inflow (Qin) in the deep and total outflow (Qout) in the shallow could not explain the occurrence of flowing wells. By representing wells using the MNW2 Package imbedded in MODFLOW 2005, we obtain the exchange of groundwater between the aquifer and the well in the discharge area of 3D unit basins and identify three scenarios: Qin = Qout, Qin > Qout > 0 and Qin > Qout = 0. The relationship of Qin > Qout well explains why flowing wells only develop in a limited part of the discharge area. Sensitivity analysis shows that well location, water table undulation, and basin length–depth ratio do not change the profile of the ratio of cumulative flow rate in a flowing well to total inflow (Qv/Qin) versus the relative elevation in the inflow segment, zin*, but could significantly change the length of the inflow segment; well depth could change both the length of the inflow segment and the profile of Qv/Qin versus zin*. Based on the numerical results in homogeneous and isotropic basins with different basin length–depth ratios in the current study, the ratio of inflow in the lower half part of a flowing well to the total inflow is found to be at least 67% and could be close to 100%, indicating that water at the outlets of flowing wells with long open sections is mainly from the deep part of the well. [ABSTRACT FROM AUTHOR]

Abstract: A new solution of transient confined–unconfined flow driven by a pumping well is developed and compared to previous approximate solutions of Moench and Prickett [Moench AF, Prickett TA. Radial flow in an infinite aquifer undergoing conversion from artesian to water table conditions. Water Resour Res 1972;8:494–9] and Hu and Chen [Hu L, Chen C. Analytical methods for transient flow to a well in a confined–unconfined aquifer. Ground Water 2008;46(4):642–6]. The problem is rewritten in dimensionless form with the Boltzmann transform. The nonlinear equation for flow in the unconfined zone is solved with the Runge–Kutta method. Position of the conversion interface is determined with an iteration scheme. This study shows that the confined–unconfined flow depends on three dimensionless parameters that represent the confined–unconfined storativity ratio (a D ), the ratio of the initial hydraulic head over the aquifer thickness (f i ), and the dimensionless pumping rate (q D ). The rate of expansion of the unconfined zone increases with q D , but decreases with a D and f i . Differences between the two previous approximate solutions and the new solution of this study are observable in the estimated position of the conversion interface and the drawdown–time curves. The new solution can be applied to estimate the time for confined–unconfined conversion to occur (critical conversion time), and the time when the pumping well becomes dry (critical drying time). The critical conversion time is found to be very sensitive to the initial hydraulic head. The critical drying time is often much larger than the critical conversion time and may never be observed during a finite pumping period. [Copyright &y& Elsevier]

• We summarize field observations of 14C age showing obvious GW aging in two basins. • Pumping causes not only rejuvenation, but also aging of GW in layered basins. • Although GW rejuvenation is larger in area, aging is greater in magnitude. • We identify two new mechanisms of pumping-induced GW aging in layered basins. • Topography-driven flow is the main cause of GW aging in the discharge area. As an indicator of the renewability of a groundwater basin, the spatial distribution of groundwater age has attracted intensive research. As reported in the literature, 3H measurements in 74 US aquifers showed pumping-induced groundwater rejuvenation, whereas 14C ages in two heavily developed aquifers in France and China indicated pumping-induced groundwater aging. Unfortunately, the mechanisms leading to the seemingly opposite behavior of pumping-induced aging and rejuvenation remain unclear. In this study, we conducted numerical simulations of groundwater age in layered aquifer-aquitard systems under different pumping scenarios and aquitard permeability to show that pumping-induced groundwater rejuvenation and aging always occur simultaneously. Groundwater rejuvenation widely occurs throughout the basin wherever the well is located, whereas groundwater aging mainly occurs to the downgradient parts of the basin. Although groundwater rejuvenation is larger in area, groundwater aging is greater in magnitude, and the magnitude of aging increases with pumping intensity and duration. Moreover, when pumping occurs in the discharge area, aging could occur in almost the entire aquifer to the downgradient of the pumping well, and the magnitude of aging is larger. The large area and large magnitude of aging in the discharge area is a direct result of upward leakage of old groundwater from underlying aquitards/aquifers, which is controlled by the topography-driven groundwater flow. By demonstrating the importance of topography-driven groundwater flow and underlying formations in groundwater age redistribution, this study has implications for determining the bottom of a groundwater flow and solute transport model and interpretation of various tracers. [ABSTRACT FROM AUTHOR]

• We separate the roles of upward discharge and aquitards in flowing artesian zone. • Basinal flow are coupled with well hydraulics to examine hydraulics of flowing wells. • A low K aquitard causes a large flowing artesian zone but low stable discharge rate. • Baseflow reduction is the main source of groundwater derived in flowing wells. Causes and water sources of flowing artesian wells attracted the interest of many hydrogeologists throughout history, however, a quantitative model that satisfactorily considers the roles of topography, groundwater recharge/discharge and aquitards on hydraulics of flowing wells is still lacking. In this study, a three-layer river-valley basin with a recharge boundary is used to obtain the basinal flow field, spatial distribution and transient discharge rates of flowing wells. Even if there is an aquitard separating the unconfined and the confined aquifers, groundwater discharge to the river still plays a critical role in occurrence of flowing wells. An aquitard with lower permeability would enlarge the zone with flowing wells, implying that flowing wells are more likely to occur in basins with continuous aquitards. By comparing flowpaths and inter-aquifer leakages before and after numerically installing a flowing well in a three-layer basin, we find decreased upward leakage in the discharge limb plays a much larger role than increased downward leakage in the recharge limb in the stable discharge rate of a flowing well. This recognition is different from previous studies where the stable discharge rate is supported only by downward leakage from the overlying aquifer. By considering well hydraulics in a basinal background flow field, this study improves understanding of the mechanisms of flowing artesian wells and the interaction of surface water and deep groundwater, indicating that flowing wells are valuable for sampling deep groundwater representative of deep and old groundwater in the discharge area of a basin. [ABSTRACT FROM AUTHOR]

• Three alternative conceptual models were constructed. • Each model was calibrated automatically and evaluated statistically. • Hydrogeological interpretations are the basis for the conceptualisation. • Modelling objectives determine the complexity of the conceptual models. Three alternative groundwater flow models were evaluated for Beijing Plain, China. The first model (AM1) was constructed with the "thin layer approach" in which all 9 model layers, including five aquifers separated by four aquitards, are continuously present in the same model area. The second model (AM2) was constructed with the "quasi-3D approach" in which the hydrogeological formations were classified into five aquifer units consisting of mixed permeable and semi-permeable layers at different depth ranges. The third model (AM3) was constructed with the "true layer approach" in which aquifers and aquitards were defined according to hydrostratigraphic properties, and model layers are absent in the area where corresponding hydrogeological formations intersect bedrocks. All 3 models were calibrated with the parameter optimization method under the steady state flow condition with the same hydrological stresses and observation data. All three models fit to observations well with the similar calibration criteria values. Furthermore, AIC and BIC information criteria could not distinguish three alternative models. Only KIC could identify AM3 as the best model. Major differences of the three alternative models were identified from a hydrogeological perspective. The AM1 model depicted an illusion through contour maps that groundwater was present everywhere in the deep aquifers. The model computed larger vertical leakages because more abstraction rates were assigned improperly in deep aquifers. The AM2 model was able to compute regional groundwater balances and depicted spatial groundwater level variations. However, the AM2 model computed longer groundwater travel times around the wellfield and should not be used for the delineation of the well field protection zones and contaminant transport simulation. The AM3 model could not only compute the regional groundwater balances and describe spatial groundwater distribution in deep confined aquifers, but also delineate the capture zone of the wellfield. It can therefore be used for simulating contaminant transport. Furthermore, the AM3 model is suitable to construct a coupled regional and local flow model for simulating a managed aquifer recharge scheme. [ABSTRACT FROM AUTHOR]