Watershed Sediment and Its Effect on Storage Capacity: Case Study of Dokan Dam Reservoir

Dokan is a multipurpose dam located on the Lesser Zab River in the Iraq/Kurdistan region. The dam has operated since 1959, and it drains an area of 11,690 km2. All reservoirs in the world suffer from sediment deposition. It is one of the main problems for reservoir life sustainability. Sustainable reservoir sediment-management practices enable the reservoir to function for a longer period of time by reducing reservoir sedimentation. This study aims to assess the annual runoff and sediment loads of the Dokan Dam watershed using the soil and water assessment tool (SWAT) model to evaluate the relative contributions in comparison with the total values delivered from both watershed and Lesser Zab River and to identify the basins with a high sediment load per unit area. These help in the process of developing a plan and strategy to manage sediment inflow and deposition. The SUFI-2 program was applied for a model calibrated based on the available field measurements of the adjacent Derbendekhan Dam watershed, which has similar geological formations, characteristics and weather. For the calibration period (1961–1968), the considered statistical criteria of determination coefficients and Nash–Sutcliffe model efficiency were 0.75 and 0.64 for runoff while the coefficients were 0.65 and 0.63 for sediment load, respectively. The regionalization technique for parameter transformation from Derbendekhan to Dokan watershed was applied. Furthermore, the model was validated based on transformed parameters and the available observed flow at the Dokan watershed for the period (1961–1964); they gave reasonable results for the determination coefficients and Nash–Sutcliffe model efficiency, which were 0.68 and 0.64, respectively. The results of SWAT project simulation for Dokan watershed for the period (1959–2014) indicated that the average annual runoff volume which entered the reservoir was about 2100 million cubic meters (MCM). The total sediment delivered to the reservoir was about 72 MCM over the 56 years of dam life, which is equivalent to 10% of the reservoir dead storage. Two regression formulas were presented to correlate the annual runoff volume and sediment load with annual rain depth for the studied area. In addition, a spatial distribution of average annual sediment load was constructed to identify the sub basin of the high contribution of sediment load.


Introduction
Most of the dams, storage schemes, and different hydraulic structures around the world suffer from sedimentation problems.For dams and reservoirs, this effect is mainly concerned with the design capacity and operation schedule.The main source of reservoir sediments is the main river flow in and incorporating both sensitivity and uncertainty analysis.A number of previous studies [14][15][16][17] were applied using the SWAT model to estimate runoff, sediment yield and/or other soluble materials for ungagged watersheds based on neighboring or similar property watersheds.Another technique is used for flow and erosion, and sediment transport is the distributed mode.Juez et al. [18] simulates the hydro-sedimentary response of the Western Mediterranean catchment to a representative rainfall storm.The simulation combines the distributed flow surface model with the empirical model for infiltration, the Soil Conservation Services Model (SCS), and the erosion model, which is the Hillslope Erosion Model (HEM), considering water depth and flow as a 2D model.The present model is a tool for analyzing the hydro-sedimentary process at a temporal and special scale.
Most countries in the Middle East suffer from water shortage problems, where the annual allocation per capita does not exceed 500 m 3 [19].For this reason, water is essential to life, socioeconomic development, and political stability in this region.Future prospects are negative; therefore, this problem is expected to be more chronic and severe in future [20].Iraq used to be considered to be a relatively rich country in terms of its water resources, until the mid-1970s, because of the presence of the Tigris and Euphrates Rivers [21].Due to regional and internal problems in Iraq, the estimation of the overall water required is about 75-81 billion cubic meters (BCM) [22], while the available quantity is 59-75 BCM and will drop to 17.61 BCM in 2025 [23].In view of this situation, it is very important to know the actual storage capacity of the reservoirs-which are unknown now-so that prudent water resources planning can be done.The sedimentation rate of several reservoirs was recently investigated in Mosul and Dohuk.This is the third reservoir to be dealt with in Iraq.The bed of the Dokan Dam reservoir (located in the northeast of Iraq) is surveyed by Hassan et al. [24], and this studied the bed sediment using 32-bed samples distributed spatially over the reservoir.The results indicate that the bed sediments of the reservoir are composed of silt (48%), clay (23%), gravel (15%) and sand (14%).
All reservoirs in the world suffer from sediment deposition.This is one of the main problems for reservoir life sustainability.Sustainable reservoir sediment-management practices enable continued reservoir functioning for a longer period by reducing reservoir sedimentation.Iraq suffers from water shortage problems, especially after the construction of a series of storage reservoirs in source countries (Turkey, Syria and Iran), so the evaluation of the actual live storage capacity of dams is important for the prudent management of the operation schedule.The aim of this study is to assess the annual runoff and sediment loads of the Dokan Dam watershed (ungagged area) using a SWAT model set-up based on the parameter transformation technique of the modeling-gauged Derbendekhan watershed to learn the hydrological behavior of the area and to assess its contribution to the total values pouring into the reservoir.Moreover, the set-up model helps us to find the spatial distribution of erosion and annual sediment yield for the sub basins.This will help us to find the sub basins with a high sediment yield and evaluate effective factors for them.These assessments help in the process of developing a plan and a strategy for managing the sediment inflow and deposition.

