Expanding the application of the Eu-oxybarometer to the lherzolitic shergottites and nakhlites: Implications for the oxidation state heterogeneity of the Martian interior

Experimentally rehomogenized melt inclusions from the nakhlite Miller Range 03346 (MIL 03346) and the lherzolitic shergottite Allan Hills 77005 (ALH 77005) have been analyzed for their rare earth element (REE) concentrations in order to characterize the early melt compositions of these Martian meteorites and to calculate the oxygen fugacity conditions they crystallized under. D(Eu/Sm)pyroxene/melt values were measured at 0.77 and 1.05 for ALH 77005 and MIL 03346, respectively. These melts and their associated whole rock compositions have similar REE patterns, suggesting that whole rock REE values are representative of those of the early melts and can be used as input into the pyroxene Eu-oxybarometer for the nakhlites and lherzolitic shergottites. Crystallization fO_2 values of IW + 1.1 (ALH 77005) and IW + 3.2 (MIL 03346) were calculated. Whole rock data from other nakhlites and lherzolitic shergottites was input into the Eu-oxybarometer to determine their crystallization fO_2 values. The lherzolitic shergottites and nakhlites have fO_2 values that range from IW + 0.4 to 1.6 and from IW + 1.1 to 3.2, respectively. These values are consistent with some previously determined fO_2 estimates and expand the known range of fO_2 values of the Martian interior to four orders of magnitude. The origins of this range are not well constrained. Possible mechanisms for producing this spread in fO_2 values include mineral/melt fractionation, assimilation, shock effects, and magma ocean crystallization processes. Mineral/melt partitioning can result in changes in fO_2 from the start to the finish of crystallization of 2 orders of magnitude. In addition, crystallization of a Martian magma ocean with reasonable initial water content results in oxidized, water-rich, late-stage cumulates. Sampling of these oxidized cumulates or interactions between reduced melts and the oxidized material can potentially account for the range of fO_2 values observed in the Martian meteorites

Intra- and Intercrystalline Oxygen Isotope Variations in Minerals from Basalts and Peridotites

Igneous phenocrysts commonly exhibit zoning in major and trace element composition, reflecting (and potentially constraining) the differentiation and/or mixing histories of their parent melts. To date, little work has been done characterizing zonation of oxygen isotopes in minerals from mafic and ultramafic rocks. We present 259 ion probe (CAMECA ims-1280) measurements of δ^(18)O in 34 natural magmatic and mantle olivines and pyroxenes from five hand samples from diverse igneous environments. We compare δ^(18)O variations with zonation in other elements [especially P; analyzed by electron microprobe analysis (EMPA) and nano-secondary ionization mass spectrometry (nanoSIMS)]. There is generally a close (average within ~0·1–0·2 ‰) agreement between average δ^(18)O values of olivines measured by SIMS (standardized against San Carlos olivine) and independently known values for bulk separates from the same samples measured by laser fluorination. These data demonstrate that current ion microprobe techniques are not only precise but also accurate enough for study of sub-per-mil oxygen isotope variations in silicates (within ~0·2 ‰), provided samples are prepared and analyzed following strict guidelines. All but one of the 34 studied grains are homogeneous in δ^(18)O within a small multiple of analytical precision [estimated ±0·2‰, 1σ for most data; poorer for a subset of measurements made on small (~5 µm) spots]. This population of isotopically homogeneous grains includes some with oscillatory micrometer-scale P banding. The lack of δ^(18)O variations suggests that whatever factors lead to this common mode of trace element zonation have no detectable effect on melt–crystal partitioning of oxygen isotopes. Large (2‰) oxygen isotope variations are observed in one olivine glomerocryst from Mauna Kea, Hawaii. This glomerocryst contains P-rich domains that are either equant or skeletal or feathery in outline, and these P-rich domains are systematically low in δ^(18)O compared with adjacent, later-grown, P-poor olivine. This unusual oxygen isotope zonation pattern might reflect a kinetic fractionation during nucleation and growth of the cores of some olivine phenocrysts. We tested this hypothesis through measurements of δ^(18)O distributions in synthetic olivines grown at a range of rates and exhibiting diverse patterns of P zoning. These synthetic olivines are homogeneous in δ^(18)O, within the limits of our analyses (± 0·3–0·4‰ in this case) and show no connection between P zonation and oxygen isotope heterogeneity. We therefore think it more plausible that unusual O isotope zonation in the Mauna Kea glomerocryst reflects addition of a low-δ^(18)O component to some Hawaiian magmas just before nucleation of olivine. More generally, this study demonstrates the utility of modern SIMS techniques for in situ study of the subtle (~1‰ range) oxygen isotope variations characteristic of common mafic and ultramafic rocks

