Abstract. In recent decades, several lake models of varying complexity have been developed and incorporated into numerical weather prediction systems and climate models. To foster enhanced forecasting ability and verification, improvement of these lake models remains essential. This especially applies to the limited simulation capabilities of biogeochemical processes in lakes and greenhouse gas exchanges with the atmosphere. Here we present multi-model simulations of physical variables and dissolved gas dynamics in a temperate lake (Harp Lake, Canada). The five models (ALBM, FLake, LAKE, LAKEoneD, MTCR-1) considered within this most recent round of the Lake Model Intercomparison Project (LakeMIP) all captured the seasonal temperature variability well. In contrast, none of the models is able to reproduce the exact dates of ice cover and ice off, leading to considerable errors in the simulation of eddy diffusivity around those dates. We then conducted an additional modeling experiment with a diffusing passive tracer to isolate the effect of the eddy diffusivity on gas concentration. Remarkably, sophisticated k−ε models do not demonstrate a significant difference in the vertical diffusion of a passive tracer compared to models with much simpler turbulence closures. All the models simulate less intensive spring overturn compared to autumn. Reduced mixing in the models consequently leads to the accumulation of the passive tracer distribution in the water column. The lake models with a comprehensive biogeochemical module, such as the ALBM and LAKE, predict dissolved oxygen dynamics adequate to the observed data. However, for the surface carbon dioxide concentration the correlation between modeled (ALBM, LAKE) and observed data is weak (∼0.3). Overall our results indicate the need to improve the representation of physical and biogeochemical processes in lake models, thereby contributing to enhanced weather prediction and climate projection capabilities.

Abstract. In recent decades, several lake models of varying complexity have been developed and incorporated in numerical weather prediction systems and climate models. To foster enhanced forecasting ability and verification, improvement of these lake models remains essential. This especially applies to the limited simulation capabilities of biogeochemical processes in lakes and greenhouse gas exchanges with the atmosphere. Here we present multi-model simulations of physical variables and dissolved gas dynamics in a temperate lake (Harp Lake, Canada). The five models (ALBM, FLake, LAKE, LAKEoneD, MTCR-1) considered within this most recent round of the Lake Model Intercomparison Project (LakeMIP) all captured the seasonal temperature variability well. In contrast, none of the models is able to reproduce the exact dates of ice-cover and ice-off, leading to considerable errors in the simulation of eddy diffusivity around those dates. We then conducted an additional modeling experiment with a diffusing passive tracer to isolate the effect of the eddy diffusivity on gas concentration. Remarkably, sophisticated k − ϵ models do not demonstrate a significant difference in the vertical diffusion of a passive tracer compared to models with much simpler turbulence closures. All models simulate less intensive spring overturn compared to autumn. Reduced mixing in the models consequently leads to the accumulation of the passive tracer distribution in the water column. The lake models with a comprehensive biogeochemical module, such as ALBM and LAKE, predict dissolved oxygen dynamics adequate to the observed data. However, for the surface carbon dioxide concentration the correlation between modeled (ALBM, LAKE) and observed data is weak (∼ 0.3). Overall our results indicate the need to improve the representation of physical and biogeochemical processes in lake models, thereby contributing to enhanced weather prediction and climate projection capabilities.

Inland waters, such as lakes, reservoirs and rivers, are important sources of climate forcing trace gases. A key parameter that regulates the gas exchange between water and the atmosphere is the gas transfer velocity, which itself is controlled by near-surface turbulence in the water. While in lakes and reservoirs, near-surface turbulence is mainly driven by atmospheric forcing, in shallow rivers and streams it is generated by bottom friction of gravity-forced flow. Large rivers represent a transition between these two cases. Near-surface turbulence has rarely been measured in rivers and the drivers of turbulence have not been quantified. We analyzed continuous measurements of flow velocity and quantified turbulence as the rate of dissipation of turbulent kinetic energy over the ice-free season in a large regulated river in Northern Finland. Measured dissipation rates agreed with predictions from bulk parameters, including mean flow velocity, wind speed, surface heat flux, and with a one-dimensional numerical turbulence model. Values ranged from ∼10 −10 m 2 s −3 to 10 −5 m 2 s −3 . Atmospheric forcing or gravity was the dominant driver of near-surface turbulence for similar fraction of the time. Large variability in near-surface dissipation rate occurred at diel time scales, when the flow velocity was strongly affected by downstream dam operation. By combining scaling relations for boundary-layer turbulence at the river bed and at the air-water interface, we derived a simple model for estimating the relative contributions of wind speed and bottom friction of river flow as a function of depth. Plain Language SummaryInland water bodies such as lakes, reservoirs and rivers are an important source of climate forcing trace gases to the atmosphere. Gas exchange between water and the atmosphere is regulated by the gas transfer velocity and the concentration difference between the water surface and the atmosphere. The gas transfer velocity depends on near-surface turbulence, but robust formulations have not been developed for river systems. Their surface area is sufficiently large for meteorological forcing to cause turbulence, as in lakes and reservoirs, but turbulence generated from bed and internal friction of gravity-driven flows is also expected to contribute. Here we quantify nearsurface turbulence using data from continuous air and water side measurements conducted over the ice-free season in a large subarctic regulated river in Finland. We find that turbulence, quantified as the dissipation rate of turbulent kinetic energy, is well described using equations for predicting turbulence from meteorological data for sufficiently high wind speeds whereas the contribution from bottom shear dominated at higher flow velocities. A one-dimensional river model successfully captured these processes.

No abstract