Search USGS⌕ Search

Geology topics

Research about Oklahoma

Source-linked reports with geographic coverage including Oklahoma.

At least 55 records · Page 3Linked to original sources

The 180-km-long Meers-Willow Fault System in the Southern Oklahoma Aulacogen: A potential U.S. mid-continent seismic hazard

We integrate new high-resolution aeromagnetic data with seismic reflection data, well logs, satellite remote sensing, and field observations to provide a regional view of buried and exposed structures in the Southern Oklahoma Aulacogen and to assess their potential for future seismicity. Trends ranging from NW−SE to ∼E−W, peaking at 330° ± 4.5° and 280° ± 3°, dominate the magnetic lineaments of the Southern Oklahoma Aulacogen, reflecting basement contacts, dikes, and faults, including a previously unmapped ∼100-km-long basement fault, which is herein referred to as the Willow fault. The fault disrupts, truncates, and vertically offsets basement-related seismic reflectors and overlying Paleozoic strata up through the Permian reflectors. Surface deformation along the trend includes fault-parallel monoclinal folds, pervasive fractures, and fracture-hosted mud dikes in Permian evaporite units. These structures indicate a Permian or post-Permian reactivation of the fault. Along-strike, the Willow fault connects to the NW-trending, seismically active Meers Fault to comprise the ∼180-km-long Meers-Willow fault system, which potentially represents a major seismic hazard along the Southern Oklahoma Aulacogen. Fault slip potential analyses of the mapped potential fault traces show that seismic hazards are elevated where faults have steeper dips. Given some uncertainty in the regional stress state, we also show that hazards along the NW−SE to E−W trending faults vary considerably within the uncertainty range. We propose that the Meers-Willow fault system originated as a Cambrian aulacogen-scale, basement-rooted fault that was later reactivated as a left-lateral strike-slip fault (with ∼40 km displacement) during the late Paleozoic Ancestral Rocky Mountain orogeny, highlighting that lateral offset accommodated a major component of deformation during the orogen.

Oklahoma↗

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↗

Oklahoma and Landsat

Oklahoma benefits from a varied landscape abundant in resources. Mountains, grasslands, reservoirs, rivers, fields, and forests offer employment and enjoyment in a State that epitomizes the transition from north to south and east to west. Wheat grows in northern Oklahoma; cotton grows in the south. Wetter deciduous forest lands in the southeast contrast with drier mesas in the northwest. Among the many lakes and reservoirs, the Great Salt Plains Lake on the Arkansas River has thousands of acres of salt flats with unique hourglass-shaped selenite crystals. Grassland and grazing permeate the Sooner State, which is the second largest cattle producer in the country. Oil, natural gas, and coal have long played a key role in Oklahoma’s economy and energy industry, whereas wind farms in western Oklahoma are a more recent addition. The Landsat series of Earth observation satellites has allowed analysts to search for oil in Oklahoma from hundreds of miles above the Earth. The satellite program also helps Oklahoma’s agriculture industry keep an eye on the condition of crops and grazing land. Here are just a few ways Landsat has benefited Oklahoma.

Oklahoma↗

Critical thermal maximum of stream fishes including distinct populations of Smallmouth Bass

Understanding the thermal tolerances of stream fishes, including sport fishes, is important for assessing thermal stressors that are common across the landscape. Our study objectives were to determine the thermal tolerances of 17 stream fishes (15 species and 2 genetically distinct populations of juvenile Smallmouth Bass Micropterus dolomieu : the Neosho subspecies M. dolomieu velox and the Ouachita strain M . sp. cf. dolomieu velox ). Fish were collected from the field and acclimated to laboratory conditions at 20°C or 25°C, with dissolved oxygen maintained above 6 mg/L. We determined the critical thermal maximum (CTM) using an incomplete block design with 9–11 replications for each species. During trials, we increased the water temperature at a rate of 2°C per hour until fish experienced loss of equilibrium. The estimated CTM ranged from 32.43°C to 38.26°C among species. The CTM values differed significantly between taxonomic groups and species, including the genetically distinct populations of Smallmouth Bass. The Neosho subspecies of Smallmouth Bass had a significantly lower thermal tolerance than the Ouachita strain at both acclimation temperatures; however, the magnitude of the difference was about 0.5°C greater at the higher acclimation temperature. Closely related species, including the Bigeye Shiner Notropis boops and Kiamichi Shiner N . ortenburgeri , had significantly different thermal tolerances despite occupying similar riverine locations. Our results suggest that our perceptions of a species’ thermal tolerance based on that of closely related species or that of species using similar habitat may be incorrect. Moreover, the differences in thermal tolerances among populations may be an important consideration for conservation and management actions, such as stocking decisions. Laboratory data such as those provided in this study can be integrated with field data to better assess thermal responses of fishes in a changing environment.

