Introduction
The world's population is expected to rise to 8.3 billion in 2030 and to 9.3 billion in 2050, despite predictions from the UN that growth will dramatically slow down. Despite this slowdown in the population increase, there is a rising need for food, feed, fuel and fibre (Ahmad et al., Reference Ahmad, Ashfaq, Rasul, Wajid, Ahmad, Khaliq, Nasir, Rasul, Riaz, Ahmad, Ahmad, Baig, Valdivia, Hoogenboom, Rosenzweig, Mutter and Contreras2021). Cotton (Gossypium hirustum L.) is commonly cultivated and then traded as fibre, yarn, fabric or finished commodities (Ballester et al., Reference Ballester, Hornbuckle, Brinkhoff and Quayle2021; Van der Sluijs, Reference Van der Sluijs2022; Ahmad et al., Reference Ahmad, Ahmad, Ahmad, Ahmad, Wajid, Khaliq, Abbas, Wilkerson and Hoogenboom2023). Customers prefer cotton fibre over synthetic fibre because it is a natural fibre that is soft to the touch, absorbs moisture well, and has a variety of uses (Wakelyn et al., Reference Wakelyn, Bertoniere and French2006; Thiry, Reference Thiry2011; Van der Sluijs and Johnson, Reference van der Sluijs and Johnson2011; Ahmad et al., Reference Ahmad, Abbas, Abbas, Fatima, Atique-ur-Rehman, Naz, Younis, Khan, Nasim, Rehman, Ahmad, Rasul, Khan and Hasanuzzaman2017, Reference Ahmad, Ashfaq, Rasul, Wajid, Ahmad, Khaliq, Nasir, Rasul, Riaz, Ahmad, Ahmad, Baig, Valdivia, Hoogenboom, Rosenzweig, Mutter and Contreras2021). Because cotton is cultivated both as a fibre and as an oilseed crop, it has great social and economic importance. It is currently one of the top ten cash crops in the agricultural sector worldwide (de Oliveira Araujo et al., Reference de Oliveira Araujo, Antonio Camacho and Vincensi2013; Khan et al., Reference Khan, Tan, Munsif, Afridi, Shah, Wei, Fahad and Zhou2017; Tariq et al., Reference Tariq, Yasmeen, Ahmad, Hussain, Afzal and Hasanuzzaman2017, Reference Tariq, Afzal, Muhammad, Ahmad, Shahzad, Kiran and Wakeel2018, Reference Tariq, Abbas, Yasmeen, Ahmad, Ahmad and Hasanuzzaman2020, Reference Tariq, Fatima, Iqbal, Nahar, Ahmad and Hasanuzzaman2021) and is a source of livelihood for 24.2 million people (International Cotton Advisory Committee, 2022). During 2021, it was grown on 32.6 million ha with a total production of 25 million tons of lint and 43 million tons of cotton seed for oil (International Cotton Advisory Committee, 2022). More than 70 countries grow cotton, with Asia contributing more than half of the total world production (Ahmad and Raza, Reference Ahmad and Raza2014; Tariq et al., Reference Tariq, Afzal, Muhammad, Ahmad, Shahzad, Kiran and Wakeel2018, Reference Tariq, Khan, Muhammad, Ahmad and Ahmed2022; Ahmad and Hasanuzzaman, Reference Ahmad and Hasanuzzaman2020; Ahmad et al., Reference Ahmad, Ahmad, Ahmad, Ahmad, Wajid, Khaliq, Abbas, Wilkerson and Hoogenboom2023).
Pakistan is a leading cotton-producing country, and cotton is Pakistan's most important cash crop (Tariq et al., Reference Tariq, Afzal, Muhammad, Ahmad, Shahzad, Kiran and Wakeel2018, Reference Tariq, Khan, Muhammad, Ahmad and Ahmed2022; Ahmad et al., Reference Ahmad, Ahmad, Ahmad, Ahmad, Wajid, Khaliq, Abbas, Wilkerson and Hoogenboom2023). Cotton produces both edible oil and yarn (Matloob et al., Reference Matloob, Aslam, Haseeb Ur Rehman, Khaliq, Ahmad, Yasmeen, Ahmad and Hasanuzzaman2020; Munir et al., Reference Munir, Rasul, Ahmad, Sajid, Ayub, Arif, Iqbal, Khan, Fatima, Ahmad, Khan, Ahmad and Hasanuzzaman2020; Rahman et al., Reference Rahman, Ahmad, Ghaffar, Haider, Ahmad, Ahmad, Tariq, Nasim, Rasul, Fahad, Ahmad, Hoogenboom, Ahmad and Hasanuzzaman2020). Exports of lint and other value-added cotton goods account for 60% of Pakistan's overall foreign exchange profits, with cotton accounting for 0.7% of GDP and 2.9% of value-added agriculture (Govt. of Pakistan, 2024). Due to the cotton leaf curl virus and floods in 2021 and 2022, the cotton cropping area decreased to 1.94 million ha with a total production of 8.33 million bales. There are other factors responsible for this decrease, which include climate change, and improper crop nutrition (Rinehardt et al., Reference Rinehardt, Edmisten, Wells and Faircloth2004; Rafi et al., Reference Rafi, Allah, Abdul, Safdar, Muhammad, Javeed and Amjed2015; Tariq et al., Reference Tariq, Afzal, Muhammad, Ahmad, Shahzad, Kiran and Wakeel2018, Reference Tariq, Khan, Muhammad, Ahmad and Ahmed2022; Afzal et al., Reference Afzal, Tariq, Ahmad, Mubeen, Khan, Afzal and Ahmad2019; Raphael et al., Reference Raphael, Echer and Rosolem2019; Mubeen et al., Reference Mubeen, Afzal, Tariq, Ahmad, Muhammad, Shehzad and Yonas2021; Ahmad et al., Reference Ahmad, Ahmad, Ahmad, Ahmad, Wajid, Khaliq, Abbas, Wilkerson and Hoogenboom2023). Despite significant genetic improvement, global yield has been stagnant at around 750 kg lint/ha since 2010 (International Cotton Advisory Committee, 2022).
The most important nutrient for crop productivity is N, and the N fertilizer application recommendations are directly related to yield targets (Snider et al., Reference Snider, Harris, Roberts, Meeks, Chastain, Bange and Virk2021; Ahmad et al., Reference Ahmad, Ahmad, Ahmad, Ahmad, Wajid, Khaliq, Abbas, Wilkerson and Hoogenboom2023). The N availability affects cotton growth, development, physiological processes and ultimately yields. Inadequate N leads to a series of deficiency symptoms that are generally recognized as N stress (Radin and Mauney, Reference Radin, Mauney, Mauney and Stewart1986; Howard et al., Reference Howard, Gwathmey, Essington, Roberts and Mullen2001; Cochran et al., Reference Cochran, Roberts, Larson and Tyler2007; Snyder et al., Reference Snyder, Bruulsema, Jensen and Fixen2009; Chen et al., Reference Chen, Hou, Wu, Liang and Wei2010; Srivastava et al., Reference Srivastava, Panda, Chakraborty and Halder2018; Dhakal et al., Reference Dhakal, Lange, Parajulee and Segarra2019). Low N exposure causes plants to grow more slowly, produce fewer fruiting sites, have lower leaf area and mature earlier (Radin and Mauney, Reference Radin, Mauney, Mauney and Stewart1986; Reddy et al., Reference Reddy, Reddy and Wang1997). The N deficiency also lead to shedding of fruit and thus fewer bolls and a decrease in yield (Boquet and Breitenbeck, Reference Boquet and Breitenbeck2000; Bondada and Oosterhuis, Reference Bondada and Oosterhuis2001).
Previous studies revealed that N application rates to cotton crops resulted in variable N accumulation in different plant parts (Hou et al., Reference Hou, Xiang, Fan, Zhang, Hu, Yan, Guo, Xiao, Li, Cheng and Li2021), accumulation in the root structure (Chen et al., Reference Chen, Liu, Wang, Zhang, Sun, Song, Bai, Lu and Li2020a, Reference Chen, Liu, Zhang, Ge, Shen, Ma, Zhang, Dong, Zhao and Pang2020b), and regulation of auxin (Krouk et al., Reference Krouk, Lacombe, Bielach, Perrine-Walker, Malinska, Mounier, Hoyerova, Tillard, Leon, Ljung, Zazimalova, Benkova, Nacry and Gogon2010), abscisic acid, and salicylic acid (Chen et al., Reference Chen, Wang, Liu, Zhang, Ge, Shen, Ma, Zhang, Dong, Zhao, Liu and Pang2021). The role of N availability in the root–shoot relationship was elaborated by Chen et al. (Reference Chen, Liu, Wang, Zhang, Sun, Song, Bai, Lu and Li2020a, Reference Chen, Liu, Zhang, Ge, Shen, Ma, Zhang, Dong, Zhao and Pang2020b), showing that the growth of both the roots and the shoots was controlled by N. Genotype response to N supply is variable due to differences in traits such as root growth and architecture (Xu et al., Reference Xu, Fan and Miller2012; Jiang et al., Reference Jiang, Sun, Tian, Hu, Michel, Gao, Jiang, Cao and Dai2017). Conversely, increasing application rates of N can delay crop maturity but enhance vegetative growth and increase the retention of poor fruit at lower nodes of the main stem (Boquet and Breitenbeck, Reference Boquet and Breitenbeck2000; Hons et al., Reference Hons, McFarland, Lemon, Nichols, Boman, Saladino, Mazac, Jahn and Stapper2004). Because of these variations, measurements of crop growth and maturity should be included in assessments of cotton response to different N application rates. Measuring plant height or the number of main stem nodes over time, as well as using sigmoidal growth curves to determine a higher value for every parameter as well as the rate of growth for any given period in the season, are common methods for assessing plant growth response (Zhao et al., Reference Zhao, Reddy, Kakani, Read and Sullivan2003). When N is applied to cotton before flowering, both the leaf area index (LAI) and blooming improve (Borowski, Reference Borowski2001). Higher N application rates promote vegetative growth, which delays crop maturity (Hons et al., Reference Hons, McFarland, Lemon, Nichols, Boman, Saladino, Mazac, Jahn and Stapper2004).
Experiments were conducted in the conventional cotton belt of southern Punjab, Pakistan, during 2007–2009 and 2018–2020 using a wide range of N application rates. The overall goal of this study was to determine the most efficient N levels for time-series growth, N status, seed cotton yield and time-series biomass for cotton grown in a semi-arid environment.
Materials and methods
Experimental site and cotton cultivars
Cotton has been sown in research fields of the Central Cotton Research Institute (CCRI) since its establishment in 1970. Two field experiments were conducted at the experimental station of CCRI, Multan (30.19°N, 71.47°E and 122 m above sea level), in southern Punjab, Pakistan for three years from 2007 to 2009 and for three years from 2018 to 2020. The soil is sandy clay loam in nature with fractions of sand, silt and clay at 0.39, 0.44 and 0.17, respectively (Table 1). Composite soil samples from variable soil depths were drawn for various physical and chemical properties of the soil at the experimental site (Table 1). The data regarding soil properties showed that the soils were alkaline and low in fertility levels. The chemical properties of the study site were determined using the Bao (Reference Bao2000) methodology. The following techniques, i.e., potassium dichromate-volumetric, alkaline hydrolysis diffusion, Mo–Sb colorimetric, NH4OAc extraction were used for determining soil organic matter, alkali-hydrolysed N, soil available P, and soil available K. The climate of the area is semi-arid with long, hot and dry summers and short, very cold winters. Daily maximum and minimum temperature, rainfall and solar radiation for the cotton growing seasons from 2007 to 2009 (Exp.-I) and 2018 to 2020 (Exp.-II) are presented in Figs 1 and 2, respectively. The cotton cultivar MNH-886 was used for planting during both field experiments (Exp.-I and Exp.-II). This cultivar is commercially grown by farmers in the cotton belt of Pakistan.
Table 1. Physical and chemical properties of study soil


