A 2-year field trial was conducted at regional research stations (RRS) of Punjab Agricultural University located at Faridkot and Abohar during the summer seasons (April–November) of 2019 and 2020. These Research Stations are located, 96 km apart, in two distinct agro-climatic locations of the southwestern cotton belt of Punjab. The experimental site of RRS Faridkot lies in agro-climatic zone IV at an altitude of 211 m above mean sea level (AMSL) and latitude of 30°40′ N and longitude of 74°44′ E, whereas the experimental site of RRS Abohar lies in agro-climatic zone V at an altitude of 186 m AMSL and latitude of 30°08′ N and longitude of 74°12′ E. Both of the experimental locations fall in the Trans Gangetic plain zone of India with the climate characterized by sub-tropical and semi-arid conditions experiencing a dry and hot summer from mid-April to June and a cold winter season from November to January, thus truly representing the cotton belt of northwestern India. The average annual rainfall of both sites varies from 300 to 400 mm, with 75% of the total precipitation mainly occurring from July to September.

The experiment having eight treatments has been laid out in a randomized complete block design with three replications ( Figure 1 ). It was comprised of six treatment combinations from three levels of SSDI [60% crop evapotranspiration (ETc), 80% ETc, and 100% ETc of PAN evapotranspiration] and two N fertigation levels [75% RDN (84.4 kg N ha-1) and 100% RDN (112.5 kg N ha-1)]. In addition, there were two extra absolute Control treatments, i.e., SF method of irrigation with 100% RDN through broadcasting of urea (Control 1) and SDI at 0.8 ETc coupled with 112.5 kg N ha-1 (Control 2).

Various treatment combinations included the following: T₁—SSDI at 60% ETc with DF of 75% RDN, T₂—SSDI at 60% ETc with DF of 100% RDN; T₃—SSDI at 80% ETc with DF of 75% RDN; T₄—SSDI at 80% ETc with DF of 100% RDN; T₅—SSDI at 100% ETc with DF of 75% RDN; T₆—SSDI at 100% ETc with DF of 100% RDN; and T₇—SF irrigation with 100% RDN through broadcasting of urea–Control 1; and T₈—SDI at 80% ETc with DF of 100% RDN (Control 2).

After harvesting of the wheat crop, the field was ploughed twice, and pre-sowing irrigation (rauni) was applied at both study locations. Bt cotton hybrid cultivar (RCH773 BGII) was sown in the first fortnight of May during both years (2018 and 2019) at Faridkot and Abohar. The sowing was done manually by hand dibbling method using two seeds per hill at a uniform depth of 5 cm. An inter-row spacing of 67.5 cm and intra-row spacing of 75 cm within plants was uniformly maintained. Each experimental plot (60.75 m²) had a total of 125 plants in five rows, while each row length was 18 m.

The sub-surface laterals for SSDI were installed using a tractor-operated machine at a uniform depth of 20 cm, having an inter-row spacing of 67.5 cm with the emitter discharging 2.2 L of water per hour. However, for Control 2, SDI drip laterals with a similar discharge were manually laid out on the ground surface. The Control 1 plots receiving irrigation through SF were delivered with water through dedicated PVC pipes. All the experimental plots received regulated and measured water supply through a water meter and were surrounded by strong bunding all around, with a sufficient buffer area (2.0 m) so that variation owing to water application among different treatments could be minimized. The bulk density of soil has been measured by using the Core method (Black and Hartage, 1986), while soil moisture was determined by employing time domain refractometry (PR2/6 Profile Probe, Delta-T Devices Ltd. UK) after calibrating it with gravimetric methods (Richards and Weaver, 1943).

Growth- and yield-attributing parameters like plant height, sympodial branches per plant, bolls per plant, and boll weight have been recorded from 10 randomly selected plants per plot. Biomass accumulation was measured by cutting the whole aboveground portion of the plant after harvest and then weighing after thorough drying under the sun. SCY constituted the total weight of both the hand pickings and is expressed in kilograms per hectare (Mishra et al., 2021).

