|
|
||||||||
a USDA-ARS, Bldg. 3702, Curtin Rd., University Park, PA 16802
b Division of Plant Science, Univ. of Missouri, Columbia, MO 65211
c Dep. of Crop and Soil Sciences, Pennsylvania State Univ., 116 ASI Bldg., University Park, PA 16802
* Corresponding author (john.schmidt{at}ars.usda.gov)
Received for publication June 27, 2006.
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: EONR, economic optimum nitrogen rate
| INTRODUCTION |
|---|
|
|
|---|
Corn yield variability within a field or between fields has been well documented with precision agriculture technologies (i.e., yield maps), and yield ranging from <1.0 Mg ha1 to >12 Mg ha1 within one field during one growing season is common (Taylor et al., 2001). However, yield variability does not automatically translate into variability in EONR (Vanotti and Bundy, 1994), so developing site-specific N fertilizer recommendations should not automatically follow from a variable yield map. Schmidt et al. (2002) illustrated that irrigated corn yield response functions to N fertilizer for several within-field locations were the same (i.e., maximum yield was obtained with the same N rate), although maximum yield for these same locations ranged from 6.4 to 10.6 Mg ha1 during one growing season. In another irrigated corn field from the same study, maximum yield for several within-field locations was obtained with N rates ranging from 56 to 182 kg N ha1. Fox and Piekielek (1995) indicated that there was no relationship between maximum yield and EONR (r2 = 0.08) in their evaluation of 57 site-years in Pennsylvania between 1982 and 1994, despite maximum yields ranging from 6.7 to 12.4 Mg ha1 and EONR ranging from 67 to 212 kg N ha1. Successful site-specific N management for corn depends on determining the spatial distribution of EONR across a field or between fields within the geographic region of a farm.
Field and sub-fieldscale variability in EONR for corn has been documented recently (Mamo et al., 2003; Scharf et al., 2005), but soil characteristics or other causal factors that could be used to develop a spatial EONR map have not been identified. Katsvairo et al. (2003) concluded that site-specific N management in the Northeast requires more spatial information than is provided by a late spring soil NO3N test and/or yield maps from previous growing seasons. Field-to-field (within and across years) variability in EONR ranged from 22 to 203 kg N ha1 for a continuous corn rotation at 11 site-years in Pennsylvania (Fox and Piekielek, 1983). In this same study, EONR varied from 0 to 215 kg N ha1 across all site-years, including various crop rotations and histories of manure application. Despite considerable research in Pennsylvania that demonstrates field-to-field (within and across years) variability in EONR (Fox and Piekielek, 1983, 1995), there has been little or no research that explores whether there might be sufficient within-field spatial variability in EONR to justify site-specific N applications in the Northeast USA.
The hillslope is a typical agricultural landscape unit in the Northeast, with potentially EONR-dependent soil characteristic variability, such as soil water content. Previous research has suggested that water availability may affect EONR. For example, Fox and Piekielek (1998) noted that maximum grain yield in corn was linearly related (r2 = 0.69) to July precipitation for 15 yr of results from a N fertilizer rate study at Rock Springs, PA. In another unpublished report, Fox and Piekielek demonstrated that maximum corn yield was linearly related to July rainfall for rainfall less than 9.4 cm, then reached a maximum yield with a linear-plateau relationship (r2 = 0.62) for 20 yr of results from the Rock Springs farm (the latter report presumably includes data from the 1998 publication). Although not presented by Fox and Piekielek (unpublished report, 2001), EONR from their study was also linearly related (r2 = 0.5; P > F = 0.001) to July rainfall but did not result in a plateau-limiting response. In this example, July rainfall was a simple indicator of water availability during a period of rapid vegetative growth (usually 8-leaf growth stage to tasseling) and high water demand by corn; however, spatial variability in soil water availability may have provided a better indicator for yield or EONR.
Variability in soil water content along a hillslope is not simply a function of elevation and rainfall, but, as Famiglietti et al. (1998) demonstrated for a hillslope near Austin, TX, surface soil water content (05 cm) depends on soil porosity and hydraulic conductivity during wet soil conditions and on relative elevation, aspect, and clay content during dry conditions. Ridolfi et al. (2003) underscored the complexity of soil moisture dynamics along a hillslope by identifying 10 different phenomena that contribute to soil moisture spatial variability. Pachepsky et al. (2001) used topographic features to explain variability in soil water content and suggested that topographic variability had a potential use for interpreting field-scale variability in precision agriculture.
The objective of this study was to characterize the spatial variability in EONR for corn along a 300-m hillslope, considering the impact of soil water availability on EONR and the potential for using this landscape and soil characteristic for site-specific N management in the northeast USA.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Before applying N treatments, soil samples were collected for routine soil analyses at Locations 2, 6, and 10, representing the toe slope, mid slope, and head slope positions, respectively (Fig. 1a). Soil samples (0- to 15-, 15- to 30-, and 30- to 60-cm depths) consisted of three subsamples collected with an open-faced bucket auger (5-cm i.d.) and composited for each depth. Analyses were completed at The Pennsylvania State University Agricultural Analytical Service Laboratory, except for inorganic N, which was determined by flow injection analysis of 2 M KCl extracts (QuikChem Methods; Lachat Instruments, Loveland, CO). Results are summarized in Table 1. Surface soil pH ranged from 5.3 to 5.7, which represented slightly less-than-optimum to optimum conditions for growing corn. Surface soil Mehlich 3 soil test P (Wolf and Beegle, 1995) was slightly less than optimum (3050 mg kg1) at Location 10 (24 mg kg1) and optimum to slightly greater than optimum at Locations 2 and 6 (46 and 77 mg kg1, respectively). Mehlich 3 soil test K (Wolf and Beegle, 1995) (015 cm) for these soils was slightly less than optimum (100150 mg kg1) to slightly greater than optimum (Table 1). Surface soil organic matter content by loss on ignition (Schulte, 1995) was similar among sites, ranging from 24.3 to 28.7 g kg1 for the highest to lowest landscape positions represented with these three locations. Inorganic soil N (015 cm) was considerably greater at the lowest landscape position at 66.2 mg kg1, compared with 14.7 and 12.1 mg kg1 for the higher landscape positions. General soil nutrient characteristics, except for inorganic N, were similar among soils along this toposequence, representing slightly less than to slightly greater than optimum categories, and represented typical growing conditions in central PA.
|
Corn grain was harvested from three of the four inside rows with a combine modified for plot work. Corn grain yield was adjusted to a moisture content of 155 g kg1. We fitted a quadratic, linear-plateau, quadratic-plateau, and exponential function to evaluate yield response to N treatment at each location. We chose these functions based on the literature (Cerrato and Blackmer, 1990; Schmidt et al., 2002) and the observed shapes of the scatterplots of yield versus N treatment. We assessed the goodness-of-fit of these functions based on the significance of an F test and the magnitude, randomness, and normality of the model-fit residuals. An ideal model would have the smallest residuals that exhibit a random pattern and are normally distributed with significant treatment effects. The EONR at each location was determined for selected yield response functions using a N fertilizer cost of $0.66 kg1 ($0.3 lb1) and a corn price of $0.078 kg1 ($2 bu1), equating the first derivative of the response equation to the fertilizer/corn price ratio and solving for X (Cerrato and Blackmer, 1990).
Nitrogen treatment effects on grain yield were determined with PROC GLM (SAS Institute, 1999). Regression analyses were completed using PROC REG (SAS Institute, 1999) for linear and quadratic functions and PROC NLIN (SAS Institute, 1999) for exponential, linear-plateau and quadratic-plateau functions.
Volumetric soil water content was determined using a factory-calibrated time domain reflectometry moisture meter (TRIME-FM3) with a cylindrical probe (T3 probe) (both from Imko GmbH, Ettlingen, Germany). The TRIME-FM3 provides an effective way to obtain profile soil water content at multiple landscape positions, and the performance of this instrument has been evaluated by Laurent et al. (2005) (RMSE = 0.0662 when compared with neutron probe measurements). Soil water access tubes were installed at Locations 1 through 10 located in the same row (between the first and second plot along the length of the hillslope) adjacent to and immediately in front (within 5 m on the downhill side) of each location (Fig. 1). A hydraulic soil probe was used to remove a soil core to 1.1-m depth. One PVC tube (5.0-cm i.d.) was fitted snugly into each hole. Access tubes were not placed in the plots receiving varying N fertilizer rates but were placed in "alley" areas that received a uniform 200 kg N ha1 immediately after planting. We sampled soil water content inside the PVC access tubes at 0- to 20-, 10- to 30-, 30- to 50-, 50- to 70-, and 70- to 90-cm depths at approximately weekly intervals between 5 June and 2 Sept. 2005 and after significant rain events. This period corresponded to approximately the 5-leaf growth stage to grain fill. Equivalent depth of soil water, Wp, was calculated for the top 90 cm of the soil profile using Eq. [1]:
![]() | [1] |
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
The quadratic-plateau model provided the smallest residual sum of squares (SS) for 6 of 10 locations (Locations 1, 2, 4, 6, 8, and 9) (Table 2). At Locations 5, 7, and 10, the exponential model provided the smallest residual SS; however, the quadratic-plateau model performed almost as well at these locations, with the residual SS within 5% of those observed for the exponential model. None of the regression models considered at Location 3 was significant. The residual SS for every model at this location was more than twice as great as observed for any response functions at any other location (Table 2). Mean grain yield response to N fertilizer across all 10 locations was best described with the quadratic-plateau function. The quadratic-plateau model is often selected in the literature to describe corn yield response to N fertilizer, especially when summarizing data from multiple fields or years (Derby et al., 2005) or data from strip plots representing the length of a field (Scharf et al., 2005). Selecting the quadratic-plateau model is also consistent with the model selected by Cerrato and Blackmer (1990) in describing corn yield response to N fertilizer at 12 site-years in Iowa. Fox and Piekielek (1995) used the quadratic-plateau model to describe corn yield response to N fertilizer at 57 site-years in Pennsylvania. Because the quadratic-plateau model most often provided the smallest residual SS across all locations, provided the smallest residual SS for the mean grain yield response for this field, and is consistent with examples cited from the literature, this model was selected as the best model from which to proceed in this study in determining EONR. However, because the exponential model provided overall conclusions similar to the quadratic-plateau model and some modelers prefer continuously differentiable response functions, parameter coefficients are also provided for the exponential model.
Parameter coefficients and model statistics for the quadratic-plateau and exponential response functions are provided in Table 3, and the responses are depicted in Fig. 2 . The primary difference (besides residual SS) between selecting one of these models in favor of the other is that EONR for the exponential model was generally slightly greater than EONR for the quadratic-response function (Fig. 2), but conclusions based on results from either model were the same.
|
|
If the EONR based on the mean yield response (137 kg N ha1) were selected as the most appropriate N application on this hillslope, 44% of the field would have received N application within 20 kg N ha1 of the observed EONR. The EONR for another 44% of the field deviated from field-mean EONR by 40 to 50 kg N ha1 (Table 3). At Location 8, representing 11% of the field, observed EONR was 90 kg N ha1 less than field-mean EONR. These results are similar to those observed by Scharf et al. (2005), who indicated that EONR for more than 50% of the field deviated from the median EONR by at least 34 kg N ha1. Variability in EONR along this hillslope implicates the potential for improving N management in this field and in the overall Chesapeake Bay Watershed through site-specific technologies.
The objective of N recommendations developed for corn should be to have producers apply N fertilizer at the EONR, avoiding the economic risk associated with less N and the environmental risks associated with more N. Nitrogen recommendations for corn in many states, which have been designed to be implemented on a field-to-field basis and for large geographic regions, are a linear function of yield (or yield goal). Some examples in this category include Pennsylvania (Beegle, 2004), Colorado (Mortvedt et al., 1996), and Nebraska (Shapiro et al., 2003); whereas some states have developed N recommendations that do not include yield goal (Iowa State University Extension, 1997) or place less emphasis on yield goal (e.g., Minnesota) (Randall et al., 2003). Whether any of these N recommendations designed for large geographic areas can be successfully implemented on a site-specific basis is undetermined; consequently, the type of research presented here is important to improving N recommendations and improving the scientific approach to developing N recommendations based on additional information accessible to producers through new technologies.
Observed grain yield at EONR in this study ranged from 11.1 Mg ha1 at Location 10 (the highest position in this toposequence) to 13.5 Mg ha1 at Location 1 (the toe slope position) (Table 3). The relationship between grain yield at EONR (Mg ha1) and local relief (m) (omitting Location 3) along this toposequence suggests that local relief could be used as an indicator of yield potential (yield = 13.30.167 x local relief; r2 = 0.59; P > F = 0.02); however, using local relief as an indicator for varying N applications was not supported by the relationship between EONR and local relief. Although greater yield was obtained on lower positions in this landscape, EONR was unrelated to local relief (r2 = 0.14; P > F = 0.31). These two relationships seem to suggest that EONR was not related to grain yield, and the relationship between EONR and grain yield at EONR (r2 = 0.43; P > F = 0.18) (Fig. 3 ) only marginally supported the concept of making N recommendations based on yield goal. Based on this relationship (Fig. 3), EONR increased from 68 to 141 kg N ha1 as grain yield observed at EONR increased from 11.1 to 13.1 Mg ha1 (Fig. 3). This corresponds to a mean increase of 36.6 kg N Mg1 (2.05 lb N bu1) between 11.1 and 13.1 Mg ha1. Although this yield response does not constitute overwhelming evidence in support of N recommendations that are a linear function of yield goal, these results are important when considering the implication for site-specific N management. Should site-specific N recommendations for this field be based on yield goal (i.e., an extension of traditional N recommendations), a spatial map of which can be easily obtained by a producer with a yield monitor?
|
Although irrigation was not available to supplement rainfall in this study, soil water content (020 cm and 090 cm) was slightly greater between late July and mid-August (Fig. 4 ; coinciding with tasseling to grain fill period of crop development) than earlier in the growing season (e.g., 30 June and 8-leaf stage). Consequently, maximum yield, or grain yield at EONR (11.113.5 Mg ha1) (Table 3), met or exceeded general expectations for this central Pennsylvania production field. The below-normal early season rainfall (Fig. 4) seemed to be masked by the JulyAugust rainfall observed in 2005. The July rainfall (12.4 cm) exceeded the yield-limiting July rainfall (9.4 cm) determined by Fox and Piekielek (unpublished report, 2001), and grain yield observed in 2005 exceeded maximum grain yield (10.4 Mg ha1) observed by Fox and Piekielek. However, the amount of water available to a growing crop is not simply a function of rainfall, as illustrated by the 10 processes affecting hillslope hydrology that were identified by Ridolfi et al. (2003).
|
Soil water content along this hillslope was less variable on 30 June, corresponding to dry soil conditions, than on 25 July and 2 September, corresponding to wetter soil conditions (Fig. 1b and 1c). Famiglietti et al. (1998) observed similar behavior in water variability in the surface soil (05 cm) along a 200-m hillslope near Austin, TX, attributing spatial heterogeneity in soil moisture under wet conditions to soil variability, and noted that joint influences of topography and soil properties contributing to more similar soil water content along the hillslope under drier conditions. Although we have not conducted a detailed analysis to describe the processes affecting water redistribution for this hillslope, the changes observed in soil water content (Fig. 1b and 1c) are suggestive of processes described by Famiglietti et al. (1998).
Perhaps less than 100% of rainfall in early July infiltrates at Location 8, whereas infiltration may equal rainfall at Location 5 and/or the soil at Location 5 is a beneficiary of subsurface redistribution of rainfall resulting in a change in profile soil water content almost equal to rainfall, despite evapotranspiration by a growing corn crop. Greater available water in the 0- to 90-cm soil profile during a growing season with less-than-normal precipitation translated to greater grain yield at EONR (Fig. 5 ). As the net change in soil profile water content between 30 June and 25 July increased from 0.5 cm to 7.3 cm, grain yield increased from slightly more than 11 Mg ha1 to slightly more than 13 Mg ha1 (r2 = 0.66; P > F = 0.04) (Fig. 5). Relationships between grain yield at EONR and the change in soil profile water content (Fig. 5) and between EONR and grain yield at EONR (Fig. 3) suggest that there might be a significant relationship between EONR and the change in soil water profile content.
|
|
These results underscore the importance in understanding the spatial variability in EONR, not simply understanding or capturing yield variability. Although there was evidence of a relationship between EONR and grain yield at EONR (Fig. 3) and a significant relationship between grain yield at EONR and the change in soil profile water content (Fig. 5), the defining relationship seems to be the one between EONR and the change in soil profile water content between 30 June and 25 July (Fig. 6). The driest part of the growing season with respect to soil water content (30 June, Fig. 4a) coincided with the period of rapid vegetative growth in corn and a very high water demand by the crop. The quite variable EONR among these 10 within-field locations (Table 3) identifies an opportunity for site-specific N management in the Northeast USA. If, as the results from this study and previous studies at this same research farm suggest, soil water availability during July is an important indicator of EONR, the practical question is: How will a producer obtain this information? Adequate soil characteristic data might be obtained through soil electrical conductivity maps and selected soil sampling and analyses to verify the electrical conductivity map. Fine-resolution topographic maps are becoming increasingly available, and Schmidt et al. (2003) demonstrated how a topographic map might be obtained with repeated passes with a typical agricultural global positioning receiver. Distributed hydrologic models could then be used to estimate soil moisture variability along a hillslope (Famiglietti et al., 1998). Although these various tools are not available to producers for site-specific N management, future research should emphasize an evaluation of site-specific EONR and provide the direction for the development of the appropriate producer tools. Development of appropriate response models will not only depend on the correlative relationships observed in this study, but also on an improved understanding of the causal relationships between soil physical and/or landscape characteristics and EONR. Eliciting these types of causal relationships with EONR, which may be unique to individual fields or perhaps to a geographic region, is the challenge for site-specific N management research.
| NOTES |
|---|
|
|
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. E. Dellinger, J. P. Schmidt, and D. B. Beegle Developing Nitrogen Fertilizer Recommendations for Corn Using an Active Sensor Agron. J., October 21, 2008; 100(6): 1546 - 1552. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. A. Hernandez and D. J. Mulla Estimating Uncertainty of Economically Optimum Fertilizer Rates Agron. J., September 1, 2008; 100(5): 1221 - 1229. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. L. Ping, R. B. Ferguson, and A. Dobermann Site-Specific Nitrogen and Plant Density Management in Irrigated Maize Agron. J., June 23, 2008; 100(4): 1193 - 1204. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Crop Science | Vadose Zone Journal | |||
| Journal of Natural Resources and Life Sciences Education |
Soil Science Society of America Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||