Oklahoma↗

Germanium redistribution during weathering of Zn mine wastes: Implications for environmental mobility and recovery of a critical mineral

Germanium (Ge) is a metal used in emerging energy technologies, communications, and defense, and has been deemed critical by the United States due to its essential applications and scarce supply. Germanium is recovered as a byproduct of zinc (Zn) sulfides, and mining and processing of these materials lead to waste that could act both as a source of extractable Ge and a source for exposure to humans and ecosystems. Yet the distribution, speciation, and mineral hosts of Ge in mining-impacted areas are poorly understood. The Tar Creek Superfund Site, a former Zn mining area and Ge producer, is a natural laboratory to understand the environmental behavior and economic implications of Ge in mine wastes. We studied the distribution and behavior of Ge in solid wastes at the Tar Creek Superfund Site using bulk and microanalytical techniques. In wastes at this site we find that Ge has been redistributed from its original host, sphalerite (ZnS), to the fine-grained weathering product hemimorphite (Zn4Si2O7(OH)2·H2O), which impacts germanium's mobility, bioaccessibility, and potential for recovery. We provide chemical and mineralogical evidence of this redistribution, along with an evaluation of the oxidation state and molecular-scale substitution of Ge into sphalerite, hemimorphite, and quartz. Geochemical modeling shows that hemimorphite is more stable than sphalerite in waste piles and provides a stable secondary repository for Ge. However, hemimorphite is fine-grained, and if ingested or inhaled is readily soluble, with the potential to release Ge. Lastly, we discuss other sites internationally where similar behavior may be important. This study shows that weathering can have a significant impact on the distribution, speciation, and mineral hosts of Ge in mine wastes; directly influence mobilization from waste piles and subsequent availability to humans and ecosystems; and dictate metallurgical strategies to target Ge for recovery.

Oklahoma↗

Mean annual runoff and annual runoff variability map for Oklahoma, 1940–2007

Hydrologic records used to create previously published maps depicting mean annual runoff are biased to a relatively dry period in Oklahoma history that was dominated by droughts. Therefore, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, developed an updated mean annual runoff and annual runoff variability map for Oklahoma and parts of adjacent States. The updated map, which is based on mean-annual-streamflow regression equations developed from available streamgage data through 2007, is assumed to be representative of the long-term mean annual runoff conditions. The map covers all 69 8-digit hydrologic units with at least 1 square mile of area in Oklahoma; those 8-digit hydrologic units contain 2,870 12-digit hydrologic units that provided the geographic framework for the analysis described in this report. Although parts of adjacent States are included in the study area, this report is primarily focused on providing a map of mean annual runoff and annual runoff variability for Oklahoma. The mean annual runoff increased from less than 0.25 inch per year in the Panhandle of northwestern Oklahoma to more than 30 inches per year in the mountainous terrain of southeastern Oklahoma. The orientation and pattern of mean annual runoff contours in this report were comparable to those of previously published map reports. The annual runoff variability, or the difference between the 80-percent and 20-percent streamflow-duration statistics, increased from less than 0.25 inch per year in the Panhandle of northwestern Oklahoma to more than 40 inches per year in the mountainous terrain of southeastern Oklahoma. The annual runoff variability data were similar in orientation and pattern to the mean annual runoff contours; annual runoff variability generally increased proportionally with increasing mean annual runoff. The annual runoff variability was also greatest, therefore, in the mountainous terrain of southeastern Oklahoma. The mean annual runoff and annual runoff variability were calculated at sampled points representing the outlets of 12-digit hydrologic units, so the map in this report is most representative of runoff conditions in rural, unregulated drainage basins at the 12-digit hydrologic-unit scale. The map was developed by using regression equations formulated on streamgage data for the entire period of record through 2007, but those equations are biased to the period 1940–2007 when streamgages became more numerous and distributed across Oklahoma. Therefore, the map is likely most representative of runoff conditions during the period 1940–2007. Because runoff is a function of climate variables that can change over time, caution is warranted when using the information in this report to project mean annual runoff and annual runoff variability conditions beyond 2007.

