scholarly journals Adjoint inversion of the thermal structure of Southeastern Australia

2019 ◽  
Vol 219 (3) ◽  
pp. 1648-1659 ◽  
Author(s):  
B Mather ◽  
L Moresi ◽  
P Rayner
Keyword(s):  
Thermal Properties ◽  
Heat Flow ◽  
Heat Production ◽  
Thermal Structure ◽  
Surface Heat ◽  
Flow Data ◽  
Data Set ◽  
Surface Heat Flow ◽  
Southeastern Australia ◽  
Curie Depth

SUMMARY The variation of temperature in the crust is difficult to quantify due to the sparsity of surface heat flow observations and lack of measurements on the thermal properties of rocks at depth. We examine the degree to which the thermal structure of the crust can be constrained from the Curie depth and surface heat flow data in Southeastern Australia. We cast the inverse problem of heat conduction within a Bayesian framework and derive its adjoint so that we can efficiently find the optimal model that best reproduces the data and prior information on the thermal properties of the crust. Efficiency gains obtained from the adjoint method facilitate a detailed exploration of thermal structure in SE Australia, where we predict high temperatures within Precambrian rocks of 650 °C due to relatively high rates of heat production (0.9–1.4 μW m−3). In contrast, temperatures within dominantly Phanerozoic crust reach only 520 °C at the Moho due to the low rates of heat production in Cambrian mafic volcanics. A combination of the Curie depth and heat flow data is required to constrain the uncertainty of lower crustal temperatures to ±73 °C. We also show that parts of the crust are unconstrained if either data set is omitted from the inversion.

Towards a global lithospheric thermal model

2020 ◽  
Author(s):  
Sheona Masterton ◽  
Samuel Cheyney ◽  
Chris Green ◽  
Peter Webb
Keyword(s):  
Heat Flow ◽  
Curie Temperature ◽  
Thermal Structure ◽  
Magnetic Layer ◽  
Surface Heat ◽  
Local Context ◽  
Magnetic Data ◽  
Flow Data ◽  
Surface Heat Flow ◽  
Temperature Depth

<p>Temperature and heat flow are key parameters for understanding the potential for source rock maturation in sedimentary basins. Knowledge of the thermal structure of the lithosphere in both a regional and local context can provide important constraints for modelling basin evolution through time.</p><p>In recent years, global coverage of heat flow data constraints have enhanced scientific understanding of the thermal state of the lithosphere. However, sample bias and variability in sampling methods continues to be a major obstacle to heat flow-derived isotherm prediction, particularly in frontier areas where data are often sparse or poorly constrained. Consideration and integration of alternative approaches to predict temperature at depth may allow interpolation of surface heat flow in such data poor areas.   </p><p>We have attempted to integrate three independent approaches to modelling temperature with depth. The first approach is based on heat flow observations, in which a 1D steady-state model of the lithosphere is constructed from quality-assessed surface heat flow data, crustal thickness estimates and associated lithospheric thermal properties. The second approach is based on terrestrial (airborne, ground and shipborne) magnetic data, in which the maximum depth of magnetisation within the lithosphere is estimated using a de-fractal method and used as a proxy for Curie temperature depth. The third approach is based on satellite magnetic data and estimates the thickness of the magnetic layer within the lithosphere based on the varying amplitudes of satellite magnetic data, accounting for global variations in crustal magnetisation. Curie temperature depth results from each of these approaches have been integrated into a single global grid, then used to calculate temperature-depth variations through the crust.</p><p>We have evaluated our isotherm predictions by comparing them with temperature-depth control points and undertook qualitative and quantitative analyses of discrepancies that exist between different modelling approaches; this has provided insights into the origin of such discrepancies that can be integrated into our models to generate a better controlled global temperature-depth result.  </p><p>We present details of our methodology and the results of our integrated studies. We demonstrate areas where the independent results are in good agreement, providing vital information for high-level basin screening. We also highlight areas of disagreement and suggest possible causes for these discrepancies and potential resolutions.</p>

Temperature structure of the Andean subduction zone as derived from the 3D geometry of crustal and upper mantle discontinuities

