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.     

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.     

碳酸盐岩裂缝与岩溶作用研究

从构造与非构造作用角度出发,突出研究断层对裂缝的控制作用、断裂与岩溶的关系和裂缝与岩溶之间的关系,揭示碳酸盐岩裂缝与岩溶的发育规律。研究表明：构造裂缝为岩溶发育提供岩溶作用进行的有利通道,非构造裂缝则可以加速溶蚀作用的进行,促进岩溶的发育;裂缝与岩溶之间的关系始终是相辅相成、相互促进的关系。     

Wenchuan Earthquake： Way of thinking is changed

On 12 May, 2008, at 14：28：04 local time, an Ms 8.0 earthquake struck the Wenchuan County, Sichuan Province, in southwest China （Figure 1）. This quake left about 70,000 people dead, 18,000 missing and over 370,000 injured. In the 75 days after the main shock, 20,000 aftershocks were recorded, which included 241 with Ms≥4.0, 205 with 4.0≤M≤4.9, 30 with 5.0≤Ms≤5.9, and 6 with Ms≥6.0. The largest aftershock took place 13 days after the main shock, in Qingchuan County, 110 km northeast of the epicenter of the quake （Chen, 2008）. The losses from the quaketriggered geological disasters accounted for over a third of the total quake losses, which is extremely rare in the history of earthquake disasters. Secondary disasters, such as landslides, rockslides, landfalls, debris flows, etc., were widely distributed. Large numbers of buildings collapsed, including ones with steel reinforcing in the area of the epicenter. Infrastructure was badly damaged. In many areas, communications, and power and water supplies were cut off. The earthquake left over 35 large quake-dammed lakes, which threatened further disasters in the catchment areas.     

河南洛宁段河金矿流体包裹体研究和矿床成因

河南省洛宁县段河石英脉型金矿主要包括石寨沟和岭东两个矿区,分别由3～4条含金石英脉构成。矿化过程从早到晚包括石英-黄铁矿、石英-多金属硫化物和石英-碳酸盐等3个阶段．其中中阶段金矿化最强,次为早阶段。各阶段石英中流体包裹体以气液两相包裹体为主．次为纯液体包裹体。激光拉曼测试表明,气液两相包裹体的液相为H2O,气相主要为Ho和CO2混合、纯H2O,次为纯CO2;纯液体包裹体为纯H2O。石寨沟矿区包裹体均一温度从早到晚依次为240．9～315．9℃．188．7～304．5℃,137．3～259．3℃：流体盐度变化依次为（6．74～12．85）wt%NaCl．eq,（2．41～8．68）wt％NaCl．eq,（2．24-7．86）wt%NaCl．eq。岭东矿区均一温度从早到晚依次为303．7-343．1℃,251．8-325．4℃,305．7～355．0℃：流体盐度变化依次为（5．11～11．70）wt％NaCl．eq,（2．74-10．11）wt％NaCl．eq,（0．53-6．74）wt％NaCl．eq。两矿区主成矿期流体均为中温、低盐度,早阶段流体为改造热液和变质热液的混合体,含一定量CO2,且流体CO2含量和盐度从早到晚逐渐降低。石寨沟矿区包裹体均一温度逐渐降低,而岭东矿区包裹体均一温度先降后升,加之岭东矿区各阶段成矿温度均高于石寨沟矿区．表明成矿流体系统主要受岩浆热驱动,岭东矿区更靠近岩体,且在晚阶段又有脉动性的岩浆加热．段河金矿区南部存在隐伏岩体。