USGS Science⌕ Search

SEARCH · USGS Science

Results for “Bulletin, Southern California Academy of Sciences”

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

Airborne geophysical analysis to decipher salinization for coastal Louisiana

Coastal Louisiana is known for saltwater intrusion that threatens wetlands, aquifers, and rivers. However, the extent of saltwater intrusion is not well understood. This study develops an innovative framework with airborne electromagnetic (AEM) data to map chloride concentration distributions for wetlands in the Mississippi River deltaic plain and Chenier plain as well as for the Mississippi River Valley alluvial aquifer (MRVA) and Chicot aquifer. Moreover, the framework maps chloride concentrations along the Mississippi River and Atchafalaya River. Key components in the framework include the establishment of resistivity-to-chloride concentration transformation, 3D resistivity architecture building through geostatistics, and the employment of a lithologic model. The transformation functions correlate AEM resistivity data with porewater salinity measurements and groundwater and river chloride samples. The results show that AEM data reliably infers soil water chloride concentrations and correlates well with the distribution of various marsh types. AEM data reveals extensive saltwater presence at depth and near the coast, originating from salt domes and the Gulf of Mexico, respectively. The saltwater upconing pattern in the Chicot aquifer is likely due to excessive groundwater withdrawals. The AEM data also confirms a distinct tongue of saltwater intruding into the Atchafalaya Basin from the Gulf. The AEM data helps to identify faults that are obscured or eroded at the surface, which appear as leaky barriers in the subsurface where dramatic changes in chloride concentration are apparent. Finally, this study uses the AEM data to infer the presence of an extensive seawater wedge in the Mississippi River and Atchafalaya River.

Louisiana↗
USGS2024

Perspectives on transportable array Alaska background noise levels

Background seismic noise fundamentally sets a lower bound on our ability to record signals arising from earthquakes. The background noise spectrum at a station is a combination of cultural noise, ocean-generated microseism noise, intrinsic instrument self-noise, and the sensitivity of the instrument to nonseismic noise sources. The USArray-Transportable Array Alaska deployed 195 stations across Alaska and parts of Canada (Yukon, British Columbia, and Northwest Territories). These stations were all installed using similar techniques and made use of instruments with similar self-noise levels. As such, this network provides an opportunity to look at how geographic location influences seismic background. Using these broadband stations, we report background noise levels from 0.2 to 75 s period in six discrete bands. By constructing “noise maps,” we depict both spatial and temporal changes in the background noise field. Using these maps, combined with targeted analysis, we infer sources and contributing factors to noise levels in these different period bands. These include cultural noise, the formation of sea ice, seasonal changes in permafrost and wave activity in the Gulf of Alaska, and magnetic field variability. We use this study as an opportunity to review several previous studies examining seismic noise in Arctic regions.

Book chapter↗
USGS2025

Groundwater structures fish growth and production across a riverscape

Landscapes are composed of habitat patches and conditions that vary across space and time. While habitat variability and complexity can support important ecological processes and ecosystem services, the dynamic nature of habitats can also constrain organismal growth and production as optimal conditions are fleeting. In riverine ecosystems, groundwater discharge to streams stabilises water temperature and flow regimes, thus mediating how habitat complexity is expressed. Yet, how stable habitats structure growth and production within the broader landscape matrix is not well understood. In this study, we explored the effects of groundwater on spatiotemporal variation in growth and production for juvenile Yellowstone cutthroat trout ( Oncorhynchus virginalis bouvieri ) across the upper Snake River catchment, Wyoming, USA. We combined machine learning techniques and remotely sensed landscape data to estimate groundwater availability across the river network, which we linked to stream temperature regimes and conspecific density. We then used Bayesian hierarchical models to quantify the effects of temperature, density and groundwater on spatiotemporal variation in fish growth and production in 52 focal reaches. Finally, we predicted body size trajectories and trends in total production continuously over both space and time to understand the effect of groundwater at the riverscape scale. Groundwater discharged to streams where topography changes abruptly in valley-bottom areas underlain by coarse glacial deposits. Groundwater stabilised temperature regimes and was associated with high trout densities. Temperature and density, in turn, interacted to influence growth rates: growth increased strongly with temperature, but this effect was reduced when density was high. Accordingly, variation in groundwater availability among stream reaches diversified growth and production regimes. In reaches with low groundwater availability, growth and production declined over time from summer maxima. In contrast, in reaches with high groundwater availability, temporal trends in growth and production were hump-shaped—peaking in autumn—and mean production was greater. At the riverscape scale, temporal asynchrony in growth rates generated convergent spatial variation in growth capacity, but—when combined with density—led to the formation of distinct hotspots of production. Our results demonstrate how groundwater, an important driver of aquatic ecosystem heterogeneity, structures trout growth and production across space and time. Importantly, rare, but stable habitats may disproportionately affect ecological processes and serve as key sources of population diversity at larger spatial scales.

