USGS Science⌕ Search

SEARCH · USGS Science

Results for “Environmental Monitoring and Assessment”

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,662 records · Page 77Linked to original sources
USGS2025

Near-surface material and topography generate anomalous high-frequency ground motion amplification in Chugiak, Alaska

An ∼3 km long nodal array oriented approximately east–west was deployed in Chugiak, Alaska, by the U.S. Geological Survey during 2021. The array intersects with the permanent NetQuakes station NP.ARTY, where peak ground acceleration (PGA) value of 1.98 g was recorded during the 2018 M w 7.1 Anchorage, Alaska, earthquake, in sharp contrast to the PGA of ∼0.3 g at a site just 4 km to the west. Seismic data for M w 1.8–4.3 aftershocks from the M w 7.1 event recorded by the nodal array confirm the anomalously large ground motions obtained at NP.ARTY as well as similar amplifications at nodes within ∼1 km to the east. Here, we performed 0–10 Hz 3D finite‐difference simulations, including high‐resolution surface topography, to explore the cause of the unexpectedly large amplification. As expected, the simulations computed with a regional 3D tomography velocity model severely underpredict the 0–10 Hz acceleration records at almost all sites. Adding a near‐surface low‐velocity taper to 300 m depth amplifies the accelerations by up to a factor of 5 and enables a reasonable match between the nodal data and simulations at sites to the west of NP.ARTY. However, this model still underpredicts the spectral energy in the area covered by glacial sediments by up to an order of magnitude. The addition of a till layer using a depth‐dependent shear‐wave velocity ( ⁠⁠ V s ) profile along with a homogeneous, 8 m thick low‐velocity layer with V s = 250 m/s representing the kame terraces improves the fit to data to within a factor of 2 at nodes located on top of the glacial sediments. Our study shows that the anomalously large high‐frequency amplification recorded at and near NP.ARTY can be explained by a combination of topographic effects and near‐surface low‐velocity material with amplification effects on the high‐frequency ground motion by up to about 40% and an order of magnitude, respectively.

Alaska↗
USGS2024

Effects of stochastically-simulated near-fault ground motions on soil liquefaction

The scarcity of historically recorded near-fault ground motions poses a challenge to systematically understanding the influence of near-fault effects on various types of seismic demands for engineering purposes. In particular, the current state of knowledge of the influence of ground-shaking intensity on soil liquefaction and its consequences does not specifically account for the effects of near-fault ground motion characteristics. In this study, the influence of near-fault ground motions on liquefaction triggering and lateral spreading are investigated using non-linear modeling of a hypothetical liquefiable soil column in the finite-element computational platform OpenSees subjected to simulated ground motion time series that represent strong earthquake shaking in the near field. The simulated ground motion time series and resulting datasets are based on a parametric stochastic model and are developed for a range of source and path parameters to represent a realistic variability of ground motion characteristics. Dependencies between ground motion intensity measures (IMs) and liquefaction demand parameters are investigated for near-fault pulse and nonpulse-like ground motion sets. Evolutionary IMs, such as cumulative absolute velocity (CAV) and the time-varying magnitude-adjusted peak ground acceleration (PGAM), are considered in developing liquefaction triggering probability density functions. Post-liquefaction triggering responses such as lateral spreading displacements are examined in relation to PGA M and CAV. The ground motion simulations are validated by comparing their liquefaction-capacity PGA M fragilities and post-triggering CAV vulnerability relationships to historical records from the 1994 Northridge earthquake in California, USA. Finally, a path forward for future studies that includes finding systematic differences in the IM-liquefaction demand relationships between near-fault and far-field stochastic ground motion sets is outlined.

Conference Paper↗
USGS2000

Regional evaluation of evapotranspiration in the Everglades

