Atypical soil behavior during the 2011 Tohoku earthquake (<Emphasis Type="Italic">М</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> = 9)

To understand physical mechanisms of generation of abnormally high peak ground acceleration (PGA; >1g) during the Tohoku earthquake, models of nonlinear soil behavior in the strong motion were constructed for 27 KiK-net stations located in the near-fault zones to the south of FKSH17. The method of data processing used was developed by Pavlenko and Irikura, Pure Appl Geophys 160:2365–2379, 2003 and previously applied for studying soil behavior at vertical array sites during the 1995 Kobe (М _w ?=?6.8) and 2000 Tottori (М _w ?=?6.7) earthquakes. During the Tohoku earthquake, we did not observe a widespread nonlinearity of soft soils and reduction at the beginning of strong motion and recovery at the end of strong motion of shear moduli in soil layers, as usually observed during strong earthquakes. Manifestations of soil nonlinearity and reduction of shear moduli during strong motion were observed at sites located close to the source, in coastal areas. At remote sites, where abnormally high PGAs were recorded, shear moduli in soil layers increased and reached their maxima at the moments of the highest intensity of the strong motion, indicating soil hardening. Then, shear moduli reduced with decreasing the intensity of the strong motion. At soft-soil sites, the reduction of shear moduli was accompanied by a step-like decrease of the predominant frequencies of motion. Evidently, the observed soil hardening at the moments of the highest intensity of the strong motion contributed to the occurrence of abnormally high PGA, recorded during the Tohoku earthquake.     

Simulating CO<Subscript>2</Subscript> transport into the ocean from a CO<Subscript>2</Subscript> lake at the seafloor using a <Emphasis Type="Italic">z</Emphasis>- and a <Emphasis Type="Italic">σ</Emphasis>-coordinate model

The ocean takes up approximately 2 GT carbon per year due to the enhanced CO₂ concentrations in the atmosphere. Several options have been suggested in order to reduce the emissions of CO₂ into the atmosphere, and among these are CO₂ storage in the deep ocean. Topographic effects of dissolution and transport from a CO₂ lake located at 3,000-m depth have been studied using the z-coordinate model Massachusetts Institute of Technology general circulation model (MITgcm) and the σ-coordinate model Bergen ocean model (BOM). Both models have been coupled with the general ocean turbulence model (GOTM) in order to account for vertical subgrid processes. The chosen vertical turbulence mixing scheme includes the damping effect from stable stratification on the turbulence intensity. Three different topographic scenarios are presented: a flat bottom and the CO₂ lake placed within a trench with depths of 10 and 20 m. The flat case scenario gives good correlation with previous numerical studies of dissolution from a CO₂ lake. When topography is introduced, it is shown that the z-coordinate model and the σ-coordinate model give different circulation patterns in the trench. This leads to different dissolution rates, 0.1 μmol cm^− 2 s^− 1 for the scenario of a 20-m-deep trench using BOM and 0.005–0.02 μmol cm^− 2 s^− 1 for the same scenario using the MITgcm. The study is also relevant for leakages of CO₂ stored in geological formations and to the ocean.     

Strong motion characteristics of the <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> 6.6 Lushan earthquake,Sichuan, China — an insight into the spatial difference of a typical thrust fault earthquake

Near-field strong ground motions are useful for engineering seismology studies and seismic design, but dense observation networks of damaging earthquakes are still rare. In this study, based on the strong-motion data from the M w 6.6 Lushan earthquake, the ground motion parameters in different spatial regions are systematically analyzed, and the contributions from different effects, like the hanging-wall effect, directivity effect, and attenuation effect are separated to the extent possible. Different engineering parameters from the observed ground motions are compared with the local design response spectra and a new attenuation relation of Western China. General results indicate that the high frequency ground motion, like the peak ground acceleration, on two sides of the fault plane is sensitive to the hanging-wall effect, whereas the low frequency ground motion, like the long period spectral acceleration, in the rupture propagation direction is affected by the directivity effect. Moreover, although the M w 6.6 Lushan earthquake is not a large magnitude event, the spatial difference of ground motion is still obvious; thus, for a thrust faulting earthquake, in addition to the hanging effect, the directivity effect should also be considered.     

