Oblique rifting is investigated through centrifuge experiments that reproduce extension of a continental lithosphere containing a preexisting weakness zone. During extension, this weakness localizes deformation, and different rift obliquity is obtained by varying its trend with respect to the stretching direction. Model results show that deformation is mostly controlled by the obliquity angle a (defined as the angle between the orthogonal to the rift trend and the extension direction). For low obliquity (alpha < 45 degrees), rifting is initially characterized by activation of large, en echelon boundary faults bordering a subsiding rift depression, with no deformation affecting the rift floor. Increasing extension results in the abandonment of the boundary faults and the development of new faults within the rift depression. These faults are orthogonal to the direction of extension and arranged in two en echelon segments linked by a complex transfer zones, characterized by strike-slip component of motion. In these models, a strong strain partitioning is observed between the rift margins, where the boundary fault systems have an oblique-slip motion, and the valley floor that away from the transfer zones is affected by a pure extension. Moderate obliquity (alpha = 45 degrees) still results in a two-phase rift evolution, although boundary fault activity is strongly reduced, and deformation is soon transferred to the rift depression. The fault pattern is similar to that of low-obliquity models, although internal faults become slightly oblique to the orthogonal to the direction of extension. Deformation partitioning between the rift margins and the valley floor is still observed but is less developed than for low-obliquity rifting. For high obliquity (alpha > 45 degrees), no boundary faults form, and the extensional deformation affects the rift depression since early stages of extension. Dominance of the strike-slip motion over extension leads to the development of oblique-slip and nearly pure strike-slip faults, oblique to both the rift trend and the orthogonal to the extension direction, with no strain partitioning between the margins and the rift floor. These results suggest that oblique reactivation of preexisting weaknesses plays a major role in controlling rift evolution, architecture, and strain partitioning, findings that have a significant relevance for natural oblique rifts.

P-wave velocities (Vp) have been measured in the laboratory and calculated using thermodynamic modeling for seven representative rock samples from the lower crust to mantle section of the Kohistan paleo-island arc. Lower crustal rocks comprise plagioclase-rich gabbro, garnet-bearing gabbro, and hornblendite; mantle rocks comprise garnetite, pyroxenite, websterite, and dunite. Measurements were performed at confining pressures up to 0.5GPa and temperatures up to 1200 degrees C. Vp were also calculated using rock major element chemistry with the Perple_X software package. Calculated Vp match closely the laboratory measurements. At depths representative for the arc root, Vp of upper mantle rocks vary from 7.7-8.1km/s, whereas the lower crustal rocks have velocities between 6.9-7.5km/s. P-wave anisotropy is small, with exceptions of sheared gabbros. Measured and calculated seismic properties are consistent with, and complement a growing database of published seismic properties from the Kohistan arc. In the light of such data, we discuss seismic imaging of present-day island arcs. Intermediate Vp (7.4-7.7km/s) in arc roots can be explained by pyroxenites and garnet-bearing mafic rocks. Strong seismic reflectors may be related to garnetites (8.0-8.2km/s).

Extensive regions of yellow and white ("bleached") sandstones within the terrestrial Jurassic red bed deposits of the Colorado Plateau reflect widespread interaction with subsurface reduced fluids which resulted in the dissolution of iron-oxide grain coatings. Reduced fluids such as hydrocarbons, CO2, and organic acids have been proposed as bleaching agents. In this study, we characterize an altered section of the Slick Rock member of the Jurassic Entrada Sandstone that exposes bleached sandstone with bitumen-saturated pore spaces. We observe differences in texture, porosity, mineralogy, and geochemistry between red, pink, yellow, and gray facies. In the bleached yellow facies we observe quartz overgrowths, partially dissolved K-feldspar, calcite cement, fine-grained illite, TiO2-minerals, and pyrite concretions. Clay mineral content is highest at the margins of the bleached section. Fe2O3 concentrations are reduced up to 3x from the red to gray facies but enriched up to 50x in iron-oxide concretions. Metals such as Zn, Pb, and rare-earth elements are significantly enriched in the concretions. Supported by a batch geochemical model, we conclude the interaction of red sandstones with reduced hydrocarbon-bearing fluids caused iron-oxide and K-feldspar dissolution, and precipitation of quartz, calcite, clay, and pyrite. Localized redistribution of iron into concretions can account for most of the iron removed during bleaching. Pyrite and carbonate stable isotopic data suggest the hydrocarbons were sourced from the Pennsylvanian Paradox Formation. Bitumen in pore spaces and pyrite precipitation formed a reductant trap required to produce Cu, U, and V enrichment in all altered facies by younger, oxidized saline brines.

Fluid flow through the axial hydrothermal system at fast spreading ridges is investigated using the Sr-isotopic composition of upper crustal samples recovered from a tectonic window at Pito Deep (NE Easter microplate). Samples from the sheeted dike complex collected away from macroscopic evidence of channelized fluid flow, such as faults and centimeter-scale hydrothermal veins, show a range of ⁸⁷Sr/⁸⁶Sr from 0.7025 to 0.7030 averaging 0.70276 relative to a protolith with ⁸⁷Sr/⁸⁶Sr of ∼0.7024. There is no systematic variation in ⁸⁷Sr/⁸⁶Sr with depth in the sheeted dike complex. Comparison of these new data with the two other localities that similar data sets exist for (ODP Hole 504B and the Hess Deep tectonic window) reveals that the extent of Sr-isotope exchange is similar in all of these locations. Models that assume that fluid-rock reaction occurs during one-dimensional (recharge) flow lead to significant decreases in the predicted extent of isotopic modification of the rock with depth in the crust. These model results show systematic misfits when compared with the data that can only be avoided if the fluid flow is assumed to be focused in isolated channels with very slow fluid-rock exchange. In this scenario the fluid at the base of the crust is little modified in ⁸⁷Sr/⁸⁶Sr from seawater and thus unlike vent fluids. Additionally, this model predicts that some rocks should show no change from the fresh-rock ⁸⁷Sr/⁸⁶Sr, but this is not observed. Alternatively, models in which fluid-rock reaction occurs during upflow (discharge) as well as downflow, or in which fluids are recirculated within the hydrothermal system, can reproduce the observed lack of variation in ⁸⁷Sr/⁸⁶Sr with depth in the crust. Minimum time-integrated fluid fluxes, calculated from mass balance, are between 1.5 and 2.6 × 10⁶ kg m⁻² for all areas studied to date. However, new evidence from both the rocks and a compilation of vent fluid compositions demonstrates that some Sr is leached from the crust. Because this leaching lowers the fluid ⁸⁷Sr/⁸⁶Sr without changing the rock ⁸⁷Sr/⁸⁶Sr, these mass balance models must underestimate the time-integrated fluid flux. Additionally, these values do not account for fluid flow that is channelized within the crust.

[1] We develop a theory that relates gas hydrate saturation in marine sediments to the depth of the sulfate‐ methane transition (SMT) zone below the seafloor using steady state, analytical expressions. These expres- sions are valid for systems in which all methane transported into the gas hydrate stability zone (GHSZ) comes from deeper external sources (i.e., advective systems). This advective constraint causes anaerobic oxidation of methane to be the only sulfate sink, allowing us to link SMT depth to net methane flux. We also develop analytical expressions that define the gas hydrate saturation profile based on SMT depth and site‐specific parameters such as sedimentation rate, methane solubility, and porosity. We evaluate our analytical model at four drill sites along the Cascadia Margin where methane sources from depth dominate. With our model, we calculate average gas hydrate saturations across GHSZ and the top occurrence of gas hydrate at these sites as 0.4% and 120 mbsf (Site 889), 1.9% and 70 mbsf (Site U1325), 4.7% and 40 mbsf (Site U1326), and 0% (Site U1329), mbsf being meters below seafloor. These values compare favorably with average saturations and top occurrences computed from resistivity log and chloride data. The analyt- ical expressions thus provide a fast and convenient method to calculate gas hydrate saturation and first‐ order occurrence at a given geologic setting where vertically upward advection dominates the methane flux. 

