Abstract
Introduction:
Water scarcity is a major constraint limiting agricultural productivity in semi-arid regions. Developing sustainable production systems that improve water-use efficiency while maintaining soil physicochemical quality is therefore essential. This study evaluated the combined effects of cropping systems and deficit irrigation on instantaneous water-use efficiency (iWUE), soil physicochemical properties, nutrient dynamics, crop growth, and productivity in dry bean (Phaseolus vulgaris L.)-based intercropping systems with Cucumis myriocarpus and Cleome gynandra.
Methods:
A field experiment was conducted over two growing seasons (2023/24 and 2024/25) using three cropping systems and three irrigation regimes supplying 100% (full irrigation, FI), 75% (moderate deficit irrigation, DI), and 50% (severe DI) of crop water requirements. Measurements included iWUE, soil physicochemical properties, exchangeable cations, crop growth, biomass production, and yield.
Results:
Irrigation level had a greater influence on iWUE than cropping system. Severe deficit irrigation (50% DI) significantly reduced crop growth, biomass accumulation, and yield, whereas full irrigation (100% FI) generally produced the highest performance. Moderate deficit irrigation (75% DI) maintained satisfactory productivity while reducing water input, demonstrating its potential as a water-saving strategy. Both cropping system and irrigation regime significantly affected soil physicochemical properties and exchangeable cations (Ca, Mg, and K). Intercropping improved soil conditions by reducing bulk density and enhancing soil pH and nutrient cycling compared with sole cropping. However, dry bean-based systems maintained higher exchangeable cation concentrations and cation exchange capacity than C. gynandra-based systems. Soil nutrient availability declined progressively with increasing water deficit. Sole dry bean under full irrigation recorded the highest growth and yield. Although intercropping reduced dry bean yield compared with sole cropping, the dry bean + C. myriocarpus intercrop consistently outperformed the dry bean + C. gynandra intercrop under water-limited conditions.
Discussion:
The findings demonstrate that irrigation management plays a more decisive role than cropping system in determining water-use efficiency and crop productivity under semi-arid conditions. Moderate deficit irrigation (75% DI), when combined with compatible intercropping systems, particularly dry bean + C. myriocarpus, offers a practical strategy for conserving irrigation water while sustaining crop productivity and improving soil physicochemical quality. These results support the adoption of integrated water and cropping management practices to enhance the resilience and sustainability of crop production in water-limited environments.
1 Introduction
Agricultural production systems worldwide are increasingly under pressure to meet rising food demand while simultaneously conserving natural resources. In response to growing food security needs, the intensive use of natural resources, particularly fertile soils, have intensified across production systems (). This has led to widespread soil degradation, including accelerated erosion, nutrient depletion, and the deterioration of key soil physicochemical properties (). These challenges are especially severe in arid and semi-arid regions, where fragile ecosystems are highly vulnerable to inappropriate land management practices. Consequently, there is an urgent need to adopt sustainable agricultural strategies, such as intercropping, to restore soil physicochemical properties and improve water-use efficiency (; Tefera et al., 2025).
Climate-smart agriculture promotes diversified cropping systems that integrate complementary crop species to maximise productivity while conserving natural resources and sustaining ecosystem services (). Within these systems, intercropping, particularly the combination of legumes and cereals, has received considerable attention due to its potential to improve soil physicochemical properties, enhance nutrient cycling, increase water-use efficiency, and boost overall system productivity (; ). Moreover, intercropping is particularly beneficial for smallholder farmers, as it offers both agronomic and nutritional advantages while reducing vulnerability to environmental and economic stresses (Toker et al., 2024).
In addition to declining soil quality, water scarcity has emerged as one of the most significant constraints to agricultural production. Water deficit stress impairs plant growth by restricting leaf expansion, reducing photosynthesis, water-use efficiency and ultimately lowering crop yields (; ). These challenges are expected to intensify under climate change, which is projected to increase temperatures, alter precipitation patterns, and heighten the frequency and severity of drought events (). In sub-Saharan Africa, water scarcity is already a major constraint to agricultural production, particularly for subsistence and smallholder farmers who rely heavily on rain-fed farming systems (). Existing limitations in water availability and access continue to constrain crop productivity, and climate change projections indicate that these pressures are likely to intensify in the future (Touch et al., 2024). These challenges underscore the urgent need to develop and implement efficient water management strategies, particularly in water-scarce regions experiencing rapid population growth. As water resources and arable land continue to decline, evaluating irrigation strategies that can sustain agricultural productivity while improving water use efficiency under conditions of limited water availability has become increasingly important. Deficit irrigation scheduling has emerged as a promising approach for enhancing water use efficiency by optimizing the use of limited water supplies while maintaining acceptable crop yields. However, despite extensive research demonstrating its potential, the adoption and practical implementation of deficit irrigation by smallholder farmers remain limited, highlighting the need for further research to improve its feasibility, accessibility, and field-level application.
Soil physicochemical properties are fundamental indicators of soil quality and play a central role in ecosystem functioning and agricultural productivity. Inadequate management of these properties exacerbates soil degradation by creating unfavorable conditions in arable lands, leading to the depletion of soil organic matter and essential nutrients (; ). Consequently, declining soil quality results in significant reductions in crop productivity (). Healthy soils support plant, animal, and human life by linking agricultural productivity with effective soil management and sustainable food production systems (). Soil physical properties form the basis for chemical and biological processes by regulating the movement of air, water, and dissolved nutrients within the soil profile, thereby influencing seed germination, root development, and overall plant growth (). For instance, increased soil bulk density can reduce pore space and water infiltration due to compaction, ultimately limiting nutrient cycling and plant productivity (). Soil chemical properties, including the concentration and availability of cations and anions, are equally important for nutrient management and fertilisation planning. Among these, soil pH and electrical conductivity (EC) are particularly critical, as they strongly influence nutrient availability and serve as key indicators of soil degradation ().
Synthetic fertilisers have played a major role in increasing crop yields since the twentieth century and have contributed substantially to global food production (Tang et al., 2022). However, their excessive and prolonged use can lead to nutrient imbalances, soil acidification, reduced microbial diversity, and increased greenhouse gas emissions (Tripathi et al., 2020). These challenges are particularly severe for subsistence and smallholder farmers in sub-Saharan Africa, many of whom lack the financial capacity to purchase sufficient fertilisers to replenish nutrients removed through crop harvesting (; Tully et al., 2015). Consequently, there is a growing need to adopt sustainable agricultural technologies, including intercropping and integrated soil physicochemical management, to maintain and enhance crop productivity in smallholder farming systems (Verma et al., 2019).
Cropping systems refer to the spatial and temporal arrangement of crops within a farming system and include practices such as monocropping, intercropping, crop rotation, and mixed cropping (Yang et al., 2020). These systems are essential for optimising resource use, improving soil physicochemical properties, and enhancing resilience to environmental stresses, particularly in semi-arid regions where agricultural production is constrained by limited water availability and high climatic variability (). In such environments, well-designed cropping systems can improve water-use efficiency (WUE), reduce production risks, and stabilise crop yields (Yang H. et al., 2021).
Within legume-based intercropping systems, nutrient acquisition is enhanced through the ability of legumes to fix atmospheric nitrogen via biological nitrogen fixation. This process facilitates interspecific nitrogen transfer, improves nutrient availability, and increases overall nutrient-use efficiency within the system (; ) Legumes therefore play a key role in sustainable agriculture by improving soil fertility and supporting long-term soil productivity (). Intercropping legumes with cereal crops, such as maize, have been widely recognized for strengthening food production systems in developing countries, particularly under water-limited conditions. In addition to increasing productivity, these systems enhance soil health and nutrient cycling through biological nitrogen fixation ().
South Africa has a rich diversity of indigenous plant species with significant nutritional, medicinal, and economic value, making their cultivation increasingly important for sustainable agricultural development (). Many indigenous and medicinal plants have been shown to produce acceptable yields and generate substantial economic returns under water-scarce conditions (). Among these, African leafy vegetables have gained increasing attention due to their contributions to food and nutrition security, dietary diversity, and rural livelihoods (). Their importance continues to grow because of their nutritional benefits and strong adaptability to marginal environments, including water-deficit conditions. Many of these species can maintain productivity and economic value under limited water availability, making them suitable for climate-resilient agricultural systems ().
Among the most important indigenous leafy vegetables is Cleome gynandra, widely cultivated and consumed across Africa and parts of Asia (). It is an important source of vitamins, minerals, and other essential nutrients, particularly for resource-constrained communities (). Another notable species is Cucumis myriocarpus, native to several southern African countries, including South Africa, Zambia, Botswana, Mozambique, and Lesotho (). Species within the Cucumis genus are generally drought tolerant and capable of surviving under minimal water supply (). Their physiological and morphological adaptations enable them to withstand environmental stress, remain dormant during adverse conditions, and exhibit moderate resistance to diseases such as downy mildew, scab, and mosaic virus (; ). Furthermore, these crops perform well in arid environments and can thrive in soil with low water-holding capacity ().
Despite increasing research on intercropping systems, limited information is available on the integration of dry beans with African leafy vegetables, particularly under deficit irrigation conditions. The drought tolerance and nutritional value of African leafy vegetables, combined with the soil fertility benefits and nitrogen-fixing ability of legumes such as dry beans, suggest strong potential for developing resilient and sustainable cropping systems. Such systems could enhance resource-use efficiency, improve soil physicochemical properties, increase water-use efficiency, and strengthen food and nutritional security in resource-constrained environments. However, the combined effects of deficit irrigation and dry bean–African leafy vegetable intercropping on soil physicochemical properties, water-use efficiency, and crop productivity remain poorly understood. Addressing this gap could provide valuable insights into the development of climate-smart and sustainable production systems that support smallholder farmers in water-limited regions.
2 Materials and methods
2.1 Study area
The study was conducted at the University of Limpopo Experimental Farm (UL Farm), South Africa (23° 50’ 36.86’’ S, 29° 40’ 54.99’’ E). The planting season was from November 2023 to March 2024 (season 1) and from November 2024 to March 2025 (season 2). The UL Farm is situated in a semi-arid region within the Polokwane Municipality, Limpopo Province, South Africa. The area experiences winter temperatures ranging from 16 °C to 18 °C (minimum) and 20 °C to 30 °C (maximum), while summer temperatures range from 18 °C to 22 °C (minimum) and 28 °C to 38 °C (maximum). Relative humidity ranges from 30–40% (minimum) to 85–95% (maximum) ().
The soils at the study site are classified as calcareous with an effective rooting depth of approximately 40 cm, according to the South African Soil Classification System (). The field had previously been cultivated under a maize–legume intercropping system. The orthic topsoil is about 70.5 cm deep and is characterised by a dark reddish-brown color (5YR 3/4) when dry and a dark brown color (7.5YR 3/3) when wet. It has a clay content ranging from 15% to 35%, classifying it as sandy clay loam. The subsoil is classified as a red apedal B horizon with a depth of 33.5 cm. It appears brownish yellow (10YR 6/6) when dry and yellowish red (5YR 4/6) when wet, with a clay content of 0–5%.
2.2 Study design
Data were analyzed using a linear mixed-effects model in which irrigation level, cropping system, growth stage, and their interactions were treated as fixed effects, while block and subplot effects were treated as random effects to account for the split-plot experimental design. This approach appropriately accommodated the hierarchical structure of the experiment and enabled accurate estimation of treatment effects and standard errors. In the split-plot design, irrigation levels were assigned to main plots, while cropping systems were assigned to subplots. The model was specified as follows:
where:
, observed response variable,
, overall mean,
, fixed effect of irrigation level,
, fixed effect of cropping system,
, irrigation × cropping system interaction,
, random effect of block,
, residual error.
For soil physicochemical analysis, a 5 × 3 split-plot design arranged in a Randomized Complete Block Design (RCBD) with three replications was used. The main-plot factor consisted of five cropping systems: (i) monocropped dry bean, (ii) monocropped Cucumis myriocarpus, (iii) monocropped Cleome gynandra, (iv) dry bean intercropped with C. myriocarpus, and (v) dry bean intercropped with C. gynandra. The subplot factor comprised three irrigation levels based on crop evapotranspiration (ETc): full irrigation (100% ETc), 75% deficit irrigation, and 50% deficit irrigation.
For dry bean performance, the experiment was analyzed as a 3 × 3 factorial arrangement in an RCBD with three replications. The first factor included three cropping systems: (i) monocropped dry bean, (ii) dry bean intercropped with C. myriocarpus, and (iii) dry bean intercropped with C. gynandra. The second factor consisted of the same three irrigation levels (100% ETc, 75% deficit irrigation, and 50% deficit irrigation).
For C. myriocarpus and C. gynandra, separate analyses were conducted using a 2 × 3 factorial arrangement in an RCBD with three replications. The first factor was cropping system, consisting of monocropping and intercropping with dry bean, while the second factor was irrigation level (100% ETc, 75% deficit irrigation, and 50% deficit irrigation). Each species was analyzed independently to account for differences in crop response.
2.3 Cultural practices
2.3.1 Soil preparation and characterization
The study area showed no topographical variation; therefore, simple random sampling was applied within a 56 × 54 m experimental field. Spacing between growth stage treatments was 11.5 m, irrigation treatments were separated by 4 m, and cropping treatments by 1.5 m (Supplementary Figure 1). Composite soil samples were randomly collected from each plot at a depth of 0–30 cm using an auger before planting to assess baseline soil fertility across four phenological stages. This sampling depth corresponds to the typical rooting zones of dry bean (Phaseolus vulgaris L.), Cucumis myriocarpus, and Cleome gynandra. The experimental field was prepared through ploughing, harrowing, and levelling prior to crop establishment (Supplementary Figure 2).
2.3.2 Irrigation
Soil was irrigated to 100% field capacity, and soil moisture was monitored using the gravimetric method under different regulated deficit irrigation regimes across the four phenological stages. Both gravimetric and volumetric water contents were used to determine soil moisture levels.
Gravimetric method: W, (Uthaiwat et al., 2024).
Volumetric water content: (Zamora-Ledezma et al., 2025).
For, Bulk density: Pw, ().
Regulated deficit irrigation was implemented using reference crop evapotranspiration (ET0) derived from high-resolution data obtained from the automated weather station (AWS) at the University of Limpopo (Figures 1, 2). These data form part of comprehensive South African climate datasets managed by the Agricultural Research Council. Irrigation scheduling was based on ET0 and crop-specific coefficients (Kc) to estimate crop evapotranspiration (ETc) and determine precise irrigation requirements.
Figure 1
Figure 2
The figures illustrate weekly meteorological conditions during the 2023/24 (Figure 1) and 2024/25 (Figure 2) growing seasons, including rainfall, temperature, relative humidity, and reference evapotranspiration (ET0). Both seasons were characterised by low and erratic rainfall, with only occasional rainfall peaks, highlighting the dependence of crop production on supplemental irrigation. Air temperatures fluctuated throughout both seasons, with slightly greater variability observed during 2024/25, while minimum temperatures remained relatively stable. Relative humidity was generally high; however, minimum relative humidity values exhibited greater fluctuations, indicating intermittent dry periods. Reference evapotranspiration (ET0) remained moderate across both seasons, suggesting relatively stable and manageable atmospheric water demand.
Daily ET0 was calculated using the Penman–Monteith method (). Crop water requirements (ETc) were subsequently estimated using the standard relationship ETc , ET0 × Kc, where ETc is crop evapotranspiration, ET0 is reference evapotranspiration, and Kc is the crop coefficient.
For the African leafy vegetables Cleome myriocarpa and Cleome gynandra, crop coefficient values developed for spinach and similar leafy vegetables were used to estimate ETc, with Kc values of 0.40 (initial stage), 1.52 (mid-season), and 0.35 (late season). Crop coefficient values for C. myriocarpa and C. gynandra are not yet well established in the literature. Therefore, the use of spinach-based coefficients was considered appropriate due to the physiological and morphological similarities between these species, including comparable growth habits, leaf morphology, canopy development, rooting characteristics, and their cultivation for vegetative biomass production (; ; ). sing these crop coefficients in conjunction with ET0 data obtained from the automatic weather station (AWS) (Figures 1, 2), crop water requirements were determined for each crop and growth stage.
Instantaneous Water Use Efficiency (iWUE) was measured during the four phenological stages as;
().
Where:
iWUE, Instantaneous Water Use Efficiency (μmol CO2 mmol-¹ H2O)
A, Net photosynthetic rate (μmol CO2 m-² s-¹)
E, Transpiration rate (mmol H2O m-² s-¹)
2.3.3 Fertilization
Fertilizer requirements were determined according to the recommendations of () and were further guided by soil analysis results (Table 1) and were further refined based on the soil analysis results (Table 1). The soil analysis indicated that phosphorus (P) and potassium (K) concentrations were within the optimal range; therefore, no additional P or K fertilizers were required to correct soil nutrient deficiencies.
Table 1
| Soil property | Value |
|---|---|
| Soil Texture | Loamy sand (82% sand, 10.4% silt, 7.2% clay) |
| Soil bulk density | 1.74 g cm-3 |
| Soil pH (H2O) | 7.86 |
| EC (EC) | 178 mS/m |
| Calcium (Ca) | 116 mg/kg |
| Potassium (K) | 203 mg/kg |
| Phosphorus (P) | 41 mg/kg |
| Nitrogen (N) | 2.3 mg/kg |
| Magnesium Mg | 147 mg/kg |
Physiochemical properties of the soil prior to planting.
Fertiliser applications were then tailored to the nutritional requirements of each crop. Cucumis myriocarpus received a basal application of 2:3:2 (30) fertiliser at 100 kg ha-¹, followed by a top-dressing of urea (46% N) at 60 kg ha-¹. Cleome gynandra received 2:3:2 (30) at 150 kg ha-¹ as a basal application, followed by a side-dressing of urea at 90 kg ha-¹. Phaseolus vulgaris L. (dry bean) received 2:3:2 (30) at 150 kg ha-¹ as a basal application and urea at 50 kg ha-¹ as a booster application. Potassium chloride (KCl; 50% K) was applied at 60 kg ha-¹ to support pod development and yield formation, following the crop-specific fertiliser.
2.4 Data collection
2.4.1 Soil parameters
Composite soil samples were randomly collected from the 0–30 cm depth. The samples were air-dried, gently crushed, and passed through a 2 mm sieve prior to analysis. Selected physical and chemical soil properties were then determined pre-planting and towards the end of the growth of the crops. Soil texture was analyzed using the hydrometer method (), while soil bulk density was calculated as the ratio of oven-dry soil mass to core volume (). Soil chemical properties, including pH, electrical conductivity (EC), phosphorus (P), potassium (K), calcium (Ca), magnesium (Mg), and sodium (Na), were also determined, as these parameters play a key role in the productivity of dry bean (Phaseolus vulgaris L.) and indigenous crops such as Cucumis myriocarpus and Cleome gynandra. Soil pH and EC were measured in a 1:2.5 soil-to-water suspension following (; ). Exchangeable cations (K, Ca, Mg, and Na) were extracted using the ammonium acetate method (M. Yang et al., 2024; ), while available phosphorus was determined using the Bray-1 extraction method (Venu et al., 2023). Total effective cation exchange capacity (CEC) was calculated as the sum of K + Ca + Mg + Na (cmol/kg) ().
2.4.2 Plant growth and yield parameters
Dry bean: Growth and yield components, including plant height, root diameter, stem diameter, stem length, and pod length, were measured using a measuring tape. Biomass, root weight, and pod weight were determined using a precision weighing scale.
African leafy crops: Growth and yield of Cucumis myriocarpus and Cleome gynandra were assessed by harvesting each plant as a bunch per plant. Leaves were cut above the growing point, following standard leaf harvest practices used for baby spinach and Swiss chard.
2.5 Data analysis
Data was subjected to analysis of variance (ANOVA) using using Statistix 10 (Analytical Software, Tallahassee, FL, USA, and means were separated using Tukey’s Honest significant difference. For all analyses, assumptions of normality and homogeneity of variance were evaluated using diagnostic residual plots. Where significant treatment effects were detected (P < 0.05), means were separated using Tukey’s Honest Significant Difference (HSD) test at the 5% significance level. Principal component analysis (PCA) was used to explore multivariate relationships among selected soil physicochemical properties and to assess treatment groupings under different irrigation levels and cropping systems. PCA was conducted using standardized variables, and the first two principal components were used to visualise associations among soil properties and treatment combinations.
3 Results
3.1 Soil water content across the four phenological growth stages
Figure 3 illustrates soil moisture dynamics under three irrigation regimes: 50% deficit irrigation (DI), 75% DI, and 100% full irrigation (FI), highlighting the effects of regulated deficit irrigation across four key phenological growth stages. Across all observation periods, soil water content was consistently highest under FI, intermediate under 75% DI, and lowest under 50% DI. Noticeable peaks around weeks 7 and 14 likely reflect recent irrigation or rainfall events, whereas subsequent declines corre-spond to progressive soil water depletion driven by crop uptake and drainage processes (Figure 3a).
Figure 3
Differences among treatments became more pronounced after flowering, particularly under 50% DI, indicating increasing water stress under reduced irrigation (Figure 3b). A general decline in soil moisture was observed across all treatments towards physiological maturity, mainly due to reduced irrigation input combined with sustained crop water demand (Figure 3d).
Overall, decreasing irrigation from 100% to 50% ETc consistently reduced soil water availability, with the strongest effects occurring after flowering and during the late growth stages. Although intermittent rainfall temporarily reduced differences among treatments, lower irrigation levels consistently resulted in reduced soil moisture throughout the crop cycle.
3.2 Instantaneous water use efficiency during four phenological growth stages
No significant interaction effects (p > 0.05) between cropping system and deficit irrigation were observed for instantaneous water use efficiency (WUE) across most growth stages in both seasons, except at the early pod stage in 2023/24, where sole dry bean under 75% DI showed a significant increase (Table 2).
Table 2
| Phenological growth stages | 50% flowering stage | Early-pod stage | Mid-pod stage | Maturity stage | ||||
|---|---|---|---|---|---|---|---|---|
| Cropping treatments | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 |
| Sole DB | 0,0748a | 0,0448a | 0,0703a | 0,1076a | 0,0408a | 0,0569a | 0,0181a | 0,0152a |
| DB + CM | 0,0483ab | 0,0387ab | 0,0371ab | 0,0406a | 0,03ab | 0,0494ab | 0,0178a | 0,0122a |
| Sole CM | 0,0271bc | 0,0257b | 0,0301ab | 0,0382a | 0,0255ab | 0,0287b | 0,0168a | 0,0104a |
| DB + CG | 0,0214bc | 0,0253b | 0,0281ab | 0,0212a | 0,0184b | 0,0267b | 0,0153a | 0,0104a |
| Sole CG | 0,02c | 0,0231b | 0,0143b | 0,0203a | 0,0122b | 0,0258b | 0,0146a | 0,0082a |
| HSD (Critical value for comparison) | 0,0272 | 0,0272 | 0,0361 | 0,0361 | 0,0501 | 0,0501 | 0,0238 | 0,0238 |
| Response | * | * | * | ns | * | * | ns | ns |
| Irrigation levels | ||||||||
| 100%FI | 0,0361b | 0,0215a | 0,0457a | 0,0904a | 0,0308a | 0,0267a | 0,0127a | 0,0127a |
| 75%DI | 0,0546a | 0,0206a | 0,0326a | 0,0279a | 0,0256ab | 0,0258a | 0,0218a | 0,0094b |
| 50%DI | 0,0242b | 0,0158a | 0,0297a | 0,0184a | 0,0197b | 0,0226a | 0,02a | 0,0084b |
| HSD (Critical value for comparison) | 0,0179 | 0,0188 | 0,0549 | 0,0549 | 0,0388 | 0,0388 | 0,0124 | 0,0124 |
| Response | * | ns | ns | ns | * | ns | ns | * |
Effect of cropping and irrigation on water use efficiency in different phenological growth stages.
DB, Dry beans; CM, Cucumis myriocarpus; CG, Cleome gynandra; FI, Full irrigation; DI, Deficit irrigation; EC, Electrical conductivity. Column means followed by the same letter were not significantly different (p < 0.05) according to Tukey’s HSD (Honestly Significant Difference) Test. Response **(Highly significant p ≤ 0.01), *(Significant p ≤ 0.05), ns (not significant). Bold values are mean values for honest significant difference.
Cropping system significantly influenced WUE at the 50% flowering, early pod, and mid pod stages in 2023/24, and at the 50% flowering and mid pod stages in 2024/25. At 50% flowering, sole dry bean (DB) consistently recorded the highest WUE in both seasons (0.0748 and 0.0448), significantly exceeding sole Cleome gynandra (CG) and most intercrop treatments (Table 2). A similar trend was observed at the early pod stage in 2023/24, where sole DB recorded the highest value (0.0703), while sole CG showed the lowest (0.0143).
At the mid pod stage, sole DB consistently outperformed all other treatments in both seasons, with values of 0.0408 and 0.0569 in 2023/24 and 2024/25, respectively. The lowest values were generally recorded under sole CG and the DB + CG intercrop (Table 2). However, at physiological maturity, differences among cropping systems were no longer significant in either season, indicating convergence of treatment effects as crops reached maturity.
Irrigation level also affected WUE at specific growth stages. During 50% flowering in 2023/24, 75% DI produced the highest WUE (0.0546), significantly higher than both FI and 50% DI. At the mid pod stage in the same season, FI recorded the highest value (0.0308), while 50% DI recorded the lowest (0.0197). At maturity in 2024/25, FI again resulted in significantly higher WUE (0.0127) compared with both deficit irrigation treatments (Table 2).
In contrast, irrigation effects were not significant at most growth stages in 2024/25, suggesting a weaker response of WUE to irrigation level under those seasonal conditions. Overall, sole dry bean consistently exhibited higher WUE than intercrop and sole Cleome gynandra treatments, particularly during vegetative and reproductive stages. While irrigation effects were variable across seasons, moderate deficit irrigation (75% DI) occasionally improved WUE compared with full irrigation.
3.3 Effect of cropping treatments and irrigation levels on selected soil physical and chemical properties
3.3.1 Main effects of cropping systems and irrigation levels on selected soil physicochemical properties
Cropping systems significantly influenced most soil chemical properties in both seasons, although phosphorus (P) concentrations did not differ significantly among treatments. During the 2024/25 season, sole dry bean (DB) and the dry bean + Cucumis myriocarpus (DB + CM) intercrop maintained relatively higher P concentrations, whereas the dry bean + Cleome gynandra (DB + CG) intercrop and sole C. gynandra recorded lower values (Table 3). The relatively stable P concentrations across cropping systems may be attributed to the low mobility of phosphorus in soil, where diffusion rather than mass flow is the dominant transport mechanism. Consequently, changes in crop composition may have had a limited effect on overall soil P availability over the relatively short duration of the experiment. Similar findings have been reported in legume-based intercropping systems, where changes in residual soil P are often less pronounced than changes in more mobile nutrients such as K and Na ().
Table 3
| Parameter | Bulkdensity (g/cm³) | EC (µS/cm) | pH H2O | pH KCL | ||||
|---|---|---|---|---|---|---|---|---|
| Seasons | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 |
| Cropping systems | ||||||||
| Sole DB | 1,347a | 1,334a | 151,28abc | 139,79a | 7,794a | 7,165bc | 7,43a | 6,52a |
| DB + CM | 1,312a | 1,311a | 153,23a | 140,16a | 7,812a | 7,221a | 7,438a | 6,53a |
| Sole CM | 1,323a | 1,306a | 149,3bc | 138,01ab | 7,73bc | 7,155bc | 7,379bc | 6,466b |
| DB + CG | 1,279a | 1,323a | 1,29a | 140,1a | 7,749b | 7,187ab | 7,418ab | 6,521a |
| Sole CG | 1,368a | 1,323a | 148,19c | 137,28b | 7,705c | 7,125c | 7,365c | 6,458b |
| HSD (Critical value for comparison) | 0,107 | 0,053 | 3,602 | 2,365 | 0,035 | 0,05 | 0,043 | 0,023 |
| Response | ns | ns | ** | * | * | * | * | ns |
| Irrigation levels | ||||||||
| 100% FI | 1,256b | 1,331a | 145,42b | 141,28a | 7,633c | 7,206b | 7,337b | 6,468b |
| 75% DI | 1,359a | 1,329a | 153,21a | 134,22b | 7,772b | 7,049c | 7,343b | 6,452c |
| 50% DI | 1,363a | 1,278b | 153,91a | 141,7a | 7,875a | 7,257a | 7,542a | 6,578a |
| HSD (Critical value for comparison) | 0,065 | 0,032 | 2,691 | 1,692 | 0,022 | 0,027 | 0,027 | 0,014 |
| Response | ns | ns | * | * | * | * | * | * |
Effect of cropping treatments and irrigation levels on selected soil physicochemical properties during 2023/24 season (Experiment 1) and 2024/25 season (experiment 2).
DB, Dry beans; CM, Cucumis myriocarpus; CG, Cleome gynandra, FI, Full irrigation; DI, Deficit irrigation; EC, Electrical conductivity. Column means followed by the same letter were not significantly different (p < 0.05) according to Tukey’s HSD (Honestly Significant Difference) Test. Response **(Highly significant p ≤ 0.01), *(Significant p ≤ 0.05), ns (not significant). Bold values are mean values for honest significant difference.
In contrast, potassium (K), calcium (Ca), magnesium (Mg), and sodium (Na) were significantly affected by cropping system in both seasons. The DB + CG intercrop consistently recorded the highest K concentrations, while sole C. gynandra produced the lowest values across both seasons (Table 3). The elevated K concentrations under DB + CG may indicate enhanced nutrient mobilization within the rhizosphere, resulting from complementary root interactions and greater biological activity associated with intercropping (). Legume roots can alter rhizosphere chemistry through proton release, root exudation, and stimulation of microbial activity, thereby increasing nutrient solubilisation and availability (). Furthermore, differences in rooting depth and nutrient acquisition strategies between dry bean and C. gynandra may have reduced direct competition for some nutrients while improving overall nutrient recovery from the soil profile. Similar improvements in nutrient availability under legume-based intercrops have been reported by several authors, who attributed these responses to niche complementarity and enhanced nutrient cycling (; ; ).
Conversely, sole dry bean maintained the highest Ca concentrations and relatively high Mg levels, whereas DB + CM generally promoted greater Mg retention than the other cropping systems. These results suggest that nutrient uptake patterns differed among species and that interspecific competition may have altered nutrient extraction from the soil. The greater retention of Ca and Mg under sole dry bean may reflect lower overall nutrient demand compared with intercrop treatments, resulting in higher residual concentrations remaining in the soil.
Irrigation level significantly influenced K, Ca, Mg, and Na during the 2023/24 season, although these effects were less pronounced in 2024/25. Full irrigation (100% FI) generally resulted in the highest concentrations of K, Mg, and Na, whereas the lowest concentrations were recorded under 50% deficit irrigation (DI) (Table 4). The decline in exchangeable cations under severe water deficit is likely linked to reduced nutrient mobility and slower mineralisation processes under dry soil conditions (). Adequate soil moisture promotes microbial activity and nutrient cycling while facilitating nutrient transport to roots through both diffusion and mass flow (). Under deficit irrigation, reduced soil water content restricts these processes, limiting nutrient availability and altering nutrient distribution within the root zone.
Table 4
| Parameter | Bulk density (g/cm³) | EC (µS/cm) | pH H2O | pH KCL | ||||
|---|---|---|---|---|---|---|---|---|
| Seasons | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 |
| Sole DB x 100%FI | 1,312abc | 1,386a | 141,83fgh | 127,03h | 7,595h | 6,968a | 7,248hi | 6,45a |
| Sole DB x 75%DI | 1,305abc | 1,37ab | 140,04gh | 142,64cde | 7,509i | 7,152a | 7,243i | 6,41a |
| Sole DB x 50%DI | 1,3abc | 1,329abc | 137,59hi | 142,59cde | 7,871bc | 7,1a | 7,547bc | 6,57a |
| DB + CM x 100%FI | 1,326abc | 1,294abc | 131,05i | 131,19gh | 7,653gh | 7,19a | 7,319ghi | 6,35a |
| DB + CM x 75%DI | 1,325abc | 1,286abc | 168,3a | 129,33gh | 7,648gh | 7,12a | 7,249hi | 6,35a |
| DB + CM x 50%DI | 1,45ab | 1,386a | 164,72ab | 150,81a | 7,929b | 7,12a | 7,603ab | 6,63a |
| Sole CM x 100%FI | 1,288abc | 1,227c | 151,17de | 133,9fg | 7,795de | 7,12a | 7,387efg | 6,51a |
| Sole CM x 75%DI | 1,272abc | 1,325abc | 149,57ef | 132,13gh | 7,76ef | 7,06a | 7,372fg | 6,51a |
| Sole CM x 50%DI | 1,24bc | 1,313abc | 148,28efg | 146,3abc | 7,696fg | 7,12a | 7,367fg | 6,5a |
| DB + CG x 100%FI | 1,153c | 1,3abc | 142,42fgh | 138,89ef | 7,693fg | 7,02a | 7,34g | 6,48a |
| DB +CG x 75%DI | 1,44ab | 1,298abc | 160,96abc | 137,93ef | 7,677g | 7,19a | 7,34gh | 6,47a |
| DB + CG x 50%DI | 1,359abc | 1,375ab | 159,14bcd | 144,67bcd | 8,013a | 7,121a | 7,65a | 6,60a |
| Sole CG x 100%FI | 1,35abc | 1,302abc | 158,69bcd | 140,24de | 7,829cde | 7,09a | 7,47cde | 6,53a |
| Sole CG x 75%DI | 1,472a | 1,312abc | 154,48cde | 140,01de | 7,867bcd | 7,12a | 7,44def | 6,52a |
| Sole CG x 50%DI | 1,45ab | 1,306abc | 154,44cde | 148,36ab | 7,861bcd | 7,19a | 7,517bcd | 6,567a |
| Critical value for comparison (HSD) | 0,215 | 0,107 | 8,812 | 5,542 | 0,020 | 0,097 | 0,090 | 0,048 |
| Response | * | ** | ** | ** | ** | ns | ** | ns |
Interactive effect of cropping systems x irrigation levels on selected soil physicochemical properties during 2023/24 (experiment 1) and 2024/25 (experiment 2) seasons.
DB, Dry beans; CM, Cucumis myriocarpus; CG, Cleome gynandra; FI, Full irrigation; DI, Deficit irrigation; EC, Electrical conductivity. Column means followed by the same letter were not significantly different (p < 0.05) according to Tukey’s HSD (Honestly Significant Difference) Test. Response **(Highly significant p ≤ 0.01), *(Significant p ≤ 0.05), ns (not significant). Bold values are mean values for honest significant difference.
3.3.2 Interaction effects of cropping systems × irrigation levels on selected soil physicochemical properties
The interaction between cropping systems and irrigation levels significantly influenced bulk density and electrical conductivity (EC) in both seasons, indicating that the response to irrigation was strongly dependent on the cropping arrangement. The highest EC values were consistently observed under the dry bean + Cucumis myriocarpus intercrop (DB + CM) combined with deficit irrigation, particularly under DB + CM × 75% DI in 2023/24 and DB + CM × 50% DI in 2024/25. This suggests that the combined effect of intercropping and reduced irrigation enhanced the accumulation of soluble salts in the soil (Table 4).
These findings are consistent with (M. ), who reported that while elevated EC may constrain early plant growth, moderate increases at later stages can improve nutrient availability and uptake through enhanced osmotic adjustment. In contrast, the lowest EC values were generally recorded under fully irrigated sole dry bean treatments (Table 4). However, all EC values remained below 4 dS m-¹, indicating non-saline conditions according to USDA standards (). This confirms that although irrigation and cropping system influenced salt dynamics, salinity levels were not high enough to adversely affect crop productivity.
Bulk density also showed significant variation across treatment combinations. The lowest values were recorded under DB + Cleome gynandra × 100% FI, suggesting that intercropping combined with adequate soil moisture may enhance soil aggregation and reduce compaction (Table 4). In contrast, higher bulk density values observed during the early pod stage, particularly under 50% and 75% deficit irrigation, are likely associated with reduced soil moisture, which limits soil swelling and biological activity and increases compaction (). At the maturity stage, lower bulk density under full irrigation (100% FI) further suggests that adequate moisture promotes improved soil structure through enhanced aggregation and reduced compaction. Overall, consistently higher bulk density under 50% DI indicates that water stress may restrict root development and microbial activity, thereby reducing soil porosity, aeration, and structural stability.
Significant interaction effects were also observed for soil pH (H2O) and pH (KCl) during the 2023/24 season, although these effects were not significant in 2024/25 (Table 4). The highest pH values were recorded under DB + CG × 50% DI and DB + CM × 50% DI, indicating that intercropping under severe water deficit conditions promoted the accumulation of basic cations and increased soil alkalinity (Table 4). Conversely, sole dry bean treatments under full or moderate irrigation generally exhibited lower pH values.
The absence of significant interaction effects in the second season suggests that seasonal variability, environmental conditions, or soil buffering capacity may have moderated treatment responses over time. Overall, the results demonstrate that the influence of irrigation on soil physicochemical properties is strongly dependent on cropping system composition, with intercropping under deficit irrigation producing the most pronounced changes in soil chemical characteristics.
3.4 Effect of cropping treatments and irrigation levels on selected soil chemical properties
3.4.1 Main effects of cropping systems and irrigation levels on selected soil chemical properties
Cropping systems significantly influenced most soil cation properties in both seasons, although phosphorus (P) concentrations did not differ significantly among treatments. During the 2024/25 season, sole dry bean (DB) and the dry bean + Cucumis myriocarpus (DB + CM) intercrop maintained relatively higher P concentrations, whereas the dry bean + Cleome gynandra (DB + CG) intercrop and sole C. gynandra recorded lower values (Table 5). The absence of significant differences in P among cropping systems suggests that phosphorus availability was relatively stable across treatments, despite variations in crop composition and nutrient uptake patterns.
Table 5
| Parameter | P (mg/kg) | K (mg/kg) | Ca (mg/kg) | Mg (mg/kg) | Na (mg/kg) | |||||
|---|---|---|---|---|---|---|---|---|---|---|
| Seasons | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 |
| Cropping treatments | ||||||||||
| Sole DB | 31,508a | 36,84a | 147,92ab | 178,92ab | 108,36a | 129,97a | 137,19ab | 148,72a | 38,422ab | 41,094ab |
| DB + CM | 31,406a | 36,82a | 147,31ab | 168,39bc | 105,97ab | 122,42ab | 141,36a | 149,69a | 37,611ab | 39,301ab |
| Sole CM | 30,853a | 36,542a | 147,17ab | 153,75c | 105,06ab | 122,33ab | 136,69b | 145,78ab | 37,092b | 37,744b |
| DB + CG | 30,797a | 32,887b | 156,86a | 185,36a | 104,64b | 114,36b | 135,75b | 138,85bc | 39,603a | 46,014a |
| Sole CG | 29,919a | 31,708b | 144,81b | 184,75a | 103,72b | 112,08b | 135,69b | 134,31c | 38,567ab | 44,319ab |
| HSD comparison | 1,974 | 3,235 | 10,447 | 16,08 | 3,66 | 11,756 | 4,263 | 7,298 | 2,136 | 7,555 |
| Response | ns | ** | ** | ** | ** | ** | ** | ** | ** | ** |
| Irrigation levels | ||||||||||
| 100% FI | 32,113a | 35,45a | 153,82a | 178,55a | 107,93a | 120,65a | 141,77a | 147,28a | 39,75a | 44,595a |
| 75% DI | 31,65a | 34,88a | 147,73ab | 173,57a | 109,25a | 120,58a | 136,6b | 143,92a | 38,177b | 40,863a |
| 50% DI | 28,927b | 34,552a | 144,88b | 170,58a | 99,47b | 119,47a | 133,65c | 139,22b | 36,85b | 39,625a |
| HSD comparison | 1,373 | 1,924 | 6,477 | 11,021 | 2,427 | 7,77 | 2,7 | 4,58 | 1,383 | 5,519 |
| Response | ns | ns | ** | ns | ** | ns | ** | * | ** | ns |
Effect of cropping treatments and irrigation levels on selected soil chemical properties during 2023/24 season (experiment 1) and 2024/25 season (experiment 2).
DB, Dry beans; CM, Cucumis myriocarpus; CG, Cleome gynandra; FI, Full irrigation; DI; Deficit irrigation; P, Phosphorus; K, Potassium; Ca, Calcium; Mg, Magnesium; Na, Sodium. Column means followed by the same letter were not significantly different (p < 0.05) according to Tukey’s HSD (Honestly Significant Difference) Test. Response **(Highly significant p ≤ 0.01), *(Significant p ≤ 0.05), ns (not significant p ≥ 0.05). Bold values are mean values for honest significant difference.
In contrast, potassium (K), calcium (Ca), magnesium (Mg), and sodium (Na) were significantly affected by cropping system in both seasons. The DB + CG intercrop consistently recorded the highest K concentrations, while sole C. gynandra produced the lowest K levels across both seasons (Table 5). Conversely, sole dry bean maintained the highest Ca concentrations and relatively high Mg levels, whereas DB + CM generally promoted greater Mg retention than the other cropping systems. These findings indicate that crop species composition influenced nutrient accumulation and cycling through differences in root architecture, nutrient uptake efficiency, residue quality, and rhizosphere interactions. The elevated K concentrations under DB + CG suggest enhanced nutrient mobilization and recycling within the intercrop system, while the greater Ca and Mg concentrations under sole dry bean indicate improved retention of these base cations in the soil.
Irrigation level also significantly influenced K, Ca, Mg, and Na during the 2023/24 season, although its effects were less pronounced in 2024/25. Full irrigation (100% FI) generally resulted in the highest concentrations of K, Mg, and Na, whereas the lowest concentrations were recorded under 50% deficit irrigation (DI) (Table 5). A similar trend was observed for Ca, with concentrations declining as irrigation water decreased and 50% DI producing significantly lower values than the fully irrigated treatment during the 2023/24 season.
The decline in exchangeable cation concentrations under severe water deficit may be attributed to reduced nutrient mineralization, restricted nutrient transport to plant roots through mass flow, and increased physiological stress associated with limited soil moisture. Reduced soil water availability can also limit nutrient diffusion and root activity, thereby decreasing nutrient uptake efficiency and altering nutrient distribution within the soil profile. However, the lack of significant irrigation effects on most cations during the 2024/25 season suggests that seasonal environmental conditions, residual nutrient accumulation, or soil buffering processes may have moderated treatment responses.
3.4.2 Interaction effects of cropping systems × irrigation levels on soil chemical properties
The interaction between cropping systems and irrigation levels had a highly significant influence on phosphorus (P), potassium (K), calcium (Ca), magnesium (Mg), sodium (Na), and total exchangeable cation exchange capacity (TEC–CEC) in most cases, indicating that soil nutrient dynamics were strongly governed by the combined effects of crop arrangement and water availability (Tables 4, 6).
Table 6
| Parameter | P (mg/kg) | K (mg/kg) | Ca (mg/kg) | Mg (mg/kg) | Na (mg/kg) | |||||
|---|---|---|---|---|---|---|---|---|---|---|
| Seasons | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 |
| Sole DB x 100%FI | 34,642a | 44,242a | 165,17a | 202,67a | 117a | 133,92a | 155,5a | 162a | 41,667a | 58,875a |
| Sole DBx 75% DI | 34,517a | 41,217ab | 164ab | 199a | 116,33a | 131,75a | 142,58b | 157,92ab | 41,5a | 47,583ab |
| Sole DB x 50%DI | 34,392a | 39,033abc | 157,08abc | 186,5ab | 113,17ab | 126,42a | 142,25bc | 156,08abc | 41,167ab | 40,942b |
| DB + CM x 100%FI | 34a | 37,917abcd | 151,17abc | 182,58ab | 111,92abc | 126,33a | 142,25bc | 153abc | 40,692ab | 45,217ab |
| DB + CM x 75%DI | 33,375ab | 37,683abcd | 150abc | 182,17ab | 110,42abcd | 124,75a | 141,75bc | 144,17bcd | 40,558abc | 44,433ab |
| DB + CM x 50%DI | 32,167abc | 36,933bcde | 148,58abc | 181,67ab | 107,92bcde | 124,25a | 141,58bc | 144bcd | 40,267abc | 40,758b |
| Sole CM x 100%FI | 32,083abc | 34,708bcdef | 148,25abc | 180,08ab | 105,83bcde | 123,67a | 138,17bcd | 144bcd | 40abc | 40,358b |
| Sole CM x 75%DI | 31,758abcd | 34,683bcdef | 145,92abc | 173,08ab | 105cde | 117,97a | 137,58bcd | 143,08bcde | 38,5abcd | 40,317b |
| Sole CM x 50%DI | 31,417abcd | 33,658cdefg | 144,92abc | 171,75ab | 104,58cdef | 117,25a | 135,83bcde | 142,83bcde | 38abcd | 39,042b |
| DB + CG x 100%FI | 29,217bcd | 33,425cdefg | 144,83abc | 171,42ab | 103,5defg | 116,03a | 135,25bcde | 141,5cde | 37,917abcde | 39,025b |
| DB +CG x 75%DI | 28,442cd | 31,617defg | 143,42abc | 170,17ab | 101,08efgh | 115,83a | 133,25cdef | 135,17de | 36,583bcdef | 38,85b |
| DB + CG x 50%DI | 28,05cde | 31,583defg | 143,25bc | 155,33b | 100,17efgh | 115,25a | 131,75def | 135,08de | 36,033cdef | 38,55b |
| Sole CG x 100%FI | 27,942cde | 30,842efg | 143,25bc | 153,42b | 96,83fgh | 110,85a | 129,58def | 132,67de | 35def | 37,667b |
| Sole CG x 75%DI | 27,533de | 29,575fg | 141,42c | 152,5b | 95,75gh | 109,88a | 127ef | 132,17de | 33,333ef | 37,017b |
| Sole CG x 50%DI | 23,917e | 27,319g | 140,92c | 151,17b | 93,75h | 109,36a | 125,75f | 128,39e | 32,667f | 36,786b |
| Critical value for comparison (HSD) | 4,4975 | 6,3033 | 21,21 | 36,09 | 7,948 | 25,444 | 8,8415 | 15,012 | 4,5303 | 18,076 |
| Response | ** | ** | ** | ** | ** | ns | ** | ** | ** | ** |
Interactive effect of cropping systems x irrigation levels on selected soil chemical properties during 2023/24 (experiment 1) and 2024/25 (experiment 2) seasons.
DB, Dry bean; CM, Cucumis myriocarpus; CG, Cleome gynandra; FI, Full irrigation; DI, Deficit irrigation; P, Phosphorus; K, Potassium; Ca, calcium; Mg, magnesium; Na, sodium. Column means followed by the same letter were not significantly different (p<0.05) according to Tukey’s HSD (High Significant Difference) Test. Response **(Highly significant P ≤ 0.01), *(Significant P ≤ 0.05), ns (not significant P ≥ 0.05). Bold values are mean values for honest significant difference.
Across both seasons, the highest concentrations of P, K, Ca, Mg, and Na were consistently observed under sole dry bean (Phaseolus vulgaris L.) combined with full irrigation (100% FI), whereas the lowest values were generally recorded under sole Cleome gynandra subjected to 50% deficit irrigation (Table 4). For example, sole DB × 100% FI produced the highest K concentrations (165.17 and 202.67 mg kg-¹) and Ca concentrations (117.00 and 133.92 mg kg-¹) across the two seasons (Table 6). These results suggest that a legume-based system combined with adequate water supply enhances nutrient availability and retention, likely through improved root activity, nutrient cycling, and soil biological processes.
A general decline in nutrient concentrations with increasing irrigation deficit was observed across most cropping systems, particularly for K, Mg, and Na. Concentrations progressively decreased from full irrigation to 50% DI in both sole and intercropped treatments (Table 6). Although DB + CG maintained relatively higher K levels than most treatments, severe water deficit still resulted in substantial nutrient depletion. Similarly, sole CG consistently recorded the lowest concentrations of most cations across irrigation regimes, suggesting lower nutrient recycling efficiency and/or greater nutrient depletion (Table 6).
Higher P levels under moderate deficit irrigation (75% DI) compared with full irrigation indicate that water stress may reduce plant uptake, leading to nutrient accumulation in the soil. This contrasts with (), who reported increased P availability under water deficit in quinoa. The discrepancy may be attributed to reduced root growth and limited P mobility under dry conditions, since phosphorus transport is strongly dependent on soil moisture and diffusion processes.
Intercropping also influenced nutrient dynamics through increased interspecific competition. Lower P levels in intercropping systems likely resulted from overlapping root zones, which enhanced nutrient uptake and reduced residual soil () This agrees with (Xu et al., 2021), who reported higher soil P under monocropping compared with intercropping due to reduced competition.
Potassium dynamics were similarly affected by the interaction between cropping system and irrigation. Higher K levels under sole cropping and under 100% and 75% FI suggest that adequate soil moisture supports K availability and reduces plant stress, thereby maintaining higher residual soil K (). Conversely, the lowest K concentrations under 50% DI, particularly in intercropping systems, reflect increased nutrient depletion driven by combined water stress and interspecific competition (W. Tan et al., 2020). The generally higher nutrient levels observed in 2024/25 compared with 2023/24 further indicate seasonal variability, likely driven by differences in rainfall distribution, soil moisture conditions, and nutrient uptake efficiency ().
Calcium dynamics showed weaker interaction effects in 2024/25, indicating greater stability of Ca across treatments in that season. Declining Ca levels at maturity, particularly under intercropping and sole C. gynandra, suggest increased plant uptake during reproductive stages. Seasonal factors such as rainfall, irrigation scheduling, and soil acidity may also have contributed to Ca leaching and redistribution (; ). Higher Ca levels under 75% DI suggest reduced uptake under moderate stress, whereas lower values under 50% DI may reflect greater early-stage uptake driven by transpiration demand.
The decline in Ca in 2024/25, particularly under sole C. gynandra, may also reflect cumulative nutrient depletion due to continuous cropping and rainfall-driven leaching losses. Magnesium showed a contrasting pattern, with higher values under full irrigation, especially in sole cropping systems, indicating that adequate soil moisture enhances Mg availability and retention. This agrees with (), who reported higher Mg concentrations under wetter conditions. In contrast, Mg declined under deficit irrigation and intercropping, likely due to reduced moisture and increased nutrient competition. Notably, Mg concentrations often exceeded Ca levels, particularly under 100% and 75% irrigation, possibly due to differences in leaching behavior and irrigation water composition ().
Interaction effects also significantly influenced TEC–CEC in both seasons. In 2023/24, sole dry bean under full irrigation recorded the highest TEC–CEC value (2.33 cmol kg-¹), significantly exceeding most other treatments (Figure 4). This was followed by sole DB × 75% DI and DB + CM under 100% and 75% FI, indicating that dry bean-based systems generally maintained higher exchangeable cation levels. This improvement is likely attributable to enhanced nutrient cycling through legume root activity and residue decomposition, as well as improved nutrient retention under adequate soil moisture.
Figure 4
A similar trend was observed in 2024/25, where TEC–CEC generally declined with increasing water deficit across cropping systems (Figure 5). The reduction was most pronounced under sole DB and DB + CM, where values decreased markedly under 50% DI compared with full irrigation. For example, TEC–CEC declined from approximately 2.80 cmol kg-¹ under sole DB × 100% FI to about 2.55 cmol kg-¹ under sole DB × 50% DI (Figure 6). Likewise, DB + CM × 50% DI recorded significantly lower values (2.18 cmol kg-¹) compared with higher irrigation treatments (Figure 5).
Figure 5
Figure 6
Overall, sole dry bean under full irrigation consistently maintained the highest TEC–CEC values, indicating improved retention of exchangeable nutrient cations, while Cleome-based systems under deficit irrigation recorded the lowest values, reflecting reduced soil fertility status. These results demonstrate that severe water stress reduces soil cation retention capacity, whereas dry bean-based systems help sustain higher levels of exchangeable nutrients under variable irrigation conditions.
3.4.3 Principal component analysis on soil physicochemical chemical properties under deficit irrigation and intercropping
Principal component analysis (PCA) was used to evaluate the relationships among selected soil physicochemical properties under three irrigation levels [100%, 75%, and 50% irrigation deficit (DI)]. The first principal component (F1) accounted for 48.08% of the total variance, while the second principal component (F2) explained 15.60%, resulting in a cumulative variance of 63.68% explained by the first two components (Figure 7). Treatments located close to one another on the biplot exhibited similar soil characteristics, whereas treatments positioned further apart displayed distinct soil property profiles. Irrigation treatments located farther from the origin and in the direction of specific soil variables indicated a strong association with those properties. Calcium (Ca), potassium (K), magnesium (Mg), and soil pH were represented by long vectors oriented in a similar direction within the third quadrant, indicating a strong positive correlation among these variables. Conversely, phosphorus (P), bulk density (BD), electrical conductivity (EC), and sodium (Na) were clustered closer to the origin in the fourth quadrant, suggesting weaker correlations and lower contributions to the overall variation explained by the first two principal components. The proximity of these variables to the origin further indicates that their distribution may be influenced by additional factors not fully captured by F1 and F2.
Figure 7
Principal component analysis was also conducted to assess the relationships among selected soil properties under five cropping systems: sole-planted dry beans, dry beans × Cleome myriocarpus, sole-planted C. myriocarpus, dry beans × Cleome gynandra, and sole-planted C. gynandra (Figure 8). The PCA revealed a strong positive association among Ca, Mg, and K, as evidenced by their closely aligned vectors and positive loadings along the F2 axis. These variables contributed substantially to the separation of treatments along the vertical dimension of the biplot. In contrast, Na, EC, and pH were more strongly associated with the F1 axis and were positioned opposite to Ca, Mg, and K, indicating contrasting patterns of variation among the soil chemical properties. The separation of variables across the principal component axes suggests that the cropping systems differentially influenced nutrient availability and soil chemical characteristics. Cropping systems positioned in the direction of Ca, Mg, and K were associated with improved soil fertility attributes, whereas those aligned with Na and EC reflected relatively greater salinity-related characteristics.
Figure 8
3.5 Effect of growth stages, cropping treatments and effect irrigation levels on dry bean pod yield
3.5.1 Main effects of cropping systems and irrigation levels on dry bean pod yield
Cropping treatments had limited effects on most dry bean pod yield components during both seasons. The number of pods per plant, pod length, and pod weight were generally not significantly affected by cropping system, indicating that intercropping dry beans with Cucumis myriocarpus (CM) or Cleome gynandra (CG) did not substantially reduce pod production compared with sole dry bean cultivation (Table 6). However, significant differences were observed for pod length in 2024/25 and number of seeds per pod in 2023/24. Sole dry beans consistently produced the longest pods and the highest number of seeds per pod, while the dry bean–CG intercrop recorded lower values (Table 6). In particular, the number of seeds per pod declined from 6.36 in sole dry beans to 5.86 in the dry bean–CG intercrop during 2023/24 (Table 7). These results suggest that although intercropping had minimal effects on pod formation, competition for light, water, and nutrients in the dry bean–CG system may have slightly reduced reproductive development and seed set.
Table 7
| Parameter | No. of pods/plant | Pod length (cm) | No. of seeds/pod | Pod weight (g) | ||||
|---|---|---|---|---|---|---|---|---|
| Season | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 |
| Cropping treatments | ||||||||
| Sole dry beans | 12,278a | 10,639a | 11,022a | 10,564a | 6,361a | 7,333a | 2,718a | 2,039a |
| Dry beans x C. myriocarpus | 12,028a | 10,583a | 10,781a | 10,369ab | 6,305a | 7,236a | 2,641a | 2,002a |
| Dry beans x C. gynandra | 11,944a | 10,361a | 10,706a | 9,731b | 5,861b | 7,097a | 2,613a | 1,91a |
| HSD (Critical value for comparison) | 0,724 | 0,712 | 0,745 | 0,706 | 0,438 | 0,755 | 0,322 | 0,234 |
| Response | ns | ns | ns | * | * | ns | ns | ns |
| Irrigation levels | ||||||||
| 100% FI | 13,694a | 11,569a | 12,25a | 11,792a | 6,888a | 8,291a | 2,937a | 2,129a |
| 75% DI | 12,917a | 11,528a | 11,964a | 11,386a | 6,555a | 7,055ab | 2,908a | 1,913a |
| 50% DI | 9,639b | 8,486b | 8,295b | 7,486b | 5,083b | 6,319b | 2,128b | 1,909a |
| HSD (Critical value for comparison) | 0,936 | 1,112 | 0,996 | 0,954 | 0,665 | 1,254 | 0,288 | 0,298 |
| Response | * | * | * | * | * | * | ns | ns |
Effect of cropping treatments and effect irrigation levels on dry bean pod yield during 2023/24 season (experiment 1) and 2024/25 season (experiment 2).
DB, Dry beans; CM, Cucumis myriocarpus; CG, Cleome gynandra; FI, Full irrigation; DI, Deficit irrigation. Column means followed by the same letter were not significantly different (p < 0.05) according to Tukey’s HSD (Honestly Significant Difference) Test. Response **(Highly significant p ≤ 0.01), *(Significant p ≤ 0.05), ns (not significant p ≥ 0.05). Bold values are mean values for honest significant difference.
Irrigation level significantly influenced most pod yield components in both seasons. Full irrigation (100% FI) produced the highest number of pods per plant, longest pods, and greatest number of seeds per pod, whereas severe water deficit (50% DI) resulted in the lowest values. For example, pod number declined from 13.69 to 9.64 pods plant-¹ in 2023/24 and from 11.57 to 8.49 pods plant-¹ in 2024/25 as irrigation decreased from 100% FI to 50% DI (Table 7). Similarly, pod length and seed number per pod decreased significantly under severe deficit irrigation. These reductions can be attributed to limited soil moisture reducing photosynthesis, flower retention, pod development, and assimilate translocation to reproductive organs. Pod weight, however, was not significantly affected by irrigation in either season, indicating that dry beans maintained individual pod mass despite reductions in pod number and seed production under water stress. Overall, adequate irrigation was essential for maximizing pod yield components, while severe water deficits negatively affected reproductive performance.
3.5.1 Interaction effects of cropping systems × irrigation levels on dry bean pod yield
The interaction between cropping systems and irrigation levels significantly affected pod length (2024/25), number of seeds per pod (2023/24), and pod weight (2023/24), while no significant interaction was observed for pod number per plant. Across both seasons, sole dry beans under full irrigation generally produced the highest values for most yield components. In 2024/25, sole dry bean × 100% FI recorded the longest pods (12.24 cm), whereas the shortest pods were observed under dry bean + CG × 50% DI (7.21 cm) (Table 8). Similarly, during 2023/24, sole dry bean × 100% FI produced the highest number of seeds per pod (7.0), while dry bean+CG × 50% DI recorded the lowest value (4.42) (Table 8). These findings indicate that the negative effects of water stress on reproductive development were amplified when dry beans were grown in association with C. gynandra.
Table 8
| Parameter | No. of pods/plant | Pod length (cm) | No. of seeds/pod | Pod weight (g) | ||||
|---|---|---|---|---|---|---|---|---|
| Season | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 |
| Sole dry beans x 100% FI | 14,083a | 11,958a | 12,4a | 12,242a | 7a | 8,541a | 2,988a | 2,14a |
| Sole dry beans x 75% DI | 13,667a | 11,792a | 12,283a | 11,992a | 6,916ab | 7,833a | 2,975a | 2,136a |
| Sole dry beans x 50% DI | 13,333a | 11,667a | 12,192a | 11,992a | 6,75ab | 6,458a | 2,960a | 2,111a |
| DB + CM x 100% FI | 13,167a | 11,625a | 12,158a | 11,658a | 6,667abc | 8,25a | 2,890ab | 1,994a |
| DB + CM x 75% DI | 13,083a | 11,292ab | 11,833a | 11,142a | 6,583abc | 7a | 2,864ab | 1,994a |
| DB + CM x 50% DI | 12,5a | 10,958abc | 11,775a | 10,508a | 6,416abc | 6,875a | 2,858ab | 1,987a |
| DB + CG x 100% FI | 9,917b | 8,833bcd | 8,592b | 7,708b | 5,583bcd | 8,083a | 2,192bc | 1,902a |
| DB + CG x 75% DI | 9,583b | 8,458cd | 8,351b | 7,542b | 5,25cd | 6,333a | 2,106c | 1,849a |
| DB + CG x 50% DI | 9,417b | 8,167d | 7,942b | 7,208b | 4,416d | 5,62a | 2,085c | 1,740a |
| HSD (Critical value for comparison) | 2,177 | 2,585 | 2,317 | 2,217 | 1,546 | 2,916 | 0,671 | 0,693 |
| Response | ns | * | ns | ns | * | ns | * | ns |
Interactive effect of cropping systems x irrigation levels on dry bean pod yield during 2023/24 season (experiment 1) and 2024/25 season (experiment 2).
DB, Dry beans; CM, Cucumis myriocarpus; CG, Cleome gynandra; FI, Full irrigation; DI, Deficit irrigation. Column means followed by the same letter were not significantly different (p < 0.05) according to Tukey’s HSD (Honestly Significant Difference) Test. Response **(Highly significant p ≤ 0.01), *(Significant p ≤ 0.05), ns (not significant p ≥ 0.05). Bold values are mean values for honest significant difference.
The interaction results further demonstrate that deficit irrigation reduced yield components within each cropping system, although the magnitude of reduction varied among intercrops. The dry bean–CM intercrop-maintained pod yield characteristics relatively close to those of sole dry beans across irrigation levels, suggesting lower interspecific competition and greater compatibility between the two species (Table 7). In contrast, the dry bean + CG intercrop consistently produced lower pod numbers, shorter pods, fewer seeds per pod, and lower pod weights, particularly under 50% DI (Table 8). The significant reductions observed under dry bean + CG × 50% DI indicate that combining severe water deficit with interspecific competition imposed greater stress on dry bean reproductive growth. Overall, the interaction results suggest that while dry beans can be successfully intercropped without substantial yield penalties, maintaining adequate irrigation is critical, and the dry bean–CM intercrop appears more resilient to water stress than the dry bean + CG system.
3.6 Effect of growth stages, cropping treatments and effect irrigation levels on C. myriocarpus and C. gynandra growth and yield
3.6.1 Main effects of cropping treatments and irrigation levels on Cucumis myriocarpus and Cleome gynandra growth and yield
Cropping treatments significantly influenced the growth and biomass production of both Cucumis myriocarpus and Cleome gynandra, although the magnitude and direction of the responses varied between the crops. For C. myriocarpus, sole cropping consistently resulted in greater root weight, shoot weight, and biomass accumulation than intercropping with dry beans, indicating that interspecific competition reduced overall growth performance (Table 9). Similarly, sole C. gynandra produced greater biomass than the intercrop treatment, despite intercropped plants exhibiting higher root and shoot weights in some instances (Table 10). These findings suggest that when grown alone, both indigenous vegetables had unrestricted access to water, nutrients, and light, enabling greater biomass accumulation. In contrast, intercropping with dry beans introduced competition for resources, resulting in reduced dry matter production. However, stem length and plant height were generally unaffected by cropping treatments in both species, indicating that plant architecture was less sensitive to intercropping than biomass-related traits (Tables 9, 11).
Table 9
| Parameter | Root weight (g) | Shoot weight (g) | Biomass weight(g) | Stem length (cm) | Plant height (cm) | |||||
|---|---|---|---|---|---|---|---|---|---|---|
| Season experiment | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 |
| Growth stages | ||||||||||
| 50% flowering stage | 6,955b | 4,889a | 9,394a | 6,877ab | 200,56ab | 208,5b | 3,807a | 2,805a | 40,5ab | 39,33a |
| Early pod stage | 7,5a | 4,157b | 9,366a | 8,472a | 215,11a | 240,56a | 3,644a | 2,483a | 42,333a | 37,5ab |
| Mid-pod stage | 6,494c | 3,971bc | 9,205a | 6,516ab | 189,78b | 184,06c | 3,405ab | 2,415a | 37,889ab | 34,88ab |
| Maturity stage | 6,25c | 3,653c | 6,827b | 5,283b | 150,56c | 159,83d | 2,966b | 1,835b | 35,556b | 32,5b |
| HSD (Critical value for comparison) | 0,347 | 0,444 | 0,831 | 2,291 | 17,103 | 14,528 | 0,511 | 0,560 | 6,294 | 5,984 |
| Response | ** | ** | ns | ** | ** | ** | ns | ns | ** | ** |
| Cropping treatments | ||||||||||
| Sole C. myriocarpus | 7,080a | 4,567a | 8,933a | 7,122a | 195,22a | 201,64a | 3,560a | 2,49a | 39,917a | 36,917a |
| C. myriocarpus + dry beans | 6,519b | 3,768b | 8,463b | 6,452b | 182,78b | 194,83b | 3,351a | 2,27a | 38,222a | 35,194a |
| HSD (Critical value for comparison) | 0,1855 | 0,214 | 0,4437 | 0,349 | 9,1253 | 5,5611 | 0,273 | 0,2217 | 3,3585 | 3,3673 |
| Response | * | * | * | * | * | * | ns | ns | ns | ns |
| Irrigation levels | ||||||||||
| 100% FI | 7,137a | 4,462a | 9,241a | 7,537a | 207,54a | 205,79a | 3,736a | 2,7a | 40,875a | 37,83a |
| 75% DI | 7,020a | 4,324a | 9,029a | 6,616ab | 182,04b | 196,42b | 3,466ab | 2,31a | 39,125a | 36,125a |
| 50% DI | 6,241b | 3,716b | 7,825b | 6,208b | 177,42b | 192,5b | 3,1642b | 2,14a | 37,208a | 34,208a |
| HSD (Critical value for comparison) | 0,273 | 0,311 | 0,654 | 1,155 | 13,459 | 7,604 | 0,402 | 0,564 | 4,953 | 5,886 |
| Response | * | * | * | * | * | * | * | ns | ns | ns |
Effect of cropping treatments and effect irrigation levels on C. myriocarpus growth and yield during 2023/24 season (experiment 1) and 2024/25 season (experiment 2).
CM, Cucumis myriocarpus; DB, Dry beans; FI, Full irrigation; DI, Deficit irrigation. Column means followed by the same letter were not significantly different (p < 0.05) according to Tukey’s HSD (Honestly Significant Difference) Test. Response **(Highly significant p ≤ 0.01), *(Significant p ≤ 0.05), ns (not significant p ≥ 0.05). Bold values are mean values for honest significant difference.
Table 10
| Parameter | Root weight (g) | Shoot weight (g) | Biomass weight (g) | Stem length (cm) | Plant height (cm) | |||||
|---|---|---|---|---|---|---|---|---|---|---|
| Season | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 |
| Sole dry beans x 100% FI | 7,583a | 5,349a | 9,575a | 5,925a | 209,92a | 205,92a | 4,173a | 2,88a | 42,33a | 39,25a |
| Sole dry beans x 75% DI | 7,032b | 4,510b | 9,4a | 4,467b | 205,17a | 205,67a | 3,569ab | 2,521ab | 41,33a | 38,33a |
| Sole dry beans x 50% DI | 7,008b | 4,138bc | 9,083ab | 4,4bc | 192,33ab | 202,42ab | 3,364b | 2,388ab | 39,417a | 36,417a |
| DB + CM x 100% FI | 6,691b | 3,841cd | 8,483abc | 3,867bcd | 188,17abc | 196,83abc | 3,3b | 2,238ab | 39a | 36a |
| DB + CM x 75% DI | 6,65b | 3,591cd | 8,141bc | 3,75cd | 171,75bc | 190,42bc | 3,185b | 2,220ab | 36,917a | 33,917a |
| DB + CM x 50% DI | 5,833b | 3,575d | 7,508c | 3,75d | 166,67c | 188,17c | 3,143b | 2,063b | 35,417a | 32,417a |
| HSD (Critical value for comparison) | 0,475 | 0,556 | 1,136 | 0,623 | 23,361 | 14,467 | 0,698 | 0,576 | 8,597 | 8,759 |
| Response | ns | * | * | * | * | * | ns | * | ns | ns |
Interactive effect of cropping systems x irrigation levels on C. myriocarpus growth and yield during 2023/24 season (experiment 1) and 2024/25 season (experiment 2).
CM, Cucumis myriocarpus; DB, Dry beans; FI, Full irrigation; DI, Deficit irrigation. Column means followed by the same letter were not significantly different (p < 0.05) according to Tukey’s HSD (Honestly Significant Difference) Test. Response **(Highly significant p ≤ 0.01), *(Significant p ≤ 0.05), ns (not significant p ≥ 0.05). Bold values are mean values for honest significant difference.
Table 11
| Parameter | Sub-stomata (vpm) | Stomatal conductance(mol m-2 s-1) | CO2 (vpm) | Photosynthesis rate(µmol m-2 s-1) | Transpiration rate(mmol m-2 s-1) | Chlorophyll content(SPAD) | ||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Season | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 |
| Cropping treatments | ||||||||||||
| Sole C. gynandra | 36,5a | 47,625a | 0,359a | 0,840a | 424,79a | 423,58a | 11,322a | 11,576b | 4,334a | 2,465a | 17,778a | 14,389a |
| C. gynandra + Dry beans | 43,375a | 46,306a | 0,327a | 0,220a | 423,58a | 423,39a | 14,012a | 14,239a | 3,841a | 2,032b | 16,194b | 13,972a |
| HSD (Critical value for comparison) | 9,903 | 2,211 | 0,132 | 1,217 | 2,057 | 3,870 | 3,238 | 0,666 | 0,551 | 0,119 | 0,810 | 1,210 |
| Response | ns | ns | ns | ns | ns | ns | ns | * | ns | * | * | ns |
| Irrigation levels | ||||||||||||
| 100% FI | 45,313a | 51,813a | 0,496a | 1,159a | 425a | 424,44a | 13,947a | 13,834a | 3,779b | 2,315a | 19,25a | 15,75a |
| 75% DI | 41,688a | 44,958b | 0,3725a | 0,225a | 424,69a | 423,69a | 13,502a | 12,975a | 5,0894a | 2,228a | 18,792a | 14,667a |
| 50% DI | 32,813a | 44,125b | 0,1607b | 0,205a | 422,87a | 422,33a | 10,551a | 11,913b | 3,395b | 2,202a | 12,917b | 12,125b |
| HSD (Critical value for comparison) | 14,607 | 3,262 | 0,195 | 1,795 | 3,034 | 5,708 | 4,776 | 0,983 | 0,813 | 0,176 | 1,194 | 1,785 |
| Response | ns | * | * | ns | ns | ns | ns | * | * | ns | * | * |
Effect of cropping treatments and effect irrigation levels on C. gynandra physiology during 2023/24 season (experiment 1) and 2024/25 season (experiment 2).
CG, Cleome gynandra, DB, Dry beans; FI, Full irrigation; DI, Deficit irrigation; CO2, Carbon dioxide; BD, bulk density; EC, Electrical conductivity. Column means followed by the same letter were not significantly different (p < 0.05) according to Tukey’s HSD (Honestly Significant Difference) Test. Response **(Highly significant p ≤ 0.01), *(Significant p ≤ 0.05), ns (not significant p ≥ 0.05). Bold values are mean values for honest significant difference.
Irrigation level significantly affected most growth and yield parameters of both crops. Full irrigation (100% FI) consistently produced the highest root weight, shoot weight, biomass accumulation, and stem length, whereas 50% deficit irrigation (DI) resulted in the lowest values (Tables 9, 11). In both species, reductions in biomass under severe water deficit reflected the negative effects of moisture stress on photosynthesis, nutrient uptake, cell expansion, and assimilate production. Although growth parameters declined with decreasing irrigation levels, plant height remained relatively stable across irrigation treatments, suggesting that both species maintained vertical growth despite reduced biomass accumulation. Overall, the results demonstrate that adequate soil moisture is essential for maximizing growth and biomass production in both C. myriocarpus and C. gynandra, while severe water deficits significantly constrain vegetative development and productivity.
3.6.2 Interaction effects of cropping systems × irrigation levels on Cucumis myriocarpus and Cleome gynandra growth and yield
The interaction between cropping systems and irrigation levels significantly influenced most growth and biomass parameters of both species, demonstrating that crop performance depended on the combined effects of water availability and cropping arrangement. Across both seasons, the highest values for root weight, shoot weight, biomass accumulation, and stem length were generally recorded under sole cropping combined with full irrigation, while the lowest values occurred under intercropping subjected to severe water deficit (50% DI). For both C. myriocarpus and C. gynandra, biomass production declined progressively as irrigation levels decreased, with the greatest reductions occurring when plants were grown together with dry beans under water-limited conditions (Tables 10, 12). These findings indicate that the benefits of adequate irrigation were maximized under sole cropping systems where competition for resources was minimal.
Table 12
| Parameter | Sub-stomata (vpm) | Stomatal conductance (mol m-2 s-1) | CO2 (vpm) | Photosynthesis rate(µmol m-2 s-1) | Transpiration rate(mmol m-2 s-1) | Chlorophyll content(SPAD) | ||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Season | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 | 2023/24 | 2024/25 |
| Sole SG x 100% FI | 50,25a | 45,25bc | 0,7175a | 2,086a | 427,12a | 427,04a | 15,278a | 15,758a | 6,202a | 2,376a | 18,833ab | 17,417a |
| Sole SG x 75% DI | 41,75a | 44,667bc | 0,5613ab | 0,233a | 425,12a | 425,58a | 14,482a | 14,192ab | 4,223b | 2,254ab | 18,75b | 15,167ab |
| Sole SG x 50% DI | 41,625a | 44,25c | 0,1756c | 0,230a | 424,25a | 422,75a | 13,411a | 12,768bc | 3,977b | 1,979bc | 17,667b | 12,667bc |
| SG + DB x 100% FI | 40,375a | 53,625a | 0,275bc | 0,220a | 423a | 421,92a | 12,275a | 11,91cd | 3,466b | 2,542a | 13,75c | 14,167bc |
| SG + DB x 75% DI | 38,25a | 50ab | 0,183c | 0,206a | 422,87a | 421,83a | 11,726a | 11,758cd | 3,33b | 2,477a | 20,833a | 14,083bc |
| SG + DB x 50% DI | 27,375a | 44c | 0,145c | 0,204a | 422,75a | 421,79a | 8,828a | 11,058d | 3,324b | 1,863c | 12,083c | 11,583c |
| HSD (Critical value for comparison) | 25,353 | 5,661 | 0,338 | 3,117 | 5,266 | 9,908 | 8,290 | 1,706 | 1,412 | 0,305 | 2,074 | 3,099 |
| Response | ns | * | * | ns | ns | ns | ns | ** | ns | * | ** | * |
Interactive effect of cropping systems x irrigation levels on C. gynandra physiology during 2023/24 season (experiment 1) and 2024/25 season (experiment 2).
SG, Cleome gynandra; DB, Dry bean; FI, Full irrigation; DI, Deficit irrigation; CO2, Carbon dioxide; BD, bulk density; EC, Electrical conductivity. Column means followed by the same letter were not significantly different (p < 0.05) according to Tukey’s HSD (Honestly Significant Difference) Test. Response **(Highly significant p ≤ 0.01), *(Significant p ≤ 0.05), ns (not significant p ≥ 0.05). Bold values are mean values for honest significant difference.
The interaction effects further revealed that water stress intensified the competitive effects associated with intercropping. Although reductions in growth occurred under deficit irrigation in all treatments, the decline was generally greater in intercropped plants than in sole crops, indicating that limited water availability increased competition for soil moisture and nutrients. Stem length also decreased significantly under combined intercropping and deficit irrigation, reflecting restricted vegetative growth under resource-limited conditions. However, plant height remained relatively stable across treatment combinations for both species, suggesting that height growth was less sensitive to competition and moisture stress than biomass-related traits (Tables 10, 12). Overall, the results demonstrate that sole cropping under full irrigation provided the most favorable conditions for growth and biomass production of both C. myriocarpus (Table 9) and C. gynandra (Table 12), whereas the combination of intercropping and severe water deficit imposed the greatest constraints on plant development.
4 Discussion
4.1 Water use efficiency and soil moisture
The observed variation in instantaneous water-use efficiency (iWUE) across growth stages highlights the strong dependence of plant water-use regulation on phenological development, crop type, and seasonal environmental conditions. The absence of significant interaction effects between cropping system and irrigation in most growth stages suggests that iWUE was primarily governed by direct physiological responses to water availability and inherent crop-specific traits, rather than strong synergistic effects between cropping arrangement and irrigation regime.
The influence of irrigation on iWUE varied markedly across developmental stages, with the mid-pod stage emerging as the most sensitive period. This stage corresponds to peak reproductive development, when sink demand is high and assimilate partitioning to developing pods is maximised (). Under adequate water supply, sustained stomatal conductance and efficient mesophyll CO2 diffusion enhance photosynthetic carbon gain relative to transpiration, thereby improving iWUE (). In contrast, the seed-filling or maturity stage showed a weaker response to irrigation, particularly in the 2023/24 season, because most structural and physiological components are already fully developed, limiting the capacity of additional water to further enhance photosynthetic activity. Similar stage-dependent responses have been widely reported, with reproductive stages consistently identified as the most water-sensitive phases in legumes due to tight coupling between carbon assimilation and transpiration processes (; ).
Seasonal variation further modified the irrigation response, with the 2023/24 season showing comparatively higher iWUE during seed filling, while the 2024/25 season exhibited stronger responses during the mid-pod stage. These differences are likely driven by interannual variation in environmental drivers such as temperature, solar radiation, vapour pressure deficit, and rainfall distribution, all of which regulate stomatal behavior and plant water status (Xiong and Nadal, 2020). Under higher evaporative demand, stomatal regulation becomes more sensitive to soil moisture availability, thereby amplifying treatment differences, whereas under milder atmospheric demand these differences are often buffered ().
Across studies, it is widely reported that iWUE typically increases under moderate water deficit due to partial stomatal closure, which reduces transpiration more strongly than photosynthetic carbon assimilation (Yang Y. J. et al., 2021; ). However, under severe water stress, reductions in stomatal conductance, mesophyll CO2 diffusion, and photosynthetic enzyme activity collectively limit carbon assimilation, resulting in reduced iWUE (). This non-linear response has been consistently observed in both legumes and cereal systems, confirming that optimal water-use efficiency occurs under moderate, non-lethal stress conditions rather than under either full irrigation or severe deficit (; ). Adequate irrigation during reproductive stages is therefore essential to maintain stomatal function and mesophyll CO2 diffusion, ensuring efficient carbon gain per unit water transpired ().
Cropping system effects on iWUE reported in the literature are generally more variable and strongly dependent on species compatibility and resource partitioning. Some intercropping systems improve water-use efficiency through complementary rooting depth, enhanced soil moisture extraction, and improved canopy light distribution, whereas others reduce efficiency due to intensified competition for water, nutrients, and light (; Wen et al., 2022). In legume-based systems, dry beans often exhibit higher iWUE than leafy vegetables due to more conservative water-use strategies and stronger stomatal regulation under stress conditions (Wen et al., 2022). However, when interspecific competition is strong, especially under water deficit, stomatal closure may occur earlier, reducing carbon assimilation and limiting iWUE gains ().
The convergence of iWUE values among dry bean treatments at physiological maturity reflects the natural decline in physiological activity as crops approach senescence (; de ). Reduced leaf area, declining chlorophyll content, and lower stomatal conductance collectively reduce carbon assimilation and transpiration, diminishing physiological differences among treatments (Wang et al., 2026). Therefore, results aligned with broader physiological evidence that iWUE is primarily regulated by phenological stage and atmospheric demand, with irrigation acting as a modulating factor and cropping system effects being secondary and highly context dependent. The mid-pod stage represents a critical window for optimizing water-use efficiency, where adequate but not excessive irrigation is essential to sustain the balance between carbon assimilation and water loss.
4.2 Soil physicochemical properties
The present study demonstrated that both cropping system and irrigation level significantly influenced soil physicochemical properties, nutrient availability, and cation dynamics. The significant interaction effects observed between these factors indicate that soil responses were regulated not only by water availability but also by plant-mediated processes occurring within the rhizosphere (B. Tan et al., 2021). In particular, dry bean-based intercropping systems improved several soil chemical properties relative to sole cropping systems. These improvements may be attributed to greater belowground functional diversity, which enhances nutrient cycling through complementary resource acquisition, increased rhizodeposition, and stimulation of microbial activity (). Diverse rooting systems create heterogeneous rhizosphere environments that support a wider range of microbial communities and nutrient transformation pathways, thereby improving nutrient availability and soil fertility ().
Bulk density is a key indicator of soil physical quality as it influences root penetration, water infiltration, aeration, and microbial activity (). The lower bulk density observed under intercropping, particularly under 75% deficit irrigation, suggests that moderate water stress combined with diverse rooting systems improved soil structural stability. Intercropped species often develop complementary root architectures that create biopores and improve pore continuity throughout the soil profile (). In addition, greater belowground biomass and root turnover contribute organic inputs that stimulate microbial activity and the production of extracellular polysaccharides, which act as binding agents that stabilize soil aggregates and improve soil structure (Xu et al., 2021). This biological reinforcement of soil aggregates improves moisture retention and buffering capacity under fluctuating irrigation regimes, helping to reduce the severity of water stress during dry periods. Similar processes have been reported in cereal–legume intercropping systems, where enhanced root–soil–microbe interactions improved soil aggregation, reduced bulk density, and supported more stable yield performance under water-limited conditions (). This contrasts with systems dominated by single crops, where lower diversity of residue inputs can limit aggregate stability and slow structural recovery. Importantly, these soil physical improvements help explain why intercropping systems, particularly dry bean + Cucumis myriocarpus, maintained relatively comparable yield and iWUE performance under moderate deficit irrigation, despite increased interspecific competition. In semi-arid environments, where water availability is the primary limiting factor, improved soil structure enhances infiltration, water storage, and root access to deeper moisture reserves, thereby buffering plants against drought stress and stabilising physiological processes (). In contrast, soils with higher bulk density restrict root expansion and reduce hydraulic conductivity, thereby intensifying drought stress and limiting stomatal regulation and carbon assimilation ().
The higher soil pH observed under intercropping systems suggests that species interactions modified rhizosphere biogeochemical processes regulating soil acid–base balance. Legumes such as dry bean influence soil pH through multiple mechanisms, including biological nitrogen fixation, selective nutrient uptake, and associated proton exchange reactions in the rhizosphere (). During nitrogen fixation and nitrate uptake, legumes may release or consume H+ ions depending on the dominant nitrogen form and uptake pathway, thereby directly influencing soil acidity (). In addition, the decomposition of diverse crop residues in intercropping systems enhances the release of basic cations such as Ca, Mg, and K, which contribute to buffering soil acidity and stabilising pH levels (Țopa et al., 2025). Greater crop diversity also tends to stimulate microbial activity, accelerating organic matter turnover and nutrient mineralisation, which further influences soil acid–base dynamics through the production and consumption of organic acids and base-forming ions (Usman et al., 2020). However, it is important to note that reported pH responses in intercropping systems are not always consistent across studies. For example (), reported higher soil pH in monocropped legumes compared with intercrops, indicating that the direction and magnitude of pH change are highly context dependent. Such contrasting findings are typically attributed to differences in soil buffering capacity, crop species composition, residue chemistry (e.g., C:N ratio and lignin content), fertiliser inputs, and prevailing environmental conditions, all of which regulate proton dynamics and cation exchange processes in soil systems (; T. Yang et al., 2020). In the present study, the observed increase in pH under dry bean–African leafy vegetable intercrops suggests that complementary nutrient cycling and enhanced residue turnover promoted a more buffered and chemically stable rhizosphere environment, particularly under conditions of variable irrigation ().
Electrical conductivity (EC) responded strongly to both irrigation level and cropping system, with generally higher values observed under deficit irrigation and intercropping treatments. This pattern is consistent with studies reporting that reduced irrigation enhances salt accumulation in the upper soil profile due to limited leaching and increased evapotranspiration under high evaporative demand (). Similar findings have been reported by (; ), who attributed increased EC under deficit irrigation to restricted percolation and greater retention of soluble salts within the root zone. Likewise, (H. N. ) highlighted that under full irrigation, continuous downward water flux promotes salt displacement beyond the rooting depth, thereby reducing EC in the topsoil. These consistent findings across studies reinforce the hydrological control of EC dynamics in irrigated agroecosystems (M. Wang et al., 2024).
However, the magnitude and ecological implications of EC increases vary considerably among studies. In some long-term deficit irrigation experiments, EC has been reported to reach levels associated with salinity stress and yield suppression, particularly under arid conditions with poor drainage (Toumi et al., 2024; Soni et al., 2021). In contrast, the present study showed that although EC increased under deficit irrigation, values remained below salinity thresholds associated with crop growth inhibition. This agrees with findings from semi-arid intercropping systems where short- to medium-term water deficits led to nutrient concentration effects without triggering true salinisation, suggesting that EC elevation in such contexts may primarily reflect transient solute concentration rather than harmful salt accumulation (; Yan-Jun et al., 2025).
Cropping system effects observed in this study also align with and partly extend previous research. Similarly intercropping was associated with higher EC, likely due to enhanced root density, increased rhizodeposition, and accelerated residue decomposition, all of which contribute to greater release of soluble ions into the soil solution. In agreement, other studies on legume-based intercropping systems have shown that increased biological activity and nutrient cycling can elevate EC through intensified mineralisation processes (). However, contrasting evidence also exists, where intercropping reduced EC due to improved nutrient uptake efficiency and reduced residual ion accumulation in soils with highly competitive root systems (). These differences highlight that the direction of EC response depends strongly on species combination, rooting depth complementarity, and water availability (Song et al., 2025).
Moderate increases in electrical conductivity (EC) can enhance nutrient availability by increasing the concentration of soluble ions in the soil solution, thereby improving nutrient accessibility to plants. However, when EC exceeds critical salinity thresholds, excessive salt accumulation raises the osmotic potential of the soil solution, making water uptake by roots more energetically demanding and potentially reducing plant growth and productivity (; ). In the present study, deficit irrigation promoted EC through reduced leaching and the concentration of dissolved ions within the root zone, a response commonly observed in semi-arid irrigated systems.
Water stress may also suppress root elongation, root hair development, and phosphatase activity, further restricting phosphorus acquisition. Consequently, phosphorus uptake declines and greater quantities remain in the soil (). Similar findings were reported by (), who observed increased residual soil phosphorus under moderate water deficit in quinoa due to reduced plant uptake. Likewise (Suriyagoda et al., 2014), demonstrated that drought conditions significantly limit phosphorus mobility and uptake because diffusion processes are highly dependent on soil moisture. In contrast, other studies have reported lower available phosphorus under water deficit conditions due to reduced microbial activity, slower mineralisation rates, and decreased phosphorus solubilisation in dry soils (Zhang et al., 2020; ). For example (Zhang et al., 2020), found that prolonged drought reduced available phosphorus in maize systems by suppressing microbial-mediated nutrient cycling. These contrasting findings suggest that the response of soil phosphorus to deficit irrigation depends on the balance between reduced plant uptake and reduced phosphorus release from soil organic and mineral pools.
The lower residual phosphorus concentrations observed under intercropping systems suggest more efficient phosphorus acquisition and utilisation. This finding agrees with (), who reported enhanced phosphorus uptake in cereal–legume intercrops due to complementary rooting patterns and rhizosphere facilitation. Similarly (Xu et al., 2021), found lower residual phosphorus concentrations in intercropped systems because diverse root systems exploited a greater soil volume and improved nutrient acquisition. Intercropped species often differ in rooting depth, root architecture, and rhizosphere chemistry, enabling them to access phosphorus from different soil microsites (Wang et al., 2023b). Furthermore, root exudates such as organic acids can mobilise sparingly soluble phosphorus fractions, while microbial associations enhance phosphorus solubilisation and uptake (; ). However, not all studies report positive effects of intercropping on phosphorus acquisition. Some researchers have found that interspecific competition may exceed facilitative effects, particularly under severe water stress or nutrient-poor conditions, leading to lower phosphorus uptake efficiency than expected (Zhu et al., 2022; ). These contrasting responses indicate that the balance between facilitation and competition is highly dependent on species compatibility, soil fertility, and water availability.
The observed differences in exchangeable cation concentrations further highlight the importance of plant–soil–water interactions. Dry bean-based systems generally maintained higher concentrations of K, Ca, and Mg than Cleome-based systems, indicating differences in nutrient uptake patterns and nutrient cycling efficiency. Similar findings were reported by (M. ), who found that legume-based systems improved nutrient retention through biological nitrogen fixation and enhanced microbial activity. Likewise (), showed that legume residues decompose more rapidly than non-leguminous residues because of their lower C:N ratio, accelerating nutrient mineralisation and replenishment of exchangeable nutrient pools. This suggests that the net effect of legumes on soil nutrient status depends on the balance between nutrient inputs from residue decomposition and nutrient removal through plant uptake.
The decline in exchangeable cation concentrations under deficit irrigation can be explained by reduced nutrient mobility in dry soils. Potassium moves primarily through diffusion and mass flow, while calcium and magnesium are transported mainly through mass flow associated with water movement toward roots (R. H. ). As soil moisture declines, nutrient diffusion coefficients decrease and nutrient replenishment within the rhizosphere becomes increasingly restricted (Sardans and Peñuelas 2021). Similar reductions in K, Ca, and Mg availability under drought stress have been reported by (), who observed that reduced soil moisture limited nutrient transport and root interception in several field crops. Likewise (), reported that water deficits reduced potassium mobility and uptake because of restricted diffusion pathways. However, some studies have reported increased concentrations of exchangeable cations under deficit irrigation due to concentration effects arising from reduced leaching losses and lower nutrient uptake (Wang et al., 2023a). For example (), found that moderate water deficits increased soil concentrations of certain cations because less water moved nutrients beyond the root zone. The higher cation concentrations observed under full irrigation in the present study suggest that the beneficial effects of improved nutrient mobility and root uptake exceeded any concentration effects associated with reduced leaching.
The interaction between cropping system and irrigation further suggests that facilitation and competition operated simultaneously within the intercrop systems. Under moderate water deficits, complementary rooting patterns may have improved resource capture and nutrient-use efficiency by reducing direct competition for soil resources. This is consistent with the findings of (), who demonstrated that niche differentiation in intercrops can enhance resource acquisition through spatial and temporal complementarity (). Similarly (), reported that species diversity often improves nutrient capture through facilitative rhizosphere interactions and greater belowground functional diversity (). However, as water availability declined further, competition for limited soil moisture and nutrients likely intensified, reducing the benefits of facilitation. The reduced nutrient concentrations and yield performance recorded in the dry bean + Cleome gynandra intercrop under severe deficit irrigation likely reflect a shift from facilitative to competitive interactions, whereby increasing resource limitation intensified competition for water and nutrients, thereby constraining crop productivity ().
Potassium concentrations declined under increasing water deficit, reflecting the strong dependence of K mobility on soil moisture. Potassium reaches roots through both diffusion and mass flow, processes that are substantially restricted under dry soil conditions (). Reduced water content decreases ion mobility and limits replenishment of potassium in the rhizosphere, thereby reducing nutrient availability (). Similar declines in soil and plant K concentrations under drought conditions have been reported by (). The higher K concentrations observed in the dry bean–Cleome intercrop may indicate enhanced nutrient mobilisation through complementary rooting patterns and increased rhizosphere activity, although interspecific competition for nutrients may become more pronounced under severe water deficits.
Calcium and magnesium exhibited similar responses, with concentrations generally declining under greater water limitation. Because both nutrients are transported primarily through mass flow, reduced soil moisture limits their movement towards roots and restricts nutrient uptake. Seasonal declines in Ca, particularly under intercropping and sole Cleome treatments, likely reflect increased nutrient removal by crops as well as possible leaching and redistribution processes associated with seasonal rainfall patterns (). The lower concentrations observed under severe deficit irrigation may also indicate greater competition for available nutrient pools among neighboring plants. Similar moisture-dependent reductions in Ca and Mg availability have been reported in drought-prone agroecosystems ().
Sodium dynamics differed somewhat from those of other cations. Higher Na concentrations under full irrigation suggest that irrigation water may have contributed additional soluble salts to the soil profile over time. This finding contrasts with (), who reported greater Na accumulation under deficit irrigation due to reduced leaching. The discrepancy may reflect differences in irrigation water quality, soil texture, drainage conditions, and climatic factors among studies (). Such contrasting findings highlight the importance of site-specific factors in determining salinity responses to irrigation management. Nevertheless, Na concentrations remained below levels associated with sodicity-related soil degradation, indicating that sodium accumulation was not sufficient to adversely affect soil quality during the study period.
The PCA results further confirmed that irrigation level and cropping system jointly influenced soil chemical dynamics. Calcium, magnesium, potassium, and pH were the principal variables driving treatment separation under different irrigation regimes, whereas the major distinction among cropping systems reflected contrasts between nutrient-related variables (Ca, Mg, and K) and salinity-related indicators (Na and EC). These patterns suggest that irrigation primarily regulated nutrient mobility and salt accumulation processes, while cropping systems influenced nutrient cycling and rhizosphere-mediated nutrient transformations.
Overall, the results demonstrate that soil nutrient dynamics were regulated by complex interactions among water availability, plant species composition, and rhizosphere processes. Deficit irrigation altered nutrient mobility and promoted the concentration of soluble salts through reduced leaching and evapoconcentration, whereas intercropping modified nutrient cycling through complementary rooting strategies, enhanced microbial activity, and improved resource capture. The superior performance of dry bean-based systems highlights the important role of legumes in maintaining soil fertility and nutrient availability under semi-arid conditions. These findings are consistent with previous studies on legume-based intercropping systems but also demonstrate that the magnitude and direction of soil responses depend strongly on irrigation management and local environmental conditions.
4.3 Crop growth and yield
The present study demonstrated that both cropping system and irrigation level significantly influenced the growth, biomass production, and yield performance of dry bean, Cucumis myriocarpus, and Cleome gynandra. The magnitude of these effects was largely dependent on their interaction, indicating that crop responses to water availability were strongly mediated by cropping arrangement. Intercropping dry bean with C. myriocarpus and C. gynandra generally reduced pod yield and associated yield components relative to sole dry bean. These reductions were most pronounced under severe deficit irrigation, suggesting that competition for soil moisture, nutrients, and light intensified as resource availability declined (; ). The dry bean + C. gynandra intercrop was particularly sensitive to water stress, exhibiting lower pod length, seed number, pod weight, and overall productivity than the dry bean + C. myriocarpus system. This suggests that under severe water limitation, competitive interactions outweighed any potential facilitative benefits of intercropping. The greater reduction in productivity observed in the dry bean + C. gynandra intercrop may be associated with overlap in rooting zones and simultaneous demand for soil water and nutrients, which likely intensified resource competition during reproductive development. In contrast, the relatively better performance of the dry bean + C. myriocarpus intercrop suggests a greater degree of niche complementarity, whereby differences in root distribution, growth habit, or resource acquisition strategies reduced direct competition and improved overall resource-use efficiency. Similar findings were reported by (), who observed greater shoot growth, root development, and seed production in sole-planted chickpea compared with intercropped systems. However, other studies have reported yield advantages in legume-based intercrops when species exhibit strong spatial and temporal complementarity in resource use (; ).
Irrigation level was a dominant factor influencing crop performance across all treatments. Full irrigation and moderate deficit irrigation generally supported higher biomass production, improved reproductive development, and greater yield components, whereas severe deficit irrigation substantially reduced growth and productivity. These responses are associated with well-documented physiological constraints imposed by water stress, including reduced stomatal conductance, lower photosynthetic rates, impaired nutrient uptake, restricted cell expansion, and reduced assimilate translocation to developing reproductive organs (). Water deficit also limits nutrient mobility within the soil by reducing mass flow and diffusion processes, thereby constraining nutrient acquisition and further suppressing growth. As soil moisture declines, plants increasingly allocate resources towards survival rather than biomass production, resulting in reduced leaf area, lower canopy development, and diminished reproductive output. Although some increases in plant height were observed under moderate water deficit, this response likely reflects an adaptive elongation strategy aimed at maintaining light interception and reducing heat stress rather than enhanced productivity, as biomass accumulation and yield were generally reduced under these conditions. Similar trends were reported by (), who observed approximately 40% reductions in plant height under 50% deficit irrigation compared with full irrigation. Likewise (), reported significantly lower dry bean growth under 50% ETc than under full irrigation. In contrast, some studies have shown that moderate deficit irrigation can improve water-use efficiency without substantial yield penalties, particularly when water stress is imposed outside critical reproductive stages. The differences among studies likely reflect variation in crop species, timing and severity of water stress, soil characteristics, and climatic conditions.
The indigenous leafy vegetables responded similarly to irrigation and cropping treatments. Sole-cropped C. myriocarpus and C. gynandra consistently produced greater biomass than intercropped plants, highlighting the suppressive effect of interspecific competition on vegetative growth. The interaction between cropping system and irrigation further showed that water deficit intensified competitive effects, resulting in greater biomass reductions under intercropping than under sole cropping systems. Under water-limited conditions, reduced soil moisture not only restricts physiological processes but also increases competition for belowground resources because nutrient mobility declines and roots must exploit a smaller effective soil volume. These findings underscore the importance of adequate soil moisture in sustaining both individual crop growth and intercrop productivity. Similar responses were reported by (Zahedi et al., 2025; Yavas et al., 2024), who noted that water-stressed plants often allocate resources towards stem elongation at the expense of leaf expansion as a drought-avoidance strategy. However, reduced leaf area limits photosynthetic capacity and assimilate production, ultimately constraining dry matter accumulation and yield (). While some researchers have reported that intercropping can partially alleviate drought stress through complementary water extraction and improved resource-use efficiency, the present results indicate that such facilitative effects were insufficient to offset competition under severe water deficit, particularly in the dry bean + C. gynandra system ().
This study was primarily limited to the assessment of soil physicochemical properties, crop water-use efficiency, and crop performance under different cropping systems and irrigation regimes. Soil biological properties, such as microbial biomass, enzyme activity, and microbial community composition, were not included in the analysis. This exclusion was mainly due to the scope and objectives of the study, which focused on soil fertility indicators directly linked to plant water availability, nutrient status, and crop productivity under deficit irrigation conditions. Despite these limitations, the selected soil physicochemical indicators provide reliable and widely accepted measures of soil fertility and environmental response to irrigation and cropping systems. However, future research should integrate soil biological, chemical, and physical indicators to provide a more comprehensive understanding of soil health and ecosystem functioning under intercropping and deficit irrigation systems.
5 Conclusions
This study demonstrated that both cropping system and irrigation level significantly influenced crop growth, biomass production, yield performance, water-use efficiency, and soil physicochemical properties. The interaction between these factors played an important role in regulating crop productivity and soil nutrient dynamics, highlighting the importance of integrating appropriate water management and cropping strategies in water-limited production systems. Although moderate deficit irrigation (75% DI) resulted in some reductions in productivity relative to full irrigation, it maintained acceptable crop growth, yield, water-use efficiency, and soil physicochemical conditions, demonstrating its potential as a water-conserving irrigation strategy.
Cropping system significantly influenced both crop performance and soil physicochemical properties. Sole dry bean consistently produced the greatest growth and yield; however, the performance of intercropping systems varied according to the companion crop. Among the intercrop combinations evaluated, the dry bean + Cucumis myriocarpus intercrop exhibited greater compatibility and tolerance to water deficit than the dry bean + Cleome gynandra intercrop, maintaining productivity levels closer to those of sole dry bean under reduced irrigation. The study further showed that dry bean-based systems generally maintained higher concentrations of exchangeable cations and greater cation exchange capacity than systems dominated by Cleome gynandra. Based on these findings, 75% DI represents a viable water-saving strategy, as it maintains relatively stable soil moisture conditions that sustain nutrient mobility and reduce excessive salinity buildup while conserving irrigation water. The dry bean + Cucumis myriocarpus intercrop appears better adapted to water-limited conditions due to stronger functional complementarity and lower interspecific competition compared with the dry bean + Cleome gynandra system.
5.1 Recommendations
Future research should the long-term impacts of deficit irrigation and intercropping on soil physicochemical and biological properties, water productivity, and economic viability across diverse agroecological zones to further evaluate the sustainability and scalability of these production systems. Investigate the long-term dynamics of salt accumulation under repeated deficit irrigation, particularly its interaction with soil texture and rainfall patterns, as well as quantify the balance between facilitation and competition in intercropping systems using detailed root trait and rhizosphere process measurements. Integrating soil biological indicators would also improve understanding of the microbial mechanisms driving nutrient cycling under these management practices.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
BL: Conceptualization, Software, Visualization, Investigation, Writing – review & editing, Writing – original draft, Methodology, Data curation, Formal analysis. PK: Funding acquisition, Conceptualization, Writing – review & editing, Resources, Project administration, Supervision, Validation, Formal analysis. PS: Formal analysis, Project administration, Supervision, Conceptualization, Resources, Visualization, Funding acquisition, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This study was funded by the National Research Foundation (NRF) and the German Academic Exchange Service (DAAD).
Acknowledgments
Many thanks are extended to Tshwane University of Technology, Central University of Technology, and the University of Limpopo for their valuable collaboration and contribution to the successful completion of this study.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fagro.2026.1867943/full#supplementary-material
References
1
AguilosM.StahlC.BurbanB.HéraultB.CourtoisE.CosteS.et al. (2019). Interannual and seasonal variations in ecosystem transpiration and water use efficiency in a tropical rainforest. Forests10 (1). doi: 10.3390/f10010014
2
AkchayaK.ParasuramanP.PandianK.VijayakumarS.ThirukumaranK.MustaffaM. R. A. F.et al. (2025). Boosting resource use efficiency, soil fertility, food security, ecosystem services, and climate resilience with legume intercropping: a review. Front. Sustain. Food Syst.9, 1527256. doi: 10.3389/fsufs.2025.1527256
3
AkinolaF. F.FasinmirinJ. T.OlanrewajuO. O.EwuloB. S.OlorunfemiI. E. (2023). Changes in soil physical and mechanical properties under different tillage and cropping systems in alfisol soil of southwestern Nigeria. Farming Syst.1 (3), 100050. doi: 10.1016/j.farsys.2023.100050
4
AlamA. (2014). Soil degradation: a challenge to sustainable agriculture. Int. J. Sci. Res. Agric. Sci.1 (4), 50–55.
5
AloriE. T.GlickB. R.BabalolaO. O. (2017). Microbial phosphorus solubilization and its potential for use in sustainable agriculture. Front. Microbiol.8, 971. doi: 10.3389/fmicb.2017.00971
6
Al-ShammaryA. A. G.KouzaniA. Z.KaynakA.KhooS. Y.NortonM.GatesW. (2018). Soil bulk density estimation methods: a review. Pedosphere28 (4), 581–596. doi: 10.1016/S1002-0160(18)60034-7
7
AnanthiT.AmanullahM. M.Al-TawahaA. R. M. S. (2017). A review on maize-legume intercropping for enhancing the productivity and soil fertility for sustainable agriculture in India. Adv. Environ. Biol.11 (5), 49–64.
8
AndrocioliL. G.ZeffaD. M.AlvesD. S.TomazJ. P.Moda-CirinoV. (2020). Effect of water deficit on morphoagronomic and physiological traits of common bean genotypes with contrasting drought tolerance. Water12 (1), 217. doi: 10.3390/w12010217
9
ArayaT.OchsnerT. E.MnkeniP. N.HounkpatinK. O. L.AmelungW. (2024). Challenges and constraints of conservation agriculture adoption in smallholder farms in sub-Saharan Africa: a review. Int. Soil Water Conserv. Res.12 (4), 828–843. doi: 10.1016/j.iswcr.2024.03.001
10
AshineE. T.BedaneM. T.MengeshaA. A. (2024). Response of carrot (Daucus carota L.) to supplementary irrigation under rain-fed agriculture at Jimma and Gera, Jimma Zone, South West Ethiopia. Adv. Agric.2024 (1), 3976619. doi: 10.1155/2024/3976619
11
BaderB. R.TabanS. K.FahmiA. H.AboodM. A.HamdiG. J. (2021). Potassium availability in soil amended with organic matter and phosphorous fertiliser under water stress during maize (Zea mays L) growth. J. Saudi Soc. Agric. Sci.20 (6), 390–394. doi: 10.1016/j.jssas.2021.04.006
12
BaiW.KongL.GuoA. (2013). Effects of physical properties on electrical conductivity of compacted lateritic soil. J. Rock Mech. Geotech. Eng.5 (5), 406–411. doi: 10.1016/j.jrmge.2013.07.003
13
BeauclaireQ.Vanden BrandeF.LongdozB. (2024). Key role played by mesophyll conductance in limiting carbon assimilation and transpiration of potato under soil water stress. Front. Plant Sci.15. doi: 10.3389/fpls.2024.1500624
14
BechemE. E.OjongA. N.EtchuK. A. (2018). The effects of intercropping and plant densities on growth and yield of maize (Zea mays L.) and soybean (Glycine max) in the humid forest zone of Mount Cameroon. Afr. J. Agric. Res.13 (12), 574–587.
15
BhattacharyaA. (2021). “ Effect of soil water deficit on growth and development of plants: a review.” in Soil Water Deficit and Physiological Issues in Plants, (Singapore: Springer Singapore (Springer Nature Singapore Pte Ltd.)), 393–488. doi: 10.1007/978-981-33-6276-5_5
16
BuckleyT. N.MagneyT. S.Berny Mier y TeranJ. C.MillsC.PalkovicA.ParkerT. A.et al. (2025). Diversity in stomatal and hydraulic responses to post-flowering drought in common (Phaseolus vulgaris) and tepary (P. acutifolius) beans. Plant Cell Environ.48 (1), 51–64. doi: 10.1111/pce.15106
17
BuelvasR. M.AdamchukV. I.LeksonoE.TikaszP.LefsrudM.HoloszkiewiczJ. (2019). Biomass estimation from canopy measurements for leafy vegetables based on ultrasonic and laser sensors. Comput. Electron. Agric.164, 104896. doi: 10.1016/j.compag.2019.104896
18
ButcherK. R.NastoM. K.NortonJ. M.StarkJ. M. (2020). Physical mechanisms for soil moisture effects on microbial carbon-use efficiency in a sandy loam soil in the western United States. Soil Biol. Biochem.150, 107969. doi: 10.1016/j.soilbio.2020.107969
19
Cárceles RodríguezB.Durán-ZuazoV. H.Soriano RodríguezM.García-TejeroI. F.Gálvez RuizB.Cuadros TaviraS. (2022). Conservation agriculture as a sustainable system for soil health: a review. Soil Syst.6 (4), 87. doi: 10.3390/soilsystems6040087
20
CescoS.ZuluagaM. Y. A.CavaniL.BorrusoL.LaudicinaA. V.MazzettoF.et al. (2026). Roots and the rhizosphere: a perspective on the hidden engine of regenerative, antifragile, and digitally enabled agriculture. Farming Syst.4, 100199. doi: 10.1016/j.farsys.2026.100199
21
ChamkhiI.ChetoS.GeistlingerJ.ZeroualY.KouisniL.BargazA.et al. (2022). Legume-based intercropping systems promote beneficial rhizobacterial community and crop yield under stressing conditions. Ind. Crops Prod.183, 114958. doi: 10.1016/j.indcrop.2022.114958
22
ChengM.WangH.FanJ.ZhangF.WangX. (2021). Effects of soil water deficit at different growth stages on maize growth, yield, and water use efficiency under alternate partial root-zone irrigation. Water13 (2), 148. doi: 10.3390/w13020148
23
CraineJ. M.DybzinskiR. (2013). Mechanisms of plant competition for nutrients, water and light. Funct. Ecol.27, 833–840. doi: 10.1111/1365-2435.12081
24
DerpschR.KassamA.ReicoskyD.FriedrichT.CalegariA.BaschG.et al (2024). Nature's laws of declining soil productivity and conservation agriculture. Soil Secur.14, 100127. doi: 10.1016/j.soisec.2024.100127
25
Di BertS.BenardP.JiaR.WankmüllerF. J.AzadS.KaestnerA.et al. (2026). Early signals of water limitations begin at the root–soil interface: linking rhizosphere drying to water uptake decline. New Phytol.249 (6), 2909–2918. doi: 10.1111/nph.70879
26
DohrmannR. (2006). Cation exchange capacity methodology I: an efficient model for the detection of incorrect cation exchange capacity and exchangeable cation results. Appl. Clay Sci.34 (1-4), 31–37. doi: 10.1016/j.clay.2005.12.006
27
DuJ.GaoQ.SunF.LiuB.JiaoY.LiuQ. (2025). Agricultural soil microbiomes at the climate frontier: nutrient-mediated adaptation strategies for sustainable farming. Ecotoxicol. Environ. Saf.295, 118161. doi: 10.1016/j.ecoenv.2025.118161
28
DucheneO.VianJ. F.CeletteF. (2017). Intercropping with legume for agroecological cropping systems: complementarity and facilitation processes and the importance of soil microorganisms. a review. Agric. Ecosyst. Environ.240, 148–161. doi: 10.1016/j.agee.2017.02.019
29
EhrmannJ.RitzK. (2014). Plant: soil interactions in temperate multi-cropping production systems. Plant Soil376 (1-2), 1–29. doi: 10.1007/s11104-013-1921-8
30
Ellawala KankanamgeC.KodithuwakkuH. (2017). Effect of interspecific competition on the growth and nutrient uptake of three macrophytes in nutrient-rich water. Aquat. Ecol.51. doi: 10.1007/s10452-017-9640-5
31
El SabaghA.IslamM. S.SkalickyM.Ali RazaM.SinghK.Anwar HossainM.et al. (2021). Salinity stress in wheat (Triticum aestivum L.) in the changing climate: adaptation and management strategies. Front. Agron.3, 661932. doi: 10.3389/fagro.2021.661932
32
EngbersenN.BrookerR. W.StefanL.StuderB.SchöbC. (2021). Temporal differentiation of resource capture and biomass accumulation as a driver of yield increase in intercropping. Front. Plant Sci.12. doi: 10.3389/fpls.2021.668803
33
Eng GiapS. G.AhmadM. F. (2025). Factors affecting soil bulk density: a conceptual model. J. Soil Environ. Agroecology1 (1), 27–45.
34
FadlM. E.SayedY. A.El-DesokyA. I.ShamsE. M.ZekariM.AbdelsamieE. A.et al (2024). Irrigation practices and their effects on soil quality and soil characteristics in arid lands: a comprehensive geomatic analysis. Soil Syst.8 (2), 52. doi: 10.3390/soilsystems8020052
35
FanH.MiaoR.GuoC.BaoX.HeW.SunY.et al. (2025). Research on the effect of diversified cropping on crop quality: a review. Agriculture15 (5), 456. doi: 10.3390/agriculture15050456
36
FlexasJ.Díaz-EspejoA.ConesaM. A.CoopmanR. E.DoutheC.GagoJ.et al (2016). Mesophyll conductance to CO2 and Rubisco as targets for improving intrinsic water use efficiency in C3 plants. Plant Cell Environ.39 (5), 965–982. doi: 10.1111/pce.12622
37
Franco-NavarroJ. D.PadillaY. G.álvarezS.Calatayudá.Colmenero-FloresJ. M.Gómez-BellotM. J.et al. (2025). Advancements in water-saving strategies and crop adaptation to drought: a comprehensive review. Physiol. Plant.177 (4), e70332. doi: 10.1111/ppl.70332
38
GaoY.LiangY.FuY.SiZ.HamaniA. K. M. (2022). Interactive effects of intraspecific competition and drought on stomatal conductance and hormone concentrations in different tomato genotypes. Horticulturae8 (1), 45. doi: 10.3390/horticulturae8010045
39
GoetG.SonkarI.KumarS.Hari PrasadK. S.OjhaC. S. P. (2024). Effect of salinity on crop growth and soil moisture dynamics: a study with root water uptake model. J. Hazard. Toxic Radioact. Waste28 (3), 04024009. doi: 10.1061/JHTRBP.HZENG-1299
40
GowdaB.DeepakT. N.SumaT. C.YadahalliG. S.MacchaS. I.PrashanthS. M. (2020). Effect of intercropping and fertilizer dose on seed quality and physiological properties of chickpea (Cicer arietinum L.). Int. J. Curr. Microbiol. Appl. Sci.9. doi: 10.20546/ijcmas.2020.901.064
41
GuerchiA.MnafguiW.GarijoG.MahjoubA.LudidiN.GonzalezE. M.et al. (2025). Impact of intercropping on agronomic and metabolic responses of Medicago sativa and Hordeum marinum under nutrient deficiency and drought stress. Front. Agron.7, 1636363. doi: 10.3389/fagro.2025.1636363
42
HeJ.HeJ.LiH.YuY.QianL.TangL.et al. (2024). Continuous intercropping increases the depletion of soil available and non-labile phosphorus. Agronomy14 (6), 1121. doi: 10.3390/agronomy14061121
43
HeuscherS. A.BrandtC. C.JardineP. M. (2005). Using soil physical and chemical properties to estimate bulk density. Soil Sci. Soc. Am. J.69. doi: 10.2136/sssaj2005.0051a
44
HolzM.ZarebanadkoukiM.BenardP.HoffmannM.DubbertM. (2024). Root and rhizosphere traits for enhanced water and nutrients uptake efficiency in dynamic environments. Front. Plant Sci.15, 1383373. doi: 10.3389/fpls.2024.1383373
45
HuangY.QiJ.MezbahuddinS.DyckM.MacKenzieD.FaramarziM. (2026). Effects of soil moisture and winter hydrologic processes on soil phosphorous accumulation and loss in canola croplands of cold-region watersheds. Sci. Total Environ.1010. doi: 10.1016/j.scitotenv.2025.181051
46
HussenA.WorkuW.ZewdieM. (2019). Effects of deficit irrigation and phosphorus levels on growth, yield, yield components and water use efficiency of mung bean (Vigna radiata (L.) Wilczek) at Alage, Central Rift Valley of Ethiopia. Agric. Res. Technol. Open Access J.21 (3). doi: 10.19080/ARTOAJ.2019.21.556170
47
Intergovernmental Panel on Climate Change (IPCC) (2023). Climate Change 2022 – Impacts, Adaptation and Vulnerability, (Cambridge: Cambridge University Press). doi: 10.1017/9781009325844
48
JiangP.WangY.ZhangY.FeiJ.RongX.PengJ.et al (2024). Intercropping enhances maize growth and nutrient uptake by driving the link between rhizosphere metabolites and microbiomes. New Phytol.243 (4), 1506–1521. doi: 10.1111/nph.19906
49
KabatoW.HailegnawN.MutumL.MolnarZ. (2025). Managing soil health for climate resilience and crop productivity in a changing environment. Sci. Total Environ.1000, 180460. doi: 10.1016/j.scitotenv.2025.180460
50
KarkanisP. G.AuK.SchaaljeG. B. (1991). Comparison of four measurement schedules for determination of soil particle-size distribution by the hydrometer method. Can. Agric. Eng.33 (2), 211–216.
51
KatariaP.SinghJ.SinghG.MaviM. S.SillanpääM.Al-FarrajS. (2024). Residual phosphorus availability after soil application of different organic waste in varying soil P status soils. Heliyon10. doi: 10.1016/j.heliyon.2024.e25732
52
KebedeH.AbbasH. K.FisherD. K.BellalouiN. (2012). Relationship between aflatoxin contamination and physiological responses of corn plants under drought and heat stress. Toxins4 (11), 1385–1403.
53
KgopaP. M.MashelaP. W.ManyevereA. (2018). Suitability of treated wastewater with respect to pH, electrical conductivity, selected cations and sodium adsorption ratio for irrigation in a semi-arid region. Water SA44 (4), 551–556. doi: 10.4314/wsa.v44i4.04
54
KimH. N.ParkJ. H. (2024). Monitoring of soil EC for the prediction of soil nutrient regime under different soil water and organic matter contents. Appl. Biol. Chem.67 (1), 1. doi: 10.1186/s13765-023-00797-2
55
KimM.RabbaniM. G.ZebroM.ChoiK. Y. (2025). Effect of nutrient solution electrical conductivity on cucumber growth and yield in controlled pot soil cultivation. Notulae Botanicae Horti Agrobotanici Cluj-Napoca53 (3), 14548. doi: 10.15835/nbha53314548
56
KumarS.SonkarI.GuptaV.Hari PrasadK. S.OjhaC. S. P. (2021). Effect of salinity on moisture flow and root water uptake in sandy loam soil. J. Hazard. Toxic Radioact. Waste25 (3). doi: 10.1061/(ASCE)HZ.2153-5515.0000582
57
KumawatA.BamboriyaS. D.MeenaR. S.YadavD.KumarA.KumarS.et al. (2022). “ Legume-based inter-cropping to achieve the crop, soil, and environmental health security.” in Advances in Legumes for Sustainable Intensification, (Cambridge, Massachusetts: Academic Press), 307–328. doi: 10.1016/B978-0-323-85797-0.00005-7
58
KwartengA. O.AbogoomJ.AmoahR. A.NyadanuD.NyamC. K.GhunneyT.et al. (2018). Phenomic characterization of twenty-four accessions of spider plant (Cleome gynandra L.) the Upper East Region of Ghana. Sci. Hortic.235, 124–131. doi: 10.1016/j.scienta.2018.02.046
59
LadhaJ. K.PeoplesM. B.ReddyP. M.BiswasJ. C.BennettA.JatM. L.et al. (2022). Biological nitrogen fixation and prospects for ecological intensification in cereal-based cropping systems. Field Crops Res.283, 108541. doi: 10.1016/j.fcr.2022.108541
60
LambersH.OliveiraR. S. (2019). “ Photosynthesis, respiration, and long-distance transport: photosynthesis.” in Plant Physiological Ecology, (Cham: Springer International Publishing AG), 11–114. doi: 10.1007/978-3-030-29639-1_2
61
LehmannJ.BossioD. A.Kögel-KnabnerI.RilligM. C. (2020). The concept and future prospects of soil health. Nat. Rev. Earth Environ.1 (10), 544–553. doi: 10.1038/s43017-020-0080-8
62
LiuH.GaoX.LiC.CaiY.SongX.ZhaoX. (2025). Intercropping increases plant water availability and water use efficiency: A synthesis. Agric. Ecosyst. Environ.379, 109360. doi: 10.1016/j.agee.2024.109360
63
LuD.LiuL.BaiY.AnQ.ChengY.HuangG. (2025). Deficit irrigation alleviates the increase in soil salinity content in saline-alkali regions of China and improves irrigation water productivity: A meta-analysis. Agric. Water Manage.320, 109872. doi: 10.1016/j.agwat.2025.109872
64
LuoC.WangR.DouX.ZhengC. (2023). Soil moisture, nutrients, and plant growths under various irrigation and fertilization regimes during the crop replacement period in an alley intercropping system on the Loess Plateau of China. Forests14 (11), 2153. doi: 10.3390/f14112153
65
MaQ.ChenH.YangY.ZhouB. (2025). Unlocking phosphorus resources: Phosphate-solubilizing microorganisms as a green strategy for activating phosphorus in acidic red soils and promoting crop growth. Front. Microbiol.16. doi: 10.3389/fmicb.2025.1630650. Preprint.
66
MahmoodM.WangJ.MehmoodS.AhmedW.AyyoubA.SeleimanM. F.et al. (2025). Influence of drought stress on phosphorus dynamics and maize growth in tropical ecosystems. BMC Plant Biol.25 (1), 62. doi: 10.1186/s12870-025-06092-x
67
Maltais-LandryG. (2015). Legumes have a greater effect on rhizosphere properties (pH, organic acids and enzyme activity) but a smaller impact on soil P compared to other cover crops. Plant Soil394, 139–154. doi: 10.1007/s11104-015-2518-1
68
MarifA. A.EsmailA. O. (2025). Assessing the impact of seasonal variability on irrigation water quality and suitability for agricultural use in wet and dry conditions. Acta Agric. Slov.121 (4), 1–13. doi: 10.14720/aas.2025.121.4.21860
69
MasekoI.NcubeB.TesfayS.FessehazionM.ModiA. T.MabhaudhiT. (2020). Productivity of selected African leafy vegetables under varying water regimes. Agronomy10 (6), 916. doi: 10.3390/agronomy10060916
70
MashamaiteC. V.NgcoboB. L.ManyevereA.BertlingI.FawoleO. A. (2022). Assessing the usefulness of Moringa oleifera leaf extract as a biostimulant to supplement synthetic fertilizers: a review. Plants11 (17), 2214. doi: 10.3390/plants11172214
71
MashelaP. W.PofuK. M. (2012). Interactive effects of Meloidogyne incognita race 2, Bradyrhizobium japonicum and crude extracts of Cucumis myriocarpus fruit on Vigna unguiculata. Crop Prot.42, 124–127. doi: 10.1016/j.cropro.2012.06.013
72
MasondoN. A.AdetunjiT. L.MofokengM. M.Naidoo-MaharajD.KoetleM. J.OlisahC.et al (2025). Medicinal plant value chain development in South Africa – A semicentennial appraisal and research trends. S. Afr. J. Bot.186, 88–110. doi: 10.1016/j.sajb.2025.08.018
73
MateuL. A. B.ValladaresD. L.RodríguezM. B.RojasS. F.GallardoL. E.PalmaM. A. M.et al. (2019). Obesity status of children from 4 to 7 years of age in the commune of Arica: A comparative study according to gender. Stud. Ethno-Med.13 (1), 1–7.
74
MavimbelaS. S.van RensburgL. D. (2019). Estimating soil hydraulic parameters characterizing rainwater infiltration and runoff properties of dryland floodplains. Comput. Water Energy Environ. Eng.8, 11–40. doi: 10.4236/cweee.2019.81002
75
MegerssaG. R.BekereY. B. (2019). Causes, consequences and coping strategies of land degradation: evidence from Ethiopia. J. Degrad. Min. Lands Manage.7 (1), 1953. doi: 10.15243/jdmlm.2019.071.1953
76
MjokweS.ElephantD. E.ManyevereA.MashamaiteC. V. (2025). Unlocking soil phosphorus resources in semi-arid regions using arbuscular mycorrhizal fungi and phosphorus solubilising bacteria: a systematic review. Discover Soil2 (1), 105. doi: 10.1007/s44378-025-00137-1
77
MokgoloM. J.MzezewaJ.ZerizghyM. G. (2024). Investigating the synergistic effect of tillage system and manure application rates on selected properties of two soil types in Limpopo Province, South Africa. Sustainability16 (20), 8941. doi: 10.3390/su16208941
78
MorrisM. L. (2007). Fertilizer Use in African Agriculture: Lessons Learned and Good Practice Guidelines (Washington, DC: World Bank). doi: 10.1596/978-0-8213-6880-0
79
MudauF. N.ChimonyoV. G. P.ModiA. T.MabhaudhiT. (2022). Neglected and underutilised crops: a systematic review of their potential as food and herbal medicinal crops in South Africa. Front. Pharmacol.12, 809866. doi: 10.3389/fphar.2021.809866
80
NavegantesP. C. A.Lopes FilhoW. R. L.RodriguesF. H. S.MonteiroG. G. T. N.ChavesR. P. F.Oliveira NetoC. F.et al. (2024). Trocas gasosas e relações hídricas em duas cultivares de açaí submetidas ao déficit hídrico. Rev. Bras. Frutic.46, e716. doi: 10.1590/0100-29452024716
81
NguyenT. B. A.LefoulonC.NguyenT. H.BlattM. R.CarrollW. (2023). Engineering stomata for enhanced carbon capture and water-use efficiency. Trends Plant Sci.28 (11), 1290–1309. doi: 10.1016/j.tplants.2023.06.002
82
NkgapeleR. J.MphosiM. S.MashelaP. W. (2011). Comparing Cucumis africanus and Cucumis myriocarpus using a classical growth analysis tool under five irrigation frequencies. WIT Trans. Ecol. Environ.145, 565–570. doi: 10.2495/WRM110501
83
NongS.ChenH.ChenZ.LinZ.SuS.LeiX.et al. (2025). Effects of seasonal rainfall changes on N, P, and K stoichiometric characteristics in leaves and soil of tropical coastal shelterbelt forests. Forests16 (4), 600. doi: 10.3390/f16040600
84
NwaforI.NwaforC.MandunaI. (2021). Constraints to cultivation of medicinal plants by smallholder farmers in South Africa. Horticulturae7 (12), 531. doi: 10.3390/horticulturae7120531
85
OgorekL. L. P.GaoY.FarrarE.PandeyB. K. (2025). Soil compaction sensing mechanisms and root responses. Trends Plant Sci.30 (5), 565–575. doi: 10.1016/j.tplants.2024.10.014
86
OliveiraM.RebacD.CoutinhoJ.FerreiraL. M. M.TrindadeH. (2020). Nitrogen mineralization of legume residues: interactions between species, temperature and placement in soil. Span. J. Agric. Res.18 (1), 1101. doi: 10.5424/sjar/2020181-15174
87
OliveiraR. H.RosolemC. A.TrigueiroR. M. (2004). “ Importance of mass flow and diffusion on the potassium supply to cotton plants as affected by soil water and potassium,” in Revista Brasileira De Ciencia do Solo, vol. 28. . Preprint. doi: 10.1590/s0100-06832004000300005
88
Oliveira NetoS. S.PereiraF. F. S.ZozT.OliveiraC. E. S.Moda CirinoV. (2023). Effect of water deficit on morphoagronomic traits of black common bean genotypes (Phaseolus vulgaris L.) with contrasting drought tolerance. J. Agron. Crop Sci.209 (1), 83–93. doi: 10.1111/jac.12594
89
OmutoC. T.KomeG. K.RamakhannaS. J.MuziraN. M.RuleyJ. A.JayeobaO. J.et al. (2024). Trend of soil salinization in Africa and implications for agro-chemical use in semi-arid croplands. Sci. Total Environ.951, 175503. doi: 10.1016/j.scitotenv.2024.175503
90
PabuayonI. L. B.BicaldoJ. J. B.AlemarZ. A.PabuayonI. M.RitchieG. L. (2025). Dynamics of source-sink relationships in crops under marginal environments. Curr. Plant Biol.44, 100563. doi: 10.1016/j.cpb.2025.100563
91
PanH.SongQ.YiL.LiB.BiR. (2025). Effects of different irrigation amounts on soil structure in newly cultivated land on the Loess Plateau. PLoS One20 (5), e0323065. doi: 10.1371/journal.pone.0323065
92
PanagosP.De RosaD.LiakosL.LabouyrieM.BorrelliP.BallabioC. (2024). Soil bulk density assessment in Europe. Agric. Ecosyst. Environ.364, 108907. doi: 10.1016/j.agee.2024.108907
93
PereiraL. S.MotaM.RazieiT.ParedesP. (2024). Water requirements and crop coefficients of edible, spicy and medicinal herbs and vegetables; a review aimed at supporting plant and water management. Irrig. Sci.42 (6), 1199–1228. doi: 10.1007/s00271-024-00960-4
94
PiccinniG.KoJ.MarekT.LeskovarD. I. (2009). Crop coefficients specific to multiple phenological stages for evapotranspiration-based irrigation management of onion and spinach. HortScience44 (2), 421–425. doi: 10.21273/HORTSCI.44.2.421
95
PozzaL. E.FieldD. J. (2020). The science of soil security and food security. Soil Secur.1, 100002. doi: 10.1016/j.soisec.2020.100002
96
QadirM.SchubertS.OsterJ. D.SpositoG.MinhasP. S.CheraghiS. A.et al. (2018). High-magnesium waters and soils: emerging environmental and food security constraints. Sci. Total Environ.642, 1108–1117. doi: 10.1016/j.scitotenv.2018.06.090
97
QiaoM.HongC.JiaoY.HouS.GaoH. (2024). Impacts of drought on photosynthesis in major food crops and the related mechanisms of plant responses to drought. Plants13 (13), 1808. doi: 10.3390/plants13131808
98
RamatsitsiM. N.DubeZ. P. (2020). Post-infectional resistance in traditional leafy vegetable infected with root-knot nematodes. S. Afr. J. Bot.131, 169–173. doi: 10.1016/j.sajb.2020.01.023
99
RazaM. A.DinA. M. U.YasinH. S.GulH.ShahG. A.ZhiqiW.et al (2025). Maize/legume intercropping increases nutrient uptake, crop yields, land productivity, and economic profits in resource-intensive arid-irrigated areas. J. Soil Sci. Plant Nutr.25 (3), 8316–8332. doi: 10.1007/s42729-025-02675-7
100
RezaeiM.BazarganK.ShahbaziK.MarziM.CheraghiM. (2025). Modelling phosphorus and potassium dynamics in drip-irrigated potato systems using coupled agro-hydrological model. Agric. Water Manage.321, 109920. doi: 10.1016/j.agwat.2025.109920
101
RhoadesJ. D.ManteghiN. A.ShouseP. J.AlvesW. J. (1989). Soil electrical conductivity and soil salinity: New formulations and calibrations. Soil Sci. Soc. Am. J.53 (2), 433–439. doi: 10.2136/sssaj1989.03615995005300020020x
102
SadequeA.AhmedS.IslamM. A.MuktadirM. A.RahmanM. Z.IslamM. S.et al. (2026). Enhancing water use efficiency in grain legumes: from agronomic practices to breeding strategies. J. Plant Nutr.49 (2), 203–245. doi: 10.1080/01904167.2025.2544701
103
SahuG.RoutP. P.MohapatraS.DasS. P.PradhanP. P. (2020). Climate smart agriculture: a new approach for sustainable intensification. Curr. J. Appl. Sci. Technol.39 (23), 138–147. doi: 10.9734/cjast/2020/v39i2330862
104
SainjuU. M.LiptzinD.DangiS. M. (2022). Enzyme activities as soil health indicators in relation to soil characteristics and crop production. Agrosystems Geosci. Environ.5, 20297. doi: 10.1002/agg2.20297
105
SaleemM.PervaizZ. H.ContrerasJ.LindenbergerJ. H.HuppB. M.ChenD.et al. (2020). Cover crop diversity improves multiple soil properties via altering root architectural traits. Rhizosphere16, 100248. doi: 10.1016/j.rhisph.2020.100248
106
SanandajiD. S.HeidariG.FathiP.KhodaverdilooH.SharifiZ. (2024). Soil properties and quinoa yield as affected by irrigation and vermicompost application under a drip-tape irrigation system. Front. Sustain. Food Syst.8. doi: 10.3389/fsufs.2024.1277934
107
SanatuA. M.AbubakariA. (2022). Assessment of cowpea (Vigna unguiculata (L.) Walp) F1 lines response to drought tolerance. Ghana J. Sci. Technol. Dev.8 (1), 48–60.
108
SardansJ.PeñuelasJ. (2021). Potassium control of plant functions: Ecological and agricultural implications. Plants10 (2), 419. doi: 10.3390/plants10020419
109
SemchenkoM.XueP.LeighT. (2021). Functional diversity and identity of plant genotypes regulate rhizodeposition and soil microbial activity. New Phytol.232 (2), 776–787. doi: 10.1111/nph.17604
110
ShahA.SultanH.KhanM. N.AliW.AbbasH. M. M.KhanM. N.et al. (2026). Iron and zinc-modified nano-biochar alleviate salinity stress in paddy soil by modulating nutrient availability and microbial communities. Front. Microbiol.16, 1729271. doi: 10.3389/fmicb.2025.1729271
111
ShaikR. S.BurrowsG. E.UrwinN. A.GopurenkoD.LepschiB. J.WestonL. A. (2017). The biology and management of prickly paddy melon (Cucumis myriocarpus L.), an important summer annual weed in Australia. Crop Prot.92, 29–40. doi: 10.1016/j.cropro.2016.10.005
112
ShenX.LiuY.LiX.LiL. (2024). Effect of intercropping with legumes at different rates on the yield and soil physicochemical properties of Cyperus esculentus L. in arid land. Front. Plant Sci.15. doi: 10.3389/fpls.2024.1351843
113
SilbergT. R.ChimonyoV. G. P.RichardsonR. B.SnappS. S.RennerK. (2019). Legume diversification and weed management in African cereal-based systems. Agric. Syst.174, 83–94. doi: 10.1016/j.agsy.2019.05.004
114
SimonN. S.MmatekoK. P.OchandaO. J. B. (2024). Maize-legume intercropping system improved soil quality under irrigated and rainfed conditions. Res. Crops25 (4), 553–559.
115
SinghN. R.SinghA.DeviN. P.KumarY. B.SangmaR. H. C.PhilanimW. S.et al. (2024). “ Agroforestry for soil health.” in Agroforestry, (Cham: Springer Nature), 255–283. doi: 10.1002/9781394231164.ch9
116
SitholeM. Z.MorepjeM. T.MokoenaT. I. (2025). Unpacking water scarcity adaptation strategies for sustainable food production systems in sub-Saharan Africa. Sustainability17 (23), 10627. doi: 10.3390/su172310627
117
SongY.SunJ.LiZ. (2025). Intercropping increases land use efficiency and reduces soil salinity in salt-affected soils, but has minimal effect on pH. Field Crops Res.332, 110028. doi: 10.1016/j.fcr.2025.110028
118
SoniP. G.BasakN.RaiA. K.SundhaP.NarjaryB.KumarP.et al. (2021). Deficit saline water irrigation under reduced tillage and residue mulch improves soil health in sorghum-wheat cropping system in semi-arid region. Sci. Rep.11 (1), 1880. doi: 10.1038/s41598-020-80364-4
119
SuriyagodaL. D. B.RyanM. H.RentonM.LambersH. (2014). Plant responses to limited moisture and phosphorus availability: a meta-analysis. Adv. Agron.124, 143–200. doi: 10.1016/B978-0-12-800138-7.00004-8
120
TanB.LiY.LiuT.KongX.WangY.XiaoJ.et al. (2021). Response of plant rhizosphere microenvironment to water management in soil- and substrate-based controlled environment agriculture (CEA) systems: a review. Front. Plant Sci.12, 691651. doi: 10.3389/fpls.2021.691651
121
TanW.SunL.ZhangX.SongC. (2020). The effects of water levels and interspecific competition on two Carex species in a temperate wetland of northeast China. Sustainability12 (24), 10654. doi: 10.3390/su122410654
122
TangQ.CottonA.WeiZ.XiaY.DaniellT.YanX. (2022). How does partial substitution of chemical fertiliser with organic forms increase sustainability of agricultural production?. Sci. Total Environ.803, 149933. doi: 10.1016/j.scitotenv.2021.149933
123
TeferaM. L.SeddaiuG.CarlettiA.AwadaH. (2025). Rainfall variability and drought in West Africa: challenges and implications for rainfed agriculture. Theor. Appl. Climatol.156 (1), 41. doi: 10.1007/s00704-024-05278-1
124
TokerP.CanciH.TurhanI.IsciA.ScherzingerM.KordrostamiM.et al (2024). The advantages of intercropping to improve productivity in food and forage production–a review. Plant Prod. Sci.27 (3), 155–169. doi: 10.1080/1343943X.2024.2372878
125
ŢopaD. C.Cap?unaS.CalistruA. E.AilincaiC. (2025). Sustainable practices for enhancing soil health and crop quality in modern agriculture: a review. Agriculture15 (9), 998. doi: 10.3390/agriculture15090998
126
TouchV.TanD. K. Y.CookB. R.Li LiuD.CrossR.TranT. A.et al. (2024). Smallholder farmers' challenges and opportunities: Implications for agricultural production, environment and food security. J. Environ. Manage.370, 122536. doi: 10.1016/j.jenvman.2024.122536
127
ToumiI.GhrabM.ZarroukO.NagazK. (2024). Impact of deficit irrigation strategies using saline water on soil and peach tree yield in an arid region of Tunisia. Agriculture14 (3), 377. doi: 10.3390/agriculture14030377
128
TripathiS.SrivastavaP.DeviR. S.BhadouriaR. (2020). “ Influence of synthetic fertilizers and pesticides on soil health and soil microbiology” in Agrochemicals Detection, Treatment and Remediation, (Winnipeg, Manitoba: Butterworth-Heinemann), 25–54.
129
TullyK.SullivanC.WeilR.SanchezP. (2015). The state of soil degradation in Sub-Saharan Africa: baselines, trajectories, and solutions. Sustainability7 (6), 6523–6552. doi: 10.3390/su7066523
130
UsmanM.FarooqM.WakeelA.NawazA.CheemaS. A.ur RehmanH.et al. (2020). Nanotechnology in agriculture: current status, challenges and future opportunities. Sci. Total Environ.721, 137778. doi: 10.1016/j.scitotenv.2020.137778
131
UthaiwatP.LeeudomwongT.SanponputeT. (2024). The comparison of flow rate calibration methods for high-performance liquid chromatography (HPLC) pump. Accredit. Qual. Assur.29 (3), 205–214. doi: 10.1007/s00769-024-01580-3
132
VenuN.DeviN. S.DasH.DeviL. N.AnushaM. (2023). Inference of rock phosphate, phosphorus solubilizing bacteria and lime on phosphorus content and economic yield of green gram. Int. J. Plant Soil Sci.35 (16), 185–195. doi: 10.9734/ijpss/2023/v35i163144
133
VermaV. K.ChandraB.KumarS. (2019). Values and ascribed responsibility to predict consumers’ attitude and concern towards green hotel visit intention. J. Business Res.96. doi: 10.1016/j.jbusres.2018.11.021
134
WangY.ChenJ.SunY.JiaoY.YangY.YuanX.et al. (2023a). Zeolite reduces N leaching and runoff loss while increasing rice yields under alternate wetting and drying irrigation regime. Agric. Water Manage.277, 108130. doi: 10.1016/j.agwat.2022.108130
135
WangY.LuoD.XiongZ.WangZ.GaoM. (2023b). Changes in rhizosphere phosphorus fractions and phosphate-mineralizing microbial populations in acid soil as influenced by organic acid exudation. Soil Tillage Res.225, 105543. doi: 10.1016/j.still.2022.105543
136
WangZ.SlotM.WangC. (2026). Decoupling of stomatal conductance, transpiration and photosynthesis in terrestrial plants under elevated temperature: a meta-analysis. Nat. Commun.17, 1528. doi: 10.1038/s41467-025-68250-x
137
WangM.WangR.SunQ.LiY.XuL.WangY. (2024). Integrated drip irrigation regulates soil water–salt movement to improve water use efficiency and maize yield in saline–alkali soil. Water (Switzerland)16. doi: 10.3390/w16172509
138
WenZ.WuJ.YangY.LiR.OuyangZ.ZhengH. (2022). Implementing intercropping maintains soil water balance while enhancing multiple ecosystem services. Catena217. doi: 10.1016/j.catena.2022.106426
139
XiongD.NadalM. (2020). Linking water relations and hydraulics with photosynthesis. Plant J.101 (4), 800–815. doi: 10.1111/tpj.14595
140
XuQ.XiongK.ChiY.SongS. (2021). Effects of crop and grass intercropping on the soil environment in the karst area. Sustainability13 (10), 5484. doi: 10.3390/su13105484
141
Yan-JunW.Gao-XiangQ.Na-NaW.Hong-YunD.ZhangY.LuH.et al (2025). The improvement effects of intercropping systems on saline-alkali soils and their impact on microbial communities. Microorganisms13 (7), 1436. doi: 10.3390/microorganisms13071436
142
YangY. J.BiM. H.NieZ. F.JiangH.LiuX. D.FangX. W.et al. (2021). Evolution of stomatal closure to optimize water-use efficiency in response to dehydration in ferns and seed plants. New Phytol.230 (5), 2001–2010. doi: 10.1111/nph.17278
143
YangT.SiddiqueK. H. M.LiuK. (2020). “ Cropping systems in agriculture and their impact on soil health-a review,” in Global Ecology and Conservation, vol. 23. Preprint. doi: 10.1016/j.gecco.2020.e01118
144
YangH.ZhangW.LiL. (2021). Intercropping: Feed more people and build more sustainable agroecosystems. Front. Agric. Sci. Eng.8. doi: 10.15302/J-FASE-2021398
145
YangM.ZhouD.HangH.ChenS.LiuH.SuJ.et al. (2024). Effects of balancing exchangeable cations Ca, Mg, and K on the growth of tomato seedlings (Solanum lycopersicum L.) based on increased soil cation exchange capacity. Agronomy14 (3), 629. doi: 10.3390/agronomy14030629
146
YavasI.JamalM. A.DinK. U.AliS.HussainS.FarooqM. (2024). Drought-induced changes in leaf morphology and anatomy: overview, implications and perspectives. Pol. J. Environ. Stud.33 (2), 1517–1530. doi: 10.15244/pjoes/174476
147
ZahediS. M.KarimiM.VendittiA.ZahraN.SiddiqueK. H.FarooqM. (2025). Plant adaptation to drought stress: the role of anatomical and morphological characteristics in maintaining the water status. J. Soil Sci. Plant Nutr.25 (1), 409–427. doi: 10.1007/s42729-024-02141-w
148
Zamora-LedezmaE.ProM. M.CastroE. J.VélezJ. V.SaltosR. B.VelásquezJ. V.et al. (2025). Advancing carbon quantification: a comparative evaluation of gravimetric and volumetric methods for soil carbon assessment in tropical ecosystems. Results Eng.25, 104141. doi: 10.1016/j.rineng.2024.104141
149
ZhangH.ShiL.LuH.ShaoY.LiuS.FuS. (2020). Drought promotes soil phosphorus transformation and reduces phosphorus bioavailability in a temperate forest. Sci. Total Environ.732, 139295. doi: 10.1016/j.scitotenv.2020.139295
150
ZhuS. G.ChengZ. G.BatoolA.WangY. B.WangJ.ZhouR.et al (2022). Plant facilitation shifts along with soil moisture and phosphorus gradients via rhizosphere interaction in the maize-grass pea intercropping system. Ecol. Indic.139, 108901. doi: 10.1016/j.ecolind.2022.108901
Summary
Keywords
indigenous leafy vegetables, intercropping systems, phaseolus vulgaris, regulated deficit irrigation, soil moisture content, soil physicochemical properties
Citation
Lekgoathi BP, Kgopa PM and Soundy P (2026) Intercropping dry bean (Phaseolus vulgaris L.) with Cucumis myriocarpus and Cleome gynandra under deficit irrigation: impacts on soil physicochemical properties, water-use efficiency and crop growth. Front. Agron. 8:1867943. doi: 10.3389/fagro.2026.1867943
Received
28 April 2026
Revised
17 July 2026
Accepted
20 July 2026
Published
13 August 2026
Volume
8 - 2026
Edited by
Naser A. Anjum, Aligarh Muslim University, India
Reviewed by
Haider Sultan, Jianghan University, China
Arcângelo Loss, Federal University of Santa Catarina, Brazil
Updates
Copyright
© 2026 Lekgoathi, Kgopa and Soundy.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Boitumelo Patience Lekgoathi, blekgoathi@cut.ac.za; Pholosho Mmateko Kgopa, pholosho.kgopa@ul.ac.za
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.