Identifying cryptotephra units using correlated rapid, nondestructive methods: VSWIR spectroscopy, X-ray fluorescence, and magnetic susceptibility

Understanding the frequency, magnitude, and nature of explosive volcanic eruptions is essential for hazard planning and risk mitigation. Terrestrial stratigraphic tephra records can be patchy and incomplete due to subsequent erosion and burial processes. In contrast, the marine sedimentary record commonly preserves a more complete historical record of volcanic activity as individual events are archived within continually accumulating background sediments. While larger tephra layers are often identifiable by changes in sediment color and/or texture, smaller fallout layers may also be present that are not visible to the naked eye. These cryptotephra are commonly more difficult to identify and often require time-consuming and destructive point counting, petrography, and microscopy work. Here we present several rapid, nondestructive, and quantitative core scanning methodologies (magnetic susceptibility, visible to shortwave infrared spectroscopy, and XRF core scanning) which, when combined, can be used to identify the presence of increased volcaniclastic components (interpreted to be cryptotephra) in the sedimentary record. We develop a new spectral parameter (BDI1000VIS) that exploits the absorption of the 1 µm near-infrared band in tephra. Using predetermined mixtures, BDI1000VIS can accurately identify tephra layers in concentrations >15–20%. When applied to the upper ∼270 kyr record of IODP core U1396C from the Caribbean Sea, and verified by traditional point counting, 29 potential cryptotephra layers were identified as originating from eruptions of the Lesser Antilles Volcanic Arc. Application of these methods in future coring endeavors can be used to minimize the need for physical disaggregation of valuable drill core material and allow for near-real-time recognition of tephra units, both visible and cryptotephra.Keywords: spectroscopy, tephrochronology, magnetic susceptibility, identification, cryptotephraKeywords: spectroscopy, tephrochronology, magnetic susceptibility, identification, cryptotephr

The relationship between eruptive activity, flank collapse, and sea level at volcanic islands: A long-term (>1 Ma) record offshore Montserrat, Lesser Antilles

Hole U1395B, drilled southeast of Montserrat during Integrated Ocean Drilling Program Expedition 340, provides a long (>1 Ma) and detailed record of eruptive and mass-wasting events (>130 discrete events). This record can be used to explore the temporal evolution in volcanic activity and landslides at an arc volcano. Analysis of tephra fall and volcaniclastic turbidite deposits in the drill cores reveals three heightened periods of volcanic activity on the island of Montserrat (?930 ka to ?900 ka, ?810 ka to ?760 ka, and ?190 ka to ?120 ka) that coincide with periods of increased volcano instability and mass-wasting. The youngest of these periods marks the peak in activity at the Soufrière Hills volcano. The largest flank collapse of this volcano (?130 ka) occurred towards the end of this period, and two younger landslides also occurred during a period of relatively elevated volcanism. These three landslides represent the only large (>0.3 km3) flank collapses of the Soufrière Hills edifice, and their timing also coincides with periods of rapid sea-level rise (>5 m/ka). Available age data from other island arc volcanoes suggests a general correlation between the timing of large landslides and periods of rapid sea-level rise, but this is not observed for volcanoes in intra-plate ocean settings. We thus infer that rapid sea-level rise may modulate the timing of collapse at island arc volcanoes, but not in larger ocean-island settings

Heat flow in the Lesser Antilles island arc and adjacent back arc Grenada basin

