Models are simplified descriptions of reality and are intrinsically limited by the assumptions that have been introduced in their formulation. With the development of automatic calibration toolboxes, finding optimal parameters that suit the environmental system has become more convenient. Here, we explore how optimization toolboxes can be applied innovatively to uncover flaws in the physical formulations of models. We illustrate this approach by evaluating the effect of simplifications embedded in the formulation of a widely used hydro-thermodynamic model. We calibrate a Delft3D model based on temperature profiles for a case study, Lake Morat (Switzerland), through the DYNO-PODS optimization tool. The results show that higher values of the light extinction coefficient can compensate for neglecting the fraction β of short-wave radiation absorbed at the surface of the water. This leads to unrealistic values of the light extinction coefficient, as the optimization pushes its value toward the limit of no transparency, consistent with the need to reproduce a significant absorption at the surface. Although it is well known that β is significantly larger than zero, its absence from the model was never noticed as critical. Automatic calibration tools provide valuable diagnostic insights into the physical robustness of models, enabling more precise evaluation of their structural integrity and performance.

The sloping boundaries of stratified aquatic systems, such as lakes, are crucial environmental dynamic zones. While the role of sloping boundaries as energy dissipation hotspots is well established, their contribution to triggering large-scale motions has received less attention. This review delves into the development of thermally driven cross-shore flows on sloping boundaries under weak wind conditions. We specifically examine 'thermal siphons' (TS), a dynamical process that occurs when local free convection transforms into a horizontal circulation over sloping boundaries. Thermal siphons result from bathymetrically induced temperature (i.e. density) gradients when a lake experiences a uniform surface buoyancy flux, also known as differential cooling or heating. In the most common case of differential cooling of waters above the temperature of maximum density, TS lead to an overturning circulation characterised by a downslope density current and a surface return flow within a convective environment. Field observations, laboratory experiments and high-fidelity simulations of TS provide insights into their temporal occurrence, formation mechanisms, water transport dynamics and cross-shore pathways, addressing pivotal questions from an aquatic system perspective. Fluid mechanics is a fundamental tool in addressing such environmental questions and thereby serves as the central theme in this review.

Turbidity currents, gravity-driven sediment-laden flows, govern material transport and shape underwater landscapes across diverse environments. While extensively studied in marine settings, their dynamics in freshwater systems remain underexplored. We present multi-temporal, multi-scale observations from the Aare Delta of Lake Brienz (Switzerland), combining Acoustic Doppler Current Profiler measurements with high-resolution repeated bathymetric surveys to capture turbidity currents and resulting morphological changes. A single turbidity current reached near-bed velocities of 0.65 m/s and induced upslope migration of a cyclic step by at least 4 m within minutes. Characterisation of overall bedform migration across the channel network over five years highlights the cumulative imprint of turbidity currents and frames these short-term changes within the broader context of channel reworking. Together, our results show striking similarities in overall flow structure, bedform morphologies, and dynamics between freshwater and marine systems, underscoring the value of lacustrine settings as accessible, scaled-down environments for advancing turbidity current research.

Lake Cadagno is a meromictic alpine lake characterized by permanent stratification, which creates a permanent anoxic environment that supports the growth of anoxygenic phototrophic sulfur bacteria. The seasonality and interseasonality of these microorganisms were examined over a three-year period (2019–2021) through regular monitoring of the water column. A variety of physical–chemical parameters, including temperature, conductivity, light, oxygen and sulfide concentrations, and the community composition of anoxygenic phototrophic sulfur bacteria in the chemocline were recorded, to observe potential influence of external weather conditions. Our findings indicate that, despite the lake’s consistent physical and chemical stratification, the composition of the anoxygenic phototrophic sulfur bacteria community exhibited notable variations in response to external environmental factors, including changes in rainfall and light irradiance. Specifically, we observed different growth dynamics in the purple (PSB) and green (GSB) sulfur bacteria communities over the three years of monitoring. These variations underscore the complexity of biogeochemical cycles in meromictic lakes and the impact of external environmental factors on this ancestral microbial community dynamics. The results provide valuable insights into the stability of redox-stratified environments, offering a modern analog for ancient aquatic ecosystems. This research emphasizes the importance of long-term regular monitoring to capture interannual dynamics and assess the implications of climate change on such unique ecosystems.

Ocean Bottom Seismometers (OBS) are primarily designed to record seismic signals in marine environments. However, the fundamental principles of their operation as broadband instruments enable a wider range of applications. Here, we demonstrate that OBS systems can also effectively monitor the impact of sudden strong wind events on lakes. We conducted an experiment in Lake Lucerne (Switzerland), at the Muota Delta, using five OBS equipped with hydrophones. In addition, an Acoustic Doppler Current Profiler (ADCP) was moored to record current speed and direction. Comprehensive data analysis of weather measurements reveals a link between a strong wind and seismic signature. A foehn event that occurred on October 10th, 2023 in Switzerland caused a wind-lake interaction that triggered an internal current in the lake, with intensified deep near-bed velocity flowing in the upwind direction. Seismic signals induced by foehn winds correlate strongly with recorded increased sediment load, likely from resuspension or erosion, demonstrating how ocean-bottom seismometer networks can enhance environmental monitoring.

Swiss geosciences are at the heart of major national and global challenges that our society is facing, including climate change and weather extremes, long-term environmental trends, natural hazards (e.g. landslides, floods, seismic and volcanic events), the energy transition towards a green and blue economy, sustainable management of natural resources (e.g. water, soils, ecosystems, critical metals, building materials), and the increasing interest towards extra-terrestrial life and resources. In order to remain at the forefront of geosciences research and provide the fertile environment necessary for groundbreaking findings addressing these societal challenges, the Swiss Geosciences Community proposes five large infrastructures, organised into four pillars built on a common datadriven research foundation. [...]

The pathway of dense river inflows into lakes, which affects the lake water quality, is not accurately predicted by existing models. The pathway of a dense riverine inflow in a lake with a submerged canyon is analyzed based on measurements during a 4-month period of weakening lake stratification and weakening density excess between river and epilimnion. In line with models, the dense riverine inflow plunges upon entering the lake, continues as an underflow on the sloping lake bottom, and finally intrudes at its level of neutral buoyancy. Underflow and interflow velocities are O(0.1 m s⁻¹). The river inflow is finally trapped in the pycnocline most of the time, even when the river's density excess and the lake's stratification are marginal. This trapping in the pycnocline is explained by the reduction of the inflow density excess due to the intense plunging mixing, which is an order of magnitude larger than that obtained in confined laboratory flumes. The pathway of the dense riverine inflow is affected by interactions of the underflow with the lake bottom and sedimentary processes. A canyon carved by the underflows confines and accelerates the underflow, which enhances its capacity to entrain and carry sediment. The entrainment of sediment that was previously deposited on the canyon bottom accelerates the underflow. Due to both effects, the underflow can temporarily break through the pycnocline and reach the hypolimnion. Existing models explain these observations qualitatively, but a quantitative prediction would require better parameterizations of the plunging mixing and the sedimentary processes.

The extent of littoral influence on lake gas dynamics remains debated in the aquatic science community due to the lack of direct quantification of lateral gas transport. The prevalent assumption of diffusive horizontal transport in gas budgets fails to explain anomalies observed in pelagic gas concentrations. Here, we demonstrate through high-frequency measurements in a eutrophic lake that daily convective horizontal circulation generates littoral-pelagic advective gas fluxes one order of magnitude larger than typical horizontal fluxes used in gas budgets. These lateral fluxes are sufficient to redistribute gases at the basin-scale and generate concentration anomalies reported in other lakes. Our observations also contrast the hypothesis of pure, nocturnal littoral-to-pelagic exchange by showing that convective circulation transports gases such as oxygen and methane toward both the pelagic and littoral zones during the daytime. This study challenges the traditional pelagic-centered models of aquatic systems by showing that convective circulation represents a fundamental lateral transport mechanism to be integrated into gas budgets.

L'état global de la qualité des lacs en Suisse a tendance à s'améliorer, mais de nombreux problèmes d'oxygénation dans le fond des lacs persistent. Ce phénomène, lié aux concentrations actuelles et passées de phosphore, demeure préoccupant. La modélisation au coeur de cette étude vise à guider le choix d'approches simples mais appropriées pour la gestion des lacs. Deux exemples d'applications pour les lacs de Schiffenen et de Baldegg sont illustrés.

Freshwater algae exhibit complex dynamics, particularly in meso-oligotrophic lakes with sudden and dramatic increases in algal biomass following long periods of low background concentration. While the fundamental prerequisites for algal blooms, namely light and nutrient availability, are well-known, their specific causation involves an intricate chain of conditions. Here we examine a recent massive Uroglena bloom in Lake Geneva (Switzerland/France). We show that a certain sequence of meteorological conditions triggered this specific algal bloom event: heavy rainfall promoting excessive organic matter and nutrients loading, followed by wind-induced coastal upwelling, and a prolonged period of warm, calm weather. The combination of satellite remote sensing, in-situ measurements, ad-hoc biogeochemical analyses, and three-dimensional modeling proved invaluable in unraveling the complex dynamics of algal blooms highlighting the substantial role of littoral-pelagic connectivities in large low-nutrient lakes. These findings underscore the advantages of state-of-the-art multidisciplinary approaches for an improved understanding of dynamic systems as a whole.

In large and stratified lakes, substantial methane stocks are often observed within the metalimnion. The origin of the methane (CH₄) accumulated in the metalimnion during stratification, which can sustain significant emissions during convective mixing, is still widely debated. While commonly attributed to the transport of methane produced anaerobically ex situ, recent evidence suggests that oxic in situ methane production could also contribute to metalimnetic methane peaks. Here, we assessed the origin, that is, pelagic CH₄ production or transport of sublittoral CH₄ through the interflow, of metalimnetic methane in Lake Geneva, the largest lake in Western Europe. Microbial diversity data do not support the hypothesis of oxic methane production in the metalimnion. In contrast, both spatial and temporal surveys of methane show that maxima occur at depths and sites most affected by the Rhône River inflow. Methane δ¹³C values point to an anaerobic sublittoral methane source, within a biogeochemical hotspot close to the river delta region, and an efficient transport across several kilometers in a vertically well-constrained metalimnion. Our current findings emphasize the indirect role of river interflows for the long-range transport of CH₄ produced in sediment biogeochemical hotspots, even for large lakes where sublittoral habitats represent a fairly limited fraction of the lake volume.

Methane (CH₄) accumulation in the well-oxygenated lake epilimnion enhances the diffusive atmospheric CH₄ emission. Both lateral transport and in situ oxic methane production (OMP) have been suggested as potential sources. While the latter has been recently supported by increasing evidence, quantifying the exact contribution of OMP to atmospheric emissions remains challenging. Based on a large high-resolution field data set collected during 2019-2020 in the deep stratified Lake Stechlin and on three-dimensional hydrodynamic modeling, we improved existing CH₄ budgets by resolving each component of the mass balance model at a seasonal scale and therefore better constrained the residual OMP. All terms in our model showed a large temporal variability at scales from intraday to seasonal, and the modeled OMP was most sensitive to the surface CH₄ flux estimates. Future efforts are needed to reduce the uncertainties in estimating OMP rates using the mass balance approach by increasing the frequency of atmospheric CH₄ flux measurements.

Lakes emit substantial amounts of carbon dioxide (CO₂) into the atmosphere, but why they do remains debated. The long-standing vision of lakes as solely respirators of the organic matter leaking from the soils has been challenged by evidence that inorganic carbon produced by weathering of the catchment bedrock could also support lake CO₂ emissions. How inorganic carbon inputs ultimately generate lake CO₂ outgassing remains a blind spot. We develop and introduce a calcite module in a coupled one-dimensional physical-biogeochemical model that we use to simulate the carbon cycle of the large Lake Geneva over the past 40 years. We mechanistically demonstrate how the so-far neglected process of calcite precipitation boosts net CO₂ emissions at the annual scale. Far from being anecdotal, we show that calcite precipitation could explain CO₂ outgassing across various lakes globally, including some of the largest lakes in the world.

Climate change is contributing to rapid changes in lake ice cover across the Northern Hemisphere, thereby impacting local communities and ecosystems. Using lake ice cover time-series spanning over 87 yr for 43 lakes across the Northern Hemisphere, we found that the interannual variability in ice duration, measured as standard deviation, significantly increased in only half of our studied lakes. We observed that the interannual variability in ice duration peaked when lakes were, on average, covered by ice for about 1 month, while both longer and shorter long-term mean ice cover duration resulted in lower interannual variability in ice duration. These results demonstrate that the ice cover duration can become so short that the interannual variability rapidly declines. The interannual variability in ice duration showed a strong dependency on global temperature anomalies and teleconnections, such as the North Atlantic Oscillation and El Niño–Southern Oscillation. We conclude that many lakes across the Northern Hemisphere will experience a decline in interannual ice cover variability and shift to open water during the winter under a continued global warming trend which will affect lake biological, cultural, and economic processes.

Three-dimensional lake hydrodynamic model is a powerful tool widely used to assess hydrological condition changes of lake. However, its computational cost becomes problematic when forecasting the state of large lakes or using high-resolution simulation in small-to-medium size lakes. One possible solution is to employ a data-driven emulator, such as a deep learning (DL) based emulator, to replace the original model for fast computing. However, existing DL-based emulators are often black-box and data-dependent models, causing poor interpretability and generalizability in practical applications. In this study, a data-driven emulator is established using deep neural network (DNN) to replace the original model for fast computing of three-dimensional lake hydrodynamics. Then, the Koopman operator and transfer learning (TL) are employed to enhance the interpretability and generalizability of the emulator. Finally, the generalizability of DL-based emulators is comprehensively analyzed through linear regression and correlation analysis. These methods are tested against an existing hydrodynamic model of Lake Zurich (Switzerland) whose data was provided by an open-source web-based platform called Meteolakes/Alplakes. According to the results, (1) The DLEDMD offers better interpretability than DNN because its Koopman operator reveals the linear structure behind the hydrodynamics; (2) The generalization of the DL-based emulators in three-dimensional lake hydrodynamics are influenced by the similarity between the training and testing data; (3) TL effectively improves the generalizability of the DL-based emulators.

