A three-dimensional, groundwater flow model (MODFLOW-NWT) was developed to examine groundwater storage changes in the Quincy Basin, Washington. The model was calibrated to conditions from 1920 to 2013. The model was used to (1) determine the change in groundwater storage from 1920 to 2013 , and (2) simulate the potential effects of increases in pumping, decrease in irrigation recharge, and increases in streamflow in Crab Creek by 100 cubic feet per second and 500 cubic feet per second. 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/sir20185162).

This groundwater-flow model archive contains all of the input and output files for an inset MODFLOW-NWT model extracted from the northern (Wisconsin) half of a published USGS steady-state regional model of the Upper Fox River Basin in the U.S. Upper Midwest. The construction and details of the published USGS steady-state model of the Upper Fox River Basin is outlined in the U.S. Geological Survey Scientific Investigations Report 2018-5038 (https://doi.org/10.3133/sir20185038). The regional model is archived in the data release at https://doi.org/10.5066/F76D5R5V. The extracted model was used to demonstrate an innovative new method for delinating fen distribution and discharge using the MODFLOW UZF package. The extracted model incorporates the Mukwonago River Basin, a 10-digit hydrologic unit code (HUC10) basin occupying 86.2 mi2 (223 km2) in southeastern Wisconsin. The extracted model was used to demonstrate how regional and local flow patterns can be enhanced by adding a version of the UZF file that automatically inserts “seepage drains” in cells where the water table is near the land surface (within the “undulation depth”). Details on the extracted model construction and calibration, including preparation of the “stripped-down” UZF file central to the proposed fen delineation method can be found in the supporting information of the journal article in Groundwater (https://doi.org/10.1111/gwat.12931). This USGS data release contains all of the input and output files for the simulations described in the journal article in Groundwater (https://doi.org/10.1111/gwat.12931).

This groundwater-flow model archive contains all of the input and output files for an inset MODFLOW-NWT model extracted from the northern (Wisconsin) half of a published USGS steady-state regional model of the Upper Fox River Basin in the U.S. Upper Midwest. The construction and details of the published USGS steady-state model of the Upper Fox River Basin is outlined in the U.S. Geological Survey Scientific Investigations Report 2018-5038 (https://doi.org/10.3133/sir20185038). The regional model is archived in the data release at https://doi.org/10.5066/F76D5R5V. The extracted model was used to demonstrate an innovative new method for delinating fen distribution and discharge using the MODFLOW UZF package. The extracted model incorporates the Mukwonago River Basin, a 10-digit hydrologic unit code (HUC10) basin occupying 86.2 mi2 (223 km2) in southeastern Wisconsin. The extracted model was used to demonstrate how regional and local flow patterns can be enhanced by adding a version of the UZF file that automatically inserts “seepage drains” in cells where the water table is near the land surface (within the “undulation depth”). Details on the extracted model construction and calibration, including preparation of the “stripped-down” UZF file central to the proposed fen delineation method can be found in the supporting information of the journal article in Groundwater (https://doi.org/10.1111/gwat.12931). This USGS data release contains all of the input and output files for the simulations described in the journal article in Groundwater (https://doi.org/10.1111/gwat.12931).

The U.S. Geological Survey (USGS), in cooperation with the Oklahoma Water Resources Board (OWRB), constructed a finite-difference numerical groundwater-flow model of the Washita River aquifer by using MODFLOW-2005 (Harbaugh, 2005) with the Newton formulation solver (MODFLOW-NWT). The 1973 Oklahoma Groundwater Law requires that the OWRB conduct hydrologic investigations of the State’s aquifers to determine the maximum annual yield (MAY) for each groundwater basin. The MAY is defined as the total amount of fresh groundwater that can be annually withdrawn while allowing a minimum 20-year life of that groundwater basin. For alluvium and terrace groundwater basins, the life requirement is satisfied if, after 20 years of MAY withdrawals, 50 percent of the groundwater basin (hereinafter referred to as an “aquifer”) retains a saturated thickness of at least 5 ft. Once a MAY has been established, the amount of land owned or leased by a groundwater-use permit applicant determines the annual volume of water allocated to that groundwater-use permit applicant. The annual volume of groundwater allocated per acre of land is known as the equal-proportionate-share (EPS) pumping rate. The OWRB issued a final order on November 13, 1990, that established the MAY (81,840 and 46,935 acre-feet per year [acre-ft/yr]) and EPS pumping rate (1.5 and 1.0 acre-foot per acre per year) for reaches 3 and 4, respectively, of the Washita River aquifer in southern Oklahoma. Because more than 20 years have elapsed since the final order was issued, the USGS, in cooperation with the OWRB, conducted an updated hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater flow and availability in the Washita River aquifer in southern Oklahoma. Reach 3 extends from near Anadarko, Okla., to Alex, Okla., and reach 4 extends from near Alex to south of Davis, Okla. Twenty-four simulations are included in this data release: a simulation for the calibrated numerical groundwater-flow model, 18 scenario simulations to evaluate the EPS pumping rate, 4 scenario simulations to evaluate groundwater storage over a 50-year period, and 1 scenario simulation to evaluate effects of a hypothetical drought. 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/sir20235072).

The U.S. Geological Survey (USGS), in cooperation with the Oklahoma Water Resources Board (OWRB), constructed a finite-difference numerical groundwater-flow model of the Washita River aquifer by using MODFLOW-2005 (Harbaugh, 2005) with the Newton formulation solver (MODFLOW-NWT). The 1973 Oklahoma Groundwater Law requires that the OWRB conduct hydrologic investigations of the State’s aquifers to determine the maximum annual yield (MAY) for each groundwater basin. The MAY is defined as the total amount of fresh groundwater that can be annually withdrawn while allowing a minimum 20-year life of that groundwater basin. For alluvium and terrace groundwater basins, the life requirement is satisfied if, after 20 years of MAY withdrawals, 50 percent of the groundwater basin (hereinafter referred to as an “aquifer”) retains a saturated thickness of at least 5 ft. Once a MAY has been established, the amount of land owned or leased by a groundwater-use permit applicant determines the annual volume of water allocated to that groundwater-use permit applicant. The annual volume of groundwater allocated per acre of land is known as the equal-proportionate-share (EPS) pumping rate. The OWRB issued a final order on November 13, 1990, that established the MAY (81,840 and 46,935 acre-feet per year [acre-ft/yr]) and EPS pumping rate (1.5 and 1.0 acre-foot per acre per year) for reaches 3 and 4, respectively, of the Washita River aquifer in southern Oklahoma. Because more than 20 years have elapsed since the final order was issued, the USGS, in cooperation with the OWRB, conducted an updated hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater flow and availability in the Washita River aquifer in southern Oklahoma. Reach 3 extends from near Anadarko, Okla., to Alex, Okla., and reach 4 extends from near Alex to south of Davis, Okla. Twenty-four simulations are included in this data release: a simulation for the calibrated numerical groundwater-flow model, 18 scenario simulations to evaluate the EPS pumping rate, 4 scenario simulations to evaluate groundwater storage over a 50-year period, and 1 scenario simulation to evaluate effects of a hypothetical drought. 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/sir20235072).

