USGS Science⌕ Search

SEARCH · USGS Science

Results for “Scientific Data”

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.

939 records · Page 49Linked to original sources
USGS2025

An empirical Green’s function approach for isolating directivity effects in earthquake ground-motion amplitudes

In this study, we apply an empirical Green’s function (eGf) method within a ground‐motion modeling framework to mitigate trade‐offs between source, path, and site effects. Many physical processes contribute to spatial variations in observed ground motions, including earthquake radiation pattern, directivity, variable path attenuation, and site effects. Current nonergodic ground‐motion models use spatially varying coefficients for path and site effects, but they do not address trade‐offs with complex earthquake source effects. To quantify the influence of directivity on ground‐motion amplitudes, we use records from multiple smaller earthquakes with epicenters near that of a larger event. We use these small magnitude events as eGfs and estimate repeatable path and site effects at individual stations, assuming that the average adjustments are not controlled by directivity. We adjust residuals from the larger earthquake using the eGf terms, isolating effects related to the rupture. This method clearly enhances the observed broadband directivity observed in the 2022 M 5.1 and 2007 M 5.4 Alum Rock earthquake ground motions, reinforcing the conclusion that their ruptures were unilateral. For the 2004 M 6.0 Parkfield earthquake, we find a bilateral rupture model better fits the data because variations in rupture velocity, slip rate, and slip distribution seem to have a stronger effect on the ground motions than rupture direction alone. Applying eGf adjustments reduces the standard deviation of the rupture models over the three earthquakes by 32% on average and by up to 57% for the 2022 Alum Rock earthquake, confirming we have effectively removed repeatable effects related to the wave propagation path and site response. We propose a novel measure of the frequency‐dependent directivity amplification strength as the reduction in ground‐motion residual variability gained by fitting a directivity model; for the three earthquakes considered, this parameter varies between 25% and 75%, indicating that directivity can strongly influence ground motions and should be considered in ground‐motion modeling.

California↗
USGS2025

Alteration mapping in granitic gneiss using handheld geophysical and geochemical instruments: Implications for iron oxide-apatite and rare earth elements exploration

The Adirondack Mountains of New York, U.S.A. contain iron oxide-apatite (IOA) mineral deposits with variable concentrations of rare earth elements (REE). The IOA mineral deposits are typically hosted in the Lyon Mountain Granite Gneiss and are spatially correlated with extensive Na metasomatism (albitization) of the surrounding country rocks, although some mineral deposits also occur in metagabbro, paragneiss, and anorthosite hosts. The location of albitization zones is key to finding new IOA mineral deposits and better understanding the mineralization processes associated with their genesis. However, the Na alteration zones are generally not visibly identifiable in outcrop or hand specimens because the color and textural changes are subtle and are thus difficult to map. Here, we discuss the results from testing two spectroscopic techniques (handheld, field gamma-ray spectroscopy, and portable X-ray fluorescence (pXRF)) to characterize albitization zones surrounding IOA mineral deposits. The gamma-ray and pXRF spectrometers successfully distinguished the relative intensity of alteration in the Lyon Mountain Granite Gneiss based on the K abundance. The measured K content decreases towards the IOA mineral deposits, and the magnitude and width of the gradient are similar along the strike of each mineral deposit. Elevated Th and U values are present in host rocks adjacent to REE-bearing IOA mineral deposits. The pXRF and gamma-ray spectrometer K values are remarkably consistent with laboratory-based, whole-rock XRF compositional data and, therefore, useful for semi-quantitative analysis. Notably, albitization aureoles are consistent around REE-rich and REE-poor mineral deposits, suggesting that saline fluids are essential for the petrogenesis of IOA mineral deposits but may not be critical to REE mineralization. Ultimately, this study demonstrates the utility of handheld gamma-ray and pXRF spectrometry for identifying otherwise cryptic albitization gradients associated with IOA mineral deposits in granitic gneiss.

New York↗
USGS2025

Magnitude, depth and methodological variations of spectral stress drop within the SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest Earthquake Sequence