Accumulation of methane in oxic waters of lakes and the ocean has been widely reported. Despite the importance for the greenhouse gas budget, mechanistic controls of such "methane paradox" remain elusive. Here, we use a combination of CH₄ concentration and isotopic (δ¹³C_CH4, δD_H2O and δ¹⁸O_H2O) measurements, plankton incubations and microbial community assessments to demonstrate the existence of the methane paradox in oxygenated waters of a meromictic lake (Lake Cadagno, Switzerland). Following mass dynamics using water isotopes, we exclude the possibility that the accumulation of CH₄ at the thermocline results solely from lateral transport. Interannual variability in the magnitude of the methane paradox (between 0.5 and 5 μmol L^-1) is associated to stratification patterns, changes in zooplankton biomass and planktonic detritus accumulation along density gradients, as well as fluctuating microbial cell numbers. The links between hydrodynamic conditions, aggregation of planktonic detritus and its microbiome, as well as the accumulation of CH₄ in the water column are further supported by high-resolution echosounder revealing backscatter maxima at the top of the thermocline, where detritus is effectively trapped, and by oxic incubations showing that CH₄ is produced in zooplankton detritus (0.046 nmol L^-1 to 0.095 CH₄ mg dry mass L^-1 d^-1). Our results also show that detritus-hosted methanogenesis is stimulated through the addition of methylphosphonate, suggesting that zooplankton-associated microbiomes exploit organic phosphorus compounds to release CH₄. Understanding the variability of the methane paradox in relation to changing hydrodynamics and plankton communities will be crucial to predict the future role of lakes in the global methane budget.

Long-term effects of climate change on lakes globally will include a substantial modification in the thermal regime and the oxygen solubility of lakes, resulting in the alteration of ecosystem processes, habitats, and concentrations of critical substances. Recent efforts have led to the development of long-term model projections of climate change effects on lake thermal regimes and oxygen solubility. However, such projections are hardly ever confronted with observations extending over multiple decades. Furthermore, global-scale forcing parameters in lake models present several limitations, such as the need of significant downscaling. In this study, the effects of climate change on thermal regime and oxygen solubility were analyzed in the four largest French peri-alpine lakes over 1850–2100. We tested several one-dimensional (1D) lake models' robustness for long-term variations based on up to 63 years of limnological data collected by the French Observatory of LAkes (OLA). Here, we evaluate the possibility of forcing mechanistic models by following the long-term evolution of shortwave radiation and air temperature while providing realistic seasonal trends for the other variables for which local-scale downscaling often lacks accuracy. Based on this approach, MyLake, forced by air temperatures and shortwave radiations, predicted accurately the variations in the lake thermal regime over the last 4 to 6 decades, with RMSE < 1.95 ^∘C. Over the previous 3 decades, water temperatures have increased by 0.46 ^∘C per decade (±0.02 ^∘C) in the epilimnion and 0.33 ^∘C per decade (±0.06 ^∘C) in the hypolimnion. Concomitantly and due to thermal change, O₂ solubility has decreased by −0.104 mg L⁻¹ per decade (±0.005 mg L⁻¹) and −0.096 mg L⁻¹ per decade (±0.011 mg L⁻¹) in the epilimnion and hypolimnion, respectively. Based on the shared socio-economic pathway SSP370 of the Intergovernmental Panel on Climate Change (IPCC), peri-alpine lakes could face an increase of 3.80 ^∘C (±0.20 ^∘C) in the next 70 years, accompanied by a decline of 1.0 mg L⁻¹ (±0.1 mg L⁻¹) of O₂ solubility. Together, these results highlight a critical alteration in lake thermal and oxygen conditions in the coming decades, and a need for a better integration of long-term lake observatories data and lake models to anticipate climate effects on lake thermal regimes and habitats.

Introduction: Bioconvection, a phenomenon characterized by the collective upward swimming of motile microorganisms, has mainly been investigated within controlled laboratory settings, leaving a knowledge gap regarding its ecological implications in natural aquatic environments. This study aims to address this question by investigating the influence of bioconvection on the eco-physiology of the anoxygenic phototrophic sulfur bacteria community of meromictic Lake Cadagno.
Methods: Here we comprehensively explore its effects by comparing the physicochemical profiles of the water column and the physiological traits of the main populations of the bacterial layer (BL). The search for eco-physiological effects of bioconvection involved a comparative analysis between two time points during the warm season, one featuring bioconvection (July) and the other without it (September).
Results: A prominent distinction in the physicochemical profiles of the water column centers on light availability, which is significantly higher in July. This minimum threshold of light intensity is essential for sustaining the physiological CO₂ fixation activity of Chromatium okenii, the microorganism responsible for bioconvection. Furthermore, the turbulence generated by bioconvection redistributes sulfides to the upper region of the BL and displaces other microorganisms from their optimal ecological niches.
Conclusion: The findings underscore the influence of bioconvection on the physiology of C. okenii and demonstrate its functional role in improving its metabolic advantage over coexisting phototrophic sulfur bacteria. However, additional research is necessary to confirm these results and to unravel the multiscale processes activated by C. okenii’s motility mechanisms.

Gravity currents contribute to the transport of heat and mass in atmospheric and aquatic environments. In aquatic systems subject to daily surface cooling, gravity currents propagate through turbulent convective surroundings. Yet, the effects of thermal convection on aquatic gravity currents remain to be quantified. This paper demonstrates how the interaction between penetrative convection and downslope gravity currents impacts the fluid dynamics and transport across littoral aquatic systems. We performed field experiments in a wind-sheltered lake experiencing differential cooling to resolve the dynamics of thermally driven gravity currents in convective environments. Our in situ observations reveal that convective plumes penetrate gravity currents, generating large vertical fluctuations that foster the erosion of the stratified layer. This enhanced vertical mixing destroys the stratified downslope flow and limits the basin-scale transport. Our results demonstrate that the interaction between penetrative convection and downslope gravity currents controls the littoral-pelagic connectivity in aquatic ecosystems.

In ice-covered lakes, near-bottom oxygen concentration decreases for most of the wintertime, sometimes down to the point that bottom waters become hypoxic. Studies insofar have reached divergent conclusions on whether climate change limits or reinforces the extent and duration of hypoxia under ice, raising the need for a comprehensive understanding of the drivers of the dissolved oxygen (DO) dynamics under lake ice. Using high-temporal resolution time series of DO concentration and temperature across 14 mountain lakes, we showed that the duration of bottom hypoxia under ice varies from 0 to 236 d within lakes and among years. The variability of hypoxia duration was primarily explained by changes in the decay rate of DO above the lake bottom rather than by differences in DO concentration at the ice onset or in the ice-cover duration. We observed that the DO decay rate was primarily linked to physical controls (i.e., deep-water warming) rather than biogeochemical drivers (i.e., proxies for lake or catchment productivity). Using a simple numerical model, we provided a proof-of-concept that the near-bed stratification can be the mechanism tying the DO decay rate to the sediment heat release under the ice. We ultimately showed that the DO decay rate and hypoxia duration are driven by the summer light climate, with faster oxygen decline found under the ice of clearer cryostratified alpine lakes. We derived a framework theorizing how the hypoxia duration might change under the ice of alpine lakes in a warmer climate.

In alkaline freshwater systems, the apparent absence of carbon limitation to gross primary production (GPP) at low CO₂ concentrations suggests that bicarbonates can support GPP. However, the contribution of bicarbonates to GPP has never been quantified in lakes along the seasons. To detect the origin of the inorganic carbon maintaining GPP, we analyze the daily stoichiometric ratios of CO₂–O₂ and alkalinity–O₂ in a deep hardwater lake. Results show that aquatic primary production withdraws bicarbonate from the alkalinity pool for two-thirds of the year. Alkalinity rather than CO₂ is the dominant inorganic carbon source for GPP throughout the stratified period in both the littoral and pelagic environments. This study sheds light on the neglected role of alkalinity in the freshwater carbon cycle throughout an annual cycle.

Recent efforts using microstructure turbulence measurements have contributed to our understanding of the overall energy budget in lakes and linkages to vertical fluxes. A paucity of lake-wide turbulence measurements hinders our ability to assess how representative such budgets are at the basin scale. Using an autonomous underwater glider equipped with a microstructure payload, we explored the spatial variability of turbulence in pelagic and near-shore regions of Lake Geneva. Dissipation rates of kinetic energy and thermal variance were estimated by fitting temperature gradient fluctuations spectra to the Batchelor spectrum. In deep waters, turbulent dissipation rates in the surface and thermocline were mild (∼10^-8 W kg^-1) and weakened towards the hypolimnion (∼10^-11 - 10^-10 W kg^-1). The seasonal thermocline exhibited inhibited interior mixing, with extremely low values of mixing efficiency (Ri_f ≪ 0.1). In contrast, in the slope zone, a band of significantly enhanced energy dissipation (∼5 × 10^-8 W kg^-1) extended well above the bottom boundary layer and was associated with strong, efficient mixing (Ri_f > 0.17). The resulting contribution of the slope region to basin-scale mixing was large, with 90% of the basin-wide mixing - and only 30% energy dissipation - occurring within 4 km of the shoreline. This boundary mixing will modify overturning circulation and the transport pathways of dissolved compounds exchanged with the sediments. The dynamics responsible for the shift in the mixing regime, which appears crucial for the mixing budget of lakes, could not be fully unraveled with the collected observations. Additional model data analyses hint at the role of submesoscale instabilities.

Severe deterioration of water quality in lakes, characterized by overabundance of algae and declining dissolved oxygen in the deep lake (DO_B), was one of the ecological crises of the 20th century. Even with large reductions in phosphorus loading, termed "reoligotrophication," DO_B and chlorophyll (CHL) have often not returned to their expected pre-20th-century levels. Concurrently, management of lake health has been confounded by possible consequences of climate change, particularly since the effects of climate are not neatly separable from the effects of eutrophication. Here, using Lake Geneva as an iconic example, we demonstrate a complementary alternative to parametric models for understanding and managing lake systems. This involves establishing an empirically-driven baseline that uses supervised machine learning to capture the changing interdependencies among biogeochemical variables and then combining the empirical model with a more conventional equation-based model of lake physics to predict DO_B over decadal time-scales. The hybrid model not only leads to substantially better forecasts, but also to a more actionable description of the emergent rates and processes (biogeochemical, ecological, etc.) that drive water quality. Notably, the hybrid model suggests that the impact of a moderate 3°C air temperature increase on water quality would be on the same order as the eutrophication of the previous century. The study provides a template and a practical path forward to cope with shifts in ecology to manage environmental systems for non-analogue futures.

When lakes experience surface cooling, the shallow littoral region cools faster than the deep pelagic waters. The lateral density gradient resulting from this differential cooling can trigger a cold downslope density current that intrudes at the base of the mixed layer during stratified conditions. This process is known as a thermal siphon (TS). TSs flush the littoral region and increase water exchange between nearshore and pelagic zones; thus, they may potentially impact the lake ecosystem. Past observations of TSs in lakes are limited to specific cooling events. Here, we focus on the seasonality of TS-induced lateral transport and investigate how seasonally varying forcing conditions control the occurrence and intensity of TSs. This research interprets 1-year-long TS observations from Rotsee (Switzerland), a small wind-sheltered temperate lake with an elongated shallow region. We demonstrate that TSs occur for more than 50 % of the days from late summer to winter and efficiently flush the littoral region within ∼10 h. We further quantify the occurrence, intensity, and timing of TSs over seasonal timescales. The conditions for TS formation become optimal in autumn when the duration of the cooling phase is longer than the time necessary to initiate a TS. The decrease in surface cooling by 1 order of magnitude from summer to winter reduces the lateral transport by a factor of 2. We interpret this transport seasonality with scaling relationships relating the daily averaged cross-shore velocity, unit-width discharge, and flushing timescale to the surface buoyancy flux, mixed-layer depth, and lake bathymetry. The timing and duration of diurnal flushing by TSs relate to daily heating and cooling phases. The longer cooling phase in autumn increases the flushing duration and delays the time of maximal flushing relative to the summer diurnal cycle. Given their scalability, the results reported here can be used to assess the relevance of TSs in other lakes and reservoirs.

Small-scale turbulent mixing plays a pivotal role in shaping ocean circulation and a broad range of physical and biogeochemical processes. Despite advances in our understanding of the geophysical processes responsible for this mixing, the nature and importance of biomixing - turbulent mixing caused by marine biota - are still debated. A major source of uncertainty pertains to the efficiency of biomixing (the fraction of the turbulent energy produced through swimming that is spent in mixing the ocean vertically), which the few in situ observations available suggest to be much lower than that of geophysical turbulence. Here we shed light on this problem by analysing 14 days of continuous measurements of centimetre-scale turbulence in an area of coastal upwelling. We show that turbulent dissipation is elevated 10- to 100-fold (reaching 10^-6-10^-5 W kg^-1) every night of the survey due to the swimming activity of large aggregations of anchovies that gather regularly over the spawning season. Turbulent mixing is invigorated concurrently with dissipation, and occurs with an efficiency comparable to that of geophysical turbulence. Our results demonstrate that biologically driven turbulence can be a highly effective mixing agent, and call for a re-examination of its impacts on productive upper ocean regions.

