Spatiotemporal decoupling of littoral and lacustrine geosmin dynamics: Implications for early warning in drinking water reservoirs
# These authors contributed equally to this work.
a College of Environmental Science and Engineering, Ocean University of China, Qingdao 266100, China.
b Key Laboratory of Environmental Aquatic Chemistry, State Key Laboratory of Regional Environment and Sustainability, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China.
c Tianjin Hydraulic Research Institute, Tianjin 300061, China.
d Tianjin Yuqiao Reservoir Authority, Tianjin 301900, China.
e Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430070, China.
f University of Chinese Academy of Sciences, Beijing 100049, China.
* Corresponding to: Ming Su (mingsu@rcees.ac.cn)
Abstract
The relationship between cyanobacterial niche characteristics and the transport dynamics of harmful metabolites to drinking water intakes remains poorly understood. This study integrated a national survey with a five-year high-frequency monitoring program to characterize these dynamics, focusing on the potent odorant geosmin. The national investigation revealed that 14% of surveyed sites exceeded the odor threshold of 10 ng L-1, indicating a non-negligible risk. In the YQ Reservoir, Planktothrix agardhii was identified as a primary producer. Monitoring revealed a distinct spatiotemporal decoupling: shallow littoral zones functioned as production centers where P. agardhii biomass peaked 8 days prior to the lacustrine intake. Time-lagged correlation analysis indicated that littoral biomass predicts intake geosmin concentrations with a 5-week lead time (R2 = 0.41). Ammonium was identified as the key regulatory factor, exhibiting its strongest correlation with geosmin in littoral zones (R2 = 0.37), though this linkage attenuated during transport. This proposed mechanistic transport model and tiered framework shift surveillance from reactive intake sampling to proactive littoral sentinel stations, establishing a critical predictive window for preventive intervention in reservoir-dependent water supplies.
Keywords: Geosmin, taste and odor, Early warning, Planktothrix agardhii, Spatiotemporal decoupling
Introduction
The increasing frequency and intensity of harmful cyanobacterial blooms (CyanoHABs) in surface waters, driven by climate change and intensified anthropogenic activities, represent a significant global threat to water security (Feng et al., 2021; Heisler et al., 2008; Huisman et al., 2018; O’neil et al., 2012; Paerl et al., 2001; Paerl and Paul, 2012; Woolway et al., 2020). While CyanoHABs in eutrophic shallow lakes have been extensively characterized (Kosten et al., 2011; Paerl et al., 2011; Reinl et al., 2022; Rousso et al., 2020; Xu et al., 2009; Zhang et al., 2023), research in reservoir systems remains comparatively limited. In these critical drinking water resources, CyanoHABs directly jeopardize water safety through the production of toxins (Carmichael et al., 2001; Merel et al., 2013; Plaas and Paerl, 2020) and taste-and-odor (T&O) compounds (Burlingame et al., 2017; Lin et al., 2018; Su et al., 2015; Yang et al., 2008). Approximately half of China’s source-water reservoirs experience recurrent algal-derived T&O episodes (Sun et al., 2014). Unlike chemical contaminants, cyanobacterial populations exhibit exponential growth under favorable conditions (Asato, 2003; Mehnert et al., 2010), resulting in rapid metabolite release (Ai et al., 2025; Singh et al., 2005) that often exceeds the treatment capacity of conventional water utilities (Cerón-Vivas et al., 2022; Mustapha et al., 2021; Srinivasan and Sorial, 2011).
Managing CyanoHABs in reservoirs involves distinct challenges due to spatial and hydraulic complexity, which results in less predictable ecological regimes than those in natural lakes (Dordoni et al., 2023; Thornton, 1984). Conventional monitoring strategies, typically concentrated at lacustrine intakes near dams, often fail to capture the heterogeneous distribution of algal biomass and metabolites. Furthermore, the limited hydraulic residence time between the intake and the treatment plant provides an insufficient operational window for intervention. Littoral zones, with their shallow depth, generally warm earlier in spring and maintain higher light availability throughout the water column compared to deeper offshore areas (Bracchini et al., 2009). These nearshore regions also exhibit stronger sediment–water coupling, evident in diffusive nutrient fluxes and episodic wind-driven resuspension, which significantly enhance internal nutrient supply (Akbarzadeh et al., 2025; Qin et al., 2006). Together, these physicochemical factors (Cefalì et al., 2016; Zohary and Gasith, 2014) create an optimal ecological niche that supports the early proliferation of cyanobacteria. Despite their potential as sentinel sites, these areas are frequently under-monitored, and the mechanistic linkages between littoral zone dynamics and downstream water quality at intakes remain insufficiently quantified.
Geosmin and 2-methylisoborneol (MIB) are the primary terpenoid T&O compounds synthesized by cyanobacteria (Jüttner and Watson, 2007; Watson, 2003). Their low molecular weight and chemical stability allow them to penetrate conventional treatment barriers and persist in distribution networks (Lin et al., 2018; Zamyadi et al., 2015). While MIB episodes in China have been well studied and managed through hydraulic controls (Lu et al., 2023; Su et al., 2015), geosmin episodes are increasing in frequency across Northern Hemisphere reservoirs (Qiu et al., 2023; Wu et al., 2022). The high phylogenetic diversity of geosmin-producing taxa (AWWA, 2010) complicates the identification of biological sources and the development of targeted management strategies, emphasizing the need for research into the ecological niches and growth dynamics of these organisms in littoral zones.
This study analyzes recent national survey data to characterize the occurrence of geosmin across China’s drinking water sources and presents a five‑year, spatially resolved investigation of a representative reservoir. By integrating water chemistry, molecular profiling, and statistical modeling, this research examines the spatiotemporal decoupling between cyanobacterial communities in littoral zones and concentrations at the lacustrine intake. The objectives are to identify the primary cyanobacterial taxa responsible for geosmin production and to demonstrate how population dynamics within littoral hotspots determine downstream odor levels. This work provides an evidence‑based framework for shifting from reactive treatment toward proactive, risk‑based management of cyanobacterial T&O threats.
Materials and methods
Study area and sampling design
National-scale geosmin survey. This study compiles geosmin concentration data from drinking water sources across China, representing an extensive national assessment of this odorant. The dataset integrates long-term monitoring results with high-quality measurements from established literature, including a survey of 33 reservoirs in Guangdong province (Zhou et al., 2024). To improve data comparability, only quantitative geosmin data from raw water samples analyzed using GC–MS-based methods were included. For risk assessment, a uniform exceedance threshold of 10 ng L-1 was applied across all sites, and values reported below method detection limits (MDLs) were treated as non-exceedance records in the exceedance-rate analysis. In total, geosmin data from 214 distinct sampling sites were analyzed to characterize spatial heterogeneity and nationwide occurrence patterns.
Intensive monitoring at a reservoir in northern China. YQ Reservoir is a large, river-type reservoir located in northern China, with a total storage capacity of approximately 1.559 billion m3, and serves as an important drinking water source for a major metropolitan area in northern China (Xu et al., 2014; Zhai et al., 2017). In recent years, the reservoir has experienced recurrent seasonal taste-and-odor episodes (Cai et al., 2017; Qiu et al., 2021). Based on bathymetric and hydrodynamic gradients, seven sampling sites (YQ01–YQ07) were established along a littoral-to-lacustrine transect. Site YQ01, located in the deep lacustrine zone (maximum depth of approximately 12 m), was monitored using integrated water column sampling. Sites YQ02–YQ07 were situated in shallower littoral and transitional zones (water depth 2–4 m), where surface water samples were collected because the shallow depth and frequent mixing were unlikely to produce persistent vertical stratification (Holgerson et al., 2022). High-frequency sampling (2–3 times per week) was conducted during the annual odor season from 2019 to 2024. Samples for geosmin analysis, water chemistry, and phytoplankton assessment were collected concurrently, stored at 4 °C in the dark, and processed within 24 hours.
Analytical methods
Geosmin analysis. Geosmin concentrations were determined via headspace solid-phase microextraction (HS-SPME) coupled with gas chromatography–mass spectrometry (GC-MS; Agilent 6890GC-5973MSD) following established protocols (Su et al., 2015). The method detection limit (MDL) for geosmin was 1 ng L-1.
Physicochemical parameters. Water temperature, pH, dissolved oxygen (DO), and turbidity were recorded in-situ using a multiparameter sonde (EXO2, Xylem, USA). Concentrations of ammonium (NH4+-N), total nitrogen (TN), and total phosphorus (TP) were determined spectrophotometrically following Chinese national standard protocols. Specifically, TN and TP were quantified after alkaline potassium persulfate digestion (HJ 636-2012 and GB/T 11893-1989, respectively). NH4+-N was measured using the Nessler’s reagent method (HJ 535-2009).
Phytoplankton analysis. At each sampling station, a 1 L water sample was collected and immediately fixed with Lugol’s iodine solution (1% final concentration). After 48 h of static sedimentation, the samples were concentrated to a final volume of 50 mL (Hawkins et al., 2005). For phytoplankton identification and enumeration, a 0.1 mL aliquot of the concentrate was transferred to a Sedgewick–Rafter counting chamber. Phytoplankton were identified and counted under an upright microscope (Olympus CX23, Tokyo, Japan) at 400× magnification based on morphological characteristics according to (Hu and Wei, 2006). Multiple microscopic fields were examined for each sample. The phytoplankton cell density (\(C\), cells L\(^{-1}\)) was calculated using the following equation:
\[ C = \frac{n \times V_{con}}{V_{count} \times V_{sample}} \tag{1}\]
where \(n\) is the number of cells counted, \(V_{con}\) is the concentrated volume (50 mL), \(V_{count}\) is the volume of the counting chamber (0.1 mL), and \(V_{sample}\) is the initial sample volume (1 L).
Molecular characterization and bioinformatics
Molecular identification of geosmin producers. Biomass for molecular analysis was collected on 1.2-μm polycarbonate filters (Millipore) and stored at –80°C. Genomic DNA was extracted using the FastDNA™ Spin Kit for Soil (MP Biomedicals). The geoA gene, which encodes geosmin synthase in cyanobacteria, was amplified using functional primers 173AF/173AR (Giglio et al., 2008; Su et al., 2013). Taxonomic identification was performed via PCR amplification of the 16S rRNA gene V4 region (primers 515F/806R). Gene fusion was carried out to combine the (geoA) and 16S rRNA gene amplicons using the primers (173AF/806R). Amplicons were verified via agarose gel electrophoresis and subjected to high-throughput sequencing. Representative sequences were taxonomically assigned via BLAST searches, and a maximum-likelihood phylogenetic tree was constructed to confirm the identity of the dominant producer.
Bioinformatics pipeline. The sequencing was performed on the Illumina MiSeq PE300 platform (Illumina Inc., San Diego, USA) (Ravi et al., 2018). The raw paired-end sequences were first merged using FLASH (v1.2.11) and processed for quality filtering within the QIIME2 pipeline (v2022.2) (Bolyen et al., 2019). The UPARSE algorithm (v7.0.1090) was used to remove singletons and chimeras, and the sequences were clustered into Operational Taxonomic Units (OTUs) at a 97% similarity threshold (Edgar, 2013). For the 16S rRNA gene fragments, representative sequences were taxonomically assigned against the SILVA rRNA reference database (release 138) (Quast et al., 2012).
Data processing and statistical analysis
Temporal aggregation and lagged correlation. Data were aggregated to weekly means to align time series and reduce noise. To investigate time-lagged relationships, weekly geosmin concentrations were shifted forward by i weeks (i = 0–7) and paired with concurrent Planktothrix agardhii biomass. A log–log linear model was applied:
\[\log_{10}(G) = \alpha + \beta \log_{10}(P) + \varepsilon \tag{2}\]
where \(G\) is the geosmin concentration (ng L-1) and \(P\) is \(P. agardhii\) cell density (cells L-1). To minimize low-value interference, only paired observations with geosmin > 2 ng L-1 and positive biomass were included. Analysis was restricted to sites with ≥20 valid weekly observations. The coefficient of determination (\(R^2\)) was used to identify the optimal lag period (i.e., the early-warning window reported in the main text), and slope significance was evaluated via t-test.
Multivariate and spatial analysis. Principal component analysis (PCA) was performed on standardized variables to explore multivariate relationships. Zone-specific regressions (littoral, transitional, lacustrine) compared the strength of the ammonium–geosmin linkage. Bathymetric data were spatially interpolated via ordinary kriging to illustrate the reservoir’s depth gradient. All analyses were conducted in R (version 4.4.2) using the dplyr, lubridate, stats, and ggplot2 packages.
Results
National survey of geosmin occurrence in drinking water sources
We conducted a nationwide survey of 214 sites across 26 provinces to assess geosmin levels in China’s major drinking water sources, including reservoirs and key surface water bodies. Risk potential at each site was evaluated using peak recorded concentrations (Fig. S1). Values varied by more than four orders of magnitude, ranging from non-detectable levels to a maximum of 150 ng L-1 in East Coastal China. Across the study area, 14% of sites exceeded the 10 ng L-1 odor threshold (Fig. S2). The geospatial distribution showed marked spatial variability; elevated geosmin levels were primarily concentrated in eastern and southern China—specifically the Yangtze and Pearl River Deltas and parts of North China. In contrast, concentrations remained low or below detection limits throughout the northwest (e.g., Tibet, Xinjiang, and Qinghai) and certain central provinces, indicating a clear longitudinal gradient. These results identify geosmin as a pervasive odorant affecting drinking water quality across diverse climatic and geographic zones in China.
Seasonal and interannual dynamics in YQ Reservoir
High-frequency monitoring was conducted at the primary intake (site YQ01) of YQ Reservoir from 2019 to 2024 to characterize temporal geosmin variability (Fig. S3). The longitudinal record demonstrates a recurrent seasonal pattern (Fig. 2A, Fig. S4A). Concentrations typically remained below 5 ng L-1 during the spring (April–May), followed by an abrupt increase in June. Primary peaks occurred during the July–August period, with secondary elevations persisting through autumn and early winter (October–December) before a gradual decline. The 10 ng L-1 odor threshold was frequently exceeded during the study period. Monthly median values surpassed this threshold in 40% of summer months (July–September) and 100% of autumn/early winter months (October–December), indicating a recurrent odor risk during the warmer half of the year.
Interannual analysis (Fig. 2B) revealed a reduction in geosmin levels following the 2019–2020 monitoring period, which recorded a median concentration of 2.5 ng L-1. Between 2020 and 2024, annual medians stabilized within a range of 5.2 to 10 ng L-1. Statistical differences between years were confirmed by a Kruskal-Wallis test (χ2(4) = 64.3, p < 0.001). These data suggest a shift in baseline conditions or source inputs after the initial year, although seasonal periodicity remained a consistent feature.
Phytoplankton community dynamics and key cyanobacterial genera
Phytoplankton composition data from YQ Reservoir (site YQ01, 2019–2024) were analyzed to assess the relationship between community dynamics and geosmin occurrence (Fig. 2B). Total phytoplankton biomass (log10-transformed cell density) demonstrated seasonal variability, with minimum values in winter (January–March, ~107 cells L-1) followed by an increase starting in April. Biomass peaked during summer (July–September, >108 cells L-1) and remained elevated through autumn (October–December) before declining in late winter, a trajectory consistent with observed geosmin fluctuations.
Cyanobacteria dominated the assemblage throughout most of the monitoring period, with peak biomass recorded from June to October. Within this group, filamentous genera associated with geosmin production followed distinct successional patterns. Planktothrix agardhii and Aphanizomenon gracile were present year-round but reached maxima in summer and autumn. In contrast, Pseudanabaena catenata was primarily restricted to spring (April–May). Several Dolichospermum species (e.g., D. circinale, D. flos-aquae, D. planctonicum) and Aphanizomenon flos-aquae emerged in late spring, proliferated during summer, and declined by late autumn. Other taxa, including Trichormus variabilis, Microcoleus amoenus, and Pseudanabaena limnetica, occurred sporadically, with peaks confined to specific intervals between June and August.
The synchronization of elevated geosmin concentrations with high cyanobacterial biomass—particularly during blooms of Planktothrix, Dolichospermum and Aphanizomenon—points to a biological origin for the odor episodes. The temporal succession of these specific taxa likely drives the characteristic bimodal geosmin peaks observed during the summer and autumn months.
Identification and spatiotemporal distribution of the primary geosmin producer
Screening of 16S rRNA and geoA gene sequences identified the genera Planktothrix, Dolichospermum, and Aphanizomenon as the primary geosmin producers in YQ Reservoir (Fig. 2D-2G). Due to the high diversity and dynamic succession of filamentous cyanobacteria observed over the five-year monitoring period, taxonomic identification at the species level remains inconclusive based on current evidence. For subsequent analysis, P. agardhii was selected as the representative taxon for subsequent analysis because it showed a more persistent temporal pattern than other candidate producers across the monitoring period (Fig. 2B–C) and exhibited the strongest correlation with geosmin concentrations (Fig. S4B).
Spatial monitoring in 2023 demonstrated a clear biomass gradient; littoral sites (YQ03, YQ04, YQ06, and YQ07, Fig. S3) maintained significantly higher cell densities than the lacustrine site (YQ01). Maximum biomass reached 5.6 × 107 cells L-1 at littoral locations, nearly double the peak of 3.9 × 107 cells L-1 recorded at the lacustrine abstraction point. Biomass exceeded the 2 × 107 cells L-1 threshold in 8.7% of littoral samples, compared to 8.6% at YQ01. Unlike the ephemeral peaks of other cyanobacteria, P. agardhii persisted throughout the warm season (Fig. 2C).
The progression of bloom events (≥2 × 107 cells L-1) followed a spatial sequence from the littoral zone toward the dam (Fig. 2H). For example, YQ05 reached this threshold 8 days before YQ01. This temporal lag suggests that littoral zones function as primary proliferation centers, with subsequent advection driving the geosmin peaks observed at the lacustrine intake.
Distinct spatiotemporal patterns of two geosmin episodes
Two distinct odor episodes were identified in YQ Reservoir during 2023 (April–August and August–December). Spatial mapping of monthly average geosmin concentrations at seven representative sites revealed contrasting trajectories between littoral and lacustrine regions (Fig. 3).
During the first episode (April–August), littoral sites (YQ02, YQ03, YQ04, YQ06, and YQ07) recorded elevated concentrations in spring (April–June), ranging from 0 to 21.2 ng L-1. These values declined by August to between 0.6 and 5.8 ng L-1. In contrast, the lacustrine site (YQ01) and the transitional site (YQ05) exhibited a delayed response; concentrations increased progressively from April minima to reach a peak in August at 8.5 ng L-1.
The second episode (August–December) showed a reversed spatiotemporal gradient. Littoral sites reached peak concentrations in October, ranging from 0 to 36.4 ng L-1, before decreasing by December. Conversely, the lacustrine site (YQ01) exhibited a monotonic increase throughout the period, reaching 24.2 ng L-1 in December. Site YQ05 demonstrated an intermediate transition, with concentrations declining in November.
These divergent trajectories (Fig. 3) indicate a systematic spatiotemporal decoupling: littoral zones demonstrate earlier onset and decline during both episodes, while the lacustrine region shows delayed and prolonged elevations. This pattern reflects the integrated effects of initial production in littoral hotspots and subsequent advective transport toward the central reservoir.
Predictive relationship between Planktothrix agardhii dynamics and geosmin episodes
Time-lagged correlation analyses were performed between Planktothrix agardhii biomass and geosmin concentrations at the near-littoral site YQ05 and the lacustrine site YQ01 (Fig. 4A, Fig. S8). The analysis demonstrated divergent temporal coupling between the two locations. At YQ05, the strongest correlation (R2 = 0.5, p < 0.001) occurred at zero lag, indicating that geosmin levels fluctuated synchronously with local P. agardhii density. In contrast, the correlation at the lacustrine site YQ01 was weakest at zero lag but increased with time, peaking at a 5-week lag (R2 = 0.41, p < 0.001). This offset indicates that P. agardhii biomass functions as a concurrent indicator in the littoral production zone but acts as a leading indicator for geosmin risk at the distal lacustrine intake.
Regression models were established using the optimal lag periods for each site (Fig. 4B). For YQ05 (zero lag), the relationship followed the equation: log10(geosmin) = 0.244 + 0.288 × log10(P. agardhii biomass). For YQ01 (5-week lag), the model was: log10(geosmin) = 0.558 + 0.196 × log10(P. agardhii biomass). The higher slope observed at YQ01 suggests a potential amplification of geosmin concentrations during transport from littoral zones to the reservoir center.
Bi-weekly time series further illustrate this temporal phasing (Fig. 4C). At YQ05, P. agardhii and geosmin peaks were closely aligned, whereas at YQ01, biomass peaks preceded major odor episodes by approximately 6.7 weeks. This consistent lag serves as a predictive window; when P. agardhii biomass at YQ05 exceeded 5.0 × 107 cells L-1, geosmin at the lacustrine intake surpassed the 10 ng L-1 threshold within 4–6 weeks in 61.9% of observations. These results support the use of near-littoral monitoring for forecasting odor risks at central abstraction points.
Environmental drivers of geosmin and proposed transport pathway
Clear spatial gradient in nutrient loading and algal response was observed in YQ Reservoir ( Fig. S9). The littoral sites exhibit the highest vulnerability to organic enrichment, particularly during the summer. In this period, the littoral zone reaches a peak mean Chl-a of 0.072 mg L-1, which is significantly higher than the 0.061 mg L-1 observed in the lacustrine zone. This spatial disparity is mirrored in the Chemical Oxygen Demand (COD) and turbidity levels, with littoral sites showing higher COD (5.87 mg L-1) and lower Secchi Disk transparency (50.6 cm) than the deeper water sites. In contrast, the lacustrine zone displays higher stability regarding immediate organic surges but acts as a major sink for nitrogen. During winter, the Total Nitrogen (TN) in the lacustrine zone reaches a maximum of 3.56 mg L-1, a level higher than those found in the littoral or transitional areas during the same season. The transitional sites serve as a buffer, with water chemistry parameters like pH and TP typically falling between the littoral peaks and the lacustrine lows.
Principal component analysis (PCA) identified ammonium (NH4+-N) as the variable most strongly associated with geosmin variability in the PC1–PC2 plane (Fig. 5A). The first two principal components accounted for 60.6% of total variance. PC1 (43.5%) represented a thermal-nutrient productivity axis, characterized by high loadings for temperature (0.342), total phosphorus (TP, 0.39), and chlorophyll-a (Chl-a, 0.428). PC2 (17.1%) represented a nitrogen availability axis, dominated by ammonium (-0.449) and total nitrogen (TN, 0.432). Geosmin exhibited moderate positive loadings on both PC1 (0.122) and PC2 (0.229), suggesting that variability is driven by a combination of thermal-nutrient status and nitrogen availability.
Log-log regressions between ammonium and geosmin were performed across three reservoir zones: littoral (YQ02–YQ04, YQ06–YQ07), transitional (YQ05), and lacustrine (YQ01) (Fig. 5B). A clear spatial gradient in correlation strength was observed, with the highest R2 at littoral sites (0.365), followed by YQ05 (0.203), and the lacustrine site YQ01 (0.035). This attenuation suggests that the ammonium-geosmin linkage is most direct in production zones and undergoes progressive decoupling during transport toward the lacustrine intake.
Based on these results and the observed spatiotemporal dynamics of Planktothrix agardhii, a transport pathway for odor episodes is proposed (Fig. 5C). Shallow littoral areas (depths < 4–5 m) serve as primary sites for ammonium accumulation, P. agardhii proliferation, and geosmin production. Advective transport then moves cells and dissolved geosmin from these littoral zones through the transitional region (YQ05) to the lacustrine intake (YQ01). This framework integrates the spatial gradient in nutrient-geosmin coupling, the observed temporal lags between zones, and the bathymetric influence on water residence times. These elements provide a mechanistic basis for the origin and delayed manifestation of geosmin episodes at the reservoir’s primary abstraction point.
Discussion
Geosmin prevalence in Chinese drinking water sources
National survey data indicate that geosmin contamination is widespread across China’s primary drinking water sources, with 14% of sites exceeding the 10 ng L-1 odor threshold ((AWWA, 2010)). These results identify geosmin as a pervasive water quality challenge, consistent with observations in other climatological regions (Devi et al., 2021; Lin et al., 2018). The observed spatial heterogeneity—specifically the higher concentrations within the Yangtze and Pearl River Deltas—suggests that regional climatic conditions, watershed land use, and nutrient loading are primary determinants of geosmin risk. Given its low odor threshold and resistance to conventional treatment processes (Cook et al., 2001; Mustapha et al., 2021; Zamyadi et al., 2015), geosmin poses significant operational challenges for water utilities (Srinivasan and Sorial, 2011) and impacts compliance with water safety objectives such as Sustainable Development Goal 6 (Arora and Mishra, 2022). Consequently, these risks drive substantial investment in advanced treatment technologies and proactive source management (Yuan and Hofmann, 2022).
Littoral zones as production centers and spatiotemporal decoupling
While actinomycetes are established sources of geosmin in freshwater sediments (Gerber, 1979; Gerber and Lechevalier, 1965; Guo et al., 2024; Jensen et al., 1994), the synchronous peaks of Planktothrix agardhii biomass and geosmin at littoral sites (R2 = 0.5), combined with the molecular identification of the geoA gene (Fig. 2D, F), confirm a cyanobacterial origin in YQ Reservoir. This supports the consensus that cyanobacteria are the primary biotic source of odorants in well-oxygenated reservoir water columns (Chislock et al., 2021; Hooper et al., 2023; Izaguirre, 1992; Pham et al., 2020; Su et al., 2013; Suurnäkki et al., 2015). Detection of geoA-harboring cyanobacteria (Giglio et al., 2011; Giglio et al., 2008; Kutovaya and Watson, 2014) in this oxygenated environment indicates that geosmin is primarily derived from cyanobacterial activity, aligning with studies linking specific blooms to odor episodes globally (Churro et al., 2020; Hayashi et al., 2019; John et al., 2018).
Monitoring data identify littoral zones as critical ecological niches for P. agardhii proliferation. The filamentous morphology and low-light adaptation of P. agardhii facilitate its dominance in shallow, intermittently mixed littoral environments (Cao et al., 2025; Jia et al., 2019; Komárek and Johansen, 2015; Su et al., 2019). These traits allow for optimal light harvesting under the variable regimes typical of nearshore areas (Scott and Marcarelli, 2012; Zohary and Gasith, 2014). The high correlation between ammonium and geosmin specifically in littoral zones (R2 = 0.37) points to nitrogen availability from shoreline inputs as a key regulator. The positive correlation observed in situ in the reservoir is consistent with physiological evidence from laboratory studies on geosmin-producing cyanobacteria, which have shown that elevated ammonium concentrations can stimulate cyanobacterial growth and thereby increase overall geosmin production capacity (Canizales et al., 2021; Flores and Herrero, 2005; Saadoun et al., 2001; Yang et al., 2023). Thus, these findings collectively suggest that ammonium availability may serve as an indicator of geosmin production potential. Together, these results support the view that littoral zones act as the primary production centers for geosmin in YQ Reservoir.
The subsequent weakening of ammonium–geosmin correlation at the lacustrine intake (R2 = 0.03) reflects the decoupling of biological production from observed concentrations during transport. This spatial decoupling is quantified by the 8-day lag in bloom development between littoral sites and the dam, representing a period of physical advection and continued metabolic accumulation. During this process, the observed delay at the lacustrine intake may partly reflect intracellular geosmin release associated with cell lysis during downstream transport. As P. agardhii cells are transported toward the dam, where conditions become less favorable for growth, some cells may undergo lysis and release intracellular geosmin into the water column (Dubourg et al., 2015; Jia et al., 2019). In addition, geosmin concentrations at the lacustrine intake may be further affected by multiple physical and biogeochemical attenuation processes, including dilution through mixing with geosmin-poor water, microbial degradation (Clercin et al., 2021; Hammond et al., 2021; Hoefel et al., 2009), volatilization across the air–water interface (Soyluoglu et al., 2022), and photochemical transformation (e.g., UV degradation) (Kim et al., 2016; Mustapha et al., 2021).
Implications for monitoring and management
The conceptual mechanism depicted in Fig. 6 illustrates that maximum geosmin production is localized in littoral zones before advective transport through the limnetic zone. This sequence demonstrates a limitation in intake-centric monitoring: by the time geosmin reaches the lacustrine zone, the window for preventive intervention is minimal (Mohanty et al., 2024; Zamyadi et al., 2015). The identified 4–6 week lag between littoral proliferation and lacustrine contamination establishes a quantifiable window for risk mitigation. This estimate is supported by both the 5-week optimal lag identified from weekly lagged-correlation analysis and the 6.7-week delay inferred from the observed time-series pattern, which together indicate a comparable multi-week early-warning interval despite differences in method. By deploying littoral sentinel sites, utilities can move beyond reactive crisis management—such as emergency high-dose activated carbon application—toward preventive catchment-based strategies. This lead time is further supported by the lagged correlation at YQ01 (R2 = 0.41 at 5 weeks). This interval allows for optimized treatment costs and proactive resource allocation to intercept risks during the incipient growth phase, aligning with modern water safety frameworks (Recknagel et al., 2017).
This monitoring framework is most applicable to reservoirs where odor production is dominated by littoral hotspots (Izydorczyk et al., 2008). Because the warning lag may vary among systems, it should be calibrated according to reservoir-specific characteristics, including hydrodynamic conditions, shoreline-to-intake distance, water residence time, and the spatial distribution of odor-producing cyanobacteria(Chong et al., 2018; Shin et al., 2022; Song et al., 2023). Systems dominated by deep-water or benthic producers may require modified sentinel-station strategies (Gaget et al., 2020; Izaguirre and Taylor, 2007; Su et al., 2017). In addition, the biomass threshold of \(5.0 \times 10^7\) cells L\(^{-1}\) at YQ05 should be regarded as a preliminary indicator rather than a fixed decision threshold (Almuhtaram et al., 2021; Recknagel et al., 2017). Accordingly, it is better used as a practical reference for sentinel monitoring and early intervention. Exceedance at littoral stations may serve as an early-warning signal to intensify surveillance and, where feasible, support operational measures such as selective withdrawal or hydrodynamic control to limit the transport of geosmin-rich water toward critical intake depths (Lehman et al., 2009; Song et al., 2023, 2023).
While the current dataset captures robust seasonal and interannual trends, longer-term observations are required to evaluate these dynamics under shifting climate scenarios. Furthermore, while P. agardhii is the dominant producer in this system, other reservoirs may be driven by different taxa, such as Dolichospermum or Aphanizomenon (Choo et al., 2025; Izaguirre and Taylor, 2004; Su et al., 2013).
Future research should focus on: (1) integrating molecular tools to characterize the diversity and expression of geosmin synthesis genes, (2) coupling ecological monitoring with high-resolution hydrodynamic modeling to refine transport time estimates, and (3) conducting cost-benefit analyses to evaluate the economic feasibility of expanded littoral monitoring. By establishing quantifiable early-warning indicators based on the littoral-to-lacustrine transport mechanism, this work provides a technical foundation for improving drinking water security.
Conclusions
This study demonstrates that geosmin risk in drinking water reservoirs is governed by a systematic spatiotemporal decoupling, whereby shallow littoral zones function as primary production hotspots, with odor manifestation at the lacustrine intake lagging by several weeks. This quantifiable lag transforms the monitoring paradigm, establishing that proactive surveillance of littoral cyanobacteria provides a critical early-warning window, thereby enabling a shift from reactive treatment toward predictive, source-water management for ensuring drinking water security.
Acknowledgements
This work was financially supported by the National Natural Science Foundation of China (W2412156, 52030002), and Science and Technology Plan Project of Tianjin Water Bureau (JSGJ202031).