Quantifying coastal ecosystem trophic state at a macroscale using a Bayesian analytical framework.

One of the goals of coastal ecological research is to describe, quantify and predict human effects on coastal ecosystems. Broad cross-systems assessments to classify ecosystem status or condition have been developed, but are not updated frequently, likely because a lot of information and effort is needed to implement them. Such assessments could be more useful if the probability of being in a class indicating status or condition could be predicted using widely available data and information, providing a useful way to interpret changes in underlying predictors by considering their expected impact on ecosystem condition. To illustrate a possible approach, we used chlorophyll-a as an indicator of condition, in place of the intended comprehensive condition assessment. We demonstrated a predictive approach starting with a random forest model to inform variable selection, then used a Bayesian multilevel ordered categorical regression to quantify a coastal trophic state index and predict system status. We initially fit the model using non-informative priors to water quality data (total nitrogen and phosphorus, dissolved inorganic nitrogen and phosphorus, secchi depth) from 2010 and a regional factor. We then updated the model using prior distributions based on posterior parameter distributions from the initial fit and data from 2015. The Bayesian model demonstrates an intuitive way to update a model or analysis with new data while retaining the benefit of prior knowledge and maintaining flexibility to consider new kinds of information. To illustrate how the model could be used, we applied our developed trophic state index and classification to a time series of water quality data from Boston Harbor, a coastal ecosystem that has undergone significant changes in nutrient inputs. The analysis shows how water quality status and trends in Boston Harbor can be understood in the comparative ecological context provided by data from estuaries around the continental US and illustrates how the analytical approach could be used as an interpretive tool by non-practitioners of Bayesian statistics as well as a framework for further model development and analysis.

Rethinking the lake trophic state index.

Lake trophic state classifications provide information about the condition of lentic ecosystems and are indicative of both ecosystem services (e.g., clean water, recreational opportunities, and aesthetics) and disservices (e.g., cyanobacteria blooms). The current classification schemes have been criticized for developing indices that are single-variable based (vs. a complex aggregate of multi-variables), discrete (vs. a continuous), and/or deterministic (vs. an inherently random). We present an updated lake trophic classification model using a Bayesian multilevel ordered categorical regression. The model consists of a proportional odds logistic regression (POLR) that models ordered, categorical, lake trophic state using Secchi disk depth, elevation, nitrogen concentration (N), and phosphorus concentration (P). The overall accuracy, when compared to existing classifications of trophic state index (TSI), for the POLR model was 0.68 and the balanced accuracy ranged between 0.72 and 0.93. This work delivers an index that is multi-variable based, continuous, and classifies lakes in probabilistic terms. While our model addresses aforementioned limitations of the current approach to lake trophic classification, the addition of uncertainty quantification is important, because the trophic state response to predictors varies among lakes. Our model successfully addresses concerns with the current approach and performs well across trophic states in a large spatial extent.

A Random Forest approach to predict the spatial distribution of sediment pollution in an estuarine system.

Modeling the magnitude and distribution of sediment-bound pollutants in estuaries is often limited by incomplete knowledge of the site and inadequate sample density. To address these modeling limitations, a decision-support tool framework was conceived that predicts sediment contamination from the sub-estuary to broader estuary extent. For this study, a Random Forest (RF) model was implemented to predict the distribution of a model contaminant, triclosan (5-chloro-2-(2,4-dichlorophenoxy)phenol) (TCS), in Narragansett Bay, Rhode Island, USA. TCS is an unregulated contaminant used in many personal care products. The RF explanatory variables were associated with TCS transport and fate (proxies) and direct and indirect environmental entry. The continuous RF TCS concentration predictions were discretized into three levels of contamination (low, medium, and high) for three different quantile thresholds. The RF model explained 63% of the variance with a minimum number of variables. Total organic carbon (TOC) (transport and fate proxy) was a strong predictor of TCS contamination causing a mean squared error increase of 59% when compared to permutations of randomized values of TOC. Additionally, combined sewer overflow discharge (environmental entry) and sand (transport and fate proxy) were strong predictors. The discretization models identified a TCS area of greatest concern in the northern reach of Narragansett Bay (Providence River sub-estuary), which was validated with independent test samples. This decision-support tool performed well at the sub-estuary extent and provided the means to identify areas of concern and prioritize bay-wide sampling.

Associations between chlorophyll a and various microcystin health advisory concentrations.