Oklahoma↗

Cimarron River alluvial aquifer hydrogeologic framework, water budget, and implications for future water availability in the Pawnee Nation Tribal jurisdictional area, Payne County, Oklahoma, 2016–18

The Cimarron River is a free-flowing river and is a major source of water as it flows across Oklahoma. Increased demand for water resources within the Cimarron River alluvial aquifer in north-central Oklahoma (primarily in Payne County) has led to increases in groundwater withdrawals for agriculture, public, irrigation, industrial, and domestic supply purposes. The Pawnee Nation of Oklahoma (Pawnee Nation) is particularly concerned about the sustainability of the Cimarron River alluvial aquifer and whether the aquifer will continue to be a viable water resource for future generations of Tribal members and residents. To better understand current (2021) water resources and possible future water availability in the Pawnee Nation Tribal jurisdictional area, the U.S. Geological Survey, in cooperation with the Bureau of Indian Affairs and the Pawnee Nation of Oklahoma, compiled available hydrogeologic data and developed conceptual and numerical groundwater-flow models for the Cimarron River alluvial aquifer in Payne County, north-central Oklahoma, including a focus area in the Pawnee Nation Tribal jurisdictional area for the 2016–18 study period. A conceptual water budget was created to establish estimates of groundwater fluxes into and out of the aquifer through hydrologic boundaries and groundwater withdrawals for use in the numerical groundwater-flow model. The conceptual water budget focuses on the alluvial aquifer, meaning that inflows include sources of water to the aquifer and that outflows include sources of water out of the aquifer, such as base-flow contributions to the Cimarron River. The conceptual water budget was constructed by using data from 2017 (the most complete year of record for each data type included in the model) for the Pawnee Nation subdomain of the Cimarron River alluvial aquifer model extent (Pawnee Nation subdomain). Groundwater withdrawals were estimated from groundwater-withdrawal rate information for permanent and temporary permitted wells that was obtained from the Oklahoma Water Resources Board. One-half of each annual permitted groundwater-withdrawal rate allotted was used as the estimated annual groundwater-withdrawal amount. Halving the permitted groundwater-withdrawal rate was done because permitted withdrawal rates are the maximum permitted rate and actual groundwater withdrawals are generally appreciably lower than the maximum permitted rate. Total groundwater withdrawals were estimated as 1,300 acre-feet per year for the Pawnee Nation subdomain. Various hydrogeologic data were measured to assist with model development, including depth to bedrock and water-table altitude data. In support of the model development, analyses pertaining to groundwater flow, groundwater/surface-water interactions, base flows in the Cimarron River, and lithological interpretations in the Pawnee Nation Tribal jurisdictional area were used to compute a conceptual water budget applicable to the 2016–18 study period. A numerical groundwater-flow model was developed using the hydrogeologic framework of the Cimarron River alluvial aquifer and the conceptual water budget. The numerical model consists of a single layer representing alluvium and terrace deposits within the alluvial aquifer model area. Hydraulic conductivities were estimated and modeled for the alluvium and terrace deposits in the alluvial aquifer. Base-flow values were estimated using the base-flow index from streamflow data collected at U.S. Geological Survey streamgages. Stream seepage values were derived from the mean 2017 base-flow index between certain streamgages. Hydraulic conductivities were specified an initial (before calibration) value of 120 feet per day for the alluvium deposits and 16 feet per day for the terrace deposits. The simulated inflows in the numerical groundwater-flow model of the Pawnee Nation subdomain were higher than the inflows of conceptual water budget, and the simulated outflows were lower than the outflows of the conceptual water budget. Overall, simulated base flows matched closely to observed base flows for the 2016 and 2017 stress periods. Simulated streamflow tended to match better with the observed streamflow for 2017, which was the period with the most data for the Cimarron River alluvial aquifer model. Streamflow capture analysis was applied to the steady-state simulation to identify areas of the aquifer where base flows in the Cimarron River were most sensitive to groundwater withdrawals. The initial base-flow value was assigned the value obtained from streamflow-routing software used to simulate stream outflow for the calibrated steady-state base model. Subsequent simulations were run in each active cell in the Pawnee Nation subdomain for a specified groundwater-withdrawal rate of 180,000 cubic feet per day. The study area that includes the Pawnee Nation subdomain is in the upper Arkansas River Basin. A groundwater-withdrawal rate of 180,000 cubic feet per second per day represents a 34 percent increase compared to the highest permitted groundwater-withdrawal rate for the study area, which corresponds to the estimated 34 percent increase in groundwater withdrawals predicted by 2060 for the upper Arkansas River Basin. Simulated streamflow capture was highest in the alluvium deposits adjacent to the Cimarron River; that is, base flow in the Cimarron River decreased the most for simulated groundwater withdrawals in the alluvium deposits adjacent to the Cimarron River. Streamflow capture increased as the distance of a well from the Cimarron River decreased in the simulation. The northeastern part of the Pawnee Nation subdomain showed greater streamflow capture in a broader area; streamflow in that part of the Pawnee Nation subdomain is likely more sensitive to groundwater withdrawals compared to other parts of the Pawnee Nation subdomain.

