Spatial Arrangement and Biofertilizers Enhance the Performance of Legume—Millet Intercropping System in Rainfed Areas of Southern India

Intercropping is a well-established practice to enhance the yield in low-input agriculture, and beneficial microbes such as arbuscular mycorrhizal fungi (AMF) combined with plant growth promoting rhizobacteria are being used as an effective and sustainable measure to improve yields. In this study, we tested if biofertilizers can not only enhance the yield of crops in monoculture as has previously been demonstrated but can also enhance the yield of intercropping systems. We hypothesized that because AMF can form common mycorrhizal networks (CMN) that can transfer nutrients and water between different plant species, biofertilization can balance belowground competition between crop species and promote thus overall yields in intercropping systems. In our study, we used a pigeon pea (PP)—finger millet (FM) intercropping system that we grew for two consecutive growing seasons (2016/17 and 2017/18) at two contrasting sites in Bengaluru and Kolli Hills, India. We also tested if the spatial arrangement (i.e., different arrangement of component plants with similar plant density in intercropping system) of intercropped plants, using either a row-wise or a mosaic design, influences the effect of biofertilizers on yield and water relations of the PP-FM intercropping system. Our results demonstrate that intercropping can improve the straw and grain yield of PP and FM compared to the respective monocultures and that intercropping effects vary depending on the site characteristic such as climate and soil type. The spatial arrangement of component plants affected the total, straw, and grain biomass in intercropping treatments, but this effect also varied across sites. Most importantly, the results from the 2017/18 growing season clearly demonstrated a positive effect of biofertilizer on biomass yield, and this effect was irrespective of site, spatial arrangement, mixed or monoculture. Our study therefore shows that yield increase in intercropping systems can further be improved through the application of biofertilizers.

Intercropping is a well-established practice to enhance the yield in low-input agriculture, and beneficial microbes such as arbuscular mycorrhizal fungi (AMF) combined with plant growth promoting rhizobacteria are being used as an effective and sustainable measure to improve yields. In this study, we tested if biofertilizers can not only enhance the yield of crops in monoculture as has previously been demonstrated but can also enhance the yield of intercropping systems. We hypothesized that because AMF can form common mycorrhizal networks (CMN) that can transfer nutrients and water between different plant species, biofertilization can balance belowground competition between crop species and promote thus overall yields in intercropping systems. In our study, we used a pigeon pea (PP)-finger millet (FM) intercropping system that we grew for two consecutive growing seasons (2016/17 and 2017/18) at two contrasting sites in Bengaluru and Kolli Hills, India. We also tested if the spatial arrangement (i.e., different arrangement of component plants with similar plant density in intercropping system) of intercropped plants, using either a row-wise or a mosaic design, influences the effect of biofertilizers on yield and water relations of the PP-FM intercropping system. Our results demonstrate that intercropping can improve the straw and grain yield of PP and FM compared to the respective monocultures and that intercropping effects vary depending on the site characteristic such as climate and soil type. The spatial arrangement of component plants affected the total, straw, and grain biomass in intercropping treatments, but this effect also varied across sites. Most importantly, the results from the 2017/18 growing season clearly demonstrated a positive effect of biofertilizer on biomass

