Climate extremes intensify global lake eutrophication by increasing the stress resistance of harmful bloom-forming algae

Domain: chemistry
Published: 2026-02-17
Editorial outcome: minor_revision

Abstract
Lake eutrophication is increasingly shaped by climate extremes, but how short-lived events yield persistent blooms remains unclear. Using global satellite records and laboratory and field experiments, here we show that heatwaves and extreme precipitation drive bloom dynamics beyond gradual warming. In harmful bloom-forming algae (HBFA), heatwaves trigger oxidative stress that rapidly induces intracellular polyphosphate and stabilisomes, dense polyphosphate-rich organelles. As intracellular ballast, stabilisomes drive downward migration, enabling access to sediment-derived phosphorus while avoiding thermal stress; CO2 depletion elevates pH and reinforces a thermo-alkaline cascade amplification effect. Extreme precipitation delivers pulsed phosphorus inputs that are stored as intracellular polyphosphate, creating long-lived phosphorus reserves that prime later heatwaves. When precipitation pulses precede heatwaves, stabilisome formation, vertical migration, and bloom expansion are amplified, even in oligotrophic lakes. Compound climate extremes thus convert episodic disturbances into sustained ecological advantages, challenging nutrient-centric models and redefining bloom-risk prediction and management.

Introduction
Eutrophication has become one of the most pervasive threats to fresh water ecosystems, driven historically by excessive nutrient loading and increasingly reshaped by a rapidly changing climate1. Although nutrient enrichment has long been viewed as the dominant force governing algal proliferation, global observations now reveal that climate-related disturbances, particularly short-lived but intense events such as heatwaves and extreme precipitation, recurrently disrupt lake biogeochemistry and trigger abrupt shifts in bloom trajectories. Heatwaves can accelerate algal metabolism2, enhance nitrogen assimilation3, and favor thermophilic and alkaline-tolerant taxa by intensifying CO2 depletion and elevating water pH4. Extreme precipitation, in turn, delivers episodic pulses of nutrient-rich runoff that restructure lake trophic conditions over timescales of hours to days5,6. These climatic convulsions increasingly co-occur and appear to generate bloom trajectories that cannot be explained by gradual warming or nutrient inputs alone, suggesting the existence of rapid physiological responses that enable harmful bloom-forming algae (HBFA) to capitalize on environmental shocks.

Among the cellular strategies that may underlie this capacity, the intracellular accumulation of inorganic polyphosphate (polyP) has emerged as a candidate mechanism linking climatic perturbations to algal resilience. PolyP buffers oxidative stress, regulates energy balance, and mediates phosphorus storage under fluctuating environmental conditions7,8. It is sequestered within dense acidic organelles known as stabilisomes, which have been associated with cellular homeostasis and stress adaptation9. Yet despite their ubiquity across bloom-forming algae, the ecological significance of polyP and stabilisomes remains poorly resolved, particularly in the context of extreme climatic events. Whether heatwaves induce stabilisome formation that alters algal buoyancy and nutrient acquisition, and whether precipitation-driven phosphorus pulses generate a biochemical memory in phosphorus storage that primes HBFA for subsequent heat stress, are questions that remain unanswered. More broadly, it is unknown how the climate extremes interact with polyP-mediated processes to regulate bloom persistence across lakes of different trophic status.

These gaps highlight a fundamental challenge for contemporary eutrophication research: the mechanisms by which short-term climatic shocks are transformed into sustained ecological advantages for HBFA remain unresolved, limiting our ability to forecast bloom risk under future climate scenarios. Current conceptual frameworks predominantly emphasize nutrient supply or long-term warming and treat extreme events as transient disturbances. However, emerging evidence suggests that climate extremes may act as catalytic triggers that activate physiological pathways4, potentially involving polyP accumulation and stabilisome formation, that restructure algal behavior, stress tolerance, and vertical mobility on timescales much longer than the disturbances themselves. A mechanistic, cross-scale understanding of these processes is urgently needed. Addressing these unresolved questions requires a framework that links short-term climatic shocks to the physiological pathways that sustain harmful algal persistence.

Here, we combine two decades of global satellite observations with controlled laboratory assays and field-derived incubations to quantify how heatwaves and extreme precipitation regulate polyP metabolism, stabilisome formation, and vertical migration in HBFA. We show that heatwaves induce oxidative stress that triggers rapid polyP accumulation and the formation of dense stabilisomes that function as intracellular ballast, enhancing downward migration, thermal avoidance, and access to sediment-derived phosphorus. We further demonstrate that precipitation-driven phosphorus pulses are rapidly assimilated and stored as polyP, creating a biochemical memory in phosphorus storage that primes HBFA for subsequent heatwaves. Together, these mechanisms reveal how compound climate extremes convert transient disturbances into sustained physiological and ecological advantages. These findings provide a mechanistic framework that reconciles global eutrophication trends with climate extremes and highlight the need to incorporate compound-event dynamics into bloom prediction and lake-management strategies.

Results
Climate extremes reshape the global trajectory of lake eutrophication

As shown in Fig. 1a, b, global lake algal bloom scale, indexed by global median bloom frequency, exhibited distinctly irregular fluctuations from 2003 to 2022 across more than 600 lakes worldwide with bloom occurrences persisting for over a decade. Despite this variability, its long-term increase aligns with the rising trajectories of major anthropogenic and climatic drivers, including population growth, surface warming, and intensifying precipitation (Fig. 1b). Although ~47% of global lake degradation has been attributed to intensified land use and nutrient emissions10, such gradually increasing pressures alone cannot account for the pronounced interannual variability in algal bloom. Consistent with this discrepancy, previous machine-learning analyses identified temperature anomalies and extreme precipitation as the dominant predictors of future algal bloom trajectories11, suggesting that short-lived climatic perturbations exert disproportionately strong influence. Extending these findings, our generalized additive model (GAM) analysis further shows that temperature anomalies and extreme precipitation metrics retain consistently high importance across alternative model formulations, whereas gradual climatic indicators exhibit greater sensitivity to model structure (Fig. 1c). Together with the documented rise in compound lake-extreme events12, these findings indicate that episodic climatic disturbances rather than long-term mean conditions primarily govern the interannual variability and amplification of contemporary eutrophication.

Although climate extremes dominate sharp interannual fluctuations in global median bloom frequency, long-term demographic pressure and gradual climate change continue to shape the baseline trophic state of lakes. Model-based importance analysis indicates that human population, temperature, and precipitation remain non-negligible contributors to global chlorophyll-a (Chl-a) variability (Fig. 1c), supporting the role of long-term drivers in shaping the baseline eutrophication trajectory13. Consistent with this pattern, analyses of 547 globally distributed bloom-prone lakes showed that Chl-a increased with watershed population density and with aqueous total nitrogen (TN) and total phosphorus (TP) concentrations (Fig. S1a), reflecting the foundational role of nutrient enrichment in supporting algal blooms14.

