USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geophysical Journal”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,423 records · Page 79Linked to original sources

Dynamics of marsh-derived sediments in lagoon-type estuaries

Salt marshes are valuable ecosystems that must trap sediments and accrete in order to counteract the deleterious effect of sea‐level rise. Previous studies have shown that the capacity of marshes to build up vertically depends on both autogenous and exogenous processes including eco‐geomorphic feedbacks and sediment supply from in‐land and coastal ocean. There have been numerous efforts to quantify the role played by the sediments coming from marsh edge erosion on the resistance of salt marshes to sea‐level rise. However, the majority of existing studies investigating the interplay between lateral and vertical dynamics use simplified modelling approaches and they do not consider that marsh retreat can affect the regional scale hydrodynamics and sediment retention in back‐barrier basins. In this study, we evaluated the fate of the sediments originating from marsh lateral loss by using high‐resolution numerical model simulations of Jamaica Bay, a small lagoonal estuary located in New York City. Our findings show that up to 42% of the sediments released during marsh edge erosion deposits on the shallow areas of the basin and over the vegetated marsh platforms, contributing positively to the sediment budget of the remaining salt marshes. Furthermore, we demonstrate that with the present‐day sediment supply from the ocean the system cannot keep pace with sea‐level rise even accounting for the sediment liberated in the bay through marsh degradation. Our study highlights the relevance of multiple sediment sources for the maintenance of the marsh complex.

New York↗

Unsaturated flow processes and the onset of seasonal deformation in slow-moving landslides

Predicting rainfall-induced landslide motion is challenging because shallow groundwater flow is extremely sensitive to the preexisting moisture content in the ground. Here, we use groundwater hydrology theory and numerical modeling combined with five years of field monitoring to illustrate how unsaturated groundwater flow processes modulate the seasonal pore water pressure rise and therefore the onset of motion for slow-moving landslides. The onset of landslide motion at Oak Ridge earthflow in California’s Diablo Range occurs after an abrupt water table rise to near the landslide surface 52–129 days after seasonal rainfall commences. Model results and theory suggest that this abrupt rise occurs from the advection of a nearly saturated wetting front, which marks the leading edge of the integrated downward flux of seasonal rainfall, to the water table. Prior to this abrupt rise, we observe little measured pore water pressure response within the landslide due to rainfall. However, once the wetting front reaches the water table, we observe nearly instantaneous pore water pressure transmission within the landslide body that is accompanied by landslide acceleration. We cast the timescale to reach a critical pore water pressure threshold using a simple mass balance model that considers variable moisture storage with depth and explains the onset of seasonal landslide motion with a rainfall intensity-duration threshold. Our model shows that the seasonal response time of slow-moving landslides is controlled by the dry season vadose zone depth rather than the total landslide thickness.

California↗

Identifying geomorphic process domains in the synthetic landscapes of West Virginia, USA

Human activities such as mining, agriculture, and urbanization have resulted in severe, large‐scale alteration to landform organization and associated geomorphic processes. The mountaintop mining (MTM) region of West Virginia, USA has experienced dramatic topographic alteration, by removing steep slopes and introducing plateau‐like areas at ridgelines and benches on valley fills. The resulting engineered landforms create synthetic landscapes, disconnected from previous geomorphic processes. Invoking the process domain concept, we compare differences in slope‐area relations, cumulative area distributions (CADs), elevation, slope, upslope accumulated area, and a slope*area product before and after mining to adjacent unmined sub‐catchments in five study basins. Differences in the slope‐area relation include a 42% slope reduction in low drainage areas, corresponding to hillslopes, unchanneled valleys, and debris flow dominated channels, which may fall below thresholds required for debris flow processes. The curved slope‐area relation that represents valley incision by debris flows is replaced by slope‐area relations that resemble basins where gullying and the stream power law dominate. Extremely high chemical weathering of unconsolidated valley fills materials may facilitate process domain shifts from debris flows to gullying and fluvial erosion. The characteristic power law scaling break in CADs that represents the headward limit of the channelized network is subdued in post‐mined sites and may reflect headward channelized network extension in mined basins. Slope‐area relations and CADs present a unique topographic signature of MTM activity, potentially providing an analytical approach to assess impacts on underlying geomorphic processes for other synthetic landscapes such as cities or large‐scale agricultural production.

West Virginia↗

Extreme precipitation across adjacent burned and unburned watersheds reveals impacts of low severity wildfire on debris-flow processes

