An ocean circulation model in <Emphasis Type="Italic">σ</Emphasis>S-<Emphasis Type="Italic">z</Emphasis>-<Emphasis Type="Italic">σ</Emphasis>B hybrid coordinate and its validation

A 3D, two-time-level, σS-z-σB hybrid-coordinate Marine Science and Numerical Modeling numerical ocean circulation model (HyMOM) is developed in this paper. In HyMOM, the σ coordinate is employed in the surface and bottom regions, and the z coordinate is used in the intermediate layers. This method can overcome problems with vanishing surface cells and minimize the unwanted deviation in representing bottom topography. The connection between the σ and z layers vertically includes an expanded “ghost” method and the linear interpolation. The governing equations in the σS-z-σB hybrid coordinate based on the complete Reynolds-averaged Navier-Stokes equations are derived in detail. The two-level time staggered and Eulerian forward and backward schemes, which are of second-order of accuracy, are adopted for the temporal difference in internal and external mode, respectively. The computation of the baroclinic gradient force is tested in an analytic test problem; the errors for two methods in HyMOM, which are relatively large only in the bottom layers, are obviously smaller than those in the pure σ and z models in almost all of the vertical layers. A quasi-global climatologic numerical experiment is constructed to test the simulation performance of HyMOM. With the monthly mean Levitus climatology data as reference, the HyMOM can improve the simulating accuracy compared with its pure z or σ coordinate implementation.     

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.     

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.     

Limitations of <Emphasis Type="Italic">H</Emphasis>-<Emphasis Type="Italic">κ</Emphasis> stacking: ambiguous results caused by crustal layering

Over the past decade, the H-κ stacking technique of Zhu and Kanamori (J Geophys Res 105:2969–2980, 2000) has become a standard tool to determine the crustal thickness H and the bulk crustal v_P/v_S ratio κ from teleseismic receiver functions. It is obvious that unfavorable noise conditions as well as a complex 3D velocity structure can severely hamper the interpretation of receiver-function data. However, we observe that ambiguities can even arise from a simple 1D layered velocity structure which raises a high potential for misinterpretations. To analyze the feasibility and basic limitations of the H-κ stacking method, we conduct a series of tests based on synthetic data. The impact of different given elementary parameters, related either to the velocity structure or to the data processing, is evaluated in a series of eight individual tests. We deliberately exclude complications such as 3D structural variations and/or noise to show that even a simple 1D velocity structure, involving, e.g., an additional inter-crustal discontinuity, can have significant consequences for the interpretation of the results. However, our modeling suggests that more complex crustal structures may lead to even less reliable results. Additionally, our tests illustrate that time shifts of the maxima in the H-κ domain due to the superposition and merging of individual phases can lead to significantly overestimated v_P/v_S ratios. In general, the depth to the Moho (or other discontinuities of interest) is less significantly affected. Our tests indicate the necessity to accurately check delay times derived from the maxima of the H-κ stacks against corresponding phases in the receiver functions. Repeating the stacking with varied weighting factors and filter ranges can help to reduce the ambiguities and to avoid possible misinterpretation.     

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.     

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.     

Improved random noise attenuation using <Emphasis Type="Italic">f</Emphasis>−<Emphasis Type="Italic">x</Emphasis> empirical mode decomposition and local similarity

Conventional f?x empirical mode decomposition (EMD) is an effective random noise attenuation method for use with seismic profiles mainly containing horizontal events. However, when a seismic event is not horizontal, the use of f?x EMD is harmful to most useful signals. Based on the framework of f?x EMD, this study proposes an improved denoising approach that retrieves lost useful signals by detecting effective signal points in a noise section using local similarity and then designing a weighting operator for retrieving signals. Compared with conventional f?x EMD, f?x predictive filtering, and f?x empirical mode decomposition predictive filtering, the new approach can preserve more useful signals and obtain a relatively cleaner denoised image. Synthetic and field data examples are shown as test performances of the proposed approach, thereby verifying the effectiveness of this method.     

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.     

Shape of Empirical and Synthetic Isoseismals: Comparison for Italian <Emphasis Type="Italic">M</Emphasis> ≤ 6 Earthquakes

