Magmatism controls global oceanic transform fault topography.

Oceanic transform faults play an essential role in plate tectonics. Yet to date, there is no unifying explanation for the global trend in broad-scale transform fault topography, ranging from deep valleys to shallow topographic highs. Using three-dimensional numerical models, we find that spreading-rate dependent magmatism within the transform domain exerts a first-order control on the observed spectrum of transform fault depths. Low-rate magmatism results in deep transform valleys caused by transform-parallel tectonic stretching; intermediate-rate magmatism fully accommodates far-field stretching, but strike-slip motion induces across-transform tension, producing transform strength dependent shallow valleys; high-rate magmatism produces elevated transform zones due to local compression. Our models also address the observation that fracture zones are consistently shallower than their adjacent transform fault zones. These results suggest that plate motion change is not a necessary condition for reproducing oceanic transform topography and that oceanic transform faults are not simple conservative strike-slip plate boundaries.

Deglaciation-enhanced mantle CO2 fluxes at Yellowstone imply positive climate feedback.

Mantle melt generation in response to glacial unloading has been linked to enhanced magmatic volatile release in Iceland and global eruptive records. It is unclear whether this process is important in systems lacking evidence of enhanced eruptions. The deglaciation of the Yellowstone ice cap did not observably enhance volcanism, yet Yellowstone emits large volumes of CO2 due to melt crystallization at depth. Here we model mantle melting and CO2 release during the deglaciation of Yellowstone (using Iceland as a benchmark). We find mantle melting is enhanced 19-fold during deglaciation, generating an additional 250-620 km3. These melts segregate an additional 18-79 Gt of CO2 from the mantle, representing a ~3-15% increase in the global volcanic CO2 flux (if degassed immediately). We suggest deglaciation-enhanced mantle melting is important in continental settings with partially molten mantle - including Greenland and West Antarctica - potentially implying positive feedbacks between deglaciation and climate warming.

How transform fault shear influences where detachment faults form near mid-ocean ridges.

Oceanic detachment faults represent an end-member form of seafloor creation, associated with relatively weak magmatism at slow-spreading mid-ocean ridges. We use 3-D numerical models to investigate the underlying mechanisms for why detachment faults predominantly form on the transform side (inside corner) of a ridge-transform intersection as opposed to the fracture zone side (outside corner). One hypothesis for this behavior is that the slipping, and hence weaker, transform fault allows for the detachment fault to form on the inside corner, and a stronger fracture zone prevents the detachment fault from forming on the outside corner. However, the results of our numerical models, which simulate different frictional strengths in the transform and fracture zone, do not support the first hypothesis. Instead, the model results, combined with evidence from rock physics experiments, suggest that shear-stress on transform fault generates excess lithospheric tension that promotes detachment faulting on the inside corner.

Sensitivity of rift tectonics to global variability in the efficiency of river erosion.

SignificanceThe efficiency of erosion in leveling relief mainly depends on climate and strength of exposed rocks. However, whether erosion is sufficiently efficient to influence the architecture of a tectonic plate boundary remains a topic of debate. Here, we analyze continental rift landscapes reworked by river incision to assess a globally representative range of fluvial erosion efficiency. We then simulate crustal extension exposed to surface processes acting within this documented range. We find that more efficient erosion favors the growth of half-grabens over horsts, which can explain contrasting tectonic styles across the Basin and Range province and the East African Rift. This suggests that variability in Earth's geological structures partly reflects variability in hydrological conditions and associated surface processes.

High 3He/4He in central Panama reveals a distal connection to the Galápagos plume.