Figure 1. Maximum and minimum temperatures (a, b and c), rainfall and solar radiation (d, e and f) for Exp.-I at the study site during 2007, 2008 and 2009.

Figure 2. Maximum and minimum temperatures (a, b and c), solar radiation and rainfall (d, e and f) for Exp.-II during 2018 (a and d), 2019 (b and e) and 2020 (c and f) growing seasons.
Experimental design, and management practices
A randomized complete block design (RCBD) with three replications was used for both experiments. The treatments in Exp.-I that were conducted from 2007 to 2009 were, i.e., 0, 40, 80, 120, 160, 200 and 240 kg N/ha, while the treatments in Exp.-II conducted from 2018 to 2020 were 0, 70, 140, 210 and 280 kg N/ha. The N levels (as per treatment in both experiments in the form of urea) and application methodology were consistent with other cotton-growing areas in the region and phosphorus and potassium were applied at the rate of 100 kg phosphorus, 90 kg potassium per hectare in the form of di-ammonium phosphate (DAP), and sulphate of potash, respectively. The total 10–12 irrigations were applied as per need of the crop based on soil, crop and environmental conditions. However, the source of water was tube well. Each experimental unit contained six cotton rows with a 75 cm row-to-row distance. Thinning was carried out at the two true-leaf stages to maintain the optimum plant population. Further field management practices, i.e., irrigation, weeding, and insect and disease treatments were uniform and carried out according to recommendations of the extension wing of the Punjab Agricultural Extension Department.
Sampling and measurements
Seed cotton was picked manually from an area of 1 m2 marked from the central rows of each experimental unit. The seed cotton yield of every picking/harvest was sun-dried and weighed to attain seed cotton yield. Eight times plants from an area of 1 m2 from each plot were harvested manually at 15-day intervals during various vegetative and reproductive phases. Biomass was calculated from the sum of plant parts from an area of 1 m2, which were separated into their parts, i.e., leaf, stem and cotton seed. The plant material was packed in craft bags for drying at 105°C for 30 min, subsequently dried at 70°C to obtain constant weight, and then weighed. The separated leaves from an area of 1 m2 were used to calculate leaf area using a leaf area meter, and afterward, LAI was calculated as the ratio of leaf area to ground surface area.
Total plant nitrogen uptake
Leaf, stem and cotton seed N were determined at different growth stages. A pair of scissors was used to separate plant parts (leaf, stem and cotton seed) during both the vegetative and reproductive phases. These components were oven-dried at 75°C, and the oven-dried plant samples were milled and screened through a 0.5 mm sieve. The N concentration was estimated according to the micro-Kjeldahl technique (Bremmer and Mulvaney, Reference Bremmer, Mulvaney, Page, Miller and Keeney1982) and presented on a dry weight basis. The N concentration of leaf, stem and cotton seed was determined individually. For the calculation of N uptake, the N content (%) was multiplied by the biomass of respective plant components (leaf, stem and cotton seed) using the following equation:

Statistical analysis
Cotton crop parameters, i.e., LAI, biomass, seed cotton yield, leaf N, stem N and cotton seed, N data sets were statistically analysed using analysis of variance (ANOVA). To optimize the appropriate N level for cotton under a semi-arid environment, Statistics 8.1 software for RCBD with three replications was used. The main and interaction effects for N treatment means were compared for significance using the least-significant difference (LSD) test at P < 0.01.
Results (Experiments I and II)
LAI
The LAI in cotton increased considerably over time in all treatments in Exp.-I and Exp.-II. The diverse N regimes significantly affected LAI during both sets of field experiments conducted under semi-arid environmental conditions. For Exp.-I (2007–2009) and Exp.-II (2018–2020), LAI increased linearly for all N rates, and the cotton crop fully covered the ground area at 98 days after sowing (DAS) when it reached its peak value, and then decreased progressively (Figs 3–6). For LAI, the highest values were recorded at 240 kg N/ha in Exp.-I and 280 kg N/ha in Exp.-II, while the lowest were recorded at 0 kg N/ha for both experiments (Table 2).

Figure 3. Effect of different nitrogen fertilizer application rates on the LAI and total biomass of cotton for Exp.-I at the study site during 2007 at 0 kg N/ha (a), 40 kg N/ha (b), 80 kg N/ha (c), 120 kg N/ha (d), 160 kg N/ha (e) 200 kg N/ha (f) and 240 kg N/ha (g). Bars represent standard errors.

