The study was undertaken in two consecutive rainy seasons from 2014 to 2015 to 2015 to 2016 at a research farm belonging to the Indian Agricultural Research Institute (Pusa Institute), New Delhi (28.38° latitude and 77.09° longitude). The soil of the experimental unit had sandy clay loam texture in the upper 30-cm layer with a nearly level to gently sloping topography. The initial soil was slightly alkaline (pH 7.8), low in soil organic carbon (0.39%), low in available nitrogen (N) (155.4 kg ha⁻¹), high in available potassium (K) (311.4 kg ha⁻¹), and had a medium level of available phosphorus (P) (14.2 kg ha⁻¹). The climate of the study site is semi-arid to subtropical, with maximum temperatures ranging from 40 to 45°C during the cropping season. Rainfall of 390.8 mm (2014–2015) and 633.1 mm (2015–2016) was received during the crop growth period.

The study consisted of three crop establishment methods [Conventional (45 cm × 10 cm), SCI: 45 cm × 45 cm, and SCI: 30 cm × 30 cm] assisted in the horizontal strips, and four soybean varieties (Pusa-9712, PS-1347, DS-12-13, and DS-12-5) assisted in the vertical strips under a strip plot design with three replications. After obtaining a proper tilth in the soil, one plowing was carried out using a tractor-drawn double disc followed by harrowing. At the time of sowing, the recommended doses for nutrients for soybean (N, 25 kg; P₂O₅, 60 kg; and K₂O, 40 kg ha⁻¹), were applied through FYM (contained 0.50% N, 0.23% P₂O₅, and 0.56% K₂O) at 5.0 t/ha treated with Trichoderma (2.5 kg t⁻¹), and the remaining potassium and phosphorus doses were applied through muriate of potash and single super phosphate (SSP), respectively. The two sprouted seeds per hill were carefully sown at a spacing of 45 × 45 cm and 30 cm × 30 cm spacing without damaging the seed coat, and the respective seed rates of the desired row spacing were 15 and 25 kg ha⁻¹. Seed treatment was carried out according to the SCI protocol (Plate 1). During both study years, the soybean crop was established using the standard agronomic practices listed in Table 1.

Five soybean plants were randomly selected at intervals of 30, 60, and 90 days after planting (DAP) for the measurement of leaf area index (LAI), crop growth rate (CGR), and net assimilation rate (NAR) from the plant dry matter. Soybean was harvested manually using sickles during the second fortnight of October in both the years of study. After the harvesting bundles were left in the field for 2–3 days for drying, threshing was carried out using a Pullman thresher. Seed and straw yields were calculated from a 15.0 m² (5.0 × 3.0 m) plot at crop maturity. Yields were reported on a dry weight basis (12% moisture w/w), and the moisture content of the grain and straw was evaluated on an oven dry basis (70°C). The oil content of the soybean seed was estimated by using a grain analyzer (FOSS Infratec™ 1241) based on the technology of near-infrared transmittance using a non-destructive method of oil estimation. Leaf area (cm²) was measured with the help of a leaf area meter (LI-3100C) and was further converted into LAI using the following formula (Equation 1):

The computation of mean growth rate was performed using Equation 2, as provided by Watson (1952).

Where, W₁ and W₂ are the dry weight (g) of plants at time interval T₁ and T₂, respectively; S is land area (m²) covered by the plants.

The computation of the mean net assimilation rate was performed using Equation 3, as provided by Watson (1958):

Where,W₁ and W₂ are the dry weight (g) of plants at time T₁ and T₂, respectively; Ln is the natural logarithm; T₂- T₁ is the interval of time in days; LA₁ and LA₂ are the leaf area (m²) covered by plants at time T₁ and T₂, respectively.

