镉在中性和碱性土柱中不同深度处的运移模拟

Human health has been potentially threatened by cadmium (Cd) contained in sewage irrigation water.Previous studies of Cd transport in soils were mainly conducted using small soil cores with pH values less than 6.The objectives of this study were to determine the parameters of the convection-dispersion equation (CDE) for Cd transport in relatively larger columns with neutral and alkaline soils,and to investigate the parameters' variability with depth.The soil columns were 50 cm in length and 12.5 cm in diameter.Ceramic suction lysimeters were buried at depths of 2.5,7.5,17.5,27.5,and 37.5 cm to abstract soil solution.Cd concentration in the soil solution samples were subsequently analyzed to obtain breakthrough curves (BTCs).Equilibrium and nonequilibrium models in CXTFIT program were used to estimate parameters of the CDE.The results suggested that both equilibrium and non-equilibrium models performed well in modeling Cd transport.The hydrodynamic dispersion coefficient (D) ranged from 0.18 to 10.70 cm 2 h 1,showing large differences among different depths.The retardation factor (R d) ranged from 25.4 to 54.7 and the standard deviation of R d value was lower than 30% of the mean value.Precipitation coefficient (R p) decreased consistently with increasing depth,varying from 1.000 × 10 10 to 0.661 h 1.Sensitivity tests showed that D was less sensitive than R d.These results would be helpful in understanding the transport and retention of Cd in non-acidic soils.     

石灰土壤原状土柱上硼和锌的迁移

Leaching of boron （B） and zinc （Zn） can be significant in some pedomorphic conditions, which can cause contamination of shallow groundwater and economic losses. Boron and Zn adsorption and transport was studied using 8.4 cm diameter × 28 cm long intact columns from two calcareous soil series with differing clay contents and vadose zone structures： Lyallpur soil series, clay loam （fine-silty, mixed, hyperthermic Ustalfic Haplargid）, and Sultanpur soil series, sandy loam （coarse-silty, mixed, hyperthermic Ustollic Camborthid）. The adsorption isotherms were developed by equilibrating soil with 0.01 tool L^-1 CaCl2 aqueous solution containing varying amounts of B and Zn and were fitted to the Langmuir equation. The B and Zn breakthrough curves were fitted to the two-domain convective-dispersive equation. At the end of the leaching experiment, 0.11 L 10 g L^-1 blue dye solution was also applied to each column to mark the flow paths. The Lyallpur soil columns had a slightly greater adsorption partition coefficient both for B and Zn than the Sultanpur soil columns. In the Lyallpur soil columns, B arrival was immediate but the peak concentration ratio （the concentration in solution at equilibrium/concentration applied） was lower than that in the Sultanpur soil columns. The breakthrough of B in the Sultanpur soil columns occurred after about 10 cm of cumulative drainage in both the columns; the rise in effluent concentration was fast and the peak concentration ratio was almost 1. Zinc leaching through the soil columns was very limited as only one column from the Lyallpur soil series showed Zn breakthrough in the effluent where the peak concentration ratio was only 0.05. This study demonstrates the effect of soil structure on B transport and has implications for the nutrient management in field soils.     

Measured and predicted transport of two S-triazine herbicides through soil columns

