The use of the PEARL model for assessing the migration of metribuzin in soil

The use of the PEARL mathematical model for describing the transfer of metribuzin in the soil profile and its migration to ground water was studied. The model was adjusted using experimental data on the water regime and the lysimetric runoff of soddy-podzolic soils in large lysimeters of Moscow State University. Without adjustment, the model successfully predicted the residual metribuzin concentrations in the topsoil but underestimated its content in the lower soil layers. After adjustment, the model more adequately predicted the herbicide concentrations in the lower horizons of the soil profile. The model adjustment improved the accuracy of predicting the accumulative removal of metribuzin beyond the soil profile by three orders of magnitude. The validation of the model with an independent data set showed satisfactory accuracy for forecasting the lysimetric runoff in soddy-podzolic soil and a rather low accuracy for predicting the distribution of the pesticide in the soil profile.     

Modeling pesticide leaching and dissipation in a mediterranean littoral greenhouse

The one-dimensional chromatographic flow model PEARL was used to simulate the movement of the insecticide imidacloprid and the fungicide procymidone through a greenhouse soil. The model was parametrized using measured and literature values of soil hydrological parameters. Soil water movement and soil temperature were reasonably well described by the model. The ability of PEARL to simulate the fate of imidacloprid and procymidone following four applications of each compound was evaluated against greenhouse data. Simulated imidacloprid residues in the 0-10 cm layer were in good agreement with measured data. Below 10 cm, the model overestimated imidacloprid remaining following the spray applications, whereas simulated residues following the chemigation applications were in reasonable agreement with measured data. Simulated residues of procymidone in the 0-10 cm layer were in general agreement with measured values. In the 10-20 cm layer, peaks in simulated concentrations occurred later than observed in the greenhouse. Predictions of procymidone residues below 20 cm were poor and underestimated compared to the measured data. For both pesticides, discrepancies between modeled and measured data in the 10-20 cm layer were attributed to the drip irrigation method used in the greenhouse. The model was unable to satisfactorily predict pesticide movement from the soil surface by irrigation water unless the scenario was modified to reflect the localized pattern of water application. Scenario analysis indicated that air boundary layer thickness is a key parameter for readily volatilized pesticides such as procymidone. This is of particular relevance to the greenhouse environment, where the boundary layer thickness may be greater than that in outdoor conditions.     

Quantitative Estimate of the Heterogeneity of Solute Fluxes Using the Dispersivity Length for Mathematical Models of Pesticide Migration in Soils

