Geology topics

Joseph D. Hughes

Publications and source records attributed to Joseph D. Hughes.

At least 37 records · Page 2Linked to original sources

USGS2016

Scripting MODFLOW model development using Python and FloPy

Graphical user interfaces (GUIs) are commonly used to construct and postprocess numerical groundwater flow and transport models. Scripting model development with the programming language Python is presented here as an alternative approach. One advantage of Python is that there are many packages available to facilitate the model development process, including packages for plotting, array manipulation, optimization, and data analysis. For MODFLOW-based models, the FloPy package was developed by the authors to construct model input files, run the model, and read and plot simulation results. Use of Python with the available scientific packages and FloPy facilitates data exploration, alternative model evaluations, and model analyses that can be difficult to perform with GUIs. Furthermore, Python scripts are a complete, transparent, and repeatable record of the modeling process. The approach is introduced with a simple FloPy example to create and postprocess a MODFLOW model. A more complicated capture-fraction analysis with a real-world model is presented to demonstrate the types of analyses that can be performed using Python and FloPy.

Groundwater↗

USGS2016

Potential effects of alterations to the hydrologic system on the distribution of salinity in the Biscayne aquifer in Broward County, Florida

To address concerns about the effects of water-resource management practices and rising sea level on saltwater intrusion, the U.S. Geological Survey in cooperation with the Broward County Environmental Planning and Community Resilience Division, initiated a study to examine causes of saltwater intrusion and predict the effects of future alterations to the hydrologic system on salinity distribution in eastern Broward County, Florida. A three-dimensional, variable-density solute-transport model was calibrated to conditions from 1970 to 2012, the period for which data are most complete and reliable, and was used to simulate historical conditions from 1950 to 2012. These types of models are typically difficult to calibrate by matching to observed groundwater salinities because of spatial variability in aquifer properties that are unknown, and natural and anthropogenic processes that are complex and unknown; therefore, the primary goal was to reproduce major trends and locally generalized distributions of salinity in the Biscayne aquifer. The methods used in this study are relatively new, and results will provide transferable techniques for protecting groundwater resources and maximizing groundwater availability in coastal areas. The model was used to (1) evaluate the sensitivity of the salinity distribution in groundwater to sea-level rise and groundwater pumping, and (2) simulate the potential effects of increases in pumping, variable rates of sea-level rise, movement of a salinity control structure, and use of drainage recharge wells on the future distribution of salinity in the aquifer. Results from the simulation of historical conditions indicate that the model generally represents the observed greater westward extent of elevated salinity in the central part of the intruded area relative to the northern and southernmost parts of the intruded area. Results of sensitivity testing indicate that the extent of elevated salinity is most sensitive to pumping in areas where the source of saltwater is largely offshore, from the Atlantic Ocean, and is most sensitive to sea-level rise in areas where the source of salinity is downward leakage of brackish water from canals. Simulations of future scenarios indicate that increases in pumping near the existing interface may cause the interface to advance and decreases in pumping may cause it to retreat. Climatic effects, such as periods of prolonged drought or high precipitation, may augment or counteract long-term effects of changes in pumping on aquifer salinity at well fields. With increasing rates of sea-level rise, the freshwater-saltwater interface advances progressively inland, and flow-averaged salinities at well fields near the existing interface increase commensurately. Hypothetical southeastward (downstream) re-positioning of the existing G–54 salinity-control structure may prevent the interface from moving northwestward along and near the North New River canal, but beneficial effects are localized. Implementation of freshwater recharge wells in the city of Hallandale Beach may also have only a localized freshening effect in the aquifer and little appreciable effect on the freshwater-saltwater interface or on concentrations of salinity at well fields. Model accuracy and use are limited by uncertainty in the physical properties and boundary conditions of the system, uncertainty in historical and future conditions, and generalizations made in the mathematical relationships used to describe the physical processes of groundwater flow and transport. Because of these limitations, model results should be considered in relative rather than absolute terms. Nonetheless, model results do provide useful information on the relative scale of response of the system to changes in pumping distribution, sea-level rise, and mitigation activities.

Florida↗

USGS2015

