Geology Reports⌕ Search

Geology topics

Devin McPhillips

Publications and source records attributed to Devin McPhillips.

At least 19 recordsLinked to original sources
USGS2026

3D semantic mapping of surface geological features

Semantic mapping in 3D is fundamental to a wide range of geoscientific studies and applications, including geomorphology, hazard assessment, and environmental monitoring. However, automatically segmenting geological features from large-scale photogrammetric datasets remains a significant challenge. We present a methodology to address this gap. Using overlapping images collected over environments of interest, Structure-from-Motion (SfM) produces georeferenced point clouds and estimates camera poses. Existing large vision models, such as Segment Anything Model, segment objects in the images, generating pixel-segmentation associations. To produce pixel-point associations, we project the points back onto the camera image planes. As objects are independently segmented across multiple images with different perspectives, we develop a segmentation mosaicking algorithm to build probabilistic point-segmentation associations that combines the pixel-segmentation associations and pixel-point associations. Our methodology is validated using both synthetic data generated by Kubric and real-world UAV-SfM data. The implementation is designed to be compatible with existing SfM software, including Agisoft and OpenDroneMap, for photogrammetry mapping in geoscience studies. As a case study, we apply our method to the semantic mapping of precariously balanced rocks (PBRs), which provide upper-bound constraints on historical ground motion shaking intensity. To support object-level identification of PBRs, we additionally integrated Grounding DINO, enabling text-prompted segmentation of features of interest within UAV imagery. This case study demonstrates the effectiveness of our method in generating a 3D semantic map of PBRs, enabling spatial distribution of PBR fragility for earthquake hazard analysis.

Computers & Geosciences↗
USGS2026

Characterizing precariously balanced rocks (PBRs) in the eastern United States for estimating maximum past earthquake ground motions

Fragile geologic features (FGFs) can be used to estimate an upper bound to the strength of earthquake ground shaking in the past thousands to tens of thousands of years. FGFs could be especially useful in the eastern United States, where few active faults are available for paleoseismic studies and where low seismic attenuation means that an earthquake could damage FGFs at large distances. This article describes a pilot study of one form of FGF, precariously balanced rocks (PBRs), in the Blue Ridge of Virginia, ∼160 km from the central Virginia seismic zone (CVSZ) and ∼100 km from Virginia’s Giles County seismic zone (GCSZ). The main goal is to develop and refine methods for studying PBRs in the eastern United States. Photogrammetric 3D models of the rocks were used to locate their centers of mass and delineate their basal contacts, with the latter corroborated by wrapping a wire around the bases to verify the circumferences. Relations between hill geometry and amplification were used to estimate topographic amplification. Seismometers on the PBRs provided the frequencies, amplitudes, and directions of rocking after giving the rocks gentle pushes. The levels of ground shaking likely to topple the rocks were compared with the median hazard curves from the U.S. Geological Survey’s National Seismic Hazard Model and were used to map magnitudes of earthquakes in the nearby region likely to topple the rocks. Results indicate that an M w ≥ 7.5 earthquake is needed in the CVSZ or near the 2020 Sparta, North Carolina, earthquake to topple the two most fragile PBRs and that an earthquake of ~ M w 7–7.5 in the GCSZ is needed. The fragilities and estimated ages of the PBRs (20–40 ka) are thus consistent with the absence of such events in the historical record in the nearby region.

eastern United States↗
USGS2026

From collision to subduction: Thermal-kinematic inversions constrain plate boundary structure and dip-slip activity in southwestern New Zealand

