Lars Rüpke

Deep Roots of the Messinian Salinity Crisis

The Messinian salinity crisis—the desiccation of the Mediterranean Sea between 5.96 and 5.33 million years (Myr) ago—was one of the most dramatic events on Earth during the Cenozoic era. It resulted from the closure of marine gateways between the Atlantic Ocean and the Mediterranean Sea, the causes of which remain enigmatic. Here we use the age and composition of volcanic rocks to reconstruct the geodynamic evolution of the westernmost Mediterranean from the Middle Miocene epoch to the Pleistocene epoch (about 12.1–0.65 Myr ago). Our data show that a marked shift in the geochemistry of mantle-derived volcanic rocks, reflecting a change from subduction-related to intraplate-type volcanism, occurred between 6.3 and 4.8 Myr ago, largely synchronous with the Messinian salinity crisis. Using a thermomechanical model, we show that westward roll back of subducted Tethys oceanic lithosphere and associated asthenospheric upwelling provides a plausible mechanism for producing the shift in magma chemistry and the necessary uplift (∼1 km) along the African and Iberian continental margins to close the Miocene marine gateways, thereby causing the Messinian salinity crisis.

Are the regional variations in Central American arc lavas due to differing basaltic versus peridotitic slab sources of fluids

Central American arc volcanism shows strong regional trends in lava chemistry that result from differing slab contributions to arc melting. However, the mechanism that transfers slab-derived trace elements into the mantle wedge remains largely unknown. By using a dynamic model for mantle flow and fluid release, we model the fate of three different slab-fluid sources: sediment, ocean crust, and serpentinized mantle. In the open subarc system, sediments lose almost all their highly fluid mobile elements by ∼50 km depth, so other fluid sources are necessary to explain the slab signal in arc-lava compositions. The well-documented transition from lavas with a strong geochemical slab signature (i.e., high Ba/La ratios) found in Nicaragua to lavas with a weaker slab signature (i.e., low Ba/La ratios) erupted in Costa Rica seems easiest to produce by a higher fraction of serpentine-hosted fluids released from the deeply faulted, highly serpentinized lithosphere subducting beneath Nicaragua than from the less deeply faulted, thicker, amphibolitic oceanic-crust and oceanic-plateau lithosphere subducting beneath Costa Rica.

Rising Arctic Ocean temperatures cause gas hydrate destabilization and ocean acidification

Vast amounts of methane hydrates are potentially stored in sediments along the continental margins, owing their stability to low temperature – high pressure conditions. Global warming could destabilize these hydrates and cause a release of methane (CH 4) into the water column and possibly the atmosphere. Since the Arctic has and will be warmed considerably, Arctic bottom water temperatures and their future evolution projected by a climate model were analyzed. The resulting warming is spatially inhomogeneous, with the strongest impact on shallow regions affected by Atlantic inflow. Within the next 100 years, the warming affects 25% of shallow and mid-depth regions containing methane hydrates. Release of methane from melting hydrates in these areas could enhance ocean acidification and oxygen depletion in the water column. The impact of methane release on global warming, however, would not be significant within the considered time span.

Modeling the fate of methane hydrates under global warming

Large amounts of methane hydrate locked up within marine sediments are vulnerable to climate change. Changes in bottom water temperatures may lead to their destabilization and the release of methane into the water column or even the atmosphere. In a multimodel approach, the possible impact of destabilizing methane hydrates onto global climate within the next century is evaluated. The focus is set on changing bottom water temperatures to infer the response of the global methane hydrate inventory to future climate change. Present and future bottom water temperatures are evaluated by the combined use of hindcast high-resolution ocean circulation simulations and climate modeling for the next century. The changing global hydrate inventory is computed using the parameterized transfer function recently proposed by Wallmann et al. (2012). We find that the present-day worlds total marine methane hydrate inventory is estimated to be 1146Gt of methane carbon. Within the next 100years this global inventory may be reduced by ∼0.03% (releasing ∼473Mt methane from the seafloor). Compared to the present-day annual emissions of anthropogenic methane, the amount of methane released from melting hydrates by 2100 is small and will not have a major impact on the global climate. On a regional scale, ocean bottom warming over the next 100years will result in a relatively large decrease in the methane hydrate deposits, with the Arctic and Blake Ridge region, offshore South Carolina, being most affected.

