High Regional Heat Flow Research Articles

Geothermal parameter assessment in the Southwestern Sokoto Basin, Nigeria, using spectral analysis (centroid method)

This study utilized spectral analysis (centroid method) to assess geothermal parameters in the southwestern part of the Sokoto Basin, Nigeria. The high-resolution airborne data comprised forty-nine (49) overlapping blocks, and each block was divided into 55x55 km to evaluate essential parameters such as depth to the top boundary (Zt ), centroiddepth (Zo), and magnetic source bottom (Zb = 2Zo-Zt ). The analysis revealed variable Curie points depths (CPD), ranging from 3.89 km to 26.56 km. The lowest CPD is primarily associated with basement rocks within anomalies A, B, C, D, E, F, G, H, and I, with an average CPD of 9.16 km. Furthermore, the thermal gradients ranged from 21.84 ◦C/km to 149.10 ◦C/km, with an average thermal gradient of 73.30 ◦C/km. The heat flow exhibited variations between 54.81 mW/m² and 374.24 mW/m², with average heat flow of 180.4 mW/m², indicating significant geothermal potential zones. The high thermal gradients and heat flow regions were identified, around anomalies A, B, and C. Additionally, temperature gradients identified at shallow depths ranged from 110 ◦C/km to 150 ◦C/km. The results reveal the presence of high-temperature points and anomalous geothermal potentials, particularly within anomalies A, B, and C, thus requiring further for sustainable geothermal energy generation in the study area.

Two Sets of High-Conductivity Systems with Different Scales Reveal the Seismogenic Mechanism of Earthquakes in the Songyuan Area, Northeastern China

To reveal the deep seismogenic environment and mechanism of earthquakes in Songyuan City, Northeastern China, 59 broadband magnetotelluric sites in the Songyuan area were arranged in this study at a spacing of 5 km. In addition, two intersecting magnetotelluric profiles, with a total of 23 measuring sites and a spacing of 2 km, were established near the Ningjiang earthquake swarm. Using a nonlinear conjugate gradient (NLCG) algorithm, resistivity structures in the lithosphere were obtained at different scales using three-dimensional (3D) inversion. The research results show that: a deep high-conductivity system (<10 Ω·m) was identified at 25–85 km depth in the lithosphere under Songyuan, corresponding closely to a region of high heat flow. It is inferred to be the molten material of mantle upwelling. In addition, a shallow high-conductivity system (<10 Ω·m) was identified beneath the Ningjiang earthquake swarm, which is interpreted to correspond to the Fuyu North fault. It is the main seismo-controlling structure of the Ningjiang earthquake swarm. The deep seismogenic environment and seismogenic mechanism of the Ningjiang earthquake swarm can be described as a deep upwelling of molten mantle material, which provides the power source. The deep magma intruded into the lower crust and accumulated, then intruded along faults and fissures, resulting in the activation of the North Fuyu fault and triggering the Ningjiang earthquake. It is attributed to the activation of shallow faults caused by the upwelling of molten mantle material.

Curie point depth from airborne magnetic data analysis for geothermal exploration in southern Thailand

Southern Thailand is a high heat flow region, which implies to high-potential area of geothermal resources. Direct heat flow measurements are only available from petroleum wells in Gulf of Thailand. Therefore, indirect method of Curie Point Depth (CPD) calculation has been used to evaluate the heat flow on land of southern Thailand. CPD is the theoretical depth at which subsurface temperature reaches the value of the Curie point of magnetite, which is 578 °C. CPD can be calculated by deriving the depth of the bottom of the magnetic body to the surface using spectral analysis from airborne magnetic data, and the calculated values can produce the CPD map of southern Thailand. The results show that CPD ranges from 7.88 to 19.75 km with an average of 13.74 km, and heat flow ranges from 61.54 to 154.25 mW/m2 with an average of 90.36 mW/m2. This study suggests that the high heat flow in southern Thailand, which is higher than the regular continental crust heat flow of 65 mW/m2, is possibly caused by a high heat anomaly from mantle upwelling and expressed through weak zones of faults and sutures in extensional tectonic regimes. The high heat from mantle has elevated by five zones: (1) Ranong faults zone, (2) Khlong Marui faults zone, (3) Coastline of Surat Thani and Nakhon Si Thammarat, (4) Trang and Satun zone, and (5) Bentong-Roab suture. These five zones are considered to be high-potential areas for future geothermal resources exploration.

