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.     

A new method for testing pile by single-impact energy and <Emphasis Type="Italic">P-S</Emphasis> curve

By studying the pile-formula and stress-wave methods (e.g., CASE method), the authors propose a new method for testing piles using the single-impact energy and P-S curves. The vibration and wave figures are recorded, and the dynamic and static displacements are measured by different transducers near the top of piles when the pile is impacted by a heavy hammer or micro-rocket. By observing the transformation coefficient of driving energy (total energy), the consumed energy of wave motion and vibration and so on, the vertical bearing capacity for single pile is measured and calculated. Then, using the vibration wave diagram, the dynamic relation curves between the force (P) and the displacement (S) is calculated and the yield points are determined. Using the static-loading test, the dynamic results are checked and the relative constants of dynamic-static P-S curves are determined. Then the subsidence quantity corresponding to the bearing capacity is determined. Moreover, the shaped quality of the pile body can be judged from the formation of P-S curves.     

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.     

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.     

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.     

<Emphasis Type="Italic">R/S</Emphasis> analysis of the <Emphasis Type="Italic">Dst</Emphasis> index

The paper addresses estimation of the Hurst exponent for time series of the hourly values of the Dst index for the period from 1957 to 2011. It is found that the Hurst exponent is 0.79–0.94 for yearly intervals and 0.8–1.0 for monthly intervals. Based on R/S graphs, the Dst cycles are identified; they range from 3–4 months to 2.2 years and from 8.5 to 22 years in length. It is shown that a Dst time series can be quite satisfactorily described by an α-stable Levy process.     

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.     

A new method of testing pile using dynamic <Emphasis Type="Italic">P-S</Emphasis>-curve made by amplitude of wave train

A new method of detecting the vertical bearing capacity for single-pile with high strain is discussed in this paper. A heavy hammer or a small type of rocket is used to strike the pile top and the detectors are used to record vibration graphs. An expression of higher degree of strain (deformation force) is introduced. It is testified theoretically that the displacement, velocity and acceleration cannot be obtained by simple integral acceleration and differential velocity when long displacement and high strain exist, namely when the pile phase generates a whole slip relative to the soil body. That is to say that there are non-linear relations between them. It is educed accordingly that the force P and displacement S are calculated from the amplitude of wave train and (dynamic) P-S curve is drew so as to determine the yield points. Further, a method of determining the vertical bearing capacity for single-pile is discussed. A static load test is utilized to check the result of dynamic test and determine the correlative constants of dynamic-static P(Q)-S curve.     

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.     

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.     

<Emphasis Type="Italic">Q</Emphasis><Subscript>c</Subscript> and <Emphasis Type="Italic">Q</Emphasis><Subscript>S</Subscript> wave attenuation of South African earthquakes

Quality factor Q, which describes the attenuation of seismic waves with distance, was determined for South Africa using data recorded by the South African National Seismograph Network. Because of an objective paucity of seismicity in South Africa and modernisation of the seismograph network only in 2007, I carried out a coda wave decay analysis on only 13 tectonic earthquakes and 7 mine-related events for the magnitude range 3.6?≤?M _L ?≤?4.4. Up to five seismograph stations were utilised to determine Q _c for frequencies at 2, 4, 8 and 16 Hz resulting in 84 individual measurements. The constants Q ₀ and α were determined for the attenuation relation Q _c(f)?=?Q ₀ f ^α . The result was Q ₀?=?396?±?29 and α?=?0.72?±?0.04 for a lapse time of 1.9*(t _s???t ₀) (time from origin time t ₀ to the start of coda analysis window is 1.9 times the S-travel time, t _s) and a coda window length of 80 s. This lapse time and coda window length were found to fit the most individual frequencies for a signal-to-noise ratio of at least 3 and a minimum absolute correlation coefficient for the envelope of 0.5. For a positive correlation coefficient, the envelope amplitude increases with time and Q _c was not calculated. The derived Q _c was verified using the spectral ratio method on a smaller data set consisting of nine earthquakes and one mine-related event recorded by up to four seismograph stations. Since the spectral ratio method requires absolute amplitudes in its calculations, site response tests were performed to select four appropriate stations without soil amplification and/or signal distortion. The result obtained for Q _S was Q ₀?=?391?±?130 and α?=?0.60?±?0.16, which agrees well with the coda Q _c result.     

Studying the spread <Emphasis Type="Italic">F</Emphasis> effect before earthquakes

The observations of spread F during the nighttime hours (0000–0500 LT) have been statistically analyzed based on data of Tokyo, Akita, Wakkanai, and Yamagawa Japan vertical ionospheric sounding stations for the time intervals a month before and a month after an earthquake. The disturbances in the probability of spread F appearance before an earthquake are revealed against a background of the variations depending on season, solar activity cycle, geomagnetic and solar disturbances. The days with increased solar (Wolf number W > 100) and geomagnetic (ΣK > 30) activity are excluded from the analysis. The spread F effects are considered for more than a hundred earthquakes with magnitude M > 5 and epicenter depth h < 80 km at distances of R < 1000 km from epicenters to the vertical sounding station. An average decrease in the spread F occurrence probability one-two weeks before an earthquake has been revealed using the superposed epoch method (the probability was minimal approximately ten days before the event and then increased until the earthquake onset). Similar results are obtained for all four stations. The reliability of the effect has been estimated. The dependence of the detected effect on the magnitude and distance has been studied.     

