Drought prediction is crucial, especially where the rainfall regime is irregular, such as in Mediterranean countries. A new combined drought indicator (CDI) integrating rainfall, soil moisture and vegetation dynamics is proposed. Standardized precipitation index (SPI) is used for evaluating rainfall trends. A bucket-type soil moisture model is employed for keeping track of soil moisture and calculating anomalies, and, finally, satellite-based normalized difference vegetation index (NDVI) data are used for monitoring vegetation response. The proposed CDI has four levels, at increasing degrees of severity: watch, warning, alert type I and alert type II.
This CDI was thus applied over the period 2003–2013 to five study sites, representative of the main grain-growing areas of SW Spain. The performance of the CDI levels was assessed by comparison with observed crop damage data.
Observations show a good match between crop damage and the CDI. Important crop drought events in 2004–2005 and 2011–2012, distinguished by crop damage in between 70 % and 95 % of the total insured area, were correctly predicted by the proposed CDI in all five areas.
Drought is a recurrent phenomenon on the Earth's surface. It is triggered by lack of water, or “an extended imbalance between supply and demand” in the precise expression of Hobbins et al. (2016), and may have economic, social and environmental impacts (Wilhite, 2000). Drought is one of the most important natural disasters threatening our society. In spite of its relevance, there is no proper definition of drought. Tannehill (1947) called drought “the creeping phenomenon”, given the complexity of accurately delimiting its start time and end time and of adequately demarcating the spatial extent of its effects.
Wilhite and Glantz (1985) distinguished four main types of droughts according to how the effects were noticed: (i) meteorological, due to the scarcity of rainfall; (ii) hydrological, detected by low streamflow; (iii) agricultural, when soil water is not sufficient to maintain a crop; and (iv) socioeconomic, when it affects the normal functioning of society.
Drought occurs worldwide but it is especially frequent in the Mediterranean
region. In a recent analysis of a tree-ring-based reconstruction of the
summer season, the Palmer drought severity index (PDSI) (Keyantash and
Dracup, 2002) for the period from 1100 to 2012, Cook et al. (2015, 2016)
detected the gravity of recent events in the area, which were apparently induced by
anthropogenic activity. Combining two drought indices, one meteorological,
the Standardized Precipitation Index (SPI), for water supply, and the other
hydrological, the standardized precipitation–evapotranspiration
index (SPEI), for water loss tendency, Stagge et al. (2017) observed, for
the European continent in the period 1958–2014, that droughts were mainly
driven by a temperature rise with the inherent increase in the
evapotranspiration rate, whereas rainfall did not change appreciably. In
the southwestern United States, Ting et al. (2018) found that, under a
One additional problem of drought is that it can spread towards other regions, as Herrera-Estrada et al. (2017) discovered in their Lagrangian analysis in several Earth regions. Andreadis et al. (2005) have elaborated on severity-area-duration maps, modifying an earlier proposal of Dalezios et al. (2000) for severity-duration-frequency maps. Therefore, drought is a present-day risk at least for a part of our society.
Drought characterization depends on the perspective of the user. The meteorological drought is possibly the simplest type to evaluate since it is reduced to a mere consideration of the rainfall. The two main meteorological drought indices are those mentioned above, the PDSI and SPI. Hydrological drought requires the conversion of rainfall into runoff, which can be done with the help of a hydrological model; the SPEI, for instance, is a widely used hydrological drought index. Nevertheless, Van Loon and Van Lanen (2012) have explored in depth the definition of hydrological drought, starting from the time perspective of the phenomenon, and distinguishing several types in terms of the sequences rain to snow, wet to dry, cold snow and warm snow seasons and what they denominated as classical rain deficit. The use of a simple hydrological model and the establishment of some threshold values allow Van Loon and Van Lanen (2012) to determine the drought occurrence in several regions with distinct climate types. Drought severity is a function of the available water storage units, as Van Loon and Laaha (2015) explained in the review of an Australian dataset. Hobbins et al. (2016) have modified the SPEI index by representing the potential evapotranspiration and the atmospheric evaporative demand on a proper physical basis, rather than on the air temperature as a proxy of it. Their evaporative demand drought index (EDDI) is a useful indicator of drought extent, as was shown by McEvoy et al. (2016) in the conterminous US. The estimation of the agricultural drought index is somewhat similar to that of the hydrological drought one, with the additional complexity of crop behavior. Several models have been proposed for the agricultural drought index estimation. As Perrin et al. (2001) warned, and Orth et al. (2015) later confirmed, the models set up to describe soil water evolution for this purpose must be very simple and limited to soil water balance. Hunt et al. (2009), Khare et al. (2013) and Sohrabi et al. (2015) proposed reasonable soil water balance models, differing only in their characterizations of rainfall infiltration, in order to prevent the generation of excess rain, deep percolation and actual evapotranspiration rate.
The different drought indices represent distinct aspects of drought. Therefore, to gain a wider perspective, Kao and Govindaraju (2010) introduced the use of copulas in a new drought indicator denominated the joint deficit index (JDI), based on the SPI for both precipitation and streamflow. Hao and AghaKouchak (2013) formulated another copula, the multivariate standardized drought index (MSDI), consisting of the SPI and a standardized soil moisture index (SSI). This index was very useful for detecting the drought onset and duration. Alternatively, Zarch et al. (2015) used two separate indices to assess droughts, the SPI and the reconnaissance drought index (RDI). A different approach was suggested by Hao et al. (2017) with a categorical drought prediction model, the U.S. Drought Monitor (USDM), which proved to be highly adequate for early warning. Azmi et al. (2016) developed a data fusion-based drought index, grouping different indices with a clustering method.
The impact of drought on vegetation can be by means of several indices. Kogan (1995) proposed a vegetation condition index (VCI) based on the normalized difference vegetation index (NDVI), which is a good indicator of vegetation status, by combining the radiance of the visible and infrared wavelengths to assess the drought effects. Some other indices have been suggested, since NDVI is sometimes influenced by other environmental factors (Quiring and Ganesh, 2010). The normalized difference water index (NDWI) was introduced by Gao (1996), and, using radiances in a higher wavelength range than that of NDVI, it is less affected than the latter by atmospheric conditions; it is also more sensitive to drought than other indices (Gulágsi and Kovács, 2015). The Joint Research Centre of the European Comission uses the fraction of absorbed photosynthetically active radiation (fAPAR) generated from the signals acquired by the Project for On-Board Autonomy - Vegetation (PROBA-V) sensor.
The abovementioned methods can be used to evaluate the impact of drought on agricultural productivity in regions worldwide, as Sepulcre-Cantó et al. (2012) have shown for Europe. These authors proposed a combined drought indicator using SPI, fAPAR and soil moisture calculated from a regional hydrological model. For the management of local policy and mitigation actions, such as farm-scale insurance schemes, smaller spatial scales than those used by Sepulcre-Cantó et al. (2012) are required.
The main objective of this work is to assess agricultural drought by means of a combined drought indicator (CDI), based on SPI and anomalies in soil moisture and NDVI. This new CDI is thus related to crop damage data in rainfed wheat-producing regions in southern Spain at the agricultural province level, which corresponds to the most important item of available yield data. It is expected that this new CDI will be useful at the local policy level and for planning farm-scale insurance schemes.
This study was made in Andalusia, southern Spain, during the 10-year period between 2003 and 2013.
Andalusia has a Mediterranean climate with dry, hot summers (Köppen-Geiger climate Csa, Peel et al., 2007). Since the main source of water is rain caused by the western and southwestern winds carrying moist air from the Atlantic Ocean, the distribution of precipitation is conditioned by the orography of the region, with a main decreasing gradient from west to east.
Location of the study area (grey) and selected representative points (blue dots) within the areas cultivated with cereal (white).
The effect of drought on agricultural production was evaluated in five
representative areas, in each of which, four representative locations were
selected in a two-step procedure. First, the distribution of the land use
class “non irrigated arable land” within the study area was analyzed, as
shown in Fig. 1. This land use distribution is derived from the regional
land use map (SIOSE: Soil Occupation Information System of Spain) applied
to Andalusia, 2005, equivalent to the European CORINE database, on a scale
of It contains homogeneous land use of rainfed wheat within each
pixel (with no other land uses present in it); It lacks external landscape elements, such as ponds, roads,
canals, houses or natural vegetation patches that could distort the NDVI
signal; It has continuous wheat cultivation during the study period (no fallow
period).
The SPI expresses the deviation of rainfall from its long-term mean. SPI is calculated by fitting the precipitation data to a gamma distribution, after which it is transformed into a normal distribution. The SPI values can then be interpreted as representing the number of standard deviations by which the observed anomaly deviates from the long-term mean. SPI was calculated over 1-, 3- and 6-month periods, using precipitation series of between 42 and 69 years, namely SPI-1, SPI-3 and SPI-6.
SPI-1 is theoretically best related to meteorological drought, together with short-term soil moisture stress, especially in periods when crop growth is sensitive to them (Guttman, 1999). SPI-3 has been shown to reflect short to medium seasonal precipitation trends (Guttman, 1999). Bussay et al. (1999) and Szalai and Szinell (2000) evaluated the relationship between SPI and agricultural drought through soil moisture and found that SPI-2 and SPI-3 yielded the best results. Other authors (Ji and Peters, 2003; Rossi and Niemeyer, 2012) have reported a high correlation between SPI-3 and vegetation response and, therefore, deemed this index to be the best suited for evaluating agricultural drought. They deemed SPI-6 to be the best one for identifying longer-term or seasonal drought trends.
The program “SPI_SL_6.EXE”, developed by the
National Drought Mitigation Center, University of Nebraska-Lincoln, was used
to calculate SPI. Details of this method can be found in McKee et al. (1993)
and Lloyd-Hughes and Saunders (2002). The same classification used by McKee
et al. (1993) was used (Table 1), and a threshold value for defining a
drought of SPI >
Classification of droughts according to SPI and their probability of occurrence following McKee et al. (1993)
SPI values were calculated for each of the five agricultural regions selected: Campiña de Cádiz (Cádiz), Campiña Baja (Córdoba), Pedroches (Córdoba), Norte/Antequera (Málaga) and La Campiña (Sevilla). The climate series selected in each region was the one at their particular weather station that had the longest available series.
The deviation of the soil moisture from its long-term mean was expressed as
a soil moisture anomaly index (SMAI). SMAI values were calculated for each
of the five selected agricultural regions, similar to those of the SPI. To
obtain this index, we first calculated soil moisture dynamics by means of
the simple water balance model of Brocca et al. (2008). The long-term mean
soil moisture was taken as the 10-year mean in the study period (2003–2013).
In this water balance model, the water depth in the soil profile,
Parameters for the water balance model used in this study.
The soil moisture anomaly index (SMAI) is then given by Eq. (5):
Different agricultural drought studies have used satellite-based vegetation indices as their main advantage is their spatial and temporal resolution. NDVI values represent the plant chlorophyll content, which is why they are highly suitable for identification of agricultural drought. Limitations in its use are related to the fact that NDVI may reflect non-drought- related stress conditions, such as plant disease, and that soil properties can induce a bias in its response. Therefore, it is important to use NDVI-based drought evaluation in combination with other indices based on precipitation or soil water, as is the case here. NDVI anomalies express deviations in NDVI from its long-term mean, and these were evaluated on a monthly basis but only taken into account from November to April, which is the normal growing season for rainfed winter cereal in Andalusia. Only during this period can NDVI and its anomalies be expected to transmit information on rain-fed cereal growth. The long-term mean NDVI was taken as the 10-year mean in the study period (2003–2013).
Thanks to its spatial continuity, NDVI trends could be analyzed for 20
different points; i.e., four points or pixels were analyzed in each of the
five agricultural regions selected. This analysis yielded a total of 20
spatially different NDVI anomaly indices. The NDVI anomaly index was
calculated using Eq. (6):
The main idea behind the combined drought indicator (CDI) for identifying agricultural drought is an idealized cause–effect relation between water deficit and yield. There are different phases in this relationship: a precipitation deficit (phase 1) leads initially to soil water deficit (phase 2), which, if prolonged over time, will result in crop water stress and be reflected in the NDVI observed (phase 3), which finally generates a reduction in cereal yields (phase 4).
In its simplest form, this CDI would allow us to identify which cause–effect relationship phase the agricultural system has reached in the event of a drought. This indicator would then allow the establishment of a series of drought warnings, depending on that phase. The CDI should be seen as a first step towards designing that warning system.
This study proposes a CDI that combines three indices:
SPI-3 to identify the first level of precipitation deficit (phase 1) SMAI to identify anomalies in the soil moisture (phase 2) NDVI anomalies to characterize the subsequent effect of soil water stress on
crops (phase 3).
The warning levels suggested for the CDI proposed are given in Table 3. They aim and are expected to help policy makers to prepare and take action
in the case of droughts.
Classification of the combined drought indicator (CDI).
The CDI uses three different levels; the first two, watch and warning, indicate that a drought could be imminent. The highest level of the CDI is “alert”. The two types of alert include those cases in which a meteorological drought results in a rapid yield decrease. The type I alert can occur even without a previous anomaly in soil moisture values, which could be related to intense droughts occurring during sensitive phenological phases of the crop. Therefore, a type I alert depends on only two indicators, SPI-3 and NDVI. The type II alert is based on all three indicators composing the CDI (SPI-3, SMAI and NDVI) so that these give firmer evidence for the existence of an agricultural drought.
The insurance area data and those of areas affected by drought per agricultural season for rainfed cereal were given by Agroseguro, the Spanish agricultural insurance provider. These data were disaggregated for each area of the five under study and each agricultural season, from 2002–2003 to 2011–2012. Note that data for the last year of the study, 2012–2013, were not provided. Crop intensity damage is expressed as the percentage of surface area that was filed for damage with respect to the total insured area and is available on an agricultural region scale. Crop damage of close to 100 % indicates important losses during that year.
The SPI values calculated over a 3-month period (SPI-3) reflected short–medium term moisture conditions and provided an estimate of the seasonal precipitation that was useful for agricultural purposes. In our study area, defined in Sect. 2.1, SPI-3 values at the end of April revealed the precipitation trends during the plant reproduction stage and the grain development. SPI-3 at the end of December showed moisture conditions at the start of the growing season.
Figure 2 gives the trends in SPI-3 for all five selected agricultural
regions. The trends are similar in all regions, with SPI-3 values moving
periodically around the long-term mean or 0 value. In the driest years, one
can observe the highest negative peaks. For example, during the agricultural
year 2004–2005, which was very dry, negative values of up to
Variation of the standardized precipitation index over 3 months
(SPI-3) during the period studied (2003–2013) in the five selected
agricultural regions:
Figure 3 shows the variation in the SMAI over the period studied and for
each of the five agricultural regions. The two main dry periods of 2004–2005
and 2011–2012 are not consistently apparent. Generally, only two regions at
that time dipped below the
Variation of the soil moisture anomaly index (SMAI) during the
period studied (2003–2013) in the five selected agricultural regions:
Figure 4 shows a map indicating the spatial and temporal variability in NDVI values over Andalusia for the year 2004. Figure 4a indicates NDVI in April, right in the growing season, while Fig. 4b shows the same area after the cereal has been harvested. The color red indicates low values of NDVI, while green represents maxima of between 0.96 in April and 0.92 in June. When comparing the distribution of the main cereal-growing regions in the area in Fig. 1, these areas present the most important variation between the two images, with high values in April and low red ones in June.
NDVI values all over Andalusia in
Figure 5 shows the monthly variation in the NDVI anomaly for the four
selected pixels within the Campiña agricultural region. The pixels in
the other four agricultural regions are not shown, but their trend is
similar. There is, of course, an important spatial variability within the
area, such that some differences appear between the four study locations. This
can be attributed to different planting dates, crop varieties or soil
properties between the locations. Over the study period, however, the same
general temporal trends appear. Important negative deviations from the mean
indicate periods of high plant stress. Values of NDVI anomaly below
Variation of the monthly NDVI anomaly for the four selected
locations within the region “La Campiña” over the study period. Red lines
indicate values below the defined threshold of
Figure 6 shows the monthly evolution of CDI between 2003 and 2013 and compares its levels against crop damage data derived from agricultural insurance information. This occurs twice during the studied period on a regionalized scale, indicating the effects of a drought. The first time is during the agricultural year 2004–2005, with losses of between 73 % and 99 % in the five agricultural regions studied. Also, for the years 2011–2012, there was considerable crop damage of between 71 % and 92 %. A third season, 2009–2010, had medium to high losses, of between 44 % and 89 %. However, crop damage during this period is, rather, due to the effects of excessive precipitation, leading to water stagnation and erosion damage. This can be seen when comparing the annual precipitation values. For example, in the Córdoba agricultural region, with a mean long-term precipitation of 600 mm, the values for 2004–2005, 2009–2010 and 2011–2012 are, respectively, 423, 1179 and 433 mm.
Evolution of the combined drought indicator (CDI) from 2003–2013
and comparison with agricultural crop damage intensity (blue lines) for the
five agricultural regions studied:
The CDI accurately captured these two important drought periods. For the first area, Campiña de Cádiz (Fig. 6a), a series of drought warning levels were issued early in the agricultural year 2004–2005, followed by a type I alert in January. There was another type I and II alert in May–June. In other words, since the seeding and during the first months of crop growth, there was a continued series of drought warnings or alerts. In that particular year, 90 % of the insured area was reported as being damaged. In 2005–2006, the CDI registered another warning indication, but it did not lead to any damage to the crop. In September 2005 there was a type II alert, but that month is outside the cereal growth period and when the crop was seeded two months later, the situation had gone back to normal. In May 2006 another warning was issued due to a precipitation and a soil moisture deficit. However, the crop was already at the moment in its cycle when it was close to harvesting and it was therefore not affected so much. In 2009–2010, characterized by considerable crop damage, 89 % of the total insured area, there was only one alert, in November. As mentioned before, crop damage during that season was probably due to precipitation excess rather than drought. For the dry period of 2011–2012, the CDI accurately indicated that critical situation with a warning, followed by type I and II alerts in the period of February–April.
For the Campiña Baja region (Fig. 6b) the dry period of 2004–2005 was characterized by a continuous series of type I and II alerts from January to June, with two more alerts during the summer, outside the cereal growing period. In this region, the insured area damaged that year was also very extensive (95 %). In 2008–2009 a warning was issued that did not cause any yield losses, as only 15 % of the insured area was damaged. This can be explained by the fact that this situation did not occur at a time when the crop was sensitive. In another dry year, 2011–2012, a series of warnings were issued, from January to March, followed by, respectively, type II and I alerts in April and May. These all occurred at times when the crop was highly sensitive, so that it was seriously damaged in 90 % of the area.
In the Pedroches region (Fig. 6c), the two main dry periods were well predicted. The year 2004–2005 was distinguished by a series of type II alerts in January, February, March and May and a type I alert in June. This sequence of critical CDI levels was reflected in an insured crop area with 73 % of damage. In 2005–2006, although there were two types of stress situations, warnings and type II alerts from November to February, the damage rate was not a high one, only 15 % of the insured area. It is difficult to understand the underlying reasons for the good performance of the crop that year. For example, during the years 2008–2009, the incidents were clearly late in the year (May to July), a period when grain growth is not sensitive. The second dry period of 2011–2012 is marked by a number of type II alerts issued from February to April, at a time when the cereal is highly vulnerable. This is reflected in a 71 % damaged insured area.
In the Comarca Norte/Antequera region (Fig. 6d), the dry period of 2004–2005 was determined by several incidents early on, with a watch issued in November and a type II alert in January, the latter being the period of cereal nascence and other sensitive growth stages. That year, the damaged insured area was 88 %. In 2007–2008 there were two warnings and a type II alert, from December to February, but these did not lead to crop damage, as the damaged insured area was only 11 %. Again, the reason could be found in those droughts occurring during a period when the cereal was not too sensitive. During the second main dry period of 2011–2012, a number of type I and II alerts were issued between February and April. These corresponded to highly sensitive moments of the crop cycle, and damaged insured areas were consequently high that year, amounting to 83 %.
The last region, La Campiña (Fig. 6e), showed a similar trend, with 2004–2005 being identified as having an extremely high damaged insured area of 99 %. The CDI worked well in predicting this, as there were multiple and continued alerts; i.e., from January to June there was a continued type II alert, except in March when it was a type I alert. In 2009–2010 there was a watch in November, and the damaged area was 72 %. However, as mentioned before, the absence of any further drought watches during that year and the high total annual rainfall indicate that the damage was likely to have been caused by excess precipitation. In the second main dry period of 2011–2012, the situation was worse, with a number of warnings from January to March, a type II alert in April and again a type I alert in May. That year the damaged insured area was high, up to 90 %.
The results led to the conclusion that the performance of the newly proposed CDI is adequate (Fig. 6). The periods of high crop damage – between 70 and 95 % – in the two important dry periods of 2004–2005 and 2011–2012 were accompanied by watches, warnings and type I or II alerts of CDI in the five agricultural regions studied. This combined indicator has several advantages over using a single one, as is evidenced by the trends in precipitation, soil moisture and vegetation alone. Soil moisture, for example, did not include the two main dry periods, 2004–2005 and 2011–2012, in all of the areas. The soil moisture anomaly index only indicated drought in two out of five regions for each of these dry periods, and this could probably be improved by measurements of in situ soil moisture. Krueger et al. (2017), for example, showed how in situ soil moisture measurements explained wildfire incidence much better than the widely used Keetch–Byram drought index (KBDI). Like our SMAI, the KBDI is a drought index calculated on a daily scale, but it only considers daily temperature and precipitation in calculating soil moisture. Whereas our SMAI uses a more advanced soil water balance algorithm (using variable infiltration rates and refining the estimation of the actual evapotranspiration rate from the potential rate computed by the FAO-Penman Monteith equation), and it is clear that future studies should focus on site-specific calibrations of soil moisture dynamics against field data or by observations from remote sensing. Martínez-Fernández et al. (2015) successfully applied in situ soil moisture measurements to predict agricultural droughts in northern Spain. Other studies, like that of Kędzior and Zawadzki (2017), have used SMOS-derived soil moisture anomalies. They concluded that these were suitable for calculating agricultural drought risk in the Vistula river catchment. Another possibility for improving drought prediction based on soil moisture is to combine different models. Cammalleri et al. (2016) used joint means from three different models, LISFLOOD, CLM and TESSEL, and were able to increase the correlation with observations and reduce the number of false drought alarms.
In any case, our results corroborate previous studies using combined indicators that also concluded that they yielded good results for agricultural drought prediction. Sepulcre-Canto et al. (2012), for example, use a similar CDI, based on SPI, soil moisture and photosynthetically active radiation (fAPAR). They evaluate this indicator on the continental scale and assess its performance against annual cereal yield at the regional level. They conclude that their indicator is successful in predicting drought periods and lower yields. While our indicator is similar in conception, there are notable differences with the CDI proposed in this study, firstly in the way soil moisture anomalies are calculated and secondly by using NDVI instead of fAPAR. Gouveia et al. (2009), comparing a soil water index against NDVI response in Portugal, found a good correlation between NDVI and soil water content under different land use conditions. They concluded that NDVI values of arable land were more sensitive to drought compared to forests, which suggests that NDVI is particularly well suited in this study of cereal growing areas.
Future studies could focus on improving this combined indicator, for example
by using other probability density functions rather than the gamma function
used for calculating the SPI. Sienz et al. (2012) obtained a better fit to
precipitation data of several world regions with the Weibull rather than with the
gamma probability distribution function. Carrão et al. (2016) selected
an empirical standardized soil moisture index, which was highly correlated
(
This study has presented a new combined drought index (CDI) for the assessment of agricultural drought. This CDI uses a combination of anomalies in precipitation (SPI-3), soil moisture and NDVI. The alert results are classified in four levels ranging from watch, warning to alert (type I and II). The CDI dynamics have been assessed for a 10-year period between 2003 and 2013, characterized by two important drought periods (2004–2005 and 2011–2012), in the five main rainfed cereal-growing regions of SW Spain. Comparison with yield data shows that both dry periods, characterized by a high crop damage extent of between 70 % and 95 %, were correctly identified by different critical CDI levels in all five study regions. This demonstrates the potential of this CDI. Further research should focus on a better representation of soil moisture data, either by improving data input from in situ measurements or by remote sensing, or by using model ensembles. Also, phenological information could be used to improve the performance of this indicator.
Data are freely available upon request by contacting the corresponding author by email.
JVG and AT conceptualized the research goals. MdPJD collected the remote sensing data and performed the data analysis. MdPJD developed the paper with contributions from all coauthors.
The authors declare that they have no conflict of interest.
This article is part of the special issue “Remote sensing, modelling-based hazard and risk assessment, and management of agro-forested ecosystems”. It is not associated with a conference.
María del Pilar Jiménez-Donaire was supported by a CEIGRAM fellowship. The insurance data were kindly supplied by Agroseguro. Ana Tarquis is grateful to the Community of Madrid (Spain) and Structural Funds 2014-2020 (ERDF and ESF) for the financial support (project AGRISOST-CM S2018/BAA-4330 and EU project 821964, BEACON).
This research has been supported by the European Commission, Horizon 2020 Framework Programme (BEACON (grant no. 821964)) and the Ministerio de Ciencia e Innovación (grant no. AGRISOST-CM S2018/BAA-4330).
This paper was edited by Jonathan Rizzi and reviewed by three anonymous referees.