INTRODUCTION
Intercropping has been considered a sustainable way to utilize and share natural resources among different crop species and to improve and stabilize crop yield (Brooker et al., 2015;Martin-Guay et al., 2018). In intercropping systems two or more crop species are grown together (Vandermeer, 1989). Crop yield in intercropping systems are often higher than in monocropping systems because resources such as soil moisture and nutrients are utilized more efficiently in intercropping than in monocropping (Lithourgidis et al., 2007;Dahmardeh et al., 2009;Martin-Guay et al., 2018). This is because interspecific competition between intercropping partners is often lower than the intraspecific competition in monocropping, resulting in a yield advantage (Davis and Woolley, 1993). In addition, beneficial effects of intercropping can come from resource facilitation: As an example, legume-cereal intercropping systems have been widely used in areas with poor soil quality (Li et al., 2007), where legumes fix nitrogen (N) and solubilize phosphorus (P) which is then used by both intercropping partners (Hinsinger et al., 2011). In return, cereals can support legumes in two ways, by preventing nitrate-N accumulation in soil which inhibits N fixation by legumes, and by increasing iron availability which enhances N fixation (Zuo et al., 2004;Schipanski and Drinkwater, 2012).
In rainfed areas of the arid and semiarid tropics, intercropping has also been suggested to enhance the water availability of shallow-rooted crops via the facilitation of water supply by deeprooted plants through hydraulic lift (HL) (Xu et al., 2008;Mao et al., 2012). The water released from deep-rooted plants due to HL into topsoil layer becomes available to neighboring shallowrooted plants, a process termed "bioirrigation" (Burgess, 2011). The functionality of bioirrigation in intercropping systems has only been tested in a few studies-mainly under controlled conditions in the greenhouse. Sekiya and Yano (2002) showed in a field experiment that pigeon pea (a deep-rooted legume) has the potential to perform HL and could supply deep water to shallowrooted maize. In another study, Sekiya et al. (2011) showed that plants with deep roots are ideal for intercropping with shallowrooted crops in water limited agricultural fields and that this kind of intercropping system allows shallow-rooted plant to access deep soil moisture indirectly without having deep roots. Other studies have also shown the transfer of hydraulically lifted water (HLW) from a deep-rooted plant to neighboring shallow-rooted plants (Caldwell and Richards, 1989;Moreira et al., 2003;Brooks et al., 2006;Bogie et al., 2018). While these experiments have suggested that bioirrigation could be an important mechanism for drought stress avoidance of intercropped field crops, evidence for the efficiency of this mechanism in the field is still lacking.
The success of an intercropping system in the field depends on the avoidance of competitive growth inhibition among the intercropping partners. This requires appropriate spacing of the intercropping partners so that competitive, complementary, and facilitative interactions are well balanced and that yield improvements can be achieved. In particular, for bioirrigation to be effective, it seems that an ideal spacing between the intercropping partners is essential. On the one side, intercropping partners have to be arranged with sufficient space among each other in order to avoid competition. On the other side, plants need to be spaced (Burgess, 2011;Prieto et al., 2011) to allow the transfer of bioirrigated water from one rhizosphere to another.
In addition to intercropping approaches, "biofertilization" such as inoculation with arbuscular mycorrhizal fungi (AMF), combined with plant growth promoting rhizobacteria (PGPR), are beginning to become established as an effective and sustainable measure to improve yields (Mäder et al., 2011;Schütz et al., 2018;Mathimaran et al., 2020). The role of AMF for the uptake and transfer of nutrients and water to host plants has been well demonstrated (Augé et al., 2001;Querejeta et al., 2003). Biofertilization might have a particular potential to boost the yield of intercropping systems because AMF can form a common mycorrhizal network (CMN) that can transfer nutrients between two plants and balance as such belowground competition (Smith and Read, 2008). In addition, a CMN between the roots of two plants can also constitute a pathway for the transfer of water. Egerton-Warburton et al. (2007) have demonstrated that arbuscular mycorrhizal hyphae provide a potential pathway for the transfer of HLW between two plants. Our recent work has shown that CMN plays a key role in facilitating the transfer of water between the rhizospheres of two intercropped partners and can in turn improve the water relations of shallow rooted crops during soil drying in a greenhouse experimental setup under controlled conditions (Singh et al., 2019). However, a further experiment with bigger pots (50 L) than in the previous experiment did not show an effect of the CMN on water relations but treatments with CMN had lower foliar damage than treatments without CMN during drought (Singh et al., 2020).
The effects of biofertilizers on stabilizing and improving the yields in intercropping systems by improving water relations via bioirrigation have not yet been tested under field conditions, though recent greenhouse studies have shown evidence of facilitation of bioirrigation by AMF and PGPRs (Saharan et al., 2018;Singh et al., 2019). Furthermore, it is unclear to what extent beneficial effects of biofertilizers in intercropping systems depend on an appropriate spacing of the crops and if-given the appropriate spatial arrangement of crops-the establishment of a CMN can indeed facilitate bioirrigation and improve the water relations of shallow-rooted crops in intercropping systems in dryland agriculture. In this study, we investigated the effects of biofertilization on the yield of a legume-cereal intercropping system and tested different spatial arrangements of the plants in combination with biofertilizer treatments. We used pigeon pea (Cajanus cajan) (PP) as a deep-rooted plant and finger millet (Eleusine coracana) (FM) as shallow-rooted plant to investigate the following research questions: (i) Does the spatial arrangement (i.e., different geometric arrangement of component crops using row-wise and mosaic pattern) of intercropping partners affect straw and grain yield in a FM-PP intercropping system compared to monocropping of the same crops? (ii) Does the application of biofertilizers further improve yields in spatially differently arranged intercropping systems? (iii) Can intercropping in conjunction with a CMN lead to an improvement of the water relations of shallow-rooted crops?

Selection of Field Experiment Site and Crop Varieties
To test the influence of spatial arrangement and biofertilizers on crop yields of PP and FM, field trials were carried out at two different locations during the growing seasons 2016/17 and 2017/18. One experimental site was located at the research field of the University of Agricultural Sciences, Gandhi Krishi Vigyana Kendra Campus (GKVK), Bengaluru, Karnataka, India. The other site was located at the research field of M S Swaminathan Research Foundation (MSSRF), Kolli Hills, Tamil Nadu, India. At both sites two field trials were conducted at different plots; at GKVK, coordinates for 2016/17 and 2017/18 are 13 • 05' 14" N 77 • 34' 13" E and 13 • 05' 13" N 77 • 34' 15" E, respectively. For the Kolli Hills site, coordinates for 2016/17 and 2017/18 are 11 • 16' 36.0" N 78 • 23' 58.6" and 11 • 16' 34.0" N 78 • 23' 58.6", respectively. Both experimental sites were selected because farmers have already adopted a cereal-legume intercropping system there and have been cultivating PP and FM as one of their main crops, and to cover broad growth conditions as both sites differ in their soil characteristics and average annual rainfall. The soil type at the Bengaluru site is an alfisol consisting of 67.8% sand, 7.7% silt and 25.2% clay with pH of 5.1, while Kolli Hills has a vertisol soil type with 33.2% sand, 30% silt and 36.8% clay with pH of 5.2. Based on farmers practice in the region and recommendations from local agronomists, we selected FM (GPU-28) and PP (BRG-2) for the field experiment at the Bengaluru site and PP (SA-1) and FM (Suruttai kelvaragu) at the Kolli Hills site .

