The role of recharge and evapotranspiration as hydraulic drivers of ion concentrations in shallow groundwater on Everglades tree islands,Florida (USA)

Recently, evapotranspiration has been hypothesized to promote the secondary formation of calcium carbonate year‐round on tree islands in the Everglades by influencing groundwater ions concentrations. However, the role of recharge and evapotranspiration as drivers of shallow groundwater ion accumulation has not been investigated. The goal of this study is to develop a hydrologic model that predicts the chloride concentrations of shallow tree island groundwater and to determine the influence of overlying biomass and underlying geologic material on these concentrations. Groundwater and surface water levels and chloride concentrations were monitored on eight constructed tree islands at the Loxahatchee Impoundment Landscape Assessment (LILA) from 2007 to 2010. The tree islands at LILA were constructed predominately of peat, or of peat and limestone, and were planted with saplings of native tree species in 2006 and 2007. The model predicted low shallow groundwater chloride concentrations when inputs of regional groundwater and evapotranspiration‐to‐recharge rates were elevated, while low evapotranspiration‐to‐recharge rates resulted in a substantial increase of the chloride concentrations of the shallow groundwater. Modeling results indicated that evapotranspiration typically exceeded recharge on the older tree islands and those with a limestone lithology, which resulted in greater inputs of regional groundwater. A sensitivity analysis indicated the shallow groundwater chloride concentrations were most sensitive to alterations in specific yield during the wet season and hydraulic conductivity in the dry season. In conclusion, the inputs of rainfall, underlying hydrologic properties of tree islands sediments and forest structure may explain the variation in ion concentration seen across Everglades tree islands. Copyright © 2012 John Wiley & Sons, Ltd.     

Smoothed particle hydrodynamics pore-scale simulations of unstable immiscible flow in porous media

We have conducted a series of high-resolution numerical experiments using the Pair-Wise Force Smoothed Particle Hydrodynamics (PF-SPH) multiphase flow model. First, we derived analytical expressions relating parameters in the PF-SPH model to the surface tension and static contact angle. Next, we used the model to study viscous fingering, capillary fingering, and stable displacement of immiscible fluids in porous media for a wide range of capillary numbers and viscosity ratios. We demonstrated that the steady state saturation profiles and the boundaries of viscous fingering, capillary fingering, and stable displacement regions compare favorably with micromodel laboratory experimental results. For a displacing fluid with low viscosity, we observed that the displacement pattern changes from viscous fingering to stable displacement with increasing injection rate. When a high viscosity fluid is injected, transition behavior from capillary fingering to stable displacement occurred as the flow rate was increased. These observations are also in agreement with the results of the micromodel laboratory experiments.     

Flow path oscillations in transient ground-water simulations of large peatland systems

Transient numerical simulations of the Glacial Lake Agassiz Peatland near the Red Lakes in Northern Minnesota were constructed to evaluate observed reversals in vertical ground-water flow. Seasonal weather changes were introduced to a ground-water flow model by varying evapotranspiration and recharge over time. Vertical hydraulic reversals, driven by changes in recharge and evapotranspiration were produced in the simulated peat layer. These simulations indicate that the high specific storage associated with the peat is an important control on hydraulic reversals. Seasonally driven vertical flow is on the order of centimeters in the deep peat, suggesting that seasonal vertical advective fluxes are not significant and that ground-water flow into the deep peat likely occurs on decadal or longer time scales. Particles tracked within the ground-water flow model oscillate over time, suggesting that seasonal flow reversals will enhance vertical mixing in the peat column. The amplitude of flow path oscillations increased with increasing peat storativity, with amplitudes of about 5 cm occurring when peat specific storativity was set to about 0.05 m⁻¹.     

Modelling solute leaching during fingered flow by integrating and expanding various theoretical and empirical concepts

Abstract

