Estimating seismic moments and Lg <Emphasis Type="Italic">Q</Emphasis> using Lg spectra

A multi-event and multi-station inverse method is presented in the paper to simultaneously estimate the seismic moments (M ₀) and source corner frequencies (f _c) of several Jiashi (Xinjiang, China) earthquakes, as well as the apparent Lg Q models for the paths from Jiashi to eight seismic stations (WMQ, AAK, TLG, MAKZ, KUR, VOS, ZRN and CHK) in Central Asia. The resultant seismic moments correlate well with the M ₀ values obtained by Harvard University using the centroid moment tensor (CMT) inversion and the surface-wave magnitudes as well. After the correction by a typical value of average radiation coefficient for regional SV waves, the M ₀ values from Lg spectral inversion are still close to the corresponding values obtained from CMT inversion. The obtained apparent Q _0Lg values (Lg Q at 1 Hz) are consistent with the tectonic features of corresponding propagation paths. The Q _0Lg values are 351±87, 349±86 and 300±27 for the paths from Jiashi to AAK, TLG and MAKZ, respectively. They are smaller than Q _0Lg values for the paths to KUR, VOS, ZRN and CHK, which are 553±72, 569±58, 550±57 and 603±65, respectively. These results agree with the condition that the paths to AAK, TLG and MAKZ mainly propagate through the mountainous Tianshan area where relatively strong seismic activities and large variations of topography are exhibited, while the paths to KUR, VOS, ZRN and CHK mainly propagate through the stable area of Kazak platform. The Q _0Lg value for the path to WMQ is 462±56. This is also in agreement with the condition that the path to WMQ is basically along the border area between Tianshan Mountain and Tarim Basin, and along this path the variations of topography and crustal thickness are moderate in comparison with that along the path to MAKZ.     

A study of Guptkashi,Uttarakhand earthquake of 6 February 2017 (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript>5.3) in the Himalayan arc and implications for ground motion estimation

The 2017 Guptkashi earthquake occurred in a segment of the Himalayan arc with high potential for a strong earthquake in the near future. In this context, a careful analysis of the earthquake is important as it may shed light on source and ground motion characteristics during future earthquakes. Using the earthquake recording on a single broadband strong-motion seismograph installed at the epicenter, we estimate the earthquake’s location (30.546° N, 79.063° E), depth (H?=?19 km), the seismic moment (M₀?=?1.12×10¹⁷ Nm, M _w 5.3), the focal mechanism (φ?=?280°, δ?=?14°, λ?=?84°), the source radius (a?=?1.3 km), and the static stress drop (Δσ _s ~22 MPa). The event occurred just above the Main Himalayan Thrust. S-wave spectra of the earthquake at hard sites in the arc are well approximated (assuming ω^?2 source model) by attenuation parameters Q(f)?=?500f^0.9, κ?=?0.04 s, and f_max?=?infinite, and a stress drop of Δσ?=?70 MPa. Observed and computed peak ground motions, using stochastic method along with parameters inferred from spectral analysis, agree well with each other. These attenuation parameters are also reasonable for the observed spectra and/or peak ground motion parameters in the arc at distances ≤?200 km during five other earthquakes in the region (4.6?≤?M _w ?≤?6.9). The estimated stress drop of the six events ranges from 20 to 120 MPa. Our analysis suggests that attenuation parameters given above may be used for ground motion estimation at hard sites in the Himalayan arc via the stochastic method.     

Mechanism of formation of <Emphasis Type="Italic">Q</Emphasis>-disturbances in the <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> region of the equatorial ionosphere

Using model simulations, the morphological picture (revealed earlier) of the disturbances in the F ₂ region of the equatorial ionosphere under quiet geomagnetic conditions (Q-disturbances) is interpreted. It is shown that the observed variations in the velocity of the vertical E × B plasma drift, related to the zonal E _y component of the electric field, are responsible for the formation of Q-disturbances. The plasma recombination at altitudes of the lower part of the F ₂ region and the dependence of the rate of this process on heliogeophysical conditions compose the mechanism of Q-disturbance formation at night. The daytime positive Q-disturbances are caused exclusively by a decrease in the upward E × B drift, and this type of disturbances could be related to the known phenomenon of counter electrojet. Possible causes of formation of the daytime negative Q-disturbances are discussed.     

