Search USGS⌕ Search

Geology topics

Christian D. Langevin

Publications and source records attributed to Christian D. Langevin.

At least 37 records · Page 2Linked to original sources
USGS2014

MODFLOW-based coupled surface water routing and groundwater-flow simulation

In this paper, we present a flexible approach for simulating one- and two-dimensional routing of surface water using a numerical surface water routing (SWR) code implicitly coupled to the groundwater-flow process in MODFLOW. Surface water routing in SWR can be simulated using a diffusive-wave approximation of the Saint-Venant equations and/or a simplified level-pool approach. SWR can account for surface water flow controlled by backwater conditions caused by small water-surface gradients or surface water control structures. A number of typical surface water control structures, such as culverts, weirs, and gates, can be represented, and it is possible to implement operational rules to manage surface water stages and streamflow. The nonlinear system of surface water flow equations formulated in SWR is solved by using Newton methods and direct or iterative solvers. SWR was tested by simulating the (1) Lal axisymmetric overland flow, (2) V-catchment, and (3) modified Pinder-Sauer problems. Simulated results for these problems compare well with other published results and indicate that SWR provides accurate results for surface water-only and coupled surface water/groundwater problems. Results for an application of SWR and MODFLOW to the Snapper Creek area of Miami-Dade County, Florida, USA are also presented and demonstrate the value of coupled surface water and groundwater simulation in managed, low-relief coastal settings.

Groundwater↗
USGS2013

Estimation of capture zones and drawdown at the Northwest and West Well Fields, Miami-Dade County, Florida, using an unconstrained Monte Carlo analysis: recent (2004) and proposed conditions

Travel-time capture zones and drawdown for two production well fields, used for drinking-water supply in Miami-Dade County, southeastern Florida, were delineated by the U.S Geological Survey using an unconstrained Monte Carlo analysis. The well fields, designed to supply a combined total of approximately 250 million gallons of water per day, pump from the highly transmissive Biscayne aquifer in the urban corridor between the Everglades and Biscayne Bay. A transient groundwater flow model was developed and calibrated to field data to ensure an acceptable match between simulated and observed values for aquifer heads and net exchange of water between the aquifer and canals. Steady-state conditions were imposed on the transient model and a post-processing backward particle-tracking approach was implemented. Multiple stochastic realizations of horizontal hydraulic conductivity, conductance of canals, and effective porosity were simulated for steady-state conditions representative of dry, average and wet hydrologic conditions to calculate travel-time capture zones of potential source areas of the well fields. Quarry lakes, formed as a product of rock-mining activities, whose effects have previously not been considered in estimation of capture zones, were represented using high hydraulic-conductivity, high-porosity cells, with the bulk hydraulic conductivity of each cell calculated based on estimates of aquifer hydraulic conductivity, lake depths and aquifer thicknesses. A post-processing adjustment, based on calculated residence times using lake outflows and known lake volumes, was utilized to adjust particle endpoints to account for an estimate of residence-time-based mixing of lakes. Drawdown contours of 0.1 and 0.25 foot were delineated for the dry, average, and wet hydrologic conditions as well. In addition, 95-percent confidence intervals (CIs) were calculated for the capture zones and drawdown contours to delineate a zone of uncertainty about the median estimates. Results of the Monte Carlo simulations indicate particle travel distances at the Northwest Well Field (NWWF) and West Well Field (WWF) are greatest to the west, towards the Everglades. The man-made quarry lakes substantially affect particle travel distances. In general near the NWWF, the capture zones in areas with lakes were smaller in areal extent than capture zones in areas without lakes. It is possible that contamination could reach the well fields quickly, within 10 days in some cases, if it were introduced into lakes nearest to supply wells, with one of the lakes being only approximately 650 feet from the nearest supply well. In addition to estimating drawdown and travel-time capture zones of 10, 30, 100, and 210 days for the NWWF and the WWF under more recent conditions, two proposed scenarios were evaluated with Monte Carlo simulations: the potential hydrologic effects of proposed Everglades groundwater seepage mitigation and quarry-lake expansion. The seepage mitigation scenario included the addition of two proposed anthropogenic features to the model: (1) an impermeable horizontal flow barrier east of the L-31N canal along the western model boundary between the Everglades and the urban areas of Miami-Dade County, and (2) a recharge canal along the Dade-Broward Levee near the NWWF. Capture zones and drawdown for the WWF were substantially affected by the addition of the barrier, which eliminates flow from the western boundary into the active model domain, shifting the predominant capture zone source area from the west more to the north and south. The 95-percent CI for the 210-day capture zone moved slightly in the NWWF as a result of the recharge canal. The lake-expansion scenario incorporated a proposed increase in the number and surface area of lakes by an additional 25 square miles. This scenario represents a 150-percent increase from the 2004 lake surface area near both well fields, but with the majority of increase proposed near the NWWF. The lake-expansion scenario substantially decreased the extent of the 210-day capture zone of the NWWF, which is limited to the lakes nearest the well field under proposed conditions.

