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

<p>Curie depth offers a valuable constraint on the thermal structure of the lithosphere, based on its interpretation as the depth to 580&thinsp;<span class="inline-formula"><sup>∘</sup></span>C, but current methods underestimate the range of uncertainty. W...

Full description

Bibliographic Details
Main Authors: B. Mather, J. Fullea
Format: Article
Language:English
Published: Copernicus Publications 2019-06-01
Series:Solid Earth
Online Access:https://www.solid-earth.net/10/839/2019/se-10-839-2019.pdf
_version_ 1818926043046608896
author B. Mather
J. Fullea
author_facet B. Mather
J. Fullea
author_sort B. Mather
collection DOAJ
description <p>Curie depth offers a valuable constraint on the thermal structure of the lithosphere, based on its interpretation as the depth to 580&thinsp;<span class="inline-formula"><sup>∘</sup></span>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 (<span class="inline-formula">22.2±5.3</span>&thinsp;km) to Avalonian crust (<span class="inline-formula">31.2±9.2</span>&thinsp;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 <span class="inline-formula">±10</span>&thinsp;mW&thinsp;m<span class="inline-formula"><sup>−2</sup></span> for shallow Curie depths at <span class="inline-formula">∼15</span>&thinsp;km to <span class="inline-formula">±80</span>&thinsp;mW&thinsp;m<span class="inline-formula"><sup>−2</sup></span> for Curie depths <span class="inline-formula">&gt;40</span>&thinsp;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 <span class="inline-formula">∼25</span>&thinsp;mW&thinsp;m<span class="inline-formula"><sup>−2</sup></span> as a result of their high rates of heat production (4–6&thinsp;<span class="inline-formula">µ</span>W&thinsp;m<span class="inline-formula"><sup>−3</sup></span>). Conversely, heat refraction around thermally resistive mafic volcanics and thick sedimentary layers induce a negative heat flow misfit of a similar magnitude. A northward thinning of the lithosphere is supported by shallower Curie depths on the northern side of the Iapetus Suture, which separates Laurentian and Avalonian terranes. Cenozoic volcanism in Northern Britain and Ireland has previously been attributed to a lateral branch of the proto-Icelandic mantle plume. Our results show that high surface heat flow (<span class="inline-formula">&gt;90</span>&thinsp;mW&thinsp;m<span class="inline-formula"><sup>−2</sup></span>) and shallow Curie depth (<span class="inline-formula">∼15</span>&thinsp;km) occur within the same region, which supports the hypothesis that lithospheric thinning occurred due to the influence of a mantle plume. The fact that the uncertainty is only <span class="inline-formula">±3</span>–8&thinsp;km in this region demonstrates that Curie depths are more reliable in hotter regions of the crust where the magnetic base is shallow.</p>
first_indexed 2024-12-20T02:50:50Z
format Article
id doaj.art-af05e45ba89e435fa3de4e00af7804e8
institution Directory Open Access Journal
issn 1869-9510
1869-9529
language English
last_indexed 2024-12-20T02:50:50Z
publishDate 2019-06-01
publisher Copernicus Publications
record_format Article
series Solid Earth
spelling doaj.art-af05e45ba89e435fa3de4e00af7804e82022-12-21T19:56:03ZengCopernicus PublicationsSolid Earth1869-95101869-95292019-06-011083985010.5194/se-10-839-2019Constraining the geotherm beneath the British Isles from Bayesian inversion of Curie depth: integrated modelling of magnetic, geothermal, and seismic dataB. MatherJ. Fullea<p>Curie depth offers a valuable constraint on the thermal structure of the lithosphere, based on its interpretation as the depth to 580&thinsp;<span class="inline-formula"><sup>∘</sup></span>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 (<span class="inline-formula">22.2±5.3</span>&thinsp;km) to Avalonian crust (<span class="inline-formula">31.2±9.2</span>&thinsp;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 <span class="inline-formula">±10</span>&thinsp;mW&thinsp;m<span class="inline-formula"><sup>−2</sup></span> for shallow Curie depths at <span class="inline-formula">∼15</span>&thinsp;km to <span class="inline-formula">±80</span>&thinsp;mW&thinsp;m<span class="inline-formula"><sup>−2</sup></span> for Curie depths <span class="inline-formula">&gt;40</span>&thinsp;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 <span class="inline-formula">∼25</span>&thinsp;mW&thinsp;m<span class="inline-formula"><sup>−2</sup></span> as a result of their high rates of heat production (4–6&thinsp;<span class="inline-formula">µ</span>W&thinsp;m<span class="inline-formula"><sup>−3</sup></span>). Conversely, heat refraction around thermally resistive mafic volcanics and thick sedimentary layers induce a negative heat flow misfit of a similar magnitude. A northward thinning of the lithosphere is supported by shallower Curie depths on the northern side of the Iapetus Suture, which separates Laurentian and Avalonian terranes. Cenozoic volcanism in Northern Britain and Ireland has previously been attributed to a lateral branch of the proto-Icelandic mantle plume. Our results show that high surface heat flow (<span class="inline-formula">&gt;90</span>&thinsp;mW&thinsp;m<span class="inline-formula"><sup>−2</sup></span>) and shallow Curie depth (<span class="inline-formula">∼15</span>&thinsp;km) occur within the same region, which supports the hypothesis that lithospheric thinning occurred due to the influence of a mantle plume. The fact that the uncertainty is only <span class="inline-formula">±3</span>–8&thinsp;km in this region demonstrates that Curie depths are more reliable in hotter regions of the crust where the magnetic base is shallow.</p>https://www.solid-earth.net/10/839/2019/se-10-839-2019.pdf
spellingShingle B. Mather
J. Fullea
Constraining the geotherm beneath the British Isles from Bayesian inversion of Curie depth: integrated modelling of magnetic, geothermal, and seismic data
Solid Earth
title Constraining the geotherm beneath the British Isles from Bayesian inversion of Curie depth: integrated modelling of magnetic, geothermal, and seismic data
title_full Constraining the geotherm beneath the British Isles from Bayesian inversion of Curie depth: integrated modelling of magnetic, geothermal, and seismic data
title_fullStr Constraining the geotherm beneath the British Isles from Bayesian inversion of Curie depth: integrated modelling of magnetic, geothermal, and seismic data
title_full_unstemmed Constraining the geotherm beneath the British Isles from Bayesian inversion of Curie depth: integrated modelling of magnetic, geothermal, and seismic data
title_short Constraining the geotherm beneath the British Isles from Bayesian inversion of Curie depth: integrated modelling of magnetic, geothermal, and seismic data
title_sort constraining the geotherm beneath the british isles from bayesian inversion of curie depth integrated modelling of magnetic geothermal and seismic data
url https://www.solid-earth.net/10/839/2019/se-10-839-2019.pdf
work_keys_str_mv AT bmather constrainingthegeothermbeneaththebritishislesfrombayesianinversionofcuriedepthintegratedmodellingofmagneticgeothermalandseismicdata
AT jfullea constrainingthegeothermbeneaththebritishislesfrombayesianinversionofcuriedepthintegratedmodellingofmagneticgeothermalandseismicdata