土壤水非均匀流动的碘-淀粉染色示踪研究

在壤土和粉质黏土条件下采用碘-淀粉染色示踪方法分别开展了6个和4个不同土壤初始含水率和入渗水量的示踪试验,通过分析入渗深度和染色面积分布模式来研究入渗条件(土壤初始含水率、入渗水量和土壤质地)对土壤水非均匀流动特征的影响。研究表明:①土壤初始含水率升高时,土壤水入渗深度增大,但土壤初始含水率对染色面积分布模式无明显影响;②土壤水入渗深度随入渗水量的增大而增大,但入渗水量的增大对土壤水入渗深度增大的影响有限;③当入渗水量较小时,土壤水流运动以优先流为主;随着入渗水量的增大,土壤水流运动逐渐转变成以基质流为主的流动形态;④土壤水在细质地土壤(粉质黏土)中运动具有更高的非均匀程度,非均匀流通道总宽度较窄。     

土壤优先流运动的活动流场模型模拟和敏感性分析

在壤土和砂土条件下分别采用碘-淀粉染色示踪方法和亮蓝染色示踪方法各开展了2个染色示踪试验,分别采用活动流场模型和二域模型模拟计算了各试验入渗后染色区内的土壤含水率和溶质浓度分布,通过相对均方根误差分析评价了两个模型模拟预测优先流发展的有效性;此外,通过敏感性分析研究了不同入渗条件（土壤质地、入渗水量和土壤初始含水率）下活动流场模型模拟预测结果（入渗深度）对活动流场模型分形特征参数变化的敏感度。模型检验分析结果显示活动流场模型对土壤水入渗深度、入渗后染色区内土壤含水率和溶质浓度分布的预测精度要明显高于二域模型的模拟预测精度;活动流场模型较好的捕捉到了优先流运动整体的非均匀特征。敏感性分析结果显示,当降雨入渗水量和土壤初始含水率相同时,入渗深度对活动流场模型分形特征参数（γ）的敏感度随着γ的增大而增大;相同活动流场模型分形特征参数（γ）值条件下（即流动非均匀程度相同）,入渗深度对活动流场模型分形特征参数（γ）的敏感度随着入渗水量的增大和土壤初始含水率的升高而减小。     

入渗水量和试验尺度对土壤水非均匀流动的影响

[目的]研究入渗水量和试验尺度对土壤水非均匀流动(最大入渗深度和平均入渗深度、土壤水流运动非均匀程度和土壤优先流尺寸的影响)的影响,为农业灌溉施肥的高效利用、地下水环境保护以及流域水文过程预测等方面提供理论支撑。[方法]采用亮蓝染色示踪方法在砂土条件下开展了9组18个(每组2个重复)不同入渗水量(分别为1.25,2.5和5.0cm)和试验尺度(分别为0.25m×0.25m,0.5m×0.5m和1.0m×1.0m)的入渗试验。[结果](1)相同试验尺度条件下,土壤水最大入渗深度和平均入渗深度均随着入渗水量的增大而增大;土壤水流运动非均匀程度随入渗水量的增大而先增大后减小。(2)相同入渗水量条件下,当入渗水量较小时,土壤水最大入渗深度和土壤水流运动非均匀程度随着试验尺度的增大而增大;当入渗水量很大时,试验尺度的变化对最大入渗深度和土壤水流运动非均匀程度的影响不明显。(3)入渗水量的增加和试验尺度的增大都将使得优先流通道尺寸增大。[结论]入渗水量和试验尺度对土壤优先流运动有重要影响。     

土壤优先流运动的活动流场模型分形特征参数计算

活动流场模型分形特征参数控制着优先流的产生和发展,因此准确获得活动流场模型分形特征参数值对提高模型的模拟预测精度具有重要意义。该研究采用染色示踪方法,将优先流流场从流动背景中显示出来,通过数字图像分析和采样分析获得优先流流场和流场内土壤含水率的分布模式,根据活动流场模型本构方程拟合活动流场模型分形特征参数值;针对3种常见的入渗后染色区内土壤含水率分布模式,分别提供了相应的活动流场模型分形特征参数的计算方法。研究结果显示,1）由于入渗后土壤水重分布的影响,活动流场区域和染色区域在整个入渗深度范围并不完全重合,因此仅可选择活动流场与染色区域相重合深度范围内受土壤初始含水率影响较小的数据来拟合活动流场模型分形特征参数值;2）土壤质地对入渗后染色区内土壤含水率的分布模式有显著影响,细质地土壤中入渗后染色区土壤含水率沿入渗方向逐渐减小,粗质地土壤中入渗后染色区土壤含水率沿入渗方向先增大后减小。     

