Geology topics

Samuel Johnstone

Publications and source records attributed to Samuel Johnstone.

At least 19 recordsLinked to original sources

USGS2026

Impacts of transient critical zone evolution on landslide susceptibility and hazard Across the Utuado batholith, Puerto Rico

Landscapes in humid climates often experience rapid chemical and mechanical weathering that can change groundwater fluxes, drainage pathways, and frictional strength via alterations in subsurface mineral composition and matrix porosity. These changes to hillslope critical zone structure and strength should control landslide susceptibility, size, and failure style. Here we test the impact of differential weathering on landslide failure processes by exploiting a natural experiment along the dioritic Utuado batholith in Puerto Rico's Cordillera Central, where a sharp weathering gradient follows an upstream-propagating knickzone into an uplifted shoreline platform that was affected by Hurricane Maria in 2017. Through a combination of landslide mapping, field geophysical and geotechnical surveys, topographic analysis, calibrated stream profile inversion, 3D back-analysis of hillslope strength, and debris-flow inundation assessments, we show that differential weathering induces changes to critical zone structure and consequently landslide susceptibility and hazard. Our mapping reveals that the hurricane-induced landslides preferentially cluster along less weathered hillslopes despite a marked increase in frictional strength. Jointly inverted seismic refraction and electrical resistivity profiles show that the frictionally weak, heavily weathered hillslopes are highly porous and dry, which we posit acts to stabilize slopes through limiting deeper hydrologic connectivity that can produce shallow water tables with destabilizing lateral seepage forces during storms. Debris-flow inundation mapping and analysis reveals that major contributors to debris-flow growth, channel steepness and landslide coalescence, scale linearly with hillslope friction angle. This leads to higher-mobility debris flows along knickzone reaches, which disproportionately affects communities at the base of the plateau.

Puerto Rico↗

USGS2026

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

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

Nebraska↗

USGS2025

A synthesis engine for constructing geologic maps of the United States

The geologic history of the United States is cataloged in thousands of geologic maps produced during many decades. However, the disparate nature of these individual maps makes it challenging to assess resources, research geologic histories, or characterize natural hazards holistically across the Nation. The U.S. House of Representatives 2020 appropriations bill for the U.S. Department of the Interior (H.R. 116-100) requires the U.S. Geological Survey to “bring together detailed national and continental-resolution [two-dimensional] and [three-dimensional] information produced throughout the Survey and by [F]ederal and [S]tate partners.” In response to this directive, this report presents a compilation and synthesis of geologic maps across the United States in the form of a relational database. The synthesis database includes thematic maps that synthesize the Nation’s geology, and retains the original input maps as well as linkages to standardized vocabularies to aid the discoverability of geologic information. Specifically, the synthesis database is targeted toward producing four National-resolution maps for the conterminous United States: Quaternary geology, the geology at the Earth’s surface, pre-Quaternary geology, and Precambrian geology. In addition, the synthesis database includes the infrastructure necessary to expand to additional resolutions in the future.

conterminous United States↗

USGS2022

The National Cooperative Geologic Mapping Program’s U.S. GeoFramework Initiative — Delivering a digital database of geologic map information by 2030

The U.S. Geological Survey's National Cooperative Geologic Mapping Program's U.S. GeoFramework Initiative responds to Congressional direction to implement the National Geologic Map Database Phase Three. The U.S. GeoFramework Initiative aims to create and deliver a new digital geologic framework of the United States by 2030. This Fact Sheet outlines the collaborations needed and plans taken to achieve that goal.

Fact Sheet↗

USGS2022

The Seamless Integrated Geologic Mapping (SIGMa) extension to the Geologic Map Schema (GeMS)

