Testing the reference Moon model in respect of the thermal regime and chemical composition of the mantle: Thermodynamics versus seismology

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, Al₂O₃, 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 V_S) 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.     

Influence of Solar Wind Plasma Parameters on the Intensity of Isolated Magnetospheric Substorms

Parameters of the interplanetary magnetic field and solar wind plasma during periods of 163 isolated substorms have been studied. It is shown that the solar wind velocity V and plasma density N remain approximately constant for at least 3 h before substorm onset Т _o and 1 h after Т _o . On average, the velocity of the solar wind exhibits a stable trend toward anticorrelation with its density over the whole data array. However, the situation is different if the values of V and N are considered with respect to the intensity of substorms observed during that period. With the growth of substorm intensity, quantified as the maximum absolute value of AL index, an increase in both the solar wind plasma velocity and density, at which these substorms appear, is obsreved. It has been found that the magnitude of the solar wind dynamic pressure P is closely related to the magnetosphere energy load defined as averaged values of the Kan–Lee electric field E_KL and Newell parameter dΦ/dt averaged for 1 h interval before Т _o . The growth of the dynamic pressure is accompanied by an increase in the load energy necessary for substorm generation. This interrelation between P and values of EKL and dΦ/dt is absent in other, arbitrarily chosen periods. It is believed that the processes accompanying increasing dynamic pressure of the solar wind result in the formation of magnetosphere conditions that increasingly impede substorm generation. Thus, the larger is P, the more solar wind energy must enter the Earth’s magnetosphere during the period of the growth phase for substorm generation. This energy is later released during the period of the substorm expansion phase and creates even more intense magnetic bays.     

Homologous temperature of olivine: Implications for creep of the upper mantle and fabric transitions in olivine

The homologues temperature of a crystalline material is defined as T/T_m, where T is temperature and T_m is the melting(solidus) temperature in Kelvin. It has been widely used to compare the creep strength of crystalline materials. The melting temperature of olivine system,(Mg,Fe)_2SiO_4, decreases with increasing iron content and water content, and increases with confining pressure. At high pressure, phase transition will lead to a sharp change in the melting curve of olivine. After calibrating previous melting experiments on fayalite(Fe_2SiO_4), the triple point of fayalite-Fe_2SiO_4 spinel-liquid is determined to be at 6.4 GPa and 1793 K. Using the generalized means, the solidus and liquidus of dry olivine are described as a function of iron content and pressure up to 6.4 GPa. The change of T/T_m of olivine with depth allows us to compare the strength of the upper mantle with different thermal states and olivine composition. The transition from semi-brittle to ductile deformation in the upper mantle occurs at a depth where T/T_m of olivine equals 0.5. The lithospheric mantle beneath cratons shows much smaller T/T_m of olivine than orogens and extensional basins until the lithosphere-asthenosphere boundary where T/T_m 0.66, suggesting a stronger lithosphere beneath cratons. In addition, T/T_m is used to analyze deformation experiments on olivine. The results indicate that the effect of water on fabric transitions in olivine is closely related with pressure. The hydrogen-weakening effect and its relationship with T/T_m of olivine need further investigation. Below 6.4 GPa(200 km), T/T_m of olivine controls the transition of dislocation glide from [100] slip to [001] slip. Under the strain rate of 10~(-12)–10~(-15) s~(-1) and low stress in the upper mantle, the [100](010) slip system(A-type fabric) becomes dominant when T/T_m 0.55–0.60. When T/T_m 0.55–0.60, [001] slip is easier and low T/T_m favors the operation of [001](100) slip system(C-type fabric). This is consistent with the widely observed A-type olivine fabric in naturally deformed peridotites, and the C-type olivine fabric in peridotites that experienced deep subduction in ultrahigh-pressure metamorphic terranes. However, the B-type fabric will develop under high stress and relatively low T/T_m. Therefore, the homologues temperature of olivine established a bridge to extrapolate deformation experiments to rheology of the upper mantle. Seismic anisotropy of the upper mantle beneath cratons should be simulated using a four-layer model with the relic A-type fabric in the upper lithospheric mantle, the B-type fabric in the middle layer, the newly formed A- or B-type fabric near the lithosphere-asthenosphere boundary, and the asthenosphere dominated by diffusion creep below the Lehmann discontinuity. Knowledge about transition mechanisms of olivine fabrics is critical for tracing the water distribution and mantle flow from seismic anisotropy.     