Three plants were selected and pulled from the soil for root nodule examination at the 50% flowering stage of the soybean crop. The rhizosphere soil was then gently shaken from the root systems into a sterilized bag. The number of nodules per plant was then determined after carefully washing the roots with distilled water. Soybean roots were sampled at the 50% flowering stage and root samples were collected from the top 15 cm of soil using a root auger of 8.0 cm in diameter and 15 cm in length (core volume = 754.3 cm³). The roots were placed in a container with sieves to remove the soil debris and were then stored in a refrigerator at 4°C for preservation. Scanning and image analysis using the WIN-RHIZO system measured the root characteristics, such as root length and root volume (Costa et al., 2000).

An infrared gas analyzer (IRGA) was used to measure physiological parameters, such as photosynthetic rate (photo, mole CO₂ m² s⁻¹), stomatal conductance (cond, mole H₂O m² s⁻¹), transpiration rate (Tr, m. mol, H₂O m² s⁻¹), and intercellular CO₂ concentration (Ci, mole CO₂ mol⁻¹). The physiological indicators were measured at the flowering stage of soybean. Two CO₂ and H₂O infrared gas analyzers each are used in the LI-COR 6400, which provide estimates using the method described by Pandey et al. (2017). Photon (quantum) flux in radiant energy between 400 and 700 nm is used to describe photosynthetic active radiation (PAR). The PAR was determined using a LI-COR Line Quantum Sensor (1 m long) connected to a LI-1000 data logger. PAR was recorded at 1,200- and 1,300-h standard time on a sunny day.

Grain and straw samples were processed for nutrient analyses. First, samples were oven dried (70°C) until a constant weight was achieved followed by grinding using a Willey mill with stainless steel blades. A modified version of the Kjeldahl method was used to estimate nitrogen (%) in grain and straw (Prasad et al., 2006). Nitrogen (N) uptake was determined using Equations 4 and 5:

Grain and straw phosphorus (P) content was assessed using the vanado-molybdophosphoric acid yellow color method, as suggested by Prasad et al. (2006). Total P uptake (kg ha⁻¹) was determined using the following equations:

Grain and straw K content was estimated using a flame photometer (Prasad et al., 2006) and uptake (grain/straw and total) was computed using the following equations:

To estimate soil microbial dynamics, soil samples obtained from a depth of 0–15 cm were collected from each experimental plot using a core sampler. Collected soil samples were processed (air-dried, powdered, passed through a 2 mm mesh sieve, etc.) to estimate microbial parameters, such as soil microbial biomass carbon (MBC), enzymatic activities (dehydrogenase, alkaline phosphatase, and nitrogenase), and total polysaccharides. Soil MBC was evaluated using the method described by Nunan et al. (1998) with 70 mL of 0.5 M potassium sulfate (K₂SO₄). Dehydrogenase enzyme activity was assessed using the processes described by Casida et al. (1964) and a spectrophotometer at a wavelength of 485 nm. Alkaline phosphate activity, which represents free enzymes, was estimated using the technique suggested by Tabatabai and Bremner (1969) and 1.0 mL of p-ntro-phenyl phosphate disodium and 0.25 mL of toluene. Soil nitrogenase activity was estimated using the method suggested by Prasanna et al. (2003) and gas chromatography (GC) (Hewlett Packard 5890 series II) (using ethylene gas) and expressed as acetylene reducing activity (ARA). ARA was calculated using Equation (10):

Soil chlorophyll was examined using the techniques suggested by Nayak et al. (2004) with pre-weighed soil-cores (0–20 cm soil depth) and acetone: dimethyl sulfoxide (DMSO) in a 1:1 ratio with 4 m L g⁻¹ of soil. For soil chlorophyll determination, optical densities at 630, 645, 663, and 775 nm were taken (Equation 11).

The depth of irrigation water from each plot was measured using a meter scale. Crop water productivity (CWP) and water use efficiency (WUE) were calculated with the help of the formula published by Rajanna et al. (2019).

Soybean seed oil content was estimated using a grain analyzer (FOSS Infratec™ 1241) based on near-infrared transmittance technology using a non-destructive method of oil estimation. It is expressed as a percentage. Oil yield was calculated using Equation (15):

