Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Earth and Space Science”

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

Induced earthquakes are generally not tidally triggered in Oklahoma and Kansas

Human-induced earthquakes occur along critically stressed faults as injected wastewater simultaneously heightens fluid pressure and pushes faults to failure. We investigate the possibility that small stresses imposed by Earth tides could trigger earthquakes in the induced seismicity region of Oklahoma and Kansas from 2011 to 2018. We decluster a catalog consisting of ∼110,000 earthquakes using three methods (Reasenberg, nearest-neighbor distance, and phase-bin). We find no significant tidal earthquake triggering using Schuster's p -value test for the declustered catalogs as a whole. We search for localized triggering using discretized space-time cells and find ∼0–6% of cells have significant tidal triggering which is close to what is randomly expected (5%) and indicates there is an insignificant amount of tidal triggering for the full study region. One area that has significant p -values across multiple time windows, ∼2014–2016 is ∼15 km from a region of large wastewater injection volume. It is possible that localized tidal triggering occurs for this time and area because faults remain critically stressed and are particularly susceptible to slip under the small stress load from the semidiurnal tide. Possible explanations for the lack of tidal triggering in our broader study are that the pre-seismic stressing rate in the earthquake nucleation area is faster than the tidal stressing rate (∼3 kPa/day), faults are not close enough to critically stressed to be affected by tidal forcing, and that nucleation occurs over longer periods than the tides considered in this study (∼1, ∼14 days). Fluid injection could be the source of a higher pre-seismic stress rate.

Kansas, Oklahoma↗
USGS2025

U.S. Geological Survey—Department of the Interior, Region 11, Alaska—2023–24 biennial science report

Introduction U.S. Geological Survey (USGS) Mission—The USGS national mission is to monitor, analyze, and predict the current and evolving dynamics of complex human and natural Earth-system interactions and to deliver actionable information at scales and timeframes relevant to decision makers. Consistent with the national mission, the USGS in Alaska provides timely and objective scientific information to help address issues and inform management decisions across five interconnected focus areas: Energy and Minerals; Geospatial Mapping; Natural Hazards; Water Quality, Streamflow, and Ice Dynamics; and Ecosystems. The USGS in Alaska consists of approximately 350 scientists and support staff working in 3 Alaska-based science centers. USGS science activities are also initiated by the Cooperative Research Unit and USGS centers outside Alaska. In the last 5 years, USGS research in Alaska has produced many scientific benefits resulting from more than 900 publications. Publications relevant to Alaska can be conveniently searched by keyword through the USGS Publications Warehouse at https://pubs.usgs.gov/.

Alaska↗
USGS2025

Land-based nutrient flux to a fringing reef: Insights from Ofu Island, American Samoa

Submarine groundwater discharge (SGD) is a critical driver of nutrient transport in coral reef ecosystems, shaping water quality, primary productivity, and overall reef health. This study quantifies SGD fluxes and associated nutrient dynamics in two reef flat pools within the Ofu Unit of the National Park of American Samoa: Papaloloa and Fatuana. A multi-method approach integrating unoccupied aerial system-based thermal infrared (UAS-TIR) surveys, radon-based SGD measurements, multichannel electrical resistivity tomography (ERT), and discrete water sampling was used to assess SGD rates and nutrient contributions. UAS-TIR imagery revealed cooler sea surface temperatures in both pools, indicative of SGD, with the higher fluxes observed in Papaloloa. Radon measurements revealed a strong inverse correlation between SGD rates and tidal stage, with a more immediate SGD response at Papaloloa due to its highly permeable calcareous sand and gravel substrate. In contrast, a 2–3-hour lag in SGD response at Fatuana suggests discharge from a more inland aquifer that has lower diffusivity. Nutrient concentrations correlated with temperature and salinity, confirming SGD as the dominant nutrient transport mechanism, whereas isotopic analyses indicated inputs from both groundwater and potential anthropogenic sources. Despite lower SGD flux at Fatuana, higher algal cover suggests additional factors influencing algal proliferation, including substrate availability and hydrodynamic conditions. Excess nutrient inputs from SGD may contribute to algal overgrowth, which threatens Ofu’s thermally tolerant corals by increasing competition for space and light. These findings underscore the complexity of SGD-mediated nutrient dynamics in reef environments and emphasize the need for integrated hydrological and ecological assessments to support effective reef conservation and management strategies.

American Samoa, Ofu Island, Olosega Island↗
USGS2024

Probabilistic assessment of postfire debris-flow inundation in response to forecast rainfall