Wyoming↗
USGS2025

An unexplained tsunami: Was there megathrust slip during the 2020 Mw7.6 Sand Point, Alaska, earthquake?

On October 19, 2020, the M w 7.6 Sand Point earthquake struck south of the Shumagin Islands in Alaska. Moment tensors indicate the earthquake was primarily strike-slip, yet the event produced an enigmatic tsunami that was larger and more widespread than expected for an earthquake of that magnitude and mechanism. Using a suite of hydrodynamic, seismic, and geodetic modeling techniques, we explore plausible causes of the tsunami. We find that strike-slip models consistent with the moment tensor orientation cannot produce the observed tsunami. Hydrodynamic inversion of sea surface deformation from deep ocean and tide gauge data suggest seafloor deformation more closely matches a megathrust, rather than a strike-slip, source. Static slip inversions, using sea level and Global Navigation Satellite System data, allow for a portion of co-seismic megathrust slip that can explain tsunamigenesis. Combining all available geophysical datasets to model the kinematic rupture, we show that considerable, relatively slow, megathrust slip is allowable in the Shumagin segment, concurrent with strike-slip faulting. We hypothesize that the slow megathrust rupture does not contribute much seismic radiation allowing it to previously go unnoticed with traditional seismic monitoring.

Alaska↗
USGS2025

Hydrogeology and groundwater quality in the Snake River alluvial aquifer at Jackson Hole Airport, Wyoming, 2011–20

The Snake River alluvial aquifer underlying the Jackson Hole Airport (JHA) in northwest Wyoming is an important source of water used for domestic, commercial, and irrigation purposes by the airport and nearby residents. The U.S. Geological Survey, in response to previously identified water-quality concerns in the area, monitored and evaluated changes in hydrogeologic characteristics and groundwater-quality conditions of the alluvial aquifer during 2011–20. During that period, the Jackson Hole Airport made several changes that potentially improved water quality at and downgradient from the airport. Well, water level, and hydrogeologic data were collected from the alluvial aquifer to identify hydrogeologic characteristic and groundwater quality changes. Additionally, results of statistical tests were applied to water-quality results to evaluate trends in selected physical properties and constituent concentrations with time. The trends of those data show that water quality did improve overall during the study period compared to previously collected data. Presumably, these trends are in response to the changes in the aircraft deicing/anti-icing fluid (ADAF) formulation used by the JHA, the many JHA infrastructure improvements made during 2011–20, the degradation of existing ADAFs in subsurface soils and groundwater, or some combination of these possibilities.

Wyoming↗
USGS2025

U.S. Geological Survey National Groundwater Climate Response Network

What is the U.S. Geological Survey National Groundwater Climate Response Network? The U.S. Geological Survey’s (USGS) National Water Monitoring Network is a network of networks that includes the Climate Response Network (CRN). The CRN is a network of wells selected to monitor the effects of climate variability, including droughts, on groundwater systems within the United States, Puerto Rico, and the U.S. Virgin Islands. The primary purpose of the CRN is to monitor these effects on groundwater levels in unconfined or near-surface confined aquifers that are minimally affected by pumping or other human-induced stresses.

Fact Sheet↗
USGS2024

Stream nitrate dynamics driven primarily by discharge and watershed physical and soil characteristics at intensively monitored sites: Insights from deep learning

