Articles | Volume 26, issue 2
Research article
24 Jan 2022
Research article |  | 24 Jan 2022

Seasonality of density currents induced by differential cooling

Tomy Doda, Cintia L. Ramón, Hugo N. Ulloa, Alfred Wüest, and Damien Bouffard

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.

1 Introduction

As lateral boundaries in lakes, littoral regions exchange matter and energy with the surrounding watershed. The littoral region receives both sediment and chemicals from the terrestrial environment, including anthropogenic contaminants (Rao and Schwab, 2007). Particulate matter accumulates in the littoral zone (Cyr, 2017), and this biologically active region can host unique and intensified biogeochemical reactions (Wetzel, 2001). Cross-shore flows, which transport water laterally, connect a lake's littoral region to its pelagic region, thereby facilitating critical ecosystem processes (MacIntyre and Melack, 1995). They control the residence time of nearshore waters and redistribute heat as well as dissolved and particulate compounds throughout the lake. Representing lakes as isolated vertical water columns fails to account for the effects of cross-shore flows on lake biogeochemistry (Effler et al., 2010; Hofmann, 2013).

Figure 1Data-based schematic of the cooling-driven thermal siphon. The schematic shows the plateau, littoral, and mixed regions; the seasonality of the forcing; and the variables used for the transport scaling. Here, the littoral zone is the region upslope of the MT mooring, where the current velocity is measured and transport variables are calculated. The cross-shore temperature field is linearly interpolated from a transect of conductivity–temperature–depth (CTD) profiles collected on the morning of 22 August 2019 (08:20–08:50 UTC), from x=225 m (white dashed line) to x=714 m. Black dots on the x axis show the location of the profiles. The green dashed line in the seasonality diagram corresponds to the transition period between the mixing period (winter) and the stratified period (summer), when the mixed layer is not well defined.


Various mechanisms drive horizontal exchange flows. The pioneer limnologists have already elucidated the effect of large-scale surface and internal gravity waves on horizontal transport (Forel, 1895; Wedderburn, 1907; Mortimer, 1952). In addition to serving as the primary energy source for gravity waves, wind is also responsible for fueling lateral flows via wind-driven cross-shore circulation, cross-shore Ekman transport in rotating environments, and differential deepening (Imberger and Patterson, 1989; Brink, 2016). Buoyancy-driven lateral flows can be induced by horizontal temperature gradients due to differential heating or cooling. Such gradients could form when heat fluxes are spatially heterogeneous, due to large-scale differences in meteorological forcing (Verburg et al., 2011), variable geothermal heating (Roget et al., 1993), shading from vegetation (Lövstedt and Bengtsson, 2008), wind sheltering (Schlatter et al., 1997), or variation in turbidity (MacIntyre et al., 2002). Horizontal temperature gradients also occur when waters with variable depth experience a spatially uniform heat flux (Horsch and Stefan, 1988; Monismith et al., 1990). For a given surface area, the volume of water in the shallow littoral region is less than its offshore counterpart. Thus, littoral waters will cool or warm faster than pelagic waters. Monismith et al. (1990) designated cross-shore circulation resulting from a bathymetrically induced temperature gradient as a thermal siphon. The fact that all lakes have sloping boundaries makes thermal siphons potentially ubiquitous.

This study focuses on thermal siphons driven by differential cooling, wherein the lake temperature exceeds the temperature of maximum density (Fig. 1). The study does not consider heating-driven thermal siphons because induced lateral exchange remains weaker and more sensitive to wind disturbance than cooling-driven thermal siphons (Monismith et al., 1990; James et al., 1994). Atmospheric forcing affects lakes at diurnal, seasonal, or synoptic scales. This includes periods of net cooling when the net surface heat flux is directed towards the atmosphere (Bouffard and Wüest, 2019). This diel or seasonal surface cooling mixes and deepens the surface layer by penetrative convection. Once the shallow littoral region becomes fully mixed, it cools down faster than the deep pelagic region. Littoral waters become denser and plunge as a downslope density current, which intrudes at the base of the mixed layer during stratified conditions. At the surface, a reverse flow brings water from the pelagic to the littoral region (Horsch and Stefan, 1988; Monismith et al., 1990). We will hereafter refer to the density current induced by differential cooling simply as a thermal siphon (TS).

TSs have often been evoked to explain cold anomalies observed in lakes; these include tilted isotherms above the sloping sides (Talling, 1963; Eccles, 1974; MacIntyre and Melack, 1995; Woodward et al., 2017) or cold water intrusions (Wells and Sherman, 2001; Peeters et al., 2003; Rueda et al., 2007; Forrest et al., 2008; Ambrosetti et al., 2010). However, only a few studies have provided detailed in situ descriptions of the process. The first extensive studies of TSs assessed power plant heat disposal within cooling lakes (Adams and Wells, 1984). Warm water discharged into these reservoirs flows onshore at the surface, cools, and induces a TS in sidearms. However, these observations are not representative of naturally formed TSs. Cooling lakes constantly receive heat from power plants, which increases the intensity and duration of surface cooling compared with that experienced in other lakes and reservoirs (Horsch and Stefan, 1988). An extensive review of the literature suggests that Monismith et al. (1990) provided the first in situ observations of a diurnal cycle associated with a naturally formed TS. The authors studied a sidearm of Wellington Reservoir (Australia) in summer and described the complete diurnal cycle of daytime differential heating and nighttime differential cooling. In both cases, they observed cross-shore flows at velocities of ∼2 cm s−1. Cross-shore circulation was characterized by inertia with a delay of several hours between the shift in forcing conditions and the flow reversal. Other studies have reported the presence of TSs along the thalweg of narrow reservoir embayments (James and Barko, 1991a, b; James et al., 1994; Rogowski et al., 2019) but also along the sloping sides of lakes (Sturman et al., 1999; Thorpe et al., 1999; Fer et al., 2002a, b; Pálmarsson and Schladow, 2008) and between different basins (Roget et al., 1993). In most of these examples, TS characteristics comported with the description given in Monismith et al. (1990). These included thermal stratification above the sloping bottom, inertial effects, and a cross-shore velocity on the order of 1 cm s−1. The work by Fer et al. (2002b) on the sloping sides of Lake Geneva (Switzerland) represents one of the most extensive TS studies. The authors collected continuous measurements over two winters and identified density currents initiated in the middle of the night and lasting until the late morning. The research detected increasing discharge with distance downslope and found that one TS event could flush the littoral region almost twice.

The lake stratification and the intensity and duration of the daily cooling period vary seasonally (Bouffard and Wüest, 2019). These changes in forcing conditions may affect the occurrence and magnitude of TSs over time. Previous studies have focused on specific TS events and did not monitor the process over different seasons. A 3-year-long dataset on Lake Anna (VA, USA; Adams and Wells, 1984) addressed a cooling lake and, thus, did not reflect the seasonality of naturally formed TSs. Although the time series collected by Roget et al. (1993) between two lobes of Lake Banyoles (Spain) spans 8 months (October–May), these authors do not address seasonal variability in the observed TS. Geothermal fluxes between the lobes enhance differential cooling in Lake Banyoles, which may modify TS seasonality. James and Barko (1991b) performed six tracer experiments between June and September. These authors obtained estimates of TS velocity that varied between the experiments and were higher for days with stronger lateral temperature gradients. While the experiments did not identify a clear seasonal trend, thermistor arrays detected an increase in the occurrence of differential cooling from May to August, suggesting a higher frequency of TS events in late summer. This hypothesis was not verified due to the lack of long-term velocity measurements. Thus, the seasonality of TS occurrence and intensity has not been systematically investigated in lakes.

Dimensional analysis and laboratory experiments have provided relationships linking forcing and background conditions (surface cooling, stratification) to cross-shore transport by TSs (Harashima and Watanabe, 1986), but these lack validation in the field. The comprehension of the TS seasonality would improve the prediction of TS events in lakes and their contribution to exchanges between littoral and pelagic regions. The present research monitored TSs in a small temperate lake over a 1-year period. The study establishes three characteristics of TS seasonality including its occurrence (Sect. 3.3), intensity (Sect. 3.4), and diurnal dynamics (Sect. 3.5). The observed TS seasonality clearly relates to forcing parameters and, thus, pinpoints key parameters for predicting TS events in lakes (Sect. 4). The observation and interpretation of TS events across seasonal timescales provides a comprehensive understanding of the conditions required to form TSs in lakes.

2 Methods

2.1 Study site and field measurements

Rotsee is a small peri-alpine monomictic lake near Lucerne, Switzerland (Fig. 2a). It is located at 419 m a.s.l. (above sea level) and has a surface area of 0.5 km2 with a mean and maximum depth of 9 and 16 m, respectively. The main in- and outflows are located at the southwestern and northeastern ends of the lake, respectively, and have a low discharge of ∼0.1 m3 s−1. Rotsee is an elongated lake (2.5 km long, 0.2 km wide), with steep sides (slope of 15) across its longitudinal axis and more gradual slopes (1.5) at its two ends. The northeastern end of the lake finishes with a 200 m long plateau region of approximately 1 m depth that can trigger TSs. Moreover, Rotsee is located in a depression and sheltered from wind; hence, it is famous for international rowing championships and an ideal field-scale laboratory where convective processes are distinctively observable.

Figure 2(a) Bathymetry of Rotsee indicating the location of the two moorings MB and MT (red crosses), the cross-shore transect of CTD profiles (white dots), the weather station (green star), and the length of the plateau and littoral regions (Lp and Llit, respectively). Blue arrows correspond to the main inflow and outflow. The location of Rotsee is shown on the map of Switzerland using a black dot (source: Federal Office of Topography). (b) Schematic of the mooring MT. A similar setup is used for the thermistor array at MB. The detailed setup is provided in Appendix A.