Using temperature gradients measured in 10 holes at 6 sites, we generate the first high fidelity heat flow measurements from Integrated Ocean Drilling Program drill holes across the northern and central Lesser Antilles arc and back arc Grenada basin. The implied heat flow, after correcting for bathymetry and sedimentation effects, ranges from about 0.1 W/m² on the crest of the arc, midway between the volcanic islands of Montserrat and Guadeloupe, to 15 km from the crest in the back arc direction. Combined with previous measurements, we find that the magnitude and spatial pattern of heat flow are similar to those at continental arcs. The heat flow in the Grenada basin to the west of the active arc is 0.06 W/m², a factor of 2 lower than that found in the previous and most recent study. There is no thermal evidence for significant shallow fluid advection at any of these sites. Present-day volcanism is confined to the region with the highest heat flow.American Geophysical Union – Geochemistry, Geophysics, Geosystems. This is the publisher’s final pdf. The published article is copyrighted by the American Geophysical Union and can be found at: http://www.agu.org/journals/gc/.Keywords: volcanic arc., Lesser Antilles, IODP, heat flow, back arc, Grenada basinKeywords: volcanic arc., Lesser Antilles, IODP, heat flow, back arc, Grenada basi

Permeability and pressure measurements in Lesser Antilles submarine slides: Evidence for pressure-driven slow-slip failure

Recent studies hypothesize that some submarine slides fail via pressure-driven slow-slip deformation. To test this hypothesis, this study derives pore pressures in failed and adjacent unfailed deep marine sediments by integrating rock physics models, physical property measurements on recovered sediment core, and wireline logs. Two drill sites (U1394 and U1399) drilled through interpreted slide debris; a third (U1395) drilled into normal marine sediment. Near-hydrostatic fluid pressure exists in sediments at site U1395. In contrast, results at both sites U1394 and U1399 indicate elevated pore fluid pressures in some sediment. We suggest that high pore pressure at the base of a submarine slide deposit at site U1394 results from slide shearing. High pore pressure exists throughout much of site U1399, and Mohr circle analysis suggests that only slight changes in the stress regime will trigger motion. Consolidation tests and permeability measurements indicate moderately low (~10⁻¹⁶–10⁻¹⁷ m²) permeability and overconsolidation in fine-grained slide debris, implying that these sediments act as seals. Three mechanisms, in isolation or in combination, may produce the observed elevated pore fluid pressures at site U1399: (1) rapid sedimentation, (2) lateral fluid flow, and (3) shearing that causes sediments to contract, increasing pore pressure. Our preferred hypothesis is this third mechanism because it explains both elevated fluid pressure and sediment overconsolidation without requiring high sedimentation rates. Our combined analysis of subsurface pore pressures, drilling data, and regional seismic images indicates that slope failure offshore Martinique is perhaps an ongoing, creep-like process where small stress changes trigger motion

Coring disturbances in IODP piston cores with implications for offshore record of volcanic events and the Missoula megafloods

Piston cores collected from IODP drilling platforms (and its predecessors) provide the best long-term geological and climatic record of marine sediments worldwide. Coring disturbances affecting the original sediment texture have been recognized since the early days of coring, and include deformation resulting from shear of sediment against the core barrel, basal flow-in due to partial stroke, loss of stratigraphy, fall-in, sediment loss through core catchers, and structures formed during core recovery and on-deck transport. The most severe disturbances occur in non-cohesive (sandy) facies, which are particularly common in volcanogenic environments and submarine fans. Although all of these types of coring disturbances have been recognized previously, our contribution is novel because it provides an easily accessible summary of methods for their identification. This contribution gives two specific examples on the importance of these coring disturbances. We show how suck-in of sediments during coring artificially created very thick volcaniclastic sand layers in cores offshore Montserrat and Martinique (Lesser Antilles). We then analyze very thick, structureless sand layers from the Escanaba Trough inferred to be a record of the Missoula mega-floods. These sand layers tend to coincide with the base of core sections, and their facies suggest coring disturbance by basal flow-in, destroying the original structure and texture of the beds. We conclude by outlining and supporting IODP-led initiatives to further reduce and identify coring disturbances, and acknowledge their recent successes in drilling challenging sand-rich settings, such as during IODP Expedition 340

Compositional Characterization of Glassy Volcanic Material From VNIR and MIR Spectra Using Partial Least Squares Regression Models

<p>This is supporting data for the paper titled "Compositional Characterization of Glassy Volcanic Material From VNIR and MIR Spectra Using Partial Least Squares Regression Models" by Leight et al. (submitted to JGR-P 11/23). Table S1 lists each spectrum used to train PLS models, its source, and which training datasets the spectrum was included in. Zip files contain the MIR and VNIR PLS model files. Model files are .asc, and can be run using the code at Ytsma, (2022), https://doi.org/10.5281/zenodo.7347345. </p&gt