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.

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 () 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 becomes zero and unsaturated hydraulic conductivity ceases [15]. Furthermore, some SWRC models represent near the wet end of the range of soil water as a continuous function, whereas other SWRC models represent 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: where is volumetric water content, is matric potential, is saturated volumetric water content, is residual soil water content when is infinitely small, is a curve-fitting parameter, known as the air entry matric suction, and is a pore-size distribution factor. The restriction is applied for 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: where is the air entry matric suction, and is a curve-fitting parameter representing pore-size distribution where , and is an empirically derived constant. The Campbell equation is similar to the Brooks-Corey model but with . The van Genuchten (VG5) model [15] is an expression originally based on five parameters: where is a curve-fitting parameter related to and the slope of the curve at inflection, and and 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 [15]. The Tani (TA) model [24] is a three-parameter model using the expression: where 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: where is related to the width of the pore-size distribution and is interpreted as the inverse of the air entry matric suction, and 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 were taken from measurements at saturation. Three derivations of were explored to determine the influence of 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 derived from measurements at −1.5 MPa (the permanent wilting point) were designated . For practical purposes, 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 was modified by a logarithmic equation describing the adsorption of water on soil in the dry range of the retention curve [27, 28]: where is a curve-fitting parameter representing the volumetric water content when , and is the matric suction at oven dryness, which is generally accepted to be 103 MPa. The modified form of is assigned . The third derivation of was estimated as one of the fitted parameters, referred to , which contained the restriction . In all cases, the optimization results were improved by setting the measured value of as a constant parameter. Previous SWRC estimates have shown improvement by using known constant rather than optimized values of [29]. The BC and CA models were analysed as discontinuous models, whereas the VG, TA, and RU models were analysed as continuous 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: where 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 and Akaike weights [31, 33]. Delta AIC is a measure of the AIC differences of each model relative to the best model derived by: where is the AIC value for model , and 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 [31]. Akaike weights provide an effective way to scale and interpret the values by comparing the ratio of each model to the best model relative to the entire set of candidate models by: given a set of models being evaluated. Akaike weights compare models on a scale of one (the sum of equals one) indicating the weight of evidence that model is superior among the set of 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 , where is the estimated best model in the set, and 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 and models [31].

2.5. Weed Seed Germination Sensitivity to Soil Water Potential

The best-fitting curve according to Akaike weights from the 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]: where is the mean water content (cm3/cm3), and is the SWRC fitted to the measured water contents (cm3/cm3). The integration boundaries and define the integration interval over which the mean water content was derived. The log-normal distribution of was calculated using to avoid overweighting of the dry end of the SWRC.

Mean water contents for the best-fitting 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 (Table 1). 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.

3.2. Soil Water Retention Model Evaluation

The respective SWRC models (excluding CA) all exhibited better fit with modified or estimated forms of over that of measured as indicated by lower values of AIC (Table 2). This may be an indication that modified or estimated forms of were better represented by values other than those measured at −1.5 MPa. Estimating or modifying as a free parameter provided better fit particularly in the dry range of the curve compared to values (Figures 13). Models containing the modified form of 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 was modified rather than estimated. The TA and RU models, however, had considerably improved fit with compared to mainly in the dry range of the curve (Figures 2 and 3).

The 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 , which was the best-fitting model based on derived from average AIC values (Table 3). An evidence ratio of 3 or less in relation to another model provides little evidence that model is superior [31]. Accordingly, BC with and , VG4 with and , and CA models would serve nearly equally well in approximating the water retention information.

The BC and VG5 models represented the SWRC similarly except at the wetter end of the curve, where 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 , where .

The BC and VG4 models with 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 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 , 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 was used.

The CA and BC models may provide unrealistic description of water retention at the wet end of the SWRC where 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 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).

The VG4 model was superior to the VG5 model no matter which derivation of   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 and 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 on 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 (Table 4). 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 derivations was retained to evaluate water potential intervals for seed germination. Models evaluated were CA with no parameter, and BC, VG5, VG4, TA and RU models with modified .

The mean water content was averaged across hillslopes because water retention was not different between hillslopes (Table 5). 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.

The interactions for hillslope position by SWRC model are shown in Figure 4. 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.

The CA model underestimated water retention at a potential less than −0.5 MPa for the summit and backslope positions (Figures 4(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 (Figure 4). 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 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 were similar but superior to models with 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 derivations showed that all models were best represented by specifying 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 parameter is necessary for modelling seed germination by avoiding deviations in estimated water potential and predicted seed germination.

Acknowledgments

This paper was funded by the Manitoba Rural Adaptation Council, the Natural Sciences and Engineering Research Council of Canada, and the Canadian Wheat Board.