Despite their global significance, the impacts of climate change on water
resources and associated ecosystem services in the Congo River basin (CRB)
have been understudied. Of particular need for decision makers is the
availability of spatial and temporal variability of runoff projections.
Here, with the aid of a spatially explicit hydrological model forced with
precipitation and temperature projections from 25 global climate models (GCMs)
under two greenhouse gas emission scenarios, we explore the
variability in modeled runoff in the near future (2016–2035) and mid-century (2046–2065).
We find that total runoff from the CRB is projected to increase by
5 % [
Sustainable management of water resources for food production, supply of safe drinking water and provision of adequate sanitation present immense challenges in many countries of central Africa where the Congo River basin (CRB) is located (IPCC, 2014; UNEP, 2011; World Food Program, 2014). The economies of the nine countries that share the waters of the CRB are agriculture based (World Bank Group, 2014) and therefore are vulnerable to the impacts of climate change. Despite the abundant water and land resources and favorable climates, the basin countries are net importers of staple food grains and are far behind in achieving Millennium Development Goals (Bruinsma, 2003; Molden, 2007; UNEP, 2011). Appropriation of freshwater resources is expected to grow in the future as the CRB countries develop and expand their economies. At the same time, climate-change-related risks associated with water resources will also increase significantly (IPCC, 2014).
Historical, present and near-future greenhouse gas emissions in the CRB countries constitute a small fraction of global emissions; however, the impacts of climate change on water resources are expected to be severe due to the region's heavy reliance on natural resources (e.g., agriculture and forestry) (Collier et al., 2008; DeFries and Rosenzweig, 2010; Niang et al., 2014). The limited adaptation capacity in the CRB region is expected to cause water and food security challenges, which, in turn, can lead to ecosystem degradation and increased greenhouse gas emissions (Gibbs et al., 2010; IPCC, 2014; Malhi and Grace, 2000).
Strategies for addressing stresses on CRB water resources, including revival of rural economies (largely agriculture based), achieving millennium development goals and environmental conservation, would benefit from detailed information on the spatial and temporal variability of water balance components under different climate projection pathways. The effect of climate change on water resources can be investigated by incorporating climate change projections (e.g., precipitation and temperature) in simulation models that reliably represent the spatial and temporal variability of the CRB's hydrology. Such a framework could be applied to project changes in storage and runoff, and hence freshwater availability, under different socioeconomic pathways that affect climate trajectories.
A predictive framework of the CRB's hydrology is hindered by insufficient data and too few evaluations of models against available data (Beighley et al., 2011; Wohl et al., 2012). Basin-scale water budgets estimated from land-based and satellite-derived precipitation datasets reveal significantly different results, and modeled runoff shows only qualitative agreement with corresponding observations (Alsdorf et al., 2016; Beighley et al., 2011; Lee et al., 2011; Schuol et al., 2008). Tshimanga and Hughes (2012, 2014) recently developed a semi-distributed hydrologic model capable of simulating runoff in CRB. This work crucially identified approaches suitable for approximating runoff generation at the basin scale, although the spatial resolution of the model predictions is rather coarse for supporting regional water management and regional planning efforts. These regional planning efforts must take into account variability and uncertainties stemming from climate model selection and projected greenhouse gas emissions, but, with respect to freshwater runoff projections for the CRB, these issues have been inadequately addressed.
The goals of this study are to (i) develop a spatially explicit hydrology model that uses downscaled output from general circulation models (GCMs) and is suitable for simulating the spatiotemporal variability of runoff in the CRB; (ii) test the ability of the hydrological model to reproduce historical data on CRB river discharges using both observed and GCM-simulated climate fields; (iii) quantify the sensitivity and uncertainty of modeled runoff projections to GCM selection; (iv) use the hydrologic model with individual GCMs and multi-GCM ensembles to project near-term (2016–2035) and midterm (2046–2065) changes in runoff for two greenhouse gas emission scenarios. We focus on the runoff projections because streams and rivers will serve as the primary sources of freshwater targeted for human appropriation (Burney et al., 2013; Molden, 2007).
Congo River basin: the river basin boundary, the extent of the rainforest, locations of lakes and wetlands, and the locations of streamflow gages are shown. The “all other vegetation” category includes grasslands and savanna ecosystems, and all managed areas. Bartholomé and Belward (2005) provide further details on land cover in the Congo River basin.
The Congo River basin, with a drainage area of 3.7 million km
In order to compare regional patterns in precipitation and runoff, we divided the basin into four regions: (i) northern Congo (NC), (ii) equatorial Congo (EQ), (iii) southwestern Congo (SW) and (iv) southeastern Congo (SE). The EQ region covers most of the rainforest. The SE region consists of numerous interconnected lakes and wetlands. Most of the CRB's population is concentrated in the NC, SE and SW regions (Center for International Earth Science Information Network (CIESIN) Columbia University et al., 2005).
We used the Soil Water Assessment Tool (SWAT), a physically based, semi-distributed watershed-scale model that operates at a daily time step (Arnold et al., 1998; Neitsch et al., 2011). The hydrological processes simulated include evapotranspiration, infiltration, surface and subsurface flows, streamflow routing and groundwater recharge. The model has been successfully employed to simulate river basin hydrology under a wide variety of conditions and to investigate climate change effects on water resources (Faramarzi et al., 2013; Krysanova and White, 2015; Schuol et al., 2008; Trambauer et al., 2013; van Griensven et al., 2012).
We delineated 1575 watersheds within the CRB based on topography (Lehner et al., 2008).
Watershed elevations varied between 15 and 2700 m with a mean value of
680 m a.m.s.l. (above mean sea level). Each watershed consisted of one stream section,
where near-surface groundwater flow and overland flow accumulated before
being transmitted through the stream channel to the watershed outlet.
Watersheds were further divided into hydrologic response units (HRUs) based
on land cover (16 classes; Bartholomé and Belward, 2005), soils (150 types; FAO/IIASA, 2009) and
topography. The runoff generated within each watershed was routed through
the stream network using the variable storage routing method. The average
watershed size and the number of HRUs within each watershed were 2300 km
Simulated runoff, estimated for each HRU and aggregated at the watershed level, was generated via three pathways: overland flow, lateral subsurface flow through the soil zone and release from shallow groundwater storage. The curve number and a kinematic storage routing method were used to simulate overland and lateral subsurface flows, and a nonlinear storage–discharge relationship was used to simulate groundwater contribution (see Arnold et al., 1998; Neitsch et al., 2011 and the Supplement). A power-law relationship was employed to simulate the lake area–volume discharge (see the Supplement and Neitsch et al., 2011). The potential evapotranspiration was estimated using the temperature-based Hargreaves method (Neitsch et al., 2011). The actual evapotranspiration was estimated based on available soil moisture and the evaporative demand (i.e., potential evapotranspiration) for the day. Additional details on model development and calibration are provided in the Supplement.
We ran the hydrology model for the period 1950–2008. Estimates of observed daily precipitation, and minimum and maximum temperatures needed to calculate potential evapotranspiration were obtained from the Land Surface Hydrology Group at Princeton University (Sheffield et al., 2006). In addition, measured monthly streamflows were obtained at 30 gage locations (Fig. 1) that had at least 10 years of records (Global Runoff Data Center, 2011; Lempicka, 1971; Vorosmarty et al., 1998).
The model was calibrated using observed streamflows for the period 1950–1957 at 20 locations. The locations of streamflow gages and time period were chosen such that they adequately captured climatic, land cover and topographic variability within the CRB. The number of model parameters estimated by calibration varied from 10 to 13, depending on the location of flow gages (e.g., gages with lakes within their catchment area have more parameters). The calibration involved minimizing an objective function defined as the sum of squared errors between observed and simulated monthly average total discharge, baseflows (estimated by the baseflow separation method of Nathan and McMahon, 1990) and water yield. The Gauss–Marquardt–Levenberg algorithm, as implemented in a model-independent parameter estimation tool (Doherty, 2004), was used to adjust the fitted parameters and minimize the objective function. Parameter estimation was done in two stages. First, parameters for the watersheds in the upstream gages were estimated. Then, the parameters for the downstream gages were estimated. To test the calibrated model, simulated streamflows were compared to streamflows measured at the same 20 locations but during a period outside of calibration (i.e., 1958–2008), as well as at 10 additional locations that were not used in the calibration.
Historical climate simulations for the period 1950–2005 and climate
projections to 2065 for two greenhouse gas emission scenarios
(Representative Concentration Pathways – RCPs), mid-range mitigation emission (RCP4.5)
and high emission (RCP8.5), were used to drive the hydrologic
model. The RCP4.5 scenario employs a range of technologies and policies that
reduce greenhouse gas emissions and stabilize radiative forcing at 4.5 W m
Global climate models whose outputs are used in this study. Further details about comparison of model outputs and key references for GCMs are given in Aloysius et al. (2016).
GCM outputs may exhibit biases in simulating regional climate. These biases, which are attributable to inadequate representation of physical processes by the models, prevent the direct use of GCM output in climate change studies (Randall et al., 2007; Salathé Jr. et al., 2007; Wood et al., 2004). Hydrological assessments that use GCM computations as input inherit the biases (Salathé Jr. et al., 2007; Teutschbein and Seibert, 2012). To mitigate this problem, we implemented a statistical method (Li et al., 2010) to bias correct the monthly historical precipitation and temperature data. In brief, the method employs a quantile-based mapping of cumulative probability density functions for monthly GCM outputs onto those of gridded observations in the historical period. The bias correction is extended to future projections as well. The observed data used in the modeling and bias correction have some limitations. That is, the number of precipitation gages decreased over the period from 1950 to 1990, and the density of the gages is sparse compared to the size of the river basin (see Sect. 3.4 and the Supplement). However, we assumed that the available ground-based observations combined with satellite-based and reanalysis data adequately captured the spatiotemporal variability in precipitation. Studies by Munzimi et al. (2014) and Nicholson (2000) draw similar conclusions.
The simulated monthly precipitation and temperature values were temporally
downscaled to daily values for use in the CRB hydrology model. We used the
3-hourly and monthly observed historical data developed for the Global
Land Data Assimilation System (Rodell et al., 2004; Sheffield et al., 2006) and the bias-corrected
monthly simulations to generate 3-hourly precipitation and temperature
data, which were subsequently aggregated to obtain daily values (see
the Supplement). The hydrological model was forced with the bias-corrected and
downscaled daily climate for the period 1950–2065. Due to the lack of
information on the effect of CO
A total of 50 projections (25 RCP4.5 and 25 RCP8.5 projections; see Table 1) were compiled and analyzed. Results of individual and multimodel means (unweighted average of all models (MMs) and an average of select models – SMs) for the near-term (2016–2035) and midterm (2046–2065) projections are presented.
Accessible flows (AFs), which exclude surface runoff associated the storm events, were estimated by applying a baseflow separation method described in Nathan and McMahon (1990).
Historical observations of average annual precipitation vary from 1100 mm
in the southeastern portion of the CRB to 1600 mm in the CRB's equatorial
region. We compared the GCM-simulated annual precipitation and its
interannual variability during the historical period with observations from
30 locations within the CRB (Fig. 2). The simulated interannual
variability among the climate models (vertical bars in Fig. 2) lies within
the range of the observed variability (horizontal bars in Fig. 2). The
linear regression slope of 1.16 (
Comparison of observed and bias-corrected GCM-simulated average
annual precipitation for 30 catchments with streamflow gages (shown in
Fig. 1) in the historical period (1950–2005). The
Comparison of observed and simulated annual runoff at the 30 streamflow
gage locations (shown in Fig. 1).
We compared the simulated monthly runoff at 30 locations with observations
(Fig. 3a and Table S3). The colored points compare observed mean annual
runoff at the 30 gage locations with historical simulations (hydrological
model forced with observed climate), while the vertical and horizontal bars
show the modeled and observed interannual variability, respectively. The
shades of colors (from light green to yellow and red) reveal the model's
skill in simulating the monthly flows in the historical period. The
Nash–Sutcliffe coefficient of efficiency (NSE), a measure of relative
magnitude of residual variance compared to the monthly observed streamflow
variance (Legates and McCabe Jr., 1999; Nash and Sutcliffe, 1970), varies between 0.01 and 0.86 (color scale in
Fig. 3a). The NSE can range from negative infinity to 1, with values
between 0.5 and 1 considered satisfactory (Moriasi et al., 2007). A total of 17 of the
30 gages show NSE greater than or equal to 0.5. Higher NSE values at locations
on both sides of the Equator, particularly at major tributaries (NSE
Comparison of modeled runoff forced with GCM-simulated and observed climate (Fig. 3b) reveals generally acceptable runoff simulations in the CRB. The black dots and red (blue) vertical bars in Fig. 3b show multimodel mean and maximum (minimum) range of interannual variability in the 25 historical GCM simulations. The results suggest that model–data agreement in precipitation translates to similarly acceptable runoff simulations.
Runoff patterns reflect seasonal rainfall that varies asymmetrically on either side of the Equator (see Fig. S1). For example, the observed peak runoff at streamflow gages 2 and 6 located north and south of the Equator (see Fig. 1) occur near the end of the rainy seasons – during September–October and March–April, respectively (Fig. 4). Augmented by flows from northern and southern tributaries (e.g., gages 1, 2, 4 and 6) and by high precipitation in the tropical equatorial watersheds during the two wet seasons (MAM and SON), the main river flows (downstream of gage 3 in Fig. 1) show low variability (Fig. 4). Differences in streamflow variability between the main river and its tributaries are illustrated through comparison of the coefficient of variation, which equals only 0.23 at the basin outlet (gage 8), but is 0.77 and 0.40 along the northern tributary (gage 2) and southern tributary (gage 4), respectively.
Runoff in the northern (NC) and southern (SW and SE) watersheds is strongly seasonal with long dry seasons, but this is not the case in the equatorial region (Fig. 5). Average watershed runoff varies between 20 and 70 mm during dry seasons to 100–140 mm during wet seasons in the NC, SW and SE. In the equatorial region, seasonal runoff varies between 100 and 150 mm with the highest values in SON. Overall, the precipitation–runoff ratio is about 0.30 in the CRB. The AF that can be appropriated for human use, and hence excludes runoff associated with flood events, is about 70 % of the total runoff.
Aloysius et al. (2016) showed that GCM projections of temperature generally increase under both
emission scenarios, in line with the historical upward trend for Africa
(Hulme et al., 2001); however, precipitation projections contain large uncertainties.
The modeled near-term (2016–2035) precipitation projections in the CRB vary
between
MM of projected changes in precipitation (%) in the four regions within the Congo River basin (see Fig. 1) for the near term (2016–2035) and the midterm (2046–2065) relative to the reference period of 1986–2005. The regions are identified in Fig. 1. The standard deviation values across the 25 GCM simulations are provided in parenthesis. DJF: December–January–February, MAM: March–April–May, JJA: June–July–August and SON: September–October–November.
Mean monthly flows at selected tributaries in the CRB. Flows are
in m
Seasonal variation in runoff in the
Number of climate model outputs projecting an increase in
precipitation in the
The GCMs project considerable spatial and seasonal variations in precipitation (Table 2 and Fig. 6). However, the standard deviation of annual and seasonal projections within the four regions exceed or are equal to the MM, indicating little agreement on the direction of change. The spatial patterns (Fig. 6), on the other hand, show regions where modeled projections strongly agree on increasing or decreasing precipitation. For example, decreasing precipitation is projected in most of the headwater catchments in the southern CRB and parts of northern CRB.
In general, the GCMs project decreasing precipitation in the driest parts of the southern CRB (mostly in the southeastern CRB but in portions of southwestern as well). Under the RCP8.5 scenario, parts of northeastern CRB also experience a reduction in precipitation in the near term (regions in Fig. 6 with fewer GCMs projecting an increase in precipitation). The areas of decreased precipitation shrink in the southeast and southwest in the midterm; however, drying expands in parts of northern CRB under the two emission scenarios. Most GCMs (14–20) project a precipitation increase outside of southeastern CRB.
Intermodel variability in precipitation projections are sensitive to seasons and climate region (Fig. 7a–d). At monthly scale, the northern and southern regions receive less than 50 mm of precipitation for at least 3 months, which persist in the future under both emission scenarios. The dry season is more prolonged in the southeast compared to the rest of the CRB. The intermodel variability is larger in the rainy seasons under RCP8.5 compared to RCP4.5. Larger variability under RCP8.5 highlights that GCMs may have limited skill in simulating precipitation under high greenhouse gas emissions.
In general, modeled runoff increases, and its interannual variability within GCMs is larger during high flow periods compared to low flow periods, except in the equatorial region (Fig. 7e–h; see Fig. 1 for regions). The model projection uncertainty increases towards the middle of century, particularly under the RCP8.5 emission scenario. The temporal patterns of runoff in the near term and midterm are similar to the precipitation patterns, but with a time lag. As with precipitation, the monthly runoff shows prolonged periods of low values in the northern and southern CRB in both projection periods. Parts of northern, southeastern and southwestern CRB also show reduced runoff projections relative to the reference period under both RCPs; these reductions are predominantly in the areas where fewer GCMs agree on the increase in modeled precipitation (see Fig. 6 and Tables S4 and S5). The area of decreasing runoff expands in the northern CRB under both emission scenarios in the midterm (see Fig. 6, which shows that more models agree on decreasing precipitation in parts of northern CRB that subsequently results in decreasing runoff). Although the northern and equatorial CRB show an overall increase in precipitation, the decrease in runoff in certain parts in the northern and equatorial CRB is caused by reduction in seasonal precipitation (e.g., JJA and SON; see Table S4). A larger reduction – up to 15 % – in the southeastern CRB covering most of northern Zambia is due to an overall decrease in precipitation simulated by more than half of the GCMs (see Fig. 6).
The multimodel mean of total runoff from the CRB shows an increase of 5 %
(
Monthly variation of precipitation
Seasonal runoff projections (as percent relative to the reference
period 1986–2005) for the near-term (2016–2035) and midterm (2046-2065)
projection periods for northern
Runoff that can be appropriated for human use is generated mostly in the
northern, southeastern and southwestern CRB, which at present varies from
130 mm yr
MM of projected changes in runoff (%) in the four regions within
the Congo River basin for the near term (2016–2035) and the
midterm (2046–2065) relative to the reference period of 1986–2005. The
regions are identified in Fig. 1. The standard deviation values across the
25 GCM simulations are provided in parenthesis. The asterisks (
Accessible streamflow hydrographs in the near term at selected locations shown in Fig. 1a. Blue and red bars (RCP4.5 and RCP8.5, respectively) show the intermodel variability. The dotted black line shows the hydrograph in the reference period (1986–2005). Plot numbers 1–8 coincide with the gage numbers in Fig. 1.
Only part of the runoff may be appropriated for human use. In the CRB, the accessible runoff (AF), excluding runoff associated with flood events, is about 70 %. The AF is largely underutilized, but its appropriation is expected to increase in the future, mostly in the populated areas of northern, southwestern and southeastern CRB. We present the uncertainty associated with GCM and scenario selection by quantifying seasonal and intermodel variability in AF at eight major tributaries (identified in Fig. 1) that drain watersheds across a range of climatic regions on both sides of the Equator (Fig. 9). Modeled AF exhibits substantial intermodel spread in the near term and widens in the midterm (Fig. S4). The intermodel variability is larger during high flow periods compared to low flow periods.
Following the general pattern of increasing precipitation and runoff in the northern and southwestern watersheds, we find that AF increases with greater model agreement in tributaries that drain these watersheds (e.g., gages 1, 2 and 6 in Fig. 9). A closer look at tributaries in the northern and southwestern CRB reveals better agreement of increased AF during low flow periods compared to high flow periods (compare gages 1, 2, 6 and 7 in Fig. 9). In contrast, tributaries that drain southeastern watersheds exhibit greater variability in modeled AF with majority of the ensembles projecting a reduction (e.g., gages 4 and 5 in Fig. 7). Overall, the AF in the main tributary (gages 3 and 8) is projected to increase, partly due to the contributions from the northern and southwestern tributaries. The decrease in modeled precipitation and AF in the southeastern CRB appears to have a marginal effect on downstream flows in the main river.
The spatial and temporal variations in the projected AF have consequences for water resources development and management. For example, projections of increased AF near the proposed Grand Inga Hydropower Project (near gage 8; Showers, 2009) are robust compared to the large variations near the proposed transboundary water diversion in the southeast (near gage 5, Lund et al., 2007). Reductions in high and low flows in streams in the southeastern region will have implications for aquatic life, channel maintenance and lake and wetland flooding.
The sources of uncertainty encountered in this work can be broadly
categorized into (i) observational uncertainty, particularly the sparse and
declining network of precipitation and streamflow gages and (ii) model
uncertainty, which in the GCMs includes model structure, model
initialization, parameterization and climate sensitivity (i.e., the response
of global temperature to a doubling of CO
The climate data used for bias correction and for historical hydrologic simulations have their own uncertainties. Gage-based, satellite-derived data and reanalysis outputs are used to develop the historical observations (Sheffield et al., 2006). Precipitation gages were more numerous at the beginning of the simulation period and declined in number toward the end of the 20th century (Mitchell and Jones, 2005; Washington et al., 2013). Available gage data varied both spatially and temporally (Figs. S5 and S6). For example, the equatorial region – nearly a third of CRB – had about 70 rain gages through the early 1990s, but only 10% of these were functioning by 2005 (Fig. S5). The southeastern and parts of northern CRB also had good rainfall-gage coverage, which has similarly decreased since the 1990s (Mitchell and Jones, 2005). However, satellite-based and sparsely distributed gage data have been used to demonstrate that spatiotemporal distribution of precipitation can be sufficiently described in the CRB region (Munzimi et al., 2014; Nicholson, 2000; Samba et al., 2008). We assume that, even with these limitations, the available historical data are adequate to model the hydrology of the CRB.
In addition to climate data, observed runoff data are another limitation that could restrict proper validation of hydrological models. However, we utilized a time period (1950–1959) when the CRB had maximum coverage of both precipitation and runoff data to calibrate and test the hydrology model (for example, see evidence in L'vovich, 1979). Where available, we used additional runoff data to further test model outputs during the historical period. The runoff gage locations are distributed within the CRB (see Fig. 1) such that they adequately capture climatic, land cover and topographic variability.
For future projections, the largest sources of uncertainty arise from the GCMs and emission scenarios. GCMs do not consistently capture observed rainfall seasonality and heavy rainfall in regions of the central CRB and in most cases do not show key features such as seasonality and heavy rainfall in regions of the central CRB (Aloysius et al., 2016; Washington et al., 2013). The biases in the GCM-simulated precipitation, particularly in the tropical regions, have been attributed to multiple factors including poorly resolved physical processes such as the mesoscale convection systems, inadequately resolved topography due to the coarse horizontal resolution and inadequate observations to constrain parameterization schemes. These limitations are unavoidable in the current set of CMIP5 projections. We assume that the combination of GCM outputs used in our work and the bias-correction method, which maintains key statistical properties in the original GCM outputs (see Aloysius et al., 2016; Li et al., 2010), adequately captures the uncertainties in GCM and emission scenarios. Based on monthly precipitation climatology, Aloysius et al. (2016) found no significant shift in seasonality in modeled future precipitation projections.
The range of projections presented here for the two emission scenarios also highlight the uncertainties planners would encounter when making climate-related decisions. For example, broader agreement on increase in runoff in parts of the CRB would help make robust decisions, whereas weaker agreement in the southern CRB calls for greater scrutiny of regional climate. Generally, the MM approach reduces the uncertainty because averaging tends to offset errors across models. However, one could also ask whether this approach would work with fewer models.
Washington et al. (2013) and Siam et al. (2013) presented evidence that evaluating atmospheric moisture flux (which is modulated by wind patterns and humidity) and soil water balance is a better way to diagnose GCM performance in data-scarce regions like the CRB. Balas et al. (2007), Hirst and Hastenrath (1983) and Nicholson and Dezfuli (2013) have shown that sea surface temperature (SST) anomalies in the Atlantic and Indian ocean sectors could partly explain precipitation in the CRB region. Along the same lines, Aloysius et al. (2016) identified five models as suitable candidates. We examined this subset of GCM projections (M6, M7, M18, M23 and M24), which we refer to as the select model average, or SM (see Giorgetta et al., 2013; Good et al., 2012; Jungclaus et al., 2013; Meehl et al., 2013; Siam et al., 2013; Voldoire et al., 2012; Yukimoto et al., 2006 and Aloysius et al., 2016 for further comparison of GCM performance). By evaluating seasonal atmospheric moisture and soil water balance in 11 CMIP5 GCMs in the CRB and Nile River basin regions, Siam et al. (2013) identified M7, M18 and M24 as good candidates for climate change assessment.
Focusing on the northern, southeastern and southwestern CRB, where human appropriation of runoff is expected to increase, we find that the projected increase of annual runoff in SM is more than that of MM (20 to 50 % higher in the SM compared to MM). In addition, the extent of reduction in runoff in the south is concentrated in the southeastern upstream watersheds in both MM and SM, although the magnitude of decrease is smaller in SM (Tables S4 and S5).
From the viewpoint of water resources for human appropriation, the changes
by seasons are also important. Future changes and uncertainties in modeled
seasonal runoff averaged over the four regions are presented in Fig. 8. In
comparison with the CRB projections, the uncertainties in subregions are
larger. Nearly all the MM and SM projections show an increase in runoff in
all the four seasons; however, there is substantial intermodel variability.
The uncertainties increase under the high emission RCP8.5 scenario during
the mid-century. Considering the southeastern region as an example, under
the RCP8.5 emission scenario, uncertainties reported as 1 intermodel standard
deviation in the midterm are
From the point of view of climate change adaptation related to water resources, agriculture and ecosystem management, the challenge faced by CRB countries is recognizing the value of making timely decisions in the absence of complete knowledge. In some settings, climate change presents opportunities as well as threats in the CRB. The projected increases in accessible runoff imply new opportunities to meet increasing demands (e.g., drinking water, food production and sanitation), while the enhanced flood runoff would pose new challenges (e.g., flood protection and erosion control). On the other hand, water managers could face different challenges in the southeast where precipitation and runoff are projected to decrease.
GCM-related variability in regional climate projections could be constrained by a subset of models based on attributes that modulate large-scale circulations (see Knutti and Sedlacek, 2013; Masson and Knutti, 2011). This approach is particularly useful because regions like the CRB lack complete coverage of observational data but the mechanisms that moderate the climate system, particularly precipitation, are fairly well understood (Hastenrath, 1984; Nicholson and Grist, 2003; Washington et al., 2013). Yet, the span in rainfall predictions among the MM, SM and individual GCMs suggests that, despite the advances in climate modeling, significant uncertainties in precipitation projections for CRB persist.
Rather than providing a narrow pathway for decision-making, our results, for the first time for CRB, provide a framework to (i) assess implications under various climate model assumptions and uncertainties, (ii) characterize and expose vulnerabilities and (iii) provide ways to guide the search for impact-oriented and actionable policy alternatives, as emphasized by Weaver et al. (2013). Projections and associated uncertainties vary widely by region within the CRB, and therefore diverse but robust planning strategies might be advisable within the river basin. We emphasize that projections provided here could be considered as part of the process of incorporating multiple stressors into climate change adaptation and engaging stakeholders in the decision-making process.
Historical observed climate data, described in Sheffield et al. (2006), are available
at
The authors declare that they have no conflict of interest.
We would like to thank Nadine Laporte, Innocent Liengola, Peter Umunay, Greg Fiske and Melanie Burr for help with data and literature search. We acknowledge the World Climate Research Program's Working Group on Coupled Modeling, which is responsible for CMIP, and we thank the climate modeling groups (Table 1) for producing and making their model output available. For CMIP, the US Department of Energy's Program for Climate Model Diagnosis and Intercomparison provided coordinating support and led development of software infrastructure in partnership with the Global Organization for Earth System Science Portals. We gratefully acknowledge the efforts of three anonymous reviewers who made thoughtful comments that substantially improved the manuscript. This work was supported in part by the facilities and staff of the Yale University Faculty of Arts and Sciences High Performance Computing Center and by the National Science Foundation under grant CNS 08-21132 that partially funded acquisition of the facilities. Edited by: Lixin Wang Reviewed by: three anonymous referees