舟山沿海土壤机油污染的电阻率变化特征及物理影响因素研究

为探索舟山沿海土壤在机油污染情况下电阻率的变化规律及影响因素,通过室内配制污染土,采用单因素分析法研究机油污染砂土、砂质壤土和粉砂质壤土的电阻率变化与土壤质地类型、含水率、含油率及时间的关系。结果表明,土壤构成及含水率与含机油率均对土壤电阻率产生影响。未污染和污染土的电阻率随着含水率的增加呈幂指数下降趋势,其规律符合Archie公式,同时含水率对土壤电阻率的影响大于含油率的影响,并且在一定程度上机油的存在会削弱含水率的影响。5%含水率时,砂质壤土电阻率随含油率增加变化幅度最大,砂土次之,粉砂质壤土最小;15%含水率时,电阻率变化幅度为砂土>砂质壤土>粉砂质壤土。另外,低含水率(5%水)下机油污染粉砂质壤土电阻率随时间的变化波动减小,而高含水率(15%)下则是波动增加,即在经过18~20天后出现小幅下降随后继续增加。根据污染土壤电阻率变化特征可以判断是否存在机油污染,监测效果砂土和砂质壤土优于粉砂质壤土,低含水率条件监测效果优于高含水率条件。研究成果可为土壤机油污染的电阻率法监测提供理论支持。     

喷灌压力、喷嘴直径和喷洒历时对土壤入渗性能的影响

选取适宜的允许喷灌强度是避免喷灌条件下水土流失的重要措施之一。通过田间试验,测定了不同喷灌压力与喷嘴直径条件下土壤结皮厚度、表层土壤容重以及土壤入渗性能随喷洒历时的动态变化。试验处理包括2种质地土壤(壤质砂土和粉质黏壤土)、3个喷灌压力(103,138,172 kPa)和3个喷嘴直径(3.97,5.95,7.94 mm)。结果表明,喷灌压力越大,土壤结皮厚度和表层容重越小,土壤稳定入渗率越大。喷嘴直径越大,土壤结皮厚度和表层容重越大,土壤稳定入渗率越小。土壤结皮厚度和表层容重随着喷洒历时的增加显著增大,而土壤稳定入渗率则随之对数减小。喷洒动能强度是一个描述喷洒水滴对土壤入渗率影响的优选参数。研究结果为允许喷灌强度的确定提供科学依据。     

土壤初始含水率对优先流的影响简

土壤中优先流的存在会造成肥料利用率降低以及土壤深层甚至地下水污染.为确定土壤初始含水率对优先流的影响,通过室内含大孔隙土柱定水头入渗实验,以长武黑垆土及渭河砂土为研究对象,分析不同土壤初始含水率对优先流的影响.结果表明：1）相同质地,当土壤初始含水率大于其斥水性的峰值含水率时,湿润锋运移深度随着土壤初始含水率的增大而增大,反之湿润锋运移深度随着初始含水率的增加而减少,不同质地土壤湿润锋运移速度随土壤粉粒质量分数升高而减慢;2）累积入渗量受到土壤储水性及斥水性的双重影响,导致与湿润锋的运移趋势并不一致,黑垆土对照CK组及优先流O-C-30组（除7.04％含水率外）累积入渗量随着土壤初始含水率的增大而减小,渭河砂土CK组及O-C-30组均呈现2阶段特征,累积入渗量先随着土壤初始含水率的增加而减少,之后随土壤初始含水率的增加而增加;3）Kostiakov入渗模型能很好地拟合累积入渗量随时间的变化过程,适用于存在优先流的土壤入渗模拟.该研究结果对土壤优先流水分模拟具有一定的指导意义.     

一维马尔可夫链模拟黑河中游流域土壤质地垂向变异

