首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
In the first part of the paper [Molodenskii, 2011], we considered the problem of ambiguity in the solution of the inverse problem of retrieval of density distribution in the Earth’s core and mantle and determination of the Q factors in the mantle from the entire set of modern data on seismic velocities (V p and V S ), the frequencies f i and quality factors Q i of free oscillations of the Earth, and the amplitudes and phases of its forced nutations. We have constructed the model distributions of these parameters, in which the root-meansquared (rms) deviations of all observed values from the predicted ones are much smaller than in the PREM model. Below, we compare the observed amplitudes of the forced nutation with the values predicted by our model. In order to understand how rigid are the constraints imposed by the amplitudes of forced nutation, we not only calculate the deviations of the observed amplitudes of nutation from the predictions by our model but also estimate the changes in these deviations caused by small variations in several parameters of the model. To the parameters to be varied we refer those which have no or barely any effect on the periods and damping constants of free oscillations but have a pronounced effect on the amplitudes of forced nutation. These parameters include (1) the rheological properties of the mantle in the interval of periods from an hour to a day; (2) the dynamical flattening of the liquid core; (3) the dynamic flattening of the solid inner core; (4) the viscosity of the liquid core; and (5) the moment of inertia of the solid inner core. In addition, we estimate the effects of variations in the moment of inertia of the liquid core to be small (±0.2%) and not to affect, within the observation error, the periods of free oscillations. We show that the uncertainty in the model depth distributions of density considerably decreases when the new data on the amplitudes and phases of the forced nutation of the Earth are taken into account. With these data, it is possible to estimate the creep function for the lower mantle in a wide range of periods from a second to a day.  相似文献   

2.
The results of solving the inverse problem of forced nutations and free oscillations of the Earth by decomposing the Q-factor and small depth variations in density in a system of orthogonal functions are considered. These functions are determined by orthogonalization of the functional derivatives of the observed parameters with respect to the depth distributions of the sought parameters (assuming there are no distributions of the velocities of body seismic waves V p and V S with depth and unchanged total mass M and inertia moments I of the Earth). The examples are presented to illustrate the numerical solution of the inverse problem on finding the density distributions in the mantle and core of the Earth using orthogonalization of the integral constraints for the probable depth distributions of density describing the conditions of unchanged M and I, as well as the constraints posed by the data on the periods of the free low-order oscillations of the Earth.  相似文献   

3.
As was shown in [Molodensky, 2004a, 2004b], modern very long base interferometer (VLBI) data on the amplitudes and phases of the Earth’s forced nutation can provide significantly more rigid constraints on possible values of the quality factor of the lower mantle Q μ and on the dynamic flattening of the liquid core e lc as compared with seismic evidence and data on damping of the free oscillations of the Earth. On the other hand, the accuracy of modern tidal gravity data (obtained from twenty-year series of observations with a cryogenic gravimeter) is also very high and these data must be taken into account while estimating the parameters Q μ and e lc . The paper presents comparative estimates of the determination accuracy of the parameters Q μ and the dynamic flattening of the liquid core from VLBI and the aforementioned tidal gravity data.  相似文献   

4.
We employ a niching genetic algorithm to invert ∼30,000 differential ScS/S attenuation values for a new spherically symmetric radial model of shear quality factor (Qμ) with high sensitivity to the lower mantle. The new radial Qμ model, QLM9, possesses greater sensitivity to Qμ at large mantle depths than previous studies. On average, lower mantle Qμ increases with depth, which supports models of increasing viscosity with depth [B.M. Steinberger, A.R. Calderwood. Mineral physics constraints on viscous flow models of mantle flow, J. Conf. Abs., 6, 2001., 2001.]. There are two higher-Qμ regions at ∼1000 and ∼2500 km depth, which roughly correspond to high-viscosity regions observed by Forte and Mitrovica [A.M. Forte and J.X. Mitrovica, Deep-mantle high-viscosity flow and thermochemical structure inferred from seismic and geodynamic data, Nature 410, 1049–1056, 2001.]. There is a lower-Qμ layer at the core–mantle boundary and a relatively low-Qμ region in the mid-lower mantle. With several caveats, we infer a divergence of the solidus and geotherm in the lower mantle and a convergence within Dʺ by relating Qμ to homologous temperature.  相似文献   

