Mediterranean-Scale Drought: Regional Datasets for Exceptional Meteorological Drought Events during 1975–2019

: Drought is one of the most complex climate-related phenomena and is expected to pro-gressively affect our lives by causing very serious environmental and socioeconomic damage by the end of the 21st century. In this study, we have extracted a dataset of exceptional meteorological drought events between 1975 and 2019 at the country and subregional scales. Each drought event was described by its start and end date, intensity, severity, duration, areal extent, peak month and peak area. To define such drought events and their characteristics, separate analyses based on three drought indices were performed at 12-month timescale: the Standardized Precipitation Index (SPI), the Standardized Precipitation Evapotranspiration Index (SPEI), and the Reconnaissance Drought Index (RDI). A multivariate combined drought index (DXI) was developed by merging the previous three indices for more understanding of droughts’ features at the country and subregional levels. Principal component analysis (PCA) was used to identify five different drought subregions based on DXI-12 values for 312 Mediterranean stations and a new special score was defined to classify the multi-subregional exceptional drought events across the Mediterranean Basin (MED). The results indicated that extensive drought events occurred more frequently since the late 1990s, showing several drought hotspots in the last decades in the southeastern Mediterranean and northwest Africa. In addition, the results showed that the most severe events were more detected when more than single drought index was used. The highest percentage area under drought was also observed through combining the variations of three drought indices. Furthermore, the drought area in both dry and humid areas in the MED has also experienced a remarkable increase since the late 1990s. Based on a comparison of the drought events during the two periods—1975–1996 and 1997–2019— we find that the current dry conditions in the MED are more severe, intense, and frequent than the earlier period; moreover, the strongest dry conditions occurred in last two decades. The SPEI-12 and RDI-12 have a higher capacity in providing a more comprehensive description of the dry conditions because of the inclusion of temperature or atmospheric evaporative demand in their scheme. A complex range of atmospheric circulation patterns, particularly the Western Mediterranean Oscillation (WeMO) and East Atlantic/West Russia (EATL/WRUS), appear to play an important role in severe, intense and region-wide droughts, including the two most severe droughts, 1999–2001 and 2007–2012, with lesser influence of the NAO, ULMO and SCAND.

thereby, sustainable human and environmental resource development at increased risk [1]. The intensive climate extremes are projected to result in a change in water availability, exacerbate carbon losses in terrestrial ecosystems, and further raise water vapor content in the atmosphere, thus amplifying the warming impacts and increasing mega-droughts worldwide [2][3][4]. Drought is one such recurring, widespread, complicated, costliest and disruptive climate extreme that frequently occurs as a result of differential responses to climate warming in different climatic regimes [5][6][7]. Recently, drought has begun to be seen not as a purely natural hazard, but partially as a result of human action, which has altered many of the characteristics of drought [1,8]. One-fifth of the global destruction caused by natural disasters can be attributed to droughts [9], where the dramatic population and economic growth has resulted in increasing demand for water that substantially intensified the frequency of global drought events [10]. Pandemics can potentially exacerbate drought impacts [11], as it was found that drought is among the most obvious reasons for displacement and disturbance during the coronavirus pandemic [12,13]. Droughts in the 21st century are identified as multifaceted, challenging natural disasters [14,15] and are characterised by longer duration, higher severity, larger spatial extent and hotter temperature that can potentially result in many deleterious impacts on ecological security, with non-linear alterations in ecosystem functions and resilience [3,4]. Systematic reviews [16,17] have indicated an accelerated transition to a more arid climate over several areas as the result of an increased tendency of frequent and intense droughts. Meanwhile, human-induced global warming and rapidly expanding human populations have already influenced water availability and storage, increasing pressure on water supplies, affecting the long-term ecosystem and increasing their sensitivity to droughts [4,18]. As a result of 1.5 to 2 °C global warming, two-thirds of the world population will experience increasing droughts; the drought magnitude is likely to double in 30% of the global land by the end of this century [19,20]. An important part of drought problem is that its current definitions refer to the drought only from the perspective of human dimensions, focusing primarily on meteorological, agricultural, socioeconomic and cultural impacts without addressing its ecological dimensions [21]. The occurrence of widespread drought in developing and developed societies has underscored the sensitivity of all communities to this natural hazard. It is not easy to know whether the drought frequency is increasing, or rather the community exposure to it. This tendency appears to be accelerating as a result of the increasing demand on both local and regional water resources [22]. The socioeconomically devastating impacts of frequent drought events have recently resulted in several global assessments of future drought conditions to better support populations and improve management plans in order to reduce direct and indirect cascading drought impacts [23][24][25]. Combining the natural and human dimensions of drought is one of the most fundamental steps in addressing the increased risk of drought in the 21st century [18,26].

Drought, Climate Change and the Mediterranean Context
The MED is identified as one of the most responsive hotspots of the climate system in the face of increased dryness at global warming levels beyond 1.5 °C [27]. The IPCC 1.5 degree Special Report (SR15) has pointed out that the increased anthropogenic warming in the MED has contributed to increased drying in the northern hemisphere mid-latitude areas including the MED [28,29]. Furthermore, with global warming of 2 °C, desertification is also predicted to occur in the MED by the end of the 21st century, resulting in an expansion of areas with significant decreases in water availability, with an accompanying increase in aridity, driving irreversible terrestrial biodiversity loss and affecting the Mediterranean ecosystem carbon storage in the coming decades [30]. The MED experienced more frequent and severe meteorological and hydrological droughts in recent decades, which is in line with the expected trend towards high frequencies of drought periods in a future warmer climate [31,32]. Above 2 °C, the MED could become more vulnerable to drought, calling for the development of a variety of adaptation mechanisms and the pursuit of drastic adaptive responses to cope with these extreme climate events, including mitigation strategies in addition to radical changes in the social structure and human communities [33,34]. The complexity of the Mediterranean climate, with its high rainfall variability and its unequal seasonal distribution in addition to the observed dynamical and physical atmospheric process complexities [35,36], produced the conditions that have led to the high sensibility and vulnerability of this region to droughts [37]. Droughts are not spatially coherent in the MED [38] demonstrating different spatial patterns even at the regional scale [39,40].
The main goal of our study is the identify and characterise meteorological drought events over the MED by constructing a robust list of the most relevant drought events that occurred in the MED between 1975 and 2019, and analysing these events. The severity, intensity, spatial extent, peak month, area involved at peak month and the frequency of drought occurrence are the essential characteristics that we have investigated. Compared to the existing drought datasets, we have provided some important developments and novelties that include a new and detailed database of drought episodes over the MED. The importance of this dataset is that it deals in details with the drought events in the Eastern Mediterranean and Middle East (EMME), where drought studies are still relatively few. Based on quality-checked and homogenised data, three drought indices were applied in order to obtain detailed information and better understanding of drought events' features over the MED in addition to formulating a strong dataset of these events, taking into account the impacts of both precipitation (PP) and temperature (TEMP). The obtained drought events records have specific entries at both the country and subregional scales, which can then be used to extract the main drought characteristics for each subregion.
The remainder of this paper is structured into five main sections. In Section 2, we present the input data, study area and the selected drought indices. The main characteristics of different drought events and the new applied system for classifying the drought events are also demonstrated in this part. Section 3 displays the analysis results and discussion and describes, separately: the database structure; the collected records at the country and subregional scales; the biggest constructed drought events; and the drought dynamics in terms of large-scale teleconnection patterns. Finally, Section 4 derives some conclusions from the main results and includes some critical research questions for evaluating future studies of drought in the region.

