Search

Boninite samples from several intervals within Hole U1439C, recovered during IODP Expedition 352, show highly variable mineral chemistries that imply complex crystallization histories. Small pyroxene grains show oscillatory zoning with cores and zones ranging from pigeonite to augite. Late crystallizing augite has highly variable Al2O3 contents (1.9–13.7 wt%) and Ca-Tschermak component contents (3–13 mol%), which reflect disequilibrium conditions. Large, euhedral, low-Ca pyroxene (i.e., enstatite/clinoenstatite) crystals exhibit complex sector and oscillatory zoning patterns. Cr-rich spinel is found as inclusions both in olivine and low-Ca pyroxene. Early crystallized olivine phenocrysts have embayed and reacted margins, and some early crystallized olivines exhibit zoning. A few olivines have multiple zones, with both normal and reverse zoning between Fo86–92. Olivine xenocrysts also have embayed and reacted margins; however, xenocrysts do not exhibit chemical zoning patterns and have consistent Fo88 compositions. Disequilibrium crystallization of pyroxene rims reflects rapid cooling during an eruption. Sector zoning in low-Ca pyroxenes is the result of crystallization during periods of moderate undercooling between mixing events. Oscillatory, normal, and reverse zoning in olivine and pyroxene appears to have formed in response to multi-stage magma mingling or mixing processes, which introduced additional Ca, Fe, Ti, and Al into parental boninitic melts. The presence of olivine xenocrysts and orthopyroxene indicate equilibrium at 2–4 kbar (Whattam et al. 2020) indicates that boninite magma mixing events likely occurred within shallow magma chambers containing olivine-rich cumulate piles. Large mixing events probably destabilized the magma chamber, resulting in devolatilization and eruption. In contrast, small mixing events lacked the energy to destabilize the chamber, resulting in repeated compositional oscillations in minerals affected by multiple small mixing events.

Central aims of IODP Expedition 352 were to delineate and characterize the magmatic stratigraphy in the Bonin forearc to define key magmatic processes associated with subduction initiation and their potential links to ophiolites. Expedition 352 penetrated 1.2 km of magmatic basement at four sites and recovered three principal lithologies: tholeiitic forearc basalt (FAB), high-Mg andesite, and boninite, with subordinate andesite. Boninites are subdivided into basaltic, low-Si, and high-Si varieties. The purpose of this study is to determine conditions of crystal growth and differentiation for Expedition 352 lavas and compare and contrast these conditions with those recorded in lavas from mid-ocean ridges, forearcs, and ophiolites. Cr# (cationic Cr/Cr+Al) vs. TiO2 relations in spinel and clinopyroxene demonstrate a trend of source depletion with time for the Expedition 352 forearc basalt to boninite sequence that is similar to sequences in the Oman and other suprasubduction zone ophiolites. Clinopyroxene thermobarometry results indicate that FAB crystallized at temperatures (1142–1190 °C) within the range of MORB (1133–1240 °C). When taking into consideration liquid lines of descent of boninite, orthopyroxene barometry and olivine thermometry of Expedition 352 boninites demonstrate that they crystallized at temperatures marginally lower than those of FAB, between ~1119 and ~1202 °C and at relatively lower pressure (~0.2–0.4 vs. 0.5–4.6 kbar for FAB). Elevated temperatures of boninite orthopyroxene (~1214 °C for low-Si boninite and 1231–1264 °C for high-Si boninite) may suggest latent heat produced by the rapid crystallization of orthopyroxene. The lower pressure of crystallization of the boninite may be explained by their lower density and hence higher ascent rate, and shorter distance of travel from place of magma formation to site of crystallization, which allowed the more buoyant and faster ascending boninites to rise to shallower levels before crystallizing, thus preserving their high temperatures.