IR band of O<Subscript>2</Subscript> at 1.27 μm as the tracer of O<Subscript>3</Subscript> in the mesosphere and lower thermosphere: Correction of the method

The problem of systematic overestimation (20–50%) of the retrieved ozone concentrations in the altitude range of 60–80 km in the TIMED–SABER satellite experiment in the daytime has been solved. The reason for overestimation is the neglect of the electronic vibrational kinetics of photolysis products of ozone and molecular oxygen O₂(b¹Σ_g ⁺, ν) and O₂(a¹Δ_g, ν). The IR emission band of O₂(a¹Δ_g, ν = 0) at 1.27 μm can be correctly used in remote sensing in order to obtain the ozone altitude profile in the altitude range of 50–88 km only with the use of a complete model of electronic vibrational kinetics of O₂ and O₃ photolysis products (YM2011) in the Earth’s mesosphere and lower thermosphere. Alternative ozone tracers have been considered, and an optimum tracer in the altitude range of 50–100 km such as O₂(b¹Σ_g ⁺, ν = 1) molecule emissions has been proposed.     

New characteristics of intensity assessment of Sichuan Lushan “4.20” <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>7.0 earthquake

The post-earthquake rapid accurate assessment of macro influence of seismic ground motion is of significance for earthquake emergency relief, post-earthquake reconstruction and scientific research. The seismic intensity distribution map released by the Lushan earthquake field team of the China Earthquake Administration (CEA) five days after the strong earthquake (M7.0) occurred in Lushan County of Sichuan Ya’an City at 8:02 on April 20, 2013 provides a scientific basis for emergency relief, economic loss assessment and post-earthquake reconstruction. In this paper, the means for blind estimation of macroscopic intensity, field estimation of macro intensity, and review of intensity, as well as corresponding problems are discussed in detail, and the intensity distribution characteristics of the Lushan “4.20” M7.0 earthquake and its influential factors are analyzed, providing a reference for future seismic intensity assessments.     

Ionospheric <Emphasis Type="Italic">N</Emphasis><Subscript>e</Subscript> disturbances before 2007 Pu’er,Yunnan, China,earthquake

Selecting three half orbits near the epicenter of Pu’er earthquake, we analyzed the N_e data recorded in their revisited orbits during a year before this earthquake, and extracted N_e precursors. The results show that: ① There are significant seasonal variations of ionospheric N_e in night time, which exhibit different shapes respectively in four seasons; ② There are three main shapes of N_e: single-peak, saddle-shaped and even-shaped, all of which may occur in four seasons, but each season with its typical shape relatively; ③ Spatial images of N_e showed high values near the epicenter in 30 days before the earthquake, and there is a good correlation between anomaly and distribution of earthquake in space and time, which reflects that these spatial anomalies were indeed concerned with the earthquake; ④ There shows a certain similarity of the N_e curves among revisited orbits, which can provide background information for distinguishing and identification of seismic anomaly.     

Hovsgol earthquake 5 December 2014, <Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript> = 4.9: seismic and acoustic effects

A moderate shallow earthquake occurred on 5 December 2014 (M _W = 4.9) in the north of Lake Hovsgol (northern Mongolia). The infrasonic signal with duration 140 s was recorded for this earthquake by the “Tory” infrasound array (Institute of Solar-Terrestrial Physics of the Siberian Branch of the Russian Academy of Science, Russia). Source parameters of the earthquake (seismic moment, geometrical sizes, displacement amplitudes in the focus) were determined using spectral analysis of direct body P and S waves. The spectral analysis of seismograms and amplitude variations of the surface waves allows to determine the effect of the propagation of the rupture in the earthquake focus, the azimuth of the rupture propagation direction and the velocity of displacement in the earthquake focus. The results of modelling of the surface displacements caused by the Hovsgol earthquake and high effective velocity of propagation of infrasound signal (~ 625 m/s) indicate that its occurrence is not caused by the downward movement of the Earth’s surface in the epicentral region but by the effect of the secondary source. The position of the secondary source of infrasound signal is defined on the northern slopes of the Khamar-Daban ridge according to the data on the azimuth and time of arrival of acoustic wave at the Tory station. The interaction of surface waves with the regional topography is proposed as the most probable mechanism of formation of the infrasound signal.     