Examples of deformation-dependent flow simulations of conjunctive use with MF-OWHM

The dependency of surface- and groundwater flows and aquifer hydraulic properties on deformation induced by changes in aquifer head is not accounted for in the standard version of MODFLOW. A new USGS integrated hydrologic model, MODFLOW-OWHM, incorporates this dependency by linking subsidence and mesh deformation with changes in aquifer transmissivity and storage coefficient, and with flows that also depend on aquifer characteristics and land-surface geometry. This new deformation-dependent approach is being used for the further development of the integrated Central Valley hydrologic model (CVHM) in California. Preliminary results from this application and from hypothetical test cases of similar systems show that changes in canal flows, stream seepage, and evapotranspiration from groundwater (ET gw ) are sensitive to deformation. Deformation feedback has been shown to also have an indirect effect on conjunctive surface- and groundwater use components with increased stream seepage and streamflows influencing surface-water deliveries and return flows. In the Central Valley model, land subsidence may significantly degrade the ability of the major canals to deliver surface water from the Delta to the San Joaquin and Tulare basins. Subsidence can also affect irrigation demand and ET gw , which, along with altered surface-water supplies, causes a feedback response resulting in changed estimates of groundwater pumping for irrigation. This modeling feature also may improve the impact assessment of dewatering-induced land subsidence/uplift (following irrigation pumping or coal-seam gas extraction) on surface receptors, inter-basin transfers, and surface infrastructure integrity.

Proceedings of the International Association of Hy↗

USGS2015

Coupling geophysical investigation with hydrothermal modeling to constrain the enthalpy classification of a potential geothermal resource.

An appreciable challenge in volcanology and geothermal resource development is to understand the relationships between volcanic systems and low-enthalpy geothermal resources. The enthalpy of an undeveloped geothermal resource in the Karckar region of Armenia is investigated by coupling geophysical and hydrothermal modeling. The results of 3-dimensional inversion of gravity data provide key inputs into a hydrothermal circulation model of the system and associated hot springs, which is used to evaluate possible geothermal system configurations. Hydraulic and thermal properties are specified using maximum a priori estimates. Limited constraints provided by temperature data collected from an existing down-gradient borehole indicate that the geothermal system can most likely be classified as low-enthalpy and liquid dominated. We find the heat source for the system is likely cooling quartz monzonite intrusions in the shallow subsurface and that meteoric recharge in the pull-apart basin circulates to depth, rises along basin-bounding faults and discharges at the hot springs. While other combinations of subsurface properties and geothermal system configurations may fit the temperature distribution equally well, we demonstrate that the low-enthalpy system is reasonably explained based largely on interpretation of surface geophysical data and relatively simple models.

Karckar Region↗

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↗

USGS2014

Quantifying the predictive consequences of model error with linear subspace analysis

All computer models are simplified and imperfect simulators of complex natural systems. The discrepancy arising from simplification induces bias in model predictions, which may be amplified by the process of model calibration. This paper presents a new method to identify and quantify the predictive consequences of calibrating a simplified computer model. The method is based on linear theory, and it scales efficiently to the large numbers of parameters and observations characteristic of groundwater and petroleum reservoir models. The method is applied to a range of predictions made with a synthetic integrated surface-water/groundwater model with thousands of parameters. Several different observation processing strategies and parameterization/regularization approaches are examined in detail, including use of the Karhunen-Loève parameter transformation. Predictive bias arising from model error is shown to be prediction specific and often invisible to the modeler. The amount of calibration-induced bias is influenced by several factors, including how expert knowledge is applied in the design of parameterization schemes, the number of parameters adjusted during calibration, how observations and model-generated counterparts are processed, and the level of fit with observations achieved through calibration. Failure to properly implement any of these factors in a prediction-specific manner may increase the potential for predictive bias in ways that are not visible to the calibration and uncertainty analysis process.

Water Resources Research↗

USGS2014

Feedback of land subsidence on the movement and conjunctive use of water resources

