Carob-Thyme Intercropping Systems Can Improve Yield Efficiency and Environmental Footprint Compared to Conservation Tillage

1. Introduction

As climate change related aridity is increasing worldwide [1], drylands are predicted to cover half of the global land surface by the end of century [2]. Given that 41 percent of the Earth’s surface are already drylands, protecting these soils from desertification is critical for ecosystems and human wellbeing [3]. Agroecosystems in Mediterranean drylands are considered particularly vulnerable to drought [4]. Average temperatures in the region will likely increase by 20 percent, precipitation is expected to decline between 4-22 percent, while droughts are expected to become 5-10 times more frequent as a result of climate change [5]. Such conditions can lead to an up to 18 percent increase of irrigation requirements [4], while negatively affecting the quality and yield of agricultural crops [5]. However, soil degradation can also result from intensive agricultural practices such as tillage, monocultures and application of synthetic pesticides and fertilizers. Such unsustainable practices can lead to soil degradation, soil erosion, impairment of water infiltration and loss of biodiversity among others [4,6,7,8].

On the other hand, regenerative agricultural practices focusing on soil health, such as application of intercropping systems and reducing tillage to zero, are promising nature-based farming solutions [9] shown to increase carbon sequestration, enhance biodiversity, facilitate nutrient cycling, prevent soil erosion and increase water availability [10,11,12,13]. Such nature-based solutions can, thus, improve agricultural production and resilience, address climate change mitigation and protect biodiversity [14,15].

Intercropping, the practice of concurrently growing two or more crop species on the same agricultural plot [16] is a nature-based farming solution gaining increased attention. Intercropped agroecosystems include annual and perennial cover crops (grains, grasses, leguminous plants) or rows of perennial herbaceous plants planted amongst productive trees. Diversification through intercropping, has been found to be one of the most sustainable agricultural practices [17] providing farmers with diverse income streams [16,18]. This method can be an alternative to soil tillage or herbicides in combating weeds [17,19,20], while improving biodiversity [21,22] and promoting overall ecosystem sustainability. Comparing intercropping and monoculture systems, higher yields can be achieved for the first under conditions of reduced agricultural inputs [23,24]. Intercropping systems can obtain the same crop yields on smaller plots of land compared to monoculture, leading to a more efficient use of land [25].

Annual or perennial living plants intercropped with a main crop, can provide permanent soil cover and improve its properties. This can be achieved by increasing water infiltration, by facilitating nutrient cycling, by reducing the risk of erosion, and by reducing pesticides, herbicides and machinery use [6,11,26]. Intercropping with nitrogen-fixing leguminous plants, such as Trifolium sp., can increase soil fertility and enhance yields for the main crop [6,27]. Furthermore, establishing intercropping systems with drought-adapted perennial herbs can be an efficient way to mitigate climate change and aridification in agroecosystems. Incorporating aromatic plants in tree-based agroforestry, can achieve better pest and weed control [17], improve soil properties and water content [28,29,30] and increase yields [31]. Intercropping systems of fruit trees and aromatic plants, can either maintain [30] or significantly increase the main crop’s yields [17,32]. Thyme alley-cropped with rainfed almond trees, significantly increased soil organic carbon and moisture with no negative effect on the main crop’s yield [30]. Overall, aromatic-based intercropping systems can diversify income and provide economic returns through increased land use efficiency [17,33].

However, intercropping systems composed of fruit tree and aromatic plants have very sparsely been studied in the Mediterranean, with most studies conducted coming from China, India and Iran [17]. Moreover, in drier climates, it has been reported that intercropping methods can increase water and nutrient competition with the main crop [31,34] leading to negative impacts on the main crop’s productivity [28]. One solution could be to strategically plant living mulch systems on the south (sun) exposed soil side (SES) of trees, to reduce via shading soil temperature/evaporation without fueling completion for water and nutrients with the fruit tree. Niche differentiation between intercropped species, such as taking advantage of shading effects between species, can lead to reduced evaporation and increased soil moisture resulting to higher overall yields for the intercropped system [17].

More research is needed on dryland intercropping systems to identify under which conditions they can benefit agroecosystems and farmland productivity. This is particularly relevant when compared to conservation tillage options (1-2 tillage applications per year at the end of the wet season) that appear more attractive to farmers for combating weeds, but are still energy demanding and do not promote soil regeneration [35,36,37].

Carob (Ceratonia siliqua L.) is an evergreen species that belongs to the Fabaceae family. Since antiquity, carob trees have been wide-spread in most countries of the Mediterranean basin usually growing in mild and dry places having marginal and calcareous soils [38]. Today carob fruit has found a diverse number of uses in the food and pharmaceutical industry including gluten-free health products, stabilizers, thickeners and anti-cancer and anti-diabetes drugs and supplements [38,39,40,41]. Traditionally carob trees were interplanted along other staple crops such as olives, grapes and almonds and were part of low-intensity agro-silvopastural systems alongside grazing animals, and cereal cropping [38]. These systems have been traditionally applied in areas like Portugal (montado) and Spain (dehesa). Moreover, carob-based agroforestry has been wide-spread in the arid/semi-arid Mediterranean region including Tunisia [42] and Greece [43]. In Cyprus, carob trees have traditionally been part of sheep-grazing and cereal silvopasture systems [44].

Carob is increasingly recognized as a strategic crop for sustainable agricultural development in arid and semi-arid regions, due to its drought tolerance, low nutrient requirements, and capacity to stabilize degraded soils. These characteristics make it particularly suitable for addressing the challenges of water scarcity and land degradation in Mediterranean and other dryland environments [45]. However, ensuring the sustainability of carob production systems requires a comprehensive evaluation of key environmental indicators, including water footprint (WF), energy intensity (EI), and carbon footprint (CF). These metrics are essential for assessing the efficiency of water use, energy demand, and greenhouse gas emissions associated with production practices [46]. Recent studies have emphasized that integrating such sustainability indicators can support the development of low-input, climate-resilient farming systems, providing valuable insights for both producers and policymakers seeking to optimize resource management in arid areas [47]. Furthermore, life cycle assessment approaches applying these indicators have proven effective in identifying trade-offs between productivity and environmental impacts, enabling more informed decision-making for the sustainable management of carob-based agroecosystems [48].

Thymbra capitata [(L.) Cav.], [synonym of Thymus capitatus (L.) Hoffmanns. et Link.] is an aromatic plant adapted to extreme drought stress found in the thermo-Mediterranean vegetation belt and is typical of garrigue or phrygana vegetation. T. capitata makes very conservative use of soil moisture [49] enabling it to grow on marginal and saline soils [50]. Such properties make it ideal for dryland agroforestry systems as its competition with the tree crop is likely to be low. Thyme, when inter-cropped with fruit trees, has been found to increase soil organic carbon [30,51], improve soil aggregate stability and water availability in both top and subsoil [52], provide efficient weed management as a living-mulch [22] and provide for better soil cover compared to other aromatic plants due to its cushion life-form [53].

In this study rainfed intercropping systems of carob with thyme and clover were tested under arid and semi-arid conditions to address the following research questions:

(1) Can rainfed (i) thyme (Thymbra capitata; T-System) or (ii) clover (Trifolium squarrosum; C-System) intercropped with carob at SES of the tree, improve carob leaf-water content and tree productivity (yield efficiency and wood production) compared to rainfed conservation tillage carob plots (TLG_dry) in both semi-arid/low-weed and arid/high-weed competition conditions?

Our first hypothesis is that in semi-arid low-weed competition conditions, both thyme- and clover-based intercropping systems, when strategically planted in the SES of carobs, may reduce soil temperature/evaporation and retain higher levels of soil moisture during the dry season, without exercising too much competition pressure on carob trees. This may alleviate drought stress on carob trees and improved productivity, compared to rainfed conservation tillage systems.

Our second hypothesis is that the T-System may perform better than the C-System, due to the thymes’ capacity to consume less water, while providing SES cover all year round, thus lowering evapotranspiration rates compared to the C-System.

Our third hypothesis is that in arid high-weed competition conditions, we expect higher evapotranspiration and lower SES establishment rates for the T- and the C-System due to weed competition. This in turn may not benefit soil moisture availability for carob, leading to reduced tree productivity compared to conservation tillage.

(2) What are the pros and cons of the T-System compared to TLG_dry and to irrigated conservation tillage carob plots (TLG_irr) in terms of yields and environmental performance (CO₂ emissions, energy and water demand), in semi-arid low weed competition conditions?

We hypothesize that the T-System may outperform TLG systems in terms of environmental performance, while securing satisfactory yields, as a low input nature-based solution for adaptation of dryland agroforestry to climatic change.

2. Materials and Methods

2.1. Experimental Sites

The field experiments were conducted under rainfed conditions from 2019 till 2022 covering three dry seasons. Selected fields were from two sites on the island of Cyprus in the southeast Mediterranean that had been planted with carob trees at 7 m distance from each other (7x7 m planting density: 204 trees/ha) (Figure 1). The first site, named Skarinou (SKR), included two neighbouring agricultural plots in Larnaca district [(i) 34°48’37.68”N 33°21’20.57”E, (ii) 34°48’40.75”N 33°21’14.62”E] at ca. 189 m elevation, while having a semi-arid climate (Table 1). The second site named Vrysoulles (VRY), was in Ammochostos district (35°04’42.2”N 33°52’46.7”E) at 32 m elevation, experiencing an arid climate. Both regions have been classified as prone to desertification (Figure 1; see also Table 1). At site SKR 22 carob trees were selected of similar size and vigor at the age of 20. Since planting they had been drip irrigated at the base of their stem during the dry months up to summer 2019. Thereafter, irrigation stopped so as not to affect the experiments. At site VRY the 20 carob trees selected were also of similar vigor and size and all 18-years old. They had only been irrigated during the first two years of establishment. Prior to the experiment, conventional soil and weed management in both locations was performed, with 3 to 4 ploughing applications every year depending on seasonal rainfall and growth of weeds.

2.2. Intercropping Design

Around each of the 42 selected carob trees, a ca. 22-25 m² rectangular surface plot was delimitated. Within each plot, at the South-sun exposed soil side (SES) of each tree, a smaller surface area of ca. 3.6 m² (3 x 1.2 m) was selected (Figure 2) for application of Thyme (Thymbra capitata [(L.) Cav.]) or Clover (Trifolium squarrosum L.) to create a carob-thyme intercropping system (T-System; 19 plots in both sites) or a carob-clover intercropping system (C-System; 24 plots in total). The SES plot area was ca. 1.5 m away from the trunk (no midday shading from the tree on the intercropped crop). The remaining plot surface was dominated by wild annual vegetation (see Matsi et al., 2024 for floristic composition) and in less than 20% by vetch Lathyrus ochrus (L.) DC. No weeding took place in the plots throughout the duration of the experiment.

2.2.1. Carob-Thyme Intercropping System (T-System)

In each SES area, 60 thymes (ca. 3-5 cm tall) at a density of ca. 16 seedlings per m², were planted in late December 2019 till early January 2020, adding up to a total of 1140 T. capitata plants in 19 plots (8 plots in SKR; 11 plots in VRY). Prior to planting each thyme the soil within each planting hole was ecologically engineered by using a beneficial microbe complex and a water retaining mineral as described in [55].

2.2.2. Carob-Clover Intercropping System (C-System)

In addition to the 19 T-System plots, another 24 plots were planted with thyme (as described in section 2.2.1) and sown annually with T. squarrosum (55 g/m² of seed) in late autumn as described in (Figure S1). The initial aim was to produce a Carob-Thyme-Clover intercropping system. However, after the first dry season thyme survival rate was below 2% (Table 2). It appears that the ecologically engineered soil combined with the very wet climate of Y1 produced a clover overgrowth that eradicated the thymes. Therefore, these 24 plots were considered as carob-clover plots (C-system, 12 plots per site), since any actual effect of thyme on carob was negligible. However, clover’s establishment the following years highly depended on climate conditions. During Y2, the very dry winter (December 2020-February 2021; Annual rainfall: SKR= 425.3mm, VRY= 331.1mm, approx. 50% less than normal precipitation, Table 1) did not permit clover to sprout in any plot. As a compensation, a layer of clover hay was applied at a biomass density of 320-400 g/m² at the SES of each tree, to retain a clover mulch layer for the following dry season (Figure S1, Table 2). Moreover, in Y3, not all plots hosted clover since in some cases weeds appeared to outcompete it. Thus, the C-System had two variations. The CCC C-System, were plots with clover-based living mulch was available to the carob trees in all 3 monitoring years. In CCW S-System plots clover-based living mulch was available in Y1 and Y2, while in Y3 weeds dominated the CCW SES plot area. No weeding or management of clover, such as termination of the crop through cutting or ploughing, took place throughout the monitoring period.

2.3. Control Plots

The T- and C-System performance was compared to two types of control plots. In Control 1 plots conservation tillage was applied within the entire plot once a year without irrigation (TLG_dry, 8 plots in total; 4 per site). Control 2 included all remaining trees within site SKR (213 trees) where irrigation and conservation tillage were applied (TLG_irr).

2.3.1. Conservation Tillage Control Plots Without Irrigation (TLG_dry)

Wild vegetation in TLG_dry plots was managed by ploughing the plots down to a maximum of 15-20 cm once, during springtime for each experiment year. The aim was to reduce vegetation’s soil cover to below ca. 30% during summertime (average soil cover in TLG_dry plots for all three monitoring summer seasons was 21% in SKR site and 33% for VRY site). No irrigation was provided to TLG_dry plots throughout the duration of the experiment.

2.3.2. Irrigated Conservation Tillage Carob Trees (TLG_irr)

Wild vegetation in TLG_irr was managed by ploughing the trees down to a maximum of 15-20 cm once during springtime for each experiment year. Supplemented irrigation was provided to the trees during the dry months (May-October) throughout the experiment. In total 213 trees at SKR site were monitored.

2.4. Experimental Monitoring

2.4.1. Thyme Survival Rate

All living thyme plants (baring green leaves) were counted in each T- and C-System plot, measured three times during the duration of the experiment (February 2021, May 2022 and December 2022) to cover plant recovery from the dry seasons of 2020, 2021 and 2022. The thyme plant survival rate (TSR) was estimated as the percentage of living thyme plants compared to the total number of thyme plants planted during establishment. Note that TSR was below 2% after the first dry season (February 2021 assessment) for the C-System plots.