Empirical evidence demonstrates that lakes and reservoirs are warming across the globe. Consequently, there is an increased need to project future changes in lake thermal structure and resulting changes in lake biogeochemistry in order to plan for the likely impacts. Previous studies of the impacts of climate change on lakes have often relied on a single model forced with limited scenario-driven projections of future climate for a relatively small number of lakes. As a result, our understanding of the effects of climate change on lakes is fragmentary, based on scattered studies using different data sources and modelling protocols, and mainly focused on individual lakes or lake regions. This has precluded identification of the main impacts of climate change on lakes at global and regional scales and has likely contributed to the lack of lake water quality considerations in policy-relevant documents, such as the Assessment Reports of the Intergovernmental Panel on Climate Change (IPCC). Here, we describe a simulation protocol developed by the Lake Sector of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP) for simulating climate change impacts on lakes using an ensemble of lake models and climate change scenarios for ISIMIP phases 2 and 3. The protocol prescribes lake simulations driven by climate forcing from gridded observations and different Earth system models under various representative greenhouse gas concentration pathways (RCPs), all consistently bias-corrected on a 0.5^∘ × 0.5^∘ global grid. In ISIMIP phase 2, 11 lake models were forced with these data to project the thermal structure of 62 well-studied lakes where data were available for calibration under historical conditions, and using uncalibrated models for 17 500 lakes defined for all global grid cells containing lakes. In ISIMIP phase 3, this approach was expanded to consider more lakes, more models, and more processes. The ISIMIP Lake Sector is the largest international effort to project future water temperature, thermal structure, and ice phenology of lakes at local and global scales and paves the way for future simulations of the impacts of climate change on water quality and biogeochemistry in lakes.

Chromium stable isotope composition (δ⁵³Cr) is a promising tracer for redox conditions throughout Earth's history; however, the geochemical controls of δ⁵³Cr have not been assessed in modern redox-stratified basins. We present new chromium (Cr) concentration and δ⁵³Cr data in dissolved, sinking particulate, and sediment samples from the redox-stratified Lake Cadagno (Switzerland), a modern Proterozoic ocean analog. These data demonstrate isotope fractionation during incomplete (non-quantitative) reduction and removal of Cr above the chemocline, driving isotopically light Cr accumulation in euxinic deep waters. Sediment authigenic Cr is isotopically distinct from overlying waters but comparable to average continental crust. New and published data from other redox-stratified basins show analogous patterns. This challenges assumptions from δ⁵³Cr paleoredox applications that quantitative Cr reduction and removal limits isotope fractionation. Instead, fractionation from non-quantitative Cr removal leads to sedimentary records offset from overlying waters and not reflecting high δ⁵³Cr from oxidative continental weathering.

Propelled by the rapid development of equipment, technology and computational power, the monitoring and simulation of the hydrodynamics in lakes have steadily advanced. In contrast, water quality simulations are more difficult to implement, due to the difficulty in obtaining large-scale, spatially resolved field observations for model validation and the number of interacting processes to be parameterized. Here we demonstrate that remote sensing data can be used to inform Lagrangian particle tracking in a large lake, and vice versa. We used total suspended matter (TSM) as a parameter that can be both estimated from the backscattering in satellite images and modelled in terms of particle abundance. Specifically, we compared TSM concentrations in Lake Geneva deduced from images taken by Sentinel-2 and Sentinel-3 satellites to those estimated from Delft3D hydrodynamic and particle tracking models. TSM concentrations obtained from both methods were compared over a time span of up to 5 days in several scenario studies, including instantaneous and continuous point sources and large-scale TSM simulations. The results demonstrate that remote sensing images can serve to calibrate and validate particle tracking models with independent observations. The model was able to capture both the position of a TSM cloud arising 5 days after an instantaneous point source release, and the direction of particle transport and TSM plume size resulting from a continuous source. Even when simulating the whole lake domain, model results closely approximated the satellite-derived TSM concentrations along lake transects within 9%. In return, the particle tracking model was able to complete partially impaired satellite images, and fill in a four-day image gaps between satellite revisits. The synergy of remote sensing techniques and particle tracking modelling allows a rapid, continuous and more accurate analysis on solute transport in lakes.

The interaction of a uniform cooling rate at the lake surface with sloping bathymetry efficiently drives cross-shore water exchanges between the shallow littoral and deep interior regions. The faster cooling rate of the shallows results in the formation of density-driven currents, known as thermal siphons, that flow downslope until they intrude horizontally at the base of the surface mixed layer. Existing parameterizations of the resulting buoyancy-driven cross-shore transport assume calm wind conditions, which are rarely observed in lakes and thereby restrict their applicability. Here we examine how moderate winds (≲ 5 m s^-1) affect this convective cross-shore transport. We derive simple analytical solutions that we further test against realistic three-dimensional numerical hydrodynamic simulations of an enclosed stratified basin subject to uniform and steady surface cooling rate and cross-shore winds. We show cross-shore winds modify the convective circulation, stopping or even reversing it in the upwind littoral region and enhancing the cross-shore exchange in the downwind region. The analytical parameterization satisfactorily predicted the magnitude of the simulated offshore unit-width discharges in the upwind and downwind littoral regions. Our scaling expands the previous formulation to a regime where both wind and buoyancy forces drive cross-shore discharges of similar magnitude. This range is defined by the non-dimensional Monin-Obukhov length scale, χ_MO: 0.1 ≲ χ_MO ≲ 0.5. The information needed to evaluate the scaling formula can be readily obtained from a traditional set of in-situ observations.

We present a Bayesian inference for a three-dimensional hydrodynamic model of Lake Geneva with stochastic weather forcing and high-frequency observational datasets. This is achieved by coupling a Bayesian inference package, SPUX, with a hydrodynamics package, MITgcm, into a single framework, SPUX-MITgcm. To mitigate uncertainty in the atmospheric forcing, we use a smoothed particle Markov chain Monte Carlo method, where the intermediate model state posteriors are resampled in accordance with their respective observational likelihoods. To improve the uncertainty quantification in the particle filter, we develop a bi-directional long short-term memory (BiLSTM) neural network to estimate lake skin temperature from a history of hydrodynamic bulk temperature predictions and atmospheric data. This study analyzes the benefit and costs of such a state-of-the-art computationally expensive calibration and assimilation method for lakes.

In recent decades, lakes have experienced unprecedented ice loss with widespread ramifications for winter ecological processes. The rapid loss of ice, resurgence of winter biology, and proliferation of remote sensing technologies, presents a unique opportunity to integrate disciplines to further understand the broad spatial and temporal patterns in ice loss and its consequences. Here, we summarize ice phenology records for 78 lakes in 12 countries across North America, Europe, and Asia to permit the inclusion and harmonization of in situ ice phenology observations in future interdisciplinary studies. These ice records represent some of the longest climate observations directly collected by people. We highlight the importance of applying the same definition of ice-on and ice-off within a lake across the time-series, regardless of how the ice is observed, to broaden our understanding of ice loss across vast spatial and temporal scales.

Gravitational convection plays a significant role in the ventilation, heat and mass distribution of aquatic systems. This study investigates thermally driven convection resulting from heat loss at the air-water interface during cooling periods in freshwater environments. In the littoral zone, where the water depth increases from the shoreline to interior waters (pelagic zone), uniform heat loss at the surface generates differential cooling between shallow and deep regions. If the latter process occurs for a long enough time, the density-induced cross-shore pressure gradient may drive an overturning circulation across littoral waters, known as "thermal siphon". This paper examines the conditions under which a thermal siphon develops in natural water bodies and its associated convective regimes. For the above setting, we derive time and velocity scales associated with the transition from Rayleigh-Bénard type convection to horizontal overturning circulation across sloping basin regions. The above transition in the convective regime is characterised by a three-way horizontal momentum balance between the cross-shore pressure gradient and inertia before reaching a quasi-steady regime. Our theoretical scaling expressions are supported by high-fidelity numerical simulations and field-scale experiments, and they provide a robust conceptual framework to characterise convective flows induced by night-time or seasonal surface cooling in nearshore aquatic systems, such as lakes, reservoirs, and coastal seas.

Cooling the surface of freshwater bodies, whose temperatures are above the temperature of maximum density, can generate differential cooling between shallow and deep regions. When surface cooling occurs over a long enough period, the thermally induced cross-shore pressure gradient may drive an overturning circulation, a phenomenon called ‘thermal siphon’. However, the conditions under which this process begins are not yet fully characterised. Here, we examine the development of thermal siphons driven by a uniform loss of heat at the air–water interface in sloping, stratified basins. For a two-dimensional framework, we derive theoretical time and velocity scales associated with the transition from Rayleigh–Bénard type convection to a horizontal overturning circulation across the shallower sloping basin. This transition is characterised by a three-way horizontal momentum balance, in which the cross-shore pressure gradient balances the inertial terms before reaching a quasi-steady regime. We performed numerical and field experiments to test and show the robustness of the analytical scaling, describe the convective regimes and quantify the cross-shore transport induced by thermal siphons. Our results are relevant for understanding the nearshore fluid dynamics induced by nighttime or seasonal surface cooling in lakes and reservoirs.

In stratified lakes, methane oxidizing bacteria are critical methane converters that significantly reduce emissions of this greenhouse gas to the atmosphere. Efforts to better understand their ecology uncovered a surprising diversity, vertical structure, and seasonal succession. It is an open question how this diversity has to be considered in models of microbial methane oxidation. Likewise, it is unclear to what extent simple microbial traits related to the kinetics of the oxidation process and temperature optimum, suggested by previous studies, suffice to understand the observed ecology of methane oxidizing bacteria. Here we incorporate niche partitioning in a mechanistic model of seasonal lake mixing and microbial methane oxidation in a stratified lake. Can we model MOB diversity and niche partitioning based on differences in methane oxidation kinetics and temperature adaptation? We found that our model approach can closely reproduce diversity and niche preference patterns of methanotrophs that were observed in seasonally stratified lakes. We show that the combination of trait values resulting in coexisting methanotroph communities is limited to very confined regions within the parameter space of potential trait combinations. However, our model also indicates that the sequence of community assembly, and variations in the stratification and mixing behavior of the lake result in different stable combinations. A scenario analysis introducing variable mixing conditions showed that annual weather conditions and the pre-existing species also affect the developing stable methanotrophic species composition of the lake. Both, effect of pre-existing species and the environmental impact suggest that the MOB community in lakes may differ from year to year, and a stable community may never truly occur. The model further shows that there are always better-adapted species in the trait parameter space that would destabilize and replace an existing stable community. Thus, natural selection may drive trait values into the specific configurations observed in nature based on physiological limits and tradeoffs between traits.

Water quality measures for inland and coastal waters are available as discrete samples from professional and volunteer water quality monitoring programs and higher-frequency, near-continuous data from automated in situ sensors. Water quality parameters also are estimated from model outputs and remote sensing. The integration of these data, via data assimilation, can result in a more holistic characterization of these highly dynamic ecosystems, and consequently improve water resource management. It is becoming common to see combinations of these data applied to answer relevant scientific questions. Yet, methods for scaling water quality data across regions and beyond, to provide actionable knowledge for stakeholders, have emerged only recently, particularly with the availability of satellite data now providing global coverage at high spatial resolution. In this paper, data sources and existing data integration frameworks are reviewed to give an overview of the present status and identify the gaps in existing frameworks. We propose an integration framework to provide information to user communities through the the Group on Earth Observations (GEO) AquaWatch Initiative. This aims to develop and build the global capacity and utility of water quality data, products, and information to support equitable and inclusive access for water resource management, policy and decision making.

Harmful blooms of the filamentous cyanobacteria Planktothrix rubescens have become common in many lakes as they have recovered from eutrophication over the last decades. These cyanobacteria, capable of regulating their vertical position, often flourish at the thermocline to form a deep chlorophyll maximum. In Lake Zurich (Switzerland), they accumulate during stratified season (May-October) as a persistent metalimnetic thin layer (~2 m wide). This study investigated the role of turbulent mixing in springtime layer formation, its persistence over the summer, and its breakdown in autumn. We characterised seasonal variation of turbulence in Lake Zurich with four surveys conducted in April, July and October of 2018 and September of 2019. Surveys included microstructure profiles and high-resolution mooring measurements. In July and October, the thin layer occurred within a strong thermocline (N ≳ 0.05 s^-1) and withstood significant turbulence, observed as turbulent kinetic energy dissipation rates (ε ≈ 10^-8 W kg^-1). Vertical turbulent overturns -monitored by the Thorpe scale- went mostly undetected and on average fell below those estimated by the Ozmidov scale (L_O ≈ 1 cm). Consistently, vertical diffusivity was close to molecular values, indicating negligible turbulent fluxes. This reduced metalimnetic mixing explains the persistence of the thin layer, which disappears with the deepening of the surface mixed layer in autumn. Bi-weekly temperature profiles in 2018 and a nighttime microstructure sampling in September 2019 showed that nighttime convection serves as the main mechanism driving the breakdown of the cyanobacterial layer in autumn. These results highlight the importance of light winds and convective mixing in the seasonal cycling of P. rubescens communities within a strongly stratified medium-sized lake.

Turbulent mixing controls the vertical transfer of heat, gases and nutrients in stratified water bodies, shaping their response to environmental forcing. Nevertheless, due to technical limitations, the redistribution of wind-derived energy fuelling turbulence within stratified lakes has only been mapped over short (sub-annual) timescales. Here we present a year-round observational record of energy fluxes in the large Lake Geneva. Contrary to the standing view, we show that the benthic layers are the main locus for turbulent mixing only during winter. Instead, most turbulent mixing occurs in the water-column interior during the stratified summer season, when the co-occurrence of thermal stability and lighter winds weakens near-sediment currents. Since stratified conditions are becoming more prevalent -possibly reducing turbulent fluxes in deep benthic environments-, these results contribute to the ongoing efforts to anticipate the effects of climate change on freshwater quality and ecosystem services in large lakes.

We simulated bottom resuspension events in Lake Erie, using a coupled three-dimensional hydrodynamic and water quality model. Key parameters in the model, including critical bottom shear stress (τ_cr) and resuspension rate (α) were calibrated and validated by comparing the model output to observations. These included total suspended solid (TSS) concentrations in the bottom boundary layer (RMSE = 0.74 mgL^-1) and water column (RMSE = 0.81 mgL^-1), and to time series of acoustic backscatter signal (R² > 0.8) and turbidity (R² ≈ 0.4) from long-term moorings near the lakebed in 2008–09 and 2013. Signals from phytoplankton, in spring and summer, caused discrepancies between modeled TSS and the observed turbidity data. Although common practice, we show that literature-based or field-observed critical shear stress should not be directly applied in large-scale Reynolds-averaged sediment model as this will likely underestimate resuspension. In agreement with the literature, the model reproduced more frequent and intensive surface-wave driven resuspension in the shallow regions (< ~20 m), particularly in the western basin, compared to the deeper central and eastern basins, where bottom stresses induced by mean currents (τ_c) were comparable with those due to surface waves (τ_w). However, on the north-shore of the eastern basin, τ_c often predominated over τ_w. We simulated thermocline motion, including up- and down-welling events and swashing of the internal Poincaré wave, to contribute to τ_c in the central basin and form nepheloid layers.

