Next Article in Journal
Simulation of Groundwater-Surface Water Interactions under Different Land Use Scenarios in the Bulang Catchment, Northwest China
Next Article in Special Issue
Estimating Evapotranspiration from an Improved Two-Source Energy Balance Model Using ASTER Satellite Imagery
Previous Article in Journal
Spatial and Temporal Correlates of Greenhouse Gas Diffusion from a Hydropower Reservoir in the Southern United States
Previous Article in Special Issue
Spatio-Temporal Patterns of Water Table and Vegetation Status of a Deserted Area
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Dynamic, Multivariate Sustainability Measure for Robust Analysis of Water Management under Climate and Demand Uncertainty in an Arid Environment

by
Christian Hunter
1,2,3,*,
Jorge Gironás
1,2,4,5,†,
Diogo Bolster
3,† and
Christos A. Karavitis
6,†
1
Departamento de Ingeniería Hidráulica y Ambiental, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna, Santiago 4860, Chile
2
Centro de Desarrollo Urbano Sustentable (CEDEUS), Av. Vicuña Mackenna, Santiago 4860, Chile
3
Civil Engineering, Environmental and Earth Sciences Department (CEEES), University of Notre Dame, Notre Dame, IN 46556, USA
4
Centro Interdisciplinario de Cambio Global, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna, Santiago 4860, Chile
5
Centro de Investigación para la Gestión Integrada de Desastres Naturales (CIGIDEN), Av. Vicuña Mackenna, Santiago 4860, Chile
6
Department of Natural Resources Management & Agricultural Engineering, Agricultural University of Athens, Iera Odos 75, Athens 118 55, Greece
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Water 2015, 7(11), 5928-5958; https://doi.org/10.3390/w7115928
Submission received: 31 August 2015 / Revised: 10 October 2015 / Accepted: 19 October 2015 / Published: 30 October 2015

Abstract

:
Considering water resource scarcity and uncertainty in climate and demand futures, decision-makers require techniques for sustainability analysis in resource management. Through unclear definitions of “sustainability”, however, traditional indices for resource evaluation propose options of limited flexibility by adopting static climate and demand scenarios, limiting analysis variables to a particular water-use group and time. This work proposes a robust, multivariate, dynamic sustainability evaluation technique and corresponding performance indicator called Measure of Sustainability (MoS) for resource management that is more adapted to withstand future parameter variation. The range of potential future climate and demand scenarios is simulated through a calibrated hydrological model of Copiapó, Chile, a case study example of an arid watershed under extreme natural and anthropogenic water stresses. Comparing MoS and cost rankings of proposed water management schemes, this paper determines that the traditional evaluation method not only underestimates future water deficits, but also espouses solutions without considering uncertainties in supply and demand. Given the uncertainty of the future and the dependence of resources upon climate and market trajectories, the MoS methodology proposes solutions that, while perhaps are not the most optimal, are robust to variations in future parameter values and are thus the best water management options in a stochastic natural world.

Graphical Abstract

1. Introduction

In the 20th century, the world population has tripled, but water use has increased six-fold, indicating augmented stresses on limited resources [1]. In light of this, politicians and decision-makers seek solutions to meet demand [2]. One approach is a shift to include both supply and demand management through Integrated Water Resources Management (IWRM) (e.g., [3,4,5,6,7,8,9]). Of growing importance in IWRM, one such method has been to improve resource efficiency through watershed-scale application of water reuse and recycling (WRR) technologies [10,11,12].
Equitable water use as proposed by IWRM faces challenges due to continuous structural transformations in society and the environment [13], yielding problems of high complexity and uncertainty [14,15,16]. Climate change will lead to uncertain future water supplies [17], while population growth, technological innovation, composition of future economies, and dynamic stakeholder goals add further demand-side uncertainties [18]. As such, methodological tools are required that enable the assessment of sustainability under an uncertain future.
It is well known, however, that defining sustainability in a rigorous and meaningful manner is challenging at best [19,20]. The first widely accepted articulation of the term sustainable development, published in the Brundtland report, calls for a “development that meets the needs of the present without compromising the ability of future generations to meet their own needs” [21]. Since its publication, close to 300 other attempts at a definition of “sustainability” have been published (e.g., [19,22,23,24]), none of which is universally accepted, but there is increasing consensus that discussions thereof should include uncertainty [25]. Under uncertainty, we must seek solutions that may not be the most optimal under present conditions, but are the most robust (i.e., less likely to fail and more immune to the effects of failure due to factors of uncertainty) [26,27]. We argue that any robust measure of water resources sustainability must consider three major sources of uncertainty under a decision support analysis as represented by the following questions [20,22]: (1) for whom is the resource being sustained; (2) at what time will sustainability be determined; (3) under what future conditions will the resources remain sustainable. From hereafter the sources of uncertainty associated with these questions will be respectively referred as to Whomever, Whenever, and No Matter What.
In line with this, for over three decades there have been attempts to create water risk indices, easy for lay-users to understand but capable of merging the intricacies of resource scarcity into a single unified number, a “synthesis of numerous factors into one given factor” [28]. Decision support tools in water resource management often aim to do this (e.g., [29,30,31,32,33,34,35]). For example, to quantify risk in a reservoir to optimize water management under uncertainty and scarcity, Hashimoto et al. proposed a landmark methodology, called RRV, that quantifies risk in terms of the frequency (Reliability, R), duration (Resilience, R), and magnitude (Vulnerability, V) of system failure [29]. RRV has been applied in numerous hydrologic settings to quantify water quantity and quality risk under varying scenarios (i.e., [13,36,37]). Loucks [31] and Sandoval-Solis et al. [35] unified the RRV factors, proposing multi-criteria analysis methodologies that consider not only the needs of various stakeholders but also incorporate multiple risk measures.
These methodologies, however, are similar to traditional analyses techniques in that, while addressing the first two sources of resource uncertainty (Whomever and Whenever) through dynamic multi-criteria analysis, their simplification of variable trajectories overlooks the resources’ supply and demand futures (No Matter What). By assuming unchanging climate parameters and resource use, these previous studies are susceptible to adopting solutions that, while optimal for present conditions, are incapable of responding to variations in these. Overlooking uncertainty in future conditions risks the possibility of espousing solutions unsustainable in the very future for which they are proposed.
The objective of this study is to propose a methodology that encompasses multiple sources of uncertainty in decision-making, including climate change, population growth, and demand trajectories defined by future socio-economic scenarios. While traditional engineering solutions operating under deterministic parameters often invoke “optimization” as the problem-solving goal, in stochastic or uncertain situations “robustness” must serve as the objective so as to maintain performance of the system when subjected to unpredictable perturbations [38]. By evaluating the performance of a dynamic and multivariate risk index in response to the range of potential supply and demand futures, this methodology seeks robust sustainable options under future uncertainties [39]. In our approach, we incorporate climate change projections, based on downscaled Global Climate Projections (GCP) following Special Report on Emission scenarios (SRES) and Representative Concentration Pathway (RCP) scenarios, as well as growth projections provided by experts in each water sector, considering conservative, average, and extreme demand scenarios. These are integrated into a Water Evaluation And Planning (WEAP)-MODFLOW coupled hydrology model that simulates the surface-subsurface system. To assess future outcomes we propose a new dynamic multidimensional evaluation methodology for sustainable water management plan ranking and selection.
This new methodology for robust sustainability analysis is presented in Section 2 of this paper, and in Section 3 is applied to the Copiapó River basin in Chile, one of the most arid regions in the world in which strong economic and population growth is taking place due to agricultural and mining industries. A recent study of this site [40] explored several water management strategies, calling for a reduction in demand of 20%–30%, a range achievable in practice through effective WRR [41]. This study, while providing useful planning information, adopted the oversimplifying assumptions common to current water management analyses of constant supply and demand futures, single variable analysis, and static rather than dynamic evaluation of sustainability (e.g., [42,43,44]). Results from the application of our new dynamic, multivariate, robust sustainability measure are compared in Section 4 with those using a static, single-variable, deterministic future (i.e., [40]) to demonstrate the importance of considering these often overlooked factors in resource decision-making. Finally, Section 5 presents the main conclusions of this work.

2. Materials and Methods

2.1. Previous Sustainability Methods—RRV and Sustainability Index

Previous studies have made important contributions to the field of sustainability analysis. Recognizing the dynamic nature of risk related to water resource management, Hashimoto et al. developed the RRV (Reliability, Resiliency, Vulnerability) methodology [29]. Reliability describes the probability a system is in a satisfactory state, resiliency is a measure of how quickly a system recovers from failure, and vulnerability is the expected magnitude of failure. For application in IWRM, consider two time series of water demand XD(t) and supply XS(t), which can be broken into n equal time increments (t = 1,2,..,n). We define a water deficit during an interval t, as D t = ( X D , t X S , t ) H ( X D , t X S , t ) , where X represents amount of water and H ( x ) is the Heaviside step function, which assigns a value of D t = 0 when supply exceeds demand, which we say represents a satisfactory state (S). Likewise, a state of failure (F) is defined as when demand exceeds supply, that is X D , t >   X S , t or Dt > 0 =   ( X D , t X S , t ) . RRV for a certain variable i, where represents the number of occurrences, are given by:
R e l i a b i l i t y i = α i = P [ X t S ] = N o ( D t = 0 ) n
R e s i l i e n c e i = ρ i = P { X t S , X t 1 F } = { N o ( D t = 0 ) | D t 1 > 0 }   N o ( D t > 0 )
V u l n e r a b i l i t y i = ε i = E [ X t | X t F ] = t = 0 t = n D t N o ( D t > 0 )
While these three metrics are individually informative, Loucks [31,45] proposed combining the RRV measures in such a way as to define a sustainability index (SI). In particular he proposed the following index S I L , given by:
S I L i = α i ρ i ( 1 ε r e l i )
Note that S I L includes the concept of relative vulnerability, ε r e l , which corresponds to a normalized measured of the maximum vulnerability within the context of all the management alternatives being compared for sustainability:
ε r e l i = ε i / ( M a x A l l ε )
S I L allows for the comparison of relative sustainability between competing management options, simplifying it down to a single number for ease of use by water resource managers [46,47]. Indeed, S I L summarizes the performance of alternative policies from the perspective of water users and the environment, and it is also a measure of a system’s adaptive capacity to reduce its vulnerability [35]. It ranges from 0 to 1 with 1 indicating complete sustainability. If any one measure has a value of 0, the sustainability will also be 0 (i.e., no complete tradeoff is possible) and its multiplicative form implicitly weighs heavily the worst performing criteria.
Sandoval-Solis et al. [35] suggest improving the original index in structure, scale, and content to make it more flexible and adjustable to the requirements of each basin. Thus, they add two new dimensionless resource performance criteria to Louck’s original SI, which are integral to resource planning: standard deviation (σ, Equation (6)), a measure of variability, and maximum deficit (Equation (7)), a measure of worst-case scenario:
σ i = [ t = 1 t = n ( X S , t i X ¯ S i ) 2 ] / ( n 1 ) X D , t i
δ m = m a x ( D A n n u a l i ) X D , t i
Moreover, Sandoval-Solis et al. [35] offer a generalized form of the SI calculation, allowing for the inclusion of other criteria depending upon the user’s perception of the necessities of each case:
S I i = [ m = 1 M C m i ] 1 M
where S I i is the sustainability index for variable i, calculated using M performance criteria C m i . This approach solves a problem of results scaling, as the geometric average of the criteria values still penalizes poor performance, but is more representative (e.g., three criteria values of 50% each would have a logical 50% final rating using Sandoval-Solis, compared to a 12.5% using Loucks).
As might be seen in the above equations, RRV analysis (and thus the SI as well) meets the uncertainty criteria of Whenever, using dynamic rather than static performance measures independent of when they are evaluated in the delimited time frame n. Also, the index may be calculated for more than one watershed variable so as to include SI calculations for a range of variables from various stakeholder sectors, thus meeting the uncertainty criteria of Whoever stated above as well. Indeed, we initially propose to consider the agriculture, industry, and potable water production sectors, although other stakeholder groups of relevance for the watershed under study may also be considered. However, the third sustainability criterion (i.e., No Matter What), which deals with uncertainty in future variable values, has been largely overlooked (No Matter What). To meet this criterion, climate change and future demand scenarios caused by uncertainties of markets, values, and population must be included. Overall, decision-makers wish to select the management option that, although possibly not optimal, is robust to uncertain future variables and thus able to withstand the unforeseen [20].

