USGS ScienceSearch

SEARCH · USGS Science

Results for “Integrated Environmental Assessment and Management”

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.

1,661 records · Page 91Linked to original sources

Temperate wetlands lose climate-cooling capacity under warming

Temperate inland wetlands have long been viewed as climate regulators, with sustained CO₂ sequestration offsetting continuous CH₄ emissions over centennial–millennial timescales. However, the net climatic effect of these opposing fluxes remains uncertain under anthropogenic warming as they differ in both magnitude and atmospheric lifetime. We combine 19,000 CH₄ chamber measurements from 202 temperate inland wetlands with carbon burial rates from 16 dated sediment cores in North America’s Prairie Pothole Region to estimate greenhouse gas dynamics over two centuries. Here we show that CH₄ emissions rise by ~1.39 kg CH₄ ha⁻¹ yr⁻¹ over the historical period and are projected to approach 950 kg CH₄ ha⁻¹ yr⁻¹ under high-forcing scenarios by 2100, while CO₂ sequestration is expected to decline by ~0.3 kg CO₂ ha⁻¹ yr⁻¹, driving a shift toward persistent net CO₂-equivalent source behaviour and a weakening of the long-term cooling effect that wetlands have provided under natural conditions.

Prairie Pothole Region

Ecohydrological response of a forested headwater catchment to a flash drought in the Southeastern U.S.

Flash droughts differ from traditionally defined droughts in their rapidity of intensification and often associated high vapor-pressure deficit. These droughts can lead to declines in streamflow and water table depth and induce water stress to vegetation at a greater rate than droughts that manifest over longer periods. However, little is known regarding the response of forested environments to flash drought because most studies of impacts have been conducted in agricultural settings. In this study we investigated water-use patterns of riparian trees using sap flow methods and examined the role of groundwater as a source of moisture over three periods that were delimited by antecedent soil moisture conditions. For a longer-term perspective we also examine monthly streamflow over the 35-year record. We observed that trees at only one monitoring plot showed a decrease in water use relative to evaporative demand during a flash drought. Total reverse sap flow (flow toward the roots rather than the canopy) greatly increased during the flash drought period, suggesting the likely occurrence of hydraulic redistribution to the excessively dry soils. Over the drought period groundwater became a more dominant source of moisture for sustaining forest water use. Monthly mean streamflow during the flash drought approached levels observed in past multiyear droughts. This is the first study, to our knowledge, to specifically investigate the response of multiple water budget components to flash drought in a humid forest. As more studies are conducted, a better understanding of the range of expected responses are likely to emerge.

Georgia

Ultralong, supershear rupture of the 2025 Mw 7.7 Mandalay earthquake reveals unaccounted risk

The 28 March 2025 moment magnitude ( M w ) 7.7 earthquake in Mandalay, Burma (Myanmar), ruptured 475 kilometers of the Sagaing Fault, which was more than twice the length predicted by magnitude scaling relationships. Kinematic slip models and observation of a Rayleigh Mach wave that passed through parts of Thailand confirmed that rupture occurred at supershear velocities of greater than 5 kilometers per second. The anomalous length exposed a vast population to violent near-fault shaking. The Mandalay earthquake is a modern analog for the M w 7.9 1906 San Francisco earthquake, another atypically long and fast rupture. Probabilistic seismic hazard analyses use scaling relations that do not account for such long ruptures at moderate magnitudes. This limitation, in conjunction with a likely increased population and infrastructure exposure for atypically long ruptures, contributes to a potential mischaracterization of seismic risk.

Mandalay

Summary of annual repeat magnetotelluric surveys of the Geysers geothermal field

As part of a multi-year project, funded by the California Energy Commission, annual magnetotelluric (MT) surveys have been collected at The Geysers geothermal field in northern California with the goal of measuring temporal changes within the steam field. The repeat surveys started in 2021 and repeated a 2017 survey (Peacock et al., 2020) with further extension to the southern part of the geothermal field. Temporal variations in the MT transfer functions are observed to be spatially coherent and compartmentalized. Mapping residual phase tensor ellipses demonstrates the direction of maximum change is often aligned with existing fracture orientations. Three dimensional inversion of the MT data, using the inversion results from previous years as the starting model, indicates that the steam reservoir has generally become more resistive over time (~10%), suggesting more steam in the field. A few pockets within the steam field have become more electrically conductive over time and are collocated with injection wells, suggesting either more fluid content in those zones, less steam, or more saline fluids.

