Geology topics

John H. Ellis

Publications and source records attributed to John H. Ellis.

7 recordsLinked to original sources

USGS2024

Use of a numerical groundwater-flow model and projected climate scenarios to simulate the effects of future climate conditions on base flow for reach 1 of the Washita River alluvial aquifer and Foss Reservoir storage, western Oklahoma

To better understand the relation between climate variability and future groundwater resources in reach 1 of the Washita River alluvial aquifer and Foss Reservoir in western Oklahoma, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation, used a previously published numerical groundwater-flow model and climate-model data to investigate changes in base flow and reservoir storage by evaluating three scenarios. The three projected climate scenarios were (1) a central-tendency scenario, (2) a warmer/drier scenario, and (3) a less-warm/wetter scenario. To estimate future base flow and groundwater availability in western Oklahoma, specifically in reach 1 of the Washita River alluvial aquifer, downscaled climate-model data from 231 Coupled Model Intercomparison Project phase 5 (CMIP5) projections coupled with a previously published numerical groundwater-flow model were used to compare the effects of different climate scenarios on the aquifer. Changes in base flow and groundwater-level elevations during a 30-year baseline scenario (1985–2014) and the three 30-year projected climate scenarios (2050–79) under central-tendency, warmer/drier, and less-warm/wetter climatic conditions were assessed by using the calibrated model. In the simulations, the amount of base flow and reservoir storage declined in the central-tendency and warmer/drier scenarios compared to the amount of base flow and reservoir storage under historical climatic conditions (baseline scenario). Mean annual change in reservoir storage decreased from the baseline scenario the most in the warmer/drier scenario, followed by the central-tendency scenario, but increased in the less-warm/wetter scenario compared to the baseline scenario. At the end of the simulation period (2079), the largest magnitude differences in groundwater-level elevations in all three projected climate scenarios relative to the baseline scenario occurred upstream from Foss Reservoir. Results from incorporating downscaled climate projections into localized numerical groundwater-flow models can highlight potential future changes in and implications for groundwater resources and availability.

Oklahoma↗

USGS2022

Treatment of the Chicot and Evangeline aquifers as a single hydrogeologic unit and use of geostatistical interpolation methods to develop gridded surfaces of water-level altitudes and water-level changes in the Chicot and Evangeline aquifers (undifferentiated) and Jasper aquifer, greater Houston area, Texas, 2021

The greater Houston area of Texas includes approximately 11,000 square miles and encompasses all or part of 11 counties (Harris, Galveston, Fort Bend, Montgomery, Brazoria, Chambers, Grimes, Liberty, San Jacinto, Walker, and Waller). From the early 1900s until the mid-1970s, groundwater withdrawn from the three primary aquifers that compose the Gulf Coast aquifer system—the Chicot, Evangeline, and Jasper aquifers—had been the primary source of water for the greater Houston area. The withdrawal of groundwater was unregulated prior to 1975, resulting in land-surface subsidence caused by large water-level declines in the greater Houston area. This report, prepared by the U.S. Geological Survey in cooperation with the Harris-Galveston Subsidence District, City of Houston, Fort Bend Subsidence District, Lone Star Groundwater Conservation District, and Brazoria County Groundwater Conservation District, describes updates to the ways in which water-level altitudes and water-level changes in the greater Houston area are presented relative to previous U.S. Geological Survey reports. The first update involves presenting water-level altitudes and water-level changes as a combined (undifferentiated) representation of the Chicot and Evangeline aquifers. The second update concerns the methods used to depict water-level altitudes and water-level changes in the greater Houston area in interpretive reports, with geostatistical interpolation methods replacing manual contouring methods. The Chicot and Evangeline aquifers have historically been described as distinct hydrogeologic units for the purpose of water-level mapping. A confining unit does not separate these two aquifers in the study area, and water-level data from colocated wells screened in these aquifers indicate that there is likely a substantial degree of hydrogeologic connection. From a groundwater-flow perspective, these two aquifer units predominantly function as a single unit. Hence, the decision was made to combine the Chicot and Evangeline aquifers into a single, undifferentiated hydrogeologic unit for the purposes of assessing water-level altitudes and water-level changes over time. The 2020 water-level altitudes for the Chicot, Evangeline, and Jasper aquifers were re-created in this report from computer algorithms of the contoured datasets as gridded surfaces to demonstrate the similarity of results from geostatistical interpolation methods to those from manual contouring methods.

Texas↗

USGS2022

Potential effects of out-of-basin groundwater transfers on spring discharge, base flow, and groundwater storage pertaining to the Rush Springs aquifer in and near the Caddo Nation of Oklahoma Tribal jurisdictional area, western Oklahoma

