Search USGS⌕ Search

Geology topics

Research about Oklahoma

Source-linked reports with geographic coverage including Oklahoma.

At least 109 records · Page 6Linked to original sources
USGS2018

Precision of four otolith techniques for estimating age of white perch from a thermally altered reservoir

The White Perch Morone americana is an invasive species in many Midwestern states and is widely distributed in reservoir systems, yet little is known about the species' age structure and population dynamics. White Perch were first observed in Sooner Reservoir, a thermally altered cooling reservoir in Oklahoma, by the Oklahoma Department of Wildlife Conservation in 2006. It is unknown how thermally altered systems like Sooner Reservoir may affect the precision of White Perch age estimates. Previous studies have found that age structures from Largemouth Bass Micropterus salmoides and Bluegills Lepomis macrochirus from thermally altered reservoirs had false annuli, which increased error when estimating ages. Our objective was to quantify the precision of White Perch age estimates using four sagittal otolith preparation techniques (whole, broken, browned, and stained). Because Sooner Reservoir is thermally altered, we also wanted to identify the best month to collect a White Perch age sample based on aging precision. Ages of 569 White Perch (20–308 mm TL) were estimated using the four techniques. Age estimates from broken, stained, and browned otoliths ranged from 0 to 8 years; whole‐view otolith age estimates ranged from 0 to 7 years. The lowest mean coefficient of variation (CV) was obtained using broken otoliths, whereas the highest CV was observed using browned otoliths. July was the most precise month (lowest mean CV) for estimating age of White Perch, whereas April was the least precise month (highest mean CV). These results underscore the importance of knowing the best method to prepare otoliths for achieving the most precise age estimates and the best time of year to obtain those samples, as these factors may affect other estimates of population dynamics.

Oklahoma↗
USGS2018

Quantifying seining detection probability for fishes of Great Plains sand‐bed rivers

Species detection error (i.e., imperfect and variable detection probability) is an essential consideration when investigators map distributions and interpret habitat associations. When fish detection error that is due to highly variable instream environments needs to be addressed, sand‐bed streams of the Great Plains represent a unique challenge. We quantified seining detection probability for diminutive Great Plains fishes across a range of sampling conditions in two sand‐bed rivers in Oklahoma. Imperfect detection resulted in underestimates of species occurrence using naïve estimates, particularly for less common fishes. Seining detection probability also varied among fishes and across sampling conditions. We observed a quadratic relationship between water depth and detection probability, in which the exact nature of the relationship was species‐specific and dependent on water clarity. Similarly, the direction of the relationship between water clarity and detection probability was species‐specific and dependent on differences in water depth. The relationship between water temperature and detection probability was also species dependent, where both the magnitude and direction of the relationship varied among fishes. We showed how ignoring detection error confounded an underlying relationship between species occurrence and water depth. Despite imperfect and heterogeneous detection, our results support that determining species absence can be accomplished with two to six spatially replicated seine hauls per 200‐m reach under average sampling conditions; however, required effort would be higher under certain conditions. Detection probability was low for the Arkansas River Shiner Notropis girardi , which is federally listed as threatened, and more than 10 seine hauls per 200‐m reach would be required to assess presence across sampling conditions. Our model allows scientists to estimate sampling effort to confidently assess species occurrence, which maximizes the use of available resources. Increased implementation of approaches that consider detection error promote ecological advancements and conservation and management decisions that are better informed.

Oklahoma↗
USGS2018

Pore-types and pore-network evolution in Upper Devonian-Lower Mississippian Woodford and Mississippian Barnett mudstones: Insights from laboratory thermal maturation and organic petrology

