The influence of different turbulence schemes on modelling primary production in a 1D coupled physical–biological model

A one-dimensional (1D) coupled physical–microbiological model has been applied to a site in the central North Sea. The impact of the choice of the turbulence closure scheme on the modelling the primary production has been investigated.The model was run with four different parameterisations of vertical mixing of heat, momentum and dissolved and suspended matters, using M2 tidal forcing and the hourly mean meteorological forcing of 1989 to reproduce the annual thermal structure and primary production. The four mixing parameterisations are: Level 2 turbulence closure scheme [Mellor, G.L., Yamada, T., 1974. A hierarchy of turbulence closure models for planetary boundary layers. J. Atmos. Sci. 31, 1791–1806; Mellor, G.L., Yamada, T., 1982. Development of a turbulence closure model for geophysical Fluid problems. Rev. Geophys. Space Phys. 20 (4) 851–875] using an explicit numerical scheme [Sharples, J., Tett, P., 1994. Modelling the effect of physical variability on the midwater chlorophyll maximum. J. Mar. Res. 52, 219–238]; a version of the Level 2.5 turbulence closure scheme [Galperin, B., Kantha, L.H., Hassid, S., Rosati, A., 1988. A quasi-equilibrium turbulent energy model for geophysical flows. J. Atmos. Sci. 45, 55–62; Ruddick, K.G., Deleersnijder, E., Luyten, P.J., Ozer, J., 1995. Haline stratification in the rhine/meuse freshwater plume: a 3D model sensitivity analysis. Cont. Shelf Res. 15 (13) 1597–1630] simplified to use an algebraic mixing length by Sharples and Simpson [Sharples, J., Simpson, J.H., 1995. Semidiurnal and longer period stability cycles in the Liverpool Bay region of freshwater influence. Cont. Shelf Res. 15, 295–313], also solved explicitly; the same simplified L2.5 scheme with an implicit numerical solution and modified vertical discretisation scheme [Annan, J.D., 1999. Numerical methods for the solution of the turbulence energy equations in the shelf seas. Int. J. Numer. Methods Fluids 29, 193–206]; and another version of the same scheme (but using a different algebraic mixing length) as described by Xing and Davies [Xing, J., Davies, A.M., 1996a. Application of turbulence energy models to the computation of tidal currents and mixing intensities in the shelf edge regions. J. Phys. Oceanogr. 26, 417–447; Xing, J., Davies, A.M., 1996b. Application of a range of turbulence models to the computation of tidal currents and mixing intensities in shelf edge regions. Cont. Shelf. Res. 16, 517–547; Xing, J., Davies, A.M., 1998. Application of a range of turbulence energy models to the computation of the internal tide. Int. J. Numer. Methods Fluids 26, 1055–1084]. Various model outputs at the sea surface and in depth profiles have been compared with data collected in 1989 as part of the North Sea Project [Huthnance, J.M., 1990. Progress on North Sea Project. NERC News, vol. 12, pp. 25–29, UK]. It is shown that the biological results are extremely sensitive to the small changes in the physical conditions, which arise due to the different turbulence schemes tested. The timing of the spring bloom and the maintenance of the midwater chlorophyll maximum all differ greatly between model runs, and the gross primary production varies by a factor of two from the highest to lowest results. The simplified Level 2.5 scheme, implemented using the numerical methods of Annan [Annan, J.D., 1999. Numerical methods for the solution of the turbulence energy equations in the shelf seas. Int. J. Numer. Methods Fluids 29, 193–206], produces results, which give the best agreement with the available data.     

Late summer stratification, internal waves, and turbulence in the Yellow Sea

Microstructure profiling measurements at two locations in the Yellow Sea (a deeper central basin and a local shelf break) were analyzed focusing on tidal and internal-wave induced turbulence near the bottom and in the pycnocline. A classical three-layer density structure consisting of weakly stratified surface and bottom boundary layers and a narrow sharp pycnocline is developed by the end of warm season. Turbulence in the surface layer was not influenced by the tidal forcing but by the diurnal cycle of buoyancy flux and wind forcing at the sea surface. The enhanced dissipation and diffusivity generated by the shear stress at the seafloor was found in the water interior at heights 10–15 m above the bottom with a phase shift of ~ 5–6 m/h. No internal waves, turbulence, or mixing were detected in the pycnocline in the central basin, in contrast to the pycnocline near the local shelf break wherein internal waves of various frequencies were observed all the time. The thickness of the surface layer near the local shelf break slightly exceeded that of the bottom layer (20 vs. 18 m). A 5–6 m high vertical displacement of the pycnocline, which emerged during the low tide, was arguably caused by the passage of an internal soliton of elevation. During this episode, the gradient Richardson number decreased below 0.25 due to enhanced vertical shear, leading to local generation of turbulence with dissipation rates exceeding the background level by an order of magnitude.     

