Fluid infiltration during contact metamorphism of interbedded marble and calc-silicate hornfels, Twin Lakes area, central Sierra Nevada, California

Abstract In the Twin Lakes area, central Sierra Nevada, California, most contact metamorphosed marbles contain calcite + dolomite + forsterite ± diopside ± phlogopite ± tremolite, and most calc-silicate hornfelses contain calcite + diopside + wollastonite + quartz ± anorthite ± K-feldspar ± grossular ± titanite. Mineral-fluid equilibria involving calcite + dolomite + tremolite + diopside + forsterite in two marble samples and wollastonite + anorthite + quartz + grossular in three hornfels samples record P± 3 kbar and T± 630° C. Various isobaric univariant assemblages record CO₂-H₂O fluid compositions of χCO₂= 0.61–0.74 in the marbles and χCO₂= 0.11 in the hornfelses. Assuming a siliceous dolomitic limestone protolith consisting of dolomite + quartz ° Calcite ± K-feldspar ± muscovite ± rutile, all plausible prograde reaction pathways were deduced for marble and hornfels on isobaric T-XCO₂ diagrams in the model system K₂O-CaO-MgO-Al₂O₃-SiO₂-H₂O-CO₂. Progress of the prograde reactions was estimated from measured modes and mass-balance calculations. Time-integrated fluxes of reactive fluid which infiltrated samples were computed for a temperature gradient of 150 °C/km along the fluid flow path, calculated fluid compositions, and estimated reaction progress using the mass-continuity equation. Marbles and hornfelses record values in the range 0.1–3.6 × 10⁴ cm³/cm² and 4.8–12.9 × 10⁴ cm³/cm², respectively. For an estimated duration of metamorphism of 10⁵ years, average in situ metamorphic rock permeabilities, calculated from Darcy's Law, are 0.1–8 × 10^?6 D in the marbles and 10–27 × 10^?6 D in the hornfelses. Reactive metamorphic fluids flowed up-temperature, and were preferentially channellized in hornfelses relative to the marbles. These results appear to give a general characterization of hydrothermal activity during contact metamorphism of small pendants and screens (dimensions ± 1 km or less) associated with emplacement of the Sierra Nevada batholith.     

Screening of polymers on selected Hawaii soils for erosion reduction and particle settling

In recent years, high‐molecular‐weight anionic polyacrylamides (PAMs) have been tested on a variety of soils, primarily in temperate climates. However, little information is available regarding the effectiveness of PAM for preventing soil loss through runoff in tropical settings. Screening tests were performed using three negatively charged PAMs and one positively charged PAM on five Hawaii soils (two Oxisols, one Vertisol, and two Aridisols) to determine erosion loss, sediment settling, and aggregate stability. A laboratory‐scale rainfall simulator was used to apply erosive rainfall at intensities from 5 to 8·5 cm h^?1 at various PAM doses applied in both dry and solution forms. Soil detachment due to splash and runoff, as well as the runoff and percolate water volumes, were measured for initial and successive storms. The impact of PAM on particle settling and aggregate stability was also evaluated for selected soil‐treatment combinations. Among the PAMs, Superfloc A‐836 was most effective, and significantly reduced runoff and splash sediment loss for the Wahiawa Oxisol and Pakini Andisol at rates varying between 10 and 50 kg ha^?1. Reduced runoff and splash sediment loss were also noted for PAM Aerotil‐D when applied in solution form to the Wahiawa Oxisol. Significant reductions in soil loss were not noted for either the Lualualei Vertisol or the Holomua Oxisol. It is believed that the high montmorillonite content of the Lualualei Vertisol and the low cation‐exchange capacity of the Holomua Oxisol diminished the effectiveness of the various PAMs tested. The polymers were also found to enhance sediment settling of all soils and helped improve their aggregate stability. This screening study shows the potential use of PAM for tropical soils for applications such as infiltration enhancement, runoff reduction, and enhanced sedimentation of detention ponds. Copyright © 2005 John Wiley & Sons, Ltd.     