California

Prospectivity modeling of the NASA VIPER landing site at Mons Mouton near the Lunar South Pole

We use a high-resolution digital elevation model and a numerical thermal model to produce a variety of inputs for a water-ice prospectivity model for the Volatiles Investigating Polar Exploration Rover (VIPER) landing site. These input data are maps of topography, surface slope, surface aspect, surface curvature, maximum temperature, depth to ice stability, permanently shadowed regions (PSRs), distance to PSRs, and PSR density. This model predicts where water ice is most likely within the top meter of regolith, assuming plausible relationships between ice concentration and the various inputs. The model is designed to be adjusted in near-real time as data are collected during the VIPER mission. As such, it is a tool for both analyzing data from the mission as well as planning operations. Since the current model, at this point, relies only on orbital remote sensing, the final version will also be a tool to extrapolate the VIPER mission results across the lunar poles.

Planetary Science Journal

Magnitude conversion relations create substantial differences in seismic hazard models

Earthquake catalogs are essential data inputs for seismic hazard modeling. Because earthquake magnitudes are reported in a variety of types (e.g., local magnitudes and moment magnitudes), magnitude conversion relationships must be used to convert the different magnitude types present in a catalog to a uniform magnitude type to avoid biases in the hazard computation. However, these conversion relationships are often uncertain and have been shown to sometimes perform poorly. Here, we investigate the sensitivity of the gridded seismicity component of the National Seismic Hazard Model (NSHM) to the catalog conversion equations in the Eastern United States. In the 2023 NSHM, magnitudes of various types were converted to moment magnitudes using equations developed by the Central and Eastern United States Seismic Source Characterization for Nuclear Facilities (CEUS‐SSCn), based on least‐squares (LS) regressions made using data from a catalog containing events up through 2008. We recompute these equations using events in the Advanced National Seismic System Comprehensive Earthquake Catalog with multiple magnitudes from 2000 to 2023. Although we prefer the use of orthogonal regressions for our datasets, LS regressions produce broadly similar results, with both approaches exhibiting large deviations from the CEUS‐SSCn conversions, especially at smaller magnitudes. We compare the spatial distribution of annual rates using three different models: (1) the 2023 NSHM conversions, (2) our updated conversions, and (3) no conversions. We find that the choice of conversions leads to substantial differences in the rate forecasts, which can greatly impact the seismic hazard model, particularly in regions with low‐seismicity rates such as the Eastern United States, where the hazard is dominated by gridded seismicity rather than a fault model.

Seismological Research Letters

Simulated hydrologic responses to proposed wastewater-returnflow scenarios in Falmouth, Massachusetts

The Cape Cod aquifer is the sole source of drinking water for communities on Cape Cod, Massachusetts, including the Town of Falmouth, where the aquifer is currently threatened by contamination from septic-system-derived nitrogen. To address this problem, the Town is proposing to sewer areas of Falmouth, treat the wastewater at the Town’s Main Wastewater Treatment Facility (a nitrogen removing/tertiary treatment facility), and discharge the treated wastewater to an ocean outfall pipe in Nantucket Sound. The U.S. Geological Survey, in cooperation with the Town of Falmouth, updated a three-dimensional steady-state groundwater flow model to represent current (defined as 2019–23) average hydrologic conditions and to simulate the long-term average freshwater hydrologic response to two wastewater-return-flow scenarios. Scenario 1 involves the sewering of all properties south of Route 28 in Falmouth, which approximates the Town’s possible sewer expansion over the next 20–30 years. Scenario 2 involves sewering of all properties in Falmouth to demonstrate the maximum potential effect of sewering on the aquifer. Overall, the simulated hydrologic response of water-table altitudes and streamflow in both scenarios was relatively small compared to fluctuations from natural recharge. In scenario 1, the water-table altitude decreased by about 0.1 feet south of Route 28, where the conversion to municipal sewers removed wastewater-return flow from onsite septic systems. The water-table altitude decreased by about 0.1–0.2 feet over a larger area in Falmouth under town-wide sewering in scenario 2. The greatest decrease in water-table altitude in both scenarios occurred near the Main Wastewater Treatment Facility, with a decrease of about 1.1 feet in scenario 1 and about 1.3 feet in scenario 2. Simulated decreases in streamflow also were estimated for six selected streams in Falmouth and Mashpee. In both scenarios, the largest simulated decreases in streamflow were at the Coonamessett River, which is the closest stream to the Main Wastewater Treatment Facility. In scenario 1, the average annual decrease in flow at the Coonamessett River was 0.1 cubic feet per second, a 1.1 percent decrease from current (2019–23) conditions. In scenario 2, streamflow at the Coonamessett River decreased by 0.6 cubic feet per second, a 5.4 percent decrease from current (2019–23) conditions.

