Next Article in Journal
Different Functional and Taxonomic Composition of the Microbiome in the Rhizosphere of Two Purslane Genotypes
Next Article in Special Issue
Optimizing Crop Water Productivity in Greenhouse Pepper
Previous Article in Journal
A Systematic Review on the Improvement of Cd Stress Tolerance in Ramie Crop, Limitations and Future Prospective
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Setting Irrigation Thresholds for Building a Platform Aimed at the Improved Management of Citrus Orchards in Coastal Syria

1
LEAF—Linking Landscape, Environment, Agriculture and Food Research Center, Associated Laboratory TERRA, Instituto Superior de Agronomia, Universidade de Lisboa, Tapada da Ajuda, 1349-017 Lisboa, Portugal
2
General Commission for Scientific Agriculture Research (GCSAR), Hejaz Station, Damascus P.O. Box 113, Syria
3
Centro de Ciência e Tecnologia do Ambiente e do Mar (MARETEC-LARSyS), Instituto Superior Técnico, Universidade de Lisboa, Av. Rovisco Pais, 1, 1049-001 Lisboa, Portugal
*
Author to whom correspondence should be addressed.
Submission received: 1 June 2023 / Revised: 26 June 2023 / Accepted: 2 July 2023 / Published: 4 July 2023
(This article belongs to the Special Issue Assessment and Mapping of Soil Water Balance)

Abstract

:
Citrus is one of the most valuable crops in Syria, with the largest production areas being in the coastal provinces of Tartus and Latakia, where this study was performed. A companion paper reported on the basal crop coefficients derived from the field water balance and on the performance assessment of various irrigation methods used in a citrus orchard located in the same region. That study evidenced the need for the improved management of irrigation water, mainly reducing water applications and increasing productivity, thus leading to the current research. The main objectives consisted of (i) providing a set of reliable basal (Kcb) and average (Kc) crop coefficients to be used in practice in the citrus orchards of the Syrian coastal area, while accounting for the diversity of characteristics observed; (ii) to estimate the seasonal consumptive use of typical orchards under different climate-demand and deficit-irrigation scenarios; and (iii) to assess possible water savings and related yield reductions. The previously calibrated water balance model SIMDualKc was used for these purposes. The computed Kcb values for the mid-season and average demand for water ranged from 0.52, when the plant density was low, to 0.84, when plant density was very high. The corresponding Kc values, which further reflected the impact of drip irrigation on controlling soil evaporation, were 0.72 and 0.97, respectively. Overall, the consumptive use of water was estimated to range from 867 to 1573 mm. The assessed water-saving scenarios consisted of adopting increased management-allowed depletion (MAD) thresholds relative to the p depletion fraction for no stress: MAD = 1.05, 1.10, 1.20, and 1.30 p. For trees under a very high climatic demand, water savings ranged from 12 to 34%, but the yield losses induced by the water deficits ranged from 8 to 48%. Although the selection of optimal strategies should be based upon economic terms, these may only be used when the Syrian economy recovers from civil war and the current crisis. The present results show the feasibility of adopting such MAD thresholds for building an irrigation management platform. The data provided by the current study are valuable because they can be efficiently used to support of the irrigation management of Syrian citrus production systems.

1. Introduction

Syria is one of the main citrus producers in the Mediterranean region, with orange production averaging 0.67 M tons year−1 and lemon and lime yielding 0.36 M tons year−1 from 2016 to 2020 [1]. Before war times, citrus production counted as an important income source for the country, representing 1.3% of the gross domestic product, 20% of the value of national fruit and vegetable exports, and 0.8% of the world global production [2]. With the war, exports drastically dropped, and citrus production has been mostly for domestic consumption [3].
The Syrian citrus cultivated area is mainly located in the Tartus and Latakia provinces, totalizing close to 42,700 ha. In the coastal region, these provinces combine favorable environmental conditions for citrus production, namely mild winters, high humidity for most of the year, annual rainfall often above 800 mm, and water availability for irrigation during the dry summer season [2]. Drip irrigation is used in 43% of the irrigated land area [4], greatly contrasting with the rest of the country, where traditional irrigation methods are still the common rule.
In an effort to rationalize agricultural water use in Syria, several studies were carried out over the last decades to rationalize agricultural water use and control the overexploitation of available water resources [5,6]; improve irrigation water management practices through the accurate estimation of crop water requirements, using the soil water balance SIMDualKc [7,8]; increase land and water productivity and farmers’ income [9,10]; and assess the water, yield, and economic performance of modernized irrigation systems [11,12,13].
The FAO56 method remains the most widely used approach for estimating crop water requirements, i.e., the crop evapotranspiration (ETc) [14]. The ETc is the product of the grass reference evapotranspiration (ETo), representing the evaporative demand of the atmosphere and defined with the FAO Penman–Monteith (PM) equation, and a crop coefficient (Kc), representing the effect of the difference in terms of ET rate between the crop under study and under pristine, eustress conditions, and the reference crop, for which the ET rate is the ETo. Kc values, which are described by a Kc or a Kcb curve, vary with crop species and crop stages [14]. When adopting the single Kc approach, an average value combining soil evaporation and crop transpiration is assumed. If ETc is partitioned into the basal transpiration coefficient (Kcb) and the soil evaporation coefficient (Ke), the dual crop coefficient approach is adopted (Kc = Kcb + Ke).
The FAO56 methodology is somewhat challenging for orchards and grape vines because these are complex systems, with heterogeneous surfaces and incomplete ground cover. Differences in the planting density, canopy height, training system, interrow management, and irrigation method influence the amount of energy available for both the transpiration and the soil evaporation processes [15]. In addition, the best crop management often corresponds to adopting controlled/regulated water deficit at given periods for water saving and enhance the quality of fruits [16], referred to as the eustress conditions. To more easily account for the variability of conditions influencing Kc in orchard systems, Allen and Pereira [17] developed an approach (hereafter denoted as the A&P approach) to compute the Kcb from observations of the fraction of the ground covered by the vegetation (fc) and plant height (h) that estimate a density coefficient, as well as by the degree of stomatal adjustment (Fr) and the opacity of the canopy (ML). Pereira et al. [18] revisited the A&P approach and tested its performance relative to a large number of annual and perennial crops, comparing the Kcb values obtained with the A&P approach with those derived from ET field research. The results, namely for trees, showed excellent predictions of the Kcb by the A&P approach. In their study, Pereira et al. [19] tabulated the parameters to be used with the A&P approach when targeting standard Kcb values. Moreover, Pereira et al. [18] also analyzed the performance of A&P to estimate Kcb when fc and h are obtained from remote sensing. This approach is used in California’s SIMS together with ETo from CIMIS to compute the ETc for irrigation advising, resulting in a small mean bias error of 6.9% for trees and vines.
The methodology based on the dual Kc approach and the A&P approach shows a good potential for application to the coastal citrus orchards of Syria, namely using ground observations instead of remote sensing. The innovation required includes assigning characteristics of orchards, namely those relative to Kcb, to the various orchards. Once this is achieved, it is possible to simulate the water balance of the citrus depending on the characteristics and the Kcb of orchards, either in real time or relative to various water-demand scenarios. The companion paper by Darouich et al. [20] followed the research referred to above and consisted of the first part of the study on improved irrigation of citrus orchards to be continued through this article. Using the software SIMDualKc, Darouich et al. [20] developed a soil water balance of a clementine orchard, derived its crop coefficients (Kc) relative to various irrigation methods (drip, bubblers, micro-sprinkler, and ring basin), and assessed the performance of those methods, but focusing on drip irrigation positive impacts on the non-beneficial uses of water in citrus orchards, i.e., by decreasing evaporation losses during the dry summer season due to the small area wetted by drippers when compared with other methods. That study also showed the impact of canopy cover, tree height, and irrigation method on crop coefficients, and the impact was confirmed later by the standard Kcb of vineyards [21], as well as for citrus and various Mediterranean crops, using the published data of more than one hundred papers reporting good-quality field research [22]. Numerous applications of the A&P approach to tree crops have been reported, e.g., by Paço et al. [23] for peach; Paço et al. [24] and Puig-Sirera et al. [25] for olive; Vinci et al. [26] for hazelnut; and Mobe et al. [27] for various apple orchards in South Africa. Applications to vineyards include Williams et al. [21], who reported on a survey of California vineyards, using A&P. The applications to citrus are also numerous [28,29,30,31,32]. These applications make us confident about the effectiveness of using the A&P approach.
The current study, which serves as a companion study to that of Darouich et al. [20], aimed at further assessing appropriate issues to implement water-saving irrigation in the coastal Tartus and Latakia citrus orchards by focusing on irrigation scheduling in combination with drip irrigation. Considering the previous calibration and validation of SIMDualKc, as well as the careful definition of the orchards, water-saving thresholds, and atmospheric demand for water, the results of the simulations performed consist of effective information tools for irrigation scheduling thresholds to be used in practice. Moreover, the SIMDualKc is a well-proved SWB software model that, in addition to support deriving Kcb and Kc values, has been previously used to search for crop and irrigation practices that could control the water demand of maize and pea for industry [33,34]. In these applications, the model was first calibrated and validated using field data, the SWB was performed, and the debilities of the current irrigation management were identified, together with possible issues for improvement. Once the model was properly calibrated, simulations were performed to assess the top issues for irrigation and crop management, namely relative to irrigation scheduling and planting dates. Similar studies have been performed with tree crops after being adapted to the orchard and crop characteristics [25,35,36,37,38]. SIMDualKc applications include Rosa [29] for lemon and pear, Peddinti and Kambhammettu [39] for citrus, Cao et al. [40] for apple, Fandiño [41] for grapes, and Vinci et al. [26] for hazelnut. Naturally, despite the fact that the objective of improving agricultural water management is basically common, the orchard characteristics and research tools and practices used were different, thus implying innovative approaches for every application.
The main objective of this study was, therefore, to further pursue the work performed by Darouich et al. [20] by developing tools to improve the management of irrigation water in citrus orchards located in the Syrian coastal region. The specific objectives were (i) to derive Kcb and Kc values for the various citrus orchards, using the A&P approach; (ii) from those Kcb and Kc, to compute the crop transpiration (Tc), crop evapotranspiration (ETc), and irrigation water needs in Syrian citrus orchards, while considering different crop densities, water-shortage conditions, and climate variability (1988–2022); and (iii) to estimate yield reductions from adopting different water-saving irrigation strategies. The immediate results of this study contribute to overcoming the existing knowledge gap in the region in terms of rational data for effective use in irrigation water management, including in real time after the war ends. The ultimate objective is, therefore, the development and implementation of water-saving irrigation practices in the Syrian coastal region, contributing to the sustainability of local soil and water resources.

