Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Hydrology: Regional Studies”

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

Statewide surficial geologic map of Nebraska underscores Quaternary landscape evolution from the High Plains to the Central Lowland

Surficial geologic mapping in Nebraska has been conducted primarily at the 1:24,000 scale since the mid-1990s, although there have also been limited efforts to map generalized Quaternary and surficial geology within the state and the region. We compiled and evaluated disparate maps (1:24,000–1:1,000,000) and datasets—including geologic, soil and soil parent material, and geomorphic maps as well as LiDAR derivatives—to produce a single 1:500,000 scale surficial geologic map of Nebraska that is the first of its kind. This new map provides a coarse-scale surficial geologic map that will be incorporated into a nationwide U.S. Geological Survey Quaternary geologic map. It also reflects the variation and uniqueness of physical landscapes in the state, where the Great Plains and Central Lowland physiographic provinces meet, further developing a richer interdisciplinary understanding of regional geomorphology in the heart of North America.

Nebraska↗
USGS2025

Controls on water quality below a reclaimed surface coal mine, southeastern Montana

Coal mining and reclamation can have a profound influence on hydrogeologic systems, with clear consequences for groundwater quality, yet their long-term influence on downgradient water quality over time following reclamation is less well documented. Geochemical trends were evaluated in water quality downgradient of a fully reclaimed landscape at the former Big Sky Mine in the Rosebud Creek watershed (southeastern Montana, USA), over a 3-year period (2020–2022), including bond release in 2022. Within 6 km downgradient from the reclaimed area, sulfate concentrations decreased from approximately 3500 to 1800 mg l −1 within the Miller Coulee alluvial aquifer. Major ions, δ 34 S SO4 values, and residence time tracers suggest that the observed decreases in sulfate concentration result from a combination of dilution by mixed-age inflows and incomplete transit of the high salinity plume from the mine boundary. Both bedrock and alluvial aquifers of the Rosebud Creek corridor contained contributions of millennia-old regional groundwater, which may serve to mitigate mine-derived high salinity waters. Rosebud Creek, which traverses the outflow zone of Miller Coulee in the study area, exhibited high sulfate concentrations during low flows and consistent downgradient increases in sulfate concentration. The possibility of plume dynamics in Miller Coulee suggests that the greatest water quality impacts may not yet have reached Rosebud Creek.

Montana↗
USGS2025

Manipulation of the Symbiodiniaceae microbiome confers multigenerational impacts on symbioses and reproductive ecology of its Exaiptasia diaphana host

Symbiodiniaceae-associated microbiota strongly influence cnidarian symbioses. We systematically reduced the bacterial and fungal communities associated with Symbiodiniaceae to study potential effects on the cnidarian holobiont Exaiptasia diaphana (Aiptasia). Clonal anemones were inoculated with xenic Breviolum minutum (SSB01) and microbiome-manipulated cultures after antibacterial or antifungal treatment. The asexual reproduction of pedal laceration allowed for three generations of clonal aposymbiotic Aiptasia to be utilised in this study, from the initial adult generation (G0), to the first (G1), and second (G2) generation. We inoculated small and large G1 Aiptasia with SSB01 algae and monitored onset of symbiosis, rate of algal proliferation, and holobiont characteristics. Sequencing the 16S and 18S rRNA gene regions identified significant differences in the bacterial and fungal communities of the G0 and G1 generations, alongside differences between the size classes of small and large G1 anemones. The microbiome of larger G1 individuals was distinct to the smaller G1 anemones, suggesting a microbiome maturation process. Control Breviolum minutum cultures exhibited a significantly greater proliferation rate in large G1 anemones when compared to antibacterial or antifungal treated cultures, whereas the opposite trend was documented in the small G1 anemones. Although no differences were observed between algal photochemical parameters, or the growth and polyp activity of G1 juveniles, we observed a significant influence in the production of G2 clones between treatments. Overall, we provide strong ecological implications of manipulating Symbiodiniaceae microbiome, not for the algae themselves, but for the maturation of the host Aiptasia, as well as for the cnidarian holobiont over multiple generations.

The ISME Journal↗
USGS2025

Assessing earthquake risks to lifeline infrastructure systems in the United States