Water-level fluctuation (WLF) is a widespread management action in lakes and reservoirs whose impacts on contaminant fate have seldom been investigated. We used near shore hourly measurements (n = 2122) of turbidity (contaminant proxy) and water velocity (sediment resuspension proxy) to track high-frequency contaminant dynamics during a 0.6 m change in water level observed in autumn 2017 in a large French lake. Simultaneously, discrete trace metal measurements highlighted that trapped sediment was more contaminated and finer than surficial sediment supporting that suspended particles (measured by turbidity) were a preferential medium for contaminant mobility. General additive models involving tensor products revealed the enhancement of wind-speed and river discharge effects on turbidity with water draw down. The decrease of the explained deviances by the models over time-lags indicated short time-scale response of turbidity to external forcing. Three of the four major turbid events occurred at the lowest water-level and were concomitant of sediment resuspension as well as precipitation events and/or river flood suggesting a complex interplay among in-lake and watershed processes at controlling sediment mobility during the WLF. These results shed in light that WLF can affect lake littoral hydrodynamic cascading up to the enhancement of contaminant mobility. Sediment resuspension may be an overlooked feature of WLF increasing contamination risk and exposure for littoral organisms with widespread ecological consequences due to the large number of water-level regulated ecosystems.

This report presents past observations and projects the future development of water temperature in Swiss lakes and rivers. Projections are made until the end of the 21st century using the CH2018 climate scenarios. Besides climate change effects on temperature, we also discuss effects on discharge for rivers, and effects on the thermal structure, and specifically the seasonal mixing regime and ice cover of lakes. [...]

Dieser Bericht beleuchtet die vergangene und zukünftige Entwicklung der Temperaturen in Schweizer Seen und Flüssen. Die Berechnungen bis zum Ende des 21. Jahrhunderts basieren auf den Klimaszenarien CH2018. Neben den Klimaauswirkungen auf die Temperaturen betrachten wir auch die Auswirkungen auf die Abflüsse der Fliessgewässer und auf die thermische Struktur von Seen, insbesondere auf die saisonale Mischung und die Eisbedeckung. [...]

Water inherent optical properties (IOPs) contain integrative information on the optical constituents of surface waters. In lakes, IOP measurements have not been traditionally collected. This study describes how high-frequency IOP profiles can be used to document short-term physical and biogeochemical processes that ultimately influence the long-term trajectory of lake ecosystems. Between October 2018 and May 2020, we collected 1373 high-resolution hyperspectral IOP profiles in the uppermost 50 m of the large mesotrophic Lake Geneva (Switzerland-France), using an autonomous profiler. A data set of this size and content does not exist for any other lake. Results showed seasonal variations in the IOPs, following the expected dynamic of phytoplankton. We found systematic diel patterns in the IOPs. Phases of these diel cycles were consistent year-round, and amplitudes correlated to the diurnal variations of dissolved oxygen, clarifying the link between IOPs and phytoplankton metabolism. Diel amplitudes were largest in spring and summer under low wind condition. Wind-driven changes in thermal stratification impacted the dynamic of the IOPs, illustrating the potential of high-frequency profiles of water optical properties to increase our understanding of carbon cycling in lake ecosystems.

Sediment oxygen demand (SOD) and hypolimnetic oxygen demand (HOD) drive deep-water dissolved oxygen (DO) depletion in lakes, yet these parameters can be difficult to be measure routinely. To address this issue, we present an empirical DO depletion model from time-series measurements of hypolimnetic DO and water-column temperature profiles to estimate hypolimnion thickness (H). The model is based on a dataset that includes 3 temperate lakes (Lake Erie, Lake Simcoe, and Eagle Lake) in Ontario, Canada, of varying size and trophic state. We report SOD (mean [standard deviation]; 0.30 [0.07] g m⁻² d⁻¹) and HOD (0.08 [0.03] g m⁻³ d⁻¹) values based on regression fits of (Formula presented.) (where t is time) and H from these lakes. The model shows that when vertical (through thermocline) and horizontal fluxes can be neglected in the DO budget during summer, SOD and HOD are the first-order parameters driving (Formula presented.). The empirical model predicted hypolimnetic DO in the 3 lakes (root-mean-square error [RMSE] of DO < 3.58 g m⁻³) and was subsequently validated against observations from a fourth lake (Little Silver Lake, Ontario; RMSE of DO = 1.07 g m⁻³). The model provides insight into the importance of physical characteristics (i.e., 1/H) in the hypolimnetic DO budget and the relative impact of physical transport versus biogeochemical sources and sinks in the DO budget. When the spring turnover DO concentration is known, the model can be used to simulate depletion of DO in the hypolimnion, including the onset of hypoxia, using routinely collected temperature profile data. We suggest that the proposed values for SOD and HOD can serve as estimates for water quality model calibration when no information is available.

1. Zooplankton can spend winter actively under the ice cover of lakes. However, dietary resources under lake ice are both quantitatively and qualitatively limited, and feeding might not be energetically rewarding for most zooplankton species. Many zooplankters are expected to fast throughout the winter, exhausting their previously accumulated fat storage. We hypothesised that only a fraction of the actively overwintering zooplankton contributes to an active food web under lake ice, leading to few trophic linkages within the planktonic community.
2. Zooplankton habitats and feeding were investigated under the ice of Lake Onego. Zooplankton habitats and migrations were studied by coupling zooplankton sampling around the clock to measurements of particle movement using an acoustic Doppler current profiler. Secondly, fatty acid-specific stable isotope compositions were used to determine whether and which zooplankton fatty acids ultimately came from the assimilation of under-ice seston.
3. The algal biomass was low under ice and mostly dominated by large diatoms. Copepods dominated the zooplankton community. Species present as late copepodite and adult instars were confined to the deeper layers, while nauplii occupied the surface layer. Diel vertical migration by Cyclops was the most tangible observation of persistent feeding under the ice. Previously accumulated fat storage represented most of zooplankton fatty acids, with few, yet detectable, fatty acids recently acquired by feeding under the ice.
4. Although some zooplankton taxa maintained feeding activity under the ice of Lake Onego, the food source available beneath the ice was not sufficiently rewarding to leave an isotopic imprint upon the dominant fatty acids of bulk zooplankton. The seston fatty acids that were passed on to zooplankton from feeding under ice were not provided by diatoms, although they made up most of the phytoplankton biovolume. Instead, the zooplankton food web was supported by mixotrophic phytoplankton (i.e. cryptophytes and chrysophytes) that represented <5% of the under-ice biovolume. Consequently, the planktonic food web under the ice of Lake Onego had few trophic linkages, and thereby low connectance.
5. Environmental conditions under the ice of Lake Onego do not depart significantly from those observed in lakes of similar latitudes. Therefore, low connectance food webs could be relatively common under ice for lakes above 60° latitude. Our methodological approach is applicable to other lakes and could thus disclose the variability of under-ice food webs. This would provide a much more complete picture of annual dynamics of food webs in lakes that ice over.

The gas transfer velocity (k) is a major source of uncertainty when assessing the magnitude of lake gas exchange with the atmosphere. For the diversity of existing empirical and process-based k models, the transfer velocity increases with the level of turbulence near the air–water interface. However, predictions for k can vary by a factor of 2 among different models. Near-surface turbulence results from the action of wind shear, surface waves, and buoyancy-driven convection. Wind shear has long been identified as a key driver, but recent lake studies have shifted the focus towards the role of convection, particularly in small lakes. In large lakes, wind fetch can, however, be long enough to generate surface waves and contribute to enhance gas transfer, as widely recognised in oceanographic studies. Here, field values for gas transfer velocity were computed in a large hard-water lake, Lake Geneva, from CO₂ fluxes measured with an automated (forced diffusion) flux chamber and CO₂ partial pressure measured with high-frequency sensors. k estimates were compared to a set of reference limnological and oceanic k models. Our analysis reveals that accounting for surface waves generated during windy events significantly improves the accuracy of k estimates in this large lake. The improved k model is then used to compute k over a 1-year time period. Results show that episodic extreme events with surface waves (6 % occurrence, significant wave height > 0.4 m) can generate more than 20 % of annual cumulative k and more than 25 % of annual net CO₂ fluxes in Lake Geneva. We conclude that for lakes whose fetch can exceed 15 km, k models need to integrate the effect of surface waves.

In late winter, solar radiation is the main driver of water motion in ice-covered lakes. The resulting circulation and mixing determine the spatial distribution of heat within the lake and affect the heat budget of the ice cover. Although under-ice lake warming is often modeled as a one-dimensional (1D) vertical process, lake bathymetry induces a relative excess heating of shallow waters, creating horizontal density gradients. This study shows that the dynamic response to these gradients depends sensitively on lake size and latitude - Earth's rotation - and is controlled by the Rossby number. In the ageostrophic limit, horizontal density gradients drive cross-shore circulation that transports excess heat to the lake interior, accelerating the under-ice warming there. In the geostrophic regime, the circulation of the near- and off-shore waters decouples, and excess heat is retained in the shallows. The flow regime controls the fate of this excess heat and its contribution to water-induced ice melt.

Anomaly detection is the process of identifying unexpected data samples in datasets. Automated anomaly detection is either performed using supervised machine learning models, which require a labelled dataset for their calibration, or unsupervised models, which do not require labels. While academic research has produced a vast array of tools and machine learning models for automated anomaly detection, the research community focused on environmental systems still lacks a comparative analysis that is simultaneously comprehensive, objective, and systematic. This knowledge gap is addressed for the first time in this study, where 15 different supervised and unsupervised anomaly detection models are evaluated on 5 different environmental datasets from engineered and natural aquatic systems. To this end, anomaly detection performance, labelling efforts, as well as the impact of model and algorithm tuning are taken into account. As a result, our analysis reveals the relative strengths and weaknesses of the different approaches in an objective manner without bias for any particular paradigm in machine learning. Most importantly, our results show that expert-based data annotation is extremely valuable for anomaly detection based on machine learning.

Studies of future 21^st century climate warming in lakes along altitudinal gradients have been partially obscured by local atmospheric phenomena unresolved in climate models. Here we forced the physical lake model Simstrat with locally downscaled climate models under three future scenarios to investigate the impact on 29 Swiss lakes, varying in size along an altitudinal gradient. Results from the worst-case scenario project substantial change at the end of the century in duration of ice-cover at mid to high altitude (−2 to −107 days), stratification duration (winter −17 to −84 days, summer −2 to 73 days), while lower and especially mid altitude (present day mean annual air temperature from 9 °C to 3 °C) dimictic lakes risk shift to monomictic regimes (seven out of the eight lakes). Analysis further indicates that for many lakes shifts in mixing regime can be avoided by adhering to the most stringent scenario.

In situ observations of biophysical interactions in natural waters typically focus on physical mechanisms influencing biological activity. Yet, biological activity can also drive physical processes in aquatic environments. A community of photoautotrophic, motile and heavy bacteria—Chromatium okenii, which requires light, sulfide, and anoxic conditions to perform anoxygenic photosynthesis, accumulates below the chemocline of the meromictic Lake Cadagno (Switzerland). Upward vertical migration drives bioconvection, which modifies the physical environment of the bacteria‐populated water to create a deep and homogeneous mixed layer of up to 1 m thickness. Continuous convection within the mixed layer and diapycnal diffusivity from its adjacent stratified surroundings determine ecologically relevant gradients. The daytime vertical migration that induce convective instabilities is well‐established. It consists in bacteria swimming upward towards light and accumulating at the upper part of the anoxic layer, leading to a locally‐unstable density excess. However, nocturnal activity has not yet been analyzed. An intensive 48‐h survey was conducted in August 2018 using standard and microstructure profilers, as well as a moored high‐resolution current meter coupled with temperature and turbidity sensors deployed across the mixed layer depth. This survey revealed a persistent mixed layer also during nighttime hours. Using a mixed layer shape model, vertical velocity observations and turbulent dissipation estimates, we conclude that photoautotrophic bacteria continue their vertical migration at night. This nocturnal activity thereby drives "dark bioconvection" and maintains the subsurface mixed bacterial layer in Lake Cadagno throughout the diel cycle.

Long-term lake ice phenological records from around the Northern Hemisphere provide unique sensitive indicators of climatic variations, even prior to the existence of physical meteorological measurement stations. Here, we updated ice phenology records for 60 lakes with time-series ranging from 107–204 years to provide the first re-assessment of Northern Hemispheric ice trends since 2004 by adding 15 additional years of ice phenology records and 40 lakes to our study. We found that, on average, ice-on was 11.0 days later, ice-off was 6.8 days earlier, and ice duration was 17.0 days shorter per century over the entire record for each lake. Trends in ice-on and ice duration were six times faster in the last 25-year period (1992–2016) than previous quarter centuries. More extreme events in recent decades, including late ice-on, early ice-off, shorter periods of ice cover, or no ice cover at all, contribute to the increasing rate of lake ice loss. Reductions in greenhouse gas emissions could limit increases in air temperature and abate losses in lake ice cover that would subsequently limit ecological, cultural, and socioeconomic consequences, such as increased evaporation rates, warmer water temperatures, degraded water quality, and the formation of toxic algal blooms.

Predicting the freezing time in lakes is achieved by means of complex mechanistic models or by simplified statistical regressions considering integral quantities. Here, we propose a minimal model (SELF) built on sound physical grounds that focuses on the pre-freezing period that goes from mixed conditions (lake temperature at 4 ^∘C) to the formation of ice (0 ^∘C at the surface) in dimictic lakes. The model is based on the energy balance involving the two main processes governing the inverse stratification dynamics: cooling of water due to heat loss and wind-driven mixing of the surface layer. They play opposite roles in determining the time required for ice formation and contribute to the large interannual variability observed in ice phenology. More intense cooling does indeed accelerate the rate of decrease of lake surface water temperature (LSWT), while stronger wind deepens the surface layer, increasing the heat capacity and thus reducing the rate of decrease of LSWT. A statistical characterization of the process is obtained with a Monte Carlo simulation considering random sequences of the energy fluxes. The results, interpreted through an approximate analytical solution of the minimal model, elucidate the general tendency of the system, suggesting a power law dependence of the pre-freezing duration on the energy fluxes. This simple yet physically based model is characterized by a single calibration parameter, the efficiency of the wind energy transfer to the change of potential energy in the lake. Thus, SELF can be used as a prognostic tool for the phenology of lake freezing.