Analytical solutions of one‐dimensional infiltration before and after ponding

This paper is a continuation of previous research, which obtained a convenient solution for arbitrary surface fluxes before ponding. By means of Fourier Transformation this has been extended to derive analytical solutions of a linearized Richards' equation for arbitrary input fluxes after surface saturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Assessing soil erosion in Mediterranean karst landscapes of Lebanon using remote sensing and GIS

Maps showing the potential for soil erosion at 1:100,000 scale are produced in a study area within Lebanon that can be used for evaluating erosion of Mediterranean karstic terrain with two different sets of impact factors built into an erosion model. The first set of factors is: soil erodibility, morphology, land cover/use and rainfall erosivity. The second is obtained by the first adding a fifth factor, rock infiltration. High infiltration can reflect high recharge, therefore decreasing the potential of surface runoff and hence the quantity of transported materials. Infiltration is derived as a function of lithology, lineament density, karstification and drainage density, all of which can be easily extracted from satellite imagery. The influence of these factors is assessed by a weight/rate approach sharing similarities between quantitative and qualitative methods and depending on pair-wise comparison matrix.The main outcome was the production of factorial maps and erosion susceptibility maps (scale 1:100,000). Spatial and attribute comparison of erosion maps indicates that the model that includes a measure of rock infiltration better represents erosion potential. Field investigation of rills and gullies shows 87.5% precision of the model with rock infiltration. This is 17.5% greater than the precision of the model without rock infiltration. These results indicate the necessity and importance of integrating information on infiltration of rock outcrops to assess soil erosion in Mediterranean karst landscapes.     

基于霍顿下渗能力曲线的流量过程线连续分割方法研究

在用前期影响雨量代替流域蓄水量的基础上,提出了基于霍顿下渗能力曲线的流量过程线的连续分割方法.选择西峡、猴子岩和东湾三个流域的降雨径流资料,采用该方法对流量过程线进行分割,并与现行的数字滤波法、非线性水库法和Boussinesq方程法相比较.结果表明,该方法具有一定的物理基础、符合产汇流基本规律,而且能够减少涨洪段流量过程线分割的主观性,对于流域时段单位线和降雨径流关系的推求均有重要的意义.     

Constraints on the application of carbon isotope thermometry in high- to ultrahigh-temperature metamorphic terranes

Nine marble horizons from the granulite facies terrane of southern India were examined in detail for stable carbon and oxygen isotopes in calcite and carbon isotopes in graphite. The marbles in Trivandrum Block show coupled lowering of δ¹³C and δ¹⁸O values in calcite and heterogeneous single crystal δ¹³C values (? 1 to ? 10‰) for graphite indicating varying carbon isotope fractionation between calcite and graphite, despite the granulite facies regional metamorphic conditions. The stable isotope patterns suggest alteration of δ¹³C and δ¹⁸O values in marbles by infiltration of low δ¹³C–δ¹⁸O‐bearing fluids, the extent of alteration being a direct function of the fluid‐rock ratio. The carbon isotope zonation preserved in graphite suggests that the graphite crystals precipitated/recrystallized in the presence of an externally derived CO₂‐rich fluid, and that the infiltration had occurred under high temperature and low f_O2 conditions during metamorphism. The onset of graphite precipitation resulted in a depletion of the carbon isotope values of the remaining fluid+calcite carbon reservoir, following a Rayleigh‐type distillation process within fluid‐rich pockets/pathways in marbles resulting in the observed zonation. The results suggest that calcite–graphite thermometry cannot be applied in marbles that are affected by external carbonic fluid infiltration. However, marble horizons in the Madurai Block, where the effect of fluid infiltration is not detected, record clear imprints of ultrahigh temperature metamorphism (800–1000 °C), with fractionations reaching <2‰. Zonation studies on graphite show a nominal rimward lowering δ¹³C on the order of 1 to 2‰. The zonation carries the imprint of fluid deficient/absent UHT metamorphism. Commonly, calculated core temperatures are > 1000 °C and would be consistent with UHT metamorphism.     