In steep landscapes, wildfire-induced changes to soil and vegetation can lead to extreme and hazardous geomorphic responses, including debris flows. The wildfire-induced mechanisms that lead to heightened geomorphic responses, however, depend on many site-specific factors including regional climate, vegetation, soil texture, and soil burn severity. As climate and land use change drive changes in fire regime, there is an increasing need to understand how fire alters geomorphic responses, particularly in areas where fire has been historically infrequent. Here, we examine differences in the initiation, magnitude, and particle-size distribution of debris flows that initiated within the area burned by the 2019 Woodbury Fire in central Arizona, USA, and those that initiated in a nearby unburned area. Despite similar rainfall intensities, unburned watersheds were less likely to produce debris flows. Debris flows in unburned areas initiated from both runoff and shallow landslides, while debris flows only initiated from runoff-related processes in the burned area. The grain-size distribution making up the matrix of debris-flow deposits within the burned area generally had a lower ratio of sand to silt relative to debris flows that initiated in the unburned area, though there were no systematic differences in the coarse fraction of debris-flow sediment between burned and unburned areas. Results help expand our ability to predict postwildfire debris-flow activity in a wider range of settings, specifically the Sonoran Desert ecoregion, and provide general insight into the impact of wildfire on geomorphic processes in steep terrain.

Arizona↗

Movement of sediment through a burned landscape: Sediment volume observations and model comparisons in the San Gabriel Mountains, California, USA

Post-wildfire changes to hydrologic and geomorphic systems can lead to widespread sediment redistribution. Understanding how sediment moves through a watershed is crucial for assessing hazards, developing debris flow inundation models, engineering sediment retention solutions, and quantifying the role that disturbances play in landscape evolution. In this study, we used terrestrial and airborne lidar to measure sediment redistribution in the 2016 Fish Fire, in the San Gabriel Mountains in southern California, USA. The lidar areas are in two adjacent watersheds, at spatial scales of 900 m 2 to 4 km 2 , respectively. Terrestrial lidar data were acquired prior to rainfall, and two subsequent surveys show erosional change after rainstorms. Two airborne lidar flights occurred (1) 7 months before, and (2) 14 months after the fire ignition, capturing the erosional effects after rainfall. We found hillslope erosion dominated the overall sediment budget in the first rainy season after wildfire. Only 7% of the total erosion came from the active channel bed and channel banks, and the remaining 93% of eroded sediment was derived from hillslopes. Within the channelized portion of the watershed erosion/deposition could be generally described with topographic metrics used in a stream power equation. Observed sediment volumes were compared with four empirical models and one process-based model. We found that the best predictions of sediment volume were obtained from an empirical model developed in the same physiographic region. Moreover, this study showed that post-wildfire erosion rates in the San Gabriel Mountains attain the same magnitude as millennial time scale bedrock erosion rates.

California↗

Vegetation affects timing and location of wetland methane emissions

Common assumptions about how vegetation affects wetland methane (CH) flux include acting as conduits for CH release, providing carbon substrates for growth and activity of methanogenic organisms, and supplying oxygen to support CH oxidation. However, these effects may change through time, especially in seasonal wetlands that experience drying and re-wetting, or change across space, dependent on proximity to vegetation. In a mesocosm study, we assessed the impacts of on CH flux using clear flux-chamber measurements directly over plants (‘whole-plant’), adjacent to plants (where roots were present but no stems; ‘plant-adjacent’), and plant-free soils (‘control’). During the establishment phase of the study (first 30-days), the whole-plant treatment had ~5-times higher CH flux rates (51.78±8.16 mg-C md) than plant-adjacent or control treatments, which was primarily due to plant-mediated transport, with little contribution from diffusive-only flux. However, high fluxes from whole-plants depleted porewater CH concentrations both directly below whole-plant and in neighboring plant-adjacent treatments, while controls accumulated a highly concentrated reservoir of porewater CH. When the water table was drawn down to simulate seasonal drying, reserve porewater CH from control soil was released as a pulse, equaling the earlier higher CH emissions from whole-plants. Plant-adjacent treatments, which had neither plant-mediated CH transport nor a concentrated reservoir of porewater CH, had low CH flux throughout the study. Our findings indicate that in seasonal wetlands, vegetation affects the timing and location of CH emissions. These results have important mechanistic and methodological implications for understanding the role of vegetation on wetland CH flux.

Journal of Geophysical Research: Biogeosciences↗

Integrating perspectives to understand lake ice dynamics in a changing world