The dependency of surface- or groundwater flows and aquifer hydraulic properties on dewatering-induced layer deformation is not available in the USGS's groundwater model MODFLOW. A new integrated hydrologic model, MODFLOW-OWHM, formulates this dependency by coupling mesh deformation with aquifer transmissivity and storage and by linking land subsidence/uplift with deformation-dependent flows that also depend on aquifer head and other flow terms. In a test example, flows most affected were stream seepage and evapotranspiration from groundwater (ET gw ). Deformation feedback also had an indirect effect on conjunctive surface- and groundwater use components: Changed stream seepage and streamflows influenced surface-water deliveries and returnflows. Changed ET gw affected irrigation demand, which jointly with altered surface-water supplies resulted in changed supplemental groundwater requirements and pumping and changed return runoff. This modeling feature will improve the impact assessment of dewatering-induced land subsidence/uplift (following irrigation pumping or coal-seam gas extraction) on surface receptors, inter-basin transfers, and surface-infrastructure integrity.

Environmental Modelling and Software↗

USGS2014

Hydrologic conditions in urban Miami-Dade County, Florida, and the effect of groundwater pumpage and increased sea level on canal leakage and regional groundwater flow

The extensive and highly managed surface-water system in southeastern Florida constructed during the 20th Century has allowed for the westward expansion of urban and agricultural activities in Miami-Dade County. In urban areas of the county, the surface-water system is used to (1) control urban flooding, (2) supply recharge to production well fields, and (3) control seawater intrusion. Previous studies in Miami-Dade County have determined that on a local scale, leakage from canals adjacent to well fields can supply a large percentage (46 to 78 percent) of the total groundwater pumpage from production well fields. Canals in the urban areas also receive seepage from the Biscayne aquifer that is derived from a combination of local rainfall and groundwater flow from Water Conservation Area 3 and Everglades National Park, which are west of urban areas of Miami-Dade County. To evaluate the effects of groundwater pumpage on canal leakage and regional groundwater flow, the U.S. Geological Survey (USGS) developed and calibrated a coupled surface-water/groundwater model of the urban areas of Miami-Dade County, Florida. The model was calibrated by using observation data collected from January 1997 through December 2004. The model calibration was verified using observation data collected from January 2005 through December 2010. A 1-year warmup period (January 1996 through December 1996) was added prior to the start of the calibration period to reduce the effects of inaccurate initial conditions on model calibration. The model is designed to simulate surface-water stage and discharge in the managed canal system and dynamic canal leakage to the Biscayne aquifer as well as seepage to the canal from the aquifer. The model was developed using USGS MODFLOW–NWT with the Surface-Water Routing (SWR1) Process to simulate surface-water stage, surface-water discharge, and surface-water/groundwater interaction and the Seawater Intrusion (SWI2) Package to simulate seawater intrusion, respectively. Automated parameter estimation software (PEST) and highly parameterized inversion techniques were used to calibrate the model to observed surface-water stage, surface-water discharge, net surface-water subbasin discharge, and groundwater level data from 1997 through 2004 by modifying hydraulic conductivity, specific storage coefficients, specific yield, evapotranspiration parameters, canal roughness coefficients (Manning’s n values), and canal leakance coefficients. Tikhonov regularization was used to produce parameter distributions that provide an acceptable fit between model outputs and observation data, while simultaneously minimizing deviations from preferred values based on field measurements and expert knowledge. Analytical and simulated water budgets for the period from 1996 through 2010 indicate that most of the water discharging through the salinity control structures is derived from within the urban parts of the study area and that, on average, the canals are draining the Biscayne aquifer. Simulated groundwater discharge from the urban areas to the coast is approximately 7 percent of the total surface-water inflow to Biscayne Bay and is consistent with previous estimates of fresh groundwater discharge to Biscayne Bay. Simulated groundwater budgets indicate that groundwater pumpage in some surface-water basins ranges from 13 to 27 percent of the sum of local sources of groundwater inflow. The largest percentage of groundwater pumpage to local sources of groundwater inflow occurs in the basins that have the highest pumping rates (C–2 and C–100 Basins). The ratio of groundwater pumpage to simulated local sources of groundwater inflow is less than values calculated in previous local-scale studies. The position of the freshwater-seawater interface at the base of the Biscayne aquifer did not change notably during the simulation period (1996–2010), consistent with the similar positions of the interface in 1984, 1995, and 2011 under similar hydrologic and groundwater pumping conditions. Landward movement of the freshwater-seawater interface above the base of the aquifer is more prone to occur during relatively dry years. The model was used to evaluate the effect of increased groundwater pumpage and (or) increased sea level on canal leakage, regional groundwater flow, and the position of the freshwater-seawater interface. Permitted groundwater pumping rates, which generally exceed historical groundwater pumping rates, were used for Miami-Dade County Water and Sewer Department groundwater pumping wells in the base-case future scenario. Base-case future and increased pumping scenario results suggest seawater intrusion may occur at the Miami-Springs well field if the Miami Springs, Hialeah, and Preston well fields are operated using current permitted groundwater pumping rates. Scenario simulations also show that, in general, the canal system limits the adverse effects of proposed groundwater pumpage increases on water-level changes and saltwater intrusion. Proposed increases (up to a 7 percent increase) in groundwater pumpage do not have a notable effect on movement of the freshwater-seawater interface. Increased groundwater pumpage increased lateral groundwater inflow into basins subject to additional groundwater pumpage; however, most (55 percent) of the additional groundwater extracted from pumping wells was supplied by changes in canal seepage and leakage in urban areas of the model. Increased sea level caused increased water-table elevations in urban areas and decreased hydraulic gradients across the system; the largest increases in water-table elevations occurred seaward of the salinity control structures. The extent of flood-prone areas and the percentage of time water-table elevations in flood-prone areas were less than 0.5 foot below land surface increased with increased sea level. Increased sea level also resulted in landward migration of the freshwater-seawater interface; the largest changes in the position of the interface occurred seaward of the salinity control structures except in parts of the model area that were inundated by increased sea level. Decreased water-table gradients reduced groundwater inflow, groundwater outflow, canal exchanges, surface-water inflow, and surface-water outflow through salinity control structures. Results for the scenario that evaluated the combination of increased groundwater pumpage and increased sea level did not differ substantially from the scenario that evaluated increased sea level alone. Groundwater inflow, groundwater outflow, and canal exchanges were reduced in urban areas of the study area as a result of decreased water-table gradients across the system, although reductions were less than those in the increased sea-level scenario. The decline in groundwater levels caused by increased groundwater pumpage was less under the increased sea-level scenario than under the increased groundwater-pumpage scenario. The largest reductions in surface-water outflow from the salinity control structures occurred with increased sea level and increased groundwater pumpage. The model was designed specifically to evaluate the effect of groundwater pumpage on canal leakage at the surface-water-basin scale and thus may not be appropriate for (1) predictions that are dependent on data not included in the calibration process (for example, subdaily simulation of high-intensity events and travel times) and (or) (2) hydrologic conditions that are substantially different from those during the calibration and verification periods. The reliability of the model is limited by the conceptual model of the surface-water and groundwater system, the spatial distribution of physical properties, the scale and discretization of the system, and specified boundary conditions. Some of the model limitations are manifested in model errors. Despite these limitations, however, the model represents the complexities of the interconnected surface-water and groundwater systems that affect how the systems respond to groundwater pumpage, sea-level rise, and other hydrologic stresses. The model also quantifies the relative effects of groundwater pumpage and sea-level rise on the surface-water and groundwater systems.