Environmental sciences depend heavily on observational data. Successful studies of ecological processes in lakes require in-situ data that cover the relevant temporal scales from milliseconds to entire seasons. Temporal and spatial coverage requirements represent a non-trivial challenge in lake sciences, which have traditionally used sampling campaigns conducted from research vessels or anchored moorings. These come with various logistical tasks and impose constraints on data coverage. An open water platform can overcome many of these limitations by providing continuous access and a wide range of analytical capabilities in direct contact with the lake environment. A consortium of five partner institutions constructed a 10 × 10 m, open-water, multipurpose platform on Lake Geneva (Switzerland/France) for a broad range of limnological research. The LéXPLORE platform, anchored since February 2019 at a position reaching 110 m depth off the lake's north-shore, provides workspace for a large number of instruments and up to 16 staff working in parallel on individual or integrated multidisciplinary projects. The safe, dry and protected floating laboratory offers direct access to the lake environment for high-sensitivity, high-throughput analyses including those which might advance sensor technology. The platform provides flexible workspace for both high-resolution measurements and investigations of larger-scale external forcing. It thus supports multidisciplinary empirical research in limnology, atmospheric sciences, and remote sensing. This article describes the platform and how it will advance aquatic sciences. The large number of projects that have already requested access to the platform demonstrate the efficacy and necessity of the LéXPLORE concept.

Many seasonally stratified lakes accumulate substantial amounts of the greenhouse gas methane in the anoxic zone. Methane oxidizing bacteria in the water column act as a converter, oxidizing methane into carbon dioxide and biomass before it reaches the atmosphere. Current observations and estimates of this methane oxidation efficiency are diverging, especially for the lake overturn period. Here, we combine a model of turbulent mixing, gas exchange, and microbial growth with a comprehensive data set for autumn mixing to quantify the relevant physical and microbial processes for a 16 m deep, wind-sheltered Swiss lake. Scenario analysis suggests that the methane converter is efficient and robust under a large range of mixing velocities and only rare events of pronounced surface cooling can trigger substantial outgassing. This case study combines in situ observation and a deterministic physically based model and suggests that the frequency of storms may strongly impact methane emissions for similar temperate lakes.

Understanding lake dynamics is crucial to provide scientifically credible information for ecosystem management. In this context, three-dimensional hydrodynamic models are a key information source to assess critical but often subtle changes in lake dynamics occurring at all spatio-temporal scales. However, those models require time-consuming calibrations, often carried out by trial-and-error. Through a new coupling of open source software, we present here a flexible and computationally inexpensive automated calibration framework. The method, tailored to the calibration data available to the user, aims at (i) reducing the time spent on calibration, and (ii) making three-dimensional lake modelling accessible to a broader range of users. It is demonstrated for two different lakes (Lake Geneva and Greifensee) with an extensive multi-variable observational dataset. Models mean absolute errors are reduced by up to ~50% over the baseline. Guidelines on heat and momentum transfer parameters are given with their dependence on the observational setup.

The understanding of physical dynamics is crucial to provide scientifically credible information on lake ecosystem management. We show how the combination of in situ observations, remote sensing data, and three-dimensional hydrodynamic (3D) numerical simulations is capable of resolving various spatiotemporal scales involved in lake dynamics. This combination is achieved through data assimilation (DA) and uncertainty quantification. In this study, we develop a flexible framework by incorporating DA into 3D hydrodynamic lake models. Using an ensemble Kalman filter, our approach accounts for model and observational uncertainties. We demonstrate the framework by assimilating in situ and satellite remote sensing temperature data into a 3D hydrodynamic model of Lake Geneva. Results show that DA effectively improves model performance over a broad range of spatiotemporal scales and physical processes. Overall, temperature errors have been reduced by 54 %. With a localization scheme, an ensemble size of 20 members is found to be sufficient to derive covariance matrices leading to satisfactory results. The entire framework has been developed with the goal of near-real-time operational systems (e.g., integration into meteolakes.ch).

Environmental management depends on high-quality monitoring and its meaningful interpretation. The combination of local weather dynamics, regional anthropogenic stresses and global environmental changes make the evaluation of monitoring information in dynamic freshwater systems a challenging task. While the lake ecosystems gather many complex biogeochemical interactions, they remain constrained by the same physical environment of mixing and transport. It is therefore crucial to obtain high-quality physical system insight. Three-dimensional hydrodynamic models are perfectly suited for providing such information. However, these models are complex to implement, and their use is often limited to modellers. Here, we aim to provide model output via a user-friendly platform to a broad audience ranging from scientists to public and governmental stakeholders.
We present a unified approach merging the apparently diverse interests through meteolakes.ch, an online platform openly disseminating lake observations and three-dimensional numerical simulations in near real-time with short-term forecasts and data assimilation. Meteolakes is scalable to a broad range of devices, modular and distributed, hence allowing its expansion to other regions and hardware infrastructures. Since 2016, the platform has continuously provided timely synoptic lake information to more than 250,000 users. This web-based system was built not only to provide guidance to scientists in the design and analysis of field experiments and to foster interdisciplinary lake studies, but also to assist governmental agencies and professionals in the long-term policy and planning of water resources management. Finally, our system aimed at promoting awareness and understanding of the complexity of lakes and providing information to the public through user-friendly interfaces. This article details the design and operation of such a platform and its products. Applications are demonstrated by examples of a recent upwelling and a storm event. Both cases illustrate how Meteolakes help scientists in their quest for process understanding as well as water professionals and civil society in providing specific warnings.

Through the artificial manipulation of hydrological residence time, hydropower operations on alpine lakes should not only alter their hydro-biogeochemical characteristics but also modulate the degree of sensitivity of the lake to atmospheric forcing. These hypotheses were addressed by means of in situ observations that were combined with a three-dimensional lake model of a high-altitude pump-storage reservoir, i.e., Lake Corne, during the ice-free period. The fitted model (Pumping scenario) was used to simulate the hydrodynamics and oxygen concentrations of the lake under natural conditions (without pump-storage). Thereafter, the lake response to direct or catchment-mediated meteorological conditions (e.g., changes in inflows during storms) was investigated with a sensitivity analysis of the forcing parameters under both the pumping and natural hydrological scenarios. The pumping operation resulted in significant changes in water temperature, heat content and water mass stability during summer, as compared to a natural scenario. The lake hydrodynamics during the ice-free season were highly responsive to changes in the summer weather conditions, and more sensitive to realistic changes in cloudiness and water transparency than changes in air temperatures for both the pumped and natural scenarios. The spatial and temporal evolution of dissolved oxygen was comparatively less responsive to summer weather conditions, in both scenarios. The pump-storage operation from Lake Corne had a comparatively smaller effect on the lake functioning than the natural meteorological variability in summer. However, the pump-storage operation decreased the sensitivity of the lake hydrodynamics to changes in water transparency and limited water mass stability during summer.

Mixing rates and biogeochemical fluxes are commonly estimated from the rate of dissipation of turbulent kinetic energy ε as measured with a single instrument and processing method. However, differences in measurements of ε between instruments/methods often vary by one order of magnitude. In an effort to identify error in computing ε, we have applied four common methods to data from the bottom boundary layer of Lake Erie. We applied the second-order structure function method (SFM) to velocity measurements from an acoustic Doppler current profiler, using both canonical and anisotropy-adjusted Kolmogorov constants, and compared the results with those computed from the law of the wall, Batchelor fitting to temperature gradient microstructure, and inertial subrange fitting to acoustic Doppler velo-cimeter data. The ε from anisotropy-adjusted constants in SFM increased by a factor of 6 or more at 0.2 m above the bed and showed a better agreement with microstructure and inertial method estimations. The maximum difference between SFM ε, computed using adjusted and canonical constants, and microstructure values was 25% and 50%, respectively. This difference was 30% and 55%, respectively, for those from inertial subrange fitting at times of high-intensity turbulence (Reynolds number at 1 m above the bed of more than 2 x 10⁴). Comparison of the SFM ε to those from law of the wall was often poor, with errors as large as one order of magnitude. From the considerable improvement in ε estimates near the bed, anisotropy-adjusted Kolmogorov constants should be applied to compute dissipation in geophysical boundary layers.

Large lakes of the world are habitats for diverse species, including endemic taxa, and are valuable resources that provide humanity with many ecosystem services. They are also sentinels of global and local change, and recent studies in limnology and paleolimnology have demonstrated disturbing evidence of their collective degradation in terms of depletion of resources (water and food), rapid warming and loss of ice, destruction of habitats and ecosystems, loss of species, and accelerating pollution. Large lakes are particularly exposed to anthropogenic and climatic stressors. The Second Warning to Humanity provides a framework to assess the dangers now threatening the world's large lake ecosystems and to evaluate pathways of sustainable development that are more respectful of their ongoing provision of services. Here we review current and emerging threats to the large lakes of the world, including iconic examples of lake management failures and successes, from which we identify priorities and approaches for future conservation efforts. The review underscores the extent of lake resource degradation, which is a result of cumulative perturbation through time by long-term human impacts combined with other emerging stressors. Decades of degradation of large lakes have resulted in major challenges for restoration and management and a legacy of ecological and economic costs for future generations. Large lakes will require more intense conservation efforts in a warmer, increasingly populated world to achieve sustainable, high-quality waters. This Warning to Humanity is also an opportunity to highlight the value of a long-term lake observatory network to monitor and report on environmental changes in large lake ecosystems.

With the increasing extraction of methane from Lake Kivu, there is a growing need to evaluate the effect of such operations on the lake's permanent density stratification. This requires understanding of the spatial structure and variability of flow velocities and constituents in Lake Kivu. In this study, we develop a 3D hydrodynamic model of Lake Kivu, set-up within DELFT3D at a 750 m grid spacing and forced by COSMO-CLM atmosphere model results at a 2.8 km grid spacing. Validation shows that the model correctly reproduces the generation and breakdown of the temperature stratification in the upper mixed layer and predicts flow velocity magnitudes and directions similar to measurements both at the surface and at greater depth. Analysis of currents reveals a surface current pattern with two clockwise circulations, one around the whole lake and a smaller one in the northern part, with velocities around 0.1 m/s. This pattern is consistently present over an (ensemble-)averaged day, both in the wet and in the dry season, while day-by-day variations are large. Time-averaged deep currents are found to be a few mm/s at maximum. However, the variations can be substantial, with standard deviations up to 2 cm/s for the currents at 220 m depth, attributed to internal seiches. The temperature stratification, present during the entire wet season, is found to first break down in the dry season in the southern part of the lake. This is explained by the spatial differences in the wind stress and the evaporation heat fluxes during the dry season.

In March 2017, we joined a multi‐disciplinary expedition on Lake Onego, in the Russian Karelia, with the primary goal to investigate zooplankton strategies to survive deep winter under ice. We intended to focus on the large zooplankters, i.e. crustaceans of tenths to a few mm of body length, through sampling by horizontal and vertical trawls using nets of relatively large mesh‐size. Delays in the delivery of sampling gear compelled us to opt, at the last minute, for a much‐smaller mesh size net for the horizontal trawls underneath ice.

High-resolution lake and reservoir bathymetric surveys can pinpoint locations that may experience underwater landslides (subaquatic sedimentary mass movements). These can pose a risk to underwater and shoreline infrastructure. This paper outlines an approach for using spatial variation in sedimentary patterns to identify areas susceptible to subaquatic mass movements in lakes and reservoirs. This study focusses on Lake Biel (Switzerland), which has experienced a protracted history of upstream alteration of river flow. Altered flow patterns increase risk of unstable sedimentary features and subaquatic mass movements. Data from sediment traps and cores, Acoustic Doppler Current Profilers and results from a 3D hydrodynamic model gave a consistent picture of spatial and temporal variation in weather-related sedimentation. Erosion caused by short-term rain events contributes the largest proportion of sediments to the lake. Strong rain events combine with typical wind patterns to drive lake circulation. The net effect results in preferential sedimentation onto a steeply sloping shelf prone to subaquatic slides. The integrated approach outlined here incorporates short- and long-term sediment dynamics to provide a systematic assessment of lake sedimentation and potential mass movement hazards. This research represents a first step in developing a risk-evaluation tool for aquatic hazard evaluation.

Pockmarks are circular depressions on the floor of oceans and lakes and constitute potential hot spots of gas ebullition and/or groundwater discharge. Marine pockmarks are well-studied, whereas lacustrine pockmarks are virtually unexplored. In Lake Neuchatel (Switzerland) four giant pockmarks of 80 to 150 m in diameter are located along the northern shore and adjacent to the karst system of the Jura Mountains. Two pockmarks have a ~60 m-deep chimney filled with mud; two are funnel-shaped 12 and 29 m deep holes. We present evidence for the presence of groundwater in the pockmark chimneys and active lacustrine groundwater discharge (LGD) at both pockmark types. Temperature, electrical conductivity and calcium concentrations of the pore water in the chimneys show values typical for karst water (maximal sampling and profiling depth: 41 m into the pockmark) and contrast strongly with the properties of the lake water. TOC and TIC indicate that the chimney mud is homogenized and composed of liquefied sediments from the entire deglacial to Holocene lacustrine sediment succession. Mini mud volcanoes apparent on the suspension surface imaged with a remotely operated vehicle (ROV) localize the groundwater exit points and confirm LGD. LGD is further corroborated by electrical conductivity anomalies detected above the lutoclines and within a funnel-shaped pockmark during the ROV survey. We conclude that the giant pockmarks in Lake Neuchatel represent a type of subaquatic springs that connect the water body of the lake with the karst system. A next essential research step will be to quantify LGD via the pockmarks in order to assess their lake-wide relevance. This study underlines the existing need for research on the connectivity of lakes and oceans with groundwater systems for completing our understanding of the hydrological cycle.

