1. Introduction

Droughts are recurrent natural disasters characterized by a significant rainfall deficit, leading to water scarcity, streamflow depletion, crop damage, and groundwater reduction. They can be classified as meteorological, hydrological, agricultural, and socio-economic droughts (Dracup et al. 1980; Wilhite, Glantz 1985). Various parameters must be considered when defining the temporal and spatial scales and the regional characteristics of droughts (Tallaksen, Lanen 2004). More than 130 published drought definitions can be categorized in several ways, as noted by Wilhite and Glantz (1985). Some of the most commonly used definitions of droughts have been listed by Tate and Gustard (2000) and Demuth and Bakenhus (1994). Droughts have negative impacts on social life, the economy, and the environment, although these impacts are difficult to detect. Several indices have been developed to identify and analyze droughts, of which the most commonly used one is the standard precipitation index (SPI) (McKee et al. 1993; Ganguli, Reddy 2014). Other indices, such as the Palmer drought severity index (PDSI) (Palmer 1965), the surface water supply index (SWSI) (Shafer, Dezman 1982), and the crop moisture index (CMI) (Palmer 1968), are also widely applied to detect drought events.

Several drought trend studies have been conducted all over the world (Piccarreta et al. 2004; Xu et al. 2011), and the hydrological time series of precipitation, temperature, soil moisture, and evapotranspiration at the catchment scale have been widely investigated (Villarini et al. 2009; Burn et al. 2010; Gocic, Trajkovic 2013; Zhang et al. 2015). As a consequence of global warming and climate change, water resources, the environment, agricultural activities, and industrial production have been significantly affected (Shi, Xu 2008). However, the detection of changes in climate and hydrological time series is complex and challenging but has become an important issue because of the role of such changes in water resources management and drought analysis. According to the report of the Intergovernmental Panel on Climate Change (Houghton et al. 1996), in the 20th century, global land precipitation has increased by 2% on average. The impact of climate change on the maximum and minimum temperatures of Gombe City in northeastern Nigeria has shown an increasing trend (Alhaji et al. 2018). Asfaw et al. (2018) investigated trends in the precipitation and temperature parameters for the Woleka Basin in northern Ethiopia and found that the number of dry years increased. Although the mean and minimum temperatures presented an increasing trend, an insignificant trend was observed for the maximum temperature values. Keskin et al. (2015) studied monthly and annual water levels in Lake Eğirdir in Turkey, using Mann-Kendall and Sen’s slope tests. Based on their results, the annual water level was reduced by approximately 0.026 m. Karabulut et al. (2008) analyzed trends in rainfall and temperature data in the Samsun region between 1931 and 2006 and noted statistically significant trends in summer temperatures. Although there are numerous investigations on trends in hydrometeorological parameters and climate change (Cannarozzo et al. 2006; Yu et al. 2006; Kampata et al. 2008; Zhong, Li 2009; Kumar, Jain 2010; Santos et al. 2011; Yuce et al. 2015, 2018), only few studies have performed trend analyses in drought events (duration and severity) (Tabari et al. 2012; Ganguli, Reddy 2014; Spinoni et al. 2014; Dashtpagerdi et al. 2015; Zhang et al. 2015; Yuce et al. 2022).

As of the beginning of 2021, most of Turkey is facing a severe drought. An increase in the number of dry days is inevitable after several low-rainfall seasons. In 2019, the summer and autumn months were largely devoid of rain, resulting in decreased reservoir water levels. The year 2020 was the driest of the last 5 years, with particularly little rainfall in the latter half of the year, according to a NASA report1. This study will investigate the presence of any systematic changes in extreme drought events, using historical hydrometeorological data in the Ceyhan River Basin, which plays a critical role in the agricultural and hydropower production of the country. If the presence of any trend is detected, its magnitude and change point in the time series will be studied in detail. Furthermore, the trends of extreme drought events, such as the annual maximum drought severity (AMDS), defined as the largest cumulative severity value for each year, and the annual maximum drought duration (AMDD), which is described as the length of the maximum drought duration for each year, will be investigated. Notably, these two parameters are independent of one another. Here, the length of the AMDS may not be the same as the length of the AMDD within 1 year. Hence, AMDS and AMDD will be evaluated as different time series in statistical analysis tests. Both parameters of the extreme events are calculated by using the SPI method for multiple-time scales of 1, 3, 6, 9, and 12 months.

2. Study area and data

The Central Taurus Mountains are located in the northern and northwestern parts of the Ceyhan River Basin, with two different mountain ranges, namely the Bolkar Mountains and the Tahtalı Mountains starting from the west. The Ceyhan River Basin is situated in the Mediterranean Climate Zone, which is characterized by dry summers and mild, wet winters. The rainiest months are December, January, and February, and the driest months are June, July, and August. The total length of the Ceyhan River is around 425 km, with an annual discharge of 82.9 m3/s and a basin yield of 10.7 L/h/km3. The maximum total annual rainfall of the catchment is recorded in Kozan with at least 842 mm, and the lowest total annual rainfall is recorded in Elbistan with 395.7 mm (Eris et al. 2019; Yuce et al. 2019). The highest average annual temperature is 19.3°C in Kozan, whereas the lowest average annual temperature is 8.9°C in Göksun (Uzunkol, Kızılelma 2016).

Fig. 1.

Locations of the meteorological stations in the Ceyhan Basin.

https://www.mhwm.pl/f/fulltexts/154573/MHWM-9-0011-g001_min.jpg