We present the first ensemble analysis of the 56 different sets of results submitted to the ongoing Community Stress Drop Validation Study using the 2019 Ridgecrest, California, earthquake sequence. Different assumptions and methods result in different estimation of the source contribution to recorded seismograms, and hence to the source parameters (principally corner frequency, f c ⁠ , spectral stress drop, Δσ, and seismic moment, M 0 ⁠ ) obtained from modeling calculated source spectra. For earthquakes smaller than magnitude (M) 2.5 there is negligible correlation between the f c values obtained by different studies, implying that no present method is reliable using available data. For larger magnitude events, correlation between f c measurements of different studies, within even a small M range is always higher than spectral ⁠Δσ , because the f c measurements simply reflect the underlying physical decrease in f c with increasing M. We model the observed trends of submitted f c with both magnitude and depth. Most methods report an increase in spectral Δσ with M, although a magnitude‐invariant spectral Δσ is within the confidence limits. The depth dependence is smaller and depends on whether a study allows attenuation to vary with source depth; a combination of depth‐dependent attenuation correction, and depth‐dependent shear‐wave velocity can compensate for reported depth trends. We model the submitted values to remove differing M and depth variation to investigate the relative interevent variability. We find consistent relative variation between individual events, and also lower relative spectral Δσ in the northwest of the aftershock sequence, and higher on the cross fault and in the region of main fault intersection. This large‐scale comparison implies that absolute spectral Δσ estimates are dependent on the methods used; studies of different regions or using different methods should not be directly compared and improved constraints on path and site corrections are needed to resolve these absolute spectral Δσ differences.

California↗
USGS2025

Assessment of western Oregon debris-flow hazards in burned and unburned environments

In the steep and mountainous environment of western Oregon, debris flows pose a considerable threat to property, infrastructure and life. Wildfire is commonly known to increase the susceptibility of steep slopes to debris flows, but the extent of this process in the western Cascades is not well understood. The US Geological Survey (USGS) currently estimates postfire debris-flow likelihood and triggering rainfall thresholds using a model calibrated to a southern California inventory of debris flows generated by excess runoff within the first year after fire. Because of a lack of available data, this model has not been tested in western Oregon, or in locations where postfire debris flows initiate via other mechanisms (e.g., shallow landslides or in-channel failures). Using repeat field observations and aerial imagery, we developed two new debris-flow inventories within and adjacent to the perimeters of five 2020 wildfires in western Oregon: Archie Creek, Holiday Farm, Beachie Creek, Lionshead and Riverside. The first inventory focuses on postfire debris flows (2020–2022); the second focuses on debris flows prior to fires (1995–2020). Our inventories of prefire and postfire debris flows were used to document initiation mechanisms in Oregon's western Cascades and to evaluate the effects of wildfire. We found that wildfire changed the distribution of debris-flow initiation mechanisms in the western Cascades. After the wildfires, annual rates of runoff-generated debris flows increased by 22% and the number of shallow landslide-initiated debris flows decreased by 17% relative to before the wildfires. Despite this shift, shallow landsliding was the dominant debris-flow initiation mechanism in both unburned and burned environments. We found the performance of the current USGS debris-flow likelihood model was degraded relative to other previously tested locations across the intermountain western United States. Our results highlight the need for improved postfire hazard assessment in western Oregon based on regional model calibration that is tuned to the dominant debris-flow initiation mechanisms.

Oregon↗
USGS2026

Hydrogeology, groundwater salinity distributions, and assessment of the effect of oil-production activities on groundwater in the Midway Valley area, western Kern County, San Joaquin Valley, California

This study seeks to determine the effects of oil field produced water disposal operations and well mechanical integrity issues on groundwater quality in oil fields in the southwest San Joaquin Valley, California. Whereas previous studies used groundwater wells to study shallow aquifers outside the oil fields, this study demonstrates that future approaches may use oil well geophysical logs to map groundwater head gradients, create salinity profiles and document changes in salinity over time in oil field areas with sparse groundwater well data and at depths greater than 330 m. We also incorporate an analysis of well histories to determine potential effects of compromised wellbore seals on changes in aquifer quality that cannot be explained by water disposal practices. Water quality in the aquifers is naturally brackish across most of the area, with better quality groundwater occurring in the eastern part. Geophysical logs are used to determine salinity variations within aquifers including the depth at which TDS exceeds 10,000 mg/L. This depth ranges from 366 m in the northwest to approximately 1,500 m in the southeast. Oil well porosity logs are used to determine water table elevations. These logs indicate the water table slopes south-southeast, showing the predominant groundwater flow direction is from oil field disposal areas toward better quality groundwater east of the oil fields. Geophysical logs show formation resistivity near some disposal facilities has decreased over time, indicating the salinity of the aquifer has increased due to disposal of saline produced water in injection wells and ponds. Oil well history analysis suggests that increased salinity over time in water-saturated sand intervals >1.5 km from disposal facilities may be caused by mechanical failures and/or incomplete borehole seals in poorly constructed or abandoned wellbores prevalent throughout the study area—particularly wells drilled prior to 1930.

