USGS Science⌕ Search

SEARCH · USGS Science

Results for “The R Journal”

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.

At least 235 records · Page 13Linked to original sources
USGS2019

Offset channels may not accurately record strike-slip fault displacement: Evidence from landscape evolution models

Slip distribution, slip rate, and slip per event for strike‐slip faults are commonly determined by correlating offset stream channels—under the assumption that they record seismic slip—but offset channels are formed by the interplay of tectonic and geomorphic processes. To constrain offset channel development under known tectonic and geomorphic conditions, we use numerical landscape evolution simulations along a theoretical strike‐slip fault with uniform and steady uplift, erosion, and diffusion. We investigate the influence of four tectonic parameters (fault zone width, earthquake recurrence interval, variance of the recurrence interval, and total slip relative to channel spacing) on offset channel development through multiple earthquake cycles. Analysis of >3,000 automatically measured offsets from >135 simulations suggests ~30% variability in individual measurements, but modeled displacement is recovered by averaging multiple measurements. However, the average of multiple offset measurements systematically underestimates modeled slip except when the fault zone is less than ~5 m wide, total slip is less than channel spacing, and offsets are measured shortly after an earthquake. In these simulations, postearthquake landscape evolution widens the geomorphic expression of the fault zone and modifies apparent channel offsets. We distinguish this “geomorphic fault zone” from the tectonic fault zone (zone of coseismic distributed deformation). This study highlights the capability of landscape evolution models to explore a range of conditions not easily defined in natural examples and the importance of averaging multiple measurements. Our results verify that paleoseismic studies must consider how geomorphic change has modified offset markers and use caution interpreting slip histories with multiple earthquakes.

Journal of Geophysical Research Solid Earth↗
USGS2021

The birth of a Hawaiian fissure eruption

Most basaltic explosive eruptions intensify abruptly, allowing little time to document processes at the start of eruption. One opportunity came with the initiation of activity from fissure 8 (F8) during the 2018 eruption on the lower East Rift Zone of Kīlauea, Hawaii. F8 erupted in four episodes. We recorded 28 min of high‐definition video during a 51‐min period, capturing the onset of the second episode on 5 May. From the videos, we were able to analyze the following in‐flight parameters: frequency and duration of explosions; ejecta heights; pyroclast exit velocities; in‐flight total mass and estimated mass eruption rates; and the in‐flight total grain size distributions. The videos record a transition from initial pulsating outgassing, via spaced, but increasingly rapid, discrete explosions, to quasisustained, unsteady fountaining. This transition accompanied waxing intensity (mass flux) of the F8 eruption. We infer that all activity was driven by a combination of the ascent of a coupled mixture of small bubbles and melt, and the buoyant rise of decoupled gas slugs and/or pockets. The balance between these two types of concurrent flow determined the exact form of the eruptive activity at any point in time, and changes to their relative contributions drove the transition we observed at early F8. Qualitative observations of other Hawaiian fountains at Kīlauea suggest that this physical model may apply more generally. This study demonstrates the value of in‐flight parameters derived from high‐resolution videos, which offer a rapid and highly time‐sensitive alternative to measurements based on sampling of deposits posteruption.

Island of Hawai'i↗
USGS2021

A rock record of complex aeolian bedforms in a Hesperian desert landscape: The Stimson formation as exposed in the Murray Buttes, Gale Crater, Mars

Lithified aeolian strata encode information about ancient planetary surface processes and the climate during deposition. Decoding these strata provides insight regarding past sediment transport processes, bedform kinematics, depositional landscape, and the prevailing climate. Deciphering these signatures requires a detailed analysis of sedimentary architecture to reconstruct dune morphology, motion, and the conditions that enabled their formation. Here, we show that a distinct sandstone unit exposed in the foothills of Mount Sharp, Gale crater, Mars, records the preserved expression of compound aeolian bedforms that accumulated in a large dune field. Analysis of Mastcam images of the Stimson formation shows that it consists of cross‐stratified sandstone beds separated by a hierarchy of erosive bounding surfaces formed during dune migration. The presence of two orders of surfaces with distinct geometrical relations reveals that the Stimson‐era landscape consisted of large dunes (draas) with smaller, superimposed dunes migrating across their lee slopes. Analysis of cross‐lamination and subset bounding surface geometries indicate a complex wind regime that transported sediment toward the north, constructing oblique dunes. This dune field was a direct product of the regional climate and the surface processes active in Gale crater during the fraction of the Hesperian Period recorded by the Stimson formation. The environment was arid, supporting a large aeolian dune field; this setting contrasts with earlier humid depositional episodes, recorded by the lacustrine sediments of the Murray formation (also Hesperian). Such fine‐scale reconstruction of landscapes on the ancient surface of Mars is important to understanding the planet’s past climate and habitability.