Irrigation under the SDI and SSDI method was applied at an interval of every 5 days. The total amount of applied irrigation water varied according to treatments, which were based on the amount of crop evapo-transpiration. In the SF method (conventional), water was applied by flood irrigation. A total of six and four SF irrigations were applied at the Faridkot and Abohar locations, respectively, during 2019, and the corresponding values were four and five during 2020. In SF plots, the first irrigation was given at 35 days after sowing (DAS) in mid-June at all locations and study years. The second, third, fourth, and fifth (last) SF irrigation, respectively, was applied during the end June, July, August, and September, respectively. The quantity of each irrigation was measured by using a water meter in all the treatment plots, and water applied under all irrigations was summed up to calculate the total irrigation water.

To work out the ETc, reference evapotranspiration (ETo) was calculated with a calculator developed by FAO from site-specific weather data for both locations. According to the user manual given by FAO, to calculate the daily crop evapotranspiration (ETc), the value of the crop coefficient (K_c) of cotton crop was 0.75 until the end of June, 1.15 during July and August, and 0.70 until boll picking. Fertilizer N in the form of urea (46% N) was applied for SDF and SSDF treatments starting from 35 DAS and delivered in 10 equal splits after every 5 days except for the SF, where 100% RDN was applied in two equal splits, i.e., half dose at thinning and the remaining half at flowering. The application of phosphorous was skipped as the recommended dose of P was applied to the preceding wheat as recommended by Punjab Agricultural University (Anonymous, 2022). Actual crop evapotranspiration (ET_a) has been worked out with the help of the soil water balance equation (Fernández et al., 2020):

where IW is amount of irrigation water applied (mm), P represents precipitation in mm, R indicates surface runoff (mm), D is deep drainage (mm), and ΔS represents soil profile moisture change (mm). Runoff was nil as ridges/buffers surrounded all the plots. Deep drainage was assumed to be zero when moisture storage in the soil profile was lesser than the field capacity and when soil moisture storage (SMS) surpassed the field capacity storage either after a rain or due to irrigation. Thereafter, deep drainage has been developed as the gap between field capacity storage and SMS plus rain/irrigation. Since the water table at both study sites was below 3 m, an upward flux from the groundwater was not considered. Two-meter wide buffers were established between various plots to eliminate water fluxes in the vicinity of root zone laterals.

Bio-physical water productivity (BWP) and economic water productivity (EWP) were computed by using the following equations (Perry et al., 2017; Fernández et al., 2020):

where SCY means seed cotton yield (kg ha^-1), ET_a is actual crop evapotranspiration (m³ ha^-1), and NR indicates net returns ($ ha^-1).

WUE was calculated by using the following equations (Perry et al., 2009; Sahoo et al., 2018):

where IW is irrigation water applied (m³ ha^-1), R represents rainfall (m³ ha^-1), and ±ΔS was change in soil profile moisture (m³ ha^-1), while ETa indicates actual crop evapotranspiration (m³ ha^-1).

Nitrogen use efficiency (NUE) is calculated using the formula given below:

Lint samples were obtained by ginning the clean and dried weighed samples of seed cotton through a single-roller electric gin, and ginning turnout (GTO) was computed by using the following formula:

Fiber samples were sent to the laboratory of ICAR-Central Institute for Research on Cotton Technology (ICAR-CIRCOT), Mumbai, for measurement of various fiber parameters like halo length, fiber strength, uniformity index, micronaire, etc. A sample of 100-gram lint was taken to measure the micronaire value by using Prectitronic Digital Mic Tester at C, Mumbai.

The total cost of cultivation incurred for raising cotton crop was calculated with the help of Enterprise budget (2021) of Kharif crops by the Department of Economics and Sociology, PAU, Ludhiana. Gross returns of different treatments were worked out by multiplying the SCY from respective treatments with the prevalent market price of $0.75 kg^-1 of seed cotton. The benefit/cost ratio (B:C) was worked out to check the economic feasibility of treatments and is calculated by dividing the net returns by the total cost of cultivation (Singh et al., 2020).

