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).

A three-dimensional groundwater flow model was developed to characterize groundwater resources the uppermost principal aquifers in the Williston structural basin in parts of Montana, North Dakota, and South Dakota in the United States and of Manitoba and Saskatchewan in Canada as part of a detailed assessment of the groundwater availability of the area. The uppermost principal aquifers are comprised of the glacial, lower Tertiary, and Upper Cretaceous aquifer systems. The model was developed as a part of the U.S. Geological Survey Water Availability and Use Science Program's effort to conduct large-scale multidisciplinary regional studies of groundwater availability. The numerical model is intended to be used to (1) simulate hydrologic scenarios of interest to groundwater managers and to advance the understanding of groundwater budgets and components including recharge, discharge, and aquifer storage for the entire system, (2) compute historical and projected system response to natural and anthropogenic stresses, and (3) evaluate potential hydrologic monitoring programs at a scale relevant to basin-wide water-management decisions. The three-dimensional groundwater-flow model was developed using the numerical modeling software, MODFLOW-NWT. The steady-state (mean) hydrological conditions included data from 1981 to 2005, and transient (temporally-varying) conditions included a combination of a steady state period with data prior to 1960, and a transient period from 1961 to 2005. The model was calibrated by attempting to match simulated and measured or estimated hydraulic heads, differences in hydraulic heads between aquifers, stream base flow, and measured flow at flowing artesian wells. Sub-regional water budgets for the model area were produced with ZONEBUDGET. This USGS data release contains all of the input and output files for the model described in the associated model documentation report (https://doi.org/10.3133/sir201755158). This data release also includes (1) MODFLOW-NWT (version 1.0.9) source code, and (2) ZONEBUDGET source code.

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 three-dimensional, groundwater flow model was developed with the numerical code MODFLOW-NWT to represent changes in groundwater pumping and aquifer recharge in the Northern Atlantic Coastal Plain aquifer system from Long Island, New York to North Carolina. The model was constructed using existing hydrogeologic and geospatial information to represent the aquifer system geometry, boundaries, and hydraulic properties of the 19 separate regional aquifers and confining units within the aquifer system. The model was calibrated using an inverse modeling parameter-estimation (PEST) technique to conditions from 1986 to 2008, the period for which data are most complete and reliable. The simulation period for this analysis spanned from predevelopment to future conditions, from 1900 to 2058. The model was used to advance the understanding of groundwater budgets and components including recharge, discharge, and aquifer storage for the entire system and for each of the statewide systems; compute historical and recent system response and project future system response to development at a scale relevant to basinwide water-management decisions; and evaluate options for hydrologic monitoring of system changes. The report ‘Documentation of a groundwater flow model developed to assess groundwater availability in the Northern Atlantic Coastal Plain aquifer system from Long Island, New York, to North Carolina: U.S. Geological Survey Scientific Investigations Report 2016–5076' (https://doi.org/10.3133/sir20165076) documents the model design and calibration, as well as several simulations to test model construction assumptions. The report 'Assessment of groundwater availability in the Northern Atlantic Coastal Plain aquifer system from Long Island, New York, to North Carolina: U.S. Geological Survey Professional Paper 1829' (https://doi.org/10.3133/pp1829) documents water-availability simulations and the resulting analysis and discussion. This USGS data release contains all of the input and output files for the simulations described in the associated reports (https://doi.org/10.3133/sir20165076) and (https://doi.org/10.3133/pp1829). This data release also includes (1) MODFLOW-NWT source code, (2) the PEST files and source code used for model calibration, and (3) the ZONEBUDGET input files and source code used for the groundwater availability analysis.

