1 Uncertainties: New Challenges for Decision Making on Fresh Water Supply Systems

Decision making and planning for fresh water supply are challenging because of a number of different reasons: First, a variety of individual measures or tactics need to be considered at both the supply and the demand side, including regulation, adaptation of infrastructure elements, subsidies or fines, operational control and coordination, etc., and appropriate combinations of these need to be made in the process of policy design. Second, different scales need to be considered, such as the local, regional, national and river basin scales, each having their associated governmental bodies. Third, a range of sometimes conflicting interests needs to be considered. Fourth, various uncertainties come into play. Traditionally, ample attention has been given to natural climate variability, for which empirical data are available and probabilistic approaches can be used to estimate the likelihood and consequences of various climate scenario’s, and to assess the effectiveness of alternative policies. It has been argued that under the rapidly changing global climate, the assumption of stationarity is no longer valid, meaning that historic return frequencies of extreme hydrological and weather conditions can no longer be relied on in water management: here the past is no longer the key to the future (Milly et al. 2008). Long-term decision making, involving, for example, decisions regarding water supply infrastructures having long lifetimes, brings a new class of important uncertainties. Climate change, itself being uncertain, may have significant – yet uncertain – impacts on precipitation and evaporation patterns, salt intrusion, salinization, and seasonal variation. In addition, a variety of uncertain but possible developments in land use and agricultural conditions, food markets, spatial developments, economic developments, future technology, and future societal values may also affect the future impacts and valuation of policies that have to be decided in the short term.

Partly but not exclusively triggered by the need to explore the impacts of climate change, during the last decade, the literature on dealing with various sorts of uncertainties has been developing rapidly (e.g., Maxim and Van der Sluijs 2011; Van Bree and Van der Sluijs 2014). Therefore, in this paper, we will first give a brief overview of some of the major developments in the uncertainty literature. Next, three examples are presented of application of one of the many approaches under development in the field of fresh water supply. We conclude with a few general observations on the state of the art, and recommendations for further work in the field.

2 Developments in Dealing with Uncertainties

Knowledge available for long-term policy making is unavoidably incomplete and often controversial. In the science-policy interface, strategic use of uncertainties seems to dominate: uncertainties are either downplayed to promote political decisions or they are overemphasized to prevent political action. Both are problematic and lead to policy strategies that are not fit for meeting the challenges posed by the uncertainties and complexities faced. While the scientific community has not ignored the problems of uncertainty, most efforts thus far have been associated with analysing and dealing with those uncertainties that can be described quantitatively (e.g., probabilistically, like in risk analysis), and with efforts at reducing uncertainties through additional research and/or refined modelling.

Not all uncertainties encountered in decision making can be quantified. In fact, unquantifiable uncertainties can sometimes be more relevant than the part where we do have enough knowledge to quantify uncertainty. Rather than using such unquantifiable uncertainties in a strategic way, or ignoring them, there is a growing consensus among scholars that deliberate policy strategies in the face of often irreducible uncertainties should be developed and implemented (Van der Sluijs et al 2008), and a wide variety of approaches and methods have emerged in the literature on the subject. As a result, the field is quite fragmented, and different key authors and groups have developed their own terminologies. A certain degree of conceptual confusion results, as approaches partially overlap, and different authors use the same terms in different ways. For example, approaches like ‘Adaptive Management’, ‘Adaptation Policy Framework’, and ‘Robust Decision making’ strongly (but not completely) overlap in their guiding principles.

Despite the fragmented character of the field, there seems to be a general agreement that different approaches are suitable under different degrees and types of uncertainty, that in real world situations often different types of uncertainty are present at the same time, and that combinations of approaches need to be fine-tuned to the characteristics of the situation. We will now briefly elaborate on these aspects.

2.1 Types and Levels of Uncertainty

In an attempt to harmonize diverse literatures in as far as they relate to model-based decision support, Walker et al. (2003) presented an uncertainty framework. The starting point for this framework is the distinction between an analyst’s perspective on uncertainty and a decision maker’s perspective on uncertainty (van Asselt 2000). The framework focuses on the analyst’s perspective on uncertainty. A central idea of this framework is that uncertainty is a multi-dimensional concept. Ironically, this framework in turn had its own variants, motivated by various perceived shortcomings of the original typology or perceived need to tailor it to a specific field of study. More recently, Kwakkel et al. (2011) presented a review of the variants and what motivated them, and presented a new framework that integrated and harmonized the variants. Both the original and the new framework categorize uncertainties by their location, nature, and level. The location dimension focuses on where the uncertainty is located. The nature dimension focuses on the character of the uncertainty, i.e., whether the uncertainty is knowledge-based, a direct consequence of inherent variability, or related to human interpretation (ambiguity) or human values now and in the future. The level dimension focuses on the severity of the uncertainty.

