Geology topics

Jonathan P. Traylor

Publications and source records attributed to Jonathan P. Traylor.

11 recordsLinked to original sources

USGS2026

Scenarios to assess the future water availability in the Mississippi River Valley Alluvial Aquifer for the Cache River and Grand Prairie Regions of Arkansas

The U.S. Geological Survey, as part of the Arkansas Groundwater Initiative, developed forecast scenarios using previously calibrated MODFLOW 6 groundwater models that focused on the Cache and Grand Prairie Critical Groundwater Areas to assess the impact of future climate and water management strategies on the Mississippi River Valley alluvial aquifer. A Soil Water Balance model was used to forecast recharge and irrigation water use. The forecast scenario period was from January 1, 2019, through December 31, 2055, with monthly stress periods. Twenty scenarios were simulated and included seven alternate climate forecasts, five 13 general groundwater pumping reduction scenarios (round 1), and groundwater pumping reduction scenarios by crop type and for the Bayou Meto Water Management Project and Grand Prairie Area Demonstration Project (round 2). Declines in saturated thickness within the Cache Critical Groundwater Area were larger for 18 of the 20 scenarios as compared to outside of the Critical Groundwater Area. The largest average increase in saturated thickness inside the Critical Groundwater Area was 6.4 m which occurred for the round 1, 50 percent reduction scenario. Automatic reductions in groundwater pumping by MODFLOW 6 in the Cache simulation ranged from 0.02 to 13.1 percent of total groundwater pumping. For the Grand Prairie model domain, the average change in saturated thickness of the Mississippi River Valley alluvial aquifer inside the Critical Groundwater Area for the forecast period ranged between -6.6 to 1.7 m. The average saturated thickness of the Mississippi River Valley alluvial aquifer inside the Grand Prairie Critical Groundwater Area declined for 16 of the 20 scenarios. The average reduction in requested groundwater pumping for all scenarios inside the Grand Prairie Critical Groundwater Area was 25.1 percent, and the largest reduction was 46.5 percent.

ESS Open Archive↗

USGS2025

An inset groundwater-flow model to evaluate the effects of layering configuration on model calibration and assess managed aquifer recharge near Shellmound, Mississippi

The U.S. Geological Survey has developed a high-resolution inset groundwater-flow model in the Mississippi Delta as part of an interdisciplinary collaboration coordinated by the Mississippi Alluvial Plain project to provide a tool that stakeholders can use to support water-resource management decisions. Groundwater withdrawals from the Mississippi River Valley alluvial (MRVA) aquifer have been vital to support agricultural production in the region, but substantial groundwater-level declines near Shellmound, Mississippi, have caused concerns for long-term sustainability of the aquifer. To better understand the subsurface and try to mitigate the long-term groundwater-level declines, stakeholders have undertaken actions including a Groundwater Transfer and Injection Pilot (GTIP) project using a riverbank filtration-based managed aquifer recharge approach. The pilot project consisted of extracting groundwater near the Tallahatchie River and reinjecting it into the aquifer 3 kilometers west where water levels have substantially declined. A high-resolution airborne electromagnetic (AEM) survey was also completed to collect electrical resistivity data to support the GTIP project and the development of the groundwater model. The inset groundwater-flow model was developed to (1) integrate the AEM data into the optimal layering configuration of the MRVA aquifer that the available observation data can support through calibration, and (2) assess the potential effect of the GTIP project on the groundwater levels. The AEM data were processed into three different layering configurations leading to the development of model A (18 layers), model B (16 layers), and model C (8 layers), all at a 100- x 100-meter cell spatial resolution using the U.S. Geological Survey modular finite-difference flow model 6 code with Newton-Raphson formulation. The model development process integrated recent advances in modeling, such as the incorporation of AEM data, the use of outputs from the soil-water-balance (SWB) model, and the Aquaculture and Irrigation Water-Use Model, and was facilitated by robust automation using the open-source python packages Modflow-setup and SFRmaker. Using Parameter Estimation ++ Iterative Ensemble Smoother, the three numerical groundwater-flow models (models A, B, and C) were calibrated against a set of observations, which included aquifer groundwater levels, streamflows, stream stage, and aquifer transmissivity. Results indicate that the detailed representation of MRVA aquifer layers in model A produced the best calibrated model by history matching, and the integration of data representing surficial connectivity played a key role in improving groundwater recharge and enhancing the ability of the model to match groundwater levels in the cone of depression. A forecast model simulated the managed aquifer recharge approach, and the results indicated that, given average irrigation and recharge conditions (2010–15), the GTIP project has the potential to induce groundwater-level increases of as much as 3 meters around the injection site, but a sustained increase would require repetition in subsequent years of water transfer at 2022 rates or above.