A numerical groundwater flow model of the Hualapai Valley Basin, using MODFLOW-NWT, was developed to assist water-resource managers in understanding the potential effects of projected groundwater withdrawals on groundwater levels in the basin. Hualapai Valley Basin is a broad, internally drained, intermountain desert basin in Mohave County, northwestern Arizona. Basin-fill aquifers are the primary groundwater source for many desert communities, and the residents, commerce, and agriculture in and near to the Hualapai Valley Basin must rely on such groundwater to meet water needs. As in many parts of the western United States, population growth in this part of Arizona is substantial. From 2000 to 2018 the population of the City of Kingman, Arizona, grew from 20,069 to 30,314, an increase of 51 percent, whereas the population of Mohave County grew from 155,062 to 209,550, an increase of 35 percent. Water managers in Mohave County have raised concern about the potential for future groundwater development and additional stresses on the groundwater system in the Hualapai Valley Basin. In particular, the City of Kingman, Ariz., water supply is primarily groundwater withdrawn from the Kingman subbasin of the Hualapai Valley Basin, northeast of the city. The potential effects of future water development on the City of Kingman well field have become a top concern to regional water-resource managers. To address these concerns the Hualapai Valley Hydrologic Model (HVHM) simulates the hydrologic system for the years 1935 through 2219, including future withdrawal scenarios that simulate large-scale agricultural expansion with and without enhanced groundwater recharge from potential new infiltration basin projects. HVHM is a highly parameterized model (75,586 adjustable parameters) capable of simulating grid-scale variability in aquifer properties (for example, conductivity, specific yield, and specific storage) and system stresses (for instance, natural recharge and groundwater withdrawals). System stresses were partially adopted from a previously-published groundwater model (Tillman and others, 2013). Parameter estimation and uncertainty quantification were performed using an iterative ensemble smoother software (PESTPP-IES) to produce an ensemble of models fit to historical data. Two future scenarios were simulated with a subset of the posterior parameter ensemble comprising the 40 best-fit realizations. In scenario 1, future pumping was simulated to increase linearly from 2019 through 2029 and then held constant through 2219. Scenario 2 includes the same specified future pumping, but also simulates enhanced recharge at proposed infiltration basins throughout the Kingman subbasin beginning in 2019. 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/sir20215077).