Each of the three dimensions is relevant when selecting an appropriate approach for handling the uncertainty. However, the level dimension plays the most important role. Broadly speaking, the level of uncertainty is the assignment of likelihood to things or events. In some cases the likelihood or plausibility of these things, events or situations can be expressed using numbers, but in other cases more imprecise labels are used, such as more likely, less likely, or equally likely. Overall, the level of uncertainty ranges from complete certainty to absolute ignorance. We define five intermediate levels of uncertainty: we speak of level 1 uncertainty, or recognized uncertainty, when one admits that one is not absolutely certain but when one is not willing to measure the degree of uncertainty in any explicit way (Hillier and Lieberman 2001). We speak of Level 2 uncertainty, or shallow uncertainty, when one is able to enumerate multiple possibilities and is able to provide probabilities. We speak of Level 3 uncertainty, or medium uncertainty, when one is able to enumerate multiple possibilities and able to rank order them in terms of perceived likelihood. However, how much more likely or unlikely one possibility is compared to another cannot be specified. We speak of Level 4 uncertainty, or deep uncertainty, when one is able to enumerate multiple possibilities without being able or willing to rank order the possibilities in terms of how likely or plausible they are judged to be. Finally, we speak of Level 5 uncertainty, or recognized ignorance when one is unable to enumerate multiple possibilities, while admitting the possibility of being surprised (Kwakkel et al. 2011). Table 1 summarizes the five levels.

Table 1 The five levels of uncertainty (adapted from Kwakkel et al. 2010b)

2.2 A New Paradigm for Handling Uncertainty

Uncertainties pose a significant challenge to planning and decision making. The dominant approach in many fields has been to ignore the uncertainties, to quantify them into error margins, to try and reduce them, or to deal with only those uncertainties that can be easily quantified (Quade 1982; Dempsey et al. 1997; Marchau et al. 2009; Van Geenhuizen et al. 2007; van Geenhuizen and Thissen 2007; McDaniel and Driebe 2005). However, such approaches suffer from the problem that they focus on those uncertainties that are “among the least of our worries; their effects are swamped by uncertainties about the state of the world and human factors for which we know absolutely nothing about probability distributions and little more about the possible outcomes” (Quade 1982). Similarly, Goodwin and Wright (2010) (p. 355) demonstrate that “all the extant forecasting methods – including the use of expert judgment, statistical forecasting, Delphi and prediction markets – contain fundamental weaknesses.” And Popper et al. (2009) state that the traditional methods “all founder on the same shoals: an inability to grapple with the long-term’s multiplicity of plausible futures.” In response to this, various new planning approaches have been put forward (e.g., Lempert et al. 2003; Walker et al. 2001; de Neufville 2000, 2003; Dewar 2002; Dewar et al. 1993; Holling 1973; Lempert 2002). These approaches generally contain three types of elements:

  • analytic methods and tools to capture and assess the (consequences of) uncertainties,

  • specific action or policy strategies in light of irreducible uncertainties,

  • specific, often participatory processes to involve stakeholders and decisionmakers

Analytic approaches (Swanson et al. 2010) generally emphasize the need for a more thorough integrated forward-looking analysis of the uncertainties through techniques such as exploratory modelling and analysis (Agusdinata 2008; Lempert et al. 2003), bounce casting (Kahan et al. 2004), and scenarios in various forms (Bradfield et al. 2005; Varum and Melo 2010).

The action or policy strategies can be divided in three partly overlapping categories: First, so-called top-down or robust strategies build on analytic approaches in order to assess the vulnerability of alternative strategies, and to select those policies that do well under a wide variety of possible futures. Second, so-called bottom-up, resilience-based strategies recognize the limited capability of analytic approaches for anticipating rare events, and focus on the system-to-be-governed or to-be-designed itself, and try to enhance the system’s inherent ability to cope with disturbances and its self-organizational capacity in a comprehensive way (e.g., Wardekker et al. 2010). Third, there is a growing interest in flexibility and adaptability in plans in which a strategic vision of the future is combined with short-term actions and a framework that can guide future actions (Albrechts 2004; Walker et al. 2001; Walker et al. 2013).

Together with a focus on involving stakeholders and decisionmakers in both the identification of uncertainties and the assessment of alternative strategies, these new approaches can be considered to form a new emerging paradigm for handling uncertainty differently. They are based on accepting the uncertainties, and focusing on the question what can best be done now, given that we don’t know what the future will bring.

For a more complete overview of concepts, approaches and methods, the reader is referred to recent review reports and articles (Dessai and Van der Sluijs 2007, 2011; Wardekker 2011; Walker et al 2013; Lourenço et al 2014a; van Bree and van der Sluijs 2014).

2.3 Applications in the Field of Fresh Water Supply

As explained above, fresh water supply planning faces a variety of uncertainties, part of which cannot adequately be dealt with in a (purely) probabilistic way. Recently, a number of example applications of the type of approaches mentioned above have been reported in the literature. For example, Groves, et al. (2007; 2012) have reported on a number of studies on fresh water supply in California and Colorado River, using the so-called ‘robust decisionmaking’ approach that is strongly based on analytical and computational techniques. Matrosov et al (2013a; 2013b) apply among others Robust Decision Making to a fresh water supply case in the United Kingdom. Groot et al. (2014) report, among other things, on several practical cases of water supply planning in Portugal, the UK, and Hungary. These studies mostly were related to dealing with uncertainties of the scenario-type, and used a combination of multiple methods such as expert elicitation, stakeholder involvement and model-based sensitivity analysis. The studies demonstrate, on the one hand, the complexity of the challenge, and, on the other, the relevance of more explicitly dealing with uncertainties, and ‘a shift towards a flexible, robust and no-regret approach’ (Groot et al. 2014, p 67). Lourenço et al. (2014b) also emphasize a preference for ‘options that contribute to enhance resilience and adaptive capacity’.

