Geology topics

Thomas W. Sisson

Publications and source records attributed to Thomas W. Sisson.

At least 19 recordsLinked to original sources

USGS2026

Cenozoic distributed volcanism of the Arabia Plate—A review

Cenozoic volcanic rocks of the Arabia Plate cover about 140,000 square kilometers across a distance of about 3,000 kilometers from southern Yemen to southeastern Turkey. The majority of volcanic products are alkali basalts that erupted in restricted areas, commonly over periods of a million or more years, building mafic lava fields, each known in Arabic as a “harrat.” Harrat volcanism commenced following the Oligocene flood-lava effusions that blanketed the (now) Ethiopian highlands, southern Sudan, and western Yemen, and overlapped the latest Oligocene to early Miocene initial riftings of the Red Sea and Gulf of Aden, but the majority of harrat volcanism has been since approximately 13–10 million years ago. Persistent harrat magmatism in restricted locations led to the development of intermediate and evolved magmas (hawaiites, mugearites, benmoreites, trachytes, comendites, and phonolites) mainly through intracrustal crystallization-differentiation. Most of these intermediate and evolved magmas erupt at sites of the greatest aggregate volcanic relief, reflecting sites of the greatest overall magmatic fluxes. Production of fractionated magmas at these sites negates “monogenetic” as an appropriate descriptor of harrat volcanism. This chapter summarizes the geologic, eruptive, and tectonic history and aspects of the petrogenesis of the Cenozoic Arabian alkalic province. Particular emphasis is placed on results of a joint study of Ḩarrat Rahat adjacent to the city of Al Madīnah al Munawwarah, Kingdom of Saudi Arabia, published as U.S. Geological Survey Professional Paper 1862 and Saudi Geological Survey Special Report SGS–SP–2021–1. A goal of this chapter is to provide an introduction to those unfamiliar with this vast, enigmatic, and fascinating region of distributed continental volcanism.

Arabia Plate↗

USGS2026

Shallow differentiation of primitive arc magmas at the Jurassic Emigrant Gap mafic complex, Sierra Nevada, California

The Emigrant Gap composite pluton exposes ultramafic to silicic intrusive rocks that preserve the chemical evolution of primitive mafic arc magmas and their open-system interactions in the upper crust during mid-Jurassic growth of the Sierra Nevada batholith (California). We present field and petrographic observations and mineral and whole-rock chemistry of the ~35-km 2 ultramafic to dioritic Emigrant Gap mafic complex and an adjacent penecontemporaneous ~90-km 2 granodiorite that together make up the composite pluton. In the Emigrant Gap mafic complex, four roughly central masses of dunite, wehrlite, and olivine clinopyroxenite are surrounded by weakly layered gabbronorite and non-layered diorite. The ultramafic rocks are cumulates formed from near-liquidus minerals of primitive arc magmas that accumulated in steep feeder zones with substantial modification by melt–mush reaction as primitive liquids repeatedly transited the mush-filled conduits. The dominant gabbronoritic rocks are the variably accumulative products of more advanced crystallization–differentiation of arc tholeiitic basalts and basaltic andesites. The adjacent granodiorite intrusion originated separately and preserves field and geochemical evidence for assimilation of metasedimentary rocks. Open-system hybridization between the gabbronoritic mushes and the granodioritic magma produced an intervening body of two-pyroxene diorite. We infer that the ultramafic rocks and gabbronorite of the Emigrant Gap mafic complex crystallized from near-primitive arc basaltic to basaltic andesitic magmas at ~0.15–0.3 GPa, with estimated f O 2 of ≥FMQ +1 and dissolved H 2 O concentrations of only ~0.5–2 wt %. Notably, the Emigrant Gap composite pluton is distinct from other Mesozoic plutons in the Sierra Nevada batholith because of (1) its abundance of mafic and ultramafic rocks that crystallized from relatively primitive mafic melts and (2) the low inferred H 2 O concentrations of its parental magmas, indicated by a near absence of igneous amphibole and by the intermediate rather than calcic compositions of plagioclase. A Jurassic regional extension event probably accounts for the formation of relatively dry primitive arc magmas, as well as for their ascending to the upper crust.

California↗

USGS2025

Igneous differentiation in arcs