The research data recorded during the present study was evaluated using the “analysis of variance” (ANOVA) of the strip plot design (SPD) approach, as suggested by Gomez and Gomez (1984). The “F” test was used to determine the significance of treatments (variance ratio). When the “F” test revealed significant differences between means, differences in the treatment means were evaluated using Duncan’s multiple range test (DMRT) and least significant difference (LSD) at a 5% probability level. Standard error (SE) was used (in Figures 1–6) to measure the degree of variability between the individual data values.

Growth indices, such as LAI, CGR, and NAR, were significantly influenced by the different cultivation methods and soybean genotypes (Table 2). At 30, 60, and 90 days after sowing (DAS), the conventional method of soybean planting exhibited a significantly higher LAI under both systems of intensification methods (SCI). Similarly, with the conventional method of sowing at 45 × 10 cm, CGR and NAR were considerably greater than both SCI methods from the 0–30-day period to the 30–60-day period. However, for 60–90 days, CGR and NAR were significantly higher in SCI at a spacing of 30 × 30 cm compared with conventional methods and SCI at 45 cm × 45 cm. Among the genotypes, in all phases of development, DS 12–13 had a higher LAI than all the other tested genotypes, followed by DS 12–5 and Pusa 9712. Similarly, DS 12–13 soybean genotypes had the highest CGR across all the time periods (0–30, 30–60, and 60–90 days), while NAR was the highest for PS 1347 and Pusa 9,712 between 30–60 and 60–90 days, respectively.

Cultivation methods and soybean genotypes significantly influenced flowering-stage stomatal conductance, intercellular CO₂ concentration, and net photosynthetic rates during 2015–16 (Table 3). The highest intercellular CO₂ concentration among cultivation methods was produced by SCI at 30 cm × 30 cm (260.1 μmol CO₂ mol⁻¹), which was significantly greater than the conventional method of soybean cultivation. Likewise, SCI at 45 cm × 45 cm exhibited significantly higher stomatal conductance (0.211 mol H₂O m⁻² s⁻¹) and transpiration rate (7.8 m mol H₂O m⁻² s⁻¹) than other cultivation techniques, whereas the net photosynthetic rate was substantially higher with the conventional method (398 μmol CO₂ m⁻² s⁻¹) than with SCI techniques. Among the soybean genotypes, DS 12–13 had a higher stomatal conductance (0.22 mol H₂O m⁻² s⁻¹), intercellular CO₂ concentration (250.3 μ mol CO₂ mol⁻¹), and rate of transpiration (7.5 mmol H₂O m⁻² s⁻¹) than other soybean genotypes. There was no significant difference between DS 12–13 and DS 12–5. However, the Pusa 9712 genotype had a significantly higher rate of net photosynthesis (353 μmol CO₂ m⁻² s⁻¹) than the other genotypes.

Intercepted photosynthetically active radiation (PAR) of the soybean genotypes was also significantly influenced by the cultivation method and genotype. The adoption of SCI at 30 cm × 30 cm soybean establishment resulted in a significantly higher intercepted PAR (1,249 μ mol m⁻² s⁻¹) than other SCI methods (45 cm × 45 cm) and the conventional method of crop establishment. Increases were 14.1 and 8.5% with SCI at 30 cm × 30 cm compared with the conventional method and SCI at 45 cm × 45 cm, respectively. As a result, the PAR interception rate was 83.7% higher in SCI at 30 cm × 30 cm than the other methods. However, soybean genotypes have little impact/no significant effect on the interception of PAR. Comparatively, DS 12–13 intercepted more PAR than DS 12–5, PS 1347, and Pusa 9712.

