The Productivity and Nutrient Use Efficiency of Rice–Rice–Black Gram Cropping Sequence Are Inﬂuenced by Location Speciﬁc Nutrient Management

: Nutrient management is critical for rice farming because the crop is grown under diverse conditions, and in most cases, the existing nutrient management practices fail to achieve an attainable yield target. During recent years site speciﬁc nutrient management gained importance for a target yield with maximum nutrient use efﬁciency. Sufﬁcient research work has not been carried out in this direction so far in the rice–rice–pulse (black gram) sequential cropping system under the red and lateritic belt of West Bengal, India. A multi-locational ﬁeld experiment was conducted from July 2013 to June 2015 at three different locations, namely, Guskara (Burdwan district) and Benuriya (Birbhum district) villages in farmers’ ﬁelds and at the university farm of Visva-Bharati, Sriniketan, West Bengal, India. The performance of nutrients was tested by providing ample doses of N, P, K, S, and Zn compared to the omission of these nutrients. The growth parameters, such as crop biomass production, leaf area index, and number of tillers, and yield attributes and yield were inﬂuenced by nutrient management treatments. Application of 100% of N, P, K, S, and Zn resulted in its superiority to other nutrient management options studied, and a similar trend was also noted with the treatment in the expression of nutrient use efﬁciency (NUE) and nutrient response (NR). The available N, P, K, S, and Zn contents in soil increased steadily due to the increase in fertilizer application. The study concluded that optimization of NPK in the rice–rice–pulse cropping system on target yield along with need-based S and Zn application was beneﬁcial for higher productivity. m − 2 and grains panicle − 1 ), and grain yield of kharif and boro rice were enhanced with balance application ample dose of nutrients, and the residual effects of the treatments were also pronounced in terms of productivity of summer black gram. Further, the agronomic efﬁciency (AE), physiological efﬁciency (PE), and nutrient response (NR) were maximum with the application of an ample dose of nutrients. The study concludes in favor of the requirement of SSNM for optimization of nutrients in the rice–rice–pulse cropping system for a target yield of rice and need-based S and Zn application for higher productivity.


Introduction
Nutrient management is critical for profitable rice farming in Asia because fertilizers incur the second-highest input cost after labor. A study in seven major irrigated rice areas in six Asian countries showed that fertilizers represented 11 to 28% of the annual costs for farmers producing rice [1]. Fertilizers must be applied by adopting the 4R rule, that is, plan related to this rice-rice-pulse (pre-kharif black gram) cropping system. Further, the nutrient use efficiency has recently gained more attention with rising fertilizer costs and continued concern over environmental impairment. Nutrient use efficiency can be viewed from different perspectives based on yield and removal of nutrients. The most common expression of the efficiency of fertilizer nutrient is defined as the percentage of fertilizer recovered in aboveground plant biomass during the growing season. In general, N, P, and K are considered to estimate nutrient use efficiency. However, P and K can be viewed over the long term, while N efficiency is regarded on a short-term basis because of its transient nature. But where there is potential for building soil reserves, long-term N efficiency is appropriate because soil balance also affects N balance. Hence, a multi-location (both the on-station and on-farm) experiment was carried out to optimize nutrients (N, P, K, S, and Zn) with consideration of nutrient balance and nutrient use efficiency for sustaining productivity of the system and generation of awareness among the farmers about improved nutrient management practice in the rice-rice-pulse cropping system.

