青藏高原东北部观赏地被植物马蔺根系固土护坡效应研究

为了进一步丰富青藏高原东北部地区植被护坡工程植物选择,推动区内植被护坡工程实用性和观赏性兼顾发展。该项研究以青藏高原东北部观赏地被植物马蔺为研究对象,以青海省大通县朔北藏族乡至东峡镇之间区域为研究区,在野外植物分布特征和生长特征调查的基础上,结合单根拉伸试验和重塑根-土复合体三轴压缩试验,定量评价了野生马蔺植株和根系生长特征、单根抗拉特性以及根系对土体抗剪强度的增强效应。试验结果表明：马蔺主要分布在研究区内洪积扇、坡积裙和河漫滩等地貌单元,立地条件包括上述地貌单元内的林下坡地、田间、道旁、湿地、荒地等;通过选取两处试验区开展调查可知,马蔺平均株高、冠幅和植株密度分别为75.22±11.40 cm、106.09±25.62 cm和1.51±0.55 株/m2;马蔺标准株根系分布深度可达50 cm,主要分布在0~30 cm深度范围内,根幅约40~50 cm;马蔺根径为0.20～0.70 mm,单根抗拉力、抗拉强度和延伸率平均值分别为7.94±2.91 N,46.30± 11.06 MPa和54.17%±17.08%,随着根径的增大,单根抗拉力和单根抗拉强度分别呈幂函数增大和幂函数降低关系,单根极限延伸率随着根径的增大呈逐渐增大变化趋势,但二者之间未表现出相对显著性的关系;马蔺根系对土体黏聚力和内摩擦角均有增强作用,可分别显著增强0~30 cm和0~10 cm深度范围内土体黏聚力和内摩擦角,且对土体黏聚力的增强作用较为显著。该项研究成果对于选用观赏地被植物马蔺进行研究区坡面水土流失、浅层滑坡等地质灾害现象的生态防护工作,具有理论研究价值和实际指导意义。     

干旱区梨园不同覆盖条件下土壤环境因子综合性评价研究

在干旱沿黄灌区开展不同覆盖材料配合滴灌的灌溉保墒方式下,进一步研究了不同覆盖方式对梨园全生育期的土壤理化性状和养分含量变化影响,分析了土壤温度、水分、pH等和矿质营养元素间的相互关系,对覆盖后的土壤质量进行了综合性评价。试验处理分为无覆盖对照（T1）、园艺地布覆盖（T2）、玉米秸秆覆盖（T3）和黑地膜覆盖（T4）4个处理;试验设计为随机区组设计,每个处理小区均为167株梨树（约占地667 m²）,重复3次;各小区土样分0~20 cm、20~40 cm、40~60 cm土层取样。结果表明：（1） 园艺地布和黑地膜覆盖处理有一定的增温效应,而玉米秸秆覆盖有较好的降温和稳温效应,且增墒效应明显,土壤含水率比其他3种处理显著提高了1.0%~2.7%。（2） 相较无覆盖处理,黑地膜覆盖可提升表层土壤pH,并加速表层土壤有机质的分解,有机质含量较无覆盖处理下降33.1%;玉米秸秆覆盖可显著降低各土层的土壤pH,降低范围为1.8%~4.6%,并促进0~20 cm土层有机质的提升,土壤有机质含量增加12.2%;园艺地布覆盖下0~40 cm土层内土壤有机质含量和全盐量均有降低。（3） 黑地膜覆盖下0~20 cm和20~40 cm土层的碱解氮含量分别为73.00 mg·g^-1和64.53 mg·g^-1,均显著地高于无覆盖处理,无覆盖条件下土壤碱解氮在深层（40~60 cm土层）积累较多,显著地高于玉米秸秆和黑地膜覆盖;各处理0~20 cm和20~40 cm土层速效磷含量差异显著,大小顺序均为T4>T3>T2>T1,玉米秸秆覆盖可提升土壤速效钾和速效铁的含量。（4） 进行主成分分析表明不同覆盖方式对梨园浅层土壤环境因子的影响要明显大于深层土壤,在0~40 cm土层内各覆盖处理效果均好于无覆盖处理,其中玉米秸秆覆盖在0~20 cm和20~40 cm土层综合得分分别为1.189和0.326,覆盖效果最佳。     

1981—2018年内蒙古典型草原季节性冻土对气候变化的响应

