Landform‐oriented flow‐routing algorithm for the dual‐structure loess terrain based on digital elevation models

The loess landform in the Loess Plateau of China is with typical dual structure, namely, the upper smooth positive terrain and the lower cliffy negative terrain (P–N terrain for short). Obvious differences in their morphological feature, geomorphological mechanism, and hydrological process could be found in the both areas. Based on the differences, a flow‐routing algorithm that separately addresses the dual‐structure terrain would be necessary to encompass this spatial variation in their hydrological behaviour. This paper proposes a mixed flow‐routing algorithm to address aforementioned problems. First, the loess landform surface is divided into P–N terrains based on digital elevation models. Then, specific catchment area is calculated with the new algorithm to simulate the water flows in both positive and negative terrain areas. The mixed algorithm consists of the multiple flow‐routing algorithm (multiple‐flow direction) for positive areas and the D8 algorithm for negative areas, respectively. The approach is validated in two typical geomorphologic areas with low hills and dense gullies in the northern Shaanxi Loess Plateau. Four indices are used to examine the results, which show that the new algorithm is more suitable for loess terrain in simulating the spatial distribution of water accumulation, as well as in modeling the flow characteristics of the true surface by considering the morphological structures of the terrain. Copyright © 2013 John Wiley & Sons, Ltd.     

Synchronous pattern of moisture availability on the southern Tibetan Plateau since 17.5 cal. ka BP – the Tangra Yumco lake sediment record

A possible asynchronicity of the spatial and temporal moisture availability on the Tibetan Plateau has been a controversial subject of discussion in recent years. Here we present the first attempt to systematically investigate possible spatial and temporal variations in moisture availability by examining two lakes, Tangra Yumco and Nam Co, on an east–west transect on the southern Tibetan Plateau using identical proxies for palaeoenvironmental reconstruction. In this study, an independent record from Tangra Yumco was analysed applying a multi‐proxy approach to reconstruct variations in moisture availability since the Lateglacial. Results were subsequently compared with previously published records from Nam Co and additional records from Tso Moriri (northwestern Himalaya) and Naleng Co (southeastern Tibetan Plateau). Our results show that Tangra Yumco was at least partially ice covered prior to 17.1 cal. ka BP. A temperature rise after 17.1 cal. ka BP probably resulted in thawing of the permafrost. At 16.0 cal. ka BP moisture availability increased, representing an initial monsoonal intensification. Warmer conditions between 13.0 and 12.4 cal. ka BP and cooler conditions between 12.4 cal. ka BP and the onset of the Holocene reflect the Bølling‐Allerød and Younger Dryas. At the onset of the Holocene moisture availability rapidly increased, with moisture highest prior to 8.5 cal. ka BP when temperatures were also highest. After 8.5 cal. ka BP the moisture availability gradually decreased and showed only minor amplitude variations. These findings are consistent with the records from large lakes like Nam Co, Tso Moriri, and Naleng Co, revealing a synchronous pattern of moisture availability on the southern Tibetan Plateau.     

Application of the numerical manifold method for stress wave propagation across rock masses

In this paper, the numerical manifold method (NMM) is extended to study wave propagation across rock masses. First, improvements to the system equations, contact treatment, and boundary conditions of the NMM are performed, where new system equations are derived based on the Newmark assumption of the space–time relationship, the edge‐to‐edge contact treatment is further developed for the NMM to handle stress wave propagation across discontinuities, and the viscous non‐reflection boundary condition is derived based on the energy minimisation principle. After the modification, numerical comparisons between the original and improved NMM are presented. The results show that the original system equations result in artificial numerical damping, which can be overcome by the Newmark system equations. Meanwhile, the original contact scheme suffers some calculation problems when modelling stress wave propagation across a discontinuity, which can be solved by the proposed edge‐to‐edge contact scheme. Subsequently, the influence of the mesh size and time step on the improved NMM for stress wave propagation is studied. Finally, 2D wave propagation is modelled, and the model's results are in good agreement with the analytical solution. Copyright © 2013 John Wiley & Sons, Ltd.     

Sensitivity of advective travel time of contaminants to correlated formation porosities