Wetting front instability (fingered flow) accelerates solute transport through the unsaturated zone to the groundwater table. Whether fingers widen or dissipate close to the groundwater is unclear. Water flow in a two-dimensional artificial capillary fringe below a dry layer exhibiting fingered flow was investigated. The flow diverged strongly in the wet soil, suggesting that fingers dissipate. Expressions for the finger radius in dry and wet soil were combined and adapted to a soil hydraulic property parameterization popular in numerical modelling. The modified equation provided finger radii for soils in humid and arid climates. The fingers in the arid soil were excessively wide. The finger radii were used to model solute transport, assuming fingers dissipated in the subsoil. Modelling was cumbersome for the arid climate. One shower may often be insufficient to trigger fingering in arid regions with short, heavy showers. In soils with shallow groundwater, the diverging subsoil flow determines solute leaching.     

Thermohaline convection in flowing groundwater

Geothermal activity creates destabilising temperature gradients which are significant in some aquifers. Usually, in such aquifers stabilising salinity gradients also exist. The combination of temperature and salinity distribution in the aquifer may induce various types of hydrodynamic instabilities which were identified in a previous article. The present article concerns the effect of anisotropic characteristics of the hydrodynamic dispersion on the growth of instabilities in the aquifer. Three different mechanisms may lead to instability of the flow field: (a) buoyancy forces may induce convection currents; if the difference between the convection velocity of salt, due to the hydraulic gradients, and that of heat is negligible, then this mechanism is generally most effective in planes parallel to the hydraulic velocity of the fluid (velocity due to the hydraulic gradient); (b) the difference between heat and salt effective diffusivities may lead to overstability; this mechanism is most effective in planes perpendicular to the hydraulic velocity; (c) the difference between the convection velocity of salt and that of heat may induce oscillations which are most effective in planes parallel to the hydraulic velocity. The growth of instabilities in an aquifer of unlimited length is different from their growth in an aquifer of limited length. In the latter thermohaline convection develops in planes perpendicular to the hydraulic velocity, whereas in the former it develops in planes forming an angle θ with the hydraulic gradient. The development of convection cells in the flow field is identified by numerical experiments. These experiments identify the convection cell length and the angle formed between the thermohaline convection plane and the hydraulic gradient.     

Numerical modeling of geothermal groundwater flow in karst aquifer system in eastern Weibei,Shaanxi Province,China

The quantitative assessment of geothermal water resources is important to the exploitation and utilization of geothermal resources. In the geothermal water systems the density of groundwater changes with the temperature, therefore the variations in hydraulic heads and temperatures are very complicated. A three-dimensional density-dependent model coupling the groundwater flow and heat transport is established and used to simulate the geothermal water flow in the karst aquifers in eastern Weibei, Shaanxi Province, China. The multilayered karst aquifer system in the study area is cut by some major faults which control the regional groundwater flow. In order to calibrate and simulate the effect of the major faults, each fault is discretized as a belt of elements with special hydrological parameters in the numerical model. The groundwater dating data are used to be integrated with the groundwater flow pattern and calibrate the model. Simulation results show that the calculated hydraulic heads and temperature fit with the observed data well.

     

Numerical modeling of geothermal groundwater flow in karst aquifer system in eastern Weibei, Shaanxi Province, China

The quantitative assessment of geothermal water resources is important to the exploitation and utilization of geothermal resources. In the geothermal water systems the density of groundwater changes with the temperature, therefore the variations in hydraulic heads and temperatures are very complicated. A three-dimensional density-dependent model coupling the groundwater flow and heat transport is established and used to simulate the geothermal water flow in the karst aquifers in eastern Weibei, Shaanxi Province, China. The multilayered karst aquifer system in the study area is cut by some major faults which control the regional groundwater flow. In order to calibrate and simulate the effect of the major faults, each fault is discretized as a belt of elements with special hydrological parameters in the numerical model. The groundwater dating data are used to be integrated with the groundwater flow pattern and calibrate the model. Simulation results show that the calculated hydraulic heads and temperature fit with the observed data well.     