Pore-evolution models from immature organic-matter (OM) -rich Barnett (0.42%R o ) and Woodford (0.49%R o ) mudstones were compared with models previously developed from low-maturity OM-lean Boquillas (Eagle Ford-equivalent) mudstones to investigate whether (1) different mineralogy (siliceous vs. calcareous) exerts different catalytic and sorption effects and influences OM-pore origin and evolution; and (2) different types of macerals show different OM pore evolution history. Laboratory gold-tube pyrolysis, scanning electron microscopy (SEM) and thin-section petrography, organic petrography, and geochemical characterization were used to investigate the role of bulk mineralogy, maceral type, and thermal maturation on OM-pore evolution. Results suggest that mineralogy has little impact on OM-pore development and evolution. Macerals, identified using both SEM (platy OM, particulate OM, organic–mineral admixtures, Tasmanites ) and organic petrology (vitrinite, inertinite, amorphous organic matter [AOM]/bituminite, telalginite [ Leiosphaeridia , Tasmanites ]), do affect the origin and evolution of OM pores owing to differences in chemical compositions, generation kinetics, and activation-energy distributions between Tasmanites , matrix bituminite, and other types of macerals. Leiosphaeridia and Tasmanites in Woodford mudstone samples exhibit a delay in onset and a shorter period of petroleum generation and pore development compared to the matrix bituminite in the Barnett and Woodford mudstone samples. Pre-oil solid bitumen was observed to have migrated into initial primary mineral pore networks at the bitumen generation stage in both Barnett and Woodford samples. At higher levels of thermal maturation, the volume of primary mineral pores decreases and the pore volume composed of modified mineral pores and OM pores becomes greater. Pore evolution and pore-type heterogeneity in these mudstones is a function of the initial mineral pore network, types of kerogen and macerals, and generation kinetics of individual macerals upon thermal maturation.

Oklahoma↗
USGS2018

Groundwater flux estimation in streams: A thermal equilibrium approach

Stream and groundwater interactions play an essential role in regulating flow, temperature, and water quality for stream ecosystems. Temperature gradients have been used to quantify vertical water movement in the streambed since the 1960s, but advancements in thermal methods are still possible. Seepage runs are a method commonly used to quantify exchange rates through a series of streamflow measurements but can be labor and time intensive. The objective of this study was to develop and evaluate a thermal equilibrium method as a technique for quantifying groundwater flux using monitored stream water temperature at a single point and readily available hydrological and atmospheric data. Our primary assumption was that stream water temperature at the monitored point was at thermal equilibrium with the combination of all heat transfer processes, including mixing with groundwater. By expanding the monitored stream point into a hypothetical, horizontal one-dimensional thermal modeling domain, we were able to simulate the thermal equilibrium achieved with known atmospheric variables at the point and quantify unknown groundwater flux by calibrating the model to the resulting temperature signature. Stream water temperatures were monitored at single points at nine streams in the Ozark Highland ecoregion and five reaches of the Kiamichi River to estimate groundwater fluxes using the thermal equilibrium method. When validated by comparison with seepage runs performed at the same time and reach, estimates from the two methods agreed with each other with an R 2 of 0.94, a root mean squared error (RMSE) of 0.08 (m/d) and a Nash–Sutcliffe efficiency (NSE) of 0.93. In conclusion, the thermal equilibrium method was a suitable technique for quantifying groundwater flux with minimal cost and simple field installation given that suitable atmospheric and hydrological data were readily available.

Oklahoma↗
USGS2018

Induced earthquake and liquefaction hazards in Oklahoma, USA: Constraints from InSAR

Oklahoma experienced three earthquakes of M w 5.0 or greater in 2016: the 13-Feb. Fairview earthquake (M w 5.1), the 03-Sep. Pawnee earthquake (M w 5.8), and the 07-Nov. Cushing earthquake (M w 5.0). These events are the first earthquakes in the state exceeding M w 5.0 since the 2011 M w 5.7 Prague earthquake and likely result from wide-scale deep fluid-injection. We use interferometric synthetic aperture radar (InSAR) observations to quantify the magnitude and location of surface deformation associated with these three events, determine the depth ranges of fault slip, and assess the spatial relationship between fault slip and well-calibrated mainshock and aftershock locations. We also include newly reported, calibrated event locations for the Cushing earthquake. We find that the Pawnee earthquake ruptured within the crystalline basement with the shallowest slip occurring at depths of 3.1–4.3 km. We find a similar, though shallower, crystalline basement source for the Cushing earthquake with a minimum depth to slip of 1.6–2.3 km. Despite the smaller magnitude of the Cushing earthquake, it generated anomalously high ground motions and damage compared to the larger Pawnee and Fairview earthquakes. We postulate that the shallow source of the Cushing earthquakes provides one explanation for the higher than expected ground motions. The Fairview earthquake generated no detectable co-seismic displacements, which is consistent with a relatively deep earthquake source (~8.5 km). We do, however, identify a 16 km stretch of floodplain where widespread liquefaction occurred in response to the Fairview earthquake, and where 30 gas production wells were exposed to surface displacements exceeding 5 cm. Consequently, the depth to crystalline basement, which limits the depth of injection-induced earthquakes in Oklahoma, and the potential for liquefaction are important factors in assessing shaking risk in the central United States.