This chapter assesses crystallization-differentiation of convergent-margin magmas both from the standpoints of natural rock suites and from laboratory phase-equilibrium experiments. Crystallization-differentiation of dry subalkaline and alkaline basalts are summarized at the beginning for reference and context. Experiments crystallizing arc basalts through andesites at middle to upper crustal pressures (100–400 MPa), saturated either with H 2 O or mixed H 2 O-CO 2 , produce liquid lines of descent similar to arc igneous suites, with higher f O 2 s giving results closer to calc-alkaline magma series. Early enrichment of melts in SiO 2 with little or no enrichment in FeO* (total Fe as FeO) comes about by the combined effects of H 2 O suppressing crystallization of plagioclase relative to olivine, shifting plagioclase to more calcic – therefore lower SiO 2 – compositions, and suppressing silicate minerals to greater extents than Fe-rich spinel. Also contributing to SiO 2 enrichment is crystallization near the liquidus of higher proportions of olivine relative to calcic clinopyroxene, compared with MORB, but it is unclear if this is brought about by H 2 O or elevated pressure. Exhumed arc terranes, however, generally fail to expose evidence for primitive liquids having reached and differentiated at middle to upper crustal depths, at least in any abundance. Instead, most differentiation in arcs appears to take place in the lower crust or near the crust-mantle transition, and although hydrous experiments at such pressures produce liquids broadly similar to arc suites in their SiO 2 enrichment, the experimental liquids differ in becoming peraluminous at early stages of differentiation while melts have basaltic, basaltic andesitic, or andesitic SiO 2 concentrations. Crystallization of abundant calcic clinopyroxene close to the liquidus causes the high-pressure experimental liquids to become peraluminous with only small increases in SiO 2 concentration. In contrast, natural arc volcanic and plutonic suites become gradually less metaluminous with increasing SiO 2 , crossing to peraluminous at dacitic to rhyolitic SiO 2 concentrations. Hybridism, encompassing magma-mixing, assimilation of old crust, and incorporation of near-solidus melts from only shortly older deep, antecedent arc intrusions may accompany natural crystallization-differentiation, preventing melts from attaining peraluminous compositions at low-SiO 2 concentrations. If so, then nearly all arc suites are to substantial extents composites, rather than preserving simple, progressive liquid lines of descent. Some comments and assessments are also made on the participation of garnet in arc differentiation and on the production of arc batholiths.

Book chapter↗

USGS2024

Shallow storage of the explosive Earthquake Flat Pyroclastics magma body, Okataina Volcanic Center, Taupo Volcanic Zone, New Zealand: Evidence from phase-equilibria experiments

Rhyolitic tuffs range widely in their crystal contents from nearly aphyric to crystal-rich, and their crystal cargoes inform concepts of upper crustal magma reservoirs. The Earthquake Flat pyroclastics (Okataina Volcanic Center, Taupo Volcanic Zone, New Zealand) are 10 km 3 of rhyolitic tuffs with abundant (~ 40 vol.%) plagioclase and quartz, minor biotite, hornblende, and orthopyroxene, and accessory Fe-Ti oxides, apatite, and zircon, set in high-silica rhyolitic glass. Major minerals form large, euhedral phenocrysts and abundant glomerocrysts with few disequilibrium textures excepting some faintly resorbed quartz. Plagioclase phenocrysts have thick rims of nearly constant composition near An 30 , and hornblende is weakly zoned or unzoned. The abundant and texturally complex mineral assemblage contrasts with the nearby (~ 25 km), nearly synchronous, but more voluminous and crystal-moderate rhyolite tuffs from Rotoiti caldera. New H 2 O-saturated phase-equilibria results on the erupted Earthquake Flat melt (glass) determine its co-saturation with the partial phenocryst assemblage of plagioclase, quartz, biotite, and Fe-Ti oxides at: 140 MPa, 755 ºC. These closely approximate the conditions of the pre-eruptive magma body assuming it was saturated with nearly pure H 2 O and at an f O 2 of ~ Ni–NiO. Absence of hornblende and orthopyroxene from the synthesized assemblages may result from those minerals being in a peritectic reaction relation with melt to produce biotite, so they would not grow from the liquid used as starting material. Experimental results on Rotoiti rhyolite (Nicholls et al. 1992) show that the two bodies resided at similar pressures, temperatures, and f O 2 s. Lower crystal abundance of the Rotoiti tuffs may result from slight compositional differences. We interpret that the Earthquake Flat pyroclastics were sourced from the crystal-rich periphery of a mushy reservoir system with the Rotoiti occupying a more melt-rich central location. Uncertain is whether this was a single intrusion zoned continuously in crystallinity, or discrete adjacent intrusions, but our results illustrate and quantify complexities of magma storage across relatively short distances.

