6.2 Geothermal gradients and heat flow measurement

Geothermal gradients measure how Earth's temperature changes with depth. This concept is central to understanding heat flow from the planet's interior to its surface, and it directly informs geothermal energy exploration, hydrocarbon maturation studies, and models of Earth's thermal evolution.

Measuring geothermal gradients involves recording temperatures in boreholes and analyzing the thermal properties of the rocks those boreholes penetrate. Combining these two datasets through Fourier's law yields heat flow values that reveal subsurface thermal patterns.

Geothermal Gradient and Its Significance

Definition and Measurement Units

The geothermal gradient is the rate at which Earth's temperature increases with depth. It's typically expressed in °C/km (or sometimes °F/100 ft in older literature). A steeper gradient means temperature rises faster as you go deeper.

Average Gradient and Variability

The average crustal geothermal gradient is roughly 25–30 °C/km, but this number varies a lot depending on location. Tectonic setting, rock type, and subsurface fluid circulation all play a role. In volcanic regions like Iceland, gradients can exceed 100 °C/km, while stable continental interiors may sit well below the average.

Driving Factors and Heat Sources

Two primary sources supply Earth's internal heat:

• Primordial heat from planetary accretion and core formation, left over from Earth's formation ~4.5 billion years ago
• Radiogenic heat from the decay of radioactive isotopes, primarily uranium ($U$), thorium ($Th$), and potassium ($^{40}K$)

Radiogenic heat production is especially significant in felsic crustal rocks. Granitic intrusions, for example, can have elevated heat production rates (on the order of 2–5 μW/m³), which raises local geothermal gradients above the regional background.

Importance in Geophysical Applications

Geothermal gradient data feeds into several applied problems:

• Geothermal energy exploration: identifying high-temperature reservoirs suitable for power generation
• Hydrocarbon maturation studies: determining whether source rocks have reached the thermal window for oil or gas generation
• Tectonic and volcanic analysis: mapping lithospheric thermal structure and identifying active heat sources
• Basin modeling: constraining thermal histories and calibrating numerical simulations of sedimentary basin evolution

Measuring Geothermal Gradients and Heat Flow

Borehole Temperature Measurements

Direct temperature measurement in boreholes is the most reliable way to determine the geothermal gradient. Here's how it works:

1. Drill or access a borehole that penetrates to the depth of interest.
2. Allow the borehole to thermally equilibrate. Drilling circulates fluid that disturbs the natural temperature profile, so the well needs time (often weeks to months) to return to equilibrium.
3. Lower a temperature logging tool (typically a thermistor or thermocouple sensor) down the borehole, recording temperature at regular depth intervals.
4. Plot temperature vs. depth. The slope of this profile gives you the geothermal gradient.

Differential temperature logging is a variant that measures the temperature difference between two closely spaced sensors. This approach is more sensitive to small gradient changes and helps identify interval-specific variations that a single-sensor log might smooth over.

Thermal Conductivity Measurements

Knowing the gradient alone isn't enough to calculate heat flow. You also need to know how easily heat moves through the rock, which is its thermal conductivity ($k$).

Thermal conductivity is measured in the lab on core samples or cuttings using methods such as:

• Divided bar method: a steady-state technique where a rock disc is sandwiched between materials of known conductivity, and heat flow through the stack is measured
• Needle probe method: a transient technique where a heated needle is inserted into the sample and the rate of temperature rise is recorded

Once you have both gradient and conductivity, you calculate heat flow using Fourier's law of heat conduction:

$Q = -k \frac{dT}{dz}$

• $Q$ = heat flow (mW/m²)
• $k$ = thermal conductivity (W/m·K)
• $\frac{dT}{dz}$ = geothermal gradient (°C/m or K/m)

The negative sign indicates heat flows in the direction of decreasing temperature (from hot to cold, i.e., upward toward the surface).

Estimations from Well Logs

When direct temperature measurements aren't available, geothermal gradients can be estimated indirectly from well logs. Resistivity logs and acoustic (sonic) logs are both sensitive to temperature because temperature affects pore fluid properties and rock matrix behavior. These estimates are less precise than direct measurements but useful when equilibrated temperature logs don't exist.

Heat Flow Measurements

Heat flow values, reported in mW/m², combine gradient and conductivity data into a single quantity that characterizes the thermal regime of a region. Global continental heat flow averages around 65 mW/m², while oceanic heat flow averages about 100 mW/m² (though it varies strongly with crustal age).

Heat flow maps are used to:

• Identify thermal anomalies that may indicate geothermal resources or active tectonic processes
• Constrain regional thermal models
• Compare thermal regimes across different geological provinces

Factors Influencing Geothermal Gradients

Lithology and Thermal Properties

