Modelling the fate and transport of negatively buoyant storm-river water in small multi-basin lakes

The dynamics of negatively buoyant river plumes in a small multi-basin kettle lake with steep bathymetry (Toolik Lake, AK) are simulated using a Cartesian hydrodynamic model based on the solution of the three-dimensional shallow water equations. To validate the model, model predictions are compared with results from previous analytical and laboratory studies and with field observations. The grid resolution adopted for the Toolik Lake model is 0.5m (=@Dz) in the vertical and 20m (=@Dx) in the horizontal, so that the ratio of the bottom slope S"0 to @Dz/@Dx is lower than 4 in 99% of the computational domain. With that resolution, the model represents correctly the rate of mixing between lake and river water and the speed of propagation of downslope gravity currents. The model provides accurate predictions of the temperature structure (RMSE=0.25^oC) and of eddy diffusivities at the depths of the intrusions of incoming water. Measurements and modelling show similar fractions and depth distribution of river water on a cross-basin transect, which suggests that the mixing dynamics of the plume as it transits between basins are well resolved. Thus, the stage is set to quantify the ecological consequences of storm events in small lakes with several interconnected basins using coupled biological measurements and 3D modelling.

[1]  M. Gregg,et al.  Diurnal restratification and turbulence in the oceanic surface mixed layer: 1. Observations , 1993 .

[2]  R. Stocker,et al.  Modeling circulation in lakes: Spatial and temporal variations , 2003 .

[3]  G. Mellor,et al.  A Hierarchy of Turbulence Closure Models for Planetary Boundary Layers. , 1974 .

[4]  J. Antenucci,et al.  Influence of wind and lake morphometry on the interaction between two rivers entering a stratified lake , 2008 .

[5]  J. Imberger,et al.  Time-Varying Underflow into a Continuous Stratification with Bottom Slope , 2008 .

[6]  Kenneth Falconer,et al.  Unsolved Problems In Geometry , 1991 .

[7]  Rex Britter,et al.  The motion of the front of a gravity current travelling down an incline , 1980, Journal of Fluid Mechanics.

[8]  B. Peterson,et al.  Water and sediment export of the upper Kuparuk River drainage of the North Slope of Alaska , 1992, Hydrobiologia.

[9]  Waterways Experiment Station.,et al.  An assessment of reservoir mixing processes , 1986 .

[10]  Peter Yeates,et al.  Pseudo two‐dimensional simulations of internal and boundary fluxes in stratified lakes and reservoirs , 2003 .

[11]  P. Bournet,et al.  Numerical Investigation of Plunging Density Current , 1999 .

[12]  S. MacIntyre,et al.  Describing fluxes within lakes using temperature arrays and surface meteorology , 2008 .

[13]  Jörg Imberger,et al.  A DYNAMIC RESERVOIR SIMULATION MODEL - DYRESM: 5 , 1981 .

[14]  A. Wüest,et al.  Application of k-ϵ turbulence models to enclosed basins: The role of internal seiches , 2002 .

[15]  H. Stefan,et al.  Plunging Flow into a Reservoir: Theory , 1984 .

[16]  F. Schanz,et al.  The influence of a flood event on phytoplankton succession , 1990, Aquatic Sciences.

[17]  K. Jin,et al.  Application of three-dimensional hydrodynamic model for lake Okeechobee , 2000 .

[18]  J. Downing,et al.  The global abundance and size distribution of lakes, ponds, and impoundments , 2006 .

[19]  M. J. Hall,et al.  Rainfall-Runoff Modelling , 2004 .

[20]  A. Wüest,et al.  Effects of upstream hydropower operation on riverine particle transport and turbidity in downstream lakes , 2006 .

[21]  H. Fischer.  Mixing in Inland and Coastal Waters , 1979 .

[22]  J. Imberger,et al.  Collie River Underflow into the Wellington Reservoir , 1979 .

[23]  S. Monismith,et al.  On the effects of topography on wind and the generation of currents in a large multi-basin lake , 2004, Hydrobiologia.

[24]  K. Ashida,et al.  BASIC STUDY OF TURBIDITY CURRENTS , 1975 .

[25]  L. Hinzman,et al.  Hydrograph separations in an arctic watershed using mixing model and graphical techniques , 1997 .

[26]  H. Fischer,et al.  Observations of transport to surface waters from a plunging inflow to Lake Mead1 , 1983 .

[27]  Sanjoy Banerjee,et al.  The Air-Water Interface: Turbulence and Scalar Exchange , 2007 .

[28]  G. Lawrence,et al.  Natural dispersion in a small lake , 1995 .

[29]  R. Gu,et al.  Two-Dimensional Simulations of Contaminant Currents in Stratified Reservoir , 1998 .

[30]  S. MacIntyre,et al.  Physical pathways of nutrient supply in a small, ultraoligotrophic arctic lake during summer stratification , 2006 .

[31]  Akira Okubo,et al.  Oceanic diffusion diagrams , 1971 .

[32]  Robert L. Haney,et al.  On the Pressure Gradient Force over Steep Topography in Sigma Coordinate Ocean Models , 1991 .

