Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “magma”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 685 records · Page 38

The Stop-and-Go Mechanism: Towards an Integrated Approach to Model Seismicity, Outgassing, Deformation, and Thermal Unrest at Active Volcanoes

Connecting the geophysical and geochemical signals recorded at and above the surface of volcanoes with source mechanisms is fundamental to understand transitions from quiescence to eruption, and to integrate the behavior of volcanoes with their regional seismotectonic and hydrological context. Some of these signals include shallow volcanic tremor, a long-lasting (minutes-to-years) ground vibration detected in volcanic areas during unrest; volcanic outgassing, which shows intriguing periodicities over multiple scales, as revealed by space-borne instruments and by ground-based, high-frequency (~1 Hz), sampling methods; inflation and deflation of volcanic edifices, commonly associated with subsurface volume changes; and large-scale (from a few to tens of km2), low-temperature, thermal anomalies, a newly-discovered signal that has been observed to emerge on volcanic flanks from months-to-years prior to gas- and magma-driven eruptions. Traditionally, these geophysical and geochemical signals have been modelled independently, which limits our interpretation of subsurface processes and thus our assessment of unrest. In this work, we show through lumped-parameter models, numerical simulations, and preliminary laboratory experiments that many of the signals recorded around volcanoes can be integrated into a common “stop-and-go” mechanism. In particular, we found that they can arise spontaneously when taking into account the feedbacks between the temporary accumulation of gas and/or magma in the crust (“stop”), and the diffusive transfer of gas and heat towards the surface (“go”). For example, we find that the “stop-and-go” mechanism can explain: (i) the emergence of monochromatic, broadband, and harmonic tremor; (ii) the periodic components commonly observed in volcanic outgassing time series; (iii) the link between large-scale thermal anomalies and shallow hydrothermal systems; (iv) the reported time lags between large-scale thermal anomalies and deformation; and (v) the breathing of potentially destructive volcanic calderas. The “stop-and-go” mechanism provides a consistent and realistic framework to link multiple geophysical and geochemical signals with the processes leading to volcanic unrest and eruption.

volcanoes↗

From Source to Surface: An Investigation of Magmatic Lunar Volatiles

Lunar missions led by NASA discovered water ice on the Moon’s surface [1-3], revolutionizing views of the abundance, distribution, and potential sources of H2O and other volatiles. Volatiles play an important role in affecting the rheological properties of minerals and melts and influencing magma eruption processes. In order to determine the indigenous volatile inventory of the Moon, it is vital to identify the magmatic and secondary processes that may have affected the volatile contents in lunar minerals [4]. As part of the Apollo Next Generation Sample Analysis (ANGSA) program we are investigating the petrogenesis of a set of four Apollo lunar basalts collected from the rim of Steno Crater 71035, 71037, 71055, and specially curated sample 71036. Our team is studying the petrology in 2D and 3D (this work), chronology (Pomeroy et al., LPSC 2022), oxidation state of sulfur in apatite (Brounce et al., LPSC 2022), and volatile inventory of the sample set in order to understand the age and genetic relationships among the basalts, their magma ascent conditions and eruptive signatures, and their degassing and surface histories. Here we present the first detailed study of the 2D and 3D mineralogy, textures, 3D vesiculation, and chemistry of these basalts to shed light on their magmatic, volcanic, and post-eruptive histories. Note: Extended abstract on document.

Z. E. Wilbur↗

Pink Spinel in Apollo Impact Melt Rock 68815: Implications for Mg-Suite Magmatism