Appraisal of geothermal potentials of some parts of the Abakaliki Anticlinorium and adjoining areas (Southeast Nigeria) using magnetic data

The Abakaliki Anticlinorium and its adjoining areas were appraised with the object of delineating high geothermal potential zones. Spectral depth analysis involving an improved centroid technique was used to analyze high-quality magnetic data. The obtained geothermal parameters were gridded to map various geothermal features within the investigated area. The obtained results varied from 4.99–9.35 km, 2.31–6.15 km, 6.11–16.28 km, 35.63°C–94.93°C/km, and 89.07–237.32 mW/m2 for centroid depth, top depth, Curie point depth, geothermal gradient, and heat flow values, respectively. The delineated semioval structure in the central zone of the investigated region characterized by a shallow Curie point depth (< 8.5 km) correlates with the location of the high-heat flow (>191.0 mW/m2) and geothermal gradient (>74.0°C/km) region. The high geothermal potential of the region is triggered by the massive post-rift tectonic event of the Santonian period related to the Abakaliki Anticlinorium. Further geophysical exploration programs should be carried out before exploitation activities at anomalous geothermal regions.

Modeling the maturation history of the stacked petroleum systems of the Williston Basin, USA

A three-dimensional petroleum systems model was built to support U.S. Geological Survey assessments of undiscovered oil and gas resources in the Williston Basin of North Dakota, Montana, and South Dakota. Numerous Paleozoic source rocks have been proven or postulated in the basin, of which five were the focus of maturation and migration modeling: the Ordovician Icebox Formation, the kukersite beds of the Ordovician Red River Formation, the shales of the Devonian–Mississippian Bakken Formation, the Mississippian Madison Group, and the Pennsylvanian Tyler Formation. Calibration of the three-dimensional model to present-day temperature data indicates the existence of a north-south trend of high heat flow in western North Dakota, along with a region of high heat flow in eastern Montana. These high heat flow trends strongly control the maturity of all studied source intervals. A Bakken-specific hydrocarbon generation kinetic model was developed to match the calibrated time-temperature history of the basin to spatial trends in hydrogen index from programmed pyrolysis data. Generation of hydrocarbons occurred in the Cretaceous through Paleogene due to increased burial. Subsequent uplift and erosion in the Neogene cooled the basin, ending hydrocarbon generation for all source rocks. The cumulative volume of hydrocarbons generated by each of the source rocks was calculated and used to compare their relative robustness. The shales of the Bakken Formation are estimated to have generated approximately 460 billion barrels of oil equivalent (BBOE), while the Red River Formation generated approximately 130 BBOE, the Tyler Formation 94 BBOE, the Madison Group 44 BBOE, and the Icebox Formation 28 BBOE. Gross migration trends were analyzed with respect to historical oil and gas production in the basin and generally indicate segregation of petroleum systems throughout the stratigraphic column. However, most modeled scenarios indicated significant loss of Bakken oil to the Madison Group, suggesting that mixing of Madison and Bakken oils may be more prevalent than has recently been recognized in the U.S. portion of the Williston Basin and is particularly likely in fractured regions of the basin.

Imaging the magmatic plumbing of the Clear Lake Volcanic Field using 3-D gravity inversions

