Geology topics

Mark E. Stelten

Publications and source records attributed to Mark E. Stelten.

At least 19 recordsLinked to original sources

USGS2026

Ignimbrites of the Sour Creek dome, Yellowstone volcano, USA: Implications for the Lava Creek Tuff eruption and Yellowstone Caldera

Field observations plus sanidine compositional and 40 Ar/ 39 Ar age data collected on rocks from the Sour Creek dome (SCD) area of Yellowstone volcano (United States) serve to identify five new ignimbrite packages with ages (average 631.5 ± 0.9 ka) that are analytically indistinguishable from those of the Lava Creek Tuff (LCT), challenging the current framework of two LCT eruptive members. Newly recognized package 1 occurs as LCT-aged recycled clasts within a proximal lag breccia, reflecting a major time break in the LCT eruption sequence. Conformably above, the scoria-bearing ignimbrite of package 2 is extensively found in the dome, having been previously mis-mapped as Huck-leberry Ridge Tuff and LCT member A. To the south, package 3 is a crystal-rich, densely welded ignimbrite (inferred to be younger than package 2) that is conformably overlain by package 4, a fine-grained ignimbrite that contains abundant clasts of recycled tuff. Lastly, the topographically highest parts of the dome are made up of package 5, a multicomponent ignimbrite containing scoria, recycled clasts, and pumices, potentially correlative with package 2 and/or package 4. Cathodoluminescence imaging and zoning assessment was integrated with sanidine geochemistry from all five packages, including their recycled clasts. Four of the five packages were sourced from separate magmatic bodies within a larger plumbing system, each defined by their sanidine Ba distributions and componentry. Although the SCD was originally mapped as a resurgent dome and is actively deforming today, faults cutting the region show comparable offsets for NW-SE “resurgent faulting” as for orthogonal NE-SW faulting, the latter of which also cuts the 162 ka Elephant Back lava flow. We infer that faulting in the SCD well postdates eruption and that the SCD instead represents a domical accumulation of ignimbrite, not a resurgent structure. The abundance of recycled materials (and absence of Absaroka country rocks) in the SCD ignimbrites implies that vent locations were well inboard of the mapped caldera boundary, with some potentially along the line of the Yellowstone River, requiring a revision of the currently accepted LCT caldera outline. The LCT eruptive events were thus far more complex than currently recognized, with multiple deposits fed from multiple magma bodies over a time scale that cannot be resolved by radiometric dating.

Wyoming↗

USGS2026

Dismembering the Lava Creek Tuff, Yellowstone volcano

The youngest large explosive eruption of Yellowstone volcano generated the 0.63 Ma Lava Creek Tuff (LCT) ignimbrite and a voluminous, widespread fall deposit accompanying formation of Yellowstone Caldera. The tuff is conventionally divided into members A (LCT-A) and B (LCT-B), defined by a reversal in density/welding profiles and the presence of phenocrystic amphibole in LCT-A. Remapping within the northeastern caldera sector has, however, documented multiple ignimbrites that yield LCT eruption ages but are distinct from documented members A and B. Thus, the framework by which the LCT was originally subdivided is re-examined here using field observations, petrologic analyses, and sanidine geochemistry from LCT-A and LCT-B ignimbrite (as previously mapped) at reference localities around the caldera. Mineral modal abundances in both LCT-A and LCT-B outcrops are dominated by quartz, followed by sanidine, with sizes and abundances that generally increase up-section, as previously documented. Amphibole is found by us, however, only in western and southern sectors, in deposits mapped as both LCT-A and LCT-B, and only as a trace phase. Further, a locality reported as LCT-A with amphibole present yields a sanidine 40 Ar/ 39 Ar date of 0.775 ± 0.002 Ma, precluding any LCT correlation. Ignimbrite samples from mapped LCT-A and LCT-B exhibit variable crystal sizes, sanidine cathodoluminescence image populations and Ba concentrations, with no clear distinction between the mapped members. Collectively this research found no evidence for a definable LCT-A versus LCT-B, meaning that the traditional A and B subdivision does not accurately represent the eruptive architecture of the Lava Creek Tuff and warrants revision.

Idaho, Montana, Wyoming↗

USGS2025

Refined chronology of late Quaternary eruptions at Harrat Khaybar, Saudi Arabia, with implications for magma dynamics and regional volcanic history