Region-specific <Emphasis Type="Italic">S</Emphasis>-wave attenuation for earthquakes in northwestern Iran

In this study, continuous wavelet transform is applied to estimate the frequency-dependent quality factor of shear waves, Q _S , in northwestern Iran. The dataset used in this study includes velocigrams of more than 50 events with magnitudes between 4.0 and 6.5, which have occurred in the study area. The CWT-based method shows a high-resolution technique for the estimation of S-wave frequency-dependent attenuation. The quality factor values are determined in the form of a power law as Q _S (f)?=?(147?±?16)f ^0.71?±?0.02 and (126?±?12)f ^0.73?±?0.02 for vertical and horizontal components, respectively, where f is between 0.9 and 12 Hz. Furthermore, in order to verify the reliability of the suggested Q _S estimator method, an additional test is performed by using accelerograms of Ahar-Varzaghan dual earthquakes on August 11, 2012, of moment magnitudes 6.4 and 6.3 and their aftershocks. Results indicate that the estimated Q _S values from CWT-based method are not very sensitive to the numbers and types of waveforms used (velocity or acceleration).     

Attenuation coefficients of Rayleigh and <Emphasis Type="Italic">Lg</Emphasis> waves

Analysis of the frequency dependence of the attenuation coefficient leads to significant changes in interpretation of seismic attenuation data. Here, several published surface-wave attenuation studies are revisited from a uniform viewpoint of the temporal attenuation coefficient, denoted by χ. Theoretically, χ( f) is expected to be linear in frequency, with a generally non-zero intercept γ?=?χ(0) related to the variations of geometrical spreading, and slope dχ/df = π/Q _e caused by the effective attenuation of the medium. This phenomenological model allows a simple classification of χ( f) dependences as combinations of linear segments within several frequency bands. Such linear patterns are indeed observed for Rayleigh waves at 500–100-s and 100–10-s periods, and also for Lg from ~2 s to ~1.5 Hz. The Lg χ( f) branch overlaps with similar linear branches of body, Pn, and coda waves, which were described earlier and extend to ~100 Hz. For surface waves shorter than ~100 s, γ values recorded in areas of stable and active tectonics are separated by the levels of \(\gamma _{D} \approx 0.2 \times 10^{-3}\) s^???1 (for Rayleigh waves) and 8 ×10^???3 s^???1 (for Lg). The recently recognized discrepancy between the values of Q measured from long-period surface waves and normal-mode oscillations could also be explained by a slight positive bias in the geometrical spreading of surface waves. Similarly to the apparent χ, the corresponding linear variation with frequency is inferred for the intrinsic attenuation coefficient, χ _i , which combines the effects of geometrical spreading and dissipation within the medium. Frequency-dependent rheological or scattering Q is not required for explaining any of the attenuation observations considered in this study. The often-interpreted increase of Q with frequency may be apparent and caused by using the Q-based model of attenuation and following preferred Q( f) dependences while ignoring the true χ( f) trends within the individual frequency bands.     

Distribution of Lg coda <Emphasis Type="Italic">Q</Emphasis> in Xinjiang and its adjacent regions

In this study, we collected 1 156 broadband vertical components records at 22 digital seismic stations in Xinjiang region, Ürümqi station, and 7 stations in the adjacent regions during the period of 1999–2003. The records were firstly processed by the stacked spectral ratio method to obtain Q ₀ (Q at 1 Hz) and the frequency correlation factor η corresponding to each path. Based on the results, the distribution images of Q ₀ and η in 1°×1° grids for Xinjiang region were gained by the back-projection technique. The results indicate that Q ₀ is high (300–450) in the Tarim platform and marginal Siberian platform, while Q ₀ is low (150–250) in the southern regions as west Kunlun fold system and Songpan-Ganzi fold system. In the northern regions as Junggar fold system and Tianshan fold system, Q ₀ is also low (250–300) and η varies between 0.5 and 0.9.     