Cyanobacteria harmful algal blooms (cHABs) are associated with a wide range of adverse health effects that stem mostly from the presence of cyanotoxins. To help protect against these impacts, several health advisory levels have been set for some toxins. In particular, one of the more common toxins, microcystin, has several advisory levels set for drinking water and recreational use. However, compared to other water quality measures, field measurements of microcystin are not commonly available due to cost and advanced understanding required to interpret results. Addressing these issues will take time and resources. Thus, there is utility in finding indicators of microcystin that are already widely available, can be estimated quickly and in situ, and used as a first defense against high levels of microcystin. Chlorophyll a is commonly measured, can be estimated in situ, and has been shown to be positively associated with microcystin. In this paper, we use this association to provide estimates of chlorophyll a concentrations that are indicative of a higher probability of exceeding select health advisory concentrations for microcystin. Using the 2007 National Lakes Assessment and a conditional probability approach, we identify chlorophyll a concentrations that are more likely than not to be associated with an exceedance of a microcystin health advisory level. We look at the recent US EPA health advisories for drinking water as well as the World Health Organization levels for drinking water and recreational use and identify a range of chlorophyll a thresholds. A 50% chance of exceeding one of the specific advisory microcystin concentrations of 0.3, 1, 1.6, and 2 µg/L is associated with chlorophyll a concentration thresholds of 23, 68, 84, and 104 µg/L, respectively. When managing for these various microcystin levels, exceeding these reported chlorophyll a concentrations should be a trigger for further testing and possible management action.

Developing a wintering waterfowl community baseline for environmental monitoring of Narragansett Bay, Rhode Island.

In 2004, the Atlantic Ecology Division of the US Environmental Protection Agency's Office of Research and Development began an annual winter waterfowl survey of Rhode Island's Narragansett Bay. Herein, we explore the survey data gathered from 2004 to 2011 in order to establish a benchmark understanding of our waterfowl communities and to establish a statistical framework for future environmental monitoring. The abundance and diversity of wintering waterfowl were relatively stable during the initial years of this survey, except in 2010 when there was a large spike in abundance and a reciprocal fall in diversity. There was no significant change in ranked abundance of most waterfowl species, with only Bufflehead ( Bucephala albeola) and Hooded Merganser ( Lophodytes cucllatus) showing a slight yet significant upward trend during the course of our survey period. Nonmetric multidimensional scaling (NMDS) was used to examine the community structure of wintering waterfowl. The results of the NMDS indicate that there is a spatial structure to the waterfowl communities of Narragansett Bay and this structure has remained relatively stable since the survey began. Our NMDS analysis helps to solidify what is known anecdotally about the bay's waterfowl ecology, and provides a formalized benchmark for long-term monitoring of Narragansett Bay's waterfowl communities. Birds, including waterfowl, are preferred bioindicators and we propose using our multivariate approach to monitor the future health of the bay. While this research focuses on a specific area of New England, these methods can be easily applied to novel areas of concern and provide a straightforward nonparametric approach to community-level monitoring. The methods provide a statistic test to examine potential drivers of community turnover and well-suited visualization tools.

Predicting and mapping potential Whooping Crane stopover habitat to guide site selection for wind energy projects.

Migratory stopover habitats are often not part of planning for conservation or new development projects. We identified potential stopover habitats within an avian migratory flyway and demonstrated how this information can guide the site-selection process for new development. We used the random forests modeling approach to map the distribution of predicted stopover habitat for the Whooping Crane (Grus americana), an endangered species whose migratory flyway overlaps with an area where wind energy development is expected to become increasingly important. We then used this information to identify areas for potential wind power development in a U.S. state within the flyway (Nebraska) that minimize conflicts between Whooping Crane stopover habitat and the development of clean, renewable energy sources. Up to 54% of our study area was predicted to be unsuitable as Whooping Crane stopover habitat and could be considered relatively low risk for conflicts between Whooping Cranes and wind energy development. We suggest that this type of analysis be incorporated into the habitat conservation planning process in areas where incidental take permits are being considered for Whooping Cranes or other species of concern. Field surveys should always be conducted prior to construction to verify model predictions and understand baseline conditions.

Enhanced migratory waterfowl distribution modeling by inclusion of depth to water table data.

In addition to being used as a tool for ecological understanding, management and conservation of migratory waterfowl rely heavily on distribution models; yet these models have poor accuracy when compared to models of other bird groups. The goal of this study is to offer methods to enhance our ability to accurately model the spatial distributions of six migratory waterfowl species. This goal is accomplished by creating models based on species-specific annual cycles and introducing a depth to water table (DWT) data set. The DWT data set, a wetland proxy, is a simulated long-term measure of the point either at or below the surface where climate and geological/topographic water fluxes balance. For species occurrences, the USGS' banding bird data for six relatively common species was used. Distribution models are constructed using Random Forest and MaxEnt. Random Forest classification of habitat and non-habitat provided a measure of DWT variable importance, which indicated that DWT is as important, and often more important, to model accuracy as temperature, precipitation, elevation, and an alternative wetland measure. MaxEnt models that included DWT in addition to traditional predictor variables had a considerable increase in classification accuracy. Also, MaxEnt models created with DWT often had higher accuracy when compared with models created with an alternative measure of wetland habitat. By comparing maps of predicted probability of occurrence and response curves, it is possible to explore how different species respond to water table depth and how a species responds in different seasons. The results of this analysis also illustrate that, as expected, all waterfowl species are tightly affiliated with shallow water table habitat. However, this study illustrates that the intensity of affiliation is not constant between seasons for a species, nor is it consistent between species.