Lateral density anomaly distribution in the lower mantle obtained from satellite low order geopotential coefficients

In this paper, we attempt to use satellite gravity data and a new inversion method to study the lateral density anomaly distribution in the mantle. First, density difference Δρ(τ,θ,φ) is expanded in terms of a three—dimensional orthogonal function system, the coefficients of the expansion are to be determined. Then, a set of observation equations is established from the relationship between density anomaly and disturbing geopotential. In the equations the unknown vector contains the coefficients of density anomaly expansion, the observational vector is obtained by computing geopotential perturbations using the potential coefficients of GEM10B, and a filtering process is done for the observational values by properly selecting the harmonic degrees of geopotentical. Finally, the lateral density variations in the lower mantle (670 km toCM boundary) are investigated. In this case, the degrees of disturbing geopotential are selected as 2–11, the truncated degrees of density anomaly expansion are taken asL=6 andK=4, and the damping least squares method is used to solve the observation equations. The resulting model shows the high level of density perturbations at 670 km discontinuity and core — mantle boundary, a high — density zone circumscribing the Pacific and a lower—density region under the center of Pacific. These features are in agreement with the three—dimensional seismic velocity variation features by Dziewonski (1984). In the Antarctic region and some parts of Atlantic and Indian Ocean, however, the resulting density anomalies are negatively correlated with the seismic velocity anomalies, the cause resulting in these phenomena is preliminarily analysed in this paper. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 53–65, 1991. The principle and method represented in this paper can also be suitable to study the lateral density anomaly distribution in the earth’s crust and the upper mantle.     

Inversion of seismic and gravity data for the composition and core sizes of the Moon

We model the internal structure of the Moon, initially homogeneous and later differentiated due to partial melting. The chemical composition and the internal structure of the Moon are retrieved by the Monte-Carlo inversion of the gravity (the mass and the moment of inertia), seismic (compressional and shear velocities), and petrological (balance equations) data. For the computation of phase equilibrium relations and physical properties, we have used a method of minimization of the Gibbs free energy combined with a Mie-Gr@uneisen equation of state within the CaO-FeO-MgO-Al₂O₃-SiO₂ system. The lunar models with a different degree of constraints on the solution are considered. For all models, the geophysically and geochemically permissible ranges of seismic velocities and concentrations in three mantle zones and the sizes of Fe-10%S core are estimated. The lunar mantle is chemically stratified; different mantle zones, where orthopyroxene is the dominant phase, have different concentrations of FeO, Al₂O₃, and CaO. The silicate portion of the Moon (crust + mantle) may contain 3.5–5.5% Al₂O₃ and 10.5–12.5% FeO. The chemical boundary between the middle and the lower mantle lies at a depth of 620–750 km. The lunar models with and without a chemical boundary at a depth of 250–300 km are both possible. The main parameters of the crust, the mantle, and the core of the Moon are estimated. At the depths of the lower mantle, the P and S velocities range from 7.88 to 8.10 km/s and from 4.40 to 4.55 km/s, respectively. The radius of a Fe-10%S core is 340 ± 30 km.     

The structure of 2D mantle convection and stress fields: Effects of viscosity distribution

Spatial fields of temperature, velocity, overlithostatic pressure, and horizontal stresses in the Earth’s mantle are studied in two-dimensional (2D) numerical Cartesian models of mantle convection with variable viscosity. The calculations are carried out for three different patterns of the viscosity distribution in the mantle: (a) an isoviscous model, (b) a four-layer viscosity model, and (c) a temperature- and pressure-dependent viscosity model. The pattern of flows, the stresses, and the surface heat flow are strongly controlled by the viscosity distribution. This is connected with the formation of a cold highly viscous layer on the surface, which is analogous to the oceanic lithosphere and impedes the heat transfer. For the Rayleigh number Ra = 10⁷, the Nusselt number, which characterizes the heat transfer, is Nu = 34, 28, and 15 in models with constant, four-layered, and p, T-dependent viscosity, respectively. In all three models, the values of overlithostatic pressure and horizontal stresses σ _xx in a vast central region of the mantle, which occupies the bulk of the entire volume of the computation domain, are approximately similar, varying within ±5 MPa (±50 bar). This follows from the fact that the dimensionless mantle viscosity averaged over volume is almost similar in all these models. In the case of temperature- and pressure-dependent viscosity, the overlithostatic pressure and stress σ _xx fields exhibit much stronger concentration towards the horizontal boundaries of the computation domain compared to the isoviscous model. This effect occurs because the upwellings and downwellings in a highly viscous region experience strong variations in both amplitude and direction of flow velocity near the horizontal boundaries. In the models considered with the parameters used, the stresses in the upper and lower mantle are approximately identical, that is, there is no denser concentration of stresses in the upper or lower mantle. In contrast to the overlithostatic pressure field, the fields of horizontal stresses σ _xx in all models do not exhibit deep roots of highly viscous downwelling flows.     