5.
We analyze the present-day data on the periods of free oscillations and amplitudes of the forced nutations of the Earth for evaluating the admissible range of the mass and moment of inertia for the liquid core. The initial model for this study is taken in the form of the model distribution of density and mechanical Q parameters of the mantle suggested in (Molodenskii, 2010; 2011a; 2011b). This model was constructed by the steepest descent method in the space of 64 parameters, which determine the distribution of density and parameters of mechanical Q in the mantle, liquid outer core, and solid inner core of the Earth. We assumed the Q parameter of the mantle and inner solid core to be constant and sought for the density variations for the simplest two-parameter model of the piecewise-linear functions with the jumps on the boundary between the liquid core and the mantle and at the olivine-spinel phase transition at a depth of 670 km in the mantle. After this, the computations were repeated for the other distributions of Q (which were also assumed to be unchanged) that correspond to their limiting admissible values. Using this approach, we managed to find the most probable values of the mass and moment of inertia of the liquid core and determine the admissible range of their values. According to our estimates, the ratios of the mass and moments of inertia of the liquid core to the mass and moment of inertia of the whole Earth fall in the intervals 0.317996 ± 0.00065 and 0.110319 ± 0.00022, respectively. These values are lower than the corresponding values for the PREM model (0.322757 and 0.112297) by (1.48 ± 0.30)% and (1.76 ± 0.35)%, respectively. The interpretation of these results requires the revision and thorough analysis of the data on the admissible temperature range of the liquid core and (or) its chemical composition.  相似文献   

6.
Preliminary reference Earth model   总被引:29,自引:0,他引:29  
A large data set consisting of about 1000 normal mode periods, 500 summary travel time observations, 100 normal mode Q values, mass and moment of inertia have been inverted to obtain the radial distribution of elastic properties, Q values and density in the Earth's interior. The data set was supplemented with a special study of 12 years of ISC phase data which yielded an additional 1.75 × 106 travel time observations for P and S waves. In order to obtain satisfactory agreement with the entire data set we were required to take into account anelastic dispersion. The introduction of transverse isotropy into the outer 220 km of the mantle was required in order to satisfy the shorter period fundamental toroidal and spheroidal modes. This anisotropy also improved the fit of the larger data set. The horizontal and vertical velocities in the upper mantle differ by 2–4%, both for P and S waves. The mantle below 220 km is not required to be anisotropic. Mantle Rayleigh waves are surprisingly sensitive to compressional velocity in the upper mantle. High Sn velocities, low Pn velocities and a pronounced low-velocity zone are features of most global inversion models that are suppressed when anisotropy is allowed for in the inversion.The Preliminary Reference Earth Model, PREM, and auxiliary tables showing fits to the data are presented.  相似文献   

7.
The models of the mechanical Q factor of the inner solid core of the Earth are reconstructed from the amplitudes and phases of forced nutation and the periods and damping constants of the high-order overtones of free radial modes. The admissible range of the Q-factor in the solid core is estimated and the stability of the obtained distributions is analyzed. The real accuracy of the obtained model distributions is estimated on the basis of the previous conclusions concerning the uncertainty in the solution of the inverse problem of reconstructing the internal structure of the Earth in the low-frequency range.  相似文献   