Ice cover plays a critical role in physical, biogeochemical, and ecological processes in lakes. Despite its importance, winter limnology remains relatively understudied. Here, we provide a primer on the predominant drivers of freshwater lake ice cover and the current methodologies used to study lake ice, including in situ and remote sensing observations, physical based models, and experiments. We highlight opportunities for future research by integrating these four disciplines to address key knowledge gaps in our understanding of lake ice dynamics in changing winters. Advances in technology, data integration, and interdisciplinary collaboration will allow the field to move toward developing global forecasts of lake ice cover for small to large lakes across broad spatial and temporal scales, quantifying ice quality and ice thickness, moving from binary to continuous ice records, and determining how winter ice conditions and quality impact ecosystem processes in lakes over winter. Ultimately, integrating disciplines will improve our ability to understand the impacts of changing winters on lake ice.

Journal of Geophysical Research: Biogeosciences↗

Getting to the root of plant‐mediated methane emissions and oxidation in a thermokarst bog

Vascular plants are important in the wetland methane cycle, but their effect on production, oxidation, and transport has high uncertainty, limiting our ability to predict emissions. In a permafrost‐thaw bog in Interior Alaska, we used plant manipulation treatments, field‐deployed planar optical oxygen sensors, direct measurements of methane oxidation, and microbial DNA analyses to disentangle mechanisms by which vascular vegetation affect methane emissions. Vegetation operated on top of baseline methane emissions, which varied with proximity to the thawing permafrost margin. Emissions from vegetated plots increased over the season, resulting in cumulative seasonal methane emissions that were 4.1–5.2 g m −2 season −1 greater than unvegetated plots. Mass balance calculations signify these greater emissions were due to increased methane production (3.0–3.5 g m −2 season −1 ) and decreased methane oxidation (1.1–1.6 g m −2 season −1 ). Minimal oxidation occurred along the plant‐transport pathway, and oxidation was suppressed outside the plant pathway. Our data indicate suppression of methane oxidation was stimulated by root exudates fueling competition among microbes for electron acceptors. This contention is supported by the fact that methane oxidation and relative abundance of methanotrophs decreased over the season in the presence of vegetation, but methane oxidation remained steady in unvegetated treatments; oxygen was not detected around plant roots but was detected around silicone tubes mimicking aerenchyma; and oxygen injection experiments suggested that oxygen consumption was faster in the presence of vascular vegetation. Root exudates are known to fuel methane production, and our work provides evidence they also decrease methane oxidation.

Alaska↗

Investigating vegetation responses to underground nuclear explosions through integrated analyses

Vegetation has the potential to respond to underground nuclear explosions, yet these links have not been fully explored. Given the lack of previously described signatures, the changes in vegetation are possibly subtle. The integration of multiple different data streams is potentially a useful approach to improve signal detection. Here, we investigate whether semi-arid vegetation growth patterns responded to eight legacy underground nuclear tests at the Nevada National Security Site in southern Nevada, USA. We tested for spatial and temporal changes in vegetation cover, tree growth patterns, and tree leaf spectral properties using ground-based measurements, including those from tree-rings and hyperspectral surface vegetation reflectance, as well as space-based measurements of Normalized Difference Vegetation Index (NDVI) from Landsat. Multiple data streams suggest a localized (<1.2 km) spatial pattern whereby tree growth is enhanced closer to the source of the underground test relative to sites further away. We also observed a more regional (>1.2–9 km) pattern whereby tree growth is suppressed coincident with a drought beginning 1 year before the 1989 tests, but continuing in the 5 years following the tests, which is anomalous relative to what is expected based on the response of tree growth to previous droughts. Quantification of the relative effects of the tests on vegetation remains a challenge due to the coincident drought and the potential for other disturbances to have impacted tree growth at this time, but the integration of these data reveals a more nuanced growth response than any other one data set indicates alone.

Nevada↗

Long-term carbon sinks in marsh soils of coastal Louisiana are at risk to wetland loss