Rainfall
The total annual precipitation at the Bengaluru site was 694 mm in 2016 and 1,104 mm in the year 2017. At Kolli Hills, the total annual precipitation was 281 mm in 2016 and 1,690 mm in 2017. Rainfall data recorded during the experimental period indicate that the Kolli Hills site received less rain than the Bengaluru site (Supplementary Figure 1). Both sites received the maximum amount of rain during the months of May, June and July. The Bengaluru site received up to 40-60 mm rain during September, October, and December, while the Kolli Hills site was completely dry after July during 2016. During the summer 2016, both research sites, received very low precipitation, as recorded by our weather station at the field sites. For the comparison and visulaization, precipitation data from nearest sites as recorded by the Climate Research Unit (Harris et al., 2020)

Intercrop Field Design With Different Spatial Arrangement of PP and FM
The plot size for each treatment was 7.2 × 3.6 m (width × length) with a net plot area of 3.6 × 1.8 m (Figure 1). The net plot area defines the central part of each plot, where all physiological, growth and yield parameters were assessed. The field experiments had six treatments: FM monoculture (T1), PP monoculture (T2), 2:8 (PP:FM) row-wise intercropping (T3), 1:4 (PP:FM) rowwise intercropping (T4), 100% mosaic (T5) and 50% mosaic (T6) ( Figure 1B). Each treatment was replicated four times in a randomized block design.
In monocultures, the density of FM was 48 plants per m 2 , and the density of PP was 6 plants per m 2 . We planted 8 times more individuals for FM than for PP per area and the total number of plants for FM in monoculture (T1) was 1,152 per plot and 288 plants in the net plot area. For PP monocrop (T2), the total number of plants was 144 in the total plot and 36 plants in the net plot area. The spacing between FM rows was 30 cm and the distance between FM plants within a row was 7.5 cm. The spacing between PP rows was 60 cm and the distance between PP plants within a row was 30 cm. In intercropping treatments, the spacing between PP and FM rows was 45 cm.
Intercropping systems were based on FM monocultures, where eight FM plants were substituted by one PP plant. Rowwise intercropping systems (treatment T3 and T4) were based on previous investigations under rainfed conditions in Karnataka, India (Ashok et al., 2010;Padhi et al., 2010;Mathimaran et al., 2020). For T3 (2:8 PP:FM row-wise arrangement), each replicate had thus 48 PP (12 plants × 4 rows) and 768 FM (48 plants ×12 rows in each total plot area). T4 (1:4 PP:FM row-wise arrangement) had the identical number of PP and FM plants as T3 but it differed in row arrangement where one row of PP was planted after four rows of FM. Treatment T5 (100% mosaic) consisted of identical numbers of PP and FM plants as T3 and T4, but PP and FM plants were planted within the same row in a mosaic design ( Figure 1B). In treatment T6 (50% mosaic), the number of PP was reduced by 50 % and replaced by FM plants. It consisted of 24 PP plants (2 plants × 12 rows) and 960 FM plants. In the 2017/18 field trial at the Bengaluru site, FM plants in T5 were not substituted by PP but PP was accidentally added into the mosaic design. Therefore, plant density of FM was higher than in the other treatments.
We established the same treatments in the years 2016/17 and 2017/18 except for T6, which was not established in 2017/18 based on results from 2016/17 field trial. While field trials during year 2016/17 had only treatments with biofertilizers, field trials during the year 2017/2018 included treatments with and without biofertilizers ( Table 1).

