Long Term Historic Changes in the Flow of Lesser Zab River , Iraq

The assessment of trends in river flows has become of interest to the scientific community in order to understand the changing characteristics of flow due to climate change. In this study, the trends in river flow of Dukan Dam located in the northern part of Iraq were assessed. The assessment was carried out for the period 1964 to 2013 using Sen’s slope and the Mann–Kendall test. Sen’s slope was used to assess the magnitude of change while the Mann–Kendall trend test was used to confirm the significance of trends. The results of the study showed that there was a decreasing trend in river flow both annually and for all individual months. The highest decreasing trend of −5.08846 m3/month was noticed in April, while the lowest change of −1.06022 m3/month was noticed in November. The annual flow also showed a significant decrease at a rate of −1.912 m3/year at a 95% level of confidence. Additionally, the findings of the study also confirmed that a decrease in precipitation and the construction of hydraulic structures reduced the flow in the river. The findings of the study suggest that decreasing trends may cause a water-scarce situation in the future if proper adaptation measures are not taken.


Introduction
High concentrations of greenhouse gases are causing an increase in global temperature [1][2][3][4].The increase in temperature is playing an adverse role in altering the components of the global hydrological cycle [5][6][7][8][9].Precipitation is considered to be the most vigorous component of the global hydrological cycle, which has reportedly altered over several regions around the world [10][11][12][13].The alteration of precipitation characteristics can have severe impacts on society in the form of floods and droughts, which can have an adverse impact on people's socio-economic situation [14][15][16][17].Changes in streamflow or river flow are often associated with changes in precipitation.Cigizoglu,et al. [18] reported that stream flow changes are very sensitive to even small changes in precipitation.Therefore, it is very important and informative to investigate trends in river flow.
Trend assessments of river flow provide basic information related to water [18][19][20][21][22][23][24][25].Several trend assessments have shown that changes in the characteristics of temperature and precipitation have direct impacts on river flows across the world [26][27][28][29][30]. Iraq, located in Asia, has also experienced changes in precipitation and temperature in the last fifty years [31][32][33].Therefore, it can be remarked that climate change may have impacted the water resources of Iraq.Additionally, Shubbar, Salman and Lee [33] reported that the temperature in Iraq will increase in future.Future temperature increases will modify precipitation patterns, which will affect the major water reservoirs.Thus, it is important to evaluate the trends in the major water reservoirs of Iraq.Knowledge of the changes in river flow in Iraq is very limited due to the availability and quality of data [31].However, it is one of the most water stressed countries in the region.The gradual decline of water resources due to a decline in rainfall has been reported in Iraq [32].This reinforces the need for a quantification of the changes in river flow in the country.Such analysis should be conducted with robust statistical methods considering the quality of data and missing values.The present study considered non-parametric methods that can provide trends and their significance when data quality is doubtful.
A number of dams have been constructed in the northern region of Iraq, such as in Mosul, Duhok, and Dukan.Dukan Dam, constructed over the lower Zab catchment is playing a vital role in supplying water to a large population in Iraq [34,35].However, in recent years changes in the characteristics of precipitation have been reported over the region, which is expected to have an impact on the flow of the river [32,36].Therefore, it is important to evaluate trends in the flow of Dukan Dam in order to avoid major crises regarding water supply in the region in the future.A review of the literature revealed that no studies has evaluated the river flow trends in the Dukan Dam flow.Thus the major objective of the study was to assess the long-term historic trends in the flow of the Lesser Zab river.
A number of trend assessment procedures are available in the literature, such as linear regression, Bayesian, and Mann-Kendall tests [37,38].The Mann-Kendall test is considered to be a robust method for trend assessment and is also recommended by the World Meteorological Organization [39][40][41].Thus, in this study the Mann-Kendall trend test was used for the confirmation of trends in the river flow (monthly, mean annual, annual maximum and annual minimum flows).The least square and Sen's slope tests are most widely used for estimating the magnitude of a trend.However, Sen's slope is not influenced by extreme values and thus is more robust compared to other methods.Thus, in this study, Sen's slope was used to estimate the magnitude of the change in trends.The assessment was carried out on data from fifty years ranging from 1964 to 2013.

