各向异性软土基坑抗隆起稳定极限平衡分析

土体的各向异性与沉积环境和应力历史密切相关。以能够考虑土体特性的各向异性不排水抗剪强度公式为基础,利用圆弧滑动方法,并考虑破坏面上强度的折减（ ）,推导各向异性软土基坑抗隆起稳定极限平衡解。对软土基坑抗隆起稳定进行研究,系统分析了开挖深度H、围护结构插入深度D、最下道支撑与坑底相对距离He/H等几何因素及土体有效内摩擦角 、各向异性强度比k及强度沿深度非均质变化的斜率 等强度参数的影响。最后结合工程实例,验证了公式的适用性,为软土基坑工程抗隆起分析提供了新的思路。     

Undrained Shear Strength in Cohesive Soils Estimated by Directional Modes of In-Situ Shear Wave Velocity

The estimated undrained shear strength (s_u) is often not a unique value because it can be evaluated by various test types and/or procedures, such as different failure modes, shear strain rates, and boundary conditions. This study explores (1) the relationship between reference undrained shear strength and in situ shear wave velocity in terms of the effective overburden stress, and (2) the independent relationships to evaluate the undrained shear strength with special consideration of different directional and polarization modes (VH, HV, HH shear waves), which has not been reported. This evaluation is done via a worldwide database compiled from 43 well-documented geotechnical test sites associated with soft ground. Finally, new correlation models are proposed to estimate the undrained shear strength based on the in situ shear wave velocity as well as the plasticity index or the overconsolidation ratio. The application of the shear wave velocity–undrained shear strength relation is illustrated through two independent case studies. The proposed relationships are expected to contribute to reasonable estimates of undrained shear strength as well as offer practical guidance on even extrapolation beyond the data that is available to geotechnical engineers.     

Interpretation of shear strength uncertainty and reliability analyses of slopes

The variance of the friction angle or friction coefficient (tan ?) is often considered in geotechnical reliability analyses, which implies that the variance of the shear strength as defined by a Mohr-Coulomb envelope increases as the normal stress on the shearing surface increases. However, shear strength data sometimes has approximately constant variance, and most simple regression techniques assume constant variance. Four effective stress shear strength data sets are evaluated using both the constant variance (homoscedastic) and constant coefficient of variation (heteroscedastic) interpretations. The impact of the variance interpretation on slope stability is evaluated using infinite slope, homogeneous dam, and zoned dam examples. For relatively shallow infinite slope surfaces, the reliability index for the heteroscedastic interpretation of shear strength variance was about twice the reliability index obtained using the homoscedastic approach. In the dam examples, the difference in the reliability indices resulting from the heteroscedastic and homoscedastic interpretations was about one, indicating a tenfold increase in the probability of failure. The typical assumption of constant coefficient of variation of shear strength may result in unconservative estimates of the reliability of shallow failure surfaces and overly conservative results for deeper failure surfaces.     

Garnet–orthopyroxene (GOPX) geothermometer: a comparative study

The garnet–orthopyroxene pairs are commonly found in the assemblages of basic granulites/charnockite and hence are suitable for estimating equilibrium temperature of these metamorphic rocks. At present, there are many calibrations of garnet–orthopyroxene thermometer that may confuse geologists in choosing a reliable thermometer. To test the accuracy of the garnet–orthopyroxene thermometers, we have applied 14 models formulated by a number of workers since 1980 to date. We have collated 51 samples from the literature all over the world, which has been processed through the “Gt-Opx.EXE” software. Based on the present study, we have identified a set of the best among all the 14 models which were considered under this comparative study. We have concluded that the five garnet–orthopyroxene (Gt-Opx) thermometers are the most valid and reliable of this kind of thermometer: Aranovich and Berman (Am Mineral 82:345–353, 1997), Raith et al. (Earth Sci 73:211–244, 1983), Harley (Contrib Mineral Petrol 86:359–373, 1984), Nimis and Grütter (Contrib Mineral Petrol 159:411–427, 2010), and Sen and Bhattacharya (Contrib Mineral Petrol 88:64–71, 1984).     

