Spherical harmonic analysis of earth’s conductive heat flow

A reappraisal of the international heat flow database has been carried out and the corrected data set was employed in spherical harmonic analysis of the conductive component of global heat flow. Procedures used prior to harmonic analysis include analysis of the heat flow data and determination of representative mean values for a set of discretized area elements of the surface of the earth. Estimated heat flow values were assigned to area elements for which experimental data are not available. However, no corrections were made to account for the hypothetical effects of regional-scale convection heat transfer in areas of oceanic crust. New sets of coefficients for 12° spherical harmonic expansion were calculated on the basis of the revised and homogenized data set. Maps derived on the basis of these coefficients reveal several new features in the global heat flow distribution. The magnitudes of heat flow anomalies of the ocean ridge segments are found to have mean values of less than 150 mW/m². Also, the mean global heat flow values for the raw and binned data are found to fall in the range of 56–67 mW/m², down by nearly 25% compared to the previous estimate of 1993, but similar to earlier assessments based on raw data alone. To improve the spatial resolution of the heat flow anomalies, the spherical harmonic expansions have been extended to higher degrees. Maps derived using coefficients for 36° harmonic expansion have allowed identification of new features in regional heat flow fields of several oceanic and continental segments. For example, lateral extensions of heat flow anomalies of active spreading centers have been outlined with better resolution than was possible in earlier studies. Also, the characteristics of heat flow variations in oceanic crust away from ridge systems are found to be typical of conductive cooling of the lithosphere, there being little need to invoke the hypothesis of unconfined hydrothermal circulation on regional scales. Calculations of global conductive heat loss, compatible with the observational data set, are found to fall in the range of 29–34 TW, nearly 25% less than the 1993 estimate, which rely on one-dimensional conductive cooling models.     

Reply to comments by Henry N. Pollack and David S. Chapman on “Spherical harmonic analysis of Earth’s conductive heat flow”

Pollack and Chapman, hereafter referred to as P&C, argue that: (1) errors arising from lack of quality control in the IHFC database are not important and not properly documented, (2) resolution of spatial patterns in global heat flux distribution should not be represented by spherical harmonics and (3) heat flow in young oceanic crust and global heat loss are better represented by a contested 1-D cooling model than by the data. We disagree and provide additional information that may help clear up such misunderstandings. We also mention briefly the results of a new improved thermal model of the lithosphere that satisfactorily reproduces the main features identified in observational data sets of heat flow and ocean floor bathymetry. Thus, there is no reason to invoke the ad hoc hypothesis of large-scale hydrothermal circulation in the ocean crust.     

Spectral harmonic analysis and synthesis of Earth’s crust gravity field

We developed and applied a novel numerical scheme for a gravimetric forward modelling of the Earth’s crustal density structures based entirely on methods for a spherical analysis and synthesis of the gravitational field. This numerical scheme utilises expressions for the gravitational potentials and their radial derivatives generated by the homogeneous or laterally varying mass density layers with a variable height/depth and thickness given in terms of spherical harmonics. We used these expressions to compute globally the complete crust-corrected Earth’s gravity field and its contribution generated by the Earth’s crust. The gravimetric forward modelling of large known mass density structures within the Earth’s crust is realised by using global models of the Earth’s gravity field (EGM2008), topography/bathymetry (DTM2006.0), continental ice-thickness (ICE-5G), and crustal density structures (CRUST2.0). The crust-corrected gravity field is obtained after modelling and subtracting the gravitational contribution of the Earth’s crust from the EGM2008 gravity data. These refined gravity data mainly comprise information on the Moho interface and mantle lithosphere. Numerical results also reveal that the gravitational contribution of the Earth’s crust varies globally from 1,843 to 12,010 mGal. This gravitational signal is strongly correlated with the crustal thickness with its maxima in mountainous regions (Himalayas, Tibetan Plateau and Andes) with the presence of large isostatic compensation. The corresponding minima over the open oceans are due to the thin and heavier oceanic crust.     

Reply to D. F. Lascelles. Discussion on ‘Microplaty hematite—its varied nature and genesis’ by R. C. Morris and ‘Genesis modelling for the Hamersley BIF-hosted iron ores of Western Australia: a critical review’ by R. C. Morris and M. Kneeshaw

Abstract