Estimation of crustal Q<Subscript>β</Subscript> in the NW Himalaya using teleseismic broadband SH waveforms of the 8 October 2005 South Asian earthquake

The average value of the intrinsic shear wave quality factor, Q _β , for the 15 km of the upper crust in the NW Himalaya is estimated. Thirty-two teleseismic broadband SH and sSH waveforms of 8 October 2005 South Asian earthquake (Mw = 7.6), from sixteen Global Seismographic Network stations of the National Earthquake Information Center network have been used. The selected windows of the direct and depth phases have been Fourier transformed and smoothed using the Hamming and Tuckey technique. Ratio of the smoothed spectra of depth and direct phases are obtained in the frequency range 0.2–1.5 Hz. A straight line is fitted in the least-square sense to the spectral ratio versus frequency. The value of Q _β is estimated from the slope of the line. The frequency independent average value of Q _β is estimated to be 218±56.     

Multi-dipole observation system and study on the abnormal variation of the geoelectric field observed at Capital Circle area before the Wen’an <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>5.1 earthquake

On July 4, 2006, a M_S5.1 earthquake occurred in Wen’an county of Hebei Province of which the epicenter is near the Beijing city. The six geoelectric field monitoring stations have been in operation for several years around the Beijing area to examine the relationship between electric field changes and earthquake. This paper firstly explains the principle of the eliminating noise method by using multi-dipole observation system of the geoelectric field. Then the data observed at the stations are studied and a lot of abnormal signals preceding the Wen’an earthquake are selected, of which five abnormal signals of the geoelectric field are finally recognized as the precursory signals. The result shows that ? there probably exists the precursory signals of the geoelectric field preceding the Wen’an earthquake; ? there are sensitive sites in the spatial distribution of the abnormal variation of the geoelectric field before the quack; ? the anomalous signals do not appear synchronously, and their durations are not same at different stations; ? the amplitudes of the abnormal signals recorded at Baodi station are small, but large at Changli station, while the latter is farther from the epicentral area than the former.     

Source mechanism of strong aftershocks (<Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>⩾5.6) of the 2008/05/12 Wenchuan earthquake and the implication for seismotectonics

Dozens of >M5, hundreds of >M4, and much more >M3 aftershocks occurred after the 2008/05/12 Wenchuan earthquake, which were well recorded by permanent and portable seismic stations. After relocated with P arrival, the >M3 aftershocks show two trends of distribution, with most of the aftershocks located along the north-east strike consistent with Longmenshan fault system, yet there is a north-west trend around the epicenter. It seems that substantially more aftershocks occur in regions with crystalline bedrocks. Then we collected waveform data from National Digital Seismograph Network and regional seismograph network of China, and employed “Cut and Paste” method to obtain focal mechanisms and depths of the big aftershocks (M⩾5.6). While most of those aftershocks show thrust mechanism, there are some strike slip earthquakes in the northern-most end of the rupture. Focal mechanisms show that the events located on the southern part of central Beichuan-Yingxiu Fault (BY) are mainly thrust earthquakes, which is consistent with initial mechanism of the main shock rupture. In the north part the aftershocks along the BY are also dominated by thrust slip, which is quite different from the right slip rupture of the main shock. Around Qingchuan-Pingwu Fault, the focal mechanisms are dominated by right-slip rupture with large depths (∼18 km). So we suspected that in the north part the main shock might rupture on two faults: Beichuan Fault and Qingchuan-Pingwu Fault. The complex pattern of aftershock mechanisms argues for presence of a complicated fault system in the Longmenshan area. Supported by Knowledge Innovation Project of Chinese Academy of Sciences (Grant Nos. KZCX3-SW-153, KZCX2-YW-116-1), National Natural Science Foundation of China (Grant No. 40604004), and National Basic Technology R & D Program (Grant No. 2006BAC01B02-01-02).     