Florida↗

USGS2014

One-Water Hydrologic Flow Model (MODFLOW-OWHM)

The One-Water Hydrologic Flow Model (MF-OWHM) is a MODFLOW-based integrated hydrologic flow model (IHM) that is the most complete version, to date, of the MODFLOW family of hydrologic simulators needed for the analysis of a broad range of conjunctive-use issues. Conjunctive use is the combined use of groundwater and surface water. MF-OWHM allows the simulation, analysis, and management of nearly all components of human and natural water movement and use in a physically-based supply-and-demand framework. MF-OWHM is based on the Farm Process for MODFLOW-2005 (MF-FMP2) combined with Local Grid Refinement (LGR) for embedded models to allow use of the Farm Process (FMP) and Streamflow Routing (SFR) within embedded grids. MF-OWHM also includes new features such as the Surface-water Routing Process (SWR), Seawater Intrusion (SWI), and Riparian Evapotrasnpiration (RIP-ET), and new solvers such as Newton-Raphson (NWT) and nonlinear preconditioned conjugate gradient (PCGN). This IHM also includes new connectivities to expand the linkages for deformation-, flow-, and head-dependent flows. Deformation-dependent flows are simulated through the optional linkage to simulated land subsidence with a vertically deforming mesh. Flow-dependent flows now include linkages between the new SWR with SFR and FMP, as well as connectivity with embedded models for SFR and FMP through LGR. Head-dependent flows now include a modified Hydrologic Flow Barrier Package (HFB) that allows optional transient HFB capabilities, and the flow between any two layers that are adjacent along a depositional or erosional boundary or displaced along a fault. MF-OWHM represents a complete operational hydrologic model that fully links the movement and use of groundwater, surface water, and imported water for consumption by irrigated agriculture, but also of water used in urban areas and by natural vegetation. Supply and demand components of water use are analyzed under demand-driven and supply-constrained conditions. From large- to small-scale settings, MF-OWHM has the unique set of capabilities to simulate and analyze historical, present, and future conjunctive-use conditions. MF-OWHM is especially useful for the analysis of agricultural water use where few data are available for pumpage, land use, or agricultural information. The features presented in this IHM include additional linkages with SFR, SWR, Drain-Return (DRT), Multi-Node Wells (MNW1 and MNW2), and Unsaturated-Zone Flow (UZF). Thus, MF-OWHM helps to reduce the loss of water during simulation of the hydrosphere and helps to account for “all of the water everywhere and all of the time.” In addition to groundwater, surface-water, and landscape budgets, MF-OWHM provides more options for observations of land subsidence, hydraulic properties, and evapotranspiration (ET) than previous models. Detailed landscape budgets combined with output of estimates of actual evapotranspiration facilitates linkage to remotely sensed observations as input or as additional observations for parameter estimation or water-use analysis. The features of FMP have been extended to allow for temporally variable water-accounting units (farms) that can be linked to land-use models and the specification of both surface-water and groundwater allotments to facilitate sustainability analysis and connectivity to the Groundwater Management Process (GWM). An example model described in this report demonstrates the application of MF-OWHM with the addition of land subsidence and a vertically deforming mesh, delayed recharge through an unsaturated zone, rejected infiltration in a riparian area, changes in demand caused by deficiency in supply, and changes in multi-aquifer pumpage caused by constraints imposed through the Farm Process and the MNW2 Package, and changes in surface water such as runoff, streamflow, and canal flows through SFR and SWR linkages.