The security and economic stability of the United States rely heavily on robust lifeline infrastructure systems and yet the risks to such systems are seldom quantified at the national scale. For example, while earthquake risks to buildings in the United States have been investigated at the national scale regularly, such risks to gas pipelines have rarely been investigated nationally. In this paper, we use examples from two critical infrastructure sectors to demonstrate (1) the nature of earthquake risks to lifeline infrastructure systems, (2) complexities involved in regional seismic risk assessments, and (3) how such risks change with time. We found that bridge risks can be underestimated by at least 64 % when viewed from repair costs instead of traffic demands and that regional risks can be underestimated by 19 % when spatial correlations of ground motion are ignored. Further, exceedance of traffic demand can be 50 times more likely to occur when viewed at the regional scale than when viewed at an individual bridge. Similarly, exceedance of repairs can be 180 times more likely to occur when viewed at the pipeline network level than at a segment-specific level. Finally, sensitivity analyses with the 2018 and 2023 USGS National Seismic Hazard Models indicate an increase in bridge risk of at least 24 % and an increase in exposed gas pipeline mileage of 43 %. The evolution of risks, complexities involved in assessments, and limited resources jointly underscore the need for more routine updates to nationwide seismic risk assessments of lifeline systems in the United States.

International Journal of Critical Infrastructure P↗
USGS2026

Opportunities for the U.S. Geological Survey’s National Seismic Hazard Model to improve seismic risk assessment of critical infrastructure.

As fragility and risk modeling techniques and computational capabilities evolve, complemented by moving toward more routine and systematic seismic risk assessment of all buildings and critical infrastructure, the authors pose a few critical questions to investigate how the U.S. Geological Survey (USGS) National Seismic Hazard Models (NSHMs) can be used and enhanced further to serve such issues. In this paper, we use three examples from multiple sectors to (1) identify the role of USGS NSHMs in evaluating seismic risks to critical infrastructure, (2) quantify potential impacts from NSHM enhancements (i.e., [i] hazard curves for the vertical component of ground motion, [ii] stochastic event sets, and [iii] maps of probabilistic ground failure hazards), and (3) clarify the feasibility of relevant NSHM improvements. We illuminate that NSHMs are commonly used in location-specific performance assessments, whereas earthquake effects on critical infrastructure can be widespread across large geospatial regions. Further, we found that without the NSHM extensions considered here, risk can be severely underestimated, e.g., neglecting ground failure hazards can underestimate regional loss by a factor of two or more. Although many challenges remain, we developed example prototypes to clarify the feasibility of the NSHM extensions, which can facilitate improved management of risks to critical infrastructure.

Earthquake Spectra Journal↗
USGS2026

Surface-wave relocation and characterization of the 2024–2025 dyking episode along the Fentale–Dofen segment of the Ethiopian rift

From late-December 2024 to mid-March 2025, a 50-km-long dyke intrusion triggered over 300 earthquakes (magnitude 4 to 5.9) between Fentale and Dofen volcanoes along the Northern Main Ethiopian Rift. Dyke intrusions periodically occur along the Fentale–Dofen magmatic segment and are an expression of ongoing rift extension. Preliminary analyses using interferometric synthetic aperture radar revealed extensive ground deformation (up to 60 cm), which closely matched the temporal and spatial evolution of surface manifestations and earthquake locations from global catalogues. While global catalogues are critical for real-time monitoring, the precision of locations in remote and or sparsely instrumented regions can be low. In this investigation, we present surface-wave relocation results of the dyking episode that began near Fentale volcano in December 2024. We estimate relative locations using differential traveltimes measured from regional-to-teleseismic distance surface-wave observations of earthquakes reported by the U.S. Geological Survey. Relative relocations reduce the initial region of diffuse seismicity to a 50-km-long narrow band bounding the strike of surface manifestations and the zone of maximum surface deformation. We demonstrate the precision of surface-wave relocations over incremental time periods, capturing the progression of dyking from seismic onset through seismic migration and caldera subsidence. Results showcase the utility of surface-wave relocations in the characterization of dyking episodes and provide complementary insights into the current understanding of the Fentale–Dofen volcanic plumbing system.

Geophysical Journal International↗
USGS2026

Ringtail (Bassariscus astutus) survival in southwestern Oregon

