Validity of Resolving the 785 km Discontinuity in the Lower Mantle with P′P′ Precursors?

2020 ◽  
Vol 91 (6) ◽  
pp. 3278-3285
Author(s):  
Baolong Zhang ◽  
Xiangfang Zeng ◽  
Jun Xie ◽  
Vernon F. Cormier

Abstract P ′ P ′ precursors have been used to detect discontinuities in the lower mantle of the Earth, but some seismic phases propagating along asymmetric ray paths or scattered waves could be misinterpreted as reflections from mantle discontinuities. By forward modeling in standard 1D Earth models, we demonstrate that the frequency content, slowness, and decay with distance of precursors about 180 s before P′P′ arrival are consistent with those of the PKPPdiff phase (or PdiffPKP) at epicentral distances around 78° rather than a reflection from a lower mantle interface. Furthermore, a beamforming technique applied to waveform data recorded at the USArray demonstrates that PKPPdiff can be commonly observed from numerous earthquakes. Hence, a reference 1D Earth model without lower mantle discontinuities can explain many of the observed P′P′ precursors signals if they are interpreted as PKPPdiff, instead of P′785P′. However, this study does not exclude the possibility of 785 km interface beneath the Africa. If this interface indeed exists, P′P′ precursors at distances around 78° would better not be used for its detection to avoid interference from PKPPdiff. Indeed, it could be detected with P′P′ precursors at epicentral distances less than 76° or with other seismic phases such as backscattered PKP·PKP waves.

1998 ◽  
Vol 132 (3) ◽  
pp. 603-642 ◽  
Author(s):  
C. Denis ◽  
M. Amalvict ◽  
Y. Rogister ◽  
S. Tomecka-Suchoń

SUMMARY After general comments (Section 1) on using variational procedures to compute the oblateness of internal strata in the Earth and slowly rotating planets, we recall briefly some basic concepts about barotropic equilibrium figures (Section 2), and then proceed to discuss several accurate methods to derive the internal flattening. The algorithms given in Section 3 are based on the internal gravity field theory of Clairaut, Laplace and Lyapunov. They make explicit use of the concept of a level surface. The general formulation given here leads to a number of formulae which are of both theoretical and practical use in studying the Earth's structure, dynamics and rotational evolution. We provide exact solutions for the figure functions of three Earth models, and apply the formalism to yield curves for the internal flattening as a function of the spin frequency. Two more methods, which use the general deformation equations, are discussed in Section 4. The latter do not rely explicitly on the existence of level surfaces. They offer an alternative to the classical first-order internal field theory, and can actually be used to compute changes of the flattening on short timescales produced by variations in the LOD. For short durations, the Earth behaves elastically rather than hydrostatically. We discuss in some detail static deformations and Longman's static core paradox (Section 5), and demonstrate that in general no static solution exists for a realistic Earth model. In Section 6 we deal briefly with differential rotation occurring in cylindrical shells, and show why differential rotation of the inner core such as has been advocated recently is incompatible with the concept of level surfaces. In Section 7 we discuss first-order hydrostatic theory in relation to Earth structure, and show how to derive a consistent reference Earth model which is more suitable for geodynamical modelling than are modern Earth models such as 1066-A, PREM or CORE11. An important result is that a consistent application of hydrostatic theory leads to an inertia factor of about 0.332 instead of the value 0.3308 used until now. This change automatically brings ‘hydrostatic’ values of the flattening, the dynamic shape factor and the precessional constant into much better agreement with their observed counterparts than has been assumed hitherto. Of course, we do not imply that non-hydrostatic effects are unimportant in modelling geodynamic processes. Finally, we discuss (Sections 7–8) some implications of our way of looking at things for Earth structure and some current problems of geodynamics. We suggest very significant changes for the structure of the core, in particular a strong reduction of the density jump at the inner core boundary. The theoretical value of the free core nutation period, which may be computed by means of our hydrostatic Earth models CGGM or PREMM, is in somewhat better agreement with the observed value than that based on PREM or 1066-A, although a significant residue remains. We attribute the latter to inadequate modelling of the deformation, and hence of the change in the inertia tensor, because the static deformation equations were used. We argue that non-hydrostatic effects, though present, cannot explain the large observed discrepancy of about 30 days.


