Multi-Temporal Analysis and Trends of the Drought Based on MODIS Data in Agricultural Areas, Romania

: The aim of this study is to analyze the performance of the Drought Severity Index (DSI) in Romania and its validation based on other data sources (meteorological data, soil moisture content (SMC), agricultural production). Also, it is to assess the drought based on a multi-temporal analysis and trends of the DSI obtained from Terra MODIS satellite images. DSI is a standardized product based on evapotranspiration (ET) and the Normalized Di ﬀ erence Vegetation Index (NDVI), highlighting the di ﬀ erences over a certain period of time compared to the average. The study areas are located in Romania: three important agricultural lands (Oltenia Plain, Baragan Plain and Banat Plain), which have di ﬀ erent environmental characteristics. MODIS products have been used over a period of 19 years (2001–2019) during the vegetation season of the agricultural crops (April–September). The results point out that those agricultural areas from the Baragan Plain and Oltenia Plain were more a ﬀ ected by drought than those from Banat Plain, especially in the years 2002, 2007 and 2012. Also, the drought intensity and the agricultural surfaces a ﬀ ected by drought decreased in the ﬁrst part of the vegetation season (March–May) and increased in the last part (August–September) in all three study areas analyzed. All these results are conﬁrmed by those of the Standardized Precipitation Evapotranspiration Index (SPEI) and Soil Moisture Anomaly (SMA) indices.


Introduction
Drought is a very complex phenomenon without a generally accepted definition, but most of them highlight the deficit of precipitation in a certain period of time. Because of the complexity of the natural hazard that affects many aspects of society and environment [1,2], drought is commonly classified into four types: meteorological, agricultural, hydrologic, and socio-economic [3,4]. According to the World Meteorological Organization (WMO), meteorological drought depends on the atmospheric conditions and is defined as: "prolonged absence or marked deficiency of precipitation" [5]. Meteorological drought is analyzed in terms of the intensity and duration of the dry period observed based on meteorological data [4]. The intensity of meteorological drought can lead to agricultural drought, associated with soil moisture deficits, higher soil temperatures and increasing evapotranspiration. Agricultural drought happens when the soil moisture available to plants has dropped to such a low level that the crop yield is reduced or even totally compromised. The hydrological drought is analyzed and eastern parts of the country [27]. According to Carrão et al. [28] and Vogt et al. [2] those regions records a higher drought risk for agricultural production compared to the rest of the country and even to other areas of Europe. The average yields of various crops represented only 35-60% of the potential yields, during the extremely droughty years [27].
In the last two decades, drought analysis based on in-situ measurements reveals the most severe droughts were recorded: i. in the Southeastern part of Romania (Baragan Plain) in the years 2000,2003,2007 and 2012 [29]; ii. in the Southwest (Oltenia Plain) in the years 2000,2002,2007 and 2012 [30]; iii. in the Western part of Romania (Banat Plain) in the years 2000, 2003 and 2007 and less intense [31,32]. Significant agricultural areas were affected in 2000 and 2003 in Banat Plain [33], and in Eastern, South-Eastern and Southern parts of the country in 2007 [32,34], according to the analysis of drought indices derived from satellite images. Drought not only affects agricultural crops, but also causes vegetation fires [35].
However, indices obtained from satellite images are used not very often in Romania to compensate the lack of meteorological data in areas where there are no measuring stations, as well as the lack of information on the affected areas and the intensity of the phenomenon. In general, the coverage of these areas is done by interpolating existing data. So far, remote sensing indices focused on monitoring and assessing drought are not well known, with the exception of few common indices such as Normalized Difference Drought Index (NDDI) [34,35] and also vegetation indices such as NDVI [30,33] and Leaf Area Index (LAI) [34], which are used to identify the impact of drought on vegetation and also certain climatic parameters such as Land Surface Temperature (LST) used for the analysis of high temperatures [33,35].
The Drought Severity Index (DSI) developed by researchers from NTSG University of Montana [36] is based on the Evapotranspiration (ET)-Potential Evapotranspiration (PET) ratio to compensate for the lack of direct meteorological data input, ET being directly related to for water, carbon and energy cycles of the land surface [37]. The evapotranspiration is an important indicator for the identification and monitoring of drought. Also, for assessing the vegetation response to drought, NDVI is used as a very good indicator of vegetation drought or water stress in plants [38]. This approach is set to provide high quality temporal and spatial information over vegetated surfaces but also to get rid of the uncertainties of meteorological data and modelling [37].
DSI is a standardized index, where the values theoretically range from unlimited negative values (drier than normal) to unlimited positive values (wetter than normal) [39]. The negative values of the DSI show how much the vegetation has been affected by meteorological drought, in other words-it highlights the agricultural drought.
The main purpose of this paper is to analyze the performance of the DSI in Romania and its validation based on other data sources (meteorological data, soil moisture content (SMC), agricultural production). The DSI is used for the first time over Romania, while other remote sensing indices (NDVI, NDWI, TCI, WHI etc.) have been used in limited studies [40,41]. Once to validate the performance of this index, there can be increase the confidence in using it.

