Simplified multi-block constitutive model predicting earthquake-induced landslide triggering and displacement along slip surfaces of saturated sand

Slopes consisting of saturated sand have recently moved down-slope tens or hundreds of meters under the action of earthquakes. This paper presents a simplified but accurate method predicting the triggering and displacement of such landslides. For this purpose, a simplified constitutive model simulating soil response of saturated sands along slip surfaces is proposed and validated. Then, this constitutive model is coupled with the multi-block sliding system model to predict the triggering and displacement of such slides. The multi-block model considers a general mass sliding on a trajectory which consists of n linear segments. The steps needed to apply this method are described in detail. The method was applied successfully to predict the triggering, the motion and the final configuration of the well-documented (a) Higashi Takezawa, (b) Donghekou and (c) Nikawa earthquake-induced slides.     

Three-dimensional active earth pressure from cohesive backfills with tensile strength cutoff

Most of previous analyses on the active earth pressure were performed in two-dimensional cases using the Mohr-Coulomb (M-C) failure function to describe the soil strength. However, all failures of retained slopes indicate a somewhat three-dimensional (3D) feature, and the M-C function is found to overestimate the tensile strength of cohesive soil. In this work, a kinematic limit analysis–based approach is developed for computing the 3D active earth pressure resulting from cohesive backfills. The concept of tensile strength cutoff is adopted to implement the reduction or elimination of tensile strength from the strength envelope. An extended 3D horn failure mechanism that is associated with the modified strength envelope is developed to characterize the collapse of retained slopes. The resultant of active earth pressure is evaluated from the work rate balance equation and expressed as an unfactored coefficient. The obtained results indicate that less support provided by the wall is required when allowing the existence of soil cohesion and 3D effects and that eliminating the tensile strength can observably increase the active earth pressure, especially for the backfill with a great level of cohesion.     

Experimental determination of REE partition coefficients between zircon,garnet and melt: a key to understanding high‐T crustal processes

The partitioning of rare earth elements (REE) between zircon, garnet and silicate melt was determined using synthetic compositions designed to represent partial melts formed in the lower crust during anatexis. The experiments, performed using internally heated gas pressure vessels at 7 kbar and 900–1000 °C, represent equilibrium partitioning of the middle to heavy REE between zircon and garnet during high‐grade metamorphism in the mid to lower crust. The D_REE (zircon/garnet) values show a clear partitioning signature close to unity from Gd to Lu. Because the light REE have low concentrations in both minerals, values are calculated from strain modelling of the middle to heavy REE experimental data; these results show that zircon is favoured over garnet by up to two orders of magnitude. The resulting general concave‐up shape to the partitioning pattern across the REE reflects the preferential incorporation of middle REE into garnet, with D_Gd (zircon/garnet) ranging from 0.7 to 1.1, D_Ho (zircon/garnet) from 0.4 to 0.7 and D_Lu (zircon/garnet) from 0.6 to 1.3. There is no significant temperature dependence in the zircon–garnet REE partitioning at 7 kbar and 900–1000 °C, suggesting that these values can be applied to the interpretation of zircon–garnet equilibrium and timing relationships in the ultrahigh‐T metamorphism of low‐Ca pelitic and aluminous granulites.     

Numerical simulation of the response of a reinforced wall to a high speed train passage

The dynamic response of a mechanically stabilized earth wall to the passing of a high‐speed train is modelled using the finite element method. A three‐dimensional analysis is carried out, using a specific framework that allows performing the analysis with a moderate computational effort. In the first place, a so‐called multiphase approach is used to take into account the reinforcing strips. The moving load is taken into account by performing the calculation in a mobile referential using the properties of symmetry of the train cars and a simplifying assumption of periodicity for the whole train. We also assume a steady state. A partial validation of the approach is obtained by means of a comparison with an analytical solution. The quick increase in displacements induced by the train passing when the speed comes close to the celerity of Rayleigh waves clearly appears in the results. The vertical displacements, vertical stresses in the backfill, tensile forces in the strips and the influence of the stiffness of the soil are discussed. Copyright © 2017 John Wiley & Sons, Ltd.     

