Extreme weather events, such as droughts, have been increasingly affecting the agricultural sector, causing several socio-economic consequences. The growing economy requires improved assessments of drought-related impacts in agriculture, particularly under a climate that is getting drier and warmer. This work proposes a probabilistic model that is intended to contribute to the agricultural drought risk management in rainfed cropping systems. Our methodology is based on a bivariate copula approach using elliptical and Archimedean copulas, the application of which is quite recent in agrometeorological studies. In this work we use copulas to model joint probability distributions describing the amount of dependence between drought conditions and crop yield anomalies. Afterwards, we use the established copula models to simulate pairs of yield anomalies and drought hazard, preserving their dependence structure to further estimate the probability of crop loss. In the first step, we analyse the probability of crop loss without distinguishing the class of drought, and in the second step we compare the probability of crop loss under drought and non-drought conditions. The results indicate that, in general, Archimedean copulas provide the best statistical fits of the joint probability distributions, suggesting a dependence among extreme values of rainfed cereal yield anomalies and drought indicators. Moreover, the estimated conditional probabilities suggest that when drought conditions are below moderate thresholds, the risk of crop loss increases between 32.53 % (cluster 1) and 32.6 % (cluster 2) in the case of wheat and between 31.63 % (cluster 2) and 55.55 % (cluster 2) in the case of barley. From an operational point of view, the results aim to contribute to the decision-making process in agricultural practices.
Agriculture is one of the activities most directly influenced by climate variability (Lesk et al., 2016; IPCC, 2012) and by extreme weather events in particular (IPCC, 2012). The latter are a major source of risk in agricultural systems, often entailing substantial crop yield losses (Bokusheva et al., 2016; Kogan et al., 2015; Saadi et al., 2015). Despite the constraints associated with the application of certain governmental policies in agricultural systems, the success of mitigating the consequences of climate extremes is largely dependent on the development of appropriate risk management strategies (Paredes et al., 2014; Quiroga et al., 2011). For this purpose, probabilistic information of the agricultural risk associated with certain meteorological conditions is currently a major requirement (Madadgar et al., 2017; Jayanthi et al., 2014; Iglesias and Quiroga, 2007), particularly within the scope of projected climate changes (Hernández-Barrera and Rodríguez-Puebla, 2017; Ferrise et al., 2011).
From both researcher's and stakeholder's perspectives, the management of agricultural drought risk has been a challenging task for decades, mainly in regions dominated by high precipitation variability and recurrent dry and warm episodes, such as the Mediterranean region and the Iberian Peninsula (IP) in particular (Vicente-Serrano et al., 2014; Sousa et al., 2011; Martin-Vide and Lopez-Bustins, 2006). Recent work has found significant negative trends of drought indexes in the IP, based on long-term time series including the entire 20th century, particularly in southern regions (Páscoa et al., 2017a; Sousa et al., 2011), and the expected declining of crop yields due to future warming conditions has also been pointed out (Hernández-Barrera and Rodríguez-Puebla, 2017; Ferrise et al., 2011).
The assessment of yield variability based on crop and meteorological information is crucial for a more stable farmer income and management (Reidsma et al., 2010). The recently developed drought index SPEI (Standardized Precipitation Evapotranspiration Index; Vicente-Serrano et al., 2010) is found to be particularly suitable for agricultural drought applications in Mediterranean regions (Zampieri et al., 2017) and shows significant correlations with crop yields in the IP (Ribeiro et al., 2019a; Páscoa et al., 2017b). On the other hand, crop models describing the biological processes are one of the existing tools used to assess crop productivity, e.g. CERES (Crop Environment REsource Synthesis) models (Capa-Morocho et al., 2016; Hlavinka et al., 2010) and AquaCrop (Paredes et al., 2016; Vergni et al., 2015). These crop models are important tools in agrometeorological studies as they are able to compute irrigation requirements and yield simulations, and they have been particularly useful for assessing the impacts of climate change on agricultural productions (Leng and Hall, 2019; Hlavinka et al., 2010). However, such models are limited in their ability to quantify the impact of climate variability on crop yields over larger scales (Estes et al., 2013), and the detailed representation of crop's biophysical interactions requires demanding parameterization settings and input data (Giménez et al., 2016; Paredes et al., 2014, 2016). Thus, empirical modelling constitutes an alternative for representing the large-scale impacts of drought conditions in the agricultural sector (Bokusheva et al., 2016; Kogan et al., 2015; Matsumura et al., 2015; Vicente-Serrano et al., 2006), requiring lower computation costs than mechanistic modelling (Estes et al., 2013; Ferrise et al., 2011).
In addition, the use of satellite-based data is increasing for agricultural purposes (Kogan et al., 2015; Rojas et al., 2011) and considerable correlations between remote sensing of vegetation and crop yield are found in the IP (Ribeiro et al., 2019a; Gouveia and Trigo, 2008; Vicente-Serrano et al., 2006). Some studies have considered the use of different remote sensing drought indicators to account for different crop sensitivities to drought, such as to moisture and thermal conditions over the vegetative cycle (Ribeiro et al., 2019a; Bokusheva et al., 2016; Zarei et al., 2013; Kogan, 2001). Moreover, the establishment of models for estimating crop yield under drought influence, using the combination of different drought indicators and different timescales of drought occurrence, have shown an added value in the performance of the crop yield simulations over the IP (Ribeiro et al., 2019a; Hernandez-Barrera et al., 2017; Vicente-Serrano et al., 2006).
The statistical modelling of crop yield variability under drought conditions has been previously done to estimate drought-related crop losses (Ribeiro et al., 2019a; Zampieri et al., 2017; Kogan et al., 2015). Some authors have estimated crop yield probability distribution functions to find crop-specific risk levels and have applied Monte Carlo methods to generate large sample sizes of yield distributions over Mediterranean areas (Resco et al., 2010; Iglesias and Quiroga, 2007). At the country level in Europe, Naumann et al. (2015) have developed drought damage functions using a single power law dependence between drought severity and the associated damage. At a regional level in the IP, regression techniques (Ribeiro et al., 2019a; Hernandez-Barrera et al., 2017; Hernández-Barrera and Rodríguez-Puebla, 2017) and artificial neural network (ANN) models (Ribeiro et al., 2019a) have been used to model the response of rainfed winter cereal yields to drought conditions. A major conclusion in Ribeiro et al. (2019a) was that there are stronger relationships between remote sensing indices and cereal yield in the northern sector of the IP and between SPEI and cereal yield in the southern sector of the IP. This character of the response of crop yields to climate conditions highlights how it varies according to the location, type of crop, moment of the vegetative cycle, drought indicator and temporal scale.
More recently, copula-based models have been applied for agricultural purposes, to model the dependence structures between crop yields and environmental conditions using joint distributions (Ribeiro et al., 2019b; Madadgar et al., 2017; Bokusheva et al., 2016; Li et al., 2015). The concept of copulas is quite popular in financial risk modelling and has been becoming a valuable tool to model the risks associated with climate hazards, such as droughts (Ganguli and Reddy, 2012; Mirabbasi et al., 2012; Serinaldi et al., 2009). Based on the Sklar's theorem (Sklar, 1959) a copula approach “joins” the probability of drought occurrence and the probability of crop losses caused by the drought event. A detailed description about the use of copulas is provided by Nelsen (2006).
A major advantage of copula methods is the generation of joint distributions independently of their marginal distribution functions (Maity, 2018; Nelsen, 2006). Copula functions show a great flexibility in modelling the dependence between individual variables (such as crop yield and drought indicators) with complex relationships without making significant assumptions. In addition, copula functions are adequate for modelling rare events in multivariate distributions and generating large samples, allowing us to find the probability that individual variables will not exceed a certain extreme (tailed) value (Madadgar et al., 2017). A recent study by Madadgar et al. (2017) has produced probability distributions of rainfed crop yields in Australia under drought impacts based on copula-based techniques, using the Standardized Precipitation Index (SPI) and the Standardized Soil moisture Index (SSI). For crop insurance purposes at the farm level in Kazakhstan, Bokusheva et al. (2015) modelled the joint distributions of wheat yields and two satellite-based drought indices (vegetation condition index, VCI, and temperature condition index, TCI). At the global scale, Leng and Hall (2019) have also used copulas to assess the likelihood of yield loss in response to droughts based on SPI for the a historical (1961–2016) and a future period (2071–2100) under the RCP8.5 emission scenario to investigate future changes in yield loss risk. The authors found that global wheat is more vulnerable to droughts than maize, rice, and soybeans and that global warming is expected to amplify drought-driven yield loss risk.
In this study, a copula-based approach is adopted to model the joint probability density function of crop yield and the drought conditions for probabilistic yield assessment, based on the data and empirical analysis previously considered in Ribeiro et al. (2019a). This method allows us to estimate the dependence structures between the probability distributions of crop yield and drought indicators using copula functions. The novelty and interest of this approach relates to the fact that this methodology will allow us to estimate the likelihood of crop loss and compare the expected losses under drought conditions and non-drought conditions in the IP. This key question is posed based on the current demand, of the most interest to stakeholders such as farmers and insurance companies, to mitigate agricultural drought risk over the major agricultural areas in the IP.
The exposure analysis performed by Ribeiro et al. (2019a) allowed the
identification of two clusters of provinces in the IP dominated by rainfed
agricultural practices (Fig. 1), located approximately in the regions of
Castilla and Léon (cluster 1 – northern region) and Castilla–La Mancha
(cluster 2 – southern region). Given the suitability of using these two
clusters for an agricultural drought analysis at the regional level, here we
have considered the same area selection criteria: provinces with more than
50 % of the territory occupied by agricultural areas and more than 50 %
of rainfed crops according to the CORINE Land Cover (2012) (for more details please see Ribeiro et al., 2019a). Considering previous requirements, and for sequential purposes, the crop and drought hazard data used in Ribeiro et al. (2019a) have been incorporated in the present study to analyse the distributions of probabilities. Spatial averages of annual yield anomalies (t ha
Provinces with more than 50 % of the territory occupied by agricultural areas and more than 50 % of rainfed crops (yellow), according to CLC (2012), and selected clusters of provinces. Cluster 1 provinces: Zamora (I), Valladolid (II), Palencia (III), Burgos (IV) and Segovia (V). Cluster 2 provinces: Toledo (VI), Cuenca (VII), Ciudad Real (VIII) and Albacete (IX).
Drought conditions were investigated using the hydro-meteorological drought
indicator SPEI and three satellite-based indices obtained from NOAA-AVHRR
since 1981, namely the VCI (Kogan, 1990), the TCI (Kogan, 1995) and the Vegetation Health Index (VHI) (Kogan, 1995). The monthly drought index SPEI gridded values, with a spatial resolution of 0.5
Considering the vegetative cycle of wheat and barley, and in accordance with the results obtained by Ribeiro et al. (2019a), the data of VCI, TCI and VHI used in this work covered the period from week 35 (early September) to week 25 (late June) and data of SPEI covered January to June. Spatial averages of all of these indicators were computed for each provincial cluster and used for further modelling of the joint probability between the drought hazard and cereal yield anomalies over the period 1986–2012. Stepwise regression models (95 % confidence level) were established to select the timescales and months of SPEI, together with the weeks of VCI, TCI and VHI, better related with wheat and barley annual yield (Ribeiro et al., 2019a). The selection of the most relevant drought indicator for each cereal and cluster was performed based on the largest absolute value of the standardized regression coefficients from the models developed in Ribeiro et al. (2019a), in order to constitute pairs of cereal yield anomalies and drought indicators. Afterwards, for each cereal time series, the joint probability with drought conditions was estimated using one drought indicator.
Copula functions are powerful tools used to estimate the joint distribution
between variables (Madadgar et al., 2017; Bokusheva et al., 2016; Zhang et al., 2011). The concept of copula was firstly introduced by Sklar (1959) to decompose a joint cumulative distribution function
There is a range of copula families described in the literature that are
able to estimate the dependence between the univariate variables (Nelsen, 2006). The most commonly used copula families focus on the Archimedean and elliptical classes (Maity, 2018). There are three Archimedean copulas
that are particularly popular, given their simple functional form and their different
patterns of dependence captures, i.e. Clayton, Gumbel and Frank, while there
the two most popular elliptical copulas are derived from elliptical distributions, i.e. Gaussian and
Equations of the copula functions, where
An important concept for studying extreme events is the tail dependence,
whose importance is more critical than the overall dependence structure for
risk analysis (Bokusheva, 2014). The joint tail behaviour describes the amount of dependence in the corners of upper-right and lower-left quadrants (i.e. joint extreme events) and its representation depends on the type of copula (Nelsen, 2006). The Frank, Gaussian and
The estimation of the copula parameters can be performed using different methods based on maximum likelihood, such as maximum likelihood estimate (MLE), inference functions for margins (IFM) or canonical maximum likelihood (CML) (Maity, 2018). With MLE, both individual margins and copula parameters are estimated together, whereas with IFM the marginal parameters are first estimated individually. In this study the statistical inference of the copula functions is performed with the CML method, which stands for a nonparametric estimation of the margins. In this way, the individual variables were first transformed to the unit scale (pseudo-observations) using the kernel density estimator of the cumulated distribution function (CDF) without making assumptions about the marginal distributions (Fig. 2). The drawback of the shorter sample size is surpassed by the nonparametric estimation of the margins, which avoids significant assumptions about their distributions, even when the available sample is rather small (Fahr, 2017; Corder and Foreman, 2011). The fitting of the bivariate copula functions was then applied to the pseudo-observations, and the dependence parameters were estimated by means of maximum likelihood (Fig. 2). Figure 2 summarizes the main steps of the copula-based approach adopted in the present study. For a detailed description on fitting methods please see Maity (2018).
Scheme of the copula-based approach adopted in the present study.
Akaike's Information Criterion (AIC) is frequently employed as a model
selection tool in copula modelling (Li et al., 2015; Mirabbasi et al., 2012). Therefore, the selection of the best copula function for each pair of cereal and drought indicators was made based on the evaluation of AIC values, calculated as AIC
After the estimation of the copula parameters, the established models are
used to simulate 1000 pairs of uniformly distributed data (Fig. 2). In the
present study, let
First the copula simulations were used to estimate the risk of crop loss in
terms of the probability of not exceeding a threshold value of yield, i.e.
probability of non-exceedance (PNE) (Fig. 2). In this study we considered
the threshold of minus one standard deviation (
Afterwards we partitioned the simulated data points of
In sum, first we describe the joint probability of drought hazard and yield anomalies and simulate pairs of data preserving their dependence structure. After that, probability of crop loss (PNE) and conditional probability of crop loss (CPNE) are estimated, addressing whether the probability of crop loss under drought conditions is higher than during non-drought conditions and if distinguishing drought severity is important. The probability distributions (based on a normal kernel function) of the generated yield anomalies are also analysed for graphical visualization of the area corresponding to crop loss.
The estimates of the dependence between the yield anomalies and drought hazard were performed using the selected drought indicators outlined in Table 2. This selection of drought indicators highlights that the response of crop yields to climate conditions varies according to the location, type of crop, moment of the vegetative cycle and chosen temporal scale. While annual yield anomalies in cluster 1 are better characterized by short-term responses to the drought conditions based on the weekly values of TCI and VHI, the annual yield anomalies in cluster 2 are better characterized by the monthly response to the dry conditions based on the SPEI. In terms of predictability, the effects of temperature (TCI) and vegetation health (VHI) during late growth stages (weeks 23 and 22 correspond approximately to end of May and beginning of June, respectively, for wheat and barley) are the most influential conditions in the northern cluster. On the other hand, the yields in cluster 2 are influenced by drought conditions described by SPEI much earlier, in the beginning of the intermediate growth stages (February and April with 5 and 1 month timescales, respectively, for wheat and barley). In this way, the importance of including multiple drought response timescales is evidenced for predictability purposes and assessment of drought-related crop losses.
Variables used for copula application. In the first column, the numbers 1 and 2 correspond to the respective provincial cluster (clusters 1 and 2). In the second column, the numbers correspond to the selected weeks in the case of the remote sensing indices and to the selected months and timescales (in months) in the case of SPEI. The values of the standardized regression coefficients were determined by Ribeiro et al. (2019a).
Empirical cumulative distribution functions (ECDF, blue points), kernel density estimation of the CDF (red line), crop loss and drought thresholds (dotted black vertical line), respective marginal probabilities of crop loss and drought occurrence (dotted black horizontal line), and pseudo-observations (scatter) of the margins on the interval [0, 1].
Figure 3 shows the non-parametric estimations of the CDF of the individual
variables from Table 2, here used to transform the variables to the unit
scale (pseudo-observations) for the copula modelling. A good agreement with
the ECDF is suggested (Fig. 3) and the crop loss and drought thresholds used
in this study (
The estimates of the dependence between the yield anomalies and drought
indicators were performed using the copula functions from Table 1 (Gaussian,
Copula dependence parameter estimates (
The selected copula functions (Table 3) suggest that, in general, the relationship between yield and drought conditions is described by an asymmetric dependence in the tails of the joint distributions, except in the case of barley in cluster 1. This feature is illustrated in Fig. 4, showing the different shapes and contours of the selected copula densities. While wheat in cluster 1 and 2 shows a stronger dependence in the upper tail of the joint distributions based on Gumbel copulas (suggesting higher probability of observing a higher value of yield anomalies given a high value of the drought indicators), barley in cluster 2 shows stronger dependence in the lower-left tail based on a Clayton copula, suggesting higher probability of finding a lower value of yield anomalies given a low value of the drought indicators. The randomly generated yield and drought data were transformed back to the original scales (Fig. 4, bottom row) and the respective scatter plots indicate that more extreme values are generated using the joint distribution relationships. In general, the modelling of the joint distributions leads to results close to the real observations (Fig. 4, bottom panel).
Selected joint probability distribution functions (PDFs) where
After estimating the joint distribution functions and simulating pairs of data preserving the modelled dependence structures, we evaluate and compare the probability of non-exceedance (PNE) and conditional probability of non-exceedance (CPNE) as a function of the crop loss threshold. In this way, we evaluate if distinguishing drought severity leads to different risk values of crop loss in comparison to disregarding a drought threshold (using only simulations of yield) and compare the probability of crop loss under drought and non-drought conditions (by means of both simulations of yield and respective drought indicator). One of the key advantages of estimating the values of PNE and CPNE by means of the copula simulations is the use of larger samples that are comprised of more joint extreme values based on the joint behaviour of crop yields and drought hazard.
Probability of non-exceedance (PNE) function (%) of yield anomalies (top row) in both clusters based on the derived simulations from the estimated copulas and respective probability density estimates (bottom row). In the bottom row, the red values indicate the probability of crop loss, which is also indicated in the top row by the intersected dashed lines, indicating the threshold of crop loss and respective PNE value.
Figure 5 shows the PNE curves and the distributions of the simulations of
yield anomalies, with the respective crop loss area correspondent to the
probability (%) of the yield anomaly not exceeding
The following target was to compare the likelihood of crop loss under drought and non-drought conditions. Figure 6 shows the simulated crop yield anomalies during drought (orange boxplots) and non-drought (blue boxplots) events. As expected, the boxplots show lower (and negative in average) values of yield anomalies during drought events in comparison with non-drought episodes. Although the number of samples simulated under drought conditions is smaller than under non-drought conditions (Fig. 6), the use of copula simulations enhances the amount of simulated joint low extremes (i.e. co-occurrence of crop loss and drought events).
Wheat and barley yield simulations, differentiating drought (orange) and non-drought conditions (blue) according to the respective drought indicator denoted in parenthesis in the
The differences in terms of crop losses between cereals and regions is much evident when differentiating the climatic conditions (Fig. 7), particularly during drought conditions. Figure 7 shows that the values of CPNE under drought (non-drought) conditions are above (below) the values of PNE illustrated in Fig. 5. In comparison with the distributions of yield simulations without conditioning to specific thresholds of the drought indicators shown in Fig. 5, in Fig. 7 the distributions of the yield simulations during drought events show a shift to the left towards negative values of yield anomalies, while the distributions of yield simulations during non-drought events show a shift to the right towards positive values of yield anomalies (Fig. 7). The case of barley in cluster 1 is quite distinct exhibiting a drought-related barley loss almost 3 times higher than the value illustrated in Fig. 5 (19.8 %), supporting the importance of conditional probabilities for agricultural drought risk purposes. The conditional probability of wheat loss (Fig. 7), is also higher when focusing on drought conditions, although it is less than 2 times the values shown in Fig. 5.
Conditional probability of non-exceedance (CPNE) function (%) based on the derived copula simulations (top row) and respective probability density estimates (bottom row) under drought (orange) and non-drought conditions (blue). In the bottom row, the orange and blue values indicate the probability of crop loss under the different climatic conditions, which is also indicated in the top row by the intersected dashed lines, which indicate the threshold of crop loss and respective CPNE value.
Regarding the drought-related barley loss, the distribution of barley in cluster 1 is more shifted to negative yield anomalies, stressing that the drought risk of barley loss is higher on cluster 1 (59.2 %) than on cluster 2 (39.4 %), while it is quite similar on both clusters in Fig. 5. While barley suggests higher conditional probabilities of crop loss under drought conditions in cluster 1, wheat suggests higher conditional probabilities of crop loss under drought conditions in cluster 2 (46.7 %) in comparison to cluster 1 (36.5 %). Among all the cases, the highest level of drought-related crop loss is 59.2 %, observed in the case of barley in cluster 1, followed by wheat in cluster 2 with 46.7 % chance of crop loss under dry conditions.
The theoretical CPNE based on Eqs. (5) and (6) (Table 4) agrees quite well with the estimates of the CPNE in Fig. 7, thus corroborating the
representativeness of the copula experiment using 1000 simulations.
Nevertheless, the use of simulations allows us to increase the sample size and
to generate more joint extreme values based on the dependence structures
characterized by the selected copulas. In addition, the effect of the copula
parameter (
Theoretical CPNE (%) during drought and non-drought conditions (Eqs. 9 and 10) and respective lower and upper bounds of the 95 % confidence interval, where
The results show that CPNE based on simulations (Fig. 7) and theoretical
equations (Table 4) indicate that the probabilities of crop loss increase
when drought conditions occur, even considering the two-sided confidence
bound values of the copula parameters. Moreover, the results indicate that
the CPNE using the simulations (Fig. 7) lies within the estimates of CPNE
using the two-sided confidence bound values of the copula parameter at the
95 % level of confidence (Table 4). The only exception is the case of
barley in cluster 2 considering the lower bound of
This study investigated the usefulness of copula methods in estimating the likelihood of drought risk in wheat and barley cropping systems, when applied to two regions in the IP. Here we proposed modelling the joint probability of yield and drought hazard using copulas, based on a prior analysis of the association between drought and crop loss (Ribeiro et al., 2019a). The advantage of using a probabilistic approach is to meet the ambitious challenge of helping farmers and stakeholders in managing their operations by identifying the probability of crop loss under specific drought conditions. Hernández-Barrera and Rodríguez-Puebla (2017) and Ferrise et al. (2011) have shown that the projected warmer and drier climate will lead to wheat yield shortfall over the IP and Mediterranean, respectively, highlighting the importance of establishing novel statistical approaches for agricultural drought risk analysis. Other crops rather than rainfed cereals are also expecting significant losses during the next century in the IP (Saadi et al., 2015; Resco et al., 2010; Quiroga and Iglesias, 2009), and the here-proposed crop-specific approach could be applied to other agricultural systems under drought conditions for different regions around the world.
The novelty of the presented models, in comparison to other works addressing climate risk in the IP (e.g. Ribeiro et al., 2019a; Resco et al., 2010; Iglesias and Quiroga, 2007), is the focus on the impacts associated with droughts and on the joint probability of rainfed yield anomalies and drought hazards. Previous works using copulas in hydro-climatology studies have tended to focus on the joint distribution of different characteristics of the hazardous events, such as frequency, intensity, severity, and duration, among others (Li et al., 2015; Chen et al., 2013; Mirabbasi et al., 2012). Moreover, the restriction to the bivariate case allowed for a simpler interpretation of the results, in contrast to higher-dimension copulas (Afshar et al., 2016; Ganguli and Reddy, 2013), for instance by adding other factors influencing crop yield beyond drought as copula variables.
More recently, copulas have been applied to estimate the joint behaviour of drought conditions and the associated impacts in agricultural systems (Leng and Hall, 2019; Ribeiro et al., 2019b; Madadgar et al., 2017; Bokusheva et al., 2016), instead of using drought information only. We have adopted a similar approach to reproduce time-, regional- and crop-specific dependence of drought conditions, and the probability distribution of crop yield anomalies under drought conditions was estimated for risk analysis. In addition, the use of different drought indicators in this study represents an advantage since crops react differently to several factors at distinct moments and locations, highlighting the importance of quantifying the contributions of different drought indices on a regional scale (Peña-Gallardo et al., 2019; Zarei et al., 2013). A recent study by Peña-Gallardo et al. (2019) focused on the responses of wheat and barley cropping systems to different drought indices over Spain, have shown the different efficacies of several drought indices, stressing the importance of the multiscalar character of droughts, in particular of the SPEI. Similarly, and in accordance to previous work by the authors (Ribeiro et al., 2019a), the present study shows the adequacy of SPEI for the assessment of the agricultural risks associated with droughts in the IP and advances the added value of using the remote sensing of vegetation.
Overall, the results of the estimated copula functions have shown that Archimedean copulas are suitable to model the joint behaviour of yield anomalies and droughts, suggesting a dependence between extreme values of rainfed cereal yield anomalies and drought indicators, and the subsequent simulated distributions of crop yield anomalies are quite consistent with the observations. The results highlighted that the use of copulas for probabilistic assessment allows the estimation of the dependence in the tails of the distribution and was able to give the likelihood of crop loss under drought conditions. This feature is of the most interest in risk analysis given that it models the joint probability of occurrence of crop loss and drier events. Moreover, this study suggests the relevance of impact-centric approaches (also referred to the literature as “bottom-up” approaches; Zscheischler et al., 2018) to identify and characterize the hazards that lead to the larger impacts.
Moreover, it is important to stress that crop anomalies decline much more when drought conditions are below the mild or moderate threshold, suggesting a high agricultural drought risk level of wheat and barley in both clusters. While values of PNE in the crop loss threshold were low and similar for wheat in cluster 1 and barley in cluster 2, the values of CPNE in the crop loss threshold during drought years are considerably larger. The higher probability of crop loss obtained when analysing only drought conditions agrees with Páscoa et al. (2017b), which show a very high agreement between low wheat yield anomalies and drought conditions in the IP, even on provinces where the linear correlation is not significant.
Although there is a greater risk of crop loss during drought conditions, some losses can still be expected during non-drought events, particularly in cluster 2 (14.1 % and 7.77 % in the cases of wheat and barley, respectively). In the northern sector (cluster 1) the probabilities of crop loss under non-drought conditions have the lower values, displaying 3.97 % in the case of wheat and 3.65 % in the case of barley. Some studies point to crop damage attributable to excessively wet soils (Zampieri et al., 2017; Rosenzweig et al., 2002), due to delayed planting or later harvest, nutrient runoff and development of pests and diseases, among other factors, highlighting the complexity of quantifying agricultural risk levels for management purposes and the non-linear relation between crop yield and climate conditions. The lower values of CPNE under non-drought conditions in cluster 1 support the fact that the slightly high values of PNE in cluster 1 are mainly dominated by drought conditions.
With the present study it is not possible to establish sharp conclusions about the adequacy of the copula models to a specific type of drought indicator (remote sensing or hydro-meteorological), since only one type of drought indicator was considered for each cereal. In contrast, Bokusheva et al. (2016) have found that Gumbel copulas provided better fits representing the joint distribution of VCI and wheat, while Frank copulas better described the dependence between TCI and wheat yields in Kazakhstan. Madadgar et al. (2017) modelled the conditional probability density functions of crop yields under wet and dry conditions using SPI and SSI and found that a Clayton copula was the best function to model the dependence structures. Similarly, Leng and Hall (2019) have also used the same copula families and found that from 10 countries studied 5 used Clayton copulas to fit the joint distribution between wheat production and SPI. However, the referred studies were somehow more restrictive as they do not take advantage of using both remote sensing and hydro-meteorological drought indicators and do not select the most important one a priori.
To further the research, the application of SPEI methodology to climate projections of precipitation and temperature holds an added value to the estimation of drought risk levels for the next century. Likewise, the use of seasonal drought forecasts is also quite plausible in the approach presented in this study. Nevertheless, the presented results indicated the likelihood of crop loss based on drought conditions observed much earlier than the harvest time, particularly in cluster 2 using SPEI (February and April with 5 and 1 month timescales). Hence, given the uncertainty associated with the seasonal forecasts for regional drought predictability in the IP, the use of past information for predictability studies is still successfully used (Pires and Ribeiro, 2016) and continues to be a source of information from an operational point of view. Another potential use of this methodology for future research is the evaluation of its suitability at the province level and the assessment of whether other hazards (such as heat waves) are amplifying the impact of droughts on crop harvests.
The agricultural drought risk levels estimated in the present work aimed to
improve the effectiveness of the agricultural management of rainfed cereals
in the major agricultural areas of the IP. The main findings of this study
are summarized below.
The dependence structure between crop yield anomalies and drought conditions is mainly asymmetrical, suggesting the existence of dependence among extreme values of yield anomalies and drought indicators. The differences between the unconditional and the conditional probability suggest that the risk of wheat loss and barley loss can be underestimated without conditioning the probabilities of non-exceedance crop thresholds to specific drought levels. The conditional probabilities of non-exceedance suggest that the risk of wheat loss and barley loss increases when drought events aggravate from normal or wet to moderate or severe conditions. The values of conditional probabilities of crop loss under dry conditions suggest that the risk of drought-related barley loss is more likely to occur in the northern sector, while the risk of drought-related wheat loss is more likely in the southern sector, suggesting that sowing in cluster 1 (cluster 2) could be more focused on wheat (barley). The overall results show the importance of the concept of conditional
probability for distinguishing different meteorological settings associated with
crop losses and the applicability of the copula theory. The use of copula
simulations for the analysis of the co-occurrence of dry and low-yield
extreme events has shown the additional value of this methodology for the
estimation of drought-related crop failure. Nevertheless, minor wheat and barley losses can still be expected during normal or wet conditions, stressing the complexity of the interactions between the agricultural systems and the climate. Particularly under the current climate change context, further high-impact-centric analysis are required, involving the cascading effects of different climate hazards.
The Portuguese crop yield database is published by the Portuguese National Statistics Institute and is available at
AFSR, AR and CMG conceived the idea of the present study. AFSR analysed the data, performed the statistical analysis, produced the figures and drafted the manuscript. AR and CMG supervised the work. CALP verified the analytical methods and designed the estimation of the conditional probabilities' uncertainty. AR performed the computations of SPEI. PP instructed the acquisition and analysis of the crop yield data. CMG instructed the acquisition and analysis of the remote sensing data. All the authors provided helpful insight in the discussion of the results and contributed to the design of the research and to the final paper.
The authors declare that they have no conflict of interest.
This article is part of the special issue “Hydroclimatic extremes and impacts at catchment to regional scales”. It is not associated with a conference.
This work was partially supported by national funds through the FCT (Fundação para a Ciência e a Tecnologia, Portugal) under the projects CLMALERT (ERA4CS/0005/2016) and IMPECAF (PTDC/CTA-CLI/28902/2017). Ana Russo and Andreia Ribeiro thank FCT for grants SFRH/BPD/99757/2014 and PD/BD/114481/2016, respectively. The authors are also sincerely thankful to Ricardo Trigo (Instituto Dom Luiz) for his valuable suggestions and support.
This research has been supported by the FCT (Fundação para a Ciência e a Tecnologia, Portugal) (grant nos. SFRH/BPD/99757/2014, PD/BD/114481/2016, ERA4CS/0005/2016, WaterJPI/0004/2014 and PTDC/CTA-CLI/28902/2017).
This paper was edited by Fernando Domínguez-Castro and reviewed by three anonymous referees.