Introduction: Magnesian rocks from the lunar highlands are collectively termed the Mg-suite. Characterized by high (>60) Mg# (molar 100×Mg/[Mg+Fe]) in mafic phases and calcic plagioclase, these rocks are plutonic to hypabyssal in origin, and include a range of bulk mineralogies such as troctolites, dunites, norites, gabbronorites, and spinel troctolites [1]. These Mg-suite lithologies have distinct trace element concentrations and ratios that differentiate them from other lunar rock types. These rocks are ancient, generally dated to between 4.5 and 4.1 Ga, although it is unknown if this represents the full range of Mg-suite ages [1–2]. Among the Mg-suite lithologies, the spinel troctolites are relatively rare, to date only found in polymict breccias [3]. Spinel troctolites, as their name suggests, consist of calcic plagioclase and forsteritic olivine, with minor amounts of spinel (MgAl2O4), ± pyroxene and cordierite [1,4]. This form of spinel is often called ‘pink’ spinel because of its appearance in thin section under plane polarized light (PPL; Fig. 1), due to minor amounts of Cr. Spinel troctolites are generally plutonic or hypabyssal in origin (subsequently exhumed and incorporated into polymict breccias), or formed through impact processes (e.g., crystalline impact melt) [5]. A spinel-rich lithology has also been found in the Moscoviense region of the Moon via the Moon Mineralogy Mapper (M3) and lacks other mafic phases [6]. Finally, while the Mg-suite sampled thus far consists of plutonic (or hypabyssal) rocks, the question remains if such magmas could have erupted on the surface of the Moon [7]. These magmas have much lower density than mare basalts, but little sample or remote sensing evidence has been found to support the idea that extrusive Mg-suite volcanism occurred [7]. Here, we present a coordinated microanalytical study of spinel-bearing lithic and mineral clasts found in Apollo sample 68815. These data will be used to understand their petrogenesis (magmatic or impact) and modification histories, and to shed light on the existence of volcanic Mg-suite rocks. Sample Description: Apollo sample 68815 is a polymict impact melt breccia containing a variety of lithic and mineral fragments embedded in devitrified impact melt. This sample was chipped off the top of a boulder at Station 8 during the Apollo 16 mission and had an original weight of nearly 1.8 kg. In this study, we investigated two polished thin sections of 68815: 68815,17 and 68815,148, both containing spinel. Methods: The thin sections of 68815 were studied using optical light microscopy (PPL, cross-polarized light, and reflected light) with a Keyence VHX-7100 Digital Microscope. Each section was then X-ray mapped for 13–14 elements using a Cameca SX100 electron probe microanalyzer (EPMA) located in the Kuiper Materials Imaging and Characterization Facility (KMICF) at the University of Arizona. We have obtained geochemical information about the phases (olivine, plagioclase, spinel, pyroxene) in the thin sections also using the EPMA. In addition, we have used ThermoScientific Helios NanoLab 660 Focused-Ion-Beam Scanning-Electron Microscope (FIB-SEM) and a Hitachi S-4800 SEM (both in KMICF) to obtain backscattered electron (BSE) images and energy dispersive Xray spectrometry (EDS) maps of areas of interest. Using a JEOL 7900F SEM at the Astromaterials Research & Exploration Science (ARES) at NASA Johnson Space Center (JSC), we have obtained electron backscatter diffraction (EBSD) maps of the spinel-bearing portions of the thin sections. The EBSD data were collected under beam conditions of 20 kV, and ~90 μA, with step sizes varying from 0.05 to 2 μm. Following EBSD data collection, we processed the data using AZtecCrystal and MTEX, a free MATLAB toolbox. Results: We have found clasts with subophitic textures, that consist of primarily olivine and plagioclase, with minor amounts of pink Mg-Al spinel and pyroxene (Fig. 1). These clasts are up to ~1 mm in length and contain spinels up to 50 μm across. We have additionally identified pink Mg-Al spinels within the impact melt (i.e., not contained in lithic clasts) in both thin sections. In one instance, a single spinel grain is approximately 300 μm across (Fig. 1b, 2). The spinel fragments embedded in impact melt have varying compositions, typically distinct from the compositions of spinels in the lithic clasts. Spinel-Bearing Clasts: Ten lithic clasts with similar textures and mineral compositions were identified between 68815,17 (two clasts) and ,148 (eight clasts). These clasts fall into two groups. The first has skeletal olivine with intergranular plagioclase, with minor amounts of pyroxene and spinel (Fig. 1a, 2c, 2d). The spinel in these clasts are found amid the plagioclase. The second group have an intergranular texture of olivine and plagioclase, again with minor spinel and pyroxene. The second group may contain spinels surrounded by plagioclase, and spinels enclosed in olivine. Spinels located within both clast types range from no apparent Cr zoning, to reverse zoning (Cr-enrichment inward; Fig. 2d), to normal zoning (Cr-enrichment outward). In the clasts thus far investigated with EPMA, plagioclase compositions range from An# (molar 100×Ca/[Ca+Na+K]) 92–96. Olivine Mg# ranged from 77 to 94, while pyroxene had Mg# from 54–84. Spinel in the clasts have Cr# (molar 100×Cr/[Cr+Al]) 2–4 and Mg# 88–91, which is within the range of pristine and plutonic spinel troctolites [8]. Isolated Spinels: These crystals are generally euhedral to subhedral, and can exhibit reverse Cr zoning (Cr enrichment inward) or no apparent Cr zoning. The spinels thus far investigated via EPMA have Cr# 9–14 and Mg# 65–82. The Cr# for these spinels is within the range reported by [8], but have lower Mg#. Future Work: We will continue to process the EBSD data for these lithic and mineral clasts. We will also continue to characterize these clasts using EPMA and SEM. By thoroughly characterizing the various spinels and spinel-bearing clasts, we aim to constrain the petrogenesis of these minerals and rock fragments. Acknowledgments: We thank NASA for the loan of these thin sections. We thank Ken Domanik and Jerry Chang for their support with data collection. Work was supported by a University of Arizona RII Core Facilities Pilot Program grant and start-up funds to JJB. We acknowledge support from NASA’s Planetary Science Research program for analysis performed at JSC. References: [1] Shearer C. K. et al. (2015) Am. Min. 100, 294–325. [2] Borg L. E. et al. (2020) GCA 290, 312–332. [3] Warren P. H. (1993) Am. Min. 78, 360–376. [4] Dymek R. F. et al. (1976) LPS VII, 2335–2378. [5] Treiman et al. (2019) Am. Min. 104, 370–384. [6] Pieters et al. (2011) JGR: Plan. 116:E00G08. [7] Prissel et al. (2016) Icarus 277, 319–329. [8] Prissel et al. (2016) Am. Min. 101, 1624–1635.