Intergranular diffusion rates from the analysis of garnet surfaces: Implications for metamorphic equilibration

Novel approaches to garnet analysis have been used to assess rates of intergranular diffusion between different matrix phases and garnet porphyroblasts in a regionally metamorphosed staurolite‐mica‐schist from the Barrovian‐type area in Scotland. X‐ray maps and chemical traverses of planar porphyroblast surfaces reveal chemical heterogeneity of the garnet grain boundary linked to the nature of the adjacent matrix phase. The garnet preserves evidence of low temperature retrograde exchange with matrix minerals and diffusion profiles documenting cation movement along the garnet boundaries. Garnet–quartz and garnet–plagioclase boundaries preserve evidence of sluggish Mg, Mn and Fe diffusion at comparable rates to volume diffusion in garnet, whereas diffusion along garnet–biotite interfaces is much more effective. Evidence of particularly slow Al transport, probably coupled to Fe³⁺ exchange, is locally preserved on garnet surfaces adjacent to Fe‐oxide phases. The Ca distribution on the garnet surface shows the most complex behaviour, with long‐wavelength heterogeneities apparently unrelated to the matrix grain boundaries. This implies that the Ca content of garnet is controlled by local availability and is thought likely to reflect disequilibrium established during garnet growth. Geochemical anomalies on the garnet surfaces are also linked to the location of triple junctions between the porphyroblasts and the matrix phases, and imply enhanced transport along these channels. The slow rates of intergranular diffusion and the characteristics of different boundary types may explain many features associated with the prograde growth of garnet porphyroblasts. Thus, minerals such as quartz, Fe‐oxides and plagioclase whose boundaries with garnet are characterized by slow intergranular diffusion rates appear to be preferentially trapped as inclusions within porphyroblasts. As such grain boundary diffusion rates may be a significant kinetic impediment to metamorphic equilibrium and garnet may struggle to maintain chemical and textural equilibrium during growth in pelites.     

Numerical and theoretical research on stress distribution in the loosening zone of the trapdoor problem

The traditional theory of soil arching effect was developed on the assumption that stress distribution in the loosening zone is uniform. However, because of the deflection of principal stress' direction, the stress distribution in the loosening zone is actually ununiform. For the evaluation of principal stress axis deflection and stress redistribution, a discrete element method numerical model of trapdoor problem is established for the simulation of soil arching effect. Based on the numerical results, an arc shape of major principal stress trajectory and uniform horizontal stress distribution at the same depth of the loosening zone are adopted. An analytical model is raised to estimate the average loosening earth pressure acting on the trapdoor and stress distribution in the loosening zone at a limit state. In addition, comparison studies are carried out between the predictions of the proposed solutions and discrete element method numerical results as well as available model test results, thereby validating the accuracy of the proposed theoretical model. Both numerical and theoretical results indicate that the vertical stress distribution in the loosening zone is obviously ununiform. The load acting in the middle of loosening zone is transferred toward two sides so that the vertical stress distribution in loosening zone is concave.     

Assessing the effectiveness of imperviousness on stormwater runoff in micro urban catchments by model simulation

Imperviousness, considered as a critical indicator of the hydrologic impacts of urbanization, has gained increasing attention both in the research field and in practice. However, the effectiveness of imperviousness on rainfall–runoff dynamics has not been fully determined in a fine spatiotemporal scale. In this study, 69 drainage subareas <1 ha of a typical residential catchment in Beijing were selected to evaluate the hydrologic impacts of imperviousness, under a typical storm event with a 3‐year return period. Two metrics, total impervious area (TIA) and effective impervious area (EIA), were identified to represent the impervious characteristics of the selected subareas. Three runoff variables, total runoff depth (TR), peak runoff depth (PR), and lag time (LT), were simulated by using a validated hydrologic model. Regression analyses were developed to explore the quantitative associations between imperviousness and runoff variables. Then, three scenarios were established to test the applicability of the results in considering the different infiltration conditions. Our results showed that runoff variables are significantly related to imperviousness. However, the hydrologic performances of TIA and EIA were scale dependent. Specifically, with finer spatial scale and the condition heavy rainfall, TIA rather than EIA was found to contribute more to TR and PR. EIA tended to have a greater impact on LT and showed a negative relationship. Moreover, the relative significance of TIA and EIA was maintained under the different infiltration conditions. These findings may provide potential implications for landscape and drainage design in urban areas, which help to mitigate the runoff risk. Copyright © 2015 John Wiley & Sons, Ltd.     

