# Modelling Soil Water Retention for Weed Seed Germination Sensitivity to Water Potential

Authors: W. John Bullied; Paul R. Bullock; Rene C. Van Acker
Journal: Applied and Environmental Soil Science
 (2012)

---

## Abstract

Soil water retention is important for the study of water availability to germinating weed seeds. Six soil water retention models (Campbell, Brooks-Corey, four- and five-parameter van Genuchten, Tani, and Russo) with residual soil water parameter derivations were evaluated to describe water retention for weed seed germination at minimum threshold soil water potential for three hillslope positions. The Campbell, Brooks-Corey, and four-parameter van Genuchten model with modified or estimated forms of the residual parameter had superior but similar data fit. The Campbell model underestimated water retention at a potential less than −0.5 MPa for the upper hillslope that could result in underestimating seed germination. The Tani and Russo models overestimated water retention at a potential less than −0.1 MPa for all hillslope positions. Model selection and residual parameter specification are important for weed seed germination by representing water retention at the level of minimum threshold water potential for germination. Weed seed germination models driven by the hydrothermal soil environment rely on the best-fitting soil water retention model to produce dynamic predictions of seed germination.


---

## Body

## 1.  Introduction

The soil water retention characteristic (SWRC) is a basic hydrophysical property of the soil that relates the water content of soil water to its energy state [1]. The water content-potential function θ(ψ) is fundamental to the characterization of water holding capacity, water retention, and water flow in soil [2, 3]. The SWRC is necessary for modelling fluxes in soil water and is needed for germination studies where soil water is measured on a content basis.The timing of seed germination is a function of soil water potential [4]. As the soil dries, soil water potential is reduced, and it becomes increasingly difficult for seeds to imbibe water. At the minimum threshold (base) water potential, seeds do not imbibe sufficient water to initiate embryo growth and complete the germination process. The minimum threshold water potential at which germination ceases to occur in many agricultural weeds ranges from −0.1 to −1.5 MPa [5, 6]. Accurate representation of the SWRC over a wide range of water potential minimum thresholds is required for predictive modelling of seed germination.One of the greatest challenges in characterizing the SWRC for the shallow depth of the seedling recruitment zone (soil layer from which seeds germinate and emerge) across field topography is obtaining the parameters of the soil hydrological property. Determining the SWRC by direct measurement is time-consuming to obtain sufficient representation of a field due to spatial variability of soil properties [7, 8]. The ability of the soil to retain and transmit water is specific to soil type, and is affected by texture, pore-size and pore arrangement [1].Diverse SWRC often exist in soils along a hillslope as a result of variability in texture and pore-size distribution [9, 10]. Soil physical properties have been shown to reflect hillslope position with much greater variability between positions than within a position [11, 12]. Soil topographical variability in texture is related to differences between hillslope positions [13].The SWRC is expressed as a nonlinear function in which water content decreases from saturation to dryness in a sigmoidal manner with decreasing (more negative) water potential [14, 15]. Predicting θ(ψ) is particularly important for the minimum threshold water potential for seed germination within the range of soil water potential of −0.1 to −1.5 MPa. Anchoring the residual water content (θr) may represent the dry end of the soil water potential range to specify the point at which seed germination ceases to occur. The residual soil water content is the water content at the dry end of the SWRC at which the gradient dθ/dψ becomes zero and unsaturated hydraulic conductivity ceases [15]. Furthermore, some SWRC models represent dθ/dψ near the wet end of the range of soil water as a continuous function, whereas other SWRC models represent dθ/dψ as discontinuous. The objectives of this study were to (1) compare SWRC models to determine the best-fitting model for the measured soil water data across categorical hillslope positions (summit, backslope, and toeslope) and (2) to evaluate the SWRC models to describe weed seed germination sensitivity to soil water potential within specified water potential intervals for the shallow seedling recruitment zone at the categorical hillslope positions.


## 2.  Materials and Methods

### 2.1.  Experimental Site

The SWRC in the shallow seedling recruitment zone was evaluated in three hillslope positions (summit, backslope, and toeslope) on two hillslopes having opposing aspects in an annually cropped agricultural field at Graysville, MB, Canada. Each hillslope contained six replications that were arranged perpendicular to the hillslope gradient to maximize homogeneous soil conditions at each hillslope position.


### 2.2.  Soil Physical Properties

Soil samples were extracted from the 25–50 mm soil depth in three replicates by pressing 50 mm diameter by 25 mm deep rings into the soil. The 25–50 mm soil depth is representative of average weed seedling recruitment depth in a conventionally tilled field [16]. Furthermore, the 25–50 mm soil depth has been shown to have optimal hydrothermal properties for seed germination over a wide range of minimum threshold water potentials compared to shallower or deeper soil depths [17]. The 25–50 mm soil depth was also selected for analysis as part of a larger study on the soil water retention characteristics of the seedling recruitment zone [18].Subsamples of soil weighing 12 g were placed in shallow Plexiglas cylinders and saturated for 24 h. Volumetric water content at saturation was determined by placing the soil samples on saturated porous plates in covered ceramic suction cups to equilibrate saturation water content over 48 h. Volumetric water content was determined on desaturation at pressures of 0.01, 0.03, 0.05, and 0.1 MPa with a model 1600 0.5-MPa pressure plate extractor with a 0.1-MPa porous ceramic pressure plate and determined at pressures of 0.5 and 1.5 MPa with a model 1500 1.5-MPa pressure plate extractor with 0.5-MPa and 1.5-MPa porous ceramic pressure plates (Soil Moisture Equipment Corp., Santa Barbara, CA, USA) according to methods by [19]. At each pressure equilibrium, soils were removed from the ceramic plates and oven-dried at 105°C for 48 h.Soil particle size was analysed by the hydrometer method [20]. Organic matter content (OM) was determined by the loss on ignition procedure [21]. Bulk density (BD) was determined by averaging gravimetric soil water samples of a known volume in six replications from each hillslope position. Soil physical properties were analysed with the mixed procedure using statistical analysis systems (SAS) [22]. Mean separation of soil physical properties was analysed across all hillslope positions using Fisher's protected LSD.


### 2.3.  Soil Water Retention Models

