Search

Understanding the retention of per‐ and polyfluoroalkyl substances (PFAS) in the vadose zone is vital to the management of impacted sites. This paper examines PFAS retention in the unsaturated zone above the water table using a mathematical model, MODFLOW‐USG‐Transport PFAS or "USGT‐PFAS." The USGT‐PFAS model incorporates adsorption onto air‐water interfaces, providing a more comprehensive understanding of PFAS retention near the water table and release to groundwater. Modeling of a hypothetical perfluorooctane sulfonic acid (PFOS) site under various idealized site conditions illustrated that the impacts on PFOS retention from smallest to largest were water table fluctuations, low episodic recharge, constant recharge, moderate episodic recharge, constant recharge with water table fluctuations, and high episodic recharge. PFOS retention also varied by sand type, with greater retention occurring in simulations incorporating coarse sand with low capillary potential versus fine sand with high capillary potential. PFAS management strategies were also explored, including the adaptation of gas sparging, a method traditionally used for volatile organic compounds. Gas sparging can concentrate PFAS in groundwater and the vadose zone around the water table, facilitating retention or removal. Model simulations for a simplified hypothetical site demonstrated that PFAS can be substantially retained in the unsaturated zone once gas sparging results in an upward concentration of PFAS in groundwater and the unsaturated zone near the water table. Modeling can aid in understanding PFAS behavior but requires simulation of multiple interrelated processes to correctly predict PFAS fate and transport in subsurface conditions. [ABSTRACT FROM AUTHOR]

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]

Simulating the interaction of groundwater with surface water networks using traditional boundary packages available with MODFLOW‐USG can be challenging for complex systems. Often several package types are required as they are typically purpose built. Moreover, these packages generally do not interact with one another which complicates accounting of groundwater discharge at different points within the system. Here, we demonstrate that the connected linear network (CLN) package of MODFLOW‐USG, and advances therein in USG‐Transport, can be used to simulate groundwater interaction with a complex surface water network comprised of creeks, ponds, wetlands, and springs, in a manner that is comparable with these other packages, but with additional benefits, including explicit routing of water between the features. Article impact statement: We demonstrate the Connected Linear Network package can be used to simulate the exchange of groundwater with complex surface water networks. [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

The χMD matrix solver package is incorporated into USGS groundwater modeling software, such as MODFLOW-NWT, MODFLOW-USG, and MT3D. The solver is used to solve matrices assembled through numerical discretization of the groundwater flow equation, and solute transport equations. χMD has demonstrated its higher robustness, faster execution speed, and more efficient memory usage compared to the existing solvers for many types of groundwater flow problems. χMD uses preconditioned iterative Krylov-subspace methods and consists of preconditioning and acceleration modules. Because the solver package uses a variety of preconditioning features including level-based incomplete lower-upper (ILU) factorization method with a drop tolerance scheme, users must choose optimal preconditioning parameters to improve execution speed and robustness. In order to examine how the preconditioning parameters, ILU factorization level, and drop tolerance values affect the overall performance of the matrix solver, we evaluated five different groundwater model applications using MODFLOWUSG that include different numerical complexities. For those five cases, the number of discretization nodes varied from 10,000 cells to 730,300 cells. From the analysis, we found that the preconditioning parameters greatly affect execution times and memory usage of the preconditioning and acceleration procedures. In addition, a combination of the ILU level between five to seven and the drop tolerance value between 10-2 and 10-3 usually resulted in shorter overall execution time. Our study suggests that the users can elicit higher performance and robustness of the χMD matrix solver using this combination of the parameters and enhance computational efficiency of solving groundwater and solute transport problems. [ABSTRACT FROM AUTHOR]

Density‐dependent flow and transport solutions for coastal saltwater intrusion investigations, analyses of fluid injection into deep brines, and studies of convective fingering and instabilities of denser fluids moving through less dense fluids typically formulate the groundwater flow equation in terms of pressure or equivalent freshwater head. A formulation of the flow equation in terms of hydraulic head is presented here as an alternative. The hydraulic‐head formulation can facilitate adaptation of existing constant‐density groundwater flow codes to include density‐driven flow by avoiding the need to convert between freshwater head and hydraulic head within the code and by incorporating density‐dependent terms as a compartmentalized "correction" to constant‐density calculations already performed by the code. The hydraulic‐head formulation also accommodates complexities such as unconfined groundwater flow and Newton‐Raphson solution schemes more readily than the freshwater‐head formulation. Simulation results are presented for four example problems solved using an implementation of the hydraulic‐head formulation in MODFLOW. [ABSTRACT FROM AUTHOR]

Axisymmetric groundwater models are used for simulating radially symmetric conditions. Groundwater simulators built specifically to model axisymmetric conditions are most commonly used for simulating aquifer tests. Although some numerical models capable of simulating flow and solute transport that are developed in the cartesian coordinate system framework offer flexibility to simulate axisymmetric conditions, most of the numerical groundwater models, such as the MODFLOW family of codes, are based on structured grids in which axisymmetric flows cannot be directly simulated. Researchers in the past have provided methods to manipulate aquifer properties to mimic axisymmetric conditions. This study presents a methodology that takes advantage of the unstructured grids of MODFLOW‐USG to simulate axisymmetric models within the MODFLOW framework. To develop axisymmetric models, the intercell interface area arrays of MODFLOW‐USG were calculated to accurately represent coaxial cylindrical model cells. Three examples are presented to demonstrate the application of MODFLOW‐USG for axisymmetric modeling: a pumping well with delayed yield effects, a vadose zone flow model simulating an infiltration basin, and a density‐dependent saltwater intrusion problem for a circular island. Results were verified against analytical solutions and published numerical codes. Article impact statement: This study presents a new methodology to create axisymmetric models using MODFLOW‐USG. [ABSTRACT FROM AUTHOR]

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]