Study Areas
The study was carried out for the agricultural land of Romania, which is about 13.3 million ha [42,43], with the focus set on three main agricultural areas: Baragan Plain, Oltenia Plain and Banat Plain, which represents about 17% of the total. Their agricultural area is around 685,000 ha in Oltenia Plain, 664,000 ha in Banat Plain and 926,000 ha in Baragan Plain, according to Corine Land Cover Data (2018) [43].
Located in areas with different climatic conditions, due to the natural obstacle of the Carpathian Arch but also to the Eastern atmospheric circulation, rainfall does not fall evenly in all three field areas.
The Banat Plain is located in the Western part of Romania, Oltenia Plain-In the Southwest of Romania and the Baragan Plain-in the Southeastern part of the country, Figure 1. The climate of Oltenia Plain is temperate-continental with Sub-Mediterranean influences characterized by mild winters and a significant precipitation regime (especially in autumn) [44]. The average annual temperatures were 11.5 • C, about one degree higher in the South and South-West than in the North and North-East of the plain, according to the climatological period 1981-2010, Table S1. Also, the precipitation amount recorded 532 mm, increasing from South to North by about 100 mm. The correlation of soil type (sandy) with the deforestation of forest protection curtains, plots fragmentation (land restitution measures) and the irrigation system degradation lead to the occurrence of deflation [45]. As a consequence, those changes lead to the emergence of drought and finally to decreasing yields [30].
The climate of the Baragan Plain is also temperate-continental, but the continental character is more enhanced than in the other two regions [44]. It is a warm region (11.2 • C) with rich thermal and radiative resources and low precipitation (480 mm/year), Table S1. It is a dry region with a significant water deficit during the active crop vegetation period (April-October) [44,46]. The expansion of dry areas has been facilitated by the same factors as in Oltenia Plain (the destruction and abandonment of irrigation systems, deforestation, excessive land fragmentation and inadequate farming practices) [47]. However, according to the values of the aridity index, climatic conditions are more severe in this region than in the other regions, which enhance the risk of agricultural damage [44,48].
In the Banat Plain, the climate is temperate-continental with influences of Atlantic Cyclones, which bring significant precipitation and moderate temperatures [44], resulting in higher annual precipitation (601 mm) and a mean annual air temperature similar with that in the other two regions (11 • C), Table S1. While the first two are vulnerable to drought, the Banat Plain, on the other hand, is vulnerable to floods [49], with drought occurring less frequently than in the other analyzed regions: 37.4% in Banat-Crisana region (including the Banat Plain) between 1961-2008, compared with 41.6% in Oltenia region or 58.4% in Muntenia region (including the Baragan Plain) [46].
The irrigation system of Romania covers an area representing 31% of the total arable land (9.4 mil. ha), according to the National Institute of Statistics (NIS). However, most parts of the system were decommissioned post-1989 (post-socialist period), only 284,000 ha (3%) being irrigated at least once in 2019. The state of irrigation in all three test areas is presented in Table 1 and is based on data from every county encompassed by each of the plains [50]. Romania has the 6th widest agriculture area of the EU countries and should have a high level of crop productivity. Even though it is one of the top ten exporters of wheat and maize worldwide, Romania has yields affected by land fragmentation, poor land reclamation systems, insufficient mechanization and human resources. The main crops grown over the three test areas are common wheat, maize, sun flower, rape, sugar beet, etc.