Coastal marshes are essential habitats for soil carbon accumulation and burial, which can influence the global carbon budget. Coastal Louisiana has extensive marsh habitats (fresh, intermediate, brackish, and saline) where soil cores were collected to a depth of 100 cm at 24 sites to assess long-term carbon accumulation and burial rates. Select soil depth intervals were analyzed for bulk density, total carbon, and radionuclide ( 137 Cs and 210 Pb) dating. Marsh habitat maps (years 1949–2013) were also used to determine the most frequently occurring habitat at each field site. Over 5 decades, half of the sites transitioned between marsh habitats at least once. Saline marshes tended to have lower mean total carbon density (0.04 ± 0.002 g cm −3 ) and lower mean long-term total carbon accumulation rates (211 ± 46 g TC m −2 yr −1 , n = 5, based on 210 Pb) compared to the other marsh habitats. Using marsh habitat specific accumulation rates and area, the total carbon burial rate for coastal Louisiana in year 2013 was estimated at 4.3 Tg TC yr −1 which accounts for about 5%–21% of the estimated tidal wetland burial rate globally. Historically, about 1.0 Tg TC yr −1 was lost due to reduced marsh area from 1949 to 2013. With no coastal restoration activity, the predicted wetland loss over the next 50 years could reduce carbon burial in coastal Louisiana to 2.1 Tg TC yr −1 , a reduction of about 50% from the year 2013 rate, with potential to significantly alter the global carbon budget.

Louisiana↗

Carbon fluxes and microbial activities from boreal peatlands experiencing permafrost thaw

Permafrost thaw in northern ecosystems may cause large quantities of carbon (C) to move from soil to atmospheric pools. Because soil microbial communities play a critical role in regulating C fluxes from soils, we examined microbial activity and greenhouse gas production soon after permafrost thaw and ground collapse (into collapse-scar bogs), relative to the permafrost plateau or older thaw features. Using multiple field and laboratory-based assays at a field site in interior Alaska, we show that the youngest collapse-scar bog had the highest CH 4 production potential from soil incubations, and, based upon temporal changes in porewater concentrations and 13 C-CH 4 and 13 C-CO 2 , had greater summer in situ rates of respiration, methanogenesis, and surface CH 4 oxidation. These patterns could be explained by greater C and N availability in the young bog, while alternative terminal electron accepting processes did not play a significant role. Field diffusive CH 4 fluxes from the young bog were 4.1 times greater in the shoulder season and 1.7–7.2 times greater in winter relative to older bogs, but not during summer. Greater relative CH 4 flux rates in the shoulder season and winter could be due to reduced CH 4 oxidation relative to summer, magnifying the importance of differences in production. Both the permafrost plateau and collapse-scar bogs were sources of C to the atmosphere due in large part to winter C fluxes. In collapse scar bogs, winter is a critical period when differences in thermokarst age translates to differences in surface fluxes.

Alaska↗

Lateral carbon exports from drained peatlands: An understudied carbon pathway in the Sacramento-San Joaquin Delta, California

Degradation of peatlands via drainage is increasing globally and destabilizing peat carbon (C) stores. The effects of drainage on the timing and magnitude of lateral C losses from degraded peatlands remains understudied. We measured spatial and temporal variability in lateral C exports from three drained peat islands in the Sacramento‐San Joaquin Delta in California across the 2017 and 2018 water years using measurements of dissolved inorganic C (DIC), dissolved organic C (DOC), and suspended particulate organic C (POC) concentration combined with discharge. These measurements were supplemented with stable isotope data (δ 13 C‐DIC, δ 13 C‐POC, δ 15 N‐PON, and δ 2 H‐H 2 O values) to provide insight into hydrological and biogeochemical controls on lateral C exports from drained peatlands. Drainage DOC and DIC concentrations were seasonally variable with the highest values in the winter rainy season, when discharge was also elevated. Seasonal differences in the mobilization of dissolved C appeared to result from changing water sources and water table levels. Peat island drainage C contributions to surrounding waterways were also greatest during the winter. Although temporal variability in C cycling processes and trends were generally similar across islands, baseline drainage DIC, DOC, and POC concentrations were spatially variable, likely a result of sub‐island‐scale differences in soil organic matter content and hydrology. This spatial variability complicates system‐wide assessments of C budgets. Net lateral C exports were water year dependent and comparable to previously published vertical C emission rates for this system. This work highlights the importance of including lateral C exports from drained peatlands in local and regional C budgets.

California↗

Soil respiration response to rainfall modulated by plant phenology in a montane meadow, East River, Colorado, USA

