Shifting institutional culture to develop climate solutions with Open Science.

To address our climate emergency, "we must rapidly, radically reshape society"-Johnson & Wilkinson, All We Can Save. In science, reshaping requires formidable technical (cloud, coding, reproducibility) and cultural shifts (mindsets, hybrid collaboration, inclusion). We are a group of cross-government and academic scientists that are exploring better ways of working and not being too entrenched in our bureaucracies to do better science, support colleagues, and change the culture at our organizations. We share much-needed success stories and action for what we can all do to reshape science as part of the Open Science movement and 2023 Year of Open Science.

nsink: An R package for flow path nitrogen removal estimation.

The nsink package estimates cumulative nitrogen (N) removal along a specified flow path and is based on methodologies outlined in Kellogg et al. (2010). For a user-specified watershed (i.e., hydrologic unit code (HUC)), nsink downloads all required datasets from public datasets in the United States, prepares data for use, summarizes N removal along a flow path and creates several static maps. The results of an nsink analysis may be exported to standard geospatial files for use in other applications.

Socioeconomic factors affecting local support for black bear recovery strategies.

There is global interest in recovering locally extirpated carnivore species. Successful efforts to recover Louisiana black bear in Louisiana have prompted interest in recovery throughout the species' historical range. We evaluated support for three potential black bear recovery strategies prior to public release of a black bear conservation and management plan for eastern Texas, United States. Data were collected from 1,006 residents living in proximity to potential recovery locations, particularly Big Thicket National Preserve. In addition to traditional logistic regression analysis, we used conditional probability analysis to statistically and visually evaluate probabilities of public support for potential black bear recovery strategies based on socioeconomic characteristics. Allowing black bears to repopulate the region on their own (i.e., without active reintroduction) was the recovery strategy with the greatest probability of acceptance. Recovery strategy acceptance was influenced by many socioeconomic factors. Older and long-time local residents were most likely to want to exclude black bears from the area. Concern about the problems that black bears may cause was the only variable significantly related to support or non-support across all strategies. Lack of personal knowledge about black bears was the most frequent reason for uncertainty about preferred strategy. In order to reduce local uncertainty about possible recovery strategies, we suggest that wildlife managers focus outreach efforts on providing local residents with general information about black bears, as well as information pertinent to minimizing the potential for human-black bear conflict.

A process for comparing and interpreting differences in two benthic indices in New York Harbor.

Often when various estuarine benthic indices disagree in their assessments of benthic condition, they are reflecting different aspects of benthic condition. We describe a process to screen indices for associations and, after identifying candidate metrics, evaluate metrics individually against the indices. We utilize radar plots as a multi-metric visualization tool, and conditional probability plots and receiver operating characteristic curves to evaluate associations seen in the plots. We investigated differences in two indices, the US EPA Environmental Monitoring and Assessment Program's benthic index for the Virginian Province and the New York Harbor benthic index of biotic integrity using data collected in New York Harbor and evaluated overall agreement of the indices and associations between each index and measures of habitat and sediment contamination. The indices agreed in approximately 78% of the cases. The New York Harbor benthic index of biotic integrity showed stronger associations with sediment metal contamination and grain size.

Beyond data management: how ecoinformatics can benefit environmental monitoring programs.

We review ways in which the new discipline of ecoinformatics is changing how environmental monitoring data are managed, synthesized, and analyzed. Rapid improvements in information technology and strong interest in biodiversity and sustainable ecosystems are driving a vigorous phase of development in ecological databases. Emerging data standards and protocols enable these data to be shared in ways that have previously been difficult. We use the U.S. Environmental Protection Agency's National Coastal Assessment (NCA) as an example. The NCA has collected biological, chemical, and physical data from thousands of stations around the U.S. coasts since 1990. NCA data that were collected primarily to assess the ecological condition of the U.S. coasts can be used in innovative ways, such as biogeographical studies to analyze species invasions. NCA application of ecoinformatics tools leads to new possibilities for integrating the hundreds of thousands of NCA species records with other databases to address broad-scale and long-term questions such as environmental impacts, global climate change, and species invasions.

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.

CProb: a computational tool for conducting conditional probability analysis.

