Melt Lens Research Articles

Pleistocene to recent evolution of Mocho-Choshuenco volcano during growth and retreat of the Patagonian Ice Sheet

Mocho-Choshuenco volcano (39.9°S, 72.0°W) produced ∼75 explosive eruptions following retreat of the >1.5-km-thick Patagonian Ice Sheet associated with the local Last Glacial Maximum (LGM, from 35 to 18 ka). Here, we extend this record of volcanic evolution to include pre- and syn-LGM lavas that erupted during the Pleistocene. We establish a long-term chronology of magmatic and volcanic evolution and evaluate the relationship between volcanism and loading/unloading of the Patagonian Ice Sheet via twenty-four 40Ar/39Ar and two 3He age determinations integrated with stratigraphy and whole-rock compositions of lava flows and glass compositions of tephra. Our findings reveal that the edifice is much younger than previously thought and preserves 106 km3 of eruptive products, of which 50% were emplaced immediately following the end of the penultimate glaciation and 20% after the end of the LGM. A period of volcanic inactivity between 37 and 26 ka, when glaciers expanded, was followed by the eruption of incompatible element-rich basaltic andesites. Several of these syn-LGM lavas dated between 26 and 16 ka, which crop out at 1500−1700 m above sea level, show ice contact features that are consistent with emplacement against a 1400- to 1600-m-thick Patagonian Ice Sheet. Small volume dacitic eruptions and two explosive rhyolitic eruptions dominate the volcanic output from 18 to 8 ka, when the Patagonian Ice Sheet began to retreat rapidly. We hypothesize that increased lithostatic loading as the Patagonian Ice Sheet grew prohibited dike propagation, thus stalling the ascent of magma, promoting growth of at least three discrete magma reservoirs, and enhancing minor crustal assimilation to generate incompatible element-rich basaltic andesitic to dacitic magmas that erupted between 26 and 17 ka. From an adjacent reservoir, incompatible element-poor dacites erupted from 17 to 12 ka. These lava flows were followed by the caldera-forming eruption at 11.5 ka of 5.3 km3 of rhyolite from a deeper reservoir atop which a silicic melt lens had formed and expanded. Subsequent eruptions of oxidized dacitic magmas from the Choshuenco cone from 11.5 to 8 ka were followed by andesitic to dacitic eruptions at the more southerly Mocho cone, as well as small flank vent eruptions of basaltic andesite at 2.5 and 0.5 ka. This complex history reflects a multi-reservoir plumbing system beneath Mocho-Choshuenco, which is characterized by depths of magma storage, oxidation states, and trace element compositions that vary over short periods of time (<2 k.y.).

Blind magmatism abets nonvolcanic continental rifting

Tectonic forces alone cannot drive rifting in old and thick continental lithosphere. Geodynamic models suggest that thermal weakening is critical for lithospheric extension, yet many active rifts lack volcanism, seeming to preclude this process. We focus on one such rift, the Tanganyika-Rukwa segment of the East African Rift System, where we analyze local seismicity for shear wave anisotropy and couple the results with numerical modeling. The strongest splitting measurements are from earthquakes with paths sampling lower crustal regions of high compressional-to-shear wave velocity ratios and have fast polarization directions parallel to the local mantle flow, implying the existence of oriented melt lenses. This lower crustal magmatism and observed high surface heat flow are consistent with substantial lithospheric weakening and explain the enigmatic relief and increasing strain accommodation along the rift axis. We conclude that progressive nonvolcanic rifting is assisted by deep crustal melts yet to breach the surface.

Organic Input to Titan's Subsurface Ocean Through Impact Cratering.