The effective conservation and management of small carnivore populations requires understanding species’ life-history traits and identifying important vital rates that drive population trajectories. However, many of these species are rare or elusive and of state or federal conservation concern, and demographic information is often lacking and difficult to obtain. At the northern limit of their range in Oregon, USA, ringtail ( Bassariscus astutus ) occupy mid-elevation forests and are a species of conservation concern because of their limited distribution and suspected low density. We initiated a radio-telemetry study in 2020 to estimate monthly and annual survival of ringtail in southwest Oregon. We monitored 26 ringtail from November 2020 to October 2022 and estimated survival rates using a known-fate framework and Program MARK. Model-averaged monthly survival estimates ranged from lows of 0.963 (SE = 0.022, 95% CI = 0.887–0.988) in spring (Feb–May) of 2021 to highs of 0.980 (SE = 0.018, 95% CI = 0.889–0.997) in summer (Jun–Oct) of 2022. Model-averaged estimates of annual survival were 0.695 (SE = 0.176, 95% CI = 0.310–0.920) during 2020–2021 and 0.728 (SE = 0.168, 95% CI = 0.336–0.934) during 2021–2022. Predation was the leading cause of mortality, but notably, no mortalities could be attributed to avian predation. Survival rates were much higher within forested landscapes in Oregon than for populations in the southwestern United States, which may have implications for regional conservation and management strategies.

Oregon↗
USGS2026

A review and synthesis of post-wildfire shifts in hydrologic processes and streamflow generation mechanisms

Critical water supply watersheds in the western United States (WUS) are impacted by wildfires, with potential negative effects on water quality and quantity. Scientific understanding is currently insufficient to deliver estimates of wildfire consequences for water quantity that are regionally accurate. Regional variability in the directionality and magnitude of post-wildfire shifts in streamflow generation fuels uncertainty in estimates of wildfire effects on water supply. In this work we provide a narrative review of wildfire effects on hydrologic processes and the resulting changes in streamflow generation mechanisms with a focus on the WUS, incorporating other global regions when pertinent. A conceptual model summary of wildfire effects on streamflow generation emphasizes: (1) precipitation seasonality, (2) synchrony of precipitation and potential evapotranspiration, (3) net shifts in interception, evaporation, and transpiration relative to total annual precipitation, (4) vegetation changes, including compensatory uptake and type conversion, (5) degree of overlap in rainfall rates and infiltration, (6) fire extent and severity, (7) burn scar positioning (e.g. in headwaters or proximal to watershed outlet), (8) scale-dependent groundwater leakage, (9) near-surface water storage reduction, and (10) soil to groundwater connectivity. Ongoing gaps and challenges include separating the influences of precipitation variability, water withdrawals, and post-fire land management; compound and overlapping disturbances; and lack of pre-fire data. Notable future opportunities include: harnessing ever-improving gridded and remotely sensed precipitation and fire-effects data; linking geophysical, isotopic tracer, and geochemical signatures to diagnose hydrologic changes; leveraging physically based and data-driven model advancements; and analyzing streamflow generation recovery trajectories across diverse watersheds.

western United States↗
USGS2025

Water supply in the conterminous United States, Alaska, Hawaii, and Puerto Rico, water years 2010–20

We present an assessment of water supply across the conterminous United States (CONUS), Alaska, Hawaii, and Puerto Rico covering water years 2010–20. Our analysis drew on two national hydrologic models, the National Hydrologic Model Precipitation-Runoff Modeling System and the Weather Research and Forecasting model hydrologic modeling system. Both models produced estimates of streamflow, evapotranspiration, soil moisture, snow water equivalent, and other hydrologic states and fluxes. The models were driven by the bias-adjusted 4-kilometer-resolution, long-term regional hydroclimate simulation over the conterminous United States dataset (CONUS404). We assessed spatial and temporal error distributions by comparing monthly simulations at the 12-digit hydrologic unit code and regional scale from both models against external benchmarking datasets. Results showed that average annual rainfall across the CONUS was 857 millimeters per year for the period of analysis, with water year 2012 the driest year (729 millimeters) and water year 2019 the wettest year (995 millimeters). Key interannual variability results included the following: (1) the California–Nevada hydrologic region had the highest variability in precipitation and snow accumulation, and (2) the Texas hydrologic region was among hydrologic regions with the highest variability in precipitation. We related interannual variability in precipitation to storage volumes in soil moisture, snow water equivalent, and lakes and reservoirs to highlight areas with little storage and large year-to-year variability in precipitation. These areas included the Southern High Plains, Central High Plains, Texas, Souris–Red–Rainy, Mississippi Embayment, and Midwest regions. Our analysis of groundwater-level data showed that several of these areas overlap aquifers where groundwater levels were considerably lower than historical averages, including the Colorado Plateaus aquifers, the Rio Grande aquifer system, and the Central and Southern regions of the High Plains aquifer. Many of these lowered groundwater levels are continuations of decades-long declines from overpumping that started well before the assessment period. The resulting water budgets and their analyses provide a high-resolution foundational assessment of the mean state and variability of the terrestrial hydrologic cycle across the CONUS and Alaska, Hawaii, and Puerto Rico to support a wide range of water resource management applications.

