Thaw-induced slope failures and stability analyses in permafrost regions of the Qinghai-Tibet Plateau,China

The distribution of permafrost-related slope failures along the Qinghai-Tibet Highway from Wuddaoliang to Fenghuoshan correlates with ice content, slope gradient, and ground temperature. Slope failures are of two types. (1) Retrogressive thaw slumps result from icy permafrost being exposed by either man-induced excavation or fluvial-thermal erosion and undercutting of basal slopes. (2) Active-layer-detachment failures are caused by thaw of icy permafrost at the active layer-permafrost interface. After initial failure, active-layer-detachment failures can lead to retrogressive thaw-slumping and localized surficial landslide. Common trigger mechanisms for failure include high summer air temperatures and heavy summer precipitation. A third possible trigger mechanism for slope failure is earthquake occurrence. A geotechnical slope stability analysis was undertaken for an active-layer-detachment failure that had progressed into a retrogressive thaw slump. A safety factor (F _s ) of 1.24 for the natural slope was determined using in situ tested strength parameters. However, the slope would lose stability when either the groundwater level over the permafrost table exceeded 1.42 m or seismic acceleration reached, or exceeded, 0.03?g.     

Geometrical characteristics of earthquake-induced landslides and correlations with control factors: a case study of the 2013 Minxian,Gansu, China, <Emphasis Type="Italic">Mw</Emphasis> 5.9 event

Geometric parameters are useful for characterizing earthquake-triggered landslides. This paper presents a detailed statistical analysis on this issue using the landslide inventory of the 2013, Minxian, China Mw 5.9 earthquake. Based on GIS software and a 5-m resolution DEM, geometric parameters of 635 coseismic landslides (with areas larger than 500 m²) were obtained, including height, length, width, reach angle (arc tangent of the height-length ratio), and aspect ratio (length-width ratio). The fitting relationship of height and length from these data is H = 0.6164L + 0.4589, with an average reach angle of 31.65°. The landslide aspect ratios concentrate in the range of 1.4～2.6, with an average of 2.11. According to the plane geometric shapes and aspect ratios, the landslides are classified into four categories: transverse landslide (LA1, L/W ≤ 0.8), isometric landslide (LA2, 0.8 < L/W ≤ 1.2), longitudinal landslide (LA3, 1.2 < L/W ≤ 3), and elongated landslide (LA4, L/W > 3). Statistics of these four types of landslides versus ten classified control factors (elevation, slope angle, slope aspect, curvature, slope position, distance to drainages, lithology, seismic intensity, peak ground acceleration, and distance to seismogenic fault) are used to examine their possible correlations and the landslide-prone areas, which would be helpful to the landslide disaster mitigation in the affected area.     

Shear wave velocity imaging of landslide debris deposited on an erodible bed and possible movement mechanism for a loess landslide in Jingyang,Xi’an,China

The South Jingyang Plateau, with a total area of 70 km², is located in Shaanxi Province, China. Since 1976, more than 50 landslides of different types have occurred repeatedly on the edge slopes of the plateau due to the start of diversion irrigation on the plateau, resulting in great loss of lives and property. To better understand the initiation and movement mechanisms of these loess landslides, we surveyed them and carried out a detailed investigation of a large landslide in the Xihetan area. Our field survey results revealed that although most of these landslides had a long runout with high mobility, most of the landslide materials originating from the edge slopes may have been in an unsaturated state when the landslide occurred. This suggests that the materials at the toe of the edge slope as well as on the travel path along the river terrace might have played a key role in landslide movement. To examine how the materials on the travel path were involved in the landsliding, we used a multichannel surface wave technique and surveyed shear wave velocity (V _s ) profiles of the landslide deposits. We also examined the internal geometry of the deposits that outcropped on the right-side slope of the landslide foot. The longitudinal profile of V _s along the direction of movement showed that terrace deposits near the toe of the edge slope may have been sheared upward, indicating that at the toe, the surface of rupture might be located inside the terrace deposits. The V _s contours showed an A-shaped fold within the landslide deposits in the middle part of the travel path and became greater in the most distal toe part. The V _s profile across the deposits showed a U-shaped belt, in which the soil layers have smaller V _s . This belt may be the boundary between the sliding landslide debris and terrace deposits. The observed internal geometry of the landslide deposits indicates that a sliding surface developed within the sandy layer underlying the gravel layer. Therefore, we inferred that after failure, the displaced landslide materials overrode and sheared the terrace deposits along its main sliding direction, resulting in the formation of thrust folds within the terrace deposits, and greater V _s on the distal toe part of the landslide.     