Florida↗
USGS2013

MODFLOW–USG version 1: An unstructured grid version of MODFLOW for simulating groundwater flow and tightly coupled processes using a control volume finite-difference formulation

A new version of MODFLOW, called MODFLOW–USG (for UnStructured Grid), was developed to support a wide variety of structured and unstructured grid types, including nested grids and grids based on prismatic triangles, rectangles, hexagons, and other cell shapes. Flexibility in grid design can be used to focus resolution along rivers and around wells, for example, or to subdiscretize individual layers to better represent hydrostratigraphic units. MODFLOW–USG is based on an underlying control volume finite difference (CVFD) formulation in which a cell can be connected to an arbitrary number of adjacent cells. To improve accuracy of the CVFD formulation for irregular grid-cell geometries or nested grids, a generalized Ghost Node Correction (GNC) Package was developed, which uses interpolated heads in the flow calculation between adjacent connected cells. MODFLOW–USG includes a Groundwater Flow (GWF) Process, based on the GWF Process in MODFLOW–2005, as well as a new Connected Linear Network (CLN) Process to simulate the effects of multi-node wells, karst conduits, and tile drains, for example. The CLN Process is tightly coupled with the GWF Process in that the equations from both processes are formulated into one matrix equation and solved simultaneously. This robustness results from using an unstructured grid with unstructured matrix storage and solution schemes. MODFLOW–USG also contains an optional Newton-Raphson formulation, based on the formulation in MODFLOW–NWT, for improving solution convergence and avoiding problems with the drying and rewetting of cells. Because the existing MODFLOW solvers were developed for structured and symmetric matrices, they were replaced with a new Sparse Matrix Solver (SMS) Package developed specifically for MODFLOW–USG. The SMS Package provides several methods for resolving nonlinearities and multiple symmetric and asymmetric linear solution schemes to solve the matrix arising from the flow equations and the Newton-Raphson formulation, respectively.

Techniques and Methods↗
USGS2013

Documentation of the seawater intrusion (SWI2) package for MODFLOW

The SWI2 Package is the latest release of the Seawater Intrusion (SWI) Package for MODFLOW. The SWI2 Package allows three-dimensional vertically integrated variable-density groundwater flow and seawater intrusion in coastal multiaquifer systems to be simulated using MODFLOW-2005. Vertically integrated variable-density groundwater flow is based on the Dupuit approximation in which an aquifer is vertically discretized into zones of differing densities, separated from each other by defined surfaces representing interfaces or density isosurfaces. The numerical approach used in the SWI2 Package does not account for diffusion and dispersion and should not be used where these processes are important. The resulting differential equations are equivalent in form to the groundwater flow equation for uniform-density flow. The approach implemented in the SWI2 Package allows density effects to be incorporated into MODFLOW-2005 through the addition of pseudo-source terms to the groundwater flow equation without the need to solve a separate advective-dispersive transport equation. Vertical and horizontal movement of defined density surfaces is calculated separately using a combination of fluxes calculated through solution of the groundwater flow equation and a simple tip and toe tracking algorithm. Use of the SWI2 Package in MODFLOW-2005 only requires the addition of a single additional input file and modification of boundary heads to freshwater heads referenced to the top of the aquifer. Fluid density within model layers can be represented using zones of constant density (stratified flow) or continuously varying density (piecewise linear in the vertical direction) in the SWI2 Package. The main advantage of using the SWI2 Package instead of variable-density groundwater flow and dispersive solute transport codes, such as SEAWAT and SUTRA, is that fewer model cells are required for simulations using the SWI2 Package because every aquifer can be represented by a single layer of cells. This reduction in number of required model cells and the elimination of the need to solve the advective-dispersive transport equation results in substantial model run-time savings, which can be large for regional aquifers. The accuracy and use of the SWI2 Package is demonstrated through comparison with existing exact solutions and numerical solutions with SEAWAT. Results for an unconfined aquifer are also presented to demonstrate application of the SWI2 Package to a large-scale regional problem.