The U.S. Geological Survey (USGS), in cooperation with the Caddo Nation of Oklahoma and Bureau of Indian Affairs, assessed four groundwater-withdrawal scenarios and their potential effects on the Rush Springs aquifer in and near the Caddo Nation of Oklahoma Tribal jurisdictional area in western Oklahoma. Increases in industrial and public water supply needs have led to increased development of water resources within the Rush Springs aquifer. As new areas within the aquifer are developed, increased water withdrawals may result in decreases in available groundwater resources and conflicts among water users. For this study, a previously published numerical groundwater-flow model of the Rush Springs aquifer was modified to simulate the potential effects of four groundwater withdrawal scenarios. For the previously published calibrated model, groundwater flow was simulated from 1979 through 2015. In this study, groundwater flow simulations were extended through 2035. The period from 2016 through 2035 is referred to as the “20-year projection.” Four groundwater withdrawal scenarios starting in 2007 and continuing through 2035 were evaluated. Scenario 1 simulated no groundwater withdrawals; scenario 2 simulated no withdrawals allocated for out-of-basin water-use transfers; scenario 3 simulated withdrawals based on reported withdrawals during the 2007–15 simulation period and compounded annual increases in groundwater use during the subsequent 20-year projection; and scenario 4 simulated maximum permitted withdrawals for allocation to out-of-basin water-use transfers. Out-of-basin water transfers were classified as withdrawals that are not returned back to the aquifer. At the springs of interest, changes in water-level altitudes in response to different groundwater withdrawal scenarios were simulated by comparing the results from different model cells. Between 2007 and 2015, scenarios 2–4 yielded similar simulated water-level altitudes in the model cells containing springs of interest, with water-level altitudes decreasing to below the land surface altitude at 13 of the total 25 springs of interest, whereas under scenario 1 there were only two model cells containing springs of interest where the simulated water-level altitudes of a spring decreased to below land surface altitude. For the 20-year projection, water-level altitudes at springs simulated in model cells in scenarios 2–4 decreased to below land surface altitude for 13 of the total 25 model cells containing springs of interest, whereas under scenario 1 there were only two model cells containing springs of interest where the simulated water-level altitudes of a spring decreased to below land surface altitude. The potential effects of groundwater withdrawals were evaluated by comparing changes in groundwater storage between the four scenarios. The 2007–15 groundwater withdrawal scenarios were used to simulate the potential effects of groundwater withdrawal rates on groundwater storage of the Rush Springs aquifer. The simulated groundwater storage change in the Rush Springs aquifer ranged from an increase of 2.8 percent for scenario 1 to an increase of 1.0 percent for scenario 4. Projected 20-year groundwater withdrawal scenarios were used to simulate the potential effects of selected groundwater withdrawal rates on groundwater storage of the Rush Springs aquifer. Simulated groundwater storage changes ranged from a decrease of 0.5 percent for scenario 1 to a decrease of 0.7 percent for scenario 4.

Oklahoma↗

USGS2021

Hydrogeology and model-simulated groundwater availability in the Salt Fork Red River aquifer, southwestern Oklahoma, 1980–2015