We developed a suite of models using deep learning to make hindcast predictions of the 7‐day average backward‐looking nitrate concentration at 46 predominantly agricultural sites across the midwestern and eastern United States. The models used daily observations of discharge and meteorological variables and watershed attributes describing anthropogenic modification to hydrology, nitrogen application, climate, groundwater, land use, watershed physiographic attributes, and soils. Across all sites, discharge and watershed soil and physiographic attributes showed a strong influence on model performance. Analysis of drivers across sites revealed considerable regional differences related to controlling processes such as groundwater contributions. We tested several ways to pool data across sites to develop accurate models and make the most effective use of available data. Single‐site models, in which models are trained and tested at a single location, showed generally strong predictive performance (median Kling‐Gupta Efficiency = 0.66), and accuracy at poorly performing sites could be improved by grouping sites with similar characteristics. Developing a single model for all sites reduced performance at several locations with distinct characteristics, suggesting that there is a threshold of dissimilarity beyond which more data does not improve the model. While many deep learning studies have shown that national or even global models can outperform local models, it is not clear that this is true for water quality constituents. This study demonstrates how data can be combined effectively, using deep learning to develop accurate and interpretable models of instream nitrate at sites where varying processes are responsible for changes in nitrate concentration.

Water Resources Research↗
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↗
USGS2024

Improving crop-specific groundwater use estimation in the Mississippi Alluvial Plain: Implications for integrated remote sensing and machine learning approaches in data-scarce regions

Study region The Mississippi Alluvial Plain (MAP) in the United States (US). Study focus Understanding local-scale groundwater use, a critical component of the water budget, is necessary for implementing sustainable water management practices. The MAP is one of the most productive agricultural regions in the US and extracts more than 11 km 3 /year for irrigation activities. Consequently, groundwater-level declines in the MAP region pose a substantial challenge to water sustainability, and hence, we need reliable groundwater pumping monitoring solutions to manage this resource appropriately. New hydrological insights for the region We incorporate remote sensing datasets and machine learning to improve an existing lookup table-based model of groundwater use previously developed by the U.S. Geological Survey (USGS). Here, we employ Distributed Random Forests, an ensemble machine learning algorithm to predict annual and monthly groundwater use (2014–2020) throughout this region at 1-km resolution, using pumping data from existing flowmeters in the Mississippi Delta. Our model compares favorably with the existing USGS model, with higher R 2 (0.51 compared to 0.42 in the previous model), and lower root mean square error (RMSE) and mean absolute error (MAE)— 0.14 m and 0.09 m, respectively in our model, compared to 0.15 m and 0.1 m in the previous model. Therefore, this work advances our ability to predict groundwater use in regions with scarce or limited in-situ groundwater withdrawal data availability.

Journal of Hydrology Regional Studies↗
USGS2025

Light absorbing particles deposited to snow cover across the Upper Colorado River Basin, Colorado Rocky Mountains, 2013-16: Interannual variations from multiple natural and anthropogenic sources

Atmospheric particulate matter (PM) as light-absorbing particles (LAPs) deposited to snow cover can result in early onset and rapid snow melting, challenging management of downstream water resources. We identified LAPs in 38 snow samples (water years 2013–2016) from the mountainous Upper Colorado River basin by comparing among laboratory-measured spectral reflectance, chemical, physical, and magnetic properties. Dust sample reflectance, averaged over the wavelength range of 0.35–2.50 μm, varied by a factor of 1.9 (range, 0.2300–0.4444) and was suppressed mainly by three components: (a) carbonaceous matter measured as total organic carbon (1.6–22.5 wt. %) including inferred black carbon, natural organic matter, and carbon-based synthetic, black road-tire-wear particles, (b) dark rock and mineral particles, indicated by amounts of magnetite (0.11–0.37 wt. %) as their proxy, and (c) ferric oxide minerals identified by reflectance spectroscopy and magnetic properties. Fundamental compositional differences were associated with different iron oxide groups defined by dominant hematite, goethite, or magnetite. These differences in iron oxide mineralogy are attributed to temporally varying source-area contributions implying strong interannual changes in regional source behavior, dust-storm frequency, and (or) transport tracks. Observations of dust-storm activity in the western U.S. and particle-size averages for all samples (median, 25 μm) indicated that regional dust from deserts dominated mineral-dust masses. Fugitive contaminants, nevertheless, contributed important amounts of LAPs from many types of anthropogenic sources.

