Numerical modeling of rainstorm-induced shallow landslides in saturated and unsaturated soils

For the assessment of shallow landslides triggered by rainfall, the physically based model coupling the infinite slope stability analysis with the hydrological modeling in nearly saturated soil has commonly been used due to its simplicity. However, in that model the rainfall infiltration in unsaturated soil could not be reliably simulated because a linear diffusion-type Richards’ equation rather than the complete Richards’ equation was used. In addition, the effect of matric suction on the shear strength of soil was not actually considered. Therefore, except the shallow landslide in saturated soil due to groundwater table rise, the shallow landslide induced by the loss in unsaturated shear strength due to the dissipation of matric suction could not be reliably assessed. In this study, a physically based model capable of assessing shallow landslides in variably saturated soils is developed by adopting the complete Richards’ equation with the effect of slope angle in the rainfall infiltration modeling and using the extended Mohr–Coulomb failure criterion to describe the unsaturated shear strength in the soil failure modeling. The influence of rainfall intensity and duration on shallow landslide is investigated using the developed model. The result shows that the rainfall intensity and duration seem to have similar influence on shallow landslides respectively triggered by the increase of positive pore water pressure in saturated soil and induced by the dissipation of matric suction in unsaturated soil. The rainfall duration threshold decreases with the increase in rainfall intensity, but remains constant for large rainfall intensity.     

Timing of UHP metamorphism in the Hong’an area,western Dabie Mountains,China: evidence from zircon U–Pb age,trace element and Hf isotope composition

The Hong’an area (western Dabie Mountains) is the westernmost terrane in the Qinling-Dabie-Sulu orogen that preserves UHP eclogites. The ages of the UHP metamorphism have not been well constrained, and thus hinder our understanding of the tectonic evolution of this area. LA-ICPMS U–Pb age, trace element and Hf isotope compositions of zircons of a granitic gneiss and an eclogite from the Xinxian UHP unit in the Hong’an area were analyzed to constrain the age of the UHP metamorphism. Most zircons are unzoned or show sector zoning. They have low trace element concentrations, without significant negative Eu anomalies. These metamorphic zircons can be further subdivided into two groups according to their U–Pb ages, and trace element and Lu–Hf isotope compositions. One group with an average age of 239 ± 2 Ma show relatively high and variable HREE contents (527 ≥ Lu_N ≥ 14) and ¹⁷⁶Lu/¹⁷⁷Hf ratios (0.00008–0.000931), indicating their growth prior to a great deal of garnet growth in the late stage of continental subduction. The other group yields an average age of 227 ± 2 Ma, and shows consistent low HREE contents and ¹⁷⁶Lu/¹⁷⁷Hf ratios, suggesting their growth with concurrent garnet crystallization and/or recrystallization. These two groups of age are taken as recording the time of prograde HP to UHP and retrograde UHP–HP stages, respectively. A few cores have high Th/U ratios, high trace element contents, and a clear negative Eu anomaly. These features support a magmatic origin of these zircon cores. The upper intercept ages of 771 ± 86 and 752 ± 70 Ma for the granitic gneiss and eclogite, respectively, indicate that their protoliths probably formed as a bimodal suite in rifting zones in the northern margin of the Yangtze Block. Young Hf model ages (T _DM1) of magmatic cores indicate juvenile (mantle-derived) materials were involved in their protolith formation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Metallographic cooling rates and origin of IVA iron meteorites