Mississippi↗

USGS2024

Inset groundwater-flow models for the Cache and Grand Prairie Critical Groundwater Areas, northeastern Arkansas

The water resources in the Mississippi alluvial plain, located in parts of Missouri, Kentucky, Tennessee, Mississippi, Louisiana, and Arkansas, supports a multibillion-dollar agricultural industry that relies heavily on pumping of groundwater for irrigation of crops and aquaculture. The primary source of groundwater for agricultural-related pumping is the Mississippi River Valley alluvial aquifer, which has declined in storage for decades; secondary groundwater sources include the middle Claiborne aquifer and Wilcox aquifer system. Two areas in northeastern Arkansas that lie within the Mississippi alluvial plain, part of the Cache and Grand Prairie regions, have been designated as Critical Groundwater Areas owing to decades of groundwater declines that resulted from past and current water use. The multidisciplinary Mississippi Alluvial Plain project, led by the U.S. Geological Survey, and funded by their Water Availability and Use Science Program, included objectives to develop numerical groundwater models in focus regions, including the part of the Cache and Grand Prairie regions of northeastern Arkansas. Two inset models were developed using the child model capabilities of MODFLOW 6, the U.S. Geological Survey’s Modular Hydrologic Model simulation software. Both models, called the Cache model and Grand Prairie model, simulated the groundwater system and surface-water/groundwater interactions for the Mississippi River Valley alluvial aquifer and underlying Tertiary-age aquifers and confining units to the Midway confining unit. Each model was spatially discretized into 500-meter x 500-meter orthogonal cells on a grid with 5-meter constant-thickness vertical layers that represented the Mississippi River Valley alluvial aquifer and increasing thickness layers for the aquifers and confining units below the alluvial aquifer. The Cache and Grand Prairie models were calibrated with the PEST++ iterative ensemble smoother Version 5 and employed high dimensional parameterization schemes of 13,740 and 30,436 parameters, respectively. The Cache mean absolute residual for groundwater-level observations within each model domain for the priority well was 1.58 meters. Grand Prairie mean absolute residuals for the alluvial aquifer and middle Claiborne aquifer groundwater-level observations were 2.71 and 10.78 meters, respectively. The groundwater budgets for the Cache and Grand Prairie models were characterized by substantial outflows to irrigation wells, which constituted about 52 and 54 percent of all outflows, with the primary source of water to those wells being releases from unconfined aquifer storage.

Arkansas↗

USGS2024

Simulation of groundwater-flow dynamics in the U.S. Northern High Plains driven by multi-model estimates of surficial aquifer recharge