Application of Fertilizers and Biofertilizers
Mineral fertilizers were applied at 50% of the recommended dose of fertilizer (RDF) to farmers for all plots as basal application, i.e., at rate of 12.5:25:12.5 NPK kg ha −1 for PP and 25:20:12.5 NPK kg ha −1 for FM. Nitrogen (N) fertilizer was given in the form of Urea (46% N-0P 2 O 5 -0K 2 O, SPIC India Fertilizer Company), Phosphate (P) fertilizer was given in the form of Single Super Phosphate (SSP, 0N-16% P 2 O 5 -0K 2 O, SPIC India Fertilizer Company), and Potash (K) fertilizer was given in the form of Muriate of Potash (MOP, 0N-0P 2 O 5 -60% K 2 O, SPIC India Fertilizer Company). The remaining 50% of the fertilizer was not added as we wanted to test the effectiveness of biofertilizer at only 50% of RDF, with the objective to minimize the input cost of farmers and potentially achieve long-term soil sustainability by reducing the use of mineral fertilizers.
Biofertilizers PGPRs (Pseudomonas sp. MSSRFD41 and Rhizobium liquid formulation) and AMF (Rhizophagus fasciculatus and Ambispora leptoticha) spores in a vermiculite substrate were applied to PP and FM respective treatments as described by Mathimaran et al. (2020). In the biofertilizer treatments, FM seeds coated with PGPR Pseudomonas sp. MSSRFD41 consisting of 1x 10 9 CFU per ml at the rate of 10 ml kg −1 seed and the PP seeds coated with both PGPR Pseudomonas sp. MSSRFD41 and Rhizobium consisting of 1x 10 9 CFU per ml at the rate of 10 ml kg −1 seed were sowed, and for no biofertilizer treatments untreated FM and PP were used. AMF A. leptoticha was applied at the rate of 5 g per PP seedling hole (germinated in a polybag, see below) with a dosage of 278 kg ha −1 . Similarly, in the rows of FM R. fasciculatus, consisting of 15 spores g −1 of substrate was applied as a band application at the rate of ca. 444 kg ha −1 . Additionally, a band application of Pseudomonas sp. MSSRFD41 were applied in the planting rows of biofertilizer treatment plots at the rate of 49.5 L ha −1 with a carrier of 7.5 t ha −1 farmyard manure (FYM) and the treatment details were mentioned in the table 1.0 and described in our earlier work Mathimaran et al. (2020). The PGPR strains were multiplied in a liquid formulation of King's B medium.

Pre-germination, Sowing of Seeds Into Field, Growth Period and Harvest
Based on an established practice in the area, PP seeds were pre-germinated in polybags (15 × 10 cm) filled with 1.6 kg of a mixture of field soil: FYM:sand (ratio of 15:1:1), and a seed hole of 4 × 1 cm was made at the top . The bottom layer of the seed hole was filled with A. leptoticha in vermiculite, two PP seeds coated with rhizobia and PGPR strains were kept above the vermiculite layer and field soil was filled on the top. The seeds were allowed to germinate and grow for 35-45 days. Later, healthy seedlings from these polybags were transplanted into the field on the third week of July 2016

Growth and Yield Parameters
Plant growth parameters such as plant height, number of pods, pod weight per plant, number of panicles, grain weight per panicle, straw and grain biomass (both sun dried and oven dried), weight of 1,000 FM seeds and 100 seeds of PP were measured after harvesting the plant material in the net plot area. For biomass, plants were harvested row-wise in the net plot area and straw and grains were separated. The sun-dried biomass was determined after drying the straw under the sun for 15 and 20 days for FM and PP, respectively. Grains were dried under sun for 10 days for PP and FM. A subsample of the sun-dried straw and grain material was oven dried at 80 • C for 24 h for calculating the dry matter per row, moisture content of oven dried samples were not measured. Biomass per plant was calculated by dividing the row biomass by the number of plants in each row; biomass in ton per ha was obtained by multiplying the row biomass with the number of rows per ha.

Land Equivalent Ratio (LER)
The facilitative and competitive interactions between PP and FM in response to the different treatments were calculated using the LER. The LER indicates the efficacy of an intercropping system for using natural resources compared with monoculture (Willey and Osiru, 1972). The baseline for LER is one. If the LER is greater than one, intercropping favors growth and yield of plants, and when it is lower than one, intercropping negatively affects the growth and yield of plants. The LER was calculated as Where Y FM and Y PP are yield of PP and FM in its monoculture, Y FM,PP is yield of finger millet in intercropping, and Y PP,FM is yield of pigeon pea in intercropping.

Measurement of Physiological Parameters
To test if different spatial arrangements of FM and PP, and the application of biofertilizers affect the water relations and growth of FM, we determined FM leaf water potential at predawn (04:00 to 05:00 h) and mid-day (12:30 to 13:30 h) toward the end of the field trial during first three weeks in November during 2016/17 and 2017/18. Due to limitation in resources, particularly manpower and time, these measurements were only performed at the Bengaluru site. Leaf water potential (LWP) was measured using a pressure chamber (model 1000, Pressure Chamber Instrument Company, USA). For predawn measurements, leaf samples were collected between 04:00 and 05:00 h and for midday measurements, leaves were sampled between 12:30 and 13:30 h. After sampling, leaves were packed into airtight Ziploc bags to avoid water loss. Bags were kept in the dark and leaf water potential was measured within 1-2 h after sampling.