Cultivation methods and genotypes significantly influenced nodules per plant⁻¹, total nodule dry weight, and root characteristics of soybeans at the blooming stage (Table 4). The SCI method of cultivation at 45 cm × 45 cm produced significantly more nodules per plant (41.1) and higher nodule dry weight (176.5 mg plant⁻¹) than the conventional method of cultivation (45 cm × 10 cm); however, the SCI method at 30 cm × 30 cm produced a similar number of nodules. Among the genotypes, DS 12–13 produced the maximum nodules per plant⁻¹ (35.8) and nodule dry weight (179.1 mg plant⁻¹), making it substantially more productive than the PS 1347 genotype. However, it was on a par with the Pusa 9712 and DS 12–5 genotypes. DS 12–13 had the most characteristics with the highest values, followed by Pusa 9712 and DS 12–5, all of which were comparable with one another.

Similarly, density of root length (1.58 cm cm⁻³), root volume (4.07 mm³ cm⁻³), surface area (0.118 cm² cm⁻³), and dry root mass (2.30 mg cm⁻³) were significantly higher in the SCI at 45 cm × 45 cm than the conventional method of cultivation. However, differences between both of the SCI methods were not significant. Among the genotypes, DS 12–13 produced a significantly higher density of root length (1.52 cm cm⁻³), root volume (3.66 mm³ cm⁻³), surface area (0.104 cm² cm⁻³), and dry root mass (2.30 mg cm⁻³) than Pusa 9712 and PS 1347 but it was on a par with DS 12–5.

During the 2014–2015 and 2015–2016 cropping seasons, the method of cultivation and varieties had a significant effect on soybean seed yields, as well as aboveground biomass yields and harvest index (Table 5). Among the cultivation methods, SCI at 30 cm × 30 cm had a considerably greater seed yield (1.99–2.13 t ha⁻¹), above-ground biomass yield (5.82–6.07 t ha⁻¹), and harvest index (34.2–35.1%) than the conventional method of cultivation (45 × 10 cm). The magnitude of increase in seed yields ranged from 9.6 to 13.3%, and the increases in biomass yields ranged from 8.2 to 10.7%. However, compared with the conventional method of cultivation, yields from both SCI methods were almost identical. Among the genotypes, DS 12–13 produced a significantly higher seed yield (1.97–2.15 t ha⁻¹) and aboveground biomass yield (5.86–6.15 t ha⁻¹) than the PS 1347 and Pusa 9,712 genotypes.

The interaction effect of crop establishment methods and genotypes significantly influenced soybean seed and biomass yields (Table 5). The combination of the DS 12–13 genotype with SCI at 30 cm × 30 cm produced a significantly higher seed yield (2.1–2.24 t ha⁻¹) and aboveground biomass yield (6.04–6.31 t ha⁻¹) than other combinations. Interestingly, all the studied genotypes produced significantly higher seed and biomass yields under an SCI at 30 cm × 30 cm than other SCI methods and conventional methods of cultivation.

Irrigation water productivity (IWP) and water usage efficiency (WUE) were considerably influenced by the different cultivation techniques and genotypes (Table 5). SCI at 30 cm × 30 cm yielded significantly higher WUE (3.06–3.34 kg ha⁻¹ mm⁻¹) and IWP (1.77–2.22 kg m⁻³) than conventional practice. Among the genotypes, DS 12–13 recorded significantly higher WUE (3.09–3.30 kg ha⁻¹ mm⁻¹) and IWP (1.64–2.01 kg m⁻³) than Pusa 9,712 and PS 1347. Among the interactions, the combination of DS 12–13 with an SCI at 30 cm × 30 cm produced significantly higher WUE (3.22–3.51 t ha⁻¹) and IWP (1.87–2.33 t ha⁻¹) than other combinations. All genotypes produced significantly higher WUE and IWP with an SCI at 30 cm × 30 cm.

