Basaltic Eruptions Research Articles

Magmatic and phreatomagmatic contributions on the ash-dominated basaltic eruptions: Insights from the April and November–December 2005 paroxysmal events at Karthala volcano, Comoros

Basaltic eruptions are commonly associated with lava emissions and relatively weak explosive activities, but they can sometimes produce strong explosive eruptive phases. In April and November 2005, two paroxysmal eruptive events occurred within the summit crater of Karthala basaltic shield volcano (Grande Comore Island, Comoros), which hosted a water lake before each of these events. Both 2005 ash plumes spread across the Comoros Archipelago and heavily impacted the whole Grande Comore Island. Associated deposits on the volcano summit are extremely fine-grained (up to 50 wt% of fine ash <63 μm for some analyzed layers) and rich in millimeter-sized rounded accretionary lapilli aggregates. Field observations, as well as textural and chemical analyses performed on both coarse- and fine-grained pyroclasts permit to identify juvenile and non-juvenile components and quantify their peculiar characteristics. Coarse ash (710–1000 μm) mainly consists of juvenile pumice particles (vesicle number density NV = 4.5 104 mm−3 and gas to melt ratio VG/VL = 1.5, on average), characterized by glassy groundmasses and representative of magma portions ascending quickly within the eruptive conduits (up to 10 m s−1). A relatively low amount of juvenile scoria particles are also observed in the coarse ash fractions, which are characterized by magma degassing (NV = 4.7 104 mm−3 and VG/VL = 0.5 on average) and associated crystallization (occurrence of dendritic microlites). Non-juvenile fragments (from blocks to coarse ash) are dense lava or intrusive fragments. Their amount decreases exponentially towards the fine ash fractions, which are mainly composed of juvenile, blocky, dense and glassy particles that are characterized by unambiguous textural signs of brittle fragmentation (hackle lines, stepped features and cracks). We support that Molten Fuel-Coolant Interactions between highly porous fast ascending basaltic magmas and external waters occurred during the paroxysmal phases of the studied eruptions, leading to a brittle-dominant and efficient regime of magma fragmentation. Variable but large amount of fine ash grains through the stratigraphic depth of the deposits can be ascribed to the brittle failure of the vesicle walls of the initial porous magma. Concurrently, thermohydraulic explosions caused the host rock fragmentation at shallow level, generating the relatively coarse non-juvenile particles. A short-lived episode of intense lava fountaining associated with steam explosions eventually occurred at the end of the November 2005 paroxysm, forming the last and relatively coarse tephra layer at the top of the studied eruptive sequence. Each paroxysmal phase lasted about a day as each associated water lake and shallow water table progressively vaporized and dried away. Both eruptions ended with lava pond and weak lava fountaining activities confined within the summit crater. We conclude that the contributions of both magmatic processes and phreatomagmatic interaction mechanisms ultimately generated the grain size, grain component and grain texture variabilities observed within the paroxysmal deposits. This work contributes to a better understanding of the generation of unusual fine ash from basaltic explosions as well as their eruptive dynamics and associated mechanisms, from magma ascent in the conduit to the fragmentation level and the interaction with intra-crateric lake waters.

Linking exhumation, paleo-relief, and rift formation to magmatic processes in the western Snake River Plain, Idaho, using apatite (U-Th)/He thermochronology

Abstract The western Snake River Plain (WSRP) in southwest Idaho has been characterized as an intracontinental rift basin but differs markedly in topography and style from other Cordilleran extensional structures and structurally from the down-warped lava plain of the eastern Snake River Plain. To investigate mechanisms driving extension and topographic evolution, we sampled granitoid bedrock from Cretaceous and Eocene-aged plutons from the mountainous flanks of the WSRP to detail their exhumation history with apatite (U-Th)/He (AHe) thermochronometry. AHe cooling dates from seventeen samples range from 7.9 ± 1.4 Ma to 55 ± 10 Ma. Most cooling dates from Cretaceous plutons adjacent to the WSRP are Eocene, while Eocene intrusions from within the Middle Fork Boise River canyon ~35 km NE of the WSRP yield Miocene cooling dates. The AHe dates provide evidence of exhumation of the Idaho batholith during the Eocene, supporting a high relief landscape at that time, followed by decreasing relief. The Miocene AHe dates show rapid cooling along the Middle Fork Boise River that we take to indicate focused river incision due to base level fall in the WSRP. Eocene AHe dates limit magnitudes of exhumation and extension on the flanks of the WSRP during Miocene rift formation. This suggests extension was accommodated by magmatic intrusions and intrabasin faults rather than basin-bounding faults. We favor a model where WSRP extension was related to Columbia River Flood Basalt eruption and enhanced by later eruption of the Bruneau-Jarbidge and Twin Falls volcanic fields, explaining the apparent difference with other Cordilleran extensional structures.