North Island, Taupo Volcanic Zone↗

USGS2023

The Saudi Geological Survey-U.S. Geological Survey northern Harrat Rahat project—Styles, rates, causes, and hazards of volcanism near Al Madīnah al Munawwarah, Kingdom of Saudi Arabia

Active volcanic systems pose serious hazards to people and property including inundation and incineration by lava, blanketing by tephra (volcanic ash), exposure to noxious volcanic gases, and damage from shallow earthquakes triggered by ascending molten material (magma). To improve understanding of volcanism and associated seismicity on the western Arabia Plate, the Saudi Geological Survey and the U.S. Geological Survey conducted a multi-year investigation of the northern Harrat Rahat volcanic field adjacent to the city of Al Madīnah al Munawwarah, Kingdom of Saudi Arabia. Project components included creation of a high-resolution digital topographic base; interpretation of eruptive history supported by detailed geologic mapping, paleomagnetism, and abundant high-precision geochronology of volcanic deposits; assessments of eruptive styles and volcanic hazards by physical volcanology; investigation of the origins of magmas in the mantle and of their differentiation in the crust revealed by chemical and isotopic petrology; gravity and magnetotelluric surveys to reveal crustal structures and to search for magma reservoirs; and regional and local seismic tomography and analyses of seismic hazards. Project results are presented in this Professional Paper as chapters written for technical scientific audiences. This initial chapter introduces the project and briefly summarizes results in plain language for readers who have more general backgrounds.

Al Madīnah al Munawwarah↗

USGS2023

Cenozoic tectonics of the western Arabia Plate related to harrat magmatism near Al Madīnah, Kingdom of Saudi Arabia

Sprawling volcanic fields, or harrats, in western Saudi Arabia have been emplaced during the past 30 million years following effusions of flood basalts in Ethiopia and Yemen. Although broadly associated with volcanism in three rift valleys (Red Sea, Gulf of Aden, and East African Rift Zone) radiating from the Afar depression, harrat abundance on the Arabian Peninsula indicates that volcanism is distinct from processes that opened the Red Sea and drive its magmatism. Harrats primarily lie unconformably on the Arabian Shield and locally on Paleozoic platform deposits, but several are conformable upon or interbedded with Paleogene supracrustal strata. Harrats erupted from linear chains of vents; some are oriented parallel to the Red Sea (N. 30° W.), but most trend northerly or north-northeasterly. Harrats consist mainly of weakly alkalic basalts and lesser amounts of their differentiation products, but some contain strongly alkalic rocks such as basanites and phonolites. Early harrats predate initial opening of the Red Sea 24–20 million years ago (Ma) and uplift of its flank 20–11 Ma, but younger harrats postdate the uplift and are larger and more abundant. Seismic studies identify substantially thinned continental lithosphere beneath the most voluminous harrats, but little crustal thinning is apparent at the surface, an indication that the lithosphere was eroded from below. The combination of plume impingement, rift valley formation, and erosion of the lithosphere beneath the Arabian Shield since 30 Ma indicates that the Afar plume thermally weakened the region, enabling the Arabia Plate to separate from the Africa Plate along the Red Sea and Gulf of Aden. The plume also thinned interior parts of the newly formed Arabia Plate, not just its margins, promoting mantle melting and volcanism.

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

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

Garnet stability in arc basalt, andesite, and dacite—An experimental study