Statistical Analysis
Analysis of yield data and LWP from field trials was carried out using GraphPad Prism software (version 7.0 for Mac OS X, GraphPad Software, La Jolla California USA). Data are expressed as mean ± standard error of mean (SEM). Tukey's test was used for post hoc multiple treatment comparison following one-way ANOVA or multifactor ANOVA using general linear models. The criterion for significance was p < 0.05.

Total Biomass, Straw, and Grain Yield per Hectare and LER
Intercropping and the spatial arrangement of the intercropping partners had a significant effect on the total biomass yield per hectare at the Bengaluru site in 2016/17 (Figure 2  (Supplementary Table 1). In summary, in 2016/17 we found a strong positive intercropping effect for total biomass yield, straw yield and grain yield at the Bengaluru site, where the intercropping effect were strongest in the 8:2 row-wise spacing.
In contrast, no yield improvements by intercropping irrespective of the spatial arrangement were observed at the Kolli Hills site. These observations are also reflected in LER values. At the Bengaluru site, LER values for total biomass were greater than one for T3+, T4+ and T5+ and where T3+ had the highest LER value. Similarly, for straw biomass, T3+ had higher LER values than T4+, T5+ and T6+. For grain biomass LER values were greater than one for the T3+ and T4+ treatment, equal to one for T5+ and less than one for T6+ (Figure 3). At Kolli Hills LER values for all treatments were less than one (Figure 3).
In 2017/18, intercropping and the spatial arrangement of the intercropping partners also had a strong and significant effect on the total biomass yield, straw yield, and grain yield at the Bengaluru site (Figure 4). As in 2016/17 the treatment T3-and T3+ produced significantly more biomass per hectare than monocultures of the constitutive crops or other spatial arrangements when compared to the respective treatments with and without biofertilizer. Importantly, the application of biofertilizers enhanced the total biomass yield, straw yield and grain yield in all treatments, and this effect was consistent irrespective of experiment site, monoor intercropping (Supplementary Table 2). At Kolli Hills, we also found significant treatment effects (Figure 4). However, intercropping treatments did not produce higher yields for total biomass and straw than any of the other treatments with or without biofertilizer. Yet, treatment T5+ was equal in total biomass yield to the most productive monoculture (T2+). For grain yield FM monoculture exceeded the productivity of PP ( Figure 4F) and in intercropping T3-, T3+ and T5+ grain yield was similar to monoculture of FM with or without biofertilizer. The effects of biofertilizers on total biomass yield, straw yield, and grain yield that we detected at the Bengaluru site were also observed at the Kolli Hills site and this effect was again consistent across all treatments (Figure 4; Supplementary Table 2). We did not find a significant interaction between treatment and biofertilizers nor a significant three-way interaction between treatment, biofertilizers, and site. However, as indicated above, the effects of biofertilizers at Kolli Hills resulted in total biomass yield, straw yield and grain yield that were of the same magnitude in some intercropping treatments as the highest yield in the corresponding monocultures (e.g., T5+ for total biomass yield and straw yield, and T3+ and T5+ for grain yield, Figure 4). In summary, in 2017/18 we found a strong positive intercropping effect for total biomass yield, straw yield and grain yield at the Bengaluru site. In Kolli Hills, no such intercropping effect was found. Importantly, biofertilizers improved the yields of crops at both sites and independently of treatment. Despite the nonsignificant biofertilization-treatment interaction, intercropping treatments at Kolli Hills yet showed a trend to be more enhanced through biofertilizers than monocultures, to an extent that they produced similar yields as the most productive monoculture, which we did not observe without biofertilizers.
These observations were confirmed by LER values for 2017/18 at both sites ( Figure 5). LER was greater than one at the Bengaluru site for all treatments. Also, LER values at the Bengaluru site were largest for T3 and declined in the other treatments. Biofertilizers had a negative effect on LER values in all spatial arrangements at the Bengaluru site. At Kolli Hills, LER values in treatments without biofertilizers were either equal to or less than one. Biofertilizers increased, however, the LER values in all spatial arrangements to values of one or greater than one and the largest values were observed for T3+ and T5+.

Per Plant Biomass Yield of PP and FM
We found a significant effect of the intercropping treatments on total biomass per plant, total straw yield per plant and total grain yield per plant of PP and FM at the Bengaluru site but not at Kolli Hills in 2016/17 (Figure 6; Supplementary Tables 3,  4). At Bengaluru, total biomass per plant in FM was highest in the monoculture (T1+), the 2:8 treatment (T3+) and the 1:4 treatment (T4+). The biomass of individual plant was significantly reduced in the mosaic treatments (T5+ and T6+) compared to monoculture (T1+) and row-wise intercropping (T3+ and T4+, Figure 6A). PP showed highest total biomass in the mosaic treatment T6+, followed by other intercropping treatments and lowest biomass in the monoculture T2+ (Figure 6C). At Kolli Hills, total biomass per plant in PP and FM did not differ significantly among treatments (Figures 6B,D). However, the trend was similar to the Bengaluru site where FM showed a reduction in biomass in mosaic treatments while PP showed an increase in biomass in mosaic treatments.
In 2017/18, we also found a significant treatment effect on the total biomass, straw yield and grain yield of FM and PP at the Bengaluru site but only for PP at Kolli Hills (Figure 7;  Supplementary Tables 5, 6). At the Bengaluru site, total biomass of FM plants in T3+ was significantly larger than total biomass of plants in treatments T1-, T1+ and T4-. Total biomass of PP plants was largest in T3+ and T4+ compared to T2-, T2+, and T4- (Figure 7C). At Kolli Hills total biomass per plant in FM did not show any significant difference among intercropping and monoculture. For PP, in contrast, total biomass per plant was largest in treatments T4+ and T5+ compared to T2-and T2+ ( Figure 7D).
A two-way ANOVA analysis was performed to test the effects of spatial arrangement and biofertilization on per plant yield (Supplementary Tables 7, 8). At both sites in 2017/18 FM yield did not show any significant effect of biofertilizer application. However, PP showed a strong significant effect of biofertilization at the Bengaluru site, and at the Kolli Hills site the effect was marginally significant. At both sites, the effect of biofertilization did not differ among treatments with different spatial arrangements of the component plants of the intercropping system.

Water Relations of PP and FM in Intercropping Treatments
Measurements of predawn leaf water potential (LWP) were done on FM leaves at the Bengaluru site to evaluate the effects of spatial arrangement and biofertilizer application on In 2017/18 at week 1 (1st week of November 2017), predawn LWP of FM in monoculture with biofertilizer (T1+) had values of −0.32 MPa which is significantly more positive than FM in monoculture without biofertilizer (T1-) (−0.60 MPa) or any of the intercropping treatments ( Figure 8B). Later, FM did not show any significant difference in LWP compared to the other intercropping treatments. Interestingly, treatments without biofertilizer showed generally lower values for LWP as compared to the respective treatments with biofertilizer. The biofertilizer application did not have a significant effect on LWP of FM, but intercropping treatments showed a strong significant effect (Supplementary Table 9). There was a significant interaction between the effect of biofertilizers and the intercropping treatments. As shown in Figure 8, treatments T1-, T1+, T5-, and T5+ consistently showed a large difference in LWP of FM with or without biofertilizer.

DISCUSSION
The results obtained from the field trials during 2016/17 and 2017/18 showed that intercropping can improve the straw and grain yield in PP-FM intercropping compared to the respective monocultures but that intercropping effects vary depending on the site characteristic such as climate and soil type as well as crop variety. The spatial arrangement of component plants also affected the total, straw, and grain biomass in intercropping treatments, and this effect also varied across sites. The results from 2017/18 clearly demonstrated a positive effect of biofertilizer on biomass yield, and this effect was irrespective of site, spatial arrangement, mixed or monoculture. Despite the positive effect of intercropping and biofertilization on FM and PP yield, water relations of FM were not enhanced in the intercropping treatments or by biofertilizers. Most likely this was due to interspecific competition between PP and FM for soil moisture in the topsoil layer. On the basis of these results, we propose that intercropping and the application of biofertilizer both enhance the yield of cropping systems and effects on yield are strongest if intercropping and biofertilization are applied in combination. However, the spatial arrangement of component crops is a key factor that affects the productivity of the involved intercropping partners and needs to be considered in the design of intercropping systems.
Is PP-FM Intercropping Beneficial Over Monocropping?
The yield advantage in intercropping systems is typically assigned to resource sharing and facilitation (Loreau and Hector, 2001;Li et al., 2014;Duchene et al., 2017). Resource complementarity reduces the niche overlap and competition between two species and allows crops to take up a greater amount of resources than the sole crops. Ghanbari et al. (2010) reported resource complementarity in maize-cowpea intercropping systems, where intercropping increased the light interception, reduced evaporation, and improved soil moisture conservation compared to maize sole crops. In most cases, facilitation occurs through increased availability of soil resources such as water and nutrients (Jensen et al., 2020). Intercropping systems with legume species (such as PP in this study) can increase agricultural productivity by providing increased nitrogen availability through N 2 fixation, and are therefore used very frequently (Hauggaard-Nielsen and Jensen, 2005;Altieri et al., 2012). Non-legumes (such as cereals) in an intercropping system with legumes may obtain additional N released by legumes into the soil, which can contribute up to 15% of their N content (Zuo et al., 2004;Li et al., 2007).
According to our expectation, we found that at the Bangalore research site, the intercropping treatments (T3+ and T4+) produced higher yields (Figures 2, 4) than monocultures in both growing seasons. In contrast, at Kolli Hills, there was no significant effect of intercropping in the 2016/17 season. Also in 2017/18, we did not observe strong intercropping effects but yields in some intercropping treatments (T5+) were as high as the highest yields in the monocrop. Accordingly, LER values were above one at Bangalore in both years but below one at Kolli Hills in 2016/17 and near zero or above in 2017/18. This illustrates that intercropping effects depend on the weather of the growing season and soil type at the experiment site. The total rainfall at Kolli Hills in 2017 was 1690 mm, while the 2016 growing season was shaped by a severe drought with a total rainfall of only 281.7 mm. Additionally, both locations differ in their soil properties. At the Bangalore site, the soil is an Alfisol with 67.8 % sand, 7.7 % silt, 25.2 % clay, C org 0.5 % and a pH of 4.8. At the Kolli Hills site, the soil type is a Vertisol with 33.2 % sand, 30.0 % silt, 36.8 % clay, C org 0.8 % and a pH of 5.2 (see Mathimaran et al., 2020, supplimentary data). The relationship between crop yield and soil depends on complex interactions between physico-chemical properties of soil and other climatic factors (Stenberg, 1998). Juhos et al. (2015), using a multivariate statistical approach, show that in years of drought, the sodification, salinization, soil texture, and nutrient content determined the yield, while in humid years soil organic matter and nutrient content were the main determining factors for crop yields. Our results indicate that the low amount of rainfall and inherent soil properties could have caused different intercropping effects at the two sites and between the two growing seasons (Figure 2).

Effect of Spatial Arrangement on Yield in PP-FM Intercropping
At the Bengaluru site, straw and grain yield (per hectare) was higher in the row-wise than in the mosaic intercropping treatments during 2016/17 and 2017/18 field trial. The results from Kolli Hills were inconsistent, perhaps because of rainfall and soil properties as indicated above. Effects of the spatial arrangements can be explained by intra-and inter-specific competition, as illustrated when data are expressed per plant biomass (Figures 6, 7). Results from Bengaluru site clearly indicate that PP benefits in terms of per plant biomass in intercropping treatments likely due to reduction in intra-specific competition that PP faces in monoculture. In contrast, FM faces higher inter-specific competition in mosaic treatments, which leads to a reduction in per plant biomass in mosaic treatments (T5+ and T6+). The field trial results from Kolli Hills, however, did not show any significant effect of spatial arrangement on per plant biomass in PP and FM during 2016/17 trial (Figures 6B,D). During 2017/18, only PP showed a significant increase in per plant biomass in intercropping treatments T4+ and T5+ compared to monoculture treatments T2-and T2+. The effect of spatial arrangement on FM per plant biomass was not significant and it was consistent during both years at Kolli Hills.
The results of this study consistently show that PP growth is favored in intercropping systems due to a reduction in the intraspecific competition, while FM faces a higher interspecific competition in the mosaic intercropping than in row-wise intercropping. This effect is modulated by the variety of intercropped PP (different varieties of PP were grown at the Bangalore and Kolli Hills research sites), soil quality and local weather. There are several factors, such as light, soil moisture and nutrient, that affect the yield of each component crop in intercropping (Bedoussac et al., 2015). The difference in penetration of light into canopy is considered to be a key factor affecting photosynthesis and ultimately growth and yield (Gwathmey and Clement, 2010;Kaggwa-Asiimwe et al., 2013). In our study, the reduction in light availability to relatively short FM plants standing next to taller PP plants in the mosaic intercropping treatments T5 & T6 (see supplementary data) could be a factor impacting growth, since in all row-wise intercropping designs PP and FM rows were well spaced to avoid a shading effect, which is not the case in the mosaic design. Similar results were reported by Martin and Snaydon (1982) and Dubey et al. (1995), who reported higher yield for barley-beans and sorghum-soybean in row-wise intercropping than mosaic (mixed within rows), respectively.
The intercropping designs tested in this study illustrate that the row-wise intercropping treatment T3+ (2:8 with biofertilizer) performed consistently better than the other arrangements, most probably due to the release of intraspecific competition. Effects of the spatial arrangement of component plants in intercropping have been shown to be species specific. Chen et al. (2004), Lauk and Lauk (2008) and Aynehband et al. (2010) have shown mixing of component plants within rows (mosaic pattern) to be the best arrangement for barley-peas, maize-soybean, and maizeamaranth, respectively. In contrast, Martin and Snaydon (1982) and Dubey et al. (1995) reported higher yields for barley-beans and sorghum-soybean sown in alternate rows than mixed within rows, respectively. Inter-specific competition could occur when two species are planted together and such competition could lead to a decrease in plant growth and yield (Jensen, 1996). In a cereal-legume intercropping system there is a significant number of days of overlapping growth period, and inter-specific competition between component crops could lead to a decrease in yield (Clément et al., 1992;Oljaca et al., 2000;Karasawa and Takebe, 2012). Therefore, the spatial arrangement between the plants needs to be carefully optimized. In this study, PP had a head start of 45 days (polybag transplantation) compared to FM, which provided PP a competitive advantage to acquire more resources (light, nutrients, and water) through its wellestablished root network, and FM may face, additionally, shading effects due to tall PP plants.

Effects of Biofertilizers
In the 2017/18 field trial, at both experimental sites, the effect of biofertilizer application was positive and showed an increase in total yield (Figure 7). The positive effect of biofertilization did not differ among intercropping treatments with different spatial arrangements (Supplementary Tables 7, 8). The effect of biofertilization was, however, specific to each component plants in the PP-FM intercropping system. Total biomass and straw yield per plant in FM was not significantly affected by biofertilization, but grain yield was significantly increased (Supplementary Table 5) similar to the observations made by Mathimaran et al. (2020). In the case of PP, the effect of biofertilization was significant on total biomass, straw, and grain yield per plant. The results of this study are in agreement with findings of Mäder et al. (2011) who reported that combined application of AMF and PGPR improves grain yield. Previous studies (Reddy, 2012;Mathimaran et al., 2017) have shown that better phosphorus uptake and crop tolerance to biotic and abiotic stresses via PGPR are among the most common mechanisms through which biofertilizers improve crop growth. The increase in grain yield in both component plants (FM and PP) in intercropping was the result of an increased number of panicle and grain weight per panicle in FM and number of pod and pod weight per plant in PP (see supplementary data). Since the process of pod and panicle formation is influenced by light availability, nutrients and soil moisture (Härdter and Horst, 1991;Thorsted et al., 2006;Raza et al., 2020), the yield improvement in row-wise intercropping could be attributed to efficient utilization of nutrients through the applied biofertilizers.

Effect of Intercropping and Biofertilizers on Water Relations of FM
In this study, the water relations (predawn LWP) of FM decreased significantly in mosaic treatments as compared to row-wise and monoculture treatments (Figures 8A,B). The trend in predawn LWP (Figures 8A,B) can also be compared with the trend in biomass production per plant (Figures 6A, 7A), therefore, competition for water could be the limiting factor here which influenced the yield and effectiveness of intercropping treatments at the Bengaluru site. Our results suggest that there exists an important degree of below-ground competition for water between PP and FM, and the facilitative effect of bioirrigation is suppressed. Similar results have been reported by Ludwig et al. (2004), They found that HL performing trees extracted a significant amount of water from the topsoil layer that resulted in lower LWP in understorey grasses. However, grasses were able to absorb soil moisture released by tree due to HL.
One of our objectives was to find out if CMN can facilitate the transfer of bioirrigated water from PP to FM and improve the water relations of FM in intercropping treatments. The results from the 2017/18 field trial showed that CMN did not affect the water relations (predawn LWP) of FM in intercropping treatments. However, at week 1 and 2 (first and second week of November 2017) FM in T3+ had higher (less negative), but not significant, LWP than T3-. Similarly, FM in monoculture treatment showed a higher LWP with CMN than without CMN ( Figure 8B). Since, we observed similar effects of CMN in both monoculture and 2:8 row-wise intercropping, we cannot assign this to bioirrigation. The effect of CMN changed over time, and at week 3 (third week of November 2017) treatments T1+, T3+, and T5+ (with CMN) had a lower LWP than T1-, T3-, and T5-(without CMN). The effect of different treatments, biofertilization and times (weekly measurement) had significant interaction with each other (Supplementary Table 9).
In this study, we could not find out if the positive intercropping effect by CMN was due to bioirrigation. This is in contrast to two greenhouse studies that we performed earlier (Singh et al., 2019(Singh et al., , 2020 where such an effect of bioirrigation was detectable. The discrepancy among these studies illustrated the difficultly to upscale plant-plant interactions obtained in controlled green-house studies to the field. The average hyphal spread rate of Glomus species is 0.7-0.8 mm per day (Jakobsen et al., 1992;Harinikumar and Bagyaraj, 1995), We did not check for the spread of CMN between PP and FM, but it is possible that the AMF introduced with the biofertilizer could not cover the distance of 45 cm between PP and FM in intercropping treatments and thus, a potential facilitative effect of bioirrigation through CMN was not observed in the field.

CONCLUSIONS
This study shows that intercropping has a positive effect on total yield of PP and FM, but that this effect varies across the sites, likely based on site characteristics such as soil type and weather. Importantly, our experiment demonstrates for the first time that the positive effect of intercropping on yield can be enhanced by the application of biofertilizers. In conclusion, the answers to our three research questions are therefore as follows: (i) the spatial arrangement of intercropping partners does affect the straw and grain yield in a FM-PP intercropping system, and the optimal spatial arrangement for PP-FM intercropping system depends on geographic location (local weather conditions) and plant variety. In general, the row-wise treatment (T3+) resulted in better yields than the mosaic treatments at the Bengaluru site, while at the Kolli Hills site in 2017/18, both the row-wise treatment (T3+) and the mosaic treatment (T5+) performed equally well. In addition, (ii) we show that the application of biofertilizer promotes yield in intercropping systems, and the spatial arrangement of component plants does not affect the effect of biofertilization. The effect of biofertilization is mainly due to the promotion of PP. We further show that (iii) the spatial arrangement of plants is a key factor that affects the competition for topsoil moisture between PP and FM.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.