Geology topics

Research about Kansas, Oklahoma

Source-linked reports with geographic coverage including Kansas, Oklahoma.

At least 19 recordsLinked to original sources

USGS2026

Bathymetric maps, surface areas, and storage capacities of Council Grove Lake and Marion Reservoir, Kansas, and Pine Creek Lake, Oklahoma, 2024

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, completed high-resolution multibeam bathymetric surveys to compute new elevation-area-capacity tables for Council Grove Lake and Marion Reservoir, Kansas, and Pine Creek Lake, Oklahoma. Elevation-area-capacity tables identify the relation between the water-surface elevation, surface area, and storage capacity of the lake. The surface areas and storage capacities of each lake were computed from bathymetric surfaces combining multibeam echo sounder data collected in 2024 and light detection and ranging point-cloud data collected in 2016 and 2018.

Kansas, Oklahoma↗

USGS2025

Induced earthquakes are generally not tidally triggered in Oklahoma and Kansas

Human-induced earthquakes occur along critically stressed faults as injected wastewater simultaneously heightens fluid pressure and pushes faults to failure. We investigate the possibility that small stresses imposed by Earth tides could trigger earthquakes in the induced seismicity region of Oklahoma and Kansas from 2011 to 2018. We decluster a catalog consisting of ∼110,000 earthquakes using three methods (Reasenberg, nearest-neighbor distance, and phase-bin). We find no significant tidal earthquake triggering using Schuster's p -value test for the declustered catalogs as a whole. We search for localized triggering using discretized space-time cells and find ∼0–6% of cells have significant tidal triggering which is close to what is randomly expected (5%) and indicates there is an insignificant amount of tidal triggering for the full study region. One area that has significant p -values across multiple time windows, ∼2014–2016 is ∼15 km from a region of large wastewater injection volume. It is possible that localized tidal triggering occurs for this time and area because faults remain critically stressed and are particularly susceptible to slip under the small stress load from the semidiurnal tide. Possible explanations for the lack of tidal triggering in our broader study are that the pre-seismic stressing rate in the earthquake nucleation area is faster than the tidal stressing rate (∼3 kPa/day), faults are not close enough to critically stressed to be affected by tidal forcing, and that nucleation occurs over longer periods than the tides considered in this study (∼1, ∼14 days). Fluid injection could be the source of a higher pre-seismic stress rate.

Kansas, Oklahoma↗

USGS2024

Recovery of working grasslands following a megafire in the southern mixed-grass prairie

While fire is a necessary ecological driver for grassland systems, Great Plains grasslands have undergone extensive land use change following European settlement (conversion, fragmentation, fire suppression, intensive grazing, etc.). Recent studies have documented the benefits of re-introducing fire to grasslands, but work has largely focused on small-scale, low-intensity fire, often at a pasture scale (i.e., prescribed fire). Over the last 30–40 years, the size and frequency of wildfires in the Great Plains has increased due to long-term fire suppression, woody encroachment, and climate change, resulting in megafires (wildfire >40,000 ha). While there is a wealth of information regarding effects of prescribed fire on Great Plains grasslands, knowledge of how large megafire events affect modern, fragmented working grasslands (i.e., grazed grasslands) is lacking and needed in the face of increasing megafire activity. To assess grassland response and recovery following a 2017 megafire (~254,000 ha), we compared vegetation characteristics pre- (2014–2015) and post-fire (2018–2019) in the mixed-grass prairie of Kansas, USA. We examined linkages between vegetation characteristics and a metric of burn severity (differenced normalized burn ratio [dNBR]) and evaluated megafire effects on limiting reproductive habitat for a declining grassland species, the lesser prairie-chicken ( Tympanuchus pallidicinctus ). One-year post-fire, we documented increased bare ground (+59%) and decreased visual obstruction (−39%), litter depth (−31%), and forb cover (−35%). Decreased visual obstruction and increased bare ground led to an 81% decrease of lesser prairie-chicken nest habitat in the first-year post-fire, but grassland structure, functional group cover, and available nest habitat largely recovered 2.5 years post-fire. Recovery to pre-fire conditions 2.5 years post-fire was primarily due to elevated growing season precipitation received in the years following the fire (>700 mm/year). Percent cover of grass (β = 0.38) and bare ground (β = −0.36) exhibited the strongest relationships with burn severity pre-fire, but overall grassland structure and functional group cover were not strongly influenced dNBR burn severity post-fire. While our results suggested that recovering grasslands were more homogenous due to the large size of megafire, working grasslands in the mixed-grass prairie appeared largely resilient to the effects of megafire.