Oklahoma↗

The sensitivity of a unionid mussel (Lampsilis siliquoidea) to a permitted effluent and elevated potassium in the effluent

Freshwater mussels are one of the most imperiled groups of animals in the world and are among the most sensitive species to a variety of chemicals. However, little is known about the sensitivity of freshwater mussels to wastewater effluents. The objectives of the present study were to (1) assess the toxicity of a permitted effluent, which entered the Deep Fork River, Oklahoma (USA), to a unionid mussel ( Lampsilis siliquoidea ) and to two standard test species (cladoceran Ceriodaphnia dubia ; and fathead minnow Pimephales promelas ) in short-term 7-day effluent tests; (2) evaluate the relative sensitivities of the three species to potassium (K), an elevated major ion in the effluent, using 7-day toxicity tests with KCl spiked into a Deep Fork River upstream reference water; (3) determine the potential influences of background water characteristics on the acute K toxicity to the mussel (96-h exposures) and cladoceran (48-h exposure) in four reconstituted waters that mimicked the hardness and ionic composition ranges of the Deep Fork River; and (4) determine the potential influence of temperature on acute K toxicity to the mussel. The effluent was found to be toxic to mussels and cladocerans, and it contained elevated concentrations of major cations and anions relative to the upstream Deep Fork River reference water. The K concentration in the effluent was 48-fold greater than in the upstream water. Compared with the standard species, the mussel was more than 4-fold more sensitive to the effluent in the 7-day effluent tests and more than 8-fold more sensitive to K in the 7-day K toxicity tests. The acute K toxicity to the mussel decreased by a factor of 2 when the water hardness was increased from soft (42 mg/L as CaCO 3 ) to very hard (314 mg/L as CaCO 3 ), whereas the acute K toxicity to the cladoceran remained almost the same as hardness increased from 84 to 307 mg/L as CaCO 3 . Acute K toxicity to the mussel at 23 °C was similar to the toxicity at an elevated temperature of 28 °C. The overall results indicate that the two standard test species may not represent the sensitivity of the tested mussel to both the effluent and K, and the toxicity of K was influenced by the hardness in test waters, but by a limited magnitude.