Estimation of seismic stress drop from the peak velocity of ground motion

ＥｓｔｉｍａｔｉｏｎｏｆｓｅｉｓｍｉｃｓｔｒｅｓｄｒｏｐｆｒｏｍｔｈｅｐｅａｋｖｅｌｏｃｉｔｙｏｆｇｒｏｕｎｄｍｏｔｉｏｎＪＩＡＺＨＥＮＧＱＩＮ（秦嘉政）ＺＵＹＩＮＬＩＵ（刘祖荫）ＸＩＡＯＤＯＮＧＱＩＡＮ（钱晓东）ＱＩＮＧＹＩＮＸＩＥ（谢庆...     

Nonlinear free thermal convection in a spherical shell:A variable viscosity model

Introduction The velocity field of surface plate motion can be split into a poloidal and a toroidal parts.At the Earth′s surface,the toroidal component is manifested by the existence of transform faults,and the poloidal component by the presence of convergence and divergence,i.e.spreading and subduc-tion zones.They have coupled each other and completely depicted the characteristics of plate tec-tonic motions.The mechanism of poloidal field has been studied fairly clearly which is related to …     

Relationship between the ionospheric potential and the ground-level electric field in the southern polar cap

The relationships between the average hourly values of the vertical ground-level electric field measured at the Vostok Antarctic station and the ionospheric potential above the station have been obtained. The ΔEz and Uext parts of both parameters controlled by the solar wind were considered. Convection models (Weimer, 1995; Lukianova and Christiansen, 2006) and a model based on the SuperDARN radar system were used to determine the ionospheric potential. An analysis has been performed for isolated days and the entire sample in 1998–2000 (including 170 days of “fine weather”). For an isolated day, the best correlation coefficients (R) between ΔEz and Uext obtained using the three models were 0.81, 0.80, and 0.88, respectively. The total correlation coefficient for the entire data set was R = 0.24−0.32. The R value was larger during daytime (R ≈ 0.4) and smaller at night (R ≈ 0.1) and slightly increased in the early morning hours. The specific features of daily variations in R apparently indicate that it is possible to adequately describe the structure of the ionospheric electric field equipotentials by using large-scale stationary convection models. The R value varies complexly, depending on the IMF orientation, but it generally tends to increase from IMF By < 0 to By > 0, which is explained by the asymmetric convection patterns for opposite By signs.     

Geochemical constraints on the model of the composition and thermal conditions of the Moon according to seismic data

A new method of reconstruction of the temperature profile in the lunar mantle from the velocities of seismic P- and S-waves for different models of chemical composition is developed. The procedure of the solution of an inverse problem is realized with the help of the minimization of the Gibbs free energy and the equations of state of a mantle substance, taking into account phase transformations, anharmonicity, and the effects of inelasticity. The geophysical and geochemical constraints on composition and temperature distribution in Moon’s mantle are established. The upper mantle can be composed of olivine pyroxenite, depleted by low-volatile oxides (∼2 wt % of CaO and Al₂O₃). On the contrary, the lower mantle must be enriched by low-volatile oxides (∼4–6 wt % of CaO and Al₂O₃). Its composition can be represented by a mineral association of the olivine + clinopyroxene + garnet or olivine + orthopyroxene + clinopyroxene + garnet type, which is close in composition to pyrolite. The temperature distribution at depths 50–1000 km are approximated by the equation: T(°C) = 351 + 1718[1–exp (−0.00082H)]. The constraints inferred make it possible to conclude that the published values of the velocities of P- and S-waves for the lunar mantle, obtained by processing the data of seismic experiments of the Apollo lunar mission are inconsistent with each other at depths below 300 km. Otherwise, the variations in the velocities of P- and S-waves disturb the symmetry between the petrological model (composition), the temperature profile, and the seismic profile.     

Sr-Nd-Pb isotopes of the Early Paleozoic maficultramafic dykes and basalts from South Qinling belt and their implications for mantle composition

