ALBERT

Description: Abstract Tectonic extension of continental lithosphere creates accommodation space in which sediments are deposited. Climate‐driven processes provide the mechanism by which mass is detached from hillslopes and sediments are transported into this accommodation space. These two forcings, climate and tectonics, act together to create either endorheic (internally drained) or exorheic (externally drained) rift basins. Here we use a large‐scale dynamic landscape evolution‐tectonics model to understand the contribution of tectonic processes in endorheic‐exorheic transitions. In the model, extension results in opening of an asymmetric half‐graben along a listric normal fault. Rift opening occurs in the models in wet, temperate, or semi‐arid climates where runoff and evapotranspiration are varied. Our numerical experiments show that slow rift‐opening rates, a slowing‐down of rift opening, or increase of headwater topography (e.g., upstream epeirogenic uplift), are tectonic situations that can cause a transition from an endorheic to an exorheic drainage state in a rift basin. Our results also show that wet climate conditions lead to a permanent exorheism that persists regardless of rift opening rates. In semi‐arid climates, endorheic conditions are favored, and may last for the duration of rifting except for when rift opening is very slow. These results form an interpretive framework to study endorheic and exorheic drainage systems in natural continental rifts. In the slow‐opening Rio Grande rift, the endorheic‐exorheic transition may have occurred without dramatic climate changes. Lake‐level variations in East African rift basins are predicted by our models to result from variations in climate.

Description: Abstract Crustal extension is commonly thought to be accommodated by faults that strike orthogonal and obliquely to the regional trend of the minimum compressive stress (σ3). Activation of oblique faults can, however, be conceptually problematic as under Andersonian faulting, it requires preexisting crustal weaknesses, high fluid pressures, and/or stress rotations. Furthermore, measurements of incremental fault displacements, which are typically used to identify oblique faulting, do not necessarily reflect regional stresses. Here, we assess oblique faulting by calculating the stress ratio (σ3/σ1, where σ1 is the maximum compressive stress), slip tendency, and effective coefficient of friction (μs′) required to reactivate variably striking normal faults under different trends of σ3. We apply this analysis to NW and NNE striking active faults at the southern end of the Malawi Rift, where NE‐SW, ENE‐WSW, E‐W, and SE‐NW σ3 trends have previously been proposed. A uniform σ3 trend is inferred for this region as recent joints sets do not rotate along the rift. With a NE‐SW trending σ3, NW‐striking faults are well oriented, however, NNE‐striking faults require μs′ 〈 0.6 to reactivate. This is inconsistent with a lack of frictionally weak phyllosilicates detected in the fault zone rocks. With an ENE‐WSW to E‐W trending σ3, all faults can reactivate at μs′ 〉 0.55. These σ3 trends are also comparable to a focal mechanism stress inversion, regional joint orientations, and previously reported geodetically derived extension directions. We therefore conclude that unlike typical models of oblique rifting, the southern Malawi Rift consists of faults that all strike slightly oblique to σ3.

Description: Abstract Ahyi seamount, a shallow submarine volcano in the Northern Mariana Islands, began erupting on 23 April 2014. Hydroacoustic eruption signals were observed on the regional Mariana seismic network and on distant hydrophones, and National Oceanic and Atmospheric Administration (NOAA) scuba divers working in the area soon after the eruption began heard and felt underwater explosion sounds. The NOAA crew observed yellow‐orange bubble mats along the shore of neighboring Farallon de Pájaros Island, but no other surface manifestations of the eruption were reported by the crew or observed in satellite data. Here, we detail the eruption chronology and its morphologic impacts through analysis of seismic and hydroacoustic recordings and repeat bathymetric mapping. Throughout the 2‐week‐long eruption, Ahyi produced several thousand short, impulsive hydroacoustic signals that we interpret as underwater explosions as well as tremor near the beginning and end of the sequence. The initial tremor, which occurred for 2 hr, is interpreted as small phreatomagmatic explosions. This tremor was followed by a 90‐min pause before the characteristic impulsive signals began. Occasional tremor (lasting up to a few minutes) during the last 1.5 days of the eruption is interpreted as more sustained eruptive activity. Bathymetric changes show that a new crater, about 150 m deep, formed near the former summit and a large landslide chute formed on the southeastern flank. Comparing to other geophysically detected submarine eruptions, we find that the signals from the 2014 Ahyi eruption were more similar to those from other shallow or at‐surface submarine eruptions than those at deep (〉500 m) eruptions.

Description: Abstract Seismic anisotropy records past and present tectonic deformations and provides important constraints for understanding the structure and dynamics of the Earth's interior. In this work, we use tremendous amounts of high‐quality P wave arrival times from local and regional earthquakes to determine a high‐resolution tomographic model of 3‐D P wave azimuthal anisotropy down to 1,000‐km depth beneath East Asia. Our results show that trench‐parallel fast‐velocity directions (FVDs) are visible in the shallow portion of the subducting Pacific slab (〈80 km), whereas the deeper portion of the Pacific slab mainly exhibits trench‐normal FVDs, except for the stagnant slab in the mantle transition zone (MTZ) where obvious NE‐SW FVDs are revealed. The FVDs in the subslab mantle change from a subduction‐parallel trend at depths of 80–400 km to a subduction‐normal trend in the MTZ. Large‐scale low‐velocity anomalies are revealed beneath the Philippine Sea plate where the FVD is NE‐SW. The FVDs along the Izu‐Bonin arc and in a slab gap exhibit a striking anticlockwise toroidal trend. All these features may reflect complex 3‐D flows in the mantle wedge due to tearing and dehydration processes of the subducting Pacific slab. The subducting Pacific slab is split at ~300‐km depth under the Bonin arc and then penetrates into the lower mantle, whereas under East Asia the Pacific slab becomes stagnant in the MTZ and reaches the North‐South Gravity Lineament in China. The intraplate volcanoes in East Asia are caused by hot and wet upwelling flows in the big mantle wedge above the stagnant Pacific slab.

Description: Abstract We determine mass transport and structural properties of binary liquid iron alloys over a wide density (5.055–11.735 g·cm−3) and temperature range (2,500–6,500 K) using first‐principles molecular dynamics. Compositions consist of 96 at% Fe and 4 at% ϕ, where ϕ = H, C, N, O, Mg, Si, S, or Ni. Self‐diffusion coefficients (D) of Fe and ϕ range from 3.5·10−9 to 1.9·10−7 m2·s−1. Results show a relation between mean atomic radius and diffusivity ratio for the alloying element and iron: Si and Ni are “iron‐like” with similar atomic radii and D compared with those of Fe; H, C, N, O, and S are “small non‐iron‐like” with smaller atomic radii and larger D; and Mg transitions from “large non‐iron‐like” with a larger atomic radius and smaller D at low density to iron‐like under conditions of the Earth's core. The effect of pressure on D for C, N, and O is negligible for densities below ~8 g·cm−3, accompanied by an increase in average coordination numbers to ~6, and an increase in mean atomic radii. For densities above ~8 g·cm−3, diffusivities and atomic radii of these elements decrease monotonically with pressure, which is typical for the iron‐like alloying elements as well as for H, Mg, and S over the whole compression range. While atomic radius ratios move toward unity with compression, diffusivity ratios for the alloying element relative to iron tend to increase for the “non‐iron‐like” elements with density.

Description: Abstract Grain size is one of the most fundamental properties of sediments. It is frequently used in paleoclimate, paleoceanographic and paleoenvironmental research as a proxy for river discharge, current and wind strength, and to identify mass flow deposits. Measuring grain‐size is, however, time‐consuming and destructive. Given the strong influence of grain size on sediment inorganic geochemistry, single element variations measured by e.g., X‐ray fluorescence (XRF) core scanning are increasingly used to estimate grain‐size variations at high resolution in sediment cores. This approach is however limited to a narrow grain‐size range since individual elements only monotonically relate to grain size over a narrow size range. Here, we present a simple, code‐free, multi‐element method based on Partial Least Square regression to predict sediment mean grain size from inorganic geochemical data over the range of sizes commonly encountered in sedimentary basins (clay to sand). The method was first tested using river sediment samples separated in eleven grain‐size fractions and it was later successfully applied to two sediment cores from the Chilean fjords. Our method only requires measuring grain size on a limited number (around ten) of selected training samples, and it allows to predict mean grain size at XRF core scanner resolution. This method has the potential to be applied to any lake or marine sediment core, provided sediment provenance, weathering, and diagenesis remain relatively stable through time, and we anticipate that it will result in a significant increase in the resolution of sediment proxy records of climate and environmental change.

Description: Abstract Dunite, harzburgite, and clinopyroxenite xenoliths from Kharchinsky volcano, Kamchatka have abundances and ratios of incompatible trace elements similar to those in arc volcanic rocks (elevated Ba/Th, La/Yb, Nd/Hf, Sr/Y). All orthopyroxenes and some clinopyroxenes in the peridotites have U‐shaped rare‐earth element patterns. Negative Ce anomalies are present in orthopyroxenes with Ce/Ce* as low as 0.01 and down to 0.22 in whole‐rock peridotite data. Ce anomaly growth is linked to increasing La/Sm and enrichments in Rb, U, Pb, and Ba over La and Ce. Isotopes (Pb, Sr, Nd, Hf) indicate pelagic sediment and hydrothermal crusts play no role in Ce anomaly development. Instead, Ce anomalies appear to be products of fluid transport and elemental scavenging under oxidizing conditions beneath the deep forearc. Textures and compositions of aluminous green spinels indicate most of the peridotites were partially melted and recrystallized at depth. Veins and pockets of amphibole reflect impregnation late in the petrogenesis of the rocks by melts similar to Kamchatka basalts. Orthopyroxenite xenoliths are fine‐grained with fibrous orthopyroxene that has high‐Mg/Mg+Fe (up to 0.96) and generally lower CaO and Al2O3 compared to peridotite orthopyroxenes, and perhaps formed by reaction of siliceous fluids with olivine. Kharchinsky xenoliths have Pb, Sr and Nd isotopes similar to Kamchatka volcanic rocks, but Hf isotopes in clinopyroxenites and gabbros are more radiogenic by 1‐3 epsilon units. Patterns in isotopic data indicate a compositional change in the source of Kamchatka volcanism within the past 20 million years.

Description: Abstract A major goal in Earth Science has been to understand how geochemical characteristics of lavas at the Earth's surface relate to the location and formation history of specific regions in the Earth's interior. For example, some of the strongest evidence for the preservation of primitive material comes from low 4He/3He ratios in ocean island basalts, but the location of the primitive helium reservoir(s) remains unknown. Here we combine whole‐mantle seismic tomography, simulations of mantle flow, and a global compilation of new and existing measurements of the 4He/3He ratios in ocean island basalts to constrain the source location of primitive 4He/3He material. Our geodynamic simulations predict the present‐day surface expression of plumes to be laterally offset from their lower mantle source locations. When this lateral offset is accounted for, a strong relationship emerges between minimum 4He/3He ratios in oceanic basalts and seismically slow regions, which are generally located within the two large low shear‐wave velocity provinces (LLSVPs). Conversely, no significant relationship is observed between maximum 208Pb*/206Pb* ratios and seismically slow regions in the lowermost mantle. These results indicate that primitive materials are geographically restricted to LLSVPs, while recycled materials are more broadly distributed across the lower mantle. The primitive nature of the LLSVPs indicates these regions are not composed entirely of recycled slabs, while complementary xenon and tungsten isotopic anomalies require the primitive portion of the LLSVPs to have formed during Earth's accretion, survived the Moon‐forming giant impact, and remained relatively unmixed during the subsequent 4.5 billion years of mantle convection.

Description: Abstract A cave monitoring study in Hatchet Bay Cave on the island of Eleuthera, Bahamas, has examined the origins of variations in oxygen and carbon isotopic and minor element composition in cave calcites. Every 3 to 8 months, between 2012 and 2016, temperature, humidity, cave air (δ13CCO2), dripwaters (δ18O and δ2H values, and Ca, Sr, and Mg concentrations), and the chemical composition of precipitating calcite (δ18O and δ13C values, and Ca, Sr, and Mg concentrations) were analyzed in two rooms in the cave. Results from the elemental analyses show that throughout the cave prior calcite precipitation (PCP) was a driver of the elemental chemistry of the precipitated calcites. In addition, cave calcites show that δ13C and δ18O values were positively correlated with Mg/Ca ratios. The Mg/Ca ratios were also positively correlated with lower calcite precipitation rates. Therefore, water/rock interactions may also influence δ13C and δ18O values and Mg/Ca ratios of the calcite. Differences were observed between the two rooms, with the Main Room of the cave exhibiting increased PCP, more ventilation, lower calcite precipitation rates, and δ18O values which were farther from equilibrium when compared to the more isolated portion of the cave. These results also validated previous interpretations from Pleistocene stalagmites collected from a nearby Bahamian cave suggesting that a positive covariation between Mg/Ca and δ13C values reflects water/rock interactions.

Description: Abstract The ferrimagnetic properties of soils are used to quantitatively reconstruct paleomonsoon precipitation from Chinese loess. Numerous magneto‐climofunctions have been established based on the magnetic proxies that are selectively sensitive to neoformation of fine‐grained superparamagnetic (SP) or single‐domain (SD) ferrimagnetic particles. Accumulating evidence has indicated that maghemite is the final product of the ferrimagnetic phases during pedogenesis in loessic soils. Quantitative estimates of abundance of maghemite of both SP and SD grains is therefore still required in developing magneto‐climofunctions. Here, we present detailed measurements on a suite of modern soil samples from the Chinese Loess Plateau to determine pedogenic ferrimagnetic mineralogy and to develop a new magneto‐climofunction based on a new parameter derived from the high‐temperature‐dependent magnetic susceptibility. Particle‐size fractionation processes combined with magnetic measurements indicate that fine‐grained SP and SD maghemite is the dominant pedogenic ferrimagnetic phases. High‐temperature dependent susceptibility measurements show that the thermally‐induced susceptibility drops between ~230°C and ~400°C during heating mainly result from the conversion of maghemite to hematite. We proposed a new parameter quantifying changes in the temperature dependence of magnetic susceptibility between 230‐400°C, "χtd ", that captures the concentration of pedogenically formed maghemite. Results show that χtd has a strong correlation with known quantities of maghemite in synthetic standard samples, and that χtd of modern soils correlates with modern mean annual precipitation (MAP) quite well (R2=0.82, n=24). The established χtd‐MAP climofunction provides a new approach to reconstructing paleorainfall during past warm interglacials from paleosols in Chinese loess.

Description: Abstract Statistical analysis of geomagnetic paleosecular variation (PSV) and time averaged field (TAF) has been largely based on global compilations of paleomagnetic data from lava flows. These show different trends in the averaged inclination anomaly (ΔI) between the two hemispheres, with small positive (〈2°) anomalies in mid‐southern latitudes and large negative (〉‐5°) anomalies in mid‐northern latitudes. To inspect the large ΔI between 20°N‐40°N we augment the global data with a new paleomagnetic dataset from the Golan‐Heights (GH), a Plio‐Pleistocene volcanic plateau in northeast Israel, located at 32°N‐33°N. The GH dataset consists of 91 lava flows sites: 40 sites obtained in the 1990s and 51 obtained in this study. The chronology of the flows is constrained by 57 40Ar/39Ar ages: 39 from previous studies and 18 from this study, which together cover most of the GH plateau. We show that the 1990s dataset might be affected by block rotations and does not fully sample PSV. The Plio‐Pleistocene pole (86.3°N, 120.8°E, N=44, k=25, α95=4.4°), calculated after applying selection criteria with Fisher precision parameter (k) ≥ 100 and number of specimens per site (n) ≥ 5 is consistent with a geocentric axial dipole field and shows smaller inclination anomaly (ΔI=‐0.4°) than predicted by global compilations and PSV models. Re‐examination of the inclination anomaly in the global compilation using different calculation methods and selection criteria suggests that inclination anomaly values are affected by: (1) inclusion of poor quality data, (2) averaging data by latitude bins and (3) the way the inclination anomaly is calculated.

Description: No abstract is available for this article.

Description: Abstract We propose an explanation to the enigmatic synrift erosional unconformities reported along the distal domain of several magma‐poor rifted margins. Using thermomechanical numerical modeling, we show that transient emersion of (future) distal domains following a phase of subsidence can be explained by asynchronous necking of first the upper mantle and subsequently the crust, without the need of prominent normal faulting caused by strain softening, mantle phase transitions, or magmatic processes. When the upper crust and upper mantle are mechanically decoupled by a weak lower crust and, in the absence of any prominent rheological heterogeneity, upper mantle, necking starts first because of the higher deviatoric stresses associated with its larger effective viscosity. Consequently, the ductile lower crustal material flows toward the necked mantle domain, delaying thinning of the overlying crust. Once the necked lithospheric mantle has locally lost most of its strength, the overdeepened Moho moves upward toward an isostatically compensated depth. This flexural rebound causes uplift and emersion of distal parts of the rift system that are composed of still relatively thick crust and triggers the necking of the overlying crust. Early necking of the upper mantle causes a transient heating event with temperatures up to 750 °C at the base of the crust in the (future) distal domain. The onset of this thermal event slightly predates emersion of the (future) distal domain. These results are consistent with field observations and thermochronological data from the fossil Alpine Tethys margins, as well as with seismic observations from several present‐day rifted margins.

