Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Space Science Reviews”

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.

373 records · Page 20Linked to original sources
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↗
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↗
USGS1999

Concentrations and distribution of manmade organic compounds in the Lake Tahoe Basin, Nevada and California, 1997-99

The U.S. Geological Survey, in cooperation with the Tahoe Regional Planning Agency and the Lahontan Regional Water-Quality Control Board, sampled Lake Tahoe, major tributary streams to Lake Tahoe, and several other lakes in the Lake Tahoe Basin for manmade organic compounds during 1997-99. Gasoline components were found in all samples collected from Lake Tahoe during the summer boating season. Methyl tert -butyl ether (MTBE), benzene, toluene, ethylbenzene, and xylenes (BTEX) were the commonly detected compounds in these samples. Most samples from tributary streams and lakes with no motorized boating had no detectable concentrations of gasoline components. Motorized boating activity appears to be directly linked in space and time to the occurrence of these gasoline components. Other sources of gasoline components to Lake Tahoe, such as the atmosphere, surface runoff, and subsurface flow, are minor compared to the input by motorized boating. Water sampled from Lake Tahoe during mid-winter, when motorized boating activity is low, had no MTBE and only one sample had any detectable BTEX compounds. Soluble pesticides rarely were detected in water samples from the Lake Tahoe Basin. The only detectable concentrations of these compounds were in samples from Blackwood and Taylor Creeks collected during spring runoff. Concentrations found in these samples were low, in the 1 to 4 nanograms per liter range. Organochlorine compounds were detected in samples collected from semipermeable membrane devices (SPMD's) collected from Lake Tahoe, tributary streams, and Upper Angora Lake. In Lake Tahoe, SPMD samples collected offshore from urbanized areas contained the largest number and highest concentrations of organochlorine compounds. The most commonly detected organochlorine compounds were cis- and trans-chlordane, p , p' -DDE, and hexachlorobenzene. In tributary streams, SPMD samples collected during spring runoff generally had higher combined concentrations of organochlorine compounds than those collected during baseflow conditions. Upper Angora Lake had the fewest number of organochlorine compounds detected of all lake samples. Dioxins and furans were not detected in SPMD samples from two sites in Lake Tahoe or from two tributary streams. The number of polycyclic aromatic hydrocarbon (PAH) compounds and their combined concentrations generally were higher in samples from Lake Tahoe than those from tributary streams. Areas of high-motorized boating activity at Lake Tahoe had the largest number and highest concentrations of PAH's. PAH compounds were detected in samples from SPMD's in four of six tributary streams during spring runoff, all tributary streams during baseflow conditions, and at all lake sites. The most commonly detected PAH's in tributary streams during spring runoff were phenanthrene, fluoranthene, pyrene, and chrysene, and during baseflow conditions were phenanthrene, 1-methylphenanthrene, diethylnaphthalene, and pyrene. Upper Truckee River, which has an urban area in its drainage basin, had the largest number and highest combined concentration of PAH's of all stream samples. Bottom-sediment from Lake Tahoe had detectable concentrations of p-cresol, a phenol, in all but one sample. A sample collected near Chambers Lodge contained phenol at an estimated concentration of 4 micrograms per kilogram (µg/kg). Bottom-sediment samples from tributary streams had no detectable concentrations of organochlorine or PAH compounds. Several compounds were detected in bottom sediment from Upper Angora Lake at high concentrations. These compounds and their concentrations were p , p' -DDD (10 µg/kg), p , p' -DDE (7.4 µg/kg), 2,6-dimethylnaphthalene (estimated at 190 µg/kg), pentachlorophenol (3,000 µg/kg), and p-cresol (4,400 µg/kg).

California, Nevada↗
USGS2024

U.S.-Mexico Borderland & vegetation community map

People on both sides of the United States-Mexico border need a high-resolution, binational vegetation community map that spans the entire United States-Mexico borderlands. Traditionally, mapping efforts in this region were impeded by complex logistics related to the international border, differing national needs and plans, and resource allocations and priorities. To address this need, scientists from the U.S. Geological Survey (USGS) Southwest Biological Science Center partnered with the Sonoran Joint Venture, the U.S. Fish and Wildlife Service (FWS) Migratory Bird Program, data engineers from the Department of Biosystems Engineering at the University of Arizona, and collaborators from the Wildlands Network, the Borderlands Program to produce the first prototype land cover map within the overlapping Mojave Desert, Sonoran Desert, and the North American Bird Conservation Initiative’s Bird Conservation Region 33 (BCR33) using Landsat satellite data . BCR33 is an area of high biodiversity, providing habitat for bird species of concern and other wildlife. The land cover map supports FWS recovery plan efforts related to conservation planning activities for many species, including Yellow-billed Cuckoo ( Coccyzus americanus ), Cactus Ferruginous Pygmy-Owl ( Glaucidium brasilianum cactorum ), Southwestern Willow Flycatcher ( Empidonax traillii extimus ), Yuma Ridgway’s Rail ( Rallus obsoletus yumanensis ), Bendire’s thrasher ( Toxostoma bendirei ), LeConte’s thrasher ( Toxostoma lecontei ), Masked Bobwhite ( Colinus virginianus ridgwayi ), jaguar ( Panthera onca ), and endangered plants such as Bartram’s stonecrop ( Graptopetalum bartramii ) and the Pima pineapple cactus ( Coryphantha robustispina ssp. robustispina ). In 2024, a Phase-II map for the full BCR33 region was completed, increasing the understanding of the binational nature of natural communities. The published map and associated paper can be found here .

Borderland↗
USGS2025

Altitude of the potentiometric surface and depth to water in the Mississippi River Valley alluvial aquifer, spring 2022

Potentiometric-surface and depth-to-water maps for spring 2022 were created for the Mississippi River Valley alluvial aquifer (MRVA) using groundwater-altitude data from 1,136 wells completed in the MRVA and from the altitude of the top of the water surface in area rivers from 160 streamgages. The potentiometric-surface and depth-to-water maps for 2022 were created to support investigations to characterize the MRVA as part of the U.S. Geological Survey Water Availability and Use Science Program. Sufficient data were available to map the potentiometric surface and depth to water of the MRVA for spring 2022 for about 83 percent of the aquifer area. The potentiometric contours ranged from 0 to 340 feet (ft) above the North American Vertical Datum of 1988. The regional direction of groundwater gradient was generally to the south-southwest, except in areas of groundwater-altitude depressions, where the groundwater gradient was into the depression, and near rivers, where the groundwater gradient can be from aquifer to the river or from the river into the aquifer. There are large depressions in the potentiometric-surface map in the lower one-half of the Cache region and in much of the Grand Prairie and Delta regions. Depth to water in the MRVA, spring 2022, by well ranged from 5.00 ft above land surface to 145.66 ft below land surface.

Mississippi River Valley alluvial aquifer↗
USGS2024

Conceptual hydrogeologic framework and groundwater budget near the southeastern part of Puget Sound, Washington

More than 1 million people live within the active model area (AMA) in the southeastern part of the lowlands surrounding Puget Sound, or Puget Lowland, Washington, and groundwater is the source for approximately one-half of their public, domestic, and irrigation water demands. The 887-square-mile AMA, located in King and Pierce Counties, represents the area of analysis for the conceptual hydrogeologic framework and numerical groundwater-flow models within the study area and includes the Puyallup River and Chambers-Clover Creek watersheds. To assess the potential hydrologic and anthropogenic impacts to groundwater and the connected surface-water resources, conceptual and numerical groundwater-flow models of groundwater flow were developed by the U.S. Geological Survey Washington Water Science Center in close cooperation with 18 water-resource agencies and stakeholders. This report presents information used to characterize the groundwater-flow system and the development of a numerical model in the AMA. Included are descriptions of the geology and conceptual hydrogeologic framework, groundwater levels and flow directions, groundwater recharge and discharge, numerical groundwater-flow model construction and results, and model limitations. The study area encompasses the western part of Pierce County and the southwestern part of King County, Washington. The study area extends south to the Nisqually River, southwest to Tanwax Creek, northeast to the Green River, and north through the valley near Auburn and adjacent uplands. It is bounded on the east by foothills of the Cascade Range, and on the northwest by Puget Sound.

Washington↗
USGS2026

Earthquake rupture arrest from depth‐dependent frictional stability on the Pütürge segment of the East Anatolian Fault Zone, Türkiye (Turkey)

Determining why earthquake ruptures stop is a central challenge in earthquake science and seismic hazard assessment. The Pütürge segment of the East Anatolian Fault Zone, Türkiye, exhibits shallow creep (<∼3 km depth) yet hosts M > 6.5 earthquake ruptures at greater depth. Here, we evaluate whether variations in frictional stability along this segment aided arrest of the 2020 M 6.7 Elaziğ and 2023 M 7.8 Pazarcık earthquake ruptures. Analysis of Sentinel‐1 Synthetic Aperture Radar imagery indicates the 2023 M 7.8 rupture propagated laterally into a metamorphic massif within the Pütürge segment, where slip rapidly decayed below detection limits. Creepmeters along this segment recorded no significant surface afterslip (<3 mm) in the following year. To investigate this fault‐slip behavior, we conducted triaxial friction experiments on Pütürge fault gouge sampled from an outcrop exposure. The gouge, composed primarily of muscovite, quartz, and calcite, is velocity strengthening at conditions approximating 0–2.5 km depth and velocity weakening at 4–5 km depth. This transition to velocity‐weakening friction is associated with enhanced comminution and shear localization observed microstructurally. Our results suggest that depth‐dependent frictional stability of the Pütürge fault segment facilitates rupture nucleation and propagation at depth while maintaining rupture‐arresting behavior near Earth’s surface.

Pütürge segment of the East Anatolian Fault Zone↗
USGS2025

A probabilistic assessment methodology for the evaluation of geologic energy storage capacity—Natural gas storage in depleted hydrocarbon reservoirs

The need for energy storage, particularly underground, where capacity and duration may far exceed battery storage technologies, is especially relevant given the increasing demands for reliable power alongside the development of intermittent renewable electricity sources. Geologic energy storage facilities already exist, and expanded use would enable storing gases such as methane and hydrogen. In 2018, a National Academies of Sciences, Engineering, and Medicine report, “Future Directions for the U.S. Geological Survey's Energy Resources Program,” recommended that the U.S. Geological Survey (USGS) prioritize assessing underground energy storage in geologic formations in the United States. The U.S. Geological Survey has since developed a methodology for assessing natural gas storage capacities in depleted hydrocarbon reservoirs on a national scale. The methodology introduced in this report prescribes three approaches for calculating gas storage capacity. This methodology relies on the availability of input data, including cumulative hydrocarbon production records, reservoir petrophysical properties, and reservoir pressure data. Assessment inputs can be obtained from public, State-level databases and propriety national-scale databases, although the use of analogs could be warranted for estimating input parameters. Probabilistic assessment results are aggregated to play, petroleum province, regional, and national scales. The steps defined in this report are demonstrated on the Michigan Basin Province, which includes the Mississippian Sandstone Gas Play and the Clinton Structural Play. This methodology could be used to systematically and consistently assess hydrocarbon plays and provinces for natural gas storage capacity across the United States.

Michigan↗
USGS2024

U.S. Geological Survey Karst Interest Group proceedings, Nashville, Tennessee, October 22-24, 2024

Karst hydrogeologic systems represent challenging and unique conditions to scientists studying groundwater flow and contaminant transport. Karst terrains are characterized by distinct and beautiful landscapes, caverns, and springs, and many of the exceptional karst areas are designated as national or state parks. The range and complexity of landforms and groundwater flow systems associated with karst terrains are enormous, perhaps more than any other aquifer type. The U.S. Geological Survey (USGS) Karst Interest Group (KIG), formed in 2000, is a loosely knit, grassroots organization of USGS and non-USGS scientists and researchers devoted to fostering better communication among scientists working on, or interested in, karst aquifers. The primary mission of the KIG is to encourage and support interdisciplinary collaboration and technology transfer among scientists working in karst areas. To accomplish its mission, the KIG has organized a series of workshops. To date (2024), nine KIG workshops, including the workshop documented in this report, have been held. The abstracts and extended abstracts provide a snapshot in time of past and current karst related studies. The USGS Water Availability and Use Science Program funded the workshop and proceedings. The planning committee for the ninth workshop includes Thomas D. Byl (USGS and Tennessee State University), Allan K. Clark (USGS), Laura M. DeMott (USGS), Eve L. Kuniansky (USGS, Emeritus), Benjamin V. Miller (USGS), and Lawrence E. Spangler (USGS, Emeritus). The workshop proceedings are edited by Eve L. Kuniansky and Lawrence E. Spangler. The field trip guide was produced by Benjamin V. Miller and Brian Ham (Tennessee Department of Environment and Conservation) and included in the proceedings from the KIG’s 2021 virtual workshop to be used on the optional field trip held on Thursday, October 24, 2024.

Open-File Report↗
USGS2024

Long-term communication of aftershock forecasts: The Canterbury earthquake sequence in New Zealand

On 14 February 2016, a magnitude (M)5.7 earthquake struck in Christchurch New Zealand (Aotearoa in the Maori language). The shaking caused damage to historic facades, power outages, cliff collapses, rock falls, and liquefaction but no reported injuries or fatalities. This Valentine’s Day earthquake was an aftershock in the Canterbury earthquake sequence (CES), which began on 4 September 2010 with the M7.1 Darfield Earthquake and included the destructive and fatal M6.2 Christchurch aftershock on 22 February 2011. This study, eight months after the Valentine’s Day earthquake and six years after the initiation of the CES, is the first to explore long-term aftershock forecast information and communication needs. The exploratory study also aimed to gather feedback on aftershock scenarios, an alternative form for communicating the forecast. The qualitative study involved workshops with emergency managers, public health officials, and members of the public in Christchurch. Key findings for long-term communication throughout an earthquake sequence include: 1. divergent earthquake experiences affect aftershock communication response and information needs; 2. understanding aftershock sequence behavior is foundational to sense-making when large aftershocks occur; 3. strategic earthquake sequence updates from the trusted science agency and local agencies could serve as important reminders for earthquake preparedness; 4. communication of aftershock forecast uncertainty could aid with both the credibility of the information and living with uncertainty, and 5. inclusion of impact information and preparedness advice into aftershock forecast scenarios could provide links to actionable information. The paper derives implications for research and practice of long-term communications during an aftershock sequence.

Christchurch↗
USGS2025

Mapping bedrock outcrops in the Sierra Nevada Mountains (California, USA) using machine learning

Accurate, high-resolution maps of bedrock outcrops can be valuable for applications such as models of land–atmosphere interactions, mineral assessments, ecosystem mapping, and hazard mapping. The increasing availability of high-resolution imagery can be coupled with machine learning techniques to improve regional bedrock outcrop maps. In the United States, the existing 30 m U.S. Geological Survey (USGS) National Land Cover Database (NLCD) tends to misestimate extents of barren land, which includes bedrock outcrops. This impacts many calculations beyond bedrock mapping, including soil carbon storage, hydrologic modeling, and erosion susceptibility. Here, we tested if a machine learning (ML) model could more accurately map exposed bedrock than NLCD across the entire Sierra Nevada Mountains (California, USA). The ML model was trained to identify pixels that are likely bedrock from 0.6 m imagery from the National Agriculture Imagery Program (NAIP). First, we labeled exposed bedrock at twenty sites covering more than 83 km 2 (0.13%) of the Sierra Nevada region. These labels were then used to train and test the model, which gave 83% precision and 78% recall, with a 90% overall accuracy of correctly predicting bedrock. We used the trained model to map bedrock outcrops across the entire Sierra Nevada region and compared the ML map with the NLCD map. At the twenty labeled sites, we found the NLCD barren land class, even though it includes more than just bedrock outcrops, accounted for only 41% and 40% of mapped bedrock from our labels and ML predictions, respectively. This substantial difference illustrates that ML bedrock models can have a role in improving land-cover maps, like NLCD, for a range of science applications.

California↗
USGS2024

An orbital comparison of a late mantling unit on Aeolis Mons with other erosion-resistant strata explored by MSL in Gale crater, Mars

Previous investigations along Curiosity 's traverse in Gale crater have explored the relationship between orbital and in situ observations. This work aims to better understand the geologic environment of units only observable from orbit and compare them to the properties of units examined by Curiosity . Here, we map an erosion-resistant dark-toned mantling unit that overlies the modern topography of Aeolis Mons (informally known as Mt. Sharp) and compare this unit to two other previously mapped dark-toned resistant units, the marker band and the mound skirting unit (MSU), that have been inferred to represent different geologic environments (lacustrine and aeolian, respectively). Visible to short wave infrared spectra from the Compact Reconnaissance Imaging Spectrometer for Mars and visual images from the High Resolution Imaging Science Experiment and Context Cameras aboard the Mars Reconnaissance Orbiter are used for this comparison. Spectral data suggest a mafic composition with minor alteration, although the composition varies more with location around Mt. Sharp rather than between units. Morphologically, the mantling unit has strong similarities to the marker band based on their consistent low-albedo, erosion-resistance, and smooth appearance, contrasting with the highly variable surface texture of the MSU. We hypothesize that all three units had a similar sediment source but experienced aqueous alteration at different times: early ubiquitous cementation in a surface aqueous environment in the mantling unit and marker band versus patchy late diagenesis in the MSU. If true, these results suggest that water activity continued within the Gale crater long after the erosion of Mt. Sharp.

JGR Planets↗
USGS2024

Hydrogeologic framework and extent of saltwater intrusion in Kings, Queens, and Nassau Counties, Long Island, New York

In 2016, the U.S. Geological Survey began a multiyear cooperative study with the New York State Department of Environmental Conservation to evaluate the sustainability of Long Island’s sole-source aquifer system through hydrogeologic mapping, compilation of groundwater chloride concentrations, and groundwater flow modeling. In the initial phase of the islandwide study, the hydrogeologic framework and extent of saltwater intrusion in aquifers in Kings, Queens, and Nassau Counties on western Long Island, N.Y., were investigated. The aquifer system underlying western Long Island has been under stress from pumping of public, irrigation (golf course), and industrial supply wells. Saltwater intrusion has occurred from surrounding embayments (East River, Long Island Sound, Jamaica Bay, and the Atlantic Ocean) due to pumping. Eighteen boreholes were drilled and cores taken during 2019–21 to collect hydrogeologic, geochemical, and geophysical data to delineate the complex subsurface hydrogeology and extent of saltwater intrusion within the study area. Evaluation of the new cores, reexamination of legacy core descriptions, and analysis of borehole geophysical logs was used to refine the previously published hydrogeologic framework of Pleistocene and Cretaceous unconsolidated sediments in the area, including delineation of a previously undefined hydrogeologic unit between the Magothy aquifer and the Raritan confining unit, herein named the “upper Raritan aquifer.” The upper Raritan aquifer was first recognized in southeastern Nassau County from an analysis of about 50 closely spaced boreholes with high-resolution core descriptions and gamma-ray (gamma) logs. Further analysis of borehole logs across the study area indicated that the upper Raritan aquifer was also present in Kings and Queens Counties. Nuclear magnetic resonance (NMR) logging was used for the first time on Long Island to provide estimates of the hydraulic properties of the major aquifer and confining units. Unlike other geophysical logs that record responses to the rock matrix and fluid properties and are strongly dependent on mineralogy, NMR logs record responses to the presence of hydrogen protons in the formation fluid to determine water fraction and pore-size distribution. NMR log analysis provided estimates of the clay-bound, capillary-bound, and mobile water fractions and hydraulic conductivity of aquifers and confining units penetrated by five wells in Nassau County. Pumpage for public-supply and industrial wells on Long Island began in the 1870s with small, localized suppliers of populated areas in Kings and Queens Counties. By 1904–16, pumpage for public water supply in Kings County averaged 21 million gallons per day, and averaged 37 million gallons per day in Queens County, mostly from the upper glacial aquifer. Saltwater intrusion was reported as early as the beginning of the 20th century and included the upper glacial-Jameco-Magothy and Lloyd-North Shore aquifer systems. By 1936, pumping in central Kings County created a major cone of depression in the water table extending to the south shore of much of Kings County and into southwestern Queens County. Saltwater intrusion has caused the shutdown of public-supply wells in Kings, Queens, and Nassau Counties. A large saltwater intrusion wedge in the Lloyd aquifer was indicated in southern Queens County in the early part of the 20th century, and the saltwater interface may have been onshore predevelopment. Most of Kings and Queens Counties are intruded with saltwater in both the upper glacial-Jameco-Magothy and Lloyd-North Shore aquifers systems. Saltwater increased during the 20th century and continues to increase to the present (2023) in the Lloyd-North Shore aquifer system in Great Neck and Manhasset Neck in northern Nassau County. A major wedge of saltwater intrusion in the upper glacial-Jameco-Magothy aquifer in southwestern Nassau County appears to be increasing.

New York↗
USGS2025

Comparisons of shoreline positions from satellite-derived and traditional field- and remote-sensing techniques

Satellite-derived shorelines (SDS) have the potential to help researchers answer critical coastal science questions and support work to predict coastal change by filling in the spatial and temporal gaps present in current field-based and remote-sensing data collection methods. The U.S. Geological Survey conducted comparison analyses of traditionally sourced shorelines and SDS in diverse coastal landscapes to determine how SDS could be used in ongoing and future work across varied coastal environments and provided some initial findings that could be used for implementation. Using CoastSeg, a browser-based program for SDS detection and mapping, SDS for the period 1984–2023 for multiple locations across the United States were compared to shoreline positions from traditionally sourced shoreline data. In this report, the authors present these comparisons alongside lessons learned and challenges encountered when building SDS workflows in different coastal locations. Results show that individual SDS have larger uncertainty and yet produced similar linear trends to sparser, traditionally sourced shoreline data; because SDS methods provide orders of magnitude more data than traditional shoreline-detection methods, they can be used to evaluate shoreline behaviors. Refining average scalar slopes used in tidal corrections did not result in substantial decreases in uncertainty. Using lessons from this work to outline needs for regional implementation, initial setup time would be considerable, being on the order of weeks. However, once complete, shoreline detections and analyses are fast (on the order of minutes to hours) and achievable using a desktop computer.

Alaska, Florida, Massachusetts, Washington↗
USGS2026

Altitude of the potentiometric surface and depth to water in the Mississippi River Valley alluvial aquifer, spring 2024

Potentiometric-surface and depth-to-water maps for spring 2024 were created for the Mississippi River Valley alluvial aquifer using groundwater-altitude data from 1,151 wells completed in the Mississippi River Valley alluvial aquifer and from the altitude of the top of the water surface in area rivers from 160 U.S. Geological Survey and U.S. Army Corps of Engineers streamgages. The potentiometric-surface and depth-to-water maps for 2024 were created to support investigations to characterize the Mississippi River Valley alluvial aquifer as part of the U.S. Geological Survey Water Availability and Use Science Program. Sufficient data were available to map the potentiometric surface and depth to water of the Mississippi River Valley alluvial aquifer for spring 2024 for about 83 percent of the aquifer area. The potentiometric contours ranged from 0 to 330 feet (ft) above the North American Vertical Datum of 1988. The regional direction of the groundwater gradient was generally to the south-southwest, except in areas of groundwater-altitude depressions, where the groundwater gradient direction was into the depression, and near rivers, where the groundwater gradient direction was either from the aquifer to the river or from the river into the aquifer. Groundwater depressions in the potentiometric-surface map are in the lower one-half of the Cache region and in most of the Grand Prairie and Delta regions. Depth to water by well in the Mississippi River Valley alluvial aquifer in spring 2024 ranged from 0.27 ft above land surface to 145.65 ft below land surface.

Arkansas, Illinois, Kentucky, Louisiana, Mississip↗
USGS2025

A scientific vision and roadmap for earthquake rupture forecast developments, a USGS perspective

We articulate a scientific vision and roadmap for the development of improved Earthquake Rupture Forecast models, which are one of the two main modeling components used in modern seismic hazard and risk analysis. One primary future objective is to provide fully time-dependent models that include both elastic rebound and spatiotemporal clustering nationwide, which is particularly important for shorter-term hazard and risk considerations (e.g., earthquake insurance products). We also discuss the importance and perennial challenges associated with quantifying epistemic uncertainties, including those associated with deformation-model slip rates, un-quantified sampling errors with respect to off-fault seismicity, and any spatial covariances. The need for more physics-based approaches is also emphasized, as is the benefit of adding model valuation (quantifying usefulness) to our verification and validation protocols. Given the multidisciplinary and system-level nature of this activity, modular design is critical. Future updates will also draw from best-available science by both the United States Geological Survey and the external community. The primary goal of this paper is to highlight plans that guide research and facilitate community engagement with model development, especially with respect to lowering the entry barrier for early career scientists and engineers. The paper is written so readers can focus on the sections that interest them most (see table of contents), with the Introduction and Discussion providing a stand-alone overview and summary.

Bulletin of the Seismological Society of America↗
USGS2024

Numerical model of the groundwater-flow system near the southeastern part of Puget Sound, Washington

Groundwater flow in the active model area (AMA) was simulated using a groundwater-flow model. A steady-state model version of the model simulates equilibrium conditions, and a transient model version simulates monthly variability. The model corresponds to the physical and temporal dimensions of the conceptual model and groundwater budget. The steady-state model version represents average conditions for an 11-year period (January 1, 2005–December 31, 2015), and the transient model represents monthly hydrologic variability within that period. The 13-layer model was constructed using MODFLOW-NWT with a uniformly spaced grid consisting of 416 rows, 433 columns, and cells with a horizontal dimension of 500 feet (ft) on a side. The model was calibrated to measured values of water levels in wells and lakes and estimated base flow for selected streamflow measurement stations, commonly referred to as streamgages. Model calibration was accomplished using a combination of manual and automatic methods, including the Model-Independent Parameter Estimation (PEST) program that adjusted model input parameters with the aim of minimizing the difference between estimated and model-simulated values of hydraulic head and base flow. Model boundary conditions consist of all simulated groundwater inflow to and outflow from the AMA. For example, a stream reach that simulates a gain from or loss to groundwater is a boundary condition that allows water to exit or enter, respectively, the groundwater system. Other boundary conditions include springs, seeps, precipitation recharge, groundwater exchange with lakes and Puget Sound, and groundwater pumping. A comparison of the estimated groundwater budget to that simulated by the steady-state model version indicates that the relative percentages of total inflow or total outflow for six major categories of boundary conditions are similar for the two budgets. The model was used to simulate three suites of scenarios of potential drought and water-use changes. Scenario 1 suite consisted of the steady-state model version that was run with 0, 15, 20, and 25 percent reduction of precipitation recharge to assess the corresponding reductions in base flow with decreasing recharge. The last simulation for the scenario 1 suite consisted of the transient model version simulating 3 years of consecutive seasonal drought, defined by the months of May through September, to assess the corresponding base-flow reductions. Scenario 2 suite consisted of the steady-state model version with all simulated groundwater use removed, compared with a simulation that includes current groundwater use to evaluate changes to potentiometric surfaces and base flows. Scenario 3 suite consisted of a transient model version of the model that simulated pumping increases for four different categories of water-supply wells (compared to no pumping increases) to evaluate resulting reductions in base flow. Although, these scenarios provide examples of model applications and useful insights, many other scenarios could be simulated. A description of how to download the model is described in the body of this report. Uncertainty is associated with most model inputs. Groundwater levels, lake levels, and land-surface altitudes are relatively certain; other model inputs are far less certain, including precipitation recharge, base flow, hydraulic properties, water use, and the three-dimensional structure of subsurface hydrogeologic units. Models are useful not because of high levels of accuracy of all model inputs, but because they combine the best information and estimates available, thereby providing the best predictions available related to physical processes. The model described in this report simulates groundwater flow on a regional scale, which has inherent limitations for simulating hydrologic scenarios at local scales. Model structures and inputs were generalized to be consistent with this regional scale. For example, the actual groundwater system has much greater heterogeneity of hydraulic conductivity than is possible within the model’s degrees of freedom. Variations in hydraulic gradients over distances less than 500 ft cannot be simulated. The distances between model features, such as a pumping well and a stream, must be placed at 500-ft intervals and are co-located if both features are within the same model cell.

Washington↗
USGS2025

Characterizing Meteor Crater impact melts through geochemistry and textural analysis

The U.S. Geological Survey Astrogeology Science Center houses the Meteor Crater sample collection, an assemblage of over 2,500 meters of cuttings from 161 drill holes into Meteor Crater’s rim, flanks, and ejecta blanket. We have utilized this unique collection to study the composition and spatial distribution of impact-generated materials from within the ejecta blanket. Meteor Crater has historically been known to have generated only a relatively small amount of impact melt compared to other terrestrial craters of similar size. A detailed compositional and textural dataset of impact-derived melts from this impact can therefore be a useful asset in improving our understanding of crater formation, and in particular impact melt formation. We have characterized 42 impact-melt particles from Meteor Crater using a scanning electron microscope and an electron microprobe for textural and compositional analysis. We analyzed samples from six drill holes in the ejecta blanket, situated to the northwest, southeast, south, and southwest of the crater (ejecta northeast of the crater is devoid of impact melts). Impact melts were collected from drill cuttings at various depths within the ejecta blanket, ranging from a few centimeters below the surface down to ~6.5 meters. Backscattered electron (BSE) images were acquired for each analyzed impact-melt particle. To characterize the various textures and phases present in each impact melt, we also took many detailed BSE images. Our geochemical analyses include full spectral profiles using energy dispersive X-ray spectrometry and well-calibrated wavelength dispersive spectrometry for a number of phases, including minerals (olivine, pyroxene, and so on), pristine glass, and metallic inclusions. The full dataset is available in ScienceBase as a data release (Gullikson and others, 2024), accessible at https://doi.org/10.5066/P9OGAJ8P . Our goal for this Open-File Report is to provide a summary of this immense dataset, details on data collection, descriptions of the different phases observed within impact-melt particles (both geochemically and texturally), and observable trends.

Arizona↗