Mantle Transition Zone Discontinuities beneath the Contiguous United States

Using over 310,000 high-quality radial receiver functions recorded by the USArray and other seismic stations in the contiguous United States, the depths of the 410 km and 660 km discontinuities (d410 and d660) are mapped in over 1,000 consecutive overlapping circles with a radius of 1⁰. The average mantle transition zone (MTZ) thickness for both the western and central/eastern U.S. is within 3 km from the global average of 250 km, suggesting an overall normal MTZ temperature beneath both areas. The Pacific Coast Ranges and the southern Basin and Range Province are underlain by a depressed d410, indicating higher-than-normal temperature in the upper MTZ. The proposed Yellowstone and Raton hot spots are not associated with clear undulations of the MTZ discontinuities, but d410 beneath another proposed hot spot, Bermuda, is depressed significantly and d660 has a normal depth. Low-temperature regions are found in the upper MTZ associated with the subducted Juan de Fuca slab beneath the northern Rocky Mountains and in two circular areas beneath the northern Basin and Range Province and the southern Colorado Plateau. Part of the Great Plains is characterized by a depressed d660. This observation, when combined with results from seismic tomography, suggests the existence of a cold region in the lower MTZ, probably associated with subducted Farallon slab segments

First principles simulations of liquid Fe-S under Earth's core conditions

First principles electronic structure calculations, based upon density functional theory within the generalized gradient approximation and ultra-soft Vanderbilt pseudopotentials, have been used to simulate a liquid alloy of iron and sulfur at Earth's core conditions. We have used a sulfur concentration of $\approx 12$wt, in line with the maximum recent estimates of the sulfur abundance in the Earth's outer core. The analysis of the structural, dynamical and electronic structure properties has been used to report on the effect of the sulfur impurities on the behavior of the liquid. Although pure sulfur is known to form chains in the liquid phase, we have not found any tendency towards polymerization in our liquid simulation. Rather, a net S-S repulsion is evident, and we propose an explanation for this effect in terms of the electronic structure. The inspection of the dynamical properties of the system suggests that the sulfur impurities have a negligible effect on the viscosity of Earth's liquid core.Comment: 24 pages (including 8 figures

Lithospheric Structure and Evolution of Southern Africa: Constraints from Joint Inversion of Rayleigh Wave Dispersion and Receiver Functions

We conduct a joint inversion of teleseismic receiver functions and Rayleigh wave phase velocity dispersion from both ambient noise and earthquakes using data from 79 seismic stations in southern Africa, which is home to some of the world\u27s oldest cratons and orogenic belts. The area has experienced two of the largest igneous activities in the world (the Okavango dyke swarm and Bushveld mafic intrusion) and thus is an ideal locale for investigating continental formation and evolution. The resulting 3-D shear wave velocities for the depth range of 0—100Â km and crustal thickness measurements show a clear spatial correspondence with known geological features observed on the surface. Higher than normal mantle velocities found beneath the southern part of the Kaapvaal craton are consistent with the basalt removal model for the formation of cratonic lithosphere. In contrast, the Bushveld complex situated within the northern part of the craton is characterized by a thicker crust and higher crustal Vp/Vs but lower mantle velocities, which are indicative of crustal underplating of mafic materials and lithospheric refertilization by the world\u27s largest layered mafic igneous intrusion. The thickened crust and relatively low elevation observed in the Limpopo belt, which is a late Archean collisional zone between the Kaapvaal and Zimbabwe cratons, can be explained by eclogitization of the basaltic lower crust. The study also finds evidence for the presence of a stalled segment of oceanic lithosphere beneath the southern margin of the Proterozoic Namaqua-Natal mobile belt

Structural and elastic anisotropy of crystals at high pressures and temperatures from quantum mechanical methods: The case of Mg2SiO4 forsterite

We report accurate ab initio theoretical predictions of the elastic, seismic, and structural anisotropy of the orthorhombic Mg2SiO4 forsterite crystal at high pressures (up to 20 GPa) and temperatures (up to its melting point, 2163 K), which constitute earth’s upper mantle conditions. Single-crystal elastic stiffness constants are evaluated up to 20 GPa and their first- and second-order pressure derivatives reported. Christoffel’s equation is solved at several pressures: directional seismic wave velocities and related properties (azimuthal and polarization seismic anisotropies) discussed. Thermal structural and average elastic properties, as computed within the quasi-harmonic approximation of the lattice potential, are predicted at high pressures and temperatures: directional thermal expansion coefficients, first- and second-order pressure derivatives of the isothermal bulk modulus, and P-V-T equation-of-state. The effect on computed properties of five different functionals, belonging to three different classes of approximations, of the density functional theory is explicitly investigated

Mid-mantle deformation inferred from seismic anisotropy

With time, convective processes in the Earth's mantle will tend to align crystals, grains and inclusions. This mantle fabric is detectable seismologically, as it produces an anisotropy in material properties—in particular, a directional dependence in seismic-wave velocity. This alignment is enhanced at the boundaries of the mantle where there are rapid changes in the direction and magnitude of mantle flow, and therefore most observations of anisotropy are confined to the uppermost mantle or lithosphere and the lowermost-mantle analogue of the lithosphere, the D" region. Here we present evidence from shear-wave splitting measurements for mid-mantle anisotropy in the vicinity of the 660-km discontinuity, the boundary between the upper and lower mantle. Deep-focus earthquakes in the Tonga–Kermadec and New Hebrides subduction zones recorded at Australian seismograph stations record some of the largest values of shear-wave splitting hitherto reported. The results suggest that, at least locally, there may exist a mid-mantle boundary layer, which could indicate the impediment of flow between the upper and lower mantle in this region

