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%.     

A broad‐scale assessment of the effect of improved soil management on catchment baseflow index

Although the importance of sustainable soil management is recognized, there are many threats to soils including widespread soil structural degradation. This reduces infiltration through the soil surface and/or the percolation of water through the soil profile, with important consequences for crop yields, nutrient cycling and the hydrological response of catchments. This article describes a broad‐scale modelling approach to assess the potential effect that improved agricultural soil management, through reduced soil structural degradation, may have on the baseflow index (BFI) of catchments across England and Wales. A daily soil–water balance model was used to simulate the indicative BFI of 45 696 thirty‐year model runs for different combinations of soil type, soil/field condition, land cover class and climate which encapsulate the variability across England and Wales. The indicative BFI of catchments was then calculated by upscaling the results by spatial weighting. WaSim model outputs of indicative BFI were within the 95% confidence intervals of the national‐average BFI values given for the Hydrology of Soil Type (HOST ? ) classes for 26 of the 28 classes. At the catchment scale, the concordance correlation coefficient between the BFI from the WaSim model outputs and those derived from HOST was 0·83. Plausible improvements in agricultural soil/field condition produced modest simulated increases of up to 10% in the indicative BFI in most catchments across England and Wales, although for much of southern and northern England the increases were less than 5%. The results suggest that improved soil management might partially mitigate the expected adverse effects of climate change on baseflow to rivers. Healthy, well‐functioning soils produce many additional benefits such as better agricultural yields and reduced pollutant movement, so improved soil management should provide win‐win opportunities for society, agricultural systems and the environment and provide resilience to some of the expected environmental impacts of climate change. Copyright © 2011 John Wiley & Sons, Ltd.     

Implementing quality control on a random number stream to improve a stochastic weather generator

For decades, stochastic modellers have used computerized random number generators to produce random numeric sequences fitting a specified statistical distribution. Unfortunately, none of the random number generators we tested satisfactorily produced the target distribution. The result is generated distributions whose mean even diverges from the mean used to generate them, regardless of the length of run. Non‐uniform distributions from short sequences of random numbers are a major problem in stochastic climate generation, because truly uniform distributions are required to produce the intended climate parameter distributions. In order to ensure generation of a representative climate with the stochastic weather generator CLIGEN within a 30‐year run, we tested the climate output resulting from various random number generators. The resulting distributions of climate parameters showed significant departures from the target distributions in all cases. We traced this failure back to the uniform random number generators themselves. This paper proposes a quality control approach to select only those numbers that conform to the expected distribution being retained for subsequent use. The approach is based on goodness‐of‐fit analysis applied to the random numbers generated. Normally distributed deviates are further tested with confidence interval tests on their means and standard deviations. The positive effect of the new approach on the climate characteristics generated and the subsequent deterministic process‐based hydrology and soil erosion modelling are illustrated for four climatologically diverse sites. Copyright © 2007 John Wiley & Sons, Ltd.