Les lacs fournissent de nombreux services écosystémiques, mais sont menacés à la fois à l'échelle globale et locale. De nouvelles approches, permettant de comprendre et d'anticiper leurs réactions aux pressions externes (effets anthropogéniques locaux et globaux), sont nécessaires. Cet article présente un nouveau modèle hydrodynamique 3D du lac de Zurich avec la plateforme en ligne: meteolakes.ch. Il offre ainsi une nouvelle perspective sur la variabilité spatio-temporelle du lac et son management.

Hydrophysical studies conducted in Petrozavodsk Bay of Lake Onego under ice-covered conditions in March 2016 and 2017 detected radiatively driven convection and revealed specific structural and dynamic parameters for the convectively mixed layer (CML). Analysis of time series, spectral energy distribution, and vertical velocity profiles indicated the presence of a mean current, seiches, and convective motion. Because of their similar spatial and temporal scales, these processes were investigated using progressive-vector diagrams (PVDs). Despite low water velocities, the CML hydrodynamic regime remained close to that of fully developed turbulence, and convective cells at a range of different scales were expected. Signal resolution constraints limited detection to only the largest cells. We investigated the horizontal structure of the CML using individual and combined observations from 3 acoustic velocity profilers located within a radius of a few tens of meters. This novel setup collected data indicating that the CML's large-scale horizontal flow structure consists of a continuum of quasi-deterministic cells. Cell parameters necessary for estimating turbulent transfer were derived from PVD and hodograph curves.

Lakes and other confined water bodies are not exposed to tides, and their wind forcing is usually much weaker compared to ocean basins and estuaries. Hence, convective processes are often the dominant drivers for shaping mixing and stratification structures in inland waters. Due to the diverse environments of lakes—defined by local morphological, geochemical, and meteorological conditions, among others—a fascinating variety of convective processes can develop with remarkably unique signatures. Whereas the classical cooling-induced and shear-induced convections are well-known phenomena due to their dominant roles in ocean basins, other convective processes are specific to lakes and often overlooked, for example, sidearm, under-ice, and double-diffusive convection or thermobaric instability and bioconvection. Additionally, the peculiar properties of the density function at low salinities/temperatures leave distinctive traces. In this review, we present these various processes and connect observations with theories and model results.

Climate change and other anthropogenic factors affect lakes in many ways that can lead to important effects on freshwater ecosystems. For this reason, the Swiss Federal Office for the Environment (FOEN) is currently evaluating options to add a nation-wide lake temperature monitoring program to their monitoring activities, which already includes a network of river temperature monitoring.

In this work, we discuss the benefits and limitations of temperature data obtained by a monitoring program depending on the spatial and temporal resolution of the measurements. We conclude that for accurately observing trends in the thermal structure of lakes caused by climate change, a high time resolution is required that cannot be provided by traditional temperature profiling. We therefore recommend installing moorings where temperature can be continuously monitored.

Individual lakes react differently to climate change according to internal and external parameters, such as their size, trophic state, or altitude. These different reactions should be considered when selecting lakes to be included in the monitoring network. We therefore discuss the factors that influence a lake’s reaction to climate change, and present four different scenarios for a monitoring network. The first scenario includes the modelling of lakes that are currently monitored, without adding new monitoring sites. The second scenario represents the installation of moorings in all lakes above a certain size. The third scenario is the monitoring of lakes at different altitudes. Scenario 4 is a combination of Scenario 2 and Scenario 3 with the addition of other lakes with various characteristics to broaden the perspective. The advantages and disadvantages of the four scenarios are assessed.

Finally, we discuss the opportunity of measuring other parameters (e.g., meteorological forcing, chemical parameters…) that would increase the coherence and impacts of the monitoring program. We also shortly discuss options for data management and data quality control.

We investigated radiatively driven under-ice convection in Lake Onego (Russia) during 3 consecutive late winters. In ice-covered lakes, where the temperature of water is below the temperature of maximum density, radiatively driven heating in the upper water column induces unstable density distributions leading to gravitational convection. In this work, we quantified the key parameters to characterise the radiatively driven under-ice convection: (1) the effective buoyancy flux, B_* (driver), and its vertical distribution; (2) the convective mixed-layer thickness, h_CML (depth scale); and (3) the convective velocity,w_*(kinematic scale). We compared analytical w_* scaling estimates to in situ observations from high-resolution acoustic Doppler current profilers. The results show a robust correlation between w_* and the direct observations, except during the onset and decay of the solar radiation. Our results highlight the importance of accurately defining the upper limit of h_CML in highly turbid water and the need for spectrally resolving solar radiation measurements and their attenuation for accurate B_* estimates. Uncertainties in the different parameters were also investigated. We finally examined the implications of under-ice convection for the growth rate of nonmotile phytoplankton and provide a simple heuristic model as a function of easily measurable parameters.