Communities downstream of burned steep lands face increases in debris-flow hazards due to fire effects on soil and vegetation. Rapid postfire hazard assessments have traditionally focused on quantifying spatial variations in debris-flow likelihood and volume in response to design rainstorms. However, a methodology that provides estimates of debris-flow inundation downstream of burned areas based on forecast rainfall would provide decision-makers with information that directly addresses the potential for downstream impacts. We introduce a framework that integrates a 24 h lead-time ensemble precipitation forecast with debris-flow likelihood, volume, and runout models to produce probabilistic maps of debris-flow inundation. We applied this framework to simulate debris-flow inundation associated with the 9 January 2018 debris-flow event in Montecito, California, USA. When the observed debris-flow volumes were used to drive the probabilistic forecast model, analysis of the simulated inundation probabilities demonstrates that the model is both reliable and sharp. In the fully predictive model, however, in which debris-flow likelihood and volume were computed from the atmospheric model ensemble's predictions of peak 15 min rainfall intensity, I 15 , the model generally under-forecasted the inundation area. The observed peak I 15 lies in the upper tail of the atmospheric model ensemble spread; thus a large fraction of ensemble members forecast lower I 15 than observed. Using these I 15 values as input to the inundation model resulted in lower-than-observed flow volumes which translated into under-forecasting of the inundation area. Even so, approximately 94 % of the observed inundated area was forecast to have an inundation probability greater than 1 %, demonstrating that the observed extent of inundation was generally captured within the range of outcomes predicted by the model. Sensitivity analyses indicate that debris-flow volume and two parameters associated with debris-flow mobility exert significant influence on inundation predictions, but reducing uncertainty in postfire debris-flow volume predictions will have the largest impact on reducing inundation outcome uncertainty. This study represents a first step toward a near-real-time hazard assessment product that includes probabilistic estimates of debris-flow inundation and provides guidance for future improvements to this and similar model frameworks by identifying key sources of uncertainty.

California↗
USGS2026

Seasonal and hydrologic variation influences habitat and functional structure of stream fish assemblages

Introduction: Hydrologic variability is a key driver of ecological structure in lotic systems, shaping habitat conditions, taxonomic diversity, and the functional traits that mediate species’ persistence and performance (e.g., reproductive success). While many studies examine taxonomic responses to variation in flows, few evaluate how spatiotemporal hydrologic variation influences the functional organization within stream fish communities. Methods: We quantified seasonal habitat structure and functional trait diversity of fish assemblages across six Ozark Plateau headwater streams representing two contrasting flow regimes: Groundwater Flashy and Runoff/Intermittent Flashy. Fish and habitat data were collected seasonally during a dry year (2002) and a wet year (2003). Functional space was constructed using PCoA of morphological, ecological, and life-history traits, and functional diversity was measured using community weighted means (CWMs), functional richness (FRic), functional evenness (FEve), and functional divergence (FDiv). Results: We found that habitat structure differed strongly by flow regime and season, with Runoff/Intermittent streams exhibiting pronounced reductions in depth, area, and velocity, while groundwater streams remained structurally stable. Functional identity of assemblages was similar across flow regimes, dominated by benthic, hydrodynamic taxa with opportunistic and periodic life-history strategies. However, functional structure differed significantly: FEve and FDiv were consistently lower in Runoff/Intermittent Flashy streams in both years, indicating assemblage dominance of species with similar trait combinations and reduced trait partitioning under variable flow. FRic and taxonomic richness remained stable across seasons and flow regimes, suggesting high functional redundancy despite species turnover. Discussion: Together, results show that flow regime mediates both habitat structural stability and functional organization. As climatic warming and extreme drought increase hydrologic instability in headwaters, functional trait approaches provide a sensitive tool for detecting losses of functional roles that may not be evident by using taxonomic metrics alone.

Arkansas, Kansas, Missouri, Oklahoma↗
USGS2024

The EnMAP spaceborne imaging spectroscopy mission: Initial scientific results two years after launch