Hybrid shallow on-axis and deep off-axis hydrothermal circulation at fast-spreading ridges

Hydrothermal flow at oceanic spreading centres accounts for about ten per cent of all heat flux in the oceans and controls the thermal structure of young oceanic plates. It also influences ocean and crustal chemistry, provides a basis for chemosynthetic ecosystems, and has formed massive sulphide ore deposits throughout Earth’s history. Despite this, how and under what conditions heat is extracted, in particular from the lower crust, remains largely unclear. Here we present high-resolution, whole-crust, two- and three-dimensional simulations of hydrothermal flow beneath fast-spreading ridges that predict the existence of two interacting flow components, controlled by different physical mechanisms, that merge above the melt lens to feed ridge-centred vent sites. Shallow on-axis flow structures develop owing to the thermodynamic properties of water, whereas deeper off-axis flow is strongly shaped by crustal permeability, particularly the brittle–ductile transition. About 60 per cent of the discharging fluid mass is replenished on-axis by warm (up to 300 degrees Celsius) recharge flow surrounding the hot thermal plumes, and the remaining 40 per cent or so occurs as colder and broader recharge up to several kilometres away from the axis that feeds hot (500–700 degrees Celsius) deep-rooted off-axis flow towards the ridge. Despite its lower contribution to the total mass flux, this deep off-axis flow carries about 70 per cent of the thermal energy released at the ridge axis. This combination of two flow components explains the seismically determined thermal structure of the crust and reconciles previously incompatible models favouring either shallower on-axis or deeper off-axis hydrothermal circulation.

Automated thermotectonostratigraphic basin reconstruction: Viking Graben case study

We present a generic algorithm for automating sedimentary basin reconstruction. Automation is achieved through the coupling of a two-dimensional thermotectonostratigraphic forward model to an inverse scheme that updates the model parameters until the input stratigraphy is fitted to a desired accuracy. The forward model solves for lithospheric thinning, flexural isostasy, sediment deposition, and transient heat flow. The inverse model updates the crustal- and mantle-thinning factors and paleowater depth. Both models combined allow for automated forward modeling of the structural and thermal evolution of extensional sedimentary basins. The potential and robustness of this method is demonstrated through a reconstruction case study of the northern Viking Graben in the North Sea. This reconstruction fits present stratigraphy, borehole temperatures, vitrinite reflectance data, and paleowater depth. The predictive power of the model is illustrated through the successful identification of possible targets along the transect, where the principal source rocks are in the oil and gas windows. These locations coincide well with known oil and gas occurrences. The key benefits of the presented algorithm are as follows: (1) only standard input data are required, (2) crustal- and mantle-thinning factors and paleowater depth are automatically computed, and (3) sedimentary basin reconstruction is greatly facilitated and can thus be more easily integrated into basin analysis and exploration risk assessment.

Modeling fluid flow in sedimentary basins with sill intrusions: Implications for hydrothermal venting and climate change

Large volumes of magma emplaced within sedimentary basins have been linked to multiple climate change events due to release of greenhouse gases such as CH4. Basin-scale estimates of thermogenic methane generation show that this process alone could generate enough greenhouse gases to trigger global incidents. However, the rates at which these gases are transported and released into the atmosphere are quantitatively unknown. We use a 2D, hybrid FEM/FVM model that solves for fully compressible fluid flow to quantify the thermogenic release and transport of methane and to evaluate flow patterns within these systems. Our results show that the methane generation potential in systems with fluid flow does not significantly differ from that estimated in diffusive systems. The values diverge when vigorous convection occurs with a maximum variation of about 50%. The fluid migration pattern around a cooling, impermeable sill alone generates hydrothermal plumes without the need for other processes such as boiling and/or explosive degassing. These fluid pathways are rooted at the edges of the outer sills consistent with seismic imaging. Methane venting at the surface occurs in three distinct stages and can last for hundreds of thousands of years. Our simulations suggest that although the quantity of methane potentially generated within the contact aureole can cause catastrophic climate change, the rate at which this methane is released into the atmosphere is too slow to trigger, by itself, some of the negative δ13C excursions observed in the fossil record over short time scales (< 10,000 years).