For this study, we used long-term monthly precipitation data. To evaluate the SPI index, the precipitation data of 11 out of 23 meteorological stations were obtained from the General Directorate of Meteorology of Turkey (MGM), whereas the data of the other 12 stations were obtained from the General Directorate of the State Hydraulic Works of Turkey (DSI). The locations of 23 meteorological stations are displayed in Figure 1. The selected monthly precipitation time series with recorded lengths ranging from 15 years (station 7767) to 54 years (station 17255) for the period from 1963-2016 were employed in the analysis. This is a significant amount of data for extreme drought investigation in this basin. The statistical characteristics determined from the monthly precipitation time series for each gauging station are presented in Table 1. The mean, standard deviation, coefficient of variation (Cv), coefficient of skewness (Cs), and lag-one autocorrelation coefficient (r1) were evaluated by using the observed monthly precipitation time series.

Table 1

Statistical parameters of observed rainfall stations; S.D., standard deviation; Cv, coefficient of variance; Cs, skewness; r1, kurtosis.

StationEarliest record yearLatest record yearLatitudeLongitudeMean (mm)St. Dev.CvCsr1
77671998201237.357536.090765.6455.530.851.060.26
82751964199537.090136.309485.7583.870.981.280.39
172551963201637.576036.915060.3965.571.091.270.50
173551986201637.102136.253969.1063.750.921.120.33
178661963201138.024036.482351.1946.290.901.070.39
178681970201138.240536.919035.3131.160.880.990.34
178701963201138.203837.198233.3028.870.870.970.33
178711993201137.786737.653260.1566.401.101.550.46
179081963201137.433735.818870.2458.790.841.110.23
179601964201137.015335.795559.0058.640.991.430.33
179791964201136.768735.790367.5470.991.051.550.33
D20M0011963200938.300037.583329.1326.700.921.120.31
D20M0021980200537.545836.774758.4066.961.151.510.42
D20M0041970198837.383337.166641.1541.771.020.960.42
D20M0061963197437.166635.700053.9860.211.121.850.30
D20M0091963200937.883336.850078.8981.621.031.280.44
D20M0111963199638.616636.933332.3229.730.921.280.27
D20M0131968201537.466637.250050.2854.041.071.360.41
D20M0141973201537.335236.005690.2778.250.871.150.25
D20M0151976201537.050036.066676.1566.450.871.120.28
D20M0161977201537.500035.833378.2769.640.891.250.25
D20M0171980200338.152036.465465.4865.481.001.430.33
D20M0181982200238.291836.760425.0928.131.122.290.19

3. Methodology

3.1. Standardized Precipitation Index

The SPI method, developed by McKee et al. (1993), is used to evaluate the level of deficiency in precipitation on different time scales ranging from 1 to 48 months. It is the most well-known and used index among many other indices to identify meteorological drought events and is based on fitting precipitation data to a probability distribution function (PDF). Gamma is the most suitable probability density distribution function for climatological data (Ganguli, Reddy 2014; Zhang et al. 2015). The dry period is considered to be the period in which the index is less than or equal to –1. The SPI, which is classified based on values presented in Table 2, is determined by Equation 1:

(1)
SPI=XijXimσ

where Xij is the monthly precipitation amount and Xim and σ are the mean and standard deviation of precipitation evaluated from the all-monthly time series, respectively.

Table 2

SPI classification (McKee et al. 1993); SPI = Standardized Precipitation Index.

SPI valueCategory
≥ 2Extremely wet
1.50 to 1.99Very wet
1.0 to 1.49Moderately wet
0.99 to 0Normal
0 to 0.99Near normal
–1.00 to –1.49Moderately dry
–1.50 to –1.99Severely dry
≤ –2Extremely dry

3.2. Homogeneity test

3.2.1. Wallis-Moore and Wald-Wolfowitz

The Wallis and Moore phase frequency test (Wallis, Moore 1941) is used to detect deviations of time series for randomness in the sequence of values. The test is based on sign differences (- or +), while the first and last phases named sequence of signs are not accounted thus the number of phases is identified. If n ≥10 and continuity correction is applied, a fairly good test may be based on the hypothesis that the data are normally distributed, when n ≥25, the correction is not used (Wallis, Moore 1941). The z test statistic is calculated as follows:

(2)
E[H]=13(2n7)
(3)
var[H]=190(16n29)
(4)
z=(HE[H])/(190(16n29))

where h is the number of phases, although the first and last phases are not taken into account. The z-statistic is normally distributed. A continuity correction of –0.5 is added to the denominator for n ≤ 30.

The Wald-Wolfowitz test, which is also known as Runs test, measures the randomness of the data and examines whether an observation influences the subsequent observation. Time series data are cut from a certain level that could be the mean, medium, or mode to determine if each value in the series is lower or higher than this level. The number of passes from one data to the other above or below a certain level is called the run number, it is small if it is below or above long periods. Such series may not have homogeneity (Wald, Wolfowitz 1940). The result of the test is z, the number of data N, the number of runs r, the number of values below the medium level Na, the number of values above the medium level Nu; the equations are as follows:

(5)
E[r]=2Na+Nu(Na+Nu)+1
(6)
V[r]=2Na+Nu(2Na+NuN)(N1)N2
(7)
z=r1E[r]V[r]

3.3. Trend Analysis

3.3.1. Mann-Kendall Test

The Mann-Kendall test is a non-parametric test to determine the trend in the hydrometeorological variable in a time series (Mann 1945; Kendall 1975; Helsel et al. 2002). It is the most preferred statistical method for time series since it is not fitted to any distribution. The application of the Mann-Kendall test statistic Z is expressed as follows:

(8)
S=k=1n1j=k+1nsgn(xjxk)

where n is the number of the data, xj and xk are the data point in years j and k (j > k), and ti is the length of the tied rank group.