2. Materials and Methods

2.1. Study Area

The study area is the Tartus and Latakia provinces (Figure 1), which concentrate, respectively, 82% and 17% of the citrus production in Syria [42]. In all, 60% of the total production corresponds to orange (Citrus sinensis (L.) Osbeck), 12% to lemon (Citrus limon (L.) Osbeck), and 22% to clementine (Citrus clementina Hort.) and other species [4]. The climate in the region is hot-summer Mediterranean (Csa) [43]. Considering the period from 1998 to 2020, the mean annual air temperature is 19.3 °C. The mean monthly values vary from 11.5 °C in January to 27.0 °C in August. The mean annual precipitation is 930 mm and occurs mostly between October and May. The annual reference evapotranspiration (ETo) averages 1363 mm. The dominant soil reference groups are Vertisols, Cambisols, and Luvisols [44]. Irrigated agricultural land covers 29,100 ha in the Tartus district and 38,000 ha in the Latakia district [4]. Irrigated water is mainly withdrawn from surface water resources. The water-table depth varies from 10 to 20 m.

2.2. The A&P Approach

Allen and Pereira [17] proposed an empirical procedure for estimating Kcb values for the different stages of development of natural vegetation, tree orchards, and landscape systems based on the amount of vegetation present in those systems and background soil evaporation. In this approach, the transpiration component of the evapotranspiration process, represented by the Kcb, is related to the amount of vegetation of the crop density coefficient (Kd):
K c b = K c   m i n + K d K c b   f u l l K c   m i n
where Kd is the density coefficient that represents the impacts of plant density and/or leaf area, Kcb full is the estimated basal Kcb for plant-growth conditions having nearly full ground cover (or LAI > 3), and Kc min is the minimum Kc for bare soil (in the absence of vegetation).
When natural vegetation or grass covering the row and interrow of orchards is observed, Equation (1) is modified to account for the impact of those plants on the total evapotranspiration of the orchard system [17,46]:
K c b = K c b   c o v e r + K d   m a x K c b   f u l l K c b   c o v e r , K c b   f u l l K c b   c o v e r 2
where Kcb cover is the Kcb of the ground cover in the absence of tree foliage, Kd is the density coefficient, and Kcb full is the basal Kcb anticipated for the crop under full-cover conditions and corrected for climate. The second term of the max function, which accounts for the effects of shading by the active ground cover, reduces the estimated Kcb by half the difference between Kcb full and Kcb cover when this difference is negative. The value for Kcb cover in Equation (2) should represent the Kcb of the surface cover in the absence of tree cover; it should therefore reflect the density and vigor of the surface cover, as in areas exposed to sunlight.
The density coefficient Kd is estimated from observations of the fraction of the ground covered by vegetation (fc) and plant height (h) and describes the increase in Kc with increases in the amount of vegetation. Kd is estimated as follows [18]:
K d = m i n 1 , M L   f c   e f f , f c   e f f 1 1 + h
where fc eff is the effective fraction of the ground covered or shaded by vegetation (0.01–1) near solar noon, ML is a multiplier on fc eff describing the effect of canopy density on shading and maximum relative ET per fraction of shaded ground [1.0–2.0], and h is the mean height of vegetation (m).
The Kcb full value represents the upper limit of the system and is estimated as follows:
K c b   f u l l = F r m i n 1.0 + k h   h , 1.20 + 0.04 u 2 2 0.004 R H m i n 45 h 3 0.3
where u2 is the mean daily wind speed at a 2 m height (m s−1) during the growth period, RHmin (%) is the mean daily minimum relative humidity during the growth period, and h is the mean plant height (m) during mid-season. Before climatic adjustment, the upper limit for Kcb full is 1.20. The effect of crop height is considered through the sum (1 + kh h), with kh = 0.1 for tree and vine crops [18]. Higher Kcb full values are expected for taller crops and when the local climate is drier and/or windier than the standard climate conditions (RHmin = 45% and u2 = 2 m s−1) adopted in FAO 56. When the vegetation shows greater stomatal adjustment upon transpiration, parameter Fr applies an empirical adjustment (Fr < 1.0); otherwise, Fr = 1.0. For trees and vines, Fr is closer to 1.0 when crops exhibit great vegetative vigor; Fr decreases with limited water supply and due to pruning and training when the crop is stressed, and stomatal adjustment occurs. The Fr is defined as follows [18]:
F r = + γ 1 + 0.34   u 2 + γ 1 + 0.34   u 2 r l r t y p
where rl and rtyp are, respectively, the estimated actual mean leaf resistance and the typical leaf resistance (s m−1) for the vegetation in question. Moreover, Δ is the slope of the saturation vapor pressure vs. air temperature curve (kPa °C−1), and γ is the psychrometric constant (kPa °C−1), both relative to the period when Kcb full is computed. When standard Kcb values are considered, e.g., as initial values of Kcb for calibration purposes, Fr = 1.0 is assumed. Differently, when searching for actual Kcb values, Fr < 1.0 may be estimated with support of the tabulated values in Pereira et al. [18].
The application of this approach does not require calibration/validation if using the tabulated parameters in Pereira et al. [19]. Nevertheless, when field data are available, a validation may be performed by comparing the Kcb values derived from the A&P approach with those computed from field data, as already performed for the study area by Darouich et al. [20].

2.3. Management Scenarios

The soil water balance was computed for different orchard scenarios related to planting density, climate demand, and irrigation strategies. For each scenario, the crop evapotranspiration was computed by combining the A&P approach for assessment of the transpiration component according to scenario characteristics, and simulations of the SIMDualKc model for assessment of the soil evaporation component. Both methodologies were previously calibrated/validated in the companion study [20], therefore providing reliable estimates for both components in the study area.
Citrus orchards in the Syrian coastal area were divided into three major groups based on flowering and harvesting dates:
  • Group 1 (G1) includes clementine and mandarin species, with flowering by March and harvesting from September to October [20].
  • Group 2 (G2) refers to sweet orange, with flowering in May and the harvest from November to January [47,48].
  • Group 3 (G3) includes lemon and lime species. Three flowering periods were considered (April, May, and November), and the harvest was throughout the year [29,49].