黑河中游地区土壤剖面砂、黏层次相间排列的特点及部分区域剖面中出现的不透水层对土壤中的水分运动和溶质迁移具有重要影响。该研究调查了黑河中游100km2区域内土壤剖面的质地分层情况,运用一维嵌入马尔可夫链模型模拟该地区的土壤质地剖面。研究表明,研究区土壤剖面共出现砂土、壤质砂土、砂质壤土、壤土、黏质壤土、粉黏壤土和粉黏土7种质地类型,层次厚度呈对数正态分布。粉黏土未在表层土壤出现,而砂土在表层出现的概率明显高于其余质地类型,剖面某一质地层之下多出现细粒含量比其高的下一质地类型。剖面上相邻两质地层之间的转移具有明显的马尔可夫链特征(简称马氏性),且马氏链是平稳的。一维嵌入马尔可夫链模型能够较好地描述研究区土壤质地层次的垂向变化,剖面主要的质地层次组合为:砂—壤、壤—砂、壤—黏、黏—壤。对土壤质地剖面的定量模拟,可以为开展该区域土壤水循环、转化和溶质迁移等相关研究提供参考。     

聚丙烯酸钠对土壤水分入渗影响的模拟研究

采用一维垂直积水入渗法,分析了聚丙烯酸钠(SA)4种浓度(0,0.08%,0.2%与0.5%)处理、不同质地土壤间(砂土、壤土和黏土)累积人渗能力差异性和不同入渗模型的拟合效果.研究结果表明,初始入渗阶段3种土壤累积入渗能力较强,累积入渗量增加很快,随着人渗过程的进行累积入渗量曲线上升趋势变缓;3种土壤对照处理的累积人渗能力顺序为:壤土>砂土>黏土.3种土壤累积人渗能力都随着SA浓度的加大而减弱;随着聚丙烯酸钠浓度的增大(0.08%～0.5%),人渗能力都不同程度的降低,砂土累积人渗能力降低幅度最大.Horton模型对拟合各个处理的累积人渗过程效果最佳,Kostiakov模型和Philip模型拟合效果次之.     

土壤初始含水率对膜孔灌湿润体特征的影响

为了研究膜孔灌中土壤初始含水率对湿润体特征及累积入渗量的影响,首先通过室内试验验证HYDRUS模拟西安粉壤土膜孔灌湿润体形状以及含水率分布的可靠性,然后基于HYDRUS模型模拟在不同初始含水率条件下膜孔灌湿润体的变化过程。结果表明:基于HYDRUS模型模拟的累积入渗量和湿润锋运移距离与室内试验结果的R~2均接近1,标准偏差绝对值均小于10%,拟合良好,表明HYDRUS模型模拟入渗过程的可靠性。膜孔灌湿润锋形状可采用椭圆方程表示。当初始含水率较小(不大于0.1 cm~3/cm~3)时,湿润体半径的含水率分布可采用椭圆方程表示;从膜孔中心到湿润锋表面,随着初始含水率的增大,湿润体内的含水率梯度减小,湿润体半径的含水率分布曲线由椭圆曲线逐渐转变为平缓曲线。基于湿润体含水率分布规律建立了考虑初始含水率的累积入渗量模型,累积入渗量与湿润体半径的三次方呈正比,湿润体半径可表示为湿润锋水平运移距离和垂向运移距离的几何平均值;对于不同的膜孔半径(1～5cm),模型计算累积入渗量与HYDRUS模拟值的R~2为0.99,标准偏差绝对值小于10%;对于粉土、粉壤土和壤土,当初始体积含水率不大于0.25 cm~3/cm~3算累积入渗量与HYDRUS模拟值的R~2为0.99,标准偏差绝对值小于10%,结果表明该模型对不同土壤质地和膜孔半径的适用性良好;该模型在计算作物灌水需求量方面优于Kostiacov模型等传统的经验模型。该研究揭示了不同初始含水率下的膜孔灌湿润体特征,并建立了累积入渗量模型,可为膜孔灌灌溉水量的计算提供参考依据。     

黄淮海平原三种土壤中优势流现象的试验研究

在黄淮海平原选择三种不同质地类型的土壤（砂壤土、淤土、风沙土）,各设１．５ｍ×１．５ｍ两个小区,其中一个小区在试验开始前两癸灌水约５６ｋｇ,以获得不同的初始含水量。将１００ｍｍ含有染色剂亮蓝的水灌入小区,一天之后,开挖１ｍ×１ｍ的剖面拍照,进行图像分析处理,计算剖面中染色百分比随深度的变异,研究所形成的优势流的状况与土壤类型、土壤初始含水量的关系。结果表明：不同土壤发生优势流的程度不同,结构发育好     

Effect of non-linear adsorption on the transport behaviour of Brilliant Blue in a field soil