Late Early Paleozoic mafic-ultramafic dykes and volcanic rocks from the South Qinling belt are char- acterized by εNd( t ) = 3.28― 5.02, (87Sr/86Sr)i= 0.70341― 0.70555, (206Pb/204Pb)i = 17.256― 18.993, (207Pb/204Pb)i= 15.505―15.642, (208Pb/204Pb)i=37.125―38.968, ?8/4=21.18―774.43, ?7/4=8.11―18.82. These charac- teristics suggest that they derived from a Middle Neoproterozoic mantle with isotopic compositions of mixed HIMU, EMII and minor EMI components. We interpret that these rocks were melting products of depleted mantle modified by subducted ancient oceanic crust and continental margin sediments along the northern margin of Yangtze block during Early Neoproterozoic.     

Helium and argon isotopes of the Tertiary basic igneous rocks from Shandong Peninsula and implications for the magma origin

Helium (He) and Argon (Ar) isotopic compositions of the Tertiary basic igneous rocks were determined by the high temperature melting extraction method. The selected samples for the studies included al-kaline basalts and diabases from the Jiyang basin,and the surrounding Shanwang and Qixia outcrops in the Shandong Peninsula,eastern China. The results show that the Paleogene basalts and diabases from the Jiyang basin yielded a wide range of P4 PHe abundance of (73.70-804.16)×10 P-8 Pcm P3 P STP·g P-1 P,with P3 PHe/ P4 PHe ratios of 0.374-2.959 Ra,which was lower than the MORB but evidently higher than the con-tinental crust value. The Neogene alkaline basalts from the Jiyang basin,Shanwang and Qixia outcrops have variable P4 PHe abundances ((42.34-286.72)×10-8 Pcm P3 P STP·g-1 P),and "continental crust-like" P3 PHe/ P4 PHe ratios (0.013-0.074 Ra). All of them contain atmospheric-like P40 PAr/ P36 PAr ratio (395.4-1312.7),reflecting the mantle sources with air components. Their low P3 PHe/ P4 PHe ratios are interpreted as the enrichment of the radiogenic P4 PHe mainly inherited from the mantle. He and Ar systematics show the mixing of MORB-type,air and a P4 PHe enriched member in the mantle source,suggesting that these igneous rocks originated from the depleted asthenospheric mantle mixed with an EMI component. Therefore,the present He and Ar isotopes do not support the viewpoints that the Cenozoic igneous rocks of Eastern North China were the products of mantle plume(s) activities.     

Stress Changes due to Recent Seismic Events in the Central Apennines (Italy)

The Central Apennines, Italy, are characterized by moderate seismic activity on normal faults, oriented in directions parallel to the Apenninic chain. The subject of this study is the Umbria-Marche Apennines, a segment approximately 200-km long, where three main seismic events occurred in the last three decades. The 1979 Norcia earthquake was a M_w = 5.8 event, taking place at the south end of the considered segment. The 1984 Gubbio earthquake was a M_w = 5.6 event which took place at the north end. The 1997-1998 Colfiorito sequence constituted 8 main shocks with magnitudes M_w between 5 and 6 and epicenters located between the Gubbio and the Norcia earthquake areas. A model made of an elastic half-space is considered, in which the seismic sources are represented by rectangular dislocations which have the appropriate values of source parameters, and in which the static stress field produced by each event is calculated. The analysis of the Coulomb stress change (ΔC) as a function of time shows that the coseismic stress transfer and fault interaction played an important role in the region during the past three decades: 7 earthquakes of the 9 considered took place where ΔC>0. Such an interaction has been confirmed by the analysis of the aftershocks in the Colfiorito zone post September 26, 1997: about the 61% of the aftershocks considered took place where ΔC>0. The comparison between the ΔC^s due to the coseismic stress transfer and the rate ΔĊ^t due to the tectonic stress allows us to quantify the time advance of the earthquakes. The ΔC^s pattern shows positive values in two areas that can be regarded as historical seismic gaps.     

Study of calibration function for surface wave magnitude of DK1 seismographs

ＳｔｕｄｙｏｆｃａｌｉｂｒａｔｉｏｎｆｕｎｃｔｉｏｎｆｏｒｓｕｒｆａｃｅｗａｖｅｍａｇｎｉｔｕｄｅｏｆＤＫ１ｓｅｉｓｍｏｇｒａｐｈｓＦＥＮＧＸＵＥ（薛峰）ＹＯＮＧＺＨＡＯ（赵永）ＣｅｎｔｅｒｆｏｒＡｎａｌｙｓｉｓａｎｄＰｒｅｄｉｃｔｉｏｎ，Ｓｔａｔ...     