Effluent concentration of chloride and two pesticides (prometon and atrazine) were measured during column displacement experiments at two water flow rates. A constant suction of approximately 1300 Pa was maintained in the packed soil columns which were positioned vertically on top of a vacuum chamber enclosing an automatic fraction collector. Measured breakthrough curves (BTC's) were analyzed in terms of two solute transport models: the standard two-parameter convection-dispersion equation (CDE), and a four-parameter two-site/two-region nonequilibrium model (TRM). Calculations obtained with the TRM model were found to be in better agreement with measured BTC's than predictions using the CDE model. Column retardation factors for prometon and artrazine calculated from equilibrium batch sorption coefficients were comparable to those estimated from the observed BTC's only when the nonequilibrium TRM model was used.     

Mobility of the herbicide pyrithiobac through intact soil columns

The relative mobility of pyrithiobac [sodium 2-chloro-6-(4, 6-dimethoxypyrimidin-2-ylthio)benzoate], a new herbicide used for postemergence control of broadleaf weeds in cotton (Gossypium hirsutum), was evaluated and compared against that of bromide (Br(-)) tracer on four soils representative of cotton-growing regions using intact soil columns under saturated flow conditions. Pyrithiobac breakthrough curves were asymmetrical in shape with significant tailing and displaced to the left of 1 pore volume in the Houston Black clay (fine, montmorillonitic, thermic Udic Pellustert), Orelia fine sandy clay loam (fine-loamy, mixed, hyperthermic Typic Ochraqualfs), and Ships silty clay (very-fine, mixed, thermic Udic Chromustert) soils. Breakthrough of pyrithiobac in the Hidalgo sandy loam soil (fine-loamy, mixed, hyperthermic Typic Calciustoll) was delayed and more symmetrical, with peak pyrithiobac concentration reached after 1.2 pore volumes. The immobile pore water (IPW) fractions estimated from the Br(-) breakthrough curves ranged from 20 to 87% of total pore water. The IPW values demonstrated that soils with the greatest amount of IPW (Ships with IPW = 87.3%) exhibited the most rapid movement of pyrithiobac (peak concentration after 0.04 pore volume). The experimentally determined pyrithiobac breakthrough curves confirmed the high mobility of this herbicide in these alkaline and predominantly smectitic soils. These results indicate that pyrithiobac mobility was influenced by soil type and preferential flow processes when leached through intact soil columns.     

利用两区和两流区模型模拟氯离子在黑垆土中的迁移过程

以原状长武黑垆土为研究对象,针对不同输入方式、不同孔隙水流速条件下的氯离子溶质迁移过程,分析原状土壤中大孔隙存在条件下溶质迁移的特征,对两区模型和两流区模型的各拟合参数进行对比和分析。结果表明:大小脉冲输入条件下,氯离子浓度在试验土柱中均随穿透试验的进行呈现先增加而后降低的趋势,其峰值大小及其出现位置随孔隙水流速变化差异显著,且小脉冲输入时溶质穿透曲线较大脉冲时的不对称和拖尾现象更为显著;两流区模型相对于两区模型可较好描述氯离子的穿透曲线;两模型所得模拟参数随孔隙水流速的变化差异性较大。两模型拟合所得弥散度(λ)随孔隙水流速的变化趋势相同,但两流区模型拟合的弥散度大两区模型拟合值2个数量级。两模型拟合所得可动水含量(β)均随孔隙水流速的增加相应增加,而质量交换系数(α)值均极小。     

太湖地区水稻土的硝态氮穿透曲线及影响因素

用模拟土柱的方法对太湖地区三种主要水稻土的硝态氮垂直运移穿透曲线和影响因素进行了研究。结果表明:对同一种土壤而言,土层由上往下土层的出流速度减慢,硝态氮的穿透过程推迟。对这三种土壤而言,乌栅土的各土层的平均流量最小,穿透时间最长。在水分饱和和充分淋洗的条件下,硝态氮能完全穿透土柱,硝态氮的垂直运移过程主要受土壤粘粒含量的影响,随着粘粒含量的增加出流速度减小,硝态氮的穿透开始时间推迟,硝态氮的穿透曲线峰值变低、峰面分布变宽,穿透曲线变化平缓。所有的穿透曲线都显得不对称,粘粒含量越高的土壤越有明显的拖尾现象。土壤容重和有机质含量对硝态氮的垂直运移有很大的影响;容重越小、有机质含量越高,出流的流量就越大,穿透开始的时间也就越提前;有机质与穿透总时间有极显著的相关性;有机质含量越高,穿透总时间就越短。     

大肠杆菌在饱和砂质壤土中非平衡运移的CDE数学模型模拟

稳态、饱和条件下，以速度u=0．214cm min^-1和0．470cm min^-1进行混合置换实验，研究大肠杆菌在砂质壤土中的运移，并根据平衡与非平衡假设下，对流一弥散方程数学模型进行数值模拟。结果表明，含非平衡的双点吸附、且含有不可逆滞留项的模型能够较好地模拟大肠杆菌在砂质壤土中的运移。模拟和实验结果均表明，大肠杆菌的BTC（Breakthrough Curve穿透曲线）与示踪剂相比峰值明显降低，拖尾明显，且出现延迟，总的流出量也明显少于示踪剂；当水流速度由0．214cm min^-1增至0．470cm min^-1时。大肠杆菌BTC峰值由0．05增至0．2，且随速度的增加，滞留系数减小。     

A semi‐analytical model for solute transport in layered dual‐porosity media

The vulnerability of groundwater from chemical leaching through soil is a concern at some locations. Because measurements are laborious, time‐consuming, and expensive, simulation models are frequently used to assess leaching risks. But the significance of simulated solute movement through a layered soil is questionable if vertical homogeneity of physical soil properties has been assumed. In the present study, a semi‐analytical model for solute leaching in soils is presented. The model is relatively simple, but it does account for soil layers having different physical properties. The model includes the mobile‐immobile model (MIM) to describe one‐dimensional (1‐D) nonequilibrium, transient solute transport under steady‐state flow conditions. The MIM is rewritten as a second‐order differential equation and solved by a numerical scheme. Differing from fully analytical or fully numerical solutions, the new approach solves the differential equation numerically with respect to time and analytically with respect to distance. Numerical experiments for a single layered soil profile show that the semi‐analytical solution (SA‐MIM) is numerically stable for a wide range of parameter values. The accuracy of SA‐MIM predictions is comparable to that of analytical solutions. Numerical experiments for a multilayered profile indicate that the model correctly predicts effluent curves from finite layered soil profiles under steady‐state flow conditions. The SA‐MIM simulations with typical parameter values suggest that neglecting vertical heterogeneity of flow paths in a layered soil can lead to inaccurate prediction of soil‐solute leaching. The quality of predictions is generally improved if parameter estimates for the different soil layers are considered. However, the mobile‐immobile‐parameter estimates obtained in a number of previous studies may not be transferable to a field situation that is characterized by a slow and steady flow of water. Further field experiments to determine mobile‐immobile parameters under such conditions are desirable.     

土壤胶体在不同饱和度土壤介质中的释放与淋溶行为研究

采用饱和与非饱和填充土柱纵向淋溶研究方法,结合对流弥散模型方程（CDE）对穿透曲线的拟合计算,全面考察了土壤介质水饱和度、土壤水pH/离子强度、土壤孔隙水流速和土壤胶体颗粒大小对天然土壤胶体在实际土壤介质中释放、沉积迁移行为的影响。分别获取胶体扩散系数和阻滞因子值,定量说明实验中水化学、水动力学等条件的作用影响力。结果显示,介质不饱和条件不利于胶体的释放和淋溶;高pH和低离子强度条件对土壤胶体释放与迁移有利;淋溶过程的间断干扰,可以促使土壤胶体的增量淋溶释放;淋溶强度及胶体颗粒粒径大小,能够影响胶体穿透时间和穿透浓度峰值大小。     

Laboratory method for determining immobile soil water content and mass exchange coefficient

Preferential flow in soil can enhance the leaching of agricultural chemicals. In a number of studies it has been shown that the mobile‐immobile solute transport model (MIM) is a useful tool to characterize preferential flow. In the present study, a new laboratory method for determining the MIM parameters θ_m and θ_im (mobile and immobile water content), as well as α (mass transfer coefficient), is developed. The computations are uncomplicated and the method requires only simple equipment. It is applied to short, undisturbed soil columns. Measured values ranged from 0.11 to 0.27 for θ_im θ^—1 and from 0.015 h^—1 to 0.034 h^—1 for α for an Iowan soil (Nicollet silt loam). For two sandy Eutric Gleysols from Germany, low values for θ_im θ^—1 from 0.04 to 0.07 and from 0.001 h^—1 to 0.008 h^—1 for α were determined. Although the new method is a flow‐interruption technique, values for the Nicollet silt loam compare well with those from conventional leaching experiments. Values for the Eutric Gleysols agree with the observation that these soils were poorly structured. Because the new method does not assume negligible dispersion, it is applicable to a wider range of soils and boundary conditions than comparable approaches. We conclude that the new method provides parameter values that are suited to describe non‐equilibrium solute transport.     

Simulation of water flow and solute transport in free-drainage lysimeters and field soils with heterogeneous structures

Lysimeters are valuable for studying the fate and transport of chemicals in soil. Large‐scale field lysimeters are used to assess pesticide behaviour and radionuclide transport, and are assumed to represent natural field conditions better than laboratory columns. Field lysimeters are usually characterized by a free‐draining lower boundary. As a result, the hydraulic gradient is disrupted, and leachate cannot be collected until the bottom of the lysimeter becomes saturated. We compared heterogeneously structured, free‐drainage lysimeters and field soils with respect to water flow and solute transport. Numerical simulations were carried out in a two‐dimensional heterogeneous sandy soil under unsaturated water flow conditions with the CHAIN_2D code. Three different soil structures (isotropic, horizontal, and vertical) were generated, and Miller–Miller similitude was used to scale the hydraulic properties of the soil. The results showed that ponding occurs at the bottom of the lysimeter for the three soil structures and that it occurred faster and was more pronounced with the vertical structure (preferential flow effect). Breakthrough curves of a conservative solute (bromide) showed that solutes are moving faster in the field than in the lysimeters. Fewer differences between lysimeters and field soils were found with the horizontal soil structure than with the isotropic and vertical structures.     

An analytical solution for advance of solute front in soils

Based on the assumption that solute transport in a semi-infinite soil column or in a field soil profile can be described by the boundary-layer method, an analytical solution is presented for the advance of a solute front with time. The traditional convection-dispersion equation (CDE) subjected to two boundary conditions: 1) at the soil surface (or inlet boundary) and 2) at the solute front, was solved using a Laplace transformation. A comparison of resident concentrations using a boundary-layer method and an exact solution (in a semi-infinite-domain) showed that both were in good agreement within the range between the two boundaries. This led to a new method for estimating solute transport parameters in soils, requiring only observation of advance of the solute front with time. This may be corroborated visually using a tracer solution with marking-dye or measured utilizing time domain reflectometry (TDR). This method is applicable to both laboratory soil columns and field soils. Thus, it could be a step forward for modeling solute transport in field soils and for better understanding of the transport processes in soils.     