Data Collection
The Terra MODIS collection 6 dataset has been used: MOD16A2, level 4 and MOD09A1, level 3, which are 8 days composite syntheses produced at a 500-m pixel resolution and MOD11A2, level 3, 8 days composite products with a 1 km spatial resolution [51][52][53]. The dataset was downloaded from EARTHDATA [54].
The surface reflectance product (MOD09A1) was used to compute the NDVI for drought analysis, while the Land Surface Temperature (MOD11A2) was used for the mean land surface temperature and cloud mask during the 8-day synthesis to eliminate errors. The MODIS surface spectral reflectance products are corrected for atmospheric conditions such as gasses, aerosols, and Rayleigh scattering [55].
The Evapotranspiration and Potential Evapotranspiration products (MOD16A2) were used to compute DSI. Each layer (ET and PET) is the sum of all of the eight days within the synthesis period and the algorithm used for these products is based on the logic of the Penman-Monteith equation [56,57].
Drought analysis was performed between the 89th and 265th days of the year (DOY), i.e., the period between 29 March and 30 September, which is considered the vegetation season (a total of 23 synthesis products per year), over 2001-2019; Table S2. It should be noted that there was a lack of MODIS satellite data for the synthesis product/DOY 169 (18th-25th of June) and DOY 177 (26th of June-3rd of July) in 2001.
An auxiliary dataset used to identify agriculture areas, was CORINE Land Cover (2018), with 100 m spatial resolution (Minimum Mapping Unit: 25 ha), [43]. The SPEI was used for validation of the DSI results and is calculated from in-situ data (measured precipitation and air temperature data from weather stations). Volumetric surface soil moisture products represent the content of liquid water in a surface soil layer of 2 to 5 cm depth expressed as m 3 water per m 3 soil. The ESA CCI Soil Moisture data products [58] were used for the validation procedure (2001-2019).
The data sets used in this paper were synthesized in Table 2 and the workflow of data used to compute the DSI is summarized in the flowchart below, Figure 2.

Methods
The Drought Severity Index was computed according to the methodology developed by Mu et al. [39]-see Equations (1) and (2) below-resulting in a DSI 8 days composite at a 500 m resolution for the 2001-2019 period. The negative values of DSI were classified in three classes of drought intensity according to Table 3, which are equivalent to those used by EDO [19].
where:   [59] was computed based on surface reflectance (Red and NIR bands) extracted from MOD09A1 product. It has valid values between −1 and 1, but in this study a NDVI threshold value greater than 0.2 was used for the vegetation season analysis. Values below 0.2 designate barren soil (0-0.2 interval), as a result of works performed on the agricultural land, such as plowing or harvesting.