Atmospheric mineral dust is a poorly constrained yet extremely important component of the climate system. Provenance studies from geologic dust archives are crucial to understand the drivers of the dust cycle over long time scales. Our multi-technique provenance analysis of a rare Paleogene (35–27 Ma) eolian dust sequence from Ulantatal, ∼400 km northwest of the Chinese Loess Plateau (CLP), shows that Paleogene dust transporting winds generally varied between northwesterly and westerly, the same as those in the late Neogene-Quaternary bipolar icehouse. We propose that, as today, westerly wind circulation patterns would have been modulated by an Arctic Oscillation (AO)-like situation, and that the warm Eocene favored a long-term negative phase of AO, leading to meridional westerly circulation and the dominance of a northwesterly dust transport pathway. After the Eocene-Oligocene transition (EOT), long-term positive phase of AO-like conditions initiated, leading to stronger and more zonal westerlies. The Siberian High (SH) also formed or strengthened at the EOT and started to control dust storm activity along the northwesterly transport pathway. We argue that increased Paleogene Northern Hemisphere (NH) ice volume was the ultimate driver of this modern-type dust transport regime in the Ulantatal region, possibly also controlling initial Ulantatal dust sequence formation via the development of the SH and modern-type eolian regime. The similarity between the Ulantatal and late Neogene northern CLP dust provenance signals suggests that the increased NH ice volume, via its control on the northwesterly dust transport, could have promoted increased loess formation also in the late Miocene.

Here we present the first four dimensional (time and three dimensional space resolved) experiment on a strongly deformed geological material. Results show that even complicated microstructures with large continuous and discontinuous changes in crystallographic orientation can be resolved quantitatively. The details that can be resolved are unprecedented and therefore the presented technique promises to become influential in a wide range of geoscientific investigations. Grain and subgrain scale processes are fundamental to mineral deformation and associated Earth Dynamics, and time resolved observation of these processes is vital for establishing an in-depth understanding of the latter. However, until recently, in situ experiments were restricted to observations of two dimensional surfaces. We compared experimental results from two dynamic, in situ annealing experiments on a single halite crystal; a 2D experiment conducted inside the scanning electron microscope and a 3D X-ray diffraction experiment. This allowed us to evaluate the possible effects of the free surface on grain and subgrain processes. The extent to which surface effects cause experimental artifacts in 2D studies has long been questioned. Our study shows that, although the nature of recovery processes are the same, the area swept by subgrain boundaries is up to 5 times larger in the volume than observed on the surface. We suggest this discrepancy is due to enhanced drag force on subgrain boundaries by thermal surface grooving. Our results show that while it is problematic to derive absolute mobilities from 2D experiments, derived relative mobilities between boundaries with different misorientation angles can be used.

Recent seismic unrest and a persistent Holocene eruption record at Katla volcano, Iceland indicate that a near-future eruption is possible. Previous petrological investigations suggest that Katla is supplied by a simple plumbing system that delivers magma directly from depth, while seismic and geodetic data also point toward the existence of upper-crustal magma storage. To characterize Katla's recent plumbing system, we established mineral-melt equilibrium crystallization pressures from four age-constrained Katla tephras spanning from 8 kyr BP to 1918. The results point to persistent shallow- (≤8 km depth) as well as deep-crustal (ca. 10 – 25 km depth) magma storage beneath Katla throughout the last 8 kyr. The presence of multiple magma storage regions implies that mafic magma from the deeper reservoir system may become gas-rich during ascent and storage in the shallow crust and erupt explosively. Alternatively, it might intersect evolved magma pockets in the shallow-level storage region, and so increase the potential for explosive mixed-magma ash eruptions.

Sheet intrusions (inclined sheets and dykes) in the deeply eroded volcanoes of Geitafell and Dyrfjoll, eastern Iceland, were studied at the surface to identify the location, depth, and size of their magmatic source(s). For this purpose, the measured orientations of inclined sheets were projected in three dimensions to produce models of sheet swarm geometries. For the Geitafell Volcano, the majority of sheets converge toward a common focal area with a diameter of at least 4 to 7 km, the location of which coincides with several gabbro bodies exposed at the surface. Assuming that these gabbros represent part of the magma chamber feeding the inclined sheets, a source depth of 2 to 4 km below the paleoland surface is derived. A second, younger swarm of steeply dipping sheets crosscuts this gabbro and members of the first swarm. The source of this second swarm is estimated to be located to the SE of the source of Swarm 1, below the present-day level of exposure and deeper than the source of the first swarm. For the Dyrfjoll Volcano, we show that the sheets can be divided into seven different subsets, three of which can be interpreted as swarms. The most prominent swarm, the Njardvik Sheet Swarm, converges toward a several kilometers wide area in the Njardvik Valley at a depth of 1.5 to 4 km below the paleoland surface. Two additional magmatic sources are postulated to be located to the northeast and southwest of the main source. Crosscutting relationships indicate contemporaneous, as well as successive activity of different magma chambers, but without a resolvable spatial trend. The Dyrfjoll Volcano is thus part of a complex volcanic cluster that extends far beyond the study area and can serve as fossil analog for nested volcanoes such as Askja, whereas in Geitafell, the sheet swarms seem to have originated from a single focus at one time, thus defining a single central volcanic complex, such as Krafla Volcano.

Sheet intrusions (inclined sheets and dykes) in the deeply eroded volcanoes of Geitafell and Dyrfjöll,eastern Iceland, were studied at the surface to identify the location, depth, and size of their magmaticsource(s). For this purpose, the measured orientations of inclined sheets were projected in three dimensionsto produce models of sheet swarm geometries. For the Geitafell Volcano, the majority of sheetsconverge toward a common focal area with a diameter of at least 4 to 7 km, the location of which coincideswith several gabbro bodies exposed at the surface. Assuming that these gabbros represent part of the magmachamber feeding the inclined sheets, a source depth of 2 to 4 km below the paleoland surface is derived.A second, younger swarm of steeply dipping sheets crosscuts this gabbro and members of the first swarm.The source of this second swarm is estimated to be located to the SE of the source of Swarm 1, below thepresent‐day level of exposure and deeper than the source of the first swarm. For the Dyrfjöll Volcano,we show that the sheets can be divided into seven different subsets, three of which can be interpretedas swarms. The most prominent swarm, the Njardvik Sheet Swarm, converges toward a several kilometerswide area in the Njardvik Valley at a depth of 1.5 to 4 km below the paleoland surface. Two additionalmagmatic sources are postulated to be located to the northeast and southwest of the main source. Crosscuttingrelationships indicate contemporaneous, as well as successive activity of different magma chambers,but without a resolvable spatial trend. The Dyrfjöll Volcano is thus part of a complex volcanic cluster thatextends far beyond the study area and can serve as fossil analog for nested volcanoes such as Askja, whereasin Geitafell, the sheet swarms seem to have originated from a single focus at one time, thus defining a singlecentral volcanic complex, such as Krafla Volcano.

Marine sediments host large amounts of methane (CH4), which is a potent greenhouse gas. Quantitative estimates for methane release from marine sediments are scarce, and a poorly constrained temporal variability leads to large uncertainties in methane emission scenarios. Here, we use 2-D and 3-D seismic reflection, multibeam bathymetric, geochemical, and sedimentological data to (I) map and describe pockmarks in the Witch Ground Basin (central North Sea), (II) characterize associated sedimentological and fluid migration structures, and (III) analyze the related methane release. More than 1,500 pockmarks of two distinct morphological classes spread over an area of 225 km(2). The two classes form independently from another and are corresponding to at least two different sources of fluids. Class 1 pockmarks are large in size (> 6 m deep, > 250 m long, and > 75 m wide), show active venting, and are located above vertical fluid conduits that hydraulically connect the seafloor with deep methane sources. Class 2 pockmarks, which comprise 99.5% of all pockmarks, are smaller (0.9-3.1 m deep, 26-140 m long, and 14-57 m wide) and are limited to the soft, fine-grained sediments of the Witch Ground Formation and possibly sourced by compaction-related dewatering. Buried pockmarks within the Witch Ground Formation document distinct phases of pockmark formation, likely triggered by external forces related to environmental changes after deglaciation. Thus, greenhouse gas emissions from pockmark fields cannot be based on pockmark numbers and present-day fluxes but require an analysis of the pockmark forming processes through geological time. Plain Language Summary Marine sediments host large amounts of methane (CH4), which is a potent greenhouse gas. The amount of methane released into the atmosphere is, however, largely unknown making it difficult to implement this methane source in climate models. Here we use geophysical, geochemical, and sedimentological data to map the distribution of fluid escape structures in the central North Sea. More than 1,500 pockmarks, which are circular to semicircular depressions of the seafloor, indicate fluid flow from the subsurface. There are two distinct morphological classes of pockmarks corresponding to at least two different fluid sources. Class 1 pockmarks are large, show active venting, and are located above vertical fluid conduits in the subsurface, which feed fluids from deeper strata. Class 2 pockmarks, which comprise 99.5% of all pockmarks, are smaller and limited to the soft sediments directly below the seafloor. Older pockmarks in the subsurface document distinct phases of pockmark formation, likely triggered by external forces after the retreat of ice in the North Sea. The amount of methane released from natural geological sources based on pockmark numbers may be wrong as these do not take into account the origin and composition of released fluids.