Solar System Processes Underlying Planetary Formation, Geodynamics, and the Georeactor

Only three processes, operant during the formation of the Solar System, are responsible for the diversity of matter in the Solar System and are directly responsible for planetary internal-structures, including planetocentric nuclear fission reactors, and for dynamical processes, including and especially, geodynamics. These processes are: (i) Low-pressure, low-temperature condensation from solar matter in the remote reaches of the Solar System or in the interstellar medium; (ii) High-pressure, high-temperature condensation from solar matter associated with planetary-formation by raining out from the interiors of giant-gaseous protoplanets, and; (iii) Stripping of the primordial volatile components from the inner portion of the Solar System by super-intense solar wind associated with T-Tauri phase mass-ejections, presumably during the thermonuclear ignition of the Sun. As described herein, these processes lead logically, in a causally related manner, to a coherent vision of planetary formation with profound implications including, but not limited to, (a) Earth formation as a giant gaseous Jupiter-like planet with vast amounts of stored energy of protoplanetary compression in its rock-plus-alloy kernel; (b) Removal of approximately 300 Earth-masses of primordial gases from the Earth, which began Earth's decompression process, making available the stored energy of protoplanetary compression for driving geodynamic processes, which I have described by the new whole-Earth decompression dynamics and which is responsible for emplacing heat at the mantle-crust-interface at the base of the crust through the process I have described, called mantle decompression thermal-tsunami; and, (c)Uranium accumulations at the planetary centers capable of self-sustained nuclear fission chain reactions.Comment: Invited paper for the Special Issue of Earth, Moon and Planets entitled Neutrino Geophysics Added final corrections for publicatio

Iron under Earth's core conditions: Liquid-state thermodynamics and high-pressure melting curve

{\em Ab initio} techniques based on density functional theory in the projector-augmented-wave implementation are used to calculate the free energy and a range of other thermodynamic properties of liquid iron at high pressures and temperatures relevant to the Earth's core. The {\em ab initio} free energy is obtained by using thermodynamic integration to calculate the change of free energy on going from a simple reference system to the {\em ab initio} system, with thermal averages computed by {\em ab initio} molecular dynamics simulation. The reference system consists of the inverse-power pair-potential model used in previous work. The liquid-state free energy is combined with the free energy of hexagonal close packed Fe calculated earlier using identical {\em ab initio} techniques to obtain the melting curve and volume and entropy of melting. Comparisons of the calculated melting properties with experimental measurement and with other recent {\em ab initio} predictions are presented. Experiment-theory comparisons are also presented for the pressures at which the solid and liquid Hugoniot curves cross the melting line, and the sound speed and Gr\"{u}neisen parameter along the Hugoniot. Additional comparisons are made with a commonly used equation of state for high-pressure/high-temperature Fe based on experimental data.Comment: 16 pages including 6 figures and 5 table

Off-Diagonal Elements of the DeWitt Expansion from the Quantum Mechanical Path Integral

The DeWitt expansion of the matrix element M_{xy} = \left\langle x \right| \exp -[\case{1}{2} (p-A)^2 + V]t \left| y \right\rangle, $(p=-i\partial)$ in powers of $t$ can be made in a number of ways. For $x=y$ (the case of interest when doing one-loop calculations) numerous approaches have been employed to determine this expansion to very high order; when $x \neq y$ (relevant for doing calculations beyond one-loop) there appear to be but two examples of performing the DeWitt expansion. In this paper we compute the off-diagonal elements of the DeWitt expansion coefficients using the Fock-Schwinger gauge. Our technique is based on representing $M_{xy}$ by a quantum mechanical path integral. We also generalize our method to the case of curved space, allowing us to determine the DeWitt expansion of \tilde M_{xy} = \langle x| \exp \case{1}{2} [\case{1}{\sqrt {g}} (\partial_\mu - i A_\mu)g^{\mu\nu}{\sqrt{g}}(\partial_\nu - i A_\nu) ] t| y \rangle by use of normal coordinates. By comparison with results for the DeWitt expansion of this matrix element obtained by the iterative solution of the diffusion equation, the relative merit of different approaches to the representation of $\tilde M_{xy}$ as a quantum mechanical path integral can be assessed. Furthermore, the exact dependence of $\tilde M_{xy}$ on some geometric scalars can be determined. In two appendices, we discuss boundary effects in the one-dimensional quantum mechanical path integral, and the curved space generalization of the Fock-Schwinger gauge.Comment: 16pp, REVTeX. One additional appendix concerning end-point effects for finite proper-time intervals; inclusion of these effects seem to make our results consistent with those from explicit heat-kernel method

Revenue optimization for the Ocean Grazer wave energy converter through storage utilization

Increased penetration of renewable energy generation motivates a change of paradigm in the way power systems are structured and operated, as advocated by the smart grid concept. Accordingly, in this paper we investigate the lossless storage capabilities of the Ocean Grazer wave energy converter (WEC), which could facilitate the aforementioned paradigm shift. This specific WEC exhibits both adaptability with respect to the incoming waves and significant lossless storage capabilities. We propose a model predictive control (MPC) strategy based on a lumped dynamical model in order to mitigate power imbalances in the power grid and maximize the revenue of the WEC. Furthermore, we illustrate that the proposed strategy exploits the WEC energy storage capabilities and we show the economic benefits it brings. Lastly, the proposed strategy is compared with a heuristic approach and a setting without storage