This study focuses on quantifying TSs originating from the northeastern plateau region (Fig. 2a). Because of the elongated shape of Rotsee, we use the 2D (xz) framework shown in Fig. 1 by orienting the x axis along the thalweg. We assume that TSs flow preferentially along the x axis, and we neglect the influence of perpendicular flows. We will now refer to the northeastern end of the lake as the “shore” and call the direction of the x axis the “offshore direction”.

We monitored the background stratification at the deepest location (“background mooring” MB, approx. 16 m deep) as well as the dynamics of TSs off the plateau region (“TS mooring” MT, approx. 4 m deep) from March 2019 to March 2020 (Fig. 2a). Mooring MT was located along the thalweg, at the beginning of the sloping region. This shallow water column is already vertically mixed in summer by the action of surface cooling. We briefly describe the moorings below and provide more detailed information on the specifications and setup of the instruments in Appendix A. Each mooring included an array of thermistors (Fig. 2b). The evolution of the background temperature was monitored at MB with Vemco Minilog II-T loggers (accuracy: 0.1 C). High-resolution thermistors RBR TR-1050 and RBRduet T.D (accuracy: 0.002 C) were used at MT to track finer-scale temperature fluctuations near the sloping boundaries. The thermistors were installed from 1 to 13 m above the bottom with a 1 m vertical resolution at MB and from 0.25 to 2 m above the bottom with a 0.25 m resolution at MT. The evolution of the near-surface temperature (0.2 to 0.3 m below the surface) was monitored with Vemco loggers attached to surface buoys at both locations. An upward-looking acoustic Doppler current profiler (ADCP, Nortek Aquadopp profiler 1 MHz) was deployed at MT for fine-scale measurements of the bottom currents (Fig. 2b). It collected velocity profiles with a vertical bin resolution of 0.05 m from 0.25 to 3 m above the bottom. Velocity was obtained by averaging 512 burst measurements every 15 min. A week-long maintenance of the moorings was conducted every 40 d on average. The measurement periods of the ADCP were shorter in December and January due to low battery (Appendix A). The monitoring was stopped in May 2019 to avoid any hazard during the European Rowing Championships.

A WxPRO Campbell weather station was installed at the lake shore, near the plateau region, in September 2019 (Fig. 2a). It measured atmospheric pressure, air temperature, wind speed, incoming shortwave radiation, and relative humidity every 10 min. In addition, incoming longwave radiation was obtained from the Lucerne weather station (470212′′ N, 81804′′ E; source MeteoSwiss), ∼4 km from Rotsee. The Lucerne station also provided meteorological forcing when the reference Rotsee weather station was not measuring (data gaps and period before September 2019). In these cases, the Lucerne dataset was adjusted based on the Rotsee observations from September 2019 to June 2020. Pressure, air temperature, and shortwave radiation were linearly corrected (R2>0.9). This simple method did not allow reliable correction to the relative humidity and wind speed. Instead, we applied a neural network fitting approach to correct these two variables. The shallow neural network consisted of one hidden layer with 50 neurons and was trained by a Levenberg–Marquardt algorithm. A total of 70 % of the reference period from September 2019 to June 2020 was used to train the neural network, and the remaining percentage was equally distributed for validation and testing. All of the meteorological variables measured in Lucerne were used as inputs to the network, and the procedure was repeated for each of the two target variables. The network performance was satisfactory, with a coefficient of determination R2 and a root-mean-square error ERMS between observations and estimates of R2≈0.72 and ERMS≈0.67 m s−1 for wind speed and of R2≈0.94 and ERMS≈4.7 % for relative humidity. The uncertainty in wind speed and relative humidity estimates leads to an average uncertainty of 5 % and 3 % of the surface heat fluxes (Sect. 2.2), respectively.

We captured the spatial variability in TSs on the morning of 22 August 2019 (Fig. 1) by collecting 11 conductivity–temperature–depth (CTD) profiles (CTD60M, Sea & Sun Technology) along the x axis, between 2 and 14 m depth (Fig. 2a).

2.2 Heat and buoyancy fluxes

Heat fluxes were estimated from meteorological data and lake surface temperature at MB. We assumed that the meteorological conditions and the heat fluxes were spatially uniform over the lake surface (0.5 km2). Heat fluxes were defined positive in the upward direction (lake cooling) and negative in the downward direction (lake heating). The non-penetrative surface heat flux HQ0 is defined as

(1) H Q 0 = H C + H E + H LW , in + H LW , out W m - 2 ,

where HC and HE are the respective sensible and latent heat fluxes, and HLW,in and HLW,out are the respective incoming and outgoing longwave radiation. The net surface heat flux H0,net includes the shortwave radiative flux at the surface HSW,0:

(2) H 0 , net = H Q 0 + H SW , 0 W m - 2 .

The different terms of the heat budget were calculated similarly to Fink et al. (2014). Cloudiness modulates the proportion of direct and diffuse shortwave radiation and was derived from daily clear-sky irradiance as in Meyers and Dale (1983). We modified the empirical wind and air temperature-based calibration function, f, in the sensible and latent heat flux estimates of Fink et al. (2014) to take the lake fetch into account, as proposed in McJannet et al. (2012).

The non-penetrative surface buoyancy flux B0 was inferred from HQ0 under the assumption that only heat fluxes modify the potential energy near the surface (Bouffard and Wüest, 2019):

(3) B 0 = α g H Q 0 ρ C p , w W kg - 1 .

In Eq. (3), α=-(1/ρ)ρ/T is the thermal expansivity of water [C−1], g is the gravitational acceleration [m s−2], ρ is the water density [kg m−3], and Cp,w is the specific heat of water [J C−1 kg−1]. The thermal expansivity was estimated from the surface temperature using the relationship reported in Bouffard and Wüest (2019). To calculate water density, we used the equation of state from Chen and Millero (1986), with a measured average salinity of S=0.2 g kg−1 for Rotsee. Our analysis assumed a constant salinity over the year, as the seasonal changes in surface density are controlled by temperature fluctuations rather than variations in salinity.

Similarly, a radiative surface buoyancy flux BSW,0 was obtained as follows:

(4) B SW , 0 = α g H SW , 0 ρ C p , w W kg - 1 .

The net buoyancy flux at the surface is B0,net=B0+BSW,0. B0,net>0 indicates a destabilizing buoyancy flux (net cooling), whereas B0,net<0 indicates a stabilizing buoyancy flux (net heating). We used B0,net to identify the daily cooling and heating phases (Sect. 3.1), and we averaged B0 over the daily cooling phase to quantify the intensity of convection (Sect. 2.6).

2.3 Mixed-layer depth

We used hourly averaged temperature data from MB to calculate the thermocline depth (ht) and the mixed-layer depth (hML). Based on the review of Gray et al. (2020), we estimated hML using the temperature-gradient method with a threshold of dT/dz=0.05C m−1. The water column was defined as fully mixed when dT/dz<0.05C m−1 at all depths. In this case, hML was set to the depth of the bottom temperature sensor. For stratified conditions, the thermocline depth was defined as the depth of maximum gradient in dT/dz. Starting from ht and moving upward, the lower end of the mixed layer was reached when the local temperature gradient dropped to dT/dz<0.05C m−1. When dT/dz>0.05C m−1 at all depths, the entire column was stratified and hML was zero.

2.4 Occurrence of thermal siphons

We used the velocity data from the Aquadopp profiler at MT to estimate the cross-shore transport over the bottom 3 m of the sloping region. Quality control of the velocity was performed by discarding values with beam correlation below 50 % or signal amplitude weaker than 6 dB above noise floor. The horizontal velocity was projected onto the x axis (angle of 56 from north), which crosses the isobath at MT perpendicularly (Fig. 2a). Following the 2D framework of Fig. 1, we will now call Ux the “cross-shore velocity”.

The time series of Ux was then divided into 24 h subsets, starting and ending at 17:00 UTC. Each 24 h subset was analyzed separately to identify TS events. We decided to focus on periods when MT was located in the sloping mixed region, with TSs flowing at the bottom of the water column (Fig. 1). This condition allowed us to relate our current measurements to scaling formulae of downslope density currents (Sect. 3.4). A downslope TS event was detected on a specific day if the three following conditions were met:

  1. There was significant cross-shore flow, where the depth-averaged velocity of the current Ux>0.5 cm s−1 for at least 2 h;

  2. There was weak wind during the selected event, where |LMO,avg|/hML,avg<0.5 (LMO,avg and hML,avg are the averaged Monin–Obukhov length scale (Wüest and Lorke, 2003) and mixed-layer depth during the cross-shore flow event, respectively);

  3. There was a mixed water column at MT before the onset of the flow, where dT/dz<0.05C m−1 at all depths.

Cross-shore flows that did not meet conditions no. 2 and no. 3 were noted as “wind-driven circulation” and “stratified flows”, respectively. The term “stratified flows” refers here to any cross-shore flow that occurred when the mixed-layer depth was shallower than the lake depth at MT. This includes basin-scale internal waves and intruding TSs. Further justifications of the above criteria are provided in Appendix D.

2.5 Cross-shore transport by thermal siphons

The cross-shore transport was calculated for each identified TS event (Sect. 2.4). For each event, we defined subperiods with continuous positive Ux (further details in Appendix D). We calculated the unit-width volume of water transported over each of these subperiods:

(5) V x , i = t 0 , i t f , i q x ( t ) d t m 2 ,

where t0,i and tf,i are the initial and final times of subperiod i, respectively, and qx is the unit-width discharge defined as

(6) q x ( t ) = z bot ( t ) z top ( t ) U x ( z , t ) d z m 2 s - 1 .