8.
The existence of uncoupled shear (S) and compression (P) wave velocity variations in Earth's mantle is a characteristic that might only be explained by the presence of significant chemical and/or phase heterogeneity, with important implications for the dynamics and evolution of Earth's interior. While making a one-to-one comparison between tomographic models for P and S velocity (VP and VS) variations for a particular geographic region is ill-posed, their global statistical distributions reveal several robust characteristics indicative of the nature of uncoupled VP and VS in the deep mantle. We find that all of the VP and VS model distributions at a given depth are Gaussian-like throughout the lowermost mantle. However, a distinct low velocity feature is present in VS distributions below ≈ 2200 km depth that is not present or is relatively weak in VP models. The presence of anomalously low VS material cannot be explained as an artifact, nor can the absence of a similarly strong feature in P models be ascribed to under-resolution. We propose that this feature can be partly explained by laterally variable occurrences of post-perovskite (pPv) lenses in the D″ layer, however, the persistence of significantly slow VS regions at heights up to ≈ 700 km or more above the core–mantle boundary is likely to be incompatible with a pPv origin and might only be explained by the presence of a laterally discontinuous layer of chemically distinct material and/or some other kind of phase heterogeneity. There also exist significant discrepancies between tomographic models with respect to the width of the distributions as well as differences between the modeled peak values. We propose a scheme for comparison between different seismic models in which the widths of the dominant features in their statistical distributions is exploited.  相似文献   

9.
In the first part of the paper, we obtained the refined estimates for the periods and Q-factors of the fundamental modes and overtones of spherical and toroidal oscillations with periods longer than 3 min from the data on the free oscillations of the Earth, which were excited by the earthquakes with magnitude 9 that occurred in Sumatra, Japan, and the Sea of Okhotsk. In (Molodenskii et al., 2013), we analyzed the limits of the admissible density distributions in the mantle and liquid core of the Earth, using the data on the amplitudes and phases of the forced nutations, as well as the periods and attenuation factors of the fundamental modes of the free spheroidal and toroidal oscillations of the Earth. These studies were conducted with the fixed values of the total mass and total moment of inertia of the Earth and the fixed distributions of the body seismic waves in the mantle and in the core. The solution was obtained by orthogonalizing the kernels of the integral equations that link the residuals of the observed frequencies and attenuation factors of the free oscillations, as well as the amplitudes and phases of the forced nutations, with the sought densities and Q-factors of the mantle and liquid core. Below, we present the solution of the same problem with allowance for the results obtained in the first part of this paper, namely, the new data on the periods and attenuation factors of the fundamental modes of free oscillations of the Earth and on the periods of the first four overtones of the free spheroidal and toroidal oscillations. Despite the involvement of the new data on the overtones, which have not been considered in our calculations, the weighted root mean square deviations of the theoretical predictions from the observed periods and attenuation factors of the free oscillations, as well as the amplitudes and phases of the forced nutations, have significantly decreased. This is due to (1) the noticeable reduction of the real errors in estimating the parameters of the free oscillations described in the first part of the paper and (2) the inclusion of the quantities determining the depth- and frequency dependences of the Q-factor in the mantle in the set of the independently varied parameters.  相似文献   

10.
To study the amount of heat generated by radioactive decay in the continental crust, the usual practice in the literature is to fit to the heat flow and radioactivity data a relationship of the form: Q = Qr + D · A where Q and A are the observed heat flow and radiogenic heat production. Qr is the “reduced” heat flow and D is a depth scale. This procedure implicitly assumes that uranium, thorium and potassium have identical distributions in the crust. We suggest that significant information may be lost as the three radioelements may in fact be affected by processes operating over different depths.Data published for four heat flow provinces throughout the world are used to estimate the distributions of uranium, thorium and potassium in the continental crust. These distributions are characterized by a depth scales defined as follows: Di =∫0h Ci(z)Ci(0)dz where h is the thickness of the layer containing the bulk of radioactivity and Ci(z) the concentration of element i at depth z. Three depth scales are computed from a least-squares fit to the following relationship: Q = Qr + DU · AU + DT · AT + DK · AT where Q is the observed heat flow and Qr some constant (a reduced heat flow). Ai is the heat generation rate due to the radioactive decay of element i, and Di is the corresponding depth scale.The analysis suggests that the three distributions are different and that they have the same basic features in the four provinces considered. The depth scale for potassium is large in granitic areas, that for thorium is small and that for uranium lies between the other two.We propose a simple model according to which each radioelement essentially provides a record for one process. Potassium gives a depth scale for the primary differentiation of the crust. Thorium gives the depth scale of magmatic or metamorphic fluid circulation. Finally, the uranium distribution reflects the late effects of alteration due to meteoric water. We show that the heat flow and radioactivity data are compatible with this model.Our analysis and numerical results are supported by data from deep boreholes and by geochemical evidence, such as detailed investigations of plutonic series and studies of U-Th-Pb systematics.  相似文献   