Negative pressure effect on the electrical conductivity of San Carlos olivine and its implication to the electrical structure in the upper mantle

Pressure effect on the electrical conductivity of San Carlos olivine was investigated by the newly installed electrical conductivity measurement system at China University of Geosciences. Electrical conductivity of San Carlos olivine aggregates was measured up to 12 GPa and 1475 K using the Walker-type multi-anvil apparatus equipped with eight WC cubes as the second-stage anvils. The pressure generation against applied load for the experimental assemblage was examined by phase transition of Bi,quartz, forsterite under different P-T conditions. To check the data validity of this new system, electrical conductivities of the serpentinites and talc samples were measured. The results are consistent with the published data of the same samples. Electrical conductivity(σ) of the San Carlos olivine aggregates and temperature(T) satisfy the Arrhenian formula: σ=σ0exp[.(ΔE+PΔV)/kT].The pre-exponential factor(σ0), activation energy(ΔE) and activation volume(ΔV) yield value of 7.74 S/m, 0.85 eV and 0.94cm3/mol, respectively. Electrical conductivities of the San Carlos olivine aggregates decline with increasing pressure at same temperatures. The negative pressure effect can be interpreted by strain energy model of defect energy together with the lattice deformation. In addition, the electrical conductivity-depth 1-D profile of the upper mantle was constructed based on our results and some assumptions. The calculated profile is concordant with the geophysical observation at the depth of 180–350 km beneath Europe, which indicates that the upper mantle beneath Europe might be dry.     

Characteristics of Seismicity in the Areas of Large Water Reservoirs and Waterfalls: The Role of Effects from Additional Load and Permanent Vibration

The characteristics of seismicity in the near vicinity of five large water reservoirs and three large waterfalls from different regions of the Earth are considered. It is found that in some cases induced seismicity manifests itself during the filling of reservoirs at quite large depths: in the lower crust and even in the upper mantle. There is negative correlation between the maximum magnitudes Мmax of the earthquakes recorded near water reservoirs and waterfalls and the water discharge in these objects (V _p ). The largest values of М_max are characteristic of earthquakes that occurred near Sarez Lake (Tajikistan) and the Koyna Reservoir (India), which have the lowest V _p ; in contrast, the smallest magnitudes are reported for earthquakes in the areas of the Khone Falls (Laos) and Niagara Falls (United States, Canada), where there are no large artificial water reservoirs, but huge water discharge takes place. The available data indicate that permanent vibration caused by falling water reduces the level of seismicity.     

Intersector disbalance of the large-scale open solar magnetic field fluxes,active and passive boundaries,and solar-terrestrial extrastorms

The dynamics (from rotation to rotation) of the absolute values of the large-scale open solar magnetic field fluxes in the four-sector field structure has been considered for the first time, using CRs 2032–2035 in July–October 2005 as examples. An important role of the ratio of the fluxes at the eastern and western sector boundaries (Φ _E /Φ _W ) is confirmed. As in the cases of the two-sector structure, Φ _E /Φ _W > 1 is typical of active rigidly corotating boundaries with intense sunspot formation, flares, and interplanetary and geomagnetic disturbances. A remarkable property of the considered structure was the presence of a rapidly increasing flux in an initially narrow sector and the flux interaction with a stable rigidly corotating sector in the zone of the main active longitudes, which caused an unexpectedly strong geoeffective long-range action of flares near the corresponding active boundary.     

Limits of transversely isotropic elastic parameters for the existence of classical Rayleigh waves

