Satellites quantify the spatial extent of cyanobacterial blooms across the United States at multiple scales.

Previous studies indicate that cyanobacterial harmful algal bloom (cyanoHAB) frequency, extent, and magnitude have increased globally over the past few decades. However, little quantitative capability is available to assess these metrics of cyanoHABs across broad geographic scales and at regular intervals. Here, the spatial extent was quantified from a cyanobacteria algorithm applied to two European Space Agency satellite platforms-the MEdium Resolution Imaging Spectrometer (MERIS) onboard Envisat and the Ocean and Land Colour Instrument (OLCI) onboard Sentinel-3. CyanoHAB spatial extent was defined for each geographic area as the percentage of valid satellite pixels that exhibited cyanobacteria above the detection limit of the satellite sensor. This study quantified cyanoHAB spatial extent for over 2,000 large lakes and reservoirs across the contiguous United States (CONUS) during two time periods: 2008-2011 via MERIS and 2017-2020 via OLCI when cloud-, ice-, and snow-free imagery was available. Approximately 56% of resolvable lakes were glaciated, 13% were headwater, isolated, or terminal lakes, and the rest were primarily drainage lakes. Results were summarized at national-, regional-, state-, and lake-scales, where regions were defined as nine climate regions which represent climatically consistent states. As measured by satellite, changes in national cyanoHAB extent did have a strong increase of 6.9% from 2017 to 2020 (|Kendall's tau (τ)| = 0.56; gamma (γ) = 2.87 years), but had negligible change (|τ| = 0.03) from 2008 to 2011. Two of the nine regions had moderate (0.3 ≤ |τ| < 0.5) increases in spatial extent from 2017 to 2020, and eight of nine regions had negligible (|τ| < 0.2) change from 2008 to 2011. Twelve states had a strong or moderate increase from 2017 to 2020 (|τ| ≥ 0.3), while only one state had a moderate increase and two states had a moderate decrease from 2008 to 2011. A decrease, or no change, in cyanoHAB spatial extent did not indicate a lack of issues related to cyanoHABs. Sensitivity results of randomly omitted daily CONUS scenes confirm that even with reduced data availability during a short four-year temporal assessment, the direction and strength of the changes in spatial extent remained consistent. We present the first set of national maps of lake cyanoHAB spatial extent across CONUS and demonstrate an approach for quantifying past and future changes at multiple spatial scales. Results presented here provide water quality managers information regarding current cyanoHAB spatial extent and quantify rates of change.

Merging of the Case 2 Regional Coast Colour and Maximum-Peak Height chlorophyll-a algorithms: validation and demonstration of satellite-derived retrievals across US lakes.

Water quality monitoring is relevant for protecting the designated, or beneficial uses, of water such as drinking, aquatic life, recreation, irrigation, and food supply that support the economy, human well-being, and aquatic ecosystem health. Managing finite water resources to support these designated uses requires information on water quality so that managers can make sustainable decisions. Chlorophyll-a (chl-a, µg L-1) concentration can serve as a proxy for phytoplankton biomass and may be used as an indicator of increased anthropogenic nutrient stress. Satellite remote sensing may present a complement to in situ measures for assessments of water quality through the retrieval of chl-a with in-water algorithms. Validation of chl-a algorithms across US lakes improves algorithm maturity relevant for monitoring applications. This study compares performance of the Case 2 Regional Coast Colour (C2RCC) chl-a retrieval algorithm, a revised version of the Maximum-Peak Height (MPH(P)) algorithm, and three scenarios merging these two approaches. Satellite data were retrieved from the MEdium Resolution Imaging Spectrometer (MERIS) and the Ocean and Land Colour Instrument (OLCI), while field observations were obtained from 181 lakes matched with U.S. Water Quality Portal chl-a data. The best performance based on mean absolute multiplicative error (MAEmult) was demonstrated by the merged algorithm referred to as C15-M10 (MAEmult = 1.8, biasmult = 0.97, n = 836). In the C15-M10 algorithm, the MPH(P) chl-a value was retained if it was > 10 µg L-1; if the MPH(P) value was ≤ 10 µg L-1, the C2RCC value was selected, as long as that value was < 15 µg L-1. Time-series and lake-wide gradients compared against independent assessments from Lake Champlain and long-term ecological research stations in Wisconsin were used as complementary examples supporting water quality reporting requirements. Trophic state assessments for Wisconsin lakes provided examples in support of inland water quality monitoring applications. This study presents and assesses merged adaptations of chl-a algorithms previously reported independently. Additionally, it contributes to the transition of chl-a algorithm maturity by quantifying error statistics for a number of locations and times.

Satellite-derived cyanobacteria frequency and magnitude in headwaters & near-dam reservoir surface waters of the Southern U.S.

Reservoirs are dominant features of the modern hydrologic landscape and provide vital services. However, the unique morphology of reservoirs can create suitable conditions for excessive algae growth and associated cyanobacteria blooms in shallow in-flow reservoir locations by providing warm water environments with relatively high nutrient inputs, deposition, and nutrient storage. Cyanobacteria harmful algal blooms (cyanoHAB) are costly water management issues and bloom recurrence is associated with economic costs and negative impacts to human, animal, and environmental health. As cyanoHAB occurrence varies substantially within different regions of a water body, understanding in-lake cyanoHAB spatial dynamics is essential to guide reservoir monitoring and mitigate potential public exposure to cyanotoxins. Cloud-based computational processing power and high temporal frequency of satellites enables advanced pixel-based spatial analysis of cyanoHAB frequency and quantitative assessment of reservoir headwater in-flows compared to near-dam surface waters of individual reservoirs. Additionally, extensive spatial coverage of satellite imagery allows for evaluation of spatial trends across many dozens of reservoir sites. Surface water cyanobacteria concentrations for sixty reservoirs in the southern U.S. were estimated using 300 m resolution European Space Agency (ESA) Ocean and Land Colour Instrument (OLCI) satellite sensor for a five year period (May 2016-April 2021). Of the reservoirs studied, spatial analysis of OLCI data revealed 98% had more frequent cyanoHAB occurrence above the concentration of >100,000 cells/mL in headwaters compared to near-dam surface waters (P < 0.001). Headwaters exhibited greater seasonal variability with more frequent and higher magnitude cyanoHABs occurring mid-summer to fall. Examination of reservoirs identified extremely high concentration cyanobacteria events (>1,000,000 cells/mL) occurring in 70% of headwater locations while only 30% of near-dam locations exceeded this threshold. Wilcoxon signed-rank tests of cyanoHAB magnitudes using paired-observations (dates with observations in both a reservoir's headwater and near-dam locations) confirmed significantly higher concentrations in headwater versus near-dam locations (p < 0.001).

Acute health effects associated with satellite-determined cyanobacterial blooms in a drinking water source in Massachusetts.

BACKGROUND: The occurrence of cyanobacterial blooms in freshwater presents a threat to human health. However, epidemiological studies on the association between cyanobacterial blooms in drinking water sources and human health outcomes are scarce. The objective of this study was to evaluate if cyanobacterial blooms were associated with increased emergency room visits for gastrointestinal (GI), respiratory and dermal illnesses. METHODS: Satellite-derived cyanobacteria cell concentrations were estimated in the source of drinking water for the Greater Boston area, during 2008-2011. Daily counts of hospital emergency room visits for GI, respiratory and dermal illnesses among drinking water recipients were obtained from an administrative record database. A two-stage model was used to analyze time-series data for an association between cyanobacterial blooms and the occurrence of illnesses. At the first stage, predictive autoregressive generalized additive models for Poisson-distributed outcomes were fitted to daily illness count data and daily predictive variables. At the second stage, residuals from the first stage models were regressed against lagged categorized cyanobacteria concentration estimates. RESULTS: The highest cyanobacteria concentration (above the 75th percentile) was associated with an additional 4.3 cases of respiratory illness (95% confidence interval: 0.7, 8.0, p = 0.02, n = 268) compared to cyanobacteria concentrations below the 50th percentile in a two-day lag. There were no significant associations between satellite derived cyanobacterial concentrations and lagged data on GI or dermal illnesses. CONCLUSION: The study demonstrated a significant positive association between satellite-derived cyanobacteria concentrations in source water and respiratory illness occurring 2 days later. Future studies will require direct measures of cyanotoxins and health effects associated with exposure to cyanobacteria-impacted drinking water sources.

Assessing cyanobacterial frequency and abundance at surface waters near drinking water intakes across the United States.

This study presents the first large-scale assessment of cyanobacterial frequency and abundance of surface water near drinking water intakes across the United States. Public water systems serve drinking water to nearly 90% of the United States population. Cyanobacteria and their toxins may degrade the quality of finished drinking water and can lead to negative health consequences. Satellite imagery can serve as a cost-effective and consistent monitoring technique for surface cyanobacterial blooms in source waters and can provide drinking water treatment operators information for managing their systems. This study uses satellite imagery from the European Space Agency's Ocean and Land Colour Instrument (OLCI) spanning June 2016 through April 2020. At 300-m spatial resolution, OLCI imagery can be used to monitor cyanobacteria in 685 drinking water sources across 285 lakes in 44 states, referred to here as resolvable drinking water sources. First, a subset of satellite data was compared to a subset of responses (n = 84) submitted as part of the U.S. Environmental Protection Agency's fourth Unregulated Contaminant Monitoring Rule (UCMR 4). These UCMR 4 qualitative responses included visual observations of algal bloom presence and absence near drinking water intakes from March 2018 through November 2019. Overall agreement between satellite imagery and UCMR 4 qualitative responses was 94% with a Kappa coefficient of 0.70. Next, temporal frequency of cyanobacterial blooms at all resolvable drinking water sources was assessed. In 2019, bloom frequency averaged 2% and peaked at 100%, where 100% indicated a bloom was always present at the source waters when satellite imagery was available. Monthly cyanobacterial abundances were used to assess short-term trends across all resolvable drinking water sources and effect size was computed to provide insight on the number of years of data that must be obtained to increase confidence in an observed change. Generally, 2016 through 2020 was an insufficient time period for confidently observing changes at these source waters; on average, a decade of satellite imagery would be required for observed environmental trends to outweigh variability in the data. However, five source waters did demonstrate a sustained short-term trend, with one increasing in cyanobacterial abundance from June 2016 to April 2020 and four decreasing.

Satellite remote sensing to assess cyanobacterial bloom frequency across the United States at multiple spatial scales.

Cyanobacterial blooms can have negative effects on human health and local ecosystems. Field monitoring of cyanobacterial blooms can be costly, but satellite remote sensing has shown utility for more efficient spatial and temporal monitoring across the United States. Here, satellite imagery was used to assess the annual frequency of surface cyanobacterial blooms, defined for each satellite pixel as the percentage of images for that pixel throughout the year exhibiting detectable cyanobacteria. Cyanobacterial frequency was assessed across 2,196 large lakes in 46 states across the continental United States (CONUS) using imagery from the European Space Agency's Ocean and Land Colour Instrument for the years 2017 through 2019. In 2019, across all satellite pixels considered, annual bloom frequency had a median value of 4% and a maximum value of 100%, the latter indicating that for those satellite pixels, a cyanobacterial bloom was detected by the satellite sensor for every satellite image considered. In addition to annual pixel-scale cyanobacterial frequency, results were summarized at the lake- and state-scales by averaging annual pixel-scale results across each lake and state. For 2019, average annual lake-scale frequencies also had a maximum value of 100%, and Oregon and Ohio had the highest average annual state-scale frequencies at 65% and 52%. Pixel-scale frequency results can assist in identifying portions of a lake that are more prone to cyanobacterial blooms, while lake- and state-scale frequency results can assist in the prioritization of sampling resources and mitigation efforts. Satellite imagery is limited by the presence of snow and ice, as imagery collected in these conditions are quality flagged and discarded. Thus, annual bloom frequencies within nine climate regions were investigated to determine whether missing data biased results in climate regions more prone to snow and ice, given that their annual summaries would be weighted toward the summer months when cyanobacterial blooms tend to occur. Results were unbiased by the time period selected in most climate regions, but a large bias was observed for the Northwest Rockies and Plains climate region. Moderate biases were observed for the Ohio Valley and the Southeast climate regions. Finally, a clustering analysis was used to identify areas of high and low cyanobacterial frequency across CONUS based on average annual lake-scale cyanobacterial frequencies for 2019. Several clusters were identified that transcended state, watershed, and eco-regional boundaries. Combined with additional data, results from the clustering analysis may offer insight regarding large-scale drivers of cyanobacterial blooms.

Spatio-Temporal Modeling for Forecasting High-Risk Freshwater Cyanobacterial Harmful Algal Blooms in Florida.

Due to the occurrence of more frequent and widespread toxic cyanobacteria events, the ability to predict freshwater cyanobacteria harmful algal blooms (cyanoHAB) is of critical importance for the management of drinking and recreational waters. Lake system specific geographic variation of cyanoHABs has been reported, but regional and state level variation is infrequently examined. A spatio-temporal modeling approach can be applied, via the computationally efficient Integrated Nested Laplace Approximation (INLA), to high-risk cyanoHAB exceedance rates to explore spatio-temporal variations across statewide geographic scales. We explore the potential for using satellite-derived data and environmental determinants to develop a short-term forecasting tool for cyanobacteria presence at varying space-time domains for the state of Florida. Weekly cyanobacteria abundance data were obtained using Sentinel-3 Ocean Land Color Imagery (OLCI), for a period of May 2016-June 2019. Time and space varying covariates include surface water temperature, ambient temperature, precipitation, and lake geomorphology. The hierarchical Bayesian spatio-temporal modeling approach in R-INLA represents a potential forecasting tool useful for water managers and associated public health applications for predicting near future high-risk cyanoHAB occurrence given the spatio-temporal characteristics of these events in the recent past. This method is robust to missing data and unbalanced sampling between waterbodies, both common issues in water quality datasets.

Quantifying national and regional cyanobacterial occurrence in US lakes using satellite remote sensing.

Cyanobacterial harmful algal blooms are the most common form of harmful algal blooms in freshwater systems throughout the world. However, in situ sampling of cyanobacteria in inland lakes is limited both spatially and temporally. Satellite data has proven to be an effective tool to monitor cyanobacteria in freshwater lakes across the United States. This study uses data from the European Space Agency Envisat MEdium Resolution Imaging Spectrometer and the Sentinel-3 Ocean and Land Color Instrument to provide a national overview of the percentage of lakes experiencing a cyanobacterial bloom on a weekly basis for 2008-2011, 2017, and 2018. A total of 2321 lakes across the contiguous United States were included in the analysis. We examined four different thresholds to define when a waterbody is classified as experiencing a bloom. Across these four thresholds, we explored variability in bloom percentage with changes in seasonality and lake size. As a validation of algorithm performance, we analyzed the agreement between satellite observations and previously established ecological patterns, although data availability in the wintertime limited these comparisons on a year-round basis. Changes in cyanobacterial bloom percentage at the national scale followed the well-known temporal pattern of freshwater blooms. The percentage of lakes experiencing a bloom increased throughout the year, reached a maximum in fall, and decreased through the winter. Wintertime data, particularly in northern regions, were consistently limited due to snow and ice cover. With the exception of the Southeast and South, regional patterns mimicked patterns found at the national scale. The Southeast and South exhibited an unexpected pattern as cyanobacterial bloom percentage reached a maximum in the winter rather than the summer. Lake Jesup in Florida was used as a case study to validate this observed pattern against field observations of chlorophyll a. Results from this research establish a baseline of annual occurrence of cyanobacterial blooms in inland lakes across the United States. In addition, methods presented in this study can be tailored to fit the specific requirements of an individual system or region.

Envisat MERIS and Sentinel-3 OLCI satellite lake biophysical water quality flag dataset for the contiguous United States.

Monitoring lake biophysical water quality is a global challenge. Satellite remote sensing offers a technology for continuous water quality information in data poor regions throughout the United States. Quality assurance flag data are provided for the presence of snow/ice, land-adjacency, and unresolvable waterbodies supporting water quality derived measures from Envisat MEdium Resolution Imaging Spectrometer and Sentinel-3 Ocean and Land Colour Instrument for the continental United States. In addition, an updated Waterbody Data mask that contains valid waterbody and coastal ocean delineation is provided. The quality assurance flag datasets can benefit the scientific community in processing lake water quality throughout the contiguous United States by addressing errors from snow/ice, land adjacency, and land masking. The dataset presented here will be used in the development of national scale metrics for derived biophysical water quality in the US.

Forecasting Seasonal Vibrio parahaemolyticus Concentrations in New England Shellfish.

Seafood-borne Vibrio parahaemolyticus illness is a global public health issue facing resource managers and the seafood industry. The recent increase in shellfish-borne illnesses in the Northeast United States has resulted in the application of intensive management practices based on a limited understanding of when and where risks are present. We aim to determine the contribution of factors that affect V. parahaemolyticus concentrations in oysters (Crassostrea virginica) using ten years of surveillance data for environmental and climate conditions in the Great Bay Estuary of New Hampshire from 2007 to 2016. A time series analysis was applied to analyze V. parahaemolyticus concentrations and local environmental predictors and develop predictive models. Whereas many environmental variables correlated with V. parahaemolyticus concentrations, only a few retained significance in capturing trends, seasonality and data variability. The optimal predictive model contained water temperature and pH, photoperiod, and the calendar day of study. The model enabled relatively accurate seasonality-based prediction of V. parahaemolyticus concentrations for 2014-2016 based on the 2007-2013 dataset and captured the increasing trend in extreme values of V. parahaemolyticus concentrations. The developed method enables the informative tracking of V. parahaemolyticus concentrations in coastal ecosystems and presents a useful platform for developing area-specific risk forecasting models.

Mobile device application for monitoring cyanobacteria harmful algal blooms using Sentinel-3 satellite Ocean and Land Colour Instruments.

Cyanobacterial harmful algal blooms (cyanoHAB) cause human and ecological health problems in lakes worldwide. The timely distribution of satellite-derived cyanoHAB data is necessary for adaptive water quality management and for targeted deployment of water quality monitoring resources. Software platforms that permit timely, useful, and cost-effective delivery of information from satellites are required to help managers respond to cyanoHABs. The Cyanobacteria Assessment Network (CyAN) mobile device application (app) uses data from the European Space Agency Copernicus Sentinel-3 satellite Ocean and Land Colour Instrument (OLCI) in near realtime to make initial water quality assessments and quickly alert managers to potential problems and emerging threats related to cyanobacteria. App functionality and satellite data were validated with 25 state health advisories issued in 2017. The CyAN app provides water quality managers with a user-friendly platform that reduces the complexities associated with accessing satellite data to allow fast, efficient, initial assessments across lakes.

An initial validation of Landsat 5 and 7 derived surface water temperature for U.S. lakes, reservoirs, and estuaries.

The United States Harmful Algal Bloom and Hypoxia Research Control Act of 2014 identified the need for forecasting and monitoring harmful algal blooms (HAB) in lakes, reservoirs, and estuaries across the nation. Temperature is a driver in HAB forecasting models that affects both HAB growth rates and toxin production. Therefore, temperature data derived from the U.S. Geological Survey Landsat 5 Thematic Mapper and Landsat 7 Enhanced Thematic Mapper Plus thermal band products were validated across 35 lakes and reservoirs, and 24 estuaries. In situ data from the Water Quality Portal (WQP) were used for validation. The WQP serves data collected by state, federal, and tribal groups. Discrete in situ temperature data included measurements at 11,910 U.S. lakes and reservoirs from 1980 through 2015. Landsat temperature measurements could include 170,240 lakes and reservoirs once an operational product is achieved. The Landsat-derived temperature mean absolute error was 1.34°C in lake pixels >180 m from land, 4.89°C at the land-water boundary, and 1.11°C in estuaries based on comparison against discrete surface in situ measurements. This is the first study to quantify Landsat resolvable U.S. lakes and reservoirs, and large-scale validation of an operational satellite provisional temperature climate data record algorithm. Due to the high performance of open water pixels, Landsat satellite data may supplement traditional in situ sampling by providing data for most U.S. lakes, reservoirs, and estuaries over consistent seasonal intervals (even with cloud cover) for an extended period of record of more than 35 years.

A method for examining temporal changes in cyanobacterial harmful algal bloom spatial extent using satellite remote sensing.

Cyanobacterial harmful algal blooms (CyanoHAB) are thought to be increasing globally over the past few decades, but relatively little quantitative information is available about the spatial extent of blooms. Satellite remote sensing provides a potential technology for identifying cyanoHABs in multiple water bodies and across geo-political boundaries. An assessment method was developed using MEdium Resolution Imaging Spectrometer (MERIS) imagery to quantify cyanoHAB surface area extent, transferable to different spatial areas, in Florida, Ohio, and California for the test period of 2008 to 2012. Temporal assessment was used to evaluate changes in satellite resolvable inland waterbodies for each state of interest. To further assess cyanoHAB risk within the states, the World Health Organization's (WHO) recreational guidance level thresholds were used to categorize surface area of cyanoHABs into three risk categories: low, moderate, and high-risk bloom area. Results showed that in Florida, the area of cyanoHABs increased largely due to observed increases in high-risk bloom area. California exhibited a slight decrease in cyanoHAB extent, primarily attributed to decreases in Northern California. In Ohio (excluding Lake Erie), little change in cyanoHAB surface area was observed. This study uses satellite remote sensing to quantify changes in inland cyanoHAB surface area across numerous water bodies within an entire state. The temporal assessment method developed here will be relevant into the future as it is transferable to the Ocean Land Colour Instrument (OLCI) on Sentinel-3A/3B missions.

Satellite monitoring of cyanobacterial harmful algal bloom frequency in recreational waters and drinking source waters.

Cyanobacterial harmful algal blooms (cyanoHAB) cause extensive problems in lakes worldwide, including human and ecological health risks, anoxia and fish kills, and taste and odor problems. CyanoHABs are a particular concern in both recreational waters and drinking source waters because of their dense biomass and the risk of exposure to toxins. Successful cyanoHAB assessment using satellites may provide an indicator for human and ecological health protection, In this study, methods were developed to assess the utility of satellite technology for detecting cyanoHAB frequency of occurrence at locations of potential management interest. The European Space Agency's MEdium Resolution Imaging Spectrometer (MERIS) was evaluated to prepare for the equivalent series of Sentine1-3 Ocean and Land Colour Imagers (OLCI) launched in 2016 as part of the Copernicus program. Based on the 2012 National Lakes Assessment site evaluation guidelines and National Hydrography Dataset, the continental United States contains 275,897 lakes and reservoirs >1 hectare in area. Results from this study show that 5.6 % of waterbodies were resolvable by satellites with 300 m single-pixel resolution and 0.7 % of waterbodies were resolvable when a three by three pixel (3×3-pixel) array was applied based on minimum Euclidian distance from shore. Satellite data were spatially joined to U.S. public water surface intake (PWSI) locations, where single-pixel resolution resolved 57% of the PWSI locations and a 3×3-pixel array resolved 33% of the PWSI locations. Recreational and drinking water sources in Florida and Ohio were ranked from 2008 through 2011 by cyanoHAB frequency above the World Health Organization's (WHO) high threshold for risk of 100,000 cells mL-1. The ranking identified waterbodies with values above the WHO high threshold, where Lake Apopka, FL (99.1 %) and Grand Lake St. Marys, OH (83 %) had the highest observed bloom frequencies per region. The method presented here may indicate locations with high exposure to cyanoHABs and therefore can be used to assist in prioritizing management resources and actions for recreational and drinking water sources.

Environmental Conditions Associated with Elevated Vibrio parahaemolyticus Concentrations in Great Bay Estuary, New Hampshire.

Reports from state health departments and the Centers for Disease Control and Prevention indicate that the annual number of reported human vibriosis cases in New England has increased in the past decade. Concurrently, there has been a shift in both the spatial distribution and seasonal detection of Vibrio spp. throughout the region based on limited monitoring data. To determine environmental factors that may underlie these emerging conditions, this study focuses on a long-term database of Vibrio parahaemolyticus concentrations in oyster samples generated from data collected from the Great Bay Estuary, New Hampshire over a period of seven consecutive years. Oyster samples from two distinct sites were analyzed for V. parahaemolyticus abundance, noting significant relationships with various biotic and abiotic factors measured during the same period of study. We developed a predictive modeling tool capable of estimating the likelihood of V. parahaemolyticus presence in coastal New Hampshire oysters. Results show that the inclusion of chlorophyll a concentration to an empirical model otherwise employing only temperature and salinity variables, offers improved predictive capability for modeling the likelihood of V. parahaemolyticus in the Great Bay Estuary.

Uncertainty in model predictions of Vibrio vulnificus response to climate variability and change: a Chesapeake Bay case study.

The effect that climate change and variability will have on waterborne bacteria is a topic of increasing concern for coastal ecosystems, including the Chesapeake Bay. Surface water temperature trends in the Bay indicate a warming pattern of roughly 0.3-0.4°C per decade over the past 30 years. It is unclear what impact future warming will have on pathogens currently found in the Bay, including Vibrio spp. Using historical environmental data, combined with three different statistical models of Vibrio vulnificus probability, we explore the relationship between environmental change and predicted Vibrio vulnificus presence in the upper Chesapeake Bay. We find that the predicted response of V. vulnificus probability to high temperatures in the Bay differs systematically between models of differing structure. As existing publicly available datasets are inadequate to determine which model structure is most appropriate, the impact of climatic change on the probability of V. vulnificus presence in the Chesapeake Bay remains uncertain. This result points to the challenge of characterizing climate sensitivity of ecological systems in which data are sparse and only statistical models of ecological sensitivity exist.