Field monitoring of groundwater responses to heavy rainfalls and the early warning of the Kualiangzi landslide in Sichuan Basin,southwestern China

The Kualiangzi landslide was triggered by heavy rainfalls in the “red beds” area of Sichuan Basin in southwestern China. Differing from other bedrock landslides, the movement of the Kualiangzi landslide was controlled by the subvertical cracks and a subhorizontal bedding plane (dip angle < 10°). The ingress of rainwater in the cracks formed a unique groundwater environment in the slope. Field measurement for rainfall, groundwater movement, and slope displacement has been made for the Kualiangzi landslide since 2013. The field monitoring system consists of two rainfall gauges, seven piezometers, five water-level gauges, and two GPS data loggers. The equipments are embedded near a longitudinal section of the landslide, where severe deformation has been observed in the past 3 years. The groundwater responses to four heavy rainfall events were analyzed between June 16 and July 24 in 2013 coincided with the flood season in Sichuan. Results showed that both of the water level and the pore-water pressure increased after each rainfall event with delay in the response time with respect to the precipitation. The maximum time lag reached 35 h occurred in a heavy rainfall event with cumulative precipitation of 127 mm; such lag effect was significantly weakened in the subsequent heavy rainfall events. In each presented rainfall event, longer infiltration period in the bedrock in the upper slope increased the response time of groundwater, compared to that of in the gravels in the lower slope. A translational landslide conceptual model was built for the Kualiangzi landslide, and the time lag was attributed to the gradual formation of the uplift pressure on the slip surface and the softening of soils at the slip surface. Another important observation is the effect on the slope movement which was caused by the water level (H _w) in the transverse tension trough developed at the rear edge of the landslide. Significant negative correlation was found for H _w and the slope stability factor (F _s), in particular for the last two heavy rainfall events, of which the drastic increase of water level caused significant deterioration in the slope stability. The rapid drop (Δ?=?22.5 kPa) of pore-water pressure in the deep bedrock within 1 h and the large increase (Δ?=?87.3 mm) of surficial displacement were both monitored in the same period. In the end, a four-level early warning system is established through utilizing H _w and the displacement rate D _r as the warning indicators. When the large deformation occurred in flood season, the habitants at the leading edge of the landslide can be evacuated in time.     

Coseismic landslides triggered by the 8th August 2017 <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript> 7.0 Jiuzhaigou earthquake (Sichuan,China): factors controlling their spatial distribution and implications for the seismogenic blind fault identification

On 8th August 2017, a magnitude M_s 7.0 earthquake struck the County of Jiuzhaigou, in Sichuan Province, China. It was the third M_s ≥?7.0 earthquake in the Longmenshan area in the last decade, after the 2008 M_s 8.0 Wenchuan earthquake and the 2013 M_s 7.0 Lushan earthquake. The event did not produce any evident surface rupture but triggered significant mass wasting. Based on a large set of pre- and post-earthquake high-resolution satellite images (SPOT-5, Gaofen-1 and Gaofen-2) as well as on 0.2-m-resolution UAV photographs, a polygon-based interpretation of the coseismic landslides was carried out. In total, 1883 landslides were identified, covering an area of 8.11 km², with an estimated total volume in the order of 25–30?×?10⁶ m³. The total landslide area was lower than that produced by other earthquakes of similar magnitude with strike-slip motion, possibly because of the limited surface rupture. The spatial distribution of the landslides was correlated statistically to a number of seismic, terrain and geological factors, to evaluate the landslide susceptibility at regional scale and to identify the most typical characteristics of the coseismic failures. The landslides, mainly small-scale rockfalls and rock/debris slides, occurred mostly along two NE-SW-oriented valleys near the epicentre. Comparatively, high landslide density was found at locations where the landform evolves from upper, broad valleys to lower, deep-cut gorges. The spatial distribution of the coseismic landslides did not seem correlated to the location of any known active faults. On the contrary, it revealed that a previously-unknown blind fault segment—which is possibly the north-western extension of the Huya fault—is the plausible seismogenic fault. This finding is consistent with what hypothesised on the basis of field observations and ground displacements.     