Massachusetts

GST-1: A high-resolution global sediment thickness model

Global Sediment Thickness 1 (GST-1) is a high-resolution sedimentary thickness model calculated on a 0.125° x 0.125° grid. It modifies the sediment thickness of the 1° x 1° Earth Crustal Model 1 (ECM1) by means of 3D inversions of free air gravity anomalies. GST-1 is calculated by performing structural inversions on high-density contrasts across two crustal boundaries: the sediment – basement interface and the crystalline crust – upper mantle interface. The inversions are calculated in each of ten overlapping 3D models that span the globe. These ten models are merged to obtain the GST-1 global model, providing an eight-fold increase in lateral spatial resolution in comparison with ECM1 and CRUST 1.0. Our sediment thickness model exploits the nearly continuous sampling of gravity data when compared to the irregular, sparse sampling of seismic refraction data. Sediment thickness values in GST-1 are in excellent agreement with independently derived cross sections from well-studied sedimentary basins, and within expected resolution limits of seismic refraction data. GST-1 offers a robust, high resolution global model of sedimentary thickness to support studies of sedimentary basins.

Tectonophysics

New constraints on northeast Seattle basin structure from converted seismic waves

The Seattle basin is a deep sedimentary basin in the Seattle–Bellevue, Washington metropolitan area within the Puget Lowland of Washington State. We determine the structure of a portion of the basin and the underlying basement using analysis of P waves converted from direct S incident from below. A deep local crustal event beneath Monroe, about 35 km northeast of Seattle, was recorded by a 100‐station nodal array deployed in 2019. The event produced a variety of coherent seismic phases, including converted waves from the sediment—basement boundary, internal structure within the basin, and additional crustal discontinuities. Using observed Sp converted waves, we apply an adjoint‐based full waveform inversion (FWI) method to determine the amplitude and extent of seismic discontinuities at depth. We find the strongest source of converted waves for this event lies ∼6 to 7 km depth below northern Lake Washington, interpreted to be the local depth to basement rock. The newly imaged shallow basement structure may be part of a deformation zone associated with the Siletzia eastern boundary. Our results highlight the utility of converted seismic waves recorded by a dense array, combined with an FWI method, to illuminate crustal structure.

Washington

A 481 m-high landslide-tsunami in a cruise ship-frequented Alaska fjord

Early in the morning of 10 August 2025, a >64 × 10 6 –cubic meter landslide struck Tracy Arm fjord in Alaska. The landslide was preconditioned by glacial retreat caused by climate change. The resulting 481-meter runup megatsunami followed an initial 100-meter-high breaking wave traveling at >70 meters per second. The landslide was preceded by several days of microseismicity, which increased in rate and magnitude until ~1 hour before failure. The landslide produced globally observed long-period seismic waves equivalent in size to a moment magnitude 5.4 earthquake. A long-period (~66 second) global seismic signal, produced by a landslide-induced seiche trapped within the fjord, persisted for up to 36 hours, the second time a days-long seiche had thus been observed. With fjord regions increasingly visited by cruise ships, and climate change making similar events more likely, this unanticipated, near-miss event highlights the growing risk from landslides and tsunamis in coastal environments.

Alaska

Uncertainty in ground-motion-to-intensity conversions significantly affects earthquake early warning alert regions