Lg magnitude determination for the eastern six Provinces of China (part 2)

Any pair ofm _{L gz} ,m _{L g} b,mmxz andmm x h in the eastern six provinces of China show good linear relation this makes it easy to convert from one scale to another.mm x h may replaceM _L (ECH) as the local earthquake magnitude scale for the eastern six provinces, yet retaining the level ofM _L (ECH). The scalemm x h is 0.36 magnitude units higher thanM _L (SC) for southern California. By comparingRm a x (Δ) withR ₃ (Δ) of Yunnan, it is found that the γ-value ofMAX phase of about 1 section in Yunnan region is half as large as that in the eastern six provinces of China. Observation indicates thatmm x h≈mb for magnitude 4–5. It follows therefore thatmb (ECH)≈mb (WUS). This is supported by the result that the attenuation coefficient, the γ-value in eastern China is the same as in western United States.     

Experimental simulations of gas-driven eruptions: kinetics of bubble growth and effect of geometry

 Simulated gas-driven eruptions using CO₂–water-polymer systems are reported. Eruptions are initiated by rapidly decompressing CO₂–saturated water containing up to 1.0 wt.% CO₂. Both cylindrical test cells and a flask test cell were used to examine the effect of magma chamber/conduit geometry on eruption dynamics. Bubble-growth kinetics are examined quantitatively in experiments using cylindrical test cells. Uninhibited bubble growth can be roughly expressed as dr/dt≈λD(β-1)/(γt ^1/3) for a CO₂–water-polymer system at 0–22  °C and with viscosities up to 5 Pa·s, where r is the radius of bubbles, λ and D are the Ostwald solubility coefficient and diffusivity of the gas in the liquid, β is the degree of saturation (decompression ratio), and γ characterizes how the boundary layer thickness increases with time and is roughly 1.0×10^–5 m/s^1/3 in this system. Unlike the radius of cylindrical test cells, which does not affect the eruption threshold and dynamics, the shape of the test cells (flask vs cylindrical) affects the dynamics but not the threshold of eruptions. For cylindrical test cells, the front motion is characterized by constant acceleration with both Δh (the height increase) and ΔV (the volume increase) being proportional to t ²; for the flask test cell, however, neither Δh nor ΔV is proportional to t ² as the conduit radius varies. Test-cell geometry also affects foam stability. In the flask test cell, as it moves from the wider base chamber into the narrower conduit, the bubbly flow becomes fragmented, affecting the eruption dynamics. The fragmentation may be caused by a sudden increase in acceleration induced by conduit-shape change, or by the presence of obstacles to the bubbly flow. This result may help explain the range in vesicularities of pumice and reticulite. Received: 16 May 1997 / Accepted: 11 October 1997     

Estimation of the gravitational potential energy of the earth based on different density models

The estimation of the Earth’s gravitational potential energy E was obtained for different density distributions and rests on the expression E = − (W_min + ΔW) derived from the conventional relationship for E. The first component W_min expresses minimum amount of the work W and the second component ΔW represents a deviation from W_min interpreted in terms of Dirichlet’s integral applied on the internal potential. Relationships between the internal potential and E were developed for continuous and piecewise continuous density distributions. The global 3D density model inside an ellipsoid of revolution was chosen as a combined solution of the 3D continuous distribution and the reference PREM radial piecewise continuous profile. All the estimates of E were obtained for the spherical Earth since the estimated (from error propagation rule) accuracy σ_E of the energy E is at least two orders greater than the ellipsoidal reduction and the contribution of lateral density inhomogeneities of the 3D global density model. The energy E contained in the 2nd degree Stokes coefficients was determined. A good agreement between E = E_Gauss derived from Gaussian distribution and other E, in particular for E = E_PREM based on the PREM piecewise continuous density model and E-estimates derived from simplest Legendre-Laplace, Roche, Bullard and Gauss models separated into core and mantle only, suggests the Gaussian distribution as a basic radial model when information about density jumps is absent or incomplete.     

Wave disturbances in the ionosphere during a lasting solar activity minimum