Colorado↗
USGS2025

Flow heterogeneity controls dissolution dynamics in topologically complex rocks

Rock dissolution is a common subsurface geochemical reaction affecting pore space properties, crucial for reservoir stimulation, carbon storage, and geothermal energy. Predictive models for dissolution remain limited due to incomplete understanding of the mechanisms involved. We examine the influence of flow, transport, and reaction regimes on mineral dissolution using 29 time-resolved data from 3D rocks. We find that initial pore structure significantly influences the dissolution pattern, with reaction rates up to two orders of magnitude lower than batch conditions, given solute and fluid-solid boundary constraints. Flow unevenness determines the location and rate of dissolution. We propose two models describing expected dissolution patterns and effective reaction rates based on dimensionless metrics for flow, transport, and reaction. Finally, we analyze feedback between evolving flow and pore structure to understand conditions that regulate/reinforce dissolution hotspots. Our findings underscore the major impact of flow arrangement on reaction-front propagation and provide a foundation for controlling dissolution hotspots.

Geophysical Research Letters↗
USGS2025

Dynamic feedbacks between river meandering and landsliding in northwestern Washington glacial terraces

Landsliding in river valleys poses unique risks for cascading hazards and can damage infrastructure and cause fatalities. In postglacial valleys, many landslides are posited to occur in relation to lateral river erosion, but the dynamics of fluvial-hillslope interactions are not well understood. Here, we investigate a section of the Nooksack River in western Washington State where the channel is flanked by landslide-prone glacial terraces similar to those that failed in the 2014 State Route 530 “Oso” landslide. We map 216 landslides through time across 17 aerial imagery data sets (1933–2022) and analyze them in relation to river meandering and curvature. We observe dynamic feedbacks between lateral river meandering and valley-adjacent landsliding. Terrace lateral retreat rates of up to 25 m/year owing to combined fluvial erosion and slope failure occur on pinned, outer meander bends immediately downstream from peaks in river curvature (>0.0075 1/m); these locations are predisposed to both shallow and deep-seated landslides. Deep-seated landslides extending 17%–32% of the active valley width into the floodplain can displace the river away from the floodplain margin and change the channel planform. River-displacing landslides relocate meanders up- or downstream, thereby conditioning the location of subsequent landslides. This conceptual model of coupled landslide-driven meander displacement and valley-adjacent landsliding is exemplified across western Washington river systems. The distance between up- and downstream valley-adjacent landsliding scales with valley width, meander wavelength, and terrace height. Our results can advance our understanding of the river-hillslope interface in landscape evolution and can be used to inform hazard management in river corridors.

Washington↗
USGS2026

Magnitude conversion relations create substantial differences in seismic hazard models

Earthquake catalogs are essential data inputs for seismic hazard modeling. Because earthquake magnitudes are reported in a variety of types (e.g., local magnitudes and moment magnitudes), magnitude conversion relationships must be used to convert the different magnitude types present in a catalog to a uniform magnitude type to avoid biases in the hazard computation. However, these conversion relationships are often uncertain and have been shown to sometimes perform poorly. Here, we investigate the sensitivity of the gridded seismicity component of the National Seismic Hazard Model (NSHM) to the catalog conversion equations in the Eastern United States. In the 2023 NSHM, magnitudes of various types were converted to moment magnitudes using equations developed by the Central and Eastern United States Seismic Source Characterization for Nuclear Facilities (CEUS‐SSCn), based on least‐squares (LS) regressions made using data from a catalog containing events up through 2008. We recompute these equations using events in the Advanced National Seismic System Comprehensive Earthquake Catalog with multiple magnitudes from 2000 to 2023. Although we prefer the use of orthogonal regressions for our datasets, LS regressions produce broadly similar results, with both approaches exhibiting large deviations from the CEUS‐SSCn conversions, especially at smaller magnitudes. We compare the spatial distribution of annual rates using three different models: (1) the 2023 NSHM conversions, (2) our updated conversions, and (3) no conversions. We find that the choice of conversions leads to substantial differences in the rate forecasts, which can greatly impact the seismic hazard model, particularly in regions with low‐seismicity rates such as the Eastern United States, where the hazard is dominated by gridded seismicity rather than a fault model.