Journal of Geophysical Research - Planets↗
USGS2022

Evidence for fluctuating wind in shaping an ancient Martian dune field: The Stimson formation at the Greenheugh pediment, Gale crater

Temporal fluctuations of wind strength and direction can influence aeolian bedform morphology and orientation, which can be encoded into the architecture of aeolian deposits. These strata represent a direct record of atmospheric processes and can be used to understand ancient Martian atmospheric processes as well as those on Earth. The strata can: give insight to ancient atmospheric circulation, how the atmosphere evolved in response to global changes in habitability, and how ancient processes differ from modern processes. The Stimson formation at the Greenheugh pediment (Gale crater) records evidence of fluctuating wind across multiple temporal scales. The strata can be subdivided into three intervals–Gleann Beag, Ladder, and Edinburgh intervals. Internally, the intervals record changes of dune morphology and orientation, correlatable to wind fluctuations at multiple temporal scales. The basal Gleann Beag interval comprises compound cross-strata, deposited by oblique compound dunes. These dunes record a bimodal wind regime, resulting in net sediment transport toward the north. The Ladder interval records a reversal of sediment transport to the south, where straight-crested simple-dunes shaped by a seasonally variable winds formed. Finally, the Edinburgh interval records sediment transport to the west, where a unimodal wind formed sinuous-crested simple dunes. These observations demonstrate active and variable atmospheric circulation in Gale crater during the accumulation of the Stimson dune field, at multiple temporal scales from seasonally driven winds to much longer time-frames, during the Hesperian. These observations can be used to further understand ancient atmospheric conditions and processes, at a high temporal resolution on Mars.

Journal of Geophysical Research - Planets↗
USGS2022

Atmospheric river storm flooding influences tidal marsh elevation building processes

Disturbances are a key component of ecological processes in coastal ecosystems. Investigating factors that affect tidal marsh accretion and elevation change is important, largely due to accelerating sea-level rise and the ecological and economic value of wetlands. Sediment accumulation rates, elevation change, and flooding were examined at five marshes along a riverine-tidal gradient in the northern San Francisco Bay-Delta, California, USA during an Atmospheric River storm event in 2017 using Surface Elevation Tables (SETs), feldspar marker horizons (MH), and continuous water-level sensors. Our results showed that localized marsh flooding increased during the storm event, but not evenly across sites. Marsh surface elevation increased the most at the tidal freshwater marsh site in response to the storms, with an average surface elevation gain of 45.6 ± 13.1 mm, and the least at a tidal saline marsh with an average surface elevation gain of 4.0 ± 1.2 mm. A marsh located on the large embayment did not exhibit an immediate response to the storm but had a surface elevation gain of 21.5 ± 13.7 mm 6 months after the storm. During the storm period, marsh distance to the bay was the strongest predictor of elevation change, followed by SET-MH elevations. Conversely, during non-storm periods, SET-MH elevation was a relatively strong predictor of elevation change. Atmospheric Rivers appear to be a major factor affecting short-term spatial and temporal variability in flooding and sedimentation rates in tidal marsh systems. Incorporating information about storms into monitoring could increase our understanding of how episodic storms can impact marshes.

Journal of Geophysical Research: Biogeosciences↗
USGS2022

Quantifying large-scale surface change using SAR amplitude images: Crater morphology changes during the 2019-2020 Shishaldin Volcano eruption