Site Characteristics
The field experiment was conducted at three different locations, namely, Guskara (Burdwan, 23 • 54 N and 87 • 67 E) and Benuriya (Birbhum district, 23 • 66 N and 87 • 62 E) villages in farmers' fields and at a university farm of Visva-Bharati, Sriniketan, (23 • 66 N and 87 • 65 E) West Bengal, India under the typical Ultisols of the red and lateritic belt. The climate is sub-tropical, greatly influenced by hot, dry summers and cold winters. It falls in the region of the southwest monsoon, and monsoon rains generally start from the end of June and continue up to the middle of October.
The mean average annual rainfall is 1377 mm, out of which nearly 80-90% is received between June and October. The meteorological data related to the weather conditions prevailing during crop seasons (July 2013 to June 2015) with respect to rainfall, relative humidity, and temperature obtained from the agro-meteorological advisory services of the university is presented in Figure 1. The initial soil fertility status, particularly pH, organic carbon content (%), available N (kg ha −1 ), available P 2 O 5 (kg ha −1 ), available K 2 O (kg ha −1 ), available S (kg ha −1 ), and available Zn (mg kg −1 ) was estimated in the laboratory at the beginning of the experiment (Tables 1 and 2). with rising fertilizer costs and continued concern over environmental impairment. Nutrient use efficiency can be viewed from different perspectives based on yield and removal of nutrients. The most common expression of the efficiency of fertilizer nutrient is defined as the percentage of fertilizer recovered in aboveground plant biomass during the growing season. In general, N, P, and K are considered to estimate nutrient use efficiency. However, P and K can be viewed over the long term, while N efficiency is regarded on a shortterm basis because of its transient nature. But where there is potential for building soil reserves, long-term N efficiency is appropriate because soil balance also affects N balance. Hence, a multi-location (both the on-station and on-farm) experiment was carried out to optimize nutrients (N, P, K, S, and Zn) with consideration of nutrient balance and nutrient use efficiency for sustaining productivity of the system and generation of awareness among the farmers about improved nutrient management practice in the rice-rice-pulse cropping system.

Site Characteristics
The field experiment was conducted at three different locations, namely, Guskara (Burdwan, 23°54′ N and 87°67′ E) and Benuriya (Birbhum district, 23°66′ N and 87°62′ E) villages in farmers' fields and at a university farm of Visva-Bharati, Sriniketan, (23°66′ N and 87°65′ E) West Bengal, India under the typical Ultisols of the red and lateritic belt. The climate is sub-tropical, greatly influenced by hot, dry summers and cold winters. It falls in the region of the southwest monsoon, and monsoon rains generally start from the end of June and continue up to the middle of October. The mean average annual rainfall is 1377 mm, out of which nearly 80-90% is received between June and October. The meteorological data related to the weather conditions prevailing during crop seasons (July 2013 to June 2015) with respect to rainfall, relative humidity, and temperature obtained from the agro-meteorological advisory services of the university is presented in Figure 1. The initial soil fertility status, particularly pH, organic carbon content (%), available N (kg ha −1 ), available P2O5 (kg ha −1 ), available K2O (kg ha −1 ), available S (kg ha −1 ), and available Zn (mg kg −1 ) was estimated in the laboratory at the beginning of the experiment (Tables 1  and 2).

Experimental Treatments and Design
The experiment was conducted for two years (six cropping seasons) during 2013-2014 and 2014-2015. The experiment was laid out in a Randomized Block Design with nine treatments (10 m × 10 m plots) and replicated thrice. The following treatments were set: T 1 : ample dose of N, P, K, S, and Zn; T 2 : P, K, S, and Zn, T 3 : N, K, S, and Zn; T 4 : N, P, S, and Zn; T 5 : N, P, K, and Zn; T 6 : N, P, K, and S, T 7 : local variety under the unfertilized check; T 8 : local variety under ample fertilizer; T 9 : control (without any fertilizer). The dose of nutrients in T 1 was 150:60:70 kg ha −1 of N:P 2 O 5 :K 2 O and 200:80:90 kg ha −1 of N:P 2 O 5 :K 2 O in the kharif and boro seasons, respectively. In the case of T 2 : P, K, S, and Zn were applied with an ample dose, and N omitted. In T 3 : N, K, S, and Zn were applied with an ample dose, P omitted, and the same manner was applied for remaining treatments up to T 6 . But, in T 7 (local variety 'Kabirajsal' in was taken in kharif and HYV 'IR 36 in boro) as unfertilized control (without any fertilizer) and in T 9 (where HYV and hybrid were considered as 'MTU 7029 and 'Arize Tej' during the kharif and boro seasons, respectively) with ample dose of fertilizer as applied in T1. However, in T 8 , local variety 'Kabirajsal' in kharif and HYV 'IR 36 in boro were considered with ample dose of fertilizer as per the recommendations (T 1 ) respectively). The ample dose of N, P, and K were calculated based on yield targets that were 5 t ha −1 and 7 t ha −1 in the kharif and boro seasons, respectively.