Statistical analysis of various recorded parameters and calculated indices was performed to evaluate the effect of various treatments (SSDI regimes and N fertigation levels), in comparison with control treatments, using SAS Proc GLM (SAS software 9.3, SAS Institute Ltd., USA). Significant mean differences were compared using Fisher’s least significant difference test at a probability level of 5%. The variance of data in both sites was homogeneous according to Bartlett’s test (p ≤ 0.05), so the data of both sites and years was pooled and analyzed.

The data pertaining to weather-related parameters has been recorded from the agro-meteorological observatories of RRS, Faridkot, and RRS, Abohar and presented in Figures 2 , 3 .

Among the different treatments, the actual crop evapotranspiration (ET_a) ranged from 345.8 to 588.8 mm at both locations over the years. During 2019, except for the SF, negative soil moisture confirmed the necessity of underground water withdrawal for meeting irrigation demand that was substantially improved in the later year for each treatment at both experimental sites. The effects of the available moisture in the soil profile witnessed the heavy deep drainage during 2020 compared with the year 2019. Similarly, a huge surface runoff, i.e., 86.4 and 88.2 mm at Faridkot and 92.3 and 93.3 mm at Abohar, was recorded for the surface flood method which was tremendously reduced under drip treatments. The data on the physio-chemical characteristics of soil at the experimental sites is given in Table 1 . The layer-wise moisture retention capacity (at field capacity and permanent wilting point, PWP) of soil at both experimental sites for soil profile (0–150 cm) has been worked out and presented to work out various water indices ( Table 2 ). Layer wise bulk density of both experimental sites is presented in Table 3 . The detailed information on cultivar, planting geometry, sowing and harvesting dates, etc., is summarized ( Table 4 ).

Different irrigation and fertigation treatments exerted a significant effect on growth parameters like plant height and biomass accumulation ( Table 5 ). Among the tested irrigation regimes, taller plants were recorded under SDI (170.0 cm), closely followed by a SSDI level of 80% ETc (164.3 cm) and 100% ETc (167.4 cm), while a significant reduction in plant height was observed under 60% ETc (148.1 cm) and SF (156.1 cm). Biomass accumulation has been highest under a SSDI level of 100% ETc (2,627 g m^-2) compared with 80% ETc (2439 g m^-2) and 60% ETc (1978 g m^-2). The reduced plant height and biomass accumulation under SSDI level of 60% ETc might be due to the fact that here least water was supplied, which, in turn, failed to maintain optimal crop growth. The better plant vigor under SDI over SF in the present findings is in close proximity with those of Yadav and Chauhan (2016), who observed taller plants under SDI over SF. Furthermore, SF not only recorded shorter plants in line with Prajapati and Subbaiah (2018) but also lesser biomass accumulation (2,138 g m^-2) in comparison with various SDF/SSDF treatments under 80%/100% ETc, elucidating that applied water and nitrogen could not be efficiently utilized ( Table 5 ).

Among N fertigation schedules, higher plant height and biomass accumulation was evident under 100% RDN (2,420 g m^-2) over the 75% RDN (2276 g m^-2). Nevertheless, both SSDF levels of 75% RDN and 100% RDN revealed better plant height and biomass accumulation over the SF method (2,138 g m^-2), which could be attributed to the optimum availability of water and nitrogen to plants under SSDF, in agreement with Ayyadurai and Manickasundaram (2014) who recorded 36% higher biomass under DF of 100% RDN over soil application under SF. These findings established water and nitrogen to be among the essential growth factors as evident from better height and biomass accumulation under their increased supply (Brar et al., 2021). Higher biomass accumulation under SSDI over SF method is also supported by Sampthkumar et al. (2006) who recorded improved biomass accumulation under DI over border strip and SF irrigation.