Intracaldera explosions and lava emissions during the 2007 caldera collapse of Piton de la Fournaise, La R\xe9union Island

The March–April 2007 Piton de la Fournaise basaltic eruption was the most significant eruption on La Réunion Island in historical times. On 2 April, a fissure opened on the southeastern flank of the volcano. Vigorous fountains fed lavas that rapidly reached the coast. Three days later, on the 5–6 April, major caldera collapse occurred at the summit, affecting the floor and walls of Dolomieu caldera. Monitoring records, primarily webcam images, have been analysed and integrated with geophysical data to reconstruct the chronology of events at the summit during caldera collapse. Those events included progressive subsidence of the former caldera floor, landslides, explosions, lava emissions and steam fumaroles, and lasted until 19 April though diminished greatly in frequency after 7 April. For two days after the main caldera collapse on 5 April, subsidence increments, intracaldera lava emission and explosions were closely associated in time and in location. Abundant steam and wet talus on the caldera walls imply that the shallow hydrothermal system and/or groundwater were exposed by subsidence. The presence of juvenile components in ash deposited at the summit during caldera collapse and the close link between intracaldera lavas and explosions suggest that many explosions were phreatomagmatic. Although caldera collapse was related to magma withdrawal from beneath the summit via the flank vent activity, numerous intracaldera lava emissions indicated that magma was nevertheless present at the summit during caldera collapse. We infer that the lava emissions were fed by an intrusion emplaced at the end of March into the eastern summit region and that the intrusion was being actively recharged during caldera collapse. Caldera collapse involved a complex combination of magma withdrawal and magma replenishment at separate summit reservoirs.

A red bole zircon record of cryptic silicic volcanism in the Deccan Traps, India

Abstract Silicic magmas within large igneous provinces (LIPs) are understudied relative to volumetrically dominant mafic magmas despite their prevalence and possible contribution to LIP-induced environmental degradation. In the 66 Ma Deccan LIP (India), evolved magmatism is documented, but its geographic distribution, duration, and significance remain poorly understood. Zircons deposited in weathered Deccan lava flow tops (“red boles”) offer a means of indirectly studying potentially widespread, silicic, explosive volcanism spanning the entire period of flood basalt eruptions. We explored this record through analysis of trace elements and Hf isotopes in zircon crystals previously dated by U–Pb geochronology. Our results show that zircon populations within individual red boles fingerprint distinct volcanic sources that likely developed in an intraplate setting on cratonic Indian lithosphere. However, our red bole zircon geochemical and isotopic characteristics do not match those from previously studied silicic magmatic centers, indicating that they must derive from yet undiscovered or understudied volcanic centers associated with the Deccan LIP.

Geochronology and petrogenesis of the Shuizi pyroxenite at Danba, western Yangtze Craton

扬子克拉通西缘在~260Ma发生短期内大规模峨眉山玄武岩溢流喷发。攀西地区发育的镁铁-超镁铁质岩被广泛认为是峨眉山大火成岩省的产物，但在北端松潘-甘孜岩区一直缺乏该类岩石的报道。本文首次报道扬子西缘丹巴水子乡单斜辉石岩的准确年龄，其锆石LA-ICP-MS U-Pb加权平均年龄为260.7±3.3Ma，表明其为峨眉山大火成岩省北端松潘-甘孜岩区镁铁-超镁铁质岩的组成部分。通过与攀枝花钒钛磁铁矿含矿岩体边缘相带苦橄岩和上部相带浅色辉长岩进行锆石微量元素对比显示，水子乡单斜辉石岩具有相近的高氧逸度，其ΔQFM为0~3，Ce_N/Ce_N^*平均为~30，该性质可能同样源自扬子西缘洋壳板片俯冲交代形成的较高氧逸度地幔源区。尽管如此，水子乡辉石岩体并未因高氧逸度而有明显的含钛磁铁矿饱和结晶，可能由其较低结晶分异程度造成。相比之下，攀枝花岩体经历了更高程度的含钛磁铁矿和斜长石分离结晶作用，伴随大规模的钒钛磁铁矿成矿。;Large-scale eruption of Emeishan flood basalt occurred in the western margin of the Yangtze Craton in a short period during~260Ma. Mafic-ultramafic rocks found in the Panxi area are widely considered to be products of the Emeishan Large Igneous Province (ELIP). However, there is a lack of reports on this kind of rocks in the Songpan-Ganzi area in the northern end of the ELIP. This paper first documents a precise age of the Shuizi pyroxenite in the Danba County. The zircon LA-ICP-MS U-Pb weighted mean age is 260.7±3.3Ma, indicating a representative mafic-ultramafic constituent of the northern end of the ELIP in the Songpan-Ganzi area. Comparing zircon trace elements of the Shuizi pyroxenite with those of picrite in the marginal zone and gabbro in the upper zone of the Panzhihua Fe-Ti-V oxide deposit were conducted. It is found that the Shuizi pyroxenite has similarly high magmatic oxygen fugacity, ΔQFM of 0~3 and average Ce_N/Ce_N^* of~30, to the Panzhihua mafic rocks. The high oxygen fugacity for both is considered to be derived from a mantle lithosphere source which had been metasomatized by subducted oceanic slab materials in the western margin of the Yangtze Craton. Although the Shuizi pyroxenite has high oxygen fugacity, there is less obvious Fe-Ti oxide saturation and precipitation here. This could be attributed to its relatively low degree of crystallization and differentiation. In contrast, the Panzhihua mafic-ultramafic intrusion experienced higher degree of titanomagnetite and plagioclase fractional crystallization, accompanied by a massive Fe-Ti-V oxide mineralization.