基于变固结系数的静压桩周土体力学特性研究

以静压沉桩后桩周土体的应力状态为初始条件,根据桩周土体孔隙比、渗透系数和有效应力之间的相关性,在考虑固结系数随固结时间变化的条件下改进了轴对称条件下的太沙基固结控制方程。随后,采用分离变量法和离散化分析推导得出了桩周超孔隙水压力消散的半解析半数值解,并将解答与实测数据进行对比验证。在此基础上,采用空间滑动面理论改进的修正剑桥模型（SMP-MCC）定义土体三维不排水抗剪强度,研究了静压桩周土体强度、剪切模量随固结时间的变化规律。研究结果表明：由于解答考虑了固结系数随固结时间的变化,因而与实测结果吻合良好;土体压缩指数与渗透指数之比对土体固结系数和孔压消散速率影响较大;当土体压缩指数与渗透指数之比为1时,土体固结系数保持不变,解答退化为经典的太沙基轴对称固结方程;土体强度和剪切模量随固结时间的增长而逐步增加,固结完成后其值超越了土体原位强度和原位剪切模量。     

An enhanced SMA based SCS-CN inspired model for watershed runoff prediction

Incorporation of initial soil moisture (V ₀) in the Soil Conservation Service Curve Number (SCS-CN) methodology helps to avoid the sudden jumps in Curve Number (CN) and, in turn, in computed runoff. It invoked the development of an enhanced (yet simple) Soil Moisture Accounting (SMA) procedure-based-SCS-CN inspired model, by incorporating initial moisture (V ₀). Its performance is tested using a dataset of 152 small to large watersheds of USDA (total 38,169 storm events), and compared with original SCS-CN method, Mishra and Singh (Acta Geophys Polon 50(3):457–477, 2002), Michel et al. (Water Resour Res 41(2):W02011, 2005) and Singh et al. (Water Resour Manag 29(11): 4111–4127, 2015) model using four statistical indices (RMSE, R ², PBIAS and NSE) and rank grading system (RGS). The proposed model scores highest (= 691 marks out of maximum 2280 marks) (Rank I) followed by Singh et al. (Water Resour Manag 29(11):4111–4127, 2015) model with 642 marks (Rank II), Michel et al. (Water Resour Res 41(2):W02011, 2005) model with 376 marks (Rank III) and Mishra and Singh model with 362 marks (= Rank IV). The original SCS-CN model, however, performs the poorest of all with 209 marks (Rank V).     

Qanat: a resource or a hazard? Rebuttal to “Qanat is not a hazard” by Parise (Environ Earth Sci 2016 75: 1476)

A reply essay is presented on the rebuttal article by Parise (Environ Earth Sci 75(23):1476, 2016) suggesting that qanat is not a hazard. It is presented as a refutation on the paper by Abbasnejad et al. (Environ Earth Sci 75:1306, 2016) in which the authors have explained the etiology, the impacts and remedies of a qanat hazard in Iranian urban areas. Since the majority of qanats in Iranian urban areas are abandoned, useless and threatening, according to definition and in comparison with similar features introduced as a hazard, they are considered as a hazard too. However, this does not mean, and Abbasnejad et al. (in Environ Earth Sci 75:1306, 2016) have not claimed, that all qanats are hazardous. In addition, the authors who have studied qanats in Iranian urban areas, before Abbasnejad et al. (in Environ Earth Sci 75:1306, 2016), have also considered qanats as a hazard.     

The vulnerability of the European air traffic network to spatial hazards