The Quaternary Clear Lake Volcanic Field (CLVF) in the Northern California Coast Range is the youngest of a string of northward-younging volcanic centers in the state. The CLVF is located within the broad San Andreas Transform Fault System and has been active intermittently for ∼2 million years. Heat beneath the CLVF supports The Geysers, one of the largest producing geothermal fields in the world.Previous geophysical studies proposed the existence of a magma reservoir beneath Mount Hannah, which is northeast of The Geysers, near the geographic center of the CLVF. The lateral extent, depth, and presence of melt within this reservoir are poorly constrained, as is the relationship between this body and the broader magmatic plumbing of the CLVF. To gain a clearer and more comprehensive picture of the CLVF magma source region, a gravity dataset was compiled and the first 3-D gravity inversions of the CLVF were completed.Field and synthetic model inversions from the current study both indicate that the gravity low roughly centered on Mount Hannah is not accurately explained by a 5–7 km thick lens of Mesozoic Great Valley Sequence (ρ=2.58g/cm3) as proposed by Stanley et al. (1998). The observed gravity low is more accurately described by one or more silicic, partial melt bodies between The Geysers and Mount Hannah. Although our inversions cannot constrain the exact depth and geometry of these bodies, the recovered models indicate the existence of a partial melt zone between 6 and 13 km depth.The prolonged eruption history of the CLVF, coupled with the compositional variation of erupted rocks over time and space, is consistent with the existence of several, potentially ephemeral, melt-bearing bodies as opposed to one large melt body. Given the density and location of the recovered anomaly, rhyolite-MELTS thermodynamic modeling suggests the existence of 10–30% rhyodacitic melt within the proposed silicic magma reservoir at about 700 °C and 8 km depth (210 MPa). Independent petrologic, geochemical, and seismic evidence indicates that this silicic partial melt zone is underlain by basaltic melt in the lower to middle crust (13 to 21 km depth), which is fed by a mantle source.Eruptions in the past ∼8.5–13.5 thousand years; high regional heat flow; 3He enrichment of hydrothermal fluids; and our modeling, which suggests the presence of a mid-crustal, silicic partial melt zone, point to a still-active CLVF. The relatively low estimates of partial melt (10–30%) predicted by thermodynamic modeling indicates that an injection of new magma into the imaged partial melt zone is needed to generate sufficient melt to incite future eruptions. Despite the low percent melt estimates within the proposed silicic partial melt zone the potential for future volcanic eruption remains. Due to the proximity of the CLVF to cities surrounding Clear Lake and the densely populated San Francisco Bay Area, continued research and monitoring of the volcanic field are warranted. The geophysical and petrologic modeling presented here improves our understanding of the CLVF magma plumbing system and allows us to better characterize its associated volcanic hazards.

Anomalously High Heat Flow Regions Beneath the Transantarctic Mountains and Wilkes Subglacial Basin in East Antarctica Inferred From Curie Depth

AbstractThe Transantarctic Mountains (TAMs) separate the warmer lithosphere of the Cretaceous‐Tertiary West Antarctic rift system and the colder and older provinces of East Antarctica. Low velocity zones beneath the TAM imaged in recent seismological studies have been interpreted as warm low‐density mantle material, suggesting a strong contribution of thermal support to the uplift of the TAM. We present new Curie Point Depth (CPD) and geothermal heat flow (GHF) maps of the northern TAM and adjacent Wilkes Subglacial Basin (WSB) based exclusively on high resolution magnetic airborne measurements. We find shallow CPD and high GHF beneath the northern TAM, reinforcing the hypothesis of thermal support of the topography of the mountain range. Additionally, this study demonstrates, that limiting spectral analysis to areas with a high density of aeromagnetic measurements increases the resolution of CPD estimates revealing localized shallow CPD and associated high heat flow in the Central Basin of the WSB and the Rennick Graben (RG). Across the study area the CPD ranges from 15 to 35 km and the GHF values range from 30 to 110 mW/m2. The recovered CPD range is compatible with recent Moho depth estimates, as the CPD predominantly lies within the crust, rather than in the magnetite‐poor mantle. GHF estimates, based on the CPD estimates, show a good agreement to sparse in situ GHF measurements and the location of active volcanoes. Comparison to existing continent‐wide GHF estimates shows strong differences from magnetically‐derived heat flow estimates, while seismologically‐derived heat flow estimates show the best agreement to our results.