Study Area and Data
The Mediterranean basin, with latitude boundaries between 30° and 49° N and longitude boundaries between 10° W and 45° E, is the study area in this work. This study area reaches portions of three continents: Africa, Asia, and Europe. Climate is very variable within this area and characterised by warm to hot, dry summers and mild to cool, wet winters [40]. The daily meteorological datasets  used in this study are compiled by both the European Climate Assessment & Dataset (ECA & D; [41]) and the Global Historical Climatology Network (GHCN; [42]) Data for Syria, Lebanon, Algeria, and Jordan are directly obtained from national meteorological and hydrological services and some regional research projects. Main statistical characteristics of total annual precipitation and mean annual temperature for all stations are shown in Appendix A. After an initial quality check, the mean monthly maximum and minimum temperature (Tmax, °C), (Tmin, °C) were computed as an arithmetic average of daily maximum and minimum temperatures, respectively. The mean temperature values for any given month were considered if there were no more three missing values in this month. Precipitation amounts (P) were transformed from daily values into monthly sums if no more than two daily values were missing. Monthly temperature and precipitation series have also undergone further quality checks and data errors and inhomogeneities are eliminated using the MASH V3.03 (Multiple Analysis of Series for Homogenization [43,44]) homogenization procedures and 312 stations were finally selected for our study. About 20% of total stations failed the homogeneity tests, more than 90% of which are located in Greece, Italy, Algeria and Turkey. Depending on these results, we discarded 78 stations and completed this study with 312 meteorological stations distributed evenly located at the west-east axis of the Mediterranean basin ( Figure 1). The values NOA, EATL/WRUS and SCAND indices used in the study were obtained through the NOAA CPC website: https://www.cpc.ncep.noaa.gov/data/teledoc/telecontents.shtml, accessed on 21 July 2021 while the following websites have been used to obtain the values of the WeMO and ULMO, respectively: https://crudata.uea.ac.uk/cru/data/moi/, http://www.ub.edu/gc/data_base/, accessed on 21 July 2021.

Selection of Drought Indicators: The SPI, SPEI and RDI
Drought indices are important and effective elements for providing a comprehensive picture of drought to plan, monitor, and quantify drought and its risk management [45,46]. Multiscale drought meteorological drought indices are commonly used as proxies for characterising soil moisture at different layers [47,48]. However, due to drought complexity, using a single index is not sufficient for precisely identifying the drought dimensions [49], and the combined use of several different drought indices would give an actual and clearer picture of drought conditions [50]. Drought is defined for each station using the Standardized Precipitation Index (SPI; McKee et al. [51][52][53]), which is based on precipitation data only and measures the normalised anomalies in accumulated precipitation through a given number of months. The SPI was recommended by the World Meteorological Organization (WMO) in 2012 to characterise the meteorological droughts and quantify the precipitation deficit at multiple time scales [54]. The later Standardized Precipitation Evapotranspiration Index (SPEI; [55,56]) is also applied. The SPEI is designed to represent climatic water balance (the difference between precipitation and potential evapotranspiration, PET), which captures the main effects on water requirement as a result of evapotranspiration. Using the ratio between cumulative precipitation and PET, the Reconnaissance Drought Index (RDI; [57,58]) is also used. The SPEI and RDI provide a more comprehensive representation of drought characteristics due to their inclusion of precipitation along with potential evapotranspiration that can be effectively used to compare the drought events under different climate conditions [59,60]. Droughts in the MED have very different dimensions and affect several sectors, but currently, the most important impacts relate to the hydrological one, which is largely affected by frequent droughts [61]. In this paper, the drought indices values are generated using a 12-month aggregation period because of its great importance in estimating the severity of hydrological droughts, determining the possibility of drought persistence conditions for long periods of time [62] and considering the management of water supplies [40,63]. This paper focuses on long term drought indices (12-month timescale) because of their great importance in estimating the severity of hydrological droughts and determining the possibility of drought persistence conditions for long periods of time [62]. The input data are fitted to a two-parameter gamma probability density function for SPI-12 and RDI-12, whereas the SPEI-12 follows the log-logistic probability distributions. The PET is derived using the Hargreaves-Samani equation [64,65] which indirectly estimates extra-terrestrial radiation using Tmax and Tmin. Although the estimated evapotranspiration using temperature-based models may magnify the impacts of the current global warming [66], it has been indicated that drought indices provide reliable estimates of the true global drought trend when using the temperature-based models [67]. The temperature-based approaches, such as Thornthwaite, have produced reliable results [68]. The Hargreaves-Samani outputs show high correlation with Penman-Monteith method [56] and confirm high reliability in arid, semi-arid, Mediterranean, and very humid areas [69]. SPI-12 and SPEI-12 values are obtained using the SPEI package developed by Beguería and Vicente-Serrano (2017), available in R (Version.3.6.3, Development Core Team 2020) [70].

Drought Concepts and Characterisation, a Multivariate Perspective for a New Combined Index
A single index usually considers particular processes while ignoring others, leading to it not capturing an accurate picture of hydrological changes, and resulting in an incorrect interpretation [71]. To overcome the single-index problem, the three drought indices were integrated into one. For each station, a new combined drought index, DXI-12, was constructed by averaging the monthly values of the SPI-12, the SPEI-12, and the RDI-12 as a multivariate drought index to incorporate multiple drought-related values. The combination of multiple drought-related variables was successfully applied for drought monitoring [72], detecting [73,74], and prediction [75,76]. For example, results based on [77] have suggested a combined use of meteorological drought indices (SPI and SPEI) in the case of lack of streamflow measurements to provide the best performances in monitoring hydrological droughts. The development of a combined drought indicator gives a synoptic and synthetic overview and provides an integrated approach that fully describes the meteorological drought conditions across the MED [72]. The selected time scale relates to the time required for the effects of drought to appear on water resources and other different life sectors [78]. Recently, the use of multivariate indices has concretely increased [79], especially in cases that take into consideration the multiple characteristics of drought events. The combined index describes droughts in more detail, showing a set of drought characteristics such as intensity, severity, duration, start and end drought period, and areal extent for each region. Applicability of the combination of several drought indices is based on the fact that drought, as a natural phenomenon, results from a combination of different factors such as lack of precipitation, prolonged periods with below average rainfall, and increased evaporative demand as a result of temperature excess [80]. For each constructed drought event, we have identified the essential aspects of drought characteristics. In this paper, seven parameters (severity, intensity, duration, areal extent, peak area, peak month, area involved with peak month) were used to describe the constructed drought events. Details about the definition of these parameters can be found in Table 1. The degree of dryness (DDI) was also assessed as an indicator to evaluate the drought intensity for different drought categories in each area and for each subregion. More details about DDI calculation can be found in Salehnia et al. [81]. Table 1. Parameters used for a drought event evaluation as listed in the database.

