Modelling climate change scenarios and their impact on microalgae communities in the Coorong

2026

Daniel Chilton, Justin Brookes, Matthew Hipsey, Sophie C. Leterme
The Coorong, a Ramsar‑listed coastal lagoon of ecological, cultural and economic importance, is undergoing change as a result of climate‑driven shifts in freshwater inflows, temperature, sea level and extreme weather events. Phytoplankton are important bioindicators of these changes due to their rapid turnover rates and are foundational to the functioning of the Coorong’s food web. This report synthesises 16 years of phytoplankton and environmental data (2007–2023) to examine how climate‑mediated flow variation shapes the Coorong’s microalgal communities, supporting a predictive framework to understand future environmental trajectories and phytoplankton responses under climate change. Through an integrated analysis combining long‑term monitoring, statistical modelling, machine learning and habitat‑based ecological characterisation, this study delivers a comprehensive assessment of how climate variability affects phytoplankton community dynamics and harmful algal bloom (HAB) risk in the Coorong.
A full range of hydrological conditions were captured in this study, from the Millennium Drought (extreme dry) to a large River Murray flood (extreme wet). These contrasting climate periods produced dramatic shifts in barrage flows, water levels, salinity, temperature and nutrient concentrations, which in turn restructured phytoplankton abundance, diversity and community composition.
The Coorong’s environmental conditions significantly differed across flow periods. During the dry and extreme dry periods, freshwater inputs declined lowering water levels and increasing salinity to levels exceeding 190 PSU in the South Lagoon. Wet and extreme wet periods reversed these conditions, producing rapid freshening, increased hydrological flushing that removed nutrients and salt from the system and increased water levels, that further diluted salt and nutrients. Phytoplankton communities demonstrated strong, statistically significant responses to these changes, with significantly lower abundances under both extreme wet and dry conditions, but higher richness under wetter conditions and higher evenness during drought.
Across an ensemble of analytical approaches, including multiple correlation analysis, ordination, distance‑based linear modelling and generalised additive models, salinity and water level emerged as the dominant drivers of phytoplankton abundance, diversity and community structure. These variables are tightly coupled and both controlled by freshwater inflows. Higher water levels and lower salinities promoted Cyanobacteria and Chlorophyte abundance. Conversely, increasing salinity altered species composition across all taxonomic groups, with increasing dominance of Diatoms and Dinoflagellates. Other important secondary drivers of phytoplankton community dynamics included temperature, distance from the Murray Mouth (representing hydrodynamic and salinity gradients) and nutrient concentrations (particularly NOX and PO₄), although their influence was weaker and more variable across statistical models.
To characterise responses to hydrological change, phytoplankton communities were classified into six salinity‑derived habitats spanning freshwater (H1) to extreme hypersalinity (H6). Community composition, abundance and diversity significantly differed between habitats. Key patterns included:
· Freshwater and low‑brackish habitats (H1–H2) had the highest richness but lower evenness, driven by Cyanobacteria and Chlorophyte dominance, particularly under wetter periods. However, H2 demonstrated the most equal representation of all taxonomic groups,
· High‑brackish to marine conditions (H3) displayed moderate diversity and high evenness and were dominated by Diatoms and Dinoflagellates,
· Hypersaline habitats (H4–H5) supported very high picophytoplankton abundance and Diatom dominance, and
· Extreme hypersaline conditions (H6) had the lowest richness but highest evenness and were dominated by Diatoms and Dinoflagellates with limited representation from other groups.
These habitat‑based community descriptions provide a strong foundation for forecasting future ecological states using hydrological model outputs that predict salinity.
A total of 106 potentially harmful algal bloom forming species (HABs) were identified across the study period, including 82 toxin-producing species. Most HABs occurred in low‑salinity habitats (H1–H2) in the North Lagoon, especially during the post‑drought wet period in 2012. Although HAB events were generally infrequent and short‑lived, two notable bloom events occurred:
1. A Cyanobacteria bloom from January to June 2012 in the North Lagoon consisting of multiple species, though dominated by Aphanocapsa sp. and exceeding medium to high risk thresholds, and
2. A Dinoflagellate bloom during the extreme dry period, where Alexandrium spp. reached action‑level abundances across all sampling sites.
Despite overall low HAB prevalence, the action-level abundances of Alexandrium species capable of producing toxins warrants ongoing vigilance, particularly under future climate scenarios that may increase the frequency and intensity of low‑flow, high‑salinity conditions or sharp nutrient pulses.
An extreme gradient boost (XGBoost) machine learning model was used to predict HAB species abundance from environmental variables. The model performed strongly across much of the dataset, explaining 77% of the variation in HAB abundance, with a normalised root mean square error of 5.7%. The model’s predictive performance declined under very high observed HAB abundances (> 1.0 x 108 cells L-1), suggesting that extreme bloom events are harder to predict due to their rarity and non‑linear ecological triggers. SHapley Additive exPlanations (SHAP) values were used to quantify each environmental predictor variable’s contribution to the variation explained by the model. Salinity was revealed as the most influential predictor, aligning with results from other analyses, with higher salinity strongly suppressing HAB formation. pH was the second most important variable, exerting mostly negative effects on HAB abundance. NOX and distance from the Murray Mouth were moderate predictors linked to nutrient supply and spatial-chemical gradients, with NOX exerting a positive influence on HAB formation, while HAB abundance declined with distance from Murray Mouth. The XGBoost model was developed as a tool to better understand the broad context of harmful algal proliferation in the Coorong, but not to be used as a primary method for risk assessment. However, this model has the potential to be further developed to predict the abundances of individual species that have defined cell count risk thresholds and pose a significant threat to the Coorong.
Climate change projections for the Coorong predict continued declines in freshwater inflows, increased air and water temperatures, higher variability in droughts and floods and increasing marine influence with sea-level rise. Under future conditions, the Coorong is expected to experience more frequent and prolonged hypersalinity, particularly in the South Lagoon, resulting in a decline in low‑salinity habitats (H1–H2) that currently support higher species richness and Cyanobacteria and Chlorophyte diversity and abundance. Picophytoplankton, Diatoms and Dinoflagellate abundances are predicted to increase, with decreasing species richness as hypersaline habitats (H4–H6) expand. Reducing species richness and the decline of several taxonomic groups will impact the food web, affecting trophic transfer efficiency. Increasing picophytoplankton dominance will increase turbidity, reducing submerged macrophyte habitat with further consequences for the Coorong’s biodiversity. Increasing frequency and intensity of drought increases the vulnerability of the Coorong to toxin-producing Dinoflagellate HABs. Without sufficient freshwater inflows, the system risks declining to a permanent degraded state that resembles the end of the Millennium Drought (2007–2010).
To maintain a healthy phytoplankton community and Coorong ecosystem, the following recommendations are outlined:
1. Maintain and enhance freshwater flows to preserve low‑salinity habitats critical for phytoplankton diversity, HAB mitigation and hydrological flushing to limit picophytoplankton dominance and maintain water clarity,
2. Implement regular and continuing phytoplankton monitoring to track changes in phytoplankton communities across varying flow periods, detect early HAB signals and improve predictive modelling of phytoplankton dynamics. It is recommended that this monitoring includes testing algal-produced toxin levels,
3. Incorporate phytoplankton community responses to coupled hydrodynamic–biogeochemical-ecological models predicting salinity, water level, temperature and nutrient inputs to capture the spatiotemporal dynamics of the Coorong,
4. Refine machine learning models predicting HABs by further refining model hyperparameters and incorporating additional variables and species‑specific predictors to improve bloom detection, particularly under extreme abundances,
5. Investigate species interactions and grazing dynamics, including co‑occurrence networks, which strongly influence phytoplankton community composition, an
Strengthen HAB risk frameworks by developing species‑specific thresholds for the Coorong and integrating outputs into management response protocols.