A new <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> estimation parameter for use in earthquake early warning systems

We propose a method that employs the squared displacement integral (ID2) to estimate earthquake magnitudes in real time for use in earthquake early warning (EEW) systems. Moreover, using τ _c and P _d for comparison, we establish formulas for estimating the moment magnitudes of these three parameters based on the selected aftershocks (4.0 ≤ M _s  ≤ 6.5) of the 2008 Wenchuan earthquake. In this comparison, the proposed ID2 method displays the highest accuracy. Furthermore, we investigate the applicability of the initial parameters to large earthquakes by estimating the magnitude of the Wenchuan M _s 8.0 mainshock using a 3-s time window. Although these three parameters all display problems with saturation, the proposed ID2 parameter is relatively accurate. The evolutionary estimation of ID2 as a function of the time window shows that the estimation equation established with ID2 ^Ref determined from the first 8-s of P wave data can be directly applicable to predicate the magnitudes of 8.0. Therefore, the proposed ID2 parameter provides a robust estimator of earthquake moment magnitudes and can be used for EEW purposes.     

Sampling uncertainties and source <Emphasis Type="Italic">b</Emphasis> likelihood for the Gutenberg-Richter <Emphasis Type="Italic">b</Emphasis> value from the Aki-Utsu method

The Aki-Utsu method of Gutenberg-Richter (G-R) b value estimation is often misapplied so that estimations not using the G-R histogram are often meaningless because they are not based on adequate samples. We propose a method to estimate the likelihood Pr(b?b _m , N, M ₁, M ₂) that an observed b _m estimate, based on a sample of N magnitudes within an [M ₁????≤?ΔM/2,?M ₂?+?ΔM/2) range, where ΔM?=?0.1 is the usual rounding applied to magnitudes, is due to a “true” source b value, b, and use these likelihoods to estimate source b ranges corresponding to various confidence levels. As an example of application of the method, we estimate the b values before and after the occurrence of a 7.4-magnitude earthquake in the Mexican subduction zone, and find a difference of 0.82 between them with 100% confidence that the b values are different.     

Dependence of geomagnetic disturbances on extreme values of the solar wind <Emphasis Type="Italic">E</Emphasis><Subscript><Emphasis Type="Italic">y</Emphasis></Subscript> component

The dependence of the zonal geomagnetic indices (AE, Ap, Kp, Kn, and Dst) on the solar wind parameters (the electric field E _y component, dynamic pressure P _d and IMF irregularity σB) has been studied for two types of events: magnetic clouds and high-speed streams. Based on the empirical relationships, it has been established that the AE, Ap, Kp, and Kn indices are directly proportional to the E _y value at E _y < 12 mV m^?1 and are inversely proportional to this value at E _y > 12 mV m^?1 for the first-type events. On the contrary, the dependence of Dst on E _y is monotonous nonlinear. A linear dependence of all geomagnetic indices on E _y is typical of the second-type events. It has been indicated that the specific features of geoeffectiveness of magnetic clouds and high-speed solar wind streams are caused by the dependence of the electric field potential across the polar cap on the electric field, solar wind dynamic pressure, and IMF fluctuations.     

The <Emphasis Type="Italic">Q</Emphasis>-factor estimates for the crust and upper mantle in the vicinity of Sochi and Anapa (North Caucasus)