Soil respiration is a primary component of the terrestrial carbon cycle. However, predicting the response of soil respiration to climate change remains a challenge due to the complex interactions between environmental drivers, especially plant phenology, temperature, and soil moisture. In this study, we use a 1‐D diffusion‐reaction model to calculate depth‐resolved CO 2 production rates from soil CO 2 concentrations and surface efflux observations in a subalpine meadow in the East River watershed, CO. Modeled rates are compared to in situ soil temperature and moisture conditions and MODIS satellite enhanced vegetation index (EVI) representing plant phenology across three hydrologically distinct growing seasons from 2016–2018. While soil respiration correlated with temperature on diel timescales ( p < 0.05), seasonal variability was dominated by soil moisture and plant phenology ( p < 0.05). We observed significant respiration increases in response to precipitation events; however, magnitude and duration were significantly higher in 2017 than 2016 despite similar wetting characteristics. Based on MODIS EVI, we suggest that the respiration response to rainfall is controlled by plant phenology, which in turn reflects the capacity of plants to respond to precipitation via increased photosynthesis and autotrophic respiration, behavior that is not captured in typical soil respiration pulse models. Projected changes in montane climate such as earlier snowmelt and prolonged fore‐summer drought may decrease soil respiration fluxes by decreasing the overlap between peak productivity and the summer monsoon. Finally, we observed significant late season CO 2 fluxes from the deep subsoil (>165 cm) that support growing evidence for the importance of subsoil processes in driving integrated respiration fluxes.

Colorado↗

Seasonal precipitation and soil moisture relationships across forests and woodlands in the southwestern United States