The Izu–Bonin–Mariana (IBM) forearc contains a record of a widespread and short-lived subduction initiation event in the early Eocene. We present new high-precision ages to determine the rates and length scales of ocean crust production following subduction initiation along the IBM convergent margin. We also explore the implications of these ages for potential causes of subduction initiation and its far-field effects. Zircons extracted from gabbros collected by Shinkai 6500 diving within the ophiolitic sequence exposed in the outer Bonin and Mariana forearcs give crystallization ages (±2σ) of 51.94 ± 0.13 Ma and 51.81 ± 0.03 Ma, respectively, consistent with a near-synchronous seafloor spreading after initiation of subduction along the entire length of the IBM system. 40Ar/39Ar ages of forearc basalts (51.34 ± 0.78 Ma) and boninites (51.27 ± 0.09 Ma to 50.33 ± 0.55 Ma) cored upslope of the Bonin sample gabbro during IODP Expedition 352 further reveal that the outer forearc comprises oceanic lithosphere with MORB to boninite compositions created by seafloor spreading during a remarkably short time period of 0.6–2 Myr at a full-rate of c. 70 mm/yr. This time period is consistent with those of some well-dated ophiolites. Thus, near-trench seafloor spreading to form basalt then boninite lithosphere is an extensive, rapid, but short-lived process after subduction initiation.

International audience

This article is a companion to Shervais et al. (2020), https://doi.org/10.1029/2020GC009093 .

The subduction of seamounts and ridge features at convergent plate boundaries plays an important role in the deformation of the overriding plate and influences geochemical cycling and associated biological processes. Active serpentinization of forearc mantle and serpentinite mud volcanism on the Mariana forearc (between the trench and active volcanic arc) provides windows on subduction processes. Here, we present (1) the first observation of an extensive exposure of an undeformed Cretaceous seamount currently being subducted at the Mariana Trench inner slope; (2) vertical deformation of the forearc region related to subduction of Pacific Plate seamounts and thickened crust; (3) recovered Ocean Drilling Program and International Ocean Discovery Program cores of serpentinite mudflows that confirm exhumation of various Pacific Plate lithologies, including subducted reef limestone; (4) petrologic, geochemical and paleontological data from the cores that show that Pacific Plate seamount exhumation covers greater spatial and temporal extents; (5) the inference that microbial communities associated with serpentinite mud volcanism may also be exhumed from the subducted plate seafloor and/or seamounts; and (6) the implications for effects of these processes with regard to evolution of life. This article is part of a discussion meeting issue ‘Serpentine in the Earth system’.

Tectonic discrimination diagrams are a key tool for understanding ancient volcanic rock origins. In this contribution we compile over 15,000 whole rock compositions to re-evaluate the Ti-V discrimination diagram and compare it to another commonly used tool, the Nb/Yb-Th/Yb diagram. We have reformulated the Ti-V diagram into a log–log plot to distinguish samples more clearly at lower concentrations. The compilation shows that MORB are dominated by Ti/V = 20–43, whereas juvenile arc tholeiites and boninites are characterized by Ti/V 43, although there may be significant overlap with MORB-like ratios for ridge-centered OIB and for some oceanic plateaux. About 56% of alkaline OIB have Ti/V > 70. Back-arc basins are dominantly MORB-like. Melt models show that MORB and juvenile arc volcanics most likely formed under different ƒO2 conditions, but are permissive of similar ƒO2 if the arc rocks form by much higher melt fractions. The Nb/Yb vs. Th/Yb plot clearly distinguishes most oceanic basalts (MORB, plateaux, OIB) from subduction-related volcanic rocks (boninite, juvenile arc tholeiite, calc-alkaline) and from flood basalts. We propose here a new two-proxy diagram of Ti/V vs. Th/Nb, which incorporates the advantages of both.