Sympodial branches, bolls per plant, boll weight, and SCY varied significantly under various irrigation regimes ( Table 5 ). Higher sympodial branches per plant were recorded under SDI (24.2), closely followed by SSDI of 100% ETc (23.1), 80% ETc (22.2), and SF method (20.5), while the number was least under SSDI of 60% ETc (17.6). The observation on improved sympodial branches per plant under SSDI is in conformity with that of Ali et al. (2017) who reported 15% higher sympodial branches under SDI over the furrow method. Higher bolls per plant was revealed under a SSDI of 100% ETc (60.2) and 80% ETc (57.2) over SSDI of 60% ETc (46.7) and SF method (46.5), which was in agreement with Prajapati and Subbaiah (2018) who recorded improved bolls per plant under SDI over the furrow method. A significant reduction by 23% and 29% for bolls per plant under the SF method over a SSDI of 80% ETc and 100% ETc was evident. Higher boll weight under SDI (4.11 g) was closely followed by SSDI at 100% ETc (4.04 g) and 80% ETc (3.95 g), while SSDI at 60% ETc recorded the significantly lowest value (3.64 g) primarily due to reduced water supply (Singh et al., 2018).

Among SSDI levels, SCY was maximum at 100% ETc (3,305 kg ha^-1), closely followed by 80% ETc (3,240 kg ha^-1) and with the least value under 60% ETc (2619 kg ha^-1), where it was 26% and 23.7% lower compared with 80% ETc and 100% ETc, respectively. Nevertheless, SDI recorded statistically at par SCY (3,300 kg ha^-1) with SSDI of 80% ETc and 100% ETc. However, SSDI of 80% and 100% ETc resulted in 18.7% and 21% higher SCY over SF (2,728 kg ha^-1) due to better yield parameters (Aladakatti et al., 2012). This was primarily due to the improved boll count per plant which was 23% and 29.4% higher under SSDI level of 80% ETc and 100% ETc, respectively, over the SF method. Singh et al. (2018) also reported higher SCY by 19% and 23% over the conventional SF method under SDI of 100% ETc and 80% ETc, owing to the improved boll count. Better SCY and yield attributes such as higher bolls per plant under SSDI of 100% ETc and 80% ETc over SF irrigation are also supported by Neelakanth et al. (2019).

A fertigation level of 100% RDN elucidated better sympodial branches per plant (21.9), boll weight (3.92 g), and bolls per plant (57.2) over 75% RDN (20.0, 3.83 g, and 52.2, respectively) in conformity with Singh and Bhati (2018) who observed higher bolls per plant under SDI of 100% RDN (50.8) over the 75% RDN (46.2). Among N levels, 100% RDN revealed significantly higher SCY by 11.9% and 18.2% over the 75% RDN (2,882 kg ha^-1) and broadcasting method (2,728 kg ha^-1), respectively. Furthermore, the data elucidated that SSDF of 100% RDN either at 80% (3,455 kg ha^-1) or 100% ETc (3,477 kg ha^-1) improved SCY by 4.6% and 5.3%, respectively, over SDF.

Seed cotton yield revealed a positive and linear relationship ( Figure 4 ) with crop evapotranspiration (ET_a), while its relationship with irrigation water applied (IWA) followed a second-order polynomial trend ( Figure 5 ). Furthermore, the R ² values for IWA and SCY for 2019 ( Figure 5A ), 2020 ( Figure 5B ), averaged over the years ( Figure 5C ), and averaged over the locations and years were significant ( Figure 5D ). The regression equations for individual year and location and also when averaged over the years and locations clearly elucidate SCY to be dependent upon ET_a ( Figures 4A–D ) and IWA ( Figures 5A–D ). A correlation heat map among various traits of cotton is given in Figure 6 . This further signifies that efficient usage of water and nutrients can be made under SSDF as evident from better yield realization. These results clearly established that SSDF could play a pivotal role in improving the yield attributes and SCY compared with SF and the soil application of nutrients. Nevertheless, improved boll weight and bolls per plant under SSDF have remained as the primary reasons for the higher SCY. The boll count per plant was 12.2% and 23.0% higher under the SSDF of 75% RDN and 100% RDN, respectively, over broadcasting of 100% RDN in SF. These findings reveal that both water and nitrogen greatly govern the yield as envisaged from the fact that, by increasing the water and N fertigation level, crop yield tends to improve (Dar et al., 2017; Thanappan et al., 2020).