Lithospheric strength and stress revisited: Pruning the Christmas tree

Whether lithospheric stress can reach the maximum level predicted by the Christmas-tree strength envelope is a fundamental question but with controversial answers. There is little controversy that a deforming lithosphere in high heat flow regions is likely critically stressed, i.e., at full yield at all depths, as described by the envelope. But different conceptual frameworks offer opposite views for very cold lithosphere, either at full yield or far below yield. Here, we use simple numerical models to investigate stresses in end-member cold cratonic lithosphere (e.g., Canadian Craton) in comparison with end-member warm plate-boundary lithosphere (e.g., Canadian Cordillera). The two key elements of our modelling are (1) that lithospheric stress builds up elastically with horizontal tectonic loading not only in the elastic–frictional brittle regime but also in the viscoelastic ductile regime, and (2) that the stress level is limited by the available tectonic force. In a cratonic lithosphere, the limiting tectonic force is sustained by competent rock material over a large depth range, represented by the competent thickness Tc that exceeds 90 km. The lithosphere undergoes mostly elastic deformation at a stress level of a few tens of MPa. While weakly stressed strong lithosphere can still produce limited earthquakes at shallow depths due to structural and stress heterogeneity, the lithospheric stress under horizontal tectonic loading is theoretically predicted to be orders of magnitude lower than predicted by the Christmas-tree envelope. Stresses in a real lithosphere may substantially deviate from this theoretical level because of spatiotemporal variations in rheology and structure. For example, the stress memory of past loading history in cold lithosphere may or may not be erased by more recent tectonic stresses. Because much of previous scientific debates on lithospheric stress levels and comparison with seismicity were focused on topographically induced flexural stress, we also investigate the effect of vertical loading. We show that the effective elastic thickness Te derived from the flexural response is a reasonable proxy for Tc derived from horizontal tectonic loading; a very large Te such as > 80 km is generally associated with very low tectonic stress far below yield. However, the flexure-induced bending stress is not directly comparable with seismicity because it may either enhance or suppress seismogenic stress in the crust.

Wet roots of high elevation in the western United States

Lithospheric heat transfer strongly modulates the distributions of melt, rock strength, and buoyancy that are responsible for volcanism, seismicity, and elevation. Geotherm models often extrapolate shallow (<3 km) heat flow measurements but are complicated by near-surface hydrology, poorly-known crustal thermal properties, and deep advective transfer by melts and volatiles. Here we compare temperatures estimated from P-wave velocities in the uppermost mantle to those modeled from surface heat flow in the western United States. We show that U.S. Cordilleran regions of high heat flow and high elevation have deep temperatures much lower than predicted by steady-state conductive cooling models. We hypothesize that the discrepancy reflects reaction thermodynamics and advection by migration of volatiles and melts up the lithospheric column. Hydration of the mantle and lower crust by Farallon subduction consumes garnet into melts that absorb latent heat in the lower crust and upper mantle, while hydration reaction enthalpy heats the subduction back-arc mantle, and advection amplifies surface heat flow. This process would increase elevation both by raising temperatures and by converting dense constituents of the mineral assemblage to more buoyant intrusions in the middle and upper crust. The results imply hydration reaction enthalpy changes play a significant role in the dynamics of Cordilleran regions.

Catalog of Enhanced Geothermal Systems based on Heat Sources