In Eq. (6), zbot and ztop are the respective bottom and top interfaces of the region over which Ux>0 at time t. For each day, we then defined the TS flushing period as the subperiod with the largest volume transported Vx,i.

Four daily quantities characterizing the cross-shore transport were calculated over the flushing period of each day (t0ttf): the average cross-shore velocity Uavg and unit-width discharge qavg, the maximum cross-shore velocity Umax (both in time and depth), and the flushing timescale τF=Vlit/qavg, with Vlit≈500 m2 representing the unit-width volume of the littoral region upslope of MT (Fig. 1). Uavg and qavg were obtained as follows:


2.6 Scaling formulae

We relate the transport quantities introduced in Sect. 2.5 to the forcing parameters using scaling formulae from theoretical and laboratory studies.

The horizontal velocity scale of TSs is U(B0L)1/3, where B0 is the destabilizing surface buoyancy flux, and L is a horizontal length scale (Phillips, 1966). Following Wells and Sherman (2001), we used the horizontal length scale LML based on the mixed-layer depth:

(9) U = c U B 0 L ML 1 / 3 m s - 1 ,

where cU is a proportionality coefficient. For each value of hML, we calculated LML as the distance between the northeastern edge of the lake and the location where the mixed layer intersects the lake sloping bottom (Fig. 1).

The unit-width discharge is scaled by q(B0L)1/3h, where h is a vertical length scale. We used the average depth of the littoral region upslope of MT (hlit=Vlit/Llit, Fig. 1) as h:

(10) q = c q B 0 L ML 1 / 3 h lit m 2 s - 1 ,

where cq is a proportionality coefficient. We defined the flushing timescale from Eq. (10) as

(11) τ F = V lit q = c F L lit B 0 L ML 1 / 3 [ s ] ,

where cF=1/cq. The flushing timescale represents the time that a TS needs to flush the littoral region of volume Vlit.

The daily cooling phase of duration τc is defined by B0,net>0. After the onset of surface cooling, a TS does not form immediately. Vertical convective mixing takes place before the horizontal density gradient starts to drive a TS (Wells and Sherman, 2001; Ulloa et al., 2022). The period between the start of the cooling phase and the TS formation is quantified by an initiation timescale (τini). If the littoral region is initially mixed vertically, τini equals the transition timescale τt derived by Ulloa et al. (2022) under constant surface cooling:

(12) τ t = 2 L ML - L p 2 / 3 B 0 1 / 3 1 - h p / h ML 1 / 3 [ s ] ,

where Lp≈170 m and hp≈1 m, representing the length and depth of the plateau region of Rotsee, respectively (Fig. 1). This timescale is based on a three-way momentum balance between the lateral pressure gradient due to differential cooling and the inertial terms (Finnigan and Ivey, 1999).

In Eqs. (9)-(12), B0 and LML are averaged over τc every day. We discuss the choice of the length scales LML and hlit in more detail in Sect. 4.3. The key parameters used for the scaling formulae and TS dynamics are listed in Table 1.

Table 1Key parameters related to the forcing parameters, bathymetry, and TS dynamics. For seasonally varying parameters, the given ranges of values are daily averages.

Download Print Version | Download XLSX

3 Results

3.1 Diurnal cycle

The magnitude of the net surface buoyancy flux B0,net oscillates daily, and the intensity of the diel fluctuation is modulated at a seasonal timescale (Fig. 3a). We illustrate the diurnal cycle with an example on 9–10 September 2019 (Fig. 3b, c, d). This summertime period is characterized by a strong day–night variability, which is ideal to elucidate how changes in forcing conditions affect the formation and destruction of TSs. The net surface buoyancy flux B0,net varies from -2×10-7 W kg−1 during daytime to 1×10-7 W kg−1 during nighttime (Fig. 3b). In contrast to B0,net, B0 remains positive over 24 h, indicating a continuous cooling at the air–water interface. In late summer, B0 is 10-7 W kg−1 with a limited variability of ±5×10-8 W kg−1 between day and night. The radiative forcing BSW,0 is the main driver of the diurnal variability, with diel variations that are 1 order of magnitude larger than B0 during cloud-free summer days like 9–10 September 2019. We observe differential cooling at night, with lateral temperature gradients of 5×10-4C m−1 in late summer (Fig. 1).

Figure 3Diurnal cycle of a TS. (a) Net surface buoyancy flux as a function of time during the day, represented over 1 year (from 13 March 2019 to 13 March 2020). DOY stands for “day of the year”. A 30 d moving average has been applied to smooth the time series. The 45 d long period without measurement in May 2019 is shown in gray. The contour lines have a spacing of 1×10-7 W kg−1. The dashed line corresponds to the diurnal cycle shown in panels (b)(c), and (d) (9–10 September 2019; DOY 252–253). (b) Time series of surface buoyancy fluxes defining the cooling and heating phases on 9–10 September 2019. (c) Vertical velocity measured every 15 min at MT (depth of ∼4 m) as a function of time and height above the sediment. Positive (purple) values correspond to a flow moving upward. Strong vertical movements are the signature of convective plumes. (d) Cross-shore velocity measured every 15 min at MT as a function of time and height above the sediment. Positive (purple) values correspond to a flow moving offshore (southwestward flow). Black lines are 0.05 C spaced isotherms that are linearly interpolated between each thermistor. The TS detected by the algorithm and used for the transport calculation is depicted in blue. The vertical resolution of the thermistors is shown by the horizontal black ticks on the right axis. Horizontal arrows identify three periods: (1) deepening of the mixed layer, (2) transition period before the initiation of the TS, and (3) flushing period.


We divide the diurnal cycle into a net cooling phase during the night (B0>|BSW,0|) and a net heating phase during the day (B0<|BSW,0|). Starting the period at 17:00 UTC, the water column is initially thermally stratified by the radiative heating of the daytime period (Fig. 3d). During the first part of the night, convective cooling mixes and deepens the surface layer (period no. 1 in Fig. 3d). The water column becomes entirely mixed at MT around 00:00 UTC (see the vertical isotherms starting around 00:00 UTC in Fig. 3d). The convective plumes intensify in the second part of the night to reach vertical velocities of ∼5 mm s−1 (Fig. 3c). A total of 3 h after complete vertical mixing at MT, a cross-shore circulation is initiated (period no. 3 in Fig. 3d), with a downslope flow near the bottom (positive Ux) and an opposite flow above (negative Ux). The velocity magnitude of ∼1 cm s−1 is typical for TSs. We will hereafter refer to the period during which a TS flows as the flushing period. The interface between the two flows oscillates vertically by 1–2 m (∼100 % of the TS thickness), suggesting that the transport is limited by the vertical mixing of convective plumes. At the end of the cooling phase, the water column is distinctly stratified near the bottom, indicating the presence of persistent colder water flowing downslope. A striking observation is the intensification of the TS in phase with the weakening of vertical convection at the onset of radiative heating. The cross-shore velocity increases above 2 cm s−1, and the induced bottom stratification reaches ∼0.1C m−1 at the beginning of the heating phase. This strong flushing lasts for several hours, while radiative heating increases and re-stratifies the water above the density current. The cross-shore flow weakens in the late morning but continues until noon.

The example detailed above presents all of the characteristics of a typical TS, with different phases and inertia between the cooling phase and the flushing period. Other examples, provided in Appendix B (Fig. B2), show that the timing and magnitude of the flow change with the forcing conditions.

3.2 Seasonal variability in the forcing parameters

Two key parameters for TS formation are B0 and hML (Sect. 2.6), which are indirectly related to B0,net. Hereafter, we refer to these parameters as “forcing parameters” of TSs. The daily intensity of convective cooling is estimated based on the convective velocity scale w*=(B0hML)1/3 (Deardorff, 1970), where B0 and hML are daily averages during the cooling phase at MB. These two parameters are both seasonally dependent, and their relationship is characterized by an annual hysteresis cycle (Fig. 4a). Although the surface heat flux HQ0 increases from winter to summer, the strong seasonal variability in B0 (1 order of magnitude larger in summer than in winter) is mostly due to seasonal changes in the thermal expansivity α.

Figure 4Seasonal variability in the forcing parameters over 1 year (from March 2019 to March 2020). (a) Daily averages of mixed-layer depth as a function of surface buoyancy flux during cooling periods. Gray dashed lines indicate the corresponding convective velocity scale. The black dashed line with arrows is a qualitative representation of the annual cycle. (b) Monthly averages of convective velocity scale w* during the cooling phase (green line), duration of the cooling phase τc (blue line), and net surface buoyancy flux B0,net (linearly interpolated color map). Error bars represent the standard deviation of w* and τc.


From spring to summer, hML remains shallow (<2 m) and B0 increases from 10-9 W kg−1 in March to 10-7 W kg−1 in August. The latter results in an increase in the convective velocity scale from w*3.8±0.7 to w*6.7±0.9 mm s−1 over the same period (Fig. 4a, b). The lake undergoes net daily heating over this period (B0,net<0), with a cooling phase that remains shorter than 15 h (Fig. 4b). In August, B0,net0 and daily heating is balanced by daily cooling. At this period, the intensity of convective mixing reaches its maximum, and the diurnal heating–cooling cycle is pronounced (Fig. 3a).