Techniques and Methods↗
USGS2012

Effects of sea-level rise on salt water intrusion near a coastal well field in southeastern Florida

A variable-density groundwater flow and dispersive solute transport model was developed for the shallow coastal aquifer system near a municipal supply well field in southeastern Florida. The model was calibrated for a 105-year period (1900 to 2005). An analysis with the model suggests that well-field withdrawals were the dominant cause of salt water intrusion near the well field, and that historical sea-level rise, which is similar to lower-bound projections of future sea-level rise, exacerbated the extent of salt water intrusion. Average 2005 hydrologic conditions were used for 100-year sensitivity simulations aimed at quantifying the effect of projected rises in sea level on fresh coastal groundwater resources near the well field. Use of average 2005 hydrologic conditions and a constant sea level result in total dissolved solids (TDS) concentration of the well field exceeding drinking water standards after 70 years. When sea-level rise is included in the simulations, drinking water standards are exceeded 10 to 21 years earlier, depending on the specified rate of sea-level rise.

Florida↗
USGS2012

Modeling variably saturated subsurface solute transport with MODFLOW-UZF and MT3DMS

The MT3DMS groundwater solute transport model was modified to simulate solute transport in the unsaturated zone by incorporating the unsaturated-zone flow (UZF1) package developed for MODFLOW. The modified MT3DMS code uses a volume-averaged approach in which Lagrangian-based UZF1 fluid fluxes and storage changes are mapped onto a fixed grid. Referred to as UZF-MT3DMS, the linked model was tested against published benchmarks solved analytically as well as against other published codes, most frequently the U.S. Geological Survey's Variably-Saturated Two-Dimensional Flow and Transport Model. Results from a suite of test cases demonstrate that the modified code accurately simulates solute advection, dispersion, and reaction in the unsaturated zone. Two- and three-dimensional simulations also were investigated to ensure unsaturated-saturated zone interaction was simulated correctly. Because the UZF1 solution is analytical, large-scale flow and transport investigations can be performed free from the computational and data burdens required by numerical solutions to Richards' equation. Results demonstrate that significant simulation runtime savings can be achieved with UZF-MT3DMS, an important development when hundreds or thousands of model runs are required during parameter estimation and uncertainty analysis. Three-dimensional variably saturated flow and transport simulations revealed UZF-MT3DMS to have runtimes that are less than one tenth of the time required by models that rely on Richards' equation. Given its accuracy and efficiency, and the wide-spread use of both MODFLOW and MT3DMS, the added capability of unsaturated-zone transport in this familiar modeling framework stands to benefit a broad user-ship.

Ground Water↗
USGS2012

Future of groundwater modeling

With an increasing need to better manage water resources, the future of groundwater modeling is bright and exciting. However, while the past can be described and the present is known, the future of groundwater modeling, just like a groundwater model result, is highly uncertain and any prediction is probably not going to be entirely representative. Thus we acknowledge this as we present our vision of where groundwater modeling may be headed.

Ground Water↗
USGS2012

Improving sub-grid scale accuracy of boundary features in regional finite-difference models

As an alternative to grid refinement, the concept of a ghost node, which was developed for nested grid applications, has been extended towards improving sub-grid scale accuracy of flow to conduits, wells, rivers or other boundary features that interact with a finite-difference groundwater flow model. The formulation is presented for correcting the regular finite-difference groundwater flow equations for confined and unconfined cases, with or without Newton Raphson linearization of the nonlinearities, to include the Ghost Node Correction (GNC) for location displacement. The correction may be applied on the right-hand side vector for a symmetric finite-difference Picard implementation, or on the left-hand side matrix for an implicit but asymmetric implementation. The finite-difference matrix connectivity structure may be maintained for an implicit implementation by only selecting contributing nodes that are a part of the finite-difference connectivity. Proof of concept example problems are provided to demonstrate the improved accuracy that may be achieved through sub-grid scale corrections using the GNC schemes.