One-dimensional hydrodynamic models are nowadays widely recognized as key tools for lake studies. They offer the possibility to analyze processes at high frequency, here referring to hourly timescales, to investigate scenarios and test hypotheses. Yet, simulation outputs are mainly used by the modellers themselves and often not easily reachable for the outside community. We have developed an open-access web-based platform for visualization and promotion of easy access to lake model output data updated in near-real time (http://simstrat.eawag.ch, last access: 29 August 2019). This platform was developed for 54 lakes in Switzerland with potential for adaptation to other regions or at global scale using appropriate forcing input data. The benefit of this data platform is practically illustrated with two examples. First, we show that the output data allows for assessing the long-term effects of past climate change on the thermal structure of a lake. The study confirms the need to not only evaluate changes in all atmospheric forcing but also changes in the watershed or throughflow heat energy and changes in light penetration to assess the lake thermal structure. Then, we show how the data platform can be used to study and compare the role of episodic strong wind events for different lakes on a regional scale and especially how their thermal structure is temporarily destabilized. With this open-access data platform, we demonstrate a new path forward for scientists and practitioners promoting a cross exchange of expertise through openly sharing in situ and model data.

Remote sensing helps foster our understanding of inland water processes allowing a synoptic view of water quality parameters. In the context of global monitoring of inland waters, we demonstrate the benefit of combining in-situ water analysis, hydrodynamic modelling and remote sensing for investigating biogeochemical processes. This methodology has the potential to be used at global scales. We take the example of four Landsat-8 scenes acquired by the OLI sensor and MODIS-Aqua imagery over Lake Geneva (France—Switzerland) from spring to early summer 2014. Remotely sensed data suggest a strong temporal and spatial variability during this period. We show that combining the complementary spatial, spectral and temporal resolutions of these sensors allows for a comprehensive characterization of estuarine, littoral and pelagic near-surface features. Moreover, by combining in-situ measurements, biogeochemical analysis and hydrodynamic modelling with remote sensing data, we can link these features to river intrusion and calcite precipitation processes, which regularly occur in late spring or early summer. In this context, we propose a procedure that can be used to monitor whiting events in temperate lakes worldwide.

Ice-covered periods might significantly contribute to lake emissions at ice-melt, yet a comprehensive understanding of under-ice carbon dioxide (CO₂) dynamics is still lacking. This study investigated the processes driving spatiotemporal patterns of under-ice CO₂ in large Lake Onego. In March 2015 and 2016, under-ice CO₂, dissolved inorganic carbon (DIC), and dissolved organic carbon (DOC) distributions were measured along a river to an open-lake transect. CO₂ decreased from 120/129 μmol L⁻¹ in the river to 51/98 μmol L⁻¹ in the bay, and 34/36 μmol L⁻¹ in the open lake, while DOC decreased from 1.18/1.55 mmol L⁻¹ in the river to 0.67/1.04 mmol L⁻¹ in the bay in 2015 and 2016, respectively. These decreases in concentrations with increasing distance from the river mouth indicate that river discharge modulates spatial patterns of under-ice CO₂. The variability between the 2 years was mainly driven by river discharge and ice transparency affecting the extent of under-ice convection. Higher discharge during winter 2016 resulted in higher CO₂ concentrations in the bay. By contrast, intensive under-ice convection led to lower, more homogeneously distributed CO₂ in 2015. In conclusion, the river-to-bay transition zone is characterized by strong CO₂ variability and is therefore an important zone to consider when assessing the CO₂ budget of large lakes.

In stratified natural waters, convective processes tend to form nearly homogeneous mixed layers. However, shear‐driven turbulence generated by large‐scale background flow often rapidly smooths them through mixing with the stratified surroundings. Here we studied the effect of background turbulence on convectively driven mixed layers for the case of bioconvection in Lake Cadagno, Switzerland. Along with microstructure measurements, a diffusive‐shape model for the mixed layers allowed us to define (i) mixed layer thickness and (ii) diffusive transition length. Further microstructure analysis was performed allowing estimation of convective turbulence in the mixed layer and shear‐driven turbulence quantified by eddy diffusion in their surroundings. Based upon these results, we propose a Péclet number scaling that relates mixed layer shape to the opposing effects of convection and diffusion. We further validate this quantitative approach by applying it to two other distinct convective systems representative of double‐diffusive convection and radiatively driven under‐ice convection.

Lake productivity is fundamental to biogeochemical budgets as well as estimating ecological state and predicting future development. Combining modelling with Earth Observation data facilitates a new perspective for studying lake primary production. In this study, primary production was modelled in the large Lake Geneva using the MEdium Resolution Imaging Spectrometer (MERIS) image archive for 2002–2012. We used a semi-empirical model that estimates primary production as a function of photosynthetically absorbed radiation and quantum yield of carbon fixation. The necessary input parameters of the model—concentration of chlorophyll a, downwelling irradiance, and the diffuse attenuation coefficient—were obtained from MERIS products. The primary production maps allow us to study decennial temporal (with daily frequency) and spatial changes in this lake that a single sample point cannot provide. Modelled estimates agreed with in situ results (R² = 0.68) and showed a decreasing trend (∼27%) in production in Lake Geneva for the selected decade. Yet, in situ monitoring measurements missed the general increase of productivity near the incoming Rhône River. We show that the temporal and spatial resolution provided by satellite observations allows estimates of primary production at the basin-scale. The phytoplankton annual primary production was estimated as ∼302 (SD 20) g C m⁻² yr⁻¹ for Lake Geneva for 2003 to 2011. This study demonstrates that maps of primary production can be obtained even with reduced resolution (1200  m) MERIS data and relatively simple methods, and thereby calls for deeper integration of remote sensing products into conventional in situ observation approaches.

This study deals with the impact of spatio-temporal heterogeneities on the assessment of lake ecological status according to the European water framework directive (WFD). A method, based on three-dimensional coupled hydrodynamic and ecological modeling, is presented to assess the variability of lake ecological status, and to locate the most representative sampling station of Lake Geneva (France/Switzerland). Five variables used in the lake ecological status evaluation were simulated by using the free software Delft3D. The numerical simulation results showed that the simulated ecological status based on chlorophyll a and total phosphorus concentrations measured at the regulatory monitoring station depend on the choice of the sampling date. Results also indicated a strong spatial heterogeneity in ecological status that varies from "poor" to "good" along an East-West gradient. Finally, the numerical simulation results showed that the most representative point of a mean theoretical ecological quality for Lake Geneva would be located in the center of the upper basin, close to the historical sampling station.

Cette étude traite de l'impact des hétérogénéités spatio-temporelles sur l'évaluation de l'état écologique des lacs conformément à la directive-cadre européenne sur l'eau (DCE). Une méthode, basée sur la modélisation tridimensionnelle couplée hydrodynamique et écologique, est présentée pour évaluer la variabilité de l'état écologique et pour localiser la station d'échantillonnage la plus représentative du lac Léman. Cinq variables utilisées dans l'évaluation de l'état écologique ont été simulées en utilisant le logiciel libre Delft3D. Les résultats des simulations numériques ont montré que l'état écologique simulé basé sur les concentrations en chlorophylle a et en phosphore total mesurées à la station de suivi réglementaire dépend du choix de la date d'échantillonnage. Les résultats ont aussi indiqué des fortes hétérogénéités spatiales de l'état écologique allant de «médiocre» à «bon» le long d'un gradient est-ouest. Finalement, les résultats de simulations numériques ont montré que le point le plus représentatif d'un état écologique théorique moyen du lac Léman serait situé au centre du grand bassin, proche de la station d'échantillonnage historique.

Various lake observables, including lake ice, are related to climate and climate change and provide a good opportunity for long-term monitoring. Lakes (and as part of them lake ice) is therefore considered an Essential Climate Variable (ECV) (WMO, 2018) of the Global Climate Observing System (GCOS). In Switzerland, the implementation of GCOS is coordinated by the Swiss GCOS Office at the Federal Office for Meteorology and Climatology MeteoSwiss. In 2007 (updated in 2018), MeteoSwiss published the first national inventory of the most important climate observations in Switzerland (MeteoSwiss, 2018). For each ECV, including observations of lake ice subsumed under the ECV "Lakes", and international centre, the inventory report identified possible gaps regarding the legal basis, definition of responsibilities, and availability of financial resources for the continuation of observations and operation, respectively. In 2008, considering the findings of this report, the Federal Council approved a financial contribution to secure the continuation of time series and international centres at risk of discontinuation, as a long-term contribution to GCOS. Concerning the ECV "Lakes", in the global National Snow and Ice Data Center (NSIDC), Boulder, Colorado database on lake ice-on/off dates, currently only Lake St. Moritz is partly included (till 2012). Further, existing observations and data from local authorities and publications are not systematic and come from different, uncoordinated and not secured sources. Traditionally, on-shore observers collected the information of lake ice recording the visible frozen-edge. Within the past two decades due to lack of budget and/or human resources, the number of field observations declined and mostly totally stopped. At the same time, the potential of different remote sensing sensors covering varying time periods and spatial coverage to measure the occurrence of lake ice was demonstrated by several investigations. Thus, following the need for an integrated multi-temporal monitoring of lake ice in Switzerland, MeteoSwiss in the framework of GCOS Switzerland supported this 2-year project to explore not only the use of satellite images but also the possibilities of Webcams and in-situ measurements. [...]

In ice‐covered lakes, penetrative radiation warms fluid beneath a diffusive boundary layer, thereby increasing its density and providing energy for convection in a diurnally active, deepening mixed layer. Shallow regions are differentially heated to warmer temperatures, driving turbulent gravity currents that transport warm water downslope and into the basin interior. We examine the energetics of these processes, focusing on the rate at which penetrative radiation supplies energy that is available to drive fluid motion. Using numerical simulations that resolve convective plumes, gravity currents, and the secondary instabilities leading to entrainment, we show that advective fluxes due to differential heating contribute to the evolution of the mixed layer in waterbodies with significant shallow areas. A heat balance is used to assess the relative importance of differential heating to the one‐dimensional effects of radiative heating and diffusive cooling at the ice‐water interface in lakes of varying morphologies.

We report novel results of a numerical experiment designed for examining the basin-scale hydrodynamics that control the mass, momentum, and energy distribution in a daily wind-forced, small thermally-stratified basin. For this purpose, the 3-D Boussinesq equations of motion were numerically solved using large-eddy simulation (LES) in a simplified (trapezoidal) stratified basin to compute the flow driven by a periodic wind shear stress working at the free surface along the principal axis. The domain and flow parameters of the LES experiment were chosen based on the conditions observed during summer in Lake Alpnach, Switzerland. We examine the diurnal circulation once the flow becomes quasi-periodic. First, the LES results show good agreement with available observations of internal seiching, boundary layer currents, vertical distribution of kinetic energy dissipation and effective diffusivity. Second, we investigated the wind-driven baroclinic cross-shore exchange. Results reveal that a near-resonant regime, arising from the coupling of the periodic wind-forcing (T = 24 h) and the V2H1 basin-scale internal seiche (T_V2H1 ≈ 24 h), leads to an active cross-shore circulation that can fully renew near-bottom waters at diurnal scale. Finally, we estimated the bulk mixing efficiency, Γ, of relevant zones, finding high spatial variability both for the turbulence intensity and the rate of mixing (10^–3 ≤ Γ ≤ 10^–1). In particular, significant temporal variability along the slopes of the basin was controlled by the periodic along-slope currents resulting from the V2H1 internal seiche.

Nonstationary convection forced by distributed buoyancy sources is a fundamental environmental fluid mechanics process, particularly in ice-covered freshwater waterbodies. In this paper, we present novel field-based results that characterise the diurnal evolution of the main energetics of radiatively-driven convection in ice-covered lakes that is the radiatively-induced buoyancy flux, B, and the kinetic energy dissipation rate, ε. To estimate the spatiotemporal distribution of ε, we applied scale similarity of the velocity structure functions to identify the fine turbulence scales from high-frequency velocity measurements. The field study was carried out at Lake Vendyurskoe, Russia, in April 2016. Small-scale velocity fluctuations were measured using acoustic Doppler current profiler in a 2 m layer beneath the ice cover. The method was proven to be valid for low-energy convection without mean shear. The inertial subrange, covering order of magnitude in the spatial domain, was identified by fitting the ²/₃ scaling power law to the structure function method, thus confirming the regime of fully developed turbulence. The calculated rate of dissipation of turbulent kinetic energy ε reaches values up to 3 × 10^–9 m² s^–3. Although a strong correlation between ε and B was observed, ε picks up about 1 h later after the onset of the heating-phase. This delay roughly corresponds to the turnover time of the energy containing eddies. We finally observed a decay of ε at night, during the relaxation-phase, but, interestingly, the level remained above the statistical error.

We derive the mechanical energy budget for shallow, ice‐covered lakes energized by penetrative solar radiation. Radiation increases the available and background components of the potential energy at different rates. Available potential energy drives under‐ice motion, including diurnally active turbulence in a near‐surface convective mixing layer. Heat loss at the ice‐water interface depletes background potential energy at a rate that depends on the available potential energy dynamics. Expressions for relative energy transfer rates show that the pathway for solar energy is sensitive to the convective mixing layer temperature through the nonlinear equation of state. Finally, we show that measurements of light penetration, temperature profiles resolving the diffusive boundary layer, and an estimate of the kinetic energy dissipation rate can be combined to estimate the forcing rate, the rate of heat loss to the ice, and efficiencies of the energy pathways for radiatively driven flows.

Basin scale seiches in lakes are important elements of the total energy budget and are a driver of fluxes of important ecological parameters, such as oxygen, nutrients, and sediments. At present, the extraction of the damping ratios of surface seiches, which are directly related to the capacity of seiches to drive these fluxes through the increased mixing of the water column, is reliant on spectral analysis which may be heavily influenced by the transformation of water level records from the time domain to the frequency domain, and which are sensitive to the level of noise present within the data. Existing spectral-based methods struggle to extract the periods of surface seiches which are of similar magnitude due to the overlap between their spectral responses. In this study, the principles of operational modal analysis, through the random decrement technique (RDT), currently used primarily in the analysis of high rise structures and in the aeronautical industry and not previously applied within the fields of limnology or ecology, are applied to barotropic seiches through the analysis of water level data for Lake Geneva, Switzerland, and Lake Tahoe, USA. Using this technique, the autocorrelation of the measurements is estimated using the RDT and modal analysis can then be carried out on this time-domain signal to estimate periods of the dominant surface seiches and the corresponding damping ratios. The estimated periods show good agreement with experimental results obtained through conventional spectral techniques and consistent damping ratios are obtained for the dominant surface seiche of Lake Tahoe. The effect of input parameters is discussed, using data for the two lakes, alongside discussion of the application of RDT to the study of internal seiches and current barriers to its application. RDT has great potential for the analysis of both surface and internal seiches, offering a method through which accurate damping ratios of seiche oscillations may be obtained using readily available data without necessitating spectral analysis.

Wind driven oscillations of water surfaces (surface seiches) are a key driver of ecological fluxes, such as of oxygen and nutrients, within lakes and reservoirs. These seiches exhibit low frequency oscillations with periods typically in the range from minutes to days, closely-spaced modes and a low signal to noise ratio, which makes modal analysis by simple spectral techniques challenging. Time domain operational modal analysis, specifically the Random Decrement Technique, was applied to extract modal parameters for the dominant surface seiches of Lake Geneva, Switzerland, and Lake Tahoe, USA. The results obtained show good agreement with existing literature on the periods and damping ratios of the seiches, historically obtained in this field by visual selection of free-decay sections of the time series or half-power peak picking. These analyses illustrate the potential for the use of techniques developed for operational modal analysis in areas outside the realm of mechanical and civil structures dynamics.

Freshwater lakes and impoundments are an important source of methane into the atmosphere. Methane is the second most important greenhouse gas after carbon dioxide and accounts for about 75% of the climate impact of lakes. The only process that reduces methane emissions from lakes is the oxidation by methanotrophic bacteria. This diverse group of bacteria has the unique ability to use methane and oxygen to gain energy and produce biomass.
Aim of this dissertation was to improve our understanding of the ecology of these bacteria and their importance to prevent methane emissions from freshwater lakes. From the outset, the interdisciplinary project was designed to complement field studies with modelling. The focus was in particular on the methane dynamics during the autumnal overturn of stratified lakes. During this critical phase, large amounts of methane that accumulate in the anoxic hypolimnion during stratification are transported to the surface. [...]

Phytoplankton growth depends on various factors, and primarily on nutrient availability, light and water temperature, whose distributions are largely controlled by hydrodynamics. Our main objective is to analyse the link between spatial and temporal variability of surface water temperature and algal concentration in a large lake by means of remote sensing and hydrodynamic modelling. We compare ten years of satellite images showing chlorophyll concentrations and surface water temperature of Lake Geneva. Our observations suggest different correlations depending on the season. Elevated chlorophyll concentrations in spring are correlated with warmer zones. But, in summer, higher chlorophyll concentrations are observed in colder zones. We show with a three-dimensional hydrodynamic model that the spatial variability of the surface water temperature reflects the upwelling and downwelling zones resulting from wind forcing. In springtime, nearshore downwellings induce locally increased surface temperature and stratification, which are associated with high chlorophyll concentration. In summertime, colder surface temperature area, often interpreted as transient upwellings, represents the thermal surface signature of wind-induced basin-scale internal waves, bringing either nutrients or phytoplankton from deeper layers to the surface. Our study suggests the latter to be the dominant process, with the basin-scale internal wave activity and associated transient summertime upwellings and downwellings having little net effects on the algal concentration. This study finally demonstrates the necessity to connect remote sensing retrievals and three-dimensional hydrodynamic modelling to properly understand the dynamic of the lake ecosystems.

Aquatic physics in inland water is a crucial subject for studying aquatic ecosystems.Transport and mixing are of tremendous importance for the pace at which chemical and biological processes develop. Recent observations allow to distinguish mixing and transport processes in stratified lakes and reservoirs. The surface and bottom boundary layer are turbulent while the lake interior remains comparatively quiescent.

We examine near-shore and pelagic current variability in Lake Geneva, a large and deep lake in western Europe, using observations from several measurement locations and a three-dimensional numerical model for the period 2014–2016. Linear internal seiche modes excited by wind forcing clearly appear as peaks in the energy spectra for measurements in off-shore locations. In contrast, spectra from the near-shore data, where currents interact with the lake bed, reveal a negligible contribution of internal seiches to the total kinetic energy. A similar contrast is seen in the spectra obtained from the numerical model at the same locations. Comparing the contribution of the different terms in the vertically-averaged momentum equation from the modeling results shows that the nonlinear advective term dominates in the near-shore boundary layer. Its contribution decays with distance from shore. The width of this near-shore boundary layer, which may extend for several kilometers, seems to be mainly determined by local topography. Both field measurements and modeling results indicate that nonlinear dynamics are of primary importance in the near-shore boundary layer.

Human impacts on biodiversity are well recognized, but uncertainties remain regarding patterns of diversity change at different spatial and temporal scales. Changes in microbial assemblages are, in particular, not well understood, partly due to the lack of community composition data over relevant scales of space and time. Here, we investigate biodiversity patterns in cyanobacterial assemblages over one century of eutrophication and climate change by sequencing DNA preserved in the sediments of ten European peri-Alpine lakes. We found species losses and gains at the lake scale, while species richness increased at the regional scale over approximately the past 100 years. Our data show a clear signal for beta diversity loss, with the composition and phylogenetic structure of assemblages becoming more similar across sites in the most recent decades, as have the general environmental conditions in and around the lakes. We attribute patterns of change in community composition to raised temperatures affecting the strength of the thermal stratification and, as a consequence, nutrient fluctuations, which favoured cyanobacterial taxa able to regulate buoyancy. Our results reinforce previous reports of human-induced homogenization of natural communities and reveal how potentially toxic and bloom-forming cyanobacteria have widened their geographic distribution in the European temperate region.

We investigate the impact on remote sensing reflectance by the vertical non-uniformities of water constituents. Reflectance simulated for 210 pairs of in situ measured chlorophyll-a and turbidity profiles (z = 0–20 m) from Lake Geneva are compared to simulations for uniform constituent gradients and non-uniform profiles approximated by Gaussian curves, orthogonal layers and steady gradients. Relevant concentration ranges are between 0 and 17 mg m⁻³ for chlorophyll-a and 0 and 4.6 g m⁻³ for total suspended matter within the photic layer. Our results show that mesotrophic lakes are specifically sensitive to non-uniformities with 20% of the 210 samples used in this study showing deviations of the spectral angle > 5° between a uniform assumption and observations which mostly occur for deeper-laying water constituents. By stressing the different use of blue and red parts of the spectrum, we argue further that algorithms are affected by variable vertical structures of algal and inorganic particles. Finally, we demonstrate that approximation models of the vertical structure of water constituents are a good solution to better account for non-uniformities in the development of invertible bio-optical models.

Extreme weather events may be just as important as gradual trends for the long-term trajectories of ecosystems. For alpine lakes, which are exposed to both exacerbated atmospheric warming and intense episodic weather events, future conditions might not be appropriately forecast by only climate change trends, i.e. warming, if extreme events have the potential to deflect their thermal and metabolic states from their seasonal ranges. We used high‐frequency monitoring data over three open‐water seasons with a one‐dimensional hydrodynamic model of the high‐altitude Lake Muzelle (France) to show that rainstorms or windstorms, notwithstanding their intensity, did not trigger long‐lasting consequences to the lake characteristics when light penetration into the lake was not modified. In contrast, storms associated with high turbidity input from the watershed ("turbid storms") strongly modified the lacustrine hydrodynamics and metabolism for the rest of the open‐water season through reduced light penetration. The long‐lasting effects of turbid storms were related to the inputs and in‐lake persistence of very light glacial suspensoids from the watershed. The occurrence of the observed turbid storms was not related to the wind or rain intensities during the events. Instead, the turbid storms occurred after dry and atypically warm spells, i.e. meteorological conditions expected to be more frequent in this alpine region in the upcoming decades. Consequently, storm events, notwithstanding their intensity, are expected to strongly imprint the future ecological status of alpine lakes under climate warming.

Thermal responses of inland waters to climate change varies on global and regional scales. The extent of warming is determined by system-specific characteristics such as fluvial input. Here we examine the impact of ongoing climate change on two alpine tributaries, the Aare River and the Rhône River, and their respective downstream peri-alpine lakes: Lake Biel and Lake Geneva. We propagate regional atmospheric temperature effects into river discharge projections. These, together with anthropogenic heat sources, are in turn incorporated into simple and efficient deterministic models that predict future water temperatures, river-borne suspended sediment concentration (SSC), lake stratification and river intrusion depth/volume in the lakes. Climate-induced shifts in river discharge regimes, including seasonal flow variations, act as positive and negative feedbacks in influencing river water temperature and SSC. Differences in temperature and heating regimes between rivers and lakes in turn result in large seasonal shifts in warming of downstream lakes. The extent of this repressive effect on warming is controlled by the lakes hydraulic residence time. Previous studies suggest that climate change will diminish deep-water oxygen renewal in lakes. We find that climate-related seasonal variations in river temperatures and SSC shift deep penetrating river intrusions from summer towards winter. Thus potentially counteracting the otherwise negative effects associated with climate change on deep-water oxygen content. Our findings provide a template for evaluating the response of similar hydrologic systems to on-going climate change.

Low oxygen concentrations in lakes and reservoirs are an ongoing environmental concern, particularly in light of increasing anthropogenic activity and climate change. Oxygen depletion processes in lakes are still not completely understood and a variety of models have been proposed based on limited field observations. Here, we present field measurements of oxygen depletion processes in a deep lake, Lake Geneva (Switzerland). The aim of this study was to quantify three basic processes controlling hypolimnetic oxygen depletion and their relative contribution to the total oxygen depletion (TOD) rate. Sediment oxygen uptake (SOU) and the flux of reduced substances were estimated based on oxygen microprofile measurements and sediment core data of reduced substances. Acoustic Doppler current profiler measurements and hydrodynamic modeling were used to ensure that SOU was measured under typical hydrodynamic conditions. Comparison with long-term monitoring data allowed for an estimate of the relative importance of SOU and water column mineralization (WCM). Results show a decrease in both SOU and WCM down to mid-depth which could not be explained by changes in hydrodynamic conditions or temperature. Below mid-depth, TOD increased due to an enhanced sediment area to water volume ratio (α). This vertical pattern of oxygen depletion is driven by (1) lake morphometry paired with increasing α, and (2) decreasing organic matter mineralization in the water column with depth. The findings are explained by a model which separates the oxygen depletion into an exponentially decreasing component, representing the fast-decaying fraction of the organic matter, and a constant background component.

Lake biological parameters show important spatio-temporal heterogeneities. This is why explaining the spatial patchiness of phytoplankton abundance has been a recurrent ecological issue and is an essential prerequisite for objectively assessing, protecting and restoring freshwater ecosystems. The drivers of these heterogeneities can be identified by modeling their dynamics. This approach is useful for theoretical and applied limnology. In this study, a 3D hydrodynamic model of Lake Geneva (France/Switzerland) was created. It is based on the Delft3D suite software and includes the main tributary (Rhône River) and two-dimensional high-resolution meteorological forcing. It provides 3D maps of water temperature and current velocities with a 1 h time step on a 1 km horizontal grid size and with a vertical resolution of 1 m near the surface to 7 m at the bottom of the lake. The dynamics and the drivers of phytoplankton heterogeneities were assessed by combining the outputs of the model and chlorophyll-a concentration (Chl-a) data from MERIS satellite images between 2008 and 2012. Results highlight physical mechanisms responsible for the occurrence of seasonal hot-spots in phytoplankton abundance in the lake. At the beginning of spring, Chl-a heterogeneities are usually caused by an earlier onset of phytoplankton growth in the shallowest and more sheltered areas; spatial differences in the timing of phytoplankton growth can be explained by spatial variability in thermal stratification dynamics. In summer, transient and locally higher phytoplankton abundances are observed in relation to the impact of basin-scale upwelling.

Ice-covered waterbodies are far from being quiescent systems. In this paper, we investigate ice-covered freshwater basins heated by solar radiation that penetrates across waters with temperatures below or near the temperature of maximum density. In this scenario, solar radiation sets a radiative buoyancy flux, Φ_r, that forces increments of temperature/density in the upper fluid volume, which can become gravitationally unstable and drive convection. The goal of this study is twofold. We first focus on formulating the mechanical energy budget, putting emphasis on the conversion of Φ_r to available potential energy, E_a. We find that E_a results from a competition among Φ_r and the irreversible mixing controlled by the diapycnal and the laminar mixing rates, respectively. Secondly, and based on the above result, we introduce an integral formulation of the mixing efficiency to quantify the rate of mixing over the relevant time scale τ , ɳ_c ≡ ∆E_{b, τ}/E_{r, τ}, where ∆E_{b, τ} and E_{r, τ} are the change of background potential energy and the time-integrated Φ_r over τ. The above definition is applied to estimate ɳ_c for the first time, finding an approximate value of ɳ_c ≈ 0.65 . This result suggests that radiatively heated ice-covered waterbodies might be subject to high mixing rates. Overall, the present work provides a framework to examine energetics and mixing in ice-covered waters.

The Swiss Federal Office for the Environment (FOEN) is currently evaluating options to add a nation-wide lake temperature monitoring program to their monitoring activities, which already includes a network of river temperature monitoring. The objective of this short report is to provide a preliminary assessment about the following questions from the perspective of a research institute:
• Why do we need a lake temperature monitoring?
• What should be considered in selecting monitoring locations?
• What should be considered in designing a monitoring system?
This report should thus form the basis for a more detailed assessment that will have to be conducted during the development of a strategy for the lake monitoring.

This paper presents an improvement of a one-dimensional lake hydrodynamic model (Simstrat) to characterize the vertical thermal structure of deep lakes. Using physically based arguments, we refine the transfer of wind energy to basin-scale internal waves (BSIWs). We consider the properties of the basin, the characteristics of the wind time series and the stability of the water column to filter and thereby optimize the magnitude of wind energy transferred to BSIWs. We show that this filtering procedure can significantly improve the accuracy of modelled temperatures, especially in the deep water of lakes such as Lake Geneva, for which the root mean square error between observed and simulated temperatures was reduced by up to 40 %. The modification, tested on four different lakes, increases model accuracy and contributes to a significantly better reproduction of seasonal deep convective mixing, a fundamental parameter for biogeochemical processes such as oxygen depletion. It also improves modelling over long time series for the purpose of climate change studies.

Turbulent mixing produced by breaking of internal waves plays an important role in setting the patterns of downwelling and upwelling of deep dense waters and thereby helps sustain the global deep ocean overturning circulation. A key parameter used to characterize turbulent mixing is its efficiency, defined here as the fraction of the energy available to turbulence that is invested in mixing. Efficiency is conventionally approximated by a constant value near one sixth. Here we show that efficiency varies significantly in the abyssal ocean and can be as large as approximately one third in density stratified regions near topographic features. Our results indicate that variations in efficiency exert a first-order control over the rate of overturning of the lower branch of the meridional overturning circulation.

Numerical simulations were carried out to investigate gyres within open lacustrine embayments subjected to parallel-to-shore currents. In such embayments, gyre formation occurs due to flow separation at the embayment’s upstream edge. High momentum fluid from the mixing layer between the embayment and offshore flows into the embayment and produces recirculating flow. Systematic numerical experiments using different synthetic embayment configurations were used to examine the impact of embayment geometry. Geometries included embayments with different aspect ratios, depths and embayment corner angles. The magnitudes of the recirculation and turbulent kinetic energy (TKE) in the embayment vary significantly for angles in the range 40°–55°. Embayments with corner angles less than 50° have much stronger recirculation and TKE, other parameters remaining the same. The numerical findings are consistent with gyre formation observed in two embayments located in Lake Geneva, Switzerland, and thus help explain flow patterns recorded in lacustrine shoreline regions.

Anthropogenic heat emissions into inland waters influence water temperature and affect stratification, heat and nutrient fluxes, deep water renewal, and biota. Given the increased thermal stress on these systems by growing cooling demands of riparian/coastal infrastructures in combination with climate warming, the question arises on how to best monitor and manage these systems. In this study, we investigate local and system-wide physical effects on the medium-sized perialpine Lake Biel (Switzerland), influenced by point-source cooling water emission from an upstream nuclear power plant (heat emission ∼700 MW, ∼18 W m⁻² lake wide). We use one-dimensional (SIMSTRAT) and three-dimensional (Delft3D-Flow) hydrodynamic numerical simulations and provide model resolution guidelines for future studies of thermal pollution. The effects on Lake Biel by the emitted excess heat are summarized as: (i) clear seasonal trend in temperature increase, locally up to 3.4°C and system-wide volume mean ∼0.3°C, which corresponds to one decade of regional surface water climate warming; (ii) the majority of supplied thermal pollution (∼60%) leaves this short residence time (∼58 days) system via the main outlet, whereas the remaining heat exits to the atmosphere; (iii) increased length of stratified period due to the stabilizing effects of additional heat; (iv) system-wide effects such as warmer temperature, prolonged stratified period, and river-caused epilimnion flushing are resolved by both models whereas local raised temperature and river short circuiting was only identifiable with the three-dimensional model approach. This model-based method provides an ideal tool to assess man-made impacts on lakes and their downstream outflows.

Der Energie Service Biel plant das Seewasserwerk Ipsach und die die dazugehörige Seewasserfassung im Bielersee zu ersetzen. Neben den üblichen Parametern der Wasserqualität (wie Temperatur, Sauerstoff, Trübung, Geruch, Geschmack) spielt im Bielersee auch die Hangstabilität der Sedimentablagerungen eine zentrale Rolle für die Wahl der Position der Seewasserfassung.

Dissolved oxygen microprofiles at the sediment-water interface of Lake Geneva were measured concurrently with velocities 0.25 to 2 m above the sediment. The measurements and scaling analyses indicate dissolved oxygen fluctuations and turbulent fluxes in exceedance of molecular diffusion in the proximity of the sediment-water interface. The measurements allowed the parameterization of the turbulent diffusion as a function of the dimensionless height above the sediment and the turbulence above the sediment-water interface. Turbulent diffusion depended strongly on the friction velocity and differed from formulations reported in the literature that are based on concepts of turbulent and developed wall-bounded flows. The dissolved oxygen microprofiles and proposed parameterization of turbulent diffusion enable a foundation for the similarity scaling of oxygen microprofiles in proximity to the sediment. The proposed scaling allows the estimation of diffusive boundary layer thickness, oxygen flux, and oxygen microprofile distribution in the near-sediment boundary layer.

Lake surface water temperatures (LSWTs) are sensitive to atmospheric warming and have previously been shown to respond to regional changes in the climate. Using a combination of in situ and simulated surface temperatures from 20 Central European lakes, with data spanning between 50 and ∼100 years, we investigate the long-term increase in annually averaged LSWT. We demonstrate that Central European lakes are warming most in spring and experience a seasonal variation in LSWT trends. We calculate significant LSWT warming during the past few decades and illustrate, using a sequential t test analysis of regime shifts, a substantial increase in annually averaged LSWT during the late 1980s, in response to an abrupt shift in the climate. Surface air temperature measurements from 122 meteorological stations situated throughout Central Europe demonstrate similar increases at this time. Climatic modification of LSWT has numerous consequences for water quality and lake ecosystems. Quantifying the response of LSWT increase to large-scale and abrupt climatic shifts is essential to understand how lakes will respond in the future.

Early-spring under-ice convection in the Petrozavodsk Bay of Lake Onega (Russia) was investigated as part of an interdisciplinary research project conducted during March 2015. Measurements performed using a thermistor chain and vertical profiling sensors were used to examine temperature dynamics in the convectively mixed and stratified layers of the lake. Radiative transfer through the ice was high leading to a large convective mixed layer (up to 20 m deep) during daytime. Convective velocity was evaluated using two different methods. It is shown that convective velocity (a maximum value of ~7.4 mm s⁻¹, and daytime average of 3.9 mm s⁻¹) is completely damped during the restratifying night hours. We observed internal waves in the thermocline below the convective mixed layer with intriguing variations between night and day. Maximum of internal wave energy was found to start in the afternoon and continue long after the end of solar radiation forcing. Our analysis indicates that local convective processes are key forcing mechanisms for the generation of internal waves in ice-covered lakes. We also hypothesize that spatial differential heating between the nearshore regions and the centre of the bay (e.g. density current intruding the thermocline) could be a source of internal waves in ice-covered lakes.

Low concentrations of dissolved oxygen remain a global concern regarding the ecological health of lakes and reservoirs. In addition to high nutrient loads, climate-induced changes in lake stratification and mixing represent additional anthropogenic menace resulting in decreased deepwater oxygen levels. The analysis of 43 years of monitoring data from Lake Geneva shows no decreasing trend neither in the areal hypolimnetic mineralization rate nor in the extent of hypoxia. Instead, hypoxic conditions are predominantly controlled by deep mixing in winter and much less by the trophic variations over the past decades. To reproduce winter mixing, the one-dimensional hydrodynamic model SIMSTRAT was specially adapted to deep lakes and run for several climate scenarios. The simulations predicted a decrease in the maximum winter mixing depth from an average of ∼172 m for 1981–2012 to ∼136 m and ∼127 m in response to predicted atmospheric temperatures between 2045–2076 and 2070–2101 according to Intergovernmental Panel on Climate Change scenarios. Concurrently, events with complete homogenization of temperature and oxygen in winter will decrease by ∼50%. Consequently, the hypolimnetic oxygen concentrations will significantly decrease. These results demonstrate that changes in deep mixing can have stronger impact than eutrophication on the deepwater oxygen levels of oligomictic lakes.

Directives and legislations worldwide aim at representatively and continuously monitoring the ecological status of surface waters. In many countries, chlorophyll-α concentrations (CHL) are used as an indicator of phytoplankton abundance and the trophic level of lakes or reservoirs. In-situ measurements of water quality parameters, however, are time-consuming, costly and of unknown but naturally limited spatial representativeness. In addition, the variety of the involved lab and field measurement methods and instruments complicates comparability and reproducibility.
Taking Lake Geneva as an example, 1234 satellite images from the MERIS sensor on the Envisat satellite from 2002 to 2012 are used to quantify the spatial and temporal variations of CHL concentrations. Based on histograms of spring, summer and autumn CHL estimates, the spatial representativeness of two existing in-situ measurement locations is analysed. Appropriate sampling frequencies to capture CHL peaks are examined by means of statistical resampling. The approaches proposed allow determining optimal in-situ sampling locations and frequencies. Their generic nature allows for adaptation to other lakes, especially to establish new survey programmes where no previous records are available.

Le laboratoire de Physique des Systèmes Aquatiques (APHYS-EPFL) a été mandaté par le Domaine Immobilier et Infrastructures (DII-EPFL) afin de mieux comprendre l’évolution de la température de l’eau ainsi que ces fluctuations observées à la station de pompage et de proposer des solutions pour l’implantation d’une nouvelle pompe d’eau du lac pour le campus.