We present results from a comparative analysis of empirical and synthetic shapes for isoseismals of low intensity (I = IV–VI on the MCS scale) for six Italian earthquakes of ML = 4.5–6. Our modeling of isoseismals is based on a plane-stratified earth model and on the double–couple point source approximation to calculate seismograms in the frequency range f 1 Hz. With a minimum of parameter variation we demonstrate that the low intensity isoseismals provide information on source geometry. We strive to avoid subjectivity in isoseismal constructions by using the new Diffuse Boundary method, which visualizes isoseismals with their uncertainty. Similar results in this direction are known for large earthquakes (M_L 6 or greater) with extended sources and for the higher isoseismals (I VI on the MM scale). The latter studies disregard the earth structure, use a greater number of parameters, and therefore have greater possibilities for fitting the data than our approach.     

Statistical analysis of ionospheric TEC anomalies before global <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> ≥ 7.0 earthquakes using data of CODE GIM

Based on Center of Orbit Determination in Europe (CODE) global ionospheric map (GIM) data, a statistical analysis of local total electron content (TEC) anomalies before 121 low-depth (D ≤ 100 km) strong (M _w ≥ 7.0) earthquakes has been made using the sliding median differential calculation method combining with a new approach of image processing technique. The results show that significant local TEC anomalies could be observed 0–6 days before 80 earthquakes, about 66.1% out of the total. The positive anomalies occur more often than negative ones. For 26 cases, both positive and negative anomalies are observed before the shock. The pre-earthquake TEC anomalies show local time recurrence for 38 earthquakes, which occur around the same local time on different days. The local time distribution of the pre-earthquake TEC anomalies mainly concentrates between 19 and 06 LT, roughly from the sunset to sunrise. Most of the pre-earthquake TEC anomalies do not locate above the epicenter but shift to the south. The pre-earthquake TEC anomalies could be extracted near the magnetic conjugate point of the epicenter for 40 events, which is 50% out of the total 80 cases with significant local TEC anomalies. In general, the signs of the anomalies around epicenter and its conjugate point are the same, but the abnormal magnitude and lasting time are not.     

Flow over a rounded backward-facing step,using a <Emphasis Type="Italic">z</Emphasis>-coordinate model and a <Emphasis Type="Italic">σ</Emphasis>-coordinate model

Homogeneous, nonrotating flow over a backward-facing rounded step is simulated using the 2D vertical version of two general circulation models, a z-coordinate model—the Massachusetts Institute of Technology general circulation model (MITgcm)—and a σ-coordinate model—the Bergen Ocean Model (BOM). The backward-facing step is a well-known testcase since it is geometrically simple but still embodies important flow characteristics such as separation point, reattachment length, and recirculation of the flow. The study compares the core of the two models and uses constant eddy viscosities and diffusivities. The Reynolds numbers ranges from 2·10² to 2·10⁶. The results correspond with previously published results having a relatively stationary separation point and a fluctuating reattachment length due to downslope propagating eddies released from the reattachment zone for Reynolds numbers higher than or equal to 2 · 10⁴. For Reynolds number within the laminar regime, the flow is stationary. The discrepancies between the models increase by enhancing Reynolds numbers. The σ-coordinate model experiences a reduction in eddy sizes with increasing resolution and Reynolds numbers in correspondence with published experiments, while the size of the eddies are independent of the Reynolds number using the MITgcm. Due to mixing generated by the staircase topography, the z-coordinate model gives a better convergence of the separation point and reattachment length compared with the BOM; however, this conclusion might change with the inclusion of a relevant turbulence scheme.     

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.     

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.     

The spatiotemporal analysis of the minimum magnitude of completeness <Emphasis Type="Italic">Mc</Emphasis> and the Gutenberg–Richter law <Emphasis Type="Italic">b</Emphasis>-value parameter using the earthquake catalog of Greece

Spatiotemporal mapping the minimum magnitude of completeness Mc and b-value of the Gutenberg–Richter law is conducted for the earthquake catalog data of Greece. The data were recorded by the seismic network of the Institute of Geodynamics of the National Observatory of Athens (GINOA) in 1970–2010 and by the Hellenic Unified Seismic Network (HUSN) in 2011–2014. It is shown that with the beginning of the measurements at HUSN, the number of the recorded events more than quintupled. The magnitude of completeness Mc of the earthquake catalog for 1970–2010 varies within 2.7 to 3.5, whereas starting from April 2011 it decreases to 1.5–1.8 in the central part of the region and fluctuates around the average of 2.0 in the study region overall. The magnitude of completeness Mc and b-value for the catalogs of the earthquakes recorded by the old (GINOA) and new (HUSN) seismic networks are compared. It is hypothesized that the magnitude of completeness Mc may affect the b-value estimates. The spatial distribution of the b-value determined from the HUSN catalog data generally agrees with the main geotectonic features of the studied territory. It is shown that the b-value is below 1 in the zones of compression and is larger than or equal to 1 in the zones dominated by extension. The established depth dependence of the b-value is pretty much consistent with the hypothesis of a brittle–ductile transition zone existing in the Earth’s crust. It is assumed that the source depth of a strong earthquake can probably be estimated from the depth distribution of the b-value, which can be used for seismic hazard assessment.     