From September to January, the lake undergoes a daily net cooling (B0,net>0) and the mixed layer deepens by convection at an average rate of ∼0.1 m d−1, leading to a complete mixing in mid-December. Over the same period, B0 continuously decreases due to the drop in α with colder temperatures. This decrease in B0 balances the deepening of hML, and w* remains constant around w*7±1 mm s−1 from August to November. The duration of the cooling phase increases in autumn to reach its maximum in November (τc21±2 h). Convective cooling is occurring almost continuously at that time. In winter, the convective velocity w* drops to w*3.5±1 mm s−1, and the effect of the cooling phase is again balanced by the heating phase (B0,net0, τc15±3 h in February). The differences between the cooling and heating phases of the diurnal cycle are reduced in winter, with B0,net0 over the entire day (Fig. 3a). Note that winter 2019–2020 was warm, and Rotsee remained mixed from December to February (no inverse stratification, surface temperature T>4C).

3.3 Seasonal occurrence of thermal siphons

We apply our detection algorithm for cross-shore flows (Sect. 2.4) on 227 d with continuous measurements at MT from March 2019 to February 2020. We identify 156 d with a significant cross-shore flow (69 % of the days with observations), over which 85 d (37 %) were reported as TS events. The remaining days with identified cross-shore flow were either associated with wind-driven circulation for 13 d (6 %) or stratified flows for 58 d (26 %). We estimate the percentage of occurrence of TSs relative to the number of days with measurements (pTS) for each month (Fig. 5). Before late July, the mixed layer is usually shallower than 4 m (Fig. 4a) and stratified flows prevail at MT (Fig. 5). The first events identified as TSs occur on 28 July and 31 July (pTS=12.5 % in July). However, TSs do not reach a significant percentage of occurrence (∼50 %) until August. The occurrence increases in autumn to reach 87 % in November. In winter, TSs are less common: they only occur on 38 % of the days in December, 2 d in early January (not shown in Fig. 5 due to the few days measured during this month), and are not observed later in winter.

Figure 5Monthly occurrence of TSs and other cross-shore (southwestward) flows at MT, expressed as a percentage of the days with measurements. The other cross-shore flows are divided between wind-driven circulation (|LMO|/hML>0.5) and stratified flows at MT (dT/dz>0.05C m−1). The red dashed line indicates the 15 d moving average of the percentage of occurrence of TSs. The number of days with measurements (n) is shown for each month. Months with less than 10 d of measurements have been removed (May 2019 and January 2020).


The 15 d moving average in Fig. 5 reveals the short-term variability due to synoptic changes in the meteorological conditions that naturally modulate the seasonal pattern. Although the monthly averaged occurrence increases from August to November, the biweekly percentage of occurrence drops over shorter periods. These periods are associated with strong heating that re-stratifies the water column at MT, reduces surface cooling at night, and prevents TS formation (Appendix B, Fig. B1). The 15 d averaged occurrence reaches 100 % in November, when a TS is present over 25 consecutive days.

Although TSs become more frequent in autumn, the percentage of occurrence of cross-shore flows (pCS) remains constant at around 80 %–90 % from July to November. Lateral flows observed at MT in spring and summer primarily result from direct wind forcing or subsequent stratified mixing. These flows may also include TSs intruding before reaching MT and, therefore, not counted as downslope TSs (further details in Appendix D). These different processes become less frequent in autumn due to the deeper mixed layer and are replaced by downslope TSs. Overall, there seems to be an effect of TSs on the number of cross-shore flows if we compare the period from July to December (pCS80±14 %) with the rest of the year (pCS47±12 %), suggesting that TSs are the main mechanism connecting the littoral and pelagic regions in Rotsee.

3.4 Scaling the cross-shore transport

We compare our field-based transport estimates for the TS events identified in Fig. 5 with the laboratory-based scaling formulae introduced in Sect. 2.6 (Fig. 6). The daily average and maximum cross-shore velocities are linearly correlated with the horizontal velocity scale from Eq. (9) (R2>0.5, pval,F<0.01; Fig. 6a, b). The daily average velocity Uavg0.9±0.2 cm s−1 is 3 times lower than the daily maximum velocity Umax2.6±0.7 cm s−1 that occurs in the morning (Fig. 3). The best linear fits (with 0 intercept) are given by Uavg=0.33(B0LML)1/3 and Umax=0.99(B0LML)1/3. The 95 % confidence interval of the proportionality coefficient is 0.32cU0.34 for Uavg and 0.95cU1.02 for Umax. Despite the natural variability, a seasonal trend is distinguishable, with a decrease in Uavg and Umax by a factor of 2 from July to December. This is consistent with the weakening of convective cooling observed in Fig. 4.

Figure 6(a) Daily averaged cross-shore velocity, (b) daily maximum cross-shore velocity, (c) daily averaged unit-width discharge, and (d) flushing timescale as a function of scaling formulae (Sect. 2.6). The equation of the linear regressions with 0 intercept, the coefficient of determination (R2), and the p value of an F test (pval,F) are indicated.


The unit-width discharge qavg and flushing timescale τF are also well predicted by the scaling formulae (Eqs. 10 and 11, respectively), despite a larger scatter than for the cross-shore velocity (R2≈0.2–0.3, pval,F<0.01; Fig. 6c, d). The unit-width discharge qavg0.015±0.004 m2 s−1 corresponds to an average thickness of hTS,avg=qavg/Uavg1.8±0.2 m and a flushing timescale of τF9.7±2.7 h. The best relationships (with 0 intercept) are given by qavg=0.34(B0LML)1/3hlit (R2=0.27) and τF=2.99Llit/(B0LML)1/3 (R2=0.18). The 95 % confidence interval of the proportionality coefficient is 0.33cq0.36 for qavg and 2.83cF3.15 for τF. The strong daily variability is explained by the fluctuating TS thickness during the convective period (Fig. 3c, d), which affects the calculation of the discharge (further details in Appendix D). As observed for the cross-shore velocity, the cross-shore transport weakens from summer to winter: qavg decreases and τF increases by a factor of 2. The flushing timescale is τF≈7 h in summer but reaches τF≈20 h in winter. However, due to longer TS events in winter than in summer, the volume flushed by a single TS event is independent of the season and remains larger than the volume of the littoral region (Appendix C).

3.5 Flushing period

The diurnal cycle described in Sect. 3.1 varies at the seasonal scale, as the forcing parameters change. To assess the effects on the flushing period (period no. 3 in Fig. 3d), we average the diurnal cycle monthly between August and December (Fig. 7). The cooling and heating phases are illustrated with the diurnal cycle of B0,net. The duration of the cooling phase (blue shading in Fig. 7) markedly increases from August to December, as already observed in Figs. 3a and 4b. The histograms of |B0,netmax| show that the magnitude of and variability in B0,net decrease over the same period (Fig. 7b, d, f, h, j).

Figure 7Seasonal variability in the diurnal cycle. (a, c, e, g, i) Monthly average of the diurnal cycle of the unit-width discharge qx and FG parameter. Only days with observed TSs between August and December have been averaged. The shaded white area corresponds to average discharge ± standard deviation. The monthly averaged net buoyancy flux represented as a color map is normalized by the maximum absolute value of each month. The vertical black dashed line indicates the time of maximal qx. The transition timescale τt defines the starting time of the flushing, represented by a pink vertical dashed line. (b, d, f, h, j) Histograms of the maximum absolute value of the net buoyancy flux |B0,netmax| for each month. The red vertical line corresponds to |B0,netmax| of the monthly averaged diurnal cycles.


These changes in the forcing parameters have a direct effect on the cross-shore transport. The unit-width discharge in Fig. 7 (white dotted line) is obtained by monthly averaging the 24 h time series qx(t). The discharge qx(t) is calculated from Eq. (6) during the flushing period t0ttf and qx(t)=0 for t<t0 or t>tf. In August and September (Fig. 7a, b), qx increases at night to reach its maximum at the beginning of the heating phase. The flow generally stops in the late morning, depicted by a drop in qx at around 10:00 UTC. This cycle corresponds to the example in Fig. 3. The daily peak in flushing is also observed in October and November, but it is more spread with a continuous increase in cross-shore transport at night, which can already be seen in the evening (Fig. 7e, g). Thus, the average flushing duration of TSs is longer in late autumn (see, for example, the 24 h long TS event in Appendix B, Fig. B2). In December, the peak in discharge is reduced and the cross-shore transport is nearly continuous over 24 h (Fig. 7i). The time of maximal qx (vertical black dashed line in Fig. 7) is delayed from August to December. We associate this maximal qx with the heating phase, as it always occurs ∼1–2 h after B0,net becomes negative (except in December, during which time the peak is less pronounced).

The diurnal cycle of the cross-shore transport is also described by the flow geometry parameter introduced by Ulloa et al. (2022) as the ratio between the root-mean-square (RMS) of the cross-shore velocity and the RMS of the vertical velocity:

(13) F G ( t ) = U x ( t ) 2 U z ( t ) 2 ,

where < … > denotes a depth average. FG provides information on the relative magnitude of TSs with respect to the surface convection. We expect FG to be close to unity when surface convection is the only process acting on the water column. A deviation from unity towards larger values of FG (i.e., FG≫1) implies that the flow is anisotropic, with a stronger cross-shore component, as a TS develops. A sharp peak in the monthly averaged FG is visible at the beginning of the heating phase, from August to November (Fig. 7). It matches the time of maximal qx and results from the combined effect of TS-induced flushing (increase in Ux) and convection weakening (decrease in Uz).

To give a theoretical estimation of the starting time of the flushing, we estimate the initiation timescale of TSs from the monthly averaged transition timescale τiniτt (Eq. 12). The assumption of constant surface cooling is reasonable for nighttime cooling in Rotsee (Fig. 3b). The starting time, depicted by a vertical pink dashed line in Fig. 7, corresponds to (tsunset+τt), where tsunset is the time of sunset when BSW,0 drops to zero. The increase in τt from τt≈3 h in late summer to τt≈20 h in winter leads to a later theoretical onset of TSs in winter, despite the earlier time of sunset. Based on this seasonal variability in τt, TSs are expected to start in the evening in summer but in the morning in winter. This theoretical starting time, however, can be improved by considering the initial stratification at the beginning of the cooling phase, as further discussed in Sect. 4.2 and 4.4.

