A Basin-Wide Survey of Coastal Wetlands of the Laurentian Great Lakes: Development and Comparison of Water Quality Indices

Coastal wetlands of the Laurentian Great Lakes are vital habitats for biota of ecological and economic importance. These habitats are susceptible to water quality impairments driven by runoff from the landscape due to their location along the shoreline. Monitoring of the overall status of biotic and abiotic conditions of coastal wetlands within the Great Lakes has been ongoing for over a decade. Here, we utilize measurements of aquatic physicochemical and land cover variables from 877 vegetation zones in 511 coastal wetland sites spanning the US and Canadian shorelines of the entire Great Lakes basin. Our objective is to develop water quality indices based on physicochemical measures (Chem-Rank), land use/land cover (LULC-Rank), and their combined effects (Sum-Rank and Simplified Sum-Rank), for both vegetation zones and wetland sites. We found that water quality differed among wetland vegetation types and among Great Lakes regions, corroborating previous findings that human land use alters coastal wetland water quality. Future monitoring can use these straightforward, easy-to-calculate indices to assess the abiotic condition of aquatic habitats. Our data support the need for management efforts focused on reducing nutrient and pollution loads that stem from human activities, particularly in the developed southern portions of the Great Lakes basin.


Introduction
In recent decades, aquatic resource managers have increasingly used biological indicators to assess ecosystem health, in place of chemical and physical sampling of habitats, inferring that biota integrate conditions over time, whereas chemical and physical samples represent only a snapshot of the ecosystem. Published references to biological indicators have increased over the last 20 years, with many new biological indices developed for aquatic ecosystems (Ruaro and Gubiani 2013). Biological indices are developed based on the underlying assumption that the suite of anthropogenic impacts on an ecosystem can be best represented by the biological community capable of persisting within that ecosystem (Karr 1981). Ecosystems with taxa sensitive to physical and chemical stressors are assumed to have experienced fewer Electronic supplementary material The online version of this article (https://doi.org/10.1007/s13157-019-01198-z) contains supplementary material, which is available to authorized users.
* Anna M. Harrison harri25a@cmich.edu 1 anthropogenic impacts and thus be of higher quality. Although this index of biotic integrity (IBI) approach is wellestablished, it remains difficult to differentiate the relative influence of anthropogenic disturbances in biological communities from effects of natural disturbances (Suter 2001), although some watershed modeling approaches have assessed these relationships (Wiley et al. 2003). In coastal wetlands of the Laurentian Great Lakes (hereafter, Great Lakes), for example, anthropogenic disturbances can mask or interact with the effects of natural disturbances, such as wave exposure  or water-level fluctuations (Wilcox et al. 2002), on community structure. As a contrast to the IBI approach, physicochemical measurements in water (which include both chemical and physical variables) can be used to infer anthropogenic disturbances on an ecosystem more directly. Unfortunately, physicochemical measurements often fail to detect habitat alterations, and individual water chemistry or physical metrics may not reveal specific anthropogenic disturbances. Further, the effect of a disturbance on the entire ecosystem can persist longer than detectable changes in water quality. Although acute disturbances such as pulses of point-source contaminant discharge may not be detectable over long periods, chronic anthropogenic disturbances (e.g., increased nutrient and sediment loadings that accompany agricultural land use and urban development) have been associated with altered water quality across Great Lakes coastal wetlands (Chow-Fraser 2006; Morrice et al. 2008). These anthropogenic disturbances also influence the scores of macroinvertebrate-and fish-based IBIs within Great Lakes coastal wetlands (Uzarski et al. , 2005Cooper et al. 2018).
In addition to physicochemical measurements, the land use surrounding aquatic ecosystems often provides an indicator of anthropogenic stress, and multiple indices of land use stress within the Great Lakes basin have been developed. For example, Danz et al. (2007) developed a shoreline stress index for the U.S. shoreline of the Great Lakes by combining five classes of stressors (i.e., agriculture, atmospheric deposition, human population, land cover, and point-source pollution). Host et al. (2011) developed a stress index for the entirety of Lake Superior using land use/land cover, population density, point-source discharge, and road density. This stress index was spatially expanded to create a watershed stress index for the entire Great Lakes Basin, as part of the State of the Lakes Ecosystem Conference . Further research assessed the influences of these watershed-based stressors on established responses of various species assemblages in coastal areas of the Great Lakes (Host et al. 2019).
Efforts to combine landscape variables with water chemistry and physical variables typically focus on a driver-response relationship (e.g., Chow-Fraser 2006), but these two types of indices are rarely combined in a monitoring framework. Combined water quality and land use indices can be useful in understanding habitat quality in Great Lakes coastal wetlands because these two approaches allow for direct estimation of current in situ water quality (physical and chemical conditions in water) and inference of the longer term surrogate measure of chronic pollution (land use). Water quality data are also easy to collect, can require less specialized collection expertise and equipment than an IBI, and are often available online through a variety of monitoring programs (e.g., governmental agencies, non-profit organizations, or universities).
The habitat condition and water quality of Great Lakes coastal wetlands can influence the capacity of coastal wetlands to provide beneficial functions. Coastal wetlands are biodiverse ecosystems that serve as a major component of nearshore lake nutrient cycling. They also filter pollutants and other contaminants running off the landscape that would otherwise imperil the pelagic biota of the Great Lakes (Uzarski et al. 2017 and references therein). Despite their ecological importance, nearly half of all coastal wetland areas in the Great Lakes have been converted to other land uses (e.g., agriculture, development) since European settlement, and many remaining wetlands are threatened by multiple anthropogenic stressors (Maynard and Wilcox 1997; Environment Canada and the U.S. Environmental Protection Agency 2014).
Since 2011, a basin-wide monitoring effort has been underway to quantify the overall ecological status of coastal wetlands in the Laurentian Great Lakes using standardized, quantitative methods. This monitoring program has provided data on the abiotic and biotic conditions of coastal wetlands to a range of stakeholders (Uzarski et al. 2017). Our primary objective in this study is to use physicochemical data collected from the coastal wetland monitoring project (see Uzarski et al. 2017), coupled with landscape metrics, to develop an overall view of coastal wetland water quality along U.S. and Canadian shorelines throughout the entire Great Lakes basin. We previously proposed an index of anthropogenic disturbance combining physicochemical and landscape stressors (Uzarski et al. 2005), but this previous index was limited in its spatial scale. Taken individually, specific indices of water quality or landscape stress have the potential to miss or overlook either acute or chronic signatures of anthropogenic impairments. Thus, we predict that the combination of landscape-level and physicochemical metrics is likely to capture a wider range of anthropogenic disturbances than consideration of either stressor type individually.
We also anticipate spatial variability within each indicator metric because water quality will vary both locally among vegetation types and regionally across the Great Lakes basin. Wetland vegetation taxa have variable tolerances for anthropogenic disturbances, including water quality (Lougheed et al. 2001), to the extent that vegetation communities can be used as indicators of wetland health (Albert and Minc 2004;Croft and Chow-Fraser 2007;Grabas et al. 2012). Specifically, we expect to see lower water quality scores in Phrgamites and Typha vegetation zones (Albert and Minc 2004;Zedler and Kercher 2004) compared to vegetation zones dominated by native vegetation, especially those with lower tolerances to anthropogenic disturbance (e.g., Schoenoplectus) (Croft and Chow-Fraser 2007).

Site Selection and Description
We performed field sampling between mid-June and early September from 2011 through 2015 as part of a coordinated, basin-wide, coastal wetland monitoring program (greatlakeswetlands.org). We used a stratified-random site-selection procedure to select a subset of all accessible coastal wetlands identified throughout the Great Lakes to sample each year (see Uzarski et al. 2017 for a detailed description of site identification and sampling approach). We measured a suite of physicochemical (physical and chemical) and biological variables at each site using a variety of sampling methods (Uzarski et al. 2017 and references therein). Here, we focus on aquatic physicochemical variables and landscape metrics to understand the drivers of overall water quality in Great Lakes coastal wetlands.
Per the program-wide sampling protocols, at each site we visually identified all mono-dominant vegetation habitat types consisting of at least 75% of a single category of wetland vegetation. Each site could have multiple mono-dominant habitat types, which we refer to as vegetation zones. Emergent vegetation zones were typically mono-dominant, whereas SAV and wet meadow were generally mixtures of different taxa. Submerged aquatic vegetation (SAV) zones were established where SAV was present, but there was not any floating or emergent vegetation (Uzarski et al. 2017). Over the 5-year period, we sampled 1114 vegetation zones from 518 unique sites, but for this analysis, we only include vegetation types that accounted for >1% of all vegetation zones sampled (i.e., we excluded extremely rare vegetation types). Ten vegetation types were included in this analysis: wet meadow (dominated by sedges; Carex, Juncus, and Eleocharis), Typha, Phragmites, mixed emergent (no dominant emergent vegetation type), Peltandra/Sagittaria/ Pontederia (PSP, combined due to similar growth form), Sparganium, Zizania (wild rice), Schoenoplectus (bulrush), Nuphar-Nymphaea (lily), and submerged aquatic vegetation (SAV) ( Table 1). Where data were missing from a location (i.e., the 12% of zones missing one or two physicochemical measurements from the rank indices), we substituted the mean value of the variable estimated from values measured within all other zones. This procedure allowed us to maintain the same sample size for the rank calculations without influencing the final rank of the zone or site. Ultimately, we included 877 unique vegetation zones from 511 unique sites (Table 1).

Landscape Characteristics
We used the North American Land Change Monitoring System (NALCMS) 2010 land cover of North America at 30 m to provide land cover data for each site (CEC 2015). To assess both local and regional land cover, we extracted land cover estimates from 1-km and 20-km buffers around each wetland site. Distance buffers allowed for a standardized land cover estimate across different wetland classes (lacustrine, barrier-protected, and riverine). We extracted four land cover categories from the NALCMS classes, including %-natural vegetation (forest + grassland + shrub), %-developed (urban), %-agricultural (cropland), and %-wetland.

Aquatic Physicochemical Characteristics
At each vegetation zone, we measured a suite of physical and chemical variables in situ at the middle of the water column. Using YSI (or similar) multi-parameter sondes, we measured in situ water temperature, pH, specific conductance (μS cm −1 ), and dissolved oxygen (DO concentration; mg L −1 and % saturation). We calibrated data sondes prior to sampling according to manufacturer protocols and a stringent quality assurance project plan that is followed by all field staff for the coastal wetland monitoring program (Uzarski et al. 2017;greatlakeswetlands.org). We also collected a composite water sample (of 3 sub-samples) within each vegetation zone in an acid-washed carboy and stored samples on ice until returning from the field. We measured water clarity in the field using a sub-sample of this water sample with a 100-cm transparency tube (Anderson and Davie 2004).
In the laboratory, we analyzed the water sample from each vegetation zone for a range of analyses using standard methods. We transferred two unfiltered water samples into 250-mL, acid-washed polypropylene bottles and stored these samples frozen for later analysis of total nitrogen (TN; mg L −1 ) and total phosphorus (TP; mg L −1 ) using standard methods (APHA 2005). An additional subsample was filtered through a Whatman GF/C glass fiber filter and frozen for later analysis. We then thawed and analyzed the glass fiber filter for chlorophyll a (chl a; μg L −1 ) using standard methods (APHA 2005). We further filtered the filtrate from the chl a sample through an acid-and deionized water (DI)-rinsed 0.45-μm Millipore membrane filter into an acid-washed, DI-rinsed polypropylene bottle and subsequently froze the sample. We later thawed and analyzed this sample for soluble reactive phosphorus (SRP, mg L −1 ), ammonium-N (NH 4 + -N, mg L −1 ), and [nitrate + nitrite]-N (NO 3 − -N, mg L −1 ) via standard methods (APHA 2005).

Index Development
Four water quality indices were calculated and compared in this study: (1) Chem-Rank Index, an index based only on physicochemical variables, (2) LULC-Rank Index, an index based only on land use/land cover data, (3) Simplified-Sum Rank Index (SS-Rank), a simplified index comprised of both physicochemical and land cover variables, and (4) the original Sum Rank Index (Sum-Rank), which includes a more complex combination of physicochemical and land cover variables than SS-Rank (Uzarski et al. 2005;Uzarski et al. 2017). The Chem-Rank index and LULC-Rank index highlight the relative contributions of each to the Sum-Rank and SS-Rank indices.
Because many of the physicochemical and land cover variables were non-normally distributed and expressed on different scales, we rank-transformed each variable prior to index development. Ranks were assigned to variables such that rank order was proportionate to the inferred degree of anthropogenic disturbance. Variables with high values indicative of low disturbance (i.e., water clarity, % natural vegetation and % wetland land cover) were ranked ascendingly, with high ranks given to high values. For example, at the site level, a water clarity rank of 511 was given to the site with the clearest water, and a rank of 1 was given to the site with the lowest water clarity. High values of variables indicative of high disturbance (i.e., specific conductance, NO 3 − -N, NH 4 + -N, SRP, TN, TP, chl a, % developed land cover, % agricultural land cover) were scaled such that higher values received lower ranks. For example, a SRP rank of 1 was given to the site with the highest SRP concentration (or highest disturbance), whereas a rank of 511 was given to the site with the lowest SRP concentration. We acknowledge that nutrients occur naturally in these systems (i.e., without human inputs); however, the majority of nitrogen and phosphorous inputs come from anthropogenic sources.
Rank values for all the variables at both the vegetation zone level and site level were then summed for index calculation. We scaled each summed index value from 0 to 100, using Eq. 1, where 'Rank i ' is the summed rank value for each site or vegetation zone and 'Rank min ' and 'Rank max ' are the corresponding lowest and highest ranked site or zone within each index. Separate scales were derived for vegetation zone (based on 877 observations) and for site (based on 511 observations).
Chem-Rank was calculated by summing individual ranks for eight physicochemical variables (specific conductance, NO 3 − -N, NH 4 + -N, SRP, TN, TP, chl a, and water clarity). The Chem-Rank index was measured at both the vegetation zone and site level. Although chloride is an important indicator of anthropogenic influences in water quality, it was not included in the indices because it was not measured at all sites in our dataset. Values of specific conductance and chloride were strongly positively correlated within vegetation zones where both variables were measured (r = 0.82). Specific conductance serves as a valuable surrogate for anthropogenic inputs of chloride, and thus is a useful metric of anthropogenic The ten regions were defined by first by country and ecoregion (Southern U.S., Northern U.S., and Canada) and then divided by lake: two regions in Lake Erie (LE North and South), three regions in Lake Huron (LH Northeast, Southwest and West), two regions in Lake Michigan (LM North and South), two regions in Lake Ontario (LO North and South), and the entirety of Lake Superior (LS) disturbance in our indices. We included both total (TN, TP) and dissolved (NO 3 − -N, NH 4 + -N, SRP) nutrients in our water quality indices because each form of inorganic (or organic) nutrient can affect an ecosystem differently. For nitrogen, NH 4 + -N is typically more biologically available than NO 3 − -N, whereas NO 3 − -N represents the majority of dissolved inorganic nitrogen and can readily alleviate N limitation. For phosphorous, SRP is bioavailable, whereas TP is tied up in biological material and can correlate with suspended sediments (Grayson et al. 1996;Baustian et al. 2018).
We calculated a land use/land cover only index, LULC-Rank, to assess the contribution of only land cover to potential effects on water quality. The LULC-Rank was calculated by summing the ranks of all eight LULC variables calculated from the NALCMS (1-km and 20-km buffers for % developed, % natural vegetation, % agriculture, % wetland). We calculated LULC-Rank at the site level only because different vegetation zones at the same site have identical LULC values.
The SS-Rank index was developed as a combination of in situ physicochemical variables and land cover. It is a simplified index, developed from the original and more complex Sum-Rank (described below). SS-Rank summed the ranks of the eight physicochemical variables in Chem-Rank with the ranks of the eight land cover ranks from LULC-Rank, giving each index equal weight in SS-Rank. SS-Rank was calculated at both the vegetation zone and site levels.
Our original water quality index, Sum-Rank, contained all the same variables as the SS-Rank index described above, but included two additional categories of variables. The first category consisted of variables ranked based on distance from the median (i.e., more extreme values were indicative of disturbance). This extreme-value approach was used to rank NO 3 − -N, DO, and pH in Sum-Rank, but was removed from SS-Rank to simplify the index calculations and because the median nutrient concentrations from on our dataset are generally higher than what is expected from a natural system without anthropogenic sources. The second additional category in Sum-Rank was a multivariate ordination, Principal Component Analysis (PCA), calculated from raw values of all physicochemical variables and land cover proportions, and then scaled so that higher PC1 values indicated higher ranking (Uzarski et al. 2005). The PCA was also omitted from SS-Rank, to minimize calculation steps required for the simplified index.

Analyses
Indices developed at the vegetation zone level were calculated from values measured within each individual vegetation zone at each site (n = 877). Vegetation type averages for water quality are available in Online Resource 1. We used linear mixed effects models to test for statistically significant differences in water quality across vegetation type (Table 1) and included region as a random effect for each model to account for regional spatial variability of vegetation types. To define the ten regions, the Great Lakes basin was divided into three broad regions by country and ecoregion (Bailey and Cushwa 1981;Omernik 1987): Southern U.S., Northern U.S., and Canada. Each region was then divided by lake, with the entirety of Lake Superior as one region because there were so few wetlands on the northern side of Lake Superior. Vegetation type group means were compared using a Tukey post-hoc multiple comparisons test. Differences among vegetation types were calculated for only three indices (SS-Rank, Sum-Rank, and Chem-Rank) because LULC was only calculated at the site level. The lme4 package (Bates et al. 2015) was used for linear mixed effects models and the emmeans package (Lenth 2018) for subsequent Tukey multiple comparison tests. All statistical analyses and figure preparations were performed in R (version 3.5.0; R Development Core Team 2018).
For site-level indices, we averaged each variable's raw values for all vegetation zones sampled within an individual site to calculate site-level values and ranks (n = 511). Regional average values for water quality variables and land cover are available in Online Resources 2 and 3. We then used analysis of variance (ANOVA) to test whether water quality differed across Great Lakes regions (Table 1) using the site-level SS-Rank, Sum-Rank, Chem-Rank, and LULC-Rank indices and compared groups using a Tukey post-hoc test.
At the site level, we also quantified the strength of relationships between the six pairs of four indices (SS-Rank, Sum-Rank, Chem-Rank, and LULC-Rank) using Pearson's correlation. We then compared individual variables included in Chem-Rank or LULC-Rank with the alternate index using Pearson's correlation with a Bonferroni correction for multiple comparisons. For example, we quantified the correlation between water clarity and LULC-Rank, as well as the correlation between developed land cover in the 20-km buffer with Chem-Rank. We did not compare individual variables with their associated index (e.g., NO 3 − -N~Chem-Rank) due to autocorrelation. For correlation analyses, we log 10 -transformed all physicochemical variables that fit a log-normal distribution (all variables except clarity), and we arcsinesquare root transformed land cover variables due to nonnormal distributions.
All four indices varied at the site level across Great Lakes regions. Each index had the lowest average value in northern Lake Erie (LE North) and the highest average in western Lake Huron (LH West), with the exception of the LULC-Rank, w h i c h w a s h i g h e s t i n L a k e S u p e r i o r ( F i g . 2 , Online Resource 6). SS-Rank values averaged from 20.2 (± 17.5 SD) to 75.8 (± 10.9 SD) and varied among regions (ANOVA; F 9, 501 = 72.59; p < 0.001; Fig. 2a). Sum-Rank values ranged, on average, from 22.5 (± 16.9 SD) to 74.5 (± 10.5 SD) across sites and differed among regions (ANOVA; F 9, 501 = 71.66; p < 0.001; Fig. 2b). Chem-Rank averages ranged across Great Lakes regions from 28.7 (± 22.6 SD) to 72.8 (± 12.0 SD), and values varied among regions (ANOVA; F 9, 501 = 29.55; p < 0.001; Fig. 2c). Regional averages for the land cover index (LULC-Rank) ranged from 14.6 (± 13.1 SD) in LE North to 74.1 (± 15.0 SD) in Lake Superior, and index values differed among regions (ANOVA; F 9, 501 = 98.35; p < 0.001; Fig. 2d). Overall trends across all rank indices show the lowest water quality and land cover quality at sites in Lake Erie, southern Lake Michigan, and southwestern Lake Huron. Lake Superior, northern Lake Michigan, and western and northeastern Lake Huron had the highest ranked scores for both water quality and land cover across all four indices (Fig. 2, Fig. 3, Online Resources 6, 7). Chem-Rank had the greatest variability among regions compared to the other three rank indices (Fig. 2, Online Resource 6, 7).
Chem-Rank was correlated with seven of the eight landcover variables included in LULC-Rank, had positive correlations with natural vegetation and 20 km wetland land cover, and had negative correlations with both developed and agricultural land cover (Table 2). Of the LULC variables, Chem-Rank was most negatively correlated with the 20-km agricultural land cover (r = −0.51) and most positively correlated with 20-km natural vegetation (r = 0.55). LULC-Rank was negatively correlated with all the physicochemical variables included in Chem-Rank, except water clarity with which it was positively correlated ( Table 2). The strongest LULC-Rank correlation with a physicochemical variable was with specific conductance (r = −0.64).
Unsurprisingly, the four indices were correlated with each other (Fig. 4). The SS-Rank and Sum-Rank were highly correlated with one another (r = 0.98; Fig. 4a), and SS-Rank had slightly higher correlations with Chem-Rank (r = 0.88; Fig.  4c) and LULC-Rank (r = 0.90; Fig. 4e) compared with Sum-Rank correlations with Chem-Rank (r = 0.85; Fig. 4d) and LULC-Rank (r = 0.89; Fig. 4f). These results suggest that SS-Rank may be a suitable and easier to calculate and interpret alternative to the Sum-Rank.
There was a fair amount of variability in the relationship between Chem-Rank and LULC-Rank (r = 0.58), and the variability in Chem-Rank values was greater at low LULC-Rank values (< 40) compared to higher LULC-Rank values (Fig.  4b). This variability shows that these two indices are not explaining identical water quality impairments. It also demonstrates that, although land use/land cover represents the potential for stress, it is not necessarily indicative of actual stressors in the ecosystems. In addition, the variability may be illustrating differences in temporal and spatial stresses, where the water quality (one point in time) does not match the land cover (large area, relatively unchanging in the long-term). Therefore, the development of SS-Rank and Sum-Rank from Chem-Rank and LULC-Rank provides additional information that may not be evident from only one of the individual indices.

Basin-Wide Water Quality Patterns
Water quality impairment has been identified as an important factor contributing to wetland degradation (Morrice et al. 2008;Cvetkovic and Chow-Fraser 2011). There is interest among natural resource managers in protecting and restoring coastal wetlands in the Great Lakes, including the chemical and physical nature of these systems (Mitsch and Wang 2000), with hopes that improving coastal wetland water quality may both enhance habitat for biota and reduce pollution loads into open water areas of the adjacent lakes. Our analyses of water quality in 877 vegetation zones from 511 unique coastal wetlands represents the largest inventory of water quality monitoring data from Great Lakes coastal wetlands. Our results confirm previous findings of poor water quality in coastal wetlands in portions of the more heavily-populated southern Great Lakes basin, with generally higher water quality in the northern portions of the basin (Chow-Fraser 2006;Cvetkovic and Chow-Fraser 2011;Uzarski et al. 2017).
The spatial distribution of vegetation types within coastal wetlands (i.e., within sites) also supports the trend of higher water quality in the northern Great Lakes. The vegetation types with highest index scores (bulrush and wet meadow) were found predominantly (82% and 99%, respectively) in the four northern regions with the highest SS-Rank index scores (northeastern and western Lake Huron, northern Lake Michigan, and Lake Superior). It is unclear, however, if these factors are mechanistically related or if they reflect the influence of another extrinsic driver not included in the present analysis. Bulrush occur throughout the Great Lakes basin; thus, this trend likely reflects lower anthropogenic disturbance (water quality impairments among other potential disturbances) in the northern regions compared with southern regions. Bulrush zones may also be prevalent in higher quality northern sites because they are outside the area of Phragmites invasion at northern latitudes and because bulrushes tolerate the low nutrient concentrations that inhibit both Phragmites and Typha dominance (Albert and Minc 2004;Lemein et al. 2017). This latitudinal trend may explain the lower overall variability in rank values within the bulrush zones compared to the high variability in the five vegetation zones with the lowest average rank values: Lily, PSP, Phragmites, SAV, and Typha. This indicates that while vegetation zone presence is related to region, these relationships are likely a product of the plant community response to the amount of anthropogenic disturbance at a site, which generally trends along a latitudinal gradient in the Great Lakes.
Despite this general latitudinal trend, water quality in coastal wetlands near anthropogenically impacted areas in the northern basin, such as Duluth, MN, was also poor (Fig. 3). This pattern of poor water quality in wetlands located closer to human activity is further supported by the negative correlations found between Chem-Rank and agricultural and developed land cover using both 1-km and 20-km buffers ( Table 2). The 20-km land cover buffers have stronger correlations with the water quality endpoints than the 1-km buffers. This is possibly a product of the land use directly adjacent to the coast (1-km buffer) being non-representative of the surrounding land use, and the 20-km buffer likely better represents what is loading into the wetlands. For example, in Saginaw Bay, some coastal wetlands have narrow buffers of forested land use around the wetland, but the majority of the land use in the 20-km buffer is agriculture. This influence of agricultural and urban land use on nutrient loading from tributaries has been modeled at broad scales in the Great Lakes basin (Bosch et al. 2011;Robertson and Saad 2011), and empirical studies have consistently shown a negative relationship between water quality and land use (Chow-Fraser 2006;Trebitz et al. 2007;Morrice et al. 2008). Our study provides additional, complementary indices of water quality to assess Great Lakes coastal wetlands based on a large sample of sites along both the U.S. and Canadian shorelines. Most of the physicochemical variables measured correlated negatively with LULC-Rank values (all except water clarity), indicating that wetlands surrounded by higher proportions of agricultural and developed land cover had lower water quality. Specific conductance was most strongly (negatively) correlated with LULC-Rank values, whereas various dissolved nutrient metrics (i.e., NO 3 − , NH 4 + , SRP) were less correlated with LULC-Rank values (Table 2). Although N and P concentrations were high in many wetlands, coastal wetland algal biofilms are frequently nutrient-limited, and N-limitation appears to be more common than P-limitation (Cooper et al. 2016). If N-limitation is common across coastal wetlands, then internal cycling may explain the fairly weak relationship between dissolved inorganic N and LULC-Rank values. Dissolved inorganic nutrients are highly labile, and internal cycling can significantly alter dissolved nutrient concentrations in coastal wetlands (Krieger 2003;Morrice et al. 2004). Total nutrient concentrations (TN or TP) showed stronger negative relationships with LULC-Rank values than dissolved inorganic forms, further reflecting the role of internal biogeochemical cycling of nutrients within wetlands.

Strengths of the Index Approach
Assessments of ecological condition of aquatic habitats can address needs that may range from individual management questions related to specific drivers of environmental quality to broad-scale surveys intended to provide a general overview of a range of environmental conditions. Our goal in deriving the SS-Rank index was to create a simple metric to characterize water quality, based on easily measured variables, that is intuitive to formulate and easy to understand and interpret. Previously developed environmental indices tend to rely primarily on either water quality (Chow-Fraser 2006) or landscape data (Danz et al. 2007). Our results highlight the explanatory value of incorporating water quality monitoring with land use/land cover data for water quality index development. Although we found significant correlations between our physicochemical-and landscape-based sub-indices, there was considerable residual variation in the relationship (Fig.  4b), suggesting that combining both factors provides more information than either one individually.
The combination of both water quality variables and land cover allowed our indices to distinguish both regional (e.g., southern vs. northern Green Bay, eastern vs. western Lake Ontario) and basin-wide (e.g., northern vs. southern latitudes) gradients in wetland water quality (Fig. 3) and to identify areas where land cover is less predictive of water quality. For example, the average Chem-Rank index score was significantly higher in northern Lake Ontario (59.5 ± 17.7 SD) than in four other regions (northern and southern Lake Erie, southern Lake Michigan, and southern Lake Ontario), whereas the LULC-Rank index score for northern Lake Ontario was much lower (27.9 ± 9.0 SD) than the Chem-Rank value, and the LULC-Rank was not significantly different from the lowest ranked regions (Fig. 2). This suggests that despite northern Lake High index scores on the y-axis indicate high quality sites, compared to low y-axis values that indicate sites with poorer water quality. Letters denote differences among regions from Tukey HSD post-hoc tests. Regions were broken into ten geographic areas (labeled on the x-axis), described in Table 1. Values for means and standard deviations of ranks by region are available in Online Resource 6 Ontario having a higher proportion of developed and agricultural land cover, the effects of these land uses are not reflected as strongly in physicochemical variables (i.e., Chem-Rank) compared with other regions. Conversely, the Duluth area of western Lake Superior has higher LULC-Rank values than Chem-Rank values, which also indicates that land use/land cover alone is not necessarily a good predictor of water quality (Fig. 3c, d). This area has high shipping traffic and associated anthropogenic disturbances that influence water quality yet are not accounted for in the land cover buffers. It is also possible that not all the upstream land use influencing these wetlands is accounted for in the 20-km land use buffer in this particular area.
An additional strength of our SS-Rank approach is its relative simplicity. In contrast to previous indices put forth by our research group (e.g., Uzarski et al. 2005; the original Sum-Rank index described in this paper), the SS-Rank index is based solely on the rank-order of variables and not multivariate ordinations or median-centered values. Although multivariate approaches may capture patterns in environmental gradients that univariate approaches do not, the use of ordinations can complicate the development and interpretation of an index and make it difficult to establish which individual factors drive patterns within a given index. The removal of mediancentered values also simplifies the SS-Rank index, compared to the Sum-Rank index. This was justified after univariate regressions found that relationships with physicochemical variables were not improved by transforming their values based on distance from the median. Dissolved oxygen and pH can also vary throughout the day (Nelson et al. 2009), which further justifies their removal in SS-Rank index. Therefore, to maximize the simplicity of the SS-Rank index, we ranked absolute values of NO 3 − -N (rather than their distance from the median) and removed DO and pH from the SS-Rank. Fig. 3 Spatial distribution of index scores were mapped across the entire Great Lakes basin (SS-Rank; a, Sum-Rank; b, Chem-Rank; c and LULC-Rank; d). Each point represents a unique coastal wetland site. Rank values were equally divided into five sub-classes to represent water quality status. Water quality ranges from light yellow-colored points representing the poorest water quality (low numbers) to dark bluecolored points representing the highest water quality (high numbers) Our simplified indices are still consistent with values of the ordination-based Water Quality Index (WQI) proposed by Chow-Fraser (2006), as we found strong correlations between the WQI and both SS-Rank and Chem-Rank using our sitelevel data (r = 0.81 and r = 0.86, respectively). The Chem-Rank index had a higher correlation than SS-Rank with the WQI, likely because the WQI is calculated from water quality variables only, and SS-Rank also incorporates land cover. Our SS-Rank builds upon the WQI by including land-cover proportions in the index calculation in addition to physicochemical variables.

Applications of Rank Indices
As a broadly adaptable tool, the SS-Rank index could be applied to assess habitat quality in other systems and regions, rather than just coastal wetlands in the Laurentian Great Lakes. Because our index is based on univariate rankings and then scaled from zero to 100, a subset of variables could be used to calculate the index rather than requiring measurement of all physicochemical variables included in the full index. To assess the quality of freshwater coastal wetland ecosystems using this index, we recommend that, at minimum, specific conductance, TN, TP, and chlorophyll a (measured via standard methods, APHA 2005) be included, along with all land cover data. These variables are all strongly associated with anthropogenic disturbance and have previously been shown to drive water quality of coastal wetlands (Crosbie and Chow-Fraser 1999;Trebitz et al. 2007;Morrice et al. 2008). Although chloride was not included in the rank indices in this paper, we recommend its inclusion, if possible, but suggest that specific conductance is also an excellent metric for anthropogenic disturbances, as it is easy to measure and correlates with multiple disturbance metrics.
In combining physicochemical and land cover variables across a large spatial scale, the SS-Rank index also provides government agencies in Canada and the U.S. with an updated tool to assess progress towards achieving the objectives agreed upon in the Great Lakes Water Quality Agreement (GLWQA), "to restore and maintain the chemical, physical, and biological integrity of the waters of the Great Lakes Basin Ecosystem" (United States and Canada 1978). The SS-Rank index could be used as a tool to meet obligations outlined in the GLWQA annexes, as each annex contains specific requirements, and the SS-Rank could be applied to the goals of several annexes (e.g., Annex 2, Lakewide Management; Annex 7, Habitat and Species; and Annex 10, Science). For example, Annex 10, Science calls for signatories to establish, maintain, and update comprehensive, science-based ecosystem indicators to assess the state of the Great Lakes and measure progress toward achievement of the objectives of the GLWQA. The SS-Rank index could provide managers with an indicator of coastal wetland health that could be used to assess changes over time in wetland health throughout the region.

Conclusions
The Laurentian Great Lakes basin spans a vast geographic area, including a range of climates, geology, and land use. The densest human land use (agricultural + developed) generally occurs in the southern portion of the basin, but all five lakes have areas that have been substantially altered by human activity. In our survey, which encompasses the largest number of Great Lakes coastal wetlands surveyed to date and is comprised of wetlands along both the U.S. and Canadian shorelines, we revealed patterns of coastal wetland water quality associated with land use and land cover alterations that complement previous studies. Such water quality impairments have been associated with changes to biological communities inhabiting Great Lakes coastal wetlands, including macrophytes (Lougheed et al. 2001;Croft and Chow-Fraser 2007), macroinvertebrates Cooper et al. 2007), and fishes (Uzarski et al. 2005;Seilheimer and Chow-Fraser 2006;Cooper et al. 2018). A focus on improving water quality in coastal wetlands is vital to enhance habitat for biota and reduce pollution loads to open-water habitats of the Great Lakes. Our results, together with previous studies, indicate that reducing nutrient loading from human activities on the landscape (e.g., agriculture and The p values reflect Bonferroni correction for multiple comparisons a Log 10 transformed b Arcsine-square root transformed urbanization) is critical for improving water quality. Best management practices designed to minimize nutrient application rates (e.g., soil testing) and retain nutrients and other pollutants on the land prior to downstream export (e.g., riparian buffers) are essential to protecting and improving abiotic and biotic conditions of coastal wetlands, as well as the overall ecological functioning of the Great Lakes.