Imaging spectroscopy has been a recognized and established remote sensing technology since the 1980s, mainly using airborne and field-based platforms to identify and quantify key bio- and geo-chemical surface and atmospheric compounds, based on characteristic spectral reflectance features in the visible-near infrared (VNIR) and short-wave infrared (SWIR). Spaceborne missions, a leap in technology, were sparse, starting with the CHRIS/PROBA and EO1/Hyperion missions in the early 2000s, and providing spectroscopy data with limited spectral coverage and/or low data quality in the SWIR. Since 2019, several countries and agencies have successfully launched a number of spaceborne imaging spectroscopy systems into orbit or deployed them on the International Space Station (ISS) such as DESIS, PRISMA, HISUI, GF-5, EnMAP and EMIT. Among these recent missions, the German Environmental Mapping and Analysis Program (EnMAP) stands for its long-term development, sophisticated design with on-board calibration, high data quality requirements, and extensive accompanying science program. EnMAP was launched in April 2022 and, following a successful commissioning phase, started its operational activities in November 2022. The EnMAP mission encompasses global coverage from 80° N to 80° S through on-demand data acquisitions. Data are free and open access with 30 m spatial resolution, a high spectral resolution with a spectral sampling distance of 6.5 nm and 10 nm in the VNIR and SWIR regions respectively, and a high signal-to-noise ratio. In this paper, we aim to present the mission's current status, coverage, science capabilities and performance two years after launch. We show the potential of EnMAP for space-based imaging spectroscopy to operate in various environments, including high and low light levels, dense forests, Antarctic glaciers, and arid agricultural areas. EnMAP enables various applications in fields such as agriculture and forestry, soil compositional, raw materials, and methane mapping, as well as water quality assessment, and snow and ice properties. The results show that EnMAP's performance exceeds the mission requirements, and highlights the significant potential for contribution to scientific exploitation in various geo- and biochemical sciences. EnMAP is also expected to serve as a key tool for the development and testing of data processing algorithms for upcoming global operational missions.

Remote Sensing of Environment↗
USGS2025

Estimation of dynamic geologic CO2 storage resources in the Illinois Basin, including effects of brine extraction, anisotropy, and hydrogeologic heterogeneity

Since the vast majority of carbon dioxide (CO 2 ) storage resources in the United States are in deep saline aquifers, optimizing the use of these saline storage resources could be crucial for efficient development of geologic CO 2 storage (GCS) resources and basin- or larger-scale deployment of GCS in the country. Maximum CO 2 injection rates can be enhanced by extracting brine from the CO 2 storage unit. However, disposal of the extracted brine is both a technological and economic challenge. The lowest-cost option would likely be reinjection of the extracted brine into another formation above or below the CO 2 storage unit. Therefore, it is important to estimate brine injectivity as it will constrain the potential to increase CO 2 injectivity at an injection site that has access to multiple geologic storage units where either CO 2 or brine can be injected. Using a simulation-optimization framework, coupled with a non-isothermal, multiphase CO 2 -water-salt equation-of-state module, we developed a computationally efficient method for evaluating optimization of simultaneous CO 2 injection, brine extraction, and brine (re)injection at hypothetical injection sites deployed across a geologic basin. The Illinois basin is ideal for testing our methodology because it contains multiple geologic storage units with seals in between them to isolate injection of CO 2 in one unit from interfering with the injection of either brine or CO 2 in another unit above or below it. In addition, we investigated the relative effects of variation in key geologic parameters as well as two reservoir structures (hydrogeologic heterogeneity/anisotropy and homogeneity/isotropy) on CO 2 injectivities and enhancement of CO 2 injectivity through extracting brine. Results suggest that permeability, depth, and especially thickness of the storage unit could be the most influential parameters determining CO 2 injectivity. They also suggest that only injecting CO 2 into the storage unit with the greatest injectivity, enhancing that unit’s injectivity by extracting brine, and disposing of the produced brine in other suitable units could maximize total CO 2 injectivity in limited regions of the basin. At the majority of simulated injection sites, however, we found that injecting CO 2 into all of the accessible and suitable storage units was more likely to maximize the CO 2 storage resource.

Illinois, Kentucky, Indiana↗
USGS2026

Cajon Pass and the southern San Andreas Fault System: Earthquake cycle stress accumulation and present-day loading

With over a century since the last major rupture affecting the wider Los Angeles region, tectonic stress has steadily built along the southern San Andreas and San Jacinto fault systems, raising concerns of an imminent large earthquake. Cajon Pass, located at the junction of these faults, represents a critical site for potential through-going ruptures in Southern California. We constructed new 4D earthquake cycle simulations using a 1000-year paleoseismic rupture history of the San Andreas Fault System (SAFS) to assess spatial and temporal variations in stress. A semi-analytic Fourier transform model was used to compute stress from 3D dislocations in an elastic plate overlying a Maxwell viscoelastic half-space, assuming a complete coseismic reset of resolved shear stress on ruptured elements. Results show highest stress accumulation north of Cajon Pass (∼1.8 MPa/100 years) due to greater slip rates, and lower rates south of Cajon Pass (∼1.0–1.5 MPa/100 years). By 2025, Coulomb stress is estimated at 2.8 MPa on the Mojave South (MOS) segment, 1.8 MPa on the North San Bernardino (NSB1) segment and 3.6 MPa on the San Jacinto Bernardino (SJB) segment. Segments accumulate stresses with characteristic ranges of pre-event stress interpreted as failure thresholds: 1.2–2.7 MPa for MOS, 0.4–1.6 MPa for NSB1, and 1.2–2.9 MPa for SJB. When the stress disparity between segments SJB and MOS narrows, the faults appear to rupture jointly, suggesting that stress levels may control how Cajon Pass acts as an earthquake gate. These results may inform seismic hazard assessments by linking stress evolution to fault interactions.