It is well established that mantle plumes are the main conduits for upwelling geochemically enriched material from Earth's deep interior. The fashion and extent to which lateral flow processes at shallow depths may disperse enriched mantle material far (>1,000 km) from vertical plume conduits, however, remain poorly constrained. Here, we report He and C isotope data from 65 hydrothermal fluids from the southern Central America Margin (CAM) which reveal strikingly high 3He/4He (up to 8.9RA) in low-temperature (≤50 °C) geothermal springs of central Panama that are not associated with active volcanism. Following radiogenic correction, these data imply a mantle source 3He/4He >10.3RA (and potentially up to 26RA, similar to Galápagos hotspot lavas) markedly greater than the upper mantle range (8 ± 1RA). Lava geochemistry (Pb isotopes, Nb/U, and Ce/Pb) and geophysical constraints show that high 3He/4He values in central Panama are likely derived from the infiltration of a Galápagos plume-like mantle through a slab window that opened ∼8 Mya. Two potential transport mechanisms can explain the connection between the Galápagos plume and the slab window: 1) sublithospheric transport of Galápagos plume material channeled by lithosphere thinning along the Panama Fracture Zone or 2) active upwelling of Galápagos plume material blown by a "mantle wind" toward the CAM. We present a model of global mantle flow that supports the second mechanism, whereby most of the eastward transport of Galápagos plume material occurs in the shallow asthenosphere. These findings underscore the potential for lateral mantle flow to transport mantle geochemical heterogeneities thousands of kilometers away from plume conduits.

Hydraulic transmissivity inferred from ice-sheet relaxation following Greenland supraglacial lake drainages.

Surface meltwater reaching the base of the Greenland Ice Sheet transits through drainage networks, modulating the flow of the ice sheet. Dye and gas-tracing studies conducted in the western margin sector of the ice sheet have directly observed drainage efficiency to evolve seasonally along the drainage pathway. However, the local evolution of drainage systems further inland, where ice thicknesses exceed 1000 m, remains largely unknown. Here, we infer drainage system transmissivity based on surface uplift relaxation following rapid lake drainage events. Combining field observations of five lake drainage events with a mathematical model and laboratory experiments, we show that the surface uplift decreases exponentially with time, as the water in the blister formed beneath the drained lake permeates through the subglacial drainage system. This deflation obeys a universal relaxation law with a timescale that reveals hydraulic transmissivity and indicates a two-order-of-magnitude increase in subglacial transmissivity (from 0.8 ± 0.3 [Formula: see text] to 215 ± 90.2 [Formula: see text]) as the melt season progresses, suggesting significant changes in basal hydrology beneath the lakes driven by seasonal meltwater input.

Aseismic transient slip on the Gofar transform fault, East Pacific Rise.

Oceanic transform faults display a unique combination of seismic and aseismic slip behavior, including a large globally averaged seismic deficit, and the local occurrence of repeating magnitude (M) [Formula: see text] earthquakes with abundant foreshocks and seismic swarms, as on the Gofar transform of the East Pacific Rise and the Blanco Ridge in the northeast Pacific Ocean. However, the underlying mechanisms that govern the partitioning between seismic and aseismic slip and their interaction remain unclear. Here we present a numerical modeling study of earthquake sequences and aseismic transient slip on oceanic transform faults. In the model, strong dilatancy strengthening, supported by seismic imaging that indicates enhanced fluid-filled porosity and possible hydrothermal circulation down to the brittle-ductile transition, effectively stabilizes along-strike seismic rupture propagation and results in rupture barriers where aseismic transients arise episodically. The modeled slow slip migrates along the barrier zones at speeds ∼10 to 600 m/h, spatiotemporally correlated with the observed migration of seismic swarms on the Gofar transform. Our model thus suggests the possible prevalence of episodic aseismic transients in M [Formula: see text] rupture barrier zones that host active swarms on oceanic transform faults and provides candidates for future seafloor geodesy experiments to verify the relation between aseismic fault slip, earthquake swarms, and fault zone hydromechanical properties.

Grain-size dynamics beneath mid-ocean ridges: Implications for permeability and melt extraction.

Grain size is an important control on mantle viscosity and permeability, but is difficult or impossible to measure in situ. We construct a two-dimensional, single phase model for the steady state mean grain size beneath a mid-ocean ridge. The mantle rheology is modeled as a composite of diffusion creep, dislocation creep, dislocation accommodated grain boundary sliding, and a plastic stress limiter. The mean grain size is calculated by the paleowattmeter relationship of Austin and Evans (2007). We investigate the sensitivity of our model to global variations in grain growth exponent, potential temperature, spreading-rate, and mantle hydration. We interpret the mean grain-size field in terms of its permeability to melt transport. The permeability structure due to mean grain size may be approximated as a high permeability region beneath a low permeability region. The transition between high and low permeability regions occurs across a boundary that is steeply inclined toward the ridge axis. We hypothesize that such a permeability structure generated from the variability of the mean grain size may focus melt toward the ridge axis, analogous to Sparks and Parmentier (1991)-type focusing. This focusing may, in turn, constrain the region where significant melt fractions are observed by seismic or magnetotelluric surveys. This interpretation of melt focusing via the grain-size permeability structure is consistent with MT observation of the asthenosphere beneath the East Pacific Rise. KEY POINTS: The grain-size field beneath MORs can vary over orders of magnitude The grain-size field affects the rheology and permeability of the asthenosphere The grain-size field may focus melt toward the ridge axis.

Greenland supraglacial lake drainages triggered by hydrologically induced basal slip.

Water-driven fracture propagation beneath supraglacial lakes rapidly transports large volumes of surface meltwater to the base of the Greenland Ice Sheet. These drainage events drive transient ice-sheet acceleration and establish conduits for additional surface-to-bed meltwater transport for the remainder of the melt season. Although it is well established that cracks must remain water-filled to propagate to the bed, the precise mechanisms that initiate hydro-fracture events beneath lakes are unknown. Here we show that, for a lake on the western Greenland Ice Sheet, drainage events are preceded by a 6-12 hour period of ice-sheet uplift and/or enhanced basal slip. Our observations from a dense Global Positioning System (GPS) network allow us to determine the distribution of meltwater at the ice-sheet bed before, during, and after three rapid drainages in 2011-2013, each of which generates tensile stresses that promote hydro-fracture beneath the lake. We hypothesize that these precursors are associated with the introduction of meltwater to the bed through neighbouring moulin systems (vertical conduits connecting the surface and base of the ice sheet). Our results imply that as lakes form in less crevassed, interior regions of the ice sheet, where water at the bed is currently less pervasive, the creation of new surface-to-bed conduits caused by lake-draining hydro-fractures may be limited.

Foundering of lower island-arc crust as an explanation for the origin of the continental Moho.

A long-standing theory for the genesis of continental crust is that it is formed in subduction zones. However, the observed seismic properties of lower crust and upper mantle in oceanic island arcs differ significantly from those in the continental crust. Accordingly, significant modifications of lower arc crust must occur, if continental crust is indeed formed from island arcs. Here we investigate how the seismic characteristics of arc crust are transformed into those of the continental crust by calculating the density and seismic structure of two exposed sections of island arc (Kohistan and Talkeetna). The Kohistan crustal section is negatively buoyant with respect to the underlying depleted upper mantle at depths exceeding 40 kilometres and is characterized by a steady increase in seismic velocity similar to that observed in active arcs. In contrast, the lower Talkeetna crust is density sorted, preserving only relicts (about ten to a hundred metres thick) of rock with density exceeding that of the underlying mantle. Specifically, the foundering of the lower Talkeetna crust resulted in the replacement of dense mafic and ultramafic cumulates by residual upper mantle, producing a sharp seismic discontinuity at depths of around 38 to 42 kilometres, characteristic of the continental Mohorovicic discontinuity (the Moho). Dynamic calculations indicate that foundering is an episodic process that occurs in most arcs with a periodicity of half a million to five million years. Moreover, because foundering will continue after arc magmatism ceases, this process ultimately results in the formation of the continental Moho.

Fracture propagation to the base of the Greenland Ice Sheet during supraglacial lake drainage.

Surface meltwater that reaches the base of an ice sheet creates a mechanism for the rapid response of ice flow to climate change. The process whereby such a pathway is created through thick, cold ice has not, however, been previously observed. We describe the rapid (<2 hours) drainage of a large supraglacial lake down 980 meters through to the bed of the Greenland Ice Sheet initiated by water-driven fracture propagation evolving into moulin flow. Drainage coincided with increased seismicity, transient acceleration, ice-sheet uplift, and horizontal displacement. Subsidence and deceleration occurred over the subsequent 24 hours. The short-lived dynamic response suggests that an efficient drainage system dispersed the meltwater subglacially. The integrated effect of multiple lake drainages could explain the observed net regional summer ice speedup.

