Search

Managed aquifer recharge has become a standard water resources management practice to promote the development of locally sustainable water supplies and combat water scarcity. However, installation of injection wells for replenishment purposes in urban areas with complex hydrogeology faces many challenges, such as limited land availability, potential impacts on municipal production wells and known subsurface contamination plumes, and complex spatially variable hydraulic connections between aquifer units. To assess the feasibility and cost‐effectiveness of injecting advanced treated water (ATW) into a complex urban aquifer system, a Simulation‐Optimization (SO) model was developed to automate a systematic search for the most cost‐effective locations to install new wells for injecting various quantities of ATW, if feasible. The generalized workflow presented here uses an existing MODFLOW groundwater model—along with advanced optimization routines that are publicly available—to flexibly accommodate a multiobjective function, complex constraints, and specific project requirements. The model successfully placed wells for injection of 1 to 4 MGD of ATW in aquifers underlying the study area. The injection well placement was primarily constrained by avoiding excessive impact on environmental sites with underlying groundwater plumes. The largest costs were for well installation and piping to the wells from the existing ATW pipes. This workflow is readily adaptable to other sites with different complexities, decision variables, or constraints. [ABSTRACT FROM AUTHOR]

MODFLOW 6 is the latest in a line of six "core" versions of MODFLOW released by the U.S. Geological Survey. The MODFLOW 6 architecture supports incorporation of additional hydrologic processes, in addition to groundwater flow, and allows interaction between processes. The architecture supports multiple model instances and multiple types of models within a single simulation, a flexible approach to formulating and solving the equations that represent hydrologic processes, and recent advances in interoperability, which allow MODFLOW to be accessed and controlled by external programs. The present version of MODFLOW 6 consolidates popular capabilities available in MODFLOW variants, such as the unstructured grid support in MODFLOW‐USG, the Newton‐Raphson formulation in MODFLOW‐NWT, and the support for partitioned stress boundaries in MODFLOW‐CDSS. The flexible multi‐model capability allows users to configure MODFLOW 6 simulations to represent the local‐grid refinement (LGR) capabilities available in MODFLOW‐LGR, the multi‐species transport capabilities in MT3DMS, and the coupled variable‐density capabilities available in SEAWAT. This paper provides a new, holistic and integrated overview of simulation capabilities made possible by the MODFLOW 6 architecture, and describes how ongoing and future development can take advantage of the program architecture to integrate new capabilities in a way that is minimally invasive and automatically compatible with the existing MODFLOW 6 code. [ABSTRACT FROM AUTHOR]

The sustainability of irrigated agriculture in many arid and semiarid areas of the world is at risk because of a combination of several interrelated factors, including lack of fresh water, lack of drainage, the presence of high water tables, and salinization of soil and groundwater resources. Nowhere in the United States are these issues more apparent than in the San Joaquin Valley of California. A solid understanding of salinization processes at regional spatial and decadal time scales is required to evaluate the sustainability of irrigated agriculture. A hydro-salinity model was developed to integrate subsurface hydrology with reactive salt transport for a 1,400-k[m.sup.2] study area in the San Joaquin Valley. The model was used to reconstruct historical changes in salt storage by irrigated agriculture over the past 60 years. We show that patterns in soil and groundwater salinity were caused by spatial variations in soil hydrology, the change from local groundwater to snowmelt water as the main irrigation water supply, and by occasional droughts. Gypsum dissolution was a critical component of the regional salt balance. Although results show that the total salt input and output were about equal for the past 20 years, the model also predicts salinization of the deeper aquifers, thereby questioning the sustainability of irrigated agriculture. regional hydrology | salinization | vadose zone

In this paper, we develop a numerical solution to equations derived in the first part of this study to simulate the thaw subsidence of an unsaturated frozen soil due to surface heating. Solutions obtained in terms of pore pressure, temperature, degree of unfrozen water content, ice content, and porosity provide the spatial and temporal variations of these quantities in a thawing soil column. A one-dimensional form of governing equations is discretized by employing a Crank-Nicolson finite-difference scheme. These highly nonlinear equations are then solved by an iterative procedure at every time step for predescribed positive surface temperatures that produce thawing. The numerical solutions are obtained for a 1-m silt column. Results show that, above the frozen region, unfrozen water saturation increases due to low hydraulic conductivity, which prevents the melting ice from draining, thus causing accumulation of water above it. It has also been observed that most of the settlement at the surface is due to thaw consolidation rather than settlement due to drainage. The former is completed in seven days of a 12-day simulation period. Furthermore, a sensitivity analysis was carried out to determine the response of the simulation to deviations in the values of certain parameters. It is found that the settlement estimations are very sensitive to the thaw-settlement parameter and phase-composition curve, thus requiring accurate determination of these parameters for a reliable simulation of the thaw-consolidation process.

Construction and maintenance of engineering structures in arctic regions demands an understanding of and the ability to cope with environmental problems produced by permafrost. Where thawing of the permafrost cannot be prevented, the range, amount, and extent of degradation are factors that must be taken into account in the design of engineering structures. Predictions of these factors are facilitated by the use of a mathematical model designed to predict the physical system under consideration. This study addresses the problem by developing a predictive model to describe permafrost thaw subsidence based upon a complete formulation of the problem. This problem is one of multiphase (air, ice, water, and solid) transport and deformation in a porous medium. The conservation of mass equations for liquid water, ice with phase change, and deforming soil solids are developed. Darcy's law is extended to include moisture movement due to thermal gradients. The quasi-static equilibrium equations and stress-strain relations that assume a perfectly elastic solid matrix are also employed. Simultaneous variations of moisture retention and phase-composition curves with temperature and pore pressure are also incorporated in the model. The energy-conservation equation, which includes terms due to viscous dissipation and compression effects, along with the other equations expressed, is found to satisfactorily represent the physical processes occurring in thawing permafrost soils. Approximations to simulate the thawing of a one-dimensional soil column are then developed from the general model.