11.
We discuss the problem of the ambiguity of gravity inversion, i.e., finding the depth distribution of density and the depth and frequency dependences of the Q-factor from the entire set of the present-day seismic and astrometric data on the travel times of seismic waves, the periods and attenuation factors of the free oscillations of the Earth, as well as the amplitudes and phases of the forced nutations. In the first part of the paper, we present the new and more accurate determinations of the periods and Q-factors for the fundamental tones and overtones of the spheroidal and toroidal oscillations of the Earth, which have periods longer than 3 min. These determinations are based on analyzing the signals from the Sumatra, Tohoku, and Okhotsk earthquakes of magnitude 9, which were recorded by the stations of the Global Seismographic Network (GSN) in Obninsk and Kazakhstan. It is shown that, although the Okhotsk earthquake had a lower magnitude than the other quakes analyzed (since its seismic source was extremely deep (about 600 km)), the amplitudes of the overtones excited by this event are significantly higher than the amplitudes of the overtones caused by the Sumatra and Tohoku events of magnitude 9. Moreover, the amplitudes of the overtones from the Okhotsk earthquake exceed the amplitudes of the overtones of the free oscillations caused by the other seismic events of magnitude 9 that occurred in the second half of the 20th century. Due to this, the data on the Okhotsk Sea earthquake are of utmost importance for the solution of the inverse problems of reconstructing the vertical profiles of Q-factor in the ultra-low frequency (ULF) range and for reconstructing the vertical distribution of density. Based on the new empirical data, we obtained new and more accurate estimates for the periods and attenuation factors of the free oscillations of the Earth.  相似文献   

12.
The VPREMOON seismic reference Moon model (Garcia et al., 2011) has been tested with respect to the thermal regime and chemical composition of the mantle. Based on a self-consistent thermodynamic approach and petrological models of the lunar mantle covering a wide range of concentrations of CaO, Al2O3, and FeO, we convert the P- and S-wave velocity profiles to the temperature–depth profiles. The solution procedure relies on the method of the Gibbs free energy minimization and the equations of state for the mantle material which take into account the effects of phase transformations, anharmonicity, and anelasticity. We find that regardless of the chemical composition, the positive P- and S-wave velocity gradient in the lunar mantle leads to a negative temperature gradient, which has no physical basis. For adequate mantle temperatures, the P- and S-wave velocities should remain almost constant or slightly decrease with depth (especially VS) as a result of the effects of the temperature, which grows faster than pressure. These findings underscore the importance of the relationship of the thermodynamics and physics of minerals with seismology.  相似文献   

13.
Based on the long period surface wave data recorded by the China Digital Seismograph Network (CDSN), theQ R of fundamental mode Rayleigh wave with periods from 10 s to 146 s is determined for the eastern Sino-Korean paraplatform in this paper. TheQ β models of the crust and upper mantle are respectively obtained for the 4 paths, with the aid of stochastic inverse method. It shows that in the eastern Sino-Korean paraplatform, the average crustalQ β is about 200, and that there exists a weak attenuation layer in the middle crust (about 10–20 km deep) which is possibly related to earthquake-prone layer. A strong attenuation layer (lowQ) of 70 km thick extensively exists in the uppermost mantle, with the buried depth about 80 km. The averageQ R of fundamental mode Rayleigh wave is between the value of stable tectonic region and that of active tectonic region, and much close to the latter.  相似文献   

14.
The motion of the Earth’s geophysical pole is modeled by numerical integration of the Liouville equations. The dependence of variable parameters of these equations (the period of free nutation T and the quality factor of the shell Q) is examined as a function of time and initial data sets used for constructing geophysical perturbing functions. In particular, we used oceanic perturbing functions constructed from TOPEX/POSEIDON altimetry data. The annual and Chandler components of the pole motion were also modeled. Apart from the subtraction of quadratic trends, initial data were not filtered or averaged. The performed analysis provided model values of T = 425–440 days and Q = 20–60, and Q was found to be somewhat unstable with time.  相似文献   

