Loading…

Refining the thermal structure of the European lithosphere by inversion of subsurface temperature data

We present an updated thermal model of the European lithosphere based on a new stochastic modeling work flow. We developed this work flow to estimate subsurface temperatures from site- to regional-scale, up to the depth of the lithosphere-asthenosphere boundary (LAB). Our model is composed of four l...

Full description

Saved in:
Bibliographic Details
Published in:Global and planetary change 2018-12, Vol.171, p.18-47
Main Authors: Limberger, Jon, van Wees, Jan-Diederik, Tesauro, Magdala, Smit, Jeroen, Bonté, Damien, Békési, Eszter, Pluymaekers, Maarten, Struijk, Maartje, Vrijlandt, Mark, Beekman, Fred, Cloetingh, Sierd
Format: Article
Language:English
Subjects:
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:We present an updated thermal model of the European lithosphere based on a new stochastic modeling work flow. We developed this work flow to estimate subsurface temperatures from site- to regional-scale, up to the depth of the lithosphere-asthenosphere boundary (LAB). Our model is composed of four layers, consisting of sediments, upper crust, lower crust and lithospheric mantle. We asigned thermal properties, including radiogenic heat production and temperature- and pressure-dependent bulk thermal conductivity on the base of broad-scale lithological variation within the European crust. We corrected thermal properties with a 1D steady-state temperature approximation, assuming only vertical heat flow. Using these corrected thermal properties, we calculated the 3D thermal field with a conjugate-gradient method, assuming fixed temperatures at the surface and at the base of the lithosphere. To obtain more robust results for our thermal model, we applied data assimilation, aiming at consistency between temperature and heat flow observations and tectonic model predictions. We used an Ensemble Smoother Multiple with Data Assimilation (ES-MDA) method to update prior estimates of thermal properties and the thermal field with temperature data. We calibrated our European thermal model with available regional thermal models due to the present lack of a unified dataset with public borehole temperature measurements. A large dichotomy is observed in the model along the Trans-European Suture Zone (TESZ). Northeast of the TESZ, geothermal gradients up to 10 km depth are mainly below 20 °C km−1. Southwest of the TESZ, gradients range from 20 °C km−1 near the Adriatic coast to more than 50 °C km−1 in volcanically and tectonically active regions. We show that the large scale thermal structure is locally and regionally perturbed by non-conductive heat transfer and affected by transient effects. •An improved 3D thermal model of the European lithosphere based on a new stochastic modeling work flow.•Physics-based reference model for thermal properties and boundary conditions at EU scale for thermal characterization.•Model calibration with 35533 temperature points derived from regional thermal models.•Updated range of thermal properties using the ensemble smoother with multiple data assimilation method (ES-MDA).
ISSN:0921-8181
1872-6364
DOI:10.1016/j.gloplacha.2018.07.009