Marine sediments host large amounts of methane (CH4), which is a potent greenhouse gas. Quantitative estimates for methane release from marine sediments are scarce, and a poorly constrained temporal variability leads to large uncertainties in methane emission scenarios. Here, we use 2-D and 3-D seismic reflection, multibeam bathymetric, geochemical, and sedimentological data to (I) map and describe pockmarks in the Witch Ground Basin (central North Sea), (II) characterize associated sedimentological and fluid migration structures, and (III) analyze the related methane release. More than 1,500 pockmarks of two distinct morphological classes spread over an area of 225 km2. The two classes form independently from another and are corresponding to at least two different sources of fluids. Class 1 pockmarks are large in size (>6 m deep, >250 m long, and >75 m wide), show active venting, and are located above vertical fluid conduits that hydraulically connect the seafloor with deep methane sources. Class 2 pockmarks, which comprise 99.5% of all pockmarks, are smaller (0.9–3.1 m deep, 26–140 m long, and 14–57 m wide) and are limited to the soft, fine-grained sediments of the Witch Ground Formation and possibly sourced by compaction-related dewatering. Buried pockmarks within the Witch Ground Formation document distinct phases of pockmark formation, likely triggered by external forces related to environmental changes after deglaciation. Thus, greenhouse gas emissions from pockmark fields cannot be based on pockmark numbers and present-day fluxes but require an analysis of the pockmark forming processes through geological time.

Mangaia hosts the most radiogenic Pb-isotopic compositions observed in ocean island basalts and represents the HIMU (high µ = ²³⁸U/²⁰⁴Pb) mantle end-member, thought to result from recycled oceanic crust. Complete geochemical characterization of the HIMU mantle end-member has been inhibited due to a lack of deep submarine glass samples from HIMU localities. We homogenized olivine-hosted melt inclusions separated from Mangaia lavas and the resulting glassy inclusions made possible the first volatile abundances to be obtained from the HIMU mantle end-member. We also report major and trace element abundances and Pb-isotopic ratios on the inclusions, which have HIMU isotopic fingerprints. We evaluate the samples for processes that could modify the volatile and trace element abundances postmantle melting, including diffusive Fe and H₂O loss, degassing, and assimilation. H₂O/Ce ratios vary from 119 to 245 in the most pristine Mangaia inclusions; excluding an inclusion that shows evidence for assimilation, the primary magmatic H₂O/Ce ratios vary up to ∼200, and are consistent with significant dehydration of oceanic crust during subduction and long-term storage in the mantle. CO₂ concentrations range up to 2346 ppm CO₂ in the inclusions. Relatively high CO₂ in the inclusions, combined with previous observations of carbonate blebs in other Mangaia melt inclusions, highlight the importance of CO₂ for the generation of the HIMU mantle. F/Nd ratios in the inclusions (30 ± 9; 2σ standard deviation) are higher than the canonical ratio observed in oceanic lavas, and Cl/K ratios (0.079 ± 0.028) fall in the range of pristine mantle (0.02–0.08).

The 2021 La Palma eruption (Tajogaite) was unprecedented in magnitude, duration, and degree of monitoring compared to historical volcanism on La Palma. Here, we provide data on melt inclusions in samples from the beginning and end of the eruption to compare the utility of both melt and fluid inclusions as recorders of magma storage. We also investigated compositional heterogeneities within the magmatic plumbing system. We found two populations of olivine crystals: a low Mg# (78–82) population present at the beginning and end of eruption, recording the maximum volatile contents (2.5 wt % H₂O, 1,800 ppm F, 700 ppm Cl, 3,800 ppm S) and a higher Mg# (83–86) population sampled toward the end of the eruption, with lower volatile contents. Despite their host composition, melt inclusions share the same maximum range of CO₂ concentrations (1.2–1.4 wt %), indicating olivine growth and inclusion capture at similar depths. Overall, both melt and fluid inclusions record similar pressures (450–850 MPa, ∼15–30 km), and when hosted in the same olivine crystal pressures are indistinguishable within error. At these mantle pressures, CO₂ is expected to be an exsolved phase explaining the similar range of CO₂ between the two samples, but other volatile species (F, Cl, S) behave incompatibly, and thus, the increase between the two olivine populations can be explained by fractional crystallization prior to eruption. Finally, based on our new data, we provide estimates on the total volatile emission of the eruption.

The Cretaceous High Arctic Large Igneous Province (HALIP) in Canada involved extrusion of continental flood basalts (CFBs) at 130-120 Ma and 100-95 Ma and emplacement of an extensive sill and dike network that intersected the Carboniferous to Paleogene Sverdrup Basin. In this paper, we present new Ar-40/Ar-39 ages, major and trace elements, and Sr-Nd-Pb isotope ratios for HALIP lava, dikes, and sills from Bukken Fiord, NW Axel Heiberg Island, Canadian Arctic Islands. Our best constrained Ar-40/(39) ages yield a weighted average of 124.1 +/- 1 (2 sigma) Ma, coincident with the first pulse of tholeiitic CFB magmatism in the Arctic-wide HALIP as exemplified by Isachsen Formation flood basalts on Axel Heiberg Island. The Bukken Fiord samples are plagioclase and clinopyroxene-phyric tholeiitic basalts, are relatively evolved (3.2-6.5 wt% MgO), and share similar major and trace element compositions to typical HALIP tholeiites. Initial Nd-143/Nd-144 ranges from 0.51260 to 0.51291 and initial Sr-87/Sr-86 ranges from 0.70362 to 0.70776, while measured Pb-206/Pb-204, Pb-207/Pb-204, and Pb-208/Pb-204 range from 18.614 to 19.199, 15.534 to 15.630, and 38.404 to 39.054, respectively. The most primitive sample in this study has Sr-Nd-Pb isotope signatures that suggest an enriched plume-derived mantle source for HALIP tholeiites. Most samples, however, possess relatively radiogenic isotope signatures that can be explained by moderate degrees of assimilation of Sverdrup Basin sedimentary rocks. Magma-crust interaction in the HALIP plumbing system was likely widespread and may have increased the environmental impact of the HALIP, particularly if crustal carbon was volatilized.

The Cretaceous High Arctic Large Igneous Province (HALIP) in Canada involved extrusion of continental flood basalts (CFBs) at 130-120 Ma and 100-95 Ma and emplacement of an extensive sill and dike network that intersected the Carboniferous to Paleogene Sverdrup Basin. In this paper, we present new ⁴⁰Ar/³⁹Ar ages, major and trace elements, and Sr-Nd-Pb isotope ratios for HALIP lava, dikes, and sills from Bukken Fiord, NW Axel Heiberg Island, Canadian Arctic Islands. Our best constrained 40Ar/39Ar ages yield a weighted average of 124.1 ± 1 (2σ) Ma, coincident with the first pulse of tholeiitic CFB magmatism in the Arctic-wide HALIP as exemplified by Isachsen Formation flood basalts on Axel Heiberg Island. The Bukken Fiord samples are plagioclase and clinopyroxene-phyric tholeiitic basalts, are relatively evolved (3.2-6.5 wt% MgO), and share similar major and trace element compositions to typical HALIP tholeiites. Initial ¹⁴³Nd/¹⁴⁴Nd ranges from 0.51260 to 0.51291 and initial ⁸⁷Sr/⁸⁶Sr ranges from 0.70362 to 0.70776, while measured ²⁰⁶Pb/²⁰⁴Pb, ²⁰⁷Pb/²⁰⁴Pb, and ²⁰⁸Pb/²⁰⁴Pb range from 18.614 to 19.199, 15.534 to 15.630, and 38.404 to 39.054, respectively. The most primitive sample in this study has Sr-Nd-Pb isotope signatures that suggest an enriched plume-derived mantle source for HALIP tholeiites. Most samples, however, possess relatively radiogenic isotope signatures that can be explained by moderate degrees of assimilation of Sverdrup Basin sedimentary rocks. Magma-crust interaction in the HALIP plumbing system was likely widespread and may have increased the environmental impact of the HALIP, particularly if crustal carbon was volatilized.