Orthopyroxene is especially suited to decode and testify to the behavior of highly immobile elements during hydrous mantle melting. Laser ablation ICP-MS analyses from orthopyroxene hosted within peridotite from the Coast Range ophiolite (CRO) demonstrates that Group A peridotites (lherzolites) have similar compositions to mid-ocean-ridge abyssal peridotite, whereas other peridotites (Groups B and C; harzburgites) retain depleted signatures, but display ‘spoon-shaped’ enrichments for the light-REE. These patterns are consistent with variable degrees of partial melting of MORB-source asthenosphere initiated within the garnet stability field (

Project HOTSPOT, the Snake River Scientific Drilling Project (International Continental Scientific Drilling Program), tested for deep geothermal resources and examined the petrology of volcanic rocks with three drillholes in the central and western Snake River Plain (western USA). The MH-2 drillhole targeted fractured crystalline and hydrothermally altered basalt in the area of the Mountain Home Air Force Base (Idaho) to a total depth of 1821 m. At 1745 m depth the drillhole encountered flowing artesian hydrothermal fluids of at least 150 °C. We integrate geological analyses of core, image log, and borehole geophysical data, and in situ stress analyses to describe the structural environment that produces permeability for artesian flow. The rocks in the lower 540 m of the drillhole consist of basalt flows as much as 30 m thick, altered basalt, and thin sedimentary horizons. The mechanical stratigraphy is defined by nine mechanical horizons that are in three ranges of rock strength on the basis of experimentally determined strength data, core logging, and geophysical log signatures. Hydrothermal alteration products and mineralization in the core are associated with three highly faulted sections; the lowermost section is associated with the zone of flowing thermal water. Shear slip indicators on faults observed in core indicate slip ranging from pure strike slip to normal failure mechanisms in the stronger horizons. The borehole breakouts indicate that the maximum horizontal stress, SH, is oriented 047° ± 7°, and drilling-induced tensile fractures indicate that SH is oriented at 67° ± 21°. The in situ stress orientations exhibit little variation over the depth of the measured interval, but the SH magnitude varies with depth, and is best explained by an oblique normal fault stress regime. The geomechanical model indicates that if pore pressures at depth are elevated above the normal hydrostatic gradient, as observed here, the system has the potential to deform by mixed normal and strike-slip failure. Our observations and interpretations suggest that the MH-2 borehole was drilled into oblique normal faults that intersect a buried 300°-trending fault block masked by the basaltic volcanic complex. These data indicate that the transition from the central to western Snake River Plain is characterized by complex structures developed in response to a transitional stress state related to Snake River Plain and western Basin and Range stress regimes. The western Basin and Range stress and tectonic regime may extend from northern Nevada into western Idaho and may enhance the potential for geothermal resources by creating interconnected fracture and fault-related permeability at depth.

International audience; nternational Ocean Discovery Program (IODP) Expedition 352 recovered a high-fidelity record of volcanism related to subduction initiation in the Bonin fore-arc. Two sites (U1440 and U1441) located in deep water nearer to the trench recovered basalts and related rocks; two sites (U1439 and U1442) located in shallower water further from the trench recovered boninites and related rocks. Drilling in both areas ended in dolerites inferred to be sheeted intrusive rocks. The basalts apparently erupted immediately after subduction initiation and have compositions similar to those of the most depleted basalts generated by rapid sea-floor spreading at mid-ocean ridges, with little or no slab input. Subsequent melting to generate boninites involved more depleted mantle and hotter and deeper subducted components as subduction progressed and volcanism migrated away from the trench. This volcanic sequence is akin to that recorded by many ophiolites, supporting a direct link between subduction initiation, fore-arc spreading, and ophiolite genesis.

Volcanic glasses are rarely preserved in the rock record, and the quality of preservation generally declines with increasing age. Records preserved in ancient basaltic glasses therefore provide important links between processes operating in the distant past, and those that are active on the Earth today. Microbial colonization has been linked to the formation of characteristic structures in basaltic glass, including tubules and granule-filled tubules, which are thought to be produced by microbially mediated glass dissolution. Structures of similar occurrence and morphology but filled almost entirely with fine-grained titanite have been documented in some ancient metabasalts. It has been suggested that the ancient titanite-mineralized structures are mineralized equivalents of hollow tubules in modern glassy basaltic rocks, but a direct link has not been firmly established. We report the discovery of tubular bioalteration structures in fresh and minimally altered basaltic glasses of middle Jurassic (164 Ma) age from the Stonyford Volcanic Complex (SFVC), Coast Range Ophiolite, California. Tubular structures hosted in unaltered basaltic glass are typically hollow, whilst those in zones of zeolitic alteration are mineralized by titanite. Tubules are continuous across zeolite-glass interfaces, which mark an abrupt change from titanite-filled to hollow tubules, demonstrating that titanite growth occurs preferentially within pre-existing tubular structures. Titanite mineralization in the SFVC represent a link between tubular structures in modern basaltic glass and titanite-mineralized features of similar morphology and spatial distribution in ancient metabasalts. Our observations support a link between textures in modern glassy basaltic rocks and some of the oldest-known putative ichnofossils.

International audience; Subduction initiation is one of the least understood aspects of plate tectonics. In an effort to obtain the first in situ magmatic record of subduction initiation, the International Ocean Discovery Program Expedition 352 drilled at four sites in the inner trench wall of the Bonin Trench to recover 1.22 km of oceanic upper crust accreted within a few m.y. of subduction initiation. The two sites nearer to the trench (U1440 and U1441) yielded axial and off-axis fore-arc basalts (FAB), while those c. 15 km further from the trench (U1439 and U1442) yielded axial low-silica boninites and high-Mg andesites overlain by off-axis high-silica boninites. This study uses Hf–Nd–Sr–Pb isotope analyses from c. 50 stratigraphically representative core samples to trace the evolution of the mantle source during the brief period of FAB-through-boninite magmatism immediately following subduction initiation. Results show that: 1) the FAB have high εHf relative to εNd and were derived from variably depleted mantle of ‘Indian’ provenance with no detectable subduction input; 2) the axial boninites follow mixing trends between a residual FAB mantle source and a subduction component derived from shallow (amphibolite facies) melting of oceanic crust of ‘Pacific’ provenance; and 3) the off-axis boninites define mixing trends between a hybrid mantle wedge (residual mantle + slab melt) and an additional subduction component with lower εNd and higher 207Pb/204Pb that requires a significant contribution from pelagic sediment. This incoming of pelagic sediments may signify a change from an accretionary to non-accretionary margin as subduction evolves. The results thus indicate a rapidly evolving system in terms of geodynamics, magma genesis and crustal accretion immediately following subduction initiation.

