Seasonality of density currents induced by differential cooling
- 1Eawag, Swiss Federal Institute of Aquatic Science and Technology, Surface Waters – Research and Management, 6047 Kastanienbaum, Switzerland
- 2Physics of Aquatic Systems Laboratory, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
- 3Department of Civil Engineering, University of Granada, 18003 Granada, Spain
- 4Department of Earth and Environmental Science, University of Pennsylvania, Philadelphia, PA 19104-6316, USA
Correspondence: Tomy Doda (firstname.lastname@example.org)
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.
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).
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.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.
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 (x, z) 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 (47∘02′12′′ N, 8∘18′04′′ 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 is defined as
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:
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 under the assumption that only heat fluxes modify the potential energy near the surface (Bouffard and Wüest, 2019):
In Eq. (3), 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:
The net buoyancy flux at the surface is . indicates a destabilizing buoyancy flux (net cooling), whereas 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 ∘C m−1. The water column was defined as fully mixed when ∘C 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 . Starting from ht and moving upward, the lower end of the mixed layer was reached when the local temperature gradient dropped to ∘C m−1. When ∘C 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:
There was significant cross-shore flow, where the depth-averaged velocity of the current cm s−1 for at least 2 h;
There was weak wind during the selected event, where (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);
There was a mixed water column at MT before the onset of the flow, where ∘C 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:
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
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 (): 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 , 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 , 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:
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 , where h is a vertical length scale. We used the average depth of the littoral region upslope of MT (, Fig. 1) as h:
where cq is a proportionality coefficient. We defined the flushing timescale from Eq. (10) as
where . 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 . 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:
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.
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 W kg−1 during daytime to 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 W kg−1 with a limited variability of 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 ∘C m−1 in late summer (Fig. 1).
We divide the diurnal cycle into a net cooling phase during the night () and a net heating phase during the day (). 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.1 ∘C 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 (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 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 α.
From spring to summer, hML remains shallow (<2 m) and B0 increases from W kg−1 in March to W kg−1 in August. The latter results in an increase in the convective velocity scale from to mm s−1 over the same period (Fig. 4a, b). The lake undergoes net daily heating over this period (), with a cooling phase that remains shorter than 15 h (Fig. 4b). In August, 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 ( 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 mm s−1 from August to November. The duration of the cooling phase increases in autumn to reach its maximum in November ( h). Convective cooling is occurring almost continuously at that time. In winter, the convective velocity w* drops to mm s−1, and the effect of the cooling phase is again balanced by the heating phase (, h in February). The differences between the cooling and heating phases of the diurnal cycle are reduced in winter, with 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>4 ∘C).
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.
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 ( %) with the rest of the year ( %), 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, ; Fig. 6a, b). The daily average velocity cm s−1 is 3 times lower than the daily maximum velocity cm s−1 that occurs in the morning (Fig. 3). The best linear fits (with 0 intercept) are given by and . The 95 % confidence interval of the proportionality coefficient is for Uavg and 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.
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, ; Fig. 6c, d). The unit-width discharge m2 s−1 corresponds to an average thickness of m and a flushing timescale of h. The best relationships (with 0 intercept) are given by (R2=0.27) and (R2=0.18). The 95 % confidence interval of the proportionality coefficient is for qavg and 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 show that the magnitude of and variability in B0,net decrease over the same period (Fig. 7b, d, f, h, j).
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 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:
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.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.
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 , 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:
where [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:
In Eq. (15), the buoyancy frequency squared is approximated as . 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 ( W kg−1) lead to a short transition timescale of 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 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 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, , 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.
We further compare our scaling for Ux, qx, and τF with previous studies on buoyancy-driven flows. The velocity scale 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 , where m2 s−1 is the kinematic viscosity of water. The empirical relationship 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.
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 , , and , 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.
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 (), for each 24 h subset starting at 17:00 UTC. 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 between days with and without TSs (Fig. B1). Days with TSs are characterized by higher 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).
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.
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 ( 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., cm s−1 for at least 2 h) is lake-specific. This criterion is based on typical examples of TSs in Rotsee, where 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., 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 ( 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.
The filter used to discard days with wind-driven circulation is based on the criterion , where κ=0.41 is the von Kármán constant, and u* the friction velocity. The ratio 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 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).
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).
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 (https://doi.org/10.25678/00057K, Doda et al., 2021).
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.
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.
This study was financed by the Swiss National Science Foundation (“Buoyancy driven nearshore transport in lakes” project; HYPOlimnetic THErmal SIphonS, HYPOTHESIS, grant no. 175919).
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, https://doi.org/10.1061/(ASCE)0733-9429(1984)110:6(773), 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, https://doi.org/10.4081/jlimnol.2010.3, 2010.
Bouffard, D. and Wüest, A.: Convection in lakes, Annu. Rev. Fluid Mech., 51, 189–215, https://doi.org/10.1146/annurev-fluid-010518-040506, 2019.
Brink, K. H.: Cross-shelf exchange, Annu. Rev. Mar. Sci., 8, 59–78, https://doi.org/10.1146/annurev-marine-010814-015717, 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, https://doi.org/10.4319/lo.1986.31.3.0657, 1986.
Cyr, H.: Winds and the distribution of nearshore phytoplankton in a stratified lake, Water Res., 122, 114–127, https://doi.org/10.1016/j.watres.2017.05.066, 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, https://doi.org/10.1175/1520-0469(1970)027<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], https://doi.org/10.25678/00057K, 2021.
Eccles, D. H.: An outline of the physical limnology of Lake Malawi (Lake Nyasa), Limnol. Oceanogr., 19, 730–742, https://doi.org/10.4319/lo.1974.19.5.0730, 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, https://doi.org/10.1127/1863-9135/2010/0178-0001, 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, https://doi.org/10.4319/lo.2002.47.2.0576, 2002a.
Fer, I., Lemmin, U., and Thorpe, S. A.: Winter cascading of cold water in Lake Geneva, J. Geophys. Res., 107, 3060, https://doi.org/10.1029/2001JC000828, 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, https://doi.org/10.1002/2013WR014448, 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, https://doi.org/10.1017/S0022112098003437, 1999.
Forel, F. A.: Le Léman – Monographie Limnologique. Tome II, Editions Rouge, Lausanne, Switzerland, https://doi.org/10.5962/bhl.title.124608, 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, https://doi.org/10.4319/lo.2008.53.5_part_2.2321, 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, https://doi.org/10.1016/j.watres.2019.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, https://doi.org/10.1029/JC091iC11p13056, 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, https://doi.org/10.1002/grl.50453, 2013.
Horsch, G. M. and Stefan, H. G.: Convective circulation in littoral water due to surface cooling, Limnol. Oceanogr., 33, 1068–1083, https://doi.org/10.4319/lo.19126.96.36.1998, 1988.
Imberger, J. and Patterson, J. C.: Physical Limnology, edited by J. W. Hutchinson and T. Y. Wu, Adv. Appl. Mech., 27, 303–475, https://doi.org/10.1016/S0065-2156(08)70199-6, 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, https://doi.org/10.4319/lo.1991.36.1.0179, 1991a.
James, W. F. and Barko, J. W.: Littoral-pelagic phosphorus dynamics during nighttime convective circulation, Limnol. Oceanogr., 36, 949–960, https://doi.org/10.4319/lo.1991.36.5.0949, 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, https://doi.org/10.1007/BF00016857, 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, https://doi.org/10.1029/2008WR006949, 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, https://doi.org/10.2307/1467544, 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, https://doi.org/10.4319/lo.2002.47.3.0656, 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, https://doi.org/10.1016/j.envsoft.2011.11.017, 2012.
Meyers, T. and Dale, R.: Predicting daily insolation with hourly cloud height and coverage, J. Clim. Appl. Meteorol., 22, 537–545, https://doi.org/10.1175/1520-0450(1983)022<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, https://doi.org/10.4319/lo.19188.8.131.526, 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, https://doi.org/10.1098/rstb.1952.0005, 1952.
Pálmarsson, S. Ó. and Schladow, S. G.: Exchange flow in a shallow lake embayment, Ecol. Appl., 18, A89–A106, https://doi.org/10.1890/06-1618.1, 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, https://doi.org/10.4319/lo.2003.48.4.1419, 2003.
Phillips, O. M.: On turbulent convection currents and the circulation of the Red Sea, Deep-Sea Res. Oceanogr. Abstr., 13, 1149–1160, https://doi.org/10.1016/0011-7471(66)90706-6, 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, https://doi.org/10.3394/0380-1330(2007)33[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, https://doi.org/10.1007/BF00877450, 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, https://doi.org/10.1002/lom3.10326, 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, https://doi.org/10.1007/s00027-006-0823-8, 2007.
Schlatter, J. W., Wüest, A., and Imboden, D. M.: Hypolimnetic density currents traced by sulphur hexafluoride (SF6), Aquat. Sci., 59, 225–242, https://doi.org/10.1007/BF02523275, 1997.
Sturman, J. J. and Ivey, G. N.: Unsteady convective exchange flows in cavities, J. Fluid Mech., 368, 127–153, https://doi.org/10.1017/S002211209800175X, 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, https://doi.org/10.1017/S0022112096001735, 1996.
Sturman, J. J., Oldham, C. E., and Ivey, G. N.: Steady convective exchange flows down slopes, Aquat. Sci., 61, 260–278, https://doi.org/10.1007/s000270050065, 1999.
Talling, J. F.: Origin of stratification in an African Rift lake, Limnol. Oceanogr., 8, 68–78, https://doi.org/10.4319/lo.1963.8.1.0068, 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, https://doi.org/10.4319/lo.19184.108.40.2065, 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, https://doi.org/10.1017/jfm.2021.883, 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, https://doi.org/10.4319/lo.2011.56.3.0910, 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, https://doi.org/10.1017/S0080456800022791, 1907.
Wells, M. G. and Sherman, B.: Stratification produced by surface cooling in lakes with significant shallow regions, Limnol. Oceanogr., 46, 1747–1759, https://doi.org/10.4319/lo.2001.46.7.1747, 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, https://doi.org/10.1002/lno.10522, 2017.
Wüest, A. and Lorke, A.: Small-scale hydrodynamics in lakes, Annu. Rev. Fluid Mech., 35, 373–412, https://doi.org/10.1146/annurev.fluid.35.101101.161220, 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, https://doi.org/10.1002/lol2.10209, 2021.