2.4.2. Thyme Biomass/Thyme Soil Cover/Total Soil Cover

Thyme biomass (TB) and thyme soil cover (TSC) were assessed following each dry season (in May of 2021, 2022 and 2023). For TB, an estimation was provided based on sampling thymes within the SES plots, aiming not to remove biomass amounts that could significantly disrupt any possible effects of thyme on soil microclimate. Thus, the minimum number of thyme plants e.g. 3 plants per SES plot at each site was randomly selected (57 plants in total) and their biomass above 4-5 cm from the ground was cut, dried and weighed (with a 2 decimal digital scale) in all T-System plots. The thyme dry biomass per site was averaged and calculated per m² (multiplying by 16 i.e. the number of planted thyme plants per m²). Since not all thymes survived, a correction was applied by multiplying the above biomass by the average TSR % per site during the corresponding monitoring period. After thyme cutting, the remaining shrub volume still on the ground (mostly xylem of the stems’ bases) was 40% (Y1), 30% (Y2) and 20% (Y3), based on the entire shrub’s volume assessment before the cutting. This biomass was estimated by multiplying the above % with the biomass measured from the remaining shrub and then added to that measured to produce the final TB per site.

TSC was estimated by measuring the elliptic cross section area of each thyme’s canopy as projected onto the soil for all thyme plants within a 1 m² quadrat. The quadrat was randomly placed within each of the T-System SES plots in each site location. To estimate the TSC, the distance between thyme branches furthest extending from the central stem (D1) as well as the distance perpendicular to D1 (D2), were used for calculating the elliptic surface of thyme canopy projected on the soil based on the equation: Thyme soil cover per plant = π*D1*D2 / 4. TSC per site was averaged and multiplied by the site’s TSR during that same monitoring period. During each thyme canopy monitoring session, the % area covered by wild vegetation was also estimated in each quadrat. This permitted to estimate the total soil cover by vegetation in each quadrat; i.e. the soil covered by both thyme and wild vegetation (TWSC).

2.4.3. Clover Biomass and Soil Cover

During each dry season (in June 2020, 2021 and 2022) clover biomass (CB) and total biomass from clover and weeds (CWB) present in a 1 m² quadrat randomly placed in SES plots were measured for 4-6 C-System plots. The CB/CWB ratio (in %) was calculated (Table 2). The total soil cover by clover and weeds on soil (CWSC) was also measured using a 1m² quadrat for 3-5 SES plots at each site.

2.4.4. Water Potential

The carob trees’ leaf water potential (Ψ) was monitored using a Schölander pressure chamber. Twig with leaf samples were collected between at midday from the southwestern (sun exposed) side of each sampled tree at 1.70-2 m from the ground. At least two samples were taken from each monitored tree (N sampled trees: SKR=15-17, VRY=13-15 each season). Water potential was monitored once at the beginning of the dry season (i, mid-June) and once at the end of the dry season and prior to the rainfall period (ii, beginning-mid October). A baseline measurement was also conducted prior to experimental establishment (October 2019). Recorded values (in bar) were transformed to MPa using the equation: Ψ (MPa) = Ψ (bar) * (-0.1). As a measure of dry season water stress on the carob trees, the leaf water potential difference ΔΨ between the beginning (i) and the end of the dry season (ii) (ΔΨ= Ψii - Ψi) for each monitored year was calculated.

2.4.5. Wood Production

To monitor wood production (secondary xylem growth) of each carob tree, the circumference of each experimental tree was recorded at a distance of 40cm perpendicular from the ground. Circumference (C, cm) was transformed to basal area (BA, cm²) using the equation: BA (cm²) = C² / 4 * 3.14. To estimate the percentage of wood production over the span of the three-year experiment (ΔBA_22-19), the difference of tree basal area between the third (BA₂₂: December 2022) and the baseline year (Prior to experimental establishment - BA₁₉: October 2019) was calculated and then divided by the baseline value based on the equation: ΔBA_22-19 (%) = (BA₂₂ – BA₁₉ ) * 100 / BA₁₉. Values were then normalized based on the site’s average ΔΒA (multiplied by 100 and divided by the site’s average) for comparison purposes between the two sites.

2.4.6. Carob Yield Efficiency

Carob fruit (pod) from each plot was collected and weighed at the end of August of each experiment year. The total carob yield per tree for all 3 years of production was then divided by canopy volume (CV) of the particular tree to estimate the trees’ yield efficiency in kg/m³ [56]. This correction was necessary to account for any differences in tree size. Trees were at a productive phase during the first 2 years of the experiments only in site SKR. Carob pod productivity at site VRY appeared to have been delayed, since these trees had been growing in more arid conditions and without irrigation apart from their first 2 years of establishment. Canopy volume (CV) was calculated based on the following equation: CV (m³) = π*D₁*D₂*[(Ht₁-Hs₁) + (Ht₂-Hs₂)]/12 [57]; with D1 being the distance of tree branches furthest extending from the trunk; D2 the distance perpendicular to D₁; Ht tree height and Hs the lowest point of canopy’s height from the ground. Measurements were made with the Nikon Forestry Pro II laser meter.

2.4.7. Soil Moisture Measurements

Soil moisture was monitored at the end of the third-year dry season (October 2022) down to 5 cm of soil using a EC5 Soil Moisture Smart Sensor with a HOBO USB Micro Station Data Logger (Onset). Soil moisture measurements were taken between 7:00am-12:30pm from 3 to 4 sun exposed spots positioned within the SES plots (N sampled trees: SKR=16, VRY=11). Measurements per spot were recorded every 10 seconds for a total of 5 minutes. Values were recorded in m³/m³ and were multiplied by a factor of a 100 to be expressed in %.

2.4.8. Surface Temperature Measurements

Surface infrared radiation was recorded during the peak of the third dry season (July 2022) using a Grade Thermal Camera (FLIR ONE Pro LT Pro) with a thermal resolution of 80x60 (4800px) and thermal pixel size of 17μm. One photo of the south facing side of each monitored plot was recorded between the hours 11:30-13:30. Thermal images were later processed using the FLIR Thermal Studio Starter Software. As experimental plot surface temperature was considered the average temperature of a 350px (width) by 97px (height) box frame (33,950px area). The box frame was always placed in the SES plot covering at least 50% of its surface, while avoiding any non-treatment related objects (e.g. large rocks).

2.5. Climate Data

As described in Matsi et al. (2024) climate data were used from the three meteorological stations closest to SKR and VRY sites. They belonged to the network of the Cyprus Department of Meteorology and were <5 km away from each experimental site. The 3 stations were Skarinou (station number: 628, 34°49’17.22”, 33°21’19.45”), Dipotamos (station number: 633, 34°51’08.24”, 33°21’40.47”), and Frenaros (station number: 1845, 35°02’32.64”, 33°55’03.414”) at 180, 175 and 78 m above sea level (a.s.l.). Due to missing data in Dipotamos station, temperature data for SKR site was estimated based on regression analysis (R² = 0.997) in 2022 using temperature data from the closest station (ca. 2 km distance from Dipotamos station, namely Kakoratzia station; number: 1634, 34°51’40.55”, 33°22’57.10”, at 231 m a.s.l).