Among the cultivation methods, an SCI at 30 cm × 30 cm had a considerably higher oil yield of 384.1 kg ha⁻¹ (2015–2016) and 413.2 kg ha⁻¹ (2014–2015) than the conventional method of cultivation (45 cm × 10 cm), whereas soybean genotypes significantly influenced soybean oil content and yield. Cultivation of the DS 12–13 soybean variety resulted in a statistically higher oil yield (394–430 kg ha⁻¹) than that of PS 1347 and Pusa 9712, but it was statistically on a par with DS 12–5, while soybean DS 12–5 produced a significantly higher oil content (20.4–20.5%) than other genotypes, but it was on a par with DS 12–13.

Perusal of data revealed that NPK concentration and their uptake in grain and straw of soybean influenced statistically due to cultivation methods during the study period (Table 6). However, the genotypes did not have any significant effect on NPK content. Statistically higher N content in grain (5.94%) and straw (2.61%) was noted in SCI at 45 cm × 45 cm compared with other methods of cultivation. Similarly, grain, straw, and total nitrogen uptake was statistically higher with SCI at 30 cm × 30 cm (108.1, 90.59, and 198.7 kg ha⁻¹, respectively) over the rest of the treatments. By contrast, crop establishment methods had no significant effect on P and K concentrations in soybean plant parts, whereas P and K uptake was significantly influenced by cultivation techniques and genotypes. SCI at 30 cm × 30 cm had a significantly greater P uptake in grain, straw, and total P uptake (9.38, 5.72, and 15.10 kg ha⁻¹, respectively) and K uptake in grain, straw, and total P uptake (15.5, 70.3, and 85.8 kg ha⁻¹, respectively) than other cultivation methods. Among the genotypes, the highest P uptake in grain and straw and total P uptake was recorded with DS 12–13 (9.33, 5.79, and 15.12 kg ha⁻¹, respectively), as was K uptake in grain and straw and total P uptake (15.6, 72.2, and 87.8 kg ha⁻¹, respectively) compared with the other genotypes.

Soil biological properties, such as dehydrogenase activity (DHA) (Figure 1), alkaline phosphatase activity (APA) (Figure 2), acetylene reduction assay (ARA) (Figure 3), total polysaccharides (Figure 4), microbial biomass carbon (MBC) (Figure 5), and soil chlorophyll (Figure 6), were significantly influenced by different methods of cultivation. Significantly higher soil DHA, APA, ARA, total polysaccharides, MBC, and soil chlorophyll were recorded in the SCI at 45 cm × 45 cm (13.63 μg TPF g⁻¹ soil hr.⁻¹, 93.2 μg p–nitro phenol g⁻¹ soil hr.⁻¹, 25.5 n moles ethylene g soil hr.⁻¹, 443.7 mg kg⁻¹ soil, 216.5 mg kg⁻¹ soil, and 0.43 mg g⁻¹ soil, respectively) compared with the conventional method. Among the genotypes, APA and ARA were considerably impacted by soybean genotypes over the study seasons. However, genotypes had no apparent impact on the levels of DHA, total polysaccharides, MBC, and soil chlorophyll. Among the soybean genotypes, DS 12–13 exhibited significantly higher APA (90.0 μg p–nitro phenol g⁻¹ soil hr.⁻¹) and ARA (18.8 n moles ethylene g soil hr.⁻¹) than the other genotypes.

A linear increase in growth indices, such as LAI, CGR, and NAR, was observed with the advancement in the growth stage of the soybean crop under the different cultivation methods and varieties. SCI either at 30 cm × 30 cm or 45 cm × 45 cm exhibited significantly higher LAI, CGR, and NAR over the conventional method (45 cm × 10 cm) due to less competition for growth factors, such as solar radiation interception, root growth, soil aeration, and microbial activity. Ball et al. (2000) reported that an optimum plant population is required for the maximum interception of light for a higher crop growth rate, LAI, NAR, and dry matter accumulation (Berger-Doyle et al., 2014) in the crop.