Description: Abstract Observations of seismic anisotropy can provide direct constraints on the character of mantle flow in subduction zones, critical for our broader understanding of subduction dynamics. Here we present over 750 new SKS splitting measurements in the vicinity of Mount St. Helens in the Cascadia subduction zone using a combination of stations from the iMUSH broadband array and Cascades Volcano Observatory network. This provides the highest density of splitting measurements yet available in Cascadia, acting as a focused “telescope” for seismic anisotropy in the subduction zone. We retrieve spatially consistent splitting parameters (mean fast direction Φ: 74°, mean delay time ∂t: 1.0 s) with the azimuthal occurrence of nulls in agreement with the fast direction of splitting. When averaged across the array, a 90° periodicity in splitting parameters as a function of back azimuth is revealed, which has not been recovered previously with single‐station observations. The periodicity is characterized by a sawtooth pattern in Φ with a clearly defined 45° trend. We present new equations that reproduce this behavior based upon known systematic errors when calculating shear wave splitting from data with realistic seismic noise. The corrected results suggest a single layer of anisotropy with an ENE‐WSW fast axis parallel to the motion of the subducting Juan de Fuca plate; in agreement with predictions for entrained subslab mantle flow. The splitting pattern is consistent with that seen throughout Cascadia, suggesting that entrainment of the underlying asthenosphere with the subducting slab is coherent and widespread.

Description: Abstract Recent seismic studies indicate the presence of seismic anisotropy near subducted slabs in the transition zone and uppermost lower mantle (mid‐mantle). In this study, we investigate the origin of radial anisotropy in the mid‐mantle using 3‐D geodynamic subduction models combined with mantle fabric simulations. These calculations are compared with seismic tomography images to constrain the range of possible causes of the observed anisotropy. We consider three subduction scenarios: (i) slab stagnation at the bottom of the transition zone; (ii) slab trapped in the uppermost lower mantle; and (iii) slab penetration into the deep lower mantle. For each scenario, we consider a range of parameters, including several slip systems of bridgmanite and its grain‐boundary mobility. Modeling of lattice‐preferred orientation shows that the upper transition zone is characterized by fast‐SV radial anisotropy anomalies up to −1.5%. For the stagnating and trapped slab scenarios, the uppermost lower mantle is characterized by two fast‐SH radial anisotropy anomalies of ∼+2% beneath the slab's tip and hinge. On the other hand, the penetrating slab is associated with fast‐SH radial anisotropy anomalies of up to ∼+1.3% down to a depth of 2,000 km. Four possible easy slip systems of bridgmanite lead to a good consistency between the mantle modeling and the seismic tomography images: [100](010), [010](100), [001](100), and . The anisotropy anomalies obtained from shape‐preferred orientation calculations do not fit seismic tomography images in the mid‐mantle as well as lattice‐preferred orientation calculations, especially for slabs penetrating into the deep lower mantle.

Description: Abstract Building on an earlier study that confirmed the stability of the 405‐kyr eccentricity climate cycle and the timing of the Newark‐Hartford astrochronostratigraphic polarity time scale (N‐H APTS) back to 215 Ma, we extend the magnetochronology of the Late Triassic Chinle Formation to its basal unconformity in scientific drill core PFNP‐1A from Petrified Forest National Park (Arizona, USA). The 335‐m‐thick Chinle section is imprinted with paleomagnetic polarity zones PF1r to PF10n, which we correlate to chrons E17r to E9n (~209 to 224 Ma) of the N‐H APTS. A sediment accumulation rate of ~34 m/Myr can be extended down to ~270 meters, close to the base of the Sonsela Member and the base of magnetozone PF5n, which we correlate to chron E14n that onsets at 216.16 Ma. Magnetozones PF5r to PF10n in the underlying 65‐m‐thick section of the mudstone‐dominated Blue Mesa and Mesa Redondo members plausibly correlate to chrons E13r to E9n, indicating a sediment accumulation rate of only ~10 m/Myr. Published high precision U‐Pb detrital zircon dates from the lower Chinle tend to be several million years older than the magnetochronological age model. The source of this discrepancy is unclear but may be due to sporadic introduction of juvenile zircons that get recycled. The new magnetochronological constraint on the base of the Sonsela Member brings the apparent timing of the included Adamanian‐Revueltian land vertebrate faunal zone boundary and the Zone II to Zone III palynofloral transition closer to the temporal range of the ~215 Ma Manicouagan impact structure in Canada.

Description: Abstract Using 8‐25s period Rayleigh and Love wave phase velocity dispersion data extracted from seismic ambient noise, we (i) model the 3D shear wave velocity structure of the West Antarctic crust and (ii) map variations in crustal radial anisotropy. Enhanced regional resolution is offered by the UK Antarctic Seismic Network. In the West Antarctic Rift System (WARS), a ridge of crust ~26‐30km thick extending south from Marie Byrd Land separates domains of more extended crust (~22km thick) in the Ross and Amundsen Sea Embayments, suggesting along‐strike variability in the Cenozoic evolution of the WARS. The southern margin of the WARS is defined along the southern Transantarctic Mountains (TAM) and Haag Nunataks‐Ellsworth Whitmore Mountains (HEW) block by a sharp crustal thickness gradient. Crust ~35‐40km is modelled beneath the Haag Nunataks‐Ellsworth Mountains, decreasing to ~30‐32km km thick beneath the Whitmore Mountains, reflecting distinct structural domains within the composite HEW block. Our analysis suggests that the lower crust and potentially the mid crust is positively radially anisotropic (VSH 〉 VSV) across West Antarctica. The strongest anisotropic signature is observed in the HEW block, emphasising its unique provenance amongst West Antarctica's crustal units, and conceivably reflects a ~13km thick metasedimentary succession atop Precambrian metamorphic basement. Positive radial anisotropy in the WARS crust is consistent with observations in extensional settings, and likely reflects the lattice‐preferred orientation of minerals such as mica and amphibole by extensional deformation. Our observations support a contention that anisotropy may be ubiquitous in continental crust.

Description: Abstract Quantitative estimates of natural climate variability are required to detect anthropogenic climate trends in the tropical Pacific, however instrumental records from this region are too short and scarce. Coral oxygen isotopic (δ18O) and strontium to calcium (Sr/Ca) records are often used to extend instrumental observations, however differences in the mean Sr/Ca and δ18O values of Porites spp. colonies from the same reef can introduce large uncertainties in coral‐based climate reconstructions. To quantify intercolony variability at Palmyra Atoll, we generate monthly‐resolved Sr/Ca and δ18O timeseries from five Porites spp. colonies that grew between 1980‐2010. Monthly to interannual variability in Sr/Ca and δ18O is well‐reproduced among different colonies, however we document intercolony offsets in mean Sr/Ca of ±0.09 mmol/mol (1σ) or ~1°C, and in mean δ18O of ±0.12‰ (1σ) or ~0.1°C. The sensitivity of each proxy to climate also varies across colonies, with Sr/Ca‐SST slopes ranging from ranging from ‐0.06 to ‐0.1 mmol mol‐1 °C‐1 and δ18O‐SST slopes ranging from ‐0.25 to ‐0.35 ‰°C‐1. Intercolony variability in both coral Sr/Ca and δ18O reduces the reproducibility of coral‐based δ18Osw reconstructions across overlapping colonies. Accounting for both intercolony variability and slope error suggests that SST reconstructions using Sr/Ca from a single Palmyra coral have an uncertainty of ±1.3°C (1σ), however replicating Sr/Ca records across multiple colonies can greatly reduce this uncertainty. A composite Sr/Ca record built using five modern cores, for example, offers a reduced error of ±0.6°C (1σ) in mean SST reconstructions, ~2.5 times smaller than errors associated with reconstructions from single corals.

Description: Abstract Hydrocarbon systems with declining or viscous oil production are often stimulated using enhanced oil recovery (EOR) techniques, such as the injection of water, steam and CO2, in order to increase oil and gas production. As EOR and other methods of enhancing production such as hydraulic fracturing have become more prevalent, environmental concerns about the impact of both new and historical hydrocarbon production on overlying shallow aquifers have increased. Noble gas isotopes are powerful tracers of subsurface fluid provenance and can be used to understand the impact of EOR on hydrocarbon systems and potentially overlying aquifers. In oil systems, produced fluids can consist of a mixture of oil, water and gas. Noble gases are typically measured in the gas phase; however, it is not always possible to collect gases and therefore produced fluids (which are water, oil and gas mixtures) must be analyzed. We outline a new technique to separate and analyze noble gases in multi‐phase hydrocarbon‐associated fluid samples. An offline double capillary method has been developed to quantitatively isolate noble gases into a transfer vessel, while effectively removing all water, oil, and less volatile hydrocarbons. The gases are then cleaned and analyzed using standard techniques. Air‐saturated water reference materials (n=24) were analyzed and results show a method reproducibility of 2.9% for 4He, 3.8% for 20Ne, 4.5% for 36Ar, 5.3% for 84Kr and 5.7% for 132Xe. This new technique was used to measure the noble gas isotopic compositions in six produced fluid samples from the Fruitvale Oil Field, Bakersfield, California.

Description: Abstract Geophysical data acquisition in oceanic domains is challenging, implying measurements with low and/or non‐homogeneous spatial resolution. The evolution of satellite gravimetry and altimetry techniques allows testing 3D density models of the lithosphere, taking advantage of the high spatial resolution and homogeneous coverage of satellites. However, it is not trivial to discretise the source of the gravity field at different depths. Here, we propose a new method for inferring tectonic boundaries at the crustal level. As a novelty, instead of modelling the gravity anomalies and assuming a flat Earth approximation, we model the Vertical Gravity Gradients (VGG) in spherical coordinates, which are especially sensitive to density contrasts in the upper layers of the Earth. To validate the methodology, the complex oceanic domain of the Caribbean region is studied, which includes different crustal domains with a tectonic history since Late Jurassic time. After defining a lithospheric starting model constrained by up‐to‐date geophysical datasets, we tested several a‐priory density distributions and selected the model with the minimum misfits with respect to the VGG calculated from the EIGEN‐6C4 dataset. Additionally, the density of the crystalline crust was inferred by inverting the VGG field. Our methodology enabled us not only to refine, confirm and/or propose tectonic boundaries in the study area, but also to identify a new anomalous buoyant body, located in the South Lesser Antilles subduction zone, and high density bodies along the Greater, Lesser and Leeward Antilles forearcs.

Description: Abstract Studying spatial and temporal trends in volcanic gas compositions and fluxes is crucial to both volcano monitoring and to constrain the origin and recycling efficiency of volatiles at active convergent margins. New compositions and fluxes are here reported for Nevado del Ruiz, Galeras and Purace, the most persistently degassing volcanoes in the Colombian Arc Segment (CAS) of the Northern Volcanic Zone (NVZ). At Nevado del Ruiz, from 2014 to 2017, plume emissions showed an average molar CO2/ST ratio of 3.9 ± 1.6 (ST is total sulfur, S). Contemporary, fumarolic chemistry at Galeras progressively shifted towards low‐temperature, S‐depleted gas discharges with an average CO2/ST ratio in excess of 10 (6.0 – 46.0, 2014‐2017). This shift in volcanic gas compositions was accompanied by a concurrent decrease in SO2 emissions, confirmed on the 21 March 2017 by high‐resolution UV camera‐based SO2 fluxes of ~2.5 kg s‐1 (~213t d‐1). For comparison, SO2 emissions remained high at Nevado del Ruiz (weighted average of 8 kg s‐1) between 2014 and 2017, while Puracé maintained rather low emission levels (〈1 kg s‐1 of SO2, CO2/SO2 ≈ 14). We here estimate carbon dioxide fluxes for Nevado del Ruiz, Galeras and Puracé of ~23, 30 and 1 kg s‐1, respectively. These, combined with recent CO2 flux estimates for Nevado del Huila of ~10 kg s‐1 (~860 t d‐1), imply that this arc segment contributes about 50% to the total subaerial CO2 budget of the Andean Volcanic Belt. Furthermore, our work highlights the northward increase in carbon‐rich sediment input into the mantle wedge via slab fluids and melts that is reflected in magmatic CO2/ST values far higher than those reported for Southern Volcanic Zone (SVZ) and Central Volcanic Zone (CVZ) volcanoes. We estimate that about 20% (~1.3 Mt C/y) of the C being subducted (~6.19 Mt C/y) gets resurfaced through subaerial volcanic gas emissions in Colombia (Nevado del Ruiz ~0.7 Mt C/y). As global volcanic volatile fluxes continue to be quantified and refined, the contribution from this arc segment should not be underestimated.

Description: Abstract Five model compounds with representative chemical structures were selected for use in simulation experiments of pyrolytic gas production. The gas production and isotopic fractionation characteristics were observed and analyzed. Then, the factors affecting carbon isotope fractionation during natural gas generation were discussed, and a fractionation model was established and calibrated. We concluded that the final hydrocarbon gas (C1‐5) yield of octadecane, octadecylamine, octadecanoic acid, decahydronaphthalene and 9‐phenylanthracene decreased in turn with the effective hydrogen content. Compared with linear alkanes or alkyl compounds, cycloalkanes have higher thermal stability and generate gas later. The variation in the carbon isotopic composition of natural gas is primarily controlled by the following three factors. a) The thermal evolution of organic matter (OM) results in a gradually heavier isotopic composition for the main gas production stage. b) Gas inherits the isotopic composition of its parent material, and this effect is evident when the chemical structure and gas generation mechanism between parent materials are similar. c) The structure of OM determines the reaction mechanism of gas generation, which has a significant influence on the range and trend of carbon isotope fractionation in the process of methane generation. An improved chemical kinetic model can accurately characterize carbon isotope fractionation during gas generation.

Description: Abstract Benthic foraminiferal assemblages and geochemical tracers (δ18O, δ13C and 14C) have been analyzed on benthic and planktonic foraminifera from core MD77‐176, located in the northern Bay of Bengal (BoB), in order to reconstruct the evolution of intermediate circulation in the northern Indian Ocean since the last glaciation. Results indicate that during the Last Glacial Maximum (LGM), Southern Sourced Water (SSW) masses were dominant at the core site. A high relative abundance of intermediate and deep infaunal species during the LGM reflects low oxygen concentration and/or meso‐ to eutrophic deep water conditions, associated with depleted benthic δ13C values. During the Holocene, benthic foraminiferal assemblages indicate an oligo– to mesotrophic environment with well‐ventilated bottom water conditions compared with LGM. Higher values for benthic foraminifera δ13C and B‐P 14C age offsets suggest an increased contribution of North Atlantic Deep Water (NADW) to the Northern BoB during the Late Holocene compared to the LGM. Millennial‐scale events punctuated the last deglaciation, with a shift in the δ13C and the ɛNd values coincident with low B‐P 14C age offsets, providing strong evidence for an increased contribution of Antarctic Intermediate Water (AAIW) at the studied site. This was associated with enhanced upwelling in the Southern Ocean, reflecting a strong sea‐atmospheric CO2 exchange through Southern Ocean ventilation during the last deglaciation.

Description: Abstract The East Taiwan Ophiolite (ETO) occurs as blocks and thrust sheets associated with the Lichi Mélange in the Coastal Range of eastern Taiwan. The blocks consist of serpentinized harzburgite, serpentinite breccia, gabbro, dikes of dolerite and plagiogranite, pillow basalts, and red clay within a mud‐ and serpentinite‐rich mélange matrix. New U‐Pb zircon dating of a pegmatite gabbro yields a weighted mean age of 16.65±0.20 Ma. This age is earlier than the North Luzon Arc, but overlaps with the late–stage spreading of the South China Sea. ETO glassy basalt has low K2O, MgO and high CaO contents, similar to MORB. REE and trace element patterns show both N‐MORB patterns with LREE depletion and E‐MORB patterns with slight LREE enrichment. A few samples show slight depletion in Nb‐Ta, and Ti and enrichment in Rb, Ba, U and Sr, indicating a hint of subduction influence. Most ETO basalt plots within the overlapping fields of N‐MORB and BABB on Ti‐V, Cr‐Y, Nb/Yb‐ Th/Yb and Hf/3‐Th‐Ta discrimination diagrams. These geochemical compositions are emblematic of mid‐ocean ridge or back‐arc lava, like South China Sea basalt. We interpret ETO basalt and gabbro as fragments of the subducted South China Sea basement that were scrapped off and accreted to the Luzon forearc during the process of subduction initiation along the Manila Trench. Blocks of mantle material in the mélange may originate from the upper plate of the arc‐continent collision and were mixed with lower plate crustal material in a subduction channel now represented by the Lichi Mélange.