Dynamics of carbon fluxes with responses to vegetation,meteorological and terrain factors in the south-eastern Tibetan Plateau

Tibetan Plateau (TP) is the highest and most extensive plateau in the world and has been known as the roof of the world, and it is sensitive to climate change. The researches of CO₂ fluxes (F _C) in the TP region play a significant role in understanding regional and global carbon balance and climate change. Eddy covariance flux measurements were conducted at three sites of south-eastern TP comprising Dali (DL, cropland ecosystem), LinZhi (LZ, alpine meadow ecosystem) and Wenjiang (WJ, cropland ecosystem); amongst those DL and LZ are located in plateau region, while WJ is in plain region. Dynamics of F _C and influences of vegetation, meteorological (air temperature, photosynthetically active radiation, soil temperature and soil water content) and terrain factors (altitude) were analysed on the basis of data taken during 2008. The results showed that, in the cool sub-season (March, April, October and December), carbon sink appeared even in December with fluxes of (?0.021 to ?0.05) mg CO₂ m^?2 s^?1 and carbon source only in October (0.03 ± 0.0048) mg CO₂ m^?2 s^?1 in DL and WJ site. In LZ site, carbon sink was observed in April: (?0.036 ± 0.0023) mg CO₂ m^?2 s^?1 and carbon sources in December and March (0.008–0.010 mg CO₂ m^?2 s^?1). In the hot sub-season (May–August), carbon source was observed only in May with (0.011 ± 0.0022), (0.104 ± 0.0029) and (0.036 ± 0.0017) fluxes in LZ, DL and WJ site, respectively, while carbon sinks with (?0.021 ± 0.0041), (?0.213 ± 0.0007) and (?0.110 ± 0.0015) mg CO₂ m^?2 s^?1 fluxes in LZ, DL, and WJ, respectively. Comparing with plain region (WJ), carbon sinks in plateau region (DL and LZ) lasted for a longer time, and the absorption sum was large and up to (–357.718 ± 0.0054) and (?371.111 ± 0.0039) g C m^?2 year^?1, respectively. The LZ site had the weakest carbon sink with (?178.547 ± 0.0070) g C m^?2 year^?1. Multivariate analysis of covariance showed that altitude (AL) as an independent factor explained 39.5 % of F _C (P < 0.026). F _C had a quadratic relationship with Normalized difference vegetation index (NDVI) (R ² ranges from 0.485 to 0.640 for three sites), an exponential relationship with soil temperature at 5-cm depth (ST ₅) at night time and a quadratic relationship with air temperature (T _a) at day time. Path analysis indicated that photosynthetically active radiation (PAR), sensible heat fluxes (H) and other factors all had direct or indirect effects on F _C in all of the three tested sites around the south-eastern TP.     

Precursory seismicity rate changes prior to the large and major earthquakes along the Sagaing fault zone,Central Myanmar

In this study, the seismicity rate changes that can represent an earthquake precursor were investigated along the Sagaing Fault Zone (SFZ), Central Myanmar, using the Z value technique. After statistical improvement of the existing seismicity data (the instrumental earthquake records) by removal of the foreshocks and aftershocks and man-made seismicity changes and standardization of the reported magnitude scales, 3574 earthquake events with a M _w ≥ 4.2 reported during 1977–2015 were found to directly represent the seismotectonic activities of the SFZ. To find the characteristic parameters specifically suitable for the SFZ, seven known events of M _w ≥ 6.0 earthquakes were recognized and used for retrospective tests. As a result, utilizing the conditions of 25 fixed earthquake events considered (N) and a 2-year time window (T _w), a significantly high Z value was found to precede most of the M _w ≥ 6.0 earthquakes. Therefore, to evaluate the prospective areas of upcoming earthquakes, these conditions (N = 25 and T _w = 2) were applied with the most up-to-date seismicity data of 2010–2015. The results illustrate that the vicinity of Myitkyina and Naypyidaw (Z = 4.2–5.1) cities might be subject to strong or major earthquakes in the future.     

Impact of pore water conductivity and porosity on the electrical conductivity of kaolinite