spinel↗

VESIcal: An Open-source Thermodynamic Model Engine for Mixed Volatile (H2O-CO2) Solubility in Silicate Melts

Modeling the solubility of volatiles in silicate melts is fundamental to the interpretation of volcanic systems and has implications for magma dynamics, eruption style, and material transport between the mantle, crust, and atmosphere. Recent advancements in computational capabilities and access to computing tools has outpaced the functionality and extensibility of previously available modeling platforms. Here we present VESIcal (Volatile Equilibria and Saturation Index calculator), the first comprehensive modeling tool for H2O, CO2, and mixed (H2O-CO2) solubility in silicate melts that: a) allows users access to seven popular models, with easy inter-comparison between models; b) provides universal functionality for all models (e.g., functions for calculating saturation pressures, degassing paths, etc.); c) can process large datasets (1,000’s of samples) automatically; d) can output computed data into an Excel spreadsheet or CSV file for post-modeling analysis; e) integrates plotting capabilities directly within the tool; and f) provides all of this within the framework of a python library, making the tool extensible by the user and allowing any of the model functions to be incorporated into any other code capable of calling python.Here we will provide a demonstration of VESIcal and its capabilities with applications to various volcanic processes affected by volatiles. VESIcal represents the first tool capable of directly comparing multiple solubility models and equations of state. We find that commonly used models predict surprisingly different volatile solubilities, particularly for pure CO2 or mixed CO2-H2O fluids. Even for melt compositions that are well represented in the calibration datasets of multiple models (e.g., MORBs), calculated solubilities for pure CO2 and pure H2O can deviate from one another by factors of >2 leading to 2x deviations in calculated saturation pressures (e.g., 5 to 10 kbar). The solubility of CO2 predicted by different rhyolitic models also differs substantially, overwhelming other sources of uncertainty such as analytical errors on measurements of volatile contents or uncertainties in crustal density profiles. This highlights the importance of model choice when drawing geological conclusions based on volatiles in magmas.

Kayla Iacovino↗

Comparing the Volatile Contents of Basaltic Rocks Through the Inner Solar System