A significant reduction in yield attributes and SCY is evident under the least level of water and nitrogen (i.e., SSDI of 60% ETc along with fertigation of 75% RDN), where the supply of both water and nitrogen was minimal ( Table 5 ). SCY improved significantly by 23.7%, while the increasing irrigation level from 60% ETc to 80% ETc was in conformity with Singh et al. (2018). However, there was only a marginal increase of 2% while moving from 80% ETc to 100% ETc. The present data manifested that, under 100% RDN, SSDI—at either 80% ETc or 100% ETc—remained on par for SCY, though it increased by 708 and 730 kg ha^-1, respectively, over the SSDF at 60% ETc. This revealed that SCY exhibited a tendency to increase between the irrigation levels of 80% ETc and 100% ETc but was optimized at 80% ETc (Kumar et al., 2021).

The total irrigation water applied has been highest under SF (Control 1), i.e., 550 and 400 mm during 2019, with a value of 400 and 475 mm during 2020 for Faridkot and Abohar, respectively ( Tables 6 , 7 ). However, actual crop evapotranspiration (ET_a) remained higher for SSDI at 100% ETc with 100% RDN (580.6 and 462.8 mm for 2019; 582.6 and 588.8 mm for 2020) for Faridkot and Abohar, respectively. A general reduction for ET_a was observed under SF (Control 1) compared with SSDI treatments across the locations. The ETa followed a linear relationship with SCY among the studied treatments over the locations during both study years ( Figures 4A–D ). During 2019, ETa was highest under T₆ (i.e., 580.6 and 462.8 mm for Faridkot and Abohar, respectively), while during 2020, T₄ (584.8 _mm) and T₈ (590.2 mm) revealed higher ETa values for Faridkot and Abohar, respectively. These findings clearly elucidated a saving of 36–54.3% and 52.4-61.4% of irrigation water under SSDI for Faridkot and Abohar, respectively, during 2019 with increased ET_a over SF (Cetin and Kara, 2019). During 2020, the corresponding values ranged from 19.8–41.8% to 30.0–50.5%. The presented results on savings of water under SSDI are well supported by Patil et al. (2009) who observed 45.6% water saving under DI at 100% ETc over the SF method. Roopashree et al. (2016) also reported a saving of 21.7% irrigation water under SSDI over SF.

Among the studied irrigation regimes, improved BWP values ( Table 8 ) were recorded under SSDI levels of 100% ETc (0.627 kg m^-3), 80% ETc (0.620 kg m^-3), and SDI (0.637 kg m^-3), while it was significantly reduced under 60% ETc (0.551 kg m^-3) and SF (0.601 kg m^-3). The reduced BWP by 12.5% under SSDI at 60% ETc compared with 80% ETc was clearly indicative of its poor efficiency (Brar et al., 2021). The SF method revealed least WUE (57%), while values were considerably improved under a SSDI level of 0.8 ETc (80.3%), closely followed by SDI (80.9%), SSDI level of 60% ETc (76.5%), and 100% ETc (78.8%). The better WUE under SDI/SSDI at 0.8 ETc was in conformity with Shruti and Aladakatti (2017), who elucidated higher values under DI applied at 80% ETc.