Environmental changes in the western Black Sea related to anthropogenic and natural conditions

The Black Sea northwestern shelf (NWS) ecosystem has been subjected to the strongest anthropogenic pressure of the entire Black Sea as about 80% of the freshwater runoff is discharged there.This paper presents a review of the global environmental changes related mainly to increased eutrophication in the western Black Sea basin.A case study (CoMSBlack 92a cruise) attempts to highlight the interaction of some natural and anthropogenic factors responsible for specific chemical and biological features in the western Black Sea environment. The significance of processes located near the Danube river mouth, which, due to their relatively small space scale, have not been recognized before, is clarified. Hydrological processess of specific importance to the distribution and transformation of antropogenic inputs are river plume dynamics, coastal upwelling and mixing and downwelling over the shelf break and slope in this part of the basin. A layer of low hydrological variability (“conservative layer”) appears to be a natural feature of the area, hypothesized to precondition specific biological and chemical processes there.The results suggest that the interplay between the Danube anthropogenic nutrient load with the natural hydrological fronts and gradients provides opportunities for enhanced biological activity thus contributing to the global environmental changes in the Black Sea NWS.     

One-dimensional modelling of the plankton ecosystem of the north-western Corsican coastal area in relation to meteorological constraints

In order to study the influence of wind mixing on the spring variability of the plankton production of the north western Corsican coastal area, a one-dimensional (1D), vertical, coupled hydrodynamic/biological model (ECOHYDROMV) is used. A hydrodynamic 1D model of the water column with a k–l turbulent closure is applied. The biological model comprises six state variables, representing the plankton ecosystem in the spring period: phytoplankton, copepods, nitrate, ammonium, particulate organic matter of phytoplanktonic origin and particulate organic matter of zooplanktonic origin. The system is influenced by turbulence (expressed by the vertical eddy diffusivity), temperature and irradiance. The model takes into account momentum and heat surface fluxes computed from meteorological data in order to simulate a typical spring atmospheric forcing for the considered area. Results show that primary production vertical structure is characterised by a subsurface maximum which deepens with time and is regulated by the opposite gradients of nitrate concentration and irradiance. Surface plankton productivity is mainly controlled by turbulent vertical transport of nutrients into the mixed layer. The short time scale variability of turbulent mixing generated by the wind appears to be responsible for the plurimodal shape of plankton blooms, observed in the considered area. Furthermore, the model is applied to the study of the spring evolution of the plankton communities off the bay of Calvi (Corsica) for the years 1986 and 1988. In order to initiate and validate the model, time series of hydrological, chemical and biological data have been used. The model reproduces accurately the spring evolution of the phytoplankton biomass measured in situ and illustrates that its strong variability in those years was in close relation to the variability of the wind intensity.     

The structure of dissipation in the western Irish Sea front

We report on an intensive campaign in the summer of 2006 to observe turbulent energy dissipation in the vicinity of a tidal mixing front which separates well mixed and seasonally stratified regimes in the western Irish Sea. The rate of turbulent dissipation ε was observed on a section across the front by a combination of vertical profiles with the FLY dissipation profiler and horizontal profiles by shear sensors mounted on an AUV (Autosub). Mean flow conditions and stratification were obtained from a bed mounted ADCP and a vertical chain of thermistors on a mooring. During an Autosub mission of 60 h, the vehicle, moving at a speed of ~ 1.2 m s^− 1, completed 10 useable frontal crossings between end points which were allowed to move with the mean flow. The results were combined with parallel measurements of the vertical profile of ε which were made using FLY for periods of up to 13 h at positions along the Autosub track. The two data sets, which show a satisfactory degree of consistency, were combined to elucidate the space–time variation of dissipation in the frontal zone. Using harmonic analysis, the spatial structure of dissipation was separated from the strong time dependent signal at the M₄ tidal frequency to yield a picture of the cross-frontal distribution of energy dissipation. A complementary picture of the frontal velocity field was obtained from a moored ADCP and estimates of the mean velocity derived from the thermal wind using the observed density distribution. which indicated the presence of a strong (0.2 m s^− 1) jet-like flow in the high gradient region of the front. Under neap tidal conditions, mean dissipation varied across the section by 3 orders of magnitude exceeding 10^− 2 W m^− 3 near the seabed in the mixed regime and decreasing to 10^− 5 W m^− 3. in the strongly stratified interior regime. The spatial pattern of dissipation is consistent in general form with the predictions of models of tidal mixing and does not reflect any strong influence by the frontal jet.     