Using the Kharkov incoherent scatter radar, observations of wave disturbances in electron concentration N in the ionosphere at heights of 120–600 km are conducted. The measurements were carried out in the periods of the spring and fall equinoxes and winter and summer solstices. The height-time dependences of the absolute ΔN and relative ΔN/N amplitudes of wave disturbances, as well as their spectral composition, were analyzed. It is shown that wave disturbances in the ionosphere with periods of 10–180 min were present at almost any time of the day and in all seasons. Their absolute and relative amplitudes varied from 6 × 10⁹ to 6 × 10¹⁰ m⁻³ and from 0.01 to 0.5, respectively. The maximum values of ΔN and ΔN/N were observed at a height of ∼200 km. The passage of the solar terminator changed substantially the wave disturbance parameters.     

The revised m Lg scale in six eastern provinces of China

The influence of the four different methods of measuring Lg amplitude, and the selection of different geometrical attenuation coefficient ζ-values (=5/6 or 1) on the determination ofγ value of Lg wave are discussed.γ=0.0034±0.0001km⁻¹ (when ζ=5/6) for six eastern provinces is redetermined. The revised magnitude calibration function ofm _Lg (mxh),q _E (Δ)=(5/6)logΔ+0.00147Δ+1.81 is deduced. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 171–178, 1991. Projects sponsored by the Chinese Joint Seismological Science Foundation.     

Geochemistry and Sr-Nd-Pb isotopic characteristics of the Mugouriwang Cenozoic volcanic rocks from Tibetan Plateau:Constraints on mantle source of the underplated basic magma

The Mugouriwang Cenozoic volcanic rocks exposed in the north Qiangtang Block of Tibetan Plateau are mainly composed of basalt and andesitic-basalt,both characterized by the lower SiO2 (51%―54%),high refractory elements (i.e. Mg,Cr,Ni) as well as the moderate enrichment in light rare earth elements (LREE) relative to a slight depleted in Eu and high strength field elements (HFSE,i.e. Nb,Ta,Ti). Be-sides,the fairly low Sm/Yb value (3.07―4.35) could signify that the rocks should be derived directly from partial melting of the spinel lherzolite at the upper part of the asthenosphere. These rocks have radiogenic Sr and Pb (87Sr/86Sr = 0.705339 to 0.705667; 208Pb/204Pb = 38.8192 to 38.8937; 207Pb/204Pb = 15.6093 to 15.6245; 206Pb/204Pb = 18.6246 to 18.6383),and non-radiogenic Nd (143Nd/144Nd = 0.512604 to 0.512639; εNd = 0.02 to -0.66) in agreement with those values of the BSE mantle reservoir. The DUPAL anomaly of the rocks can be evidently attested by the △8/4Pb = 66.82 to 74.53 ,△7/4Pb = 9.88 to 11.42,△Sr>50,implying that the Mugouriwang volcanic rock is likely to be generated by partial melting of a Gondwana-bearing asthenospheric mantle ever matasomatised by the fluid from subduction zone. Depending on the previous study on the high-K calc-alkaline intermediate-felsic volcanics in the study area,this paper proposed that the fluids derived from the subducted Lhasa Block metasomatised the asthenosphere beneath the Qiangtang Block,and induced its partial melting,and then the melt under-plated the thickened Qiangtang lithosphere and caused the generation of the Cenozoic adakite-like felsic magmas in the Qiangtang region.     

He,Ne and Ar isotopic composition of Fe-Mn crusts from the western and central Pacific Ocean and implications for their genesis