Techniques and Methods↗

USGS2014

Utilizing dimensional analysis with observed data to determine the significance of hydrodynamic solutions in coastal hydrology

In this paper, the authors present an analysis of the magnitude of the temporal and spatial acceleration (inertial) terms in the surface-water flow equations and determine the conditions under which these inertial terms have sufficient magnitude to be required in the computations. Data from two South Florida field sites are examined and the relative magnitudes of temporal acceleration, spatial acceleration, and the gravity and friction terms are compared. Parameters are derived by using dimensionless numbers and applied to quantify the significance of the hydrodynamic effects. The time series of the ratio of the inertial and gravity terms from field sites are presented and compared with both a simplified indicator parameter and a more complex parameter called the Hydrodynamic Significance Number (HSN). Two test-case models were developed by using the SWIFT2D hydrodynamic simulator to examine flow behavior with and without the inertial terms and compute the HSN. The first model represented one of the previously-mentioned field sites during gate operations of a structure-managed coastal canal. The second model was a synthetic test case illustrating the drainage of water down a sloped surface from an initial stage while under constant flow. The analyses indicate that the times of substantial hydrodynamic effects are sporadic but significant. The simplified indicator parameter correlates much better with the hydrodynamic effect magnitude for a constant width channel such as Miami Canal than at the non-uniform North River. Higher HSN values indicate flow situations where the inertial terms are large and need to be taken into account.

Florida↗