California↗
USGS2025

Spatiotemporal variations in strain release and seismic rupture in multifault systems: An example from Panamint Valley, southeastern California

Geometrically complex, multifault ruptures have been observed in recent, damaging earthquakes in southeastern California, sparking renewed efforts to identify physical conditions that promote or inhibit fault discontinuity-spanning coseismic ruptures. The likelihood of ruptures propagating across fault discontinuities is thought to be partly controlled by fault geometries, rupture direction, and the history of strain release. However, these parameters vary in space and time over multiple earthquake cycles, making it difficult to forecast the likelihood that an earthquake on one fault will trigger rupture on a nearby fault. Here we use tectono-geomorphic mapping of a geometrically complex fault zone in Panamint Valley, southeastern California, to assess spatiotemporal variations of paleo-rupture patterns and geometries of fault discontinuities over multiple earthquake cycles. First, we identify ten generations of late Pleistocene to Holocene alluvium using geomorphic parameters and luminescence dating to constrain ages of alluvium and bracket late Holocene earthquake timing. Then, we quantify slip kinematics using high-resolution structure from motion digital surface models. We find the Panamint Valley transtensional relay (PVTR) hosted four late Holocene earthquakes, bracketed to ~5.8–3.4 ka, ~3.8–2.2 ka, ~2.4–0.6 ka, and ~0.64–0.16 ka, with ~0.6–1.1 m of slip per event, correlative to Mw ≈ 6.7–6.9 earthquakes. Additionally, we find similarities in earthquake timing on the Ash Hill, PVTR, and Panamint Valley faults and similarities in the slip magnitude and slip kinematics between the Ash Hill and PVTR faults, implying that the PVTR may co-rupture with nearby faults. Paleo-rupture patterns indicate that seismogenic strain transfer may occur through the PVTR, along different combinations of fault segments and jump distances, over multiple earthquake cycles. These data highlight the utility of tectono-geomorphic mapping in evaluating paleo-rupture patterns and suggest that the PVTR may act to propagate and/or arrest rupture between the Ash Hill and Panamint Valley faults.

California↗
USGS2026

Vegetation cover and composition in environments surrounding uranium mines in the Grand Canyon ecosystem, Northern Arizona

Mining uranium from breccia-pipe deposits in the greater Grand Canyon region has occurred since the mid-1900s. However, possible ecosystem contamination with harmful levels of radionuclides may have occurred due to mining activities in the 21st century. In response, a 20-year Federal moratorium on new mining claims in the Grand Canyon watershed was initiated in 2012, to allow time to evaluate the potential effects of uranium exploration and mining on human health, wildlife, and water resources. This moratorium, nor the 2023 designation of the “Baaj Nwaavjo I’tah Kukveni–Ancestral Footprints of the Grand Canyon National Monument,” precludes operation or development of mining claims predating 2012. Vegetation is a core ecosystem component that may be affected by uranium mining (for instance, through uptake and storage of radionuclides from the air or soil) or may act as a vector of exposure to wildlife, livestock, and humans (for instance, via their consumption of contaminated plant tissues). To provide baseline information about the plant communities associated with uranium mines in the Grand Canyon region, the U.S. Geological Survey surveyed an approximately 200-meter-wide buffer surrounding four breccia-pipe deposits, each in a unique stage of mine development, and at one reference area (a livestock water tank) that underwent ground disturbance but contains no mineral deposits. We sectioned the buffer zones into 0.65–4.52 hectare plots, within which we (1) inventoried all plant species, (2) measured percent cover of plant species, plant functional groups, and ground surface types (dark cyanobacteria, lichen, moss, bedrock, rock, embedded litter, duff, plant bases, and bare soil) using line-point intercept, and (3) measured length and frequency of gaps between perennial plant canopies using canopy gap intercept. We found that plant composition at the mines and the reference area differed from one another but were all characteristic of expected regional vegetation patterns. We provide this data summary as potential baseline information for future research and management efforts.