Constraints on three-dimensional fault system geometry and kinematics are required to understand both long-term tectonics and contemporary seismic hazard. In many cases, these parameters must be inferred from a combination of different data or analyses. Here, we invert low-temperature thermochronological data to quantitatively parameterize both the crustal-scale geometry and late Cenozoic dip-slip activity of the plate boundary system in southwest New Zealand. We find that the modern transpressive plate boundary system has accommodated the bulk of Australian-Pacific Plate convergence since establishment in the middle Miocene-early Pliocene. Shortening has occurred via (a) dextral-reverse slip on the moderate angle, listric central Alpine Fault décollement, and (b) reverse-slip on either the sharply curved, anti-listric northern Puysegur subduction interface or low-angle offshore thrust fault system that splays from the high-angle southern Alpine Fault at depth. We resolve late Cenozoic shortening rates on these structures that are indistinguishable, and interpret these to reflect the transfer of convergent deformation outboard of the predominantly strike-slip southern Alpine Fault. Our work helps elucidate the transition from strain localized oblique continental collision to strain partitioned subduction. We also provide an estimate of the across-strike width of the central-southern Alpine Fault section boundary at seismogenic depths. This is the first quantitative constraint on the geometry of a structural feature inferred to have arrested several major prehistoric earthquakes. This study furthers understanding of fault structure and dip-slip kinematics along the greatest sources of seismic hazard in and offshore of southwest New Zealand.

South Island↗
USGS2025

Seismic response of rock towers at the Trona Pinnacles (U.S.A.) to the 2019 Ridgecrest earthquake sequence: Theory, observations, and models

We analyze the seismic response of a class of fragile geologic features (FGFs), referred to as rock towers (RTs) at the Trona Pinnacles, a group of RTs in southern California that suffered strong shaking during the 2019 Ridgecrest earthquake sequence. FGFs, including RTs, may provide maximum constraints on past earthquake shaking intensity, and thereby support probabilistic seismic hazard assessments (PSHAs). In a rare case study, we explore the hypothesis that RT structural integrity is time dependent, as damage accumulates progressively. We develop finite‐element method (FEM) models of the RTs using photogrammetric shape models. We validate the models by comparing numerical simulations of their response to broadband ground shaking with low‐intensity seismic recordings obtained at the Pinnacles. Results of our simulations are in good agreement with the seismic recordings of actual earthquake aftershocks. We next use the results of the FEM models to analyze the response and evolution of RTs. Our analyses elucidate the influence of geometry over their seismic response, providing a rationale that may explain the rarity of slender RTs at Trona: high‐aspect‐ratio structures that respond in bending develop detrimental tensile stresses that crack the rock, whereas low‐aspect‐ratio ones’ response also includes shearing, which does not compromise material integrity as much as tension. Field measurements with a rebound hammer support this finding, suggesting that the material around the base of slender rocks has been weakened relative to other parts of the RT. We also study how to define simplified mechanical models (“archetypes”) to predict the natural frequencies of RTs. Results from our work illuminate the fundamental mechanisms of seismic response and progressive failure of RTs, and open new avenues of research to potentially incorporate these geologic features as long‐return period constraints on PSHA, in ways analogous to those of the widely used precariously balanced rocks.

California↗
USGS2025

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

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

Bulletin of the Seismological Society of America↗
USGS2025

Strong shaking from past Cascadia Subduction Zone earthquakes encoded in coastal landforms

Strong earthquakes along subduction zones are often devastating events, but sparse records along some tectonic margins limit our understanding of seismic hazards. Constraining shaking intensities is critical, especially in subduction zones with infrequent but large-magnitude earthquakes like the Cascadia Subduction Zone (CSZ), where the lack of recorded ground motions has led to uncertainty in the severity and potential impacts of future earthquakes. Here we fill this observational gap with a novel inventory of quantitative estimates of past shaking intensities from geotechnical modeling of coastal landforms. One hundred fifty-four deep-seated landslides and 65 fragile geologic features constrain minimum and maximum peak ground accelerations, respectively. These estimates are broadly consistent with model predictions of M9 ruptures, suggesting strong shaking of 0.4–0.8 g during past CSZ earthquakes. Local discrepancies between our geologic shaking constraints and earthquake simulations may inform past rupture behavior, leading to better predictions of shaking intensity for future earthquakes.

California, Oregon, Washington↗
USGS2024

Precariously balanced rocks in northern New York and Vermont, U.S.A.: Ground-motion constraints and implications for fault sources