The purpose of this study is to quantify the magnitudes of surface conduction and pore water conduction from the measured electrical conductivity of kaolinite, with the ultimate goal of estimating the electrical conductivity of kaolinite with a wide range of pore water conductivities (σ _w = 0.013–3.356 S/m) and porosities (n = 0.368–1.0). Therefore, the theoretical background of the electrical conductivity in soils was reviewed, and electrical conductivity measurements on kaolinite were performed using both slurry and consolidation tests in this study. The results of this study demonstrate that the variations of measured electrical conductivity (σ _mix) with n are debatable according to the values of σ _w, because a decrease in n results in both an increase in surface conduction (K _s) and a decrease in pore water conduction (K _w); this causes the relative magnitude of K _s compared to that of K _w to vary with σ _w and n. Consequently, this study develops the relation between the porosity-normalized K _s/K _w and 1/σ _w. Additionally, the surface conductivity of the tested kaolinite is back-calculated and compared with the previous relationship between K _s and zeta potential of kaolinite. The measured and estimated σ _mix values are compared with the varying pore water conductivity and porosity values.     

On the crystal chemistry and elastic behavior of a phlogopite 3<Emphasis Type="Italic">T</Emphasis>

The crystal chemistry and the elastic behavior under isothermal conditions up to 9 GPa of a natural, and extremely rare, 3T-phlogopite from Traversella (Valchiusella, Turin, Western Alps) [(K_0.99Na_0.05Ba_0.01)(Mg_2.60Al_0.20Fe _0.21 ²⁺ )[Si_2.71Al_1.29O₁₀](OH)₂, space group P3₁12, with a = 5.3167(4), c = 30.440(2) Å, and V = 745.16(9) ų] have been investigated by electron microprobe analysis in wavelength dispersion mode, single-crystal X-ray diffraction at 100 K, and in situ high-pressure synchrotron radiation powder diffraction (at room temperature) with a diamond anvil cell. The single-crystal refinement confirms the general structure features expected for trioctahedral micas, with the inter-layer site partially occupied by potassium and sodium, iron almost homogeneously distributed over the three independent octahedral sites, and the average bond distances of the two unique tetrahedra suggesting a disordered Si/Al-distribution (i.e., 〈T1-O〉 ~ 1.658 and 〈T2-O〉 ~ 1.656 Å). The location of the H-site confirms the orientation of the O–H vector nearly perpendicular to (0001). The refinement converged with R ₁(F) = 0.0382, 846 unique reflections with F _O > 4σ(F _O) and 61 refined parameters, and not significant residuals in the final difference-Fourier map of the electron density (+0.77/?0.37 e ^?/ų). The high-pressure experiments showed no phase transition within the pressure range investigated. The P–V data were fitted with a Murnaghan (M-EoS) and a third-order Birch-Murnaghan equation of state (BM-EoS), yielding: (1) M-EoS, V ₀ = 747.0(3) ų, K _T0 = 44.5(24) GPa, and K′ = 8.0(9); (2) BM-EoS, V ₀ = 747.0(3) ų, K _T0 = 42.8(29) GPa, and K′ = 9.9(17). A comparison between the elastic behavior in response to pressure observed in 1M- and 3T-phlogopite is made.     

Evaluating of the earthquake hazard parameters with Bayesian method for the different seismic source regions of the North Anatolian Fault Zone

The earthquake hazard parameters and earthquake occurrence probabilities are computed for the different regions of the North Anatolia Fault Zone (NAFZ) using Bayesian method. A homogenous earthquake catalog for M _S magnitude which is equal or larger than 4.0 is used for a time period between 1900 and 2015. Only two historical earthquakes (1766, M _S = 7. 3 and 1897, M _S = 7. 0) are included in Region 2 (Marmara Region) where a large earthquake is expected in the near future since no large earthquake has been observed for the instrumental period. In order to evaluate earthquake hazard parameters for next 5, 10, 20, 50, 100 years, M _max (maximum regional magnitude), β value, λ (seismic activity or density) are computed for the different regions of NAFZ. The computed M _max values are changed between 7.11 and 7.89. While the highest magnitude value is calculated in the Region 9 related to Tokat-Erzincan, the lowest value in the Region 10 including the eastern of Erzincan. The “quantiles” of “apparent” and “true” magnitudes of future time intervals of 5, 10, 20, 50, and 100 years are calculated for confidence limits of probability levels of 50, 70 and 90 % of the 10 different seismic source regions. The region between Tokat and Erzincan has earthquake hazard level according to the determined parameters. In this region the expected maximum earthquake size is 7.8 with 90 % occurrence probability in next 100 years. While the regional M _max value of Marmara Region is computed as 7.61, expected maximum earthquake size is 7.37 with 90 % occurrence probability in next 100 years.     