2.6. Environmental Indicators

The environmental performance indicators 1) water footprint (WF - L of water per kg of carobs), 2) energy intensity (EI - MJ energy used per kg of carobs) and 3) carbon footprint (CF - kg CO₂ – equivalent per kg of carobs) were calculated, based on experimental data. The calculations were made for the SKR site and the treatments TLG_irr, TLG_dry and T-system, for two years (2021, 2022). In the experiment, the inputs that are related to quantifiable indicators are tillage and irrigation, where energy is consumed in the form of fuel and electricity to power the irrigation pump. In the case of thyme-based living mulch these inputs were zero, as tillage and irrigation were not applied. Additionally, the C stored in the wood biomass and the soil C balance were not considered. In the case of TLG_irr, data were obtained from the producer, related to irrigation schedule, network, water application and electricity consumption by the motor pump employed. The irrigation is typically applied every 2 weeks from May to October, and the applied irrigation dose is 100-200 L per tree, based on the seasonal variation of rainfall, which is also accounted for in the WF estimations, as the effective rainfall, which is 30% of the annual precipitation. For the EI, electricity and fuel use were considered. Total values were estimated for C emissions, Water and Energy consumption and, then, expressed per kg of yield. For electricity use, the emission factor of the electricity production in Cyprus, using the country’s energy mix (e.g., renewable, crude oil) was used, equal to 0.627 kg CO₂ per kWh (data from Electricity Authority of Cyprus). For the energy use and emissions related to tillage, data from Ecoinvent were used, adjusted to the equipment and time needed for soil cultivation in SKR site. Therefore, per tree, 4.404 MJ of energy are related to soil cultivation, which corresponds to emissions of 0.571 kg CO₂ eq. Table S1, provides the data used for the estimation of the environmental impact indicators.

2.7. Statistical Analysis

To test for normality of the data sets the One-Sample Kolmogorov-Smirnov Test was used. When data were non-normal, the Kruskal-Wallis test was conducted to compare medians. To compare means for normally distributed data one-way ANOVA was used followed by Tukey’s Honestly Significant Difference test (HSD) and Bonferroni test to analyse multiple comparisons between homogenous groups. The level of significance used was set at p < 0.05. All statistical analysis was carried out using SPSS Statistics 29 (SPSS, Chicago, IL, USA).

3. Results

Our experiment spanned over a three-year monitoring period. SKR site was classified as having a semi-arid climate under normal conditions with average Oct-Sep rainfall of 425mm during 1981-2010. However, plants in SKR experienced very wet conditions during Y1 i.e. 2019-20 (+54% than normal rainfall; 655 mm Oct-Sep rainfall, De Martone’s Aridity Index; DMAI= 21.29; Table 1). This was followed by a very dry climate in 2020-21 (Y2; -27% than normal rainfall; 311 mm Oct-Sep rainfall, DMAI= 9.96) creating arid conditions. In Y3 (2021-22) climate was close to normal (+17% than normal rainfall; 496 mm Oct-Sep rainfall, DMAI= 15.77). Climate, at the classified as arid site of VRY average Oct-Sep rainfall of 331mm during 1981-2010), was likewise very wet during the hydrological year 2019-20 (Y1; +55% than normal rainfall; 513 mm Oct-Sep rainfall, DMAI= 16.9), but arid for the two years that followed (Y2: 2020-21; -13% than normal; 289 mm Oct-Sep rainfall, DMAI= 9.4; and Y3: 2021-22; dry year, -15% than normal; 283 mm Oct-Sep rainfall, DMAI= 9.3; Table 1). Preceding our experiment, very wet conditions prevailed at both sites (SKR: +60% than normal rainfall; 678 mm Oct-Sep rainfall, DMAI= 22.28. VRY: +42% than normal rainfall; 469mm Oct-Sep rainfall, DMAI=15.35; Table 1).

Table 1. Annual climate data for SKR and VRY site by hydrological year (Oct-Sep). Exp.= Experiment, DMAI = De Martonne’s Aridity Index, r = rainfall in bold/italics for below/above the 1981-2010 mean, Ta = mean temperature.

Exp. : Hydrological year	r (mm)		% change to previous year		% of normal r		Ta (°C)		DMAI
	SKR	VRY	SKR	VRY	SKR	VRY	SKR	VRY	SKR	VRY
Y0: 2018-2019	678	469			160	142	20.4	20.5	22.28	15.35
Y1: 2019-2020	655	513	-<1%	+<1%	154	155	20.8	20.4	21.29	16.88
Y2: 2020-2021	311	289	-53%	-44%	73	87	21.2	20.8	9.96	9.40
Y3: 2021-2022	496	283	59%	-2%	117	85	20	20.3	15.77	9.33

Table 1. Annual climate data for SKR and VRY site by hydrological year (Oct-Sep). Exp.= Experiment, DMAI = De Martonne’s Aridity Index, r = rainfall in bold/italics for below/above the 1981-2010 mean, Ta = mean temperature.

Exp. : Hydrological year	r (mm)		% change to previous year		% of normal r		Ta (°C)		DMAI
	SKR	VRY	SKR	VRY	SKR	VRY	SKR	VRY	SKR	VRY
Y0: 2018-2019	678	469			160	142	20.4	20.5	22.28	15.35
Y1: 2019-2020	655	513	-<1%	+<1%	154	155	20.8	20.4	21.29	16.88
Y2: 2020-2021	311	289	-53%	-44%	73	87	21.2	20.8	9.96	9.40
Y3: 2021-2022	496	283	59%	-2%	117	85	20	20.3	15.77	9.33

The T-System at SKR site maintained relatively the same TSR during the three-year experiment, starting from a 97% TSR (Y1) and ending with a 94% TSR (Y3) (Table 2; Figure 3). At VRY, TSR dropped from 72 to 40%, three years after the T-System establishment (Table 2; Figure 3; [55]). Thyme biomass for the T- System in SKR, increased from 43 g/m² (Y1), and 335 g/m² (Y2) to 420 g/m² (Y3). Soil cover within the SES plots included thyme and unmanaged weeds (TWSC). It measured >40% (Y1), 55% (Y2) and 85% (Y3). TSC within the SES plots increased from 15% in Y1 and 35% in Y2, to 83% in Y3.

The T-System at VRY site, yielded a thyme biomass of 28 g/m² in Y1, increasing to 45 g/m² in Y3, reaching only at ca. 10% of that in site SKR. Soil cover of both thyme plants and wild vegetation (TWSC) in SES plots in VRY reached >97% in Y3, from >55% in Y1 (Table 2). TSC within the SES plots increased from 11% in Y1 and 9% in Y2, to 25% in Y3. Thus, unlike site SKR, wild vegetation (mostly Avena sp.) dominated SES plots in VRY site (Matsi et al. 2024). For the C-System in site SKR, the average biomass for 3 years in CCC plots (clover + weeds) was 288 g/m² (Table 2: Y1= 215 g/m², Y2=329 g/m², Y3>320 g/m²), with 70% of biomass therein produced by clover in Y1 and 55% during Y3. Soil cover (clover & weeds; CWSC) remained over 90% in CCC plots for all years of experiment (Table 2: Y1>90, Y2=97, Y3=100%). Total plant biomass (clover & weeds) in the C-System plots at VRY (CCW plots), remained higher than 300 g/m², during all experiment years (Table 2; Y1= 305 g/m², Y2>400 g/m² and Y3=316 g/m²), with clover biomass therein dropping from 50% in Y1, to 7% in Y3. In CCW plots soil cover (clover & weeds; CWSC) remained higher than 95% during all three years of experiment (Table 2; Y1=97, Y2>95, Y3=100%).