Nine sites in the Florida Everglades were selected and instrumented for collection of data necessary for evapotranspiration-determination using the Bowen-ratio energy-budget method. The sites were selected to represent the sawgrass or cattail marshes, wet prairie, and open-water areas that constitute most of the natural Everglades system. At each site, measurements necessary for evapotranspiration (ET) calculation and modeling were automatically made and stored on-site at 15- or 30-minute intervals. Data collected included air temperature and humidity at two heights, wind speed and direction, incoming solar radiation, net solar radiation, water level and temperature, soil moisture content, soil temperature, soil heat flux, and rainfall. Data summarized in this report were collected from January 1996 through December 1997, and the development of site-specific and regional models of ET for this period is described. Latent heat flux is the energy flux density equivalent of the ET rate. Modified Priestley-Taylor models of latent heat flux as a function of selected independent variables were developed at each site. These models were used to fill in periods of missing latent heat flux measurement, and to develop regional models of the entire Everglades region. The regional models may be used to estimate ET in wet prairie, sawgrass or cattail marsh, and open-water portions of the natural Everglades system. The models are not applicable to forested areas or to the brackish areas adjacent to Florida Bay. Two types of regional models were developed. One type of model uses measurements of available energy at a site, together with incoming solar energy and water depth, to estimate hourly ET. This available-energy model requires site data for net radiation, water heat storage, and soil heat flux, as well as data for incoming solar radiation and water depth. The other type of model requires only incoming solar energy, air temperature, and water depth data to provide estimates of hourly ET. The second model thus uses data that are more readily available than the data required for the available-energy model. Computed ET mean annual totals for all nine sites for the 1996-97 period ranged from 42.4 inches per year at a site where the water level is below land surface for several months each year to 57.4 inches per year at an open-water site with no emergent vegetation. Although the density of photosynthetically-active plant leaves has been shown to relate directly to ET in some studies, it does not appear to relate directly to ET in the Everglades, based on comparison of annual ET data with leaf-area index, defined as the Normalized Difference Vegetation Index (NDVI), data from satellite imagery. NDVI and ET appear to be inversely related in the Everglades. The greatest ET rates occurred at open-water sites where the NDVI data indicated the lowest leaf-area index. Among the remaining vegetated sites, there is no clear relation between ET and NDVI, though the highest ET rate corresponded to the lowest NDVI and one of the lowest ET rates corresponded to the highest NDVI value. The variation in ET follows a seasonal pattern, with lowest monthly ET totals occurring in December through February, and highest ET occurring in May through August. The greatest range in monthly ET among all nine sites for the 2-year period occurred at site 3: from 1.81 inches in December 1997 to 6.84 inches in July 1996. A study to compare the Bowen-ratio/energy balance method of ET measurement with the eddy-correlation method was done at one site from June 22, 1998, through September 28, 1998. This comparison indicated that both methods gave comparable values of the Bowen ratio, but there was a considerable difference in available energy measured by the two methods. The mean of all 30-minute measured turbulent heat fluxes from the eddy-correlation apparatus for June 22 through September 29, 1998, was 137.4 watts per square meter, and the mean of the corresponding measured energy was 163.6 watts per square meter, or about 20 percent greater. The disagreement in mean energy fluxes measured by the two methods is problematical and is not fully understood. Although the difference seems to be related to friction velocity, and is practically non-existent at values of friction velocity greater than 0.3 meter per second, the "correctness" of either method cannot be determined with the data available.

Florida↗
USGS2025

A comprehensive geologic framework of the National Crustal Model for seismic hazard studies in the conterminous United States

A three-dimensional (3D) geologic framework has been developed for the conterminous United States (U.S.) as part of the U.S. Geological Survey National Crustal Model to enhance seismic hazard modeling. The geologic framework is created from geologic maps and multiple subsurface geologic unit boundaries including the base of the Miocene, Cenozoic, Phanerozoic, and the Mohorovičić discontinuity. Modifications are made to surficial geologic maps to remove discontinuities across state and country borders. The subsurface distribution of rock type and age is extrapolated from the surface, seeded with subsurface geologic information, and constrained by a map of basement geology. The framework provides the basis for estimates of subsurface seismic velocity and density that is needed to improve estimates of earthquake ground shaking and seismic hazard. The present framework greatly expands and updates a previously published 3D geologic framework of the western part of the U.S. that was itself a first-of-its-kind digital 3D portrayal of the nation.

conterminous United States↗
USGS2024

California State Waters Map Series—Benthic habitat characterization in the region offshore Humboldt Bay, California