Garnet’s stability in arc magmas and its influences on their differentiation were explored experimentally in a typical basalt, andesite, and dacite at conditions of 0.9–1.67 GPa, 800–1300 °C, with 2–9 wt.% added H 2 O, and with oxygen fugacity buffered near Re + O 2 = ReO 2 (~ Ni-NiO + 1.7 log 10 bars). Garnet did not grow at 0.9 GPa in any of the compositions, even with garnet seeds added to facilitate nucleation. At 1.0–1.2 GPa, garnet grew as thin rims (< 5 µm) on introduced garnet seeds coexisting with dacitic to rhyodacitic liquids at temperatures ≤ 1000 °C. At 1.3 GPa, garnet grew readily with no seeds from 900 to 1100 °C coexisting with liquids ranging from peraluminous basaltic andesite to rhyodacite, and at 1.46 GPa, garnet was stable as hot as 1150 °C in metaluminous basaltic liquid. Garnet grew as a liquidus phase only in the dacite, a composition similar to the average upper continental crust. Inverse experiments on the dacite determined a liquidus multiple-saturation point with garnet, plagioclase, orthopyroxene, calcic clinopyroxene, and amphibole at 975 °C, 1.46 GPa, with 7 wt.% dissolved H 2 O. Such dacitic and more evolved melts can be products of peritectic reactions that with decreasing temperature consume garnet, calcic clinopyroxene, and melt components, producing amphibole and less abundant but more evolved melts. For this reason, experiments on product melts need not produce reactant minerals, accounting for some disparities in published experimental results on the apparent stability of garnet in intermediate-to-evolved arc magmas. Results on more mafic compositions are more reliable guides and show that liquids of arc dacitic composition, and more evolved compositions, would coexist stably with garnet only in the deepest portions of continental-margin arc crust with average thickness and density (~ 43 km, ~ 1.2 GPa) or in the underlying shallow mantle. Metaluminous arc basaltic, basaltic andesitic, and many andesitic liquids would not coexist stably with garnet at pressures ranging from the crust to at least the midpoint of the mantle wedge, but results in the literature allow that some andesitic liquids with higher Fe/Mg than common in arcs may also saturate with garnet in the deeper portions of average-thickness continental arc crust.

Contributions to Mineralogy and Petrology↗

USGS2022

Geologic field-trip guide to volcanism and its interaction with snow and ice at Mount Rainier, Washington

Mount Rainier is the Pacific Northwest’s iconic volcano. At 4,393 meters and situated in the south-central Cascade Range of Washington State, it towers over cities of the Puget Lowland. As the highest summit in the Cascade Range, Mount Rainier hosts 26 glaciers and numerous permanent snow fields covering 87 square kilometers and having a snow and ice volume of about 3.8 cubic kilometers. It remains by far the most heavily glacier-clad mountain in the conterminous United States despite having lost about 14 percent of its ice volume between 1970 and 2008. Five major rivers head at Mount Rainier—the White, Carbon, Puyallup, Nisqually, and Cowlitz Rivers. Because Mount Rainier is situated west of the Cascade Range crest, all of these rivers eventually turn and drain westward. The Puget Lowland, situated west to northwest of Mount Rainier, is the Pacific Northwest’s most densely populated area, including Seattle, Tacoma, and Olympia. The Puget Lowland is now home to a population of more than 4.5 million and a vibrant economy. Mount Rainier is one of the most hazardous volcanoes in the United States, not so much because of its explosivity, but rather because of its frequent eruptions, its propensity to produce voluminous far-traveled lahars, and its proximity to large population centers of the Puget Lowland. Steep-sided, glacially carved valleys serve as lahar conduits, and even mild eruptions commonly produced large lahars that traveled into areas now populated by hundreds of thousands of people. This guide describes a five-day field trip to view the geology of Mount Rainier as it relates to volcanism and its interaction with snow and ice. Day 1 will focus on lahars in the White River valley. We will drive to Enumclaw, Washington, to begin the day then work our way back upvalley toward Mount Rainier. Day 2 concentrates on geology of the Sunrise-Glacier Basin area within Mount Rainier National Park. As part of day 2 activities, we will hike about 10 miles from Sunrise to the top of Burroughs Mountain, down into Glacier Basin, and be picked up at White River Campground. On day 3 we will pack up and move to Paradise, stopping to examine geology along Stevens Canyon Road. We will hike from Paradise along the Golden Gate Trail and eventually eastward to the former Paradise Ice Caves area (the ice caves have melted out). Day 4 involves hiking from Comet Falls trailhead to Mildred Point and return (~7 miles; 11 km), examining geology along the way. During the first half of day 5, we will visit sites on the south side of Mount Rainier to study lahar deposits, then return to the tour origin.

Washington↗

USGS2021

Constraints on deep, CO2-rich degassing at arc volcanoes from solubility experiments on hydrous basaltic andesite of Pavlof Volcano, Alaska Peninsula, at 300 to 1200 MPa