Figure 4. Effect of different nitrogen fertilizer application rates on the LAI and total biomass of cotton for Exp.-I at the study site during 2008 at 0 kg N/ha (a), 40 kg N/ha (b), 80 kg N/ha (c), 120 kg N/ha (d), 160 kg N/ha (e) 200 kg N/ha (f) and 240 kg N/ha (g). Bars represent standard errors.

Figure 5. Effect of different nitrogen fertilizer application rates on the LAI and total biomass of cotton for Exp.-I at the study site during 2009 at 0 kg N/ha (a), 40 kg N/ha (b), 80 kg N/ha (c), 120 kg N/ha (d), 160 kg N/ha (e) 200 kg N/ha (f) and 240 kg N/ha (g). Bars represent standard errors.

Figure 6. Effect of different nitrogen fertilizer application rates on the LAI of cotton for Exp.-II during the 2018 (a–e), 2019 (f–j) and 2020 (k–o) growing seasons at 0 kg N/ha (a, f and k), 70 kg N/ha (b, g, and l), 140 kg N/ha (c, h, and m), 210 kg N/ha (d, i, and n) and 280 kg N/ha (e, j and o). Bars represent standard errors.
Table 2. Effect of different N fertilizer levels on the maximum LAI, seed cotton yield, lint yield and final biomass for Exp.-I and Exp.-II