Solutions of P-SV equations of motion in a homogeneous transversely isotropic elastic layer contain a factor exp(±ν _j z), where z is the vertical coordinate and j?=?1, 2. For computing Rayleigh wave dispersion in a multi-layered half space, ν _j is computed at each layer. For a given phase velocity (c), ν _j becomes complex depending on the transversely isotropic parameters. When ν _j is complex, classical Rayleigh waves do not exist and generalised Rayleigh waves propagate along a path inclined to the interface. We use transversely isotropic parameters as α _H , β _V , ξ, ? and η and find their limits beyond which ν _j becomes complex. It is seen that ν _j depends on ? and η, but does not depend on ξ. The complex ν _j occurs when ? is small and η is large. For a given c/β _V , the region of complex ν _j in a ? -η plane increases with the increase of α _H /β _V . Further, for a given α _H /β _V , the complex region of ν _j increases significantly with the decrease of c/β _V . This study is useful to compute dispersion parameters of Rayleigh waves in a layered medium.     

Crustal Structure Along Sunda-Banda Arc Transition Zone from Teleseismic Receiver Functions

We analyzed receiver function of teleseismic events recorded at twelve Indonesian-GEOFON (IA-GE) broadband stations using nonlinear Neighbourhood Algorithm (NA) inversion and H-k stacking methods to estimate crustal thickness, V _p /V _s ratios and S-wave velocity structure along Sunda-Banda arc transition zone. We observed crustal thickness of 34–37 km in Timor Island, which is consistent with the previous works. The thick crust (> 30 km) is also found beneath Sumba and Flores Islands, which might be related to the arc-continent collision causing the thickened crust. In Timor and Sumba Islands, we observed high V _p /V _s ratio (> 1.84) with low velocity zone that might be associated with the presence of mafic and ultramafic materials and fluid filled fracture zone. The high V _p /V _s ratio observed at Sumbawa and Flores volcanic Islands might be an indication of partial melt related to the upwelling of hot asthenosphere material through the subducted slab.     

Site Specific Ground Motion Modeling and Seismic Response Analysis for Microzonation of Baku,Azerbaijan

We investigated ground response for Baku (Azerbaijan) from two earthquakes of magnitude M6.3 occurred in Caspian Sea (characterized as a near event) and M7.5 in Shamakhi (characterized as a remote extreme event). S-wave velocity with the average shear wave velocity over the topmost 30 m of soil is obtained by experimental method from the V _P values measured for the soils. The downtown part of Baku city is characterized by low V_S30 values (< 250 m/s), related to sand, water-saturated sand, gravel-pebble, and limestone with clay. High surface PGA of 240 gal for the M7.5 event and of about 190 gal for the M6.3 event, and hence a high ground motion amplification, is observed in the shoreline area, through downtown, in the north-west, and in the east parts of Baku city with soft clays, loamy sands, gravel, sediments.     

Shear wave velocity-based evaluation and design of stone column improved ground for liquefaction mitigation

The evaluation and design of stone column improvement ground for liquefaction mitigation is a challenging issue for the state of practice. In this paper, a shear wave velocity-based approach is proposed based on the well-defined correlations of liquefaction resistance (CRR)-shear wave velocity (V _s)-void ratio (e) of sandy soils, and the values of parameters in this approach are recommended for preliminary design purpose when site specific values are not available. The detailed procedures of pre- and post-improvement liquefaction evaluations and stone column design are given. According to this approach, the required level of ground improvement will be met once the target V _s of soil is raised high enough (i.e., no less than the critical velocity) to resist the given earthquake loading according to the CRR-V _s relationship, and then this requirement is transferred to the control of target void ratio (i.e., the critical e) according to the V _s-e relationship. As this approach relies on the densification of the surrounding soil instead of the whole improved ground and is conservative by nature, specific considerations of the densification mechanism and effect are given, and the effects of drainage and reinforcement of stone columns are also discussed. A case study of a thermal power plant in Indonesia is introduced, where the effectiveness of stone column improved ground was evaluated by the proposed V _s-based method and compared with the SPT-based evaluation. This improved ground performed well and experienced no liquefaction during subsequent strong earthquakes.     

Mantle upwelling beneath Madagascar: evidence from receiver function analysis and shear wave splitting