Geologic maps are the fundamental building blocks of surface and subsurface three-dimensional geologic framework models of the Earth’s crust. However, as the production and availability of geologic map databases continues to increase, inconsistent data models and the lack of synthesized, national geologic map data at scales appropriate for informed decision making negatively affect the functional integration of geologic map data with other national datasets. The Geologic Map Schema (GeMS) is the publication and archive database standard for geologic map data funded by the U.S. Geological Survey National Cooperative Geologic Mapping Program, and standardizes the organization and content of a single map database. However, synthesizing multiple databases into a seamless geologic map database creates a different set of challenges and database needs than GeMS was designed to accommodate. The Seamless Integrated Geologic Mapping (SIGMa) extension is designed to expand the capabilities of GeMS by enabling integration of map-based geoscience data. In particular, the SIGMa extension enables capturing a diverse and ever-changing set of map units, produced by many contributors operating independently, and by incremental and noncontiguous assembly and publication. Feature-level metadata fields allow data sources and digital compilation methods to be attributed separately and a relational structure is designed to support the link between data sources and features attributed with multiple data sources. Instead of paragraph-style map-unit descriptions that can be highly inconsistent, SIGMa parses fundamental map-unit attributes, including material, genetic process, and age, into thematically specific fields. The SIGMa extension uses a hierarchical map-unit organization to facilitate a dynamic and evolving, formation-level stratigraphic framework. The hierarchy is developed around geologic provinces that represent temporally restricted geologic events, processes, and settings. Geologic provinces can include magmatic events, depositional settings associated with tectonic processes or stable continental margins, and processes that are actively shaping the modern landscape. A geologic province hierarchy places map units into a geologic context at subregional to continental scales and provides the flexibility to support incremental assembly of the stratigraphy.

Scientific Investigations Report↗

USGS2022

The Pondosa fault zone: A distributed dextral-normal-oblique fault system in northeastern California, USA

The tectonic domains of Basin and Range extension, Cascadia subduction zone contraction, and Walker Lane dextral transtension converge in the Mushroom Rock region of northeastern California, USA. We combined analysis of high-resolution topographic data, bedrock mapping, 40 Ar/ 39 Ar geochronology, low-temperature thermochronology, and existing geologic and fault mapping to characterize an extensive dextral-normal-oblique fault system called the Pondosa fault zone. This fault zone extends north-northwest from the Pit River east of Soldier Mountain, California, into moderately high-relief volcanic topography as far north as the Bartle (California) townsite with normal and dextral offset apparent in geomorphology and fault exposures. New and existing 40 Ar/ 39 Ar and radiocarbon dating of offset lava flows provides ages of 12.4 ka to 9.6 Ma for late Cenozoic stratigraphic units. Scarp morphology and geomorphic expression indicate that the fault system was active in the late Pleistocene. The Pondosa fault zone may represent a dextral-oblique accommodation zone between north-south–oriented Basin and Range extensional fault systems and/or part of the Sierra Nevada–Oregon Coast block microplate boundary.

California↗

USGS2022

Mountain rivers reveal the earthquake hazard of geologic faults in Silicon Valley

The 1989, M w = 6.9 Loma Prieta earthquake resulted in tens of lives lost and cost California almost 3% of its gross domestic product. Despite widespread damage, the earthquake did not clearly rupture the surface, challenging the identification and characterization of these hidden hazards. Here, we show that they can be illuminated by inverting fluvial topography for slip-and moment accrual-rates—fundamental components in earthquake hazard assessments—along relief-generating geologic faults. We applied this technique to thrust faults bounding the mountains along the western side of Silicon Valley in the San Francisco Bay Area, and discovered that these structures may be capable of generating a M w = 6.9 earthquake every 250–300 years based on moment accrual rates. This method may be deployed broadly to evaluate seismic hazard in developing regions with limited geological and geophysical information.

California↗

USGS2022

Continental shelves as detrital mixers: U-Pb and Lu-Hf detrital zircon provenance of the Pleistocene–Holocene Bering Sea and its margins

