Hydrographic changes in the Labrador Sea, 1960–2005

The Labrador Sea has exhibited significant temperature and salinity variations over the past five decades. The whole basin was extremely warm and salty between the mid-1960s and early 1970s, and fresh and cold between the late 1980s and mid-1990s. The full column salinity change observed between these periods is equivalent to mixing a 6 m thick freshwater layer into the water column of the early 1970s. The freshening and cooling trends reversed in 1994 starting a new phase of heat and salt accumulation in the Labrador Sea sustained throughout the subsequent years. It took only a decade for the whole water column to lose most of its excessive freshwater, reinstate stratification and accumulate enough salt and heat to approach its record high salt and heat contents observed between the late 1960s and the early 1970s. If the recent tendencies persist, the basin’s storages of salt and heat will fairly soon, likely by 2008, exceed their historic highs.The main process responsible for the net cooling and freshening of the Labrador Sea between 1987 and 1994 was deep winter convection, which during this period progressively developed to its record depths. It was caused by the recurrence of severe winters during these years and in its turn produced the deepest, densest and most voluminous Labrador Sea Water (LSW_1987–1994) ever observed. The estimated annual production of this water during the period of 1987–1994 is equivalent to the average volume flux of about 4.5 Sv with some individual annual rates exceeding 7.0 Sv. Once winter convection had lost its strength in the winter of 1994–1995, the deep LSW_1987–1994 layer lost “communication” with the mixed layer above, consequently losing its volume, while gaining heat and salt from the intermediate waters outside the Labrador Sea.While the 1000–2000 m layer was steadily becoming warmer and saltier between 1994 and 2005, the upper 1000 m layer experienced another episode of cooling caused by an abrupt increase in the air-sea heat fluxes in the winter of 1999–2000. This change in the atmospheric forcing resulted in fairly intense convective mixing sufficient to produce a new prominent LSW class (LSW₂₀₀₀) penetrating deeper than 1300 m. This layer was steadily sinking or deepening over the years following its production and is presently overlain by even warmer and apparently less dense water mass, implying that LSW₂₀₀₀ is likely to follow the fate of its deeper precursor, LSW_1987–1994. The increasing stratification of the intermediate layer implies intensification in the baroclinic component of the boundary currents around the mid-depth perimeter of the Labrador Sea.The near-bottom waters, originating from the Denmark Strait overflow, exhibit strong interannual variability featuring distinct short-term basin-scale events or pulses of anomalously cold and fresh water, separated by warm and salty overflow modifications. Regardless of their sign these anomalies pass through the abyss of the Labrador Sea, first appearing at the Greenland side and then, about a year later, at the Labrador side and in the central Labrador Basin.The Northeast Atlantic Deep Water (2500–3200 m), originating from the Iceland–Scotland Overflow Water, reached its historically freshest state in the 2000–2001 period and has been steadily becoming saltier since then. It is argued that LSW_1987–1994 significantly contributed to the freshening, density decrease and volume loss experienced by this water mass between the late 1960s and the mid 1990s via the increased entrainment of freshening LSW, the hydrostatic adjustment to expanding LSW, or both.     

Modification of the damping function in the k–ε model to analyse oscillatory boundary layers

A simple relationship has been developed between the wall coordinate y⁺ and Kolmogorov's length scale using direct numerical simulation (DNS) data for a steady boundary layer. This relationship is then utilized to modify two popular versions of low Reynolds number k–ε model. The modified models are used to analyse a transitional oscillatory boundary layer. A detailed comparison has been made by virtue of velocity profile, turbulent kinetic energy, Reynolds stress and wall shear stress with the available DNS data. It is observed that the low Reynolds number models used in the present study can predict the boundary layer properties in an excellent manner.     

A spatial ecosystem and populations dynamics model (SEAPODYM) – Modeling of tuna and tuna-like populations

An enhanced version of the spatial ecosystem and population dynamics model SEAPODYM is presented to describe spatial dynamics of tuna and tuna-like species in the Pacific Ocean at monthly resolution over 1° grid-boxes. The simulations are driven by a bio-physical environment predicted from a coupled ocean physical–biogeochemical model. This new version of SEAPODYM includes expanded definitions of habitat indices, movements, and natural mortality based on empirical evidences. A thermal habitat of tuna species is derived from an individual heat budget model. The feeding habitat is computed according to the accessibility of tuna predator cohorts to different vertically migrating and non-migrating micronekton (mid-trophic) functional groups. The spawning habitat is based on temperature and the coincidence of spawning fish with presence or absence of predators and food for larvae. The successful larval recruitment is linked to spawning stock biomass. Larvae drift with currents, while immature and adult tuna can move of their own volition, in addition to being advected by currents. A food requirement index is computed to adjust locally the natural mortality of cohorts based on food demand and accessibility to available forage components. Together these mechanisms induce bottom-up and top-down effects, and intra- (i.e. between cohorts) and inter-species interactions. The model is now fully operational for running multi-species, multi-fisheries simulations, and the structure of the model allows a validation from multiple data sources. An application with two tuna species showing different biological characteristics, skipjack (Katsuwonus pelamis) and bigeye (Thunnus obesus), is presented to illustrate the capacity of the model to capture many important features of spatial dynamics of these two different tuna species in the Pacific Ocean. The actual validation is presented in a companion paper describing the approach to have a rigorous mathematical parameter optimization [Senina, I., Sibert, J., Lehodey, P., 2008. Parameter estimation for basin-scale ecosystem-linked population models of large pelagic predators: application to skipjack tuna. Progress in Oceanography]. Once this evaluation and parameterization is complete, it may be possible to use the model for management of tuna stocks in the context of climate and ecosystem variability, and to investigate potential changes due to anthropogenic activities including global warming and fisheries pressures and management scenarios.