Determining accurate and precise ages for Quaternary volcanic centers is essential for reconstructing volcanic field histories, understanding magmatic processes, and assessing potential hazards or risk. Harrat Khaybar, western Saudi Arabia, is one of the youngest and potentially most active volcanic fields on the Arabian plate, has been active since ca. 1.7 Ma, and is characterized by a spectrum of rock compositions ranging from predominantly alkalic basalt to trachyte and comendite. Previous work in Harrat Khaybar utilizing 40 Ar/ 39 Ar incremental heating geochronology to constrain morphological preservation and superpositional relationships bracketed the volcanic activity into broad age groups in intervals of ∼150 k.y., and the youngest and most compositionally evolved volcanoes, including Jabal Abyad, Jabal Bayda, and Jabal Qidr, were assigned to the age groups between ca. 300 ka and present. Herein, we establish a detailed chronology of prehistoric silicic and historical basaltic eruptions at central Harrat Khaybar using four independent eruption age determination techniques: zircon double-dating (ZDD), which combines 238 U- 230 Th disequilibrium or U-Pb with (U-Th)/He dating; zircon U-Pb dating; cosmogenic 3 He dating; and cosmogenic 36 Cl geochronology. These were employed to accurately date six volcanic centers, including the comenditic Jabal Abyad, Jabal Bayda, Jabal Ibayl, and Jabal Alhayyirah, the trachytic Jabal Aluthmor, and the basaltic Jabal Qidr. Additionally, our previously published 40 Ar/ 39 Ar ages have been recalculated using isochron intercept (nonatmospheric) 40 Ar/ 36 Ar for the trapped Ar component. Our new results reveal that zircon rims from Jabal Abyad and Jabal Bayda define isochron 238 U- 230 Th crystallization ages of 125 ± 4 ka and 144 ± 6 ka, concordant with ZDD eruption ages of 132 ± 4 ka and 149 ± 5 ka, respectively. Zircon U-Pb crystallization and (U-Th)/He eruption ages from Jabal Alhayyirah are concordant at 471 ± 14 ka and 458 ± 18 ka, respectively. Finally, zircons from the nearby Jabal Ibayl yield a U-Pb weighted mean crystallization age of 566 ± 16 ka concordant with the corresponding (U-Th)/He eruption age of 554 ± 12 ka, both of which are notably older than the previously proposed eruption ages of 300−150 ka. Recalculations of published 40 Ar/ 39 Ar ages for the youngest volcanoes at central Harrat Khaybar are now in excellent agreement with new geochronological data. Our new age data reveal several new insights into the development of Harrat Khaybar. It is now clear that the comenditic eruptions do not belong to the same eruptive phase and indicate an extended history during which comendites have episodically punctuated the basaltic volcanism since at least 600 ka. The data indicate that Jabal Ibayl and Jabal Alhayyirah represent separate older volcanic events, whereas the younger Jabal Abyad and Jabal Bayda volcanoes appear to be coeval, and their spatial proximity implies that they share a magmatic lineage and maybe a common plumbing system. Zircon age spectra of the comendites reveal obvious xenocrysts and antecrysts indicating assimilation of basement and plutonic progenitors, but otherwise they define broad unimodal populations of crystallization ages that overlap within error with the respective eruption ages. We interpret this to indicate that zircon crystallization continued up to the time of eruption. The cosmogenic ages of Jabal Aluthmor and Jabal Qidr reveal that these centers erupted as recently as ca. 2000 and 760 years ago, respectively. Broadly coeval young silicic and basaltic eruptions at northern Harrat Rahat, the harrat immediately south of Harrat Khaybar, may imply a shared geodynamic forcing between the two adjacent volcanic fields.

Harrat Khaybar↗

USGS2025

Origin and evolution of mafic volcanism associated with 3 m.y. of andesite production at the Goat Rocks volcanic cluster, southern Washington Cascade Range