The Snake River Plain of central Idaho represents the world’s best example of a mantle hotspot track impinging upon continental crust and provides a record of bimodal volcanism extending over 12 Ma to the present. Project Hotspot recovered almost 2 km of continuous drill core from the Kimama borehole, located in central Idaho on the axial volcanic zone. The Kimama drill core represents the most complete record of mafic volcanism along the Yellowstone–Snake River Plain hotspot track. A total of 432 basalt flow units, representing 183 basalt flows, 78 basalt flow groups, and 34 super groups, along with 42 sediment interbeds are recognized using volcanic facies observations, stratigraphic relationships, borehole geophysical logs, and paleosecular variation in magnetostratigraphy. Rhyolite and other non-basaltic volcanic materials were not encountered in the drill core. Ages for six basalt lava flows were determined by 40Ar/39Ar using incremental heating experiments. Paleomagnetic inclination was measured on over 1200 samples collected at roughly 2-m-depth intervals, yielding mean values of paleosecular variation between ±50° to ±70° in Kimama flow groups, close to the expected 61° axial dipole average for the Kimama borehole location. Twenty-three magnetic reversals were identified and correlated to dated geomagnetic chrons and subchrons and compared with the 40Ar/39Ar radiometric ages. A linear fit to 40Ar/39Ar dates, geomagnetic chron and subchron boundaries, and volcanogenic zircon U-Pb ages defines a mean accumulation rate of ∼320 m/m.y. and extrapolates to a bottom hole age of 6.3 Ma. Average thicknesses of lithologic units increase from 2.7 m (sediment), 4 m (flow units), 10 m (flows), 23 m (flow groups), to 53 m (super groups). On average, one lava flow inundated the Kimama borehole location every 33 k.y. Intercalated sediments, ranging from 0.06 to 24.5 m thick, make up roughly 6% of the drill core and indicate lulls in local volcanic activity that may have lasted up to 77 k.y. Neutron and gamma-ray logs supplement observations from the drill cores: neutron logs document individual flow units through the contrast between massive flow interiors and more porous flow surfaces, and gamma-ray logs document the depth and thickness of sedimentary interbeds and high–K-Fe basalts. The 5.8 m.y. duration of basaltic volcanism in the Kimama drill core implies a steady rate of volcanism, indicating a relatively stable rate of mantle upflow along the lithosphere-mantle boundary in the wake of Yellowstone–Snake River Plain plume volcanism.