Satellite orbits contracting under the influence of air drag experience 15th-order resonance when the track over the Earth repeats after 15 revolutions. If the orbital decay rate is slow enough, an orbit passing through the resonance is appreciably perturbed by the effects of 15th-order harmonics in the geopotential. We have used the observed perturbations in 23 resonant orbits, at various inclinations to the equator, to determine the harmonic coefficients of order 15 and degree 15, 16, 17,... 35. Analysis of the changes in orbital inclination on the 23 orbits gives the harmonics of odd degree, while those of even degree are found from the changes in eccentricity on 16 of the orbits. The values derived are given in tables 6 and 8. The coefficients of degrees 15, 16, 17,... 23, should be more accurate than any previously obtained; their average s. d. is 1.4 x 10 -9 , equivalent to 1 cm in geoid height. Comparisons with comprehensive Earth models show the Goddard Earth Model 10B to be the best, and a standard deviation of about 3 x 10 -9 in the GEM 10B 15th-order coefficients is indicated.


1992 ◽  
Vol 82 (3) ◽  
pp. 1453-1474
Author(s):  
J. F. Cassidy

Abstract The use of broadband receiver function analysis to estimate the fine-scale S-velocity structure of the lithosphere is becoming increasingly popular. A series of numerical experiments shows several important aspects of this technique, with emphasis on estimation of dipping interfaces. The recent modification introduced to the receiver function analysis technique that preserves absolute amplitudes (Ammon, 1991) is more robust than the previous technique of modeling receiver functions that were normalized to unit amplitude. Using the latter method, shallow (e.g., depths less than ∼2 km) high-velocity contrast interfaces may alter the apparent amplitudes of Ps phases and produce inaccuracies in the Earth model developed. The use of absolute amplitudes minimizes this potential for error. When research targets include deep dipping structure, tight stacking bounds (e.g., ≦ 10° in backazimuth (BAZ) and epicentral distance (Δ)) should be applied to avoid attenuating Ps phases and to aid in the identification of reverberations or scattered energy. Reverberations sample a relatively large lateral range about the recording site (e.g., a radius of 1 to 1.5 times the depth of the reflecting interface) and in the presence of dipping interfaces exhibit drastic variations in amplitude and arrival time as a function of BAZ and Δ. Thus, they cannot readily be used to provide constraints on the Earth structure. Formal inversion techniques, which attempt to match all arrivals in the waveform, must be used with caution when modeling receiver functions from complex regions. Only those phases whose amplitude and arrival-time variations as a function of BAZ and Δ are consistent with those of Ps conversions should be modeled. Forward modeling may resolve, depending upon the data quality and noise level, S-velocity contrasts greater than ∼ 0.2 to 0.4 km / sec. Layers of thickness 2 to 5 km may be accurately imaged, and transition zones may be examined by considering various frequency bands of the data. In order to better understand the resolving power of the data, the averaging functions associated with the receiver functions may be calculated from the observed data and, if desired, used in the forward modeling process.


Geophysics ◽  
1987 ◽  
Vol 52 (10) ◽  
pp. 1412-1417 ◽  
Author(s):  
Sri Niwas ◽  
M. Israil

A computational method has been developed based on a linear relationship between apparent resistivity and the kernel function through a matrix operator. Forward modeling of dipole apparent resistivity values over a one‐dimensional layered earth model as well as the inversion of electrical sounding data can be easily accomplished using this matrix operator. For a few two‐layer, three‐layer, and four‐layer earth models, the apparent resistivity values can be calculated quickly and accurately. It is shown that the technique is adaptable to automatic interpretation in both the resistivity and the kernel domains.


2010 ◽  
Vol 63 (3) ◽  
pp. 395-417 ◽  
Author(s):  
Kenneth Gade

Position calculations, e.g. adding, subtracting, interpolating, and averaging positions, depend on the representation used, both with respect to simplicity of the written code and accuracy of the result. The latitude/longitude representation is widely used, but near the pole singularities, this representation has several complex properties, such as error in latitude leading to error in longitude. Longitude also has a discontinuity at ±180°. These properties may lead to large errors in many standard algorithms. Using an ellipsoidal Earth model also makes latitude/longitude calculations complex or approximate. Other common representations of horizontal position include UTM and local Cartesian ‘flat Earth’ approximations, but these usually only give approximate answers, and are complex to use over larger distances. The normal vector to the Earth ellipsoid (called n-vector) is a non-singular position representation that turns out to be very convenient for practical position calculations. This paper presents this representation, and compares it with other alternatives, showing that n-vector is simpler to use and gives exact answers for all global positions, and all distances, for both ellipsoidal and spherical Earth models. In addition, two functions based on n-vector are presented, that further simplify most practical position calculations, while ensuring full accuracy.