More than 3 m.y. of mafic volcanism near the Goat Rocks volcanic cluster in the southern Washington Cascade Range, USA, lends insight into the evolution of basalts and the subarc mantle at a long-lived, major arc volcanic locus. We contribute field observations, 40 Ar/ 39 Ar dates, paleomagnetic directions, and bulk rock and mineral compositions to characterize nine mafic units that erupted in association with the Goat Rocks volcanic cluster. The time frame of mafic volcanism, ca. 3.6 Ma to 60 ka, encompasses the lifespan of the central volcanic cluster (3.1 Ma to 115 ka), with a lull from ca. 2.7 Ma to 1.4 Ma. A climactic period of voluminous mafic activity and far-traveled lava flows, including construction of the Hogback Mountain shield volcano, coincided with voluminous andesite eruptions from the central volcanic cluster. The basaltic rocks in the Goat Rocks area are calc-alkaline to barely tholeiitic and have high field strength element depletion relative to large-ion lithophile elements characteristic of calc-alkaline basalts (CAB) of the Cascade volcanic arc. Unlike at neighboring andesitic volcanic centers (Mounts Adams, St. Helens, and Rainier), no other mafic end members such as high-aluminum olivine tholeiite (HAOT) or intraplate-type basalt (IPB) are present at or near the Goat Rocks volcanic cluster, although some of the calc-alkaline basalts in this study have IPB-like affinities. The Goat Rocks mafic units exhibit two main temporal trends in composition: (1) the most primitive basalts erupted earlier, compared to less primitive and more evolved compositions later, and (2) high field strength element concentrations are higher in the younger basalt units relative to the oldest two. In contrast to these temporal trends, the mafic units define two compositional groups that recur through time, a low-Sr and a high-Sr group, each with distinct trace element and Sr and Nd isotope ratios. Although radiogenic isotope ratios are generally aligned with High Cascades CAB and HAOT, some extend toward IPB of Mount Adams and Simcoe Mountains volcanic field. Olivine-dominated crystal fractionation at shallow pressure from a small range of parent magma compositions accounts for much of the variation among the basalts and basaltic andesites. A high-pressure fractionation model is plausible for only one of the youngest basalt units (basalt of Walupt Lake volcano). Mafic recharge and crustal assimilation accounts for the incompatible-element enriched composition of basaltic andesites erupted during construction of the largest andesitic centers, further supporting sustained basalt mass flux and thermal energy driving andesite genesis. We model the most primitive members of the Goat Rocks mafic units as partial melts of successively less depleted mantle in time. Variable degrees of fluxing with fluids and melts from subduction explain the distinction between high-Sr and low-Sr groups. We propose that mantle metasomatism by ancestral subduction and fluid-flux melting is heterogeneously distributed through the local subarc mantle and played a greater role in the genesis of the high-Sr basalt group. The limited range of primitive basalt types around the Goat Rocks volcanic cluster contrasts with the much greater diversity of basalts throughout the southern Washington to northern Oregon Cascade arc. On the other hand, the central volcanic cluster encompasses nearly the entire diversity observed at neighboring composite volcanoes. In the case of the Goat Rocks area at least, and perhaps attributable to the entire region, this means that the genesis of diverse intermediate magmas is independent from and does not require vastly different parental basalt compositions.

California, Oregon, Washington↗

USGS2025

New 40Ar/39Ar eruption ages reveal an important temporal relationship between mafic and silicic volcanism in the Yellowstone Plateau volcanic field

The chronology of mafic eruptions and their temporal relation to rhyolitic volcanism in the Yellowstone Plateau volcanic field are poorly known, thereby limiting our understanding of the way(s) in which mafic magmatism drives rhyolitic activity. To address this, we measured 40 Ar/ 39 Ar eruption ages on 13 mafic samples collected from Henrys Fork Caldera (eastern Idaho, western United States), which represents a region of known volcanic activity immediately west of Yellowstone caldera for which the relationship to Yellowstone volcano’s most recent caldera-forming cycle remains unclear. Our new ages indicate that mafic activity was occurring throughout the Henrys Fork Caldera both leading up to and following the emplacement of the Lava Creek Tuff. Furthermore, these ages reveal that mafic volcanism in the Henrys Fork Caldera region occurred concurrently with second- and third-cycle rhyolite volcanism in and around Yellowstone caldera. Our new ages therefore provide unique and definitive evidence that the mafic magmatism of Henrys Fork Caldera played a critical role in the development of shallow-crustal rhyolitic magma chambers that ultimately fueled the large caldera-forming eruptions within the Yellowstone volcanic system.

Idaho↗

USGS2025

The progression of basaltic–rhyolitic melt storage at Yellowstone Caldera

Yellowstone Caldera is one of the largest volcanic systems on Earth, hosting three major caldera-forming eruptions in the past two million years, interspersed with periods of less explosive, smaller-volume eruptions 1 . Caldera-forming eruptions at Yellowstone are sourced by rhyolitic melts stored within the mid- to upper crust. Seismic tomography studies have suggested that a broad region of rhyolitic melt extends beneath Yellowstone Caldera, with an estimated melt volume that is one to four times greater than the eruptive volume of the largest past caldera-forming eruption, and an estimated melt fraction of 6–28 per cent 2 , 3 , 4 , 5 . Seismic velocity is strongly influenced by temperature, pressure and melt; however, magnetotelluric data are primarily sensitive to the presence of melt, making these data ideal for constraining volcanic systems. Here we utilize magnetotelluric data to model the resistivity structure of Yellowstone Caldera’s crustal magma reservoir and constrain the region’s potential for producing major volcanic eruptions. We find that rhyolitic melts are stored in segregated regions beneath the caldera with low melt fractions, indicating that the reservoirs are not eruptible. Typically, these regions have melt volumes equivalent to small-volume post-caldera Yellowstone eruptions. The largest region of rhyolitic melt storage, concentrated beneath northeast Yellowstone Caldera, has a storage volume similar to the eruptive volume of Yellowstone’s smallest caldera-forming eruption. We identify regions of basalt migrating from the lower crust, merging with and supplying heat to the northeast region of rhyolitic melt storage. On the basis of our analysis, we suggest that the locus of future rhyolitic volcanism has shifted to northeast Yellowstone Caldera.