The 1973 Oklahoma Water Law (82 OK Stat § 82-1020.5) requires that the Oklahoma Water Resources Board (OWRB) conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin (hereinafter referred to as an “aquifer”). The maximum annual yield allocated per acre of land is known as the equal-proportionate-share (EPS) pumping rate. At present (2021), the OWRB has not yet established a maximum annual yield and EPS pumping rate for the Salt Fork Red River aquifer. To provide updated information to the OWRB that could support evaluation and determination of an appropriate maximum annual yield, the U.S. Geological Survey (USGS), in cooperation with the OWRB, conducted a hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater availability in the Salt Fork Red River aquifer. The Salt Fork Red River aquifer in Greer, Harmon, and Jackson Counties of southwestern Oklahoma is composed of about 274.5 square miles of alluvium and terrace deposits associated with the Salt Fork Red River. The mean annual recharge rate to the Salt Fork Red River aquifer for the period 1980–2015 was estimated to be about 2.94 inches per year, or 10.0 percent of the mean annual precipitation for the same period (29.4 inches per year). This 1980–2015 mean annual recharge rate is equivalent to a mean annual recharge rate of about 38,000 acre-feet per year (acre-ft/yr) for the Salt Fork Red River aquifer excluding about 19,764 acres comprising the Mulberry Creek and Horse Creek terraces. The mean annual recharge rates upgradient and downgradient from USGS streamgage 07300500 Salt Fork Red River at Mangum, Okla. (hereinafter referred to as the “Mangum gage”), apportioned by aquifer area (41.5 and 58.5 percent, respectively), were about 16,000 and 22,000 acre-ft/yr, respectively. Mean annual groundwater use for the study period (1980–2015) was 3,532.7 acre-ft/yr; about 77 percent of that groundwater use was for irrigation, and about 23 percent was for public supply. Most groundwater use for irrigation was associated with wells in the Martha terrace. A hydrogeologic framework was developed for the Salt Fork Red River aquifer and included a definition of the aquifer extent and potentiometric surface, as well as a description of the textural and hydraulic properties of aquifer materials. The hydrogeologic framework was used in the construction of the numerical groundwater-flow model of the Salt Fork Red River aquifer described in this report. A conceptual model for the Salt Fork Red River aquifer that reasonably represents the groundwater-flow system was developed to constrain the construction and calibration of the numerical model. The conceptual-model water budget estimated mean annual inflows to, and outflows from, the Salt Fork Red River aquifer for the period 1980–2015 and included a subaccounting of mean annual inflows and outflows for the portions of the aquifer that were upgradient and downgradient from the Mangum gage. The numerical groundwater-flow model of the Salt Fork Red River aquifer was constructed by using MODFLOW-2005 with the Newton formulation solver. The model of the Salt Fork Red River aquifer was spatially discretized into 1,050 rows, 1,125 columns, about 170,000 active cells measuring 200 by 200 feet (ft), and a single convertible layer. The model was temporally discretized into 432 monthly transient stress periods (each with two time steps to improve model stability). An initial steady-state stress period 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 as well as some groundwater-availability scenarios. The model was calibrated to water-table-altitude observations at selected wells and base-flow observations at selected streamgages. The simulated saturated thickness of the Salt Fork Red River aquifer was determined by subtracting the altitude of the aquifer base from the simulated water-table altitude at the end of the numerical-model period (2015). The simulated saturated thickness was more than 75 ft in a paleochannel in the Dodson terrace near the Texas border. The mean aquifer thickness (sum of saturated and unsaturated) was 49.62 ft, and the mean saturated thickness was 28.55 ft. A simulated mean transmissivity of 1,024 feet squared per day was computed from the calibrated hydraulic conductivity and saturated thickness of each cell. The simulated available water in storage at the end of the numerical-model period (2015) was 526,117 acre-feet (acre-ft); about 42 percent of that total was available upgradient from the Mangum gage, and about 58 percent of that total was available downgradient from the Mangum gage (including the Mangum terrace). Three types of groundwater-availability scenarios were run using the calibrated numerical model. These scenarios were used to (1) estimate the EPS pumping rate that ensures a minimum 20-, 40-, and 50-year life of the aquifer, (2) quantify the potential effects of projected well withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical 10-year drought on base flow and groundwater storage. The 20-, 40-, and 50-year EPS pumping rates under normal recharge conditions were about 0.51, 0.48, and 0.48 acre-foot per acre per year, respectively. Given the 155,929-acre modeled aquifer area, these rates correspond to annual yields of about 78,800, 74,900, and 74,700 acre-ft/yr, respectively. For the 20-year EPS scenario, decreasing and increasing recharge by 10 percent resulted in a 6-percent change in the EPS pumping rate in both cases; for the 40- and 50-year EPS scenarios, decreasing and increasing recharge by 10 percent resulted in a 7-percent change in the EPS pumping rate in both cases. Projected 50-year pumping scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage of the Salt Fork Red River aquifer and base flows in the Salt Fork Red River. The effects of well withdrawals were evaluated by quantifying differences in groundwater storage and base flow in four 50-year scenarios, which applied (1) no groundwater pumping, (2) mean pumping rates for the study period (1980–2015), (3) 2015 pumping rates, and (4) increasing demand pumping rates at simulated wells. The increasing demand pumping rates assumed a cumulative 20.4-percent increase in pumping over 50 years based on 2010–60 demand projections for southwestern Oklahoma. Groundwater storage after 50 years with no pumping was 535,000 acre-ft, or 8,900 acre-ft (1.7 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean water-table-altitude increase of 0.48 ft. Groundwater storage after 50 years of pumping at the mean rate for the study period (1980–2015) was 519,900 acre-ft, or 6,200 acre-ft (1.2 percent) less than the initial groundwater storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.34 ft. Groundwater storage at the end of the 50-year period with 2015 pumping rates was 513,100 acre-ft, or 13,000 acre-ft (2.5 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.71 ft. Groundwater storage at the end of the 50-year period with increasing demand pumping rates was 509,700 acre-ft, or 16,500 acre-ft (3.1 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.89 ft. A hypothetical 10-year drought scenario was used to simulate the effects of a prolonged period of reduced recharge on groundwater storage. The period January 1983–December 1992 was chosen as the simulated drought period. Drought effects were quantified by comparing the results of the drought scenario to those of the calibrated numerical model (no drought) at the end of the simulated drought period (1992). To simulate the hypothetical drought, recharge in the calibrated numerical model was reduced by 50 percent during the simulated drought period (1983–92). Upstream inflows from the Salt Fork Red River, Turkey Creek, and Bitter Creek were reduced by 75 percent. Groundwater storage at the end of the drought period (1992) was 479,200 acre-ft, or 53,200 acre-ft (10.0 percent) less than the groundwater storage of the calibrated numerical model at the end of the drought period. This decrease in groundwater storage is equivalent to a mean water-table-altitude decline of 2.9 ft. At the end of the 10-year hypothetical drought period, simulated base flows at the Mangum gage and USGS streamgage 07301110 Salt Fork Red River near Elmer, Okla., had decreased by about 80 and 70 percent, respectively.

Oklahoma↗

USGS2020

Hydrogeology, numerical simulation of groundwater flow, and effects of future water use and drought for reach 1 of the Washita River alluvial aquifer, Roger Mills and Custer Counties, western Oklahoma, 1980–2015

The Washita River alluvial aquifer is a valley-fill and terrace alluvial aquifer along the valley of the Washita River in western Oklahoma that provides a productive source of groundwater for agricultural irrigation and water supply. The Oklahoma Water Resources Board (OWRB) has designated the westernmost section of the aquifer in Roger Mills and Custer Counties, Okla., as reach 1 of the Washita River alluvial aquifer; reach 1 is the focus of this report. The OWRB issued an order on November 13, 1990, that established the maximum annual yield (MAY; 120,320 acre-feet per year [acre-ft/yr]) and equal-proportionate-share (EPS) pumping rate (2.0 acre-feet per acre per year [(acre-ft/acre)/yr]) for reach 1 of the Washita River alluvial aquifer. The MAY and EPS were based on hydrologic investigations that evaluated the effects of potential groundwater withdrawals on groundwater availability in the Washita River alluvial aquifer. Every 20 years, the OWRB is statutorily required to update the hydrologic investigation on which the MAY and EPS were based. Because 30 years have elapsed since the last order was issued, the U.S. Geological Survey, in cooperation with the OWRB, conducted a new hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater flow and availability in the Washita River alluvial aquifer. The Washita River is the primary source of inflow to Foss Reservoir, a Bureau of Reclamation reservoir constructed in 1961 for flood control, water supply, and recreation. Foss Reservoir provides water for Bessie, Clinton, New Cordell, and Hobart, Okla. Nearly 98 percent of the total groundwater use from the Washita River alluvial aquifer during 1967 to 2015 was for irrigation; other uses of groundwater in the study area include public supply, mining, and agriculture. A hydrogeologic framework was developed for the Washita River alluvial aquifer and included the physical characteristics of the aquifer, the geologic setting, the hydraulic properties of hydrogeologic units, the potentiometric surface (water table), and groundwater-flow directions at a scale that captures the regional controls on groundwater flow. The Washita River alluvial aquifer consists of alluvium and terrace deposits that were transported primarily by water and range from clay to gravel in size. The terrace includes windblown deposits of silt size and, in some cases, contains gravel laid down at several levels along former courses of present-day rivers. A conceptual flow model is a simplified description of the aquifer system that includes hydrologic boundaries, major inflow and outflow sources of the groundwater-flow system, and a conceptual water budget with the estimated mean flows between those hydrologic boundaries. During the study period 1980–2015, mean annual groundwater withdrawals, predominantly used for agricultural irrigation, totaled 5,502 acre-ft/yr, or 14 percent of aquifer outflows. When applied across the 132-square-mile aquifer area used for modeling purposes (84,366 acres), mean annual recharge of 3.15 inches per year corresponds to a mean annual recharge volume of 22,169 acre-ft/yr, or 56 percent of aquifer inflows. The annual saturated-zone evapotranspiration outflow was 11,828 acre-ft/yr for the Washita River alluvial aquifer, or about 30 percent of aquifer outflows. For the Washita River alluvial aquifer, lateral flow was 17,157 acre-ft/yr, or 44 percent of the aquifer inflows. The conceptual flow model and hydrogeologic framework were used to conceptualize, design, and build the numerical groundwater-flow model. A numerical groundwater-flow model of the Washita River alluvial aquifer was constructed by using MODFLOW-2005. The Washita River alluvial aquifer groundwater-model grid was spatially discretized into 350-foot (ft) cells and two layers. Layer 1 represented the undifferentiated alluvium and terrace deposits of Quaternary age, and layer 2 represented the bedrock of Permian age, which was given a uniform nominal thickness of 100 ft. The groundwater-simulation period was temporally discretized into 433 monthly transient stress periods, representing January 1980 to December 2015. An initial 365-day steady-state stress period was configured to represent mean annual inflows and outflows from the Washita River alluvial aquifer for the study period. The groundwater-flow model was calibrated manually and by automated adjustment of model inputs by using PEST++. Calibration targets for the Washita River alluvial aquifer model included groundwater-level observations and reservoir-stage observations, as well as base-flow and stream-seepage estimates. Three groundwater-availability scenarios were used in the calibrated groundwater model to (1) estimate the EPS pumping rate that retains the saturated thickness that meets the minimum 20-year life of the aquifer, (2) quantify the effects of projected pumping rates on groundwater storage over a 50-year period, and (3) evaluate how projected pumping rates extended 50 years into the future and sustained hypothetical drought conditions over a 10-year period affect base flow and groundwater in storage. The results of the groundwater-availability scenarios could be used by the OWRB to reevaluate the established MAY of groundwater from the Washita River alluvial aquifer. EPS scenarios for the Washita River alluvial aquifer were run for periods of 20, 40, and 50 years. The 20-, 40-, and 50-year EPS pumping rates under normal recharge conditions were 1.7, 1.6, and 1.6 (acre-ft/acre)/yr, respectively. Given the aquifer area used for modeling purposes (84,366 acres), these rates correspond to annual yields of 142,579, 134,986, and 134,986 acre-ft/yr, respectively. Groundwater storage at the end of the 20-year EPS scenario was about 281,000 acre-feet (acre-ft), or about 306,000 acre-ft (52 percent) less than the starting storage. Considering the land-surface area of the Washita River alluvial aquifer and using a specific yield of 0.12, this decrease in storage was equivalent to a mean groundwater-level decline of about 30 ft. The Washita River downstream from Foss Reservoir and most of the streams in the study area were dry at the end of the 20-year EPS scenario. Foss Reservoir stage was below the dead-pool stage of 1,597 ft after about 7 years of pumping in the 20-year EPS scenario. Four projected 50-year groundwater-use scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage in the Washita River alluvial aquifer. These four scenarios used (1) no groundwater use, (2) groundwater use at the 2015 pumping rate, (3) mean groundwater use for the simulation period, and (4) increasing groundwater use. Groundwater storage after 50 years with no groundwater use was 545,249 acre-ft, or 693 acre-ft (0.1 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean groundwater-level increase of 0.1 ft. Groundwater storage at the end of the 50-year period with 2015 pumping rates was 543,831 acre-ft, or 723 acre-ft (0.1 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean groundwater-level decrease of 0.1 ft. Groundwater storage after 50 years with the mean pumping rate for the study period was 543,202 acre-ft, or 1,349 acre-ft (0.2 percent) less than the initial groundwater storage; this groundwater storage decrease is equivalent to a mean groundwater-level decrease of 0.1 ft. Groundwater storage at the end of the 50-year period with an increasing demand groundwater-pumping rate, which was 38 percent greater than the 2015 groundwater-pumping rate, was 542,584 acre-ft, or 1,967 acre-ft (0.4 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean groundwater-level decrease of 0.2 ft. A hypothetical 10-year-drought scenario was used to simulate the effects of a prolonged period of reduced recharge on groundwater storage in the Washita River alluvial aquifer and Foss Reservoir stage and storage. To simulate the hypothetical drought, recharge in the calibrated model was reduced by 50 percent during the simulated drought period (1983–1992). Groundwater storage at the end of the drought period in December 1992 was 562,000 acre-ft, or 36,000 acre-ft (6 percent) less than the groundwater storage of the calibrated groundwater model (598,000 acre-ft). At the end of the hypothetical drought, the largest changes in saturated thickness (as great as 43.5 ft) were in the area upgradient from Foss Reservoir, particularly in the terrace at the model boundary. Substantial decreases in the Foss Reservoir stage began during the fall of 1985 in conjunction with base-flow decreases of up to 100 percent at U.S. Geological Survey streamgage 07324200 Washita River near Hammon, Okla. These lake-stage declines outpaced groundwater-level declines in the surrounding aquifer. The minimum Foss Reservoir storage simulated during the drought period was 77,954 acre-ft, which was a decrease of 46 percent from the nondrought storage.

Oklahoma↗

USGS2017

Hydrogeology and simulated groundwater flow and availability in the North Fork Red River aquifer, southwest Oklahoma, 1980–2013

On September 8, 1981, the Oklahoma Water Resources Board established regulatory limits on the maximum annual yield of groundwater (343,042 acre-feet per year) and equal-proportionate-share (EPS) pumping rate (1.0 acre-foot per acre per year) for the North Fork Red River aquifer. The maximum annual yield and EPS were based on a hydrologic investigation that used a numerical groundwater-flow model to evaluate the effects of potential groundwater withdrawals on groundwater availability in the North Fork Red River aquifer. The Oklahoma Water Resources Board is statutorily required (every 20 years) to update the hydrologic investigation on which the maximum annual yield and EPS were based. Because 20 years have elapsed since the final order was issued, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, conducted an updated hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater flow and availability in the North Fork Red River aquifer in Oklahoma. This report describes a hydrologic investigation of the North Fork Red River aquifer that includes an updated summary of the aquifer hydrogeology. As part of this investigation, groundwater flow and availability were simulated by using a numerical groundwater-flow model. The North Fork Red River aquifer in Beckham, Greer, Jackson, Kiowa, and Roger Mills Counties in Oklahoma is composed of about 777 square miles (497,582 acres) of alluvium and terrace deposits along the North Fork Red River and tributaries, including Sweetwater Creek, Elk Creek, Otter Creek, and Elm Fork Red River. The North Fork Red River is the primary source of surface-water inflow to Lake Altus, which overlies the North Fork Red River aquifer. Lake Altus is a U.S. Bureau of Reclamation reservoir with the primary purpose of supplying irrigation water to the Lugert-Altus Irrigation District. A hydrogeologic framework was developed for the North Fork Red River aquifer and included a definition of the aquifer extent and potentiometric surface, as well as a description of the textural and hydraulic properties of aquifer materials. The hydrogeologic framework was used in the construction of a numerical groundwater-flow model of the North Fork Red River aquifer described in this report. A conceptual model of aquifer inflows and outflows was developed for the North Fork Red River aquifer to constrain the construction and calibration of a numerical groundwater-flow model that reasonably represented the groundwater-flow system. The conceptual-model water budget estimated mean annual inflows to and outflows from the North Fork Red River aquifer for the period 1980–2013 and included a sub-accounting of mean annual inflows and outflows for the portions of the aquifer that were upgradient and downgradient from Lake Altus. The numerical groundwater-flow model simulated the period 1980–2013 and was calibrated to water-table-altitude observations at selected wells, monthly base flow at selected streamgages, net streambed seepage as estimated for the conceptual model, and Lake Altus stage. Groundwater-availability scenarios were performed by using the calibrated numerical groundwater-flow model to (1) estimate the EPS pumping rate that guarantees a minimum 20-, 40-, and 50-year life of the aquifer, (2) quantify the potential effects of projected well withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical (10-year) drought on base flow and groundwater storage. The results of the groundwater-availability scenarios could be used by the Oklahoma Water Resources Board to reevaluate the maximum annual yield of groundwater from the North Fork Red River aquifer. EPS scenarios for the North Fork Red River aquifer were run for periods of 20, 40, and 50 years. The 20-, 40-, and 50-year EPS pumping rates under normal recharge conditions were 0.59, 0.52, and 0.52 acre-foot per acre per year, respectively. Given the 497,582-acre aquifer area, these rates correspond to annual yields of about 294,000, 259,000, and 259,000 acre-feet per year, respectively. Groundwater storage at the end of the 20-year EPS scenario was about 951,000 acre-feet, or about 1,317,000 acre-feet (58 percent) less than the starting EPS scenario storage. This decrease in storage was equivalent to a mean water-level decline of about 22 feet. Most areas of the active alluvium near the North Fork Red River, Elk Creek, and Elm Fork Red River remained partially saturated through the end of the EPS scenario because of streambed seepage. Lake Altus storage was reduced to zero after 6–7 years of EPS pumping in each scenario. Projected 50-year pumping scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage of the North Fork Red River aquifer and base flows in the North Fork Red River upstream from Lake Altus. The effects of well withdrawals were evaluated by comparing changes in groundwater storage and base flow between four 50-year scenarios using (1) no groundwater pumping, (2) mean pumping rates for the study period (1980–2013), (3) 2013 pumping rates, and (4) increasing demand pumping rates. The increasing demand pumping rates assumed a 20.4-percent increase in pumping over 50 years based on 2010–60 demand projections for southwest Oklahoma. Groundwater storage after 50 years with no pumping was about 2,606,000 acre-feet, or 137,000 acre-feet (5.5 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean water-level increase of 2.3 feet. Groundwater storage after 50 years with the mean pumping rate for the study period (1980–2013) was about 2,476,000 acre-feet, or about 7,000 acre-feet (0.3 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean water-level increase of 0.1 foot. Groundwater storage at the end of the 50-year period with 2013 pumping rates was about 2,398,000 acre-feet, or about 70,000 acre-feet (2.8 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-level decline of 1.2 feet. Groundwater storage at the end of the 50-year period with increasing demand pumping rates was about 2,361,000 acre-feet, or about 107,000 acre-feet (4.3 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-level decline of 1.8 feet. Mean annual base flow simulated at the Carter streamgage (07301500) on North Fork Red River increased by about 4,000 acre-feet (10 percent) after 50 years with no pumping and decreased by about 5,400 acre-feet (13 percent) after 50 years with increasing demand pumping rates. Mean annual base flow simulated at the North Fork Red River inflow to Lake Altus increased by about 7,400 acre-feet (15 percent) after 50 years with no pumping and decreased by about 5,800 acre-feet (12 percent) after 50 years with increasing demand pumping rates. A hypothetical 10-year drought scenario was used to simulate the effects of a prolonged period of reduced recharge on groundwater storage and Lake Altus stage and storage. Drought effects were quantified by comparing the results of the drought scenario to those of the calibrated numerical model (no drought). To simulate the hypothetical drought, recharge in the calibrated numerical model was reduced by 50 percent during the simulated drought period (1984–1993). Groundwater storage at the end of the drought period was about 2,271,000 acre-feet, or about 426,000 acre-feet (15.8 percent) less than the groundwater storage of the calibrated numerical model. This decrease in groundwater storage is equivalent to a mean water-table-altitude decline of 7.1 feet. At the end of the 10-year hypothetical drought period, base flows at the Sweetwater (07301420), Carter (07301500), Headrick (07305000), and Snyder (07307010) streamgages had decreased by about 37, 61, 44, and 45 percent, respectively. The minimum Lake Altus storage simulated during the drought period was 403 acre-feet, which was a decline of 92 percent from the nondrought storage. Reduced base flows in the North Fork Red River were the primary cause of Lake Altus storage declines.

Oklahoma↗

USGS2017

Hydrogeology and simulation of groundwater flow and analysis of projected water use for the Canadian River alluvial aquifer, western and central Oklahoma

This report describes a study of the hydrogeology and simulation of groundwater flow for the Canadian River alluvial aquifer in western and central Oklahoma conducted by the U.S. Geological Survey in cooperation with the Oklahoma Water Resources Board. The report (1) quantifies the groundwater resources of the Canadian River alluvial aquifer by developing a conceptual model, (2) summarizes the general water quality of the Canadian River alluvial aquifer groundwater by using data collected during August and September 2013, (3) evaluates the effects of estimated equal proportionate share (EPS) on aquifer storage and streamflow for time periods of 20, 40, and 50 years into the future by using numerical groundwater-flow models, and (4) evaluates the effects of present-day groundwater pumping over a 50-year period and sustained hypothetical drought conditions over a 10-year period on stream base flow and groundwater in storage by using numerical flow models. The Canadian River alluvial aquifer is a Quaternary-age alluvial and terrace unit consisting of beds of clay, silt, sand, and fine gravel sediments unconformably overlying Tertiary-, Permian-, and Pennsylvanian-age sedimentary rocks. For groundwater-flow modeling purposes, the Canadian River was divided into Reach I, extending from the Texas border to the Canadian River at the Bridgeport, Okla., streamgage (07228500), and Reach II, extending downstream from the Canadian River at the Bridgeport, Okla., streamgage (07228500), to the confluence of the river with Eufaula Lake. The Canadian River alluvial aquifer spans multiple climate divisions, ranging from semiarid in the west to humid subtropical in the east. The average annual precipitation in the study area from 1896 to 2014 was 34.4 inches per year (in/yr). A hydrogeologic framework of the Canadian River alluvial aquifer was developed that includes the areal and vertical extent of the aquifer and the distribution, texture variability, and hydraulic properties of aquifer materials. The aquifer areal extent ranged from less than 0.2 to 8.5 miles wide. The maximum aquifer thickness was 120 feet (ft), and the average aquifer thickness was 50 ft. Average horizontal hydraulic conductivity for the Canadian River alluvial aquifer was calculated to be 39 feet per day, and the maximum horizontal hydraulic conductivity was calculated to be 100 feet per day. Recharge rates to the Canadian River alluvial aquifer were estimated by using a soil-water-balance code to estimate the spatial distribution of groundwater recharge and a water-table fluctuation method to estimate localized recharge rates. By using daily precipitation and temperature data from 39 climate stations, recharge was estimated to average 3.4 in/yr, which corresponds to 8.7 percent of precipitation as recharge for the Canadian River alluvial aquifer from 1981 to 2013. The water-table fluctuation method was used at one site where continuous water-level observation data were available to estimate the percentage of precipitation that becomes groundwater recharge. Estimated annual recharge at that site was 9.7 in/yr during 2014. Groundwater flow in the Canadian River alluvial aquifer was identified and quantified by a conceptual flow model for the period 1981–2013. Inflows to the Canadian River alluvial aquifer include recharge to the water table from precipitation, lateral flow from the surrounding bedrock, and flow from the Canadian River, whereas outflows include flow to the Canadian River (base-flow gain), evapotranspiration, and groundwater use. Total annual recharge inflows estimated by the soil-water-balance code were multiplied by the area of each reach and then averaged over the simulated period to produce an annual average of 28,919 acre-feet per year (acre-ft/yr) for Reach I and 82,006 acre-ft/yr for Reach II. Stream base flow to the Canadian River was estimated to be the largest outflow of groundwater from the aquifer, measured at four streamgages, along with evapotranspiration and groundwater use, which were relatively minor discharge components. Objectives for the numerical groundwater-flow models included simulating groundwater flow in the Canadian River alluvial aquifer from 1981 to 2013 to address groundwater use and drought scenarios, including calculation of the EPS pumping rates. The EPS for the alluvial and terrace aquifers is defined by the Oklahoma Water Resources Board as the amount of fresh water that each landowner is allowed per year per acre of owned land to maintain a saturated thickness of at least 5 ft in at least 50 percent of the overlying land of the groundwater basin for a minimum of 20 years. The groundwater-flow models were calibrated to water-table altitude observations, streamgage base flows, and base-flow gain to the Canadian River. The Reach I water-table altitude observation root-mean-square error was 6.1 ft, and 75 percent of residuals were within ±6.7 ft of observed measurements. The average simulated stream base-flow residual at the Bridgeport streamgage (07228500) was 8.8 cubic feet per second (ft 3 /s), and 75 percent of residuals were within ±30 ft 3 /s of observed measurements. Simulated base-flow gain in Reach I was 8.8 ft 3 /s lower than estimated base-flow gain. The Reach II water-table altitude observation root-mean-square error was 4 ft, and 75 percent of residuals were within ±4.3 ft of the observations. The average simulated stream base-flow residual in Reach II was between 35 and 132 ft 3 /s. The average simulated base-flow gain residual in Reach II was between 11.3 and 61.1 ft 3 /s. Several future predictive scenarios were run, including estimating the EPS pumping rate for 20-, 40-, and 50-year life of basin scenarios, determining the effects of current groundwater use over a 50-year period into the future, and evaluating the effects of a sustained drought on water availability for both reaches. The EPS pumping rate was determined to be 1.35 acre-feet per acre per year ([acre-ft/acre]/yr) in Reach I and 3.08 (acre-ft/acre)/yr in Reach II for a 20-year period. For the 40- and 50-year periods, the EPS rate was determined to be 1.34 (acre-ft/acre)/yr in Reach I and 3.08 (acre-ft/acre)/yr in Reach II. Storage changes decreased in tandem with simulated groundwater pumping and were minimal after the first 15 simulated years for Reach I and the first 8 simulated years for Reach II. Groundwater pumping at year 2013 rates for a period of 50 years resulted in a 0.2-percent decrease in groundwater-storage volumes in Reach I and a 0.6-percent decrease in the groundwater-storage volumes in Reach II. The small changes in storage are due to groundwater use by pumping, which composes a small percentage of the total groundwater-flow model budgets for Reaches I and II. A sustained drought scenario was used to evaluate the effects of a hypothetical 10-year drought on water availability. A 10-year period was chosen where the effects of drought conditions would be simulated by decreasing recharge by 75 percent. In Reach I, average simulated stream base flow at the Bridgeport streamgage (07228500) decreased by 58 percent during the hypothetical 10-year drought compared to average simulated stream base flow during the nondrought period. In Reach II, average simulated stream base flows at the Purcell streamgage (07229200) and Calvin streamgage (07231500) decreased by 64 percent and 54 percent, respectively. In Reach I, the groundwater-storage drought scenario resulted in a storage decline of 30 thousand acre-feet, or an average decline in the water table of 1.2 ft. In Reach II, the groundwater-storage drought scenario resulted in a storage decline of 71 thousand acre-feet, or an average decline in the water table of 2.0 ft.

Oklahoma↗