Start event, month (ST)
The first month with an index value below a certain threshold End event, month (EN) The last month with an index value below a certain threshold

Drought duration (DD)
Number of consecutive months between drought start and drought end (length of drought event)

Drought severity (DRS)
Sum of the absolute values of drought indices during a drought event Drought intensity (DRI) The mean value of drought indices during a drought event (DRS/DD)

Drought frequency (DRF)
Number of drought events during specific interval, usually expressed as number of events every decade.

Areal Extent (AE)
Percentage of the region with indicator values below a certain threshold

Peak Area (PA)
The largest area exposed to drought during the event (%). Usually the PA is presented with its corresponding date in month and year (AMY)

Peak month and year (PM)
The month with the lowest value of the index over the drought event. This parameter is usually presented with severity value of this month.

Principal Component Analysis (PCA)
PCA is a multivariate technique that is extensively used in meteorological, climatological, hydrological studies, and drought regionalisation as a dimensionality reduction statistical method for extracting structural information in a dataset, and allowing a clearer interpretation of large datasets of dependent variables without losing much of the original information [82]. In this paper, the PCA was used to capture the major spatial patterns of drought across the MED by reducing dimensionality in a group of the DXI-12. Specifically, the S-mode PCA with the varimax orthogonal rotation was applied to the DXI-12 at 312 stations for identification of the regions within the MED that have homogeneous drought features and conditions. The S-mode PCA enables the analysis to describe the areas with close temporal fluctuations [83]. Using the REdaS R package developed in 2015, the Kaiser-Meyer-Olkin (KMO) test [84] was utilized to examine the quality of the principal components of the time series of all drought indices before the application of the PCA method, and thus, to check the suitability of the data for factor analysis. Both the scree plot of eigenvalues and the rule of thumb [85] were used to make the decision on how many loadings to retain for rotation. The PCA was performed using the "prcomp" function in R.

Definition of the Constructed Meteorological Drought Event Database
According to Henriques and Santos [86], a drought event occurs when the value of the variable falls below a fixed threshold. Many thresholds were selected for classifying drought events. McKee et al. [51] established that a drought event starts when the SPI value reaches −1 and ends when it becomes positive. Spinoni et al. [76] indicated that a drought spell requires at least two consecutive months with a negative drought indicator measured. Spinoni et al. [28,87] assumed that a drought event starts when the indicator values fall below a certain value corresponding to a given negative standard deviation (X = −1s) for at least two consecutive months and ends when this value turns positive. In this study, drought durations for each country/subregion were constructed by creating sequences of the number of two or more consecutive months over which the value of indicator persists below a given threshold (X = −1s). The time series of drought durations for each station was constructed by creating a succession of the number of consecutive months, at least two, over which the value of the drought indicator persists below this threshold. The drought event ends when the indicator value turns positive. Each drought event within the constructed database has its own property which reflects the drought dimensions and characteristics. We have investigated the drought events at the national and subregional level to determine the biggest droughts that occurred in each region over the MED during last four decades. Monthly series for DXI-12 during 1975-2019 were ob-tained for each station, and then the corresponding series for each country were constructed. Finally, the subregional series was derived for each subregion. We have also assembled a list of the overall most relevant Mediterranean drought events that took place between 1975 and 2019 by also considering the events that involved more than one region (multi-subregional events).

Developing a New System for Classifying the Drought Events According to Three Different Attributes
Spinoni et al. [28] compared the drought events at different scales according to a special scoring system that was derived from an aggregation of six different parameters. To further identify the most exceptional drought events over the past 40 years, a similar approach will be applied using only the three most important parameters (severity, intensity and the areal extent) to describe and classify the drought events in different areas over the MED. Since the drought is a three dimensional interdisciplinary and multi-scalar phenomenon [73,88], our new approach to classification of the drought events is modified to incorporate these three drought parameters. The new classification is modified to evaluate the drought events depending on the previous grouped parameters. The main objective of this method is to identify coherent time and coverage structures through a three-dimensional array (severity, intensity, and areal extent) of drought period identified by the DXI-12 in order to transfer this qualitative classification into a quantitative one. In other words, a given parameter received grades from 1 to 6, relying on the event itself compared to others in the series, at a regional scale. While this approach has been tested for two different accumulation periods (3 and 12-month timescale), in this paper, we only present the results for 12-month accumulations. The percentile is used as a way to classify the status of each parameter and different percentiles are examined to assess the drought event condition through the combination of the three aforementioned indicators. Table 2 shows the classes of the three selected parameters and their respective grades, and this adopted drought events classification system has established the categories of "extreme drought events" for total grades from 18 to 15, "severe drought events" for total grades from 14 to 11, and "moderate drought events" for grades from 10 to 1 as a new attempt to classify the drought periods according to these attributes. The new classification is applied for each constructed drought event through the DXI-12. Table 2. Drought parameters classes allocated to each drought event and their competent grades.

Grade
DP Status 6 DP > 95th percentile over the recorded events at the subregional and regional scale 5 90th percentile < DP ≤ 95th percentile over the recorded events at the subregional and regional scale 4 70th percentile < DP ≤ 90th percentile over the recorded events at the subregional and regional scale 3 50th percentile < DP ≤ 70th percentile over the recorded events at the subregional and regional scale 2 30th percentile < DP ≤ 50th percentile over the recorded events at the subregional and regional scale 1 DP < 30th percentile over the recorded events at the subregional and regional scale

Atmospheric Circulation-Drought Conditions Analyses
The MED is prone to the effects of the large-scale teleconnection patterns which considerably vary in their strength within the region [89]. Several studies have considered the association between fluctuations in the atmospheric teleconnection patterns (ATPs) and droughts [90][91][92] and their different impacts on the hydrological regime in the MED [91,93]. Here, we focus on the North Atlantic Oscillation (NAO) [94], the Upper-Level Mediterranean Oscillation index (ULMO) [95], the Western Mediterranean Oscillation (WeMO) [96], the East Atlantic/West Russia (EATL/WRUS) [97] and the Scandinavia pattern (SCAND) [98] because a preliminary analysis has showed that they are the best candidates among existing teleconnection indices which describe the drought conditions in the region of interest. Relationships between atmospheric circulation and the percent area of drought based on SPI-12, SPEI-12 and RDI-12 were explored using Pearson's correlation analysis.

