USGS ScienceSearch

SEARCH · USGS Science

Results for “Earth Surface Processes and Landforms”

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,665 records · Page 91Linked to original sources

Horizontal transport of Picture Gorge Basalt magma through the Monument Dike Swarm determined by magnetic fabric

Flood basalts of the mid-Miocene Columbia River Basalt Group (CRBG) cover 210,000 km 2 of Washington, Oregon, and Idaho. The source of CRBG melt is debated; widely spaced feeder dike swarms can be projected toward hypothetical sources near the Oregon-Idaho border. In this study, we use anisotropy of magnetic susceptibility (AMS) to track magma flow in the Monument dike swarm (MDS), the feeder dikes of the Picture Gorge Basalt (PGB). This small formation of the main-phase CRBG eruptions allows us to explore in detail the localized dynamics of a large igneous province feeder system, with implications for the larger CRBG picture. We measured the magnetic fabric of 205 oriented paleomagnetic specimens subsampled from 97 samples collected from 15 dikes of the MDS. Thermal demagnetization and hysteresis loops show that the magnetic minerals are a mixture of single domain and multidomain sized titanomagnetites. At three dikes, the paleodepth of sampling was determined to be shallow (<350 m). Magma flowing through dikes has been shown—in most cases— to acquire an anisotropic magnetic fabric with an AMS ellipsoid minimum axis perpendicular to the wall and maximum axis aligned in the direction of flow. Of 15 dikes, 12 show horizontal flow directions in the plane of the dike. Only one dike displayed imbricated fabrics, showing westward flow away from the Oregon-Idaho border. We conclude that magma flow in the MDS was sub-horizontal from a distal source.

Oregon

Variable pH habitats could help prepare stone crabs for coastal acidification

Coastal acidification is being exacerbated by terrestrial organic inputs, especially after high precipitation events. Florida’s rainy season coincides with stone crab ( Menippe mercenaria) reproduction, and pH extremes could limit future harvests by reducing reproductive output. Populations that experience pH variability can serve as “natural laboratories” for estimating the tolerance of a species to coastal acidification. Here, we conducted a series of experiments to determine if ovigerous stone crabs conditioned in more variable pH habitats (seagrass) would result in faster embryonic development, greater hatching success, and higher larval survival relative to crabs conditioned in habitats with lower pH variability (sandy habitats). After field conditioning, crabs were transported to the laboratory and randomly acclimated to either a control pH (pH = ~ 7.90) or a reduced pH condition (pH = ~ 7.60) until larval release. The rate of embryo development was slower in the laboratory reduced pH treatment, however, there were no observable field condition effects on embryo development rate. Crabs conditioned in the more pH variable seagrass habitat did have greater hatching success and higher larval survival than crabs in less pH variable sandy habitats; however, larval survival was low across all treatments. These results suggest that the pH variability experienced in seagrass habitats during brooding may serve as a mechanism for stone crabs to acclimatize to extremes in seawater pH.

Marine Biology

Annual NLCD (National Land Cover Database)—The next generation of land cover mapping

Introduction The widely used National Land Cover Database (NLCD) has long been the foundational land cover source for scientists, resource managers, and decision makers across the United States. In 2024, a reinvention as Annual NLCD added the key improvement of annual time steps to show decades of change at a higher frequency than the intervals of 2–3 years used in the legacy NLCD. Annual NLCD was derived primarily from the long Landsat satellite data record, and it includes data from other sources. The first release in 2024 of Annual NLCD provides Collection 1.0 of products encompassing land cover and land change from 1985 through 2023 for the conterminous United States (CONUS). The Annual NLCD Collection 1.0 consists of six operational products that map the unique characteristics of land cover. A map created from Annual NLCD shows 16 land cover classes for the CONUS in 2023.

conterminous United States

Predictions of groundwater PFAS occurrence at drinking water supply depths in the United States

Per- and polyfluoroalkyl substances (PFAS), known colloquially as “forever chemicals”, have been associated with adverse human health effects and have contaminated drinking water supplies across the United States owing to their long-term and widespread use. People in the United States may unknowingly be drinking water that contains PFAS because of a lack of systematic analysis, particularly in domestic water supplies. We present an extreme gradient boosting model for predicting the occurrence of PFAS in groundwater at the depths of drinking water supply for the conterminous United States. Our model results indicate that 71 to 95 million people in the conterminous United States potentially rely on groundwater with detectable concentrations of PFAS for their drinking-water supplies prior to any treatment.

