1. Introduction
Potassium (K) is an important nutrient for upland cotton (Gossypium hirsutum L.) but can be difficult to manage over time (Howard et al., Reference Howard, Essington, Hayes and Percell2001). In the late 1980s and early 1990s, there were frequent reports of late-season K deficiencies in the U.S. Cotton Belt (Maples, Thompson, and Varvil, Reference Maples, Thompson and Varvil1988; Mullins, Burmester, and Reeves, Reference Mullins, Burmester and Reeves1997), resulting in numerous agronomic studies on cotton response to K while considering K levels in the soil that are readily available for plant consumption (soil K) from the previous production year (i.e., carryover) (Essington et al., Reference Essington, Howard, Savoy and Lessman2002; Howard et al., Reference Howard, Gwathmey, Roberts and Lessman1998; Mullins, Schwab, and Burmester, Reference Mullins, Schwab and Burmester1999). Findings from these studies triggered extension personnel in several southeastern states to recalibrate K recommendations to maintain adequate soil K levels for cotton production and circumvent yield loss because of late-season K deficiencies. Methods to analyze soil K levels were developed more than 60 years ago (Mehlich, Reference Mehlich1953), but gathering this information using soil tests and applying the information to make K applications is still a growing practice among cotton producers (Zhou et al., Reference Zhou, English, Boyer, Roberts, Larson, Lambert and Velandia2015).
If a profit-maximizing producer considers information on soil K levels before applying a K rate, the producer's decision framework changes from maximizing net returns in a given year to maximizing the present value of net returns (NPV) over a planning horizon. The switch in the producer's decision framework occurs because application decisions in a given year are based on their application rates from the previous year. Economists developed dynamic programming models to determine nutrient application rates that maximize NPV with nutrient carryover (Fuller, Reference Fuller1965; Heady and Dillon, Reference Heady and Dillon1961; Kennedy, Reference Kennedy1986; Kennedy et al., Reference Kennedy, Whan, Jackson and Dillon1973; Stauber, Burt, and Linse, Reference Stauber, Burt and Linse1975). The dynamic programming approach separates the producer's planning horizon into discrete sequential maximization problems, solved by deriving optimality conditions for each period (Bellman, Reference Bellman1957). The single period application rates are conditioned on some knowledge of soil carryover levels prior to planting.
These dynamic models have been adapted to different crops and nutrients (Ackello-Ogutu, Paris, and Williams, Reference Ackello-Ogutu, Paris and Williams1985; Jomini et al., 2001; Kennedy et al., Reference Kennedy, Whan, Jackson and Dillon1973; Lambert, Lowenberg-DeBoer, and Malzer, Reference Lambert, Lowenberg-DeBoer and Malzer2007; Lanzer and Paris, Reference Lanzer and Paris1981; Park et al., Reference Park, Stocker, Hattey and Turner2007; Schnitkey, Hopkins, and Tweeten, Reference Schnitkey, Hopkins and Tweeten1996; Segarra et al., Reference Segarra, Ethridge, Deussen and Onken1989; Watkins, Lu, and Huang, Reference Watkins, Lu and Huang1998). However, Harper et al. (Reference Harper, Lambert, Larson and Gwathmey2012) were first to use a dynamic programming model to determine K application rates that maximize NPV for upland cotton production. They developed an application for valuing the information from soil testing in cotton production by considering multiple information scenarios in the dynamic programming framework. Harper et al. (Reference Harper, Lambert, Larson and Gwathmey2012) used 3 years of data on K soil fertility and cotton lint yields in Tennessee to determine the applied K rate that maximized NPV when soil K levels were considered and were not considered in the applied K rate decision. They found that using the soil K level information annually would increase NPV relative to not considering the soil K level information over a 5-year horizon. Additionally, using soil K carryover knowledge reduced the amount of annual K that was applied and the soil K carryover levels, which may be helpful for reducing off-site K leaching.
A common assumption in the existing dynamic programming literature is that producers annually soil test to update their soil nutrient information. However, gathering information about soil nutrient variability on an annual basis might not increase the NPV (e.g., by reducing fertilizer costs) enough to pay for the cost of gathering soil nutrient information, especially for K because it is immobile in the soil profile (Walworth, Reference Walworth2011). Currently, state extension recommendations encourage cotton producers to update K soil test information from annually to every 3 years in the Southeast (Kissel and Sonon, Reference Kissel and Sonon2011; Mylavarapu, Reference Mylavarapu1997; Savoy and Joines, 2013). These recommendations vary across states because production factors such as cropping intensity, soil type, tillage practices, and weather conditions play an important role in determining the length of time to wait until retesting soil nutrient levels. Nonetheless, the common assumption in the literature of soil testing annually may not be appropriate for all nutrients and crops to maximize NPV.
Lambert et al. (Reference Lambert, English, Harper, Larkin, Larson, Mooney, Roberts, Velandia and Reeves2014) used survey data of cotton producers in 13 southern states to determine the factors affecting the length of time between updating soil test information (temporal frequency) for precision soil sampling. They found that farm size, land ownership, farm location, and farming experience were correlated with the temporal frequency that producers tested soils. Overall, cotton producers who adopted precision soil sampling indicated they retested soils on average every 2.5 years, which is within the range encouraged by Southeast extension agronomists. Moreover, Lambert et al.’s (Reference Lambert, English, Harper, Larkin, Larson, Mooney, Roberts, Velandia and Reeves2014) research indicates that annual soil testing might not be a profit-maximizing frequency.
The objective of this research was to determine the K application and temporal frequency for obtaining K soil test information that maximizes NPV to K in cotton production over a 10-year planning horizon. Optimal K rates and NPV were determined ex ante for five soil testing schedules of varying temporal frequencies using Kennedy's (Reference Kennedy1986) dynamic programming framework. The conceptual modeling of this study extends the literature by considering the temporal frequency of obtaining soil test information in a dynamic programming framework to determine optimal K rates and expected returns for cotton production. The results can guide producers and extension personnel on optimal K rates and length of time between soil tests for upland cotton growers.
2. Empirical Framework
2.1. Dynamic Programming Model
A risk-neutral, profit-maximizing cotton producer chooses an amount of K to apply (At ) at the beginning of each production year (t = 1,. . .,T [T = 10]), conditioned on some knowledge of soil K carryover (Ct ), that maximizes the NPV of returns to K over a planning horizon (Kennedy, Reference Kennedy1986; Kennedy et al., Reference Kennedy, Whan, Jackson and Dillon1973). This producer also selects the optimal temporal frequency j of soil testing by choosing some discrete number of years between retesting soils. Five soil testing schedules of varying temporal frequency from annually to every fifth year (j = 1,. . ., 5) were considered, where j is the number of years between soil testing.
At the beginning of each production year, the producer had knowledge of current elemental K prices, an expectation of cotton lint prices at harvest, cotton lint yield response to total available K, and information on soil K levels from the most recent soil test. For each soil testing schedule, the producer used this a priori knowledge to apply a profit-maximizing amount of K. The optimal temporal frequency was the soil testing schedule that provides the greatest NPV over a 10-year planning period. Therefore, the maximized NPV is calculated as follows:
$$\begin{equation}
\begin{array}{@{}l@{}} \mathop {\max }\limits_{{A_{t,j}}...,{A_{T,j}}} NP{V_j} = \sum\limits_{t = 1}^T {{\delta ^{t - 1}}N{R_{t,j}}} \\[12pt] {\rm{Subject\; to{:}}}\\[3pt] {A_{t,j}},{C_{t,j}} \ge 0\\ {C_{t + 1,j}} = {a_0} + {a_1}({A_{t,j}} + {C_{t,j}})\\ {C_{t,j}} = {\lambda _j}{Q_{t,j}} + (1 - {\lambda _j}){Q_{t - 1,j}}\\ {Q_0}\;{\rm{given,}} \end{array}
\end{equation}$$
where NPVj is the sum of discounted net returns ($/ac.) over T years for a producer following soil testing schedule j; At,j is applied K (lb./ac.); NRt,j is the net returns ($/ac.) to K for cotton production; δ is a discount factor reflecting the time value of money 1/(1 + r) t , where r is the discount rate; Ct,j is the producer's knowledge of carryover K(lb./ac.) in time period t; Qt,j is the actual carryover K (lb./ac.) in time period t; Ct +1, j is the producer's knowledge of carryover K(lb./ac.) prior to planting in year t + 1, which is a function of applied K and the producer's knowledge of soil K carryover (i.e., total K available [lb./ac.]) in year t; λj is an indicator variable that is equal to 1 in the year a producer updates his or her knowledge of soil K by soil testing, and 0 otherwise; a 0 and a 1 are estimated parameters for the linear carryover function; and Q 0 is the actual soil K level before K is applied in the first production period. Partial budgeting was used to calculate the single period net returns for a risk-neutral profit-maximizing producer, where single period net returns (NR) are the following:
$$\begin{equation}
N{R_{t,j}} = \delta p_t^c{y_{t,j}}({A_{t,j}} + {C_{t,j}}) - p_t^K{A_{t,j}} - {\lambda _j}s,
\end{equation}$$
where pc t is the producer's expectation of the price of cotton lint at harvest ($/ac.);pK t is the producer's observed price of K at the time of application ($/ac.); yt,j is cotton lint yield (lb./ac.) in period t; and s is the cost of obtaining soil test information prior to applying K ($/ac.), which only occurs in years when the producer tests the soil (λj = 1), or otherwise s = 0.
Assume the producer conducts a soil test prior to production in the first year; therefore, the producer knows the soil K level at the beginning of production year 1 and uses this information in selecting an optimal application rate in year 1 (A 1,j ). After the first production year, the producer chooses to update his or her knowledge of soil K at the beginning of each period, following five soil test schedules j defined as follows:
$$\begin{equation}
j = 1,\,\,\,\,\,\lambda = \left\{ {\begin{array}{l@{\quad}l} 1 & \hbox{\it if}\quad t = 1,...,10\\ 0 & otherwise \end{array}} \right.,
\end{equation}$$
$$\begin{equation}
j = 2,\,\,\,\,\,\lambda = \left\{ {\begin{array}{l@{\quad}l} 1& \hbox{\it if}\quad t = 1,\,3,\,5,\,7,\,9\,\\ 0 & otherwise \end{array}} \right.,
\end{equation}$$
$$\begin{equation}
j = 3,\,\,\,\,\,\lambda = \left\{ {\begin{array}{l@{\quad}l} 1 & \hbox{\it if}\quad t = 1,\,4,\,7,\,10 \\ 0 & otherwise \end{array}} \right.,
\end{equation}$$
$$\begin{equation}
j = 4,\,\,\,\,\,\lambda = \left\{ {\begin{array}{l@{\quad}l} 1 & \hbox{\it if}\quad t = 1,\,5,\,9\\ 0 & otherwise \end{array}} \right.,
\end{equation}$$
$$\begin{equation}
j = 5,\,\,\,\,\,\lambda = \left\{ {\begin{array}{l@{\quad}l} 1 & \hbox{\it if}\quad t = 1,\,6\,\\ 0 & otherwise \end{array}} \right..
\end{equation}$$
When a producer did not test soil in a given year, the producer's knowledge of carryover K(Ct,j ) was assumed to be the actual carryover K obtained from the most recent soil test (Qt- 1,j ). However, when a producer did soil test, the producer's knowledge of carryover K(Ct,j ) was still assumed to be the actual carryover K obtained from the most recent soil test (Qt,j ). Table 1 shows how the producer's knowledge of K carryover is updated with actual soil test information by temporal frequency.
Table 1. Producer's Knowledge of K Carryover by Temporal Frequency of Soil Testing and Year
Note: This table shows how the constraint C t, j = λ j Q t, j + (1 − λ j )Q t − 1, j is updated in equation (1).
When maximizing NPV, the economic optimality principle of marginal value product (MVP) equals marginal factor cost (MFC) is complicated by intertemporal factors such as the time value of money (opportunity cost) and fertilizer carryover (Kennedy, Reference Kennedy1986; Kennedy et al., Reference Kennedy, Whan, Jackson and Dillon1973). Therefore, a dynamic optimization technique was required to determine optimal total K levels in each period when using soil test information. This study extends Kennedy's (Reference Kennedy1986) dynamic programing framework to determine optimal total available K levels in each period:
$$\begin{equation}
\begin{array}{@{}l@{}} {V_{t,j}}\{ {C_{t,j}}\} = \mathop {{\rm{max}}}\limits_{{A_{t,j}}} [N{R_{t,j}} + \delta {V_{t + 1,j}}\{ {C_{t + 1,j}}\} ]\\ {\rm{Subject\; to{:}}}\\[3pt] {A_{t,j}},{C_{t,j}} \ge 0\\ {C_{t + 1,j}} = {a_0} + {a_1}({A_{t,j}} + {C_{t,j}})\\ {C_{t,j}} = {\lambda _j}{Q_{t,j}} + (1 - {\lambda _j}){Q_{t - 1,j}}\\ {V_{T + 1,j}}\{ {C_{T + 1,j}}\} = 0\\ {Q_0} = {\rm{ given,}} \end{array}
\end{equation}$$
where Vt,j {Ct,j } is the present value of net returns ($/ac.) from applying the profit-maximizing K application in year t; and VT+ 1,j {CT +1,j } = 0 is the terminal condition stating that the producer does not receive any economic value from the available K remaining in the soil after the last period of the planning horizon because the producer will not be able to utilize these available soil K levels (Chiang, Reference Chiang1992).
The optimal single period applied K level was determined using Bellman's (Reference Bellman1957) recursive equation. The optimality conditions were solved by differentiating equation (8) with respect to the decision variable At,j as follows:
$$\begin{equation}
\frac{{\partial {V_{t,j}}}}{{\partial {A_{t,j}}}} = \delta p_t^c\frac{{\partial {y_{t,j}}}}{{\partial {A_{t,j}}}} - p_t^K + \delta \frac{{d{V_{t + 1,j}}}}{{d{C_{t + 1,j}}}}{a_1} = 0,
\end{equation}$$
which can be rearranged as
$$\begin{equation}
\delta p_t^c\frac{{\partial {y_{t,j}}}}{{\partial {A_{t,j}}}} = p_t^K - \delta \frac{{d{V_{t + 1,j}}}}{{d{C_{t + 1,j}}}}{a_1}.
\end{equation}$$
Using the envelope theorem (Léonard and Van Long, Reference Léonard and Van Long1992), differentiating equation (8) with respect to the state variable Ct,j gives the following:
$$\begin{equation}
\frac{{\partial {V_{t,j}}}}{{\partial {C_{t,j}}}} = \delta p_t^c\frac{{\partial {y_{t,j}}}}{{\partial {A_{t,j}}}} + \delta \frac{{d{V_{t + 1,j}}}}{{d{C_{t + 1,j}}}}{a_1}.
\end{equation}$$
Substituting equation (10) into equation (11) and simplifying gives
$$\begin{equation}
\frac{{\partial {V_{t,j}}}}{{\partial {C_{t,j}}}} = p_t^K,
\end{equation}$$
which indicates soil carryover at the beginning of year t is valued at the price of K in year t. This result can be updated to year t + 1:
$$\begin{equation}
\frac{{\partial {V_{t + 1,j}}}}{{\partial {C_{t + 1,j}}}} = p_{t + 1}^K,
\end{equation}$$
which can be substituted into equation (10) to obtain the optimality condition for single period K application rate when using soil test information:
$$\begin{equation}
\delta p_t^c\frac{{\partial {y_{t,j}}}}{{\partial {A_{t,j}}}} = p_t^K - \delta p_{t + 1}^K{a_1}.
\end{equation}$$
Equation (14) indicates the current period optimal total available K level occurs where the MVP (left-hand side of equation 14) was equal to the MFC less the discounted savings associated with K carried over to the next year (right-hand side of equation 14).
2.2. K Carryover Function
Soil K carryover was estimated as a linear function of total K available (i.e., applied K and actual carryover K), which is a commonly used functional form (Harper et al., Reference Harper, Lambert, Larson and Gwathmey2012; Jomini et al., Reference Jomini, Deuson, Lowenberg-DeBoer and Bationo1991; Lanzer and Paris, Reference Lanzer and Paris1981; Segarra et al., Reference Segarra, Ethridge, Deussen and Onken1989). Parameter estimates for the carryover function were obtained using the actual measured total K available. However, depending on the producer's temporal frequency of soil testing, the producer's knowledge of carryover K in the dynamic programming model was updated in each time period with the actual soil K level or the soil K level from the previous soil test. A year random effect was included in the intercept:
$$\begin{equation}
{C_{t + 1,i}} = {a_0} + {a_1}({A_{t,i}} + {Q_{t,i}}) + {\tau _t} + \,{u_{t,i}},
\end{equation}$$
where τt ~ N(0, στ 2) is a random effect capturing the variation in carryover levels across years, and ut,i ~ N(0, uτ 2) is a random error term capturing the variation in carryover levels because of unobserved factors in plot i. The two error terms were assumed to be independent. The intercept, a 0, represents some constant amount of available K that remains in the soil over the planning horizon; the slope, a 1, is the proportion of total K from the current year readily available to the next crop estimated using observed carryover K. Soil K from previous applications accumulates into current-period soil K levels; thus, the only relevant soil K carryover level is for the current period for each of the schedules. Maximum likelihood parameter estimates for equation (15) were obtained using the MIXED procedure in SAS 9.3 (SAS Institute Inc., 2011).
2.3. Yield Response Function
The selection of a functional form to characterize cotton lint yield response to K is important for determining application rates that maximize NPV (Ackello-Ogutu, Paris, and Williams, Reference Ackello-Ogutu, Paris and Williams1985). Plateau-type response functions, such as the linear response plateau or quadratic plus plateau, have been suggested to be more appropriate for characterizing yield response to fertilizer than polynomial or other nonlinear specifications (Ackello-Ogutu, Paris, and Williams, Reference Ackello-Ogutu, Paris and Williams1985; Bullock and Bullock, Reference Bullock and Bullock1994; Cerrato and Blackmer, Reference Cerrato and Blackmer1990). Plateau-type functional forms assume yield responds to an input in a linear or polynomial manner until it reaches a plateau, beyond which the input no longer limits yield. Tembo et al. (Reference Tembo, Brorsen, Epplin and Tostão2008) extended the linear response plateau by including a normally distributed random effect in the plateau to capture variation in the plateau from stochastic events such as insects, weather, and disease. The linear response stochastic plateau (LRSP) developed by Tembo et al. (Reference Tembo, Brorsen, Epplin and Tostão2008) has been found to be more appropriate than similar deterministic functional forms to model yield response to nutrient applications for several crops and provide more accurate economically optimal nutrient rates (Biermacher et al., Reference Biermacher, Brorsen, Epplin, Solie and Raun2009; Boyer et al., Reference Boyer, Larson, Roberts, McClure, Tyler and Zhou2013; Tumusiime et al., Reference Tumusiime, Brorsen, Mosali, Johnson, Locke and Biermacher2011). Furthermore, Harmon et al. (Reference Harmon, Boyer, Lambert, Larson and Gwathmey2016) found the LRSP to be more appropriate than a deterministic linear plateau function in the dynamic programming framework. Therefore, cotton yield response to K applied and actual measured soil K was estimated using the LRSP function:
$$\begin{equation}
{y_{t,i}} = \min [{\beta _0} + {\beta _1}({A_{t,i}} + {Q_{t,i}}),\,\mu + {v_t}] + {w_t} + \,{\varepsilon _{t,i}},
\end{equation}$$
where β 0 and β 1 are the yield response parameters estimated using observed yields; µ is the expected plateau yield parameter (lb./ac.); vt ~ N(0, σv 2) is a normally distributed plateau random effect; wt ~ N(0, σw 2) is the intercept random effect isolating variation in yields from year to year; and εt,i ~ N(0, σε 2) is the random error term accounting for variation in yields from unexplained factors. Independence is assumed across the three random effects. Maximum likelihood parameter estimates for equation (16) were obtained using the NLMIXED procedure in SAS 9.3 (SAS Institute Inc., 2011).
To solve for the optimal applied K rate, the optimality condition (equation 14) was updated with the first-order condition for the LRSP yield response function with respect to applied K, defined as
$$\begin{equation}
\delta p_t^c[{\beta _1}(1 - \Phi )] = p_t^K - \delta p_{t + 1}^K{a_1},
\end{equation}$$
where Φ = Φ{[β 0 + β 1(At,j + Ct,j ) – μ]/σv } is the standard normal cumulative distribution function (Tembo et al., Reference Tembo, Brorsen, Epplin and Tostão2008). Equation (17) is rearranged to obtain the optimal K application rate for period t:
$$\begin{equation}
A_{t,j}^* = \frac{{\left[ {{\Phi ^{ - 1}}\left( {\left. {1 - \frac{{p_t^K - \delta p_{t + 1}^K{a_1}}}{{\delta p_t^c{\beta _1}}}} \right)} \right.} \right]{\sigma _v} + \mu - {\beta _0}}}{{{\beta _1}}} - {C_{t,j}},
\end{equation}$$
which can be simplified to
$$\begin{equation}
A_{t,j}^* = \frac{{{z_m}{\sigma _v} + \mu - {\beta _0}}}{{{\beta _1}}} - {C_{t,j}},
\end{equation}$$
where zm is a standard normal z-score,Footnote 1 and m is the observed probability in the right-hand tail of the standard normal distribution (Tembo et al., Reference Tembo, Brorsen, Epplin and Tostão2008). Thus, the optimal application decision depends on the deviation from the expected plateau yield (zmσv ), which is conditioned on the ratio of per unit K cost and cotton lint price (Tembo et al., Reference Tembo, Brorsen, Epplin and Tostão2008). Additionally, the optimal application rate will depend on the temporal frequency of soil testing and the accuracy of the producer's knowledge of their soil K levels. For the producer who soil tests in every year (j = 1), the producer's knowledge of carryover K is equal to the actual K carryover in each year. However, for a producer who soil tests every other year (j = 2), the producer's knowledge of carryover K is equal to the actual soil K levels in years when a soil test occurs, but the producer's knowledge of carryover K may be higher or lower than the actual K carryover level in periods when the producer does not update soil testing information. Therefore, the producer who updates soil K levels less frequently may over- or underapply K depending on the variability of soil K carryover levels between years. Given that yield responds to the amount of applied and carryover K, the producer's decision of which soil testing schedule to follow will affect yield and the subsequent returns achieved in each period. The NPV for the different soil testing schedules indicates how often a producer needs to update his or her information on soil K to maximize NPV.
2.4. Monte Carlo Simulation
A Monte Carlo simulation was used to introduce uncertainty into the dynamic programming model. One thousand iterations of a 10-year planning period were simulated to generate output distributions of NPV for each of the five scenarios. The prices of K and cotton lint yield, as well as the yield response and carryover coefficients, were assumed to be stochastic, providing an ex ante analysis of the NPV of returns to K for each scenario. Figure 1 summarizes the general process used to solve the dynamic programming model. Shaded boxes correspond with stochastic parameters in the model.
Figure 1. Flow Chart of the Dynamic Programming Model and Simulation Process of Solving for Optimal K Rates
Uncertainty surrounding the prices of cotton lint and K were introduced into the model by bootstrapping the observed real average annual prices of cotton lint and K for each period of the 10-year planning horizon. To introduce uncertainty in the expected yield response, the yield response coefficients were simulated as multivariate normal (MVN) random variables:
$$\begin{eqnarray}
\left[{\begin{array}{@{}*{1}{c}@{}} {\beta _0^*}\\ {\beta _1^*}\\ {\mu _{}^*}\\ {\sigma _v^{2\;*}} \end{array}} \right] \sim MVN\left( {\left[ {\begin{array}{@{}*{1}{c}@{}} {\beta _0^*}\\ {\beta _1^*}\\ \mu \\ {\sigma _v^2} \end{array}} \right]} ,\left[ {\begin{array}{l@{\;\;\,}l@{\;\;\,}l@{\;\;\,}l} {\sigma _{{\beta _0}}^2}&{{\rho _{{\beta _0},{\beta _1}}}{\sigma _{{\beta _0}}}{\sigma _{{\beta _1}}}}&{{\rho _{{\beta _0},\mu }}{\sigma _{{\beta _0}}}{\sigma _\mu }}&{{\rho _{{\beta _0},\sigma _v^2}}{\sigma _{{\beta _0}}}{\sigma _{\sigma _v^2}}}\\ {{\rho _{{\beta _1},{\beta _0}}}{\sigma _{{\beta _1}}}{\sigma _{{\beta _0}}}}&{\sigma _{{\beta _1}}^2}&{{\rho _{{\beta _1},\mu }}{\sigma _{{\beta _1}}}{\sigma _\mu }}&{{\rho _{{\beta _1},\sigma _v^2}}{\sigma _{{\beta _1}}}{\sigma _{\sigma _v^2}}}\\ {{\rho _{\mu ,{\beta _0}}}{\sigma _\mu }{\sigma _{{\beta _0}}}}&{{\rho _{\mu ,{\beta _1}}}{\sigma _\mu }{\sigma _{{\beta _1}}}}&{\sigma _\mu ^2}&{{\rho _{\mu ,\sigma _v^2}}{\sigma _\mu }{\sigma _{\sigma _v^2}}}\\ {{\rho _{\sigma _v^2,{\beta _0}}}{\sigma _{\sigma _v^2}}{\sigma _{{\beta _0}}}}&{{\rho _{\sigma _v^2,{\beta _1}}}{\sigma _{\sigma _v^2}}{\sigma _{{\beta _1}}}}&{{\rho _{\sigma _v^2,\mu }}{\sigma _{\sigma _v^2}}{\sigma _\mu }}&{\sigma _{\sigma _v^2}^2} \end{array}} \right] \right) , \hspace*{-8pt} \nonumber\\
\end{eqnarray}$$
where the mean of the distribution was a vector of the estimated coefficients for the yield response function (equation 16); the variance of the distribution was a four-by-four matrix of the robust covariance estimator of the parameter estimates, where ρ is the correlation coefficient; and the asterisk denotes a randomly drawn coefficient for the simulation (Cuvaca et al., Reference Cuvaca, Lambert, Walker, Marake and Eash2015). The preplanting carryover levels after the initial year were estimated by the linear carryover function (equation 15), where the carryover coefficients followed a MVN distribution:
$$\begin{equation}
\left[ {\left. {\begin{array}{@{}*{1}{c}@{}} {a_{_0}^*}\\ {a_{_1}^*} \end{array}} \right]} \right.\sim MVN\left( {\left[ {\left. {\begin{array}{@{}*{1}{c}@{}} {{a_0}}\\ {{a_1}} \end{array}} \right]} \right.} \right.,\left[ {\left. {\left. {\begin{array}{l@{\quad}l} {\sigma _{{a_0}}^2}&{{\rho _{{a_0},{a_1}}}{\sigma _{{a_0}}}{\sigma _{{a_1}}}}\\ {{\rho _{{a_1},{a_0}}}{\sigma _{{a_1}}}{\sigma _{{a_0}}}}&{\sigma _{{a_1}}^2} \end{array}} \right]} \right)} \right..
\end{equation}$$
For each iteration, new coefficients and prices were randomly sampled to determine, ex ante, the total available K, yield, and NPV.
Uncertainty surrounding the initial carryover level (Q 0) was introduced into the model by bootstrapping the observed carryover levels. In year 1, prices of K and cotton were randomly drawn along with parameter estimates for the carryover and yield response function. The producer was assumed to soil test in year 1, so the preplanting carryover level was determined using the initial carryover level. The yield and carryover parameter estimates, prices of K and cotton, and the estimated preplanting K carryover level were substituted into equation (18) to obtain the optimal application rate in year 1. Subsequent yield and single period net returns were calculated at the optimal application rate and the carryover level. In time period 2, new prices of K and cotton were randomly drawn. After the first period of production, the producer's decision to soil test in each period was determined by the soil testing temporal frequency he or she follows. If the producer tested soil K levels, the optimal K rate was determined using the actual carryover K rate in time period 2. If the producer did not soil test, the producer's knowledge of carryover K was used to determine the optimal K application rate. The subsequent yield and NPV were determined using the application rate and actual carryover levels, and the soil testing decision process was repeated for the remaining years of the planning period for each soil testing schedule. Therefore, after the first year of production, K applications, yield, and NPV for years t = 2,. . .,10 were influenced by the temporal frequency with which producers updated soil test information.
For each scenario, output distributions were generated for the annual and 10-year average applied K, carryover K, lint yield, and NPV for each scenario. The Monte Carlo simulation was conducted using @Risk (Palisade Corporation, 2014). The expected NPVs of each scenario were compared to determine the soil testing temporal frequency that provided the greatest NPV of returns to K.
3. Data
Data on cotton yield response and soil K fertility levels were collected from a 9-year field study (2000 to 2008) conducted at the University of Tennessee, West Tennessee Research and Education Center at Jackson (35.63° N, 88.85° W). The soil type was Loring-Calloway silt loam (thermic Oxyaquic Fragiudal and thermic Typic Fragiaqualf). The plots were not tilled. Each year, K fertilizer (muriate of potash, 0–0–60) was broadcast by hand to individual plots prior to planting at rates of 0, 25, 50, 75, 100, 125, and 149 lb./ac. of elemental K. These treatments were applied to the same plots each year, beginning 5 years prior to the start of the study (2000) through the last year (2008). Plots were arranged in a randomized complete block design, with five or six replications of the fertilizer treatments.
Cotton was planted between April 30 and May 15 of each year using a four-row John Deere MaxEmerge planter. From 2000 to 2002, the cultivar ‘PM1218BG/RR’ was planted on all plots. From 2003 to 2008, two contrasting cultivars were planted in a factorial arrangement relative to the K-fertility plots. The cultivars ‘PM1218BG/RR’ and DP555BG/RR’ were planted from 2003 to 2005, the cultivars ‘FM960BR’ and DP555BG/RR’ were planted from 2006 to 2007, and the cultivars ‘ST455B2RF’ and ‘ST5327B2RF’ were planted in 2008.Footnote 2 Plots were 30 by 12 feet, each containing four rows spaced 38 inches apart. Shortly before or after planting each year, nitrogen fertilizer (ammonium nitrate, 34–0–0) was uniformly drop-spread to all plots at a rate of 80 lb./ac. University of Tennessee Extension Service recommendations were followed for lime and phosphorus applications (Savoy and Joines, Reference Savoy and Joines2001). Supplemental irrigation was used during dry spells in all years except 2002 and 2003. Monthly growing season rainfall for Jackson, Tennessee, is summarized in Table 2 (National Oceanic and Atmospheric Administration, National Climatic Data Center, 2016). All other production practices followed the University of Tennessee Agricultural Extension Service (2001) guidelines for cotton production.
Table 2. Total Monthly Precipitation Levels for the Growing Season of Upland Cotton in Jackson, Tennessee, 2000–2008
Source: National Oceanic and Atmospheric Administration, National Climatic Data Center (2016).
Seed cotton was harvested from the two interior rows of each plot twice each year using a modified John Deere 9930 spindle picker. First harvest occurred from September 7 to October 8, with a second harvest occurring 14 to 28 days later. Lint yields were calculated using seed cotton weights, gin turnouts, and plot areas harvested. Yield response functions were estimated using observed lint yields from 2000 to 2008. Average annual lint yields by K rate are displayed in Table 3. Yields may have increased over time because of improved biotechnology from different cultivars. Therefore, cotton lint yields were tested with a deterministic quadratic time response function (Just and Weninger, Reference Just and Weninger1999). Similar to cotton yields in Oklahoma (Boyer, Brorsen, and Tumusiime, Reference Boyer, Brorsen and Tumusiime2015), a time trend was not present.
Table 3. Average Annual Cotton Lint Yield and Postharvest K Carryover Level by K Application Rate in Jackson, Tennessee, 2000–2008
a Preplanting K carryover levels were not measured for the year 2000.
Within 6 weeks after harvest in each year, soil samples were collected from all plots at the 0- to 6-inch depth using the Mehlich I extraction method (Howard et al., Reference Howard, Essington, Hayes and Percell2001). The samples were tested at the University of Tennessee Soil and Forage Test Laboratory in Nashville, Tennessee. Preplanting soil test levels from 2001 to 2009 were used to estimate the carryover function. The average soil test level for the experiment was characterized by the medium soil fertility range (Savoy and Joines, Reference Savoy and Joines2001) (Table 2). Soil test levels were corrected for heteroskedasticity across years.
Average annual cotton lint and elemental K prices ($/lb.) from 1994 to 2013 were used to determine the K fertilization rates that maximized NPV over a 10-year planning horizon. The Federal Reserve implicit price deflator (Federal Reserve Bank of St. Louis, 2016) was used to adjust nominal prices to reflect real prices in 2013. From 1994 to 2013, real average annual cotton prices varied from $0.38 to $1.07/lb., and real average annual elemental K prices varied from $0.20 to $0.91/lb. (U.S. Department of Agriculture, Economic Research Service, 2013, 2014). Real cotton and K prices were not correlated over time. The real cost of soil testing included the cost of obtaining the soil sample and the chemical analysis. The cost of obtaining the soil sample was $7.27/ac. per year, which was based on the University of Tennessee Custom Rate Survey (Reference BowlingBowling, 2013). The cost of the chemical analysis was $0.70/ac. per year, which assumes a producer soil tests on a 10-acre grid, following University of Tennessee recommendations for soil testing (Savoy and Joines, 2013). A 5% discount rate was used to represent the opportunity cost of land in cotton production, similar to previous dynamic programming literature (Harper et al., Reference Harper, Lambert, Larson and Gwathmey2012; Kennedy et al., Reference Kennedy, Whan, Jackson and Dillon1973; Park et al., Reference Park, Stocker, Hattey and Turner2007; Segarra et al., Reference Segarra, Ethridge, Deussen and Onken1989; Watkins, Lu, and Huang, Reference Watkins, Lu and Huang1998).
4. Results
4.1. Yield Response and Carryover
The parameter estimates for the yield response and carryover functions are presented in Table 4. The intercept of the LRSP function was insignificant and negative, indicating that the intercept was not statistically different than zero. Total available K was always observed to be greater than zero (Table 3); thus, a negative yield would be outside the range of the data. Similarly, Watkins, Lu, and Huang (Reference Watkins, Lu and Huang1998) and Stauber, Burt, and Linse (Reference Stauber, Burt and Linse1975) found insignificant negative estimated intercepts in their yield response to nitrogen when carryover was considered. The remaining parameter estimates had the expected positive signs and were significant at the 1% level. Estimated yields from the LRSP function were plotted against observed yield response to total available K in Figure 2.
Table 4. Parameter Estimates for the Linear Response Stochastic Plateau Yield Response to Total Available K and the Linear Carryover Function
a Single, double, and triple asterisks (*, **, ***) represent significance at the 10%, 5%, and 1% levels, respectively.
b Standard errors are in parentheses.
c Carryover data were corrected for heteroskedasticity.
d Units are reported in lb./ac.
Figure 2. Visualizing the Fit of the Linear Response Stochastic Plateau (LRSP) Functional Form on the Observed Cotton Lint Yield Data
The estimated K carryover function had positive estimates for the intercept and slope (P ≤ 0.01). The intercept indicated that 25.46 lb. K/ac. of soil K did not come from the amount of total K in the previous year but remains available to the plant over the planning period (Lanzer and Paris, Reference Lanzer and Paris1981). The estimated slope indicated that 73% of the total K available in the current period (t) will be carried over to the next period (t + 1). The carryover coefficient was similar to Harper et al.’s (Reference Harper, Lambert, Larson and Gwathmey2012) estimated K carryover coefficient of 0.72.
4.2. Simulation
Monte Carlo simulation results for the annual and 10-year average K application, K carryover, and yield are presented in Table 5. The 10-year average profit-maximizing K application rates for all temporal frequencies of soil testing varied from 29 to 31 lb./ac. per year. However, the range of optimal K application rates across the years varied by temporal frequency of soil testing. The optimal K application in years when a producer did not update his or her soil K information with a soil test was lower than the optimal K application rate when a producer did update his or her information on soil K carryover. These lower applications were offset by applying a higher rate in years when soil testing information was updated to correct for K deficiencies in soil carryover levels. For example, in year 2, the producer that soil tested every year applied a higher K rate to maximize NPV than the producer that soil tested every other year. However, the optimal K application rate for the producer that soil tested every other year was higher in year 3 than the optimal K rate for the producer that soil tested every year, which was necessary to rebuild total available K levels to maximize NPV. The longer the producer waited to update his or her information on soil K carryover with a soil test, the greater the range of optimal K application rate increased.
Table 5. Monte Carlo Simulation Results for the Optimal K Application Rate, Potassium Carryover, and Yield by Simulation Year for a 10-Period Planning Horizon
The 10-year average optimal K carryover levels were 204 lb./ac. per year when a producer updated soil testing information annually. However, when a producer soil tested every other year and every third year, average soil carryover levels were reduced by 2 lb./ac. per year, and average carryover levels were 3 and 4 lb./ac. per year lower than annual soil testing when a producer waited 4 and 5 years to update soil testing information, respectively. However, as temporal frequency of soil testing decreased, the lower bound of the range of annual K carryover level decreased. For instance, when a producer waited 5 years to soil test, the carryover K level dropped from 238 lb./ac. at the beginning of period 1 to 159 lb./ac. when the producer updated soil test information at the beginning of period 6, which was 26 lb./ac. lower than the K carryover level for a producer who soil tests annually. To correct this deficiency, the producer applied a higher rate of K, relative to annual soil testing, to increase the soil K levels. This result shows that when a producer does not know or use soil K carryover information in managing K applications, the optimal strategy was to draw down soil K levels. In turn, the drawing down of soil K levels requires producers to apply higher rates of K in years that soil test information is updated. This finding also resembles what southeastern cotton producers experienced in the late 1980s and early 1990s with K deficiencies (Maples, Thompson, and Varvil, Reference Maples, Thompson and Varvil1988; Mullins, Burmester, and Reeves, Reference Mullins, Burmester and Reeves1997), which reiterates the importance of soil testing for K levels in cotton production in the southeastern United States.
The 10-year average K carryover levels reported in this study would be classified in the medium soil fertility range for each soil testing schedule according to the University of Tennessee Extension Service guidelines (Savoy and Joines, Reference Savoy and Joines2001). However, the University of Tennessee Extension Service recommends K application rates of 90 lb./ac. for medium testing soils in cotton production, which are higher than the optimal rates reported in this study. Moreover, Harper et al. (Reference Harper, Lambert, Larson and Gwathmey2012) reported optimal total K levels of 516 lb./ac., which is higher than the optimal total K levels reported in this study. However, Harper et al. (Reference Harper, Lambert, Larson and Gwathmey2012) used data from a shorter time series on a different soil type and a different yield response function than used in this study, which might explain the differences in results. Our results suggest that recommended K application rates in cotton production based on soil test levels could be decreased.
Annual soil testing produced the highest 10-year average yield of 1,384 lb./ac. per year, and the 10-year average lint yields decreased as the producer's temporal frequency of soil testing decreased. By waiting until every fifth year to soil test, the producer decreased his or her yield by 32 lb./ac. per year relative to the producer that soil tests annually. The longer a producer waited to update soil testing information, the lower the annual yields decreased, which might be attributed to the deficient soil carryover levels limiting yield.
4.3. Optimal Temporal Frequency
The expected NPV increased as temporal frequency increased from soil testing every fifth year ($7,436/ac.) to soil testing every other year ($7,580/ac.) (Figure 3). This indicates that the additional information on soil K carryover had a greater value than the cost of soil testing. The expected NPV increased $83/ac. when a producer went from soil testing every 5 years to every 4 years, $48/ac. when a producer went from soil testing every 4 years to every 3 years, and $13/ac. when a producer went from soil testing every 3 years to every other year (Figure 3). Thus, the findings further indicate that the value of soil testing increased at a decreasing rate when the temporal frequency of soil testing increased. Annual soil testing provided the producer with the most accurate knowledge of soil K carryover variability over time. However, results from this study suggest that the additional value from having the most accurate information on soil K carryover was less than the cost of soil testing to gather that information. When a producer soil tested annually, the expected NPV decreased by $12/ac. from soil testing every other year (Figure 3).
Figure 3. Net Present Value from Applying Optimal K Rates over a 10-Year Planning Period for Five Soil Testing Schedules (different letters indicate a significant difference at the 0.05 level; least significant difference = 58)
Pairwise comparisons were made between the distributions of simulated NPVs for all temporal frequencies of soil testing. The NPV of soil testing every other year was statistically different than the NPVs from soil testing every fourth (P value = 0.08) and fifth year (P value = <0.0001) at the 0.05 level. However, the NPV of soil testing every other year was not statistically different than the NPVs from soil testing annually and every third year. Although producers who use soil testing information to manage K in cotton production can maximize their expected NPV by updating soil testing information every other year, the optimum soil testing temporal frequency may be found soil testing annually to every 3 years.Footnote 3 The conclusion about optimal temporal frequency of soil testing was supported by what Lambert et al. (Reference Lambert, English, Harper, Larkin, Larson, Mooney, Roberts, Velandia and Reeves2014) found in their survey of cotton producers in the southern United States.
5. Conclusion
The objective of this study was to determine the temporal frequency of soil testing for K that maximized NPV of returns to K in cotton production. Cotton lint yield and soil testing data were obtained from a 9-year experiment in Jackson, Tennessee. Cotton lint yield was characterized by an LRSP yield response function, and soil testing information was characterized by a linear carryover function. The producer's objective was to choose the temporal frequency of soil testing and the subsequent annual K application rates that maximized his or her NPV of returns to K over a 10-year planning period. Kennedy's (Reference Kennedy1986) dynamic programming framework was applied, and a Monte Carlo simulation was used to introduce uncertainty into the model.
Previous studies that investigate optimal nutrient application to maximize NPV assume that producers soil test on an annual basis; however, this assumption does not appear to follow producer practices (Lambert et al., Reference Lambert, English, Harper, Larkin, Larson, Mooney, Roberts, Velandia and Reeves2014) or extension recommendations in the Southeast (Kissel and Sonon, Reference Kissel and Sonon2011; Mylavarapu, Reference Mylavarapu1997; Savoy and Joines, 2013). This study expands previous literature by incorporating the temporal frequency of soil testing into a dynamic programming model, and results will assist producers in making better-informed decisions regarding their use of soil testing as a tool to manage K in cotton production.
On average over the 10-year production horizon, the profit-maximizing K application rates for all temporal frequencies of soil testing varied slightly from 29 to 31 lb./ac. per year. The longer producers waited to update their information on soil K carryover with a soil test, the greater the range of optimal K application rate increased. As producers decreased their temporal frequency of soil testing, the lower bounds of the carryover levels and yields decreased because of yield-limiting levels of K. The expected NPV increased at a decreasing rate as temporal frequency increased from soil testing every fifth year to soil testing every other year, indicating that the value of the additional information from soil testing was greater than the cost of soil testing. However, the expected NPV decreased when producers increased their temporal frequency from every other year to annually; thus, the additional value from having the most accurate information on soil K carryover was less than the cost of soil testing to gather that information. The findings of this study are limited to the production region and soil type of the experiment. Production practices, weather, and soil type may produce an optimal soil testing temporal frequency that differs from the results of this study. Further research is needed on the optimal temporal frequency in other cotton-producing regions of the United States.
This study used an LRSP functional form assuming the expected plateau yield is normally distributed, but nonnormality in the expected plateau yield could increase or decrease the optimal K rates (Boyer, Brorsen, and Tumusiime, Reference Boyer, Brorsen and Tumusiime2015). Nonnormality in the expected plateau yield may give producers an incentive to update information more frequently to sustain higher total K levels. Future research could investigate if skewness in the distribution of expected plateau lint yields has an impact on the optimal temporal frequency of soil testing or optimal K application rates for cotton.
Open access
