Next Article in Journal
Impacts of Legal and Institutional Changes on Irrigation Management Performance: A Case of the Gezira Irrigation Scheme, Sudan
Next Article in Special Issue
Landscape Drivers and Social Dynamics Shaping Microbial Contamination Risk in Three Maya Communities in Southern Belize, Central America
Previous Article in Journal
Analysis of the Public Flood Risk Perception in a Flood-Prone City: The Case of Jingdezhen City in China
Previous Article in Special Issue
Modeling Water Yield: Assessing the Role of Site and Region-Specific Attributes in Determining Model Performance of the InVEST Seasonal Water Yield Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Socioeconomic and Environmental Proxies for Comparing Freshwater Ecosystem Service Threats across International Sites: A Diagnostic Approach

by
Thomas C. Harmon
1,*,
Robyn L. Smyth
2,
Sudeep Chandra
3,
Daniel Conde
4,18,
Ramesh Dhungel
5,
Jaime Escobar
6,7,
Natalia Hoyos
8,
Juan Pablo Lozoya
4,
Mariana Nin
4,
Gerardo M.E. Perillo
9,
Stephanie Pincetl
10,
M. Cintia Piccolo
9,
Brian Reid
11,
James A. Rusak
12,13,
Facundo Scordo
9,
Maria I. Velez
14,
Sandra R. Villamizar
15,
Beverley Wemple
16 and
Mariana Zilio
17
1
Environmental Systems Graduate Program, Department of Civil & Environmental Engineering, University of California Merced, Merced, CA 94343, USA
2
Environmental and Urban Studies, Bard College, Annandale-on-Hudson, New York, NY 12504, USA
3
Department of Biology, University of Nevada, Reno, NV 89557, USA
4
Centro de Manejo Costero (CURE), Universidad de la República, Maldonado 20000, Uruguay
5
Department of Agronomy, Kansas State University, Manhattan, KS 66506, USA
6
Department of Civil and Environmental Engineering, Universidad del Norte, Barranquilla 080001, Colombia
7
Center for Tropical Paleoecology and Archaeology, Smithsonian Tropical Research Institute, Panama City 32401, Panama
8
Department of History and Social Sciences, Universidad del Norte, Barranquilla 080001, Colombia
9
Instituto Argentino de Oceanografía, CONICET & Universidad Nacional del Sur, Bahía Blanca 8000, Argentina
10
Institute of the Environment, University of California Los Angeles, Los Angeles, CA 90095, USA
11
Centro de Investigaciones en Ecosistemas de la Patagonia, Universidad Austral de Chile, Coyhaique 5951601, Chile
12
Department of Biology, Queen’s University, Kingston, ON K7L 3N6, Canada
13
Dorset Environmental Science Centre, Dorset, ON P0A 1E0, Canada
14
Department of Geology, University of Regina, Regina, SK S4S 4P5, Canada
15
Department of Civil Engineering, Universidad Pontificia Bolivariana, Bucaramanga 681007, Colombia
16
Department of Geography, University of Vermont, Burlington, VT 05405, USA
17
Instituto de Investigaciones Económicas y Sociales del Sur—CONICET & Universidad Nacional del Sur Argentina, Bahía Blanca 8000, Argentina
18
Sección Limnología, IECA, Facultad de Ciencias, Universidad de la República, Montevideo 11400, Uruguay
*
Author to whom correspondence should be addressed.
Water 2018, 10(11), 1578; https://doi.org/10.3390/w10111578
Submission received: 11 September 2018 / Revised: 24 October 2018 / Accepted: 2 November 2018 / Published: 4 November 2018
(This article belongs to the Special Issue Current and Emerging Issues Surrounding Water in the Americas )

Abstract

:
In this work, we develop and test proxy-based diagnostic tools for comparing freshwater ecosystem services (FWES) risks across an international array of freshwater ecosystems. FWES threats are increasing rapidly under pressure from population, climate change, pollution, land use change, and other factors. We identified spatially explicit FWES threats estimates (referred to as threat benchmarks) and extracted watershed-specific values for an array of aquatic ecosystems in the Western Hemisphere (Ramsar sites). We compared these benchmark values to values extracted for sites associated with an international FWES threat investigation. The resulting benchmark threats appeared to provide a meaningful context for the diagnostic assessment of study site selection by revealing gaps in coverage of the underlying socio-environmental problem. In an effort to simplify the method, we tested regularly updated environmental and socioeconomic metrics as potential proxies for the benchmark threats using regression analysis. Three category proxies, aggregated from (i) external (global to regional, climate-related), (ii) internal (watershed management-related), and (iii) socioeconomic and governance related proxies produced strong relationships with water supply threat benchmarks, but only weak relationships with biodiversity-related and nutrient regulation benchmark threats. Our results demonstrate the utility of advancing global FWES status and threat benchmarks for organizing coordinated research efforts and prioritizing decisions with regard to international socio-environmental problems.

1. Introduction

Freshwater ecosystems provide our global society with provisioning (e.g., water and food supply), regulating (e.g., pollutant attenuation), cultural (e.g., iconic fishery conservation), and supporting services [1,2,3]. While critical to human and aquatic life and threatened globally, many freshwater ecosystem services (FWES) are complex and difficult to quantify [4,5,6]. Therefore, most ecosystem services (ES) assessments are site-specific, limiting our ability to compare across sites and make informed decisions about the allocation of research and management resources. To enhance our understanding of threats to FWES and make sound, transferable recommendations for mitigation, it is important to consider a broad range of socio-environmental conditions. Coordinated research networks provide new opportunities to develop understanding at larger spatial scales in the context of socio-environmental problems like threats to FWES [7,8,9,10]. However, optimizing site selection in such networks is also a challenging problem that has received little attention in the socio-environmental and ES literature.
Evaluating a network of sites for the study of FWES threats requires at least some knowledge of these threats and their spatial variation. Some knowledge may be gleaned from environmental quality and performance, water scarcity, biodiversity and other ES-related states, which are increasingly being assessed and mapped using socioeconomic, environmental, and human development data [11,12,13,14,15,16]. However, advancing toward understanding and mitigating FWES threats requires further investigation, as this is a “wicked problem” related to natural and anthropogenic factors, and to society’s capacity to manage those factors in terms of governance, values, and perceptions [17,18]. Most relevant to this work are the Riverthreat.net [14] and Aqueduct Water Risk Atlas [16] efforts to integrate available environmental, agronomic and human development metrics and simulation outputs to obtain FWES threat-related products at the watershed scale, which is the scale of interest most relevant to FWES decision-makers. The Riverthreat.net products include global maps of threat drivers (e.g., catchment disturbance, pollution, fishing pressure) and the estimated incident threat to human water security (HWS) and biodiversity. Both threats are relatable to FWES threats, as are several of the drivers. For example, Vörösmarty et al. [14] classify nitrogen loading as a driver for the threat to biodiversity, but nitrogen loading is also connected to nutrient regulation FWES. These and related approaches that take more FWES-relevant metrics and simulations into account show great promise in supporting the understanding of the distribution of FWES threats [19,20,21].
Approaches integrating hydro-climatic, agronomic, and other information to map environmental states have advanced markedly in the past decade and research is ongoing to provide reliable FWES and FWES threat estimates [2,4,6,14,15,16]. A key issue to address is that most such approaches fail to capture social dynamics that would help them advance toward forecasting environmental conditions or performances at the watershed scale and identify high threat/low capacity areas in need of attention. For example, diagnostic investigations and modifications may be warranted for advanced freshwater ecosystem monitoring and management programs, such as the European Union’s Water Framework Directive (WFD) [22,23,24,25]. Coordinated inquiries across socio-environmental gradients can help to address this shortcoming and, despite their limitations, current FEWS-related products may be adequate to support the design of coordinated research activity on global FWES threats and adaptation strategies. For a first approximation, societal conditions and aquatic state (hydroclimate, watershed management, etc.) should be useful for predicting FWES threat. It also follows, and others have noted, that some FWES threats are easier to characterize than others [5,6]. For example, provisioning services like water supply are universally valued and straightforward to assess. More complex FWES, such as nutrient regulation, are affected by social, hydro-climatic, watershed management and governance factors and more difficult to quantify. Some FWES, such as biodiversity-related ES, may be hampered by unanswered ecological questions rendering them much more difficult to assess.
In this paper, we develop and test proxy-based diagnostic tools for comparing FWES risks across an international (Pan-American) array of freshwater ecosystems. Our approach was to first identify meaningful benchmarks for threats to FWES, and then to construct and test more easily attainable proxies for these benchmarks by combining readily available and regularly updated indicators. We hypothesized that the reliability of FWES threat proxies will decrease from provisioning services (most reliable), to regulating services, to cultural and supporting services (least reliable). In addressing this hypothesis, we explored two supporting research questions: (1) Do the benchmark threats provide useful diagnostics in the context of a coordinated international research network? (2) Do meaningful benchmark threat proxy relationships exist that could offer similar threat diagnostics in the absence of the more effort-intensive benchmarks that are not necessarily updated over time?