以1981—2018年内蒙古典型草原季节性冻土为研究对象,通过气候倾向率、Mann-Kendall法、多元线性回归等方法,分析最大冻土时空分布特征、年际、年代际变化,研究影响最大冻土深度变化的气象因子。结果表明：（1） 内蒙古典型草原季节性冻土冻结初日在9—11月,终日在4—6月,年内最大冻土深度出现在2—3月,深度在100~280 cm之间。（2） 最大冻土深度年际变化分为下开口抛物线型、上开口抛物线型、正弦曲线型,从最大冻土深度气候倾向率看呈现减小趋势的站点有68%。（3） 最大冻土深度年代际变化分为逐年代递减、减-增型和无明显变化规律,50%的站点在1989年以后最大冻土深度发生突变。（4） 多元线性回归表明气温冻结指数、年平均风速、年极端最低气温对最大冻土深度产生显著影响。该研究揭示了最大冻土深度存在退化的事实,为草原应对气候变化提供指导,为陆地土壤和大气碳循环交换的研究给出提示。     

基于包络检波和STFT谱分析的探地雷达土壤分层信息识别

土壤分层信息,特别是表土层结构,对土地生产力具有重要影响,是评价土壤质量的一个重要指标。为了快速、准确地获取土壤分层信息,本文利用探地雷达对分层土壤进行了回波信号采集,并分别在时域和频域分析土壤层位置和层厚信息。首先在信号预处理的基础上,借助包络检波方法确定在土壤分层界面在时域上的位置;然后获取电磁波速度,得到土壤分层厚度。考虑到土壤介电常数与电磁波在土壤中传播速度的相关性,采用短时傅里叶变换方法（Short-time Fourier Transform,STFT）获取各土壤层时频域特征值,并利用回归分析建立特征值与介电常数之间的数学关系,实现对各土壤层的介电常数估算,从而计算出电磁波传播速度,进而确定土壤各层厚度。为验证算法的有效性,分别对理想模拟实验环境和农田环境进行了探地雷达实验,结果表明利用包络检波对探地雷达回波信息进行分析,土壤层检出率达到94.5%,借助STFT谱分析进行探地雷达回波速度估计,对于70 cm深度以上土层厚度计算误差大都保持在10%以下,但随着土壤深度的增加,误差变大。总体来说,本方法能有效识别浅层土壤的分层信息,可应用于实际生产中耕层厚度的估测。     

1980—2020年中国九大流域蒸散发及其组分时空评估

蒸散发是地表陆气水分交换的纽带,准确量化蒸散发的时空演变格局对于水资源规划与管理至关重要。本文基于GLEAM模型的蒸散发及其组分数据集,借助7个通量观测站数据、120个流域的流域水量平衡及PML_V2蒸散发产品,在中国九大流域系统评估了GLEAM-ET产品,分析了植被恢复背景下,蒸散发（ET）及其组分（植被蒸腾Ec,截留蒸发Ei,土壤蒸发Es）在1980—2020年的时空演变格局。本文主要得到以下结论：① GLEAM-ET产品在中国九大流域具有较好的适用性,其性能与气候类型有关,干旱区效果优于湿润区。此外,GLEAM与PML_V2模型在九大流域相关性较好（R>0.7）,分布格局与变化趋势整体保持一致。② 全国尺度上,ET均值为416.88 mm,增长速率为1.21 mm/a。Ec与ET均呈自东南向西北递减的分布格局,而Es与其相反。Ec、ET在九大流域均呈显著增加趋势（p<0.001）。Ei、Es在季风区流域分别呈显著增加和显著减小趋势;在内陆区流域呈不显著减小（p>0.05）和显著增加趋势。③在植被恢复背景下,ET组分比例发生了变化。Ec占比变化存在南北差异,南方流域Ec占比均减小,北方流域均增加。Ei占比在各流域均增加,Es占比均减小。黄河流域ET组分对植被恢复的响应最为明显,Ec占比增加了5.21%,Es占比减小了5.56%。     

Factors and processes of permanent gully evolution in a Mediterranean marly environment (Cape Bon,Tunisia)

Abstract

Gully erosion is considered to be one of the most important soil erosion processes in Mediterranean marly environments, but its actual contribution to total soil loss is still under discussion. The objectives of this paper are: (a) to acquire the distributed value of erosion rate in a permanent gully developed on a marly substratum in a Mediterranean environment; and (b) to quantify the key factors responsible for the spatial and temporal differences in erosion rates observed within the gully. A permanent gully located in Cap Bon (northeastern Tunisia) has been intensively and regularly monitored over a 7-year period with electronic survey equipment (total station) to give five field topographic surveys, as well as hydrological measurements at the gully outlet. The net soil loss for the 7-year period comprised a denudation of 51 m³ of sediment on the gully bank slopes, which corresponds to a mean soil loss of 61 m³ ha^?1 year^?1 or 6.1 mm year^?1. Denudation was observed on bed units with a slope gradient greater than 20%, while the remainder showed deposition. By confirming the factors involved in gully evolution, and by refining the statistical link between factors and erosion rates within the gully, the results provide important information to predict gully erosion rates in Mediterranean marly environments.

