James D. Muirhead

Lateral magma flow in mafic sill complexes

The structure of upper crustal magma plumbing systems controls the distribution of volcanism and influences tectonic processes. However, delineating the structure and volume of plumbing systems is difficult because (1) active intrusion networks cannot be directly accessed; (2) field outcrops are commonly limited; and (3) geophysical data imaging the subsurface are restricted in areal extent and resolution. This has led to models involving the vertical transfer of magma via dikes, extending from a melt source to overlying reservoirs and eruption sites, being favored in the volcanic literature. However, while there is a wealth of evidence to support the occurrence of dike-dominated systems, we synthesize field- and seismic reflection–based observations and highlight that extensive lateral magma transport (as much as 4100 km) may occur within mafic sill complexes. Most of these mafic sill complexes occur in sedimentary basins (e.g., the Karoo Basin, South Africa), although some intrude crystalline continental crust (e.g., the Yilgarn craton, Australia), and consist of interconnected sills and inclined sheets. Sill complex emplacement is largely controlled by host-rock lithology and structure and the state of stress. We argue that plumbing systems need not be dominated by dikes and that magma can be transported within widespread sill complexes, promoting the development of volcanoes that do not overlie the melt source. However, the extent to which active volcanic systems and rifted margins are underlain by sill complexes remains poorly constrained, despite important implications for elucidating magmatic processes, melt volumes, and melt sources.

Interconnected sills and inclined sheet intrusions control shallow magma transport in the Ferrar large igneous province, Antarctica