2020 ◽  
Author(s):  
Andres Tassara ◽  
Joaquín Julve ◽  
Iñigo Echeverría ◽  
Ingo Stotz
Keyword(s):  
Heat Flow ◽  
Heat Production ◽  
Thermal Structure ◽  
Surface Heat ◽  
Thermal Parameters ◽  
Continental Lithosphere ◽  
Flow Measurements ◽  
Study Region ◽  
Surface Heat Flow ◽  
Input Model

<p>The distribution of temperature inside active continental margins plays a fundamental role on regulating first order geodynamic processes as the isostatic balance, rheologic behavior of crust and mantle, magmagenesis, volcanism and seismogenesis. In spite of these major implications, well-constrained 3D thermal models are known for few regions of the world (Europe, Western USA, China) where large geophysical databases have been integrated into compositional and structural models of crust and lithospheric mantle from which a thermal model is derived. Here we present a three-dimensional representation of the distribution of temperature underneath the Andean active margin of South America (10°-45°S) that is based on a geophysically-constrained model for the geometry of the subducted slab, continental lithosphere-asthenosphere boundary (LAB), Moho discontinuity and an intracrustal discontinuity (ICD). This input model was constructed by forward modelling the satellite gravity anomaly under the constraint of most of the seismic information published for this region. We use analytical expressions of 1D conductive continental geotherms with adequate boundary conditions that consider the compositional stratification of crust and mantle included in the input model, and the advective thermal effect of slab subduction. The 1D geotherms are assembled into a 3D volume defining the thermal structure of the study region. We test the influence of several thermal parameters and structural configurations of the Andean lithosphere by comparing the resulting surface heat flow distribution of these different models against a database containing heat flow measurements that we compile from the literature. Our results show that the thermal structure and derived surface heat flow is dominantly controlled by the geometry of the thermal boundary layer at the base of the lithosphere, i.e. the slab upper surface below the forearc and LAB inland. Variations on the modeled configuration of the continental lithosphere (i.e. the way on which the geometry of the continental Moho and ICD are considered into the definition of a space-variable thermal conductivity or the length scale for radiogenic heat production) have an effect on surface heat flow that is lower than the average uncertainty of the measurements and therefore can be considered as second-order. The simplicity of our analytical approach allows us to compute hundreds of different models in order to test the sensitivity of results to changes on thermal parameters (conductivity, heat production, mantle potential temperature, etc), which provides a tool for discussing their possible range of values in the context of a subduction margin. We will also show how variations of these models impact on the Moho temperature and therefore in the expected mechanical behavior of crust and mantle in this geotectonic context</p>

Heat production and thermal conductivity of rocks from the Pikwitonei–Sachigo continental cross section, central Manitoba: implications for the thermal structure of Archean crust

1987 ◽  
Vol 24 (8) ◽  
pp. 1583-1594 ◽  
Author(s):  
David M. Fountain ◽  
Matthew H. Salisbury ◽  
Kevin P. Furlong
Keyword(s):  
Thermal Conductivity ◽  
Heat Flow ◽  
Cross Section ◽  
Heat Production ◽  
Thermal Structure ◽  
Granulite Facies ◽  
Surface Heat ◽  
Gamma Ray Spectrometry ◽  
Surface Heat Flow ◽  
Greenstone Belts

The Pikwitonei and Sachigo subprovinces of central Manitoba provide a cross-sectional view of the Superior Province crust. In cross section, the upper to mid-level crust is composed of synformal greenstone belts surrounded by tonalitic gneisses, both of which are intruded by granitoid plutons. This crustal structure persists downward into the granulite facies, where keels of the greenstone belts can be found. To constrain thermal models of the crust, we measured heat production and thermal conductivity in 60 rocks from this terrain using standard gamma-ray spectrometry and divided bar techniques. Large vertical and lateral heterogeneities in heat production in the upper crust are evident; heat production is high in granites and metasedimentary rocks, intermediate in tonalite gneisses, and low in the portions of greenstone belts dominated by mafic meta-igneous rocks. In the deeper granulite facies rocks, heat production decreases by a factor of two in the tonalitic gneisses and remains low in the high-grade mafic rocks. When applied to the Pikwitonei–Sachigo crust cross section, the laboratory data here do not support step function or exponential models of the variation of heat production with depth. However, estimates of surface heat flow and surface heat production for various sites in the crustal model yield the well-known linear relationship between surface heat production and surface heat flow observed for heat-flow provinces for both one- and two-dimensional models. This demonstrates that determinations of heat production with depth based on inversion of the linear heat-production–heat-flow relationship are nonunique.