Kansas, Oklahoma↗

USGS2022

Demographic effects of a megafire on a declining prairie grouse in the mixed-grass prairie

Recent studies have documented benefits of small, prescribed fire and wildfire for grassland-dependent wildlife, such as lesser prairie-chickens ( Tympanuchus pallidicintus ), but wildlife demographic response to the scale and intensity of megafire (wildfire >40,000 ha) in modern, fragmented grasslands remains unknown. Limited available grassland habitat makes it imperative to understand if increasing frequency of megafires could further reduce already declining lesser prairie-chicken populations, or if historical evolutionary interactions with fire make lesser prairie-chickens resilient. To evaluate lesser prairie-chicken demographic response to megafires, we compared lek counts, nest density, and survival rates of adults, nests, and chicks before (2014–2016) and after (2018–2020) a 2017 megafire in the mixed-grass prairie of Kansas, USA (Starbuck fire ~254,000 ha). There was a 67% decline in attending males on leks post-fire and a 57% decline in occupied leks post-fire. Despite population declines as indicated by lek counts, adult female breeding season survival ( Ŝ ) was similar pre- ( Ŝ = 0.65 ± 0.08 [SE]) and post-fire (0.61 ± 0.08), as was chick survival (pre-fire: 0.23 ± 0.07; post-fire: 0.27 ± 0.11). Nest survival appeared lower post-fire (pre-fire: 0.38 ± 0.06; post-fire: 0.20 ± 0.06), but did not differ at the 95% confidence interval. Nest density of marked females declined 73% in areas burned by megafire. Although lesser prairie-chickens persisted in the study area and we documented minimal effects on most demographic rates, reduced lesser prairie-chicken abundance and reproductive output suggests full recovery may take >3 years. Increased propensity for megafire resulting from suppression of smaller fires, compounded by climate change and woody encroachment, may impose a short-term (3–5 year) threat to already declining lesser prairie-chicken populations.

Kansas, Oklahoma↗

USGS2021

A big problem for small earthquakes: Benchmarking routine magnitudes and conversion relationships with coda-envelope-derived Mw in southern Kansas and northern Oklahoma

Earthquake magnitudes are widely relied upon measures of earthquake size. Although moment magnitude ( ⁠ M w "> M w Mw ⁠ ) has become the established standard for moderate and large earthquakes, difficulty in reliably measuring seismic moments for small (generally M w < 4 "> M w < 4 Mw<4 ⁠ ) earthquakes has meant that magnitudes for these events remain plagued by a patchwork of inconsistent measurement scales. Because of this, magnitudes of small earthquakes and statistics derived from them can be biased. Furthermore, because small earthquakes are much more numerous than large ones, many applications, such as seismic hazard modeling, depend critically on analysis of events characterized by magnitudes other than M w "> M w Mw ⁠ . To assess this problem, we apply coda envelope analysis to reliably determine moment magnitudes for a case study of small earthquakes from northern Oklahoma and southern Kansas. Not surprisingly, we find significant differences among M L "> M L ML ⁠ , m b L g "> m b L g mbLg ⁠ , and M w "> M w Mw for M ∼2–4 earthquakes examined here. More troublingly, we find that relations designed to convert other magnitudes to M w "> M w Mw ⁠ , which are relied upon for important applications such as seismic hazard analysis, often increase rather than decrease this bias for our dataset. In our case study, we find that converted magnitudes can result in a systematic bias sometimes exceeding 0.5 magnitude units, a difference that typically corresponds to a factor of ∼3 in seismicity rate. Moreover, we find a correspondingly large bias in Gutenberg–Richter b ‐values, controlled primarily by inaccurate magnitude scaling in the conversion relationships. Although this study focuses on a relatively small geographic area, we can expect that similar issues exist with varying severity in other regions. Therefore, magnitudes of small earthquakes and their associated statistics, including seismicity rates and b ‐values, should be treated with caution.