Each group was divided into 5 subgroups (Table 1), accounting for the density of citrus orchards, which were characterized based on representative fc and h values observed in the region, as well as the Fr values tabulated in Pereira et al. [19]. This information was used to estimate the Kcb values for each orchard following the A&P approach.
The net irrigation water requirements were estimated by computing the daily soil water balance using the SIMDualKc model [50], which was used in the companion paper by Darouich et al. [20]. The model adopts the FAO56 dual Kc approach, thus partitioning evapotranspiration fluxes into its components, namely crop transpiration (Tc) and soil evaporation (Es) [14,51]. The SIMDualKc model has proved to be quite precise for computing crop evapotranspiration and performing the water balance to assess crop irrigation requirements, as analyzed by Pereira et al. [52] and referred to in the Introduction. The model was calibrated and validated for citrus, as reported in the companion paper (Darouich et al. [20]), where a complete description of the model and of its performance is provided. Therefore, because the companion paper is open access, a full description of the model and of its application is not given here.
In the current study, while the Kcb values were provided by the A&P approach, allowing the direct estimate of the Tc component, the Es component required the daily computation of the soil water balance of the soil evaporative layer, which was performed with the SIMDualKc model. The parametrization of irrigation followed observations reported in the companion paper [20]. Drip irrigation was used. The wetted fraction (fw) was small, up to 0.25. In all scenarios, the net irrigation depths were fixed at 5 mm per event, except for the case of high-density and tall trees, for which the application depths were fixed to 10 mm to compensate for the higher irrigation demand.
Simulations were performed for the period 1998–2022 (25 years). Daily meteorological data were taken from the Zahid station (34°41′37″ N, 35°59′16″ E; 12 m a.s.l.), which was considered to be representative of the Syrian coastal area. The data included the minimum and maximum air temperatures (Tmin and Tmax; °C), precipitation (P; mm), minimum and maximum relative humidity (RHmin and RHmax; %), number of sunshine hours (Isun; h), and wind speed at 2 m height (u2; m s−1). As in Darouich et al. [8], missing u2 values in 2011–2013 and 2015 were filled with wind data from the closest weather station (Trípoli, Lebanon), following the recommendations in FAO56. Figure 2 presents a brief characterization of the weather data used, as well as the annual precipitation and ETo variability throughout the study period.
To build the climate scenarios, data were then ordered from the lowest to the highest based on atmospheric demand (annual ETo). The years with the probability of occurrence of 20%, 50, 80%, and 95%, corresponding to low (2012), medium (2011), high (2010), and very high (2022) atmospheric demand, were then selected for the analysis.
The dates of the crop stages were defined based on the growing-degree days reported in the companion paper [20]. The dates of the initial stage were set following an analysis of the daily air temperatures at the beginning of the year. Table 2 presents the dates of the crop stages for every citrus orchard and climate scenario considered in this study. The parametrization of soil properties also followed observations reported in the companion paper [20]. The soil at the Zahid station was assumed to represent the study region, namely relative to the total and readily available water (TAW and RAW; mm). The soil water balance simulations were performed using the calibrated parameters of runoff (CN) and deep percolation functions (aD and bD), the soil water depletion fraction for no stress (p), the total and readily evaporable water (TEW and REW; mm), and the thickness of the soil evaporation layer (Ze; mm). The model parameters are given in Table 3.
Deficit-irrigation/water-saving scenarios were defined to improve irrigation water use while further considering the expected increasing scarcity of water in the region. Irrigation schedules were defined by imposing a water deficit of 5%, 10%, 20%, and 30% of p, the soil water depletion fraction for no stress. These values define the irrigation trigger thresholds, i.e., the management-allowed depletion (MAD) = 1.05 p, 1.10 p, 1.20 p, and 1.30 p, respectively.

2.4. Evaluation of Yield Decline in Relation to Deficit-Irrigation Management

For each scenario, the linear crop water-yield function (Doorenbos and Kassam [53]) was used to assess the impact of the proposed deficit-irrigation schedules on crop yields, as follows:
1 Y a Y m = K y 1 E T c   a c t E T c
where Ya and Ym are the actual and maximum crop yields, respectively; ETc act and ETc are the corresponding actual and potential seasonal crop evapotranspiration, respectively; and Ky is the yield response factor describing the reduction in relative yield due to the reduction in ET caused by the soil water shortage. The Ky values were defined for the mid-season stage. For the citrus species in G1 (clementine and mandarin) and G2 (sweet orange), Ky = 1.1 was assumed, considering their relatively high sensitivity to water stress [53,54,55]. The species in G3 (lemon) were assumed to exhibit a greater tolerance to water stress; thus, Ky = 0.9 was adopted [29,49]. Following Manssur [42], the Ym values were set to 30, 40, and 30 tons ha−1 for the crops in G1, G2, and G3, respectively. However, it is well-known that yields change with variety, crop management, environmental conditions, irrigation management, pests and diseases, and other factors.

3. Results and Discussion

3.1. Crop Coefficients for Citrus Orchards

The basal crop coefficients (Kcb) computed for the different crop stages using the A&P approach are shown in Table 4.
A single average value is presented for each stage following the single Fr values tabulated in Pereira et al. [19], which are not different among the citrus species; further field studies will provide different Fr values. Table 4 further includes the mean crop coefficient (Kc) for the different crop stages, thus also including the evaporation coefficient (Ke) derived from the computation of the daily soil water balance using the SIMDualKc model.
Estimates of the Kcb values provided by the A&P approach were already validated in Darouich et al. [20], returning the same or very similar values as those estimated from changes in soil water storage and the SIMDualKc model. As expected, Kcb values were very similar between the different crop stages and increased with the increase of tree density, i.e., with higher fc and h values. In Darouich et al. [20], the relationship between the Kcb and fc and h had already been revealed, with 10–14-year-old citrus and lower fc and h values returning lower Kcb values than 18–20-year-old trees having larger fc and h values. The Kcb values also changed with climate scenarios, from low to very high evapotranspiration demand, namely considering the inverse relation between Kcb full and RHmin (Equation (4)). The Kcb values for scenarios with a high-climate and very high climate demand were in close agreement with those tabulated in Pereira et al. [19] and obtained using the A&P approach. For low- and medium-climate demand scenarios, the estimated Kcb values in Table 4 were inferior to the tabulated ones in Pereira et al. [19].
The Kcb estimates for low-density orchards agreed with those for orange in Vilallobos et al. [56], for which an fc of 0.42 corresponded to Kcb mid = 0.58; and with El-Raki et al. [57], for which an fc of 0.30 corresponded to Kcb mid = 0.50. They were higher than those for orange, clementine, and mandarin in Ramos et al. [58] because of the lower tree height in these orchards with 5–6-year-old trees. Inversely, they were lower than those for clementine in Ballester et al. [59], for which an fc of 0.38 resulted in a Kcb mid = 0.65. Generally, the present results agreed with those reported by Pereira et al. [22]. The same happened with the medium-density orchards, for which the estimated Kcb values corresponded with those for clementine in Maestre-Valero et al. [60], with fc = 0.66 and a Kcb mid ranging from 0.50 to 0.65; and for orange reported by Taylor et al. [61], with an fc from 0.60 to 0.63 and a Kcb mid from 0.55 to 0.65. Relative to the Kcb estimates for high-density orchards, the values also agreed with those for lemon in Rosa [29], with fc = 0.75 and Kcb = 0.67 for the mid- and end-season; and for orange in Taylor et al. [60], Peddinti and Kambhammettu [39], Jamshidi et al. [62], and Jafari et al. [63], with the fc in orchards ranging from 0.70 to 0.88 and the Kcb mid from 0.70 to 0.85.
The Kc values in Table 4 showed similar dynamics as those already reported in Darouich et al. [20], i.e., lower values during the dry irrigation season because of the small wetted area below drippers, and higher values during the rainy season when the entire soil surface was wetted and subjected to soil evaporation. Again, the reported Kc values find agreement with several studies performed in Mediterranean citrus orchards [57,64,65], as well as with those now reported by Pereira et al. [22]. As such, the Kcb and Kc values in Table 4 are found to be adequate for estimating the water requirements of Syrian citrus orchards.

3.2. Consumptive Use

The consumptive use of citrus species in G1 (clementine and mandarin), G2 (sweet orange), and G3 (lemon and lime) for climate scenarios of low, average, high, and very high demand and for the water-saving/deficit-irrigation scenarios relative to the threshold MAD of 1.05, 1.10, 1.20, and 1.30 is presented in Table 5, Table 6 and Table 7.
The consumptive use represents the water taken up by plants in the evapotranspiration process (use) and that evaporated from the soil and transpired by the plants, i.e., water that is not reusable. Because Kc values were assumed to be the same among citrus species (Table 4), as with FAO56 and in Rallo et al. [15], the potential consumptive use of citrus trees ended quite similarly among all the groups, i.e., G1, G2, and G3, with minor differences resulting from different dates and lengths of the crop-growth stages (Table 2).
The most important factors affecting the potential consumptive use (ETc) were the density of trees in the citrus orchards and the tree heights. The larger the trees’ canopy, the greater the consumptive use, resulting also in higher Kcb values being estimated with the A&P approach. Large differences were naturally found among climate scenarios, with seasonal ETc values ranging from 867 mm under low demand to 1573 mm under very high climatic demand. The beneficial water use, i.e., the Tc component, represented 65.4% to 84.3% of the total consumption use and was higher for the climate scenarios of very high demand because the transpiration surface of trees was larger then, while the evaporative soil surface (few) was smaller since it was reduced to the exposed area wetted by drippers.
Darouich et al. [20] have already shown that drip irrigation could reduce the non-beneficial use of water, i.e., the Es component, when compared with other irrigation methods (bubblers, micro-sprinkler, and ring basin). Jovanovic et al. [66] listed several other strategies to further reduce the non-beneficial use of water in orchard systems, mainly deficit-irrigation strategies. Figure 3 presents the consumption-use savings (ETc–ETc act) in G1, G2, and G3, resulting from deficit irrigation. The consumption-use savings increased with climate demand and with deficit irrigation because less water was applied, thus reducing the ETc act values. Naturally, such strategies impacted yields, as discussed in Section 3.4.