4 Discussion

4.1 Seasonality of thermal siphons

Our 1-year-long measurements show that both the occurrence and intensity of TSs vary seasonally in lakes with shallow littoral zones comparable to Rotsee. TSs occur regularly from late summer to early winter, with a maximum frequency in November, and are absent the rest of the year (Sect. 3.3). While the frequency of TSs increases in autumn, the intensity of the net cross-shore transport decreases compared with the summer period (Sect. 3.4). However, the daily averaged volume flushed from the littoral region increases from summer to autumn (Appendix C). The diurnal cycle is well defined in summer and divided into a cooling phase at night and a heating phase during daytime (Sect. 3.1). TSs form during the second part of the night and last until late morning, reaching a maximal flushing at the beginning of the heating phase. In autumn, the flushing duration increases, and the time of maximal flushing is shifted later in the day (Sect. 3.5).

To explain the TS seasonality, we need to relate it to the seasonal variability in the forcing parameters. The duration of the cooling phase (τc) is a key parameter to understand the TS occurrence and the diurnal dynamics, while the intensity of surface cooling (B0) and the stratification (hML) parametrize the cross-shore transport. We discuss the effects of the forcing parameters on the TS seasonality in the following sections.

Figure 8Timescales determining the occurrence of TSs. (a) Schematic of the three periods of the cooling phase at MT (Fig. 3d) parameterized by the mixing timescale τmix, transition timescale τt, and cooling duration τc. The mixed-layer depth is expressed as the relative depth hML=hML-hMT, with respect to hMT=4 m. The mixed layer deepens during the first period (0<t<τmix), until the complete mixing of the water column. The transition period (τmix<t<τini) is dominated by convection. TSs occur during the flushing period (τini<t<τc). (b) Effects of the seasonality of τmix, τt, τc, and hML on the TS occurrence. Monthly averages are represented, with shaded areas (τc, τt, and τmix) and error bars (hML) indicating the monthly standard deviation. Note the log scale for the axis of timescales. The gray shaded period corresponds to optimal conditions for the TS occurrence.


4.2 Effects of the forcing parameters on the occurrence of thermal siphons

TS occurrence can be predicted by estimating the initiation timescale τini and comparing it with the duration of the cooling phase (τc). A TS will occur if the cooling phase is long enough, i.e., τc>τini. In initially stratified surface waters, the initiation period is divided into (1) a stratified period, during which the mixed layer deepens, and (2) a transition period, during which the littoral region is vertically mixed (Figs. 3d, 8a). We express the initiation timescale as τini=τmix+τt, where τmix is the time needed to vertically mix the littoral region, and τt is the transition timescale given by Eq. (12). To estimate τmix, we use the deepening rate of the mixed layer (Zilitinkevĭc, 1991) as follows:

(14) d h ML d t = ( 1 + 2 A ) B 0 h ML N 2 m s - 1 ,

where N2=-(g/ρ)(dρ/dz) [s−2] is the squared buoyancy frequency below the mixed layer, ρ is the depth-averaged water density, and A≈0.2 is an empirical coefficient. The model of Eq. (14) assumes that the mixed layer deepens by convection only, without any wind contribution. This assumption is valid for calm conditions, which prevail in Rotsee due to wind sheltering (Zimmermann et al., 2021). The average duration required for mixing the water column at MT, assuming a constant surface cooling B0 and an initial mixed-layer depth hML,ini, can be derived from Eq. (14) as follows:

(15) τ mix = h MT 2 - h ML , ini 2 N 2 2 B 0 ( 1 + 2 A ) [ s ] .

In Eq. (15), the buoyancy frequency squared is approximated as N2=-(g/ρ)[ρ(z=-hML,ini)-ρ(z=-hMT)]/[hMT-hML,ini]. The depth hMT corresponds to the maximum depth of the littoral region upslope of MT (Fig. 1), which has to be mixed to observe downslope TSs at MT.

From the seasonality of the forcing parameters (Fig. 4), we predict the optimal period for the TS occurrence based on τini (Fig. 8b). In summer, the shallow mixed layer (hML≈2 m) and strong nighttime surface cooling (B010-7 W kg−1) lead to a short transition timescale of τt2±0.5 h. However, the strong stratification near the surface limits the deepening of the mixed layer and causes τmix to be too large for the water column to mix at MT (orange period in Fig. 8b). The littoral region remains stratified for most of the nights (hML<hMT), and the occurrence of downslope TSs at MT is low (τini>τc on average). Downslope TSs start to be observed in late summer when the mixed layer reaches hMT before the end of the night. Starting in October, the monthly averaged mixed layer is deeper than hMT, and τmix≈0. At the same time, τt increases because of the deepening of the mixed layer and the decrease in B0. The transition timescale becomes longer than τc in winter (purple period in Fig. 8b), reaching its maximum of τt23±7 h in January. This prevents TSs from forming in winter, except for days with continuous surface cooling (i.e., τc>24 h). The transition timescale decreases in spring, but τmix increases simultaneously due to the re-stratification of the littoral region, which prevents TSs from occurring. Thus, the conditions to observe downslope TSs at MT are optimal between September and December, when τiniτc (gray shaded period in Fig. 8b). This period coincides well with our frequent observations of TSs in autumn. The limiting factor for the formation of TSs in spring and summer is τmix, whereas it is τt in winter.

During the optimal period of occurrence, the variability in forcing conditions (B0,net) between days can modify the occurrence condition τini<τc at a shorter timescale than the seasonal scale. Days with TSs are characterized by higher B0,net than days without TSs (Appendix B, Fig. B1). This results in the short-term variability in the biweekly averaged occurrence in Fig. 5. Other mechanisms affect the TS occurrence on a daily basis, such as wind events. Direct wind forcing either enhances or suppresses TSs. We identified several wind events from September to December that prevented any flow in the offshore direction (i.e., no significant flow detected by the algorithm) or, conversely, generated a strong cross-shore circulation (wind-driven circulation in Fig. 5).

Although the TS occurrence over short periods may vary from one year to another, we expect the seasonal trend previously described to be repeated every year in Rotsee. We also suggest that a similar trend should be common in other temperate lakes. Most of the TS events reported for other systems were observed in late summer or autumn, corresponding to the occurrence period for Rotsee (Monismith et al., 1990; James and Barko, 1991a, b; James et al., 1994; Rogowski et al., 2019). James and Barko (1991b) showed that the percentage of days with differential cooling in Eau Galle Reservoir (WI, USA) increased from 31 % in May to 95 % in August. However, these high values do not correspond to the percentage of occurrence of TSs, as differential cooling can happen without necessarily forming a cross-shore circulation. Pálmarsson and Schladow (2008) also found that differential cooling was rare in May and June in Clear Lake (CA, USA). Only 2 nights were cold enough to result in cooling-driven TSs. In winter, TS events are not frequent in Rotsee, but they have been reported for larger systems such as Lake Geneva in January (Thorpe et al., 1999; Fer et al., 2002a, b). A different bathymetry modifies τt (Eq. 12), which becomes shorter if hp/hML is small or the slope is high, leading to a more frequent TS occurrence in winter. In addition, convective cooling stops earlier in small lakes like Rotsee that reach the temperature of maximum density before larger and deeper lakes.

4.3 Effects of the forcing parameters on cross-shore transport

The seasonality of the cross-shore transport induced by TSs can be examined using the scaling for the cross-shore velocity Ux (Eq. 9) and the unit-width discharge qx (Eq. 10). The timescale τF is also a key parameter to assess the seasonality of the littoral region flushing (Eq. 11). We showed that these scaling formulae reproduce the observed seasonality in Rotsee (Fig. 6). The weakening of the transport from summer to winter is explained by the decrease in B0 by a factor of 10, which overcomes the increase in LML due to the deepening of the mixed layer (Fig. 4). However, τF remains shorter than the duration of the cooling phase (Fig. 8b), and the littoral region upslope of MT (Vlit≈500 m2) is entirely flushed by a single event (Appendix C). The use of the constant length scale hlit in Eq. (10) implies that the daily averaged TS thickness, hTS,avg=qavg/Uavg, should not show any seasonal variability. We observe the following in Rotsee: hTS,avg remains close to its yearly average of 1.8±0.2 m, without any seasonal trend.

Table 2Comparison of the transport scaling formula q=cqh(B0L)1/3 between different studies on sloping basins, where h is a vertical length scale, and L is a horizontal length scale. The different notations can be found in Fig. 1 and Table 1. The littoral region is based on the location of the discharge measurement, as in Fig. 1. Lforc is the length over which the destabilizing forcing B0 applies.

Download Print Version | Download XLSX

