Research Article

Nitrogen Requirements of Promising Cotton Cultivars in Arid Climate of Multan

Aftab Wajid1, Ashfaq Ahmad1, Muhammad Awais2, Muhammad Habib-ur-Rahman3*, M. Aown Sammar Raza2, Usman Bashir2, Muhammad Naveed Arshad1, Sana Ullah4, Muhammad Irfan5 and Umair Gull1

1Agro-Climatology Labortary Department of Agronomy, University of Agriculture, Faisalabad, 38040, Pakistan; 2Department of Agronomy, University Collage Environmental Sciences, The Islamia University of Bahawalpur, Pakistan; 3Muhammad Nawaz Shareef University of Agriculture, Multan-Pujnab-Pakistan; 4Institute of Soil and Environmental Sciences University of Agriculture, Faisalabad, 38040, Pakistan; 5Department of Crop Physiology, University of Agriculture, Faisalabad, 38040, Pakistan.

Received | October 06, 2017; Accepted | July 15, 2017; Published | August 15, 2017

*Correspondence | Muhammad Habib-ur-Rahman, Muhammad Nawaz Shareef University of Agriculture, Multan-Pujnab-Pakistan; Email: habibagri@hotmail.com

Citation | Wajid, A., A. Ahmad, M. Awais, M. Habib-ur-Rahman, A. Sammar, U. Bashir, M.N. Arshad, Sana Ullah, M. Irfan and U. Gull. 2017. Nitrogen Requirements of Promising Cotton Cultivars in Arid Climate of Multan. Sarhad Journal of Agriculture, 33(3): 397-405.

DOI | http://dx.doi.org/10.17582/journal.sja/2017/33.3.397.405

Keywords | Cotton cultivars adaptability, LAI, LAD, TDM, Seed cotton yield empirical model

Introduction

Cotton, being a cash crop, plays an important role in the development of Pakistan’s economy. It accounts for 1.5 and 7.1% in GDP and agriculture value addition, respectively. The textile industry of Pakistan earned a huge amount of foreign exchange (10.22 billion US$) during July to March of 2014-15. Cotton crop was planted on an area of 2961 thousand hectares during 2014-15 revealing a 5.5% increase as compared to area (2806 thousand hectares) of 2013-14. Similarly cotton production (13.983 million bales) was also 9.5% more as compared to last years (12.769 million bales) (GoP, 2014-15).

Nitrogen (N) plays an important role in cotton productivity (Ali et al., 2003) and its deficiency is one of the main yields limiting factor for cotton cultivars (Milroy and Bange, 2004; Nicholos et al., 2004). Over application of N can lead to excessive vegetative growth, delay in maturity and harvest, promotion of boll shedding, reduction of fiber quality and yield, and increased water pollution (Boquet and Breitenbeck, 2000). Zhao and Oosterhuis (2000) noted a reduction in leaf chlorophyll contents, photosynthetic rate and finally leaf area of cotton crop with N deficiency after vegetative growth. N deficiency promotes the fruit abscission and decrease crop productivity. Excessive N application also reduces crop performance (Howard et al., 2001). Studies show that older cultivars have much lower yield and also poor response to N fertilization than modern cultivars. Aslam et al. (2013) reported that cotton cultivar CIM-506 produced 2568 kg ha-1 of seed cotton with the application of 200 kg N ha-1. Dar et al. (2004) noted that application of 200 kg N ha-1 to cotton cultivar CIM-499 enhanced seed cotton yield (2978 kg ha-1) over all other nitrogen treatments (50, 100, 150 and 200 kg ha-1). The selection of crop variety according to the prevailing climatic conditions is of primary importance (Hussain et al., 2007). Bange et al. (2003) concluded that yield losses of different cultivars were due to reduced boll number associated with less dry matter production. Cotton cultivars differed significantly for total dry matter (TDM) production (Bange et al., 2003) number of bolls per plant (Anwar et al., 2002; Rahman et al., 2016) susceptibility for bollworm (Lisheng, 2005) and economic yield (Hofs et al., 2006; Ehsaan et al., 2008, Wajid et al., 2014). The present study was conducted with the hypothesis that cotton cultivars differ in their response to different nitrogen fertilizer application. The objective of the study was to identify optimum nitrogen level and the best cultivar x nitrogen combination that would lead to maximum economic cotton yield in the cotton growing area of Cotton-Wheat cropping system of Punjab.

Materials and Methods

