A Coupled Solid-Fluid Method for Modeling Subduction
International audience; We present a novel dynamic approach for solid/fluid coupling by joining two different numerical methods: Boundary Element Method (BEM) and Finite Element Method (FEM). FEM results describe the thermo-mechanical evolution of the solid while the fluid is solved with the BEM. The bidirectional feedback between the two domains evolves along a Lagrangian interface where the FEM domain is embedded inside the BEM domain. The feedback between the two codes is based on the calculation of a specific drag tensor for each boundary/finite element. The approach is presented here to solve the complex problem of the descent of a cold subducting oceanic plate into a hot fluid like ma…
A smeared seismicity constitutive model
The classical application of rate and state dependent frictional constitutive laws has involved the instabilities developed between two sliding surfaces. In such a situation, the behaviour and evolution of asperities is the controlling mechanism of velocity weakening. However, most faults have a substantial thickness and it would appear that it is the bulk behaviour of the fault gouge, at whatever scale, that is important. The purpose of this paper is to explore how bulk frictional sliding behaviour may be described. We explore here the consequences of applying the rate and state framework initially developed to describe the frictional behaviour at the interface between two interacting slid…
Fluid reservoirs in the crust and mechanical coupling between the upper and lower crust
An important observation associated with seismic activity on the Nagamachi-Rifu Fault is the existence of tabular, fluid rich zones at mid-crustal levels. These zones resemble the “bright spots” seen in many seismic images of the crust worldwide. The aim of this paper is to develop the mechanical foundations for the formation of such zones. To do so requires an understanding of the distribution of pore fluid pressure in a deforming crust. In a hydrostatically stressed porous material, the pore fluid pressure should equal the mean stress in order to keep the pores from collapsing. Past discussions of this subject imply very high pore fluid pressures, two to three times lithostatic. Considera…
Water and Geodynamics
Hydrogen is the most abundant element (Fig. 1⇓) in the galaxy and our solar system (Lodders 2003). Therefore it is not astonishing that hydrogen is a key player in the geodynamic evolution of planets. Its fate in the early Earth, after condensation of the solar nebula, the accretion of our planet and hydrogen reprocessing through early asteroidal and cometary bombardment (Dauphas et al. 2000) and segregation of a proto-Earth into iron core and silicate mantle is described elsewhere in this volume (Marty and Yokochi 2006). Figure 1. Abundance of elements in the solar system in numbers of atoms per 106 atoms of silicon (Lodders 2003). This chapter concerns itself with the geodynamics of the m…
Heat generation associated with collision of two plates: the Himalayan geothermal belt
Abstract An analysis of heat discharged by thermal fluids along the c. 3000-km-long Himalayan geothermal belt (HGB) shows that heat transfer is concentrated along 30- to 50-km-wide `heat bands' which are associated with at least 600 geothermal systems. The bands have been interpreted as segments of major, concentric slip lines caused by plastic deformation of the ductile crust within the Asian plate resulting from plate collision. Assuming that this crust behaves like an ideal plastic medium, the heat transfer within and along a slip line can be estimated. It amounts to c. 55 mW/m2 for a 40-km-wide band. Estimates of present-day heat discharges point to 20–35 mW/m2 for convective, and 10–30…
Continental extension: From core complexes to rigid block faulting
Extension of overthickened continental crust is commonly characterized by an early core complex stage of extension followed by a later stage of crustal-scale rigid block faulting. These two stages are clearly recognized during the extensional destruction of the Alpine orogen in northeast Corsica, where rigid block faulting overprinting core complex formation eventually led to crustal separation and the formation of a new oceanic backarc basin (the Ligurian Sea). Here we investigate the geodynamic evolution of continental extension by using a novel, fully coupled thermomechanical numerical model of the continental crust. We consider that the dynamic evolution is governed by fault weakening, …
Analysis of thermally induced flows in the laboratory by geoelectrical 3-D tomography
[1] Many natural bodies as well as materials inside industrial installations, such as the Earth's mantle and the glass inside melting furnaces, exchange matter through convection. These processes result from differences in temperature, density, and chemical concentration. In this analysis, we focus on the visualization of thermally driven flows in the laboratory. In nature and in industrial installations, it is difficult to measure the temperature inside the object of interest directly. We benchmark a new DC-geoelectrical 3-D tomography method for temperature measurements that allows obtaining temperature values without influencing the flow pattern. For verification of the method, we use di…
Dilatant plasticity applied to Alpine collision: ductile void growth in the intraplate area beneath the Eifel volcanic field
The Eifel is located in the middle of the European plate far away from any active plate boundary, yet it appears to be a maximum of intraplate tectonic activity. A map of intraplate seismic energy flow shows that the Eifel is linked to the Alpine collisional belt via a narrow seismoactive shear zone. Two parallel Quaternary volcanic zones (the East Eifel Volcanic Zone EEVZ and the West Eifel Volcanic Zone WEVZ) line up with the seismogenic shear zone. Xenoliths ejected from these volcanic lineaments indicate upper mantle shearing by dynamic recrystallization textures and metasomatized chemistry. Important CO2-dominated mantle degassing observed in mineral springs, lakes or dry degassing sug…
The effect of energy feedbacks on continental strength
The classical strength profile of continents is derived from a quasi-static view of their rheological response to stress--one that does not consider dynamic interactions between brittle and ductile layers. Such interactions result in complexities of failure in the brittle-ductile transition and the need to couple energy to understand strain localization. Here we investigate continental deformation by solving the fully coupled energy, momentum and continuum equations. We show that this approach produces unexpected feedback processes, leading to a significantly weaker dynamic strength evolution. In our model, stress localization focused on the brittle-ductile transition leads to the spontaneo…
Rapid conversion of elastic energy into plastic shear heating during incipient necking of the lithosphere
An important and novel mechanism for ductile failure of the lithosphere is identified here, which is intrinsic to the thermal-mechanical feedback in a temperature dependent plastic body with coupled elastic fields. Both a temperature-dependent power-law visco-elasto-plastic rheology and a temperature-dependent elasto-plastic rheology are employed to study in a self-consistent fashion the deformation of the lithosphere subject to extension by means of a two-dimensional, finite-element code. A structural perturbation initially localizes elasto-plastic deformation only in its immediate vicinity. However, after 800,000 years have elapsed the localized zone of deformation takes off in a ‘crack-l…
On the thermodynamics of listric faults
We investigate a novel fully coupled thermal-mechanical numerical model of the crust in order to trace the physics of interaction of its brittle and ductile layers. In a unified approach these layers develop in a natural transition as a function of the state variables pressure, deviatoric stress, temperature and strain-rate. We find that the main storage of elastic energy lies in the domain where brittle and ductile strain-rates overlap so that shear zones are attracted to this zone of maximum energy dissipation. This dissipation appears as a local heat source (shear heating). The brittle-ductile transition zone evolves through extreme weakening by thermo-mechanical feedback. The physics of…