In another study in Hungary (Malatinszky et al. 2013), the focus was on climate adaptation strategies for wetlands and grasslands. The authors emphasize the importance and various benefits of stakeholder involvement in such processes, which include “enhanced awareness, willingness to taking action, inclusion of local knowledge, information exchange among affected parties, identification of win-win-solutions for land users and nature conservation, and building trust in authorities”.

Clearly, while there seems to be agreement on the general tendencies of preferred approaches, much remains to be developed and learned, both with respect to methods and tools and their potential and limitations in practical applications and with respect to what strategies fit what uncertainty situations.

In the following sections, we will therefore, in order to learn more about their applicability, illustrate the application of three novel but different approaches to two different real-world fresh water supply cases in The Netherlands. First, the resilience approach is applied to an area of peat grassland in the west of the Netherlands. Second, a system robustness approach is applied to a low-lying polder area mainly used for agricultural purposes. Third, an exploratory modeling approach is applied to the same case study area.

3 A Resilience Approach: the Venen-Vechtstreek Region

3.1 Resilience

Resilience is a concept that emerged in ecology research in the 1960s. It was used in relation to the stability of ecosystems in the face of a wide range of external perturbations: “the amount of disturbance that a system can absorb and still remain within the same state or domain of attraction” (Holling 1973; Walker and Salt 2006). The concept has since been adopted by various disciplines, ranging from engineering to psychology, disaster studies, and climate change adaptation. Walker et al (2004) define resilience in relation to ‘social-ecological systems’ as: “the capacity of a system to absorb disturbance and reorganize while undergoing change so as to still retain essentially the same function, structure, identity, and feedbacks”. Resilience focusses on coping with disturbances and preventing collapse into a qualitatively different state, rather than predicting and preventing disturbances.

In relation to climate change adaptation, resilience can be understood as a bottom-up approach, which aims at strengthening the system, removing vulnerabilities, and enhancing the capacity to cope with and recover from climate change-related disturbances and surprises. The case-study examined the consequences of various proposed measures and developments on the climate-resilience of the freshwater systems in the Venen-Vechtstreek region. These measures are plans to manage the area in general and are not specifically proposed to adapt to climate change.

Resilience analysis can take several steps. The first step is to determine ‘resilience of what?’. This includes an overview of the area, the trends and challenges it faces, the way it is used, and the key characteristics and functions that make the system for what it is (i.e., if these are not retained, the system has shifted to a different state). The second step is to address the question of ‘resilience to what?’. In this study, an inventory of climate change impacts was made. The third step is to assess the resilience of the system, measures that will enhance resilience, or the resilience consequences of developments. This case study concerned the latter. We applied a set of ‘resilience principles’ describing different mechanisms by which systems can absorb disturbances and retain identity (Watt and Craig 1986; Barnett 2001; Wardekker et al. 2010):

  • Homeostasis: multiple feedback loops counteract disturbances and stabilize the system.

  • Omnivory: vulnerability is reduced by diversification of resources and means.

  • High flux: a fast rate of movement of resources (e.g., information, money, expertise, goods for emergency relief, fuel) through the system ensures fast mobilization of these resources to enable coping with perturbations through for instance reduced feedback delay. If information is proactively and rapidly exchanged, parties can act on this faster when there is a problem.

  • Flatness: the hierarchical levels relative to the base should not be top-heavy. Overly hierarchical systems are too inflexible and too slow to cope with surprise. The more hierarchical levels exists through which stabilizing feedback signals must operate, the lower the signal to noise ratio and the lower the speed of signal transmission, and thus of action to counteract the perturbation or remedy the damage.

  • Buffering: essential capacities are over-dimensioned such that critical thresholds are less likely to be crossed.

  • Redundancy: overlapping functions; if one fails, others can take over.

3.2 Case Study Area and Trends

The Venen-Vechtstreek region is an area of peat grassland in the west of the Netherlands. It is known for its wide views, cow-filled meadows, and long and narrow stretches of land separated by water (Fig. 1). The appearance of the region has changed little over the centuries. The case study focused on the specific area of Groot Wilnis – Vinkeveen. The region is part of a range of wetlands forming a ‘robust ecological corridor’ between areas protected under Natura 2000, the European ecological network of protected areas. The corridor is intended to provide a connection and shelter that allows animals and plants to migrate between nature areas. The area is situated at approximately 2.5 m below sea level and is artificially drained.

Fig. 1
figure 1

Satellite view of the case-study area (source: Google Maps)

Due to the artificial drainage, peat comes into contact with oxygen and decomposes. Consequently, the area suffers from soil subsidence. The most sensitive parts of the area have seen subsidence of up to 12 mm per year. Agricultural areas suffer more subsidence than nature areas, due to higher levels of drainage, resulting in large local differences in water levels. The decomposing peat also releases nutrients into the surface water and greenhouse gasses into the atmosphere.

The area is valuable for both ecological and agricultural reasons. The area’s management covenant (Stichting Ontwikkeling De Venen 2010) specifies that the area should be preserved as an open landscape in which the dairy sector can continue to develop. Four key functions can be distinguished: fresh water supply, nature, agriculture, and recreation. Stakeholders aim to strengthen these functions, as well as to reduce the amount of soil subsidence. The functions have different implications for the management of the area. Connecting them through integrated measures will be an important step in enhancing resilience.