California↗
USGS2026

Numerical modeling of Late Pleistocene to Holocene earthquake-induced progressive rock slope damage in the paraglacial Serpentine valley, Prince William Sound, Alaska

Landslides in deglaciating fjords pose a potential tsunami threat to nearby communities; however, processes contributing to long-term progressive rock damage and landslide conditioning remain poorly constrained in paraglacial settings. Here, we analyse the role of earthquake-induced rock mass damage over late Pleistocene to Holocene time scales, as a conditioning factor for modern landslides, using distinct element numerical modelling to assess spatial and temporal patterns of fracture propagation influenced by varying glacier thickness. Conceptualised numerical models were parameterised by in situ rock mass, glacial and topographic conditions in Serpentine valley, located in Prince William Sound, Alaska, where several large landslides are actively developing along the western valley wall. Results show that although rigid glacier buttressing reduces co-seismic rock mass damage, it does not suppress it completely, and damage occurs both above and below the glacier surface elevation. We demonstrate that topography, and especially steep slopes with topographic convexity, as well as preexisting damage of joint networks and faults inherited from tectonic and exhumation induced stresses, exert primary control on the location of new co-seismic damage. Simulations representing a simplified deglaciation sequence over the past 25 ka, with a series of 10 evenly spaced earthquakes, generated rock mass damage patterns that qualitatively match in situ landslide structural and kinematic observations at one instability in Serpentine valley. Our conceptual study helps clarify the role of repeated seismicity over glacial timescales as a long-term conditioning process for paraglacial rock slope failure, highlighting spatial and temporal patterns of progressive damage accumulation, with outcomes relevant for modern landslide hazard assessment.

Alaska↗
USGS2024

A benchmark dataset and workflow for landslide susceptibility zonation

Landslide susceptibility shows the spatial likelihood of landslide occurrence in a specific geographical area and is a relevant tool for mitigating the impact of landslides worldwide. As such, it is the subject of countless scientific studies. Many methods exist for generating a susceptibility map, mostly falling under the definition of statistical or machine learning. These models try to solve a classification problem: given a collection of spatial variables, and their combination associated with landslide presence or absence, a model should be trained, tested to reproduce the target outcome, and eventually applied to unseen data. Contrary to many fields of science that use machine learning for specific tasks, no reference data exist to assess the performance of a given method for landslide susceptibility. Here, we propose a benchmark dataset consisting of 7360 slope units encompassing an area of about 4,100 km 2 "> 4,100 km 2 in Central Italy. Using the dataset, we tried to answer two open questions in landslide research: (1) what effect does the human variability have in creating susceptibility models; (2) how can we develop a reproducible workflow for allowing meaningful model comparisons within the landslide susceptibility research community. With these questions in mind, we released a preliminary version of the dataset, along with a “call for collaboration,” aimed at collecting different calculations using the proposed data, and leaving the freedom of implementation to the respondents. Contributions were different in many respects, including classification methods, use of predictors, implementation of training/validation, and performance assessment. That feedback suggested refining the initial dataset, and constraining the implementation workflow. This resulted in a final benchmark dataset and landslide susceptibility maps obtained with many classification methods. Values of area under the receiver operating characteristic curve obtained with the final benchmark dataset were rather similar, as an effect of constraints on training, cross–validation, and use of data. Brier score results show larger variability, instead, ascribed to different model predictive abilities. Correlation plots show similarities between results of different methods applied by the same group, ascribed to a residual implementation dependence. We stress that the experiment did not intend to select the “best” method but only to establish a first benchmark dataset and workflow, that may be useful as a standard reference for calculations by other scholars. The experiment, to our knowledge, is the first of its kind for landslide susceptibility modeling. The data and workflow presented here comparatively assess the performance of independent methods for landslide susceptibility and we suggest the benchmark approach as a best practice for quantitative research in geosciences.

Earth-Science Reviews↗
USGS2024