AbstractIt is common sense that a deeper well implies higher temperature in the exploration of deep geothermal resources, especially with hot dry rock (HDR) geothermal resources, which are generally exploited in terms of enhanced geothermal systems (EGS). However, temperature is always different even at the same depth in the upper crust due to different heat sources. This paper summarizes the heat sources and classifies them into two types and five sub‐types: crust‐origin (partial melting, non‐magma‐generated tectonic events and radiogenic heat production), and mantle‐origin (magma and heat conducted from the mantle). A review of global EGS sites is presented related to the five sub‐types of heat sources. According to our new catalog, 71% of EGS sites host mantle‐origin heat sources. The temperature logging curves indicate that EGS sites which host mantle‐origin magma heat sources have the highest temperature. Therefore, high heat flow (&gt;100 mW/m2) regions with mantle‐origin magma heat sources should be highlighted for the future exploration of EGS. The principle to identify the heat source is elucidated by applying geophysical and geochemical methods including noble gas isotope geochemistry and lithospheric thermal structure analysis. This analytical work will be helpful for the future exploration and assessment of HDR geothermal resources.

Ambient Noise Tomography of Upper Crustal Structures and Quaternary Faults in the Seoul Metropolitan Area and Its Geological Implications

AbstractWe investigate the upper‐crustal seismic velocity structure in the Seoul metropolitan area, where about 20 million people live. The Chugaryeong fault zone (CFZ) is placed in this area, but the seismic hazard potential remains unclear. We conducted ambient noise tomography to illuminate the high‐resolution upper‐crustal structure in the Seoul metropolitan area. We analyzed continuous vertical seismic records for ∼5 months from a dense seismic array with 77 broadband stations. Group velocity dispersion curves and tomographic maps were extracted between 0.5 and 10 s periods. We inverted 3‐D group velocity tomography models up to a depth of ∼10 km from the group velocity maps. The shear‐wave velocity model is consistent with the geological features. High‐velocity anomalies at shallow depths are correlated with the surface topography and geology. The CFZ is located at a low velocity below the 5 km depth and presented as the simplified model. The large VS contrast regions are located beneath NS‐trending faults. The cross‐sections coincide with the near‐vertical strike‐slip faults in this area. In the southern region of the Seoul metropolitan area, low‐velocity anomalies correlate with high heat flow regions. Our results effectively suggest high resolution upper‐crustal structures and subsurface hidden faults in the urban area.

Ancient deep roots for Mesozoic world-class gold deposits in the north China craton: An integrated genetic perspective

The North China Craton (NCC) hosts some of the world-class gold deposits that formed more than 2 billion years after the major orogenic cycles and cratonization. The diverse models for the genesis of these deposits remain equivocal, and mostly focused on the craton margin examples, although synchronous deposits formed in the interior domains. Here we adopt an integrated geological and geophysical perspective to evaluate the possible factors that contributed to the formation of the major gold deposits in the NCC. In the Archean tectonic framework of the NCC, the locations of the major gold deposits fall within or adjacent to greenstone belts or the margins of micro-continents. In the Paleoproterozoic framework, they are markedly aligned along two major collisional sutures – the Trans North China Orogen and the Jiao-Liao-Ji Belt. Since the Mesozoic intrusions hosting these deposits do not carry adequate signals for the source of gold, we explore the deep roots based on available geophysical data. We show that the gold deposits are preferentially distributed above zones of uplifted MOHO and shallow LAB corresponding to thinned crust and eroded sub-lithospheric mantle, and that the mineralization is located above regions of high heat flow representing mantle upwelling. The NCC was at the center of a multi-convergent regime during the Mesozoic which intensely churned the mantle and significantly enriched it. The geophysical data on Moho and LAB upwarp from the centre towards east of the craton is more consistent with paleo-Pacific slab subduction from the east exerting the dominant control on lithospheric thinning. Based on these results, and together with an evaluation of the geochemical and isotopic features of the Mesozoic magmatic intrusions hosting the gold mineralization, we propose a genetic model that invokes reworking of ancient Au archives preserved in the lower crust and metasomatised upper mantle and which were generated through multiple subduction, underplating and cumulation events associated with cratonization of the NCC as well as the subduction-collision of Yangtze Craton with the NCC. The heat and material input along zones of heterogeneously thinned lithosphere from a rising turbulent mantle triggered by Mesozoic convergent margins surrounding the craton aided in reworking the deep roots of the ancient Au reservoirs, leading to the major gold metallogeny along craton margins as well as in the interior of the NCC.