A lightweight thermal modelling tool for physics-based continental heat flow interpolation

2020 ◽  
Author(s):  
Alberto Pastorutti ◽  
Carla Braitenberg
Keyword(s):  
Temperature Field ◽  
Heat Flow ◽  
Central Europe ◽  
Heat Production ◽  
Surface Heat ◽  
Gravity Models ◽  
Forward Modelling ◽  
Flow Data ◽  
Satellite Gravity ◽  
Surface Heat Flow

<p><em>Both energy applications, such as assessing one of the controlling factors of conductive geothermal plays, and geodynamics modelling, are influenced by the large uncertainties arising from uneven sampling of the direct observable of the Earth's thermal state, surface heat flow. Heterogeneity in structure and composition of the continental lithosphere complicate the temperature field even in stable provinces in thermal equilibrium. The measurements deviate from what simple relationships with geological and geophysical data predict, requiring more sophisticated schemes such as those based on multivariate inversion (e.g. Mather et al. 2018) and geostatistics (e.g. the similarity method employed by Lucazeau, 2019).</em></p><p><em>Recently, we aimed at assessing the performance of satellite-gravity-constrained modelling of surface heat flow [1], with the aim of employing the unparalleled spatial uniformity of global gravity models in the fill-in of sparsely sampled surface heat flow data. The model we obtained, in a test area in Central Europe, provided additional information on the lithospheric structure and revealed a satisfactory coherence with the geological features in the area and their controlling effect on the conductive heat transport. That test was based on a fit of radioactive heat production to available heat flow data, based on a misfit linearization and substitution strategy, which we have shown to be independently consistent with available heat production relationships (e.g. Hasterok and Webb, 2017). Furthermore, model validation techniques provide additional metrics on the predictability in areas devoid of heat flow measurements.</em></p><p><em>T</em><em>o reach those objectives, we developed a finite-difference based solver for the heat equation in conductive, stable lithosphere, relying on the assumption of steady state, 3-D heat conduction from the thermal base of the lithosphere to surface. It allows for non-homogeneous heat production and thermal conductivity, and non-flat upper and bottom boundaries. Concurrent joint forward modelling of the gravity field is also possible.<br>Through compromise between complexity and approximation, it was designed favouring easy and fast forward modelling, such as in assessing parameter sensitivity and performing grid searches or parameter fitting. Geological models and parameters can be defined using an user-friendly plain text layer-wise definition, which is then turned into a volume, on a rectangular mesh.<br>Computational requirements are lean: a 75 × 75 × 104 node model such as the one employed in [1] can be forward-modelled on an ordinary workstation in 135 seconds. A direct solver is employed to solve the FD system of linear equations: the Matlab built-in Cholesky decomposition for sparse arrays (Davis, 2006).</em></p><p><em>Albeit initially developed as an ad-hoc tool for a proof of concept, its ease of use and versatility suggest its potential in other applications. We therefore present the solver and the accompanying tool set, both openly available, along with a set of promising examples.<br><br>[1] Pastorutti, A., Braitenberg, C. (2019) "A geothermal application for GOCE satellite gravity data: modelling the crustal heat production and lithospheric temperature field in Central Europe." Geophysical Journal International, doi:10.1093/gji/ggz344</em></p>

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