Coastal and Marine Ecological Classification Standard (CMECS) geoform, substrate, and biotic component geographic information system (GIS) products were developed for the California State Waters of northern California in the region offshore of Humboldt Bay. The study was motivated by interest in development of offshore wind-energy capacity and infrastructure in Federal waters offshore. This project, carried out by the U.S. Geological Survey (USGS), resulted in four data releases for individual map blocks that are part of the “California State Waters Map Series”: (1) Offshore of Arcata, (2) Offshore of Eureka, (3) Offshore of the Eel River, and (4) Offshore of Cape Mendocino. The study area consists of 436 square kilometers of multibeam echo sounder (MBES) data acquired by Fugro Pelagos, Inc., in 2007. Towed camera-sled video was acquired in 2009 and 2010 to supervise the classification of the MBES data into habitats, and single channel sparker data were collected to calculate sediment thickness above the transgressive unconformity. Using video observations of habitat as ground truth, derivatives of the MBES data were classified into 3 seafloor character types (hard-rugose, hard-flat, and soft-flat), 26 induration-slope-depth groups, and 15 geoforms. The study area substrate is predominantly soft-flat sediment (mud and fine sand) covering 73.6 percent of the area. Hard-flat substrate areas, predominantly coarse sediment in scour depressions, cover 5.4 percent of the study area. The hard-rugose substrate areas are primarily outcrops of layered sedimentary bedrock and constitute 20.9 percent of the study area. Fifteen geoforms were identified in the analysis. The predominant geoforms mirror the seafloor character results, shelf geoforms, rock outcrop geoforms, and scour depression geoforms. Rock and scour areas are restricted to the southern portion of the study area off Cape Mendocino where uplift has exposed bedrock. On the flat shelf area post-transgressive sediment varies in thickness from 1.7 meters (m) nearshore to 28.1 m offshore.

California↗
USGS2025

An evaluation of the effects of different deicing salt application rates on three watersheds in Essex County, New York

The U.S. Geological Survey, in cooperation with the New York State Department of Transportation, evaluated the effects of different deicing salt application rates on surface water, groundwater, and highway runoff quality near State highways in northern New York. Three reaches of State highways were tested with different deicing treatments between October 2019 and November 2022: a salt-sand mixture (Treatment A), a salt mixture applied at a lower rate (Treatment B), and a control mixture consistent with typical deicing salt amounts and application rates. Data on pavement conditions and the quality of surface water, highway runoff, and groundwater were collected. Surface electromagnetic data were also collected. Surface-water and groundwater quality downgradient from the State highways were compared with water quality at upgradient locations. The percentage of snow or ice coverage was used to evaluate the effectiveness of the salt applications. This report provides an overview of the transport of deicing salt. The Treatment B watershed had deicing mixture applied more frequently than other highway reaches, which caused it to have the highest annual total chloride application. Despite differences in chloride application, flow-weighted mean chloride concentrations in highway runoff were comparable across treatments. Chloride concentrations were elevated in surface water and groundwater downgradient from highways relative to chloride concentrations upgradient from highways. A chloride mass balance, calculated for one treatment watershed, indicated that groundwater affected by legacy deicing practices may be contributing additional chloride to surface water. Spatial patterns from electromagnetic surveys show a shallow saline plume alongside the highway in that area. Differences in winter severity and pavement-surface conditions drove deicing salt applications in the treatment areas. This study found that several factors affect chloride loads in the watersheds, including variable winter conditions, adaptive snow and ice management, legacy management practices, and area-specific aquifer and groundwater conditions.

New York↗
USGS2024

Bedrock geologic map of the Woodstock quadrangle, Grafton County, New Hampshire

The bedrock geology of the Woodstock 7.5-minute quadrangle consists of highly deformed metasedimentary rocks of the Central Maine trough, including the Silurian Rangeley and Perry Mountain Formations and the Devonian Littleton Formation. The central, northern, and eastern parts of the quadrangle are underlain by the oldest rocks in the area, the Rangeley Formation. In the southwest and south-central part of the quadrangle, metaturbidites of the Perry Mountain Formation and subsequent Littleton Formation overly the Rangeley Formation in a deformed F 1 synform, herein informally called the Bagley Brook basin. The metasedimentary rocks were intruded by widespread syn- to post-tectonic granitoids of the Devonian New Hampshire Plutonic Suite and minor post-metamorphic Jurassic-Cretaceous mafic dikes of the White Mountain Plutonic-Volcanic Suite. The metasedimentary rocks were affected by at least two episodes of deformation in the Devonian Acadian orogeny. The dominant regional foliation is second-generation (S 2 /D 2 ) and formed during the development of sillimanite-muscovite mineral assemblages. Large bodies of the Early Devonian Kinsman Granodiorite intruded the metasedimentary rocks semi-concordantly during D 2 deformation. Dikes of the Late Devonian Concord Granite cut the Kinsman Granodiorite and the metasedimentary rocks and were emplaced either syn- or post-D 2 . The map pattern in the Rangeley Formation is dominated by northeast to northwest trending, moderately to steeply north-dipping F 2 and F 3 folds. Map-scale F 1 folds are defined by the Bagley Brook basin. Previous division of Rangeley Formation stratigraphy in this region into “upper” and “lower” parts was not corroborated by 1:24,000-scale mapping of lithodemic units, and rocks previously mapped as part of the Smalls Falls and Madrid Formations are here reassigned to the Rangeley Formation. Some rocks previously mapped as the lower part of the Littleton Formation are now assigned to the Perry Mountain Formation. The Littleton Formation on this map is approximately equivalent to rocks previously mapped as the upper part of the same formation. Steeply dipping fractures in the quadrangle show a preferred northeast orientation, consistent with subsurface fracture orientations in the well fields near Mirror Lake. Jurassic-Cretaceous mafic dikes and normal faults show preferred northeast orientations, similar to the fractures, suggesting that the extensional stress field that controlled dike orientation during the Mesozoic also produced the dominant brittle fabrics in the area.