Climatic and topographic limits to the abundance of bog pools

On patterned peatlands, open water pools develop within a matrix of terrestrial vegetation (‘ridges’). Regional patterns in the distribution of ridge–pool complexes suggest that the relative cover of these two surface types is controlled in part by climate wetness, but landscape topography must also be an important controlling factor. In this paper, a functional model that relates relative cover of ridges and pools to climate and surface gradient was developed and tested. The model was formulated in terms of a water budget, based on the differential effects of ridges and pools on losses by evapotranspiration and subsurface flow. It predicts a positive relationship between surface gradient and ridge proportion, with a linear effect related to water supply and ridge hydraulic conductivity, modified at high ridge proportion by differences in evapotranspiration between ridges and pools. The limit to patterned peatland distribution occurs where the surface is completely covered by ridges. The model may be sensitive or insensitive to climate differences between localities, depending on whether hydraulic characteristics of ridge peat co‐vary with water supply. To distinguish between these alternative hypotheses, surface gradient and ridge proportion were surveyed along 20 transects in each of three localities in Scotland that differ threefold in net precipitation to pools. The results of the field survey served to reject the climate‐sensitive hypothesis, but were consistent with the climate‐insensitive hypothesis. Analysis of the residuals suggested that variation within localities was related more to topographic control of water supply than to ridge hydraulic conductivity or developmental stage. Hence, within this maritime climate region, the distribution of ridge–pool complexes and the relative abundance of pools are controlled mainly by topographic variables. Field surveys across both maritime and continental regions are required to confirm a subtle climatic effect that allows pools to occur on higher gradients in drier climates than in wetter climates. Further development and testing of the functional model will provide a stronger basis for assessing potential feedback between climate change, peatland surface structure and methane emission from pools. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical modeling of geothermal groundwater flow in karst aquifer system in eastern Weibei,Shaanxi Province,China

The quantitative assessment of geothermal water resources is important to the exploitation and utilization of geothermal resources. In the geothermal water systems the density of groundwater changes with the temperature, therefore the variations in hydraulic heads and temperatures are very complicated. A three-dimensional density-dependent model coupling the groundwater flow and heat transport is established and used to simulate the geothermal water flow in the karst aquifers in eastern Weibei, Shaanxi Province, China. The multilayered karst aquifer system in the study area is cut by some major faults which control the regional groundwater flow. In order to calibrate and simulate the effect of the major faults, each fault is discretized as a belt of elements with special hydrological parameters in the numerical model. The groundwater dating data are used to be integrated with the groundwater flow pattern and calibrate the model. Simulation results show that the calculated hydraulic heads and temperature fit with the observed data well.     

Identification of hydraulic conductivity distributions in density dependent flow fields of submarine groundwater discharge modeling using adjoint-state sensitivities

A mathematical optimal control method is developed to identify a hydraulic conductivity distribution in a density dependent flow field. Using a variational method, the adjoint partial differential equations are obtained for the density- dependent state equations used for the saline aquifer water flow. The adjoint equations are numerically solved in through a finite difference method. The developed method is applied to identify the hydraulic conductivity distribution through the numerical solution of an optimal control problem. To demonstrate the effectiveness of the optimal control method, three numerical experiments are conducted with artificial observation data. The results indicate that the developed method has the potential to accurately identify the hydraulic conductivity distribution in a saline water aquifer flow system.     

The effect of density-driven flow on the transport of high-concentration solutes in the hyporheic zone

In this study, both laboratory experiments and numerical simulations were conducted to investigate the effect of density-driven flow on the transport of high-concentration pollutants in the hyporheic zone. The results show that the density gradient can change the flow of pore water and the strong density-driven flow can lead to an unstable flow, which increases the effect of preferential flow and thus causes the appearance of solute fingers in the hyporheic zone. Notably, these solute fingers become more obvious with the increase of depth. The appearance of solute fingers depends on the relative strength of the pumping exchange and density gradient, which are represented by the dimensionless number M* and N* respectively. Finger flows appear near the interface when M^* is less than 0.5 N^*. This study may contribute to better understanding the transport and destination of solutes and thus may provide some insights into the assessment on pollution incidents.     