Precariously balanced rocks (PBRs) and other fragile geologic features have the potential to constrain the maximum intensity of earthquake ground shaking over millennia. Such constraints may be particularly useful in the eastern United States (U.S.), where few earthquake‐source faults are reliably identified, and moderate earthquakes can be felt at great distances due to low seismic attenuation. We describe five PBRs in northern New York and Vermont—a region of elevated seismic hazard associated with historical seismicity. These boulders appear to be among the most fragile PBRs in the region, based on reports from hobbyists. The PBRs are glacial erratics, best evidenced by glacial striations on bedrock pedestals. The pedestals themselves are locally high knobs, often situated on regionally high topography; this setting limits soil development and indicates that any outwash deposits were likely ephemeral. As a result, PBR ages can be reliably established by the retreat of the last continental ice sheet, ∼15–13 ka. To quantify the fragility of the PBRs, we surveyed them with ground‐based light detection and ranging and calculated geometric parameters from the point clouds, field observations, and seismic responses. Preliminary validation of the 2023 time‐independent U.S. National Seismic Hazard Model (NSHM) shows that the existence of PBRs is generally consistent with the median site‐specific hazard curves. Only the Blue Ridge Road site suggests a modest reduction in hazard. To visualize the ensemble of data, we mapped the minimum permissible distance to potential source faults around each PBR site as a function of source magnitude by using the ground‐motion models from the 2023 NSHM. Viewed in this manner, our data are consistent with potential M ∼6.5 earthquake‐source faults in many parts of the Lake Champlain Valley and northern Adirondack Mountains. Our work illustrates a potential pathway for better constraining earthquake‐source faults in regions of cryptic faults.

Bulletin of the Seismological Society of America↗
USGS2024

Inbuilt age, residence time, and inherited age from radiocarbon dates of modern fires and late Holocene deposits, Western Transverse Ranges, California

Radiocarbon dates of sedimentary deposits include the elapsed time between formation of the organic material and deposition at the sample site, known as the inherited age. Long inherited ages reduce the accuracy of estimates of the timing of depositional events used to infer paleoclimate change, fire histories, and paleoearthquake timing. An inherited age distribution combines the inbuilt age distribution, which reflects the age composition of the vegetation of the source area, and the residence time distribution, which includes transport and interim storage prior to final deposition. Differentiating residence time and inbuilt age is difficult given typical dispersion of ages in a sedimentary deposit. We address this problem by comparing charcoal dates from two modern fires in southern California, the 2020 Bobcat and the 2013 Grand Fire, with a well-dated late Holocene deposit in the Pallett Creek watershed. The modern fire deposits have negligible transport time (<1 year), and 56 radiocarbon dates indicate a median age of 25 years (300-year 95% range) provides an estimate of inbuilt age for the San Gabriel Mountains. The inherited age calculated from the paleodeposits is older with a median age of ~90 years and has a positive skew (850-year 95% range). A modeled inherited age, calculated by applying the pre-bomb radiocarbon calibration curve to the modern fire age distribution, is shorter than the paleodeposit inherited age by only 21 years, indicating samples with long residence times are not common in the deposit. Comparison of inherited ages calculated from organic-rich and clastic paleodeposits indicate a slight facies dependence that may reflect longer residence time in clastic deposits. The results provide insight into the transport of charcoal through the landscape are useful for refining estimates of past environmental and tectonic events.

California↗
USGS2024

Fault activity in the San Gabriel Mountains, southern California, USA: Insights from landscape morphometrics, erosion rates, and fault-slip rates