Titan has an organic-rich atmosphere and surface with a subsurface liquid water ocean that may represent a habitable environment. In this work, we determined the amount of organic material that can be delivered from Titan's surface to its ocean through impact cratering. We assumed that Titan's craters produce impact melt deposits composed of liquid water that can founder in its lower-density ice crust and estimated the amount of organic molecules that could be incorporated into these melt lenses. We used known yields for HCN and Titan haze hydrolysis to determine the amount of glycine produced in the melt lenses and found a range of possible flux rates of glycine from the surface to the subsurface ocean. These ranged from 0 to 1011 mol/Gyr for HCN hydrolysis and from 0 to 1014 mol/Gyr for haze hydrolysis. These fluxes suggest an upper limit for biomass productivity of ∼103 kgC/year from a glycine fermentation metabolism. This upper limit is significantly less than recent estimates of the hypothetical biomass production supported by Enceladus's subsurface ocean. Unless biologically available compounds can be sourced from Titan's interior, or be delivered from the surface by other mechanisms, our calculations suggest that even the most organic-rich ocean world in the Solar System may not be able to support a large biosphere.

Subduction of solidifying melt lens beneath ultraslow-spreading ridges revealed by (Grt-)amphibolites and dolerites in the Xiugugabu ophiolitic mélange, southwestern Tibet

Ridge subduction is an inevitable process during the evolution history of a vanished ocean. Due to its special characteristics, such as buoyant density, high heat fluxes and active magmatism, subduction of a ridge segment would leave specific fingerprints of high geothermal gradients and unique metamorphism. Therefore, investigations on ridge subduction could provide invaluable details of the evolution of oceanic lithosphere, and also help to understand the emplacement and preservation mechanism of worldwide ophiolites. This study mainly focuses on the mafic blocks of garnet-(Grt-)free and Grt-bearing amphibolites, foliated dolerites and dolerite dykes intruding into other blocks outcropped in the Xiugugabu ophiolitic mélanges of the west segment of the Yarlung-Tsangpo Suture Zone (YTSZ). On the basis of the detailed field, petrographical and geochemical investigations, combined with previous research on the other YTSZ ophiolites, this study highlights that the (grt-)amphibolites in the ophiolitic mélange are metamorphosed products under conditions of high geothermal gradients with protoliths to be subducted cumulate layered gabbros, which were likely to be newly crystallized in solidifying melt lenses beneath ridges. Hence, the generation of the XGGB ophiolite was contemporary with Neo-Tethyan ridge subduction. Intruding dolerite dykes show identical UPb ages to their intruded metamorphic soles, indicating that these amphibolites had gone through extremely rapid exhumation after peak conditions. Finally, scenarios of ridge destruction and subduction were proposed to complete the evolution history of the Neo-Tethyan Ocean and to provide an explanation for ophiolite emplacement and extremely rapid exhumation of metamorphic soles.

Beyond spreading rate: Controls on the thermal regime of mid-ocean ridges.

The thermal state of mid-ocean ridges exerts a crucial modulation on seafloor spreading processes that shape ~2/3 of our planet's surface. Standard thermal models treat the ridge axis as a steady-state boundary layer between the hydrosphere and asthenosphere, whose thermal structure primarily reflects the local spreading rate. This framework explains the deepening of axial melt lenses (AMLs)-a proxy for the basaltic solidus isotherm-from ~1 to ~3 km from fast- to intermediate-spreading ridges but fails to account for shallow crustal AMLs documented at slow-ultraslow spreading ridges. Here, we show that these can be explained by a numerical model that decouples the potentially transient ridge magma supply from spreading rate, captures the essential physics of hydrothermal convection, and considers multiple modes of melt emplacement. Our simulations show that melt flux is a better thermal predictor than spreading rate. While multiple combinations of melt/dike emplacement modes, permeability structure, and temporal fluctuations of melt supply can explain shallow crustal AMLs at slow-ultraslow ridges, they all require elevated melt fluxes compared to most ridge sections of comparable spreading rates. This highlights the importance of along-axis melt focusing at slow-ultraslow ridges and sheds light on the natural variability of their thermal regimes.

Stable Brine Layers beneath Europa’s Chaos