Idaho, Montana, Wyoming↗

USGS2024

A far-traveled basalt lava flow in north-central Oregon, USA

Widely separated basalt lava-flow outcrops in north-central Oregon, USA, expose products of a single eruptive episode. A Pliocene lava flow, here informally termed the Tetherow basalt, issued from vents near Redmond, in the Deschutes basin of Oregon, as a plains-forming basalt now exposed in continuous outcrops northward for 60 km. A similar basalt crops out 47 km farther north, near Maupin, within what was then a slightly incised ancestral Deschutes River canyon. The northernmost outcrops of this lava flow lie on Fulton Ridge, in the Dalles basin, near the confluence of the Deschutes and Columbia Rivers. Complementary lines of evidence confirm these rocks are all from the same volcanic eruption. Outcrops in the Deschutes and Dalles basins are chemically similar high-titanium basalts, petrographically similar to each other and distinct from other lava flows in the area. Paleomagnetic directions from 11 scattered sites are similar and indistinguishable by various tests for a common mean. Three new 40 Ar/ 39 Ar ages indicate the Tetherow basalt eruption occurred between 5.5 Ma and 5.0 Ma, likely at ca. 5.2 Ma. The widely separated outcrops of this lava flow span 160–180 km along the ancestral Deschutes River and downstream Columbia River. The lava flow’s length and erupted volume of 15–20 km 3 are extraordinarily large in a non-flood-basalt setting. This lava flow provides a datum with which to describe regional physiographic history, assess incision rates, and infer tectonic history. Spanning different depositional basins, the Tetherow basalt is a useful chronologic and stratigraphic marker bed.

Oregon↗

USGS2023

Eruptive history within the vicinity of Al Madīnah in northern Harrat Rahat, Kingdom of Saudi Arabia

The northernmost part of the Harrat Rahat volcanic field contains early Pleistocene to Holocene mafic eruptive products within the vicinity of the city of Al Madīnah, Kingdom of Saudi Arabia. A detailed geologic investigation into the eruptive history of a 570 square kilometer (km 2 ) area covering Al Madīnah and the surrounding area has yielded 33 mapped Quaternary volcanic units consisting of lava flows, scoria cones, and shield volcanoes. These eruptive products consist of continental, intraplate alkalic and minor transitional basalts, hawaiites, and a single mugearite that were emplaced from at least 1,014±14 thousand years ago (ka) to a single Holocene eruption in 1256 C.E. Lava flows are generally 10 to 15 kilometers (km) long (but can reach 23 km long), 1 to 3 km wide, and at least 10 meters thick. Most of the mapped units erupted episodically between 400 and 340 ka and 180 and 100 ka. Despite small individual volumes (less than 1 cubic kilometers dense rock equivalent), each unit represents eruption of a distinct magma batch that was strongly influenced by clinopyroxene, olivine, and plagioclase fractionation. Some of these units are interpreted to have undergone magma mixing pre- and (or) syneruptively. Integrating eruption ages, geochemistry, and paleomagnetic data yields evidence that some eruptions were temporally and (or) spatially clustered. Aligned scoria cones and elongate vent edifices were constructed atop fissure vent systems that reflect the local stress field, which controls dike ascent through the middle and upper crust.

Professional Paper↗

USGS2023

Eruptive history of northern Harrat Rahat—Volume, timing, and composition of volcanism over the past 1.2 million years

Harrat Rahat, one of several large, basalt-dominated volcanic fields in the western part of the Kingdom of Saudi Arabia, is a prime example of continental, intraplate volcanism. Excellent exposure makes this an outstanding site to investigate changing volcanic flux and composition through time. We present 93 40 Ar/ 39 Ar ages and 6 36 Cl surface-exposure ages for volcanic deposits throughout northern Harrat Rahat that, integrated with a new geologic map, define 12 eruptive stages. Exposed volcanic deposits in the study area erupted less than 1.2 million years ago (Ma), and 214 of 234 identified eruptions occurred less than 570 thousand years ago (ka). Two eruptions were in the Holocene, including a historically described basaltic eruption in 1256 C.E. and a trachyte eruption newly recognized as Holocene (4.2±5.2 ka). An estimated approximately 82 cubic kilometers (km 3 ; dense rock equivalent) of volcanic products can be documented as having erupted since 1.2 Ma, though this is a lower limit because of concealment of deposits older than 570 ka. Over the last 570 thousand years (k.y.), the average eruption rate was 0.14 cubic kilometers per thousand years (km 3 /k.y.), but volcanism was episodic with periods alternating between low (0.04–0.06 km 3 /k.y.) and high (0.1–0.3 km 3 /k.y.) effusion rates. Before 180 ka, eruptions vented from the volcanic field’s dominant eastern vent axis and from a subsidiary, diffuse, western vent axis. After 180 ka, volcanism focused along the eastern vent axis, and the composition of volcanism varied systematically along its length from basalt dominated in the north to trachyte dominated in the south. We hypothesize that these compositional variations younger than 180 k.y. reflect the growth of a mafic intrusive complex beneath the southern part of the vent axis, which led to the development of evolved magmas. Lastly, these new age data allow for a reassessment of the volcanic recurrence interval at northern Harrat Rahat. Based on available data, volcanism in northern Harrat Rahat over the last 180 k.y. is poorly described using a Poisson distribution with a single recurrence interval. Instead, data for northern Harrat Rahat are better described using a mixed exponential distribution that is applicable for volcanic systems characterized by two different eruptive states, where one state with a longer recurrence interval corresponding to periods of low eruption frequency and one state with a shorter recurrence interval corresponding to periods of high eruption frequency. The preferred model for northern Harrat Rahat over the last 180 k.y. uses a long recurrence interval of 4.0 k.y. and a short recurrence interval of 0.22 k.y.

