| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|
|
|
|||||||
|
|||||||||
a ERS-USDA, 1800 M St., N.W., Washington, DC 20036
b Dep. of Agricultural Economics, Purdue Univ., 403 West State St., West Lafayette, IN 47907-20566
c Dep. of Soil, Water, and Climate, Univ. of Minnesota, 1991 Upper Buford Circle, St. Paul, MN 55108
* Corresponding author (dlambert{at}ers.usda.gov)
Received for publication January 3, 2005.
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONMATERIALS AND METHODSRESULTS AND DISCUSSIONCONCLUSIONSREFERENCES |
|---|
Abbreviations: EOR, economically optimal rate • EONR, economically optimal rates for N • EORP, economically optimal rates for P • GEO, Geostatistic Regression Model • LM, Lagrange Multiplier • NPV, net present value • PSNT, Pre-Sidedress Nitrogen Test • SSM, site specific management • UNI, uniform application • VRA, variable rate application • VRN, variable rate nitrogen • VRT, variable rate technology • WF, whole field
|
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTS AND DISCUSSIONCONCLUSIONSREFERENCES |
|---|
Until recently, multi-year on-farm production data from controlled variable rate technology (VRT) experiments have been limited. Studies that have looked at site-specific yield response to specific inputs over two or more growing seasons are more frequent because georeferenced panel production data are increasingly available. Finck's (1998) article reported results from an on-farm VRT-NP experiment with corn and soybean over three seasons in Illinois, USA. Bongiovanni and Lowenberg-DeBoer (2002) analyzed variable rate nitrogen (VRN) profitability and site-specific crop response stability for corn over two growing seasons in Argentina. Kasper et al. (2003) used 6 yr of corn data to estimate corn response with respect to topography and elevation. Mamo et al. (2003) evaluated a VRN trial for corn and an N loss inhibitor over three corn production cycles in Minnesota, USA. Like Kasper et al., they found that corn yield response to N varied spatially and was temporally unstable. Swinton et al. (2002) analyzed a VRN trial on 14 farms in Michigan, USA, between 1999 and 2001 (in all, a total of 24 site years). They also found that corn yield response to N was not stable over time.
The objectives of this paper are twofold: (i) to test whether yield response of corn to N and P and soybean response to P exhibit stable temporal and spatial patterns, and (ii) to evaluate the profitability of a variable rate nitrogen and phosphorus fertility management strategy over a 5-yr corn–soybean rotation using this response information. The VRT study spans 5 yr of production data from a corn–soybean rotation in Minnesota, USA. The analysis focuses on the temporal and spatial variability of crop response rather than the temporal and spatial stability of yield itself. The spatial and temporal analysis of crop response to variable rate N and P identifies landscape positions where corn and soybean response to inputs is similar between production seasons. Variable rate application profitability is compared with a uniform (UNI) N and P management strategy with an ex post, multi-season net present value (NPV) comparison. Spatial correlation between yield observations is modeled using a geostatistically weighted regression (Cressie, 1993). The site-specific management (SSM) hypothesis (Whelan and McBratney, 2000) is tested over the production surface and between growing seasons. The null hypothesis states that the optimal risk management strategy for producers is a uniform input management strategy, given temporal variations in weather patterns. The findings have implications with respect to N and P fertilizer management, and identification of areas where yield response is temporally stable and spatially homogenous. The results may interest corn and soybean producers, researchers, crop consultants, and the persons who advise them.
|
MATERIALS AND METHODS |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTS AND DISCUSSIONCONCLUSIONSREFERENCES |
|---|
The experimental design was 3 replications of 13 treatments in a split plot arrangement of a randomized complete block design where P was the main plot. Nitrogen rates were randomized within the P treatments (Fig. 1 ). The P treatments were randomized within each replication. Treatments were applied at constant rates in strips across the entire field. All treatments were repeated three times (a total of 39 strips). Fertilizer N treatments were 0, 67, 112, 157, and 202 kg ha–1 of N as either anhydrous ammonia (fall 1996) or urea (fall 1998 and 2000), and P treatments were 0, 56, and 112 kg ha–1 of P, as triple super phosphate. In the fall of 1998 and 2000, P fertilizer was applied to the same areas of the field that had received the treatments in 1996. The experiment was designed to accommodate extra N treatment strips that would ensure zero N rate treatments were not placed in the same location as previous zero N rate treatments. Other N treatments were the same each corn year.
|
The experimental site was partitioned into 69 sub-blocks for purposes of spatial yield response analysis (Fig. 1). Previous applications of this methodology are found in Mamo et al. (2003), Hurley et al. (2003), and Malzer et al. (1996). For each corn year, 897 observations were available. Therefore, there were 13 observations available for estimating corn response to N and P treatments in each sub-block. Only 621 observations were collected with the research combine during the soybean years. Nine observations were available to estimate soybean response to P in each sub-block.
Because a yield response is estimated for each sub-block, economically optimal rates for N (EONR) and economically optimal rates for P (EORP) can be estimated for each sub-block, along with corresponding profit maximizing yields and profits. This facilitates ex post comparisons of variable rate fertilizer strategies with uniform fertilizer management strategies.
Corn and Soybean Yield Response Model
Because inputs are economically quantifiable, response functions facilitate comparison between input changes and the cost of making those changes. A quadratic response function is used to estimate corn response to N and P and soybean response to P. Quadratic functions have theoretical appeal because they are interpretable as second-order approximations of any function in general, and crop response functions to inputs in particular (Hurley et al., 2003). Alternatively, a quadratic response and plateau model could be considered. However, this is computationally onerous given the multiple response zones and semivariogram estimation. Most of the recent work estimating site-specific crop response functions used single, quadratic equations to estimate plant response to inputs over space (e.g., Bruulsema et al., 1996; Malzer et al., 1996; Bongiovanni, 2002; Swinton et al., 2002; Hurley et al.,2003; Mamo et al., 2003; Anselin et al., 2004; Lambert et al., 2004). For crop response economic analyses, quadratic functions are suitable because they are concave functions with f'> 0 and f''< 0, which permits diminishing marginal returns to inputs. The estimation of economically optimal input rates is also tractable because a closed-form solution to f' exists.
For each crop year, the response equations are:
 |
 | [1] |
t(s)
(0, 
t) is a disturbance term that may be spatially correlated; N and P are applied inputs at position s; and ds,z is a dummy variable identifying sub-blocks z = 2...69, with the first block the reference group. Parameters of the sub-block dummy variables are constrained as 
z = 269 
i,z = 0 to capture differences between site-specific crop response functions and the whole field (WF), average crop response. This constraint tests Whelan and McBratney's (2000) SSM hypothesis because significant response differences are a precondition for VRT profitability. Significant t tests for these coefficients indicate whether the response in a sub-block is different from the average WF response in a given year. Rejection of the null hypothesis that a sub-block is not different from the WF response to N or P is evidence in favor of VRA management. The subscript t indicates which crop is grown in year t, t = 1997...2001, with corn starting in 1997 and soybean rotating in the even years thereafter. Because N was not the focus of the soybean study, P was the only fertilizer variable included in the soybean regression equation. Equation [1] was estimated separately for each year because the number of corn and soybean observations was not the same for each year.
Modeling Spatial Covariance
A geostatistic regression model (Cressie, 1993) was used to estimate yield response when spatial correlation between observations was significant. Because yield observations were not equally dispersed, the assumption that equal weights could be placed on neighboring yield points was not tenable. Therefore, a spatial econometric approach using a spatial weighting design based on a lattice contiguity matrix was not appropriate (for example, Anselin et al., 2004). Hurley et al. (2003) and Lambert et al. (2004) have used the GEO method to estimate crop response to variable rate inputs. Lark and Wheeler (2003) used a variant of this model to estimate yield response using yield monitor data. Robust empirical semivariograms (Cressie, 1993) were estimated using the corn or soybean residuals from Eq. [1]. The functional form used to model covariance between observations is 
(s) = c0 + 
1 – exp (–[||s||/a]
); where ||s|| is the Euclidean distance between observations; and c0, , and a are the nugget, sill, and range parameters, respectively. The innovation of this model is the shape parameter, . To ensure admissibility of this function, is restricted to be positive. Experience with the Windom data has shown that empirical semivariograms are sometimes intermediate between Gaussian and exponential functional forms.
Regression Estimation and Diagnostics
A Breusch-Pagan (Breusch and Pagan, 1980) Lagrange Multiplier (LM) test for heteroskedasticity was conducted for each crop response equation, each season. When the null hypothesis of homogenous variance between observations was rejected, Eq. [1] was re-estimated using General Method of Moments (GMM) with White's robust standard errors. The GMM procedure is useful when the exact form heteroskedasticity is unknown (SAS Institute, 1999). Furthermore, no assumptions about the error distribution are required except that errors are independent and identically distributed (Mittelhammer et al., 2000).
Cressie's (1993) weighted nonlinear least squares procedure was used to fit (s) to the regression residual semivariograms. The F test of each semivariogram regression indicates if there is significant spatial structure in the residuals of each response equation. When spatial structure was evident, Eq. [1] was estimated using geostatistically weighted GMM.
Ex Post Marginal Analysis of Management Strategies
Uniform and VRT profitability is compared using a partial budget (Boehlje and Eidman, 1984), and marginal analysis is used to estimate net returns from applied N and P (Beattie and Taylor, 1985). Marginal analysis shows that profit is maximized when the value of the increased yield from added N or P equals the cost of applying an additional unit of fertilizer, or when the marginal value product equals the marginal factor cost.
The profitability analysis is ex post because it assumes that the producer knows, in hindsight, yield response and optimal N and P applications (Bullock and Bullock, 2000). Ex post analyses are useful as starting points in economic assessments of solutions to management problems (Anselin et al., 2004). While ideally such evaluations should be based on ex ante decision-making, ex post estimates provide insight into possible outcomes of VRT-UNI comparisons. If ex post results are not profitable, then ex ante results are unlikely to be profitable. Alternatively, if the ex post approach does show profitable results there still remains the question of how well it would do in an ex ante decision making context.
A spatial optimization framework proposed by Bullock et al. (2002) is modified to accommodate the multi-period corn–soybean rotation. Over the course of the experiment, N and P were only applied in the fall before corn was planted (except for the 2001 year when bad weather conditions postponed fertilizer N application). Therefore, the P-for-soybean decision is made when fertilizer amounts are determined for corn.
The producer's maximization problem in the mth corn–soybean rotation and zth location is
 |
)–m is a discounting factor with an opportunity cost of capital of 7.5% (Chiarella et al., 2002); yc,zm (ys,zm) is the estimated corn(soybean) yield from Eq. [1] in sub-block z, rotation m; xz,m = [Pz,m, Nz,m]' is a k x 1 input vector applied at the beginning of rotation m on sub-block z; F are fixed costs; g and v are quasi-fixed, site-specific information and variable rate application costs, respectively; r is a k x 1 vector of variable costs for the k inputs applied at the beginning of the corn year in rotation m; and I(m)(•) is an indicator variable equal to one when g and v costs are incurred by the producer, zero otherwise. For example, soil tests or maps (g) may only be needed every 4 yr or more, or v may be incurred only during corn years.
The sum of NPVs over the m corn–soybean rotations and z locations is the cumulative expected NPV over the course of the experiment. Here, total expected NPV is E[NPVTOTAL] = z = 169
z[E(NPVz,1997–1998) + E(NPVz,1999–2000) + E(NPVz,2001)], with z the portions of the field covered by sub-block z (z = 1/69). When the NPV from the UNI management strategy is compared with VRT cash flow, fixed costs (F) disappear because they are the same for both. Analytical derivatives solve optimal N and P factor demand levels for each sub-block in each crop year.
Price Data Used in the Ex Post Profitability Analysis
The average market prices (1997–2001; NASS, 2002) for corn and soybean observed in Cottonwood County, Minnesota, were used as output prices ($0.0787 and $0.175 kg–1, corn and soybean, respectively). Input prices were $0.077 and $0.118 kg–1 for N and P, respectively. Information costs for SSM included costs for map making and management zone development ($2.96 per map), and soil sampling, including lab and collection fees ($13.59 ha–1). These are the costs reported in the 1997 Akridge and Whipker dealership survey. Variable rate application costs ($13.21 ha–1) were the average of the VRT costs reported in Akridge and Whipker (1997, 1999) and Whipker and Akridge (2001). A uniform application cost of $9.88 ha–1 is assumed for the UNI budget analysis (Aghib and Lowenberg-DeBoer, 1999). The input application costs for VRT and UNI are doubled because P and N are usually applied at different times. Because the UNI strategy uses soil test information to follow extension recommendations for P, a whole field soil test fee of $0.82 ha–1 is charged under the uniform strategy. Mapping and soil test fees are charged once, in the 1997 year. These products are assumed to have a useful life of approximately 4 yr (Swinton and Lowenberg-DeBoer, 1998). Uniform and VRT fees are only charged before the corn production years. The NPV of the VRT strategy is compared with the NPV of returns from the UNI N and P fertilization strategy.
|
RESULTS AND DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTS AND DISCUSSIONCONCLUSIONSREFERENCES |
|---|
|
|
(s). The 1999 and 2001 corn years were the only years where spatial error dependence was significant (F test, F = 2.37, P = 0.08 and F = 4.76, P = 0.005, respectively). This is not surprising because spatial error dependence is usually considered an omitted variable problem (Anselin, 1988). The sub-block regressions are essentially fixed effects models because block-specific dummy variables identify intercept, linear, quadratic, and interaction response terms in each sub-block. Allowing each sub-block its own response explained most of the spatial error dependence in the corn 1997 and soybean 1998, 2000 years. Therefore, the GEO method was only used to estimate corn response in 1999 and 2001.
Spatial and Temporal Response Variability and Stability
A graphical approach is used to convey spatial and temporal relationships of corn response to N and P and soybean response to P (Figs. 3
and 4)
. The t test statistics associated with the linear N and P coefficients for corn and soybean in Eq. [1] determine whether crop response to N and P was significantly different from the average whole-field response. The bars associated with each point estimate are 90% confidence intervals. The complete set of response coefficients are found in Lambert (2005). Quadratic and interaction terms are important with respect to economic analysis of crop response, but only linear response terms are highlighted here. The linear term is the plant response to the first kilogram of fertilizer and, as such, is the key to understanding profitability. In profit maximization, the interpretation becomes how much profit is increased when an extra kilogram of input is applied. If the first kilogram of fertilizer is not profitable, then it is unlikely that any fertilizer will be profitable. The quadratic term is important for determining the rate at which response diminishes at higher application rates and, consequently, the peak yield. Linear and quadratic terms were highly negatively correlated for N and P in all years for both crops (Pearson's r > 0.93). This means examination of the quadratic coefficients shows similar patterns, but in a mirror image.
|
|
5%). These findings were expected because soil high in organic matter has a high potential to supply N through mineralization (Mamo et al., 2003). Spatial patterns where corn response to P was significantly different from the WF average response were clearly discernable (Fig. 3), especially in the bare hilltop area dominated by a Clarion-Swanlake loam series soil with a 6 to 12% slope (Sections 2, 3, Rows 10–17), the high organic matter portion (Section 2, Rows 1–6), and in Section 1 in the low-lying seepage areas of the field (Rows 1–8, Jeffers variant clay loams series) for the 1997 and 1999 corn years. In the 1997 and 1999 corn years, 32 and 40% of the sub-blocks exhibited P response levels significantly different from the WF average (Table 1). In the 2001 season, P response heterogeneity for corn increased to 49%. Locations where corn response to P was lowest were in areas where soil organic matter or P soil test levels were relatively high in Sections 2 and 3 (Fig. 2).
Corn response to P in Section 3, Rows 19 to 23 was temporally and spatially stable. The point estimates of all years rest in the same confidence interval band and were highly significant (Fig. 3). On the other hand, corn response to P was not temporally stable in Section 2, Row 5 because the confidence intervals of those point estimates do not overlap. Likewise, 1999 corn response to P in Section 2, Row 4 was significantly different from corn response to P in Section 2, Rows 20 to 23 in the same year with respect to location. Corn response profiles along the P treatment sections suggest that although yield response was significantly different between crop rotations, similar patterns were evident in physically distinct portions of the field (Fig. 3).
In the low-lying seepage area dominated by Jeffers clay loam and Clarion loam (2–4% slopes) series soils (Site A) and hillside seepage areas (Sites B, C, Clarion loams, 2–4% slope and Webster-Delft clay loam soils, respectively) in P treatment Section 1, the corn response pattern to applied P in 2001 was highly variable compared with the 1997 and 1999 seasons (Fig. 3). The high variability associated with P in Section 1 was most likely due to the fact that only three P levels were used in the treatment block, and that this portion of the field has poor drainage capability. In Section 3, soil in the portion of the field where the old dairy entrance was located (Site G) had a high P content. Corn response to P in this landscape position was remarkably similar across all corn years. The south end of Section 2 is dominated by a Canisteo clay loam series soil (fine-loamy mixed [calcareous] Typic Haplaquolls, Rows 20–23), and corn response to P was also stable across all years. But in Section 1, the corn response pattern to applied P in 2001 was highly variable compared with the 1997 and 1999 seasons. This area of the site corresponds with high soil organic levels (>4.5%) and soil test P levels of 5 to 11 mg kg–1 [ppm]. In P treatment Sections 2 and 3, a prominent physical feature is the eroded hilltop area composed of a Clarion-Swanlake loam series soil with 6 to 12% slopes (Sites E2 and F2; between Rows 10 and 17). Hillside seepage areas dominated by a Jeffers variant clay loam series soil (2–4% slope) are also present on the downward slope of these transects (Site E1, F1, Rows 10–14). During the 2001 corn year corn response to P was highly variable in these spots, which are associated with low organic matter and low soil P. The organic sink in the north end of Section 2 (Site D) also exhibited stable and relatively small P response profiles for corn where soil organic matter ranged between 5 and 9% and soil test P levels were between 8 and 11 mg kg–1 [ppm].
When spatial P response variability was considered across the corn years, corn response to P was similar in 51% of the blocks. Here, response similarity means that the SSM null hypothesis for a given sub-block was equally rejected across all years. For soybean, 57% of the blocks responded similarly to P between growing seasons. But for corn response to N, only 26% of the sub-blocks responded similarly between growing seasons. The relative instability of corn response to P in 2001 was likely caused by the interaction between the number of P treatments and the poor drainage capability of some parts of the field. The relative stability of P response was likely the result of its mobility in soils and the physiological effects of P on root system growth.
Similarities in corn response to P can be expected due to P soil carryover dynamics and placement, but consistent spatial response patterns are less discernable for corn response to N. Nitrogen is more susceptible to hydrologic conditions that are affected by annual precipitation and topography interactions. These findings for N correspond with the Mamo et al. (2003) N study; temporal inconsistencies in response can be attributed to organic matter mineralization, denitrification, and water availability, which in turn are influenced by topography. These findings for N are also consistent with the Fiez et al. (1994) conclusions about N response and topography relations.
The most striking temporal response trends are observed in the soybean response profiles (Fig. 4). The curves of soybean response trends within transects are remarkably similar. Although the magnitude of soybean yield response to P was different between years, the same underlying effects of P response dynamics are present. These relations were strongest in Sections 2 and 3 (Fig. 4). Soybean response to P was different from the WF average in 33% of the sub-blocks in 1998, and different from the WF average in 62% of the sub-blocks in the 2000 season. Soybean response to P was significantly different from the WF average in the high organic matter portions of the field and the eroded hilltop area for the 2000 season. For both soybean years, response to P was significantly different from the WF average along the entire transect of Section 3, excluding the southeastern upper portion of the field where the entrance to an old dairy barn was once located (Site G). Because of the localized distribution of manure from the barn, P soil test values were very high in this position. These findings are consistent with Wittry and Mallarino's (2004) VRT-P study for corn and soybean. They found that plant response to P was only significant in areas of the field recording <15 mg kg–1 [ppm] P. For P, 33 and 62% of the sub-blocks reported significant response to the first kilogram of P for soybean. As a response sensitivity test, N was included in the soybean response equations. Soybean marginal response to the first kilogram of N was significant in only 4 and 14% of the sub-blocks for the 1998 and 2000 soybean years, respectively.
A side-by-side comparison of corn and soybean response patterns to N and P is possible by standardizing the response coefficients for both crops across all years as N(0,1) with their respective averages and standard deviations (Fig. 5
). In the third P section, similar response patterns to P are observed in corn and soybean. Likewise, similar response patterns are also discernable in the second P section (Site E1). In the first section, comovement of response patterns is less noticeable but still present, particularly in the areas between the hillside seepage areas (Sites B, C, Rows 12–19). Comovement of N across corn years is less discernable, except perhaps in Section 2, Rows 1 to 15. By these measures the spatial patterns of temporal (in)stability become clear. For corn response to N, each year is different and response patterns appear to be random. For P, moving west to east, P response over time for corn and soybean converges toward a stable trend over landscape position and cropping years.
|
Economically Optimal Nitrogen and Phosphorus Rates
In S1, the field average level of applied N under the VRA program was higher than the extension recommendation amount for all years by 15, 44, and 27% for 1997, 1999, and 2001, respectively (Table 2). The EORP was lower for the 1997 and 2001 years by 2 and 8%, but only slightly higher in the 1999–2000 rotation by 1%. On average, EONR was higher than the extension N recommendation (mean ± SD, 115 kg ha–1 ± 13), but EORP was slightly lower (64 kg ha–1 ± 5) than extension recommended P amounts for Scenarios 1 and 2 (Table 2). In the third corn year, 13% of the cells had zero EONRs. One location for these recommendations was in the lower corner section of Section 1 (Rows 21–23), which is dominated by a Webster-Delft clay loam series with poor drainage capacity. Other zero N levels occurred on the margins of the eroded hilltop area, and the downward slope just north of this position (Section 2, Row 8 and 15). This indicates that in these areas the mineralization potential of N by organic matter is substantial. In the control strip, corn yield was 7843 kg ha–1 in this section, indicating the fertility of the position. Other zero N rates occurred along the slopes of the eroded hilltop area (Sites E1, F1). This landscape position is prone to seepage problems during wet years and is dominated by a Jeffers clay loam soil series relatively high in organic matter. In these positions, the mineralization potential of N by organic matter content is also high.
|
The 1997 EONR rate in S4 is much lower than the S1 optimization EONR by, on average, 45%. This difference is also due to the N x P interaction term. The cause of this difference is the negative sign on the N x P interaction term in the 1997 corn year in S4 (Table 1). When P (or N) is not allowed to vary, the first-order conditions of the optimization problem are quite different than those of the joint N and P optimization problem. For N, holding P constant, solving the FOC for the optimal N rate is
with B the discount factor, rN the nitrogen price, pc the corn price, 
the phosphorus extension recommendation, and the subscripts on the betas indicating N linear (N) and quadratic (N x N), and N x P interaction terms. When ßNxP < 0, ßNxN < 0, and ßN > 0, N* will decrease. This only happens in the 1997 year of S4, and is why the EONR is much lower than the S1 optimization scenario. For the N-UNI, variable P scenario, the FOC is slightly different because of the additional soybean term. Assuming that soybean linear (quadratic) terms are positive (negative), P increases because of the revenue boost from soybean, while N is held constant.
Marginal Analysis: Ex Post VRT/UNI Yield and Profitability Comparison
The SSM null hypothesis was rejected in the first and second scenarios (Table 2). The NPV of the VRA-NP strategy was significantly greater than the UNI strategy by $28.38 ha–1 in the S1 (baseline) comparison (Paired t test, T = 3.06, P = 0.003). The largest negative cashflow for VRA occurred in the 1997 corn season ($12.56 ha–1). Cashflows were lower for this season in part because of initial capital outlay to maps, management zone identification, and soil test costs, combined with the VRA-P application fee. In a cost sensitivity analysis, when the initial costs were annualized over the 5 yr, the returns for VRA exceeded those of the uniform application every season. Differences between UNI and VRA-NP cashflows were significant only for the corn years. In the extremely wet year of 2001, the cashflow margin between VRA-NP and UNI was even greater than the previous seasons.
It is widely known that information and application costs have constrained widespread adoption of VRA technologies (Griffin et al., 2004). It is also generally recognized that applying the right amount of inputs in the right locations will produce higher yields and increases in returns over input costs compared with uniform application of inputs. But increased returns may not cover information and VRA costs with current technology. In S2, the SSM null hypothesis is rejected because the expected cumulative NPV of the VRA-NP strategy was significantly higher than returns to the UNI-NP strategy by $68.47 ha–1 (T = 7.38, P < 0.0001) (Table 2). When VRA does not require higher information or application costs, all rotation NPVs were higher under the VRA-NP strategy compared with the UNI strategy.
The third (S3) and fourth (S4) scenarios compare VRA-P or VRA-N strategies combined with constant N or P amounts applied at the extension recommendation, respectively (Table 2). The SSM null hypothesis is not rejected in Scenario 3. Returns to VRA were still higher than the UNI strategy ($0.25 ha–1), but not significantly (T = 0.08, P = 0.93). In S4, the SSM null hypothesis is not rejected because the profit margin significantly favored the UNI strategy by $19.81 ha–1 (T = 13.27, P < 0.0001). The lack of VRA-N profitability is likely due to the lower corn yield response to site-specific N in 1997 and 2001. The profitability of site-specific P is increased by returns in soybean years. However, the VRA-P, uniform N, net present value would be negative were it not for the revenue gained from the soybean years. Additionally, N is more volatile than P, making it more difficult to spatially manage over multiple growing seasons without recourse to address spatial and temporal variability of soil mineralization or weather effects on response.
Because of the unbalanced nature of the experimental design (3 yr of corn and 2 yr of soybean), the partial budget results may be biased in favor of VRA fertilizer application. To test this assumption, the marginal analysis was conducted using the corn years only. Results showed that although the nature of site-specific relations was limited with soybean, they do not change the results of S1 or S2: VRA profits were greater than UNI by about $11 and $18 ha–1, respectively. However, the results do change in S3: VRA was no longer profitable compared with the uniform recommendation (difference = –$6.76 ha–1). When considering the corn years only, applied P decreased, respectively, in the 1997 and 1999 years by 10 and 15% in S1. Optimal P applications decreased by 51 and 44% for the 1997 and 1999 years, respectively, in Scenario 3. This occurred because the first-order condition terms for soybean P drop out of the corn–soybean rotation optimization problem. When expected revenue from soybean is not jointly considered with corn in the planning horizon, conclusions about VRT as applied to multiple corn–soybean rotations may be underestimated.
|
CONCLUSIONS |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTS AND DISCUSSIONCONCLUSIONSREFERENCES |
|---|
In an ex post cash flow analysis, discounted returns and corn–soybean yields from a VRA-NP fertilizer program were $28 ha–1 greater than returns and yields from a UNI management strategy. Corn and soybean yields also increased under the VRA-P, uniform N strategy (8354 and 2939 kg ha–1, respectively) compared with the UNI strategy (8296 and 2881 kg ha–1, respectively). However, the $0.25 ha–1 profit margin between the VRA-P and UNI strategy was not significant. In this study, returns to a VRA-N strategy were significantly less than the UNI strategy, indicating that spatial management of this input over multiple growing seasons is more difficult than spatial management of P.
Ignoring the opportunities gained from knowing timely nutrient level information may exaggerate the temporal level of yield response risk, as well as interactions between temporal and spatial yield response to risk. With respect to corn yield response, N clearly drives the profitability results. Today, presidedress nitrate tests (PSNT) are commercially available. The PSNT is just one example of many that might be used to minimize temporal risk. Although these tests are designed to predict plant N needs, they do not always work as well as most producers would like. The relative success of these tests is in part due to climatic and soil variability from one region to the next. But when PSNT tests are on the mark, they may substantially reduce the temporal risks of N management. Later-season N applications would be another strategy to deal with the risk caused by N instability over multiple growing seasons. For farmers whose income depended on revenue from corn production only, appropriate timing and placement of N would play more of a role with respect to risk management. But, if response functions are relatively flat (as they are in many parts of the field in this study), then variable rate N application will not substantially reduce risk, even if these parameters could be determined.
Future research could uncover the temporal and spatial mechanisms regulating crop nutrient removal and nutrient carryover dynamics. This next step would require estimation of nutrient carryover functions jointly with site-specific estimation of crop response functions. Results of this analysis are a starting point for management zone delineation, based not only on yield patterns or soil characteristics, but also on response variability and fertilizer carryover dynamics.
|
ACKNOWLEDGMENTS |
|---|
|
REFERENCES |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTS AND DISCUSSIONCONCLUSIONSREFERENCES |
|---|
This article has been cited by other articles:
|
|
 |
|
  P. M. Kyveryga, A. M. Blackmer, and T. F. Morris Disaggregating Model Bias and Variability when Calculating Economic Optimum Rates of Nitrogen Fertilization for Corn Agron. J., June 5, 2007; 99(4): 1048 - 1056. [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 | |||