Loading Patterns of SPI-12, RDI-12, SPEI-12 and DXI-12 Derived from PCA
The SPI-12, SPEI-12, RDI-12 and DXI-12 were used to identify drought patterns for the 1975-2019 period. For all drought indices, five principal components (PCs) were retained identifying five regions that are well-structured and spatially coherent and show the highest values of loadings and explain the highest variance. The obtained regions are characterised by different drought variability due to different precipitation and temperature regimes in these areas.
High values obtained in the KMO test (>0.71) suggest that all proposed drought indices are adequate for the PCA, which also indicates that we can apply the PCA method to all selected time series. The KMO measures of sample adequacy applied to the drought indices are, respectively, 0.71, 0.95, 0.69, and 0.92 for the SPI-12, SPEI-12, RDI-12 and DXI-12, suggesting that all are adequate for PCA analysis (KMO test > 0.50). Appling the varimax rotation to the PCs has allowed the identification of subregions within the MED that have homogeneous drought conditions. Generally, there is a spatial coherence in the distribution of stations associated with each PC. Some stations have shown the highest correlation with more than one component identified in the PCA analysis. Generally, the main spatial patterns of all drought indices' variability characterise the MED uniformly. The patterns were harmonious and had a coherent spatial structure that was consistent with climate influences over the MED. Although the general drought patterns existing over the MED are captured by all indices, the proportions of variance by the leading PCs are relatively higher in the case of the SPEI-12, RDI-12, and DXI-12, indicating the significant role of PET and the large response of drought to temperature variations. This means that the sensitivity of a purely precipitation-based drought index to drought variability differs from the precipitation and temperature-based indices.
The percentages of the total variance explained by the first retained loadings for the previous drought indices, which explain 10.34%, 15.53%, 17.78%, and 17.58% of total variance, respectively, clearly exhibit consistent spatial structures as well as temporal variation. The second rotated loadings account for 10.08%, 10.5%, 10.49%, and 10.99% of the total variance, respectively, with a mean cumulative variance of almost 59.9% for the five PCs for all indices. Thus, the rotated loadings for SPI-12, SPEI-12, RDI-12 and DXI-12 successfully reflect a coherent structure with five distinct subregions that are characterised by different drought variability levels, which are possibly related to the different rainfall regimes in these areas. The spatial patterns of the first five principal components are almost the same for SPI-12, SPEI-12 and RDI-12 but, in the case of the SPEI-12 and RDI-12, the loadings are higher and correspond to a wider extent than SPI. Therefore, our study is based on the analysis of the drought-event records, which are constructed according to DXI-12 for five subregions spread over 20 countries along the Mediterranean area ( Figure  2 and Table 3).

Meteorological Drought Events Identified by DXI-12 and Their Characteristics at Country and Subregional Scales
Based on the analysis of the DXI-12, a new dataset of meteorological drought events was constructed during the period between 1975 and 2019. The overall number of drought events is 375, 169 of which occurred after 1996 at the country level. Any drought event that involved more than one country was allocated to represent its subregion. At the subregional scale, 19 drought events are detected in IBE-NWA, whereas the dataset has reported 20, 13, 13 and 15 drought events in EMME, AEG, BLK-FR and NA, respectively. The spatiotemporal characteristics of each drought event, including its duration, severity, intensity, areal extent and area involved with peak were examined at the country and subregional scales. Considering that the area affected by drought is a very important drought characteristic to evaluate the drought risk, we first reported the widest drought events in each region that hit more than one country with a peak regional area that exceeded 50% ( Table 4). The corresponding Degree of Dryness (DDI) values for the peak year for each drought event at the subregional scale are also presented.
The areal extent of annual droughts (Table 4) shows many wide droughts over the MED that have struck more than 40% of the area. The exceptional multi-year drought event (1999-2001) might be marked as critical in all subregions, as the MED was exposed to drought episodes which affected the wider area by more than 50% in the IBW-NWA, EMME and AEG, while 37% of the area of NA was affected by this exception event. In terms of individual drought events and their other main characteristics, the differences were also significant between subregions. The meteorological drought events in the MED usually exhibit spiky behavior, inherently related to the high temporal variability of precipitation in this region. Generally, the most pronounced drought events in the last two decades experienced larger areas than other periods, compared with the late 1970s and 1980s. Sub regional drought area reached its highest level in several years after 2015 with severe and large coverage drought episodes during 2015-2018. Considering both the average of involved area and DDI values, the 1999-2001 period is the period of the broadest drought spell in the MED as a whole and the driest period in the EMME, with the drought also extended until 2002 in the NA, at which point it had the largest severity and highest intensity ( Table 4). These findings reinforce what is already established by Hoell et al. [99] and Barlow et al. [100] who considered the catastrophic 1999-2001 droughts to be one of the worst drought periods in the EMME and the period with the widest drought impacts in the Southeastern Mediterranean. Furthermore, it appears that drought distribution within a multi-year is one of the most important factors in determining the severity of drought episodes. In Europe as a whole, a very spatially extensive event was shown between 2017 and 2018 that covered at least 35% of the total area and peaked in December 2017 with 74.3% in the Iberian subregion. The highest severities were also recorded during the 2017-2018 drought incidents in the IBE-NWA, as shown in Table 4. These results are in agreement with other studies [101] that ranked the 2017 drought period as one of the most severe drought events in Europe since 1979. In the IBE-NWA, an extensive coverage drought period (52.7%) was also notable in the mid-2000s, typified by high severe conditions (DDI = 11.9), and peaked at 74.2% in August 2005. The same degree of dryness was identical for the AEG in 1990, with a peak area of 77.6%. Severe drought events with high percentage areal extent were experienced in the AGE through the 1990s (59.2% in 1990), 2000s (65.3 in 2001) and mid 2010s (50% in 2014), with an almost complete lack of wide drought events in 1980s. Table 4. List of the widest drought events detected by the DXI-12 and their corresponding areal extent, the area involved by the peak month, and the Degree of Dryness (DII) of widest area at the subregional scale during 1975-2019. A wide drought event was recorded in the BLK-FR in the early 1990s, covering more than 50% of this area and followed by multiple periods of large drought extent in the mid-2000s with more 45% of spatial coverage that were not repeated until the late-2010s. Moreover, large sectors of the Adriatic and west central Mediterranean (BLK-FR) suffered from a persistent spring and summer drought in 2003 affecting 31% of this area and reaching its largest coverage in December 2003 (50%). This intensive drought episode was accompanied with an exceptional heatwave that broke the temperature record, particularly cross southern and central France [102,103], indicating a very good example of a compound drought and heat event. In the NA, drought events have been recorded in different periods and there was a succession of wide drought events since the 1980s but the most severe and widest drought event was recorded in the late 1990s, a period characterised by extremely severe long-lasting droughts. The latter half of the twentieth century also emerged as the driest of the last nine hundred years in the NA [104]. Furthermore, the multi-year drought (1999)(2000)(2001)(2002) in the NA was reported as the most severe of the last five centuries [105].