Oklahoma↗
USGS2018

Long-Term Movement and Estimated Age of a Paddlefish (Polyodon spathula) in the Arkansas River Basin of Oklahoma

We report the age and distance moved for an individual paddlefish ( Polyodon spathula ) that was tagged March 1998 in the Cimarron River Arm of Keystone Lake, Oklahoma, and snagged by an angler in April 2016 downstream of Eufaula Dam, Oklahoma. The fish was part of a cohort spawned in 1995. At the time of initial capture, the fish measured 795 mm eye–fork length, was estimated to be 3 y old, and 18 y had elapsed before its recapture by an angler in 2016, indicating this fish was 21 y old at recapture. Although paddlefish as old as 27 have been estimated in the Grand River basin of Oklahoma, this is the oldest fish known in the Arkansas River basin of Oklahoma. At the place of its recapture, this fish would have traveled approximately 235 km, passing downstream through three dams before moving upstream to Eufaula Dam.

Oklahoma↗
USGS2018

Assessing wild juvenile trout ecology in the lower Mountain Fork

Reservoir tailwaters can be valuable fisheries for Rainbow Trout ( Oncorhynchus mykiss ), which is commonly stocked as mitigation for the altered habitat because it performs well as a put-and-take species in these thermally depressed systems. These fisheries are usually sustained by stocking due to flow fluctuations and lack of suitable spawning habitat that may limit natural reproduction. The Lower Mountain Fork River (LMFR) below Broken Bow Dam in southeastern Oklahoma is one of two year-round trout fisheries in the state and wild, juvenile Rainbow Trout were documented beginning in 2006, prompting speculation about the potential for a self-sustaining population. To determine this potential, we searched several sites over two years throughout the LMFR for wild, juvenile Rainbow Trout to estimate several parameters related to their population status (e.g., age, growth, time of spawning and hatching, and prey use). We also assessed the availability of macroinvertebrate prey to determine how food resources may affect trout sustainability. We found wild, juvenile Rainbow Trout each year, but only at sites within the first 4.5 km of the 19-km tailwater. Juvenile trout were the result of spawning that took place from late January through mid-April. Growth and body condition were variable between years, but similar to other systems. Weekly survival estimates using catch curves were low (<80%), suggesting limited potential for recruitment; however, declining catchability of larger juvenile fish likely biased these estimates. Wild, juvenile Rainbow Trout ate a variety of food items, but selected for Amphipoda and Diplostraca and against Trichoptera. Overlap in diet with adult Rainbow Trout was low (Bray-Curtis dissimilarity = 0.70). Macroinvertebrate prey resources available to trout varied among management zones, being most abundant in Zone 1 and Zone 3. Potential for a wild fishery may exist in the upper portion of the LMFR but additional research on recruitment to adulthood would be required to provide a more definitive answer.

Oklahoma↗
USGS2017

Wastewater disposal and the earthquake sequences during 2016 near Fairview, Pawnee, and Cushing, Oklahoma

Each of the three earthquake sequences in Oklahoma in 2016—Fairview, Pawnee, and Cushing—appears to have been induced by high-volume wastewater disposal within 10 km. The Fairview M 5.1 main shock was part of a 2 year sequence of more than 150 events of M 3, or greater; the main shock accounted for about half of the total moment. The foreshocks and aftershocks of the M 5.8 Pawnee earthquake were too small and too few to contribute significantly to the cumulative moment; instead, nearly all of the moment induced by wastewater injection was focused on the main shock. The M 5.0 Cushing event is part of a sequence that includes 48 earthquakes of M 3, or greater, that are mostly foreshocks. The cumulative moment for each of the three sequences during 2016, as well as that for the 2011 Prague, Oklahoma, and nine other sequences representing a broad range of injected volume, are all limited by the total volumes of wastewater injected locally.

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

Low stress drops observed for aftershocks of the 2011 Mw 5.7 Prague, Oklahoma, earthquake