Volatiles in planetary interiors play important roles in magma genesis, crust formation, eruption style, and even habitability of planetary surfaces. Allan Treiman’s career has touched on volatiles in the Moon, Mars, HED-parent body, the Earth, Venus, and other planetary bodies through studies of igneous rocks, volatile-bearing igneous minerals within these rocks, and the secondary alteration products produced from these rocks and minerals. His seminal work in many of these fields has inspired our careers and many others in our field. Here, we will build on and summarize what we know of the volatile content of basaltic rocks in the inner solar system. Specifically, this talk will focus on using volatile-bearing minerals (amphibole and apatite) along with bulk compositions to constrain the pre-eruptive, and potentially predegassed, volatile content of the magma and their source regions. We will then use these estimates to compare the volatile content of the Earth, Moon, Mars, Venus, and Mercury and what open questions there are for each planetary body.

Justin Filiberto↗

Terrestrial Aqueously Altered Magmatic Dike Forming Sulfate-Rich Hydrothermal Fluids to Constrain Martian Habitability

Intrusion of magma into a water-bearing crust can create a habitable hydrothermal environment, as seen on Earth. As the rock surrounding the intrusion is subjected to heat it can 1) release water, and different chemical species, from pre-existing minerals, and 2) interact with volatiles released from the cooling magma, which can create a geothermal brine with a distinct composition suitable for a range of chemotrophs. Hence, studying such environments on Earth is important in the search for life beyond our planet. The geothermal brine’s nature is dependent on the host rock and the magmatic properties. On Mars basaltic and sulfur-rich soils are dominant and therefore investigating basaltic magmatic interaction with sulfate-rich sediments and their associated habitats is key for recognizing potential habitats on Mars. The NASA Perseverance Rover at Jezero Crater has detected aqueously altered igneous rocks, potentially intrusive in origin, with secondary carbonates, sodium perchlorate and sulfates. Such an environment may have been habitable. Therefore, we investigate the intrusion of a magmatic dike into Jurassic sulfate-rich sediments from the San Rafael Swell, Utah, to reconstruct the geothermal fluids and assess their putative habitability.

B Baharier↗

Mantle Plume Magmatism in Elysium Planitia as Constrained by InSight Seismic Observations

Based on cratering statistics, there has been active volcanism in the last 100 Ma and possibly the last 1 to 10 Ma in both central Elysium Planitia and in the Cerberus Fossae region of south-eastern Elysium. Seismic measurements by the InSight mission show a strong concentration of seismicity in Cerberus Fossae which are interpreted as due to magma transport along dike systems. Together, these observations strongly support the existence of geologically recent and possibly of currently active magmatism in Elysium. InSight seismic observations were also used to construct seismic velocity models for the martian mantle. These models are effectively a regional aver-age of the seismic velocity structure for the Elysium region, where both the InSight lander and most of the measured seismic events are located. These seismic velocity models have been interpreted in terms of the lithospheric thickness and mantle potential temperature. Khan et al. found a lithospheric thickness of 400 to 600 km and a mantle potential temperature of 1325 to 1425 °C. Durán et al. found a lithospheric thickness of ~450 km and a mantle potential tempera-ture of 1375-1475 °C. Drilleau et al. found a thermal lithosphere of 420-660 km and a mantle potential temperature of 1380-1560 °C. In this study, we combine these InSight seismic observations with finite element mantle plume magma production models to place new constraints on the properties of the mantle beneath Elysium Planitia.

Walter S Kiefer↗

Gravitational Inferences on Lunar Igneous History Beyond Mare Volcanism

The source, distribution, and volume of igneous material on the Moon sheds light into its thermal evolution. Previous work has focused on extrusive mare volcanism, but intrusions and non-mare compositions are equally as important. Constraints on the ratio of magma production to volcanic output (intrusive/extrusive ratio) inform mantle melting rates and the conditions that allowed transport of magma through the crust. Additionally, the lithological diversity of lunar volcanic products sheds light on the timeline of lunar cooling and fracturing of the crust.

K Izquierdo↗

Timing, Abundance, and Spatial Extent of Initial Magmatism on the Moon Explained By Cumulate Mantle Overturn