The solubility of CO 2 in hydrous basaltic andesite was examined in f O2 -controlled experiments at a temperature of 1125 °C and pressures between 310–1200 MPa. Concentrations of dissolved H 2 O and CO 2 in experimental glasses were determined by ion microprobe calibrated on a subset of run glasses analyzed by high-temperature vacuum manometry. Assuming that the solubility of H 2 O in mafic melt is relatively well known, estimates of 𝑋H2Ofluid and 𝑃H2Ofluid in the saturating fluid were modeled, and by difference, values for 𝑋CO2fluid and 𝑃CO2fluid were obtained ( X CO2 ~0.5–0.9); f CO2 could be then calculated from the fluid composition, temperature, and pressure. Dissolved H 2 O over a range of 2.3–5.5 wt% had no unequivocal influence on the dissolution of CO 2 at the pressures and fluid compositions examined. For these H 2 O concentrations, dissolved CO 2 increases with f CO2 following an empirical power-law relation: dissolved CO 2 (ppmw) = 14.9−3.5+4.5 [ f CO2 (MPa)] 0.7±0.03 . The highest-pressure results plot farthest from this equation but are within its 1 standard-error uncertainty envelope. We compare our experimental data with three recent CO 2 -H 2 O solubility models: Papale et al. (2006) ; Iacono-Marziano et al. (2012) ; and Ghiorso and Gualda (2015) . The Papale et al. (2006) and Iacono-Marizano et al. (2012) models give similar results, both over-predicting the solubility of CO 2 in a melt of the Pavlof basaltic andesite composition across the f CO2 range, whereas the Ghiorso and Gualda (2015) model under-predicts CO 2 solubility. All three solubility models would indicate a strong enhancement of CO 2 solubility with increasing dissolved H 2 O not apparent in our results. We also examine our results in the context of previous high-pressure CO 2 solubility experiments on basaltic melts. Dissolved CO 2 correlates positively with mole fraction (Na+K+Ca)/Al across a compositional spectrum of trachybasalt-alkali basalt-tholeiite-icelandite-basaltic andesite. Shortcomings of current solubility models for a widespread arc magma type indicate that our understanding of degassing in the deep crust and uppermost mantle remains semi-quantitative. Experimental studies systematically varying concentrations of melt components (Mg, Ca, Na, K, Al, Si) may be necessary to identify solubility reactions, quantify their equilibrium constants, and thereby build an accurate and generally applicable solubility model.

Alaska↗

USGS2020

Felsic melt and gas mobilisation during magma solidification: An experimental study at 1.1 kbar

Melt and gas transfer processes are essential to the formation and growth of the Earth’s crust and for sustaining volcanic activity. These processes also play a major role in magma fractionation at shallow depths (<10 km) where magmas stall rheologically and solidify. In this scenario, the conditions of melt and gas mobilization during progressive cooling of crystal mushes down to their solidus remain poorly understood. We present experimental data (at 1.1 kbar) showing how a combination of temperature and crystal content control the ability of melt and gas to escape from cooling and solidifying hydrous silicic magmas with initial crystal volume fractions (Φ) of 0.6, 0.7, and 0.8, and for temperature snapshots of 850, 800, and 750°. Microstructural observations and chemical data show that the amount of extracted melt increases by 70% from 850 to 750° and by 40% from Φ = 0.6 to 0.8 at 750°, due to the formation of interconnected crystal frameworks, gas expansion in constricted pore space, and filter pressing during cooling. As a result, our experiments suggest that melt and gas extraction from cooling mushes increases in proximity to their solidus and can operate efficiently at 0.6 < Φ <0.93. These observations shed light on maximum estimates of the segregation of gas-rich, crystal-poor magmas (0.02 m/year at 850° to 9 m/year at 750°) to form felsic dykes or eruptible systems feeding volcanoes.

Frontiers in Earth Science↗

USGS2020

Geologic map of the northern Harrat Rahat volcanic field, Kingdom of Saudi Arabia