Description: Abstract Magnetic properties from the Reinfjord Ultramafic Complex, in northern Norway, which formed as part of a deep magmatic conduit system, have been investigated to determine the magnetic signature of ultramafic rocks now exposed at the surface and deeper in the lower crust. The dominant carriers in these ultramafic rocks are a chrome spinel with Fe‐rich exsolution blebs, and exsolution lamellae of magnetite in clinopyroxene. Except locally, in a fault zone and in discrete small fractures, these rocks show only minor to no alteration. We infer that the magnetic oxides characterized here are representative of pristine magnetic carriers in similar rocks deeper in the crust. These oxides can be stable in lower crustal, possibly upper mantle, depths when temperatures are below the Curie temperature of magnetite, taking into account pressure effects. These ultramafic rocks are candidates for potential sources of long‐wavelength anomalies.

Description: Abstract The most recent eruption of Mt. Fuji (Japan), the VEI 5 Hōei plinian eruption (CE 1707) heavily impacted Lake Yamanaka, a shallow lake located at the foot of Mt. Fuji. Here, we discuss the influence of the Hōei eruption on the lacustrine sedimentation of Lake Yamanaka using high resolution geophysical and geochemical measurements on gravity cores. Hōei scoria fall‐out had two major impacts on Lake Yamanaka: (i) reduction of the sedimentation rate (from ~0.16 cm/yr to ~0.09 cm/yr); and (ii) the increase of in‐situ lake productivity. Sedimentation rates after the eruption were relatively low due to the thick scoria layer, trapping underlying sediments in the catchment. The lacustrine system took over more than ~170 years to begin to recover from the Hōei eruption: sedimentation recovery have been accelerated by changes in land use. Since the beginning of the 20th Century, vegetated strips delimited cultivated parcels, trapping sediment and minimizing the anthropogenic impacts on the sedimentation rate. Over the last decade, the decline of agriculture and the increase of other human activities led to an increase in the sedimentation rate (~1 cm/yr). This study highlights the effect of the grainsize of the volcanic ejecta on the sedimentation rate following a volcanic eruption. Coarse‐grained tephra are difficult to erode. Therefore, their erosion and remobilization is largely limited to intense typhoons when porous scoria deposits are saturated by heavy rains. Moreover, this study suggests that recent anthropogenic modifications of the catchment had a greater impact on the sedimentation rate than the Hōei eruption.

Description: Abstract The advantages in provenance research of U‐Pb dating different detrital minerals along with simultaneously analyzing trace elements is demonstrated in a study of sand from the mouth of the Merrimack River in New England, USA. Zircon ages record episodes of magmatism in the Early Paleozoic, peaking in the Early Devonian, followed by quiescence through the remainder of the Paleozoic and additional magmatic episodes in the Jurassic and Cretaceous. Simultaneous measurement of trace elements in zircons reveals a shift from arc magmatism to crustal melting associated with terrane collision in the Early Devonian while many Jurassic grains are clearly derived from A‐type granites. Detrital monazites and rutiles have Devonian and Permian ages. Many of the older monazites have trace element characteristics suggestive of igneous origin while Permian monazites are clearly metamorphic and record orogenesis that is absent from the detrital zircon record. Rutile grains have trace element chemistry indicative of mostly metasedimentary source rocks, and Zr thermometry indicates growth under amphibolite facies conditions. Age offsets between monazite and rutile populations provide information about the region's cooling history. Titanite grains have trace element chemistry mostly consistent with igneous origin and U‐Pb ages lining up with minor zircon age populations in the Ordovician‐Silurian and the Middle Devonian, suggesting that these magmatic episodes produced metaluminous compositions. These results show that combing trace element fingerprinting with dating and analyzing multiple detrital mineral species provide a more complete portrait of the geologic history of the sediment source region than U‐Pb dating of zircon alone.

Description: Abstract The thermal evolution of a solid planet is governed by mantle convection and therefore the dependence of viscosity on temperature. In this study, over a span of five billion years, we investigate the effect of viscosity clipping (i.e., limiting the maximum value of the viscosity) on the thermal evolution of lunar‐sized initially hot bodies featuring decaying internal heat sources. Models with a decreasing viscosity contrast resulting from limiting the maximum viscosity to 105.5 times the initial viscosity at the core‐mantle boundary were first examined. At times determined by the initial internal heating rate, rapid cooling sets in as a result of a convective regime change from stagnant‐lid to mobile‐lid convection, followed by gradual cooling to a weakly convecting and eventually nearly conductive state. Subsequently, we employ a dynamic clipping viscosity of 105.5 times the viscosity at the core‐mantle boundary, throughout the planet's evolution. In this case, stagnant‐lid convection is the only convective regime observed. Finally, convection with an initially large viscosity contrast (1010) is modeled in both 2D and 3D spherical geometry and we find strong agreement in the thermal evolution when compared with the dynamic clipping model. Our findings show that convective regime changes due to secular cooling can occur due to implementing a fixed viscosity contrast that becomes sub‐critical with respect to obtaining a stagnant‐lid. To avoid spurious convective regime changes, the specification of a dynamic clipping viscosity can be used to emulate much higher viscosity contrasts.

Description: Abstract Sand‐shale mélanges from the Kodiak Accretionary Complex and Shimanto Belt of Japan record deformation during underthrusting along a paleosubduction interface in the range 150 to 350°C. We use observations from these mélanges to construct a simple kinetic model that estimates the maximum time required to seal a single fracture as a measure of the rate of fault zone healing. Crack sealing involves diffusive redistribution of Si from mudstones with scaly fabric to undersaturated fluid‐filled cracks in sandstone blocks. Two driving forces are considered for the chemical potential gradient that drives crack sealing: 1) a transient drop in fluid pressure ∆Pf , and 2) a difference in mean stress between scaly slip surfaces in mudstones and cracks in stronger sandstone blocks. Sealing times are more sensitive to mean stress than ∆Pf, with up to four orders of magnitude faster sealing. Sealing durations are dependent on crack‐spacing, silica diffusion kinetics, and magnitude of the strength contrast between block and matrix, each of which are loosely constrained for conditions relevant to the seismogenic zone. We apply the model to three active subduction zones and find that sealing rates are fastest along Cascadia and several orders of magnitude slower for a given depth along Nicaragua and Tohoku slab‐top geotherms. The model provides: 1) a framework for geochemical processes that influence subduction mechanics via crack sealing and shear fabric development and 2) demonstration that kinetically‐driven mass redistribution during the interseismic period is a plausible mechanism for creating asperities along smooth, sediment‐dominated convergent margins.

Description: Abstract We developed a fully automated magnetic field scanner dedicated to uniaxial magnetic field measurements to determine remanent magnetization intensities and their variations in sedimentary U‐channels. A fluxgate magnetometer located as close as possible to the sedimentary section is used to perform uniaxial measurements of magnetic fields generated by the isothermal remanent magnetization of the sediment. This artificial magnetization, which is known to be a powerful proxy in environmental magnetism, is produced perpendicular to the U‐channel long axis, and parallel to the fluxgate axis, using a Halbach cylinder prior to the measurement. The present magnetic scanner offers a maximal spatial resolution of 5.8 mm for point sources. A spatial resolution of 14 mm is obtained for U channel samples. The magnetic scanner presents a reliable magnetic field range over about 3 orders of magnitude allowing measurement of magnetizations that saturate the Superconducting Rock Magnetometer in its classical configuration. The estimation of remanent magnetization intensities along the U‐channel is based on a modeling approach that uses successive uniformly magnetized prisms. In lacustrine laminated sections, comparison between modeling results based on prisms of a constant thickness, on prisms determined from sedimentary facies and on prisms determined from XRF (X‐Ray Fluorescence) data helps to understand the detrital vs. diagenetic history of the sedimentary succession.

Description: Abstract Ocean bottom seismometers (OBS) commonly record short duration events (SDEs), that could be described by all of these characteristics: (i) duration 〈 1 s, (ii) one single‐wave train with no identified P‐ nor S‐wave arrivals and (iii) a dominant frequency usually between 4 Hz and 30 Hz. In many areas, SDEs have been associated with gas or fluid‐related processes near cold seeps or hydrothermal vents, although fish bumps, instrumental or current‐generated noise have been proposed as possible sources. In order to address some remaining issues, this study presents results from in situ and laboratory experiments combined with observations from 2 contrasting areas, the Sea of Marmara (Turkey) and the Chilean subduction zone. The in situ experiment was conducted at the EMSO‐Molène submarine observatory (near Brest, France) and consisted in continuously monitoring two OBSs with a camera. The images revealed that no fish regularly bumped into the instruments. Laboratory experiments aimed at reproducing SDEs’ waveforms by injecting air or water in a tank filled by sand and sea‐water and monitored with an OBS. Injecting air in the sediments produced waveforms very similar to the observed SDEs, while injecting air in the water column did not, constraining the source of SDEs in the seafloor sediments. Finally, the systematic analysis of two real data sets revealed that it is possible to discriminate gas‐related SDEs from biological or sea‐state related noise from simple source parameters, such as the temporal mode of occurrence, the back azimuth and the dominant frequency.

Description: Abstract The variety of coral reefs morphologies highlights their sensitivities to several forcings; fossil reefs stack in sequences that are accordingly diverse. In order to understand their genesis and architectures, we devised a numerical approach, accounting for Quaternary sea‐level oscillations, vertical land motion, initial slope, wave erosion, and reef growth. We first test our model on the subsiding sequence of Hawaii, and on the uplifting sequence of Wangi‐Wangi (Sulawesi) that bears active barriers. We then construct a parametric study, that we analyse based on a comprehensive yet compact description of sequences as barcodes, that depict the vertical distribution of a few geometrical characteristics (number, width and height of the terraces, barriers). We find that geological factors suffice to explain the variety of architectures of reefal sequences at first order, regardless of additional ecosystemic processes. Vertical land motion and foundation slopes are the prime players, while reef growth rates only play a minor role. Barriers may develop both in uplift and subsidence mode, and their preservation attests for the erosional power. Last, we reappraise the genesis of sequences and find that sequences do not fingerprint discrete events of sea‐level oscillations but a continuous process harrowed by stochastic events: Major sea‐level fluctuations can be over‐represented by several terraces, or conversely absent; reoccupations may yield composite terraces representing multiple events. Overall, sequences shall not be regarded as stacks of reef bodies forming during sea‐level highstands, which implies that the commonly assumed bijective relationship between sea‐level highstands and terraces shall be abandoned.

Description: ABSTRACT Ion‐microprobe 206Pb/238U geochronology and trace element geochemistry of the unpolished rims and sectioned interiors of zircons from Yellowstone caldera's oldest post‐caldera lavas provide insight into the magmatic system during the prelude and aftermath of the caldera‐forming Lava Creek supereruption. The post‐caldera lavas compose the Upper Basin Member of the Plateau Rhyolite, and fall into two groups based on zircon crystallization age: early lavas with zircon ages between ~750‐550 ka and late lavas with zircon ages between ~350‐250 ka. Zircons from the early‐erupted East Biscuit Basin flow yield U‐Pb dates and trace element compositions, which when considered with the Pb isotopic compositions of their coexisting feldspars and pyroxenes, point to an isotopically distinct parental melt present during crystallization of the Lava Creek magma but untapped by the supereruption. Distinct zircon crystallization ages and Pb‐isotope compositions of major minerals between the early and late Upper Basin Member groups suggest contrasting sources in the magma reservoir. As proxies for melt evolution, the zircons indicate that Yellowstone's post‐caldera rhyolites became more evolved between mid‐ to late‐Pleistocene time, during the same interval that melting of hydrothermally‐altered wall rock and recharge by new silicic magmas changed in their relative roles. The results from this study indicate that discrete and ephemeral bodies of silicic magma, at times within a mush dominated reservoir and including during the prelude to the Lava Creek eruption, have characterized Yellowstone's subvolcanic reservoir.

Description: Abstract Studying offshore slow‐slip events (SSEs) along subduction zone interfaces is important for constraining the overall slip budget and potential for seismic slip, and the relationship with large megathrust earthquakes. Models using only onshore data increasingly lack model resolution the further from the shore the SSE occurs. In this study, we combine data from the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip (HOBITSS) seafloor absolute pressure gauge (APG) network with daily position timeseries from New Zealand's GeoNet to create time‐dependent models of slip during the 2014 Gisborne, New Zealand SSE using the Network Inversion Filter (NIF). We compare models assuming heterogeneous vs. homogenous elastic properties to explore their influence on our models. The time‐dependent results show that slip uncertainties under the APGs drop by about 23%. We also find that the peak value of slip increases with heterogeneous elastic properties as compared to homogenous models. The inclusion of the offshore APG data in our models places more slip near the trench and detects a more defined migration of slip, especially in the heterogeneous model. These differences are important for interpreting the relationship between the SSE and associated tremor, which occurs after the peak SSE slip‐rate. Additionally, we use a static “potency bounding” technique in order to gauge the range of models that can fit the data. This analysis demonstrates that the inclusion of offshore data helps to substantially narrow the range of plausible slip models.

Description: Abstract For ~82 million years, the Hawaiian‐Emperor chain volcanoes have sampled the Pacific mantle via the Hawaiian mantle plume, providing evidence that its composition varies on a range of temporal and spatial scales. Hawaiian volcanoes from 0 to 2 Ma are divided into southwestern (Loa) and northeastern (Kea) geographic and geochemical trends that are interpreted to reflect the bilateral chemical structure of the underlying plume and its corresponding deep mantle sources. Older Hawaiian volcanoes that formed between 8 and 3 Ma record a geochemical transition between the Kea‐dominated Northwest Hawaiian Ridge (8 to 49 Ma) and the bilateral trends of the younger Hawaiian Islands. High‐precision Pb isotopic analyses conducted on 55 new shield‐stage samples from two of these key volcanoes, Kauaʻi and Waiʻanae, show that Loa‐like Pb isotopic ratios (e.g., elevated 208Pb*/206Pb*) gradually increase with decreasing age among the northern Hawaiian volcanoes and dominate for over two million years prior to the onset of the bilateral Loa and Kea geochemical trends. Distinct isotopic groups are observed across Kauaʻi and the distribution of Loa and Kea isotopic compositions is rotated relative to that observed on the younger Hawaiian Islands. Protracted Loa compositions and the atypical Loa‐Kea trend on Kauaʻi are accounted for by: 1) the arrival of a voluminous, Loa mantle heterogeneity possibly associated with anchoring of the Hawaiian plume to the Pacific Large Low Shear Velocity Province, and 2) a different orientation of the Pacific plate relative to the Loa‐Kea compositional boundary prior to 2 Ma.

Description: Abstract Experiments have been conducted in which CO2 gases with varying C and O isotopic compositions and with stochastic and non‐stochastic Δ47 values have been allowed to equilibrate with phosphoric acid of two concentrations in reaction vessels of varying dimensions at temperatures of 25 and 90oC. Rates of 13C18O and 18O exchange between the CO2 and the phosphoric acid varied as a function of the length of exposure, volume of reaction vessel, acid strength, and difference of the initial Δ47 and δ18O values of the CO2 from theoretical equilibrium values. The Δ47 values were also altered by heated stainless steel surfaces such as those found within the Kiel Device and other preparation systems. These results have been used to explain variations in the differences in the fractionation between 25 and 90oC reported for calcite by different workers as well as differences in the slopes between temperature and Δ47 values produced by reacting samples at different temperatures (25 and 90oC).

Description: Abstract The Louisville Seamounts display age progressive volcanism and are thought to have formed as the Pacific plate moved over a long‐lived primary hotspot. Here we present 70 new 40Ar/39Ar age results from the Integrated Ocean Drilling Program (IODP) Sites U1372, U1375, U1376 and U1377 drilled and cored during Expedition 330 to the northern, older end of the Louisville Seamounts. The five seamounts drilled are flat‐topped guyots with ages ranging from ∼74 Ma (Canopus Guyot) to ∼51 Ma (Hadar Guyot) recovering up to ∼510 m of basaltic material beneath thin sediment interfaces. Our 40Ar/39Ar measurements reveal that throughout each drill hole most dates are constant within ∼500‐900 kyr at the 2σ confidence interval. In this study we use the new Louisville age information to compare against the Hawaiian‐Emperor trail on the Pacific Plate. The Louisville hotspot trail is low volume, dominantly alkali basalt, and seamounts have a relatively short lifespan up to ∼4 Myr, whereas the Hawaiian hotspot trail is high volume, has a tholeiitic shield‐building stage capped by an alkalic post‐shield stage, and have a lifespan up to ∼6‐7 Myr. Here we show a new approach to estimating a seamount's inception age based on the known cumulative age distributions for seamounts in the Louisville and Hawaii‐Emperor seamount trails. Based on our new 40Ar/39Ar analyses we conclude that existing absolute plate motion models misrepresent the age progression of Louisville seamounts and that the timing of the Hawaiian‐Emperor Bend and the 169°W bend in the Louisville seamounts are asynchronous by ~3.7 Myr.