Region
Additionally, it is apparent that the late 1980s drought episode in the IBW-NWA, AGE, and NA droughts are not such extreme events compared with other drought episodes. The results suggest they display considerably lower severity and intensity with respect to the multi-year 1999-2001, 2005-2007 and 2015-2018 droughts. Table 5 and Figure 3 list the top five drought events that struck each subregion ranked by duration, severity, intensity, areal extent and degree of dryness. Several drought hotspots were indicated during last two decades along the MED. The most severe and intense droughts of the 1990s and late 2010s are clearly visible for the region as a whole. The 1990s and 2000s also stand out in terms of spatial area in drought. The most spatially extensive events were two of the five long-lasting droughts (2008)(2009)) in the southeastern Mediterranean across the EMME, and in 2002 in NA. The multi-year drought in 2015-2018 displayed a highly unusual spatial pattern that affected huge parts of Europe including the AEG and BLK-FR, as shown in Table 5. Dry conditions were reported in over 90% of central-western Europe, which hit record-breaking values in more than 25% of the area [100]. Furthermore, the spatial coverage of the top drought events in all subregions has displayed large affected areas with the highest coverage in the EMME (≈76%) in the late 2000s. For the same period, this area also displayed its most long-lasting drought event, lasting 47 months, which was the longest drought episode in the last 900 years in the EMME [38]. This prolonged drought episode (2007-2012) was one of the worst long-term droughts in modern history, and was associated with the steepest decline in crop productivity in thousands of years [100]. With this continuous harsh drought phenomenon, hundreds of thousands of people in have been displaced in the Euphrates and Tigris catchment basins, leading to an increased risk of regional conflicts, and thus contributing to the outbreak of war and helping to spark the Syrian Civil War [106,107]. The most intense drought events were recorded in the late 1990s and mid-2000s in the IBE-NWA, EMME, and NA, whereas the late 2010s experienced the most intense events in the AEG. Results have revealed top drought clusters in specific decades such as the mid-2000s and late 2010s, and more widespread events were most notable during the 2000s, with the highest degrees of dryness also shown after 1990 in the region as a whole. The highest degree of dryness was recorded in the 1990s in the EMME, NA and AEG, while the other subregions recorded their highest scores in the mid-2000s. These results are in line with other previous ones which confirmed that there is clear evidence of a strong drying trend in the EMME starting in the early 1960s, with the lowest precipitation amounts associated with the driest years in the late 1990s [99,108]. Very similar findings have been obtained by computing the DDI to describe the area affected by the droughts in the AEG. This result was also reported by Korner et al., 2005 [109], who indicated that the most crop yield declined in the Aegean area, namely in Greece, during the 1990s due to the consecutive severe droughts that hit this area. Varol and Ertugrul, 2015 [110] have also reported an increase of 14.75 times in the burned forest area in the Aegean region after 1990 due to increased droughts in this area. Finally, the results suggest the MED as a region which is characterised by a marked and high spatial variability of drought events, even for the most extreme and exceptional episodes recorded in the region, which is consistent with the results of previous observation [111]. The higher values of the DDIs during the last two decades also represent the trends toward drier conditions in all subregions ( Figure 4).  The results demonstrate that the MED has exhibited a high drought tendency during the past four decades with more frequent and severe droughts occurring after the 1990s and in the 2000s, and fewer droughts in the 1980s. Based on Figure 4, the DDIs change in the MED exhibits a large subregional increase over the past four decades, with a remarkable upward (drying) shift since the late 1990s. Generally, the MED displays the highest tendency toward extreme dry conditions, represented by higher DDIs values (>10) in the IBE-NWA, EMME and AEG, where the highest DDI was recorded in 1999 in the EMME (14.3). Large dry conditions were also present for the other four subregions in the 2000s, while relatively small dry conditions were present from the early-1980s into the mid-1990s.