2.2. Proposed Sustainability Methodology—Measure of Sustainability (MoS)

We propose a sustainability measure that continues to address the first two uncertainty questions, maintaining the multisectorial approach proposed by Loucks to address the environmental, political, social, and economic impacts of water resources [31]. In addition, we retain the five-factor RRV sustainability index presented by Sandoval-Solis et al. [35], but this research advances the state of present sustainability evaluation in that it addresses uncertainty in future parameters under which decisions will be forced to operate. By considering a range of possible supply futures, and using climate change models as well as a cone of potential demand trajectories, this research proposes a robust methodology for water management decision-making. To differentiate from previous indices, we shall refer to this method as the Measure of Sustainability (MoS).

2.2.1. Future Water Demand

Given growth observed in recent years not only in population but in per capita consumption and in economies, static projections can be both erroneous and lead to inadequate long-term planning. While future conditions might experience a wide range of potential trajectories, econometrics in resource planning and sociological literature lead us to consider the three different logical growth options presented in Table 1 to estimate the population in year t, Pt, given the initial population in the reference year 0, P0. Throughout its life cycle, any socio-economic sector may experience all three of these different growth patterns [48]. During the early stages of sectorial growth, rapid change is modeled according to a uniform percentage growth rate r (exponential/extreme growth). Later, this rate steadies to a constant K (linear/average growth). Finally, as the sector nears its saturation level, the rate of growth continues to decrease with limit at 0 (cubic/conservative growth). Given the difficulty of determining in which of these three stages present conditions lie, we consider all three possibilities simultaneously through a “demand cone” (Figure 1) which, starting with present values, expands into the future to consider the range of possibilities. Demand futures in potable, agricultural, and industrial/mining water will be considered separately, each using all three of the demand cone possible scenarios: extreme, average, and conservative growth (Figure 1, Table 1) [49,50].
Table 1. Generic future demand models according to life cycle stage.
Table 1. Generic future demand models according to life cycle stage.
Life Cycle StageGrowth ScenarioModelEquation
EarlyExtreme (Ext)Exponential   P t = P 0 e r t
MiddleAverage (Avg)Linear P t = P 0 + K t
LateConservative (Con)Declining Growth P t = P 0 + α 3 t 3 + α 2 t 2 + α 1 t
Figure 1. Generic “demand cone” of future trajectories; extreme (Ext), average (Avg), and conservative (Cons) futures.
Figure 1. Generic “demand cone” of future trajectories; extreme (Ext), average (Avg), and conservative (Cons) futures.
Water 07 05928 g001

2.2.2. Future Climates

Future supply variation must also be considered in the analysis, as it might provide an even greater source of uncertainty. Both Loucks [31] and Sandoval-Solis et al. [35] recognize the need for a sustainability analysis that encompasses climate change (e.g., [51]), but neither attempts it. Given documented potential changes in global and regional climate patterns that affect both temperature and precipitation distribution, it is crucial for resource planning to consider climate change.
Given the large grid size (over 1° in latitude and longitude) of current GCP´s compared with watershed scales, GCM projections must be downscaled for use in hydrologic supply analysis. For smaller watersheds such as that of Copiapó, Chile, a single-climate model grid cell might be larger than an entire watershed, covering elevations between 0 and 6000 m, with average annual precipitations ranging from 0 to 500 mm. Thus, the predictions of the climate models must be adjusted so as to be applied to discrete weather stations located within the watershed. Following the quantile mapping and bias correction methodology of Wood et al. [52] that is used by Meza et al. [53], the downscaling process adopted in this study first attempts to correct for a spatial bias, using the inverse distance squared method [54] between the four grid points nearest the water station of interest to determine a geometric average. Then, comparing the probability distribution of the model average for the control period (2000–2010) with the probability distribution of historical data, the methodology adjusts for elevation and seasonal biases by shifting the Probability Distribution Function (PDF) of the model vertically and stretching horizontally to align with the mean and standard deviation of the historical data. This same shift, then, is applied to the climate projection data for future temperature and precipitation. While this methodology does assume that the bias between the historical and control data shall remain constant in the future, the factors causing the bias (location, elevation, seasons) typically change slowly enough not to be a concern [52].
Dozens of climate models exist, which in turn consider different emissions scenarios, resulting in hundreds of potential hydrologic supply futures. Users of this methodology must first select a period for evaluation; we have chosen 30 years as a typical period of capital investment. Calculating the 30-year changes in average annual temperature ( Δ T , in °C) and annual precipitation ( Δ P , in %) for different climate models under a variety of emission scenarios for the watershed, results are plotted as may be observed in Figure 2. Given that we are not specifically interested in the performance of the system under each and every climate model (blue points, Figure 2), but rather in the range of possible futures, we choose to limit the volume of data by selecting the four extreme climate models (dry-hot, dry-cool, wet-hot, wet-cool) as well as an average model to represent the variety of possible climates (red points, Figure 2). The dotted vertical line in Figure 2 represents 2.5 standard deviations above the mean for percent change in precipitation. While we have chosen to exclude outliers that fall to the right of this line (p = 0.006), this determination is made by the user of the methodology. Indeed, the election of the four extreme climate scenarios is subjective, and if desired, confidence intervals could be used depending on the decision-maker’s perspective. For the final MoS calculation, results for these five models X i = ( X i , P , X i , T ) will be weighted based upon the percentage of all climate models M j = ( M j , P , M j , T ) closest to each on the Δ T   v s .   Δ P plot (e.g., Figure 2 below for the Copiapó, Chile example). For each climate model M j , Euclidean distances are calculated between the model and each of the five representative models X i . The model M j will be represented by the closest X i .
Figure 2. Shown is Δ T   v s .   Δ P , 30-year average percent change in climate (Copiapó, Chile example). Blue points indicate all climate models with emission scenarios while red points are the chosen representative scenarios. Vertical dotted line indicates cut-off for outliers (>2.5 standard deviations above mean).
Figure 2. Shown is Δ T   v s .   Δ P , 30-year average percent change in climate (Copiapó, Chile example). Blue points indicate all climate models with emission scenarios while red points are the chosen representative scenarios. Vertical dotted line indicates cut-off for outliers (>2.5 standard deviations above mean).
Water 07 05928 g002

2.2.3. Computation of the MoS Index

For robust analysis that addresses the three previously mentioned major questions of uncertainty, each management plan is evaluated using the RRV-based five-factor SI presented by Sandoval-Solis et al. [35] (Equation (9)) for the four different water use dimensions (environmental, economic, social, political) under five different climate models and three different demand permutations, yielding one final MoS score (Figure 3).
Following is a list of generalized steps for the calculation of the MoS, applicable to any study site (Figure 3):
  • Identify the five climate models and three demand trajectories using the process described above for a total of 15 combinations of potential futures;
  • For each of these 15 futures, collect time-series data for variables representing the four water use dimensions (environmental, political, social, economic), giving 60 variables to be considered over a predetermined time frame;
  • For each of these 60 variables i, calculate sustainability according to:
    S S S i = [ α i ρ i ( 1 ε r e l i ) ( 1 δ M i ) ( 1 σ i ) ] 1 / 5
  • For each potential future, combine the S I S S i scores for the four water use dimensions through a weighted sum by means of user-defined weights of importance to give a future demand score (default: equal importance [31]);
  • For each climate model, combine the future demand scores from Step 4 for each of the three demand trajectories by using user-defined weights of potentiality (default: equally probable) (red diamond in Figure 3);
  • Calculate the final index score (MoS) by combining the five S I S S i scores from Step 5, weighted according to the method described in the previous section, based on the percentage of all climate models closet to each simulated future on the Δ T   v s .   Δ P graph. (Figure 2).