New Hampshire↗
USGS2026

Origin of the Pd/Pt ratio of the J-M Reef, Stillwater Complex, Montana, USA

The J-M reef of the Stillwater Complex is characterized by a high Pd/Pt ratio (mean ~3.8 with a standard error of 0.03) with a homogeneous geospatial distribution at the deposit scale. In this contribution, we demonstrate that the Pd/Pt ratio of the reef is the product of equilibration of an immiscible sulfide liquid with a silicate melt rich in Pd relative to Pt. Despite the high tenors of the J-M reef sulfides (avg 2,700 ppm Pt and 770 ppm Pt), numerical modeling shows that the parental melts did not have extraordinary Pd and Pt concentrations. Instead, the initial composition of a plausible parental silicate melt can have Pd and Pt contents well within the expected range of a normal, mantle-derived partial melt (i.e., ~10–20 ppb for both Pd and Pt with Pd/Pt of ~1). The relative differences in the partitioning behavior of Pt and Pd between sulfide liquid and silicate melt are unlikely to produce a consistent Pd/Pt ratio across a wide range of silicate melt to sulfide liquid mass ratios (i.e., R factors). Instead, the pre-emplacement fractionation of Pt alloy from S-undersaturated silicate magma accounts for the homogeneous and high Pd/Pt ratio of the J-M reef. We show that batch equilibration of sulfide liquid with silicate melt can produce the high Pd/Pt ratios of the reef if the partition coefficients between sulfide liquid and silicate melt for Pd and Pt are extremely high (>10 6 ). In an alternative model, Pd enrichment could be achieved by sulfide upgrading in resident footwall mush even if the partition coefficients between sulfide liquid and silicate melt are relatively small (between 10 4 and 10 6 ) because the instantaneous mass ratio of silicate melt to sulfide liquid is small (R ≈ 100–700), so the partitioning behavior of Pt and Pd has little impact on the composition of sulfide liquid.

Montana↗
USGS2026

Subsurface vertical connectivity shapes solute transport to montane streams: Insights from watershed-scale geophysics-informed modeling

Headwater catchments provide essential water and nutrients to downstream ecosystems. The timing and evolution of their transport is shaped in part by the “invisible” subsurface structure but have been largely unexplored at the watershed scale. This raises a fundamental question: how and to what extent does subsurface structure influence streamflow and solute transport in montane watersheds? Here, we test controls of three-dimensionally resolved subsurface structure on solute transport, which reflects source waters and water residence time. We integrate the geophysically mapped architecture of Coal Creek, a headwater catchment in the Colorado River Basin, into numerical flow-transport models. Our results reveal greater vertical connectivity elongates flow paths and enhances deep groundwater contributions to streams, therefore substantially influencing the timing and magnitude of solute transport. In contrast, watershed-integrated outlet discharge can be predicted without high-resolution subsurface information. These findings demonstrate incorporating subsurface architecture is critical for predicting ecosystem health and contaminant transport.

Colorado↗
USGS2024

Geophysical modeling of a possible blind geothermal system near Battle Mountain, NV