Description: Abstract Ciomadul is the youngest volcano in the Carpathian‐Pannonian Region, Eastern‐Central Europe, which last erupted 30 ka. This volcano is considered to be inactive, however, combined evidence from petrologic and magnetotelluric data, as well as seismic tomography studies, suggests the existence of a subvolcanic crystal mush with variable melt content. The volcanic area is characterized by high CO2 gas output rate, with a minimum of 8.7 × 103 t/year. We investigated 31 gas emissions at Ciomadul to constrain the origin of the volatiles. The δ13C–CO2 and 3He/4He compositions suggest the outgassing of a significant component of mantle‐derived fluids. The He isotope signature in the outgassing fluids (up to 3.10 Ra) is lower than the values in the peridotite xenoliths of the nearby alkaline basalt volcanic field (R/Ra 5.95 Ra ± 0.01), which are representative of a continental lithospheric mantle and significantly lower than MORB values. Considering the chemical characteristics of the Ciomadul dacite, including trace element and Sr–Nd and O isotope compositions, an upper crustal contamination is less probable, whereas the primary magmas could have been derived from an enriched mantle source. The low He isotopic ratios could indicate a strongly metasomatized mantle lithosphere. This could be due to infiltration of subduction‐related fluids and postmetasomatic ingrowth of radiogenic He. The metasomatic fluids are inferred to have contained subducted carbonate material resulting in a heavier carbon isotope composition (δ13C is in the range of −1.4‰ to −4.6‰) and an increase of CO2/3He ratio. Our study shows the magmatic contribution to the emitted gases.

Description: Abstract The Talamanca Cordillera in the Central America Arc (Costa Rica‐Panama) preserves the record of the geochemical evolution from an intraoceanic arc to a juvenile continental arc in an active subduction zone, making it a testbed to explore processes that resulted in juvenile continental crust formation and explore potential mechanisms of early continental crust generation. Here we present a comprehensive set of geochronological, geochemical, and petrological data from the Talamanca Cordillera that tracks the key turning point (12–8 Ma) from the evolution of an oceanic arc depleted in incompatible elements to a juvenile continent. Most plutonic rocks from this transition and postintrusive rocks share striking similarities with average upper continental crust and Archean tonalite, trondhjemite, and granodiorite. We complement these data with seismic studies across the arc. Seismic velocities within the Caribbean Plate (basement of the arc) show a relatively uniform lateral structure consistent with a thick mafic large igneous province. Comparisons of seismic velocity profiles in the middle and lower crust beneath the active arc and remnant Miocene arc suggest a transition toward more felsic compositions as the volcanic center migrated toward the location of the modern arc. Seismic velocities along the modern arc in Costa Rica compared with other active arcs and average continental crust suggest an intermediate composition beneath the active arc in Costa Rica closer to average crust. Our geochemical modeling and radiogenic isotopes systematics suggest that input components from melting of the subducting Galapagos hotspot tracks are required for this compositional change.

Description: Abstract Several recent studies have employed variations in the concentration and isotopic composition of molybdenum as tracers of igneous processes. In this study we present new Mo concentration and δ98/95Mo data on the peculiar subduction‐related potassic magmas of the Central‐Southern Italian peninsula; the leucite‐free (lamproite‐like) rocks of the Tuscan Magmatic Province and the leucite‐bearing rocks of Mt. Vesuvius. These rocks display exotic and distinctive geochemical and isotopic features due to differences in the lithology of the subducted material in their respective mantle sources. We examine the elemental and isotopic systematics of Mo in the context of these geochemical variations. The two different associations of magmas display significantly different Ce/Mo values but surprisingly similar δ98/95Mo values (0.10–0.26‰ for Vesuvius and 0.07–0.24‰ for Tuscan Magmatic Province), which are significantly heavier than typical mid‐ocean ridge basalts. While the δ98/95Mo implicate an isotopically heavy sedimentary component recycled into their respective mantle sources, the different Ce/Mo ratios reflect contrasting elemental fractionation during sediment melting related to the lithology and consequent residual mineralogy (sulfides vs. epidote) of the subducted sedimentary material undergoing melting (Ca poor vs. Ca rich). This indicates that the heavy Mo isotopic signature of these magmas is independent of the lithology of the recycled material, which instead controls the elemental fractionation of Mo.

Description: Abstract Magma transfer from the mantle to the crust in arcs is an important step in the global cycling of elements and volatiles from Earth's interior to the atmosphere. Arc intrusive rocks dominate the total magma mass budget over extrusive rocks. However, their total volume and rate of addition is still poorly constrained, especially in continental arcs. We present lateral (forearc to backarc) and depth‐dependent (volcanics to deep crust) magma volume additions and arc‐wide magma addition rates (MARs) calculated from three continental arc crustal sections preserving magma flare‐up periods. We observe an increase in volume addition with depth and less magma added in the forearc (~15%) and backarc (~10% to 30%) compared to the main arc. Crustal‐wide MARs for each section are remarkably similar and around 0.7–0.9 km3/km2/Ma. MARs can be used to estimate CO2 fluxes from continental arcs. With initial magma CO2 contents of 1.5 wt.%, global continental arc lengths, and MARs, we calculate changes in C (Mt/year) released from continental arcs since 750 Ma. Calculated present‐day global C fluxes are similar to values constrained by other methods. Throughout the Phanerozoic, assuming equal durations of flare‐up and lull magmatism, calculated continental CO2 flux rates vary between 4 and 18 Mt C/year with highest values in the Mesozoic. These fluxes are considered minima since the intake of mantle and/or crustal carbon is not considered. Magmatic episodicity in continental arcs and changes in arc thickness and width are critical to consider when calculating MARs through time.

Description: Abstract First‐order variations in sea level exhibit amplitudes of ∼200 m over periods that coincide with those of supercontinental cycles (∼300–500 Myr). Proposed mechanisms for this sea level change include processes that change the container volume of the ocean basins and the relative elevation of continents. Here we investigate how unbalanced rates of water exchange between Earth's surface and mantle interior, resulting from fluctuations in tectonic rates, can cause sea level changes. Previous modeling studies of subduction water fluxes suggest that the amount of water that reaches sub‐arc depths is well correlated with the velocity and age of the subducting plate. We use these models to calibrate a parameterization of the deep subduction water flux, which we together with a parameterization of mid‐ocean ridge outgassing, then apply to reconstructions of Earth's tectonic history. This allows us to estimate the global water fluxes between the oceans and mantle for the past 230 Myr and compute the associated sea level change. Our model suggests that a sea level drop of up to 130 m is possible over this period and that it was partly caused by the ∼150Ma rift pulse that opened the Atlantic and forced rapid subduction of old oceanic lithosphere. This indicates that deep water cycling may be one of the more important sea level changing mechanisms on supercontinental time scales and provides a more complete picture of the dynamic interplay between tectonics and sea level change.

Description: Abstract Micrometer‐scale maps of authigenic microstructures in submarine basaltic tuff from a 1979 Surtsey volcano, Iceland, drill core acquired 15 years after eruptions terminated describe the initial alteration of oceanic basalt in a low temperature hydrothermal system. An integrative investigative approach uses synchrotron source X‐ray microdiffraction (μXRD), microfluoresence (μXRF), micro‐computed tomography (μCT), and scanning transmission electron microscopy (S/TEM) coupled with Raman spectroscopy to create finely resolved spatial frameworks that record a continuum of alteration in glass and olivine. Micro‐analytical maps of vesicular and fractured lapilli in specimens from 157.1, 137.9, and 102.6 m depth, and borehole temperatures of 83, 93.9 and 141.3 °C measured in 1980, respectively, describe the production of nanocrystalline clay mineral, zeolites, and Al‐tobermorite in diverse microenvironments. Irregular alteration fronts at 157.1 m depth resemble microchannels associated with biological activity in older basalts. By contrast, linear microstructures with little resemblance to previously described alteration features have nanocrystalline clay mineral (nontronite) and zeolite (amicite) texture. The crystallographic preferred orientation rotates around an axis parallel to the linear feature. Raman spectra indicating degraded and poorly‐ordered carbonaceous matter of possible biological origin are associated with nanocrystalline clay mineral in a crystallographically‐oriented linear microstructure in altered olivine at 102.6 m and with sub‐circular nanoscale cavities in altered glass at 137.9 m depth. Although evidence for biotic processes is inconclusive, the integrated analyses describe the complex organization of previously unrecognized mineral texture in very young basalt. They provide a foundational mineralogical reference for longitudinal, time‐lapse characterizations of palagonitized basalt in oceanic environments.

Description: Abstract A certain type of deep‐sea sediment exhibits very high content of rare‐earth elements and yttrium (REY) and is therefore expected to serve as a novel resource for these industrially essential metals. In this paper, we statistically analyzed the bulk chemical composition of deep‐sea sediments collected from the western North Pacific Ocean. By applying independent component analysis (ICA) to the multi‐elemental dataset, we extracted three independent components (ICs) that can be interpreted as the influence of Mn‐oxides (IC1), REY‐enriched biogenic calcium phosphate (BCP) (IC2), and possibly a diagenetic effect involving Cu‐enrichment (IC3) on bulk sediment geochemistry. Subsequently, we selected representative samples based on the ICA result, and implemented Sr–Nd–Pb isotopic analyses of bulk sediments. The results indicate that the extremely REY‐rich mud characterized by IC2 inherits the geochemical signature of deep Pacific seawater, whereas the non‐REY‐rich mud with less diagenetic alterations, characterized by IC3, implies an influence of terrigenous dust probably from the Taklimakan Desert–Chinese loess plateau. IC1 may reflect the variation in sedimentation rates. Characteristic downhole variations of IC1 and IC3 scores imply the presence of hiatus and/or erosion of the sediment column across the REY content peak. The putative cause is an enhanced bottom current, which can physically separate coarse BCP grains with very high REY content and thus produce an extremely REY‐enriched sediment layer.

Description: Abstract Carbon dioxide draw‐down resulting from enhanced chemical weathering during orogenesis has been invoked to explain late Cenozoic global cooling. Establishing chemical weathering records from the India–Asia collision zone is important to test this hypothesis because uplift of the Tibetan Plateau is thought to be responsible for Cenozoic cooling. However, proxies for the intensity of chemical weathering can be affected by additional factors, such as sediment grain size and provenance. Here we report major element compositions and calculated chemical weathering intensity records of three size fractions (0–5, 5–20, 20–63 μm) from the Dahonggou section of the Qaidam Basin in the northeastern Tibetan Plateau, and compare those records with published provenance data from the same section. Results show that the indices of the fine (0–5 μm) fraction fractions vary in coordination with provenance shifts, but variations of the 5–20 μm and 20–63 μm fractions are less affected by provenance variations. Comparison with Upper Continental Crust reveals that some labile elements are not leached, but instead are enriched in the fine fraction, indicating that it do not faithfully record chemical weathering intensity. In addition, weathering can result in clay mineral transformation instead of elemental variations, complicating the relationship between element‐based parameters and weathering intensity. This work suggests that changes in sediment provenance must be accounted for when inferring variations in chemical weathering intensity on the basis of element‐based weathering intensity indices of the clay fraction.

Description: Abstract African basin‐and‐swell morphology is often attributed to the planform of sub‐plate mantle convection. Across North Africa, the coincidence of Neogene and Quaternary (i.e. 〈23 Ma) magmatism, topographic swells, long wavelength gravity anomalies, and slow shear wave velocity anomalies within the asthenosphere provides observational constraints for this hypothesis. Admittance analysis of topographic and gravity fields corroborates the existence of sub‐plate support. To investigate quantitative relationships between intraplate magmatism, shear wave velocity, and asthenospheric temperature, we collected and analyzed a suite of 224 lava samples from Tibesti, Jabal Eghei, Haruj, Sawda/Hasawinah and Gharyan volcanic centers of Libya and Chad. Forward and inverse modeling of major, trace, and Rare Earth elements were used for thermobarometric studies and to determine melt fraction as a function of depth. At each center, mafic magmatism is modeled by assuming adiabatic decompression of dry peridotite with asthenospheric potential temperatures of 1300–1360°C. Surprisingly, the highest temperatures are associated with the low‐lying Haruj volcanic center rather than with the more prominent Tibesti swell. Our results are consistent with earthquake tomographic models which show that the slowest shear wave anomalies within the upper mantle occur directly beneath the Haruj center. This inference is corroborated by converting observed velocities into potential temperatures, which are in good agreement with those determined by geochemical inverse modeling. Our results suggest that North African volcanic swells are primarily generated by thermal anomalies located beneath thinned lithosphere.

Description: Abstract The Yermak Plateau is one of several regions in the Arctic Ocean where paleomagnetism yields controversial results. Despite low sedimentation rates, late Pleistocene paleomagnetic excursions have been reconstructed from many cores in the region, but they are characterized by considerably longer durations when compared to established ones. Self‐reversal during maghemitization of (titano‐)magnetite has been proposed as one explanation. Rock magnetic, 14C dating, sedimentological and stable isotope (δ18O) methods were employed to three new sediment cores to put paleomagnetic results in the context of the regional stratigraphy and chronology. Coherence of lithological parameters and δ18O variations validated the ratio of anhysteretic remanent susceptibility to bulk magnetic susceptibility (κARM/κ) as a parameter for cross‐core correlation. As established by earlier studies, we use the link between glacial/interglacial cycles and κARM/κ to tune our records to a global δ18O stack, which provides age models that are independent of radicarbon ages and paleomagnetic data. Our results show that zones of negative magnetic inclination are asynchronous across the plateau. Alternating field demagnetization data revealed that negative inclinations are contained in a medium‐high coercivity (〉25‐35 mT) magnetic phase that may be the result of post‐depositional alteration of (titano‐)magnetite. We note a positive relationship between water depth and excursion duration, which may be driven by changes in water mass circulation on glacial/interglacial timescales.

Description: Abstract We have performed an extensive characterization by transmission electron microscopy (including precession electron diffraction tomography and ab initio electron diffraction refinement as well as electron energy loss spectroscopy) of anhydrous phase B (Anh‐B) formed directly from olivine at 14 GPa, 1400 °C. We show that Anh‐B, which can be considered as a superstructure of olivine, exhibits strong topotactic relationships with it. This lowers the interfacial energy between the two phases and the energy barrier for nucleation of Anh‐B which can form as a metastable phase. We have calculated the elastic and seismic properties of Anh‐B. From the elastic point of view, Anh‐B appears to be more isotropic than olivine. Anh‐B displays only a moderate seismic anisotropy quite similar to the one of wadsleyite.

Description: Abstract Slope‐gradient maps of swath bathymetry around the Hawaiian islands locate 12 shield slope breaks associated with former shorelines that are now submerged, ranging in age back to 5 Ma. The age of their drowning correlates with the waning of tholeiitic shield volcanism that ceases to repave the shoreline sometime between the beginning and the end of the late‐shield stage. Slope breaks on Mauna Loa's northeast and southwest rifts are consistent with it being in the waning stages of shield building. Superposition of shore breaks shows that the Hualālai shield is older than Mauna Kea's. We find evidence for three volcanic shields forming Ka‘ena Ridge, for a simultaneous waning of Maui‐Lāna‘i‐Kaho‘olawe shield building and the initiation of Kohala shield building by 1.3 Ma, and for Mahukona growing to just above sea level about 0.6 Ma. A contiguous shoreline slope break formed ~1.8 Ma at the end of the West Moloka‘i late‐shield stage and the beginning of the East Moloka‘i late‐shield stage. At its western end the shield dips obliquely away from, rather than toward, this slope break. Similarly, the short Hualālai slope break plunges north, which is opposite to that expected for volcanic loading to the south. Shields locally sloping away from paleo‐shorelines may be related to landslides causing flexural back‐tilting away from their unloaded footwall. The Nu‘uanu and Wailau slides and Pololū slump were mostly shield flank events that cut marginally into their former shorelines, enough to preferentially nucleate headward erosion by streams.

Description: Abstract Active rifts release large amounts of gases from deep sources to the atmosphere by advection and diffusion processes along permeable fracture zones. The objective of this study is to develop geothermal exploration concepts for areas with little or no hydrothermal surface expressions suitable for fluid sampling and analyses (e.g., hot springs, geysers, and fumaroles). In such areas, soil gas surveys can complement established geophysical and geochemical exploration. We report CO2, 222Rn (Radon) and 220Rn (Thoron) emission data and ground temperatures from the Aluto volcanic complex in the Main Ethiopian Rift to improve understanding of tectonic and volcanic controls on the existing geothermal system. This suite of gas emission measurements allows us to identify major, deep‐rooted permeable structures with active fluid circulation and identify suitable drilling targets for geothermal production wells on Aluto. We show that significant differences in gas signatures (i.e., efflux and spatial pattern) can be used to identify predominantly volcanically and/or tectonically influenced compartments. Major gas emissions indicate significant fluid circulation at depth, which is typical for magmatic systems. Such high gas emissions have been observed in areas affected by major tectonic structures interacting with magmatic bodies at depth (tectono‐volcanic). Predominantly fault‐controlled sectors also show hydrothermal fluid circulation, but to a lower extent compared to tectono‐volcanic sectors. Within the Aluto volcanic complex geothermal production wells mainly target such fault‐controlled domains, whereas results of the study indicate strongest fluid circulation in tectono‐volcanic sectors. This result should be considered for the future exploration and development strategy of the site.