Professional Paper↗
USGS2025

Central Valley Hydrologic Model version 2 (CVHM2): Decision support tool for groundwater and land subsidence management

The San Joaquin Valley (SJV) of California is one of the world’s most productive agricultural regions. Reliance on groundwater has led to some of the greatest rates of human-induced land subsidence in the world in the 20th century, as well as more recently. The United States Geological Survey (USGS) has recently developed an integrated surface–subsurface hydrologic model, the Central Valley Hydrologic Model 2 (CVHM2), that represents the major components of the hydrologic system of California’s Central Valley. In this study, CVHM2 was applied as a decision support tool while simulating various management strategies to mitigate the land subsidence caused by the extraction of groundwater. CVHM2 was extended through to 2073 and applied to simulate management scenarios in terms of three primary drivers and their impact on subsidence along the Delta–Mendota Canal (DMC), a critical piece of infrastructure in the western SJV. The drivers considered were agricultural water demands, managed aquifer recharge (MAR), and changes in future climate. The results show that future subsidence is most sensitive to water demands, second most sensitive to future changes in climate, and relatively insensitive to MAR when it is applied as a surface application in the western SJV. However, we demonstrate via proof-of-concept scenarios that the MAR is capable of arresting subsidence when implemented via injection below the Corcoran Clay Member of the Tulare Formation instead of as a surface application. We also examine the uncertainty that is the result of climate variability and how to use the tool to identify the most appropriate strategies to constrain future subsidence to acceptable levels.

California↗
USGS2025

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

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

Planetary Science Journal↗
USGS2026

Northern bobwhites select for native grasses on working grazing land

Northern bobwhite ( Colinus virginianus ; bobwhite) populations have experienced an 85% decline across most of its range since the 1960s. The most drastic decreases have been in the southeastern United States where biologists attribute the decline to widespread habitat loss, including conversion of native grasslands to tall fescue ( Lolium arundinaceum ) and other exotic cool-season grass (CSG) pasture and hayfields. Including agricultural lands in conservation programs could improve habitat conditions on a regional scale. One working-lands conservation strategy involves the use of native warm-season grasses (NWSGs) rather than cool-season exotic grasses as economically viable cattle forage. To evaluate this management practice for creating bobwhite habitat, we conducted a field study on a bobwhite population in eastern Kentucky, USA, where grazed NWSG and burned NWSG fields were adjacent to exotic grazed and exotic hayed CSG fields. Between April 2019 and October 2022, we used radio telemetry to track bobwhites to evaluate resource selection at second- and third-order scales. Bobwhites used traditionally managed CSGs less than would be expected by chance at the second-order scale in both breeding and non-breeding seasons. Bobwhite use of grazed NWSG and burned NWSG was greater than would be expected by chance at the second-order scale in both breeding and non-breeding seasons. Similarly, at the third-order scale, bobwhites used CSGs less than would be expected by chance, whereas NWSGs and woody stems were used more than by chance alone. Our study suggests that bobwhites select NWSGs at multiple spatial scales despite broad-scale availability of managed exotic CSGs. Under a working-lands framework, the integration of NWSGs into working lands could create bobwhite habitat in the southeastern United States, especially if woody stems are present.

Kentucky↗
USGS2026

Continental-scale prediction of hydrologic signatures and processes