Gridded maps of geological methane emissions and their isotopic signature

Abstract. Methane (CH4) is a powerful greenhouse gas, whose natural and anthropogenic emissions contribute ∼20 % to global radiative forcing. Its atmospheric budget (sources and sinks), however, has large uncertainties. Inverse modelling, using atmospheric CH4 trends, spatial gradients and isotopic source signatures, has recently improved the major source estimates and their spatial–temporal variation. Nevertheless, isotopic data lack CH4 source representativeness for many sources, and their isotopic signatures are affected by incomplete knowledge of the spatial distribution of some sources, especially those related to fossil (radiocarbon-free) and microbial gas. This gap is particularly wide for geological CH4 (geo-CH4) seepage, i.e. the natural degassing of hydrocarbons from the Earth's crust. While geological seepage is widely considered a major source of atmospheric CH4, it has been largely neglected in 3-D inverse CH4 budget studies given the lack of detailed a priori gridded emission maps. Here, we report for the first time global gridded maps of geological CH4 sources, including emission and isotopic data. The 1∘×1∘ maps include the four main categories of natural geo-CH4 emission: (a) onshore hydrocarbon macro-seeps, including mud volcanoes, (b) submarine (offshore) seeps, (c) diffuse microseepage and (d) geothermal manifestations. An inventory of point sources and area sources was developed for each category, defining areal distribution (activity), CH4 fluxes (emission factors) and its stable C isotope composition (δ13C-CH4). These parameters were determined considering geological factors that control methane origin and seepage (e.g. petroleum fields, sedimentary basins, high heat flow regions, faults, seismicity). The global geo-source map reveals that the regions with the highest CH4 emissions are all located in the Northern Hemisphere, in North America, in the Caspian region, in Europe and in the East Siberian Arctic Shelf. The globally gridded CH4 emission estimate (37 Tg yr−1 exclusively based on data and modelling specifically targeted for gridding, and 43–50 Tg yr−1 when extrapolated to also account for onshore and submarine seeps with no location specific measurements available) is compatible with published ranges derived using top-down and bottom-up procedures. Improved activity and emission factor data allowed previously published mud volcanoes and microseepage emission estimates to be refined. The emission-weighted global mean δ13C-CH4 source signature of all geo-CH4 source categories is about −49 ‰. This value is significantly lower than those attributed so far in inverse studies to fossil fuel sources (−44 ‰) and geological seepage (−38 ‰). It is expected that using this updated, more 13C-depleted, isotopic signature in atmospheric modelling will increase the top-down estimate of the geological CH4 source. The geo-CH4 emission grid maps can now be used to improve atmospheric CH4 modelling, thereby improving the accuracy of the fossil fuel and microbial components. Grid csv (comma-separated values) files are available at https://doi.org/10.25925/4j3f-he27.

Terrestrial heat flow and crustal thermal structure of the Gonghe-Guide area, northeastern Qinghai-Tibetan plateau

The Gonghe-Guide area, which is located at the northeastern edge of the Qinghai-Tibetan Plateau, has the greatest hot dry rock (HDR) geothermal resources exploration and development potential in China. Separated by the Waligong tectonomagmatic belt, the study area includes the Gonghe and Guide basins, which are characterized by distinctly different thermal state and surface thermal manifestations. Although a number of geophysical explorations have been undertaken, studies on fundamental geothermal theory remain scarce, including a lack of high-quality heat flow determinations and thermal structure studies. In this study, we obtained continuous steady-state temperature logs from three deep boreholes (GR2 and DR3 in the Gonghe basin and ZR1 in the Guide basin). The thermal conductivities and radiogenic heat production rates of the outcrop samples gathered from the two basins were measured. The heat-flow values were determined to be 116.3–122.3mWm−2 for the Gonghe basin, yielding a mean of 119.3mWm−2. Two existed heat flow sites (R2 and R3 boreholes) gave an average heat flow of 75mWm−2 for the Guide basin, which is approximately that of the regional background heat flow level of the northeastern Tibetan Plateau. The 1-D theoretical crustal thermal structures of the Gonghe and Guide basins were constructed by solving steady-state heat conduction equation. The results show that the Moho temperature and mantle heat flow in the Guide basin were 1016°C and 27mWm−2, respectively. However, the temperature field in the Gonghe basin indicated partial melting in the crust. The high temperatures and thus heat flow in the Gonghe basin may suggest the existence of anomalous heat source body in crust. Successful inversions of vertical and horizontal temperature variations of boreholes in the Gonghe basin were obtained by numerically simulating a cooling rectangular magma chamber with a buried depth of 5–6km and horizontal dimensions of 12–14km. The study emphasizes that the anomalous heat flow in the Gonghe basin may be the joint thermal effect of: (1) the high regional background heat flow level of the northeastern Tibetan Plateau, which is mainly due to the radiogenic heat of thickened crust; (2) the cooling of a shallow magma chamber, which may be caused by the uplift and denudation of the plateau.