Editor Z.W. Kundzewicz; Associate editor G. Mahé

Citation El Khalili, A., Raclot, D., Habaieb, H., and Lamachère, J.M., 2013. Factors and processes of permanent gully evolution in a Mediterranean marly environment (Cape Bon, Tunisia). Hydrological Sciences Journal, 58 (7), 1519–1531.     

An iterative pressure‐stabilized fractional step algorithm in saturated soil dynamics

This paper presents an iterative, incremental pressure‐stabilized fractional step algorithm for coupled hydro‐mechanical problems with mixed formulations of the displacement–pressure ( u –p) model in saturated soil dynamics that allows the use of finite elements with equal low order of interpolation approximation of u and p. In comparison with the original fractional step algorithm, the distinct features of the proposed algorithm lie in its enhanced stability owing to the introduction of both an iteration procedure and a finite increment calculus (FIC) process into the algorithm. The introduction of the iterative procedure makes the velocity term satisfy the momentum conservation equation in an implicit sense and allows much larger time step sizes to be used than those limited in existing explicit and semi‐implicit versions of the algorithm. The introduction of the FIC process removes the dependence of the stability of the proposed algorithm on the time step size, as a result it allows to using the incremental version of the algorithm and evades the minimum time step size requirement presented in the existing versions of the fractional step algorithm that restricts the application of the algorithm to saturated soil dynamics problems with high frequencies. Numerical experiments demonstrate the effectiveness and improved performance of the proposed iterative pressure‐stabilized fractional step algorithm. Copyright © 2009 John Wiley & Sons, Ltd.     

Time‐scale analysis in unsaturated porous media under external wave loads

Extreme waves caused by tsunamis and storm surges can lead to soil failures in the near‐shore region, which may have severe impact on coastal environments and communities. Multiphase flows in deformable porous media involve several coupled processes and multiple time scales, which are challenging for numerical simulations. The objective of this study is to investigate the roles of the various processes and their interactions in multiphase flows in unsaturated soils under external wave loading, via theoretical time‐scale analysis and numerical simulations. A coupled geomechanics–multiphase flow model based on conservation laws is used. Theoretical analysis based on coupled and decoupled models demonstrates that transient and steady‐state responses are governed by pore pressure diffusion and saturation front propagation, respectively, and that the two processes are essentially decoupled. Numerical simulations suggest that the compressibility of the pore fluids and the deformation of the soil skeleton are important when the transient responses of the media are of concern, while the steady‐state responses are not sensitive to these factors. The responses obtained from the fully coupled numerical simulations are explained by a simplified time‐scale analysis based on coupled and decoupled models. Copyright © 2010 John Wiley & Sons, Ltd.     

Effects of the Soil Heat Flux Estimates on Surface Energy BalanceClosure over a Semi-Arid Grassland

Soil heat flux is important for surface energy balance (SEB), and inaccurate estimation of soil heat flux often leads to surface energy imbalance. In this paper, by using observations of surface radiation fluxes and soil temperature gradients at a semi-arid grassland in Xilingguole, Inner Mongolia, China from June to September 2008, the characters of the SEB for the semi-arid grassland were analyzed. Firstly, monthly averaged diurnal variations of SEB components were revealed. A 30-min forward phase displacement of soil heat flux (G) observed by a fluxplate at the depth of 5-cm below the soil surface was conducted and its effect on the SEB was studied. Secondly, the surface soil heat flux (Gs) was computed by using harmonic analysis and the effect of the soil heat storage between the surface and the fluxplate on the SEB was examined. The results show that with the 30-min forward phase displacement of observed G, the slope of the ordinary linear regression (OLR) of turbulent fluxes (H+LE) against available energy (Rn–G) increased from 0.835 to 0.842, i.e., the closure ratio of SEB increased by 0.7%, yet energy imclosure of 15.8% still existed in the SEB. When Gs, instead of G was used in the SEB equation, the slope of corresponding OLR of (H+LE) against (Rn–Gs) reached 0.979, thereby the imclosure ratio of SEB was reduced to only 2.1%.