Shear wave velocity as function of cone penetration resistance and grain size for Holocene-age uncemented soils: a new perspective

For feasibility studies and preliminary design estimates, field measurements of shear wave velocity, V _s, may not be economically adequate and empirical correlations between V _s and more available penetration measurements such as cone penetration test, CPT, data turn out to be potentially valuable at least for initial evaluation of the small-strain stiffness of soils. These types of correlations between geophysical (Vs) and geotechnical (N-SPT, q _c-CPT) measurements are also of utmost importance where a great precision in the calculation of the deposit response is required such as in liquefaction evaluation or earthquake ground response analyses. In this study, the stress-normalized shear wave velocity V _s1 (in m/s) is defined as statistical functions of the normalized dimensionless resistance, Q _tn-CPT, and the mean effective diameter, D ₅₀ (in mm), using a data set of different uncemented soils of Holocene age accumulated at various sites in North America, Europe, and Asia. The V _s1–Q _tn data exhibit different trends with respect to grain sizes. For soils with mean grain size (D ₅₀) < 0.2 mm, the V _s1/Q _tn ^0.25 ratio undergoes a significant reduction with the increase in D ₅₀ of the soil. This trend is completely reversed with further increase in D ₅₀ (D ₅₀ > 0.2 mm). These results corroborate earlier results that stressed the use of different CPT-based correlations with different soil types, and those emphasized the need to impose particle-size limits on the validity of the majority of available correlations.     

Probabilistic seismic hazard assessment in the Constantine region,Northeast of Algeria

This paper presents a seismic hazard evaluation and develops an earthquake catalogue for the Constantine region over the period from 1357 to 2014. The study contributes to the improvement of seismic risk management by evaluating the seismic hazards in Northeast Algeria. A regional seismicity analysis was conducted based on reliable earthquake data obtained from various agencies (CRAAG, IGN, USGS and ISC). All magnitudes (M _l, m _b) and intensities (I ₀, I _MM, I _MSK and I _EMS) were converted to M _s magnitudes using the appropriate relationships. Earthquake hazard maps were created for the Constantine region. These maps were estimated in terms of spectral acceleration (SA) at periods of 0.1, 0.2, 0.5, 0.7, 0.9, 1.0, 1.5 and 2.0 s. Five seismogenic zones are proposed. This new method differs from the conventional method because it incorporates earthquake magnitude uncertainty and mixed datasets containing large historical events and recent data. The method can be used to estimate the b value of the Gutenberg-Richter relationship, annual activity rate λ(M) of an event and maximum possible magnitude M _max using incomplete and heterogeneous data files. In addition, an earthquake is considered a Poisson with an annual activity rate λ and with a doubly truncated exponential earthquake magnitude distribution. Map of seismic hazard and an earthquake catalogue, graphs and maps were created using geographic information systems (GIS), the Z-map code version 6 and Crisis software 2012.     

The Influence of Benthic Macrofauna on the Erodibility of Intertidal Sediments with Varying mud Content in Three New Zealand Estuaries

Fine sediment inputs can alter estuarine ecosystem structure and function. However, natural variations in the processes that regulate sediment transport make it difficult to predict their fate. In this study, sediments were sampled at different times (2011–2012) from 45 points across intertidal sandflat transects in three New Zealand estuaries (Whitford, Whangamata, and Kawhia) encompassing a wide range in mud (≤63 μm) content (0–56 %) and macrofaunal community structure. Using a core-based erosion measurement device (EROMES), we calculated three distinct measures of sediment erosion potential: erosion threshold (? _c ; N m^?2), erosion rate (ER; g m^?2 s^?1), and change in erosion rate with increasing bed shear stress (m _e ; g N^?1 s^?1). Collectively, these measures characterized surface (? _c and ER) and sub-surface (m _e ) erosion. Benthic macrofauna were grouped by functional traits (size and motility) and data pooled across estuaries to determine relationships between abiotic (mud content, mean grain size) and biotic (benthic macrofauna, microbial biomass) variables and erosion measures. Results indicated that small bioturbating macrofauna (predominantly freely motile species <5 mm in size) destabilized surface sediments, explaining 23 % of the variation in ? _c (p ≤ 0.01) and 59 % of the variation in ER (p ≤ 0.01). Alternatively, mud content and mean grain size cumulatively explained 61 % of the variation in m _e (p ≤ 0.01), where increasing mud and grain size stabilized sub-surface sediments. These results highlight that the importance of biotic and abiotic predictors vary with erosion stage and that functional group classifications are a useful way to determine the impact of benthic macrofauna on sediment erodibility across communities with different species composition.     