There is growing interest in incorporating higher-resolution groundwater modeling within the framework of large-scale land surface models (LSMs), including processes such as three-dimensional flow, variable soil saturation, and surface water/groundwater interactions. Conversely, complex groundwater models (e.g., the U.S. Geological Survey Groundwater-Flow Model, MODFLOW) often use simpler representations of land surface dynamics (e.g., surface vegetation, evapotranspiration, recharge) and may benefit from higher process fidelity and temporal resolutions in these inputs. This study investigates the potential of improving groundwater representation in LSMs and land surface dynamics in MODFLOW through forcing MODFLOW with recharge from LSMs. Groundwater simulations build on an existing and well-calibrated MODFLOW model of the U.S. Northern High Plains aquifer, a hydrologically complex basin under the dual impacts of conversion of native vegetation to intense irrigated agricultural fields and climate change. Simulated groundwater recharge from four different land models are used to drive MODFLOW groundwater simulations. Results show relatively large discrepancies between recharge estimates among simulations. Forcing MODFLOW using recharge simulated by some of the LSMs in place of a simple water balance model marginally improves MODFLOW groundwater simulation. Further, our results support the efficacy of coupling LSMs to a sophisticated groundwater model such as MODFLOW. The coupling results in notable improvements in matching the historical groundwater levels through reduction of the skewness coefficient in percent bias histogram (from 1.50 and 1.41 in original LSMs to 0.44 and 0.27, respectively, when MODFLOW is forced by groundwater recharge from LSMs) and reduction of bias. This modeling effort seeks to identify the best compromise between comprehensive land surface processes from global LSMs and advanced representation of groundwater from regional models.

Journal of Hydrology↗

USGS2023

Effects of auto-adaptive localization on a model calibration using ensemble methods

Simulations of the natural systems for environmental decision-making typically benefit from a highly parameterized approach (Hunt et al. 2007; Doherty and Hunt 2010), which enhances the flow of information contained in state observations to the parameters and improves application to decision support. However, parameter estimation (PE) with highly parameterized environmental models using traditional approaches (e.g., Doherty and Hunt 2010) is computationally intensive. Attempts at addressing the computational burden include improved computing approaches (e.g., Schreüder 2009; Hunt et al. 2010) and advances in algorithmic approaches (e.g., Tonkin and Doherty 2005; Welter et al. 2012, 2015). Recently, the iterative ensemble smoother (IES) approach (Chen and Oliver 2013; White 2018; White et al. 2020a) has greatly improved the efficiency of the PE calibration process compared to previous algorithms while concurrently providing nonlinear estimates of uncertainty (Hunt et al. 2021).

Groundwater↗

USGS2023

An integrated hydrologic model to support the Central Platte Natural Resources District Groundwater Management Plan, central Nebraska

The groundwater and surface-water supply of the Central Platte Natural Resources District supports a large agricultural economy from the High Plains aquifer and Platte River, respectively. This study provided the Central Platte Natural Resources District with an advanced numerical modeling tool to assist with the update of their Groundwater Management Plan. An integrated hydrologic model, called the Central Platte Integrated Hydrologic Model, was constructed using the MODFLOW-One-Water Hydrologic Model code with the Newton solver. This code integrates climate, landscape, surface water, and groundwater-flow processes in a fully coupled approach. Model framework included 163 rows; 327 columns; 2,640 feet cell sides; and 3 vertical layers. A predevelopment model simulated steady-state hydrologic conditions prior to April 30, 1895, and a development period model discretized into 610 stress periods simulated transient hydrologic conditions from May 1, 1895, to December 31, 2016, using 170 biannual stress periods from 1895 to 1980, and monthly stress periods from May 1, 1980, to December 31, 2016. Calibration of the Central Platte Integrated Hydrologic Model involved two phases: a manual adjustment of parameters, followed by the automated calibration completed using BeoPEST that was facilitated by the employment of the singular value decomposition-assist features of PEST that specified 50 super parameters assembled from the 435 adjustable parameters and Tikhonov regularization. The average absolute groundwater-level residuals for model layers one, two, and three were 6.1, 12.4, and 7.4 feet, respectively. Calibrated horizontal hydraulic conductivity was about 70, 32, and 35 feet per day for layers 1, 2, and 3, respectively. The largest development period inflow to groundwater was recharge from deep percolation past the root zone, averaging 1,122,257 acre-feet per year (2.7 inches per year), and the largest outflow was to irrigation wells, averaging 693,171 acre-feet per year (10.2 inches per year for the Central Platte Natural Resources District). Other substantial groundwater outflows included evapotranspiration and base flow. For the total development period, there was a net change in storage of −122,393 acre-feet per year (−0.3 inch per year). The calibrated Central Platte Integrated Hydrologic Model was used to simulate eight different potential future climate and irrigation pumping conditions from January 1, 2017, to December 31, 2049. Simulated future groundwater levels within the Central Platte Natural Resources District varied significantly between scenarios and locally, from 13.8 feet below to 7.6 feet above baseline 1982 groundwater levels. Most areas exhibited groundwater-level declines for the drought scenarios and rises for the alternate irrigation scenarios. Changes in scenario groundwater levels correlated with the relations between farm net recharge and irrigation pumping. Linear “first order second moment” techniques indicated that the uncertainty in projected groundwater altitudes was reduced by 15.33 feet through model calibration.