Solid Earth ◽  
2019 ◽  
Vol 10 (3) ◽  
pp. 839-850 ◽  
Author(s):  
Ben Mather ◽  
Javier Fullea
Keyword(s):  
Heat Flow ◽  
Flow Field ◽  
Mantle Plume ◽  
Heat Production ◽  
Surface Heat ◽  
Integrated Modelling ◽  
Moho Depth ◽  
British Isles ◽  
Surface Heat Flow ◽  
Curie Depth

Abstract. Curie depth offers a valuable constraint on the thermal structure of the lithosphere, based on its interpretation as the depth to 580 ∘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 ∼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 ∼25 mW m−2 as a result of their high rates of heat production (4–6 µW m−3). 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 (>90 mW m−2) and shallow Curie depth (∼15 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 ±3–8 km in this region demonstrates that Curie depths are more reliable in hotter regions of the crust where the magnetic base is shallow.

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

2019 ◽  
Author(s):  
Ben Mather ◽  
Javier Fullea
Keyword(s):  
Heat Flow ◽  
Flow Field ◽  
Mantle Plume ◽  
Heat Production ◽  
Surface Heat ◽  
Integrated Modelling ◽  
Moho Depth ◽  
British Isles ◽  
Surface Heat Flow ◽  
Curie Depth

Abstract. Curie depth offers a valuable constraint on the thermal structure of the lithosphere, based on its interpretation as the depth to 580 °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 windows 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 realisations 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 data sets. The uncertainty of the simulated heat flow field substantially increases from ± 10 mW m−2 for shallow Curie depths ~ 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 ∼ 25 mW m−2 as a result of their high rates of heat production (4–6 μW m−3). 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 (> 90 mW m−2) and shallow Curie depth (∼ 15 km) occur within the same region, which supports the hypothesis that lithospheric thinning occurred due to the influence of a mantle plume. That the uncertainty is only ± 3–8 km in this region, demonstrates that Curie depths are more reliable in hotter regions of the crust where the magnetic base is shallow.

scholarly journals Surface heat flow and lithosphere thermal structure of the Rhenohercynian Zone in the greater Luxembourg region

Geothermics ◽  
2015 ◽  
Vol 56 ◽  
pp. 93-109 ◽  
Author(s):  
Tom Schintgen ◽  
Andrea Förster ◽  
Hans-Jürgen Förster ◽  
Ben Norden
Keyword(s):  
Heat Flow ◽  
Thermal Structure ◽  
Surface Heat ◽  
Surface Heat Flow ◽  
Rhenohercynian Zone

scholarly journals Crustal heat production and estimate of terrestrial heat flow in central East Antarctica, with implications for thermal input to the East Antarctic ice sheet

2018 ◽  
Vol 12 (2) ◽  
pp. 491-504 ◽  
Author(s):  
John W. Goodge
Keyword(s):  
Heat Flow ◽  
Heat Production ◽  
Ice Sheet ◽  
Surface Heat ◽  
East Antarctica ◽  
Antarctic Ice Sheet ◽  
Surface Heat Flow ◽  
Terrestrial Heat Flow ◽  
Production Values ◽  
East Antarctic Ice Sheet

Abstract. Terrestrial heat flow is a critical first-order factor governing the thermal condition and, therefore, mechanical stability of Antarctic ice sheets, yet heat flow across Antarctica is poorly known. Previous estimates of terrestrial heat flow in East Antarctica come from inversion of seismic and magnetic geophysical data, by modeling temperature profiles in ice boreholes, and by calculation from heat production values reported for exposed bedrock. Although accurate estimates of surface heat flow are important as an input parameter for ice-sheet growth and stability models, there are no direct measurements of terrestrial heat flow in East Antarctica coupled to either subglacial sediment or bedrock. As has been done with bedrock exposed along coastal margins and in rare inland outcrops, valuable estimates of heat flow in central East Antarctica can be extrapolated from heat production determined by the geochemical composition of glacial rock clasts eroded from the continental interior. In this study, U, Th, and K concentrations in a suite of Proterozoic (1.2–2.0 Ga) granitoids sourced within the Byrd and Nimrod glacial drainages of central East Antarctica indicate average upper crustal heat production (Ho) of about 2.6  ±  1.9 µW m−3. Assuming typical mantle and lower crustal heat flux for stable continental shields, and a length scale for the distribution of heat production in the upper crust, the heat production values determined for individual samples yield estimates of surface heat flow (qo) ranging from 33 to 84 mW m−2 and an average of 48.0  ±  13.6 mW m−2. Estimates of heat production obtained for this suite of glacially sourced granitoids therefore indicate that the interior of the East Antarctic ice sheet is underlain in part by Proterozoic continental lithosphere with an average surface heat flow, providing constraints on both geodynamic history and ice-sheet stability. The ages and geothermal characteristics of the granites indicate that crust in central East Antarctica resembles that in the Proterozoic Arunta and Tennant Creek inliers of Australia but is dissimilar to other areas like the Central Australian Heat Flow Province that are characterized by anomalously high heat flow. Age variation within the sample suite indicates that central East Antarctic lithosphere is heterogeneous, yet the average heat production and heat flow of four age subgroups cluster around the group mean, indicating minor variation in the thermal contribution to the overlying ice sheet from upper crustal heat production. Despite these minor differences, ice-sheet models may favor a geologically realistic input of crustal heat flow represented by the distribution of ages and geothermal characteristics found in these glacial clasts.

scholarly journals Crustal heat production and estimate of terrestrial heat flow in central East Antarctica, with implications for thermal input to the East Antarctic ice sheet

2017 ◽  
Author(s):  
John W. Goodge
Keyword(s):  
Heat Flow ◽  
Heat Production ◽  
Ice Sheet ◽  
Surface Heat ◽  
East Antarctica ◽  
Antarctic Ice Sheet ◽  
Surface Heat Flow ◽  
Terrestrial Heat Flow ◽  
Production Values ◽  
East Antarctic Ice Sheet

Abstract. Terrestrial heat flow is a critical first-order factor governing the thermal condition and, therefore, mechanical stability of Antarctic ice sheets, yet heat flow across Antarctica is poorly known. Previous estimates of terrestrial heat flow come from inversion of seismic and magnetic geophysical data, by modeling temperature profiles in ice boreholes, and by calculation from heat production values reported for exposed bedrock. Although accurate estimates of surface heat flow are important as an input parameter for ice-sheet growth and stability models, there are no direct measurements of terrestrial heat flow in East Antarctica coupled to either subglacial sediment or bedrock. Bedrock outcrop is limited to coastal margins and rare inland exposures, yet valuable estimates of heat flow in central East Antarctica can be extrapolated from heat production determined by the geochemical composition of glacial rock clasts eroded from the continental interior. In this study, U, Th and K concentrations in a suite of Proterozoic (1.2–2.0 Ga) granitoids sourced within the Byrd and Nimrod glacial drainages of central East Antarctica indicate average upper crustal heat production (Ho) of about 2.6 ± 1.9 μW m-3. Assuming typical mantle and lower crustal heat flux for stable continental shields, and a length scale for the distribution of heat production in the upper crust, the heat production values determined for individual samples yield estimates of surface heat flow (qo) ranging from 33–84 mW m-2 and an average of 48.0 ± 13.6 mW m-2. Estimates of heat production obtained for this suite of glacially-sourced granitoids therefore indicate that the interior of the East Antarctic ice sheet is underlain in part by Proterozoic continental lithosphere with average surface heat flow, providing constraints on both geodynamic history and ice-sheet stability. The ages and geothermal characteristics of the granites indicate that crust in central East Antarctica resembles that in the Proterozoic Arunta and Tenant Creek inliers of Australia, but is dissimilar to other areas characterized by anomalously high heat flow in the Central Australian Heat Flow Province. Age variation within the sample suite indicates that central East Antarctic lithosphere is heterogeneous, yet the average heat production and heat flow of four age subgroups cluster around the group mean, indicating minor variation in thermal contribution to the overlying ice sheet from upper crustal heat production. Despite their minor differences, ice-sheet models may favor a geologically realistic model of crustal heat flow represented by such a distribution of ages and geothermal characteristics.

Sign in / Sign up

Export Citation Format

Share Document