In the current study, growth indices were higher with the conventional method up to 60 DAS; however, with the advancement of growth stages, LAI, CGR, and NAR increased with the SCI at 30 cm × 30 cm compared with the other cultivation methods. Greater leaf area per plant was observed due to the higher growth performance of the plants that were more widely spaced. A higher population density with the conventional method caused the ceiling effect in the plant canopy, resulting in greater stem elongation leading to increased plant height (Rahman et al., 2004). Rehman et al. (2014) concluded that narrow row spacing (30–40 cm) significantly increased NAR, CGR, RGR, and LAI of soybean compared with wider spacing (70 cm). Among the genotypes, growth attributes were higher in DS 12–13 and DS 12–5 than in other tested cultivars. The growth attributes were significantly influenced by the genetic makeup of plants.

Leaf is the prime plant part to receive incident solar radiation (PAR) and transform it into photosynthates. Therefore, better soybean growth attributes resulted in better physiological parameters. In the current study, the highest intercellular CO₂ concentration was produced by an SCI at 30 cm × 30 cm (260.1 μmol CO₂ mol⁻¹), and the SCI at 45 cm × 45 cm exhibited a significantly higher stomatal conductance (0.211 mol H₂Om⁻² s⁻¹) and rate of transpiration (7.8 mmol H₂O m⁻² s⁻¹) than the other cultivation techniques. The net photosynthetic rate for soybean was substantially higher with the conventional method (398 μmol CO₂ m⁻² s⁻¹) than with the SCI techniques. A higher transpiration rate and stomatal conductance was observed in the SCI due to the greater leaf surface area and higher temperature (26–36°C). An increased transpiration rate leads to high stomatal conductance, thereby increasing CO₂ influx into the chloroplasts, possibly leading to a higher net photosynthetic rate. These findings are corroborated by Moreira et al. (2015) who reported that low planting densities and N rates increased the greenness index (SPAD value), concentration of intercellular carbon dioxide (Ci), and inherent WUE.

Among the soybean genotypes, DS 12–13 had the highest intercellular CO₂ concentration (250.3 μmol CO₂ mol⁻¹), rate of transpiration (7.5 mmol H₂O m⁻² s⁻¹), and stomatal conductance (0.220 mol H₂O m⁻² s⁻¹). There was no significant difference between DS 12–13 and DS 12–5. However, the Pusa 9,712 genotype had a significantly higher rate of net photosynthesis (353 μmol CO₂ m⁻² s⁻¹) than other genotypes. An increased transpiration rate leads to high stomatal conductance, thereby increasing CO₂ influx into the chloroplasts. A significantly higher chlorophyll content, net photosynthetic rate, and stomatal conductance in JS 95–60 than in JS 97–52 was also reported by Vyas and Khandwe (2014).

The PAR interception and transmission by soybean in SCI at 30 cm × 30 cm was 14.05 and 8.5%, respectively, and it was greater than the values obtained with the conventional method and the SCI at 45 cm × 45 cm (Table 5). PAR radiation was intercepted more with the SCI due to the green leaf coverage per unit area. Thereby, SCI at 30 cm × 30 cm intercepted 83.65% of incident PAR compared with 73.35% with the conventional method. High foliage surface aids a higher transpiration rate with a high maximum temperature (22–26°C) at the anthesis stage of soybean; furthermore, an increased transpiration rate leads to high stomatal conductance, thus increasing CO₂ entry into the chloroplasts, which might be the reason behind the higher net photosynthetic rate with the SCI method. The density of plants significantly influenced the economic yield in soybean by altering leaf area and therefore light interception and photosynthesis (Wells, 1991) and increasing plant spacing and decreasing row spacing, chlorophyll content (SPAD value), LAI, and the rate of photosynthesis (Jiang et al., 2015). With narrow plant spacing and wider row SPAD, LAI and Pn were higher. The use of wheat straw mulch and wider spacing by Dass and Bhattacharyya (2017) improved leaf SPAD values, PAR interception, the net photosynthetic rate, stomatal conductance, and WUE of soybean in dryland areas. There was a non-significant variation noted among all soybean cultivars with respect to intercepted PAR.