Morphological processes often induce meter-scale elevation changes. When a volcano erupts, tracking such processes provides insights into the style and evolution of eruptive activity and related hazards. Compared to optical remote-sensing products, synthetic aperture radar (SAR) observes surface change during inclement weather and at night. Differential SAR interferometry estimates phase change between SAR acquisitions and is commonly applied to quantify deformation. However, large deformation or other coherence loss can limit its use. We develop a new approach applicable when repeated digital elevation models (DEMs) cannot be otherwise retrieved. Assuming an isotropic radar cross-section, we estimate meter-scale vertical morphological change directly from SAR amplitude images via an optimization method that utilizes a high-quality DEM. We verify our implementation through simulation of a collapse feature that we modulate onto topography. We simulate radar effects and recover the simulated collapse. To validate our method, we estimate elevation changes from TerraSAR-X stripmap images for the 2011–2012 eruption of Mount Cleveland. Our results reproduce those from two previous studies; one that used the same dataset, and another based on thermal satellite data. By applying this method to the 2019–2020 eruption of Shishaldin Volcano, Alaska, we generate elevation change time series from dozens of co-registered TerraSAR-X high-resolution spotlight images. Our results quantify previously unresolved cone growth in November 2019, collapses associated with explosions in December–January, and further changes in crater elevations into spring 2020. This method can be used to track meter-scale morphology changes for ongoing eruptions with low latency as SAR imagery becomes available.

Alaska↗
USGS2022

Toward next-generation lava flow forecasting: Development of a fast, physics-based lava propagation model

During effusive volcanic crises, the eruption and propagation of lava flows pose a significant hazard to nearby populations, homes, and infrastructure. Consequently, timely lava flow forecasts are a critical need for volcano observatory and emergency management operations. Previous lava flow modeling tools are typically either too slow to produce timely forecasts, or are fast, but lack critical aspects of lava physics or important forecasting outputs. In particular, the strong thermal stratification present in laminar, high-Prandtl number flows has generally been neglected. Bulk rheological changes have previously been computed from cell-averaged temperatures, assuming that the flow is thermally mixed. Here, we detail the development and initial testing of Lava2d, a new two-dimensional depth-averaged finite volume model of lava flow propagation over natural terrain which accounts for bulk rheological changes due to thermorheological stratification. We use a novel approach to energy conservation based on tracking cooling and solidifying at the flow base and at the moving flow surface, allowing for the estimation of more realistic vertical thermorheological profiles, while maintaining computational efficiency, producing very timely model runs. We validate our approach with three examples: comparison with theoretical propagation of crust-dominated lava flows, comparison with a large-scale molten basalt experiment from the Syracuse University Lava Project, and efficiency testing and comparison with the initial phase of the 1984 Mauna Loa lava flows. Our model is shown to produce rapid, realistic forecasts, making it a good candidate for operationalization in active volcanic regions such as in Hawai'i.

Journal of Geophysical Research - Solid Earth↗
USGS2023

Using corrected and imputed polarity measurements to improve focal mechanisms in a regional earthquake catalog near the Mt. Lewis Fault Zone, California

We utilized relative polarity measurements and machine learning techniques to better resolve focal mechanisms and stress orientations considering a catalog of ∼29,000 relocated earthquakes that occurred during 1984–2021 in the southeastern San Francisco Bay Area. Earthquake focal mechanisms are commonly produced using P wave first motion polarities, which traditionally requires events to be well-recorded across a seismic network with good focal sphere coverage. We adapted recently developed approaches that are less dependent on high signal-to-noise records and exploit similar waveforms to produce relative polarity and amplitude measurements between earthquake pairs. These techniques were previously only applied on localized earthquake sequences, and we further developed these approaches so that they can be utilized for regional catalogs. We validated or corrected manually identified polarities by performing polarity consensuses using earthquake pairs. Missing and unreliable polarity measurements were imputed using iterative random forests, an unsupervised ensemble machine learning method. Relative P and S wave amplitude measurements were made between earthquakes, constraining S / P ratios for low signal-to-noise waveforms. Using these techniques, we were able to reduce focal mechanism uncertainties by an average of ∼13° and produced well-constrained focal mechanisms for ∼6 times as many earthquakes than those produced using only the traditionally derived polarities. We performed stress inversions using the focal mechanisms by grouping the focal mechanism results into a quadtree structure. Our stress results are consistent with previous work, albeit at a higher spatial resolution, and demonstrate these techniques can aid our understanding of fault structures and kinematics in more detail than was previously possible.

California↗
USGS2022

Sedimentological and geochemical perspectives on a marginal lake environment recorded in the Hartmann’s Valley and Karasburg members of the Murray formation, Gale crater, Mars