Steady state flow with hydraulic conductivity change around large diameter wells

Darcian flow law in aquifers assumes that the aquifer hydraulic conductivity is constant and the groundwater movement is due only to the piezometric level changes through hydraulic gradient. In practice, after the well development the aquifer just around the well has comparatively larger hydraulic conductivity and gradient. Patchy aquifer solutions in the literature consider sudden hydraulic conductivity changes with distance for the steady state flow. The change of transmissivity is demonstrated by the application of slope‐matching procedure to actual field data. It is the main purpose of this paper to derive simple analytical expressions for aquifer parameter evaluations with steadily decreasing hydraulic conductivity around the well. Spatial nonlinear hydraulic conductivity changes around a large‐diameter well within the depression cone of a confined aquifer are considered as exponentially decreasing functions of the radial distance. Copyright © 2010 John Wiley & Sons, Ltd.     

Radionuclide Transport Behavior in a Generic Geological Radioactive Waste Repository

We performed numerical simulations of groundwater flow and radionuclide transport to study the influence of several factors, including the ambient hydraulic gradient, groundwater pressure anomalies, and the properties of the excavation damaged zone (EDZ), on the prevailing transport mechanism (i.e., advection or molecular diffusion) in a generic nuclear waste repository within a clay‐rich geological formation. By comparing simulation results, we show that the EDZ plays a major role as a preferential flowpath for radionuclide transport. When the EDZ is not taken into account, transport is dominated by molecular diffusion in almost the totality of the simulated domain, and transport velocity is about 40% slower. Modeling results also show that a reduction in hydraulic gradient leads to a greater predominance of diffusive transport, slowing down radionuclide transport by about 30% with respect to a scenario assuming a unit gradient. In addition, inward flow caused by negative pressure anomalies in the clay‐rich formation further reduces transport velocity, enhancing the ability of the geological barrier to contain the radioactive waste. On the other hand, local high gradients associated with positive pressure anomalies can speed up radionuclide transport with respect to steady‐state flow systems having the same regional hydraulic gradients. Transport behavior was also found to be sensitive to both geometrical and hydrogeological parameters of the EDZ. Results from this work can provide useful knowledge toward correctly assessing the post‐closure safety of a geological disposal system.     

Numerical simulation of groundwater flow patterns using flux as upper boundary

Since the 1960s, most of the studies on groundwater flow systems by analytical and numerical modelling have been based on given‐head upper boundaries. The disadvantage of the given‐head approach is that the recharge into and discharge from a basin vary with changes in hydraulic conductivity and/or basin geometry. Consequently, flow patterns simulated with given‐head boundaries but with different hydraulic conductivities and/or basin geometry may not reflect the effects of these variables. We conducted, therefore, numerical simulations of groundwater flow in theoretical drainage basins using flux as the upper boundary and realistically positioned fluid‐potential sinks while changing the infiltration intensity, hydraulic conductivities, and geometric configuration of the basin. The simulated results demonstrate that these variables are dominant factors controlling the flow pattern in a laterally closed drainage basin. The ratio of infiltration intensity to hydraulic conductivity (R_ic) has been shown to be an integrated pattern‐parameter in a basin with a given geometric configuration and possible fluid‐potential‐sink distribution. Successively, the changes in flow patterns induced by stepwise reductions in R_ic are identical, regardless of whether the reductions are due to a decrease in infiltration intensity or an increase in hydraulic conductivity. The calculated examples show five sequential flow patterns containing (i) only local, (ii) local–intermediate, (iii) local–intermediate–regional, (iv) local–regional, and (v) just regional flow systems. The R_ic was found to determine also whether a particular sink is active or not as a site of discharge. Flux upper boundary is preferable for numerical simulation when discussing the flow patterns affected by a change of infiltration, the hydraulic conductivity, or the geometry of a basin. Copyright © 2012 John Wiley & Sons, Ltd.     