A numerical groundwater flow model of the Hualapai Valley Basin, using MODFLOW-NWT, was developed to assist water-resource managers in understanding the potential effects of projected groundwater withdrawals on groundwater levels in the basin. Hualapai Valley Basin is a broad, internally drained, intermountain desert basin in Mohave County, northwestern Arizona. Basin-fill aquifers are the primary groundwater source for many desert communities, and the residents, commerce, and agriculture in and near to the Hualapai Valley Basin must rely on such groundwater to meet water needs. As in many parts of the western United States, population growth in this part of Arizona is substantial. From 2000 to 2018 the population of the City of Kingman, Arizona, grew from 20,069 to 30,314, an increase of 51 percent, whereas the population of Mohave County grew from 155,062 to 209,550, an increase of 35 percent. Water managers in Mohave County have raised concern about the potential for future groundwater development and additional stresses on the groundwater system in the Hualapai Valley Basin. In particular, the City of Kingman, Ariz., water supply is primarily groundwater withdrawn from the Kingman subbasin of the Hualapai Valley Basin, northeast of the city. The potential effects of future water development on the City of Kingman well field have become a top concern to regional water-resource managers. To address these concerns the Hualapai Valley Hydrologic Model (HVHM) simulates the hydrologic system for the years 1935 through 2219, including future withdrawal scenarios that simulate large-scale agricultural expansion with and without enhanced groundwater recharge from potential new infiltration basin projects. HVHM is a highly parameterized model (75,586 adjustable parameters) capable of simulating grid-scale variability in aquifer properties (for example, conductivity, specific yield, and specific storage) and system stresses (for instance, natural recharge and groundwater withdrawals). System stresses were partially adopted from a previously-published groundwater model (Tillman and others, 2013). Parameter estimation and uncertainty quantification were performed using an iterative ensemble smoother software (PESTPP-IES) to produce an ensemble of models fit to historical data. Two future scenarios were simulated with a subset of the posterior parameter ensemble comprising the 40 best-fit realizations. In scenario 1, future pumping was simulated to increase linearly from 2019 through 2029 and then held constant through 2219. Scenario 2 includes the same specified future pumping, but also simulates enhanced recharge at proposed infiltration basins throughout the Kingman subbasin beginning in 2019. 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/sir20215077).

A three-dimensional groundwater flow model, constructed in MODFLOW-NWT, was developed to evaluate the groundwater flow system of the Kitsap Peninsula, west-central Washington. A transient model was constructed to simulate groundwater flow for January 1985–December 2012 using annual stress periods for 1985–2004 and monthly stress periods for 2005–2012. The model was used to simulate six hydrologic scenarios, including simulations of a steady-state system, no-pumping and return flows, 15-percent increase in current withdrawals in all wells, 80-percent decrease in outdoor water to simulate effects of conservation efforts, 15-percent decrease in recharge from precipitation to simulate a drought, and particle tracking to determine flow paths. This data release contains the input and output files for the simulations described in the associated model documentation report (https://doi.org/10.3133/sir20165052).

A three-dimensional groundwater flow model, constructed in MODFLOW-NWT, was developed to evaluate the groundwater flow system of the Kitsap Peninsula, west-central Washington. A transient model was constructed to simulate groundwater flow for January 1985–December 2012 using annual stress periods for 1985–2004 and monthly stress periods for 2005–2012. The model was used to simulate six hydrologic scenarios, including simulations of a steady-state system, no-pumping and return flows, 15-percent increase in current withdrawals in all wells, 80-percent decrease in outdoor water to simulate effects of conservation efforts, 15-percent decrease in recharge from precipitation to simulate a drought, and particle tracking to determine flow paths. This data release contains the input and output files for the simulations described in the associated model documentation report (https://doi.org/10.3133/sir20165052).

A three-dimensional groundwater flow model (MODFLOW-NWT) of the Columbia Plateau Regional aquifer (CPRAS) in Washington, Oregon, and Idaho was developed to provide an integrated understanding of the hydrologic system to implement effective water-resource management strategies. The U.S. Geological Survey (USGS) Groundwater Resources Program assessed the groundwater availability as part of a national study of regional systems (https://pubs.usgs.gov/circ/1323/). The CPRAS assessment includes the status of groundwater resources, how these resources have changed over time, and development and application of tools to estimate system responses to stresses from future uses and climate variability and change. A major product of this assessment is a numerical groundwater-flow model of the system. Two models were constructed to simulate groundwater flow in the CPRAS: a steady-state predevelopment model representing conditions before large-scale pumping and irrigation altered the system, and a transient model representing the period 1900–2007. Construction of the model, development and calibration of the steady-state and transient simulations, as well as, a commingling scenario, is documented in the Scientific Investigations Report 2014-5127 (https://doi.org/10.3133/sir20145127). Two additional scenarios were completed and documented in the U.S. Geological Survey Professional Paper 1817 (https://doi.org/10.3133/pp1817). One scenario represents long-term equilibrium under 2007 conditions and the other modified the equilibrium conditions to account for potential increased pumpage under projected temperature increases with climate change. The model is a useful tool for investigating water supply, water demand, management strategies, groundwater-surface water exchanges, and potential effects of changing climate on the hydrologic system. This USGS data release contains all of the input and output files for the simulations described in the associated model documentation reports: U.S. Geological Survey Professional Paper 1817 (https://doi.org/10.3133/pp1817) and the modeling report that documents the construction and calibration of the model, Scientific Investigations Report 2014-5127 (https://doi.org/10.3133/sir20145127).