A three-dimensional, numerical groundwater flow model of the Hualapai Plateau and Truxton Basin was developed to assist water-resource managers in understanding the potential effects of projected groundwater withdrawals on groundwater levels and storage in the basin. The Truxton Basin Hydrologic Model (TBHM) was developed using previously published data, as well as, new geophysical data collected as part of this investigation (Ball, 2020; Kennedy, 2020). TBHM is a transient model that simulates the hydrologic system between the years 1976–2140, including three hypothetical groundwater withdrawal scenarios between 2020–2140. The predevelopment (before 1976) groundwater system is simulated using a steady-state model. The TBHM study area encompasses 1,200 square miles (mi2), including the Truxton basin (approximately 75 mi2) and much of the surrounding Hualapai Plateau. Potential future pumping scenarios are simulated in Truxton Basin at the existing Truxton well field and at a hypothetical new pumping location in the northern portion of the basin. 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/sir20205017E).

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

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

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

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