Continental shelves serve as critical transfer zones in sediment-routing systems, linking the terrestrial erosional and deep-water depositional domains. The degree to which clastic sediment is mixed and homogenized during transfer across broad shelves has important implications for understanding deep-sea detrital records. Wide continental shelves are thought to act as capacitors characterized by transient sediment storage during sea level rise and sediment remobilization during sea level fall. This study attempts to test the hypothesis that sea level lowstand yields more efficient and direct sediment transfer from fluvial sources to deep-sea sinks compared to highstand when sediment is sequestered and mixed on the shelf. We test this by evaluating U-Pb and Lu-Hf detrital zircon provenance trends along the vast Bering Sea shelf and deep-marine Beringian continental margin. We present 5884 U-Pb ages and 402 Lu-Hf analyses from 30 samples to characterize the provenance of modern to Pleistocene sediment across the Bering Sea region. We used both forward and inverse numerical mixture modeling to estimate the abundance of distinct fluvial sources in shelfal and deep-water deposits. These results demonstrate that sediment in the Bering Sea is derived from a mixture of regional fluvial sources, but that the Yukon River is the primary detrital source for sediment throughout the region. Although Yukon River signatures are abundant in all basin samples, the relative proportions of Yukon vs other sources vary spatially across the shelf. A comparison of Holocene and surficial sediment with Pleistocene deposits shows that sediment across the shelf and in the deep-sea remains well-mixed between climate states. Thus, detrital provenance signatures in deep-marine deposits outward of broad transfer zones are likely to represent mixtures of fluvial sources regardless of sea level.

Bering Sea↗

USGS2021

Detrital signals of coastal erosion and fluvial sediment supply during glacio-eustatic sea-level rise, Southern California, USA

Coastal erosion, including sea-cliff retreat, represents both an important component of some sediment budgets and a significant threat to coastal communities in the face of rising sea level. Despite the importance of predicting future rates of coastal erosion, few prehistoric constraints exist on the relative importance of sediment supplied by coastal erosion versus rivers with respect to past sea-level change. We used detrital zircon U-Pb geochronology as a provenance tracer of river and deep-sea fan deposits from the Southern California Borderland (United States) to estimate relative sediment contributions from rivers and coastal erosion from late Pleistocene to present. Mixture modeling of submarine canyon and fan samples indicates that detrital zircon was dominantly (55%–86%) supplied from coastal erosion during latest Pleistocene (ca. 13 ka) sea-level rise, with lesser contributions from rivers, on the basis of unique U-Pb age modes relative to local Peninsular Ranges bedrock sources. However, sediment that was deposited when sea level was stable at its highest and lowest points since the Last Glacial Maximum was dominantly supplied by rivers, suggesting decreased coastal erosion during periods of sea-level stability. We find that relative sediment supply from coastal erosion is strongly dependent on climate state, corroborating predictions of enhanced coastal erosion during future sea-level rise.

California↗

USGS2021

A Bayesian nonparametric approach to unmixing detrital geochronologic data

Sedimentary deposits constitute the primary record of changing environmental conditions that have acted on Earth’s surface over geologic time. Clastic material is eroded from source locations (parents) in sediment routing systems and deposited at sink locations (children). Both parents and children have characteristics that vary across many different dimensions, including grain size, chemical composition, and the geochronologic age of constituent detrital minerals. During transport, sediment from different parents is mixed together to form a child, which in turn may serve as the parent for other sediment farther down-system or later in time when buried sediment is exhumed. The distribution of detrital mineral ages observed in parent and child sediments allows for investigation of the proportion of each parent in the child sediment, which reflects the properties of the sediment routing system. To model the proportion of dates in a child sample that comes from each of the parent distributions, we use a Bayesian mixture of Dirichlet processes. This model enables us to estimate the mixing proportions with associated uncertainty while making minimal assumptions. We also present an extension to the model whereby we reconstruct unobserved parent distributions from multiple observed child distributions using mixtures of Dirichlet processes. The model accounts for uncertainty in both the number of mineral formation events that constitute each parent distribution and the mixing proportions of each parent distribution that constitutes a child distribution. To demonstrate the model, we perform analyses using simulated data where the true age distribution is known as well as using a real-world case study from the coast of central California, USA.

Mathematical Geosciences↗

USGS2021

Identifying elusive piercing points along the North American transform margin using mixture modeling of detrital zircon data from sedimentary units and their crystalline sources