Many studies use landscape form to determine spatial patterns of tectonic deformation, and these are particularly effective when paired with independent measures of rock uplift and erosion. Here, we use morphometric analyses and 10 Be catchment-averaged erosion rates, together with reverse slip rates from the Sierra Madre−Cucamonga fault zone, to reveal patterns in uplift, erosion, and fault activity in the range front of the San Gabriel Mountains in southern California, USA. Our analysis tests two prevailing hypotheses: (1) the range front of the San Gabriel Mountains is at steady state, in which rock uplift balances erosion and topographic elevations are stable throughout time, and (2) that west-to-east increases in elevation, relief, erosion rate, and stream-channel steepness across the interior of the massif reflect a parallel reverse-slip rate gradient on the range-bounding Sierra Madre−Cucamonga fault zone. We show that although deviations from steady state occur, the range-front hillslopes and stream channels are typically both well-connected and adjusted to patterns in Quaternary uplift driven by motion on the range-front fault network. Accordingly, landscape morphometrics, 10 Be erosion rates, and model erosion rates effectively image spatial and temporal patterns in uplift. Interpreted jointly, these data reveal comparable peak slip rates on the Sierra Madre−Cucamonga fault zone and show that they do not monotonically increase from west to east. Thus, the eastward-increasing gradients developed within the interior of the massif are not solely related to reverse slip on the range-front faults. Evaluated on shorter length scales (<10 km), morphometric data corroborate earlier descriptions of the Sierra Madre−Cucamonga fault zone as multiple individual faults or fault sections, with slip rates tapering toward fault tips. We infer that these patterns imply the predominance of independent fault or fault section ruptures throughout the Quaternary, though data cannot rule out the possibility of large, connected Sierra Madre−Cucamonga fault zone ruptures. Deeper in the hanging wall of the Sierra Madre−Cucamonga fault zone, secondary faults accommodate range-front uplift. Motion on these faults may contribute to active uplift of the highest topography within the massif, in addition to partly reconciling differences between geologic and geodetic Sierra Madre−Cucamonga fault zone reverse-slip rates. This study provides a new, unified perspective on tectonics and landscape evolution in the San Gabriel Mountains.

California↗
USGS2023

The 2023 US 50-State National Seismic Hazard Model: Overview and implications

The US National Seismic Hazard Model (NSHM) was updated in 2023 for all 50 states using new science on seismicity, fault ruptures, ground motions, and probabilistic techniques to produce a standard of practice for public policy and other engineering applications (defined for return periods greater than ∼475 or less than ∼10,000 years). Changes in 2023 time-independent seismic hazard (both increases and decreases compared to previous NSHMs) are substantial because the new model considers more data and updated earthquake rupture forecasts and ground-motion components. In developing the 2023 model, we tried to apply best available or applicable science based on advice of co-authors, more than 50 reviewers, and hundreds of hazard scientists and end-users, who attended public workshops and provided technical inputs. The hazard assessment incorporates new catalogs, declustering algorithms, gridded seismicity models, magnitude-scaling equations, fault-based structural and deformation models, multi-fault earthquake rupture forecast models, semi-empirical and simulation-based ground-motion models, and site amplification models conditioned on shear-wave velocities of the upper 30 m of soil and deeper sedimentary basin structures. Seismic hazard calculations yield hazard curves at hundreds of thousands of sites, ground-motion maps, uniform-hazard response spectra, and disaggregations developed for pseudo-spectral accelerations at 21 oscillator periods and two peak parameters, Modified Mercalli Intensity, and 8 site classes required by building codes and other public policy applications. Tests show the new model is consistent with past ShakeMap intensity observations. Sensitivity and uncertainty assessments ensure resulting ground motions are compatible with known hazard information and highlight the range and causes of variability in ground motions. We produce several impact products including building seismic design criteria, intensity maps, planning scenarios, and engineering risk assessments showing the potential physical and social impacts. These applications provide a basis for assessing, planning, and mitigating the effects of future earthquakes.

Earthquake Spectra↗
USGS2023

The USGS 2023 Conterminous U.S. time‐independent earthquake rupture forecast

We present the 2023 U.S. Geological Survey time‐independent earthquake rupture forecast for the conterminous United States, which gives authoritative estimates of the magnitude, location, and time‐averaged frequency of potentially damaging earthquakes throughout the region. In addition to updating virtually all model components, a major focus has been to provide a better representation of epistemic uncertainties. For example, we have improved the representation of multifault ruptures, both in terms of allowing more and less fault connectivity than in the previous models, and in sweeping over a broader range of viable models. An unprecedented level of diagnostic information has been provided for assessing the model, and the development was overseen by a 19‐member participatory review panel. Although we believe the new model embodies significant improvements and represents the best available science, we also discuss potential model limitations, including the applicability of logic tree branch weights with respect different types of hazard and risk metrics. Future improvements are also discussed, with deformation model enhancements being particularly worthy of pursuit, as well as better representation of sampling errors in the gridded seismicity components. We also plan to add time‐dependent components, and assess implications with a wider range of hazard and risk metrics.

Bulletin of the Seismological Society of America↗
USGS2023