Structural analysis and interpretation of the surface deformation associated with the Ning’er,Yunnan Province,China <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>6.4 earthquake of June 3, 2007

The seismogenic fault and the dynamic mechanism of the Ning’er, Yunnan Province M_S6.4 earthquake of June 3, 2007 are studied on the basis of the observation data of the surface fissures, sand blow and water eruption, landslide and collapse associated with the earthquake, incorporating with the data of geologic structures, focal mechanism solutions and aftershock distribution for the earthquake area. The observation of the surface fissures reveals that the Banhai segment of the NW-trending Ning’er fault is dominated by right-lateral strike-slip, while the NNE-trending fault is dominated by left-lateral strike-slip. The seismo-geologic hazards are concentrated mainly within a 330°-extending zone of 13.5 km in length and 4 km in width. The major axis of the isoseismal is also oriented in 330° direction, and the major axis of the seismic intensity VIII area is 13.5 km long. The focal mechanism solutions indicate that the NW-trending nodal plane of the Ning’er M_S6.4 earthquake is dominated by right-lateral slip, while the NE-trending nodal plane is dominated by left-lateral slip. The preferred distribution orientation of the aftershocks of M_S≥2 is 330°, and the focal depths are within the range of 3～12 km, predominantly within 3～10 km. The distribution of the aftershocks is consistent with the distribution zone of the seismo-geologic hazards. All the above-mentioned data indicate that the Banhai segment of the Ning’er fault is the seismogenic fault of this earthquake. Moreover, the driving force of the Ning’er earthquake is discussed in the light of the active block theory. It is believed that the northward pushing of the Indian plate has caused the eastward slipping of the Qinghai-Tibetan Plateau, which has been transformed into the southeastern-southernward squeezing of the southwest Yunnan region. As a result, the NW-trending faults in the vicinity of the Ning’er area are dominated by right-lateral strike-slip, while the NE-trending faults are dominated by left-lateral strike-slip. This tectonic framework might be the main cause of the frequent occurrence of M_S6.0～6.9 earthquakes in the area.     

Temporal and spatial variations of <Emphasis Type="Italic">δ</Emphasis><Superscript>15</Superscript>N and <Emphasis Type="Italic">δ</Emphasis><Superscript>18</Superscript>O for atmospheric N<Subscript>2</Subscript>O above the oceanic surface from Shanghai to Antarctica

During the 22nd Chinese Antarctic Research Expedition (CHINARE-22), the atmospheric gas samples above the oceanic surface and near the surface were collected on the track for the scientific ship “Xuelong” and on Millor Peninsula of eastern Antarctica, respectively, using the Tedlar gas bags. Every day the sampling times were 10:00 and 22:00 (local time), respectively. In the laboratory, high-precision measurement of the isotopic compositions for N₂O in these gas samples was conducted using Thermo Finnigan MAT-253 Isotopic Mass Spectrometer with a fully automated interface for the pre-GC concentration (PreCon) of trace gases. The temporal and spatial variations of δ ¹⁵N and δ ¹⁸O in atmospheric N₂O were analyzed. The mean δ ¹⁵N and δ ¹⁸O-N₂O values above the oceanic surface were (7.21±0.50)‰ and (44.52±0.52)‰, respectively. From 30°N to Antarctica, the δ ¹⁵N (6.05‰–7.88‰) linearly increased with the rate of about 0.01‰ with the latitude while the δ ¹⁸O (43.05‰–48.78‰) showed a large fluctuation. The δ ¹⁵N negatively correlated with air temperature and N₂O concentration, and slightly positively correlated with δ ¹⁸O. The summertime variations of δ ¹⁵N and δ ¹⁸O-N₂O appeared the same trend on Millor Peninsula of eastern Antarctica. They significantly positively correlated with each other and negatively with N₂O concentration. The δ ¹⁵N and δ ¹⁸O-N₂O at different sites averaged (7.42±0.35)‰ and (44.69±0.49)‰, respectively, slightly higher than those above the oceanic surface, significantly higher than those of atmospheric N₂O in the low-latitude regions of Northern Hemisphere. The predominant factors affecting the spatial variations of δ ¹⁵N and δ ¹⁸O values were also discussed. The isotopic data given in this study can help to investigate the global and regional N₂O budgets. Supported by the National Natural Science Foundation of China (Grant Nos. 40676005 and 40406001)     