Figure 3. Flow chart of Steps 1–5 of MoS methodology for Avg. climate scenario. Step 6 combines all climate scenarios using weights from Δ T   v s .   Δ P plot for final MoS score. Step numbers correspond to the list of steps outlined below.
Figure 3. Flow chart of Steps 1–5 of MoS methodology for Avg. climate scenario. Step 6 combines all climate scenarios using weights from Δ T   v s .   Δ P plot for final MoS score. Step numbers correspond to the list of steps outlined below.
Water 07 05928 g003
To measure performance in the water use dimensions mentioned in Step 2 of the MoS methodology above, decision-makers may select any dynamic variable represented by a time-series that they feel exemplifies the state of water in relation to the dimension. This paper, however, proposes the variables observed in Table 2, with their respective failure thresholds, as simple but succinct measures of water dimension “health”. The economic dimension will be gauged by the sum of the percentage of water demand met for each productive stakeholder, weighted by that member’s percent participation in the watershed economy. Economic failure shall be considered below 100% of weighted water demand, as healthy economies generally do not allow for deflation. Given the crucial importance of domestic potable water for health, comfort, and hygiene, the threshold for the social variable of the percentage of potable water demand met will be set at 100%; this threshold is a metric common to studies of potable water supply [55]. The political impact of water shall be considered in a state of failure when less than 75% of stakeholders have their water rights met. With more than 25% of water rights left unsatisfied, growing discontent among voters begins to make political re-election more difficult for a single decision-maker in a democratic watershed (inclusion threshold of competitor) [56].
Table 2. Sustainability variables and failure thresholds for water use dimensions.
Table 2. Sustainability variables and failure thresholds for water use dimensions.
Water Use DimensionVariableFailure Threshold
EnvironmentalAquifer levelSlope (current-critical)
Economic% demand economic weightLess than 100%
Social% of potable demand metLess than 100%
Political% of all demands metLess than 75%
Special mention must be made for the environmental water variable, aquifer volume. While the other water use dimensions are evaluated according to variations on demand met, and thus not dependent on whether a zone is arid or not, the variable chosen to evaluate the environmental dimension (aquifer volume) is particular to arid basins. In more water-rich environments, a suitable variable to measure the health of the environmental dimension might be stage height of a river or total volume of a reservoir; such measures, however, are often not appropriate for arid zone hydrology. In these areas, groundwater often provides a large proportion of the resource and thus aquifer volume serves as a reasonable evaluation variable. Nevertheless, the general methodology proposed would be equally applicable under more humid scenarios given appropriately defined evaluation variables.
For arid basins, while a certain depth to groundwater could be set as the failure threshold, due to overuse of groundwater compared with the relatively slow recharge process in water-scarce zones, true sustainable use of these resources (that is, maintaining resources as constant) is quite simply not feasible. Recovery once surpassing a static threshold of depth, will most likely not occur unless demand is reduced by unmanageable amounts. “Sustainability” in this case shall be defined as the continual presence of water resources during the modeled time period in more than half of the current wells in the watershed. Using the example of Copiapó (Chile) presented later, this threshold shall be surpassed when the aquifer level falls more than 60 m below the land surface, or rather when the total volume drops below 7373 million m3. In some sense this approach defines sustainability as maintaining a sufficiently slow decline in resources with respect to time, in the hope that over several decades alternative unidentified solutions will emerge that can address declining resources effectively. While optimistic, it is likely the only currently realistic option. Thus, for similar cases in arid regions where withdrawal greatly exceeds recharge, the environmental dimension shall be evaluated according to a sloping linear threshold between present volume levels and the critical volume expected at the end of the simulation period (i.e., a slope of −42.54 million m3/year for the Copiapó case study to be illustrated below). While the relative variables and failure thresholds for these four water use dimensions are all interchangeable according to the values and opinion of the decision-maker, we believe that these proposed variables and failure thresholds represent the logical functioning of a sustainable watershed.

3. Application to a Case Study

To demonstrate the application of this generalized methodology for a robust, dynamic, multi-criteria sustainability index, we utilize the example of the Copiapó River Basin (18,538 km2), located in the Atacama Desert of Chile (Figure 4, divided into six administrative sectors) and whose main urban center is the city of Copiapó (elev. 391 m, pop. 158,000, sector 4a). This snowmelt-dominated watershed, situated between 27° S and 29º S latitude, represents a basin under extreme water stress. Positioned in an arid region (28 mm average annual precipitation over the entire watershed, but with large inter-annual variation due to El Niño and large spatial variation due to elevation change) in which constantly growing demands increasingly strain scarce water resources, the Copiapó River Basin experiences an overexploitation of water resources by 30%–40%. Mining, agriculture, water utility, and tourism all compete for access to limited resources; indeed, water rights already granted are estimated to represent almost 500% of the groundwater recharge rate. The mining sector (copper, gold, iron) owns 17% of these rights while the agricultural sector holds 56% (mostly grapes, olives, and vegetable crops) [40]. The recent gifting of extraction rights by the Candelaria mine to the Colla indigenous community has skewed these percentages even further towards agriculture [57]. Growth in both activities, however, is considered fundamental to regional and national economic development with 33 new mining projects (est. $50 billion USD), and a 20% increase in agricultural surface area projected over the next five years [58]. Additionally, economic development has led to an increase of local potable consumption (8.8% of water rights) [40]. This water stress has left the Lautaro Reservoir empty during most months and is forcing decision-makers to look elsewhere for solutions. Further physical and operational details of the Copiapó watershed may be found in Suarez et al. [40].
Figure 4. The Copiapó River, with an average flow of ~1.5 m3/s, largely from snowmelt, drops from 1230 m to sea level over a distance of approximately 162 km. Here it is divided into six administrative sectors (Source: [59]).
Figure 4. The Copiapó River, with an average flow of ~1.5 m3/s, largely from snowmelt, drops from 1230 m to sea level over a distance of approximately 162 km. Here it is divided into six administrative sectors (Source: [59]).
Water 07 05928 g004

3.1. Watershed Modeling in WEAP-MODFLOW

The literature review, government reports, and fieldwork data utilized to construct the original AQUATOOL model presented by Suarez et al. [40] served as the basis for the construction of our hydrologic model, built with the Water Evaluation and Planning System (WEAP) water accounting platform [60]. A “user-friendly software tool that takes an integrated approach to water resources planning” (www.weap21.org), WEAP allows for the visualization of water supply and demand, and permits the input of necessary resource parameters (rainfall, evaporation, demand, infiltration, etc.) for integrated watershed simulation with monthly time steps. This platform was chosen over the previous AQUATOOL model because of its growing popularity in the literature, its ability to incorporate existing MODFLOW groundwater models into a coupled surface/groundwater management scheme (such linkage was not considered by Suarez et al. [40]), its ability to simultaneously compare various scenarios, and its free distribution for academic purposes in developing nations.
Using Geographical Information System (GIS) layers created with information from the DGA Santiago, Chile, Chilean Water Authority), SERNAGEOMIN (Santiago, Chile, Chilean National Geologic and Mining Service), the Junta de Vigilancia del Copiapó (Copiapó, Chile, Water Management Committee), Aguas Chañar (Copiapó, Chile, The Water Utility), DICTUC (Santiago, Chile, Engineering Firm Associated with the Universidad Católica de Chile), and local experts, a 33 catchment WEAP model with 13 mining and potable water demand sites was crafted and coupled with MODFLOW to incorporate groundwater flows. This MODFLOW model, originally provided by SERNAGEOMIN (the Chilean National Geologic and Mining Service), was cross-referenced with field well data to assure accuracy [59]. Rainfall data was acquired from the DGA for 13 different rain gauges with more than 10 years of data located throughout the watershed, covering an elevation range between 15 and 2250 m: Canto de Agua, Conay en Albericoque, Copiapó, Iglesia Colorada, Lautaro Embalse, Los Loros, Caldera, El Totral, Elibor Campamento, Jorquera en la Guardia, Las Vegas, Hacienda Manflas, and Pastos Grandes. Of these, the first six also provided historical temperature information for the calibration period (1993–2012). Due to the greater spatial variability in precipitation compared with that of the temperature, fewer weather stations are needed for an appropriate representation of the latter. Reservoir and canal operation parameters were provided by the Junta de Vigilancia. Historical demand information for the industrial/mining and potable water sectors were obtained from DICTUC [59] and cross-referenced with information provided by COCHILCO [61] and Aguas Chañar, respectively (A. Sola, personal communication, 11 June 2012). Agricultural demand is calculated directly by WEAP using the crop coefficient method multiplied by the amount of land sown [62]. This method requires the input of variables pertaining to evapotranspiration including crop coefficients, effective precipitation, agricultural acreage, and reference evapotranspiration or ET0 [59,63].
While most parameters for the proper calibration of WEAP may be directly input from data, this last factor (ET0) requires additional analysis due to changes in future climates. To maintain parsimony for the model, we compute ET0 using the Thornthwaite model [64,65,66], which is one of the simplest and yet most widespread empirical methods for evaporation estimation [67]:
E T 0 = { 16 ( L 12 ) ( N 30 ) ( 10 T I ) α   0 T 26 o r 415.85 + 34.24 T 0.43 T 2   T >   26
where
I = n = 1 n = 12 ( T n / 5 ) 1.514
and
α = ( 6.75 × 10 7 ) I 3 ( 7.71 × 10 5 ) I 2 + ( 1.792 × 10 2 ) I + 0.49239
In these Equations, L is the average daylength in (h) for the month being calculated, N the number of days in the month, and T the average daily temperature ( ).
The constant underestimation of E T 0 , observed in Table 3 when compared to pan evaporation ( E p a n ) using the Thornthwaite Equations (10)–(12), we believe to be explained by Pruitt’s [68] “oasis factor” correction of 1.76 (calculated as E p a n / E T 0 ). This factor accounts for desert air surrounding a source of evaporation and acting as a near-infinite sink for humidity ([69,70,71,72]). Although the average calculated “oasis factor” for the Copiapó basin (1.7) was very similar to Pruitt´s proposed 1.76, some extreme values were as large as 2.2. Large values such as these, however, are acceptable for the Atacama desert as a previous study in the region reported temp/pan evaporation pairings with an average oasis factor as high as 2.48 [73]. Monthly adjustment factors were calculated using Thornthwaite (Equations (10)–(12)), as were average monthly values and observed monthly E p a n for each of the six temperature gauge stations in the Copiapó basin. These were verified by comparing “oasis factor”-adjusted Thornthwaite estimates E T O a s i s with historical observations, yielding a relatively high correlation (i.e., a Nash-Sutcliffe coefficient (Equation (14)) of 0.52, and a coefficient of determination of 0.72). We thus believe that the adjusted Thornthwaite evaporation values E T O a s i s ) provide a simple and rather accurate method for forecasting evaporation and ET0 for use as a WEAP variable. Hence, while historical precipitation values were used for calibration, Thornthwaite evaporation values calculated using historical temperature and “oasis factors” were used for the calibration period so as to maintain continuity between calibration and experimental periods.
Table 3. Observed and calculated evaporation, with “oasis factor” values at latitude (28º S) for the Lautaro reservoir station.
Table 3. Observed and calculated evaporation, with “oasis factor” values at latitude (28º S) for the Lautaro reservoir station.
MonthJanuaryFebruaryMarchAprilMayJuneJulyAugustSeptemberOctonerNovemberDecember
L (h)13.512.912.111.310.610.310.41111.712.613.313.7
T ( ) 20.420.419.517.315.113.412.51415.116.317.519.1
Epan (mm)186154133786447517295131156188
ET0 (mm)100.686.983.861.847.736.433.943.451.064.173.991.2
“Oasis Factor”1.851.771.591.261.341.291.501.661.862.042.112.06
Daily ETOasis (mm)6.05.54.32.62.11.61.62.33.24.25.26.1
Another factor difficult to enter directly for future-casting, the loss to groundwater in the Lautaro reservoir clearly depends upon the volume of water contained therein. WEAP-MODFLOW, therefore, must be programmed to estimate this value based upon reservoir volume. Using bathymetry of volume vs. stage height provided by the Junta de Vigilancia, a cubic regression ( R 2 = 0.995 ) was calculated using historic infiltration data, resulting in the following relationship between stage height h (m) and reservoir infiltration RI (m3):
R I = 485.03 h 3 23018 h 2 + 411710 h 2174040
Transmission canal leakage and losses, a major “user” of water resources, must be estimated by the WEAP-MODFLOW user. A 10%–40% volume loss during transmission is reported as typical in arid regions [74,75]. The Chilean Sanitary Service (SISS) [76] reported that 40.4% of the water treated and delivered by the water utility company in Copiapó, Aguas Chañar, did not reach final users (with 35.4% being the Chilean national average). While the same report defines an efficiently operating system as one losing only 20% of water during transmission, the historic and continual growth of losses (only 29.5% in 1998) combined with the outlook of the SISS considering improvement possibilities lead us to propose calibrating using the 40% loss rate but simulating the future using the national average of 35% to represent proposed improvements in efficiency.
We calibrated the model, using the Nash-Sutcliffe coefficient (Equation (14)), considering both surface and groundwater variables [77].
E = 1 = 1 N ( O i P i ) 2 i = 1 N ( O i O ¯ ) 2
In this equation, O ¯   represents the observed mean of the variable in question, O i   the observed value, P i   the value predicted by the model being tested, and N the number of observations. Considering E values at the 10 different river gauge stations to evaluate surface water modeling performance and E of total aquifer volume for groundwater performance, the WEAP-MODFLOW model was compared with that created by AQUATOOL, presented in Suarez et al. [40]. Of the 10 river sectors, nine performed better in the WEAP-MODFLOW model, with the lower accuracy at the Copiapó city station where river withdrawals are insignificant to non-existent (Table 4). The only river sector in WEAP-MODFLOW exhibiting poor performance (Canal Mal Paso), still outperformed the AQUATOOL model; the difficulty in modeling this particular canal may be attributed to known undocumented (illegal) water withdrawals and inconsistent management. For groundwater volume simulation, the WEAP-MODFLOW model exhibited highly accurate predictive power (Table 4). Hence, the WEAP-MODFLOW model more accurately predicts performance and historical time dependence of both the river and groundwater, and thus may be considered at least as appropriate as, if not more than, the previous study (i.e., [40]) for our purposes of future simulation for sustainability analysis [78]. Moreover, the new model allows a more physically based representation of the groundwater as the previous model used a more conceptual approach.
Table 4. Comparative calibration, Nash-Sutcliffe, and Q ¯ ; stream flow and aquifer volume.
Table 4. Comparative calibration, Nash-Sutcliffe, and Q ¯ ; stream flow and aquifer volume.
River GaugeObs. FlowWEAP-MODFLOWAQUATOOL
Q ¯ obs (m3/s)Nash (E) Q ¯ cal (m3/s)Nash (E) Q ¯ cal (m3/s)
Manflas1.241.001.150.991.18
Jorquera0.620.970.670.970.60
Pulido1.600.991.550.991.55
Pastillo2.440.922.350.902.26
Emb. Lautaro2.440.742.33−1.450.42
La Puerta2.200.672.120.491.91
Mal Paso0.440.470.870.030.44
Canal Mal Paso0.70−0.330.74−7.730.44
Copiapó0.290.490.350.570.31
Angostura0.290.360.090.070.03
Aquifer Volume0.98-0.44-