Kansas, Oklahoma↗

USGS2021

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

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

Kansas, Oklahoma↗

USGS2021

Fluid-earthquake and earthquake-earthquake interactions in southern Kansas, USA

An increase in injection activity associated with energy production in southern Kansas starting in 2013 has been linked to the occurrence of more than 130,000 earthquakes (M −1.5 to 4.9) between 2014 and 2017. Studies suggest that the dramatic increase in seismicity rate is related to wastewater injection into the highly permeable Arbuckle formation. Most of the seismicity is located in the underlying crystalline basement, for which hydrological properties and specific fault geometries are unknown. Additionally, some earthquake clusters occurred relatively far (tens of kilometers) from the main injection wells. Therefore, the effect of pore pressure diffusion may be insufficient to explain the relation between the volume of injected fluids and the spatiotemporal evolution of seismicity. Combining physical models (static stress and poroelasticity) and a statistical cluster analysis applied to a high-resolution relocated catalog, we analyze the evolution of seismicity in southern Kansas. We find that pore pressure changes (Δ p ) and Coulomb stress changes (ΔCFS) due to fluid diffusion smaller than 0.1 MPa are enough to initiate seismic sequences, which then evolve depending on their distance from the major injection wells. However, we find that earthquake sequences have different seismogenic responses to Δ p and ΔCFS in terms of triggering threshold. In regions located close to disposal wells (Harper area) our cluster analysis suggests that both earthquake interactions and fluid diffusion control the evolution of seismicity. On the other hand, at greater distances (Milan area), where clustering behavior suggests greater earthquake interactions, we find that coseismic ΔCFS are larger than Δ p .

Kansas, Oklahoma↗

USGS2021

Forecasting induced earthquake hazard using a hydromechanical earthquake nucleation model

In response to the dramatic increase in earthquake rates in the central United States, the U.S Geological Survey began releasing 1 yr earthquake hazard models for induced earthquakes in 2016. Although these models have been shown to accurately forecast earthquake hazard, they rely purely on earthquake statistics because there was no precedent for forecasting induced earthquakes based upon wastewater injection data. Since the publication of these hazard models, multiple physics‐based methods have been proposed to forecast earthquake rates using injection data. Here, we use one of these methods to generate earthquake hazard forecasts. Our earthquake hazard forecasts are more accurate than statistics‐based hazard forecasts. These results imply that fluid injection data, where and when available, and the physical implications of fluid injection should be included in future induced earthquake hazard forecasts.

Kansas, Oklahoma↗

USGS2019

The LArge-n Seismic Survey in Oklahoma (LASSO) experiment

In 2016, the U.S. Geological Survey deployed more than 1,800 vertical-component nodal seismometers in Grant County, Oklahoma to study induced seismic activity associated with production of the Mississippi Limestone Play. The LArge-n Seismic Survey in Oklahoma (LASSO) array operated for approximately one month, covering a 25-km-by-32-km region with a nominal station spacing of ~400 m. Primary goals of the deployment were to detect microseismic events not captured by the sparser regional network stations and to provide nearly unaliased records of the seismic wavefield. A more complete record of earthquakes allows us to map the spatiotemporal evolution of induced event sequences and illuminates the structures on which the events occur. Dense records of the seismic wavefield also provide improved measurements of the earthquake source, including focal mechanisms and stress drops. Taken together, we can use these findings to glean insights into the processes that induce earthquakes. Here, we describe the array layout, features of the nodal sensors, data recording configurations, and the field deployment. We also provide examples of earthquake waveforms recorded by the array to illustrate data quality and initial observations. LASSO array data provide a significant resource for understanding the occurrence of earthquakes induced by wastewater disposal.