Wide space sowing decreases aboveground competition and ground competitiveness between roots for the absorption of water and nutrients and provides the adequate arena for nodulation and root growth of the crop plants. Therefore, in the current study, wider spacing at 45 cm × 45 cm with an SCI method of cultivation recorded statistically higher nodule counts, higher nodule dry weight (176.5 mg plant⁻¹), better root length, volume, and mass density than the conventional method of cultivation (45 cm × 10 cm). Seed inoculation using Rhizobium culture might have enhanced the population of the effective and healthy Rhizobium strain in the root nodules, which had a better ability to fix atmospheric N. Improved nodulation with the SCI was due to the slow release of nitrogen from FYM, which tends to form more nodules to meet the nitrogen demand of the crop (Salvagiotti et al., 2008). Mondal et al. (2014) showed that plant population had a significant effect on root length, root surface area, and the number of lateral roots/plant at a broader spacing of 20 cm × 30 cm followed by 15 cm × 30 cm.

Among the genotypes, DS 12–13 genotype produced a significantly higher density of root length (1.52 cm cm⁻³), root volume (3.66 mm³ cm⁻³), surface area (0.104 cm² cm⁻³), and dry root mass (2.30 mg cm⁻³) Pusa 9712 and PS 1347, but it was on a par with DS 12–5. Varietal root growth performance depends upon the cultivation practice. Additionally, leaf area/plant might be due to the higher growth of individual plants at a wider arrangement than the narrow arrangement. Imtiyaz et al. (2014) observed that DS 9712 had a greater reduction in root dry weight and fresh weight; however, the greatest reduction in shoot fresh and dry weight was noted in PS-1347 and DS-9712, respectively.

The above- and below-ground part of the soybean plant proliferated with the SCI because it provides broader spacing and square planting. Broader spacing leads to the improved growth of specific plants due to the availability of proper space, light, and moisture, which ultimately resulted in higher productivity of the crop (Dass et al., 2015). In the current study, SCI at 30 cm × 30 cm considerably increased seed yields from 9.6 to 13.3% and biomass yields from 8.2 to 10.7% compared with the conventional method of ‘cultivation and the SCI at 45 cm × 45 cm during 2014–2015 and 2015–2016. The increase in yield attributing parameters due to better plant growth and vigor brings out higher yield. More widely spaced soybean sowing resulted in higher seed and straw yield than with narrow row spacing (Pandya et al., 2005). However, compared with the conventional method of cultivation, the yields from both SCI methods were almost identical. Among the genotypes, DS 12–13 produced a significantly higher seed yield (1.97–2.15 t ha⁻¹) and aboveground biomass yield (5.86–6.15 t ha⁻¹) than PS 1347 and Pusa 9712 during both the study seasons. Different genotypes of soybean failed to affect the harvest index significantly during 2014–2015 and 2015–2016. The combination of the DS 12–13 genotype with the SCI at 30 cm × 30 cm produced a significantly higher seed yield (2.10–2.24 t ha⁻¹) and aboveground biomass yield (6.04–6.31 t ha⁻¹) than other combinations. Interestingly, all the studied genotypes produced significantly higher seed and biomass yields with the SCI at 30 cm × 30 cm than other SCI methods and the conventional method of cultivation. The SCI offers wider spacing, square planting, and microbe-loaded organic sources of nutrients, which provide better options for the development of aboveground parts of plants. A similar finding was reported by Abbas et al. (1994) who reported that an optimal plant population boosts the plant to improve yield and yield-attributing characteristics in plants (Abbas et al., 1994).