2. Methods

Our approach involved three main tasks. First, we identified the best available FWES threat estimates for a large number of well-known aquatic ecosystems distributed throughout the Americas, including sites from an existing research network. Second, with respect to experimental design, we used the resulting threat estimates to assess whether the research network site as a subset of the greater FWES threat space adequately spanned that space. These tasks allowed us to address our first research question. Third, to address our second research question, we assembled potential threat proxies from the published literature on environmental and socioeconomic metrics and indicators, and used linear regression analysis to explore potential relationships between threat proxies and benchmarks. Regression approaches such as this one are commonly used to explore global socio-environmental problems and connect regularly available and updated metrics (e.g., Gross Domestic Product) with less available and less frequently updated environmental quality or performance indicators, e.g., [15,16].

2.1. Study Sites

To explore the variation in threats to FWES across the Americas we selected 32 sites in 23 nations from the international Convention on Wetlands (www.ramsar.org, Figure 1 and Table A1), including a mix of lakes, rivers and wetlands. The 32 sites included three sites each from larger nations with heterogeneous climate and population distributions (Canada, USA, Brazil), two sites from medium-sized nations (Argentina, Chile, Mexico), and one site from each of 17 smaller nations in Central and South America. Several smaller Caribbean nations were omitted because they lacked Ramsar sites, or indicators needed for the analysis were unavailable. We used the Ramsar sites to extract threat estimates from a global FWES threat dataset (Section 2.2), then used those results to diagnostically assess the international coordinated research effort “Sensing the Americas’ Freshwater Ecosystem Risk from Climate change” (SAFER Project, www.safer.conicet.gob.ar). The SAFER research network includes seven aquatic sites (rivers, lakes, and coastal lagoons) located in six nations (Figure 1). SAFER sites were selected mainly because they were sites of ongoing investigations by SAFER researchers and allowed the SAFER network to leverage prior work and additional funding. In this work, we examine their potential for exploring socio-environmental gradients to improve the understanding of FWES threats.
As part of the SAFER project, investigators held workshops with local experts and in some cases stakeholders to classify and rank key ecosystem services (ES) and identify threats to those services (Table 1). See Smyth et al. (this issue) for more details on the SAFER sites, particularly with respect to stakeholder engagement. For the purposes of this work, we selected three FWES common to all of the SAFER sites: (i) Water Provisioning, (ii) Biodiversity-related ES (broadly defined), and (iii) Nutrient regulation which reflect water quantity, ecological condition, and water quality, respectively.

2.2. Identifying Freshwater Ecosystem Service Threat Benchmarks

We identified the global Riverthreat.net spatial 0.5° gridded data products [14] as the best available FWES threat metrics (hereafter referred to as threat benchmarks). Other prospective products provided similar results but were available only at coarser resolutions [26,27,28,29]. We designated these as our benchmarks because (i) they emphasize freshwater ecosystems (primarily rivers) and are organized at the watershed scale, and (ii) they incorporate a comprehensive range of stressors related to climate and human influences. Specific to the three FWES targeted in this work (water supply, biodiversity conservation, nutrient regulation), we selected the adjusted Human Water Security threat (aHWS), incident Biodiversity threat (iBD), and the nitrogen loading (NL) driver as our focal threat benchmarks. The aHWS and iBD metrics are intended to summarize threats of inadequate human water supply and biodiversity loss in a watershed. These threats are based on accumulation of 23 weighted drivers in four thematic areas: catchment disturbance (four drivers); pollution (nine); water resource development (six); and biotic factors (four) [14]. The adjustment in the HWS score (leading to aHWS) is calculated from water resources development drivers associated with risk-ameliorating water infrastructure and management investments. The nitrogen loading (NL) score was categorized as one of the 23 drivers in the Riverthreat.net data product (nutrient regulation threats per se were not available), and we chose to employ it separately as a threat benchmark.
We extracted aHWS, iBD and NL benchmark values from the global dataset for each Ramsar and SAFER study site by averaging the five highest pixel values in contact with each site’s water body. Two of the SAFER site water bodies (La Salada and Laguna de Rocha) were small relative to the data set spatial resolution (0.5°). In these cases, we averaged only the two or three pixels in contact with the water bodies. We used the resulting Ramsar benchmark values to diagnostically assess the SAFER sites as a coordinated research network in terms of its coverage of the range of FWES threats observed in the broader set of sites.

2.3. Identifying Proxies for Threat Benchmarks

Because the benchmarks threats are (to date) one-time estimates, we tested a wide array of readily available hydrologic, environmental performance, and socioeconomic indicators as potential proxies for the benchmarks. If available, such proxies could facilitate site selection and prioritization in coordinated research networks. Prospective proxies included (Table 2): (i) indicators of water use, water stress/scarcity, and water vulnerability, (ii) scores from the Yale Environmental Performance Index (EPI) and its components, (iii) access to clean water and improved sanitation indicators, (iv) World Governance indicators, and (v) wealth as GDP per capita. All of these proxies are available only at the national scale with the exception of some of the hydroclimatic variables (WRI 2014 proxies; [14]), which are watershed-based values.
After preliminary assessment of a broader array of data sources, using methods described below, we settled on the most promising and reliable proxies (their readily available and regularly updated). We categorized the selected proxies in terms of (1) external threats related to regional hydroclimatic changes, (2) internal threats from human activities within a watershed, or (3) social threats, i.e., threats due to lack of resources, poor governance, and other socioeconomic factors (Table 2). Many of these indicators are associated with more than one of these categories, but we categorized them with respect to their primary influence. For example, upstream storage (STOR) is categorized here as an internal indicator because it is primarily associated with geomorphology and major watershed manipulations (e.g., reservoirs), but is also clearly driven by climate and adaptive capacity (e.g., arid climates require more water storage infrastructure to enable agricultural enterprises and require government policy to construct).
We normalized all non-normalized proxy values using their maximum theoretical value, resulting in a range from 0 to 1 (low to high threat, strong to poor performance, etc.). Proxies scored with the opposite convention were inverted for our purposes (1-proxy value in Table 2). The directionality of the performance or threat scale was unclear for normalized GDP per capita (GDPP). While increasing national wealth (GDP) is expected to be associated with lower threats to some services (e.g., water supply), it may be associated with elevated threats to other services (e.g., biodiversity-related threats). For this proxy, we made no assumption about directionality. The directionality of the threat scale is sometimes dependent of the ES in question. For example, increasing amounts of upstream storage in a watershed (STOR, Table 2) typically reduces the threat of loss of water supply services. However, storage is also likely to increase the threat of loss of biodiversity-related ES.
We used a three-step approach to search for FWES threat proxies that were a combination of those from the three categories (Figure 2). First, within each category, we considered the proxy values both individually and in combination with the other proxies in the category. Combinations were created by simple averaging and all combinations were tested. Second, we aggregated the three categories using both weighted means (Pwm) and geometric means (Pgm) as follows:
P w m = w 1 T e + w 2 T i + w 3 T s .
P g m = T e T i T s 3 ,
where Pwm and Pgm are referred to as three-category threat proxies based on the intermediate proxies for external, internal, and social categories (Te, Ti, Ts,) and w1, w2, w3 are weighting factors (valued 0 to 1 and summing to 1). Third, we used regression analysis to identify the strongest linear relationships (highest R2) between the threat benchmarks as the independent variables (aHWS, iBD, NL) and the 3-category proxy as the dependent variable (Pwm or Pgm). As with the benchmarks, we used the resulting regression relationships to diagnostically assess the SAFER network of sites in terms of its coverage of the proxy-benchmark threat space.