Cooler Equatorial Climate in the Late Lopingian Estimated from Paleosols Developed on Emeishan Basalts

Paleosols that developed on large igneous province flood basalts provide a valuable opportunity to quantitatively reconstruct the terrestrial paleoclimate during a given period. This study uses a newly constructed basaltic soil–climate function to obtain terrestrial paleoclimatic data following the termination of the Emeishan basalt eruption in the late Lopingian (∼253 Ma). The constructed relationships are mean annual temperature of 32.81e−0.4166X (R2=0.55, standard error [SE]=4.6°C) and mean annual precipitation of −875.5ln(X)+1792 (R2=0.69, SE=328 mm), where X=SiO2/(Al2O3+Fe2O3) (the Saf index). We collected samples from nine paleosol profiles developed on the last-phase Emeishan basalts and conducted field observations as well as micromorphological, mineral composition, and geochemical analyses. On the basis of these paleosols, we reconstructed the regional mean temperature and precipitation to be 18.9°C±4.6°C and 1546±328 mm, respectively. The regional climate was cooler than those in modern equatorial regions of similar latitudes. The reduced CO2 concentration in the atmosphere resulting from rapid weathering of Emeishan basalts in a tropical climate could have played a dominant role in causing the regional cool terrestrial paleoclimate following the termination of the Emeishan basalt eruption.

Soil Carbon Dioxide Planetary Thermostat.

Biological regulation of planetary temperature has been explained with the Daisyworld model, in which reflective-cooling white daises balance absorbing-warming black daisies. This article advances the proposition that cooling "daisies" of Daisyworld represent carbon sequestration and consumption by productive soils and ecosystems, such as grasslands expanding into deserts and tropical forests migrating toward the poles. On the other hand, warming "daisies" represent continued CO2 emissions from volcanoes and springs allowed by unproductive frigid and desert ecosystems. Greenhouse spikes of CO2 in deep time from large perturbations, such as flood basalt eruptions and asteroid impacts, did not continue as lethal runaway greenhouses, such as Venus, nor did low CO2 of ice ages decline to a sterile global snowball, such as Mars. These hypotheses are quantified and tested by new global soil maps derived from paleosols of the last extremes of atmospheric CO2: middle Miocene (16 Ma) and last glacial maximum (20 ka), when CO2 levels were 588 ± 72 and 180 ppm, respectively. Observed expansion of productive soils curbed large atmospheric injections of CO2 in deep time and observed expansion of unproductive soils during ice ages of low CO2 was thwarted by continued metamorphic and volcanic degassing. This short-term Soilworld thermostat of biogeographic redistribution of ecosystems supplemented long-term evolution of terrestrial carbon sequestration curbing solar radiation increases over billions of years. Similar agricultural management of ecosystems has potential for short-term carbon sequestration.

Global trends in novel stable isotopes in basalts: Theory and observations