3.2. Input of Future Scenarios

For future simulation, WEAP requires the input of future supply and demand conditions. Extreme, average, and conservative demand projections for the three different water use sectors (i.e., agriculture, mining/industry, and potable water) were estimated based on the evolution of historical data. Table 5 presents the analytical expressions as functions of time t for these projections while Figure 5, Figure 6 and Figure 7 graphically show these projections for each sector with irrigated surface area serving as a representative proxy of the growth of the agricultural sector; demand growth in mining/industry and potable water is expressed in m3/s.
Table 5. Demand scenario future models where t is the calendar year; response variables are million hectares for agriculture, m3/s for mining/industry, and population for potable water.
Table 5. Demand scenario future models where t is the calendar year; response variables are million hectares for agriculture, m3/s for mining/industry, and population for potable water.
Demand SectorExtremeAverageConservative
Agriculture 9 × 10 20 e 0.0293 t 68677 t 10 8 8.01 t 3 + 48017 t 2 10 8 t + 6 × 10 10
Mining/Industry 10 98 e 0.1124 t 0.0227 t 44.837 0.000004 t 3 0.0256 t 2 + 53.12 t 36738
Potable Water 269.2 × 10 13 e 0.0206 t 134.6 × ( 3275 t 6000000 ) 134.6 × ( 15.9 t 2 + 67035 t 70000000 )
In the agricultural sector, demand projections assume that the percentages of each crop planted do not change, evaluating the growth in demand by the growth in area dedicated to production in millions of hectares (MHa) (Figure 5). Projections were made using data from INE (Chilean National Statistics Institute) censuses from 1976, 1997, and 2007, as well as the 2015 estimation provided by experts to DICTUC [59,79]. Space limitations, given the concentration of already productive lands within the Copiapó river valley and the lack of arable land outside the valley, would indicate a curtailing and conservative future estimate. To address the uncertainty in future agricultural markets, however, the MoS methodology incorporates all three of the different trajectories observed in the demand cone in Figure 5.
Figure 5. Agricultural growth projections with historic data.
Figure 5. Agricultural growth projections with historic data.
Water 07 05928 g005
In Copiapó, the mining sector dominates industrial production. Given the historical highs seen in copper prices during the past 15 years, mining projects in the north of Chile have grown exponentially. In 2012, the Chilean Copper Commission (COCHILCO) published maximum, optimistic, and pessimistic projections for the growth of fresh water demand in m3/s, as seen in Figure 6 [61]. Providing future foreseeable demand considering optimistic and pessimistic predictions about technological efficiency, as well as maximum extractions (m3/s) considering no technological improvements, COCHILCO clearly anticipates an exponentially growing demand function. Fitting historical data to cubic, linear, and exponential regressions, we observe that although the conservative estimate best models historical data (R2 = 0.89), it is the extreme future that meets industry expectations.
Figure 6. Mining growth projections with historic data and industry estimates.
Figure 6. Mining growth projections with historic data and industry estimates.
Water 07 05928 g006
The population of Copiapó and the surrounding towns has experienced considerable growth over the past 100 years, from 27,000 (1907) to 200,000 (2014), with increased growth observed since the early 1990s when large-scale copper mining began. According to the INE in 2010, however, this population growth will not mimic the exponential growth of the mining industry, but rather should begin to curtail as it reaches a saturation point. All three potential futures are conceivable given the past 22 years of data and should be considered for robust analysis (Figure 7). Future water demands are estimated using the three population projections after assuming a constant per capita water use during the 30-year generation simulated. The average water production by Aguas Chañar during the last 20 years has been 18.7 m3/capita/month which, considering transmission losses (40%), results in the consumption of 134.6 m3/capita/year [76].
Figure 7. Population growth projections with historic data and expert projections.
Figure 7. Population growth projections with historic data and expert projections.
Water 07 05928 g007
In addition to the possible future demand trajectories, robust risk analysis would require also considering the water supply changes expected. Using the climate model downscaling methodology previously described, we analyzed the climate predictions made using 72 different Global Climate Models under some of the SRES (Special Report on Emissions Scenarios: A1, A2, B1, B2) and RCP (Representative Concentration Pathways: 8.5, 6, 4.5, 2.6) emissions scenarios for a total of 183 different climate projections. The 30-year average percent changes in temperature and annual precipitation were calculated and plotted on a Δ T   v s .   Δ P graph (Figure 2). All models predict temperature increases, while the future trajectory of precipitation is unclear.
Eliminating models that fall more than 2.5 standard deviations above the mean in precipitation (nearly 100% increase) as outliers, we selected the four models whose extreme predictions encompass the range of climate possibilities. Csiromk360r7i1p1_rcp85 [80] foresees a hotter and drier future; ncarpcm1_b1 [81] predicts one that is hotter but also wetter; inmcm4_rcp85 [82] predicts less temperature increase but a rise in precipitation; and inmcm3_b1 [82] expects drier but slightly warmer conditions than are presently experienced in the area. We also selected the model that best represents the average change in precipitation and temperature estimated by all models, noresm1m_rcp45 [83], which forecasts a 9% decrease in precipitation and a 9% increase in temperature, both of which would further exacerbate already water-scarce conditions. Weighting factors for the simulation results using these five climate models were computed as the percentage of all 183 climate projections that are best represented by each of the five, as determined by the Pythagorean numeric distance from each (Table 6).
Table 6. Representative climate change models.
Table 6. Representative climate change models.
ScenarioModelWeight
Dry-Hotcsiromk360r7i1p1_rcp850.077
Dry-Coolinmcm3_b10.317
Wet-Hotncarpcm1_b10.093
Wet-Coolinmcm4_rcp850.093
Averagenoresm1m_rcp450.421

4. Results and Discussion

In order to demonstrate the importance of utilizing a dynamic, multivariate, robust analysis methodology in resource decision support, results from our WEAP-MODFLOW simulations were compared with those from Suarez et al. [40] who employed the static, single-variable analytics common to traditional water management. This previous research selected five general management schemes (Table 7) as an experimental manipulation of time, space, and quantity variables. While our work supports this attempt to address these variables, Suarez et al. [40] neglected the uncertain factors previously mentioned (i.e., global change, time-independent risk measures) which are necessary for robust resource planning. Utilizing the before-mentioned methodology for future planning (climate and demand) in our WEAP-MODFLOW simulation, as well as the MoS risk analysis methodology outlined in Section 2, we evaluated the same management scenarios as Suarez et al. [40] and found a nearly identical ranking of plans despite our methodological refinements (Table 7).
Table 7. Original plan results comparison between Suarez et al. [40] and MoS.
Table 7. Original plan results comparison between Suarez et al. [40] and MoS.
Water Management PlanSuarez et al., RankMoS Rank
  • Current status.