The 2010 eruption of the Eyjafjallajökull volcano had a devastating effect on the European air traffic network, preventing air travel throughout most of Europe for 6 days (Oroian in ProEnvironment 3:5–8, 2010). The severity of the disruption was surprising as previous research suggests that this type of network should be tolerant to random hazard (Albert et al. in Nature 406(6794):378–382, 2000; Strogatz in Nature 410(6825):268–276, 2001). The source of this hazard tolerance lies in the degree distribution of the network which, for many real-world networks, has been shown to follow a power law (Albert et al. in Nature 401(6749):130–131, 1999; Albert et al. in Nature 406(6794):378–382, 2000). In this paper, we demonstrate that the ash cloud was unexpectedly disruptive because it was spatially coherent rather than uniformly random. We analyse the spatial dependence in air traffic networks and demonstrate how the combination of their geographical distribution and their network architectures jeopardises their inherent hazard tolerance.     

Miocene echinoids from the Sadat area,south Gebel Ataqa,NW Gulf of Suez,Egypt: systematics,palaeobiogeography

The echinoid fauna from the Miocene sedimentary succession cropping out south Wadi Tweirig, and Wadi Hommath, south Gebel Ataqa, NW Gulf of Suez, has been examined with the aim to known their stratigraphic and paleogeographic distribution. The Miocene succession includes two formations: Sadat Formation, unconformably overlying the middle/upper Eocene rocks at the base and Hommath Formation at the top. Twenty-eight echinoid species (8 regular and 20 irregular) belonging to 18 genera, 13 families, and 7 orders have been identified, systematically described, and illustrated in this work. Eleven species are recorded for the first time from Egypt: ten of these came from the Hommath Formation (Schizechinus cf. serresii Desor (1856), Schizechinus pentagonus Kier 1972, Clypeaster cf. martini des Moulins 1837, Scutella checchiae occidentalis Desio 1934, Scutella melitensis Airaghi 1902, Echinodiscus desori Duncan and Sladen 1883, Echinolampas cf. zeitensis Fourtau 1920, Schizaster lovisatoi Cotteau 1895, Agassizia (Agassizia) powersi Kier 1972, and Hemipatagus ocellatus Defrance (1827)), and one from the Sadat Formation (Clypeaster campanulatus Schlotheim (1820)). The identified fauna shows a strong affinity with the Mediterranean bio-province.     