Harrat Rahat, in the west-central part of the Kingdom of Saudi Arabia, is the largest of 15 Cenozoic harrats (Arabic for “volcanic field”) distributed on the Arabian plate. It extends more than 300 km north-south and 50 to 75 km east-west, and it covers an area of approximately 20,000 km 2 , has a volume of approximately 2,000 km 3 , and encompasses more than 900 observable vents. Volcanism commenced around 10 Ma and has continued into historic time, the most recent eruption occurring in 1256 C.E. Volcanic products are dominated by alkali basalt and hawaiite lava flows, with subordinate mugearite lava flows, as well as benmoreite and trachyte lava flows, domes, and pyroclastic flows. This geologic map distinguishes 239 eruptive units that cover an area of 3,340 km 2 in northern Harrat Rahat and the adjacent city of Al-Madinah. Results are presented as a geologic map of the study area at 1:75,000 scale and of smaller regions of particular interest at 1:25,000 scale, along with interpretive text. Most units are basaltic lava flows that erupted from the broadly north-northwest-trending main vent axis that constructed the topographic crest of the volcanic field. This 300- to 400-m-high vent axis, which has a width of 6 to 10 km, lies in the eastern one-third of northern Harrat Rahat. Basalt and hawaiite lava flows can extend as far as 27 km from their vents, but most are 10 to 15 km long. Evolved products such as mugearites, benmoreites, and trachytes are less extensive; the trachytic pyroclastic flows extend as far as 9 km from their source vents, although most only reach 4 to 6 km. Vents of the evolved products are restricted to the main vent axis or its flanks. No volcanic rocks older than 1.2 Ma are exposed in the map area, and about 90 percent of the exposed volcanic rocks erupted during the past 570 thousand years. As depicted on the geologic maps, eruption ages and field relations define 12 eruptive stages for northern Harrat Rahat for the past 1.2 million years. Other important geochronological findings include (1) several late Pleistocene lava flows near Al-Madinah, which previously were interpreted as Holocene from archeological evidence; (2) the eruption age of a cluster of cinder cones and small lava flows in the western outskirts of Al-Madinah (previously ascribed to an eruption in 641 C.E.) is actually 13.3±1.9 ka, close to the Pleistocene-Holocene boundary; and (3) only two Holocene eruptions have been identified in the map area, those of the historically described basalt of Al Labah in 1256 C.E. and the dome and pyroclastic flows of the trachyte of Um Rgaibah at 4.2±5.2 ka.

Scientific Investigations Map↗

USGS2019

Age of the dacite of Sunset Amphitheater, a voluminous Pleistocene tephra from Mount Rainier (USA), and implications for Cascade glacial stratigraphy

The dacite of Sunset Amphitheater, Mount Rainier (USA), illustrates the difficulties in establishing accurate ages of Pleistocene tephra eruptions. Nearly uniform whole-rock, glass, and mineral compositions, texture, and phenocryst assemblage establish that certain conspicuous dissected pumice exposures scattered from Mount Rainier to southern Puget Sound are products of the same Pleistocene Plinian eruption . Deposit thicknesses and pumice sizes support an eruption on the order of low Volcanic Explosivity Index (VEI) 5, atypically explosive for dominantly lava-producing Mount Rainier. Statistically permissible 40 Ar/ 39 Ar plateau ages of plagioclase phenocryst separates are 138 ± 20 ka and 101 ± 11 ka (2σ). A previously published result of 206 ± 11 ka is herein shown to result from a sample selection error. Zircon from the pumice yields a U-Th crystallization age of 147 ± 8 ka if the isochron is required to pass through the tephra U-Th isotopic composition . In contrast, pooled (U-Th)/He measurements on the zircon yield an age of 85 ± 6 ka (2σ), which accords with well-behaved 40 Ar/ 39 Ar ages of stratigraphically associated lavas high on Mount Rainier, and is the best estimate of the pumice's true eruption age. Inclusions of undegassed melt (glass) in the plagioclase separates are proposed as biasing apparent 40 Ar/ 39 Ar plateau ages to old values through coupling of undegassed magmatic excess Ar with radiogenic Ar that accumulated post-eruptively from relatively K-rich glass. U-Th ages record zircon growth prior to eruption, consistent with a possible complex history of advanced solidification followed by remobilization . The ca. 85 ka eruption age confirms that bracketing glacial tills on the flanks of Mount Rainier were products of the Penultimate Glaciation (MIS 6) and Last Glacial Maximum (MIS 2). This eruption age also provides an important time marker for glacial and other sedimentary deposits in southern Puget Sound lowland that, excepting the Vashon Drift (MIS 2), generally lack reliable age determinations.

Washington↗