Three-Dimensional Identification and Characterisation of Meteorological Drought Events at Subregional and Regional Scales
According to the new classification system that we previously described, multiple drought events have been recorded between 1975 and 2019 at the national and regional scales. Following the new classification, the meteorological drought events were firstly investigated at the country scale (Table 6), as the impacts of drought, corresponding damages and other economic losses are generally recorded at the national scale for each drought event [112]. Then, all events at the country scale which involved more than one country have been grouped to create a new list of drought events at the subregional scale (Table 7). Only drought events with high grade (>13) during 1975-2019 have been considered for both the national and regional scales. In order to prove the effectiveness of our approach in extracting drought periods, we have listed some scientific literature and scientific text-based reports that discussed single or several drought case studies across the MED. Generally, the literature is not as extensive about all drought events in the MED, particularly in southern and southeastern areas. However, all the drought events reported in Tables 6 and 7 has been previously documented. The spatial patterns of the obtained drought events in this study are in good agreement with most of the major drought events presented in the international database of text-based reports and add a series of other drought events, particularly in the Eastern Mediterranean and Maghreb countries. As shown in Table 6, the largest numbers of droughts occurred from the 1990s up to the present and are considered to be exceptional events. All countries have experienced extreme drought events becoming more severe and more intense in the last two decades. The areal extent of annual droughts (Table 6) shows that many wide droughts, in all Mediterranean countries, have struck more than 50% of the area. The period after the late 1990s might be marked as critical, as the MED was exposed to drought episodes which affected more than 90% of countries in the wider area. For example, the event in 1989-1991 was the most severe one in Greece, Italy, and Croatia over the last four decades and the area under these very dry conditions has dramatically expanded to reach its maximum (97%) in Croatia in this exceptional drought spell. The late 1990s were characterised by very exceptional events that hit the Southern Levant countries (Syria, Lebanon, Israel, Jordan, Palestine, and Turkey), two of them ending in the early 2000s (Syria and Turkey). The 1999 widespread drought also affected neighboring Jordan and Israel, which were hit by their worst drought over the past 100 years in 1999 [113]. The highest drought severity ever recorded at the country level was detected in Algeria in the late 1990s, namely between 1999 and 2002 ( Table 6). The recent drought (1999)(2000)(2001)(2002) in the Mediterranean parts of Algeria and Tunisia was classified as the longest period of consecutive drought years since at least the middle of the 15th century and the worst drought event that struck these two countries since at least the middle of the 15th century [105,114]. The period between 2000 and 2001 included very severe and prolonged droughts in France, Italy, and Greece, which started in 1999 in Greece. This period is seen as the most serious event in Italy since the 1920s [115]. The extensive and persistent droughts in the mid-2000s significantly exceeded the other events in their spatial extent to include all the Mediterranean countries. Table 6 shows many wide droughts that struck more than 40% of the area, such as the event which spanned between 2003 and 2004 in Croatia and affected 74% of the entire country. The results display that around 34% of the area that was hit by drought in 2006 and substantial parts of Morocco and Spain (50 % and 70%, respectively) were also drought-prone between 2004 and 2005. Very dry conditions occupied several western and eastern countries in the mid and late-2010s as one of the most recent exceptional drought periods was consistent with a major drought in Spain, Italy, France, Algeria, and Cyprus and severe drought from 2016 to 2018 in Turkey. Based on these results, the two decades of 2000-2010 can be classified as a period with the most severe and long-lasting droughts during the last 50 years that appear in all Mediterranean countries [87,116]. More specifically, the droughts of the late 1990s, mid-2000s, and the multi-year drought series around 2015, appear to be the most serious and large spatial extent drought episodes, which covered almost all Mediterranean countries and reaching high severity levels. It should be noted, however, that recent droughts at particular locations (e.g., Algeria and Tunisia) are more severe than in other Mediterranean countries, as Algeria and Tunisia recorded the highest drought severity in 2002 and 1981, respectively. Although the other countries have a broader spatial extent, the NA experienced droughts in the late 20th century that were exceptional in the context of the prior 500 years [117]. It should be noted, however, that recent droughts at particular locations (e.g., Algeria and Tunisia) are more severe than other sites, although the others have broader spatial extent as the NA experienced droughts in the late 20th century that were exceptional in the context of the prior 500 years [117]. Compared to previous works on global-scale drought change [28,66], we provide, through this paper, a more comprehensive and detailed drought-event dataset for the MED by considering multiple meteorological drought indices to identify drought in this region. Table 7 reports a corresponding compiled list of the extraordinary, extreme and severe drought events, with scores exceeding 13 at the subregional scale between 1975 and 2019 identified by the three-dimensional array system, and presents twenty-four multiregion severe events detected by the DXI-12. Moreover, Table 7 is based on input data up to December 2019, but a few drought events were still ongoing at the end of 2019. These updates will be included in the coming dataset in future research. A notable uneven temporal distribution of the drought events is detected as most events occurred after the 2000s, which have been characterised by increasing and more severe and extreme rainfall events mixed with prolonged dry conditions [40,143]. In the 2000s, 17% of the global population (1.1 billion people), who mostly live in the Middle East, North Africa, and South and East Asia, experienced critical water deficiency and high water stress [33] with 61-89% of the population under water shortage. Severe, intense, wide, and long-lasting drought events are frequent in all subregions across the MED over the last decades and were more pronounced in the Middle East and some parts of the Eastern Mediterranean. The largest drought severity was recorded in the EMME (>53) with very high intensity (>2.1) between 2014 and 2019. The multi-drought event in the mid-2000s also experienced high severity (>40) with very high intensity (2.2). The common pattern of drought occurrence displays marked spatial variability of the Mediterranean drought events, as droughts in this region seldom affect the area as a whole; the matter could be explained through the east-west climate dipole as the dominant component of variability between western and eastern sites of the MED [144]. Moreover, the EMME reflects transition between summer subtropical in the north and semi-arid climates in the south, making it more vulnerable to climatic sensitivity [141]. Combining the results reported in Tables 6 and 7 and Figure 4 makes it possible to clearly identify the "hotspots" hit by more severe, more intense and wider meteorological drought events at the national and regional scales during 1975-2019. The recent drought in the Levant region between 2008 and 2012 was probably the worst drought of the past nine centuries. It was 50% drier than the worst drought period of the previous 500 years and 20% drier than the driest spell over the past 900 years [38]. A very severe drought period struck the northern parts of Africa between 2004 and 2005, leading to a reduction in the economic growth rate by 2.2% [142]. In 2005, severe droughts were recorded in the west of the Iberian Peninsula, southeast France, which is consistent with the results supported by Quintana-Seguí et al. [141].
In 2017, one of the worst drought periods in 150 years hit large parts of Europe and threatened the Balkans, eastern and southern Spain, and most of the Italian peninsula. Table 7 shows that exceptional droughts occurred in the 2010s and encompassed five subregions with grades higher than 13, while the AEG experienced the most exceptional drought events in the early 1990s and the mid-2000s with grade (16).  (Table 8). Table 8 illustrates the number of drought events at the regional scale and their average duration, severity, intensity and their spatial extent identified by the DXI-12. The average degree of dryness was also reported in Table 8. In order to evaluate the spatial and temporal patterns of drought events, we have compared the average severity and intensity ( Figure 5), frequency and spatial extent through the two subperiods. The drought periods identified by the DXI-12 were longer, wider, most severe, most intense and most frequent over 1997-2019 (Table 8 and Figure 5). The DDI values have also exerted visible upward tendency in all subregions. The Table 8 shows the frequency of the drought events in all subregions has increased from 1975-1996 to 1997-2019. The same applies to other parameters, as the drought intensity, severity and spatial coverage experienced a remarkable change in all subregions. Additionally, the DXI-12 identified dryer conditions and an increase in the area affected by droughts over all subregions. A more remarkable change is shown across the EMME, which indicates the highest drought exposure in the second period. The arid and semi-arid areas over the EMME are generally characterised by high dry conditions induced by increases in temperature and evaporation [38,107]. The highest increase in the different drought parameters was demonstrated in this subregion, while the change towards largest number of drought events mainly occurred in the IBE-NWA. Furthermore, the dry climate of EMME has not only tended to experience more intense drought events over the last two decades, as shown in this research, but also more prolonged and long-lasting events. This was confirmed by Tabari and Willems (2018) [145], who found evidence that the number of dry days is increasing and the drought episodes will be longer (up to 90%) over 80% of the Middle East. The most relevant decrease in severity was clearly evident in some parts of BLK-FR area and NA, mainly in Croatia, southeast France and Spain, and Tunisia. The results have shown a clear correspondence between the severity and intensity of drought events, particularly in the Levant area and Turkey. Based on the DXI-12, a significant change towards larger frequencies of moderate, severe and extreme drought events was evident, in all subregions, in 75%, 81% and 82% of the total stations, respectively ( Figure 6).