Recent studies have shown that biosignatures of ancient microbial life exist in mineral coatings in deep bedrock fractures of Precambrian cratons, but such surveys have been few and far between. Here, we report results from southwestern Sweden in an area of 1.6–1.5 Ga Paleoproterozoic rocks heavily reworked by the 1.14–0.96 Ga Sveconorwegian orogeny, a terrane previously scarcely explored for ancient microbial biosignatures. Calcite-pyrite-adularia-illite-coated fractures were analyzed for stable isotopes via Secondary Ion Mass Spectrometry (δ13C, δ18O, δ34S) and in situ Rb/Sr geochronology via Laser-ablation inductively coupled plasma mass spectrometry. The Rb/Sr ages for calcite-adularia and calcite-illite show that several fluid flow events can be discerned (797 ± 18–769 ± 7, 391 ± 5–387 ± 6, 356 ± 5–347 ± 4, and 301 ± 7 Ma). The δ13C, δ18O and 87Sr/86Sr values of different calcite growth zones further confirmed episodic fluid flow. Pyrite δ34S values down to −49.9‰V-CDT, together with systematically increased δ34S from crystal core to rim, suggest formation following microbial sulfate reduction under semi-closed conditions. Assemblages involving MSR-related pyrite generally have Devonian to Permian Rb/Sr ages, indicating an association to extension-related fracturing and fluid mixing during foreland-basin formation linked to Caledonian orogeny in the northwest. An assemblage with an age of 301 ± 7 Ma is potentially related to Oslo Rift extension, whereas the Neo-Proterozoic ages relate to post-Sveconorwegian extensional tectonics. Remnants of short-chained fatty acids in the youngest calcite coatings further indicate a biogenic origin, while the absence of organic molecules in older calcite is in line with thermal degradation, potentially related to heating during Caledonian foreland basin burial.

Recent studies have shown that biosignatures of ancient microbial life exist in mineral coatings in deep bedrock fractures of Precambrian cratons, but such surveys have been few and far between. Here, we report results from southwestern Sweden in an area of 1.6–1.5 Ga Paleoproterozoic rocks heavily reworked by the 1.14–0.96 Ga Sveconorwegian orogeny, a terrane previously scarcely explored for ancient microbial biosignatures. Calcite-pyrite-adularia-illite-coated fractures were analyzed for stable isotopes via Secondary Ion Mass Spectrometry (δ13C, δ18O, δ34S) and in situ Rb/Sr geochronology via Laser-ablation inductively coupled plasma mass spectrometry. The Rb/Sr ages for calcite-adularia and calcite-illite show that several fluid flow events can be discerned (797 ± 18–769 ± 7, 391 ± 5–387 ± 6, 356 ± 5–347 ± 4, and 301 ± 7 Ma). The δ13C, δ18O and 87Sr/86Sr values of different calcite growth zones further confirmed episodic fluid flow. Pyrite δ34S values down to −49.9‰V-CDT, together with systematically increased δ34S from crystal core to rim, suggest formation following microbial sulfate reduction under semi-closed conditions. Assemblages involving MSR-related pyrite generally have Devonian to Permian Rb/Sr ages, indicating an association to extension-related fracturing and fluid mixing during foreland-basin formation linked to Caledonian orogeny in the northwest. An assemblage with an age of 301 ± 7 Ma is potentially related to Oslo Rift extension, whereas the Neo-Proterozoic ages relate to post-Sveconorwegian extensional tectonics. Remnants of short-chained fatty acids in the youngest calcite coatings further indicate a biogenic origin, while the absence of organic molecules in older calcite is in line with thermal degradation, potentially related to heating during Caledonian foreland basin burial.

Continental collision causes deformation in the crust along shear zones. However, the physical and chemical conditions at which these zones operate and the deformation processes that enable up to hundreds of km of tectonic transport are still unclear because of the depth at which they occur and the challenges in imaging them. Ancient exhumed collision zones allow us to investigate these processes much better, for example at the COSC‐1 borehole in the central Scandinavian Caledonides. This study combines data from the COSC‐1 borehole with different seismic measurements to provide constraints on the spatial lithological and textural configuration of the Seve Nappe Complex. This is one of the few studies that shows that core‐log‐seismic integration in metamorphic rocks allows to identify the spatial distribution of major lithological units. Especially gamma ray logs in combination with density data are powerful tools to distinguish between mafic and felsic lithologies in log‐core correlation. Our results indicate that reflections along the borehole are primarily caused by compositional rather than textural changes. Reflections in the Seve Nappe Complex are not as distinct as in greater depths but continuous and several of them can be linked to magmatic intrusions, which have been metamorphically overprinted. Their setting indicates that the Seve Nappe Complex consists of the remnants of a volcanic continental margin. Our results suggest that ductile‐deformed middle crustal reflectivity is primarily a function of pre‐orogenic lithological variations which has to be considered when deciphering mountain building processes.

Plain Language Summary

Areas where continents collide experience different kind of deformation. However, these processes and the conditions at which they take place are difficult to study because of the great depth at which they occur. Former collision zones that are closer to the surface these days allow the investigation of these processes much better, for example at the COSC‐1 borehole in the central Scandinavian Caledonides. The challenge remains to image the remnant of these processes in high detail but at the same time over a large area. This study combines data from the COSC‐1 borehole with different geophysical measurements to better understand the lithology and structure of the Seve Nappe Complex. We show that the combination of these data allows us to distinguish between rocks from mafic and sedimentary origin. Our results indicate that the geophysical data along the borehole image the change of the composition of the rocks which probably originates from magmatic intrusions and have been overprinted by geological processes, rather than from fracture zones.

The 2014-2015 Holuhraun eruption on Iceland was located within the Askja fissure swarm but was accompanied by caldera subsidence in the Baroarbunga central volcano 45 km to the southwest. Geophysical monitoring of the eruption identified a seismic swarm that migrated from Baroarbunga to the Holuhraun eruption site over the course of two weeks. In order to better understand this lateral connection between Baroarbunga and Holuhraun, we present mineral textures and compositions, mineral-melt-equilibrium calculations, whole rock and trace element data, and oxygen isotope ratios for selected Holuhraun samples. The Holuhraun lavas are compositionally similar to recorded historical eruptions from the Baroarbunga volcanic system but are distinct from the historical eruption products of the nearby Askja system. Thermobarometry calculations indicate a polybaric magma plumbing system for the Holuhraun eruption, wherein clinopyroxene and plagioclase crystallized at average depths of approximate to 17 km and approximate to 5 km, respectively. Crystal resorption textures and oxygen isotope variations imply that this multilevel plumbing system facilitated magma mixing and assimilation of low-O-18 Icelandic crust prior to eruption. In conjunction with the existing geophysical evidence for lateral migration, our results support a model of initial vertical magma ascent within the Baroarbunga plumbing system followed by lateral transport of aggregated magma batches within the upper crust to the Holuhraun eruption site.