An experiment was conducted to study the influence of various N rates on growth and yield of cotton cultivars at the Central Cotton Research Institute, Multan during cotton growing season in 2014. The treatments that include cotton cultivars (CIM-496 and NIAB-111) and nitrogen rates (50, 100, 150 and 200 kg ha-1) were randomly allocated to the main and sub plots, respectively utilizing randomized complete block design (RCBD) with split plot arrangement. Used cultivars have potential to survive better in arid climatic conditions than others. Performance of these cultivars were found reasonable good under semi-arid to arid climate so these cultivars were selected for this study. Each subplot has a length and width of 10 m and 3 m, respectively. The seeds of tested verities were sown on May 16, 2014 utilizing hand drill and maintained 75 cm row spacing while plant spacing was maintained at 25 cm by thinning (20 days after sowing). Phosphorus (Diammonium phosphate) and potassium (potassium sulphate) were applied at the rate of 60 and 50 kg ha-1 in all the plots. Complete amount of phosphorus and 1/3rd nitrogen (Urea) was applied at seeding, 1/3rd N at flowering and remaining 1/3rd at boll formation. The crop was irrigated 9 times at different stages (monopodial branches, squaring, flowering, boll formation and boll opening initiation phases) and insecticides (actara 25WP @ 24 g, mospilan 20SP @ 125 mg, polo 500SC @ 250 ml, karate 2.5EC @ 400 ml and talstar 10EC @ 250 ml) were used to control the sucking insects and bollworms. Final yield was determined from a net plot size of 1.5 m × 6.0 m and rest area of the plot was utilized to collect sample for crop growth. Randomly three plants from each plot were harvested with an interval of 30 days. After recoding weight of all samples, these samples were placed in an oven for total dry matter (TDM) calculation. Leaf area index (LAI) (Leaf area / Land area) and Leaf area duration (LAD) = [(LAI1 + LAI2) /2 × {(T2 – T1)] were calculated as proposed by Hunt (1978). At maturity, an area of 1.5 m x 6 m which had 48 plants was harvested for final yield and various yield components. Sigma Plot was used to develop empirical model for the prediction of seed cotton yield on the basis of nitrogen application in the field.

Weather: Mean monthly weather data for the crop growth period is presented in Table 1. The maximum average air temperature was recorded in the month of June (34.7°C) and minimum average temperature was recorded in the month of November (19.7°C). The average relative humidity increased from May to December (maximum 87.3%). The maximum average rainfall was recorded in the month of July (134.3 mm) while no rainfall was occurred in the months of June, August, October and November. The total rainfall throughout the cotton growing season was 155.3 mm.

Table 1: Meteorological data for cropping season 2014.

LAI has prime significance in increasing the crop yield (Awais et al., 2013). Development of canopy during the growing season in response to cultivars and nitrogen applications (Figure 1 a and b). LAI increased rapidly, in all the treatments, up to 140 days after sowing (DAS) and then decreased gradually with the passage of time due to leaf senescence. Cultivar CIM-496 developed its canopy at much faster rate as compared to the other cultivar (NIAB-111). At final harvest CIM-496 gained more LAI (2.08) than NIAB-111 (1.80). Wajid et al. (2010) reported that cultivar NIAB-111 produced more LAI (3.5) as compared to CIM-496 (3.2) which did not confirm the findings of our study. Among the nitrogen levels, 200 kg N ha-1 produced maximum LAI (4.98) at 140 DAS. Decline in LAI was more obvious at lower N rates than high N rates due to senescence of leaves during the early growth season. Higher nitrogen levels increased LAI due to increase in leaf expansion. Leaf greenness was increased; radiation absorption was improved while reflectance was decreased at higher use of nitrogen. Nitrogen application boosted the photosynthetic rate, leaf expansion and leaf persistency. As leaf area directly relates to the rate as well as duration of leaf expansion, so LAI finds to be sensitive to N availability to crop plants. Wajid et al. (2010) reported that less LAI (2.8) was observed at 80 kg N ha-1 and value of 3.9 was attained from the application of 185 kg N ha-1.

Total dry matter

Nitrogen application improves crop development and growth (Kumbhar et al., 2008; Awais et al., 2015).

Table 2: Performance of two cotton cultivars under different rates of nitrogen.

The response of TDM to cultivars and nitrogen doses at different harvests during the growth season of the crop is indicated in Figure 2 a and b. TDM continued to increase linearly and gradually throughout the growth period of crop up to maturity. Five harvests for TDM were recorded with an interval of 30 days starting from 50 days after sowing (DAS) to maturity. Among varieties CIM-496 produced more TDM (868.49 g m-2) than NIAB-111 (731.77 g m-2). In case of nitrogen levels, 200 kg N ha-1 produced maximum TDM (1025.50 g m-2) followed by 150 kg N ha-1 which produced TDM of 818.00 g m-2. And lowest TDM (645.60 g m-2) was produced by the 50 kg N ha-1 treatment. Nitrogen takes part in the production of energy rich compound, adenosine triphosphate, (ATP) that control photosynthesis and ultimately assimilate crop biomass (Sawan et al., 2009). Wajid et al. (2010) concluded that each increment of nitrogen application had a positive effect on TDM and it ranged between 874-1308 g m-2.

Data indicates that CIM-496 remains green for longer period (223 days) as compared to NIAB-111 (170 days). Data showed that each increment of nitrogen enhanced the capacity of leaves to remain photosynthetically active (Table 2). At higher rates of N (200 kg ha-1) leaf tissues remained green for more days (248.40 days) as compared to lower rate of nitrogen. Leaf area duration has positive association with total dry matter (R2 = 0.96) while TDM has also positive relationship with seed cotton yield (Figure 3). An increase in leaf area duration of cotton cultivars in response to different rates of nitrogen was also reported by Wajid et al. (2010).

Plant height

Data reveal that two varieties responded differently in attaining plant height. It was noted that variety CIM-496 attained higher plant height (134.61 cm) as compared to NIAB-111 (118.37 cm). Taller plants in case of CIM-496 may be due to better genetic potential and efficiency for N uptake. Plant height depends on N application in cotton (Rochester et al., 2001). Copur (2006) also reported considerable differences in plant height of different cultivars. Each increment of nitrogen application had significant effect on plant height and it ranged from 121 to 132 cm. The tallest plants (131.78 cm) were observed with the 200 kg N ha-1 application. The increase in plant height due to increasing levels of nitrogen application could be due to increased cell division and cell elongation. Our results are supported by earlier findings that plant height was significantly affected by the genetic make up (Wankhade et al., 2002; Jatt et al., 2007), cultivars (Anwar et al., 2002) and nitrogen application rate (Rochester et al., 2001).

Sympodial branches

Higher amount of sympodial branches in a plant is an indication of higher potential of cotton crop for high production of seed cotton because these are considered the boll bearing branches (Hussain et al., 2007). Cotton variety CIM-496 produced more number of sympodial branches per plant (22.0) as compared to NIAB-111 (19.60). These significant differences may be due to high sink capacity and potency for withdrawing the food materials in case of variety CIM-496. Khan et al. (2007) also reported that sympodial branches per plant may vary among different cotton cultivars. Nitrogen levels significantly increased sympodial branches per plant. Maximum number of sympodial branches per plant (22.34) was observed with the application of 200 kg N ha-1. The interaction between cultivars and nitrogen levels was also significant as indicated in Table 2. The cultivar CIM-496 at nitrogen level of 200 kg N ha-1 gave the largest number of sympodial branches per plant (23.97) and was statistically different from all other combinations, while the variety NIAB-111 at nitrogen level of 50 kg N ha-1 produced the least number of sympodial branches per plant. Nitrogen boosted vegetative growth of crop plants. The sympodial branches, being a component of vegetative growth, indicated a positive response to each N increment. A gradual increase in sympodial branches with the subsequent increases in the nitrogen application rates was also reported by Kumbhar et al. (2008).

Bolls per plant

More number of bolls per plant ultimately produced higher seed cotton yield. Variety CIM-496 produced more bolls (30.42) per plant than NIAB-111, which

expressed 25.17 bolls per plant. The difference in number of bolls per plant among verities was due to difference in number of sympodial branches per plant and number of flowers formed. Similar results were revealed from the study of Rahman et al. (2016) for different cotton cultivars shows that genetic potential does matter and contribute differently to number of bolls and ultimately to cotton yield. Difference in number of bolls per plant among different cotton cultivars was also reported by Khan et al. (2007). Maximum number of bolls per plant was produced by higher nitrogen application (Seilsepour and Rashidi, 2011). The same was recorded in our study, however interaction between cultivar and nitrogen was found to have non significant effect on number of bolls per plant. The highest boll number (31.00) was recorded with 200 kg ha-1 N rate which was statistically different from all other nitrogen rates. The minimum boll number (24.83) was recorded with 50 kg N ha-1. The similar results were also reported by several researchers (Ali and El-Sayed, 2001; Kumbhar et al., 2008; Rashidi and Gholami 2011).

Average weight per boll

Both size and weight of seed cotton per boll are directly related to the final yield. Cotton cultivars indicate diverse potential for producing light or heavy weight bolls (Hofs et al., 2006). In case of average boll weight cultivar CIM-496 gave heavy bolls (3.49 g) than the cultivar NIAB-111, which gained 3.02 g weight per boll. Hofs et al. (2006) also described the differences in boll weight among different varieties of cotton. Boll weight increased by increasing application of nitrogen to cotton crop (Rashidi and Gholami, 2011). As far as the results for the effect of nitrogen rates on bolls weight are concerned, more weight per boll (3.48 g) was recorded with the application of 200 kg N ha-1 followed by 150 kg N ha-1. Nitrogen enhanced the photosynthetic rate, so more photosynthates were produced and stored in different plant parts (Evans, 1989). So, nitrogen positively contributed to growth of cotton plant that finally led to higher boll weight. Minimum boll weight (3.04 g) was produced by the 50 kg N ha-1 treatment. Numerous studies have shown the positive effect of nitrogen application on boll weight (Sawan et al., 2006; Saleem et al., 2010; Seilsepour and Rashidi, 2011).

The cultivars showed significant differences for cotton productivity (Ehsan et al., 2008), which is the interaction of entire yield components (Khan et al., 2007). Seed cotton yield was found significantly different between cultivars (Table 2). Cultivar CIM-496 produced more seed cotton yield (2.17 t ha-1) as compared to cultivar NIAB-111 (1.96 t ha-1). Cultivars vary significantly in seed cotton yield (Copur, 2006; Hofs et al., 2006), this may be due to differences in fruiting branches, productive bolls and seed cotton weight per boll. Seed cotton yield boosted gradually with each N fertilizer increment. The considerable higher seed cotton yield (2.42 t ha-1) was achieved from the treatment where 200 kg N ha-1 was applied. Application of 200 kg N ha-1 enhanced yield by 29% over low dose of N (50 kg ha-1). Leaf area duration (days) and fruit bearing branches (sympodial) have direct and positive association with seed cotton yield with reasonably well high values of coefficient of regression 0.91 and 0.88, respectively (Figure 4). Nitrogen promoted biomass by accumulating more photosynthetically active radiation (PAR) on larger canopy promoted by higher nitrogen application compared with lower noes while positive and strong association was found between PAR and total dry mater with R2 of 0.99 (Figure 5). Number of bolls and boll weight have also positive relationship with seed cotton yield as reasonable good association was found with reasonably good values of coefficient of regression, 0.86 and 0.77, respectively (Figure 6). Empirical model was developed on the basis of observed findings of seed cotton yield in relation with nitrogen application, higher coefficient of regression (0.99) and lower standard error (0.0142) was found for this model (Figure 7). This model was used for the prediction of seed cotton yield with the application of nitrogen beyond 200 kg N ha-1 , it is forecasted that seed cotton yield increased but with lower increase in yield for each unit addition in nitrogen beyond 200 kg N ha-1 (Figure 7). Nitrogen application promoted photosynthetic rate, assimilates production and accumulation that ultimately boosted final seed cotton yield (Kumbhar et al., 2008). Role of N fertilize in improving seed cotton yield was also reported by previous studies (Felix et al., 2003; Wajid et al., 2010; Rashidi and Gholami, 2011).

Conclusions

From the above discussion it was indicated that variety CIM-496 produced higher plant height, more number of sympodial branches per plant, higher boll weight, maximum number of bolls per plant and maximum seed cotton yield. Among nitrogen increments, nitrogen application of 200 kg ha-1 was best for attaining maximum seed cotton yield under Multan climatic conditions in Cotton-Wheat cropping system of Punjab.

Authors Contribution

Aftab Wajid and Ashfaq Ahmad designed the experiment and supervised all work. M. Awais, M. Habib ur Rahman, Muhammad Naveed Arshad, M. Irfan and Umair Gull collect the data and analyzed. Usman Bashir, M. Awais, M. Aown Sammar Raza write the article while it was improved by Muhammad Irfan and M. Habib ur Rahman. It was corresponded by Dr. M. Habib ur Rahman.

References

Ali, L., M. Ali and Q. Mohy-ud-din. 2003. Assesment of optimum nitrogen requirement and economics of cotton crop for seed yield. Int. J. Agric. Biol. 5:493-495.

Ali, S.A. and El-Sayed. 2001. Effect of sowing dates and nitrogen levels on growth, earliness and yield of Egyptian cotton cultivar Giza 88. Egypt. J. Agric. Res. 79:221-232.

Aslam, M., A.H. Sanghi, L. Khalid and S. Javed. 2013. Effect of different doses of nitrogen on the yield of different cotton varieties. Asian J. Agric. Biol. 1(4):179-182.

Anwar, A.M., M.I. Gill, D. Muhammad and M.N. Afzal. 2002. Evaluation of cotton varities at different doses of nitrogen fertilizer. The Pak. Cottons. 46:35-41.

Awais, M., A. Wajid, A. Ahmad and A. Bakhsh. 2013. Narrow plant spacing and nitrogen application enhances sunflower (Helianthus annuusL.) productivity. Pak. J. Agric. Sci. 50(4):689-697.

Awais, M., A. Wajid, A. Ahmad, M.F. Saleem, M.U. Bashir, U. Saeed, J. Hussain and M.H. Rahman. 2015. Nitrogen fertilization and narrow plant spacing stimulates sunflower productivity. Turk. J. Field Crops. 20:99-108. https://doi.org/10.17557/.17374

Bange, M.P., S.P. Milroy and P. Thongbai. 2003. Growth analysis of short and long season cotton cultivars. Aust. Cotton Co-operative Res. Centre, Australia. Field Crops Res. 88:129-142.

Boquet, D.J. and G.A. Breitenbec. 2000. Nitrogen rate effect on partitioning of nitrogen and dry matter by cotton. Crop Sci. 40:1685-1593.

Copur, O. 2006. Determination of yield and yield components of some cotton cultivars in semi-arid conditions. Pak. J. Biol. Sci. 9(14):2572-2578. https://doi.org/10.3923/pjbs.2006.2572.2578

Dar, J.S. and B. Khan. 2004. Yield Variations of CIM-499, CIM-511 AND CIM-707 Cotton Varieties as Affected by Different Nitrogen Levels. Pak. J. life Soc. Sci. 3(1):178-181

Ehsan, F., A. Ali, M.A. Nadeem, M. Tahir and A. Majeed. 2008. Comparative yield performance of new varieties of cotton (Gossypium hirsutum L.). Pak. Life Soc. Sci. 6: 1-3.

Evans, J.R. 1989. Photosynthesis and nitrogen relationships in leaves of C 3 plants. Oecologia. 78(1): 9-19. https://doi.org/10.1007/BF00377192

Felix, B., B.A. Robert, R.L. Travis, D.W. Rains, and R.B. Hutmacher. 2003. Response of irrigated Acala and Pima cotton to Nitrogen fertilization, growth, Dry matter partitioning and yield. Agron. J. 95:133-146. https://doi.org/10.2134/agronj2003.0133

Government of Pakistan. 2014-2015. Economic Survey of Pakistan. Ministry of Food, Agriculture and Livestock, Finance Division, Economic Advisor Wing, Islamabad, Pakistan. pp: 26.

Hofs, J.L., B. Hau and D. Marais. 2006. Boll distribution patterns in Bt and non-Bt cotton cultivars: I. Study on commercial irrigated farming systems in South Africa. Field Crops Res. 98:203-209. https://doi.org/10.1016/j.fcr.2006.01.007

Howard, D.D., C.O. Gwathmey, M.E. Essington, R.K. Roberts and M.D. Mullen. 2001. Nitrogen fertilization of no-till cotton on loess-derived soils. Agron. J. 93:157-163. https://doi.org/10.2134/agronj2001.931157x

Hunt, R. 1978. Plant growth analysis. The institute of Biology’s studies in Biology. Edward Arnold (Pub.) Ltd., 96:8-38.

Hussain, M., A. Ahmad and S.I. Zamir. 2007. Evaluation of Agro-Qualitative Characters of five cotton cultivars (Gossypium hirsutum L.) grown under Toba Tek Singh conditions. Pak. J. Agric. Sci. 44:575-580.

Jatt, T., H. Abro, A.S. Larik and Z.A. Soomro. 2007. Performance of different cotton avrities under the climatic conditions of Jamshoro. Pak. J. Bot. 39:2427-2430.

Khan, N.U., H.U. Khan, K. Usman, H. Khan and S. Alam. 2007. Growth and yield related characterstics of improved commercial cotton cultivars. Sarhad J. Agric. 23:261-263.

Kumbhar, A.M., U.A. Buriro, S. Junejo, F.C. Oad, G.H. Jamro, B.A. Kumbhar and S.A. Kumbhar. 2008. Impact of different nitrogen levels on cotton growth, yield and N-uptake planted in legume rotation. Pak. J. Bot. 40:767-778.

Lisheng, B. 2005. Comprehensive evaluation of new varieties of insect resistant hybrid cotton by grey association degree. Acta Agric. Shanghai.21: 12-14.

Milroy, S.P and M.P. Bange. 2004. Nitrogen and light responses of cotton photosynthesis and implications for crop growth. Australian Cotton Co-operative Res. Centre, Narrabri NSW 2390, Australia. Aust. J. Crop Sci. 43:904-913. https://doi.org/10.2135/cropsci2003.9040

Nicholos, S.P., C.E. Snipes and M.A. Jones. 2004. Cotton growth, lint yield, and fiber quality as affected by row spacing and cultivar. J. Cotton Sci. 8:1-12.

Rahman, M. H., A. Ahmad, A. Wajid, M. Hussain, J. Akhtar, G. Hoogenboom. 2016. Estimation of temporal variation resilience in cotton varieties using statistical models. Pak. J. Agric. Sci. 53(4):169-186. https://doi.org/10.21162/PAKJAS/16.4549

Rashidi, M. and M. Gholami. 2011. Nitrogen and boron effects on yield and quality of cotton (Gossypium hirsutum L.). Inter. Res. J. Agric. Sci. Soil Sci. 1:118-125.

Rochester, I.J., M.B. Peoples, N.R. Hulugalle, R.R. Gault and G.A. Constable. 2001. Using legumes to enhance nitrogen fertility and improve soil condition in cotton cropping systems. Field Crops Res. 70:27-41. https://doi.org/10.1016/S0378-4290(00)00151-9

Saleem, M.F., M.F. Bilal, M. Awais, M.Q. Shahid and S.A. Anjum. 2010. Effect of nitrogen on seed cotton yield and fiber qualities of cotton (Gossypium hirsutum L.) cultivars. J. Anim. Plant Sci. 20:23-27.

Sawan, Z.M., A.H. Fahmy and S.E. Yousaf. 2009. Direct and residual effect of nitrogen fertilization, foliar application of potassium and plant growth retardant on Egyptian cotton growth, seed yield, seed viability and seed vigor. Acta Ecological Sinica, 29:116-123. https://doi.org/10.1016/j.chnaes.2009.05.008

Sawan, Z.M., M.H. Mahmoud and A.H.E. Guibali. 2006. Response of yield, yield components and fibre properties of Egyptian cotton (Gossypium barbadense L.) to nitrogen fertilization and foliar-applied potassium and mepiquat chloride. J. Cotton Sci. 10:224-234.

Seilsepour, M. and M. Rashidi. 2011. Effect of Different Application Rates of Nitrogen on Yield and Quality of Cotton (Gossypium hirsutum). American-Eurasian J. Agric. Environ. Sci. 10:366-370.

Wajid, A., A. Ahmad, T. Khaliq, S. Alam, A. Hussain, K. Hussain, W. Naseem, M. Usman and S. Ahmad. 2010. Quantification of growth, yield and radiation use efficiency of promising cotton cultivars at varying nitrogen levels. Pak. J. Bot. 42:1703-1711.

Wajid, A., A. Ahmad, M. Hussain, M.H. Rahman, T. Khaliq, M. Mubeen, F. Rasul, U. Bashir, M. Awais, J. Iqbal, S.R. Sultana and G. Hoogenboom. 2014. Modeling growth, development and seed-cotton yield for varying nitrogen increments and planting dates using DSSAT. Pak. J. Agri. Sci. 51:639-647.

Wankhade, S.T. and A.K. Gobble. 2002. Performance of different cotton (Gossypium hirsutum L.) cultivars at different nitrogen rates. Annals Plant Phys. 1:569-571.

Zhao, D. and D.M. Oosterhuis. 2000. Nitrogen application effect on leaf photosynthesis, nonstructural carbohydrate concentrations and yield of field-grown cotton. Special Report - Arkansas Agric. Exp. Station. 198: 69-71.