Author: Torabi, Benjamin; Soltani, Afshin; Galeshi, Serollah; Zeinali, Ebrahim
Date published: July 1, 2011
Nitrogen (N) is often considered as the most important limiting factor for biomass production in natural ecosystems, after water deficit. N fertilization practices can provide a sufficient N supply for plants to achieve the potential yield allowed by the actual climatic conditions (Lemaire et al., 2008). But because of climatic variability, quantities of applied N fertilizer by farmers often are unsuitable for achieving maximum yield. Passioura (2002) showed that for a typical region limited by water availability, poor nitrogen conditions have noticeably limited yield. Therefore, it could be hypothesized that part of the gap between attainable and potential yield might be partially covered by increases in the availability of N, independently of the occurrence of water deficits (Abeledo et al., 2008). It is essential to apply N fertilizers on adequate time and rate. The economically optimum rate of N fertilizer for crops may vary spatially due to variation in soil characteristics and temporally due to the interactions of environmental factors (Mamo et al., 2003; Katsvario et al., 2003; Subedi and Ma, 2007). The rate of applied N fertilizer depends on whether factors influencing soil temperature (Westerman et al., 1999), soil N status, cultivar (Kalra et al., 2007), solar radiation and precipitation (Abeledo et al., 2008). The presence of long-term field experiments to determine optimum N fertilizer and analyze the yield gap under different management and environmental conditions demand great time and cost. Therefore, simulation models for this purpose seem to be helpful. Systems analysis and crop growth simulation have been used in conjunction with surveys of farmer practices, supplemented by measurements of soil properties and crop performance in yield gap analysis (Lobell et al., 2005; Boling et al., 2010a,b; Calvino and Sadras, 2002; Sadras et al., 2002). In recent years, several processes based on dynamic crop simulation models have been developed to predict crop growth, development, and yield using systems approach that integrate knowledge of the underlying processes and interaction of different components of crop production (Boote et al., 1996). The information from long-term simulation experiments is being increasingly used in assessing crop management options and the yield gap analysis by assessing the non-limiting water potential, limiting water potential or limiting nutrient potential yields for a particular region (Aggarwal and Kalra, 1994; Lansigan et al., 1996; Naab et al., 2004; Bhatia et al., 2008). For example, Abeledo et al. (2008) used CERES-Wheat model for assessing yield gap resulted from nitrogen fertilizer. They showed that in the high-yielding potential years, the main restriction for growth was water shortage, and fertilizing only slightly reduced the gap. Conversely, in rainy years characterized by low potential yields and mild water stress, N management may constitute a simple tool for effectively reducing yield gap under rain-fed conditions. Asseng et al. (2008) showed grain yields are often low due to low inputs as results of large uncertainties of rainfall, particularly during the latter part of the season. In addition, they showed soils with higher plant-available water-holding capacity responded more to N applications and thus at high N applications achieved higher yields. Boling et al. (2010a) showed a simulated yield gap of 1.76 t ha-1 (41%) currently exists in rainfed rice farmers' fields. Yield gaps could be substantially reduced by 1.48 t ha-1 (34%) through improved N-management practices. Also in other study, Boling et al. (2010b) analyzed yield-limiting factors (water, N) on rice yields using data from on-farm experiments. Potential, water-limited, and N-limited yields were simulated using the ORYZA2000 crop growth model. Farmers' fields showed large spatial and temporal variation in hydrology and fertilizer doses. The range of yield gap caused by water limitations was 0-28% and that caused by N limitations was 35-63%. Therefore, the yield gap is estimated as the yield difference between optimum and suboptimal factor(s) conditions. Crop simulation models can be used to estimate crop yields under these conditions. However, with respect to the importance of N nutrition in wheat yield production and unawareness of farmers about quantities of N fertilizer application, the aim of the present study was to (i) determine the optimum N fertilizer rate for achieving maximum attainable yield, (ii) estimate the attainable yields under suboptimal N fertilizer conditions and (iii) assess role of N fertilizer on the yield gap in both simulation and farm conditions.
Results and discussion
The studied region is in the north of Iran and it has a sub-humid weather with short and intermittent water stress that can impose wheat yield during the crop growth cycle (Soltani et al., 2000). The long-term mean of annual rainfall is 607 mm. Rainfall is not distributed evenly during the year and usually a relatively wet and cold season (autumn, winter, and early spring) is followed by a dry and warm season (Fig. 1). The amount of annual precipitation was fluctuated from 260-668 mm during a period of 45 years. Monthly maximum and minimum rainfall occurs in March (76.5 mm) and July (20.2 mm), respectively (Fig. 1). Based on the long-term means, the maximum (22.4 MJ m2 day-1) and minimum (8.2 MJ m2 day-1) solar radiation occur in June and December, respectively (Fig. 1). The soil is usually deep and characterized to have relatively high silt and clay and low sand percentage. The soil has a loess nature. The dominant soil textures are loam, silt loam, silty clay loam, and clay loam.
Assessing yield and applied N fertilizer on farms
The amount of yield on farms ranged from 2.5 to 6.3 t ha-1 for 2007 and from 3.4 to 5.9 t ha-1 for 2008 (Fig. 2). The mean yields were 4.7 and 4.6 t ha-1 for 2007 and 2008, respectively. The cumulative distribution functions showed the yield with 50% probability is 4.7 t ha-1 in both years, but the difference in yields between 75 and 25% probability levels was 1.1 and 0.6 t ha-1 for 2007 and 2008, respectively. This indicated that yield variations between farms were more in 2007 than 2008 (Fig. 2). One of the reasons for yield variations between farms is due to the variation of N fertilizer management such as splitting, timing, and amount of N fertilizer application. On-farm assessment showed that the amount of N fertilizer application varied between 45 and 172.5 kg ha-1 in 2007 and between 65 and 175.5 kg ha-1 in 2008 (Fig 3). The mean N fertilizer application was 112 and 94 kg ha-1 for 2007 and 2008, respectively. The cumulative probability distribution showed the mean N fertilizer application with 50% probability is 109 and 87.5 kg ha-1 for 2007 and 2008 respectively. The difference in N fertilizer application between 75 and 25% probability levels was 43 and 26 kg ha-1 in 2007 and 2008, respectively. This indicated N fertilizer application variations between farms were more in 2007 than 2008 (Fig. 3). Variation of N fertilizer application can be explained by the amount of soil's initial N, soil texture, soil moisture, previous crop and N fertilizer availability for farmers (Kalra et al., 2007).
We simulated yield under different N fertilizer application over 1969-2008. The results showed that the mean long-term simulated yield increased between 0 and 200 kg N ha-1. There was no increase in the mean long-term simulated yield over 200 kg N ha-1 application (Table 1). The changes of mean yield over long-term simulation versus N fertilizer were describable using a non-linear, segmented regression model. The segmented model consists of two intersecting lines, a sloping line for the linear increase in yield versus N fertilizer application, and a horizontal line, which determines attainable maximum yield. Mathematically, the segmented model may be expressed as (Soltani et al., 2004, 2005; Fig. 4):
y = a + bx for x < x^sub o^,
y = a + bx^sub o^ for x > x^sub o^
where, y is the mean yield over long-term simulation for each N treatment (t ha-1), x is N fertilizer rate (kg ha-1), a is the intercept with the vertical axis, b is the rate of linear increase in grain yield (t per kg N), xo is the optimum N fertilizer rate for achieving maximum yield. The slope of yield increase was about 0.02 t per kg N fertilizer in hectare. The measure of this slope became zero in 171 kg N ha-1 and above it. In other words, N fertilizer application in more than 171 kg ha-1 did not increase the yield. Therefore, the rate of optimum N fertilizer for achieving maximum attainable yield was estimated 171 kg ha-1. This amount of optimum N fertilizer calculated during the mean long-term simulation can be a suitable criterion for assessing yield gap. Kalra et al. (2007) showed that the higher N fertilizer application results in increasing N-uptake and LAI, albeit these relations had a quadratic behavior. The gain in LAI results in rapid canopy closure. Positive effects of rapid canopy closure are as reducing soil evaporation, inhibition of weeds, the greater CO2 fixation per unit of water transpired, high light interception, and consequently yield increase (Soltani et al., 2001; Soltani and Galeshi, 2007).
Assessing the gap between attainable and maximum attainable yield
To measure the gap between maximum attainable yield under optimum N fertilizer and attainable yield under different N fertilizer, the yield simulated under 171 kg N ha-1 fertilizer. The mean of maximum attainable yield under condition of 171 kg N ha-1 was 5.9 t ha-1 over long-term simulation. The yield gap between maximum attainable yield and attainable yield under N0, N50, N75, N100 and N150 was 3.8, 2.7, 2.1, 1.6 and 0.5 t ha-1, respectively (Table 1). The yield gap was zero in 171 kg N ha-1 and above it. Therefore, in the lower N fertilizer treatments, the addition of N fertilizer increases the attainable yield and consequently yield gap would be decreased (Table 1). The slope of yield gap decrease against N fertilizer application was about 0.02 t ha-1 per kg N fertilizer in hectare (Fig. 5). Yield gap decrease continued until the measure of yield gap became zero in 171 kg N ha-1 and above it (Fig. 5). Grassini et al. (2011) showed with a higher N fertilizer rates, N fertilizer efficiency and yield would be increased. Abeldo et al. (2008) showed that increasing N availability declines yield gap in wheat but the rate of N availability depends on rainfall rate. In high rainfall years with removing water stress, the response of yield to fertilization increased.
Assessing yield gap on farms
The results showed that yield gap between maximum attainable yield and farms' yield ranged from -0.4 to 3.4 t ha-1 for 2007 and 0 to 2.5 t ha-1 for 2008 (Fig. 5). Also, on-farm assessment showed that the amount of N fertilizer application varied between 45 and 172.5 kg ha-1 in 2007 and between 65 and 175.5 kg ha-1 in 2008 (Fig. 3). With respect to the optimum N fertilizer obtained in long-term simulation, it was identified that in 2007 and 2008 nearly all of the farmers (except one farm for each year) used N fertilizers lower than the optimum N fertilizer (Fig. 3). Therefore, a decrease in yield gap between maximum attainable yield and farms' yield was expected with increasing N fertilizer application on farms. Kalra et al. (2007) showed that the application of nitrogen could explain the variation in grain yield. N fertilization generally improves TE through an increase in green area index and radiation use efficiency (Korentager and Berliner, 1987; Heitholt, 1989, cited by Debaeke and Aboudrare, 2004). Application of fertilizer may be increasing the depth of water extraction, or the amount extracted from specific soil layers or both (Brown, 1971, cited by Debaeke and Aboudrare, 2004). In addition, with respect to linear relation between yield gap and variation range of N fertilizer application on farms, it was expected that: (i) range of yield gap in studied farms to be 0-2.7 and 0-2.3 t ha-1 in 2007 and 2008, respectively; (ii) the measure of yield gap in each level of applied N fertilizer on farm to be on regression line. Unlike above points, variation range of yield gap on farms was -0.4 to 3.4 t ha-1 for 2007 and 0 to 2.5 t ha-1 for 2008 and majority of yield gaps were not on linear regression model. These differences between yield gap of simulation and farms are due to different management condition, such as planting date, variety of wheat, N fertilizer splitting, time of N fertilizer application, soil N status , other fertilizers
application, plant density, and soil type of farms (results not shown). Farms yields gap above and below regression line are indicating management condition on farms have been more unsuitable and suitable than it in simulation, respectively. Negative yield gap indicates actual yield on farm has been higher than maximum attainable yield due to the better management condition on farm to simulation. For example, shifting N application from fall to spring or at planting and greater use of split N fertilizer during the growing season, rather than a single large N application, represent options to achieve better congruence (Cassman et al., 2002).
Material and method
The on-farm assessment was carried out in five sites namely Jelin, Mohammad Abad, Karim Abad, Varsan and Nodijeh in the eastern part of the coast of the Caspian Sea in Gorgan, Iran, during 2007 and 2008 years. For assessing crop management, 95 farms with different management levels were selected in these sites. In the first year 50 farms and in the second year 45 farms were selected. Data on crop management such as planting date, cultivar, fertilizer, irrigation, harvest date, and incidence of weeds and disease were collected from these farms through interview and phone. For this purpose, a list of all the agronomic operations, from seedbed preparation to harvesting, was prepared. Then farmers were asked questions to find out how they performed the operations. In this list actual yield harvested by farmer was recorded. To ensure the accuracy of farmers' answers, supervisor engineers of these farms reviewed these completed lists. Data associated with N fertilizer application were used to consider influence of N fertilizer on yield in this study. Cumulative distribution functions were used to evaluate variations of actual yield and applied N fertilizer on the studied farms.
Simulation and statistics analysis
CropSyst model (Stockle et al., 2003) was used for crop growth simulation under different management. Evaluation of CropSyst model under Gorgan condition showed that the model predicted growth and yield of wheat reasonably well (Soltani et al., 2010). Four input data series are needed to run CropSyst: weather statistics, physicochemical soil properties, crop properties, and management practices. Weather statistics, including daily temperature, rain, and radiation, were collected over 1969-2008 at closest meteorological station. Soil, crop, and management practices inputs were collected through farms survey in 2007 and 2008 in this study and also in study of Soltani et al. (2010). However, simulations of estimated attainable yield in different levels of N fertilizer (0, 50, 75 100, 150, 200, 250 and 300 kg N ha-1; shown with N0 to N300, respectively) were made for wheat using CropSyst model over 1969-2008 in Gorgan. In this simulation, each level of N fertilizer was used in three stages: 25 kg N ha-1 before planting and the rest was divided into 60 and 90 days after planting. In all simulations, irrigation was assumed to be implemented in two stages of growth cycle (at flowering an amount of 60 mm and 20 days after flowering an amount 70 mm) that is a common practice in the area. Simulations were performed at a standard sowing date (10 December) and density (450 plant m-2) which is correct for most of the farms in considerable parts of the area. In addition, Soil texture in these simulations was silty clay loam that its depth was 1.2 meter. There were these properties in the majority of studied area soils.
Yield gap analysis
For considering effect of N fertilizer on yield gap, we estimated optimum N fertilizer in which we obtained maximum attainable yield. Optimum N fertilizer obtained using non-linear regression fitted to mean long-term simulated yield in each level of N fertilizer by NLIN procedure of SAS software (Soltani, 2007). The difference between mean maximum attainable yield and mean attainable yields under suboptimal N fertilizers were calculated as yield gap. In addition, the difference between maximum attainable yield and actual yield in farms was calculated in 2007 and 2008. For assessing yield gap and effect of N fertilizer on yield gap, we compared the variation range of yield gap and applied N fertilizer in both simulation and on-farm conditions.
In generally, in this study optimum N fertilizer for achieving maximum attainable yield was estimated about 171 kg ha-1. Farms survey showed that nearly all of the farmers use N fertilizer less than optimum N fertilizer. It concluded that deficiency of applied N fertilizer can be a yield-limiting factor. Thus, farmers should be using higher rate of N fertilizer. In addition, the results showed that in a given level of applied N fertilizer for both simulation and farm conditions, estimated yield gap was different. This difference is related to other agronomic practices. This study suggests considering influence of these practices on yield in the next studies.
Abeledo LG, Savin R, Slafer GA (2008) Wheat productivity in the Mediterranean Ebro Valley: Analyzing the gap between attainable and potential yield with a simulation model. Europ J Agron 28: 541-550.
Aggarwal PK, Kalra N (1994) Analyzing the limitations set by climatic factors, genotype, and water and nitrogen availability on productivity of wheat II. Climatically potential yields and management strategies. Field Crops Res 38: 93-103.
Asseng S, Milroy SP, Poole ML (2008) Systems analysis of wheat production on low water-holding soils in a Mediterranean-type environment: I. Yield potential and quality. Field Crops Res 105: 97-106.
Bhatia VS, Singh P, Wani SP, Chauhan GS, Kesava Rao AVR, Mishra AK, Srinivas K (2008) Analysis of potential yields and yield gaps of rainfed soybean in India using CROPGRO-Soybean model. Agric Forest Meteorol 148: 1252-1265.
Boling AA, Bouman BAM, Tuong TP, Konboon Y, Harnpichitvitaya D (2010a) Yield gap analysis and the effect of nitrogen and water on photoperiod-sensitive Jasmine rice in north-east Thailand. NJAS - Wageningen J of Life Sci doi:10.1016/j.njas.2010.05.001.
Boling AA, Tuong TP, van Keulen H, Bouman BAM, Suganda H, Spiertz JHJ (2010b) Yield gap of rainfed rice in farmers' fields in Central Java, Indonesia. Agric Sys 103: 307-315.
Boote KJ, Jones JW, Pickering NB (1996) Potential uses and limitations of crop models. Crop Sci 88: 704-716.
Brown PL (1971) Water use and soil water depletion by dryland winter wheat as affected by nitrogen fertilization. Agron J 63: 43-46.
Calvino P, Sadras V (2002) On-farm assessment of constraints to wheat yield in the south-eastern Pampas. Field Crops Res 74: 1-11.
Cassman KG, Dobermann A, Walters DT (2002) Agroecosystems, nitrogen-use efficiency, and nitrogen management. Ambio 31: 132-140.
Debaeke P, Aboudrare A (2004) Adaptation of crop management to water-limited environments. Europ J Agron 21: 433-446.
Grassini P, Thorburn J, Burr C, Cassman KG (2011) High-yield irrigated maize in the Western U.S. Corn Belt: I. On-farm yield, yield potential, and impact of agronomic practices. Field Crops Res 120: 142-150.
Heitholt JJ (1989) Water use efficiency and dry matter distribution in nitrogen and water-stressed winter wheat. Agron J 81: 464-469.
Kalra N, Chakraborty D, Kumar PR, Jolly M, Sharma PK (2007) An approach to bridging yield gaps, combining response to water and other resource inputs for wheat in northern India, using research trials and farmers' fields data. Agric Water Manag 93: 54-64.
Katsvario TW, Cox WJ, Van Es HM, Glos M (2003) Spatial yield response of two corn hybrids at two nitrogen levels. Agron J 95: 1012-1022.
Korentager L, Berliner PR (1987) Effects of moisture stress on nitrogen fertilizer response in dryland area. Agron J 80: 977-981.
Lansigan FP, Bouman BAM, Aggarwal PK (1996) Yield gaps in selected rice-producing areas in the Philippines. In: Aggarwal PK, Lansigan FP, Thiyagarajan TM, Rubia EG (Eds) Towards Integration of Models in Rice Research. SAARP Research Proceedings Wegeningen and Los Bafios, 1996.
Lemaire G, Jeuffroy MH, Gastal F (2008) Diagnosis tool for plant and crop N status in vegetative stage: Theory and practices for crop N management. Europ J Agron 28: 614-624.
Lobell DB, Ortiz-Monasterio JI, Asner GP, Naylor RL, Falcon WP (2005) Combining Field Surveys, Remote Sensing, and Regression Trees to Understand Yield Variations in an Irrigated Wheat Landscape. Agron J 97: 241-249.
Mamo M, Malzer GL, Mulla DJ, Huggins DR, Strock J (2003) Spatial and temporal variation in economically optimum nitrogen rate for corn. Agron J 95: 958-964.
Naab JB, Singh P, Boote KJ, Jones JW, Marfo KO (2004) Using CROPGRO-Peanut model to quantify yield gaps in the Guinean Savanna zone of Ghana. Agron J 96: 1231-1242.
Passioura JB (2002) Environmental biology and crop improvement. Funct. Plant Biol 29: 537-546.
Sadras V, Roget D, O'Leary G (2002) On-farm assessment of environmental and management constraints to wheat yield and efficiency in the use of rainfall in the Mallee. Aust J Agric Res 53: 587-598.
Soltani A (2007) Application of SAS in Statistical Analysis. JDM Press, Mashhad, Iran,182p.
Soltani A, Galeshi S (2002) Importance of rapid canopy closure for wheat production in a temperature sub-humid environment: experimentation and simulation. Field Crops Res 77: 17-30.
Soltani A, Galeshi S, Attarbashi MR, Taheri AH (2004) Comparison of two methods for estimating parameters of harvest index increase during seed growth. Field Crops Res 89: 369-378.
Soltani A, Khooie FR, Ghassemi-Golezani K, Moghaddam M (2001) A simulation study of chickpea crop response to limited irrigation in a semiarid environment. Agric Water Manag 49: 225-237.
Soltani A, Khooie FR, Ghassemi-Golezani K, Moghaddam M (2000) Thresholds for chickpea leaf expansion and transpiration response to soil water deficit. Field Crops Res 68: 205-210.
Soltani A, Mahroo-Kashani AH, Dastmalchi A, Maddah V, Zeinali E, Kamkar B (2010) Simulating wheat growth and development using DSSAT, APSIM and CropSyst models under Gorgan and Gonbad conditions (Research Report). Gorgan University of Agricultural Sciences and Natural Resources. 66 Pp.
Soltani A, Torabi B, Zarei H (2005) Modeling crop yield using a modified harvest index-based approach: application in chickpea. Field Crops Res 91: 273-285.
Stockle CO, Donatelli M, Nelson R (2003) CropSyst, a cropping systems simulation model. Europ J Agron 18: 289-307.
Subedi KD, Ma BL (2009) Assessment of some major yield-limiting factors on maize production in a humid temperate environment. Field Crops Res 110: 21-26.
Westerman RL, Raun WL, Johnson GV (1999) Nutrient and water use efficiency. In: Sommers ME (Ed) Handbook of Soil Science, CRC Press, Boca Raton, FL, D, Pp: 175-189.
Benjamin Torabi*, Afshin Soltani, Serollah Galeshi, Ebrahim Zeinali
Department of Agronomy, Gorgan University of Agricultural Sciences, Gorgan, Iran
*Corresponding author : Ben_Torabi@yahoo.com