The regularities in the radiation and propagation of seismic waves in the regions of the North Caucasus are analyzed for estimating the ground motion parameters during the probable future strong earthquakes. Based on the records of the regional earthquakes with magnitudes M_W ~ 3.9–5.6 within epicentral distances up to ~300 km obtained during the period of digital measurements at the Sochi and Anapa seismic stations, the Q-factors in the vicinities of these sites are estimated at ~55 f^0.9 and ~90f^0.7, respectively. The estimates were obtained by the coda normalization method developed by Aki, Rautian, and other authors. This method is based on the phenomenon of suppression of the earthquake (source) effects and local (site) responses by coda waves in the S-wave spectra. The obtained Q-factor estimates can be used for forecasting the ground shaking parameters for the future probable strong earthquakes in the North Caucasus in the vicinities of Sochi and Anapa.     

Variation of <Emphasis Type="Italic">Q</Emphasis> value before and after the 1999 Xiuyan,Liaoning Province, <Emphasis Type="Italic">M</Emphasis>=5.4 earthquake on the basis of analysis on attenuation dispersion of P waves

A method for determining medium quality factor is developed on the basis of analyzing the attenuation dispersion of the arrived first period P wave. In order to enhance signal to noise ratio, improve the resolution in measurement and reduce systematic error we applied the data resampling technique. The group velocity delay of P wave was derived by using an improved multi-filtering method. Based on a linear viscoelastic relaxation model we deduced the medium quality factor Q _m, and associated error with 95% confidence level. Applying the method to the seismic record of the Xiuyan M=5.4 earthquake sequences we obtained the following result: (1) High Q _m started to appear from Nov. 9, 1999. The events giving the deduced high Q _m value clustered in a region with their epicenter distances being between 32 and 46 km to the Yingkou station. This Q _m versus distance observation obviously deviates from the normal trend of Q _m linearly increasing with distance. (2) The average Q _m before the 29 Dec. 1999 M=5.4 earthquake is 460, while the average Q _m between the M=5.4 event and the 12 Jan. 2000 M=5.1 earthquake is 391, and the average Q _m after the M=5.1 event is 204.     

Magnitude <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> in metropolitan France

The recent seismicity catalogue of metropolitan France Sismicité Instrumentale de l’Hexagone (SI-Hex) covers the period 1962–2009. It is the outcome of a multipartner project conducted between 2010 and 2013. In this catalogue, moment magnitudes (M _w) are mainly determined from short-period velocimetric records, the same records as those used by the Laboratoire de Détection Géophysique (LDG) for issuing local magnitudes (M _L) since 1962. Two distinct procedures are used, whether M _L-LDG is larger or smaller than 4. For M _L-LDG >4, M _w is computed by fitting the coda-wave amplitude on the raw records. Station corrections and regional properties of coda-wave attenuation are taken into account in the computations. For M _L-LDG ≤4, M _w is converted from M _L-LDG through linear regression rules. In the smallest magnitude range M _L-LDG <3.1, special attention is paid to the non-unity slope of the relation between the local magnitudes and M _w. All M _w determined during the SI-Hex project is calibrated according to reference M _w of recent events. As for some small events, no M _L-LDG has been determined; local magnitudes issued by other French networks or LDG duration magnitude (M _D) are first converted into M _L-LDG before applying the conversion rules. This paper shows how the different sources of information and the different magnitude ranges are combined in order to determine an unbiased set of M _w for the whole 38,027 events of the catalogue.     

Tomographic inversion of near-surface <Emphasis Type="Italic">Q</Emphasis> factor by combining surface and cross-hole seismic surveys

The estimation of the quality factor Q plays a fundamental role in enhancing seismic resolution via absorption compensation in the near-surface layer. We present a new geometry that can be used to acquire field data by combining surface and cross-hole surveys to decrease the effect of geophone coupling on Q estimation. In this study, we drilled number of receiver holes around the source hole, each hole has different depth and each geophone is placed geophones into the bottom of each receiver hole to avoid the effect of geophone coupling with the borehole wall on Q estimation in conventional cross-hole seismic surveys. We also propose a novel tomographic inversion of the Q factor without the effect of the source signature, and examine its stability and reliability using synthetic data. We estimate the Q factors of the near-surface layer in two different frequency bands using field data acquired in the Dagang Oilfield. The results show that seismic absorption in the near-surface layer is much greater than that in the subsurface strata. Thus, it is of critical practical importance to enhance the seismic solution by compensating for near-surface absorption. In addition, we derive different Q factors from two frequency bands, which can be treated, to some extent, as evidence of a frequency-dependent Q.     