This study utilizes instruments from the Curiosity rover payload to develop an integrated paleoenvironmental and compositional reconstruction for the 65-m thick interval of stratigraphy comprising the Hartmann's Valley and Karasburg members of the Murray formation, Gale crater, Mars. The stratigraphy consists of cross-stratified sandstone (Facies 1), planar-laminated sandstone (Facies 2), and planar-laminated mudstone (Facies 3). Facies 1 is composed of sandstone showing truncated sets of concave-curvilinear laminae stacked into cosets. Sets are estimated to be meter-to sub-meter-scale, consistent with low-height dunes. Thin stratigraphic intervals of Facies 1 and stacking patterns with Facies 2 and 3 support a wet aeolian dune interpretation. Meter-thick packages of planar-laminated sandstone (Facies 2) are interpreted to represent interfingering dune-interdune strata. Facies 3 consists of meter-thick packages of planar-laminated mudstone interpreted to represent lacustrine deposition with persistent standing water. Integration of geochemistry with each facies reveals some compositional control based on the depositional process. Models for source rock composition from Alpha Particle X-Ray Spectrometer measurements show that facies derived from a basaltic source. Alteration indices and geochemical trends provide evidence that moderate chemical weathering occurred before compositional changes due to diagenesis. Differences in wt% FeO (T) and TiO 2 between facies are minimal, though trends point to sediment sorting in transport. Comparisons to terrestrial basaltic sedimentary systems indicate that the Hartmann's Valley and Karasburg facies reflect deposition in an environment where diverse subaqueous and subaerial facies persisted adjacent to a long-lived body of water.

Journal of Geophysical Research Planets↗
USGS2023

Variation in leaf reflectance spectra across the California flora partitioned by evolutionary history, geographic origin, and deep time

Evolutionary relatedness underlies patterns of functional diversity in the natural world. Hyperspectral remote sensing has the potential to detect these patterns in plants through inherited patterns of leaf reflectance spectra. We collected leaf reflectance data from across the California flora from plants grown in a common garden. Regions of the reflectance spectra vary in the depth and strength of phylogenetic signal. We also show that these differences are much greater than variation due to the geographic origin of the plant. At the phylogenetic extent of the California flora, spectral variation explained by the combination of ecotypic variation (divergent evolution) and convergent evolution of disparate lineages was minimal (3 to 7 %) but statistically significant. Interestingly, at the extent of a single genus ( Arctostaphylos ) no unique variation could be attributed to geographic origin. However, up to 18% of the spectral variation among Arctostaphylos individuals was shared between phylogeny and intraspecific variation stemming from ecotypic differences (i.e., geographic origin). Future studies could conduct more structured experiments (e.g., transplants or observations along environmental gradients) to disentangle these sources of variation and include other intraspecific variation (e.g., plasticity). We constrain broad scale spectral variability due to ecotypic sources, providing further support for the idea that phylogenetic clusters of species might be detectable through remote sensing. Phylogenetic clusters could represent a valuable dimension of biodiversity monitoring and detection.

California↗
USGS2025

Advancing current understanding of Martian impact-generated hydrothermal systems through novel coupled modeling: Insights from Gale, Jezero, and other craters

Impact-generated hydrothermal systems are considered potentially habitable environments on Mars, Earth, and other planetary bodies for microbial life. However, there is an ongoing debate regarding what geological features on Mars provide definitive evidence for such systems. Although earlier studies have modeled hydrothermal processes in Martian craters, they often lacked integration with shock physics hydrocodes to constrain initial impact conditions. The importance of this two-code coupling was demonstrated by successfully replicating alteration signatures in the Earth's Haughton impact structure. In this study, we use a similar two-code approach, combining the iSALE hydrocode with the HYDROTHERM hydrothermal model to simulate the full evolution of impact-generated hydrothermal systems. We apply this method to craters the size of Jezero (∼50 km) and Gale (∼154 km) in diameter. Although Jezero's interior is largely buried, our results align with hypothesized hydrothermal vents and alteration minerals near central uplifts in similarly sized exposed craters, such as Toro and Auki. Furthermore, our models correspond to alteration patterns observed by the Curiosity in the lower layers of Mount Sharp, which may represent remnants of impact-driven hydrothermal activity. A key finding is that these systems may persist much longer than previously estimated. Our simulations suggest that a Jezero-sized system could remain habitable for thermophiles for approximately 720,000 years, whereas a Gale-sized system could persist for nearly 2 million years. Additionally, simulations under unsaturated crustal conditions reveal that air-dominated near-surface layers can suppress vertical fluid flow, enabling deep subsurface alteration without producing detectable mineral signatures at the surface.

Journal of Geophysical Research - Planets↗
USGS2024

Constraining mean landslide occurrence rates for non-temporal landslide inventories using high-resolution elevation data