Arizona↗
USGS2024

Limited evidence of late Quaternary tectonic surface deformation in the eastern Tennessee seismic zone, USA

The ~300-km-long eastern Tennessee seismic zone (ETSZ), USA, is the second-most seismically active region east of the Rocky Mountains. Seismicity generally occurs below the Paleozoic fold-and-thrust belt within the Mesoproterozoic basement, at depths of 5–26 km, and earthquake magnitudes during the instrumental record have been moment magnitude ( M w )≤4.8. Evidence of surface deformation may not exist or be difficult to detect because of the vegetated and soil-mantled landscape, landslides, locally steep topography, anthropogenic landscape modification, or long, irregular recurrence intervals between surface-rupturing earthquakes. Despite the deep seismicity, analog models indicate that accumulation of strike-slip or oblique-slip displacement at depth could be expected to propagate upward through the Paleozoic section, producing a detectable surficial signal of distributed faulting. To identify potential surface deformation, we interrogated the landscape at different spatial scales. We evaluated morphotectonic and channel metrics, such as channel sinuosity and catchment-scale hypsometry. Additionally, we mapped possible fault-related topographic features on 1-m lidar. Finally, we integrated our observations with available bedrock and Quaternary surficial mapping and subsurface geophysical data. At a regional scale, most morphotectonic and channel metrics have a strong lithologic control. Within smaller regions of similar lithology, we observe changes in landscape metrics like channel sinuosity and catchment-scale hypsometry that spatially correlate with new lineaments identified in this study and previously mapped east–west Cenozoic faults. These faults have apparent left-lateral offsets, are optimally oriented to slip in the current stress field, and match kinematics from recent focal mechanisms, but do not clearly preserve evidence of late Pleistocene or Holocene tectonic surface deformation. Most newly mapped lineaments might be explained by either tectonic or non-tectonic origins, such as fluvial or karst processes. We also re-evaluated a previously described paleoseismic site and interpret that the exposure does not record evidence of late Pleistocene faulting but instead is explained by fluvial stratigraphy.

Tennessee↗
USGS2025

Broadband stochastic simulation of earthquake ground motions with multiple strong phases with an application to the 2023 Kahramanmaraş, Turkey (Türkiye), earthquake

Stochastic ground motion simulation models are often less accurate at lower frequencies than at higher frequencies when fitting recorded data unless supplemented by a deterministic forward directivity velocity pulse model. Moreover, time-modulated stochastic models, which adjust ground motion amplitudes over time, typically use functions that fail to capture multiple strong-motion phases. The February 2023 Turkey (Türkiye) earthquake exhibited diverse recordings, including near-fault and far-field motions with pulse-like and non-pulse-like characteristics, along with single and multiple strong-motion phases. To better represent such a diverse set of recordings, this study enhances a fully non-stationary site-based stochastic model without combining it with a deterministic model. Improvements include a new band-pass filter with upper- and lower-frequency limits, which refines the representation of the low-frequency content. Moreover, a time-modulating function that can represent energy arrival in multiple strong phases is introduced. The reference model’s parameters are identified by fitting to the energy content, zero-level crossings, and cumulative counts of positive-minima and negative-maxima of a target accelerogram. This fitting procedure is modified to address the increased number of parameters. These improvements broaden the reference model’s applicability while preserving its simplicity, a key aspect appealing to engineering practitioners. The improved model’s applicability is demonstrated by simulating a dataset from the February 2023 Türkiye earthquake, and the accuracy is tested using a pulse-like Next Generation Attenuation Relationships for Western United States dataset. Validations are performed based on total energy, zero-level crossings, Fourier amplitude spectrum, elastic response spectra, and peak ground motion parameters. Validations are performed schematically in the time and frequency domains and quantitatively using goodness-of-fit scores, various validation-metrics errors, and inter-period correlations. Overall, the improved stochastic model can effectively simulate a set of diverse ground motion recordings, including near-fault pulse-like records, records with multiple strong phases, and far-field motions across a broad frequency range.

Kahramanmaraş↗
USGS2026

Bathymetric survey and storage capacity of Upper Lake Mary near Flagstaff, Arizona in 2024