Among the fertigation levels, 100% RDN recorded a numerically higher BWP of 0.614 kg m^-3 over 75% RDN and SF. However, SSDF of 75% and 100% RDN exhibited better WUE over broadcasting fertilizer (Jayakumar et al., 2015). These results clearly established that SSDF of either 80 or 100% ETc along with 100% RDN recorded better BWP and WUE over the SF method despite the limited application of water. Thus, SSDF can become an effective tool for increasing the cotton productivity as well as water use efficiency without sacrificing SCY (Singh et al., 2022). Further details are elaborated in the following subsections.

Among the studied treatments, the mean evapotranspiration (Eta) at Faridkot during 2019 was 504.6, 530.2, 553.9, 576.8, 554.6, 580.6, 485.6, and 555.2 mm for T₁ to T₈, respectively. During 2020, the respective ETa values were 515.3, 534.7, 572.3, 584.8, 573.4, 582.6, 486.7, and 562.6 mm. At Abohar, the ETa values during 2019 were 380.1, 385.6, 404.6, 451.1, 403.7, 462.8, 345.8, and 410.8 mm for T₁ to T₈, respectively. The respective ETa value during 2020 was 517.4, 541.2, 530.5, 585.3, 543.4, 588.8, 544.6, and 590.2 mm. Reduction in ETa at Abohar during 2019 compared with a later year was observed, while at Faridkot the values were akin during both study years ( Tables 6 , 7 ). Higher ETa values for treatments receiving more water and fertilizer under drip fertigation resulted in vigorous cotton growth, which might have led to variability in the uptake of water and its distribution. The higher ETa in 2020 might be owing to more rainfall which was primarily accountable for continuous soil wetting, thus leading to a higher evaporative loss (Dar et al., 2017). The variation in ETa over the seasons and locations was owing to the huge variation in the amount of rainfall and its distribution ( Figure 2 ), which could have resulted into distinct effects on soil wetting followed by evaporative loss and crop water uptake. The variation among irrigation regimes is primarily owing to the variable quantity of IWA ( Tables 6 , 7 ). The higher ETa was observed under T₆, which received the maximum quantity of irrigation through SSDF, as a linear relationship has been observed between SCY and ETa ( Figure 4 ) with the R ² value ranging from 0.723 to 0.983 for individual and pooled data over the seasons and locations, respectively. The increase in ETa with IWA is in agreement with Irmak et al. (2016).

The higher BWP was recorded in T₄ treatment (0.638 kg m^-3), while the lowest (0.533 kg m^-3) was observed in T₁ treatment ( Table 8 ). The BWP declined owing to a huge reduction in SCY ( Table 5 ) compared with the reduced crop ETa ( Tables 6 , 7 ). Among drip irrigation regimes, the BWP values were 0.551, 0.620, and 0.627 kg m^-3 for 60% ETc, 80% ETc, and 100% ETc, respectively ( Table 8 ). The improvement in BWP with any hike in irrigation level could be owing to the reason that IWA just fulfilled the required soil water deficit. This means that treatments which received more irrigation water experienced lesser water stress compared with treatments receiving lesser irrigation. As a result, the SCY was statistically better under 80% ETc and 100% ETc treatments having a high IWA compared with 60% ETc and resulted in higher BWP (Sun et al., 2006; Dar et al., 2017). When averaged over the years, a linear relationship between ETa and SCY was observed, with R ² value of 0.723 at Abohar and 0.983 for Faridkot and R ² of 0.936 ( Figure 4C ) and R ² of 0.891 when averaged over the locations and years. This clearly elucidated that the efficient use of water by the cotton plants increased the SCY. Therefore, irrigating cotton with an adequate amount of water at regular intervals through drip might enhance the SCY and water productivity with additional saving of water compared with surface flood method.

The WUE in SSDF treatments ranged from 75.4 to 83.5%, compared with SF (57.0%). The highest WUE of 83.5% was recorded in T₄ treatment, i.e., drip irrigation at 80% ETc along with 100% RDN ( Table 8 ), which elucidated that a huge water can be saved by adopting novel water-savvy techniques. During 2019, applied water under SF has been 56% and 110% higher compared with what was applied in SSDI at 100% ETc for Faridkot and Abohar, respectively. The corresponding values were higher by 24.6% and 42.5% during 2020, which was indicative of huge water saving.