(9)
sgn(xjxk)={1(xjxk)>00(xjxk)=01(xjxk)<0
(10)
Var(S)=n(n1)(2n+5)irti(ti1)(2ti+5)18
(11)
Z={s1Var(S)S>00S=0s+1Var(S)S<0

A positive Z value indicates an increasing trend, whereas a negative value indicates a decreasing trend. Critical test statistical values are 1.645, 1.97, and 2.57 for different significance levels of 90%, 95%, and 99%, respectively (Yu et al. 1993).

3.3.2. Spearman’s Rho Method

Spearman’s rho test is a commonly employed nonparametric method to investigate the presence of trends (Lehmann, D’Abrera 1975; Sneyers 1990). Its statistics rs and z (standardized test statistics) are computed as follows:

(12)
rs=1[6i=1n(Rxii)2)]n(n21)
(13)
z=rsn21rs2

where Rxi (rank statistic) is obtained by sorting the data, and n is the length of the time series. Negative z values show decreasing trends, whereas positive z values indicate increasing trends. At the 90% significance level, for z > ±1.645, the null hypothesis of no trend is rejected.

3.4. Trend Slope

3.4.1. Linear Regression Method

Regression analysis is a statistical method used to predict the relationships between one dependent variable and one or more independent variables. It forms the basis of complex estimation methods and is used to evaluate the strength of the relationship between variables and to model the relationship between them. The dependent variable is usually represented by Y, and the independent variables are denoted by X (Gocic, Trajkovic 2013). The linear regression formula is given as follows:

(14)
Y=a+bX

where a is the intercept, and b is the slope of the line. A positive b value indicates increasing trends, and a negative a value indicates decreasing trends.

3.4.2. Sen’s Slope Estimator

Sen’s slope test (Sen 1968) is a non-parametric test that calculates the slope of the trend in a data set; it is used for equiponderant time series. For each data point, the slope difference is calculated per changing time. The slope of the trend can be estimated by the median of all slopes between data pairs in the same season (Helsel et al. 2002). All slope pairs are ranked from smallest to largest, and if the calculated number of slopes (n) is odd, the median slope gives the slope S. If n is even, the two median slopes are averaged. Here, Q represents data, n is the length of the data, and T is time. The slope of n pair of data is predicted as follows:

(15)
S=Q2Q1T2T1
(16)
S={Sn+12n=evenSn2n=odd

3.5. Tests for Change Point Detection

Considering the importance of climatic variability in terms of water availability, increasing irrigation demand, crop yields, and other factors, it is crucial to identify and evaluate the trends in the monthly, seasonal, and annual historical data series. Change point detection is substantial to evaluate the period in which a significant change occurs in a time series. In this study, the standard normal homogeneity test (SNHT), Buishand’s range test (BRT), and Pettitt’s test (PT) were employed to detect change points in the time series.

3.5.1. Standard Normal Homogeneity Test

The statistic of the SNHT (Alexandersson 1986), (Tk), is used to compare the average of the first n year with the average of the last (nk) year with n data points (Vezzoli et al., 2012; Jaiswal et al., 2015). The Tk equation is as follows:

(17)
Tk=kZ12+(nk)Z22

Z1 and Z2 can be calculated as follows:

(18)
Z1=1ki=1k(XiX¯)σX
(19)
Z2=1nki=k+1k(XiX¯)σX

Here, X¯ and σX are mean and standard deviation, respectively. The year in which Tk reaches the maximum value is considered as the point of change. The significant critical values for SNHT are given in Table 3.

Table 3

Critical values for different change point detection test statistics.

Number of ObservationsBuishand Range TestPettitt TestStandard Normal Homogeneity Test
N1%5%10%1%5%1%5%10%
201.61.431.171579.566.956.07
301.71.51.1213310710.457.656.6
401.741.531.1320816711.018.16.92
501.781.551.1429323511.388.457.15
701.811.591.1748839311.898.87.47
1001.861.621.2284167712.329.157.77

3.5.2. Buishand´s Range Test

The adjusted partial sum, Sk, is computed as follows:

(20)
Sk=i=1k(XiX¯)

When the series is homogeneous, the value of Sk fluctuates around zero. The test captures data showing sensitivity to break in the middle of distorting the homogeneity of the time series. If there is a change in K year, it reaches the maximum or minimum value around the year k = K (Buishand 1982). The significant change can be computed by evaluating the rescale-adjusted range R:

(21)
R=Max(Sk)Min(Sk)X¯

3.5.3. Pettitt’s Test

PT is a non-parametric method (Pettitt 1979) to determine change detection in a time series. It is commonly used to calculate the occurrence of an unexpected change in hydrometeorological records (Sneyers 1990; Smadi, Zghoul 2006). The null hypothesis indicates that the series has an independent and random distribution, whereas the alternative hypothesis indicates a sudden change. The Uk test statistic is expressed as follows:

(22)
Uk=2i=1krik(n+1)k=1,n

The Uk test results are presented graphically. In case of a break in k = K year, the test statistic is at maximum or minimum values. Significant critical values are given in Table 3.

(23)
K=max|Uk|

4. Analysis of Results

4.1. Homogeneity Analyses

The precipitation data of 23 meteorological stations in the Ceyhan River Basin were analyzed by the SPI drought index for time scales of 1, 3, 6, 9, and 12 months. The AMDS and AMDD time series for station 17355 are illustrated in Figure 2. These graphs do not show any statistically significant trend in the data and provide only a rough view of the variation of the data over time.

The Wallis and Moore phase-frequency test and the Wald-Wolfowitz test were employed for the analysis of extreme drought events (AMDS and AMDD), verifying the homogeneity for different time scales. The results of these two tests were compared to achieve an accurate homogeneity and randomness in the time series. If the outcomes of both tests suggest homogeneity in the time series, the presence of homogeneity is accepted. When only one of the methods suggests homogeneity of the time series, the results are defined as doubtful; if both methods suggest no homogeneity in the data set, it is labeled as non-homogeneous. The homogeneity test results for the AMDS and the AMDD time series of 23 meteorological stations for time scales of 1, 3, 6, 9, and 12-month are presented in Table 4.

In the AMDS time series, for SPI 1-month time scale, the data of 19 out of 23 stations indicated homogeneity by both tests, whereas stations 17871, 17979, and D20A011 presented homogeneity by one of the methods. Only station 17908 showed non-homogeneity by both tests for the significance level of 90%. For the same time scale, stations 8275, 17355, 17870, 17908, and D20A016 presented homogeneity by one of the methods, and only station D20A006 illustrated non-homogeneity by both methods; the data for the remaining stations were homogeneous in the AMDD time series. The AMDS time series were homogenous in all stations for the SPI 3-month time scale, apart from station 17960, which presented homogeneity by the Wald-Wolfowitz test. Although the AMDS time series presented homogeneity by one of the methods at 10 stations (17866, 17870, 17960, 17979, D20A002, D20A009, D20A011, D20A013, D20A015, and D20A016), the data at only one station (8275) showed non-homogeneous characteristics. The data at the other stations showed homogeneous features.

Fig. 2.

The extreme drought events´ both AMDS and AMDD time series of station 17355 at different time scales.

https://www.mhwm.pl/f/fulltexts/154573/MHWM-9-0011-g002_min.jpg

The SPI 6-month time scale AMDS time series presented homogeneity by one of the methods at stations 8275, 17868, and D20A002, whereas the other stations indicated homogenous data sets. The AMDD time series displayed non-homogenous features at stations 8275 and D20A013, homogeneous features by one of the methods at stations 17868, D20A016, and D20A017, and the data sets at the other stations were homogeneous by the two tests. With the increase in the length of the time scales, the number of non-homogenous stations also increased. For instance, three and five stations were non-homogenous for the SPI 9-month time scale and the SPI 12-month time scale, respectively. Station 17868 in both AMDS and AMDD as well as stations D20A001 and D20A016 in the AMDD time series were non-homogenous for the SPI 9-month time scale. Station 17868 in AMDS, stations D20A011 and D20A013 in both AMDS and AMDD, and stations D20A009 and D20A016 in the AMDD time series were non-homogenous for the SPI 12-month time scale.

Table 4

Homogeneity test results for extreme drought events at different time scales.

StationTestsSPI 1SPI 3SPI 6SPI 9SPI 12
AMDSAMDDAMDSAMDDAMDSAMDDAMDSAMDDAMDSAMDD
7767WM-0.1090.762-0.1091.4150.7620.1090.1091.4150.7620.762
WW-0.0231.0230.601-0.3691.6071.5132.2461.8152.1611.714
ResultHHHHHHDDDD
8275WM1.7273.4530.8633.4532.5901.7271.7271.7271.2950.432
WW0.3871.2380.6712.8871.3311.7100.1030.837-0.0110.608
ResultHDHN-HDN-HDDHH
17255WM0.8750.5470.4380.4380.2190.4380.2190.2190.8750.875
WW0.8750.5470.4380.4380.2190.4380.2190.2190.8750.875
ResultHHHHHHHHHH
17355WM0.5852.3410.7320.1460.2931.0240.1461.0240.2931.499
WW0.4550.0270.3090.4380.8581.0390.0800.450-0.4760.759
ResultHDHHHHHHHH
17866WM1.9571.1511.2662.1870.1151.4960.8060.8060.4601.496
WW-1.505-1.232-0.1840.6651.3811.2441.2831.9301.4951.289
ResultHHHDHHHDHH
17868WM0.4990.6241.3721.7460.6241.3722.1202.1202.1201.746
WW0.0990.3121.2360.9773.1921.9613.3022.0383.3481.585
ResultHHHHDDN-HN-HN-HD
17870WM1.2662.5320.5752.8770.8061.4960.8060.8061.8411.841
WW-0.988-0.854-0.2930.0480.2170.8390.3480.6030.9561.540
ResultHDHDHHHHDD
17871WM0.0950.0950.0950.4740.667-0.0950.667-0.0950.4770.477
WW-2.218-1.781-1.531-1.103-0.3790.6030.9100.7472.0952.133
ResultDDHHHHHHDD
17908WM2.3020.9060.8061.4960.2300.2300.4601.1510.4600.230
WW-2.217-1.924-1.293-1.062-0.1660.0880.4520.5110.7420.336
ResultN-HDHHHHHHHH
17960WM0.4650.2331.8612.1600.1160.1160.1160.5820.9310.233
WW-1.060-0.731-1.440-0.9410.1160.1161.1370.9892.2421.343
ResultHHDDHHHHDH
17979WM1.8610.2330.8142.1550.4650.4650.1160.9310.4651.603
WW-1.273-1.320-1.346-0.7900.376-0.4710.3580.2291.8291.483
ResultDHHDHHHHDH
D20A001WM1.1631.6290.4650.9310.2331.6291.6292.3271.6291.214
WW-1.294-1.493-1.297-0.574-0.2761.2310.1661.9451.0982.059
ResultHHHHHHHN-HHD
D20A002WM0.7230.7230.7231.6880.7230.7230.7230.7230.2411.741
WW-1.070-0.054-0.1420.2241.8440.1991.8671.1892.3660.952
ResultHHHDDHDHDD
D20A004WM0.4770.4770.0950.4770.4771.6211.6211.6210.0950.477
WW0.299-0.257-0.0970.1290.9701.2651.5961.6850.9270.960
ResultHHHHHHHDHH
D20A006WM-0.1242.8490.1241.6100.1241.6101.6101.6100.1240.124
WW-1.027-1.976-0.9510.3600.104-0.4910.1420.0530.091-0.351
ResultHN-HHHHHHHHH
D20A009WM0.3530.3531.0591.7640.3531.0591.7641.2821.0592.470
WW0.7070.005-0.6290.2440.4841.6100.6752.1660.7632.414
ResultHHHDHHDDHN-H
D20A011WM1.8120.5580.9752.4430.9750.9750.1391.6651.8122.648
WW-0.1900.403-1.091-1.2330.3820.5390.8221.1631.7192.297
ResultDHHDHHHDN-HN-H
D20A013WM0.2330.9310.5821.9780.5822.6760.5821.9781.9781.978
WW-0.6880.043-0.4910.2531.1801.7291.7381.2732.0212.414
ResultHHHDHN-HDDN-HN-H
D20A014WM0.5621.6010.6160.9851.3550.6160.8621.6010.6160.862
WW-0.977-0.156-0.549-0.235-0.3710.349-0.2011.510-0.0611.018
ResultHHHHHHHHHH
D20A015WM1.0241.2790.6400.1280.6400.1280.1280.6400.6400.128
WW-0.954-0.985-0.952-1.950-0.1890.732-0.1921.006-0.1031.576
ResultHHHDHHHHHH
D20A016WM0.5191.8151.0371.8150.6481.4261.4262.2041.4262.204
WW-0.636-0.255-0.496-0.5441.6281.6991.7532.1841.4822.095
ResultHDHDHDDN-HHN-H
D20A017WM0.4201.5940.4200.0841.0912.6020.0840.0840.0841.091
WW-0.0440.614-0.176-0.1870.1871.160-0.375-0.120-0.3580.290
ResultHHHHHDHHHH
D20A018WM0.0901.1730.4511.1730.4510.0901.1730.0901.715-0.090
WW-0.095-0.103-0.592-0.373-0.755-0.3620.4100.0561.1510.282
ResultHHHHHHHHDH

H: Homogeneous, D: Doubtful, N-H: Non-homogeneous

4.2. Trend Analyses

Trend detection analysis of extreme drought event time series under the influence of climate change is critical to propose and make use of management strategies. Table 5 summarizes the outcomes of the Spearman Rho and Mann-Kendall statistical tests used to detect any trends in the SPI 1-,3-,6-,9-, and 12-month time series. The null hypothesis is not rejected at the 90% significance level where –1.64 < z < 1.64. The Mann-Kendall and Spearman Rho methods performed similarly in the analysis of trend detection in the extreme drought event time series (Table 5). For the SPI 1-month time scale, a significant trend was not detected by both tests in the extreme drought events, except in the data of stations 17868 and D20A014. Although two tests showed a statistically significant increasing trend in the AMDS series of station17868, the AMDD series of station D20A014 showed a statistically significant decreasing trend. In the SPI 3-month time scale, the statistically significant increasing trend was evaluated only in the AMDS series of station 7767 and in the AMDD series of station 8275. In the case of the SPI 6-month time scale, the time series at four stations showed statistically significant trends in both test results. The AMDS and AMDD time series at station 7767 showed increasing and decreasing trends, respectively. Statistically significant increasing trends were observed in the AMDS time series of stations 8275, 17866, and D20A016, whereas decreasing trends were noted in the AMDS time series of station D20A016, with a significance level of 90%. As in the SPI 3-month period, similar results were clearly seen for the SPI 9-month time scale of stations 7767, 8275, and D20A016. In the time series at five stations, statistically significant trends were determined for the SPI 12-month time scale. Increasing trends were perceived in the AMDS time series of stations 7767 and D20A018, whereas decreasing trends were detected in the AMDS time series of stations 8275 and D20A016. Furthermore, both tests demonstrated a statistically significant increasing trend in the AMDD time series of stations 8275 and D20A009.

Table 5

Trend detection results for extreme drought events.

StationTestsSPI 1SPI 3SPI 6SPI 9SPI 12
AMDSAMDDAMDSAMDDAMDSAMDDAMDSAMDDAMDSAMDD
7767MK1.683-1.0931.683-0.8552.474-2.4482.428-1.9102.138-1.706
SR1.537-1.1311.684-0.9742.499-2.4972.160-1.9351.918-1.628
ResultNNINIDIDIN
8275MK-0.2760.685-0.5841.755-1.6741.947-1.3841.704-1.8601.874
SR-0.2970.697-0.5341.691-1.5951.787-1.4091.693-1.7621.957
ResultNNNININIDI
17255MK0.149-0.2650.2610.0080.351-0.1960.441-0.5650.211-0.479
SR0.143-0.2730.4110.0790.425-0.2590.478-0.5060.179-0.419
ResultNNNNNNNNNN
17355MK-1.1050.283-0.9521.211-0.3910.242-0.2550.377-0.3070.309
SR-1.2030.369-1.0141.123-0.3700.066-0.2340.345-0.4740.251
ResultNNNNNNNNNN
17866MK-0.0600.894-0.6810.755-0.5341.712-1.0361.050-1.5291.432
SR-0.0630.897-0.6740.787-0.6631.901-1.1511.086-1.6161.440
ResultNNNNNINNNN
17868MK1.767-1.0190.900-0.3730.2930.000-0.1630.0870.043-0.515
SR1.667-1.0121.009-0.5120.2820.232-0.1790.153-0.020-0.547
ResultINNNNNNNNN
17870MK0.4480.3230.509-0.7310.5520.409-0.2760.7720.0000.529
SR0.4590.2430.413-0.6670.4980.426-0.3330.8370.626-0.984
ResultNNNNNNNNNN
17871MK1.469-0.1840.3500.6060.0000.748-0.7710.634-0.7030.958
SR1.236-0.2570.0550.582-0.0450.816-0.7600.823-0.8001.065
ResultNNNNNNNNNN
17908MK0.620-0.2541.121-1.2430.9920.452-0.1730.564-0.4760.923
SR0.6200.1111.228-1.3230.9460.469-0.1040.531-0.4630.836
ResultNNNNNNNNNN
17960MK-0.4530.0090.604-0.5790.871-0.7290.792-0.6260.561-1.433
SR-0.5130.0450.053-0.4930.791-0.7410.696-0.5670.562-1.403
ResultNNNNNNNNNN
17979MK-1.1290.256-0.044-0.4690.836-1.0440.818-1.1170.873-0.777
SR-1.1370.2640.200-0.4450.852-1.0000.655-1.1850.938-0.689
ResultNNNNNNNNNN
D20A001MK0.756-0.3470.0530.0000.2310.072-0.044-0.134-0.2500.467
SR0.831-0.3250.1590.1550.2160.126-0.010-0.185-0.2210.400
ResultNNNNNNNNNN
D20A002MK-1.058-0.206-1.236-0.246-0.243-0.1780.133-0.621-0.4410.155
SR-1.250-0.180-1.280-0.203-0.456-0.107-0.175-0.472-0.0240.155
ResultNNNNNNNNNN
D20A004MK1.295-0.6160.9100.1460.560-0.5430.070-0.0350.5960.035
SR1.396-0.6440.9680.0350.514-0.6220.246-0.1800.775-0.120
ResultNNNNNNNNNN
D20A006MK0.4800.287-0.069-0.071-0.2060.972-0.7541.038-0.6301.196
SR0.2550.349-0.162-0.250-0.1860.971-0.7420.864-0.4491.031
ResultNNNNNNNNNN
D20A009MK-0.5780.973-0.6040.046-0.5690.805-0.9911.519-1.5841.656
SR-0.5410.978-0.2050.123-0.3960.702-0.9651.502-1.5751.706
ResultNNNNNNNNNI
D20A011MK-0.030-0.076-0.0300.406-0.4450.569-0.7420.313-0.3880.165
SR-0.180-0.109-0.0200.422-0.4420.577-0.8830.449-0.5600.158
ResultNNNNNNNNNN
D20A013MK-0.764-0.009-0.0980.451-0.0180.537-0.4630.789-0.6020.976
SR-0.680-0.083-0.0900.4460.0170.580-0.5690.881-0.5910.897
ResultNNNNNNNNNN
D20A014MK1.026-1.9011.256-1.063-0.3450.453-1.0571.381-0.7441.157
SR0.949-1.9101.181-1.058-0.4630.520-1.0641.204-0.7691.084
ResultNDNNNNNNNN
D20A015MK-0.350-0.3030.105-0.1420.245-0.6470.478-0.844-0.1170.459
SR-0.304-0.316-0.013-0.0990.305-0.6840.471-0.787-0.1780.407
ResultNNNNNNNNNN
D20A016MK-1.4760.224-1.6941.583-2.7622.421-2.4072.002-1.9471.783
SR-1.4240.225-1.5991.550-2.5862.301-2.3351.917-1.6961.514
ResultNNNNDIDIDN
D20A017MK-0.819-0.0270.050-1.236-0.2730.505-0.8930.526-0.821-1.706
SR-0.776-0.1380.123-1.205-0.2250.521-0.9130.481-0.8440.785
ResultNNNNNNNNNN
D20A018MK-1.3000.415-0.1510.093-0.091-0.092-0.1510.731-0.9671.798
SR-1.3770.418-0.2790.340-0.0870.029-0.2960.921-1.0781.714
ResultNNNNNNNNNI

N: No Trend, D: Decreasing, I: Increasing, MK: Mann-Kendall, SR: Spearman Rho

4.3. Trend Magnitude

Sen’s slope estimator and linear regression analysis were employed to determine the slope of trends in the AMDS and AMDD time series. The signs of the slopes were in line with the results of the Spearman Rho and Mann-Kendall tests. The slopes of the trend lines are indicated in bold in the time series where the trend was detected (Table 6). The highest slopes of the trend line, ranging from 0.6 to 0.8 for both tests, were detected in the AMDS time series of station 7767 for the SPI 9-month time scale. The other evaluated slopes of the trend were 0.067 in AMDD time series of station 8275 for the SPI 3-month time scale, 0.179 in the AMDD time series of station 8275 for the SPI 6-month period, 0.106 in the AMDD time series of station 8275 for the SPI 9-month period, 0.083 in AMDS time series of station 8275 for the SPI 12-month period, and 0.138 in the AMDD time series of station 8275 for the SPI 12-month period.

Table 6

Trend Magnitude Test results for all stations.

StationTestsSPI 1SPI 3SPI 6SPI 9SPI 12
AMDSAMDDAMDSAMDDAMDSAMDDAMDSAMDDAMDSAMDD
7767SS0.1320.0000.340-0.1250.495-0.3330.768-0.6000.710-0.500
R0.150-0.0790.297-0.1180.480-0.3790.759-0.4640.798-0.432
8275SS-0.0070.000-0.0320.067-0.0910.179-0.0760.106-0.0710.083
R-0.0180.024-0.0450.073-0.1290.138-0.1490.145-0.1890.138
17255SS0.0020.0000.0060.0000.0100.0000.0050.0000.0000.000
R-0.003-0.0030.0000.003-0.0140.001-0.010-0.0200.002-0.027
17355SS-0.0250.000-0.0500.050-0.0250.000-0.0180.000-0.0070.000
R0.001-0.003-0.0120.057-0.0330.000-0.0380.019-0.0520.001
17866SS0.0000.000-0.0130.000-0.0160.042-0.0280.032-0.0230.036
R-0.0020.017-0.0210.021-0.0250.059-0.0330.034-0.0500.060
17868SS0.0240.0000.0320.0000.0160.000-0.0020.0000.0000.000
R0.018-0.0210.032-0.0120.068-0.0070.0610.0050.079-0.037
17870SS0.0050.0000.0130.0000.0180.000-0.0050.010-0.0090.000
R0.019-0.0100.021-0.0120.0240.0120.0050.023-0.0080.023
17871SS0.0810.0000.0330.0000.0170.100-0.0810.125-0.0730.143
R0.042-0.0040.0800.042-0.0070.096-0.0780.109-0.1240.175
17908SS0.0060.0000.0300.0000.0310.000-0.0010.000-0.0010.000
R0.005-0.0060.031-0.0230.0260.013-0.0030.020-0.0180.026
17960SS-0.0060.0000.0140.0000.0280.0000.0190.0000.010-0.048
R-0.0150.0040.0040.0010.029-0.0220.061-0.0450.080-0.072
17979SS-0.0100.000-0.0010.0000.0290.0000.022-0.0280.0190.000
R-0.0110.0030.003-0.0080.043-0.0260.062-0.0370.098-0.046
D20A001SS0.0080.0000.0010.0000.0080.0000.0000.0000.0000.000
R-0.010-0.001-0.0410.022-0.0090.002-0.005-0.017-0.0210.149
D20A002SS-0.0380.000-0.0110.000-0.0230.000-0.003-0.0770.0000.000
R-0.0250.057-0.029-0.009-0.1160.021-0.102-0.049-0.123-0.014
D20A004SS0.0620.0000.1060.0000.1080.0000.0160.0000.0380.000
R0.064-0.0470.144-0.0720.175-0.1110.190-0.0210.1370.044
D20A006SS0.0800.000-0.0190.000-0.0620.367-0.2560.438-0.3450.500
R0.0120.042-0.0050.038-0.0830.245-0.1010.329-0.1340.378
D20A009SS-0.0050.000-0.0020.000-0.0210.000-0.0230.050-0.0250.037
R-0.0100.016-0.006-0.002-0.0030.016-0.0170.060-0.0200.066
D20A011SS-0.0010.000-0.0010.000-0.0270.000-0.0290.0000.0000.000
R-0.0140.006-0.0080.008-0.0130.017-0.0330.024-0.0970.028
D20A013SS-0.0110.000-0.0050.0000.0000.000-0.0060.0000.0000.000
R-0.010-0.0030.0030.0110.0210.0100.0200.0160.0120.032
D20A014SS0.0140.0000.0420.000-0.0090.000-0.0240.054-0.0120.000
R0.024-0.0420.051-0.0350.0170.018-0.0100.051-0.0270.057
D20A015SS-0.0060.0000.0050.0000.0070.0000.017-0.0340.0000.000
R0.007-0.005-0.011-0.006-0.001-0.0250.016-0.038-0.0010.013
D20A016SS-0.0330.000-0.0720.043-0.1560.133-0.0960.100-0.0710.091
R-0.0340.011-0.0810.050-0.1630.125-0.1840.111-0.1980.115
D20A017SS-0.0340.0000.002-0.080-0.0340.000-0.0420.024-0.0340.079
R-0.008-0.0160.034-0.084-0.0230.020-0.0730.047-0.0590.065
D20A018SS-0.0500.000-0.0260.000-0.0510.000-0.0180.100-0.1270.286
R-0.0560.031-0.0110.025-0.019-0.014-0.0980.106-0.2030.227

4.4. Change Point Detection

After determining the presence of a trend as well as its slope in the time series of extreme drought events, SNHT, BRT, and PT were applied to define the point where the change had started. The change point is identified as the point where the change point occurs in SNHT and BRT, which can be depicted graphically. The outcomes of the three methods were compared to assess a precise change point. As two of the three tests presented the same results, that specific point was chosen to be the change point for the particular time series. However, if the outcomes of all three tests were inconsistent, the change point selection was denoted as suspicious (S). The change points for extreme drought events of the AMDS and AMDD time series are presented in Table 7.

Table 7

Change point results of extreme drought events.

StationTestsTest criteriaSPI 1SPI 3SPI 6SPI 9SPI 12
AMDSAMDDAMDSAMDDAMDSAMDDAMDSAMDDAMDSAMDD
7767BRTR/sqrt(n)min0.911.140.990.831.481.091.291.171.341.21
CPNCNC2000NC20012008200120082001NC
SNHTT3.154.976.083.807.086.057.625.687.254.94
CPNCNC2000NC20012008200120082001NC
PTU28315928626158656834
CPNCNC2000NC20012004200120082008NC
ResultNCNC2000NC20012008200120082001NC
8275BRTR/sqrt(n)min0.990.931.071.491.21.221.061.51.051.4
CPNCNCNC1979NC1980NC196919801980
SNHTT3.923.12.316.914.46.73.556.36.26.4
CPNCNCNC1979NC1969NC196919821969
PTU6565631235910592102113105
CPNCNCNC1979NC1979NC196919801969
ResultNCNCNC1979NCSCNC196919801969
17866BRTR/sqrt(n)min0.90.930.840.90.941.120.960.951.041.07
CPNCNCNCNCNC1984NCNCNCNC
SNHTT1.33.112.092.382.177.12.172.222.843.49
CPNCNCNCNCNC1981NCNCNCNC
PTU102127122135148225185157228195
CPNCNCNCNCNC1981NCNCNCNC
ResultNCNCNCNCNC1981NCNCNCNC
D20A009BRTR/sqrt(n)min0.831.660.790.820.80.870.820.90.831.14
CPNCNCNCNCNCNCNCNCNC1982
SNHTT3.955.771.762.822.651.892.383.882.017.3
CPNCNCNCNCNCNCNCNCNC1982
PTU143130112121115126148130123205
CPNCNCNCNCNCNCNCNCNC1982
ResultNCNCNCNCNCNCNCNCNC1982
D20A014BRTR/sqrt(n)min0.951.141.0680.961.051.531.141.871.231.64
CPNC1980NCNCNCNCNCNCNCNC
SNHTT2.783.43.425.532.173.892.75.53.175.32
CPNC1975NCNCNCNCNCNCNCNC
PTU128159134134136157141128116142
CPNC1999NCNCNCNCNCNCNCNC
ResultNCSCNCNCNCNCNCNCNCNC
D20A016BRTR / sqrt(n)min0.820.971.040.881.221.261.31.651.340.9
CPNCNCNCNC19881988198819881996NC
SNHTT3.955.435.163.487.027.286.937.76.955.3
CPNCNCNCNC19881984198819841988NC
PTU11669144120188181193188169158
CPNCNCNCNC19881988198819881988NC
ResultNCNCNCNC19881988198819881988NC
D20A018BRTR/sqrt(n)min0.770.850.940.820.610.670.830.841.011.04
CPNCNCNCNCNCNCNCNCNC1988
SNHTT3.421.821.682.911.421.641.621.942.996.98
CPNCNCNCNCNCNCNCNCNC1984
PTU45343631222328424463
CPNCNCNCNCNCNCNCNCNC1988
ResultNCNCNCNCNCNCNCNCNC1988

NC: No change

All three tests were performed with a significance level of 90%. According to test results, no abrupt changes were found at station 7767 for the SPI 1-month period in the AMDS and AMDD time series, the SPI 3-month period in the AMDD time series, and the SPI 12-month period in the AMDD time series. Changes in the AMDS time series were noted in the year 2000 for SPI 3-, in 2001 for SPI 6-, 9-, and 12-month time scales. Changes in the AMDD series were found in 2008 for SPI 6- and SPI 9- at station 7767. We detected no changes in both the AMDS and AMDD time series at station D20A014. Changes were observed in 1981 at station 17866 for the SPI 6- in the AMDD time series, in 1982 at station D20A009 for the SPI 12- in the AMDD time series, and in 1988 at station D20A018 for the SPI 12-AMDD time series (Fig. 3).

In general, numerous change points were detected at stations 7767, 8275, and D20A016, with a significance level of 90%. Although there were no changes for the SPI 1-month time scales and the SPI 3-month time scales at station 8275, the years 1979, 1969, and 1980 were critical for the SPI 6-AMDS time series, the SPI 9-AMDD time series, and the SPI 12-AMDS time series, respectively. At station D20A016, the change point was detected in 1988, which was a critical year for the SPI 6-AMDS and AMDD time series, for the SPI 9-AMDS and AMDD time series, and for the SPI 12-AMDS time series.

5. Discussion and Conclusions

The AMDS and the AMDD trends in extreme drought events were evaluated in this study. The SPI drought index was used to assess the extreme events for 23 meteorological stations in Ceyhan Basin over several time scales of 1, 3, 6, 9, and 12 months. The AMDS and AMDD time series are important extreme hydrological drought events in disaster preparedness and the regional management of water supplies. The homogeneity of the extreme drought occurrence time series was investigated using the Wallis-Moore and Wald-Wolfowitz methods, and Mann-Kendal and Spearman Rho trend analyses were carried out. Sen’s slope estimator and linear regression analysis were used to assess the magnitude of trends, and the SNHT, BRT, and PT were used to identify the change point that denotes the start of the trend.

The number of stations with non-homogeneous data increased as the SPI time scales increased, especially in the AMDS time series. Most of the stations seemed to have homogeneous time series; however, the data in a few stations were homogeneous by one of the methods and non-homogeneous by the other one, with a significance level of 90%.

In the AMDS series, an increasing trend was observed only in the station 7767 data for 3-, 6-, 9-, and 12-month time scales, whereas decreasing trends were detected in station 8275 for the 12-month time scale and in station D20A016 for the 6-, 9-, and 12-month time scales. Nevertheless, in the AMDD time series, increasing trends were noted in station 8275 for the 3-, 6-, 9-,and 12-month time scales, station 17866 for the 6-month time scale, station D20A009 for the 12-month time scale, station D20A016 for the 6- and 9-month time scales, and station D20A018 for the 12-month time scale, whereas decreasing trends were perceived in the station 7767 for the 6- and 9-month time scales and in station D20A014 for the 1-month time scale.

Fig. 3.

Change point results for different time scales of SPI 3, 6, 9, and 12 months.

https://www.mhwm.pl/f/fulltexts/154573/MHWM-9-0011-g003_min.jpg

These results indicate that the AMDD time series with higher time scales revealed an increasing trend, whereas the AMDS time series showed a decreasing trend, which was statistically significant at the 90% confidence level.

Possible change points started after 1988 for D20A016 and D20A018, whereas the data in stations 7767, 8275, 17866, and D20A009 shifted after late 2000, 1969, 1981, and 1982, respectively. Change point detection tests captured different years in the data of some stations, most likely because of climate type, topography, and hydrological conditions, among others. Investigations on drought and climate change at local scales need to be periodically repeated to assess the impacts of precipitation, temperature, evapotranspiration, and groundwater level on agriculture, water resources and food security, energy production and security, economic growth, ecosystem sustainability, environmental health, and health security. The results of this study are consistent with previous findings. For instance, Di̇ki̇ci̇ and Aksel, (2021) indicated that precipitation is generally decreasing, considering a 50-year period. According to their calculations, the basin faces a severe drought threat in the 20-year return period. Tosunoglu and Kisi (2016) investigated the trends of maximum hydrologic drought variables, annual maximum duration (AMD), and annual maximum severity (AMS) using MK, modified MK, and the recently developed Sen’s innovative trend analysis (ITA) techniques. Increasing trends in the severity of the drought at the stations 2316 and 2323 were captured by the modified MK and ITA techniques, respectively. The possibility of future droughts or water shortages in these stations should be considered. The results of this paper provide valuable information to water resource management decision-makers in the Ceyhan River Basin for evaluating the effects of droughts and preparing for drought mitigation measures to avoid future drought risks.