The Mojave section of the San Andreas fault (California), 1: Shaping the terrace stratigraphy of Littlerock Creek through the competition between rapid strike-slip faulting and lateral stream erosion over the last 40ka.

To determine the post-40 ka slip-rate along the Mojave section of the San Andreas Fault (MSAF) we re-analyze the sedimentary record preserved where Little Rock (LR) Creek flows across the fault. At this location, interaction between the northeast-flowing stream and right-lateral fault has resulted in the abandonment and preservation of 11 strath terraces and one paleo-floodplain in the downstream trailing corner of the river, two of which are also preserved upstream to provide cross-fault matches. A new model of fault-induced river deflection, together with standard terrace riser restoration, yields strike-slip displacements of 1,140 ± 160 m for the older terrace and 360 ± 70 m for the younger one. When combined with new 10 Be dating and reinterpretation of prior measurements the displaced terraces yield right-lateral slip-rates of 27.7 +6.9/−3.5 and 26.8 +3.4/−3.0 mm/yr over the last 23 k.y. and last 40 k.y., where uncertainties are at 95% credible intervals. These new rate determinations are consistent with independent late Holocene estimates, indicating that the long-term rate of strain accumulation along the MSAF is relatively fast and does not vary significantly when averaged over timescales of 15–20 k.y. Using our new model of stream deflection, we find that the fluvial sequence was emplaced in two distinct periods, each characterized by a temporally stable but markedly different deflected river geometry. Each period coincides with a distinct stage of erosive power along LR Creek determined from independent paleoclimate proxies. Importantly, application of the new river-deflection model allows strike-slip displacements to be determined in the absence of upstream piercing points.

California↗
USGS2023

Spatial and temporal analysis of geologic slip rates, Cucamonga Fault, California, USA: Implications for along-strike applications and multi-fault rupture

To constrain fault processes and hazard, fault slip rates may be extrapolated over different fault lengths or time intervals. Here, we investigate slip rates for the Cucamonga Fault (CF). The CF is located at the junction of the Transverse Range fault system with the San Andreas and San Jacinto Faults, and it is hypothesized to connect with these faults, promoting the propagation of large, multi-fault earthquakes. Previous work has shown that CF displacements on late Quaternary alluvial fan surfaces are highly variable along strike. We present two new 10 Be surface exposure ages from depth profiles on the alluvial fans. Slip rates are consistent with a rate of 1.4 ± 0.3 m/kyr over time intervals of ∼20, ∼30, and ∼40 kyr. If the CF participates in multi-fault ruptures, then these earthquakes occur either rarely or with sufficient regularity to maintain apparently steady rates over multiple intervals. We also explore along-strike fault displacement variability using a calibrated morphological model. The model successfully reproduces scarp profiles and indicates that fault displacement variability can be explained in part by scarp age but not uplift rate. We infer that both erosion by ephemeral gullying and distributed deformation contribute to fault displacement variability, although both are difficult to detect confidently without excavations across the scarp. These investigations show that better characterization of cumulative-slip variability along strike may improve accuracy and precision of slip rates. Slip rates that do not consider epistemic uncertainties may not be suitable for extrapolation over longer fault sections.

California↗
USGS2022

Revised earthquake recurrence intervals in California, USA: New paleoseismic sites and application of event likelihoods

Recurrence intervals for ground rupturing earthquakes are critical data for assessing seismic hazard. Recurrence intervals are presented here for 38 paleoseismic sites in California. Eleven of these include new or updated data; the remainder use data previously included in the Unified California Earthquake Rupture Forecast Version 3 (UCERF3). The methods and results are consistent with UCERF3. In addition, revised recurrence intervals are presented at every site. The revised recurrence intervals incorporate uncertainty in the interpretation of paleoseismic evidence, which is expressed as event likelihood. Event likelihood is the probability that the evidence has been correctly interpreted as a unique earthquake. Event likelihoods are estimated here for 85 inferred past earthquakes at eight paleoseismic sites in California, using a single, consistent methodology. The average event likelihood is 0.85. The revised recurrence intervals are 16% longer, on average, than conventional estimates, and their confidence intervals are disproportionately wider. These recurrence intervals are suitable for inclusion in a “grand inversion” rupture forecast, and they may be important for addressing a systematic misfit in the UCERF3 grand inversion. The revised recurrence intervals may also be important for assessing the unusually long earthquake hiatus in California. Other applications may not need to consider event likelihoods because the effects are small relative to typical uncertainties.