The precipitation infiltration and runoff recharging experiment observation in the Taklimakan Desert

１　ＩＮＴＲＯＤＵＣＴＩＯＮＴｈｅｌｏｃａｌｐｒｅｃｉｐｉｔａｔｉｏｎｉｓｍｏｓｔｉｍｐｏｒｔａｎｔａｖａｉｌａｂｌｅｆｒｅｓｈｗａｔｅｒｓｏｕｒｃｅｉｎｄｅｓｅｒｔｓ．ＴｈｅＴａｋｌｉｍａｋａｎＤｅｓｅｒｔｂｅｌｏｎｇｓｔｏａｎｅｘｔｒｅｍｅａｒｉｄｂｅｌｔ．Ａｃｃｏｒｄｉｎｇｔｏａｆｅｗｙｅａｒｓ’ｌｉｍｉｔｅｄｏｂｓｅｒｖａｔｉｏｎｄａｔａ,ｉｔｉｓｅｓｔｉｍａｔｅｄｔｈａｔｔｈｅｐｒｅｃｉｐｉｔａｔｉｏｎｓｈｏｕｌｄｂｅｉｎｔｈｅｒａｎｇｅｏｆ３０－５０ｍｍａｎｎｕａｌｌｙ．Ｈｏｗｅｖ…     

降雨条件下考虑饱和渗透系数变异性的边坡可靠度分析

土体饱和渗透系数表现为天然的变异性,为此基于Green-Ampt模型建立了考虑饱和渗透系数变异性的降雨入渗物理模型,并藉此模型确定了坡体湿润锋深度和含水率分布。然后结合无限长非饱和土边坡稳定模型得到解析形式的反映边坡稳定性的极限状态函数。采用Monte Carlo法对饱和渗透系数进行随机抽样并最终建立降雨条件下考虑饱和渗透系数变异性的边坡概率分析框架。针对一假想边坡,探讨了饱和渗透系数的变异系数、降雨持时和降雨强度对边坡破坏概率以及破坏发生时间概率分布的影响,结果表明：在降雨初期,边坡的破坏概率随饱和渗透系数变异性的增强而逐渐增加,但随着降雨的持续,破坏概率开始随变异性的增强而显著降低;滑坡最可能发生时间的大小并不受饱和渗透系数变异性的影响,而是直接取决于降雨强度;滑坡最可能发生时间所对应的概率却随变异性的增强而逐渐减小。     

考虑径流影响的滑坡降雨入渗二维有限元模拟及应用

降雨是滑坡失稳的常见诱发因素之一,现有数值模拟方法由于不能考虑来自基岩边坡的径流对滑体渗流的补给而低估了真实的降雨入渗水量,导致在评价降雨对滑坡稳定性特别是深层滑动影响时存在一定不足。以Richards方程和有限元法为基础,忽略降雨对渗透性极低的滑床和滑带的影响,将滑坡渗流计算域缩小为滑体,以避免因滑体与滑床（滑带）渗透性差异巨大引起的数值计算困难;依据基岩边坡水平长度和滑体降雨入渗边界饱和情况,修正降雨入渗边界,实现了考虑径流补给的滑坡降雨入渗简化数值模拟。算例表明,该方法所得渗流场更加符合实际情况。以三峡库区某滑坡为例,模拟了2006年10月～2009年12月间滑坡在库水升降和降雨条件下的渗流场和安全系数的演化过程;计算结果表明,考虑径流补给时滑坡后部的渗流场饱和区域明显较大,稳定系数降低较多,与位移监测资料显示的佐证较为吻合;若不考虑径流补给,则降雨对滑坡稳定性影响不大。