Mathematical models of solute movement in soils are used as management means for assessing the migration of agrochemicals and calculating environmental risks. In these models, the soil block includes the soil values of dispersivity length or solute diffusivity as one of the main parameters. Under laboratory conditions, this parameter was determined experimentally in soil columns by recording the dynamics of the effluent concentration and solving inverse problems. A direct experimental method of field determination based on the movement of marker solution was tested. For the prediction, risk calculation, and management of pesticide application using physically based mathematical models, the following stepwise procedure is recommended: (1) model parameterization based on experimental soil properties; (2) use of the field dispersivity length, which was 3–11 cm for the 0- to 40-cm layer, exceeded 12 cm for deeper fissured BEL and BT (40- to 60-cm) horizons with prismatic structure, and had a median value of 21 cm in layers below 60 cm in experiments on an agrosoddy-podzolic soil (Eutric Albic Glossic Retisol (Abruptic, Loamic, Aric) (WRB, 2014) from Moscow oblast.     

Review of modelling crop growth, movement of water and chemicals in relation to topsoil and subsoil compaction

Soil compaction influences crop growth, movement of water and chemicals in numerous ways. Mathematical modelling contributes to better understanding of the complex and variable effects. This paper reviews models for simulating topsoil and subsoil compaction effects. The need for including both topsoil and subsoil compaction results from still increasing compactive effect of vehicular pressure which penetrates more and more into the subsoil and which is very persistent. The models vary widely in their conceptual approach, degree of complexity, input parameters and output presentation. Mechanistic and deterministic models were most frequently used. To characterise soil compactness, the models use bulk density and/or penetration resistance and water content data. In most models root growth is predicted as a function of mechanical impedance and water status of soil and crop yield—from interactions of soil water and plant transpiration and assimilation. Models for predicting movement of water and chemicals are based on the Darcy/Richards one-dimensional flow equation. The effect of soil compaction is considered by changing hydraulic conductivity, water retention and root growth. The models available allow assessment of the effects of topsoil and subsoil compaction on crop yield, vertical root distribution, chemical movement and soil erosion. The performance of some models was improved by considering macro-porosity and strength discontinuity (spatial and temporal variability of material parameters). Scarcity of experimental data on the heterogeneity is a constraint in modelling the effects of soil compaction. Suitability of most models was determined under given site conditions. Few of the models (i.e. SIBIL and SIMWASER) were found to be satisfactory in modelling the effect of soil compaction on soil water dynamics and crop growth under different climate and soil conditions.     

田间原位试验分析长期机械作业下稻麦轮作地块土壤入渗性能

田间原位不同深度入渗试验是表达土壤分层状态、展示土层物理分异以及定量土壤剖面水功能变化的关键。为了探究不同深度水稻土的入渗能力及保水作用,该研究以华东稻麦轮作区小农户长期机械化耕整模式下代表地块的土层分异为目标,设计田间原位不同深度入渗试验。在试验地块内开挖7个不同深度的入渗坑并在坑底进行入渗试验,然后渗透48 h分层测取土壤含水率,研究不同坑底深度（坑深）土壤的入渗能力和入渗后各土层含水率的变化。结果表明,不同深度入渗试验准确表达了不同坑深土壤的水分入渗及土层持水分异,同时也能清晰地鉴别出犁底层所在位置和厚度,犁底层始于15 cm深,且耕作层与犁底层分异明显,耕作层平均紧实度为1 005.79 kPa,犁底层平均紧实度为1 910.73 kPa;土壤剖面分析表明,耕作层土壤形态疏松,根系分布稠密,犁底层土壤容重大,孔隙度小,透水性差,心土层土壤铁锰斑点较多,结构性差;土壤入渗参数随坑深的增加而减少,其中0～15 cm坑深范围内平均的平均入渗速率和累计入渗量分别为>20～30 cm的17.04倍和18.06倍;通过对比初始含水率和渗透48 h后含水率,得到坑深在15 cm范围内的...     

基于PEST的土壤-作物系统模型参数优化及灵敏度分析

农业生产管理系统模型输入参数多,参数率定过程十分耗时费力,大大限制了其推广应用。该研究以华北平原2 a的冬小麦-夏玉米田间试验观测数据为基础,使用PEST(parameter estimation)参数自动优化工具对土壤-作物-大气系统水热碳氮过程藕合模型(soil water heat carbon and nitrogen simulator,WHCNS)的土壤水力学参数、氮素转化参数和作物遗传参数进行自动寻优,同时计算分析模型参数的相对综合敏感度,并将优化结果与土壤实测水力学参数和试错法的模拟结果进行比较。参数敏感度分析结果表明,18个模型参数的相对综合敏感度较高,其中土壤水力学参数普遍具有较高的敏感度,以饱和含水率敏感度最高;作物参数中,作物生长发育总积温和最大比叶面积具有较高的综合敏感度;而氮素转化参数的敏感度远低于土壤水力学参数和作物参数。评价模型模拟效果的统计性指标(均方根误差、模型效率系数和一致性指数)表明,PEST法比实测水力学参数的模拟精度有所提高,其中土壤含水率、土壤硝态氮含量、作物产量和叶面积指数的均方根误差分别降低了61.8%、23.5%、73.6%和23.3%。同时PEST法比试错法对土壤水分和作物产量的模拟精度也有较大提高,但对土壤氮素和叶面积指数的模拟精度提高不明显。由于该方法大大节约了模型校准时间,在较短的时间内获得了明显高于试错法的模拟精度,因此PEST软件在WHCNS模型参数自动优化中是一个值得推广的工具。     

不同含盐土壤圆盘入渗特征试验

不同含盐土壤水分入渗特征是获得准确的土壤水力参数的基础。该文通过圆盘入渗试验,分析了4种土壤在5个（-1、-3、-6、-9和-12 cm）负水头下的入渗特征。结果表明,随着水头的减小,4种土壤的吸湿率线性减小,稳定入渗率和非饱和导水率呈不同程度减小。随土壤含盐量增加稳定入渗率和导水率呈增大规律。根据实测资料确定了不同负水头下非饱和导水率的Gardner指数模型参数,为盐渍化土壤水力参数的确定提供理论参考。     

利用一维代数模型分析微咸水入渗特征

由于淡水资源短缺,微咸水的开发利用成为缓解农业水资源短缺的重要措施,为了研究微咸水条件下的土壤入渗过程,进行了室内微咸水垂直一维入渗试验,分析了不同矿化度的微咸水累积入渗水量随时间的变化过程,并利用一维代数入渗模型模拟计算了土壤含水率剖面。结果表明：含水率模拟计算结果与实测值相吻合。同时利用一维代数入渗模型和Green-Ampt模型对试验资料进行处理,结果表明两模型均可以比较精确地描述微咸水入渗过程,因此两模型可以互相计算参数,但长历时入渗用Green-Ampt模型计算更加精确。     

秸秆覆盖对夏玉米农田土壤水分与热量影响的模拟研究

利用试验实测资料确定作物根系吸水模型和水热运移模型的参数,分别利用有限差分法对水热运移模型进行离散和迭代法对模型进行求解,通过模拟研究秸秆覆盖对夏玉米农田土壤水分和热量的影响。研究表明该模型能描述不同覆盖率条件下夏玉米生育期农田土壤水分和热量的变化规律,土壤含水率在0～40 cm逐渐减少,在40～100 cm逐渐增大;土壤温度在0～20 cm逐渐减少。不同覆盖率下土壤含水率和土壤温度的变化规律均为覆盖率150%＞覆盖率100%＞覆盖率50%＞不覆盖。可见秸秆覆盖节水和控温效果明显,为节水条件下的作物高产提     

水蚀风蚀交错带土壤剖面水力学性质变异

土壤剖面水力学性质的确定是土壤水分动态预测的基础。该文在水蚀风蚀交错区六道沟流域分别对居于坡中和坡上两块样地160 cm土层不同深度未扰动土壤的水分特征曲线进行了测定，将Van Genuchtens水分特征曲线模式与Mualem导水模式相结合，确定了两样地土壤剖面的水力学参数，对水力学参数在剖面的变化进行了分析。结果表明，土壤剖面饱和含水率、滞留含水率、进气吸力倒数和孔隙大小分布因子沿剖面变化不大，滞留含水率、进气吸力倒数属于中等程度变异，饱和含水率和孔隙大小分布指标属于弱变异，但经方差检验均不显著，说明该地区160 cm土壤剖面可以处理成均质剖面。     

Simulation of water permeability processes in chernozems of the Kamennaya Steppe

The physical and hydrophysical properties of ordinary chernozems from the Kamennaya Steppe were studied for the experimental support of a soil water infiltration model. Increased fissuring and density of the plow horizon were noted because of the dry weather conditions, which caused high infiltration values. The soil water retention curve parameters were calculated from the experimentally determined particle-size distribution, the bulk density of the soil, the solid phase density, the organic matter content, and the physicomechanical and soil-hydrological constants using the known pedotransfer functions and the Voronin secant method. A model experiment with different input parameters was performed. It was shown that the best experimental support of the model included parameters calculated by the Voronin method with the following approximation by the van Genuchten function.     

崩岗不同土层土壤水力学特性差异性分析

为研究崩岗不同土层土壤水力学特性的差异性,采用离心法测定不同土层土壤水分特征曲线,筛选出适合的土壤水分特征曲线拟合模型,结合统计模型,推求土壤的当量孔径分布、比水容量、非饱和导水率和扩散率,分析崩岗不同土层土壤水力学参数的变化规律。结果表明,崩岗土层从红土层到砂土层的变化过程中,土壤质地由黏土向砂土变化;Fredlund&Xing模型对崩岗土壤土水特征曲线拟合效果最好;参数θs、α、n随着质地变黏重逐渐减小;随着土层深度的增加,土壤的持水性能降低;土壤比水容量、非饱和导水率和扩散率受土壤质地和基质吸力的共同影响。在低吸力阶段,3个指标随基质吸力变化比较平缓,砂土层土壤比水容量和非饱和导水率最大,扩散率最小;而在高吸力阶段,砂土层土壤的这些指标降低较快,且低于其他土层,各层土壤间导水率和扩散率差异随着基质吸力的增加而增大。     

Comparison of measured and simulated changes in base cation amounts using a one-layer and a multi-layer soil acidification model

Biogeochemical processes describing nutrient cycling strongly affect the development of different soil horizons and soil acidification. However, in many soil acidification models these processes are not directly simulated, and spatially distributed soil properties are lumped into average values. This paper describes the extension of a one-layer soil acidification model (SMART) to a multi-layer model (SMARTEX) by including nutrient cycling, and the performance of both models regarding their ability for the assessment of the effects of acid deposition on forest soil. This was done by comparing simulation results of the two models with respect to changes in base cation amounts to historical data from sites in south Sweden. The results showed that the multi-layer model gives better results for the top soil, whereas the one-layer model gives better results for the rootzone and the total profile. Since the top soil is most important, the extended model is more suited to relate acid deposition to forest growth than the more general results of the one-layer model. However, the sensitivity of the multi-layer model to some parameters is large. A change in these parameters had strong effects on the separate layers, whereas the effects were cancelled out over the total profile. This sensitivity requires a careful parameterization of the model which might be troublesome to lack of data.     

作物根系层土壤水分模拟的经验-机理模型

To predict soil water variation in the crop root zone,a general exponential recession (GER) model was developed to depict the recession process of soil water storage.Incorporating the GER model into the mass balance model for soil water,a GER-based physicoempirical (PE-GER) model was proposed for simulating soil water variation in the crop root zone.The PE-GER model was calibrated and validated with experimental data of winter wheat in North China.Simulation results agreed well with the field experiment results,as well as were consistent with the simulation results from a more thoroughly developed soil water balance model which required more detailed parameters and inputs.Compared with a previously developed simple exponential recession (SER) based physicoempirical (PE-SER) model,PE-GER was more suitable for application in a broad range of soil texture,from light soil to heavy soil.Practical application of PE-GER showed that PE-GER could provide a convenient way to simulate and predict the variation of soil water storage in the crop root zone,especially in case of insufficient data for conceptual or hydrodynamic models.     

基于模型平均法的表层土壤含水率多模型综合反演

土壤水分是农作物和牧草生长的关键因素,也是影响全球气候变化和水循环的重要因子,因此准确监测土壤水分对促进农牧业可持续管理至关重要。基于Landsat 8 OLI多光谱影像,在表面生物物理特性和地形参数的基础上,结合野外实测土壤含水率,采用经验模型法分别构建了反距离加权法（Inverse Distance Weighting,IDW）、多元线性回归（Multivariable Linear Regression,MLR）、随机森林（Random Forest,RF）的土壤含水率反演模型,并采用模型平均法（Granger-Ramanathan,GR）进行组合反演,结果显示,通过模型平均法组合单一模型后在该研究区的土壤含水率反演中表现出了更好的适用性,模型平均法GR对于土壤含水率的估算效果佳,建模集与验证集R2≥0.88,均方根误差均不大于1.42%,且模型性能远优于反距离加权法、多元线性回归和随机森林。     

土壤持水曲线van Genuchten模型求参的Matlab实现

土壤持水曲线是研究土壤水力学性质必不可少的 ,在已经建立的众多数学模型中 ,vanGenuchten模型是目前运用最广泛的模型 ,而运用该模型的关键是其 4个参数的求解。为此 ,本文对同一组东北褐土的土壤水吸力和对应的土壤含水量数据较详细地介绍了Matlab软件的非线性拟合和非线性回归函数的运用 ,得出了该土壤vanGenuchten模型的 4个参数值 ,分别建立了该土壤的vanGenuchten模型 ,并利用Mat lab强大的绘图功能对它们进行了直观比较。最后运用方差分析和残差分析对该模型的计算值与实测数据进行了分析 ,结果表明 :非线性拟合和非线性回归函数求参结果的显著水平均达到p<0 0 0 0 1,残差平方和均小于 0 0 0 0 5 ,其中非线性回归函数的求参结果较非线性拟合好。因此 ,运用Matlab软件的非线性拟合和非线性回归函数对土壤持水曲线的vanGenuchten模型进行求参是切实可行的 ,从而为土壤学工作者寻求出了一条运用数值计算方法的新途径。     

利用PEARL模型评价农药渗透对地下水的污染

 运用荷兰开发的PEARL(PesticideEmissionAssessmentofRegionalandLocalArea)模型,模拟农药在作物—土壤中的运动情况,评价在不同浓度标准要求下,哪些农药会对地下水造成威胁,哪些地下水是安全的。针对我国农药施用现状,提出了关于农药污染地下水的研究,在农药登记中加入农药在地下水中的限量指标等建议。     

Physically based soil water characteristic curves function for soils with inner porosity

Many unimodal models have been developed to describe the soil water characteristic curves (SWCC). However, for soils with inner porosity the independent draining of the structural and matrix pores frequently results in two distinct air-entry values, which any single unimodal function does not reproduce adequately. In the presented study, a bimodal lognormal (BLN) SWCCs function have been developed for soils with inner porosity. The model involved two overlaying continua pore systems, the pore probability density functions of systems were assumed as lognormal distribution and them can be superposed to obtain the overall probability density function of the soil. The experimental SWCCs data were used to verify the proposed method. The fitted result showed that the BLN approach resulted in the best agreement between measurement and simulation (R² > 0.99, RSME < 0.66). Another advantage of proposed function was capable of simulating bimodal SWCC using parameters which can be related to physical properties of the structural soils. These physically based parameters could more intuitively analyse the measured data. The BLN approach can potentially be used as effective tool for identifying mechanical and hydraulic porosities in the structural mediums.