California↗
USGS2022

Toppling of a Trona Pinnacles Spire following the M5.5 RidgecrestaAftershock of June 2020

The 2019 M w 7.1 Ridgecrest California earthquake rupture passed within 4 km of the Trona Pinnacles, a large group of tufa rock pillars. Reconnaissance following the Ridgecrest mainshock documented fresh damage to several of the Pinnacles. Repeated aerial photogrammetric surveys also documented damage during subsequent aftershocks. Here, we describe the photogrammetric data with emphasis on a specific rock spire that toppled during an M w 5.5 aftershock. We calculate the volumes of the intact spire and of its subsequent debris. To explore the utility of the pinnacles as fragile geologic features for constraining past earthquake shaking intensity, we calculate the quasi static, horizontal acceleration required to break the spire at its base. We also examine the response of this feature to observed shaking using a dynamic model of the spire. In this case, we find that the quasi‐static estimate provides a conservative maximum constraint on fragility. The dynamic model of the spire suggests that shaking during the M w 7.1 mainshock likely generated tensile stresses in excess of the spire’s bulk strength, thereby making it vulnerable to collapse in subsequent aftershocks.

California↗
USGS2021

Survey of fragile geologic features and their quasi-static earthquake ground-motion constraints, southern Oregon

Fragile geologic features (FGFs), which are extant on the landscape but vulnerable to earthquake ground shaking, may provide geological constraints on the intensity of prior shaking. These empirical constraints are particularly important in regions such as the Pacific Northwest that have not experienced a megathrust earthquake in written history. Here, we describe our field survey of FGFs in southern Oregon. We documented 58 features with fragile geometric characteristics, as determined from field measurements of size and strength, historical photographs, and light detection and ranging point clouds. Among the surveyed FGFs, sea stacks have particular advantages for use as ground‐motion constraints: (1) they are frequently tall and thin; (2) they are widely distributed parallel to the coast, proximal to the trench and the likely megathrust rupture surface; and (3) they are formed by sea cliff retreat, meaning that their ages may be coarsely estimated as a function of distance from the coast. About 40% of the surveyed sea stacks appear to have survived multiple Cascadia megathrust earthquakes. Using a quasi‐static analysis, we estimate the minimum horizontal ground accelerations that could fracture the rock pillars. We provide context for the quasi‐static results by comparing them with predictions from kinematic simulations and ground‐motion prediction equations. Among the sea stacks old enough to have survived multiple megathrust earthquakes ( n = 16), eight yield breaking accelerations lower than the predictions, although they generally overlap within uncertainty. FGFs with the lowest breaking accelerations are distributed uniformly over 130 km of coastline. Results for inland features, such as speleothems, are in close agreement with the predictions. We conclude that FGFs show promise for investigating both past earthquake shaking and its spatial variability along the coasts of Oregon and Washington, where sea stacks are often prevalent. Future work can refine our understanding of FGF age and evolution.

Oregon↗
USGS2021

Holocene depositional history inferred from single-grain luminescence ages in southern California, North America

Significant sediment flux and deposition in a sedimentary system are influenced by climate changes, tectonics, lithology, and the sedimentary system's internal dynamics. Identifying the timing of depositional periods from stratigraphic records is a first step to critically evaluate the controls of sediment flux and deposition. Here, we show that ages of single-grain K-feldspar luminescence subpopulations may provide information on the timing of previous major depositional periods. We analyzed 754 K-feldspar single-grains from 17 samples from the surface to ∼9 m-depth in a trench located downstream of the Mission Creek catchment. Single-grain luminescence subpopulation ages significantly overlap at least eight times since ∼12.0 ka indicating a common depositional history. These depositional periods correspond reasonably well with the Holocene intervals of wetter than average climate conditions based on hydroclimatic proxies from nearby locations. Our findings imply a first-order climatic control on sediment depositional history in southern California on a millennial timescale.

California↗