Seismotectonic studies of the pleistoseist area of the <Emphasis Type="Italic">M</Emphasis> <Subscript> <Emphasis Type="Italic">s</Emphasis> </Subscript> = 6.9 Ilin-Tass earthquake in northeast Yakutia

The complex seismotectonic studies of the pleistoseist area of the Ilin-Tas earthquake (M_s = 6.9), one of the strongest seismic events ever recorded by the regional seismic network in northeastern Russia, are carried out. The structural tectonic position, morphotectonic features of present-day topography, active faults, and types of Cenozoic deformations of the epicentral zone are analyzed. The data of the instrumental observations are summarized, and the manifestations of the strong seismic events in the Yana–Indigirka segment of the Cherskii seismotectonic zone are considered. The explanation is suggested for the dynamical tectonic setting responsible for the Andrei-Tas seismic maximum. This setting is created by the influence of the Kolyma–Omolon indenter, which intrudes into the Cherskii seismotectonic zone from the region of the North American lithospheric plate and forms the main seismogenic structures of the Yana–Indigirka segment in the frontal zone (the Ilin-Tas anticlinorium). The highest seismic potential is noted in the Andrei- Tas block—the focus of the main tectonic impacts from the Kolyma–Omolon superterrane. The general trend of this block coincides with the orientation of the major axis of isoseismal ellipses (azimuth 50°–85°), which were determined from the observations of macroseismic effects on the ground after the Uyandina (M_s = 5.6), Andrei-Tas (M_s = 6.1), and Ilin-Tas (M_s = 6.9) earthquakes.     

Differential physiological responses of the coastal cyanobacterium <Emphasis Type="Italic">Synechococcus</Emphasis> sp. PCC7002 to elevated <Emphasis Type="Italic">p</Emphasis>CO<Subscript>2</Subscript> at lag,exponential, and stationary growth phases

We studied the effects of expected end-of-the-century pCO₂ (1000 ppm) on the photosynthetic performance of a coastal marine cyanobacterium Synechococcus sp. PCC7002 during the lag, exponential, and stationary growth phases. Elevated pCO₂ significantly stimulated growth, and enhanced the maximum cell density during the stationary phase. Under ambient pCO₂ conditions, the lag phase lasted for 6 days, while elevated pCO₂ shortened the lag phase to two days and extended the exponential phase by four days. The elevated pCO₂ increased photosynthesis levels during the lag and exponential phases, but reduced them during the stationary phase. Moreover, the elevated pCO₂ reduced the saturated growth light (Ik) and increased the light utilization efficiency (α) during the exponential and stationary phases, and elevated the phycobilisome:chlorophyll a (Chl a) ratio. Furthermore, the elevated pCO₂ reduced the particulate organic carbon (POC):Chl a and particulate organic nitrogen (PON):Chl a ratios during the lag and stationary phases, but enhanced them during the exponential phase. Overall, Synechococcus showed differential physiological responses to elevated pCO₂ during different growth phases, thus providing insight into previous studies that focused on only the exponential phase, which may have biased the results relative to the effects of elevated pCO₂ in ecology or aquaculture.     

Precursor times of abnormal <Emphasis Type="Italic">b</Emphasis>-values prior to mainshocks