A three-dimensional groundwater flow model was developed to characterize groundwater resources and the interaction of groundwater with streams and other hydrologic features in the Northern High Plains aquifer. The Northern High Plains aquifer is generally unconfined; most of the water withdrawn from the aquifer is used for irrigation. A modified version of a previously published soil-water-balance (SWB) model estimates recharge using climatic, soils, land cover data, in addition to data for groundwater withdrawals for irrigation. The SWB output was adjusted in areas where surface water is used for irrigation. The groundwater flow model results were calibrated using parameter estimation to measured groundwater levels and estimated stream base flows. The model was designed as a tool for regional evaluations of groundwater resources and of groundwater interactions with streams and other hydrologic features resulting from current or forecasted conditions. This USGS data release contains all of the input and output files for the model described in the associated model documentation report (https://doi.org/10.3133/sir20165153). This data release also includes (1) MODFLOW-NWT (version 1.0.9) source code, and (2) SWB source code in two formats.

A previously developed model (https://doi.org/10.3133/sir20175098) was coupled with downscaled climate model data to determine the impact of climate variability on base flow and groundwater storage in the North Fork Red River aquifer, Oklahoma. The North Fork Red River aquifer is an alluvial aquifer that discharges groundwater to the North Fork Red River, which provides inflow to Lake Altus, an important water source for the surrounding communities. The impact of climate variability on hydrologic systems and the resulting effects on basins has become an important topic in assessing future water resources. Global climate projections from general circulation models, including the Coupled Model Intercomparison Project Phase 5 (CMIP5), have been developed to improve the understanding of climate science and forecast future climatic conditions. Due to the impact of climate variations on groundwater resources, it is important to communicate the ranges of results with water resource managers. To approximate a range in future base flow conditions and flow into Lake Altus, the Coupled Model Intercomparison Project Phase 5 climate data was downscaled to watershed scale using monthly Bias-Correction Spatial Disaggregation techniques. A time-series of scaling factors were developed and interpolated for three climate scenarios (central tendency, warmer/drier, and less warm-wetter) representing a range of future climate conditions for the period 2045–2074. These scaling factors were then applied to an existing soil-water-balance model dataset with climate data for the baseline period 1980–2009 to produce recharge and evapotranspiration estimations for this future period. The downscaled climate data was applied to the finite-difference numerical groundwater-flow model of the North Fork Red River aquifer using MODFLOW-2005 with the Newton formulation solver (MODFLOW-NWT) which was temporally discretized into 373 monthly transient stress periods representing the period 1980–2010. Three climate scenarios (central tendency, warmer/drier, and less warm/wetter) representing a range of future climate conditions for the period 2045–2074 were simulated. This USGS data release contains all of the input and output files for the simulations described in the associated journal article (http://doi.org/10.1007/s10040-020-02230-x).

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 Salt Fork Arkansas River and Chikaskia River alluvial aquifers by using MODFLOW-2005 with the Newton formulation solver (MODFLOW-NWT). The model included the Chikaskia River alluvial aquifer, which is classified as a minor aquifer by the OWRB and is hydrologically connected to the Salt Fork Arkansas River alluvial aquifer. 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 feet. 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. At the time of this publication (2025), a hydrologic investigation and determination of the MAY for the Salt Fork Arkansas River alluvial aquifer had not been completed. The U.S. Geological Survey, in cooperation with the OWRB, conducted a hydrologic investigation and evaluated the simulated effects of potential groundwater withdrawals on groundwater flow and availability in the Salt Fork Arkansas River alluvial aquifer in northern Oklahoma for a study period spanning 1980–2020. Fifteen simulations are included in this data release: a simulation for the calibrated numerical groundwater-flow model, 9 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 the input and output files for the simulations described in the associated model documentation report (https://doi.org/10.3133/sir20255043).

A previously calibrated MODFLOW-NWT groundwater-flow model (https://doi.org/10.3133/sir20205118) was used to determine the effects of climate variability under a range of future climatic conditions on groundwater resources in the reach 1 of the Washita River alluvial aquifer in western Oklahoma. The study area focuses on reach 1 of the Washita River alluvial aquifer; the entire Washita River alluvial aquifer consists of four administrative sections, or reaches, that are designated as reaches 1–4 by the Oklahoma Water Resources Board (OWRB, 2012). To approximate a range in future base-flow conditions in reach 1 of the Washita River alluvial aquifer and base-flow into Foss Reservoir, the Coupled Model Intercomparison Project Phase 5 Global Climate Model climate data were downscaled to watershed scale using monthly Bias-Correction Spatial Disaggregation techniques. A time series of scaling factors was developed and spatially interpolated for three climate scenarios (central tendency, warmer/drier, and less warm-wetter) representing a range of future climate conditions for the period 2050–2079. These scaling factors were then applied to an existing soil-water-balance model (https://doi.org/10.3133/sir20205118) with climate data for the baseline period 1985–2014 to produce recharge and evapotranspiration estimations for this future period. The downscaled climate data were applied to the groundwater-flow model of the reach 1 of the Washita River alluvial aquifer using MODFLOW-NWT. This data release contains the input and output files for the scenarios described in the associated model documentation report (https://doi.org/10.3133/sir20245082).

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 Boone and Roubidoux aquifers in northeastern Oklahoma by using MODFLOW-NWT (version 1.1.4) with the Newton formulation solver to simulate groundwater flow and account for the drying and rewetting of cells within the groundwater-flow model. The numerical groundwater-flow model was discretized into four layers consisting of 354 rows by 261 columns with a 2,000-feet by 2,000-feet cell size. The model layers were used to simulate the Western Interior Plains confining system, the Boone aquifer, the Ozark confining unit, and the Roubidoux aquifer. The model was temporally discretized into one steady-state stress period followed by 456 monthly transient stress periods spanning from January 1980 to December 2017. The steady-state stress period typically consisted of mean annual inputs from January 1980 to December 2017, but inputs from 1979 were included for some of the simulations.

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 Boone and Roubidoux aquifers in northeastern Oklahoma by using MODFLOW-NWT (version 1.1.4) with the Newton formulation solver to simulate groundwater flow and account for the drying and rewetting of cells within the groundwater-flow model. The numerical groundwater-flow model was discretized into four layers consisting of 354 rows by 261 columns with a 2,000-feet by 2,000-feet cell size. The model layers were used to simulate the Western Interior Plains confining system, the Boone aquifer, the Ozark confining unit, and the Roubidoux aquifer. The model was temporally discretized into one steady-state stress period followed by 456 monthly transient stress periods spanning from January 1980 to December 2017. The steady-state stress period typically consisted of mean annual inputs from January 1980 to December 2017, but inputs from 1979 were included for some of the simulations.