Ross W. Griffiths
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Research Interests
Ross W. Griffiths has a total of 108 co-authors. Here are the top co-authors by number of shared publications:
- Graham O. Hughes
Imperial College London
23 shared publications - Bishakhdatta Gayen
Centre for Atmospheric and Oceanic Sciences Indian Institute of Science Bengaluru India
11 shared publications - Ian H. Campbell
Research School of Earth Sciences, The Australian National University, Canberra, Australia
10 shared publications Jonathan H. Fink
8 shared publications- Ross C. Kerr
Research School of Earth Sciences, The Australian National University, Canberra, ACT 2601, Australia
8 shared publications - Andrew McC. Hogg
Research School of Earth Sciences and ARC Centre of Excellence for Climate Extremes Australian National University Canberra ACT Australia
7 shared publications
Latest Publications
- 26 February 2019
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 865, 681-719
Transport by deep convection in basin-scale geostrophic circulation: turbulence-resolving simulations
Direct numerical simulations are used to investigate the nature of fully resolved small-scale convection and its role in large-scale circulation in a rotating $f$ -plane rectangular basin with imposed surface temperature difference. The large-scale circulation has a horizontal geostrophic component and a deep vertical overturning. This paper focuses on convective circulation with no wind stress, and buoyancy forcing sufficiently strong to ensure turbulent convection within the thermal boundary layer (horizontal Rayleigh numbers $Ra\approx 10^{12}{-}10^{13}$ ). The dynamics are found to depend on the value of a convective Rossby number, $Ro_{\unicode[STIX]{x0394}T}$ , which represents the strength of buoyancy forcing relative to Coriolis forces. Vertical convection shifts from a mean endwall plume under weak rotation ( $Ro_{\unicode[STIX]{x0394}T}>10^{-1}$ ) to ‘open ocean’ chimney convection plus mean vertical plumes at the side boundaries under strong rotation ( $Ro_{\unicode[STIX]{x0394}T}<10^{-1}$ ). The overall heat throughput, horizontal gyre transport and zonally integrated overturning transport are then consistent with scaling predictions for flow constrained by thermal wind balance in the thermal boundary layer coupled to vertical advection–diffusion balance in the boundary layer. For small Rossby numbers relevant to circulation in an ocean basin, vertical heat transport from the surface layer into the deep interior occurs mostly in ‘open ocean’ chimney convection while most vertical mass transport is against the side boundaries. Both heat throughput and the mean circulation (in geostrophic gyres, boundary currents and overturning) are reduced by geostrophic constraints.- 26 February 2019
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 865, 681-719
Transport by deep convection in basin-scale geostrophic circulation: turbulence-resolving simulations
Direct numerical simulations are used to investigate the nature of fully resolved small-scale convection and its role in large-scale circulation in a rotating $f$ -plane rectangular basin with imposed surface temperature difference. The large-scale circulation has a horizontal geostrophic component and a deep vertical overturning. This paper focuses on convective circulation with no wind stress, and buoyancy forcing sufficiently strong to ensure turbulent convection within the thermal boundary layer (horizontal Rayleigh numbers $Ra\approx 10^{12}{-}10^{13}$ ). The dynamics are found to depend on the value of a convective Rossby number, $Ro_{\unicode[STIX]{x0394}T}$ , which represents the strength of buoyancy forcing relative to Coriolis forces. Vertical convection shifts from a mean endwall plume under weak rotation ( $Ro_{\unicode[STIX]{x0394}T}>10^{-1}$ ) to ‘open ocean’ chimney convection plus mean vertical plumes at the side boundaries under strong rotation ( $Ro_{\unicode[STIX]{x0394}T}<10^{-1}$ ). The overall heat throughput, horizontal gyre transport and zonally integrated overturning transport are then consistent with scaling predictions for flow constrained by thermal wind balance in the thermal boundary layer coupled to vertical advection–diffusion balance in the boundary layer. For small Rossby numbers relevant to circulation in an ocean basin, vertical heat transport from the surface layer into the deep interior occurs mostly in ‘open ocean’ chimney convection while most vertical mass transport is against the side boundaries. Both heat throughput and the mean circulation (in geostrophic gyres, boundary currents and overturning) are reduced by geostrophic constraints.- 28 January 2019
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 863, 545-571
Ablation of sloping ice faces into polar seawater
The effects of the slope of an ice–seawater interface on the mechanisms and rate of ablation of the ice by natural convection are examined using turbulence-resolving simulations. Solutions are obtained for ice slopes $\unicode[STIX]{x1D703}=2^{\circ }{-}90^{\circ }$ , at a fixed ambient salinity and temperature, chosen to represent common Antarctic ocean conditions. For laminar boundary layers the ablation rate decreases with height, whereas in the turbulent regime the ablation rate is found to be height independent. The simulated laminar ablation rates scale with $(\sin \unicode[STIX]{x1D703})^{1/4}$ , whereas in the turbulent regime it follows a $(\sin \unicode[STIX]{x1D703})^{2/3}$ scaling, both consistent with the theoretical predictions developed here. The reduction in the ablation rate with shallower slopes arises as a result of the development of stable density stratification beneath the ice face, which reduces turbulent buoyancy fluxes to the ice. The turbulent kinetic energy budget of the flow shows that, for very steep slopes, both buoyancy and shear production are drivers of turbulence, whereas for shallower slopes shear production becomes the dominant mechanism for sustaining turbulence in the convective boundary layer.- 28 January 2019
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 863, 545-571
Ablation of sloping ice faces into polar seawater
The effects of the slope of an ice–seawater interface on the mechanisms and rate of ablation of the ice by natural convection are examined using turbulence-resolving simulations. Solutions are obtained for ice slopes $\unicode[STIX]{x1D703}=2^{\circ }{-}90^{\circ }$ , at a fixed ambient salinity and temperature, chosen to represent common Antarctic ocean conditions. For laminar boundary layers the ablation rate decreases with height, whereas in the turbulent regime the ablation rate is found to be height independent. The simulated laminar ablation rates scale with $(\sin \unicode[STIX]{x1D703})^{1/4}$ , whereas in the turbulent regime it follows a $(\sin \unicode[STIX]{x1D703})^{2/3}$ scaling, both consistent with the theoretical predictions developed here. The reduction in the ablation rate with shallower slopes arises as a result of the development of stable density stratification beneath the ice face, which reduces turbulent buoyancy fluxes to the ice. The turbulent kinetic energy budget of the flow shows that, for very steep slopes, both buoyancy and shear production are drivers of turbulence, whereas for shallower slopes shear production becomes the dominant mechanism for sustaining turbulence in the convective boundary layer.- 13 October 2017
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 831, 491-523
Turbulent horizontal convection under spatially periodic forcing: a regime governed by interior inertia
Differential heating applied at a single horizontal boundary forces ‘horizontal convection’, even when there is no net heat flux through the boundary. However, almost all studies of horizontal convection have been limited to a special class of problem in which temperature or heat flux differences were applied in only one direction and over the horizontal length of a box (the Rossby problem; Rossby, Deep-Sea Res., vol. 12, 1965, pp. 9–16). These conditions strongly constrain the flow. Here we report laboratory experiments and direct numerical simulations (DNS) extending the results of Griffiths & Gayen (Phys. Rev. Lett., vol. 115, 2015, 204301) for horizontal convection forced by boundary conditions imposed in a two-dimensional periodic array at a horizontal boundary. The experiments use saline and freshwater fluxes at a permeable base with the imposed boundary salinity having a horizontal length scale one quarter of the width of the box. The flow reaches a state in which the net boundary buoyancy flux vanishes and the bulk of the fluid shows an inertial range of turbulence length scales. A regime transition is seen for increasing water depth, from an array of individual coherent plumes on the forcing scale to convection dominated by emergent larger scales of overturning. The DNS explore the analogous thermally forced case with sinusoidal boundary temperature of wavenumber $n=4$ , and are used to examine the Rayleigh number ( $Ra$ ) dependence for shallow- and deep-water cases. For shallow water the flow transitions with increasing $Ra$ from laminar to turbulent boundary layer regimes that are familiar from the Rossby problem and which have normalised heat transport scaling as $Nu\sim Ra^{1/5}$ and $Nu\sim (Ra\,Pr)^{1/5}$ , with $Nu$ the Nusselt number and $Pr$ the Prandtl number, in this case maintaining a stable array of coherent turbulent plumes. For deep-water and large $Ra$ the laminar scaling transitions to $Nu\sim (Ra\,Pr)^{1/4}$ , with the scales of turbulence extending to the dimensions of the box. The $1/4$ power law regime is explained in terms of the momentum of symmetric, inviscid large scales of motion in the interior coupled to diffusive loss of heat through stabilised parts of the boundary layer. The turbulence production is predominantly by shear instability rather than convection, with viscous dissipation distributed throughout the bulk of the fluid. These conditions are not seen in the highly asymmetric flow in the Rossby problem even at Rayleigh numbers up to six orders of magnitude greater than the transition found here. The new inertial interior regime has the rate of supply of available potential energy, and its removal by mixing of density, increasing as $Ra^{5/4}$ , which is faster than $Ra^{6/5}$ in the Rossby problem. Irreversible mixing is confined close to the forcing boundary and is very much larger than the viscous dissipation, which is proportional to $Ra$ .- 13 October 2017
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 831, 491-523
Turbulent horizontal convection under spatially periodic forcing: a regime governed by interior inertia
Differential heating applied at a single horizontal boundary forces ‘horizontal convection’, even when there is no net heat flux through the boundary. However, almost all studies of horizontal convection have been limited to a special class of problem in which temperature or heat flux differences were applied in only one direction and over the horizontal length of a box (the Rossby problem; Rossby, Deep-Sea Res., vol. 12, 1965, pp. 9–16). These conditions strongly constrain the flow. Here we report laboratory experiments and direct numerical simulations (DNS) extending the results of Griffiths & Gayen (Phys. Rev. Lett., vol. 115, 2015, 204301) for horizontal convection forced by boundary conditions imposed in a two-dimensional periodic array at a horizontal boundary. The experiments use saline and freshwater fluxes at a permeable base with the imposed boundary salinity having a horizontal length scale one quarter of the width of the box. The flow reaches a state in which the net boundary buoyancy flux vanishes and the bulk of the fluid shows an inertial range of turbulence length scales. A regime transition is seen for increasing water depth, from an array of individual coherent plumes on the forcing scale to convection dominated by emergent larger scales of overturning. The DNS explore the analogous thermally forced case with sinusoidal boundary temperature of wavenumber $n=4$ , and are used to examine the Rayleigh number ( $Ra$ ) dependence for shallow- and deep-water cases. For shallow water the flow transitions with increasing $Ra$ from laminar to turbulent boundary layer regimes that are familiar from the Rossby problem and which have normalised heat transport scaling as $Nu\sim Ra^{1/5}$ and $Nu\sim (Ra\,Pr)^{1/5}$ , with $Nu$ the Nusselt number and $Pr$ the Prandtl number, in this case maintaining a stable array of coherent turbulent plumes. For deep-water and large $Ra$ the laminar scaling transitions to $Nu\sim (Ra\,Pr)^{1/4}$ , with the scales of turbulence extending to the dimensions of the box. The $1/4$ power law regime is explained in terms of the momentum of symmetric, inviscid large scales of motion in the interior coupled to diffusive loss of heat through stabilised parts of the boundary layer. The turbulence production is predominantly by shear instability rather than convection, with viscous dissipation distributed throughout the bulk of the fluid. These conditions are not seen in the highly asymmetric flow in the Rossby problem even at Rayleigh numbers up to six orders of magnitude greater than the transition found here. The new inertial interior regime has the rate of supply of available potential energy, and its removal by mixing of density, increasing as $Ra^{5/4}$ , which is faster than $Ra^{6/5}$ in the Rossby problem. Irreversible mixing is confined close to the forcing boundary and is very much larger than the viscous dissipation, which is proportional to $Ra$ .- 15 June 2017
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 823, 57-99
Geostrophic and chimney regimes in rotating horizontal convection with imposed heat flux
Convection in a rotating rectangular basin with differential thermal forcing at one horizontal boundary is examined using laboratory experiments. The experiments have an imposed heat flux boundary condition, are at large values of the flux Rayleigh number ( $Ra_{F}\sim O(10^{13}{-}10^{14})$ based on the box length $L$ ), use water with Prandtl number $Pr\approx 4$ and have a small depth to length aspect ratio. The results show the conditions for transition from non-rotating horizontal convection governed by an inertial–buoyancy balance in the thermal boundary layer, to circulation governed by geostrophic flow in the boundary layer. The geostrophic balance constrains mean flow and reduces the heat transport as Nusselt number $Nu\sim (Ra_{F}Ro)^{1/6}$ , where $Ro=B^{1/2}/f^{3/2}L$ is the convective Rossby number, $B$ is the imposed buoyancy flux and $f$ is the Coriolis parameter. Thus flow in the geostrophic boundary layer regime is governed by the relative roles of horizontal convective accelerations and Coriolis accelerations, or buoyancy and rotation, in the boundary layer. Experimental evidence suggests that for more rapid rotation there is another transition to a regime in which the momentum budget is dominated by fluctuating vertical accelerations in a region of vortical plumes, which we refer to as a ‘chimney’ following related discussion of regions of deep convection in the ocean. Coupling of the chimney convection in the region of destabilising boundary flux to the diffusive boundary layer of horizontal convection in the region of stabilising boundary flux gives heat transport independent of rotation in this ‘inertial chimney’ regime, and the new scaling $Nu\sim Ra_{F}^{1/4}$ . Scaling analysis predicts the transition conditions observed in the experiments, as well as a further ‘geostrophic chimney’ regime in which the vertical plumes are controlled by local geostrophy. When $Ro<10^{-1}$ , the convection is also observed to produce a set of large basin-scale gyres at all depths in the time-averaged flow.- 15 June 2017
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 823, 57-99
Geostrophic and chimney regimes in rotating horizontal convection with imposed heat flux
Convection in a rotating rectangular basin with differential thermal forcing at one horizontal boundary is examined using laboratory experiments. The experiments have an imposed heat flux boundary condition, are at large values of the flux Rayleigh number ( $Ra_{F}\sim O(10^{13}{-}10^{14})$ based on the box length $L$ ), use water with Prandtl number $Pr\approx 4$ and have a small depth to length aspect ratio. The results show the conditions for transition from non-rotating horizontal convection governed by an inertial–buoyancy balance in the thermal boundary layer, to circulation governed by geostrophic flow in the boundary layer. The geostrophic balance constrains mean flow and reduces the heat transport as Nusselt number $Nu\sim (Ra_{F}Ro)^{1/6}$ , where $Ro=B^{1/2}/f^{3/2}L$ is the convective Rossby number, $B$ is the imposed buoyancy flux and $f$ is the Coriolis parameter. Thus flow in the geostrophic boundary layer regime is governed by the relative roles of horizontal convective accelerations and Coriolis accelerations, or buoyancy and rotation, in the boundary layer. Experimental evidence suggests that for more rapid rotation there is another transition to a regime in which the momentum budget is dominated by fluctuating vertical accelerations in a region of vortical plumes, which we refer to as a ‘chimney’ following related discussion of regions of deep convection in the ocean. Coupling of the chimney convection in the region of destabilising boundary flux to the diffusive boundary layer of horizontal convection in the region of stabilising boundary flux gives heat transport independent of rotation in this ‘inertial chimney’ regime, and the new scaling $Nu\sim Ra_{F}^{1/4}$ . Scaling analysis predicts the transition conditions observed in the experiments, as well as a further ‘geostrophic chimney’ regime in which the vertical plumes are controlled by local geostrophy. When $Ro<10^{-1}$ , the convection is also observed to produce a set of large basin-scale gyres at all depths in the time-averaged flow.- 31 May 2016
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 798, 284-298
Simulation of convection at a vertical ice face dissolving into saline water
We investigate the convection and dissolution rate generated when a wall of ice dissolves into seawater under Antarctic Ocean conditions. In direct numerical simulations three coupled interface equations are used to solve for interface temperature, salinity and ablation velocity, along with the boundary layer flow and transport. The main focus is on ambient water temperatures between $-1\,^{\circ }\text{C}$ and $6\,^{\circ }\text{C}$ and salinities around 35 ‰, where diffusion of salt to the ice–water interface depresses the freezing point and enhances heat diffusion to the ice. We show that fluxes of both heat and salt to the interface are significant in governing the dissolution of ice, and the ablation velocity agrees well with experiments and a recent theoretical prediction. The same turbulent flow dynamics and ablation rate are expected to apply at any depth in a deeper ocean water column (after choosing the relevant pressure coefficient for the liquidus temperature). At Grashof numbers currently accessible by direct numerical simulation, turbulence is generated both directly from buoyancy flux and from shear production in the buoyancy-driven boundary layer flow, whereas shear production by the convective flow is expected to be more important at geophysical scales. The momentum balance in the boundary layer is dominated by buoyancy forcing and wall stress, with the latter characterised by a large drag coefficient.- 31 May 2016
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 798, 284-298
Simulation of convection at a vertical ice face dissolving into saline water
We investigate the convection and dissolution rate generated when a wall of ice dissolves into seawater under Antarctic Ocean conditions. In direct numerical simulations three coupled interface equations are used to solve for interface temperature, salinity and ablation velocity, along with the boundary layer flow and transport. The main focus is on ambient water temperatures between $-1\,^{\circ }\text{C}$ and $6\,^{\circ }\text{C}$ and salinities around 35 ‰, where diffusion of salt to the ice–water interface depresses the freezing point and enhances heat diffusion to the ice. We show that fluxes of both heat and salt to the interface are significant in governing the dissolution of ice, and the ablation velocity agrees well with experiments and a recent theoretical prediction. The same turbulent flow dynamics and ablation rate are expected to apply at any depth in a deeper ocean water column (after choosing the relevant pressure coefficient for the liquidus temperature). At Grashof numbers currently accessible by direct numerical simulation, turbulence is generated both directly from buoyancy flux and from shear production in the buoyancy-driven boundary layer flow, whereas shear production by the convective flow is expected to be more important at geophysical scales. The momentum balance in the boundary layer is dominated by buoyancy forcing and wall stress, with the latter characterised by a large drag coefficient.- 1 March 2016
- Published by American Meteorological Society in Journal of Physical Oceanography
- Vol. 46 (3) , 731-747
Adjustment of the Meridional Overturning Circulation and Its Dependence on Depth of Mixing
The relative roles of advective processes and mixing on the temporal adjustment of the meridional overturning circulation are examined, in particular the effects of mixing in either the abyssal or upper ocean. Laboratory experiments with convectively driven overturning and imposed stirring rates show that the circulation adjusts toward an equilibrium state on time scales governed by mixing in the upper boundary layer region but independent of the mixing rate in the bulk of the interior. The equilibrium state of the stratification is dependent only on the rate of mixing in the boundary layer. An idealized high-resolution ocean model shows adjustment (of a two-cell circulation) dominated primarily by the advective ventilation time scale, consistent with a view of the circulation determined by water mass transformation occurring primarily near the surface. Both the experiments and the model results indicate that adjustments of the circulation are controlled by surface buoyancy uptake (or rejection) and that the nonequilibrium circulation is dominated by advective processes, especially if the average abyssal ocean diffusivity is less than 3 × 10−5 m2 s−1.
- 1 March 2016
- Published by American Meteorological Society in Journal of Physical Oceanography
- Vol. 46 (3) , 731-747
Adjustment of the Meridional Overturning Circulation and Its Dependence on Depth of Mixing
The relative roles of advective processes and mixing on the temporal adjustment of the meridional overturning circulation are examined, in particular the effects of mixing in either the abyssal or upper ocean. Laboratory experiments with convectively driven overturning and imposed stirring rates show that the circulation adjusts toward an equilibrium state on time scales governed by mixing in the upper boundary layer region but independent of the mixing rate in the bulk of the interior. The equilibrium state of the stratification is dependent only on the rate of mixing in the boundary layer. An idealized high-resolution ocean model shows adjustment (of a two-cell circulation) dominated primarily by the advective ventilation time scale, consistent with a view of the circulation determined by water mass transformation occurring primarily near the surface. Both the experiments and the model results indicate that adjustments of the circulation are controlled by surface buoyancy uptake (or rejection) and that the nonequilibrium circulation is dominated by advective processes, especially if the average abyssal ocean diffusivity is less than 3 × 10−5 m2 s−1.
- 9 November 2015
- Published by American Physical Society (APS) in Physical Review Letters
- Vol. 115 (20) , 204301
Turbulent Convection Insights from Small-Scale Thermal Forcing with Zero Net Heat Flux at a Horizontal Boundary
A large-scale circulation, a turbulent boundary layer, and a turbulent plume are noted features of convection at large Rayleigh numbers under differential heating on a single horizontal boundary. These might be attributed to the forcing, which in all studies has been limited to a unidirectional gradient over the domain scale. We instead apply forcing on a length scale smaller than the domain, and with variation in both horizontal directions. Direct numerical simulations show turbulence throughout the domain, a regime transition to a dominant domain-scale circulation, and a region of logarithmic velocity in the boundary layer, despite zero net heat flux. The results show significant similarities to Rayleigh-Bénard convection, demonstrate the significance of plume merging, support the hypothesis that the key driver of convection is the production of available potential energy without necessarily supplying total potential energy, and imply that contributions to domain-scale circulation in the oceans need not be solely from the large-scale gradients of forcing.
- 9 November 2015
- Published by American Physical Society (APS) in Physical Review Letters
- Vol. 115 (20) , 204301
Turbulent Convection Insights from Small-Scale Thermal Forcing with Zero Net Heat Flux at a Horizontal Boundary
A large-scale circulation, a turbulent boundary layer, and a turbulent plume are noted features of convection at large Rayleigh numbers under differential heating on a single horizontal boundary. These might be attributed to the forcing, which in all studies has been limited to a unidirectional gradient over the domain scale. We instead apply forcing on a length scale smaller than the domain, and with variation in both horizontal directions. Direct numerical simulations show turbulence throughout the domain, a regime transition to a dominant domain-scale circulation, and a region of logarithmic velocity in the boundary layer, despite zero net heat flux. The results show significant similarities to Rayleigh-Bénard convection, demonstrate the significance of plume merging, support the hypothesis that the key driver of convection is the production of available potential energy without necessarily supplying total potential energy, and imply that contributions to domain-scale circulation in the oceans need not be solely from the large-scale gradients of forcing.
- Open Access
- 1 January 2015
- Published by Elsevier in Procedia IUTAM
- Vol. 15, 165-171
Coupling of Cooling, Solidification and Gravity-driven Flow
- 1 January 2015
- Published by Elsevier in Procedia IUTAM
- Vol. 15, 165-171
Open AccessCoupling of Cooling, Solidification and Gravity-driven Flow
- Open Access
- 1 January 2015
- Published by Elsevier in Procedia IUTAM
- Vol. 15, 78-85
Melting Driven Convection at the Ice-seawater Interface
- 1 January 2015
- Published by Elsevier in Procedia IUTAM
- Vol. 15, 78-85
Open AccessMelting Driven Convection at the Ice-seawater Interface
- 12 October 2014
- Published by Wiley in Limnology and Oceanography
- Vol. 59 (6) , 2121-2132
Interannual variability of circulation under spring ice in a boreal lake
A small range (∼ 1°C) of under‐ice water temperature is shown to result in remarkably different circulation regimes under spring ice in a deep, oligotrophic boreal lake. With the water column at < 4°C, melting of snow led to deepening vertical convection before ice break and a final depth of convection inversely correlated with earlier deep‐water temperature. We attribute that to the nonlinear dependence of water density on temperature, albeit further affected by stochastic weather factors. In four of nine study years, convection led to complete under‐ice overturn of the lake, indicating that this may not be uncommon in similar lakes with steep topography. River inflow and more intensive warming of water in the littoral zone also created a horizontal density differential, convection that involved flow down the sloping bottom and a lateral intrusion of this sinking water at a depth between the vertical convection and the quiescent deep‐water layers. The vertical and horizontal convection together produced a profile of temperature slightly increasing from the surface to the bottom of the convection layer. The contribution of horizontal convection to under‐ice mixing was interannually variable, and in one of the study years it eventually dominated under‐ice mixing. A thermal bar circulation regime developed occasionally and only in the open water between ice and shoreline. We identified five different under‐ice mixing regimes that form an interannually variable continuum of behavior during the ice melting period. The dependence on a narrow temperature range likely makes the circulation regime sensitive to a warming climate.
- 12 October 2014
- Published by Wiley in Limnology and Oceanography
- Vol. 59 (6) , 2121-2132
Interannual variability of circulation under spring ice in a boreal lake
A small range (∼ 1°C) of under‐ice water temperature is shown to result in remarkably different circulation regimes under spring ice in a deep, oligotrophic boreal lake. With the water column at < 4°C, melting of snow led to deepening vertical convection before ice break and a final depth of convection inversely correlated with earlier deep‐water temperature. We attribute that to the nonlinear dependence of water density on temperature, albeit further affected by stochastic weather factors. In four of nine study years, convection led to complete under‐ice overturn of the lake, indicating that this may not be uncommon in similar lakes with steep topography. River inflow and more intensive warming of water in the littoral zone also created a horizontal density differential, convection that involved flow down the sloping bottom and a lateral intrusion of this sinking water at a depth between the vertical convection and the quiescent deep‐water layers. The vertical and horizontal convection together produced a profile of temperature slightly increasing from the surface to the bottom of the convection layer. The contribution of horizontal convection to under‐ice mixing was interannually variable, and in one of the study years it eventually dominated under‐ice mixing. A thermal bar circulation regime developed occasionally and only in the open water between ice and shoreline. We identified five different under‐ice mixing regimes that form an interannually variable continuum of behavior during the ice melting period. The dependence on a narrow temperature range likely makes the circulation regime sensitive to a warming climate.