Understanding how dominant hydrologic processes and their drivers vary across diverse continental-scale landscapes is critical for hydrologic modeling and water management applications. Our research addresses this question by synthesizing large-sample watershed datasets, Caravan and GAGES-II, and developing random forest models to identify patterns in hydrologic function. We assessed dominant processes by examining hydrologic signatures – summary indicators of watershed function derived from hydroclimatic time series and random forest models across 14 146 gauged United States watersheds. The results reveal clear continental-scale gradients in hydrologic processes, including baseflow, overland flow, storage, and water balance losses. Our map of dominant processes highlights, for example, the transition from baseflow to fast responses and back to baseflow along the elevation gradient from the Appalachian spine, through the Piedmont, to the Eastern Coastal Plain; a distinct outer ring around the Great Lakes region; and sharp contrasts between coastal and inland processes in the West. Variable importance analysis from random forest models show that processes in the western U.S. are primarily controlled by climate, whereas in the eastern U.S., soil, geology, and topography play larger roles, with distinct human influences apparent in urban areas. Our approach of estimating dominant processes and their drivers facilitates extending process knowledge from research watersheds to the continental scale, assessing current hydrological understanding, and evaluating hydrological model structures.

conterminous United States↗
USGS2026

Competition between invasive round goby and native slimy sculpin in a coldwater laboratory setting

Slimy sculpin ( Cottus cognatus ) abundance has substantially declined in many regions of the Laurentian Great Lakes coincident with the arrival and range expansions of invasive round goby ( Neogobius melanostomus ). Previous studies have shown that round goby can outcompete native mottled sculpin ( C. bairdii ) for shelter and negatively affect native spoonhead sculpin ( C. ricei ) growth, but direct behavioral interactions of slimy sculpin and round goby have not been observed. We conducted a laboratory experiment in which a pair of slimy sculpin were placed in an aquarium that contained an artificial spawning shelter and one round goby (treatment, n = 10) or no other fish (control, n = 10) for at least 30 days at temperatures reflective of their potential offshore overlap. We measured slimy sculpin egg deposition, aggressive interactions, artificial shelter occupancy, growth, and survival. Slimy sculpin were less likely to occupy shelters and experienced increased intraspecific competition when a round goby was present. In contrast, round goby rarely chased slimy sculpin, slimy sculpin were more likely to chase round goby, and slimy sculpin growth, egg deposition, and survival were not significantly different between groups. Differences in density between treatments and the relative size of fish could confound our observations. However, the dominant slimy sculpin occupied the shelter at similar rates between experimental groups if the round goby was smaller. Our results indicate that larger round goby overlapping with slimy sculpin could affect slimy sculpin shelter and space use, but the Great Lakes environment may create more stress for slimy sculpin than we observed.

Lake Michigan, Lake Superior↗
USGS2026

Thermal detectability of subsurface water ice on Mars: A comparative analysis for the Subsurface Water Ice Mapping (SWIM) Project

We have developed a new global map of shallow ground ice distribution, SWIM23, based on Mars Global Surveyor Thermal Emission Spectrometer data and made systematic comparisons between this new map and two similar, previously developed data products. We have explored the origin of differences between the three ice maps by detailing technical and procedural differences in their development, by making global pixel-by-pixel comparisons, and by carrying out a series of one-dimensional thermal simulations to explore fundamental physical limitations of thermal ice-detection techniques. These efforts and the production of a composite thermal ice-consistency map supported integration of multiple geophysical data products relevant to ice detection in the upper meter of the Martian regolith into the larger Mars Subsurface Water Ice Mapping project. Our work also highlights fundamental physical limitations to thermal ice detection as a technique, particularly the rapid fall-off in ice detection sensitivity at depths >30 cm, which produces maximum uncertainty in the presence and depth of ice within regions preferred for potential human landing sites. A future Mars orbiter mission designed to detect ice and support crewed landing site selection in the midlatitude region should give payload priority to an instrument capable of probing the 1–5 m depth range (i.e., a high-frequency radar), over a next-generation thermal spectrometer, which is unlikely to offer clarity on ice table depths or lateral continuity of the ice table in the locations of highest interest.

Planetary Science Journal↗
USGS2026

Sedimentological and geochemical characterization of lacustrine deposits of the Babouri-Figuil basin, northern Cameroon: Implications for source rocks distribution and petroleum exploration