Description: Abstract The number of heat flow measurements at the Earth surface has significantly increased since the last global analysis Pollack et al.1993, and the most recent of them provide insights into key locations. This paper presents a new compilation, which includes approximately 70,000 measurements. Continental heat flow (67 mWm‐2) does not change significantly, but the differences are more important for oceanic heat flow. The divergence with conductive cooling models is reduced significantly for young ages of the seafloor, since the most recent measurements (92 mWm‐2) are significantly higher on average than the older ones (79 mWm‐2). This is related to a better quality and a better sampling of measurements in regions affected by hydrothermal circulation. The total Earth heat loss derived from these most recent measurements is estimated to ~40‐42 TW and represents only 3‐5 TW less than with a conductive cooling model (~45‐47 TW). Hydrothermal heat loss in the oceanic domain is estimated with a new method based on the ruggedness of the seafloor, and represents ~1.5 TW more than previous estimates. The heat flow variability on continents is so large that defining a trend with stratigraphic or tectono‐thermal age is difficult and makes extrapolation from age a poor predictor. On the other hand, additional geological and geophysical information can be combined with age for better predictions. A generalized similarity method was used here to predict heat flow on a global 0.5° × 0.5 ° grid. The agreement with local measurements is generally good and increases with the number and quality of proxies.

Description: Abstract Matsushima Bay was less affected by the 2011 Tohoku Earthquake Tsunami than other Pacific coastal areas of northeast Japan because of the sheltering effects of islands in the bay mouth. To understand the recovery from the 2011 tsunami of benthic environments in Matsushima Bay, we conducted surveys of the geochemical properties of surface sediments over five years following the earthquake and compared them to sedimentary data before the earthquake. Before the tsunami, the mud fraction (〈 63 μm size) proportion of the surface sediments varied over a wide range, whereas after the tsunami, the mud content range was relatively narrow. During 2012–2015, the mud content was linearly correlated with the total organic carbon (TOC) content, but the slope of the relationship differed from that before the tsunami. The tightly coupled mud–TOC relationship and the almost constant C/N ratios in the surface sediments suggest that the bay sediments were resuspended and transported by the tsunami. In addition, loss of some organic matter sources in the bay may partly account for the uniform C/N ratios. By 2016, the slope of the mud–TOC relationship was almost the same as the pre‐tsunami value, but the C/N ratios remained constant. These results suggest that ecosystems and biogeochemical cycles in the bay had not yet fully recovered their pre‐tsunami state. Trace element compositions of core samples indicated that sediment sources were little changed by the tsunami, probably because islands in the mouth of the bay reduced sediment transport into the bay from distant sources.

Description: Abstract Serpentinization of exhumed mantle peridotite is thought to influence the mechanics and evolution of hyperextended rifted margins. The Galicia margin is an example of where this may have occurred. New constraints on the 3‐D seismic velocities in the upper mantle beneath the Galicia S‐reflector detachment fault suggest heterogeneously distributed serpentinization (0‐90%). To the first order, areas with high degrees of serpentinization correspond with the intersections of crust‐cutting normal faults and the S‐reflector detachment, suggesting that these faults served as conduits transporting fluids to the upper mantle. The role of temperature on alteration was also investigated, as it is a key condition for serpentinization. Estimated paleo‐temperatures ranged from 55 to 140oC along the S‐reflector detachment. This temperature range may have played a role in the heterogeneity of alteration. These broader patterns of both high and low mantle alteration correspond to areas of thick and thin overburden immediately post‐rift, respectively. Several areas that are estimated to lie within the favorable paleo‐temperature range, and that have crust‐cutting fault intersections, exhibit unusually low degrees of serpentinization. This is possibly due to locally reduced fault permeability. The preservation of relationships between the degree of serpentinization and the occurrence of crust‐cutting fault intersections and post‐rift overburden thickness, implies that volumetrically significant amounts of mantle serpentinization occurred after movement ceased on the detachment. Serpentinization‐related volume increases may have caused post‐rift uplift of the detachment fault surface. Thus, the current morphology of this detachment surface may be due largely to heterogenous post‐kinematic processes, rather than primary fault‐related processes.

Description: Abstract Natural gas hydrate is often found in marine sediment in heterogeneous distributions in different sediment types. Diffusion may be a dominant mechanism for methane migration and affect hydrate distribution. We use a 1D advection‐diffusion‐reaction model to understand hydrate distribution in and surrounding thin coarse‐grained layers to examine the sensitivity of four controlling factors in a diffusion‐dominant gas hydrate system. These factors are the particulate organic carbon (POC) content at seafloor, the microbial reaction rate constant, the sediment grading pattern and the cementation factor of the coarse‐grained layer. We use available data at Walker Ridge 313 in the northern Gulf of Mexico where two ~3 m‐thick hydrate‐bearing coarse‐grained layers were observed at different depths. The results show that the hydrate volume and the total amount of methane within thin, coarse‐grained layers is most sensitive to the POC of fine‐grained sediments when deposited at the seafloor. The thickness of fine‐grained hydrate free zones (HFZs) surrounding the coarse‐grained layers are most sensitive to the microbial reaction rate constant. Moreover, it may be possible to estimate microbial reaction rate constants at other locations by studying the thickness of the HFZs using the Damköhler number. In addition, we note that sediment grading patterns have a strong influence on gas hydrate occurrence within coarse‐grained layers.

Description: Abstract Subduction is a key component of Earth's long‐term sulfur cycle; however, the mechanisms that drive sulfur from subducting slabs remain elusive. Isotopes are a sensitive indicator of the speciation of sulfur in fluids, sulfide dissolution‐precipitation reactions, and inferring fluid sources. To investigate these processes, we report δ34S values determined by secondary ion mass spectroscopy (SIMS) in sulfides from a global suite of exhumed high‐pressure rocks. Sulfides are classified into two petrogenetic groups: 1. Metamorphic, which represent closed‐system (re)crystallization from protoliths‐inherited sulfur, and 2. Metasomatic, which formed during open system processes, such as an influx of oxidized sulfur. δ34S values for metamorphic sulfides tend to reflect their precursor compositions: −4.3 to +13.5 ‰ for metabasic rocks, and −32.4 to −11.0 ‰ for metasediments. Metasomatic sulfides exhibit a range of δ34S from −21.7 to +13.9 ‰. We suggest that sluggish sulfur‐self diffusion prevents isotopic fractionation during sulfide breakdown, and that slab fluids inherit the isotopic composition of their source. We estimate a composition of −11 to +8 ‰ for slab fluids, a significantly smaller range than observed for metasomatic sulfides. Large fractionations during metasomatic sulfide precipitation from sulfate‐bearing fluids, and an evolving fluid composition during reactive transport may account for the entire ~36 ‰ range of metasomatic sulfide compositions. Thus, we suggest that sulfate is likely the dominant sulfur species in slab‐derived fluids.

Description: Abstract Submarine groundwater discharge (SGD) in reefs can be controlled regionally by structures. Since such structures are poorly characterized, their general importance to SGD is likely under‐appreciated. This study helps fill this gap by analyzing the effects and implications of faults and associated fractures on offshore hydrogeologic processes in coral reefs. The study area in Bolinao, northwestern Philippines is an extensive reef flat complex fringing a small, limestone island. The highly linear morphology of the island's western reef front is thought to indicate a fault which may serve as a pathway for SGD. Focusing on its northwestern (NW) side, we integrated measurements of seawater 222Rn activity, salinity and subsurface electrical resistivity (ER) to assess the nature and magnitude of offshore fault‐controlled regional SGD. In situ 222Rn activities comparable to previously reported reef flat values with corresponding bottom water freshening were observed at several NW reef front sites. ER features consistent with freshwater‐saturated sediments were also detected underlying these sites. These “fresh” subsurface ER features occur continuously along the NW reef front and discretely in the northern reef front, coinciding with the orientation of suspected faults in each of these reef areas. The multi‐faceted observations support the claim that faults across the area provide long and connected pathways for SGD in the reef front, which is in contrast with the typically lower SGD influence in offshore marine environments. The regional fault‐controlled SGD analyzed here is likely to be important for other similar coastal areas.

Description: Abstract A low‐viscosity crustal layer (LVCL) due to partial melt and/or hydration has been detected lying within the continental crust (e.g., in South China block), but its roles in the development of hyper‐extended margins are still not entirely understood. Using 2D thermo‐mechanical modeling, we simulate the lithospheric extension with a LVCL embedded within the continental crust. Results show that the ductile layer determines the width of highly thinned crust and the crustal thermal distribution. Detailed effects are concluded by varying the initial effective viscosity and thickness of the layer: (1) a rheologically weaker and/or thicker LVCL accommodates more deformation and thus results in a slow and distributed crustal thinning; and (2) preferential removal of the lithospheric mantle makes the crust in direct contact with the hot upwelling asthenosphere, promoting the development of highly thinned and hot continental crust. Both the widths of highly thinned crust and the thermal values increase when the rheology of the ductile layer decreases or the thickness increases. Furthermore, we obtain the crustal thickness and geothermal characteristics of the northern South China Sea (SCS) margin using reprocessed geophysical data and newly complied heat flow data. A detailed comparison of our modeling results with the northern SCS margin suggests that variation in the LVCL thickness may be a possible mechanism for the differential width of highly thinned crust and high heat flow zone along strike of the northern SCS margin.

Description: Abstract Detrital zircon (DZ) U‐Pb geochronology is a widely used provenance tool that leverages bedrock age signatures of hinterland source terranes. However, complex sediment recycling of multicycle zircon and hinterland provinces with nondiagnostic U‐Pb ages represent possible pitfalls for provenance reconstructions. Additional biases pertain to source rock zircon fertilities and insensitivity to low‐ and medium‐grade tectonothermal events that do not result in zircon generation. To bridge these inherent biases and gaps in DZ U‐Pb provenance data sets and derive more comprehensive tectonic reconstruction, this study combined U‐Pb and trace element analyses on DZ, apatite, and rutile as well as zircon (U‐Th)/He analyses from the Proto‐Zagros foreland basin in western Iran to shed light on Late Cretaceous tectonic accretion along Arabia. Integrated multimineral, multimethod data sets record formation of a 110‐ to 85‐Ma island arc within Triassic mid‐oceanic ridge crust of the Neotethys ocean, simultaneous obduction starting in Santonian‐Early Campanian times, and inversion of the Arabian rift margin. Overall, integration of these techniques constrains provenance based on multiple independent criteria including crystallization age, cooling history, and petrogenic‐geodynamic environment. This approach not only more completely described provenance signatures but also helped avoid significant pitfalls. For example, while Triassic DZ U‐Pb ages might be mistaken as input from Eurasia, zircon trace element analysis reveals a MORB signature and attributes these DZ to Neotethyan oceanic rather than Eurasian continental origin, having fundamentally different paleogeographic/tectonic implications for the Arabia‐Eurasia collision. Moreover, detrital rutile and apatite resolve and characterize multiple Paleozoic‐Mesozoic thermotectonic events not recorded by DZ U‐Pb due to their largely amagmatic nature.

Description: Abstract Three parameters: temperature, hydrogen content, and the presence of partial melt, are the dominant controls on the rheology of the convecting upper mantle. As such, they determine the dynamics that control plate tectonics and continental evolution. Since hydrogen depresses the peridotite solidus temperature, these parameters are strongly linked petrologically. We have developed a genetic algorithm code to statistically assess the likelihood that a section of upper mantle contains partial melt. This code uses magnetotelluric observations and petrological constraints on composition and solidus temperatures and allows for uncertainties in the geotherm and the electrical conductivity structure. We have applied this code to the convecting upper mantle beneath (1) a stable continent (the Superior Craton), (2) a hotspot (Tristan da Cunha), (3) stable, old oceanic lithosphere (the north‐west Pacific Ocean), and (4) young oceanic lithosphere (adjacent to the East Pacific Rise). Results show that the volume of melt in the convecting upper mantle is heterogeneous. The highest melt proportions are beneath the hotspot while little to no melt is required in the other regions. All regions show low water contents (generally 〈50 wt ppm in olivine) in the shallow convecting upper mantle, making it unlikely that water causes a large or sharp viscosity contrast between the lithosphere and the convecting mantle. Results differ significantly for different experimental olivine hydrogen conductivity models, highlighting the importance of reconciling these experimental constraints.

Description: Abstract Application of novel proxies, such as the stable isotope compositions and noble gas concentrations of fossil drip water trapped as inclusions in stalagmites, have the potential to provide unique constraints on past hydroclimate states and surface temperatures. Geochemical analysis of inclusion waters, however, requires an understanding of the three‐dimensional spatial distribution of dominantly liquid‐ versus air‐filled inclusions in a given stalagmite. Here we couple neutron computed tomography and medium‐ to high‐resolution X‐ray computed tomography to map out the three‐dimensional calcite density and distribution of liquid‐ versus air‐filled inclusions within a Sierra Nevada stalagmite (ML‐1), which formed during the last deglaciation (18.5 to 11.7 ka). Comparison of coupled neutron computed tomography‐X‐ray computed tomography results with a time series of stalagmite calcite fabrics indicates that although highest density calcite contains abundant liquid (fluid)‐filled inclusions, calcite density and fabric overall were secondary controls on the liquid inclusion distribution (LID). Furthermore, a multistatistical evaluation of the stalagmite time series indicates a significant relationship at the multicentury‐ to millennial‐scale between LID and calcite δ18O and δ13C that suggests a potential link between LID and water availability to the cave.

Description: Abstract Oceanic Anoxic Event‐2 (OAE‐2; Cenomanian‐Turonian) is characterized by extensive deposition of organic carbon‐rich deposits (black shales) in ocean basins worldwide as result of a major perturbation of the global carbon cycle. While the sedimentological, geochemical, and paleontological aspects of deep water expressions of OAE‐2 have been intensively studied in the last few decades, much less attention has been given to the coeval shallow water deposits. In this study, we present the results of a detailed facies and petrographic (optical microscope and scanning electron microscopy) and geochemical studies (δ13Ccarb, δ 13C org, δ 15Nbulk, TOC, and Rock‐Eval pyrolysis) on two key shallow marine sections from the Apennine Carbonate Platform (ACP; Italy). Here a continuous record of shallow water carbonates through the OAE‐2 interval is preserved, offering the unique opportunity to document the archive of paleoenvironmental changes in a neritic setting, at a tropical latitude and far from the influence of a large continental block. Two conspicuous intervals are characterized by abundant and closely spaced “dark” microbial laminites found at correlative stratigraphic horizons in the two stratigraphic sections. These laminites contain elevated concentrations of TOC (up to 1%) relative to microbial capping cycles laminites stratigraphically above and below. The organic matter preserved in these fine‐grained laminites is dominated by cyanobacteria remains, which accounted for most of the organic matter produced on the ACP in these intervals. Our study suggests that Tethyan carbonate platforms experienced significant biological changes during OAE‐2, alternating, in few kiloyears, between eutrophic phases dominated by microbial communities and mesotrophic/oligotrophic conditions favoring “normal” carbonate production/sedimentation. The synchronous occurrence of microbialite facies at different locations across the ACP and on other platforms worldwide (e.g., Mexico and Croatia) suggests a causal link between Large Igneous Province volcanism and the environmental conditions necessary to trigger cyanobacterial proliferation on shallow carbonate platforms.

Description: ABSTRACT Spectral variations of gravity and topography data show non‐isostatic compensation of elevation in the Southeast Carpathians that is consistent with our new residual topography calculations. Multi‐dimensional thermo‐mechanical models are created to consider a possible mantle flow component of the surface topography with varying temperature models and crustal configurations. The temperature anomalies derived from large‐scale P‐wave velocity anomalies and combined temperature model are estimated as input for the parametric numerical modeling experiments. Model results show that a surface uplift (up to ~ 1 km) develops in the southeastern Carpathians due to the rising hot upper mantle beneath the region. The upwelling flow is also in concordance with observations such as volcanism and enhanced surface heat flow in the region. Subsidence in the eastern part (0.5‐1 km over the Vrancea zone and the Focsani Basin) forms in relation to the down‐going Vrancea slab. The findings indicate the significant role of upper mantle structures on the present‐day anomalous topography related to proposed post‐collisional tectonics in the Southeast Carpathians.