The northeastern portion of the Reese River basin in north-central Nevada is the focus of detailed geophysical and geological studies as part of the INGENIOUS project, which aims to identify new, commercially viable hidden geothermal systems in the Great Basin region of the western U.S. This location, herein referred to as Argenta Rise, occupies a broad (~15km wide) left-step between major range-front fault systems along the northwestern edge of the Shoshone Range and Argenta Rim, with numerous ENE-striking intra-basin faults presumably accommodating sinistral-normal oblique slip across the step-over. Four discrete regions have been identified within the study area that have favorable structural settings for hosting a blind hydrothermal system. However, with no definitive or extensive surface manifestations of an active hydrothermal system (e.g., geysers, steam vents, sinter, etc.), detailed geophysical studies are necessary to resolve subsurface geology and structure, and identify zones of enhanced structural complexity that may promote hydrothermal fluid flow. Hence, we collected high-resolution gravity, MT, and rock property data (density, magnetic susceptibility), and analyzed the recently acquired GeoDAWN aeromagnetic data to characterize potential geothermal resources in this region. Using the new geophysical datasets, we jointly modeled gravity and magnetic data along a series of intersecting 2D profiles that integrated information from recent, local-scale fault mapping. Rock property measurements performed on outcrops and hand samples throughout the study area constrained the models. The MT data were used to construct a 3D resistivity model that highlights the location of inferred alteration and fluids in the subsurface. Combined MT and potential field results reveal which structures may be most important for controlling hydrothermal fluid migration, as well as which geologic units may host hydrothermal fluids. Our gravity derived depth to basement surface coincides well with the base of shallow conductive anomalies, suggesting hydrothermal fluids may be confined to basin fill sediments and volcanics. This work supports our development of 3D geophysical and geologic models that are focused along the western flank of the northern Shoshone Range and aids the process of selecting sites for temperature gradient drilling.

Nevada↗
USGS2024

Global geologic map of Europa

Discovered by Galileo Galilei more than 400 years ago and imaged in detail by the Voyager 2 Galileo spacecraft, Jupiter’s icy moon Europa has been a source of intrigue. A range of science investigations indicate that it contains the key ingredients for habitability, notably energy, chemistry, and liquid water. Europa’s surface is geologically complex and, based on the dearth of impact craters, interpreted to be as young as ~60 Ma. The array of geologic features that characterize the surface include extensive ridged plains, regions of broad disruption termed chaotic terrain, long, quasi-linear ridges that span thousands of kilometers, and bands as much as 60 kilometers wide and that extend 100s of kilometers. These features, along with other geophysical measurements, indicate the presence of a global briny liquid water ocean beneath the ice shell. It was not until the arrival of the Galileo spacecraft in 1995 that the true nature and level of complexity of the surface was revealed. Although image data returned by Galileo provided insight into the structure of a variety of regions, the entire satellite has yet to be observed at a regional scale (less than 250 meters per pixel) and the detailed geologic nature of much of its surface remains a mystery. Establishing the global context of the distribution and timing of Europan geologic units forms a basis to understand regional and local scale processes, serves as a tool for the planning of future missions, and most of all is essential to gaining insight into the potential habitability of this icy world.

Scientific Investigations Map↗
USGS2026

Postglacial ecosystem development of a hydrothermal landscape in Yellowstone National Park

The Yellowstone geo-ecosystem has been the subject of much research, but the ecological history of the Yellowstone Plateau volcanic field and its iconic geyser basins is less known. In this investigation, paleoenvironmental analyses of sediment cores from lakes in Lower Geyser Basin were compared with regional records and paleoclimate model simulations to reconstruct the vegetation, wildfire, limnology, hydrothermal dynamics, and climate drivers since deglaciation, 15,000 to 14,000 years ago. Pollen data from Lower Geyser Basin lakes reveal the strong influence of infertile rhyolitic soils on vegetation history: an initial late-glacial steppe was replaced by lodgepole pine forest from 12.8 to 11.0 ka, with little change in forest composition or cover thereafter despite changing climate. This stability contrasts with the more dynamic vegetation response on nonrhyolite substrates in the Yellowstone region where nutrient and moisture availability is greater. Highest wildfire activity and low lake nutrient levels in Lower Geyser Basin occurred from 12 to 4 ka, when summers were substantially warmer and drier and fire-inducing vapor pressure deficits were 29 to 56% higher. The hydrothermal history, inferred from sedimentary arsenic and cesium abundances, was spatially and temporally variable but lake-forming hydrothermal events align with periods of abundant moisture. Thus, long-term changes in wildfire, limnology, and, to some extent, hydrothermal activity were governed by insolation-driven climate variations, whereas the vegetation response was muted and constrained by geologic processes. These findings suggest that warmer, drier conditions in the future could result in less hydrothermal activity yet little change in forest cover across the Yellowstone Plateau volcanic field despite more wildfires.

