No Arabic abstract
In a first approximation, the Earths interior has an isotropic structure with a spherical symmetry. Over the last decades the geophysical observations have revealed, at different spatial scales, the existence of several perturbations from this basic structure. In this paper we discuss the hemispheric perturbations induced to this basic structure if the inner core is displaced from the center of mass of the Earth. Using numerical simulations of the observed hemispheric asymmetry of the seismic waves traveling through the upper inner core, with faster arrival times and higher attenuation in the Eastern Hemisphere, we estimate that the present position of the inner core is shifted by tens of kilometers from the Earths center eastward in the equatorial plane. If the only forces acting on the inner core were the gravitational forces, then its equilibrium position would be at the Earths center and the estimated displacement would not be possible. We conjecture that, due to interactions with the flow and the magnetic field inside the outer core, the inner core is in a permanent chaotic motion. To support this hypothesis we analyze more than ten different geophysical phenomena consistent with an inner core motion dominated by time scales from hundreds to thousands of years.
It has long been assumed the Earths solid inner core started to grow when molten iron cooled to its melting point. However, the nucleation mechanism, which is a necessary step of crystallization, has not been well understood. Recent studies found it requires an unrealistic degree of undercooling to nucleate the stable hexagonal close-packed (hcp) phase of iron, which can never be reached under the actual Earths core conditions. This contradiction leads to the inner core nucleation paradox [1]. Here, using a persistent-embryo method and molecular dynamics simulations, we demonstrate that the metastable body-centered cubic (bcc) phase of iron has a much higher nucleation rate than the hcp phase under inner-core conditions. Thus, the bcc nucleation is likely to be the first step of inner core formation instead of direct nucleation of the hcp phase. This mechanism reduces the required undercooling of iron nucleation, which provides a key factor to solve the inner-core nucleation paradox. The two-step nucleation scenario of the inner core also opens a new avenue for understanding the structure and anisotropy of the present inner core.
In a first approximation the Earths interior has an isotropic structure with a spherical symmetry. Over the last decades the geophysical observations have revealed, at different spatial scales, the existence of several perturbations from this basic structure. Some of them are situated in the neighborhood of the inner core boundary (ICB). One of the best documented perturbations is the asymmetry at the top of the inner core (ATIC) characterized by faster seismic wave velocity in the eastern hemisphere than in the western hemisphere. All existing explanations are based on a hemispheric variation of the material properties near ICB inside the inner core. Using numerical simulations of the seismic ray propagation, we show that the ATIC can be explained as well by the displacement of the inner core towards east in the equatorial plane tens of kilometers from the Earths center, without modifying the spherical symmetry in the upper inner core. The hypothesis of a displaced inner core is also sustained by other observed hemispheric asymmetries at the top of the inner core and at the bottom of the outer core. A displaced inner core would have major implications for many mechanical, thermal, and magnetic phenomena in the Earths interior.
Ensembles of geophysical models improve projection accuracy and express uncertainties. We develop a novel data-driven ensembling strategy for combining geophysical models using Bayesian Neural Networks, which infers spatiotemporally varying model weights and bias while accounting for heteroscedastic uncertainties in the observations. This produces more accurate and uncertainty-aware projections without sacrificing interpretability. Applied to the prediction of total column ozone from an ensemble of 15 chemistry-climate models, we find that the Bayesian neural network ensemble (BayNNE) outperforms existing ensembling methods, achieving a 49.4% reduction in RMSE for temporal extrapolation, and a 67.4% reduction in RMSE for polar data voids, compared to a weighted mean. Uncertainty is also well-characterized, with 90.6% of the data points in our extrapolation validation dataset lying within 2 standard deviations and 98.5% within 3 standard deviations.
The geoid is the true physical figure of the Earth, a particular equipotential surface of the gravity field of the Earth that accounts for the effect of all subsurface density variations. Its shape approximates best (in the sense of least squares) the mean level of oceans, but the geoid is more difficult to determine over continents. Satellite missions carry out distance measurements and derive the gravity field to provide geoid maps over the entire globe. However, they require calibration and extensive computations including integration, which is a non-unique operation. Here we propose a direct method and a new tool that directly measures geopotential differences on continents using atomic clocks. General Relativity Theory predicts constant clock rate at sea level, and faster (resp. slower) clock rate above (resp. below) sea level. The technology of atomic clocks is on the doorstep of reaching an accuracy level in clock rate that is equivalent to 1 cm in determining equipotential surface (including geoid) height. We discuss the value and future applicability of such measurements including direct geoid mapping on continents, and joint gravity and geopotential surveying to invert for subsurface density anomalies. Our synthetic calculations show that the geoid perturbation caused by a 1.5 km radius sphere with 20% density anomaly buried at 2 km depth in the crust of the Earth is already detectable by atomic clocks of achievable accuracy. Therefore atomic clock geopotential surveys, used together with relative gravity data to benefit from their different depth sensitivities, can become a useful tool in mapping density anomalies within the Earth.
Regional characterization of the continental crust has classically been performed through either geologic mapping, geochemical sampling, or geophysical surveys. Rarely are these techniques fully integrated, due to limits of data coverage, quality, and/or incompatible datasets. We combine geologic observations, geochemical sampling, and geophysical surveys to create a coherent 3-D geologic model of a 50 x 50 km upper crustal region surrounding the SNOLAB underground physics laboratory in Canada, which includes the Southern Province, the Superior Province, the Sudbury Structure and the Grenville Front Tectonic Zone. Nine representative aggregate units of exposed lithologies are geologically characterized, geophysically constrained, and probed with 109 rock samples supported by compiled geochemical databases. A detailed study of the lognormal distributions of U and Th abundances and of their correlation permits a bivariate analysis for a robust treatment of the uncertainties. A downloadable 3D numerical model of U and Th distribution defines an average heat production of 1.5$^{+1.4}_{-0.7}$$mu$W/m$^{3}$, and predicts a contribution of 7.7$^{+7.7}_{-3.0}$TNU (a Terrestrial Neutrino Unit is one geoneutrino event per 10$^{32}$ target protons per year) out of a crustal geoneutrino signal of 31.1$^{+8.0}_{-4.5}$TNU. The relatively high local crust geoneutrino signal together with its large variability strongly restrict the SNO+ capability of experimentally discriminating among BSE compositional models of the mantle. Future work to constrain the crustal heat production and the geoneutrino signal at SNO+ will be inefficient without more detailed geophysical characterization of the 3D structure of the heterogeneous Huronian Supergroup, which contributes the largest uncertainty to the calculation.