Mapping capable faulting hazard in a moderate-seismicity, high heat-flow environment: The Tuscia province (southern Tuscany-northern Latium, Italy)

This study is focused on potential evidence of capable faulting affecting Southern Tuscany-Northern Latium, a region of high heat flow and moderate to low seismicity. Here we refer to this region as Tuscia, the Latin name for the land of the Etruscans.The review of the evidence for active tectonics in this area has allowed us to identify several potentially active/capable faults, mapped using a combination of analysis of digital terrain models with remote sensing and local geological information derived from bibliographic data and different scale maps.Active and capable faulting is certainly more evident in the eastern sector of Tuscia (Chianti-Cetona Ridge, Pienza, Amiata and Bolsena caldera). However, even in the western Tuscia (e.g. Colline Metallifere and Tyrrhenian coast), capable faulting may be hypothesized along numerous structural elements.The characterization of active faulting potential, in terms of geometry, slip types and rates, in such a high heat-flow environment, beyond its usefulness in improving the assessment of fault displacement and earthquake hazard, also helps to improve knowledge of the stress field and structural setting of the many geothermal fields in the area. This allows new useful insights into geothermal exploitation and hazard (e.g., induced-triggered seismicity). Moreover, this new knowledge of active tectonics and capable faults of the hinterland sector of the Northern Apennines may also provide useful insights into the siting processes of highly hazardous and strategic structures and, more generally, in land use planning and management of infrastructures, especially during emergencies.

Detection of earthquake induced radon precursors by Hilbert Huang Transform

Continuous measurement of radon-222 concentration in soil was carried out across duration of one year at a geologically faulted area having high regional heat flow to detect anomalies caused by seismic activities. The data reveals a range of periodicities present in the radon time series. To identify seismic induced radon changes we treat the time series data through various filtering methods to remove inherent periodicities. The Ensemble Empirical Mode Decomposition (EEMD) is deployed to decompose the signal into its characteristic modes. Hilbert Huang Transform (HHT) is applied for the first time on the physically significant modes obtained by EEMD to represent time–energy–frequency of the recorded soil radon time series. After removing the periodic and quasi-periodic constituents from the original time series, the simulated result shows a forceful correlation in recorded radon-222 anomalies with regional and local seismic events.

Crust beneath the northwestern Deccan Volcanic Province, India: Evidence for uplift and magmatic underplating

AbstractThe northwestern Deccan Volcanic Province of India and its pericratonic rift basins were reactivated during different stages after the breakup of India from the Gondwanaland and collision with the Asian plate. In this study, we present results of crustal thickness and average crustal Vp/Vs ratios beneath this plume‐affected region using common conversion point imaging and H‐k stacking analysis of 6893 receiver functions using data from a network comprising 58 broadband seismic stations sited on diverse tectonic terrains. We find large variations in crustal thickness, with the Moho depths varying from 28 to 43 km in the Kachchh rift, 28 to 38 km in the Cambay rift, 39.5–41.5 km in the north and eastern parts of the Cambay rift, and 29 to 39 km in the Saurashtra region and South Gujarat. A Moho upwarp of 6 to 7 km in the Saurashtra region can be attributed to positive buoyancy and uplift due to thermal influx affected by the Reunion plume. High crustal Vp/Vs ratios beneath the Kachchh rift (1.8 to 2.05), coastal areas of Saurashtra (1.75 to 2.06), and North Gujarat (1.81 to 1.85) indicate dominance of a mafic/ultramafic crust. High regional heat flow, high electrical conductivity, large intracrustal S wave velocity reduction, and high average crustal Poisson's ratios are consistent with partial melt related to the process of magmatic underplating in the lower crust. At other stations, the crust appears to be felsic with Vp/Vs ratios in the range of 1.57 to 1.76.

Silica Geothermometry Applications in the Taiwan Orogenic Belt

The silica concentrations from 130 hot spring measurements were calculated and used to indicate the heat flow distribution in the active Taiwan orgenic belt. The heat flow spatial distribution seems to be controlled by the Taiwan arc-continent collision zone tectonic development. Two abnormally high silica heat flow regions, 160 and 190 mW m^(-2), are identified in northeastern and southeastern Taiwan, respectively. In northeastern Taiwan, the Chingshui area, the data suggests that the anomalous heat flow distribution may be controlled mainly by hot fluids and faults. The highest heat flow in southeastern Taiwan may be generated by heat advection caused by a rapidly uplifting hot crust with a decreasing mature arc-continent collision zone exhumation rate. In addition, we estimated the exhumation rate in different tectonic zones in Southern Taiwan by applying fission track ages and silica heat flow values. The exhumation rates are 1.33, 1.72-3.87, and 0.51-1.74 mm yr^(-1) in the tectonic zones around Taiwan for advanced, initial arc-continent collision and accretionary wedge deformation, respectively.

Estimation of depth to top of magnetic sources using the local-wavenumber approach in an area of shallow Moho and Curie depth — The Red Sea

The local-wavenumber method estimates the depth to a magnetic source based on the spectral content of a single anomaly assuming that the base of the magnetic body is at infinite depth. However, the “infinite-depth” assumption can lead to significant underestimation of the depth to the top of magnetic bodies, especially in areas where the depth to the bottom of the magnetic layer is not large compared to the depth to the top, as would occur in high heat-flow regions and thinned continental margins. Such underestimation of depths has been demonstrated in model studies and using real data with seismic and well control. We evaluated a modification to the local-wavenumber approach to estimate the depth to the top of magnetic sources assuming that the depth to the bottom of the magnetic sources is controlled by the Curie temperature or crustal thickness. We applied this new method to a simple model of a continental margin and to magnetic survey data over the central Red Sea where the Curie isotherm is shallow. The effective structural index of this finite depth extent model is found to increase continuously from the continent to the ocean as the depth to the magnetic basement increases and the depth to the bottom of the magnetic layer decreases. We have also discovered in this study that the local-wavenumber maxima correlate well with major seafloor spreading magnetic reversal epochs in the central Red Sea segment.

Heat Flow, Heat Generation, and the Thermal State of the Lithosphere

The thermal structure of the lithosphere controls many properties and processes of Earth's crust. The total ∼47-TW heat loss of Earth is key to understanding and modeling this thermal structure, as is partitioning the various sources of that heat into heat entering the base of the lithosphere, heat generated within the lithosphere by radioactive decay (primarily within the continental crust), and secular cooling of the mantle lithosphere (primarily in oceanic lithosphere). A set of framework geotherms for the continental lithosphere explains deep crustal melting in high heat flow regions, metamorphic pressure-temperature (P-T) space in the crust, partial melting at the base of the lithosphere to produce an S-wave low-velocity zone in Phanerozoic and younger terranes, and the P-T fields inferred from mantle xenoliths. Important perturbations to a standard thermal state are produced by orogenic overprints, transient thermal regimes, and exhumation.