On the dense water spreading off the Ross Sea shelf (Southern Ocean)

In this study, current meter and hydrological data obtained during the X Italian Expedition in the Ross Sea (CLIMA Project) are analyzed. Our data show a nice agreement with previous data referring to the water masses present in this area and their dynamics. Here, they are used to further analyze the mixing and deepening processes of Deep Ice Shelf Water (DISW) over the northern shelf break of the Ross Sea. In more detail, our work is focused on the elementary mechanisms that are the most efficient in removing dense water from the shelf: either classical mixing effects or density currents that interact with some topographic irregularity in order to drop to deeper levels, or also the variability of the Antarctic Circumpolar Current (ACC) which, in its meandering, can push the dense water off the shelf, thus interrupting its geostrophic flow. We also discuss in detail the (partial) evidence of dramatic interactions of the dense water with bottom particulate, of geological or biological origin, thus generating impulsive or quasi-steady density-turbidity currents. This complex interaction allows one to consider bottom particular and dense water as a unique self-interacting system. In synthesis, this is a first tentative analysis of the effect of bottom particulate on the dense water dynamics in the Ross Sea.     

Upwelling mechanisms in the northwestern Alboran Sea

From April 1996 to July 1997, a series of hydrographic surveys were carried out in the Northwestern part of the Alboran Sea to investigate the upwelling that is an almost permanent feature in this area. Simultaneously a mooring line was deployed in the north part of the eastern section of the Strait of Gibraltar to monitor the variability of the Atlantic Jet (AJ). Two mechanisms are shown to be relevant for the upwelling dynamic in the region: the southward drifting of the AJ and wind stress. A linear relation between the angle under which the Jet enters the Alboran Sea and the distance from the coastline to the front associated with the Jet has been found. This angle that has been estimated from the low passed time series of current velocity measured by the uppermost instrument of the moored line has been then used to identify the onshore–offshore excursions of the Jet. Both upwelling mechanisms are identified from hydrographic data, because each of them has associated a different type of water mass, and they take place in different locations. Wind-driven upwelling dominates in coastal zones, on the shelf, while upwelling associated with southward drifting of the AJ prevails further offshore. The amount of sub-surface water brought up to the surface by each one is of the same order. However, wind-driven upwelling contributes to the fertilization of this region in a major extent because water upwelled by wind is richer in nutrient concentration.     

A two-dimensional response to a tropical storm on the Gulf of Mexico shelf

Surface current data from drifting buoys and remotely sensed wind data recorded over the continental shelf in the northeastern Gulf of Mexico during the passage of tropical storm Josephine in October 1996 are examined. Drifter data show the existence of a strong surface jet (velocities reaching 1 m s⁻¹) that moves up the west Florida shelf and westward along the Louisiana–Texas shelf, and lasts for nearly 1 week. The coastal jet occurs during an intense synoptic scale wind event where wind speeds reach 15 m s⁻¹. A simple force balance and statistical analysis are performed to assess the role of strong wind forcing. The primary balance shows an Ekman-type current. The role of local acceleration is greatest when winds are directed along bathymetry. A simple two-dimensional strongly forced shelf response model developed from the linear steady-state momentum equations also indicates larger along-shore currents due to both Ekman-type forcing by cross-shore winds and a cross-shore pressure gradient arising from conservation of mass. Model parameters fit empirically are within 15% of theoretical values. The simple model explains 30% and 46% of the variance in the observed along-shore and cross-shore surface currents, respectively.     

白令海潮汐能通量和底边界能耗散

The spatial distribution of the energy flux, bottom boundary layer (BBL) energy dissipation, surface elevation amplitude and current magnitude of the major semidiurnal tidal constituents in the Bering Sea are examined in detail. These distributions are obtained from the results of a three-dimensional numerical simulation model (POM). Compared with observation data from seven stations, the root mean square errors of tidal height are 2.6 cm and 1.2 cm for M₂ and N₂ respectively, and those of phase-lag are 21.8° and 15.8° respectively. The majority of the tidal energy flux off the deep basin is along the shelf edge, although some of this flux crosses the shelf edge, especially in the southeast of the shelf break. The total M₂ energy dissipation in the Bering Sea is 30.43 GW, which is about 10 times of that of N₂ and S₂. The semidiurnal tidal energy enters mainly to the Bering Sea by Samalga Pass, Amukta Pass and Seguam Pass, accounting more than 60% of the total energy entering the Being Sea from the Pacific.     