Nebraska↗

USGS2023

System-scale airborne electromagnetic surveys in the lower Mississippi River Valley support multidisciplinary applications

The lower Mississippi River Valley spans over 200,000 square kilometres in parts of seven states, encompassing areas of critical groundwater supplies, natural hazards, infrastructure, and low-lying coastal regions. From 2018 - 2022, the U.S. Geological Survey acquired over 82,000 line-kilometres of airborne electromagnetic, radiometric, and magnetic data over this region to provide comprehensive and systematic information about subsurface geologic and hydrologic properties that support multiple scientific and societal interests. Most of the data were acquired on a regional grid of west-east flight lines separated by 3 - 6 kilometres; however, several high-resolution inset grids with line spacing as close as 200 m were acquired in targeted areas of interest. Approximately 8,000 line-kilometres were acquired along streams and rivers to characterise the potential for surface water-groundwater connection, and another 6,000 line-kilometres were acquired along the Mississippi and Arkansas River levees to characterise this critical infrastructure. Here, we present a summary of the data along with several examples of how they are being used to inform regional groundwater model development, inferences of groundwater salinity, identification of faults in the New Madrid seismic zone, and levee infrastructure.

lower Mississippi River Valley↗

USGS2021

Assessment of water availability in the Osage Nation using an integrated hydrologic-flow model

