This data release contains a three-dimensional groundwater flow model with example applications using MODFLOW-2000. The calibrated model is able to simulate observed decadal-scale climate-driven fluctuations in the groundwater system as well as observed shorter-term pumping-related fluctuations. Example model simulations show that the timing and location of the effects of groundwater pumping vary markedly depending on the pumping location. The complete description for the models in Gannett et al., 2012.

An existing three-dimensional model (MODFLOW-2000) by Petkewich and Campbell (2007) (https://pubs.usgs.gov/sir/2007/5126/) was updated to simulate six predictive water-management scenarios that were created to simulate potential changes in groundwater flow and groundwater-level conditions in the Mount Pleasant, South Carolina area. The model was recalibrated to conditions from 1900 to 2015. Simulations included six scenarios: (1) maximize Mount Pleasant Waterworks reverse-osmosis plant capacity by increasing groundwater withdrawals from 3.9 million gallons per day (Mgal/d) in 2015 to 8.6 Mgal/d from the Middendorf aquifer; (2) same as Scenario 1, but with the addition of a 0.5 Mgal/d supply well in the Middendorf aquifer near Moncks Corner, SC; (3) same as Scenario 1, but with the addition of a 1.5 Mgal/d supply well in the Middendorf aquifer near Moncks Corner, SC; (4) maximize Mount Pleasant Waterworks well capacity by increasing withdrawals from the Middendorf aquifer from 3.9 Mgal/d in 2015 to 10.2 Mgal/d (5) minimizing Mount Pleasant Waterworks surface-water purchase from the Charleston Water System by adding supply wells and increasing withdrawals from the Middendorf aquifer from 3.9 Mgal/d in 2015 to 12.2 Mgal/d; and (6) same as Scenario 1, but with the addition of quarterly model stress periods to simulate seasonal variations in the groundwater withdrawals. This USGS data release contains all of the input and output files for the simulations described in the associated model documentation report (https://doi.org/10.3133/sir20175128).

An existing three-dimensional model (MODFLOW-2000) by Petkewich and Campbell (2007) (https://pubs.usgs.gov/sir/2007/5126/) was updated to simulate six predictive water-management scenarios that were created to simulate potential changes in groundwater flow and groundwater-level conditions in the Mount Pleasant, South Carolina area. The model was recalibrated to conditions from 1900 to 2015. Simulations included six scenarios: (1) maximize Mount Pleasant Waterworks reverse-osmosis plant capacity by increasing groundwater withdrawals from 3.9 million gallons per day (Mgal/d) in 2015 to 8.6 Mgal/d from the Middendorf aquifer; (2) same as Scenario 1, but with the addition of a 0.5 Mgal/d supply well in the Middendorf aquifer near Moncks Corner, SC; (3) same as Scenario 1, but with the addition of a 1.5 Mgal/d supply well in the Middendorf aquifer near Moncks Corner, SC; (4) maximize Mount Pleasant Waterworks well capacity by increasing withdrawals from the Middendorf aquifer from 3.9 Mgal/d in 2015 to 10.2 Mgal/d (5) minimizing Mount Pleasant Waterworks surface-water purchase from the Charleston Water System by adding supply wells and increasing withdrawals from the Middendorf aquifer from 3.9 Mgal/d in 2015 to 12.2 Mgal/d; and (6) same as Scenario 1, but with the addition of quarterly model stress periods to simulate seasonal variations in the groundwater withdrawals. This USGS data release contains all of the input and output files for the simulations described in the associated model documentation report (https://doi.org/10.3133/sir20175128).

A three-dimensional groundwater flow model, MODFLOW-OWHM, was developed to provide a better understanding of the hydrogeology of the Lucerne Valley Groundwater Basin, California. The model was used to investigate the historical groundwater storage loss and subsidence associated with anthropogenic groundwater demands. The model was calibrated to 1942 through 2016 conditions. This USGS data release contains all of the input and output files for the simulation described in the associated model documentation report https://doi.org/10.3133/sir20225048

A three-dimensional groundwater flow model, MODFLOW-OWHM, was developed to provide a better understanding of the hydrogeology of the Lucerne Valley Groundwater Basin, California. The model was used to investigate the historical groundwater storage loss and subsidence associated with anthropogenic groundwater demands. The model was calibrated to 1942 through 2016 conditions. This USGS data release contains all of the input and output files for the simulation described in the associated model documentation report https://doi.org/10.3133/sir20225048

This model archive data release includes all models used to characterize the magnitude, spatial distribution and timing of groundwater (GW) flow through lakebed sediments to Upper Klamath Lake (UKL), Oregon, described in the associated journal article (https://doi.org/10.1016/j.scitotenv.2020.142768). One-dimensional vertical models of GW flow (MODFLOW-2005) and solute transport (MT3D-USGS) were calibrated (UCODE) to 2014 observed dissolved silica (Si, 0.2-micron filtered) porewater concentrations in the upper 0.1 m of lakebed sediment to estimate GW flow and Si exchange across the lakebed interface. The Si-based calibrated GW flow rates were then used in conjunction with observed dissolved phosphate-phosphorus (PP) porewater concentrations in the upper 0.1 m of lakebed sediment to estimate the amount of PP reacted during upward flow through the lakebed sediment and the PP discharge to the lake. One-dimensional, vertical GW flow and heat transport models (VS2DH) were calibrated (UCODE) to 2015 and 2017 observed lakebed temperatures to provide estimates of GW-inflow rates at multiple UKL locations. Calibrated GW inflows were greatest in the spring and decreased through the summer. The magnitude and timing of the GW-lake water exchange estimates obtained from these methods were compared to rates obtained from a generalized cross-sectional GW flow model (MODFLOW-NWT) with time-varying recharge. The cross-sectional GW flow model demonstrated that snow-melt GW recharge could be transported rapidly to the lake due to the relatively high permeability and low specific storage of the surrounding volcanic rocks explaining the greater GW discharge to the lake in the spring. This USGS data release contains all the input and output files for the simulations described in the associated journal article (https://doi.org/10.1016/j.scitotenv.2020.142768).

This model archive data release includes all models used to characterize the magnitude, spatial distribution and timing of groundwater (GW) flow through lakebed sediments to Upper Klamath Lake (UKL), Oregon, described in the associated journal article (https://doi.org/10.1016/j.scitotenv.2020.142768). One-dimensional vertical models of GW flow (MODFLOW-2005) and solute transport (MT3D-USGS) were calibrated (UCODE) to 2014 observed dissolved silica (Si, 0.2-micron filtered) porewater concentrations in the upper 0.1 m of lakebed sediment to estimate GW flow and Si exchange across the lakebed interface. The Si-based calibrated GW flow rates were then used in conjunction with observed dissolved phosphate-phosphorus (PP) porewater concentrations in the upper 0.1 m of lakebed sediment to estimate the amount of PP reacted during upward flow through the lakebed sediment and the PP discharge to the lake. One-dimensional, vertical GW flow and heat transport models (VS2DH) were calibrated (UCODE) to 2015 and 2017 observed lakebed temperatures to provide estimates of GW-inflow rates at multiple UKL locations. Calibrated GW inflows were greatest in the spring and decreased through the summer. The magnitude and timing of the GW-lake water exchange estimates obtained from these methods were compared to rates obtained from a generalized cross-sectional GW flow model (MODFLOW-NWT) with time-varying recharge. The cross-sectional GW flow model demonstrated that snow-melt GW recharge could be transported rapidly to the lake due to the relatively high permeability and low specific storage of the surrounding volcanic rocks explaining the greater GW discharge to the lake in the spring. This USGS data release contains all the input and output files for the simulations described in the associated journal article (https://doi.org/10.1016/j.scitotenv.2020.142768).

This model archive provides the necessary documentation of the numerical models developed for the Central Sands Lake study in central Wisconsin and will be included as a technical appendix (Appendix C) in the report to the Wisconsin State Legislature by the Wisconsin Department of Natural Resources (WDNR) in response to 2017 Wisconsin Act 10. This legislation directed DNR to determine whether existing and potential groundwater withdrawals are causing or are likely to cause significant reduction of mean seasonal water levels at Pleasant Lake, Long Lake, and Plainfield Lake (s. 281.34(7m)(2)(b), Wis. Stats.) in Waushara County, Wisconsin. To evaluate the potential hydrologic connection between groundwater withdrawals and the nearby study lakes, hydrologic models were created that focused on the lakes of interest and yet were large enough to cover a broad enough region to extend to the major hydrologic boundaries of the natural flow system. The areas near the lakes require finer-scale grid discretization (or spacing) to better represent the lakes and streams in the model, but also need to cover a large enough area to include the groundwater withdrawal locations that have the potential to cause reduction in water levels in the lakes. To accomplish these goals, three groundwater models were created: a regional model extending to major hydrologic boundaries; and two inset models, inheriting boundaries from the regional model but focused near the lakes. Each of the inset models, in turn, included a detailed area close to the lakes surrounded by an area at the same spatial scale as the regional model. To support WDNR in evaluating the connection between groundwater withdrawals and lake levels, a representative time period was required over which to compare land use with and without irrigated agriculture and for WDNR to evaluate potential lake stage and flux changes related to irrigated agriculture. WDNR chose the climate period of 1981-2018 to be representative of a typical period and provided two land use scenarios—one with no irrigated agriculture and one with assumed crop rotations similar to current conditions—to simulate with groundwater models to, then, compare lake responses with. As a result, simulations over this climate record are not intended to recreate the history of 1981-2018 because land use changed over that time. These runs are, instead, intended to provide a basis on which to compare land use with and without irrigation-related groundwater withdrawals based on the current arrangement of land use and a varied climatic record. Groundwater withdrawals focused on irrigated-agriculture-related water use because greater than 95% of groundwater withdrawal in the two inset models around the study lakes is for irrigated agriculture water use. The period of 2012-2018 was used for parameter estimation (synonymously referred to as “history matching”) for the groundwater models. This time period was chosen because it includes the most complete water use records to simulate groundwater withdrawals. History matching was performed using groundwater elevations, lake stages, and streamflow observations over the 2012-2018 time period and processed observations derived from those raw data. Climatic data were incorporated into the model using a soil-water balance approach. A soil water balance model (Westenbroek and others, 2021) was constructed at the scale of the regional groundwater model to both calculate recharge based on land use and climate, and in the long-term climate-period runs, to estimate water use required by irrigated agriculture to apply as well boundary conditions in the groundwater model in the absence of reported water use values over that period. The model archive presents all the inputs needed to run the models, the model software, information on history matching to estimate parameters of the model, model scenario files, and model outputs that the user should be able to recreate using the model files in this archive.

This model archive provides the necessary documentation of the numerical models developed for the Central Sands Lake study in central Wisconsin and will be included as a technical appendix (Appendix C) in the report to the Wisconsin State Legislature by the Wisconsin Department of Natural Resources (WDNR) in response to 2017 Wisconsin Act 10. This legislation directed DNR to determine whether existing and potential groundwater withdrawals are causing or are likely to cause significant reduction of mean seasonal water levels at Pleasant Lake, Long Lake, and Plainfield Lake (s. 281.34(7m)(2)(b), Wis. Stats.) in Waushara County, Wisconsin. To evaluate the potential hydrologic connection between groundwater withdrawals and the nearby study lakes, hydrologic models were created that focused on the lakes of interest and yet were large enough to cover a broad enough region to extend to the major hydrologic boundaries of the natural flow system. The areas near the lakes require finer-scale grid discretization (or spacing) to better represent the lakes and streams in the model, but also need to cover a large enough area to include the groundwater withdrawal locations that have the potential to cause reduction in water levels in the lakes. To accomplish these goals, three groundwater models were created: a regional model extending to major hydrologic boundaries; and two inset models, inheriting boundaries from the regional model but focused near the lakes. Each of the inset models, in turn, included a detailed area close to the lakes surrounded by an area at the same spatial scale as the regional model. To support WDNR in evaluating the connection between groundwater withdrawals and lake levels, a representative time period was required over which to compare land use with and without irrigated agriculture and for WDNR to evaluate potential lake stage and flux changes related to irrigated agriculture. WDNR chose the climate period of 1981-2018 to be representative of a typical period and provided two land use scenarios—one with no irrigated agriculture and one with assumed crop rotations similar to current conditions—to simulate with groundwater models to, then, compare lake responses with. As a result, simulations over this climate record are not intended to recreate the history of 1981-2018 because land use changed over that time. These runs are, instead, intended to provide a basis on which to compare land use with and without irrigation-related groundwater withdrawals based on the current arrangement of land use and a varied climatic record. Groundwater withdrawals focused on irrigated-agriculture-related water use because greater than 95% of groundwater withdrawal in the two inset models around the study lakes is for irrigated agriculture water use. The period of 2012-2018 was used for parameter estimation (synonymously referred to as “history matching”) for the groundwater models. This time period was chosen because it includes the most complete water use records to simulate groundwater withdrawals. History matching was performed using groundwater elevations, lake stages, and streamflow observations over the 2012-2018 time period and processed observations derived from those raw data. Climatic data were incorporated into the model using a soil-water balance approach. A soil water balance model (Westenbroek and others, 2021) was constructed at the scale of the regional groundwater model to both calculate recharge based on land use and climate, and in the long-term climate-period runs, to estimate water use required by irrigated agriculture to apply as well boundary conditions in the groundwater model in the absence of reported water use values over that period. The model archive presents all the inputs needed to run the models, the model software, information on history matching to estimate parameters of the model, model scenario files, and model outputs that the user should be able to recreate using the model files in this archive.

A three-dimensional groundwater flow model (MODFLOW200-FMP1_1) of the Central Valley in California was developed to aid water managers in understanding how water moves through the aquifer system, to predict water-supply scenarios, and to address issues related to water competition. The USGS Groundwater Resources Program made a detailed assessment of groundwater availability of the Central Valley aquifer system, which includes: (1) the present status of groundwater resources; (2) how these resources have changed over time; and (3) tools to assess system responses to stresses from future human uses and climate variability and change. This effort builds on previous investigations, such as the USGS Central Valley Regional Aquifer System and Analysis (CV-RASA) project and several other groundwater studies in the Valley completed by Federal, State and local agencies at differing scales. The principal product of this new assessment is a tool referred to as the Central Valley Hydrologic Model (CVHM) that accounts for integrated, variable water supply and demand, and simulates surface-water and groundwater-flow across the entire Central Valley system. The current model was extended to incorporate a slightly larger geographic area, has a finer spatial and temporal discretization, uses a more-detailed depiction of subsurface geology. In addition, the model utilizes a modified version of MODFLOW2000 (version 1.15.03) to include an updated and refined Farm Process (FMP1) to simulate groundwater and surface-water flow, irrigated agriculture, land subsidence, and other key processes in the Central Valley on a monthly basis for April 1961 through September 2003. This USGS data release contains all of the input and output files for the simulation and calibration of the CVHM described in the associated model documentation report (https://pubs.er.usgs.gov/publication/pp1766).