LAI, leaf area index; N, nitrogen; LSD, least-significant difference; **P ⩽ 0.01.
Seed cotton yield
The results presented in Table 2 for Exp.-I (2007–2009) and Exp.-II (2018–2020) showed that for both field experiments conducted under semi-arid environments, the seed cotton yield (SCY) was significantly affected by the diverse N regimes. The highest SCY was achieved at 240 kg N/ha in Exp.-I and 280 kg N/ha in Exp.-II and the lowest SCY was achieved at 0 kg N/ha for both experiments.
Lint yield
The results of the two field experiments conducted under semi-arid conditions depicted in Table 2 (Exp.-I; 2007–2009, and Exp.-II; 2018–2020) showed that LY was significantly affected by the diverse N rates. The highest LY was recorded at 240 kg N/ha in Exp.-I and 280 kg N/ha in Exp.-II, while the lowest LY was recorded at 0 kg N/ha for both field studies.
Above-ground biomass
The cotton crop biomass accumulation increased continuously over time for all treatments. The different N regimes significantly affected biomass accumulation for Exp.-I and Exp.-II (Table 2). In Exp.-I (2007–2009) and Exp.-II (2018–2020), biomass accumulation increased progressively for all N rates, and reached its peak values at 160 DAS (Figs 3–5 and 7). The highest biomass accumulation values were recorded at the highest N rates 240 kg N/ha in Exp.-I and at 280 kg N/ha in Exp.-II, while the lowest were recorded at 0 kg N/ha for both experiments (Table 2).

Figure 7. Effect of different nitrogen fertilizer application rates on total biomass (left-side) and seed cotton yield (right-side) for Exp.-II during 2018 (a–e), 2019 (f–j) and 2020 (k–o) cotton growing seasons at 0 kg N/ha (a, f and k), 70 kg N/ha (b, g and l), 140 kg N/ha (c, h and m), 210 kg N/ha (d, i and n), and 280 kg N/ha (e, j and o) at study site. Bars represent standard errors.
Canopy height
The results revealed that for Exp.-II (2018–2020) the canopy height was significantly affected by the N treatments (Table 3). The highest canopy height was recorded at 280 kg N/ha and the lowest canopy height was recorded at 0 kg N/ha. On average across treatments and years, the canopy height ranged from 1.53 to 2.05 m during this study.
Table 3. Effect of different N levels on canopy height, leaf N, stem N and cotton seed N for Exp.-II

N, nitrogen; LSD, least-significant difference; **P ⩽ 0.01.
Leaf nitrogen
The leaf-N was significantly affected by N treatments in Exp.-II (2018–2020). For leaf-N, a linearly increasing trend was observed at all N regimes, values peaked at 81 DAS, and afterward a decreasing trend was noticed (Fig. 8). The results revealed that the year-wise highest leaf-N values (65, 69 and 72 kg/ha; during 2018, 2019, 2020, respectively) were recorded at 210 kg N/ha, while the lowest (9.6, 11.5 and 11 kg/ha; during 2018, 2019, 2020, respectively; Table 3) were recorded at 0 kg N/ha.