Commonly used SWRC equations having a low number of parameters were evaluated for their ability to fit the data. The Brooks-Corey (BC) model [14] is a four-parameter equation in which water content is expressed as a power function of soil water pressure using the expression:(1)θ={θr+(θs-θr)[ψψe]-λforψ<ψe,θsforψ≥ψe,
where θ is volumetric water content, ψ is matric potential, θs is saturated volumetric water content, θr is residual soil water content when ψ is infinitely small, ψe is a curve-fitting parameter, known as the air entry matric suction, and λ is a pore-size distribution factor. The restriction θ=θs is applied for ψ≥ψe to prevent overestimation of water content at low matric suctions. The Campbell (CA) model [23] is a three-parameter power function to obtain water content from soil water potential:(2)θ={θs[ψψe]λforψ<ψe,θ=θsforψ≥ψe,
where ψe is the air entry matric suction, and λ is a curve-fitting parameter representing pore-size distribution where λ=-1/β, and β is an empirically derived constant. The Campbell equation is similar to the Brooks-Corey model but with θr=0. The van Genuchten (VG5) model [15] is an expression originally based on five parameters:(3)θ=θr+(θs-θr)[1+(αψ)n]-m,
where α is a curve-fitting parameter related to ψ-1 and the slope of the curve at inflection, and n and m are dimensionless curve-fitting parameters. The parameters are contingent on the pore-size distribution. The four-parameter van Genuchten (VG4) model is similar to the VG5 model except for m=1-1/n [15]. The Tani (TA) model [24] is a three-parameter model using the expression:(4)θ=θr+(θs-θr)[1+ψψo]exp⁡(-ψψo),
where ψo is the soil water potential at the inflection point on the curve. The Russo (RU) model [25] is a four-parameter model that produces the water conductivity-capillary potential relationship when integrated into Mualem’s model [26] model for relative hydraulic conductivity:(5)θ=θr+(θs-θr)[(1+0.5α|ψ|)exp⁡(-0.5α|ψ|)]2/(m+2),
where α is related to the width of the pore-size distribution and is interpreted as the inverse of the air entry matric suction, and m accounts for the dependence of the tortuosity and the correlation factors on the soil water content.


### 2.4.  Soil Water Retention Model Evaluation

The fitted SWRC for the 25–50 mm soil depth of the hillslope positions was determined by nonlinear analysis with the likelihood-based NLMIXED procedure using iterative optimization to compute the parameter estimates [22]. The values for θs were taken from measurements at saturation. Three derivations of θr were explored to determine the influence of θr on the dry range of the SWRC, since the minimum threshold water potential for seed germination of most weed species is within the dry range of water potential values (from −0.1 to −1.5 MPa). Values of θr derived from measurements at −1.5 MPa (the permanent wilting point) were designated θr(W). For practical purposes, θr can sufficiently be defined as the water content at a large negative value such as the permanent wilting point, despite further desorption of water with increasing suction [15]. The second derivation of θr was modified by a logarithmic equation describing the adsorption of water on soil in the dry range of the retention curve [27, 28]:(6)θ=[1-ln⁡(ψ)ln⁡(ψm)]θa,
where θa is a curve-fitting parameter representing the volumetric water content when ψ=1, and ψm is the matric suction at oven dryness, which is generally accepted to be 103 MPa. The modified form of θr is assigned θr(M). The third derivation of θr was estimated as one of the fitted parameters, referred to θr(E), which contained the restriction θr≥0. In all cases, the optimization results were improved by setting the measured value of θs as a constant parameter. Previous SWRC estimates have shown improvement by using known constant rather than optimized values of θs [29]. The BC and CA models were analysed as discontinuous dθ/dψ models, whereas the VG, TA, and RU models were analysed as continuous dθ/dψ models.The Akaike information criterion (AIC) [30] was used to evaluate the SWRC models. The AIC selected for the most parsimonious model [31, 32]. The AIC is a likelihood-based comparison of the data for model discrimination expressed as:(7)AIC=-2log⁡(likelihood)+2k,
where k is the number of estimable model parameters. The AIC selects for a model that fits well and has a low number of parameters. The AIC compares the data to a probability index, with a lower value of AIC indicating a better fitting model. The AIC is a relative ranking statistic; therefore, values are interpreted in terms of the magnitude of their differences among all models being considered.Model comparisons were facilitated using delta AIC(Δi) and Akaike weights (wi) [31, 33]. Delta AIC is a measure of the AIC differences of each model relative to the best model derived by:(8)Δi=AICi-AICmin⁡,
where AICi is the AIC value for model i, and AICmin is the smallest value of AIC in the set of candidate models. The relative likelihood of each SWRC model was realized by relative scaling of the models with wi [31]. Akaike weights provide an effective way to scale and interpret the Δi values by comparing the ratio of each model to the best model relative to the entire set of candidate models by:(9)wi=exp⁡(-Δi/2)∑r=1Rexp⁡(-Δr/2),
given a set of R models being evaluated. Akaike weights compare models on a scale of one (the sum of wi equals one) indicating the weight of evidence that model i is superior among the set of R models [31].Model inference was based on evidence ratios that evaluate the relative likelihood of model pairs. Evidence ratios are calculated as the ratio of Akaike weights(wi/wj), where i is the estimated best model in the set, and j indexes the remaining models in the set. Evidence ratios provide support for the fitted models as to which one is superior, by comparing models in a pairwise approach that is invariant to all models in the set except the i and j models [31].


### 2.5.  Weed Seed Germination Sensitivity to Soil Water Potential

The best-fitting curve according to Akaike weights from theθr derivations for each SWRC model was retained to explore differences in weed seed germination sensitivity to soil water potential across the hillslope positions. The range of measured water potential values (−0.01 to −1.5 MPa) was divided into four intervals (A: −0.01 to −0.1 MPa; B: −0.1 to −0.5 MPa; C: −0.5 to −1.0 MPa; D: −1.0 to −1.5 MPa) to represent groups of minimum threshold water potential for germination. The mean water content for the full range of measured water potential values and water potential intervals was calculated by numerical quadrature with the following integral using the average value of a function [34]:(10)θ̅=1b-a∫abθmdlog⁡ψ,
where θ̅ is the mean water content (cm3/cm3), and θm is the SWRC fitted to the measured water contents (cm3/cm3). The integration boundaries a and b define the integration interval over which the mean water content was derived. The log-normal distribution of ψ was calculated using log10ψ to avoid overweighting of the dry end of the SWRC.Mean water contents for the best-fittingθr derivation of each SWRC model were compared across hillslope positions by mean separation with Fisher’s LSD using mixed model analysis. Hillslopes were combined, and hillslope position by SWRC model interactive effect was explored. Mean water content of each SWRC model was compared to that of the best-fitting model.


## 2.1.  Experimental Site

The SWRC in the shallow seedling recruitment zone was evaluated in three hillslope positions (summit, backslope, and toeslope) on two hillslopes having opposing aspects in an annually cropped agricultural field at Graysville, MB, Canada. Each hillslope contained six replications that were arranged perpendicular to the hillslope gradient to maximize homogeneous soil conditions at each hillslope position.


## 2.2.  Soil Physical Properties

Soil samples were extracted from the 25–50 mm soil depth in three replicates by pressing 50 mm diameter by 25 mm deep rings into the soil. The 25–50 mm soil depth is representative of average weed seedling recruitment depth in a conventionally tilled field [16]. Furthermore, the 25–50 mm soil depth has been shown to have optimal hydrothermal properties for seed germination over a wide range of minimum threshold water potentials compared to shallower or deeper soil depths [17]. The 25–50 mm soil depth was also selected for analysis as part of a larger study on the soil water retention characteristics of the seedling recruitment zone [18].Subsamples of soil weighing 12 g were placed in shallow Plexiglas cylinders and saturated for 24 h. Volumetric water content at saturation was determined by placing the soil samples on saturated porous plates in covered ceramic suction cups to equilibrate saturation water content over 48 h. Volumetric water content was determined on desaturation at pressures of 0.01, 0.03, 0.05, and 0.1 MPa with a model 1600 0.5-MPa pressure plate extractor with a 0.1-MPa porous ceramic pressure plate and determined at pressures of 0.5 and 1.5 MPa with a model 1500 1.5-MPa pressure plate extractor with 0.5-MPa and 1.5-MPa porous ceramic pressure plates (Soil Moisture Equipment Corp., Santa Barbara, CA, USA) according to methods by [19]. At each pressure equilibrium, soils were removed from the ceramic plates and oven-dried at 105°C for 48 h.Soil particle size was analysed by the hydrometer method [20]. Organic matter content (OM) was determined by the loss on ignition procedure [21]. Bulk density (BD) was determined by averaging gravimetric soil water samples of a known volume in six replications from each hillslope position. Soil physical properties were analysed with the mixed procedure using statistical analysis systems (SAS) [22]. Mean separation of soil physical properties was analysed across all hillslope positions using Fisher's protected LSD.


## 2.3.  Soil Water Retention Models

Commonly used SWRC equations having a low number of parameters were evaluated for their ability to fit the data. The Brooks-Corey (BC) model [14] is a four-parameter equation in which water content is expressed as a power function of soil water pressure using the expression:(1)θ={θr+(θs-θr)[ψψe]-λforψ<ψe,θsforψ≥ψe,
where θ is volumetric water content, ψ is matric potential, θs is saturated volumetric water content, θr is residual soil water content when ψ is infinitely small, ψe is a curve-fitting parameter, known as the air entry matric suction, and λ is a pore-size distribution factor. The restriction θ=θs is applied for ψ≥ψe to prevent overestimation of water content at low matric suctions. The Campbell (CA) model [23] is a three-parameter power function to obtain water content from soil water potential:(2)θ={θs[ψψe]λforψ<ψe,θ=θsforψ≥ψe,
where ψe is the air entry matric suction, and λ is a curve-fitting parameter representing pore-size distribution where λ=-1/β, and β is an empirically derived constant. The Campbell equation is similar to the Brooks-Corey model but with θr=0. The van Genuchten (VG5) model [15] is an expression originally based on five parameters:(3)θ=θr+(θs-θr)[1+(αψ)n]-m,
where α is a curve-fitting parameter related to ψ-1 and the slope of the curve at inflection, and n and m are dimensionless curve-fitting parameters. The parameters are contingent on the pore-size distribution. The four-parameter van Genuchten (VG4) model is similar to the VG5 model except for m=1-1/n [15]. The Tani (TA) model [24] is a three-parameter model using the expression:(4)θ=θr+(θs-θr)[1+ψψo]exp⁡(-ψψo),
where ψo is the soil water potential at the inflection point on the curve. The Russo (RU) model [25] is a four-parameter model that produces the water conductivity-capillary potential relationship when integrated into Mualem’s model [26] model for relative hydraulic conductivity:(5)θ=θr+(θs-θr)[(1+0.5α|ψ|)exp⁡(-0.5α|ψ|)]2/(m+2),
where α is related to the width of the pore-size distribution and is interpreted as the inverse of the air entry matric suction, and m accounts for the dependence of the tortuosity and the correlation factors on the soil water content.


## 2.4.  Soil Water Retention Model Evaluation

The fitted SWRC for the 25–50 mm soil depth of the hillslope positions was determined by nonlinear analysis with the likelihood-based NLMIXED procedure using iterative optimization to compute the parameter estimates [22]. The values for θs were taken from measurements at saturation. Three derivations of θr were explored to determine the influence of θr on the dry range of the SWRC, since the minimum threshold water potential for seed germination of most weed species is within the dry range of water potential values (from −0.1 to −1.5 MPa). Values of θr derived from measurements at −1.5 MPa (the permanent wilting point) were designated θr(W). For practical purposes, θr can sufficiently be defined as the water content at a large negative value such as the permanent wilting point, despite further desorption of water with increasing suction [15]. The second derivation of θr was modified by a logarithmic equation describing the adsorption of water on soil in the dry range of the retention curve [27, 28]:(6)θ=[1-ln⁡(ψ)ln⁡(ψm)]θa,
where θa is a curve-fitting parameter representing the volumetric water content when ψ=1, and ψm is the matric suction at oven dryness, which is generally accepted to be 103 MPa. The modified form of θr is assigned θr(M). The third derivation of θr was estimated as one of the fitted parameters, referred to θr(E), which contained the restriction θr≥0. In all cases, the optimization results were improved by setting the measured value of θs as a constant parameter. Previous SWRC estimates have shown improvement by using known constant rather than optimized values of θs [29]. The BC and CA models were analysed as discontinuous dθ/dψ models, whereas the VG, TA, and RU models were analysed as continuous dθ/dψ models.The Akaike information criterion (AIC) [30] was used to evaluate the SWRC models. The AIC selected for the most parsimonious model [31, 32]. The AIC is a likelihood-based comparison of the data for model discrimination expressed as:(7)AIC=-2log⁡(likelihood)+2k,
where k is the number of estimable model parameters. The AIC selects for a model that fits well and has a low number of parameters. The AIC compares the data to a probability index, with a lower value of AIC indicating a better fitting model. The AIC is a relative ranking statistic; therefore, values are interpreted in terms of the magnitude of their differences among all models being considered.Model comparisons were facilitated using delta AIC(Δi) and Akaike weights (wi) [31, 33]. Delta AIC is a measure of the AIC differences of each model relative to the best model derived by:(8)Δi=AICi-AICmin⁡,
where AICi is the AIC value for model i, and AICmin is the smallest value of AIC in the set of candidate models. The relative likelihood of each SWRC model was realized by relative scaling of the models with wi [31]. Akaike weights provide an effective way to scale and interpret the Δi values by comparing the ratio of each model to the best model relative to the entire set of candidate models by:(9)wi=exp⁡(-Δi/2)∑r=1Rexp⁡(-Δr/2),
given a set of R models being evaluated. Akaike weights compare models on a scale of one (the sum of wi equals one) indicating the weight of evidence that model i is superior among the set of R models [31].Model inference was based on evidence ratios that evaluate the relative likelihood of model pairs. Evidence ratios are calculated as the ratio of Akaike weights(wi/wj), where i is the estimated best model in the set, and j indexes the remaining models in the set. Evidence ratios provide support for the fitted models as to which one is superior, by comparing models in a pairwise approach that is invariant to all models in the set except the i and j models [31].


## 2.5.  Weed Seed Germination Sensitivity to Soil Water Potential

The best-fitting curve according to Akaike weights from theθr derivations for each SWRC model was retained to explore differences in weed seed germination sensitivity to soil water potential across the hillslope positions. The range of measured water potential values (−0.01 to −1.5 MPa) was divided into four intervals (A: −0.01 to −0.1 MPa; B: −0.1 to −0.5 MPa; C: −0.5 to −1.0 MPa; D: −1.0 to −1.5 MPa) to represent groups of minimum threshold water potential for germination. The mean water content for the full range of measured water potential values and water potential intervals was calculated by numerical quadrature with the following integral using the average value of a function [34]:(10)θ̅=1b-a∫abθmdlog⁡ψ,
where θ̅ is the mean water content (cm3/cm3), and θm is the SWRC fitted to the measured water contents (cm3/cm3). The integration boundaries a and b define the integration interval over which the mean water content was derived. The log-normal distribution of ψ was calculated using log10ψ to avoid overweighting of the dry end of the SWRC.Mean water contents for the best-fittingθr derivation of each SWRC model were compared across hillslope positions by mean separation with Fisher’s LSD using mixed model analysis. Hillslopes were combined, and hillslope position by SWRC model interactive effect was explored. Mean water content of each SWRC model was compared to that of the best-fitting model.


## 3.  Results and Discussion

### 3.1.  Soil Physical Properties

The hillslope positions (summit, backslope, and toeslope) represented a categorical range in soil physical properties along the hillslope. Surface soil texture ranged from loamy fine sand to silty clay (Table1). Greater clay content in the toeslope was likely due to sorting of soil particles by hydrologic processes that would, over time, result in smaller soil particles, predominantly clay, to accumulate downslope [35]. Greater OM levels in the toeslope positions could be attributed to downhill transport and accumulation of OM. Bulk density did not vary greatly across the hillslope positions probably as a result of rotary tillage used for the accompanying experiment. The exception to relatively homogeneous BD among the hillslope positions was the SW toeslope which had higher OM that contributed to lower BD measured at that position.Table 1
Average measured soil physical properties for the shallow seedling recruitment zone of the hillslope positions.
Hillslope positionTextural classClaySiltSandOrganic matterBulk density(g/kg)(g/cm3)SW summitSilt loam56.7d746.7a196.7c35.3d1.12aSW backslopeSilt loam96.7c616.7c286.7b38.7cd1.04aSW toeslopeClay loam286.7a323.3e390.0a60.6a0.86bNE summitSilt loam83.3cd690.0b226.7c44.0bc1.09aNE backslopeSilt loam93.3c640.0bc266.7b34.7d1.13aNE toeslopeLoam156.7b463.3d380.0a46.0b1.05aSoil property means within a column, followed by different letters are significantly different atP≤0.05 according to Fisher’s protected LSD.


### 3.2.  Soil Water Retention Model Evaluation

The respective SWRC models (excluding CA) all exhibited better fit with modified(θr(M)) or estimated (θr(E)) forms of θr over that of measured (θr(W)) as indicated by lower values of AIC (Table 2). This may be an indication that modified or estimated forms of θr were better represented by values other than those measured at −1.5 MPa. Estimating or modifying θr as a free parameter provided better fit particularly in the dry range of the curve compared to θr(W) values (Figures 1–3). Models containing the modified form of θr retain the form of the original model in the wet range of the curve and the form of the logarithmic adsorption equation in the dry range [28]. The BC and VG models showed marginal improvement of fit when θr was modified rather than estimated. The TA and RU models, however, had considerably improved fit with θr(M) compared to θr(E) mainly in the dry range of the curve (Figures 2 and 3).Table 2
Average fitted parameter and AIC values for the soil water retention models across hillslopes.
Hillslope positionAverage parameter values (SE)AICgθr  derivationModelθsaθrbP1cP2dP3eP4f(cm3/cm3)(MPa)(dimensionless)Summit0.379AbsentCA0.000 (—)−0.0022 (0.0006)—3.29 (0.34)—−27.6MeasuredBC0.065 (—)−0.0042 (0.0006)—0.60 (0.05)—−28.5VG50.065 (—)—−0.0042 (0.0006)11.14 (0.04)0.054 (0.005)−26.5VG40.065 (—)—−0.0053 (0.0010)1.65 (0.07)—−28.1TA0.065 (—)—−0.0115 (0.0035)——−11.1RU0.065 (—)−0.0011 (<.0001)—16.66 (5.07)—−13.5ModifiedBC0.064 (0.011)−0.0042 (0.0008)—0.69 (0.16)—−30.1VG50.064 (0.008)—−0.0042 (<.0001)11.83 (<.01)0.059 (0.006)−28.1VG40.066 (0.011)—−0.0053 (0.0012)1.79 (0.20)—−29.4TA0.088 (0.005)—−0.0055 (0.0007)——−24.1RU0.087 (0.005)−0.0011 (<.0001)—5.43 (1.40)—−24.5EstimatedBC0.051 (0.015)−0.0037 (0.0007)—0.51 (0.10)—−29.4VG50.051 (0.011)—−0.0037 (<.0001)11.05 (<.01)0.046 (0.004)−27.4VG40.055 (0.015)—−0.0048 (0.0012)1.58 (0.13)—−28.5TA0.106 (0.014)—−0.0070 (0.0019)——−15.6RU0.098 (0.015)−0.0011 (<.0001)—9.93 (4.11)—−16.8Backslope0.391AbsentCA0.000 (—)−0.0023 (0.0004)—3.61 (0.23)—−32.1MeasuredBC0.073 (—)−0.0040 (0.0006)—0.54 (0.04)—−28.8VG50.073 (—)—−0.0040 (0.0006)21.21 (0.05)0.025 (0.002)−26.8VG40.073 (—)—−0.0053 (0.0009)1.59 (0.06)—−29.4TA0.073 (—)—−0.0140 (0.0008)——−10.2RU0.073 (—)−0.0011 (<.0001)—20.87 (6.05)—−12.8ModifiedBC0.061 (0.014)−0.0034 (0.0006)—0.51 (0.10)—−33.6VG50.061 (0.010)—−0.0034 (<.0001)9.68 (<.01)0.052 (0.005)−31.6VG40.066 (0.009)—−0.0044 (<.0001)1.58 (0.06)—−32.9TA0.102 (0.006)—−0.0055 (0.0009)——−21.9RU0.100 (0.007)−0.0011 (<.0001)—5.57 (1.94)—−22.5EstimatedBC0.042 (0.015)−0.0032 (0.0005)—0.39 (0.06)—−33.3VG50.042 (0.011)—−0.0032 (<.0001)9.38 (<.01)0.042 (0.003)−31.3VG40.049 (0.011)—−0.0042 (<.0001)1.45 (0.03)—−32.3TA0.120 (0.018)—−0.0078 (0.0031)——−13.6RU0.108 (0.019)−0.0011 (<.0001)—13.00 (5.66)—−15.4Toeslope0.417AbsentCA0.000 (—)−0.0014 (0.0003)—5.63 (0.30)—−34.5MeasuredBC0.120 (—)−0.0030 (0.0009)—0.41 (0.06)—−23.3VG50.120 (—)—−0.0030 (0.0009)19.63 (0.33)0.021 (0.003)−21.3VG40.120 (—)—−0.0039 (0.0013)1.44 (0.07)—−24.0TA0.120 (—)—−0.0177 (0.0057)——−8.1RU0.120 (—)−0.0011 (<.0001)—29.02 (10.26)—−10.5ModifiedBC0.000 (—)−0.0014 (0.0003)—0.18 (0.01)—−32.5VG50.000 (—)—−0.0014 (<.0001)14.34 (0.24)0.012 (<.001)−30.5VG40.000 (—)—−0.0016 (<.0001)1.18 (<.01)—−32.8TA0.149 (0.005)—−0.0039 (0.0008)——−23.4RU0.149 (0.005)−0.0011 (<.0001)—2.35 (1.23)—−23.4EstimatedBC0.000 (—)−0.0014 (0.0003)—0.18 (0.01)—−32.5VG50.000 (—)—−0.0014 (<.0001)14.45 (1.64)0.012 (0.001)−30.5VG40.000 (—)—−0.0016 (<.0001)1.18 (<.01)—−32.8TA0.182 (0.020)—−0.0065 (0.0019)——−11.3RU0.162 (0.028)−0.0011 (<.0001)—15.39 (10.95)—−12.4aθs is measured saturation water content.bθr is residual water content.cP1 is the air entry matric suction (ψe,ψe, and α-1 of the CA, BC, and RU models, resp.).dP2 is the capillary pressure at the inflection point on the water retention curve (α-1, α-1, and ψo  of the VG5, VG4, and TA models, resp.).eP3 is a dimensionless parameter (β,λ,n,n, and m of the CA, BC, VG5, VG4, and RU models, resp.).fP4 is a dimensionless parameter (m of the VG5 model).gA lower value of Akaike Information Criterion (AIC) indicates a better model fit.Fitted SWRC for soil water contents(θ) at soil matric potentials (ψ) for the (a) summit, (b) backslope, and (c) toeslope positions averaged across two hillslopes using the  Campbell (CA) model with absent θr parameter, and Brooks-Corey (BC), five- and four-parameter van Genuchten (VG5 and VG4), Tani (TA), and Russo (RU) models with θr constrained to measured values at −1.5 MPa. Error bars are SE for six water retention measurements.
(a)
Summit(b)
Backslope(c)
ToeslopeFitted SWRC for soil water contents(θ) at soil matric potentials (ψ) for the (a) summit, (b) backslope, and (c) toeslope positions averaged across two hillslopes using the Campbell (CA) model with absent θr parameter, and Brooks-Corey (BC), five- and four-parameter van Genuchten (VG5 and VG4), Tani (TA), and Russo (RU) models with θr modified by a logarithmic equation describing the adsorption of water on soil in the dry range of the retention curve. Error bars are SE for six water retention measurements.
(a)
Summit(b)
Backslope(c)
ToeslopeFitted SWRC for soil water contents(θ) at soil matric potentials (ψ) for the (a) summit, (b) backslope, and (c) toeslope positions averaged across two hillslopes using the Campbell (CA) model with absent θr parameter, and Brooks-Corey (BC), five- and four-parameter van Genuchten (VG5 and VG4), Tani (TA), and Russo (RU) models with θr estimated as one of the fitted parameters. Error bars are SE for six water retention measurements.
(a)
Summit(b)
Backslope(c)
ToeslopeThe evidence ratio provided a discrete comparison of the water retention models in a pairwise manner. The water retention models were compared to the BC model withθr(M), which was the best-fitting model based on wi derived from average AIC values (Table 3). An evidence ratio of 3 or less in relation to another model provides little evidence that model i is superior [31]. Accordingly, BC with θr(M) and θr(E), VG4 with θr(M) and θr(E), and CA models would serve nearly equally well in approximating the water retention information.Table 3
Akaike weights and evidence ratios for the soil water retention models.
θr derivation ModelAkaike weightaEvidence ratiobRankAbsentCA0.1431.44MeasuredBC0.01513.29VG50.00635.910VG40.01811.48TA0.00067, 045.116RU0.00019, 531.615ModifiedBC0.2031.01VG50.0752.76VG40.1691.23TA0.00285.312RU0.00372.211EstimatedBC0.1721.22VG50.0633.17VG40.1321.55TA0.00010, 542.014RU0.0005, 323.013Akaikea weight (wi) is based on average AIC value across all hillslope positions.bEvidence ratio (wi/wj) is based on wi from average AIC values, where model i is the best-fitting model (BC model with modified θr), and j indexes the remaining models in the set.The BC and VG5 models represented the SWRC similarly except at the wetter end of the curve, where dθ/dψ for the BC model was discontinuous (Figure 1). The BC model had lower AIC values (Table 2), which reflected parsimony in the number of parameters for the BC model. However, the CA and BC models represent the SWRC with additional complexity due to discontinuous dθ/dψ, where ψ≥ψe.The BC and VG4 models with θr(E) fit the SWRC nearly identical to one another between −0.01 and −1.5 MPa (Figure 3). The CA model was considered among the best models, due at least in part to its parsimonious equation. However, the CA model generally underestimated water retention less than approximately −0.5 MPa for the summit and backslope positions, compared to the BC and VG models.The VG5 model with modified and estimated forms of θr would be only slightly less adequate in describing the water retention relationship compared to the highest ranking five models (Table 3). The TA and RU models consistently had relatively high evidence ratios, particularly for the measured and estimated derivations of θr, indicating low prediction ability for the SWRC compared to the other models. Nonetheless, the TA and RU models were generally adequate to represent θ(ψ) at field capacity (−0.033 MPa). The TA and RU models also estimated water retention fairly well at −1.5 MPa, but only when the measured form of θr was used.The CA and BC models may provide unrealistic description of water retention at the wet end of the SWRC wheredθ/dψ is modelled as a discontinuous function [36, 37]. The continuous VG models contain an inflection point enabling better representation of water retention near saturation [36]. However, the ability of the VG4 model to predict dθ/dψ as a continuous function did not provide superior fit compared to the CA and BC models (Table 3). Regardless, the minimum threshold water potential for cessation of germination in most weed species is less than −0.1 MPa (Table 4).Table 4
Minimum threshold water potential for germination of weed species grouped by intervals of germination sensitivity to water potential.
Sensitivity interval (MPa) Minimum threshold (MPa)Common nameScientific nameReferenceB (−0.1 to −0.5)−0.10Common cockleburXanthium strumarium L.[6]−0.20Redroot pigweedAmaranthus retroflexus L.[38]C (−0.5 to −1.0)−0.50Common groundselSenecio vulgaris L.[39]−0.50Large crabgrassDigitaria sanguinalis (L.) Scop.[40]−0.60Wild oatAvena fatua L.[41]−0.64Common lambsquartersChenopodium album L.[42]−0.64VelvetleafAbutilon theophrasti Medik.[43]−0.69Yellow foxtailSetaria glauca (L.) Beauv.[5]−0.70Green foxtailSetaria viridis (L.) Beauv.[5]−0.80GoosegrassEleusine indica (L.) Gaertn.[44]−0.80Common ragweedAmbrosia artemisiifolia L.[45]−0.83Large crabgrassDigitaria sanguinalis (L.) Scop.[5]D (−1.0 to −1.5)−1.13Common chickweedStellaria media (L.) Vill.[46]−1.15Rigid ryegrassLolium rigidum Gaudin[47]−1.20Wild oatAvena fatua L.[48]−1.21GoosegrassEleusine indica (L.) Gaertn.[5]−1.50Wild oatAvena fatua L.[49]−1.50Perennial ryegrassLolium perenne L.[50]−1.53BlackgrassAlopecurus myosuroides Huds.[51]The VG4 model was superior to the VG5 model no matter which derivation of θr  was used (Table 3). The VG4 model has one less parameter making it more parsimonious than the VG5 model. The VG4 model may also provide a superior estimation of the water retention relationships compared to the VG5 model due to independence of the n and m parameters in the VG5 model, that can lead to distinctiveness problems in the estimation process, resulting in a less accurate description of the SWRC in the dry range [52].The SWRC for the coarse soils in the summit and backslope positions were best described by the BC model, based onwi values. The CA model had the best fit to the data for the soils with greater clay content in the toeslope position (Table 2).


### 3.3.  Weed Seed Germination Sensitivity to Soil Water Potential

The measured water potential interval (−0.01 to −1.5 MPa) was divided into four subintervals to represent sensitivity groups of minimum threshold water potential for germination of weed species (Table4). No species from the literature had a minimum threshold water potential within interval A (−0.01 to −0.1 MPa). The minimum threshold water potential for germination ranged from high sensitivity (interval B) to low sensitivity (interval D).The best-fitting SWRC for each model from the variousθr derivations was retained to evaluate water potential intervals for seed germination. Models evaluated were CA with no θr parameter, and BC, VG5,  VG4, TA and RU models with modified θr.The mean water content was averaged across hillslopes because water retention was not different between hillslopes (Table5). Similar mean water contents for the water potential intervals for the summit and backslope positions occurred as a result of few differences in the underlying soil physical properties between these two hillslope positions (Table 1). The upper hillslope positions differed from the toeslope as a result of differences in soil physical properties that influence the SWRC. The summit and backslope positions contained silt loam soils, whereas the toeslopes had loam and clay loam soils with greater OM content.Table 5
Mean water content for water potential intervals based on weed seed germination sensitivity to water potential.
Mean water contentWater potential interval (MPa)Main effectFullABCD−0.01 to −1.5−0.01 to −0.1−0.1 to −0.5−0.5 to −1.0−1.0 to −1.5(cm3/cm3)HillslopeSW0.1580.2020.1320.1080.098NE0.1500.1900.1260.1020.092Hillslope positionSummit0.124b0.164b0.099b0.080b0.072bBackslope0.141b0.185b0.114b0.091b0.082bToeslope0.196a0.239a0.175a0.145a0.132aModelCA0.153c0.202a0.125d0.095c0.083cBC0.153bc0.198b0.126c0.100b0.090bVG50.153bc0.199b0.126c0.100b0.089bVG40.153b0.198b0.126c0.101b0.091bTA0.156a0.189d0.137a0.118a0.109aRU0.156a0.191c0.135b0.116a0.108aSource of variationP>FHillslope (H)0.2160.1890.3130.2430.222Position (P)<.001<.001<.001<.001<.001H×P0.3600.4140.1610.0950.078Model (M)<.001<.001<.001<.001<.001H×M0.001<.001<.0010.0030.013P×M<.0010.0100.429<.001<.001H×P×M<.0010.4310.6820.4110.309Means within a column and main effect, followed by different letters are significantly different atP≤0.05 according to Fisher’s protected LSD.The interactions for hillslope position by SWRC model are shown in Figure4. The VG4 models did not differ from the BC model for any water potential interval in either the upper or lower hillslope. The CA model overestimated water retention between −0.01 to −0.1 MPa for the upper hillslope (Figures 4(a) and 4(b)). Since interval A does not contain minimum threshold water potential for weed seed germination, deviations from the measured water retention data in interval A would not be expected to influence seed germination of any of the species listed in Table 4. Deviations in water retention by the CA model were confined to the upper hillslope.Fitted SWRC for soil water contents(θ) at soil matric potentials (ψ) for the (a) summit, (b) backslope, and (c) toeslope positions using the Campbell (CA) model with absent θr parameter, and Brooks-Corey (BC), five- and four-parameter van Genuchten (VG5 and VG4), Tani (TA), and Russo (RU) models with θr modified as one of the fitted parameters. Hillslopes are averaged. Error bars are SE for six water retention measurements. Model* represents SWRC models that have significantly different water retention within a water potential interval compared to the BC model. Water potential (ψ) intervals A: −0.01 to −0.1 MPa; B: −0.1 to −0.5 MPa; C: −0.5 to −1.0 MPa; D: −1.0 to −1.5 MPa represent groups of minimum threshold water potential sensitivity for germination.
(a)
Summit(b)
Backslope(c)
ToeslopeThe CA model underestimated water retention at a potential less than −0.5 MPa for the summit and backslope positions (Figures4(a) and 4(b)). An underestimation in water retention less than −0.5 MPa would result in an underestimation in predicted germination for weeds with minimum threshold water potential less than −0.5 MPa. However, because soil water is generally measured in the field on a content basis, and SWRC are used to translate water content to a potential basis, the CA model would predict higher water potential for a given field water content compared to the BC or VG models. This case would overestimate germination of weeds that have minimum threshold water potential less than −0.5 MPa. Thus, germination would be predicted for species that would otherwise not occur in the upper hillslope. For example, at the summit, soil with a water content of 0.07 cm3/cm3 would have water potentials of −1.08 MPa and −0.57 MPa predicted by the BC and CA models, respectively. The BC model would predict germination of weeds only for interval D in Table 4, whereas the CA model would predict germination of weeds for interval D and almost all of interval C.The mean water content of the TA and RU models differed from that of the BC model across all water potential intervals for both upper and lower hillslope positions (Figure4). Water retention was underestimated by the TA and RU models for the −0.01 to −0.1 MPa interval, which would not be expected to influence modelling of germination for weeds listed in Table 4. However, overestimation of water retention by the TA and RU models for soils drier than −0.1 MPa would result in overestimation of seed germination at specific potentials. On a field water content basis, this would result in underestimation of seed germination by all species listed in Table 4. For instance, at the toeslope position, a soil water content of 0.15 cm3/cm3 would have a predicted water potential of −0.45 MPa with the BC model and −0.94 MPa with the TA model (Figure 4(c)). The BC model would, therefore, predict germination of weeds for intervals C and D, whereas the TA model would predict germination of weeds for only interval D. Predictions of weed seed germination from a range of minimum threshold water potentials are clearly dependent on the fit of the SWRC model.


## 3.1.  Soil Physical Properties

The hillslope positions (summit, backslope, and toeslope) represented a categorical range in soil physical properties along the hillslope. Surface soil texture ranged from loamy fine sand to silty clay (Table1). Greater clay content in the toeslope was likely due to sorting of soil particles by hydrologic processes that would, over time, result in smaller soil particles, predominantly clay, to accumulate downslope [35]. Greater OM levels in the toeslope positions could be attributed to downhill transport and accumulation of OM. Bulk density did not vary greatly across the hillslope positions probably as a result of rotary tillage used for the accompanying experiment. The exception to relatively homogeneous BD among the hillslope positions was the SW toeslope which had higher OM that contributed to lower BD measured at that position.Table 1
Average measured soil physical properties for the shallow seedling recruitment zone of the hillslope positions.
Hillslope positionTextural classClaySiltSandOrganic matterBulk density(g/kg)(g/cm3)SW summitSilt loam56.7d746.7a196.7c35.3d1.12aSW backslopeSilt loam96.7c616.7c286.7b38.7cd1.04aSW toeslopeClay loam286.7a323.3e390.0a60.6a0.86bNE summitSilt loam83.3cd690.0b226.7c44.0bc1.09aNE backslopeSilt loam93.3c640.0bc266.7b34.7d1.13aNE toeslopeLoam156.7b463.3d380.0a46.0b1.05aSoil property means within a column, followed by different letters are significantly different atP≤0.05 according to Fisher’s protected LSD.


## 3.2.  Soil Water Retention Model Evaluation

The respective SWRC models (excluding CA) all exhibited better fit with modified(θr(M)) or estimated (θr(E)) forms of θr over that of measured (θr(W)) as indicated by lower values of AIC (Table 2). This may be an indication that modified or estimated forms of θr were better represented by values other than those measured at −1.5 MPa. Estimating or modifying θr as a free parameter provided better fit particularly in the dry range of the curve compared to θr(W) values (Figures 1–3). Models containing the modified form of θr retain the form of the original model in the wet range of the curve and the form of the logarithmic adsorption equation in the dry range [28]. The BC and VG models showed marginal improvement of fit when θr was modified rather than estimated. The TA and RU models, however, had considerably improved fit with θr(M) compared to θr(E) mainly in the dry range of the curve (Figures 2 and 3).Table 2
Average fitted parameter and AIC values for the soil water retention models across hillslopes.
Hillslope positionAverage parameter values (SE)AICgθr  derivationModelθsaθrbP1cP2dP3eP4f(cm3/cm3)(MPa)(dimensionless)Summit0.379AbsentCA0.000 (—)−0.0022 (0.0006)—3.29 (0.34)—−27.6MeasuredBC0.065 (—)−0.0042 (0.0006)—0.60 (0.05)—−28.5VG50.065 (—)—−0.0042 (0.0006)11.14 (0.04)0.054 (0.005)−26.5VG40.065 (—)—−0.0053 (0.0010)1.65 (0.07)—−28.1TA0.065 (—)—−0.0115 (0.0035)——−11.1RU0.065 (—)−0.0011 (<.0001)—16.66 (5.07)—−13.5ModifiedBC0.064 (0.011)−0.0042 (0.0008)—0.69 (0.16)—−30.1VG50.064 (0.008)—−0.0042 (<.0001)11.83 (<.01)0.059 (0.006)−28.1VG40.066 (0.011)—−0.0053 (0.0012)1.79 (0.20)—−29.4TA0.088 (0.005)—−0.0055 (0.0007)——−24.1RU0.087 (0.005)−0.0011 (<.0001)—5.43 (1.40)—−24.5EstimatedBC0.051 (0.015)−0.0037 (0.0007)—0.51 (0.10)—−29.4VG50.051 (0.011)—−0.0037 (<.0001)11.05 (<.01)0.046 (0.004)−27.4VG40.055 (0.015)—−0.0048 (0.0012)1.58 (0.13)—−28.5TA0.106 (0.014)—−0.0070 (0.0019)——−15.6RU0.098 (0.015)−0.0011 (<.0001)—9.93 (4.11)—−16.8Backslope0.391AbsentCA0.000 (—)−0.0023 (0.0004)—3.61 (0.23)—−32.1MeasuredBC0.073 (—)−0.0040 (0.0006)—0.54 (0.04)—−28.8VG50.073 (—)—−0.0040 (0.0006)21.21 (0.05)0.025 (0.002)−26.8VG40.073 (—)—−0.0053 (0.0009)1.59 (0.06)—−29.4TA0.073 (—)—−0.0140 (0.0008)——−10.2RU0.073 (—)−0.0011 (<.0001)—20.87 (6.05)—−12.8ModifiedBC0.061 (0.014)−0.0034 (0.0006)—0.51 (0.10)—−33.6VG50.061 (0.010)—−0.0034 (<.0001)9.68 (<.01)0.052 (0.005)−31.6VG40.066 (0.009)—−0.0044 (<.0001)1.58 (0.06)—−32.9TA0.102 (0.006)—−0.0055 (0.0009)——−21.9RU0.100 (0.007)−0.0011 (<.0001)—5.57 (1.94)—−22.5EstimatedBC0.042 (0.015)−0.0032 (0.0005)—0.39 (0.06)—−33.3VG50.042 (0.011)—−0.0032 (<.0001)9.38 (<.01)0.042 (0.003)−31.3VG40.049 (0.011)—−0.0042 (<.0001)1.45 (0.03)—−32.3TA0.120 (0.018)—−0.0078 (0.0031)——−13.6RU0.108 (0.019)−0.0011 (<.0001)—13.00 (5.66)—−15.4Toeslope0.417AbsentCA0.000 (—)−0.0014 (0.0003)—5.63 (0.30)—−34.5MeasuredBC0.120 (—)−0.0030 (0.0009)—0.41 (0.06)—−23.3VG50.120 (—)—−0.0030 (0.0009)19.63 (0.33)0.021 (0.003)−21.3VG40.120 (—)—−0.0039 (0.0013)1.44 (0.07)—−24.0TA0.120 (—)—−0.0177 (0.0057)——−8.1RU0.120 (—)−0.0011 (<.0001)—29.02 (10.26)—−10.5ModifiedBC0.000 (—)−0.0014 (0.0003)—0.18 (0.01)—−32.5VG50.000 (—)—−0.0014 (<.0001)14.34 (0.24)0.012 (<.001)−30.5VG40.000 (—)—−0.0016 (<.0001)1.18 (<.01)—−32.8TA0.149 (0.005)—−0.0039 (0.0008)——−23.4RU0.149 (0.005)−0.0011 (<.0001)—2.35 (1.23)—−23.4EstimatedBC0.000 (—)−0.0014 (0.0003)—0.18 (0.01)—−32.5VG50.000 (—)—−0.0014 (<.0001)14.45 (1.64)0.012 (0.001)−30.5VG40.000 (—)—−0.0016 (<.0001)1.18 (<.01)—−32.8TA0.182 (0.020)—−0.0065 (0.0019)——−11.3RU0.162 (0.028)−0.0011 (<.0001)—15.39 (10.95)—−12.4aθs is measured saturation water content.bθr is residual water content.cP1 is the air entry matric suction (ψe,ψe, and α-1 of the CA, BC, and RU models, resp.).dP2 is the capillary pressure at the inflection point on the water retention curve (α-1, α-1, and ψo  of the VG5, VG4, and TA models, resp.).eP3 is a dimensionless parameter (β,λ,n,n, and m of the CA, BC, VG5, VG4, and RU models, resp.).fP4 is a dimensionless parameter (m of the VG5 model).gA lower value of Akaike Information Criterion (AIC) indicates a better model fit.Fitted SWRC for soil water contents(θ) at soil matric potentials (ψ) for the (a) summit, (b) backslope, and (c) toeslope positions averaged across two hillslopes using the  Campbell (CA) model with absent θr parameter, and Brooks-Corey (BC), five- and four-parameter van Genuchten (VG5 and VG4), Tani (TA), and Russo (RU) models with θr constrained to measured values at −1.5 MPa. Error bars are SE for six water retention measurements.
(a)
Summit(b)
Backslope(c)
ToeslopeFitted SWRC for soil water contents(θ) at soil matric potentials (ψ) for the (a) summit, (b) backslope, and (c) toeslope positions averaged across two hillslopes using the Campbell (CA) model with absent θr parameter, and Brooks-Corey (BC), five- and four-parameter van Genuchten (VG5 and VG4), Tani (TA), and Russo (RU) models with θr modified by a logarithmic equation describing the adsorption of water on soil in the dry range of the retention curve. Error bars are SE for six water retention measurements.
(a)
Summit(b)
Backslope(c)
ToeslopeFitted SWRC for soil water contents(θ) at soil matric potentials (ψ) for the (a) summit, (b) backslope, and (c) toeslope positions averaged across two hillslopes using the Campbell (CA) model with absent θr parameter, and Brooks-Corey (BC), five- and four-parameter van Genuchten (VG5 and VG4), Tani (TA), and Russo (RU) models with θr estimated as one of the fitted parameters. Error bars are SE for six water retention measurements.
(a)
Summit(b)
Backslope(c)
ToeslopeThe evidence ratio provided a discrete comparison of the water retention models in a pairwise manner. The water retention models were compared to the BC model withθr(M), which was the best-fitting model based on wi derived from average AIC values (Table 3). An evidence ratio of 3 or less in relation to another model provides little evidence that model i is superior [31]. Accordingly, BC with θr(M) and θr(E), VG4 with θr(M) and θr(E), and CA models would serve nearly equally well in approximating the water retention information.Table 3
Akaike weights and evidence ratios for the soil water retention models.
θr derivation ModelAkaike weightaEvidence ratiobRankAbsentCA0.1431.44MeasuredBC0.01513.29VG50.00635.910VG40.01811.48TA0.00067, 045.116RU0.00019, 531.615ModifiedBC0.2031.01VG50.0752.76VG40.1691.23TA0.00285.312RU0.00372.211EstimatedBC0.1721.22VG50.0633.17VG40.1321.55TA0.00010, 542.014RU0.0005, 323.013Akaikea weight (wi) is based on average AIC value across all hillslope positions.bEvidence ratio (wi/wj) is based on wi from average AIC values, where model i is the best-fitting model (BC model with modified θr), and j indexes the remaining models in the set.The BC and VG5 models represented the SWRC similarly except at the wetter end of the curve, where dθ/dψ for the BC model was discontinuous (Figure 1). The BC model had lower AIC values (Table 2), which reflected parsimony in the number of parameters for the BC model. However, the CA and BC models represent the SWRC with additional complexity due to discontinuous dθ/dψ, where ψ≥ψe.The BC and VG4 models with θr(E) fit the SWRC nearly identical to one another between −0.01 and −1.5 MPa (Figure 3). The CA model was considered among the best models, due at least in part to its parsimonious equation. However, the CA model generally underestimated water retention less than approximately −0.5 MPa for the summit and backslope positions, compared to the BC and VG models.The VG5 model with modified and estimated forms of θr would be only slightly less adequate in describing the water retention relationship compared to the highest ranking five models (Table 3). The TA and RU models consistently had relatively high evidence ratios, particularly for the measured and estimated derivations of θr, indicating low prediction ability for the SWRC compared to the other models. Nonetheless, the TA and RU models were generally adequate to represent θ(ψ) at field capacity (−0.033 MPa). The TA and RU models also estimated water retention fairly well at −1.5 MPa, but only when the measured form of θr was used.The CA and BC models may provide unrealistic description of water retention at the wet end of the SWRC wheredθ/dψ is modelled as a discontinuous function [36, 37]. The continuous VG models contain an inflection point enabling better representation of water retention near saturation [36]. However, the ability of the VG4 model to predict dθ/dψ as a continuous function did not provide superior fit compared to the CA and BC models (Table 3). Regardless, the minimum threshold water potential for cessation of germination in most weed species is less than −0.1 MPa (Table 4).Table 4
Minimum threshold water potential for germination of weed species grouped by intervals of germination sensitivity to water potential.
Sensitivity interval (MPa) Minimum threshold (MPa)Common nameScientific nameReferenceB (−0.1 to −0.5)−0.10Common cockleburXanthium strumarium L.[6]−0.20Redroot pigweedAmaranthus retroflexus L.[38]C (−0.5 to −1.0)−0.50Common groundselSenecio vulgaris L.[39]−0.50Large crabgrassDigitaria sanguinalis (L.) Scop.[40]−0.60Wild oatAvena fatua L.[41]−0.64Common lambsquartersChenopodium album L.[42]−0.64VelvetleafAbutilon theophrasti Medik.[43]−0.69Yellow foxtailSetaria glauca (L.) Beauv.[5]−0.70Green foxtailSetaria viridis (L.) Beauv.[5]−0.80GoosegrassEleusine indica (L.) Gaertn.[44]−0.80Common ragweedAmbrosia artemisiifolia L.[45]−0.83Large crabgrassDigitaria sanguinalis (L.) Scop.[5]D (−1.0 to −1.5)−1.13Common chickweedStellaria media (L.) Vill.[46]−1.15Rigid ryegrassLolium rigidum Gaudin[47]−1.20Wild oatAvena fatua L.[48]−1.21GoosegrassEleusine indica (L.) Gaertn.[5]−1.50Wild oatAvena fatua L.[49]−1.50Perennial ryegrassLolium perenne L.[50]−1.53BlackgrassAlopecurus myosuroides Huds.[51]The VG4 model was superior to the VG5 model no matter which derivation of θr  was used (Table 3). The VG4 model has one less parameter making it more parsimonious than the VG5 model. The VG4 model may also provide a superior estimation of the water retention relationships compared to the VG5 model due to independence of the n and m parameters in the VG5 model, that can lead to distinctiveness problems in the estimation process, resulting in a less accurate description of the SWRC in the dry range [52].The SWRC for the coarse soils in the summit and backslope positions were best described by the BC model, based onwi values. The CA model had the best fit to the data for the soils with greater clay content in the toeslope position (Table 2).


## 3.3.  Weed Seed Germination Sensitivity to Soil Water Potential

The measured water potential interval (−0.01 to −1.5 MPa) was divided into four subintervals to represent sensitivity groups of minimum threshold water potential for germination of weed species (Table4). No species from the literature had a minimum threshold water potential within interval A (−0.01 to −0.1 MPa). The minimum threshold water potential for germination ranged from high sensitivity (interval B) to low sensitivity (interval D).The best-fitting SWRC for each model from the variousθr derivations was retained to evaluate water potential intervals for seed germination. Models evaluated were CA with no θr parameter, and BC, VG5,  VG4, TA and RU models with modified θr.The mean water content was averaged across hillslopes because water retention was not different between hillslopes (Table5). Similar mean water contents for the water potential intervals for the summit and backslope positions occurred as a result of few differences in the underlying soil physical properties between these two hillslope positions (Table 1). The upper hillslope positions differed from the toeslope as a result of differences in soil physical properties that influence the SWRC. The summit and backslope positions contained silt loam soils, whereas the toeslopes had loam and clay loam soils with greater OM content.Table 5
Mean water content for water potential intervals based on weed seed germination sensitivity to water potential.
Mean water contentWater potential interval (MPa)Main effectFullABCD−0.01 to −1.5−0.01 to −0.1−0.1 to −0.5−0.5 to −1.0−1.0 to −1.5(cm3/cm3)HillslopeSW0.1580.2020.1320.1080.098NE0.1500.1900.1260.1020.092Hillslope positionSummit0.124b0.164b0.099b0.080b0.072bBackslope0.141b0.185b0.114b0.091b0.082bToeslope0.196a0.239a0.175a0.145a0.132aModelCA0.153c0.202a0.125d0.095c0.083cBC0.153bc0.198b0.126c0.100b0.090bVG50.153bc0.199b0.126c0.100b0.089bVG40.153b0.198b0.126c0.101b0.091bTA0.156a0.189d0.137a0.118a0.109aRU0.156a0.191c0.135b0.116a0.108aSource of variationP>FHillslope (H)0.2160.1890.3130.2430.222Position (P)<.001<.001<.001<.001<.001H×P0.3600.4140.1610.0950.078Model (M)<.001<.001<.001<.001<.001H×M0.001<.001<.0010.0030.013P×M<.0010.0100.429<.001<.001H×P×M<.0010.4310.6820.4110.309Means within a column and main effect, followed by different letters are significantly different atP≤0.05 according to Fisher’s protected LSD.The interactions for hillslope position by SWRC model are shown in Figure4. The VG4 models did not differ from the BC model for any water potential interval in either the upper or lower hillslope. The CA model overestimated water retention between −0.01 to −0.1 MPa for the upper hillslope (Figures 4(a) and 4(b)). Since interval A does not contain minimum threshold water potential for weed seed germination, deviations from the measured water retention data in interval A would not be expected to influence seed germination of any of the species listed in Table 4. Deviations in water retention by the CA model were confined to the upper hillslope.Fitted SWRC for soil water contents(θ) at soil matric potentials (ψ) for the (a) summit, (b) backslope, and (c) toeslope positions using the Campbell (CA) model with absent θr parameter, and Brooks-Corey (BC), five- and four-parameter van Genuchten (VG5 and VG4), Tani (TA), and Russo (RU) models with θr modified as one of the fitted parameters. Hillslopes are averaged. Error bars are SE for six water retention measurements. Model* represents SWRC models that have significantly different water retention within a water potential interval compared to the BC model. Water potential (ψ) intervals A: −0.01 to −0.1 MPa; B: −0.1 to −0.5 MPa; C: −0.5 to −1.0 MPa; D: −1.0 to −1.5 MPa represent groups of minimum threshold water potential sensitivity for germination.
(a)
Summit(b)
Backslope(c)
ToeslopeThe CA model underestimated water retention at a potential less than −0.5 MPa for the summit and backslope positions (Figures4(a) and 4(b)). An underestimation in water retention less than −0.5 MPa would result in an underestimation in predicted germination for weeds with minimum threshold water potential less than −0.5 MPa. However, because soil water is generally measured in the field on a content basis, and SWRC are used to translate water content to a potential basis, the CA model would predict higher water potential for a given field water content compared to the BC or VG models. This case would overestimate germination of weeds that have minimum threshold water potential less than −0.5 MPa. Thus, germination would be predicted for species that would otherwise not occur in the upper hillslope. For example, at the summit, soil with a water content of 0.07 cm3/cm3 would have water potentials of −1.08 MPa and −0.57 MPa predicted by the BC and CA models, respectively. The BC model would predict germination of weeds only for interval D in Table 4, whereas the CA model would predict germination of weeds for interval D and almost all of interval C.The mean water content of the TA and RU models differed from that of the BC model across all water potential intervals for both upper and lower hillslope positions (Figure4). Water retention was underestimated by the TA and RU models for the −0.01 to −0.1 MPa interval, which would not be expected to influence modelling of germination for weeds listed in Table 4. However, overestimation of water retention by the TA and RU models for soils drier than −0.1 MPa would result in overestimation of seed germination at specific potentials. On a field water content basis, this would result in underestimation of seed germination by all species listed in Table 4. For instance, at the toeslope position, a soil water content of 0.15 cm3/cm3 would have a predicted water potential of −0.45 MPa with the BC model and −0.94 MPa with the TA model (Figure 4(c)). The BC model would, therefore, predict germination of weeds for intervals C and D, whereas the TA model would predict germination of weeds for only interval D. Predictions of weed seed germination from a range of minimum threshold water potentials are clearly dependent on the fit of the SWRC model.


## 4.  Conclusion

Models of the SWRC with derivations ofθr were developed to describe the water retention properties of the shallow seedling recruitment zone. Soils along the hillslopes in this study had diverse underlying soil physical properties that influence the SWRC. Soil topographical relationships for the physical properties that affect the SWRC allowed hillslope positions to be used as categorical units for describing the SWRC.Evaluation of analytical models to describe the SWRC of the shallow seedling recruitment zone indicated that the BC, VG4, and CA models with modified and estimated forms of θr were similar but superior to models with θr designated at −1.5 MPa. However, when considering intervals of water potential, the CA model underestimated water retention for the upper hillslope at the dry end of the SWRC in the range for minimum threshold water potential for seed germination. Setting the residual water content at zero as in the case of the CA model will inevitably underestimate water content. The TA and RU models estimated water retention reasonably well at field capacity, but tended to overestimate water retention for the dry end of the SWRC.Weed seed germination is generally predicted by models using hydrothermal properties of the soil environment, of which water potential is a component. Because the best-fitting SWRC is needed for a dynamic model to predict weed germination within the range of minimum threshold water potentials for germination, water retention models that specifically deal with the SWRC at the dry end of the curve need further investigation. Future efforts could include analysis of the SWRC below −1.5 MPa to include species with very low minimum threshold water potentials for germination. As well, future research could evaluate the SWRC for a range of soil textures beyond that examined in the current study.The evaluatedθr derivations showed that all models were best represented by specifying θr as a modified or estimated parameter. For modelling weed seed germination, it is most important to adequately describe the drier end of the SWRC below −0.1 MPa. Accurate representation of the SWRC over the range of minimum threshold water potential for seed germination by modifying or estimating the θr parameter is necessary for modelling seed germination by avoiding deviations in estimated water potential and predicted seed germination.



---

*Source: 812561-2012-05-06.xml*