2.1 Experimental site description
The experiment was established on farmers’ fields in the Bogura district (Sherpur sub-district at 25°01'0.12" N and 89°19'0.12" E), Cumilla district (Chandina sub-district at 23°45'0.57" N and 90° 99'0.77" E), and Mymensingh district (Muktagacha sub-district at 24°45'53.42" N and 90°15'25.13" E) of Bangladesh (Fig. 1). The Mymensingh site is located in Agro-ecological Zone 28, which is comprised of Madhupur Tract loam soil (Eutric Gleysols) and contains 49% sand, 40% silt, and 11% clay, and has a bulk density of 1.2 g/cm3, pH of 5.3, total organic C of 8.12 g/kg, total N of 0.73 g/kg, available P of 5.9 mg/kg, exchangeable K of 43 mg/kg, and available S of 12.8 mg/kg. In contrast, the Bogura site with a Grey Karatoya Bangali Floodplain soil (Aeric Haplaquept) is located in Agro-ecological Zone 4, which contains sand 32%, silt 53%, and clay 15%, with a bulk density of 1.26 g/cm3, pH of 5.5, total organic C of 8.1 g/kg, total N of 1 g/kg, available P of 55.5 mg/kg, exchangeable K of 46.9 mg/kg, and available S of 16.4 mg/kg. The Cumilla site comprised of Meghna Estuarine Floodplain silt loam soil (Aeric Haplaquept) situated in Agro-ecological Zone 19 and contains sand 15%, silt 73%, and clay 13%, with a bulk density of 1.32 g/cm3, pH of 6.3, total organic C of 10.3 g/kg, total N of 1.1 g/kg, available P of 56.9 mg/kg, exchangeable K of 58.7 mg/kg, and available S of 16.3 mg/kg. All the sites are under sub-tropical climatic conditions and the maximum and minimum rainfall data are presented in Fig. 2.
2.2 Experimental design and crop management
A 2 x 4 split-plot experiment in randomized complete design with three replications was conducted during the 2020–2021 and 2021–2022 crop growing seasons at all sites. Crop residue retention was assigned as no residue (NR) and 30% residue (CR) in the main plots, and N fertilization rates were assigned as control (CL), farmers' practice (FP), recommended N dose (RD), and 125% of RD (125RD) in the sub-plots.
The cropping sequences were Potato-Boro-T., Aus-T. and Aman rice (P-R-R-R) in Bogura; Boro-T., Aus-T. and Aman rice (R-R-R) in Cumilla; and Maize-T. and Aman rice (M-R) in Mymensingh. For potatoes, CL (0), FP @ 138, RD @ 135, and 125RD @ 169 kg N/ha were applied. In Boro rice, CL (0), FP @ 120, RD @ 150, and 125RD @ 187 kg N/ha were applied. In T. Aus rice, CL (0), FP @ 77, RD @ 74, and 125RD @ 92 kg/N ha were applied, whereas T. Aman rice had CL (0), FP @ 98, RD @ 90, and 125RD @ 122 kg/N. Maize received CL (0), FP @ 152, RD @ 225, and 125RD @ 281 kg/N. While FP rates were assessed (n = 60) from a survey, the RD rates were based on Bangladesh national fertilizer recommendations (FRG, 2018).
Potatoes were grown at the Bogura site from October to February, while maize was grown at the Mymensingh site from December to May (midwinter to pre-monsoon season). Both T. Aus and T. Aman rice crops were grown from May to July (monsoon) and August to November (post monsoon to prewinter), respectively, in 2020-21 (First year) and 2021-22 (Second year). Moreover, transplanted Boro rice was grown from January to May (post-winter to pre-monsoon) at both the Bogura and Cumilla sites.
The N, P, K, S, Mg, Zn, and B fertilizer sources were urea, triple super phosphate, muriate of potash, gypsum, magnesium sulphate, zinc sulphate, and boric acid, respectively. During the final stage of land preparation, P, K, S, Mg, Zn, and B were incorporated into the soil as basal application following the Bangladesh national fertilizer recommendation guide. For maize, the N fertilizer was applied in three equal splits at 15, 35, and 55 days following seed sowing, while for potatoes, it was applied in two equal splits at 15 and 35 days after planting. For Boro rice, N was applied at 12, 37, and 57 days after transplantation. In contrast, N fertilizer was applied at 12 and 38 days in T. Aus and 12, 27, and 35 days in T. Aman rice after transplantation.
Potato (Diamond, cv.) and maize (BARI hybrid maize-14) seeds were planted/sown at the rate of 2 and 30 kg/ha, respectively. The seed-to-seed distance in each line of maize seeds was 30 cm, and the lines were spaced 60 cm apart. Potato seeds were sown on ridges at 60 cm spacing and 20 cm apart. Boro, T. Aus, and T. Aman rice were transplanted at the rate of 40 kg/ha with 20 cm spacing. The rice cultivars were BRRI-49, BINA-19, BRRI-71 in Boro, T. Aus, and T. Aman, respectively, at all the sites.
Three days prior to final land preparation, the field was sprayed with non-selective herbicide glyphosate (Round up®) at the rate of 1.85 kg/ha. Furthermore, Pretilachlor (Superhit®), a post-emergence herbicide, was applied at the rate of 450 g/ha seven days after rice seedlings transplantation. When necessary, Brifer 5G and Cidial 5G (ACI Bangladesh Ltd.) were implemented.
2.3 Crop yield and N-use efficiency
For all rice crops, a 4 m2 area in the middle of each replicated plot was harvested. For maize, a 2 m2 area was harvested. Grain yield (GY) of Boro, T. Aus, T. Aman, maize, and tuber yield of potatoes was calculated in a dry-matter basis by using Eq. 1 (Akter et al., 2023).
GY or tuber yield (ton/ha) = \(\frac{\text{G}\text{Y}, \text{t}\text{u}\text{b}\text{e}\text{r} \text{y}\text{i}\text{e}\text{l}\text{d} \left(\text{k}\text{g}\right) \times 10}{\text{A}\text{r}\text{e}\text{a} \text{o}\text{f} \text{t}\text{h}\text{e} \text{p}\text{l}\text{o}\text{t} \left({\text{m}}^{2}\right)}\) Eq. 1
The rice and maize grain and straw, potato tuber, and haulm N uptake were determined at oven-dried moisture basis.
Total plant N was determined by the Kjeldahl method as described by Jackson (1973) and expressed as kg/ha. The Eq. 2 was used to calculate the NUE (Moll et al., 1982):
NUE =\(\frac{{\text{N} \text{u}\text{p}\text{t}\text{a}\text{k}\text{e} }_{\text{F}\text{e}\text{r}\text{t}\text{i}\text{l}\text{i}\text{z}\text{e}\text{d} \text{p}\text{l}\text{o}\text{t}} \left(\text{k}\text{g}/\text{h}\text{a}\right)-{\text{N} \text{u}\text{p}\text{t}\text{a}\text{k}\text{e} }_{\text{C}\text{o}\text{n}\text{t}\text{r}\text{o}\text{l} \text{p}\text{l}\text{o}\text{t}} \left(\text{k}\text{g}/\text{h}\text{a}\right) }{\text{T}\text{o}\text{t}\text{a}\text{l} \text{N} \text{i}\text{n}\text{p}\text{u}\text{t} (\text{k}\text{g}/\text{h}\text{a})}\) Eq. 2
Here, N uptake (kg/ha) was calculated as per the Eq. 3 (Assefa et al., 2021).
N uptake (kg/ha) = \(\frac{\text{N}\text{u}\text{t}\text{r}\text{i}\text{e}\text{n}\text{t} \text{c}\text{o}\text{n}\text{c}\text{e}\text{n}\text{t}\text{r}\text{a}\text{t}\text{i}\text{o}\text{n} \left(\text{%}\right)\times \text{D}\text{r}\text{y} \text{B}\text{i}\text{o}\text{l}\text{o}\text{g}\text{i}\text{c}\text{a}\text{l} \text{Y}\text{i}\text{e}\text{l}\text{d} \left(\text{k}\text{g}/\text{h}\text{a}\right) }{100}\) Eq. 3
According to Pythagorean's Theorem for finding a triangle's hypotenuse, the mathematical definition of SNMI is as follows (Eq. 4) (Zhang & Davidson, 2022). The definition of SNMI based on normalized (NUE*) and N yield (NYield*).
SNMI = \(\sqrt{{\left(1-\text{N}\text{Y}\text{i}\text{l}\text{e}\text{d}\text{*}\right)}^{2}+{\left(1-\text{N}\text{U}\text{E}\text{*}\right)}^{2}}\) Eq. 4
Here, NYiled* was calculated from the Eq. 5 or 6
if NYieldref < NYield; then NYield* = 1 Eq. 6
Here, NYieldref = 90 kg ha− 1yr− 1 (Zhang & Davidson, 2022)
and if NUE ≤ NUEref ; then NUE* was estimated from Eq. 7. In this experiment, it was NUE ≤ NUEref.
The N response index (NRI) was calculated from Eq. 8 (Johnson and Raun 2003).
NRI = \(\frac{\text{M}\text{a}\text{x}\text{i}\text{m}\text{u}\text{m} {\text{Y}\text{i}\text{e}\text{l}\text{d}}_{\text{N} \text{F}\text{e}\text{r}\text{t}\text{i}\text{l}\text{i}\text{z}\text{e}\text{d} \text{c}\text{r}\text{o}\text{p}} \left(\text{t}\text{o}\text{n}/\text{h}\text{a}\right)}{{\text{Y}\text{i}\text{e}\text{l}\text{d}}_{\text{C}\text{o}\text{n}\text{r}\text{o}\text{l} \text{p}\text{l}\text{o}\text{t} \left(\text{t}\text{o}\text{n}/\text{h}\text{a}\right)}}\) Eq. 8
Kg N/ton yield was calculated as the Eq. 9, where kg N/ton yield represents the requirement of N (kg) to produce 1 ton yield.
Kg N/ton = \(\frac{\text{T}\text{o}\text{t}\text{a}\text{l} \text{N} \text{i}\text{n}\text{p}\text{u}\text{t} \text{b}\text{y} \text{f}\text{e}\text{r}\text{t}\text{i}\text{l}\text{i}\text{z}\text{e}\text{r} \left(\text{k}\text{g}/\text{h}\text{a}\right)}{\text{T}\text{o}\text{t}\text{a}\text{l} \text{y}\text{i}\text{e}\text{l}\text{d} \left(\text{t}\text{o}\text{n}/\text{h}\text{a}\right)}\) Eq. 9
2.4 N-use efficiency indices, N fertilization, and crop yield model to calculate N optimum
Statistically maximum NUE, grain yield, kg N/ton of yield, and minimum SNMI were taken into account and marked with red markers in an Excel scatter X, Y plot against N rate, where maximum NRI and minimum N rate were considered for all types of crops (potato, Boro, T. Aus, T. Aman, and maize) in order to develop the model to achieve a potentially optimal N rate. Data from Bogura and Mymensingh sites over the past two years were utilized for potatoes (n = 48) and maize (n = 48), respectively. However, for the Boro (n = 96) and T. Aus (n = 96) rice models, the average of two years data sets from two locations (Bogura and Cumilla) were taken into consideration, while for the T. Aman (n = 144)) rice model, the average of two years’ data from all three sites (Bogura, Cumilla, and Mymensingh) were utilized.
2.5 Soil analysis
Soil pH was determined (dry soil: water, 1:2.5) by electrode method (Jackson, 1973). Total organic C was determined in finely-ground (< 0.125 µg) soil by K2Cr2O7 and concentrated H2SO4 wet oxidation method (Walkley & Black 1934). Total N was measured in finely-ground (< 0.125 µg) soil by a micro-kjeldahl digestion and distillation method (Jackson, 1973). Soil available P was measured using the Olsen method (1954). Exchangeable K in soil was extracted with 1 M neutral ammonium acetate followed by ammonia distillation using the method described by Knudsen et al. (1983). Soil available S was determined by following Hunt (1980).
Bulk density was determined by standard core method using soil mass over volume ratio (Blake & Hartge, 1986). Soil particle size analysis was performed by the standard Bouyoucos hydrometer method after oxidizing the soil organic matter with 5% H2O2 followed by overnight dispersion with 10% Na-hexametaphosphate solution (Gee & Bauder, 1979).
2.6 Statistical Analysis
A two-way analysis of variance (ANOVA) was performed using residue management and N fertilization as fixed independent factors. The data distribution for normality was performed before the ANOVA test. The data were statistically examined to find out the variations between residue management and N fertilization rates simple and interactive effects. The Tukey-Kramer Multiple Comparisons Test was used in post-hoc analyses to distinguish the significant differences among main- and sub-plots variations and their interactions at p ≤ 0.05, unless otherwise mentioned. All the statistical analyses were performed using Statistix-10® and the graphic illustrations were created using Microsoft Excel.