The 2014–2015 Holuhraun eruption on Iceland was located within the Askja fissure swarm butwas accompanied by caldera subsidence in the Barðarbunga central volcano 45 km to the southwest. Geophysicalmonitoring of the eruption identified a seismic swarm that migrated from Barðarbunga to theHoluhraun eruption site over the course of two weeks. In order to better understand this lateral connectionbetween Barðarbunga and Holuhraun, we present mineral textures and compositions, mineral-meltequilibriumcalculations, whole rock and trace element data, and oxygen isotope ratios for selected Holuhraunsamples. The Holuhraun lavas are compositionally similar to recorded historical eruptions from theBarðarbunga volcanic system but are distinct from the historical eruption products of the nearby Askja system.Thermobarometry calculations indicate a polybaric magma plumbing system for the Holuhraun eruption,wherein clinopyroxene and plagioclase crystallized at average depths of 17 km and 5 km,respectively. Crystal resorption textures and oxygen isotope variations imply that this multilevel plumbingsystem facilitated magma mixing and assimilation of low-d18O Icelandic crust prior to eruption. In conjunctionwith the existing geophysical evidence for lateral migration, our results support a model of initial verticalmagma ascent within the Barðarbunga plumbing system followed by lateral transport of aggregatedmagma batches within the upper crust to the Holuhraun eruption site.

Analysis of sediment samples in the visible to mid infrared (IR) region requires small amounts of sample material and enables rapid and cost efficient geochemical analysis of mineral and organic sediment components. Here we use geochemical properties (total organic and inorganic carbon, biogenic silica, total nitrogen) from the ICDP deep drilling project PASADO to compare three different IR spectroscopy techniques: Diffuse Reflectance Fourier Transform IR Spectrometry (DRIFTS), Attenuated Total Reflectance Fourier Transform IR Spectroscopy (ATR-FTIRS) and Visible Near IR Spectroscopy (VNIRS). ATR-FTIRS and VNIRS are more rapid techniques compared to DRIFTS. Results show that calibration models developed using DRIFTS are most robust (correlation coefficient: R = 0.92 for TIC, R = 0.84 for BSi, R = 0.97 for TOC, R = 0.95 for TN). However, good statistical performance was also obtained by using ATR-FTIRS and VNIRS. When time and costs are limiting factors, these tools may be given preference for rapid biogeochemical screening.

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, microfluoresence, micro‐computed tomography, and scanning transmission electron microscopy coupled with Raman spectroscopy to create finely resolved spatial frameworks that record a continuum of alteration in glass and olivine. Microanalytical maps of vesicular and fractured lapilli in specimens from 157.1‐, 137.9‐, and 102.6‐m depths 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 subcircular 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.

[1] This paper presents new results from Optically Stimulated Luminescence (OSL) dating on a sediment core raised from the crest of the Lomonosov Ridge in the central Arctic Ocean. There has been much debate about dating sediment cores from the central Arctic Ocean and by using an independent absolute dating technique we aim to test whether or not relatively fast, cm-scale/ka, sedimentation rates were typical of Arctic’s Pleistocene depositional mode. On the basis of mainly paleomagnetic reversal stratigraphy, many previous studies suggest mm-scale/ka sedimentation rates. A common feature in these studies is that the first down core paleomagnetic negative inclination is consistently interpreted as the Brunhes/Matuyama boundary at about 780 ka. Our OSL dating results indicate that this assumption is not generally valid, and that the first encountered negative inclination represents younger age excursions within the Brunhes Chron, implying reinterpretation of many published core studies where paleoenvironmental reconstructions have been made for the central Arctic Ocean. Our dating results furthermore corroborates a correlation of the uppermost 2-3 m of the Lomonosov Ridge cores to a well-dated core located off the Barents-Kara Sea margin that in turn is correlated to cores in the Fram Strait. Valuable information on the paleoceanographical evolution in the Arctic Ocean from MIS 6 to the Holocene is given through this correlation of records from the central Arctic Ocean to records off the Eurasian continental margin.

[1] This paper presents an analysis of the Arctic Ocean and its constituent seas for seafloor area distribution versus depth and ocean volume. The bathymetry from the International Bathymetric Chart of the Arctic Ocean (IBCAO) is used together with limits defining this ocean and its constituent seas from the International Hydrographic Organization (IHO) as well as redefined limits constructed to confine the seas to the shallow shelves. IBCAO is a bathymetric grid model with a resolution of 2.5 2.5 km, which significantly improved the portrayal of the Arctic Ocean seafloor through incorporation of newly released bathymetric data including echo soundings from U.S. and British navies, scientific nuclear submarine cruises, and icebreaker cruises. This analysis of seafloor area and ocean volume is the first for the Arctic Ocean based on this new and improved portrayal of the seafloor as represented by IBCAO. The seafloor area and volume are calculated for different depths starting from the present sea level and progressing in increments of 10 m to a depth of 500 m and in increments of 50 m from 550 m down to the deepest depth within each of the analyzed seas. Hypsometric curves expressed as simple histograms of the frequencies in different depth bins and depth plotted against cumulative area for each of the analyzed seas are presented. The area and volume calculations show that the entire IHO-defined Arctic Ocean makes up 4.3% of the total ocean area but only 1.4% of the volume. Furthermore, the IHO Arctic Ocean is the shallowest (mean depth 1201 m) of all the major oceans and their adjacent seas. The continental shelf area, from the coasts out to the shelf break, make up as much as 52.9% of the total area in the Arctic Ocean, defined in this work as consisting of the oceanic deep Arctic Ocean Basin; the broad continental shelves of the Barents, Kara, Laptev, East Siberian, Chukchi, and Beaufort Seas; the White Sea; and the narrow continental shelf off both the Canadian Arctic Archipelago and northern Greenland. This result indicates that the Arctic Ocean has significantly larger continental shelves compared with all the other oceans, where previous studies show that the proportion of shelves, from the coasts out to the foot of the continental slopes, only ranges between about 9.1 and 17.7%. Furthermore, the derived hypsometric curves show that most of the Arctic Ocean shelf seas besides the Barents Sea, Beaufort Sea, and the shelf off northern Greenland have a similar shape, with the largest seafloor area between 0 and 50 m. The East Siberian and Laptev seas, in particular, show area distributions concentrated in this shallow depth range, and together with the Chukchi Sea they form a large flat shallow shelf province composing as much as 22% of the entire Arctic Ocean area but only 1% of the volume. This implies that the circulation in the Arctic Ocean might be very sensitive to eustatic sea level changes. One of the aims with this work is to make up-to-date high-resolution area and volume calculations for the Arctic Ocean at various depths available for download.

The Macquarie Ridge Complex (MRC) on the Australia‐Pacific plate boundary south of New Zealand is an extinct mid‐ocean ridge that has experienced a complex tectonic history and produced highly heterogeneous mid‐ocean ridge basalts (MORBs). When and how seafloor spreading ceased along the proto‐Macquarie mid‐ocean ridge remain elusive, and it is unclear how the mantle source of MORBs is affected by the gradual cessation of seafloor spreading at mid‐ocean ridges. To constrain the tectonic evolution of the MRC, the mantle source variations for MORBs at dying mid‐ocean ridges, and the mechanisms of mantle enrichment and asthenospheric heterogeneities, we report 11 pyroxene, plagioclase, basaltic glass, groundmass, and sericite ⁴⁰Ar/³⁹Ar and one zircon U‐Pb ages for the MRC MORBs. Our data reveal that basalts from the MRC seamounts were erupted between 25.9 and 1.6 Ma and Macquarie Island at ∼10 Ma. Combined age and plate reconstruction results reveal that the cessation of seafloor spreading at the MRC generally propagated from south to north along the ridge. Basalts produced by the then dying Macquarie mid‐ocean ridge at different times on different seamounts/island show a large variation in isotopic compositions and there is no clear correlation between ages and isotopic ratios. The heterogeneity of mantle source for MORBs from the proto‐Macquarie mid‐ocean ridge suggests that the upper asthenospheric mantle is heterogeneous, and such heterogeneity becomes most obvious at dying mid‐ocean ridges where the degrees of partial melting are low and a large range of melt compositions are produced.

The response of Arctic Ocean biogeochemistry to subsurface flow driven by permafrost thaw is poorly understood. We present dissolved chloride and water isotopic data from the Chukchi Sea Shelf sediments that reveal the presence of a meteoric subsurface flow enriched in cations with a radiogenic Sr fingerprint. This subsurface fluid is also enriched in dissolved inorganic carbon and methane that bear isotopic compositions indicative of a carbon reservoir modified by reactions in a closed system. Such fluid characteristics are in stark contrast with those from other sites in the Chukchi Sea where the pore water composition shows no sign of meteoric input, but reflect typical biogeochemical reactions associated with early diagenetic sequences in marine sediment. The most likely source of the observed subsurface flow at the Chukchi Sea Shelf is from the degradation of permafrost that had extended to the shelf region during the Last Glacial Maximum. Our data suggest that the permafrost-driven subsurface flow most likely took place during the 2-3 degrees C warming in the Early Holocene Thermal Maximum. This time scale is supported by numerical simulation of pore water profiles, which indicate that a minimum of several thousand years must have passed since the cessation of the subsurface methane-bearing fluid flow.

The Himalaya and the Tibetan Plateau are uplifted by the ongoing northward underthrusting of the Indian continental lithosphere below Tibet resulting in lithospheric stacking. The layered structure of the Tibetan upper mantle is imaged by seismic methods, most detailed with the receiver function method. Tibet is considered as a place where the development of a future craton is currently under way. Here we study the upper mantle from Germany to northern Sweden with seismic S receiver functions and compare the structure below Scandinavia with that below Tibet. Below Proterozoic Scandinavia, we found two low-velocity zones on top of each other, separated by a high-velocity zone. The top of the upper low-velocity zone at about 100 km depth extends from Germany to Archaean northern Sweden. It agrees with the lithosphere-asthenosphere boundary (LAB) below Germany and Denmark. Below Sweden it is known as the 8 degrees discontinuity, or as a mid-lithospheric discontinuity (MLD), similar to observations in North America. Seismic tomography places the LAB near 200 km in Scandinavia, which is close to the top of our deeper low-velocity zone. We also observed the bottom of the asthenosphere (the Lehmann discontinuity) deepening from 180 km in Germany to 260 km below Sweden. Remnants of old subduction in the upper about 100 km below Scandinavia and Finland are known from controlled source seismic experiments and local earthquake studies. Recent tomographic studies indicate delamination of the lithosphere below southern Scandinavia and northern Germany. We are suggesting that the large-scale layered structure in the Scandinavian upper mantle may be caused by processes similar to the ongoing lithospheric stacking in Tibet.

Day‐Dunlop plots are widely used in paleomagnetic and environmental studies as a tool to determine the magnetic domain state of magnetite, i.e., superparamagnetic (SP), stable single‐domain (SD), pseudosingle‐domain (PSD), multidomain (MD), and their mixtures. The few experimental studies that have examined hysteresis properties of SD‐SP mixtures of magnetite found that the ratios of saturation remanent magnetization to saturation magnetization and the coercivity of remanence to coercivity are low, when compared to expected theoretical mixing trends based on Langevin theory. This study reexamines Day‐Dunlop plots using experimentally controlled mixtures of SD and SP magnetite grains. End‐members include magnetotactic bacteria (MSR‐1) as the SD source, and a commercial ferrofluid or magnetotactic bacteria (ΔA12) as the SP source. Each SP‐component was added incrementally to a SD sample. Experimental results from these mixing series show that the magnetization and coercivity ratios are lower than the theoretical prediction for bulk SP magnetic size. Although steric repulsion was present between the particles, we cannot rule out interaction in the ferrofluid for higher concentrations. The SP bacteria are noninteracting as the magnetite was enclosed by an organic bilipid membrane. Our results demonstrate that the magnetization and coercivity ratios of SD‐SP mixtures can lie in the PSD range, and that an unambiguous interpretation of particle size can only be made with information about the magnetic properties of the end‐members.

Large-scale, catastrophic mass wasting is a major process contributing to the dismantling of oceanic intraplate volcanoes. Recent studies, however, have highlighted a possible feedback relationship between flank collapse, or incipient instability, and subsequent episodes of structural rearrangement and/or renewed volcano growth. The Teno massif, located in northwestern Tenerife (Canary Islands), is a deeply eroded Miocene shield volcano that was built in four major eruptive phases punctuated by two lateral collapses, each removing >20–25 km³ of the volcano's north flank. In this paper, we use detailed field observations and petrological and geochemical data to evaluate possible links between large-scale landslides and subsequent volcanism/magmatism during Teno's evolution. Inspection of key stratigraphic sequences reveals that steep angular unconformities, relics of paleolandslide scars, are marked by polymict breccias. Near their base, these deposits typically include abundant juvenile pyroclastic material, otherwise scarce in the region. While some of Teno's most evolved, low-density magmas were produced just before flank collapses, early postlandslide lava sequences are characterized by anomalously high proportions of dense ankaramite flows, extremely rich in clinopyroxene and olivine crystals. A detailed sampling profile shows transitions from low-Mg # lavas relatively rich in SiO₂ to lavas with low silica content and comparatively high Mg # after both landslides. Long-term variations in Zr/Nb, normative nepheline, and La/Lu are coupled but do not show a systematic correlation with stratigraphic boundaries. We propose that whereas loading of the growing precollapse volcano promoted magma stagnation and differentiation, the successive giant landslides modified the shallow volcano-tectonic stress field at Teno, resulting in widespread pyroclastic eruptions and shallow magma reservoir drainage. This rapid unloading of several tens of km³ of near-surface rocks appears to have upset magma differentiation processes, while facilitating the remobilization and tapping of denser ankaramite magmas that were stored in the uppermost mantle. Degrees of mantle melting coincidently reached a maximum in the short time interval between the two landslides and declined shortly after, probably reflecting intrinsic plume processes rather than a collapse-induced influence on mantle melting. Our study of Teno volcano bears implications for other oceanic volcanoes where short-term compositional variations may also directly relate to major flank collapse events.

Deciphering the polymorphic nature and the degree of iron lattice-preferred orientation in the Earth's inner core holds a key to understanding the present status and evolution of the inner core. A multiphase lattice-preferred orientation pattern is obtained for the top 350 km of the inner core by means of the ab initio based Candy Wrapper Velocity Model coupled to a Monte Carlo phase discrimination scheme. The achieved geographic distribution of lattice alignment is characterized by two regions of freezing, namely within South America and the Western Central Pacific, that exhibit an uncommon high degree of lattice orientation. In contrast, widespread regions of melting of relatively weak lattice ordering permeate the rest of the inner core. The obtained multiphase lattice-preferred orientation pattern is in line with mantle-constrained geodynamo simulations and allows to setup an ad hoc mineral physics scenario for the complex Earth's inner core. It is found that the cubic phase of iron is the dominating iron polymorph in the outermost part of the inner core.

Granitic magma bodies form in the ephemeral part of magma mush systems and are emplaced by a variety of mechanisms in different tectonic settings. This study investigates how granitic magma emplacement processes and tectonomagmatic interactions assert control over the architecture of mush state pluton-scale magma transport pathways. The 1.45 Ga shallow-crustal Gotemar pluton is a 4.5 km diameter circular pluton that consists of three granite units: a coarse-grained red granite, a medium-grained pale to red granite, and fine-grained pale microgranite sheets. We employed geological mapping supported by Anisotropy of Magnetic Susceptibility (AMS) to examine the magmatic and regional tectonic controls on late-stage magma transport in the Gotemar granitic magma mush system. Multiple parallel arcuate subhorizontal microgranite and medium-grained granite sheets (from 0.1 to 10s of meters thick) were mapped within the pluton. The arcuate sheets pinch out from the northern part of the pluton toward the SE inferring magma propagation direction. A dominant set of vertical granitic sheets within the granite body strikes NW-SE. The AMS fabrics are contact-parallel in the main medium-grained granite body and indicate inflation. Within the microgranite sheets, the AMS fabrics are parallel to the sheet strike and support a sheet propagation direction to the SE. The Gotemar pluton displays a clear link between arcuate (concentric) magma-transporting sheets and concentric strain-partitioning related to the intrusion of medium-grained granite magma. The vertical magma sheet orientations are consistent with an NE-SW extensional stress field that is associated with the extensional back-arc stress regime of the contemporary Hallandian Orogen. The eruptive products of volcanoes are thought to be stored in pockets of melt in crystal-dominated magmatic systems called crystal mushes prior to volcanic eruptions. An understanding of where magma is stored and how it is transported in mush systems is important in order to predict the eruptive behavior of the volcanic system. This contribution investigates the magma transport pathways in the Gotemar granite in Sweden and its relationship to local magmatic deformation and regional deformation related to the Hallandian mountain building event. We show that magma is transported in vertical sheets parallel to the front of the Hallandian Orogen and laterally in sub-horizontal arcuate sheets that reflect the circular shape of the granite pluton. Our study highlights the importance of understanding the shape and the formation of the magmatic granite body for deciphering the melt transport in the magma mush system under volcanoes. Subhorizontal magma transport in a granitic magma mush controlled by magma emplacement structuresVertical magma sheets in granitic magma mush controlled by the regional stress fieldMagma transported in concentric magma fingers in the circular Gotemar granite pluton

How subsurface microbial life changed at the bottom of the kilometers-deep (hypo) Critical Zone in response to evolving surface conditions over geologic time is an open question. This study investigates the burial and exhumation, biodegradation, and fluid circulation history of hydrocarbon reservoirs across the Colorado Plateau as a window into the hypo-Critical Zone. Hydrocarbon reservoirs, in the Paradox and Uinta basins, were deeply buried starting ca. 100 to 60 Ma, reaching temperatures >80–140°C, likely sterilizing microbial communities present since the deposition of sediments. High salinities associated with evaporites may have further limited microbial activity. Upward migration of hydrocarbons from shale source rocks into shallower reservoirs during maximum burial set the stage for microbial re-introduction by creating organic-rich “hot spots.” Denudation related to the incision of the Colorado River over the past few million years brought reservoirs closer to the surface under cooler temperatures, enhanced deep meteoric water circulation and flushing of saline fluids, and likely re-inoculated more permeable sediments up to several km depth. Modern- to paleo-hydrocarbon reservoirs show molecular and isotopic evidence of anaerobic oxidation of hydrocarbons coupled to bacterial sulfate reduction in areas with relatively high SO4-fluxes. Anaerobic oil biodegradation rates are high enough to explain the removal of at least some portion of postulated “supergiant oil fields” across the Colorado Plateau by microbial activity over the past several million years. Results from this study help constrain the lower limits of the hypo-Critical Zone and how it evolved over geologic time, in response to changing geologic, hydrologic, and biologic forcings.

While basaltic volcanism is dominant during rifting and continental breakup, felsic magmatism may be a significant component of some rift margins. During International Ocean Discovery Program (IODP) Expedition 396 on the continental margin of Norway, a graphite-garnet-cordierite bearing dacitic unit (the Mimir dacite) was recovered in two holes within early Eocene sediments on Mimir High (Site U1570), a marginal high on the Vøring Transform Margin. Here, we present a comprehensive textural, petrological, and geochemical study of the Mimir dacite in order to assess its origin and discuss the geodynamic implications. The major mineral phases (garnet, cordierite, quartz, plagioclase, alkali feldspar) are hosted in a fresh rhyolitic, vesicular, glassy matrix that is locally mingled with sediments. The major element chemistry of garnet and cordierite, the presence of zircon inclusions with inherited cores, and thermobarometric calculations all support an upper crustal metapelitic origin. While most magma-rich margin models favor crustal anatexis in the lower crust, thermobarometric calculations performed here show that the Mimir dacite was produced at upper-crustal depths (<5 kbar, 18 km depth) and high temperature (750–800°C) with up to 3 wt% water content. In situ U-Pb analyses on zircon inclusions give a magmatic crystallization age of 54.6 ± 1.1 Ma, consistent with emplacement that post-dates the Paleocene-Eocene Thermal Maximum. Our results suggest that the opening of the Northeast Atlantic was associated with a phase of low-pressure, high-temperature crustal anatexis preceding the main phase of magmatism.

Laminated, organic-rich silts and clays with high dissolved gas content characterize sediments at IODP Site M0063 in the Landsort Deep, which at 459 m is the deepest basin in the Baltic Sea. Cores recovered from Hole M0063A experienced significant expansion as gas was released during the recovery process, resulting in high sediment loss. Therefore, during operations at subsequent holes, penetration was reduced to 2 m per 3.3 m core, permitting expansion into 1.3 m of initially empty liner. Fully filled liners were recovered from Holes B through E, indicating that the length of recovered intervals exceeded the penetrated distance by a factor of > 1.5. A typical down-core logarithmic trend in gamma density profiles, with anomalously low-density values within the upper similar to 1 m of each core, suggests that expansion primarily occurred in this upper interval. Thus, we suggest that a simple linear correction is inappropriate. This interpretation is supported by anisotropy of magnetic susceptibility data that indicate vertical stretching in the upper similar to 1.5 m of expanded cores. Based on the mean gamma density profiles of cores from Holes M0063C and D, we obtain an expansion function that is used to adjust the depth of each core to conform to its known penetration. The variance in these profiles allows for quantification of uncertainty in the adjusted depth scale. Using a number of bulk C-14 dates, we explore how the presence of multiple carbon source pathways leads to poorly constrained radiocarbon reservoir age variability that significantly affects age and sedimentation rate calculations.

Laminated, organic-rich silts and clays with high dissolved gas content characterize sediments at IODP Site M0063 in the Landsort Deep, which at 459 m is the deepest basin in the Baltic Sea. Cores recovered from Hole M0063A experienced significant expansion as gas was released during the recovery process, resulting in high sediment loss. Therefore, during operations at subsequent holes, penetration was reduced to 2 m per 3.3 m core, permitting expansion into 1.3 m of initially empty liner. Fully filled liners were recovered from Holes B through E, indicating that the length of recovered intervals exceeded the penetrated distance by a factor of >1.5. A typical down-core logarithmic trend in gamma density profiles, with anomalously low-density values within the upper ∼1 m of each core, suggests that expansion primarily occurred in this upper interval. Thus, we suggest that a simple linear correction is inappropriate. This interpretation is supported by anisotropy of magnetic susceptibility data that indicate vertical stretching in the upper ∼1.5 m of expanded cores. Based on the mean gamma density profiles of cores from Holes M0063C and D, we obtain an expansion function that is used to adjust the depth of each core to conform to its known penetration. The variance in these profiles allows for quantification of uncertainty in the adjusted depth scale. Using a number of bulk 14C dates, we explore how the presence of multiple carbon source pathways leads to poorly constrained radiocarbon reservoir age variability that significantly affects age and sedimentation rate calculations.

Surface heat flow data in the Arctic Ocean are needed to assess hydrocarbon and methane hydrate distributions, and provide constraints into the tectonic origins and nature of underlying crust. However, across broad areas of the Arctic, few published measurements exist. This is true for the outer continental shelf and slope of the East Siberian Sea, and the adjoining deep water ridges and basins. Here we present 21 new surface heat flow measurements from this region of the Arctic Ocean. On the Southern Lomonosov Ridge, the average measured heat flow, uncorrected for effects of sedimentation and topography, is 57±4 mW/m² (n=4). On the outer continental shelf and slope of the East Siberian Sea (ESS), the average is 57±10 mW/m² (n=16). An anomalously high heat flow of 203±28 mW/m² was measured at a single station in the Herald Canyon. With the exception of this high heat flow, the new data from the ESS are consistent with predictions for thermally equilibrated lithosphere of continental origin that was last affected by thermotectonic processes in the Cretaceous to early Cenozoic. Variability within the data likely arises from differences in radiogenic heat production within the continental crust and overlying sediments. This can be further explored by comparing the data with geophysical constraints on sediment and crustal thicknesses.

Surface heat flow data in the Arctic Ocean are needed to assess hydrocarbon and methane hydrate distributions, and provide constraints into the tectonic origins and nature of underlying crust. However, across broad areas of the Arctic, few published measurements exist. This is true for the outer continental shelf and slope of the East Siberian Sea, and the adjoining deep water ridges and basins. Here we present 21 new surface heat flow measurements from this region of the Arctic Ocean. On the Southern Lomonosov Ridge, the average measured heat flow, uncorrected for effects of sedimentation and topography, is 574 mW/m(2) (n=4). On the outer continental shelf and slope of the East Siberian Sea (ESS), the average is 5710 mW/m(2) (n=16). An anomalously high heat flow of 20328 mW/m(2) was measured at a single station in the Herald Canyon. With the exception of this high heat flow, the new data from the ESS are consistent with predictions for thermally equilibrated lithosphere of continental origin that was last affected by thermotectonic processes in the Cretaceous to early Cenozoic. Variability within the data likely arises from differences in radiogenic heat production within the continental crust and overlying sediments. This can be further explored by comparing the data with geophysical constraints on sediment and crustal thicknesses.

We use high-definition seafloor digital imagery and multibeam bathymetric data acquired during the 2007 Arctic Gakkel Vents Expedition (AGAVE) to evaluate the volcanic characteristics of the 85 degrees E segment of the ultraslow spreading Gakkel Ridge (9 mm yr(-1) full rate). Our seafloor imagery reveals that the axial valley is covered by numerous, small-volume (order similar to 1000 m(3)) lava flows displaying a range of ages and morphologies as well as unconsolidated volcaniclastic deposits with thicknesses up to 10 cm. The valley floor contains two prominent volcanic lineaments made up of axis-parallel ridges and small, cratered volcanic cones. The lava flows appear to have erupted from a number of distinct source vents within the similar to 12-15 km-wide axial valley. Only a few of these flows are fresh enough to have potentially erupted during the 1999 seismic swarm at this site, and these are associated with the Oden and Loke volcanic cones. We model the widespread volcaniclastic deposits we observed on the seafloor as having been generated by the explosive discharge of CO2 that accumulated in (possibly deep) crustal melt reservoirs. The energy released during explosive discharge, combined with the buoyant rise of hot fluid, lofted fragmented clasts of rapidly cooling magma into the water column, and they subsequently settled onto the seafloor as fall deposits surrounding the source vent.

We use high-definition seafloor digital imagery and multibeam bathymetric data acquired during the 2007 Arctic Gakkel Vents Expedition (AGAVE) to evaluate the volcanic characteristics of the 85 degrees E segment of the ultraslow spreading Gakkel Ridge (9 mm yr(-1) full rate). Our seafloor imagery reveals that the axial valley is covered by numerous, small-volume (order similar to 1000 m(3)) lava flows displaying a range of ages and morphologies as well as unconsolidated volcaniclastic deposits with thicknesses up to 10 cm. The valley floor contains two prominent volcanic lineaments made up of axis-parallel ridges and small, cratered volcanic cones. The lava flows appear to have erupted from a number of distinct source vents within the similar to 12-15 km-wide axial valley. Only a few of these flows are fresh enough to have potentially erupted during the 1999 seismic swarm at this site, and these are associated with the Oden and Loke volcanic cones. We model the widespread volcaniclastic deposits we observed on the seafloor as having been generated by the explosive discharge of CO2 that accumulated in (possibly deep) crustal melt reservoirs. The energy released during explosive discharge, combined with the buoyant rise of hot fluid, lofted fragmented clasts of rapidly cooling magma into the water column, and they subsequently settled onto the seafloor as fall deposits surrounding the source vent.

Computed tomography (CT) of sediment cores allows for high-resolution images, three-dimensional volumes, and down core profiles. These quantitative data are generated through the attenuation of X-rays, which are sensitive to sediment density and atomic number, and are stored in pixels as relative gray scale values or Hounsfield units (HU). We present a suite of MATLAB^TM tools specifically designed for routine sediment core analysis as a means to standardize and better quantify the products of CT data collected on medical CT scanners. SedCT uses a graphical interface to process Digital Imaging and Communications in Medicine (DICOM) files, stitch overlapping scanned intervals, and create down core HU profiles in a manner robust to normal coring imperfections. Utilizing a random sampling technique, SedCT reduces data size and allows for quick processing on typical laptop computers. SedCTimage uses a graphical interface to create quality tiff files of CT slices that are scaled to a user-defined HU range, preserving the quantitative nature of CT images and easily allowing for comparison between sediment cores with different HU means and variance. These tools are presented along with examples from lacustrine and marine sediment cores to highlight the robustness and quantitative nature of this method.

Although it is widely accepted that shallow silicic magma reservoirs exist, and can feed eruptions, their dynamics and longevity are a topic of debate. Here, we use field mapping, geochemistry, 3D pluton reconstruction and a thermal model to investigate the assembly and eruptive history of the shallow Reyðarártindur Pluton, southeast Iceland. Primarily, the exposed pluton is constructed of a single rock unit, the Main Granite (69.9–77.7 wt.% SiO₂). Two further units are locally exposed as enclaves at the base of the exposure, the Granite Enclaves (67.4–70.2 wt.% SiO₂), and the Quartz Monzonite Enclaves (61.8–67.3 wt.% SiO₂). Geochemically, the units are related and were likely derived from the same source reservoir. In 3D, the pluton has a shape characterized by flat roof segments that are vertically offset and a volume of >2.5 km³. The pluton roof is intruded by dikes from the pluton, and in two locations displays depressions associated with large dikes. Within these particular dikes the rock is partially to wholly tuffisitic, and rock compositions range from quartz monzonite to granite. We interpret these zones as eruption-feeding conduits from the pluton. A lack of cooling contacts throughout the pluton indicates rapid magma emplacement and a thermal model calculates the top 75 m would have rheologically locked up within 1,000 years. Hence, we argue that the Reyðarártindur Pluton was an ephemeral part of the wider plumbing system that feeds a volcano, and that timeframes from emplacement to eruption were rapid.

Although it is widely accepted that shallow silicic magma reservoirs exist, and can feed eruptions, their dynamics and longevity are a topic of debate. Here, we use field mapping, geochemistry, 3D pluton reconstruction and a thermal model to investigate the assembly and eruptive history of the shallow Reyoarartindur Pluton, southeast Iceland. Primarily, the exposed pluton is constructed of a single rock unit, the Main Granite (69.9-77.7 wt.% SiO2). Two further units are locally exposed as enclaves at the base of the exposure, the Granite Enclaves (67.4-70.2 wt.% SiO2), and the Quartz Monzonite Enclaves (61.8-67.3 wt.% SiO2). Geochemically, the units are related and were likely derived from the same source reservoir. In 3D, the pluton has a shape characterized by flat roof segments that are vertically offset and a volume of >2.5 km(3). The pluton roof is intruded by dikes from the pluton, and in two locations displays depressions associated with large dikes. Within these particular dikes the rock is partially to wholly tuffisitic, and rock compositions range from quartz monzonite to granite. We interpret these zones as eruption-feeding conduits from the pluton. A lack of cooling contacts throughout the pluton indicates rapid magma emplacement and a thermal model calculates the top 75 m would have rheologically locked up within 1,000 years. Hence, we argue that the Reyoarartindur Pluton was an ephemeral part of the wider plumbing system that feeds a volcano, and that timeframes from emplacement to eruption were rapid.

The Atlantic-Mediterranean exchange of water at Gibraltar represents a significant heat and freshwater sink for the North Atlantic and is a major control on the heat, salt and freshwater budgets of the Mediterranean Sea. Consequently, an understanding of the response of the exchange system to external changes is vital to a full comprehension of the hydrographic responses in both ocean basins. Here, we use a synthesis of empirical (oxygen isotope, planktonic foraminiferal assemblage) and modeling (analytical and general circulation) approaches to investigate the response of the Gibraltar Exchange system to Atlantic freshening during Heinrich Stadials (HSs). HSs display relatively flat W–E surface hydrographic gradients more comparable to the Late Holocene than the Last Glacial Maximum. This is significant, as it implies a similar state of surface circulation during these periods and a different state during the Last Glacial Maximum. During HS1, the gradient may have collapsed altogether, implying very strong water column stratification and a single thermal and δ¹⁸O_water condition in surface water extending from southern Portugal to the eastern Alboran Sea. Together, these observations imply that inflow of Atlantic water into the Mediterranean was significantly increased during HS periods compared to background glacial conditions. Modeling efforts confirm that this is a predictable consequence of freshening North Atlantic surface water with iceberg meltwater and indicate that the enhanced exchange condition would last until the cessation of anomalous freshwater supply into to the northern North Atlantic. The close coupling of dynamics at Gibraltar Exchange with the Atlantic freshwater system provides an explanation for observations of increased Mediterranean Outflow activity during HS periods and also during the last deglaciation. This coupling is also significant to global ocean dynamics, as it causes density enhancement of the Atlantic water column via the Gibraltar Exchange to be inversely related to North Atlantic surface salinity. Consequently, Mediterranean enhancement of the Atlantic Meridional Overturning Circulation will be greatest when the overturning itself is at its weakest, a potentially critical negative feedback to Atlantic buoyancy change during times of ice sheet collapse.