We have determined the metallographic cooling rates for 13 IVA irons using the most recent and most accurate metallographic cooling rate model. Group IVA irons have cooling rates that vary from 6600 °C/Myr at the low-Ni end of the group to 100 °C/Myr at the high-Ni end of the group. This large cooling rate range is totally incompatible with cooling in a mantled core which should have a uniform cooling rate. Thermal and fractional crystallization models have been used to describe the cooling and solidification of the IVA asteroid. The thermal model indicates that a metallic body of 150 ± 50 km in radius with less than 1 km of silicate on the outside of the body has a range of cooling rates that match the metallographic cooling rates in IVA irons in the temperature range 700-400 °C where the Widmanstätten pattern formed. The fractional crystallization model for Ni with initial S contents between 3 and 9 wt% is consistent with the measured variation of cooling rate with bulk Ni and the thermal model. New models for impacts in the early solar system and the evolution of the primordial asteroid belt allow us to propose that the IVA irons crystallized and cooled in a metallic body that was derived from a differentiated protoplanet during a grazing impact. Other large magmatic iron groups, IIAB, IIIAB, and IVB, also show significant cooling rate ranges and are very likely to share a similar history.     

Isotopic studies of leaf water. Part 1: A physically based two-dimensional model for pine needles

Understanding of isotopic variations in leaf water is important for reconstruction of paleoclimate and assessment of global biochemical processes. We report here a study of isotopic distributions within a single needle of two pine species, Pinus resinosa Ait and Pinus strobes L., with the objective of understanding how isotopic compositions of leaf water are controlled by environmental and physiological variables. A 2D model was developed to simulate along-leaf isotopic variations and bulk leaf water isotopic compositions. In addition to variables common to all leaf water isotopic models, this 2D model also takes into account the specific geometry and dimensions of pine needles and the isotopic transport in xylem and mesophyll. The model can successfully simulate oxygen isotopic variations along a single needle and averaged over a leaf (bulk leaf water). The simulations suggest that isotopic composition of the bulk leaf water does not always depend only upon the average transpiration rate, which in turn raises questions about using leaf water isotopic values to estimate transpiration rates. An unsuccessful attempt to simulate along-needle hydrogen isotopic variations suggests that certain unknown biological process(es) may not have been incorporated into our 2D model, and if so, it calls for a reevaluation of all other models for hydrogen isotopic simulations of leaf water since they too lack these processes.Existing leaf water isotopic models are reviewed in this work. In particular, we evaluate the most frequently used model, the stomatal boundary layer model (also referred to as the Craig-Gordon model). We point out that discrepancy between the boundary layer model and the measured bulk leaf water seems to depend upon relative humidity. Using our 2D model, we show that this humidity dependency is a result of an interplay between environmental and physiological conditions: if the transpiration rate of plant leaves decreases with increasing relative humidity, our 2D model can reproduce the pattern of isotopic discrepancy between boundary layer model predictions and observations, enabling us to understand better the reason behind this discrepancy.     

Impact of lightning on organic matter-rich soils: influence of soil grain size and organic matter content on underground fires

A detailed study was carried out on a piece of land that had been struck by lightning during the violent rainstorm that raged over the Island of S?o Miguel (Azores Archipelago) in late October 2006. Temperature and gas measurements (CO₂, CO, H₂S and CH₄) were performed in four study trenches, dug in an area of ∼3 m², where an underground fire had been initiated by the impact with a lightning stroke, followed by the emission of a column of gases and smoke. The soil under study was originally a well-pedogenized about 80 cm thick bed, made of volcanic clayey to silty tephra fallouts and contained 5.5–9.7% of organic matter. The underground fire was monitored for one week and revealed a peak release of 404 ppm CO and 3.4% CO₂ originating from a horizon located about 45 cm under the soil surface. Measurements of temperature, performed one week after the impact, indicated a maximum value of 326°C inside the soil, while 516.5°C were measured on the surface of a lava block interred about 20 cm under the surface. Subsequently, a stratigraphic and sedimentologic study proved the role of the grain-size of the soil and of the organic matter content of the different horizons of the impact area, in determining the ratio between anoxic/oxidised combustion conditions and in the progress of the process itself. It was also noticed that combustion was not total all over in the soil bed and that the process had slightly migrated toward SW during the observation period. The combustion process went on for about ten days, in spite of several other violent rainstorms, until it was artificially extinguished through the excavations made to obtain study trenches. This particular circumstance evidenced the potential natural hazard represented by this kind of atmospheric event, especially in a land where the volcanic nature of the soil may easily mislead inexperienced observers and, consequently, delay proper action.     