The performance of time series models is assessed using synthetic head series simulated with a numerical model that solves Richards' equation for variably saturated flow. Heads were simulated in a homogeneous unconfined aquifer between two parallel canals; measured daily precipitation and potential evaporation are specified at the land surface and root water uptake is simulated. The head response to a precipitation event is nonlinear and depends on the saturation degree and rainfall before and after the precipitation event while evaporation reduction occurs during summers. Synthetic series were generated for 27 years and three different soil types; the unsaturated zone thickness varies between 0 and >5 m. The synthetic head series were simulated with a linear and nonlinear time series model. Performance of a linear time series model with four parameters, using a scaled Gamma response, gave R2 values ranging from 0.67 to 0.96. The nonlinear time series model with five parameters simulates recharge using a root zone reservoir after which the head response to recharge is simulated with a scaled Gamma response function. The nonlinear time series model was able to simulate all synthetic head series very well with R2 values above 0.9 for almost all models. The head response of the nonlinear model to a single precipitation event compares well to the response of the variably saturated groundwater model. The provided scripts may be used to simulate synthetic head series for other climates or for systems with additional complexity to assess the performance of other data‐driven models. [ABSTRACT FROM AUTHOR]

Abstract: A groundwater flow model is typically used to provide the flow field for conducting groundwater solute transport simulations. The advection term of the mass conserved formulation for groundwater transport assumes that the flow field is perfectly balanced and that all water flowing into a numerical grid cell is exactly balanced by outflows after accounting for sources/sinks or internal storage. However, in many complicated regional or site‐scale models, there may be localized flow balance errors that may be difficult to eliminate through tighter flow convergence tolerances due to simulation time constraints or numerical limits on convergence tolerances. Thus, if water is erroneously gained or lost within a grid cell during the flow computation, the solutes within it will also be numerically affected in the associated transport simulation. Transport solutions neglect this error in groundwater flow as the transport equations that are solved assume no error in flow. This flow imbalance error can however have consequences on the transport solution ranging from unnoticeable errors in the resulting concentrations to spurious oscillations that can grow in time and hinder further solution. An approach has been suggested here, to explicitly handle these flow imbalances during mass conserved advective transport computations and report them in the corresponding transport mass balance output, as corrections that are needed to handle errors originating in the flow solution. Example problems are provided to explain the concepts and demonstrate the impacts. [ABSTRACT FROM AUTHOR]

Infiltration along ephemeral channels plays an important role in groundwater recharge in arid regions. A model is presented for estimating spatial variability of seepage due to streambed heterogeneity along channels based on measurements of streamflow-front velocities in initially dry channels. The diffusion-wave approximation to the Saint-Venant equations, coupled with Philip's equation for infiltration, is connected to the groundwater model MODFLOW and is calibrated by adjusting the saturated hydraulic conductivity of the channel bed. The model is applied to portions of two large water delivery canals, which serve as proxies for natural ephemeral streams. Estimated seepage rates compare well with previously published values. Possible sources of error stem from uncertainty in Manning's roughness coefficients, soil hydraulic properties and channel geometry. Model performance would be most improved through more frequent longitudinal estimates of channel geometry and thalweg elevation, and with measurements of stream stage over time to constrain wave timing and shape. This model is a potentially valuable tool for estimating spatial variability in longitudinal seepage along intermittent and ephemeral channels over a wide range of bed slopes and the influence of seepage rates on groundwater levels. Copyright © 2012 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

Guha, Hillol and Sorab Panday, 2012. Impact of Sea Level Rise on Groundwater Salinity in a Coastal Community of South Florida. Journal of the American Water Resources Association (JAWRA) 48(3): 510-529. DOI: 10.1111/j.1752-1688.2011.00630.x Abstract: Freshwater resources of coastal communities in the United States and world over are threatened by the rate of sea level rise. According to recent estimates by various governmental agencies and climate researchers, the global sea level rise is likely to be between 0.6 and 2.1 m by the year 2100. South Florida is a coastal community and much of its coastline is subject to sea level rise and potential impacts to wetlands and the water resources of the area. To understand what the impact of sea level rise would cause to the groundwater level and salinity intrusion, an integrated groundwater and surface water model was developed for North Miami-Dade and Broward Counties of South Florida. The model was calibrated against daily groundwater heads, base flows in canals, and chloride concentrations for a period of one year and six months. Three separate sensitivity analyses were conducted by increasing the sea level by 0.6, 0.9, and 1.22 m. Results of the simulations shows increase of groundwater heads in some areas from 4 to 15%; whereas the average relative chloride concentrations increased significantly by 100-600% in some wells. The increase in groundwater elevations and chloride concentrations varies from location of the wells and its proximity to the coast. The model results indicate that even a 0.6 m increase in sea level (which is the conservative estimate) is likely to impair the vital freshwater resources in many parts of South Florida. [ABSTRACT FROM AUTHOR]