Study Area and Datasets
Dukan Dam (Figure 1) is a multi-purpose cylindrical concrete arch dam located in the northern regions of Iraq (latitude: 35 • 57'15" N and longitude: 44 • 57'10" E) [42].The dam is 360 m long, 116.5 m high, and 34.3 m and 6.2 m wide at the base and crest, respectively.The dam has a catchment size of 11,690 km 2 , with a surface area of 270 km 2 .The normal operational level of the dam is 511 m above the mean sea level.The dam was constructed across the Lesser Zab in 1954 to ensure water supply to the population for domestic, industrial and agricultural purposes [42,43].The dam is now also used for controlling floods and power generation [44].The dam has the capacity to impound 6,100,000,000 m 3 of water with two spillways, i.e., a tunnel and a bell-mouth.The tunnel is used to discharge excess water while the bell-mouth is operated during emergency periods.One of the major advantages of the reservoir is its capacity and gravel type bed, which helps to increase the water level in the area.
The study area experiences two climatic seasons, summer and winter.Summer begins in June and lasts until September, while winter begins in November and lasts until March.Around 90% of the rainfall occurs from November to April, and the rest of the months are mostly dry [32].Rainfall in the area varies from 200 mm to 1000 mm [35].The land use/land cover of the study area is mainly characterized by winter plants and pastures (82.7%), forest (15.6%), vegetables (1.6%) and urban areas (0.1%) [42].
The monthly river flow data recorded at two gauging stations were acquired from the Directorate of Dukan Dam.The data was collected for the period 1964 to 2013.The percentage of missing data was found to be 0.67% (upstream side) and 2.83% (downstream side) of the dam.Missing values in any dataset can be estimated by various statistical methods, for example, arithmetic averaging [45], inverse distance interpolation [46], multiple regression analysis [47], closest station method [48], artificial neural network [49,50] and expectation maximization (EM) algorithm [51].In this study, missing data were filled using the EM approach [52].EM has been widely used in recent years for filling missing data in hydro-climatological records [53,54].The completed data series were arranged in Excel for further analysis.The monthly flow of each year was used to compute the mean annual, annual maximum and annual minimum flow.
arithmetic averaging [45], inverse distance interpolation [46], multiple regression analysis [47], closest station method [48], artificial neural network [49,50] and expectation maximization (EM) algorithm [51].In this study, missing data were filled using the EM approach [52].EM has been widely used in recent years for filling missing data in hydro-climatological records [53,54].The completed data series were arranged in Excel for further analysis.The monthly flow of each year was used to compute the mean annual, annual maximum and annual minimum flow.

Distributional Analysis
Distributional analyses are frequently performed on hydrological studies to understand the nature of the data [55].A number of frequency distribution methods are available in the literature such as normal, log normal, uniform, generalized extreme value, binomial, gamma, etc.In the present study, gamma, generalized extreme value and normal distributions, which are widely used in hydro-climatic studies, were applied to the datasets of Dukan Dam.The gamma and generalized extreme value are widely used for fitting skewed data, which is often noticed in hydrological data [55].In addition, normal distribution was used as monthly river flow was often found to follow such a distribution.The goodness of fit for the tests was assessed using the Kolmogorov-Smirnov (KS) test.The KS test is found to be more powerful than the Anderson-Darling and Chi-square test [55].Therefore, the test was applied for each month, mean annual flow, annual maximum flow and annual minimum flow at a 95% level of confidence with a null-hypothesis.Therefore, the test was applied for each month, mean annual flow, annual maximum flow and annual minimum flow at a 95% level of confidence with a null-hypothesis of the data followed the particular distribution (gamma, generalized extreme value and normal distribution).

Estimation of Magnitude of Change
In this study, Sen's slope [56] method was applied to calculate the intensity of the change in river flow.The method is non-parametric and calculates the slope with a pair of equal data.The

Distributional Analysis
Distributional analyses are frequently performed on hydrological studies to understand the nature of the data [55].A number of frequency distribution methods are available in the literature such as normal, log normal, uniform, generalized extreme value, binomial, gamma, etc.In the present study, gamma, generalized extreme value and normal distributions, which are widely used in hydro-climatic studies, were applied to the datasets of Dukan Dam.The gamma and generalized extreme value are widely used for fitting skewed data, which is often noticed in hydrological data [55].In addition, normal distribution was used as monthly river flow was often found to follow such a distribution.The goodness of fit for the tests was assessed using the Kolmogorov-Smirnov (KS) test.The KS test is found to be more powerful than the Anderson-Darling and Chi-square test [55].Therefore, the test was applied for each month, mean annual flow, annual maximum flow and annual minimum flow at a 95% level of confidence with a null-hypothesis.Therefore, the test was applied for each month, mean annual flow, annual maximum flow and annual minimum flow at a 95% level of confidence with a null-hypothesis of the data followed the particular distribution (gamma, generalized extreme value and normal distribution).

Estimation of Magnitude of Change
In this study, Sen's slope [56] method was applied to calculate the intensity of the change in river flow.The method is non-parametric and calculates the slope with a pair of equal data.The method is well-known and widely used for the calculation of trends [23,24,57].The slope of a trend can be calculated using the following equation: where x j and x k are the data values at times j and k, (j > k), respectively.The median of these N values of Q is Sen's slope.If there is only one data in each time period, then: where n is the number of periods.The median of the N estimated slopes is obtained in the usual way, the N values of Q are ranked by

Mann-Kendall Trend Tests
The Mann-Kendall test was used to reveal the trends in monthly and mean annual, annual maximum and annual minimum river flow data [18, 58,59].The Mann-Kendall test (MK) is frequently used for trend assessment and is a nonparametric test [60].The test is also endorsed by the World Meteorological Organization (WMO) for hydro-meteorological trend assessments.The MK test has the ability to handle outliers and it is distribution free [16,38,61].Mann-Kendall test statistics (S) for river flow series can be calculated as: In Equation ( 4), n is data length, and x i and x j are sequential data in series, while the sign is calculated as below: The variance of the statistic can be calculated as: In Equation ( 6), t p refers to ties of the p th value, and q refers to a tied values number.Standardized test statistics for the Mann-Kendall test can be calculated as: Positive and negative µ 1 values indicate that the direction of the trend exists in a time series.

Distributional Analysis
Gamma, generalized extreme value and normal distribution were individually applied to the data for each month, the mean annual, annual maximum, and annual minimum flow.Tests were applied to the data of both gauges located on the upstream and downstream side of the dam at 95% level of confidence.The obtained results are presented in Table 1.Bold values in the table show that the hypothesis is rejected at the 95% level of confidence, while values in italic indicate the best fit.As can be seen in Table 1, the null hypothesis was accepted for most of the months in all tests.In the upstream side gauge, all data was found to follow the gamma and generalized extreme value distributions, while the null hypothesis was rejected in the month of March for the normal distribution.On the other hand, for the downstream side gauge all data was found to follow the generalized extreme value distribution.In the gamma distribution, the months of March, April, May and the annual maximum flow were found to show a rejected null hypothesis.Similarly, in the normal distribution, the months of January, March, April, May and the annual maximum flow were rejected.Overall, the generalized extreme value was found to have the best fit for the Dukan Dam data.

Seasonal Distribution of Flow
The seasonal distribution of monthly flow in the Dukan Dam is shown in Figure 2. Figure 2 shows the monthly distribution of the changes recorded at the station located before the Dukan reservoir.The figure clearly shows that flow was high in the months of December to May and relatively low in the months of June to November.The highest flow can be seen in the month of July, while the lowest was recorded in the month of September.The month of April showed a flow of around 450 m 3 /sec, while the flow was around 50 m 3 /sec during September.It is also clear from Figure 2 that the months of August to October had more or less the same flow, which was around 50 m 3 /sec.It is also clear from Figure 2 that the flow increased from December to April and then started to decline.This change from increase to decrease was due to the rainfall over the region during that period.The seasonal distribution of the monthly flow recorded after Dukan Dam is shown in Figure 3.The figure shows that flow was highest in the month of August and lowest in April.Overall, most of the months had a flow of more than 100 m 3 /sec.It can be noted that flow started to rise from the month of July and onwards, which is contradictory to the flow before the dam, where the flow started decreasing from the month of July.The increases in flow were due to the surplus water from the reservoir and also in order to maintain a minimum flow in the river.It is also clear from Figures 2  and 3 that the flow was high from the months of January to May upstream and low from June to December on the downstream side of the dam.

Trends in Monthly and Annual Flow
The magnitude of trends in river flow was assessed using Sen's slope while the significance was confirmed using the Mann-Kendall trend test.Sen's slope and the Mann-Kendall trend test were individually applied to each month and the mean annual flow for the period 1964 to 2013.The obtained results for the upstream and downstream side gauges are presented in Table 2. Table 2 shows the Z statistic, significance at the 90% and 95% levels of confidence and the change in each month and mean year.Table 2 shows that there were significant increases and decreases in trends on The seasonal distribution of the monthly flow recorded after Dukan Dam is shown in Figure 3.The figure shows that flow was highest in the month of August and lowest in April.Overall, most of the months had a flow of more than 100 m 3 /sec.It can be noted that flow started to rise from the month of July and onwards, which is contradictory to the flow before the dam, where the flow started decreasing from the month of July.The increases in flow were due to the surplus water from the reservoir and also in order to maintain a minimum flow in the river.It is also clear from Figures 2 and 3 that the flow was high from the months of January to May upstream and low from June to December on the downstream side of the dam.The seasonal distribution of the monthly flow recorded after Dukan Dam is shown in Figure 3.The figure shows that flow was highest in the month of August and lowest in April.Overall, most of the months had a flow of more than 100 m 3 /sec.It can be noted that flow started to rise from the month of July and onwards, which is contradictory to the flow before the dam, where the flow started decreasing from the month of July.The increases in flow were due to the surplus water from the reservoir and also in order to maintain a minimum flow in the river.It is also clear from Figures 2  and 3 that the flow was high from the months of January to May upstream and low from June to December on the downstream side of the dam.

Trends in Monthly and Annual Flow
The magnitude of trends in river flow was assessed using Sen's slope while the significance was confirmed using the Mann-Kendall trend test.Sen's slope and the Mann-Kendall trend test were individually applied to each month and the mean annual flow for the period 1964 to 2013.The obtained results for the upstream and downstream side gauges are presented in Table 2. Table 2 shows the Z statistic, significance at the 90% and 95% levels of confidence and the change in each month and mean year.Table 2 shows that there were significant increases and decreases in trends on

Trends in Monthly and Annual Flow
The magnitude of trends in river flow was assessed using Sen's slope while the significance was confirmed using the Mann-Kendall trend test.Sen's slope and the Mann-Kendall trend test were individually applied to each month and the mean annual flow for the period 1964 to 2013.The obtained results for the upstream and downstream side gauges are presented in Table 2. Table 2 shows the Z statistic, significance at the 90% and 95% levels of confidence and the change in each month and mean year.Table 2 shows that there were significant increases and decreases in trends on both sides of the dam.In the upstream gauge, the months of April, May, November and December showed significant trends.The trends in April, May, November were significant at the 95% level of confidence, while the trends in December were found to be significant at the 99% level of confidence.It is important to mention here that annual flow also had a negative significant trend at the 90% level of confidence.The highest significant trend of −5.08846 m 3 /year was observed for the month of April, followed by −2.70121 m 3 /year in the month of May.On the other hand, the trend in mean annual river flow showed a change of −1.912 m 3 /year at the 90% level of confidence.It also can be noted that none of the month or mean annual flows have shown an increasing trend.A decrease in river flow trends has also been observed in other parts of the world, such as Turkey [18] and Canada [22].
Similar to the upstream side gauge, the downstream side gauge station also showed a decrease in flow in all months except June, where a change of 0.73651 m 3 /sec was recorded.It can be noted that the highest decrease (−4.3828 m 3 /sec) was recorded in the month of September, followed by October (−3.7975 m 3 /sec).It is also important to note that the changes were not significant in both months.A significant change was observed in the months of January, February, March, April, August, November and December.Most of the months were significant at the 90% level of confidence.The trends in the month of August and the mean annual flow were found to be significant at the 95% level of confidence.It is interesting to note that the flow for the month of December significantly decreased on both sides of the dam at the 99% level of confidence.

Trends in Maximum Annual Flow
The trends in maximum flow assessed at the upstream (a) and downstream (b) side gauges are shown in Figure 4.The trends were assessed based on the maximum flow recorded for each month.The figure shows that flow had a decreasing trend on both sides of the dam.The trend estimated using Sen's slope showed a negative change of −4.073 m 3 /sec/year on the upstream side, which was also found to be significant at the 99% confidence level.On the other side of the dam, the flow was also found to have decreased.The decrease rate was recorded as −2.540 m 3 /sec and was significant at the 90% level of confidence.It is important to note that there was a decrease in the maximum flow on the upstream side of the dam, which indicates that the decrease in precipitation resulted in a decrease in the maximum flow.A decrease in precipitation was also reported in the study by Salman, et al. [62] and the study by Shubbar, Salman and Lee [33].Besides the reduction in precipitation, the media also reported that the construction of Garan Dam in Iran, which is near Dukan Dam, has also reduced the flow.Additionally, Hassan, et al. [63] reported that the construction of hydrological schemes in the catchment has reduced the flow in Dukan Dam.Overall, the decrease in maximum flow was relatively higher in recent years on both sides of the dam.

Trends in Minimum Annual Flow
The trends in the minimum flow on both sides of the dam were also assessed.The assessment was performed based on the minimum values recorded each year.The obtained results are displayed in Figure 5.The figure clearly shows that there was a decrease in flow on both sides of the dam.However, the decrease was relatively higher on the upstream side of the dam.The Sen's slope estimated a negative change of −1.082 m 3 /sec, while the Mann-Kendall trend test did not show any significance at the 90%, 95% or 99% confidence level.Similar to the upstream side, the downstream side showed a negative change of −0.948 m 3 /sec.The negative change in the downstream side of the dam was found to be significant at the 95% confidence level.It is also important to note the minimum flow was higher on the downstream side compared with upstream.This could be due to the dependence of flow or discharge from the reservoir.

Trends in Minimum Annual Flow
The trends in the minimum flow on both sides of the dam were also assessed.The assessment was performed based on the minimum values recorded each year.The obtained results are displayed in Figure 5.The figure clearly shows that there was a decrease in flow on both sides of the dam.However, the decrease was relatively higher on the upstream side of the dam.The Sen's slope estimated a negative change of −1.082 m 3 /sec, while the Mann-Kendall trend test did not show any significance at the 90%, 95% or 99% confidence level.Similar to the upstream side, the downstream side showed a negative change of −0.948 m 3 /sec.The negative change in the downstream side of the dam was found to be significant at the 95% confidence level.It is also important to note the minimum flow was higher on the downstream side compared with upstream.This could be due to the dependence of flow or discharge from the reservoir.

Conclusions
This study presented the distributional analysis results, and trends in the monthly, mean annual, annual maximum and annual minimum flow of Dukan Dam located in the northern region of Iraq.The trends were assessed for the period 1964 to 2013.The rate of change in the river flow was assessed using Sen's slope, and the Mann-Kendall trend test was used to confirm the significance of the trends.The results of the study showed that the generalized extreme value trend fit best with the flow data of Dukan Dam.The results of the trend analysis showed that river flow was significantly decreased in all months on both sides of the dam.The month of April showed the highest rate of change upstream, while September showed the highest rate of change on the downstream side of the dam.The month of November and July showed the lowest change in the upstream and downstream sides of the dam, respectively.Significant changes were noticed in the months of April, May, November and December.Similarly, changes in the annual river flow were also found to be significant.Negative significant changes were also noticed in the maximum and minimum flow on both sides of the dam.The maximum flow was found to have decreased relatively significantly on the upstream side of the dam, which indicates a decrease in precipitation.In addition, the construction of hydraulic structures can be a major cause of the reduction in flow.It is expected that the findings of this study will help in the planning and management of the agriculture and water resources of Dukan region.
It can be noted that the results presented in this study are based on the standard Mann-Kendall trend test.In recent years, it has been reported that long-term persistence or the Hurst phenomenon influence the variance of test statistics and can alter the significance of the results acquired using the assumptions of independence and short-term persistence.Therefore, in the future, modified

Conclusions
This study presented the distributional analysis results, and trends in the monthly, mean annual, annual maximum and annual minimum flow of Dukan Dam located in the northern region of Iraq.The trends were assessed for the period 1964 to 2013.The rate of change in the river flow was assessed using Sen's slope, and the Mann-Kendall trend test was used to confirm the significance of the trends.The results of the study showed that the generalized extreme value trend fit best with the flow data of Dukan Dam.The results of the trend analysis showed that river flow was significantly decreased in all months on both sides of the dam.The month of April showed the highest rate of change upstream, while September showed the highest rate of change on the downstream side of the dam.The month of November and July showed the lowest change in the upstream and downstream sides of the dam, respectively.Significant changes were noticed in the months of April, May, November and December.Similarly, changes in the annual river flow were also found to be significant.Negative significant changes were also noticed in the maximum and minimum flow on both sides of the dam.The maximum flow was found to have decreased relatively significantly on the upstream side of the dam, which indicates a decrease in precipitation.In addition, the construction of hydraulic structures can be a major cause of the reduction in flow.It is expected that the findings of this study will help in the planning and management of the agriculture and water resources of Dukan region.
It can be noted that the results presented in this study are based on the standard Mann-Kendall trend test.In recent years, it has been reported that long-term persistence or the Hurst phenomenon influence the variance of test statistics and can alter the significance of the results acquired using the assumptions of independence and short-term persistence.Therefore, in the future, modified versions

Figure 1 .
Figure 1.Location map of the Dukan Dam Reservoir.(a) Iraq map shows the main rivers in Iraq the location of the Dukan Dam Reservoir.(b) A photo taken from Google Earth showing the geographic location, the main rivers feeding the reservoir and its division into two sub-reservoirs (adapted from [43]).

Figure 1 .
Figure 1.Location map of the Dukan Dam Reservoir.(a) Iraq map shows the main rivers in Iraq the location of the Dukan Dam Reservoir.(b) A photo taken from Google Earth showing the geographic location, the main rivers feeding the reservoir and its division into two sub-reservoirs (adapted from [43]).

Figure 2 .
Figure 2. Seasonal distribution of flow recorded before Dukan Dam.

Figure 3 .
Figure 3. Seasonal distribution of flow recorded after Dukan Dam.

Figure 2 .
Figure 2. Seasonal distribution of flow recorded before Dukan Dam.

Figure 2 .
Figure 2. Seasonal distribution of flow recorded before Dukan Dam.

Figure 3 .
Figure 3. Seasonal distribution of flow recorded after Dukan Dam.

Figure 3 .
Figure 3. Seasonal distribution of flow recorded after Dukan Dam.

Hydrology 2019, 6 , 13 Figure 4 .
Figure 4. Trends in the maximum annual river flow on the (a) upstream and (b) downstream sides of Dukan dam.

Figure 4 .
Figure 4. Trends in the maximum annual river flow on the (a) upstream and (b) downstream sides of Dukan dam.

Figure 5 .
Figure 5. Trends in the minimum annual river flow on the (a) upstream and (b) downstream sides of Dukan dam.

Figure 5 .
Figure 5. Trends in the minimum annual river flow on the (a) upstream and (b) downstream sides of Dukan dam.

Table 1 .
Kolmogorov-Smirnov test statistics obtained during the fitting distribution of river flow on the upstream and downstream side gauge data of the dam.

Table 2 .
Trends in monthly and mean annual river flow obtained using Sen's slope and Mann-Kendall trend test located up-and down-stream of Dukan Dam.