Numerical study of circulation and temperature-salinity distributions in the Bras d’Or Lakes

The Bras d’Or Lakes (BdOL) are a large, complex and virtually land-locked estuary in central Cape Breton Island of Nova Scotia and one of Canada’s charismatic ecosystems, sustaining ecological and cultural communities unique in many aspects. The BdOL comprise two major basins, many deep and shallow bays, several narrow channels and straits and a large, geologically complex watershed. Predictive knowledge of the water movement within the estuary is a key requirement for effective management and sustainable development of the BdOL ecosystem. A three-dimensional (3D) primitive-equation ocean circulation model is used to examine the estuary’s response to tides, winds and buoyancy forcing associated with freshwater runoff in a series of numerical experiments validated with empirical data. The model results generate intense, jet-like tidal flows of about 1 m s−1 in the channels between the basins and connecting them to the ocean and relatively weak tidal currents in other regions, which agrees well with previous observations and numerical results. Wind forcing and buoyancy forcing associated with river runoff play important roles in generating the significant sub-tidal circulations in the estuary, including narrow channels, deep basins and shallow bays. The circulation model is also used to reconstruct the 3D circulation and temperature-salinity distributions in the summer months of 1974, when current and hydrographic measurements were made at several locations. The sub-tidal circulation in the estuary produced by the model is characterised by wind and barometric set-up and set-down in different sections of the system, and a classic two-layer estuarine circulation in which brackish, near-surface waters flow seaward from the estuary into the Atlantic Ocean, and deep salty waters flow landward through the major channel. The model results reproduce reasonably well the overall features of observed circulation and temperature-salinity fields made in the BdOL in 1974 but generally underestimate the observed currents and density stratification. The model discrepancies reflect the use of spatially mean wind forcing and spatially and monthly mean surface heat flux and the inability of the coarse model horizontal resolution (∼500 m) to resolve narrow channels and straits.