The West and Central African Rift System (WCARS) refers to the series of Cretaceous rift basins where commercial hydrocarbon accumulations have been discovered. Some of the WCARS frontier basins are currently being investigated to increase our understanding of these basins in light of new commercial discoveries. The present study was performed in the Babouri-Figuil Basin (BFB), which is genetically related to the WCARS and constitutes an area of interest in terms of petroleum prospecting, where the distribution of petroleum source rocks and potential targets for petroleum exploration across the entire basin is poorly understood. For the current study, an integrated facies analysis along with organic and inorganic geochemical techniques were applied to the basin's Cretaceous deposits with the aim of reconstructing the paleodepositional environment, assessing factors that triggered the input of organic matter, and providing a spatial overview of the organic matter accumulation in the basin based on outcrop samples. An alluvial fan-lacustrine-braided river system is inferred from the facies analysis of the stratigraphic sequence consisting of conglomerate, sandstone, siltstone, limestone, marlstone, and claystone. Bulk analysis of organic matter reveals that black shale and massive claystone are the main prospective petroleum source rocks in the basin. Inorganic geochemical analyses reveal the influence of anoxic conditions, moderate to high primary productivity, and low terrigenous inputs in organic matter enrichment. The formations rich in organic matter are predominantly concentrated in the western and eastern parts of the basin which may represent areas with depressions, characterized by high accommodation space. In terms of the regional context of the WCARS rift basins, typical hydrocarbon exploration in the BFB may target basal-conglomerate, sandstone beds situated directly above or/and interbedded with the Lower Cretaceous source rocks, and the Upper Cretaceous sandstone beds. Basement rocks (granite, granodiorite, and gneisses) and oil shale deposits may represent potential unconventional hydrocarbon exploration. The current integrated study provides an insight that should guide future hydrocarbon exploration campaigns in the basin.

Babouri-Figuil basin↗
USGS2026

Subsurface water ice mapping on Mars: A probabilistic approach

Subsurface water ice deposits on Mars are an important resource for potential future human exploration. They are also an indicator of the planet’s past climate. However, the distribution of subsurface water ice in Mars’s midlatitudes is uncertain because spacecraft imagery cannot directly observe subsurface ice in most cases. Various spacecraft remote sensing instruments are sensitive to subsurface water ice, including thermal imaging spectrometers, radar sounders, and neutron spectrometers. Geomorphic analyses of images can also implicate subsurface ice. Building upon the data products from the Mars Subsurface Water Ice Mapping project, we provide a probabilistic framework to jointly interpret existing data and estimate the likelihood of subsurface water ice in the Martian midlatitudes between 60 ∘ S and 60 ∘ N with uncertainty. Broadly, we find that near-surface ice is likely present poleward of ∼45 ∘ in both the northern and southern hemispheres. However, closer to the equator, existing remote sensing data cannot uniquely constrain the presence of subsurface water ice. Our probabilistic results provide a framework for quantifying the abundance of ice on Mars, and our uncertainty estimates allow future analysis and exploration to target regions of high uncertainty.

Planetary Science Journal↗
USGS2025

Identification of representative earthquakes for probabilistic tsunami hazard analysis (PTHA) using earthquake rupture forecasts and machine learning

As probabilistic tsunami hazard analysis (PTHA) focuses more on assessments for localized, populous regions, techniques are needed to identify a subsample of representative earthquake ruptures to make the computational requirements for producing high-resolution hazard maps tractable. Moreover, the greatest epistemic uncertainty in seismic PTHA is related to source characterization, which is often poorly defined and subjective. We address these two salient issues by applying streamlined earthquake rupture forecasts (ERFs), based on combinatorial optimization methods, to an unsupervised machine learning workflow for identifying representative ruptures. ERFs determine the optimal distribution of a millennia-scale sample of earthquakes by inverting the observed slip rate on major faults. We use two previously developed combinatorial optimization ERFs, integer programming and greedy sequential, to produce the optimal location of ruptures with seismic moments sampled from a regional Gutenberg–Richter magnitude–frequency distribution. These ruptures in turn are used to calculate peak nearshore tsunami amplitude, using computationally efficient tsunami Green's functions. An unsupervised machine learning workflow is then used to identify a small subsample of the earthquakes input to ERFs for onshore PTHA analysis. We eliminate epistemic uncertainty related to source distribution under traditional PTHA analysis; in its place, a quantifiable, less subjective and generally smaller uncertainty related to the input to ERFs is included. The Nankai subduction zone is used as a test case, where previous ERFs have been conducted. Results indicate that the locations of representative earthquakes are sensitive to choice of magnitude–area relation and to whether a minimum cumulative stress objective is imposed on the fault. In general, incorporating ERFs into PTHA provide a physically self-consistent method to incorporate fault slip information in determining representative earthquakes for onshore PTHA, eliminating a major source of epistemic uncertainty.

Nankai subduction zone↗