3. Results

3.1. Freshwater Ecosystem Service Threat Benchmarks

The FWES threat benchmarks for water supply (aHWS), biodiversity-related services (iBD), and nutrient regulation (N-load) ranged from relatively minor (~0.3) to major threats (~1.0) in all three cases (Figure 3). For the Ramsar site aHWS values (Figure 3a), we found few low-threat sites (aHWS < 0.4) and a relatively high mean threat (0.67). The aHWS values were roughly normally distributed (Figure 3a inset graph), though slightly skewed toward lower values. The sites causing the tailing at the low end are all highly developed nations (CA and US). With respect to SAFER project diagnostics, the SAFER site aHWS values ranged from 0.35 to 0.75, encompassing the lower two-thirds of the range exhibited by the Ramsar sites, as evidenced by the lower mean value for the SAFER sites (0.53). There was a lack of coverage of the higher water supply threat range by the SAFER sites.
The biodiversity-related services threat benchmark (iBD) for the Ramsar sites (Figure 3b) ranged from 0.39 to 0.93, with an average (0.66) suggesting moderately high threat levels at most sites. Similar to the aHWS case, the Ramsar iBD values were roughly normally distributed and skewed toward lower values. In this case, the sites associated with lower iBD values were mainly in remote, undeveloped areas. In contrast to aHWS, the higher end of iBD range included most of the Ramsar sites in highly developed nations. For the iBD-related threat benchmark, the SAFER site threat values represented good coverage of the entire range presented by the Ramsar sites.
The nutrient regulation threat benchmark (N-load) for the Ramsar sites (Figure 3c) ranged from 0.29 to 0.85, with a slightly lower average (0.58) than the other two benchmarks. The Ramsar N-load values were normally distributed (skewness = −0.02). The SAFER site values were skewed toward the intermediate and high threat portions of the range, but fail to appear in the lower portion.

3.2. Proxies for Freshwater Ecosystem Service Threat Benchmarks

Multiple combinations of the external, internal, and social proxies resulted in strong relationships with the water supply stress benchmark (aHWS). In the simplest case, national wealth (GDPP) as an independent proxy accounted for 72% of the variation in aHWS across the Ramsar sites (Figure 4a). The majority of the Ramsar sites are from developing nations, which led to the grouping of sites between about 0.1 and 0.5 (x-axis, PaHWS in Figure 3a), and a large gap in coverage between 0.5 and 0.8, with highly-developed nations (Canada and USA) falling above 0.8. Notably, aggregation of prospective governance indicators with GDPP did not increase the regression coefficient in this case. The strongest relationship with the Ramsar aHWS values (R2 = 0.92) was provided by the following three-category proxy (Figure 4b):
P a H W S = 0.2 D R O + 0.5 ( 1 W R I ) , W A T S U P , A C S A T ¯ + 0.3 G D P P ,
where DRO is the drought frequency proxy (external), ( 1 W R I ) , W A T S U P , A C S A T ¯ is the average of the water return index, access to clean water, and access to sanitation (internal) proxies, and GDPP is the gross domestic product per capita (social) proxy.
Five of the seven SAFER sites were within the 95% confidence bands developed from the Ramsar site data for GDPP-only proxy case. For the more complex proxy there was a much tighter confidence interval encompassing four SAFER sites, with the remaining three just outside the lower confidence band. These results suggest that this benchmark-proxy relationship could be useful for informing the design of a coordinated research network. The SAFER sites displayed the same grouping behavior mentioned above for development level. This outcome suggested a potential limitation of the resulting relationships in these equations due to an overabundance of sites characterized by roughly the same GDPP.
Viable relationships between the benchmark threats to biodiversity services (iBIO) and the prospective proxies were more challenging to identify relative to those for water supply. We identified one relatively weak (R2 = 0.31) but significant relationship for a three-category threat proxy (PiBIO) defined by (Figure 4):
P i B D = 1 H F O , 1 D R O , 1 W S V , 1 S V ¯ × S T O R , W R I ¯ × R L , G D P P ¯ 3 ,
where 1 H F O , 1 D R O , 1 W S V , 1 S V ¯ is the average of all four of the hydro-climate (external) proxies, S T O R , W R I ¯ is the average of the watershed storage and water return (internal) proxies, and R L , G D P P ¯ is the average of the Rule of Law and GDPP (social) proxies. The regression slope is positive, meaning that any increase in PiBD predicts an increase in the threat to biodiversity.
The utility of PiBD as a comparative metric for improving the SAFER project is questionable relative to the case for PaHWS (Figure 3b). Only two of the seven SAFER sites (San Joaquin River, USA and La Salada, AR) fall within the confidence bands. Furthermore, the SAFER sites are narrowly clustered (0.3 < PiBD < 0.5) and appear to be aligned along a steeper slope than the Ramsar sites. Lastly, the scale of the optimal proxy is compressed here, achieving the maximum threat (iBD = 1.0) at a PiBD value of 0.7. This compression could be an artifact of the iBD indicator since biodiversity threat classification is itself a challenging topic [36].
Relationships between the benchmark threats to nutrient regulation services (NL) and proxies (PNL) were the most challenging to identify. One significant but weak relationship (R2 = 0.22) proxy for nitrogen loading (PNL) was found (Figure 6):
P N L = ( 1 W S V ) , ( 1 S V ) ¯ × N U E , N B a l ¯   × R L , G D P P ¯ 3 ,
where ( 1 W S V ) , ( 1 S V ) ¯ is the average of the inter-annual and seasonal precipitation variability (external) proxies, N U E , N B a l ¯ is the average of the nitrogen use efficiency and nitrogen balance (internal) proxies, and R L , G D P P ¯ is the average of the Rule of Law and GDPP (social) proxies. In this case, the response between the proxy and benchmark threats was weak (slope = 0.35) relative to those for iBD (0.67) and aHWS (−1.23). The 95% confidence bands for the N-loading regression were broad due to the relatively weak relationship, and hence it is not surprising that most of the SAFER sites (six of seven) fell within the bands. The SAFER site proxy values provided good coverage of the upper half of the proxy space (0.5 < PiBD < 0.9) but were absent in the lower range.

4. Discussion

4.1. Freshwater Ecosystem Service Threat Benchmarks

Our first research question concerned whether the benchmarks can provide useful diagnostics with respect to a coordinated research network. The threat benchmarks extracted for the Ramsar sites demonstrate the promise of this approach for identifying meaningful comparative FWES threats metrics across an international array of sites. Superposition of the SAFER sites benchmark values on the spectrum of Ramsar site values (Figure 3a–c) reveals that the SAFER network covers some parts of the more general range of threat behavior better than others. For water supply threats (aHWS), the results point to the need for study sites in higher threat areas (Figure 3a), which mainly coincide with the less well-developed nations in our study (e.g., Haiti, Ecuador, El Salvador, Dominican Republic). There are logistical and other challenges associated with including such nations, and these must be overcome if the goal is to test a wider range of threats. In contrast, for nutrient regulation threats (NL, Figure 3c) the SAFER sites provide better coverage of the high and intermediate threat levels, but fail to cover the lowest levels. Improving coverage in this regard may present the challenge of identifying funding sources for study sites where threats to FWES are relatively low.
As noted at the outset of this paper, identifying the real threat to FWES is a complex and site-specific socio-environmental problem. The outcomes here are considered first approximations for the purpose of better aligning future studies with the diversity of threats present. That said, our benchmarking approach could provide diagnostic support for decisions regarding inclusion of additional sites to a research network should research capacity and local interests allow it.

4.2. Proxies for Threat Benchmarks