The geochemistry of global mantle melts suggests that both mid-ocean ridge basalts (MORB) and ocean island basalts (OIB) sample lithological and temperature heterogeneities originating in both the upper and lower mantle. Recently, non-traditional stable isotopes have been suggested as a new tool to complement existing tracers of mantle heterogeneity (e.g., major and trace elements, radiogenic isotopes), because mineral- and redox-specific equilibrium stable isotope fractionation effects can link the stable isotope ratios of melts to their source mineralogy and melting degree. Here, we investigate five stable isotope systems (Mg-Ca-Fe-V-Cr) that have shown promise in models or natural samples as tracers of mantle temperature and/or lithological heterogeneity. We use a quantitative model, combining thermodynamically self-consistent mantle melting and equilibrium isotope fractionation models, to explore the behaviour of the isotope ratios of these elements during melting of three mantle lithologies (peridotite, and silica-excess and silica-deficient pyroxenites), responding to changes in mantle mineralogy, oxygen fugacity, temperature and pressure. We find that, given current analytical precision, the stable isotope systems examined here are not predicted to be sensitive to mantle potential temperature variations through equilibrium isotope fractionation processes. By contrast, source lithological heterogeneity is predicted to be detectable in some cases in the stable isotope ratios of erupted basalts, although generally only at proportions of >10% MORB-like pyroxenite in the mantle source, given current analytical precision. Magnesium and Ca stable isotopes show most sensitivity to a garnet-bearing source lithology, and Fe and Cr stable isotopes are potentially sensitive to the presence of MORB-like pyroxenite in the mantle source, although the behaviour of Cr isotopes is comparatively under-constrained and requires further work to be applied with confidence to mantle melts. When comparing the magnitude and direction of predicted equilibrium isotopic fractionation of peridotite and pyroxenite melts to natural MORB and OIB data, we find that aspects of the natural data (including the mean Mg-Ca-Fe-V isotopic composition of MORB, the range of Mg-Ca isotopic compositions seen in MORB data, the mean Mg-Ca-Cr isotopic composition of OIB, and the range of Mg-V-Cr isotopic compositions in OIB data) can be matched by equilibrium isotope fractionation during partial melting of peridotite and pyroxenite sources – with pyroxenite required even for some MORB data. However, even when considering analytical uncertainty on natural sample measurements, the range in stable isotope compositions seen across the global MORB and OIB datasets suggests that kinetic isotope fractionation, or processes modifying the isotopic composition of recycled crustal material such that it is distinct from MORB, may be required to explain all the natural data. We conclude that the five stable isotope systems considered here have potential to be powerful complementary tracers to other geochemical tracers of the source lithology of erupted basalts. However, continued improvements in analytical precision in conjunction with experimental and theoretical predictions of isotopic fractionation between mantle minerals and melts are required before these novel stable isotopes can be unambiguously used to understand source heterogeneity in erupted basalts.

New Estimates on the Basalt Volume of the Tarim (Not So Large) Igneous Province, NW China

AbstractThe early Permian basalts of the Tarim Basin of NW China are geographically widespread, covering over 200,000 km2 and are in places up to 1‐km thick and hence frequently referred to represent a Large Igneous Province (LIP). Available volume estimates are based on average lava thickness multiplied by area, using these values results in volumes exceeding 100,000 km3, and is commonly cited as being evidence of LIP status. Using a database of 58 boreholes, field sections, and seismic data we calculate a range of basalt isopach maps, and use these to estimate new, more refined overall eruption volumes. These volumes are significantly smaller than previously reported estimates. We demonstrate that the thickness distribution of erupted basalts is nonuniform and strongly skewed, indicating that volume estimates using an arithmetically averaged thickness of input data will result in overestimation of total volumes. We find that even with the addition of an estimated nonerupted volume component, representing cumulates and intrusions, that the total volume of the Tarim volcanics is significantly below the minimum threshold of a LIP. Using a simple Monte Carlo model we demonstrate that even while accounting for uncertainty and error in the input data, the probability of the total volume exceeding 100,000 km3 is very low. The new volume estimates imply a revision of the climate and environmental impact of the Tarim volcanism. We therefore suggest that the term Tarim Igneous Province should be used instead of Tarim Large Igneous Province.

A Model for Buoyant Tephra Plumes Coupled to Lava Fountains With an Application to the 29th of August 2011 Paroxysmal Eruption at Mount Etna, Italy