The SCI at 30 cm × 30 cm yielded significantly higher WUE (3.06–3.34 kg ha⁻¹ mm⁻¹) and IWP (1.77–2.22 kg m⁻³) than the conventional method (Table 5). Among the genotypes, DS 12–13 recorded significantly higher WUE (3.09–3.30 kg ha⁻¹ mm⁻¹) and IWP (1.64–2.01 kg m⁻³) than Pusa 9712 and PS 1347. Among the interactions, the combination of DS 12–13 with the SCI at 30 cm × 30 cm produced significantly higher WUE (3.22–3.51 t ha⁻¹) and IWP (1.87–2.33 t ha⁻¹) than the other combinations.

Oil yield is the function of its content and seed production. Among the cultivation methods, SCI at 30 cm × 30 cm had a considerably higher oil yield than the conventional method of cultivation (45 cm × 10 cm) and the SCI at 45 cm × 45 cm. Relatively higher oil yields of 19.38 and 19.27% were obtained with the SCI at 30 cm × 30 cm, followed by the SCI at 45 cm × 45 cm with 19.29 and 19.21%, compared with the conventional method during 2014 and 2015, respectively (Table 5). Improvements in oil content and oil yields with the SCI might be due to increased N and S content in grain on account of enhanced soil nitrogen and S through FYM and single super phosphate. Significantly higher oil content and yield were recorded in DS 12–5 (20.5 and 20.4%) and DS 12–13 (431 and 394 kg/ha) compared with other varieties in both years. A similar finding was reported by Moreira et al. (2015) who reported that plant population, row spacing, and nitrogen rates did not influence oil content in soybean seeds. Diep et al. (2016) reported that biofertilizer application in soybean significantly increased seed oil content and oil yield compared with the control, but remained on a par with 400 kg ha⁻¹ NPK 15-15-15.

Cultivation methods and genotypes significantly affected NPK content in grain and straw, and NPK uptake. Statistically higher NPK content in grain and straw was reported in SCI at 45 cm × 45 cm compared with other methods of cultivation. Similarly, NPK grain and straw and total uptake was statistically higher with SCI at 30 cm × 30 cm. The higher concentration of nutrients in SCI might be due to the microbial breakdown of organic manures, resulting in a release of nutrients into the soil (Prajapati et al., 2014) and, therefore, faster absorption and translocation into plants. Nutrient uptake is a function of nutrient concentration and biomass, and an increase in yield coupled with increased nutrient content resulted in a higher total uptake (Upadhyay et al., 2019) of nutrients with SCI at 30 × 30 cm. The greater microbial and enzymatic activities observed in the present investigation (Table 6) might also be responsible for the release of nutrients, ensuring an adequate supply to the plants. Among the genotypes, DS 12–13 had the highest NPK uptake in grain and straw and total P.

Soil biological properties, such as dehydrogenase activity (DHA) (Figure 1), alkaline phosphatase activity (APA) (Figure 2), ARA (Figure 3), total polysaccharides (Figure 4), MBC (Figure 5), and soil chlorophyll (Figure 6), were significantly higher with the SCI at 45 cm × 45 cm than with the conventional method of cultivation. Improved soil biological properties with the SCI might be due to the application of FYM and split application of vermicompost, which increased soil moisture availability, and the addition of nutrient and organic carbon through the decomposition of organic residue. The availability of easily decomposable carbon (C) sources in soil allows a significant relationship between microbial and enzyme activity and improves the soil water status (Shen et al., 2016). However, carbon is positively correlated with the microbial biomass carbon (Almeida et al., 2011). Hoeing practice in the SCI created a soil mulch, which reduced evaporation and improved soil moisture and, thus, enhanced soil enzymatic activities compared with the conventional method. Among the genotypes, APA and ARA were considerably impacted by soybean genotype. However, genotypes had no apparent impact on the levels of DHA, total polysaccharides, MBC, or soil chlorophyll. Thus, DS12-13 exhibited significantly higher APA (90.0 μg p–nitro phenol g⁻¹ soil hr.⁻¹) and ARA (18.8 n moles ethylene g soil hr.⁻¹) than the other genotypes.