Experimental Procedure
The fertilizers were applied in the plots after layout as per treatments. The sources of fertilizers for nitrogen, phosphorus, and potash were urea, diammonium phosphate (DAP), and muriate of potash (MOP), respectively. The S was applied at 30 kg ha −1 through bentonite-S (90% S) and Zn at 3 kg ha −1 through Zn-Ethylenediamine tetraacetic acid (EDTA) (12% Zn) in rice as the basal dose during both seasons. After basal application, the top dressing of the remaining nitrogen was applied in two equal splits at maximum tillering and panicle initiation stages during both seasons. The HYV of kharif rice was 'MTU 7029 and hybrid of the boro season was 'Arize Tej' for the treatments T1 to T6 and T9, but for T7 and T8, the local variety of rice 'Kabirajsal' was taken during the kharif season and HYV IR 36 was chosen during the boro season. The black gram cultivar 'Sarada' was grown with residual soil fertility, after the harvest of hybrid boro rice ( Table 3). The crops were irrigated as and when required. The weeds were removed manually at 20 and 40 days after transplanting (DAT). Plant protection measures were taken as per recommended protocols of the university. The second rows from the border of each side of a plot were destructively sampled to record biometric observations, such as dry matter accumulation, number of tillers, and leaf area index (LAI). To determine yield attributes, five plants of each plot were randomly selected and marked. At maturity, these five plants were harvested, dried and their yield attributes were recorded. To determine dry matter accumulation, rice plants were cut at ground level from each plot randomly as destructive samples. Plants of each plot were dried in a hot air oven, kept at 65 • C for 48 h until constant weights were obtained. The dry weights were recorded and used for the determination of dry matter accumulation. The representative green leaves were taken randomly from destructive samples, and their areas were recorded by leaf area meter. The leaves were then dried in a hot air oven at 80 • C for 10 h until constant weights were obtained, and then weights were recorded. The ratio of leaf area weight −1 of these leaves was used to measure the leaf area indices. Since LAI is the area of leaf surface per unit of the land surface [23], it (LAI) was obtained by multiplying this ratio of area and weight with the dry weight of green leaf produced per unit area (square meter) of the land surface. The Leaf Area Index (LAI) during the period of 60 DAT was determined with the following formula (Equation (1)) [24].
Leaf area index = leaf area ground area (1)

Collection and Analysis of Plant and Soil Samples
The crops were harvested manually from the whole plot (net area of 100 m 2 ). Grain and straw sub-samples were drawn from each treatment plot after recording the yields. Further, plant samples were chopped into small pieces and washed in acidified detergent solution, and finally rinsed three times in deionized water. These plant samples were dried in a forced-air circulation oven at 65 • C to bring a constant weight. The samples were then pulverized in a wiring blender, which was cleaned with a hairbrush after each sample, and digested in di-acid (9:4 v/v) of nitric acid (HNO 3 )/perchloric acid (HClO 4 ). The nutrient (N, P, K, S, and Zn) concentration in plant samples was determined by the methods (Table 1) as described by Subbiah and Asija [18]. For N content, grain and straw samples were digested in concentrated H 2 SO 4 -salicylic acid mixture and digestion mixture (potassium sulfate 400 parts, copper sulfate 20 parts, mercuric oxide 3 parts, selenium powder 1 part) and measured by following the micro-Kjeldahl method [18]. For estimation of P and K in plant samples, the vanadate molybdate method and flame photometer method [17] were used. The filtered sample was used for zinc (Zn) estimation with the help of atomic absorption spectrophotometer by adopting DTPA extractable Zn estimation method [21]. Sulfur in whole plant parts and seeds was estimated by the Turbid Metric Method by using a spectrophotometer at 490 λnm, and the concentration was expressed in percentage [22].
Five soil sub-samples (0-15 cm) from each plot were collected before the growing of rice during kharif 2013 (initial soil sample) and after the harvest of the crop in the rabi season of 2015 (post-harvest) using stainless steel tube augers. Initial soil samples were taken from 15cm topsoil. After uniform mixing of sub-samples, representative soil samples were air dried and ground to pass through a 2-mm stainless steel sieve for determination of different soil parameters by standard methods as mentioned in Table 1. The initial soil fertility has been mentioned in Table 2.

Calculation of Nutrient Indices
The agronomic efficiency (AE), physiological efficiency (PE), and nutrient response (NR) were calculated using the equations described by Ray et al. [25] as given below: where Gf is the economic (grain) yield of the fertilized plot (kg) and Go is the economic yield of the nutrient omitted plot (kg), and Na is the quantity of nutrient applied (kg).
where Ntf is the nutrient accumulation by straw/stover and economic product in the fertilized plot (kg) and Nto is the nutrient accumulation by straw/stover and economic product in the nutrient omitted plot (kg).