Table 2. Living mulch performance in the SES for T-Systems and C-Systems (CCC and CCW) at SKR and VRY sites during the three-year monitoring period (Y1=2019-20, Y2=2020-21, Y3=2021-22). Thyme and clover biomass (TB, CB). Clover and weed biomass (CWB). Soil cover (thyme without and with weeds, clover without and with weeds; TSC, TWSC, CSC, CWSC). Thyme survival rate (TSR). Clover to clover plus weeds ratio (CR). Exp.= Experiment. Weed comp.=Weed competition *Representing values for the SES layer of clover-based hay applied during Y2.

SKR site
			T system				C-System: CCC
Exp. year	Climate year	Weed comp.	TB (g/m2) mean	TSC (%) mean	TWSC (%) mean	TSR (%) median	CB (g/m2) mean	TB (g/m2) mean	TSC (%) mean	CWB (g/m2) mean	CR (%) mean	CWSC (%) mean	TSR (%) median
Y1	Very wet	Low	43	15	>40	97	157	1	<0.3	215	70	>90	2
Y2	Very dry	Low	335	35	55	97				329*		97*
Y3	Normal	Low	420	83	85	94	>180			>320	55	100
VRY site
			T system				C-System: CCW
Y1	Very wet	High	28	11	>55	72	171	0	0	305	50	97	0
Y2	Very dry	High	39	9	>80	40				>400*		>95*
Y3	Very dry	High	45	25	>97	40	22			316	7	100

Table 2. Living mulch performance in the SES for T-Systems and C-Systems (CCC and CCW) at SKR and VRY sites during the three-year monitoring period (Y1=2019-20, Y2=2020-21, Y3=2021-22). Thyme and clover biomass (TB, CB). Clover and weed biomass (CWB). Soil cover (thyme without and with weeds, clover without and with weeds; TSC, TWSC, CSC, CWSC). Thyme survival rate (TSR). Clover to clover plus weeds ratio (CR). Exp.= Experiment. Weed comp.=Weed competition *Representing values for the SES layer of clover-based hay applied during Y2.

SKR site
			T system				C-System: CCC
Exp. year	Climate year	Weed comp.	TB (g/m2) mean	TSC (%) mean	TWSC (%) mean	TSR (%) median	CB (g/m2) mean	TB (g/m2) mean	TSC (%) mean	CWB (g/m2) mean	CR (%) mean	CWSC (%) mean	TSR (%) median
Y1	Very wet	Low	43	15	>40	97	157	1	<0.3	215	70	>90	2
Y2	Very dry	Low	335	35	55	97				329*		97*
Y3	Normal	Low	420	83	85	94	>180			>320	55	100
VRY site
			T system				C-System: CCW
Y1	Very wet	High	28	11	>55	72	171	0	0	305	50	97	0
Y2	Very dry	High	39	9	>80	40				>400*		>95*
Y3	Very dry	High	45	25	>97	40	22			316	7	100

In our three year experiment one heavy (Y1), one poor (Y2) and one medium (Y3) fruit-bearing year was recorded (Table S2). During the harvest of 2020 at the end of a very wet hydrological year (Y1; Table 1), carob yield efficiency was greater for the TLG_dry (1.27 kg/m³) and the T-System carob trees (1.11 kg/m³), with irrigated conservation tillage plots (TLG_irr) coming last (0.95 kg/m³). However, for the harvests of 2021 and 2022, towards the end of a very dry (Y2) and a normal hydrological year (Y3), the carob yield efficiency of TLG_irr appeared related to the water supply (Table S3: TLG_irr; Yield efficiency: 2021=0.92 kg/m³, 2022=0.59 kg/m³ , Irrigation water: 2021= 43.23mm, 2022= 62.27mm) Yield efficiency for the T-System during the very dry Y2 was 0.26 kg/m³, amounting to 28% of TLG_irr plots (Table S3). Nevertheless, at Y3, when climate was normal, the yield efficiency of the T-System reached 0.33 kg/m³, which equaled to 56% of TLG_irr carob tree productivity. Note, that during Y2 and Y3, carob pod production was practically zero for TLG_dry plots (Table S3; Yield efficiency 2021< 0.01 kg/m³, 2022 no production).

Prior to treatment initiation, midday leaf water potential of the trees belonging to the 3 carob tree groups tested (T-System, C-System and TLG_dry) was measured at the end of the dry season (October 2019 – Y0). No significant differences were identified in tree water status among the 3 groups at both sites (Figure 4a and b). However, in Y3 (October 2022) at SKR site, carob trees of the T- and C- System plots had significantly improved tree water content, compared to TLG_dry carob plots (Figure 4c, p < 0.001). At VRY site, no signigicant differences in the trees’ water content were identified between the intercropping systems and the TLG_dry plots (Figure 4d). Similar results were produced when the same trees were compared during Y0 and Y3 monitoring periods (N for SKR: T=3, CCC=3, TLG_dry=3 & VRY: T=4, CCW=4, TLG_dry=3).

For the wet (Y1) and the normal (Y3) climate year, midday leaf water potential differences, ΔΨ, between the beginning and the end of the dry season for T-System carob trees at SKR site, was signigicantly higher compared to TLG_dry i.e. the unirrigated trees under conservation tillage (Figure 5a p < 0.01; Figure 5c p < 0.001). During the very dry year (Y2), no difference in the trees’ water content was identified between the two systems (Figure 5b, p > 0.05). Carob tree yield efficiency added up for Y1, Y2 and Y3 at SKR, was significantly greater for the T-System compared to TLG_dry, with productivity increase ranging from 27% to 71% (based on standard error margins; mean 47%; Figure 5d, p < 0.05). It concided with 12.5% higher carob wood production over the 3 year period (based on carob stem basal area change) compared to TLG_dry, although not statistically different (Figure 5e, p > 0.05).

At SKR, midday leaf water potential differences, ΔΨ, between the beginning and the end of the dry season for C-System carob trees growing under the effect of clover-based living mulch in all 3 years of monitoring (CCC), was signigicantly higher (Figure 6c, p < 0.05) on Y3, while having no negative impact for Y1 and Y2 compared to TLG_dry (Figure 6a and b, p > 0.05). CCC carob trees appeared to produce a tendency towards improved yield efficiency compared to TLG_dry considering all 3 years of monitoring, although not statistically significant (Figure 6d, p > 0.05). Carob wood production for both systems was again not statistically different over the 3-year monitoring period (Figure 6e, p > 0.05).

At VRY site, no significant differences in midday leaf water potential, ΔΨ, between the beginning and the end of the dry season were recorded for the T-System carob trees compared to TLG_dry , for all 3 years of monitoring (Figure 7 a-c, p > 0.05), as was also the case for carob wood production (Figure 7d, p > 0.05). C-System carob trees growing under the effect of clover-based living mulch in 2 out of the 3 years of monitoring (CCW) exhibited the same pattern as above (Figure 8, p > 0.05).