We further compare our scaling for Ux, qx, and τF with previous studies on buoyancy-driven flows. The velocity scale U(B0L)1/3 from Phillips (1966) is commonly used, with L defined as the horizontal length scale along which a lateral density gradient is set. This length scale varies with the experimental configuration and basin geometry (Table 2). In laboratory experiments with a localized buoyancy loss, L is the length along which the destabilizing buoyancy flux applies (Harashima and Watanabe, 1986; Sturman et al., 1996; Sturman and Ivey, 1998). For more realistic systems undergoing uniform surface cooling, L is the length of the littoral region, which can be defined as the plateau zone (Ulloa et al., 2022) or the region vertically mixed by convection (Wells and Sherman, 2001). In the first case, L=Lp is a constant, whereas L=LML varies with the stratification in the second case. We chose to use LML because it corresponds to the region affected by differential cooling. We expect that a deeper mixed layer will increase the temperature difference ΔT between the littoral and pelagic regions, resulting in a higher velocity Ux. In Eqs. (10) and (11), the transport is expressed as a function of the velocity scale and the size of the littoral region of length Llit and depth hlit (Fig. 1). We defined the littoral region based on the location of MT because we measured the discharge out of this region. This definition of hlit in Eq. (10) is consistent with the fact that qx varies spatially, increasing with distance from the shore (Fer et al., 2002b). The vertical length scale h can also be chosen as the depth of the plateau (hp) (Sturman and Ivey, 1998; Wells and Sherman, 2001; Fer et al., 2002b; Ulloa et al., 2022) to predict the discharge from the plateau region only. Despite these different choices of length scales, the coefficient cq≈0.34 is close to other estimates from the literature (Table 2). Harashima and Watanabe (1986) found that cq increases with the flux Reynolds number defined as Ref=B01/3hlit2/LML2/3ν, where ν1.5×10-6 m2 s−1 is the kinematic viscosity of water. The empirical relationship cq0.43-1.67/Ref1/2 was obtained for Ref>50. For the days where TSs are observed in Rotsee, Ref varies from 250 in summer to 50 in winter, with an average value of ∼140 which leads to cq≈0.29, close to our estimate of cq≈0.34.

4.4 Effects of the forcing parameters on the flushing period

The seasonality of the forcing parameters does not only affect the occurrence and magnitude of TSs but also the diurnal dynamics of the littoral flushing (Fig. 7). The initiation of TSs always occurs several hours after the beginning of the cooling phase. This delay is modulated by τini and is consistent with previous studies reporting the formation of a TS at night (Monismith et al., 1990; James et al., 1994; Fer et al., 2002b; Pálmarsson and Schladow, 2008) or in the morning (Rogowski et al., 2019; Sturman et al., 1999). τt is a good estimate of the initiation period when the water column is initially mixed (τmix=0). In summer and on warm autumn days, however, the initiation of TSs occurs later than predicted by τt, due to the initial stratification at the beginning of the cooling phase (τmix≠0). This leads to a delay in the increase in qx and FG (Fig. 7a, c). In autumn, TS events lasting for more than 24 h lead to monthly averaged qx>0 before the starting time predicted by τt (Fig. 7e, g). In winter, τt→24 h implies that the rare TS events are all characterized by a continuous flow (Fig. 7i), without any diurnal cycle of flushing.

When a diurnal cycle is present in summer and autumn, the flushing period ends several hours after the beginning of the heating phase (Figs. 3, 7). This inertia has already been reported by Monismith et al. (1990), who observed TSs until the mid-afternoon, despite the opposite pressure gradients due to differential heating. In addition to the delay between the end of the cooling and the end of the flow, we also showed that a peak in flushing occurred a few hours after the beginning of the heating phase (Fig. 3). The time of maximal flushing is shifted later in the day from summer to winter, as the heating phase starts later (Fig. 7). The maximal flushing seems related to the transition from the cooling phase to the heating phase. This finding suggests that convective plumes eroding the flow control the weak initial transport during the cooling phase. A weakening of convection reduces this vertical mixing and finally enables the flushing to reach its maximum intensity. Further investigations are required to better understand this process.

4.5 Practical recommendations to predict and measure thermal siphons in other lakes

Rotsee is an ideal field-scale laboratory to investigate TSs, due to wind sheltering and its elongated shape minimizing complex recirculation. To assess the TS seasonality in other lakes, long-term velocity and temperature measurements in the sloping region are required. The developed algorithm to detect TS events (Sect. 2.4) and calculate the cross-shore transport (Sect. 2.5) fulfilled our requirements but remained lake-specific. The limitations of the algorithm are discussed in more detail in Appendix D. The general structure of the algorithm (Sect. 2.4) can serve as a basis for detecting TSs in lakes, by adapting the lake-specific criteria to other systems. The 2D framework of TSs requires specific validation in more complex nearshore systems and large lakes, where the topography, large-scale circulation, and Coriolis may also affect the TS dynamics (Fer et al., 2002b). In these systems, the alongshore velocity component of TSs must be considered in the cross-shore transport analysis. Further development is needed to build a more robust algorithm with lake-independent physically based criteria.

The effects of a different bathymetry can be predicted from the scaling discussed in the previous sections. A shallower nearshore plateau region or a steeper sloping region would decrease the transition timescale τt (Eq. 12), causing the initiation of TSs to happen earlier and more often. Higher slopes would also decrease the length of the mixed region LML and reduce the horizontal velocity of TSs (Eq. 9). Past observations of TSs have reported horizontal velocities ranging from ∼0.1 cm s−1 (James and Barko, 1991a, b) to ∼10 cm s−1 (Roget et al., 1993; Fer et al., 2002b), with Ux∼1 cm s−1 in most cases (Monismith et al., 1990; Sturman et al., 1999; Pálmarsson and Schladow, 2008; Rogowski et al., 2019). The cross-shore transport qx and flushing timescale τF are strongly dependent on the size and depth of the littoral region considered (Eqs. 10, 11). A deeper littoral region (larger hlit) would lead to a stronger discharge, and a longer littoral region (larger Llit) would take more time to flush. Values of TS thickness and discharge can be 1 order of magnitude larger in lakes deeper than Rotsee (Thorpe et al., 1999; Fer et al., 2002a; Rogowski et al., 2019). The effects of more complex bathymetries, departing from our 2D framework, could be further investigated with 3D numerical simulations.

Information on lake bathymetry and thermal structure are needed to predict the occurrence (τini<τc) and intensity of TSs (Ux, qx, τF) under specific forcing. This approach does not consider wind effects, which deserve further investigation. Windier conditions would hinder TSs by locally enhancing and reducing the cross-shore transport. In this case, the scaling formulae (Eqs. 9–11) should be modified to take wind shear into account. The littoral region would also be mixed faster under windy conditions, and the effects of wind mixing on τmix should be included in Eq. (15). These questions could be addressed by performing numerical simulations with varying wind speed.

The TS seasonality may evolve in a changing climate, which also needs to be investigated. Changes in heat fluxes, summer stratification, and surface temperature would affect both the intensity (B0) and occurrence (τini) of TSs.

5 Conclusions

The flushing of the littoral region by cross-shore flows increases the exchange between nearshore and pelagic waters. In this study, we investigated one of the processes that enhance the renewal of littoral waters, the so-called thermal siphon (TS) driven by differential cooling. From a 1-year-long monitoring of TSs in a small, temperate, wind-sheltered lake, we quantified the seasonality of the cross-shore transport induced by TSs. This seasonality is related to the intensity of surface cooling (surface buoyancy flux B0) and the lake stratification (mixed-layer depth hML). Three aspects of the TS seasonality are highlighted. First, TSs are a recurring process from late summer to winter (when τc>τini), occurring on ∼80 % of autumn days. Second, the seasonal changes in the TS-induced transport are well reproduced by the scaling Uavg0.33(B0LML)1/3, qavg0.34(B0LML)1/3hlit, and τF=Vlit/qavg2.99Llit/(B0LML)1/3, where LML is the length of the region mixed by convection, and Vlit=Llithlit is the unit-width volume of the littoral region. This study provides a field validation of this laboratory and theoretically based scaling. Third, the diurnal dynamics of the flushing follow the seasonal changes in the cooling and heating phases, by evolving from a well-defined diurnal cycle in summer to a more continuous flow in winter.

Our results demonstrate that TSs significantly contribute to the flushing of the nearshore waters. This process occurs frequently during the cooling season (>50 % of the time), each time flushing the entire littoral region. We stress that this buoyancy-driven transport is, perhaps counterintuitively, stronger in summer and in the morning. Such a timing has implications for the transport of dissolved compounds, with, for instance, maximal exchange between littoral and pelagic waters at a time of high primary production (summer and daytime). Overall, this study provides a solid framework to integrate the role of TSs in the lake ecosystem dynamics.

Appendix A: Instrumentation

Table A1Specifications and setup of the sensors from the two moorings and the meteorological station. For the ADCP, only the complete days of measurements used for the transport estimates are taken into account in the measurement periods.

Download Print Version | Download XLSX

Appendix B: Diurnal variability in the occurrence and duration of thermal siphons

In addition to the seasonality, TSs also vary between days. A change in forcing conditions (B0,net) between days may enhance or prevent TS formation, by modifying the occurrence condition τini<τc at a short timescale. We average the net buoyancy flux over 48 h (B0,net2d), for each 24 h subset starting at 17:00 UTC. B0,net2d is the average of B0,net between 00:00 UTC on the first day of the 24 h subset and 23:59 UTC on the following day. This quantity includes the effects of the forcing on (1) the initial stratification, which affects τini (day preceding the TS event), and (2) the duration of the cooling phase τc (following day). We compare the values of B0,net2d between days with and without TSs (Fig. B1). Days with TSs are characterized by higher B0,net2d than days without TSs, which is confirmed by a t test at a 5 % significance level, from August to December. The variability in the forcing conditions between days also affects the flushing duration of TSs (Fig. B2).

Figure B1Box plots of the 2 d averaged net surface buoyancy flux for each month, depending on the occurrence of TS events.


Figure B2Three examples of TS events, with (a–c) buoyancy fluxes, (d–f) vertical velocity, and (g–i) cross-shore velocity, as in Fig. 3. (a, d, g) A long TS event on 21–22 August 2019 due to short heating phases on both days. (b, e, h) A short TS event on 26–27 October 2019 due to strong heating on both days (re-stratification). (c, f, i) Continuous flushing on 20–21 November 2019 due to continuous net cooling. See the caption of Fig. 3 for more details about each panel. The gray shaded area in panels (a), (d), and (g) indicates the time of the CTD transect shown in Fig. 1.