Calculations and Statistical Analysis
The experimental data were analyzed statistically by using analysis of variance (ANOVA). The standard error of means (SEm±) and critical difference at 5% probability level of significance (CD, p ≤ 0.05) [26]. The Excel software (Microsoft Office Home and Student version 2019-en-us, Microsoft Inc., Redmond, Washington, (USA) was used for statistical analysis and drawing graphs and figures.

Growth Parameters of Crops under Rice-Rice-Black Gram Cropping System
The pooled data of two years study are presented in Table 4 to show the impact of nutrient management on biomass production (gm −2 ) of kharif and boro rice. The observations showed a positive and favorable effect of nutrient management on improving biomass production of rice in different treatments under the study at various locations. The mean data of multi-locational trial indicated that treatments T 1 and T 8 with ample doses of nutrients (N, P, K, S, and Zn) registered significantly higher biomass yield compared to the same varieties/hybrid receiving omission nutrient doses (omission -N, -P, -K, -Zn, -S, and control) during different seasons. The maximum effect of omission was observed in N omission treatment (T 2 ) followed by P (T 3 ), K (T 4 ), and other nutrients.
The data on tiller number (m −2 ) and leaf area index (LAI) of kharif and boro rice revealed that the nutrient management practices influenced the above growth parameters (Tables 5 and 6). The treatments T 1 and T 8 with ample doses of nutrients (N, P, K, S, and Zn) registered significantly higher values of the number of tillers at 90 days after transplanting (DAT) and LAI at 60 DAT than treatments with the omission of N, P, K, Zn, and S and unfertilized control. Application of an ample dose of nutrients assured proper nutrition through uptake and metabolism, which were reflected in the enhancement of growth characters, such as biomass production, number of tillers, and leaf area index. Earlier, Buresh et al. [27] and Shankar et al. [28] noticed that increasing the level of nutrient application influenced the growth parameters of rice significantly. Table 4. Effect of nutrient management on growth attributes of in rice-rice cropping system in multiple locations (pooled data of two years).

Yield Attributes of Crops under Rice-Rice-Black Gram Cropping System
The pooled data on yield attributes of rice grown during both seasons clearly indicated that nutrient management treatments influenced panicle m −2 , number of grains panicle −1 , and test weight (Tables 7-9). The data of multi-locational trials revealed that treatments T 1 and T 8 with an ample dose of nutrients (N, P, K, S, and Zn) registered significantly more panicle m −2 than the same varieties/hybrid grown with omission plots (omission -N, -P, -K, -Zn, -S, and control) during different seasons. A similar type of observation was recorded with the number of grains panicle −1 . The maximum effect of omission was observed in the treatment where N was omitted (T 2 ), which was closely followed by the omission of P (T 3 ) and K (T 4 ) and other nutrients (namely S and Zn). In the case of test weight of rice varieties/hybrid grown in different seasons, much variation was not noted within treatments, and it was nonsignificant. The test weight is a genetic character that may or may not be influenced by nutrient levels. Similar observations on yield attributes of rice were noticed by earlier researchers [7,12,28].

Yield of Crops under Rice-Rice-Black Gram Cropping System
The yield data of the rice-rice-pulse (Pre-kharif black gram) cropping system was recorded for two years from different locations, and pooled data are presented in Table 10 and Figures 2-4. The yield of the rice-based cropping system was positively affected by ample doses of N, P, K, S, and Zn for the target yield of 5 t ha −1 and 7 t ha −1 during the kharif and boro seasons of rice, respectively. Further, the residual effect of nutrients on Pre-kharif pulse as black gram was also satisfactory in terms of productivity where ample doses of nutrients were applied. With the application of an ample dose of N, P, K, S, and Zn in kharif rice, the high yielding variety 'MTU 7029 (T 1 ) yielded much higher than local rice variety 'Kabirajsal' (T 8 ). A similar result was observed between the hybrid rice 'Arize Tej' and HYV IR 36 during the boro season. In both the years, higher available nutrient was received to the rice crop due to the application of 100% N, P, K, S, and Zn (T 1 and T 8 ). Both the treatments expressed higher growth parameters, namely, dry matter production, leaf area index and tillers, and yield attributes (particularly, panicles m −2 and number of grains panicle −1 ), and the impact of these characters was reflected in the productivity of rice during kharif and rabi seasons. Table 10. Effect of nutrient management on grain yield (t/ha) of rice and seed yield of black gram (kg/ha) (pooled of two years). Treatment details are available in Table 4.

Pre-kharif
Sustainability 2021, 13,3222 rice variety 'Kabirajsal' (T8). A similar result was observed between the hybrid rice 'Arize Tej' and HYV IR 36 during the boro season. In both the years, higher available nutrient was received to the rice crop due to the application of 100% N, P, K, S, and Zn (T1 and T8). Both the treatments expressed higher growth parameters, namely, dry matter production, leaf area index and tillers, and yield attributes (particularly, panicles m −2 and number of grains panicle −1 ), and the impact of these characters was reflected in the productivity of rice during kharif and rabi seasons.  Table 4. Standard error mean (SEm±) bar was calculated from three replicates for each treatment.  Table 4. Standard error mean (SEm±) bar was calculated from three replicates for each treatment.  Table 4. Standard error mean (SEm±) bar was calculated from three replicates for each treatment.  Table 4. Standard error mean (SEm±) bar was calculated from three replicates for each treatment.   T1  T2  T3  T4  T5  T6  T7  T8 T1  T2  T3  T4  T5  T6  T7  T8 Table 4. Standard error mean (SEm±) bar was calculated from three replicates for each treatment.  Table 4. Standard error mean (SEm±) bar was calculated from three replicates for each treatment.  Table 4. Standard error mean (SEm±) bar was calculated from three replicates for each treatment.  Table 4. Standard error mean (SEm±) bar was calculated from three replicates for each treatment.   T1  T2  T3  T4  T5  T6  T7  T8 T1  T2  T3  T4  T5  T6  T7  T8 Table 4. Standard error mean (SEm±) bar was calculated from three replicates for each treatment.
On the other hand, the omission of nutrients reduced the availability of essential nutrients in the kharif season and more especially for the succeeding boro rice. The omission of N was more crucial as the mean data showed there was a yield reduction in rice by 43 and 45 percent during the kharif and boro seasons, respectively. Similarly, the omission of P and K reduced kharif rice grain yield by 9.7 and 3.8 percent, respectively; however, in boro, rice grain yield reduction was 11.2 and 8.6 percent, respectively. As the soil is rich in K, the omission of K reduced the grain yield of rice comparatively less than N and P during both seasons. The omission of S and Zn also did not reduce much yield of rice during both seasons. Moreover, pre-kharif black gram raised with residual nutrients also showed a similar trend in terms of productivity as noted in rice.
The omission of N decreased black gram yield by 27%, but the omission of P, K, S, and Zn decreased black gram productivity by 7.4, 2.9, 2.8, and 1.7%, respectively. On the other hand, the omission of nutrients reduced the availability of essential nutrients in the kharif season and more especially for the succeeding boro rice. The omission of N was more crucial as the mean data showed there was a yield reduction in rice by 43 and 45 percent during the kharif and boro seasons, respectively. Similarly, the omission of P and K reduced kharif rice grain yield by 9.7 and 3.8 percent, respectively; however, in boro, rice grain yield reduction was 11.2 and 8.6 percent, respectively. As the soil is rich in K, the omission of K reduced the grain yield of rice comparatively less than N and P during both seasons. The omission of S and Zn also did not reduce much yield of rice during both seasons. Moreover, pre-kharif black gram raised with residual nutrients also showed a similar trend in terms of productivity as noted in rice.
The omission of N decreased black gram yield by 27%, but the omission of P, K, S, and Zn decreased black gram productivity by 7.4, 2.9, 2.8, and 1.7%, respectively.    On the other hand, the omission of nutrients reduced the availability of essential nutrients in the kharif season and more especially for the succeeding boro rice. The omission of N was more crucial as the mean data showed there was a yield reduction in rice by 43 and 45 percent during the kharif and boro seasons, respectively. Similarly, the omission of P and K reduced kharif rice grain yield by 9.7 and 3.8 percent, respectively; however, in boro, rice grain yield reduction was 11.2 and 8.6 percent, respectively. As the soil is rich in K, the omission of K reduced the grain yield of rice comparatively less than N and P during both seasons. The omission of S and Zn also did not reduce much yield of rice during both seasons. Moreover, pre-kharif black gram raised with residual nutrients also showed a similar trend in terms of productivity as noted in rice.
The omission of N decreased black gram yield by 27%, but the omission of P, K, S, and Zn decreased black gram productivity by 7.4, 2.9, 2.8, and 1.7%, respectively.

Nutrient Use Efficiency under Rice-Rice-Black Gram Cropping System
The pooled data of two years data showed the efficiency of applied nutrient in terms of NUE for grain production and presented in Tables 11 and 12. The study highlighted that the balanced nutrient management practices showed a positive effect on agronomic efficiency (AE) and physiological efficiency (PE) in the rice-based sequential cropping system. Table 11. Effect of nutrient management on agronomic efficiency (AE) of rice-rice cropping system (pooled data of two years).   Table 4.

Treatment
The AE and PE during both the seasons were greater in the treatment where crop received an ample dose of N, P, K, S, and Zn through chemical fertilizers than other fertility treatments, namely, the omission -N, -P, -K, -Zn, -S, and control. The results clearly showed that crop nutrition through an ample dose of chemical fertilizers was beneficial

Nutrient Use Efficiency under Rice-Rice-Black Gram Cropping System
The pooled data of two years data showed the efficiency of applied nutrient in terms of NUE for grain production and presented in Tables 11 and 12. The study highlighted that the balanced nutrient management practices showed a positive effect on agronomic efficiency (AE) and physiological efficiency (PE) in the rice-based sequential cropping system. Table 11. Effect of nutrient management on agronomic efficiency (AE) of rice-rice cropping system (pooled data of two years).   Table 4.

Treatment
The AE and PE during both the seasons were greater in the treatment where crop received an ample dose of N, P, K, S, and Zn through chemical fertilizers than other fertility treatments, namely, the omission -N, -P, -K, -Zn, -S, and control. The results clearly showed that crop nutrition through an ample dose of chemical fertilizers was beneficial

Nutrient Use Efficiency under Rice-Rice-Black Gram Cropping System
The pooled data of two years data showed the efficiency of applied nutrient in terms of NUE for grain production and presented in Tables 11 and 12. The study highlighted that the balanced nutrient management practices showed a positive effect on agronomic efficiency (AE) and physiological efficiency (PE) in the rice-based sequential cropping system. Table 11. Effect of nutrient management on agronomic efficiency (AE) of rice-rice cropping system (pooled data of two years).
The AE and PE during both the seasons were greater in the treatment where crop received an ample dose of N, P, K, S, and Zn through chemical fertilizers than other fertility treatments, namely, the omission -N, -P, -K, -Zn, -S, and control. The results clearly showed that crop nutrition through an ample dose of chemical fertilizers was beneficial for improving productivity and nutrient use efficiency. The low response of crops to nitrogenous fertilizers was due to various nitrogen loss mechanisms, namely, ammonia volatilization, leaching, and denitrification. An enhancement of AE and PE was recorded with an ample dose of nutrients because of the proper relationship between source and sink, which ultimately resulted in a positive impact on the productivity of a rice-based cropping system. A similar finding was also reported by Singh and Bansal [32] and Xu et al. [33].

Nutrient Responses of Crops under Rice-Rice-Black Gram Cropping System
The mean data of the multi-location trial clearly revealed that the maximum nutrient response (Table 13) during the kharif season was noted when the crop received an ample dose of N, P, K, S, and Zn through chemical fertilizers, and it was closely followed by the omission of K and S. Similarly, in boro rice, the highest value of nutrient response was recorded with 100% application of N, P, K, S, and Zn, followed by the omission of sulfur and zinc. In the kharif season, the HYV 'MTU 7029 gave a better nutrient response over 'Kabirajsal', while the hybrid rice 'Arize Tej' proved more responsive over 'IR 36 during the boro season.
During the pre-kharif season, the residual effect of nutrients was tested by growing black gram, and the treatment with ample dose of N, P, K, S, and Zn resulted in the maximum nutrient response. The available N, P, K, S, and Zn contents in soil increased steadily due to increased fertilizers application. The highest available N, P, K, S, and Zn contents in soil were recorded from the treatment having an ample dose of N, P, K, S, and Zn, which were markedly higher than all other treatments. The increase in the nutrient of the soil in available form is necessary for sustaining crop productivity. Thus, optimization of NPK in rice-rice-pulse based cropping system on target yield basis along with need-based S and Zn application should be recommended for higher productivity and profitability, as mentioned by [33]. Early research evidence mentioned that interaction among essential nutrients influenced productivity of a rice-based cropping system, and it is important to make fertilizer recommendations for farmers of a locality [34].