Different rock types conduct heat at different rates, and this directly affects the geothermal gradient. The relationship is inverse: for a given heat flow, rocks with lower thermal conductivity will have higher geothermal gradients, and vice versa.

• Sedimentary rocks (shales, mudstones) generally have low thermal conductivities (1–2 W/m·K), producing steeper gradients in sedimentary basins.
• Salt and quartzite are highly conductive (5–7 W/m·K for salt), so they transmit heat efficiently and produce lower local gradients. Salt diapirs, for instance, act as thermal chimneys, funneling heat upward and creating distinctive thermal anomalies.
• Igneous and metamorphic rocks typically fall between these extremes, though composition matters (felsic rocks tend to be less conductive than mafic rocks).

Tectonic Setting

The tectonic environment exerts first-order control on heat flow and geothermal gradients:

• Extensional settings (continental rifts, back-arc basins) tend to have elevated gradients because lithospheric thinning brings hot asthenospheric material closer to the surface. The Basin and Range Province in the western U.S. is a classic example, with heat flow values commonly exceeding 80–100 mW/m².
• Convergent margins display complex thermal patterns. The forearc above a subducting slab is typically cold (the descending plate depresses isotherms), while the volcanic arc behind it shows elevated heat flow from magmatic activity.
• Stable cratons generally have low, uniform heat flow (40–50 mW/m²) due to their thick, old lithosphere.

Fluid Circulation Effects

Subsurface fluid movement can redistribute heat far more efficiently than conduction alone, and this is one of the biggest complications in interpreting geothermal data:

• Upward fluid migration (e.g., hydrothermal convection) transports heat toward the surface, creating localized hot spots and steepening the near-surface gradient.
• Downward fluid flow (e.g., meteoric water recharge) carries cool surface water to depth, depressing the gradient.
• Convective heat transfer by circulating fluids can dominate over conduction in fractured or permeable zones, making it difficult to extract the purely conductive thermal signal.

Recognizing the signature of advective (fluid-driven) vs. conductive heat transfer is critical when interpreting borehole temperature profiles.

Radiogenic Heat Production

Rocks enriched in $U$, $Th$, and $^{40}K$ generate heat internally through radioactive decay. This effect is most pronounced in:

• Upper-crustal granitic bodies
• Felsic gneisses
• Some sedimentary formations enriched in organic matter or heavy minerals

In regions with thick, radiogenic upper crust, the contribution to surface heat flow can be substantial (sometimes 30–50% of the total).

Interpreting Geothermal Data for Subsurface Characterization

Identifying Thermal Anomalies

Mapping geothermal gradients and heat flow across a region reveals areas that deviate from the background thermal field. Positive anomalies (higher than expected) may indicate:

• Shallow magmatic intrusions or hydrothermal systems
• Geothermal reservoirs suitable for energy extraction
• Upward fluid migration pathways

Negative anomalies can point to downward fluid flow, recent sedimentation (which buries and insulates the surface), or the presence of highly conductive lithologies redistributing heat laterally.

Constraining Basin Thermal History

Geothermal data integrates with other geological and geophysical datasets (vitrinite reflectance, fission track thermochronology, burial history curves) to reconstruct how a sedimentary basin's temperature field has evolved over geologic time. This thermal history is essential for predicting:

• When and where source rocks entered the oil and gas generation windows
• The timing and pathways of hydrocarbon migration
• The degree of diagenetic alteration in reservoir rocks

Estimating Depth to the Brittle-Ductile Transition

The brittle-ductile transition is the depth at which crustal rocks shift from fracturing under stress to flowing plastically. Its depth depends strongly on the geothermal gradient because rock strength decreases with temperature. A steeper gradient pushes this transition shallower, which has implications for:

• Maximum depth of seismicity (earthquakes generally nucleate in the brittle zone)
• Fluid flow pathways (fracture permeability dominates above the transition)
• Crustal deformation style (faulting vs. ductile shear)

Geothermal Resource Assessment

Heat flow and gradient data are the starting point for evaluating geothermal energy potential. High gradients mean you can reach economically useful temperatures (typically >150 °C for electricity generation) at shallower, more drillable depths. Lower-temperature resources (50–100 °C) can still be valuable for direct-use applications like district heating, greenhouse agriculture, and aquaculture.

Calibrating Basin and Thermal Models

Measured geothermal gradients and heat flow values serve as ground truth for numerical models of basin evolution and lithospheric thermal structure. These measurements help:

• Set appropriate thermal boundary conditions (e.g., basal heat flow at the base of the sedimentary section)
• Validate model predictions against observed temperature-depth profiles
• Reduce uncertainty in forward models used for resource exploration and tectonic reconstruction