The Makhool Dam has brought much attention to the Mosul-Makhool Basin (MMB) in Iraq. Dam construction needs comprehensive studies on the impact of climate change on streamflow and sediment yield in watersheds that supply dams; however, there have been no such investigations conducted on the MMB. An ensemble of three General Circulation Models from Coupled Model Intercomparison Project Phase 6 (CMIP6) with two Shared Socioeconomic Pathways scenarios (SSP2–4.5 and SSP5–8.5) was used to predict the impact of climate change. A distribution mapping downscaling method was utilized to improve the biased climate data. The results indicate that precipitation decreased by 9.5 and 18.7% under SSP2–4.5 and SSP5–8.5 scenarios at the end of the 21st century. The average maximum and minimum temperatures are projected to rise by 2.4 and 1.8 °C in SSP2–4.5 and rise by 3.5 and 2.8 °C in SSP5–8.5 during the study period. The streamflow is expected to decrease by 36.6 and 45.9%, and sediment yield will reduce by 46 and 55% within the same scenarios. Understanding the consequences of climate change helps to face abrupt climate changes to manage dam construction and treat water resources successfully. The results of this research are expected to contribute to improving water management strategies in the study region.

  • The research introduces pre-dam construction investigation.

  • No previous studies explain the impact of climate change on streamflow and sediment yield in the MMB.

  • A CMIP6 project was used to predict the impact of climate change.

  • The study explains the impact of climate change on precipitation and temperature.

  • A guide is provided for decision-makers on managing the water resources.

  • Decision-makers can take steps to mitigate potential challenges that may arise in the future regarding water resources.

Over the past century, the global climate has undergone rapid and notable alteration in the global climate, which is characterized by rising temperatures and an increased frequency of extreme precipitation events. These changes, which are primarily attributed to global warming resulting from excessive carbon emissions, have profoundly impacted hydrological processes (Parmesan et al. 2022). The Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC) revealed that the average temperature of the Earth's surface experienced an increase of approximately 0.85 °C between the years 1880 and 2012 (Stocker et al. 2014). Based on the findings of Reddy et al. (2023) and the IPCC report from 2021 (Allan et al. 2021), the severity and frequency of climate extremes are expected to increase as temperatures continue to rise.
Figure 1

The location of the MMB.

Figure 1

The location of the MMB.

Close modal