Kansas, Oklahoma↗

USGS2019

Combining dynamic rupture simulations with ground motion data to characterize seismic hazard from Mw 3-5.8 earthquakes in Oklahoma and Kansas

Many seismically active areas suffer from a lack of near‐source ground‐motion recordings, making ground‐motion prediction difficult at distances within ∼ 40    km "> ∼ 40 km ∼40 km from an earthquake. We aim to aid the development of near‐source ground‐motion prediction equations (GMPEs) by generating synthetic ground‐motion data via simulation. Building on previous work using point‐source moment tensor sources to simulate small ( ⁠ M w "> M w Mw 3–4) earthquakes for a target region encompassing north central Oklahoma and south central Kansas, we perform dynamic rupture simulations of earthquakes up to M w "> M w Mw 5.8. We introduce complexity into the rupture process by adding stochastically generated heterogeneity to initial stress conditions. Our simulated ground‐motion data are added to a catalog of recorded ground‐motion data to construct a composite recorded‐simulated ground‐motion catalog that we use to develop a GMPE for the target region. This procedure can be generalized and used for GMPE development in other regions with near‐source ground‐motion data scarcity, which could directly benefit critical applications such as the National Seismic Hazard Maps produced by the U.S. Geological Survey (USGS).

Kansas, Oklahoma↗

USGS2019

Amphibian conservation in the Anthropocene

Research is necessary to identify patterns in nature, to understand how a system functions, and to make predictions about the future state of an ecosystem. Applied research in conservation biology can identify effective strategies to maintain biodiversity, though many papers end with the conclusion that more research is needed. However, more research does not necessarily lead to solutions. We use the ongoing global decline of amphibians as a salient example to highlight limitations in current conservation research, and to focus on finding solutions which are directly relevant for conservation. While research has been conducted since declines were first detected in the 1990s, outside a few specific examples, little progress in conservation has been achieved. We suggest that the case of amphibian declines is relevant to conservation science in general, as the current paradigm for conservation is that management is planned after research is completed; research and management are not effectively (and not directly) connected. This disconnect illustrates the knowledge-action divide which has been identified as a serious deficiency in conservation. Accordingly, we use this introductory paper to the Special Issue (Amphibian conservation in the Anthropocene: Progress and challenges) to describe amphibians as a conservation dilemma, and to make the case for a different, more pragmatic, and more solutions-focused view of conservation research.

Kansas, Oklahoma↗

USGS2019

Grassland bird and butterfly responses to Sericea lespedeza control via late-season grazing pressure

Sericea lespedeza ( Lespedeza cuneata ) is a high-tannin, late-season invasive forb species that reduces biodiversity in tallgrass prairie ecosystems. The largest tallgrass prairie remnant exists in the Flint Hills of Kansas and Oklahoma, where the most common grazing management practice involves prescribed fire in early spring followed by intensive stocking with yearling beef cattle from April to July. Sericea has continued to spread under this management regime. From 2013 to 2016, in Kansas Flint Hills tallgrass prairie, we tested the effects of using spring burning with early-season steer grazing, followed by late-season sheep grazing (Steer+Sheep) compared to spring burning followed by steer grazing only (Steer) on sericea vigor, grassland birds, and pollinators. Density and nest success of Grasshopper Sparrows (Ammodramus savannarum ) and Eastern Meadowlarks ( Sturnella magna ) were not negatively affected by Steer+Sheep relative to Steer treatments, whereas there was evidence of a negative effect in these same metrics for Dickcissels ( Spiza americana ). Abundance of butterflies and their nectar sources were similar between treatments but abundance of grassland specialist butterfly species was low, overall. Comprehensively, Steer+Sheep effectively controls the spread of sericea but may not create habitat for all tallgrass prairie wildlife species.