Precipitation [P: mm] controls forest and woodland dynamics in the southwestern United States (SWUS) by altering soil moisture [ θ : mm 3 mm −3 ] availability, but the influence of P on θ is complex, varying across space and time. We evaluated seasonal P and θ relationships at shallow (0‐20 cm) and intermediate (50 cm) soil depths for 9 semiarid forest and woodland sites (56 total years), which comprised 3 elevation gradients in the SWUS. We developed time‐varying definitions of winter (snow accumulation), spring (moisture recharge), and summer (moisture deficit), and determined how these sites exhibited similar P influence on θ across depths in the soil profile, between seasons, and in seasons with above‐ and below‐average P. Higher elevation sites ( > 2800 m) experienced greater winter P, longer springs, and shorter summers compared to lower elevation sites ( < 2500 m). Seasons with above‐ and below‐average P reduced elevation‐associated differences. θ at 0‐20 cm was generally decoupled from θ at 50 cm in seasons with average and below‐average P, imparted by differences in spring and summer rainfall versus winter snowfall. Notably, across‐season influence of θ (e.g. a season's similarity to subsequent seasons) was high when the first season experienced above‐ or below‐average P, and the subsequent season experienced average P, illustrating an important temporal connection initiated by wet and dry conditions. These results illustrate similarities in P‐ θ relationships across widely differing ecosystems in the SWUS, and elucidate how these relationships may be altered in a changing climate.

Journal of Geophysical Research Biogeosciences↗

Warmer winters increase the biomass of phytoplankton in a large floodplain river

Winters are changing rapidly across the globe but the implications for aquatic productivity and food webs are not well understood. In addition, the degree to which winter dynamics in aquatic systems respond to large-scale climate versus ecosystem-level factors is unclear but important for understanding and managing potential changes. We used a unique winter data set from the Upper Mississippi River System to explore spatial and temporal patterns in phytoplankton biomass (chlorophyll a , CHL) and associated environmental covariates across 25 years and ∼1,500 river km. To assess the role of regional climate versus site-specific drivers of winter CHL, we evaluated whether there were coherent long-term CHL dynamics from north to south and across lotic-lentic areas. We then estimated the degree to which these patterns were associated with climate variability (i.e., the Multivariate El Nino-Southern Oscillation Index), winter severity (freezing degree days), river discharge, or site-specific environmental variables (ice depth, snow depth, and nutrient concentrations). We found that winter CHL was typically highest in ice-free reaches and backwater lakes, occasionally exceeding summer values. We did not find highly synchronous CHL dynamics across the basin, but instead show that temporal trends were independent among river reaches and lotic-lentic areas of the river. Moreover, after accounting for these spatial dynamics, we found that CHL was most responsive to winter air temperature, being consistently higher in years with warmer winters across the basin. These results indicate that although productivity dynamics are highly dynamic within large river ecosystems, changes in the duration and severity of winter may uniformly increase wintertime productivity.

Minnesota, Missouri, Illinois, Iowa, Wisconsin↗

Quantifying eruptive and background seismicity, deformation, degassing, and thermal emissions at volcanoes in the United States during 1978–2020

An important aspect of volcanic hazard assessment is determination of the level and character of background activity at a volcano so that deviations from background (called unrest) can be identified. Here, we compile the instrumentally recorded eruptive and noneruptive activity for 161 US volcanoes between 1978 and 2020. We combine monitoring data from four techniques: seismicity, ground deformation, degassing, and thermal emissions. To previous work, we add the first comprehensive survey of US volcanoes using medium-spatial resolution satellite thermal observations, newly available field surveys of degassing, and new compilations of seismic and deformation data. We report previously undocumented thermal activity at 30 volcanoes using data from the spaceborne ASTER sensor during 2000–2020. To facilitate comparison of activity levels for all US volcanoes, we assign a numerical classification of the Activity Intensity Level for each monitoring technique, with the highest ranking corresponding to an eruption. There are 96 US volcanoes (59%) with at least one type of detected activity, but this represents a lower bound: For example, there are 12 volcanoes where degassing has been observed but has not yet been quantified. We identify dozens of volcanoes where volcanic activity is only measured by satellite (45% of all thermal observations), and other volcanoes where only ground-based sensors have detected activity (e.g., all seismic and 62% of measured degassing observations). Our compilation provides a baseline against which future measurements can be compared, demonstrates the need for both ground-based and remote observations, and serves as a guide for prioritizing future monitoring efforts.

Journal of Geophysical Research↗

Exploring GPS observations of postseismic deformation following the 2012 MW7.8 Haida Gwaii and 2013 MW7.5 Craig, Alaska Earthquakes: Implications for viscoelastic Earth structure

The Queen Charlotte-Fairweather Fault (QC-FF) system off the coast of British Columbia and southeast Alaska is a highly active dextral strike-slip plate boundary that accommodates ∼50 mm/yr of relative motion between the Pacific and North America plates. Nine M W ≥ 6.7 earthquakes have occurred along the QC-FF system since 1910, including a M S (G-R) 8.1 event in 1949. Two recent earthquakes, the October 28, 2012 Haida Gwaii ( M W 7.8) and January 5, 2013 Craig, Alaska ( M W 7.5) events, produced postseismic transient deformation that was recorded in the motions of 25 nearby continuous Global Positioning System (cGPS) stations. Here, we use 5+ yr of cGPS measurements to characterize the underlying mechanisms of postseismic deformation and to constrain the viscosity structure of the upper mantle surrounding the QC-FF. We construct forward models of viscoelastic deformation driven by coseismic stress changes from these two earthquakes and explore a large set of laterally heterogeneous viscosity structures that incorporate a relatively weak back-arc domain; we then evaluate each model based on its fit to the postseismic signals in our cGPS data. In determining best-fit model structures, we additionally incorporate the effects of afterslip following the 2012 event. Our results indicate the occurrence of a combination of temporally decaying afterslip and vigorous viscoelastic relaxation of the mantle asthenosphere. In addition, our best-fit viscosity structure (transient viscosity of 1.4–2.0 × 10 18 Pa s; steady-state viscosity of 10 19 Pa s) is consistent with the range of upper mantle viscosities determined in previous studies of glacial isostatic rebound and postseismic deformation.

Alaska, British Columbia↗

The proliferation of induced seismicity in the Permian Basin, Texas

The Permian Basin has a long history of induced earthquakes, but the seismicity rates have increased dramatically over the past two decades and included a M W 5.0 likely induced by wastewater disposal (WD) in March 2020. A detailed characterization of the proliferation of seismicity in the Permian Basin throughout this time period is needed for improving the scientific understanding of the mechanisms responsible and for mitigating future seismic hazard. Due to a sparse regional seismic network before the advent of Texas Seismological Network in 2017, we characterize seismicity using the 10-station TXAR array that is 100s of km away from most of the seismicity, with the objective of improving upon the substantial contributions from previous work. By exploiting the nature of waveform similarity, we detect events with template matching, performing a quantitative analysis of spatially varying detection capabilities throughout the study area. From an initial catalog of 10,753 events, we identify 45,009 earthquakes and 10,208 quarry blasts. Using our catalog of earthquakes, we improve epicentral locations, compare relative magnitude techniques, and associate earthquakes to WD or hydraulic stimulations. We further use our earthquake catalog to investigate the relationship between seismicity and human activities near the city of Pecos, Texas. Through a comparison of our earthquake catalog with industrial records, we determine that the vast majority seismicity near Pecos, Texas, since 2000 is likely induced by an increase of WD at wells injecting at depths greater than 1.5 km.

Texas↗