Evaluation of Different Drought Events Parameters during Two Subperiods
The increases in the number of droughts over the MED have been detected between the subperiods 1975-1996 and 1997-2019 in all subregions at different categories. During the second period (1997-2019), the increasing trend in drought frequency based on DXI-12 was experienced over the majority of the MED. The results revealed a clear signal concerning the increase in the average number of drought events over the MED by 1.36 events per decade (not shown). The EMME demonstrates the highest increase in drought occurrence by more than one event per decade (+1.9 event/per decade), where Syria and Israel show the highest increase in drought frequency (more than two events per decade). On a spatial basis and based on the DXI-12, the frequency of moderate droughts shows similar patterns for all subregions, but the moderate drought events display a larger increase in the last two decades than both severe and extreme ones. The EMME retained the largest share of the increase in all extreme, severe and moderate events, with increases of almost one moderate event, 0.7 severe events, and 0.3 extreme events per decade ( Figure 6). The AEG, for example, experienced a pronounced increase in severe drought events by 0.  [28,76] following different approaches ( Figure 6). Thus, drought occurrence has shown an increasing trend in addition to more prolonged and frequent periods, in last two decades, at the regional level ( Table 8). The change in drought event frequency was also assessed according to the SPI, SPEI-12 and RDI-12 between the subperiods 1975-1996 and 1997-2019 (Figure 7). An overall increase in the frequency of dry spells appears over the MED for all the combinations of SPI-12, SPEI-12 and RDI-12 during 1997-2019, with clear spatial patterns, suggesting that droughts were more frequent in the last two decades. Based on the SPI-12, SPEI-12 and RDI-12, the results indicated that 70%, 96% and 84% of the total stations have reported an increase in the frequency of drought periods during 1996-2019, respectively (Figure 7). Although there is a coherent change in the characteristics of the original drought events, which was identified by all indicators, there is a general increase in the frequency of the drought events that were studied using SPEI-12 and RDI-12 compared with those obtained by the SPI-12. Therefore, their capacity to detect drought conditions is potentially higher than the SPI-12 due to their incorporating of the influence of temperature on multitemporal drought. The representation of the full water balance by SPEI and RDI makes them superior to SPI in terms of the accurate assessment of meteorological droughts with increasing climate [146,147,148]. Most parts of the Iberian Peninsula have also displayed consistent characteristics and experienced a drought frequency increase identified by the different drought indicators, and these were more evident when measured by the SPEI-12 and RDI-12. The results further suggest that drought events could exhibit entirely different characteristics when we include the influence of the temperature in the drought indices.

Change in Precipitation (PP) and Potential Evapotranspiration (PET) between 1975-1996 and 1997-2019
Changes in precipitation amounts and PET values between 1975-1996 and 1997-2019 were analysed to explore whether the dry and wet tendency reflected by the SPI-12, SPEI-12 and RDI-12 was driven by the precipitation or temperature anomalies-that were included in PET-, or both (Figure 8a,b). It is not easy to understand the actual roles of both temperature and precipitation variations in drought changes over the MED during the last two decades. This is due to nonlinear interactions among temperature, precipitation, and droughts. Although an extended period of reduced precipitation plays an important role in producing drought conditions over the MED, the effect of increased temperature is notable as a significant factor in enhancing drought by enhancing the evapotranspiration rate. According to SPI-12, 62% of total stations experienced a drying tendency. In the last two decades, some areas across west Turkey, south Spain, and the Adriatic experienced a wetting tendency. The more extreme negative values for both the SPEI-12 and RDI-12 compared with the SPI-12 (not shown), particularly in the semi-arid areas in the southeastern Mediterranean, indicate the impacts of temperature that may initiate or increase drought via the enhancement of PET. Although the RDI-12 identified dryer conditions in the last two decades compared with the SPI-12, the SPEI-12 displays the highest contribution of PET to drought conditions compared with the SPI-12 and RDI-12, particularly in the Balkans, north-east Tunisia, and west Turkey. It was concluded that the DXI-12 was able to provide a physically sound, temporally flexible index, which can be directly combined with all possible variants and linked to the climate conditions, that allowed us to list hundreds drought periods and characterise several exceptional drought events between 1975 and 2019. In this sense, DXI-12 can evaluate the actual drought conditions and effectively lessen overestimation of any drought conditions by SPEI-12 in terms of drought duration and intensity. All drought events detected by the DXI-12 were previously confirmed and documented at both the country and subregional scales, while some drought accidents, especially the moderate and short ones, could not be detected by SPI-12 and RDI-12 as well. According to the SPEI-12, the dry conditions across the MED have been enhanced in recent years as a result of increases in temperature and decreases in precipitation in approximately 63% of total stations. The findings also show that approximately 28% of the total stations experienced drier conditions which resulted only from positive temperature anomalies over the last two decades, where the SPI-12 showed an increasing trend, while the rest of the indicators showed a significant decreasing trend (Figure 8). There is a detectable and significant decreasing trend of precipitation and increasing temperature over the Mediterranean and Middle Eastern countries [33,100,107]. Thus, this result indicates that the positive trend of precipitation is not strong enough to outweigh the temperature increase in this area, and that the drought exhibits a larger response to temperature variations, which is more than its response to the precipitation anomalies in these areas. This proportion decreases to approximately 15% for RDI-12. Generally, the MED as a whole displays a clear example of the worst drying combination (decrease in precipitation amount associated with increase in the PET).

Drought Dynamics, Large Scale Atmospheric Controls of Widespread Drought Events over the MED
The correlation of the annual percent of drought area extent based on the SPI-12, SPEI-12, RDI-12 and DXI-12 indices with large-scale atmospheric fields over the whole study period (1975-2019) is investigated. The correlation analysis gives a more detailed evaluation of the processes linking large-scale atmospheric variation to Mediterranean meteorological droughts. In addition to the general relationship between the percentage of MED experiencing meteorological droughts and ATPs (Table 9), the widespread droughts and their atmospheric drivers were also investigated. The spatial pattern of correlation is virtually identical for the SPEI-12, RDI-12 and RDI-12, with positive correlations with NAO and ULMO and negative ones with WeMO, EATL/WRUS and SCAND (Table  9). Table 9. Correlation coefficients between the percent areas in drought identified by SPI-12, SPEI-12 and RDI-12 and ATPs (only coefficients significant at 95% confidence level are shown).