[33]  K. Beven.  Rainfall-Runoff Modelling: The Primer , 2012 .

[34]  J. Imberger,et al.  Modeling basin‐scale internal waves in a stratified lake , 2000 .

[35]  P. Liu,et al.  Advances in Coastal and Ocean Engineering , 1999 .

[36]  S. MacIntyre,et al.  Vertical mixing in a shallow, eutrophic lake: Possible consequences for the light climate of phytoplankton , 1993 .

[37]  J. Lumley,et al.  A First Course in Turbulence , 1972 .

[38]  A. Jassby,et al.  Vertical patterns of eddy diffusion during stratification in Castle Lake, California1 , 1975 .

[39]  S. McCutcheon,et al.  Modeling Reservoir Density Underflow And Interflow From A Chemical Spill , 1996 .

[40]  J. Turner,et al.  Turbulent entrainment in stratified flows , 1959, Journal of Fluid Mechanics.

[41]  Alan F. Blumberg,et al.  Three-Dimensional Hydrothermal Model of Onondaga Lake, New York , 1999 .

[42]  Anthony J. Jakeman,et al.  Ten iterative steps in development and evaluation of environmental models , 2006, Environ. Model. Softw..

[43]  Marc C. Johnson,et al.  An Assessment of Reservoir Density Currents and Inflow Processes , 1983 .

[44]  W. Vincent,et al.  Eutrophication processes regulated by a plunging river inflow , 1991, Hydrobiologia.

[45]  S. Schladow,et al.  Basin-scale internal wave dynamics during a winter cooling period in a large lake , 2003 .

[46]  Herbert E. Huppert,et al.  Entrainment into two-dimensional and axisymmetric turbulent gravity currents , 1996, Journal of Fluid Mechanics.

[47]  J. Imran,et al.  Three-Dimensional Modeling of Negatively Buoyant Flow in Diverging Channels , 2003 .

[48]  J. Barko,et al.  The effects of disturbance events on phytoplankton community structure in a small temperate reservoir , 1999 .

[49]  Robert Hallberg,et al.  Simulation of Density-Driven Frictional Downslope Flow in Z-Coordinate Ocean Models , 1998 .

[50]  H. Stefan,et al.  Negatively buoyant flow in a diverging channel. I: Flow regimes , 1987 .

[51]  E. Cowen,et al.  Residence time of a freshwater embayment connected to a large lake , 2005 .

[52]  C. R. Murthy,et al.  Horizontal Diffusion Characteristics in Lake Ontario , 1976 .

[53]  Dieter M. Imboden,et al.  Horizontal mixing in lakes , 1996 .

[54]  S. Savage,et al.  Analysis Of Plunging Phenomena In Water Reservoirs , 1975 .

[55]  H. Nepf,et al.  The production of chemical heterogeneity in Upper Mystic Lake , 2000 .

[56]  J. Wettlaufer,et al.  The long-term circulation driven by density currents in a two-layer stratified basin , 2007, Journal of Fluid Mechanics.

[57]  R. Döscher,et al.  A Method for Improved Representation of Dense Water Spreading over Topography in Geopotential-Coordinate Models , 1997 .

[58]  D. Durran.  Numerical methods for wave equations in geophysical fluid dynamics , 1999 .

[59]  Vahid Alavian.  Investigation Of Plunging Density Currents Generated By Inflows In Lakes , 2010 .

[60]  J. Imberger,et al.  Coupling an Underflow Model to a Three-Dimensional Hydrodynamic Model , 2003 .

[61]  J. R. Romero,et al.  Spatial‐temporal variability in surface layer deepening and lateral advection in an embayment of Lake Victoria, East Africa , 2002 .

[62]  K. Loague.  Rainfall-Runoff Modelling , 2010 .

[63]  L. Kantha,et al.  An improved mixed layer model for geophysical applications , 1994 .

[64]  J. R. Romero,et al.  Boundary mixing and nutrient fluxes in Mono Lake, California , 1999 .

[65]  David P. Hamilton,et al.  Ten steps applied to development and evaluation of process-based biogeochemical models of estuaries , 2008, Environ. Model. Softw..

[66]  J. Imberger,et al.  Relative Buoyancy Dominates Thermal-Like Flow Interaction along an Incline , 2008 .

[67]  Jonathan P. Fram,et al.  Climate‐related variations in mixing dynamics in an Alaskan arctic lake , 2009 .

[68]  Julie D. Pietrzak,et al.  On the pressure gradient error in sigma coordinate ocean models: A comparison with a laboratory experiment , 1999 .

[69]  G. Kling,et al.  Bacterioplankton Community Shifts in an Arctic Lake Correlate with Seasonal Changes in Organic Matter Source , 2003, Applied and Environmental Microbiology.

[70]  J. Imberger,et al.  Entrainment and turbulence in saline underflow in Lake Ogawara , 2001 .

[71]  P. Killworth,et al.  A Turbulent Bottom Boundary Layer Code for Use in Numerical Ocean Models , 1999 .