Editorial: From cold seeps to hydrothermal vents: Geology, chemistry, microbiology, and ecology in marine and coastal environments

This Research Topic compiles contemporary studies on cold seeps, hydrothermal vents, mud volcanoes, and related seafloor features that are associated with focused fluid emissions and the transfer of carbon, other chemical species, and sometimes heat from the geosphere to the ocean. Because these features sometimes tap fluids and gas originating kilometers below the seafloor, they provide an important window into deep processes that are otherwise inaccessible to scientists. At the shallow portion of their journey, migrating fluids nearing the seafloor contribute to a range of unique biological, physical, and chemical processes within the sediments themselves and at the sediment-water interface. Seafloor fluid emissions play a critical role in global biogeochemical cycles, ocean chemistry, and possibly even climate change. Seafloor leakage points often emit hydrocarbon gases (especially methane and CO 2 ) and are sometimes the loci for deposition of seafloor minerals that have economic value. A burgeoning area of research focuses on natural products generated at these features, seeking compounds with potential pharmaceutical or other applications. Multidisciplinary studies have become routine for characterization of seafloor fluid emission sites, attesting to the inseparability of geologic, physical, chemical, and biological processes in these settings. It is increasingly common for researchers to combine in a single research cruise: subbottom imaging and seafloor mapping; porewater and water column geochemistry and gas sampling; sediment retrieval for lithologic, biostratigraphic, and solid phase analyses; and studies of benthic and subseafloor communities at the microbial to macrofaunal scales. This multidisciplinary approach has the advantage of ensuring the spatial and temporal coincidence of surveys and samples, an important factor at highly dynamic seafloor fluid emission sites. In addition, researchers often use remotely operated vehicles (ROVs), autonomous underwater vehicles (AUVs), or human-occupied vehicles (HOVs) to record video of the seafloor, compile photomosaics, collect targeted samples, and survey with high-resolution geophysical near-seafloor systems, providing a degree of detail about seafloor fluid emission sites that is unprecedented compared to most areas of the deep ocean. While rarer, long-term cabled observatories or shorter-term deployments of portable observatories are also used at some loci for seafloor fluid flux and are particularly helpful for capturing temporal variations at these dynamic features. Here we summarize the Research Topic’s contribution to multidisciplinary seafloor emission studies in the categories of cold seeps, mud volcanoes, and hydrothermal vents. Figure 1 shows the geographic distribution of the studies in this Research Topic and key features referred to in this Introduction.

Frontiers in Earth Science↗
USGS2026

Evaluating groundwater quality influences from oil field operations and other anthropogenic activities in an urban setting, Santa Fe Springs, California

Groundwater quality is often affected by anthropogenic activities in urban settings. This study examines groundwater quality in and around the Santa Fe Springs Oil Field in Los Angeles County, California, where oil and gas production commonly intersects with high density industrial, commercial and residential land uses. Utilizing a combination of new and historical data, we evaluated potential pathways that would allow for oil field formation fluids to migrate into groundwater and whether mixing may have occurred based on the distribution of groundwater and oil field formation fluid tracers in samples. Samples were analyzed for a wide array of constituents including volatile organic compounds, light hydrocarbons, major ions, and various isotopic compositions. Despite evidence of oil field infrastructure providing potential pathways of migration via uncemented annular spaces, casing breaches and historical disposal of oil field formation water in surface ponds, the distribution and occurrence of stable isotopes of water, chloride, boron, and total dissolved solids do not indicate mixing of oil field formation water and groundwater. However, methane isotopic signatures and the presence of heavier alkanes suggest gas from oil-bearing formations have migrated from depth via oil field well infrastructure. Volatile organic compound detections were mainly from manufactured compounds unrelated to oil and gas production, with a relatively limited number of petroleum hydrocarbons also detected. Volatile organic compounds were generally found in wells tapping shallow, modern aged groundwater, indicating anthropogenic activities occurring at or near land surface as the source. Study results suggest that while oil field infrastructure provides migration pathways for oil field formation fluids to be introduced into groundwater, urban land uses not related to oil and gas production are the primary drivers of groundwater quality degradation.

California↗
USGS2026

Guidelines for producing integrated 210Pb and 14C age-models