Appendix C: Seasonality of the flushing of the littoral region

The unit-width volume flushed by each TS event is calculated by integrating the discharge qx over the flushing period (Eq. 5). We estimate the daily averaged flushed volume Vflush,day and the averaged flushed volume of TS events Vflush,TS by dividing the total volume of water flushed every month by the number of days with measurements and by the number of days with TSs, respectively. The volume Vflush,day includes the effect of the occurrence of TS events on the daily flushing, whereas Vflush,TS only depends on the intensity and duration of a TS event. We express both flushed volumes as fractions of the littoral region of volume Vlit (Fig. C1). The seasonal trend of Vflush,day is similar to the TS occurrence (Fig. 5), with an increase in flushing from July to November (Fig. C1a). This indicates that the increase in the TS occurrence from summer to autumn overcomes the weakening of the transport over the same period. The flushed volume of TS events does not vary significantly with season (Fig. C1b). A single TS event flushes the littoral region more than once on average. Overall, the TS occurrence is the primary factor controlling the seasonality of the littoral flushing in Rotsee.

Figure C1Fraction of the littoral region flushed by TSs in 1 d for each month. Months with less than 10 d of measurements have been removed as in Fig. 5 (May 2019 and January 2020). (a) The daily averaged flushed fraction including all days with measurements. (b) The daily averaged flushed fraction for days with TSs only.


Appendix D: Identification of thermal siphons by the algorithm

The developed algorithm used to detect TS events (Sect. 2.4) and calculate the cross-shore transport (Sect. 2.5) aimed at automatizing the identification of TSs and, thus, at limiting subjective bias in the characterization of the process. We manually assessed the performance of the algorithm over different days. We further tested the validity of the identified TS events with the fact that a cross-shore transport resulting from a TS should be associated with a decrease in water temperature (Fig. 1). We calculated the correlation r(Ux,T) between the cross-shore velocity Ux and the temperature T from the thermistor array at MT, at each depth and during the cross-shore flow events. The negative correlation during 93 % of the identified TS events (r(Ux,T)<-0.48 according to the 95 % confidence interval from a t test) and the lack of clear negative correlation (at a 5 % significance level) during the other days justified the skills of the detection algorithm. The general structure of the algorithm can be used in other systems. However, several criteria are lake-specific and must be adapted to the system of interest. We discuss the limitations of the algorithm below and provide suggestions for improvement.

The definition of a significant cross-shore flow by the algorithm (i.e., Ux>0.5 cm s−1 for at least 2 h) is lake-specific. This criterion is based on typical examples of TSs in Rotsee, where Ux is above 0.5 cm s−1 over most of the flushing period. Theoretical estimates of the period between the onset time and the end of the cooling phase (e.g., τc-τini10 h in autumn) show that most TS events last more than 2 h and justify the chosen threshold for the duration. From a sensitivity analysis (Fig. D1), we notice that the TS occurrence does not change significantly if the duration threshold is increased by a few hours. Shorter cross-shore flow events must be discarded because they are not consistent with the gravitational adjustment triggered by differential cooling; they could be wind-driven or result from free surface convection. In Rotsee, cross-shore flows driven by internal waves are expected to be short and to be discarded by the algorithm (TV1H1/22.5 h in summer, where TV1H1 is the period of V1H1 internal waves). Other types of data could be included to detect TSs, such as lateral temperature gradients, near-surface velocity (return flow), or vertical velocity (convective plumes). Using other techniques, like machine learning algorithms, could also be a useful approach to better identify TS events.

Figure D1Sensitivity analysis for the criteria used to define cross-shore flows. The effects of modifying the depth-averaged cross-shore velocity Ux and the averaging duration Δt are shown as changes in the percentage of occurrence of TSs (pTS) for each month. ΔpTS is the percentage difference with respect to the occurrence shown in Fig. 5, obtained with Ux=0.5 cm s−1 and Δt=2 h (reference point shown by the red cross).


The filter used to discard days with wind-driven circulation is based on the criterion |LMO|/hML=κ-1|u*/w*|3>0.5, where κ=0.41 is the von Kármán constant, and u* the friction velocity. The ratio u*/w* expresses the relative importance of wind shear compared with convection. The threshold value of 0.5 was successfully tested on wind-driven circulation events observed in Rotsee in winter. This filter discards the strongest wind events, but cross-shore flows with |LMO|/hML<0.5 can still be affected by a wind peak during the flushing period (e.g., interaction between wind and a TS). Being more conservative and decreasing the threshold value from 0.5 to 0.1 would discard TS events mainly in December (when w* is low) but would not significantly modify the TS occurrence during the other periods (Fig. D2). Additional filters could be implemented to distinguish between TSs and wind-driven cross-shore flows, based, for example, on high-resolution vertical velocity measurements, observed oscillations of the thermocline (e.g., wavelet analysis), estimates of the period of internal waves, and identification of upwelling events (e.g., Wedderburn and Lake numbers) (Imberger and Patterson, 1989).

Figure D2Filtering of the 227 d analyzed, shown as a time series of LMO/hML from March 2019 to March 2020. Days without a significant cross-shore flow or without a mixed layer are depicted as bottom ticks. Stratified flows are cross-shore flows without complete mixing at MT. The wind-driven circulation threshold LMO/hML=0.5 is shown by the horizontal dashed line.


By discarding days with stratified conditions at MT, the second filter does not consider TSs that intrude before reaching MT and appear as interflows. Thus, TS occurrence in summer is underestimated. In Fig. 5, the higher percentage of occurrence of stratified flows in July (81 %) compared with spring (around 40 %–50 %) suggests a possible contribution of intruding TSs. From visual inspection of each event, we estimated that intruding TSs could represent almost half of the stratified flows in July. Including downslope TSs only was relevant, as the transport properties of interflows can be different from downslope gravity currents, for which the scaling formulae were derived. In addition, intruding TS events could be confounded with other baroclinic flows, and further filtering steps would have been required to detect them correctly.

The transport quantities of TSs are averaged over the flushing period (Sect. 2.5). The onset of the flushing is challenging to define because of the interaction between TSs and convective plumes (Fig. 3c). The cross-shore velocity field can be very variable over time, with vertical fluctuations of the region with positive Ux (Fig. 3d). We decided to include only the period over which the region with positive Ux remained at the same depths (i.e., the limits of the flushing period correspond either to vertical displacements of the region of positive Ux or to Ux≤0 at all depths). This approach might sometimes include only a part of the TS event in the transport calculations. Daily averaged velocities should not be affected. The average thickness and discharge, however, vary between days depending whether the vertical oscillations of the TS interface are included or not. The latter is a possible reason for the larger variability observed for qavg and τF (Fig. 6).

Code and data availability

The data (raw and processed data as well as figures) and the scripts used for all of the analyses are available for download from the Eawag Research Data Institutional Collection (, Doda et al., 2021).

Author contributions

TD and DB designed the field experiments. TD led the data collection and analysis with help from DB, HNU, and CLR. TD wrote the initial draft of the manuscript, and all co-authors commented on and edited the text.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


We would like to sincerely thank our technician, Michael Plüss, for organizing the field campaigns and helping to set up and maintain the different instruments. We are grateful to the Canton of Luzern, the municipalities of Luzern (Lucerne) and Ebikon, the Rowing Centre Lucerne-Rotsee, the Quartierverein Maihof and Pro Natura associations, and the Rotsee-Badi for their support of our measurements. We are also indebted to Bieito Fernández-Castro for his help with estimating the heat fluxes. Moreover, we acknowledge Love Råman Vinnå, Edgar Hédouin, Josquin Dami, and Alois Zwyssig for their assistance in the field. Discussions with Mathew Wells, Martin Schmid, Oscar Sepúlveda Steiner, and Love Råman Vinnå helped to improve the data analysis and the quality of the paper. The meteorological data from the Lucerne weather station were provided by MeteoSwiss, the Swiss Federal Office of Meteorology and Climatology. We gratefully acknowledge the three anonymous reviewers for their helpful comments and suggestions.

Financial support

This study was financed by the Swiss National Science Foundation (“Buoyancy driven nearshore transport in lakes” project; HYPOlimnetic THErmal SIphonS, HYPOTHESIS, grant no. 175919).

Review statement

This paper was edited by Matthew Hipsey and reviewed by three anonymous referees.


Adams, E. E. and Wells, S. A.: Field measurements on side arms of Lake Anna, Va., J. Hydraul. Eng., 110, 773–793,, 1984. 

Ambrosetti, W., Barbanti, L., and Carrara, E. A.: Mechanisms of hypolimnion erosion in a deep lake (Lago Maggiore, N. Italy), J. Limnol., 69, 3–14,, 2010. 

Bouffard, D. and Wüest, A.: Convection in lakes, Annu. Rev. Fluid Mech., 51, 189–215,, 2019. 

Brink, K. H.: Cross-shelf exchange, Annu. Rev. Mar. Sci., 8, 59–78,, 2016. 

Chen, C.-T. A. and Millero, F. J.: Precise thermodynamic properties for natural waters covering only the limnological range, Limnol. Oceanogr., 31, 657–662,, 1986. 

Cyr, H.: Winds and the distribution of nearshore phytoplankton in a stratified lake, Water Res., 122, 114–127,, 2017. 

Deardorff, J. W.: Convective velocity and temperature scales for the unstable planetary boundary layer and for Rayleigh convection, J. Atmos. Sci., 27, 1211–1213,<1211:CVATSF>2.0.CO;2, 1970. 