Constraining landslide occurrence rates can help to generate landslide hazard models that predict the spatial and temporal occurrence of landslides. However, most landslide inventories do not include any temporal data due to the difficulties of dating landslide deposits. Here we introduce a method for estimating the mean landslide occurrence rate of deep-seated rotational and translational slides derived solely from high-resolution (≤3 m) elevation data and globally available estimates of the diffusion coefficient for sediment flux. The method applies a linear diffusion model to the roughest landslide deposits until they reach a representative non-landslide roughness distribution. This estimates the time for a landslide deposit to be unrecognizable in high-resolution digital elevation data, which we term the mean lifetime of the landslide. Using the mean lifetime and number of landslides within an area of interest, we can estimate the mean occurrence rate of landslides over that domain. We validate this approach using a comprehensive temporal inventory of landslides in western Oregon created using age-roughness curves that are calibrated with high-resolution elevation data and radiocarbon data. We find good agreement between our diffusion method and the existing age-roughness-derived estimates, producing mean lifetimes of 4500 and 5200 years (4% difference), respectively. Hazard maps produced using the two methodologies generally agree, with the maximum differences in landslide probability reaching 0.1. Due to the relative abundance of high-resolution elevation data compared with age-dated landslides, our method could help constrain landslide occurrence rates in areas previously considered unfeasible.

Oregon↗
USGS2012

Fractal properties of tremor and gas piston events observed at Kilauea Volcano, Hawaii

Studies the fractal properties of shallow volcanic tremor and gas piston events associated with magma degassing at Kilauea Volcano, Hawaii, using data from two dense short-baseline arrays of seismographs deployed near the active crater of Puu Oo on the east rift of the volcano. The existence of a categorically stable attractor characterizing both types of activities strongly suggests that the excitation mechanism of tremor is similar to that of gas piston events, which in turn are correlated with visual observations at the volcanic vent. Low values derived for the dimension of the attractor in phase space points to significant self-organization in the process of generation of tremor and offers general constraints on the dimensionality of attractors derived from models of acoustic emission associated with magma flow, vesiculation, and degassing. -from Authors

Hawaii↗
USGS2012

Limitations of quantitative analysis of deep crustal seismic reflection data: Examples from GLIMPCE

Amplitude preservation in seismic reflection data can be obtained by a relative true amplitude (RTA) processing technique in which the relative strength of reflection amplitudes is preserved vertically as well as horizontally, after compensating for amplitude distortion by near-surface effects and propagation effects. Quantitative analysis of relative true amplitudes of the Great Lakes International Multidisciplinary Program on Crustal Evolution seismic data is hampered by large uncertainties in estimates of the water bottom reflection coefficient and the vertical amplitude correction and by inadequate noise suppression. Processing techniques such as deconvolution, F-K filtering, and migration significantly change the overall shape of amplitude curves and hence calculation of reflection coefficients and average reflectance. Thus lithological interpretation of deep crustal seismic data based on the absolute value of estimated reflection strength alone is meaningless. The relative strength of individual events, however, is preserved on curves generated at different stages in the processing. We suggest that qualitative comparisons of relative strength, if used carefully, provide a meaningful measure of variations in reflectivity. Simple theoretical models indicate that peg-leg multiples rather than water bottom multiples are the most severe source of noise contamination. These multiples are extremely difficult to remove when the water bottom reflection coefficient is large (>0.6), a condition that exists beneath parts of Lake Superior and most of Lake Huron.

Great Lakes↗
USGS2012

Interseismic deformation at the Nankai Trough, Japan, subduction zone

Uplift along the coast of southwestern Japan following the 1944 Tonankai ( M s = 8.0) and 1946 Nankaido ( M s = 8.2) earthquakes has been inferred from the 1950–1985 tide gage records. Although uplift rates vary with geographic position, the temporal dependence at each site can be described as the superposition of an exponentially decaying (time constant ∼5 years) transient, significant only in the first decade postseismic, and a steady interseismic rate. The steady interseismic uplift rate is defined only over the 1956–1985 interval, but the linearity in time ofthat uplift record is such that uplift over the remainder of the interseismic cycle (∼80 years) might reasonably be expected to occur at the same rate. Moreover, the spatial pattern and magnitude of the steady uplift rate are reasonably close to those predicted by a simple dislocation model of subduction. However, the currently observed steady uplift rate is not consistent with the uplift rate inferred from geodetic leveling for the final half of the previous earthquake cycle. Thus, either the hypothesis of a steady interseismic uplift rate after the first decade postseismic is wrong or the strain accumulation pattern for one earthquake need not be the same as for the preceding earthquake on the same segment.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Paleomagnetism and geochronology of 23 Ma gabbroic intrusions in the Keku Strait, Alaska, and implications for the Alexander Terrane