conterminous United States

Deep syntectonic burial of the Anthracite belt, Eastern Pennsylvania

Fluid inclusion microthermometry and Raman spectroscopy of fluid inclusions in quartz veins from the Pennsylvanian rocks of the Anthracite belt, eastern Pennsylvania support a deep burial model of coalification in favor of focused orogenic hot fluid flow. High-temperature (250 to 255 °C) trapping of CH 4 ± CO 2 saturated aqueous fluids and CH 4 ± CO 2 inclusions indicate fluid trapping at depths of 11.5 to 13.4 km under a cover of Pennsylvanian to Permian(?) syntectonic load. In the folded rocks to the south of the Anthracite belt, CH 4 ± CO 2 fluid inclusions indicate a sediment load that was up to 16.3 km thick. Re -equilibrated aqueous fluid inclusions from veins in Silurian through Devonian rocks give the same range of trapping conditions but a wide range of fluid salinities suggesting that folding, fracturing, and meteoric recharge resulted in the intermixing of fluids from throughout the stratigraphic succession.

Pennsylvania

topoBuilder quick start guide

TopoBuilder is a public web application from the National Geospatial Program that enables anyone to create customized digital U.S. Geological Survey (USGS) topographic maps, called OnDemand Topos, with the best available, most up-to-date data from The National Map (nationalmap.gov). OnDemand Topos can be made at different scales or quadrangles and can cover anywhere within the United States and its territories.

Fact Sheet

Exploratory water well, St. George Island, Florida

A 1,026-foot well was drilled into the Floridan aquifer on St. George Island in search of a supply of water for residential development in 1971. This is the only well on the island on record to have been drilled below 200 feet for water. The well was abandoned; water flowing from the well was saline. It had a specific conductance of 7,500 micromhos per centimeter at 25° Celsius. The specific conductances of samples taken at drilling depths 217-1,026 feet ranged from 3,140 to 39,600, with the maximum value at 470 ft. A composite geophysical and lithologic log of the well shows spontaneous potential, single-point resistance, gamma-ray, neutron, gamma-gamma, and caliper logs. Well-site notes, a written description of cuttings samples, a graph showing the relation of specific conductance of water samples to the depth of the well at time of sampling, and four chemical analyses of the water are included in the report.

Florida

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

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

Modeling future groundwater depletion to evaluate sustainability goals set under the Sustainable Groundwater Management Act in the critically overdrafted basins of the Central Valley, California, USA (2020–2070)

In 2014, California's Sustainable Groundwater Management Act (SGMA) mandated local agencies to devise and implement groundwater sustainability plans to address critically overdrafted conditions throughout the state's aquifers. However, the feasibility of these agencies' sustainability goals has not previously been assessed through a regional-scale, integrative lens. Here, we develop and analyze a novel, basin-wide database of 936 sustainability indicator wells located within Central Valley subbasins designated as critically overdrafted, most of which lie in the San Joaquin Valley. Our database shows 2040 groundwater elevation goals vary widely from 60 m above to 80 m below 2020 levels, with variability within and between adjacent subbasins. To evaluate the feasibility of achieving these goals, we coupled the database with a regional hydrologic model (Central Valley Hydrologic Model version 2) and simulated multiple future pumping scenarios. Results show that under increased groundwater demand, 60%–70% of indicator wells may fail to meet their 2040 goals. Even a 50% reduction from 2020 demand levels leaves nearly 40% of wells failing to meet their sustainability thresholds by 2040. Baseline models show that by 2070, up to 70% of wells could fail to meet their goals due to large-scale, spatially connected regions of groundwater depletion. This integrated framework, linking the first region-wide compilation of SGMA indicator wells with a regional groundwater model, demonstrates that many local sustainability goals may be unattainable with substantial (up to 50%) reductions in pumping. Additional management interventions, such as expanded recharge or coordinated demand reductions, may help achieve sustainability goals.

California

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

The U.S. Geological Survey National Water Quality Network—Groundwater—2023

The U.S. Geological Survey (USGS) operates a National Water Quality Network (NWQN) to monitor trends in groundwater quality and assess emerging contaminants of concern. It is a “network of networks” with 81 subnetworks being sampled on a decadal time scale. Each year, eight of the subnetworks are sampled. Subnetworks have 20–30 wells each and include studies of domestic supply wells or shallow groundwater (20–50 feet deep) underlying urban land use or agricultural land use. Currently there are 2,114 wells in the network.

conterminous United States

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