Numerical simulations of seasonal and interannual variability of the Black Sea thermohaline circulation

The Black Sea general circulation is simulated by a primitive equation model with active free surface. The forcing is seasonally variable and is based on available climatic data. The model reproduces the main features of the Black Sea circulation, including the river discharge effects on the mean sea level and the Bosphorus outflow. Model results show that the simulated sea surface elevation increases in spring over the whole sea, reaching a maximum in the Danube delta area. In the same region, a minimum is observed in winter. The amplitude of the seasonal oscillations (about 8–12 cm over the whole basin) is of the same order of magnitude as the maximum horizontal variations (about 15–18 cm between the coastal areas and the basin interior). This strong signal formed mostly by river discharges, along with the seasonal variability in the other forcing functions and the local dynamics creates a well-pronounced interannual variability. The performance of the model in simulating the seasonal and interannual variability is critically analyzed, with a special attention on the cold intermediate water formation and the circulation in the upper 150 m. The simulations demonstrate that the source of intermediate waters is on the shelf, and that the water mass in the core of cold intermediate layer changes with time as a response to the periodic forcing at sea surface. This type of variability is characterized by pronounced interannual changes, proving that important differences could exist between water mass structure in different years, even when using identical atmospheric forcings each year.     

Primary productivity and export fluxes on the Canadian shelf of the Beaufort Sea: A modelling study

We present a coupled sea ice–ocean-biological (including ice algae) model in the Arctic Ocean. The 1D model was developed and implemented on the Canadian Beaufort Sea shelf to examine the importance of different physical processes in controlling the timing and magnitude of primary production and biogenic particle export over an annual cycle (1987). Our results show that the snow and sea ice cover melt and/or break-up controls the timing of the phytoplankton bloom but primary producers (ice algae and phytoplankton) on the outer shelf are essentially nutrient limited. The total annual primary production (22.7 to 27.7 g-C m^? 2) is thus controlled by nutrient “pre-conditioning” in the previous fall and winter and by the depth of wind mixing that is controlled in part by the supply of fresh water at the end of spring (ice melt or runoff). The spring bloom represents about 40% of the total annual primary production and occurs in a period of the year when sampling is often lacking. Time interpolation of observed values to obtain total annual primary production, as done in many studies, was shown to lead to an underestimation of the actual production. Our simulated ratios of export to primary production vary between 0.42 and 0.44.     

An analysis of the characteristics of chlorophyll in the Sulu Sea

The seasonal and spatial variations of chlorophyll concentrations, Sea Surface Temperature (SST), wind fields and wind-induced Ekman pumping in the Sulu Sea are investigated using a set of new remote sensing measurements from October 1997 to December 2004. The results show the seasonality of chlorophyll, wind fields and SST and reveal the phytoplankton blooming events in the Sulu Sea basin during the northeast monsoon season. In summer, chlorophyll concentrations were relatively low (< 0.2 mg/m³) and distributed uniformly throughout the basin with a narrow belt of high chlorophyll concentrations along the coastal waters, particularly the coasts of Borneo and of the Sulu Archipelago. In winter, chlorophyll concentrations increased (> 0.2 mg/m³) throughout the entire basin, and phytoplankton bloomed southward to the central basin, while chlorophyll concentrations reached high levels (1 mg/m³) in the center of the blooms. One peak was observed during the northeast monsoon season each year. SSTs have significant negative correlations with chlorophyll concentrations; i.e., high and uniformly distributed in summer but lower with an obvious tongue of cold waters southward to the central basin in winter. The seasonal variation of chlorophyll concentrations and SST distribution were associated with the seasonally reversing monsoon. The winter phytoplankton blooming and the tongue of the cold waters were correlated to the vertical upwelling cold and nutrient-rich waters drawn by the northeast wind, with the center of the blooms and the location of cold tongues coinciding with the maximum of the wind speeds and the Ekman pumping velocities.     

Mesozooplankton species distribution in the NW and N Iberian shelf during spring 2004: Relationship with frontal structures