The third (Es₃) member of the Eocene Shahejie Formation is the main hydrocarbon source rock interval in the Raoyang Sag of the Bohai Bay Basin. The lower Es₃ (Es₃^L) and upper Es₃ (Es₃^U) submembers display clear differences in source rock quality. The formation mechanisms of the source rocks are investigated via geochemical methods using 60 samples from three wells to reconstruct the paleo-environment during deposition. The major element parameters exhibit changes in paleoclimate from humid to arid. Indicators such as Sr/Ba, B/Ga, the gammacerane index and isotopic data suggest fresh–brackish and hydrological open lakes with unstable water column stratification likely occurred during deposition of the Es₃^L submember and saline and hydrologically closed lakes with stable water column stratification likely during deposition of the Es₃^U submember. Carbon isotope values of organic matter, trace elements and biomarker parameters suggest that the Es₃^L submember had moderate productivity, with a significant contribution from terrigenous organic matter whereas the Es₃^U submember had slightly enhanced productivity, with no or minor contributions from terrigenous organic matter. Furthermore, the pristane/phytane ratio and the enrichment of Mo and U indicate that the euxinic bottom-water conditions (sulfidic) of the lakes during deposition of the Es₃^U submember were best for preserving organic matter. Comparison of the models of source rock deposition of the Es₃^L and Es₃^U submembers indicates that the redox conditions play an important role in the formation of organic-rich source rocks in the Raoyang Sag.     

Model or measurements? A discussion of the key issue in Chapman and Pollack’s critique of Hamza et al.’s re-evaluation of oceanic heat flux and the global power

Chapman and Pollack (C and P)[2007, Int J Earth Sci] criticize Hamza et al. [2007, Int J Earth Sci] for using actual heat flux measurements in young oceanic crust instead of values from 1-D cooling models. The rationalization of C and P and previous authors is that hydrothermal circulation causes the discrepancy between model and measurement. However, the discrepancy between model values and measured heat flux exists over the entire ocean floor and is opposite to the perturbations that hydrothermal circulation would superimpose on a conductive system [Hofmeister and Criss (2005) Tectonophysics 409:199–203]. The error lies in force-fitting a 1-D cooling model to the 3-D oceanic crust [Hofmeister and Criss (2005) Tectonophysics 395:159–177]. Shortcomings of the 1-D model include mathematical errors, such as use of volumetric rather than linear thermal expansivity to describe contraction which, by assumption, is limited only to the Z -direction [Hofmeister and Criss (2006) Tectonophysics]. This 3× error, traceable to McKenzie and Sclater [1969, Bull Vocanol 33–1:101–118], accidentally provides good agreement of model values with globally averaged seafloor depths for young, but not old ages, and is the sole rationale for using the simplistic cooling model. There is no justification for selective substitution of erroneous 1-D model values for measurements only for the younger half of the 3-D oceanic crust, as stridently and arbitrarily promoted by C and P. Hamza et al. [2007, Int J Earth Sci], in contrast, use the scientific method, which calls for discarding models that do not well describe physical phenomena. The remainder of this report summarizes the shortcomings of cooling models, particularly the half-space cooling (HSC) model touted by C and P, and explains how hydrothermal circulation affects heat flux. We focus on the basics, as these have been misunderstood. With the key issues of C and P being erroneous, it is not necessary to address their remaining comments, many of which enumerate the vote for an imagined, gargantuan circulation of hot fluid through oceanic basins that is somehow warmed without removing heat from the rocks. The use of “consensus” to belittle valid challenge is the enemy of the scientific method.     

The nature,origin and modification of insoluble organic matter in chondrites,the major source of Earth’s C and N

All chondrites accreted ∼3.5 wt.% C in their matrices, the bulk of which was in a macromolecular solvent and acid insoluble organic material (IOM). Similar material to IOM is found in interplanetary dust particles (IDPs) and comets. The IOM accounts for almost all of the C and N in chondrites, and a significant fraction of the H. Chondrites and, to a lesser extent, comets were probably the major sources of volatiles for the Earth and the other terrestrial planets. Hence, IOM was both the major source of Earth’s volatiles and a potential source of complex prebiotic molecules.Large enrichments in D and ¹⁵N, relative to the bulk solar isotopic compositions, suggest that IOM or its precursors formed in very cold, radiation-rich environments. Whether these environments were in the interstellar medium (ISM) or the outer Solar System is unresolved. Nevertheless, the elemental and isotopic compositions and functional group chemistry of IOM provide important clues to the origin(s) of organic matter in protoplanetary disks. IOM is modified relatively easily by thermal and aqueous processes, so that it can also be used to constrain the conditions in the solar nebula prior to chondrite accretion and the conditions in the chondrite parent bodies after accretion.Here we review what is known about the abundances, compositions and physical nature of IOM in the most primitive chondrites. We also discuss how the IOM has been modified by thermal metamorphism and aqueous alteration in the chondrite parent bodies, and how these changes may be used both as petrologic indicators of the intensity of parent body processing and as tools for classification. Finally, we critically assess the various proposed mechanisms for the formation of IOM in the ISM or Solar System.     

Exchange rate of mercury between atmosphere and different kinds of Earth’s surfaces on the east slope of Mt. Gongga

Seasonal mercury flux measurements (three sampling periods: 2005 autumn, 2005 winter, and 2006 spring) were conducted at two soil sites, two grasslands and a cornfield on the east slope of Mt. Gongga, located in the southwest of Sichuan Province, China. T…