AbstractExplosive basaltic eruptions pose significant threats to local communities, regional infrastructures and international airspace. They produce tephra plumes that are often associated with a lava fountain, complicating their dynamics. Consequently, source parameters cannot be easily constrained using traditional formulations. Particularly, mass flow rates (MFRs) derived from height observations frequently differ from field deposit‐derived MFRs. Here, we investigate this discrepancy using a novel integral plume model that explicitly accounts for a lava fountain, which is represented as a hot, coarse‐grained inner plume co‐flowing with a finer‐grained outer plume. The new model shows that a plume associated with a lava fountain has higher variability in rise height than a standard plume for the same initial MFR depending on initial conditions. The initial grain‐size distribution and the relative size of the lava fountain compared to the surrounding plume are primary controls on the final plume height as they determine the strength of coupling between the two plumes. We apply the new model to the August 29, 2011 paroxysmal eruption of Mount Etna, Italy. The modeled MFR profile indicates that the field‐derived MFR does not correspond to that at the vent, but rather the MFR just above the lava fountain top. High fallout from the lava fountain results in much of the erupted solid material not reaching the top of the plume. This material deposits to form the proximal cone rather than dispersing in the atmosphere. With our novel model, discrepancies between the two types of observation‐derived MFR can be investigated and understood.

Miocene to sub-Recent magmatism at the intersection between the Dead Sea Transform and the Ash Shaam volcanic field: evidence from the Yarmouk River gorge and vicinity

AbstractThe Yarmouk River gorge extends along the Israel–Jordan–Syria border junction. It marks the southern bound of the Irbid–Azraq rift and Harrat Ash Shaam volcanic field at their intersection with the younger Dead Sea Transform plate boundary. During the last ∼13 Ma, the gorge has repeatedly accumulated basaltic units, chronologically named the Lower, Cover, Yarmouk and Raqqad Basalt formations. We examined their origin and distribution through aerial photos, and geological and geophysical evidence. Our results define a southern Golan magmatic province, which includes exposed Miocene (∼13 Ma) basalts, gabbro–diabase intrusions below the gorge and the adjacent Dead Sea Transform valley, and numerous Pliocene–Pleistocene volcanic sources along the gorge. Cover Basalt (∼5.0–4.3 Ma) eruptions formed two adjacent 0–100 m thick plateaus on the transform shoulder before flowing downslope to fill the topographically lower Dead Sea Transform valley with ∼700 m thick basalts. Later incision of the Yarmouk River and displacement along its associated fault divided the plateaus and formed the gorge. The younger Yarmouk (0.8–0.6 Ma) and Raqqad (0.2–0.1 Ma) basalts erupted in the upper part of the gorge from volcanos reported here, and flowed downstream toward the Dead Sea Transform valley. Consequently, eruptions from six phreatic volcanic vents altered the Yarmouk River morphology from sinuous to meandering. Our results associate the ∼13 Ma long southern Golan volcanism with the proposed SW-trending extensional Yarmouk Fault, located east of the Dead Sea Transform. Hence, the Yarmouk volcanism is associated with the ongoing Harrat Ash Shaam activity, which is not directly linked to the displacement along the Dead Sea Transform.

Rapid eruption of the Emeishan continental flood basalts: New paleomagnetic and geochronologic constraints

The duration of the eruption of the Emeishan large igneous province is hotly debated. We conducted a magnetostratigraphic and geochronological study of the core area of the large igneous province in the Binchuan area of Yunnan Province, southwestern China, in order to constrain the duration of the eruption. The results of detailed thermal demagnetization experiments revealed two remanent magnetic components from the volcanic rocks of 11 composite sections. A low-temperature component separated below 300 °C is interpreted as a recent viscous remanence. Additionally, reliable characteristic remanent magnetizations were revealed above 400 °C, with unblocking temperatures up to 580−680 °C, which passed the fold test and record three magnetozones. Zircons from the felsic ignimbrites exposed in the final stage of the mafic volcanism are dated to 258.2 ± 0.7 Ma (n = 15; mean square of weighted deviates = 1.3) by sensitive high-resolution ion microprobe. Stratigraphic and magnetostratigraphic correlations of the Emeishan basalts in the Binchuan sections indicate that the eruption of the mafic rocks of the Emeishan large igneous province can be clearly divided into early (reverse polarity subzone), middle (normal polarity subzone), and late (reverse polarity subzone) stages, with a total duration of less than 1.7 m.y. (260.8−259.1 Ma). However, by combining this chronology with previously reported conodont biostratigraphic results from locations around the Emeishan large igneous province, and comparing the dominant normal-reverse polarity sequence in the Emeishan large igneous province with the geomagnetic polarity time scale, we obtain a much shorter duration of the main eruptive stage of &amp;lt;0.8 m.y. (260.4−259.6 Ma). About three quarters of the basalts of the Emeishan large igneous province record have a normal polarity and erupted within 0.4 m.y., while the other quarter, mainly distributed in the central zone, shows a reverse polarity and much shorter duration. Given the short duration of the eruption, gas volatiles would have been released into the atmosphere at high rates, which might provide a causal link between the rapid eruption and the end-Guadalupian mass extinction. Before the mantle plume eruption, localized eruptions probably occurred. After eruption of the mafic Emeishan flood basalts, an acid volcanic eruption occurred in the early Wuchiapingian, which was sporadically distributed in the Emeishan large igneous province.

The Diagenetic Alteration of the Carbonate Rocks from the Permian Qixia Formation as Response to Two Periods of Hydrothermal Fluids Charging in the Central Uplift of Sichuan Basin, SW China

The hydrothermal fluid–carbonate rock reaction is frequently regarded to occur in deep-burial diagenesis, and the hydrothermal dissolution is usually distributed and takes place along the faults. Previous studies have suggested that there was hydrothermal fluid activity locally in the Permian Qixia Formation in Sichuan Basin, likely related to the Emeishan basalt eruption. However, the effect of hydrothermal fluids on the carbonate rocks of the Qixia Formation in the central uplift of Sichuan Basin is still unclear. Based on the characteristics and geochemical parameters of the diagenetic minerals, this study aims to reveal the diagenetic alteration related to the hydrothermal fluid–rock reaction in the Qixia Formation and reestablish the diagenetic evolution by using the timing of diagenetic mineral precipitation. The methods include petrographic observation; trace and rare earth element (REE) analysis; C, O and Sr isotope measurement; fluid inclusion temperature measurement and cathodoluminescence analysis. According to the petrographic characteristics, the dolostones are mainly of crystalline structure, namely fine-medium crystalline dolostone, meso-coarse crystalline dolostone, and coarse crystalline dolostone, with the cathodoluminescence color becoming brighter in that order. The limestones from the Qixia Formation are of the bioclastic limestone type, with no cathodoluminescence color. Compared with dolostones, limestones have higher Sr content, lower Mn content, and heavier oxygen isotopes. With the crystalline size of dolostone becoming coarser, the oxygen isotopes of dolostones tend to become lighter. The meso-coarse crystalline dolostone has the highest Mn content and negative carbon isotope. Both limestones and dolostones have an obvious positive Eu anomaly in the Qixia Formation. However, the REE patterns of fine-medium crystalline dolostones are very different from those of meso-coarse crystalline dolostones. It is credible that there were two periods of hydrothermal fluid charging, with different chemical compositions. The first period of hydrothermal fluids could laterally migrate along the sequence boundary. Fine-medium crystalline dolostones were almost completely distributed below the sequence boundary and were dolomitized during the shallow burial period. As products of the hydrothermal fluid–dolostone reaction, the saddle-shaped dolomites in the meso-coarse crystalline dolostones were the evidence of the second period of hydrothermal fluids. As a result, the dolomitization model was established according to the timing of diagenetic mineral precipitation, which can improve that the geological understanding of the effect of hydrothermal fluid activities on the carbonate rocks in the Qixia Formation.

Timescales of magma transport in the Columbia River flood basalts, determined by paleomagnetic data

Flood basalts represent major events in Earth History, in part because they are linked to large climate perturbations and mass extinctions. However, the durations of individual flood basalt eruptions, which directly impact potential environmental crises, are poorly constrained. Here we use a combination of paleomagnetic data and thermal modeling to create a magnetic geothermometer (MGT) that can constrain the active transport lifetime of magmatic conduits and intrusions. We apply the MGT technique to eight feeder dike segments of the Columbia River basalts (CRB), demonstrating that some dike segments were actively heating host rocks for less than one month, while other segments may have been active for several years. Results suggest that eruption rates, localized spatially along-strike of dike segments, were as high as 1–8 km3day−1. These results help contextualize field evidence for contrasting CRB eruption durations and suggest a pathway for constraining the tempo of global flood basalt magmatism that is beyond the resolution of geochronology.

Mayotte seismic crisis: building knowledge in near real-time by combining land and ocean-bottom seismometers, first results

SUMMARYThe brutal onset of seismicity offshore Mayotte island North of the Mozambique Channel, Indian Ocean, that occurred in May 2018 caught the population, authorities and scientific community off guard. Around 20 potentially felt earthquakes were recorded in the first 5 d, up to magnitude Mw 5.9. The scientific community had little pre-existing knowledge of the seismic activity in the region due to poor seismic network coverage. During 2018 and 2019, the MAYOBS/REVOSIMA seismology group was progressively built between four French research institutions to improve instrumentation and data sets to monitor what we know now as an on-going exceptional submarine basaltic eruption. After the addition of 3 medium-band stations on Mayotte island and 1 on Grande Glorieuse island in early 2019, the data recovered from the Ocean Bottom Seismometers were regularly processed by the group to improve the location of the earthquakes detected daily by the land network. We first built a new local 1-D velocity model and established specific data processing procedures. The local 1.66 low VP/VS ratio we estimated is compatible with a volcanic island context. We manually picked about 125 000 P and S phases on land and sea bottom stations to locate more than 5000 events between February 2019 and May 2020. The earthquakes outline two separate seismic clusters offshore that we named Proximal and Distal. The Proximal cluster, located 10 km offshore Mayotte eastern coastlines, is 20–50 km deep and has a cylindrical shape. The Distal cluster start 5 km to the east of the Proximal cluster and extends below Mayotte's new volcanic edifice, from 50 to 25 km depth. The two clusters appear seismically separated, however our data set is insufficient to firmly demonstrate this.

Melting conditions and mantle source composition from probabilistic joint inversion of major and rare earth element concentrations

The chemical composition of erupted basalts provides a record of the thermo-chemical state of their source region and the melting conditions that lead to their formation. Here we present the first probabilistic inversion framework capable of inverting both trace and major element data of mafic volcanic rocks to constrain mantle potential temperature, depth of melting, and major and trace element source composition. The inversion strategy is based on the combination of (i) a two-phase multi-component reactive transport model, (ii) a thermodynamic solver for the evolution of major elements and mineral/liquid phases, (iii) a disequilibrium model of trace element partitioning and (iv) an adaptive Markov chain Monte Carlo algorithm. The mechanical and chemical evolution of melt and solid residue are therefore modelled in an internally- and thermodynamically-consistent manner.We illustrate the inversion approach and its sensitivity to relevant model parameters with a series of numerical experiments with increasing level of complexity. We show the benefits and limitations of using major and trace element compositions separately before demonstrating the advantages of a joint inversion. We show that such joint inversion has great sensitivity to mantle temperature, pressure range of melting and composition of the source, even when realistic uncertainties are assigned to both data and predictions. We further test the reliability of the approach on a real dataset from a well-characterised region: the Rio Grande Rift in western North America. We obtain estimates of mantle potential temperature (∼1340 °C), lithospheric thickness (∼60 km) and source composition that are in excellent agreement with numerous independent geochemical and geophysical estimates. In particular, this study suggests that the basalts in this region originated from a moderately hot upwelling and include the contribution from a slightly depleted source that experienced a small degree of melt or fluid metasomatism. This component is likely associated with partial melting of the lower portions of the lithosphere. The flexibility of both the melting model and inversion scheme developed here makes the approach widely applicable to assessing the thermo-chemical structure and evolution of the lithosphere-asthenosphere system and paves the way for truly joint geochemical-geophysical inversions.

GEOSCCM Columbia River Flood Basalt Eruption

GEOSCCM Columbia River Flood Basalt Eruption

Degassing and gas percolation in basaltic magmas

Due to their generally low eruptive melt viscosities and concomitant high diffusivities of volatiles, basaltic magmas degas relatively efficiently. This relative efficiency, combined with variations in style, extent, timing and length scales of degassing govern the range of eruptive styles observed at basaltic volcanoes. The result is a surprising complexity of degassing regimes and products in basaltic volcanism. In particular, the transition between closed- and open-system degassing at low pressure at the percolation threshold may strongly affect the type of eruption. Here we aim to better understand degassing and gas percolation processes in basaltic magmas and their implications for eruptive style. Combining new and literature data, we present a database of vesicle metrics in basaltic rocks including vesicularity, vesicle number density, vesicle size distribution (and its polydispersivity), vesicle connectivity and permeability. We combine these textural and petrophysical data with a numerical model of percolation for systems having polydisperse vesicle size distributions. Using this model, we also evaluate different definitions of vesicle connectivity inherent to different measurement techniques. Our results show that polydispersivity exerts a strong control on the percolation threshold of basaltic magmas and consequently on eruptive style. Intermediate to highly polydisperse bubble networks are more typical of Hawaiian activity and are characterized by higher values of percolation threshold. This results in delayed coalescence and an increase in magma vesicularity hindering the formation of large decoupled and buoyant bubbles, which in turn can promote magma acceleration, fragmentation by inertia below the percolation threshold and sustained fountaining activity. Bubble populations with lower polydispersivity, typical of Strombolian eruptions, promote early coalescence prior to fragmentation, which in turn may lead to the formation of large decoupled slugs or gas pockets and/or plugs at the surface via outgassing. Further, we discuss the implications of our findings for Plinian, violent Strombolian, Surtseyan, deep submarine and effusive basaltic eruptions.