To adapt to and mitigate the impacts of climate change, communities require reliable climate data. General Circulation Models (GCMs) serve as a valuable spatial and temporal climate data resource. Essential information about the climate system can be obtained from the GCMs. This information would be difficult to derive from observational data alone (Hamed et al. 2022a; Reddy et al. 2023). Earlier studies utilized GCM datasets to construct the representative concentration pathways (RCPs) scenarios of the CMIP5 project. However, more recent data from Phase Six of the Coupled Model Intercomparison Project (CMIP6) have become available to the global scientific community, significantly advancing the analysis of global climate change. Within CMIP6, the Scenario Model Intercomparison Project (ScenarioMIP) has developed five narratives representing the predicted worldwide social development by the end of the 21st century. These narratives, along with the new RCPs, have produced eight comprehensive scenarios addressing both climate change and society. Collectively, these scenarios are referred to as Shared Socioeconomic Pathways (SSPs) (O'Neill et al. 2017).

The description of SSP1 and SSP5 portrays a future that is characterized by a positive outlook. These scenarios envision the presence of effective institutions, notable acceleration in economic growth, and significant investments in health and education. However, SSP1 presumes a transition toward sustainable energy consumption, whereas SSP5 anticipates a continued reliance on fossil-based energy sources. In contrast, SSP3 and SSP4 depict a less optimistic trajectory, characterized by the limited allocation of resources toward education and health, resulting in heightened increased disparities and population growth. Consequently, these communities are particularly susceptible to the effects of climate change. SSP4 recognizes the presence of severe inter- and intranational inequalities, whereas SSP3 emphasizes regional security as the highest priority. The SSP2 scenario represents the primary pathway, assuming continuing historical trends without significant deviations (O'Neill et al. 2017). The superior performance of CMIP6 compared to CMIP5 has been documented in numerous worldwide studies. Among these, studies are Zamani et al. (2020) in Iran, Kamruzzaman et al. (2021) in Bangladesh, Pour et al. (2022) in Peninsular Malaysia, Hamed et al. (2022a) in Southeast Asia, and Hamed et al. (2022b) in the Middle East and North Africa (MENA). Hamed et al. (2022b) studied the accuracy of the CMIP5 and CMIP6 models in predicting precipitation and temperature patterns in the MENA region. The study results revealed that CMIP6 models demonstrated enhanced simulation capabilities compared to CMIP5 models. Based on the results of Rivera & Arnould (2020), it was observed that most GCMs in the CMIP6 precipitation simulations could effectively replicate climate attributes and long-term precipitation fluctuations. Nevertheless, the utilization of CMIP6 for evaluating streamflow and sediment yield is infrequent, necessitating further investigation to improve the accuracy of simulations in this context.

Successful water resource management, including drought and flood prediction and sustainable agriculture promotion, requires an accurate evaluation of climate change's effects on streamflow (Li & Fang 2021). Numerous studies have focused on investigating the impact of climate change on water resources (Mohseni et al. 2023), revealing that a significant proportion of world's 200 major rivers have experienced notable alterations in streamflow patterns since the 1950s. Rising greenhouse gas concentrations significantly increase the potential risks associated with climate change's consequences on freshwater availability. Many semi-arid and arid regions expect their sustainable surface water and groundwater resources to decrease considerably over the 21st century (Li & Fang 2021). Semi-arid regions exhibit heightened susceptibility to anthropogenic alterations in terrestrial ecosystems and environmental conditions (Stocker et al. 2014). Roughly 15% of the global population resides in these regions, making precipitation an essential resource for satisfying domestic, industrial, and agricultural demands (Schwinning et al. 2004). The influence may be strongly felt in arid and semi-arid regions, especially in western Asia (Salman et al. 2018). The dry climatic system is particularly vulnerable to even minor fluctuations in the current climate due to the fragility of its ecosystems (Salman et al. 2017; Ahmed et al. 2019). The intricate relationship between precipitation and sea surface temperatures in the region gives rise to irregular precipitation patterns and recurring drought conditions (Barlow et al. 2021). Previous research has demonstrated that even minor fluctuations in precipitation levels can lead to a substantial increase in extreme weather events associated with precipitation (Nashwan & Shahid 2020; Salman et al. 2022). Therefore, significant alterations in climatic factors, particularly rainfall, can result in the emergence of climate change-related challenges, such as hydro-meteorological hazards, within the western region of Asia (Salman et al. 2017, 2019; Salehie et al. 2022). Iraq, located in the western region of Asia, is globally recognized as one of the countries with the lowest level of climate resilience. The country has observed a substantial temperature increase, surpassing the global average temperature rise rate over the past 40 years (Salman et al. 2017). Furthermore, it has experienced a significant decrease in rainfall over time by 7.8% (Salman et al. 2022). These changes have substantially impacted various economic sectors and individuals' means of survival.

Hydrological conditions and water availability will change on regional and global scales due to climate change. These changes will influence rainfall, evaporation, water yield, and land surface runoff (Abbas et al. 2022). Droughts and floods threaten various sectors, particularly in developing countries where rain-fed agriculture substantially contributes to the economy (Oguntunde & Abiodun 2013). Multiple studies have provided evidence suggesting that climate change will have profound and significant effects on water resources worldwide, necessitating adaptation strategies (Lehner et al. 2019; Schilling et al. 2020; Li & Fang 2021). Effective adaptation plans require a thorough understanding of global climate change projections and their effects on water resources. Catchment area consequences from climate change are commonly evaluated using rainfall-runoff models and global climate change scenarios. The results highlight the importance of precise climate projections and robust hydrological models capable of withstanding changing climatic conditions (Ruelland et al. 2015). The Soil and Water Assessment Tool (SWAT) is a hydrological, basin-scale tool based on physics. Due to its effectiveness in assessing the impacts of global warming on runoff and sediment yield, it has gained broad acceptance (Zhang et al. 2016).

The future effects of global warming on streamflow and sediment yield have been the focus of much recent research using the SWAT for simulation. The investigation by Azari et al. (2016) focused on the Gorganroud River watershed in northern Iran. The study aimed to determine the impact of global warming on sediment yield and streamflow in this area; the climate change data were derived mostly from the CMIP5 model. De Oliveira et al. (2019) in the Upper Paranabá River Basin, Southeastern Brazil, Tian et al. (2020) in the Tibet Lhasa River Basin, and Sharafati et al. (2020) in the Iran Dehbar River Basin used the CMIP5 model for future climate data under two emission scenarios (RCP4.5 and RCP8.5) to evaluate the effects of global warming on sediment yield and streamflow. Climate data from CMIP5 were also used by Shrestha et al. (2022) in the Songkhram River Basin, Thailand, and dos Santos et al. (2021) in the Caatinga/Atlantic forest ecotone, Brazil to study climate change effects on streamflow and sediment yield.

Numerous studies indicate that variations in precipitation and temperature may have a considerable influence on watershed processes, while the magnitude and nature of this impact can vary widely from region to region. Given the anticipated global negative effects of climate change across all regions, it is crucial to quantify the hydrological impacts to enhance our understanding and predictive capabilities concerning sediment yield and streamflow processes. This is particularly important as part of a pre-dam construction study. Therefore, the objective of this study is to assess the potential consequences of climate change on streamflow and sediment yield in the Tigris River Basin, with specific emphasis on the Mosul-Makhool Basin (MMB). To achieve this, we will utilize the SWAT model and CMIP6 data from three GCMs, analyzing the time period from 2021 to 2100. The analysis will be conducted under two emissions scenarios, namely SSP2–4.5 and SSP5–8.5, and is divided into four time intervals: 2021–2040, 2041–2060, 2061–2080, and 2081–2100. The use of 20-year intervals serves as a practical and scientifically justifiable approach, aligning with the requirements of dam construction planning, risk management, and adaptation strategies. This approach allows us to capture the evolving nature of climate change and its impacts on streamflow and sediment yield, which are critical considerations for the proposed dam project. The findings from this study will significantly contribute to decision-making processes, particularly in relation to the proposed Makhool Dam. Moreover, the results will promote sustainable practices in water resource management in arid and semi-arid regions facing the challenges of climate change.

Study area

The MMB is a sub-basin of the Tigris River Basin, which is located in the northeastern region of Iraq. It borders Iran and Turkey to the east and the north, respectively (Figure 1). The basin spans 4.51 × 104 km2 at the location of the Makhool Dam at 140 meters above mean sea level (m.a.s.l.). The basin has various climatic conditions, from low-lying parts at 140 m.a.s.l. to the high-altitude Zagros Mountains at 4,168 m.a.s.l. (Osman et al. 2019). The climatic conditions in the northern and northeastern regions of Iraq are generally typical of a Mediterranean climate. The annual precipitation exhibits significant variability, with the mountainous areas in the northeast experiencing precipitation levels exceeding 1,500 mm, while the location of the Makhool Dam receives less than 350 mm. The predominant precipitation in the river basin is observed during the spring and winter seasons, mostly as rain and snowmelt in Turkey and Iran's mountainous headwater regions. Winter temperatures in the basin can drop below freezing (0 °C) in January, while temperatures in February rise above 18 °C. Conversely, the summer season consistently witnesses extremely high temperatures, frequently surpassing 45 °C (Al-Hasani 2019). The average annual streamflow in the Baji station is 779 m3/s, varying from 8,947 m3/s in spring to 409 m3/s in autumn. The geological composition of the basin's mountainous regions is primarily distinguished by the prevalence of carbonate rocks and clastic rocks. The western part of the watershed mainly consists of silt, clay, rare sand, and pebbles (Sissakian et al. 2018). The MMB is currently facing significant challenges related to soil erosion, sediment yields, floods, and debris flow. These problems stem from rainfall-induced erosion in sensitive soils, aggravated by extensive changes in land use, specifically the change of range lands and forests into arid areas. These factors have caused increased runoff, leading to significant soil erosion and sediment accumulation within the watershed.

Data collection

Hydro-meteorological data

The daily precipitation data used in the study were sourced from the Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS) dataset, which has a spatial resolution of 0.05° × 0.05°. CHIRPS is a comprehensive precipitation dataset that offers high-resolution information from 1981 to the present. This dataset represents a gridded time series of rainfall information derived from a combination of satellite data and data from in situ rain gauge stations. By incorporating both sources of information, this approach offers a comprehensive and precise depiction of precipitation patterns (Katsanos et al. 2016). The CHIRPS product and the related data are freely available on the official website: http://chg.geog.ucsb.edu/data/chirps/. Numerous studies have underscored the superior performance of CHIRPS compared to other datasets in simulating streamflow within the SWAT framework (Duan et al. 2019; Dhanesh et al. 2020; Mengistu et al. 2022).

The daily maximum and minimum temperature (Tmax and Tmin) data for this study were obtained from the reanalysis products provided by the National Aeronautics and Space Administration Prediction of Worldwide Energy Resource (NASA POWER). These datasets offer extensive information on daily, monthly, and annual Tmin and Tmax at a spatial resolution of 0.5° × 0.5° (Natumanya et al. 2022). The NASA POWER datasets can be accessed through the data access viewer at the following URL: https://power.larc.nasa.gov/data-access-viewer/. Given the limited number of ground observation stations and the associated data quality concerns, 45 meteorological stations were used in the MMB, which were provided by the CHIRPS and NASA POWER.

The outflow data from the Mosul Dam and the Dokan Dam, covering the period from 1994 to 2020, were sourced from the National Center for Water Resources, an entity operating under the Iraqi Ministry of Water Resources. The streamflow observations used in this study were of high quality, and there were no missing data points in the dataset. An interpolation method was employed if minor gaps occurred to ensure data continuity and reliability.

Topographic, soil, and land use/land cover data

This research utilized 30 m spatial resolution data of the Digital Elevation Model (DEM). The DEM was created by using Shuttle Radar Topography Mission (SRTM) data and obtained from the official United States Geological Survey (USGS) website at https://earthexplorer.usgs.gov/.

The study utilized the 1995 global soil map from the Food and Agriculture Organization of the United Nations for soil data (http://www.fao.org). The land use/land cover (LULC) data were sourced from the Environmental Systems Research Institute (ESRI), with a spatial resolution of 10 m. The official data source provided by the ESRI can be accessed at https://livingatlas.arcgis.com/landcover/. Throughout the model simulation process, the LULC remained constant, as the primary objective was to evaluate climate change's effect on streamflow and sediment yield.

SWAT model

The United States Department of Agriculture (USDA) developed the SWAT hydrological model. Arnold et al. (1998) described it as a process-based and semi-distributed model that operates continuously over time to simulate various hydrological processes within watersheds. It is valuable for simulating and modeling the quantity and quality of groundwater and surface water and their responses to multiple land management techniques (Li & Fang 2021). The effective operation of the SWAT requires hydrological data, daily precipitation and temperature, DEM, soil information, and LULC data. Figure 2 shows a schematic representation of the SWAT model. The datasets partition a basin into sub-basins and hydrologic response units (HRUs) using predetermined thresholds associated with soil type, LULC, and slope. This partitioning procedure permits a more accurate representation of the basin's characteristics and improves the accuracy of hydrological modeling by operating at a higher spatial resolution (Arnold et al. 2012). This research subdivided the MMB into 113 sub-basins and 1,526 HRUs.
Figure 2

Schematic representation of the SWAT model.

Figure 2

Schematic representation of the SWAT model.

Close modal
Within the SWAT, the hydrological processes within each HRU are governed by the water balance equation. This equation considers various daily factors such as runoff, precipitation, return flow, evapotranspiration, and percolation. To estimate surface runoff in the SWAT, two different methods are employed: the Soil Conservation Service Curve Number (SCS-CN) method and the Green and Ampt infiltration method. The SCS-CN method estimates the runoff quantity by considering the variability in soil types and LULC conditions. In contrast, the Green and Ampt equation estimates infiltration, assuming that water is abundant at the surface during the entire simulation period. The SWAT can simulate soil erosion and sediment yield from in-stream channel processes and hillslopes. The Modified Universal Soil Loss Equation estimates soil degradation caused by runoff and rainfall (Neitsch et al. 2011). The Muskingum method or the variable storage coefficient method is used for river flow routing (Arnold et al. 1998). Figure 3 provides the methodology used in this investigation.
Figure 3

A schematic diagram representation of the methodology used in this study.

Figure 3

A schematic diagram representation of the methodology used in this study.

Close modal

Calibration and verification of models

Sensitivity analysis, calibration, and validation of models are essential to accurately predicting the potential impact of climate change on sediment yield and runoff. By reviewing the model's performance statistics within acceptable ranges, these procedures establish its acceptability.

This study used the SWAT Calibration and Uncertainty Program (SWAT-CUP) to calibrate and validate the hydrological model. The effectiveness of the model is assessed using various performance metrics. Nash–Sutcliffe efficiency (NSE), percent bias (PBIAS), coefficient of determination (R2), and root mean squared error-observation deviation ratio (RSR) are all utilized in this investigation. Detailed information regarding these performance indicators can be found in Moriasi et al. (2007). These metrics comprehensively assess the model's performance and capacity to accurately simulate streamflow and sediment yield.

Among the options available in the SWAT-CUP, the Sequential Uncertainty Fitting (SUFI-2) algorithm was selected for this investigation. Previous studies have demonstrated that the SUFI-2 algorithm exhibits fewer iterations when optimizing parameters than alternative methods like particle swarm optimization and Generalized Likelihood Uncertainty Estimation. The SUFI-2 algorithm has gained recognition for its notable efficacy in optimizing parameters and analyzing uncertainty within the SWAT model (Abbaspour et al. 2017). These findings indicate that the SUFI-2 algorithm offers a more efficient approach for parameter optimization within the context of SWAT-CUP.

Twenty parameters were selected to calibrate the model of the MMB (Table 1). The sensitivity parameters were evaluated utilizing p-values and t-stats. Parameters with lower p-values and higher t-stats indicate a heightened sensitivity to the optimization function, suggesting their substantial influence on the model's performance. This sensitivity analysis helps identify the crucial parameters that play a key role in improving the accuracy and reliability of the model through the calibration process (Bhatta et al. 2019). After completing sensitivity analysis, calibration and validation were carried out. The calibration phase utilized data spanning from 1994 to 2011, while the subsequent validation period extended from 2012 to 2020.

Table 1

The sensitive parameters for streamflow of the MMB with their maximum, minimum, and fitted values

No.ParameterInitial values
Final values
Min_valueMax_valuet-statp-valueMin_valueMax_valueFitted value
ALPHA_BNK.rte −0.5 0.5 8.84 0.00 0.00 0.400 0.268 
CH_K2.rte −5 150 −5.41 0.00 1.00 200.00 69.257 
CH_N2.rte 0.2 −4.98 0.00 0.058 0.167 0.088 
CN2.mgt −0.5 0.5 −3.80 0.00 −0.280 0.359 −0.097 
ESCO.hru −0.5 0.6 1.85 0.06 0.00 0.798 0.357 
SOL_K(..).sol −60 500 −1.62 0.11 −324.763 836.326 −258.581 
SOL_AWC(..).sol 0.5 0.8 1.23 0.22 0.653 0.987 0.717 
SLSUBBSN.hru 15 1.22 0.22 10.00 20.00 16.290 
SMTMP.bsn −5 15 1.00 0.32 −2.474 5.312 3.887 
10 SFTMP.bsn −10 10 0.95 0.34 −0.203 5.693 1.690 
11 REVAPMN.gw −50 300 0.90 0.37 −274.970 316.050 −111.257 
12 GWQMN.gw 10 0.87 0.38 3.726 15.799 12.745 
13 SOL_BD(..).sol 2.5 −0.70 0.48 0.143 2.582 0.950 
14 RCHRG_DP.gw 0.5 0.49 0.62 0.00 0.500 0.119 
15 GW_REVAP.gw 0.25 0.65 −0.45 0.66 0.249 0.435 0.364 
16 HRU_SLP.hru 0.5 0.38 0.71 0.800 1.00 0.813 
17 GW_DELAY.gw −50 375 0.32 0.75 −296.962 451.885 418.187 
18 ALPHA_BF.gw 0.5 −0.21 0.83 0.00 0.511 0.252 
19 OV_N.hru −0.1 0.1 −0.20 0.84 0.010 0.100 0.015 
20 TIMP.bsn 0.1 0.5 −0.20 0.85 −0.031 0.522 0.117 
No.ParameterInitial values
Final values
Min_valueMax_valuet-statp-valueMin_valueMax_valueFitted value
ALPHA_BNK.rte −0.5 0.5 8.84 0.00 0.00 0.400 0.268 
CH_K2.rte −5 150 −5.41 0.00 1.00 200.00 69.257 
CH_N2.rte 0.2 −4.98 0.00 0.058 0.167 0.088 
CN2.mgt −0.5 0.5 −3.80 0.00 −0.280 0.359 −0.097 
ESCO.hru −0.5 0.6 1.85 0.06 0.00 0.798 0.357 
SOL_K(..).sol −60 500 −1.62 0.11 −324.763 836.326 −258.581 
SOL_AWC(..).sol 0.5 0.8 1.23 0.22 0.653 0.987 0.717 
SLSUBBSN.hru 15 1.22 0.22 10.00 20.00 16.290 
SMTMP.bsn −5 15 1.00 0.32 −2.474 5.312 3.887 
10 SFTMP.bsn −10 10 0.95 0.34 −0.203 5.693 1.690 
11 REVAPMN.gw −50 300 0.90 0.37 −274.970 316.050 −111.257 
12 GWQMN.gw 10 0.87 0.38 3.726 15.799 12.745 
13 SOL_BD(..).sol 2.5 −0.70 0.48 0.143 2.582 0.950 
14 RCHRG_DP.gw 0.5 0.49 0.62 0.00 0.500 0.119 
15 GW_REVAP.gw 0.25 0.65 −0.45 0.66 0.249 0.435 0.364 
16 HRU_SLP.hru 0.5 0.38 0.71 0.800 1.00 0.813 
17 GW_DELAY.gw −50 375 0.32 0.75 −296.962 451.885 418.187 
18 ALPHA_BF.gw 0.5 −0.21 0.83 0.00 0.511 0.252 
19 OV_N.hru −0.1 0.1 −0.20 0.84 0.010 0.100 0.015 
20 TIMP.bsn 0.1 0.5 −0.20 0.85 −0.031 0.522 0.117 

Projection of future climate change

GCMs and selection of specific CMIP models

General circulation models (GCMs) are widely accepted as effective and reliable for assessing global climate change. The data from the CMIP6 were utilized to evaluate the potential future changes in temperature and precipitation. The CMIP6 was preferred due to its improved precision and resolution compared to the previous iteration (CMIP5) (Nourani et al. 2022) and the incorporation of new scenarios that integrate socioeconomic aspects of climate change. This research employed ensemble data from three global climate models to predict the future climate and spatiotemporal variations of hydromorphological changes in the study region (Table 2). The decision of using three specific GCMs in our research was informed by recent studies in the field of climate modeling, specifically Salman et al. (2022) and Hamed et al. (2023). These studies have identified these three GCMs as the most suitable choices for simulating climate conditions in the Middle East and Iraq. Their findings align with our research objectives and the geographic focus of our study. The historical period (1994–2020) and the future (2021–2100) period were used for the analysis.

Table 2

Global climate models (GCMs) considered in this study

ModelInstitutionResolutionCountry
ACCESS-CM2 Australian Research Council Centre of Excellence for Climate System Science 1.25° × 1.875° Australia 
BCC-CSM2-MR Beijing Climate Center 1.125° × 1.125° China 
MRI-ESM2 − 0 Meteorological Research Institute 1.125° × 1.125° Japan 
ModelInstitutionResolutionCountry
ACCESS-CM2 Australian Research Council Centre of Excellence for Climate System Science 1.25° × 1.875° Australia 
BCC-CSM2-MR Beijing Climate Center 1.125° × 1.125° China 
MRI-ESM2 − 0 Meteorological Research Institute 1.125° × 1.125° Japan 

The CMIP6 has replaced the RCPs with the SSPs for exploring potential interactions between socioeconomic development pathways and resulting greenhouse gas emissions and climate change. These pathways incorporate factors, including population growth, economic development, energy use, and land use changes, to generate integrated assessments of future climate change by linking with global climate models. Potential consequences of global warming on sediment yield and streamflow in the studied basin were assessed utilizing two emission scenarios: SSP2–4.5 (medium emission) and SSP5–8.5 (high emission). The scenario known as SSP2–4.5 indicates a moderate degree of social vulnerability and radiative forcing, whereas SSP5–8.5 is distinguished by the highest level of anthropogenic radiative forcing, reaching 8.5 W/m2 by the year 2100 (Zhang et al. 2019).

Bias correction of climate data

This study employed the Climate Model data for hydrologic modeling (CMhyd) tool for statistical bias correction. CMhyd is a Python-based software tool designed to downscale and bias-correct climate model data obtained from general circulation models (GCMs) for hydrologic modeling purposes. GCMs typically have coarse spatial resolution and may produce biased climate projections unsuitable for SWAT use (Kumar et al. 2022). CMhyd utilizes various statistical downscaling techniques to generate high-resolution climate projections. This paper specifically employed the distribution mapping method for bias correction. This method was chosen based on studies highlighting its effectiveness in correcting biases in climate models (Smitha et al. 2018; Enayati et al. 2021). The adjustments to precipitation and temperature were carried out using a transfer function derived from the monthly mean values of the projected GCM data. This approach aimed to enhance the accuracy of simulations by fine-tuning the GCMs' outputs and reducing the discrepancies between observed and simulated climate variables (Siabi et al. 2023). By implementing these bias correction techniques via CMhyd, the researchers intended to improve the accuracy and reliability of the hydrological model's simulations, thereby ensuring more realistic and significant findings for their study (Vaittinada Ayar et al. 2021).

Sensitivity analysis, calibration, and validation of parameters

Twenty parameters were selected to calibrate the SWAT-based model of the MMB (SWAT-MMB). These selected parameters are illustrated in Table 1. In SWAT-CUP, the calibration process involves conducting 500 simulations. The effectiveness of optimization trials depends on the results of sensitivity analyses performed on various model parameters. The sensitivity analysis of these parameters involves assessing their significance by considering p-values and t-stats. A large absolute t-stats indicates that the parameter significantly affects the model output, while a small p-value (typically less than 0.05) indicates that the parameter has a significant effect on the model output. In general, parameters with small p-values and large t-stats are considered significant contributors to the model (Abbaspour et al. 2017). Among the parameters considered, it was observed that the Alpha factor of the base flow in bank storage (ALPHA_BNK) exhibited the greatest sensitivity in the MMB. Following ALPHA_BNK, the parameters that displayed significant sensitivity were the main channel's effective hydraulic conductivity (CH_K2), Manning's value (n) for the main channel (CH_N2), SCS runoff curve number (CN2), and the compensating factor for soil evaporation (ESCO). It is worth mentioning that during the model calibration process, the parameters related to routing water flow (ALPHA_BNK, CH_K2, and CH_N2) and parameters associated with surface runoff (CN2) were found to be the most sensitive. The results agree with the studies of Abbas et al. (2016) and Tian et al. (2020). The SWAT-MMB model was found to be reasonably influenced by snow-related factors like SFTMP and SMTMP. This influence was especially pronounced in high-altitude areas where winter temperatures consistently drop below freezing, leading to an accumulation of snow cover that melts during spring months.

Through this study, the SWAT-MMB model was calibrated and validated using monthly streamflow. The concurrence between observed and simulated flow was assessed using various evaluation metrics, including NSE, R2, PBIAS, and RSR. The calibration process used a timeframe of 1994–2011, while the validation process used a timeframe of 2012–2020. As shown in Table 3, the model demonstrated significant consistency between the simulated and observed data. This suggests that the model performed ‘very well’ in replicating the monthly streamflow patterns during both the calibration and validation periods. The simulation results of the model are considered reliable when certain performance metrics meet specific criteria. According to previous studies (Moriasi et al. 2007; Abou Rafee et al. 2019), a model is deemed satisfactory if it achieves NSE and R2 values above 0.5 and 0.6, respectively, while maintaining RSR below 0.7 and PBIAS below 25%. However, obtaining a robust fit was difficult due to the complexity of hydrological processes and the substantial monthly variations in runoff. Overall, the model predicts streamflow in the MMB with sufficient accuracy, and the future runoff and sediment yield as a consequence of climate change could be simulated with the optimum model parameters (Figure 4).
Table 3

The performance of the SWAT model in simulating streamflow was evaluated during both the calibration and validation processes

Parameter usedMMB
CalibrationValidation
NSE 0.93 0.82 
R2 0.94 0.84 
RSR 0.27 0.42 
PBIAS −4.3 2.8 
Parameter usedMMB
CalibrationValidation
NSE 0.93 0.82 
R2 0.94 0.84 
RSR 0.27 0.42 
PBIAS −4.3 2.8 
Figure 4

Comparison between observed and best simulated monthly streamflow data at the Baji station during calibration and validation periods.

Figure 4

Comparison between observed and best simulated monthly streamflow data at the Baji station during calibration and validation periods.

Close modal

In this study, the validity of precipitation and temperature data obtained from CHIRPS and NASA-POWER datasets was verified by comparing them with measurements from two weather stations in the study region, such as Tikrit Station and Mosul Station. This comparison involved calculating NSE, R2, PBIAS, and RSR metrics (Table 4). The results indicated significant consistency between observed data from the weather stations and corresponding data from CHIRPS and NASA-POWER datasets, confirming the reliability and suitability of these datasets for precipitation and temperature analysis in the study area.

Table 4

The evaluation of the performance of CHIRPS and NASA-POWER datasets to obtain the precipitation and temperature data, respectively, for both the calibration and validation processes

Parameter usedCHIRPS
NASA-POWER
Tikrit Station
Mosul Station
Tikrit Station
Mosul Station
CalibrationValidationCalibrationValidationCalibrationValidationCalibrationValidation
NSE 0.85 0.80 0.78 0.74 0.71 0.67 0.61 0.64 
R2 0.90 0.85 0.92 0.95 0.74 0.73 0.70 0.65 
RSR 0.34 0.39 0.27 0.31 0.50 0.55 0.49 0.51 
PBIAS −5.2 3.4 −4.7 4.3 11.5 13.2 −15.6 16.4 
Parameter usedCHIRPS
NASA-POWER
Tikrit Station
Mosul Station
Tikrit Station
Mosul Station
CalibrationValidationCalibrationValidationCalibrationValidationCalibrationValidation
NSE 0.85 0.80 0.78 0.74 0.71 0.67 0.61 0.64 
R2 0.90 0.85 0.92 0.95 0.74 0.73 0.70 0.65 
RSR 0.34 0.39 0.27 0.31 0.50 0.55 0.49 0.51 
PBIAS −5.2 3.4 −4.7 4.3 11.5 13.2 −15.6 16.4 

Projected future temperature and precipitation

The present study employed an ensemble of three global climate models (ACCESS-CM2, BCC-CSM2-MR, and MRI-ESM2–0) to represent future climate conditions regarding precipitation and temperature (Tmax, Tmin). Two emission scenarios were considered: SSP2–4.5 (medium) and SSP5–8.5 (high). The bias-corrected data were compared with the baseline climatic data from 1994 to 2020 to identify the precise changes in precipitation and temperature. Before developing climate scenarios, an evaluation was conducted to assess the effectiveness of the distribution mapping technique for bias correction. For evaluating the global climate models, the raw data were compared with the bias-corrected data using R2, NSE, RSR, and PBIAS. Table 5 shows a significant consistency between the bias-corrected and observed data assessed by R2, NSE, RSR, and PBIAS. Figure 5 shows the comparison between observed, uncorrected, and bias-corrected data for precipitation, Tmax, and Tmin for the baseline climatic period.
Table 5

Performance of the distribution mapping bias correction method

Parameter usedPrecipitationTmaxTmin
NSE 0.75 0.96 0.97 
R2 0.97 0.997 0.99 
RSR 0.00 0.04 0.03 
PBIAS 8.00 9.04 12.50 
Parameter usedPrecipitationTmaxTmin
NSE 0.75 0.96 0.97 
R2 0.97 0.997 0.99 
RSR 0.00 0.04 0.03 
PBIAS 8.00 9.04 12.50 
Figure 5

The comparison between observed, uncorrected, and bias-corrected data for (a) precipitation, (b) Tmax, and (c) Tmin.

Figure 5

The comparison between observed, uncorrected, and bias-corrected data for (a) precipitation, (b) Tmax, and (c) Tmin.

Close modal
In the MMB region, the average annual precipitation currently stands at 508 mm, with maximum and minimum temperatures of 27.3 and 13.8 °C, respectively. Future projections indicate a rise in both average minimum and maximum temperatures. Under the SSP2–4.5 scenario, the average annual minimum temperature is expected to increase by 1.5, 2.0, 2.8, and 3.5 °C for the periods 2021–2040s, 2041–2060s, 2061–2080s, and 2081–2100s. In contrast, under the SSP5–8.5 scenario, the increases are projected to be 1.4, 2.7, 4.3, and 5.8 °C for the same time intervals. Average annual maximum temperatures are projected to rise, with increases of 1.1, 1.5, 2.1, and 2.5 °C under SSP2–4.5 and 1.0, 2.2, 3.5, and 4.7 °C under SSP5–8.5 over the same time periods. Notably, the study found that the increase in minimum temperatures exceeded that of maximum temperatures, potentially leading to consequences like reduced snow accumulation, shorter winters, earlier snowmelt in mountainous areas, and changes in ecosystems and agricultural patterns. These disparities between minimum and maximum temperatures were most pronounced during the summer months, particularly in July and August. The heightened impact during the summer months accentuates the need for adaptation and mitigation strategies to address the challenges posed by these changing climate conditions. Both emission scenarios exhibit similar temperature change patterns, with overall temperature increases expected across all seasons. However, the SSP5–8.5 scenario, which is characterized by higher radiative forcing and greenhouse gas emissions, projects more substantial temperature increases compared to SSP2–4.5. Figure 6 depicts the absolute changes in the MMB temperature relative to the baseline period. The findings of this study underscore the urgency of addressing climate change in the MMB region. With the prospect of significant temperature increases in both emission scenarios, proactive measures to adapt to changing conditions and reduce emissions are crucial to safeguarding the region's ecosystems, agriculture, and overall well-being in the face of a warming climate.
Figure 6

Absolute change concerning the baseline period (1994–2020) in the relative change in precipitation (top), maximum temperature (middle), and minimum temperature (bottom) in the MMB.

Figure 6

Absolute change concerning the baseline period (1994–2020) in the relative change in precipitation (top), maximum temperature (middle), and minimum temperature (bottom) in the MMB.

Close modal
In Iraq, the distribution of rainfall exhibits a distinct seasonality pattern. The rainy season typically unfolds during the months of January, February, March, and April, collectively contributing to over 65% of the annual precipitation in the region. During this period, the country experiences its highest rainfall rates, supporting critical agricultural and ecological processes. In contrast, the months of June, July, August, and September mark the dry season, with minimal precipitation. These months receive a mere 0.05% of the yearly rainfall, making their contribution to the overall precipitation in Iraq negligible. As a result of global warming, precipitation is expected to reduce in the future under both SSP scenarios. The average annual rainfall is projected to decrease by 4.8, 9.2, 15.4, and 4.9% during the periods 2021–2040s, 2041–2060s, 2061–2080s, and 2081–2100s, respectively, under the SSP2–4.5 scenario. Meanwhile, under the SSP5–8.5 scenario, a substantial decrease in the average rainfall is anticipated. Specifically, during the same periods, we expect a decrease of 11.9, 16.6, 17.5, and 25.5%. The anticipated changes in precipitation patterns in the study area are illustrated in Figure 6, showing the predicted changes in rainfall compared to the baseline period. According to the projections in Figure 6, a slight increase in the average precipitation is projected for March and April, particularly during the 2021–2040s and 2041–2060s under the SSP2–4.5 scenario. An increase in the average rainfall is anticipated during November and December in the years 2081–2100, as shown in Figure 6. Throughout all periods and under both scenarios, precipitation levels in June, July, August, and September are projected to approach zero. It is worth noting that, across all seasons during future periods, the rate of decrease is anticipated to be higher under the SSP5–8.5 scenario compared to the SSP2–4.5 scenario. These findings confirm those of Osman et al. (2017) and Salman et al. (2022) in the north of Iraq. Figure 7 shows the changes in the frequency of extreme precipitation events over different periods compared to the baseline period. Extreme precipitation events are defined here as those that exceed the 99th percentile of historical precipitation amounts, which was 20 mm/day. The frequency of extreme precipitation events shows a decreasing trend toward the end of the century for all periods and under both emission scenarios.
Figure 7

The changes in the precipitation frequency of extreme events over different periods compared with the historical period.

Figure 7

The changes in the precipitation frequency of extreme events over different periods compared with the historical period.

Close modal
Figures 810 demonstrate the spatial distributions of annual precipitation, mean maximum and minimum temperatures for the baseline period (1994–2020), and future periods. The inverse distance weighted interpolation technique and geographical information systems were utilized to map the spatial distribution of temperature and precipitation, as well as to analyze their trends and variability. The lowest values of Tmax and Tmin occurred in the northern region of the watershed, while the highest values occurred in the southern region. The highest annual precipitation was recorded in the middle of the region. Precipitation trends from 2020 until 2100 indicate a decrease in the northern region of 6.1 and 27.1% under both the SSP2–4.5 and SSP5–8.5 scenarios. In the southern region, precipitation is projected to decrease by 7% under SSP2–4.5 and 25% under SSP5–8.5 scenarios. Similarly, in the middle regions, there is an expected decline of 4% under the SSP2–4.5 scenario and a more substantial decrease of 25% under SSP5–8.5.
Figure 8

Spatial distribution of annual precipitation in the MMB.

Figure 8

Spatial distribution of annual precipitation in the MMB.

Close modal
Figure 9

Spatial distribution of mean Tmax in the MMB.

Figure 9

Spatial distribution of mean Tmax in the MMB.

Close modal
Figure 10

Spatial distribution of mean Tmin in the MMB.

Figure 10

Spatial distribution of mean Tmin in the MMB.

Close modal

Regarding temperature changes, Tmax is anticipated to increase by 2.5 °C under the SSP2–4.5 scenario and 4.8 °C under the SSP5–8.5 scenario in the northern region, while, in the southern region, it is projected to increase by 2.6 °C under SSP2–4.5 and 4.8 °C under SSP5–8.5. Similarly, Tmin is expected to increase by 3.4 and 6.1 °C under both the SSP2–4.5 and SSP5–8.5 scenarios in the northern region, and in the southern region, Tmin is projected to increase by 3.4 °C under SSP2–4.5 and 5.8 °C under SSP5–8.5.

Impact of climate change on the streamflow

Changes in climatic factors can have an impact on various components of the water balance within a catchment area, leading to subsequent modifications in river flow patterns. This study presents streamflow data for the MMB under two emission scenarios: SSP2–4.5 and SSP5–8.5. The projected annual average streamflow for the period 2021–2100 indicates a reduction in streamflow under both the SSP2–4.5 and SSP5–8.5 scenarios compared to the baseline period of 1994–2020 (see Figure 11(a)). Specifically, the average annual streamflow demonstrates a gradual decrease of 31.4, 37.8, 40.5, and 36.7% during the time windows of 2021–2040s, 2041–2060s, 2061–2080s, and 2081–2100s, respectively, under the SSP2–4.5 scenarios. According to the SSP5–8.5 scenario, reductions of 32.0, 47.4, 46.0, and 58.0% are evident during the same time periods. Compared to the baseline period, climate change under the SSP2–4.5 scenario observed a decrease in streamflow by 36.6% from 2021 to 2100. On the other hand, the SSP5–8.5 scenario exhibited a larger reduction, with streamflow decreasing by 45.9% over the same periods.
Figure 11

Changes in average annual (a), streamflow (b), wet season, and (c) dry season for the periods 2021–2040s, 2041–2060s, 2061–2080s, and 2081–2100s.

Figure 11

Changes in average annual (a), streamflow (b), wet season, and (c) dry season for the periods 2021–2040s, 2041–2060s, 2061–2080s, and 2081–2100s.

Close modal

Figure 11(b) and 11(c) illustrates the changes in streamflow during the wet and dry seasons across the four future periods under the SSP2–4.5 and SSP5–8.5 scenarios. The observed effects of climate change on streamflow reveal a consistent pattern for both wet and dry seasons. Specifically, in the wet season, a decrease in streamflow is projected across all future periods. Under the SSP2–4.5 scenario, percentage changes in streamflow are 14.1, 20.2, 14.1, and 6.7% during the periods of 2021–2040s, 2041–2060s, 2061–2080s, and 2081–2100s, respectively. Similarly, under the SSP5–8.5 scenario, streamflow exhibits a reduction of 17.0, 30.0, 19.9, and 32.9% during the same time periods (Figure 11(b)). In the dry season, streamflow consistently decreases across all future decades under both SSP scenarios (Figure 11(c)). This reduction in streamflow during the dry season is particularly concerning as it exacerbates challenges related to water scarcity and availability during a period when demand is often at its peak. A stepwise multiple linear regression analysis was conducted to examine the effects of temperature and precipitation on streamflow. The results revealed that during the wet season, the influence of precipitation on streamflow had greater significance than that of temperature. This finding clarifies that the decline in streamflow during the wet season is anticipated to be less severe than the reductions observed during the dry season. Temperature is projected to have a more significant impact than precipitation during the dry season in future climate change scenarios (Table 6). The insight that precipitation remains a more influential factor than temperature in driving streamflow during the wet season offers a ray of hope for water resource management strategies. While climate change is projected to bring about reductions in streamflow during this critical period, there is an opportunity for mitigation through the enhancement of precipitation harvesting and storage techniques. During the dry season, both maximum and minimum temperatures are negatively correlated with streamflow, adversely impacting the water flow in the streams. This observation suggests that evaporation rates increase as temperatures rise, surpassing the effect of precipitation and reducing streamflow (Bhatta et al. 2019). Consistent with this, Azari et al. (2016) and Li & Fang (2021) reported a reduction in streamflow during the dry season and an increase during the wet season in various regions worldwide.

Table 6

The standardized coefficients for the relationship between climate variables and streamflow according to the stepwise regression analysis

VariableStreamflow
JanFebMarAprMayJunJulAugSepOctNovDec
Baseline Prcp 0.453**         −0.490* 0.496** 0.561** 
Tmin 0.522** 0.870** 0.981*  −0.403*  −0.380* −0.512* −0.412*    
Tmax  −0.760** −1.28** −0.607**  −0.540*       
SSP2–4.5 Prcp 0.585** 0.419** 0.431** 0.335**       0.741** 0.546** 
Tmin     −0.302* −0.214* −0.274*  −0.331** −0.376**   
Tmax        −0.540**     
SSP5–8.5 Prcp 0.620** 0.644** 0.518** 0.328**       0.527** 0.465* 
Tmin     −0.448** −0.398**    −0.655** −0.337**  
Tmax       −0.702** −0.657** −0.485**    
VariableStreamflow
JanFebMarAprMayJunJulAugSepOctNovDec
Baseline Prcp 0.453**         −0.490* 0.496** 0.561** 
Tmin 0.522** 0.870** 0.981*  −0.403*  −0.380* −0.512* −0.412*    
Tmax  −0.760** −1.28** −0.607**  −0.540*       
SSP2–4.5 Prcp 0.585** 0.419** 0.431** 0.335**       0.741** 0.546** 
Tmin     −0.302* −0.214* −0.274*  −0.331** −0.376**   
Tmax        −0.540**     
SSP5–8.5 Prcp 0.620** 0.644** 0.518** 0.328**       0.527** 0.465* 
Tmin     −0.448** −0.398**    −0.655** −0.337**  
Tmax       −0.702** −0.657** −0.485**    

*Significant at the 0.05 level.

**Highly significant at the 0.01 level.

Figure 12 illustrates the average monthly streamflow changes from 2021 to 2100 under the SSP2–4.5 and SSP5–8.5 emission scenarios. The findings indicate a consistent decline in streamflow across all months under both scenarios. Moreover, the decline in streamflow is expected to be more substantial in the SSP5–8.5 scenario compared to the SSP2–4.5 scenario. The difference between the two scenarios becomes notably more distinct in April, while it remains relatively stable in the remaining months. Notably, the peak flow is projected to shift from May to April under both emission scenarios. The primary factor driving this shift in peak streamflow is the accelerated thawing of snowpacks due to global warming (Gould et al. 2016; Romshoo & Marazi 2022). Climate change poses a significant threat to water security. This is evident in the projected changes in peak streamflow magnitude and frequency under both the SSP2–4.5 and SSP5–8.5 emission scenarios. Under SSP2–4.5, the magnitude of peak streamflow is expected to increase notably during the 2021–2040, 2041–2060, and 2081–2100 periods, with the exception of the 2061–2080 period. However, the frequency of peak streamflow is expected to decrease in future periods compared to the historical period (Figure 13). Under SSP5–8.5, the magnitude and frequency of peak streamflow are both expected to decrease. Specifically, the periods of 2061–2080 and 2081–2100 show no flow exceeding the threshold of the 95th percentile of historical monthly peak streamflow (2001–2020). These findings suggest that climate change will lead to more extreme precipitation events, such as floods, as well as more severe droughts.
Figure 12

Monthly mean streamflow comparisons between the baseline period and emission scenarios for the duration of the projection period.

Figure 12

Monthly mean streamflow comparisons between the baseline period and emission scenarios for the duration of the projection period.

Close modal
Figure 13

Comparing changes in peak streamflow frequency and magnitude with historical data.

Figure 13

Comparing changes in peak streamflow frequency and magnitude with historical data.

Close modal

Impact of climate change on sediment yield

The projected annual decrease in sediment yield during the periods 2021–2040s, 2041–2060s, 2061–2080s, and 2081–2100s is estimated to be 38.9, 45.4, 52.2, and 47.6%, respectively, under the SSP2–4.5 scenario (Figure 14(a)). Similarly, the SSP5–8.5 scenario suggests a decrease in sediment yield throughout all periods. These expectations are based on extracting sediment yield values from a model calibrated using streamflow data, as sediment gauge stations are scarce. The behavior of sediments during climate changes is similar to that of streamflow, with only slight variations in proportions.
Figure 14

Changes in mean annual (a), sediment yield (b), wet season, and (c) dry season for the periods 2021–2040s, 2041–2060s, 2061–2080s, and 2081–2100.

Figure 14

Changes in mean annual (a), sediment yield (b), wet season, and (c) dry season for the periods 2021–2040s, 2041–2060s, 2061–2080s, and 2081–2100.

Close modal

Figure 14(b) and 14(c) illustrates the variations in sediment yield during wet and dry seasons over specific periods. The future projections indicate a decline in sediment during the wet season for all four future periods under the SSP5–8.5 and SSP2–4.5 scenarios, as depicted in Figure 14(b). Specifically, the sediment yield is projected to decrease by 14.5, 21.7, 21.0, and 12.1% for the periods of 2021–2040s, 2041–2060s, 2061–2080s, and 2081–2100s, respectively, under the SSP2–4.5 scenario. Similarly, under the SSP5–8.5 scenario, the decrease in sediment yield exceeds 18.5, 35.6, 24.0, and 39.7% during the same periods. However, during the dry season, the sediment yield is projected to decline even further than during the wet season throughout the entire study period under the SSP2–4.5 and SSP5–8.5 scenarios, as shown in Figure 14(c).

The smaller decrease in sediment yield during the wet season, as compared to the dry season, can be attributed to the elevated streamflow resulting from increased precipitation levels during this specific period. In contrast, during the dry season, the reduced streamflow can be attributed to the combined effects of increased Tmax and Tmin and an accelerated evaporation rate. The rise in precipitation during the wet season enhances the overall volume of water flowing through the watersheds, leading to increased streamflow. This increased flow has a greater capacity to transport sediments, resulting in higher sediment yield during this season. Conversely, the dry season experiences higher temperatures in terms of Tmax and Tmin. These increased temperatures contribute to an elevated evaporation rate, causing the water levels in watersheds to decrease. The reduced streamflow during this period limits the transport capacity of sediments, resulting in lower sediment yield. Piao et al. (2004) stated that rising temperatures can stimulate vegetation growth. An increase in vegetation coverage can subsequently lead to a reduction in sediment yield. Therefore, the contrasting patterns of sediment yield between the wet and dry seasons can be attributed to the opposite effects of precipitation and temperature-related factors, such as streamflow and evaporation rates. These findings also indicate that climate changes significantly impact sediment yield more than streamflow. The substantial reduction in sediment yield has important implications for reservoir siltation. Reduced sediment yield can lead to slower rates of reservoir infilling, potentially impacting water storage capacity and the long-term sustainability of water resources. Additionally, the reduction in sediment supply may impact habitat quality for aquatic species and alter the geomorphology of riverbeds. This significant impact on sediment yield underscores the need for adaptive strategies in managing sediment transport and considering sediment-related aspects in climate resilience planning.

Discussion of uncertainty analysis

Future climate projections using the CMIP6 model are considered to be more reliable than those of CMIP5 (Nourani et al. 2022). However, further investigation utilizing alternative models is still necessary. Al-Hasani (2019) highlighted that streamflow is more responsive to precipitation changes than the temperature in the Mediterranean catchments of the Tigris River Basin. In contrast, our study demonstrates a clear impact of temperature on streamflow, particularly during dry seasons. Furthermore, it is worth noting that Osman et al. (2019) conducted a study specifically focused on the Greater Zab watershed, which is a part of the Tigris River Basin. In this study, the obtained changes in streamflow projections differ from the present research's findings. In Osman et al. (2019), both increases and decreases in streamflow magnitudes were projected. Particularly, the effect of temperature on streamflow was ignored in this investigation. Notably, it indicated a rise in streamflow from July to November, contrasting with the streamflow changes observed in the present research. The sources of uncertainty in the projected streamflow can be attributed to various factors, including the climate models utilized, the downscaling techniques employed, and the structure of the hydrological model.

Indeed, temperature and precipitation projections play a crucial role in assessing the uncertainties related to water resources and basin hydrology. Changes in precipitation patterns and temperature levels directly influence streamflow, water availability, sediment yield, and overall basin hydrology. Accurate predictions of these climatic variables are imperative for precisely evaluating water resources and hydrological processes. Any uncertainties in these projections can significantly impact the understanding and management of water systems (Ahmadalipour et al. 2018). This research utilized three global climate models to project future climate conditions. These models are well recognized for their accuracy and reliability in Iraq's climate research. The distribution mapping downscaling method enhanced the accuracy of climate model projections. This strategy is known for its effectiveness in addressing biases in climate models (Smitha et al. 2018; Enayati et al. 2021). By utilizing this technique, the present research intended to eliminate any potential discrepancies between the observed data and climate model results, ultimately improving projection accuracy. By considering these factors and employing ACCESS-CM2, BCC-CSM2-MR, and MRI-ESM2–0 models, along with the distribution mapping downscaling method, this study sought to reduce uncertainties and provide more accurate projections of future climate conditions related to water resources and basin hydrology. Developing suitable models and utilizing statistical downscaling methods are crucial steps in reducing the uncertainties associated with projected temperature and precipitation. These approaches will help quantify the effect of natural climate variability and anthropogenic activities on streamflow. Nevertheless, due to the limitations of this study, further research is necessary to explore these dimensions. Additional analyses can enhance our understanding of specific impacts and refine projections, ultimately reducing uncertainties in assessing water resources, streamflow, and sediment yield dynamics.

Changes in land use and human activities, such as urban expansion and reservoir operations, significantly impact runoff generation, sediment yield, and water resource availability (Zhou et al. 2015). This investigation examined the potential consequences of climate change under the assumption of no changes in land use. Additionally, potential variations in dam releases resulting from modifications in the reservoir operation, which can consequently affect streamflow, were considered. The use of a calibrated streamflow model to estimate sediment yield introduces a level of uncertainty. However, employing such a model is a reasonable approach in regions where sediment yield measurements are unavailable (ungauged regions). If sediment data becomes accessible in subsequent periods, it is advisable to re-evaluate and enhance the model by incorporating it for calibration and validation. Despite these assumptions, these findings still serve as a valuable reference for understanding the potential effect of future climate change on streamflow and sediment yield.

The primary objective of this research was to investigate and evaluate how potential climate change would affect streamflow and sediment yield in the MMB that was located in northern Iraq during the 21st century. The study also aimed to evaluate the possible changes in streamflow and sediment yield under two distinct climate scenarios. A SWAT-based model was employed to simulate and assess streamflow and sediment yield in the basin. The sensitivity analysis and model performance were evaluated using SWAT-CUP software. The CMIP6 was queried for precipitation and temperature projections in the future under two scenarios (SSP2–4.5 and SSP5–8.5). A distribution mapping bias correction method was utilized to enhance the precision of climate data.

In the late 21st century, the MMB is projected to experience a decline in average annual precipitation. Specifically, the projected decrease is approximately 9.5% in the SSP2–4.5 scenario, while under the SSP5–8.5 scenario, it is anticipated to decrease by around 18.7%. The results of this study demonstrate that the rate of increase in minimum temperature exceeds that of maximum temperature. Furthermore, the changes in both minimum and maximum temperatures during the summer, particularly in July and August, were more substantial compared to the rainy and winter seasons. The combined effects of reduced precipitation and increased temperatures synergistically influence the streamflow. According to the SSP2–4.5 scenario, the study projects a significant reduction in streamflow of more than 36.6% between 2021 and 2100. However, under the more severe SSP5–8.5 scenario, the streamflow exhibited an even higher rate of decline, with a reduction of 45.9% over the same period. It is worth noting that regardless of the emission scenario, there was a notable change in the pattern of peak flow. The peak flow will shift from May to April under both emission scenarios. The sediment yield decreased by 46.0% from 2021 to 2100 during the SSP2–4.5 scenario and decreased by 55.0% during the SSP5–8.5 scenario over the same period. The monthly analysis revealed that the runoff and sediment yield increases were more noticeable during the wet seasons but decreased during the dry season. The findings also indicate that climate changes significantly impact sediment yield more than streamflow.

The projected decrease in streamflow and sediment yield due to climate change could lead to challenges in maintaining sufficient reservoir levels for consistent water supply, especially in regions relying on dams for their water needs. Reduced sediment inflow into the reservoir can influence water quality by trapping and removing nutrients, pollutants, and other substances from the water. A decrease in sediment supply may lead to changes in water quality, which can affect the ecological health of the reservoir and the suitability of the water for various uses.

The findings of this paper provide valuable insights into how future climate scenarios may impact streamflow and sediment yield in the MMB, serving as a pre-dam construction study. Decision-makers could use this information to develop and implement proactive water resource management strategies in the MMB and the Makhool Dam. This will help in mitigating potential future challenges related to water resources. Decision-makers could develop adaptive strategies to address the potential impacts of climate change on water resources of more frequent droughts and changing precipitation patterns. This paper can serve as a useful reference for future studies focusing on the design and operation of the dams. Indeed, in addition to considering climate change impacts, it is essential to account for other factors, such as LULC changes. Incorporating these factors into future research efforts will provide a more comprehensive understanding of the dynamics affecting water resources in the region.

For their helpful and encouragement, the authors would like to extend their gratitude to the Civil Engineering Department professors and employees at the University of Technology in Iraq.

This research received no external funding.

Data cannot be made publicly available; readers should contact the corresponding author for details.

The authors declare there is no conflict.

Abbas
N.
,
Wasimia
S. A.
&
Al-Ansari
N.
2016
Assessment of climate change impacts on water resources of Khabour in Kurdistan, Iraq using SWAT model
.
J. Environ. Hydrol.
24
,
1
21
.
Abbas
M.
,
Zhao
L.
&
Wang
Y.
2022
Perspective impact on water environment and hydrological regime owing to climate change: A review
.
Hydrology
9
,
203
.
https://doi.org/10.3390/hydrology9110203
.
Abou Rafee
S. A.
,
Uvo
C. B.
,
Martins
J. A.
,
Domingues
L. M.
,
Rudke
A. P.
,
Fujita
T.
&
Freitas
E. D.
2019
Large-scale hydrological modelling of the Upper Paraná River Basin
.
Water
11
,
882
.
https://doi.org/10.3390/w11050882
.
Ahmadalipour
A.
,
Moradkhani
H.
&
Rana
A.
2018
Accounting for downscaling and model uncertainty in fine-resolution seasonal climate projections over the Columbia River Basin
.
Clim. Dyn.
50
,
717
733
.
https://doi.org/10.1007/s00382-017-3639-4
.
Ahmed
K.
,
Shahid
S.
,
Wang
X.
,
Nawaz
N.
&
Khan
N.
2019
Spatiotemporal changes in aridity of Pakistan during 1901–2016
.
Hydrol. Earth Syst. Sci.
23
,
3081
3096
.
https://doi.org/10.5194/hess-23-3081-2019
.
Al-Hasani
A. A. J.
2019
Sensitivity assessment of the impacts of climate change on streamflow using climate elasticity in Tigris River Basin, Iraq
.
Int. J. Environ. Stud.
76
,
7
28
.
https://doi.org/10.1080/00207233.2018.1494924
.
Allan
R. P.
,
Hawkins
E.
,
Bellouin
N.
&
Collins
B.
2021
IPCC, 2021: Summary for Policymakers
.
https://doi.org/10.1017/9781009157896.001
.
Arnold
J. G.
,
Srinivasan
R.
,
Muttiah
R. S.
&
Williams
J. R.
1998
Large area hydrologic modeling and assessment part I: Model development
.
J. Am. Water Resour. Assoc.
34
,
73
89
.
https://doi.org/10.1111/j.1752-1688.1998.tb05961.x
.
Arnold
J. G.
,
Moriasi
D. N.
,
Gassman
P. W.
,
Abbaspour
K. C.
,
White
M. J.
,
Srinivasan
R.
,
Santhi
C.
,
Harmel
R. D.
,
Van Griensven
A.
&
Van Liew
M. W.
2012
SWAT: Model use, calibration, and validation
.
Trans. ASABE
55
,
1491
1508
.
Azari
M.
,
Moradi
H. R.
,
Saghafian
B.
&
Faramarzi
M.
2016
Climate change impacts on streamflow and sediment yield in the North of Iran
.
Hydrol. Sci. J.
61
,
123
133
.
https://doi.org/10.1080/02626667.2014.967695
.
Bhatta
B.
,
Shrestha
S.
,
Shrestha
P. K.
&
Talchabhadel
R.
2019
Evaluation and application of a SWAT model to assess the climate change impact on the hydrology of the Himalayan River Basin
.
Catena
181
,
104082
.
https://doi.org/10.1016/j.catena.2019.104082
.
de Oliveira
V. A.
,
de Mello
C. R.
,
Beskow
S.
,
Viola
M. R.
&
Srinivasan
R.
2019
Modeling the effects of climate change on hydrology and sediment load in a headwater basin in the Brazilian Cerrado biome
.
Ecol. Eng.
133
,
20
31
.
https://doi.org/10.1016/j.ecoleng.2019.04.021
.
Dhanesh
Y.
,
Bindhu
V. M.
,
Senent-Aparicio
J.
,
Brighenti
T. M.
,
Ayana
E.
,
Smitha
P. S.
,
Fei
C.
&
Srinivasan
R.
2020
A comparative evaluation of the performance of CHIRPS and CFSR data for different climate zones using the SWAT model
.
Remote Sens.
12
,
3088
.
https://doi.org/10.3390/rs12183088
.
dos Santos
J. Y. G.
,
Montenegro
S. M. G. L.
,
da Silva
R. M.
,
Santos
C. A. G.
,
Quinn
N. W.
,
Dantas
A. P. X.
&
Neto
A. R.
2021
Modeling the impacts of future LULC and climate change on runoff and sediment yield in a strategic basin in the Caatinga/Atlantic forest ecotone of Brazil
.
Catena
203
,
105308
.
https://doi.org/10.1016/j.catena.2021.105308
.
Duan
Z.
,
Tuo
Y.
,
Liu
J.
,
Gao
H.
,
Song
X.
,
Zhang
Z.
,
Yang
L.
&
Mekonnen
D. F.
2019
Hydrological evaluation of open-access precipitation and air temperature datasets using SWAT in a poorly gauged basin in Ethiopia
.
J. Hydrol.
569
,
612
626
.
https://doi.org/10.1016/j.jhydrol.2018.12.026
.
Enayati
M.
,
Bozorg-Haddad
O.
,
Bazrafshan
J.
,
Hejabi
S.
&
Chu
X.
2021
Bias correction capabilities of quantile mapping methods for rainfall and temperature variables
.
J. Water Clim. Change
12
,
401
419
.
https://doi.org/10.2166/wcc.2020.261
.
Gould
G. K.
,
Liu
M.
,
Barber
M. E.
,
Cherkauer
K. A.
,
Robichaud
P. R.
&
Adam
J. C.
2016
The effects of climate change and extreme wildfire events on runoff erosion over a mountain watershed
.
J. Hydrol.
536
,
74
91
.
https://doi.org/10.1016/j.jhydrol.2016.02.025
.
Hamed
M. M.
,
Nashwan
M. S.
,
Shahid
S.
,
Bin Ismail
T.
,
Wang
X.
,
Dewan
A.
&
Asaduzzaman
M.
2022a
Inconsistency in historical simulations and future projections of temperature and rainfall: A comparison of CMIP5 and CMIP6 models over Southeast Asia
.
Atmos. Res.
265
,
105927
.
https://doi.org/10.1016/j.atmosres.2021.105927
.
Hamed
M. M.
,
Nashwan
M. S.
,
Shiru
M. S.
&
Shahid
S.
2022b
Comparison between CMIP5 and CMIP6 models over MENA region using historical simulations and future projections
.
Sustainability
14
,
10375
.
https://doi.org/10.3390/su141610375
.
Hamed
M. M.
,
Sammen
S. S.
,
Nashwan
M. S.
&
Shahid
S.
2023
Spatiotemporal variation of drought in Iraq for shared socioeconomic pathways
.
Stoch. Environ. Res. Risk Assess.
37
,
1321
1331
.
https://doi.org/10.1007/s00477-022-02343-7
.
Kamruzzaman
M.
,
Shahid
S.
,
Islam
A. R. M. T.
,
Hwang
S.
,
Cho
J.
,
Zaman
M. A. U.
,
Ahmed
M.
,
Rahman
M. M.
&
Hossain
M. B.
2021
Comparison of CMIP6 and CMIP5 model performance in simulating historical precipitation and temperature in Bangladesh: A preliminary study
.
Theor. Appl. Climatol.
145
,
1385
1406
.
https://doi.org/10.1007/s00704-021-03691-0
.
Katsanos
D.
,
Retalis
A.
&
Michaelides
S.
2016
Validation of a high-resolution precipitation database (CHIRPS) over Cyprus for a 30-year period
.
Atmos. Res.
169
,
459
464
.
https://doi.org/10.1016/j.atmosres.2015.05.015
.
Lehner
F.
,
Wood
A. W.
,
Vano
J. A.
,
Lawrence
D. M.
,
Clark
M. P.
&
Mankin
J. S.
2019
The potential to reduce uncertainty in regional runoff projections from climate models
.
Nat. Clim. Change
9
,
926
933
.
https://doi.org/10.1038/s41558-019-0639-x
.
Mengistu
A. G.
,
Woldesenbet
T. A.
&
Dile
Y. T.
2022
Evaluation of observed and satellite-based climate products for hydrological simulation in data-scarce Baro-Akob River Basin, Ethiopia
.
Ecohydrol. Hydrobiol.
22
,
234
245
.
https://doi.org/10.1016/j.ecohyd.2021.11.006
.
Moriasi
D. N.
,
Arnold
J. G.
,
Van Liew
M. W.
,
Bingner
R. L.
,
Harmel
R. D.
&
Veith
T. L.
2007
Model evaluation guidelines for systematic quantification of accuracy in watershed simulations
.
Trans. ASABE
50
,
885
900
.
https://doi.org/10.13031/2013.23153
.
Nashwan
M. S.
&
Shahid
S.
2020
A novel framework for selecting general circulation models based on the spatial patterns of climate
.
Int. J. Climatol.
40
,
4422
4443
.
https://doi.org/10.1002/joc.6465
.
Natumanya
E.
,
Ribeiro
N.
,
Mwanjalolo
M. J. G.
&
Steinbruch
F.
2022
Using SWAT model and field data to determine potential of NASA-POWER data for modelling rainfall-runoff in Incalaue River Basin
.
Comput. Water Energy Environ. Eng.
11
,
65
83
.
https://doi.org/10.4236/cweee.2022.112004
.
Neitsch
S. L.
,
Arnold
J. G.
,
Kiniry
J. R.
&
Williams
J. R.
2011
Soil and Water Assessment Tool Theoretical Documentation Version 2009
.
Texas Water Resources Institute
,
College Station, TX
. .
Nourani
V.
,
Paknezhad
N. J.
&
Huang
J. J.
2022
Application of PPIE method to assess the uncertainty and accuracy of multi-climate model-based temperature and precipitation downscaling
.
Theor. Appl. Climatol.
147
,
1327
1343
.
https://doi.org/10.1007/s00704-021-03884-7
.
Oguntunde
P. G.
&
Abiodun
B. J.
2013
The impact of climate change on the Niger River Basin hydroclimatology, West Africa
.
Clim. Dyn.
40
,
81
94
.
https://doi.org/10.1007/s00382-012-1498-6
.
O'Neill
B. C.
,
Kriegler
E.
,
Ebi
K. L.
,
Kemp-Benedict
E.
,
Riahi
K.
,
Rothman
D. S.
,
Van Ruijven
B. J.
,
Van Vuuren
D. P.
,
Birkmann
J.
&
Kok
K.
2017
The roads ahead: Narratives for shared socioeconomic pathways describing world futures in the 21st century
.
Glob. Environ. Change
42
,
169
180
.
https://doi.org/10.1016/j.gloenvcha.2015.01.004
.
Osman
Y.
,
Abdellatif
M.
,
Al-Ansari
N.
,
Knutsson
S.
&
Jawad
S.
2017
Climate change and future precipitation in an arid environment of the Middle East: Case study of Iraq
.
J. Environ. Hydrol.
25
,
1
18
.
Osman
Y.
,
Al-Ansari
N.
&
Abdellatif
M.
2019
Climate change model as a decision support tool for water resources management in northern Iraq: A case study of Greater Zab River
.
J. Water Clim. Change
10
,
197
209
.
https://doi.org/10.2166/wcc.2017.083
.
Parmesan
C.
,
Morecroft
M. D.
&
Trisurat
Y.
2022
Climate Change 2022: Impacts, Adaptation and Vulnerability
.
IPCC
,
Geneva, Switzerland
.
https://hal.science/hal-03774939/document (accessed 29 December 2023)
.
Piao
S.
,
Fang
J.
,
Ji
W.
,
Guo
Q.
,
Ke
J.
&
Tao
S.
2004
Variation in a satellite-based vegetation index in relation to climate in China
.
J. Veg. Sci.
15
,
219
226
.
https://doi.org/10.1658/1100-9233(2004)015[0219:VIASVI]2.0.CO;2
.
Pour
S. H.
,
Shahid
S.
&
Mainuddin
M.
2022
Relative performance of CMIP5 and CMIP6 models in simulating rainfall in Peninsular Malaysia
.
Theor. Appl. Climatol.
149
,
709
725
.
https://doi.org/10.1007/s00704-022-04076-7
.
Reddy
N. M.
,
Saravanan
S.
,
Almohamad
H.
,
Al Dughairi
A. A.
&
Abdo
H. G.
2023
Effects of climate change on streamflow in the Godavari basin simulated using a conceptual model including CMIP6 dataset
.
Water
15
,
1701
.
https://doi.org/10.3390/w15091701
.
Romshoo
S. A.
&
Marazi
A.
2022
Impact of climate change on snow precipitation and streamflow in the Upper Indus Basin ending twenty-first century
.
Clim. Change
170
,
6
.
https://doi.org/10.1007/s10584-021-03297-5
.
Ruelland
E.
,
Kravets
V.
,
Derevyanchuk
M.
,
Martinec
J.
,
Zachowski
A.
&
Pokotylo
I.
2015
Role of phospholipid signalling in plant environmental responses
.
Environ. Exp. Bot.
114
,
129
143
.
https://doi.org/10.1016/j.envexpbot.2014.08.009
.
Salehie
O.
,
Hamed
M. M.
,
Ismail
T. B.
&
Shahid
S.
2022
Projection of droughts in Amu river basin for shared socioeconomic pathways CMIP6
.
Theor. Appl. Climatol.
149
,
1009
1027
.
https://doi.org/10.1007/s00704-022-04097-2
.
Salman
S. A.
,
Shahid
S.
,
Ismail
T.
,
Chung
E.-S.
&
Al-Abadi
A. M.
2017
Long-term trends in daily temperature extremes in Iraq
.
Atmos. Res.
198
,
97
107
.
https://doi.org/10.1016/j.atmosres.2017.08.011
.
Salman
S. A.
,
Shahid
S.
,
Ismail
T.
,
Ahmed
K.
&
Wang
X.-J.
2018
Selection of climate models for projection of spatiotemporal changes in temperature of Iraq with uncertainties
.
Atmos. Res.
213
,
509
522
.
https://doi.org/10.1016/j.atmosres.2018.07.008
.
Salman
S. A.
,
Shahid
S.
,
Ismail
T.
,
Al-Abadi
A. M.
,
Wang
X.
&
Chung
E.-S.
2019
Selection of gridded precipitation data for Iraq using compromise programming
.
Measurement
132
,
87
98
.
https://doi.org/10.1016/j.measurement.2018.09.047
.
Salman
S. A.
,
Hamed
M. M.
,
Shahid
S.
,
Ahmed
K.
,
Sharafati
A.
,
Asaduzzaman
M.
,
Ziarh
G. F.
,
Ismail
T.
,
Chung
E.
&
Wang
X.
2022
Projecting spatiotemporal changes of precipitation and temperature in Iraq for different shared socioeconomic pathways with selected Coupled Model Intercomparison Project Phase 6
.
Int. J. Climatol.
https://doi.org/10.1002/joc.7794
.
Schilling
J.
,
Hertig
E.
,
Tramblay
Y.
&
Scheffran
J.
2020
Climate change vulnerability, water resources and social implications in North Africa
.
Reg. Environ. Change
20
,
1
12
.
https://doi.org/10.1007/s10113-020-01597-7
.
Schwinning
S.
,
Sala
O. E.
,
Loik
M. E.
&
Ehleringer
J. R.
2004
Thresholds, memory, and seasonality: Understanding pulse dynamics in arid/semi-arid ecosystems
.
Oecologia
.
https://doi.org/10.1007/s00442-004-1683-3
.
Sharafati
A.
,
Pezeshki
E.
,
Shahid
S.
&
Motta
D.
2020
Quantification and uncertainty of the impact of climate change on river discharge and sediment yield in the Dehbar river basin in Iran
.
J. Soils Sediments
20
,
2977
2996
.
https://doi.org/10.1007/s11368-020-02632-0
.
Shrestha
S.
,
Bhatta
B.
,
Talchabhadel
R.
&
Virdis
S. G. P.
2022
Integrated assessment of the land use change and climate change impacts on the sediment yield in the Songkhram River Basin, Thailand
.
Catena
209
,
105859
.
https://doi.org/10.1016/j.catena.2021.105859
.
Siabi
E. K.
,
Awafo
E. A.
,
Kabo-bah
A. T.
,
Derkyi
N. S. A.
,
Akpoti
K.
,
Mortey
E. M.
&
Yazdanie
M.
2023
Assessment of Shared Socioeconomic Pathway (SSP) climate scenarios and its impacts on the Greater Accra region
.
Urban Clim.
49
,
101432
.
https://doi.org/10.1016/j.uclim.2023.101432
.
Sissakian
V. K.
,
Al-Ansari
N.
,
Adamo
N.
,
Laue
J.
&
Knutsson
S.
2018
Geology of the Tigris River with emphasize on the Iraqi part
.
J. Earth Sci. Geotech. Eng.
8
,
145
166
.
Smitha
P. S.
,
Narasimhan
B.
,
Sudheer
K. P.
&
Annamalai
H.
2018
An improved bias correction method of daily rainfall data using a sliding window technique for climate change impact assessment
.
J. Hydrol.
556
,
100
118
.
https://doi.org/10.1016/j.jhydrol.2017.11.010
.
Stocker
T.
,
Qin
D.
,
Plattner
G.-K.
,
Tignor
M.
,
Allen
S.
,
Boschung
J.
,
Nauels
A.
,
Xia
Y.
,
Bex
V.
&
Midgley
P.
2014
Summary for Policymakers
.
Tian
P.
,
Lu
H.
,
Feng
W.
,
Guan
Y.
&
Xue
Y.
2020
Large decrease in streamflow and sediment load of Qinghai–Tibetan Plateau driven by future climate change: A case study in Lhasa River Basin
.
Catena
187
,
104340
.
https://doi.org/10.1016/j.catena.2019.104340
.
Zamani
Y.
,
Hashemi Monfared
S. A.
,
Azhdari Moghaddam
M.
&
Hamidianpour
M.
2020
A comparison of CMIP6 and CMIP5 projections for precipitation to observational data: The case of Northeastern Iran
.
Theor. Appl. Climatol.
142
,
1613
1623
.
https://doi.org/10.1007/s00704-020-03406-x
.
Zhang
Y.
,
You
Q.
,
Chen
C.
&
Ge
J.
2016
Impacts of climate change on streamflows under RCP scenarios: A case study in Xin River Basin, China
.
Atmos. Res.
178
,
521
534
.
https://doi.org/10.1016/j.atmosres.2016.04.018
.
Zhang
L.-X.
,
Xiao-Long
C.
&
Xiao-Ge
X. I. N.
2019
Short commentary on CMIP6 scenario model intercomparison project (ScenarioMIP)
.
Adv. Clim. Change Res.
15
,
519
.
Zhou
J.
,
He
D.
,
Xie
Y.
,
Liu
Y.
,
Yang
Y.
,
Sheng
H.
,
Guo
H.
,
Zhao
L.
&
Zou
R.
2015
Integrated SWAT model and statistical downscaling for estimating streamflow response to climate change in the Lake Dianchi watershed, China
.
Stoch. Environ. Res. Risk Assess.
29
,
1193
1210
.
https://doi.org/10.1007/s00477-015-1037-1
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).