Idaho, Montana, Wyoming↗
USGS2025

A method to obtain remotely sensed grain size distributions from nonplanar granular deposits

Constraining the grain size distribution of granular deposits with complex surfaces is difficult with existing approaches. Field and laboratory techniques are time consuming and limited by the maximum grain size that laboratories can accommodate. In this study, we present a new method to identify the coarse fraction of the grain size distribution at a debris-flow fan deposit surveyed with terrestrial laser scanning (TLS) in Glenwood Canyon, Colorado, USA. This method is a novel grain segmentation algorithm developed for application to point cloud data of deposits with complex surfaces and angular grains ranging in size from centimeters to a meter. This approach combines an existing random forest machine learning method with a novel iterative clustering algorithm. We compared the grain size distribution from our algorithm with a Wolman pebble count conducted in the field, and found a root mean squared error of less than 2 cm from the 5th to 95th percentile of the grain size distribution of grains ranging from cobble to boulder sized (6.3–78 cm in our application). Finally, we compared our new algorithm with an existing open-source grain segregation algorithm, and our method outperformed the selected alternative when applied to the debris-flow deposit point cloud.

Colorado↗
USGS2024

Coastal Science Navigator companion guide—Discover the U.S. Geological Survey coastal science products you need

The Coastal Science Navigator is an online gateway to a wide variety of U.S. Geological Survey (USGS) coastal change hazards-related information, data, and tools relevant to stakeholders’ scientific and decision-making needs. The products within the Coastal Science Navigator provide data related to past, present, and future threats to our coastlines. The filter search allows users to see all available products and identify relevant options by time scale, geographic scope, coastal hazard theme, and other filters. The guided search suggests products based on users’ answers to a short series of questions. A comprehensive summary is available for each product. The idea for the Coastal Science Navigator arose in 2020 in response to stakeholder feedback identifying the need for a central source for USGS coastal science information. It was published in July 2023 and initially included 55 products. Regular updates are planned to integrate other existing and new products. This guide introduces some of the many coastal change hazards-related products available through the USGS. In it, we showcase the products included in the Coastal Science Navigator’s initial publication in July 2023. While it is not representative of all the information, tools, and data available, we hope it serves as a compelling snapshot of what the USGS has to offer and encourages you to explore the Coastal Science Navigator to discover more of the products you need. To navigate this guide, the products have been organized by the time scale they are best suited for—past, present, or future—although many products cover multiple time scales. An additional section features software, one of the many product types available as filters within the Coastal Science Navigator. Other products include downloadable data, websites, and geonarratives (web pages that combine text, images, and interactive maps into narratives you can scroll through). Featured geographic scopes are also highlighted within this guide, detailing some of the many regions in which the USGS conducts research and illustrating another way to filter products within the Coastal Science Navigator.

Circular↗
USGS2026

The digital archivist: Automating legacy macroseismic data processing using large language models

Macroseismic data are a key resource to investigate shaking and damage from preinstrumental and early instrumental eras. However, data are often stored as inconsistently formatted reports describing observed shaking and damage, making manually parsing and interpreting accounts labor‐intensive. We introduce a novel workflow using Google’s Gemini 2.5 Pro large language model (LLM) to automate the extraction and structuring of macroseismic observations from summary reports. We apply this workflow to the 22 March 1957 M 5.3 Daly City, California, earthquake as a case study. We used Gemini to extract addresses, originally assigned modified Mercalli intensity values, and descriptions from each report. To address coordinate precision limits, addresses were geocoded via Google’s Geocoding application programming interface. This workflow yielded over 2300 geocoded intensity reports for the Daly City earthquake. We use the geocoded accounts, with the original report intensity assignments, to develop a shaking intensity map that in some respects rivals modern Did You Feel It? Maps. We also extract and present data for the 9 February 1971 M L 6.7 Sylmar, California, earthquake. Our results demonstrate the potential of LLMs for reliably extracting and analyzing large, unstructured macroseismic datasets. LLMs offer a scalable solution for rapidly digitizing macroseismic archives, enabling their broader use to constrain ground‐motion models in modern seismic hazard analysis and to improve our understanding of site effects in urban areas. The concepts explored here may also be applied to the handling of other legacy seismological and earth science data.