The best proxies-benchmark relationships (Figure 4b, Figure 5 and Figure 6) provided some support for our research hypothesis. Proxies exhibited strong relationships with the water supply benchmark (aHWS), while only relatively weaker relationships were identified for proxies and the benchmarks for biodiversity-related (iBD) and nutrient regulation (NL) service threats. This outcome is likely related to the complexity of the respective ecosystem services, and also pertains to our second research question regarding the meaningfulness of our proxy benchmark relationships.
For water provisioning ES threats, our three-category proxy was strongly (R2 = 0.85) and inversely related to the benchmark aHWS (Figure 4b), indicating a relatively sharp decline in risk per unit increase in proxy. The best-fitting proxy included drought occurrence (DRO), the non-return (natural) flow ratio (1-WRI), access to clean water (WATSUP) and advanced sanitation (ACSAT), and GDP per capita (GDPP). Although it would be speculative to attribute cause and effect to these simple regression results, the constituting proxies here could be meaningful. Drought prone regions with sufficient financial resources generally develop water resources infrastructure to provide a buffer against drought impacts (often at the expense of the aquatic ecosystems). Increasing access to clean water and sanitation has long been used to document human development and could be interpreted as another indicator for improving watershed management. All of these conditions are consistent with a reduced threat to the water provisioning service.
For the biodiversity-related ES threat (iBD), the best proxy-benchmark relationship was much weaker than that for water supply (R2 = 0.31, Figure 5) and therefore less reliable. The three-category proxy identified is consistent with the rationale that threats to these services increases with decreasing hydro-climate variability, increasing water storage and return flows in the watershed, and increasing national wealth and Rule of Law. These conditions are consistent with development and likely urbanization, which are well-known to result in permanent habitat loss and other threats to biodiversity-related services [37,38]. However, given the weakness of the relationship, it is more likely the case that better and/or more consistently assessed metrics for biodiversity-related threats are needed before meaningful proxies can be identified.
The benchmark-proxy threat relationship for nutrient regulation (NL) was the weakest of the three FWES threats tested (R2 = 0.22, Figure 6). The large degree of scatter combined with the weak positive slope in this relationship reveal that nitrogen loading does not differ greatly across the spectrum of sites. Less variable climate and increased national wealth correspond to higher threat values, which may correspond to increased agriculture and fertilizer usage. The appearance of nitrogen-use efficiency (NUE) and nitrogen balance (NBal) from the internal (watershed management) proxy category appears to be counterintuitive. Higher values for these metrics correspond to better agronomic management of nitrogen on watershed farms. Thus, one would expect the nitrogen loading to decrease with increase NUE and/or NBal. If this benchmark-proxy threat relationship is valid, then it suggests that these agronomic metrics for nitrogen use may also be indicative of nitrogen-related threats to nutrient regulation, as from widespread nitrogen application in a watershed [39].

4.3. Testing for Nationally Clustered Behavior

