Mantle Heat Flow Distribution Calculator
Our planetary & earth system science calculator computes mantle heat flow distribution accurately. Enter your values for instant results.
Reviewed for accuracy by Daniel Agrici, Founder & Lead Developer
Mantle Heat Flow Distribution Calculator
Calculator
Adjust values & calculateEnter your values below. Every result is computed in your browser โ no data is sent to any server.
Formula: q_m = q_s - A x D; q(t) = C / sqrt(t); Q = q_ocean x A_ocean + q_cont x A_cont
Worked example โ Mantle heat flow 35 mW/m2 | Crustal radiogenic 30 mW/m2 | Crust supplies 46.2 percent
Formula
q_m = q_s - A x D; q(t) = C / sqrt(t); Q = q_ocean x A_ocean + q_cont x A_cont
Continental partition (Birch-Roy reduced heat flow): q_m is the mantle heat flow reaching the base of the lithosphere, q_s is the measured surface heat flow in mW/m2, A is crustal radiogenic heat production in microwatts per cubic meter, and D is the thickness of the enriched layer in km. Because 1 uW/m3 acting over 1 km equals exactly 1 mW/m2, the product A x D is already in mW/m2. Oceanic partition: conductive cooling of a semi-infinite half space gives q(t) = k(Tm - T0) / sqrt(pi x kappa x t), which collapses to C / sqrt(t) with C about 510 mW/m2 per root-Myr for k = 3.138 W/(m K), Tm - T0 = 1450 K and kappa = 8.05e-7 m2/s. Global heat loss Q is the area-weighted sum over ocean floor and continents using Earth surface area 4 pi R^2 = 5.101e14 m2 for R = 6371 km. The oceanic area fraction here is the share of the surface floored by oceanic lithosphere, about 0.607 in the Pollack, Hurter and Johnson partition, with continents plus their submerged shelves making up the remaining 0.393; it is deliberately smaller than the 0.71 of the surface merely covered by water, because the shelves are continental crust.
Worked Examples
Example 1: Continental Partition into Crustal and Mantle Heat Flow
Problem:A shield site has a measured surface heat flow of 65 mW/m2. Outcrop samples give a radiogenic heat production of 3.0 microwatts per cubic meter over an enriched layer 10 km thick. Split the surface flux into its crustal and mantle parts.
Solution:Unit identity: 1 uW/m3 x 1 km = 1e-6 W/m3 x 1e3 m = 1e-3 W/m2 = 1 mW/m2 Crustal radiogenic flux q_c = A x D = 3.0 x 10 = 30 mW/m2 Mantle (reduced) heat flow q_m = q_s - q_c = 65 - 30 = 35 mW/m2 Crustal share = 30 / 65 = 0.4615 = 46.2 percent
Result:Mantle heat flow 35 mW/m2 | Crustal radiogenic 30 mW/m2 | Crust supplies 46.2 percent
Example 2: Oceanic Heat Flow Versus Seafloor Age
Problem:Using the half-space cooling constant C = 510 mW/m2 per root-Myr, compare the conductive heat flow on 4 Myr ridge flank, 60 Myr crust, and 100 Myr crust.
Solution:q(t) = C / sqrt(t), with t in Myr and q in mW/m2 q(4) = 510 / sqrt(4) = 510 / 2 = 255.0 mW/m2 q(60) = 510 / sqrt(60) = 510 / 7.746 = 65.8 mW/m2 q(100) = 510 / sqrt(100) = 510 / 10 = 51.0 mW/m2 The flux falls by exactly a factor of 5 from 4 Myr to 100 Myr because sqrt(100/4) = 5
Result:255.0 mW/m2 at 4 Myr | 65.8 mW/m2 at 60 Myr | 51.0 mW/m2 at 100 Myr
Example 3: Global Heat Loss Integrated Over Oceans and Continents
Problem:Integrate to a planetary total using an oceanic area fraction of 0.607, an oldest seafloor age of 180 Myr, C = 510 mW/m2 per root-Myr, and a continental surface heat flow of 65 mW/m2.
Solution:Earth surface area = 4 pi R^2 = 4 pi (6.371e6 m)^2 = 5.101e14 m2 Ocean area = 0.607 x 5.101e14 = 3.096e14 m2; continent area = 0.393 x 5.101e14 = 2.005e14 m2 Area-weighted oceanic mean, triangular age distribution = (8/3) C / sqrt(T) = (8/3)(510) / sqrt(180) = 1360 / 13.416 = 101.4 mW/m2 Ocean heat loss = 0.1014 W/m2 x 3.096e14 m2 = 3.14e13 W = 31.4 TW Continent heat loss = 0.065 W/m2 x 2.005e14 m2 = 1.30e13 W = 13.0 TW Global total = 31.4 + 13.0 = 44.4 TW Global mean flux = 4.44e13 W / 5.101e14 m2 = 0.0871 W/m2 = 87.1 mW/m2
Result:Ocean 31.4 TW | Continents 13.0 TW | Global 44.4 TW at a mean flux of 87.1 mW/m2
Frequently Asked Questions
What is geothermal heat flow and how is it measured?
Geothermal heat flow is the rate at which thermal energy escapes from Earth's interior through its surface, expressed in milliwatts per square meter (mW/m2). It is measured by drilling boreholes and recording the temperature gradient with depth alongside the thermal conductivity of the rock. Heat flow equals the product of thermal conductivity and the temperature gradient, which is Fourier's law of heat conduction: q equals negative k times dT/dz. Global average surface heat flow is approximately 87 mW/m2, representing the combined contribution of radiogenic heat production in the crust and primordial heat left from Earth's formation and differentiation.
What is Fourier's law of heat conduction and how does it apply to the mantle?
Fourier's law states that the heat flux through a material is proportional to the negative of the temperature gradient and the thermal conductivity of the material: q equals negative k times dT/dx, where q is heat flux in W/m2, k is thermal conductivity in W/(m K), and dT/dx is the temperature gradient in K/m. In the context of Earth's mantle, this law governs conductive heat transport through lithospheric plates. The mantle itself transfers heat primarily by slow convection rather than conduction, but the rigid lithosphere above it conducts heat conductively to the surface. Typical mantle thermal conductivity ranges from 3 to 4 W/(m K) for peridotite at relevant pressures.
How is continental heat flow split between crustal radiogenic heat and mantle heat flow?
Surface heat flow measured on a continent combines heat generated inside the crust by radioactive decay with heat conducted into the base of the lithosphere from the mantle below. Roy, Blackwell and Birch showed in 1968 that heat flow measured across a tectonic province falls on a straight line when plotted against the radiogenic heat production of the surface rocks, giving q_s equals q_m plus A times D. Here q_m is the intercept, known as the reduced or mantle heat flow, A is the heat production of the near-surface rocks in microwatts per cubic meter, and D is the characteristic thickness of the enriched layer, typically 7 to 16 km. The units work out cleanly because one microwatt per cubic meter acting over one kilometer is exactly one milliwatt per square meter. Their original fits gave a reduced heat flow of about 33 mW/m2 with D near 7.5 km for the eastern United States and about 17 mW/m2 with D near 10 km for the Sierra Nevada.
How does oceanic heat flow depend on seafloor age?
Oceanic lithosphere is created hot at a spreading ridge and cools by conduction as it moves away, so its heat flow depends on age rather than on distance or spreading rate. Treating the cooling plate as a semi-infinite half space held at the seafloor temperature gives q(t) equal to the thermal conductivity times the temperature contrast, divided by the square root of pi times the thermal diffusivity times the age. Collecting the constants gives the compact working form q(t) equals C divided by the square root of the age in millions of years. With a thermal conductivity of 3.138 W/(m K), a mantle-to-surface temperature contrast of 1450 K and a thermal diffusivity of 8.05 x 10^-7 m2/s, C evaluates to about 510 mW/m2 per root million years, the value used in the GDH1 model of Stein and Stein. Earlier fits gave 473 in Parsons and Sclater and 480 in Lister and colleagues. The relation predicts 255 mW/m2 on 4-million-year-old crust and 51 mW/m2 on 100-million-year-old crust. Beyond roughly 80 million years the observed flux flattens near 48 mW/m2 instead of continuing to fall, which is why plate models with a fixed basal temperature replace the pure half-space solution for old seafloor.
How does continental heat flow differ from oceanic heat flow?
Oceanic crust has significantly higher average heat flow than continental crust, around 101 mW/m2 compared to approximately 65 mW/m2 for continents, though both values have wide variability. Young oceanic crust near mid-ocean ridges can exhibit heat flow exceeding 200 to 300 mW/m2 because hot mantle rock is close to the surface. As oceanic crust ages and moves away from spreading centers it cools and subsides, reducing heat flow to values around 50 mW/m2 for old ocean basins. Continental heat flow is elevated in active tectonic regions and geothermal areas but is lower in stable cratons that have not been volcanically or tectonically active for hundreds of millions of years.
What role does radioactive decay play in Earth's heat flow?
Radioactive decay is the dominant source of heat generated inside Earth's continental crust, and decay in the crust and mantle together supplies roughly 50 percent of Earth's total heat output of about 47 terawatts. Continental crust by itself accounts for only about 6 to 8 terawatts of that, because the heat-producing elements are concentrated in a thin granitic layer covering less than half the surface; the rest of the radiogenic heat comes from the far larger volume of the mantle. The primary heat-producing isotopes are uranium-238, uranium-235, thorium-232, and potassium-40, all of which undergo spontaneous decay that releases energy as heat. Granitic continental crust is enriched in these elements relative to oceanic crust and the mantle. Typical continental crustal heat production is 1 to 3 microwatts per cubic meter. The remaining heat flow comes from the cooling of the primordial Earth, the latent heat of inner core crystallization, and gravitational energy from core formation that occurred early in Earth's history.
Why is heat flow elevated near mid-ocean ridges and subduction zones?
Mid-ocean ridges exhibit high heat flow because they sit directly above upwelling mantle rock. As tectonic plates diverge, hot asthenosphere wells up to fill the gap, bringing temperatures near the melting point of peridotite to within a few kilometers of the seafloor. This creates heat flow values that can locally exceed 500 mW/m2 near active spreading centers, though hydrothermal circulation through fractured crust redistributes much of this heat. Subduction zones show elevated heat flow on the volcanic arc side due to partial melting of the mantle wedge above the subducting slab, while the slab itself brings relatively cold oceanic crust down into the mantle, creating anomalously low heat flow in the fore-arc region.
What is the global heat budget of Earth's interior?
Earth loses approximately 47 terawatts of heat to space through its surface, with about 32 terawatts escaping through the ocean floor and 15 terawatts through the continents. This total heat loss is divided between two fundamental sources: approximately 50 percent comes from the decay of long-lived radioactive isotopes distributed throughout the mantle and crust, while the other 50 percent is primordial heat inherited from accretion, differentiation, and early radioactive decay of now-extinct short-lived isotopes. This secular cooling is very slow; estimates suggest Earth's interior cools by only about 50 to 100 degrees Celsius per billion years. The ratio of radiogenic to primordial heat, called the Urey ratio, is debated and lies somewhere between 0.4 and 0.7.
How does heat flow relate to the geothermal gradient?
The geothermal gradient is the rate at which temperature increases with depth in Earth's interior, typically expressed in degrees Celsius per kilometer. Near the surface in continental crust the average gradient is about 25 to 30 degrees Celsius per kilometer, though it varies widely with local heat flow and rock type. The relationship between heat flow and gradient is given by Fourier's law: gradient equals heat flow divided by thermal conductivity. A region with high heat flow and low thermal conductivity rock will have a steep gradient, reaching temperatures capable of melting rock at relatively shallow depths. In contrast, stable cratons with low heat flow and high conductivity rocks have shallow gradients and cold, thick lithospheres.
What are the units of heat flow and how large is Earth's heat output?
Heat flow is expressed in milliwatts per square meter (mW/m2), where one milliwatt equals one thousandth of a watt. The global average surface heat flow in the 1993 compilation is approximately 87 mW/m2, which seems small but when integrated over Earth's entire surface area of about 510 million square kilometers yields a total heat output of roughly 44 terawatts, since 87 x 10^-3 W/m2 multiplied by 5.10 x 10^14 m2 gives 4.4 x 10^13 W. The denser 2010 reassessment raises the mean to about 92 mW/m2 and the total to 47 terawatts. For comparison, global human primary energy consumption is roughly 19 terawatts, so Earth's internal heat engine releases more than twice as much power as all human energy use combined. Older literature used the heat flow unit (HFU) where one HFU equals 41.86 mW/m2, but modern publications use SI units exclusively.
References
- Pollack, Hurter & Johnson (1993) - Heat Flow from the Earth's Interior: Analysis of the Global Data Set
- Davies & Davies (2010) - Earth's Surface Heat Flux (Solid Earth)
- Stein & Stein (1992) - A Model for the Global Variation in Oceanic Depth and Heat Flow with Lithospheric Age
- International Heat Flow Commission - Global Heat Flow Database
Background & Theory
History
Reviewed for accuracy by Daniel Agrici, Founder & Lead Developer ยท Editorial policy
Related Calculators
๐งฎGeothermal Heat Flow Calculator
Calculate geothermal heat flow with inputs, formulas, and instant results.
๐งฎHeat Flow From Gradient and Conductivity Calculator
Calculate heat flow from gradient and conductivity with inputs, formulas, and instant results.
๐งฎDarcyโs Law Flow Rate Calculator (geologic)
Calculate darcyโs law flow rate (geologic) with inputs, formulas, and instant results.
๐งฎFlow Duration Curve Calculator
Calculate flow duration curve with inputs, formulas, and instant results.
๐งฎFlow Resistance Converter (nโf) Calculator
Calculate flow resistance converter (nโf) with inputs, formulas, and instant results.
๐งฎGroundwater Flow Velocity Calculator
Calculate groundwater flow velocity with inputs, formulas, and instant results.
๐งฎHydraulic Gradient Calculator (natural Flow)
Calculate hydraulic gradient (natural flow) with inputs, formulas, and instant results.
๐งฎGlacial Flow Velocity (glenโs Law) Calculator
Calculate glacial flow velocity (glenโs law) with inputs, formulas, and instant results.