Climate change will affect Dutch freshwater systems through increasing air and water temperatures, sea level rise, changes in seasonal river discharge, and changes in precipitation (PBL 2013). Average precipitation may either increase or decrease in summer, and increase in winter. For the study area, this may result in an increasing potential precipitation deficit during summer. Dry periods could negatively impact water availability and quality for fresh water supply as well as for nature and agriculture. Wet periods could result in flooding, negatively impacting all four functions. They could also exacerbate the current tension in water level requirements for agriculture (relatively low) and nature (relatively high), and enhance nutrient leaching from fields into surface water, negatively impacting nature. Temperature changes can impact various plants and animals. For 2050, a scenario study projects a moderate to high terrestrial drought risk, moderate risk of aquatic eutrophication, low risk of terrestrial eutrophication, and low to high risk of aquatic salinization (Verhoeven 2012). The high salinization risk occurred in a scenario with warm, dry summers, when brackish water from outside the area would be required to manage the water levels.

3.3 Analysis

A list was drafted including measures that are currently planned for managing the study area, based on the management covenant, other policy documents, and position papers of regional and national stakeholders. This list was supplemented with options that could contribute to the area’s overall resilience. These supplementary options were collected through interviews with stakeholders and experts with experience concerning the study area. Options were categorized according to the key functions they contribute to: nature, agriculture, water quality, and recreation.

The research team performed an initial reflection on the resilience implications of the options, followed by a second set of interviews. These reflections were used in preparation of a workshop with stakeholders and experts on the area and on resilience: to select a subset of options to be discussed in the limited time available, and to seed the discussion. Workshop participants (n = 7; this number is adequate for an expert assessment; 6–12 is advisable, to cover relevant fields of expertise, but allowing for sufficient time for discussion; Knol et al. 2010) scored the shortlist of options on each of the resilience principles. They applied a five-point scale, ranging from strong decrease to strong increase of resilience. Aggregated ‘resilience scores’ were also generated. See Fig. 2.

Fig. 2
figure 2

Mean aggregated resilience scores. Options are arranged in four groups, from left to right: agriculture (planned and new) and nature (planned and new)

Nearly all of the assessed options for nature improved resilience, according to the scores of most participants (occasionally, one or two participants disagreed). The participants expected the option ‘marshland construction and capillarity’ to perform particularly well. Agricultural options scored less well; half resulted in (slightly) reduced resilience. The option ‘underwater drainage’ did receive a good score. For both types of options, the range between individual aggregated scores was generally one to two points. Scores on separate resilience scores had similar ranges. They could on occasion diverge more strongly (up to scores of -2 to +2, cf. Fig. 2, although these were exceptions), but the interquartile ranges (25–75 percentiles of scores) were two points in most cases. The group size was sufficiently large that a single expert scoring notably different did not change the interquartile ranges.

The approach to assessing resilience applied in this case study provided a useful tool to scan a variety of options on a range of impacts for a range of stakeholders on their implications. It allows for a relatively quick evaluation of options and the inclusion of stakeholders and their perspectives. The approach was also easy to combine with other decision criteria in a multicriteria assessment, allowing for a broader evaluation of policy options. Preferably, this would include multiple workshops with a broader range of expertise amongst the participants. A downside is that the assessment takes place on an ordinal scale (categories, rather than full quantification) and does not examine spatial patterns. Also differences in opinion between participants scores will remain, even with extensive discussion. These require further investigation and may hint at deeper differences in opinion (uncertainty due to ambiguity; see e.g., Brugnach et al. 2008), or different ways in which an option can be implemented. The approach can be applied for a broad range of cases, and does not require extensive computing power or special software. It does require some insight in the system dynamics of the case study (among researchers and participants), as well as experience with participatory processes (among researchers). It is possible to scan a broader set of options with the research team only, and/or include a larger set of stakeholders by means of a survey. However, discussion with, and among, local stakeholders and experts remains important to analyze local resilience in a meaningful way.

We concluded that, considering that agriculture is a key function of the area, covering most of the land use, it would be advisable to pay attention to increasing resilience in agricultural developments. Scores on separate resilience principles (i.e., not aggregated into a single resilience score) showed that very few options increase flatness. Additionally, few agricultural options address redundancy.

4 A System Robustness Analysis: Case Rijnland

4.1 System Robustness

System robustness is a concept that originates from engineering and biology, where it refers to the ability of systems to maintain system characteristics when subjected to disturbances (Carlson and Doyle 2002; Stelling et al. 2004). For example, a road network is robust to incidents when the effect on traffic congestion is prevented and/or limited (Snelder et al. 2012). According to the Merriam-Webster dictionary, the term robustness indicates that systems or networks are capable of performing without failure under a wide range of conditions. In the context of flood and drought risk management, system robustness is defined as the ability of a socio-economic system to remain functioning under a range of disturbance magnitudes (Mens et al. 2011). The term disturbance refers to an external event resulting from climate variability, for example floods or droughts. The system for which robustness is assessed comprises biophysical as well as socio-economic aspects. For example, the robustness of a flood risk system does not only depend on the flood protection system, but also on the characteristics that determine the flood impact, and on the social and economic capacity to recover from these impacts. Failure of the system thus means that impacts are too large to recover from. Such impacts may be considered a catastrophe or a disaster. Likewise, a robust drought risk system is one that can resist frequent droughts (zero impact) and that can recover from socio-economic impacts of rare droughts (limited impact).