Crustal receiver functions have been calculated from 128 events for two three-component broadband seismomenters located on the south coast (FOMA) and in the central High Plateaux (ABPO) of Madagascar. For each station, crustal thickness and V _p /V _s ratio were estimated from H- κ plots. Self-consistent receiver functions from a smaller back-azimuthal range were then selected, stacked and inverted to determine shear wave velocity structure as a function of depth. These results were corroborated by guided forward modeling and by Monte Carlo error analysis. The crust is found to be thinner (39 ± 0.7 km) beneath the highland center of Madagascar compared to the coast (44 ± 1.6 km), which is the opposite of what would be expected for crustal isostasy, suggesting that present-day long wavelength topography is maintained, at least in part, dynamically. This inference of dynamic support is corroborated by shear wave splitting analyses at the same stations, which produce an overwhelming majority of null results (>96 %), as expected for vertical mantle flow or asthenospheric upwelling beneath the island. These findings suggest a sub-plate origin for dynamic support.     

Generation of superDreicer electric fields in the solar chromosphere

The electric field generation at the front of the current pulse, which originates in a coronal magnetic loop owing to the development of the Rayleigh–Taylor magnetic instability at loop footpoints, has been considered. During the τ_A ≈ l/V _A ≈ 5?25 s time (where l is the plasma plume height entering a magnetic loop as a result of the Rayleigh–Taylor instability), a disturbance related to the magnetic field tension B _?(r,t), “escapes” the instability region with the Alfvén velocity in this case. As a result, an electric current pulse I_z(z ? V _A t), at the front of which an induction magnetic field E _z, which is directed along the magnetic tube axis and can therefore accelerate particles, starts propagating along a magnetic loop with a characteristic scale of Δξ ≈ l. In the case of sufficiently large currents, when B _? ²/8π > p, an electric current pulse propagates nonlinearly, and a relatively large longitudinal electric field originates E _z ≈ 2I _z ³ V _A/c ⁴a²B_z ²l, which can be larger than the Dreicer field, depending on the electric current value.     

The noble gas isotope geochemical composition of chert at the bottom of Cambrian in Tarim Basin,China

In the Tarim Basin, black shale series at the bottom of Cambrian is one of the important marine facies hydrocarbon source rocks. This research focuses on the analysis of the isotope of noble gas of 11 cherts. The R/R _a ratio of chert in the Keping area is 0.032–0.319, and ⁴⁰Ar/³⁶Ar is 338–430. In Quruqtagh the R/R _a ratio is 0.44–10.21, and ⁴⁰Ar/³⁶Ar is 360–765. The R/R _a ratio of chert increases with ⁴⁰Ar/³⁶Ar from the west to the east accordingly. They have evolved from the crust source area to the mantle source area in a direct proportion. Surplus argon ⁴⁰Ar_E in chert is in direct proportion to the R/R _a ratio, indicating that it has the same origin of excess argon as in fluid and mantle source helium. Comparison of the R/R _a ratios between the west and the east shows that the chert in the eastern part formed from the activity system of the bottom hydrothermal venting driven by the mantle source, where the material and energy of crust and mantle had a strong interaction in exchange; whereas in the western part, chert deposited from the floating of hydrothermal plume undersea bottom, which is far away from the centre of activities of the hydrothermal fluid of ocean bottom. In addition, from noble gas isotope composition of chert, it is suggested that the ocean anoxia incident happened at the black shale of the Cambrian bottom probably because of the large-scaled ocean volcanoes and the following hydrothermal activities.     

DHI evaluation by combining rock physics simulation and statistical techniques for fluid identification of Cambrian-to-Cretaceous clastic reservoirs in Pakistan