As noted in Section 3, in some cases benchmark values for the sites appeared to cluster with respect to national factors. We did not analyze the clustering behavior exhaustively, but we began to explore it with respect to the World Bank’s GDP estimates and governance quality indicators (Table 2). More specifically, we tested for clusters in nations using historical values (1996–2015) of annual GDP and the four governance indicators (R Package https://uc-r.github.io/kmeans_clustering). Three distinct clusters of nations emerged using the historical trends in GDPP, Rule of Law (RL) and access to advanced sanitation (ACSAT) (see Figure A1). These can be characterized as higher (Cluster 1), intermediate (Cluster 2) and lower (Cluster 3) clusters of wealth, sanitation access, and Rule of Law. Reanalyzing 32 Ramsar sites to create three benchmark-proxy correlations using these clusters revealed further insight into the threat benchmark-proxy relationships.
Many of the new benchmark-proxy relationships for the clusters span narrow ranges and hence it is important to avoid over-interpretation. With this qualification, the results suggest that there are “sub-trends” within the overall trends (Figure 7) that may warrant further investigation. These sub-trends may be related to factors like national development level and governance quality. For instance, the first and second clusters for the water supply threat (aHWS) resulted in significant negative relationships, while the third cluster failed to yield a relationship (Figure 7 left column). Furthermore, if the trends are accurate, they point to stronger negative relationship (based on the slope) for the more developed nations relative to intermediate nations. These results suggest that the connection between the external, internal, and socioeconomic/governance proxy indicators and the water provisioning service benchmark weakens with decreasing development levels in regard to national wealth, rule of law, and access to sanitation.
For the biodiversity-based benchmark threats (iBD), the first and second clusters resulted in positive slopes, similar to the overall results (Figure 4b), while the third cluster again failed to yield a significant relationship (Figure 7 center column). In this case, if accurate, the slopes for the second cluster were more positive than for the first cluster. These results suggest that the connection between the proxy indicators and the biodiversity-based services threat benchmark weakens with increasing development level in regard to national wealth, rule of law, and access to sanitation. However, the abundance of elevated biodiversity threat values for the rightmost portion of the plot also suggests that many of these sites are already experiencing elevated threats and hence further development cannot increase these threats substantially.
The nutrient regulation ES threats (NL) resulted in significant positive and negative relationships with the proxy for the first and third clusters, respectively, while cluster 2 failed to yield a significant relationship (Figure 7, right column). From cluster 1, we see that wealthier nations with relatively strong environmental performance (ACSAT) and Rule of Law appear to show an increase in nitrogen loading. With the United States, Canada and other nations in the group, increasing population and wealth led to increased agricultural production supported by the widespread fertilizer application. From cluster 3, we see nitrogen loading decreasing as PNL increases. This trend would be consistent with a stage in development in which farming pressure in the form of nitrogen loading on waterways is not yet evident, either due to lack of agricultural development or lack of data [39].

5. Conclusions

We extracted benchmark threat estimates for three FWES from a global data set using an array of Ramsar aquatic ecosystems distributed throughout the Americas. The resulting benchmark threats provided a meaningful context for a posteriori diagnostic assessment of study site selection for the SAFER project coordinated research network. To our knowledge, this is the first demonstration of the potential utility of such an approach to examine and optimize socio-environmental research site selection. With new global data sets becoming available [40,41,42,43], future work should compare them in this context.
We also explored the potential for using a wide array of readily accessible and regularly updated proxies to estimate the same three benchmark threats. We did this because the benchmark threats are currently one-time snapshot estimates, and proxies would therefore be useful moving forward to apply this now-dated large-scale initiative to current and future FWES assessments. Recognizing that estimating FWES threats is a complex socio-environmental problem, we proposed aggregating proxy values representing external (regional climatic), internal (watershed management-related), and social (wealth and governance-related) factors. We hypothesized that this approach would yield stronger relationships for relatively simple FWES than for relatively complex FWES. Our results for three FWES threats supported this hypothesis with water provisioning benchmark threats yielding a much stronger relationship with proxies than those for either biodiversity-related services or nutrient regulation.
Although the proxy approach proposed is promising, the unreliable outcomes for two of the three ES tested confirm that alternative data are needed if this approach is to apply to a broad array of FWES. Future research would involve identifying and testing new prospective data sources. Such data may include high resolution regional climate products with stronger connections to FWES threats. The amount and quality of such data varies globally, which is another challenge for applications that span countries and hemispheres. For example, we have observed that the climatic pressures on FWES have been scrutinized much more in some regions of the SAFER project (e.g., SJR, California) than others (e.g., SR, Argentina). Additional types of and higher granularity socioeconomic and governance data sources also need to be explored as our proxies were mainly national values, which may be biased toward urban conditions. In contrast, many of the aquatic ecosystems in this study are in relatively rural settings, which can lag in terms of environmental policy enforcement [44,45]. Lastly, our approach failed to capture the likely connections between governance, political factors, and local stakeholders. Understanding these connections is critical to managing socio-environmental problems and is detailed in a related SAFER project paper [46]. Disconnects between environmental policy and implementation that can bias environmental assessments have also been identified in Chile [47] and elsewhere. Thus, additional efforts aimed at validating data and indexes calculated for socio-environmental performance or problems are warranted.
Overall, our results demonstrate the utility of advancing the state of global FWES status and threat benchmarks for organizing coordinated research efforts and prioritizing management decision-making. Identifying risks to freshwater ecosystem services will require intensive coordinated research efforts into the nature of the FWES threats and how best to monitor those threats and communicate them to stakeholders. This study identifies important differences in how easily those threats can be recreated from proxies when data are lacking, and which FWES could benefit most from additional research.

Author Contributions

Conceptualization, T.C.H., R.L.S., D.C., J.P.L., M.N.; Methodology, T.C.H., R.L.S., M.Z.; Formal Analysis, R.D., T.C.H.; Writing—Original Draft Preparation, T.C.H.; Writing—Review & Editing, T.C.H., R.L.S., S.C., D.C., J.E., N.H., J.P.L., G.M.E.P., S.P., M.C.P., B.R., J.A.R., F.S., M.I.V., S.R.V., B.W., M.Z.; Visualization, T.C.H., N.H. S.R.V.; Supervision, G.M.E.P., M.C.P. T.C.H.; Project Administration, G.M.E.P.; Funding Acquisition, T.C.H., G.M.E.P., J.A.R., B.R., M.I.V., J.E.

Funding

This work was supported by the Inter-American Institute for Global Change Research (IAI) CRN3038 under U.S. National Science Foundation (NSF) Award GEO-1128040, and by the NSF Water Sustainability & Climate and Science across Virtual Institutes programs (Award CBET #1204841 and #1336839). J.E. and N.H. were partially funded by Canadian Queen Elizabeth II Diamond Jubilee Scholarships (QES), a partnership among Universities in Canada, the Rideau Hall Foundation (RHF), and Community Foundations of Canada (CFC). The QES-AS is made possible by financial support from IDRC and SSHRC. N.H. was partially supported by the Fulbright Visiting Scholar Program from the Fulbright Commission in Colombia while on research leave from Universidad del Norte.

Acknowledgments

The authors thank F.E.H. for proofreading the final version of the manuscript, and the editor and anonymous reviewers for greatly improving the final product.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analysis, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A.

Table A1. Ramsar convention sites employed (three for larger nations, two for medium-sized nations, and one for all other nations).
Table A1. Ramsar convention sites employed (three for larger nations, two for medium-sized nations, and one for all other nations).
Site NameNationCoordinatesRamsar Sites Link
Beaverhill LakeCanada53°30′ N, 113°30′ Whttps://rsis.ramsar.org/ris/370
Lake ChamplainCanada, US44°57′ N, 73°10′ Whttps://rsis.ramsar.org/ris/2200
Lac Saint-FrancoisCanada45°02′ N, 74°29′ Whttps://rsis.ramsar.org/ris/361
Laguna de Santa RosaUS38°24′ N, 122°47′ Whttps://rsis.ramsar.org/ris/1930
Caddo LakeUS32°45′ N, 94°08′ Whttps://rsis.ramsar.org/ris/633
Laguna de ZapotlánMexico19°45′ N, 103°29′ Whttps://rsis.ramsar.org/ris/1466
Manglares y Humedales de TuxpanMexico21°00′ N, 097°21′ Whttps://rsis.ramsar.org/ris/1602
Lago AtitlánGuatemala15°25′ N, 89°22′ WNot a Ramsar site
Lago YojoaHonduras14°51′ N, 88°00′ Whttps://rsis.ramsar.org/es/ris/1467
Lago ArenalCosta Rica0°30′ N, 84°51′ Whttps://rsis.ramsar.org/ris/1022
Golfo de MontijoPanama7°45′ N, 81°07′ Whttps://rsis.ramsar.org/ris/510
Laguna de OlomegaEl Salvador13°19′ N, 88°04′ Whttps://rsis.ramsar.org/ris/1899
Sistema de Humedales de San MiguelitoNicaragua11°25′ N, 84°51′ Whttps://rsis.ramsar.org/ris/1140
Laguna de LecheCuba22°19′ N, 78°29′ Whttps://rsis.ramsar.org/ris/1235
Lago EnriquilloDominican Republic18°28′ N, 71°39′ Whttps://rsis.ramsar.org/es/ris/1179
Lac AzuéiHaiti18°35′ N, 72°0′ WNot a Ramsar site
Laguna de CochaColombia01°03′ N, 77°12′ Whttps://rsis.ramsar.org/es/ris/1047
Caroni SwampTrinidad and Tobago10°34′ N, 61°27′ Whttps://rsis.ramsar.org/ris/1497
Parque Nacional CajasEcuador02°50′ N, 79°14′ Whttps://rsis.ramsar.org/ris/1203
Lago TiticacaBolivia, Peru16°10′ S, 68°52′ Whttps://rsis.ramsar.org/ris/959
Cienega de los OlivitosVenezuela10°55′ N, 71°26′ Whttps://rsis.ramsar.org/es/ris/859
Lagoa do PeixeBrazil31°14′ S, 50°57′ Whttps://rsis.ramsar.org/ris/603
Ilha do BananalBrazil10°31′ S, 50°12′ Whttps://rsis.ramsar.org/ris/624
MamirauáBrazil2°18′ S, 66°02′ Whttps://rsis.ramsar.org/ris/623
Laguna BlancaArgentina39°02′ S, 70°21′ Whttps://rsis.ramsar.org/ris/556
Laguna de LlancaneloArgentina35°45′ S, 69°08′ Whttps://rsis.ramsar.org/ris/759
Bañados del Este y Franja CosteraUruguay33°48′ S, 53°50′ Whttps://rsis.ramsar.org/ris/290
Salar de TaraChile22°56′ S, 67°15′ Whttps://rsis.ramsar.org/ris/875
Carlos Anwandter SanctuaryChile39°41′ S, 73°11′ Whttps://rsis.ramsar.org/ris/222
Lago YpoaParaguay26°30′ S, 57°33′ Whttps://rsis.ramsar.org/ris/728
Table A2. Regression outcomes for Figure 7, with nations clustered by historical trends in GDPP, access to advanced sanitation (ACSAT), and Rule of Law (RL), where the regression coefficient (trend line slope).
Table A2. Regression outcomes for Figure 7, with nations clustered by historical trends in GDPP, access to advanced sanitation (ACSAT), and Rule of Law (RL), where the regression coefficient (trend line slope).
Sites (by Nation)Benchmark IndicatorSlope (Regr. Coeff.)R2p-Value
Cluster 1Canada (3), Chile (2), Costa Rica, Trinidad and Tobago, United States (3), UruguayaHWS
iBIO
NL
−1.26
0.56
0.45
−0.85
0.61
0.62
<0.001
0.05
0.05
Cluster 2Argentina (2), Brazil (3), Colombia, Cuba, Dominican Republic, Ecuador, Honduras, Mexico (2), Panama, Paraguay, Peru, VenezuelaaHWS
iBIO
NL
−0.82
1.20
0
−0.51
0.68
0
0.05
<0.01
NS
Cluster 3Bolivia, El Salvador, Guatemala, Haiti, NicaraguaaHWS
iBIO
NL
−0.47
−1.89
−1.15
−0.49
−0.45
−0.81
NS
NS
0.05
Figure A1. Results from cluster analysis of study nations based on historical values (1996–2015) of GDPP, ACSAT, and RL (R Package https://uc-r.github.io/kmeans_clustering).
Figure A1. Results from cluster analysis of study nations based on historical values (1996–2015) of GDPP, ACSAT, and RL (R Package https://uc-r.github.io/kmeans_clustering).
Water 10 01578 g0a1

References

  1. De Groot, R.S.; Wilson, M.A.; Boumans, R.M.J. A typology for the classification, description and valuation of ecosystem functions, goods and services. Ecol. Econ. 2002, 41, 393–408. [Google Scholar] [CrossRef] [Green Version]
  2. Millennium Ecosystem Assessment Panel (MEAP). Ecosystems and Human Well Being: Synthesis; Island Press: Washington, DC, USA, 2005; ISBN 1-59726-040-1. [Google Scholar]
  3. Haines-Young, R.; Potschin, M. Common International Classification of Ecosystem Services CICES Version 4: Response to Consultation; Center for Environmental Management, University of Nottingham: Nottingham, UK, 2012. [Google Scholar]
  4. Daily, G.C.; Polasky, S.; Goldstein, J.; Kareiva, P.M.; Mooney, H.A.; Pejchar, L.; Ricketts, T.H.; Salzman, J.; Shallenberger, R. Ecosystem services in decision making: Time to deliver. Front. Ecol. Environ. 2009, 7, 21–28. [Google Scholar] [CrossRef]
  5. Bull, J.W.; Jobstvogt, N.; Böhnke-Henrichs, A.; Mascarenhas, A.; Sitas, N.; Baulcomb, C.; Lambini, C.K.; Rawlins, M.; Baral, H.; Zähringer, J.; et al. Strengths, Weaknesses, Opportunities and Threats: A SWOT analysis of the ecosystem services framework. Ecosyst. Serv. 2016, 17, 99–111. [Google Scholar] [CrossRef]
  6. Maron, M.; Mitchell, M.G.E.; Runting, R.K.; Rhodes, J.R.; Mace, G.M.; Keith, D.A.; Watson, J.E.M. Towards a threat assessment framework for ecosystem services. Trends Ecol. Evol. 2017, 32, 240–248. [Google Scholar] [CrossRef] [PubMed]
  7. Redman, C.L.; Grove, J.M.; Kuby, L.H. Integrating social science into the Long-Term Ecological Research (LTER) Network: Social dimensions of ecological change and ecological dimensions of social change. Ecosystems 2004, 7, 161–171. [Google Scholar] [CrossRef]
  8. Montgomery, J.L.; Harmon, T.; Haas, C.N.; Hooper, R.; Clesceri, N.L.; Graham, W.; Kaiser, W.; Sanderson, A.; Minsker, B.; Schnoor, J.; et al. The waters network: An integrated environmental observatory network for water research. Environ. Sci. Technol. 2007, 41, 6642–6647. [Google Scholar] [CrossRef] [PubMed]
  9. Keller, M.; Schimel, D.S.; Hargrove, W.W.; Hoffman, F.M. A continental strategy for the National Ecological Observatory Network. Front. Ecol. Environ. 2008, 6, 282–284. [Google Scholar] [CrossRef] [Green Version]
  10. Hanson, P.C.; Weathers, K.C.; Kratz, T.K. Networked lake science: How the Global Lake Ecological Observatory Network (GLEON) works to understand, predict, and communicate lake ecosystem response to global change. Inland Waters 2016, 4, 543–554. [Google Scholar] [CrossRef]
  11. Ebert, U.; Welsch, H. Meaningful environmental indices: A social choice approach. J. Environ. Econ. 2004, 47, 270–283. [Google Scholar] [CrossRef]
  12. Niemeijer, D.; de Groot, R.S. A conceptual framework for selecting environmental indicator sets. Ecol. Indic. 2008, 8, 14–25. [Google Scholar] [CrossRef]
  13. Raudsepp-Hearne, C.; Peterson, G.D.; Tengoe, M.; Bennett, M.E.; Holland, T.; Benessaiah, K.; MacDonald, G.K.; Pfeifer, L. Untangling the environmentalist’s paradox: Why is human well-being increasing as ecosystem services degrade? BioScience 2010, 60, 576–589. [Google Scholar] [CrossRef]
  14. Vörösmarty, C.J.; McIntyre, P.B.; Gessner, M.O.; Dudgeon, D.; Prusevich, A.; Green, P.; Glidden, S.; Bunn, S.E.; Sullivan, C.A.; Reidy Liermann, C.; et al. Global threats to human water security and river biodiversity. Nature 2010, 467, 555. [Google Scholar] [CrossRef] [PubMed]
  15. Hsu, A.; Emerson, J.; Levy, M.; de Sherbinin, A.; Johnson, L.; Malik, O.; Schwartz, J.; Jaiteh, M. The 2014 Environmental Performance Index. New Haven, CT: Yale Center for Environmental Law and Policy. Available online: http://www.epi.yale.edu. (accessed on 15 November 2017).
  16. Gassert, F.; Rieg, P.; Shiao, T.; Luck, M.; Landis, M. Aqueduct Global Maps 2.1 Indicators: Constructing Decision-Relevant Global Water Risk Indicators, Working Paper, 2015. Available online: http://www.wri.org/publication/aqueduct-global-maps-21-indicators (accessed on 11 May 2017).
  17. Cosens, B.; Gunderson, L.; Allen, C.; Benson, M.; Cosens, B.; Gunderson, L.; Linkov, I. Identifying legal, ecological and governance obstacles, and opportunities for adapting to climate change. Sustainability 2014, 6, 2338–2356. [Google Scholar] [CrossRef]
  18. Angeler, D.G.; Allen, C.R.; Garmestani, A.S.; Gunderson, L.H.; Linkov, I. Panarchy use in environmental science for risk and resilience planning. Environ. Syst. Decis. 2016, 36, 225–228. [Google Scholar] [CrossRef] [Green Version]
  19. Burkhard, B.; Kroll, F.; Nedkov, S.; Müller, F. Mapping ecosystem service supply, demand and budgets. Ecol. Indic. 2012, 21, 17–29. [Google Scholar] [CrossRef]
  20. Albert, C.; Bonn, A.; Burkhard, B.; Daube, S.; Dietrich, K.; Engels, B.; Frommer, J.; Götzl, M.; Grêt-Regamey, A.; Job-Hobenet, B.; et al. Towards a national set of ecosystem service indicators: Insights from Germany. Ecol. Indic. 2016, 61, 38–48. [Google Scholar] [CrossRef]
  21. Diamantini, E.; Lutz, S.R.; Mallucci, S.; Majone, B.; Merz, R.; Bellin, A. Driver detection of water quality trends in three large European river basins. Sci. Total Environ. 2018, 612, 49–62. [Google Scholar] [CrossRef] [PubMed]
  22. Liefferink, D.; Wiering, M.; Uitenboogaart, Y. The EU Water framework directive: A multi-dimensional analysis of implementation and domestic impact. Land Use Policy 2011, 28, 712–722. [Google Scholar] [CrossRef]
  23. Collins, A.; Ohandja, D.-G.; Hoare, D.; Voulvoulis, N. Implementing the water framework directive: A transition from established monitoring networks in England and Wales. Environ. Sci. Policy 2012, 17, 49–61. [Google Scholar] [CrossRef]
  24. Arle, J.; Mohaupt, V.; Kirst, I.; Arle, J.; Mohaupt, V.; Kirst, I. Monitoring of surface waters in Germany under the water framework directive—A review of approaches, methods and results. Water 2016, 8, 217. [Google Scholar] [CrossRef]
  25. Voulvoulis, N.; Arpon, K.D.; Giakoumis, T. The EU Water Framework Directive: From great expectations to problems with implementation. Sci. Total Environ. 2017, 575, 358–366. [Google Scholar] [CrossRef] [PubMed]
  26. Abell, R.; Thieme, M.L.; Revenga, C.; Bryer, M.; Kottelat, M.; Bogutskaya, N.; Coad, B.; Mandrak, N.; Balderas, S.C.; Bussing, W.; et al. Freshwater ecoregions of the world: A new map of biogeographic units for freshwater biodiversity conservation. BioScience 2008, 58, 403–414. [Google Scholar] [CrossRef]
  27. Carr, G.M.; Neary, J.P. Water Quality for Ecosystem and Human Health, 2nd ed.; UNEP/Water Programme: Burlington, ON, Canada, 2008; ISBN 92-95039-51-7. [Google Scholar]
  28. Seitzinger, S.P.; Harrison, J.A.; Dumont, E.; Beusen, A.H.W.; Bouwman, A.F. Sources and delivery of carbon, nitrogen, and phosphorus to the coastal zone: An overview of Global Nutrient Export from Watersheds (NEWS) models and their application. Glob. Biogeochem. Cycles 2005, 19. [Google Scholar] [CrossRef] [Green Version]
  29. Groombridge, B.; Jenkins, M.D. Freshwater Biodiversity: A Preliminary Global Assessment—UNEP. WCMC Biodiversity Series 8; UNEP-WCMC: Cambridge, UK, 1998; ISBN 1-899628-12-6. [Google Scholar]
  30. Center for International Earth Science Information Network (CIESIN)—Columbia University, and Centro Internacional de Agricultura Tropical (CIAT). Gridded Population of the World, Version 3 (GPWv3): Population Density Grid, Future Estimates; NASA Socioeconomic Data and Applications Center (SEDAC): Palisades, NY, USA, 2005.
  31. Dartmouth Flood Observatory. 2011. Available online: http://floodobservatory.colorado.edu/ (accessed on 20 January 2018).
  32. Sheffield, J.; Wood, E.F. Characteristics of global and regional drought, 1950–2000: Analysis of soil moisture data from off-line simulation of the terrestrial hydrologic cycle. J. Geophys. Res. Atmos. 2007, 112, 1–21. [Google Scholar] [CrossRef]
  33. FAO AQUASTAT 2012. Available online: http://www.fao.org/nr/water/aquastat/main/index.stm (accessed on 10 February 2018).
  34. World Bank, 2014 Governance. Available online: http://info.worldbank.org/governance/wgi/index.aspx#home (accessed on 30 November 2017).
  35. World Bank, 2014 National Accounts Data. Available online: https://data.worldbank.org/indicator/NY.GDP.PCAP.CD (accessed on 2 December 2017).
  36. Chiarucci, A.; Bacaro, G.; Scheiner, S.M. Old and new challenges in using species diversity for assessing biodiversity. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2011, 366, 2426–2437. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. McKinney, M.L. Effects of urbanization on species richness: A review of plants and animals. Urban Ecosyst. 2008, 11, 161–176. [Google Scholar] [CrossRef]
  38. Pauchard, A.; Aguayo, M.; Peña, E.; Urrutia, R. Multiple effects of urbanization on the biodiversity of developing countries: The case of a fast-growing metropolitan area (Concepción, Chile). Biol. Conserv. 2006, 127, 272–281. [Google Scholar] [CrossRef]
  39. Lu, C.; Tian, H. Global nitrogen and phosphorus fertilizer use for agriculture production in the past half century: shifted hot spots and nutrient imbalance. Earth Syst. Sci. Data 2017, 9, 181–192. [Google Scholar] [CrossRef] [Green Version]
  40. Seitzinger, S.P.; Mayorga, E.; Bouwman, A.F.; Kroeze, C.; Beusen, A.H.W.; Billen, G.; Drecht, G.V.; Dumont, E.; Fekete, B.M.; Garnier, J.; et al. Global river nutrient export: A scenario analysis of past and future trends. Glob. Biogeochem. Cycles 2010, 24. [Google Scholar] [CrossRef] [Green Version]
  41. Bouwman, L.; Goldewijk, K.K.; Hoek, K.W.V.D.; Beusen, A.H.W.; Vuuren, D.P.V.; Willems, J.; Rufino, M.C.; Stehfest, E. Exploring global changes in nitrogen and phosphorus cycles in agriculture induced by livestock production over the 1900–2050 period. PNAS 2013, 110, 20882–20887. [Google Scholar] [CrossRef] [PubMed]
  42. Lehner, B.; Grill, G. Global river hydrography and network routing: Baseline data and new approaches to study the world’s large river systems. Hydrol. Process. 2013, 27, 2171–2186. [Google Scholar] [CrossRef]
  43. Dallaire, C.O.; Lehner, B.; Sayre, R.; Thieme, M. A multidisciplinary framework to derive global river reach classifications at high spatial resolution. Environ. Res. Lett. 2018. [Google Scholar] [CrossRef]
  44. McGranahan, G.; Satterthwaite, D.; Tacoli, C. Rural-Urban Change, Boundary Problems and Environmental Burdens. Human Settlements Working Paper Series Rural-Urban Interactions and Livelihood Strategies No. 10; IIED: London, UK, 2004; ISBN 978-1-84369-538-7. [Google Scholar]
  45. Wang, J.; Webber, M.; Finlayson, B.; Barnett, J. Rural industries and water pollution in China. J. Environ. Manag. 2008, 86, 648–659. [Google Scholar] [CrossRef] [PubMed]
  46. Smyth, R.L.; Ratima, U.; Segarra, M.; Borre, L.; Reid, B.; Pincetl, S.; Zilio, M.I.; Astorga, A.; Huamantinco Cisneros, M.A.; Conde, D.; et al. Engaging stakeholders in water research across socio-environmental gradients in the Americas. Water 2018, in press. [Google Scholar]
  47. Barandiaran, J. Chile’s environmental assessments: Contested knowledge in an emerging democracy. Sci. Cult. 2015, 24, 251–271. [Google Scholar] [CrossRef]
Figure 1. Map showing estimated human population density [30] and the locations of Ramsar sites (black-in-white circles) and SAFER project sites (white-in-black circles) used in this study (see Ramsar site list in Table A1; see SAFER abbreviations in Table 1).
Figure 1. Map showing estimated human population density [30] and the locations of Ramsar sites (black-in-white circles) and SAFER project sites (white-in-black circles) used in this study (see Ramsar site list in Table A1; see SAFER abbreviations in Table 1).
Water 10 01578 g001
Figure 2. Method for developing three category proxies for benchmark freshwater ecosystem service threats: (1) averaging of proxies within each category; (2) aggregation of proxies from the three categories, (i) external (regional hydro-climatic), (ii) internal (land use and watershed management-related), (iii) social (wealth and governance related); and (3) building benchmark-proxy relationships through regression analysis.
Figure 2. Method for developing three category proxies for benchmark freshwater ecosystem service threats: (1) averaging of proxies within each category; (2) aggregation of proxies from the three categories, (i) external (regional hydro-climatic), (ii) internal (land use and watershed management-related), (iii) social (wealth and governance related); and (3) building benchmark-proxy relationships through regression analysis.
Water 10 01578 g002
Figure 3. Estimated normalized FWES threat benchmark values and summary statistics for water supply (aHWS, a), biodiversity-related services (iBD, b), and nutrient regulation (N-load, c) for Ramsar (black bars) and SAFER sites (red bars). Inset graphs are frequency histograms for the same values (0.05 bins).
Figure 3. Estimated normalized FWES threat benchmark values and summary statistics for water supply (aHWS, a), biodiversity-related services (iBD, b), and nutrient regulation (N-load, c) for Ramsar (black bars) and SAFER sites (red bars). Inset graphs are frequency histograms for the same values (0.05 bins).
Water 10 01578 g003
Figure 4. Regression analysis for ecosystem service water supply benchmark aHWS (adjusted Human Water Stress; [14]) as a function of proxy (PaHWS) obtained using external, internal and socioeconomic indicators (Table 2): (a) regression results for GDPP as single proxy, and (b) best regression results for three-category proxy. Regression models (solid line; dashed lines 95% confidence interval) based on Ramsar sites (black circles), with SAFER sites (red circles) plotted for comparison.
Figure 4. Regression analysis for ecosystem service water supply benchmark aHWS (adjusted Human Water Stress; [14]) as a function of proxy (PaHWS) obtained using external, internal and socioeconomic indicators (Table 2): (a) regression results for GDPP as single proxy, and (b) best regression results for three-category proxy. Regression models (solid line; dashed lines 95% confidence interval) based on Ramsar sites (black circles), with SAFER sites (red circles) plotted for comparison.
Water 10 01578 g004aWater 10 01578 g004b
Figure 5. Strongest correlation for biodiversity-related threat benchmark, iBD (incident biodiversity threat; [14]), as a function of proxy (PiBD) obtained using external, internal and socioeconomic indicators (Table 2). Regression model (solid line; dashed lines 95% confidence interval) based on Ramsar sites (black circles), with SAFER sites (red circles) plotted for comparison.
Figure 5. Strongest correlation for biodiversity-related threat benchmark, iBD (incident biodiversity threat; [14]), as a function of proxy (PiBD) obtained using external, internal and socioeconomic indicators (Table 2). Regression model (solid line; dashed lines 95% confidence interval) based on Ramsar sites (black circles), with SAFER sites (red circles) plotted for comparison.
Water 10 01578 g005
Figure 6. Regression analysis for ecosystem services provided by ecosystem nutrient regulation service threat benchmark N-load (NL; [14]) as a function of proxy (PNL) based on external, internal and socioeconomic indicators (Table 2). Regression model (solid line; dashed lines 95% confidence interval) based on Ramsar sites (black circles), with SAFER sites (red circles) plotted for comparison.
Figure 6. Regression analysis for ecosystem services provided by ecosystem nutrient regulation service threat benchmark N-load (NL; [14]) as a function of proxy (PNL) based on external, internal and socioeconomic indicators (Table 2). Regression model (solid line; dashed lines 95% confidence interval) based on Ramsar sites (black circles), with SAFER sites (red circles) plotted for comparison.
Water 10 01578 g006
Figure 7. Three-category threat benchmark proxy regression outcomes for Ramsar sites clustered by national GDPP and access to advanced sanitation. Plots are for water supply (aHWS), biodiversity-related (iBD), and nutrient regulation (NL) threats (left to right) for higher, intermediate and lower development levels (top to bottom). (Significance levels: ** p < 0.01; * p < 0.05; NS = not significant; see Table A2 for regression results).
Figure 7. Three-category threat benchmark proxy regression outcomes for Ramsar sites clustered by national GDPP and access to advanced sanitation. Plots are for water supply (aHWS), biodiversity-related (iBD), and nutrient regulation (NL) threats (left to right) for higher, intermediate and lower development levels (top to bottom). (Significance levels: ** p < 0.01; * p < 0.05; NS = not significant; see Table A2 for regression results).
Water 10 01578 g007
Table 1. Current project (SAFER) study sites, priority freshwater ecosystem services, and major threats.
Table 1. Current project (SAFER) study sites, priority freshwater ecosystem services, and major threats.
Site Name, IDNationKey Ecosystem ServicesMajor Threats to Ecosystem Services
La Paloma Complex, LPCChileLivestock water supply
Hydropower
Recreation
Forestry (artisanal wood gathering)
Livestock, grazing
Aquaculture
Hydropower development
Recreational pressure
Invasive species
La Salada, LSArgentinaWater quality regulation
Biodiversity conservation
Recreation, education
Recreational pressure
Hydro-climate change
Land use/land cover change
Overfishing
Wastewater and irrigation drainage
Water diversion
Senguer River, SRArgentinaIrrigation water supply
Drinking water supply
Water quality regulation
Water diversions
Flood control engineering
Pollution (agriculture and oil exploration)
Dewatering (agriculture and oil exploration)
Hydro-climate change
Population increase (oil discovery)
Ciénaga Grande de Santa Marta, CGSMColombiaIrrigation water supply
Water quality regulation
Biodiversity conservation
Overfishing
Wastewater discharge
Agricultural drainage
Grazing
Sea level rise
Hydro-climate change
Diversions (agriculture)
Laguna de Rocha, LdRUruguayWater quality regulation
Recreation
Cultural/aesthetic
Recreational pressure (tourism)
Hydraulic engineering
Land use/land cover change (urbanization)
Storm surge
Flooding
Coastal erosion
San Joaquin River, SJRUSAIrrigation water supply
Water quality regulation
Biodiversity conservation
Population increase
Hydro-climate change
Diversions (agriculture)
Wastewater return flow and agricultural drainage
Land use/land cover change
Flood control engineering
Groundwater overdraft
Invasive species
Muskoka River Watershed, MRWCanadaWater quality regulation
Recreation
Cultural/aesthetic
Population increase
Invasive species
Land use/land cover change
Recreational pressure
Precipitation regime change (rain on frozen ground)
Temperature increase
Eutrophication
Table 2. Environmental and demographic indicators tested as proxies for aquatic ecosystem service threat benchmarks.
Table 2. Environmental and demographic indicators tested as proxies for aquatic ecosystem service threat benchmarks.
Proxy IDProxy NameDescription (0–1 = Low–High Threat Range Unless Otherwise Noted)References
External Proxies (Hydro-Climate, Water Demand)
WSVWater supply inter-annual variabilityWater supply variation from year to year[16]
SVWater supply seasonal variabilityWater supply variation from month to month [16]
HFOFlood occurrence Number of floods in recent history (1985–2011) [31]
DRODrought severityProduct of the average drought length and drought dryness soil dryness (1901–2008) [16,32]
Internal Proxies (Water Access, Watershed Management, Water Quality, Nutrient Management)
WRIWater return indexFraction of available water previously used and discharged upstream as wastewater effluent (0–1 = high–low threat for water supply ES; 0–1 = low–high threat for biodiversity-related and nutrient regulation ES)[16,33]
AGSUBAgricultural subsidiesDegree of environmental pressure exerted by subsidizing agricultural inputs (0–1 = high–low threat)[15]
NUENitrogen use efficiencyMeasure of the appropriate management of nitrogen resources for agricultural production (0–1 = high–low threat)[15]
NBalNitrogen use balance Measure of the appropriate management of nitrogen resources for agricultural production (0–1 = high–low threat)[15]
STORUpstream storageUpstream water storage capacity relative to total water supply (0–1 = high–low threat for water supply ES; 0–1 = low-high threat for biodiversity-related and nutrient regulation ES)[16]
ECO_SUpstream protected landFraction of total water supply that originates from protected watersheds (0–1 = high–low threat)[16]
WATSUPAccess to drinking waterFraction of nation’s population with access to improved drinking water (0–1 = high–low threat to water supply)[15]
ACSATAccess to sanitationFraction of a nation’s population with access to improved sanitation (0–1 = high–low threat to water supply)[15]
WWTWastewater treatedFraction of collected wastewater that is treated (0–1 = high–low threat)[15]
TPATerrestrial protected AreasDegree to which a nation achieves target of protecting 17% of its biomes (0–1 = high–low threat)[15]
Socioeconomic and Governance Proxies
RLRule of LawCaptures perceptions of the extent to which agents have confidence in and abide by rules of society (especially quality of contract enforcement, property rights, police, courts and likelihood of crime and violence) [34]
VAVoice & AccountabilityCaptures perceptions of the extent to which a country’s citizens are able to participate in selecting their government, as well as freedom of expression, freedom of association, and a free media[34]
GEGovernment EffectivenessCaptures perceptions of a nation’s quality of public services and civil services, the degree of its independence from political pressures, the quality of policy formulation and implementation, and the credibility of the government’s commitment to such policies[34]
GDPPGross Domestic Product per CapitaValue of annual goods produced and services provided by a nation divided by its population (0–1 = low–high normalized GDPP; threat scale tested in both directions)[35]

Share and Cite

MDPI and ACS Style

Harmon, T.C.; Smyth, R.L.; Chandra, S.; Conde, D.; Dhungel, R.; Escobar, J.; Hoyos, N.; Lozoya, J.P.; Nin, M.; Perillo, G.M.E.; et al. Socioeconomic and Environmental Proxies for Comparing Freshwater Ecosystem Service Threats across International Sites: A Diagnostic Approach. Water 2018, 10, 1578. https://doi.org/10.3390/w10111578

AMA Style

Harmon TC, Smyth RL, Chandra S, Conde D, Dhungel R, Escobar J, Hoyos N, Lozoya JP, Nin M, Perillo GME, et al. Socioeconomic and Environmental Proxies for Comparing Freshwater Ecosystem Service Threats across International Sites: A Diagnostic Approach. Water. 2018; 10(11):1578. https://doi.org/10.3390/w10111578

Chicago/Turabian Style

Harmon, Thomas C., Robyn L. Smyth, Sudeep Chandra, Daniel Conde, Ramesh Dhungel, Jaime Escobar, Natalia Hoyos, Juan Pablo Lozoya, Mariana Nin, Gerardo M.E. Perillo, and et al. 2018. "Socioeconomic and Environmental Proxies for Comparing Freshwater Ecosystem Service Threats across International Sites: A Diagnostic Approach" Water 10, no. 11: 1578. https://doi.org/10.3390/w10111578

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop