Constraining the geotherm beneath the British Isles from Bayesian inversion of Curie depth: integrated modelling of magnetic, geothermal, and seismic data

Curie depth offers a valuable constraint on the thermal structure of the lithosphere, based on its interpretation as the depth to 580 degrees C, but current methods underestimate the range of uncertainty. We formulate the estimation of Curie depth within a Bayesian framework to quantify its uncertai...

Full description

Bibliographic Details
Published in:Solid Earth
Main Authors: Mather, Ben, Fullea Urchulutegui, Javier
Format: Article in Journal/Newspaper
Language:English
Published: European Geosciences Union 2019
Subjects:
Online Access:https://hdl.handle.net/20.500.14352/92957
https://doi.org/10.5194/se-10-839-2019
Description
Summary:Curie depth offers a valuable constraint on the thermal structure of the lithosphere, based on its interpretation as the depth to 580 degrees C, but current methods underestimate the range of uncertainty. We formulate the estimation of Curie depth within a Bayesian framework to quantify its uncertainty across the British Isles. Uncertainty increases exponentially with Curie depth but this can be moderated by increasing the size of the spatial window taken from the magnetic anomaly. The choice of window size needed to resolve the magnetic thickness is often ambiguous but, based on our chosen spectral method, we determine that significant gains in precision can be obtained with window sizes 15-30 times larger than the deepest magnetic source. Our Curie depth map of the British Isles includes a combination of window sizes: smaller windows are used where the magnetic base is shallow to resolve small-scale features, and larger window sizes are used where the magnetic base is deep in order to improve precision. On average, the Curie depth increases from Laurentian crust (22.2 +/- 5.3 km) to Avalonian crust (31.2 +/- 9.2 km). The temperature distribution in the crust, and associated uncertainty, was simulated from the ensemble of Curie depth realizations assigned to a lower thermal boundary condition of a crustal model (sedimentary thickness, Moho depth, heat production, thermal conductivity), constructed from various geophysical and geochemical datasets. The uncertainty in the simulated heat flow field substantially increases from +/- 10 mW m(-2) for shallow Curie depths at similar to 15 km to +/- 80 mW m(-2) for Curie depths > 40 km. Surface heat flow observations are concordant with the simulated heat flow field except in regions that contain igneous bodies. Heat flow data within large batholiths in the British Isles exceed the simulated heat flow by similar to 25 mW m(-2) as result of their high rates of heat production (4-6 mu W m(-3)). Conversely, heat refraction around thermally resistive mafic volcanics and ...