77
2.
Uniform reduction of water demand.
20%44
30%22
50%11
3.
Segmented reduction of water demand.
33
4.
Water resource management with uniform reduction of water demand and translocation of water between aquifer zones.
65
5.
Water resource management with segmented reduction of water demand and translocation of water between aquifer zones.
56
This parallel between the Suarez et al. [40] and MoS results is to be expected. These original management plans have widely varying values of demand reduction (0%–50%), thus making total water use the dominant variable. A simple analysis shows that the plan that consumes the least amount of water would be the highest rated in terms of sustainability. Indeed, the correlation between the percentage of overall demand reduction and the MoS scores is extremely high (R2 = 0.94), thus indicating a direct linear dependence of MoS upon the amount of water consumed by each plan. Therefore, to eliminate the dominating effects of demand reduction upon the MoS score in order to objectively compare the management schemes without conflicting variables, we propose four new management plan options that offer the same overall demand reduction (20%, a level achievable through water reuse technologies) while maintaining the original proportionality between water users and administrative sectors of demand reduction presented in Suarez et al. [40] defined by local stakeholders. For example, plan 5.1 (Table 8) originally called for the mining and agriculture users to reduce water demand by 20% in sectors 1 and 2, 35% in sector 3, and 50% in sectors 4 and 5. To achieve an overall demand reduction of 20%, these sectorial values were reduced proportionally by the same 4:7:10 ratio, so that in the new plan 5.1, reduction levels are 14.4% in sectors 1 and 2, 25.2% in sector 3, and 36% in sectors 4 and 5. For plan 4.1, which originally required a reduction in demand of 35% in agriculture and mining, demand reduction for these users was merely lowered to 21% so that the total reduction over all water users would be 20%. Demand transfers are not reduced as they do not affect the total demand reduction. No changes were required for plans 2.1 or 3.1 as they originally resulted in the desired 20% total demand reduction. Having eliminated demand reduction as the dominating variable in sustainability analysis, we are now able to determine the most sustainable distribution of resources between stakeholders and geographical sectors.
Following the MoS flow chart in Figure 3 to evaluate the functioning of the tested management plans outlined above, the constructed WEAP model was run 75 times considering three demand scenarios, five climate futures, and five management schemes. Utilizing the five-factor sustainability index outlined above (Equation (9)) to consider performance in the four different variable dimensions (environmental, economic, political, social), we obtained numerous results for each plan tested. These individual results were recorded in tables similar to that below (Table 9 shows results for water use plan 2.1). To obtain a unique final value for management plan comparison and ranking, the four variable dimensions and the three demand scenarios are each given equal weight, as suggested by Loucks [31]. These weights, however, are flexible and should be adjusted according to the opinions and values of the user of the MoS algorithm. Decision-makers predicting a certain demand future over another or placing more importance in one water use sector than another might weight these as they see fit. Climate, however, should remain weighted as expressed above, as these weights are derived as a general consensus of many climate models. Considering all five-factor Sustainability Index values (Equation 9) for the 60 permutations of climate (Table 6), demand (Table 5), and dimension (Table 2), and the weightings described, the final MoS is listed in the bottom right corner of the table as a single and simple number ranging from 0 to 1 (higher values represent more desirable sustainability situations). MoS scores and relative rankings for all five management schemes are presented in Table 10.
Table 8. Management schemes, original and adjusted for equal demand reduction (“sec” refers to administrative sectors, see Figure 4).
Table 8. Management schemes, original and adjusted for equal demand reduction (“sec” refers to administrative sectors, see Figure 4).
Mgmt. Plan1.12.13.14.15.1
Old plan (Suarez et al., 2014 [40])No reductionAll uses: 20% demand reductionAll uses: 20% demand reduction in sec 1 and 2, 35% in sec 3, 50% in sec 4Mining/Agriculture: 30% reduction
Potable: 50% demand transfer from 4a to 5a
Mining/Agriculture: 20% in sec 1 and 2, 35% in sec 3, 50% in sec 4 and 5
Potable: 50% demand transfer from 4a to 5a
Total demand Reduction-20%20%28%27%
New PlanNo reductionAll uses: 20% demand reductionAll uses: 20% demand reduction in sec 1 and 2, 35% in sec 3, 50% in sec 4Mining/Agriculture: 21% reduction
Potable: 50% demand transfer from 4a to 5a
Mining/Agriculture: 14.4% in sec 1 and 2, 25.2% in sec 3, 36% in sec 4 and 5
Potable: 50% demand transfer from 4a to 5a
Total demand reduction-20%20%20%20%
Table 9. Example table for the computation of MoS completed for each management option under Extreme (Ext), Average (Avg), and Conservative (Con) demand futures (using management scheme 2.1 as example).
Table 9. Example table for the computation of MoS completed for each management option under Extreme (Ext), Average (Avg), and Conservative (Con) demand futures (using management scheme 2.1 as example).
Water Use (2.1)EnvironmentalEconomicPolitical Social
Demand ModelExtAvgConExtAvgConExtAvgConExtAvgCon
Climate ScenarioDry Hot0.1980.2180.2260.0770.3440.4490.2070.4920.7320.6381.0001.000
Dry Cool0.2130.2310.2390.0800.3460.4550.2090.5060.9000.6301.0001.000
Wet Hot0.2040.2240.2310.0780.3450.4500.1990.4930.7720.6151.0001.000
Wet Cool0.1980.2300.2380.0770.3450.4530.2070.5180.8390.6301.0001.000
Average0.2080.2270.2340.0790.3450.4520.2020.4980.7800.6301.0001.000
Weighted Total0.2070.2280.2350.0790.3450.4530.2050.5020.8190.6291.0001.000
MoS0.475
Table 10. MoS Scores and rankings of adjusted management plans.
Table 10. MoS Scores and rankings of adjusted management plans.
New Mgmt Plan1.12.13.14.15.1
MoS score0.3240.4750.5470.3610.348
Ranking52134
While the purpose of the MoS was to create an index able to combine the range of supply and demand futures and their impact upon numerous variables affecting different water use dimensions into a single, easily understandable number, the analysis of the individual sustainability scores such as those in Table 9 also presents results of interest. Maintaining constant one of the two variables of future scenario uncertainty (demand and climate), the variability in sustainability introduced by the other factor may be measured for each of the four water use dimensions. Comparing the change in sustainability caused by differences in demand scenario with that produced by modification of the climate future, we observe that for all four water dimensions, demand is the more influential variable (Effect on SI of demand variance/Effect of climate variance by dimension: env., 1.89; econ., 89.5; pol., 6.9; soc., 31.4; Table 11). While much attention has been placed upon the effects of climate change on water resource availability, in arid regions where water is naturally scarce, the efficient use of this resource might perhaps be a much more important determinant of future sustainability, as indicated by these numbers.
Table 11. Effect on SI caused by future variables (second variable held constant) and relative impact.
Table 11. Effect on SI caused by future variables (second variable held constant) and relative impact.
Water Use DimensionEnvironmentalEconomicPoliticalSocial
Demand0.02720.30030.52280.3965
Climate0.01440.00340.07580.0126
Demand/Climate1.8989.496.9031.38
Higher-rated plans in the MoS index, however, are not necessarily the most practical solutions. Management plan implementation cost is an essential variable in decision-making and cannot be overlooked. A plan requiring 50% reduction in water demand, for example, would logically conserve more water, but without considering the costs of enacting such a measure we are unable to compare trade-offs. For the calculation of plan implementation costs, we assume that, without improvements in technological efficiency or philanthropy on the part of economic sectors, the reduction in fresh water demand must be replaced by other means. Following the suggestion of Suarez and current practices of the Candelaria Mine in Copiapó and other mines in arid Northern Chile, we propose SWRO (sea water reverse osmosis) as the source for replacement resources.
Currently, the cost of desalinating and transporting 0.5 m3/s to the Candelaria Mine, the principal mine in the river basin located at 675 m, is approximately $3 USD/m3. On the other hand, a lower cost of $1.78/m3 is estimated for production and transport of potable water to the city of Copiapó using SWRO [84]. Being that no study has yet to calculate the economic feasibility of utilizing desalinated water for agricultural purposes in Copiapó, we chose to utilize an average of the values for desalination already observed in Copiapó ($2.39/m3). While this value is approximately four times more expensive than the highly efficient Ashkalon desalination plant in Israel ($0.57/m3), it falls within international norms [85,86]. Based on projected demands for each stakeholder under the three different future demand scenarios and the associated demand reduction percentages proposed by each management plan, we calculated a first-order approximation of plan implementation costs under the above-stated desalination costs (Table 12, column 2).
Table 12. Cost and performance evaluations of new management plans. Costs are calculated using the average cost of all demand and climate scenarios.
Table 12. Cost and performance evaluations of new management plans. Costs are calculated using the average cost of all demand and climate scenarios.
PlanImplementation Cost (Million $$)Deficit Reduction Cost (Million $$)Total Cost (Million $$)MoS
1.1015,95315,9530.324
2.16901771814,6190.475
3.16865711513,9800.547
4.1(new)667114,93521,6060.361
5.1(new)6632770514,3370.348
The highest-performing management plan utilizing the traditional single-variable static approach seen in Suarez et al. (50% reduction in all demands) would also be the most expensive option to implement, costing an estimated $9.3 billion USD more than the next most expensive option [40]. Suarez et al. [40] agree upon the infeasibility of this option due to difficulty of implementation and instead propose plan 2.2 (30% reduction in all demands), which still happens to be the second most expensive option, costing an estimated $2.8 billion USD more than the next most expensive option. These additional costs over a more traditional method that employs static future demand scenarios result from the fact that even in the best case scenario, conservative demand projections are greater than current consumption.
Implementation, however, is not the only cost inherent in a management plan, but as Suarez et al. [40] correctly argue, each plan will result in a different future water deficit based upon the ability to meet demand after implementation. While Suarez’s cost calculations were a necessary first attempt to convert hydrological performance into an economic value, certain assumptions led to a large underestimation of deficit cost. The potential effects associated with assuming static futures in climate and markets have been addressed in this work, but the assumption that ex post facto deficit recovery involves the same cost as preventative resource production must also be addressed. While the nearly deterministic nature of their demands allows the mining and potable water sectors to prepare for deficits, the fact that additional water demand and deficits are calculated after the fact for agricultural purposes dependent upon stochastic rainfall limits the plausibility of desalination for deficit reduction (at the rate of $1.35–$2.01/m3 as employed by Suarez). Rather, the most likely costs to be incurred by farmers will come in the form of reduced crop yield. Taking the crop irrigation demand requirements for Copiapó [59] of 17,735 m3/ha/year, the average yield for table grape production in climatic zones similar to Copiapó [87] of 19,952 metric tons/ hectare, and the average table grape prices over the last five years for produce originating in Copiapó of $1.1/kg as reported by the Chilean Office of Agrarian Studies and Politics ODEPA [88], we find that average economic product lost by deficit is $14.85/m3. Clearly, when compared with the marginal cost of SWRO (~$2.39/m3), economic loss due to water deficit is much greater and should be avoided. Were it not possible to reutilize the desalination plants mentioned for management plan implementation, increasing their capacity so as to meet deficit needs, this new factor, combined with the dominance of deficits by the agricultural sector, would greatly raise the costs of deficit reduction and increase the importance of implementing preventative sustainability plans.
Independent of this decision to calculate crop losses or agricultural drought prevention, the accepted static futures methodology for resource analysis (e.g., [40]) underestimates deficit by a factor of more than 10 times by overlooking the projections of mainly increased demand and, to a minor extent, decreased supply. It should also be noted that although Suarez et al. [40] estimated that a 50% demand reduction would lead to aquifer recharge and thus an environmentally sustainable position, by excluding the potential effects of climate change and demand growth, they likely overstate their case. According to calculations in this study that consider these factors, even a 50% reduction in future demand, at a cost of at least $16.1 billion USD for implementation, would not be sufficient to reach absolute sustainability in Copiapó, and would still require additional expenditure for deficit reduction through SWRO. As might be seen in Table 12 below, at that cost it would be more economic to follow nearly any other management scheme, thus the consideration of future demand growth is integral to the decision-making process.
Combing the costs of implementation of the new management plans (equal percentage of demand reduction, Table 8) with the costs of deficits caused by the failure of these plans to meet demand requirements, we arrive at a first order estimation for the total cost of adopting the management plans for water resource sustainability. As might be seen in Table 12 below, plan 3.1 (previously ranked only third by Suarez et al. [40]) now not only stands as the best-performing option according to our MoS, but also is estimated to be the least expensive.
Merely by showing that results are different does not necessarily prove that the MoS methodology is any more appropriate than the traditional method which does not account for climate change, demand change, or multi-sectorial multivariable analysis. Being that these techniques are intended to aid in resource planning for uncertain futures, their sensitivity to variations in future variables is of prime concern. A more robust analysis methodology would be one that either eliminates variable uncertainty or includes the range of uncertainty into calculations. Maintaining constant the climate scenario (dry-cool, dry-hot, wet-cool, wet-hot, average), we calculate the range of total water deficit during the 30-year period over all three demand scenarios (conservative, average, extreme) (Table 13). By eliminating the effects of the climate variable, we observe that the range of uncertainty introduced by the demand variable is much greater (between 3200 and 5600 Mm3 depending on the management scheme) than the range between all management schemes using the traditional methodology (500 Mm3) observed in Suarez et al. [40]. The uncertainty caused by climate also results in a much greater range under the extreme demand scenario (1200–2000 Mm3) and also under conservative and average demand if we exclude the base scheme in which no water conservation measures are taken. In fact, by overlooking climate and demand changes altogether, the traditional method underestimates water deficits even in the best-case scenario (conservative demand and cool-wet climate future) in all but the base scenario by an average of 250%. These ranges of uncertainty introduced by variables that are overlooked by the traditional analysis technique observed in Suarez et al. [40], being greater than the range of deficit values between the potential management schemes being analyzed (experimental variable), demonstrate the need for a robust analytical methodology. Reliant upon uncontrolled stochastic variables rather than upon robust analysis makes the conclusions of the traditional method highly fragile for decision-making under an uncertain future. The MoS methodology, however, includes these uncertainties in its calculations by estimating the range of effects caused by potential futures rather than gambling on a single demand and climate scenario. Thus, the high range of potential water deficits resultant from uncertain futures, rather than introducing a conflicting variable, is incorporated into the MoS evaluation of management schemes.
Table 13. Range of deficits introduced by uncertain variables in million m3, considering either a single demand scenario over all climates or a single climate scenario over all demands. To the right are projected deficits using the traditional and MoS techniques.
Table 13. Range of deficits introduced by uncertain variables in million m3, considering either a single demand scenario over all climates or a single climate scenario over all demands. To the right are projected deficits using the traditional and MoS techniques.
Mgmt SchemeSingle Demand Scenario (All Climates Considered)Single Climate Scenario (All Demands Considered)All Variables Considered
Cons.Avg.Ext.Avg.Dry-CoolDry-HotWet-CoolWet-HotSuarez et al.MoS
1.11392671972539752755588530055015282311
2.152120140839823887412839134052321534
3.12791136439363814412238414038181438
4.152789140939823887412239148686541795
5.162128142839993925411839464051451571
While the entirety of the difference in the MoS score does not lie in a single category but rather is distributed throughout the different sustainability indicators and water use dimensions, a gain in 0.072 MoS as that observed between plans 3.1 and 2.1 is equivalent to 25.8 more months where water is provided to all water utility clients during the 30-year period, or it is equivalent to the average aquifer depth being 5.59 m closer to the surface. This 25.8-month difference in reliability multiplied by 239 clients without water on average per failure and 16.2 m3/client average monthly usage reported by Aguas Chañar, who charges approximately $1.31/m3, gives a total of $130,868 USD gained by choosing the better-performing option. The 5.59 m difference in average aquifer depth translates into a difference of $15.1 million USD extra expenditure in pumping costs during the 30-year period when choosing plan 2.1 over 3.1. According to these numbers it may be observed that the MoS does not have a strict monetary interpretation as the same reduction of 0.072 in the index score would correlate to incurred costs that vary over two orders of magnitude. This difference is to be expected, however, as potable water usage is nearly two orders of magnitude lower than total aquifer pumping. These translations of the MoS into monetary values are merely examples to aid in user understanding of the index. Nevertheless, decision-makers should be presented with both plans and costs to aid in decision-making.