Professional Paper↗

USGS2023

Lava flow emplacement in Harrat Rahat with implications for eruptions in mafic volcanic fields

Mafic volcanic fields are widespread, but few have erupted in historical times, providing limited observations of the magnitudes, dynamics, and timescales of lava flow emplacement in these settings. The Harrat Rahat volcanic field in western Saudi Arabia offers a good opportunity to study eruptions in such a setting, with a historical eruption in 1256 C.E. (654 in the year of the Hijra) and numerous well-preserved late Pleistocene lava flows. We combine historical observations and rheological and morphological analyses of the youngest flows with analytical models to reconstruct eruptive histories and lava flow emplacement conditions in Harrat Rahat. Petrologic analysis of samples for emplacement temperatures and crystallinities show cooling trends from vent to toe of ~1,140 to ~1,090 degrees Celsius (°C) at rates of 2 to 7 °C per kilometer, crystallinities increasing from 0.5 to 60 volume percent, and apparent viscosities increasing from 10 2 to 10 9 pascal seconds. High-resolution topographic data facilitate quantitative analysis of morphology and interpolation of preeruptive surfaces to measure flow thicknesses, channels, and levees, and enable calculation of eruptive volumes. Analytical models relating flow morphology to emplacement conditions are applied to estimate effusion rates. Within the suite of studied flows, minimum volume estimates range from 0.07 to 0.42 cubic kilometers dense rock equivalent, with effusion rates on the order of tens to hundreds of cubic meters per second and durations from 1 to 15 weeks. These integrated analyses quantify past lava flow emplacement conditions and dynamics in Harrat Rahat, improving our understanding and observations of fundamental parameters and controls of effusive eruptions in Harrat Rahat and other mafic volcanic fields.

northeastern Harrat Rahat↗

USGS2023

The duration and characteristics of magmatic differentiation from basalt to trachyte within the Matan volcanic center, northern Harrat Rahat, Kingdom of Saudi Arabia

A fundamental goal of igneous petrology is to quantify the durations required to produce evolved magmas following influx of basalt into the crust. However, in many cases, complex field relations and (or) the presence of a long-lived magmatic system make it difficult to assess how basaltic inputs relate to more evolved magmas, therefore precluding calculation of meaningful timescales. Here we present field relations, geochemistry, 40 Ar/ 39 Ar ages, and 36 Cl surface-exposure ages for volcanic rocks from the Matan volcanic center, located in the northern part of the Harrat Rahat volcanic field, in the Kingdom of Saudi Arabia. These data document a systematic and repeated temporal progression from alkali basalt to trachyte for the youngest eruptive products. From ~155–17 thousand years ago, the following eruptive sequence occurred four times: (1) alkali basalt, (2) hawaiite, mugearite, or benmoreite, and (3) trachyte. We interpret each eruptive sequence to result from injection of basalt into the crust, and its subsequent differentiation and eruption of progressively evolved magmas. We use the interval time between successive eruptions within a given sequence to calculate the timespans required to produce trachyte from alkali basalt. Differentiation from alkali basalt to intermediate compositions (hawaiite, mugearite, and benmoreite) took ≤3 thousand years (k.y.) on average. Differentiation from intermediate compositions to trachyte took a maximum of 6.7±3.6 to 22.9±1.7 k.y. Thus, the total duration of differentiation was ~10–25 k.y. Timescales presented here are independent of the processes evoked to drive differentiation because they are based solely on the ages and compositions of eruptive products from a system characterized by a simple, repeated differentiation sequence.

northern Harrat Rahat↗

USGS2023

Explosive trachyte eruptions from the Al Efairia volcanic center in northern Harrat Rahat, Kingdom of Saudi Arabia

