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

Oxygen conditions in ice-covered lakes depend on many factors, which, in turn, are influenced by a changing climate, so detection of the oxygen trend becomes difficult. Our research was based on data of long-term measurements of dissolved oxygen (2007–2020), water temperature, under-ice solar radiation, and snow-ice thickness (1995–2020) in Lake Vendyurskoe (Northwestern Russia). Changes of air temperature and precipitation in the study region during 1994–2020 and ice phenology of Lake Vendyurskoe for the same period based on field data and FLake model calculations were analyzed. The interannual variability of ice-on and ice-off dates covered wide time intervals (5 and 3 weeks, respectively), but no significant trends were revealed. In years with early ice-on, oxygen content decreased by more than 50% by the end of winter. In years with late ice-on and intermediate ice-off, the oxygen decrease was less than 40%. A significant negative trend was revealed for snow-ice cover thickness in spring. A climatic decrease of snow-ice cover thickness contributes to the rise of under-ice irradiance and earlier onset of under-ice convection. In years with early and long convection, an increase in oxygen content by 10–15% was observed at the end of the ice-covered period, presumably due to photosynthesis of phytoplankton.

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

Lake ice has a major impact on the functioning of lake ecosystems, the thermal and gas regimes of lakes, habitat conditions, socio-economic aspects of human life, local climate, etc. The multifaceted influence of lake ice makes it important to study its changes associated with global warming, including lake ice phenology, ice thickness, and the snow–ice fraction. This article presents a study of lake ice changes in different regions of Eurasia: the Arctic (Lake Imandra in the Murmansk region and Lake Kilpisjärvi in Finland), the temperate zone (six small and medium lakes in Karelia, Mozhaysk Reservoir in the Moscow region, and Lake Pääjärvi in Finland), the arid zone (Lake Ulansuhai in China), and the highlands (lakes Arpi and Sevan in Armenia). In the study regions, a statistically significant increase in winter air temperature has been recorded over the past few decades. The number of days with thaw (air temperature above 0 °C) has increased, while the number of days with severe frost (air temperature below −10 °C and −20 °C) has decreased. The share of liquid or mixed precipitation in winter increases most rapidly in the temperate zone. For two Finnish lakes, lakes Vendyurskoe and Vedlozero in Karelia, and Mozhaysk Reservoir, a decrease in the duration of the ice period was revealed, with later ice-on and earlier ice-off. The most dramatic change occurred in the large high-mountain Lake Sevan, where the water area has no longer been completely covered with ice every winter. In contrast, the small high-mountain Lake Arpi showed no significant changes in ice phenology over a 50-year period. Changes in the ice composition with an increase in the proportion of white ice and a decrease in the proportion of black ice have occurred in some lakes. In the temperate lakes Pääjärvi and Vendyurskoe, inverse dependences of the thickness of black ice on the number of days with thaw and frost in December–March for the first lake and on the amount of precipitation in the first month of ice for the second were observed. In the arid study region of China, due to the very little winter precipitation (usually less than 10 mm) only black ice occurs, and significant interannual variability in its thickness has been identified.

Until now, the phenomenon of radiatively driven convection (RDC) in ice-covered lakes has not been sufficiently studied, despite its important role in the functioning of aquatic ecosystems. There have been very few attempts to numerically simulate RDC due to the complexity of this process and the need to use powerful computing resources. The article presents the results of Large Eddy Simulations (LES) of RDC with periodic external energy pumping, which imitates the diurnal variations in solar radiation in the subglacial layer of lakes in spring. The research is aimed at numerically studying the initial stages in the formation and development of a convective mixed layer (CML). A numerical calculation was carried out for three variants of external energy pumping that differed in intensity. A diurnal acceleration and suppression of RDC due to a change in external pumping was revealed for all three variants. The results of numerical simulations provide estimates of such integral parameters of RDC development as the rate of deepening of the lower boundary of the CML, and the rate of water temperature rise within this layer. It was shown that as the cumulative heating of the CML increases over several days, daily increments in temperature and depth slowed down; that is, the dependence of the integral RDC parameters on external pumping was nonlinear. The LES results on RDC parameters were in good agreement with our observational data.

In winter, the reflective properties of ice-and-snow-covered lakes and the surrounding surface covered with snow differ little. Estimation of the albedo becomes especially important in spring when the reflectivity of the ice-covered lakes differs markedly from the surrounding snow-free surface. The measurements of the snow and ice albedo in situ are necessary to ensure adequate analysis of remote sensing data and monitoring of the state of the land ice. A wide range of surface albedo of small boreal ice-covered Lake Vendyurskoe (North-Western Russia) is calculated from observational data of solar radiation obtained in 1995–2000 and 2002–2017. A noticeable variability of albedo due to spring melting and changes in weather conditions is described. Climatic variability (North Atlantic Oscillation (NAO) index, the average monthly air temperature and sum of precipitation in December-March) is analyzed in the context of possible influence on the maximal thickness of snow and white ice, and, consequently, on the spring lake albedo dynamics. In current work to estimate the external heat flux that comes from the atmosphere to the surface of the lake and presumably determines the dynamics of albedo during spring melting, a simple parameter was used: the accumulated sum of positive temperatures, starting from March 1. The obtained prognostic equation of the cubic dependence of albedo on the named parameter demonstrates a different rate of albedo decrease in the melting of snow, white, and congelation ice.