USGS2013

Use of general purpose graphics processing units with MODFLOW

To evaluate the use of general-purpose graphics processing units (GPGPUs) to improve the performance of MODFLOW, an unstructured preconditioned conjugate gradient (UPCG) solver has been developed. The UPCG solver uses a compressed sparse row storage scheme and includes Jacobi, zero fill-in incomplete, and modified-incomplete lower-upper (LU) factorization, and generalized least-squares polynomial preconditioners. The UPCG solver also includes options for sequential and parallel solution on the central processing unit (CPU) using OpenMP. For simulations utilizing the GPGPU, all basic linear algebra operations are performed on the GPGPU; memory copies between the central processing unit CPU and GPCPU occur prior to the first iteration of the UPCG solver and after satisfying head and flow criteria or exceeding a maximum number of iterations. The efficiency of the UPCG solver for GPGPU and CPU solutions is benchmarked using simulations of a synthetic, heterogeneous unconfined aquifer with tens of thousands to millions of active grid cells. Testing indicates GPGPU speedups on the order of 2 to 8, relative to the standard MODFLOW preconditioned conjugate gradient (PCG) solver, can be achieved when (1) memory copies between the CPU and GPGPU are optimized, (2) the percentage of time performing memory copies between the CPU and GPGPU is small relative to the calculation time, (3) high-performance GPGPU cards are utilized, and (4) CPU-GPGPU combinations are used to execute sequential operations that are difficult to parallelize. Furthermore, UPCG solver testing indicates GPGPU speedups exceed parallel CPU speedups achieved using OpenMP on multicore CPUs for preconditioners that can be easily parallelized.

Ground Water↗

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

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↗

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↗

USGS2010

Effect of numerical dispersion as a source of structural noise in the calibration of a highly parameterized saltwater intrusion model

A model with a small amount of numerical dispersion was used to represent saltwater 7 intrusion in a homogeneous aquifer for a 10-year historical calibration period with one 8 groundwater withdrawal location followed by a 10-year prediction period with two groundwater 9 withdrawal locations. Time-varying groundwater concentrations at arbitrary locations in this low-10 dispersion model were then used as observations to calibrate a model with a greater amount of 11 numerical dispersion. The low-dispersion model was solved using a Total Variation Diminishing 12 numerical scheme; an implicit finite difference scheme with upstream weighting was used for 13 the calibration simulations. Calibration focused on estimating a three-dimensional hydraulic 14 conductivity field that was parameterized using a regular grid of pilot points in each layer and a 15 smoothness constraint. Other model parameters (dispersivity, porosity, recharge, etc.) were 16 fixed at the known values. The discrepancy between observed and simulated concentrations 17 (due solely to numerical dispersion) was reduced by adjusting hydraulic conductivity through the 18 calibration process. Within the transition zone, hydraulic conductivity tended to be lower than 19 the true value for the calibration runs tested. The calibration process introduced lower hydraulic 20 conductivity values to compensate for numerical dispersion and improve the match between 21 observed and simulated concentration breakthrough curves at monitoring locations. 22 Concentrations were underpredicted at both groundwater withdrawal locations during the 10-23 year prediction period.

Book↗

USGS2010

Evaluating the effect of Tikhonov regularization schemes on predictions in a variable-density groundwater model

Calibration of highly‐parameterized numerical models typically requires explicit Tikhonovtype regularization to stabilize the inversion process. This regularization can take the form of a preferred parameter values scheme or preferred relations between parameters, such as the preferred equality scheme. The resulting parameter distributions calibrate the model to a user‐defined acceptable level of model‐to‐measurement misfit, and also minimize regularization penalties on the total objective function. To evaluate the potential impact of these two regularization schemes on model predictive ability, a dataset generated from a synthetic model was used to calibrate a highly-parameterized variable‐density SEAWAT model. The key prediction is the length of time a synthetic pumping well will produce potable water. A bi‐objective Pareto analysis was used to explicitly characterize the relation between two competing objective function components: measurement error and regularization error. Results of the Pareto analysis indicate that both types of regularization schemes affect the predictive ability of the calibrated model.

Conference Paper↗