Harrat Rahat is an alkali basalt, continental, intraplate volcanic field located within the central-western part of the Kingdom of Saudi Arabia. The northern quarter of Harrat Rahat contains evolved volcanic products that achieve trachyte compositions (>60 weight percent SiO 2 ). Within the Al Efairia volcanic center, pyroclastic-flow and -surge deposits that reflect explosive trachyte volcanism (and minor exposed lava domes that reflect effusive trachyte eruptions) sit at the surface as the youngest expression of volcanic activity within this part of the Harrat Rahat volcanic field. Five trachyte deposits emplaced explosively have been identified within the Al Efairia volcanic center based on geologic mapping, petrography, geochemistry, and paleomagnetism. These units are the trachytes of Um Rgaibah, Gura 5, Gura 4, Al Efairia, and Al Qayf, in descending stratigraphic order. Here, we present 14 40 Ar/ 39 Ar analyses from four of these units, which yield eruption ages of 4.2±5.2 thousand years (ka) for the trachyte of Um Rgaibah, 79.7±1.6 ka for the trachyte of Gura 5, 84.3±1.6 ka for the trachyte of Gura 4, and 88.0±1.8 ka for the trachyte of Al Efairia. The eruption age of the trachyte of Al Qayf has been constrained to between 410.3±3.4 and 418.8±1.9 ka using paleomagnetic correlations and 40 Ar/ 39 Ar ages from overlying and underlying intermediate composition lava flows. Most of these trachytes have distinct geochemical compositions, petrographic characteristics, and directions of remanent magnetization. The exceptions are for the trachytes of Gura 4 and Gura 5, which overlap in their geochemical, petrographic, paleomagnetic, and geochronologic affinities. Based on these similarities, we interpret the trachytes of Gura 4 and Gura 5 to have erupted during a closely spaced (a few decades) time interval from the same magma batch but from craters that are >2 kilometers (km) apart. The eruption of the trachyte of Al Efairia at 88.0±1.8 ka is the result of a different magma batch that erupted a few thousand years prior to the trachytes of Gura 4 and Gura 5. The Al Efairia volcanic center is remarkably different from the Matan volcanic center located ~10 km to the north, which has also erupted young (<150 ka) trachytes. The Matan volcanic center has been shown to produce trachyte compositions only after eruption of basalt followed by intermediate lava flows, whereas only trachyte compositions have erupted within the Al Efairia volcanic center over this same time interval.

Al Efairia volcanic center, Harrat Rahat↗

USGS2023

Paleomagnetism of the Harrat Rahat volcanic field, Kingdom of Saudi Arabia—Geologic unit correlations and geomagnetic cryptochron identifications

Paleomagnetic rock samples were collected from 173 drill sites in the Quaternary alkali basaltic volcanic field of northern Harrat Rahat, Kingdom of Saudi Arabia. Laboratory measurements on these samples established that lava flows and vent complexes—identified and mapped from field characteristics, rock types, and compositions as products of single or temporally close eruptions—typically record single, or very similar, directions of remanent magnetization. Correlations defined through geologic mapping, spatial association, geochemistry, geochronology, and identical mean remanent directions indicate at least 16 brief episodes of temporally clustered eruptions. These episodes had durations of a few centuries or less. Anomalous remanent magnetic directions were found for at least 13 mapped lavas of northern Harrat Rahat, which demonstrate that they were acquired during brief geomagnetic cryptochrons during the Brunhes Normal Polarity Chron. These uncharacteristic directions enhance the opportunity to identify common eruptive episodes, and to better understand and evaluate assessments of eruption ages based on 40 Ar/ 39 Ar geochronology. Combining paleomagnetic and regional archaeomagnetic results for the youngest eruptions allows us to evaluate their historical age assignments and, in one case, refute a previously assigned provisional age.

Harrat Rahat volcanic field↗

USGS2023

Mantle origin and crustal differentiation of basalts and hawaiites of northern Harrat Rahat, Kingdom of Saudi Arabia