Seismological Research Letters↗
USGS2026

Channel morphology and large wood control postfire debris-flow erosion and deposition

Runoff-generated debris flows are a known response to wildfire, and accurately predicting the volume of these debris flows is important for estimating the magnitude of downstream hazards. Prior data collection efforts have focused on debris-flow volume measurements at catchment outlets, but few studies have considered how erosion and deposition modulate the volume of sediment arriving at catchment outlets. This study takes advantage of a high-resolution dataset to document the factors that control the total debris-flow volume reaching the catchment outlet during a fatal postfire debris flow. Using pre- and post-event airborne lidar, satellite imagery and field mapping, we found that a postfire debris flow in the Black Hollow catchment in northern Colorado eroded 136,000 ± 30,000 m 3 and redeposited 27,000 ± 7,500 m 3 in the main channel. Most of the in-channel deposition (52% by volume) occurred where a confined channel reach transitioned to an unconfined channel reach downstream, allowing the flow to widen and deposit material. Wood jams played multiple roles in the debris-flow dynamics, both nucleating deposition (25% of the deposit volume was stored behind wood jams) and exacerbating erosion (50% of the total erosion occurred downstream from a wood dam break). The remaining deposition occurred due to spatial changes in channel slope as well as deposition observed at newly formed channel bars. Using these data in this study, we identified topographic and vegetation metrics that can be used (pre-event) to anticipate where deposition may occur in channels prior to a debris flow.

Colorado↗
USGS2025

Estimating agricultural irrigation water consumption for the High Plains aquifer region with integrated energy- and water-balance evapotranspiration modeling approaches

Estimation of irrigation water use provides essential information for the management and conservation of agricultural water resources. Conventionally, water use data are created based on reports and surveys from water users, whereas manual records may not be complete due to lacking flow meters, measurement gaps, inconsistent methods across regions, and time- and cost-consuming data processing. Alternatively, spatially explicit estimation of irrigation water use can be conducted efficiently using remote sensing evapotranspiration (ET) modeling approaches. In this study, we created a gridded blue water evapotranspiration (BWET) dataset to estimate historical irrigation water consumption (1986 – 2020) in the croplands across the United States High Plains aquifer region. The BWET data were generated by integrating an energy-balance ET model [Operational Simplified Surface Energy Balance model (SSEBop)] and a water-balance ET model [Vegetation ET model (VegET)]. BWET in croplands indicates crop consumptive use of irrigation water extracted from surface water and groundwater resources. The BWET estimates were compared with reported irrigation water use data for all counties within the aquifer region. The results revealed high agreement between growing season (May – September) BWET and annual water withdrawal at county level. Specifically, correlation coefficients of volumetric BWET and water withdrawal were 0.90 and 0.96, respectively, for the entire aquifer region and western Kansas. The timeseries of BWET and water withdrawal showed similar temporal trends and high covariations. The BWET estimates were systematically lower than the water withdrawal measurements, which was primarily attributed to blue water losses in the irrigation system. The irrigation efficiency, calculated as the ratio of BWET to water withdrawal depth, was 0.57 and 0.74 for the entire aquifer region and western Kansas, respectively. This study demonstrates the capability of using satellite-based ET models (e.g., SSEBop and VegET) to efficiently estimate crop water consumption and evaluate irrigation efficiency at landscape, county, and regional scales.

Colorado, Kansas, Nebraska, New Mexico, Oklahoma, ↗
USGS2025

Hydraulic conductivity and transmissivity estimates from slug tests in wells within the Mississippi Alluvial Plain, Arkansas and Mississippi, 2020

During the spring and summer of 2020, the U.S. Geological Survey conducted single-well slug tests on selected observation wells within the Mississippi Alluvial Plain in Arkansas and Mississippi to estimate hydraulic conductivity and transmissivity values for the Mississippi River Valley alluvial and middle Claiborne aquifers. Well and aquifer data were collected from field measurements, well-construction reports, and published aquifer-thickness information. A total of 324 slug-in and slug-out tests were conducted on 48 wells by using mechanical slugs to displace the water column and submersible pressure transducers to record changes in water levels in the wells. Hydraulic conductivity of the aquifers in which the wells are screened was estimated by curve fitting the water-level-change data using aquifer test analysis software. Estimates of aquifer transmissivity were made by multiplying the estimated hydraulic conductivity value by the aquifer thickness at well locations. Mean hydraulic conductivity estimates for 44 observation wells screened in the Mississippi River Valley alluvial aquifer range from 3 to 401 feet per day, and mean transmissivity estimates range from 285 to 80,559 feet squared per day. Mean hydraulic conductivity estimates for four observation wells screened in units of the middle Claiborne aquifer range from 0.14 to 183 feet per day, and mean transmissivity estimates range from 55 to 67,913 feet squared per day. The results from these tests can be used to improve the understanding of water availability and groundwater migration, to refine groundwater models, and to ultimately provide stakeholders and decisionmakers better information for management of the groundwater resources within the Mississippi Alluvial Plain.

Arkansas, Mississippi↗
USGS2024

New methodology for assessing underground natural gas storage resources – Example from Michigan Basin, United States

Energy consumption in the United States (U.S.) and across the world is shifting away from traditional fossil fuels like coal and oil, and towards natural gas and renewable sources, including hydrogen. Because gas demand is typically greatest during cold seasons and renewable sources sometimes produce variable supplies, it is important to store energy for use when demand exceeds supply. Whereas batteries and tanks typically store energy above ground, geologic (underground) storage may be able to retain much greater quantities of energy over much longer durations (e.g., Matos et al. 2019; Buursink et al. 2023). Consequently, the U.S. Geological Survey (USGS) is developing a methodology to locate new geologic energy storage sites with an initial focus on depleted hydrocarbon reservoirs followed by estimates of associated pore space or gas storage capacity. Depleted hydrocarbon fields in Michigan Basin in the U.S. already host large quantities of injected natural gas (about 681 billion cubic feet of seasonal gas capacity in 42 gas storage facilities; U.S. Energy Information Administration 2024). Therefore, we used this basin as an example to assess potential new storage resources alongside existing storage capacity and carbon dioxide sequestration resources (e.g., Katz and Coats 1968; Haagsma et al. 2020). In Michigan Basin gas storage is taking place in reservoirs within several key hydrocarbon producing formations with different lithologies, including the informal Stray sandstone of the Michigan Formation and the pinnacle reefs of the Niagara Group (Gautier et al. 1995; Swezey et al. 2015). The assessment methodology introduced here consists of two steps relying on both well information and reservoir data.

Michigan↗
USGS2025

Structural controls on splay fault rupture dynamics during Cascadia megathrust earthquakes

Great subduction earthquakes ( M w ≥ 8.0) can generate devastating tsunamis by rapidly displacing the seafloor and overlying water column. These potentially tsunamigenic seafloor offsets result from coseismic fault slip and deformation beneath or within the accretionary wedge. The mechanics of these shallow rupture phenomena and their dependence on subduction zone properties remain unresolved, partly due to the sparsity of offshore observations of shallow megathrust earthquake deformation. Here, we analyze how offshore structure influences shallow rupture mechanics and slip partitioning using 3D dynamic earthquake simulations of the Cascadia subduction zone (CSZ) megathrust with and without variably dipping seaward- or landward-vergent splay faults in the wedge that sole into the megathrust. Resulting tradeoffs between splay and megathrust slip reveal structural controls on rupture partitioning, with greater splay slip leading to less shallow megathrust slip updip. Gently dipping and seaward-vergent splays host more slip than those with steeper, landward-vergent splays. To isolate the underlying mechanisms, we compare models with Andersonian and plunging principal stresses. Results suggest distinct static and dynamic processes control the dip- and vergence-dependence of splay rupture: static (mis)alignment relative to far-field tectonic loading favors slip on more optimally oriented, shallowly dipping splay faults. In contrast, dynamic stress interactions of an updip-propagating megathrust rupture front with the free surface and potential branch faults favor forward branching onto seaward-vergent splays and inhibit backward branching onto landward-vergent splays. Resulting seafloor displacements suggest splay fault structure may influence coseismic tsunami source processes, highlighting the importance of dynamically viable rupture scenarios in subduction hazard assessments.

Cascadia subduction zone↗