System robustness can be assessed by studying the so-called response curve: the relationship between disturbance and response, for example the relationship between precipitation deficit and crop yield deficit. Key aspects of the response curve for this application are:

  • Resistance threshold: the smallest deficit causing yield deficit;

  • Manageability: the severity of the yield deficit in comparison with a recovery threshold (i.e., the societally unacceptable level of yield deficit).

Manageability refers to how severe the impacts of a drought can become. Whether a response is considered severe depends on the social and economic capacity to recover from the response to a disturbance (see De Bruijn 2004). This capacity can be quantified by a recovery threshold, indicating the level of response from which it will be very difficult to recover. For example, when economic impacts of flooding exceed 5 % of the national GDP, it often means that aid is needed from other countries (5). Manageability indicates the extent to which the response stays far from the recovery threshold for the range of considered disturbance magnitudes.

4.2 Case Rijnland

Droughts occur as a natural temporary feature of the climate (AMS 2013), but it may have an effect on reservoir storage, groundwater resources and river discharges, thereby potentially causing socio-economic and environmental impacts. Climate change may cause a change in the frequency and severity of droughts. Many water resources studies focus on the impact of climate change on long-term average flows (Lehner et al. 2006), thereby analysing the average demand and supply balance in the long-term. However, even if water supply is sufficient to meet demand on average, droughts may cause a temporary shortage of water which may have far-reaching impacts on society. Insight is thus needed into the impact of potential drought events resulting from climate variability, now and in the future.

To understand how systems robustness analysis can aid in the management of drought risk systems, it was applied on a case: Rijnland polder areas (Mens and Klijn 2015; Mens et al. this issue). A drought risk system is a geographical area that depends on one or more water resource for several functions (e.g., drinking water supply, industry, nature and agriculture). Drought, in this case defined as a period of precipitation deficit (i.e., evapotranspiration is larger than precipitation), may cause water shortage with an impact on the area’s functions. Precipitation deficit is thus considered the disturbance to the system.

The case study area, Rijnland, is a system of low-lying polders with agriculture (potatoes, horticulture, bulbs, flowers, dairy farming) as the main land use type. This is a mixture of rain-fed and irrigated agriculture. The water supply and drainage system consists of ditches, canals, lakes, and pumps and sluices. In winter, there is an excess of water that is drained and pumped out, and in summer fresh water is let into the surface water system from a nearby river (see also Mens et al. this issue). The inlet water is used for water level control, flushing for water quality management and water supply (irrigation, drinking water, industry). Water level control guarantees the stability of canal embankments, and flushing reduces the salinity in the canals and ditches, caused by seepage of saline groundwater. During an average summer, the total external fresh water inlet amounts 40 to 60 106 m3. In a dry summer, this may increase to about 100 106 m3 (Rijnland 2009). The analysis focuses on the impact of drought on agriculture.

4.3 Model

Figure 3 presents an overview of the model structure. The water demand module generates water demands for irrigation and water level control in rural areas and is a simple two layer grid-based groundwater model with a resolution of 250 by 250 m, taking into account a limited number of land use and soil types. For each layer in each grid cell, the model calculates the water balance. First, the potential evaporation is calculated by multiplying the reference evaporation with a crop factor that is specified for each crop and 10-day period. The actual evaporation is a function of the potential evaporation, the moisture in the root zone, and the soil moisture suction (pF value). Lateral flow from groundwater to local surface water and vice versa is a function of groundwater depth relative to surface water level. Water flowing from the root zone to the subsoil (percolation) depends on the root depth, porosity, and precipitation. Capillary rise (flow from subsoil to root zone) is calculated as a function of the groundwater depth below the surface level and the root zone suction (Kabat and Beekma 1994; Oosterbaan 2001). The lower boundary condition of each plot is an annual seepage flux taken from results of the complex model for an average year. In case the root zone and subsoil are saturated, excess water is moved through surface runoff. In urban areas surface runoff is a function of the net precipitation and a runoff co- efficient of 0.8 (Urbonas and Roesner 1993). The water demand is determined from the difference between the actual and potential evaporation. The amount of water requested for maintaining the target water level in the local surface waters areas is derived from the net precipitation and the surface area of these waters. The grid cells are aggregated over a watershed area (called district).

Fig. 3
figure 3

Diagram showing the basic structure of the model (adapted from Haasnoot et al. (2014))

The salt intrusion module simulates the salt concentration at the Gouda inlet depending on river discharge and sea level. This module is based on empirical correlation between the Rhine discharge at Lobith and salt concentrations in the lower river reaches calculated using a 1D hydraulic model (SOBEK) (van den Boogaard and van Velzen 2012).

$$ \begin{array}{l} Salt=1700+\left(90-1700\right)\times \frac{e^{Fact}}{1+{e}^{Fact}}\hfill \\ {} Fact={\left(\frac{Q_{lobith}-600}{2.211}\right)}^{0.309}\hfill \end{array} $$

where Q lobith is the discharge at Lobith in cubic meter per second and Salt is the salt concentration at the Gouda inlet in milligram per liter. As discussed in Haasnoot et al. (2014), this relation will slightly underestimate the frequency of closure at Gouda.

The focus of this analysis is on the economic damages to agriculture due to drought. For this, we use Agricom (Mulder and Veldhuizen 2014) which is an agro-economic model to estimate agricultural yield losses due to water shortage, saline soil moisture and water excess. Drought is defined in terms of