A dynamic subtime-stepping method is described for solving contaminant transport problems that utilize higher-order implicit time-marching procedures with higher-accuracy nonoscillating spatial-discretization methods to resolve sharp-plume fronts in advection-dominated systems. Nonoscillating spatial-discretization methods for the advective term prevent unphysical oscillations and minimize numerical diffusion. Second-order temporal accuracy is achieved by using the Crank-Nicholson implicit scheme, however nonoscillating properties may be violated for Courant Numbers larger than one, leading to spurious oscillations. The proposed subtiming method allows use of small time-step sizes in critical portions of a domain, with larger time-step sizes in other locations. This locally limits the Courant number where required and still keeps the general solution free of time-step size restrictions. This technique makes it possible to apply higher-order nonoscillating schemes with higher-order temporal weighting for advection-dominated flows, even when the Courant number is much greater than one. Feasibility and applicability of the dynamic implicit subtiming method are demonstrated through three proof-of-concept example problems. [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]

The article discusses the uncertainty on the future of groundwater modeling. It explores works on the future of hydrology, hydrogeology, and trends in groundwater modeling software development. It asserts that the future of groundwater modeling will depend largely on the types of problems relative to water supply and water quality. It also notes the use of models in connection with optimization methods to better manage resources.

Various approaches have been proposed to manage the nonlinearities associated with the unconfined flow equation and to simulate perched groundwater conditions using the MODFLOW family of codes. The approaches comprise a variety of numerical techniques to prevent dry cells from becoming inactive and to achieve a stable solution focused on formulations of the unconfined, partially-saturated, groundwater flow equation. Keeping dry cells active avoids a discontinuous head solution which in turn improves the effectiveness of parameter estimation software that relies on continuous derivatives. Most approaches implement an upstream weighting of intercell conductance and Newton-Raphson linearization to obtain robust convergence. In this study, several published approaches were implemented in a stepwise manner into MODFLOW for comparative analysis. First, a comparative analysis of the methods is presented using synthetic examples that create convergence issues or difficulty in handling perched conditions with the more common dry-cell simulation capabilities of MODFLOW. Next, a field-scale three-dimensional simulation is presented to examine the stability and performance of the discussed approaches in larger, practical, simulation settings. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Dissolved compounds in the water phase are rejected during freezing of soils. This phenomenon occurs due to the inability of ice crystals to incorporate most alien molecules. The ice grows only by association with water molecules, and molecules of solutes are rejected into the unfrozen water. A mathematical model for salt rejection in freezing saline soils has been developed by employing the mass balance equations for water and solute in a saturated porous medium, and the conservation of energy equation. Additional relations introduced include a moisture retention curve as a function of ice-water capillary pressure, the Clapeyron equation and expressions for the heat capacity and the effective thermal conductivity of the soil. The resulting set of seven equations in terms of solute concentration, soil temperature, pore water, ice, and capillary pressures, and degrees of water and ice saturation have been solved simultaneously using a Newton-Raphson linearization, with implicit iterative treatment of nonlinearities. Results indicate the development of unfrozen high salt concentration regions trapped in a frozen zone at sufficiently high freezing rates. The numerical results are compared with experimental data obtained by some researchers. Favorable match is obtained between the theoretical and experimental results. Finally, extension of the model to situations where ice lensing and soil deformations are significant is discussed. [ABSTRACT FROM AUTHOR]