15.
The aim of this paper is to show the application of short-period surface waves recorded during deep seismic sounding experiment for constraining shallow velocity structure of the crust. Phase velocity of fundamental mode Rayleigh waves, observed along the CELEBRATION 2000 experiment profile CEL09, were obtained by a p-ω method and has been subsequently inverted for one-dimensional shear velocity models for the top 2 km. Multiple filter technique applied to one shot gather was used to carry out a joint inversion of phase and group velocity data and to provide γR data to be used for Qβ inversion. Validity of obtained VS and Qβ models was confirmed by the reflectivity method. Noticeably, no clear dispersive wawes were observed in the Tepla-Barrandian Unit. Quasi-2D model based on the individual 1D VS models is well correlated with the surface geology. Lower VS are observed in the Saxothuringian Zone in comparison to the Moldanubian Zone. In the vicinity of the Central Bohemian and Moldanubian Plutons, the near-surface VS values are relatively low, but below 1 km depth, they are higher than in surrounding areas. We interpret it as the result of the weathering and cracks within the granitoid rocks.  相似文献   

16.
The “anomalous” layer in the lowermost mantle, identified as D″ in the notation of K.E. Bullen, appears in the PREM Earth model as a 150 km-thick zone in which the gradient of incompressibility with pressure, dKdP, is almost 1.6, instead of 3.2 as in the overlying mantle. Since PREM shows no accompanying change in density or density gradient, we identify D″ as a thermal boundary layer and not as a chemically distinct zone. The anomaly in dKdP is related to the temperature gradient by the temperature dependence of Ks, for which we present a thermodynamic identity in terms of accessible quantities. This gives the numerical result (?Ks/?T)P=?1.6×107 Pa K?1 for D″ material. The corresponding temperature increment over the D″ range is 840 K. Such a layer cannot be a static feature, but must be maintained by a downward motion of the lower mantle toward the core-mantle boundary with a strong horizontal flow near the base of D″. Assuming a core heat flux of 1.6 × 1012 W, the downward speed is 0.07 mm y?1 and the temperature profile in D″, scaled to match PREM data, is approximately exponential with a scale height of 73 km. The inferred thermal conductivity is 1.2 W m?1 K?1. Using these values we develop a new analytical model of D″ which is dynamically and thermally consistent. In this model, the lower-mantle material is heated and softened as it moves down into D″ where the strong temperature dependence of viscosity concentrates the horizontal flow in a layer ~ 12 km thick and similarly ensures its removal via narrow plumes.  相似文献   

17.
Consider a plane homogeneous harmonic SH wave incident upon an interface between two anelastic half-spaces. Computing the plane wave displacement and energy-flux-based reflection and transmission coefficients correctly requires determining the proper signs of the vertical slowness components of all the reflected and transmitted waves, i.e., determining which of the two values of the square root for a given vertical slowness should be chosen. For anelastic media, this can be problematic, as unphysical results can arise. Previous research has led to a specific recommendation on how to choose the signs. However, when this recommendation is employed, it is found via numerical experimentation that for certain values of the medium parameters, the energy-flux-based transmission coefficient T can be negative for certain supercritical values of the incidence angle, whereas physical reasoning suggests it should be zero. To investigate this seemingly unphysical result, an analytical determination of the mathematical conditions under which it occurs would be useful. Such a determination is performed in this article. Letting V 1 and V 2 be the wave speeds of homogeneous SH plane waves in the incidence and transmission media respectively, and Q 1 and Q 2 be the corresponding quality factors, with Q 1, Q 2 ? 1, it is found that if V 1 < V 2 and Q 1 < Q 2 (a common situation in the Earth), then T will be negative for part of the supercritical incidence angle range if 1 < Q 2/Q 1 ?? 2 ? (V 1/V 2)2 and for all of it if 2 (V 1/V 2)2 < Q 2/Q 1.  相似文献   