Quaternary volcanic rocks of northern Harrat Rahat, Kingdom of Saudi Arabia, are chiefly alkali basalts with subordinate transitional basalts, hawaiites, mugearites, benmoreites, and trachytes. Geochemical and isotopic results indicate that crystallization-differentiation, mixing, and cumulate reassimilation within the magmatic system produced most of its compositional diversity, with only minor involvement of Neoproterozoic crust. With increasing evolution, crystal assemblages that separated from and drove basalt-to-hawaiite differentiation passed from (1) dunitic or troctolitic to (2) olivine gabbroic to (3) titanomagnetite-bearing olivine gabbroic, with typical hawaiites representing about 20 weight percent residual liquids from an estimated primary alkali basaltic parent. Crystallization-differentiation simulations for midcrustal pressures yield the closest compositional match to the basalt-hawaiite suite, and stagnation in the midcrustal area (near 20 kilometers [km] depth) may result from density trapping beneath the lower density upper continental crust. Least differentiated alkali basalts have fractionation-adjusted major-oxide compositions that are consistent with separation from the deeper parts (60–70 km) of the spinel-lherzolite stability field at pressures that are close to the local lithosphere-asthenosphere boundary (60–80 km). Mantle potential temperature estimates are strongly sensitive to modeling approach, but potential temperatures (1,345–1,390 degrees Celsius [°C]) are not discernably greater than for midocean ridge basalts (MORB; 1,350–1,410 °C) if adherence to spinel-lherzolite melting relations is required. Inversion of the trace-element concentrations of the lesser fractionated basalts indicates a depleted mantle source, similar to MORB-source estimates, but one that is enriched in Sr and includes a greater relative proportion of melting in the garnet-peridotite stability field. These geochemical and thermal relations, as well as radiogenic isotopes, point to a dominantly depleted mantle asthenospheric source for Harrat Rahat basalts, admixed with subordinate materials either from the Afar mantle plume or an enriched MORB component in the ambient asthenosphere. The lithosphere-asthenosphere boundary is shallower beneath the belt of major volcanic fields on the Arabia Plate, and restoration of rifting across the Red Sea and Gulf of Aden places the south end of this belt adjacent to the northern part of the Afar region, suggesting a once-continuous structure, possibly an arch, a weakness, or a discontinuity along the base of the lithosphere. Magma generation can be ascribed to focused upwelling and decompression melting, perhaps driven by a magmatic-feedback mechanism whereby basaltic intrusions into the deep lithosphere solidify as eclogites, causing lithospheric foundering and further asthenospheric upwelling and decompression melting in a restricted region.

northern Harrat Rahat↗

USGS2023

Isotopic and geochemical evidence for the source of volcanism at Harrat Rahat, Kingdom of Saudi Arabia