The food dye Brilliant Blue FCF (Color Index 42090) is often used as dye tracer in field studies for visualizing the flow pathways of water in soils. Batch studies confirmed findings of other researchers that non‐linear sorption is important for Brilliant Blue, especially at small concentrations (< 10 g l^?1 for our soil), and that retardation increases with decreasing concentrations as well as with increasing ionic strength of solutions. Therefore, it is not obvious if it can be used as an indicator for water flow paths as is often done. In this study, we compared the mobility of Brilliant Blue in a field soil (gleyic Luvisol) with that of bromide. Brilliant Blue and potassium bromide were simultaneously applied as a 6‐mm pulse on a small plot in the field, and the tracers were displaced with 89 mm of tracer‐free water using a constant intensity of 3.9 ± 0.2 mm hour^?1. Both tracer concentrations were determined on 144 soil cores taken from a 1 m × 1 m vertical soil profile. The transport behaviour differed in both (i) mean displacement and (ii) spatial concentration pattern. We found the retardation of Brilliant Blue could not be neglected and, in contrast to the bromide pattern, a pulse splitting was observed at the plough pan. Numerical simulations with a particle tracking code revealed that the one‐dimensional concentration profile of bromide was represented fairly well by the model, but the prediction of the double peak in the Brilliant Blue concentration profile failed. With additional assumptions, there were indications that Brilliant Blue does not follow the same flow paths as bromide. However, the question of Brilliant Blue taking the same flow pathways as bromide cannot be adequately answered by comparing both concentration distributions, because we look at two different transport distances due to the retardation of Brilliant Blue. It became obvious, however, that Brilliant Blue is not a suitable compound for tracing the travel time of water itself.     

Analysing flow patterns from dye tracer experiments in a forest soil using extreme value statistics

Preferential flow of water in soil is now recognized as a common phenomenon. It results in complex flow patterns that can be visualized by dye tracers and increases the risk of pollutants reaching greater depths. We analysed the behaviour of a risk index for vertical solute propagation based on extreme value theory. This risk index can be calculated from binary images of dye‐stained soil profiles and is defined as the form parameter of the generalized Pareto distribution. We did five tracer experiments with Brilliant Blue and iodide under changing initial (variable initial soil moisture) and experimental conditions (different irrigation rates). Our results indicate some persistence of the risk index against small changes of experimental conditions such as the irrigation rate. On the other hand, it seems to be affected by initial soil moisture. Comparisons of Brilliant Blue and iodide patterns show that the form parameter alone is not sufficient to estimate the risk of vertical solute propagation. Therefore we propose to combine the risk index with the scale parameter of the generalized Pareto distribution.     

引入分维理论评价土壤对优势流的敏感性

Food dye Brilliant Blue was introduced as the tracer in a dye-tracing experiment to obtain dye profile patterns of sandy loam soil,aeolian sandy soil,percolating paddy soil and permeable paddy soil.The dyed soil profiles were then photographed and the photos were scanned into a computer.Edited with certain software,only the dyed areas were left on the profile photos,which indicted the preferential flow paths for water and solute transport.Fractal dimensions of the dye patterns were calculated according to Arnold‘s function.Soil particle size distribution was analyzed by pipette method.The regression analysis showed that there was significant relationship between soil clay content and fractal dimension D of the dye pattern of soil profile.Based on the experiment results,the possibility of introducing fractal dimension to estimation of soil sensitivity to preferential flow is discussed.     

Dye tracer infiltration in the plough layer after straw incorporation

Straw residues may be accumulated in isolated zones in the plough layer after mouldboard ploughing. This may limit straw decomposition because of limited nitrogen availability depending on the prevalent water flow pathways induced by tillage. In this study, we used the food dye Brilliant Blue FCF (Color Index 42090) to make visible the flow pathways in the plough layer of a loamy soil after incorporation of 10 tonne of dry mass ha^− 1, and to get a qualitative and quantitative interpretation of the small-scale flow behaviour. Approximately 1.3 L of a 4 g L^− 1 Brilliant Blue solution was applied at a matric head of − 1 cm with an infiltrometer, equipped with a 25 cm-diameter disk. Horizontal cross-sections of 50 × 50 cm were photographed at 1.5 to 3.5 cm depth intervals. High-resolution spatial maps of Brilliant Blue concentration were derived from the scanned photographs, using separate calibration relationships between the measured Brilliant Blue concentrations and the color spectra and depth for the soil and the straw. Pronounced lateral dye movement was observed through the soil matrix owing to soil sorptivity. In cases of high sorptivity, the dye barely reached the depth of straw incorporation. Otherwise, enhanced preferential dye transport directed towards the incorporated straw was made visible and the stained water eventually ponded on the plough pan which induced lateral redistribution. Although the straw inclusions contributed to dye transport, they did not retain a substantial part of the applied dye mass, owing to the low density and surface area of straw. However, these inclusions may enhance the fast migration of potential pollutants such as nitrate or pesticides out of the tilled layer where much of the transformation occurs due to biological activity.     

