|Home | About | Journals | Submit | Contact Us | Français|
In this paper, we use a robust inversion algorithm, which we have tested in many regional studies, to obtain the first global model of Curie-point depth (GCDM) from magnetic anomaly inversion based on fractal magnetization. Statistically, the oceanic Curie depth mean is smaller than the continental one, but continental Curie depths are almost bimodal, showing shallow Curie points in some old cratons. Oceanic Curie depths show modifications by hydrothermal circulations in young oceanic lithosphere and thermal perturbations in old oceanic lithosphere. Oceanic Curie depths also show strong dependence on the spreading rate along active spreading centers. Curie depths and heat flow are correlated, following optimal theoretical curves of average thermal conductivities K=~2.0W(m°C)−1 for the ocean and K=~2.5W(m°C)−1 for the continent. The calculated heat flow from Curie depths and large-interval gridding of measured heat flow all indicate that the global heat flow average is about 70.0mW/m2, leading to a global heat loss ranging from ~34.6 to 36.6 TW.
Our planet Earth is under constant cooling and differentiation since its origin. Its current thermal state and heat loss can be estimated from surface heat flow measurements1,2,3,4,5,6. However, heat flow estimates have evident drawbacks; they are often measured at sparsely and irregularly distributed sites, and they are strongly affected by shallow hydrothermal circulation and therefore are limited in inferring deep thermal structure of the lithosphere. These uncertainties put quite different estimates of global heat loss. In one early study1, and in some recent calibrations minimizing the hydrothermal effects3,4, the global heat loss was estimated at about 40–47 TW (or 4.4–4.7×1013W). Other reappraisals of global heat flow database, however, concluded that global conductive heat loss falls in the range from 29 to 34 TW2,5,6. The differences between these estimates vary considerably.
An alternate and independent method of studying deep thermal structure and global heat loss is by detecting Curie depths from inversion of surface total field magnetic anomalies7,8,9,10,11,12,13,14. Global coverage of magnetic anomalies is being constantly improved in recent years15,16, and this allows it possible to map the global Curie isotherm in high-resolution, assuming that the Curie-point temperature is constant globally around 550°C17,18 and that lateral compositional variations exert neglectable influence on the Curie temperature.
There are several different technical schemes of Curie depth inversion from magnetic anomalies, and even for the same method the selected inversion parameters can also vary among scientists7,8,9,10,11,12,13,14. Furthermore, previous studies focus only on a particular region of the Earth. These facts render it impossible the global comparison of Curie depth results.
In early regional studies, we have successfully automated the centroid method based on radially-averaged amplitude spectrum to obtain high-resolution Curie depths9,10,11. This saves considerably the computational time. Now we can naturally apply this algorithm to the Earth Magnetic Anomaly Grid of 2-arc-minute resolution (EMAG2, http://geomag.org/)16 to obtain the first global reference model of Curie-point depths (GCDM).
The Earth’s surface from −75° to 75° latitudes is divided into 72 zones, and the two polar zones are treated separately. Within each zone, we transform the magnetic data from geographical coordinate to Cartesian coordinate, before gridding the data in a constant 2.6km interval using the minimum curvature method, which iteratively solves a set of differential equations to minimize the total second horizontal derivative and honor input data19. We then estimate Curie depths using three different window sizes, 98.8×98.8km2, 195.0×195.0km2, and 296.4×296.4km2. The size of 296.4×296.4km2 is large enough to capture the deepest Curie depths around 50km, and is larger than in most existing applications of the algorithm. The moving steps for these three window sizes are 49.4km, 97.5km, and 98.8km, respectively. As did before20, we take the average of Curie depths from the three windows as the final Curie depth model. Taken to be constant at 3, the scaling factor of 3D fractal magnetic sources, is defined by, in which is the 3D power spectrum of the magnetization, kx, ky, and kz are wavenumbers in x, y, and z directions, respectively, and their Euclidean norm
The datum altitude (4km) of the EMAG2 is subtracted from the estimated Curie depths to be in reference to the geoid (Fig. 1).
We applied the standard fast Fourier transform scheme in estimating spectra, but with careful consideration and treatment of boundary effects and short-wavelength tails21. We use the wavenumber band of ~0.005–0.03 1/km for estimating the centroid depths, not just the first slope for stability reasons. It is true that same spectral wavenumber range is not applicable for every window. Our algorithm can judge automatically if the first (smallest wavenumber) point is smaller than the second in the spectra and, if so, does not account this point in the calculation, because theoretically this cannot happen11. In each window, we select the best depth in the least-squares sense.
Our efforts in deriving the first global reference model of Curie depth from magnetic anomalies give us a total number of 516,772 raw Curie depth estimates using the three different window sizes. Different window sizes emphasize different scales of Curie depth variations, but their average, as the final reference model, reduces random noises and improves the resolution model data is available in the supplementary Information. The final spatial resolution of the Curie depth model should be smaller than the smallest moving step, which is 49.4km. The gridding resolution of our Curie depth model is 10minutes.
Our global Curie depth model (GCDM) shows evident shallowing along mid-ocean ridges, particularly in the Indian, north Atlantic, and Arctic Oceans (Fig. 1). This reflects the active magmatism and upwelling of the hot asthenosphere. GCDM confirms the early results in the north Atlantic11, that in the background of increasing Curie depths with oceanic crustal ages22, anomalous off-ridge upwelling and oscillations in Curie isotherm are observed in all major sea basins, and in most cases, they are associated with small-scale convection or hotspots. For example, the Hawaii seamount trail, the Bermuda Rise, and the Cape Verde Islands all show small Curie depths.
In both the Pacific and Atlantic, large-scale upwells in Curie points are observed in the age range from ~100 to 150Ma (Fig. 1). By plotting all Curie depth estimates against oceanic crustal ages, we found that this large-scale Curie depth anomaly persists globally (Fig. 2). Assuming an average lithospheric thermal diffusivity κ of 0.319mm2/s (or 10.06km2/Ma), observed average Curie depths are larger than predicted from the half-space cooling model for ages <~30Ma, but are smaller for older crustal ages. This implies that, near the mid-ocean ridges, temperatures are lower than that from the model prediction, probably due to stronger hydrothermal activities, but deep thermal perturbations exist in older lithospheres, particularly in the age range of ~100–150Ma.
We find that Curie depths within the two 1Ma isochrons of active ridges are strongly dependent on the spreading rate (Rr) of Müller et al.22. Ultra-slow spreading ridges (Rr<4mm/yr) show large Curie depths averaged at ~17.5km (Fig. 3). However, the average Curie depth decreases abruptly at slightly higher rates, and then with increasing spreading rates from about 4 to 45mm/yr (slow to medium spreading rate), the average Curie depth increases linearly, from ~12 to ~18km. Along fast spreading ridges, average Curie depths decrease very slowly with rates and keep almost constant at ~17.7km. Likewise along super-fast ridges, average Curie depths also decrease very slowly with spreading rates, but keep at a slightly smaller constant at ~15.2km and with a narrower distribution. If the fractal scaling factors of magnetization are not dependent on the spreading rate and the Curie temperatures are constant, this overall complex dependence of Curie depths on the spreading rate indicate that the near-ridge thermal structures vary considerably with spreading rate.
The following hypothesis may be perceived to interpret these observations:
The majority of the largest Curie depths are found on the continents, most noticeably South America, Africa, and Indian Shield (Fig. 1). However, small Curie depths are also found within some old and stable cratons, such as the Northern China, Siberia, and North America, where large Curie depths are expected from their old ages. This may reflect that reactivation of these old continental cratons has induced thermal perturbations at depth. North China craton, for example, may have been reactivated by the subduction of the Pacific Plate23. The zone of small Curie depths along the western margin of America is also caused by active subductions and related arc magmatism and extension. In southern Africa, a large patch of shallow Curie points could be induced by hot spot magmatism, so is the noticeable Curie depth anomaly associated with the large intraplate hot spot Ahaggar Swell in northern Africa.
Oceanic Curie depths show a good normal distribution with a smaller mean of 20.18km and a smaller standard deviation than continental Curie depths (Fig. 4). By fitting continental Curie depth distribution with a mixture of two normal distributions, we find that the distribution is nearly bimodal (Fig. 4), which shows a first peak depth of 22.43km, but a minor second peak occurs around 32.76km, giving a tail of large Curie depths.
Surface heat flow measures the heat loss at the surface of the earth; it could be modulated by hydrothermal circulation or sedimentary cover, but if evenly and densely distributed, its global average quantifies the global heat loss. Large-interval gridding also reduces local measurement biases in heat flow11,20, and hence increases the reliability of interpretation. Our heat flow model, based on a constant 1° interval gridding of all the raw data from the International Heat Flow Commission Database (http://www.heatflow.und.edu/; last updated in January 2011), captures the characteristics of main tectonic features (Fig. 5). Averaging of heat flow values is based on block averaging of (x, y, z) data by L2 norm. We estimate a mean position and heat flow value for every non-empty block in each 1°×1° grid region to avoid aliasing short wavelengths. We then apply a constant 1° interval gridding of the averaged heat flow using the minimum curvature algorithm with a tension factor of 0.5. Histograms show that gridded heat flow reduces the number of extremely low and high heat flow values caused by hydrothermal process, and retains the original statistical mean (Fig. 6a,b). At the 5% significance level, the Kolmogorov Smirnov statistical test rejected the null hypothesis that the raw and gridded heat flow are from the same continuous distribution (Fig. 6c).
If we interpolate all the raw heat flow data in a constant 100km interval, using the minimum curvature algorithm with a tension factor of 0.5, the global heat flow average is 67.8mW/m2. Alternatively, we interpolate the same data set in a constant 1° interval (Fig. 5), and then resample the grid in a constant 10minute interval. Then in each 1°×1° quadrangle we estimate the average of the gridded heat flow, and the global mean from these quadrangles is 68.3mW/m2. These values are substantially smaller than some of the recent estimates (e.g., 91.6mW/m2)24.
Previous regional studies have shown that heat flow and Curie depth are correlated, following a theoretical thermal conduction relationship10,11. By comparing the Curie depth map (Fig. 1) with heat flow map (Fig. 5), we observe that this correlation should hold for the global data, particularly in active tectonic areas such as mid-ocean ridges and subduction zones. In old and stable continental cratons, surface heat flow may deviate from deep thermal structures due to their long evolutional histories and late-stage thermal perturbations. In thick continental lithosphere, there may be long delays between changes in asthenosphere temperature and measurable effects arriving at the surface25,26. Shallow radioactive contribution to surface heat flow is also larger in the continental domains.
Figure 7 Shows a large scattering in the global correlation between heat flow and Curie depth; this is expected since there are diverse geological units and local structures. Nevertheless, good correlations between heat flow and Curie depth can be observed, i.e., high heat flow measurements tend to correlate with small Curie depths, and vice versa. It is seen that most continental points are clustered in the vicinity of the continental theoretical curve with an average thermal conductivity K of ~2.5W(m°C)−1, and most oceanic data can be best fitted with the oceanic theoretical curve of an average K=~2.0W(m°C)−1. These conductivities are compatible with those of granite and basalt27,28.
The Earth is divided into 360×180 bins of 1°×1° in size. The exact surface area (Ai) of these latitude-longitude quadrangle bins are calculated assuming a spherical ellipsoid. Calculated heat flow data from Curie depths are gridded in a 10minute interval. The surface heat flow measurements are interpolated using a constant 1° interval and the grid is also resampled in a 10minute interval. These data are then averaged, respectively, within each 1°×1° latitude-longitude quadrangle (Qsi). Quadrangles without heat flow measurements or Curie depth constraints are assumed to take the global average of heat flow. The total heat loss (Thl) is then calculated by
We find that the measured global heat loss is 3.46×1013W (or 34.6 TW) based on our interpolated heat flow. Calculated conductive heat flow from Curie depths using an average continental thermal conductivity k=2.5W(m°C)−1 and an oceanic average k=2.0W(m°C)−1, is averaged at 72.125mW/m2. Here the radioactive heat contribution is included in the estimation. The global heat loss hereby obtained is 3.66×1013W (or 36.6 TW). These values are much lower than some of the early estimates3, but slightly larger than some others2,5,6.
Our global Curie depth model (GCDM) shows that oceanic Curie depths increase with crustal ages, but deviate from the half-space cooling model in both young and old oceanic lithospheres. Measured Curie depths are larger than the model predictions in young oceanic lithosphere of age <~30Ma, but are smaller in older lithospheres. In particular, there is a significant decrease of Curie depths for crustal ages between ~100 and 150Ma, noticeable in both Atlantic and Pacific. These deviations reflect strong hydrothermal circulations in young oceanic lithospheres and thermal perturbations in aged oceanic lithospheres.
Along active mid-ocean ridges, our study shows that Curie depths are strongly dependent on the spreading rate. There is a unique Curie depth distribution for each single spreading-rate type (i.e., ultra-slow, slow, medium, fast, and super-fast ridges), indicating strong coupling between spreading rate and seafloor-spreading mechanisms, thermal structures, and hydrothermal activities.
Continental domains have the largest Curie depths, but also show nearly a bimodal distribution; small Curie depths are found in some of old continental cratons, indicating that deep thermal reactivations caused inconsistencies between Curie depth and surface heat flow.
Despite a large scatter of points, global correlations between Curie depth and gridded surface heat flow in 1° interval follow theoretical curves of an average thermal conductivity of K=2.0W(m°C)−1 for the oceanic data and an continental average K=2.5W(m°C)−1. With these optimal conductivities, the calculated global heat flow average is 72.125mW/m2. This estimate is close to measured global heat flow average from direct large-interval gridding of raw heat flow data, which is about 68.0mW/m2. We propose that the global heat flow average is around 70.0mW/m2 and the global heat loss ranges from ~34.6 to 36.6 TW.
How to cite this article: Li, C.-F. et al. A global reference model of Curie-point depths based on EMAG2. Sci. Rep. 7, 45129; doi: 10.1038/srep45129 (2017).
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This research is funded by National Science Foundation of China (Grant No. 91428309 and 91028007). We thank Dhananjay Ravat and Yaoling Niu for discussions. We also thank two anonymous reviewers who provided in-depth reviews and suggestions for improving this paper. Data processing and mapping are supported by the USGS potential field software, and by GMT.
The authors declare no competing financial interests.
Author Contributions C.-F. Li improved the algorithm, wrote the software, processed and interpreted the data, prepared figures, and wrote the main manuscript text. Y. Lu processed the data in the two polar regions. J. Wang helped in the data processing.