Figure 8. Effect of different nitrogen fertilizer application rates on the leaf-N and stem-N for Exp.-II during 2018 (a–e), 2019 (f–j) and 2020 (k–o) cotton growing seasons at 0 kg N/ha (a, f and k), 70 kg N/ha (b, g and l), 140 kg N/ha (c, h and m), 210 kg N/ha (d, i and n), and 280 kg N/ha (e, j and o) at study site. Bars represent standard errors.
Stem nitrogen
During Exp.-II (2018–2020), the stem-N was significantly affected by N regimes (Table 3). The stem-N was measured at a 15-day interval. The stem-N showed a linear increase during study, peaked at 81 DAS, and then a decreasing trend was noticed (Fig. 8). The data showed that year-wise the highest stem-N values (32.1, 32.7 and 32.4 kg/ha; during 2018, 2019, 2020, respectively) were recorded at 280 kg N/ha, while the lowest (4.6, 5.6 and 5.5 kg/ha; during 2018, 2019, 2020, respectively; Table 3) were recorded at 0 kg N/ha.
Cotton seed nitrogen
The results depicted in Table 3 reveal significant differences among N treatments for cotton seed-N. The cotton seed-N presented a sigmoid trend, and its highest values were recorded at 160 DAS (Fig. 9). The highest N rate 280 kg N/ha depicted the highest values (83.1, 101.1 and 100.5 kg/ha; during 2018, 2019, 2020, respectively), while the lowest N rate 0 kg N/ha recorded the lowest values 19, 22 and 24 kg/ha during 2018, 2019, 2020, respectively.