Earthquake probabilities and magnitude distribution (<Emphasis Type="Italic">M</Emphasis>≥6.7) along the Haiyuan fault,northwestern China

In recent years, some researchers have studied the paleoearthquake along the Haiyuan fault and revealed a lot of paleoearthquake events. All available information allows more reliable analysis of earthquake recurrence interval and earthquake rupture patterns along the Haiyuan fault. Based on this paleoseismological information, the recurrence probability and magnitude distribution for M≥6.7 earthquakes in future 100 years along the Haiyuan fault can be obtained through weighted computation by using Poisson and Brownian passage time models and considering different rupture patterns. The result shows that the recurrence probability of M _S≥6.7 earthquakes is about 0.035 in future 100 years along the Haiyuan fault.     

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

Evidence for a Tang-Song Dynasty great earthquake along the Longmen Shan Thrust Belt prior to the 2008 <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> 7.9 Wenchuan earthquake,China

The magnitude (M _w) 7.9 Wenchuan earthquake occurred on 12 May 2008 in the Longmen Shan region of China, the transition zone between the Tibetan Plateau and the Sichuan Basin, resulting in widespread damage throughout central and western China. The steep, high-relief eastern margin of the Tibetan Plateau has undergone rapid Cenozoic uplift and denudation accompanied by folding and thrusting, yet no large thrust earthquakes are known prior to the 2008 M _w 7.9 Wenchuan earthquake. Field and excavation investigations reveal that a great historical earthquake occurred in the Sichuan region that ruptured a >200-km-long thrust fault within the Longmen Shan Thrust Belt, China, which also triggered the 2008 M _w 7.9 Wenchuan earthquake. The average co-seismic slip amount produced by this historical earthquake is estimated to be 2–3 m, comparable with that caused by the 2008 Wenchuan earthquake. Paleoseismic and archaeological evidence and radiocarbon dating results show that the penultimate great earthquake occurred in the Sichuan region during the late Tang-Song Dynasty, between AD 800 and 1000, suggesting a recurrence interval of ~1,000–1,200 years for Wenchuan-magnitude (M = ~8) earthquakes in the late Holocene within the Longmen Shan Thrust Belt. This finding is in contrast with previous estimates of 2,000–10,000 years for the recurrence interval of large earthquakes within the Longmen Shan Thrust Belt, as obtained from long-term slip rates based on the Global Positioning System and geological data, thereby necessitating substantial modifications to existing seismic-hazard models for the densely populated region at the eastern marginal zone of the Tibetan Plateau.     

<Emphasis Type="Italic">Prochlorococcus</Emphasis> viruses—From biodiversity to biogeochemical cycles

As the dominant primary producer in oligotrophic oceans, the unicellular picocyanobacterium Prochlorococcus is the smallest and most abundant photosynthetic phytoplankton in the world and plays an important role in marine carbon cycling. Cyanophages that infect Prochlorococcus influence the growth, carbon fixation, diversity, evolution, and environmental adaptation of their hosts. Here, we review studies on the isolation, genomics, and phylogenetic diversity of Prochlorococcus viruses and their interactions with Prochlorococcus. We also review the potential effects of Prochlorococcus viruses on biogeochemical cycling in the ocean.     

Deep structural characteristics and seismogenesis of the <Emphasis Type="Italic">M</Emphasis>≥8.0 earthquakes in North China

Based on the results from seismogeological study, aeromagnetic inversion and deep seismic sounding (DSS), it is found that the M≥8.0 earthquakes in North China have three common deep structural characteristics, i.e., they all took place above the ultra-crustal deep faults or on the edges of the tectonic blocks with higher intensity, and there are low-velocity, low-density and high-conductive layers deep in the epicentral regions. The origins of the earthquakes are also discussed and the two possibilities of seismogenesis are proposed, i.e., tectonic movement and intracrustal explosion.