Field observations and structural data from intrusive complexes at Allan Hills and Terra Cotta Mountain, South Victoria Land, Antarctica, demonstrate that interconnected sills and inclined sheets transported magma through the shallow subsurface. These sills and sheets represent the upper-crustal (top 4 km) plumbing system of the 183 Ma Ferrar large igneous province. The sheets are short in length (

Varying styles of magmatic strain accommodation across the East African Rift

Observations of active dike intrusions provide present day snapshots of the magmatic contribution to continental rifting. However, unravelling the contributions of upper crustal dikes over the timescale of continental rift evolution is a significant challenge. To address this issue, we analyzed the morphologies and alignments of >1500 volcanic cones to infer the distribution and trends of upper crustal dikes in various rift basins across the East African Rift (EAR). Cone lineament data reveal along-axis variations in the distribution and geometries of dike intrusions as a result of changing tectonomagmatic conditions. In younger ( 10 Ma) in Ethiopia and the Kenya Rift, rift-parallel dikes accommodate upper crustal extension along the full length of the basin.

Evolution of upper crustal faulting assisted by magmatic volatile release during early-stage continental rift development in the East African Rift

During the development of continental rifts, strain accommodation shifts from border faults to intra-rift faults. This transition represents a critical process in the evolution of rift basins in the East African Rift, resulting in the focusing of strain and, ultimately, continental breakup. An analysis of fault and fluid systems in the younger than 7 Ma Natron and Magadi basins (Kenya-Tanzania border) reveals the transition as a complex interaction between plate flexure, magma emplacement, and magmatic volatile release. Rift basin development was investigated by analyzing fault systems, lava chronology, and geochemistry of spring systems. Results show that extensional strain in the 3 Ma Natron basin is primarily accommodated along the border fault, whereas results from the 7 Ma Magadi basin reveal a transition to intra-rift fault–dominated strain accommodation. The focusing of strain into a system of intra-rift faults in Magadi also occurred without oblique-style rifting, as is observed in Ethiopia, and border fault hanging-wall flexure can account for only a minor portion of faulting along the central rift axis (∼12% or less). Instead, areas of high upper crustal strain coincide with the presence of hydrothermal springs that exhibit carbon isotopes and N2-He-Ar abundances indicating mixing between mantle-derived (magmatic) fluids and air saturated water. By comparing the distribution of fault-related strain and zones of magmatic fluid release in the 3 Ma Natron and 7 Ma Magadi basins, we present a conceptual model for the evolution of early-stage rifting. In the first 3 m.y., border faults accommodate the majority of regional extension (1.24–1.78 mm yr–1 in Natron at a slip rate ranging 1.93–3.56 mm yr–1), with a significant portion of intra-rift faulting (38%–96%) driven by flexure of the border fault hanging wall. Fluids released from magma bodies ascend along the border fault and then outward into nearby faults forming in the flexing hanging wall. By 7 m.y., there is a reduction in the amount of extension accommodated along the border fault (0.40–0.66 mm yr–1 in Magadi at a slip rate ranging from 0.62 to 1.32 mm yr–1), and regional extension is primarily accommodated in the intra-rift fault population (1.34–1.60 mm yr–1), with an accompanying transition of magmatic volatile release into the rift center. The focusing of magma toward the rift center and concomitant release of magmatic fluids into the flexing hanging wall provides a previously unrecognized mechanism that may help to weaken crust and assist the transition to intra-rift dominated strain accommodation. We conclude that the flow of magmatic fluids within fault systems plays an important role in weakening lithosphere and focusing upper crustal strain in early-stage continental rift basins prior to the establishment of magmatic segments.

Initial pulse of Siberian Traps sills as the trigger of the end-Permian mass extinction

Mass extinction events are short-lived and characterized by catastrophic biosphere collapse and subsequent reorganization. Their abrupt nature necessitates a similarly short-lived trigger, and large igneous province magmatism is often implicated. However, large igneous provinces are long-lived compared to mass extinctions. Therefore, if large igneous provinces are an effective trigger, a subinterval of magmatism must be responsible for driving deleterious environmental effects. The onset of Earth’s most severe extinction, the end-Permian, coincided with an abrupt change in the emplacement style of the contemporaneous Siberian Traps large igneous province, from dominantly flood lavas to sill intrusions. Here we identify the initial emplacement pulse of laterally extensive sills as the critical deadly interval. Heat from these sills exposed untapped volatile-fertile sediments to contact metamorphism, likely liberating the massive greenhouse gas volumes needed to drive extinction. These observations suggest that large igneous provinces characterized by sill complexes are more likely to trigger catastrophic global environmental change than their flood basalt- and/or dike-dominated counterparts.Author SummaryAlthough the mass end-Permian extinction is linked to large igneous provinces, its trigger remains unclear. Here, the authors propose that the abrupt change from flood lavas to sills resulted in the heating of sediments and led to the release of large-scale greenhouse gases to drive the end-Permian extinction.

Emplacement of magma at shallow depth: insights from field relationships at Allan Hills, south Victoria Land, East Antarctica

Abstract Allan Hills nunatak, south Victoria Land, Antarctica, exposes an exceptional example of a shallow depth (< 500 m) intrusive complex formed during the evolution of the Ferrar large igneous province (LIP). Dyke distribution, geometries and relationships allow reconstruction of its history and mechanics of intrusion. Sills interconnect across host sedimentary layers, and a swarm of parallel inclined dolerite sheets is intersected by a radiating dyke-array associated with remnants of a phreatomagmatic vent, where the dolerite is locally quenched and mixed to form peperite. Intrusion geometries, and lack of dominant rift-related structures in the country rock indicate that magma overpressure, local stresses between mutually interacting dykes and vertical variations of host rock mechanical properties controlled the intrusive process throughout the thick and otherwise undeformed pile of sedimentary rocks (Victoria Group). Dolerite sills connected to one another by inclined sheets are inferred to record the preferred mode of propagation for magma-carrying cracks that represent the shallow portions of the Ferrar LIP plumbing system.

Fault-magma interactions during early continental rifting: seismicity of the Magadi-Natron-Manyara basins, Africa

Although magmatism may occur during the earliest stages of continental rifting, its role in strain accommodation remains weakly constrained by largely 2D studies. We analyze seismicity data from a 13-month, 39-station broadband seismic array to determine the role of magma intrusion on state-of-stress and strain localization, and their along-strike variations. Precise earthquake locations using cluster analyses and a new 3D velocity model reveal lower crustal earthquakes beneath the central basins and along projections of steep border faults that degas CO2. Seismicity forms several disks interpreted as sills at 6-10 km below a monogenetic cone field. The sills overlie a lower crustal magma chamber that may feed eruptions at Oldoinyo Lengai volcano. After determining a new ML scaling relation, we determine a b-value of 0.87 ± 0.03. Focal mechanisms for 65 earthquakes, and 13 from a catalogue prior to our array reveal an along-axis stress rotation of ∼60° in the magmatically active zone. New and prior mechanisms show predominantly normal slip along steep nodal planes, with extension directions ∼ N90°E north and south of an active volcanic chain consistent with geodetic data, and ∼ N150°E in the volcanic chain. The stress rotation facilitates strain transfer from border fault systems, the locus of early stage deformation, to the zone of magma intrusion in the central rift. Our seismic, structural, and geochemistry results indicate that frequent lower crustal earthquakes are promoted by elevated pore pressures from volatile degassing along border faults, and hydraulic fracture around the margins of magma bodies. Results indicate that earthquakes are largely driven by stress state around inflating magma bodies.

Shallow magma diversions during explosive diatreme-forming eruptions

The diversion of magma is an important mechanism that may lead to the relocation of a volcanic vent. Magma diversion is known to occur during explosive volcanic eruptions generating subterranean excavation and remobilization of country and volcanic rocks. However, feedbacks between explosive crater formation and intrusion processes have not been considered previously, despite their importance for understanding evolving hazards during volcanic eruptions. Here, we apply numerical modeling to test the impacts of excavation and subsequent infilling of diatreme structures on stress states and intrusion geometries during the formation of maar–diatreme complexes. Explosive excavation and infilling of diatremes affects local stress states which inhibits magma ascent and drives lateral diversion at various depths, which are expected to promote intra-diatreme explosions, host rock mixing, and vent migration. Our models demonstrate novel mechanisms explaining the generation of saucer-shaped sills, linked with magma diversion and enhanced intra-diatreme explosive fragmentation during maar-diatreme volcanism. Similar mechanisms will occur at other volcanic vents producing crater-forming eruptions.Shallow magmatic feeder systems in monogenetic volcanic fields may determine how a volcano erupts. Here, the authors use numerical modeling to show that explosive excavation and infilling of eruptive craters affects local stress states, with feedbacks controlling sites and depths of crater-forming explosions.

Lacustrine fan delta deposition alongside intrabasinal structural highs in rift basins: an example from the Early Cretaceous Jiuquan Basin, Northwestern China

Development of fan deltas alongside intrabasinal structural highs has been overlooked compared to those forming on basin margins. However, these fan deltas may provide important clues regarding the tectonic and climatic controls on deposition during rift development. This paper documents fan delta deposition alongside an intrabasinal structural high within the Early Cretaceous Xiagou Formation of the Jiuquan Basin, China, using subsurface geological and geophysical data. Deposits observed in drill core support fan delta deposition occurring almost exclusively through subaerial and subaqueous gravity flows. Subsurface mapping reveals a consistent decrease in the areal extent of fan deltas from lowstand to highstand system tracts, suggesting that deposition alongside the structural high is sensitive to lake-level changes. The temporal and spatial distribution of the fan deltas display retrogradational stacking patterns, where fan deltas exhibit a decreasing lateral extent up-sequence until fan delta deposition terminated and was replaced by deposition of fine-grained lacustrine deposits. The retrogradational stacking patterns observed alongside the intrabasinal structural high are not observed in fan deltas along the basin margin in the lower parts of the Xiagou Formation. Subsidence profiles also show differential subsidence across the basin during the earliest stages of this formation, likely resulting from border fault movements. These data suggest that non-uniform stacking patterns in the lower parts of the Xiagou Formation reflect basin-scale tectonic movements as the dominant control on synrift deposition patterns. However, later stages of Xiagou Formation deposition were characterized by uniform subsidence across the basin, and uniform retrogradational stacking patterns for fan deltas alongside the intrabasinal structural high and border fault. These observations suggest that basin-scale tectonic movements played a relatively limited role in controlling sediment deposition, and imply a potential change to regional-scale processes affecting fan delta deposition during later synrift stages. Climate change is favored here as the region-scale control on the uniform retrogradational fan delta stacking patterns. This assertion is supported by pollen assemblages, isotope signatures, and organic geochemical analyses, which collectively suggest a change from a humid to semi-arid environment during later synrift stages. We suggest that variations in stacking patterns between different fan delta systems can provide insights into the basin- and regional-scale processes that control rift basin deposition.