The Osage Nation of northeastern Oklahoma, conterminous with Osage County, covers about 2,900 square miles. The area is primarily rural with 62 percent of the land being native prairie grass, and much of the area is used for cattle ranching and extraction of petroleum and natural gas. Protection of water rights are important to the Osage Nation because of its reliance on cattle ranching and the potential for impairment of water quality by petroleum extraction. Additionally, the potential for future population increases, demands for water from neighboring areas such as the Tulsa metropolitan area, and expansion of petroleum and natural-gas extraction on water resources of this area further the need for the Osage Nation to better understand its water availability. Therefore, the U.S. Geological Survey, in cooperation with the Osage Nation, completed a hydrologic investigation to assess the status and availability of surface-water and groundwater resources in the Osage Nation. A transient integrated hydrologic-flow model was constructed using the U.S. Geological Survey fully integrated hydrologic-flow model called the MODFLOW One-Water Hydrologic Model. The integrated hydrologic-flow model, called the Osage Nation Integrated Hydrologic Model (ONIHM), was constructed and uses an orthogonal grid of 276 rows and 289 columns, and each grid cell measures 1,312.34 feet (ft; 400 meters) per side, with eight variably thick vertical layers that represented the alluvial and bedrock aquifers within the study area, including the alluvial aquifer, the Vamoosa-Ada aquifer, and the minor Pennsylvanian bedrock aquifers, and the confining units. Landscape and groundwater-flow processes were simulated for two periods: (1) the 1950–2014 period from January 1950 through September 2014 and (2) the forecast period from October 2014 through December 2099. The 1950–2014 period ONIHM simulated past conditions using measured or estimated inputs, and the forecast-period ONIHM simulated three separate potential forecast conditions under constant dry, average, or wet climate conditions using calibrated input values from the 1950–2014 period ONIHM. The 1950–2014 period ONIHM was calibrated by linking the Parameter Estimation software (PEST) with the MODFLOW One-Water Hydrologic Model. PEST uses statistical parameter estimation techniques to identify the best set of parameter values to minimize the difference between measured or estimated calibration targets and their simulated equivalent values (residuals). Tikhonov regularization and singular-value decomposition-assist features of PEST were used during the calibration process. The 1950–2014 period ONIHM was calibrated to 713 measured groundwater levels at 195 wells; 95,636 estimated monthly mean groundwater levels at 124 wells; 5,307 measured streamflows at 13 streamgages; and 8,679 simulated mean monthly streamflows at 10 streamgages extracted from a surface-water model by adjusting 231 parameters. The estimated groundwater-level observations and streamflows were included as observations to improve the spatial and temporal density of observation targets during calibration. The best set of parameter values obtained during the calibration process of the 1950–2014 model was then used as the input parameter values for the forecast model simulations. A comparison of the calibration targets to their corresponding simulated values indicated that the model adequately reproduced streamflows and groundwater levels for some streamgages and wells and underestimated streamflows and groundwater levels at other locations. Measured and simulated streamflows correlated adequately with a coefficient of determination of 0.938, as did water levels with a coefficient of determination of 0.795. The 1950–2014 period ONIHM underestimated certain groundwater levels and streamflows, but generally measured or estimated calibration targets correlated well with simulated equivalents, which indicated that the model can adequately simulate the response of the hydrologic system to stresses in the 1950–2014 and forecast periods. In the 1950–2014 period ONIHM, the calibrated mean horizontal hydraulic conductivity for layer 1 alluvial aquifer was 30.7 feet per day, and the seven lower layers had a calibrated mean horizontal hydraulic conductivity of less than 3.3 feet per day. The mean calibrated groundwater-level residual was 16.6 ft, and the mean calibrated streamflow residual of the Arkansas River at Ralston, Oklahoma, streamgage (U.S. Geological Survey station 07152500) was within 6 percent (373 cubic feet per second) of mean measured streamflow for the 1950–2014 period ONIHM. The ONIHM simulated landscape fluxes of precipitation; groundwater applied by irrigation wells; evapotranspiration from precipitation, groundwater, and irrigation; runoff from precipitation; and deep percolation from precipitation. The largest loss of water from the landscape was evapotranspiration from precipitation with a calibrated mean annual outflow of 32 inches (in.): mean annual precipitation was about 36 in. Calibrated mean annual runoff and deep percolation (recharge to the water table) rates were 4.7 inches per year (in/yr) and 0.70 in/yr, respectively, for the 1950–2014 period ONIHM. The calibrated 1950–2014 period ONIHM groundwater fluxes included net farm net recharge (calculated as the difference between the inflow of recharge to the water table and the outflow of evapotranspiration from the water table such that negative values indicate that evapotranspiration from the water table was greater than deep percolation [recharge to the water table] and vice versa). Net farm net recharge was the largest flux from the groundwater system with a mean annual net outflow of 153.4 cubic feet per second. Stream leakage was the largest flux to the groundwater system with a mean annual net inflow of 152.5 cubic feet per second, indicating that, on average, the groundwater/surface-water interaction was a “losing” system where stream water leaked into the subsurface and recharged the water table. Simulated monthly trends demonstrated that net stream leakage was the largest inflow to the groundwater-flow system for 10 of the 12 months; for the other 2 months (January and March), farm net recharge (January) and net storage (March) were the largest inflow to the groundwater-flow system. A saline groundwater interface map was created for the study and compared to the water levels from the final stress period of the 1950–2014 model to identify the presence of fresh/marginal groundwater throughout the study area. Fresh/marginal groundwater was characterized as groundwater with less than 1,500 milligrams per liter of total dissolved solids. Fresh/marginal groundwater thickness ranged from 0 to 438.2 ft within the study area. The thickest regions of fresh/marginal groundwater were in the eastern part of the study area near Sand Creek, Bird Creek, and Hominy Creek and in the Arkansas River alluvial aquifer in the region downstream from the Arkansas River at Ralston, Okla. Like the 1950–2014 model, forecast model results for the landscape indicated that transpiration from precipitation was the largest flux out of the landscape for all three forecasts, constituting 77, 73, and 58 percent of precipitation for the dry, average, and wet forecasts, respectively. The dry and average forecast landscape fluxes demonstrated similar trends and magnitudes, whereas the wet forecast landscape fluxes indicated the largest changes compared to the average forecast fluxes. Most notably, runoff increased from a mean of 1.1 and 1.6 in/yr for the dry and average forecasts, respectively, to 10 in/yr for the wet forecast. Similar changes occurred for the other wet forecast landscape fluxes. The calibrated 1950–2014 period ONIHM simulated three forecasts to assess the effects of potential climatic changes on the hydrologic system from October 2014 to December 2099. The three forecasts simulated theoretical dry, average, and wet conditions using precipitation and potential evapotranspiration datasets from selected years in the calibrated 1950–2014 period ONIHM. Annual precipitation amounts were 26.89, 35.47, and 50.73 in. for the dry, average, and wet forecasts, respectively. Groundwater-flow component forecast results indicated that stream leakage is always a net inflow to the groundwater-flow system for dry, average, and wet conditions, meaning the study area stream network is always predominantly a “losing” regime where stream water infiltrates into the underlying aquifer. Storage was only a net outflow from the groundwater-flow system and indicated a replenishment to groundwater storage that resulted in an increase in groundwater levels only during the wet forecast. Further, these gains in groundwater storage for the wet forecast occurred only during February through June. Mean fresh/marginal groundwater saturated thicknesses were 125 and 126 ft for the dry and average forecast conditions, respectively, and wet forecast average thickness was 145 ft and ranged from 0 to 443 ft. The spatial extents of fresh/marginal groundwater at the end of the dry, average, and wet forecast model periods (December 2099) did not change substantially from the end of the 1950–2014 model period (September 2014).