We have recently quantified the timing, abundance, and spatial extent of lower mantle melting induced by cumulate overturn on the Moon through a series of 3D geodynamical models. Our dynamical modeling indicates that overturn of thin (~30-50 km) and weak ilmenite-bearing cumulates (IBC) triggers a rapid, short-lived, and widespread period of lower mantle melting which reproduces the key geochronological, volume, and spatial characteristics associated with the onset of secondary magmatism on the Moon (Figs. 1,2), and without energy contributions from KREEP (potassium, rare earth elements, phosphorus, radiogenic U, Th). Within the guiding paradigms of global differentiation via magma ocean crystallization and subsequent cumulate mantle overturn, our model provides explanation for near contemporaneous primary and secondary crust production constrained by geochronology of returned lunar samples and meteorites. In this abstract, we discuss our results in context with several intricacies of lunar chronology including models of a long-lived magma ocean, the hypothesis that mantle overturn was induced by the giant South Pole-Aitken basin forming impact, and ancient lunar zircon.

T. C. Prissel↗

Preliminary Characterization of Spinel Troctolite Clast Identified in Apollo Next Generation Sample Analysis (ANGSA) Core 73002

The lunar magnesian (Mg-) suite is a diverse lithological group composed of Mg-rich mafic minerals. The Mg-suite consists primarily of troctolites, norites, gabbronorites, and spinel troctolites. These rocks exhibit contradictory geochemical traits, with mafic minerals bearing high Mg# (molar 100×Mg/[Mg+Fe]) indicating primitive parental magmas, but also plagioclase saturation and enrichments in incompatible trace elements (i.e., KREEP) indicate more evolved parental magmas. A variety of models have been invoked to explain the occurrence and petrogenesis of the Mg-suite. Here, we report on a new spinel troctolite clast found in Apollo Next Generation Sample Analysis (ANGSA) program core 73002. Sample 73002 was acquired at Station 3 during Apollo 17 as the upper part of a double-drive tube with sample 73001. During processing, the spinel troctolite clast (0.066 g) was found at the 4.0–4.5 cm depth interval.

A. C. Stadermann↗

Quantifying mass flows at Mt. Cleveland, Alaska between 2001 and 2020 using satellite photogrammetry

Measuring eruption volume provides constraints on the magma supply rate and plumbing systems and therefore is a critical component for monitoring volcanoes. We use ArcticDEM—a large collection of time-dependent digital elevation models (DEMs) derived from satellite stereo-photogrammetry — to construct a first-of-its-kind measurement of the volumes of recent mass flows at Mount Cleveland, Alaska. We quantify the volume of the products of the 2001 eruption (the largest in the past half-century) as (54.8 ± 0.5) × 10 -6 m -3 covering a total area of 5.09 km -2 . The total volume of material loss at the summit crater is (0.67 ± 0.02) × 10 -6 m -3 , which is likely caused by later explosions and collapses of the shallow magma chamber. The total eruptive volume between 2017 and 2020 is (0.086 ± 0.002) × 10 -6 m -3 . Elevation changes associated with lahars are variable. On the upper northern slopes of the volcano, the lahar channels were almost exclusively erosive, suggesting that lahars originating at the summit eroded and entrained loose materials high on the cone. In general, lahar deposits were thickest near their toes, with some reaching ~20 m thickness.

Chunli Dai↗

Review of melt inclusions in lunar rocks: constraints on melt and mantle composition and magmatic processes

Mineral-hosted melt inclusions provide a window into magmatic processes and pre-eruptive liquid compositions. Because melt inclusions are small (typically < 100 µm), the study of lunar melt inclusions is enabled by advancements of microbeam instrumental techniques. In the 1970s immediately following the Apollo and Luna missions, major and minor oxide concentrations of lunar melt inclusions were measured using electron microprobes. The data were used to understand magma evolution, and they revealed the immiscibility of two silicate liquids in the late stage of lunar magma evolution. More recently, the development of secondary ion mass spectrometry as well as laser ablation–inductively coupled plasma–mass spectrometry has enabled the measurement of key volatile elements and other trace elements in lunar melt inclusions, down to about the 0.1 ppm level. The applications of these instruments have ushered in a new wave of lunar melt inclusion studies. Recent advances have gone hand in hand with improved understanding of post-entrapment loss of volatiles. These studies have provided deep insights into pre-eruptive volatiles in lunar basalts, the abundance of volatiles in the lunar mantle, the isotopic ratios of some volatile elements, and the partition of trace elements between host olivine and melt inclusions. The recent studies of lunar melt inclusions have played a critical role in establishing a new paradigm of a fairly wet Moon with about 100 ppm H2O in the bulk silicate Moon (rather than a “bone-dry” Moon) and have been instrumental in developing an improved understanding of the origin and evolution of the Moon.