In November 2011, three M w ≥ 4.8 earthquakes and thousands of aftershocks occurred along the structurally complex Wilzetta fault system near Prague, Oklahoma. Previous studies suggest that wastewater injection induced a M w 4.8 foreshock, which subsequently triggered a M w 5.7 mainshock. We examine source properties of aftershocks with a standard Brune-type spectral model and jointly solve for seismic moment ( M 0 ), corner frequency ( f 0 ), and kappa ( κ ) with an iterative Gauss-Newton global downhill optimization method. We examine 934 earthquakes with initial moment magnitudes ( M w ) between 0.33 and 4.99 based on the pseudospectral acceleration and recover reasonable M 0 , f 0 , and κ for 87 earthquakes with M w 1.83–3.51 determined by spectral fit. We use M 0 and f 0 to estimate the Brune-type stress drop, assuming a circular fault and shear-wave velocity at the hypocentral depth of the event. Our observations suggest that stress drops range between 0.005 and 4.8 MPa with a median of 0.2 MPa (0.03–26.4 MPa with a median of 1.1 MPa for Madariaga-type), which is significantly lower than typical eastern United States intraplate events (>10 MPa). We find that stress drops correlate weakly with hypocentral depth and magnitude. Additionally, we find the stress drops increase with time after the mainshock, although temporal variation in stress drop is difficult to separate from spatial heterogeneity and changing event locations. The overall low median stress drop suggests that the fault segments may have been primed to fail as a result of high pore fluid pressures, likely related to nearby wastewater injection.

Oklahoma↗
USGS2017

Poroelastic properties of the Arbuckle Group in Oklahoma derived from well fluid level response to the 3 September 2016 M w 5.8 Pawnee and 7 November 2016 M w 5.0 Cushing earthquakes

The Arbuckle Group (Arbuckle) is a basal sedimentary unit that is the primary target for saltwater disposal in Oklahoma. Thus, the reservoir characteristics of the Arbuckle, including how the poroelastic properties change laterally and over time are of significant interest. We report observations of fluid level changes in two monitoring wells in response to the 3 September 2016 M w 5.8 Pawnee and the 7 November 2016 M w 5.0 Cushing earthquakes. We investigate the relationship between static strain resulting from these events and the fluid level changes observed in the wells. We model the fluid level response by estimating static strains from a set of earthquake source parameters and spatiotemporal poroelastic properties of the Arbuckle in the neighborhood of the monitoring wells. Results suggest that both the direction of the observed fluid level step and the amplitude can be predicted from the computed volumetric strain change and a reasonable set of poroelastic parameters. Modeling results indicate that poroelastic parameters differ at the time of the Pawnee and Cushing earthquakes, with a moderately higher Skempton’s coefficient required to fit the response to the Cushing earthquake. This may indicate that dynamic shaking resulted in physical alteration of the Arbuckle at distances up to ∼50 km from the Pawnee earthquake.

Oklahoma↗
USGS2017

Geodetic slip model of the 3 September 2016 Mw 5.8 Pawnee, Oklahoma, earthquake: Evidence for fault‐zone collapse

The 3 September 2016 M w 5.8 Pawnee earthquake in northern Oklahoma is the largest earthquake ever recorded in Oklahoma. The coseismic deformation was measured with both Interferometric Synthetic Aperture Radar and Global Positioning System (GPS), with measureable signals of order 1 cm and 1 mm, respectively. We derive a coseismic slip model from Sentinel‐1A and Radarsat 2 interferograms and GPS static offsets, dominated by distributed left‐lateral strike slip on a primary west‐northwest–east‐southeast‐trending subvertical plane, whereas strike slip is concentrated near the hypocenter (5.6 km depth), with maximum slip of ∼1 m located slightly east and down‐dip of the hypocenter. Based on systematic misfits of observed interferogram line‐of‐sight (LoS) displacements, with LoS based on shear‐dislocation models, a few decimeters of fault‐zone collapse are inferred in the hypocentral region where coseismic slip was the largest. This may represent the postseismic migration of large volumes of fluid away from the high‐slip areas, made possible by the creation of a temporary high‐permeability damage zone around the fault.

Oklahoma↗
USGS2017

Comparative precision of age estimates from two southern reservoir populations of paddlefish [Polyodon spathula (Walbaum, 1792)]

The aim of the study was to determine whether location and sex affected the age precision estimates between two southern, reservoir populations of paddlefish [ Polyodon spathula (Walbaum, 1792)]. From 589 paddlefish collected in Grand Lake and Keystone Lake, Oklahoma in 2011, ages from dentaries were estimated using three independent readers and precision was compared with coefficient of variation between locations and sexes. Ages were more precisely estimated from Grand Lake and from females.

Oklahoma↗
USGS2017