The Tolud Burst of Seismicity and the Earthquake of November 30, 2012 (<Emphasis Type="Italic">M</Emphasis><Subscript>C</Subscript> = 5.4, <Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript> = 4.8) that Accompanied the Start of the 2012‒2013 Tolbachik Eruption

This paper reports a study of the Tolud earthquake sequence; the sequence was a burst of shallow seismicity between November 28 and December 7, 2012; it accompanied the initial phase in the Tolbachik Fissure Eruption of 2012?2013. The largest earthquake (the Tolud earthquake of November 30, 2012, to be referred to as the Tolud Earthquake in what follows, with K_S = 11.3, M_L = 4.9, M_C = 5.4, and M_W = 4.8) is one of the five larger seismic events that have been recorded at depths shallower than 10 km beneath the entire Klyuchevskoi Volcanic Cluster in 1961?2015. It was found that the Tolud earthquake sequence was the foreshock–aftershock process of the Tolud Earthquake. This is one of the larger seismicity episodes ever to have occurred in the volcanic areas of Kamchatka. Data of the Kamchatka seismic stations were used to compute some parameters for the Tolud Earthquake and its largest (M_L = 4.3) aftershock; the parameters include the source parameters and mechanisms, and the moment magnitudes, since no information on these is available at the world seismological data centers. The focal mechanisms for the Tolud Earthquake and for its aftershock are consistent with seismic ruptures at a tension fault in the rift zone. Instrumental data were used to estimate the intensity of shaking due to the Tolud Earthquake. We discuss the sequence of events that was a signature of the time-dependent seismic and volcanic activity that took place in the Tolbachik zone in late November 2012 and terminated in the Tolud burst of seismicity. Based on the current ideas of the tectonics and magma sources for the Tolbachik volcanic zone, we discuss possible causes of these earthquakes.     

Was the M<Subscript>w</Subscript> 7.5 1952 Kern County,California, earthquake induced (or triggered)?

Several recent studies have presented evidence that significant induced earthquakes occurred in a number of oil-producing regions during the early and mid-twentieth century related to either production or wastewater injection. We consider whether the 21 July 1952 M_w 7.5 Kern County earthquake might have been induced by production in the Wheeler Ridge oil field. The mainshock, which was not preceded by any significant foreshocks, occurred 98 days after the initial production of oil in Eocene strata at depths reaching 3 km, within ~1 km of the White Wolf fault (WWF). Based on this spatial and temporal proximity, we explore a potential causal relationship between the earthquake and oil production. While production would have normally be expected to have reduced pore pressure, inhibiting failure on the WWF, we present an analytical model based on industry stratigraphic data and best estimates of parameters whereby an impermeable splay fault adjacent to the main WWF could plausibly have blocked direct pore pressure effects, allowing the poroelastic stress change associated with production to destabilize the WWF, promoting initial failure. This proof-of-concept model can also account for the 98-day delay between the onset of production and the earthquake. While the earthquake clearly released stored tectonic stress, any initial perturbation on or near a major fault system can trigger a larger rupture. Our proposed mechanism provides an explanation for why significant earthquakes are not commonly induced by production in proximity to major faults.     

The seasonal behavior of the <Emphasis Type="Italic">B</Emphasis><Subscript><Emphasis Type="Italic">z</Emphasis></Subscript> daily variations according to the monitoring at the Bishkek geodynamic site (1993–1995)