Seismic observations exhibit the presence of abnormal b-values prior to numerous earthquakes. The time interval from the appearance of abnormal b-values to the occurrence of mainshock is called the precursor time. There are two kinds of precursor times in use: the first one denoted by T is the time interval from the moment when the b-value starts to increase from the normal one to the abnormal one to the occurrence time of the forthcoming mainshock, and the second one denoted by T _p is the time interval from the moment when the abnormal b-value reaches the peak one to the occurrence time of the forthcoming mainshock. Let T* be the waiting time from the moment when the abnormal b-value returned to the normal one to the occurrence time of the forthcoming mainshock. The precursor time, T (usually in days), has been found to be related to the magnitude, M, of the mainshock expected in a linear form as log(T)?=?q?+?rM where q and r are the coefficient and slope, respectively. In this study, the values of T, T _p , and T* of 45 earthquakes with 3?≤?M?≤?9 occurred in various tectonic regions are compiled from or measured from the temporal variations in b-values given in numerous source materials. The relationships of T and T _p , respectively, versus M are inferred from compiled data. The difference between the values of T and T _p decreases with increasing M. In addition, the plots of T*/T versus M, T* versus T, and T* versus T-T* will be made and related equations between two quantities will be inferred from given data.     

Dynamics of the midlatitude ionospheric <Emphasis Type="Italic">F</Emphasis> region at sunrise

The behavior of the F2 layer at sunrise has been studied based on vertical-incidence ionospheric sounding data in Almaty (76°55′E, 43°15′N). Records with small amplitudes of electron density background fluctuations were selected in order to exactly estimate the onsets of a pronounced increase in the electron density at different altitudes. It has been indicated that the electron density growth rate is a function of altitude; in this case, the growth rate at the F2 layer maximum is much lower than such values at fixed altitudes of ～30–55 km below the layer maximum. The solar zenith angle (χ) and the blanketing layer thickness (h ₀) at the beginning of a pronounced increase in the electron density at altitude h are linearly related to the h value, and these quantities vary within ～90° < χ < 100° and 180 km < h ₀ < 260 km, respectively.     

Experimental examination on the heterogeneity parameter <Emphasis Type="Italic">C</Emphasis><Subscript>v</Subscript> of earthquake precursors

Two rock samples with different structures and materials were deformed under a biaxial loading system, and multipoint strain measurements were performed for each sample. The distribution of strain anomalies during the deformation and the instability process were analyzed by using C _v value put forward by WANG Xiao-qing and CHEN Xue-zhong, et al, a parameter to describe the heterogeneous distribution of earthquake precursors, so as to examine the method of C _v value and to explore its physical meaning experimentally. The result shows that the change of C _v value is correlated to the change of deformation characteristics and is an effective parameter to describe the heterogeneity of precursor distribution. C _v value increases firstly and then decreases before the instability, and the instability occurs when C _v value decreases to the level before increasing. This indicates that C _v value may be a useful parameter for earthquake prediction.     

Long-term changes in the relation between <Emphasis Type="Italic">foF</Emphasis>2 and <Emphasis Type="Italic">hmF</Emphasis>2

The change in the dependence of the F2-layer critical frequency on its height hmF2 is considered based on two sources of initial data used earlier by the authors. It is found that the slope k of the foF2 dependence on hmF2 systematically decreases from the earlier (“etalon”) period, 1958–1980, to the later periods of 1988–2010, 1998–2010, and 1998–2014. Since the foF2 value depends on the atomic oxygen concentration in the F region much more strongly than hmF2, the found decrease in k confirms the concept of a decrease in the atomic oxygen concentration in the thermosphere with time previously formulated by the authors.     

Construction of constant-<Emphasis Type="Italic">Q</Emphasis> viscoelastic model with three parameters

The popularly used viscoelastic models have some shortcomings in describing relationship between quality factor (Q) and frequency, which is not consistent with the observation data. Based on the theory of viscoelasticity, a new approach to construct constant-Q viscoelastic model in given frequency band with three parameters is developed. The designed model describes the frequency-independence feature of quality factor very well, and the effect of viscoelasticity on seismic wave field can be studied relatively accurate in theory with this model. Furthermore, the number of required parameters in this model has been reduced fewer than that of other constant-Q models, this can simplify the solution of the viscoelastic problems to some extent. At last, the accuracy and application range have been analyzed through numerical tests. The effect of viscoelasticity on wave propagation has been briefly illustrated through the change of frequency spectra and waveform in several different viscoelastic models.