Simulation des Anionen-Transports in ungestörten Bodensäulen unter stationären Fließbedingungen

Simulation of anion transport in undisturbed soil cores under steady-state flow conditions Miscible displacement experiments with undisturbed soil columns were carried out in the laboratory. Objective of the experiments was to collect information about the transport of Cl^? and NO₃^? through field soils. The experiments were carried out with a forest soil and an agricultural soil. The flow velocity of the chloride and the nitrate solution was either 1 cm/day or 0,3 cm/day. Of each soil there were 5 replicates. The effluent of each of the columns was analyzed and the collected data were used for model calculations. It was found that the breakthrough curves of Cl^? and NO₃^? were similar in shape. For the well-aggregated forest soil the apparent diffusion coefficient was much larger than for the agricultural soil. For both soils practically all of the chloride could be recovered in the effluent, but for nitrate considerable losses within the soil column were noted. It was also found that the chloride breakthrough could be described with a simple convection-dispersion equation. However a reduction of the total pore space, accounting for anion exclusion, was needed. Furthermore it was observed that for nitrate an additional sink term in the convection-dispersion equation was needed to account for the observed nitrate losses. It appears that the transport of nitrate and chloride through the soils that were studied can be described mathematically, provided the anion exclusion space and the rate of nitrate losses are known. The nature of the anion exclusion and the nitrate transformation needs further study.     