ET-Evapotranspiration
The NDVI, ET, PET [56,60] indices were clipped by the LST and by Romanian border. Because clouds affect the quality of data and could give some erroneous results or bare soil leads to some drought overestimation, the NDVI, ET and PET data were clipped by the positive value of LST (the LST values > 0 were considered valid, while values < 0 were considered invalid), which was resampled to 500 m in order to have the same spatial resolution as the other indices. Performing this process, the possibility to have negative LST pixels was eliminated, especially during the summer season (June-August). The lack of LST data marks the clouds so the vegetation indices values are excluded in the DSI computation.
The results obtained based on the DSI from MODIS were compared with monthly SPEI obtained from meteorological data and 8 days Soil Moisture Anomaly obtained from EEA soil moisture products.
The SPEI is an index developed by Serano et al. [11] for drought variability analysis on multiple time scales, suited to detecting, monitoring, and exploring the consequences of global warming on drought conditions [61,62]. As an extension of the widely used SPI, it takes into account both precipitation and PET. It can measure drought severity according to its intensity and duration, and can identify the beginning and end of drought episodes [63][64][65]. The PET can be estimated by several formulas (e.g., the Thornthwaite equation, the Penman-Monteith equation, the Hargreaves equation etc.), but SPEI is not linked to any particular one. In this study, the Thornthwaite equation [66,67] was used to calculate PET. Also, the negative values of SPEI were classified in three classes of drought intensity according to Table 3.
Since the SPEI is calculated from in-situ data, which represents measured points data from weather stations, SPEI values are interpolated at 1 km spatial resolution, by using Radial Basis interpolation. The reference period taken into account in the calculation of the SPEI was 1961-2019. This is the common period in which meteorological data are available at most stations in Romania.
Soil Moisture Anomaly was obtained from European Environmental Agency (EEA) soil moisture products. Soil Moisture dataset combine long-term (1978-2019) passive and active products retrieved from radiometrically calibrated backscatter or brightness temperature measurements [68]. The ACTIVE (active-microwave-based only) product merges data from AMI-WS and ASCAT (Metop-A and Metop-B), and the PASSIVE (passive-microwave-based only) product merges data from SMMR, SSM/I, TMI, AMSR-E, WindSat, AMSR2, and SMOS [69]. The SMA equation takes into account the multiannual average of the soil moisture as well as the standard deviation. The negative values of SMA are equivalent to those used by EDO [19] to establish drought intensity (Table 3).
In order to highlight the trend and its statistical significance, the non-parametric Mann-Kendall test and Sen's slope estimator were used [70]. These were calculated with MAKESENS-application [71]. Also, ANOVA statistical test was applied to check the hypothesis of equal mean (H 0 ) between agricultural areas affected by drought based on three drought indices (DSI, SPEI and SMA), against the alternative hypothesis (H 1 ) which shows that there are at least two data sets that are different. If the null hypothesis (H 0 ) is true, then any differences in the means from the agricultural areas affected by the three drought indices are explained by chance alone [72].    (Figure 5a). The differences can be explained by the fact that SPEI is calculated from interpolated meteorological data and does not take into account the soil moisture reserve, as DSI reflects drought in terms of vegetation.

