1. Introduction
The Mediterranean region is often identified as a climate change hotspot [
1]. More intense warming and drying are projected in this area than in other parts of the world [
2]. Even considering adaptation efforts, new climatic conditions will impact the countries and the communities of the Mediterranean region in a variety of ways, with different magnitudes and consequences [
3]. It is widely accepted that agriculture will be one of the most affected sectors with projections of decreased yields, greater irrigation needs [
3,
4], and reduced water availability [
5,
6].
Barley is one of the most resilient crops to water and temperature stress [
7,
8,
9,
10]. Thanks to its deep rooting system, it can endure drought [
8,
10], and it is believed to adapt to stressful conditions better than other cereals [
7,
8]. Barley thrives under both cold and warm climates [
10], and thanks to its relatively short growing season, it can be grown in areas with limited rainfall [
9]. Thanks to its resilience, barley thrives in the Mediterranean basin, serving multiple purposes, including livestock feed, beer and other alcohol production, and human nutrition [
9,
10]. Barley can also be a crucial stage in crop rotation, with benefits related to pest and disease control [
10].
Notwithstanding its resilience and endurance, barley production is expected to be impacted by climate change [
11,
12]. The foreseen alteration in temperature and precipitation patterns that will affect the Mediterranean, along with the increased probability of extreme events, will reduce yields and increase the risk of crop failure [
3,
10]. An increase in dry conditions has already been reported by the IPCC for the whole Mediterranean region (defined by the cited source as the area ranging from northern Italy to the north of Morocco and Tunisia, and from Spain to Lebanon), and the trend is projected to continue in future climate change scenarios, posing inevitable threats to water availability [
3]. Under the projected drier and warmer environmental conditions, even this stress-resilient crop might face crop losses [
10].
Several studies have investigated and tried to quantify the effects of climate change on barley production in the Mediterranean region. Bento et al. [
12] found that under climate change, the production of barley and wheat will decrease in the southern region of the Iberian Peninsula, mainly due to an increase in spring maximum temperature. On the contrary, an increase in yields is projected in the northern parts of the Iberian Peninsula, with the main driver being early winter warming [
12]. Cammarano et al. [
11] focused on the whole Mediterranean basin and defined three scenarios of water availability: “dry”, “mid”, and “wet”. By mid-century and under Representative Concentration Pathway (RCP) 4.5, the study projected a decrease in yields of 27% under the “dry” scenario and increases of 4% and 8% under the “mid” and “wet” scenarios. The authors stressed the importance of soil water content at the beginning of the growing season as a critical factor for barley growth [
11]. This was supported by Al-Bakri et al. [
13] when focusing on barley growth in Jordan. The same authors analysed three case studies representing semi-arid Mediterranean environments in Jordan (BSk in Köppen–Geiger classification [
14]), considering RCP4.5 and RCP8.5 at both the mid-century (2030–2050) and end-century (2080–2100) time periods [
15]. The results of their study indicate a decrease in barley yields between 5% and 30%, with an average decrease of 12% in grain yields and of 16% in biological yields due to drier environmental conditions [
15].
Due to its suitability for water-stress conditions, the FAO AquaCrop model is often used to study barley growth; however, only a few cases exist focusing on the Mediterranean area. Dhoiub et al. [
16] carried out a multicriteria analysis to evaluate barley growth in hilly agrosystems in Tunisia. Alaya et al. [
17] adopted barley as one of the sample crops to evaluate the scalability of AquaCrop results through GIS. Both of these works utilize the same parametrization, which was developed for salinity stress and irrigated conditions by El Mokh et al. [
18]. López-Urrea et al. [
19] also provided a parametrization for barely in a Mediterranean setting, but it focuses on high-yielding barley.
When interested in the applications of AquaCrop for climate change impact assessment on barley, we have to move outside the Mediterranean. In a study on Northern Serbia, Daničić et al. [
20] found that while the flowering time and the overall phenology of the crop changed, yields were largely unimpacted in the mid- and end-century time periods. Dubey and Sharma [
21] focused on the basin of the Banas River in India and showed how in the mid-century (2021–2050) time period and under RCP4.5 yields were projected to increase. Yawson et al. [
22,
23,
24] studied climate change’s impact on barley production for malting beer and for food security in the UK. They suggest possible increases in yields if the effects of climate change alone are considered [
23]. Lastly, Arce-Romero et al. [
25] applied the model to two Mexican case studies, projecting yield decreases in the future, smoothed if mulching or a change to later sowing dates was adopted.
All the aforementioned studies include, to differing extents, a calibration of the crop’s parameters. Such calibration is generally conducted through field experiments [
19,
20,
21,
26], literature data [
24,
25], algorithms [
27], or satellite imagery [
28,
29]. For experimental calibration, the process requires multiple years of data and agronomic knowledge, making the process resource- and time-intensive. After calibration, validation is required. This is normally conducted by comparing the results obtained in the model using the calibrated parameters and experimental data and through a set of statistical indicators such as Root Mean Squared Error, Index of Agreement, Nash–Sutcliffe efficiency, and Goodness of Fit [
12,
19,
26,
30].
The process of calibration enhances the accuracy of the results; however, as highlighted by Coudron et al. [
31], these can still be affected by other uncertainties within the model, such as those linked to local climate variability. When considering the accuracy of the results, it is important to remember that AquaCrop was developed for field-scale modelling [
32]. Therefore, it requires high-quality inputs to model inter-annual changes in productivity. However, such a focus on the crop parameters becomes less important in larger-scale applications of the model, as in the work by de Roos et al. [
33].
Aim of the Study
Our study focuses on the Spanish province of Almería, an area of transition from a Mediterranean to a semi-arid climate. Almería is an interesting case study as it relies on the agricultural sector for its economy, particularly on greenhouse-based intensive horticulture which uses a large amount of water [
34,
35]. In Almería, barely is the main field crop, in most cases grown under rainfed conditions and mainly planted in northern areas of the province, in which greenhouse intensive agriculture is not present [
36]. Due to Almería’s situation in southern Spain, climate change threatens to make barley cropping no longer a source of livelihoods for local farmers [
12]. To improve the understanding of the impact of climate change on agriculture in a semi-arid area as Almería, and provide tools for developing adaptation strategies, the work is guided by two central research questions:
How and to what extent will barley production change in the mid- and end-century time periods under different climate change scenarios in terms of multi-year trends?
What is the efficacy of irrigation, mulching, and changing the sowing date as adaptation strategies to address the effects of climate change in the province?
Research question 1 will provide knowledge on future trends in barley production at the provincial scale, allowing for more readily available information for local stakeholders to drive climate action. Research question 2 will provide insights on the efficacy of climate change adaptation practices that have been suggested by practitioners in the region [
37], giving a more profound understanding of possibilities to adapt barley production to climate change. It should be noted that the impact of climate change on rainfed barley production in this study is intended solely in terms of overall yields, and the agronomic aspects explaining these are beyond the scope of the current paper.
We applied the FAO AquaCrop model [
32] in its Python implementation—AquaCrop-OSPy [
38]—at the scale of Almería province and for two 30-year time horizons, mid-century (2041–2070) and end-century (2071–2100), comparing them with a baseline period (1985–2014). For each time period, three Coupled Model Intercomparison Project phase 6 (CMIP6) Tier 1 Shared Socioeconomic Pathway (SSP)-based scenarios [
39] were analysed: SSP1-2.6, SSP2-4.5, and SSP5-8.5. To account for possible losses in soil water availability at the beginning of the growing season, three levels of Initial Soil Water Content (ISWC) were considered: 30% of the Total Available Water (TAW), 20% TAW, and 10% TAW.
The focus on multi-year trends, the limits of Mediterranean parametrizations for barley [
18,
19], and the time and resource requirements for crop calibration led to the development of a framework where AquaCrop was run with the standard barley crop parameters available in Aqua-Crop-OSPy (based on an Ethiopian case study), without local recalibration. To evaluate the applicability of the proposed approach, the model’s results from the baseline period were compared through a statistical analysis with a historically observed dataset of barley yields for the target province from the Spanish Ministry of Agriculture. This allowed us to retrieve a reference ISWC and the minimum number of years needed to be studied together to retrieve meaningful information on barley yield trends.
The results obtained in this work can be used by local stakeholders and policymakers to develop long-term strategies for climate change adaptation. The overall simulation framework can constitute an example to be followed in other regions, to obtain preliminary indications of crop yield changes. Having acknowledged the lack of studies applying AquaCrop to study climate change impact on barley in the Mediterranean region, this study also suggests the possibility of extending the use of this model with such purpose in this area.
3. Results
3.1. Climate Change Projections for Almería Province
Climate change projections for the province of Almería indicate a warmer and drier future. Minimum and maximum temperatures, along with potential evapotranspiration, are projected to increase while precipitation gradually decreases over time (
Figure 3). The SSPs exhibit different trends, with SSP5-8.5 projecting the largest rate of change concerning the baseline period (1985–2014), for all four variables. Under this scenario, the increase in yearly minimum temperature can reach over 60%, with yearly maximum temperature and potential evapotranspiration increasing by 40% and 25%, respectively, by the end of the century. Under SSP5-8.5, precipitation is projected to decrease by 20% (mid-century) and 40% (end-century) when compared to the baseline period, while SSPs 1-2.6 and 2-4.5 show more moderate decreases, reaching 10% and 20%, respectively, by the end of the century.
Water Availability
Water availability trends in Almería province were quantified by calculating the yearly difference between precipitation and potential evapotranspiration, expressed as the percentage change from the mean of the baseline period (1985–2014) (
Figure 4).
Since precipitation will not mirror the increase in potential evapotranspiration (
Figure 3), water availability is expected to continue to decrease in the region until the end of the century in all three analysed scenarios. However, there are differences in the changing trend between scenarios, in line with their different expected climate forcings.
Under SSP5-8.5, there is the greatest decrease in water availability. The difference between precipitation and potential evapotranspiration reaches an average of −60% of the baseline period by the end of the century, with a more pronounced decrease in precipitation and increase in potential evapotranspiration after the 2060s.
Water availability decreases at a lower rate in SSP1-2.6 and SSP2-4.5, with the results showing a constant trend with a decrease between −10% and −30% of the baseline period average.
The uncertainty range in the model’s ensemble for scenarios SSP1-2.6 and SSP2-4.5 shows the possibility that water availability in Almería may remain equal, or even increase relative to the baseline period. However, and with the notable exception of the period up to 2060, this is not the case with SSP5-8.5, which shows a clear drying trend over the entire period.
3.2. Selection of the Reference Initial Soil Water Content (ISWC)
The modelling of barley production in AquaCrop requires inputting a reference value of Initial Soil Water Content (ISWC) in the form of the percentage of Total Available Water (TAW) in the soil. A statistical-based approach was used to better validate and match data with the real situation in the field without having to undergo field data campaigns. AquaCrop was run for the baseline period (1985–2014) with different values of ISWC, which were then compared to observed historical data of rainfed barley yields using a set of statistical analyses (paired t-test, Spearman and Pearson correlation coefficients, and error analysis). The results of this analysis are described in the following sections.
3.2.1. Paired T-Test
The AquaCrop-OSPy modelled yield distributions with the least significant mean difference to the historical observed dataset were obtained with 20% and 30% TAW (
Table 3 and
Figure S2 display instead the yield distributions obtained with the different ISWC values). Although modelled results using a TAW of 40% show a
p-value of less than 0.05 (and thus the mean difference between modelled and observed paired samples may be significantly different from 0), the results were notably higher than for all other remaining TAW values (
Table 3). For this reason, a TAW of 40% was also included as a model parameterization in our analysis.
3.2.2. Pearson and Spearman Correlation Coefficients
The relationships between the historical observed yield dataset and the modelled yield datasets for the baseline time period under all selected TAW levels (20%, 30%, and 40%) were analysed using Pearson and Spearman correlation coefficients. Statistical correlations were computed for multiple rolling average windows of yield datasets and the results are depicted in
Figure 5.
Results show that, under all TAW levels, both correlation coefficients increase with an increasing rolling average window, reaching a maximum with a 13-year rolling average. Additionally, there are no significant differences between TAW levels in either coefficient, with correlation coefficients exceeding 0.6 observed above the 10-year rolling average windows.
The statistical significance of results for each rolling average window and TAW level (20%, 30%, and 40%) were evaluated for both correlation coefficients, using a
t-test with a threshold of 0.05. In both correlation coefficients, statistical significance was verified for rolling average windows above 9 years, independently of the TAW level. Detailed results of this analysis are presented in
Table S1.
3.2.3. Error Analysis
We used Root Mean Squared Error (RMSE) and relative Root Mean Squared Error (rRMSE) to analyse predicted error levels in modelled yield datasets, accounting for the three TAW levels studied. The results show that both RMSE and rRMSE become lower with an increasing rolling average window (
Figure 6). The results include a steep decrease in both RMSE and rRMSE between 0 and 5 years of the rolling average window. Values reach low levels after that, namely below 0.4 of RMSE and 30 of rRMSE. Different error behaviours can be described depending on TAW levels, with a 30% TAW distribution reducing the error at a greater pace with the average window (going below 0.2 and 10 after the 10-year rolling average), and the curves of 20% TAW and 40% TAW plateauing at around 0.3 and 20 for each coefficient.
Based on these results, a value of 30% TAW was selected as the reference ISWC for the model parameterization of the baseline period. This value leads to a reduction in errors while maintaining a good correlation and a good similarity with the historical observed dataset. This finding means that the analysis of future scenarios of decreases in water availability and adaptation options will be based on a comparison with the 30% TAW threshold, considered as the current condition.
3.3. Climate Change Impact on Rainfed Barley Yields
Climate change impacts on rainfed barley yields as percentage changes in yield are presented in
Figure 7 for three climate scenarios (SSPs 1-2.6, 2-4.5, and 5-8.5) and three ISWC parameterizations (TAW of 30%, 20%, and 10%), by the mid- and end-century time periods.
The results indicate a decrease in yields for all the analysed cases, except for scenarios where 30% TAW is sustained, or, in other words, where the current ISWC at the sowing date is maintained. Even in this case, model uncertainty points toward the possibility of a potential loss of yield, with model results spreading to negative values in all three climate scenarios, and for both time horizons.
The largest changes in average yield reach −55.1% for a 10% TAW level at the end of the century in SSP5-8.5. For this TAW level, the minimum change observed is −37.6% at the mid-century period for SSP1-2.6.
In scenarios where the ISWC is maintained at the current level of 30% TAW, results range from an average yield increase of 14% under SSP2-4.5 at the mid-century period to a 4% increase under SSP5-8.5, for the end-century period.
Table S2 provides in-depth details on the average yield changes for each timeline and climate change scenario analysed.
In general, the ISWC parameterization affects barley yield more than climate change scenarios over both time horizons. For example, in the mid-century period, the results show minimal changes in yield for the same ISWC across climate scenarios. Nonetheless, there were significant changes for different ISWCs within the same climate scenario. This tendency remains at the end of the century, although the magnitude of change between ISWC parametrizations increases substantially compared to the previous period.
3.4. Climate Adaptation Pathways Analysis
3.4.1. Irrigation
To evaluate the climate adaptation potential of barley irrigation in Almería, each combination of climate scenarios and ISWC parameters was modelled in AquaCrop-OSPy using three different irrigation parameterizations: (a) absence of irrigation (rainfed); (b) irrigation triggered when soil water content drops below 20%; and (c) irrigation triggered when soil water content reaches 0% TAW (
Table 4).
As for barley productivity impacts, ISWC parameterization influences irrigation needs more than climate scenarios. More water available at the beginning of the growing season means lower irrigation needs regardless of the climate scenario and time horizon.
Additionally, the choice of the irrigation threshold has different impacts according to the ISWC parameterization. For example, at the 10% and 20% ISWC levels, a TAW irrigation threshold set at 0% is not useful to avoid productivity losses. However, for the same ISWC levels, a 20% TAW irrigation threshold always increases productivity, across all scenarios and time horizons.
Another example of the importance of ISWC is the use of a 20% TAW irrigation threshold across scenarios. Somewhat contra-intuitively, this threshold yields better productivity results in the 10% ISWC than in the 20% and 30% ISWC scenarios, for all climate scenarios and time horizons. The reason for such an effect is that a lower soil water content at the beginning of the growing season causes the soil TAW to drop below the defined threshold more often, thus triggering more (and earlier) irrigation events in the model, leading to an increase in productivity.
The additional water requirements for barley irrigation in Almería linked to the results of
Table 4 are presented in
Table 5.
Using the example above of a 10% ISWC level scenario, and with the aim of ensuring that irrigation is triggered when TAW reaches 20% (i.e., when enough rain has fallen to go from 10% TAW to more than 20% TAW but then the soil dries), the additional water requirement varies between 355.9 m3/ha (364.5 m3/ha) and 386.7 m3/ha (396.6 m3/ha), by the mid-century (end-century) period.
Similarly, with a 20% TAW irrigation threshold and for a 20% TAW ISWC level, water requirements are still substantial. In comparison, for the 30% ISWC level, irrigation needs drop to residual levels since irrigation is seldom triggered, on account of soil moisture being kept at today’s levels, even under warming scenarios.
An extreme example of the importance of soil water content is the already mentioned scenario of irrigating when ISWC is under 10% TAW levels. In this case, setting the irrigation threshold to 0% TAW (thus using less water) would still require between 55.4 m3/ha and 109.6 m3/ha by the end of the century to reduce losses in productivity.
As expected, irrigation needs are higher with a higher threshold and lower ISWC. The scenario of a 10% ISWC level with a 20% TAW irrigation threshold is the one requiring more additional water during the growth phase. However, this is also the scenario leading to higher positive yield changes (
Table 4). These can be explained by the triggering of (more) irrigation events at earlier stages of the plant growth phase when the above-mentioned couple of the ISWC and irrigation threshold is verified.
Conversely, the scenario where ISWC is kept at the current 30% TAW at the sowing date highlights the importance of soil water content. This scenario requires almost no additional water, either because irrigation is triggered rarely (at 20% TWA threshold) or not at all (at 0% TWA threshold) but still presents positive changes in barley productivity, even in a no-irrigation situation.
3.4.2. Mulching
The second adaptation strategy modelled was the application of mulches coupled with irrigation. The results of this strategy are reported in
Table 6. Yield change and seasonal irrigation requirements are compared against values obtained for the same irrigation conditions but without the application of mulches (as per
Table 4 and
Table 5 above).
The findings indicate that mulching is partially efficient in reducing irrigation needs while improving yields. The effect of mulches appears uniform throughout the time horizons, climate scenarios, and model parameterization. Mulching always lowers water needs while improving yields regardless of the ISWC level and irrigation threshold.
However, the effect of reducing water needs is more heterogeneous and affected by a much larger variance than increasing yields, with marked differences between climate scenarios, ISWC, and TAW irrigation thresholds. For example, the increase in yields promoted by the use of mulches never exceeds 10% (when compared to the irrigation-only strategy). Despite this, the reduction in irrigation water needs can reach up to 40% by the mid-century period under SSP5-8.5.
Mulches are more effective in reducing water needs for a 0% TAW irrigation threshold in the 10% TAW and 20% TAW ISWC parameterization and the 20% TAW irrigation threshold with a 30% TAW ISWC.
If annual inter-variability is factored in, these results open up the possibility that irrigation coupled with mulches could halt productivity losses in wetter years, for example, if the ISWC level is kept at least at 20% TAW by the sowing date. Mulching efficiency in absolute values can be found in
Table S3.
3.4.3. Changing the Sowing Date
The last adaptation option tested was the change in sowing date (anticipation or delay) from the reference date of the 10th of November. Percentage changes in rainfed barley yield obtained for different modelled sowing dates are presented in
Table 7 for all climate change scenarios, time horizons, and ISWC parameterizations.
Results do not show a clear trend of improving barley yields with earlier or later sowing dates. Changes in productivity are within the −17.8% to +15.5% range, with only the 4th of November (one-week anticipation of sowing) showing positive values for all scenarios, time horizons, and ISWC parameterization.
By the mid-century period, the results show productivity improvements with up to 2 weeks earlier sowing dates but only for higher warming levels (SSP2-4.5 and 5-9.5) and lower ISWC (10% to 20%). However, this effect is not carried through to the end of the century, where productivity gains appear connected to a later sowing date of about 1 week.
These results may point to some positive gains of an earlier sowing date (in line with a warmer climate but not necessarily with a drier one). There is no clear evidence that an earlier or later sowing date would significantly improve barley yield in Almería.
5. Limitations and Future Research
Our study has various limitations that have been made explicit in the paper.
The most significant limitation is the use of uncalibrated barley parameters in AquaCrop using the standard parameters available for a short-cycle barley variety grown in Ethiopia. This introduces uncertainties linked to the crop’s responses to climate occurring outside the modelled growing season (i.e., in late spring). Future work in this area should carry out a preliminary calibration to improve the accuracy of the results. Such results could be compared to the ones illustrated in this paper to evaluate the benefits of crop calibration. The modelling approach adopted does not capture the inter-annual variability of crop growth. While this approximation is suitable for the paper’s focus on multi-year trends, it does not allow for the analysis of finer time scales. This limits the applicability of this approach to long-term analyses. The results are therefore intended as a quantification of a range of probable future trends in rainfed barley growth in Almería. Future research can improve the modelling framework and investigate inter-annual variability in rainfed barley yields.
Secondly, the accuracy of modelling cropland management practices associated with climate adaptation in Almería could benefit from the use of field data for validation, e.g., calendar types, sowing dates, crop characteristics, or in its absence, from the interaction with local stakeholders.
Furthermore, the modelling used average values for the whole of the Almería province, thus not allowing for the possibility of studying spatial patterns in barley productivity change. Future studies could implement a modelling strategy based on pixel data, such as the one proposed by de Roos et al. [
33].
Then, in our study, adaptation strategies were kept identical across climate change scenarios. However, different scenarios also include different socio-economic, technological, and behavioural developments. Future research would benefit from modelling adaptation strategies that reflect these different narratives. This should come in tandem with an overall improvement in the accuracy of the modelling framework.
Lastly, this paper focuses on the effect of climate change on overall rainfed barely yields in Almería province, as climate was identified as the key variable of interest for local stakeholders. A deep investigation of the physiological reasons that cause such changes in yields is beyond the scope of the present study. Future studies might be interested in using the results of this study as a basis for the investigation of the agronomic aspects of such yield change.
We plan to address these limitations in future works, using the results and methods of this paper as a basis.
6. Conclusions
This paper applied the AquaCrop-OSPy model to evaluate the impact of climate change on rainfed barley cropping in the previously unstudied province of Almería, Spain. In this study, irrigation, mulching, and changing the sowing date were included in the modelling to evaluate their efficacy as climate change adaptation solutions.
Our work suggests a modelling framework for AquaCrop-OSPy without using field data or locally recalibrating the barley parameters. Such an approach was evaluated to be effective in retrieving meaningful information if at least 9 years of data were analysed. A reference Initial Soil Water Content, in this case identified as 30% TAW, allowed for finding a distribution of modelled yield that better matched the historical observed dataset.
Subsequently, AquaCrop-OSPy was run for two future time periods (2041–2070 and 2071–2100) under three climate change scenarios (SSP1-2.5, SSP2-4.5, and SSP5-8.5). To account for possible decreases in soil water availability, varying Initial Soil Water Content levels (10%, 20%, and 30% TAW) were considered for each climate change scenario and time horizon.
Results indicate a possible reduction of −55.1% at the end of the century, for a 10% TAW soil water content at sowing. In the case of an Initial Soil Water Content of 30% TAW, the results indicate yield increases of up to 14% and no decreases in average yields, for the same time horizon. Therefore, it can be concluded that the most important parameter determining future rainfed barley productivity changes is the Initial Soil Water Content (ISWC).
Irrigation emerges as an effective adaptation option in cases of reduced Initial Soil Water Content but carries the burden of exploiting already limited water resources. However, such water requirements were quantified to be slightly less than 400 m3/ha during the crop growing season, which would account for a minimal percentage of the water resources of Almería province, even considering climate change.
Two additional adaptation strategies were modelled. Mulching was effective in limiting irrigation water needs and partially improving yields. Anticipating or delaying the sowing dates of rainfed barley, on the other hand, was not shown to improve rainfed barley production. The findings, however, support the conclusion that the ideal sowing window might be reduced in the future.