During Y1 and following a very wet hydrological year (Y0=2018-2019, +60% than normal rainfall, Table 1) at site SKR, no differences in carob tree yield efficiency were recorded among all systems investigated (Figure 9a, p > 0.05). Following a second very wet hydrological year (Y1=2019-2020, +54% than normal rainfall, Table 1), during Y2 of monitoring, carob tree yield efficiency was significantly higher at T- System plots compared to C-System and TLG_dry plots (Figure 9b, p < 0.05). In Y3 and after a very dry wet season (Y2=2020-2021, -27% than normal rainfall, Table 1), a tendency for higher carob tree yield efficiency was noted for T-System and C-System CCC plots, although not statistically significant compared to TLG_dry (Figure 9c, p > 0.05). For Y1 and Y2, the C-System included trees growing under the influence of clover as a cover crop within the SES plots (Figure 9a,b). However, since no clover sprouted in 8 of the 11 plots in Y3, the C-System was seperated into two groups CCC (N=3) and CCW (N=8) (Figure 9c).

Measurements within the dry season of Y3 indicated that both the T- and the C-System CCC at SKR site can significantly reduce the plots’ midday thermal stress by 9 ^o C and 13^o C compared to TLG_dry (Figure 10a p < 0.001, Figure 10b p < 0.01). However, no such effect was detected in site VRY (Figure 11a and b, p > 0.05). At site SKR, at the end of the same dry season, soil moisture in both the T- and the C- System plots was found significantly higher by 18% and 3% respectively (Figure 10c p < 0.001, Figure 10d p < 0.001) compared to TLG_dry plots. In site VRY, the opposite was detected for plots of the C-System CCW, with soil moisture found to be significantly lower by 4% compared to TLG_dry plots (Figure 11d, p < 0.001), while no significant difference in soil moisture was recorded betweenT-System plots and TLG_dry plots (Figure 11c, p > 0.05).

For the WF the values for Y2 were 2.0, 0.8 and 38.5 x 10³ L/kg carob pods for the treatments T, TLG_irr and TLG_dry (Figure 12a). The respective values for Y3 were 0.7 and 0.8 x 10³ L/kg carob pods for T and TLGirr. Regarding EI, the values for Y2 were 0, 2 and 18 MJ/kg carob pods for the treatments T, TLG_irr and TLG_dry (Figure 12b). The respective values Y3 were 0 and 5 MJ/kg carob pods for T and TLG_irr. The WF, EI and CF for TLG_dry could not be determined for Y3, as there was no production (n.p.). Finally, the CF values for Y2 were 0, 0.4 and 2.3 kg CO₂-eq/kg carob pods, for the treatments T, TLG_irr and TLG_dry (Figure 12c). For the Y3 these were 0 and 0.8 kg CO₂-eq/kg carob pods, for the treatments T and TLG_irr. Note that in Y3 the T-System outperforms all other systems for all indicators evaluated.

4. Discussion

Although conservation tillage is often praised for its benefits compared to conventional ploughing [58], it is still based on disrupting soil structure, a practice far from sustainable when it comes to promoting soil health and to exploiting the ecosystem services of the soil food web [59]. More efficient nature-based farming solutions are currently available that help conserve water and soil resources [6,30], enhance yields [17,27] and make better use of land overall [25]. Complementarity in resource use, water and nutrient uptake, are the foundations for the benefits in intercropping systems [60,61]. The alignment of root system traits, in both space and time, for intercropped species can be key in explaining how facilitation can outweigh competition [60], similarly to how plant facilitation in unmanaged dryland ecosystems is based on both root and canopy traits, capable of combating desertification phenomena [62,63,64].

For promoting alternatives to soil tillage, we designed agroforestry systems based on perennial (thyme) and annual (clover) living mulch. They were established only on the sun-south exposed soil side (SES) of carob trees and in only ca. 15% out of a total surface plot around each monitored tree. This would reduce as much as possible any[6,30] competition effects for soil water with the carob trees, while providing soil cover that would reduce thermal stress and evaporation where it would be the highest around the tree (Figure 2 and Figure 3). Our systems would also promote wild vegetation establishment around the trees, which, if left unmanaged at least for a large part of the year, can contribute to increasing the systems’ biodiversity [65,66] compared to conservation tillage. Moreover, the more perennial cover from aromatic plants expands in the SES of carobs, the more reduced weed competition and increase in beneficial insects can be expected as shown for Mediterranean olive intercropped systems [22].

4.1. The Drivers for the Greater Carob Tree Productivity in T-System Plots

Our findings are in line with our first hypothesis that in semi-arid low-weed competition conditions, as is the case for site SKR, the thyme inter-cropping system, strategically planted at the SES of carob trees (T-System), may reduce soil temperature/evaporation and retain higher levels of soil moisture during the dry season without exercising too much competition pressure on carob trees. This may alleviate drought stress on carob trees and improve productivity, compared to rainfed conservation tillage systems (TLG_dry). However, this was not found to be the case for the clover living mulch system (C-System).

At site SKR, T. capitata at the SES of carob trees achieved very high thyme establishment rates (Table 2: TSR> 90%) due to lower weed completion and more favourable climate [55]. This living mulch appears to have contributed towards significantly improving carob yield efficiency compared to TLG_dry (Figure 5d, p < 0.05). Thyme biomass in SES plots surpassed 300g/m² in Y2, 400g/m² in Y3 and dominated plant soil cover, achieving from 35% TSC in Y2, more than 80% TSC in Y3 (Table 2). Using perennial species planted in high densities, seems to be key for the significant improvement of soil properties and water content in intercropping systems under semiarid rainfed conditions [30] in line with our findings.

However, in Y1 during the T-System’s establishment, the >40% soil cover was achieved mostly from wild vegetation (>25%) with the contribution of thyme at 15%. The dominant type of wild vegetation at site SKR in Y1 (above 10% plant population abundance; p. ab.) was composed of Avena sterilis subsp. ludoviciana with p. ab. of 33.3 % and presence in 63% of the sampled plots, Crepis foetida subsp. foetida with p. ab. of 18.1 % and presence in 50% of the sampled plots, followed by Lactuca serriola (p. ab. of 10.9 % and presence in 63% of the sampled plots) and Lolium perenne (p. ab. of 9.4 % and presence in 100% of the sampled plots; see [55] for details). Such plant species and abundances combined with T. capitata do not appear to have produced a negative effect on carob yields.