Conditional probability is the probability of observing one event given that another event has occurred. In an environmental context, conditional probability helps to assess the association between an environmental contaminant (i.e., the stressor) and the ecological condition of a resource (i.e., the response). These analyses, when combined with controlled experiments and other methodologies, show great promise in evaluating ecological conditions from observational data and in defining water quality and other environmental criteria. Current applications of conditional probability analysis (CPA) are largely done via scripts or cumbersome spreadsheet routines, which may prove daunting to end-users and do not provide access to the underlying scripts. Combining spreadsheets with scripts eases computation through a familiar interface (i.e., Microsoft Excel) and creates a transparent process through full accessibility to the scripts. With this in mind, we developed a software application, CProb, as an Add-in for Microsoft Excel with R, R(D)com Server, and Visual Basic for Applications. CProb calculates and plots scatterplots, empirical cumulative distribution functions, and conditional probability. In this short communication, we describe CPA, our motivation for developing a CPA tool, and our implementation of CPA as a Microsoft Excel Add-in. Further, we illustrate the use of our software with two examples: a water quality example and a landscape example. CProb is freely available for download at http://www.epa.gov/emap/nca/html/regions/cprob.

Predicting estuarine sediment metal concentrations and inferred ecological conditions: an information theoretic approach.

Empirically derived relationships associating sediment metal concentrations with degraded ecological conditions provide important information to assess estuarine condition. Resources limit the number, magnitude, and frequency of monitoring activities to acquire these data. Models that use available information and simple statistical relationships to predict sediment metal concentrations could provide an important tool for environmental assessment. We developed 45 predictive models for the total concentrations of copper, lead, mercury, and cadmium in estuarine sediments along the Southern New England and Mid-Atlantic regions of the United States. Using information theoretic model-averaging approaches, we found total developed land and percent silt/clay of estuarine sediment were the most important variables for predicting the presence of all four metals. Estuary area, river flow, tidal range, and total agricultural land varied in their importance. The model-averaged predictions explained 78.4, 70.5, 56.4, and 50.3% of the variation for copper, lead, mercury, and cadmium, respectively. Overall prediction accuracies of selected sediment benchmark values (i.e., effects ranges) were 83.9, 84.8, 78.6, and 92.0% for copper, lead, mercury, and cadmium, respectively. Our results further support the generally accepted conclusion that sediment metal concentrations are best described by the physical characteristics of the estuarine sediment and the total amount of urban land in the contributing watershed. We demonstrated that broad-scale predictive models built from existing monitoring data with information theoretic model-averaging approaches provide valuable predictions of estuarine sediment metal concentrations and show promise for future environmental modeling efforts in other regions.

Top-down and bottom-up controls on southern New England salt marsh crab populations.

Southern New England salt marsh vegetation and habitats are changing rapidly in response to sea-level rise. At the same time, fiddler crab (Uca spp.) distributions have expanded and purple marsh crab (Sesarma reticulatum) grazing on creekbank vegetation has increased. Sea-level rise and reduced predation pressure drive these changing crab populations but most studies focus on one species; there is a need for community-level assessments of impacts from multiple crab species. There is also a need to identify additional factors that can affect crab populations. We sampled crabs and environmental parameters in four Rhode Island salt marshes in 2014 and compiled existing data to quantify trends in crab abundance and multiple factors that potentially affect crabs. Crab communities were dominated by fiddler and green crabs (Carcinus maenas); S. reticulatum was much less abundant. Burrow sizes suggest that Uca is responsible for most burrows. On the marsh platform, burrows and Carcinus abundance were negatively correlated with elevation, soil moisture, and soil percent organic matter and positively correlated with soil bulk density. Uca abundance was negatively correlated with Spartina patens cover and height and positively correlated with Spartina alterniflora cover and soil shear strength. Creekbank burrow density increased dramatically between 1998 and 2016. During the same time, fishing effort and the abundance of birds that prey on crabs decreased, and water levels increased. Unlike in other southern New England marshes where recreational overfishing is hypothesized to drive increasing marsh crab abundance, we propose that changes in crab abundance were likely unrelated to recreational finfish over-harvest; instead, they better track sea-level rise and changing abundances of alternate predators, such as birds. We predict that marsh crab abundance will continue to expand with ongoing sea-level rise, at least until inundation thresholds for crab survival are exceeded.