The Izu‐Bonin‐Mariana (IBM) fore arc preserves igneous rock assemblages that formed during subduction initiation circa 52 Ma. International Ocean Discovery Program (IODP) Expedition 352 cored four sites in the fore arc near the Ogasawara Plateau in order to document the magmatic response to subduction initiation and the physical, petrologic, and chemical stratigraphy of a nascent subduction zone. Two of these sites (U1440 and U1441) are underlain by fore‐arc basalt (FAB). FABs have mid‐ocean ridge basalt (MORB)‐like compositions, however, FAB are consistently lower in the high‐field strength elements (TiO2, P2O5, Zr) and Ni compared to MORB, with Na2O at the low end of the MORB field and FeO* at the high end. Almost all FABs are light rare earth element depleted, with low total REE, and have low ratios of highly incompatible to less incompatible elements (Ti/V, Zr/Y, Ce/Yb, and Zr/Sm) relative to MORB. Chemostratigraphic trends in Hole U1440B are consistent with the uppermost lavas forming off axis, whereas the lower lavas formed beneath a spreading center axis. Axial magma of U1440B becomes more fractionated upsection; overlying off‐axis magmas return to more primitive compositions. Melt models require a two‐stage process, with early garnet field melts extracted prior to later spinel field melts, with up to 23% melting to form the most depleted compositions. Mantle equilibration temperatures are higher than normal MORB (1,400 °C–1,480 °C) at relatively low pressures (1–2 GPa), which may reflect an influence of the Manus plume during subduction initiation. Our data support previous models of FAB origin by decompression melting but imply a source more depleted than normal MORB source mantle., Key Points Fore‐arc basalts (FABs) formed by decompression melting in response to subduction initiation and differ from mid‐ocean ridge basaltsFABs form by two‐stage melting, with early garnet field melts extracted prior to spinel field melting, resulting in LREE/HREE depletionHighly depleted FAB may reflect both an older depletion event and higher ambient temperatures related to the Manus plume

Basalts erupted in the Snake River Plain of central Idaho and sampled in the Kimama drill core link eruptive processes to the construction of mafic intrusions over 5.5 Ma. Cyclic variations in basalt composition reveal temporal chemical heterogeneity related to fractional crystallization and the assimilation of previously-intruded mafic sills. A range of compositional types are identified within 1912 m of continuous drill core: Snake River olivine tholeiite (SROT), low K SROT, high Fe-Ti, and evolved and high K-Fe lavas similar to those erupted at Craters of the Moon National Monument. Detailed lithologic and geophysical logs document 432 flow units comprising 183 distinct lava flows and 78 flow groups. Each lava flow represents a single eruptive episode, while flow groups document chemically and temporally related flows that formed over extended periods of time. Temporal chemical variation demonstrates the importance of source heterogeneity and magma processing in basalt petrogenesis. Low-K SROT and high Fe-Ti basalts are genetically related to SROT as, respectively, hydrothermally-altered and fractionated daughters. Cyclic variations in the chemical composition of Kimama flow groups are apparent as 21 upward fractionation cycles, six recharge cycles, eight recharge-fractionation cycles, and five fractionation-recharge cycles. We propose that most Kimama basalt flows represent typical fractionation and recharge patterns, consistent with the repeated influx of primitive SROT parental magmas and extensive fractional crystallization coupled with varying degrees of assimilation of gabbroic to ferrodioritic sills at shallow to intermediate depths over short durations. Trace element models show that parental SROT basalts were generated by 5-10% partial melting of enriched mantle at shallow depths above the garnet-spinel lherzolite transition. The distinctive evolved and high K-Fe lavas are rare. Found at four depths, 319 m, 1045 m, 1078 m, and 1189 m, evolved and high K-Fe flows are compositionally unrelated to SROT magmas and represent highly fractionated basalt, probably accompanied by crustal assimilation. These evolved lavas may be sourced from the Craters of the Moon/Great Rift system to the northeast. The Kimama drill core is the longest record of geochemical variation in the central Snake River Plain and reinforces the concept of magma processing in a layered complex.