For T-System carob trees at site SKR we recorded an improvement in carob yield efficiency by at least 27% and up to 71% compared to TLG_dry control plots (based on the standard error range of yields, mean 47%; Figure 5d), when considering performance from all 3 years of monitoring. It coincided with a 12.5% higher carob wood production over the 3-year period (based on carob stem basal area change) compared to TLG_dry (Figure 5e). However, in Y1 during the T-System’s establishment and under a very wet year, no significant yield increase was evident compared to TLG_dry. The improvement occurred in Y2 and Y3 (Figure 9). To explain these results, it is critical to comprehend the phenology of carob pod production. It appears that in carob trees, flower induction occurs at the beginning of the dry season, while blooming and pollination takes place towards its end, one year in advance of the carob pods’ harvesting year (Figure S2). We consider that the water made available to the tree, the winter months before the dry season, determines the number of flowers produced for fertilization thereafter. Thus, it can be related to the number of mature pods that can be harvested the next dry season. This considering that flower production and pollination remains relatively undisturbed by other environmental factors, the soil is well supplied in nutrients [67,68] and the hydrological year just before pod harvesting provides a minimum amount of soil moisture needed for pod elongation [67]. It appears that the above criteria are met in our experiments. Our trees were regularly supplied with manure, extreme rainfall events did not occur during the dry season, whereas rainfall during the hydrological year coinciding with the year of harvest was above 300 mm (a threshold above which carob survival does not appear to be threatened). Thus, plenty of rainfall during the wet season of Y0, i.e. the year before experiments initiated, would be expected to provide abundant flowering and the potential to produce a high carob yield during the harvest of Y1. This hypothesis is fully in line with our findings.

The very high yields of Y1 (2020) harvest (Tables A2 and A3, Figure 9) can be explained by the very wet year of Y0 (on top of which irrigation also took place during the dry season). No difference was recorded between T-System and TLG_dry carob pod harvest in Y1. It coincided with no difference in leaf water potential between the two systems at the end of the dry season of Y0 (Figure 4). Considering the above carry-over effect, the Y2 harvest of 2021 is actually the first year’s harvest on which our experimental trials could have had an impact on productivity, by effecting soil moisture of Y1. The significant improvement in carob leaf-water content in T-System plots compared to TLG_dry , during the dry season of Y1, reflects the soil moisture conditions the trees experienced during the wet and dry season of Y1. It coincided with a significant carob yield efficiency increase by 0.256 kg/m³ (p < 0.05), on average, in Y2 for T-System plots compared to TLG_dry (Table S3; Figure 9b). Soil cover at the SES of T-System plots was >40% during the dry season of Y1, when considering both thyme and the type of wild vegetation present at site SKR that produces low competition pressure, as mentioned above. Moreover, an improvement in carob leaf-water content was recorded during the dry season of Y2 in T-System plots compared to TLG_dry, at the end of a very dry hydrological year. It coincided with a higher carob yield efficiency by 0.331 kg/m³ on average in Y3 (Table S3; Figure 9c). Soil cover at the SES of T-System plots was 55% during the dry season of Y2, when considering both thyme and the SKR site wild vegetation.

The above findings are in line with the significantly improved leaf-water content during the dry season, for trees growing under the effect of the T-System compared to TLG_dry in Y3 at TSC of > 80% (Figure 5)_. In Y3 the contribution of wild vegetation on plant soil cover was negligible (ca. 2%); with T. capitata practically eliminating any competition from wild vegetation in site SKR conditions (Table 2; Figure 3; Figure 5c ; [55]). Therefore, it appears that T. capitata based living mulch most likely has the potential to reduce the loss of soil moisture due to evaporation during the dry season, out of which the carob trees can benefit from.

T. capitata in its mature life form, develops into a typical dryland cushion plant. The ability of cushion plants to buffer ecological extremes is a well-known ecological principal [69,70,71,72]. In an alley cropping experiment of thyme and almond trees conducted under rainfed Mediterranean conditions, thyme plants contributed to a significant increase in soil moisture compared to tillage application [30]. In our case, thyme-based living mulch measured within the dry season of Y3 displayed its ability to significantly reduce the plots’ midday thermal stress by 15% (-9 ^o C) and improve soil moisture content by 18%, compared to TLG_dry (Figure 10a and c, p < 0.001). Except from any cooling produced from evaporation of essential oils (water transpiration being less likely due to leaf stomata closing in mid-day or leaf desiccation in summer) [73], reflectance of solar radiation by stems/brown leaves and shading of soils appears as the most likely mechanism that protects soil from extreme heat stress in summer time. This in turn most likely leads to less evaporative forcing on soil, explaining its significantly higher soil moisture content in the dry season compared to TLG_dry soil. The T-Systems performance could even be compared with that of irrigated conservation tillage plots (TLG_irr). In Y2, yield efficiency for the T-System equalled almost a third of the yield of TLG_irr (Table S3: 2021 - Yield efficiency= 0.26 kg/m³, TLG_irr= 0.91 kg/m³). As the T-System’s soil cover became further established, yield efficiency of T-System carob trees in Y3 was more than 50% of that provided by TLG_irr plots (Table S3: 2022 - Yield efficiency: T= 0.33 kg/m³, TLG_irr= 0.58 kg/m³).

However, the above effect may become even stronger as the thyme canopy further closes up, permiting more efficient soil shading, while acting as an excellent means to control weeds [17]. Our preliminary analysis of comparing soil temperature on thymes canopies vs. bare soil between thyme rows indicated lower surface temperatures by 22 % on thyme rows (40.3 °C median), compared to the 51.4 °C (median) recorded on the bare soil surface in between rows at site SKR in July 2022 (Table S4 ; Figure S4) ; a significant difference of 11.1 °C (p <0.05). At site VRY these differences were even greater with 30% lower surface temperatures (45.3 °C median) monitored on the surface of thyme plants compared to the 64.7 °C (median) recorded on bare soil in between thyme rows; a recorded difference of 19.4 °C (p <0.05; Table S4). Thus, as the canopy of the aromatic living mulch grows its ability to mitigate soil temperature extremes should increase [74,75].

In terms of explaining wood production, one needs to consider that the rooting strategy of carob permits uptake of water from deep soil layers ([76,77] and references therein). Deep moisture pools may be supplied by long integration periods of rainfall, even of multiple years, as is the case for other tree species of the same climate zone (e.g. pines ; [78,79]). One of the first dendroclimatological analysis ever performed on carob tree rings was conducted in Cyprus based on trees from sites SKR and VRY [80]. It indicated that although annual rainfall could explain part of tree growth variability in site SKR, a significant effect of accumulated precipitation of 4 consecutive years, before and including the year of tree growth, could also be found in trees of both sites. Thus, apart from annual or inter-annual precipitation [81] it could be that multiple years of past rainfall [82,83] may also influence carob productivity. This may explain why the tree growth increase observed at site SKR was not as great as the carob yield efficiency for our 3-year monitoring period; with a 4th monitoring year likely needed to capture the full effect of our experiments on carob tree growth. It could also be the case that the T-System carob trees first invest in fruit productivity, and thus less resources may be allocated for simultaneous wood growth. Research conducted on carobs by von Haselberg [84] demonstrated that increased fruit yield came at a cost for vegetative growth. Specifically, during years of heavy fruit bearing (“on” years) trunk growth was significantly lower compared to years of light carob yield (“off” years). It is not yet well known how extended periods of drought impact carob yields [67]. Results of a recent review on carob tree publications of the last 40 years, indicate the need for further research on the correlation between carob yield and climate [45].