Accurate reconstructions of past environmental changes are crucial in paleoecological research and require reliable chronologies of sedimentary archives. Establishing robust age-models and obtaining the most appropriate proxies for analysis is a complex scientific endeavor, requiring extensive resources and collaboration among specialists, including radiochronologists. Radiometric dating methods, such as 210 Pb and radiocarbon ( 14 C), are frequently employed to establish chronologies in aquatic sedimentary deposits and peat bogs. In this study, we review key aspects of sampling, analysis, and the principles underlying 210 Pb and 14 C age-models, focusing on methods for developing robust joint chronologies for paleoenvironmental research. Drawing largely from the authors' experiences and group discussions during and after a scientific workshop in 2022, we discuss important considerations for site selection, sampling strategies, and radiometric dating to construct integrated 210 Pb and 14 C age-models. Using expert consensus, this group – called Paleostats – aims to provide a set of best practices for other geochronologists with this methods paper. Among our conclusions, we emphasize the importance of accounting for site-specific factors such as prior information on sedimentation rates to establish appropriate sampling and analytical strategies. The use of appropriate coring devices can minimize disturbance to sediments and ensure the core surface remains intact and preserved until sectioning. Where excess 210 Pb is expected, sectioning at intervals of ≤1 cm provides an adequate sampling resolution for 210 Pb dating. Exceptions are possible, allowing for ~2–3 cm sections in areas with confirmed high sedimentation rates (e.g., > 1 cm yr −1 ). Recovering deeper core sections for 14 C dating with sufficient overlap allows for accounting errors in depth estimates made in the field. Special attention is advised during time intervals where validation proxies, such as the human-made radionuclides 137 Cs or post-bomb 14 C, are expected, and to determine the depth of secular equilibrium between 210 Pb and 226 Ra. Radiocarbon analyses are commonly performed by accelerator mass spectrometry, and age models are constructed mainly using Bayesian statistics with Markov Chain Monte Carlo techniques (e.g., Bacon ). A Bayesian approach ( Plum ) is now available for producing 210 Pb age-models, which infers the 210 Pb ex flux, eliminates the need for selecting an equilibrium depth, and allows dating cores with incomplete 210 Pb ex inventory. Plum offers improved chronologies by integrating raw 210 Pb and 14 C data, and these age-models can be enriched with other dating methodologies, such as identifying tephras and other well-recorded historical events. Harmonized reporting would contribute to making radiometric age-models reproducible, which would benefit from an international effort. Using 210 Pb and 14 C to produce integrated age-models may yield better insights into the interplay between natural and recent anthropogenic forcings on ecosystems. This can enhance our understanding of environmental processes and their impacts on climate change, ultimately supporting science-based assessments and decisions.

Earth-Science Reviews↗
USGS2025

Constraining landslide frequency across the United States to inform county-level risk reduction

Informative landslide hazard estimates are needed to support landslide mitigation strategies to reduce landslide risk across the United States. Whereas existing national-scale landslide susceptibility products assess where landslides are likely to occur, they do not address how often , which is a critical element of landslide hazard and risk assessments. In particular, the U.S. Federal Emergency Management Agency's National Risk Index (NRI) requires landslide frequency estimates to inform expected annual loss estimates. We present county-level landslide frequency (landslides per area per year) estimates for the 50 US states. We applied Bayesian negative binomial regression to estimate both the expected (average) reported landslide frequency and full distribution of annual landslide counts for each county. We compared a suite of models that used combinations of landslide-susceptible area, probability of potentially triggering earthquakes, frequency of potentially triggering precipitation, and ecological region as predictors. We trained our models with landslide inventory data from counties with the most comprehensive records available nationwide and used zero-inflated negative binomial distributions as an incompleteness model to correct for temporal reporting gaps. We selected a preferred frequency model to inform the NRI based on information criteria and physically plausible parameter estimates. The model showed that average annual reported landslide frequencies vary by 5 orders of magnitude across US counties, ranging from 0.002 (0.00015–0.05) landslides 1000 km −2 yr −1 in Kusilvak Census Area, Alaska, to 29 (19–46) landslides 1000 km −2 yr −1 in Lake County, California, reflecting the country's strong variations in landslide susceptibility, earthquake probability, and other factors for which ecological region serves as a proxy. Counties with estimated frequencies in the top 20 % of all counties are predominately along the West Coast of the continental United States, in mountainous regions of the Pacific Northwest and Intermountain West, in locally steep or earthquake-prone regions of the Midwest and Southeast, along the Appalachians, in southern and southeastern Alaska, and on some Hawaiian islands. By examining the number of landslides predicted in 99th percentile years for each county, we identified that 26 % of US counties likely have potential for widespread landsliding with more than 10 landslides 1000 km −2 yr −1 , even when such large events have not been reported in the training data for that county. Overall, our results better represent the range of possible landslide frequencies and spatial variations than previous national-scale estimates reported in the NRI, and our approach can inform other risk-reduction and loss-mitigation efforts across the United States and globally.