5. Conclusions

Addressing issues of oversimplification of futures observed in current sustainability analysis methods, this research proposes a robust, dynamic, and multivariate approach for decision support in water resources. While the simulation of results in multiple climate and demand scenarios increases calculation cost, by evaluating the range of watershed performance in response to these uncertain parameters avoids potentially costly underestimations produced by simpler techniques. To demonstrate its application, this new methodology (the Measure of Sustainability or “MoS”) was applied to the Copiapó, Chile, watershed, an example of an arid watershed under extreme natural and anthropogenic water stress. Results of the MoS were compared to those previously calculated by a more traditional approach in which climate and water demands are considered to be stationary.
The MoS results, coupled with first order estimated costs, provide decision-makers with a simplified manner of selecting resource management possibilities whose economic feasibility might further be analyzed by financial experts. In the analysis of Copiapó we observed that the sustainability of the original plans evaluated by a more traditional approach [40] is highly dependent upon the water reduction and thus the inherent implementation cost proposed by each. Even in the best possible scenario of a conservative demand future, the traditional method of static demand underestimates the cost to implement any management scheme while overlooking the effects of climate change also misjudges the magnitude of future water deficits. Specifically in the agricultural sector, the marginal cost of crop loss due to water deficit is much greater than the marginal cost of SWRO and thus should be avoided where possible through preparation rather than remediation. This leads us to suggest further investigation into the economic feasibility of SWRO solutions for agricultural deficit mitigation.
Independently analyzing the uncertain future variables of climate and demand, we observe that the uncertainty introduced by each of these is larger than the range of the experimental variable deficit of each management scheme in the traditional methodology. In addition, we observe that in arid regions, such as Copiapó, Chile, where water resources are naturally scarce, demand variability might introduce more risk to sustainability than potential changes in climate, leading us to suggest increased efforts in hydrologic efficiency. Combining these two variables, however, to compare four management schemes with equal water reduction but different distribution structures, we were able to find that reducing all demands by a strict 20% (plan 2.1) and reducing all demands by geographic sector (plan 3.1) are the best performing management schemes, with 3.1 scoring slightly higher in the MoS index (0.072) but also costing slightly less (approximately $639 million USD). Not only for lower economic cost but also for improved performance in the Measure of Sustainability (MoS) index, we propose the adoption of plan 3.1 over Suarez et al.’s selection of 2.1 for water management in the Copiapó basin.
Mounting concerns about the impacts of human activities, potential climatic shifts, expanding populations and water demands as well as new knowledge are part of the pressing need to develop alternative management schemes in an integrated manner, particularly in environments with scarce natural resources. In this background, it is important to emphasize the transition in theory and practice from Environmental Impact Assessment to Risk Assessment, and then to Sustainability Analysis and Vulnerability Assessment. Thus, given the uncertainty of future events and the dependence upon climate and market trajectories, we reiterate that the plans selected by this methodology are not necessarily the most optimal solutions to problems in water resource sustainability. By considering the performance of each plan under the expected range of likely futures, however, the MoS methodology proposes solutions that are robust to variations in future parameter values and are thus the best management options in a stochastic natural world. While important, sustainability and cost are not the only factors that influence environmental decision-making. Social acceptance, political feasibility, and technological capabilities also influence the opinions of stakeholders and thus add a component of uncertainty in values not directly addressed by the MoS methodology. Further research in decision support should address issues of political and social will for integrated conservation measures. In particular, due to its dominance of demand, inelasticity of infrastructure, and disjointed individualist organization, the agricultural sector’s motivation to participate in water demand management should be of pressing interest for future study. As observed by considering the range of future water supply and demand trajectories, the MoS demonstrates that sustainability management is even more vital than previously thought.

Acknowledgments

The authors acknowledge the funding from the European Union (EU) Seventh Framework Programme (FP7/2007-2013) under grant agreement NO. 283025 to Project COROADO, the Centro de Desarrollo Urbano Sustentable (CEDEUS, Conicyt/Fondap/15110020) and the Centro de Investigación para la Gestión Integrada de Desastres Naturales (CIGIDEN Conicyt/Fondap/15110017). This work was also funded by CONICYT Programa de Formación de Capital Humano Avanzado/Doctorado Nacional 21130266. We also acknowledge the Chilean Ministry of Public Works and Alvaro Sola for providing relevant information for the study, Jean-Marc Dorsaz for his work on early stages of this paper, and Centro de Cambio Global UC for providing downscaled climate change data. Finally, we also thank two referees for their detailed comments and constructive suggestions. The present work was partially developed within the framework of the Panta Rhei Research Initiative of the International Association of Hydrological Sciences.

Author Contributions