Carob’s alternate bearing is a known phenomenon [56]-, [85]. However, our results do not appear to be determined by it primarily. After a heavy fruit-baring year (Y1: “on” year) we would expect yields to decline (Y2: “off” year) and then increase again during the third monitoring period (Y3). Despite Y1 being a heavy yield year for TLG_dry plots, their yields for the following two years approached zero. These results are in line with the significant reduction in SES soil moisture and leaf-water content in TLG_dry plots compared to our living mulch experiments (Figure 5, Figure 6 and Figure 10). More specifically, the significant difference recorded in yield efficiency between the T-System and TLG_dry control plots in Y2, could be attributed to T-System plots being in “on” and TLG_dry in “off” status (Figure 9b, Table S3). If that would have been the case though, we would expect TLG_dry control plots to have an “on” year in Y3 – something we did not note (Figure 9c, Table S3). Therefore, although carob yield has been mostly explained by the status of the trees’ alternate bearing [86], environmental stress, such as intense drought stress, can also significantly reduce carob fruit yield [85]. In our case, water availability seems to be the limiting factor for carob productivity. When irrigating conservation tillage plots (TLG_irr), carobs at site SKR were able to maintain a relatively stable yield efficiency through fluctuating climatic conditions (Table S3: Yield efficiency: Y1= 0.93 kg/m³, Y2= 0.92 kg/m³, Y3= 0.58 kg/m³). Hence, under specific conditions, out of which the lack of water and nutrient uptake are key, the status of the tree’s alternate bearing can play a weaker effect on carob yield [85].

Our second hypothesis was that the T-System may perform better than the C-System, due to the thymes’ capacity to consume less water, while providing SES cover all year round, thus lowering evapotranspiration rates compared to the C-System. Although the T-System did outperform the C-System in terms of improving leaf-water content and carob tree yield efficiency, the handicap of the C-System does not appear to be the high completion for water with the carob tree. In fact, the SES plot covers a small surface for such potential effects to become evident. It is more likely that the C-System’s handicap is related to its high dependence on soil moisture in order for it to produce enough biomass/soil cover that can efficiently protect soil from high thermal stress and evaporation. In this sense the perennial T-System with its all-year capacity to provide soil cover has the advantage.

More specifically, carobs of the CCC version of the C-System at SKR site retained more leaf-water during the dry season compared to TLG_dry only during the normal climate year (Y3) (Figure 4 and Figure 6c). Competition for water between the carob and clover does not appear to be a problem, since in Y3 the highest clover biomass (CB) and clover-weed biomass (CWB) were observed (>180 g/m² for CB and >320 g/m² for CWB; Table 2). The 3-year cumulative effect led to greater carob yields compared to rainfed conservation tillage plots, but the difference was not significant (Figure 6d). No significant improvement in carob yield efficiency for the C-System could be found in Y2 nor in Y3 (Figure 9), while no increase in carob wood production could be identified (Figure 6d). Hence, T. squarosum failed to propagate autonomously and create living mulch in all climate conditions tested, while adding dry clover mulch to compensate for this deficiency did not improve carob leaf-water content during Y2, the driest of the years tested. It could be a theme of further investigation if other more drought adapted clovers or legumes, particularly perennial ones could perform better in order for the carobs or other dryland tree crop systems to be able to benefit from the legumes nitrogen fixing capacity [87,88].

Our third hypothesis was that in arid high-weed competition conditions of site VRY, we expect higher evapotranspiration and lower SES establishment rates for the T- and the C-System due to weed competition. This in turn may not benefit soil moisture availability for carob, leading to reduced tree productivity compared to conservation tillage. Although, both systems’ establishment was poor (Table 2; Figure 3), carob leaf-water content did not appear to be significantly affected (Figure 7 and Figure 8) at 25% SES soil cover and 45 g/m² T. capitata biomass nor at 7% clover abundance (Y3). No negative impacts on carob productivity, could also be identified. Overall, our interventions cannot be considered as detrimental for carob systems even in arid conditions, when it comes to somewhat improving the systems biodiversity or in providing some biologically fixed nitrogen to the carob trees. The coexistance of trifolium and unmanaged wild vegetation in the plots improves biodiversity indices compared to conservation ploughed plots [89,90]. Moreover, low-weed competition conditions seem to be key for the successful establishment of rainfed carob-thyme intercropping system [55].

There is evidence to suggest that the stress gradient hypothesis is not only limited to natural ecosystems, but also applies to agroecosystems [91], proposing that the benefits of intercropping can be further accentuated under conditions of environmental stress [92,93,94]. When comparing intercropped to monoculture systems under equal nutrient inputs, the first achieved higher resource use efficiency and yields – almost 50% more [91]. The difference, though, significantly declines as inputs become more abundant. Intercropping seems to become even more critical under water-stressed conditions [93]. In a water limited environment, an intercropped system was able to improve productivity by more than 90% compared to monoculture. In the case though of high-water availability, the same study identified only 14% improvement in the intercropped system yields in comparison to monoculture [94].

4.2. Overall Yield and Environmental Performance

We hypothesized that the T-System may outperform TLG systems in terms of environmental performance, while securing satisfactory yields, as a low input nature-based solution for adaptation of dryland agroforestry to climatic change. Our results fully support our original hypothesis. Productivity of rainfed carobs can range between 6,000 and 7,400 kg/ha under arid and semi-arid conditions [85] with as low as 400 kg/ha in drier sites [84]. Average carob pod yield efficiency in a four-year experiment conducted in rainfed carobs in Portugal under soil tillage management, ranged between 0.20 and 0.54 kg of pods per m³ canopy with less than 0.13 kg/m³ recorded in the site with the lowest annual rainfall [84]. In Y2, yield efficiency for the T-System at site SKR was 0.26 kg/m³ and increased to 0.33 kg/m³ in Y3 (Table S3), very well within the range or even above the performance of the tillage system reported above.

However, it is the environmental performance indicators of the carob production system that reveal the distinct profiles for each treatment. For the Water Footprint (WF) calculated considering both rainwater and irrigation, the TLG_dry treatment demonstrated great water waste during the harvest of Y2, with yield efficiency as low as <0.01 kg/m² (Figure 12a). The following year the T-System’s yield efficiency was the greatest, while TLG_dry had no production at all (Table S3). Therefore, for the later, WF could not be determined (Figure 12). WF data for carob trees are scarce. Yet, the WF values observed in this study for the T-System were lower than those presented in [95] for rainfed olive trees in Apulia, Italy. For Energy Efficiency (EI), the T-System consumed zero energy for its operation (Figure 12b), since no tillage and irrigation were applied and the same was observed for its Carbon Footprint (CF) measuring zero (Figure 12c). Nonetheless as the T-System has the potential to capture and store carbon into the soil, its long term carbon footprint is most likely to be negative.

The range of EI reported for grapes (2,5 - 4,2 MJ/kg) [96] is similar to that observed for TLG_irr, while the CF estimated for carob pods in the case of TLG_irr is comparable with the values estimated for indigenous grape varieties in Cyprus [96]. The TLG_dry treatment had the highest CF and EI values (Figure 12b and c), indicating a non-sustainable management system characterized by inputs but practically no yield under our experimental testing. Overall, our results highlight a clear trade-off between environmental impacts and yield stability. While irrigation increases water, energy, and carbon inputs, it supports more consistent productivity. However, irrigation water is in short supply in Cyprus as in dryland regions of the Mediterranean and beyond, comprising a significant limiting factor for carob production.

Hence, the intercropping strategy provided by the T-System should be further investigated for its promising performance in optimizing soil management and irrigation inputs, balancing resource efficiency with yield stability, parameters critical for sustainable carob and similar crop (e.g. olive) production in drylands under climate change.