Pleistocene and Holocene basalts, hawaiites, mugearites, benmoreites, and trachytes from the northern part of the Harrat Rahat volcanic field, Kingdom of Saudi Arabia, were analyzed for Sr, Nd, Hf, and Pb isotopic compositions. Evolved trachytes with Mg number <0.1 (Mg# = Mg/[Mg+Fe 2+ ], molar) have relatively radiogenic Sr isotopic compositions indicating that they were influenced by contamination probably in the upper crust. Volcanic rocks with Mg# >0.1, consisting chiefly of alkali basalts but encompassing hawaiites, mugearites, and benmoreites, show a limited range in Hf, Nd, Sr, and Pb isotopic compositions. Although the total Pb isotope variation is only 1 percent, the Pb isotope values correlate with Mg#, where the least radiogenic Pb is in samples with the lowest Mg#. The trend formed in Pb isotope space points toward an unradiogenic Pb composition that is similar to the Pb isotopic composition of lower crust of the Precambrian Arabian-Nubian Shields, as well as to feldspars and galena in the upper crust of the western Arabian Shield. This trend is interpreted as progressive but overall minor (no more than 5 weight percent) assimilation of shield rocks, or their partial melts, during fractional crystallization. Isotopic compositions of the least evolved northern Harrat Rahat magmas are most similar among analyzed Arabian harrats to depleted spreading-ridge basalts of the active Red Sea rift, but isotopic values are displaced toward those of spreading-ridge basalts of the Gulf of Aden that are proximal to the site of the Afar mantle plume. The Pb isotopic compositions very near the Northern Hemisphere Reference Line indicate no discernable lithospheric contribution to yield the parental basalts of northern Harrat Rahat, and their isotopic compositions are consistent with derivation predominantly from depleted Northern Hemisphere asthenosphere with a subordinate (20–30 weight percent) component from the Afar mantle plume. Trace-element variations show that appreciable portions of melting were in the garnet stability field, confirming the sub-lithospheric origin of the magmas, and that melting extents were low, accounting for the alkalic, trace-element-enriched character of the suite. The presence of possible Afar mantle beneath the western part of the Arabian Shield and its absence beneath the Red Sea rift may result from capture and channelized flow along high-relief structures incised into the base of the sub-continental lithosphere, as revealed by geophysical images.

northern Harrat Rahat↗

USGS2023

Spatio-temporal clustering of post-caldera eruptions at Yellowstone caldera: Implications for volcanic hazards and pre-eruptive magma reservoir configuration

At least 28 rhyolitic lava flows, domes, and tuffs erupted within Yellowstone caldera following its formation 631 ka. Understanding the timing of intracaldera eruptions is essential for characterizing natural hazards posed by Yellowstone volcano. We present 40 Ar/ 39 Ar eruption ages for the Mallard Lake Member and Central Plateau Member of the Plateau Rhyolite, which comprise products from 22 eruptions spanning ~160 to ~70 ka. These eruptions occurred along two linear vent zones within Yellowstone caldera that appear to be extensions of major extracaldera normal faults. Eruptions took place in five brief clusters with group mean ages of 160.3±1.0 ka (2σ), 149.8±4.0 ka, 111.4±0.9 ka, 104.1±0.8 ka, and 70.8±0.7 ka. All mapped products within each informal group have indistinguishable 40 Ar/ 39 Ar eruption ages implying brief eruption durations. Eruption products of each group are spatially clustered, with only one of the two linear vent zones being active during an eruption group. Using the oldest group (160 ka), we apply paleomagnetic and geochemical analyses to investigate the duration of these eruption groups and their pre-eruptive magma reservoir configuration. Paleomagnetic analyses suggest that the nine rhyolite eruptions ~160 ka spanned 400 years or less. Glass and sanidine compositions of each rhyolite are nearly identical, which is consistent with pre-eruptive storage of a large (~130 km 3 ), interconnected melt body within Yellowstone volcano’s crystal-rich magma reservoir. These characteristics differ from intracaldera eruptions that vented near the caldera’s inner ring fracture system, which tend to be compositionally disparate with less temporal clustering. These insights inform hazard mitigation scenarios.

Idaho, Wyoming↗

USGS2023

Evolution of Miocene normal and dextral faulting in the lower Colorado River region near Blythe, California, USA

The evolution of strain in nascent continental plate boundaries commonly involves distributed deformation and transitions between different styles of deformation as the plate boundary matures. Distributed NW-striking faults, many with km-scale right-lateral separation, are prevalent near Blythe, California, and have been variably interpreted to have accommodated either Middle Miocene NE-SW extension as normal faults or Late Miocene to Pliocene dextral shear as strike-slip faults. However, with poor timing and kinematic constraints, it is unclear how these faults relate to known domains of Neogene deformation and the evolution of the Pacific–North America plate boundary. We present kinematic data (n = 642 fault planes, n = 512 slickenlines) that demonstrate that these faults dominantly dip steeply northeast; ~96% of measured faults record normal, dextral, or oblique dextral-normal kinematics that likely reflect a gradational transition between normal and dextral oblique kinematic regimes. We constrain fault timing with 11.7 Ma and 7.0 Ma 40 Ar/ 39 Ar dates of rocks cut by faults, and laser ablation–inductively coupled plasma–mass spectrometry U-Pb dating of calcite mineralized during oblique dextral faulting that demonstrates fault slip at ca. 10–7 Ma and perhaps as late as ca. 4 Ma. This Late Miocene dextral oblique faulting is best compatible with a documented regional transition from Early to Middle Miocene NE-directed extension during detachment fault slip to subsequent NW-directed dextral shear. We estimate 11–38 km of cumulative dextral slip occurred across a 50-km-wide zone from the Palen to Riverside mountains, including up to 20 km of newly documented dextral shear that may partly alleviate the regional discrepancy of cumulative dextral shear along this part of the Late Miocene Pacific–North America plate boundary.

California↗

USGS2022

Physical and geochemical architecture and age of the Pliocene Bicycle Lake basalt, southeastern Fort Irwin, California

The informally named Bicycle Lake basalt forms a volcanic field in southeastern Fort Irwin, California, disrupted by three east-striking faults and linked cross-faults of the Eastern California Shear Zone, and its distribution provides a framework for assessing volcanic field development, groundwater resources, and fault offsets. Previous geochronologic studies yielded ages ranging from ~2.9-5.6 Ma, and a new cycle of dating reported herein yields a more robust age of ~4.55 Ma. No vents or dikes have been identified for the basalt. No interstratified sedimentary rocks have been identified, not even aeolianites, implying that the field was monogenic and short-lived. Thickness of the field varies from 1 m at distal edges to as much as 20-35 m (projected in cross-section) in the central part of the field. The field consists of lava flows <8 m thick, and the relative abundance of flows <2 m thick and 3-8 m thick varies across the field, as does the interstratification of <3 m thick cosets of <1 m thick flows. Most flows are basaltic andesite with minor basalt. Within stratigraphic sequences, variations in major oxides and trace element geochemistry indicate interstratification of basalt with basaltic andesite, but also different compositions of interstratified basaltic andesite. The local abundance of basalt or different compositions of basaltic andesite varies. The distribution of compositions suggests minor changes in the composition of magma erupted and minor sector partitioning of flow paths during growth of the field. The basalt has been deformed as a result of transpressional stress transferred from the Coyote Lake to Bicycle Lake faults, which formed NW- striking reverse faults and faulted kink folds along with broad arching of the basalt. Restoring separation along the Coyote Canyon and Bicycle Lake faults indicates some flows of the basalt field traveled at least 10 km from their inferred source.

California↗