Approach
The authors investigate optimal reservoir operations across hydrometeorological conditions to identify a system performance threshold under climate change uncertainty. That is, system operation is formulated as a feedback control policy, and the overall objective is to find the theoretical upper limit for a water resources system's adaptation to climate change uncertainty, i.e. the maximum operational adaptive capacity beyond which even optimal operational rules fail to achieve societal objectives for a system of interest.
The approach is a variant of the formal method referred to as "decision scaling". Testing performance of the optimal feedback control policies across changes in temperature and precipitation is a form of stress testing, and the identification of failures to meet objectives is a form of vulnerability analysis.
Feedback control operating rules are modelled using an approximation technique called "Gaussian radial basis functions" (RBFs), describing reservoir releases over time and as a function of lake level. In addition to estimating a historical operating rule, Pareto-optimal RBF policies are identified using evolutionary multiobjective "direct policy search".
Methodological steps
- Generation of an exposure space (H):
- making incremental changes to hydrological variables' attributes. A n-dimensional exposure space is created, made of selected attributes (h), such as H = [h1,...,hn].
- sampling the exposure space using a scaling of historical data to develop perturbed time series for the relevant climatic variables.
- Assessment of current system performance for each point of the exposure space:
- Running a simulation model with current operation policy with each climate input (see exploratory modelling)
- Separating "failure" from "success" conditions based on the definition of thresholds for acceptable values.
- Assessment of the optimised system performance for each point of the exposure space, i.e. with optimal operation rules for each climate input
- Superimposition of climate projections onto the exposure space for the current and optimised system performance to identify climate conditions and timing for which the system may fail and when to implement management strategies to maximise a system adaptation capacity.
Case study details
Lake Como is a regulated lake in Northern Italy with an active storage capacity of 254 Mm3, sustaining water needs for irrigation and water supply, together with ensuring flood protection and environmental flow requirements in the Adda River. Its subalpine location implies that the lake is dependent on snowmelt to be filled, which requires the definition of strategic operation rules on an annual basis to support flood control and irrigation water demand objectives. In addition, changes in climate conditions and snowmelt runoff into the lake could increase tensions in the region.
The problem is identification of operational rules to maximise the Lake Como system's adaptive capacity under climate change uncertainty
Conditions of failure for the system :
- daily flooded area over 100 m2
- daily average squared water deficit = 400 kL2
The optimal operation rule is based on a penalty function assessing reliability:
- flood reliability function: the ratio of the days without a flood event within the simulation horizon
- irrigation reliability function: ratio of available supply to the demand averaged across all time steps, with ratio values capped at 1
Climate change projections in the climate exposure space are obtained from multiple models for the years 2025, 2050, 2075 and 2100, with statistical downscaling based on quantile mapping to generate future local climate change projections.
Attributes for generating the climate exposure space:
- Mean annual precipitation (mm): Precipitation was perturbed as a percentage change to daily wet days (over 1mm) from 90% to 130% of current values.
- Mean annual temperature: Temperature was perturbed for each daily time-step from -5℃ to 15℃ of current values.
The main result is the identification of the upper limit for adaptive capacity (failure boundary) to support irrigation and flood protection objectives in Lake Como. The identification of such a limit allows for dynamic management and the setting of optimal operating policies acoording to climate changes. It was found that the Lake Como system could adapt, by modifying operating policies, to more than three times as many hydrometeorological states compared to current operating rules.
The result is conditional on certain key assumptions:
- the selection of sensitive variables to generate the exposure space for which operational rules are sensitive. For that matter, a sensitivity analysis of the inputs is recommended to identify such variables.
- accurate modelling of stakeholder demand to ensure the proper identification of conditions of failure
- the proper setting of operating policy structure and constraints for the optimisation process
- the choice of an appropriate model to represent reservoir decisions, able to change releases in response to extreme climate conditions
Future considerations to improve the ability of the framework to identify upper limits include:
- considering changes in water users' demands over time and application of a nonstationary system model within the framework
- including more than two dimensions to generate the exposure space with an increased number of uncertain variables that influence the effectiveness of the operation rules
Reference
Also see presentation for the Queensland Water Modelling Network.