Performance of National Maps of Watershed Integrity at Watershed Scales.

Watershed integrity, the capacity of a watershed to support and maintain ecological processes essential to the sustainability of services provided to society, can be influenced by a range of landscape and in-stream factors. Ecological response data from four intensively monitored case study watersheds exhibiting a range of environmental conditions and landscape characteristics across the United States were used to evaluate the performance of a national level Index of Watershed Integrity (IWI) at regional and local watershed scales. Using Pearson's correlation coefficient (r), and Spearman's rank correlation coefficient (rs ), response variables displayed highly significant relationships and were significantly correlated with IWI and ICI (Index of Catchment Integrity) values at all watersheds. Nitrogen concentration and flux-related watershed response metrics exhibited significantly strong negative correlations across case study watersheds, with absolute correlations (|r|) ranging from 0.48 to 0.97 for IWI values, and 0.31 to 0.96 for ICI values. Nitrogen-stable isotope ratios measured in chironomids and periphyton from streams and benthic organic matter from lake sediments also demonstrated strong negative correlations with IWI values, with |r| ranging from 0.47 to 0.92, and 0.35 to 0.89 for correlations with ICI values. This evaluation of the performance of national watershed and catchment integrity metrics and their strong relationship with site level responses provides weight-of-evidence support for their use in state, local and regionally focused applications.

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.

Estimating summer nutrient concentrations in Northeastern lakes from SPARROW load predictions and modeled lake depth and volume.

Global nutrient cycles have been altered by the use of fossil fuels and fertilizers resulting in increases in nutrient loads to aquatic systems. In the United States, excess nutrients have been repeatedly reported as the primary cause of lake water quality impairments. Setting nutrient criteria that are protective of a lakes ecological condition is one common solution; however, the data required to do this are not always easily available. A useful solution for this is to combine available field data (i.e., The United States Environmental Protection Agency (USEPA) National Lake Assessment (NLA)) with average annual nutrient load models (i.e., USGS SPARROW model) to estimate summer concentrations across a large number of lakes. In this paper we use this combined approach and compare the observed total nitrogen (TN) and total phosphorus (TN) concentrations in Northeastern lakes from the 2007 National Lake Assessment to those predicted by the Northeast SPARROW model. We successfully adjusted the SPARROW predictions to the NLA observations with the use of Vollenweider equations, simple input-output models that predict nutrient concentrations in lakes based on nutrient loads and hydraulic residence time. This allows us to better predict summer concentrations of TN and TP in Northeastern lakes and ponds. On average we improved our predicted concentrations of TN and TP with Vollenweider models by 18.7% for nitrogen and 19.0% for phosphorus. These improved predictions are being used in other studies to model ecosystem services (e.g., aesthetics) and dis-services (e.g. cyanobacterial blooms) for ~18,000 lakes in the Northeastern United States.

Predicting maximum lake depth from surrounding topography.

Information about lake morphometry (e.g., depth, volume, size, etc.) aids understanding of the physical and ecological dynamics of lakes, yet is often not readily available. The data needed to calculate measures of lake morphometry, particularly lake depth, are usually collected on a lake-by-lake basis and are difficult to obtain across broad regions. To span the gap between studies of individual lakes where detailed data exist and regional studies where access to useful data on lake depth is unavailable, we developed a method to predict maximum lake depth from the slope of the topography surrounding a lake. We use the National Elevation Dataset and the National Hydrography Dataset - Plus to estimate the percent slope of surrounding lakes and use this information to predict maximum lake depth. We also use field measured maximum lake depths from the US EPA's National Lakes Assessment to empirically adjust and cross-validate our predictions. We were able to predict maximum depth for ∼28,000 lakes in the Northeastern United States with an average cross-validated RMSE of 5.95 m and 5.09 m and average correlation of 0.82 and 0.69 for Hydrological Unit Code Regions 01 and 02, respectively. The depth predictions and the scripts are openly available as supplements to this manuscript.