$$ {E}_{ratio} = \frac{E{T}_{act}}{E{T}_{pot}} $$

Where ET act is the actual evatranspiration and ET pot is the potential evapotranspiration. Given E ratio and crop specific damage curves that take into account the growing season, the E ratio is translated into a loss of yield in kilogram, which in turn is monetized. In our analysis, we calculate both the potential yield in Euro, assuming perfect conditions, and the actual yield in Euro’s, allowing us to calculate the relative loss due to deficiencies in the system.

4.4 Results

The response curve of the reference situation is given in Fig. 4. Each point in the curve is a sample from climate variability according the current climate and according to the worst-case climate change scenario for 2100. Together they provide an indication of the range of drought magnitudes that may occur and the corresponding impact on agricultural yield deficit. The vertical range is explained by the timing of the precipitation deficit within the growing season. Some crops are more sensitive to drought in the beginning of the season, while others are more sensitive at the end of the season. Droughts with similar precipitation deficits may therefore occur, but when they differ in timing the impact will differ as well.

Fig. 4
figure 4

Response curve of the Rijnland system: yield deficit as a function of maximum cumulative precipitation deficit (MCPD), the indicator for drought magnitude

There are two lines drawn in the figure: one at 5 % yield deficit and one at 10 % yield deficit. Yield deficits below 5 % are considered not significant, because farmers take into account some level of business risk. The first point that exceeds this line (indicated with an arrow) indicates the resistance threshold at 250 mm precipitation deficit. The 10 % line is arbitrarily chosen to indicate the level of unacceptable yield deficit that must be avoided (recovery threshold). Only two of the 35 events exceed this 10 % line in the future climate. When it is aimed to enhance the manageability, measures should be designed that reduce the impact for these events.

In practice, the location of both the 5 % and the 10 % line is a political choice. In fact, the real recovery threshold in this system will never be exceeded, since it can be expected that the area can still economically recover from 100 % yield deficit (see Mens et al. this issue). The Rijnland system can be considered robust to droughts, because the maximum simulated yield deficit does not exceed 25 % for a large range of precipitation deficits. This can be partly explained by the storage capacity as well as the diversity of crop types that vary in when they are most sensitive to water shortage (Mens et al. this issue).

4.5 Reflection

Although the model obviously is a limited representation of reality, it suffices for the purpose of the robustness analysis. Different model parameter assumptions could slightly change the resistance threshold, but the score on proportionality and manageability is not expected to change. Thus, the conclusion on how well the system is able to deal with a range of drought events is not affected by the choice of the model.

Two types of measures are often considered in a drought risk context: those that increase the water supply and those that reduce the water demand. In traditional decision making approaches, where supply reliability is the main decision criterion, both types of measures will increase the supply reliability. However, a high reliability does not inform about the consequences of low-probability drought events. Additional criteria are thus needed that assess measures on how they affect the system in terms of avoiding large societal and economic impacts due to droughts. The two robustness criteria, resistance threshold and manageability, seem to fulfil this need. More applications are needed to explore how different types of drought risk reduction measures score on the robustness criteria.

We believe that a robust system, that can deal well with frequent as well as rare events, will also be able to deal well with changed variability under a future climate (Watts et al. 2012). Knowledge or assumptions about the future climate variability are not needed to perform robustness analysis. However, climate change scenarios are helpful in determining the range of drought events for which the system’s robustness is assessed.

5 An Exploratory Modeling Approach for Designing Adaptation Pathways

5.1 Adaptation Pathways

An emerging approach to the design of climate adaptation plans and strategies is adaptation pathways (Haasnoot et al. 2013; Wise et al. 2014). This approach combines two bodies of literature on planning under uncertainty: work on adaptive policymaking (Walker et al. 2001; Kwakkel et al. 2010; Hamarat et al. 2013) and work on adaptation tipping points and policy pathways (Kwadijk et al. 2010; Haasnoot et al. 2012; Offermans 2012). A plan is conceptualized as a series of actions taken over time. In order to come to a good plan, it is necessary to first identify candidate policy actions and their adaptation tipping points. An adaptation tipping point is the condition under which a given policy action no longer meets its objectives. The point in time at which this happens (the so-called sell-by-date) is scenario dependent. In light of an analysis of the sell-by dates of various policy actions, concatenations of options, or policy pathways, can be specified where a new policy option is activated once its predecessor is no longer able to meet the definition of success. Typically, there is a portfolio of pathways that decision-makers would like to keep open for the future. This adaptation map forms the basis for the plan. Figure 5 shows an example of such a map. For a more detailed elaboration on DAPP, see Haasnoot et al. (2013).

Fig. 5
figure 5

An example of an Adaptation Pathways map (left) and a scorecard presenting the costs and benefits of the 9 possible pathways presented in the map. In the map, starting from the current situation, targets begin to be missed after 4 years. Following the grey lines of the current plan, one can see that there are four options. Actions A and D should be able to achieve the targets for the next 100 years in all climate scenarios. If Action B is chosen after the first 4 years, a tipping point is reached within about 5 years; a shift to one of the other three actions will then be needed to achieve the targets (follow the orange lines). If Action C is chosen after the first 4 years, a shift to Action A, B, or D will be needed after approximately 85 years in the worst case scenario (follow the solid green lines). In all other scenarios, the targets will be achieved for the next 100 years (the dashed green line). The colors in the scorecard refer to the actions: A (red), B (orange), C (green), and D (blue)

