The Accuracy Assessment of the PREM and AK135-F Radial Density Models

The Earth’s synthetic density and gravitational models can be used to validate numerical methods for global (or large-scale) gravimetric forward and inverse modelling formulated either in the spatial or spectral domains. The Preliminary Reference Earth Model (PREM) density parameters can be adopted as a 1-D reference density model and further refined using more detailed 2-D or 3-D crust and mantle density models. Alternatively, the AK135-F density parameters can be used for this purpose. In this study, we investigate options for a refinement of the Earth’s synthetic density model by assessing the accuracy of available 1-D density models, specifically the PREM and AK135-F radial density parameters. First, we use density parameters from both models to estimate the Earth’s total mass and compare these estimates with published results. We then estimate the Earth’s gravity field parameters, particularly the geoidal geopotential number W₀ and the mean gravitational attraction and compare them with published values. According to our results, the Earth’s total mass from the two models (the PREM and the AK135-F) differ less than 0.02% and 0.01%, respectively, when compared with the value adopted by the International Astronomical Union (IAU). The geoidal geopotential values of the two models differ from the value adopted by the IAU by less than 0.1% and 0.04%, respectively. The values of the mean gravitational attraction of the two models differ less than 0.02% and 0.08%, respectively, when compared with the value obtained from the geocentric gravitational constant and the Earth’s mean radius. These numerical findings ascertain that the PREM and AK135-F density parameters are suitable for defining a 1-D reference density model.

Lateral Variation of Crustal Lg Attenuation in Eastern North America

We perform Q _Lg tomography for the northeastern part of North America. Vertical broadband seismograms of 473 crustal earthquakes recorded by 302 stations are processed to extract the Lg amplitude spectra. Tomographic inversions are independently conducted at 58 discrete frequencies distributed evenly in log space between 0.1 and 20.0 Hz. This relatively large dataset with good ray coverage allows us to image lateral variation of the crustal attenuation over the region. Obtained Q _Lg maps at broadband and individual frequencies provide new insights into the crustal attenuation of the region and its relationship to geological structures and past tectonic activity in the area. The Q _Lg shows more uniform values over the older, colder, and drier Canadian Shield, in contrast to higher variations in the younger margins. Results confirm the correlation of large-scale variations with crustal geological features in the area. Existence of low-velocity anomalies, thick sediments, volcanic rocks, and thin oceanic crust are potential sources of observed anomalies. The mean Q values are inversely correlated with average heat flow/generation for main geological provinces.

Electromagnetic detection of a 410-km-deep melt layer in the southwestern United States

A deep-seated melt or fluid layer on top of the 410-km-deep seismic discontinuity in Earth's upper mantle, as proposed in the transition-zone 'water filter' hypothesis, may have significant bearing on mantle dynamics and chemical differentiation. The geophysical detection of such a layer has, however, proved difficult. Magnetotelluric and geomagnetic depth sounding are geophysical methods sensitive to mantle melt. Here we use these methods to search for a distinct structure near 410-km depth. We calculate one-dimensional forward models of the response of electrical conductivity depth profiles, based on mineral physics studies of the effect of incorporating hydrogen in upper-mantle and transition-zone minerals. These models indicate that a melt layer at 410-km depth is consistent with regional magnetotelluric and geomagnetic depth sounding data from the southwestern United States (Tucson). The 410-km-deep melt layer in this model has a conductance of 3.0 x 10(4) S and an estimated thickness of 5-30 km. This is the only regional data set that we have examined for which such a melt layer structure was found, consistent with regional seismic studies. We infer that the hypothesized transition-zone water filter occurs regionally, but that such a layer is unlikely to be a global feature.

Seismic anisotropy: tracing plate dynamics in the mantle

Elastic anisotropy is present where the speed of a seismic wave depends on its direction. In Earth's mantle, elastic anisotropy is induced by minerals that are preferentially oriented in a directional flow or deformation. Earthquakes generate two seismic wave types: compressional (P) and shear (S) waves, whose coupling in anisotropic rocks leads to scattering, birefringence, and waves with hybrid polarizations. This varied behavior is helping geophysicists explore rock textures within Earth's mantle and crust, map present-day upper-mantle convection, and study the formation of lithospheric plates and the accretion of continents in Earth history.

The mantle flow field beneath western North America

Although motions at the surface of tectonic plates are well determined, the accompanying horizontal mantle flow is not. We have combined observations of surface deformation and upper mantle seismic anisotropy to estimate this flow field for western North America. We find that the mantle velocity is 5.5 +/- 1.5 centimeters per year due east in a hot spot reference frame, nearly opposite to the direction of North American plate motion (west-southwest). The flow is only weakly coupled to the motion of the surface plate, producing a small drag force. This flow field is probably due to heterogeneity in mantle density associated with the former Farallon oceanic plate beneath North America.

Formation and evolution of Archaean cratons: insights from southern Africa