Intermittent plate tectonics?

Although it is commonly assumed that subduction has operated continuously on Earth without interruption, subduction zones are routinely terminated by ocean closure and supercontinent assembly. Under certain circumstances, this could lead to a dramatic loss of subduction, globally. Closure of a Pacific-type basin, for example, would eliminate most subduction, unless this loss were compensated for by comparable subduction initiation elsewhere. Given the evidence for Pacific-type closure in Earth's past, the absence of a direct mechanism for termination/initiation compensation, and recent data supporting a minimum in subduction flux in the Mesoproterozoic, we hypothesize that dramatic reductions or temporary cessations of subduction have occurred in Earth's history. Such deviations in the continuity of plate tectonics have important consequences for Earth's thermal and continental evolution.

Trench-parallel anisotropy produced by foundering of arc lower crust.

Many volcanic arcs display fast seismic shear-wave velocities parallel to the strike of the trench. This pattern of anisotropy is inconsistent with simple models of corner flow in the mantle wedge. Although several models, including slab rollback, oblique subduction, and deformation of water-rich olivine, have been proposed to explain trench-parallel anisotropy, none of these mechanisms are consistent with all observations. Instead, small-scale convection driven by the foundering of dense arc lower crust provides an explanation for the trench-parallel anisotropy, even in settings with orthogonal convergence and no slab rollback.

Spreading rate dependence of gravity anomalies along oceanic transform faults.

Mid-ocean ridge morphology and crustal accretion are known to depend on the spreading rate of the ridge. Slow-spreading mid-ocean-ridge segments exhibit significant crustal thinning towards transform and non-transform offsets, which is thought to arise from a three-dimensional process of buoyant mantle upwelling and melt migration focused beneath the centres of ridge segments. In contrast, fast-spreading mid-ocean ridges are characterized by smaller, segment-scale variations in crustal thickness, which reflect more uniform mantle upwelling beneath the ridge axis. Here we present a systematic study of the residual mantle Bouguer gravity anomaly of 19 oceanic transform faults that reveals a strong correlation between gravity signature and spreading rate. Previous studies have shown that slow-slipping transform faults are marked by more positive gravity anomalies than their adjacent ridge segments, but our analysis reveals that intermediate and fast-slipping transform faults exhibit more negative gravity anomalies than their adjacent ridge segments. This finding indicates that there is a mass deficit at intermediate- and fast-slipping transform faults, which could reflect increased rock porosity, serpentinization of mantle peridotite, and/or crustal thickening. The most negative anomalies correspond to topographic highs flanking the transform faults, rather than to transform troughs (where deformation is probably focused and porosity and alteration are expected to be greatest), indicating that crustal thickening could be an important contributor to the negative gravity anomalies observed. This finding in turn suggests that three-dimensional magma accretion may occur near intermediate- and fast-slipping transform faults.

Automated analysis of electrospray ionization fourier transform ion cyclotron resonance mass spectra of natural organic matter.

The advent of ultra-high-resolution mass spectrometry has revolutionized the ability of aquatic biogeochemists to examine molecular-level components of complex mixtures of organic matter. The ability to accurately assess the chemical composition, elemental formulas, or both of detected compounds is critical to these studies. Here we build on previous work that uses functional group relationships between compounds to extend elemental formulas of low molecular weight compounds to those of higher molecular weight. We propose an automated compound identification algorithm (CIA) for the analysis of ultra-high-resolution mass spectra of natural organic matter acquired by electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry. This approach is benchmarked with synthetic data sets of compounds cited in the literature. The sensitivity of our results is examined for different sources of error, and CIA is applied to two previously published data sets. We find that CIA works well for data sets with high mass accuracy (<1 ppm) and can accurately determine the elemental formulas for >95% of all compounds composed of C, H, O, and N. Data with lower mass accuracy must be accompanied with additional knowledge of chemical structure, composition, or both in order to yield accurate elemental formulas.