The Fate of Cryptosporidium Parvum Oocysts in Reclaimed Water Irrigation-history and Non-history Soils Irrigated with Various Effluent Qualities

The present study aimed to look at the fate of protozoan parasite Cryptosporidium parvum oocysts applied through surface drip irrigation on reclaimed water irrigation-history and non-history sandy-loam (Hamra) soil columns. A new and simple isolation method for recovery of oocysts from soil samples was developed and used along this study. The new soil isolation method of oocysts is based on the “two phase separation method” formerly used to recover Clostridium perfringens spores from sediments and soil samples with minor modifications. The range recovery achieved by this method was 64–95% (mean 61.2?±?17.4). The objectives of the second part of this study were to investigate several physical and chemical factors governing transport and survival of C. parvum oocysts in sandy-loam soil columns by breakthrough curves. Comparison of fresh water and reclaimed water irrigation revealed that reclaimed water irrigated-history soil was more hydrophobic allowing water flow through channels with poor oocysts retention and fast flow. Examination of the organic matter effect (originating from reclaimed water irrigation) on oocysts breakthrough revealed that their soil infiltration increased. Calculations of oocysts concentration at different columns depths showed that most of the oocysts were retained in the first 5 cm of soil column. In the present study, comparing the two soil types (history and non-history of effluents irrigation) beside the surface electrostatic charge, one of the main elements found to affect oocysts infiltration and transport in soil columns was soil hydrophobicity caused by soluble organic matter originating from reclaimed water irrigation. Therefore, prior to application in soil irrigation, reclaimed water should be treated to high quality (i.e. membrane technology as the best option) to prevent enhanced transport of various pathogens through those irrigated soils.     

紫色土坡耕地土壤大孔隙流的定量评价