We examine how the choice of ground‐motion‐to‐intensity conversion equations (GMICEs) in earthquake early warning (EEW) systems affects resulting alert regions. We find that existing GMICEs can underestimate observed shaking at short rupture distances or overestimate the extent of low‐intensity shaking. Updated GMICEs that remove these biases would improve the accuracy of alert regions for the ShakeAlert EEW system for the West Coast of the United States. ShakeAlert uses ground‐motion prediction equations (GMPEs), which calculate spatial distributions of peak ground acceleration (PGA) and peak ground velocity (PGV) from earthquake source estimates, combined with GMICEs to translate GMPE output into modified Mercalli intensity (MMI). We find significant epistemic uncertainty in alert distances; near‐source MMI estimates from different GMICEs can differ by over 1 MMI unit, and MMI extents used for public EEW alerts can differ by hundreds of kilometers for larger magnitude earthquakes ( M ∼6.5+). We use a catalog of “Did You Feel It?” shaking reports to evaluate how well GMICEs predict observed shaking. Our preferred GMICE is the one that computes MMI using PGV for high intensities and transitions to using PGA for nondamaging intensities. These results motivate updating GMICE relationships more generally, including in ShakeMap applications.

The Seismic Record

Importance of fish in the diet of Adélie penguins across multiple life stages

Over recent decades, the Adélie penguin ( Pygoscelis adeliae ) population has grown across most of its range, the exception being the northern coast of the western Antarctic Peninsula (WAP). In the Ross Sea, the very large Cape Crozier colony grew to reach ~9% of the global population. Demographic factors driving the Cape Crozier trend have yet to be identified. However, low chick fledging mass - a negative influence in the northern WAP where the diet is mostly less-energy-dense krill - is not showing an effect in the Cape Crozier colony. Previously, we hypothesized that Cape Crozier fledglings, and post-breeding adults, must be finding sufficient higher-quality prey once free of the prey-depleted colony foraging area, with subsequent higher survival. Here we explore penguin diet within and outside the colony foraging area using stable isotope analysis (SIA) of feathers: those grown by post-breeding adults outside the foraging area, prior to moult, compared to those of near-fledged chicks raised on meals obtained within the foraging area. Second, we explore whether a major difference occurs between breeding and post-breeding diet of adults from a small, nearby colony (Cape Royds) that exploits a smaller, little-depleted preyscape. Finally, we compare the results to SIA analyses published by others to explore whether the pre-moult diet of Ross Island adults, mostly foraging within the north-eastern Ross Sea, compares with that of adults nesting in the southern WAP (Bellingshausen Sea coast). At the latter, populations are not decreasing, with diet containing an appreciable contribution of fish. The results confirmed that fish contributed greatly to chick diets at both Ross Island colonies (δ 15 N 10.6–11.4). Once adults were released from central-place foraging, SIA levels (δ 15 N 9.7–12.9) were similar to those of the southern WAP. Foraging on energy-dense fish in coastal Antarctic waters could be one reason Adélie penguins are among the most abundant and therefore most ecologically successful penguin species.

Antarctic Science

Observations of tear-drinking by lepidopterans on moose ( Alces alces americana ) in northeastern North America

Lepidoptera have long been known to feed on the tears of vertebrates as a presumed source of minerals or nutrients. While this unusual behavior has been observed in a variety of species, only a single previous record has been documented outside of the tropics. Here, we present the first documentation of moths visiting the eyes of a bull moose ( Alces americanus americanus ), captured via trail camera in Green Mountain National Forest, Vermont, United States. We discuss the biogeography of this behavior, how it may differ between tropical and temperate climates, and its potential impact on moose health.

Vermont

Heat-flow data from southeastern Oregon

With the exception of values from two holes drilled within 2 km of Mickey Hot Springs, 17 new heat-flow values in southeastern Oregon are within or somewhat below the range one would normally expect in non-anomalous parts of the North American Cordillera. This is not surprising for a region in which most igneous rocks on the surface are 5 m.y. old or more. There is a suggestion of a thermal anomaly associated with the very young (late Pleistocene or Holocene) Diamond Craters lava field, and the thermal regime on both sides of Steens Mountain seems to be controlled, to some degree, by lateral and vertical movement of water.

Oregon

Evidence for nonlocal sediment transport on hillslopes from fault scarp morphology