Basin modelling of a transform margin setting: structural, thermal and hydrocarbon evolution of the Tano Basin, Ghana

ABSTRACT This study explores the structural and thermal evolution of the Ghana transform margin. The main objective is to explore how the opening of the Atlantic Ocean and subsequent interaction with the Mid-Atlantic Ridge (MAR) has affected the margins structural and thermal evolution. Two representative evolution scenarios are described: a reference case that neglects the influence of continental breakup and a second scenario that accounts for a possible heat influx during the passage of the MAR as well as magmatic underplating. These two scenarios have further been analysed for the implications for the hydrocarbon potential of the region. The scenario analysis builds on a suite of 2D realizations performed with TECMOD2D, modelling software for automated basin reconstructions. As the observed stratigraphy is input, the structural and thermal evolution of the basin is automatically reconstructed. This is achieved through the coupling of a lithosphere scale forward model with an inverse algorithm for model parameter optimization. We find that lateral heat transport from the passing MAR in combination with flexure of the lithosphere can explain the observed uplift of the margin. These results were obtained for a broken plate elasticity solution with a relative large value for the effective elastic thickness (Te=15) and necking level (15 km). Lateral heat flow from oceanic lithosphere is clearly visible in elevated basement heat flow values up to 50 km away from the ocean–continent transition (OCT). This influx of heat does not seem to have affected the maturation history along the margin significantly. Only the deepest sediments close to the OCT show slightly elevated vitrinite reflectance in simulations that account for the passage of the MAR. In conclusion, it appears that that lateral heat transport from the oceanic lithosphere is instrumental in shaping the Ghana transform margin but seems to have only limited control on the maturation history.

Modeling of craton stability using a viscoelastic rheology

Archean cratons belong to the most remarkable features of our planet since they represent continental crust that has avoided reworking for several billions of years. Even more, it has become evident from both geophysical and petrological studies that cratons exhibit deep lithospheric keels which equally remained stable ever since the formation of the cratons in the Archean. Dating of inclusions in diamonds from kimberlite pipes gives Archean ages, suggesting that the Archean lithosphere must have been cold soon after its formation in the Archean (in order to allow for the existence of diamonds) and must have stayed in that state ever since. Yet, although strong evidence for the thermal stability of Archean cratonic lithosphere for billions of years is provided by diamond dating, the long-term thermal stability of cratonic keels was questioned on the basis of numerical modeling results. We devised a viscoelastic mantle convection model for exploring cratonic stability in the stagnant lid regime. Our modeling results indicate that within the limitations of the stagnant lid approach, the application of a sufficiently high temperature-dependent viscosity ratio can provide for thermal craton stability for billions of years. The comparison between simulations with viscous and viscoelastic rheology indicates no significant influence of elasticity on craton stability. Yet, a viscoelastic rheology provides a physical transition from viscously to elastically dominated regimes within the keel, thus rendering introduction of arbitrary viscosity cutoffs, as employed in viscous models, unnecessary.

Controls of bathymetric relief on hydrothermal fluid flow at mid‐ocean ridges

We present quantitative modeling results for the effects of surface relief on hydrothermal convection at ocean-spreading centers investigating how vent site locations and subsurface flow patterns are affected by bathymetry induced sub-seafloor pressure variations. The model is based on a 2-D FEM solver for fluid flow in porous media and is used to simulate hydrothermal convection systematically in 375 synthetic studies. The results of these studies show that bathymetric relief has a profound effect on hydrothermal flow: bathymetric highs induce subsurface pressure variations that can deviate upwelling zones and favor venting at structural highs. The deviation angle from vertical upwelling can be expressed by a single linear dependence relating deviation angle to bathymetric slope and depth of the heat source. These findings are confirmed in two case studies for the East Pacific Rise at 9°30′N and Lucky Strike hydrothermal fields. In both cases, it is possible to predict the observed vent field locations only if bathymetry is taken into account. Our results thereby show that bathymetric relief should be considered in simulations of submarine hydrothermal systems and plays a key role especially in focusing venting of across axis hydrothermal flow onto the ridge axis of fast spreading ridges.