Model-based decision support for the design of adaptation pathways is challenging for a variety of reasons. There exist a wide variety of candidate policies that has to be considered. These policies might be sequenced in different ways, taking into account co-dependencies between the various options. Their assessment needs to explicitly take into account the dynamics over time, rather than some endpoint. The assessment also should cover the variety of outcomes of interest the different stakeholders involved in the decision-making care about. And last but not least, the assessment should cover the variety of irreducible uncertainties that exist and the performance of candidate pathways needs to be assessed over the full range of uncertain factors.

One approach that might be amendable to supporting the design of adaptation pathways is exploratory modeling. Weaver et al. (2013) have argued that adaptive planning requires an exploratory modeling approach. In exploratory modeling, modelers account for the various unresolvable uncertain factors by conducting series of computational experiments that systematically explore the consequences of alternative realizations of the various uncertain factors (Bankes et al. 2013). These computational experiments can be used also a test bed for candidate policies (Lempert et al. 2003).

5.2 The Case, the Model and the Associated Uncertain Factors

As case, we use the Rijnland case as described in section 4. We use the same model as used there with a few small changes. First, we adapted the model to account for the following uncertain factors:

  • River runoff in the Rhine

  • Rainfall

  • Land use

  • Evatranspiration

For each of these sources of uncertainty we use transient scenarios (Haasnoot et al. 2014) to address them. So, rather than looking at the system in e.g., 2050, we consider the change over time from the current system to the state of the system in 2100. For river runoff, rainfall, and evapotranspiration we use different possible realizations of the climate scenarios W+ and G.. These are climate scnearios for the Netherlands and have been developed by the Dutch Royal metrological institute. W+ is an extreme warming scenario, while G represents a modest warming scenario (KNMI 2006). These realizations are the same as those used by Haasnoot et al. (2014), although here we only use the parts relevant to the case study area. We use 10 possible realizations of both scenarios. For land use, we use the land use as described in the four delta scenarios as developed in the Deltaprogram. These scenarios are labelled Warm, Pressure, Rest, and Steam. These delta scenarios combine the 2006 W+ and G climate scenarios (KNMI 2006) with socio-economic scenarios. Figure 6 shows the scenario logic for the scenarios. Maps were available for 2050 and 2100. We interpolated the maps in between. Given that we sum up over the region, a rather simple and crude interpolation has been used. Combining the different realizations with he delta scenarios gives 40 scenarios. To assess the role of changing land use, we include, in addition to the changing land use also a no change case, giving 60 scenarios in total (i.e., Warm, Pressure, Rest, Steam, No change W+, and no change G). Note that in the current results land use has some influence on water demand, but the irrigation maps do not evolve with land use.

Fig. 6
figure 6