Description: Abstract During expedition MSM37 on the German RV Maria S. Merian, bottom water temperature and sediment temperature profiles were measured in the vicinity of North Pond (western flank of Mid‐Atlantic Ridge) during exploratory dives with Remotely Operated Vehicle Jason II. In addition, push cores were taken at locations with high sediment temperature gradients. We could identify two locations where sediment temperature gradients exceed 1 K/m and bottom water temperatures showed an anomaly of up to 0.04 °C above background. We interpret these observations as clear indication of low‐temperature diffuse venting of fluids that have traveled through the uppermost crust. We can safely assume that the observed phenomena are widespread at ridge flank settings where sediment cover is thin or absent, and hence, we can explain the efficient heat mining on ridge flanks. Due to the difficulties of locating diffuse low‐temperature discharge sites and due to the fact that discharge can occur through thin sediment cover as well as through sediment‐free basement outcrops, it will be very difficult to quantify fluxes of energy and mass from low‐temperature diffuse venting in ridge flank settings; however, thermal anomalies may be used to locate sites of discharge for geochemical, microbial, and hydrologic characterization.

Description: Abstract Although seismic velocity and electrical conductivity are both sensitive to temperature, thermal lithosphere properties are derived almost exclusively from seismic data because conductivity is often too strongly affected by minor highly conductive phases to be a reliable indicator of temperature. However, in certain circumstances, electrical observations can provide strong constraints on mantle temperatures. In the southeastern United States (SEUS), magnetotelluric (MT) data require high resistivity values (〉300 Ωm) to at least 200‐km depth. As dry mantle mineral conduction laws provide an upper bound on temperature for an observed resistivity value, the only interpretation is that lithospheric temperatures (〈1330 °C) persist to 200 km. However, seismic tomography shows that velocities in this region are generally slightly slow with respect to references models; this observation has led to a view of relatively thin (〈150 km), eroded thermal lithosphere beneath the SEUS. We show that MT and seismic (tomography, attenuation, receiver function) results are consistent with thick (~200 km), coherent thermal lithosphere in this region. Reduced seismic velocities (relative to reference models) can be explained by considering the effect of finite grain size (anelasticity). Calculated velocity as a function of temperature is overall slower when including anelastic effects, even at reasonable grain sizes of 1 mm to 1 cm; this permits mantle temperatures that are colder than would typically be inferred. We argue for a geodynamic scenario in which the present thermal lithosphere in the SEUS formed in association with the Central Atlantic Magmatic Province and has subsequently survived intact for ~200 Ma.