Natural Hazards and Earth System Sciences↗
USGS2025

System characterization report on the Environmental Mapping and Analysis Program (EnMAP)

This report addresses system characterization of the Environmental Mapping and Analysis Program hyperspectral sensor by the DLR (German Aerospace Center, ground segment project management), GFZ (Deutsches Geoforschungszentrum, science lead) and is part of a series of system characterization reports produced and delivered by the U.S. Geological Survey Earth Resources Observation and Science Cal/Val Center of Excellence. These reports present and detail the methodology and procedures for characterization; present technical and operational information about the EnMAP hyperspectral sensor; and provide a summary of test measurements, data retention practices, data analysis results, and conclusions. The Earth Resources Observation and Science Cal/Val Center of Excellence system characterization team completed data analyses to characterize the geometric (interior and exterior), and radiometric performances of the EnMAP hyperspectral sensor. Results of these analyses indicate that the Environmental Mapping and Analysis Program has a band-to-band geometric performance in the range of −0.135 to 0.15 pixel, geometric performance relative to the Operational Land Imager in the range of −27.716 meters (−0.92 pixel) to 32.892 meters (1.09 pixels) offset in comparison to Landsat 8 Operational Land Imager, offset of a radiometric comparison in the range of −0.012 to 0.020, slope of a radiometric comparison in the range of 0.947 to 1.031.

Open-File Report↗
USGS2025

Airborne geophysics for geologic mapping of critical mineral systems in the United States southern midcontinent

The increased demand for clean energy technology and a significant reliance on foreign supply chains have given impetus to understanding critical mineral systems and locating potential resources within the United States. At least thirteen critical mineral-bearing systems have been identified throughout the U.S. southern Midcontinent (Hofstra and Kreiner, 2020) but much of the region’s geologic framework is concealed by vegetation and sedimentary cover that hinder traditional geologic mapping efforts. Airborne geophysical data provide an effective way to overcome these obstacles and to provide additional insight into the deeper structures that underlie shallow mineralization. However, legacy airborne magnetic and radiometric data were collected using now-outdated instruments and methods, inconsistent survey parameters, and large flight-line spacings resulting in low-resolution data that present challenges to regional-scale study and interpretation. Over the last decade, the U.S. Geological Survey Earth Mapping Resources Initiative (EMRI) and National Cooperative Geologic Mapping Program have conducted a series of high-resolution airborne magnetic and radiometric surveys across the southern Midcontinent (Fig. 1) as part of an effort to improve understanding of the geophysical framework and natural resource potential in the region. These surveys are designed using modern survey methods and instruments with consistent parameters for flight-line spacing and flight height relative to magnetic sources. The EMRI airborne surveys are planned in collaboration with State geological surveys based on focus areas (Dicken et al., 2022) according to the presence of or potential for critical mineral deposits. High-resolution airborne magnetic and radiometric data cover focus areas such as the southeast Missouri iron metallogenic province and South-Central iron-oxide-apatite (IOA) – iron-oxide-copper-gold (IOCG) province, the Magnet Cove alkaline-carbonatite complex, the Midwest Permian ultramafic dike district, the Illinois-Kentucky fluorspar district, and several Mississippi Valley-type lead-zinc deposits and districts (Fig. 1). These focus areas represent known deposits or prospective host systems of critical minerals including rare earth elements (REEs), platinum-group elements (PGEs), cobalt, lithium, fluorspar, niobium, titanium, vanadium, lead, zinc, gallium, germanium, and many more. Other significant geologic and geophysical features covered include the Reelfoot rift, the New Madrid seismic zone, the Illinois basin, the Arkoma basin, the South-Central magnetic lineament, and the Kentucky-Tennessee magnetic anomaly (Fig. 1). This presentation focuses on new airborne magnetic and radiometric data with continuous coverage across parts of six states, preliminary interpretations, examples of geologic mapping applications, and discussion of newly discovered magnetic anomalies and follow-up investigations.

Conference Paper↗
USGS2026

Climate-adaptive urban planning: Quantitative assessment of drought impacts and practical strategies for climate-resilient urban green spaces