Doda, T., Ramón, C. L., Ulloa, H. N., Wüest, A., and Bouffard, D.: Data for: Seasonality of density currents induced by differential cooling, Eawag Data Institutional Collection [data set],, 2021. 

Eccles, D. H.: An outline of the physical limnology of Lake Malawi (Lake Nyasa), Limnol. Oceanogr., 19, 730–742,, 1974. 

Effler, S. W., Prestigiacomo, A. R., Matthews, D. A., Gelda, R. K., Peng, F., Cowen, E. A., and Schweitzer, S. A.: Tripton, trophic state metrics, and near-shore versus pelagic zone responses to external loads in Cayuga Lake, New York, U.S.A., Fundam. Appl. Limnol., 178, 1–15,, 2010. 

Fer, I., Lemmin, U., and Thorpe, S. A.: Contribution of entrainment and vertical plumes to the winter cascading of cold shelf waters in a deep lake, Limnol. Oceanogr., 47, 576–580,, 2002a. 

Fer, I., Lemmin, U., and Thorpe, S. A.: Winter cascading of cold water in Lake Geneva, J. Geophys. Res., 107, 3060,, 2002b. 

Fink, G., Schmid, M., Wahl, B., Wolf, T., and Wüest, A.: Heat flux modifications related to climate-induced warming of large European lakes, Water Resour. Res., 50, 2072–2085,, 2014. 

Finnigan, T. D. and Ivey, G. N.: Submaximal exchange between a convectively forced basin and a large reservoir, J. Fluid Mech., 378, 357–378,, 1999. 

Forel, F. A.: Le Léman – Monographie Limnologique. Tome II, Editions Rouge, Lausanne, Switzerland,, 1895. 

Forrest, A. L., Laval, B. E., Pieters, R., and Lim, D. S. S.: Convectively driven transport in temperate lakes, Limnol. Oceanogr., 53, 2321–2332,, 2008. 

Gray, E., Mackay, E. B., Elliott, J. A., Folkard, A. M., and Jones, I. D.: Wide-spread inconsistency in estimation of lake mixed depth impacts interpretation of limnological processes, Water Res., 168, 115136,, 2020. 

Harashima, A. and Watanabe, M.: Laboratory experiments on the steady gravitational circulation excited by cooling of the water surface, J. Geophys. Res., 91, 13056–13064,, 1986. 

Hofmann, H.: Spatiotemporal distribution patterns of dissolved methane in lakes: how accurate are the current estimations of the diffusive flux path?, Geophys. Res. Lett., 40, 2779–2784,, 2013. 

Horsch, G. M. and Stefan, H. G.: Convective circulation in littoral water due to surface cooling, Limnol. Oceanogr., 33, 1068–1083,, 1988. 

Imberger, J. and Patterson, J. C.: Physical Limnology, edited by J. W. Hutchinson and T. Y. Wu, Adv. Appl. Mech., 27, 303–475,, 1989. 

James, W. F. and Barko, J. W.: Estimation of phosphorus exchange between littoral and pelagic zones during nighttime convective circulation, Limnol. Oceanogr., 36, 179–187,, 1991a. 

James, W. F. and Barko, J. W.: Littoral-pelagic phosphorus dynamics during nighttime convective circulation, Limnol. Oceanogr., 36, 949–960,, 1991b. 

James, W. F., Barko, J. W., and Eakin, H. L.: Convective water exchanges during differential cooling and heating: implications for dissolved constituent transport, Hydrobiologia, 294, 167–176,, 1994. 

Lövstedt, C. B. and Bengtsson, L.: Density-driven current between reed belts and open water in a shallow lake, Water Resour. Res., 44, W10413,, 2008. 

MacIntyre, S. and Melack, J. M.: Vertical and horizontal transport in lakes: linking littoral, benthic, and pelagic habitats, J. N. Am. Benthol. Soc., 14, 599–615,, 1995. 

MacIntyre, S., Romero, J. R., and Kling, G. W.: Spatial-temporal variability in surface layer deepening and lateral advection in an embayment of Lake Victoria, East Africa, Limnol. Oceanogr., 47, 656–671,, 2002. 

McJannet, D. L., Webster, I. T., and Cook, F. J.: An area-dependent wind function for estimating open water evaporation using land-based meteorological data, Environ. Model. Softw., 31, 76–83,, 2012. 

Meyers, T. and Dale, R.: Predicting daily insolation with hourly cloud height and coverage, J. Clim. Appl. Meteorol., 22, 537–545,<0537:PDIWHC>2.0.CO;2, 1983. 

Monismith, S. G., Imberger, J., and Morison, M. L.: Convective motions in the sidearm of a small reservoir, Limnol. Oceanogr., 35, 1676–1702,, 1990. 

Mortimer, C. H.: Water movements in lakes during summer stratification; evidence from the distribution of temperature in Windermere, Philos. T. Roy. Soc. Lond. B, 236, 355–398,, 1952. 

Pálmarsson, S. Ó. and Schladow, S. G.: Exchange flow in a shallow lake embayment, Ecol. Appl., 18, A89–A106,, 2008. 

Peeters, F., Finger, D., Hofer, M., Brennwald, M., Livingstone, D. M., and Kipfer, R.: Deep-water renewal in Lake Issyk-Kul driven by d cooling, Limnol. Oceanogr., 48, 1419–1431,, 2003. 

Phillips, O. M.: On turbulent convection currents and the circulation of the Red Sea, Deep-Sea Res. Oceanogr. Abstr., 13, 1149–1160,, 1966. 

Rao, Y. R. and Schwab, D. J.: Transport and mixing between the coastal and offshore waters in the Great Lakes: a review, J. Great Lakes Res., 33, 202–218,[202:TAMBTC]2.0.CO;2, 2007. 

Roget, E., Colomer, J., Casamitjana, X., and Llebot, J. E.: Bottom currents induced by baroclinic forcing in Lake Banyoles (Spain), Aquat. Sci., 55, 206–227,, 1993. 

Rogowski, P., Merrifield, S., Ding, L., Terrill, E., and Gesiriech, G.: Robotic mapping of mixing and dispersion of augmented surface water in a drought frequent reservoir, Limnol. Oceanogr. Meth., 17, 475–489,, 2019. 

Rueda, F., Moreno-Ostos, E., and Cruz-Pizarro, L.: Spatial and temporal scales of transport during the cooling phase of the ice-free period in a small high-mountain lake, Aquat. Sci., 69, 115–128,, 2007. 

Schlatter, J. W., Wüest, A., and Imboden, D. M.: Hypolimnetic density currents traced by sulphur hexafluoride (SF6), Aquat. Sci., 59, 225–242,, 1997. 

Sturman, J. J. and Ivey, G. N.: Unsteady convective exchange flows in cavities, J. Fluid Mech., 368, 127–153,, 1998. 

Sturman, J. J., Ivey, G. N., and Taylor, J. R.: Convection in a long box driven by heating and cooling on the horizontal boundaries, J. Fluid Mech., 310, 61–87,, 1996. 

Sturman, J. J., Oldham, C. E., and Ivey, G. N.: Steady convective exchange flows down slopes, Aquat. Sci., 61, 260–278,, 1999. 

Talling, J. F.: Origin of stratification in an African Rift lake, Limnol. Oceanogr., 8, 68–78,, 1963. 

Thorpe, S. A., Lemmin, U., Perrinjaquet, C., and Fer, I.: Observations of the thermal structure of a lake using a submarine, Limnol. Oceanogr., 44, 1575–1582,, 1999. 

Ulloa, H. N., Ramón, C. L., Doda, T., Wüest, A., and Bouffard, D.: Development of overturning circulation in sloping waterbodies due to surface cooling, J. Fluid Mech., 930, A18,, 2022. 

Verburg, P., Antenucci, J. P., and Hecky, R. E.: Differential cooling drives large-scale convective circulation in Lake Tanganyika, Limnol. Oceanogr., 56, 910–926,, 2011. 

Wedderburn, E. M.: The temperature of the fresh-water Lochs of Scotland, with special reference to Loch Ness, Trans. Roy. Soc. Edinb., 45, 407–489,, 1907. 

Wells, M. G. and Sherman, B.: Stratification produced by surface cooling in lakes with significant shallow regions, Limnol. Oceanogr., 46, 1747–1759,, 2001. 

Wetzel, R. G.: Limnology: Lake and River Ecosystems, 3rd Edn., Academic Press, San Diego, California, USA, ISBN 13 9780127447605, 2001. 

Woodward, B. L., Marti, C. L., Imberger, J., Hipsey, M. R., and Oldham, C. E.: Wind and buoyancy driven horizontal exchange in shallow embayments of a tropical reservoir: Lake Argyle, Western Australia, Limnol. Oceanogr., 62, 1636–1657,, 2017. 

Wüest, A. and Lorke, A.: Small-scale hydrodynamics in lakes, Annu. Rev. Fluid Mech., 35, 373–412,, 2003. 

Zilitinkevic, S. S.: Turbulent Penetrative Convection, edited by: Henderson-Sellers, B., Avebury, Adelshot, UK, ISBN 13 9781856281868, 1991. 

Zimmermann, M., Mayr, M. J., Bürgmann, H., Eugster, W., Steinsberger, T., Wehrli, B., Brand, A., and Bouffard, D.: Microbial methane oxidation efficiency and robustness during lake overturn, Limnol. Oceanogr. Lett., 6, 320–328,, 2021. 

Short summary
At night or during cold periods, the shallow littoral region of lakes cools faster than their deeper interior. This induces a cold downslope current that carries littoral waters offshore. From a 1-year-long database collected in a small temperate lake, we resolve the seasonality of this current and report its frequent occurrence from summer to winter. This study contributes to a better quantification of lateral exchange in lakes, with implications for the transport of dissolved compounds.