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.     

川中—川南地区上三叠统沉积相研究

综合分析岩心、露头剖面和测井资料,认为川中—川南地区上三叠统发育三角洲相和湖泊相,其中,三角洲相可进一步划分为三角洲平原、三角洲前缘和前三角洲三个亚相,湖泊相则可细分为滨湖和浅湖两个亚相。不同时期三角洲相和湖泊相发育程度不同。须一段、须三段、须五段和须六段以湖泊沉积为主,由东向西,水体逐渐加深。须二段和须四段以三角洲沉积为主,湖泊沉积分布范围较小。须二段发育时期,研究区北部构造活动强烈,沉积物供给充分,三角洲相主要由北向南延伸;须四段沉积时期,由于北部构造活动逐渐减弱,东南部沉积物供给增加,三角洲相主要由东南向西北方向延伸。随着周缘构造活动的变化,研究区沉积范围逐渐扩大。早期沉积范围仅局限于宜宾—泸州以北地区,中晚期沉积范围扩大,宜宾—泸州以南地区开始接收沉积。