Urban green spaces (UGSs) are vital for enhancing a city’s resilience and livability; however, their functionality is increasingly jeopardized by drought, particularly in water-scarce regions. This study evaluates drought impact on UGSs in Metropolitan Adelaide, Australia, a representative semi-arid urban system, using satellite-derived Normalized Difference Vegetation Index (NDVI) time-series data spanning 2000–2020. Vegetation dynamics were analyzed through Seasonal-Trend decomposition using Loess (STL), standardized anomaly assessment, lagged Pearson correlation, Ordinary Least Squares (OLS) regression, and Mann–Kendall trend analysis. To isolate climatically sensitive signals, 29 urban lawn patches were examined separately from mixed urban canopy, given their shallow root systems and direct dependence on surface moisture. NDVI declined by approximately 0.09 units during the Millennium Drought (2001–2009), with summer greenness deficits reaching 24% below the 20-year benchmark. Temperature was the dominant driver of lawn NDVI variability (r = −0.863, R 2 = 74.5%), substantially exceeding the effect of rainfall (r = 0.156, R 2 = 2.4%). El Niño–Southern Oscillation (ENSO) cycles modulated vegetation responses, with La Niña years supporting recovery and El Niño years amplifying decline. Post-drought recovery remained incomplete, with NDVI deficits of 8–20% persisting through 2020; full recovery was observed only in 2017, coinciding with the highest recorded summer rainfall. No significant directional trend was detected over the full study period (Mann–Kendall τ = 0.005, p = 0.908). These findings demonstrate that heat, rather than water limitation alone, is the primary driver of vegetation stress in urban systems, highlighting the benefits of integrated management strategies that address both warming and moisture deficits to sustain urban green infrastructure under future climate conditions. We introduce the concept of “urban greenery drought,” referring to a form of vegetation stress in managed urban landscapes where greenness is reduced primarily by elevated temperature and atmospheric demand despite water availability.

Adelaide↗
USGS2025

Refining the earthquake history of south-central Alaska through lake records

The Alaska–Aleutian subduction zone (AASZ) is one of the world's most seismically active plate boundaries and the source of the 1964 M w 9.2 Great Alaska earthquake–the second largest instrumentally recorded earthquake in the world. Understanding the nature and frequency of such earthquakes is necessary for seismic and tsunami hazard assessment, but instrumental and historical records that span less than 150 years are too short to allow a statistically reliable analysis of earthquake recurrence times. This calls for studies of evidence of past earthquakes, extending the earthquake catalog further back in time. Subduction-zone paleoseismology in south-central Alaska is predominantly based on coastal evidence of land-level changes and tsunamis generated by megathrust earthquakes and preserved in the geological record. A complementary approach is lacustrine paleoseismology, which is still a relatively young discipline in Alaska. However, globally, lake basins are well-established high-resolution and continuous recorders of paleoseismic activity along subduction zones, relying on the identification of underwater landslide deposits and turbidites generated by seismic shaking. As a result, lake basins not only register ground shaking from megathrust earthquakes, but also from intraslab and crustal earthquakes, which are typically not accompanied by significant land-level changes. In this review paper, we combine coastal and lacustrine paleoseismology approaches to refine the south-central Alaskan earthquake history by comparing the paleoseismic records from two lakes (i.e., Eklutna Lake, located in the Chugach Mountain Range, and Skilak Lake, situated on the Kenai Peninsula) with the coastal and crustal earthquake catalog in Alaska. The resulting age ranges of all known megathrust earthquakes involving the Alaskan megathrust between the Kodiak and Prince William Sound (PWS) sections are more precise and accurate for the last 1.3 kyrs BP than the previously published age ranges from coastal records. As a result, this study supports the following key conclusions: (1) The 1964 CE earthquake was an exceptionally strong and unique event in the last 2000 years, rupturing the PWS, Kenai, Barren Islands, and Kodiak sections simultaneously. (2) The high-resolution and seasonal markings of the varved lake records now disentangle for the first time closely timed earthquakes, which was not possible based on the coastal evidence alone. (3) No persistent megathrust rupture boundaries exist. So, the possibility of a full rupture of the entire eastern AASZ, from PWS to Semidi cannot be excluded. (4) The rupture pattern in the eastern AASZ reveals superimposed cycles of multi-asperity ruptures (1964 earthquake) and clustered complementary partial ruptures, or rupture cascades. (5) The PWS section hosts the largest asperity in the eastern AASZ. (6) The shaking record of megathrust earthquakes indicates a time-dependent (quasiperiodic) behavior for the study area, but the observation of complementary clusters means that the hazards will not drop to zero but instead may even increase for a neighboring section. (7) The time-independent behavior of intraplate earthquakes implies that the intraslab hazard did not decrease following the 2016 and 2018 earthquakes. This study utilizes an integrated approach for subduction zone paleoseismology as a solution for unraveling recurrence and rupture patterns in Alaska, which can be applied worldwide.

Alaska↗