Homogeneity and Trend Analysis of Hydrometeorological Data of the Eastern Black Sea Region , Turkey

Eastern Black Sea Region in northeastern part of Turkey has the highest precipitation total in the country, approaching 2500 mm per a year. It is therefore an important region as it frequently encounters with flash floods due to heavy rains. For future planning of water resources, environment and urbanization, it is important to know the expected behavior of hydrometeorological processes, mainly precipitation and flow. Due to these facts, in this study, homogeneity of long-term annual precipitation and streamflow series of the Eastern Black Sea Region, Turkey is checked using double mass curve method and trends are determined by means of the Mann-Kendall test. The data network consists of 38 precipitation gauging stations and 40 flow gauging stations across the Eastern Black Sea Region. It is found that 27 precipitation stations out of 38 are homogeneous and no trend is available. Out of the remaining stations, nine are found non-homogeneous and four with trend. For annual flow data, it is found that 22 stations out of 40 are homogeneous and no trend is available. The remaining 18 stations are found non-homogeneous, among which 5 stations have trend at the same time.


Introduction
Use of precipitation and flow data in water resources assessment, hydrological modeling, climate change studies, and urban and environmental planning of a region is a must.Network of precipitation and flow stations is not well-distributed in undeveloped/developing countries [1] and in mountainous terrains such that even a station record is valuable.On the other hand, climate change, measurement errors, extreme meteorological events affect structural characteristics of the hydrometeorological time series such as homogeneity and possible trends.In order to use as many records as available for a region, structural characteristics of the data should be derived and evaluated.
There have been many studies regarding homogeneity and particularly trend availability of various hydrometeorological data for Turkey.The homogeneity of hydrometeorological variables such as precipitation [2,3], temperature and relative humidity [3] has been previously investigated for Turkey.Similarly, trend analysis has been performed on precipitation [4][5][6][7][8], temperature [6,[9][10][11][12], solar radiation [13], streamflow [6,[14][15][16][17].The results of these studies vary depending on the number of stations, record period and statistical test used.A common point is to classify a station as homogeneous or non-homogeneous and/or to determine trend availability and characteristics of the observed data mostly for evidence of climate change.However, no comparison has been made between original (observed) and homogenized and trend-removed data.Such a comparison provides one to use as many records as available, instead of ignoring any data which have inconsistency and/or trend.
Therefore, in this study, first, the homogeneity and trend availability of the precipitation and streamflow data for the Eastern Black Sea Region, Turkey are investigated.The results of the data analysis are then compared to that of aforementioned studies.Secondly, the observed data from non-homogeneous stations are homogenized and if available, trends are removed.By this way, the difference between observed and adjusted data can be evaluated.

Study Area and Data
The coastal part of the Eastern Black Sea Region which is located in the north east of Turkey, between 40˚31' -41˚24'N and 38˚08' -41˚26'E is selected as study area.This coastal part of the region covers an area bordered with the Eastern Black Sea Mountain chain to the south and the Black Sea to the north as seen in Figure 1.These high mountain ranges run parallel to the sea coast as the north boundary of the study area, and rise to more than 3000 m above mean sea level (MSL).The Black Sea Region has a steep rocky coast with some rivers that cascade through the gorges of the coastal ranges.
In the coastal area of the Eastern Black Sea Region, mild and humid climate dominates.Snowfall may be seen in winter.Yearly average temperature is about 14˚C -15˚C in the coastline, however it decreases with increasing elevation.The average precipitation of the coastal area of this region is more than 1000 mm.
Mean annual precipitation and mean annual flow observations are used in this study.Precipitation data were taken from 38 observation stations of which 19 are located on the coastline of the area.Mean annual flow observations are obtained from 40 flow stations.Locations of the stations are shown in Figure 1 where precipitation stations are numbered from 1 to 38 and flow stations from 1 to 40 starting from west to east.Characteristics of the observation stations are shown in Table 1.Numbers in the column of "No" in Table 1 corresponds to numbers on the map in Figure 1.
The length of precipitation data ranges from 10 to 46 years between 1960 and 2005, whereas the flow record length ranges from 10 to 49 years between 1944 and 2006 with some gaps in the data.To complete the gap in any station, regression equations were developed using continuous data from the neighboring stations.The observed flow is not influenced by any upstream dam or water.

Tests Employed
The homogeneity of precipitation and flow data is checked using double mass curve.The double mass curve is a well-known simple tool for checking and adjusting inconsistencies in hydrometeorological data caused by changes in observation methods or data processing.The theory of double mass curve is based on a graph of accumulation of one quantity against the accumulation of another during the same period.If the process is homogeneous, the graph plots a straight line.A break in slope of the line indicates a change in the constant of proportionality between the two variables [18].In hydrological studies, double mass curve may give indefinite results such that it is unable to say which of the variables caused a break in the slope.In order to obtain more definite results, the accumulations of one of the variables can be plotted against the accumulations of a pattern composed of all similar records in a given area.
For determination of trends in precipitation and flow data for the region the Mann-Kendall test is employed.The Mann-Kendall test is a non-parametric rank-based method which does not require the normal distribution of data.It is robust to the effects of outliers and allows existence of missing data (as only ranks are used).The basic principle of Mann-Kendall test for detecting a trend in a time series is to examine the sign of all pair wise differences of observed values.The Mann-Kendall test has been widely used to detect trends in hydro-meteorologycal time series [5,14,[19][20][21][22][23] and is well documented in the literature [24,25].

Results and Discussion
Results of the tests employed for the precipitation data are analyzed.It is found that 29 stations out of 38 are homogeneous with no change in the slope of the double mass curves.For the remaining 9 precipitation stations, there are inconsistencies observed.To adjust inconsistency in any precipitation station, a correction factor is estimated by comparing cumulative annual mean precipitation data of non-homogenous stations with remaining stations, thus the precipitation data can be considered as homogenous after correction.An example of this process is shown in Figure 2 for Tamdere (3) precipitation station [Numbers in parenthesis show station number].Mann-Kendall trend test is then applied to the corrected precipitation data of 9 precipitation stations together with that data of 29 homogenous stations.Upward trends are found in Giresun ( 4) and Trabzon ( 15) stations whereas downward trends are observed in Tamdere (3) and Macka (14).
In a previous study [5] in which Giresun (4), Trabzon (15) and Rize (26) were used as common stations, no trend was found in Giresun (4) and Trabzon (15) whereas a trend was obtained in Rize (26).The beginning of the trend for Rize (26) was determined as 1952.Note also that the data record length in the study by [5] ranges from 1929-1993, while in this study it covers only the years between 1960 and 2005.
Comparison of the results of previous and current study are given in Figure 3 including record periods and stations used.Available upward/downward trends are removed from the precipitation data of 4 stations in order to compare adjusted data with observations (Figure 4).In this manner, in addition to determination of non-homogenous and/or trend-available stations, the difference between adjusted and observed data can clearly be recognized.
Positive and negative differences between adjusted and observed annual precipitation data are found as +16.2% and -27.6%, respectively.The difference between observed and adjusted data is shown in Figure 3 with the results of available upward/downward trends.In Figure 3, stations are arranged in ascending order based on their precipitation values.
For streamflow data, it was found that 22 stations out of 40 are homogeneous.For the remaining 18 stations, streamflow data are non-homogeneous.Streamflow data has been corrected using double mass curve in these 18 stations, Mann-Kendall trend test is then performed to all flow stations.Only downward trends are observed in Kanlipelit (14), Findikli (18), Agnas (19), Toskoy (28) and Topluca (34) stations.
In a previous study [15], no trend was found in the Eastern Black Sea Region according to the results of the parametric and Mann-Kendall tests applied on mean annual flow data of 12 flow stations operated by EIE with record lengths changing from 25 to 66 years.
In another previous study [14] on trend analysis of streamflow in Turkey, Dereli (4), Simsirli (26), Toskoy (28) and Topluca (34) were used as common stations, no trend was found.However, in this study a downward trend was found in stations Toskoy (28) and Topluca (34).Note that the data record length in [14] ranges from 1964 to 1994 while in this study, it covers a longer period between 1944 and 2006.From Table 1, one can realize that there are two stations named Toskoy (28 and 30) on the same stream.A trend was found in station (28) whereas no trend was available in (30).Trend was found when data in (28) was homogenized.In addition, data length is 38 years from 1965 to 2002 for station (28) and 10 years from 1986 to 2001 for (30).This shows the effect of data length on trend analysis and also depicts how controversial results can be obtained for the same region.
Reference [16] studied the trend detection over Turkey.For Eastern Black Sea Region, mean annual flow from the stations Agnas (19), Dereli (4), Simsirli (26), Topluca (34), and Toskoy (28) were used.Insignificant downward trend was found in the stations (19,34,28) and downward trend determined in (26) was found significant.No trend was found in (4).It should be pointed out once again that the data period is between the years of 1968-97.
After removing available trends from the data, positive and negative differences between adjusted and observed annual mean streamflow are found as +17.5% and -8.6%, respectively.The difference between observed and adjusted data is shown in Figure 5 together with the results of available downward trends.In Figure 5, stations are arranged in ascending order based on their flow values.

Conclusions
In this study, homogeneity and trend-availability of mean annual precipitation and streamflow data of the Eastern Black Sea Region, Turkey were investigated by means of double mass curve method and Mann-Kendall trend test.The results are compared to previous studies on precipitation and streamflow data network across Turkey for different time periods.This comparison shows the effect of data length on trend analysis and also depicts how different results can be obtained for the same region.
Precipitation and streamflow data which have been found non-homogeneous, were first homogenized and trend, if available, was removed.The differences between observed and adjusted mean annual precipitation and streamflow data were calculated.Instead of only detection of inconsistency and/or trend in a station record, determination of difference between adjusted data and  observations is more helpful to decide which station record should be used or ignored in any hydrological and water resources application.For undeveloped and developing countries, and particularly in mountainous regions where the station network is generally inadequate, rather than ignoring any non-homogenous station record, it is suitable to follow the homogeneity and trend analysis results.When the difference between observed and adjusted data is below a specified limit, then the data set can be used with their original observation values.

Figure 1 .
Figure 1.Study area and locations of the precipitation and flow stations.

Figure 3 .
Figure 3. Results of previous and current study and the difference between original and adjusted precipitation data in percent.

Figure 5 .
Figure 5. Results of previous and current study and the difference between original and adjusted flow data in percent.

Table 1 . Characteristics of precipitation and flow stations.
a,b DMI (State Meteorological Service), DSI (State Hydraulics Works), EIE (Electrical Power Resources Survey and Development Administration) with Turkish acronym.