Oklahoma↗

Shear-wave velocity site characterization in Oklahoma from joint inversion of multi-method surface seismic measurements: Implications for central U.S. Ground Motion Prediction

We analyze multimethod shear (SH)‐wave velocity ( ⁠ V S "> V S ) site characterization data acquired at three permanent and 25 temporary seismograph stations in Oklahoma that recorded M 4+ earthquakes within a 50 km hypocentral distance of at least one of the 2016 M 5.1 Fairview, M 5.8 Pawnee, or M 5.0 Cushing earthquakes to better constrain earthquake ground‐motion modeling in the region. We acquired active‐source seismic data for time‐averaged V S "> V S to 30 m depth ( V S "> V S V S 30 "> 30 ) at 28 seismograph stations near the Fairview, Pawnee, and Cushing epicentral areas. The SH‐wave refraction travel times coupled with Rayleigh‐ and Love‐wave phase velocity dispersion were extracted and modeled in a nonlinear least‐squares (L2) joint inversion to obtain a best‐fit 1D V S "> V S versus depth profile for each site. At a subset of sites where the preferred L2 inverse model did not optimally fit each of the Love, Rayleigh, and SH travel‐time datasets, we explore application of simulated annealing in a joint inversion to find a more global solution. V S "> V S V S 30 "> 30 values range from 262 to 807    m / s "> 807 m / s for the preferred measured (in situ) V S "> V S profiles, or National Earthquake Hazards Reduction Program (NEHRP) site class D to B, and are broadly comparable with estimates from previous data reports in the region. Site amplification estimates were calculated next from 1D SH transfer functions of the preferred V S "> V S VS profiles and then compared against observed horizontal‐to‐vertical spectral ratios (HVSRs) from nearby seismograph stations. We generally see good agreement between the predicted in situ model and the observed HVSR resonant frequencies, with nominal amplifications between 2 and 10 within the 2–15 Hz frequency band. Next, using 40 known in situ V S "> V S V S 30 "> 30 measurements in the region, we demonstrate that the in situ V S "> V S V S 30 "> 30 values improve the fit for selected suites of ground‐motion models (GMMs) for M 4+ earthquakes within a 50 km hypocentral distance when compared with proxy methods, arguing for future development of GMMs implementing in situ V S "> V S profiles.

Oklahoma↗

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↗

Estimating the effects of fish quality and size on the economic value of fishing in Oklahoma streams and rivers: A revealed preference and contingent behavior approach

Fishing in Oklahoma’s rivers and streams provides a unique experience for anglers in the state. Despite its popularity, information on total demand and economic benefits associated with stream fishing is limited in the state. Research on the role of site quality indicators, such as fish size and quantity, on recreational fishing has shown mixed results. Whether fish size or quantity plays an important role in determining fishing demand and economic value may have important management implications. We estimated the demand and economic value of fishing under varying scenarios by using anglers’ responses to hypothetical behavioral questions related to fishing in Ozark Highland streams and rivers in Oklahoma. We asked how intended number of trips might change in the future given hypothetical increases in catch rates of fish, catch rates of trophy-sized fish, and catch rates of preferred fish species, in combination with anglers’ trip-related data. Under current conditions, we estimated consumer surplus per person per trip to be $55 and aggregate value across all stream anglers in Oklahoma to be $68.51 million. Changes in marginal benefits varied among hypothetical scenarios of fish size and abundance but was maximized with a 25% increase in catch rates of trophy-sized fish. The study findings contribute to the understanding of the economic benefit of fishing in streams and suggest that fish size, rather than fish quantity, is more important to stream anglers in the area.

Oklahoma↗

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↗

Variation in black bass angler characteristics by stream size and accessibility in Oklahoma’s Ozark Highland streams

Fishing in streams and rivers is a popular outdoor recreation activity in eastern Oklahoma, where most anglers target black bass ( Micropterus ) species. Since the early 1990s, when the last assessment of black bass fishing in the region was conducted, broadscale factors such as harvesting behavior, state fishery regulations, and bass population dynamics have changed. In 2018, we conducted creel and fish tagging surveys in three tributaries of Lake Tenkiller (Caney Creek, Baron Fork, and Illinois River) that differed in size and accessibility to provide current estimates of catch, harvest, and effort directed toward black bass. We then related these estimates to angler socioeconomic characteristics. The amount of angler effort was concomitant with stream size and accessibility, being greatest in the largest stream with the most access (Illinois River). However, catch rates were highest in the medium-sized stream (Baron Fork). Harvest rates and exploitation were near zero in all systems. Anglers fishing Caney Creek, the smallest and least accessible stream, were nearly all local, coming from zip codes ~42 km away, with low median household incomes compared to anglers at the other streams who came from a broader array of more distant zip codes and had higher median household incomes. Anglers fishing the smallest stream were also more interested in harvesting fish and having higher creel limits than anglers at the other two systems. In the Oklahoma Ozark Highlands, stream size and accessibility appear to be a significant factor in angler demographics, potentially necessitating different management strategies.

Oklahoma↗

Bathymetric map, surface area, and capacity of Grand Lake O’ the Cherokees, northeastern Oklahoma, 2019

The U.S. Geological Survey (USGS), in cooperation with the Grand River Dam Authority, completed a high-resolution multibeam bathymetric survey to compute a new area and capacity table for Grand Lake O’ the Cherokees in northeastern Oklahoma. Area and capacity tables identify the relation between the elevation of the water surface and the volume of water that can be impounded at each water-surface elevation. The area and capacity of Grand Lake O’ the Cherokees were computed from a triangular irregular network surface created in Global Mapper Version 21.0.1. The triangular irregular network surface was created from three datasets: (1) a multibeam mapping system bathymetric survey of Grand Lake O’ the Cherokees completed during April–July 2019, (2) a previous bathymetric survey of the Neosho, Spring, and Elk Rivers, and (3) a 2010 USGS lidar-derived digital elevation model. The digital elevation model data were used in areas with land-surface elevations greater than 744 feet above the North American Vertical Datum of 1988 where the multibeam sonar data could not be collected. The 2019 multibeam sonar data were the predominant data used to compute the new area and capacity table for Grand Lake O’ the Cherokees.

Oklahoma↗

Survival and movement patterns of Rainbow Trout stocked in a groundwater-influenced warmwater stream

Stocking Rainbow Trout Oncorhynchus mykiss to create additional angling opportunities is common; however, the spatial and temporal dynamics of such stocking practices are unclear in groundwater-influenced, warmwater streams. Our objectives were to determine Rainbow Trout dispersal from a stocking location on Spavinaw Creek, Oklahoma and to quantify apparent survival of two cohorts of Rainbow Trout stocked in November 2018 and February 2019. Rainbow Trout were PIT-tagged prior to autumn ( N = 495) and winter ( N = 605) stocking and located poststocking using both an active floating array and two passive fixed arrays. We actively tracked a 6-km extent every 2–3 weeks poststocking and found that PIT-tagged Rainbow Trout dispersal ranged from 0 to 4 km with ~90% of detected tagged fish remaining within 1 km of the stocking location. Directional movement by stocked Rainbow Trout was evident with upstream movements related to increases in daily water discharge and downstream movements related to decreases in daily discharge. Estimated apparent weekly survival of Rainbow Trout was lowest during the first 2 weeks poststocking (91% and 75% for autumn and spring cohorts, respectively). Apparent weekly survival rates for both autumn and spring stocked fish increased through winter and spring when maximum water temperatures were below 25℃, reaching survival rates of 99% by the end of spring. Rainbow Trout persisted in Spavinaw Creek throughout the summer; however, the weekly apparent survival rate declined for both cohorts (93%) when maximum daily water temperatures exceeded 25°C. Our results can be used to guide stocking decisions in similar warmwater streams and suggest that Rainbow Trout persist through spring and summer poststocking under the conditions encountered during our study.

Oklahoma↗

Detectability and abundance of snowy plovers at Salt Plains National Wildlife Refuge, Oklahoma

In the past two decades, Salt Plains National Wildlife Refuge has been increasingly recognized as important habitat for both breeding and migratory shorebirds. North American snowy plovers Charadrius nivosus in particular rely on the nearly 5,000-ha salt flat at Salt Plains National Wildlife Refuge, which thousands use as breeding and stopover habitat. Elsewhere on the Southern Great Plains, decadal declines up to 75% within snowy plover subpopulations have been documented and attributed to vegetation encroachment, increased rates of nest predation, and decreased availability of fresh surface water. Despite many attempts to estimate this species' abundance across the continent, to date, no known attempt at distance sampling of snowy plovers has occurred. To address this paucity of data, we assessed feasibility of distance sampling methods to accurately estimate snowy plover abundance and detectability. Distance sampling surveys (2017–2018) indicated high detection probability ( P = 0.80) and the population abundance estimate across the salt flat extrapolated to 3,307 individuals. The distance-sampling population abundance estimate is lower than population abundance estimates determined by two previous studies within the past decade but far greater than 2,105 estimated for a study in 2006. Overall, distance sampling snowy plovers at Salt Plains National Wildlife Refuge proved to be an effective addition to pre-established survey protocols but further investigation is needed to compare accuracy and precision of methods used in this study, annual surveys conducted by Salt Plains National Wildlife Refuge, and other potential snowy plover surveys.

Oklahoma↗

Spectral inversion for seismic site response in central Oklahoma: Low-frequency resonances from the Great Unconformity

We investigate seismic site response by inverting seismic ground‐motion spectra for site and source spectral properties, in a region of central Oklahoma, where previous ground‐motion studies have indicated discrepancies between observations and ground‐motion models (GMMs). The inversion is constrained by a source spectral model, which we computed from regional seismic records, using aftershocks as empirical Green’s functions to deconvolve site and path effects. Site spectra across the region exhibit multiple, strong, low‐frequency ( ⁠ f <2 Hz ⁠ ) resonances. Modeling of vertically propagating SH waves reproduces the mean amplitudes and frequencies of the site spectra and requires a deep ( ⁠ ∼1–2 km ⁠ ) impedance contrast. Comparison of regional seismic velocity models and geologic profiles indicates that the seismic impedance contrast is, or is in proximity to, the Great Unconformity, which marks the interface between Precambrian basement rocks and overlying Paleozoic sedimentary rocks. Depth to Precambrian basement increases to the southwest across the study region ( ⁠ ∼1500–4500 m ⁠ ), and the fundamental frequencies of the site spectra are anticorrelated with basement depth. The first higher‐mode resonance also exhibits dependence on basement depth; although modeling suggests that the second higher mode should depend on basement depth, site spectra do not support this. The low‐frequency resonances in central Oklahoma are not represented in the GMMs used in current seismic hazard analyses for tectonic earthquakes, though approaches to account for such features are under consideration in other regions of the central and eastern United States. Given the broad spatial extent of the Great Unconformity underlying eastern North America, it is likely that similar effects on seismic site response also occur in other areas. This study highlights the impact of regional geologic structure on earthquake ground motions and reiterates the need for modeling regional effects to improve ground‐motion predictions and seismic hazard assessments.

Oklahoma↗

Using a large-n seismic array to explore the robustness of spectral estimations

Spectral analysis is widely used to estimate and refine earthquake source parameters such as source radius, seismic moment, and stress drop. This study aims to quantify the precision of the single spectra and empirical Green's function spectral ratio approach using the Large‐ n Seismic Survey in Oklahoma (LASSO) array. The dense station coverage in an area of local saltwater disposal offers a unique opportunity to observe and quantify radiation pattern effects and subsequent precision of spectral estimates of small earthquakes ( M < 3). The results suggest that the precision of source properties estimated from direct phase arrivals for arrays with less than 20 stations should be assumed to be not less than 30% and could be as high as 150% if less than five stations are used. Furthermore, we do not see clear evidence for, or against, a scaling of stress drop with magnitude of small earthquakes ( M < 3) as observed by other studies.

Oklahoma↗