Groundwater pumping is a boundary condition in groundwater models, and surface application, reuse, diversions, and water reclamation plant discharges are boundary conditions in hydrologic or hydraulic (i.e., surface) models. These processes can be simulated, however, as internal flow and solute transport components within a fully integrated groundwater-surface water model. A methodology has been developed and incorporated into the MODHMS fully integrated surface-subsurface modeling framework to enable evaluation of water distribution and application and impacts to water quality on surface and subsurface domains. This methodology allows the supply and delivery of water and associated water quality to be a part of the model solution rather than being predetermined boundary conditions to the system. Furthermore, this methodology handles temporal and spatial changes to land use and water quality and associated effects on outdoor water distribution, infiltration, runoff, and evapotranspiration properties that are related to surface cover or land use classifications. This new methodology was applied to a portion of the Upper Santa Clara River system and underlying groundwater basins of southern California to aid in the evaluation of a Cl- total maximum daily load (TMDL) study. The outcome of this modeling evaluation has led to an enhanced understanding of how Cl- is transported throughout the surface-subsurface model domain and improved collaboration among stakeholders. The model serves as a powerful tool with which to develop and evaluate various approaches for the Cl- TMDL implementation for the Upper Santa Clara River system. [ABSTRACT FROM AUTHOR]

A diverse group of problems requires quantification of the entire hydrologic cycle by the integrated simulation of water flow in the surface and subsurface regimes. In a transient integrated simulation of the water cycle, the time step size is a key factor in controlling the solution accuracy and the simulation efficiency for a given spatial discretization. In general, if the time step size is sufficiently small, the resulting solution becomes more accurate but with higher computational cost. Thus, to maintain an acceptable level of solution accuracy in the entire simulation domain, the time step size is restricted by the relatively rapid responses in the surface flow regime. As the relatively rapid responses are typically limited to a small portion of the surface domain compared with the groundwater system, a large portion of the domain tends to be temporally overdiscretized. The implicit subtime stepping approach described here can apply smaller subtime steps only to the subdomain where the accuracy requirements are needed. In this work, generalized formulations for implicit subtime stepping in the numerical solution of the nonlinear coupled surface-subsurface equations were derived and implemented into the integrated model HydroGeoSphere. Application to several problems showed that implicit subtime stepping can significantly improve the simulation efficiency with minimal loss in accuracy. The methodology was successfully applied to enhance the computational efficiency of an integrated flow simulation in the San Joaquin Valley, California, where the characteristic response time near surface drainage streams is orders of magnitude shorter than in the groundwater regime. [ABSTRACT FROM AUTHOR]

MODFLOW SURFACT is a state-of-the-art simulator that utilizes vadose zone flow and transport equations to provide practical solutions to the analysis of flow and contaminant transport at various levels of complexity and sophistication as needed for site evaluation and closure. The variably saturated flow equation can be solved with standard retention functions or with bimodal or multimodal relative permeability curves for unsaturated flow in porous and fractured systems. The equation can further be solved with pseudo-soil retention functions for confined--unconfined simulations and for use in wellbore hydraulics. Finally, the equation can be cast in terms of air phase flow to analyze subsurface air flow behavior. The variably saturated transport equation can be solved for an unsaturated medium or can be used for confined--unconfined situations. The passive phase of flow can be included in the equation to include both air and water phases in the transport situation. An immobile multicomponent nonaqueous phase liquid (NAPL) phase can further be included in the transport simulation with equilibrium partitioning providing mass transfer between phases, which adjusts NAPL saturations. Dual domain equations can be condensed into the transport equation to provide capabilities for analyzing transport in fractured media. General reaction capabilities provide analyses of complex environmental and geochemical interactions. Two examples are provided to demonstrate the value of a comprehensive simulation capability for site investigations. [ABSTRACT FROM AUTHOR]