Nitrogen is one of the most important limiting factors governing the growth and productivity of crop plants. The present findings clearly established a reduction in the NUE with each increase in N level. Within irrigation regimes, SSDI of 100% ETc recorded a higher NUE (33.9 kg SCY kg^-1 N) which was closely followed by 80% ETc (33.2 kg SCY kg^-1 N), while 60% ETc resulted in a significantly lowest value (26.9 kg SCY kg^-1 N). Furthermore, NUE under SF (Control 1) has been 3.8%, 28.1%, and 30.8% lesser compared with SSDI of 60% ETc, 80% ETc, and 100% ETc, respectively ( Table 8 ).

Among fertigation levels, SSDF of 75% RDN (34.1 kg SCY kg^-1 N) recorded a statistically better NUE by 19.2% over 100% RDN (28.6 kg SCY kg^-1 N), which was in agreement with Singh et al. (2018) who reported 27% higher NUE under 75% RDN over 100% RDN. However, both SSDF levels of 75% and 100% RDN exhibited increased NUE by 31.6% and 10.4%, respectively, over SF (25.9 kg SCY kg^-1 N). The higher NUE under SSDF has been primarily due to the application of nitrogen directly to the root zone, which, in turn, minimized volatilization and leaching losses and consequently improved the yield (Brar et al., 2022). Moreover, the application of N in 10 equal splits under SSDF might have checked the N losses in the soil compared with only two splits under the broadcasting method, which further helped to enhance NUE, in agreement with Bharathraj et al. (2015) wherein 47% higher efficiency has been observed under fertigation compared with soil application.

Ginning turnout was significantly reduced under SF (31.9%) compared with all SDI/SSDF treatments ( Table 9 ). However, fiber strength was significantly higher (30.5) under a SSDF combination of 80% ETc along with 100% of RDN, which was in conformity with Zhang et al. (2019). However, fertigation treatments could not differentiate much for either of the studied quality parameters (Magare et al., 2018; Zahid et al., 2021).

A higher gross return ($2,572.9 ha^-1), net return ($1,648.9 ha^-1), and B:C ratio (1.78) have been recorded in the treatment receiving SSDF at 100% ETc with fertigation of 100% RDN. However, SSDI of 80% ETc with 100% RDN recorded much lesser net returns by $12.0 ha^-1, but with an additional saving of 20% irrigation water ( Table 10 ). A statistically lowest gross return ($1,842.6 ha^-1), net return ($932.0 ha^-1), and B:C (1.02) was observed under SSDI at 0.6 ETc with 75% RDN, thus making it the least remunerative among the drip combinations studied. Nevertheless, SSDI at 0.8 ETc with 100% RDN recorded significantly higher net returns by 54.1% over the SF method. These findings get fair support from Aladakatti et al. (2012) and Pawar et al. (2015). Reduced B:C under the conventional practice of the SF method (1.11) further substantiated that fertilizer and water application through SSDF are more rewarding (Neelakanth et al., 2019).

Economic water productivity remained higher under the SSDI level of 100% ETc and SDI, closely followed by SSDI of 80% ETc, while the least values were exhibited under 60% ETc ($0.218 m^-3) and the SF method ($0.233 m^-3). EWP under SF was lower by 17.9% and 20%, respectively, over the SSDI levels of 80% ETc and 100% ETc, which indicated it to be inferior compared with other treatments. Among fertigation levels, SSDF of 100% RDN recorded better EWP ($0.279 m^-3) over 75% RDN. However, both SSDF levels of 75% and 100% RDN exhibited increased EWP by 7% and 17%, respectively, over Control 1. The SSDF of 80% ETc with 100% RDN and 100% ETc with 100% RDN recorded 29.6% higher EWP over SF. These findings clearly elucidated the advantage of SSDF in increasing monetary advantage over SF.