Seismological Research Letters↗
USGS2025

Did they feel it? Legacy maroseismic data illuminates an engimatic 20th century earthquake

The challenges and the importance of preserving legacy instrumental records of earthquakes are now well-recognized (e.g., Richards & Hellweg, 2020, https://doi.org/10.1785/0220200053 ). Seismologists may not be aware of parallel challenges and opportunities with legacy macroseismic data for earthquakes in the United States. For much of the 20th century, macroseismic data were collected by a series of U.S. government agencies using a standard questionnaire distributed on postcards. Published summaries of postcards provide macroseismic data akin to modern Did You Feel It? questionnaire responses. In this paper we focus on the M 6.5 Fickle Hill, California earthquake, on 21 December 1954 (Hellweg et al., 2025) as a proof-of-concept, illustrating the potential of what we dub Did They Feel It? (DTFI) data to improve our understanding of significant 20th century U.S. earthquakes for which instrumental data are sparse. Legacy macroseismic data interpreted following modern conventions can potentially constrain traditional ShakeMaps at a level of detail and accuracy that in some respects rival maps for modern earthquakes. The updated ShakeMap for the 1954 Fickle Hill earthquake, also drawing from recently published media and first-person accounts, supports the location, depth, and stress drop value estimated from available instrumental data (Hellweg et al., 2025).

California↗
USGS2025

Widespread anhydrite saturation in Laramide-age arc magmas of southwestern USA

Anhydrite is considered a rare mineral phase in magmas, with only ∼33 documented occurrences worldwide. However, anhydrite readily decomposes in the near-surface environment, making it difficult to recognize its former presence in rocks collected at or near Earth’s surface. In such samples, only small anhydrite inclusions fully shielded within other minerals can have survived. During a recent field trip to the southwestern USA, we sampled 17 Laramide-age (ca. 40−80 Ma) magma systems, most of which are associated with porphyry copper deposits. A systematic search for anhydrite inclusions preserved within apatite, amphibole, plagioclase, and quartz phenocrysts in ∼100 rock samples by optical microscopy and Raman spectroscopy revealed that each of these 17 magma systems was at least temporarily anhydrite-saturated. Also, most previously identified magmatic anhydrite-bearing intrusions are associated with porphyry copper deposits, and both intrusive and volcanic rocks containing magmatic anhydrite show high Sr/Y ratios. These observations suggest that anhydrite saturation and porphyry copper formation are linked via magma fractionation at high pressure. Compared to average arc magmas, anhydrite-bearing magmas are unusually oxidized and sulfur-rich and seem to also be unusually water-rich. Hence, our preferred interpretation is that magma generation and/or fractionation at high pressure promotes the formation of superhydrous and oxidized magmas, which in turn promotes high sulfur contents and ultimately the precipitation of anhydrite. The high mineralization potential of these magmas does not need to result from their high sulfur content but could be caused by other properties of high-pressure magmas.

Arizona, New Mexico↗
USGS2026

Black rails in Baja California

Summarizing records from 1905 to 2025, we document the presence and distribution of the California black rail ( Laterallus jamaicensis coturniculus ) in northwestern México, particularly Baja California and the Colorado River Delta. This subspecies has a restricted distribution in the western U.S. and northwestern México. We found few records over the past century: sparse detections in the early 1900s and an increase in reported records after 2000. Over 120 years, we compiled 35 documented records of California black rail occurrence in Baja California with a minimum reported count of 56 rails. Our call-broadcast survey methods reconfirmed previously occupied sites and documented the species at a previously unreported site. Surveys in the Baja California portion of the Colorado River Delta yielded 20 detections across nine survey events between 2000 and 2022. In April 2003, standardized surveys detected 12 California black rails: 7 in Bahía de San Quintín, 3 in Arroyo San Telmo, and 2 in Arroyo El Rosario. Our recent survey work (2022–2025) recorded maximum monthly counts of 31 detections at Bahía de San Quintín and 15 detections at Estero de Punta Banda, the latter representing a previously undocumented location for the subspecies. These surveys also yielded the first photographic and acoustic documentation of California black rail in México. The contribution of spatial and temporal data provided by citizen science platforms increased our understanding of the subspecies’ distribution. Enforceable protection of this subspecies’ habitat is needed to safeguard its persistence.

Baja California↗