Numerical analysis under seismic loads of soils improvement with floating stone columns

Geotechnical Engineering has developed many methods for soil improvement so far. One of these methods is the stone column method. The structure of a stone column generally refers to partial change of suitable subsurface ground through a vertical column, poor stone layers which are completely pressed. In general terms, to improve bearing capacity of problematic soft and loose soil is implemented for the resolution of many problems such as consolidation and grounding problems, to ensure filling and splitting slope stability and liquefaction that results from a dynamic load such as earthquake. In this study, stone columns method is preferred as an improvement method, and especially load transfer mechanisms and bearing capacity of floating stone column are focused. The soil model, 32 m in width and 8 m in depth, used in this study is made through Plaxis 2D finite element program. The clay having 5° internal friction angle with different cohesion coefficients (c 10, c 15, c 20 kN/m²) are used in models. In addition, stone columns used for soil improvement are modeled at different internal friction angles (? 35°, ? 40°, ? 45°) and in different s/D ranges (s/D 2, s/D 3), stone column depths (B, 2B, 3B) and diameters (D 600 mm, D 800 mm, D 1000 mm). In the study, maximum acceleration (a _max = 1.785 m/s²) was used in order to determine the seismic coefficient used. In these soil models, as maximum acceleration, maximum east–west directional acceleration value of Van Muradiye earthquake that took place in October 23, 2011 was used. As a result, it was determined that the stone column increased the bearing capacity of the soil. In addition, it is observed that the bearing capacity of soft clay soil which has been improved through stone column with both static and earthquake load effect increases as a result of increase in the diameter and depth of the stone column and decreases as a result of the increase in the ranges of stone column. In the conducted study, the bearing capacity of the soil models, which were improved with stone column without earthquake force effect, was calculated as 1.01–3.5 times more on the average, compared to the bearing capacity of the soil models without stone column. On the other hand, the bearing capacity of the soil models with stone columns, which are under the effect of earthquake force, was calculated as 1.02–3.7 times more compared to the bearing capacity of the soil models without stone column.     

Effect of 2-D Random Field Discretization on Failure Probability and Failure Mechanism in Probabilistic Slope Stability

This paper investigates, using the random field theory and Monte Carlo simulation, the effects of random field discretization on failure probability, p _f, and failure mechanism of cohesive soil slope stability. The spatial sizes of the discretized elements in random field Δx, Δy in horizontal and vertical directions, respectively, are assigned a series of combinational values in order to model the discretization accuracy. The p _f of deterministic critical slip surface (DCSS) and that of the slope system both are analyzed. The numerical simulation results have demonstrated that both the ratios of Δy/λ _y (λ _y  = scale of fluctuation in vertical direction) and Δx/λ _x (λ _x  = scale of fluctuation in horizontal direction) contribute in a similar manner to the accuracy of p _f of DCSS. The effect of random field discretization on the p _f can be negligible if both the ratios of Δx/λ _x and Δy/λ _y are no greater than 0.1. The normalized discrepancy tends to increase at a linear rate with Δy/λ _y when Δx/λ _x is larger than 0.1, and vice versa for p _f of DCSS. The random field discretization tends to have more considerable influence on the p _f of DCSS than on that of the slope system. The variation of p _f versus λ _x and λ _y may exhibit opposite trends for the cases where the limit state functions of slope failure are defined on DCSS and on the slope system as well. Finally, the p _f of slope system converges in a more rapid manner to that of DCSS than the failure mechanism does to DCSS as the spatial variability of soil property grows from significant to negligible.     

Properties of the aftershocks sequences of the 2000 and 2002 earthquakes in Aksehir-Afyon graben,west-central Anatolia,Turkey