Hotspot: The Snake River Geothermal Drilling Project was undertaken to better understand geothermal systems across the Snake River Plain volcanic province. A series of surface and borehole seismic profiles were obtained to provide insights into volcanic stratigraphy and test the capabilities of engineering-scale seismic imaging in such terranes. The Kimberly site drilled through 1.9 km of mostly rhyolite, with thin sedimentary interbeds in the upper part of the section. The Kimama site drilled through 1.9 km of mostly basalt with sedimentary interbeds at 200 m depth and 1700 m depth. The Mountain Home site contained numerous sediment and volcanic rock layers. Downhole and surface vibroseis seismic results suggest sedimentary interbeds at depth correspond with low-velocity, high-temperature zones that relate to reflections on seismic profiles. Our results suggest that eruption flow volumes can be estimated and flow boundaries can be imaged with surface seismic methods using relatively high-fold and wide-angle coverage. High-frequency attenuation is observed at all sites, and this deficit may be countered by acquisition design and a focus on signal processing steps. Separation of surface and body waves was obtained by muting, and the potential for large static effects was identified and addressed in processing. An accurate velocity model and lithology contacts derived from borehole information improved the confidence of our seismic interpretations.

Convergent plate margins are key features of our dynamic planet, and the processes taking place along such margins have had a correspondingly large impact on society and the geologic record. As wit...

A 1912-m exploration corehole was drilled along the axis of the eastern Snake River Plain, Idaho. Two temperature logs run on the corehole display an obvious inflection point at about 960 m. Such behavior is indicative of downward fluid flow in the wellbore. The geothermal gradient above 935 m is 4.5 °C/km, while the gradient is 72–75 °C/km from 980 to 1440 m. Projecting the higher gradients upward to where they intersect the lower gradient on the temperature logs places the bottom of the cold, freshwater Snake River Plain aquifer, which suppresses the geothermal gradient at this location, at least 860 m below the surface. The average heat flow for the corehole between 983 and 1550 m is 132 mW/m2. Although the maximum bottom-hole temperature extrapolated from a measured time–temperature curve was only 59.3 °C, geothermometers suggest an equilibrium temperature on the order of 125–140 °C based on a single fluid sample from 1070 m. Furthermore, below 960 m the basalt core shows obvious signs of alteration, including a distinct color change, the formation of smectite clay, and the presence of secondary minerals filling vesicles and fracture zones. This alteration boundary could act as an effective cap or seal for a hot-water geothermal system.

Coordinated drilling efforts are an important method to investigate active tectonics and magmatic processes related to faults and volcanoes. The US National Science Foundation (NSF) recently sponsored a series of workshops to define the nature of future continental drilling efforts. As part of this series, we convened a workshop to explore how continental scientific drilling can be used to better understand active tectonic and magmatic processes. The workshop, held in Park City, Utah, in May 2013, was attended by 41 investigators from seven countries. Participants were asked to define compelling scientific justifications for examining problems that can be addressed by coordinated programs of continental scientific drilling and related site investigations. They were also asked to evaluate a wide range of proposed drilling projects, based on white papers submitted prior to the workshop. Participants working on faults and fault zone processes highlighted two overarching topics with exciting potential for future scientific drilling research: (1) the seismic cycle and (2) the mechanics and architecture of fault zones. Recommended projects target fundamental mechanical processes and controls on faulting, and range from induced earthquakes and earthquake initiation to investigations of detachment fault mechanics and fluid flow in fault zones. Participants working on active volcanism identified five themes: the volcano eruption cycle; eruption sustainability, near-field stresses, and system recovery; eruption hazards; verification of geophysical models; and interactions with other Earth systems. Recommended projects address problems that are transferrable to other volcanic systems, such as improved methods for identifying eruption history and constraining the rheological structure of shallow caldera regions. Participants working on chemical geodynamics identified four major themes: large igneous provinces (LIPs), ocean islands, continental hotspot tracks and rifts, and convergent plate margins (subduction zones). This workshop brought together a diverse group of scientists with a broad range of scientific experience and interests. A particular strength was the involvement of both early-career scientists, who will initiate and carry out these new research programs, and more senior researchers with many years of experience in scientific drilling and active tectonics research. Each of the themes and questions outlined above has direct benefits to society, including improving hazard assessment, direct monitoring of active systems for early warning, renewable and non-renewable resource and energy exploitation, and predicting the environmental impacts of natural hazards, emphasizing the central role that scientific drilling will play in future scientific and societal developments.