Kansas, Oklahoma↗

USGS2018

Ground motions from induced earthquakes in Oklahoma and Kansas

Improved predictions of earthquake ground motions are critical to advancing seismic hazard analyses and earthquake response. The high seismicity rate from 2009 to 2016 in Oklahoma and Kansas provides an extensive data set for examining the ground motions from these events. We evaluate the ability of three suites of ground‐motion prediction equations (GMPEs)—appropriate for modeling tectonic earthquakes in active crustal and stable continental regions—to reproduce the observed ground motions. Mixed‐effects regressions are used to separate the ground‐motion residuals into bias, between‐event, and within‐event terms. Although the residuals depict differing accuracies in the ability of the three GMPE suites to predict the ground motions, some consistent trends emerge in the period, magnitude, and distance dependence. The trends suggest that aspects of the ground motions from these induced earthquakes are not well modeled by current tectonic GMPEs. Most important, we find evidence for relatively poor overall fit to the ground motions, by all of the GMPE suites, at periods less than about 0.2 s and above 3 s, greater‐than‐predicted magnitude scaling for small to moderate‐magnitude events ( ⁠ M ≲ 5 "> M ≲ 5 ⁠ ), higher‐than‐predicted within‐event variability above 3 s, and an apparent geometric attenuation that is stronger than average predictions at close distances ( ⁠ R ≲ 20    km "> R ≲ 20 km ) and short periods ( ⁠ T ≲ 1 s "> T ≲ 1 s ). Our results indicate that regionally appropriate GMPEs for induced earthquakes should be used for predicting ground motions in Oklahoma and Kansas.

Kansas, Oklahoma↗

USGS2018

A flatfile of ground motion intensity measurements from induced earthquakes in Oklahoma and Kansas

We have produced a uniformly processed database of orientation-independent (RotD50, RotD100) ground motion intensity measurements containing peak horizontal ground motions (accelerations and velocities) and 5-percent-damped pseudospectral accelerations (0.1–10 s) from more than 3,800 M ≥ 3 earthquakes in Oklahoma and Kansas that occurred between January 2009 and December 2016. Ground motion time series were collected from regional, national, and temporary seismic arrays out to 500 km. We relocated the majority of the earthquake hypocenters using a multiple-event relocation algorithm to produce a set of near-uniformly processed hypocentral locations. Ground motion processing followed standard methods, with the primary objective of reducing the effects of noise on the measurements. Regional wave-propagation features and the high seismicity rate required careful selection of signal windows to ensure that we captured the entire ground motion record and that contaminating signals from extraneous earthquakes did not contribute to the database. Processing was carried out with an automated scheme and resulted in a database comprising more than 174,000 records ( https://dx.doi.org/10.5066/F73B5X8N ). We anticipate that these results will be useful for improved understanding of earthquake ground motions and for seismic hazard applications.

Kansas, Oklahoma↗

USGS2018

Damage accumulation for a two-story wood-frame building in sequences of induced earthquakes

In this study, a nonlinear model of a two-story multifamily wood-frame residential structure is subjected to recordings of sequences of induced earthquakes in order to quantify changes in fragility and accumulation of damage throughout multiple earthquake loadings. Initial efforts consisting of ground motion selection, building design, numerical modeling, and preliminary results are presented. Damage is quantified through a seismic loss estimation procedure that accounts for damage to nonstructural and structural components of the building. Examining damage accumulation from sequential earthquake shaking enables us to explore how the occurrence of damage in an earthquake, even to relatively small levels, may increase a structure’s susceptibility to collapse or damage in subsequent ground shaking. Ultimately, the goal is to compare damage fragilities and seismic losses as a function of the building’s initial damage state.

Kansas, Oklahoma↗

USGS2017

Considerations in comparing the U.S. Geological Survey one‐year induced‐seismicity hazard models with “Did You Feel It?” and instrumental data

The recent steep increase in seismicity rates in Oklahoma, southern Kansas, and other parts of the central United States led the U.S. Geological Survey (USGS) to develop, for the first time, a probabilistic seismic hazard forecast for one year (2016) that incorporates induced seismicity. In this study, we explore a process to ground‐truth the hazard model by comparing it with two databases of observations: modified Mercalli intensity (MMI) data from the “Did You Feel It?” (DYFI) system and peak ground acceleration (PGA) values from instrumental data. Because the 2016 hazard model was heavily based on earthquake catalogs from 2014 to 2015, this initial comparison utilized observations from these years. Annualized exceedance rates were calculated with the DYFI and instrumental data for direct comparison with the model. These comparisons required assessment of the options for converting hazard model results and instrumental data from PGA to MMI for comparison with the DYFI data. In addition, to account for known differences that affect the comparisons, the instrumental PGA and DYFI data were declustered, and the hazard model was adjusted for local site conditions. With these adjustments, examples at sites with the most data show reasonable agreement in the exceedance rates. However, the comparisons were complicated by the spatial and temporal completeness of the instrumental and DYFI observations. Furthermore, most of the DYFI responses are in the MMI II–IV range, whereas the hazard model is oriented toward forecasts at higher ground‐motion intensities, usually above about MMI IV. Nevertheless, the study demonstrates some of the issues that arise in making these comparisons, thereby informing future efforts to ground‐truth and improve hazard modeling for induced‐seismicity applications.

Kansas, Oklahoma↗

USGS2016

Likelihood testing of seismicity-based rate forecasts of induced earthquakes in Oklahoma and Kansas

Likelihood testing of induced earthquakes in Oklahoma and Kansas has identified the parameters that optimize the forecasting ability of smoothed seismicity models and quantified the recent temporal stability of the spatial seismicity patterns. Use of the most recent 1-year period of earthquake data and use of 10–20-km smoothing distances produced the greatest likelihood. The likelihood that the locations of January–June 2015 earthquakes were consistent with optimized forecasts decayed with increasing elapsed time between the catalogs used for model development and testing. Likelihood tests with two additional sets of earthquakes from 2014 exhibit a strong sensitivity of the rate of decay to the smoothing distance. Marked reductions in likelihood are caused by the nonstationarity of the induced earthquake locations. Our results indicate a multiple-fold benefit from smoothed seismicity models in developing short-term earthquake rate forecasts for induced earthquakes in Oklahoma and Kansas, relative to the use of seismic source zones.

Kansas, Oklahoma↗

USGS2015

Climate, water use, and land surface transformation in an irrigation intensive watershed - streamflow responses from 1950 through 2010

Climatic variability and land surface change have a wide range of effects on streamflow and are often difficult to separate. We analyzed long-term records of climate, land use and land cover, and re-constructed the water budget based on precipitation, groundwater levels, and water use from 1950 through 2010 in the Cimarron–Skeleton watershed and a portion of the Cimarron–Eagle Chief watershed in Oklahoma, an irrigation-intensive agricultural watershed in the Southern Great Plains, USA. Our results show that intensive irrigation through alluvial aquifer withdrawal modifies climatic feedback and alters streamflow response to precipitation. Increase in consumptive water use was associated with decreases in annual streamflow, while returning croplands to non-irrigated grasslands was associated with increases in streamflow. Along with groundwater withdrawal, anthropogenic-induced factors and activities contributed nearly half to the observed variability of annual streamflow. Streamflow was more responsive to precipitation during the period of intensive irrigation between 1965 and 1984 than the period of relatively lower water use between 1985 and 2010. The Cimarron River is transitioning from a historically flashy river to one that is more stable with a lower frequency of both high and low flow pulses, a higher baseflow, and an increased median flow due in part to the return of cropland to grassland. These results demonstrated the interrelationship among climate, land use, groundwater withdrawal and streamflow regime and the potential to design agricultural production systems and adjust irrigation to mitigate impact of increasing climate variability on streamflow in irrigation intensive agricultural watershed.

Kansas, Oklahoma↗