A three-dimensional, steady-state groundwater-flow model representing 2003-13 mean hydrologic conditions was developed and calibrated to assess groundwater and lake-water exchanges and the effects of groundwater withdrawals and precipitation on water levels in lakes in the northeast Twin Cities Metropolitan Area, Minnesota. The USGS groundwater-flow model program MODFLOW-NWT version 1.0.8 was used to simulate groundwater flow in the approximately 1,000-square-mile area of the northeast Twin Cities Metropolitan Area and western Wisconsin. Water levels were below normal for several lakes in the northeast Twin Cities Metropolitan Area during 2003 through 2013. Previous periods of low lake-water levels generally correlate with periods of below-normal precipitation but increases in groundwater withdrawals and land-use changes have put into question whether recent lake-water-level declines are due solely to declines in precipitation. The groundwater-flow model was developed to provide a thorough understanding of regional groundwater and surface-water exchanges under different groundwater withdrawal and precipitation scenarios. This USGS data release contains all of the input and output files for the simulations described in the associated U.S. Geological Survey Scientific Investigations Report 2016-5139B (https://doi.org/10.3133/sir20165139B). This data release also contains input and output data for ancillary soil-water balance (SWB) models used to simulate runoff to a subset of lakes and areal groundwater recharge.

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

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

In 2017 the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, published a calibrated numerical groundwater-flow model and associated model documentation report that evaluated the effects of potential groundwater withdrawals on groundwater flow and availability in the North Fork Red River aquifer in southwest Oklahoma. The results of groundwater-availability scenarios run on the calibrated numerical groundwater-flow model could be used by the Oklahoma Water Resources Board to reevaluate the maximum annual yield of groundwater from the North Fork Red River aquifer in Oklahoma. The numerical groundwater-flow model was built on a hydrogeologic framework and a conceptual groundwater-flow model derived from previously published and newly collected hydrologic data. A hydrogeologic framework is a three-dimensional representation of the aquifer and the surrounding geologic units at a scale that captures the regional controls on groundwater flow. The hydrogeologic framework for the North Fork Red River aquifer included a definition of the aquifer extent and potentiometric surface, as well as a description of the textural and hydraulic properties of aquifer materials. A conceptual groundwater-flow model is a simplified description of the major inflow and outflow sources (hydrologic boundaries) of a groundwater-flow system as well as an accounting of the estimated mean flows from those sources (water budget) for a specified period of time. The hydrogeologic framework and conceptual model are necessary constraints used in the construction and calibration of a scientifically defensible numerical groundwater-flow model that reasonably represents the groundwater-flow system. A finite-difference numerical groundwater-flow model of the North Fork Red River aquifer was constructed by using MODFLOW-2005 with the Newton formulation solver (MODFLOW-NWT). Data inputs for each package were specified in machine-readable text files. The numerical model of the North Fork Red River aquifer had 385 rows, 460 columns, about 27,600 active cells of 886 by 886 ft (270 by 270 meters), and 2 convertible layers. The top layer (layer 1) represented the undifferentiated Quaternary alluvium and terrace deposits with variable thickness determined from the hydrogeologic framework, and the bottom layer (layer 2) represented the Permian bedrock with a nominal thickness of about 100 feet. The model active area was created from the North Fork Red River aquifer extent and expanded in some areas to ensure that each active cell was in connection with at least one other active cell. One terrace lobe in northern Beckham County was not included in the model active area because it was almost separated spatially and hydraulically from the rest of the North Fork Red River aquifer. The numerical model was temporally discretized into 408 monthly transient stress periods (each with 2 time steps) representing the period 1980–2013. An initial steady-state stress period, in which the groundwater-flow equation had no storage component, represented mean annual inflows to and outflows from the aquifer and produced a solution that was used as the initial condition for subsequent transient stress periods. The numerical model was constructed in units of meters and days. 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/sir20175098).