The algorithm for monitoring the magnetotelluric field has been briefly described, and the experimental data, obtained at the Bishkek geodynamic site, have been analyzed. The regularities in the seasonal behavior of the variations in the field vertical magnetic component have been revealed.     

QuakeIST<Superscript>®</Superscript> earthquake scenario simulator using interdependencies

Earthquakes are a permanent threat to urban environments worldwide. The communication of the related risk demands accurate damage model simulations and an interactive visualization of results. The aim of this paper is to provide a realistic problem-solving environment for earthquake discussions among decision makers, stakeholders, and the general public. QuakeIST^® is an integrated earthquake simulator developed by Instituto Superior Técnico (Lisbon University), oriented towards the performance of risk calculations concerning damage propagations that use the Disruption Index concept. This software imports data stored in a GIS environment, handles different ground motion scenarios, and deals with a complex situation of different soils and vulnerabilities of various layers of civil structures (buildings, lifelines, and other urban structures). It models interdependencies between several infrastructures and between infrastructures and the urban tissue. The computer programme is very versatile, written in separate modules, allowing an experimented user to incorporate new formulations. Results can be treated with any statistical application and most common GIS commercial environments can produced their geographic visualization. Current progress and new upcoming are briefly described at the end of the paper.     

The <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> = 5.8 14 August 2016 middle Sakhalin earthquake on a boundary between Okhotsk and Eurasian (Amurian) plates

An earthquake with the moment magnitude M _w ?=?5.8 occurred in the middle part of the Sakhalin Island, Russian Federation, on 14 August 2016, at 11:17 a.m. UTC. The earthquake source was located west of the Central Sakhalin Fault Zone, which is considered to mark the boundary between the Okhotsk and Eurasian (Amurian) plates. Moment tensor solution of the mainshock as well as the configuration of aftershock cloud suggests that the earthquake was caused by slip on a SW-dipping reverse fault. For the first time for Sakhalin, we have got the felt reports unified in accordance with DYFI. We also analyzed observed PGA values and, based on them, produced shaking maps.     

Evaluation of the dynamic <Emphasis Type="Italic">p</Emphasis>–<Emphasis Type="Italic">y</Emphasis><Subscript>p</Subscript> loops of pile-supported structures on sloping ground

A key issue in the design of pile-supported structures on sloping ground is soil–pile interaction, which becomes more complicated in case of dynamic loading. This study aimed to evaluate the effect of slope on the dynamic behavior of pile-supported structures by performing a series of centrifuge tests. Three models were prepared by varying the slope and soil density of dry sand grounds. The mass supported on 3 by 3 group piles was shaken applying sinusoidal wave with various amplitudes. Test results showed that the location of maximum values and distribution shape of the bending moment below the ground surface varied noticeably with the pile position in the slope case. The relationship between the soil resistance and pile deflection (p–y_p loops) was carefully evaluated by applying the piecewise cubic spline method to fit the measured bending moment curves along piles. It was found that the shape of the p–y_p loops was irregular due to the effect of slope, and immensely influenced by the movement of the unstable zone. In addition, the effect of the pile group in the horizontal case was evaluated by comparing with the previously suggested curves that represent the relationship between the soil resistance and pile–soil relative displacement (p–y curves) to propose the multiplier coefficients.     

Methods and technique of measurements of ρ<Subscript><Emphasis Type="Italic">k</Emphasis></Subscript> variations by GEOMES-R6 device

Upgrading of the electrical resistance variometer GEOMES-R6 of the Institute of Geophysics of the Poland Academy of Sciences, with enhancement of functional capabilities of the GEOMES assemblage for its application in laboratory measurements on rock samples, has been carried out. The assemblage operation principle, specifications, and theory of operation are given. Results of the assemblage test at measurements on marble sample and concrete model are presented. It is shown that the installation allows long-duration registration of variations of the value ρ _k on high-ohmic samples.