The U.S. Geological Survey (USGS), in cooperation with the city of Flagstaff, collected bathymetric, light detection and ranging (lidar), and land-survey data of Upper Lake Mary in Arizona during the months of April and October 2024. The city of Flagstaff uses a combination of groundwater from well fields throughout the Flagstaff area and surface water, mainly from Upper Lake Mary, for its potable water supply. The purpose of the survey is to update previous surveys using new technology and compare the results to previous surveys to determine if there was a decrease in storage capacity that could affect the city’s water supply. The lakebed was mapped in April 2024 using a vessel equipped with a multibeam echosounder (MBES) and mobile lidar scanner with positioning captured using a real-time kinematic global navigation satellite system (RTK GNSS) base and receivers. In October 2024, areas of the reservoir that were too shallow for the boat and shoreline that were not captured by the vessel-based lidar were surveyed on foot using hand-held RTK GNSS receivers. At full pool (spillway elevation of 6,831.82 feet above NAVD 88 [2,082.34 meters (m)], Upper Lake Mary has a storage capacity of 16,449.80 acre-feet (20,290,611.73 cubic meters) and a surface area of 953.57 acres (3,860,926.075 square meters). The reservoir is 5.7 miles (9.7 kilometers) long and varies in width from 326 feet (99.36 m) near the central, narrow portion of the reservoir to 2,613 feet (796.44 m) in the upper portion. Comparisons between this survey and the previous two surveys from the 1950s and 2006 indicate no apparent decrease in reservoir area or storage capacity. Results of the 2024 survey indicate that Upper Lake Mary’s storage capacity increased by 0.9 percent from the 2006 survey and a 1.6 percent increase in surface area from the 2006 survey.

Arizona↗
USGS1976

Availability of ground water in the area surrounding the Trident submarine construction facility, Kitsap County, Washington

This report provides general information on the water resources--with emphasis on ground-water occurrence and availability--in that part of Kitsap County (referred to as Trident Impact Area) that would be most affected by the development of the Trident submarine-construction facility at Bangor. Ground water occurs principally in a lower, confined aquifer near and below sea level and an upper, unconfined aquifer above sea level; both aquifers are extensive beneath the study area and are composed of sand and gravel layers separated by poorly permeable silt, clay, and till layers. Locally throughout the area efficiently constructed wells tapping each of the two aquifers are capable of yielding 100 to 1,300 or more gallons per minute. The estimated 1970 water use in the study area averaged about 13 million gal/day (gallons per day); of this amount about 9 million gal/day came from surface-water sources--from a large reservoir outside the study area--and about 4 million gal/day came from ground water pumped from the two aquifers. Anticipated water use soon will be about 18 to 21 million gal/day; virtually all the additional quantity required (about 5 to 8 million gal/day) above present use must come from ground-water sources. Preliminary evaluation of the aquifers suggests that an additional 1.5 million gal/day can be developed from the upper aquifer and 7 million gal/day from the lower aquifer. Existing wells tapping the lower aquifer might yield additional water and increase the total yield in the area by 3.5 million gal/day, and new wells drilled in selected areas could produce an additional 3.5 million gal/day from this aquifer. However, additional, large-scale ground-water withdrawal from the lower aquifer could induce saltwater intrusion into wells situated in coastal areas. Recommended future studies include additional data collection, monitoring of chloride concentrations in water from large-yield coastal wells, confirmation of well yields and capacities, and test-well drilling and analysis.

Washington↗
USGS2024

Fairweather transform boundary Oligocene to present orogenesis: Fairweather Range vertical extrusion and rotation of the Yakutat microplate at ca. 3 Ma

Oblique-slip along transform fault boundaries is often partitioned between a strike-slip system and thrust faults that accommodate contraction. However, topography along the Yakutat-North American transform (Fairweather fault), is asymmetric with low-terrain above active thrusts on the western, Yakutat side of the transform and high topography on the continental side with peaks >4500 m (Mount Fairweather: 4671 m) to the west of the Border Ranges fault, limited recorded earthquakes >M4, and no apparent reverse faults to generate the highest terrain. In this study we compile, for the first time, published U-Pb zircon, 40 Ar/ 39 Ar and K-Ar (hornblende, muscovite, and biotite) and U-Th/He and fission-track (zircon and apatite) bedrock ages (109) from 75 samples to investigate the exhumation history of the Fairweather Range region, complemented by a published detrital sample (ZFT and AFT) and 13 new 40 Ar/ 39 Ar (hornblende, biotite, and K-feldspar) ages on 9 bedrock samples from both sides of the Fairweather fault. Additionally, we examined published seismicity and geodetic data of the Fairweather region and assessed if plate paleo-vectors correlate with the cooling history of the Fairweather Range. Cooling age, seismic, and block-motion patterns indicate the Fairweather Range has been vertically extruded between the Fairweather and the Border Ranges faults as a coherent block since ca. 25 Ma. The pre-6 Ma Pacific plate motion (N30°W) aligns with the N33°W strike of the Fairweather Fault whereas a hypothetical pre-6 Ma Yakutat microplate paleo-vector of (N39°W) does not: indicating a post-6 Ma timing for Yakutat microplate counter-clockwise rotation (9°). We infer that rotation and impingement of the Yakutat microplate along the Fairweather fault at ca. 3 Ma led to the development of the Fairweather restraining bend and increased cooling rates. The resultant thickened Fairweather welt and the ∼30 km thick southeast end of the Yakutat microplate compounded double-indenter tectonics into Alaska's southeast convergent corner

Tectonophysics↗
USGS2024

A methodology to estimate CO2 and energy gas storage resources in depleted conventional gas reservoirs

Depleted hydrocarbon reservoirs are subsurface geological structures capable of sequestering vast quantities of carbon dioxide (CO 2 ) as well as storing other energy gases for later usage, such as natural gas, and potentially hydrogen (H 2 ). Here we outline a methodology to quantify multi-gas storage resources in depleted conventional gas reservoirs for usage in assessments by the United States Geological Survey (USGS) at the scale of sedimentary basins. The methodology consists first of quantifying accessible pore volume in a depleted reservoir for natural gas storage using up to three equations. Input data are derived from commonly reported or estimated reservoir parameters and natural gas production volumes, and equations may be combined in linear models to improve pore volume estimates. Storage estimates from these equations are tested and validated for 31 reservoirs in the Michigan Basin Province, USA that were previously converted to underground gas storage facilities and have known (federally reported) natural gas storage capacities. Secondly, natural gas storage capacities can be transformed via fluid substitution calculations to estimate the storage resources for non-native fluids, applied here for, CO 2 , H 2 , and methane-H 2 blends, accounting for molecule-specific deviations from ideal gas behavior at reservoir pressures and temperatures as well as differing storage efficiencies. Importantly, the storage of non-native fluids may not be appropriate in all depleted gas reservoir settings due to potential risks like leakage, in particular in the case of H 2 storage, requiring additional knowledge of caprock sealing capacity. Given this caveat, we demonstrate the fluid substitution method for natural gas reservoirs of the Northern Niagaran Reef and Southern Niagaran Reef USGS plays in the Michigan Basin Province, as these trends of Silurian pinnacle reefs are capped with tight-sealing evaporite facies. The deterministic equations outlined from this methodology can be incorporated into future probabilistic USGS gas storage assessments for CO 2 , H 2 , and natural gas resources in the United States.

Michigan↗
USGS2026

Repeated coral bleaching events have eroded reef accretion potential in the Florida Keys

Coral bleaching is now the primary driver of coral-reef degradation globally, threatening the persistence of coral-reef functions and the invaluable ecosystem services they support. The Florida Keys have been impacted by four severe, regional-scale coral-bleaching events, which occurred in 1997–1998, 2005, 2014–2015, and 2023. Our aim in this study was to quantify the effects of those thermal-stress events on the persistence of coral-reef structures. We leveraged long-term coral-reef monitoring data to quantify the impacts of coral bleaching on coral cover, carbonate production, and reef-accretion potential throughout the Florida Keys from 1996–2024. Our study indicates that coral bleaching drove significant declines of reef-building corals, which successively diminished regional-scale reef-accretion potential (from 1.65 ± 0.33 mm y -1 in 1996) and culminated with a regional-scale transition to net reef erosion (-0.12 ± 0.08 mm y -1 ) following the 2023 event. Coral bleaching was the primary driver of declining reef-building capacity, with the cumulative impacts of the four coral-bleaching events accounting for >75% of total declines since 1996. The 1997–1998 event alone drove ~50% of overall declines, largely due to losses of the shallow-water ecosystem engineer, Acropora palmata . More moderate impacts of thermal stress on regional reef-accretion potential in 2005 (6%) and 2014–2015 (9%) could lend some support to the hypothesis that the remaining corals were more thermally tolerant; however, the unprecedented thermal stress event in 2023 caused widespread coral-bleaching and bleaching-related mortality, and caused an additional 12% of the total decline in reef-accretion potential. Stony-coral tissue loss disease accounted for an additional 6% of the decline, with the remaining 19% unexplained. Coral bleaching has already dramatically reshaped the structure and function of coral reefs in the Florida Keys over the last three decades. With the high likelihood that thermal stress will accelerate moving forward, the persistence of Florida’s essential coral-reef functions and services is uncertain.

Florida↗
USGS2025

Mapping potential sensitivity to hydrogeomorphic change in the UMRS riverscape

In 2020 the U.S. Geological Survey (USGS), as part of the Upper Mississippi River Restoration (UMRR) Program, began a new project to characterize potential hydrogeomorphic change associated with hydrogeomorphic units (HGUs) and their catenae (units linked by their association with sediment sources and flow origins). The goal of the project was to develop a geographic information system (GIS) database of HGUs for the Upper Mississippi River System (UMRS) available to both scientists and river managers working on UMRR studies and HREP planning and design studies. The characterization was based on a hydrogeomorphic change hierarchical classification developed previously for the UMRS. The products were generated with automated techniques in a GIS using systemic datasets. Landforms were mapped from the 2015 UMRS topobathymetric dataset with geomorphon (shape-based) tools tailored for the large riverscape, valley bottom environments in the UMRS. A clustering analysis was applied to the resulting landforms to identify HGUs associated with zones of perennial low flows, bankfull flows, and overbank floodplains. Catenae were assembled based on the proximity of the units to the main channel, tributary mouths, and side channels from previously published aquatic areas (USACE, 2018) coupled with least-cost flowpath linkages between potential sediment origins and planform change units developed by Rogala, Fitzpatrick, and Henderson (2020). These GIS-based analyses were successful at identifying a range of HGUs using an automated technique with available data across the entire riverscape, with emphasis on those that have the potential for hydrogeomorphic change. Most of the resulting features are depositional, as expected in a large river system. However, this is the first attempt of linking tributary inputs, side channel erosion and levee breaches with their depositional counterparts. The approach was successfully piloted in Pools 8 and 10 in the Upper Impounded Reach and Pool 14 in the Lower Impounded Reach, with next steps for application in reaches of the unimpounded section and Illinois River. This report emphasizes results from Pool 10, which was the focus of most of our attention during the pilot phase.

Illinois, Indiana, Iowa, Minnesota, Missouri, Sout↗
USGS2025

Land cover change within wetland complexes at Dixie Meadows, Churchill County, Nevada: 2015 – 2023

Dixie Meadows, Nevada, is a system of geothermal springs and seeps that feed a complex of marshes and wetland meadows that are located within lands managed by the Bureau of Land Management (BLM) and the Department of Defense (DOD). A previous U.S. Geological Survey report documented variability in satellite imagery-based land cover classifications for seven wetland complexes at near monthly time intervals between October 2015 and January 2022. This report presents additional data, extending analysis to November 2023. Land cover classifications between October 2015 and November 2023 demonstrated an association between vegetation cover characteristics and surface moisture, with Class 1 having dry, bare soil or sparse upland vegetation, Class 2 having moist, bare soil or sparse to small vegetation, Class 3 having dense green vegetation with potentially saturated soil conditions, Class 4 having a mix of shallow surface water, saturated soil, and dense green vegetation, and Class 5 having open surface water. Most of the wetland complexes occur close to spring outflows primarily within land managed by the DOD, though portions are also within BLM lands. The intervening and surrounding landscape outside of the wetland complexes assessed in this study are managed by the BLM. As a result, Class 1 land covers had the largest areal coverage for BLM managed lands. Classes 2 and 3 land covers were primarily mapped inside the wetland complexes and thus had the largest area coverage within DOD managed lands. Class 4 was almost exclusively mapped within the wetland complexes and thus was largely contained within DOD managed lands. Class 5 (open water) was exclusively mapped in and adjacent to a single wetland complex with catchment ponds on land managed by the BLM. The distribution of these land cover classes over the study period was seasonally and annually variable. Land cover areas of Classes 1 and 2 were larger during the spring months. Conversely, land cover areas of Classes 3 and 4 tended to be greatest during the summer or fall. These patterns might be influenced by differences in seasonal water sources and phenology.

Nevada↗
USGS2024

Cold blood in warming waters: Effects of air temperature, precipitation, and groundwater on Gulf Sturgeon thermal habitats in a changing climate

Objective In a changing climate, the effects of air temperature, precipitation, and groundwater on water temperature and thermal habitat suitability for Gulf Sturgeon Acipenser desotoi , listed as threatened under the U.S. Endangered Species Act, are not well understood. Hence, we incorporated these factors into thermal habitat models to forecast how Gulf Sturgeon may be affected by wide‐ranging climate change scenarios in 2024–2074. Methods Using data from the Choctawhatchee River, Florida, we developed precipitation‐ and groundwater‐corrected air–water temperature models, compared their accuracy with that of conventional air–water temperature models used in fisheries management, and projected future Gulf Sturgeon thermal habitat suitability for normal physiological functioning and fieldwork (i.e., population sampling and telemetry surgeries) in summer (May–August) under 16 climate change scenarios. Result Precipitation‐ and groundwater‐corrected models were more accurate than conventional air–water temperature models (mean improvement in adjusted R 2 = +0.45; range = +0.09 to +0.75). Water temperature was projected to warm at widely variable rates across climate change scenarios encompassing different air temperature, precipitation, and groundwater regimes. Importantly, Gulf Sturgeon summer aggregation areas were cooler and influenced more by precipitation and groundwater and less by air temperature than were non‐aggregation areas. If precipitation and groundwater—as drivers of cooling—become warm in a changing climate, summer aggregation areas were projected to exhibit thermal habitat degradation equivalent to or greater than that of non‐aggregation areas. Conclusion Our results add hydrological context to the premise that aggregation areas provide cool water and energetic savings for Gulf Sturgeon during summer, underscoring the importance of protecting these habitats through groundwater conservation, water quality monitoring, and riparian/watershed habitat management. Our findings indicate that identifying thermally appropriate times for fieldwork activities will be increasingly important and time‐restricted as climate change intensifies. However, our research provides managers with a portfolio of water temperature models and an accurate, cost‐effective, management‐relevant approach to forecasting thermal habitat conditions for Gulf Sturgeon and other species in a changing climate.

Alabama, Florida↗
USGS2024

Methodology for inclusion of produced and stored carbon dioxide in the U.S. Geological Survey Federal lands greenhouse gas inventory

The U.S. Geological Survey (USGS) has developed two new carbon dioxide (CO2) emissions and sequestration accounting methods for use in future reports. The first method is a Federal lease-produced CO2 emissions calculation for an update of the report, “Federal Lands Greenhouse Gas Emissions and Sequestration in the United States.” The methodology to incorporate Federal lease CO2 production emissions into the updated report relies on CO2 sales royalty data from the Office of Natural Resources Revenue (ONRR). The end usage points for the gas include enhanced oil recovery with CO2 (CO2-EOR), food and beverage, and chemical production. CO2-EOR is the main end point for natural CO2 production in the United States; it accounted for 94% of usage in 2022 [1]. Federal lands emissions from this sector are estimated at 460 metric tons of CO2 in 2022, a very small amount relative to most other Federal lands emissions sector estimates. The second new method, planned for a separate report, is a calculation of the geologic storage of CO2 on Federal lands. The second method estimates the CO2 stored under Federal surface lands and documents Federal climate change mitigation efforts. Currently, there is no storage of CO2 at an industrial level on Federal lands, however multiple proposals and projects are planned. This method was developed on non-Federal lands datasets in an effort to prepare for when these activities on Federal lands will require accounting. National estimates for CO2 geologic storage using this method, but without a Federal lands filtering step, totaled 8.0 million metric tons (Mt) in 2022. The two methods described here are new benchmark methods in a collection of accounting procedures to document the current state of greenhouse gas emissions and their storage on Federal lands. These benchmarks can then be used to measure any subsequent changes in emissions from or carbon storage beneath Federal lands. While the magnitude of the values is currently non-existent to small, emissions mitigation goals established by decision makers indicate that these values will grow, and their documentation will take on greater value and use.

continental United States↗