Chl-a concentrations are generally projected to rise under future climate scenarios, with nutrients acting as the dominant amplifiers in landscapes undergoing intensified land use11. Nevertheless, reductions in watershed nutrient inputs can suppress algal bloom, demonstrating that eutrophication remains responsive to nutrient management even under warming15. Yet this expectation does not always hold. In China, for example, harmful algal blooms have intensified despite declining external nutrient loads16, suggesting that climatic disturbances may interact with large internal legacy nutrient pools to push systems across ecological thresholds17. The nonlinear relationships between Chl-a and either temperature or precipitation (Fig. S1a) further support this view, indicating that climatic forcing modulates bloom outcomes through its interactions with internal nutrient dynamics, and that the strength of these interactions varies across management regimes.

Although global human population and Chl-a exhibit broadly synchronous upward trends, multivariable importance analyses indicate that hydrometeorological extremes consistently rank among the most influential predictors of Chl-a variability once climate and anthropogenic drivers are considered together (Fig. 1c). Specifically, temperature anomalies and multiple indices of extreme precipitation consistently explain a larger share of Chl-a variability and retain stable importance across alternative model formulations, indicating that lake algal blooms are more sensitive to short-term climatic perturbations. These results support a coherent conclusion that, under contemporary climate conditions, hydrometeorological extremes exert a more direct and stronger control on bloom expansion than do gradual anthropogenic and climatic changes. Independent evidence corroborates this pattern: aquatic heatwaves alone can increase surface Chl-a by more than 20% and occur alongside bloom events at significantly above-random frequencies18. Event-attribution analyses further reveal that 24 of 29 documented summer lake heatwaves coincided with pronounced Chl-a increases, with mean enhancements exceeding 50% in low-nutrient lakes19. Together, these findings indicate that hydrometeorological extremes function as threshold-like disturbances that rapidly translate accumulated nutrient legacies into elevated algal biomass, generating eutrophication risks that exceed those imposed by gradual warming alone.

In addition to governing Chl-a variability, climate extremes also reshape nutrient availability, thereby influencing the resource environment that supports algal bloom development. To evaluate nutrient-response pathways, TN and TP dynamics were modeled across 547 lakes using GAMs (Fig. S1b). Strong positive associations with watershed population density reflected sustained anthropogenic nutrient inputs, whereas nonlinear dependencies on temperature and precipitation revealed an additional climate-sensitive mechanism (Fig. S1b). Storm events can abruptly increase external loading through runoff and sediment resuspension, especially in shallow or disturbed catchments20,21. At event-dominated timescales, stormflow pipelines deliver high-flux, short-duration nutrient pulses to lake inlets, while sediment resuspension mobilizes benthic nutrient reserves into bioavailable forms, producing bloom peaks that frequently lag rainfall18. Thus, climate extremes can transiently elevate nutrient availability for HBFA through pulsed inputs even without changes in external anthropogenic loads.

It should be noted that the bloom-frequency metric used here is derived from surface-sensitive satellite observations and therefore primarily captures surface or near-surface algal accumulations, potentially underrepresenting subsurface biomass in some lake systems. Moreover, while the use of a single, internally consistent MODIS–Aqua sensor and a validated detection framework minimizes cross-sensor and classification uncertainties, lake-specific hydrodynamics and internal nutrient cycling processes are not explicitly resolved in this global-scale analysis.

Collectively, these results demonstrate that contemporary eutrophication risk is governed by the interaction between long-term nutrient enrichment and short-duration climate extremes. Rather than diminishing the importance of anthropogenic drivers, hydrometeorological extremes act as pulse-like amplifiers that rapidly convert pre-existing nutrient legacies into bloom expansion. This pulse-sensitive coupling between episodic climatic disturbances and accumulated nutrient reservoirs defines a central feature of modern lake eutrophication dynamics, underscoring the need to complement mean-load reduction strategies with management approaches that explicitly address climate extremes and nutrient pulse regimes. Recent evidence showing that nutrient-pulse regimes govern phytoplankton functional stability provides further mechanistic support for shifting from mean-load reduction to pulse-oriented management strategies12,17,22. The mechanistic pathways through which climate extremes intensify eutrophication, including their effects on nutrient pulses and HBFA stress physiology, are examined in the following sections.

Heatwaves promote bloom persistence by enhancing the stress resistance of HBFA

Heatwaves impose pronounced oxidative stress on HBFA while simultaneously triggering physiological priming that enhances their capacity for survival and regrowth. In our experiments, exposure to 40°C increased reactive oxygen species (ROS; O2⁻, H2O2) by 35.9–215.8%, elevated antioxidant enzyme activities (superoxide dismutase [SOD] and catalase [CAT]) by 0.55–6.30 fold, and raised malondialdehyde (MDA) levels by ~0.81–1.96 fold relative to 25°C (Fig. 2a). Photosynthetic performance declined in parallel, as indicated by reductions in effective quantum yield of PSII (Y[II]), maximum quantum yield of PSII (Fv/Fm), and maximum electron transport rate (ETRmax) (Fig. S2a). These responses mirror those observed in Scenedesmus obliquus and Chlorella under heat stress23, and demonstrate that HBFA engage heat-stress priming mechanisms that contribute to the establishment of thermal memory24. The transient fall in Fv/Fm is consistent with accelerated PSII damage and impaired D1 repair25. However, cell density still increased at 40°C, suggesting activation of heat-resilient metabolic pathways, including heterotrophic carbon use and nutrient remobilization26.

Recovery assays showed that photosynthetic performance and growth were rapidly restored when the temperature returned to 25°C (Fig. S2b), demonstrating pronounced thermal plasticity. This response parallels findings in Pyropia yezoensis, where non-lethal heat priming establishes heat-stress memory and improves post-stress survival27. The enhanced tolerance in this red alga is partly attributed to memory-driven induction of nutrient transporter genes, which promotes nutrient acquisition and post-stress growth28. Consistent with these priming phenomena, HBFA in our study exhibited marked increases in cellular protectants: heat shock proteins rose (HSP) by 47.5–66.2% and polyP by 57.9–197.0% at 40°C relative to 25°C (Figs. 2a and 3a). Although HSP stabilize protein complexes during thermal stress, polyP provides the central metabolic and redox buffering capacity that underpin sustained heat tolerance29. Acting as a dynamic energy reservoir, polyP supports ATP regeneration and antioxidant defenses30, and helps maintain ATP/ADP balance under fluctuating thermal conditions31. Heatwaves also promoted nutrient acquisition, as intracellular TP and dissolved inorganic phosphorus (DIP) increased by 18.4–86.5% and 36.7–93.1%, respectively, at 40°C (Fig. 3b), consistent with heat-induced luxury phosphorus uptake through high-affinity transport systems32. Together, these responses indicate that heatwaves prime HBFA for enhanced phosphorus acquisition and structural stabilization, enabling rapid recovery once thermal stress subsides.