Archaean cratons are the stable remnants of Earth’s early continental lithosphere, and their structure, composition and survival over geological time make them unique features of the Earth’s surface. The Kaapvaal Project of southern Africa was organized around a broadly diverse scientific collaboration to investigate fundamental questions of craton formation and mantle differentiation in the early Earth. The principal aim of the project was to characterize the physical and chemical nature of the crust and mantle of the cratons of southern Africa in geological detail, and to use the 3D seismic and geochemical images of crustal and mantle heterogeneity to reconstruct the assembly history of the cratons. Seismic results confirm that the structure of crust and tectospheric mantle of the cratons differs significantly from that of post-Archaean terranes. Three-dimensional body-wave tomographic images reveal that high-velocity mantle roots extend to depths of at least 200 km, and locally to depths of 250–300 km beneath cratonic terranes. No low-velocity channel has been identified beneath the cratonic root. The Kaapvaal Craton was modified approximately 2.05 Ga by the Bushveld magmatic event, and the mantle beneath the Bushveld Province is characterized by relatively low seismic velocities. The crust beneath undisturbed Archaean craton is relatively thin (c. 35–40 km), unlayered and characterized by a strong velocity contrast across a sharp Moho, whereas post-Archaean terranes and Archaean regions disrupted by large-scale Proterozoic magmatic or tectonic events are characterized by thicker crust, complex Moho structure and higher seismic velocities in the lower crust. A review of Re-Os depletion model age determinations confirms that the mantle root beneath the cratons is Archaean in age. The data show also that there is no apparent age progression with depth in the mantle keel, indicating that its thickness has not increased over geological time. Both laboratory experiments and geochemical results from eclogite xenoliths suggest that subduction processes played a central role in the formation of Archaean crust, the melt depletion of Archaean mantle and the assembly of early continental lithosphere. Co-ordinated geochronological studies of crustal and mantle xenoliths have revealed that both crust and mantle have experienced a multi-stage history. The lower crust in particular retains a comprehensive record of the tectonothermal evolution of the lithosphere. Analysis of lower-crustal xenoliths has shown that much of the deep craton experienced a dynamic and proteracted history of tectonothermal activity that is temporally associated with events seen in the surface record. Cratonization thus occurred not as a discrete event, but in stages, with final stabilization postdating crustal formation.

Deflection of mantle plume material by cratonic keels

Lithosphere that formed in Archaean and possibly early Proterozoic time is thicker, more buoyant, and geochemically distinct from lithosphere that formed after about 2.3 Ga. Mantle xenolith and seismic data indicate that some cratonic roots, or ‘keels’, extend to depths of c. 250 km, compared with normal continental lithosphere of thickness 150 km or less; yet many cratons have experienced uplift, dyking and kimberlite emplacement, suggesting interactions with hot, rising asthenosphere referred to as mantle plumes. Plumes supply additional heat to the base of the lithospheric plates, whose base can be heated and entrained in the flow (thermal erosion). How have these cratonic keels persisted despite their interactions with mantle plumes? The geometry of cratonic keels during their interactions with mantle plumes is a critical factor controlling keel preservation. To a laterally spreading plume head, cratonic keels appear as major obstacles, and the hot, buoyant plume material ponds beneath thinner lithosphere. Our model simulations show that deep keels deflect mantle plume material and that steep gradients at the lithosphere-asthenosphere boundary between Archaean keels and ‘normal’ lithosphere will focus flow, leading to localized adiabatic decompression melting. Plume processes can lead to a reduction in the breadth of a cratonic root where the plume rises beneath the craton, regardless of the initial breadth of the craton. Where the plume rises beneath a craton the hot plume material will spread laterally beneath the keel and attain thicknesses of tens of kilometres. This transfers heat to the base of the lithosphere and could generate small volumes of melt at considerable depth, depending on the composition of the lower lithosphere. We have used model simulations of plumes beneath Africa to predict the magnitude and direction of seismic anisotropy caused by lateral flow of hot plume material beneath and around a cratonic keel. The shear-wave splitting in our models is greatest at the edge of the cratonic keel, and its azimuth is parallel to the plume flow direction.

African hot spot volcanism: small-scale convection in the upper mantle beneath cratons

Numerical models demonstrate that small-scale convection develops in the upper mantle beneath the transition of thick cratonic lithosphere and thin oceanic lithosphere. These models explain the location and geochemical characteristics of intraplate volcanos on the African and South American plates. They also explain the presence of relatively high seismic shear wave velocities (cold downwellings) in the mantle transition zone beneath the western margin of African cratons and the eastern margin of South American cratons. Small-scale, edge-driven convection is an alternative to plumes for explaining intraplate African and South American hot spot volcanism, and small-scale convection is consistent with mantle downwellings beneath the African and South American lithosphere.

Multiple seismic discontinuities near the base of the transition zone in the Earth's mantle

The seismologically defined boundary between the transition zone in the Earth's mantle (410-660 km depth) and the underlying lower mantle is generally interpreted to result from the breakdown of the gamma-spinel phase of olivine to magnesium-perovskite and magnesiowustite. Laboratory measurements of these transformations of olivine have determined that the phase boundary has a negative Clapeyron slope and does indeed occur near pressures corresponding to the base of the transition zone. But a computational study has indicated that, because of the presence of garnet minerals, multiple seismic discontinuities might exist near a depth of 660 km (ref. 4), which would alter the simple negative correlation of changes in temperature with changes in the depth of the phase boundary. In particular, garnet minerals undergo exothermic transformations near this depth, acting to complicate the phase relations and possibly effecting mantle convection processes in some regions. Here we present seismic evidence that supports the existence of such multiple transitions near a depth of 660 km beneath southern California. The observations are consistent with having been generated by garnet transformations coupling with the dissociation of the gamma-spinel phase of olivine. Temperature anomalies calculated from the imaged discontinuity depths--using Clapeyron slopes determined for the various transformations--generally match those predicted from an independent P-wave velocity model of the region.