Samples of Tertiary gabbro from 24 sites in the Keku Strait, Alaska, help constrain the displacement history of the Alexander terrane. Step heating experiments on a plagioclase separate from these previously undated intrusions indicate a discordant 40 Ar/ 39 Ar age of 23.1 ± 1.7 Ma. The characteristic magnetization resides in magnetite, is easily isolated by thermal and alternating field demagnetization, and has both normal and reversed polarities. The mean paleomagnetic pole, with no structural correction, is latitude 87.1°N, longitude 141.6°E, A 95 = 10.1°, with N = 20 sites. This pole indicates insignificant tectonic displacement (0.5° ± 8.2° southward) and rotation (0.6° ± 15.2° counterclockwise). We therefore conclude that any northward displacement or vertical axis rotation of the Alexander terrane, and/or tilting in the vicinity of the Keku Strait must have occurred before 23 Ma.

Alaska↗
USGS2012

Seasonal fluctuations in sea level on the South Carolina shelf and their relationship to the Gulf Stream

Near-bottom pressure measurements obtained over a 9-month period in 1978 on the outer continental shelf off South Carolina document that sea level rose 35 to 50 cm between July and October. Records of coastal sea level showed a similar rise. When the pressure and coastal sea level records were corrected for the effects of wind stress and temperature, the largest portion of the rise in the two variables remained. The amplitude of the change in residual sea level was largest at the shelf break, 28 cm, and decayed to 22 cm at the coast. The relatively small magnitude of the decay suggests that alongshelf pressure gradients imposed at the shelf break have minimum alongshelf wavelengths of 2000 km. The changes in residual sea level were seasonal and associated with seasonal changes in the transport of the Gulf Stream. The 3 years of residual sea level records examined here indicate that the amplitude of a typical decrease in transport of the Gulf Stream off South Carolina between July and October has a strong interannual variability. The change in transport can be 50% lower to over 100% higher than the historical average. The residual sea level records also suggest that a change in transport of the Gulf Stream is not strongly related to the onshore/offshore position of the Gulf Stream.

South Carolina↗
USGS2012

Evolution of the caldera‐forming eruption at Crater Lake, Oregon, indicated by component analysis of lithic fragments

Crater Lake caldera (8 × 10 km), formed 6845 years B. P. ( 14 C age) during the climactic eruption of the volcanic edifice known as Mount Mazama, is intermediate in size between small calderas associated with central vent eruptions and large calderas that have ring fracture vent systems. Our quantitative study of lithic fragments in the ejecta confirms the existing model of changes in vent configuration during the climactic eruption of Mount Mazama. Initial activity was from a single vent that produced a rhyodacite pumice fall from a Plinian column. Altered preexisting volcanic rocks are the predominant lithic type in the Plinian deposit, and their extensive hydrothermal alteration is considered as evidence of their relatively deep origin. The Wineglass Welded Tuff lies atop the Plinian deposit and contains a higher proportion of fresh volcanic rocks, suggesting enlargement of the single vent by slumping of its walls. This same vent enlargement caused the Plinian eruption column to collapse and feed valley‐hugging pyroclastic flows that deposited the Wineglass Welded Tuff. When enough material was erupted from the shallow magma chamber that its roof was no longer adequately supported, Mount Mazama collapsed to form the caldera, while highly energetic pyroclastic flows produced the climactic ignimbrite. A lag breccia that represents the proximal facies of the compositionally zoned climactic ignimbrite lies atop the Wineglass Welded Tuff and contains predominantly altered volcanic rocks of deeper origin, accompanied by minor granitoids from the magma chamber walls. Azimuthal differences in lithic component proportions in the lag breccia correlate well with the geology of the caldera walls, indicating that the climactic ignimbrite was ejected by multiple vents along a ring fracture system. Systematic lithic component changes within the lag breccia suggest different quarrying levels that reflect waxing and waning of the discharge rate during the caldera collapse phase of the climactic eruption. Our lithic component analysis demonstrates that calderas that may be too small to experience structural resurgence, such as Crater Lake, nevertheless may form by syneruptive subsidence along ring fractures.

Oregon↗