Bloom intensification during heatwaves substantially enhances photosynthetic CO2 uptake, which in turn elevates water pH, often exceeding 9 in eutrophic lakes33. To mechanistically resolve how HBFA withstand such thermo-alkaline conditions, we assessed cellular responses under controlled alkaline stress. Stabilisomes, which are polyP-rich, high-density granules that function both as intracellular phosphorus reservoirs and as biomechanical ballast structures, were examined as a potential mediator of resilience. Exposure to pH 10 led to sharp increases in intracellular TP, DIP, and polyP (27.1–209.1%, 81.6–156.4%, and 98.5–153.5% relative to pH 7), together with elevated activities of SOD and CAT and higher MDA levels (Figs. 2b and 3d). These responses indicate that alkalinity stimulates luxury phosphorus uptake, providing abundant substrate for polyP synthesis and thereby strengthening intracellular phosphorus storage and metabolic robustness34. Consistent with this pattern, Microcystis aeruginosa exhibits heightened ROS production, lipid peroxidation, and antioxidant activity at high pH 935. Taken together, these findings show that polyP plays dual protective roles under alkaline stress by buffering oxidative damage and maintaining intracellular phosphate homeostasis, thereby stabilizing cellular metabolism36,37. PolyP accumulation therefore emerges as a central mechanism enabling HBFA to preserve metabolic integrity and withstand oxidative pressure in alkaline environments.

The alkaline responses observed in HBFA are mirrored across diverse bloom-forming taxa and natural systems, indicating that thermo-alkaline coupling is broadly conserved. In Aphanizomenon flos-aquae, elevated pH (pH 8) activates alkaline phosphatase (APase), enabling the mobilization of organic phosphorus that sustains growth under external phosphorus limitation38. Cyanobacterial mats show a similar response, in which high APase activity is coupled with coordinated polyP accumulation to maintain growth during fluctuating phosphate availability39,40. Field observations from Lake Taihu further support this coupling: during bloom seasons, increases in water temperature (26–35°C), pH (8.6–9.2), and intracellular phosphorus content occurred simultaneously (Fig. S3c), indicating that thermo-alkaline stress is closely linked to enhanced phosphorus assimilation. Comparable alkalinization has been reported across diverse lake systems, including Lake Erie (Microcystis; pH > 9.5, episodically ~11), Upper Klamath Lake (Aphanizomenon; pH > 9), and Lake Bogoria (Arthrospira; ~10.5)41,42. Collectively, these patterns suggest that elevated pH is an intrinsic metabolic feedback of bloom development rather than an externally imposed condition.

Having established the ubiquity of thermo-alkaline coupling, we next asked whether these responses exhibit persistence or priming effects. PolyP continued to accumulate even after temperature had stabilized (Fig. S3d), and algae pre-exposed to pH 10 recovered more rapidly than those pre-conditioned at pH 8 (Fig. S2c). These patterns indicate that alkaline pre-conditioning can establish a primed state analogous to thermal acclimation. They also support the view that heatwaves and elevated pH interact to initiate a cascade of polyP-mediated stress responses in HBFA. By simultaneously enhancing phosphorus retention and oxidative stress tolerance, these responses generate a thermo-alkaline amplification effect that reinforces HBFA persistence. The resulting positive feedback loop strengthens bloom resilience and may accelerate eutrophication under increasingly frequent climate extremes.

While polyP accumulation under thermo-alkaline stress enhances metabolic stability, bloom persistence requires additional ecological strategies (e.g., buoyancy-driven vertical migration) that translate intracellular resilience into improved opportunities for nutrient acquisition and stress avoidance. Microcystis and other buoyant taxa modulate intracellular carbohydrate ballast and gas vesicles to adjust buoyancy and descend under adverse conditions32,43, a metabolically regulated avoidance behavior that promotes bloom stability and expansion44. Rising temperatures further accelerate buoyancy regulation and vertical migration in Microcystis, facilitating rapid redistribution within the water column45. Beyond providing temporary refuge from photoinhibition and oxidative damage, this vertical repositioning enables access to sediment-derived phosphorus, thereby reinforcing bloom persistence.

In the context of these thermo-alkaline and nutrient-mediated responses, our results identify stabilisome formation as an additional and critical pathway underpinning HBFA resilience. Microscopy revealed that stabilisomes are ubiquitous among bloom-forming aglae and accumulate prominently in actively growing cells exposed to either 40°C or pH 10, conditions under which cell density and downward migration were consistently enhanced (Figs. 3 and 4). Functionally, stabilisomes act as high-density intracellular ballast that promote sinking, as their bulk density (1.79–3.13 g cm−3) far exceeds that of major cellular constituents, including proteins (1.3–1.43 g cm−3), nucleic acids (2 g cm−3), carbohydrates (1.55–1.62 g cm−3), and lipids (0.91–1.01 g cm−3)46. Vertical migration assays corroborated this role, showing that heat-induced stabilisome proliferation across both unicellular and filamentous HBFA resulted in substantially greater algal accumulation at depths of 15–30 cm relative to surface layers (Fig. 4a, b). Under non-stress conditions, polyP stored within stabilisomes can be rapidly mobilized as a preferential phosphorus source to support macromolecular synthesis and restore buoyancy7, enabling reversible control of cellular density. This dynamic regulation is especially pronounced in small unicellular or colonial taxa, whereas vertical migration in larger taxa is further shaped by colony size, aggregation, and the balance between buoyancy and drag forces43,47.

Thus, stabilisome accumulation complements the canonical carbohydrate–gas vesicle buoyancy system by providing a biomechanical means for cells to alternate between deep, stress-protected refuges and light-rich surface waters. This stress-induced, metabolically powered density cycling tightly couples nutrient acquisition with stress avoidance, positioning stabilisomes as both biochemical reserves and dynamic regulators of vertical mobility that enhance HBFA growth. By enabling rapid redistribution within the water column, this mechanism confers resilience under thermal and chemical extremes and reinforces the positive feedbacks characteristic of thermo-alkaline amplification.