Drought Frequency According to the Drought Severity Index
The frequencies of agricultural drought were between 5% and 37%, during the vegetation season, according to the DSI analysis from 2001 to 2019. The 10-20% of drought frequency annually affects the same area of about 9.3 million ha (68.8% of the country's agricultural land) (Figure 6a top), of which about: 458,000 ha in the Banat Plain (71.5% of the agricultural area of the plain), 461,000 ha in the Oltenia Plain (67.3%) and 608,000 ha in the Baragan Plain (65.7%) (Figure 6b top), Table S3. Over 20% of drought frequency annually affects the same area of about 2.8 million ha, and less than 10% occupies the same 2.4 million ha of Romania's agriculture areas. The distribution of the drought frequency (over 20%) is kept at the level of the three analyzed plains, the largest affected areas being in the Baragan Plain (220,000 ha), followed by the Banat Plain (148,000 ha) and the Oltenia Plain (138,000 ha), Table S3.
Extreme drought has the lowest frequency of all types of drought, only from 5% to 16%. The largest affected areas (over 5 mil ha) in Romania are reported in the June-August period with a frequency up to 16% (Figure 6a bottom). On the other hand, extreme drought affected during the May-June period over 390,000 ha in the Oltenia Plain and over 300,000 ha in Baragan Plain, while in Banat Plain there were affected around 200,000 ha, with a frequency up to 11% (Figure 6b bottom).

Drought Variability Over the Last Two Decades
In the last 20 years, there have been droughts and extremely droughty periods that have had a negative impact on the agricultural crops.
The   In contrast with the results at the Banat Plain, in the Oltenia Plain the extreme drought affected 57% of the agricultural area (2002), and 28.5% according to the 75th percentile, Figure 9b. Also, the drought affected 33% of the agricultural area in 2007, between 21 and 28 July ( Figure S1).

Drought Severity Index Validation
The SPEI and SMA values have been used for DSI results validation. The same extent and the same period of time have been used for comparison and validation. Thus, the differences can be  gives and also the air temperature registered high values, above average. Therefore, although the potential evapotranspiration could have had high values, it did not have a negative impact on the agricultural vegetation, and the soil registered optimal humidity values. The correlation between DSI and SPEI has a Pearson coefficient of 0.72, which means that both identify the occurrence of drought, while the correlation coefficient between DSI and SMA is 0.82.
Regarding the agricultural areas affected by drought in Banat Plain, the three indices have similar results except 2011 and 2018 when a large difference is present between DSI/SMA and SPEI (Figure 11b). In many cases, the SPEI has a higher value than DSI and SMA. Therefore, even though the meteorological drought occurred in those years, soil moisture did not decrease so as to cause agricultural drought. Based on these data, the correlation between DSI and SPEI has a Pearson coefficient of 0.57, while between DSI and SMA the correlation coefficient is 0.87.
The The results of the ANOVA statistical test applied for the agricultural areas affected by drought based on DSI, SPEI and SMA shows that the null hypothesis (H 0 ) is true for all the study areas (Table S4). Therefore, any differences in the means of the agricultural areas affected by the three indices drought are explained by chance alone. The three data sets that were compared have the same mean values. There are no differences in the results of the three drought indices.
The negative effects of drought in Romania can be ascertained by the low yields of certain crops, such as wheat and maize. Thus, the average wheat and maize yields in Romania in the 2001-2019 period were 3215.9 kg/ha and 4057.6 kg/ha [54]. According to the data from the NIS, the wheat production in Romania registered a decrease of up to 55 (Figure 12d). It can be seen that the lowest productions were registered mainly in the former part of the analyzed period (2001-2010), the latter part being characterized by above average yields.

Drought Severity Index Trend According to the Mann-Kendall Test
According to the Mann-Kendall test, applied on the values of the DSI (2001-2019), the trend is positive, increasing the values of the index, between 30 March-27 July and negative between 28 July and 29 September, both at the level of the three analyzed fields and in Romania overall (Table 4). This means that in the former part of the year the drought phenomenon tends to decrease, while in the latter part it tends to increase.
The trend of DSI is for statistically significant growth at a level of 0.001 (99.9% confidence level) at the beginning of May in the Baragan Plain and at the end of April-early May for the agricultural area of Romania. Also, a statistically significant increase to a level of 0.01 (99.0% confidence level) is recorded in April and May in the Baragan Plain, end of April-early May in the Banat Plain and in May in the Oltenia Plain.
In contrast to the positive trend of DSI, the negative trend has a level of statistical significance of 0.01 (99.0% confidence level) in September and 0.1 (90.0% confidence level) in August, in the Baragan Plain. In the other analyzed areas, there were recorded statistically significant decreases of the DSI values to a level of 0.05 (95.0% confidence level) and 0.1 (90.0% confidence level), respectively.   The increasing tendency of the areas affected by the drought is statistically significant at a level of 0.001 for the moderate drought in the Baragan Plain and at a level of 0.01 and 0.1 for the severe drought, respectively extreme drought (Table S5). In contrast, in the other two plains the increasing tendency is significant at a level of 0.05 (moderate drought) and 0.1 (severe drought), the rest being a statistically insignificant trend. The decreasing tendency of the areas affected by the drought is statistically significant for the Baragan Plain and the Oltenia Plain for all types of drought, while in the Banat Plain the trend is significant only for moderate and extreme drought at a level of 0.05 and 0.1.
The trend of the SPEI for the period 2001-2019, is similar to that of the DSI, being statistically significant decrease at a level of 0.01 in the Baragan Plain and statistically insignificant in the other plains (Table 5). Also, the decreasing trend of SMA in the second part of the vegetation period is statistically significant for the three plains analyzed, at a level of 0.01, 0.05 respectively ( Table 6). Unlike DSI, soil moisture has a decreasing trend in the first part of the analyzed period.

Discussion
Two decades have passed since the MODIS sensor is being collecting daily information at the Earth's surface. This makes it possible to calculate some statistical parameters that are necessary for DSI analysis such as: multiannual average and standard deviation, frequency but also the trend.
The performance of the DSI index can be comparable to that of the SPEI and SMA indices, the intensity of the drought being identified as accurately as possible in time and space. Based on the analysis performed we were able to identify some strengths but also some limitations regarding its use. Thus, the strengths are: • The DSI takes into account the evapotranspiration, being one of the few remote sensing indices that take into account this parameter in drought analysis; • The values of DSI are standardized, therefore the results may be comparable to those from other indices; • DSI is useful to identify the spatial extent of agricultural drought, at 500 m spatial resolution, over different time intervals (weekly, monthly, seasonal etc.); • It can be applied globally and also locally.
The limitations are: • Cloud cover (especially for daily products) and bare soil (it can offer wrong information-overestimate the drought). The cloud cover is a very well-known problem of the optical images; thus, 8-days synthesis and cloud mask filter have been used.

•
Vegetation responds differently to the drought according to phenological phase, being necessary to identify the type of vegetation. The spatial resolution of MODIS data (500 m) does not allow the vegetation identification by species; therefore, limitation could not be solved. • MK test cannot be applying for DSI on the spatial scale due to existing pixels with missing value.

•
The short DSI samples (19), used in this case study, are not enough for the MK test.
The drought results, in this study, were influenced by the lack of information caused by the clouds and the state of the agricultural land before sowing and after harvesting. Therefore, the analysis of drought over the agricultural land was limited to those pixels without clouds and filtered by an NDVI threshold. According to the surface affected by clouds and filtered by the 0.2 NDVI threshold during all of the 435 synthesis images for each study area, it was noticed that a wide agricultural surface was removed from analyses: around 12.6% (116,571 ha agricultural land in Baragan Plain), 15.1% (100,506 ha agricultural land in Banat Plain) and 12% (82,227 ha agricultural land in Oltenia Plain).
At annual level (average of the 23 synthesis), in the Baragan Plain and Oltenia Plain, the highest cloud cover was registered in 2019 around 16.5%, and the lowest, in 2005 and 2007 with 8% in Baragan Plain and in 2011 in Oltenia Plain (9.1%). In the Banat Plain, the highest cloud cover was registered in 2013 (20.4%), and the lowest in 2005 (10.7%).
However, the results obtained from DSI highlight the droughts of the last 20 years. The dry periods and spatial extent identified in this study are consistent with those identified based on other indicators in previous studies [46,73]. However, Southern, Southeastern and Eastern Romania are the main agriculture areas affected by drought. The occurrence, duration and intensity of the drought are different from one region to another due to natural conditions (geographical position, relief, soil, climate etc.), agricultural measures and farm practices, as mentioned in previous studies [29,48,49].
The ANOVA test was applied to check the hypothesis of equal mean between agricultural areas affected by drought based on DSI, SPEI and SMA and the result was that there were no differences between the three data sets or there is insufficient evidence to deny the null hypothesis. Thus, the null hypothesis (H0) was confirmed. However, there are some differences between those obtained on the basis of DSI, SPEI and SMA. The spatial extent (%) of Romanian agricultural areas affected by drought, based on SPEI data, is greater than the one displayed by DSI, which favors possibility of an overestimation of drought when the meteorological index is used, to the detriment of satellite vegetation index. Higher drought estimation, obtained from SPEI data, was also observed by Um et al., 2018 [74], who used MODIS and AVHRR data in order to calculate and revise the DSI in East Asia.
The spatial droughts areas are not expected to match each other exactly each time, since the SPEI and DSI measure different types of drought. While the SPEI is used to measure meteorological drought, as a result of precipitation and potential evapotranspiration, the DSI measures vegetation, agricultural drought, based on the NDVI and actual and potential evapotranspiration [74]. On the other hand, the results from DSI are much closer to those of SMA as expected, the agricultural drought being strongly influenced by the soil moisture reserve.
Compared to DSI, the SPEI grids and SMA are not affected by clouds. Therefore, the extend of agricultural land analyzed remains the same. Despite the fact that SPEI does not depend on cloud cover, the results are considerably influenced by the number of in-situ measurements, which do not continuously monitor the drought through satellite information. Meanwhile, the results of the SMA depends on the resolution of the input data.
Regarding the DSI trend, it is similar to the SPEI and SMA indices. The trend of the drought phenomenon decreases in the first part of the year, while in the second half of the year the trend increase. This trend is confirmed by other studies conducted on Romania.
According to Busuioc et. al. (2016) [75], the increasing trend in the spring rain shower frequency in Romania, between 1961 and 2010, could be due to more instability events covering Romania, accompanied by more frequent cyclonic structures and more precipitable water. Also, the atmospheric instability is accompanied by changes in atmospheric circulation in spring (more frequent southern circulations transporting humid air masses from the Mediterranean Sea to Romania) and precipitable water in the case of summer that shows an upward trend over Romania.
Compared to the European trends of the mean air temperature that increased everywhere, in Romania decreases, according to the analyses made by Spinoni et al. 2016 [76], and precipitation shows the increasing trend in northern and the decreasing trend in the southern of Europe. As well, in northern and eastern European countries, drought severity shows a decrease linked to increase PET (excluding Romania) and precipitation.

Conclusions
The Drought Severity Index is a useful tool for monitoring and evaluating the drought phenomenon. The performance of the DSI in Romania is comparable to that of the SPEI and SMA indices. Negatively influencing the quality of life and the environment, drought can become a long-term phenomenon that can lead to arid agricultural lands and, in the worst case, desertification. The identification of the occurrence of the drought and its severity can be achieved through a complex analysis by the combined use of the data sets from in-situ measurements and from satellite images.
DSI capability was proved for the agricultural drought monitoring in Romania and especially in the three investigated agricultural areas: Banat Plain, Oltenia Plain and Baragan Plain. Thus, according to the results, the agricultural lands more affected by drought are located in the South, Southeast and East of Romania, mostly in 2002Romania, mostly in , 2003Romania, mostly in , 2007Romania, mostly in and 2012. The results obtained from DSI were confirmed by the low yields in the dry years.
Although the affected agricultural area was higher in 2007 than in 2012, the drought duration in 2012 was longer than in 2007. 25% of the period analyzed in 2012, the drought phenomenon occurred which affected over 50% of the agricultural area.
Results indicated that the Baragan Plain is the most vulnerable region to agricultural drought, the intensity and frequency of the phenomenon being strongly influenced by its climatic characteristics. After Baragan Plain, the Oltenia Plain is the second most affected and vulnerable region to drought, while the Banat Plain is the least affected by the drought effects.
The drought trend in the last two decades is decreasing in the former part of the vegetation season and increasing in the latter part. Thus, an enhancement of drought is expected at the end of the vegetation season, which may have negative effects on next year's agricultural vegetation, such as autumn wheat and rape.