Advances in Water Resources↗
USGS2012

Simulated flow of groundwater and brine from a flooded salt mine in Livingston County, New York, and effects of remedial pumping on an overlying aquifer

Two ceiling collapses in the Retsof salt mine near Geneseo in upstate New York in spring 1994 resulted in the upward propagation of two columns of rubble through 600 feet of overlying shale and carbonate bedrock. This upward propagation formed a hydraulic connection between the lower confined aquifer (LCA) and the mine and allowed water from the aquifer and bedrock fracture zones that intersected the rubble columns to flow into the mine at a rate of 18,000 gallons per minute (gal/min) . All salt mining ceased in September 1995, and the mine was completely flooded by January 1996. The flow of water from the lower confined aquifer into the mine caused widespread drawdowns, and water levels in the aquifer declined by as much as 400 feet near the collapse area and by more than 50 feet at wells 7 miles to the north and south. Within 3 to 4 weeks of the collapses, water levels in about a dozen domestic and industrial wells had declined severely, and some wells went dry. Water levels in at least 58 wells in the lower and middle confined aquifers were affected by mine flooding. Groundwater in the upper unconfined aquifer and surface water in streams were unaffected by water-level drawdown, but channels of the Genesee River and Beards Creek were altered by land subsidence related to the mine collapse. Water levels recovered from 1996 through 2006, but the mine is now filled with about 15 billion gallons of saturated halite brine. The weight of the overlying rock and sediment is expected to cause the salt beds to deform and fill the mine cavity during the next several hundred years; this in turn could displace as much as 80 percent of the brine and cause it to move upward through the rubble chimneys, rendering the LCA unusable as a source of water supply. Saline water was detected in the LCA in 2002 but was found to be derived primarily from fractures in the limestone and shale units between the mine and the LCA, rather than from the mine. In September 2006, the mine company began a brine-mitigation project that entailed pumping five wells finished in limestone and shale units within the collapse areas to alter the flow gradient and thereby prevent further movement of brine and saline water into the LCA. The pumped brine was routed to an onsite desalination plant. At the same time, the U.S. Geological Survey (USGS) began a study in cooperation with the New York State Office of the Attorney General to construct numerical models to analyze the groundwater chemistry and delineate the directions of flow. Specific objectives of the study were to: * Assess the sources of salinity within the collapse area and identify the factors that control the movement and mixing of freshwater, saline waters from fracture zones, and brine; * Evaluate the likelihood that the pumping will induce anhydrite dissolution and lead to continued land subsidence; * Construct variable-density groundwater flow models to predict the effect of remedial pumping on salinity within the LCA; * Evaluate the effectiveness of remedial pumping in preventing the movement of saline water into the LCA; and * Predict the extent of brine migration 8 years after a hypothetical shutdown of all pumping in 2008. This report (1) summarizes the hydrogeologic setting and effects of mine flooding, (2) describes the geochemical and variable-density model simulations and their principal results, (3) discusses the implications of (a) continued pumping and desalination to protect the LCA and (b) a full shutdown of pumping after 2008, and (4) suggests further research that could lead to refinement of model predictions. Additional information may be found in Yager and others (2001 and 2009). These reports can be accessed at http://pubs.usgs.gov/pp/pp1611/ and http://pubs.usgs.gov/pp/pp1767/, respectively. A summary of simulation results can be accessed at http://ny.water.usgs.gov/projects/Coram/seawat/seawat.html.

New York↗
USGS2012

Documentation of the Surface-Water Routing (SWR1) Process for modeling surface-water flow with the U.S. Geological Survey Modular Ground-Water Model (MODFLOW-2005)

A flexible Surface-Water Routing (SWR1) Process that solves the continuity equation for one-dimensional and two-dimensional surface-water flow routing has been developed for the U.S. Geological Survey three-dimensional groundwater model, MODFLOW-2005. Simple level- and tilted-pool reservoir routing and a diffusive-wave approximation of the Saint-Venant equations have been implemented. Both methods can be implemented in the same model and the solution method can be simplified to represent constant-stage elements that are functionally equivalent to the standard MODFLOW River or Drain Package boundary conditions. A generic approach has been used to represent surface-water features (reaches) and allows implementation of a variety of geometric forms. One-dimensional geometric forms include rectangular, trapezoidal, and irregular cross section reaches to simulate one-dimensional surface-water features, such as canals and streams. Two-dimensional geometric forms include reaches defined using specified stage-volume-area-perimeter (SVAP) tables and reaches covering entire finite-difference grid cells to simulate two-dimensional surface-water features, such as wetlands and lakes. Specified SVAP tables can be used to represent reaches that are smaller than the finite-difference grid cell (for example, isolated lakes), or reaches that cannot be represented accurately using the defined top of the model. Specified lateral flows (which can represent point and distributed flows) and stage-dependent rainfall and evaporation can be applied to each reach. The SWR1 Process can be used with the MODFLOW Unsaturated Zone Flow (UZF1) Package to permit dynamic simulation of runoff from the land surface to specified reaches. Surface-water/groundwater interactions in the SWR1 Process are mathematically defined to be a function of the difference between simulated stages and groundwater levels, and the specific form of the reach conductance equation used in each reach. Conductance can be specified directly or calculated as a function of the simulated wetted perimeter and defined reach bed hydraulic properties, or as a weighted combination of both reach bed hydraulic properties and horizontal hydraulic conductivity. Each reach can be explicitly coupled to a single specific groundwater-model layer or coupled to multiple groundwater-model layers based on the reach geometry and groundwater-model layer elevations in the row and column containing the reach. Surface-water flow between reservoirs is simulated using control structures. Surface-water flow between reaches, simulated by the diffusive-wave approximation, can also be simulated using control structures. A variety of control structures have been included in the SWR1 Process and include (1) excess-volume structures, (2) uncontrolled-discharge structures, (3) pumps, (4) defined stage-discharge relations, (5) culverts, (6) fixed- or movable-crest weirs, and (7) fixed or operable gated spillways. Multiple control structures can be implemented in individual reaches and are treated as composite flow structures. Solution of the continuity equation at the reach-group scale (a single reach or a user-defined collection of individual reaches) is achieved using exact Newton methods with direct solution methods or exact and inexact Newton methods with Krylov sub-space methods. Newton methods have been used in the SWR1 Process because of their ability to solve nonlinear problems. Multiple SWR1 time steps can be simulated for each MODFLOW time step, and a simple adaptive time-step algorithm, based on user-specified rainfall, stage, flow, or convergence constraints, has been implemented to better resolve surface-water response. A simple linear- or sigmoid-depth scaling approach also has been implemented to account for increased bed roughness at small surface-water depths and to increase numerical stability. A line-search algorithm also has been included to improve the quality of the Newton-step upgrade vector, if possible. The SWR1 Process has been benchmarked against one- and two-dimensional numerical solutions from existing one- and two-dimensional numerical codes that solve the dynamic-wave approximation of the Saint-Venant equations. Two-dimensional solutions test the ability of the SWR1 Process to simulate the response of a surface-water system to (1) steady flow conditions for an inclined surface (solution of Manning's equation), and (2) transient inflow and rainfall for an inclined surface. The one-dimensional solution tests the ability of the SWR1 Process to simulate a looped network with multiple upstream inflows and several control structures. The SWR1 Process also has been compared to a level-pool reservoir solution. A synthetic test problem was developed to evaluate a number of different SWR1 solution options and simulate surface-water/groundwater interaction. The solution approach used in the SWR1 Process may not be applicable for all surface-water/groundwater problems. The SWR1 Process is best suited for modeling long-term changes (days to years) in surface-water and groundwater flow. Use of the SWR1 Process is not recommended for modeling the transient exchange of water between streams and aquifers when local and convective acceleration and other secondary effects (for example, wind and Coriolis forces) are substantial. Dam break evaluations and two-dimensional evaluations of spatially extensive domains are examples where acceleration terms and secondary effects would be significant, respectively.

Techniques and Methods↗
USGS2011

Effect of tidal fluctuations on transient dispersion of simulated contaminant concentrations in coastal aquifers

Variable-density groundwater models require extensive computational resources, particularly for simulations representing short-term hydrologic variability such as tidal fluctuations. Saltwater-intrusion models usually neglect tidal fluctuations and this may introduce errors in simulated concentrations. The effects of tides on simulated concentrations in a coastal aquifer were assessed. Three analyses are reported: in the first, simulations with and without tides were compared for three different dispersivity values. Tides do not significantly affect the transfer of a hypothetical contaminant into the ocean; however, the concentration difference between tidal and non-tidal simulations could be as much as 15%. In the second analysis, the dispersivity value for the model without tides was increased in a zone near the ocean boundary. By slightly increasing dispersivity in this zone, the maximum concentration difference between the simulations with and without tides was reduced to as low as 7%. In the last analysis, an apparent dispersivity value was calculated for each model cell using the simulated velocity variations from the model with tides. Use of apparent dispersivity values in models with a constant ocean boundary seems to provide a reasonable approach for approximating tidal effects in simulations where explicit representation of tidal fluctuations is not feasible.

Hydrogeology Journal↗
USGS2011

Predictive uncertainty analysis of a saltwater intrusion model using null-space Monte Carlo

Because of the extensive computational burden and perhaps a lack of awareness of existing methods, rigorous uncertainty analyses are rarely conducted for variable-density flow and transport models. For this reason, a recently developed null-space Monte Carlo (NSMC) method for quantifying prediction uncertainty was tested for a synthetic saltwater intrusion model patterned after the Henry problem. Saltwater intrusion caused by a reduction in fresh groundwater discharge was simulated for 1000 randomly generated hydraulic conductivity distributions, representing a mildly heterogeneous aquifer. From these 1000 simulations, the hydraulic conductivity distribution giving rise to the most extreme case of saltwater intrusion was selected and was assumed to represent the "true" system. Head and salinity values from this true model were then extracted and used as observations for subsequent model calibration. Random noise was added to the observations to approximate realistic field conditions. The NSMC method was used to calculate 1000 calibration-constrained parameter fields. If the dimensionality of the solution space was set appropriately, the estimated uncertainty range from the NSMC analysis encompassed the truth. Several variants of the method were implemented to investigate their effect on the efficiency of the NSMC method. Reducing the dimensionality of the null-space for the processing of the random parameter sets did not result in any significant gains in efficiency and compromised the ability of the NSMC method to encompass the true prediction value. The addition of intrapilot point heterogeneity to the NSMC process was also tested. According to a variogram comparison, this provided the same scale of heterogeneity that was used to generate the truth. However, incorporation of intrapilot point variability did not make a noticeable difference to the uncertainty of the prediction. With this higher level of heterogeneity, however, the computational burden of generating calibration-constrained parameter fields approximately doubled. Predictive uncertainty variance computed through the NSMC method was compared with that computed through linear analysis. The results were in good agreement, with the NSMC method estimate showing a slightly smaller range of prediction uncertainty than was calculated by the linear method. Copyright 2011 by the American Geophysical Union.

Water Resources Research↗
USGS2011

Use of hydrologic and hydrodynamic modeling for ecosystem restoration

Planning and implementation of unprecedented projects for restoring the greater Everglades ecosystem are underway and the hydrologic and hydrodynamic modeling of restoration alternatives has become essential for success of restoration efforts. In view of the complex nature of the South Florida water resources system, regional-scale (system-wide) hydrologic models have been developed and used extensively for the development of the Comprehensive Everglades Restoration Plan. In addition, numerous subregional-scale hydrologic and hydrodynamic models have been developed and are being used for evaluating project-scale water management plans associated with urban, agricultural, and inland costal ecosystems. The authors provide a comprehensive summary of models of all scales, as well as the next generation models under development to meet the future needs of ecosystem restoration efforts in South Florida. The multiagency efforts to develop and apply models have allowed the agencies to understand the complex hydrologic interactions, quantify appropriate performance measures, and use new technologies in simulation algorithms, software development, and GIS/database techniques to meet the future modeling needs of the ecosystem restoration programs.

Florida↗
USGS2010

Application of SEAWAT to select variable-density and viscosity problems

SEAWAT is a combined version of MODFLOW and MT3DMS, designed to simulate three-dimensional, variable-density, saturated groundwater flow. The most recent version of the SEAWAT program, SEAWAT Version 4 (or SEAWAT_V4), supports equations of state for fluid density and viscosity. In SEAWAT_V4, fluid density can be calculated as a function of one or more MT3DMS species, and optionally, fluid pressure. Fluid viscosity is calculated as a function of one or more MT3DMS species, and the program also includes additional functions for representing the dependence of fluid viscosity on temperature. This report documents testing of and experimentation with SEAWAT_V4 with six previously published problems that include various combinations of density-dependent flow due to temperature variations and/or concentration variations of one or more species. Some of the problems also include variations in viscosity that result from temperature differences in water and oil. Comparisons between the results of SEAWAT_V4 and other published results are generally consistent with one another, with minor differences considered acceptable.

Scientific Investigations Report↗
USGS2010

Creative use of pilot points to address site and regional scale heterogeneity in a variable-density model

Pilot points for parameter estimation were creatively used to address heterogeneity at both the well field and regional scales in a variable-density groundwater flow and solute transport model designed to test multiple hypotheses for upward migration of fresh effluent injected into a highly transmissive saline carbonate aquifer. Two sets of pilot points were used within in multiple model layers, with one set of inner pilot points (totaling 158) having high spatial density to represent hydraulic conductivity at the site, while a second set of outer points (totaling 36) of lower spatial density was used to represent hydraulic conductivity further from the site. Use of a lower spatial density outside the site allowed (1) the total number of pilot points to be reduced while maintaining flexibility to accommodate heterogeneity at different scales, and (2) development of a model with greater areal extent in order to simulate proper boundary conditions that have a limited effect on the area of interest. The parameters associated with the inner pilot points were log transformed hydraulic conductivity multipliers of the conductivity field obtained by interpolation from outer pilot points. The use of this dual inner-outer scale parameterization (with inner parameters constituting multipliers for outer parameters) allowed smooth transition of hydraulic conductivity from the site scale, where greater spatial variability of hydraulic properties exists, to the regional scale where less spatial variability was necessary for model calibration. While the model is highly parameterized to accommodate potential aquifer heterogeneity, the total number of pilot points is kept at a minimum to enable reasonable calibration run times.

Book↗
USGS2010

A comparison between SWI and SEAWAT: the importance of dispersion, inversion and vertical anisotropy

SWI and SEAWAT are both computer codes designed to model variable-density systems. One of the options in SWI is to model Dupuit interface flow, where freshwater and seawater are separated by an interface. In this paper we compare seawater intrusion model results of SWI to model results of SEAWAT, which simulates full variable-density flow and transport. Results indicate that SWI is valid for many variable-density systems. For the case considered in this paper, SWI results are accurate when the simulated width of the transition zone between seawater to freshwater is 15% or less of the scale of the problem, density inversion (saltwater over freshwater) occurs over only a small part of the model domain, and the ratio of vertical to horizontal hydraulic conductivity is larger than 0.01. Results also show that the simulated interface moves further inland using SWI than for the same conditions using SEAWAT. SWI is preferable to be used in systems where run times for a fully-coupled variable-density flow and transport model would be prohibitive; for the case considered here, SWI run times were a few seconds and SEAWAT run times were almost three hours.

Conference Paper↗
USGS2010

Effect of sea-level rise on future coastal groundwater resources in southern Florida, USA

An existing variable‐density groundwater flow and solute transport model, developed for the northern part of Broward County, Florida, was used to predict the effect of sealevel rise on future coastal groundwater resources. Using average annual conditions from 2005, simulations were performed for 100 years into the future using four different rates of sea‐level rise: 0, 24, 48, and 88 centimeters per century. Results from these predictive analyses suggest that the average concentration of groundwater withdrawn at the municipal well field will exceed the potable limit after 70, 60, 55, and 49 years, respectively, for the four simulations.

Florida↗