In 2017 the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, published a calibrated numerical groundwater-flow model and associated model documentation report that evaluated the effects of potential groundwater withdrawals on groundwater flow and availability in the North Fork Red River aquifer in southwest Oklahoma. The results of groundwater-availability scenarios run on the calibrated numerical groundwater-flow model could be used by the Oklahoma Water Resources Board to reevaluate the maximum annual yield of groundwater from the North Fork Red River aquifer in Oklahoma. The numerical groundwater-flow model was built on a hydrogeologic framework and a conceptual groundwater-flow model derived from previously published and newly collected hydrologic data. A hydrogeologic framework is a three-dimensional representation of the aquifer and the surrounding geologic units at a scale that captures the regional controls on groundwater flow. The hydrogeologic framework for the North Fork Red River aquifer included a definition of the aquifer extent and potentiometric surface, as well as a description of the textural and hydraulic properties of aquifer materials. A conceptual groundwater-flow model is a simplified description of the major inflow and outflow sources (hydrologic boundaries) of a groundwater-flow system as well as an accounting of the estimated mean flows from those sources (water budget) for a specified period of time. The hydrogeologic framework and conceptual model are necessary constraints used in the construction and calibration of a scientifically defensible numerical groundwater-flow model that reasonably represents the groundwater-flow system. A finite-difference numerical groundwater-flow model of the North Fork Red River aquifer was constructed by using MODFLOW-2005 with the Newton formulation solver (MODFLOW-NWT). Data inputs for each package were specified in machine-readable text files. The numerical model of the North Fork Red River aquifer had 385 rows, 460 columns, about 27,600 active cells of 886 by 886 ft (270 by 270 meters), and 2 convertible layers. The top layer (layer 1) represented the undifferentiated Quaternary alluvium and terrace deposits with variable thickness determined from the hydrogeologic framework, and the bottom layer (layer 2) represented the Permian bedrock with a nominal thickness of about 100 feet. The model active area was created from the North Fork Red River aquifer extent and expanded in some areas to ensure that each active cell was in connection with at least one other active cell. One terrace lobe in northern Beckham County was not included in the model active area because it was almost separated spatially and hydraulically from the rest of the North Fork Red River aquifer. The numerical model was temporally discretized into 408 monthly transient stress periods (each with 2 time steps) representing the period 1980–2013. An initial steady-state stress period, in which the groundwater-flow equation had no storage component, represented mean annual inflows to and outflows from the aquifer and produced a solution that was used as the initial condition for subsequent transient stress periods. The numerical model was constructed in units of meters and days. 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/sir20175098).

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

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.

The U.S. Geological Survey constructed a steady-state numerical groundwater flow model in cooperation with Des Moines Water Works (DMWW) to simulate groundwater flow conditions in the Des Moines River alluvial aquifer (DMRA) during winter low-flow conditions typical of December 2018-2020. The Des Moines River alluvial aquifer (DMRA) is an important source of water for Des Moines Water Works (DMWW), the municipal water utility that serves residential and commercial water needs in the city of Des Moines, Iowa and surrounding municipalities. A comprehensive understanding of groundwater flow processes in the DMRA is needed for DMWW to make decisions related to the management of this water resource. A three-layered model was constructed using MODFLOW-NWT to simulate an area of about 15 square kilometers near Prospect Park in Des Moines, Iowa. The model has 130 rows and 130 columns of cells within the model boundary. Parameter ESTimation software (PEST) was used for model calibration to assess and optimize performance of individual parameters including the horizontal and vertical hydraulic conductivity of the various units, evapotranspiration rate, and recharge rate. 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/ofr20211110).

The U.S. Geological Survey constructed a steady-state numerical groundwater flow model in cooperation with Des Moines Water Works (DMWW) to simulate groundwater flow conditions in the Des Moines River alluvial aquifer (DMRA) during winter low-flow conditions typical of December 2018-2020. The Des Moines River alluvial aquifer (DMRA) is an important source of water for Des Moines Water Works (DMWW), the municipal water utility that serves residential and commercial water needs in the city of Des Moines, Iowa and surrounding municipalities. A comprehensive understanding of groundwater flow processes in the DMRA is needed for DMWW to make decisions related to the management of this water resource. A three-layered model was constructed using MODFLOW-NWT to simulate an area of about 15 square kilometers near Prospect Park in Des Moines, Iowa. The model has 130 rows and 130 columns of cells within the model boundary. Parameter ESTimation software (PEST) was used for model calibration to assess and optimize performance of individual parameters including the horizontal and vertical hydraulic conductivity of the various units, evapotranspiration rate, and recharge rate. 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/ofr20211110).