The formation mechanism of Europa’s large chaos terrain (>∼100 km diameter) and associated lenticulae (<∼10 km diameter) has been debated since their observations by the Galileo spacecraft. Their geomorphology and distribution suggest there may be reservoirs of saline liquid water 1–3 km beneath the surface—the “shallow water” model—generated by injection of ocean water or melting of the ice shell. Recent investigations on the evolution of small shallow-water bodies (≤103 km3) suggests that salts with a small effect on melting point (MgSO4) can extend the lifetime of saline bodies by ∼5% compared to freshwater reservoirs. However, sodium chloride, identified as a potential oceanic salt, has a significantly stronger impact on the freezing point, suggesting a further extension of liquid lifetimes. Moreover, the substantial volumes of liquid water (∼104 km3) beneath large chaos could be melted in situ rather than injected through a fracture, implying a distinct chemistry and formation environment. Here, we use numerical models to explore how the chemistry and disparate origins of shallow water control its evolution and lifetime. For small, injected sills, we find that NaCl can extend their liquid lifetime to ∼140 kyr—up to a ∼60% increase over freshwater sills. Saline melt lenses will last at least 175 kyr but, in contrast to sills, may persist as a stable layer of brine beneath the surface for over 500 kyr. Our results provide further support for the presence of liquid water at shallow depths within Europa’s ice shell today.

Two Distinct Magma Storage Regions at Ambrym Volcano Detected by Satellite Geodesy

AbstractThe flux of eruptible magma into a magmatic plumbing system influences eruption size and timing. If magma transfer is possible between two hydraulically‐connected magma lenses, system destabilization can tap a larger magma volume than stored in any one melt lens. This study identifies two distinct magma reservoirs beneath Ambrym, a basaltic island volcano in Vanuatu, during the time period February 2019 to January 2022. Using InSAR time series and a data assimilation approach, we estimate pressure changes within two reservoirs (located 5–7 and 4–6 km b.s.l.). Furthermore, a theoretical model demonstrates that the reservoirs may not currently be hydraulically connected, despite evidence of physical mixing of magma derived from each reservoir during the December 2018 eruption. These findings further our understanding of how magmatic plumbing systems at basaltic calderas may change after rift‐zone eruptions.

Mush Amalgamation, Short Residence, and Sparse Detectability of Eruptible Magma Before Andean Super‐Eruptions

AbstractGiant volcanic eruptions have the potential to overturn civilizations. Yet, the driving mechanism and timescale over which batholithic magma reservoirs transition from non‐eruptible crystal mush to mobile melt‐dominated stages and our capacity to detect a pending super‐eruption remain obscure. Here we show, using Sr isotope zonation in plagioclase crystals from three Andean large‐magnitude eruptions (Atana, Toconao, and Tara ignimbrites), that eruptible magma forms by amalgamation of isotopically diverse crystal populations and silicic melt without large‐scale reheating. In each case, crystals record large isotopic diversity in crystal cores, converging toward a common value in crystal rims that coincides with the composition of the rhyolitic carrier melt. Using diffusion chronometry, we show that the assembled magma resided pre‐eruptively in the crust for timescales of no more than decades to centuries for Atana and Tara, but up to several millennia for Toconao. These timescales and isotopic observations are consistent with the accumulation and destabilization of melt‐rich layers in crystal mush. While the prospect of capturing such melt lenses with most geophysical monitoring techniques is pessimistic, gravity modeling indicates that such structures are potentially resolvable. Our findings provoke a new assessment of the origin and hazards associated with large magnitude explosive eruptions.

The spreading rate dependence of the distribution of axial magma lenses along mid-ocean ridges

Seismically imaged axial melt lenses (AMLs) are seen almost everywhere along the axis of fast-spreading ridges but at only a few localized segment centers on slow-spreading ridges. Standard models assuming that AMLs form when melt percolating upward pools where freezing produces an impermeable cap do not explain this fundamental observation. To tackle this long-standing problem, we combine a crustal density model and a thermal model with a recent mechanical model for sill formation. The mechanical model predicts that AMLs form below the axial lithosphere but only if the average density of the axial brittle lithosphere is not greater than the magma density. For standard thermal models, crustal density structures inferred from seismic velocity data and normal crustal thicknesses, AMLs are found to be stable along all of a ridge segment for spreading rates greater than about 50 mm/y. To explain slow-spreading observations, we assume that a share of the melt produced by the mantle upwelling all along a segment is focused to the segment center. Some of this melt partially crystallizes, releasing latent heat, before the evolved magma flows along the axis to build the crust away from the segment center. This "extra" heat, beyond what is supplied by the magma that builds the crust near the segment center, results in the lithosphere thin enough for stable melt lenses at the segment center. Our results are consistent with observations and offer a quantitative explanation of the marked difference in the distribution of AMLs along fast- versus slow-spreading centers.

ICDP Oman Drilling Project: varitextured gabbros from the dike–gabbro transition within drill core GT3A

Abstract. The Oman ophiolite (Samail massif, Sultanate of Oman) is the largest sub-aerial exposure of oceanic lithosphere on Earth and provides the opportunity to study the accretion and alteration of oceanic lithosphere formed under fast-spreading conditions. Drill hole GT3A (23∘06′50.7′′ N, 58∘12′42.2′′ E) of the ICDP (International Continental Scientific Drilling Program) Oman Drilling Project with a length of 400 m aimed at penetrating the dike–gabbro transition of the Samail ophiolite paleocrust in order to shed light on the role of the axial melt lens (AML) during accretion of the lower plutonic crust. AMLs beneath fast-spreading mid-ocean ridges are sandwiched between the sheeted dike complex and the uppermost gabbros and are believed to feed the upper crust and, at least partially, the underlying crystal mush. Typical gabbroic rocks from dike–gabbro transitions of fast-spreading systems are the so-called “varitextured gabbros”, often showing considerable variations in mineral mode, texture and grain size, which are regarded as the frozen fillings of axial melt lenses. Here, we present a detailed petrographic, microanalytical and bulk-chemical investigation of 36 mafic rocks from the drill hole GT3A, which represent mostly varitextured gabbros, revealing a complex formation with several evolution stages. Poikilitic domains formed first, corresponding to an early crystallization stage, where only plagioclase and clinopyroxene of more primitive composition crystallized. Later, domains of granular textures containing also interstitial amphibole and Fe–Ti oxide were formed. This stage is characterized by a magma evolution that underwent crystal fractionation established by lower temperatures due to more efficient hydrothermal cooling at the margin of the AML. A last stage is characterized by pervasive hydrothermal alteration, where all primary minerals have been altered under temperature conditions, varying from the magmatic regime down to greenschist facies. A highlight of this stage is amphiboles showing noticeable compositional zoning. The observation of peculiar microgranular domains, representing relics of stoped exogenic material from the sheeted dike complex, documents the upward migration of an AML in a replenishment event, forcing the AML to burn through previously altered sheeted dikes. This process is responsible for significant assimilation of hydrothermally altered components, indicated by a marked Cl enrichment in the outer zones of magmatic amphiboles. Petrological modeling involving gabbros and basalts revealed that the GT3A rock suite followed a fractional crystallization evolution trend, with a primitive MORB as parental melt with an estimated water content of 0.2 wt % to 0.8 wt %. The modeled liquid lines of descent suggest a magmatic evolution via fractional crystallization, where the basalts correspond to frozen liquids, while the gabbros, especially the more primitive ones, show a significant cumulate component.

Temperatures and Timescales of Crystallization in the “Mush Zone” Beneath Fast‐Spreading Ocean Ridges

