Determining subsurface temperature & lithospheric structure from joint geophysical-petrological inversion: A case study from Ireland

Dec 1, 2023·
Emma L. Chambers
,
Raffaele Bonadio
,
Javier Fullea
,
Sergei Lebedev
,
Yihe Xu
,
Duygu Kiyan
,
Christopher J. Bean
,
Patrick A. Meere
,
Ben Mather
,
Brian M. O'Reilly
· 1 min read
Abstract
High quality maps of the geothermal gradient and temperature are essential when assessing geothermal potential for a region. However, determining geothermal potential is a challenge as direct measurements of in situ temperature are sparse and individual geophysical methods are sensitive to a range of parameters, not solely temperature. Here, we develop a novel approach to determine the geothermal gradient using a new joint geophysical-petrological inversion where seismic velocities and density in the mantle are related to temperature and bulk composition within a thermodynamic framework. Large datasets of phase velocities of seismic surface-waves are now incorporated into the inversion, and provide essential constraints on the lithospheric thickness and temperature, which shape the crustal geotherms to a significant extent. We also include all available measurements of the surface heat flow, radiogenic heat production (RHP) and thermal conductivity within the crust, to further constrain the temperature and geothermal gradient. We use Ireland as a case study and show how our new methodology can reproduce the results of previous work but also improve on them, thanks to the complementary sensitivities of the full range of data.
Type
Publication
Tectonophysics
publications

Plain Language Summary

Understanding the temperature deep underground is critical for assessing whether a region could be used for geothermal energy — a clean, renewable energy source. However, direct temperature measurements from boreholes are rare and expensive, so scientists need indirect methods to estimate subsurface temperatures over wide areas.

This study develops a new approach that combines multiple types of geophysical data — including seismic wave speeds, gravity measurements, surface heat flow, and rock property measurements — within a single computational framework. By jointly inverting all of these datasets together, the method produces more reliable maps of underground temperature and the structure of the Earth’s outer layer (the lithosphere) than any single dataset could provide alone.

Applied to Ireland as a test case, the results show that the thickness of the lithosphere and the crust are the primary controls on the geothermal gradient, with thinner lithosphere areas showing higher temperatures closer to the surface. In some locations, rocks rich in radioactive elements generate extra heat that boosts the geothermal gradient locally. This methodology can be applied to any region with limited direct temperature data, helping to identify promising areas for geothermal energy development.