The b-value of the Gutenberg–Richter’s frequency–magnitude relation and the p-value of the modified Omori law, which describes the decay rate of aftershock activity, were investigated for more than 500 aftershocks in the Aksehir-Afyon graben (AAG) following the 15 December 2000 Sultandagi–Aksehir and the 3 February 2002 Çay–Eber and Çobanlar earthquakes. We used the Kandilli Observatory’s catalog, which contains records of aftershocks with magnitudes ≥2.5. For the Çobanlar earthquake, the estimated b-values for three aftershock sequences are in the range 0.34 ≤  b ≤ 2.85, with the exception of the one that occurred during the first hour (4.77), while the obtained p-values are in the range 0.44 ≤ p ≤ 1.77. The aftershocks of the Sultandagi earthquake have a high p-value, indicating fast decay of the aftershock activity. A regular increase of b can be observed, with b < 1.0 after 0.208 days for the Çay–Eber earthquake. A systematic and similar increase and decrease pattern exists for the b- and p-values of the Çobanlar earthquakes during the first 5 days.     

Elastic behavior of vanadinite,Pb<Subscript>10</Subscript>(VO<Subscript>4</Subscript>)<Subscript>6</Subscript>Cl<Subscript>2</Subscript>, a microporous non-zeolitic mineral

The high-pressure behavior of a vanadinite (Pb₁₀(VO₄)₆Cl₂, a = b = 10.3254(5), c = 7.3450(4) Å, space group P6₃/m), a natural microporous mineral, has been investigated using in-situ HP-synchrotron X-ray powder diffraction up to 7.67 GPa with a diamond anvil cell under hydrostatic conditions. No phase transition has been observed within the pressure range investigated. Axial and volume isothermal Equations of State (EoS) of vanadinite were determined. Fitting the P–V data with a third-order Birch-Murnaghan (BM) EoS, using the data weighted by the uncertainties in P and V, we obtained: V ₀ = 681(1) Å³, K ₀ = 41(5) GPa, and K′ = 12.5(2.5). The evolution of the lattice constants with P shows a strong anisotropic compression pattern. The axial bulk moduli were calculated with a third-order “linearized” BM-EoS. The EoS parameters are: a ₀ = 10.3302(2) Å, K ₀(a) = 35(2) GPa and K′(a) = 10(1) for the a-axis; c ₀ = 7.3520(3) Å, K ₀(c) = 98(4) GPa, and K′(c) = 9(2) for the c-axis (K ₀(a):K ₀(c) = 1:2.80). Axial and volume Eulerian-finite strain (fe) at different normalized stress (Fe) were calculated. The weighted linear regression through the data points yields the following intercept values: Fe _a (0) = 35(2) GPa for the a-axis, Fe _c (0) = 98(4) GPa for the c-axis and Fe _V (0) = 45(2) GPa for the unit-cell volume. The slope of the regression lines gives rise to K′ values of 10(1) for the a-axis, 9(2) for the c-axis and 11(1) for the unit cell-volume. A comparison between the HP-elastic response of vanadinite and the iso-structural apatite is carried out. The possible reasons of the elastic anisotropy are discussed.     

Geological Characteristics and Stability Evaluation of Wanjia Middle School Slope in Wenchuan Earthquake Area

To investigate the formation mechanism and the stability of Wanjia middle school slope in Wenchuan Earthquake Area, the macroscopic geological characteristics and the failure process of the landslide are researched by engineering geology analysis method, limit equilibrium method, and finit element method. The results show that after the Wenchuan Earthquake, retaining walls, houses and other infrastructure on the foot of Wanjia middle school slope were severely destroyed, 10 cm wide tension fracture appeared at the trailing edge of the slope. Wanjia middle school slope is a type of medium-sized soil landslide. The area of the deformation body is about 19,314 m², the total volume of the deformation body is about 23 × 10⁴ m³. There may be two potential sliding surfaces in the unstable slope: shallow and deep landslide. The analysis results of the limit equilibrium method and the finite element method show that: under dead weight, dead weight + rainstorm, dead weight + earthquake conditions, the plastic zone occurs mainly at the middle part or the trailing edge of the slope, and it doesn’t fully cut through the deep landslide body, so the deep landslide is stable. However, under rainstorm or earthquake conditions, the plastic zone almost completely cut thorough the shallow landslide body, it shows that the shallow landslide is in the understable–basic stable state. It is found that the results of finite element method is concordant with the results of the limit equilibrium method (F _s = 1.06–1.29, the shallow landslide is in the basic stable–stable state). The calculation results show that shallow landslides are likely to occur in Wanjia middle school slope during a rainstorm or an earthquake, so monitoring and control of the slope should be strengthened. The shallow landslide should be managed by some measures, such as anti slide pile retaining structures and drainage works, and the dangerous rock bodies on the slope surface should be cleaned up.     

