Description: Highlights • Increasing influence of multiple environmental drivers produces changes on the temporal variability of species. • The intensity and hierarchy of drivers acting upon organisms within alternative regimes of variability may differ. • We identified regimes of variability of phytoplankton and depicted cascading effects of multiple drivers in each period. • The number of factors driving the response of phytoplankton increased along time and produced the erosion of productivity patterns. • The hierarchy and interactions of drivers changed over time, revealing that management policies require constant update. Abstract Estuaries are among the most valuable aquatic systems by their services to human welfare. However, increasing human activities at the watershed along with the pressure of climate change are fostering the co-occurrence of multiple environmental drivers, and warn of potential negative impacts on estuaries resources. At present, no clear understanding of how coastal ecosystems will respond to the non-stationary effect of multiple drivers. Here we analysed the temporal interaction among multiple environmental drivers and their changing priority on shaping phytoplankton response in the Bahía Blanca Estuary, SW Atlantic Ocean. The interaction among environmental drivers and the number of significant direct and indirect effects on chlorophyll concentration increased over time in concurrence with enhanced anthropogenic stress, changing winter climate and wind patterns. Over the period 1978–1993, proximal variables such as nutrients, water temperature and salinity, showed a dominant effect on chlorophyll, whereas in more recent years (1993–2009) climate signals (SAM and ENSO) boosted indirect effects through its influence on precipitation, wind, water temperature and turbidity. Turbidity emerged as the dominant driver of chlorophyll while in recent years acted synergistically with the concentration of dissolved nitrogen. As a result, chlorophyll concentration showed a significant negative trend and a loss of seasonal peaks reflecting a pronounced reorganisation of the phytoplankton community. We stress the need to account for the changing priority of drivers to understand, and eventually forecast, biological responses under projected scenarios of global anthropogenic change.

Description: Understanding phytoplankton species-specific responses to multiple biotic and abiotic stressors is fundamental to assess phenological and structural shifts at the community level. Here, we present the case of Thalassiosira curviseriata, a winter-blooming diatom in the Bahía Blanca Estuary, Argentina, which displayed a noticeable decrease in the past decade along with conspicuous changes in phenology. We compiled interannual field data to assess compound effects of environmental variations and grazing by the invasive copepod Eurytemora americana. The two species displayed opposite trends over the period examined. The diatom decreased toward the last years, mainly during the winters, and remained relatively constant over the other seasons, while the copepod increased toward the last years, with an occurrence restricted to winter and early spring. A quantitative assessment by structural equation modeling unveiled that the observed long-term trend of T. curviseriata resulted from the synergistic effects of environmental changes driven by water temperature, salinity, and grazing. These results suggest that the shift in the abundance distribution of T. curviseriata toward higher annual ranges of temperature and salinity—as displayed by habitat association curves—constitutes a functional response to avoid seasonal overlapping with its predator in late winters. The observed changes in the timing and abundance of the blooming species resulted in conspicuous shifts in primary production pulses. Our results provide insights on mechanistic processes shaping the phenology and structure of phytoplankton blooms.