3.3. Gross Irrigation Water Savings

Figure 4 presents the gross irrigation water savings in G1, G2, and G3, considering an irrigation efficiency of drip systems of 85% [67]. A full irrigation schedule was considered to be the baseline to quantify water savings in each deficit-irrigation scenario. This full irrigation schedule aimed to maintain soil moisture levels close to the RAW, without the crop entering stress. Therefore, the values shown in Figure 4 represent a minimum threshold of possible savings because substantially higher values would be obtained if the target of the full irrigation schedule was to maintain soil moisture closer to the soil field capacity.
Estimated gross (and net) crop irrigation requirements (CIRs) found a general correspondence with values in the literature, particularly when deficit-irrigation strategies were considered [58,68,69]. In all groups, the CIR increased with plant density and atmospheric demand. Additional water was needed to maintain soil moisture conditions within the proposed trigger and target thresholds of each irrigation strategy. Water savings followed the same trend, with higher values also found with increasing planting density and climate demand. Still, the most important factor contributing to saving water was deficit irrigation. This resulted in substantial savings of water, which ranged from 65 mm (low-density orchards under low demand) to 376 mm (high-density orchards and very high demand). The higher the imposed stress, the lower the amount of irrigation water applied and the larger the water saving when compared with the full irrigation schedule.

3.4. Expected Yield Decline

The expected yield declines due to the estimated deficits produced for water saving are presented in Table 8, Table 9 and Table 10 for G1, G2, and G3 crops. Yield reductions were similar in G1 and G2. G3 showed more contrasting differences due to the greater tolerance of lemon and lime to water stress.
Plant density and atmospheric demand had little influence on relative yields. On the other hand, the greater the water deficit condition imposed by the irrigation strategy, the greater the yield reduction. Considering that some of the estimated yield reduction, namely for irrigation deficit scenarios of 20% and 30%, represent a significant impact on farmers’ income and profitability, an economic analysis needs to be performed to complement the viability of those strategies. Still, results need also to be interpreted in a context of increasing scarcity that requires the implementation of policy measures to help mitigate farmers’ losses.

3.5. Improving Water Management of Citrus Orchards in Coastal Syria

The citrus sector in Syria is slowly recovering from years of civil war and embargo to economic trade. Once one of the most important sources of income in the agricultural sector, Syrian growers face enormous challenges today because of the economic crisis and limits to exportations that have caused the drop of fruit prices and increased costs of inputs such as fertilizers. At the same time, drought and increasing temperatures have built up a larger pressure on water resources, stressing the need for improved management.
This study contributes the necessary information for improving irrigation water management in Syrian citrus orchards in order to build a software platform to progressively implement farmers’ advice on the most adequate water savings, presently the best MAD irrigation scheduling thresholds. This can not only contribute to reducing water withdrawals in the region by improving crop water use and irrigation efficiency but also by increasing water and land productivity and further minimizing nutrient leaching from agricultural plots. On the one hand, Kcb and Kc values are provided to be efficiently used in irrigation management (Table 4), allowing us to accurately estimate crop water requirements in orchards with different planting densities and under different climate-demand conditions. On the other hand, MAD thresholds for the consumptive use of citrus trees and respective irrigation requirements are defined by considering the planting densities and climate conditions, as well as the deficit-irrigation/yield-decrease options. These deficit-irrigation strategies, which consisted of adopting an increased MAD of 1.05, 1.10, 1.20, and 1.30 p, shall be optimized when yield values will not be submitted to the war pressure and when production costs result from an open market and water costs will reflect an agricultural policy and not a war policy. Nevertheless, the implementation of this methodology is appropriate to help build an information platform and to help farmers improve their crop water use under scarcity.
The seasonal consumptive use of the different citrus species considered in the analysis ranged from 867 to 869 mm in orchards of low planting density under low demand and from 1549 to 1573 mm in orchards of very high planting density under very high climatic demand (Table 5, Table 6 and Table 7). In a full irrigation schedule, about 44% (425 mm) to 66% (820 mm) of the seasonal consumption use is likely to be met by irrigation. The estimated values are comparable to those in the literature, although the highest are not commonly observed in the region, and we hope that they will not be often achieved. The related gross irrigation requirements were estimated to be 15% higher than those net estimates, representing substantial water extractions that may not be available in a context of increasing scarcity. For this reason, several deficit-irrigation strategies were evaluated to reduce the seasonal consumption of water and irrigation needs of citrus orchards.
The mild deficit-irrigation strategy (MAD = 1.05 p) led to an almost nil reduction in the consumptive use of water (1–4%) but to important savings in the gross CIR (8–14%), with quite a small impact on yields (5–9% reduction). As deficit-irrigation strategies become more pronounced, ETc act may reduce by less than 10–16% (e.g., MAD = 1.10 p); thus, a small reduction in the consumptive use of water was also foreseen. However, when the gross water savings were high (30–37%), so were the impacts on yields (30–50% reduction), and their economic feasibility is then doubtful. Thus, moderate-to-high water-saving strategies may only be viable when supported by an economic analysis of the trade-offs between saved water and monetary incomes. This seems even more relevant in today’s context of low prices for marketable yields. Based on this analysis, it seems reasonable to conclude that substantial water savings can only be reached when the prices of citrus fruits recover to the levels prior to the civil war and the economic crisis, and the market of goods for crop production and water application and management will not be constrained as at present.

4. Conclusions

This study has to be considered together with the one reported in a companion paper which showed a clear need for the improvement of the irrigation water management of citrus orchards in the Syrian coastal area. The current study created data required to progressively build up an information platform that offers growers’ advice on the irrigation management and water saving required for the increased impacts of climate change, namely relative to higher water scarcity and higher temperature. The representative Kcb and Kc values for the irrigation management of citrus orchards in the region were proved, with Kcb values computed from representative fc and h by using the previously validated A&P approach. The Ke values were obtained through the computation of the soil water balance, using the SIMDualKc model, which was also previously calibrated and validated in the study area, using the observed data of soil moisture dynamics collected during eight growing seasons. Both the derived Kcb and Kc values found correspondence with those determined in Mediterranean climate orchards with similar conditions.
The consumptive use of citrus orchards was determined on the basis of a variety of characteristics observed in the Syrian coastal area, mainly tree density and plant height. Estimates further considered the variability of climate conditions in the Syrian coastal area, with water-saving thresholds set for scenarios of low, medium, high, and very high demand. The gross irrigation water savings ranged from 65 to 376 mm, with natural impacts on yields. Relationships between water savings and expected yield declines were drawn. Estimates also considered deficit-irrigation strategies based on prospects of increasing scarcity. These information issues shall be progressively updated when the platform is installed, and research using the A&P approach and SIMDualKc shall continue to better identify the debilities that need to be focused on and also implement real-time irrigation management.

Author Contributions