Extreme precipitation increases phosphorus pulse inputs and facilitates rapid HBFA assimilation

In addition to thermo-alkaline amplification, bloom initiation and persistence depend critically on short-lived nutrient pulses generated by hydrological extremes. Intense precipitation mobilizes particle-bound phosphorus through watershed runoff and sediment disturbance, producing high-concentration inputs that can substantially supplement algal nutrition48. Because such events typically occur during warm seasons when sedimentary phosphorus desorption is temperature sensitive, precipitation-driven resuspension and thermal conditions often co-occur in ways that jointly influence nutrient mobilization. Laboratory assays further showed that DIP and TP released from suspended solids and sediments increased by ~1.57–9.43 fold at 25–35°C relative to 15°C (Fig. 5b, c), indicating that precipitation-driven mixing and temperature-enhanced desorption can act synergistically to accelerate phosphorus release. Consistent with this mechanism, declines in precipitation intensity were accompanied by reduced Chl-a concentrations, reflecting diminished external nutrient delivery (Fig. S1a). Similar hydrologically induced nutrient mobilization has been reported in the North American Great Lakes, where disturbed particulate phosphorus constitutes a major source of phytoplankton nutrition49.

Once mobilized into the water column, this storm-driven phosphorus is rapidly assimilated by HBFA. Within 4 h at elevated temperature, DIP uptake increased by 201.1–2,011.1% relative to initial cellular content, far exceeding immediate metabolic requirements, whereas TP increased by 142.2–868.0% (Fig. 5d, e). This disparity reflects classic luxury uptake, whereby excess phosphorus is sequestered intracellularly as polyP, transforming brief hydrological inputs into durable biochemical reserves48,50. Field and experimental evidence further indicates that storm-mobilized particulate phosphorus represents a major component of the bioavailable pool and is rapidly intercepted by algae before it can resettle50–52. Thus, precipitation pulses not only trigger bloom onset but also facilitate metabolic “investment” into polyP reserves that enhance resilience to subsequent stressors7.

The ability to convert pulsed phosphorus inputs into polyP stores is not unique to HBFA but reflects a deeply conserved algal strategy. PolyP metabolism is widely conserved across algal groups and underpins resilience to diverse environmental disturbances. In cyanobacteria, polyP contributes to energy regulation and stress survival34, whereas in Karenia mikimotoi and Achnanthidium minutissimum, both polyP content and gene expression respond sharply to external phosphorus availability53,54. Even picoplankton in oligotrophic systems sustain high productivity through efficient polyP accumulation and recycling55. This conserved strategy is mirrored during extreme precipitation, when storm-driven phosphorus surpluses are rapidly assimilated and converted into polyP, effectively extending the pulse memory of short-lived hydrological inputs. Collectively, these findings identify polyP accumulation as a universal mechanism that enables HBFA to convert transient nutrient pulses into a sustained biochemical capacity for coping with later stress.

Once formed during nutrient pulses, these polyP reserves continue to support bloom development even in the absence of new external inputs. Beyond intracellular granules, polyP is stored across multiple ecological compartments, including cyanobacterial mats, sediment–water interface layers, and the extracellular polymeric substances of large colonies, forming a distributed internal phosphorus reservoir that buffers cells against fluctuations in external supply. Cyanobacterial mats, for example, rapidly accumulate polyP and maintain stoichiometric flexibility under shifting phosphate regimes42, while biogenic polyP in sediments contributes to seasonal phosphorus release that remains bioavailable during subsequent stratification or resuspension events56. These internal-storage pathways help explain why load-reduction programs often yield delayed or muted improvements in water quality, because internal polyP pools can continue to support biomass persistence even as external inputs decline57. Extreme precipitation further reinforces this internal reservoir by mobilizing particle-bound phosphorus that is rapidly intercepted and incorporated into polyP before resettling, effectively renewing the internal pool on each hydrological cycle. Consequently, polyP acts not only as a biochemical reserve that stabilizes HBFA under heat and alkaline stress, but also as a mechanism that converts episodic hydrological disturbances into long-lived metabolic reinforcement.

This reservoir function of polyP is not limited to eutrophic systems. In oligotrophic lakes, even modest precipitation-induced increases in phosphorus can reach thresholds sufficient for HBFA uptake. Under such transient enrichment, HBFA activate APase and mobilize organic phosphorus to meet their nutritional demand58. PolyP pathways frequently operate in parallel with these enzymatic mechanisms, providing dual safeguards against nutrient fluctuation37,59. Consequently, luxury phosphorus uptake and polyP storage confer adaptive advantages in both eutrophic and nutrient-poor systems. Despite this broad functional role, further work is needed to resolve how thermal stress, hydrological variability, and cellular physiological adjustments interact to support bloom persistence and intensification under extreme climatic conditions.

A unified framework linking climate extremes to bloom persistence

Our findings demonstrate that short-lived, high-intensity climate extremes exert disproportionate control over the intensification and persistence of HBFA. To integrate these processes, we propose a unified conceptual framework (Fig. S5), the thermo-alkaline cascade amplification model, in which cellular and ecological adjustments enable HBFA to convert episodic environmental shocks into sustained ecological advantage. Heatwaves impose strong oxidative stress that activates antioxidant defenses, stimulates phosphorus uptake, and drives substantial polyP accumulation. In parallel, HBFA deploys a stabilisome-based ballast system that facilitates downward migration into cooler, deeper waters where cells can continue proliferating using previously stored luxury phosphorus. As bloom biomass increases, enhanced photosynthetic CO2 drawdown elevates water pH, initiating a secondary alkaline loop that further enhances buoyancy regulation and intracellular polyP storage. Together, these responses recast heatwaves and elevated pH not simply as stressors but as ecological cues that activate resource-conserving strategies, positioning stress resistance as a central mechanism linking climate extremes to bloom persistence.

Extreme precipitation further strengthens this cascade by rapidly increasing the pool of bioavailable phosphorus for immediate assimilation. The resulting pulses enhance polyP storage, antioxidant capacity, and metabolic flexibility, providing resilience during subsequent thermal stress. Consistent with these patterns, incorporating precipitation into temperature-based GAMs increased the explanatory strength of temperature effects on Chl-a across both trophic states. Specifically, the deviance explained by temperature increased from 0.138 to 0.153 in oligotrophic lakes and from 0.093 to 0.146 in eutrophic lakes when precipitation was included (Fig. S4a–d), indicating that precipitation enhances temperature-driven Chl-a variability rather than acting as an independent driver.

When storm-driven phosphorus pulses recur in sequence with heatwaves, their effects compound, as phosphorus accumulated during precipitation events becomes a critical internal reserve during later heat stress, supporting continued proliferation and enabling vertical redistribution. This climate-rhythm coupling sustains bloom persistence even in oligotrophic lakes, illustrating how HBFA convert episodic environmental shocks into long-lived ecological advantage.

Collectively, our results show that episodic climate extremes mobilize nutrient pulses, accelerate their biological capture, and initiate self-reinforcing feedback loops that enhance algal resilience far more rapidly than gradual climate trends. When anthropogenic influences are removed, temperature anomalies and precipitation extremes exert particularly strong control over lake Chl-a variability, highlighting the acute sensitivity of bloom trajectories to short-term climatic perturbations. Heatwaves induce oxidative stress that stimulates polyP accumulation and antioxidant activation, while bloom-driven CO2 depletion elevates pH, forming a thermo-alkaline amplification cycle in which stabilisome accumulation accelerates downward migration and deep-water refuge. Extreme precipitation, by contrast, mobilizes sediment- and particle-bound phosphorus that is rapidly assimilated before resettling, and when such pulses coincide with heatwaves, nutrient-retention and thermal-priming mechanisms act jointly to expand blooms beyond what either stressor could induce alone (Fig. 6). These integrated mechanisms explain how climate extremes intensify global eutrophication by enabling HBFA to convert transient shocks into sustained biomass advantage. As climate extremes become more frequent, bloom trajectories will increasingly be governed by the timing, intensity, and co-occurrence of heatwaves and hydrological pulses rather than by gradual warming or mean precipitation alone.

Methods
Data sources

In this study, global lake algal bloom dynamics were characterized using a recently published global lake algal bloom dataset developed by Wang et al.60 (https://doi.org/10.1093/nsr/nwaf011). The dataset is derived from the complete MODIS–Aqua archive spanning 2003–2022 and comprises ~0.8 million satellite images with a temporal resolution of 1–2 days. It covers 1956 large freshwater lakes with surface areas greater than 50 km2, representing ~51% of the global freshwater lake area. For the GAM analyses, we further focused on algal bloom frequency in 607 large shallow lakes worldwide that experienced persistent blooms for more than 10 years during 2003–2022.

Algal bloom detection is based on Rayleigh-corrected reflectance (Rrc) and a normalized Floating Algal Index (nFAI), which is specifically optimized for detecting surface and near-surface algal accumulations. To reduce non-algal interference, CIE color space masking was applied to distinguish algal blooms from suspended sediments, and floating-leaved macrophytes were excluded using the Global Surface Water Occurrence dataset. This multi-step processing framework effectively minimizes misclassification arising from atmospheric correction errors, cloud contamination, sun glint, and shoreline mixed-pixel effects.

The bloom detection algorithm has been systematically validated across lakes from multiple climatic and geographic regions using independent visual interpretation and in situ observations, achieving an overall classification accuracy of ~81%. Moreover, the dataset is constructed from a single, internally consistent sensor (MODIS–Aqua), thereby avoiding cross-sensor inconsistencies and temporal discontinuities that may affect multi-sensor merged Chl-a products.

In this study, annual algal bloom frequency at the lake scale was defined as the ratio of bloom-positive observations to the total number of valid satellite observations within a given year. Lake-level bloom frequency metrics were subsequently aggregated to the global scale and integrated with climate extreme indices and socioeconomic variables for subsequent statistical analyses.

To complement the satellite observations and provide partial ground-truthing, we also incorporated a long-term in situ dataset comprising 547 lakes worldwide, originally compiled by Qin et al.61. Although widely used for ecological assessments, this dataset has not previously been examined in the context of population density, temperature, or precipitation. The integration of satellite and in situ records, therefore, provides a more robust and representative basis for assessing global eutrophication dynamics and their climatic and demographic drivers.

Annual anomalies in global land temperature (2003–2022), global surface air temperature increases relative to the 1850–1900 pre-industrial baseline, and global daily precipitation were obtained from the ERA5 reanalysis dataset produced by the European Center for Medium-Range Weather Forecasts and routinely disseminated through the Copernicus Climate Change Service. These datasets enabled the large-scale analysis of temperature- and precipitation-driven variation in algal bloom dynamics. From the daily precipitation fields, we computed three widely used indices of extreme rainfall: Rx1day, the annual maximum 1-day precipitation amount; R95p, the annual total precipitation from days exceeding the 95th percentile of the reference period; and PTOT, the cumulative precipitation associated with extreme events above the 95th percentile, which serves as a proxy for cumulative intense rainfall62. Global human population data were obtained from the United Nations World Population Prospects (https://population.un.org/wpp/). Lake-specific population density was derived from the WorldPop global population dataset (https://www.worldpop.org/) in combination with the GeoNames geographical database (http://www.geonames.org/), providing spatially explicit demographic information for watershed-scale analyses.

Preparation of HBFA and measurement of HBFA stress responses

The resistance mechanisms of HBFA under climate extremes were investigated using M. aeruginosa, Raphidiopsis raciborskii, Dolichospermum flos-aquae, Aphanizomenon flos-aquae, all of which are globally dominant harmful bloom-forming algae. Cultures were obtained from the Institute of Aquatic Biology, Chinese Academy of Sciences, and maintained in BG11 medium at 25°C under a 12 h light/12 h dark photoperiod with a light intensity of 33 μmol photons m−2s−1. Cells were agitated three times daily and harvested during exponential growth by centrifugation (~1790×g, 10 min). The pellet was resuspended in 2 L of BG11 medium, homogenized, and evenly distributed into fifteen 250-mL Erlenmeyer flasks. For temperature treatments, cultures were incubated at 25, 30, 35, 40, and 45°C at a constant pH of 7. For pH treatments, cultures were maintained at 25°C and adjusted to pH 7, 8, 9, or 10 using dilute HCl or NaOH prepared in BG11. All incubations were performed under a 12 h light/12 h dark cycle (33 μmol photons m−2s−1) with shaking at 150 rpm for 4 h per day to ensure culture homogeneity.

At the end of each experiment, all four HBFA species were analyzed for a suite of physiological and biochemical parameters, including phosphorus fractions, superoxide anion and hydrogen peroxide levels, antioxidant enzyme activities, heat shock protein abundance, photosynthetic performance, cell density, and activities of growth-related enzymes. To evaluate recovery from heat stress, cultures exposed to 35°C and 40°C were subsequently transferred to 25°C and maintained for an additional 4 days under the same conditions as the 25°C control. Likewise, after 4 h of alkaline exposure at pH 9 or 10, cultures were adjusted to pH 8 and incubated for 4 days under the conditions used for the pH 8 control. MDA content, SOD and CAT activities, Chl-a concentration, and cell density were quantified on Days 1 and 4 of the recovery phase.

To assess whether stabilisomes are present in natural HBFA populations, in situ bloom samples were collected during active bloom events from Lake Chaohu, Lake Hulun, Lake Tianlai, Lake Dianchi, Lake Taihu, and Lake Khanka (Fig. S3a, b, Table S1). Stabilisomes in all HBFA samples were visualized using a modified Albert staining procedure following the manufacturer’s protocol (Cat. RS4253, Lot 20250917, G-CLONE, Beijing, China).

Stabilisome-mediated vertical migration was examined using axenic cultures of A. flos-aquae, D. flos-aquae, R. raciborskii, and M. aeruginosa. Before the vertical migration assays, cultures were ultrasonicated for 1 min at 60% amplitude to disperse filaments into single trichomes and to break Microcystis colonies into unicells; dispersion was confirmed microscopically. Initial cell densities were standardized to 2.0×106, 2.4×106, 2.3×106, and 2.5×106 cells mL−1 for the four strains, respectively. Algal suspensions were transferred into cylindrical glass columns (10 cm diameter) wrapped in black paper to eliminate lateral illumination. Each column was fitted with a central vertical rod (0.8 cm diameter) holding three silicone sampling ports (3 mm diameter) positioned at depths of 2, 15, and 30 cm. The ports remained clamped until sampling to prevent unintended exchange. Columns were then incubated at either 25°C or 40°C for 1 h (n=5 biological replicates per treatment). Following incubation, 5-mL subsamples were gently withdrawn from each depth using sterile syringes, and algal cell densities were quantified.

Field in situ simulation experiment

Algal bloom material was collected from Lake Taihu using a 64-μm mesh phytoplankton net to concentrate cells on 8 September 2022. In situ water quality and ambient conditions during sampling were TP=0.076 mg L−1, TN=1.03 mg L−1, DO=8.1 mg L−1, with air temperature, light intensity, and UV radiation ranging from 26–35°C, 0–530 μmol photons m−2s−1, and 0–45.9 μW cm−2, respectively. A 20-mL aliquot of the concentrated suspension was placed between two 0.45-μm filter membranes, secured between acrylic plates, and clamped to maintain structural integrity. Six time-series groups were prepared (n=5 replicates per group), and samples were harvested at 0, 4, 8, 12, 16, and 20 h. At each time point, the material was divided into four equal subsamples and transferred into 50-mL centrifuge tubes (Tubes 1–4). Tube 1 was fixed with 15 mL of formaldehyde for algal cell density determination; Tube 2 was fixed with 15 mL of trichloroacetic acid for intracellular DIP measurement; Tube 3 was treated with 15 mL of alkaline sodium hypochlorite for intracellular polyP determination; and Tube 4 was left unfixed for subsequent intracellular TP analysis. All fixed samples were stored at 4°C in insulated containers and transported to the laboratory for chemical and microscopic analyses.

Surface sediments (0–10 cm) were collected from Lake Taihu using a Petersen grab sampler (HNM1-2), sealed immediately, kept in the dark, and transported to the laboratory, where coarse debris (e.g., branches, leaves, shell fragments) was manually removed. Suspended solids were obtained by sampling lake water 10 cm below the surface, passing it through a 200-mesh sieve to remove large particulates, and retaining the resulting suspension. Additional algal material from Lake Taihu was collected, harvested by centrifugation (50-mL tubes, ~1010×g, 15 min), washed three times with phosphorus-free BG11 medium, and resuspended in 500 mL of phosphorus-free BG11 to induce phosphorus starvation. Cultures were maintained for 8 days under standardized conditions (25°C, 41 μmol photons m–2s−1, 14 h light/10 h dark cycle, pH 8). The experiment comprised two treatment groups, sediment and suspended solids, each with three temperature regimes (15, 25, and 35°C) and three biological replicates. For the sediment treatment, 100 mL of sediment was added to 800 mL of phosphorus-free BG11 medium, whereas for the suspended-solids treatment, 100 mL of collected particulates was added to 800 mL of medium. All flasks were equilibrated to the designated temperature. Subsequently, 100 mL of the phosphorus-starved algal culture (1×108 cells mL−1), enclosed in dialysis bags with a molecular weight cutoff of 100,000 Da, was introduced into each flask (Fig. 5a). During the 250-min incubation, water samples were collected every 25 min. At each time point, 8 mL of filtrate and 6 mL of algal suspension were withdrawn to quantify TP in the overlying water and in algal cells, respectively.

Measurement of water quality parameters and HBFA

Water temperature (T), dissolved oxygen (DO), and pH were measured in situ using a HACH HQ30-D multiparameter probe (Hach Corporation, Loveland, CO, USA). TP and DIP in water were quantified colorimetrically, with or without persulfate (K2S2O8) digestion, following the protocol of Ebina et al.63. Intracellular TP and phosphate in algal cells were determined according to Zhou et al.64. Detailed step-by-step protocols for TP/DIP and intracellular TP/DIP measurements are provided in Supplementary Information. Absolute polyP was quantified by mixing 20 mL of filtered algal suspension with 2 mL of toluidine blue O (TBO) assay solution (15 mg TBO mL−1 in 0.1 N acetic acid) and incubating the mixture for 15 min at room temperature. A 1-mL subsample was subsequently hydrolyzed in 2 N HCl at 95°C for 40 min, after which the released phosphate was measured using the molybdenum blue method65. For intracellular assays, a 3.0-mL aliquot of algal culture was centrifuged at ~16,100×g for 5 min at 4°C. The pellet was resuspended in 3.0 mL of phosphate-buffered saline (PBS, 0.05 mol L−1, pH 6.8), gently mixed, and stored at −80°C. Prior to analysis, samples underwent four freeze–thaw cycles at 25°C to ensure complete cell lysis.

MDA content and the activities of SOD and CAT were quantified using commercial assay kits (Nanjing Jiancheng Bioengineering Institute, China; A003-1-2, A001-3-2, and A007-1-1, respectively) with absorbance read in 96-well plates. Detailed experimental procedures, calculation formulas, and quality controls are provided in Supplementary Information. Photosynthetic parameters, including Y(II), Fv/Fm, and ETRmax, were measured using a Phyto-PAM fluorometer (Heinz Walz, Germany). Chl-a was extracted with 90% acetone and quantified spectrophotometrically; the full extraction workflow and equations (Strickland and Parsons66) are described in Supplementary Information.

Stabilisome size and abundance were quantified by nanoparticle tracking analysis using a NanoSight NS300 (Malvern Panalytical, UK). Samples were diluted in PBS, injected into the instrument chamber, and analyzed based on Brownian motion trajectories to derive hydrodynamic diameters and particle concentrations (n=5). Dry-mass concentration was determined by freeze-drying 30 mL of the extracted stabilisome suspension to constant weight and measuring the residual mass with an analytical balance. Particle density was calculated from the total dry mass, particle size, and particle number, assuming spherical geometry for individual particles (V=4/3πr3).

The influence of stabilisome accumulation on algal vertical migration was assessed using M. aeruginosa. Two 5-mL cultures (2×106 cells mL−1) were exposed to either 40°C or 25°C for 4 h and centrifuged at their respective temperatures (~1010×g, 5 min). Pellets were gently layered onto the surface of vertically oriented glass columns (25 cm length, 4 cm diameter) filled with distilled water, and sinking time to defined depths was recorded.

Data treatment and statistical analysis

For global-scale analyses, lake Chl-a, nutrient concentrations (TN and TP), and all anthropogenic and climatic predictors were log10-transformed prior to visualization and modeling to ensure dimensional consistency and reduce scale-dependent bias. Statistical significance was assessed at p < 0.05.

Nonlinear relationships between algal bloom metrics and environmental drivers were quantified using GAMs fitted with penalized smooth functions, as detailed in the Supplementary Methods. Model diagnostics, including deviance explained, effective degrees of freedom, F statistics, p values, and sample size, are reported in the corresponding panels.

All statistical analyses and figure generation were performed in R (version 4.4.2; R Core Team, 2023), including GAM-based response curves, variable–response relationships, and all forms of line plots (e.g., temporal trends, gradient responses, and recovery trajectories). All GAM analyses were implemented in R using the mgcv package (v1.9-4). Solid lines represent fitted relationships or mean trends, and shaded regions denote 95% confidence intervals (Supplementary Code 1).

Peer review rounds
[{"round": 1, "reviewer_comments": {"reviewer_1": "The manuscript presents an interesting integration of global satellite data and laboratory experiments to explain how climate extremes may intensify lake eutrophication through algal stress resistance mechanisms. The hypothesis is compelling, and the combination of large-scale and mechanistic approaches is commendable. My comments are included below.\nI suggest including a more methodological description and adding an uncertainty evaluation (data preprocessing, model diagnostics, cross-validation).\nLimitations should be discussed more:\nThe satellite dataset (GlobColour, 2002–2024) has uncertainties, but the manuscript does not mention atmospheric correction errors, cloud contamination, or sensor variability across sensors.\nFrom what I understand, the authors conducted Lab experiments using Microcystis aeruginosa as a model organism. While they briefly mention testing other species, the data are not shown in detail; however, they generalize the results to all HBFA.\nMinor Comments\n• Provide GAM model formulas and parameter settings in Supplementary Materials.\n• Include a graphical workflow linking datasets, analyses, and experiments.", "reviewer_2": "The manuscript presents an interesting and potentially valuable study on the effects of climate extremes on lake eutrophication. It does exhibit some new findings regarding the mechanisms, which might attract widespread attention. However, I feel that the overall quality and novelty might not match the high standards of NC journal. For the global remote-sensing data analysis section, some findings are not brand new, which weakened the novelty. Additionally, it could be more solid if the authors further clarify the data processing and experimental detail.\nSome specific comments are as follows:\n1. The writing can be further improved, including the standardization of language. Additionally, it is not necessary to use abbreviation for the word “phosphorus”, it is not such a complicated word. Using “Days 1 and 4…” instead of “days 1 and 4”; Using “Tube 1… Tube 2”, instead of “tube 1… tube 2”….\n2. The manuscript relies heavily on satellite-derived data and socioeconomic data analysis. However, there is an insufficient explanation on data integration and validation methods.\n3. The experimental design to assess the effects of temperature and pH stress on Microcystis aeruginosa is generally sound, but the manuscript lacks details regarding the number of biological replicates per treatment. Additionally, why the pH treatments are designed as and the pH was adjusted to 6, 7, 8, 9, and 10? The large range of pH value definitely lead to different data results, however, are these values indeed closely simulate the real nature condition in lakes?\n4. The study involves some in-situ sampling and filed experiments, it could be better to further clarify the details and environmental conditions during the experiment, for example, in which year? What are the water quality conditions and weather conditions?"}, "author_response": "Response to Reviewers\nReviewer #1\n1. The manuscript presents an interesting integration of global satellite data and laboratory\nexperiments to explain how climate extremes may intensify lake eutrophication through algal\nstress resistance mechanisms. The hypothesis is compelling, and the combination of large-scale\nand mechanistic approaches is commendable.\nResponse: Thank you for this positive and encouraging assessment. We appreciate the\nrecognition of the integrative framework that combines global satellite observations with\ncontrolled laboratory experiments to link large-scale eutrophication patterns with underlying\nmechanisms.\nIn response, we have further strengthened the connection between global-scale patterns and\ncellular and physiological processes in the revised manuscript. Specifically, we (i) extended\nthe analyses to natural mixed phytoplankton assemblages collected from Lake Taihu and Lake\nChaohu (Figure S3), (ii) expanded mechanistic experiments to four globally dominant harmful\nbloom-forming cyanobacterial genera (Microcystis, Raphidiopsis, Dolichospermum, and\nAphanizomenon) (Figures 2, 3, 4, 5, S2), and (iii) enhanced methodological transparency and\nuncertainty assessment throughout the study (Line 403–429). These additions reinforce both\nthe robustness of the experimental evidence and the generality of the proposed mechanisms.\nThe corresponding results are presented in Figures 2, 3, 4, 5, S2, and S3, and detailed\ndescriptions of the additional experimental design, field sampling, and analytical procedures\nhave been incorporated into the Methods section (Line 453–529). All modifications have been\nclearly indicated in the revised manuscript using blue text.\n[... figures and methods text omitted for brevity ...]\n2. I suggest including a more methodological description and adding an uncertainty evaluation (data preprocessing, model diagnostics, cross-validation).\nLimitations should be discussed more: The satellite dataset (GlobColour, 2002–2024) has uncertainties, but the manuscript does not mention atmospheric correction errors, cloud contamination, or sensor variability across sensors.\nResponse: Thank you for this important suggestion regarding methodological transparency,\nuncertainty evaluation, and the discussion of limitations associated with satellite-based datasets.\nIn response, we substantially revised the satellite data processing workflow, uncertainty\nassessment, and limitations discussion in the manuscript. To directly address concerns related\nto atmospheric correction errors, cloud contamination, and cross-sensor variability in the\noriginal GlobColour chlorophyll-a product, we replaced the GlobColour dataset (2002–2024)\nwith a recently published global lake algal bloom dataset developed by Wang et al. (DOI:\n10.1093/nsr/nwaf011). This dataset is specifically designed for inland freshwater systems and\nrelies exclusively on the internally consistent MODIS–Aqua archive (2003–2022), thereby\navoiding cross-sensor inconsistencies (Line 403–429).\n[... methods text omitted ...]\n3. From what I understand, the authors conducted Lab experiments using Microcystis\naeruginosa as a model organism. While they briefly mention testing other species, the data are\nnot shown in detail; however, they generalize the results to all HBFA.\nResponse: We agree that, in the original version of the manuscript, the evidence supporting\ngeneralization across harmful bloom-forming algae (HBFA) was not presented with sufficient\nclarity. In the revised manuscript, we have substantially expanded both the experimental\nevidence and its presentation to more rigorously justify this generalization, as outlined below.\n(i) Inclusion of multiple HBFA genera in laboratory experiments\nIn addition to Microcystis aeruginosa, we conducted parallel laboratory experiments on three\nother globally dominant bloom-forming cyanobacterial genera: Raphidiopsis,\nDolichospermum, and Aphanizomenon. These taxa differ markedly in morphology, phylogeny,\nand physiological traits. Despite these differences, all four genera consistently exhibited\nconvergent responses to heat stress, including elevated oxidative stress, enhanced\npolyphosphate accumulation, stabilisome formation, increased stress tolerance, and\nstrengthened vertical migration capacity (Figure 2, 3, 4, 5, S2).\n[... figures and methods text omitted ...]\n4. Provide GAM model formulas and parameter settings in Supplementary Materials.\nResponse: Thank you for this important suggestion. We fully agree that explicit reporting of\nGAM formulations and parameter settings is essential for transparency, reproducibility, and\nproper interpretation of the results. In response, we have substantially expanded the description\nof all GAM analyses and added a dedicated subsection entitled “Supplementary Methods for\nGAMs” to the Supplementary Information. Across the manuscript, GAMs were implemented\nusing a consistent methodological framework, but distinct GAM specifications were\nconstructed for different analytical purposes, corresponding to Figure 1, Figure S1, and Figure\nS4.\n[... methods text omitted ...]\nComment 5: Include a graphical workflow linking datasets, analyses, and experiments.\nResponse: Thank you for this constructive suggestion. In response, we have added a new\nschematic workflow figure to clarify the logical integration of datasets, analytical approaches,\nand experimental validation used in this study (Figure S5).\n[... methods text omitted ...]\nReviewer #2\n1. The manuscript presents an interesting and potentially valuable study on the effects of\nclimate extremes on lake eutrophication. It does exhibit some new findings regarding the\nmechanisms, which might attract widespread attention. However, I feel that the overall quality\nand novelty might not match the high standards of NC journal. For the global remote-sensing\ndata analysis section, some findings are not brand new, which weakened the novelty.\nAdditionally, it could be more solid if the authors further clarify the data processing and\nexperimental detail.\nResponse：We sincerely thank you for the careful and candid evaluation, and we appreciate\nthe concern regarding whether the novelty and overall contribution meet the standards of\nNature Communications, particularly with respect to the global remote-sensing analyses.\nWe agree that some broad associations between climate variability, anthropogenic pressure,\nand lake eutrophication have been reported previously. However, the central novelty of this\nstudy does not lie in re-identifying such correlations. Instead, our contribution is conceptual\nand mechanistic: we demonstrate how short-lived climatic extremes are converted into\nsustained bloom persistence through stress-memory-like physiological mechanisms in HBFA.\n[... methods text omitted ...]\n2. The writing can be further improved, including the standardization of language. Additionally, it is not necessary to use abbreviation for the word “phosphorus”, it is not such a complicated word. Using “Days 1 and 4…” instead of “days 1 and 4”; Using “Tube 1… Tube 2”, instead of “tube 1… tube 2”….\nResponse：Thank you for this helpful suggestion. We agree that improvements in language\nconsistency and stylistic standardization were necessary. Accordingly, we carefully revised the\nentire manuscript to ensure uniform terminology, capitalization, and formatting, and to improve\noverall clarity and readability.\n[... methods text omitted ...]\n3. The manuscript relies heavily on satellite-derived data and socioeconomic data analysis.\nHowever, there is an insufficient explanation on data integration and validation methods.\nResponse：Thank you for raising this important point. We agree that, given the strong reliance\non satellite-derived and socioeconomic datasets, a clearer explanation of data integration and\nvalidation is essential. We have therefore revised the manuscript to more explicitly describe\nhow these data streams are validated, integrated, and interpreted within the overall study\nframework (Line 403–429).\n[... methods text omitted ...]\n4. The experimental design to assess the effects of temperature and pH stress on Microcystis\naeruginosa is generally sound, but the manuscript lacks details regarding the number of\nbiological replicates per treatment.\nResponse: Thank you for this comment. We agree that explicitly reporting the number of\nbiological replicates is essential for assessing the robustness and reproducibility of the\nexperimental results. We have therefore revised the manuscript to clearly specify the number\nof biological replicates for all temperature and pH treatments (Line 491–492, 500–503, 552–\n554, 769–770, 783).\n[... methods text omitted ...]\n5. Additionally, why the pH treatments are designed as and the pH was adjusted to 6, 7, 8, 9,\nand 10? The large range of pH value definitely lead to different data results, however, are these\nvalues indeed closely simulate the real nature condition in lakes?\nResponse: Thank you for raising this important point. We agree that the ecological relevance\nof the selected pH treatments requires clear justification. We first clarify that pH 6 is no longer\nincluded in the revised manuscript; all experimental pH treatments are restricted to pH 7, 8, 9,\nand 10. References to pH 6 in the original version have been removed to avoid confusion (see\nrevised Methods, Line 463–464).\n[... methods text omitted ...]\n6. The study involves some in-situ sampling and filed experiments, it could be better to further\nclarify the details and environmental conditions during the experiment, for example, in which\nyear? What are the water quality conditions and weather conditions?\nResponse: Thank you for this constructive suggestion. We agree that clearly reporting the\ntiming and environmental context of the in situ sampling and field experiments is essential for\ntransparency, interpretability, and reproducibility. In response, we have revised the manuscript\nto provide detailed information on the date, meteorological conditions, and water-quality\ncharacteristics associated with the field investigations.\n[... methods text omitted ...]"}, {"round": 2, "reviewer_comments": {"reviewer_2": "I have read the revised version of this mauscript, I feel the quality has been remarkably improved. I will support its publication in NC."}, "author_response": ""}]

Decision letter