RDI-12
Despite the similarity of the spatial patterns of correlation coefficients presented in Table 9, the correlation of the area in drought with the NAO is the weakest. Generally, the areas in drought obtained by the SPEI-12 and RDI-12 are better related to the ATPs in comparison to SPI-12 and detected the highest correlations with DXI-12 (Table 9). Statistically, the WeMO and EATL/WRUS are remarkably correlated with the area in drought identified by SPEI-12, RDI-12 and DXI-12. The widespread dry conditions over central and Western Europe (BLK-FR) were associated with negative phases of WeMO and SCAND and a positive phase of ULMO, as shown in Table 9. The results presented here show a major role played by the WeMO in drought within the EMME and BLK-FR. The negative phase of the WeMO generally shows very wide spatial coverage of drought over the EMME and BLK-FR. The large area under drought noted over the EMME during 2011-2012 (47%) was often linked with severe negative values of the WeMO. The severe negative phase of the WeMO (<−1.5) over 2017 was also connected with the largest percentage of the affected area in the EMME. Similarly, the ongoing negative annual phase of the WeMO between 1985 and 1991 was associated with very large area under drought in the BLK-FR. This affected area spanned over large territories in France and Croatia which were hit by severe drought episodes and involved 61.2% of the entire subregion. The consecutive low annual negative values of the WeMO between 2002 and 2018 were found to be associated with very high percentages of the area under drought in the BLK-FR. With the exception of the BLK-FR, the negative phase of EATL/WRUS was often associated with large drought spatial coverage within the MED. The persistent negative anomalies in the EATL/WRUS values in the late 2010s (July 2017-September 2018) were associated with high spatial coverage (40%) in the IBE-NWA. Furthermore, the exceptional multi-year 2007-2012 drought in the EMME was developed in a predominantly negative phase of both the WeMO and EATL/WRUS patterns. The combination of the negative phase of WeMO and EATL/WRUS over 2008 resulted in a strong influence on the drought areal extent in the EMME (76%), as shown in Table 5.
The major drought episodes that occurred between 1999 and 2002 in the NA were driven by qualitatively similar conditions with the complete domination of the negative EATL/WRUS, which had the largest mean drought area in 2002. The negative phase of the EATL/WRUS is usually linked to below-average precipitation and above-average temperatures in these areas [148,149], which clearly explains the strong negative relationship between the EATL/WRUS and the area in drought except for the BLK-FR.
The impact of the ULMO significantly appeared through both the 1994 and 2012 exceptional droughts. The very large areas of the IBE-NWA affected by droughts during 1992, 1999 and 2012 (52%, 44% and 62%, respectively) were concordant with the positive phase of ULMO. The widest affected area in in the BLK-FR during 2002 (Table 5) was remarkably associated with a large predominance of negative values of the WeMO, SCAND and positive values of the ULMO. These findings perhaps suggest that the WeMO and EATL/WRUS patterns show most potential role in explaining the spatial evolution and characteristics of large-scale drought over the MED with lesser influence of the NAO ULMO, and SCAND. In NA, the EATL/WRUS may determine the climatic water balance drought, while the NAO and ULMO play a significant role in regulating the internal rainfall variability, and thereby the spread of drought events. These results correspond with previous findings [150,151], which indicated high drought conditions during the positive phase on the NAO in southern Europe and North Africa. At the regional scale, similar correlation patterns have been detected with drought severity, intensity and area in drought. Correlation patterns of all selected ATPs indicate relationships of Mediterranean drought severity, intensity with atmospheric circulation at the regional level (Figure 9). High drought intensities and severities are significantly associated with the negative phases of both WeMO and EATL/WRUS across the MED. The results confirm the higher capability of the SPEI-12 and RDI-12, and DXI-12 compared with SPI-12 in improving climatic water deficit detection in the MED, which confirms the capacity of temperature in detecting drought conditions. The drought intensities and severities obtained by the DXI-12 showed the highest negative correlation coefficient values with EATL/WRUS and WeMO, respectively.

Conclusions
This study has provided a comprehensive dataset of exceptional meteorological drought events over the Mediterranean during the 1975-2019 period that contains hundreds of constructed drought events at the national and subregional scales. Three drought indices were computed at 12-month accumulation scales and a new combined indicator was obtained (DXI-12) as the average of these three indices to identify the meteorological drought due to both rainfall shortage and climatic water balance. The drought event databases are identified by monthly series of the DXI-12 at the national and subregional scales. Based on the drought-event dataset, we have also investigated in detail the main characteristics of the obtained drought events in different subregions across the MED, including duration, severity, intensity, spatial coverage, peak month, area involved at peak month and the frequency of drought occurrence. This new approach has proved to be effective in detecting the most relevant Mediterranean drought events at the national and subregional scales as it is based on data obtained with high-quality control, a multi-indicator method, and a strong statistical setting. The evaluated drought quantification, which includes these three indices, gives a more reasonable estimation of water demand as a result of including evapotranspiration. Generally, our newly developed DXI-12 provides better information about drought identification and more realistic characteristics of dry spells compared to a single drought index and, thereby, provides a better tool for monitoring water resources, especially in semi-arid and arid ecosystems. According to the results presented in this study, the MED is shown to be a region that is water-stressed and prone to severe and exceptional droughts. The drought episodes show that the EMME and most areas of the NA and southern Europe, including south Italy and Spain, underwent a drying phase during the study period, with the highest drought severity and intensity from the late 1990s onwards. These findings are in line with recent scientific literature [28,101]. The combination of individual indicators depicts the spatial extent of a drought situation across the MED and gives an overview of the drought drivers. The DXI-12 indicates that droughts over the MED were mainly driven by temperature increases with the inherent increase in the PET, particularly in the humid area, whereas rainfall did not change appreciably. A new grade classification approach for identifying the most relevant extreme events is identified by the three-dimensional array (severity, intensity and spatial coverage). The new classification approach and database extracted from combined indices the time series have listed twenty-four multi-region big drought events, which were previously documented. The new classification has allowed us to discover the most hotspots that have displayed the highest severity and intensity and included the largest area over the MED. The significant trends towards larger drought areas and higher drought severities and intensities in all subregions and mainly occurred in the EMME. Although the actual precipitation anomalies exhibit relatively large effects on drought occurrence, we can also conclude that the rapid warming of the MED in recent decades presents a relatively greater role in determining droughts than that of precipitation variations in several areas, especially in the Balkans and west Turkey.
In this paper, we have detected several drought hotspots that have tended to experience more severe, intense wide and more long-lasting drought events between the subperiods 1975-1996 and 1997-2019 in all subregions. Dynamically, the region lies at a crossroads of large-scale influence from the Atlantic Ocean and the Mediterranean Sea, which results in regionally severe and prolonged droughts. Dynamically, the region is a junction of regional influence from the Mediterranean Sea, Europe, and Asia, as well as large-scale influence from the Atlantic Ocean and the Mediterranean Sea. A complex range of atmospheric circulation patterns influences have been shown to have at least some influence on the droughts of the MED, both nationally and regionwide, particularly the WeMO and the EATL/WRUS.
The WeMO pattern and the EATL/WRUS appear to play the strongest role in forcing wide droughts at the subregional scale, including the two widest and most devastating droughts of the last 40 years (1999-2001 and 2007-2012). A relative role is played by the NAO and ULMO in the NA in these episodes; in addition, a relationship between the SCAND and the area in drought in western and the Adriatic is found. The drought event datasets produced by this study still need to be improved in some respects. The extension of the data to drought events to include more historical data is significant for studies of climatic dynamics and ATPs in order to understand the long-term variation of droughts and develop monitoring and forecasting tools to help to adapt to and mitigate drought. Moreover, including realistic socio-economic data will strongly improve the prospects for drought risk assessment in the future as many combined meteorological, hydrological, and social drivers contribute to drought occurrence. Finally, these constructed exceptional drought events at the national and subregional scales can help researchers and relevant stakeholders in dealing with the multitude of such drought events and assess their risk by investigating the correlation between the severity and duration of each event and documented impacts in different sectors. As an initial possible application, we are planning to exploit this dataset in exploring the synergy and combination between these exceptional droughts and heatwaves across the MED.

Conflicts of Interest:
The authors declare no conflict of interest.