AbstractQuantifying the timescales of magma solidification and eruption beneath mid‐ocean ridges is critical to understanding how melts are emplaced and crystallized, and how long these magmas and their crystals remain in a semi‐molten “mush‐like” state. Here, we present a new method to quantify the crystallization and solidification temperatures of individual gabbros in the lower crust of the Wadi Tayin massif (Samail ophiolite). Cotectic crystallization temperatures of olivine + plagioclase + clinopyroxene are constrained by our clinopyroxene Mg# thermometer, and supported by clinopyroxene‐olivine Fe‐Mg equilibration temperatures. We find that lower crustal gabbros reach the cotectic at approximately 1,160 ± 16°C, and are fully solidified between 781 and 1011°C, depending on cooling rate. Our results show that lower crustal gabbros spend, on average, 7,500–11,500 years in a partially molten state with no systematic change in duration with depth in the crust. Samples with very fast cooling rates result in solidification only a few tens of meters off‐axis from their emplacement location and likely indicate solidification near the distal edge of the on‐axis melt lens or in off‐axis melt lenses. Incorporating a uniform 1 km wide emplacement region directly on‐axis, as inferred geophysically, and a half‐spreading rate of 55 mm/yr, our results suggest that the Mush Zone at fast‐spreading ridges remains partially molten for 0–3 km away from the ridge axis, but up to 4–7 km considering all possible sources of uncertainty. Both geophysical and geochemical estimates of the Mush Zone are inconsistent with conductive cooling of the lower oceanic crust.

Accretion Processes of Oceanic Crust in Slow‐spreading Ridges: Plagiogranite Perspective of the Xigaze Ophiolite, South Tibet

AbstractStructural and petrological data suggest that the Xigaze ophiolite from the Yarlung Zangbo Suture Zone (YZSZ) in south Tibet was a typical slow‐spreading ridge. A new field, geochemical, mineral, and U‐Pb zircon dataset of plagiogranite intrusions were used to constrain the dynamic processes of oceanic accretion in this slow‐spreading ridge. Plagiogranites mainly occur as dykes or intrusions intruded into the whole sequence of the ophiolite and have a similar orientation to the dolerite dykes developed in the late stage of detachment faulting. U‐Pb zircon ages of 122–123 Ma were obtained for two types of plagiogranites and associated dolerite dykes. Detailed geochemical and mineralogical examinations suggest that the plagiogranites are the product of low‐pressure (2–3 kbar) fractional crystallization of mid‐ocean ridge basalt‐like magma and unlikely to have been derived from the partial melting of hydrous gabbroic rocks. The complex cross‐cut relationship between the plagiogranites and ophiolite sequence reflects that they are controlled by small discontinued melt lenses rather than a big magma chamber under the ridge axis and reveals multiple injections during the oceanic crust accretion. The formation of plagiogranites possibly reflects the complex characteristic of oceanic accretion at slow‐spreading ridges, time‐dependent on structural (external) and magmatic (internal) processes.

Magmatic Activity and Dynamics of Melt Supply of Volcanic Centers of Ultraslow Spreading Ridges: Hints From Local Earthquake Tomography at the Knipovich Ridge

AbstractAlong ultraslow spreading ridges melt is distributed unequally, but melt focusing guides melt away from amagmatic segments toward volcanic centers. An interplay of tectonism and magmatism is thought to control melt ascent, but the detailed process of melt extraction is not yet understood. We present a detailed image of the seismic velocity structure of the Logachev volcanic center and adjacent region along the Knipovich Ridge. With travel times of P‐ and S‐waves of 3,959 earthquakes we performed a local earthquake tomography. We simultaneously inverted for source locations, velocity structure and the Vp/Vs‐ratio. An extensive low velocity anomaly coincident with high Vp/Vs‐ratios &gt;1.9 lies underneath the volcanic center at depths of 10 km below sea level in an aseismic area. More shallow, tightly clustered earthquake swarms connect the anomaly to a shallow anomaly with high Vp/Vs‐ratio beneath the basaltic seafloor. We consider the deep low‐velocity anomaly to represent an area of partial melt from which melts ascent vertically to the surface and northwards into the adjacent segment. By comparing tomographic studies of the Logachev and Southwest Indian Ridge Segment‐8 volcano we conclude that volcanic centers of ultraslow spreading ridges host spatially confined, circular partial melt areas below 10 km depth, in contrast to the shallow extended melt lenses along fast spreading ridges. Lateral feeding over distances of 35 km is possible at orthogonal spreading segments, but limited at the obliquely spreading Knipovich Ridge.

