Behn Mark D.

No Thumbnail Available
Last Name
Behn
First Name
Mark D.
ORCID
0000-0002-2001-1335

Search Results

Now showing 1 - 7 of 7
  • Preprint
    Fracture propagation to the base of the Greenland Ice Sheet during supraglacial lake drainage
    ( 2008-02-20) Das, Sarah B. ; Joughin, Ian ; Behn, Mark D. ; Howat, Ian M. ; King, Matt A. ; Lizarralde, Daniel ; Bhatia, Maya P.
    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 m 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 following 24 hours. The short-lived dynamic response suggests an efficient drainage system dispersed the meltwater subglacially. The integrated effect of multiple lake drainages could explain the observed net regional summer ice speedup.
  • Preprint
    Mid-ocean ridge jumps associated with hotspot magmatism
    ( 2007-10-11) Mittelstaedt, Eric ; Ito, Garrett T. ; Behn, Mark D.
    Hotspot-ridge interaction produces a wide range of phenomena including excess crustal thickness, geochemical anomalies, off-axis volcanic ridges and ridge relocations or jumps. Ridges are recorded to have jumped toward many hotspots including, Iceland, Discovery, Galapagos, Kerguelen and Tristan de Cuhna. The causes of ridge jumps likely involve a number of interacting processes related to hotspots. One such process is reheating of the lithosphere as magma penetrates it to feed near-axis volcanism. We study this effect by using the hybrid, finite-element code, FLAC, to simulate two-dimensional (2-D, cross-section) viscous mantle flow, elasto-plastic deformation of the lithosphere and heat transport in a ridge setting near an off-axis hotspot. Heating due to magma transport through the lithosphere is implemented within a hotspot region of fixed width. To determine the conditions necessary to initiate a ridge jump, we vary four parameters: hotspot magmatic heating rate, spreading rate, seafloor age at the location of the hotspot and ridge migration rate. Our results indicate that the hotspot magmatic heating rate required to initiate a ridge jump increases non-linearly with increasing spreading rate and seafloor age. Models predict that magmatic heating, itself, is most likely to cause jumps at slow spreading rates such as at the Mid-Atlantic Ridge on Iceland. In contrast, despite the higher magma flux at the Galapagos hotspot, magmatic heating alone is probably insufficient to induce a ridge jump at the present-day due to the intermediate ridge spreading rate of the Galapagos Spreading Center. The time required to achieve a ridge jump, for fixed or migrating ridges, is found to be on the order of 105-106 years. Simulations that incorporate ridge migration predict that after a ridge jump occurs the hotspot and ridge migrate together for time periods that increase with magma flux. Model results also suggest a mechanism for ridge reorganizations not related to hotspots such as ridge jumps in back-arc settings and ridge segment propagation along the Mid-Atlantic Ridge.
  • Preprint
    Topographic controls on dike injection in volcanic rift zones
    ( 2006-04-03) Behn, Mark D. ; Buck, W. Roger ; Sacks, I. Selwyn
    Dike emplacement in volcanic rift zones is often associated with the injection of “bladelike” dikes, which propagate long distances parallel to the rift, but frequently remain trapped at depth and erupt only near the tip of the dike. Over geologic time, this style of dike injection implies that a greater percentage of extension is accommodated by magma accretion at depth than near the surface. In this study, we investigate the evolution of faulting, topography, and stress state in volcanic rift zones using a kinematic model for dike injection in an extending 2-D elastic-viscoplastic layer. We show that the intrusion of blade-like dikes focuses deformation at the rift axis, leading to the formation of an axial rift valley. However, flexure associated with the development of the rift topography generates compression at the base of the plate. If the magnitude of these deviatoric compressive stresses exceeds the deviatoric tensile stress associated with far-field extension, further dike injection will be inhibited. In general, this transition from tensile to compressive deviatoric stresses occurs when the rate of accretion in the lower crust is greater than 50-60% of the far-field extension rate. These results indicate that over geologic time-scales the injection of blade-like dikes is a self-limiting process in which dike-generated faulting and topography result in an efficient feedback mechanism that controls the time-averaged distribution of magma accretion within the crust.
  • Preprint
    Spatio-temporal evolution of strain accumulation derived from multi-scale observations of Late Jurassic rifting in the northern North Sea : a critical test of models for lithospheric extension
    ( 2005-01-23) Cowie, Patience A. ; Underhill, John R. ; Behn, Mark D. ; Lin, Jian ; Gill, Caroline E.
    We integrate observations of lithospheric extension over a wide range of spatial and temporal scales within the northern North Sea basin and critically review the extent to which existing theories of lithospheric deformation can account for these observations. Data obtained through a prolonged period of hydrocarbon exploration and production has yielded a dense and diverse data set over the entire Viking Graben and its flanking platform areas. These data show how syn-rift accommodation within the basin varied in space and time with sub-kilometer-scale spatial resolution and a temporal resolution of 2–3 Myr. Regional interpretations of 2D seismic reflection, refraction and gravity data for this area have also been published and provide an image of total basin wide stretching for the entire crust. These image data are combined with published strain rate inversion results obtained from tectonic subsidence patterns to constrain the spatio-temporal evolution of strain accumulation throughout the lithosphere during the 40 Myr (170–130 Ma) period of Late Jurassic extension across this basin. For the first 25–30 Myr, strain localisation dominated basin development with strain rates at the eventual rift axis increasing while strain rates over the flanking areas declined. As strain rates across the whole basin were consistently very low (< 3 × 10- 16 s- l), thermally induced strength loss cannot explain this phenomenon. The strain localisation is manifest in the near-surface by a systematic migration of fault activity. The pattern and timing of this migration are inconsistent with flexural bending stresses exerting an underlying control, especially when estimates of flexural rigidity for this area are considered. The best explanation for what is observed in this time period is a coupling between near-surface strain localisation, driven by brittle (or plastic) failure, and the evolving thermal structure of the lithosphere. We demonstrate this process using a continuum mechanics model for normal fault growth that incorporates the strain rate-dependence of frictional strength observed in laboratory studies. During the final 10 Myr of basin formation, strain accumulation was focused within the axis and strain rates declined rapidly. Replacement of weak crust by stronger mantle material plus crustal buoyancy forces can adequately explain this decline.
  • Preprint
    Diapirs as the source of the sediment signature in arc lavas
    ( 2011-05-31) Behn, Mark D. ; Kelemen, Peter B. ; Hirth, Greg ; Hacker, Bradley R. ; Massonne, Hans-Joachim
    Many arc lavas show evidence for the involvement of subducted sediment in the melting process. There is debate whether this “sediment melt” signature forms at relatively low temperature near the fluid-saturated solidus or at higher temperature beyond the breakdown of trace-element-rich accessory minerals. We present new geochemical data from high- to ultrahigh-pressure rocks that underwent subduction and show no significant depletion of key trace elements in the sediment melt component until peak metamorphic temperatures exceeded ~1050ºC from 2.7 to 5 GPa. These temperatures are higher than for the top of the subducting plate at similar pressures based on thermal models. To address this discrepancy, we use instability calculations for a non-Newtonian buoyant layer in a viscous half-space to show that, in typical subduction zones, solid-state sediment diapirs initiate at temperatures between 500–850ºC. Based on these calculations, we propose that the sediment melt component in arc magmas is produced by high degrees of dehydration melting in buoyant diapirs of metasediment that detach from the slab and rise into the hot mantle wedge. Efficient recycling of sediments into the wedge by this mechanism will alter volatile fluxes into the deep mantle compared to estimates based solely on devolatilization of the slab.
  • Preprint
    Sensitivity of seafloor bathymetry to climate-driven fluctuations in mid-ocean ridge magma supply
    ( 2015-09) Olive, Jean-Arthur L. ; Behn, Mark D. ; Ito, Garrett T. ; Buck, W. Roger ; Escartin, Javier E. ; Howell, Samuel M.
    Recent studies have proposed that the bathymetric fabric of the seafloor formed at mid-ocean ridges records rapid (23–100 kyr) fluctuations in ridge magma supply caused by sea level changes that modulate melt production in the underlying mantle. Using quantitative models of faulting and magma emplacement, we demonstrate that, in fact, seafloor-shaping processes act as a low-pass filter on variations in magma supply, strongly damping fluctuations shorter than ~100 kyr. We show that the systematic decrease in dominant seafloor wavelengths with increasing spreading rate is best explained by a model of fault growth and abandonment under a steady magma input. This provides a robust framework for deciphering the footprint of mantle melting in the fabric of abyssal hills, the most common topographic feature on Earth.
  • Preprint
    Implications of grain size evolution on the seismic structure of the oceanic upper mantle
    ( 2009-03-04) Behn, Mark D. ; Hirth, Greg ; Elsenbeck, James R.
    We construct a 1-D steady-state channel flow model for grain size evolution in the oceanic upper mantle using a composite diffusion-dislocation creep rheology. Grain size evolution is calculated assuming that grain size is controlled by a competition between dynamic recrystallization and grain growth. Applying this grain size evolution model to the oceanic upper mantle we calculate grain size as a function of depth, seafloor age, and mantle water content. The resulting grain size structure is used to predict shear wave velocity (VS) and seismic quality factor (Q). For a plate age of 60 Myr and an olivine water content of 1000 H/106Si, we find that grain size reaches a minimum of ~15 mm at ~150 km depth and then increases to ~20–30 mm at a depth of 400 km. This grain size structure produces a good fit to the low seismic shear wave velocity zone (LVZ) in oceanic upper mantle observed by surface wave studies assuming that the influence of hydrogen on anelastic behavior is similar to that observed for steady state creep. Further it predicts a viscosity of ~1019 Pa s at 150 km depth and dislocation creep to be the dominant deformation mechanism throughout the oceanic upper mantle, consistent with geophysical observations. We predict larger grain sizes than proposed in recent studies, in which the LVZ was explained by a dry mantle and a minimum grain size of 1 mm. However, we show that for a 1 mm grain size, diffusion creep is the dominant deformation mechanism above 100– 200 km depth, inconsistent with abundant observations of seismic anisotropy from surface wave studies. We therefore conclude that a combination of grain size evolution and a hydrated upper mantle is the most likely explanation for both the isotropic and anisotropic seismic structure of the oceanic upper mantle. Our results also suggest that melt extraction from the mantle will be significantly more efficient than predicted in previous modeling studies that assumed grain sizes of ~1 mm.