Discussion of “From saturated to unsaturated conditions and vice versa” by Martí Lloret-Cabot,Simon J. Wheeler,Jubert A. Pineda,Enrique Romero,and Daichao Sheng (https://doi.org/10.1007/s11440-017-0577-6)

The authors (Lloret-Cabot et al. in Acta Geotech 1–23, 2017) applied the glasgow coupled model (GCM), originally proposed by Wheeler et al. (Géotechnique 53(1):41–54, 2003), to the simulation of several experimental tests that involve transition between saturated and unsaturated states. The authors show qualitatively, but not quantitatively, predictions of GCM for shrinkage during air drying of normally consolidated samples (Fig. 13) under low mechanical stress conditions, without presenting the material parameters. The discussers, who have worked with GCM to model multilayer deposition of tailings/soft soils (Qi in Numerical investigation for slope stability of expansive soils and large strain consolidation of soft soils. Doctoral dissertation, University of Ottawa, 2017; Qi et al. in J Geotech Geoenviron Eng 143(7):04017018, 2017, Qi et al. J Geotech Geoenviron Eng 143(7):04017019, 2017), have made quantitative predictions of similar cases. Satisfactory simulations of such cases using GCM are sensitive to the selection of the coupling parameters k₁ and k₂. By considering an alternative analytical form of GCM, an analytical procedure can be derived for calibrating the coupling parameters for problems involving virgin drying.     

SWAT hydrological model as a DaaS cloud service

Earth Science community depends on the exploration, analysis and reprocessing of high volumes of data as well as the modeling and simulation of complex coupled systems on multiple scales. The main aim of this article is to introduce a new hydrological modeling service based on the Soil and Water Assessment Tool (SWAT) (Arnold et al. J American Water Resour Assoc 34(1), 73–89, 1998 ; Arnold and Fohrer Hydrol Process 19(3), 563–572, 2005) model using high efficiency, resource sharing and low cost cloud computing resources (Astsatryan et al. International Journal of Scientific & Engineering Research 1(1), 1130–1133, 2014). Such a Desktop as a Service (DaaS) approach allowing users to work from anywhere, and gives centralized desktop management and great performance. Within the Spatial Data Infrastructure (SDI) and cloud platform, the DaaS service gives secure access to the model and a centralized data storage to get a SWAT model input. The article illustrates the analyses of the implementation of the SWAT model for the Sotk watershed of Lake Sevan in Armenia (Sargsyan 2007).     

Effects of anisotropy in correlation structure on the stability of an undrained clay slope

Slopes are mainly naturally occurred deposits, so slope stability is highly affected by inherent uncertainty. In this paper, the influence of heterogeneity of undrained shear strength on the performance of a clay slope is investigated. A numerical procedure for a probabilistic slope stability analysis based on a Monte Carlo simulation that considers the spatial variability of the soil properties is presented to assess the influence of randomly distributed undrained shear strength and to compute reliability as a function of safety factor. In the proposed method, commercially available finite difference numerical code FLAC 5.0 is merged with random field theory. The results obtained in this study are useful to understand the effect of undrained shear strength variations in slope stability analysis under different slope conditions and material properties. Coefficient of variation and heterogeneity anisotropy of undrained shear strength were proven to have significant effect on the reliability of safety factor calculations. However, it is shown that anisotropy of the heterogeneity has a dual effect on reliability index depending on the level of safety factor adopted.     

Effect of gravel content on creep behavior of clayey soil at residual state: implication for its role in slow-moving landslides

Soil creep is mostly manifested in slow-moving landslides. It is often the case for active slow-moving landslide with slip zone comprising clayey soil where creep would develop in the residual condition. If gravel is presented in clayey soil, this will have considerable impact on creep behavior of clayey soil. However, knowledge about creep behavior of the clayey soil containing gravel particles is scarce. This paper discusses creep behavior of natural clayey soil with gravel at residual state through a series of creep shear tests. Soil samples for this testing program were collected from the slip zones of two large slow-moving landslides in China. The collected soil samples consisted of clayey soil containing various amounts of gravel particles. The test results show that the soil specimens underwent two different creep patterns. An attenuating creep pattern was observed when the soil specimens were subjected to creep stress less than the residual strength, and a creep with increasing strain rate, or a nonattenuating creep pattern, with no evident secondary creep was noted when creep stress was intestinally increased to a level slightly greater than the residual strength. The creep patterns of clayey soil with gravel at the residual state observed in this study were noted to be very consistent with those of the gravel-free clayey soil reported in Bhat et al. (Int J Geomater 1(1):39–43, 2011, Nat Hazards 69(3):2161–2178, 2013) and Di Miao et al. (Eng Geol 162:53–66, 2013). This gives an indication that creep patterns of clayey soils with and without gravel at the residual state are essentially the same, or in other words, that the presence of gravel does not change the creep pattern of clayey soil. However, the test results in this study illustrate that the presence of gravel does have a notable effect on creep behavior of clayey soil. Specifically, the creep stress leading to creep failure of clayey soil and the minimum ratio of the creep stress to residual strength (RCSR) increased with gravel content, and the displacement until the tertiary creep was also larger in samples containing more gravel particles. It is postulated that creep behavior of clayey soil at the residual state in this study and its relationship with gravel content may be related to strength recovery and crushing of gravel-sized particles during creep.     

Seismic risk assessment and mapping at different levels

In 2012, the damage costs of floods in Russia amounted to about €300m, and these floods have caused nearly 200 fatalities (Kotlyakov et al. in Reg Res Rus 3(1):32–39, 2013). Risk assessment is one of the most pressing scientific topics in Russia, but most of the works are devoted to natural hazards assessment. The purpose of this work is to estimate the influence of hazardous hydrological phenomena on society. The field research was conducted in the Slavyansk municipal district in the Krasnodar region (the south-western part of Russia), which is a highly populated coastal territory with a high frequency of hazardous hydrological events. Modified methods of the Ministry of the Russian Federation for Affairs for Civil Defence, Emergencies and Elimination of Consequences of Natural Disasters (EMERCOM) were used for potential economic damage calculation. The paper did not only focus on direct, tangible risks, but also included social risk (i.e. risk to life and health). Social vulnerability has been calculated directly as a percentage of vulnerable people, estimated in opinion polls, while in many recent papers the social vulnerability index was calculated as a combination of several statistical indicators. The resulting percentage of vulnerable people was converted to numbers of potential victims. Finally, the social risk was expressed by financial indicators in terms of the cost of the value of statistical life lost (Mrozek and Taylor in J Policy Anal Manag 21(2):253–270, 2002; Viscusi and Aldy in J Risk Uncertain 27(1):5–76, 2003). Social risk can be underestimated in comparison with economic risk because of a low “value of life” in Russia (no life insurance, neglecting of basic safety rules, etc.) (Guriev in Myths of economics, Alpina Business Books, Moscow, 2009).     

Multiscale model reduction for shale gas transport in fractured media

In this paper, we develop a multiscale model reduction technique that describes shale gas transport in fractured media. Due to the pore-scale heterogeneities and processes, we use upscaled models to describe the matrix. We follow our previous work (Akkutlu et al. Transp. Porous Media 107(1), 235–260, 2015), where we derived an upscaled model in the form of generalized nonlinear diffusion model to describe the effects of kerogen. To model the interaction between the matrix and the fractures, we use Generalized Multiscale Finite Element Method (Efendiev et al. J. Comput. Phys. 251, 116–135, 2013, 2015). In this approach, the matrix and the fracture interaction is modeled via local multiscale basis functions. In Efendiev et al. (2015), we developed the GMsFEM and applied for linear flows with horizontal or vertical fracture orientations aligned with a Cartesian fine grid. The approach in Efendiev et al. (2015) does not allow handling arbitrary fracture distributions. In this paper, we (1) consider arbitrary fracture distributions on an unstructured grid; (2) develop GMsFEM for nonlinear flows; and (3) develop online basis function strategies to adaptively improve the convergence. The number of multiscale basis functions in each coarse region represents the degrees of freedom needed to achieve a certain error threshold. Our approach is adaptive in a sense that the multiscale basis functions can be added in the regions of interest. Numerical results for two-dimensional problem are presented to demonstrate the efficiency of proposed approach.     

Separate-ice frost heave model for one-dimensional soil freezing process

For one-dimensional soil freezing process, a separate-ice frost heave model is established, and the coupled process of heat transfer, fluid flow and stress development is considered in the model. First, a coupled heat–fluid–stress model describing the growth of a single ice lens is developed by extending the coupled heat–fluid model presented by Zhou and Zhou (Can Geotech J 49(6):686–693, 2012). Second, the mechanism for the formation of a new ice lens in the frozen fringe is studied, and we indicate that if the total vertical disjoining pressure at certain place exceeds the sum of the external pressure and the critical pressure, a new ice lens will emerge. By combining the growth model of a single ice lens and the criterion for the formation of a new ice lens, the separate-ice frost heave model is then established. The difference between the separate-ice model and the rigid-ice model is explained, and the relations for different mathematical models which describe the soil freezing process are also discussed. Numerical analysis of the separate-ice model is conducted using the finite volume method. The freezing tests for Devon silt under no external pressure and Xuzhou silty clay under a constant external pressure are applied to verify the computational results. The consistence between the calculation and the observation validates the separate-ice frost heave model.     

Anisotropic mechanical properties of zircon and the effect of radiation damage

This study provides new insights into the relationship between radiation-dose-dependent structural damage due to natural U and Th impurities and the anisotropic mechanical properties (Poisson’s ratio, elastic modulus and hardness) of zircon. Natural zircon samples from Sri Lanka (see Muarakami et al. in Am Mineral 76:1510–1532, 1991) and synthetic samples, covering a dose range of zero up to 6.8 × 10¹⁸ α-decays/g, have been studied by nanoindentation. Measurements along the [100] crystallographic direction and calculations, based on elastic stiffness constants determined by Özkan (J Appl Phys 47:4772–4779, 1976), revealed a general radiation-induced decrease in stiffness (~54 %) and hardness (~48 %) and an increase in the Poisson’s ratio (~54 %) with increasing dose. Additional indentations on selected samples along the [001] allowed one to follow the amorphization process to the point that the mechanical properties are isotropic. This work shows that the radiation-dose-dependent changes of the mechanical properties of zircon can be directly correlated with the amorphous fraction as determined by previous investigations with local and global probes (Ríos et al. in J Phys Condens Matter 12:2401–2412, 2000a; Farnan and Salje in J Appl Phys 89:2084–2090, 2001; Zhang and Salje in J Phys Condens Matter 13:3057–3071, 2001). The excellent agreement, revealed by the different methods, indicates a large influence of structural and even local phenomena on the macroscopic mechanical properties. Therefore, this study indicates the importance of acquiring better knowledge about the mechanical long-term stability of radiation-damaged materials.     

Micro-mechanical failure analysis of wet granular matter

We employ a novel fluid–particle model to study the shearing behavior of granular soils under different saturation levels, ranging from the dry material via the capillary bridge regime to higher saturation levels with percolating clusters. The full complexity of possible liquid morphologies Scheel et al. (Nat Mater 7(3):189–193, 2008. doi: 10.1038/nmat2117) is taken into account, implying the formation of isolated arbitrary-sized liquid clusters with individual Laplace pressures that evolve by liquid exchange via films on the grain surface Melnikov et al. (Phys Rev E 92(022):206, 2015. doi: 10.1103/PhysRevE.92.022206). Liquid clusters can grow in size, shrink, merge and split, depending on local conditions, changes of accessible liquid and the pore space morphology determined by the granular phase. This phase is represented by a discrete particle model based on contact dynamics Brendel et al. (Contact dynamics for beginners. Wiley-VCH, Weinheim, 2005. doi: 10.1002/352760362X.ch14), where capillary forces exerted from a liquid phase add to the motion of spherical particles. We study the macroscopic response of the system due to an external compression force at various liquid contents with the help of triaxial shear tests. Additionally, the change in liquid cluster distributions during the compression due to the deformation of the pore space is evaluated close to the critical load.     

The formation of porphyry copper deposits

Copper is a moderately incompatible chalcophile element. Its behavior is strongly controlled by sulfides. The speciation of sulfur is controlled by oxygen fugacity. Therefore, porphyry Cu deposits are usually oxidized (with oxygen fugacities > ΔFMQ +2) (Mungall 2002; Sun et al. 2015). The problem is that while most of the magmas at convergent margins are highly oxidized, porphyry Cu deposits are very rare, suggesting that high oxygen fugacity alone is not sufficient. Partial melting of mantle peridotite even at very high oxygen fugacities forms arc magmas with initial Cu contents too low to form porphyry Cu deposits directly (Lee et al. 2012; Wilkinson 2013). Here we show that partial melting of subducted young oceanic slabs at high oxygen fugacity (>ΔFMQ +2) may form magmas with initial Cu contents up to >500 ppm, favorable for porphyry mineralization. Pre-enrichment of Cu through sulfide saturation and accumulation is not necessarily beneficial to porphyry Cu mineralization. In contrast, re-melting of porphyritic hydrothermal sulfide associated with iron oxides may have major contributions to porphyry deposits. Thick overriding continental crust reduces the “leakage” of hydrothermal fluids, thereby promoting porphyry mineralization. Nevertheless, it is also more difficult for ore forming fluids to penetrate the thick continental crust to reach the depths of 2–4 km where porphyry deposits form.