An approach is presented to quantify sensitivity of advective contaminant travel time to porosities of hydrogeologic units (HGUs) along the flowpaths when the porosities of different HGUs are correlated. The approach is an extension of the previous sensitivity analysis technique for independent input porosity cases. It is applicable in situations where a calibrated groundwater flow model exists, but a full contaminant transport model is not available. Three sensitivity indices are introduced based on the decomposition of covariance between the advective contaminant travel time and individual input porosities of HGUs. When the input HGU porosities are correlated, the three sensitivity indices quantify the total, intrinsic and correlated contributions from each individual HGU porosity, which should be considered in order to determine the relative importance to the uncertainty in advective contaminant travel time of the input HGU porosities contributing either independently or correlatedly. Two simple one-dimensional flow examples are presented to illustrate the applications of the approach to scenarios when each HGU has distinct porosity and situations of spatially variable porosity field. The approach is applicable to more complex multi-dimensional cases where advective contaminant travel time can be calculated based on simulated flow results from groundwater flow models.     

Iron isotopic systematics of UHP eclogites respond to oxidizing fluid during exhumation

The iron stable isotope compositions (δ⁵⁶Fe) and iron valence states of ultrahigh‐pressure eclogites from Bixiling in the Dabie orogen belt, China, were measured to trace the changes of geochemical conditions during vertical transportation of earth materials, for example, oxygen fugacity. The bulk Fe³⁺/ΣFe ratios of retrograde eclogites, determined by Mössbauer spectroscopy, are consistently higher than those of fresh eclogites, suggesting oxidation during retrograde metamorphism and fluid infiltration. The studied eclogites (five samples) display limited mid‐ocean ridge basalts (MORB)‐like (~0.10‰) δ⁵⁶Fe values, which are indistinguishable from their protoliths, that is, gabbro cumulates formed through differentiation of mantle‐derived basaltic magma. This suggests that Fe isotope fractionation during continental subduction is limited. Garnet separates display limited δ⁵⁶Fe variation ranging from ?0.08 ± 0.07 ‰ to 0.02 ± 0.07‰, whereas coexisting omphacite displays a large variation of δ⁵⁶Fe values from 0.15 ± 0.07‰ to 0.47 ± 0.07‰. Omphacite also has highly variable Fe³⁺/ΣFe ratios from 0.367 ± 0.025 to 0.598 ± 0.024, indicating modification after peak metamorphism. Omphacite from retrograde eclogites has elevated Fe³⁺/ΣFe ratios (0.54–0.60) compared to that from fresh eclogites (~0.37), whereas garnet displays a narrow range of ferric iron content with Fe³⁺/ΣFe ratios from 0.039 ± 0.013 to 0.065 ± 0.022. The homogenous δ⁵⁶Fe values and Fe³⁺/ΣFe ratios of garnet suggest that it survived the retrograde metamorphism and preserved its Fe‐isotopic features and ferric contents of peak metamorphism. Because of similar diffusion rates of Fe and Mg in garnet and omphacite, and constant Δ²⁶Mg_{omphacite‐garnet} values (1.14 ± 0.04‰), equilibrium iron isotope fractionation between garnet and omphacite was probably achieved during peak metamorphism. Elevated Fe³⁺/ΣFe ratios of omphacite from retrograde eclogites and variant Δ⁵⁶Fe_{omphacite‐garnet} values of the studied eclogites (0.13 ± 0.10‰ to 0.48 ± 0.10‰) indicate that oxidized geofluid infiltration resulted in the elevation of δ⁵⁶Fe values of omphacite during retrograde metamorphism.     

依兰-伊通断裂新构造活动规律分析

作为郯庐断裂带北段主干的依兰-伊通断裂, 其新构造活动性与活动规律仍然存在不同的认识.本次工作通过详细的野外调查, 发现该断裂内活断层广泛存在, 由东、西两支北东走向的主干活断层构成, 沿着古近纪地堑边界断层发育.这些活断层主要呈破碎型结构, 多为逆右行平移活动.通过对这些活断层一系列实测擦痕反演应力场, 显示它们多是在东西向挤压中活动的, 而现今应力场转变为北东东-南西西向区域性挤压.依据本次野外观察与¹⁴ C定年, 并结合前人定年结果与近代地震分布, 表明依兰-伊通西支活断层的最新活动时代为全新世与晚更新世相间, 而东支活断层的最新活动时代主要为早-中更新世.依兰-伊通断裂内活断层显示了明显的差异性活动, 表现为西支的活动强度明显大于东支, 西支的最新活动时代皆晚于东支, 沿走向上活动性强、弱相间与最新活动时代不断变化, 以及近代地震活动不均一分布.它们沿走向上的分段性、差异性活动主要是因为被一系列北西向断层切断所致.