The manuscript was primarily written by Christian K. Hunter, with all authors contributing to its preparation and review. The WEAP-MODFLOW model was developed by Christian K. Hunter under the guidance of Jorge Gironás and Diogo Bolster. Data analysis was carried out by Christian K. Hunter, discussing of results together with Jorge Gironás and Diogo Bolster.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cosgrove, W.J.; Rijsberman, F.R. World Water Vision: Making Water Everybody’s Business; Earthscan Publications: London, UK, 2000. [Google Scholar]
  2. Rosegrant, M.W.; Cai, X.; Cline, S.A. Averting an Impending Crisis, Global Water Outlook to 2025, Food Policy Report; International Water Management Institute (IWMI): Colombo, Sri Lanka, 2002. [Google Scholar]
  3. Rijsberman, F.R. Water scarcity: Fact or fiction? Agric. Water Manag. 2006, 80, 5–22. [Google Scholar] [CrossRef]
  4. Karavitis, C.A. Towards Strategic Water Resources Management Policies: Challenges and Options. In Proceedings of the HYDRORAMA 2002, 3rd International Forum, Athens, Greece, 21–22 March 2002; pp. 218–224.
  5. Katsiardi, P.; Manoli, E.; Karavitis, C.A.; Assimacopoulos, D. Scenario-based strategy development for integrated water management. Glob. NEST J. 2005, 7, 360–368. [Google Scholar]
  6. Maia, R.; Silva, C. DSS application at a river basin scale, taking into account water resources exploitation risks and associated costs: The Ribeiras do Algarve River Basin. Desalination 2008, 237, 81–91. [Google Scholar] [CrossRef]
  7. Cowling, R.M.; Egoh, B.; Knight, A.T.; O’Farrell, P.J.; Reyers, B.; Rouget, M.; Roux, D.J.; Welz, A.; Wilhelm-Rechman, A. An operational model for mainstreaming ecosystem services for implementation. Proc. Natl. Acad. Sci. USA 2008, 105, 9483–9488. [Google Scholar] [CrossRef] [PubMed]
  8. Ranganathan, J.; Raudsepp-Hearne, C.; Lucas, N.; Irwin, F.; Zurek, M.; Bennett, K.; Ash, N.; West, P. Ecosystem Services: A Guide for Decision-Makers; World Resources Institute: Washington, DC, USA, 2008. [Google Scholar]
  9. Daily, G.C.; Polasky, S.; Goldstein, J. Ecosystem services in decision making: Time to deliver. Front. Ecol. Environ. 2009, 7, 21–28. [Google Scholar] [CrossRef]
  10. Asano, T.; Levine, A.D. Wastewater reclamation recycling and reuse: Past, present, and future. Water Sci. Technol. 1996, 33, 1–4. [Google Scholar] [CrossRef]
  11. Angelakis, A.N.; Koutsoyiannis, D. Urban water resources management in ancient Greek times. In The Encyclopedia of Water Sciences; Stewart, B.A., Howell, T., Eds.; Markel Dekker Inc.: New York, NY, USA, 2003; pp. 999–1008. [Google Scholar]
  12. Angelakis, A.N.; Durham, B. Water recycling and reuse in EUREAU countries: Trends and challenges. Desalination 2006, 218, 3–12. [Google Scholar] [CrossRef]
  13. Lettenmaier, D.P.; Wood, A.W.; Palmer, R.N.; Wood, E.F.; Stakhiv, E.Z. Water resources implications of global warming: A US regional perspective. Clim. Chang. 1999, 43, 537–579. [Google Scholar] [CrossRef]
  14. Anderson, J. The environmental benefits of water recycling and reuse. Water Supply 2003, 3, 1–10. [Google Scholar]
  15. Hochstrat, R.; Wintgens, T.; Melin, T.; Jeffrey, P. Assessing the European wastewater reclamation and reuse potential—A scenario analysis. Desalination 2006, 188, 1–8. [Google Scholar] [CrossRef] [Green Version]
  16. Wintgens, T.; Salehi, F.; Hochstrat, R.; Melin, T. Emerging contaminants and treatment options in water recycling for indirect potable use. Water Sci. Technol. 2008, 57, 99–107. [Google Scholar] [CrossRef] [PubMed]
  17. Milly, P.C.D.; Betancourt, J.; Falkenmark, M.; Hirsch, R.M.; Zbigniew, W.; Lettenmaier, D.P.; Stouffer, R.J. Stationarity is dead: Whither water management? Science 2008, 319, 573–574. [Google Scholar] [CrossRef] [PubMed]
  18. Pachauri, R.K.; Allen, M.R.; Barros, V.R.; Broome, J.; Cramer, W.; Christ, R.; Church, J.A.; Clarke, L.; Dahe, Q.; Dasgupta, P.; et al. Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Intergovernmental Panel on Climate Change (IPCC): Geneva, Switzerland, 2014. [Google Scholar]
  19. Johnston, P.; Everard, M.; Santillo, D.; Robért, K. Reclaiming the definition of sustainability. Environ. Sci. Pollut. Res. Int. 2007, 14, 60–66. [Google Scholar] [PubMed]
  20. Ludwig, D. Magic, science, and religion in natural resource management. Ecol. Appl. 1993, 3, 555–558. [Google Scholar]
  21. World Commission on Environment and Development (WCED). Our Common Future (The Brundtland. Report); Oxford University Press: Oxford, UK, 1987. [Google Scholar]
  22. Constanza, R.; Patten, B.C. Defining and predicting sustainability. Ecol. Econ. 1995, 15, 193–196. [Google Scholar] [CrossRef]
  23. Brown, B.J.; Hanson, M.E.; Liverman, D.M.; Merideth, R.W., Jr. Global sustainability: Toward definition. Environ. Manag. 1987, 11, 713–719. [Google Scholar] [CrossRef]
  24. Takeuchi, K.; Hamlin, M.; Kundzewicz, Z.W.; Rosbejerg, D.; Simonovic, S.P. Sustainable Reservoir Development and Management; No. 251; International Association of Hydrological Sciences: Wallingford, UK, 1998. [Google Scholar]
  25. Gatto, M. Sustainability: Is it a well-defined concept? Ecol. Appl. 1995, 5, 1181–1183. [Google Scholar]
  26. Regan, H.M.; Ban-Haim, Y.; Langford, B.; Wilson, W.G.; Lundberg, P.; Andelman, S.J.; Burgman, M.A. Robust decision-making under severe uncertainty for conservation management. Ecol. Appl. 2005, 15, 1471–1477. [Google Scholar] [CrossRef]
  27. McDonald-Madden, E.; Baxter, P.W.J.; Possingham, H.P. Making robust decisions for conversation with restricted money and knowledge. J. Appl. Ecol. 2008, 45, 1630–1638. [Google Scholar] [CrossRef]
  28. Sainz, P. An Index of Social Welfare. Towards A New Way to Measure Development; Office of the South Commission: Caracas, Venezuela, 1989. [Google Scholar]
  29. Hashimoto, T.; Stedinger, J.R.; Loucks, D.P. Reliability, resiliency and vulnerability criteria for water resources planning. Water Resour. Res. 1982, 18, 14–21. [Google Scholar] [CrossRef]
  30. Hipel, K.W. Multiple objective decision making in water resources. J. Am. Water Resour. Assoc. 1992, 28, 3–12. [Google Scholar] [CrossRef]
  31. Loucks, D.P. Quantifying trends in system sustainability. Hydrol. Sci. J. 1997, 42, 513–530. [Google Scholar] [CrossRef]
  32. Tkach, R.J.; Simonovic, P. A new approach to multi-criteria decision making in water resrouces. J. Geogr. Inf. Decis. Anal. 1997, 1, 25–44. [Google Scholar]
  33. Howe, C.; White, S. Integrated resource planning for water and wastewater. Water Int. 1999, 24, 356–362. [Google Scholar] [CrossRef]
  34. Loomis, J. Environmental valuation techniques in water resource decision making. J. Water Resour. Plan. Manag. 2000, 126, 339–344. [Google Scholar] [CrossRef]
  35. Sandoval-Solis, S.; McKinney, D.C.; Loucks, D.P. Sustainability index for water resources planning and management. J. Water Resour. Plan. Manag. 2011, 137, 381–390. [Google Scholar] [CrossRef]
  36. Zongxue, X.; Jinno, K.; Kawanura, A.; Takesaki, S.; Ito, K. Performance risk analysis for Fukuoka water supply system. Water Resour. Manag. 1998, 12, 13–30. [Google Scholar] [CrossRef]
  37. Rodak, C.; Silliman, S.E.; Bolster, D. Time-dependent health risk from contaminated groundwater including use of reliability, resilience, and vulnerability as measures. J. Am. Water Resour. Assoc. 2014, 50, 14–28. [Google Scholar] [CrossRef]
  38. Carlson, J.M.; Doyle, J. Complexity and robustness. Proc. Natl. Acad. Sci. USA 2002, 99, 2538–2545. [Google Scholar] [CrossRef] [PubMed]
  39. Ascough, J.C.; Maier, H.R.; Ravalico, J.K.; Strudley, M.W. Future research challenges for incorporation of uncertainty in environmental and ecological decision-making. Ecol. Model. 2008, 219, 383–399. [Google Scholar] [CrossRef]
  40. Suárez, F.; Muñoz, J.F.; Fernández, B.; Dorsaz, J.M.; Hunter, C.K.; Karavitis, C.A.; Gironás, J. Integrated Water Resource Management and Energy Requirements for Water Supply in the Copiapó River Basin, Chile. Water 2014, 6, 2590–2613. [Google Scholar] [CrossRef]
  41. Shelef, G.; Azov, Y. The coming ERA of intensive wastewater reuse in the Mediterranean region. Water Sci. Technol. 1996, 33, 115–125. [Google Scholar] [CrossRef]
  42. Cai, X.; Lasdon, L.; Michelsen, A.M. Group decision making in water resources planning using multiple objective analysis. J. Water Resour. Plan. Manag. 2004, 130, 4–14. [Google Scholar] [CrossRef]
  43. Srdjevic, B.; Medeiros, Y.D.P.; Faria, A.S. An objective multi-criteria evaluation of water management scenarios. Water Resour. Manag. 2004, 18, 35–54. [Google Scholar] [CrossRef]
  44. Hadded, R.; Nouiri, I.; Alshihabi, O.; Mabmann, J.; Huber, M.; Laghouane, A.; Yahiaoui, H.; Tarhouni, J. A decision support system to manage the groundwater of the Zeuss Koutine aquifer using the WEAP-MODFLOW framework. Water Resour. Manag. 2013, 27, 1981–2000. [Google Scholar] [CrossRef]
  45. Loucks, D.P. Sustainable Water Resources Management. Water Int. 2000, 25, 3–10. [Google Scholar] [CrossRef]
  46. Kay, P.A. Measuring sustainability in Israel’s water system. Water Int. 2000, 25, 617–623. [Google Scholar] [CrossRef]
  47. Kjeldsen, T.R.; Rosbjerg, D. A framework for assessing water resources system sustainability. In Regional Management of Water Resources; Schumann, A.H., Acreman, M.C., Davis, R., Marino, M.A., Rosbjerg, D., Jun, X., Eds.; IAHS Press: Wallingford, UK, 2001; pp. 107–113. [Google Scholar]
  48. Notestein, F.W. Population: The Long View. In Food for the World; Schultze, T.W., Ed.; University of Chicago Press: Chicago, IL, USA, 1945. [Google Scholar]
  49. Lee, C. Models in Planning: An Introduction to the Use of Quantitative Models in Planning. Elsevier: Amsterdam, The Netherlands, 2013. [Google Scholar]
  50. Judge, G.G.; Griffiths, W.E.; Hill, R.C.; Lee, T.C. The Theory and Practice of Econometrics; Wiley: New York, NY, USA, 1980. [Google Scholar]
  51. Brekke, L.; White, D.; Olsen, K., Jr.; Townsley, E.; Williams, D. Addressing Climate Change in Long-Term Water Resources Planning and Management: User Needs for Improving Tools and Informatio; Government Reports Announcements & Index (GRA&I); Bureau of Reclamation: Chester, MT, USA, 2011; p. 11. [Google Scholar]
  52. Wood, A.W.; Leung, L.R.; Sridhar, V.; Lettenmaier, D.P. Hydrologic implications of dynamical and statistical approaches to downscaling climate model outputs. Clim. Chang. 2004, 62, 189–216. [Google Scholar] [CrossRef]
  53. Meza, F.J.; Vicuña, S.; Jelinek, M.; Bustos, E.; Bonelli, S. Assessing water demands and coverage sensitivity to climate change in the urban and rural sectors in central Chile. J. Water Clim. Chang. 2014, 5, 192–203. [Google Scholar] [CrossRef]
  54. Shepard, D. A two-dimensional interpolation function for irregularly-spaced data. In Proceedings of the 1968 23rd ACM National Conference, New York, NY, USA, 27–29 August 1968; pp. 517–524.
  55. Bonelli, S.; Vicuña, S.; Meza, F.J.; Gironás, J.; Barton, J. Incorporating climate change adaptation strategies in urban water supply planning: The case of central Chile. J. Water Clim. Chang. 2014, 5, 357–376. [Google Scholar] [CrossRef]
  56. Lijphart, A. Electoral Systems and Party Systems; Oxford University Press: Oxford, UK, 1994. [Google Scholar]
  57. Consejo Minero News Page, Construcción Planta Desalinizadora y Otorgamiento del uso de Derechos de Agua en la Región de Atacama. Available online: http://www.consejominero.cl/buena-practica-de-la-gran-mineria-freeport-mcmoran-candelaria-construccion-planta-desalinizadora-y-otorgamiento-del-uso-de-derechos-de-agua-en-la-region-de-atacama/ (accessed on 17 August 2015).
  58. McFarlane, D.; Norgate, T. Summary Report on Copiapó. Water Yields and Demands; Technical Report for AusAID as part of the Minerals Down Under Flagship; CSIRO: Urrbrae, Australia, 2012. [Google Scholar]
  59. División de Ingeniería Hidráulica y Ambiental, Pontificia Universidad Católica de Chile (DICTUC). Análisis Integrado de Gestión en Cuenca del río Copiapó; DICTUC: Santiago, Chile, 2010. [Google Scholar]
  60. Yates, D.; Sieber, J.; Purkey, D.R.; Huber-Lee, A. WEAP21—A Demand-, Priority-, and Preference-Driven Water Planning Model: Part 1, Model Characteristics. Water Int. 2005, 30, 487–500. [Google Scholar] [CrossRef]
  61. COCHILCO, Comisión Chilena del Cobre. Yearbook: Copper and Other Mineral. Statistics 1987–2006. Available online: http://www.cochilco.cl/desarrollo/estudios/anuario-pdf-final.pdf (accessed on 4 August 2007).
  62. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration-Guidelines for Computing Crop Water Requirements-FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998; p. 300. [Google Scholar]
  63. Dirección General de Aguas (DGA). Levantamiento Catastro de Usuarios de Aguas del Valle del Río Copiapó y sus Afluentes. III Región de Atacama. S.I.T. Nº150; Departamento de Estudios y Planificación, Dirección General de Aguas, Ministerio de Obras Públicas: Santiago, Chile, 2008. [Google Scholar]
  64. Thornthwaite, C.W. An approach toward a rational classification of climate. Geogr. Rev. 1948, 38, 55–94. [Google Scholar] [CrossRef]
  65. Thornthwaite, C.W.; Mather, J.R. The Water Balance. Publ. Climatol. 1955, 8, 1–104. [Google Scholar]
  66. Thornthwaite, C.W.; Mather, J.R. Instructions and Tables for the Computing Potential Evapotranspiration and the Water Balance. Publ. Climatol. 1957, 10, 311. [Google Scholar]
  67. Black, P.E. Revisiting the Thornthwaite and Mather water balance. J. Am. Water Resour. Assoc. 2007, 43, 1604–1605. [Google Scholar] [CrossRef]
  68. Pruitt, W.O. Relation of consumptive use of water to climate. Am. Soc. Agric. Eng. Trans. 1960, 3, 9–13. [Google Scholar] [CrossRef]
  69. Pelton, W.L.; King, K.M.; Tanner, C.B. An evaluation of the Thornthwaite and mean temperature methods for determining potential evapotranspiration. Agron. J. 1960, 52, 387–395. [Google Scholar] [CrossRef]
  70. Stanhill, G. A comparison of methods of calculating evapotranspiration from climatic data. Isr. J. Agric. Res. 1961, 11, 159–171. [Google Scholar]
  71. Hashemi, F.; Habibian, M.T. Limitations of temperature-based methods in estimating crop evapotranspiration in arid-zone agricultural projects. Agric. Meteorol. 1979, 20, 237–247. [Google Scholar] [CrossRef]
  72. Pereira, A.R.; de Camargo, Â.P. An analysis of the criticism of Thornthwaite’s equation for estimating potential evapotranspiration. Agric. For. Meteorol. 1989, 46, 149–157. [Google Scholar] [CrossRef]
  73. Garcia, M.; Raes, D.; Allen, R.; Herbas, C. Dynamic of reference evapotranspiration in the Bolivian highlands (Altiplano). Agric. For. Meteorol. 2004, 125, 67–82. [Google Scholar] [CrossRef]
  74. Bookman and Edmonston Inc. Report on Investigation of the Water Resources of Kaweah Delta Water Conservation District; Technical report for the Kaweah Delta Water Conservation District: Visalia, CA, USA, 1972. [Google Scholar]
  75. Fugro West Inc. Water Resources Investigation of the Kaweah Delta Water Conservation District Interim Report Task 4—Surface Water; Technical report for the Kaweah Delta Water Conservation District: Visalia, CA, USA, 2003. [Google Scholar]
  76. Superintendencia de Servicios Sanitarios (SISS). Informe de Gestión del Sector Sanitario. Superintendencia de Servicios Sanitarios; Ministerio de Obras Públicas: Santiago, Chile, 2014. [Google Scholar]
  77. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual models part I—A discussion of principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  78. Legates, D.R.; McCabe, G.J., Jr. Evaluating the use of “goodness-of-fit” measures in hydrologic and hydroclimatic model validation. Water Resour. Res. 1999, 35, 233–241. [Google Scholar] [CrossRef]
  79. Instituto Nacional de Estadisticas (INE). Cambios Estructurales en la Agricultura Chilena; Instituto Nacional de Estadisticas, Ministerio de Economia: Santiago, Chile, 2009. [Google Scholar]
  80. Collier, M.A.; Jeffrey, S.J.; Rotstayn, L.D.; Wong, K.K.; Dravitzki, S.M.; Moseneder, C.; Hamalainen, C.; Syktus, J.I.; Suppiah, R.; Antony, J.; et al. The CSIRO-Mk3. 6.0 Atmosphere-Ocean. GCM: Participation in CMIP5 and data publication. In Proceedings of the International Congress on Modelling and Simulation–MODSIM, Perth, Australia, 12–16 December 2011.
  81. Washington, W.M.; Weatherly, J.W.; Meehl, G.A.; Semtner, A.J., Jr.; Bettge, T.W.; Craig, A.P.; Strand, W.G., Jr.; Arblaster, J.; Wayland, V.B.; James, R.; et al. Parallel climate model (PCM) control and transient simulations. Clim. Dyn. 2000, 16, 755–774. [Google Scholar] [CrossRef]
  82. Volodin, E.M. Atmosphere-ocean general circulation model with the carbon cycle. Izv. Atmos. Ocean. Phys. 2007, 43, 266–280. [Google Scholar] [CrossRef]
  83. Bentsen, M.; Bethke, I.; Debernard, J.B.; Iversen, T.; Kirkevåg, A.; Seland, Ø.; Drange, H.; Roelandt, C.; Seierstad, I.A.; Hoose, C.; et al. The Norwegian Earth System Model, NorESM1-M—Part 1: Description and basic evaluation of the physical climate. Geosci. Model Dev. 2013, 6, 687–720. [Google Scholar] [CrossRef]
  84. Aqua Advise. Estudio de Plantas Desaladoras Para Consumo Humano Para la Provincia de Copiapó y Comuna de Chañaral.; Technical report for CCIRA: Copiapó, Chile, 2012. [Google Scholar]
  85. Dreizin, Y. Ashkelon seawater desalination project—Off-taker’s self-costs, supplied water costs, total costs and benefits. Desalination 2006, 190, 104–116. [Google Scholar] [CrossRef]
  86. WateReuse. Seawater Desalination Costs: White Paper. Available online: http://www.watereuse.org/sites/default/files/u8/WateReuse_Desal_Cost_White_Paper.pdf (accessed on 6 June 2015).
  87. University of California Cooperative Extension. Sample Costs to Establish and Produce Table Grapes: San Joaquin Vally South; University of California: Davis, CA, USA, 2007; Available online: http://coststudyfiles.ucdavis.edu/uploads/cs_public/79/08/79084aba-1425–430c-b3b84482182b6c07/grapecrimsonvs2007.pdf (accessed on 17 March 2015).
  88. ODEPA. Boletín Diario de Frutas y Hortalizas en Terminales; Oficina de Estudios y Políticas Agrarias, Ministerio de Agricultura: Santiago, Chile, 2015; Available online: http://www.ine.cl/canales/chile_estadistico/estadisticas_agropecuarias/2009/07/cambios_estructurales.pdf (accessed on 31 March 2015).