Application of ERT,Saline Tracer and Numerical Studies to Delineate Preferential Paths in Fractured Granites

Accurate quantification of in situ heterogeneity and flow processes through fractured geologic media remains elusive for hydrogeologists due to the complexity in fracture characterization and its multiscale behavior. In this research, we demonstrated the efficacy of tracer-electrical resistivity tomography (ERT) experiments combined with numerical simulations to characterize heterogeneity and delineate preferential flow paths in a fractured granite aquifer. A series of natural gradient saline tracer experiments were conducted from a depth window of 18 to 22 m in an injection well (IW) located inside the Indian Institute of Technology Hyderabad campus. Tracer migration was monitored in a time-lapse mode using two cross-sectional surface ERT profiles placed in the direction of flow gradient. ERT data quality was improved by considering stacking, reciprocal measurements, resolution indicators, and geophysical logs. Dynamic changes in subsurface electrical properties inferred via resistivity anomalies were used to highlight preferential flow paths of the study area. Temporal changes in electrical resistivity and tracer concentration were monitored along the vertical in an observation well located at 48 m to the east of the IW. ERT-derived tracer breakthrough curves were in agreement with geochemical sample measurements. Fracture geometry and hydraulic properties derived from ERT and pumping tests were further used to evaluate two mathematical conceptualizations that are relevant to fractured aquifers. Results of numerical analysis conclude that dual continuum model that combines matrix and fracture systems through a flow exchange term has outperformed equivalent continuum model in reproducing tracer concentrations at the monitoring wells (evident by a decrease in RMSE from 199 to 65 mg/L). A sensitivity analysis on model simulations conclude that spatial variability in hydraulic conductivity, local-scale dispersion, and flow exchange at fracture-matrix interface have a profound effect on model simulations.     

Transverse structure of tidal and residual flow and sediment concentration in estuaries

An analytical and a numerical model are used to understand the response of velocity and sediment distributions over Gaussian-shaped estuarine cross-sections to changes in tidal forcing and water depth. The estuaries considered here are characterized by strong mixing and a relatively weak along-channel density gradient. It is also examined under what conditions the fast, two-dimensional analytical flow model yields results that agree with those obtained with the more complex three-dimensional numerical model. The analytical model reproduces and explains the main velocity and sediment characteristics in large parts of the parameter space considered (average tidal velocity amplitude, 0.1–1 m s^− 1 and maximum water depth, 10–60 m). Its skills are lower for along-channel residual flows if nonlinearities are moderate to high (strong tides in deep estuaries) and for transverse flows and residual sediment concentrations if the Ekman number is small (weak tides in deep estuaries). An important new aspect of the analytical model is the incorporation of tidal variations in the across-channel density gradient, causing a double circulation pattern in the transverse flow during slack tides. The gradient also leads to a new tidally rectified residual flow component via net advection of along-channel tidal momentum by the density-induced transverse tidal flow. The component features landward currents in the channel and seaward currents over the slopes and is particularly effective in deeper water. It acts jointly with components induced by horizontal density differences, Coriolis-induced tidal rectification and Stokes discharge, resulting in different along-channel residual flow regimes. The residual across-channel density gradient is crucial for the residual transverse circulation and for the residual sediment concentration. The clockwise density-induced circulation traps sediment in the fresher water over the left slope (looking up-estuary in the northern hemisphere). Model results are largely consistent with available field data of well-mixed estuaries.     

LE TRITIUM DANS LES PRÉCIPITATIONS DE L'ASIE DU SUD-EST ET LE PHÉNOMÉNE DE LA MOUSSON EN 1966 ET EN 1967

Abstract