Comparison of Bayesian estimates of peak ground acceleration (<Emphasis Type="Italic">A</Emphasis>max) with PSHA in Iran

Bayesian probability theory is an appropriate and useful method for estimating parameters in seismic hazard analysis. The analysis in Bayesian approaches is based on a posterior belief, also their special ability is to take into account the uncertainty of parameters in probabilistic relations and a priori knowledge. In this study, we benefited the Bayesian approach in order to estimate maximum values of peak ground acceleration (A_max) also quantiles of the relevant probabilistic distributions are figured out in a desired future interval time in Iran. The main assumptions are Poissonian character of the seismic events flow and properties of the Gutenberg-Richter distribution law. The map of maximum possible values of A_max and also map of 90% quantile of distribution of maximum values of A_max on a future interval time 100 years is presented. According to the results, the maximum value of the A_max is estimated for Bandar Abbas as 0.3g and the minimum one is attributed to Esfahan as 0.03g. Finally, the estimated values in Bayesian approach are compared with what was presented applying probabilistic seismic hazard (PSH) methods based on the conventional Cornel (1968) method. The distribution function of A_max for future time intervals of 100 and 475 years are calculated for confidence limit of probability level of 90%.     

Enhanced nonlinear static analysis with the drift pushover procedure for tall buildings

In this study a new method for nonlinear static analysis based on the relative displacements of stories is proposed that is able to be implemented in a single stage analysis and considers the effects of an arbitrary number of higher modes. The method is called the extended drift pushover analysis procedure (EDPA). To define the lateral load pattern, values of the relative displacements of stories are calculated using the elastic modal analysis and the modal combination factors introduced. For determining the combination factors, six different approaches are examined. Buildings evaluated in this study consist of four special steel moment-resisting frames with 10–30 stories. Responses including relative displacements of stories, story shear forces and rotation of plastic hinges in each story are calculated using the proposed approaches in addition to modal pushover analysis and nonlinear dynamic time history analyses. The nonlinear dynamic analysis is implemented using ten consistent earthquake records that have been scaled with regard to ASCE7-10. Distribution of response errors of story shears and plastic hinge rotations show that a major part of error corresponds to the second half of the buildings studied. Thus, the mentioned responses are corrected systematically. The final results of this study show that implementing the EDPA procedure using the third approach of this research is able to effectively overcome the limitations of both the traditional and the modal pushover analyses methods and predict the seismic demands of tall buildings with good accuracy.     

Critical frequencies <Emphasis Type="Italic">foF</Emphasis>2 as an indicator of trends in thermospheric dynamics

The time variations in three parameters during the last decades are considered. R(foF2) is the correlation coefficient between the nighttime and daytime values of foF2 for the same day. Stable trends are found for the minimum (R(foF2)(max)) and maximum (R(foF2)(min)) values of R(foF2) during a year. The foF2(night)/foF2(day) ratio demonstrates both, negative and positive trends, and the trend sign depends on the inclination I and declination D of the magnetic field. The correlation coefficient r(h, fo) between foF2 and the 100 hP level in the stratosphere demonstrates a decrease (in the years of maximum and minimum solar activity) from the 1980s to the 1990s. The trends in all three groups of data are considered under the assumption of long-term changes in the circulation in the upper atmosphere.     

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.     

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).     

Spatial and time variations in the <Emphasis Type="Italic">foF</Emphasis>2 (night)/<Emphasis Type="Italic">foF</Emphasis>2 (day) ratio: Specification of some effects

A series of important aspects of the time and spatial variations in the critical frequency ratio fo(night)/fo(day) considered by one of the authors [Danilov, 2007, 2008] is specified. The list of the Eastern-Hemisphere stations, for which an analysis of the above indicated ratio was performed, is completed. The available stations of the Western Hemisphere are considered. It is shown that the character of the variations in the fo(night)/fo(day) ratio is independent of limitations imposed on the Ap index of magnetic activity for the analyzed days. The consideration of the long-term variations in geomagnetic activity using the annual mean value of the Ap index does not influence the principal conclusions of the study, although in some cases changes statistical characteristics of the obtained behavior of the fo(night)/fo(day) ratio after 1980. A comparison of the results, based on the interpretation of the trends of the fo(night)/fo(day) value with the changes in the stratosphere-ionosphere coupling from the 1980s to the 1990s obtained by the authors earlier, confirms the assumption that there occurs a systematic change in the zonal wind in the upper atmosphere.     

Gutenberg-Richter <Emphasis Type="Italic">b</Emphasis>-value maximum likelihood estimation and sample size

The Aki-Utsu maximum likelihood method is widely used for estimation of the Gutenberg-Richter b-value, but not all authors are conscious of the method’s limitations and implicit requirements. The Aki/Utsu method requires a representative estimate of the population mean magnitude; a requirement seldom satisfied in b-value studies, particularly in those that use data from small geographic and/or time windows, such as b-mapping and b-vs-time studies. Monte Carlo simulation methods are used to determine how large a sample is necessary to achieve representativity, particularly for rounded magnitudes. The size of a representative sample weakly depends on the actual b-value. It is shown that, for commonly used precisions, small samples give meaningless estimations of b. Our results give estimates on the probabilities of getting correct estimates of b for a given desired precision for samples of different sizes. We submit that all published studies reporting b-value estimations should include information about the size of the samples used.     

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.