The noble gas nuclide abundances and isotopic ratios of the upmost layer of Fe-Mn crusts from the western and central Pacific Ocean have been determined. The results indicate that the He and Ar nu- clide abundances and isotopic ratios can be classified into two types: low 3He/4He type and high 3He/4He type. The low 3He/4He type is characterized by high 4He abundances of 191×10-9 cm3·STP·g-1 on average, with variable 4He, 20Ne and 40Ar abundances in the range (42.8―421)×10-9 cm3·STP·g-1, (5.40―141)×10-9 cm3·STP·g-1, and (773―10976)×10-9 cm3·STP·g-1, respectively. The high 3He/4He samples are characterized by low 4He abundances of 11.7×10-9 cm3·STP·g-1 on average, with 4He, 20Ne and 40Ar abundances in the range of (7.57―17.4)×10-9 cm3·STP·g-1, (10.4―25.5)×10-9 cm3·STP·g-1 and (5354―9050)×10-9 cm3·STP·g-1, respectively. The low 3He/4He samples have 3He/4He ratios (with R/RA ratios of 2.04―2.92) which are lower than those of MORB (R/RA=8±1) and 40Ar/36Ar ratios (447―543) which are higher than those of air (295.5). The high 3He/4He samples have 3He/4He ratios (with R/RA ratios of 10.4―12.0) slightly higher than those of MORB (R/RA=8±1) and 40Ar/36Ar ratios (293―299) very similar to those of air (295.5). The Ne isotopic ratios (20Ne/22Ne and 21Ne/22Ne ratios of 10.3―10.9 and 0.02774―0.03039, respectively) and the 38Ar/36Ar ratios (0.1886―0.1963) have narrow ranges which are very similar to those of air (the 20Ne/22Ne, 21Ne/22Ne, 38Ar/36Ar ratios of 9.80, 0.029 and 0.187, respectively), and cannot be differentiated into different groups. The noble gas nuclide abundances and isotopic ratios, together with their regional variability, suggest that the noble gases in the Fe-Mn crusts originate primarily from the lower mantle. The low 3He/4He type and high 3He/4He type samples have noble gas characteristics similar to those of HIMU (High U/Pb Mantle)- and EM (Enriched Mantle)-type mantle material, respectively. The low 3He/4He type samples with HIMU-type noble gas isotopic ratios occur in the Magellan Seamounts, Marcus-Wake Seamounts, Marshall Island Chain and the Mid-Pacific Sea- mounts whereas the high 3He/4He type samples with EM-type noble gas isotopic ratios occur in the Line Island Chain. This difference in noble gas characteristics of these crust types implies that the MagellanSeamounts, Marcus-Wake Seamounts, Marshall Is- land Chain, and the Mid-Pacific Seamounts originated from HIMU-type lower mantle material whereas the Line Island Chain originated from EM-type lower mantle material. This finding is consistent with varia- tions in the Pb-isotope and trace element signatures in the seamount lavas. Differences in the mantlesource may therefore be responsible for variations in the noble gas abundances and isotopic ratios in the Fe-Mn crusts. Mantle degassing appears to be the principal factor controlling noble gas isotopic abundances in Fe-Mn crusts. Decay of radioactive isotopes has a negligible influence on the nuclide abundances and isotopic ratios of noble gases in these crusts on the timescale of their formation.     

High Frequency PKKP<Subscript>BC</Subscript> Around 2.5 Hz Recorded Globally

Anomalous high frequency PKKP_BC signals (displaying a large amount of energy around 2.5 Hz), recorded globally for deep and intermediate depth earthquakes, are compared to PKKP_AB signals. The attenuation difference t_\textAB^* - t_\textBC^* t_{\text{AB}}^{*} - t_{\text{BC}}^{*} is evaluated from spectral amplitudes in the range 96–111°, being approximately twice the results provided by full-wave theory and PREM (with no low Q_μ zone in the lowermost mantle and a nearly infinite Q_K in the outer core). Most ray paths for such recordings are piercing the D″ region in the proximity of regions where ultra-low velocity zones (ULVZ) have been previously reported beneath the North Atlantic Ocean, the Southwest Pacific and the southwestern part of South America. If BC amplitudes around 2.5 Hz and at low frequencies (0.5–1.5 Hz) are comparable, the observed attenuation difference (in the frequency range 0.2–2.5 Hz) is small (around 0.25 s) and close to the PREM value. The particle motion of the high-frequency PKKP_BC at 2.5 Hz is quite similar to that of the raw recording, suggesting a deep source. An explanation for this might be scattering of the BC branch in some very restricted areas of the lowermost mantle. Alternately, the presence of a thin layer with high attenuation in the D″ region would most likely be associated with either the ultra-low velocity zone (ULVZ) or light sediments on the underside of the core-mantle boundary (CMB). Correlated to other methods to investigate the lowermost mantle, the high-frequency PKKP_BC can be used to map lateral variations of attenuation above the CMB, possibly associated with the boundary of the superplumes, especially when PKKP_AB is observed.     

Model of vibrational level populations of Herzberg states of oxygen molecules at heights of the lower thermosphere and mesosphere

The paper presents a model of the kinetics of electronically excited O₂(c¹Σ_u⁻,v), O₂(A′³Δ_u,v), O₂(A³Σ_u⁺,v) molecules at heights of the lower thermosphere and mesosphere with allowance for electronic excitation transfer processes during molecular collisions. The model is used to calculate the relative O₂(A³Σ_u⁺,v) and O₂(A′³Δ_u,v) populations at heights of 80–110 km. The calculated populations are compared with the available literature results on experimental estimates, and good agreement is obtained. It is shown how the increase in the quenching rates of the considered states by oxygen atoms affects the calculation results.