面向WebGL的矢量数据三维绘制技术

矢量数据的叠加显示能够提高三维虚拟地球的表达效果与分析能力。受限于GPU的计算精度,在三维虚拟地球中矢量数据绘制普遍存在抖动现象和深度冲突现象。对基于WebGL的矢量数据三维绘制中计算精度问题进行了分析,提出了使用CPU RTC技术和GPU RTE技术提高顶点变换的精度,使用多视锥渲染算法和深度平面技术解决深度缓存精度问题。实验证明,这几种技术和算法可以有效缓解抖动现象和深度冲突现象导致的视觉干扰,改善了各种尺度和范围的矢量数据在三维地形上的叠加显示效果。     

Selection and quality assessment of Landsat data for the North American forest dynamics forest history maps of the US

Using the NASA Earth Exchange platform, the North American Forest Dynamics (NAFD) project mapped forest history wall-to-wall, annually for the contiguous US (1986–2010) using the Vegetation Change Tracker algorithm. As with any effort to identify real changes in remotely sensed time-series, data gaps, shifts in seasonality, misregistration, inconsistent radiometry and cloud contamination can be sources of error. We discuss the NAFD image selection and processing stream (NISPS) that was designed to minimize these sources of error. The NISPS image quality assessments highlighted issues with the Landsat archive and metadata including inadequate georegistration, unreliability of the pre-2009 L5 cloud cover assessments algorithm, missing growing-season imagery and paucity of clear views. Assessment maps of Landsat 5–7 image quantities and qualities are presented that offer novel perspectives on the growing-season archive considered for this study. Over 150,000+ Landsat images were considered for the NAFD project. Optimally, one high quality cloud-free image in each year or a total of 12,152 images would be used. However, to accommodate data gaps and cloud/shadow contamination 23,338 images were needed. In 220 specific path-row image years no acceptable images were found resulting in data gaps in the annual national map products.     

Upper Maastrichtian shallow marine environments and neoselachian assemblages in North Iberian palaeomargin (Castilian Ramp,Spain)

Upper Maastrichtian deposits formed in a nearshore subtidal environment within the Valdenoceda Formation (Castilian Ramp, North Iberian margin) are described together with two recently found selachian assemblages. Rare earth element concentrations (REE) have been used to assess the degree of taphonomic mixing and reworking, observing that it is minor or non-existent, and differences in degree of preservation and ecologic mixing can be explained by biostratinomic processes. The patterns of REE also helped to obtain a better understanding of the depositional environment, including the diagenetic history from burial to final degree of bone preservation.The fossil assemblages here described are close to that of the late Maastrichtian of Albaina (in the enclave of Condado de Treviño, Burgos), both in the Basque-Cantabrian Region, but their age may be slightly older (early late Maastrichtian). In total, the new assemblages consist of 17 taxa, assigned to 11 genera of shallow-water dwellers combined with individuals from the outer shelf. They represent cosmopolitan taxa (Squalicorax pristodontus, Serratolamna serrata and Rhombodus binkhorsti) together with local species (Rhinobatos echavei, Rhinobatos ibericus). Although there are not significant differences between Albaina and Quintanilla la Ojada faunas, the new assemblages add interesting taphonomic and geochemical information to the few existing uppermost Cretaceous deposits with fossil sharks in southwestern Europe.