The San Gabriel and Canton faults represent early stages in the development of the San Andreas fault system. However, questions of timing of initiation and magnitude of slip on these structures remain unresolved, with published estimates ranging from 42-75 km and likely starting in the Miocene. This uncertainty in slip history reflects an absence of appropriate piercing points. We attempt to better constrain the slip history on these faults by quantifying the changing proportions of source terranes contributing sediment to the Ventura Basin, California, through the Cenozoic, including refining data for a key piercing point. Ventura Basin sediments show an increase in detrital zircon U-Pb dates and mineral abundances associated with crystalline sources in the northern San Gabriel Mountains through time, which we interpret to record the basin’s northwest translation by dextral strike-slip faulting. In particular, an Oligocene unit mapped as part of the extra-regional Sespe Formation instead has greater affinity to the Vasquez Formation. Specifically, the presence of a unimodal population of ~1180 Ma zircon, high (57%) plagioclase content, and proximal alluvial fan facies indicate that the basin was adjacent to the San Gabriel anorthosite during deposition of the Vasquez Formation, requiring 35-60 km of slip on the San Gabriel-Canton fault system. Mixture modeling of detrital zircon data supported by automated mineralogy highlights the importance of this piercing point along the San Gabriel-Canton fault system and suggests that fault slip began during the late Oligocene to early Miocene, which is earlier than published models. These two lines of evidence disagree with recent models that estimate >60 km of offset, requiring a reappraisal of the slip history of an early strand of the San Andreas transform zone.

California↗

USGS2021

Detrital zircon record of magmatism and sediment dispersal across the North American Cordilleran arc system (28-48°N)

As zircon U-Pb geochronology has become a leading method in sediment provenance studies and basin analysis over the past 20 years, the volume of detrital zircon data made available in published literature has enabled researchers to go beyond source-to-sink provenance studies to explore increasingly complex geologic problems. In this review, we utilize the growing body of detrital zircon data acquired from Jurassic-Paleocene forearc and foreland basin strata of the North American Cordillera to investigate the Mesozoic to earliest Cenozoic evolution of the arc and its associated basins between 28°N and 48°N. Our compilation includes 830 detrital zircon samples (101,898 individual ages) from 70 studies published between 2000 and 2020. For comparative purposes, we also compile 1307 igneous zircon U-Pb ages that characterize the magmatic history of the arc. We place primary emphasis on detrital zircon ages between 251 and 56 Ma that we infer to be uniquely derived from magmatic sources in the arc. Informed by existing knowledge of magmatic, structural, and sedimentological processes that acted on the orogen, we investigate spatial and temporal trends in these “arc-derived zircon” to establish a detrital record of arc magmatism, investigate source-to-sink relationships between the arc and adjacent basins, and discuss controls on sediment dispersal across the orogen. Our review shows that compilations of detrital zircon data from the Cordilleran forearc and foreland basin systems are excellent proxies for arc magmatism because the basins are enriched in arc-derived zircon and compilations provide space- and time-integrated records of crystallization ages. The compiled detrital zircon data support a history of continuous arc magmatism throughout Mesozoic and earliest Cenozoic time, characterized by low-volume magmatism from Triassic-Early Jurassic time (~251–174 Ma) and episodic higher-volume magmatism from Middle Jurassic-Late Cretaceous time (~174–66 Ma). These trends elucidate the initiation and timing of magmatic events at the orogen-scale and corroborate our understanding of cyclic arc behavior. Detrital zircon distributions are spatially and temporally variable both within and across basins, which we discuss relative to topographic development of the orogen and attendant responses of sediment dispersal systems. Detrital zircon distributions in the forearc signal rapid transfer of sediment from the arc to basins dominantly via fluvial processes. In contrast, detrital zircon distributions across the foreland reflect the presence of topographic barriers in the hinterland region of the arc that effectively isolated parts of the foreland. The presence of hinterland topography in turn highlights the important role of ash-fall events in delivering arc-derived zircon to the foreland, underscoring the need to consider ash-fall processes in paleodrainage reconstructions. These broad regional trends, and in general the close linkage between orogenic process and sediment dispersal, emerge from our compilation because it averages out much of the local variability observed in studies of more limited geographic or temporal extent.

North American Cordilleran arc system↗

USGS2020

The grass is not always greener on the other side: Seasonal reversal of vegetation greenness in aspect-driven semiarid ecosystems