Conceptualization, H.D. and L.S.P.; methodology, H.D.; software, H.D.; validation, R.K.; formal analysis, H.D.; investigation, H.D. and R.K.; writing—original draft preparation, H.D. and T.B.R.; writing—review and editing, T.B.R. and L.S.P.; supervision, L.S.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by national funds through the FCT, Fundação para a Ciência e a Tecnologia, I.P., under projects UIDB/04129/2020 of LEAF (Linking Landscape, Environment, Agriculture and Food Research Unit) and UIDP/EEA/50009/2020 of LARSyS. The support of the FCT through grants attributed to T. B. Ramos (CEECIND/01152/2017) and H. Darouich (CEECIND/01153/2017) is also acknowledged.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Food and Agriculture Organization of the United Nations. Statistics. Available online: https://www.fao.org/statistics/en/ (accessed on 21 January 2022).
  2. Westlake, M. The citrus sub-sector. Syrian agriculture at the crossroads. In Syrian Agriculture at the Crossroads, FAO Agricultural Policy and Economic Development Series No. 8; Fiorillo, C., Vercueil, J., Eds.; Food and Agriculture Organization of the United Nations: Rome, Italy, 2003; Chapter 8; pp. 193–216. [Google Scholar]
  3. Amoudi, M.; Ahmed, T. Using Factorial Analysis to Study the Factors Affecting the Production of Fruit Trees in the Syrian Coast. Syr. J. Agric. Res. 2021, 8, 232–245. [Google Scholar]
  4. Central Bureau of Statistics. Agriculture, Chapter 4, Tables 7 and 10, Damascus, Syria. Available online: http://cbssyr.sy/index-EN.htm (accessed on 28 December 2021).
  5. Oweis, T.; Rodrigues, P.N.; Pereira, L.S. Simulation of supplemental irrigation strategies for wheat in Near East to cope with water scarcity. In Tools for Drought Mitigation in Mediterranean Regions; Rossi, G., Cancelliere, A., Pereira, L.S., Oweis, T., Shatanawi, M., Zairi, A., Eds.; Kluwer: Dordrecht, The Netherlands, 2003; pp. 259–272. [Google Scholar] [CrossRef]
  6. Mourad, K.A.; Alshihabi, O. Assessment of future Syrian water resources supply and demand by the WEAP model. Hydrol. Sci. J. 2016, 61, 393–401. [Google Scholar] [CrossRef] [Green Version]
  7. Darouich, H.; Karfoul, R.; Eid, H.; Ramos, T.B.; Baddour, N.; Moustafa, A.; Assaad, M.I. Modeling zucchini squash irrigation requirements in the Syrian Akkar region using the FAO56 dual-Kc approach. Agric. Water Manag. 2020, 229, 105927. [Google Scholar] [CrossRef]
  8. Darouich, H.; Karfoul, R.; Ramos, T.B.; Moustafa, A.; Shaheen, B.; Pereira, L.S. Crop water requirements and crop coefficients for jute mallow (Corchorus olitorius L.) using the SIMDualKc model and assessing irrigation strategies for the Syrian Akkar region. Agric. Water Manag. 2021, 255, 107038. [Google Scholar] [CrossRef]
  9. Abou Zakhem, B.; Al Ain, F.; Hafez, R. Assessment of field water budget components for increasing water productivity under drip irrigation in arid and semiarid areas, Syria. Irrig. Drain. 2019, 68, 452–463. [Google Scholar] [CrossRef]
  10. Fader, M.; Shi, S.; von Bloh, W.; Bondeau, A.; Cramer, W. Mediterranean irrigation under climate change: More efficient irrigation needed to compensate increases in irrigation water requirements. Hydrol. Earth Syst. Sci. 2016, 20, 953–973. [Google Scholar] [CrossRef] [Green Version]
  11. Darouich, H.; Gonçalves, J.M.; Muga, A.; Pereira, L.S. Water saving vs. farm economics in cotton surface irrigation: An application of multicriteria analysis. Agric. Water Manag. 2012, 115, 223–231. [Google Scholar] [CrossRef]
  12. Darouich, H.; Pedras, C.M.G.; Gonçalves, J.M.; Pereira, L.S. Drip vs. surface irrigation: A comparison focusing water saving and economic returns using multicriteria analysis applied to cotton. Biosyst. Eng. 2014, 122, 74–90. [Google Scholar] [CrossRef]
  13. Darouich, H.; Cameira, R.M.; Gonçalves, J.M.; Paredes, P.; Pereira, L.S. Comparing sprinkler and surface irrigation for wheat using multi-criteria analysis: Water saving vs. economic returns. Water 2017, 9, 50. [Google Scholar] [CrossRef] [Green Version]
  14. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration—Guidelines for Computing Crop Water Requirements; Irrigation & Drainage Paper 56; Food and Agriculture Organization of the United Nations (FAO): Rome, Italy, 1998. [Google Scholar]
  15. Rallo, G.; Paço, T.A.; Paredes, P.; Puig-Sirera, À.; Massai, R.; Provenzano, G.; Pereira, L.S. Updated single and dual crop coefficients for tree and vine fruit crops. Agric. Water Manag. 2021, 250, 106645. [Google Scholar] [CrossRef]
  16. López-Urrea, R.; Montoro, A.; Mañas, F.; López-Fuster, P.; Fereres, E. Evapotranspiration and crop coefficients from lysimeter measurements of mature Tempranillo wine grapes. Agric Water Manag. 2012, 112, 13–20. [Google Scholar] [CrossRef]
  17. Allen, R.G.; Pereira, L.S. Estimating crop coefficients from fraction of ground cover and height. Irrig. Sci. 2009, 28, 17–34. [Google Scholar] [CrossRef] [Green Version]
  18. Pereira, L.S.; Paredes, P.; Melton, F.; Johnson, L.; Wang, T.; López-Urrea, R.; Cancela, J.J.; Allen, R. Prediction of crop coefficients from fraction of ground cover and height. Background and validation using ground and remote sensing data. Agric. Water Manag. 2020, 240, 106197. [Google Scholar] [CrossRef]
  19. Pereira, L.S.; Paredes, P.; Melton, F.; Johnson, L.; Mota, M.; Wang, T. Prediction of crop coefficients from fraction of ground cover and height: Practical application to vegetable, field, and fruit crops with focus on parameterization. Agric. Water Manag. 2021, 252, 106663. [Google Scholar] [CrossRef]
  20. Darouich, H.; Karfoul, R.; Ramos, T.B.; Moustafa, A.; Pereira, L.S. Searching for Sustainable-Irrigation Issues of Clementine Orchards in the Syrian Akkar Plain: Effects of Irrigation Method and Canopy Size on Crop Coefficients, Transpiration, and Water Use with SIMDualKc Model. Water 2022, 14, 2052. [Google Scholar] [CrossRef]
  21. Williams, L.E.; Levin, A.D.; Fidelibus, M.W. Crop coefficients (Kc) developed from canopy shaded area in California vineyards. Agric. Water Manag. 2022, 271, 107771. [Google Scholar] [CrossRef]
  22. Pereira, L.S.; Paredes, P.; Oliveira, C.M.; Montoya, F.; López-Urrea, R.; Salman, M. Single and basal crop coefficients for estimation of water use of tree and vine woody crops with consideration of fraction of ground cover, height, and training system for Mediterranean and warm temperate fruit and leaf crops. Irrig. Sci. 2023; submitted. [Google Scholar]
  23. Paço, T.; Ferreira, M.; Rosa, R.; Paredes, P.; Rodrigues, G.; Conceição, N.; Pacheco, C.; Pereira, L. The dual crop coefficient approach using a density factor to simulate the evapotranspiration of a peach orchard: SIMDualKc model versus eddy covariance measurements. Irrig. Sci. 2012, 30, 115–126. [Google Scholar] [CrossRef]
  24. Paço, T.A.; Paredes, P.; Pereira, L.S.; Silvestre, J.; Santos, F.L. Crop coefficients and transpiration of a super intensive Arbequina olive orchard using the dual Kc approach and the Kcb computation with the fraction of ground cover and height. Water 2019, 11, 383. [Google Scholar] [CrossRef] [Green Version]
  25. Puig-Sirera, À.; Rallo, G.; Paredes, P.; Paço, T.A.; Minacapilli, M.; Provenzano, G.; Pereira, L.S. Transpiration and water use of an irrigated traditional olive grove with sap-flow observations and the FAO56 dual crop coefficient approach. Water 2021, 13, 2466. [Google Scholar] [CrossRef]
  26. Vinci, A.; Traini, C.; Portarena, S.; Farinelli, D. Assessment of the Midseason Crop Coefficient for the Evaluation of the Water Demand of Young, Grafted Hazelnut Trees in High-Density Orchards. Water 2023, 15, 1683. [Google Scholar] [CrossRef]
  27. Mobe, N.T.; Dzikiti, S.; Zirebwa, S.F.; Midgley, S.J.E.; von Loeper, W.; Mazvimavi, D.; Ntshidi, Z.; Jovanovic, N.Z. Estimating crop coefficients for apple orchards with varying canopy cover using measured data from twelve orchards in the Western Cape Province, South Africa. Agric. Water Manag. 2020, 233, 106103. [Google Scholar] [CrossRef]
  28. Taylor, N.J.; Annandale, J.G.; Vahrmeijer, J.T.; Ibraimo, N.A.; Mahohoma, W.; Gush, M.B.; Allen, R.G. Modelling water use of subtropical fruit crops: The challenges. Acta Hortic. 2017, 1160, 277–284. [Google Scholar] [CrossRef] [Green Version]
  29. Rosa, R.D. Modelação da Evapotranspiração com o Modelo SIMDualKc: Aplicação à Rega de Fruteiras, a Consociações de Culturas e a Condições Salinas, e Ligação ao SIG para Análise à Escala do Projeto de Rega. Ph.D. Thesis, Instituto Superior de Agronomia, Lisboa, Portugal, 2018. [Google Scholar]
  30. Fernández Hung, K.; Vargas Rodríguez, P.; Cueto Rodríguez, J.R.; Brown Manrique, O. Water requirements of new citrus orchards in “Jiguaní” Agricultural Enterprise. INGE CUC 2022, 18, 95–104. [Google Scholar]
  31. Ippolito, M.; De Caro, D.; Ciraolo, G.; Minacapilli, M.; Provenzano, G. Estimating crop coefficients and actual evapotranspiration in citrus orchards with sporadic cover weeds based on ground and remote sensing data. Irrig. Sci. 2023, 41, 5–22. [Google Scholar] [CrossRef]
  32. Mashabatu, M.; Ntshidi, Z.; Dzikiti, S.; Jovanovic, N.; Dube, T.; Taylor, N.J. Deriving crop coefficients for evergreen and deciduous fruit orchards in South Africa using the fraction of vegetation cover and tree height data. Agric. Water Manag. 2023, 286, 108389. [Google Scholar] [CrossRef]
  33. Paredes, P.; Rodrigues, G.C.; Alves, I.; Pereira, L.S. Partitioning evapotranspiration, yield prediction and economic returns of maize under various irrigation management strategies. Agric. Water Manag. 2014, 135, 27–39. [Google Scholar] [CrossRef]
  34. Paredes, P.; Pereira, L.S.; Rodrigues, G.C.; Botelho, N.; Torres, M.O. Using the FAO dual crop coefficient approach to model water use and productivity of processing pea (Pisum sativum L.) as influenced by irrigation strategies. Agric. Water Manag. 2017, 189, 5–18. [Google Scholar] [CrossRef]
  35. Phogat, V.; Skewes, M.A.; Cox, J.W.; Sanderson, G.; Alam, J.; Šimunek, J. Seasonal simulation of water, salinity and nitrate dynamics under drip irrigated mandarin (Citrus reticulata) and assessing management options for drainage and nitrate leaching. J. Hydrol. 2014, 513, 504–516. [Google Scholar] [CrossRef] [Green Version]
  36. Li, D.; Hendricks Franssen, H.-J.; Han, X.; Jiménez-Bello, M.A.; Martínez Alzamora, F.; Vereecken, H. Evaluation of an operational real-time irrigation scheduling scheme for drip irrigated citrus fields in Picassent, Spain. Agric. Water Manag. 2018, 208, 465–477. [Google Scholar] [CrossRef] [Green Version]
  37. Abou Ali, A.; Bouchaou, L.; Er-Raki, S.; Hssaissoune, M.; Brouziyne, Y.; Ezzahar, J.; Khabba, S.; Chakir, A.; Labbaci, A.; Chehbouni, A. Assessment of crop evapotranspiration and deep percolation in a commercial irrigated citrus orchard under semi-arid climate: Combined Eddy-Covariance measurement and soil water balance-based approach. Agric. Water Manag. 2023, 275, 107997. [Google Scholar] [CrossRef]
  38. Zhu, S.; Li, H.; Xu, W.; Zhang, T.; Liu, C.; Peng, Y. Dynamic process and model simulation of soil water content during citrus growth period. Acta Geophys. 2023, 71, 1525–1537. [Google Scholar] [CrossRef]
  39. Peddinti, S.R.; Kambhammettu, B.V.N.P. Dynamics of crop coefficients for citrus orchards of central India using water balance and eddy covariance flux partition techniques. Agric. Water Manag. 2019, 212, 68–77. [Google Scholar] [CrossRef]
  40. Cao, H.; Wang, H.; Li, Y.; Hamani, A.K.M.; Zhang, N.; Wang, X.; Gao, Y. Evapotranspiration Partition and Dual Crop Coefficients in Apple Orchard with Dwarf Stocks and Dense Planting in Arid Region, Aksu Oasis, Southern Xinjiang. Agriculture 2021, 11, 1167. [Google Scholar] [CrossRef]
  41. Fandiño, M. Necesidades de Agua e Influencia de los Sistemas de Riego en Vitis vinifera cv. Albariño. Ph.D. Thesis, Escuela de Doctorado Internacional, Universidad de Santiago de Compostela, Lugo, Spain, 2021. [Google Scholar]
  42. Manssur, F. The Agricultural Commodity Perspective n°5: The Citrus in Syria; FAO Project GCP/SYR/006/ITA; The National Agricultural Policy Center (NAPC): Damascus, Syria, 2007. [Google Scholar]
  43. Köppen, W. Die Wärmezonen der Erde, nach der Dauer der heissen, gemässigten und kalten Zeit und nach der Wirkung der Wärme auf die organische Welt betrachtet. Meteorol. Z. 1884, 1, 215–226. [Google Scholar]
  44. FAO; IIASA; ISRIC; ISS-CAS; JRC. Harmonized World Soil Database (Version 1.1); FAO: Rome, Italy; IIASA: Laxenburg, Austria, 2009. [Google Scholar]
  45. Mohamed, M.A. An assessment of forest cover change and its driving forces in the Syrian Coastal Region during a period of conflict, 2010 to 2020. Land 2021, 10, 191. [Google Scholar] [CrossRef]
  46. Allen, R.G.; Wright, J.L.; Pruitt, W.O.; Pereira, L.S.; Jensen, M.E. Water requirements. In Design and Operation of Farm Irrigation Systems, 2nd ed.; Hoffman, G.J., Evans, R.G., Jensen, M.E., Martin, D.L., Elliot, R.L., Eds.; ASABE: St. Joseph, MI, USA, 2007; pp. 208–288. [Google Scholar]
  47. Pérez-Pérez, J.G.; Romero, P.; Navarro, J.M.; Botía, P. Response of sweet orange cv ‘Lane late’to deficit-irrigation strategy in two rootstocks. II: Flowering, fruit growth, yield and fruit quality. Irrig. Sci. 2008, 26, 519–529. [Google Scholar] [CrossRef]
  48. Reykande, J.M.; Amiri, N.A.; Shahabian, M. Analyzing phenological stages of three citrus varieties at foothills, plain and shoreline areas of Sari in North of Iran. Int. J. Agric. Crop Sci. 2013, 6, 452–457. [Google Scholar]
  49. Shatanawi, M.R.; Suleiman, A.A.; Al-Bakri, J. Effect of deficit irrigation on lemon water requirements and yield in the Jordan Valley. Jordan J. Agric. Sci. 2011, 7, 564–573. [Google Scholar]
  50. Rosa, R.D.; Paredes, P.; Rodrigues, G.C.; Alves, I.; Fernando, R.M.; Pereira, L.S.; Allen, R.G. Implementing the dual crop coefficient approach in interactive software. 1. Background and computational strategy. Agric. Water Manag. 2012, 103, 8–24. [Google Scholar] [CrossRef]
  51. Allen, R.G.; Pereira, L.S.; Smith, M.; Raes, D.; Wright, J.L. FAO-56 dual crop coefficient method for estimating evaporation from soil and application extensions. J. Irrig. Drain. Eng. 2005, 131, 2–13. [Google Scholar] [CrossRef] [Green Version]
  52. Pereira, L.S.; Paredes, P.; Jovanovic, N. Soil water balance models for determining crop water and irrigation requirements and irrigation scheduling focusing on the FAO56 method and the dual Kc approach. Agric. Water Manag. 2020, 241, 106357. [Google Scholar] [CrossRef]
  53. Doorenbos, J.; Kassam, A.H. Yield Response to Water; Irrigation & Drainage Paper 33; Food and Agriculture Organization of the United Nations (FAO): Rome, Italy, 1979. [Google Scholar]
  54. González-Altozano, P.; Castel, J.R. Regulated deficit irrigation in ‘Clementina de Nules’ citrus trees. II: Vegetative growth. J. Hortic. Sci. Biotechnol. 2000, 75, 388–392. [Google Scholar] [CrossRef]
  55. García-Tejero, I.; Durán-Zuazo, V.H.; Arriaga-Sevilla, J.; Muriel Fernández, J.L. Impact of water stress on citrus yield. Agron. Sustain. Dev. 2012, 32, 651–659. [Google Scholar] [CrossRef] [Green Version]
  56. Villalobos, F.J.; Testi, L.; Orgaz, F.; García-Tejera, O.; Lopez-Bernal, A.; González-Dugo, M.V.; Ballester-Lurbe, C.; Castel, J.R.; Alarcón-Cabañero, J.J.; Nicolás-Nicolás, E.; et al. Modelling canopy conductance and transpiration of fruit trees in Mediterranean areas: A simplified approach. Agric. For. Meteorol. 2013, 171–172, 93–103. [Google Scholar] [CrossRef]
  57. Er-Raki, S.; Chehbouni, A.; Guemouria, N.; Ezzahar, J.; Khabba, S.; Boulet, G.; Hanich, L. Citrus orchard evapotranspiration: Comparison between eddy covariance measurements and the FAO-56 approach estimates. Plant Biosyst. 2009, 143, 201–208. [Google Scholar] [CrossRef] [Green Version]
  58. Ramos, T.B.; Darouich, H.; Oliveira, A.R.; Farzamian, M.; Monteiro, T.; Castanheira, N.; Paz, A.; Gonçalves, M.C.; Pereira, L.S. Water use and soil water balance of Mediterranean tree crops assessed with the SIMDualKc model in orchards of southern Portugal. Agric. Water Manag. 2023, 279, 108209. [Google Scholar] [CrossRef]
  59. Ballester, C.; Castel, J.; Introgliolo, D.S.; Castel, J.R. Response of Clementina de Nules citrus trees to summer deficit irrigation. Yield components and fruit composition. Agric. Water Manag. 2011, 98, 1027–1032. [Google Scholar] [CrossRef]
  60. Maestre-Valero, J.F.; Testi, L.; Jiménez-Bello, M.A.; Castel, J.R.; Intrigliolo, D.S. Evapotranspiration and carbon exchange in a citrus orchard using eddy covariance. Irrig Sci. 2017, 35, 397–408. [Google Scholar] [CrossRef]
  61. Taylor, N.J.; Mahohoma, W.; Vahrmeijer, J.T.; Gush, M.B.; Allen, R.G.; Annandale, J.G. Crop coefficient approaches based on fixed estimates of leaf resistance are not appropriate for estimating water use of citrus. Irrig. Sci. 2015, 33, 153–166. [Google Scholar] [CrossRef] [Green Version]
  62. Jamshidi, S.; Zand-Parsa, S.; Kamgar-Haghighi, A.A.; Shahsavar, A.R.; Niyogi, D. Evapotranspiration, crop coefficients, and physiological responses of citrus trees in semi-arid climatic conditions. Agric. Water Manag. 2020, 227, 105838. [Google Scholar] [CrossRef]
  63. Jafari, M.; Kamali, H.; Keshavarz, A.; Momeni, A. Estimation of evapotranspiration and crop coefficient of drip-irrigated orange trees under a semi-arid climate. Agric. Water Manag. 2021, 248, 106769. [Google Scholar] [CrossRef]
  64. Consoli, S.; Vanela, D. Mapping crop evapotranspiration by integrating vegetation indices into a soil water balance model. Agric. Water Manag. 2014, 143, 71–81. [Google Scholar] [CrossRef]
  65. Consoli, S.; O’Connell, N.; Snyder, R. Estimation of evapotranspiration of different-sized navel-orange tree orchards using energy balance. J. Irrig. Drain. Eng. 2006, 132, 2–8. [Google Scholar] [CrossRef]
  66. Jovanovic, N.; Pereira, L.S.; Paredes, P.; Pôças, I.; Cantore, V.; Todorovic, M. A review of strategies, methods and technologies to reduce non-beneficial consumptive water use on farms considering the FAO56 methods. Agric. Water Manag. 2020, 239, 106267. [Google Scholar] [CrossRef]
  67. Keller, J.; Bliesner, R.D. Sprinkle and Trickle Irrigation; Van Nostrand Reinhold: New York, NY, USA, 1990. [Google Scholar]
  68. Puig-Sirera, À.; Provenzano, G.; González-Altozano, P.; Intrigliolo, D.S.; Rallo, G. Irrigation water saving strategies in citrus orchards: Analysis of the combined effects of timing and severity of soil water deficit. Agric. Water Manag. 2021, 248, 106773. [Google Scholar] [CrossRef]
  69. Rallo, G.; González-Altozano, P.; Manzano-Juárez, J.; Provenzano, G. Using field measurements and FAO-56 model to assess the eco-physiological response of citrus orchards under regulated deficit irrigation. Agric. Water Manag. 2017, 180, 136–147. [Google Scholar] [CrossRef]
