On the radiogenic heat production of igneous rocks

Radiogenic heat production is a physical parameter crucial to properly estimating lithospheric temperatures and properly understanding processes related to the thermal evolution of the Earth. Yet heat production is, in general, poorly constrained by direct observation because the key radiogenic elements exist in trace amounts making them difficulty image geophysically. In this study, we advance our knowledge of heat production throughout the lithosphere by analyzing chemical analyses of 108,103 igneous rocks provided by a number of geochemical databases. We produce global estimates of the average and natural range for igneous rocks using common chemical classification systems. Heat production increases as a function of increasing felsic and alkali content with similar values for analogous plutonic and volcanic rocks. The logarithm of median heat production is negatively correlated ( r 2 = 0.98) to compositionally-based estimates of seismic velocities between 6.0 and 7.4 km s −1 , consistent with the vast majority of igneous rock compositions. Compositional variations for continent-wide models are also well-described by a log-linear correlation between heat production and seismic velocity. However, there are differences between the log-linear models for North America and Australia, that are consistent with interpretations from previous studies that suggest above average heat production across much of Australia. Similar log-linear models also perform well within individual geological provinces with ∼1000 samples. This correlation raises the prospect that this empirical method can be used to estimate average heat production and natural variance both laterally and vertically throughout the lithosphere. This correlative relationship occurs despite a direct causal relationship between these two parameters but probably arises from the process of differentiation through melting and crystallization.