含硫化物脉状矿床成矿流体的中阶段 δD亏损：实例与原因

很多含硫化物脉状金银矿床的成矿过程可以大致分为3个阶段,从早到晚分别以形成黄铁矿-石英组合、多金属硫化物和碳酸盐网脉为标志.所谓中阶段δD亏损,是指在成矿流体δ18O逐渐降低(或升高)时,中阶段δD同时低于早阶段和晚阶段的δD.列举了一些矿床的氢氧同位素研究结果,例证了这种现象的存在;根据中阶段伴随大量硫化物沉淀的事实,认为硫化物沉淀时流体中H2S或HS-的H 被金属离子置换出来,H 加入到流体系统后导致了流体的中阶段δD亏损.对不同温度H2O-H2S体系1000In ot值的概算显示,在200-650℃温度范围,H2S或HS-的δD低于H2O的δD约211‰～478%‰,H2S或HS-中H 的少量释放即可导致流体水δD的显著降低.因此建议,成矿流体从早到晚δ18O逐渐降低,中阶段δD亏损的规律,可以作为含硫化物脉状矿床氢氧同位素演化的概念性模型.     

金川镍铜矿床数学模型对深部矿化变化趋势的指示

按矿化特征将金川X矿区主矿体分为致密块状特富矿、海绵状富矿和浸染状贫矿.应用MicroMine软件对3类矿化体分别建立了数学模型,并应用距离平方反比法对模型单元块进行了镍、铜品位估值.通过分析和对比矿床数学模型900m至1300m不同水平剖面上的Ni、Cu品位及Ni/Cu比值的变化,得出在X矿区的中部西侧是一个Ni和Cu的矿化中心,并且向深部Cu的矿化强度高于Ni.     

Background concentrations of elements in surface soils and their changes as affected by agriculture use in the desert-oasis ecotone in the middle of Heihe River Basin,North-west China

The concentrations of twenty four chemical elements in the surface layer of natural desert soils and the cultivated farmland soils were measured at a desert-oasis ecotone in the middle of Heihe river basin, north-west China. Background values were estimated for (a) major elements (Si 335.3 g kg^− 1, Al 49.4 g kg^− 1, Fe 19.1 g kg^− 1, Ca 29.4 g kg^− 1, Mg 8.9 g kg^− 1, K 20.1 g kg^− 1, Na 17.5 g kg^− 1 and P 0.338 g kg^− 1), (b) heavy metals and non-metals (Cr 55.8 mg kg^− 1, Mn 404.8 mg kg^− 1, Ni 17.7 mg kg^− 1, Cu 5.1 mg kg^− 1, Zn 33.7 mg kg^− 1, Pb 15.5 mg kg^− 1 and As 5.2 mg kg^− 1) and (c) other trace elements (Ti 2.0 mg kg^− 1, V 55.3 mg kg^− 1, Co 5.7 mg kg^− 1, Rb 82.4 mg kg^− 1, Sr 232.9 mg kg^− 1, Y 14.7 mg kg^− 1, Zr 194.9 mg kg^− 1, Nb 7.8 mg kg^− 1 and Ba 720.6 mg kg^− 1). After natural desert soil was cultivated for agricultural use, significant changes in element concentrations occurred under tillage, irrigation and fertilisation management. Compared to natural soil, the for the levels of Si, K, Na, Sr, Zr and Ba decreased, and no changes were observed for Rb, while the values of the other 17 elements increase in agricultural soil from 1.2 to 3.5 times. However, their absolute concentrations are still low, suggesting that the arable soil in this region remains comparatively a clean soil. The increased silt, clay and organic carbon content, under long-term irrigation, enriched the fine-grained materials, and application of fertilisers and manure contributed to the accumulation of most elements in arable soil. The accumulation of elements in agricultural soil increased with increasing cultivation years and extent of soil development.