Spatiotemporal mapping of Karenia brevis blooms

Measurement of cell concentrations

Measurements of Karenia brevis cell abundance have been compiled by the Fish and Wildlife Research Institute, Florida Fish and Wildlife Conservation Commission, for more than seventy years. Blooms are most frequent along the southwest Florida coast, between Anna Maria Island (N lat 27.54) and Marco Island (N lat 25.80), where newly 120,000 samples have been collected. The sampling rate has varied greatly over this period:

The 1954-57 period of intense sampling, funded by the federal government, followed very severe blooms along the Southwest Florida coast in the late 1940s and early 1950s. Sampling rates were greatly decreased in following years due to lack of funds, except for the uptick in the late-1970s. This was the result of studies proving that fish kills accompanying blooms in 1976-78 were caused by brevetoxins produced by K. brevis (then known as Gymnodinium breve), not oxygen deprivation due decaying organic matter. The very severe 1994-96 bloom resulted in action by the Florida Legislature establishing a routine, proactive sampling protocol that continues today.

Data processing: separating signal from noise

Blooms of K. brevis are inherently heterogeneous; even in the midst of a severe bloom, samples taken from the same general area can have cell concentrations ranging from zero to millions of cells per liter. This is caused by current and wind forcing, vertical migration, and localized nutrient sources. As a result, time series of sample cell concentrations form a cloud of points, not a line. The first step to understanding the temporal structure of a bloom is the removal of the noise hiding that structure. This was done with a three-stage cascaded filter:

  1. simple average of the daily measurements to remove sub-daily spatial patchiness;
  2. seven-day, centered moving average acting as a notch filter to remove the weekly sampling protocol artifact; and
  3. 21-day, centered moving average filter to remove transient, short-term events while preserving the lower-frequency oscillations that represent the true temporal structure of the bloom.

Here is the result of this process:

ecause blooms are hetroA prominent feature observed across these smoothed timelines is that severe red tides are rarely singular, linear events. Instead, sustained disruptions—such as the extended 2017–2019 interval—frequently manifest as a sequence of discrete, blended population peaks occurring roughly 100 days apart.

While the ultimate magnitude and footprint of a bloom are heavily governed by extrinsic environmental factors—including riverine runoff, upwelling events, and nearshore nutrient loading—the persistent ~100-day cadence suggests an important interaction with the organism’s internal growth dynamics.

Integrating Population Kinetics with Environmental Forcing

To explore the mechanism behind this multi-wave persistence, these historical trends can be modeled using a time-delay logistic framework:

In this structural model, environmental resources and nutrient availability dictate the system’s carrying capacity (), while an intrinsic biological latency parameter ()—rooted in the K. brevis cell division cycle and reproductive constraints—introduces a delayed density-dependent feedback.

When regional nutrient conditions elevate the carrying capacity, this biological lag allows the population to temporarily overshoot immediate baseline limits before experiencing a feedback response. Mathematically, a 24-day internal delay operates near a Hopf bifurcation, naturally generating sustained oscillation periods () that closely mirror the empirical peak separations observed in the field.

Implications for Resource Management

Reframing bloom tracking through this combined lens highlights a vital metric for ecological and economic impact assessment: total bloom severity is a function of both average concentration and cumulative duration (). Recognizing the intrinsic temporal rhythms that sustain these multi-wave events can help resource managers better anticipate bloom longevity and optimize regional mitigation strategies.

 

 

 


Sep 1994 – Apr 1997

Jun 1997 – Feb 2001

May 2001 – Jun 2004

Aug 2004 – Jun 2007

Jun 2011 – Feb 2015

Jun 2015 – Jul 2017

Aug 2017 – Apr 2019

Aug 2019 – Feb 2020

Sep 2020 – Jan 2022

Aug 2022 – Jul 2023

Aug 2024 – Mar 2025