The physical and chemical evolution of the Earth is driven by geodynamic cycles that are global in scale, operating over 4.57 Ga of Earth’s history. Some processes are truly cyclic, e.g., the Wilson Cycle, while others are irreversible (e.g., core formation). Heat and mass transfer between the lowermost mantle (e.g., core-mantle boundary) and the surface drives these global geodynamic processes. Subduction of lithospheric plates transfers cool fractionated material into the lower mantle and leads indirectly to the formation of new oceanic lithosphere, while the rise of thermochemical plumes recycles the remnants of these plates back to the surface, driven by heat transfer across the core–mantle boundary. These global geodynamic cycles are responsible for hotspot volcanism, the formation of continental crust, collisional orogenies, continental rifting, subduction zone processes (arcs, accretionary prisms), and ore deposits. Each of these presents opportunities for investigation by continental scientific drilling. In addition, these cycles affect other processes that are targets of continental scientific drilling: the origin and evolution of life and an oxygenated atmosphere, the impact of large volcanic eruptions on climate, and geological hazards such as earthquakes and volcanic eruptions. In this paper, we present the scientific rationale for continental scientific drilling to study global geodynamic processes, review past successes in this realm that were sponsored in part by ICDP, and suggest potential new targets for drilling campaigns that focus on solid earth evolution. This paper builds on discussions at the 2013 ICDP Science Meeting on the future of continental scientific drilling, held in Potsdam in November 2013.

We examine the fine-scale variations in mineralogical composition, geochemical alteration, and texture of the fault-related rocks from the Phase 3 whole-rock core sampled between 3,187.4 and 3,301.4 m measured depth within the San Andreas Fault Observatory at Depth (SAFOD) borehole near Parkfield, California. This work provides insight into the physical and chemical properties, structural architecture, and fluid-rock interactions associated with the actively deforming traces of the San Andreas Fault zone at depth. Exhumed outcrops within the SAF system comprised of serpentinite-bearing protolith are examined for comparison at San Simeon, Goat Rock State Park, and Nelson Creek, California. In the Phase 3 SAFOD drillcore samples, the fault-related rocks consist of multiple juxtaposed lenses of sheared, foliated siltstone and shale with block-in-matrix fabric, black cataclasite to ultracataclasite, and sheared serpentinite-bearing, finely foliated fault gouge. Meters-wide zones of sheared rock and fault gouge correlate to the sites of active borehole casing deformation and are characterized by scaly clay fabric with multiple discrete slip surfaces or anastomosing shear zones that surround conglobulated or rounded clasts of compacted clay and/or serpentinite. The fine gouge matrix is composed of Mg-rich clays and serpentine minerals (saponite ± palygorskite, and lizardite ± chrysotile). Whole-rock geochemistry data show increases in Fe-, Mg-, Ni-, and Cr-oxides and hydroxides, Fe-sulfides, and C-rich material, with a total organic content of >1 % locally in the fault-related rocks. The faults sampled in the field are composed of meters-thick zones of cohesive to non-cohesive, serpentinite-bearing foliated clay gouge and black fine-grained fault rock derived from sheared Franciscan Formation or serpentinized Coast Range Ophiolite. X-ray diffraction of outcrop samples shows that the foliated clay gouge is composed primarily of saponite and serpentinite, with localized increases in Ni- and Cr-oxides and C-rich material over several meters. Mesoscopic and microscopic textures and deformation mechanisms interpreted from the outcrop sites are remarkably similar to those observed in the SAFOD core. Micro-scale to meso-scale fabrics observed in the SAFOD core exhibit textural characteristics that are common in deformed serpentinites and are often attributed to aseismic deformation with episodic seismic slip. The mineralogy and whole-rock geochemistry results indicate that the fault zone experienced transient fluid–rock interactions with fluids of varying chemical composition, including evidence for highly reducing, hydrocarbon-bearing fluids.