Location and Topography
The considered study area is the watershed of the Dokan Dam reservoir, situated in the northeast of Iraq (Figure 1).Dokan dam is a concrete arch dam located in the Lower Zab River, about 65 km southeast of Sulaimaniyah city and 295 km north of Baghdad, the capital city of Iraq.The dam height is about 116 m at maximum river depth, having a total storage capacity of 6.87 × 10 9 m 3 (6.14 × 10 9 m 3 live storage and 0.73 × 10 9 m 3 dead storage) at normal operation level of 511 m.a.s.l.[22].The dam has been built to serve irrigation, power generation, water supply and flood control needs.Due to the limited observed data of flow at the Dokan Dam watershed and the unavailability of sediments load data, the second watershed considered for this research is the Derbendekhan Dam watershed; it is the nearest watershed to the study area.The properties of the Dokan and Derbendikhan watersheds are shown in Table 1.The digital elevation model (DEM) with 30 m resolution is considered to identify the watershed boundary, classification of overland and channel flow, slopes and other properties.

Soil Type and Land Use
The exposed rocks at the Dokan and Derbendekhan watershed areas are mainly limestone and minor exposures of dolomitic limestone, dolomite, and Quaternary alluvial deposits [25,26].Based on the Reconnaissance Soil Map of the three Northern Governorates, Iraq [26] and the Food and Agriculture Organization of the United Nations (FAO) soil map [27], both watersheds are located on a common extended type of soil classification.Samples for different soil classes were taken depending on the soil map of the study area.A map of soil types is prepared for this study as a shape file for each watershed to be used in the SWAT model (see below for model details).The soil samples analysis includes grain size distribution in different types, organic matter content and hydraulic conductivity.The analysis of soil samples indicates that the area generally consists of four major soil types.Most of the area (85.6%) is covered by gravelly sandy mud; 6.9% is gravelly mud; and 7.5% of two types of muddy gravel (the main differences between the two types are the percent of gravel, which is 74% and 56% for types 1 and 2, respectively).Figure 2 shows the shape file of soil type considered in the SWAT model for the Dokan watershed.
The land use map for the years (1976-1979) [28] and available satellite image (NASA's Landsat GeoCover, 2007, with a spatial resolution of 14.25 m) indicates that the winter plants and pastures represent the main part of the land use map of the studied area.This depends on rainfall as a main source of irrigation water.The other parts are forests, vegetables and urban areas (villages).The land use change is limited (see Table 2).This is mainly due to the geological nature of both watersheds.In addition, the topography of the area does not enhance any changes.It is noteworthy to mention that rain is the main source of irrigation.For these reasons, the land cover did not change widely through the study period since the operation of the dam from the year 1959 to the year 2014.The Dam and

Soil Type and Land Use
The exposed rocks at the Dokan and Derbendekhan watershed areas are mainly limestone and minor exposures of dolomitic limestone, dolomite, and Quaternary alluvial deposits [25,26].Based on the Reconnaissance Soil Map of the three Northern Governorates, Iraq [26] and the Food and Agriculture Organization of the United Nations (FAO) soil map [27], both watersheds are located on a common extended type of soil classification.Samples for different soil classes were taken depending on the soil map of the study area.A map of soil types is prepared for this study as a shape file for each watershed to be used in the SWAT model (see below for model details).The soil samples analysis includes grain size distribution in different types, organic matter content and hydraulic conductivity.The analysis of soil samples indicates that the area generally consists of four major soil types.Most of the area (85.6%) is covered by gravelly sandy mud; 6.9% is gravelly mud; and 7.5% of two types of muddy gravel (the main differences between the two types are the percent of gravel, which is 74% and 56% for types 1 and 2, respectively).Figure 2 shows the shape file of soil type considered in the SWAT model for the Dokan watershed.
The land use map for the years (1976-1979) [28] and available satellite image (NASA's Landsat GeoCover, 2007, with a spatial resolution of 14.25 m) indicates that the winter plants and pastures represent the main part of the land use map of the studied area.This depends on rainfall as a main source of irrigation water.The other parts are forests, vegetables and urban areas (villages).The land use change is limited (see Table 2).This is mainly due to the geological nature of both watersheds.In addition, the topography of the area does not enhance any changes.It is noteworthy to mention that Water 2018, 10, 858 5 of 16 rain is the main source of irrigation.For these reasons, the land cover did not change widely through the study period since the operation of the dam from the year 1959 to the year 2014.The Dam and the studied watershed are located away from the main cities, so changes in the urban and rural areas are limited.Table 2 shows the percentage of different land use cover for the two periods of the available land use map.
Water 2018, 10, x FOR PEER REVIEW 5 of 16 the studied watershed are located away from the main cities, so changes in the urban and rural areas are limited.Table 2 shows the percentage of different land use cover for the two periods of the available land use map.Due to the small difference between the percentage of land cover for the two available years, a map of land use for the study area is prepared as a shape file (Figure 3).The area consists of four types of land use/cover.Winter plants (pasture) and forests of different types of trees cover the main part of the study area, while the remaining small area is planted with vegetables near the reservoir boundary and/or in urban areas (villages).Due to the small difference between the percentage of land cover for the two available years, a map of land use for the study area is prepared as a shape file (Figure 3).The area consists of four types of land use/cover.Winter plants (pasture) and forests of different types of trees cover the main part of the study area, while the remaining small area is planted with vegetables near the reservoir boundary and/or in urban areas (villages).the studied watershed are located away from the main cities, so changes in the urban and rural areas are limited.Table 2 shows the percentage of different land use cover for the two periods of the available land use map.Due to the small difference between the percentage of land cover for the two available years, a map of land use for the study area is prepared as a shape file (Figure 3).The area consists of four types of land use/cover.Winter plants (pasture) and forests of different types of trees cover the main part of the study area, while the remaining small area is planted with vegetables near the reservoir boundary and/or in urban areas (villages).

Applied Model
The soil and water assessment tool (SWAT) is a physically based continuous simulation model for short or long times that can be applied to large river basins and complex watersheds.It was developed by the US Department of Agriculture, Agricultural Research Service [6].The model is an efficient tool to estimate the flow and sediment load in addition to different chemical and nutrient materials.The model divides the watershed into sub basins based on their DEM date and hydrological response units (HRU); each unit has the same soil type, land use and land slope.
The hydrological simulation is based on the topographical terrain, soil type, land use and hydrological data of daily precipitation, maximum and minimum temperature, wind speed, relative humidity and solar radiation.The flow can be estimated based on a water balance equation; this equation is simulated in the SWAT model by different modular: the land phase and routing phase [14].For the land phase, the soil water balance calculation is based on the following form [29]: where, SW t : Water content of the soil (mm); SW 0 : Initial water content (mm); R day : Depth of precipitation (mm); Q sur : Equivalent depth of surface runoff (mm); E a : Evapotranspiration depth (mm); w seep : Depth of water seepage out of considered surface profile (mm); Q qw : Equivalent depth of return flow (mm).
The Penman-Monteith method is considered for potential evapotranspiration estimation.The required input data to estimate the potential evapotranspiration (PET) using the Penman-Monteith method are daily solar radiation, air temperature, relative humidity and wind speed.The formula of this method considers three effective factors for evapotranspiration, which are the required energy to sustain evaporation, the required strength to remove the water vapor and the aerodynamic in addition to resistance of the surface.The Penman-Monteith method is in the following form [29]: Also, the different parameters of land management are recognized based on soil type, land use and land cover.The soil water content and soil infiltration can be estimated by two methods based on the available data, either by the Green-Ampt infiltration equation or curve number methods.The Green-Ampt equation requires rainfall data of a sub daily interval, which is not available in Iraqi weather stations, so the curve number method is utilized throughout this work using the following form [29]: where, Q sur f : Equivalent depth of surface runoff (mm); R day : Rainfall depth of the considered day (mm); S: Retention parameter (mm).
The value of S can be estimated by the following equation [29]: where, CN is the curve number of that considered day.
The second process includes the estimation of soil erosion from the overland due to rainfall detachment and surface runoff in addition to channel erosion and deposition.The sediments, routing in both the overland and channel flow, are estimated based on rainfall data, soil properties, land use/land cover and topography.The maps of soil type and land use are required with the digital elevation model (DEM) data to identify the topography of the watershed and to classify it into overland and channel sediment flow.The modified universal soil loss equation (MUSLE) is considered in the following form [29]: where, sed: Yield of sediment for the considered storm or day (ton.);Q sur f : Volume of surface runoff (mm/ha); q peak : Greatest surface runoff rate (m

Runoff and Sediment Load Calibration
Although the mathematical and conceptual models are considered widely in hydrological studies to simulate different events, such as runoff flow, sediment and both suspended and dissolved material transport, they still require calibration with measured values to ensure the accuracy of the model outputs.
Two types of dataset are prepared to be applied to SWAT model.The metrological (climate) data include daily precipitation, maximum and minimum temperature, wind speed, relative humidity and solar radiation; hydrometric data are also present.The second dataset is the topography data, which includes the DEM map.In view of the limited available measurements of flow and unavailable sediment records of the Dokan Dam watershed, the available data of the Derbendekhan Dam watershed (the adjacent watershed to Dokan, Figure 1), were considered for both flow and sediment model calibration.Due to the similarity in the geological formation, soil type, land use, topographical and watershed characteristics (Table 1), the parameters of the calibrated watershed can be transformed to an ungauged watershed model [17].To calibrate the results of the SWAT model for both runoff and sediment, the SWAT-CUP software is applied.It is an efficient tool to adjust different parameters of the SWAT model to obtain optimal local results and create an uncertainty analysis of SWAT model parameters to provide an easy and quick method of calibration and standardized calibration [30].The considered software for the model is the Sequential Uncertainly Fitting version 2 (SUFI-2, Swiss Federal Institute of Aquatic Science and Technology, Zurich, Switzerland).In this program, all the uncertainty parameters can be used in the model calibration, including uncertainly in driving variables, parameters of the conceptual model and considered data [30].Different statistical criteria can be considered in the model objective function to evaluate model performance, such as the determination coefficient, Nash-Sutcliff model efficiency, root mean square error and Chi-square.The determination coefficient is considered to be effective criteria to obtain the optimal values of flow and sediment concentration between the observed and measured data.
The SWAT project for the Derbendekhan watershed is set-up including the required DEM data, soil type land use, as shown in Figures 1-3, respectively, and meteorological data based on the nearest stations to the area as shown in Figure 1.The monthly average flow rate data at Derbendekhan station are considered for model calibration.To obtain an enhanced calibration of the model and for more understanding of the model's performance, the monthly recorded flow data are separated into base flow and surface runoff [14].The recursive digital filter technique [31] is used to obtain a monthly separation based on the original daily separation technique.The separated monthly runoff from the total flow as measured values is applied in SUFI-2 to calibrate the model parameters.The statistical criterion of the determination coefficient is used as an objective function criterion.Besides this, the Nash-Sutcliffe model is employed to evaluate the model performance.For monthly runoff flow, the highest obtained values are 0.75 and 0.64 for the coefficient of determination and Nash-Sutcliffe model, respectively.The model uncertainty was measured using two factors: the P-factor reflects the percentage of measured data bracketed by 95% prediction uncertainty (95PPU).This means that one minus the P-factor represents the presence of poorly simulated values.The R-factor is another measurement of model uncertainly equal to the average thickness of the 95PPU band divided by the standard deviation of measured data.For the runoff calibration period, the P-factor is 0.78 and the R-factor is 1.27.Figure 4 shows the observed and simulated runoff at the Derbendekhan watershed outlet and the uncertainty band (95PPU) for the period (1961)(1962)(1963)(1964)(1965)(1966)(1967)(1968).
Water 2018, 10, x FOR PEER REVIEW 8 of 16 ungauged watershed model [17].To calibrate the results of the SWAT model for both runoff and sediment, the SWAT-CUP software is applied.It is an efficient tool to adjust different parameters of the SWAT model to obtain optimal local results and create an uncertainty analysis of SWAT model parameters to provide an easy and quick method of calibration and standardized calibration [30].The considered software for the model is the Sequential Uncertainly Fitting version 2 (SUFI-2, Swiss Federal Institute of Aquatic Science and Technology, Zurich, Switzerland).In this program, all the uncertainty parameters can be used in the model calibration, including uncertainly in driving variables, parameters of the conceptual model and considered data [30].Different statistical criteria can be considered in the model objective function to evaluate model performance, such as the determination coefficient, Nash-Sutcliff model efficiency, root mean square error and Chi-square.
The determination coefficient is considered to be effective criteria to obtain the optimal values of flow and sediment concentration between the observed and measured data.
The SWAT project for the Derbendekhan watershed is set-up including the required DEM data, soil type land use, as shown in Figures 1-3, respectively, and meteorological data based on the nearest stations to the area as shown in Figure 1.The monthly average flow rate data at Derbendekhan station are considered for model calibration.To obtain an enhanced calibration of the model and for more understanding of the model's performance, the monthly recorded flow data are separated into base flow and surface runoff [14].The recursive digital filter technique [31] is used to obtain a monthly separation based on the original daily separation technique.The separated monthly runoff from the total flow as measured values is applied in SUFI-2 to calibrate the model parameters.The statistical criterion of the determination coefficient is used as an objective function criterion.Besides this, the Nash-Sutcliffe model is employed to evaluate the model performance.For monthly runoff flow, the highest obtained values are 0.75 and 0.64 for the coefficient of determination and Nash-Sutcliffe model, respectively.The model uncertainty was measured using two factors: the P-factor reflects the percentage of measured data bracketed by 95% prediction uncertainty (95PPU).This means that one minus the P-factor represents the presence of poorly simulated values.The R-factor is another measurement of model uncertainly equal to the average thickness of the 95PPU band divided by the standard deviation of measured data.For the runoff calibration period, the P-factor is 0.78 and the Rfactor is 1.27.Figure 4 shows the observed and simulated runoff at the Derbendekhan watershed outlet and the uncertainty band (95PPU) for the period (1961-1968).The available measured sediment concentration data of the Diyala River at Derbendekhan station, as presented by Assad [32], were utilized to calibrate the SWAT model parameters for the same period of runoff flow calibration.Figure 5 shows the measured and optimal simulated values of sediment load concentration at Derbendekhan outlet.The same statistical criteria are implemented, and the optimal resultant values are 0.65 and 0.63 for determination coefficient and Nash-Sutcliffe model efficiency, respectively, while the P-factor and R-factor equal to 0.68 and 2.27 respectively.
Water 2018, 10, x FOR PEER REVIEW 9 of 16 The available measured sediment concentration data of the Diyala River at Derbendekhan station, as presented by Assad [32], were utilized to calibrate the SWAT model parameters for the same period of runoff flow calibration.Figure 5 shows the measured and optimal simulated values of sediment load concentration at Derbendekhan outlet.The same statistical criteria are implemented, and the optimal resultant values are 0.65 and 0.63 for determination coefficient and Nash-Sutcliffe model efficiency, respectively, while the P-factor and R-factor equal to 0.68 and 2.27 respectively.The most effective parameters are selected for the runoff and sediment model calibration as proposed by [28] in addition to other parameters.The resultant best-fitted values (optimal values) of the parameters and the considered range are shown in Table 3. Parameters listed in Table 3 are also mentioned in other previous literature [14,33].The most effective parameters are selected for the runoff and sediment model calibration as proposed by [28] in addition to other parameters.The resultant best-fitted values (optimal values) of the parameters and the considered range are shown in Table 3. Parameters listed in Table 3 are also mentioned in other previous literature [14,33].

Runoff Validation
A SWAT model project is set-up for the Dokan Dam watershed, the main part of the study area.The limited recorded data of the monthly flow rate at the outlet of the Dokan watershed and the absence of any sediment load measurements leads to utilizing the regionalization technique to transfer the effective parameters from the adjacent gauged (Derbendekhan) watershed.Due to its similarity in geological formation, soil type, land use, watershed characteristics and weather data, the effective hydrological parameters obtained from the Derbendekhan (gauged) watershed can be transformed to the Dokan (ungauged) watershed.The process of parameter transformation is called regionalization.There are a number of presented methods for the regionalization of the watershed hydrological parameters: Kokkonen [34] applies the regression approach, while Parajka et al. [35] employs kringing and a similarity approach and Heuvelmans et al. [36] investigates the application of artificial neural nets and other methods.Since the Derbendekhan watershed is adjacent to the Dokan watershed (Figure 1), and the physical, topographical properties and rainfall are similar (Table 1) along with the geological formation, land use/land cover, and soil type, both watersheds have similar flow and sediment parameters.In this case, the effective parameters can be transformed from a donor watershed to an ungauged watershed.The fitted values of the Derbendekhan parameters calibrated by the SUFI-2 program are transferred to the Dokan watershed SWAT project.The SUFI-2 program is implemented for the calibration, uncertainly analysis and regionalization of the considered parameters of the SWAT model for the Debendekhan watershed for runoff and sediment load of concentration.Here, the SUFI-2 algorithm is used for the calibration, validation and measurement of the uncertainty for input data, model and sensitive parameters.The degree of uncertainty is measured by two values: P-factor and R-factor.The percent of measured values bracketed by 95% prediction uncertainly represent (95PPU), which is the P-factor while the ratio of 95PPU thickness divided by standard deviation of measured values is equal to the R-factor.When the simulated values are exactly the measured ones, the value of P-factor equals 1 and the value of R-factor equals zero [37].
Based on the transformed parameters, the simulated runoff flows are compared with measured values at the Dokan watershed outlet after the separation of the base flow for the period 1961-1964 (Figure 6) to evaluate the effectiveness of the transformation of the hydrological parameters.For the validation period of runoff, the obtained values of the determination coefficient and Nash-Sutcliffe model efficiency are 0.68 and 0.64, respectively, indicating that the transformation process is successful.The P-factor and R-factor for the model uncertainty also indicate a reasonable model performance: the P-factor is 0.71 and R-factor is 0.97.
Water 2018, 10, x FOR PEER REVIEW 10 of 16

Runoff Validation
A SWAT model project is set-up for the Dokan Dam watershed, the main part of the study area.The limited recorded data of the monthly flow rate at the outlet of the Dokan watershed and the absence of any sediment load measurements leads to utilizing the regionalization technique to transfer the effective parameters from the adjacent gauged (Derbendekhan) watershed.Due to its similarity in geological formation, soil type, land use, watershed characteristics and weather data, the effective hydrological parameters obtained from the Derbendekhan (gauged) watershed can be transformed to the Dokan (ungauged) watershed.The process of parameter transformation is called regionalization [14].There are a number of presented methods for the regionalization of the watershed hydrological parameters: Kokkonen [34] applies the regression approach, while Parajka et al. [35] employs kringing and a similarity approach and Heuvelmans et al. [36] investigates the application of artificial neural nets and other methods.Since the Derbendekhan watershed is adjacent to the Dokan watershed (Figure 1), and the physical, topographical properties and rainfall are similar (Table 1) along with the geological formation, land use/land cover, and soil type, both watersheds have similar flow and sediment parameters.In this case, the effective parameters can be transformed from a donor watershed to an ungauged watershed.The fitted values of the Derbendekhan parameters calibrated by the SUFI-2 program are transferred to the Dokan watershed SWAT project.The SUFI-2 program is implemented for the calibration, uncertainly analysis and regionalization of the considered parameters of the SWAT model for the Debendekhan watershed for runoff and sediment load of concentration.Here, the SUFI-2 algorithm is used for the calibration, validation and measurement of the uncertainty for input data, model and sensitive parameters.The degree of uncertainty is measured by two values: P-factor and R-factor.The percent of measured values bracketed by 95% prediction uncertainly represent (95PPU), which is the P-factor while the ratio of 95PPU thickness divided by standard deviation of measured values is equal to the R-factor.When the simulated values are exactly the measured ones, the value of P-factor equals 1 and the value of Rfactor equals zero [37].
Based on the transformed parameters, the simulated runoff flows are compared with measured values at the Dokan watershed outlet after the separation of the base flow for the period 1961-1964 (Figure 6) to evaluate the effectiveness of the transformation of the hydrological parameters.For the validation period of runoff, the obtained values of the determination coefficient and Nash-Sutcliffe model efficiency are 0.68 and 0.64, respectively, indicating that the transformation process is successful.The P-factor and R-factor for the model uncertainty also indicate a reasonable model performance: the P-factor is 0.71 and R-factor is 0.97.

Results and Discussion
After Durbendekhan SWAT project calibration for runoff and sediment data, the best-fitted values of the hydrological watershed parameters are transformed by the regionalization technique to the Dokan SWAT project, the main project of the study.The sensitivity analysis is also studied for the most effective parameters on both the runoff and sediment load.The sensitivities are accomplished to identify the effective parameters on runoff and sediment load values for the watershed.The parameter sensitivity is estimated in the SUFI-2 model based on the multiple regression system presented by Abbaspour [38] to evaluate the effect of the considered parameter value (b i ) on the objective function (g); its sensitivity is in the following form: This formula calculates the average changes in the objective function due to the change in a given parameter while other parameters are changing.The comparative significance and sensitivity of each parameter are estimated based on the statistical criteria of the t-stat and p-value.The t-stat value is obtained from the coefficient of a parameter in the multiple regression analysis divided by its standard error.If the coefficient value is large in comparison to the standard errors, this mean that the parameter is sensitive.The p-value can be obtained by comparing the t-stat value with the student's distribution table [37].The p-value of each term test is the null hypothesis, in which the coefficient is not affected.If the p-value is less than 0.05, it indicates that the null hypothesis can be rejected.The t-stat and p-values of different effective parameters are shown in Table 4.The parameters are arranged from low to high sensitivity, i.e., from low t-stat value or high p-value.The result of the test indicates that the soil water content, soil curve number at normal conditions (CN2) and the soil saturated hydraulic conductivity are the most effective parameters while threshold depth in a shallow aquifer, the channel erodibility factor and soil erodibility (k) factor in USLE have the lowest effect on runoff and sediment load simulation.The model is applied for the study period to estimate the runoff and sediment load reaching the Dokan Reservoir.The considered years of simulation have begun since the operation of the dam in 1959 to the year of the bathometry survey (2014) carried out by [39].The resultant annual runoff volume that enters the Dukan Reservoir from the HRUs ranges from 300 to 4600 MCM (Figure 7), depending on rainfall intensity, depth and distribution through the rainy season.The runoff average volume from the watershed represents 35% of the live storage capacity of the dam, indicating that watershed runoff makes a significant contribution to reservoir inflow.
The SWAT model is an efficient tool to estimate the runoff hydrograph, but, in some hydrological studies, such as scheduling reservoir operations to supply the demand rate, the assessment of water resource income only is required.A regression formula is determined based upon the input-output data of SWAT model of the study area.This is a simple and quick tool to correlate the annual runoff depth with annual rainfall with good correlation (R 2 = 0.9) results without the need for more detailed input data as required in SWAT projects.The relationship used is in the following form: Ann − 1.92 R 2 = 0.90 ( f or R Ann > 16 mm) (7) where Run Ann : Annual runoff depth mm); R Ann : Annual rainfall depth (mm).
Water 2018, 10, x FOR PEER REVIEW 12 of 16 volume from the watershed represents 35% of the live storage capacity of the dam, indicating that watershed runoff makes a significant contribution to reservoir inflow.The SWAT model is an efficient tool to estimate the runoff hydrograph, but, in some hydrological studies, such as scheduling reservoir operations to supply the demand rate, the assessment of water resource income only is required.A regression formula is determined based upon the input-output data of SWAT model of the study area.This is a simple and quick tool to correlate the annual runoff depth with annual rainfall with good correlation (R 2 = 0.9) results without the need for more detailed input data as required in SWAT projects.The relationship used is in the following form: where : Annual runoff depth mm); : Annual rainfall depth (mm).The formula is suitable when the annual runoff depth is greater than 15 mm to avoid negative runoff; this value is already much lower than the minimum historical recorded value.
The sediment load delivered to the Dokan reservoir is also estimated based on MUSLE programmed into the SWAT model for each single storm.The results are presented here annually.The average annual sediment load concentration is 650 mg/ℓ.This concentration can be considered relatively low in comparison with other locations or measurements in the region as proposed by [40,41] as well as the worldwide rate [42].This is due to the nature of the rocks of the area and the effect of plant cover, such as winter pasture and plants and some forest trees throughout the region which reduce the detachment of soil particles transported with runoff flow.The estimated annual sediment load delivered to the Dokan reservoir from the watershed ranges from 3.6 to 0.16 × 10 6 ton for studied period, Figure 7.The average annual value is 1.63 × 10 6 ton.
The sediment trap efficiency of the reservoir is estimated based on the method presented by Garg V. and Jothiprakash V. [43].This depends on reservoir storage capacity and annual inflow.The trap efficiency of the Dokan reservoir changes through the study period from 1 to 0.985.
Based on the results obtained from the simulation model, the estimated sediment load volume deposited in the reservoir for the considered period is about 10% of the dead storage capacity.This value is for the watershed only, which is considered a reasonable value and does not affect the The formula is suitable when the annual runoff depth is greater than 15 mm to avoid negative runoff; this value is already much lower than the minimum historical recorded value.
The sediment load delivered to the Dokan reservoir is also estimated based on MUSLE programmed into the SWAT model for each single storm.The results are presented here annually.The average annual sediment load concentration is 650 mg/ .This concentration can be considered relatively low in comparison with other locations or measurements in the region as proposed by [40,41] as well as the worldwide rate [42].This is due to the nature of the rocks of the area and the effect of plant cover, such as winter pasture and plants and some forest trees throughout the region which reduce the detachment of soil particles transported with runoff flow.The estimated annual sediment load delivered to the Dokan reservoir from the watershed ranges from 3.6 to 0.16 × 10 6 ton for studied period, Figure 7.The average annual value is 1.63 × 10 6 ton.
The sediment trap efficiency of the reservoir is estimated based on the method presented by Garg V. and Jothiprakash V. [43].This depends on reservoir storage capacity and annual inflow.The trap efficiency of the Dokan reservoir changes through the study period from 1 to 0.985.
Based on the results obtained from the simulation model, the estimated sediment load volume deposited in the reservoir for the considered period is about 10% of the dead storage capacity.This value is for the watershed only, which is considered a reasonable value and does not affect the designed project life.However, the Lesser Zab load should be also considered to evaluate the total amount of sediment load delivered and deposited within the reservoir.The total amount of sediment deposited in the Dokan reservoir for the period (1959-2014) is 209 MCM [38].This means that the sediment load delivered into the reservoir from the watershed based on simulated results is about 34% of the total sediments deposited within the reservoir, which mean that the watershed sediment contribution is an effective value.
Due to the huge amount of data required, including topographical, metrological and hydrological data, different maps of soil, and land use to estimate the sediment load based on the applied model, a simple regression formula is used.It is based on simulated values to correlate the sediment load per unit area of the Dokan watershed with annual runoff depth in the following form: where Sed Ann : Annual sediment load per unit area, (ton/km 2 ); R Ann : Annual rainfall depth (mm).
The formula is suitable when the annual runoff depth is greater than 60 mm; this value is already lower than the minimum historical recorded value in Dokan area.
A special distribution map of average annual sediment yield per unit area of sub basins is also prepared (Figure 8a).It can be noticed that the annual sediment load contribution ranges from 13 to 950 ton/km 2 approximately.The rate of erosion and sediment yield depends on a number of factors: topography, soil type, land cover and rainfall intensity.Comparing the sub basins of different soil types and land uses, in both the annual sediment yield map (Figure 8a) and the sub basin slope map (Figure 8b), it can be noticed that the most effective factor for the sediment yield is the land slope rather than other factors.This can be clearly noticed in that some basins have the same soil type and land use but a higher slope gives higher sediment yield.Sub basins having an average slope between 25 to 45% represent the area of high sediment yields from 400 to 950 t/km 2 ; however, when the sub basin slope is less than 20%, the sediment yield per unit area is reduced to about 40 t/km 2 .
Water 2018, 10, x FOR PEER REVIEW 13 of 16 designed project life.However, the Lesser Zab load should be also considered to evaluate the total amount of sediment load delivered and deposited within the reservoir.The total amount of sediment deposited in the Dokan reservoir for the period (1959-2014) is 209 MCM [38].This means that the sediment load delivered into the reservoir from the watershed based on simulated results is about 34% of the total sediments deposited within the reservoir, which mean that the watershed sediment contribution is an effective value.Due to the huge amount of data required, including topographical, metrological and hydrological data, different maps of soil, and land use to estimate the sediment load based on the applied model, a simple regression formula is used.It is based on simulated values to correlate the sediment load per unit area of the Dokan watershed with annual runoff depth in the following form: 0.056. .7.33 R 0.97 60 (8) where : Annual sediment load per unit area, (ton/km 2 ); : Annual rainfall depth (mm).
The formula is suitable when the annual runoff depth is greater than 60 mm; this value is already lower than the minimum historical recorded value in Dokan area.
A special distribution map of average annual sediment yield per unit area of sub basins is also prepared (Figure 8a).It can be noticed that the annual sediment load contribution ranges from 13 to 950 ton/km 2 approximately.The rate of erosion and sediment yield depends on a number of factors: topography, soil type, land cover and rainfall intensity.Comparing the sub basins of different soil types and land uses, in both the annual sediment yield map (Figure 8a) and the sub basin slope map (Figure 8b), it can be noticed that the most effective factor for the sediment yield is the land slope rather than other factors.This can be clearly noticed in that some basins have the same soil type and land use but a higher slope gives higher sediment yield.Sub basins having an average slope between 25 to 45% represent the area of high sediment yields from 400 to 950 t/km 2 ; however, when the sub basin slope is less than 20%, the sediment yield per unit area is reduced to about 40 t/km 2 .
This map is a tool that can enable decision-makers to apply a suitable method to reduce the erosion load, especially from high erosion rate areas.Depending on the selected area, the treatment may include practicing strip planting, terracing, or contour forming to reduce the effect of slope on surface runoff flow velocity, erosion and sediment transport capacity.This map is a tool that can enable decision-makers to apply a suitable method to reduce the erosion load, especially from high erosion rate areas.Depending on the selected area, the treatment may include practicing strip planting, terracing, or contour forming to reduce the effect of slope on surface runoff flow velocity, erosion and sediment transport capacity.

Conclusions
The soil and water assessment tool (SWAT) model is applied to assess the runoff and sediment delivered from the Dokan Dam watershed.Due to the limited recorded flow data and the absence of sediment measurements data at a station near the inlet of Dokan Reservoir, the model is calibrated for both runoff and sediment load for the Derbendekhan watershed adjacent to the Dokan watershed.The regionalization technique is employed to transfer the calibrated parameters of the SWAT project from a gauged (Derbendekhan) to an ungauged (Dokan) watershed.The resultant monthly runoff flow for the Dokan SWAT project is based on transformed parameters which were compared with measured values to evaluate the regionalization technique and model performance.The determination coefficient (R 2 ) and Nash-Sutcliffe model efficiency (Eff.) are 0.68 and 0.64, respectively, indicating a reasonable model performance with this technique.The average watershed contribution for annual runoff represents 35% of the dam life storage; this percentage is considered effective in the dam operation schedule.The total sediment load delivered to the Dokan reservoir from the watershed for the studied period is about 72 MCM.This load forms about 10% of the dead storage capacity of the reservoir.Generally, the total sediment load delivered and deposited in the reservoir for the period of dam operation is considered acceptable within the allowed limits.The map of special distribution of annual sediment load yield per unit area of each sub basins is presented; the average slopes map reflects a good agreement with the map of annual sediment load yield in comparison to other effective factors to be considered.This indicates that the land slope is the most effective factor on erosion and sediment transport.This can be used for soil conservation treatment to reduce the erosion rate.

Figure 1 .
Figure 1.Topographic map of the watershed areas of the Dokan and Derbendekhan Dams and their locations in Iraq.

Figure 2 .
Figure 2. Soil type classification of the Dokan Reservoir watershed.

Figure 3 .
Figure 3. Land use and land cover classification of the Dokan Reservoir watershed.

Figure 2 .
Figure 2. Soil type classification of the Dokan Reservoir watershed.

Water 2018 ,
10, x FOR PEER REVIEW 5 of 16

Figure 2 .
Figure 2. Soil type classification of the Dokan Reservoir watershed.

Figure 3 .
Figure 3. Land use and land cover classification of the Dokan Reservoir watershed.

Figure 3 .
Figure 3. Land use and land cover classification of the Dokan Reservoir watershed.

Figure 4 .
Figure 4. Monthly observed and simulated runoff at the Derbendekhan outlet and 95PPU for the period 1961-1968.

Figure 4 .
Figure 4. Monthly observed and simulated runoff at the Derbendekhan outlet and 95PPU for the period 1961-1968.

Figure 6 .
Figure 6.Monthly observed and simulated runoff at the Dokan outlet and 95PPU for the period 1962-1965.

Figure 6 .
Figure 6.Monthly observed and simulated runoff at the Dokan outlet and 95PPU for the period 1962-1965.

Figure 7 .
Figure 7. Annual runoff volume and sediment load delivered to the Dokan Reservoir for the period 1959-2014.

Figure 7 .
Figure 7. Annual runoff volume and sediment load delivered to the Dokan Reservoir for the period 1959-2014.

Figure 8 .
Figure 8.(a) Spatial distribution of average annual sediment load for the Dokan watershed sub basins; (b) average slope of sub basins.

Figure 8 .
Figure 8.(a) Spatial distribution of average annual sediment load for the Dokan watershed sub basins; (b) average slope of sub basins.

Table 1 .
Properties of the Dokan and Derbendekhan watersheds.

a.s.l) Min. Elevation (m.a.s.l) Average Slope (%)
Figure 1.Topographic map of the watershed areas of the Dokan and Derbendekhan Dams and their locations in Iraq.

Table 1 .
Properties of the Dokan and Derbendekhan watersheds.

Table 2 .
Percentages of different land use types for two years of the studied period.

Table 2 .
Percentages of different land use types for two years of the studied period.

Table 2 .
Percentages of different land use types for two years of the studied period.

Table 3 .
The range, optimal (fitted) values and sensitivity analysis of considered parameters.
Note: R: Relative, V: Replace.(…) for different soil or plant type.

Table 3 .
The range, optimal (fitted) values and sensitivity analysis of considered parameters.

Table 4 .
Effective parameters arranged from low to high sensitivity based on t-stat and p-value.