Water management substantially alters natural regimes of streamflow through modifying retention time and water exchanges among different components of the terrestrial water cycle. Accurate simulation of water cycling in intensively managed watersheds, such as the Yakima River basin (YRB) in the Pacific Northwest of the US, faces challenges in reliably characterizing influences of management practices (e.g., reservoir operation and cropland irrigation) on the watershed hydrology. Using the Soil and Water Assessment Tool (SWAT) model, we evaluated streamflow simulations in the YRB based on different reservoir operation and irrigation schemes. Simulated streamflow with the reservoir operation scheme optimized by the RiverWare model better reproduced measured streamflow than the simulation using the default SWAT reservoir operation scheme. Scenarios with irrigation practices demonstrated higher water losses through evapotranspiration (ET) and matched benchmark data better than the scenario that only considered reservoir operations. Results of this study highlight the importance of reliably representing reservoir operations and irrigation management for credible modeling of watershed hydrology. The methods and findings presented here hold promise to enhance water resources assessment that can be applied to other intensively managed watersheds.
Ever-intensifying human activities have profoundly affected terrestrial water cycling across the globe (Jackson et al., 2001), particularly at the watershed scale (Vörösmarty and Sahagian, 2000; Yang et al., 2014, 2015). Water management substantially alters natural regimes of streamflow through modifying retention time and water exchanges among different components of the terrestrial water cycle (Haddeland et al., 2007). Hydrologic consequences of management activities should be explicitly investigated for effective water resource management (Siebert et al., 2010), especially for watersheds striving to maintain sustainable water supply for multiple users. Accurate simulation of water cycling in intensively managed watersheds faces challenges in reliably characterizing influences of management practices (e.g., reservoir operations and cropland irrigation) on the hydrologic cycling (Wada et al., 2017). Explicit analyses of how model representations of water impoundments and withdrawals would affect hydrologic modeling are needed to advance knowledge of water cycling in managed watersheds.
Construction of dams and reservoirs has substantial influences on the magnitude and variability of downstream runoff (Lu and Siew, 2006; Vicente-Serrano et al., 2017). For example, reservoir operations reduced 9 %–25 % of summer runoff to the Pacific Ocean in the western US and Mexico (Haddeland et al., 2007). In heavily dammed regions, reductions of streamflow following dam construction even reached 100 % (Graf, 1999). Reservoir operations affect the temporal variability of streamflow at multiple temporal scales in different regions across the globe (Huang et al., 2015; Zajac et al., 2017). Regulated streamflow from reservoirs to downstream areas contributes to attenuating flood peaks and volumes, but it could increase baseflow in dry seasons (Batalla et al., 2004).
Reliable representation of reservoir operations in hydrological models is critical for credible simulation of water cycling (Coerver et al., 2018). To characterize impacts of reservoir operations on watershed hydrology, multiple methods have been developed to simulate reservoir releases. These models include mathematical tools which optimize water release for achieving management objectives, simulation models which consider physical processes of water cycling in reservoirs to allow users to evaluate impacts of different management alternatives on reservoir storage and releases, and a combination of these two types of models for reservoir planning and management (Branets et al., 2009; Dogrul et al., 2016; Yeh, 1985). Among these models, the RiverWare model and models developed based on RiverWare consider both management policies and physical processes (Zagona et al., 2001) and have proven their capability of simulating reservoir storage and downstream flows. However, how reservoir operations affect watershed hydrology is still not explicitly examined.
In addition to reservoir operations, cropland irrigation also affects watershed hydrology. Water withdrawal for irrigation has been widely adopted to increase crop production in arid and semi-arid regions. Water redistribution through irrigation enhances water and energy fluxes between soils and the atmosphere (Rost et al., 2008; Sacks et al., 2009) and results in elevated water loss through evapotranspiration (Hao et al., 2015; Malek et al., 2017; Polo et al., 2016) and depletion of water resources (Aeschbach-Hertig and Gleeson, 2012) in different regions of the world. To better simulate impacts of irrigation, numerical models have been developed to quantify water fluxes among soils, vegetation, and water bodies induced by irrigation (Leng et al., 2013; Santhi et al., 2005). Impacts of irrigation on watershed hydrology should be further evaluated to apply this tool for effective management of water resources in basins with competing demands for water.
The Soil and Water Assessment Tool (SWAT) has been widely used to simulate water cycle dynamics in response to management practices across the watershed and regional scales (Arnold et al., 1998). Previous studies indicated that the default SWAT reservoir operation scheme which simulates water release based on target storage may either overestimate reservoir storage in no-flood seasons (Lv et al., 2016) or underestimate water releases when actual reservoir storage is lower than the target storage (Wu and Chen, 2012). SWAT simulates water withdrawal for irrigation from different water sources (e.g., reservoirs, streams, and groundwater aquifers). Multiple efforts have employed SWAT to evaluate impacts of different irrigation practices on watershed hydrology (Ahmadzadeh et al., 2016; Chen et al., 2017; Maier and Dietrich, 2016) and emphasized the importance of balancing water supply and irrigation demands in hydrologic simulations. However, applicability of SWAT in watersheds with interacting reservoir operations and irrigation has not been well studied and thus deserves further investigation to inform effective water resource management.
The Yakima River basin (YRB) in the Pacific Northwest of the US has been regulated for regional hydropower, flood control, fishery, crop cultivation, and drinking water supply. Water supply for irrigation is one of the most important water resource management objectives in the YRB (USBR, 2018). The Yakima River reservoir system supplies water to 180 000 ha of cropland through the operation of five reservoirs which store ca. 30 % of the mean annual runoff of the basin (Vano et al., 2010). Reservoir operations and cropland irrigation in the YRB altered historical streamflow regimes, resulted in severe low flow, and elevated flow events. Since the 1990s, increasing demands for irrigation, municipal water consumption, and critical environmental flow for conserving wildlife habitats in the context of climate change have challenged water resource management in the basin. Thus, there is an urgent need to reliably simulate water cycling in the basin to provide a solid basis for policy formulation and management actions which strive to achieve a balance among water demands for different purposes (Poff et al., 2003).
In recognition of the challenges in modeling hydrology in heavily managed watersheds, this study investigated impacts of water management on streamflow modeling in the YRB. Using the YRB as a test bed, we evaluated streamflow simulations with different model representations of management activities. Objectives of this study are to (1) examine how different representations of reservoir operations influence watershed streamflow simulations and (2) assess impacts of cropland irrigation on watershed hydrology. Methods and findings derived from this study hold promise to provide valuable information for improving hydrologic modeling in intensively managed basins across the globe.
The Yakima River basin (Fig. 1) is located in central Washington, US
(45.98–47.60
Location and land use of the Yakima River basin (67, 99, 160, and 171 are sub-basins used for streamflow calibration and validation). BARL: spring barley; CORN: corn; FRSD: deciduous forest; FRSE: evergreen forest; FRST: mixed forest; HAY: hay; ORCD: orchard; PAST: pasture; POTA: potato; RNGB: range bush; RNGE: range grasses; SWHT: spring wheat; URHD: residential – high density; URLD: residential – low density; URMD: residential – medium density; WATER: water; WETF: wetland – forested; WETN: wetland – non-forested; WWHT: winter wheat).
Settings of the five reservoirs, including locations, height, storage capacity, operating purpose, and surface area were compiled and added to SWAT input files (Table 1). We use three scenarios (R0, R1, and R2) to evaluate reservoir operation simulations in the YRB. Scenario R0 does not simulate reservoir operations, and we use it as a baseline scenario. Scenario R1 uses the SWAT model's built-in reservoir management schemes, which specify monthly target volumes for managed reservoirs (Neitsch et al., 2011). Under the R2 scenario, the SWAT model uses reservoir releases calculated by the RiverWare model as the outflow from these reservoirs to downstream reaches.
The SWAT model calculates water balance for a reservoir on a daily scale as
follows:
Reservoir information of the YRB's five reservoirs (locations are marked in Fig. 1).
Under the R1 scenario, the target release approach calculates reservoir
storage using the following equations:
When the target volume is determined, the reservoir outflow
(
Under the R2 scenario, outflow from a reservoir is calculated based on the
estimated daily release provided by the RiverWare model as follows:
RiverWare simulates operations and scheduling of reservoir management objectives, including hydropower production, flood control, and irrigation (Zagona et al., 2001). RiverWare can model a variety of physical processes for reservoirs with computational time steps ranging from 1 h to 1 year. In RiverWare simulations, the solver is based on operating rules or operating policies that provide instructions for operation decisions such as reservoir releases (Zagona et al., 2001). The rules are strictly prioritized, with high-priority rules requiring that reservoir release should not be less than the minimum flow for downstream reaches, whereas a low-priority rule requires that reservoir storage should fit a seasonal guide-curve value. Conflicts are resolved by giving higher priority rules precedence. This model has been applied to the YRB to simulate outflow from the reservoirs (USBR, 2012).
SWAT irrigation schemes consider multiple water sources including
reservoirs, streams, shallow aquifers, and sources outside the watershed.
Irrigation can be triggered by a water stress threshold (a fraction of
potential plant growth). In SWAT, water stress is simulated as a function of
actual and potential plant transpiration:
We used a plethora of geospatial datasets to parameterize and drive
hydrological simulations in the YRB (Table 2). Topography information was
derived from U.S. Geological Survey (USGS) National Elevation Dataset (NED;
Dataset used in the SWAT simulations. n/a – not applicable.
We quantified parameter sensitivities with a global sensitivity method
described by Abbaspour et al. (2017), which
employs model runs driven by randomly sampled parameter sets, a
multi-regression approach, and a
Parameter sensitivity analysis under various scenarios.
Continued.
Table 3 shows the ranking of parameter sensitivity under different scenarios. In general, selected parameters demonstrated similar sensitives among all scenarios, particularly for the 10 most sensitive parameters, indicating snow melting (SMFMX, SFTMP, and SMTMP), soil water dynamics (CN2, SOL_k, and SOL_Z), and water routing (CH_N2 and SLSUBBSN), which are critical for water cycling in the basin (Tables 3 and S1 in the Supplement). For all scenarios, the most sensitive parameters are CN2 and the snow factors, including SFTMP, SMTMP, SMFMX, SMFMN, and TIMP, indicating that snowmelt is the key hydrological process in the YRB. SWAT uses the Soil Conservation Service Curve Number method (SCS-CN) to predict runoff. As a result, parameter CN2 affects the partition of water between the surface runoff and infiltration and has significant impacts on streamflow estimates. We also observed that sensitivities of several parameters were different among the five scenarios. Specifically, parameters relevant to reservoir operations or irrigation management, including the RES_K and NDTARGR, played important roles in simulations with reservoir operations. The differences could be attributed to the inclusion of reservoir operation and irrigation schemes and further suggest that significant impacts of the management activities on water cycling should be considered in hydrologic modeling. Note that although the inclusion of management activities altered the sensitivity of reservoir- and irrigation-related parameters, snow melting and soil water dynamics may still play the fundamental role in water cycling, as evidenced by the high sensitivity of CN2 and SFTMP.
Without considering impacts of reservoir operations and water withdrawals on
water cycling, the R0 scenario demonstrated poor performance in streamflow
simulations (Fig. 2, Table S2). Streamflow simulations in R1 and R2 were
significantly improved when reservoir operation schemes were added to SWAT,
which further confirmed the importance of considering reservoir operations
in hydrologic modeling in the YRB. Note that reservoirs either increase or
reduce streamflow, as reservoirs could increase water release in dry
seasons or retain upstream water for flood control in wet seasons. In
addition, streamflow simulated in the R2 scenario (average correlation
coefficient of 0.59) showed a better agreement with measured flow than that
of the R1 scenario (average correlation coefficient of 0.57). R2 exhibits
better Ens in three of the four sub-basins than R1 (Table S2), indicating that
reservoir outflow estimated by RiverWare more accurately simulated water
releases than the default reservoir operation scheme in SWAT. The streamflow
simulations in sub-basins 67 and 99 were more sensitive to the different
reservoir schemes, as evidenced by greater improvements in the Ens and
Calibration and validation results in four sub-basins under the R0 scenario (baseline simulation does not consider management activities).
We also compared ET simulations of the YRB under the three scenarios (R0, R1, and R2). Specifically, ET estimates increased in May and June but decreased in winter for R1 and R2 simulations (Fig. S1 in the Supplement). In addition, annual ET increased by 7.83 % and 8.05 % for R1 and R2 simulations relative to the R0 simulation, respectively (Fig. 5). The changes could be attributed to increased evaporation from reservoirs.
Settings of scenario R2S1, which used reservoirs and streams as water sources for irrigation, are consistent with the actual irrigation practices in the YRB, where surface water is the primary irrigation water source (Fig. 6). For the R2S2 scenario, shallow groundwater was assumed to be the water source for irrigation (Fig. 7). Consequently, streamflow simulations under the scenario R2S1 matched observations better than those under R2S2. Compared with the R2 scenario, the simulated flow decreased by 24.87 % and 31.29 % in R2S1 and R2S2, respectively.
Calibration and validation results under the R1 scenario (default SWAT schemes for reservoir operations).
Calibration and validation results under the R2 scenario (RiverWare for reservoir operations).
ET is an important component of terrestrial water cycling, and this variable is used in the calculation of irrigation demand in SWAT simulations. Figure 8 compares simulated monthly ET of the irrigation scenarios (R2S1 and R2S2) with the RiverWare reservoir operation scenario (R2) which did not consider irrigation. The mean monthly ET rates of the irrigation scenarios (R2S1 and R2S2) were significantly higher (85 % and 63 % for R2S1 and R2S2, respectively) than simulations without irrigation, particularly during March–July, when irrigation was applied to support crop growth.
Annual ET simulated under reservoir operation-only scenarios (R0, R1, and R2).
We further compared the simulated annual ET in the R2S1 and R2 scenarios (Fig. 9). We observed low cropland ET in the R2 scenario relative to the R2S scenario. Specifically, when irrigation was included in our simulation, SWAT ET estimates increased by ca. 85 % at the annual scale. Monthly scale comparison showed that increases in ET mainly occurred in growing seasons (April to August, Fig. S2). The comparison demonstrated that inclusion of irrigation schemes achieved better estimates of water losses during irrigation and contributed to enhancing streamflow simulations (Fig. 6). In addition to magnitude, the irrigation scenario (R2S1) also simulated the interannual variability of ET well, as evidenced by the high coefficient of determination in the scatter plot against ET estimates based on remote sensing data (Fig. S3).
The mean annual irrigation depth for the irrigation scenarios of R2S1 and
R2S2 was 480.66 and 228.46 mm yr
Calibration and validation results under the R2S1 scenario (RiverWare for reservoir operation and surface water as the water source for irrigation).
Calibration and validation results under the R2S2 scenario (RiverWare for reservoir operation and groundwater as the water source for irrigation).
As indicated by the improved Ens and
ET in the composite scenarios (R2S1 and R2S2) was higher than the R0 scenario, which can be attributed to the elevated evaporation from reservoirs and irrigated cropland. Direct evaporation from reservoirs increased by 7 %–8 % over the study period (1980 to 2010), due to the improved simulation of reservoir surface areas in the R1 and R2 simulations relative to the R0 simulation. Irrigation practices led to more pronounced increases in ET in R2S1 and R2S2 simulations, as compared with those of R2 (Fig. 8). These results indicate that irrigation may have more pronounced impacts on ET through stimulating ET than reservoir operations in the study area.
Monthly ET simulated under the irrigation operation scenarios (R2S1 and R2S2) relative to the reservoir operation-only scenario (R2).
Comparison of ET simulations for cropland during 2000–2009 under the R2 and R2S1 scenarios.
In recent decades, water users of the YRB passed the Yakima River Basin Integrated Water Management Plan, which is a comprehensive agreement that advances water infrastructures and management (USBR, 2012). Enhanced hydrologic modeling provided by this study will provide valuable information for goals of the Integrated Plan, which requires accurate streamflow information for managing water resources to meet ecological objectives as well as for securing water supply for domestic uses.
Although previous investigations highlighted the importance of irrigation and reservoir management to water balance and availability (Hillman et al., 2012; Malek et al., 2014), joint impacts of these two water management practices on watershed hydrology have not been fully understood. In recognition of this challenge, we enhanced SWAT representations of the two critical water management activities, including reservoir operations and irrigation, to constrain uncertainties in hydrologic simulations. We achieved improved model performances through including the two activities in the SWAT modeling framework. The simulated streamflow was generally lower in simulations with management activities than the baseline simulation (R0). Without including reservoir management and irrigation, SWAT may overestimate streamflow due to the unreasonably estimated water loss through ET.
Water management activities have altered natural hydrological cycling and posed challenges to reliable simulation of watershed hydrology. The YRB is a typical watershed that is regulated to support agricultural production. Maintaining sustainable water supply in basins like the YRB calls for a sound understanding of hydrological impacts of management activities. Management schemes developed and evaluated in this study will be transferable and applicable to future SWAT and other watershed model applications for investigating the water cycling that is influenced by reservoir operations and water withdrawal for irrigation across broader spatial scales.
Reservoir operations have both direct and indirect impacts on streamflow. Water release from reservoirs directly affects the magnitude and variability of streamflow in downstream reaches. Dam and water diversion operations determine the amount and timing of water discharge to downstream river channels. As a result, reservoir operations may either attenuate flood peaks in wet seasons or increase streamflow in dry years, in compliance with minimum instream flow policies (Yoder et al., 2017). In addition, multiple hydrological processes, such as vertical flow in surface or subsurface waters, water routing, evaporation, precipitation, and microclimate, are also responsive to reservoir operations (Lv et al., 2016). Our simulations suggested that reservoir operations altered both streamflow and ET in the YRB.
Most precipitation in the YRB occurs in winter as snowfall. Snowpack serves as a water reservoir for spring and summer streamflow. Consequently, streamflow is high in spring but low in summer. As shown in Table 1, most of the reservoirs were built to support cropland irrigation. The presence of reservoirs positively contributed to water availability in dry periods. Water storage management in reservoirs is one adaptation strategy particularly applicable to snowmelt-dominant watersheds like the YRB, which experiences water scarcity during the summer irrigation season (Yoder et al., 2017), and thus alters natural flow regimes. Without representing reservoir regulations, SWAT simulations failed to reasonably reconstruct temporal variability in streamflow (R0 scenario). Results of this study indicated that reservoir algorithms based on RiverWare (R2) were relatively more realistic compared with the default reservoir operation algorithms in SWAT (R1), as evidenced by the improved model performances. Enhanced model performances in the R1 and R2 scenarios further corroborated the significant impacts of reservoir operations on seasonal patterns of streamflow (Adam et al., 2007).
Compared with the baseline scenario (R0), R1 and R2 simulations showed that the ET rates increased considerably from April to September due to reservoir operation. Direct evaporation from reservoirs increased under the R1 and R2 scenarios because of improved estimates of reservoir surface areas. The consideration of such an impact on ET in the R1 and R2 scenarios also contributed to enhanced model performances relative to the baseline scenario (R0).
Water withdrawal for irrigation has increased pressures on maintaining sustainable water resources in the YRB (Malek et al., 2017). Insufficient water supply for agricultural production, drinking water supply, and environmental flows have raised concerns on the local economy and ecosystem integrity (Hillman et al., 2012). Due to the significant impacts on soil moisture and plant growth, the amount and timing of irrigation influence ET losses and watershed hydrology (Maier and Dietrich, 2016). As a result, the irrigation impacts on streamflow should be evaluated to provide reliable estimates of streamflow in basins like the YRB, which helps balance the water supplies and demands for effective water resource management.
As reported in previous studies, most of the water for agricultural irrigation was provided by surface water, and one-third was from groundwater in the YRB (USBR, 2012). Under the R2S1 scenario, our assumption that irrigation water was from the reservoirs and streams generally agreed with the actual water uses for irrigation in the basin. The less satisfactory model performances in the R2S2 scenario may stem from the unrealistic assumption of water sources, irrigation efficiencies, and the return flow of irrigation. In addition, SWAT simulates streamflow based on water balance among multiple water pools, including shallow groundwater which is recharged by subsurface runoff (Shadkam et al., 2016). Under the R2S2 scenario, water withdrawal from the shallow renewable groundwater was used in our simulation. This simplification did not consider water withdrawal from deep nonrenewable aquifers. As a result, water availability based on shallow groundwater for irrigation and the groundwater recharge may have been unreasonably estimated and partially contributed to unsatisfactory model performances under this scenario (R2S2).
To better investigate hydrological consequences of water management, future studies should further constrain uncertainties in streamflow simulations by incorporating additional reservoir management and irrigation information. Including of observed reservoir release will help improve model representations water discharge from reservoirs. In addition, model representation of irrigation should be improved in the future. Note that model performances of the R2S1 scenario were not substantially improved relative to the R2 scenario. The irrigation operation scheme that used surface water as the single source may have introduced uncertainties to streamflow simulations, since groundwater is also an important water source for irrigation, particularly in dry years in the YRB. Future simulations need to incorporate explicit irrigation information about irrigated areas and the source, amount, and timing of groundwater withdrawals into hydrologic modeling to better simulate agricultural hydrology. We observed different seasonal patterns of ET under the five scenarios. How management activities affected water and energy exchanges between soil and the atmosphere should also be investigated in the future.
As most reservoirs were built for irrigation in the YRB, impacts of reservoirs should be assessed jointly with the accelerating development of irrigated agriculture in the basin. The presence of reservoirs positively contributed to water availability for irrigation, particularly for dry seasons. In general, the combination of reservoir operations and irrigation have reduced streamflow in the YRB when compared with the baseline scenario (R0). This is attributable to the large amounts of water loss through ET in irrigation and additional water storage in reservoirs.
Among the multiple modeling scenarios, we found that linking RiverWare reservoir model with SWAT achieved better performance than those model structures that reply on simplified reservoir operations, as evidenced by the relatively higher correlation coefficient and the Ens. However, it is worth noting that these statistical metrics are calculated based on a limited set of hydrological variables (e.g., streamflow) but cannot guarantee that other hydrological processes are well represented (Zhang et al., 2013). Therefore, we further used MODIS-estimated ET and reported irrigation water demand data to justify the favorable performance of the combined SWAT-RiverWare watershed model configuration.
Our model evaluation process follows the widely accepted procedures for model calibration and evaluation (Moriasi et al., 2007; Arnold et al., 2012). We also would like to point out that the complexity difference between the SWAT-RiverWare and other watershed model configurations was not explicitly considered in model evaluation. Previous research notes that model complexity is an important factor in selecting the most robust model configuration that can fulfill a specific purpose. For example, Höge et al. (2018) reviewed existing methods and laid the foundation for a comprehensive framework for understanding the critical role of model complexity in model selection. The lack of reliable prior knowledge of the model structure and associated model parameters makes it difficult to directly consider model complexity here. However, the framework laid out by Höge et al. (2018) deserves further exploration in comparing the performance of different watershed model configurations in the future.
Reservoir operations and irrigation have substantial impacts on water
cycling globally. Hydrologic simulation in the managed basins faces
challenges in reliably characterizing water management activities. This
study assessed the hydrological impacts of reservoir systems and irrigation
practices through numerical model experiments with SWAT. Representing
reservoir operations by coupling the RiverWare model and SWAT significantly
improved streamflow simulations. We achieved reasonable model performances
in the scenario using reservoirs and streams as the water sources for
irrigation, since these assumptions are consistent with the actual
irrigation practices in the basin. Model simulations suggested that
reservoir operations and irrigation water withdrawal generally reduced
streamflow by enhancing water loss through ET in the study area. Results of
this study demonstrated the importance of incorporating water management
activities into hydrologic modeling. Both SWAT and RiverWare are community
models that have been widely tested and applied in diverse regions across
the globe, as evidenced by the numerous peer-reviewed publications in the
fields of reservoir operation and watershed modeling
(
Sources of the geospatial and climate forcing data used to configure the SWAT model have been described in the “2 Materials and methods section”. The simulations of streamflow and evapotranspiration data are shown in the figures of the paper and are available through contacting the authors.
The supplement related to this article is available online at:
JQ and QY contributed equally to the paper. JQ, QY, XZ, and MH designed the study and conducted the watershed configuration and modeling efforts. JCA and KM contributed RiverWare reservoir outflow data. All authors contributed to result analysis and paper preparation.
The authors declare that they have no conflict of interest.
We greatly appreciate the valuable comments from four anonymous reviewers, which have much improved the quality of the paper. This research was supported by the U.S. Department of Energy, Office of Science, as part of research in Multi-Sector Dynamics, Earth and Environmental Systems Modeling program. Jennifer C. Adam's and Keyvan Malek's contributions were supported by the NSF EAR no. 1639458 and USDA no. 2017-67004-26131 awards. Edited by: Alberto Guadagnini Reviewed by: four anonymous referees