The problem of non-steady flow of water in a soil-plant system can be described by adding a sink term to the continuity equation for soil water flow. In this paper the sink term is defined in two different ways. Firstly it is considered to be dependent on the hydraulic conductivity of the soil, on the difference in pressure head between the soil and the root-soil interface and some root effectiveness function. Secondly the sink is taken to be a prescribed function of the soil water content. The partial differential equation applying to the first problem is solved by both a finite difference (FD 1) and a finite element (FE 1) technique, that applying to the second problem by a finite difference approach (FD 2). The purpose of this paper is to verify the numerical models against field measurements, to compare the results obtained by the three numerical methods and to show how the finite element method can be applied to complex but realistic two-dimensional flow situations. Two examples are given. The first concerns one-dimensional flow and it compares numerical results with those obtained experimentally in the field from water balance studies on red cabbage (Brassica oleracea L. ‘Rode Herfst’) grown on a clay soil in the presence of a water table. The second example describes two-dimensional flow in a complex field situation in the Netherlands where flow takes place under cropped field conditions through five anisotropic layers. Water is supplied to the system by infiltration from two unlined ditches and is withdrawn from the system by evapotranspiration and by leakage to an underlying pumped aquifer.     

Fluid dispersion effects on density-driven thermohaline flow and transport in porous media

This study introduces the dispersive fluid flux of total fluid mass to the density-driven flow equation to improve thermohaline modeling of salt and heat transports in porous media. The dispersive fluid flux in the flow equation is derived to account for an additional fluid flux driven by the density gradient and mechanical dispersion. The coupled flow, salt transport and heat transport governing equations are numerically solved by a fully implicit finite difference method to investigate solution changes due to the dispersive fluid flux. The numerical solutions are verified by the Henry problem and the thermal Elder problem under a moderate density effect and by the brine Elder problem under a strong density effect. It is found that increment of the maximum ratio of the dispersive fluid flux to the advective fluid flux results in increasing dispersivity for the Henry problem and the brine Elder problem. The effects of the dispersive fluid flux on salt and heat transports under high density differences and high dispersivities are more noticeable than under low density differences and low dispersivities. Values of quantitative indicators such as the Nusselt number, mass flux, salt mass stored and maximum penetration depth in the brine Elder problem show noticeable changes by the dispersive fluid flux. In the thermohaline Elder problem, the dispersive fluid flux shows a considerable effect on the shape and the number of developed fingers and makes either an upwelling or a downwelling flow in the center of the domain. In conclusion, for the general case that involves strong density-driven flow and transport modeling in porous media, the dispersive fluid flux should be considered in the flow equation.     

Hydraulic Gradient Comparison Method to Estimate Aquifer Hydraulic Parameters Under Steady-State Conditions

The hydraulic gradient comparison method is an inverse method for estimation of aquifer hydraulic conductivity (or trans-missivity) and boundary conductance for a ground water flow model under steady-state conditions. This method, following formal optimization techniques, defines its objective function to minimize differences between interpreted (observed) and simulated hydraulic gradients, which results in minimization of differences between observed and simulated hydraulic heads. The key features of this method are that (1) the derived optimality conditions have an explicit form with a clear hydrology concept that is con-sistent with Darcy's law, and (2) the derived optimality conditions are spatially independent as they are a function of only local hydraulic conductivity and local hydraulic gradient. This second feature allows a multidimensional optimization problem to be solved by many one-dimensional optimization procedures simultaneously, which results in a substantial reduction in computation time. The results of the numerical performance testing on a heterogeneous hypothetical case confirm that minimizing gradient residuals in the entire model domain leads to minimizing head residuals. Application of the method in real-world projects requires rigorous conceptual model development, use of a global calibration target, and an iterative calibration proess. The conceptual model development includes interpretation of a potentiometric surface and estimation of other hydrologic parameters. This method has been applied to a wide range of real-world modeling projects, including the Rocky Mountain Arsenal and Rocky Flats sites in Colorado, which demonstrates that the method is efficient and practical.