Description: Abstract Paleomagnetic secular variations in South Asia are not well understood due to insufficient records. In this study, we produce a new paleomagnetic secular variation curve covering the past 40,000 years by stacking four marine sediment records from the middle portion of Bay of Bengal. Age models are based on radiocarbon dating and core correlation. Stepwise alternating field demagnetization was used to isolate the primary remanent component. Rock magnetic measurements indicate that the primary magnetic remanence carriers are single domain and pseudo‐single domain magnetites with no sign of diagenetic alteration. Our Bengal relative paleointensity stack is broadly consistent with local records and the PISO‐1500 relative paleointensity stack (Channell et al., 2009; https://doi.org/10.1016/j.epsl.2009.03.012). Two geomagnetic excursions presented by paleointensity minima are recognized, which correlate to the Mono Lake (29–32 ka) and Laschamp excursion (38–40 ka), respectively. We conclude that the Mono Lake and Laschamp excursion are two separate events.

Description: Abstract Observations of the time lag between the last deglaciation and a surge in volcanic activity in Iceland constrain the average melt ascent velocity to be ≥50 m/year. Although existing theoretical work has explained why the surge in eruption rates increased 5–30‐fold from the steady‐state rates during the last deglaciation, they cannot account for large variations of Rare Earth Element (REE) concentrations in the Icelandic lavas. Lavas erupted during the last deglaciation are depleted in REEs by up to 70%; whereas existing models, which assume instantaneous melt transport, can only produce at most 20% depletion. Here we develop a numerical model with finite melt ascent velocity and show that the variations of REEs are strongly dependent on the melt ascent velocity. When the average melt ascent velocity is 100 m/year, the variation of La calculated by our model is comparable to that of the observations. In contrast, when the melt ascent velocity is 1,000 m/year or above, the model variation of La becomes significantly lower than observed, which explains why previous models with instantaneous melt transport did not reproduce the large variations. We provide the first model that takes account of the diachronous response of volcanism to deglaciation. We show by comparing our model calculations of the relative volumes of different eruption types (subglacial, finiglacial, and postglacial) and the timing of the bursts in volcanic eruptions with the observations across different volcanic zones that the Icelandic average melt ascent velocity during the last deglaciation is likely to be ∼100 m/year.

Description: Abstract Compact mass spectrometers can provide simultaneous multi‐species analysis with high sensitivity and precision. For volcanic gas monitoring in situ, instrumental mass spectrometry requires reliability and ruggedness combined with low power usage and portability. The Volcanic Gas Analytical Monitor (VGAM) is capable of quantitative molecular analysis of a variety of atmospheric and volcanic gases in a single sensor by ion trap mass spectrometry. These gases include: H2, He, H2O, N2, O2, Ar, NO, N2O, CO, CO2, H2S, SO, SO2, and CH4. Unlike previous field instruments using magnetic sector and quadrupole mass spectrometers (MS) with vacuums backed by compact turbomolecular pumps, the VGAM uses a low‐power autoresonant ion trap MS and NEG‐Ion vacuum that operates at only 25 W total power, which is often the power requirement for just the mass spectrometer. Ratio‐metric mass spectral response is combined with total pressure measurements to report absolute partial pressures. Data are generated in real time and are recorded to internal flash memory. Relatively low power (〈1 Watt) data telemetry to remote sites is possible. Analysis of volcanic plumes, fumaroles, and solfatara fields is accomplished by direct inlet of gases, after trapping as much excess water vapor in both the external and internal foreline as possible. We report herein on the VGAM auto‐run and post‐processing procedures, initial calibration results for CO2, and the results of a brief field deployment at Sulfur Banks solfatara field, Kilauea Volcano, Hawaii.

Description: Abstract The scarcity of reliable paleomagnetic secular variation (PSV) records from East Asia especially from low‐latitude regions impedes better understanding of global PSV mechanisms. Here we report on a radiocarbon‐dated Holocene PSV record from a composite ~6.7‐m‐long core collected from the high‐sedimentation‐rate Huguangyan Maar Lake (HML) in subtropical‐tropical South China. Rock magnetic results demonstrate that the natural remanent magnetization resides in single‐domain magnetite. Alternating field demagnetization experiments at 1‐cm spacing on u‐channel samples reveal six distinct inclination highs at ~7,500 BCE, ~5,100 BCE, ~4,600 BCE, ~3,600–3,400 BCE, ~1,600–1,200 BCE, and 600–800 CE; three inclination lows at ~4,800 BCE, ~600–300 BCE, and ~1,000–1,300 CE; and three eastward declination trends at ~3,600–3,200 BCE, ~2,600–2,400 BCE, and 400 BCE to 200 CE. The similarity between the HML PSV record and other independently dated records from East Asia and geomagnetic field models corroborates the robustness of our age model and Holocene PSV record. Strikingly, centennial‐ to millennial‐scale PSV features of the HML are comparable, within age uncertainties, with other Holocene records from Europe, North America, and Canada, suggesting that such directional patterns are likely to be hemispheric in scale. Although relative paleointensity data of HML are affected by environmental factors (e.g., organic matter diagenesis), the record still provides a regionally important new PSV reference curve whose conspicuous features may serve as stratigraphic markers for East Asian paleorecords.

Description: Abstract Magma emplacement in organic‐rich sedimentary basins is a main driver of past environmental crises. Using a 2‐D numerical model, we investigate the process of thermal cracking in contact aureoles of cooling sills and subsequent transport and emission of thermogenic methane by hydrothermal fluids. Our model includes a Mohr‐Coulomb failure criterion to initiate hydrofracturing and a dynamic porosity/permeability. We investigate the Karoo Basin, taking into account host‐rock material properties from borehole data, realistic total organic carbon content, and different sill geometries. Consistent with geological observations, we find that thermal plumes quickly rise at the edges of saucer‐shaped sills, guided along vertically fractured high‐permeability pathways. Contrastingly, less focused and slower plumes rise from the edges and the central part of flat‐lying sills. Using a novel upscaling method based on sill‐to‐sediment ratio, we find that degassing of the Karoo Basin occurred in two distinct phases during magma invasion. Rapid degassing triggered by sills emplaced within the top 1.5 km emitted ~1.6·103 Gt of thermogenic methane, while thermal plumes originating from deeper sills, carrying a 13‐times‐greater mass of methane, may not reach the surface. We suggest that these large quantities of methane could be remobilized by the heat provided by neighboring sills. We conclude that the Karoo large igneous province may have emitted as much as ~22.3·103 Gt of thermogenic methane in the half million years of magmatic activity, with emissions up to 3 Gt/year. This quantity of methane and the emission rates can explain the negative δ13C excursion of the Toarcian environmental crisis.

Description: Abstract Serpentinite fault rheology is fundamental to tectonic and earthquake processes, yet links between deformation textures and strength evolution during fault initiation are poorly constrained. Here I present field and petrographic microstructural observations of unsheared and sheared serpentinite that demonstrate a progression of fault development. I compliment observations with a clast‐size distribution analysis to investigate the evolution of fault rigidity and a numerical model to query the stress distribution of a clast‐in‐matrix geometry. Unsheared microstructures reveal well‐aligned, elongate serpentine in the matrix and short, randomly oriented serpentine in clasts. Sheared matrix displays cataclastic textures, discrete brittle surfaces, dissolution bands, and ductile textures defined by anastomosing networks of well‐aligned, fine‐grained serpentine. During fault initiation, matrix serpentine anisotropy promotes slip on basal planes or fiber aggregates, and clast‐on‐clast interactions drive a high bulk viscosity prone to stick‐slip behavior. As deformation progresses, clast fracturing is focused at clast tips, and smaller clasts are preferentially removed by dissolution‐precipitation processes, increasing the relative abundance of matrix and promoting fault creep. Strain is continually focused in the matrix and as clast content reduces the bulk viscosity drops. This study reveals that viscosity contrasts formed by primary serpentinization textures are a major driver for the development and strength evolution of faults. On a continental‐scale, similar processes may govern earthquake distributions, fault growth, and segmentation patterns on young, serpentinite‐hosted faults.

Description: Abstract The Punchbowl fault, California, USA, is an example of a simple fault zone that has a relatively narrow fault core with further slip localization to principal slip zones. We sampled the principal slip zone, fault core, and wall rocks, as well as conducted hydrous pyrolysis experiments to analyze biomarker thermal maturity within the Punchbowl fault. Using biomarker maturity as a proxy for temperature rise, we show that the existing principal slip zone experienced greater temperature rise than the surrounding fault core and wall rock, and therefore we infer that earthquake slip localized along this layer. Furthermore, evidence of slight thermal maturity within the ultracataclasite suggests that the fault core is made up at least in part of reworked former principal slip zones. Using a wide range of possible layer thicknesses, we find that the maximum temperature range during a single earthquake could have varied from ~460‐1060 °C at 1 m/s slip velocity. However, not all samples from within the principal slip zone show a temperature rise, indicating that either layer thickness, slip, or shear stress varied during slip. Our temperature estimate also allows us to constrain the frictional energy dissipated during the earthquake to 2.2‐25 MJ/m2. Our study demonstrates that localized slip structures can be directly linked to seismicity, and that small changes in earthquake or fault parameters can lead to changes in temperature (and likely fault strength) at small scales.

Description: Abstract The clumped isotopic composition of carbonate‐derived CO2 (denoted Δ47) is a function of carbonate formation temperature, and in natural samples can act as a recorder of paleoclimate, burial or diagenetic conditions. The absolute abundance of heavy isotopes in the universal standards VPDB and VSMOW (defined by four parameters: R13VPDB, R17VSMOW, R18VSMOW and λ) impact calculated Δ47 values. Here, we investigate whether use of updated and more accurate values for these parameters can remove observed inter‐laboratory differences in the measured T‐Δ47 relationship. Using the updated parameters, we reprocess 14 published calibration datasets measured in 11 different laboratories, representing many mineralogies, bulk compositions, sample types, reaction temperatures, and sample preparation and analysis methods. Exploiting this large composite dataset (n=1253 sample replicates), we investigate the possibility for a “universal” clumped isotope calibration. We find that applying updated parameters improves the T‐Δ47 relationship (reduces residuals) within most labs and improves overall agreement but does not eliminate all inter‐laboratory differences. We reaffirm earlier findings that different mineralogies do not require different calibration equations, and that cleaning procedures, method of pressure baseline correction, and mass spectrometer type do not affect inter‐laboratory agreement. We also present new estimates of the temperature dependence of the acid digestion fractionation for Δ47 (Δ*25‐X), based on combining reprocessed data from four studies, and new theoretical equilibrium values to be used in calculation of the empirical transfer function. Overall, we have ruled out a number of possible causes of inter‐laboratory disagreement in the T‐Δ47 relationship, but many more remain to be investigated.

Description: Abstract Mt. Etna is among the largest global volcanic outgassers with respect to carbon and sulfur, yet questions remain regarding the source of these volatiles and their systematics in the crust and mantle. The importance of heterogeneous mantle sources, mixing, crustal assimilation and disequilibrium degassing are investigated using melt inclusions erupted during the A.D. 1669 eruption of Mt. Etna, Italy. We find that the melt inclusion compositions define a mixing array between two geochemically distinct melts. One end‐member melt is depleted in light rare Earth elements (LREE) and enriched in strontium (Sr), carbon and sulfur; the other is enriched in LREE and depleted in Sr, carbon and sulfur. We infer, through modeling, that the melts may either have been generated by melting a mantle source that includes a recycled oceanic crustal component; or they may have assimilated carbonate material in the crust. The resulting LREE‐depleted, Sr‐enriched melts were also alkali‐rich, which enhanced the solubility of carbon and sulfur. The LREE‐depleted, Sr‐ and volatile‐rich melt ascended through the crust and likely became supersaturated with respect to CO2 and sulfur. The melt intruded into a LREE‐enriched, relatively degassed magma body in the shallow crust, cooled rapidly and vesiculated, likely triggering eruption. The melt inclusion array trapped by growing olivines during this intrusion process records a snapshot of incomplete mixing between the two melts. Mt. Etna is renowned for the large increases in CO2 gas fluxes shortly before and during eruption. The intrusion of supersaturated, CO2‐enhanced magmas into shallow reservoirs may be a common process at Mt. Etna.

Description: Abstract Hazard analysis at caldera volcanoes is challenging due to the wide range of eruptive and environmental conditions that can plausibly occur during renewed activity. Taupo volcano, New Zealand, is a frequently active and productive rhyolitic caldera volcano that has hosted the world's youngest known supereruption and numerous smaller explosive events. To assess ashfall hazard from future eruptions, we have simulated atmospheric ash dispersal using the Ash3d model. We consider five eruption scenarios spanning magma volumes of 0.1–500 km3 and investigate the main factors governing ash dispersal in modern atmospheric conditions. Our results are examined in the context of regional synoptic weather patterns (Kidson types) that provide a framework for assessing the variability of ashfall distribution in different wind fields. For the smallest eruptions (~0.1 km3 magma), ashfall thicknesses 〉1 cm are largely confined to the central North Island, with dispersal controlled by day‐to‐day weather and the dominance of westerly winds. With increasing eruptive volume (1–5 km3 magma), ashfall thicknesses 〉1 cm would likely reach major population centers throughout the North Island. Dispersal is less dependent on weather patterns as the formation of a radially‐expanding umbrella cloud forces ash upwind or crosswind, although strong stratospheric winds significantly restrict umbrella spreading. For large eruptions (50–500 km3 magma), powerful expansion of the umbrella cloud results in widespread ashfall at damaging thicknesses (〉10 cm) across most of the North Island and top of the South Island. Synoptic climatology may prove a useful additional technique for long‐term hazard planning at caldera volcanoes.

Description: Abstract The mixing zone between high‐temperature hydrothermal fluids and seawater produces redox gradients, promoting the development of unique ecosystems based on chemotrophy. The structure of microbial communities depends on their environment, which can vary according to space and time. Hydrothermal circulation within the oceanic crust determines the chemical composition and flow of fluids, depending on underground events (earthquakes, volcanic episodes, etc.) and impacts the development of microbial communities. This link between hydrothermal vent communities and deep geological events is the focus of the present study, the first of its kind for slow‐spreading ridge. In this study, we present a unique set of multidisciplinary data collected from 2008 to 2011 on the Eiffel Tower hydrothermal site (Lucky Strike vent field, Mid‐Atlantic Ridge, MAR). We benefit from continuous geophysical monitoring (temperature and seismicity) of the site, annual sampling of hydrothermal fluids (hot and diffuse) for geochemistry analyses, sampling of hydrothermal chimneys, and an in situ microbial colonization experiment over a year for microbial study. The high CO2 content and concentrations of major elements (Cl, Ca, and Si) and SO4 in the end‐member fluids collected in 2010 indicate that a magmatic degassing occurred between 2009 and 2010 under the Lucky Strike hydrothermal field. This is supported by the large temperature variations observed in March–April 2010. These magmatic CO2 inputs seem to have affected microbial communities colonizing the high‐temperature chimney, as well as the basalts in the more diffuse and mixed zone, promoting the development of thermophilic/anaerobic Archaea and Bacteria (Archaeoglobales, Nautiliales, and Nitratiruptoraceae).

Description: Abstract Clumped isotope studies on CO2, Δ47, i.e. excess in the isotopologue containing both 13C and 18O at mass 47, can be very useful since they can give temperature estimates independent of the bulk isotopic composition. The measurement itself can be affected by a number of items. Here we develop a data processing model to examine the effects different interferences might have on the final calculated value. It incorporates known issues, e.g. pressure baseline (PBL), 17O excess, and shifts in absolute ratios for primary reference materials and parameters used for 17O correction. We also included linearity effects as well as differences in isotopologue absolute abundances at a given m/z. What normally would be considered acceptable mass spectrometer 45R and 46R linearity can skew Δ47 results. That is, 0.04‰/V and ‐0.04‰/V linearity on 45R and 46R respectively would also cause an apparent shift in the parameters used for 17O corrections. Measurements were made on CO2(g) equilibrated with water and we were able to match up the effects seen with model results. Linearity and small differences in amplitude between sample and working reference gas affected Δ47 determination, as did apparent shifts in isotopologue abundances under different conditions. This may (partially) account for discrepancies seen in Δ47‐temperature calibrations curves between laboratories. We also developed an easy way to precisely calculate the δ13C and δ18O that works well in spreadsheets without the need for multiple iterations.

Description: Abstract A mylonitic quartzite with conjugate and synthetic shear bands was investigated by Electron Backscatter Diffraction and optical microscopy to obtain insights on recrystallization mechanisms and strain localization in quartz at plastic to semibrittle conditions close to the brittle‐ductile transition. The mylonitic quartzite deformed during late Miocene thrusting coeval with contact metamorphism in the high‐strain domains of the Calamita Schists (Elba Island, Italy). Mylonitic deformation occurred from amphibolite to lower greenschist facies conditions during cooling of the aureole. Dynamic recrystallization, dominated by the activity of dislocation creep by prism 〈a〉 slip, produced recrystallized quartz layers mantling relic large quartz porphyroclasts. Under decreasing temperature and fluid‐rich conditions, quartz porphyroclasts acted as relatively rigid bodies and fractured along synthetic and conjugate C′ shear bands. Shear bands developed along kinematically favored orientations, just locally assisted by weak crystallographic planes in quartz. Fracturing along shear bands was assisted by cataclasis and fluid infiltration enhancing fracture propagation and healing by recrystallization and authigenesis of new quartz and phyllosilicate grains. The process enhanced the propagation of and strain localization in shear bands, with the development of bands of weak phyllosilicates. Furthermore, we observed the development of a c axis preferred orientation (CPO) related to dissolution and precipitation of new grains with their c axis oriented parallel to shear bands. This study highlights the importance of the interplay between brittle and crystal‐plastic processes and fluid ingress in the semibrittle regime to understand deformation partitioning and strain localization.

Description: Abstract The East Antarctic Ice Sheet (EAIS) is underlain by a series of low‐lying subglacial sedimentary basins. The extent, geology, and basal topography of these sedimentary basins are important boundary conditions governing the dynamics of the overlying ice sheet. This is particularly pertinent for basins close to the grounding line wherein the EAIS is grounded below sea level and therefore potentially vulnerable to rapid retreat. Here we analyze newly acquired airborne geophysical data over the Pensacola‐Pole Basin (PPB), a previously unexplored sector of the EAIS. Using a combination of gravity and magnetic and ice‐penetrating radar data, we present the first detailed subglacial sedimentary basin model for the PPB. Radar data reveal that the PPB is defined by a topographic depression situated ~500 m below sea level. Gravity and magnetic depth‐to‐source modeling indicate that the southern part of the basin is underlain by a sedimentary succession 2–3 km thick. This is interpreted as an equivalent of the Beacon Supergroup and associated Ferrar dolerites that are exposed along the margin of East Antarctica. However, we find that similar rocks appear to be largely absent from the northern part of the basin, close to the present‐day grounding line. In addition, the eastern margin of the basin is characterized by a major geological boundary and a system of overdeepened subglacial troughs. We suggest that these characteristics of the basin may reflect the behavior of past ice sheets and/or exert an influence on the present‐day dynamics of the overlying EAIS.

Description: Abstract The importance of gradual erosion relative to landsliding depends predominantly on the slope angle. One factor of critical influence in landsliding along with slope angle and slope shape is the soil depth. Understanding soil depth development on steep topography is fundamental for understanding and predicting the occurrence of landsliding at threshold landscapes. We develop a model to predict soil depth that addresses both threshold and gradual processes. If erosion is a gradual process, soil depth increases until the soil production rate no longer exceeds the erosion rate, and steady state is reached. The predicted soil depth (x) is proportional to the ratio of the infiltration to the erosion rate. Identifying a predictive result for erosion as a function of slope angle (S) allows a test of both the erosion and soil production models with field observations. The same theoretical approach to soil production should be applicable when the principal erosion process is shallow landsliding. After landslides, soil recovery initially follows our predicted power law increase in time, though with increasing time background erosion processes become important. At a time equal to a landslide recurrence interval, the soil depth can exceed the steady state depth by as much as a factor 2. By comparing predicted and observed x(S) results, we show that the accessed result for erosion as a function of slope angle is accurate. Soils deeper than the depth predicted at the landslide recurrence interval are beyond the stability limit. This result suggests an important practical relevance of the new soil production function.

Description: Abstract This study investigates the impact of extreme heat wave events on long‐lived massive corals (Porites spp.) from the central Saudi Arabian Red Sea using trace element (Sr/Ca, Li/Mg, Mg/Ca, U/Ca, B/Ca, and Li/Ca) records preserved in the coral skeleton for the period between 1992 and 2012. Prior to 1998, the trace element records show strong correlations with sea surface temperature. However, during the prolonged high temperature phase associated with the 1998 El Niño event, the seasonal trace element signals were disrupted, which also coincided with a reduction in extension rates. This disruption in normally highly correlated seasonal trace element ratios‐sea surface temperature relationships was unusually long, lasting for approximately 2 years in the inner‐shelf reef site and nearly 4 years in the outer‐shelf reef site. Although the seasonal signal of trace element ratios in both cores eventually stabilized, for the inner‐shelf core the amplitude and absolute values in most trace element ratios remained significantly different compared to pre‐1998 levels. This suggests that prolonged thermal stress can induce subtle but potentially long‐lasting physiological changes that affect the elemental composition of the coral's calcifying fluid. The lack of indication of stress in the core records during later bleaching events (2003, 2005, and 2010) suggests that some of these physiological changes could have induced increased thermal tolerance, particularly for inner‐shelf corals, lending support to the capacity for corals to acclimatize.

Description: Abstract To better constrain the global carbon cycle fundamental knowledge of the role of carbon cycling on continental margins is crucial. Fjords are particularly important shelf areas for carbon burial due to relatively high sedimentation rates and high organic matter fluxes. As terrigenous organic matter is more resistant to remineralization than marine organic matter, a comprehensive knowledge of the carbon source is critical to better constrain the efficiency of organic carbon burial in fjord sediments. Here we investigated highly productive fjords in northern Norway and compare our results with both existing and new organic carbon to organic nitrogen ratios and carbon stable isotope compositions from fjords in mid‐Norway, west Svalbard, and east Greenland. The marine organic carbon contribution varies significantly between these fjords, and the contribution of marine organic carbon in Norwegian fjords is much larger than previously suggested for fjords in NW Europe and also globally. Additionally, northern Norwegian fjords show very high marine carbon burial rates (73.6 gC · m‐2 · year‐1) suggesting that these fjords are probably very distinct carbon burial hotspots. We argue that the North Atlantic Current inflow sustains these high burial rates and changes in the current strength due to ongoing climate change are likely to have a pronounced effect on carbon burial in North Atlantic fjords.

Description: Abstract Reliability of magnetic recordings of the ancient magnetic field is strongly dependent on the magnetic mineralogy of natural samples. Theoretical estimates of long‐term stability of remanence were restricted to single‐domain (SD) states, but micromagnetic models have recently demonstrated that the so‐called single‐vortex (SV) domain structure can have even higher stability that SD grains. In larger grains ( 10 μm in magnetite) the multidomain (MD) state dominates, so that large uniform magnetic domains are separated by narrow domain walls. In this paper we use a parallelized micromagnetic finite element model to provide resolutions of many millions of elements allowing us, for the first time, to examine the evolution of magnetic structure from a uniform state, through the SV state up to the development of the domain walls indicative of MD states. For a cuboctahedral grain of magnetite, we identify clear domain walls in grains as small as ∼3 μm with domain wall widths equal to that expected in large MD grains; we therefore put the SV to MD transition at ∼3 μm for magnetite and expect well‐defined, and stable, SV structures to be present until at least ∼1 μm when reducing the grain size. Reducing the size further shows critical dependence on the history of domain structures, particularly with SV states that transition through a so‐called “unstable zone” leading to the recently observed hard‐aligned SV states that proceed to unwind to SD yet remain hard aligned.

Description: Abstract Molnar and England [1990] introduced equations using a semi‐analytical approach that approximate the thermal structure of the forearc regions in subduction zones. A detailed new comparison with high resolution finite element models shows that the original equations provide robust predictions and can be improved by a few modifications that follow from the theoretical derivation. The updated approximate equations are shown to be quite accurate for a straight‐dipping slab that is warmed by heat flowing from its base and by shear heating at its top. The approximation of radiogenic heating in the crust of the overriding plate is less accurate but the overall effect of this heating mode is small. It is shown that the previous and updated approximate equations become increasingly inaccurate with decreasing thermal parameter and increasing variability of slab dip. It is also shown that the approximate equations cannot be extrapolated accurately past the brittle‐ductile transition. Conclusions in a recent paper [Kohn et al., 2018] that modest amount of shear heating can explain the thermal conditions of past subduction from the exhumed metamorphic rock record are invalid due to a number of compounding errors in the application of the Molnar and England [1990] equations past the brittle‐ductile transition. The use of the improved approximate equations is highly recommended provided their limitations are taken into account. For subduction zones with variable dip and/or low thermal parameter finite element modeling is recommended.

Description: Abstract Fluid injection into the Earth's crust can induce seismic events that cause damage to local infrastructure but also offer valuable insight into seismogenesis. The factors that influence the magnitude, location, and number of induced events remain poorly understood but include injection flow rate and pressure as well as reservoir temperature and permeability. The relationship between injection parameters and injection‐induced seismicity in high‐temperature, high‐permeability reservoirs has not been extensively studied. Here we focus on the Ngatamariki geothermal field in the central Taupō Volcanic Zone, New Zealand, where three stimulation/injection tests have occurred since 2012. We present a catalog of seismicity from 2012 to 2015 created using a matched‐filter detection technique. We analyze the stress state in the reservoir during the injection tests from first motion‐derived focal mechanisms, yielding an average direction of maximum horizontal compressive stress (SHmax) consistent with the regional NE‐SW trend. However, there is significant variation in the direction of maximum compressive stress (σ1), which may reflect geological differences between wells. We use the ratio of injection flow rate to overpressure, referred to as injectivity index, as a proxy for near‐well permeability and compare changes in injectivity index to spatiotemporal characteristics of seismicity accompanying each test. Observed increases in injectivity index are generally poorly correlated with seismicity, suggesting that the locations of microearthquakes are not coincident with the zone of stimulation (i.e., increased permeability). Our findings augment a growing body of work suggesting that aseismic opening or slip, rather than seismic shear, is the active process driving well stimulation in many environments.

Description: Abstract Dissolved major ions and Sr concentrations of the Brahmaputra River at Guwahati, India, have been investigated on weekly basis for one year to understand the seasonality in weathering pattern and relative contributions from possible solute sources. Comparison of major ion data sets from present and earlier studies for this location shows no appreciable change during last ~50 years. Elemental concentrations and Ca/Na* (Na* = Na − Cl) ratio of the Brahmaputra covary (inversely) with the water discharge; the degree of seasonality, however, is less pronounced compared to other Himalayan (e.g., Ganga and Salween) rivers. The monthly averaged Ca/Si ratios of the Brahmaputra (3.7 ± 0.2), which is lower by ~2 times than those reported earlier for the Ganga outflow (6 ± 1), show minimal (~6%) seasonal changes. Seasonal variability in Na*/K with higher ratios during non‐monsoon period has been attributed to proportionally higher Na supply from hot springs and/or alkaline salts. The silicate‐derived cations (Cats) and Sr (Srs) have been estimated using an inversion method. Although these estimated values broadly show seasonal changes, the average Cats and Srs values for the monsoon (26 ± 4% [Cats]; 26 ± 6% [Srs]) and non‐monsoon (27 ± 3% (Cats); 24 ± 4% (Srs)) seasons are statistically same. These estimates indicate a weak runoff‐weathering linkage for the Brahmaputra river. Outcomes from this study suggest that the chemical weathering intensity of this basin is more dominated by regional rapid weathering around the eastern syntaxis than the climatic (runoff) parameters.

Description: Abstract We investigated the rate‐and‐state frictional properties of simulated anhydrite‐carbonate fault gouge derived from the basal Zechstein caprock overlying the seismogenic Groningen gas reservoir in the NE Netherlands. Direct shear experiments were performed at in situ conditions of 50–150 °C and 40‐MPa effective normal stress, using sliding velocities of 0.1–10 μm/s. Reservoir pore fluid compositions were simulated using 4.4 Molar NaCl brine, as well as methane, air, and brine/gas mixtures. Brine‐saturated samples showed friction coefficients (μ) of 0.60–0.69, with little dependence on temperature, along with velocity strengthening at 50–100 °C, transitioning to velocity weakening at 120 °C and above. By contrast, gas filled, evacuated and partially brine‐saturated samples showed μ values of 0.72 ± 0.02 plus strongly velocity‐weakening behavior accompanied by stick slip at 100 °C (the only temperature investigated for gas‐bearing and dry samples). A microphysical model for gouge friction, assuming competition between dilatant granular flow and thermally activated compaction creep, captures the main trends seen in our brine‐saturated samples but offers only a qualitative explanation for our gas‐bearing and dry samples. Since the reservoir temperature is ~100 °C, our results imply high potential for seismogenic slip nucleation on faults that cross cut and juxtapose the basal Zechstein anhydrite‐carbonate caprock against the Groningen reservoir sandstone, specifically in the gas‐filled upper portion of the reservoir system.

Description: Abstract Sediment color reflectance contains important information on paleoenvironmental changes. Such data are routinely measured on sediment cores from the Southern Ocean. Their usefulness, however, is undervalued. Here for the first time, color reflectance of an ~520‐kyr sediment core retrieved off Prydz Bay (East Antarctica) was intensively investigated to correlate sediment color changes to changes in sediment compositions. Total carbon, total organic carbon, X‐ray fluorescence‐derived Fe/Ti and Br/Ti ratios, water content and mineral composition of sediment from this core, and color reflectance of another two sediment cores retrieved nearby were analyzed in this study. The raw reflectance data were transformed into their first‐order derivatives and subjected to factor analysis algorithms. A four‐factor model was selected to unmix the data set, explaining 〉80% of the total data variance. Its robustness is supported by similar factor analysis results from the other two cores. The four factors represent (from Factors 1–4) pulses of sedimentary Mn enrichment and Fe3+/Fe2+ changes in clay minerals, relative content changes of colorless components, hematite and goethite, and water, respectively. Scores of Factor 1 show high spikes at glacial terminations, recording enhanced deglacial ventilation in the abyssal Southern Ocean. Scores of Factor 2 mainly reflect biogenic components and chlorite but are further complicated by refractory terrigenous organic carbon and Mn‐oxides/‐hydroxides. Higher Scores of Factor 3 occur mainly at peak interglacial intervals, associated with enhanced transport of sediment of Antarctic origin by oceanic currents. Water content (Factor 4) in the core sediment was significantly influenced by opal content and sediment grain size compositions.

Description: Abstract Oceanic plagiogranite was first defined by Coleman in 1975 as an assemblage of felsic rocks in the ocean crust and is the sole rock type in ophiolites for age dating. However, the petrogenesis of these felsic rocks remains controversial. Some consider them as partial melting products of gabbros, while others interpret them as representing highly evolved residual melt of ocean ridge magmatism. To resolve this debate, we focus our study on felsic veins from Ocean Drilling Program Hole 735B (Leg 176) in the Southwest Indian Ridge. We carried out detailed petrography, mineral compositional analysis, bulk‐rock trace element, and Sr‐Nd‐Pb‐Hf isotope analysis and petrological modeling. These data and observations lead to the conclusion that these felsic vein lithologies (i.e., oceanic plagiogranite) are solidified residual melts after advanced extents of fractional crystallization of ocean ridge basaltic magmas, rather than partial melting products of gabbros. At the late stage of mid‐ocean ridge basaltic magma evolution, Fe‐Ti oxides appear on the liquidus to crystallize, resulting in the residual melts rapidly enriched in SiO2. Such silicic melts are buoyant and can transport through “cracking” within the lithologies of the gabbroic sections of the ocean crust in the form of felsic dykes, veins, and veinlets. The volumetrically small (~0.5%) but widespread plagiogranite veins and veinlets throughout the gabbroic sections provide evidence that multiple small melt batch intrusion and differentiation are primary mechanism and mode of ocean crust accretion at slow‐spreading ridges.

Description: Abstract How rocks deform at depth during lithospheric convergence and what are the magnitudes of stresses they experience during burial/exhumation processes constitute fundamental questions for refining our vision of short‐term (i.e. seismicity) and long‐term tectonic processes in the Earth's lithosphere. Field evidence showing the coexistence of both brittle and ductile deformation at high pressure – low temperature (HP‐LT) conditions particularly fuels this questioning. We here present 2D numerical models of eclogitic rock deformation by simple shear performed at cm‐scale. To approximate the eclogite paragenesis we considered the deformed medium as composed of two mineral phases: omphacite and garnet. We run a series of models at 2.0 GPa and 550°C for different background strain rates (from 10‐14 s‐1 to 10‐8 s‐1) and for different garnet proportions (from 0% to 55%). Results show that whole rock fracturing can occur under HP‐LT conditions for strain rates larger than ~10‐10 s‐1. This suggests that observation of brittle features in eclogites does not necessarily mean that they underwent extreme strain rate. Care should therefore be taken when linking failure of eclogitic rocks to seismic deformation. We also explore the ranges of parameters where garnet and omphacite are deforming with a different deformation style (i.e. frictional vs. viscous) and discuss our results in the light of naturally deformed eclogitic samples. This study illustrates that effective stresses sustained by rocks can be high at these P‐T conditions. They reach up to ~1 GPa for an entirely fractured eclogite and up to ~500 MPa for rocks that contain fractured garnet.

Description: Abstract To decipher sediment provenance, understand dispersal pattern and its controlling factors in the Indian Ocean, we have studied major and trace element abundances and isotope compositions of Sr (87Sr/86Sr) and Nd (εNd) in silicate fractions of 38 surface sediment samples collected at ~1.5° interval from India (~10°N) to Antarctica (~70°S). The 87Sr/86Sr ratio and εNd vary from 0.70527 and 0.1 near Crozet Island to 0.77854 and −34.8 near Antarctica, respectively. Higher 87Sr/86Sr ratio and lower εNd near Antarctica indicate source of old granitic rocks, and the opposite is observed in sediments from the Central and the Southwest Indian Ridges indicating source of seafloor spreading and mid‐ocean ridge hydrothermal activity. Chondrite normalized Rare Earth Elements (REEs) pattern reveals positive Europium anomaly for Central and Southwest Indian Ridge sediments, while positive Gadolinium anomaly is observed in Indian and Antarctic coastal sediments. Sediments in the Indian sector of Southern Ocean and coastal Antarctica are highly dispersed through deep, bottom, and Antarctic circumpolar currents. The εNd spatial map generated for the entire Indian Ocean based on compilation of present study and results from literature shows dominance of terrigenous inputs from continental erosion with minor contribution from the ridges in the Indian Ocean. The εNd map together with trace elements can be employed for reconnaissance for mineral exploration associated with hydrothermal deposits. Further, comparison of detrital εNd with its counter phases, that is, authigenic fractions and modern seawater at the sampling sites, would enhance our knowledge of benthic exchange processes, with an implication in using Nd isotopes as a water mass circulation proxy.

Description: Abstract We performed a series of extrusion experiments on partially molten samples of forsterite plus 10 vol% of an anorthite‐rich melt to investigate melt segregation in a pipe‐extrusion geometry and test the predictions of two‐phase flow theory with viscous anisotropy. The employed flow geometry has not been experimentally investigated for partially molten rocks; however, numerical solutions for a similar, pipe‐Poiseuille geometry are available. Samples were extruded from a 6‐mm diameter reservoir into a 2‐mm diameter channel under a fixed normal stress at 1350°C and 0.1 MPa. The melt distribution in the channel was subsequently mapped with optical and backscattered electron microscopy and analyzed via quantitative image analysis. Melt segregated from the center toward the outer radius of the channel. The melt fraction at the wall increased with increasing extrusion duration and with increasing shear stress. The melt fraction profiles are parabolic with the melt fraction at the wall reaching 0.17–0.66, values 2 to 16 times higher than at the channel center. Segregation of melt toward the wall of the channel is consistent with base‐state melt segregation as predicted by two‐phase flow theory with viscous anisotropy. However, melt‐rich sheets inclined at a low angle to the wall, which are anticipated from two‐phase flow theory, were not observed, indicating that the compaction length is larger than the channel diameter. The results of our experiments are a test of two‐phase flow theory that includes viscous anisotropy, an essential theoretical frame work needed for modeling large‐scale melt migration and segregation in the upper mantle.

Description: Abstract We analyze peridotites from a wide range of tectonic settings to investigate relationships between olivine crystallographic preferred orientation (CPO) and deformation conditions in naturally deformed rocks. These samples preserve the five olivine CPO types (A through E‐type) that rock deformation experiments have suggested are controlled by water content, temperature, stress magnitude, and pressure. The naturally deformed specimens newly investigated here (65 samples) and compiled from an extensive literature review (445 samples) reveal that these factors may matter less than deformation history and/or geometry. Some trends support those predicted by experimentally determined parametric dependence, but several observations disagree — namely that all CPO types are able to form at very low water contents and stresses, and that there is no clear relationship between water content and CPO type. This implies that at the low stresses typical of deformation in the mantle, CPO type more commonly varies as a function of strain geometry. Because olivine CPO is primarily responsible for seismic anisotropy in the upper mantle, the results of this study have several implications. These include (1) the many olivine CPO types recorded in samples from individual localities may explain some of the complex seismic anisotropy patterns observed in the continental mantle, and (2) B‐type CPO – where olivine's “fast axes" align perpendicular to flow direction – occurs under many more conditions than traditionally thought. This study highlights the need for more experiments, and the difficulty in using olivine CPO in naturally‐deformed peridotites to infer deformation conditions.

Description: Abstract The ArcCRUST model consists of crustal thickness and estimated crustal thinning factors grids for the High Arctic and Circum‐Arctic regions (north of 67oN). This model is derived by using 3D forward and inverse gravity modelling. Updated sedimentary thickness grid, an oceanic lithosphere age model together with inferred microcontinent rifting ages, variable crystalline crust and sediment densities, and dynamic topography models constrain this inversion. We use published high‐quality 2D seismic crustal‐scale models to create a database of Depths to Seismic Moho (DSM) profiles. To check the quality of the ArcCRUST model, we have performed a statistical analysis of misfits between the ArcCRUST Moho depths and DSM values. Systematic analysis of the misfits within the Arctic sedimentary basins provides information about tectonic processes unaccounted by the assumed model of pure‐shear lithospheric extension. In particular, our model implies a less‐dense and/or thin mantle lithosphere underneath microcontinents in the deep Arctic Ocean where the ArcCRUST depth to Moho values exceed the depth to seismic Moho. A systematically larger gravity‐derived crustal thickness (ca. 3 km) under the western and northern Greenland Sea points to a hotter upper mantle implied by the seismic tomography models in the North Atlantic.

Description: Abstract Pass‐through superconducting rock magnetometers (SRMs) enable rapid and precise remanence measurement of continuous samples and are essential for paleomagnetic studies. Due to convolution effect of the SRM sensor response, pass‐through measurements need to be deconvolved to restore accurate and high‐resolution signal. A key step toward successful deconvolution is a reliable estimate of the SRM sensor response. Here, we present new tool URESPONSE for accurate SRM sensor response estimate based on measurements of a well‐calibrated magnetic point source. URESONSE allows sensor response to be estimated for continuous samples with different cross‐section geometry. We estimate sensor responses for an old liquid helium‐cooled SRM (SRM‐old) and a new liquid helium‐free SRM (SRM‐new) at the University of Southampton and compare remanence measurement of a u‐channel on both SRMs before and after deconvolution. For each SRM, sensor response estimates based on data collected using different magnetic point source samples and/or measurement procedures generally yield small differences (std. 〈~1%), while sensor response estimates for continuous samples with different cross‐section geometry often show larger differences (std. up to ~2%). Compared with SRM‐old, SRM‐new has smaller cross‐axis responses, less negative zones, and significantly broader main axis responses. We demonstrate that normalization of data using a nine‐element “effective length” matrix calculated from sensor response estimate is necessary to minimize differences in measurements on two SRMs. Deconvolution of measurements on two SRMs using accurate sensor response estimates yields highly consistent and high‐resolution results, while deconvolution using inaccurate sensor response data can lead to significant differences especially for data from SRM‐old that has large cross‐axis responses.

Description: Abstract The Mono Craters are an overlapping chain of at least 28 domes and coulees located south of Mono Lake, east‐central California, and represents the most recent eruptions of high‐silica rhyolite magma in the Mono Lake–Long Valley volcanic region. Regionally widespread tephra‐fall deposits from the Mono Craters serve as important chronostratigraphic markers for correlations of late Quaternary terrestrial archives in the western United States. A well‐resolved eruption chronology for the Mono Craters is thus not only critical for volcanic hazard considerations, but also for paleoclimatic and paleomagnetic studies. Here, we constrain the timing of early Mono Craters volcanism using ion microprobe 238U–230Th dating of unpolished crystal faces of autocrystic zircon and allanite microphenocrysts, which samples the final stage of crystallization prior to eruption. The stratigraphically oldest domes (biotite‐bearing rhyolites) yield 238U–230Th isochron dates between ca. 42 ka and ca. 19 ka and are unambiguously linked to dated tephras in the Wilson Creek formation of ancestral Mono Lake via coupled 238U–230Th geochronology and titanomagnetite geochemistry. The second oldest set of domes (orthopyroxene‐bearing rhyolites) have overlapping 238U–230Th isochron dates that are within uncertainty of published K/Ar and 40Ar/39Ar dates for the fayalite‐bearing rhyolite domes, suggesting a period of intense effusive activity for the Mono Craters near the Pleistocene–Holocene boundary. Our new and previously published 238U–230Th dates for tephras in the Wilson Creek formation provide robust geochronologic constraints for the controversial geomagnetic excursion recorded in the Wilson Creek formation.

Description: Abstract Multibeam bathymetry compiled along fracture zones of the northern Atlantic reveals a striking morphological asymmetry. Seafloor fabric produced at western ridge‐transform intersections (RTIs) tends to consist of linear ridges, small in amplitude and regularly spaced, and oceanic core complexes (OCCs) occur infrequently. In contrast, seafloor fabric produced at eastern RTIs is more irregular and blocky and displays characteristics usually associated with melt‐poor accretion: These include more than double the occurrence of OCCs as well as greater seafloor depths, even where seafloor is younger than that on the opposite side of the fracture zone. We propose that this asymmetry is a consequence of the westward migration of the Mid‐Atlantic Ridge. Such migration is expected to result in an enhanced melt supply at leading (western) RTIs compared to trailing (eastern) RTIs. The morphological asymmetry is observed for ridge offsets of ~40 to ~200 km, a range that may be related to the width of melting regimes that supply ridge segments. At slow‐spreading ridges, contrasting melt supplies across fracture zones may therefore be best expressed in the distinct seafloor fabrics preserved beyond the dynamically maintained relief of the axial rift valley and walls rather than in the contrasting axial depths documented for faster spreading ridges. Although OCCs appear to form preferentially at spreading rates lower than 30 mm/yr, their common occurrence along the northern Mid‐Atlantic Ridge may also reflect the significantly higher rates at which it is migrating compared to the southern Mid‐Atlantic Ridge.

Description: Abstract This study reports the composition of the oceanic crust from the 16.5°N region of the Mid‐Atlantic Ridge, a spreading ridge segment characterized by a complex detachment fault system and three main oceanic core complexes (southern, central, and northern OCCs). Lithologies recovered from the core complexes include both greenschist facies and weathered pillow basalt, diabase, peridotite, and gabbro, while only weathered and fresh pillow basalt was dredged from the rift valley floor. The gabbros are compositionally bimodal, with the magmatic crust in the region formed by scattered intrusions of chemically primitive plutonic rocks (i.e., dunites and troctolites), associated with evolved oxide‐bearing gabbros. We use thermodynamic models to infer that this distribution is expected in regions where small gabbroic bodies are intruded into mantle peridotites. The occurrence of ephemeral magma chambers located in the lithospheric mantle enables large proportions of the melt to be erupted after relatively low degrees of fractionation. A large proportion of the dredged gabbros reveal evidence for deformation at high‐temperature conditions. In particular, chemical changes in response to deformation and the occurrence of very high‐temperature ultramylonites (〉1000 °C) suggest that the deformation related to the oceanic detachment commenced at near‐solidus conditions. This event was likely associated with the expulsion of interstitial, evolved magmas from the crystal mush, a mechanism that enhanced the formation of disconnected oxide‐gabbro seams or layers often associated with crystal‐plastic fabrics in the host gabbros. This granulite‐grade event was soon followed by hydrothermal alteration revealed by the formation of amphibole‐rich veins at high‐temperature conditions (~900 °C).

Description: Abstract A comprehensive characterization of gas hydrate system offshore the western Black Sea was performed through an integrated analysis of geophysical data. We detected the bottom‐simulating reflector (BSR), which marks, in this area, the base of gas hydrate stability. The observed BSR depth does not fit the theoretical steady state base of gas hydrate stability zone (BGHSZ). We show that the disparity between the BSR and predicted BGHSZ is the result of a transient state of the hydrate system due to the ongoing reequilibrium since the Last Glacial Maximum. When gas hydrates are brought outside the stability zone due to changes in temperature and sea level, their dissociation generates an increase in interstitial pore pressure. This process is favorable to the recrystallization of gas hydrates and delays the upward migration of the hydrate stability zone explaining the anomalously deep BSR. The BSR depth, which is commonly used to derive geothermal gradient values by assuming steady state conditions, is used here to derive the maximum excess pore pressure at the BGHSZ. Derived excess pore pressure values of 1–2 MPa are probably the result of the low permeability of hydrate‐bearing sediments. Higher pore pressure values derived at the location of a fault system could cause hydrofracturing enabling the free gas to cross the gas hydrate stability zone and emerge at the seafloor, forming the flares observed in close vicinity to where the shallow gas hydrates were sampled.

Description: Abstract Plagioclase‐hosted melt inclusions are infrequently used to investigate magmatic processes owing to the perception that they are less robust than olivine‐hosted inclusions. Based on a set of time series experiments ranging from 30 min to 4 days, we demonstrate that plagioclase‐hosted melt inclusions can preserve the original (primitive) magmatic signal if steps are taken to correct for post entrapment crystallization. Diffusion within plagioclase‐hosted melt inclusions is sufficiently fast that the melt inclusions can be homogenized within 30 min through heating experiments. As heating time increases, the composition of melt inclusions drifts. We attribute this longer‐term phenomenon to plastic deformation of the host plagioclase (crystal relaxation) inducing a decrease in the internal pressure in the melt inclusion combined with diffusion of elements across the host/melt inclusion interface. The rate of chemical reequilibration within melt inclusions is limited by the much slower rate of diffusion within the solid host. Indeed, the host/melt inclusion partition coefficient for MgO decreases by 25% from the 30‐min to the 4‐day experiments. Our results suggest that the primary character of plagioclase‐hosted melt inclusions can be recovered if one recognizes and corrects for the effects of the complex physical and chemical processes that occur during melt inclusion homogenization.

Description: Abstract Detailed reconstruction of Indian summer monsoons is necessary to better understand the late Quaternary climate history of the Bay of Bengal and Indian peninsula. We established a chronostratigraphy for a sediment core from Hole 19B in the western Bay of Bengal, extending to approximately 80 kyr BP and examined major and trace element compositions and clay mineral components of the sediments. Higher δ18O values, lower TiO2 contents, and weaker weathering in the sediment source area during marine isotope stages (MIS) 2 and 4 compared to MIS 1, 3, and 5 are explained by increased Indian summer monsoonal precipitation and river discharge around the western Bay of Bengal. Clay mineral and chemical components indicate a felsic sediment source, suggesting the Precambrian gneissic complex of the eastern Indian peninsula as the dominant sediment source at this site since 80 kyr. Trace element ratios (Cr/Th, Th/Sc, Th/Co, La/Cr, and Eu/Eu*) indicate increased sediment contributions from mafic rocks during MIS 2 and 4. We interpret these results as reflecting the changing influences of the eastern and western branches of the Indian summer monsoon and a greater decrease in rainfall in the eastern and northeastern parts of the Indian peninsula than in the western part during MIS 2 and 4.