The use of seismic direct hydrocarbon indicators is very common in exploration and reservoir development to minimise exploration risk and to optimise the location of production wells. DHIs can be enhanced using AVO methods to calculate seismic attributes that approximate relative elastic properties. In this study, we analyse the sensitivity to pore fluid changes of a range of elastic properties by combining rock physics studies and statistical techniques and determine which provide the best basis for DHIs. Gassmann fluid substitution is applied to the well log data and various elastic properties are evaluated by measuring the degree of separation that they achieve between gas sands and wet sands. The method has been applied successfully to well log data from proven reservoirs in three different siliciclastic environments of Cambrian, Jurassic, and Cretaceous ages. We have quantified the sensitivity of various elastic properties such as acoustic and extended elastic (EEI) impedances, elastic moduli (K _sat and K _sat–μ), lambda–mu–rho method (λρ and μρ), P-to-S-wave velocity ratio (V _P/V _S), and Poisson’s ratio (σ) at fully gas/water saturation scenarios. The results are strongly dependent on the local geological settings and our modeling demonstrates that for Cambrian and Cretaceous reservoirs, K _sat–μ, EEI, V _P/V _S, and σ are more sensitive to pore fluids (gas/water). For the Jurassic reservoir, the sensitivity of all elastic and seismic properties to pore fluid reduces due to high overburden pressure and the resultant low porosity. Fluid indicators are evaluated using two metrics: a fluid indicator coefficient based on a Gaussian model and an overlap coefficient which makes no assumptions about a distribution model. This study will provide a potential way to identify gas sand zones in future exploration.     

Composition,temperature, and thickness of the lithosphere of the Archean Kaapvaal craton

A new model is proposed for the structure of the Kaapvaal craton lithosphere. Based on chemical thermodynamics methods, profiles of the chemical composition, temperature, density, and S wave velocities are constructed for depths of 100–300 km. A solid-state zone of lower velocities is discovered on the S velocity profile in the depth interval 150–260 km. The temperature profiles are obtained from absolute values of P and S velocities, taking into account phase transformations, anharmonicity, and anelastic effects. The examination of the sensitivity of seismic models to the chemical composition showed that relatively small variations in the composition of South African xenoliths result in lateral temperature variations of ～200°C. Inversion of some seismic profiles (including IASP91) with a fixed bulk composition of garnet peridotites (the primitive mantle material) leads to a temperature inversion at depths of 200–250 km, which is physically meaningless. It is supposed that the temperature inversion can be removed by gradual fertilization of the mantle with depth. In this case, the craton lithosphere should be stratified in chemical composition. The depleted lithosphere composed by garnet peridotites exists to depths of 175–200 km. The lithospheric material at depths of 200–250 km is enriched in basaltoid components (FeO, Al₂O₃, and CaO) as compared with the material of garnet peridotites but is depleted in the same components as compared with the fertile substance of the underlying primitive mantle. The material composing the craton root at a depth of ～275 km does not differ in its physical and chemical characteristics from the composition of the normal mantle, and this allows one to estimate the thickness of the lithosphere at 275 km. The results of this work are compared with data of seismology, thermal investigations, and thermobarometry.     

Shear Wave Velocity Estimates through Combined Use of Passive Techniques in a Tectonically Active Area

We made an attempt to assess the shear wave velocity values V _S and, to a lesser extent, the V _P values from ambient noise recordings in an array configuration. Five array sites were situated in the close proximity to borehole sites. Shear wave velocity profiles were modeled at these five array sites with the aid of two computational techniques, viz. spatial autocorrelation (SPAC) and H/V ellipticity. Out of these five array sites, velocity estimates could be reliably inferred at three locations. The shear wave velocities estimated by these methods are found to be quite consistent with each other. The computed V _S values up to 30 m depth are in the range from 275 to 375 m/s in most of the sites, which implies prevalence of a low velocity zone at some pocket areas. The results were corroborated by evidence of site geology as well as geotechnical information.     

Determination of source parameters of strong earthquakes of Kamchatka,the Kurile Islands,and Japan from aftershock sequences

Aftershock sequences of some strong earthquakes of Kamchatka, the Kurile Islands, and Japan are examined. Such source parameters as the length L, along-dip width W, motion on fault D, and stress drop Δσ are determined from the aftershock sequences considered. The values of these parameters were obtained by the formal estimation of linear source parameters (lower bound estimates) and visually (upper bound estimates). The correlation dependences of the obtained parameters on the surface wave (M _S ) and seismic moment (M _W ) magnitudes are calculated.     

Influence of aquatic plants on the hydrogen isotope composition of sedimentary long-chain n-alkanes in the Lake Qinghai region,Qinghai-Tibet Plateau