Geochemical Constraints on the Structure of the Earth's Deep Mantle and the Origin of the LLSVPs

AbstractGeophysical analysis of the Earth's lower mantle has revealed the presence of two superstructures characterized by low shear wave velocities on the core‐mantle boundary. These Large Low Shear Velocity Provinces (LLSVPs) play a crucial role in the dynamics of the lower mantle and act as the source region for deep‐seated mantle plumes. However, their origin, and the characteristics of the surrounding deep mantle, remain enigmatic. Mantle plumes located above the margins of the LLSVPs display evidence for the presence of this deep‐seated, thermally and/or chemically heterogeneous mantle material ascending into the melting region. As a result, analysis of the spatial geochemical heterogeneity in ocean island basalts provides constraints on the structure of the Earth's lower mantle and the origin of the LLSVPs. In this study, we focus on the Galápagos Archipelago in the eastern Pacific, where bilateral asymmetry in the radiogenic isotopic composition of erupted basalts has been linked to the presence of LLSVP material in the underlying plume. We show, using spatial variations in the major element contents of high‐MgO basalts, that the isotopically enriched south‐western region of the Galápagos mantle—assigned to melting of LLSVP material—displays no evidence for lithological heterogeneity in the mantle source. As such, it is unlikely that the Pacific LLSVP represents a pile of subducted oceanic crust. Clear evidence for a lithologically heterogeneous mantle source is, however, found in the north‐central Galápagos, indicating that a recycled crustal component is present near the eastern margin of the Pacific LLSVP, consistent with seismic observations.

Top–down control on eruptive style at Masaya volcano inferred from melt composition

Highly explosive Plinian eruptions of basaltic magma are enigmatic because low melt viscosities should inhibit such eruptive style. Masaya volcano, Nicaragua, is a persistently active basaltic system capable of a wide range of eruptive styles, from open-conduit lava lake activity to voluminous Plinian eruptions; it is thus an ideal natural laboratory to constrain potential controls on basaltic eruption style. Here we report the major, trace, and volatile (CO2, H2O, S, Cl, F) element composition of olivine-, plagioclase- and clinopyroxene-hosted melt inclusions as well as matrix glasses from lava lake ejecta and two Plinian tephra deposits—the 2.1 ka Masaya Triple Layer and the 1.8 ka Ticuantepe Lapilli—to test whether pre-eruptive volatile contents and degassing history may be linked to eruptive style. All samples display a relatively narrow and largely overlapping basaltic–basaltic andesitic compositional range (51.7±1.0 wt.% SiO2, 4.8±0.4 wt.% MgO) with similar trace element signatures (e.g., Ce/Y=0.82±0.10, Ba/La=74±11). However, lava lake and Plinian samples show systematic differences in pre-eruptive volatile contents, forming distinct groups with mean H2O contents of 0.6±0.2 wt.% (lava lake), 1.1±0.2 wt.% (Masaya Triple Layer), and 1.9±0.3 wt.% (Ticuantepe Lapilli). Together, these groups generate broad positive correlations between S, Cl and H2O concentrations, with maximum values reaching 920 ppm, 1300 ppm and 2.3 wt.%, respectively, which are low compared to typical Central American arc magmas. Magma temperature estimates overlap and average at 1115±30°C, while volatile saturation pressures are low, mainly <100 MPa, although only lava lake samples record pressures <31 MPa. These observations reiterate the compositionally buffered state of the volcano's magmatic system highlighted by previous work and demonstrate that — regardless of eruption style — all Masaya magmas undergo variable, but extensive, pre-eruptive degassing at low pressure. Geohygrometry, gas emissions, and H2O/Ce–Ba/La systematics suggest initial, undegassed H2O contents on the order of 3.9–5.5 wt.%. Our results imply that pre-eruptive volatile contents are not the culprit for Plinian events at Masaya. Instead, we propose that the volcano's vigorous magma supply is modulated in a top–down manner to produce a wide range of eruptive styles, whereby temporary sealing of the conduit may instigate a transition to explosive behavior. In this model, rapid magma ascent is triggered when the seal eventually breaks from degassing-induced pressurization, yielding high degrees of undercooling, rapid microlite growth, and a dramatic increase in magma viscosity and explosive eruption potential. There may thus be a thin line between open-conduit conditions and Plinian eruptions at Masaya.