Share and Cite

MDPI and ACS Style

Hunter, C.; Gironás, J.; Bolster, D.; Karavitis, C.A. A Dynamic, Multivariate Sustainability Measure for Robust Analysis of Water Management under Climate and Demand Uncertainty in an Arid Environment. Water 2015, 7, 5928-5958. https://doi.org/10.3390/w7115928

AMA Style

Hunter C, Gironás J, Bolster D, Karavitis CA. A Dynamic, Multivariate Sustainability Measure for Robust Analysis of Water Management under Climate and Demand Uncertainty in an Arid Environment. Water. 2015; 7(11):5928-5958. https://doi.org/10.3390/w7115928

Chicago/Turabian Style

Hunter, Christian, Jorge Gironás, Diogo Bolster, and Christos A. Karavitis. 2015. "A Dynamic, Multivariate Sustainability Measure for Robust Analysis of Water Management under Climate and Demand Uncertainty in an Arid Environment" Water 7, no. 11: 5928-5958. https://doi.org/10.3390/w7115928

APA Style

Hunter, C., Gironás, J., Bolster, D., & Karavitis, C. A. (2015). A Dynamic, Multivariate Sustainability Measure for Robust Analysis of Water Management under Climate and Demand Uncertainty in an Arid Environment. Water, 7(11), 5928-5958. https://doi.org/10.3390/w7115928

Article Metrics

Back to TopTop