为阐明大孔隙丰富且孔径呈两极分化的紫色土坡耕地土壤大孔隙流的运移规律,通过室内土柱试验获取耕作层0~20 cm、非耕作层20~40 cm原状土柱和填装土柱的穿透曲线,分析饱和条件下土壤大孔隙流发生规律,并采用解析法CXTFIT软件拟合了水分优先运移参数,PFSP指标(大孔隙流引起的穿透曲线延展量与水动力弥散作用及两区作用引起的延展量的比值)定量评价土壤大孔隙流的贡献率。研究结果表明:1)以填装土柱水流为平衡基质流计算,耕作层0~20 cm原状土柱中大孔隙流的导水贡献率为66.2%~68.5%,而Br-累积淋出量占总淋出量的62.3%~66.1%。对于非耕作层20~40 cm,土壤大孔隙流导水贡献率为0.2%~1.7%,而Br-随大孔隙流运移的比例却达14.5%~20.5%。说明耕作层土壤中大孔隙流现象远比在非耕作层土壤中更为显著;2)PFSP值结果表明大孔隙流作用对穿透曲线延展量的贡献率最大,两区交换运移作用次之,水动力弥散作用的最小。即PFSP值越大,大孔隙流对总水流通量的贡献率越大。     

非均质饱和土壤盐分优先运移的随机模拟

在室内使用特殊形状的土柱隔板成功地填装了在水平横截面上呈“川”字型分布、由质地相差较大的两种土壤相间构成的非均质土柱。当土柱出流达到稳态后，灌入CaCl2溶液，监测土柱出流液的浓度动态。通过计算表征氯离子迁移时间随机特征的概率密度函数，对出流盐分的优先运移采用传递函数模型进行仿真，并对构成非均质土柱的两种均质土壤分别进行了条件类似的水盐入渗实验和模拟。在此基础上获得了参与氯离子输运的土壤水运移体积和可动体积以及土壤溶液中氯离子的体积平均驻留浓度。     

Modeling Surfactant LAS Influenced PAHs Migration in Soil Column

The objective of the present work was to study the transport of Polycyclic Aromatic Hydrocarbons (PAHs) in artificially contaminated and oil field soil columns while Linear Alkylbenzene Sulfonate (LAS) was used as leaching solvent. Through soil column leaching experiments, the through curves (BTCs) of the tracer (Br^?) and PAHs were obtained. The batch equilibrium experiments were conducted. The partition coefficients, and the retardation factor (R) were calculated. Nonlinear least-squares optimization approach was used to fit BTCs of Br^? and PAHs. The symmetrical BTC for Br^? was fitted using equilibrium convention-dispersion equation (CDE) model, and the physical and hydrodynamic parameter, i.e. the dispersion coefficient (D) was calculated. Based on these, the equilibrium and non-equilibrium CDE models were applied to fit the asymmetrical and tailing BTCs of PAHs. Results showed that two-site CDE model is better in fitting the observed data. The concentration distributions of PAHs with leaching time at different depth in the soil column were also estimated using two-site CDE model.     

土壤容重对红壤水分溶质运移特征影响的试验研究

通过测定南方第四纪红壤不同土壤容重条件下,土壤水分水平入渗过程和溴化钾溶液穿透曲线,分析了土壤容重对入渗速率、土壤饱和导水率的影响;根据CDE模型、运用CXTFIT软件对溴离子的穿透曲线(BTC)进行拟合.并求其运移参数.结果表明:土壤水分运移速率和土壤饱和导水率随着红壤容重增大而减少;随着容重增大,Br-的水动力弥散系数D增大,阻力因子R也增大.     

Cu/Pb/Zn/Cd在石英砂中的迁移实验及模拟

通过稳定流混合置换实验,研究了孔隙水流速和pH变化对Cu/Pb/Zn/Cd在石英砂中迁移行为的影响,获得了示踪剂Br-和Cu/Pb/Zn/Cd的穿透曲线(BTCs);并通过室内批量平衡实验获得吸附系数,进而计算出阻滞因子Rd。基于这些实验结果,借助CXTFIT2.1软件,用平衡CDE模型拟合了Br-的BTCs,得到了弥散系数D;在此基础上应用CDE非平衡模型拟合Cu/Pb/Zn/Cd在不同流速和pH条件下的BTCs,并预测了平行实验和不同埋深处Cu/Pb/Zn/Cd浓度的动态变化。结果表明,Cu/Pb/Zn/Cd的迁移能力随流速的增大而增强,而随pH的增大而降低;化学非平衡的两点模型能较好地模拟本文实验条件下Cu/Pb/Zn/Cd的迁移过程。