Youxue Zhang↗

Kilauea Magnetotelluric Dataset

In 2002 and 2003 a collaborative effort was undertaken between Lawrence Berkeley National Laboratory, Sandia National Laboratories, the USGS Menlo Park, the USGS Hawaiian Volcano Observatory, and Electromagnetic Instruments Inc. to study the Kilauea volcano in Hawaii using the magnetotelluric (MT) technique. The work was motivated by a desire to improve understanding of the magma reservoirs and conduits within Kilauea and the East and Southwest Rift zones, which has implications for understanding Kilaueas plumbing system. An improved understanding of the rift zones has implications in understanding large-scale landslides that are generated in the Hilina Slump, which produce significant impacts on coastal communities. Up to eight stations operated simultaneously, with multiple remote reference sites, and data were processed using multi-station robust processing techniques. In total, data were acquired at 70 sites over the Southwest and East rift zones. Good to excellent quality data were obtained even in the harshest conditions, such as those encountered on the fresh lava flows of the East Rift Zone (ERZ), where electrical contact resistances are on the order of 100 kOhm. This data supports the continuing efforts to increase geothermal power on the island of Hawaii. Each of the 70 EDI files are the MT impedance tensors for 1 site. There is also a description of the processing of the data and a site map showing the locations of each site.

15 GEOTHERMAL ENERGY↗

Recreating Giants Impacts in the Laboratory: Shock Compression of MgSiO 3 Bridgmanite to 14 Mbar

Understanding giant impacts requires accurate description of how extreme pressures and temperatures affect the physical properties of the constituent materials. Here, we report shock experiments on two polymorphs of MgSiO 3 : enstatite and bridgmanite (perovskite) crystals.We obtain pressure-density shock equation of state to 14 Mbar and more than 9 g/cm 3 , a 40% increase in density from previous data on MgSiO 3 . Density-functional-theory molecular dynamics (DFT-MD) simulations provide predictions for the shock Hugoniot curves for bridgmanite and enstatite and suggest that the Grüneisen parameter decreases with increasing density. The good agreement between the simulations and the experimental data, including for the shock temperature along the enstatite Hugoniot reveals that DFT-MD simulations reproduce well the behavior of dense fluid MgSiO 3 .We also reveal a high optical reflectance indicative of a metal-like electrical conductivity which supports the hypothesis that magma oceans may contribute to planetary magnetic field generation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The Viscosity and Atomic Structure of Volatile-Bearing Melilititic Melts at High Pressure and Temperature and the Transport of Deep Carbon

Understanding the viscosity of mantle-derived magmas is needed to model their migration mechanisms and ascent rate from the source rock to the surface. High pressure–temperature experimental data are now available on the viscosity of synthetic melts, pure carbonatitic to carbonate–silicate compositions, anhydrous basalts, dacites and rhyolites. However, the viscosity of volatile-bearing melilititic melts, among the most plausible carriers of deep carbon, has not been investigated. In this study, we experimentally determined the viscosity of synthetic liquids with ~31 and ~39 wt% SiO 2 , 1.60 and 1.42 wt% CO 2 and 5.7 and 1 wt% H 2 O, respectively, at pressures from 1 to 4.7 GPa and temperatures between 1265 and 1755 °C, using the falling-sphere technique combined with in situ X-ray radiography. Our results show viscosities between 0.1044 and 2.1221 Pa·s, with a clear dependence on temperature and SiO 2 content. The atomic structure of both melt compositions was also determined at high pressure and temperature, using in situ multi-angle energy-dispersive X-ray diffraction supported by ex situ microFTIR and microRaman spectroscopic measurements. Our results yield evidence that the T–T and T–O (T = Si,Al) interatomic distances of ultrabasic melts are higher than those for basaltic melts known from similar recent studies. Based on our experimental data, melilititic melts are expected to migrate at a rate ~from 2 to 57 km·yr -1 in the present-day or the Archaean mantle, respectively.