Kansas, Oklahoma↗

USGS2020

Groundwater availability of the Northern High Plains aquifer in Colorado, Kansas, Nebraska, South Dakota, and Wyoming

Executive Summary The Northern High Plains aquifer underlies about 93,000 square miles of Colorado, Kansas, Nebraska, South Dakota, and Wyoming and is the largest subregion of the nationally important High Plains aquifer. Irrigation, primarily using groundwater, has supported agricultural production since before 1940, resulting in nearly $50 billion in sales in 2012. In 2010, the High Plains aquifer had the largest groundwater withdrawals of any major aquifer system in the United States. Nearly one-half of those withdrawals were from the Northern High Plains aquifer, which has little hydrologic interaction with parts of the aquifer farther south. Land-surface elevation ranges from more than 7,400 feet (ft) near the western edge to less than 1,100 ft near the eastern edge. Major stream primarily flow west to east and include the Big Blue River, Elkhorn River, Loup River, Niobrara River, Republican River and Platte River with its two forks—the North Platte River and South Platte River. Population in the Northern High Plain aquifer area is sparse with only 2 cities having a population greater than 30,000. Droughts across much of the area from 2001 to 2007, combined with recent (2004–18) legislation, have heightened concerns regarding future groundwater availability and highlighted the need for science-based water-resource management. Groundwater models with the capability to provide forecasts of groundwater availability and related stream base flows from the Northern High Plains aquifer were published recently (2016) and were used to analyze groundwater availability. Stream base flows are generally the dominant component of total streamflow in the Northern High Plains aquifer, and total streamflows or shortages thereof define conjunctive management triggers, at least in Nebraska. Groundwater availability was evaluated through comparison of aquifer-scale water budgets compared for periods before and after major groundwater development and across selected future forecasts. Groundwater-level declines and the forecast amount of groundwater in storage in the aquifer also were examined. Major Findings Aquifer losses to irrigation withdrawals increased greatly from 1940 to 2009 and were the largest average 2000–9 outflow (49 percent of total). Basin to basin groundwater flows were not a large part of basin water budgets. Development of irrigated land and associated withdrawals were not uniform across the Northern High Plains aquifer, and different parts of the Northern High Plains aquifer responded differently to agricultural development. For the Northern High Plains aquifer, areas with high recharge and low evapotranspiration had the most streamflow, and most streams only remove water from the aquifer. Results of a baseline future forecast indicated that groundwater levels declined overall, indicating an overdraft of the aquifer when climate was about average and agricultural development was held at the same state as 2009. Results of two human stresses future forecasts indicated that increases of 13 percent or 23 percent in agricultural development, mostly near areas of previous development, caused increases in groundwater pumping of 8 percent or 11 percent, and resulted in continued groundwater-level declines, at rates 0.3 or 0.5 million acre-feet per year larger than the baseline forecast. Results of environmental stresses forecasts (generated from two downscalings of global climate model outputs) compared with the baseline forecast indicated that even though annual precipitation was nearly the same, differences in temperature and a redistribution of precipitation from the spring to the growing season (from about May 1 through September 30), created a large (12–15 percent) decrease in recharge to the aquifer. For the two environmental stresses forecasts, temperature and precipitation were distributed about the same among basins of the Northern High Plains aquifer, but the amounts were different.