Bathymetry and capacity of Shawnee Reservoir, Oklahoma, 2016

Shawnee Reservoir (locally known as Shawnee Twin Lakes) is a man-made reservoir on South Deer Creek with a drainage area of 32.7 square miles in Pottawatomie County, Oklahoma. The reservoir consists of two lakes connected by an equilibrium channel. The southern lake (Shawnee City Lake Number 1) was impounded in 1935, and the northern lake (Shawnee City Lake Number 2) was impounded in 1960. Shawnee Reservoir serves as a municipal water supply, and water is transferred about 9 miles by gravity to a water treatment plant in Shawnee, Oklahoma. Secondary uses of the reservoir are for recreation, fish and wildlife habitat, and flood control. Shawnee Reservoir has a normal-pool elevation of 1,069.0 feet (ft) above North American Vertical Datum of 1988 (NAVD 88). The auxiliary spillway, which defines the flood-pool elevation, is at an elevation of 1,075.0 ft. The U.S. Geological Survey (USGS), in cooperation with the City of Shawnee, has operated a real-time stage (water-surface elevation) gage (USGS station 07241600) at Shawnee Reservoir since 2006. For the period of record ending in 2016, this gage recorded a maximum stage of 1,078.1 ft on May 24, 2015, and a minimum stage of 1,059.1 ft on April 10–11, 2007. This gage did not report reservoir storage prior to this report (2016) because a sufficiently detailed and thoroughly documented bathymetric (reservoir-bottom elevation) survey and corresponding stage-storage relation had not been published. A 2011 bathymetric survey with contours delineated at 5-foot intervals was published in Oklahoma Water Resources Board (2016), but that publication did not include a stage-storage relation table. The USGS, in cooperation with the City of Shawnee, performed a bathymetric survey of Shawnee Reservoir in 2016 and released the bathymetric-survey data in 2017. The purposes of the bathymetric survey were to (1) develop a detailed bathymetric map of the reservoir and (2) determine the relations between stage and reservoir storage capacity and between stage and reservoir surface area. The bathymetric map may serve as a baseline to which temporal changes in storage capacity, due to sedimentation and other factors, can be compared. The stage-storage relation may be used in the reporting of real-time Shawnee Reservoir storage capacity at USGS station 07241600 to support water-resource management decisions by the City of Shawnee.

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

Oklahoma experiences largest earthquake during ongoing regional wastewater injection hazard mitigation efforts

The 3 September 2016, M w 5.8 Pawnee earthquake was the largest recorded earthquake in the state of Oklahoma. Seismic and geodetic observations of the Pawnee sequence, including precise hypocenter locations and moment tensor modeling, shows that the Pawnee earthquake occurred on a previously unknown left-lateral strike-slip basement fault that intersects the mapped right-lateral Labette fault zone. The Pawnee earthquake is part of an unprecedented increase in the earthquake rate in Oklahoma that is largely considered the result of the deep injection of waste fluids from oil and gas production. If this is, indeed, the case for the M 5.8 Pawnee earthquake, then this would be the largest event to have been induced by fluid injection. Since 2015, Oklahoma has undergone wide-scale mitigation efforts primarily aimed at reducing injection volumes. Thus far in 2016, the rate of M 3 and greater earthquakes has decreased as compared to 2015, while the cumulative moment—or energy released from earthquakes—has increased. This highlights the difficulty in earthquake hazard mitigation efforts given the poorly understood long-term diffusive effects of wastewater injection and their connection to seismicity.

Oklahoma↗
USGS2017

Geologic influence on induced seismicity: Constraints from potential field data in Oklahoma

Recent Oklahoma seismicity shows a regional correlation with increased wastewater injection activity, but local variations suggest that some areas are more likely to exhibit induced seismicity than others. We combine geophysical and drill hole data to map subsurface geologic features in the crystalline basement, where most earthquakes are occurring, and examine probable contributing factors. We find that most earthquakes are located where the crystalline basement is likely composed of fractured intrusive or metamorphic rock. Areas with extrusive rock or thick (>4 km) sedimentary cover exhibit little seismicity, even in high injection rate areas, similar to deep sedimentary basins in Michigan and western North Dakota. These differences in seismicity may be due to variations in permeability structure: within intrusive rocks, fluids can become narrowly focused in fractures and faults, causing an increase in local pore fluid pressure, whereas more distributed pore space in sedimentary and extrusive rocks may relax pore fluid pressure.

Oklahoma↗