Hillslope sediment transport processes such as bioturbation, rainsplash, and granular mechanics occur across the entire planet. Yet, it remains uncertain how these small-scale processes act together to shape landscapes. Longstanding hillslope diffusion theory posits that hillslope processes are spatially limited, whereas new concepts of nonlocal sediment transport argue otherwise. However, each theory produces subtly different, but distinct, predictions for the evolution of fault scarps. We use the topographic change of fault scarps to demonstrate that hillslope processes produce nonlocal sediment transport. Analysis of a global compilation of 340 dated single-earthquake scarp profiles reveals a statistically significant ( p < 0.05 ) relationship between scarp age and scarp asymmetry, here defined as the ratio of imaginary to real components of the Fourier transform of absolute slope. Numerical simulations show that nonlocal models predict this relationship, whereas hillslope diffusion models do not. To further investigate this result, we examined the depositional geometry of a well-exposed colluvial wedge along the Wasatch fault in central Utah, United States. Our quantitative comparison between the exposure and numerical simulations reveals better agreement with the nonlocal model. Nonlocal sediment transport theory appears to best capture the physics of how hillslope processes shape fault scarps, yet hillslope diffusion provides a useful approximation in many cases. As the processes that act on fault scarps are nearly identical to those acting on hillslopes, our results provide evidence supporting nonlocality as a generalized model of hillslope sediment transport.

Utah

Re-Os geochronology of molybdenite by LA-ICP-MS/MS

Re-Os geochronology is a powerful tool to directly date sulphide mineralisation using phases such as molybdenite, pyrite and chalcopyrite. We present a snapshot of a method in development utilizing laser ablation inductively coupled-plasma ‘triple quad’ mass spectrometry (LA-ICP-MS/MS) for Re-Os geochronology of molybdenite. While traditional LA-(MC)-ICP-MS measurements are limited by the >90% isobaric interference correction of parent radioactive 187 Re on radiogenic daughter 187 Os, the MS/MS technology minimizes this correction through careful instrument tuning and use of a (CH 4 + H 2 + He) reaction gas mixture to mass shift 187 Os away from 187 Re. We demonstrate the potential of the MS/MS technique by presenting robust Re-Os dates for three different molybdenite samples ranging in age from ~28 to 1520 Ma.

Conference Paper

A model uncertainty quantification protocol for evaluating the value of observation data

The history-matching approach to parameter estimation with models enables a powerful offshoot analysis of data worth—using the uncertainty of a model forecast as a metric for the worth of data. Adding observation data will either have no impact on forecast uncertainty or will reduce it. Removing existing data will either have no impact on forecast uncertainty or will increase it. The history-matching framework makes it possible to perform this quantitative analysis leveraging the connections among observations, model parameters, and model forecasts. We show this behavior on a specific groundwater flow model of the Mississippi Alluvial Plain and show where the analysis can be informative for considering the potential design of an observation network based on existing or potential observations.

Scientific Investigations Report

Elastic stress coupling between supraglacial lakes

Supraglacial lakes have been observed to drain within hours of each other, leading to the hypothesis that stress transmission following one drainage may be sufficient to induce hydro-fracture-driven drainages of other nearby lakes. However, available observations characterizing drainage-induced stress perturbations have been insufficient to evaluate this hypothesis. Here, we use ice-sheet surface-displacement observations from a dense global positioning system array deployed in the Greenland Ice Sheet ablation zone to investigate elastic stress transmission between three neighboring supraglacial lake basins. We find that drainage of a central lake can place neighboring basins in either tensional or compressional stress relative to their hydro-fracture scarp orientations, either promoting or inhibiting hydro-fracture initiation beneath those lakes. For two lakes located within our array that drain close in time, we identify tensional surface stresses caused by ice-sheet uplift due to basal-cavity opening as the physical explanation for these lakes' temporally clustered hydro-fracture-driven drainages and frequent triggering behavior. However, lake-drainage-induced stresses in the up-flowline direction remain low beyond the margins of the drained lakes. This short stress-coupling length scale is consistent with idealized lake-drainage scenarios for a range of lake volumes and ice-sheet thicknesses. Thus, on elastic timescales, our observations and idealized-model results support a stress-transmission hypothesis for inducing hydro-fracture-driven drainage of lakes located within the region of basal cavity opening produced by the initial drainage, but refute this hypothesis for distal lakes.

JGR Earth Surface