Our current understanding of semiarid ecosystems is that they tend to display higher vegetation greenness on polar-facing slopes (PFS) than on equatorial-facing slopes (EFS). However, recent studies have argued that higher vegetation greenness can occur on EFS during part of the year. To assess whether this seasonal reversal of aspect-driven vegetation is a common occurrence, we conducted a global scale analysis of vegetation greenness on a monthly time scale over an 18-year period (2000-2017). We examined the influence of climate seasonality on the normalised difference vegetation index (NDVI) values of PFS and EFS at 60 different catchments with aspect-controlled vegetation located across all continents except Antarctica. Our results show that an overwhelming majority of sites (70%) display seasonal reversal, associated with transitions from water-limited to energy-limited conditions during wet winters. These findings highlight the need to consider seasonal variations of aspect-driven vegetation patterns in ecohydrology, geomorphology, and earth system models.

Geophysical Research Letters↗

USGS2020

Timing of Cenozoic extension in the southern Stillwater Range and Dixie Valley, Nevada

The Dixie Valley fault bounds the east side of the Stillwater Range in west‐central Nevada and last ruptured in 1954. Offset basalts indicate that slip began more recently than ~14 Ma, and prior work has interpreted the southern segment as an active low‐angle normal fault. Oligocene igneous rocks in the southern Stillwater Range were steeply tilted during large‐magnitude extension prior to ~14 Ma. To refine the timing of early extension and the onset of slip on the Dixie Valley fault, we collected two transects of samples for apatite fission track, apatite and zircon (U‐Th)/He (AHe and ZHe), and apatite 4 He/ 3 He thermochronometry. Apatite fission track ages from the Oligocene IXL pluton indicate rapid cooling ~18–14 Ma, and AHe and ZHe ages from the Cretaceous La Plata Canyon pluton indicate rapid cooling ~16–19 Ma. We interpret these data to record cooling during rapid extension. AHe ages from the IXL pluton are ~6–8 Ma and record cooling during slip on the Dixie Valley fault. We modeled these ages and 4 He/ 3 He spectra from one sample as the result of cooling during exhumation of a tilted fault block at a constant extension rate. The model predicts slip on the Dixie Valley fault beginning ~8 Ma. Although it does not constrain the initial fault dip, the model illustrates how a low‐angle fault requires a higher extension rate to reproduce cooling ages. Consequently, we prefer a high‐angle southern Dixie Valley fault for strain compatibility with the high‐angle northern segment.

Nevada↗

USGS2020

Coexisting seismic behavior of transform faults revealed by high-resolution bathymetry

Transform faults are known to have anomalously low rates of seismicity, but no direct observations reveal why this is the case. We use new, autonomous underwater vehicle high-resolution seafloor mapping to image the morphology of and offsets along transform fault segments in the Gulf of California. Fault splays display a varied history of activation and deactivation of individual fault strands over time, not unlike those mapped onshore or imaged within the bathymetry of the Queen Charlotte-Fairweather and the Palos Verdes faults of offshore western Canada and Southern California. A series of six identically offset depositional fans evidence 21–23 meters of slip along the main transform fault, which could not have been produced by a single earthquake. Rather, the lack of smaller-magnitude offsets indicates synchronous deposition and an absence of multiple slope failure-inducing earthquakes, thus providing the first direct evidence that creep and earthquakes occur at different times in the slip history of a given transform fault segment.

Geology↗

USGS2020

Dilution and propagation of provenance trends in sand and mud: Geochemistry and detrital zircon geochronology of modern sediment from central California (U.S.A.)