18.
A second-order hydrostatic theory is developed on the assumption that the trace of the Earth's inertia tensor, its mass and mean radius are invariant under any process causing deviations from the hydrostatic state.The hydrostatic flattening and the zonal coefficients of the hydrostatic gravitational field are obtained as ??1 = 299.638, J2 = 1072.618 × 10?6 and J4 = ?2.992 × 10?6, respectively.The internal theory using the preliminary reference earth model (PREM) of Dziewonski and Anderson (1981) yields ??1 = 299.627, J2 = 1072.701 × 10?6 and J4 = ?2.992 × 10?6. The agreement between these and the hydrostatic values indicate that PREM is suitable as a reference model as it represents the spheroidal density distribution in a state of zero non-hydrostatic stress while satisfying the fundamental geodetic observations of the invariant quantities.The small discrepancy between the hydrostatic flattening and the value deduced from PREM suggests that the density is underestimated at large depths and/or it is slightly overestimated in shallow regions of the Earth.The discrepancies between the hydrostatic and observed quantities persist after the removal of the accountable effects of isostatically compensated topography, permanent tidal deformation and the present mass anomalies associated with the Late-Pleistocene deglaciation. These ‘corrected’ discrepancies point to a triaxial non-hydrostatic figure which cannot be explained by the delayed response of the Earth to tidal deceleration.  相似文献   

19.
The attenuation characteristics of Indian lithosphere and its comparison with different tectonic settings in the world are determined from the observations of the Q for Lg(QLg)-, and S(QS)-waves in the 1-30 Hz frequency range. The scattering is approximated with a Gaussian distribution of spherical scatterers. To approximate single scattering, we use Dainty's [Geophy. Res. Lett. 8 (11) (1981) 1126] model that attenuation is given by 1/Q(ω) = 1/Qi + g(ω)v/ω, where Qi is intrinsic Q due to anelastic attenuation, v is shear wave velocity, ω is angular frequency, g = ∫n(a)σ da is the total scattering coefficient for S-to-S scattering, n(a) da is the number of scattering spheres of radius a per unit volume, and σ is the scattering cross-section for the sphere. We find that if n(a) is described by a simple two parameter (a0 and c) Gaussian of amplitude c and standard deviation and mean a0, the attenuation data for different regions of the world are well approximated over the frequency band of seismic observations. Our major findings are: (1) the maximum effect of scattering on attenuation occurs at 0.84 Hz or a wavelength of 4.16 km; (2) the values of g are frequency dependent. Values of g are of the order of 10−3 km−1 at 1-30 Hz, varying from 0.0031 to 0.01 and 0.001 to 0.0083 km−1 for tectonically active and stable regions, respectively; (3) regions of active tectonics and seismicity generally have lower Qi values (1000) than that in stable regions (2000); and (4) regions of high Qi value exhibit low intensity of scattering.  相似文献   

20.
This paper deals with characteristics of the short period S-wave attenuation field in the rupture zones of 37 large and great earthquakes with M s = 7.0–8.6, as well as in low seismicity areas. We estimate the effective quality factor from Sn and Lg coda envelopes in two time intervals (Q 1 and Q 2). The quantity Q 1 is a measure of shear wave attenuation in the uppermost mantle, at depths of down to approximately 200–250 km, while Q 2 is relevant to deeper horizons of the upper mantle. We studied variations in the attenuation field in the rupture zone of the 1950 Assam earthquake. We examined the parameters Q 1, Q 2, and Q 1/Q 2 as functions of the time ΔT elapsed after a large earthquake. It is shown that the parameter Q 2 in rupture zones is practically independent of ΔT, while the quantities Q 1 and Q 1/Q 2 increase until ΔT ~ 20–25 years, especially rapidly for normal, normal-oblique, and strike-slip earthquake mechanisms. This analysis provides evidence that, as ΔT increases, so does the quality factor in the upper mantle for shear waves. It is supposed that this is related to the rise of mantle fluids to the crust. Geodynamic mechanisms are discussed that can support a comparatively rapid “drying” of the upper mantle beneath earthquake rupture zones.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号