Figure 1. Location of the Tartus and Latakia provinces (adapted from Mohamed [45]).
Figure 1. Location of the Tartus and Latakia provinces (adapted from Mohamed [45]).
Agronomy 13 01794 g001
Figure 2. Above: monthly averages of maximum and minimum air temperatures (Tmax and Tmin; °C), maximum and minimum relative humidity (RHmax and RHmin; %), number of sunshine hours (Isun; h), and wind speed at a 2 m height (u2; m s−1). Below: annual precipitation (mm) and annual grass reference evapotranspiration (ETo; mm); all data are for the period 1998–2022.
Figure 2. Above: monthly averages of maximum and minimum air temperatures (Tmax and Tmin; °C), maximum and minimum relative humidity (RHmax and RHmin; %), number of sunshine hours (Isun; h), and wind speed at a 2 m height (u2; m s−1). Below: annual precipitation (mm) and annual grass reference evapotranspiration (ETo; mm); all data are for the period 1998–2022.
Agronomy 13 01794 g002
Figure 3. Consumptive-use savings in G1 (clementine and mandarin), G2 (orange), and G3 (lemon and lime) under different scenarios of tree density, climate demand, and deficit irrigation.
Figure 3. Consumptive-use savings in G1 (clementine and mandarin), G2 (orange), and G3 (lemon and lime) under different scenarios of tree density, climate demand, and deficit irrigation.
Agronomy 13 01794 g003
Figure 4. Gross irrigation water savings in G1 (clementine and mandarin), G2 (orange), and G3 (lemon and lime) under scenarios of different climate demand and deficit irrigation.
Figure 4. Gross irrigation water savings in G1 (clementine and mandarin), G2 (orange), and G3 (lemon and lime) under scenarios of different climate demand and deficit irrigation.
Agronomy 13 01794 g004
Table 1. Scenarios of plant density of citrus orchards.
Table 1. Scenarios of plant density of citrus orchards.
Orchardh
(m)
fc
(-)
Fr
(-)
Low density, tall tree4.50.400.63
Med density, small tree3.50.650.53
Med density, tall tree4.50.650.57
High density, small tree3.50.700.61
High density, tall tree4.50.700.78
Note: h, plant height; fc, the fraction of the ground covered by vegetation; Fr, empirical parameter.
Table 2. Dates of the crop stages for each citrus group (G1–G3) and climate-demand scenario (low, medium, high, and very high demand).
Table 2. Dates of the crop stages for each citrus group (G1–G3) and climate-demand scenario (low, medium, high, and very high demand).
Demand
Scenario
Non-
Growing
InitialCrop
Development
Mid-
Season
Late-
Season
End-
Season
G1:
Low01-Jan20-Feb14-Mar15-Jun05-Oct31-Dec
Medium01-Jan01-Feb15-Feb01-Jun01-Oct31-Dec
High-01-Jan31-Jan21-May01-Oct31-Dec
Very high-01-Jan22-Feb29-May19-Sep31-Dec
G2:
Low01-Jan20-Feb21-Mar21-Jun01-Nov31-Dec
Medium01-Jan19-Jan01-Mar15-Jun24-Oct31-Dec
High01-Jan15-Jan01-Mar15-Jun01-Nov31-Dec
Very high01-Jan06-Feb29-Mar19-Jun01-Nov31-Dec
G3:
Low01-Jan19-Feb31-Mar06-Jun15-Sep31-Dec
Medium01-Jan01-Mar22-Mar27-May15-Sep31-Dec
High01-Jan14-Feb08-Mar22-May05-Sep31-Dec
Very high01-Jan12-Feb04-Apr22-May05-Sep31-Dec
Table 3. Calibrated model parameters (from Darouich et al. [20]).
Table 3. Calibrated model parameters (from Darouich et al. [20]).
ParameterSymbolValue
Depletion fraction for no stresspini0.60
pmid0.60
pend0.60
Total evaporable waterTEW (mm)40
Readily evaporable waterREW (mm)8
Depth of the soil evaporation layerZe (m)0.10
Deep percolationaD490
bD−0.02
Runoff curve numberCN80
Table 4. Basal crop coefficients (Kcb) and mean crop coefficients (Kc) for orchards under different planting densities and climate scenarios (low, medium, high, and very high demand).
Table 4. Basal crop coefficients (Kcb) and mean crop coefficients (Kc) for orchards under different planting densities and climate scenarios (low, medium, high, and very high demand).
ScenarioInitialMid-SeasonLate Season
KcbKcKcbKcKcbKc
Low Density
Low0.501.120.500.700.491.09
Average0.501.150.520.720.501.12
High0.481.120.530.730.511.16
Very high0.541.250.560.780.541.21
Med density, small tree
Low0.550.950.520.670.530.91
Average0.550.970.550.710.550.97
High0.530.930.550.720.560.98
Very high0.591.040.590.770.591.03
Med density, tall tree
Low0.601.000.560.710.580.96
Average0.601.020.590.750.601.01
High0.570.970.600.770.611.04
Very high0.641.090.640.820.651.09
High density, small tree
Low0.640.990.610.750.620.95
Average0.640.990.630.780.640.98
High0.620.970.640.790.661.03
Very high0.681.060.680.860.691.06
High density, tall tree
Low0.841.170.790.930.811.11
Average0.841.160.820.970.841.16
High0.801.150.841.000.861.20
Very high0.901.270.891.060.911.25
Table 5. Season potential and actual crop evapotranspiration and transpiration (mm) for G1 crops (clementine and mandarin) under various scenarios of plant density, climate demand, and deficit-irrigation MAD * thresholds.
Table 5. Season potential and actual crop evapotranspiration and transpiration (mm) for G1 crops (clementine and mandarin) under various scenarios of plant density, climate demand, and deficit-irrigation MAD * thresholds.
Plant-Density
Scenarios
Climate-
Demand
Scenarios
ETc
(mm)
Tc
(mm)
MAD Irrigation Thresholds
1.051.101.201.30
ETc actTc actETc actTc actETc actTc actETc actTc act
LowLow869579856561841547809520778493
Average939614927597911583880556847528
High959648939624920605879570838533
Very high12298131203782117575711247101074665
Med, small treeLow815609800591785577754549724523
Average878655863637850625819597789570
High902682880657861639822603781565
Very high1157864112682911008051050758997709
Med, tall treeLow863658846637829622795590762560
Average929705912684896670863640829609
High963743938714916693872652826610
Very high12349401200902117087411128201055766
High, small treeLow900712879689860671823636786603
Average955753936731918714881680844646
High996795967763943740896696850652
Very high12659981223952119592611348691074811
High, tall treeLow1100926107890210548791000827952781
Average11719821149957111992910728841017832
High1235104112051009117197511079141039848
Very high1562131015171263147512211396114413091061
* ETc and ETc act are potential and actual crop evapotranspiration, respectively; Tc and Tc act are potential and actual crop transpiration, respectively; MAD, management-allowed deficit.
Table 6. Season potential and actual crop evapotranspiration and transpiration (mm) for G2 crops (orange) under various scenarios of plant density, climate demand, and deficit-irrigation MAD * thresholds.
Table 6. Season potential and actual crop evapotranspiration and transpiration (mm) for G2 crops (orange) under various scenarios of plant density, climate demand, and deficit-irrigation MAD * thresholds.
Plant-Density
Scenarios
Climate-
Demand
Scenarios
ETc
(mm)
Tc
(mm)
MAD Irrigation Thresholds
1.051.101.201.30
ETc actTc actETc actTc actETc actTc actETc actTc act
LowLow869579856560841547809520778493
Average938613925595910582878555846528
High953644933619915603873566833531
Very high12278121199778117575711227091072663
Med, small treeLow826619811601795586763558733531
Average878655863637850625819597789570
High898679877655858636818600779564
Very high1157864112682911008051050758997709
Med, tall treeLow874668857647838631805599771569
Average927703910683894667859637827607
High961741936713913692869650825609
Very high12359401199900117187511148201057767
High, small treeLow903715882691863673826639788604
Average965762943738926722889688853654
High991790962759939736892692845648
Very high127810111233963120693811458781083819
High, tall treeLow1103928107890110578821004831952780
Average11709811149956111992810688801017831
High122110271192995115796110938991027836
Very high1573132215271272148612331401115213171069
* ETc and ETc act are potential and actual crop evapotranspiration, respectively; Tc and Tc act are potential and actual crop transpiration, respectively; MAD, management-allowed deficit.
Table 7. Season potential and actual crop evapotranspiration and transpiration (mm) for G3 crops (lemon and limes) under various scenarios of plant density, climate demand, and deficit-irrigation MAD * thresholds.
Table 7. Season potential and actual crop evapotranspiration and transpiration (mm) for G3 crops (lemon and limes) under various scenarios of plant density, climate demand, and deficit-irrigation MAD * thresholds.
Plant-Density
Scenarios
Climate-
Demand
Scenarios
ETc
(mm)
Tc
(mm)
MAD Irrigation Thresholds
1.051.101.201.30
ETc actTc actETc actTc actETc actTc actETc actTc act
LowLow867576854558838545807518776491
Average938613926595910581880555846527
High956646937622917603876568836532
Very high12258091198777117375411227071070661
Med, small treeLow820614806596789581759554727526
Average878655863637850625819597789570
High909689887664867645827609787571
Very high1147855111882210947991041750990703
Med, tall treeLow871666854644836628802597767566
Average927703909682893667859636827607
High962741936713915693870651826610
Very high12229281187889116286511048111047758
High, small treeLow899711878687859670822635786602
Average965762943738926722889688853654
High1002801973769949745901700854656
Very high12649971222951119492511348681073810
High, tall treeLow1101926107890010548771003829951779
Average11789891155963112893810748861021836
High1231103811991002116496911029091035845
Very high1549129415041247146312071382112813021051
* ETc and ETc act are potential and actual crop evapotranspiration, respectively; Tc and Tc act are potential and actual crop transpiration, respectively; MAD, management-allowed deficit.
Table 8. Expected yield decline (YD) versus water saving (WS) for G1 (clementine and mandarin) under deficit-irrigation (management-allowed deficit, MAD: 1.05, 1.10, 1.20, and 1.30) and climate (low, medium, high, and very high demand) scenarios.
Table 8. Expected yield decline (YD) versus water saving (WS) for G1 (clementine and mandarin) under deficit-irrigation (management-allowed deficit, MAD: 1.05, 1.10, 1.20, and 1.30) and climate (low, medium, high, and very high demand) scenarios.
Plant DensityClimatic
Demand
Deficit-Irrigation Strategies
MAD = 1.05MAD = 1.10MAD = 1.20MAD = 1.30
WS
(%)
YD
(%)
WS
(%)
YD
(%)
WS
(%)
YD
(%)
WS
(%)
YD
(%)
LowLow147181228263742
Average146191128243838
High127171426293547
Very high118151524303150
Med, small treeLow136181228253742
Average146191129223835
High127171425283546
Very high118161524303349
Med, tall treeLow136181228253742
Average147191228243738
High118161425293448
Very high108151524303350
High, small treeLow127181227263642
Average137181227253740
High128161425293448
Very high119151524303350
High, tall treeLow85131022233138
Average95141223253343
High96121222253044
Very high87131321273145
Table 9. Expected yield decline (YD) versus water saving (WS) for G2 (orange) under deficit-irrigation (management allowed deficit, MAD: 1.05, 1.10, 1.20, and 1.30) and climate (low, medium, high, and very high demand) scenarios.
Table 9. Expected yield decline (YD) versus water saving (WS) for G2 (orange) under deficit-irrigation (management allowed deficit, MAD: 1.05, 1.10, 1.20, and 1.30) and climate (low, medium, high, and very high demand) scenarios.
Plant DensityDemandDeficit-Irrigation Strategies
5%10%20%30%
WS
(%)
YD
(%)
WS
(%)
YD
(%)
WS
(%)
YD
(%)
WS
(%)
YD
(%)
Low Low147181228263742
Average146181227243738
High128161425293446
Very high109141423303249
Med, small treeLow136181127223735
Average146191129233837
High126171226253441
Very high118161424293348
Med, tall treeLow136181127223636
Average146191229253740
High117161225253442
Very high108141423293348
High, small treeLow126181127223636
Average147181228253741
High117161225253442
Very high119151424303449
High, tall treeLow9514823203233
Average95141123233338
High95141122233139
Very high86131223223145
Table 10. Expected yield decline (YD) versus water saving (WS) for G3 (lemon and limes) under deficit-irrigation (management allowed deficit, MAD: 1.05, 1.10, 1.20, and 1.30) and climate (low, medium, high, and very high demand) scenarios.
Table 10. Expected yield decline (YD) versus water saving (WS) for G3 (lemon and limes) under deficit-irrigation (management allowed deficit, MAD: 1.05, 1.10, 1.20, and 1.30) and climate (low, medium, high, and very high demand) scenarios.
Plant DensityDemandDeficit-Irrigation Strategies
5%10%20%30%
WS
(%)
YD
(%)
WS
(%)
YD
(%)
WS
(%)
YD
(%)
WS
(%)
YD
(%)
LowLow146181128233739
Average135181027193730
High137171225243440
Very high116151223243238
Med, small treeLow136181128233738
Average14519929183828
High127171225243539
Very high116151124233337
Med, tall treeLow136181126243639
Average145191029193730
High127161226243440
Very high107151224243339
High, small treeLow126181127243640
Average146181028203732
High117161225243441
Very high118151224253340
High, tall treeLow9514923213236
Average95141023213435
High95131122223237
Very high86131124223037
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Darouich, H.; Karfoul, R.; Ramos, T.B.; Pereira, L.S. Setting Irrigation Thresholds for Building a Platform Aimed at the Improved Management of Citrus Orchards in Coastal Syria. Agronomy 2023, 13, 1794. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy13071794

AMA Style

Darouich H, Karfoul R, Ramos TB, Pereira LS. Setting Irrigation Thresholds for Building a Platform Aimed at the Improved Management of Citrus Orchards in Coastal Syria. Agronomy. 2023; 13(7):1794. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy13071794

Chicago/Turabian Style

Darouich, Hanaa, Razan Karfoul, Tiago B. Ramos, and Luís S. Pereira. 2023. "Setting Irrigation Thresholds for Building a Platform Aimed at the Improved Management of Citrus Orchards in Coastal Syria" Agronomy 13, no. 7: 1794. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy13071794

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop