USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Research of the U.S. Geological Survey”

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,697 records · Page 88Linked to original sources
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↗
USGS2025

Pathways for potential exposure to onshore oil and gas wastewater: What we need to know to protect human health

Produced water is a chemically complex waste stream generated during oil and gas development. Roughly four trillion liters were generated onshore in the United States in 2021 (ALL Consulting, 2022, https://www.gwpc.org/wp-content/uploads/2021/09/2021_Produced_Water_Volumes.pdf ). Efforts are underway to expand historic uses of produced water to offset freshwater needs in water-stressed regions, avoid induced seismic activity associated with its disposal, and extract commodities. Understanding the potential exposures from current and proposed produced water uses and management practices can help to inform health-protective practices. This review summarizes what is known about potential human exposure to produced water from onshore oil and gas development in the United States. We synthesize 236 publications to create a conceptual model of potential human exposure that illustrates the current state of scientific inquiry and knowledge. Exposure to produced water can occur following its release to the environment through spills or leaks during its handling and management. Exposure can also arise from authorized releases, including permitted discharges to surface water, crop irrigation, and road treatment. Knowledge gaps include understanding the variable composition and toxicity of produced water released to the environment, the performance of treatment methods, migration pathways through the environment that can result in human exposure, and the significance of the exposures for human and ecosystem health. Reducing these uncertainties may help in realizing the benefits of produced water use while simultaneously protecting human health.

Environmental Health Perspectives↗
USGS2026

Deformed submarine terraces in Puget Sound, Pacific Northwest, indicate only one M >~7.5 earthquake on the Seattle fault zone in the past 11,000 yr

Submerged marine terraces in Puget Sound, deformed across the Seattle fault zone (SFZ), indicate that only one earthquake as large as M~7.5 has occurred in at least the past 11 kyr. Previous paleoseismic studies document a M~7.5 earthquake between 923–4 CE, which uplifted coastal marine terraces by as much as 8 m. We demonstrate that this earthquake was the only such event since ~11 ka by mapping and quantifying deformation of older marine terraces, now submerged in Puget Sound. The submerged terraces, attributed to a late-glacial sea-level lowstand, record both glacial isostatic rebound and tectonic deformation. Vertical offset of the ~11 ka terraces within the SFZ is comparable to that of the marine terraces uplifted in 923 CE, implying no additional large (M>~7.5) earthquake on the SFZ since ~11 ka. This result implies a longer recurrence interval than current hazard estimates, which assumes recurrence of M>7.1 events every 5 kyr. Our mapping of SFZ deformation since ~11 ka also supports fault segmentation and contiguous block uplift between the Seattle and Tacoma fault zones.

Washington↗
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↗
USGS2026

Riverscape heterogeneity shapes population diversity for a migratory fish

Habitat patch dynamics can scale up to influence population demography and diversity with implications for resilience to environmental stochasticity. But how the spatial arrangement and size of habitat patches interact with other components of habitat heterogeneity to shape population diversity at larger spatial scales is not well understood. For riverine fishes, there is increasing evidence that tributary streams provide critical demographic support to main stem rivers. However, the extent to which main stem rivers rely on demographic contributions from tributaries, and the factors underlying this dependence, have not been assessed. Here, we used genetic stock identification to evaluate the effect of tributaries on population diversity of Yellowstone cutthroat trout ( Oncorhynchus virginalis bouvieri ) occupying the main stem Snake River, Wyoming, USA. We found that the main stem relied almost entirely on tributaries for demographic support, but main stem composition varied spatially among river sections. Distance between habitat patches, catchment area, and groundwater availability acted in concert to determine the contribution of specific tributaries to the main stem, but contributions were ultimately modulated by habitat connectivity. We also found evidence for multi-scale spatial structure in tributary contributions, providing insight into untested drivers of main stem river population diversity. Our results demonstrate how spatially discrete and distributed riverscape attributes influence population diversity at broader spatial scales, illustrating how ecosystem resilience emerges from the dynamic, two-way exchange of individuals and energy across habitat networks. Management plans for large rivers that address the ecological contributions of tributaries may be needed to achieve optimal outcomes. Similarly, conservation strategies that exclusively focus on headwater streams may fail to capture the broader habitat requirements necessary to maintain robust cold-water fish populations and associated recreational fisheries, particularly under global environmental change.

Wyoming↗
USGS2025

Global patterns of coseismic landslide runout mobility differ from aseismic landslide trends

Coseismic landslides significantly contribute to human and economic losses during and immediately following earthquakes, yet very little data on the runout of such landslides exist. While well-established behavior of aseismic (e.g., hydrologically triggered) landslide runout mobility suggests strong correlation between landslide size and mobility, limited studies of coseismic landslide runout find conflicting mobility trends. We present a global dataset of runout lengths produced from a new automated method for estimating landslide runout, developed and validated using 1726 manually mapped landslides from five unique earthquakes. We then apply the automated runout tool to 23 global earthquake-induced landslide inventories, producing a compiled database of 73,665 measured and estimated runout lengths of coseismic landslides to assess mobility trends. We find a significant divergence between well-established aseismic mobility trends and that of coseismic landslides, with far greater scatter and more complex mobility patterns in earthquake-triggered landslides. As a function of landslide size, we observe global coseismic landslide mobility patterns are bilinear, becoming increasingly less mobile with increasing size above some threshold. This discordance between aseismic and coseismic landslide mobility may be a function of landslide type, kinematics, hydrology, and or setting that systematically differ between triggering mechanisms and should be explored in more depth to develop predictive models of these unique runout patterns. These results suggest hazard and risk models for coseismic landslides may significantly under-predict or over-predict impacts, depending on the size of triggered landslides.

Engineering Geology↗
USGS1976

Distribution functions for statistics derived from bivariate normal and bivariate two-parameter log-normal populations

The distribution functions for statistics that may be used to assess the significance of differences between sample means, standard deviations, coefficients of skewness, and coefficients of variation are obtained by Monte Carlo experiments. The distributions are expressed as functions of sample size, cross correlation, and skewness. In general, the distributions are more sensitive to cross correlation than to skewness. As sample size increases, however, the distributions tend to become more sensitive to skewness.

Open-File Report↗
USGS2025

Factors influencing landslide occurrence in low-relief formerly glaciated landscapes: Landslide inventory and susceptibility analysis in Minnesota, USA

In landscapes recently impacted by continental glaciation, landslides may occur where topographic relief has been generated by the drainage of glacial lakes and ensuing post-glacial fluvial network development into unconsolidated glacially derived sediments and exhumed bedrock. To investigate linkages among environmental variables, post-glacial landscape development, and landslides, we created a landslide inventory of nearly 10,000 landslides in five regions of the formerly glaciated low-relief state of Minnesota, USA. Multivariate logistic regression indicates the importance of slope angle, lithology, and the development of stream valleys to landslide distribution. Areas underlain by fine-grained glaciolacustrine and nearshore deposits that are incised by streams are particularly prone to shallow (<1-2 m depth) landslides. Landslides also occur in a wide range of glacial and fluvial deposits, and as rockfall in layered Paleozoic sedimentary rocks in central and southern Minnesota and Precambrian igneous and sedimentary rocks in northeastern Minnesota. Although no more than 1-2% of the studied regions are susceptible to landslides, they can pose risk to life and safety, damage infrastructure, and impact water quality. The combination of recently generated low-relief steep slopes, extensive unconsolidated sediments, and layered sedimentary bedrock make this formerly glaciated landscape more susceptible to landslides than current national-scale models indicate.

Minnesota↗
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↗
USGS2025

Factors regulating the potential for freshwater mineral soil wetlands to function as natural climate solutions

There are increasing global efforts and initiatives aiming to tackle climate change and mitigate its impacts via natural climate solutions (NCS). Wetlands have been considered effective NCS given their capacity to sequester and retain atmospheric carbon dioxide (CO 2 ) while also providing a myriad of other ecosystem functions that can assist in mitigating the impacts of climate change. However, wetlands have a dual impact on climate, influencing the atmospheric concentrations of both CO 2 and methane (CH 4 ). The cooling effect associated with wetland CO 2 sequestration can be counterbalanced by the warming effect caused by CH 4 emissions from wetlands. The relative ability of wetlands to sequester CO 2 versus emit CH 4 is dependent on a suite of interacting physical, chemical, and biological factors, making it difficult to determine if/which wetlands are considered important NCS. The fact that wetlands are embedded in landscapes with surface and subsurface hydrological connections to other wetlands (i.e., wetlandscapes) that flow over and through geochemically active soils and sediments adds a new layer of complexity and poses further challenges to understanding wetland carbon sequestration and greenhouse gas fluxes at large spatial scales. Our review demonstrates how additional scientific advances are required to understand the driving mechanisms associated with wetland carbon cycling under different environmental conditions. It is vital to understand wetland functionality at both wetland and wetlandscape scales to effectively implement wetlands as NCS to maximize ecological, social, and economic benefits.

Wetlands↗
USGS2026

HyFlood: A surrogate-model-based framework for compound coastal flooding

Compound coastal flooding is a major threat to low-lying coastal regions and is expected to intensify under future climate change projections. However, modeling the joint interaction of waves, storm surge, tides, and rainfall remains computationally demanding, limiting the development of fast and reliable forecast tools. Here we present HyFlood, a hybrid statistical-numerical downscaling framework capable of computing and mapping high-resolution compound flood hazards while substantially reducing the computational cost compared with fully process-based hydrodynamic modeling. HyFlood combines statistical sampling and selection algorithms with a cascade of reduced-complexity surrogate models that emulate nearshore wave transformation, surf-zone hydrodynamics, and coastal, fluvial, and pluvial flooding. The surrogate models employ machine-learning and regression algorithms applied to a low-dimensional representation of the flooding outputs, obtained through statistical dimensionality reduction. The framework is demonstrated in southern O'ahu, Hawai'i, a region exposed to elevated sea levels driven by tides, waves, and storm surge along with frequent precipitation-driven flash flooding. Validation of the surrogates against the physics-based model outputs demonstrates that HyFlood accurately reproduces daily maxima of spatially distributed flooding depths. This hybrid approach offers a scalable and efficient tool to better quantify how changes in flooding drivers translate into hazard and impact assessments, and to support compound-flood risk assessments and climate-change adaptation planning.

Hawaii↗
USGS2026

An entropic explanation for Gutenberg-Richter scaling

We develop a simple explanation for Gutenberg-Richter (G-R) size scaling of earthquakes on a single fault. We discretize the fault and consider all possible contiguous ruptures at that level of discretization. In this static model, we assume that slip scales with rupture length, and that the rupture rates at each point along the fault are consistent with an a priori long-term slip rate. These simple assumptions define an (under-determined) non-negative least-squares inverse problem. Each solution to this inverse problem is a set of earthquake rates that matches the slip-rate constraint. We use a Markov Chain Monte Carlo (MCMC) algorithm to uniformly sample the solution space assuming constant slip rates along the fault. At finer discretizations, deviations from G-R behavior decrease, which is consistent with an entropic pressure towards G-R solutions. When the fault is discretized into 10 or more segments, random solutions found by the MCMC algorithm have G-R size scaling, even though there are trivial solutions that, for example, have earthquakes of only one size. This is because there are simply far more solutions that have G-R scaling; as the problem size increases, the strong degeneracy of GR solutions results in other solutions becoming improbably rare. Also, the entropically favored G-R distribution has a b -value of approximately 1, which agrees with measured b -values in real earthquake catalogs.

JGR Solid Earth↗
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↗
USGS2025

Groundwater dominates snowmelt runoff and controls streamflow efficiency in the western United States

Climate change in seasonally snow-covered mountain catchments is reducing water supply and decreasing streamflow predictability. Here, we use tritium age dating to show that contrary to the common assumption that snowmelt quickly contributes to runoff, streamflow during snowmelt in western US catchments is dominated by older groundwater. The average age of streamwater during snowmelt runoff (5.7 ± 4.3 years) was intermediate to the average age of groundwater (10.4 ± 4.5 years) and recent precipitation, indicating that 58% (±34%) of snowmelt runoff was derived from groundwater. Water ages, streamflow, and groundwater storage were mediated by bedrock geology: low-permeability hard rock/shale catchments exhibited younger ages, less storage, and more efficient streamflow generation than high-permeability sandstone/clastic catchments. Our results demonstrate that snowmelt runoff is the result of multiple prior years of climate mediated by groundwater storage. Including these interactions will be crucial for predicting water resources as climate and landscape changes accelerate.

western united States↗
USGS2024

Aboveground carbon stocks across a hydrological gradient: Ghost forests to non-tidal freshwater forested wetlands

Upper estuarine forested wetlands (UEFWs) play an important role in the sequestration of atmospheric carbon (C), which is facilitated by their position at the boundary of terrestrial and maritime environments but threatened by sea level rise. This study assessed the change in aboveground C stocks along the estuarine–riverine hydrogeomorphic gradient spanning salt-impacted freshwater tidal forested wetlands to freshwater forested wetlands in seasonally tidal and nontidal landscape positions. Standing stocks of C in forested wetlands were measured along two major coastal river systems, the Winyah Bay in South Carolina and the Savannah River in Georgia (USA), replicating and expanding a previous study to allow the assessment of change over time. Aboveground C stocks on these systems averaged 172.9 Mg C ha −1 , comparable to those found in UEFWs across the globe and distinct from the terrestrial forested ecosystems they are often considered to be a part of during large-scale C inventory efforts. Groundwater salinity conditions as low as 1.3 ppt were observed in conjunction with losses of aboveground C. When viewed in context alongside expected sea level rise and corresponding saltwater intrusion estimates, these data suggest a marked decrease in aboveground C stocks in forested wetlands situated in and around tidal estuaries.

Georgia, South Carolina↗
USGS2024

Debris-flow entrainment modelling under climate change: Considering antecedent moisture conditions along the flow path

Debris-flow volumes can increase along their flow path by entraining sediment stored in the channel bed and banks, thus also increasing hazard potential. Theoretical considerations, laboratory experiments and field investigations all indicate that the saturation conditions of the sediment along the flow path can greatly influence the amount of sediment entrained. However, this process is usually not considered for practical applications. This study aims to close this gap by combining runout and hydrological models into a predictive framework that is calibrated and tested using unique observations of sediment erosion and debris-flow properties available at a Swiss debris-flow observation station (Illgraben). To this end, hourly water input to the erodible channel is predicted using a simple, process-based hydrological model, and the resulting water saturation level in the upper sediment layer of the channel is modelled based on a Hortonian infiltration concept. Debris-flow entrainment is then predicted using the RAMMS debris-flow runout model. We find a strong correlation between the modelled saturation level of the sediment on the flow path and the channel-bed erodibility for single-surge debris-flow events with distinct fronts, indicating that the modelled water content is a good predictor for erosion simulated in RAMMS. Debris-flow properties with more complex flow behaviour (e.g., multiple surges or roll waves) are not as well predicted using this procedure, indicating that more physically complete models are necessary. Finally, we demonstrate how this modelling framework can be used for climate change impact assessment and show that earlier snowmelt may shift the peak of the debris-flow season to earlier in the year. Our novel modelling framework provides a plausible approach to reproduce saturation-dependent entrainment and thus better constrain event volumes for current and future hazard assessment.

Illgraben basin↗
USGS2026

Relating glassy rind thicknesses to ambient air temperatures at the Lost Jim flow field in the Imuruk Lake volcanic field, Alaska

The Lost Jim flow field, in the Imuruk Lake volcanic field, Alaska, extends west ~34 km from a single vent, crossing subarctic tundra and currently touches several lakes and streams. The weighted mean of five 36Cl cosmogenic exposure ages from the Lost Jim pāhoehoe flow is 7.73 ± 0.37 ka, indicating this eruption occurred substantially after the eruption of the underlying Camille flow, which was emplaced at 39.7 ± 1.3 ka. Paleoclimate records indicate the period when the Lost Jim flow field was emplaced was after deglaciation, and the climate was similar to today. We propose that the emplacement of lava in these cold subarctic conditions can lead to faster cooling of the lava surface compared to lava emplaced in warmer locations such as mid- latitude cold deserts. Glass abundance in the outermost rinds at the Lost Jim flow field was on average 74 % with 6.4 mm thick rims, compared to 60 % with 2.9 mm rims for cold mid-latitude desert samples. We interpret increased glass content as a proxy for rapid cooling likely occurring partly during winter. Glassiness values varied less across vent, margin, and mid-flow locations when compared to the mid-latitude flows suggesting the Lost Jim flow field was broadly impacted by the subarctic climate as opposed to responding to local microclimates. Our results indicate that lava glassiness may be a useful environmental indicator of cooler (in this case subarctic) conditions.

Alaska↗
USGS2026

Early Miocene volcanic rocks and associated tectonics, Lava Hills and southern Bristol Mountains, California

Volcanic rocks of latest Oligocene to early Miocene age form an east-west belt across part of the central eastern Mojave Desert from the Whipple Mountains on the east to the Rosamond Hills on the west. We term this the central belt because it is separated from northern and southern belts by swaths with no volcanic rocks. Limited geochronologic data indicate that much of the belt is latest Oligocene and early Miocene in age, about 24 to 19 Ma, a finding that is consistent with these rocks being overlain by the 18.8 Ma Peach Spring Tuff in many places. We describe Miocene geology in a central area of the belt, in the Lava Hills, southern Bristol Mountains, and southern Old Dad Mountains. Sedimentary basins formed coeval with early andesite to rhyolite volcanism, progressing from fluvial and lacustrine tuffaceous sandstone to volcanic lavas, tuffs, and breccias, indicating that early basins formed proximal to volcanic edifices. Higher strata are fluvial and lacustrine with lavas punctuating the sequence. Although basins may partly have been formed within topographic lows bounded by volcanic domes, plateaus, and stratovolcanoes, consistent stratigraphic sections over wide areas indicate that tectonic basin evolution affected broad areas. The volcanic section is capped by local basalt flows and the regional Peach Spring Tuff. Limited data on normal faults support interpretations of early extensional basin development caused by northeast-southwest oriented stretching. Later extension caused stratal rotations, tilting early deposits down to the southwest. This tilted and subsequently beveled basin architecture was overlain by the youngest volcanic deposits, primarily rhyolite and basalt. The Peach Spring Tuff, 18.8 Ma, lies within this upper unit. Similar stratigraphic and structural relations are exposed in the nearby Marble Mountains and Van Winkle Mountain sections, reinforcing that a broad area underwent similar volcanism and tectonism. In our study area the upper unit is only very gently tilted except near dextral strike-slip faults of the eastern California shear zone. These late Miocene to Recent faults are represented as four main fault zones spaced about 5 km apart, representing distributed shear north of the Bristol Lake basin.

California↗