36 MATERIALS SCIENCE↗

Core-Mantle Fractionation of Carbon in Earth and Mars: the Effects of Sulfur

Constraining carbon (C) fractionation between silicate magma ocean (MO) and core-forming alloy liquid during early differentiation is essential to understand the origin and early distribution of C between reservoirs such as the crust-atmosphere, mantle, and core of Earth and other terrestrial planets. Yet experimental data at high pressure (P)-temperature (T) on the effect of other light elements such as sulfur (S) in alloy liquid on alloy-silicate partitioning of C and C solubility in Fe-alloy compositions relevant for core formation is lacking. Here we have performed multi-anvil experiments at 6–13 GPa and 1800–2000 °C to examine the effects of S and Ni on the solubility limit of C in Fe-rich alloy liquid as well as partitioning behavior of C between alloy liquid and silicate melt (𝐷𝐷(sub C)(sup alloy/silicate)). The results show that C solubility in the alloy liquid as well as 𝐷𝐷(sub C)(sup alloy/silicate) decreases with increasing in S content in the alloy liquid. Empirical regression on C solubility in alloy liquid using our new experimental data and previous experiments demonstrates that C solubility significantly increases with increasing temperature, whereas unlike in S-poor or S-free alloy compositions, there is no discernible effect of Ni on C solubility in S-rich alloy liquid. Our modelling results confirm previous findings that in order to satisfy the C budget of BSE, the bulk Earth C undergoing alloy-silicate fractionation needs to be as high as those of CI-type carbonaceous chondrite, i.e., not leaving any room for volatility-induced loss of carbon during accretion. For Mars, on the other hand, an average single-stage core formation at relatively oxidized conditions (1.0 log unit below IW buffer) with 10-16 wt.% S in the core could yield a Martian mantle with a C budget similar to that of Earth’s BSE for a bulk C content of ~0.25-0.9 wt.%. For the scenario where C was delivered to the proto-Earth by a S-rich differentiated impactor at a later stage, our model calculations predict that bulk C content in the impactor can be as low as ~0.5 wt.% for an impactor mass that lies between 9-20% of present day Earth’s mass. This value is much higher than 0.05-0.1 wt.% bulk C in the impactor predicted by Li et al. (2016) because C-solubility limit of 0.3 wt.% in a S-rich alloy predicted by their models is significantly lower than the experimentally derived C-solubility of ~1.6 wt. % for the relevant S-content in the core of the impactor.

Carbon↗

Rapid transition from primary to secondary crust building on the Moon explained by mantle overturn

Geochronology indicates a rapid transition (tens of Myrs) from primary to secondary crust building on the Moon. The processes responsible for initiating secondary magmatism, however, remain in debate. Here we test the hypothesis that the earliest secondary crust (Mg-suite) formed as a direct consequence of density-driven mantle overturn, and advance 3D mantle convection models to quantify the resulting extent of lower mantle melting. Our modeling demonstrates that overturn of thin ilmenite-bearing cumulates ≤ 100 km triggers a rapid and short-lived episode of lower mantle melting which explains the key volume, geochronological, and spatial characteristics of early secondary crust building without contributions from other energy sources, namely KREEP (potassium, rare earth elements, phosphorus, radiogenic U, Th). Observations of globally distributed Mg-suite eliminate degree-1 overturn scenarios. We propose that gravitational instabilities in magma ocean cumulate piles are major driving forces for the onset of mantle convection and secondary crust building on differentiated bodies.

Moon↗

Seismic identification of geothermal prospecting in Harrat Rahat, Northern Arabian Shield

Harrat Rahat volcanic field lies along the Makkah-Madinah-Nafud active volcanic line. It has extensive linear weak zones of wrench faults that could provide ideal conduits for upwelling magma and hydrothermal fluids as well. Here, passive seismic data recorded by seismic stations installed within and around Harrat Rahat has been processed for the detection of microseismic events. Moreover, ambient noise was also gathered and analyzed to construct a subsurface velocity model. This model indicates low velocity zone in Harrat Rahat, and two areas with coincident basement faults and possibly with enhanced subsurface temperatures were identified. It is concluded that Harrat Rahat has high geothermal potential and promising for geothermal renewable power.

58 GEOSCIENCES↗