The scenario logic of the delta scenarios (adapted from (Delta Programme 2012)

5.3 Results

A variety of analyses have been carried out. As a first experiment, we explore how the system responds if water supply is deliberately closed. To this end we simulate two situations. In the first situation, we keep the water inlet at Gouda open at all time, in the other, we close the inlet at all times. We assessed the performance for both situations by analyzing the results for 5000 experiments. A boxplot for the loss of income of farmers due to reduced yields in both situations is shown in Fig. 7. So, this is the monetized loss of yield as explained near the end of section 4.3. As can be seen, there are a few more extreme cases in terms of income loss when the Gouda inlet is closed at all times. Still, the differences are quite small and it is relevant to investigate this further.

Fig. 7
figure 7

Boxplot of percentage income loss for farmers due to a reduction in yields, comparing the situation where the fresh water inlet at Gouda is open at all times with the situation where the fresh water inlet at Gouda is closed at all times

In order to understand why there is only a minor difference in income loss between the situation where the fresh water inlet at Gouda is always open, and the situation where the fresh water inlet at Gouda is always closed, we looked at the worst-case delta scenario. In this scenario, warm, the system experiences severe climate change, combined with substantial land use developments putting even more pressure on the system. We explored the behavior of the system for all ten possible realizations. The resulting water demands per decade are shown in Fig. 8.

Fig. 8
figure 8

Water demand per decade for the Warm scenario. The left hand boxplot contains the water demand when the fresh water inlet at Gouda is closed because of salt intrusion. The right hand boxplot contains all the decades irrespective of whether the fresh water inlet at Gouda is closed or not

In Fig. 8 we compare the water demand at a decade level for only those situations where Gouda is closed on the one hand, and for all decades on the other. As can be seen, during the decades during which the fresh water inlet at Gouda is closed, the water demand is quite small. In fact, the demand is below the supply available via the small scale water supply (5 cubic meter per second). As such, no severe water shortages occur. This suggests that decades with high demand for fresh water do not coincide with low runoff off of the Rhine. In light of this analysis, it appears that no future adaptation actions are necessary. There are a few caveats with this conclusion.

First, the relation used to calculate closure of Gouda is known to underestimate closures (Haasnoot et al. 2014). As such, a closer analysis on a different time scale might be needed. The current analysis uses decades, and it might be necessary to shift to a daily analysis to get a better insight into what is happening.

Second, the model does consider land use change. This land use change affects water demand. We have not translated the changes in land use into changes in the areas being irrigated. Doing so is possible. Ter Maat et al (2013), for example, introduce several additional assumptions regarding future irrigation demands, allowing the translation of land use change into changes in water demand. We speculate that one of the reasons for the differences between the results reported here and those reported in ter Maat et al. (2013) are due to how changes in land use affect irrigation.

Third, it appears that the model at every time step resets the water levels to the norm height. Inspection of the source code of the model did not provide clarity as to where the water for this comes from. That is to say, it appears that a hidden assumption in the model is that that the full capacity of the KWA is available for maintaining water levels. Based on ter Maat et al. (2013), this assumption is questionable. If the model were modified such that the KWA is first used to maintain water levels, and only subsequently if possible for irrigation, it is plausible that there would be more severe fresh water shortages for irrigation. In turn, this would increase the income loss of farmers.

From a methodological point of view, the main conclusions of the analyses reported on are that it is possible to consider a wide range of different uncertainties simultaneously. Exploratory modeling offers an approach for systematically mapping out the consequences of various uncertainties. We also demonstrated that it is possible to do a variety of potentially useful further analyses on the exploratory modeling results. On the one hand, these analyses provide usefull information for policy-making. On the other hand, the analyses create new insight into the behavior of the models. Moreover, these kinds of analyses can also help in assessing the usefulness of models, since they stress test the model and clarify the domain over which the model provides meaningful results.

6 Conclusions and Discussion

Long term planning for fresh water supply faces a variety of sources and types of uncertainty, many of which cannot be adequately dealt with using probabilistic approaches. In response, a variety of complementary approaches has been proposed and developed to explicitly recognize and deal with these uncertainties. On the one hand, analytical methods are being developed to capture and assess uncertainties and their impacts, on the other, various action strategies are proposed in view of uncertainties, including approaches based on robustness, on resilience, and on flexibility and adaptation. We have illustrated the application of respectively a resilience approach, a robustness approach, and an exploratory modeling approach. Reflecting on these different approaches, we observe that the robustness approach requires quantification of uncertainties and computational analysis. The exploratory modeling approach builds on computational power and enables exploration of a broad set of possible assumptions about the future while relaxing requirements of probabilistic specification.

The resilience approach, on the other hand, is primarily qualitative. As a result, the effectiveness and efficiency of resilience options is difficult to assess in quantitative terms. But the resilience approach may, as illustrated, be more suitable for tailoring adaptation to local situations, and appeal to notions of stakeholder involvement and ownership, self-governance and self-relience. In contrast to the robustness and exploratory modelling approach, the resilience approach can be used for those adaptation policy challenges where the state of scientific knowledge does not (yet) allow a sufficiently adequate form of quantification of impacts at the spatial and temporal scales that are pertinent to that particular adaptation policy challenge.

Application in the Venen-Vechtstreek case illustrates that the resilience approach enables involvement of multiple stakeholder perspectives and consideration of multiple societal reactions (that in turn can influence the effectiveness of the various management options) as well as a range of other criteria relevant to the decision making challenge such as costs and co-benefits of the management options. In our workshop, the co-benefits turned out to be crucial to obtain support for resilience-based climate adaptation measures from local stakeholders such as farmers. The approach is especially useful for a quick screening of management options.

The robustness approach focusses on a partial challenge around a single metric (vulnerability of the present fresh water system from the view point of the water manager only), but enables spatially explicit detail and quantitative assessment of societal costs and benefits often required in policy making contexts to justify interventions that affect various stakeholders in various ways.

In practical terms, the resilience approach is useful in the phase where long term management and adaptation strategies for a region are being designed in multiple stakeholder processes and considering multiple impacts on multiple sectors. Robustness analysis is more useful to assess specific issues (vulnerability to a specific impact for a specific sector), and fit for situations where the relevant set of future scenario’s is known and quantifiable. As such the approaches can be highly complementary (see also Wardekker et al. 2010). Especially in cases where decisions on dimensioning of options is key, a combination of approaches can have added value: if a resilience based screening of options reveals the wish to create fresh water buffers in an area (buffering) or to increase the inlet of fresh water from outside the area (ominvory), then a specific robustness analysis can help to answer the questions how big the buffer should be or what is the required capacity of the extra fresh water inlet should be.

More generally speaking, in practical cases different types of uncertainties will play a role, and analyses and strategies at different levels will be needed. As a result, a tailored combination of different approaches may be the best way forward, combining more top-down measures with bottom-up approaches. Clearly, both the analytical approach and the choice of strategies will have to be tuned to the characteristics of the situation, and more research and experience are needed to increase our understanding of the advantages and disadvantages for various approaches, and how they relate to the various types of uncertainties.

Another lesson from the cases presented here is that the match between the choice of policy objectives (which possible impacts do we seek to reduce with the adaptation policies) and the choice of scope and system boundaries of the scientific assessment carried out to support that decision making process is crucial. Is the key concern in agricultural impacts the physical yield (in kg of crops per ha) or the economic yield (in monetary units per ha)? Should we focus on the robustness of water supply only (how can we maintain the fresh supply under future climates to meet the given future demand) or should we address the robustness of the entire fresh water system? In the latter case, adaptation of water demand to a reduced supply would also need to be considered in the scientific assessment. This relates to the question whose perspective should guide the scientific assessment and who is and is not deemed responsible for adaptation.