2013 ◽  
Vol 194 (1) ◽  
pp. 362-366 ◽  
Author(s):  
Yingjie Xia ◽  
Sidao Ni ◽  
Xiangfang Zeng

Abstract Based on studies of continuous waveform data recorded on broad-band seismograph stations in Africa, Europe and North America, we report evidences for two temporally persistent and spatially localized monochromatic vibrating sources (around 0.036 and 0.038 Hz, respectively) in the Gulf of Guinea, instead of just one source (0.038 Hz or 26 s) found 50 yr ago. The location of the 0.036 Hz source is close to the Sao Tome Volcano, therefore it may be related to volcano processes. However, the 0.038 Hz source cannot be explained with known mechanisms, such as tectonic or oceanic processes. The most likely mechanism is volcano processes, but there is no reported active volcano in source region. Such repetitive vibration sources may provide valuable tools for detecting temporal variation of crustal structure of the Earth.


1963 ◽  
Vol 53 (3) ◽  
pp. 483-501 ◽  
Author(s):  
Leonard E. Alsop

Abstract Periods of free vibrations of the spheroidal type have been calculated numerically on an IBM 7090 for the fundamental and first two shear modes for periods greater than about two hundred seconds. Calculations were made for four different earth models. Phase and group velocities were also computed and are tabulated herein for the first two shear modes. The behavior of particle motions for different modes is discussed. In particular, particle motions for the two shear modes indicate that they behave in some period ranges like Stoneley waves tied to the core-mantle interface. Calculations have been made also for a model which presumes a solid inner core and will be discussed in Part II. The two computer programs which were made for these calculations are described briefly.


1969 ◽  
Vol 59 (5) ◽  
pp. 2079-2099
Author(s):  
John S. Derr

abstract All observations of the free oscillations of the Earth published through 1968 are weighted to produce a set of means and standard errors of the means. Fundamental orders 0 to 97 for spheroidal and 2 to 99 for torsional are treated, as well as many overtones up to order 49. Statistical tests indicate that some observations are path dependent at the 99 per cent confidence level. Comparison of these means and standard errors with published Earth models indicate that they form a consistent basis for inversion of free oscillation observations to infer Earth structures.


Author(s):  
A. M. Dziewonski ◽  
D. L. Anderson
Keyword(s):  

2020 ◽  
Vol 33 (2) ◽  
pp. 443-459 ◽  
Author(s):  
E. Powell ◽  
N. Gomez ◽  
C. Hay ◽  
K. Latychev ◽  
J. X. Mitrovica

AbstractThe West Antarctic Ice Sheet (WAIS) overlies a thin, variable-thickness lithosphere and a shallow upper-mantle region of laterally varying and, in some regions, very low (~1018 Pa s) viscosity. We explore the extent to which viscous effects may affect predictions of present-day geoid and crustal deformation rates resulting from Antarctic ice mass flux over the last quarter century and project these calculations into the next half century, using viscoelastic Earth models of varying complexity. Peak deformation rates at the end of a 25-yr simulation predicted with an elastic model underestimate analogous predictions that are based on a 3D viscoelastic Earth model (with minimum viscosity below West Antarctica of 1018 Pa s) by ~15 and ~3 mm yr−1 in the vertical and horizontal directions, respectively, at sites overlying low-viscosity mantle and close to high rates of ice mass flux. The discrepancy in uplift rate can be reduced by adopting 1D Earth models tuned to the regional average viscosity profile beneath West Antarctica. In the case of horizontal crustal rates, adopting 1D regional viscosity models is no more accurate in recovering predictions that are based on 3D viscosity models than calculations that assume a purely elastic Earth. The magnitude and relative contribution of viscous relaxation to crustal deformation rates will likely increase significantly in the next several decades, and the adoption of 3D viscoelastic Earth models in analyses of geodetic datasets [e.g., Global Navigation Satellite System (GNSS); Gravity Recovery and Climate Experiment (GRACE)] will be required to accurately estimate the magnitude of Antarctic modern ice mass flux in the progressively warming world.


Sign in / Sign up

Export Citation Format

Share Document