Seismic source characteristics in Kachchh and Saurashtra regions of Western India: <Emphasis Type="Italic">b</Emphasis>-value and fractal dimension mapping of aftershock sequences

Seismic source characteristics in the Kachchh rift basin and Saurashtra horst tectonic blocks in the stable continental region (SCR) of western peninsular India are studied using the earthquake catalog data for the period 2006–2011 recorded by a 52-station broadband seismic network known as Gujarat State Network (GSNet) running by Institute of Seismological Research (ISR), Gujarat. These data are mainly the aftershock sequences of three mainshocks, the 2001 Bhuj earthquake (M _w 7.7) in the Kachchh rift basin, and the 2007 and 2011 Talala earthquakes (M _w ≥ 5.0) in the Saurashtra horst. Two important seismological parameters, the frequency–magnitude relation (b-value) and the fractal correlation dimension (D _c) of the hypocenters, are estimated. The b-value and the D _c maps indicate a difference in seismic characteristics of these two tectonic regions. The average b-value in Kachchh region is 1.2 ± 0.05 and that in the Saurashtra region 0.7 ± 0.04. The average D _c in Kachchh is 2.64 ± 0.01 and in Saurashtra 2.46 ± 0.01. The hypocenters in Kachchh rift basin cluster at a depth range 20–35 km and that in Saurashtra at 5–10 km. The b-value and D _c cross sections image the seismogenic structures that shed new light on seismotectonics of these two tectonic regions. The mainshock sources at depth are identified as lower b-value or stressed zones at the fault end. Crustal heterogeneities are well reflected in the maps as well as in the cross sections. We also find a positive correlation between b- and D _c-values in both the tectonic regions.     

High-pressure displacive phase transition of a natural Mg-rich pigeonite

A high-pressure single-crystal X-ray diffraction study has been carried out on a P2₁/c natural Mg-rich pigeonite sample with composition ca. Wo₆En₇₆Fs₁₈ using a diamond anvil-cell. The unit-cell parameters were determined at 14 different pressures to 7.14 GPa. The sudden disappearance of the b-type reflections (h + k = odd) and a strong discontinuity (about 2.8%) in the unit-cell volume indicated a first-order P2₁/c–C2/c phase transition between 4.66 and 4.88 GPa. The P(V) data of the P2₁/c phase were fitted to 4.66 GPa by a third-order Birch–Murnaghan equation of state (BM3 EoS), whereas the limited number of experimental data collected within the C2/c phase between 4.88 and 7.14 GPa were fitted using the same equation of state but with K′ constrained to the value obtained for the P2₁/c fitting. The equation of state coefficients are V ₀ = 424.66(6) ų, K _T0 = 104(2) GPa and K′ = 8(1) for the P2₁/c phase, and V ₀ = 423.6(1) ų, K _T0 = 112.4(8) GPa, and K′ fixed to 8(1) for the C2/c phase. The axial moduli for a, b, and c for the P2₁/c phase were obtained using also a BM3-EoS, while for the C2/c phase only a linear calculation could be performed, and therefore the same approach was applied for comparison also to the P2₁/c phase. In general the C2/c phase exhibits axial compressibilities (β _c > β _a >> β _b) lower than those of the P2₁/c phase (β _b > β _c ≈ β _a; similar to those found in previous studies in clinopyroxenes and orthopyroxenes). The lower compressibility of the C2/c phase compared with that of the P2₁/c could be ascribed to the greater stiffness along the b direction. A previously published relationship between P _c and M2 average cation radius (i.r.) has been updated using all the literature data on P2₁/c clinopyroxene containing large cations at M2 site and our new data. The following weighted regression was obtained: P _c (GPa) = 26(4) ? 28(5) ×  i.r (Å), R ² = 0.97. This improved equation can be used to predict the critical pressure of natural P2₁/c clinopyroxene samples just knowing the composition at M2 site.