Colorado, Kansas, Nebraska, South Dakota, Wyoming↗

USGS2016

Groundwater-flow model of the northern High Plains aquifer in Colorado, Kansas, Nebraska, South Dakota, and Wyoming

The High Plains aquifer is a nationally important water resource underlying about 175,000 square miles in parts of eight states: Colorado, Kansas, Oklahoma, Nebraska, New Mexico, South Dakota, Texas, and Wyoming. Droughts across much of the Northern High Plains from 2001 to 2007 have combined with recent (2004) legislative mandates to elevate concerns regarding future availability of groundwater and the need for additional information to support science-based water-resource management. To address these needs, the U.S. Geological Survey began the High Plains Groundwater Availability Study to provide a tool for water-resource managers and other stakeholders to assess the status and availability of groundwater resources. A transient groundwater-flow model was constructed using the U.S. Geological Survey modular three-dimensional finite-difference groundwater-flow model with Newton-Rhapson solver (MODFLOW–NWT). The model uses an orthogonal grid of 565 rows and 795 columns, and each grid cell measures 3,281 feet per side, with one variably thick vertical layer, simulated as unconfined. Groundwater flow was simulated for two distinct periods: (1) the period before substantial groundwater withdrawals, or before about 1940, and (2) the period of increasing groundwater withdrawals from May 1940 through April 2009. A soil-water-balance model was used to estimate recharge from precipitation and groundwater withdrawals for irrigation. The soil-water-balance model uses spatially distributed soil and landscape properties with daily weather data and estimated historical land-cover maps to calculate spatial and temporal variations in potential recharge. Mean annual recharge estimated for 1940–49, early in the history of groundwater development, and 2000–2009, late in the history of groundwater development, was 3.3 and 3.5 inches per year, respectively. Primary model calibration was completed using statistical techniques through parameter estimation using the parameter estimation suite of software with Tikhonov regularization. Calibration targets for the groundwater model included 343,067 groundwater levels measured in wells and 10,820 estimated monthly stream base flows at streamgages. A total of 1,312 parameters were adjusted during calibration to improve the match between calibration targets and simulated equivalents. Comparison of calibration targets to simulated equivalents indicated that, at the regional scale, the model correctly reproduced groundwater levels and stream base flows for 1940–2009. This comparison indicates that the model can be used to examine the likely response of the aquifer system to potential future stresses. Mean calibrated recharge for 1940–49 and 2000–2009 was smaller than that estimated with the soil-water-balance model. This indicated that although the general spatial patterns of recharge estimated with the soil-water-balance model were approximately correct at the regional scale of the Northern High Plains aquifer, the soil-water-balance model had overestimated recharge, and adjustments were needed to decrease recharge to improve the match of the groundwater model to calibration targets. The largest components of the simulated groundwater budgets were recharge from precipitation, recharge from canal seepage, outflows to evapotranspiration, and outflows to stream base flow. Simulated outflows to irrigation wells increased from 7 percent of total outflows in 1940–49 to 38 percent of 1970–79 total outflows and 49 percent of 2000–2009 total outflows.

Colorado, Kansas, Nebraska, South Dakota, Wyoming↗

USGS2015

Numerical simulation of groundwater flow, resource optimization, and potential effects of prolonged drought for the Citizen Potawatomi Nation Tribal Jurisdictional Area, central Oklahoma