The Snake River Plain represents 17 m.y. of volcanic activity that took place as the North American continent migrated over a relatively fixed magma source, or hotspot. We present new Pb, Sr, and Nd data for a suite of 25 basalts collected from Western and Central Snake River Plain (SRP). The new isotope data, combined with previously published data from the SRP, provide a traverse of the Wyoming craton margin, from the 87Sr/86Sr = 0.706 line boundary of western SRP with Phanerozoic accreted terranes, east through the central and eastern SRP, to the Yellowstone Plateau. Low-K basalts from the western SRP, overlain by high-K basalts, provide a temporal record of regional source variation from ∼16.8 to 0.2 Ma. Principal Component Analysis (PCA) of the new and previously published SRP basalt Pb isotopes reveals that >97% of the total variability is accounted for by mixing between three end-members and is consistent with a sublithospheric Yellowstone hotspot mantle source with a radiogenic isotope composition similar to the mantle source of the early Columbia River Basalt Group (CRBG) and two continental lithosphere end-members, heterogeneous in age and composition. We use the SRP Pb, Sr, and Nd isotope data to model the Yellowstone Hotspot–continental lithosphere interaction by three component mixing between two continental lithospheric components, Archean lithosphere (CL1) that represents older lithosphere underlying the Yellowstone Plateau in the east, and Paleoproterozoic lithosphere (CL2) representing the younger lithosphere underlying the SRP in the west near the craton margin, and a sublithospheric end-member, representing the Yellowstone hotspot (PL). The results suggest a continuous flow of PL material westward as the NA continental lithosphere migrated over the upwelling hotspot along a shoaling gradient in the sub-continental mantle lithosphere. The model shows a decrease in Total Lithosphere end-members (CL1 + CL2) and the Lithosphere Ratio (CL1/CL2), from the craton interior at Yellowstone toward its western margin, consistent with geologic and geophysical evidence that the continental lithosphere beneath the SRP decreases in age and thickness from east to west. The Lithosphere Ratio shows step-like decreases from Yellowstone in the east to the 87Sr/86Sr = 0.706 line in the west, indicating that the SRP cuts across geochemically distinct parcels of lithospheric mantle, consistent with terrane accretion models for the craton margin. In the western SRP, young high-K basalts have a lower mass fraction of Total Lithospheric compared to the underlying low-K tholeiites, but the same Lithosphere Ratio, consistent with a recent (700–900 ka) decrease in lithosphere contribution between eruption of early low- and younger high-K basalts.

[1] There are three methods of investigating the oxygen-isotope composition of oceanic crust: (1) data from drilling-obtained samples; (2) study of obducted ophiolites; and (3) from crustal-derived xenoliths; however, each approach has limitations. Drilling capabilities are limited to the upper 1–3 km of the oceanic crust. Ophiolites are segments of oceanic crust and underlying mantle exposed at the Earth's surface, but their origin and emplacement are still debated. Nonetheless, the study of ophiolite sequences has enabled many of the fundamental observations regarding ancient-oceanic crust and the oxygen-isotope signature of the upper mantle. Analyses of ophiolite sequences (e.g., Samail) reveal δ18O variations ranging from isotopically enriched in relation to MORB (δ18O > MORB) to isotopically depleted (low δ18O ≤ MORB), interpreted to show the effects of seawater interaction at different temperatures. Mantle-derived xenoliths also display compelling geochemical evidence suggesting they are remnants of subducted ocean crust, and may in fact, be more representative samples. Some kimberlite-hosted diamondiferous eclogite xenoliths, with crustal protoliths, reveal a marked oxygen-isotope disparity between oceanic crust and the upper mantle; i.e., a bimodal distribution of both mantle-like (δ18O ≈ 5.6‰) and isotopically enriched oxygen (δ18O ≈ 6.8‰). The isotopically high δ18O values are often observed in eclogite xenoliths with ultramafic protoliths, which form beyond the low-temperature alteration zone. This finding is in contrast to traditional oxygen-isotope versus depth-curves of ophiolite sequences. Here, we present a compiled data set from Siberian eclogites and pyroxenites that illustrates this bimodality, using magnesium number (Mg#) as a depth proxy, and investigate the conditions that could potentially result in the observed “abnormal” oxygen-isotope variations.