Rapid shifting of a deep magmatic source at Fagradalsfjall volcano, Iceland

Recent Icelandic rifting events have illuminated the roles of centralized crustal magma reservoirs and lateral magma transport1–4, important characteristics of mid-ocean ridge magmatism1,5. A consequence of such shallow crustal processing of magmas4,5 is the overprinting of signatures that trace the origin, evolution and transport of melts in the uppermost mantle and lowermost crust6,7. Here we present unique insights into processes occurring in this zone from integrated petrologic and geochemical studies of the 2021 Fagradalsfjall eruption on the Reykjanes Peninsula in Iceland. Geochemical analyses of basalts erupted during the first 50 days of the eruption, combined with associated gas emissions, reveal direct sourcing from a near-Moho magma storage zone. Geochemical proxies, which signify different mantle compositions and melting conditions, changed at a rate unparalleled for individual basaltic eruptions globally. Initially, the erupted lava was dominated by melts sourced from the shallowest mantle but over the following three weeks became increasingly dominated by magmas generated at a greater depth. This exceptionally rapid trend in erupted compositions provides an unprecedented temporal record of magma mixing that filters the mantle signal, consistent with processing in near-Moho melt lenses containing 107–108 m3 of basaltic magma. Exposing previously inaccessible parts of this key magma processing zone to near-real-time investigations provides new insights into the timescales and operational mode of basaltic magma systems.

Melt inclusion formation during olivine recrystallization: Evidence from stable isotopes

Melt inclusions are often used to infer melting processes or to determine source magmas that are usually overprinted in bulk rocks due to late stage mixing or near surface contamination. Here we present the first investigation of oxygen (O) isotope equilibrium between melt inclusions and their host olivines from arc samples. Olivines in all but one sample record either magma mixing or fractional crystallization. All six melt inclusions from Vulcano, 83% of seven from Sukumoyama, 44% of 21 from St Vincent, 37% of four from Iwate, and 21% of 13 from Aoba are not in isotopic equilibrium with their olivine host, despite the other major elements being in apparent equilibrium. A detailed study of some of the olivines shows that only a small volume around the melt inclusions is in equilibrium with its host. This strongly suggests that in these olivines melt inclusions are trapped in partly recrystallized olivines, highlighting the importance of magma mixing and crystal recycling in the magmatic plumbing system of these volcanoes. Oxygen isotope fractionation between melt inclusions and their host olivines, as well as phosphorus-δ18O systematics, could be used to better understand the formation of melt inclusions and crystal history. It could also provide valuable information to help characterize the magmatic plumbing system that the inclusions and their olivine hosts formed in (e.g., crystal rich-mush versus crystal poor melt lenses).

Long-lived dacitic magmatic systems and recharge dynamics in the Jemez Mountains volcanic field, western USA