The hydrogen isotopic composition(δD) of leaf wax long-chain n-alkanes(C27, C29, and C31) from lacustrine sediments has been widely applied to reconstruct terrestrial paleoclimatic and paleohydrological changes. However, few studies have addressed whether the aquatic-derived n-alkanes can affect the δD values of lake sedimentary long-chain n-alkanes, which are usually regarded as a recorder of the terrestrial hydrological signals. Here we systematically investigated δD values of long-chain n-alkanes from modern aquatic plants, both near-shore and off-shore surface sediments, surrounding terrestrial plant litters, as well as river water and lake water in Lake Qinghai and its satellite lakes on the northeastern Qinghai-Tibet Plateau. Our data showed that(i) δD values of long-chain n-alkanes from aquatic plants varied from-184‰ to-132‰ for n-C27, from-183‰ to-138‰ for n-C29, and from-189‰ to-130‰ for n-C31, respectively, with no significant differences among the three n-alkanes homologues;(ii) δD values of long-chain n-alkanes from aquatic plants were generally more positive than those from surrounding terrestrial plants, possibly because that they recorded the D-enrichment of lake water in this semi-arid region;(iii) δD values of long-chain n-alkanes from surface sediments showed significant differences among the three n-alkanes homologues, due to the larger aquatic input of n-C27 to the sedimentary lipid pool than that of n-C31, and(iv) n-C27 δD values of near-shore aquatic plants and near-shore sediments are more negative than those from off-shore as a result of lower δD values of near-shore lake water. Our findings indicate that in this region(i) the offset between sedimentary n-C27 and n-C31 δD values(ΔδDC27-C31) could potentially be used to evaluate if sedimentary long-chain n-alkanes are derived from a single source;(ii) while δD values of n-C27 may be influenced by lake water hydrological changes, sedimentary n-C31 is derived predominantly from terrestrial plants and thus its δD can serve as a relatively reliable indicator for terrestrial paleoclimatic and paleohydrological reconstructions.     

Frictional and viscous effects on the nucleation phase of an earthquake

We study the frictional and viscous effects on earthquake nucleation, especially for the nucleation phase, based on a one-degree-of-freedom spring-slider model with friction and viscosity. The frictional and viscous effects are specified by the characteristic displacement, U _c, and viscosity coefficient, η, respectively. Simulation results show that friction and viscosity can both lengthen the natural period of the system and viscosity increases the duration time of motion of the slider. Higher viscosity causes a smaller amplitude of lower velocity motion than lower viscosity. A change of either U _c (under large η) or η (under large U _c) from a large value (U _ch for U _c and η _h for η) to a small one (U _cl for U _c and η _l for η) in two stages during sliding can result in a clear nucleation phase prior to the P-wave. The differences δU _c = U _ch ? U _cl and δη = η _h ? η _l are two important factors in producing a nucleation phase. The difference between the nucleation phase and the P-wave increases with either δU _c or δη. Like seismic observations, the peak amplitude of P-wave, which is associated with the earthquake magnitude, is independent upon the duration time of nucleation phase. A mechanism specified with a change of either η or U _c from a larger value to a smaller one due to temporal variations in pore fluid pressure and temperature in the fault zone based on radiation efficiency is proposed to explain the simulation results and observations.     

Scaling relations of moment magnitude,local magnitude,and duration magnitude for earthquakes originated in northeast India

In this study, we aim to improve the scaling between the moment magnitude (M _W), local magnitude (M _L), and the duration magnitude (M _D) for 162 earthquakes in Shillong-Mikir plateau and its adjoining region of northeast India by extending the M _W estimates to lower magnitude earthquakes using spectral analysis of P-waves from vertical component seismograms. The M _W-M _L and M _W-M _D relationships are determined by linear regression analysis. It is found that, M _W values can be considered consistent with M _L and M _D, within 0.1 and 0.2 magnitude units respectively, in 90 % of the cases. The scaling relationships investigated comply well with similar relationships in other regions in the world and in other seismogenic areas in the northeast India region.