Quantifying dye tracers in soil profiles by image processing

Developing and testing models for solute transport in the field requires experimental data on the spreading of solutes in the soil. Obtaining such data is costly, and a substantial part of the total costs is in the preparation and chemical analysis of the tracing compounds in the gathered samples. We developed a cheap method to quantify the concentration of the mobile dye tracer Brilliant Blue FCF from digitized photographs of stained soil profiles, and we have tested it in the field. Soil sampling and chemical analyses were necessary only to establish a calibration relation between the dye content and the colour of the soil. The digital images were corrected for geometrical distortions, varying background brightness, and colour tinges, and then they were analysed to determine the soil colour at sampling points in the profiles. The resident concentration of the dye was modelled by polynomial regression with the primary colours red, green, blue and the soil depth as explanatory variables. Concentration maps of Brilliant Blue were then computed from the digitized images with a spatial resolution of 1 mm. Validation of the technique with independent data showed that the method predicted the concentration of the dye well, provided the corrected images contained only the colours included in the calibration.     

Fate of heavy metals in a strongly acidic shooting‐range soil: small‐scale metal distribution and its relation to preferential water flow

To assess the mobility of Pb and associated metals in a highly contaminated shooting range soil (Losone, Ticino, Switzerland), we investigated the spatial distribution of the metals and their relation to preferential water flow paths. A 2.2 m² plot located 40 m behind the stop butt was irrigated with a solution containing bromide and Brilliant Blue, a slightly sorbing dye. A soil profile 50 cm in width was sampled down to 80 cm with a spatial resolution of 2.5 cm, resulting in 626 samples. Concentrations of elements (12 ≤ Z ≤ 92) were determined by energy‐dispersive Xray fluorescence spectrometry, and Brilliant Blue concentrations were determined with a chromameter. In the acidic (pH 3), organic matter‐rich, well drained Dystric Cambisol, maximum concentrations of 80.9 g kg^‐1 Pb, 4.0 g kg^‐1 Sb, and 0.55 g kg^‐1 Cu were measured in the topsoil. Within 40 cm soil depth, however, Pb, Sb, and Cu approached background concentrations of 23 mg kg^‐1, 0.4 mg kg^‐1, and 9.4 mg kg^‐1, respectively. The even horizontal distribution and the steep gradient along soil depth indicate tight metal binding in the topsoil, and a fairly homogeneous transport front. In contrast, water flow through the profile was highly heterogeneous. In the uppermost 20 cm, preferential flow was initiated by heterogeneous infiltration at the soil surface, but had no influence on metal distribution. Below 20 cm, however, preferential flow originated from larger tree roots, and metal concentrations were significantly elevated along these macropores. Spatial distributions of Pb, Sb, and Cu were similar, suggesting that all three metals are strongly retained in the topsoil and transported along preferential water flow paths in the subsoil.     

A morphological approach to understanding preferential flow using image analysis with dye tracers and X-ray Computed Tomography

A key problem facing soil physics and hydropedology at present is some of the standard theories of water flow in soils do not fully reflect the processes at the pore scale, and thus, cannot be adequately used for prediction. As such, examination of soil structure is vital for hydropedologists. Realisation that solutes move preferentially through soil into groundwaters has meant research in this area has increased in importance. This paper describes a multi-scale approach to analyse transport mechanisms using visualisation techniques. Chloride and Brilliant Blue tracers were applied to undisturbed soil cores to examine the physical and morphological properties associated with preferential flow in a range of soil types. Following collection of serial digital images, it was possible to examine and quantify the nature of active water flow mechanisms in terms of both dye-stained pathways and spatial distribution of dye concentration, using image analysis. Preferential flow linked to water potential and soil structural discontinuity was observed in all but the coarsest textured soil which conformed to uniform flow theory. A high level of variability in flow patterns was noted between the soil types. Such information as to how a soil dynamically re-wets is key for hydropedologists involved in applications such as pollution modelling. This is especially significant when considering a wetting mechanism, such as preferential flow, that cannot be adequately described by conventional soil physics.