Abstract
The difficulty of monitoring hydrometeorological trends remains relevant because of the importance of climate change. In this study, innovative polygon trend analysis (IPTA), new graphical methods, have recently been presented, Mann-Kendall and innovative trend analysis (ITA) with significance tests are employed at monthly and annual hydrometeorological data for Ankara province in Turkey. Finally, sequential Mann-Kendall (SQ-MK), CUSUM, and standard normal homogeneity test (SNHT) are applied to detect any abrupt changes in annual time series. The results indicate that ITA with significance test, IPTA, and MK capture precipitation trends in 83, 94, and 0.2% of all months (3 stations × 12 months), temperature trends in 86, 75, and 22% of all months, relative humidity trends in 80, 80, and 30% of all months, and evapotranspiration series in 91, 80, and 47% of all months. These findings suggest that the ITA with significance test and IPTA are more sensitive than the MK test, with the precipitation series being the most sensitive and the evapotranspiration series being the least sensitive. SQ-MK, CUSUM, and SNHT tests on station 17664 were successful in detecting a change in annual evapotranspiration in 2005. The annual total evaporation series have been rising since 2005, according to the SQ-MK test.
HIGHLIGHTS
Graphical and traditional trend tests were utilized at monthly and annual hydrometeorological data in the context of climate change.
Graphical tests are most sensitive than the traditional test with the precipitation series being the most sensitive.
Evaporation in February, July, and September indicate an increasing trend, and maximum temperature in February and July observed an increasing trend.
Graphical Abstract
INTRODUCTION
Climate change, rapid urbanization and industrialization, land degradation, and growing populations are all factors that have a substantial impact on numerous hydrometeorological variables. Variability in rainfall, relative humidity, evapotranspiration, and temperature can have a negative impact on the occurrences of natural hazards, particularly floods and droughts (Lehner et al. 2018; Arab Amiri & Gocić 2021; Esit et al. 2021). By 2100, climate change is expected to increase a 0.5–1.0 m rises in mean sea level and a 1.0–5.7 °C increase in mean temperature as well as higher rainfall, wind speed, and more frequent extreme events in some areas, according to the Intergovernmental Panel on Climate Change (Allan 2021). Climate change has a variety of substantial effects in humid, semi-arid, and dry regions, according to IPCC studies (IPCC 2013). Climate change influences the hydrological cycle, water availability, quality, and services, and observations and records from all continents and oceans indicate unexpected behaviors. These unexpected hydrological cycle behaviors may be seen as either trends or sudden increases in average and variation values. They also have an impact on a variety of human activities, including groundwater recharge, water supply, agricultural activities, hydroelectric power generation, and irrigation methods (Şen 2012, 2014; Şan et al. 2021). As a result, interested institutions and organizations should investigate changes in hydrometeorological information to obtain a preliminary understanding of climate change.
One of the most useful methods for observing the effects of climate change on hydrometeorological variables is trend analysis. Most of the hydrometeorological variables are stochastic and complex (Hırca et al. 2022). While monitoring their trends can provide a useful perspective on water resources and meteorological sciences, the intrinsic characteristics of hydroclimatic variables require new technologies and methods for trend assessment (Şen 2020). The use of non-parametric tests in statistical analysis is crucial. They have been commonly employed for the analysis of hydrometeorological time series for the past two decades. Mann-Kendall (MK) test, Spearman's Rho (SR) test, regression analysis, and wavelet analysis are all frequently suggested as classic mathematical methods (Malik et al. 2020; Ivanovski et al. 2021; Yuce & Esit 2021; Shokrian et al. 2022).
The approaches for trend analysis described above, on the other hand, are pure statistical methods that do not allow for the detection of trends in low, medium, and high values in a single computation procedure. The majority of non-parametric tests can be worthless in some cases without graphical and exploratory data analyses (Cengiz et al. 2020). Onyutha (2016) proved that using a combination of graphical and statistical trend testing methods provides more useful and influential information than using only statistical methods. The majority of classical trend analysis methods show monotone trends across the whole time series (Yue et al. 2002; Shadmani et al. 2012; Eris et al. 2019; Isensee et al. 2022). However, it does not demonstrate any periodic variations (such as monthly). Hence, trend characteristics are critical when illustrating seasonal trend behavior. Because of anthropogenic effects in the atmospheric and hydrologic environments, as well as the resulting climate change, the stationarity principle is no longer applicable (Mohorji et al. 2017).
Hydrometeorological parameters (precipitation, temperature, relative humidity, and evapotranspiration) variation have a direct effect on agriculture, irrigation, and water resources system (Abdi et al. 2016; Cetin et al. 2018; Yuce et al. 2019). In particular, the change in precipitation is important in the hydrological cycle, as it has effects on both the occurrence of drought and floods. To overcome the above restrictions, innovative trend analysis (ITA) introduced by Şen (2012), was recently used to detect trends in hydrological and meteorological data in various parts of the world (Mallick et al. 2021; Singh et al. 2021; Akçay et al. 2022; Gumus et al. 2022). Although studies have demonstrated that sample size, distribution shape, and serial correlation influence the ITA approach, its benefits can be implemented in a variety of aspects (Serinaldi et al. 2020). In comparison to the MK and SR tests, Kisi (2015) discovered that the graphs created by the ITA approach can better demonstrate hidden trends in pan evaporation. In comparison to other statistical tests, the strength of ITA is its high performance of graphical visualization in identifying hidden trends in hydrometeorological time series (Cengiz et al. 2020). The ITA method is not significant if the data do not fit a normal distribution, the data length is small, or the data are serially dependent (Şen 2012, 2017). The method has been widely utilized to detect trends (Yilmaz 2019; Caloiero 2020). A novel method also proposed by Şen et al. (2019), the Innovative Polygonal Trend Analysis (IPTA), which can identify small-scale time series (daily, weekly, monthly, etc.) as well as determine trend transitions between successive sections of two equal segments of the original hydrometeorological time series, resulting in the polygon trend. A few studies were conducted to determine the applicability of IPTA in practice (Ceribasi & Ceyhunlu 2020; Achite et al. 2021; Harka et al. 2021; Ahmed et al. 2022). However, more research is needed to see if IPTA can be used to analyze hydrometeorological data and if it may be a suitable alternative trend analysis approach, such as the MK test.
The degree and slope of trend transitions between consecutive segments (e.g., months) can be evaluated using the IPTA method, and information on the magnitude and slope of the trend transitions between consecutive segments (e.g., months) can be obtained. This approach can also be used to observe the one-year behavior of time series. Internal variability that is time-dependent might be assessed in this way (Akçay et al. 2022). This approach was applied to monthly precipitation and flow data of selected stations in various parts of the world by Şen et al. (2019). The findings revealed many increasing/decreasing trends, and the slopes and lengths of these trends were determined. Furthermore, the monthly precipitation is studied using the IPTA graphical approach, which provides a clear and intuitive knowledge of precipitation fluctuations throughout the year, allowing for active and efficient response to soil erosion and vegetation degradation. Yenice & Yaqub (2022) evaluated trend analysis of temperature data by applying the IPTA method at six stations in Turkey. Şan et al. (2021) compared ITA, IPTA, and MK trend test for rainfall data in Vietnam. Their results showed that trends could be seen in almost 90% of all months in IPTA and ITA, while with the MK test, this rate was only 23%. They also noted that ITA and IPTA are more sensitive than MK in detecting a trend. Ahmed et al. (2022) observed the trend analysis of monthly precipitation in the Hindukush–Karakoram–Himalaya River basins of Pakistan using MK, Sen's slope, ITA, and IPTA methods. According to their findings, the MK method was less sensitive than the IPTA and ITA methods at determining trends. Even if 90% of the innovative methods are compatible with one another, the IPTA provides more details on trend transitions between subsequent time series segments. Wu et al. (2022) analyzed the precipitation variations in the Tai Lake Basin from 1971 to 2018 based on innovative trend approaches. Their results indicated that the detection results from the ITA method and traditional trend analysis methods (linear regression analysis, MK, and modified MK) were consistent. As an addition to ITA, IPTA can systematically identify consecutive seasons and monthly transition characteristics. The period and transition properties of a time series cannot be seen using the general trend detection method. The innovative polygon trend analysis (IPTA), which is based on ITA, can reflect the periodic properties of hydrometeorological variables and is frequently used for monthly and seasonal scales as a way to discover the transition trends between months and seasons (Ceribasi & Ceyhunlu 2020; Adelodun et al. 2022; Akçay et al. 2022; Wu et al. 2022).
Since ITA and IPTA methods are newly published, a few studies have been done. The previous studies generally focused on precipitation and temperature data. But other meteorological parameters that are important in the hydrology cycle need to be examined. Therefore, all parameters should be evaluated together. Unlike other studies, evapotranspiration and relative humidity data were examined for the first time. Also, no research employing innovative techniques on temperature or precipitation data has yet been done for the Ankara province. Therefore, it provides the opportunity to look at the region to be analyzed from a broad perspective. The aim of this study is to detect the trend of hydroclimate data (precipitation, temperature, relative humidity, and evapotranspiration) monthly and annually by comparing it with both the classical MK method and the ITA and IPTA methods. Ankara region was chosen as the study area because it is one of the most important provinces as the capital city of Turkey. To identify the trend slope, Sen's slope was used in this study. CUSUM and SQ-MK tests were employed to detect abrupt changes in hydrometeorological variables.
STUDY AREA AND DATA
Ankara is located in the northern part of the Central Anatolian Region, which constitutes the interior of Turkey. It is neighbors with Çankırı in the northeast, Kırıkkale in the east, Kırşehir and Aksaray in the southeast, Konya in the south, Eskişehir in the west, and Bolu in the northwest. In addition, the city of Ankara is located between 38°40′ and 40°45′ north latitudes and 30°50′ and 33°53′ east longitudes. There are climatic differences from place to place in the wide area of the province. In the south, there is the steppe climate, which is the distinctive feature of the Central Anatolian climate, while in the north, the temperate and rainy states of the Black Sea climate can be seen. In this region where the continental climate prevails, winter temperatures are low and summer is hot. The hottest month is July–August, and the coldest month is January. The amount of precipitation falling in the region differs in the northern and southern parts. Kızılcahamam and Çubuk districts located in the northern direction are characterized by the precipitation regime of the Black Sea Region. In the south direction, it shows the climate characteristic of the Central Anatolia Region (Çi̇çek 2003).
The average temperature in the province is 11.7 °C, and the annual average precipitation is 389.1 mm. The highest temperature value was determined as 40.8 °C and the lowest temperature as −24.9 °C. The number of days with frost is 60–117, and the number of snowy days is 30.5 days a year. The highest snow thickness was determined as 30 cm. When we look at the wind conditions of the city center and stations in general. It is seen that the prevailing wind changes depending on the terrain structure. According to this, prevailing winds are northeast in Ankara (center), Esenboğa, Çubuk, Ayaş, and Yenimahalle districts, west in Haymana (İkizce), Sincan, Dikmen district and Nallıhan districts, north in Polatlı and Şerefikoçhisar districts, southwest in Etimesgut and Elmadağ districts, southeast in Kızılcahamam district, and northeast in Beypazarı district work in the direction. The months with strong winds are March and April. The highest wind speed detected in Ankara is 29.2 m/s. The only exception to these measurements is the tornado disaster that caused great damage in the Akyurt district in 2007. According to the values measured for many years, Ankara's average pressure value is 913.1 Mb, the highest-pressure value detected is 935.0 Mb, and the lowest pressure value is 891.0 Mb (Afshar et al. 2020; Danandeh Mehr et al. 2020).
Station code . | Meteorological parameters . | Latitude . | Longitude . | Earliest record year . | Latest record year . | Mean . | SD . | Cv . | Cs . | r1 . |
---|---|---|---|---|---|---|---|---|---|---|
17128 | ATP | 40.124 | 32.999 | 1956 | 2021 | 410.503 | 84.345 | 0.205 | −0.032 | 0.244 |
AMT | 24.091 | 0.936 | 0.039 | −0.061 | 0.315 | |||||
AMRH | 65.206 | 3.359 | 0.052 | −0.317 | 0.584 | |||||
ATE | 649.695 | 51.305 | 0.079 | 0.453 | 0.417 | |||||
17664 | ATP | 40.472 | 32.644 | 1959 | 2021 | 579.602 | 113.560 | 0.196 | 0.350 | 0.012 |
AMT | 23.804 | 1.031 | 0.043 | −0.123 | 0.416 | |||||
AMRH | 66.366 | 3.121 | 0.047 | 0.091 | 0.422 | |||||
ATE | 635.094 | 43.885 | 0.069 | 0.515 | 0.367 | |||||
17130 | ATP | 39.972 | 32.863 | 1926 | 2021 | 391.051 | 79.802 | 0.204 | 0.373 | 0.214 |
AMT | 25.054 | 0.995 | 0.040 | −0.562 | 0.148 | |||||
AMRH | 60.112 | 2.880 | 0.048 | −0.385 | 0.483 | |||||
ATE | 712.428 | 52.285 | 0.073 | 0.653 | 0.400 |
Station code . | Meteorological parameters . | Latitude . | Longitude . | Earliest record year . | Latest record year . | Mean . | SD . | Cv . | Cs . | r1 . |
---|---|---|---|---|---|---|---|---|---|---|
17128 | ATP | 40.124 | 32.999 | 1956 | 2021 | 410.503 | 84.345 | 0.205 | −0.032 | 0.244 |
AMT | 24.091 | 0.936 | 0.039 | −0.061 | 0.315 | |||||
AMRH | 65.206 | 3.359 | 0.052 | −0.317 | 0.584 | |||||
ATE | 649.695 | 51.305 | 0.079 | 0.453 | 0.417 | |||||
17664 | ATP | 40.472 | 32.644 | 1959 | 2021 | 579.602 | 113.560 | 0.196 | 0.350 | 0.012 |
AMT | 23.804 | 1.031 | 0.043 | −0.123 | 0.416 | |||||
AMRH | 66.366 | 3.121 | 0.047 | 0.091 | 0.422 | |||||
ATE | 635.094 | 43.885 | 0.069 | 0.515 | 0.367 | |||||
17130 | ATP | 39.972 | 32.863 | 1926 | 2021 | 391.051 | 79.802 | 0.204 | 0.373 | 0.214 |
AMT | 25.054 | 0.995 | 0.040 | −0.562 | 0.148 | |||||
AMRH | 60.112 | 2.880 | 0.048 | −0.385 | 0.483 | |||||
ATE | 712.428 | 52.285 | 0.073 | 0.653 | 0.400 |
ATP, annual total precipitation (mm); AMT, annual maximum temperature (°C); AMRH, annual mean relative humidity (%); ATE, annual total evapotranspiration (mm).
METHODS
Mann-Kendall (MK) test
A positive Z number shows an upward trend, while a negative value shows a downward trend. Critical test statistical values at 90, 95, and 99% probability levels are 1.645, 1.96, and 2.58, respectively (Yuce & Esit 2021).
ITA with significance test
The ITA method developed by Şen (2012) is a technical analysis method that shows a possible general increase or decrease in a given time series. In this method, the available data are arranged according to the order, and then they are divided into two equal series. These series, which are divided into two, are ordered from smallest to largest. The first part of the series (Xi) is placed on the X-axis of the Cartesian coordinate system, and the second part (Xj) is placed on the Y-axis. The data are first examined for its position on the 1:1 line. Data above the 1:1 line indicate an increasing trend in the data, while data below the 1:1 line indicate a decreasing trend in the data. Secondly, it is examined in which cluster the data are among the low, medium, or high clusters. Finally, if the data are clustered around the 1:1 line, it is concluded that there is no trend (Caloiero et al. 2018; Serinaldi et al. 2020).
The value z produced from the standard normal distribution at a given confidence level is referred to as . It is considered as an increasing (decreasing) trend if the trend slope exceeds the upper (lower) confidence level. There is no statistically significant trend at a specific confidence level if these conditions are not met. The level of confidence in this study is 95%.
IPTA method
Sen's slope estimator
Standard normal homogeneity test (SNHT), CUSUM test, and sequential Mann-Kendall test (SQ-MK)
Here, and are mean and SD, respectively. The year in which reaches the maximum value is considered as the point of change.
Distribution free CUSUM (trendchange package in R) test is employed to determine the abrupt change point in a series of data. The change point in the time series will be the point at which the cumulative sum reaches its maximum value. The significance of the change point is indicated if the highest value is equal to or greater than the critical value ((Patakamuri & Das 2019).
It is well known that using the MK method to detect trends at the end of any time period does not provide a complete trend picture (structure of the trend) for the entire series. There may be changes in the trend across the investigation period, which can be noticed by performing the test for each period separately (Sneyers 1990). According to Makokha & Shisanya (2010), negative or positive trends are not always significant. However, it can be discovered using SQ-MK graphs.
The backward sequential statistic u () are calculated in the same way, but starting from the end of the series. Many scientists have used this approach to determine the starting point of trends (Rahman et al. 2017; Salehi et al. 2020; Alhathloul et al. 2021) rather than identifying the whole trend. In this study, the SQ-MK method is applied to detect an abrupt change in the time series of hydrometeorological data.
RESULTS AND DISCUSSION
Correlation analysis between meteorological variables
MK test results
Table 2 shows the results of the MK test for annual hydrometeorological variables (precipitation, temperature, relative humidity, and evapotranspiration). Confidence levels are selected to be 10, 5, and 1%. Annual total precipitation has no significant trend at stations 17128 and 17664, while an increasing trend has been found at station 17128. Stations 17128 and 17130 indicate a significant increasing trend for annual maximum Temperature, whereas there is no significant trend at station 17664. Additionally, all stations indicate no increasing/decreasing trend at annual mean relative humidity. But all stations show an increasing trend (10 and 5% confidence level) at annual total evapotranspiration.
MK and Sen's slope tests results for monthly hydroclimate variables are presented in Table 3. As shown in Table 3, the months in bold demonstrate significant trends at the 90, 95, and 99% confidence levels. According to monthly precipitation data, an increasing trend is indicated for station 17130 in January (90% confidence level), March (95% confidence level), and October (90% confidence level), while a decreasing trend is detected for station 17128 in February (90% confidence level). An increasing trend is obtained from stations 17128 and 17130 in January, February, June, July, and September for temperature data, whereas no trends are observed at station 17664.
Station No . | Meteorological Variables . | Tests . | Jan . | Feb . | Mar . | Apr . | May . | Jun . | Jul . | Aug . | Sep . | Oct . | Nov . | Dec . |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
17128 | ATP | MK | 0.094 | − 1.826* | 0.465 | −0.149 | −1.107 | 1.218 | 0.011 | 1.484 | −1.218 | 1.533 | 1.173 | −1.256 |
SS | 0.025 | −0.263 | 0.066 | −0.029 | −0.220 | 0.196 | 0.000 | 0.100 | −0.076 | 0.176 | 0.140 | −0.195 | ||
AMT | MK | 1.82* | 2.16** | 1.617 | 0.255 | 1.932* | 2.375** | 2.226** | 1.639 | 2.857*** | 1.035 | 0.736 | −0.537 | |
SS | 0.041 | 0.060 | 0.025 | 0.005 | 0.027 | 0.033 | 0.029 | 0.020 | 0.035 | 0.018 | 0.011 | −0.007 | ||
AMRH | MK | 0.083 | − 2.640 | − 2.873** | −0.974 | −1.472 | 1.273 | 0.482 | −0.487 | −1.893 | 0.315 | −0.310 | −0.011 | |
SS | 0.000 | −0.092 | −0.096 | −0.031 | −0.056 | 0.050 | 0.012 | −0.023 | −0.100 | 0.019 | −0.012 | 0.000 | ||
ATE | MK | 1.370 | 2.34** | 1.938** | 0.554 | 2.210** | 2.155** | 3.345*** | 3.244*** | 3.877*** | 2.696*** | 0.709 | 0.124 | |
SS | 0.000 | 0.015 | 0.104 | 0.033 | 0.125 | 0.143 | 0.230 | 0.261 | 0.216 | 0.157 | 0.024 | 0.000 | ||
17130 | ATP | MK | 1.64* | 0.332 | 2.415** | 1.317 | 0.237 | 1.456 | 0.256 | 0.329 | −0.196 | 1.880* | 0.282 | −0.101 |
SS | 0.144 | 0.033 | 0.189 | 0.119 | 0.027 | 0.127 | 0.009 | 0.006 | −0.007 | 0.125 | 0.018 | −0.010 | ||
AMT | MK | 2.43** | 3.036*** | 0.779 | 0.003 | 0.855 | 3.119*** | 2.014** | 1.257 | 1.779* | 1.605 | −1.150 | 0.516 | |
SS | 0.029 | 0.039 | 0.009 | 0.000 | 0.007 | 0.023 | 0.015 | 0.009 | 0.015 | 0.014 | −0.008 | 0.005 | ||
AMRH | MK | − 2.504** | − 4.110*** | − 2.384** | −0.108 | 0.247 | 2.219** | 2.532** | 3.44*** | 1.599 | 2.206* | −1.019 | − 2.289** | |
SS | −0.031 | −0.080 | −0.050 | −0.003 | 0.006 | 0.055 | 0.053 | 0.083 | 0.046 | 0.067 | −0.024 | −0.030 | ||
ATE | MK | 2.16** | 2.51** | 2 93*** | 0.874 | 0.618 | 0.805 | 2.265** | 1.568 | 2.29** | 1.600 | −0.931 | 0.279 | |
SS | 0.000 | 0.028 | 0.110 | 0.036 | 0.023 | 0.031 | 0.102 | 0.072 | 0.090 | 0.061 | −0.022 | 0.000 | ||
17664 | ATP | MK | 0.415 | −1.257 | 0.973 | 0.504 | −0.896 | 0.243 | −0.249 | −0.474 | −0.581 | 0.712 | −1.163 | −1.341 |
SS | 0.157 | −0.367 | 0.163 | 0.112 | −0.204 | 0.065 | −0.032 | −0.043 | −0.070 | 0.144 | −0.191 | −0.475 | ||
AMT | MK | −0.457 | 0.932 | −1.323 | −1.033 | 0.237 | 0.160 | 0.861 | 0.463 | −0.481 | −1.478 | −0.487 | −0.154 | |
SS | −0.009 | 0.025 | −0.022 | −0.024 | 0.003 | 0.000 | 0.013 | 0.005 | −0.007 | −0.029 | −0.009 | −0.003 | ||
AMRH | MK | 3.30*** | −0.285 | −0.570 | −1.459 | −1.263 | 0.540 | − 2.047** | − 1.934* | −0.807 | 1.436 | 0.042 | 1.317 | |
SS | 0.130 | −0.015 | −0.019 | −0.060 | −0.074 | 0.029 | −0.083 | −0.089 | −0.035 | 0.067 | 0.000 | 0.050 | ||
ATE | MK | 0.117 | 2.07** | 0.777 | −0.154 | 1.550 | 2.26** | 3.35*** | 4.09*** | 3.01*** | 1.507 | −0.196 | −0.670 | |
SS | 0.000 | 0.000 | 0.050 | 0.000 | 0.100 | 0.150 | 0.230 | 0.336 | 0.172 | 0.083 | −0.008 | 0.000 |
Station No . | Meteorological Variables . | Tests . | Jan . | Feb . | Mar . | Apr . | May . | Jun . | Jul . | Aug . | Sep . | Oct . | Nov . | Dec . |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
17128 | ATP | MK | 0.094 | − 1.826* | 0.465 | −0.149 | −1.107 | 1.218 | 0.011 | 1.484 | −1.218 | 1.533 | 1.173 | −1.256 |
SS | 0.025 | −0.263 | 0.066 | −0.029 | −0.220 | 0.196 | 0.000 | 0.100 | −0.076 | 0.176 | 0.140 | −0.195 | ||
AMT | MK | 1.82* | 2.16** | 1.617 | 0.255 | 1.932* | 2.375** | 2.226** | 1.639 | 2.857*** | 1.035 | 0.736 | −0.537 | |
SS | 0.041 | 0.060 | 0.025 | 0.005 | 0.027 | 0.033 | 0.029 | 0.020 | 0.035 | 0.018 | 0.011 | −0.007 | ||
AMRH | MK | 0.083 | − 2.640 | − 2.873** | −0.974 | −1.472 | 1.273 | 0.482 | −0.487 | −1.893 | 0.315 | −0.310 | −0.011 | |
SS | 0.000 | −0.092 | −0.096 | −0.031 | −0.056 | 0.050 | 0.012 | −0.023 | −0.100 | 0.019 | −0.012 | 0.000 | ||
ATE | MK | 1.370 | 2.34** | 1.938** | 0.554 | 2.210** | 2.155** | 3.345*** | 3.244*** | 3.877*** | 2.696*** | 0.709 | 0.124 | |
SS | 0.000 | 0.015 | 0.104 | 0.033 | 0.125 | 0.143 | 0.230 | 0.261 | 0.216 | 0.157 | 0.024 | 0.000 | ||
17130 | ATP | MK | 1.64* | 0.332 | 2.415** | 1.317 | 0.237 | 1.456 | 0.256 | 0.329 | −0.196 | 1.880* | 0.282 | −0.101 |
SS | 0.144 | 0.033 | 0.189 | 0.119 | 0.027 | 0.127 | 0.009 | 0.006 | −0.007 | 0.125 | 0.018 | −0.010 | ||
AMT | MK | 2.43** | 3.036*** | 0.779 | 0.003 | 0.855 | 3.119*** | 2.014** | 1.257 | 1.779* | 1.605 | −1.150 | 0.516 | |
SS | 0.029 | 0.039 | 0.009 | 0.000 | 0.007 | 0.023 | 0.015 | 0.009 | 0.015 | 0.014 | −0.008 | 0.005 | ||
AMRH | MK | − 2.504** | − 4.110*** | − 2.384** | −0.108 | 0.247 | 2.219** | 2.532** | 3.44*** | 1.599 | 2.206* | −1.019 | − 2.289** | |
SS | −0.031 | −0.080 | −0.050 | −0.003 | 0.006 | 0.055 | 0.053 | 0.083 | 0.046 | 0.067 | −0.024 | −0.030 | ||
ATE | MK | 2.16** | 2.51** | 2 93*** | 0.874 | 0.618 | 0.805 | 2.265** | 1.568 | 2.29** | 1.600 | −0.931 | 0.279 | |
SS | 0.000 | 0.028 | 0.110 | 0.036 | 0.023 | 0.031 | 0.102 | 0.072 | 0.090 | 0.061 | −0.022 | 0.000 | ||
17664 | ATP | MK | 0.415 | −1.257 | 0.973 | 0.504 | −0.896 | 0.243 | −0.249 | −0.474 | −0.581 | 0.712 | −1.163 | −1.341 |
SS | 0.157 | −0.367 | 0.163 | 0.112 | −0.204 | 0.065 | −0.032 | −0.043 | −0.070 | 0.144 | −0.191 | −0.475 | ||
AMT | MK | −0.457 | 0.932 | −1.323 | −1.033 | 0.237 | 0.160 | 0.861 | 0.463 | −0.481 | −1.478 | −0.487 | −0.154 | |
SS | −0.009 | 0.025 | −0.022 | −0.024 | 0.003 | 0.000 | 0.013 | 0.005 | −0.007 | −0.029 | −0.009 | −0.003 | ||
AMRH | MK | 3.30*** | −0.285 | −0.570 | −1.459 | −1.263 | 0.540 | − 2.047** | − 1.934* | −0.807 | 1.436 | 0.042 | 1.317 | |
SS | 0.130 | −0.015 | −0.019 | −0.060 | −0.074 | 0.029 | −0.083 | −0.089 | −0.035 | 0.067 | 0.000 | 0.050 | ||
ATE | MK | 0.117 | 2.07** | 0.777 | −0.154 | 1.550 | 2.26** | 3.35*** | 4.09*** | 3.01*** | 1.507 | −0.196 | −0.670 | |
SS | 0.000 | 0.000 | 0.050 | 0.000 | 0.100 | 0.150 | 0.230 | 0.336 | 0.172 | 0.083 | −0.008 | 0.000 |
*, **, *** represent the significance levels at 90, 95, 99%, respectively. See
There is a decreasing trend in February (α = %1), March (α = %1), and September (α = %10) for station 17128 considering relative humidity data. For station 17130, a decreasing trend in January–March, and December, and an increasing trend are observed in June–August and October. According to station 17664, an increasing trend is detected in January (α = %1), while there is a decreasing trend in July (α = %5) and August (α = %10). For evapotranspiration data, an increasing trend in February–March and May–October periods at station 17128, an increasing trend in the January–March, July, and September periods for station 17130, and an increasing trend are observed in February and June–September periods. No trends are captured in April and November and December (except relative humidity for station 17130) significantly for all hydrometeorological variables by the MK test.
The highest decreasing Sen's slope is indicated in February (−0.263) for station 17128 (precipitation data), whereas the highest increasing SS is observed in August (0.336) for station 17664 (evapotranspiration data).
ITA with significance test results
At station 17128, trend slopes are observed higher than the lower limits at the winter season (December–January–February) for precipitation data. It means that there is a decreasing trend for the winter season. Decreasing trends are captured in May and September, while increasing trends are detected in the summer season (June–July–August) and October and November considering precipitation data. According to temperature and evapotranspiration data, the trend slope is much higher than the upper limits for all months (except December as a decreasing trend) at significance levels. Thus, increasing trends are observed for all seasons. Relative humidity shows non-similarity compared with temperature data. The decreasing trend appears in the spring season (March–April–May), August, September, and November, while the increasing trend is just observed in October. No trend is detected in June, July, and December.
At station 17130, decreasing trends are detected in January, March, April, June, August, October, and November, while decreasing trends are observed in February, September, and December for precipitation data. Trend slopes indicate lower values than upper/lower limits in May and July. For temperature data, increasing trends in the winter (DJF) season, March, June, July, September, and October and decreasing trends are detected just in November. Evaporation in the summer season (JJA) is in the increasing trend with the January–April period and September and October, while November and December are observed a decreasing trend. According to relative humidity data, decreasing trends occur in the winter season, and notable increasing trends are observed in the summer season. There are no trends in May and November.
At station 17664, precipitation in the winter (DJF) season, as well as May and November, indicate a decreasing trend, whereas an increasing trend appears in March, April, July, September, and October. The temperature in the summer season (JJA) and January, February, and November are observed an increasing trend, notable decreasing trends are shown in October and December. Relative humidity in the spring (MAM) and summer (JJA) seasons appears a decreasing trend with September, while December and January are detected in an increasing trend. According to evapotranspiration data, the summer season, as well as February, March, May, September, and October, demonstrate an increasing trend, while a decreasing trend is captured in January and April.
IPTA results
Station No . | Variables . | Tests . | Jan–Feb . | Feb–Mar . | Mar–Apr . | Apr–May . | May–Jun . | Jun–Jul . | Jul–Aug . | Aug–Sep . | Sep–Oct . | Oct–Nov . | Nov–Dec . | Dec–Jan . |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
17128 – Ankara Airport | P (mm) | Trend length | 15.33 | 8.40 | 4.99 | 11.11 | 20.30 | 31.74 | 2.99 | 7.11 | 19.53 | 4.95 | 29.62 | 7.79 |
Trend slope | 0.95 | − 4.41 | 0.50 | 0.16 | 0.23 | 1.19 | − 0.47 | − 0.01 | 2.87 | − 0.04 | 0.56 | 1.74 | ||
T (°C) | Trend length | 4.23 | 9.57 | 6.51 | 4.74 | 5.03 | 4.71 | 0.24 | 4.50 | 7.02 | 10.44 | 8.48 | 4.27 | |
Trend slope | 1.21 | 0.88 | 0.90 | 1.31 | 0.99 | 1.15 | − 3.80 | 0.98 | 1.11 | 0.96 | 1.23 | 0.68 | ||
RH (%) | Trend length | 6.35 | 9.89 | 6.93 | 0.81 | 6.55 | 12.96 | 1.26 | 6.26 | 13.92 | 13.52 | 10.77 | 0.90 | |
Trend slope | 2.07 | 1.09 | 0.59 | − 2.28 | 0.51 | 1.02 | 1.49 | 0.62 | 1.39 | 0.86 | 1.11 | 4.79 | ||
ET (mm) | Trend length | 3.61 | 18.27 | 38.84 | 48.23 | 36.25 | 40.87 | 11.19 | 60.86 | 50.77 | 43.20 | 16.66 | 3.39 | |
Trend slope | 1.21 | 1.24 | 0.94 | 1.09 | 1.03 | 1.13 | 0.84 | 1.09 | 1.02 | 1.19 | 1.04 | 0.94 | ||
17130 − Ankara Center | P (mm) | Trend length | 9.50 | 8.19 | 4.53 | 14.63 | 23.75 | 29.45 | 3.08 | 8.00 | 16.12 | 6.00 | 19.82 | 8.41 |
Trend slope | 8.98 | 43.69 | 6.06 | 0.37 | 0.74 | 1.31 | − 0.20 | 0.29 | 4.27 | 0.23 | 0.52 | − 0.03 | ||
T (°C) | Trend length | 3.58 | 9.87 | 6.66 | 5.32 | 4.94 | 4.19 | 0.46 | 4.90 | 7.34 | 9.74 | 9.28 | 3.40 | |
Trend slope | 1.48 | 0.86 | 0.84 | 1.03 | 1.32 | 0.90 | − 0.61 | 0.85 | 0.91 | 1.30 | 0.86 | 0.80 | ||
RH (%) | Trend length | 7.03 | 11.72 | 8.86 | 2.56 | 7.06 | 10.82 | 2.06 | 6.85 | 15.22 | 15.65 | 11.72 | 0.24 | |
Trend slope | 1.89 | 0.77 | 0.61 | 1.49 | 0.61 | 0.96 | 0.19 | 0.58 | 1.11 | 0.73 | 0.85 | − 0.19 | ||
ET (mm) | Trend length | 3.80 | 18.92 | 40.69 | 53.15 | 40.06 | 41.64 | 12.31 | 64.46 | 53.73 | 42.50 | 20.48 | 4.70 | |
Trend slope | 1.91 | 1.23 | 0.88 | 0.97 | 1.02 | 1.10 | 1.17 | 0.96 | 1.02 | 1.20 | 0.87 | 0.73 | ||
17664 − Ankara − Kizilcahamam | P (mm) | Trend length | 23.46 | 10.75 | 4.98 | 10.77 | 24.71 | 31.54 | 3.89 | 3.00 | 19.23 | 15.54 | 51.81 | 12.74 |
Trend slope | 0.63 | − 1.36 | 1.06 | − 0.33 | 0.80 | 0.76 | − 8.45 | − 1.78 | 1.64 | − 0.06 | 0.94 | 2.20 | ||
T (°C) | Trend length | 3.56 | 10.03 | 6.91 | 4.99 | 4.72 | 4.97 | 0.25 | 4.28 | 7.36 | 10.28 | 9.37 | 3.68 | |
Trend slope | 0.87 | 0.82 | 1.01 | 1.36 | 1.01 | 1.23 | − 38.75 | 1.20 | 1.05 | 0.91 | 1.11 | 0.68 | ||
RH (%) | Trend length | 4.85 | 9.18 | 5.21 | 0.38 | 3.26 | 9.65 | 1.06 | 4.01 | 11.46 | 9.40 | 9.31 | 2.02 | |
Trend slope | 3.06 | 1.19 | 1.14 | −2.96 | 0.46 | 1.65 | 0.31 | 1.64 | 1.33 | 0.91 | 1.27 | 0.20 | ||
ET (mm) | Trend length | 3.14 | 18.02 | 8.12 | 47.33 | 35.71 | 38.04 | 13.02 | 57.04 | 48.00 | 41.03 | 17.25 | 4.03 | |
Trend slope | 1.46 | 0.98 | 0.95 | 1.22 | 1.00 | 1.15 | 0.78 | 1.18 | 1.03 | 1.15 | 0.99 | 1.02 |
Station No . | Variables . | Tests . | Jan–Feb . | Feb–Mar . | Mar–Apr . | Apr–May . | May–Jun . | Jun–Jul . | Jul–Aug . | Aug–Sep . | Sep–Oct . | Oct–Nov . | Nov–Dec . | Dec–Jan . |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
17128 – Ankara Airport | P (mm) | Trend length | 15.33 | 8.40 | 4.99 | 11.11 | 20.30 | 31.74 | 2.99 | 7.11 | 19.53 | 4.95 | 29.62 | 7.79 |
Trend slope | 0.95 | − 4.41 | 0.50 | 0.16 | 0.23 | 1.19 | − 0.47 | − 0.01 | 2.87 | − 0.04 | 0.56 | 1.74 | ||
T (°C) | Trend length | 4.23 | 9.57 | 6.51 | 4.74 | 5.03 | 4.71 | 0.24 | 4.50 | 7.02 | 10.44 | 8.48 | 4.27 | |
Trend slope | 1.21 | 0.88 | 0.90 | 1.31 | 0.99 | 1.15 | − 3.80 | 0.98 | 1.11 | 0.96 | 1.23 | 0.68 | ||
RH (%) | Trend length | 6.35 | 9.89 | 6.93 | 0.81 | 6.55 | 12.96 | 1.26 | 6.26 | 13.92 | 13.52 | 10.77 | 0.90 | |
Trend slope | 2.07 | 1.09 | 0.59 | − 2.28 | 0.51 | 1.02 | 1.49 | 0.62 | 1.39 | 0.86 | 1.11 | 4.79 | ||
ET (mm) | Trend length | 3.61 | 18.27 | 38.84 | 48.23 | 36.25 | 40.87 | 11.19 | 60.86 | 50.77 | 43.20 | 16.66 | 3.39 | |
Trend slope | 1.21 | 1.24 | 0.94 | 1.09 | 1.03 | 1.13 | 0.84 | 1.09 | 1.02 | 1.19 | 1.04 | 0.94 | ||
17130 − Ankara Center | P (mm) | Trend length | 9.50 | 8.19 | 4.53 | 14.63 | 23.75 | 29.45 | 3.08 | 8.00 | 16.12 | 6.00 | 19.82 | 8.41 |
Trend slope | 8.98 | 43.69 | 6.06 | 0.37 | 0.74 | 1.31 | − 0.20 | 0.29 | 4.27 | 0.23 | 0.52 | − 0.03 | ||
T (°C) | Trend length | 3.58 | 9.87 | 6.66 | 5.32 | 4.94 | 4.19 | 0.46 | 4.90 | 7.34 | 9.74 | 9.28 | 3.40 | |
Trend slope | 1.48 | 0.86 | 0.84 | 1.03 | 1.32 | 0.90 | − 0.61 | 0.85 | 0.91 | 1.30 | 0.86 | 0.80 | ||
RH (%) | Trend length | 7.03 | 11.72 | 8.86 | 2.56 | 7.06 | 10.82 | 2.06 | 6.85 | 15.22 | 15.65 | 11.72 | 0.24 | |
Trend slope | 1.89 | 0.77 | 0.61 | 1.49 | 0.61 | 0.96 | 0.19 | 0.58 | 1.11 | 0.73 | 0.85 | − 0.19 | ||
ET (mm) | Trend length | 3.80 | 18.92 | 40.69 | 53.15 | 40.06 | 41.64 | 12.31 | 64.46 | 53.73 | 42.50 | 20.48 | 4.70 | |
Trend slope | 1.91 | 1.23 | 0.88 | 0.97 | 1.02 | 1.10 | 1.17 | 0.96 | 1.02 | 1.20 | 0.87 | 0.73 | ||
17664 − Ankara − Kizilcahamam | P (mm) | Trend length | 23.46 | 10.75 | 4.98 | 10.77 | 24.71 | 31.54 | 3.89 | 3.00 | 19.23 | 15.54 | 51.81 | 12.74 |
Trend slope | 0.63 | − 1.36 | 1.06 | − 0.33 | 0.80 | 0.76 | − 8.45 | − 1.78 | 1.64 | − 0.06 | 0.94 | 2.20 | ||
T (°C) | Trend length | 3.56 | 10.03 | 6.91 | 4.99 | 4.72 | 4.97 | 0.25 | 4.28 | 7.36 | 10.28 | 9.37 | 3.68 | |
Trend slope | 0.87 | 0.82 | 1.01 | 1.36 | 1.01 | 1.23 | − 38.75 | 1.20 | 1.05 | 0.91 | 1.11 | 0.68 | ||
RH (%) | Trend length | 4.85 | 9.18 | 5.21 | 0.38 | 3.26 | 9.65 | 1.06 | 4.01 | 11.46 | 9.40 | 9.31 | 2.02 | |
Trend slope | 3.06 | 1.19 | 1.14 | −2.96 | 0.46 | 1.65 | 0.31 | 1.64 | 1.33 | 0.91 | 1.27 | 0.20 | ||
ET (mm) | Trend length | 3.14 | 18.02 | 8.12 | 47.33 | 35.71 | 38.04 | 13.02 | 57.04 | 48.00 | 41.03 | 17.25 | 4.03 | |
Trend slope | 1.46 | 0.98 | 0.95 | 1.22 | 1.00 | 1.15 | 0.78 | 1.18 | 1.03 | 1.15 | 0.99 | 1.02 |
At station 17128, a sharp transition decreasing trend area to an increasing area is observed from February to March, while an increasing trend area to the decreasing area is seen in the transition from August to September for mean precipitation data. The winter (DJF) season is in the decreasing area far away from 45° (1:1 line). October appears to be a strong increasing trend, while April, May, and September showed a significant decreasing trend. The maximum trend length for precipitation data is evaluated as 31.74 mm in the June–July transition zone, whereas the lowest trend length is calculated as 4.95 mm within the October–November transition zone. The maximum trend slope is seen in the February–March transition zone as −4.41 and the lowest trend slope is evaluated as −0.01 within the August–September transition zone. According to temperature data, all months are in an increasing transition zone except December. In other words, the fours season is almost in an upward trend with far away from 45° (1:1 line). The trend length varies from 0.23 to 10.43 (°C) and the trend slope ranges from −3.8 to 1.308. The maximum trend length is seen in the October–November transition zone. The mean value of evapotranspiration data is similar to the temperature data as a significantly increasing trend zone. The maximum trend length is observed as 50.77 in the September–October transition zone, while the lowest trend length is seen in January–February as 3.61. The trend slope changes from 0.83 (July–August) to 1.24 (February–March). A sharp transition decreasing trend area is observed for relative humidity data in the spring (MAM) season, January–February, and August–September periods, a notable increasing trend appears in October. The trend slope varies between −2.27 (April–May) and 4.78 (December–January). The trend length ranges from 0.81 (April–May) to 13.91 (September–October).
At station 17130, rainfall has a transition from September to October in the decreasing trend zone to the increasing trend zone, while January–February, August–September, and November–December are in the significantly increasing trend zone. Rainfall in May and June are clustered around 45° (1:1 line) which means no trend. October–November has seen a strong increasing trend according to Figure 6. The maximum trend length is calculated as 23.75 in May–June transition zone, while the lowest value is seen as 4.53 from March to June periods. The trend slope for rainfall data varies from −0.02 (December–January) to 43.68 (February–March). According to temperature data, the winter (DJF) season is observed in the increasing area with far away from 45° (1:1 line), a notable decreasing trend is just seen in the November–December transition zone. The temperature in April, May, and August is clustered around no trend line. According to Table 4, the maximum trend length occurred as 9.86 from February to March, while the minimum trend length is calculated as 0.46 from July to August. The highest trend slope appears as 1.48 (January–February), notable the lowest trend slope is seen as −0.61 (July–August). Evaporation data for the same station, there is a distinct transition from October to November from the significantly increasing to decreasing trend area, while no trend is observed in May, June, and December. The trend slope ranges from 1.90 (January–February) to 0.73 (December–January). The maximum trend length varies from 53.72 (September–October) to 3.80 (January–February). In addition, the minimum trend slope and trend length are calculated in January–February transition zone. According to relative humidity data, the winter (DJF) season is evaluated in the decreasing area with far away from 45° (1:1 line), whereas the summer (JJA) season shows in the increasing area with the September–October periods. The maximum trend length is seen as 15.64 from October to November, while the minimum trend length is calculated as 0.23 from December to January. Trend slope for relative humidity data changes from −0.18 (December–January) to 1.89 (January–February).
At station 17664, the winter (DJF) season is in the decreasing transition area with May–June, August, and November periods, whereas the March–April, July, and September–October periods are in the significantly increasing for precipitation data. A sharp increasing trend from an upward to downward trend zone is observed in the transition from April to May, July to August, and October to November. A larger trend is clearly seen in October–November than in another month in the upward trend zone but in a decreasing direction. The maximum trend length for precipitation is calculated as 51.81 from November to December, notable the minimum decreasing trend length is seen as 2.99 from August to September. The highest trend slope is shown as −8.45 (July–August), while the lowest trend slope is seen as −0.06 (October–November). According to temperature data, there is one transition (February–March) that appears from an increasing trend area to a decreasing trend area. In general, the remaining months are clustered at an upward trend area except for May–June, September–October, and December which are near to no trend lines. Temperature from July to August has the minimum trend length and the maximum trend slope as 0.25 and −38.75, respectively. In addition, the maximum trend length and the minimum trend slope are observed as 10.08 (October–November) and 0.6 (December–January). According to evaporation data for the same station, from April to May is also one observed transition from a downward trend area to an upward area. The summer (JJA) season is shown in a significantly increasing trend, while the winter (DJF) season are observed as no trend due to clustering around 45° (1:1 line). The maximum trend length for ET data is calculated as 57.03 from August to September, whereas the minimum trend length is seen as 3.14 from January to February. From July to August is evaluated as the minimum trend slope with 0.78, notable that the maximum trend slope is observed from January to February as 1.46. According to relative humidity data, spring (MAM) and summer (JJA) seasons are shown to a decrease trend zone significantly. The increasing trend is observed in January, October, and December at the upward trend area. The trend length ranges from 0.38 (April–May) to 11.46 (September–October). The trend slope varies from 0.195 (December–January) to 3.05 (January–February).
A general assessment of trend tests (MK, ITA, and IPTA)
The results of MK, ITA, and IPTA tests are presented in Tables 5,67–8 for each hydrometeorological variable. When the tables are examined in general, there are differences between the MK test and the innovative tests. According to Table 5, only 1 of 36 months has a significant trend (decreasing/increasing) when using rainfall data, according to the MK test. However, IPTA and ITA capture a significant trend in 30 and 34 of 36 months, respectively. Innovative approaches caught the same trends as MK captured in the months when MK captured a trend. The ITA and IPTA, on the other hand, exhibit trends for the large majority of the 35 months where MK did not find any significant trend for precipitation data. The explanation for the difference in results between MK and innovative methods is that innovative methods are more sensitive in detecting trends (Achite et al. 2021; Şan et al. 2021; Akçay et al. 2022). The hidden trend in meteorological variables is detected effectively using graphical methods rather than the classical trend method.
According to three approaches, it is obtained that there is an increasing trend in February, Jun, July, and September in temperatures and a decreasing trend in February and March in relative humidity for station 17228. Evaporation for the same station, the summer (JJA) season, February, May, September, and October indicate a significantly increasing trend considering all three methods. At station 17130, precipitation in March shows an increasing trend, while no trends are detected in May and July. Increasing trends are detected in January, February, June, and July, whereas April, May, and August show no significant trend for temperature. For relative humidity, the decreasing trend that begins in the winter season continues until the end of March. However, the summer season shows an increasing trend that continues until the end of October in terms of MK, ITA, and IPTA methods. Furthermore, there is no trend found in May and November. For evapotranspiration, an increasing trend is observed in January, February, March, July, and September, while no trend just is detected in May. According to station 17664, precipitation and temperature indicate the difference between MK and innovative approaches. In other words, the innovative method shows a significant trend, while the MK test does not capture in the same manner. Relative humidity in January is seen as an increasing trend by considering three methods. For evapotranspiration, an increasing trend is detected in the summer season, February and September. In addition, February, July, and September are observed as increasing for all three methods as well as for all stations.
Change point detection results
Station . | Meteorological variables . | SNHT . | CUSUM test . | |||
---|---|---|---|---|---|---|
Tk . | Year . | Vk . | α = %5 . | Year . | ||
17128 | Annual total precipitation (mm) | 5.0448 | – | 6 | 11.04 | – |
Annual maximum temperature (°C) | 22.331 | 1997 | 20 | 11.04 | 1993 | |
Annual mean relative humidity (%) | 29.215 | – | 10 | 11.04 | – | |
Annual total evapotranspiration (mm) | 14.881 | 2015 | 20 | 11.04 | 2011 | |
17130 | Annual total precipitation (mm) | 12.696 | 1960 | 15 | 13.35 | 1960 |
Annual maximum temperature (°C) | 17.596 | 1997 | 18 | 13.35 | 1997 | |
Annual mean relative humidity (%) | 12.584 | – | 12 | 13.35 | – | |
Annual total evapotranspiration (mm) | 39.276 | 1997 | 18 | 13.35 | 1997 | |
17664 | Annual total precipitation (mm) | 4.8957 | – | 16 | 10.79 | 2008 |
Annual maximum temperature (°C) | 12.701 | – | 48 | 10.79 | 2005 | |
Annual mean relative humidity (%) | 4.1061 | – | 7 | 10.79 | – | |
Annual total evapotranspiration (mm) | 23.283 | 2005 | 16 | 10.79 | 2005 |
Station . | Meteorological variables . | SNHT . | CUSUM test . | |||
---|---|---|---|---|---|---|
Tk . | Year . | Vk . | α = %5 . | Year . | ||
17128 | Annual total precipitation (mm) | 5.0448 | – | 6 | 11.04 | – |
Annual maximum temperature (°C) | 22.331 | 1997 | 20 | 11.04 | 1993 | |
Annual mean relative humidity (%) | 29.215 | – | 10 | 11.04 | – | |
Annual total evapotranspiration (mm) | 14.881 | 2015 | 20 | 11.04 | 2011 | |
17130 | Annual total precipitation (mm) | 12.696 | 1960 | 15 | 13.35 | 1960 |
Annual maximum temperature (°C) | 17.596 | 1997 | 18 | 13.35 | 1997 | |
Annual mean relative humidity (%) | 12.584 | – | 12 | 13.35 | – | |
Annual total evapotranspiration (mm) | 39.276 | 1997 | 18 | 13.35 | 1997 | |
17664 | Annual total precipitation (mm) | 4.8957 | – | 16 | 10.79 | 2008 |
Annual maximum temperature (°C) | 12.701 | – | 48 | 10.79 | 2005 | |
Annual mean relative humidity (%) | 4.1061 | – | 7 | 10.79 | – | |
Annual total evapotranspiration (mm) | 23.283 | 2005 | 16 | 10.79 | 2005 |
There is one intersection point of and in annual relative humidity series except station 17664. No upward/downward trends are detected as well as no change point by three tests. In general, annual total evaporation series are an upward trend after 2005 by the SQ-MK test. Similarly, an upward trend is observed in the annual maximum temperature that the start of 1997. However, there are no significant trends detected in both annual precipitation (except station 17130) or annual relative humidity series.
CONCLUSION
Ankara is one of the most critical provinces in Turkey. Therefore, this paper aims to employ new graphical trend analysis as well as the classical method. The three most important meteorological stations evenly distributed in Ankara province are selected. For this purpose, hydrometeorological variables including precipitation, temperature, relative humidity, and evapotranspiration trends are investigated at monthly and annual scales. SNHT, CUSUM, and SQ-MK methods are used to utilize any abrupt change points in the annual series. IPTA (one of the most recent trend analysis approaches reported in the literature), ITA with significance test, and MK (one of the classical approaches of trend analysis) are applied to detect the presence of trends at both monthly and annual scales. The findings of the trend analysis are also compared. While the MK test fails to detect significant trends in the majority of variables, ITA and IPTA capture different trends. This paper results are convenient with previous studies in the literature (Caloiero et al. 2018; Ceribasi & Ceyhunlu 2020; Hırca et al. 2022). In other words, ITA with significance test, IPTA, and MK capture trends in 83, 94, and 0.2% of all months (3 meteorological stations × 12 months) for precipitation, 86, 75, and 22% of all months for temperature, 80, 80 and 30% of all months for relative humidity, and 91, 80 and 47% of all months for evapotranspiration series. These findings indicate that ITA with significance test and IPTA are more sensitive than the MK test with the highest in the precipitation series and lowest in the evapotranspiration series. In addition, three methods are convenient for capturing trends in evapotranspiration series due to higher trends. Thus, the MK test is unable to capture the hidden trend in time series, especially in precipitation series. Therefore, the hidden trend in precipitation is successfully captured by ITA and IPTA graphical methods. Previous studies (Şan et al. 2021; Hırca et al. 2022) also indicate the same result for precipitation data. In addition, although the trend captured by the MK method is observed by the IPTA and ITA methods, the trends that the MK method does not show in some months are shown by the IPTA and ITA methods. The advantage of IPTA graphs shows the trend slope as well as the trend length compared with the MK test. According to IPTA, MK and ITA with significance test indicate a holistic trend. According to the results, for station 17228, it is discovered that there is an increasing tendency in temperature in February, June, July, and September, and a decreasing trend in relative humidity in February and March. Considering all three approaches, evaporation for the same station during the summer (JJA) season shows a significantly increasing trend. Precipitation in March shows an increasing trend, while no trends are detected in May and July for station 17230. The decreasing trend in relative humidity that begins in the winter season continues through the end of March. In terms of MK, ITA, and IPTA procedures, however, the summer season shows an increasing trend that will last through the end of October. In the summer season, there is an increasing trend in evapotranspiration. Furthermore, all three methodologies as well as all stations show an increasing trend in February, July, and September. According to annual time series findings, at stations 17128 and 17664, there is no apparent trend in annual total precipitation, however, at station 17128, there is an increasing trend. At annual mean relative humidity, there is no indication of an increasing or decreasing trend at any of the stations. However, at annual total evapotranspiration, all stations show an increasing trend (10 and 5% confidence levels). According to change point detection, the results reveal that for station 17128 in the annual maximum temperature and annual evapotranspiration series, the SNHT and CUSUM tests fail to detect any abrupt changes. For station 17664, three tests successfully captured a change point in annual evapotranspiration in 2005. According to the SQ-MK test, annual total evaporation series have been increasing since 2005. Similarly, from the beginning of 1997, an increasing trend in annual maximum temperature has been noticed. However, neither the annual precipitation (except station 17130) nor the annual relative humidity series show any significant changes.
FUNDING
No funding was received for conducting this study.
AUTHOR CONTRIBUTIONS
M.E. data gathering, hydrometeorological data trend analysis, interpretation of the findings, manuscript writing, and submission.
DATA AVAILABILITY STATEMENT
All relevant data are included in the paper or its Supplementary Information.
CONFLICT OF INTEREST
The authors declare there is no conflict.