|
|
||||||||
a USDA-ARS, Alternate Crops and Systems Lab., Bldg. 007, Rm. 116, BARC-West, 10300 Baltimore Ave., Beltsville, MD, 20705
b USDA-ARS Hydrology and Remote Sensing Lab., Bldg. 007, Rm. 106, BARC-W, 10300 Baltimore Ave., Beltsville, MD 20705
c College of Agricultural Sciences, P.O. Box 21360, Rio Piedras, PR 00928
d Dept. of Soil, Crop, and Atmospheric Sciences, Bradfield Hall, Cornell Univ., Ithaca, NY 14853
* Corresponding author (dtimlin{at}asrr.arsusda.gov)
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
An objective of site specific management is to apply appropriate levels of agricultural inputs to localized areas. Since areas with deeper soils that facilitate extensive root growth tend to be more productive, it would be useful to have a tool to estimate or quantify the potential productivity due to water availability related to soil depth. An estimate of the spatial distribution of potential soil productivity can be obtained from knowledge of the spatial variability of soil depth. The relationship between soil properties and fertilizer requirements should be considered for site specific management (Cahn et al., 1994). Fertilizer or seed applications could therefore be tailored for the potential response of the crop as a function of soil depth and water availability (Mathews et al., 1997; Barnhisel et al., 1996).
Crop simulation models have become a useful tool to characterize and quantify yield and available water. Paz et al. (1998), using a soybean model, showed that yield variability correlated with variability of simulated water stress. Rooting depth and soil water holding capacity were important variables. Soil depth was an important parameter in a productivity index model used by Khakural et al. (1996) to estimate the spatial variability of crop yields. Moore and Tyndale-Brisco (1999) using crop models showed that much of the variability of wheat response to nitrogen could be explained by differing soil water holding capacities. Mathews and Crosser (1997) used CERES-wheat to study variable rate nitrogen application. They reported that the optimal nitrogen distribution was to apply more fertilizer to deeper soils and less to shallow soils.
Many of the models currently applied in precision agriculture have complex input requirements and may be more detailed than necessary for certain applications. They also require some form of calibration. Timlin et al. (1986) described a simple water budget model that calculated corn grain yields in shallow soils. The model was based on the assumption that a major portion of the yield variability in an otherwise uniformly managed field, is due to variability in soil available water. This in turn would be a function of depth to some root restricting layer, pore size distribution, and weather. The model only requires knowledge of soil water holding capacity, rooting depth, crop growth stages, and weather data.
The purpose of this study was to gain further understanding of yield and available water relationships under varying seasonal weather conditions and soil rooting depths. The objectives were (i) to determine if the effect of soil depth on yield is mainly due to reduced water availability and (ii) to evaluate the use a simple water budget to describe corn grain yield variability due to variation in soil depth.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Soil depth, defined as depth to bedrock, was measured with a hand auger. Soil water characteristics were measured at selected sites on cores or clods using a pressure plate apparatus. Water content at wilting point was taken to be that at -1.5 MPa using disturbed samples. The matric potential corresponding to field capacity was the tensiometer reading at 2 d after a saturating rainfall, when evapotranspiration was minimal.
The Aurora site had been in continuous corn. The field was moldboard plowed in the spring to 0.15 m and treated before planting with atrazine [6-chloro-N-ethyl-N'-(1-methylethyl)-1,3,5-triazine-2,4-diamine] to control weeds. The field was fertilized with 30 kg ha-1 nitrogen, 60 kg ha-1 phosphorus, and 30 kg ha-1 potassium. An early season variety of corn was planted in 0.76-m rows and nitrogen was sidedressed at the rate of 70 kg ha-1 N as anhydrous ammonia.
Yields were recorded over three growing seasons. In 1983 and 1984, yields were harvested at the nodes of a grid; yields from a total of 70 plots were recorded. The grid cells measured 20 by 17 m. Each plot was two rows wide by 5.3 m long and was located at the upper left hand corner of the grid cell, except in 1983 when four rows per plot were harvested. In 1985, the grid locations were not harvested due to weed problems. Instead, five transects parallel to the slope and containing 13 to 16 plots each were harvested to determine yield from areas where weeds had been controlled.
Nine pairs of irrigated and non irrigated plots were set out along transects. The irrigation experiment was carried out over one growing season, during the summer of 1985. Each plot included four corn rows 6.1 m long and the pairs of plots were separated by distances of approximately 10 to 20 m. The irrigated and non irrigated plots were instrumented with tensiometers installed at 0.15 m depth intervals. Gypsum blocks were installed in the non irrigated plots to cover matric potential ranges more negative than about -80 to -100 kPa. A pressure transducer as described by Marthaler et al. (1983) was used to measure the air pressure inside the partially water filled tensiometers. Irrigation was carried out when the soil water tension, as measured by the tensiometers in the 0 to 0.3-m layer, reached -50 kPa. Irrigation in each plot was accomplished by connecting four parallel drip irrigation lines to two 264-L barrels. Each drip irrigation line was no more than 15 cm from the crop row. The inner two rows of each four-row plot were harvested for yield determination.
Weather data were collected at Aurora from a weather station at the farm office located about 1 km from the field. The weather data include daily max-min temperatures, precipitation, and class A pan evaporation.
Sixteen years of meteorological, and corn grain yield and crop phenological data also were obtained from the crop-climate experiments at Canton Agricultural and Technical College at Canton, NY (R.F. Lucey, unpublished data, 1984). Recommended management practices and a short season hybrid of corn (Funk G-11-A) were used. Soil water holding capacities were also available. These data were used for further testing and validation of the model.
Model Background
Corn grain yield predictions in the Soil Depth-Soil Productivity Model (SD-SPM) are based on the relationship between relative yields and relative transpiration rates (Hanks, 1974):
![]() |
Yields are estimated from a water budgeting procedure that calculates the relative transpiration ratio.
The SD-SPM requires weather, soil, and crop data. Weather data include daily evapotranspiration and rainfall throughout the growing season. Soil data required by the model are depth to a root restricting layer, soil water holding capacities, and the water content at wilting point of the soil horizons. As a first approximation, the 1.5 MPa water content is used for the wilting point water content. Most of the water held in soil pores at tensions greater than 1.5 MPa can be potentially lost through evaporation and is not readily available to field crops such as Z. mays. These soil data are generally available from standard soil characterization.
Plant data that the model requires include the maximum potential rooting depth, the depth at which the rooting density is one half that at the surface for a fully developed root system grown to the potential depth, and the distribution of corn growth stages. A 1-m depth was used as the maximum potential rooting depth for corn hybrids specific to the soils and climate of New York (Olson, 1983). Layers with a bulk density >1.8 Mg m-3 were considered restricting to root growth (Olson, 1983). In non restricting subsoils, the depth at which the root length density is one half that of the surface can be taken to be one half the potential rooting depth (i.e., 0.5 m) (Mengel and Barber, 1974). If there are soil horizons that have a bulk density >1.50 Mg m-3 with few roots at some depth <0.5 m, then the depth of the lower boundary of this horizon can be used. The limiting value of 1.5 Mg m-3 was determined from an analysis of data from Fehrenbacher and Rust (1956). The depth at which root length density is one half the root density at the soil surface was estimated from soil characterization data and profile descriptions.
Crop growth stages are input to the model. Corn growth stages used in the model are vegetative, late vegetative, silking, blister kernel, and maturity (Table 1). If crop growth stages are not known, the calendar dates when the crop reaches late vegetative and silking stages can be estimated from cumulative heat units since planting. The number of heat units required for the corn crop to reach a given growth stage can be obtained by comparing accumulated heat units and measured crop phenological data over a period of several years (Daynard, 1972).
|
In SD-SPM, a root growth function adapted from Davidson et al. (1978) was used to calculate relative root density. The root growth function was configured to calculate a root distribution that, at maturity, extends to a specified lower boundary and has one half the surface root density at a specified depth (both given as input). The purpose of this function was to obtain a reasonable relative root distribution with depth. The maximum potential corn rooting depth and the depth at which root length density is one half that at the surface for a mature root system are input parameters. The maximum potential rooting depth for a particular soil is defined as the maximum depth roots would reach in a deep well drained soil similar to that soil type. It is assumed that maximum potential rooting depth is reached at the time the crop has fully tasselled. Lateral root growth may continue until 75% of the crop has silked.
Calculated root density in the model is assumed to be independent of the presence of a root restricting layer below it. Where soil conditions restrict root growth, the potential root system is simply truncated at the depth of the restricting layer. While this may not be accurate, the purpose of the root growth routine is to only to provide an approximate measure of root distribution to model relative water uptake from different layers given simple and easily obtainable input data.
The root length density, R [cm root (cm depth)-1], at any time, t, and depth, z, is defined as:
![]() | (1) |
0) at time t,
is a term that gives the root distribution function its shape (cm-2), Z is depth, and L(t) is the vertical extent of the roots at time t. Time, t, is defined as days after planting. The shape factor,
, the depth of the roots, L, the maximum root length density near the surface, RMAX, and the depth, Z1/2, at which the root length density is one half RMAX, at time t are modeled as quadratic functions in time.
, RMAX, L, and Z1/2, are:
![]() | (2) |
![]() | (3) |
![]() | (4) |
![]() | (5) |
There are several fixed conditions: (i) the depth of the roots at 29 d (0.47 m here), (ii) the potential depth of the roots at the end of the tasseling stage (given as input), (iii) the depth at which root length density is one half that at the surface (Z1/2) at the time of 75% silking (given as input), (iv)Z1/2 is zero until 5 d after planting, and (v) the rates of change of a) vertical root growth and b) Z1/2 with time are both zero at the time of 75% silking. These given conditions provide six known values of the variables that allow us to solve for the six coefficients (ß1-ß6) in Eqs. [4] and [5]. The conditions given in 1 to 5 above were developed from data given in NaNagara et al. (1976).
In SD-SPM, water from precipitation is added to the soil profile as piston flow displacement. Actual transpiration (Ta) and potential transpiration (Tp) are estimated by a water budgeting procedure based on precipitation, pan evaporation, crop growth stage, root growth, and the water holding capacity of the soil within the rooting zone. The ratio of actual to potential transpiration (Ta/Tp) is estimated as a function of the available water in the rooting zone (Fig. 1). Available water is a function of climate and the soil properties that determine water storage capacity. Water uptake from a particular layer depends on the relative root density in that layer (Selirio and Brown, 1979).
|
![]() | (6) |
![]() |
Hiler and Clark (1971) and Shaw (1974) used this stress index to calculate corn grain yield as:
![]() | (7) |
In order to be able to represent results from different sites on a more general basis, we will redefine Eq. [2] as a relative yield equation by dividing Eq. [7] by potential yield (YP). The result is:
![]() | (8) |
Relative yields and a relative stress index were also calculated for the Canton data which were originally used to test the water stress component of the model (Timlin et al., 1986).The authors tested the water stress component by using half the Canton data (8 yr) to calibrate the model and the other half to test the predictions. The equation fit to all the data was:
![]() | (9) |
Where yield was in kg ha-1 and OCHU are Ontario Corn Heat Units (Daynard, 1972). Note that in this equation, Yp from Eq. [7] was calculated as a function of Ontario Corn Heat Units (7.51 OCHU-4441) and the water stress-response coefficient, A from Eq. [7] is 198 kg ha-1. A calibration equation developed from half the data predicted the validation grain yields well.
Potential yield calculation was separated from the water stress component in order to calculate a relative water stress response coefficient (AR in Eq. [8]) for the Canton Data. Since the rooting depth of the soil at this site was relatively deep (>0.90 m), the variation in potential yield from year to year was assumed to be a function of temperature alone and potential yield was estimated from Ontario Corn Heat Units (OCHU). In New York state, temperature and precipitation are major sources of variation in corn grain yields (Russo, 1978). Estimates of potential yields (Yp) for 16 years of data were obtained by regressing measured yields on OCHU only (setting SS to zero). A relative water stress response coefficient for the Canton data was obtained by regressing relative yields (Y/Yp) on calculated water stress indices.
Potential yield (Yp) was obtained by three methods for the Canton and Aurora data sets. These included mean irrigated yield (Aurora), OCHU in Eq. [9] (Canton), and by regressing yields measured over a range in soil depths on calculated seasonal stress indices (Aurora). The use of mean irrigated yield provides an estimate of potential yield that is independent of the data used to calculated stress indices and independent of the reduction in yield due to water stress, slope A in Eq. [8]. Potential yield calculated from OCHU can be shown to be statistically independent of the reduction of yield due to water stress.
We recognize that there may be other sources of limiting factors for crop yields and it may not be possible to completely eliminate water stress. The proposed relationship, however should be able to quantify relative differences in yields when the main source of yield variability is due to water stress.
Statistical Analysis
Correlations, regressions, and Students t-tests were carried out using SAS statistical software (SAS Institute, 1995). Analysis of the irrigation treatment effects was done by testing the significance of the relationship between the increase in grain yield due to irrigation and soil depth. The increase in yield due to irrigation was regressed on soil depth and the significance of the relationship was tested by comparing the slope to 0. The water stress response coefficients for the relative yield equations at the different sites and years were compared using Proc Reg of SAS (SAS Institute, 1995). We used indicator variables (Neter and Wasserman, 1974, p. 279338) in multiple linear regression to obtain parameters for all relative yield equations to calculate values of the Students t-statistic to make comparisons.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
Weather
The weather conditions over the three year period varied from very wet and cool to dry and warm (Fig. 2). The summers were warm and dry in 1983 and 1985. In contrast, the summer of 1984 was cool and wet. The water potential data (Fig. 3) reflect the dry conditions for 1985. The soil in the shallow root zone (<0.5 m) reached lower matric potentials than the soil with the deeper root zone and remained at these lower potentials for a long period of time.
|
|
|
|
|
|
|
Other sources of variability were most likely factors associated with lack of uniformity in management. Weeds were difficult to control and had varying effects on yield. However, in the irrigated and non irrigated plots, weeds were controlled by hand weeding. Plant populations and fertilizer applications appeared to be somewhat variable. Nevertheless, in spite of these sources of variation, the calculated stress indices do follow the trends of the data.
Relationships Between Relative Yields and Relative Stress Indexes
When Eq. [7] was converted to the relative yield basis (Eq. [8]) for the Aurora data by dividing measured yields by the calculated YP values (Table 2) and regressing on the stress indices, the resulting relative water stress response coefficients, AR, were similar for both 1983 and 1985 (Table 2).
A relative water stress response coefficient was also calculated for the Canton data. Potential yields were calculated from the OCHU by setting the seasonal stress index (SS) to zero (Eq. [9]). There was no significant interaction at this site between the transpiration ratio and temperature (as determined by regressing yield on both calculated stress index and temperature and comparing the coefficients with those obtained from single variable regressions). The relative stress index was 0.029 and was similar to those calculated for Aurora (Table 2).
The relative grain yield vs calculated seasonal stress index for the 1983 and 1985 Aurora data, and the Canton data were linear and the three data sets appear to conform to the same relationship (Fig. 8). The relative stress response coefficients were not significantly different from each other. The F values to test the differences between slopes were close to 0 for all the comparisons. Therefore the three data sets were combined and one regression equation fit to them. The fitted equation is:
![]() | (10) |
|
![]() | (11) |
The fitted equation describes the relative yield-seasonal water stress index data well (Fig. 8).
Comparison of Measured and Estimated Relative Yields Over a Range of Soil Depths
Independent estimates of relative yields for 1985 were calculated by dividing measured yields by mean irrigated yield. This gave an estimate of relative yield that was independent of the regression with stress indices. In order to have an estimate of the relative stress index that was independent of the 1985 data for comparative purposes, the slope in Eq. [8] was taken to be 0.029, that calculated from the Canton data. The comparison, shown in Fig. 9 demonstrate that the model reasonably predicts the average changes in relative yield as rooting depth varies. The mean square error was 0.017 in units of relative yield.
|
| SUMMARY AND CONCLUSIONS |
|---|
|
|
|---|
Where soil depthyield relationships existed, the measured data followed the same trend as the predicted relationships. One relative yield equation was found to predict all trends with reasonable accuracy. The equation was:
![]() |
The objective of this research was to show that the major source of variability in corn grain yields could be related to water holding capacity, a function of soil depth. Furthermore, a simple water budget model could be used to quantify this variability. Many factors other than water were probably involved in the yield variability. However, observed trends of yield with soil depth were predicted reasonably well by the model on the basis of water relations, and the irrigation experiments verified this. The water stress response coefficient for the water stress index to predict relative yield was shown to only vary slightly for three data sets from New York using a short season variety of corn. The generality of this expression needs to be investigated over a wider range of data. An independent test of the model where potential yield was obtained from irrigation experiments and the relative stress index component from the Canton data, showed that the model predicted relative yields for the 1985 Aurora data reasonably well (mean square error of 0.017). The predicted values were closest to measured, where the soil depth was less than 0.60 m. The model tended to underpredict relative yield where soil depth was greater than 0.60 m.
The simple water budget and relative yield model as described here can be used as a tool to evaluate potential productivity in a landscape as a function of available water and weather. The model can be used to quantify the effects of weather on corn yield.
| ACKNOWLEDGMENTS |
|---|
Received for publication August 17, 1999.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
P. Jiang, N. R. Kitchen, S. H. Anderson, E. J. Sadler, and K. A. Sudduth Estimating Plant-Available Water Using the Simple Inverse Yield Model for Claypan Landscapes Agron. J., May 7, 2008; 100(3): 830 - 836. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Jiang, S. H. Anderson, N. R. Kitchen, K. A. Sudduth, and E. J. Sadler Estimating Plant-Available Water Capacity for Claypan Landscapes Using Apparent Electrical Conductivity Soil Sci. Soc. Am. J., October 29, 2007; 71(6): 1902 - 1908. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. L. S. Morgan, J. M. Norman, and B. Lowery Estimating Plant-Available Water Across a Field with an Inverse Yield Model Soil Sci. Soc. Am. J., March 1, 2003; 67(2): 620 - 629. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. L. Feng, A. Meiri, and J. Letey Evaluation of a Model for Irrigation Management Under Saline Conditions: II. Salt Distribution and Rooting Pattern Effects Soil Sci. Soc. Am. J., January 1, 2003; 67(1): 77 - 80. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Crop Science | |||
| Journal of Natural Resources and Life Sciences Education |
Vadose Zone Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||