Integrated, multi-method provenance studies of siliciclastic sedimentary deposits are increasingly used to reconstruct the history of source-to-sink transport, paleogeography, and tectonics. Invariably, analysis of large-scale depositional systems must confront issues regarding how to best sample the system and adequately cope with the details of sediment mixing. Potential biases including variations in grain size, sediment flux, and zircon concentration may cause provenance tracking tools to misrepresent the contributions of source-areas that contribute to large drainage networks. We have acquired U-Pb detrital zircon data from modern sand and whole rock geochemistry from mud sampled from the Sacramento-San Joaquin drainage of central California to elucidate conditions that can skew provenance trends along the course of a major river system. This drainage network is fed by headwaters that tap the Mesozoic pluton-dominated southern Sierra Nevada, the Paleozoic-Mesozoic wallrock and volcanic-dominated northern Sierra Nevada, the ultramafic-dominated eastern Klamath Mountains, and the intermediate to mafic Cascades volcanic arc. Analysis of the results indicates that detrital zircon provenance trends effectively record source variations for the southern, granite-dominated portion of the drainage network where contrasts in lithology and inferred zircon fertility are relatively minor. In these circumstances, mixture modeling of U-Pb detrital zircon data calibrated with a measure of zircon fertility approximates relative sediment flux contributed by individual drainages. Alternatively, in the northern parts of the system, source regions underlain by ultramafic and /or volcanic rocks are poorly represented, or entirely missing, in down-stream detrital zircon records. In some cases, mud geochemistry data more faithfully represents sediment provenance trends. Sampling performed at the confluence of the Sacramento, American, Mokolumne, and San Joaquin rivers within the Sacramento Delta region yields a detrital zircon age distribution that is indistinguishable from that of an independently established database of Sierra Nevada batholith crystallization ages. However, when the combined river flows along a recently established passage to the Pacific through the San Francisco Bay region, dredged sediment is found to be significantly contaminated by locally eroded material from the Franciscan Complex and other rocks that crop out within the Coast Ranges. Large variation of Zr concentrations measured throughout the Bay area document that significant hydrodynamic fractionation impacts sediment delivery through this segment of the system. The more Sierra Nevada-like detrital zircon age distribution yielded by a piston-core sample from the continental slope may be explained by either early-stage unroofing of the Coast Ranges or more efficient sand delivery from the delta to the Pacific by a free flowing river driven by a low stand in sea level.

California↗

USGS2019

A stratigraphic approach to inferring depositional ages from detrital geochronology data

With the increasing use of detrital geochronology data for provenance analyses, we have also developed new constraints on the age of otherwise undateable sedimentary deposits. Because a deposit can be no older than its youngest mineral constituent, the youngest defensible detrital mineral age defines the maximum depositional age of the sampled bed. Defining the youngest `defensible' age in the face of uncertainty (e.g., analytical and geological uncertainty, or sample contamination) is challenging. The current standard practice of finding multiple detrital minerals with indistinguishable ages provides confidence that a given age is not an artifact; however, we show how requiring this overlap reduces the probability of identifying the true youngest component age. Barring unusual complications, the principle of superposition dictates that sedimentary deposits must get younger upsection. This fundamental constraint can be incoporated into the analysis of depositional ages in sedimentary sections through the use of Bayesian statistics, allowing for the inference of bounded estimates of true depositional ages and uncertainties from detrital geochronology so long as some minimum age constraints are present. We present two approaches for constructing a Bayesian model of deposit ages, first solving directly for the ages of deposits with the prior constraint that the ages of units must obey stratigraphic ordering, and second describing the evolution of ages with a curve that represents the sediment accumulation rate. Using synthetic examples we highlight how this method preforms in less-than-ideal circumstances. In an example from the Magallanes Basin of Patagonia, we demonstrate how introducing other age information from the stratigraphic section (e.g., fossil assemblages or radiometric dates) and formalizing the stratigraphic context of samples provides additional constraints on and information regarding depositional ages or derived quantities (e.g., sediment accumulation rates) compared to isolated analysis of individual samples.

Frontiers in Earth Science↗

USGS2019

A supervolcano and its sidekicks: A 100 ka eruptive chronology of the Fish Canyon Tuff and associated units of the La Garita magmatic system

Establishing temporal constrains on major volcanic eruptions is limited by the precision of existing geochronometers. Prior work on the La Garita caldera, created by the eruption of the Fish Canyon Tuff, failed to resolve temporal differences between pre-, syn-, and post-collapse eruptive units. Here, we report 40Ar/39Ar geochronologic data supporting a ca. 100 ka eruptive history of the La Garita caldera, and resolving the timing of the precaldera Pagosa Peak Dacite , syncaldera Fish Canyon Tuff, and postcaldera Dacite of Nutras Creek. Minimizing uncertainty in neutron fluence by rotating samples during irradiation and employing Bayesian statistical interpretation of analytical data enables resolution of the ca. 60 ka precaldera eruptive history and a hiatus of 0 - 20 ka prior to the eruption of postcaldera lavas. The improved precision demonstrated using these methods provides previously unresolvable temporal constraints on physical processes in the La Garita magmatic system and underscores the potential of unraveling other closely-spaced events in geologic time.

Geology↗