Figure 9. Effect of different nitrogen fertilizer application rates on cotton seed nitrogen for Exp.-II during 2018 (a–e), 2019 (f–j) and 2020 (k–o) cotton growing seasons at 0 kg N/ha (a, f and k), 70 kg N/ha (b, g and l), 140 kg N/ha (c, h and m), 210 kg N/ha (d, i and n), and 280 kg N/ha (e, j and o) at study site. Bars represent standard errors.
Discussion
Cotton growth (LAI and biomass)
In the present study, we evaluated growth time-series LAI and biomass (Exp.-I: 2007–2009; and Exp.-II: 2018–2020) under a semi-arid environment. The study site was selected because it is the highest-yielding area, centrally located in the cotton region, and managed by the Federal Government of Pakistan. The highest values of time course LAI and biomass were reported at 240 kg N/ha in Exp.-I, while in Exp.-II the highest values of canopy height, time-series LAI and biomass were reported at 280 kg N/ha. These optimum N rates in both studies were consistent with findings of large-scale N rate response studies conducted under semi-arid conditions by Ali et al. (Reference Ali, Afzal, Ahmad, Ahmad, Akhtar and Atif2011) and Khan et al. (Reference Khan, Khalid, Fazal, Abdus, Hayatullah, Saleem and Sadia2014). The LAI and biomass accumulation were maximized at the highest applied N rate at the study site in the semi-arid environment in southern Punjab, Pakistan (Wajid et al., Reference Wajid, Ahmad, Khaliq, Alam, Hussain, Hussain, Naseem, Usman and Ahmad2010; Afzal et al., Reference Afzal, Tariq, Ahmad, Mubeen, Khan, Afzal and Ahmad2019). However, Girma et al. (Reference Girma, Teal, Freeman, Boman and Raun2007) and Rochester (Reference Rochester2011) have documented 200 kg N/ha as the optimum rate. The higher amounts of N in the present study caused a considerable increase in time course growth and development during the entire cotton crop duration. The increased LAI and higher biomass were attained during the rapid leaf area development and biomass accumulation period (Liu et al., Reference Liu, Li, Zhang, Cui, Zhan, Xu, Zhang, Dai, Li, Nie, Yang, Li and Dong2022). This increase may also be due to higher N uptake and N accumulation by cotton crop components, i.e., leaf, stem and cotton seed, which might have contributed to the higher LAI and biomass values for the higher N application of 240 or 280 kg/ha. The N being a vital component of cotton crop plant nutrition played a significant role in growth and development. The LAI and biomass were lowest at control or 0 kg/ha N application, and our results were consistent with Tariq et al. (Reference Tariq, Fatima, Iqbal, Nahar, Ahmad and Hasanuzzaman2021), who also found that an improper N level reduces the leaf area development and biomass accumulation in cotton (Tariq et al., Reference Tariq, Yasmeen, Ahmad, Hussain, Afzal and Hasanuzzaman2017, Reference Tariq, Khan, Muhammad, Ahmad and Ahmed2022). Control and lower N rates resulted in lower N availability rates for the crop required for fulfilling the cotton nutritional requirement, thus explaining why lower LAI and biomass were produced at these N rates (Tariq et al., Reference Tariq, Afzal, Muhammad, Ahmad, Shahzad, Kiran and Wakeel2018, Reference Tariq, Khan, Muhammad, Ahmad and Ahmed2022). The lower rates of N were not sufficient to meet the full N demands for cotton growth, which principally restricted leaf area development and biomass accumulation during cotton crop growth.
In-season N status
When the appropriate N rate is applied, efficient N uptake and utilization assure the formation of crop yields. Higher N application rates or availability improved the N supply capacity of the soil, which finally improves crop N uptake and utilization under semi-arid conditions (Ahmad et al., Reference Ahmad, Atique-ur-Rehman, Ejaz, Fatima, Kan and Ahmad2018; Abbas et al., Reference Abbas, Atique-ur-Rehman, Sarwar, Hussain, Abbas, Khan, Fatima, Naz, Younis, Ahmed, Ullah, Iqbal, Iqbal, Kan and Ahmad2020, Reference Abbas, Atique-ur-Rehman, Sarwar, Fatima, Hussain, Ahmed, Raza, Kan, DoĞan, Khan and Ahmad2023; Wahab et al., Reference Wahab, Javed, Khan, Zakir, Abbas, Ahmed, Tariq, Kan, Doğan, Ahmad, Doğan and Fidan2022, Reference Wahab, Hussain, Fatima, Ahmed and Ahmad2024). The study also determined that sampling of crop components in vegetative and reproductive growth phases could be a valuable indicator of in-season crop performance along with final productivity. In both experiments, higher N application rates of 240 kg/ha (Exp.-I) and 280 kg/ha (Exp.-II) in cotton significantly improved N supply capacity at all growth stages and phases compared with control or 0 kg N/ha. Our results showed that higher N uptake by cotton crop components, i.e., leaf, stem and cotton seed, might have contributed to the greatest in-season N status at a higher N application. Rahman et al. (Reference Rahman, Ahmad, Wang, Wajid, Nasim, Hussain, Ahmad, Ahmad, Ali, Ishaque, Awais, Muddasir, Shelia, Ahmad, Fahad, Alam, Ullah and Hoogenboom2018) and Ahmad et al. (Reference Ahmad, Ashfaq, Rasul, Wajid, Ahmad, Khaliq, Nasir, Rasul, Riaz, Ahmad, Ahmad, Baig, Valdivia, Hoogenboom, Rosenzweig, Mutter and Contreras2021) described that higher N fertilizer levels were obligatory to meet critical component (leaf, stem and cotton seed) requirements and higher yields. Their findings corresponded closely with the results presented here where optimal N rates met critical a threshold for cotton growth in a semi-arid environment and hot climate (Ahmad et al., Reference Ahmad, Abbas, Abbas, Fatima, Atique-ur-Rehman, Naz, Younis, Khan, Nasim, Rehman, Ahmad, Rasul, Khan and Hasanuzzaman2017, Reference Ahmad, Ahmad, Ahmad, Ahmad, Wajid, Khaliq, Abbas, Wilkerson and Hoogenboom2023; Nouri et al., Reference Nouri, Lee, Yoder, Jagadamma, Walker, Yin and Arelli2020; Rahman et al., Reference Rahman, Ahmad, Ghaffar, Haider, Ahmad, Ahmad, Tariq, Nasim, Rasul, Fahad, Ahmad, Hoogenboom, Ahmad and Hasanuzzaman2020; Ma et al., Reference Ma, Wang, Li, Wang and Chen2022). Liu et al. (Reference Liu, Li, Zhang, Cui, Zhan, Xu, Zhang, Dai, Li, Nie, Yang, Li and Dong2022) found that occurrences of higher N concentration in cotton parts, i.e., leaf, stem and cotton seed resulted in higher growth and production. The control and lower rates of N application were insufficient to meet the requirements of the cotton crop through efficient uptake and utilization in our experiments, possibly due to lesser N fertilizer uptake and utilization and non-synchronization with cotton N demand during all growth stages and phases (Afzal et al., Reference Afzal, Tariq, Ahmad, Mubeen, Khan, Afzal and Ahmad2019; Tariq et al., Reference Tariq, Abbas, Yasmeen, Ahmad, Ahmad and Hasanuzzaman2020, Reference Tariq, Fatima, Iqbal, Nahar, Ahmad and Hasanuzzaman2021).
Seed cotton yields
In this research, we also determined seed cotton yield in the two independent experiments (Exp.-I: 2007–2009; and Exp.-II: 2018–2020) under semi-arid conditions. The highest values of seed cotton yields were recorded at 240 kg N/ha in Exp.-I, while in Exp.-II at 280 kg N/ha. Seed cotton and final biomass yields gradually increased with increasing N application rates in both experiments (I and II) during all study years. The increased physiological traits and processes at appropriate N rates in cotton plants improved leaf area development and canopy cover that contributed towards higher seed cotton yield and final biomass production (Boquet and Breitenbeck, Reference Boquet and Breitenbeck2000; Yang et al., Reference Yang, Tang, Nie and Zhang2011, Reference Yang, Geng, Huo, Lei, Lang, Li and Liu2021; Ahmad et al., Reference Ahmad, Raza, Ali, Shahzad, Atiq-ur-Rehman and Sarwar2014; Ali et al., Reference Ali, Hameed, Ahmad, Shahzad and Sarwar2014). The higher seed cotton yield was attained by the cotton crop during the rapid reproductive phase (Liu et al., Reference Liu, Li, Zhang, Cui, Zhan, Xu, Zhang, Dai, Li, Nie, Yang, Li and Dong2022). Overall, N nutrition determines the production and quality of cotton crops. Higher N application rates increased photosynthetic rates, in turn leading to a higher accumulation of photosynthates and metabolites and higher production and productivity. These two experiments demonstrated that 240 and 280 kg N/ha are optimum N levels and are adequate quantities to produce higher seed cotton yields under semi-arid environmental conditions. Gormus and El Sabagh (Reference Gormus and El Sabagh2016) found that overall improvement in crop performance and higher yield parameters was recorded with increasing N rates. Geng et al. (Reference Geng, Ma, Zhang, Li, Liu, Lyu and Zheng2015) found that cotton crop requires a continuous supply of N from the soil and that synchronizing N inputs with the needs of cotton is vital for yield formation. Geng et al. (Reference Geng, Ma, Zhang, Li, Liu, Lyu and Zheng2015) and Ghaffar et al. (Reference Ghaffar, Rahman, Ali, Haider, Ahmad, Fahad, Ahmad, Ahmad and Hasanuzzaman2020) found that cotton crop needs not only a high supply of N but, more importantly, a greater proportion of N from the first bloom stage to the initial bloom opening. The results of this study revealed that better vegetative growth is a prerequisite for better reproductive growth in cotton (Figs 3–9) under a semi-arid environment. Synchronization of timely optimum N application rates with optimum crop growth rates and application timing strongly influences time course crop productivity, as does partitioning of growth between vegetation and reproductive organs during vegetative and reproductive phases.
Way forward
Cotton as an indeterminate crop requires special N fertilizer management compared to determinate crops under semi-arid environmental conditions. Therefore, synchronization of timely optimum N application rates with crop N fertilizer requirements for optimum crop growth rates greatly influences the productivity of specific cultivars, especially for partitioning of photo-assimilates between vegetation and reproductive organs as it progresses simultaneously in cotton under a semi-arid environment. It is, therefore, recommended that under semi-arid conditions slightly higher doses of N could be safely used for obtaining higher seed cotton yields without having any negative impact on the environment through leaching or volatilization. The results of this study will be shared with farmers and extension workers of the region through outreach activities at the farmer fields. The N nutrition results for the semi-arid environment are also useful for policymakers and planner to help set priorities for the agricultural sector.
Conclusions
The present study provided updated research results on cotton crop response to different N management regimes for semi-arid environmental conditions. The results showed that there were higher values for LAI, seed cotton yield, lint yield, biomass and cotton plant components (leaf, stem and cotton seed) N-status during the growing season for the higher N application rates. The rates used in these studies were similar to the agriculture department's recommendation of 250 kg/ha for this region.
Acknowledgements
The authors acknowledge the Central Cotton Research Institute (CCRI) support for this research work. The authors also acknowledge the anonymous reviewers for their valuable comments for improving the initial draft.
Author contributions
G. A. and M. T. conceived and designed the study. Z. F. and A. A. W. conducted data gathering. M. A. performed statistical analyses. S. A. wrote the article and C. W. J. and G. H. improved the final draft. All authors wrote the final draft before submission to journal.
Funding statement
The study was financially supported by Bahauddin Zakariya University and the Higher Education Commission of Pakistan.
Competing interests
The authors declare no competing interests.
Ethical standards
Not applicable