We examine the connected history of dacite-dominant volcanic rocks of the Tschicoma Formation, erupted between 5.5 and 2 Ma from the Jemez Mountains volcanic field, western USA. Zircon samples from two separate eruptions have continuous SHRIMP U–Pb age spectra spanning 0.84–1.08 Myr duration (3.12–3.96 Ma and 3.50–4.58 Ma, respectively), following an episode of zircon crystallization 0.28–0.50 Myr earlier (at 4.46 Ma and 4.86 Ma, respectively). Zircon chemical variations, as well as ubiquitous resorption textures that commonly show large core-rim age differences (up to 720–740 kyr), suggest that they grew in separate melt lenses. Zircons were likely stored at near-solidus or even sub-solidus conditions after crystallization, but may have been reactivated in response to at least four major magma recharge events every 300–400 kyr and smaller events in between. A cycle of zircon dissolution (from heating), recrystallization (during cooling), and storage repeated in different locations in the Tschicoma mush system throughout its lifespan; each recharge-induced heating stage may last for several hundred to more than a thousand years based on calculations of zircon dissolution. We envisage the melt lenses to be distributed in a crystal mush zone, coalescing into a single magma batch as magma recharge occurs shortly before eruption. Once active, increasing magma supply rates may trigger large-scale partial melting of the pre-existing mush and caldera-forming eruptions.

A Reference Section Through Fast‐Spread Lower Oceanic Crust, Wadi Gideah, Samail Ophiolite (Sultanate of Oman): Whole Rock Geochemistry

AbstractDue to its inaccessibility, no complete and coherent data set exists for the composition of modern fast‐spreading oceanic crust. We sampled outcrops through 6,500 m of fossil oceanic crust in the Oman Ophiolite (Wadi Gideah Transect) that is regarded as best analogue of fast‐spreading crust on land. Here we report a complete set of whole‐rock major and trace element data displaying systematic and contrasting compositional trends in lower and upper gabbros being correlated with stratigraphic depth. A significant discontinuity in crystallization regime is observed at ∼3,525 m above the mantle‐crust boundary: gabbros below ∼3,525 m have in general very low incompatible element mass fractions which develop upwards in a barely noticeable way to more differentiated compositions while Mg# decreases. More pronounced trends indicating progressive fractionation of ascending melts can be observed for incompatible elements and their element ratios as a consequence of in situ crystallization. Locally, more variable compositions within narrow depth intervals testify for advanced differentiation in situ within individual sills. Gabbros above ∼3,525 m become significantly more evolved and show considerable variations in composition. Fractional crystallization and mixing processes in a transient axial melt lens control the composition of isotropic “varitextured” gabbros and sheeted dike basalts where fractionation of high field strength elemental ratios is minor. New average compositions of fast‐spread (paleo) oceanic crust are reported for major and 38 trace elements. Comparison with new data from Wadi Khafifah close to Wadi Gideah suggests robustness of crustal accretion processes in both space and time.

Melt in the Greenland EastGRIP ice core reveals Holocene warm events

Abstract. We present a record of melt events obtained from the East Greenland Ice Core Project (EastGRIP) ice core in central northeastern Greenland, covering the largest part of the Holocene. The data were acquired visually using an optical dark-field line scanner. We detect and describe melt layers and lenses, seen as bubble-free layers and lenses, throughout the ice above the bubble–clathrate transition. This transition is located at 1150 m depth in the EastGRIP ice core, corresponding to an age of 9720 years b2k. We define the brittle zone in the EastGRIP ice core as that from 650 to 950 m depth, where we count on average more than three core breaks per meter. We analyze melt layer thicknesses, correct for ice thinning, and account for missing layers due to core breaks. Our record of melt events shows a large, distinct peak around 1014 years b2k (986 CE) and a broad peak around 7000 years b2k, corresponding to the Holocene Climatic Optimum. In total, we can identify approximately 831 mm of melt (corrected for thinning) over the past 10 000 years. We find that the melt event from 986 CE is most likely a large rain event similar to that from 2012 CE, and that these two events are unprecedented throughout the Holocene. We also compare the most recent 2500 years to a tree ring composite and find an overlap between melt events and tree ring anomalies indicating warm summers. Considering the ice dynamics of the EastGRIP site resulting from the flow of the Northeast Greenland Ice Stream (NEGIS), we find that summer temperatures must have been at least 3 ± 0.6 ∘C warmer during the Early Holocene compared to today.

The Surface Deformation Signature of a Transcrustal, Crystal Mush‐Dominant Magma System

