scholarly journals Shear heating reconciles thermal models with the metamorphic rock record of subduction

2018 ◽  
Vol 115 (46) ◽  
pp. 11706-11711 ◽  
Author(s):  
Matthew J. Kohn ◽  
Adrian E. Castro ◽  
Buchanan C. Kerswell ◽  
César R. Ranero ◽  
Frank S. Spear
Keyword(s):  
Heat Flow ◽  
Subduction Zones ◽  
Thermal Structure ◽  
Metamorphic Rocks ◽  
Shear Stresses ◽  
Thermal Models ◽  
Average Value ◽  
Temperature Conditions ◽  
Mechanical Models ◽  
Shear Heating

Some commonly referenced thermal-mechanical models of current subduction zones imply temperatures that are 100–500 °C colder at 30–80-km depth than pressure–temperature conditions determined thermobarometrically from exhumed metamorphic rocks. Accurately inferring subduction zone thermal structure, whether from models or rocks, is crucial for predicting metamorphic reactions and associated fluid release, subarc melting conditions, rheologies, and fault-slip phenomena. Here, we compile surface heat flow data from subduction zones worldwide and show that values are higher than can be explained for a frictionless subduction interface often assumed for modeling. An additional heat source––likely shear heating––is required to explain these forearc heat flow values. A friction coefficient of at least 0.03 and possibly as high as 0.1 in some cases explains these data, and we recommend a provisional average value of 0.05 ± 0.015 for modeling. Even small coefficients of friction can contribute several hundred degrees of heating at depths of 30–80 km. Adding such shear stresses to thermal models quantitatively reproduces the pressure–temperature conditions recorded by exhumed metamorphic rocks. Comparatively higher temperatures generally drive rock dehydration and densification, so, at a given depth, hotter rocks are denser than colder rocks, and harder to exhume through buoyancy mechanisms. Consequently––conversely to previous proposals––exhumed metamorphic rocks might overrepresent old-cold subduction where rocks at the slab interface are wetter and more buoyant than in young-hot subduction zones.

scholarly journals The thermal structure of Israel

2011 ◽  
Vol 3 (1) ◽  
pp. 431-452 ◽  
Author(s):  
E. Shalev ◽  
V. Lyakhovsky ◽  
Y. Weinstein ◽  
Z. Ben-Avraham
Keyword(s):  
Heat Flux ◽  
Thermal Structure ◽  
Local Heating ◽  
Arabian Shield ◽  
Supporting Evidence ◽  
Hydrologic Models ◽  
Average Value ◽  
Water Wells ◽  
Shear Heating ◽  
Deep Seismicity

Abstract. Heat flux at the Arabian Shield is a significant component in reconstructing tectonic, seismic, and hydrologic models. In this paper we analyze temperature data from all the available oil and water wells in Israel. We show that the average heat flux in Israel is 40–45 mW m−2. A supporting evidence for the low heat flux is the relatively deep seismicity, extending almost to the mantle in the region. A Heat flux anomaly that exists in Northern Israel and Jordan could be attributed to groundwater flow or young magmatic activity (~100 000 years) that is common in this area. Xenoliths that yield relatively steep geothermal gradients could be the result of local heating by magmas or by lithospheric necking and shear heating. The higher Heat flux in Southern Israel and Jordan probably reflects the opening of the Red Sea and the Gulf of Eilat and does not reflect the average value of the Arabian Shield.

scholarly journals Deep decoupling in subduction zones: Observations and temperature limits

Geosphere ◽  
2020 ◽  
Vol 16 (6) ◽  
pp. 1408-1424 ◽  
Author(s):  
Geoffrey A. Abers ◽  
Peter E. van Keken ◽  
Cian R. Wilson
Keyword(s):  
Subduction Zone ◽  
Shear Zone ◽  
Subduction Zones ◽  
Thermal Structure ◽  
Shear Zones ◽  
Seismic Attenuation ◽  
Plate Interface ◽  
Brittle Ductile Transition ◽  
Wide Range ◽  
Shear Heating

Abstract The plate interface undergoes two transitions between seismogenic depths and subarc depths. A brittle-ductile transition at 20–50 km depth is followed by a transition to full viscous coupling to the overlying mantle wedge at ∼80 km depth. We review evidence for both transitions, focusing on heat-flow and seismic-attenuation constraints on the deeper transition. The intervening ductile shear zone likely weakens considerably as temperature increases, such that its rheology exerts a stronger control on subduction-zone thermal structure than does frictional shear heating. We evaluate its role through analytic approximations and two-dimensional finite-element models for both idealized subduction geometries and those resembling real subduction zones. We show that a temperature-buffering process exists in the shear zone that results in temperatures being tightly controlled by the rheological strength of that shear zone’s material for a wide range of shear-heating behaviors of the shallower brittle region. Higher temperatures result in weaker shear zones and hence less heat generation, so temperatures stop increasing and shear zones stop weakening. The net result for many rheologies are temperatures limited to ≤350–420 °C along the plate interface below the cold forearc of most subduction zones until the hot coupled mantle is approached. Very young incoming plates are the exception. This rheological buffering desensitizes subduction-zone thermal structure to many parameters and may help explain the global constancy of the 80 km coupling limit. We recalculate water fluxes to the forearc wedge and deep mantle and find that shear heating has little effect on global water circulation.

scholarly journals Bayesian inversion of surface heat flow in subduction zones: a framework to refine geodynamic models based on observational constraints

2020 ◽  
Vol 222 (1) ◽  
pp. 103-109
Author(s):  
M Morishige ◽  
T Kuwatani
Keyword(s):  
Heat Flow ◽  
Subduction Zones ◽  
Thermal Structure ◽  
Surface Heat ◽  
Dislocation Creep ◽  
Model Parameters ◽  
Limited Information ◽  
Surface Heat Flow ◽  
Trade Offs ◽  
Geodynamic Models

SUMMARY Surface heat flow has been widely used to constrain the thermal structure of subduction zones. However, the forward modelling approaches in previous geodynamic studies have only provided limited information on the model parameters controlling the thermal structure, which makes model validation difficult. Here we apply a probabilistic inversion technique based on Bayes’ theorem to surface heat flow data from Tohoku in Japan and Cascadia to simultaneously infer five model parameters that appear to have the greatest influence on the thermal structure of subduction zones. The surface heat flow is predicted via 2-D steady-state thermomechanical modelling. The Metropolis algorithm is used to obtain the posterior probability distributions. A comparison of our results with previous estimates indicates that our activation energy for the shear viscosity of dislocation creep is lower in both regions, and our radiogenic heat production rate in the upper continental crust is lower in Cascadia. These findings suggest that our geodynamic models cannot explain the surface heat flow observations with the acceptable ranges of model parameter values. We therefore need to refine the models by including, for example, the effects of recent backarc extension, vigorous thermal convection beneath the overriding plate and fluid circulation in the uppermost part of the oceanic crust. The approach presented here also allows us to determine trade-offs between the parameters. This study provides a framework to validate and refine geodynamic models based on various types of observations.

scholarly journals Advances in the thermal and petrologic modeling of subduction zones

Geosphere ◽  
2020 ◽  
Vol 16 (4) ◽  
pp. 936-952 ◽  
Author(s):  
Simon M. Peacock
Keyword(s):  
Heat Flux ◽  
Subduction Zones ◽  
Mantle Wedge ◽  
Dynamic Models ◽  
Seismic Velocity ◽  
Thermal Structure ◽  
Metamorphic Rocks ◽  
Ridge Subduction ◽  
Dehydration Reactions ◽  
Flux Measurements

Abstract In the two decades since Subduction: Top to Bottom was published in 1996, improved analytical and numerical thermal-petrologic models of subduction zones have been constructed and evaluated against new seismological and geological observations. Advances in thermal modeling include a range of new approaches to incorporating shear (frictional, viscous) heating along the subduction interface and to simulating induced flow in the mantle wedge. Forearc heat-flux measurements constrain the apparent coefficient of friction (μ′) along the plate interface to <∼0.1, but the extent to which μ′ may vary between subduction zones remains challenging to discern owing to scatter in the heat-flux measurements and uncertainties in the magnitude and distribution of radiogenic heat production in the overriding crust. Flow in the mantle wedge and the resulting thermal structure depend on the rheology of variably hydrated mantle rocks and the depth at which the subducting slab becomes coupled to the overlying mantle wedge. Advances in petrologic modeling include the incorporation of sophisticated thermodynamic software packages into thermal models and the prediction of seismic velocities from mineralogic and petrologic models. Current thermal-petrologic models show very good agreement between the predicted location of metamorphic dehydration reactions and observed intermediate-depth earthquakes, and between the predicted location of the basalt-to-eclogite transition in subducting oceanic crust and observed landward-dipping, low-seismic-velocity layers. Exhumed high-pressure, low-temperature metamorphic rocks provide insight into subduction-zone temperatures, but important thermal parameters (e.g., convergence rate) are not well constrained, and metamorphic rocks exposed at the surface today may reflect relatively warm conditions in the past associated with subduction initiation or ridge subduction. We can anticipate additional advances in our understanding of subduction zones as a result of further testing of model predictions against geologic and geophysical observations, and of evaluating the importance of advective processes, such as diapirism and subduction-channel flow, that are not captured in hybrid kinematic-dynamic models of subduction zones but are observed in fully dynamical models under certain conditions.

scholarly journals Anisotropic Differences in the Thermal Conductivity of Rocks: A Summary from Core Measurement Data in East China

Minerals ◽  
2021 ◽  
Vol 11 (10) ◽  
pp. 1135
Author(s):  
Yibo Wang ◽  
Zhuting Wang ◽  
Lin Shi ◽  
Yuwei Rong ◽  
Jie Hu ◽  
...  
Keyword(s):  
Thermal Conductivity ◽  
Heat Flow ◽  
Thermal Structure ◽  
Eastern China ◽  
Measurement Data ◽  
East China ◽  
Metamorphic Rocks ◽  
Strong Anisotropy ◽  
Conductivity Anisotropy ◽  
Development And Utilization

The study of thermal conductivity anisotropy is of great importance for more accurate heat flow calculations, geodynamic studies, development and utilization of hot dry rock, and simulation of heat transfer in geological reservoirs of nuclear waste, and so on. To study the thermal conductivity anisotropy of rocks, 1158 cores from 60 boreholes in East China were tested for thermal conductivity, including thermal conductivity values parallel to (λ∥) and perpendicular to (λ⊥) structural planes of basalt, mudstones, gneisses, sandstones, carbonates, evaporites, and metamorphic rocks. The thermal conductivity anisotropy is not obvious for sand, clay, and evaporate, and the average anisotropic factors of 1.19 ± 0.22, 1.18 ± 0.17, and 1.18 ± 0.17 for tuff/breccia, granitoid and contact metamorphic rocks, respectively, indicate that these three rocks have strong anisotropy characteristics. Finally, the effect of thermal conductivity anisotropy on heat flow is studied and discussed in detail, showing that the results of thermal conductivity tests have a significant effect on the calculation of heat flow and thermal structure, and the data show that a deviation of about 10% in thermal conductivity causes a deviation of about 11% in heat flow, which may lead to a misperception of deep thermal structure studies. The regular and anisotropic characteristics of thermal conductivity of various rocks in Eastern China obtained in this paper can provide parameter support for projects such as heat flow calculations, thermal structure studies, and geothermal resource development and utilization.

scholarly journals Heat Flow in the Cretaceous of Northwestern Kansas and Implications for Regional Hydrology

1997 ◽  
pp. 1-11
Author(s):  
Andrea Forster ◽  
Daniel F. Merriam
Keyword(s):  
Thermal Conductivity ◽  
Heat Flow ◽  
Thermal Structure ◽  
Flow Density ◽  
Thermal Gradients ◽  
Conductive Heat ◽  
Flow Estimation ◽  
Average Value ◽  
Dakota Formation ◽  
Porosity And Permeability

Temperature logs are interpreted to investigate the thermal structure of the units overlying the Kansas portion of the Cretaceous Dakota aquifer. The aim of this study is to determine if additional heat input by fluids exists and thus clarify whether the overall conductive heat flow from the basement through the sequence might be overprinted by heat advection. Although interval thermal gradients are determined for different lithologic (stratigraphic) units, the shale thermal gradients are preferred for heat-flow estimation. Shale thermal conductivity as measured in Mesozoic shales in Nebraska and South Dakota is extrapolated to the area because of the similar lithology. A few thermal-conductivity values are determined in sandstone samples of the Dakota Formation and also used in heat-flow estimation. In general, the noncalcareous, marine Cretaceous shales (Pierre, Carlile, Graneros, and Kiowa) show different thermal gradients. Gradients in the Pierre (average value 58.5°C/km) and Carlile (55.5°C/km) are slightly higher than the average gradient in the Graneros Shale (45.1°C/km) and Kiowa Formation (46.5°C/km). The higher thermal gradients are limited to the extreme northwestern corner of the study area where the Pierre and Carlile are present. The heat-flow density of 69-74 mW/m2 observed there is slightly higher than the average of 60 mW/m2 typical for central and eastern Kansas. The higher heat flow observed is in the range of data reported and mapped for northeastern Colorado and the Nebraska Panhandle on the western flank of the Chadron Arch, an area with geothermal overprint by warm fluids. Regional differences in heat flow in western Kansas seemingly are caused by the different composition, porosity, and permeability of the aquifer and the nearness to recharge areas.

The effect of temperature-dependent thermal parameters in thermal models of subduction zones

2021 ◽  
Author(s):  
Iris van Zelst ◽  
Timothy J. Craig ◽  
Cedric Thieulot
Keyword(s):  
Subduction Zones ◽  
Thermal Parameters ◽  
Dislocation Creep ◽  
Seismogenic Zone ◽  
Effect Of Temperature ◽  
Temperature Dependent ◽  
Thermal Models ◽  
Intermediate Depth ◽  
Dehydration Reactions ◽  
Model Setup

<p>The thermal structure of subduction zones plays an important role in the seismicity that occurs there with e.g., the downdip limit of the seismogenic zone associated with particular isotherms (350 °C - 450 °C) and intermediate-depth seismicity linked to dehydration reactions that occur at specific temperatures and pressures. Therefore, accurate thermal models of subduction zones that include the complexities found in laboratory studies are necessary. One of the often-ignored effects in models is the temperature-dependence of the thermal parameters such as the thermal conductivity, heat capacity, and density.<span> </span></p><p>Here, we build upon the model setup presented by Van Keken et al., 2008 by including temperature-dependent thermal parameters to an otherwise clearly constrained, simple model setup of a subducting plate. We consider a fixed kinematic slab dipping at 45° and a stationary overriding plate with a dynamic mantle wedge. Such a simple setup allows us to isolate the effect of temperature-dependent thermal parameters. We add a more complex plate cooling model for the oceanic plate for consistency with the thermal parameters.<span> </span></p><p>We test the effect of temperature-dependent thermal parameters on models with different rheologies, such as an isoviscous wedge, diffusion and dislocation creep. We find that slab temperatures can change by up to 65 °C which affects the location of isotherm depths. The downdip limit of the seismogenic zone defined by e.g., the 350 °C isotherm shifts by approximately 4 km, thereby increasing the maximum possible rupture area of the seismogenic zone. Similarly, the 600 °C isotherm is shifted approximately 30 km deeper, affecting the depth at which dehydration reactions and hence intermediate-depth seismicity occurs. Our results therefore show that temperature-dependent thermal parameters in thermal models of subduction zones cannot be ignored when studying subduction-related seismicity.<span> </span></p>

scholarly journals Present-day geothermal regime of the Uliastai Depression, Erlian Basin, North China

2018 ◽  
Vol 37 (2) ◽  
pp. 770-786 ◽  
Author(s):  
Wei Xu ◽  
Shaopeng Huang ◽  
Jiong Zhang ◽  
Ruyang Yu ◽  
Yinhui Zuo ◽  
...  
Keyword(s):  
Thermal Conductivity ◽  
Heat Flow ◽  
Thermal Structure ◽  
Rheological Characteristic ◽  
Rheological Parameter ◽  
Lithospheric Thickness ◽  
Terrestrial Heat Flow ◽  
Geothermal Regime ◽  
Mantle Heat Flow ◽  
Erlian Basin

In this study, we calculated the present-day terrestrial heat flow of the Uliastai Depression in Erlian Basin by using systematical steady-state temperature data obtained from four deep boreholes and 89 thermal conductivity measurements from 22 boreholes. Then, we calculated the lithospheric thermal structure, thermal lithospheric thickness, and lithospheric thermo-rheological structure by combining crustal structure, thermal conductivity, heat production, and rheological parameter data. Research from the Depression shows that the present-day terrestrial heat flow ( qs) is 86.3 ± 2.3 mW/m2, higher than the average of 60.4 ± 12.3 mW/m2 of the continental area of China. Mantle heat flow ( qm) in the Depression ranges from 33.7 to 39.3 mW/m2, qm/ qs ranges from 40 to 44%, show that the crust plays the dominant position in the terrestrial heat flow. The thermal thickness of the lithosphere is about 74–88 km and characterized by a “strong crust–weak mantle” rheological characteristic. The total lithospheric strength is 1.5 × 1012 N/m under wet mantle conditions. Present-day geothermal regime indicates that the Uliastai Depression has a high thermal background, the activity of the deep-seated lithosphere is relatively intense. This result differs significantly from the earlier understanding that the area belongs to a cold basin. However, a hot basin should be better consistent with the evidences from lithochemistry and geophysical observations. The results also show the melts/fluids in the study area may be related to the subduction of the Paleo-Asian Ocean. The study of the geothermal regime in the Uliastai Depression provides new geothermal evidence for the volcanic activity in the eastern part of the Central Asian Orogenic Belt and has significant implications for the geodynamic characteristics.

Sign in / Sign up

Export Citation Format

Share Document