We have analysed the mesozooplankton community structure in the southern Bay of Biscay shelf and its relationship with the hydrographic conditions during spring 2004. According to thermohaline characteristics, we observed two frontal zones of distinct origin along the shelf (around 7° and 3°W), that allowed us to differentiate three different hydrographic domains. The westernmost part of the shelf (WC), defined by the presence of relatively warm and salty water related to the presence of the Iberian Poleward Current (IPC), the easternmost region (EC), characterized by colder and fresher water and subject to the influence of freshwater inputs from the Adour river in the French coast, and a region in the Central Cantabrian Sea (CC), where thermohaline characteristics were intermediate between these two extremes. The mixing layer depth (MLD) regime in these areas was also different: the WC region was characterized by a mixed water column, whereas in the EC region the river discharges produces stratification of the upper meters of the water column (< 10 m); in the CC region, we found a distinct vertical mixing regime that separated coastal (stratification) from shelf (mixed water column) stations, giving rise to a notorious across-shelf front. We found a good match between the aforesaid hydrographic regions and the distribution of mesozooplankton species composition and community assemblages: the Mantel correlation between physical variables and mesozooplankton distribution was highly significant (n = 63, r = 0.70, α < 0.001). In the WC region, the community was dominated by Paracalanus parvus, Oithona helgolandica, Acartia clausi and Clausocalanus pergens, while in the EC region the most dominant species were Noctiluca scintillans, Oncaea media and Temora longicornis. The CC region showed similar composition of copepods than the WC region, but larvaceans (Oikopleura spp. and Fritillaria spp.) were more abundant in the CC region than in the WC region. Within each zone, the relative abundances of the dominant species differed between coastal and shelf locations.     

Assimilation of sea surface temperature in the POL Coastal Ocean Modelling System

A one-dimensional scheme is used to assimilate satellite Sea Surface Temperature data into the Proudman Oceanographic Laboratory Coastal Ocean Modelling System, set up in the Irish Sea with a fine resolution ( 1.8 km). The capabilities of the assimilation scheme are investigated using two different sets of satellite data, of lower and similar resolution to that of the model respectively. Comparison of results with independent data show that assimilation improves the modelled Sea Surface Temperature, but does not address model representation of the temperature vertical structure. It is concluded that for the Irish Sea and at the scales resolved by the model, the assimilation problem cannot be approached in a one-dimensional framework. It is also pointed out that forecast error needs to account explicitly for errors in the representation of the vertical structure of the thermal field.Three-dimensional methods that are suited for coastal systems are then suggested.     

PIV velocity field measurements of flow around a KRISO 3600TEU container ship model

The main purpose of this investigation was to demonstrate a useful application of the particle image velocimetry (PIV) method to analyze the complex flow characteristics around a ship. For a sample illustration, the KRISO 3600TEU container ship model was chosen. The flow structure in the stern and near-wake region of the model has been investigated experimentally in a circulating water channel. Instantaneous velocity fields measured by the PIV velocity field measurement technique have been ensemble-averaged to give details of flow structures such as the spatial distributions of the local mean velocity, vorticity, and turbulent kinetic energy. The free-stream velocity was fixed at U _o = 0.6m/s, and the corresponding Reynolds number based on the length between perpendiculars was about 9.0 × 10⁵. The container ship model shows a complicated three-dimensional flow structure in the stern and near-wake regions. The PIV results clearly revealed the formation of large-scale bilge vortices in the stern region and their effect on the flow in the near-wake. The results shown here provide valuable information for hull form design and the validation of viscous ship flow codes and of turbulence models.     

Tidal current variability in the Central Great Barrier Reef

With tidal data from the literature and the field, a two-dimensional (depth-averaged) numerical model was formulated to simulate the dominant semi-diurnal tidal hydrodynamics of the Central Great Barrier Reef continental shelf. Importantly, the individual mesh dimensions of the numerical scheme were set at approximately 2 × 2 km which was sufficient resolution to incorporate the topography of each reef within the matrix. The model provided a new detailed understanding of the influence of the reef matrix on the tidal currents of the outer shelf. In particular, the model demonstrated that the spatial variability in the tidal current's speed and direction exists down to the scale of the 2 km grid size. The model also demonstrated the significant tidally-induced residual currents that result from the interaction with the complex topography of the reef matrix. The advective effect of these tidal currents would be significant as the tidally-induced residual currents are of similar magnitude to the non-tidal currents of the region. Further, the spatial variability in the modelled and observed tidal currents suggests highly spatially variable advective processes operate within the reef waters.     

Interactions of wind and density driven currents in North Sea ROFIs—a model study