AbstractThe transcrustal, mush‐dominated magma storage paradigm, which posits liquid melt is heterogeneously distributed within a vertically extensive magma mush, differs significantly from classical geodetic models, where melt is stored within liquid‐dominant chambers within an elastic crust. Here, we present mechanical models consistent with transcrustal melt storage by separating the magmatic system into three domains: liquid melt lenses, surrounding crystal‐dominated poroelastic magma mush, and elastic crust. Our results indicate that pressure changes within the melt lens may induce surface displacements that approximate the displacements predicted by spheroidal pressure sources that mimic the geometry of the mush zone. Adopting constitutive parameters of the mush dependent on mush melt fraction, we show that a magma storage system will have an effective geometry inferred from surface displacements that smoothly transitions from the geometry of the melt lens to the geometry of the mush as mush melt fraction increases. This holds true across multiple storage zone geometries, including a “transcrustal” storage zone with a magma mush that extends deep in the crust. Accounting for the presence of a magma mush can lead to an increase in the estimated volume of injected or withdrawn magma (by several multiples) compared to values obtained using fully elastic models. Comparing erupted magma volumes to source volume changes allows for an estimation of magma compressibility; we show the presence of a mush can increase this estimated magma compressibility by up to approximately 50%, suggesting magmas may have higher bubble fraction than previous geodetically derived estimates.

Generation of Continental Alkalic Mafic Melts by Tholeiitic Melt–Mush Reactions: a New Perspective from Contrasting Mafic Cumulates and Dikes in Central Tibet

Abstract Alkalic basalts from continental regions have been widely used to fingerprint chemically and lithologically distinct components in the mantle. However, superimposed on mantle source effects are the effects of intra-crustal magma differentiation, which may be underestimated when exclusively considering fractional crystallization (FC) in melt-rich magma chambers. Here we present extensive mineralogical and geochemical data for mafic cumulates and dolerite dikes from central Tibet to demonstrate that tholeiitic magma differentiation occurred in crystal-rich mushes and melt-filled fractures. The mafic cumulates provide a snapshot of the percolation of water-poor (&amp;lt;0.9 wt% H2O), oxide-undersaturated, tholeiitic melts through slightly hotter mushes (50–60% crystals) with a framework of primitive olivine ± plagioclase crystals. The highly corroded olivines indicate that thermo-chemical disequilibrium between the later percolating melts and pre-existing crystal framework triggered melt–mush reactions in an open system. High-temperature (&amp;gt;1060°C) reactions were accompanied by a decreasing melt mass, causing the initial tholeiitic melts to evolve into alkalic mafic melts during slow cooling. This process involved the dissolution and reprecipitation of olivine, plagioclase, and clinopyroxene as well as the late crystallization of high-TiO2 (up to 3.4 wt%) clinopyroxene and kaersutite. In contrast, the same tholeiitic magmas (i.e. dolerite dikes) that intruded the cold country rocks along fractures produced silica-rich residual melts saturated in Fe–Ti oxides and sodic plagioclase by FC during rapid cooling. The residual melt was always in equilibrium with the instantaneous crystal during closed-system FC, so there was no reaction between the interstitial melt and adjacent crystal even after a high degree of crystallization. Alkalic mafic melts that can be formed by melt–mush reactions exhibit strong enrichment of incompatible element contents (e.g. Ti and alkali and rare earth elements) and ratios (e.g. Ce/Y), but little Si enrichment and Mg depletion, compared with the parental tholeiitic magma, which cannot be reproduced by simple FC. We argue that tholeiitic melt–mush reactions at mid- to lower-crustal pressures (≥0.5 GPa) are a new mechanism for the generation of interstitial alkalic mafic melts in continental mush reservoirs. Continental alkalic basalts could be produced by this mechanism if such interstitial melts were concentrated into melt lenses by compaction and then transported rapidly upward to the surface by focused flow.