A hydrogeological study including two numerical groundwater-flow models was completed for the Citizen Potawatomi Nation Tribal Jurisdictional Area of central Oklahoma. One numerical groundwater-flow model, the Citizen Potawatomi Nation model, encompassed the jurisdictional area and was based on the results of a regional-scale hydrogeological study and numerical groundwater flow model of the Central Oklahoma aquifer, which had a geographic extent that included the Citizen Potawatomi Nation Tribal Jurisdictional Area. The Citizen Potawatomi Nation numerical groundwater-flow model included alluvial aquifers not in the original model and improved calibration using automated parameter-estimation techniques. The Citizen Potawatomi Nation numerical groundwater-flow model was used to analyze the groundwater-flow system and the effects of drought on the volume of groundwater in storage and streamflow in the North Canadian River. A more detailed, local-scale inset model was constructed from the Citizen Potawatomi Nation model to estimate available groundwater resources for two Citizen Potawatomi Nation economic development zones near the North Canadian River, the geothermal supply area and the Iron Horse Industrial Park. Groundwater pumping rates at potential well locations were optimized using the most recent version of the U.S. Geological Survey Groundwater-Management Process for MODFLOW. The objectives of optimization were to determine if a total pumping rate of 500 gallons per minute could be pumped from 5 wells at the geothermal supply area and to maximize discharge from 16 wells at the Iron Horse Industrial Park without exceeding specified head drawdown constraints at the pumping wells and thus prevent groundwater depletion. The inset model was used to estimate North Canadian River streamflow depletion caused by optimized pumping at the Iron Horse Industrial Park because water quality was a concern, and the river may have degraded water quality compared to water in other parts of the alluvial aquifer. The fate of streamflow that infiltrates into groundwater because of pumping was not directly determined, but it was assumed that this water could end up in the well discharge, and was considered to be a maximum proportion of well discharge derived from the North Canadian River. The total optimized continuous pumping rate from five managed wells at the geothermal supply area was 638 gallons per minute, which exceeded the target pumping rate of 500 gallons per minute. The total continuous pumping rate from 16 wells at the Iron Horse Industrial Park was 1,472 gallons per minute, which induced stream infiltration of approximately 4.1 gallons per minute (approximately 0.3 percent of the total well discharge) from the North Canadian River. To estimate the effects of drought on water resources in the Citizen Potawatomi Nation Tribal Jurisdictional Area, a hypothetical 10-year drought during which precipitation would decrease by 50 percent was simulated by decreasing model groundwater recharge by the same proportion for the period 1990–2000 of the transient model. The effects of the drought were estimated by calculating the change in the volume of groundwater storage and groundwater flow to streams at the end of the drought period, and the change in simulated streamflow in the North Canadian River at the streamflow-gaging station at Shawnee, Okla., during and after the drought. The hypothetical decrease in recharge during the simulated drought caused groundwater in storage over the entire model in the study area to decrease by 361,500 acre-feet (14,100 acre-feet in the North Canadian River alluvial aquifer and 347,400 acre-feet in the Central Oklahoma aquifer), or approximately 0.2 percent of the total groundwater in storage over the drought period. This small percentage of groundwater loss showed that the Central Oklahoma aquifer as a bedrock aquifer has relatively low rates of recharge from the surface relative to the approximate storage. The budget for base flow to the North Canadian River indicated that the change in groundwater flow to the North Canadian River decreased during the 10-year drought by 386,500 acre-feet, or 37 percent. In all other parts of the Citizen Potawatomi Nation Tribal Jurisdictional Area, base flow decreased by 292,000 acre-feet, or 28 percent. Streamflow in the North Canadian River at the streamflow-gaging station at Shawnee, Okla., decreased during the hypothetical drought by as much as 28 percent, and the mean change in streamflow decreased as much as 16 percent. Streamflow at the Shawnee streamflow-gaging station did not recover to nondrought conditions until about 3 years after the simulated drought ended, during the relatively wet year of 2007.

Oklahoma↗