Three different versions of a baroclinic three-dimensional circulation model of the North Sea are used to obtain information on the wind and density interactions in the North Sea ROFIs (Regions Of Freshwater Influence): the standard version with fully prognostic treatment of salinity and temperature is compared to a barotropic model run on the same grid on the one hand and to an also fully prognostic model run on a four times coarser grid on the other hand. In order to gain knowledge on the wind and density interactions, two opposing wind directions are chosen for investigation, namely a time of strong north wind, 21st–28th April 1982, and a time of strong southwest wind, 22nd–24th May 1982. In the April case the effect of the salinity gradients on the border of the ROFIs of Rhine, Weser, Ems and Elbe, i.e. along the continental shore, is shown to lead to a clear enhancement of the mean surface currents. In May this result is partly disguised by the additional effect of the thermocline in the deeper parts of the North Sea, i.e. in the classical shelf sea regime region. Nevertheless, the same pattern of enhanced mean surface currents along the coast is detected and is of the same order of magnitude as in the April case. It is thus concluded that although the circulation in the North Sea is reversed by the wind, the density induced component of the general circulation is modified only slightly.     

Nutrient dynamics and phytoplankton development along an estuary–coastal zone continuum: A model study

This study presents a first attempt to quantify the biogeochemical transformations and fluxes of carbon and nutrients along the entire mixing zone of the shallow, tidally-dominated estuary–coastal zone continuum of the Scheldt (Belgium/The Netherlands). A fully transient, two-dimensional, nested-grid hydrodynamic model of the continuum is coupled to the biogeochemical MIRO model for the coastal zone and the CONTRASTE model for the estuary. Transient model simulations are performed with a high spatial (80–750 m) and temporal (30 min) resolution over a period of one year (January–December 1995). The high temporal resolution allows including the short-term variability triggered by the tides, the freshwater discharge and the wind stress. System scale simulations provide time series of nutrient transformations and fluxes along the entire estuary–coastal zone continuum, as well as highly resolved nutrient inventories for the estuarine and the coastal zone sub-domains. Simulation results reveal that the balance between highly variable estuarine nutrient inputs and physical constrains set by the unsteady residual transport field exert an important control on the magnitude and succession of phytoplankton blooms and the ecosystem structure in the coastal zone. In addition, they suggest that the poorly surveyed estuarine–coastal zone interface plays a central role in the continuum. In this dynamic area, marked spatial concentration gradients develop and episodically lead to a reversal of material fluxes from the coast into the estuary. During distinct episodes of the productive period, euryhaline coastal diatoms intrude far upstream into the saline estuary. This intrusion reduces the estuarine nutrient concentrations and export fluxes, thereby reinforcing the nutrient limitation in the coastal area. As a consequence, the estuarine filter does not operate independently from the processes in the coastal zone. The dynamic interplay between the two ecosystems and the intense process rates operating at their transition, therefore, strongly supports our continuum approach.     

Variability of the western Galician upwelling system (NW Spain) during an intensively sampled annual cycle. An EOF analysis approach

The key features of the western Galician shelf hydrography and dynamics are analyzed on a solid statistical and experimental basis. The results allowed us to gather together information dispersed in previous oceanographic works of the region. Empirical orthogonal functions analysis and a canonical correlation analysis were applied to a high-resolution dataset collected from 47 surveys done on a weekly frequency from May 2001 to May 2002. The main results of these analyses are summarized bellow. Salinity, temperature and the meridional component of the residual current are correlated with the relevant local forcings (the meridional coastal wind component and the continental run-off) and with a remote forcing (the meridional temperature gradient at latitude 37°N). About 80% of the salinity and temperature total variability over the shelf, and 37% of the residual meridional current total variability are explained by two EOFs for each variable. Up to 22% of the temperature total variability and 14% of the residual meridional current total variability is devoted to the set up of cross-shore gradients of the thermohaline properties caused by the wind-induced Ekman transport. Up to 11% and 10%, respectively, is related to the variability of the meridional temperature gradient at the Western Iberian Winter Front. About 30% of the temperature total variability can be explained by the development and erosion of the seasonal thermocline and by the seasonal variability of the thermohaline properties of the central waters. This thermocline presented unexpected low salinity values due to the trapping during spring and summer of the high continental inputs from the River Miño recorded in 2001. The low salinity plumes can be traced on the Galician shelf during almost all the annual cycle; they tend to be extended throughout the entire water column under downwelling conditions and concentrate in the surface layer when upwelling favourable winds blow. Our evidences point to the meridional temperature gradient acting as an important controlling factor of the central waters thermohaline properties and in the development and decay of the Iberian Poleward Current.