Coastal observing and forecasting system for the German Bight – estimates of hydrophysical states - E. V. Stanev; J. Schulz-Stellenfleth; J. Staneva; S. Grayek; J. Seemann; W. Petersen (Ocean Science, 2011)
A coastal observing system for Northern and Arctic Seas (COSYNA) aims at construction of a long-term observatory for the German part of the North Sea, elements of which will be deployed as prototype modules in Arctic coastal waters. At present a coastal prediction system deployed in the area of the German Bight integrates near real-time measurements with numerical models in a pre-operational way and provides continuously state estimates and forecasts of coastal ocean state. The measurement suite contributing to the pre-operational set up includes in situ time series from stationary stations, a High-Frequency (HF) radar system measuring surface currents, a FerryBox system and remote sensing data from satellites. The forecasting suite includes nested 3-D hydrodynamic models running in a data-assimilation mode, which are forced with up-to-date meteorological forecast data. This paper reviews the present status of the system and its recent upgrades focusing on developments in the field of coastal data assimilation. Model supported data analysis and state estimates are illustrated using HF radar and FerryBox observations as examples. A new method combining radial surface current measurements from a single HF radar with a priori information from a hydrodynamic model is presented, which optimally relates tidal ellipses parameters of the 2-D current field and the M2 phase and magnitude of the radials. The method presents a robust and helpful first step towards the implementation of a more sophisticated assimilation system and demonstrates that even using only radials from one station can substantially benefit state estimates for surface currents. Assimilation of FerryBox data based on an optimal interpolation approach using a Kalman filter with a stationary background covariance matrix derived from a preliminary model run which was validated against remote sensing and in situ data demonstrated the capabilities of the pre-operational system. Data assimilation significantly improved the performance of the model with respect to both SST and SSS and demonstrated a good skill not only in the vicinity of the Ferry track, but also over larger model areas. The examples provided in this study are considered as initial steps in establishing new coastal ocean products enhanced by the integrated COSYNA-observations and numerical modelling.

A study of the hydrographic conditions in the Adriatic Sea from numerical modelling and direct observations (2000–2008) - P. Oddo; A. Guarnieri (Ocean Science, 2011)
The inter-annual variability of Adriatic Sea hydrographic characteristics is investigated by means of numerical simulation and direct observation. The period under investigation runs from the beginning of 2000 to the end of 2008. The model used to carry out the simulation is derived from the primitive equation component of the Adriatic Forecasting System (AFS). The model is based on the Princeton Ocean Model (POM) adapted in order to reproduce the features of the Adriatic. Both numerical findings and observations agree in depicting a strong inter-annual variability in the entire Adriatic Sea and its sub-basins. Nevertheless, two model deficiencies are identified: an excessive vertical/horizontal mixing and an inaccurate representation of the thermohaline properties of the entering Mediterranean Waters. The dense water formation process has been found to be intermittent. In addition to inter-annual variability, a long-scale signal has been observed in the salinity content of the basin as a consequence of a prolonged period of reduced Po river runoff and high evaporation rates. As a result, the temperature and salinity of the northern Adriatic dense water vary considerably between the beginning and the end of the period investigated.

Silicon pool dynamics and biogenic silica export in the Southern Ocean inferred from Si-isotopes - F. Fripiat; A.-J. Cavagna; F. Dehairs; S. Speich; L. André; D. Cardinal (Ocean Science, 2011)
Silicon isotopic signatures (δ30Si) of water column silicic acid (Si(OH)4) were measured in the Southern Ocean, along a meridional transect from South Africa (Subtropical Zone) down to 57° S (northern Weddell Gyre). This provides the first reported data of a summer transect across the whole Antarctic Circumpolar Current (ACC). δ30Si variations are large in the upper 1000 m, reflecting the effect of the silica pump superimposed upon meridional water transfer across the ACC: the transport of Antarctic surface waters northward by a net Ekman drift and their convergence and mixing with warmer upper-ocean Si-depleted waters to the north. Using Si isotopic signatures, we determine different mixing interfaces: the Antarctic Surface Water (AASW), the Antarctic Intermediate Water (AAIW), and thermoclines in the low latitude areas. The residual silicic acid concentrations of end-members control the δ30Si alteration of the mixing products and with the exception of AASW, all mixing interfaces have a highly Si-depleted mixed layer end-member. These processes deplete the silicic acid AASW concentration northward, across the different interfaces, without significantly changing the AASW δ30Si composition. By comparing our new results with a previous study in the Australian sector we show that during the circumpolar transport of the ACC eastward, the δ30Si composition of the silicic acid pools is getting slightly, but significantly lighter from the Atlantic to the Australian sectors. This results either from the dissolution of biogenic silica in the deeper layers and/or from an isopycnal mixing with the deep water masses in the different oceanic basins: North Atlantic Deep Water in the Atlantic, and Indian Ocean deep water in the Indo-Australian sector. This isotopic trend is further transmitted to the subsurface waters, representing mixing interfaces between the surface and deeper layers. Through the use of δ30Si constraints, net biogenic silica production (representative of annual export), at the Greenwich Meridian is estimated to be 5.2 ± 1.3 and 1.1 ± 0.3 mol Si m−2 for the Antarctic Zone and Polar Front Zone, respectively. This is in good agreement with previous estimations. Furthermore, summertime Si-supply into the mixed layer of both zones, via vertical mixing, is estimated to be 1.6 ± 0.4 and 0.1 ± 0.5 mol Si m−2, respectively.

Coastal sea responses to atmospheric forcings at two different resolutions - Z. B. Klaić; Z. Pasarić; G. Beg Paklar; P. Oddo (Ocean Science, 2011)
We investigated coastal sea responses to three, multi-day strong wind episodes that occurred in the middle Adriatic during the Target Operational Period (TOP) of the European COastal sea OPerational observing and forecasting system (ECOOP) project. A high-resolution oceanographic model (1 km horizontal, 16 σ vertical layers) based on the modified Princeton Ocean Model (POM) was applied to a highly complex domain located in the coastal area of the eastern Adriatic Sea. The oceanographic model was nested into the Adriatic REGional model (AREG-2) covering the entire Adriatic Sea. Meteorological forcing was prepared by two atmospheric models. The coarser model was the European Centre for Medium-range Weather Forecast model (ECMWF, with horizontal and temporal resolutions of 0.25° and 6 h, respectively), and the finer one was the Aire Limitée Adaptation dynamique Développement InterNational model (ALADIN, with horizontal and temporal resolutions of 8 km and 3 h, respectively, and winds dynamically adapted to a horizontal resolution of 2 km). The results show that small-scale atmospheric features, which arise due to the orographically complex mainland and the number of islands and were not reproduced by the coarser atmospheric model, substantially affected surface currents, mass transports, sea surface temperature (SST) and surface salinity in the coastal area during strong Bora. For strong Sirocco, the atmospheric model's resolution was important for currents on the lee sides of islands.

Numerical simulation and decomposition of kinetic energy in the Central Mediterranean: insight on mesoscale circulation and energy conversion - R. Sorgente; A. Olita; P. Oddo; L. Fazioli; A. Ribotti (Ocean Science, 2011)
The spatial and temporal variability of eddy and mean kinetic energy of the Central Mediterranean region has been investigated, from January 2008 to December 2010, by mean of a numerical simulation mainly to quantify the mesoscale dynamics and their relationships with physical forcing. In order to understand the energy redistribution processes, the baroclinic energy conversion has been analysed, suggesting hypotheses about the drivers of the mesoscale activity in this area. The ocean model used is based on the Princeton Ocean Model implemented at 1/32° horizontal resolution. Surface momentum and buoyancy fluxes are interactively computed by mean of standard bulk formulae using predicted model Sea Surface Temperature and atmospheric variables provided by the European Centre for Medium Range Weather Forecast operational analyses. At its lateral boundaries the model is one-way nested within the Mediterranean Forecasting System operational products. The model domain has been subdivided in four sub-regions: Sardinia channel and southern Tyrrhenian Sea, Sicily channel, eastern Tunisian shelf and Libyan Sea. Temporal evolution of eddy and mean kinetic energy has been analysed, on each of the four sub-regions, showing different behaviours. On annual scales and within the first 5 m depth, the eddy kinetic energy represents approximately the 60 % of the total kinetic energy over the whole domain, confirming the strong mesoscale nature of the surface current flows in this area. The analyses show that the model well reproduces the path and the temporal behaviour of the main known sub-basin circulation features. New mesoscale structures have been also identified, from numerical results and direct observations, for the first time as the Pantelleria Vortex and the Medina Gyre. The classical kinetic energy decomposition (eddy and mean) allowed to depict and to quantify the permanent and fluctuating parts of the circulation in the region, and to differentiate the four sub-regions as function of relative and absolute strength of the mesoscale activity. Furthermore the Baroclinic Energy Conversion term shows that in the Sardinia Channel the mesoscale activity, due to baroclinic instabilities, is significantly larger than in the other sub-regions, while a negative sign of the energy conversion, meaning a transfer of energy from the Eddy Kinetic Energy to the Eddy Available Potential Energy, has been recorded only for the surface layers of the Sicily Channel during summer.

Tidal generation of large sub-mesoscale eddy dipoles - W. Callendar; J. M. Klymak; M. G. G. Foreman (Ocean Science, 2011)
Numerical simulations of tidal flow past Cape St. James on the south tip of Haida Gwaii (Queen Charlotte Islands) are presented that indicate mesoscale dipoles are formed from coalescing tidal eddies. Observations in this region demonstrate robust eddy generation at the Cape, with the primary process being flow separation of buoyant or wind driven outflows forming large anti-cyclonic, negative potential vorticity, Haida Eddies. However, there are other times where dipoles are observed in satellites, indicating a source of positive potential vorticity must also be present. The simulations here build on previous work that implicates oscillating tidal flow past the cape in creating the positive vorticity. Small headland eddies of alternating vorticity are created each tide. During certain tidal cycles, the headland eddies coalesce and self organize in such a way as to create large >20-km diameter eddies that then self-advect into deep water. The self advection speed is faster than the beta drift of anti-cyclones, and the propagation direction appears to be more southerly than typical Haida Eddies, though the model contains no mean wind-driven flows. These eddies are smaller than Haida Eddies, but given their tidal origin, may represent a more consistent source of coastal water that is injected into the interior of the subpolar gyre.

High frequency variability of the Atlantic meridional overturning circulation - B. Balan Sarojini; J. M. Gregory; R. Tailleux; G. R. Bigg; A. T. Blaker; D. R. Cameron; N. R. Edwards; A. P. Megann; L. C. Shaffrey; B. Sinha (Ocean Science, 2011)
We compare the variability of the Atlantic meridional overturning circulation (AMOC) as simulated by the coupled climate models of the RAPID project, which cover a wide range of resolution and complexity, and observed by the RAPID/MOCHA array at about 26° N. We analyse variability on a range of timescales, from five-daily to interannual. In models of all resolutions there is substantial variability on timescales of a few days; in most AOGCMs the amplitude of the variability is of somewhat larger magnitude than that observed by the RAPID array, while the time-mean is within about 10 % of the observational estimate. The amplitude of the simulated annual cycle is similar to observations, but the shape of the annual cycle shows a spread among the models. A dynamical decomposition shows that in the models, as in observations, the AMOC is predominantly geostrophic (driven by pressure and sea-level gradients), with both geostrophic and Ekman contributions to variability, the latter being exaggerated and the former underrepresented in models. Other ageostrophic terms, neglected in the observational estimate, are small but not negligible. The time-mean of the western boundary current near the latitude of the RAPID/MOCHA array has a much wider model spread than the AMOC does, indicating large differences among models in the simulation of the wind-driven gyre circulation, and its variability is unrealistically small in the models. In many RAPID models and in models of the Coupled Model Intercomparison Project Phase 3 (CMIP3), interannual variability of the maximum of the AMOC wherever it lies, which is a commonly used model index, is similar to interannual variability in the AMOC at 26° N. Annual volume and heat transport timeseries at the same latitude are well-correlated within 15–45° N, indicating the climatic importance of the AMOC. In the RAPID and CMIP3 models, we show that the AMOC is correlated over considerable distances in latitude, but not the whole extent of the North Atlantic; consequently interannual variability of the AMOC at 50° N, where it is particularly relevant to European climate, is not well-correlated with that of the AMOC at 26° N, where it is monitored by the RAPID/MOCHA array.

Modelling the seasonal variability of the Antarctic Slope Current - P. Mathiot; H. Goosse; T. Fichefet; B. Barnier; H. Gallée (Ocean Science, 2011)
One of the main features of the oceanic circulation along Antarctica is the Antarctic Slope Current (ASC). This circumpolar current flows westwards and contributes to communication between the three major oceanic basins around Antarctica. The ASC is not very well known due to remote location and the presence of sea ice during several months, allowing in situ studies only during summertime. Moreover, only few modelling studies of this current have been carried out. Here, we investigate the sensitivity of this simulated current to four different resolutions in a coupled ocean-sea ice model and to two different atmospheric forcing sets. Two series of simulations are conducted. For the first series, global model configurations are run at coarse (2°) to eddy-permitting (0.25°) resolutions with the same atmospheric forcing. For the second series, simulations with two different atmospheric forcings are performed using a regional circumpolar configuration (south of 30° S) at 0.5° resolution. The first atmospheric forcing is based on a global atmospheric reanalysis and satellite data, while the second is based on a downscaling of the global atmospheric reanalysis by a regional atmospheric model calibrated to Antarctic meteorological conditions. Sensitivity experiments to resolution indicate that a minimum model resolution of 0.5° is needed to capture the dynamics of the ASC in terms of water mass transport and recirculation. Sensitivity experiments to atmospheric forcing fields shows that the wind speed along the Antarctic coast strongly controls the water mass transport and the seasonal cycle of the ASC. An increase in annual mean of easterlies by about 30 % leads to an increase in the mean ASC transport by about 40 %. Similar effects are obtained on the seasonal cycle: using a wind forcing field with a larger seasonal cycle (+30 %) increases by more than 30 % the amplitude of the seasonal cycle of the ASC. To confirm the importance of wind seasonal cycle, a simulation without wind speed seasonal cycle is carried out. This simulation shows a decrease by more than 50 % of the amplitude of the ASC transport seasonal cycle without changing the mean value of ASC transport.

An ECOOP web portal for visualising and comparing distributed coastal oceanography model and in situ data - A. L. Gemmell; R. M. Barciela; J. D. Blower; K. Haines; Q. Harpham; K. Millard; M. R. Price; A. Saulter (Ocean Science, 2011)
As part of a large European coastal operational oceanography project (ECOOP), we have developed a web portal for the display and comparison of model and in situ marine data. The distributed model and in situ datasets are accessed via an Open Geospatial Consortium Web Map Service (WMS) and Web Feature Service (WFS) respectively. These services were developed independently and readily integrated for the purposes of the ECOOP project, illustrating the ease of interoperability resulting from adherence to international standards. The key feature of the portal is the ability to display co-plotted timeseries of the in situ and model data and the quantification of misfits between the two. By using standards-based web technology we allow the user to quickly and easily explore over twenty model data feeds and compare these with dozens of in situ data feeds without being concerned with the low level details of differing file formats or the physical location of the data. Scientific and operational benefits to this work include model validation, quality control of observations, data assimilation and decision support in near real time. In these areas it is essential to be able to bring different data streams together from often disparate locations.

A new tide model for the Mediterranean Sea based on altimetry and tide gauge assimilation - D. N. Arabelos; D. Z. Papazachariou; M. E. Contadakis; S. D. Spatalas (Ocean Science, 2011)
The tides for the Mediterranean Sea are described through a high resolution model (MEDI10) developed by assimilation of tide-gauge data and T/P data into a barotropic ocean tide model. Tidal parameters from 56 coastal tide-gauge stations around the Mediterranean for eight principal constituents: M2, S2, N2, K2, K1, O1, P1 and Q1 and from 20 stations for M2, S2, K1, O1 are included in the model. TOPEX/Poseidon data with all corrections applied except for the ocean tides and bathymetry from TOPO 13.1 were used for development of the model. Numerical experiments were carried out for the estimation of the friction velocity and of the decorrelation length scale. The experiments related to the friction velocity showed that the use of spatially varying friction velocity, estimated as a function of position in the model domain, gives better results than a constant value. The experiments related to the estimation of the decorrelation length suggest that the results are not sensitive for lengths close to ten times the length of the grid cell. The assessment of the model is based on ten tide-gauge observations that are not used for the assimilation. Comparisons were carried out with contemporary published global or regional models. The final solution is computed using 76 selected coastal tide-gauge stations. The comparison between the observed and the model constituents results in a Root Sum of Squares (RSS) equal to 1.3 cm.

About the seasonal and fortnightly variabilities of the Mediterranean outflow - C. Millot; J. Garcia-Lafuente (Ocean Science, 2011)
CTD time series from the HYDRO-CHANGES programme and INGRES projects have been collected simultaneously (2004–2008) on the shelf of Morocco and at the sills of Camarinal and Espartel in the strait of Gibraltar. They provide information that supports results recently obtained from the analysis of the two former time series, as well as from a reanalysis of GIBEX CTD profiles (1985–1986). The outflow of Mediterranean Waters, which does not show a clear seasonal variability before entering the strait, strongly mixes within the strait, due mainly to the internal tide, with the seasonally variable inflow of Atlantic Water. The outflow thus gets marked seasonal and fortnightly variabilities within the strait. Furthermore, since the outflowing waters entering the strait display marked spatial heterogeneity and long-term temporal variabilities, accurately predicting the characteristics of the Mediterranean outflow into the North Atlantic Ocean appears almost impossible.

Sensitivity analysis of an ocean carbon cycle model in the North Atlantic: an investigation of parameters affecting the air-sea CO 2 flux, primary production and export of detritus - V. Scott; H. Kettle; C. J. Merchant (Ocean Science, 2011)
The sensitivity of the biological parameters in a nutrient-phytoplankton-zooplankton-detritus (NPZD) model in the calculation of the air-sea CO2 flux, primary production and detrital export is analysed. We explore the effect on these outputs of variation in the values of the twenty parameters that control ocean ecosystem growth in a 1-D formulation of the UK Met Office HadOCC NPZD model used in GCMs. We use and compare the results from one-at-a-time and all-at-a-time perturbations performed at three sites in the EuroSITES European Ocean Observatory Network: the Central Irminger Sea (60° N 40° W), the Porcupine Abyssal Plain (49° N 16° W) and the European Station for Time series in the Ocean Canary Islands (29° N 15° W). Reasonable changes to the values of key parameters are shown to have a large effect on the calculation of the air-sea CO2 flux, primary production, and export of biological detritus to the deep ocean. Changes in the values of key parameters have a greater effect in more productive regions than in less productive areas. The most sensitive parameters are generally found to be those controlling well-established ocean ecosystem parameterisations widely used in many NPZD-type models. The air-sea CO2 flux is most influenced by variation in the parameters that control phytoplankton growth, detrital sinking and carbonate production by phytoplankton (the rain ratio). Primary production is most sensitive to the parameters that define the shape of the photosynthesis-irradiance curve. Export production is most sensitive to the parameters that control the rate of detrital sinking and the remineralisation of detritus.

North Atlantic 20th century multidecadal variability in coupled climate models: sea surface temperature and ocean overturning circulation - I. Medhaug; T. Furevik (Ocean Science, 2011)
Output from a total of 24 state-of-the-art Atmosphere-Ocean General Circulation Models is analyzed. The models were integrated with observed forcing for the period 1850–2000 as part of the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report. All models show enhanced variability at multi-decadal time scales in the North Atlantic sector similar to the observations, but with a large intermodel spread in amplitudes and frequencies for both the Atlantic Multidecadal Oscillation (AMO) and the Atlantic Meridional Overturning Circulation (AMOC). The models, in general, are able to reproduce the observed geographical patterns of warm and cold episodes, but not the phasing such as the early warming (1930s–1950s) and the following colder period (1960s–1980s). This indicates that the observed 20th century extreme in temperatures are due to primarily a fortuitous phasing of intrinsic climate variability and not dominated by external forcing. Most models show a realistic structure in the overturning circulation, where more than half of the available models have a mean overturning transport within the observed estimated range of 13–24 Sverdrup. Associated with a stronger than normal AMOC, the surface temperature is increased and the sea ice extent slightly reduced in the North Atlantic. Individual models show potential for decadal prediction based on the relationship between the AMO and AMOC, but the models strongly disagree both in phasing and strength of the covariability. This makes it difficult to identify common mechanisms and to assess the applicability for predictions.

The effects of biogeochemical processes on oceanic conductivity/salinity/density relationships and the characterization of real seawater - R. Pawlowicz; D. G. Wright; F. J. Millero (Ocean Science, 2011)
As seawater circulates through the global ocean, its relative composition undergoes small variations. This results in changes to the conductivity/salinity/density relationship, which is currently well-defined only for Standard Seawater obtained from a particular area in the North Atlantic. These changes are investigated here by analysis of laboratory experiments in which salts are added to seawater, by analysis of oceanic observations of density and composition anomalies, and by mathematical investigation using a model relating composition, conductivity, and density of arbitrary seawaters. Mathematical analysis shows that understanding and describing the effect of changes in relative composition on operational estimates of salinity using the Practical Salinity Scale 1978 and on density using an equation of state for Standard Seawater require the use of a number of different salinity variables and a family of haline contraction coefficients. These salinity variables include an absolute Salinity SAsoln, a density salinity SAdens, the reference salinity SR, and an added-mass salinity SAadd. In addition, a new salinity variable S∗ is defined, which represents the preformed salinity of a Standard Seawater component of real seawater to which biogeochemical processes add material. In spite of this complexity, observed correlations between different ocean biogeochemical processes allow the creation of simple formulas that can be used to convert between the different salinity and density measures, allowing for the operational reduction of routine oceanographic observations.

Comparison of SeaWiFS and MODIS time series of inherent optical properties for the Adriatic Sea - F. Mélin (Ocean Science, 2011)
Time series of inherent optical properties (IOPs) derived from SeaWiFS and MODIS are compared for the Adriatic Sea. The IOPs are outputs of the Quasi-Analytical Algorithm and include total absorption a, phytoplankton absorption aph, absorption associated with colored detrital material (CDM) acdm, and particle backscattering coefficient bbp. The average root-mean square difference Δ computed for log-transformed distributions decreases for $a$ from 0.084 at 412 nm to 0.052 at 490 nm, is higher for aph(443) (0.149) than for acdm(443) (0.071), and is approximately 0.165 for bbp at various wavelengths. The SeaWiFS a at 443 and 490 nm, aph at 443 nm and bbp are on average higher than the MODIS counterparts. Statistics show significant variations in space and time. There is an overall increasing gradient for Δ associated with the absorption terms from the open southern and central Adriatic to the northwest part of the basin, and a reversed gradient for the particulate backscattering coefficient. For time series analysis, only a(412) and acdm(443) currently present an unbiased continuity bridging the SeaWiFS and MODIS periods for the Adriatic Sea.

Mixing, heat fluxes and heat content evolution of the Arctic Ocean mixed layer - A. Sirevaag; S. de la Rosa; I. Fer; M. Nicolaus; M. Tjernström; M. G. McPhee (Ocean Science, 2011)
A comprehensive measurement program was conducted during 16 days of a 3 week long ice pack drift, from 15 August to 1 September 2008 in the central Amundsen Basin, Arctic Ocean. The data, sampled as part of the Arctic Summer Cloud Ocean Study (ASCOS), included upper ocean stratification, mixing and heat transfer as well as transmittance solar radiation through the ice. The observations give insight into the evolution of the upper layers of the Arctic Ocean in the transition period from melting to freezing. The ocean mixed layer was found to be heated from above and, for summer conditions, the net heat flux through the ice accounted for 22 % of the observed change in mixed layer heat content. Heat was mixed downward within the mixed layer and a small, downward heat flux across the base of the mixed layer accounted for the accumulated heat in the upper cold halocline during the melting season. On average, the ocean mixed layer was cooled by an ocean heat flux at the ice/ocean interface (1.2 W m−2) and heated by solar radiation through the ice (−2.6 W m−2). An abrupt change in surface conditions halfway into the drift due to freezing and snowfall showed distinct signatures in the data set and allowed for inferences and comparisons to be made for cases of contrasting forcing conditions. Transmittance of solar radiation was reduced by 59 % in the latter period. From hydrographic observations obtained earlier in the melting season, in the same region, we infer a total fresh water equivalent of 3.3 m accumulated in the upper ocean, which together with the observed saltier winter mixed layer indicates a transition towards a more seasonal ice cover in the Arctic.

Comparison between three implementations of automatic identification algorithms for the quantification and characterization of mesoscale eddies in the South Atlantic Ocean - J. M. A. C. Souza; C. de Boyer Montégut; P. Y. Le Traon (Ocean Science, 2011)
Three methods for automatic detection of mesoscale coherent structures are applied to Sea Level Anomaly (SLA) fields in the South Atlantic. The first method is based on the wavelet packet decomposition of the SLA data, the second on the estimation of the Okubo-Weiss parameter and the third on a geometric criterion using the winding-angle approach. The results provide a comprehensive picture of the mesoscale eddies over the South Atlantic Ocean, emphasizing their main characteristics: amplitude, diameter, duration and propagation velocity. Five areas of particular eddy dynamics were selected: the Brazil Current, the Agulhas eddies propagation corridor, the Agulhas Current retroflexion, the Brazil-Malvinas confluence zone and the northern branch of the Antarctic Circumpolar Current (ACC). For these areas, mean propagation velocities and amplitudes were calculated. Two regions with long duration eddies were observed, corresponding to the propagation of Agulhas and ACC eddies. Through the comparison between the identification methods, their main advantages and shortcomings were detailed. The geometric criterion presents the best performance, mainly in terms of number of detections, duration of the eddies and propagation velocities. The results are particularly good for the Agulhas Rings, which have the longest lifetimes of all South Atlantic eddies.

On the freshening of the northwestern Weddell Sea continental shelf - H. H. Hellmer; O. Huhn; D. Gomis; R. Timmermann (Ocean Science, 2011)
We analyzed hydrographic data from the northwestern Weddell Sea continental shelf of the three austral winters 1989, 1997, and 2006 and two summers following the last winter cruise. During summer a thermal front exists at ~64° S separating cold southern waters from warm northern waters that have similar characteristics as the deep waters of the central basin of the Bransfield Strait. In winter, the whole continental shelf exhibits southern characteristics with high Neon (Ne) concentrations, indicating a significant input of glacial melt water. The comparison of the winter data from the shallow shelf off the tip of the Antarctic Peninsula, spanning a period of 17 yr, shows a salinity decrease of 0.09 for the whole water column, which has a residence time of

Flow and mixing near a glacier tongue: a pilot study - C. L. Stevens; C. L. Stewart; N. J. Robinson; M. J. M. Williams; T. G. Haskell (Ocean Science, 2011)
A glacier tongue floating in the coastal ocean presents a significant obstacle to the local flow and so influences oceanic mixing and transport processes. Here acoustic Doppler current profiler and shear microstructure observations very near to a glacier tongue side-wall capture flow accelerations and associated mixing. Flow speeds reached around 40 cm s−1, twice that of the ambient tidal flow amplitude, and generated vertical velocity shear squared as large as 10−5 s−2. During the time of maximum flow, turbulent energy dissipation rates reached 10−5 m2 s−3, around three decades greater than local background levels. This is in keeping with estimates of the gradient Richardson Number which dropped to ~1 during maximum flow. Associated vertical diffusivities estimated from the shear microstructure results were substantial, reflecting the influence of the glacier on velocity gradients.

Sensitivity study of the generation of mesoscale eddies in a numerical model of Hawaii islands - M. Kersalé; A. M. Doglioli; A. A. Petrenko (Ocean Science, 2011)
The oceanic circulation around the Hawaiian archipelago is characterized by a complex circulation and the presence of mesoscale eddies west of the islands. These eddies typically develop and persist for weeks to several months in the area during persistent trade winds conditions. A series of numerical simulations on the Hawaiian region has been done in order to examine the relative importance of wind, inflow current and topographic forcing on the general circulation and the generation of eddies. Moreover, numerical cyclonic eddies are compared with the one observed during the cruise E-FLUX (Dickey et al., 2008). Our study demonstrates the need for all three forcings (wind, inflow current and topography) to reproduce the known oceanic circulation. In particular, the cumulative effect plays a key role on the generation of mesoscale eddies. The wind-stress-curl, via the Ekman pumping mechanism, has also been identified as an important mechanism upon the strength of the upwelling in the lee of the Big Island of Hawaii. In order to find the best setup of a regional ocean model, we compare more precisely numerical results obtained using two different wind databases: COADS and QuikSCAT. The main features of the ocean circulation in the area are well reproduced by our model forced by both COADS and QuickSCAT climatologies. Nevertheless, significant differences appear in the levels of kinetic energy and vorticity. The wind-forcing spatial resolution clearly affects the way in which the wind momentum feeds the mesoscale phenomena. The higher the resolution, the more realistic the ocean circulation. In particular, the simulation forced by QuikSCAT wind data reproduces well the observed energetic mesoscale structures and their hydrological characteristics and behaviors.

Numerical modelling of thermodynamics and dynamics of sea ice in the Baltic Sea - A. Herman; J. Jedrasik; M. Kowalewski (Ocean Science, 2011)
In this paper, a numerical dynamic-thermo-dynamic sea-ice model for the Baltic Sea is used to analyze the variability of ice conditions in three winter seasons. The modelling results are validated with station (water temperature) and satellite data (ice concentration) as well as by qualitative comparisons with the Swedish Meteorological and Hydrological Institute ice charts. Analysis of the results addresses two major questions. One concerns effects of meteorological forcing on the spatio-temporal distribution of ice concentration in the Baltic. Patterns of correlations between air temperature, wind speed, and ice-covered area are demonstrated to be different in larger, more open sub-basins (e.g., the Bothnian Sea) than in the smaller ones (e.g., the Bothnian Bay). Whereas the correlations with the air temperature are positive in both cases, the influence of wind is pronounced only in large basins, leading to increase/decrease of areas with small/large ice concentrations, respectively. The other question concerns the role of ice dynamics in the evolution of the ice cover. By means of simulations with the dynamic model turned on and off, the ice dynamics is shown to play a crucial role in interactions between the ice and the upper layers of the water column, especially during periods with highly varying wind speeds and directions. In particular, due to the fragmentation of the ice cover and the modified surface fluxes, the ice dynamics influences the rate of change of the total ice volume, in some cases by as much as 1 km3 per day. As opposed to most other numerical studies on the sea-ice in the Baltic Sea, this work concentrates on the short-term variability of the ice cover and its response to the synoptic-scale forcing.

Estimation of friction parameters in gravity currents by data assimilation in a model hierarchy - A. Wirth (Ocean Science, 2011)
This paper is the last in a series of three investigating the friction laws and their parametrisation in idealised gravity currents in a rotating frame. Results on the dynamics of a gravity current (Wirth, 2009) and on the estimation of friction laws by data assimilation (Wirth and Verron, 2008) are combined to estimate the friction parameters and discriminate between friction laws in non-hydrostatic numerical simulations of gravity current dynamics, using data assimilation and a reduced gravity shallow water model. I demonstrate, that friction parameters and laws in gravity currents can be estimated using data assimilation. The results clearly show that friction follows a linear Rayleigh law for small Reynolds numbers and the estimated value agrees well with the analytical value obtained for non-accelerating Ekman layers. A significant and sudden departure towards a quadratic drag law at an Ekman layer based Reynolds number of around 800 is shown, in agreement with classical laboratory experiments. The drag coefficient obtained compares well to friction values over smooth surfaces. I show that data assimilation can be used to determine friction parameters and discriminate between friction laws and that it is a powerful tool in systematically connecting models within a model hierarchy.

Comparison of the fall rate and structure of recent T-7 XBT manufactured by Sippican and TSK - S. Kizu; C. Sukigara; K. Hanawa (Ocean Science, 2011)
The fall rate of recent T-7 expendable bathythermograph (XBT; 760 m) is evaluated based on a series of concurrent measurement with a calibrated Conductivity Temperature Depth profiler (CTD) in the sea east of Japan. An emphasis is placed on comparing the fall rates of T-7 produced by the two present manufacturers, the Lockheed Martin Sippican Inc., and the Tsurumi Seiki Co. Ltd., which have been believed to be identical but had never been compared directly. It is found that the two manufacturers' T-7 fall at rates different by about 3.5%. The Sippican T-7 falls slower than given by the fall-rate equation (FRE) of Hanawa et al. (1995) by about 2.1%, and the TSK T-7 falls faster by about 1.4%. The fall-rate coefficients estimated based on the sea test by applying the equation of traditional quadratic form, d(t)=at−bt2 where d is depth in meters and t is the time elapsed, in seconds, are a=6.553 (m s−1) and b=0.00221 (m s−2) for the LMS T-7, and a=6.803 (m s−1) and b=0.00242 (m s−2) for the TSK T-7. By detail examination of the probes, we found that the two companies' T-7 have different total weight and many structural differences. Because the difference in the fall rate is about twice larger than the difference in weight (about 2%), it is inferred that the structural differences give sizable impact to the difference in their fall rates. Our results clearly show that the recent T-7 of the two companies needs to be discriminated.

Multifractal analysis of oceanic chlorophyll maps remotely sensed from space - L. de Montera; M. Jouini; S. Verrier; S. Thiria; M. Crepon (Ocean Science, 2011)
Phytoplankton patchiness has been investigated with multifractal analysis techniques. We analyzed oceanic chlorophyll maps, measured by the SeaWiFS orbiting sensor, which are considered to be good proxies for phytoplankton. The study area is the Senegalo-Mauritanian upwelling region, because it has a low cloud cover and high chlorophyll concentrations. Multifractal properties are observed, from the sub-mesoscale up to the mesoscale, and are found to be consistent with the Corssin-Obukhov scale law of passive scalars. This result indicates that, in this specific region and within this scale range, turbulent mixing would be the dominant effect leading to the observed variability of phytoplankton fields. Finally, it is shown that multifractal patchiness can be responsible for significant biases in the nonlinear source and sink terms involved in biogeochemical numerical models.

The effect of tides on dense water formation in Arctic shelf seas - C. F. Postlethwaite; M. A. Morales Maqueda; V. le Fouest; G. R. Tattersall; J. Holt; A. J. Willmott (Ocean Science, 2011)
Ocean tides are not explicitly included in many ocean general circulation models, which will therefore omit any interactions between tides and the cryosphere. We present model simulations of the wind and buoyancy driven circulation and tides of the Barents and Kara Seas, using a 25 km × 25 km 3-D ocean circulation model coupled to a dynamic and thermodynamic sea ice model. The modeled tidal amplitudes are compared with tide gauge data and sea ice extent is compared with satellite data. Including tides in the model is found to have little impact on overall sea ice extent but is found to delay freeze up and hasten the onset of melting in tidally active coastal regions. The impact that including tides in the model has on the salt budget is investigated and found to be regionally dependent. The vertically integrated salt budget is dominated by lateral advection. This increases significantly when tides are included in the model in the Pechora Sea and around Svalbard where tides are strong. Tides increase the salt flux from sea ice by 50% in the Pechora and White Seas but have little impact elsewhere. This study suggests that the interaction between ocean tides and sea ice should not be neglected when modeling the Arctic.

Upper ocean stratification and sea ice growth rates during the summer-fall transition, as revealed by Elephant seal foraging in the Adélie Depression, East Antarctica - G. D. Williams; M. Hindell; M.-N. Houssais; T. Tamura; I. C. Field (Ocean Science, 2011)
Southern elephant seals (Mirounga leonina), fitted with Conductivity-Temperature-Depth sensors at Macquarie Island in January 2005 and 2010, collected unique oceanographic observations of the Adélie and George V Land continental shelf (140–148° E) during the summer-fall transition (late February through April). This is a key region of dense shelf water formation from enhanced sea ice growth/brine rejection in the local coastal polynyas. In 2005, two seals occupied the continental shelf break near the grounded icebergs at the northern end of the Mertz Glacier Tongue for several weeks from the end of February. One of the seals migrated west to the Dibble Ice Tongue, apparently utilising the Antarctic Slope Front current near the continental shelf break. In 2010, immediately after that year's calving of the Mertz Glacier Tongue, two seals migrated to the same region but penetrated much further southwest across the Adélie Depression and sampled the Commonwealth Bay polynya from March through April. Here we present observations of the regional oceanography during the summer-fall transition, in particular (i) the zonal distribution of modified Circumpolar Deep Water exchange across the shelf break, (ii) the upper ocean stratification across the Adélie Depression, including alongside iceberg C-28 that calved from the Mertz Glacier and (iii) the convective overturning of the deep remnant seasonal mixed layer in Commonwealth Bay from sea ice growth. Heat and freshwater budgets to 200–300 m are used to estimate the ocean heat content (400→50 MJ m−2), flux (50–200 W m−2 loss) and sea ice growth rates (maximum of 7.5–12.5 cm day−1). Mean seal-derived sea ice growth rates were within the range of satellite-derived estimates from 1992–2007 using ERA-Interim data. We speculate that the continuous foraging by the seals within Commonwealth Bay during the summer/fall transition was due to favorable feeding conditions resulting from the convective overturning of the deep seasonal mixed layer and chlorophyll maximum that is a reported feature of this location.

A global comparison of Argo and satellite altimetry observations - A. -L. Dhomps; S. Guinehut; P.-Y. Le Traon; G. Larnicol (Ocean Science, 2011)
Differences, similarities and complementarities between Sea Level Anomalies (SLA) deduced from altimeter measurements and dynamic height anomalies (DHA) calculated from Argo in situ temperature (T) and salinity (S) profiles are globally analyzed. SLA and DHA agree remarkably well and, compared to previous studies, Argo dataset allows an improvement in the coherence between SLA and DHA. Indeed, Argo data provides a much better spatial coverage of all oceans and particularly the Southern Ocean, the use of an Argo mean dynamic height, the use of measured salinity profiles (versus climatological salinity), and the use of a deeper reference level (1000 m versus 700 m). The large influence of Argo salinity observations on the consistency between altimetry and hydrographic observations is particularly demonstrated with an improvement of 35% (relative to the SLA minus DHA signal) by using measured salinity profiles instead of climatological data. The availability of observations along the Argo float trajectories also provides a means to describe the sea level variability of the global ocean both for the low frequency and the mesoscale part of the circulation. Results indicate that sea level variability is dominated by baroclinic signal at seasonal to inter-annual periods for all latitudes. In the tropics, sea level variability is baroclinic for meso-scale to interannual periods and at high latitudes, sea level variability is barotropic with also deep baroclinic signals (i.e. influence of deep temperature and salinity signals) for intra seasonal and mesoscale periods. These results emphasize the need to separate the different time and space scales in order to improve the merging of the two data sets. The qualitative study of seasonal to interannual SLA minus DHA signals finally reveals signals related to deep ocean circulation variations and basin-scale barotropic signals. Future work is, however, needed to understand the observed differences and relate them to different forcing mechanisms.

Effect of tidal stream power generation on the region-wide circulation in a shallow sea - G. I. Shapiro (Ocean Science, 2011)
This paper quantifies the backward effect on the ocean currents caused by a tidal stream farm located in the open shallow sea. Recent studies in channels with 1-D models have indicated that the power potential is not given purely by the flux of kinetic energy, as has been commonly assumed. In this study, a 3-D ocean circulation model is used to estimate (i) practically extractable energy resource at different levels of rated generation capacity of the farm, (ii) changes in the strength of currents due to energy extraction, and (iii) alterations in the pattern of residual currents and the pathways of passive tracers. As well as tidal streams, the model also takes into account the wind-driven and density-driven ocean currents. Numerical modelling has been carried out for a hypothetical tidal farm located in the Celtic Sea north of Cornwall, an area known for its high level of tidal energy. Modelling results clearly indicate that the extracted power does not grow linearly with the increase in the rated capacity of the farm. For the case study covered in this paper, a 100-fold increase in the rated generation capacity of the farm results in only 7-fold increase in extracted power. In the case of a high power farm, kinetic energy of currents is altered significantly as far as 10–20 km away from the farm. At high levels of extracted energy the currents tend to avoid flowing through the farm, an effect which is not captured with 1-D models. Residual currents are altered as far as a hundred kilometres away. The magnitude of changes in the dispersion of tracers is highly sensitive to the location. Some of the passive drifters analysed in this study experience significant variations in the end-to-start distance due to energy extraction ranging from 13% to 238% while others are practically unaffected. This study shows that both energy extraction estimates and effects on region wide circulation depend on a complex combination of factors, and the specific figures given in the paper should be generally considered as first estimates.

An ocean modelling and assimilation guide to using GOCE geoid products - K. Haines; J. A. Johannessen; P. Knudsen; D. Lea; M.-H. Rio; L. Bertino; F. Davidson; F. Hernandez (Ocean Science, 2011)
We review the procedures and challenges that must be considered when using geoid data derived from the Gravity and steady-state Ocean Circulation Explorer (GOCE) mission in order to constrain the circulation and water mass representation in an ocean general circulation model. It covers the combination of the geoid information with time-mean sea level information derived from satellite altimeter data, to construct a mean dynamic topography (MDT), and considers how this complements the time-varying sea level anomaly, also available from the satellite altimeter. We particularly consider the compatibility of these different fields in their spatial scale content, their temporal representation, and in their error covariances. These considerations are very important when the resulting data are to be used to estimate ocean circulation and its corresponding errors. We describe the further steps needed for assimilating the resulting dynamic topography information into an ocean circulation model using three different operational forecasting and data assimilation systems. We look at methods used for assimilating altimeter anomaly data in the absence of a suitable geoid, and then discuss different approaches which have been tried for assimilating the additional geoid information. We review the problems that have been encountered and the lessons learned in order the help future users. Finally we present some results from the use of GRACE geoid information in the operational oceanography community and discuss the future potential gains that may be obtained from a new GOCE geoid.

Long-term spatial variations in the Baltic Sea wave fields - T. Soomere; A. Räämet (Ocean Science, 2011)
This study focuses on spatial patterns in linear trends of numerically reconstructed basic wave properties (average and extreme wave heights, wave periods) in the Baltic Sea under the assumption of no ice cover. Numerical simulations of wave conditions for 1970–2007, using the WAM wave model and adjusted geostrophic winds, revealed extensive spatial variations in long-term changes in both average and extreme wave heights in the Baltic Sea but almost no changes in the basinwide wave activity and wave periods. There has been a statistically significant decrease in the annual mean significant wave height by more than 10% between the islands of Öland and Gotland and in the southward sea area, and a substantial increase to the south-west of Bornholm, near the coast of Latvia, between the Åland Archipelago and the Swedish mainland, and between the Bothnian Sea and the Bothnian Bay. Variations in extreme wave heights (defined as the threshold for 1% of the highest waves each year) show similar patterns of changes. In several areas the trends in average and extreme wave heights are different. Such a complicated pattern of changes indicates that (i) different regions of the Baltic Sea basin have experienced widespread but essentially different changes in wind properties and (ii) many seemingly controversial trends and variations established in wave properties at different sites in the recent past may reflect the natural spatial variability in the Baltic Sea wave fields.

Measurements of bubble size spectra within leads in the Arctic summer pack ice - S. J. Norris; I. M. Brooks; G. de Leeuw; A. Sirevaag; C. Leck; B. J. Brooks; C. E. Birch; M. Tjernström (Ocean Science, 2011)
The first measurements of bubble size spectra within the near-surface waters of open leads in the central Arctic pack ice were obtained during the Arctic Summer Cloud-Ocean Study (ASCOS) in August 2008 at 87–87.6° N, 1–11° W. A significant number of small bubbles (30–100 μm diameter) were present, with concentration decreasing rapidly with size from 100–560 μm; no bubbles larger than 560 μm were observed. The bubbles were present both during periods of low wind speed (U

Phytoplankton distribution and nitrogen dynamics in the southwest indian subtropical gyre and Southern Ocean waters - S. J. Thomalla; H. N. Waldron; M. I. Lucas; J. F. Read; I. J. Ansorge; E. Pakhomov (Ocean Science, 2011)
During the 1999 Marion Island Oceanographic Survey (MIOS 4) in late austral summer, a northbound and reciprocal southbound transect were taken along the Southwest Indian and Madagascar Ridge, between the Prince Edward Islands and 31° S. The sections crossed a number of major fronts and smaller mesoscale features and covered a wide productivity spectrum from subtropical to subantarctic waters. Associated with the physical survey were measurements of size fractionated chlorophyll, nutrients and nitrogen (NO3, NH4 and urea) uptake rates. Subtropical waters were characterised by low chlorophyll concentrations (max = 0.27.3 mg m−3 dominated by pico-phytoplankton cells (> 81%) and very low f-ratios (

Corrigendum to "Influence of cross-shelf water transport on nutrients and phytoplankton in the East China Sea: a model study" published in Ocean Sci., 7, 27–43, 2010 - L. Zhao; X. Guo (Ocean Science, 2011)

Preface "Deep Ocean Exchange with the Shelf (DOES)" - J. Johnson; P. Chapman (Ocean Science, 2011)

Tracers confirm downward mixing of Tyrrhenian Sea upper waters associated with the Eastern Mediterranean Transient - W. Roether; J. E. Lupton (Ocean Science, 2011)
Observations of tritium and 3He in the Tyrrhenian Sea, 1987–2009, confirm the enhanced vertical mixing of intermediate waters into the deep waters that has been noted and associated with the Eastern Mediterranean Transient in previous studies. Our evidence for the mixing rests on increasing tracer concentrations in the Tyrrhenian deep waters, accompanied by decreases in the upper waters, which are supplied from the Eastern Mediterranean. The downward transfer is particularly evident between 1987 and 1997. Later on, information partly rests on increasing tritium-3He ages; here we correct the observed 3He for contributions released from the ocean floor. The Tyrrhenian tracer distributions are fully compatible with data upstream of the Sicily Strait and in the Western Mediterranean. The tracer data show that mixing reached to the bottom and confirm a cyclonic nature of the deep water circulation in the Tyrrhenian. They furthermore indicate that horizontal homogenization of the deep waters occurs on a time scale of roughly 5 years. Various features point to a reduced impact of Western Mediterranean Deep Water (WMDW) in the Tyrrhenian during the enhanced-mixing period. This is an important finding because it implies less upward mixing of WMDW, which has been named a major process to enable the WMDW to leave the Mediterranean via the Gibraltar Strait. On the other hand, the TDW outflow for several years represented a major influx of enhanced salinity and density waters into the deep-water range of the Western Mediterranean.

Assessment of the three dimensional temperature and salinity observational networks in the Baltic Sea and North Sea - W. Fu; J. L. Høyer; J. She (Ocean Science, 2011)
The spatial averaged correlations are presented in 1.5° × 1.5° bins for the North and Baltic Sea region. The averaged correlations are computed based on the proxy ocean data generated by the operational forecast model of Danish Meteorology Institute (DMI). It is shown that the spatial distribution of the averaged correlations could reflect the overall influence of the local atmospheric forcing, complex topography, coastlines, boundary and bottom effect, etc. Comparisons with the satellite SST data demonstrate that the proxy ocean data reproduce realistic results at the surface. Based on the spatial bin-averaged correlations, a general correlation model is assumed to approximate the spatial and temporal correlation structure. Parameters of the correlation model are obtained on the standard Levitus levels. It is found that the correlation model is not the typical Guaussian-type function. For instance, the exponents of the correlation model vary in the longitudinal direction from 0.75 at the surface to 1.33 at the depth of 250 m for temperature. For salinity, the temporal correlation can be approximated with an exponential function. Two complementary quality-indicators, effective coverage rate and "explained" variance, are defined based on the correlation models obtained above. The two indicators are able to identify the "influence area" of the information content in a given observation network and the relative importance of observations at different locations. By these indicators, the 3-D temperature and salinity observational networks are assessed in the Baltic Sea and North Sea for the period 2004–2006. It is found that the surface level is more effectively covered than the deep waters with existing networks. In addition, the Belt Sea and the Baltic Proper also show good coverage for both temperature and salinity. However, more observations are required in the Norwegian Trench and Kattegat. In the vertical, the two indicators show smaller values from 50 m to 125 m in this region, indicating the need for more observations.

Salinity-induced mixed and barrier layers in the southwestern tropical Atlantic Ocean off the northeast of Brazil - M. Araujo; C. Limongi; J. Servain; M. Silva; F. S. Leite; D. Veleda; C. A. D. Lentini (Ocean Science, 2011)
High-resolution hydrographic observations of temperature and salinity are used to analyze the formation and distribution of isothermal depth (ZT), mixed depth (ZM) and barrier layer thickness (BLT) in a section of the southwestern Atlantic (0°30´ N–14°00´ S; 31°24´–41°48´ W), adjacent to the northeastern Brazilian coast. Analyzed data consists of 279 CTD casts acquired during two cruises under the Brazilian REVIZEE Program. One occurred in late austral winter (August–October 1995) and another in austral summer (January–April 1997). Oceanic observations are compared to numerical modeling results obtained from the French Mercator-Coriolis Program. Results indicate that the intrusion of subtropical Salinity Maximum Waters (SMW) is the major process contributing to the seasonal barrier layer formation. These waters are brought by the South Equatorial Current (SEC), from the subtropical region, into the western tropical Atlantic boundary. During late austral winter southeastern trade winds are more intense and ITCZ precipitations induce lower surface salinity values near the equator. During this period a 5–90 m thick BLT (median = 15 m) is observed and BLT > 30 m is restricted to latitudes higher than 8° S, where the intrusion of salty waters between 8°–12.3° S creates shallow mixed layers over deep (ZT ≥ 90 m) isothermal layers. During austral summer, shallow isothermal and mixed layers prevail, when northeasterly winds are predominant and evaporation overcomes precipitation, causing saltier waters at the surface/subsurface layers. During that period observed BLT varies from 5 to 70 m and presents thicker median value of 35 m, when comparing to the winter. Furthermore, BLT ≥ 30 m is observed not only in the southernmost part of the study area, as verified during late winter, but in the latitude range 2°–14° S, where near-surface salty waters are transported westward by the SEC flow. These results indicate that the inclusion of salinity dynamics and its variability are necessary for studying mixed and barrier layer behaviors in the tropical Atlantic, where ocean-atmosphere coupling is known to be stronger.

Metrological traceability of oceanographic salinity measurement results - S. Seitz; R. Feistel; D. G. Wright; S. Weinreben; P. Spitzer; P. De Bièvre (Ocean Science, 2011)
Consistency of observed oceanographic salinity data is discussed with respect to contemporary metrological concepts. The claimed small uncertainty of salinity measurement results traceable to the conductivity ratio of a certified IAPSO Standard Seawater reference is not metrologically justified if results are compared on climatic time scales. This applies in particular to Practical Salinity SP, Reference Salinity SR, and the latest estimates of Absolute Salinity using the TEOS-10 formalism. On climate time scales an additional contribution to the uncertainty that is related to unknown property changes of the reference material must be accounted for. Moreover, when any of these measured or calculated quantity values is used to estimate Absolute Salinity of a seawater sample under investigation, another uncertainty contribution is required to quantify the accuracy of the equations relating the actually measured quantity to the Absolute Salinity. Without accounting for these additional uncertainties, such results cannot be used to estimate Absolute Salinity with respect to the International System of Units (SI), i.e. to the unit chosen for the mass fraction of dissolved material in the sample, which is "g kg−1". From a metrological point of view, such deficiencies in the calculations involving other quantities will produce SI-incompatible results. We outline how these problems can be overcome by linking salinity to primary SI measurement standards.

Influence of cross-shelf water transport on nutrients and phytoplankton in the East China Sea: a model study - L. Zhao; X. Guo (Ocean Science, 2011)
A three dimensional coupled biophysical model was used to examine the supply of oceanic nutrients to the shelf of the East China Sea (ECS) and its role in primary production over the shelf. The model consisted of two parts: the hydrodynamic module was based on a nested model with a horizontal resolution of 1/18 degree, whereas the biological module was a lower trophic level ecosystem model including two types of phytoplankton, three elements of nutrients, and biogenic organic material. The model results suggested that seasonal variations occurred in the distribution of nutrients and chlorophyll a over the shelf of the ECS. After comparison with available observed nutrients and chlorophyll a data, the model results were used to calculate volume and nutrients fluxes across the shelf break. The annual mean total fluxes were 1.53 Sv for volume, 9.4 kmol s−1 for DIN, 0.7 kmol s−1 for DIP, and 18.2 kmol s−1 for silicate. Two areas, northeast of Taiwan and southwest of Kyushu, were found to be major source regions of oceanic nutrients to the shelf. Although the onshore fluxes of nutrients and volume both had apparent seasonal variations, the seasonal variation of the onshore nutrient flux did not exactly follow that of the onshore volume flux. Additional calculations in which the concentration of nutrients in Kuroshio water was artificially increased suggested that the oceanic nutrients were distributed in the bottom layer from the shelf break to the region offshore of the Changjiang estuary from spring to summer and appeared in the surface layer from autumn to winter. The calculations also implied that the supply of oceanic nutrients to the shelf can change the consumption of pre-existing nutrients from rivers. The response of primary production over the shelf to the oceanic nutrients was confirmed not only in the surface layer (mainly at the outer shelf and shelf break in winter and in the region offshore of the Changjiang estuary in summer) but also in the subsurface layer over the shelf from spring to autumn.

Absolute Salinity, ''Density Salinity'' and the Reference-Composition Salinity Scale: present and future use in the seawater standard TEOS-10 - D. G. Wright; R. Pawlowicz; T. J. McDougall; R. Feistel; G. M. Marion (Ocean Science, 2011)
Salinity plays a key role in the determination of the thermodynamic properties of seawater and the new TEOS-101 standard provides a consistent and effective approach to dealing with relationships between salinity and these thermodynamic properties. However, there are a number of practical issues that arise in the application of TEOS-10, both in terms of accuracy and scope, including its use in the reduction of field data and in numerical models. First, in the TEOS-10 formulation for IAPSO Standard Seawater, the Gibbs function takes the Reference Salinity as its salinity argument, denoted SR, which provides a measure of the mass fraction of dissolved material in solution based on the Reference Composition approximation for Standard Seawater. We discuss uncertainties in both the Reference Composition and the Reference-Composition Salinity Scale on which Reference Salinity is reported. The Reference Composition provides a much-needed fixed benchmark but modified reference states will inevitably be required to improve the representation of Standard Seawater for some studies. However, the Reference-Composition Salinity Scale should remain unaltered to provide a stable representation of salinity for use with the TEOS-10 Gibbs function and in climate change detection studies. Second, when composition anomalies are present in seawater, no single salinity variable can fully represent the influence of dissolved material on the thermodynamic properties of seawater. We consider three distinct representations of salinity that have been used in previous studies and discuss the connections and distinctions between them. One of these variables provides the most accurate representation of density possible as well as improvements over Reference Salinity for the determination of other thermodynamic properties. It is referred to as "Density Salinity" and is represented by the symbol SAdens; it stands out as the most appropriate representation of salinity for use in dynamical physical oceanography. The other two salinity variables provide alternative measures of the mass fraction of dissolved material in seawater. "Solution Salinity", denoted SAsoln, is the most obvious extension of Reference Salinity to allow for composition anomalies; it provides a direct estimate of the mass fraction of dissolved material in solution. "Added-Mass Salinity", denoted SAadd, is motivated by a method used to report laboratory experiments; it represents the component of dissolved material added to Standard Seawater in terms of the mass of material before it enters solution. We also discuss a constructed conservative variable referred to as "Preformed Salinity", denoted S∗, which will be useful in process-oriented numerical modelling studies. Finally, a conceptual framework for the incorporation of composition anomalies in numerical models is presented that builds from studies in which composition anomalies are simply ignored up to studies in which the influences of composition anomalies are accounted for using the results of biogeochemical models. 1TEOS-10: international Thermodynamic Equation of Seawater 2010, http://www.teos-10.org/.

Stability and forcing of the Iceland-Faroe inflow of water, heat, and salt to the Arctic - B. Hansen; H. Hátún; R. Kristiansen; S. M. Olsen; S. Østerhus (Ocean Science, 2010)
The flow of Atlantic water across the Greenland-Scotland Ridge (Atlantic inflow) is critical for conditions in the Nordic Seas and Arctic Ocean by importing heat and salt. Here, we present a decade-long series of measurements from the Iceland-Faroe inflow branch (IF-inflow), which carries almost half the total Atlantic inflow. The observations show no significant trend in volume transport of Atlantic water, but temperature and salinity increased during the observational period. On shorter time scales, the observations show considerable variations but no statistically significant seasonal variation is observed and even weekly averaged transport values were consistently uni-directional from the Atlantic into the Nordic Seas. Combining transport time-series with sea level height from satellite altimetry and wind stress reveals that the force driving the IF-inflow across the topographic barrier of the Ridge is mainly generated by a pressure gradient that is due to a continuously maintained low sea level in the Southern Nordic Seas. This implies that the relative stability of the IF-inflow derives from the processes that lower the sea level by generating outflow from the Nordic Seas, especially the thermohaline processes that generate overflow. The IF-inflow is an important component of the system coupling the Arctic region to the North Atlantic through the thermohaline circulation, which has been predicted to weaken in the 21st century. Our observations show no indication of weakening.

Retroflection from a double-slanted coastline: a model for the Agulhas leakage variability - V. Zharkov; D. Nof; W. Weijer (Ocean Science, 2010)
The Agulhas leakage to the South Atlantic exhibits a strong anti-correlation with the mass flux of the Agulhas Current. When the Agulhas retroflection is in its normal position near Cape Agulhas, leakage is relatively high and the nearby South African coastal slant (angle of derivation from zonal) is very small and relatively invariant alongshore. During periods of strong incoming flux (low leakage), the retroflection shifts upstream to Port Elizabeth or East London, where the coastline shape has a "kink", i.e., the slant changes abruptly from small on the west side, to large (about 55°) on the east side. Here, we show that the variability of rings shedding and anti-correlation between Agulhas mass flux and leakage to the South Atlantic may be attributed to this kink. To do so, we develop a nonlinear analytical model for retroflection near a coastline that consists of two sections, a zonal western section and a strongly slanted eastern section. The principal difference between this and the model of a straight slanted coast (discussed in our earlier papers) is that, here, free purely westward propagation of eddies along the zonal coastline section is allowed. This introduces an interesting situation in which strong slant of the coast east of the kink prohibits the formation and shedding of rings, while the almost zonal coastal orientation west of the kink encourages shedding. Therefore, the kink "locks" the position of the retroflection, forcing it to occur just downstream of the kink. Rings are necessarily shed from the retroflection area in our kinked model, regardless of the degree of eastern coast slant. In contrast, a no-kink model with a coastline of intermediate slant indicates that shedding is almost completely arrested by that slant. We suggest that the observed difference in ring-shedding intensity during times of normal retroflection position and times when the retroflection is shifted eastward is due to the change in the retroflection location with respect to the kink. When the incoming flux detaches from the coast north of the kink, ring transport is small; when the flux detaches south of the kink, transport is large. Simple process-oriented numerical simulations are in fair agreement with our analytical results.

The influence of the Brazil and Malvinas Currents on the Southwestern Atlantic Shelf circulation - R. P. Matano; E. D. Palma; A. R. Piola (Ocean Science, 2010)
The oceanic circulation over the southwestern Atlantic shelf is influenced by large tidal amplitudes, substantial freshwater discharges, high wind speeds and – most importantly – by its proximity to two of the largest western boundary currents of the world ocean: the Brazil and Malvinas currents. This review article aims to discriminate the dynamical processes controlling the interaction between this extensive shelf region and the deep-ocean. The discussion is focused on two broad regions: the South Brazil Bight to the north, and Patagonia to the south. The exchanges between the Brazil Current and the South Brazil Bight are characterized by the intermittent development of eddies and meanders of the Brazil Current at the shelfbreak. However, it is argued that this is not the only – nor the most important – influence of the Brazil Current on the shelf. Numerical simulations show that the thermohaline structure of the South Brazil Bight can be entirely ascribed to steady state, bottom boundary layer interactions between the shelf and the Brazil Current. The Malvinas Current does not show the development of eddies and meanders, but its influence on the Patagonian shelf is not less important. Models and observations indicate that the Malvinas Current not only controls the shelfbreak dynamics and cross-shelf exchanges but also influences the circulation in the shelf's interior.

Thermophysical property anomalies of Baltic seawater - R. Feistel; G. M. Marion; R. Pawlowicz; D. G. Wright (Ocean Science, 2010)
While the thermodynamic properties of Standard Seawater are very well known, the quantitative effect of sea salt composition anomalies on various properties is difficult to estimate since comprehensive lab experiments with the various natural waters are scarce. Coastal and estuarine waters exhibit significant anomalies which also influence to an unknown amount the routine salinity calculation from conductivity measurements. Recent numerical models of multi-component aqueous electrolytes permit the simulation of physical chemical properties of seawater with variable solute composition. In this paper, the FREZCHEM model is used to derive a Gibbs function for Baltic seawater, and the LSEA_DELS model to provide estimates for the conductivity anomaly relative to Standard Seawater. From additional information such as direct density measurements or empirical salinity anomaly parameterisation, the quantitative deviations of properties between Baltic and Standard Seawater are calculated as functions of salinity and temperature. While several quantities show anomalies that are comparable with their measurement uncertainties and do not demand special improvement, others exhibit more significant deviations from Standard Seawater properties. In particular density and sound speed turn out to be significantly sensitive to the presence of anomalous solute. Suitable general correction methods are suggested to be applied to Baltic Sea samples with known Practical Salinity and, optionally, directly determined density.

Spatio-Temporal Complexity analysis of the Sea Surface Temperature in the Philippines - Z. T. Botin; L. T. David; R. C. H. del Rosario; L. Parrott (Ocean Science, 2010)
A spatio-temporal complexity (STC) measure which has been previously used to analyze data from terrestrial ecosystems is employed to analyse 21 years of remotely sensed sea-surface temperature (SST) data from the Philippines. STC on the Philippine wide SST showed the monsoonal variability of the Philippine waters. STC is correlated with the SST mean (R2 ≈ 0.7), and inversely correlated with the SST standard deviation (R2 ≈ 0.9). Both STC and SST are highest during the middle of the year, which coincides with the Southwest Monsoon, but with the STC values being higher towards the end of the monsoon until the start of the inter-monsoon. In order to determine if STC has the potential to define limits of bio-regions, the spatial domain was subsequently divided into six thermal regions computed via clustering of temperature means. STC and EOF of the STC values were computed for each thermal region. Our STC analysis of the SST data, and comparisons with SST values suggest that the STC measure may be useful for characterising environmental heterogeneity over space and time for many long-term remotely sensed data.

A numerical scheme to calculate temperature and salinity dependent air-water transfer velocities for any gas - M. T. Johnson (Ocean Science, 2010)
The ocean-atmosphere flux of a gas can be calculated from its measured or estimated concentration gradient across the air-sea interface and the transfer velocity (a term representing the conductivity of the layers either side of the interface with respect to the gas of interest). Traditionally the transfer velocity has been estimated from empirical relationships with wind speed, and then scaled by the Schmidt number of the gas being transferred. Complex, physically based models of transfer velocity (based on more physical forcings than wind speed alone), such as the NOAA COARE algorithm, have more recently been applied to well-studied gases such as carbon dioxide and DMS (although many studies still use the simpler approach for these gases), but there is a lack of validation of such schemes for other, more poorly studied gases. The aim of this paper is to provide a flexible numerical scheme which will allow the estimation of transfer velocity for any gas as a function of wind speed, temperature and salinity, given data on the solubility and liquid molar volume of the particular gas. New and existing parameterizations (including a novel empirical parameterization of the salinity-dependence of Henry's law solubility) are brought together into a scheme implemented as a modular, extensible program in the R computing environment which is available in the supplementary online material accompanying this paper; along with input files containing solubility and structural data for ~90 gases of general interest, enabling the calculation of their total transfer velocities and component parameters. Comparison of the scheme presented here with alternative schemes and methods for calculating air-sea flux parameters shows good agreement in general. It is intended that the various components of this numerical scheme should be applied only in the absence of experimental data providing robust values for parameters for a particular gas of interest.

A quantitative description of the Norwegian Atlantic Current by combining altimetry and hydrography - K. A. Mork; Ø. Skagseth (Ocean Science, 2010)
In this study, we quantified the mean flow and the variability in the Norwegian Atlantic Current (NwAC), including the individual branches, with associated error estimates. We accomplished this by combining repeated hydrographic data in the Svinøy section (at 62° N off the coast of Norway toward the northwest) with absolute dynamic sea surface topography data for the period from 1992–2009. The analysis shows a two-branched structure of the NwAC in this section, with calculated absolute velocities that are largely in accordance with the independent current measurements. The mean volume flux of the NwAC was estimated to be 5.1 ± 0.3 Sv (Sv = 106 m3 s−1). In terms of branches, the estimated 3.4 ± 0.3 Sv in the eastern branch is slightly below previous estimates. The estimated volume flux in the western branch is 1.7 ± 0.2 Sv, but a lack of robust estimates makes a comparison difficult. There is a significant seasonal cycle in the volume transport (the winter maximum is almost twice as large as the summer minimum) with a major contribution from the eastern branch. On the inter-annual scale, the volume flux and temperature are significantly and negatively correlated in both the western branch and in the total. Examining the results showed that increased volume flux is followed by a decrease in the stability of the upper water column, averaged over the Norwegian Sea, and a delay in the phytoplankton spring bloom at the Ocean Weather Station M (66° N, 2° E), by a lag of 1 year. The cause of this decrease in stability, whether a direct effect of the increased volume flux or a consequence of indirect effects, is however beyond the scope of this study.

Numerical simulations of spreading of the Persian Gulf outflow into the Oman Sea - M. Ezam; A. A. Bidokhti; A. H. Javid (Ocean Science, 2010)
A three dimensional numerical model namely POM (Princeton Ocean Model) and observational data are used to study the Persian Gulf outflow structure and its spreading pathways during 1992. In the model, the monthly wind speed data were taken from ICOADS (International Comprehensive Ocean-Atmosphere Data Set) and the monthly SST (sea surface temperatures) were taken from AVHRR (Advanced Very High Resolution Radiometer) with the addition of monthly net shortwave radiations from NCEP (National Center for Environmental Prediction). The mean monthly precipitation rates from NCEP data and the calculated evaporation rates are used to impose the surface salinity fluxes. At the open boundaries the temperature and salinity were prescribed from the mean monthly climatological values from WOA05 (World Ocean Atlas 2005). Also the four major components of the tide were prescribed at the open boundaries. The results show that the outflow mainly originates from two branches at different depths in the Persian Gulf. The permanent branch exists during the whole year deeper than 40 m along the Gulf axis and originates from the inner parts of the Persian Gulf. The other seasonal branch forms in the vicinity of the shallow southern coasts due to high evaporation rates during winter. Near the Strait of Hormuz the two branches join and form the main outflow source water. The results of simulations reveal that during the winter the outflow boundary current mainly detaches from the coast well before Ras Al Hamra Cape, however during summer the outflow seems to follow the coast even after this Cape. This is due to a higher density of the colder outflow that leads to more sinking near the coast in winter. Thus, the outflow moves to a deeper depth of about 500 m (for which some explanations are given) while the main part detaches and spreads at a depth of about 300 m. However in summer it all moves at a depth of about 200–250 m. During winter, the deeper, stronger and wider outflow is more affected by the steep topography, leading to separation from the coast. While during summer, the weaker and shallower outflow is less influenced by bottom topography and so continues along the boundary.

Snapshot observation of the physical structure and stratification in deep-water of the South Caspian Sea (western part) - P. Ghaffari; H. A. Lahijani; J. Azizpour (Ocean Science, 2010)
The physical parameters structures and different water masses using CTD measurements in southwestern part of the Caspian Sea (CS) adjacent to Anzali Port (AP) are investigated. CTD profiles were conducted along a transect perpendicular to the coastline on 13 stations from the coast down to 720 m on 22 January 2008. According to the results the continental shelf waters are located in the surface mixed layer. The surface mixed layer extends itself down to almost 100 m in outer areas of the continental shelf with a weak seasonal thermocline layer between 80 to 140 m Freshwaters inflow of local rivers is clearly seen outside the continental shelf at the surface layers. Investigating the dissolved oxygen reveals that winter convection is traceable down to 500 m in the lateral waters over the shelf break. Among the deeper stations that are located in continental rise and abyssal plain, 300 m seems to be a threshold for penetration of seasonal changes; therefore deeper waters tend to be impermeable against seasonal variances. Despite the small variations, stability is positive in the study area and temperature plays an important role in static stability and in triggering the lateral mixing. In view of both temperature-salinity and temperature-oxygen distributions in the southwestern part of the CS, two different water masses are separable in cold phase. Snapshot observations of physical properties in the early winter 2008, to some extent revealed that a mixing was triggered at least in the lateral waters of the study area.

Estimates of radiance reflected towards the zenith at the surface of the sea - E. Aas (Ocean Science, 2010)
Remote sensing of water colour by ship-mounted sensors represents an important tool for the validation of satellite products and the monitoring of water quality. The recorded radiance from the sea has to be corrected for the surface-reflected radiance from sun and sky in order to obtain the water-leaving radiance. Here the simple case of radiance reflected towards the zenith is studied. A set of observed sky radiance and solar irradiance data from Oslo has been used together with a Gaussian slope distribution for the sea surface in order to estimate the reflected radiance. The spectral range studied is 405–650 nm, the solar zenith angles are in the range 37°–76°, and the wind speeds are up to 10 m s−1. The analysis of the results show that the reflected radiance has to be separated into three contributions: sky radiance and sun rays reflected at the foam-free surface and irradiance reflected by whitecaps and foam. It is then demonstrated that by using four input values, namely the downward irradiance, the sky radiance from the zenith, the solar zenith angle and the wind speed, it is possible to obtain by simple expressions estimates of the reflected radiance that only differ from the former calculated values by relative errors of less than 5%. The analysis also indicates that for the spectral range studied neither the water-leaving radiance nor the surface-reflected radiance can be disregarded relative to the other one in the Case 2 waters of the Oslofjord-Skagerrak area. The results form a first step towards the study of reflected radiance in viewing angles differing from the nadir direction.

On the seasonal cycles and variability of Florida Straits, Ekman and Sverdrup transports at 26° N in the Atlantic Ocean - C. P. Atkinson; H. L. Bryden; J. J-M. Hirschi; T. Kanzow (Ocean Science, 2010)
Since April 2004 the RAPID array has made continuous measurements of the Atlantic Meridional Overturning Circulation (AMOC) at 26° N. Two key components of this system are Ekman transport zonally integrated across 26° N and western boundary current transport in the Florida Straits. Whilst measurements of the AMOC as a whole are somewhat in their infancy, this study investigates what useful information can be extracted on the variability of the Ekman and Florida Straits transports using the decadal timeseries already available. Analysis is also presented for Sverdrup transports zonally integrated across 26° N. The seasonal cycles of Florida Straits, Ekman and Sverdrup transports are quantified at 26° N using harmonic analysis of annual and semi-annual constituents. Whilst Sverdrup transport shows clear semi-annual periodicity, calculations of seasonal Florida Straits and Ekman transports show substantial interannual variability due to contamination by variability at non-seasonal frequencies; the mean seasonal cycle for these transports only emerges from decadal length observations. The Florida Straits and Ekman mean seasonal cycles project on the AMOC with a combined peak-to-peak seasonal range of 3.5 Sv. The combined seasonal range for heat transport is 0.40 PW. The Florida Straits seasonal cycle possesses a smooth annual periodicity in contrast with previous studies suggesting a more asymmetric structure. No clear evidence is found to support significant changes in the Florida Straits seasonal cycle at sub-decadal periods. Whilst evidence of wind driven Florida Straits transport variability is seen at sub-seasonal and annual periods, a model run from the 1/4° eddy-permitting ocean model NEMO is used to identify an important contribution from internal oceanic variability at sub-annual and interannual periods. The Ekman transport seasonal cycle possesses less symmetric structure, due in part to different seasonal transport regimes east and west of 50 to 60° W. Around 60% of non-seasonal Ekman transport variability occurs in phase section-wide at 26° N and is related to the NAO, whilst Sverdrup transport variability is more difficult to decompose.

Upper layer current variability in the Central Ligurian Sea - P. Picco; A. Cappelletti; S. Sparnocchia; M. E. Schiano; S. Pensieri; R. Bozzano (Ocean Science, 2010)
Long-time series of surface currents and meteorological parameters were analysed to estimate the variability of the upper layer circulation and the response to the local winds. Current meter data were collected by an upward-looking RDI Sentinel 300 kHz ADCP deployed in the Central Ligurian Sea (43°47.77' N; 9°02.85' E) near the meteo-oceanographic buoy ODAS Italia 1 for more than eight months, from 13th of September 2003 to 24th of May 2004. The ADCP sampled the upper 50 m of water column at 8 m vertical resolution and 1 h time interval; surface marine and atmospheric hourly averaged data were provided by the buoy. Currents in the sampled layer were mainly barotropic, directed North-West in accordance with the general circulation of the area, and had a mean velocity of about 18 cm/s and hourly mean peaks up to 80 cm/s. Most of the observed variability in the upper thermocline was determined by inertial currents and mesoscale activity due to the presence of the Ligurian Front. Local wind had a minor role in the near-surface circulation but induced internal waves propagating downward in the water column.

Trends in coastal upwelling intensity during the late 20th century - N. Narayan; A. Paul; S. Mulitza; M. Schulz (Ocean Science, 2010)
This study presents linear trends of coastal upwelling intensity in the later part of the 20th century (1960–2001) employing various indices of upwelling, derived from meridional wind stress and sea surface temperature. The analysis was conducted in the four major coastal upwelling regions in the world, which are off North-West Africa, Lüderitz, California and Peru. The trends in meridional wind stress showed a steady increase of intensity from 1960–2001, which was also reflected in the SST index calculated for the same time period. The steady cooling observed in the instrumental records of SST off California substantiated this observation further. It was also noted that the trends in meridional wind stress obtained from different datasets differ substantially from each other. Correlation analysis showed that basin-scale oscillations like the Atlantic Multidecadal Oscillation (AMO) and the Pacific Decadal Oscillation (PDO) could not be directly linked to the observed increase of upwelling intensity off NW Africa and California respectively. The relationship of the North Atlantic Oscillation (NAO) with coastal upwelling off NW Africa turned out to be ambiguous due to a negative correlation between the NAO index and the meridional wind stress and a lack of correlation with the SST index. Our results give additional support to the hypothesis that the coastal upwelling intensity increases globally because of raising greenhouse gas concentrations in the atmosphere and an associated increase of the land-sea pressure gradient and meridional wind stress.

Suspended particles in the Canada Basin from optical and bottle data, 2003–2008 - J. M. Jackson; S. E. Allen; E. C. Carmack; F. A. McLaughlin (Ocean Science, 2010)
It is expected that coastal erosion, upwelling, and increased river runoff from Arctic warming will increase the concentration of suspended particles in the Arctic Ocean. Here we analyze in situ transmissometer and fluorometer data from the summers of 2003 through 2008 and bottle-derived particulate organic carbon (POC) and total suspended solids (TSS) measurements sampled in the summers of 2006 and 2007 from the Canada Basin and surrounding shelves. We divided our study area into five regions to account for the significant spatial variability and found that the highest attenuation, POC and TSS values were observed along the Beaufort shelf and the lowest values were located along the eastern shelf of the Canada Basin. We then explored the correlation of POC and TSS with beam attenuation coefficients to assess the viability of estimating POC concentrations from archived transmissometer data. POC (but not TSS) and attenuation were well-correlated over the Northwind Ridge, in the Canada Basin interior, and along the eastern shelf of the Canada Basin. Neither TSS nor POC were well-correlated with attenuation along the entire Beaufort shelf. An interannual comparison of the attenuation and fluorescence data was done. We found no evidence of increasing attenuation from the summers of 2003 through 2008 and, although not statistically significant, it even appeared that attenuation decreased over time in the upper 25 m of the Northwind Ridge and in the 25–100 m layer (that includes the chlorophyll maximum) of the eastern Beaufort shelf and within the Canada Basin. In the Canada Basin interior, the subsurface chlorophyll maximum deepened at a rate of 3.2 m per year from an average of 45 m in 2003 to 61 m in 2008, an example of how changes to the Arctic climate are impacting its ecology.

Modeling the Nd isotopic composition in the North Atlantic basin using an eddy-permitting model - T. Arsouze; A. M. Treguier; S. Peronne; J.-C. Dutay; F. Lacan; C. Jeandel (Ocean Science, 2010)
Boundary Exchange (BE – exchange of elements between continental margins and the open ocean) has been emphasized as a key process in the oceanic cycle of neodymium (Nd) (Lacan and Jeandel, 2005a). Here, we use a regional eddy-permitting resolution Ocean General Circulation Model (1/4°) of the North Atlantic basin to simulate the distribution of the Nd isotopic composition, considering BE as the only source. Results show good agreement with the data, confirming previous results obtained using the same parameterization of the source in a coarse resolution global model (Arsouze et al., 2007), and therefore the major control played by the BE processes in the Nd cycle on the regional scale. We quantified the exchange rate of the BE, and found that the time needed for the continental margins to significantly imprint the chemical composition of the surrounding seawater (further referred as characteristic exchange time) is of the order of 0.2 years. However, the timescale of the BE may be subject to large variations as a very short exchange time (a few days) is needed to reproduce the highly negative values of surface waters in the Labrador Sea, whereas a longer one (up to 0.5 years) is required to simulate the radiogenic influence of basaltic margins and distinguish the negative isotopic signatures of North Atlantic Deep Water from the more radiogenic southern origin water masses. This likely represents geographical variations in erosion fluxes and the subsequent particle load onto the continental margins. Although the parameterization of the BE is the same in both configurations of the model, the characteristic exchange time in the eddy-permitting configuration is significantly lower than the previous evaluations using a low resolution configuration (6 months to 10 years), but however in agreement with the available seawater Nd isotope data. This results highlights the importance of the model dynamics in simulating the BE process.

Seasonal cycles of surface layer salinity in the Pacific Ocean - F. M. Bingham; G. R. Foltz; M. J. McPhaden (Ocean Science, 2010)
The seasonal variability of surface layer salinity (SLS) is examined in the Pacific Ocean between 40° S and 60° N using a variety of data sources. Significant seasonal cycles were found in 5 regions: 1) The western North Pacific, 2) The northeastern North Pacific and Alaska gyre, 3) the intertropical convergence zone (ITCZ), 4) an area of the central North Pacific north of the Hawaiian Islands, 5) the central South Pacific along 10–20° S. Amplitudes range from 0.1 to > 0.5. The largest amplitudes are in the tropical band and the western North Pacific. Maximum salinity is obtained in late (northern) winter in the western North Pacific, late winter and early spring in the northeastern North Pacific, early summer in the ITCZ area, late summer and early fall in the central North Pacific area and (austral) winter in the central South Pacific. Large areas of the Pacific have no significant seasonal variation in SLS. Seasonal variability of evaporation rate, precipitation rate and the difference between them (E-P) were calculated from the OAFlux and Global Precipitation Climatology Project datasets. Typical amplitudes of E-P are 0.1–1 × 10−4 kg m−2 s−1. The seasonal variability of E-P is largely dominated by variability in evaporation in the western North Pacific and precipitation elsewhere. The largest amplitudes are in areas along the edge of the western North Pacific and in the far eastern tropical Pacific around 10° N. Phases in these areas indicate maximum E-P in mid- to late winter in these areas of large amplitude. The closest correspondence between E-P and SLS is in the ITCZ. E-P was combined with seasonal variation of the mixed-layer depth to calculate the freshwater flux forcing term of the SLS balance equation. The term was found to be similar in magnitude and distribution to E-P. Some other terms of the SLS balance were calculated. Horizontal advection was found to have seasonal cycles in a region near the equator. Entrainment was found to be mostly not significant except for a small region along 2.5–7.5° N in the eastern Pacific. Averaged spatially over large areas in the western North Pacific, ITCZ, South Pacific and northern North Pacific, the seasonal cycle is mostly a balance between changes in SLS and E-P, with entrainment and advection playing relatively minor roles. This work highlights the potentially significant role of surface salinity in the hydrologic cycle and in subtropical mode water formation. It can also help to interpret measurements that will soon be available from the Aquarius and SMOS (Soil Moisture and Ocean Salinity) satellite missions.

Impact of hydrographic data assimilation on the modelled Atlantic meridional overturning circulation - G. C. Smith; K. Haines; T. Kanzow; S. Cunningham (Ocean Science, 2010)
Here we make an initial step toward the development of an ocean assimilation system that can constrain the modelled Atlantic Meridional Overturning Circulation (AMOC) to support climate predictions. A detailed comparison is presented of 1° and 1/4° resolution global model simulations with and without sequential data assimilation, to the observations and transport estimates from the RAPID mooring array across 26.5° N in the Atlantic. Comparisons of modelled water properties with the observations from the merged RAPID boundary arrays demonstrate the ability of in situ data assimilation to accurately constrain the east-west density gradient between these mooring arrays. However, the presence of an unconstrained "western boundary wedge" between Abaco Island and the RAPID mooring site WB2 (16 km offshore) leads to the intensification of an erroneous southwards flow in this region when in situ data are assimilated. The result is an overly intense southward upper mid-ocean transport (0–1100 m) as compared to the estimates derived from the RAPID array. Correction of upper layer zonal density gradients is found to compensate mostly for a weak subtropical gyre circulation in the free model run (i.e. with no assimilation). Despite the important changes to the density structure and transports in the upper layer imposed by the assimilation, very little change is found in the amplitude and sub-seasonal variability of the AMOC. This shows that assimilation of upper layer density information projects mainly on the gyre circulation with little effect on the AMOC at 26° N due to the absence of corrections to density gradients below 2000 m (the maximum depth of Argo). The sensitivity to initial conditions was explored through two additional experiments using a climatological initial condition. These experiments showed that the weak bias in gyre intensity in the control simulation (without data assimilation) develops over a period of about 6 months, but does so independently from the overturning, with no change to the AMOC. However, differences in the properties and volume transport of North Atlantic Deep Water (NADW) persisted throughout the 3 year simulations resulting in a difference of 3 Sv in AMOC intensity. The persistence of these dense water anomalies and their influence on the AMOC is promising for the development of decadal forecasting capabilities. The results suggest that the deeper waters must be accurately reproduced in order to constrain the AMOC.

Short-term impacts of enhanced Greenland freshwater fluxes in an eddy-permitting ocean model - R. Marsh; D. Desbruyères; J. L. Bamber; B. A. de Cuevas; A. C. Coward; Y. Aksenov (Ocean Science, 2010)
In a sensitivity experiment, an eddy-permitting ocean general circulation model is forced with realistic freshwater fluxes from the Greenland Ice Sheet, averaged for the period 1991–2000. The fluxes are obtained with a mass balance model for the ice sheet, forced with the ERA-40 reanalysis dataset. The freshwater flux is distributed around Greenland as an additional term in prescribed runoff, representing seasonal melting of the ice sheet and a fixed year-round iceberg calving flux, for 8.5 model years. By adding Greenland freshwater fluxes with realistic geographical distribution and seasonality, the experiment is designed to investigate the oceanic response to a sudden and spatially/temporally uniform amplification of ice sheet melting and discharge, rather than localized or gradual changes in freshwater flux. The impacts on regional hydrography and circulation are investigated by comparing the sensitivity experiment to a control experiment, without additional fluxes. By the end of the sensitivity experiment, the majority of additional fresh water has accumulated in Baffin Bay, and only a small fraction has reached the interior of the Labrador Sea, where winter mixed layer depth is sensitive to small changes in salinity. As a consequence, the impact on large-scale circulation is very slight. An indirect impact of strong freshening off the west coast of Greenland is a small anti-cyclonic component to the circulation around Greenland, which opposes the wind-driven cyclonic circulation and reduces net southward flow through the Canadian Archipelago by ~10%. Implications for the post-2000 acceleration of Greenland mass loss are discussed.

Acoustic Doppler Current Profiler observations in the southern Caspian Sea: shelf currents and flow field off Feridoonkenar Bay, Iran - P. Ghaffari; V. Chegini (Ocean Science, 2010)
The results of offshore bottom-mounted ADCP measurements and wind records carried out from August to September 2003 in the coastal waters off Feridoon-kenar Bay (FB) in the south Caspian Sea (CS) are examined in order to characterize the shelf motion, the steady current field and to determine the main driving forces of currents on the study area. Owing to closed basin and absence of the astronomical tide, the atmospheric forcing plays an important role in the flow field of the CS. The lasting regular sea breeze system is present almost throughout the year. This system performs the forcing in diurnal and semi-diurnal bands similar to tides in other regions. In general, current field in the continental shelf could be separated into two distinguishable schemes, which in cross-shelf direction is dominated by high frequencies (1 cpd and higher frequencies), and in along-shelf orientation mostly proportional to lower frequencies in synoptic weather bands. Long-period wave currents, whose velocities are much greater than those of direct wind-induced currents, dominates the current field in the continental shelf off FB. The propagation of the latter could be described in terms of shore-controlled waves that are remotely generated and travel across the shelf in the southern CS. It has also been shown that long term displacements in this area follow the classic cyclonic, circulation pattern in the southern CS.

Variability of heat and salinity content in the North Atlantic in the last decade - V. O. Ivchenko; N. C. Wells; D. L. Aleynik; A. G. P. Shaw (Ocean Science, 2010)
The analysis of the heat and salinity contents has been made for the Northern Atlantic for the decade between January 1999 and December 2008. This analysis is based on the Argo profiling data for the upper 2000 m. Basin-averaged values of heat content deviation (HCD) and salinity content deviation (SCD) are robust and stable. The HCD and SCD demonstrate positive trends in the last decade in the upper 2000 m of the North Atlantic. The linear trend of HCD and SCD are (11.14±3.17)×1020 J/yr, and (2.80±1.17)×1013 kg/yr, respectively. Both trends are significant at 95% level of significance.

Numerical implementation and oceanographic application of the thermodynamic potentials of liquid water, water vapour, ice, seawater and humid air – Part 2: The library routines - D. G. Wright; R. Feistel; J. H. Reissmann; K. Miyagawa; D. R. Jackett; W. Wagner; U. Overhoff; C. Guder; A. Feistel; G. M. Marion (Ocean Science, 2010)
The SCOR/IAPSO1 Working Group 127 on Thermodynamics and Equation of State of Seawater has prepared recommendations for new methods and algorithms for numerical estimation of the the thermophysical properties of seawater. As an outcome of this work, a new International Thermodynamic Equation of Seawater (TEOS–10) was endorsed by IOC/UNESCO2 in June 2009 as the official replacement and extension of the 1980 International Equation of State, EOS-80. As part of this new standard a source code package has been prepared that is now made freely available to users via the World Wide Web. This package includes two libraries referred to as the SIA (Sea-Ice-Air) library and the GSW (Gibbs SeaWater) library. Information on the GSW library may be found on the TEOS-10 web site (http://www.TEOS-10.org). This publication provides an introduction to the SIA library which contains routines to calculate various thermodynamic properties as discussed in the companion paper. The SIA library is very comprehensive, including routines to deal with fluid water, ice, seawater and humid air as well as equilibrium states involving various combinations of these, with equivalent code developed in different languages. The code is hierachically structured in modules that support (i) almost unlimited extension with respect to additional properties or relations, (ii) an extraction of self-contained sub-libraries, (iii) separate updating of the empirical thermodynamic potentials, and (iv) code verification on different platforms and between different languages. Error trapping is implemented to identify when one or more of the primary routines are accessed significantly beyond their established range of validity. The initial version of the SIA library is available in Visual Basic and FORTRAN as a supplement to this publication and updates will be maintained on the TEOS-10 web site. 1SCOR/IAPSO: Scientific Committee on Oceanic Research/International Association for the Physical Sciences of the Oceans 2IOC/UNESCO: Intergovernmental Oceanographic Commission/United Nations Educational, Scientific and Cultural Organization

Mixed layer sub-mesoscale parameterization – Part 1: Derivation and assessment - V. M. Canuto; M. S. Dubovikov (Ocean Science, 2010)
Several studies have shown that sub-mesoscales (SM ~1 km horizontal scale) play an important role in mixed layer dynamics. In particular, high resolution simulations have shown that in the case of strong down-front wind, the re-stratification induced by the SM is of the same order of the de-stratification induced by small scale turbulence, as well as of that induced by the Ekman velocity. These studies have further concluded that it has become necessary to include SM in ocean global circulation models (OGCMs), especially those used in climate studies. The goal of our work is to derive and assess an analytic parameterization of the vertical tracer flux under baroclinic instabilities and wind of arbitrary directions and strength. To achieve this goal, we have divided the problem into two parts: first, in this work we derive and assess a parameterization of the SM vertical flux of an arbitrary tracer for ocean codes that resolve mesoscales, M, but not sub-mesoscales, SM. In Part 2, presented elsewhere, we have used the results of this work to derive a parameterization of SM fluxes for ocean codes that do not resolve either M or SM. To carry out the first part of our work, we solve the SM dynamic equations including the non-linear terms for which we employ a closure developed and assessed in previous work. We present a detailed analysis for down-front and up-front winds with the following results: (a) down-front wind (blowing in the direction of the surface geostrophic velocity) is the most favorable condition for generating vigorous SM eddies; the de-stratifying effect of the mean flow and re-stratifying effect of SM almost cancel each other out, (b) in the up-front wind case (blowing in the direction opposite to the surface geostrophic velocity), strong winds prevents the SM generation while weak winds hinder the process but the eddies amplify the re-stratifying effect of the mean velocity, (c) wind orthogonal to the geostrophic velocity. In this case, which was not considered in numerical simulations, we show that when the wind direction coincides with that of the horizontal buoyancy gradient, SM eddies are generated and their re-stratifying effect partly cancels the de-stratifying effect of the mean velocity. The case when wind direction is opposite to that of the horizontal buoyancy gradient, is analogous to the case of up-front winds. In conclusion, the new multifaceted implications on the mixed layer stratification caused by the interplay of both strength and directions of the wind in relation to the buoyancy gradient disclosed by high resolution simulations have been reproduced by the present model. The present results can be used in OGCMs that resolve M but not SM.

Numerical implementation and oceanographic application of the thermodynamic potentials of liquid water, water vapour, ice, seawater and humid air – Part 1: Background and equations - R. Feistel; D. G. Wright; D. R. Jackett; K. Miyagawa; J. H. Reissmann; W. Wagner; U. Overhoff; C. Guder; A. Feistel; G. M. Marion (Ocean Science, 2010)
A new seawater standard referred to as the International Thermodynamic Equation of Seawater 2010 (TEOS-10) was adopted in June 2009 by UNESCO/IOC on its 25th General Assembly in Paris, as recommended by the SCOR/IAPSO Working Group 127 (WG127) on Thermodynamics and Equation of State of Seawater. To support the adoption process, WG127 has developed a comprehensive source code library for the thermodynamic properties of liquid water, water vapour, ice, seawater and humid air, referred to as the Sea-Ice-Air (SIA) library. Here we present the background information and equations required for the determination of the properties of single phases and components as well as of phase transitions and composite systems as implemented in the library. All results are based on rigorous mathematical methods applied to the Primary Standards of the constituents, formulated as empirical thermodynamic potential functions and, except for humid air, endorsed as Releases of the International Association for the Properties of Water and Steam (IAPWS). Details of the implementation in the TEOS-10 SIA library are given in a companion paper.

Seismic imaging of a thermohaline staircase in the western tropical North Atlantic - I. Fer; P. Nandi; W. S. Holbrook; R. W. Schmitt; P. Páramo (Ocean Science, 2010)
Multichannel seismic data acquired in the Lesser Antilles in the western tropical North Atlantic indicate that the seismic reflection method has imaged an oceanic thermohaline staircase. Synthetic acoustic modeling using measured density and sound speed profiles corroborates inferences from the seismic data. In a small portion of the seismic image, laterally coherent, uniform layers are present at depths ranging from 550–700 m and have a separation of ~20 m, with thicknesses increasing with depth. The reflection coefficient, a measure of the acoustic impedance contrasts across these reflective interfaces, is one order of magnitude greater than background noise. Hydrography sampled in previous surveys suggests that the layers are a permanent feature of the region. Spectral analysis of layer horizons in the thermohaline staircase indicates that internal wave activity is anomalously low, suggesting weak internal wave-induced turbulence. Results from two independent measurements, the application of a finescale parameterization to observed high-resolution velocity profiles and direct measurements of turbulent dissipation rate, confirm these low levels of turbulence. The lack of internal wave-induced turbulence may allow for the maintenance of the staircase or may be due to suppression by the double-diffusive convection within the staircase. Our observations show the potential for seismic oceanography to contribute to an improved understanding of occurrence rates and the geographical distribution of thermohaline staircases, and should thereby improve estimates of vertical mixing rates ascribable to salt fingering in the global ocean.

Technical Note: A low cost unmanned aerial vehicle for ship based science missions - E. Waugh; M. Mowlem (Ocean Science, 2010)
A low-cost Unmanned Aerial Vehicle is compared with those already available and the motivation for its development is established. It is targeted at ship-based science missions and potential applications are described including a specific science case to measure white capping in the deep ocean. The current vehicle includes a range of more than 1000 Km, carrying a payload of 2 Kg and it can be launched and recovered from a coastal research vessel. The vehicle has flown successfully in Force 4 gusting Force 6–7 wind conditions, an important requirement for operation at sea. Data analysis is performed on images captured by the vehicle to provide a measurement of wave period and white capping fraction. The next stage of the project is to develop a suitable payload and perform a demonstration science mission.

Anthropogenic carbon dynamics in the changing ocean - J. F. Tjiputra; K. Assmann; C. Heinze (Ocean Science, 2010)
The long-term response of CO2 fluxes to climate change at the ocean surface and within the ocean interior is investigated using a coupled climate-carbon cycle model. This study also presents the first attempt to quantify the evolution of lateral transport of anthropogenic carbon under future climate change. Additionally, its impact on regional carbon storage and uptake is also evaluated. For the 1850–2099 period, our climate change simulation predicts oceanic uptake of anthropogenic carbon of about 538±23 Pg C. Another simulation indicates that changes in physical climate and its associated biogeochemical feedbacks result in a release of natural carbon of about 22±30 Pg C. The natural carbon outgassing is attributed to the reduction in solubility and change in wind pattern in the Southern Hemisphere. After the anthropogenic carbon passes through the air-sea interface, it is predominantly transported along the large scale overturning circulation below the surface layer. The spatial variations in the transport patterns in turn influence the evolution of future regional carbon uptake. In the North Atlantic, a slow down in the Atlantic Meridional Overturning Circulation weakens the penetration strength of anthropogenic carbon into the deeper ocean, which leads to a reduced uptake rate in this region. In contrast, more than half of the anthropogenic carbon taken up in the high latitude Southern Ocean region (south of 58° S) is efficiently and continuously exported northward, predominantly into intermediate waters. This transport mechanism allows continuous increase in future carbon uptake in the high latitude Southern Ocean, where the annual uptake strength could reach 39.3±0.9 g C m−2 yr−1, more than twice the global mean of 16.0±0.3 g C m−2 yr−1 by the end of the 21st century. Our study further underlines the key role of the Southern Ocean in controlling long-term future carbon uptake.

Super-ensemble techniques applied to wave forecast: performance and limitations - F. Lenartz; J.-M. Beckers; J. Chiggiato; B. Mourre; C. Troupin; L. Vandenbulcke; M. Rixen (Ocean Science, 2010)
Nowadays, several operational ocean wave forecasts are available for a same region. These predictions may considerably differ, and to choose the best one is generally a difficult task. The super-ensemble approach, which consists in merging different forecasts and past observations into a single multi-model prediction system, is evaluated in this study. During the DART06 campaigns organized by the NATO Undersea Research Centre, four wave forecasting systems were simultaneously run in the Adriatic Sea, and significant wave height was measured at six stations as well as along the tracks of two remote sensors. This effort provided the necessary data set to compare the skills of various multi-model combination techniques. Our results indicate that a super-ensemble based on the Kalman Filter improves the forecast skills: The bias during both the hindcast and forecast periods is reduced, and the correlation coefficient is similar to that of the best individual model. The spatial extrapolation of local results is not straightforward and requires further investigation to be properly implemented.

Temporal energy partitions of Florida extreme sea level events as a function of Atlantic multidecadal oscillation - J. Park; J. Obeysekera; J. Barnes (Ocean Science, 2010)
An energy-conservative metric based on the discrete wavelet transform is applied to assess the relative energy distribution of extreme sea level events across different temporal scales. The metric is applied to coastal events at Key West and Pensacola Florida as a function of two Atlantic Multidecadal Oscillation (AMO) regimes. Under AMO warm conditions there is a small but significant redistribution of event energy from nearly static into more dynamic (shorter duration) timescales at Key West, while at Pensacola the AMO-dependent changes in temporal event behaviour are less pronounced. Extreme events with increased temporal dynamics might be consistent with an increase in total energy of event forcings which may be a reflection of more energetic storm events during AMO warm phases. As dynamical models mature to the point of providing regional climate index predictability, coastal planners may be able to consider such temporal change metrics in planning scenarios.

Importance of the variability of hydrographic preconditioning for deep convection in the Gulf of Lion, NW Mediterranean - L. Grignon; D. A. Smeed; H. L. Bryden; K. Schroeder (Ocean Science, 2010)
We study the variability of hydrographic preconditioning defined as the heat and salt contents in the Ligurian Sea before convection. The stratification is found to reach a maximum in the intermediate layer in December, whose causes and consequences for the interannual variability of convection are investigated. Further study of the interannual variability and correlation tests between the properties of the deep water formed and the winter surface fluxes support the description of convection as a process that transfers the heat and salt contents from the top and intermediate layers to the deep layer. A proxy for the rate of transfer is given by the final convective mixed layer depth, that is shown to depend equally on the surface fluxes and on the preconditioning. In particular, it is found that deep convection in winter 2004–2005 would have happened even with normal winter conditions, due to low pre-winter stratification.

On the numerical resolution of the bottom layer in simulations of oceanic gravity currents - N. Laanaia; A. Wirth; J. M. Molines; B. Barnier; J. Verron (Ocean Science, 2010)
The role of an increased numerical vertical resolution, leading to an explicit resolution of the bottom Ekman layer dynamics, is investigated. Using the hydrostatic ocean model NEMO-OPA9, we demonstrate that the dynamics of an idealised gravity current (on an inclined plane), is well captured when a few (around five) sigma-coordinate levels are added near the ocean floor. Such resolution allows to considerably improve the representation of the descent and transport of the gravity current and the Ekman dynamics near the ocean floor, including the important effect of Ekman veering, which is usually neglected in today's simulations of the ocean dynamics. Results from high resolution simulations (with σ and z-coordinates) are compared to simulations with a vertical resolution commonly employed in today's ocean models. The latter show a downslope transport that is reduced by almost an order of magnitude and the decrease in the along slope transport is reduced six-fold. We strongly advocate for an increase of the numerical resolution at the ocean floor, similar to the way it is done at the ocean surface and at the lower boundary in atmospheric models.

Automated gas bubble imaging at sea floor – a new method of in situ gas flux quantification - K. Thomanek; O. Zielinski; H. Sahling; G. Bohrmann (Ocean Science, 2010)
Photo-optical systems are common in marine sciences and have been extensively used in coastal and deep-sea research. However, due to technical limitations in the past photo images had to be processed manually or semi-automatically. Recent advances in technology have rapidly improved image recording, storage and processing capabilities which are used in a new concept of automated in situ gas quantification by photo-optical detection. The design for an in situ high-speed image acquisition and automated data processing system is reported ("Bubblemeter"). New strategies have been followed with regards to back-light illumination, bubble extraction, automated image processing and data management. This paper presents the design of the novel method, its validation procedures and calibration experiments. The system will be positioned and recovered from the sea floor using a remotely operated vehicle (ROV). It is able to measure bubble flux rates up to 10 L/min with a maximum error of 33% for worst case conditions. The Bubblemeter has been successfully deployed at a water depth of 1023 m at the Makran accretionary prism offshore Pakistan during a research expedition with R/V Meteor in November 2007.

The role of continental shelves in nitrogen and carbon cycling: Northwestern North Atlantic case study - K. Fennel (Ocean Science, 2010)
Continental shelves play a key role in the cycling of nitrogen and carbon. Here the physical transport and biogeochemical transformation processes affecting the fluxes into and out of continental shelf systems are reviewed, and their role in the global cycling of both elements is discussed. Uncertainties in the magnitude of organic and inorganic matter exchange between shelves and the open ocean is a major source of uncertainty in observation-based estimates of nitrogen and carbon fluxes. The shelf-open ocean exchange is hard to quantify based on observations alone, but can be inferred from biogeochemical models. Model-based nitrogen and carbon budgets are presented for the Northwestern North Atlantic continental shelf. Results indicate that shelves are an important sink for fixed nitrogen and a source of alkalinity, but are not much more efficient in exporting organic carbon to the deep ocean than the adjacent open ocean for the shelf region considered.

Statistical trend analysis and extreme distribution of significant wave height from 1958 to 1999 – an application to the Italian Seas - G. Martucci; S. Carniel; J. Chiggiato; M. Sclavo; P. Lionello; M. B. Galati (Ocean Science, 2010)
The study is a statistical analysis of sea states timeseries derived using the wave model WAM forced by the ERA-40 dataset in selected areas near the Italian coasts. For the period 1 January 1958 to 31 December 1999 the analysis yields: (i) the existence of a negative trend in the annual- and winter-averaged sea state heights; (ii) the existence of a turning-point in late 80's in the annual-averaged trend of sea state heights at a site in the Northern Adriatic Sea; (iii) the overall absence of a significant trend in the annual-averaged mean durations of sea states over thresholds; (iv) the assessment of the extreme values on a time-scale of thousand years. The analysis uses two methods to obtain samples of extremes from the independent sea states: the r-largest annual maxima and the peak-over-threshold. The two methods show statistical differences in retrieving the return values and more generally in describing the significant wave field. The r-largest annual maxima method provides more reliable predictions of the extreme values especially for small return periods (

Exchange across the shelf break at high southern latitudes - J. M. Klinck; M. S. Dinniman (Ocean Science, 2010)
Exchange of water across the Antarctic shelf break has considerable scientific and societal importance due to its effects on circulation and biology of the region, conversion of water masses as part of the global overturning circulation and basal melt of glacial ice and the consequent effect on sea level rise. The focus in this paper is the onshore transport of warm, oceanic Circumpolar Deep Water (CDW); export of dense water from these shelves is equally important, but has been the focus of other recent papers and will not be considered here. A variety of physical mechanisms are described which could play a role in this onshore flux. The relative importance of some processes are evaluated by simple calculations. A numerical model for the Ross Sea continental shelf is used as an example of a more comprehensive evaluation of the details of cross-shelf break exchange. In order for an ocean circulation model to simulate these processes at high southern latitudes, it needs to have high spatial resolution, realistic geometry and bathymetry. Grid spacing smaller than the first baroclinic radius of deformation (a few km) is required to adequately represent the circulation. Because of flow-topography interactions, bathymetry needs to be represented at these same small scales. Atmospheric conditions used to force these circulation models also need to be known at a similar small spatial resolution (a few km) in order to represent orographically controlled winds (coastal jets) and katabatic winds. Significantly, time variability of surface winds strongly influences the structure of the mixed layer. Daily, if not more frequent, surface fluxes must be imposed for a realistic surface mixed layer. Sea ice and ice shelves are important components of the coastal circulation. Ice isolates the ocean from exchange with the atmosphere, especially in the winter. Melting and freezing of both sea ice and glacial ice influence salinity and thereby the character of shelf water. These water mass conversions are known to have an important effect on export of dense water from many Antarctic coastal areas. An artificial dye, as well as temperature, is used to diagnose the flux of CDW onto the shelf. Model results for the Ross Sea show a vigorous onshore flux of oceanic water across the shelf break both at depth and at the surface as well as creation of dense water (High Salinity Shelf Water) created by coastal polynyas in the western Ross Sea.

Biofouling protection for marine environmental sensors - L. Delauney; C. Compère; M. Lehaitre (Ocean Science, 2010)
These days, many marine autonomous environment monitoring networks are set up in the world. These systems take advantage of existing superstructures such as offshore platforms, lightships, piers, breakwaters or are placed on specially designed buoys or underwater oceanographic structures. These systems commonly use various sensors to measure parameters such as dissolved oxygen, turbidity, conductivity, pH or fluorescence. Emphasis has to be put on the long term quality of measurements, yet sensors may face very short-term biofouling effects. Biofouling can disrupt the quality of the measurements, sometimes in less than a week. Many techniques to prevent biofouling on instrumentation are listed and studied by researchers and manufacturers. Very few of them are implemented on instruments and of those very few have been tested in situ on oceanographic sensors for deployment of at least one or two months. This paper presents a review of techniques used to protect against biofouling of in situ sensors and gives a short list and description of promising techniques.

Long term trends in the sea surface temperature of the Black Sea - G. I. Shapiro; D. L. Aleynik; L. D. Mee (Ocean Science, 2010)
There is growing understanding that recent deterioration of the Black Sea ecosystem was partly due to changes in the marine physical environment. This study uses high resolution 0.25° climatology to analyze sea surface temperature variability over the 20th century in two contrasting regions of the sea. Results show that the deep Black Sea was cooling during the first three quarters of the century and was warming in the last 15–20 years; on aggregate there was a statistically significant cooling trend. The SST variability over the Western shelf was more volatile and it does not show statistically significant trends. The cooling of the deep Black Sea is at variance with the general trend in the North Atlantic and may be related to the decrease of westerly winds over the Black Sea, and a greater influence of the Siberian anticyclone. The timing of the changeover from cooling to warming coincides with the regime shift in the Black Sea ecosystem.

Sensitivity of oxygen dynamics in the water column of the Baltic Sea to external forcing - S. Miladinova; A. Stips (Ocean Science, 2010)
A 1-D biogeochemical/physical model of marine systems has been applied to study the oxygen cycle in four stations of different sub-basins of the Baltic Sea, namely, in the Gotland Deep, Bornholm, Arkona and Fladen. The model consists of the biogeochemical model of Neumann et al. (2002) coupled with the 1-D General Ocean Turbulence Model (GOTM). The model has been forced with meteorological data from the ECMWF reanalysis project for the period 1998–2003, producing a six year hindcast which is validated with datasets from the Baltic Environmental Database (BED) for the same period. The vertical profiles of temperature and salinity are relaxed towards both profiles provided by 3-D simulations of General Estuarine Transport Model (GETM) and observed profiles from BED. Modifications in the parameterisation of the air-sea oxygen fluxes have led to a significant improvement of the model results in the surface and intermediate water layers. The largest mismatch with observations is found in simulating the oxygen dynamics in the Baltic Sea bottom waters. The model results demonstrate the good capability of the model to predict the time-evolution of the physical and biogeochemical variables at all different stations. Comparative analysis of the modelled oxygen concentrations with respect to observation data is performed to distinguish the relative importance of several factors on the seasonal, interannual and long-term variations of oxygen. It is found that natural physical factors, like the magnitude of the vertical turbulent mixing, wind speed and the variation of temperature and salinity fields are the major factors controlling the oxygen dynamics in the Baltic Sea. The influence of limiting nutrients is less pronounced, at least under the nutrient flux parameterisation assumed in the model.

The contribution of eastern-boundary density variations to the Atlantic meridional overturning circulation at 26.5° N - M. P. Chidichimo; T. Kanzow; S. A. Cunningham; W. E. Johns; J. Marotzke (Ocean Science, 2010)
We study the contribution of eastern-boundary density variations to sub-seasonal and seasonal anomalies of the strength and vertical structure of the Atlantic Meridional Overturning Circulation (AMOC) at 26.5° N, by means of the RAPID/MOCHA mooring array between April 2004 and October 2007. The major density anomalies are found in the upper 500 m, and they are often coherent down to 1400 m. The densities have 13-day fluctuations that are apparent down to 3500 m. The two strategies for measuring eastern-boundary density – a tall offshore mooring (EB1) and an array of moorings on the continental slope (EBH) – show little correspondence in terms of amplitude, vertical structure, and frequency distribution of the resulting basin-wide integrated transport fluctuations, implying that there are significant transport contributions between EB1 and EBH. Contrary to the original planning, measurements from EB1 cannot serve as backup or replacement for EBH: density needs to be measured directly at the continental slope to compute the full-basin density gradient. Fluctuations in density at EBH generate transport variability of 2 Sv rms in the AMOC, while the overall AMOC variability is 4.8 Sv rms. There is a pronounced deep-reaching seasonal cycle in density at the eastern boundary, which is apparent between 100 m and 1400 m, with maximum positive anomalies in spring and maximum negative anomalies in autumn. These changes drive anomalous southward upper mid-ocean flow in spring, implying maximum reduction of the AMOC, and vice-versa in autumn. The amplitude of the seasonal cycle of the AMOC arising from the eastern-boundary densities is 5.2 Sv peak-to-peak, dominating the 6.7 Sv peak-to-peak seasonal cycle of the total AMOC. Our analysis suggests that the seasonal cycle in density may be forced by the strong near-coastal seasonal cycle in wind stress curl.

The role of atmosphere and ocean physical processes in ENSO in a perturbed physics coupled climate model - S. Y. Philip; M. Collins; G. J. van Oldenborgh; B. J. J. M. van den Hurk (Ocean Science, 2010)
We examine the behaviour of the El Niño – Southern Oscillation (ENSO) in an ensemble of global climate model simulations with perturbations to parameters in the atmosphere and ocean components respectively. The influence of the uncertainty in these parametrisations on ENSO are investigated systematically. The ensemble exhibits a range of different ENSO behaviour in terms of the amplitude and spatial structure of the sea surface temperature (SST) variability. The nature of the individual feedbacks that operate within the ENSO system are diagnosed using an Intermediate Complexity Model (ICM), which has been used previously to examine the diverse ENSO behaviour of the CMIP3 multi-model ensemble. Unlike in that case, the ENSO in these perturbed physics experiments is not principally controlled by variations in the mean climate state. Rather the parameter perturbations influence the ENSO characteristics by modifying the coupling feedbacks within the cycle. The associated feedbacks that contribute most to the ensemble variations are the response of SST to local wind variability and damping, followed by the response of SST to thermocline anomalies and the response of the zonal wind stress to those SST anomalies. Atmospheric noise amplitudes and oceanic processes play a relatively minor role.

First images and orientation of fine structure from a 3-D seismic oceanography data set - T. M. Blacic; W. S. Holbrook (Ocean Science, 2010)
We present 3-D images of ocean fine structure from a unique industry-collected 3-D multichannel seismic dataset from the Gulf of Mexico that includes expendable bathythermograph casts for both swaths. 2-D processing reveals strong laterally continuous reflections throughout the upper ~800 m as well as a few weaker but still distinct reflections as deep as ~1100 m. We interpret the reflections to be caused by reversible fine structure from internal wave strains. Two bright reflections are traced across the 225-m-wide swath to produce reflection surface images that illustrate the 3-D nature of ocean fine structure. We show that the orientation of linear features in a reflection can be obtained by calculating the orientations of contours of reflection relief, or more robustly, by fitting a sinusoidal surface to the reflection. Preliminary 3-D processing further illustrates the potential of 3-D seismic data in interpreting images of oceanic features such as internal wave strains. This work demonstrates the viability of imaging oceanic fine structure in 3-D and shows that, beyond simply providing a way visualize oceanic fine structure, quantitative information such as the spatial orientation of features like fronts and solitons can be obtained from 3-D seismic images. We expect complete, optimized 3-D processing to improve both the signal to noise ratio and spatial resolution of our images resulting in increased options for analysis and interpretation.

Stochastic Heterogeneity Mapping around a Mediterranean salt lens - G. G. Buffett; C. A. Hurich; E. A. Vsemirnova; R. W. Hobbs; V. Sallarès; R. Carbonell; D. Klaeschen; B. Biescas (Ocean Science, 2010)
We present the first application of Stochastic Heterogeneity Mapping based on the band-limited von Kármán function to a seismic reflection stack of a Mediterranean water eddy (meddy), a large salt lens of Mediterranean water. This process extracts two stochastic parameters directly from the reflectivity field of the seismic data: the Hurst number, which ranges from 0 to 1, and the correlation length (scale length). Lower Hurst numbers represent a richer range of high wavenumbers and correspond to a broader range of heterogeneity in reflection events. The Hurst number estimate for the top of the meddy (0.39) compares well with recent theoretical work, which required values between 0.25 and 0.5 to model internal wave surfaces in open ocean conditions based on simulating a Garrett-Munk spectrum (GM76) slope of −2. The scale lengths obtained do not fit as well to seismic reflection events as those used in other studies to model internal waves. We suggest two explanations for this discrepancy: (1) due to the fact that the stochastic parameters are derived from the reflectivity field rather than the impedance field the estimated scale lengths may be underestimated, as has been reported; and (2) because the meddy seismic image is a two-dimensional slice of a complex and dynamic three-dimensional object, the derived scale lengths are biased to the direction of flow. Nonetheless, varying stochastic parameters, which correspond to different spectral slopes in the Garrett-Munk spectrum (horizontal wavenumber spectrum), can provide an estimate of different internal wave scales from seismic data alone. We hence introduce Stochastic Heterogeneity Mapping as a novel tool in physical oceanography.

Modeling the effects of size on patch dynamics of an inert tracer - P. Xiu; F. Chai (Ocean Science, 2010)
Mesoscale iron enrichment experiments have revealed that additional iron affects the phytoplankton productivity and carbon cycle. However, the role of initial size of fertilized patch in determining the patch evolution is poorly quantified due to the limited observational capability and complex of physical processes. Using a three-dimensional ocean circulation model, we simulated different sizes of inert tracer patches that were only regulated by physical circulation and diffusion. Model results showed that during the first few days since release of inert tracer, the calculated dilution rate was found to be a linear function with time, which was sensitive to the initial patch size with steeper slope for smaller size patch. After the initial phase of rapid decay, the relationship between dilution rate and time became an exponential function, which was also size dependent. Therefore, larger initial size patches can usually last longer and ultimately affect biogeochemical processes much stronger than smaller patches.

Structure and forcing of the overflow at the Storfjorden sill and its connection to the Arctic coastal polynya in Storfjorden - F. Geyer; I. Fer; L. H. Smedsrud (Ocean Science, 2010)
Storfjorden (Svalbard) is a sill-fjord with an active polynya and exemplifies the dense water formation process over the Arctic shelves. Here we report on our simulations of Storfjorden covering the freezing season of 1999–2000 using an eddy-permitting 3-D ocean circulation model with a fully coupled dynamic and thermodynamic sea-ice model. The model results in the polynya region and of the dense water plume flowing over the sill crest are compared to observations. The connections of the overflow at the sill to the dense water production at the polynya and to the local wind forcing are investigated. Both the overflow and the polynya dynamics are found to be sensitive to wind forcing. In response to freezing and brine rejection over the polynya, the buoyancy forcing initiates an abrupt positive density anomaly. While the ocean integrates the buoyancy forcing over several polynya events (about 25 days), the wind forcing dominates the overflow response at the sill at weather scale. In the model, the density excess is diluted in the basin and leads to a gradual build-up of dense water behind the sill. The overflow transport is typically inferred from observations using a single current profiler at the sill crest. Despite the significant variability of the plume width, we show that a constant overflow width of 15 km produces realistic estimates of the overflow volume transport. Another difficulty in monitoring the overflow is measuring the plume thickness in the absence of hydrographic profiles. Volume flux estimates assuming a constant plume width and the thickness inferred from velocity profiles explain 58% of the modelled overflow volume flux variance and agrees to within 10% when averaged over the overflow season.

Sediment 231 Pa/ 230 Th as a recorder of the rate of the Atlantic meridional overturning circulation: insights from a 2-D model - Y. Luo; R. Francois; S. E. Allen (Ocean Science, 2010)
A two dimensional scavenging model is used to investigate the patterns of sediment 231Pa/230Th generated by the Atlantic Meridional Overturning Circulation (AMOC) and further advance the application of this proxy for ocean paleocirculation studies. The scavenging parameters and the geometry of the overturning circulation cell have been chosen so that the model generates meridional sections of dissolved 230Th and 231Pa consistent with published water column profiles and an additional 12 previously unpublished profiles measured in the North and Equatorial Atlantic. The processes that generate the meridional sections of dissolved and particulate 230Th, dissolved and particulate 231Pa, dissolved and particulate 231Pa/230Th, and sediment 231Pa/230Th are discussed in detail. The results indicate that the relationship between sediment 231Pa/230Th at any given site and the overturning circulation is very complex. They clearly show that constraining past changes in the strength and geometry of the AMOC requires an extensive data set and they suggest strategies to maximize information from a limited number of samples.

Corrigendum to "The density of seawater as a function of salinity (5 to 70 g kg−1) and temperature (273.15 to 363.15 K)" published in Ocean Sci., 5, 91–100, 2009 - F. J. Millero; F. Huang (Ocean Science, 2010)

A model for predicting changes in the electrical conductivity, practical salinity, and absolute salinity of seawater due to variations in relative chemical composition - R. Pawlowicz (Ocean Science, 2010)
Salinity determination in seawater has been carried out for almost 30 years using the Practical Salinity Scale 1978. However, the numerical value of so-called practical salinity, computed from electrical conductivity, differs slightly from the true or absolute salinity, defined as the mass of dissolved solids per unit mass of seawater. The difference arises because more recent knowledge about the composition of seawater is not reflected in the definition of practical salinity, which was chosen to maintain historical continuity with previous measures, and because of spatial and temporal variations in the relative composition of seawater. Accounting for these spatial variations in density calculations requires the calculation of a correction factor δSA, which is known to range from 0 to 0.03 g kg−1 in the world oceans. Here a mathematical model relating compositional perturbations to δSA is developed, by combining a chemical model for the composition of seawater with a mathematical model for predicting the conductivity of multi-component aqueous solutions. Model calculations for this estimate of δSA, denoted δSRsoln, generally agree with estimates of δSA based on fits to direct density measurements, denoted δSRdens, and show that biogeochemical perturbations affect conductivity only weakly. However, small systematic differences between model and density-based estimates remain. These may arise for several reasons, including uncertainty about the biogeochemical processes involved in the increase in Total Alkalinity in the North Pacific, uncertainty in the carbon content of IAPSO standard seawater, and uncertainty about the haline contraction coefficient for the constituents involved in biogeochemical processes. This model may then be important in constraining these processes, as well as in future efforts to improve parameterizations for δSA.

Malvinas-slope water intrusions on the northern Patagonia continental shelf - A. R. Piola; N. Martínez Avellaneda; R. A. Guerrero; F. P. Jardón; E. D. Palma; S. I. Romero (Ocean Science, 2010)
The Patagonia continental shelf located off southeastern South America is bounded offshore by the Malvinas Current, which extends northward from northern Drake Passage (~55° S) to nearly 38° S. The transition between relatively warm-fresh shelf waters and Subantarctic Waters from the western boundary current is characterized by a thermohaline front extending nearly 2500 km. We use satellite derived sea surface temperature, and chlorophyll-a data combined with hydrographic and surface drifter data to document the intrusions of slope waters onto the continental shelf near 41° S. These intrusions create vertically coherent localized negative temperature and positive salinity anomalies extending onshore about 150 km from the shelf break. The region is associated with a center of action of the first mode of non-seasonal sea surface temperature variability and also relatively high chlorophyll-a variability, suggesting that the intrusions are important in promoting the local development of phytoplankton. The generation of slope water penetrations at this location may be triggered by the inshore excursion of the 100 m isobath, which appears to steer the Malvinas Current waters over the outer shelf.

Seasonal variability of the Caspian Sea three-dimensional circulation, sea level and air-sea interaction - R. A. Ibrayev; E. Özsoy; C. Schrum; H. İ. Sur (Ocean Science, 2010)
A three-dimensional primitive equation model including sea ice thermodynamics and air-sea interaction is used to study seasonal circulation and water mass variability in the Caspian Sea under the influence of realistic mass, momentum and heat fluxes. River discharges, precipitation, radiation and wind stress are seasonally specified in the model, based on available data sets. The evaporation rate, sensible and latent heat fluxes at the sea surface are computed interactively through an atmospheric boundary layer sub-model, using the ECMWF-ERA15 re-analysis atmospheric data and model generated sea surface temperature. The model successfully simulates sea-level changes and baroclinic circulation/mixing features with forcing specified for a selected year. The results suggest that the seasonal cycle of wind stress is crucial in producing basin circulation. Seasonal cycle of sea surface currents presents three types: cyclonic gyres in December–January; Eckman south-, south-westward drift in February–July embedded by western and eastern southward coastal currents and transition type in August–November. Western and eastern northward sub-surface coastal currents being a result of coastal local dynamics at the same time play an important role in meridional redistribution of water masses. An important part of the work is the simulation of sea surface topography, yielding verifiable results in terms of sea level. The model successfully reproduces sea level variability for four coastal points, where the observed data are available. Analyses of heat and water budgets confirm climatologic estimates of heat and moisture fluxes at the sea surface. Experiments performed with variations in external forcing suggest a sensitive response of the circulation and the water budget to atmospheric and river forcing.

Mediterranean intermediate circulation estimated from Argo data in 2003–2010 - M. Menna; P. M. Poulain (Ocean Science, 2010)
Data from 38 Argo profiling floats are used to describe the intermediate Mediterranean currents for the period October 2003–January 2010. These floats were programmed to execute 5-day cycles, to drift at a neutral parking depth of 350 m and measure temperature and salinity profiles from either 700 or 2000 m up to the surface. At the end of each cycle the floats remained at the sea surface for about 6 h, enough time to be localised and transmit the data to the Argos satellite system. The Argos positions were used to determine the float surface and intermediate displacements. At the surface, the float motion was approximated by a linear displacement and inertial motion. Intermediate velocities estimates were used to investigate the Mediterranean circulation at 350 m, to compute the pseudo-Eulerian statistics and to study the influence of bathymetry on the intermediate currents. Maximum speeds, as large as 33 cm/s, were found northeast of the Balearic Islands (western basin) and in the Ierapetra eddy (eastern basin). Typical speeds in the main along-slope currents (Liguro-Provençal-Catalan, Algerian and Libyo-Egyptian Currents) were ~20 cm/s. In the central and western part of Mediterranean basin, the pseudo-Eulerian statistics show typical intermediate circulation pathways which can be related to the motion of Levantine Intermediate Water. In general our results agree with the qualitative intermediate circulation schemes proposed in the literature, except in the southern Ionian where we found westward-flowing intermediate currents. Fluctuating currents appeared to be usually larger than the mean flow. Intermediate currents were found to be essentially parallel to the isobaths over most of the areas characterized by strong bathymetry gradients, in particular, in the vicinity of the continental slopes.

Entrainment in the Denmark Strait overflow plume by meso-scale eddies - G. Voet; D. Quadfasel (Ocean Science, 2010)
The entrainment of buoyant ambient water into the overflow plume of Denmark Strait and the associated downstream warming of the plume are estimated using time series of currents and temperature from moored instrumentation and classical hydrographic data. Warming rates are highest (0.4–0.5 K/100 km) within the first 200 km of the sill, and decrease to 0.05–0.1 K/100 km further downstream. Stirring by mesoscale eddies causes lateral heat fluxes that explain the 0.1 K/100 km warming, but in the first 200 km from the sill also vertical diapycnal fluxes, probably caused by breaking internal waves, must contribute to the entrainment.

On contribution of horizontal and intra-layer convection to the formation of the Baltic Sea cold intermediate layer - I. Chubarenko; N. Demchenko (Ocean Science, 2010)
Seasonal cascades down the coastal slopes and intra-layer convection are considered as the two additional mechanisms contributing to the Baltic Sea cold intermediate layer (CIL) formation along with conventional seasonal vertical mixing. Field measurements are presented, reporting for the first time the possibility of denser water formation and cascading from the Baltic Sea underwater slopes, which take place under fall and winter cooling conditions and deliver waters into intermediate layer of salinity stratified deep-sea area. The presence in spring within the CIL of water with temperature below that of maximum density (Tmd) and that at the local surface in winter time allows tracing its formation: it is argued that the source of the coldest waters of the Baltic CIL is early spring (March–April) cascading, arising due to heating of water before reaching the Tmd. Fast increase of the open water heat content during further spring heating indicates that horizontal exchange rather than direct solar heating is responsible for that. When the surface is covered with water, heated above the Tmd, the conditions within the CIL become favorable for intralayer convection due to the presence of waters of Tmd in intermediate layer, which can explain its well-known features – the observed increase of its salinity and deepening with time.

Impact of global ocean model resolution on sea-level variability with emphasis on interannual time scales - T. Penduff; M. Juza; L. Brodeau; G. C. Smith; B. Barnier; J.-M. Molines; A.-M. Treguier; G. Madec (Ocean Science, 2010)
Four global ocean/sea-ice simulations driven by the same realistic 47-year daily atmospheric forcing were performed by the DRAKKAR group at 2°, 1°, &frac12°, and ¼° resolutions. Simulated mean sea-surface heights (MSSH) and sea-level anomalies (SLA) are collocated over the period 1993–2004 onto the AVISO dataset. MSSH fields are compared with an inverse estimate. SLA datasets are filtered and compared over various time and space scales with AVISO regarding three characteristics: SLA standard deviations, spatial correlations between SLA variability maps, and temporal correlations between observed and simulated band-passed filtered local SLA timeseries. Beyond the 2°−1° transition whose benefits are moderate, further increases in resolution and associated changes in subgrid scale parameterizations simultaneously induce (i) strong increases in SLA standard deviations, (ii) strong improvements in the spatial distribution of SLA variability, and (iii) slight decreases in temporal correlations between observed and simulation SLA timeseries. These 3 effects are not only clear on mesoscale (14–180 days) and quasi-annual (5–18 months) fluctuations, but also on the slower (interannual), large-scale variability ultimately involved in ocean-atmosphere coupled processes. Most SLA characteristics are monotonically affected by successive resolution increases, but irregularly and with a strong dependance on frequency and latitude. Benefits of enhanced resolution are greatest in the 1°−½° and ½°−¼° transitions, in the 14–180 day range, and within eddy-active mid- and high-latitude regions. In the real ocean, most eddy-active areas are characterized by a strong SLA variability at all timescales considered here; this localized, broad-banded temporal variability is only captured at ¼° resolution.

Preface "Operational oceanography in the Mediterranean Sea: the second stage of development" - N. Pinardi; G. Coppini (Ocean Science, 2010)

Characterization of mixing errors in a coupled physical biogeochemical model of the North Atlantic: implications for nonlinear estimation using Gaussian anamorphosis - D. Béal; P. Brasseur; J.-M. Brankart; Y. Ourmières; J. Verron (Ocean Science, 2010)
In biogeochemical models coupled to ocean circulation models, vertical mixing is an important physical process which governs the nutrient supply and the plankton residence in the euphotic layer. However, vertical mixing is often poorly represented in numerical simulations because of approximate parameterizations of sub-grid scale turbulence, wind forcing errors and other mis-represented processes such as restratification by mesoscale eddies. Getting a sufficient knowledge of the nature and structure of these errors is necessary to implement appropriate data assimilation methods and to evaluate if they can be controlled by a given observation system. In this paper, Monte Carlo simulations are conducted to study mixing errors induced by approximate wind forcings in a three-dimensional coupled physical-biogeochemical model of the North Atlantic with a 1/4° horizontal resolution. An ensemble forecast involving 200 members is performed during the 1998 spring bloom, by prescribing perturbations of the wind forcing to generate mixing errors. The biogeochemical response is shown to be rather complex because of nonlinearities and threshold effects in the coupled model. The response of the surface phytoplankton depends on the region of interest and is particularly sensitive to the local stratification. In addition, the statistical relationships computed between the various physical and biogeochemical variables reflect the signature of the non-Gaussian behaviour of the system. It is shown that significant information on the ecosystem can be retrieved from observations of chlorophyll concentration or sea surface temperature if a simple nonlinear change of variables (anamorphosis) is performed by mapping separately and locally the ensemble percentiles of the distributions of each state variable on the Gaussian percentiles. The results of idealized observational updates (performed with perfect observations and neglecting horizontal correlations) indicate that the implementation of this anamorphosis method into sequential assimilation schemes can substantially improve the accuracy of the estimation with respect to classical computations based on the Gaussian assumption.

Assessment of sensor performance - C. Waldmann; M. Tamburri; R. D. Prien; P. Fietzek (Ocean Science, 2010)
There is an international commitment to develop a comprehensive, coordinated and sustained ocean observation system. However, a foundation for any observing, monitoring or research effort is effective and reliable in situ sensor technologies that accurately measure key environmental parameters. Ultimately, the data used for modelling efforts, management decisions and rapid responses to ocean hazards are only as good as the instruments that collect them. There is also a compelling need to develop and incorporate new or novel technologies to improve all aspects of existing observing systems and meet various emerging challenges. Assessment of Sensor Performance was a cross-cutting issues session at the international OceanSensors08 workshop in Warnemünde, Germany, which also has penetrated some of the papers published as a result of the workshop (Denuault, 2009; Kröger et al., 2009; Zielinski et al., 2009). The discussions were focused on how best to classify and validate the instruments required for effective and reliable ocean observations and research. The following is a summary of the discussions and conclusions drawn from this workshop, which specifically addresses the characterisation of sensor systems, technology readiness levels, verification of sensor performance and quality management of sensor systems.

Heat in the Barents Sea: transport, storage, and surface fluxes - L. H. Smedsrud; R. Ingvaldsen; J. E. Ø. Nilsen; Ø. Skagseth (Ocean Science, 2010)
A column model is set up for the Barents Sea to explore sensitivity of surface fluxes and heat storage from varying ocean heat transport. Mean monthly ocean transport and atmospheric forcing are synthesised and force the simulations. Results show that by using updated ocean transports of heat and freshwater the vertical mean hydrographic seasonal cycle can be reproduced fairly well. Our results indicate that the ~70 TW of heat transported to the Barents Sea by ocean currents is lost in the southern Barents Sea as latent, sensible, and long wave radiation, each contributing 23–39 TW to the total heat loss. Solar radiation adds 26 TW in the south, as there is no significant ice production. The northern Barents Sea receives little ocean heat transport. This leads to a mixed layer at the freezing point during winter and significant ice production. There is little net surface heat loss annually in the north. The balance is achieved by a heat loss through long wave radiation all year, removing most of the summer solar heating. During the last decade the Barents Sea has experienced an atmospheric warming and an increased ocean heat transport. The Barents Sea responds to such large changes by adjusting temperature and heat loss. Decreasing the ocean heat transport below 50 TW starts a transition towards Arctic conditions. The heat loss in the Barents Sea depend on the effective area for cooling, and an increased heat transport leads to a spreading of warm water further north.

Surface expression of Mediterranean Water dipoles and their contribution to the shelf/slope – open ocean exchange - N. Serra; I. Ambar; D. Boutov (Ocean Science, 2010)
The generation of dipolar eddies within the Mediterranean Water (MW) layers, at the Portuguese continental slope, was observed using subsurface RAFOS floats. The surface expression of these mid-depth dipoles is here characterized with remote sensing data, namely with sea surface temperature, chlorophyll concentration and sea surface height measurements. Two cases are presented demonstrating the remote detection of these underwater structures in the Gulf of Cadiz during 1998 and 2001. The presence of subsurface MW dipoles in the Iberian coastal zone is shown to influence the development of coastal upwelling filaments. The surface circulation induced by the dipoles causes the upwelling filaments to extend offshorewards and thus to enhance the transport of physical, chemical and biological properties into the open ocean. A numerical model simulation of the ocean circulation around the Iberian Peninsula forced by heat and freshwater fluxes (computed using the NCEP reanalysis atmospheric state) and by the overflow of MW at the Strait of Gibraltar, corroborates the connection between the surface and the mid-depth flows. The high-resolution numerical experiment is used to help clarifying the occurrence of the MW dipoles surface expression and the impact of these dipoles on the eddy kinetic energy of the upper ocean and on the exchange of volume and salt between the shelf/slope and the open ocean.

Variability of scaling time series in the Arctic sea-ice drift dynamics - A. Chmel; V. N. Smirnov; I. B. Sheikin (Ocean Science, 2010)
The motion of an individual ice floe in the Arctic Ocean was monitored at the Russian research station North Pole 35 established on the ice pack in 2008. The ice floe speed (V) was found to be correlated with wind speed (v) in main features, such as the positions of maxima and minima of V and v. However, the fine structure of the V-variation cannot be explained by the wind forcing alone. There were periods of time when the floe drift was affected by the interactions of ice floes between each other or by the periodical forcing due to either the Coriolis inertia effect or the tidal activity. These data were compared with the "waiting times" statistics that are the distributions of time intervals between subsequent, sufficiently strong changes in the kinetic energy of drifting ice floe. These distributions were measured in several time windows differing in the average wind speed and wind direction, and/or in the mechanical state of the ice pack. The distribution functions N (t>τ), where N is the number of successive events of energy change separated by the time interval t that exceeds τ, constructed in different time windows demonstrate fractal or a multifractal nature of the time series during motion in the consolidated ice pack but were truly random when the ice floe drifted in the highly fragmented sea ice. The latter result shows the existence of a relationship between the long-range mechanical interactions in the pack and long-term memory (time scaling behaviour) of the sea-ice motion.

The gyre-scale circulation of the North Atlantic and sea level at Brest - P. L. Woodworth; N. Pouvreau; G. Wöppelmann (Ocean Science, 2010)
The relationship between the gyre-scale circulation of the North Atlantic, represented by air pressure near to the centre of the sub-tropical gyre, and sea level measured at the eastern boundary of the ocean has been investigated using records commencing in the middle of the 18th century. These time series are twice as long as those employed in an earlier study of this relationship. Near-continuous values of annual mean sea level and mean high water from Brest, and air pressure fields for the eastern North Atlantic derived from terrestrial instrumental pressure records and ship logbook information, have been used to demonstrate that sea level on the eastern boundary does indeed appear to be related to air pressure at the centre of the gyre (subject to reservations concerning short sub-sections of data near to the ends of the records). These findings confirm the earlier conclusions but over much longer timescales. This relationship can explain at least part of the century timescale accelerations in European sea level records obtained from tide gauge and saltmarsh data. This finding has important implications for interpretation of the observed sea level rise and acceleration on the European Atlantic coast, suggesting that redistribution of water could play an important role instead of (or as well as) change in ocean volume.

The timescale and extent of thermal expansion of the global ocean due to climate change - S. Marčelja (Ocean Science, 2010)
With recently improved instrumental accuracy, the change in the heat content of the oceans and the corresponding contribution to the change of the sea level can be determined from in situ measurements of temperature variation with depth. Nevertheless, it would be favourable if the same changes could be evaluated from just the sea surface temperatures because the past record could then be reconstructed and future scenarios explored. Using a single column model we show that the average change in the heat content of the oceans and the corresponding contribution to a global change in the sea level can be evaluated from the past sea surface temperatures. The calculation is based on the time-dependent diffusion equation with the known fixed average upwelling velocity and eddy diffusivity, as determined from the steady-state limit. In this limit, the model reduces to the 1966 Munk profile of the potential temperature variation as a function of depth. There are no adjustable parameters in the calculation and the results are in good agreement with the estimates obtained from the in situ data. The method allows us to obtain relevant timescales and average temperature profiles. The evaluation of the thermosteric sea level change is extended back to the beginning of accurate sea surface temperature records. The changes in sea surface temperature from 1880 until the present time are estimated to have produced a thermosteric sea level rise of 35 mm. Application to future IPCC scenarios gives results similar to the average prediction of more complex climate models.

Ensemble perturbation smoother for optimizing tidal boundary conditions by assimilation of High-Frequency radar surface currents – application to the German Bight - A. Barth; A. Alvera-Azcárate; K.-W. Gurgel; J. Staneva; A. Port; J.-M. Beckers; E. V. Stanev (Ocean Science, 2010)
High-Frequency (HF) radars measure the ocean surface currents at various spatial and temporal scales. These include tidal currents, wind-driven circulation, density-driven circulation and Stokes drift. Sequential assimilation methods updating the model state have been proven successful to correct the density-driven currents by assimilation of observations such as sea surface height, sea surface temperature and in-situ profiles. However, the situation is different for tides in coastal models since these are not generated within the domain, but are rather propagated inside the domain through the boundary conditions. For improving the modeled tidal variability it is therefore not sufficient to update the model state via data assimilation without updating the boundary conditions. The optimization of boundary conditions to match observations inside the domain is traditionally achieved through variational assimilation methods. In this work we present an ensemble smoother to improve the tidal boundary values so that the model represents more closely the observed currents. To create an ensemble of dynamically realistic boundary conditions, a cost function is formulated which is directly related to the probability of each boundary condition perturbation. This cost function ensures that the boundary condition perturbations are spatially smooth and that the structure of the perturbations satisfies approximately the harmonic linearized shallow water equations. Based on those perturbations an ensemble simulation is carried out using the full three-dimensional General Estuarine Ocean Model (GETM). Optimized boundary values are obtained by assimilating all observations using the covariances of the ensemble simulation.

Transformation of an Agulhas eddy near the continental slope - S. Baker-Yeboah; G. R. Flierl; G. G. Sutyrin; Y. Zhang (Ocean Science, 2010)
The transformation of Agulhas eddies near the continental slope of southern Africa and their subsequent self-propagation are analyzed in both observational data and numerical simulations. Self-propagation results from a net dipole moment of a generalized heton structure consisting of a surface-intensified anticyclonic eddy and deep cyclonic pattern. Such Agulhas vortical structures can form near the retroflection region and further north along the western coast of southern Africa. We analyze nonlinear topographic wave generation, vortex deformations, and filament production as an important part in water mass exchange. Self-propagating structures provide a conduit for exchange between the deep ocean and shelf regions in the Benguela upwelling system.

Thermodynamic properties of sea air - R. Feistel; D. G. Wright; H.-J. Kretzschmar; E. Hagen; S. Herrmann; R. Span (Ocean Science, 2010)
Very accurate thermodynamic potential functions are available for fluid water, ice, seawater and humid air covering wide ranges of temperature and pressure conditions. They permit the consistent computation of all equilibrium properties as, for example, required for coupled atmosphere-ocean models or the analysis of observational or experimental data. With the exception of humid air, these potential functions are already formulated as international standards released by the International Association for the Properties of Water and Steam (IAPWS), and have been adopted in 2009 for oceanography by IOC/UNESCO. In this paper, we derive a collection of formulas for important quantities expressed in terms of the thermodynamic potentials, valid for typical phase transitions and composite systems of humid air and water/ice/seawater. Particular attention is given to equilibria between seawater and humid air, referred to as "sea air" here. In a related initiative, these formulas will soon be implemented in a source-code library for easy practical use. The library is primarily aimed at oceanographic applications but will be relevant to air-sea interaction and meteorology as well. The formulas provided are valid for any consistent set of suitable thermodynamic potential functions. Here we adopt potential functions from previous publications in which they are constructed from theoretical laws and empirical data; they are briefly summarized in the appendix. The formulas make use of the full accuracy of these thermodynamic potentials, without additional approximations or empirical coefficients. They are expressed in the temperature scale ITS-90 and the 2008 Reference-Composition Salinity Scale.

Spatiotemporal variations of f CO 2 in the North Sea - A. M. Omar; A. Olsen; T. Johannessen; M. Hoppema; H. Thomas; A. V. Borges (Ocean Science, 2010)
Data from two Voluntary Observing Ship (VOS) (2005–2007) augmented with data subsets from ten cruises (1987–2005) were used to investigate the spatiotemporal variations of the CO2 fugacity in seawater (fCO2sw) in the North Sea at seasonal and inter-annual time scales. The observed seasonal fCO2sw variations were related to variations in sea surface temperature (SST), biology plus mixing, and air-sea CO2 exchange. Over the study period, the seasonal amplitude in fCO2sw induced by SST changes was 0.4–0.75 times those resulting from variations in biology plus mixing. Along a meridional transect, fCO2sw normally decreased northwards (−12 μatm per degree latitude), but the gradient disappeared/reversed during spring as a consequence of an enhanced seasonal amplitude of fCO2sw in southern parts of the North Sea. Along a zonal transect, a weak gradient (−0.8 μatm per degree longitude) was observed in the annual mean fCO2sw. Annually and averaged over the study area, surface waters of the North Sea were CO2 undersaturated and, thus, a sink of atmospheric CO2. However, during summer, surface waters in the region 55.5–54.5° N were CO2 supersaturated and, hence, a source for atmospheric CO2. Comparison of fCO2sw data acquired within two 1°×1° regions in the northern and southern North Sea during different years (1987, 2001, 2002, and 2005–2007) revealed large interannual variations, especially during spring and summer when year-to-year fCO2sw differences (≈160–200 μatm) approached seasonal changes (≈200–250 μatm). The springtime variations resulted from changes in magnitude and timing of the phytoplankton bloom, whereas changes in SST, wind speed and total alkalinity may have contributed to the summertime interannual fCO2sw differences. The lowest interannual variation (10–50 μatm) was observed during fall and early winter. Comparison with data reported in October 1967 suggests that the fCO2sw growth rate in the central North Sea was similar to that in the atmosphere.

Comparison of global ocean colour data records - S. Djavidnia; F. Mélin; N. Hoepffner (Ocean Science, 2010)
The extending record of ocean colour derived information, an important asset for the study of marine ecosystems and biogeochemistry, presently relies on individual satellite missions launched by several space agencies with differences in sensor design, calibration strategies and algorithms. In this study we present an extensive comparative analysis of standard products obtained from operational global ocean colour sensors (SeaWiFS, MERIS, MODIS-Aqua, MODIS-Terra), on both global and regional scales. The analysis is based on monthly mean chlorophyll a (Chl-a) sea surface concentration between 2002 and 2009. Based on global statistics, the Chl-a records appear relatively consistent. The root mean square (RMS) difference Δ between (log-transformed) Chl-a from SeaWiFS and MODIS Aqua amounts to 0.137, with a bias of 0.074 (SeaWiFS Chl-a higher). The difference between these two products and MERIS Chl-a is approximately 0.15. Restricting the analysis to 2007 only, Δ between MODIS Aqua and Terra is 0.142. This global convergence is significantly modulated regionally. Statistics for biogeographic provinces representing a partition of the global ocean, show Δ values varying between 0.08 and 0.3. High latitude regions, as well as coastal and shelf provinces are generally the areas with the largest differences. Moreover, RMS differences and biases are modulated in time, with a coefficient of variation of Δ varying between 10% and 40%, with clear seasonal patterns in some provinces. The comparison of the province-averaged time series obtained from the various satellite products also shows a level of agreement that is geographically variable. Overall, the Chl-a SeaWiFS and MODIS Aqua series appear to have similar levels of variance and display high correlation coefficients, an agreement likely favoured by the common elements shared by the two missions. These results are degraded if the MERIS series is compared to either SeaWiFS or MODIS Aqua. An important outcome of the study is that the results of the inter-comparison analysis are variable with time and location, and therefore globally averaged statistics are not necessarily applicable on a seasonal or regional basis.

What are "ecogenomic sensors?" A review and thoughts for the future - C. A. Scholin (Ocean Science, 2010)
The OceanSensors08 workshop held 31 March–4 April 2008 in Warnemünde, Germany, brought together an international group of marine scientists, sensor developers and technologists with a common interest in shaping the future of ocean sensing. In preparation for that meeting a series of review papers was commissioned, one of which was meant to cover Sensors for Ocean-omics. The "ocean-omics" topic was cast very broadly. The notion was to review use of genetic techniques for assessing presence and diversity of organisms, their genomic capacity and gene expression, and to provide a prospectus of how such methods could be applied in an autonomous capacity. I chose "ecogenomic sensor" as a descriptor to convey the essence of such a system – a device that integrates genetic level sensing with larger scale environmental characterization. This phrase is derived from workshops refining the US's Ocean Observatories Initiative (OOI) and visions for instrument systems that could be deployed on such a network. But what exactly are ecogenomic sensors? A clear definition is lacking and conceptualizations far outweigh actual hardware that can be deployed in the ocean. This prospectus builds from that point. I advance a definition of "ecogenomic sensor" and outline the opportunities and challenges associated with developing such instruments. Suggestions as to how this technology may be further refined and applied are offered against the backdrop of the Autonomous Microbial Genosensor (AMG) and Environmental Sample Processor (ESP). Applications that center on detection of DNA and RNA are emphasized. The word "review" appears in the title at the request of the editors.

Indian Ocean Subtropical Mode Water: its water characteristics and spatial distribution - T. Tsubouchi; T. Suga; K. Hanawa (Ocean Science, 2010)
We have improved a basic description (water characteristics and spatial distribution) of the Indian Ocean Subtropical Mode Water (IOSTMW) using an isopycnally averaged three-dimensional hydrographic dataset. Two mode waters and corresponding wintertime mixed layer depth maxima were observed north of the subtropical front (STF) in the South Indian Ocean: IOSTMW (within 25.8–26.2 σθ) in the region of 28–45° E and another subtropical mode water in the subtropical gyre (within 26.4–26.7 σθ) in the 60–80° E longitudinal band. Through careful examination of the spatial distribution and water characteristics of a core in the layer of minimum vertical temperature gradient (LMVTG), we identified that a mass of LMVTG corresponds to IOSTMW. The average water characteristics of the IOSTMW during approximately 1960–2004 were 16.54 ± 0.49 °C, 35.51 ± 0.04 psu and 26.0 ± 0.1 σθ. The IOSTMW distribution area was estimated to be 25–50° E, 27–38° S. The formation region and approximate water characteristics of the second subtropical mode water were also estimated. Its probable formation region was 37–42° S, 60–80° E and north of the STF, with approximate water characteristics of 12.84 ± 0.57 °C, 35.17 ± 0.11 psu and 26.57 ± 0.04 σθ.

Pre-operational short-term forecasts for Mediterranean Sea biogeochemistry - P. Lazzari; A. Teruzzi; S. Salon; S. Campagna; C. Calonaci; S. Colella; M. Tonani; A. Crise (Ocean Science, 2010)
Operational prediction of the marine environment is recognised as a fundamental research issue in Europe. We present a pre-operational implementation of a biogeochemical model for the pelagic waters of the Mediterranean Sea, developed within the framework of the MERSEA-IP European project. The OPATM-BFM coupled model is the core of a fully automatic system that delivers weekly analyses and forecast maps for the Mediterranean Sea biogeochemistry. The system has been working in its current configuration since April 2007 with successful execution of the fully automatic operational chain in 87% of the cases while in the remaining cases the runs were successfully accomplished after operator intervention. A description of the system developed and also a comparison of the model results with satellite data are presented, together with a measure of the model skill evaluated by means of seasonal target diagrams. Future studies will address the implementation of a data assimilation scheme for the biogeochemical compartment in order to increase the skill of the model's performance.

Density and Absolute Salinity of the Baltic Sea 2006–2009 - R. Feistel; S. Weinreben; H. Wolf; S. Seitz; P. Spitzer; B. Adel; G. Nausch; B. Schneider; D. G. Wright (Ocean Science, 2010)
The brackish water of the Baltic Sea is a mixture of ocean water from the Atlantic/North Sea with fresh water from various rivers draining a large area of lowlands and mountain ranges. The evaporation-precipitation balance results in an additional but minor excess of fresh water. The rivers carry different loads of salts washed out of the ground, in particular calcium carbonate, which cause a composition anomaly of the salt dissolved in the Baltic Sea in comparison to Standard Seawater. Directly measured seawater density shows a related anomaly when compared to the density computed from the equation of state as a function of Practical Salinity, temperature and pressure. Samples collected from different regions of the Baltic Sea during 2006–2009 were analysed for their density anomaly. The results obtained for the river load deviate significantly from similar measurements carried out forty years ago; the reasons for this decadal variability are not yet fully understood. An empirical formula is derived which estimates Absolute from Practical Salinity of Baltic Sea water, to be used in conjunction with the new Thermodynamic Equation of Seawater 2010 (TEOS-10), endorsed by IOC/UNESCO in June 2009 as the substitute for the 1980 International Equation of State, EOS-80. Our routine measurements of the samples were accompanied by studies of additional selected properties which are reported here: conductivity, density, chloride, bromide and sulphate content, total CO2 and alkalinity.

Corrigendum to "Observed and simulated estimates of the meridional overturning circulation at 26.5° N in the Atlantic" published in Ocean Sci., 5, 575–589, 2009 - J. Baehr; S. Cunningham; H. Haak; P. Heimbach; T. Kanzow; J. Marotzke (Ocean Science, 2010)

Carbon export and sequestration in the southern Benguela upwelling system: lower and upper estimates - H. N. Waldron; P. M. S. Monteiro; N. C. Swart (Ocean Science, 2009)
Three independent studies of carbon export and sequestration in the southern Benguela upwelling system are presented. They were undertaken by Waldron (upwelling index), Monteiro (discrete upwelling centres – gate hypothesis model) and Swart (cross-shelf advection in bottom nelpheloid layers). The annual estimates were, 3.9×1013, 0.72×1013 and 8.6×1011 gC, respectively. The lowest estimate was derived from a consideration of low frequency lateral carbon export in the bottom nepheloid layer and was thought likely to be an under-estimate. Taking into account high frequency episodic events, intermediate nepheloid layers and along isopycnal export of DOC at surface and intermediate depths was thought likely to result in a substantial upward revision. The remaining two estimates were considered to be an upper and lower estimate of carbon export and sequestration due to factors inherent in the methodologies. The upper estimate presents a two-dimensional system, integrated alongshore; the lower estimate sums a series of upwelling centres in order to obtain a system flux. The former is therefore a uniform extrapolation along the coast while the latter omits upwelling between the upwelling centres.

Electrochemical techniques and sensors for ocean research - G. Denuault (Ocean Science, 2009)
This paper presents a review of applications of electrochemical methods in ocean sensing. It follows the white paper presented at the OceanSensors08 workshop held at the Leibniz-Institut für Ostseeforschung, Warnemünde, Germany, from 31 March to 4 April 2008. The principles of electrochemical techniques are briefly recalled and described. For each technique, relevant electrochemical sensors are discussed; known successful deployments of electrochemical sensors are recalled; challenges experienced when taking sensors from the research lab to the field are raised; future trends in development and applications are proposed and assessed for their potential for oceanographic applications; where possible technological readiness levels are estimated. The document is supported with references drawn from both the electrochemical and oceanographic literature.

Technical Note: Animal-borne CTD-Satellite Relay Data Loggers for real-time oceanographic data collection - L. Boehme; P. Lovell; M. Biuw; F. Roquet; J. Nicholson; S. E. Thorpe; M. P. Meredith; M. Fedak (Ocean Science, 2009)
The increasing need for continuous monitoring of the world oceans has stimulated the development of a range of autonomous sampling platforms. One novel addition to these approaches is a small, relatively inexpensive data-relaying device that can be deployed on marine mammals to provide vertical oceanographic profiles throughout the upper 2000 m of the water column. When an animal dives, the CTD-Satellite Relay Data Logger (CTD-SRDL) records vertical profiles of temperature, conductivity and pressure. Data are compressed once the animal returns to the surface where it is located by, and relays data to, the Argos satellite system. The technical challenges met in the design of the CTD-SRDL are the maximising of energy efficiency and minimising size, whilst simultaneously maintaining the reliability of an instrument that cannot be recovered and is required to survive its lifetime attached to a marine mammal. The CTD-SRDLs record temperature and salinity with an accuracy of better than 0.005 °C and 0.02 respectively. However, due to the limited availability of reference data, real-time data from remote places are often associated with slightly higher errors. The potential to collect large numbers of profiles cost-effectively makes data collection using CTD-SRDL technology particularly beneficial in regions where traditional oceanographic measurements are scarce or even absent. Depending on the CTD-SRDL configuration, it is possible to sample and transmit hydrographic profiles on a daily basis, providing valuable and often unique information for a real-time ocean observing system.

Optical tools for ocean monitoring and research - C. Moore; A. Barnard; P. Fietzek; M. R. Lewis; H. M. Sosik; S. White; O. Zielinski (Ocean Science, 2009)
Requirements for understanding the relationships between ocean color and suspended and dissolved materials within the water column, and a rapidly emerging photonics and materials technology base for performing optical based analytical techniques have generated a diverse offering of commercial sensors and research prototypes that perform optical measurements in water. Through inversion, these tools are now being used to determine a diverse set of related biogeochemical and physical parameters. Techniques engaged include measurement of the solar radiance distribution, absorption, scattering, stimulated fluorescence, flow cytometry, and various spectroscopy methods. Selective membranes and other techniques for material isolation further enhance specificity, leading to sensors for measurement of dissolved oxygen, methane, carbon dioxide, common nutrients and a variety of other parameters. Scientists are using these measurements to infer information related to an increasing set of parameters and wide range of applications over relevant scales in space and time.

Deep ocean exchange with west-European shelf seas - J. M. Huthnance; J. T. Holt; S. L. Wakelin (Ocean Science, 2009)
We review mechanisms and studies of exchange between the north-east Atlantic and the adjacent shelf seas. Well-developed summer upwelling and associated filaments off Portugal and north-west Spain give exchange O(3 m2/s per unit length of shelf). Prevailing westerly winds further north drive exchange O(1 m2/s). Poleward flow along most of the upper slope has associated secondary circulation O(1 m2/s), meanders and eddies. Eddies are shed from slope waters into the Bay of Biscay, and local exchanges occur at shelf spurs and depressions or canyons (e.g. dense-water cascading of order 1 m2/s). Tidal transports are larger, but their reversal every six hours makes exchange largely ineffective except where internal tides are large and non-linear, as in the Celtic Sea where solitons carry water with exchange O(1 m2/s). These various physical exchanges amount to an estimated 2–3 m2/s per unit length of shelf, between ocean and shelf. A numerical model estimate is comparable: 2.5×106 m3/s onto and off the shelf from Brittany to Norway. Mixing controls the seasonal thermocline, affecting primary production and hence fluxes and fate of organic matter. Specifically, CO2 take-up by primary production, settling below the thermocline before respiration, and then off-shelf transport, make an effective shelf-sea "pump" (for CO2 from the atmosphere to the deep ocean). However, knowledge of biogeochemical fluxes is generally sparse, giving scope for more measurements, model validation and estimates from models.

A review of the role of submarine canyons in deep-ocean exchange with the shelf - S. E. Allen; X. Durrieu de Madron (Ocean Science, 2009)
Cross shelf-break exchange is limited by the tendency of geostrophic flow to follow bathymetric contours, not cross them. However, small scale topography, such as canyons, can reduce the local lengthscale of the flow and increase the local Rossby number. These higher Rossby numbers mean the flow is no longer purely geostrophic and significant cross-isobath flow can occur. This cross-isobath flow includes both upwelling and downwelling due to wind-driven shelf currents and the strong cascading flows of dense shelf-water into the ocean. Tidal currents usually run primarily parallel to the shelf-break topography. Canyons cut across these flows and thus are often regions of generation of strong baroclinic tides and internal waves. Canyons can also focus internal waves. Both processes lead to greatly elevated levels of mixing. Thus, through both advection and mixing processes, canyons can enhance Deep Ocean Shelf Exchange. Here we review the state of the science describing the dynamics of the flows and suggest further areas of research, particularly into quantifying fluxes of nutrients and carbon as well as heat and salt through canyons.

Impact of data assimilation of physical variables on the spring bloom from TOPAZ operational runs in the North Atlantic - A. Samuelsen; L. Bertino; C. Hansen (Ocean Science, 2009)
A reanalysis of the North Atlantic spring bloom in 2007 was produced using the real-time analysis from the TOPAZ North Atlantic and Arctic forecasting system. The TOPAZ system uses a hybrid coordinate general circulation ocean model and assimilates physical observations: sea surface anomalies, sea surface temperatures, and sea-ice concentrations using the Ensemble Kalman Filter. This ocean model was coupled to an ecosystem model, NORWECOM (Norwegian Ecological Model System), and the TOPAZ-NORWECOM coupled model was run throughout the spring and summer of 2007. The ecosystem model was run online, restarting from analyzed physical fields (result after data assimilation) every 7 days. Biological variables were not assimilated in the model. The main purpose of the study was to investigate the impact of physical data assimilation on the ecosystem model. This was determined by comparing the results to those from a model without assimilation of physical data. The regions of focus are the North Atlantic and the Arctic Ocean. Assimilation of physical variables does not affect the results from the ecosystem model significantly. The differences between the weekly mean values of chlorophyll are normally within 5–10% during the summer months, and the maximum difference of ~20% occurs in the Arctic, also during summer. Special attention was paid to the nutrient input from the North Atlantic to the Nordic Seas and the impact of ice-assimilation on the ecosystem. The ice-assimilation increased the phytoplankton concentration: because there was less ice in the assimilation run, this increased both the mixing of nutrients during winter and the area where production could occur during summer. The forecast was also compared to remotely sensed chlorophyll, climatological nutrients, and in-situ data. The results show that the model reproduces a realistic annual cycle, but the chlorophyll concentrations tend to be between 0.1 and 1.0 mg chla/m3 too low during winter and spring and 1–2 mg chla/m3 too high during summer. Surface nutrients on the other hand are generally lower than the climatology throughout the year.

Mediterranean Forecasting System: forecast and analysis assessment through skill scores - M. Tonani; N. Pinardi; C. Fratianni; J. Pistoia; S. Dobricic; S. Pensieri; M. de Alfonso; K. Nittis (Ocean Science, 2009)
This paper describes the first evaluation of the quality of the forecast and analyses produced at the basin scale by the Mediterranean ocean Forecasting System (MFS) (http://gnoo.bo.ingv.it/mfs). The system produces short-term ocean forecasts for the following ten days. Analyses are produced weekly using a daily assimilation cycle. The analyses are compared with independent data from buoys, where available, and with the assimilated data before the data are inserted. In this work we have considered 53 ten days forecasts produced from 16 August 2005 to 15 August 2006. The forecast skill is evaluated by means of root mean square error (rmse) differences, bias and anomaly correlations at different depths for temperature and salinity, computing differences between forecast and analysis, analysis and persistence and forecast and persistence. The Skill Score (SS) is defined as the ratio of the rmse of the difference between analysis and forecast and the rmse of the difference between analysis and persistence. The SS shows that at 5 and 30 m the forecast is always better than the persistence, but at 300 m it can be worse than persistence for the first days of the forecast. This result may be related to flow adjustments introduced by the data assimilation scheme. The monthly variability of SS shows that when the system variability is high, the values of SS are higher, therefore the forecast has higher skill than persistence. We give evidence that the error growth in the surface layers is controlled by the atmospheric forcing inaccuracies, while at depth the forecast error can be interpreted as due to the data insertion procedure. The data, both in situ and satellite, are not homogeneously distributed in the basin; therefore, the quality of the analyses may be different in different areas of the basin.

On the multiple time scales of variability in the Northeast Pacific Ocean - R. Tokmakian (Ocean Science, 2009)
The spatial and temporal sea surface height energy distribution of the Northeast Pacific Ocean is described and discussed. Using an altimetric data set covering 15 years (1993–2007), the energy within the 3–9 month band is primarily located within 10° of the coast. In the Gulf of Alaska, this energy signal is on the shelf, while further south, west of the California/Oregon coast, the significant energy in this band is west of the shelf break. In both cases, it is primarily forced by the local wind. Within the 2–3 year band, the signal reflects energy generated by local changes to the wind stress from large atmospheric shifts indicated by the Pacific North American Index and by advective or propagating processes related to El Niño-Southern Oscillation. Over the two 4–6 year periods within this data set, the change is primarily due to the large scale shift in atmospheric systems north of about 30° N which also affect changes in current strengths. Based on the distribution of the energy signal and its variability, a set of three winter-time indices are suggested to characterize the distinct differences in the SSH anomalies in these areas.

Observed and simulated estimates of the meridional overturning circulation at 26.5° N in the Atlantic - J. Baehr; S. Cunnningham; H. Haak; P. Heimbach; T. Kanzow; J. Marotzke (Ocean Science, 2009)
Daily timeseries of the meridional overturning circulation (MOC) estimated from the UK/US RAPID/MOCHA array at 26.5° N in the Atlantic are used to evaluate the MOC as simulated in two global circulation models: (I) an 8-member ensemble of the coupled climate model ECHAM5/MPI-OM, and (II) the ECCO-GODAE state estimate. In ECHAM5/MPI-OM, we find that the observed and simulated MOC have a similar variability and time-mean within the 99% confidence interval. In ECCO-GODAE, we find that the observed and simulated MOC show a significant correlation within the 99% confidence interval. To investigate the contribution of the different transport components, the MOC is decomposed into Florida Current, Ekman and mid-ocean transports. In both models, the mid-ocean transport is closely approximated by the residual of the MOC minus Florida Current and Ekman transports. As the models conserve volume by definition, future comparisons of the RAPID/MOCHA mid-ocean transport should be done against the residual transport in the models. The similarity in the variance and the correlation between the RAPID/MOCHA, and respectively ECHAM5/MPI-OM and ECCO-GODAE MOC estimates at 26.5° N is encouraging in the context of estimating (natural) variability in climate simulations and its use in climate change signal-to-noise detection analyses. Enhanced confidence in simulated hydrographic and transport variability will require longer observational time series.

Sensors and instruments for oceanic dissolved carbon measurements - U. Schuster; A. Hannides; L. Mintrop; A. Körtzinger (Ocean Science, 2009)
Highly accurate and precise measurements of marine carbon components are required in the study of the marine carbon cycle, particularly when investigating the causes for its variability from seasonal to interannual timescales. This is especially true in the investigation of the consequences of anthropogenic influences. The analysis of any marine carbon component requires elaborate instrumentation, most of which is currently used onboard ships, either in manual or automated mode. Technological developments result in more and more instruments that have sufficient long-term reliability so that they can be deployed on commercial ships, surface moorings, and buoys, whilst the great technological and operational challenges mean that only few sensors have been developed that can be used for sub-surface in situ measurements on floats, robots, or gliders. There is a special need for autonomous instruments and sensors that are able to measure a combination of different components, in order to increase the spatial and temporal coverage of marine carbon data. This paper describes analytical techniques used for the measurement of the marine dissolved carbon components, both inorganic and organic: the fugacity of CO2, total dissolved inorganic carbon, pH, alkalinity, and dissolved organic carbon. By pointing out advantages, disadvantages, and/or challenges of the techniques employed in the analysis of each component, we aim to aid non-carbon marine scientists, sensor developers and technologists, in the decision of which challenges to address in further development.

Surface circulation in the Eastern Mediterranean using drifters (2005–2007) - R. Gerin; P.-M. Poulain; I. Taupier-Letage; C. Millot; S. Ben Ismail; C. Sammari (Ocean Science, 2009)
Within the framework of the EGITTO/EGYPT programs, the spatial structure and the temporal variability of the surface circulation in the eastern basin of the Mediterranean Sea was studied with satellite-tracked drifters. A total of 97 drifters drogued to 15-m nominal depth were released between September 2005 and March 2007, regularly along ship-of-opportunity routes in the Sicily Channel and within specific structures during dedicated campaigns in the Levantine sub-basin. This study complements previous ones in the former and in the Ionian, but it is the first one in the latter. After editing and low-pass filtering, the drifter trajectories were used to estimate pseudo-Eulerian statistics: mean current, variance ellipses, mean and eddy kinetic energies. A statistical analysis was also performed dividing the dataset in two extended seasons (winter and summer). A branching behaviour of the surface water after passing through the Sicily Channel is evidenced, together with a seasonal variability inducing a reversal of the circulation in the southern part of the Ionian. In the Levantine, the surface circulation describes an eastward flow along the Libyan and Egyptian slopes (the Libyo-Egyptian Current: LEC) that continues in a cyclonic circuit along the Middle East and Turkish slopes. This general alongslope circuit can be perturbed locally and temporally by the numerous anticyclonic eddies that co-exist in the Levantine, mainly created by the instability of the LEC in the south (Libyo-Egyptian Eddies: LEEs), but also by the wind (Ierapetra and Pelops), and by the topography (over the Eratosthenes Seamount and off Latakia). The most frequent perturbation is the entrainment of part of the flow seaward: the LEEs close to the slope can interact with the LEC, which then spreads more or less around them, so that a series of contiguous LEEs (paddle-wheel effect) can possibly result in an eastward offshore transport (the so-called Mid Mediterranean Jet). Additionally, when LEEs are close enough to the slope, most of the surface flow is spread seaward. Along the slope the current is then induced by the LEEs southern side which results in a westward current. Locally and temporally the circulation along the slope can thus be reversed.

A chemical ionization mass spectrometer for continuous underway shipboard analysis of dimethylsulfide in near-surface seawater - E. S. Saltzman; W. J. De Bruyn; M. J. Lawler; C. A. Marandino; C. A. McCormick (Ocean Science, 2009)
A compact, low-cost atmospheric pressure, chemical ionization mass spectrometer ("mini-CIMS") has been developed for continuous underway shipboard measurements of dimethylsulfide (DMS) in seawater. The instrument was used to analyze DMS in air equilibrated with flowing seawater across a porous Teflon membrane equilibrator. The equilibrated gas stream was diluted with air containing an isotopically-labeled internal standard. DMS is ionized at atmospheric pressure via proton transfer from water vapor, then declustered, mass filtered via quadrupole mass spectrometry, and detected with an electron multiplier. The instrument described here is based on a low-cost residual gas analyzer (Stanford Research Systems), which has been modified for use as a chemical ionization mass spectrometer. The mini-CIMS has a gas phase detection limit of 220 ppt DMS for a 1 min averaging time, which is roughly equivalent to a seawater DMS concentration of 0.1 nM DMS at 20°C. The mini-CIMS has the sensitivity, selectivity, and time response required for underway measurements of surface ocean DMS over the full range of oceanographic conditions. The simple, robust design and relatively low cost of the instrument are intended to facilitate use in process studies and surveys, with potential for long-term deployment on research vessels, ships of opportunity, and large buoys.

Sensors for observing ecosystem status - S. Kröger; E. R. Parker; J. D. Metcalfe; N. Greenwood; R. M. Forster; D. B. Sivyer; D. J. Pearce (Ocean Science, 2009)
This paper aims to review the availability and application of sensors for observing marine ecosystem status. It gives a broad overview of important ecosystem variables to be investigated, such as biogeochemical cycles, primary and secondary production, species distribution, animal movements, habitats and pollutants. Some relevant legislative drivers are listed, as they provide one context in which ecosystem studies are undertaken. In addition to literature cited within the text the paper contains some useful web links to assist the reader in making an informed instrument choice, as the authors feel that the topic is so broad, it is impossible to discuss all relevant systems or to provide appropriate detail for those discussed. It is therefore an introduction to how and why ecosystem status is currently observed, what variables are quantified, from what platforms, using remote sensing or in-situ measurements, and gives examples of useful sensor based tools. Starting with those presently available, to those under development and also highlighting sensors not yet realised but desirable for future studies.

Relating Agulhas leakage to the Agulhas Current retroflection location - E. van Sebille; C. N. Barron; A. Biastoch; P. J. van Leeuwen; F. C. Vossepoel; W. P. M. de Ruijter (Ocean Science, 2009)
The relation between the Agulhas Current retroflection location and the magnitude of Agulhas leakage, the transport of water from the Indian to the Atlantic Ocean, is investigated in a high-resolution numerical ocean model. Sudden eastward retreats of the Agulhas Current retroflection loop are linearly related to the shedding of Agulhas rings, where larger retreats generate larger rings. Using numerical Lagrangian floats a 37 year time series of the magnitude of Agulhas leakage in the model is constructed. The time series exhibits large amounts of variability, both on weekly and annual time scales. A linear relation is found between the magnitude of Agulhas leakage and the location of the Agulhas Current retroflection, both binned to three month averages. In the relation, a more westward location of the Agulhas Current retroflection corresponds to an increased transport from the Indian Ocean to the Atlantic Ocean. When this relation is used in a linear regression and applied to almost 20 years of altimetry data, it yields a best estimate of the mean magnitude of Agulhas leakage of 13.2 Sv. The early retroflection of 2000, when Agulhas leakage was probably halved, can be identified using the regression.

Problems with estimation and interpretation of oceanic heat transport – conceptual remarks for the case of Fram Strait in the Arctic Ocean - U. Schauer; A. Beszczynska-Möller (Ocean Science, 2009)
While the concept of oceanic heat transport – or rather heat transport divergence – is well known, it is sometimes applied inaccurately. Often so-called "heat transports" are computed across a partial section which means that the volume flow through such a section is not zero. In this case the "heat transports" depend entirely on the choice of the temperature scale. The consequences of such arbitrariness are demonstrated with a simple calculation exercise for the passages to the Arctic Ocean. To circumvent the arising difficulties for the Fram Strait in the Arctic we propose a stream tube concept to define a net zero volume flow section which can, with coarse assumptions, be used to determine oceanic heat transport by the portion of Atlantic water flow that passes through Fram Strait. Weaknesses of this approach and consequences for observational strategies are discussed.

Application of the Gaussian anamorphosis to assimilation in a 3-D coupled physical-ecosystem model of the North Atlantic with the EnKF: a twin experiment - E. Simon; L. Bertino (Ocean Science, 2009)
We consider the application of the Ensemble Kalman Filter (EnKF) to a coupled ocean ecosystem model (HYCOM-NORWECOM). Such models, especially the ecosystem models, are characterized by strongly non-linear interactions active in ocean blooms and present important difficulties for the use of data assimilation methods based on linear statistical analysis. Besides the non-linearity of the model, one is confronted with the model constraints, the analysis state having to be consistent with the model, especially with respect to the constraints that some of the variables have to be positive. Furthermore the non-Gaussian distributions of the biogeochemical variables break an important assumption of the linear analysis, leading to a loss of optimality of the filter. We present an extension of the EnKF dealing with these difficulties by introducing a non-linear change of variables (anamorphosis function) in order to execute the analysis step in a Gaussian space, namely a space where the distributions of the transformed variables are Gaussian. We present also the initial results of the application of this non-Gaussian extension of the EnKF to the assimilation of simulated chlorophyll surface concentration data in a North Atlantic configuration of the HYCOM-NORWECOM coupled model.

Enhancing temporal correlations in EOF expansions for the reconstruction of missing data using DINEOF - A. Alvera-Azcárate; A. Barth; D. Sirjacobs; J.-M. Beckers (Ocean Science, 2009)
DINEOF (Data Interpolating Empirical Orthogonal Functions) is an EOF-based technique for the reconstruction of missing data in geophysical fields, such as those produced by clouds in sea surface temperature satellite images. A technique to reduce spurious time variability in DINEOF reconstructions is presented. The reconstruction of these images within a long time series using DINEOF can lead to large discontinuities in the reconstruction. Filtering the temporal covariance matrix allows to reduce this spurious variability and therefore more realistic reconstructions are obtained. The approach is tested in a three years sea surface temperature data set over the Black Sea. The effect of the filter in the temporal EOFs is presented, as well as some examples of the improvement achieved with the filtering in the SST reconstruction, both compared to the DINEOF approach without filtering.

The multifractal structure of satellite sea surface temperature maps can be used to obtain global maps of streamlines - A. Turiel; V. Nieves; E. Garcia-Ladona; J. Font; M.-H. Rio; G. Larnicol (Ocean Science, 2009)
Nowadays Earth observation satellites provide information about many relevant variables of the ocean-climate system, such as temperature, moisture, aerosols, etc. However, to retrieve the velocity field, which is the most relevant dynamical variable, is still a technological challenge, specially in the case of oceans. New processing techniques, emerged from the theory of turbulent flows, have come to assist us in this task. In this paper, we show that multifractal techniques applied to new Sea Surface Temperature satellite products opens the way to build maps of ocean currents with unprecedented accuracy. With the application of singularity analysis, we show that global ocean circulation patterns can be retrieved in a daily basis. We compare these results with high-quality altimetry-derived geostrophic velocities, finding a quite good correspondence of the observed patterns both qualitatively and quantitatively; and this is done for the first time on a global basis, even for less active areas. The implications of this findings from the perspective both of theory and of operational applications are discussed.

Observations of turbulence beneath sea ice in southern McMurdo Sound, Antarctica - C. L. Stevens; N. J. Robinson; M. J. M. Williams; T. G. Haskell (Ocean Science, 2009)
The first turbulence profiler observations beneath land fast sea ice which is directly adjacent to an Antarctic ice shelf are described. The stratification in the 325 m deep water column consisted of a layer of supercooled water in the upper 40 m lying above a quasi-linearly stratified water column with a sharp step in density at mid-depth. Turbulent energy dissipation rates were on average 3×10−8 m2 s−3 with peak bin-averaged values reaching 4×10−7 m2 s−3. The local dissipation rate per unit area was estimated to be 10 m Wm−2 on average with a peak of 50 m Wm−2. These values are consistent with a moderate baroclinic response to the tides. The small-scale turbulent energetics lie on the boundary between isotropy and buoyancy-affected. This will likely influence the formation and aggregation of frazil ice crystals within the supercooled layer. The data suggest that the large crystals observed in McMurdo Sound will transition from initial growth at scales smaller than the Kolmogorov lengthscale to sizes substantially (1–2 orders of magnitude) greater than the Kolmogorov scale. An estimate of the experiment-averaged vertical diffusivity of mass Kρ yields a coefficient of around 2×10−4 m2s−1 although this increased by a factor of 2 near the surface. Combining this estimate of Kρ with available observations of average and maximum currents suggests the layer of supercooled water can persist for a distance of ~250 km from the front of the McMurdo Ice Shelf.

A nested Atlantic-Mediterranean Sea general circulation model for operational forecasting - P. Oddo; M. Adani; N. Pinardi; C. Fratianni; M. Tonani; D. Pettenuzzo (Ocean Science, 2009)
A new numerical general circulation ocean model for the Mediterranean Sea has been implemented nested within an Atlantic general circulation model within the framework of the Marine Environment and Security for the European Area project (MERSEA, Desaubies, 2006). A 4-year twin experiment was carried out from January 2004 to December 2007 with two different models to evaluate the impact on the Mediterranean Sea circulation of open lateral boundary conditions in the Atlantic Ocean. One model considers a closed lateral boundary in a large Atlantic box and the other is nested in the same box in a global ocean circulation model. Impact was observed comparing the two simulations with independent observations: ARGO for temperature and salinity profiles and tide gauges and along-track satellite observations for the sea surface height. The improvement in the nested Atlantic-Mediterranean model with respect to the closed one is particularly evident in the salinity characteristics of the Modified Atlantic Water and in the Mediterranean sea level seasonal variability.

Adjustment of the basin-scale circulation at 26° N to variations in Gulf Stream, deep western boundary current and Ekman transports as observed by the Rapid array - H. L. Bryden; A. Mujahid; S. A. Cunningham; T. Kanzow (Ocean Science, 2009)
The Rapid instrument array across the Atlantic Ocean along 26° N provides unprecedented monitoring of the basin-scale circulation. A unique feature of the Rapid array is the combination of full-depth moorings with instruments measuring temperature, salinity, pressure time series at many depths with co-located bottom pressure measurements so that dynamic pressure can be measured from surface to bottom. Bottom pressure measurements show a zonally uniform rise (and fall) of bottom pressure of 0.015 dbar on a 5 to 10 day time scale, suggesting that the Atlantic basin is filling and draining on a short time scale. After removing the zonally uniform bottom pressure fluctuations, bottom pressure variations at 4000 m depth against the western boundary compensate instantaneously for baroclinic fluctuations in the strength and structure of the deep western boundary current so there is no basin-scale mass imbalance resulting from variations in the deep western boundary current. After removing the mass compensating bottom pressure, residual bottom pressure fluctuations at the western boundary just east of the Bahamas balance variations in Gulf Stream transport. Again the compensation appears to be especially confined close to the western boundary. Thus, fluctuations in either Gulf Stream or deep western boundary current transports are compensated in a depth independent (barotropic) manner very close to the continental slope off the Bahamas. In contrast, compensation for variations in wind-driven surface Ekman transport appears to involve fluctuations in both western basin and eastern basin bottom pressures, though the bottom pressure difference fluctuations appear to be a factor of 3 too large, perhaps due to an inability to resolve small bottom pressure fluctuations after removal of larger zonal average, baroclinic, and Gulf Stream pressure components. For 4 tall moorings where time series dynamic height (geostrophic pressure) profiles can be estimated from sea surface to ocean bottom and bottom pressure can be added, there is no general correlation between surface dynamic height and bottom pressure. Dynamic height on each mooring is strongly correlated with sea surface height from satellite observations and the variability in both dynamic height and satellite sea surface height decrease sharply as the western boundary is approached.

Controlling atmospheric forcing parameters of global ocean models: sequential assimilation of sea surface Mercator-Ocean reanalysis data - C. Skandrani; J.-M. Brankart; N. Ferry; J. Verron; P. Brasseur; B. Barnier (Ocean Science, 2009)
In the context of stand alone ocean models, the atmospheric forcing is generally computed using atmospheric parameters that are derived from atmospheric reanalysis data and/or satellite products. With such a forcing, the sea surface temperature that is simulated by the ocean model is usually significantly less accurate than the synoptic maps that can be obtained from the satellite observations. This not only penalizes the realism of the ocean long-term simulations, but also the accuracy of the reanalyses or the usefulness of the short-term operational forecasts (which are key GODAE and MERSEA objectives). In order to improve the situation, partly resulting from inaccuracies in the atmospheric forcing parameters, the purpose of this paper is to investigate a way of further adjusting the state of the atmosphere (within appropriate error bars), so that an explicit ocean model can produce a sea surface temperature that better fits the available observations. This is done by performing idealized assimilation experiments in which Mercator-Ocean reanalysis data are considered as a reference simulation describing the true state of the ocean. Synthetic observation datasets for sea surface temperature and salinity are extracted from the reanalysis to be assimilated in a low resolution global ocean model. The results of these experiments show that it is possible to compute piecewise constant parameter corrections, with predefined amplitude limitations, so that long-term free model simulations become much closer to the reanalysis data, with misfit variance typically divided by a factor 3. These results are obtained by applying a Monte Carlo method to simulate the joint parameter/state prior probability distribution. A truncated Gaussian assumption is used to avoid the most extreme and non-physical parameter corrections. The general lesson of our experiments is indeed that a careful specification of the prior information on the parameters and on their associated uncertainties is a key element in the computation of realistic parameter estimates, especially if the system is affected by other potential sources of model errors.

Application of a hybrid EnKF-OI to ocean forecasting - F. Counillon; P. Sakov; L. Bertino (Ocean Science, 2009)
Data assimilation methods often use an ensemble to represent the background error covariance. Two approaches are commonly used; a simple one with a static ensemble, or a more advanced one with a dynamic ensemble. The latter is often non-practical due to its high computational requirements. Some recent studies suggested using a hybrid covariance, which is a linear combination of the covariances represented by a static and a dynamic ensemble. Here, the use of the hybrid covariance is first extensively tested with a quasi-geostrophic model and with different analysis schemes, namely the Ensemble Kalman Filter (EnKF) and the Ensemble Square Root Filter (ESRF). The hybrid covariance ESRF (ESRF-OI) is more accurate and more stable than the hybrid covariance EnKF (EnKF-OI), but the overall conclusions are similar regardless of the analysis scheme used. The benefits of using the hybrid covariance are large compared to both the static and the dynamic methods with a small dynamic ensemble. The benefits over the dynamic methods become negligible, but remain, for large dynamic ensembles. The optimal value of the hybrid blending coefficient appears to decrease exponentially with the size of the dynamic ensemble. Finally, we consider a realistic application with the assimilation of altimetry data in a hybrid coordinate ocean model (HYCOM) for the Gulf of Mexico, during the shedding of Eddy Yankee (2006). A 10-member EnKF-OI is compared to a 10-member EnKF and a static method called the Ensemble Optimal Interpolation (EnOI). While 10 members seem insufficient for running the EnKF, the 10-member EnKF-OI reduces the forecast error compared to the EnOI, and improves the positions of the fronts.

On the porosity of barrier layers - J. Mignot; C. de Boyer Montégut; M. Tomczak (Ocean Science, 2009)
Barrier layers are defined as the layer between the pycnocline and the thermocline when the latter are different as a result of salinity stratification. We present a revisited 2-degree resolution global climatology of monthly mean oceanic Barrier Layer (BL) thickness first proposed by de Boyer Montégut et al. (2007). In addition to using an extended data set, we present a modified computation method that addresses the observed porosity of BLs. We name porosity the fact that barrier layers distribution can, in some areas, be very uneven regarding the space and time scales that are considered. This implies an intermittent alteration of air-sea exchanges by the BL. Therefore, it may have important consequences for the climatic impact of BLs. Differences between the two computation methods are small for robust BLs that are formed by large-scale processes. However, the former approach can significantly underestimate the thickness of short and/or localized barrier layers. This is especially the case for barrier layers formed by mesoscale mechanisms (under the intertropical convergence zone for example and along western boundary currents) and equatorward of the sea surface salinity subtropical maxima. Complete characterisation of regional BL dynamics therefore requires a description of the robustness of BL distribution to assess the overall impact of BLs on the process of heat exchange between the ocean interior and the atmosphere.

An ensemble study of extreme storm surge related water levels in the North Sea in a changing climate - A. Sterl; H. van den Brink; H. de Vries; R. Haarsma; E. van Meijgaard (Ocean Science, 2009)
The height of storm surges is extremely important for a low-lying country like The Netherlands. By law, part of the coastal defence system has to withstand a water level that on average occurs only once every 10 000 years. The question then arises whether and how climate change affects the heights of extreme storm surges. Published research points to only small changes. However, due to the limited amount of data available results are usually limited to relatively frequent extremes like the annual 99%-ile. We here report on results from a 17-member ensemble of North Sea water levels spaning the period 1950–2100. It was created by forcing a surge model of the North Sea with meteorological output from a state-of-the-art global climate model which has been driven by greenhouse gas emissions following the SRES A1b scenario. The large ensemble size enables us to calculate 10 000 year return water levels with a low statistical uncertainty. In the one model used in this study, we find no statistically significant change in the 10 000 year return values of surge heights along the Dutch during the 21st century. Also a higher sea level resulting from global warming does not impact the height of the storm surges. As a side effect of our simulations we also obtain results on the interplay between surge and tide.

Metrics of hurricane-ocean interaction: vertically-integrated or vertically-averaged ocean temperature? - J. F. Price (Ocean Science, 2009)
The ocean thermal field is often represented in hurricane-ocean interaction by a metric termed upper Ocean Heat Content (OHC), the vertical integral of ocean temperature in excess of 26°C. High values of OHC have proven useful for identifying ocean regions that are especially favorable for hurricane intensification. Nevertheless, it is argued here that a more direct and robust metric of the ocean thermal field may be afforded by a vertical average of temperature. In the simplest version, dubbed T100, the averaging is from the surface to 100 m, a typical depth of vertical mixing by a category 3 hurricane. OHC and T100 are well correlated over the deep open ocean in the high range of OHC, ≥75 kJ cm−2. They are poorly correlated in the low range of OHC, ≤50 kJ cm−2, in part because OHC is degenerate when evaluated on cool ocean regions, ≤26°C. OHC and T100 can be qualitatively different also over shallow continental shelves: OHC will generally indicate comparatively low values regardless of the ocean temperature, while T100 will take on high values over a shelf that is warm and upwelling neutral or negative. In so far as the ocean thermal field alone is concerned, these warm, shallow continental shelves would appear to be as favorable for hurricane intensification as are warm, deep ocean regions.

Detecting marine hazardous substances and organisms: sensors for pollutants, toxins, and pathogens - O. Zielinski; J. A. Busch; A. D. Cembella; K. L. Daly; J. Engelbrektsson; A. K. Hannides; H. Schmidt (Ocean Science, 2009)
Marine environments are influenced by a wide diversity of anthropogenic and natural substances and organisms that may have adverse effects on human health and ecosystems. Real-time measurements of pollutants, toxins, and pathogens across a range of spatial scales are required to adequately monitor these hazards, manage the consequences, and to understand the processes governing their magnitude and distribution. Significant technological advancements have been made in recent years for the detection and analysis of such marine hazards. In particular, sensors deployed on a variety of mobile and fixed-point observing platforms provide a valuable means to assess hazards. In this review, we present state-of-the-art of sensor technology for the detection of harmful substances and organisms in the ocean. Sensors are classified by their adaptability to various platforms, addressing large, intermediate, or small areal scales. Current gaps and future demands are identified with an indication of the urgent need for new sensors to detect marine hazards at all scales in autonomous real-time mode. Progress in sensor technology is expected to depend on the development of small-scale sensor technologies with a high sensitivity and specificity towards target analytes or organisms. However, deployable systems must comply with platform requirements as these interconnect the three areal scales. Future developments will include the integration of existing methods into complex and operational sensing systems for a comprehensive strategy for long-term monitoring. The combination of sensor techniques on all scales will remain crucial for the demand of large spatial and temporal coverage.

Regional impacts of ocean color on tropical Pacific variability - W. Anderson; A. Gnanadesikan; A. Wittenberg (Ocean Science, 2009)
The role of the penetration length scale of shortwave radiation into the surface ocean and its impact on tropical Pacific variability is investigated with a fully coupled ocean, atmosphere, land and ice model. Previous work has shown that removal of all ocean color results in a system that tends strongly towards an El Niño state. Results from a suite of surface chlorophyll perturbation experiments show that the mean state and variability of the tropical Pacific is highly sensitive to the concentration and distribution of ocean chlorophyll. Setting the near-oligotrophic regions to contain optically pure water warms the mean state and suppresses variability in the western tropical Pacific. Doing the same above the shadow zones of the tropical Pacific also warms the mean state but enhances the variability. It is shown that increasing penetration can both deepen the pycnocline (which tends to damp El Niño) while shifting the mean circulation so that the wind response to temperature changes is altered. Depending on what region is involved this change in the wind stress can either strengthen or weaken ENSO variability.

The Kuroshio exchange with the South and East China Seas - T. Matsuno; J.-S. Lee; S. Yanao (Ocean Science, 2009)
The Kuroshio flows along the edges of the marginal East Asian seas such as the South China Sea (SCS) and East China Sea (ECS). Exchanges of materials and energy between the Kuroshio and the marginal seas partly control the environments of the marginal seas. In particular, saline water from the Kuroshio maintains certain salinity in the shelf water in the ECS. Nutrients from the subsurface of the Kuroshio may influence primary production on the shelf. We summarize how the Kuroshio comes into contact with the shelf water or marginal seas, describing phenomena related to the exchange between the Kuroshio and the ECS along with the SCS, using reports in the literature along with original data. The Kuroshio tends to intrude into the SCS through the Luzon Strait in various manners such as direct intrusion, associated with eddies and as a loop current. The Kuroshio intrusion into the shelf region of the ECS has distinct seasonal variation and the Taiwan Warm Current plays a significant role in the determination of water properties in the outer shelf associated with the Kuroshio intrusion. We then examine physical processes related to the interaction between the Kuroshio and shelf water. Interaction between the Kuroshio and the bottom topography is an important process in the control of the exchange around the shelf break. Vertical mixing and frontal eddies are also important factors that control the water exchange and formation of water masses in the outer shelf. Wind stress plays a significant role in the exchange with a rather event-like manner. To determine the source of the water masses, chemical tracers could be powerful tools and it is suggested that a significant part of the shelf water consists of Kuroshio intermediate water.

Frequency- or amplitude-dependent effects of the Atlantic meridional overturning on the tropical Pacific Ocean - G. J. van Oldenborgh; L. A. te Raa; H. A. Dijkstra; S. Y. Philip (Ocean Science, 2009)
Using the ECHAM5/MPI-OM model, we study the relation between the variations in the Atlantic meridional overturning circulation (AMOC) and both the Pacific sea surface temperature (SST) and the El Niño-Southern Oscillation (ENSO) amplitude. In a 17-member 20C3M/SRES-A1b ensemble for 1950–2100 the Pacific response to AMOC variations on different time scales and amplitudes is considered. The Pacific response to AMOC variations associated with the Atlantic Multidecadal Oscillation (AMO) is very small. In a 5-member hosing ensemble where the AMOC collapses due to a large freshwater anomaly, the Pacific SST response is large and in agreement with previous work. Our results show that the modelled connection between AMOC and ENSO depends very strongly on the frequency and/or the modelled amplitude of the AMOC variations. Interannual AMOC variations, decadal AMOC variations and an AMOC collapse lead to entirely different responses in the Pacific Ocean.

Precipitation of solid phase calcium carbonates and their effect on application of seawater S A – T – P models - G. M. Marion; F. J. Millero; R. Feistel (Ocean Science, 2009)
At the present time, little is known about how broad salinity and temperature ranges are for seawater thermodynamic models that are functions of absolute salinity (SA), temperature (T) and pressure (P). Such models rely on fixed compositional ratios of the major components (e.g., Na/Cl, Mg/Cl, Ca/Cl, SO4/Cl, etc.). As seawater evaporates or freezes, solid phases [e.g., CaCO3(s) or CaSO42H2O(s)] will eventually precipitate. This will change the compositional ratios, and these salinity models will no longer be applicable. A future complicating factor is the lowering of seawater pH as the atmospheric partial pressures of CO2 increase. A geochemical model (FREZCHEM) was used to quantify the SA−T boundaries at P=0.1 MPa and the range of these boundaries for future atmospheric CO2 increases. An omega supersaturation model for CaCO3 minerals based on pseudo-homogeneous nucleation was extended from 25–40°C to 3°C. CaCO3 minerals were the boundary defining minerals (first to precipitate) between 3°C (at SA=104 g kg−) and 40°C (at SA=66 g kg−). At 2.82°C, calcite(CaCO3) transitioned to ikaite(CaCO36H2O) as the dominant boundary defining mineral for colder temperatures, which culminated in a low temperature boundary of −4.93°C. Increasing atmospheric CO2 from 385 μatm (390 MPa) (in Year 2008) to 550 μatm (557 MPa) (in Year 2100) would increase the SA and t boundaries as much as 11 g kg−1 and 0.66°C, respectively. The model-calculated calcite-ikaite transition temperature of 2.82°C is in excellent agreement with ikaite formation in natural environments that occurs at temperatures of 3°C or lower. Furthermore, these results provide a quantitative theoretical explanation (FREZCHEM model calculation) for why ikaite is the solid phase CaCO3 mineral that precipitates during seawater freezing.

Understanding mixing efficiency in the oceans: do the nonlinearities of the equation of state for seawater matter? - R. Tailleux (Ocean Science, 2009)
There exist two central measures of turbulent mixing in turbulent stratified fluids that are both caused by molecular diffusion: 1) the dissipation rate D(APE) of available potential energy APE; 2) the turbulent rate of change Wr, turbulent of background gravitational potential energy GPEr. So far, these two quantities have often been regarded as the same energy conversion, namely the irreversible conversion of APE into GPEr, owing to the well known exact equality D(APE)=Wr, turbulent for a Boussinesq fluid with a linear equation of state. Recently, however, Tailleux (2009) pointed out that the above equality no longer holds for a thermally-stratified compressible, with the ratio ξ=Wr, turbulent/D(APE) being generally lower than unity and sometimes even negative for water or seawater, and argued that D(APE) and Wr, turbulent actually represent two distinct types of energy conversion, respectively the dissipation of APE into one particular subcomponent of internal energy called the "dead" internal energy IE0, and the conversion between GPEr and a different subcomponent of internal energy called "exergy" IEexergy. In this paper, the behaviour of the ratio ξ is examined for different stratifications having all the same buoyancy frequency N vertical profile, but different vertical profiles of the parameter Υ=α P/(ρCp), where α is the thermal expansion coefficient, P the hydrostatic pressure, ρ the density, and Cp the specific heat capacity at constant pressure, the equation of state being that for seawater for different particular constant values of salinity. It is found that ξ and Wr, turbulent depend critically on the sign and magnitude of dΥ/dz, in contrast with D(APE), which appears largely unaffected by the latter. These results have important consequences for how the mixing efficiency should be defined and measured in practice, which are discussed.

Dynamically constrained ensemble perturbations – application to tides on the West Florida Shelf - A. Barth; A. Alvera-Azcárate; J.-M. Beckers; R. H. Weisberg; L. Vandenbulcke; F. Lenartz; M. Rixen (Ocean Science, 2009)
A method is presented to create an ensemble of perturbations that satisfies linear dynamical constraints. A cost function is formulated defining the probability of each perturbation. It is shown that the perturbations created with this approach take the land-sea mask into account in a similar way as variational analysis techniques. The impact of the land-sea mask is illustrated with an idealized configuration of a barrier island. Perturbations with a spatially variable correlation length can be also created by this approach. The method is applied to a realistic configuration of the West Florida Shelf to create perturbations of the M2 tidal parameters for elevation and depth-averaged currents. The perturbations are weakly constrained to satisfy the linear shallow-water equations. Despite that the constraint is derived from an idealized assumption, it is shown that this approach is applicable to a non-linear and baroclinic model. The amplitude of spurious transient motions created by constrained perturbations of initial and boundary conditions is significantly lower compared to perturbing the variables independently or to using only the momentum equation to compute the velocity perturbations from the elevation.

Role of cabbeling in water densification in the Greenland Basin - Y. Kasajima; T. Johannessen (Ocean Science, 2009)
The effects of cabbeling mixing on water mass modification in the Greenland Sea were explored by hydrographic observations across the Greenland Basin in summer 2006. The neutral surface was chosen as a reference frame, and the strength of cabbeling mixing was quantified by the dianeutral velocity magnitude. Active cabbeling spots were detected with the criterion of the velocity magnitude >1 m/day, and four active cabbeling areas were identified; the west of Bear Island (SB), the Arctic Frontal Zone (AFZ), the central Greenland Sea (CG) and the western Greenland Sea (WG). The most vigorous cabbeling mixing was found at SB, where warm North Atlantic Water (NAW) mixed with cold water from the Barents Sea, inducing a maximum velocity of 7.5 m/day and a maximum density gain of 4.7×10−3 kg/m3. At AFZ and CG, the mixing took place between NAW, modified NAW and Arctic Intermediate Water (AIW), and the density gain at these fronts were 1.5×10−3 kg/m3 (AFZ) and 1.3×10−3 kg/m3 (CG). In the western Greenland Sea, the active cabbeling spots were widely separated and mixing appeared to be rather weak, with a maximum velocity of 2.5 m/day. The mixing source waters at WG were modified NAW, AIW and even denser water, and the density gain in this area was 0.4×10−3 kg/m3. The deepest mixing produced water whose density is equivalent to that of the dense water of the basin, indicating that cabbeling in the western Greenland Sea contributed directly to basin-scale water densification. The water mass modification rate was the highest at AFZ (about 8.0 Sv), suggesting that cabbeling may play an important role in water transformation in the Greenland Basin.

Freshwater components and transports in the Fram Strait – recent observations and changes since the late 1990s - B. Rabe; U. Schauer; A. Mackensen; M. Karcher; E. Hansen; A. Beszczynska-Möller (Ocean Science, 2009)
We present the late summer distribution and transports of freshwater components in the upper western part of the Fram Strait during 1998, 2004 and 2005. Hydrographic data and and water δ18O values are analyzed to distinguish Atlantic Water, ice melt (SIM) and freshwater removal from ice formation (IFB), and Meteoric Water (precipitation and riverine sources; MW). Concentrations of these water masses are combined with volume transport estimates from an inverse model. The average liquid freshwater transport relative to a reference salinity of 34.92, was 2500 km3/yr or 80 mSv southward, which is at the upper end of values reported in the literature. Our results indicate that not only the region of the continental slope but also parts of the East Greenland Shelf are important for freshwater transports. We estimate the average transports of of MW and IFB to be between 130 to 160 mSv (4100 to 5000 km3/yr) and 60 to 90 mSv (1900 to 2800 km3/yr) southward, respectively. The southward transport of MW was higher in 2005 than in 1998, but was compensated by a higher IFB transport. These differences in transports were associated with stronger southward velocities and the absence of northward velocities over the continental slope and the eastern East Greenland Shelf in 2005. A simulation using the North Atlantic-Arctic Ocean Sea Ice Model (NAOSIM) shows that the high transport of MW in the Fram Strait in 2005 is in agreement with the temporary storage of river water on the Siberian shelf in the mid-1990s, which reached the north of Greenland in 2003. Our results indicate that the accumulation of increased amounts of river water on the shelves is associated with enhanced ice formation.

Thermodynamic properties of standard seawater: extensions to high temperatures and pressures - J. Safarov; F. Millero; R. Feistel; A. Heintz; E. Hassel (Ocean Science, 2009)
Measurements of (p, ρ, T) properties of standard seawater with practical salinity S≈35, temperature T=(273.14 to 468.06) K and pressures, p, up to 140 MPa are reported with the reproducibility of the density measurements observed to be in the average percent deviation range Δρ/ρ=±(0.01 to 0.03)%. The measurements are made with a newly constructed vibration-tube densimeter which is calibrated using double-distilled water, methanol and aqueous NaCl solutions. Based on these and previous measurements, an empirical expression for the density of standard seawater has been developed as a function of pressure and temperature. This equation is used to calculate other volumetric properties including isothermal compressibility, isobaric thermal expansibility, differences in isobaric and isochoric heat capacities, the thermal pressure coefficient, internal pressure and the secant bulk modulus. The results can be used to extend the present equation of state of seawater to higher temperatures for pressure up to 140 MPa.

Geothermal heating, diapycnal mixing and the abyssal circulation - J. Emile-Geay; G. Madec (Ocean Science, 2009)
The dynamical role of geothermal heating in abyssal circulation is reconsidered using three independent arguments. First, we show that a uniform geothermal heat flux close to the observed average (86.4 mW m−2) supplies as much heat to near-bottom water as a diapycnal mixing rate of ~10−4 m2 s−1 – the canonical value thought to be responsible for the magnitude of the present-day abyssal circulation. This parity raises the possibility that geothermal heating could have a dynamical impact of the same order. Second, we estimate the magnitude of geothermally-induced circulation with the density-binning method (Walin, 1982), applied to the observed thermohaline structure of Levitus (1998). The method also allows to investigate the effect of realistic spatial variations of the flux obtained from heatflow measurements and classical theories of lithospheric cooling. It is found that a uniform heatflow forces a transformation of ~6 Sv at σ4=45.90, which is of the same order as current best estimates of AABW circulation. This transformation can be thought of as the geothermal circulation in the absence of mixing and is very similar for a realistic heatflow, albeit shifted towards slightly lighter density classes. Third, we use a general ocean circulation model in global configuration to perform three sets of experiments: (1) a thermally homogenous abyssal ocean with and without uniform geothermal heating; (2) a more stratified abyssal ocean subject to (i) no geothermal heating, (ii) a constant heat flux of 86.4 mW m−2, (iii) a realistic, spatially varying heat flux of identical global average; (3) experiments (i) and (iii) with enhanced vertical mixing at depth. Geothermal heating and diapycnal mixing are found to interact non-linearly through the density field, with geothermal heating eroding the deep stratification supporting a downward diffusive flux, while diapycnal mixing acts to map near-surface temperature gradients onto the bottom, thereby altering the density structure that supports a geothermal circulation. For strong vertical mixing rates, geothermal heating enhances the AABW cell by about 15% (2.5 Sv) and heats up the last 2000 m by ~0.15°C, reaching a maximum of by 0.3°C in the deep North Pacific. Prescribing a realistic spatial distribution of the heat flux acts to enhance this temperature rise at mid-depth and reduce it at great depth, producing a more modest increase in overturning than in the uniform case. In all cases, however, poleward heat transport increases by ~10% in the Southern Ocean. The three approaches converge to the conclusion that geothermal heating is an important actor of abyssal dynamics, and should no longer be neglected in oceanographic studies.

A new assessment of the error budget of global mean sea level rate estimated by satellite altimetry over 1993–2008 - M. Ablain; A. Cazenave; G. Valladeau; S. Guinehut (Ocean Science, 2009)
A new error budget assessment of the global Mean Sea Level (MSL) determined by TOPEX/Poseidon and Jason-1 altimeter satellites between January 1993 and June 2008 is presented using last altimeter standards. We discuss all potential errors affecting the calculation of the global MSL rate. We also compare altimetry-based sea level with tide gauge measurements over the altimetric period. Applying a statistical approach, this allows us to provide a realistic error budget of the MSL rise measured by satellite altimetry. These new calculations highlight a reduction in the rate of sea level rise since 2005, by ~2 mm/yr. This represents a 60% reduction compared to the 3.3 mm/yr sea level rise (glacial isostatic adjustment correction applied) measured between 1993 and 2005. Since November 2005, MSL is accurately measured by a single satellite, Jason-1. However the error analysis performed here indicates that the recent reduction in MSL rate is real.

Corrigendum to "A new method for forming approximately neutral surfaces" published in Ocean Sci., 5, 155–172, 2009 - A. Klocker; T. J. McDougall; D. R. Jackett (Ocean Science, 2009)

Evaluating two numerical advection schemes in HYCOM for eddy-resolving modelling of the Agulhas Current - B. C. Backeberg; L. Bertino; J. A. Johannessen (Ocean Science, 2009)
A 4th order advection scheme is applied in a nested eddy-resolving Hybrid Coordinate Ocean Model (HYCOM) of the greater Agulhas Current system for the purpose of testing advanced numerics as a means for improving the model simulation for eventual operational implementation. Model validation techniques comparing sea surface height variations, sea level skewness and variogram analyses to satellite altimetry measurements quantify that generally the 4th order advection scheme improves the realism of the model simulation. The most striking improvement over the standard 2nd order momentum advection scheme, is that the southern Agulhas Current is simulated as a well-defined meandering current, rather than a train of successive eddies. A better vertical structure and stronger poleward transports in the Agulhas Current core contribute toward a better southwestward penetration of the current, and its temperature field, implying a stronger Indo-Atlantic inter-ocean exchange. It is found that the transport, and hence this exchange, is sensitive to the occurrences of mesoscale features originating upstream in the Mozambique Channel and southern East Madagascar Current, and that the improved HYCOM simulation is well suited for further studies of these inter-actions.

A new method for forming approximately neutral surfaces - A. Klocker; T. J. McDougall; D. R. Jackett (Ocean Science, 2009)
We introduce a simple algorithm to improve existing density surfaces to ensure that the resulting surfaces are as close to neutral as possible. This means the slopes at any point on the surfaces are close to neutral tangent planes – the directions along which layered stirring and mixing occurs – minimizing the fictitious diapycnal diffusivity. Inverse techniques and layered models have been used for decades to understand ocean circulation. The most-used density surfaces are potential density or neutral density surfaces. Both these density surfaces and all others produce a fictitious diapycnal diffusivity to some degree due to the helical nature of neutral trajectories – with the magnitude of this artificial diffusivity in some cases being larger than the values measured in the ocean. Here we show how this error can be reduced by up to four orders of magnitude and therefore becomes insignificant compared to measured values, thus providing surfaces which would produce more accurate results when used for inverse techniques.

Deep currents in the Gulf of Guinea: along slope propagation of intraseasonal waves - C. Guiavarc'h; A. M. Treguier; A. Vangriesheim (Ocean Science, 2009)
In the Gulf of Guinea, intraseasonal variability is large at the equator and along the coast. Current data on the continental slope near 7.5° S show very energetic biweekly oscillations at 1300 m depth. A high resolution primitive equation numerical model demonstrates that this deep variability is forced by equatorial winds, through the generation of equatorial Yanai waves that propagate eastward and at depth, and then poleward as coastally-trapped waves upon reaching the coast of Africa. Intraseasonal variability is intensified along the coast of the Gulf of Guinea, especially in the 10–20 day period range and at depths between 500 and 1500 m. The kinetic energy distribution is well explained at first order by linear theory. Along the equator, eastward intensification of energy and bottom intensification are in qualitative agreement with vertically propagating Yanai waves, although the signal is influenced by the details of the bathymetry. Along the coast, baroclinic modes 3 to 5 are important close to the equator, and the signal is dominated by lower vertical modes farther south. Additional current meter data on the continental slope near 3° N display an energy profile in the 10–20 day period band that is strikingly different from the one at 7.5° S, with surface intensification rather than bottom intensification and a secondary maximum near 800 m. The model reproduces these features and explains them: the surface intensification in the north is due to the regional wind forcing, and the north-south asymmetry of the deep signal is due to the presence of the zonal African coast near 5° N. A 4 years time series of current measurements at 7.5° S shows that the biweekly oscillations are intermittent and vary from year to year. This intermittency is not well correlated with fluctuations of the equatorial winds and does not seem to be a simple linear response to the wind forcing.

Ekman layers in the Southern Ocean: spectral models and observations, vertical viscosity and boundary layer depth - S. Elipot; S. T. Gille (Ocean Science, 2009)
Spectral characteristics of the oceanic boundary-layer response to wind stress forcing are assessed by comparing surface drifter observations from the Southern Ocean to a suite of idealized models that parameterize the vertical flux of horizontal momentum using a first-order turbulence closure scheme. The models vary in their representation of vertical viscosity and boundary conditions. Each is used to derive a theoretical transfer function for the spectral linear response of the ocean to wind stress. The transfer functions are evaluated using observational data. The ageostrophic component of near-surface velocity is computed by subtracting altimeter-derived geostrophic velocities from observed drifter velocities (nominally drogued to represent motions at 15-m depth). Then the transfer function is computed to link these ageostrophic velocities to observed wind stresses. The traditional Ekman model, with infinite depth and constant vertical viscosity is among the worst of the models considered in this study. The model that most successfully describes the variability in the drifter data has a shallow layer of depth O(30–50 m), in which the viscosity is constant and O(100–1000 m2 s−1), with a no-slip bottom boundary condition. The second best model has a vertical viscosity with a surface value O(200 m2 s−1), which increases linearly with depth at a rate O(0.1–1 cm s−1) and a no-slip boundary condition at the base of the boundary layer of depth O(103 m). The best model shows little latitudinal or seasonal variability, and there is no obvious link to wind stress or climatological mixed-layer depth. In contrast, in the second best model, the linear coefficient and the boundary layer depth seem to covary with wind stress. The depth of the boundary layer for this model is found to be unphysically large at some latitudes and seasons, possibly a consequence of the inability of Ekman models to remove energy from the system by other means than shear-induced dissipation. However, the Ekman depth scale appears to scale like the climatological mixed-layer depth.

Molecular biology techniques and applications for ocean sensing - J. P. Zehr; I. Hewson; P. Moisander (Ocean Science, 2009)
The study of marine microorganisms using molecular biological techniques is now widespread in the ocean sciences. These techniques target nucleic acids which record the evolutionary history of microbes, and encode for processes which are active in the ocean today. Molecular techniques can form the basis of remote instrumentation sensing technologies for marine microbial diversity and ecological function. Here we review some of the most commonly used molecular biological techniques. These techniques include the polymerase chain reaction (PCR) and reverse-transcriptase PCR, quantitative PCR, whole assemblage "fingerprinting" approaches (based on nucleic acid sequence or length heterogeneity), oligonucleotide microarrays, and high-throughput shotgun sequencing of whole genomes and gene transcripts, which can be used to answer biological, ecological, evolutionary and biogeochemical questions in the ocean sciences. Moreover, molecular biological approaches may be deployed on ocean sensor platforms and hold promise for tracking of organisms or processes of interest in near-real time.

The density of seawater as a function of salinity (5 to 70 g kg −1 ) and temperature (273.15 to 363.15 K) - F. J. Millero; F. Huang (Ocean Science, 2009)
New seawater density measurements were made as a function of temperature T=(273.15 to 363.15) K and salinity (5 to 70 g kg−1). The measurements (N=230) from T=273.15 to 313.15 K and Practical Salinity (S) from 0 to 40 were found to be in good agreement (σ=0.0036 kg m−3) with the equation of state of seawater (Millero and Poisson, 1981) made on samples with a known chlorinity (Cl). These results indicate that the Practical Salinities (S) are in agreement to within ±0.003 kg m−3 with the values calculated from the Chlorinity, SCl=1.80655 Cl. The measurements from 298.15 to 363.15 were used to extend the equation of state to high temperatures and salinities. All the densities were made relative to pure water (ρ−ρ0 where ρ0 is the density for pure water) were fitted to equations of the form (ρ−ρ0)/(kg m−3)=A SA+B SA1.5+C SA2 where A, B, and C are functions of temperature and SA (g kg−1) is the absolute salinity, SA=(35.16504/35) S g kg−1. The fitted results from SA=0 to 50 g kg−1 and T=273.15 to 313.15 K (N=242) gave standard errors of 0.0037 kg m−3. The fitted results from 298.15 to 363.15 K (N=280) gave standard errors of 0.0063 kg m−3 and all the results (N=522) from 273.15 to 363.15 K gave standard errors of 0.0063 kg m−3. The earlier density measurements (Millero et al., 1976b; Poisson et al.,1980) used to determine the equation of state of seawater (Millero and Poisson, 1981) were combined to derive equations that are valid from 273.15 to 313.15 K and 273.15 to 363.15 K. The standard errors of these fits are respectively, 0.0038 kg m−3 (N=713) and 0.0063 kg m−3 (N=962). These new measurements expand the equation of state of seawater to a wider range of temperature (273.15 to 363.15) K and absolute salinity (0 to 70 kg m−3.

"Scale oriented" embedded modeling of the North-Western Mediterranean in the frame of MFSTEP - C. Estournel; F. Auclair; M. Lux; C. Nguyen; P. Marsaleix (Ocean Science, 2009)
An embedded forecasting system was developed for the North-Western Mediterranean at 3-km resolution. The system is based on the Symphonie hydrodynamic free surface model and on the variational initialization and forcing platform VIFOP. The regional model is initialized and forced at its open lateral boundaries by the MFS GCM and forced at the surface by the ALADIN numerical weather prediction model. Once a week, a five-day forecast is produced after a hindcast of seven days. This pre-modeling period of 7 days before beginning the forecast allows the development of the small scale features associated to the high resolution. The relevance of the 5-day forecast strategy has been examined by comparing the forecasted fields to hindcast fields (forced by meteorological and oceanic analyses) considered as a reference. Mismatches remain at a very low level indicating a good quality of the forecasted forcing fields and also possibly to the strong wind conditions which prevailed during the period. The embedded forecasts have been compared to the MFS observing system (SST and MedArgo) during 6 forecast cycles between September 2004 and February 2005. It was basically found that in the North-Western Mediterranean, the MFS basin-scale model and thus the regional model forecasts are characterized by large negative biases of salinity in the first 100 m under the surface leading thus to too light subsurface waters. The underestimation of temperature by the regional model just below the surface and its overestimation at 30m deep can be attributed to an overestimation of the turbulent mixing. The regional model allows to represent a number of processes especially those induced by the wind as coastal upwelling under stratified conditions, dense water formation over the Gulf of Lion shelf, deep mixing in the convection zone or influence on the Northern Current penetration in the Gulf of Lion.

Towards closure of regional heat budgets in the North Atlantic using Argo floats and surface flux datasets - N. C. Wells; S. A. Josey; R. E. Hadfield (Ocean Science, 2009)
The upper ocean heat budget (0–300 m) of the North Atlantic from 20°–60° N is investigated using data from Argo profiling floats for 1999–2005 and the NCEP/NCAR and NOC surface flux datasets. Estimates of the different terms in the budget (heat storage, advection, diffusion and surface exchange) are obtained using the methodology developed by Hadfield et al. (2007a, b). The method includes optimal interpolation of the individual profiles to produce gridded fields with error estimates at a 10°×10° grid box resolution. Closure of the heat budget is obtained within the error estimates for some regions – particularly the eastern subtropical Atlantic – but not for those boxes that include the Gulf Stream. Over the whole range considered, closure is obtained for 13 (9) out of 20 boxes with the NOC (NCEP/NCAR) surface fluxes. The seasonal heat budget at 20–30° N, 35–25° W is considered in detail. Here, the NCEP based budget has an annual mean residual of −55±35 Wm−2 compared with a NOC based value of −4±35 Wm−2. For this box, the net heat divergence of 36 Wm−2 (Ekman=−4 Wm−2, geostrophic=11 Wm−2, diffusion=29 Wm−2) offsets the net heating of 32 Wm−2 from the NOC surface heat fluxes. The results in this box are consistent with an earlier evaluation of the fluxes using measurements from research buoys in the subduction array which revealed biases in NCEP but good agreement of the buoy values with the NOC fields.

Turbulence closure: turbulence, waves and the wave-turbulence transition – Part 1: Vanishing mean shear - H. Z. Baumert; H. Peters (Ocean Science, 2009)
This paper extends a turbulence closure-like model for stably stratified flows into a new dynamic domain in which turbulence is generated by internal gravity waves rather than mean shear. The model turbulent kinetic energy (TKE, K) balance, its first equation, incorporates a term for the energy transfer from internal waves to turbulence. This energy source is in addition to the traditional shear production. The second variable of the new two-equation model is the turbulent enstrophy (Ω). Compared to the traditional shear-only case, the Ω-equation is modified to account for the effect of the waves on the turbulence time and space scales. This modification is based on the assumption of a non-zero constant flux Richardson number in the limit of vanishing mean shear when turbulence is produced exclusively by internal waves. This paper is part 1 of a continuing theoretical development. It accounts for mean shear- and internal wave-driven mixing only in the two limits of mean shear and no waves and waves but no mean shear, respectively. The new model reproduces the wave-turbulence transition analyzed by D'Asaro and Lien (2000b). At small energy density E of the internal wave field, the turbulent dissipation rate (ε) scales like ε~E2. This is what is observed in the deep sea. With increasing E, after the wave-turbulence transition has been passed, the scaling changes to ε~E1. This is observed, for example, in the highly energetic tidal flow near a sill in Knight Inlet. The new model further exhibits a turbulent length scale proportional to the Ozmidov scale, as observed in the ocean, and predicts the ratio between the turbulent Thorpe and Ozmidov length scales well within the range observed in the ocean.

Eastern Mediterranean biogeochemical flux model – Simulations of the pelagic ecosystem - G. Petihakis; G. Triantafyllou; K. Tsiaras; G. Korres; A. Pollani; I. Hoteit (Ocean Science, 2009)
During the second phase (2003–2006) of the Mediterranean ocean Forecasting System Project (MFS) named Toward Environmental Predictions (MFSTEP) one of the three major aims was the development of numerical forecasting systems. In this context a generic Biogeochemical Flux Model (BFM) was developed and coupled with hydrodynamic models already operating at basin scale as well as in regional areas. In the Eastern Mediterranean basin the BFM was coupled with the Aegean Levantine Eddy Resolving MOdel (ALERMO). The BFM is a generic highly complex model based on ERSEM and although a detailed description of the model and its components is beyond the scope of this work, a short overview of the main processes, laying emphasis on the parameter values used is presented. In addition the performance of the model is evaluated with some preliminary results being qualitatively compared against field observations. The model in its present form is rather promising and reproduces all important major features although there are some slight inefficiencies mostly related to primary and bacterial productivity rates.

On the time to tracer equilibrium in the global ocean - F. Primeau; E. Deleersnijder (Ocean Science, 2009)
An important issue for the interpretation of data from deep-sea cores is the time for tracers to be transported from the sea surface to the deep ocean. Global ocean circulation models can help shed light on the timescales over which a tracer comes to equilibrium in different regions of the ocean. In this note, we discuss how the most slowly decaying eigenmode of a model can be used to obtain a relevant timescale for a tracer that enters through the sea surface to become well mixed in the ocean interior. We show how this timescale depends critically on the choice between a Neumann surface boundary condition in which the flux of tracer is prescribed, a Robin surface boundary condition in which a combination of the flux and tracer concentration is prescribed or a Dirichlet surface boundary condition in which the concentration is prescribed. Explicit calculations with a 3-box model and a three-dimensional ocean circulation model show that the Dirichlet boundary condition when applied to only part of the surface ocean greatly overestimate the time needed to reach equilibrium. As a result regional-"injection" calculations which prescribe the surface concentration instead of the surface flux are not relevant for interpreting the regional disequilibrium between the Atlantic and Pacific found in paleo-tracer records from deep-sea cores. For tracers that enter the ocean through air-sea gas exchange a prescribed concentration boundary condition can be used to infer relevant timescales if the air-sea gas exchange rate is sufficiently fast, but the boundary condition must be applied over the entire ocean surface and not only to a patch of limited area. For tracers with a slow air-sea exchange rate such as 14C a Robin-type boundary condition is more relevant and for tracers such as δ18O that enter the ocean from melt water, a Neumann boundary condition is presumably more relevant. Our three-dimensional model results based on a steady-state modern circulation suggest that the relative disequilibrium between the deep Atlantic and Pacific is on the order of "only" 1200 years or less for a Neumann boundary condition and does not depend on the size and location of the patch where the tracer is injected.

The structure of the Persian Gulf outflow subjected to density variations - A. A. Bidokhti; M. Ezam (Ocean Science, 2009)
Oceanographic data and a dynamic model are used to consider the structure of Persian Gulf outflow. This outflow influences the physical properties of Oman seawater which appear in the CTD profiles of the Oman Sea. The observations show that thickness of the outflow, which is banked against the Oman coasts due to the earth rotation, is about 200 m with tongues extending east and north that may be due to the internal waves. A simple dynamical model of the outflow based on potential vorticity conservation is used to find the horizontal extension of the outflow from the coast. Typical mass transport estimate by the outflow is about 0.4 Sv, which is larger than those reported by others. This may be due to the fact the model is inviscid but the outflow is influenced by the bottom friction. Variability of the outflow structure may reflect the changing ecosystem of the Persian Gulf. Any change of the outflow source, the Persian Gulf Water (PGW), say salinity increase due to excessive evaporation (climate factor) or desalination (anthropogenic factors) of the PGW may change the outflow structure and the product waters in the Oman Sea. Hence, one can test different scenarios of changing the outflow source, the Persian Gulf Water (PGW), say by salinity increase due to excessive evaporation or desalination to estimate changes in the outflow structure and the product waters in the Oman Sea. The results of the model show that these can increase the outflow width and mass transport substantially.

Sequential assimilation of multi-mission dynamical topography into a global finite-element ocean model - S. Skachko; S. Danilov; T. Janjić; J. Schröter; D. Sidorenko; R. Savcenko; W. Bosch (Ocean Science, 2008)
This study focuses on an accurate estimation of ocean circulation via assimilation of satellite measurements of ocean dynamical topography into the global finite-element ocean model (FEOM). The dynamical topography data are derived from a complex analysis of multi-mission altimetry data combined with a referenced earth geoid. The assimilation is split into two parts. First, the mean dynamic topography is adjusted. To this end an adiabatic pressure correction method is used which reduces model divergence from the real evolution. Second, a sequential assimilation technique is applied to improve the representation of thermodynamical processes by assimilating the time varying dynamic topography. A method is used according to which the temperature and salinity are updated following the vertical structure of the first baroclinic mode. It is shown that the method leads to a partially successful assimilation approach reducing the rms difference between the model and data from 16 cm to 2 cm. This improvement of the mean state is accompanied by significant improvement of temporal variability in our analysis. However, it remains suboptimal, showing a tendency in the forecast phase of returning toward a free run without data assimilation. Both the mean difference and standard deviation of the difference between the forecast and observation data are reduced as the result of assimilation.

Retroflection from slanted coastlines-circumventing the "vorticity paradox" - V. Zharkov; D. Nof (Ocean Science, 2008)
The balance of long-shore momentum flux requires that the solution of zonally retroflecting currents involve ring shedding on the western side. An important aspect of the ring dynamics is the ring intensity α (analogous to the Rossby number), which reaches its maximum value of unity when the upstream potential vorticity (PV) is zero. Friction leads to a slow-down and a decrease in α. The main difficulty is that the solution of the system of equations for conservation of mass and momentum of zonal currents leads to the conclusion that the ratio (Φ) of the mass flux going into the rings and the total incoming mass flux is approximately 4α/(1+2α). This yields the "vorticity paradox" – only relatively weak rings (α≤1/2) could satisfy the necessary condition Φ≤1. Physically, this means, for example, that the momentum-flux of zero PV currents upstream is so high that, no matter how many rings are produced and, no matter what size they are, they cannot compensate for it. To avoid this paradox, we develop a nonlinear analytical model of retroflection from a slanted non-zonal coastline. We show that when the slant of coastline (γ) exceeds merely 150, Φ does not reach unity regardless of the value of α. Namely, the paradox disappears even for small slants. Our slowly varying nonlinear solution does not only let us circumvent the paradox. It also gives a detailed description of the rings growth rate and the mass flux going into the rings as a function of time. For example, in the case of zero PV and zero thickness of the upper layer along the coastline, the maximal values of Φ can be approximately expressed as, 1.012+0.32exp(−γ/3.41)−γ/225. Interestingly, for significant slants γ≥300), the rings reach a terminal size corresponding to a balance between the β-force and both the upstream and downstream momentum fluxes. This terminal size is unrelated to the ultimate detachment and westward drift due to β. Our analytical solutions are in satisfactory agreement with the results of a numerical model that we run.

Mutually consistent thermodynamic potentials for fluid water, ice and seawater: a new standard for oceanography - R. Feistel; D. G. Wright; K. Miyagawa; A. H. Harvey; J. Hruby; D. R. Jackett; T. J. McDougall; W. Wagner (Ocean Science, 2008)
A new seawater standard for oceanographic and engineering applications has been developed that consists of three independent thermodynamic potential functions, derived from extensive distinct sets of very accurate experimental data. The results have been formulated as Releases of the International Association for the Properties of Water and Steam, IAPWS (1996, 2006, 2008) and are expected to be adopted internationally by other organizations in subsequent years. In order to successfully perform computations such as phase equilibria from combinations of these potential functions, mutual compatibility and consistency of these independent mathematical functions must be ensured. In this article, a brief review of their separate development and ranges of validity is given. We analyse background details on the conditions specified at their reference states, the triple point and the standard ocean state, to ensure the mutual consistency of the different formulations, and the necessity and possibility of numerically evaluating metastable states of liquid water. Computed from this formulation in quadruple precision (128-bit floating point numbers), tables of numerical reference values are provided as anchor points for the consistent incorporation of additional potential functions in the future, and as unambiguous benchmarks to be used in the determination of numerical uncertainty estimates of double-precision implementations on different platforms that may be customized for special purposes.

Characterization of ASCAT measurements based on buoy and QuikSCAT wind vector observations - A. Bentamy; D. Croize-Fillon; C. Perigaud (Ocean Science, 2008)
The new scatterometer Advanced SCATterometer (ASCAT) onboard MetOp-A satellite provides surface wind speed and direction over global ocean with a spatial resolution of 25 km square over two swaths of 550 km widths. The accuracy of ASCAT wind retrievals is determined through various comparisons with moored buoys. The comparisons indicate that the remotely sensed wind speeds and directions agree well with buoy data. The root-mean-squared differences of the wind speed and direction are less than 1.72 m/s and 18°, respectively. At global scale, ASCAT winds are compared with surface winds derived from QuikSCAT scatterometer. The results confirm the buoy analyses, especially for wind speed ranging between 3 m/s and 20 m/s. For higher wind conditions, ASCAT is biased low. The ASCAT underestimation with respect to QuikSCAT winds is wind speed dependent. The comparisons based on the collocated scatterometer data collected after 17 of October 2007 indicate that there are significant improvements compared to previous periods.

Sensors for physical fluxes at the sea surface: energy, heat, water, salt - R. A. Weller; E. F. Bradley; J. B. Edson; C. W. Fairall; I. Brooks; M. J. Yelland; R. W. Pascal (Ocean Science, 2008)
The current status of meteorological sensors used aboard ships and buoys to measure the air-sea fluxes of momentum, heat, and freshwater is reviewed. Methods of flux measurement by the bulk aerodynamic, inertial dissipation and eddy-correlation methods are considered; and areas are identified where improvements are needed in measurement of the basic variables. In some cases, what is required is the transition from emergent to operational technology, in others new technologies are needed. Uncertainties in measured winds caused by flow distortion over the ship are discussed; and the possible role of computational fluid mechanics models to obtain corrections is considered. Basic studies are also needed on the influence of waves and rain on the fluxes. The issues involved in the specification of sea surface temperature are described, and the relative merits of the available sensors are discussed. The improved capability of buoy-mounted systems will depend on the emergence of low-power instruments, and/or new means of increasing the available power capacity. Other issues covered include the continuing uncertainty about the performance of rain gauges and short-wave radiometers. Also, the requirements for new instruments to extend the range of observations to extreme wind conditions are outlined, and the latest developments in the measurement of aerosol fluxes by eddy-correlation are presented.

Improving the parameterisation of horizontal density gradient in one-dimensional water column models for estuarine circulation - S. Blaise; E. Deleersnijder (Ocean Science, 2008)
A new parameterisation of horizontal density gradient for a one-dimensional water column estuarine model, inspired by the first-order finite-difference upwind scheme, is presented. This parameterisation prevents stratification from growing indefinitely, a deficiency usually referred to as "runaway stratification". It is seen that, using this upwind-like parameterisation, the salinity must remain comprised between upper and lower bounds set a priori and that any initial over- or under-shooting is progressively eliminated. Simulations of idealised and realistic estuarine regimes indicate that the new parameterisation lead to results that are devoid of the runaway stratification phenomenon, as opposed to previously used models.

Agulhas ring injection into the South Atlantic during glacials and interglacials - V. Zharkov; D. Nof (Ocean Science, 2008)
Recent proxies suggest that, at the end of the last glacial, there was a significant increase in the injection of Agulhas rings into the South Atlantic (SA). This brought about a dramatic increase in the salt-influx (from the Indian Ocean) into the SA helping re-start the then-collapsed meridional overturning cell (MOC), leading to the termination of the Younger Dryas (YD). Here, we propose a mechanism through which large variations in ring production take place. Using nonlinear analytical solutions for eddy shedding, we show that there are restricted possibilities for ring detachment when the coast is oriented in the north-south direction. We define a critical coastline angle below which there is rings shedding and above which there is almost no shedding. In the case of the Agulhas region, the particular shape of the African continent implies that rings can be produced only when the retroflection occurs beyond a specific latitude where the angle is critical. During glaciation, the wind stress curl (WSC) vanished at a latitude lower than that of the critical angle, which prohibited the retroflection from producing rings. When the latitude at which the WSC vanishes migrated poleward towards its present day position, the corresponding coastline angle decreased below the critical angle and allowed for a vigorous production of rings. Simple process-oriented numerical simulations (using the Bleck and Boudra model) are in satisfactory agreement with our results and enable us to affirm that, during the glacials, the behavior of the Agulhas Current (AC) was similar to that of the modern East Australian Current (EAC), for which the coastline slant is supercritical.

A comparison between vertical motions measured by ADCP and inferred from temperature data - H. van Haren (Ocean Science, 2008)
Combined vertical current (w) and thermistor string data demonstrate that high-, near-buoyancy frequency internal "wave" trains along a pycnocline in a flat-bottom shelf sea consist for 2 periods of a dominant mode-1 non-linear part, while thereafter mainly of linear [mode-2, quadrupled frequency] waves, to first order. In a simple [linear] heat budget the use of unfiltered temperature gradient or its time mean changes results by only 10%. The observations also demonstrate that temperature is not always adequate to estimate vertical motions using the linear 1-D heat equation. In shallow seas, tidal-w estimated from temperature data can be an order of magnitude weaker than directly observed w, and thus do not represent free internal waves. In the ocean, not too far from the main internal wave topography source, tidal motions represent linear waves and are well described by temperature-inferred w. There however, temperature-inferred w and directly observed w differ strongly near the buoyancy frequency, at which w is dominated by non-linear waves, and near [sub]inertial frequencies, at which w is dominated by eddies and gyroscopic waves.

Influence of Rossby waves on primary production from a coupled physical-biogeochemical model in the North Atlantic Ocean - G. Charria; I. Dadou; P. Cipollini; M. Drévillon; V. Garçon (Ocean Science, 2008)
Rossby waves appear to have a clear signature on surface chlorophyll concentrations which can be explained by a combination of vertical and horizontal mechanisms. In this study, we investigate the role of the different physical processes in the north Atlantic to explain the surface chlorophyll signatures and the consequences on primary production, using a 3-D coupled physical/biogeochemical model for the year 1998. The analysis at 20 given latitudes, mainly located in the subtropical gyre, where Rossby waves are strongly correlated with a surface chlorophyll signature, shows the important contribution of horizontal advection and of vertical advection and diffusion of inorganic dissolved nitrogen. The main control mechanism differs according to the biogeochemical background conditions of the area. The surface chlorophyll anomalies, induced by these physical mechanisms, have an impact on primary production. We estimate that Rossby waves induce, locally in space and time, increases (generally associated with the chlorophyll wave crest) and decreases (generally associated with the chlorophyll wave trough) in primary production, ~±20% of the estimated background primary production. This symmetrical situation suggests a net weak effect of Rossby waves on primary production.

On the Indonesian Throughflow in the OCCAM 1/4 degree ocean model - U. W. Humphries; D. J. Webb (Ocean Science, 2008)
The Indonesian Throughflow is analysed in two runs of the OCCAM 1/4 degree global ocean model, one using monthly climatological winds and one using ECMWF analysed six-hourly winds for the period 1993 to 1998. The long-term model throughflow agrees with observations and the value predicted by Godfrey's Island Rule. The Island Rule has some skill in predicting the annual signal each year but is poor at predicting year to year and shorter term variations in the total flow, especially in El Niño years. The spectra of transports in individual passages show significant differences between those connecting the region to the Pacific Ocean and those connecting with the Indian Ocean. On investigation we found that changes in the northern transports were strongly correlated with changes in the position of currents in the Celebes Sea and off Halmahera. Vertical profiles of transport are in reasonable agreement with observations but the model overestimates the near surface transport through the Lombok Strait and the dense overflow from the Pacific through the Lifamatola Strait into the deep Banda Sea. In both cases the crude representation of the passages by the model appears responsible. In the north the model shows, as expected, that the largest transport is via the Makassar Strait. However this is less than expected and instead there is significant flow via the Halmahera Sea. If Godfrey's Island Rule is correct and the throughflow is forced by the northward flow between Australia and South America, then the Halmahers Sea route should be important. It is the most southerly route around New Guinea to the Indian Ocean and there is no apparent reason why the flow should go further north in order to pass through the Makassar Strait. The model result thus raises the question of why in reality the Makassar Strait route appears to dominate the throughflow.

The low-resolution CCSM2 revisited: new adjustments and a present-day control run - M. Prange (Ocean Science, 2008)
The low-resolution (T31) version of the Community Climate System Model CCSM2.0.1 is revisited and adjusted by deepening the Greenland-Scotland ridge, changing oceanic mixing parameters, and applying a regional freshwater flux adjustment at high northern latitudes. The main purpose of these adjustments is to maintain a robust Atlantic meridional overturning circulation which collapses in the original model release. The paper describes the present-day control run of the adjusted model (referred to as "CCSM2/T31x3a") which is brought into climatic equilibrium by applying a deep-ocean acceleration technique. The accelerated integration is extended by a 100-year synchronous phase. The simulated meridional overturning circulation has a maximum of 14×106 m3 s−1 in the North Atlantic. The CCSM2/T31x3a control run is evaluated against observations and simulations with other climate models. Most shortcomings found in the CCSM2/T31x3a control run are identified as "typical problems" in global climate modelling. Finally, examples (simulation of North Atlantic hydrography, West African monsoon) are shown in which CCSM2/T31x3a has a better simulation skill than the latest low-resolution Community Climate System Model release, CCSM3/T31.

Observations of new western Mediterranean deep water formation using Argo floats 2004–2006 - R. O. Smith; H. L. Bryden; K. Stansfield (Ocean Science, 2008)
The deep convection that occurs in the western basin of the Mediterranean Sea was investigated using Argo float data over two consecutive winters in 2004–2005 and 2005–2006. The results showed deep mixed layers reaching 2000 m in surprising locations, namely the eastern Catalan subbasin (39.785° N, 4.845° E) and the western Ligurian subbasin (43.392° N, 7.765° E). Subsequently, new deep water was formed in March of 2005 and 2006 with θ=12.89–12.92°C, S=38.48–38.49 and σθ=29.113 kg m−3. The deep water produced in the Ligurian subbasin during 2006 was more saline, warmer and denser than any historical observations of western Mediterranean deep water. The results show S, θ and σθ in the western Mediterranean deep water are higher than 1990s values, with a salinity increase of 1.5×10−3 yr−1, a temperature increase of 3.6×10−3 °C yr−1 and a density increase of 4.0×10−4 kg m−3 yr−1 apparent from a dataset of western Mediterranean deep water properties spanning 1955–2006.

Descent and mixing of the overflow plume from Storfjord in Svalbard: an idealized numerical model study - I. Fer; B. Ådlandsvik (Ocean Science, 2008)
Storfjorden in the Svalbard Archipelago is a sill-fjord that produces significant volumes of dense, brine-enriched shelf water through ice formation. The dense water produced in the fjord overflows the sill and can reach deep into the Fram Strait. For conditions corresponding to a moderate ice production year, the pathway of the overflow, its descent and evolving water mass properties due to mixing are investigated for the first time using a high resolution 3-D numerical model. An idealized modeling approach forced by a typical annual cycle of buoyancy forcing due to ice production is chosen in a terrain-following vertical co-ordinate. Comparison with observational data, including hydrography, fine resolution current measurements and direct turbulence measurements using a microstructure profiler, gives confidence on the model performance. The model eddy diffusivity profiles contrasted to those inferred from the turbulence measurements give confidence on the skill of the Mellor Yamada scheme in representing sub-grid scale mixing for the Storfjorden overflow, and probably for gravity current modeling, in general. The Storfjorden overflow is characterized by low Froude number dynamics except at the shelf break where the plume narrows, accelerates with speed reaching 0.6 m s−1, yielding local Froude number in excess of unity. The volume flux of the plume increases by five-fold from the sill to downstream of the shelf-break. Rotational hydraulic control is not applicable for transport estimates at the sill using upstream basin information. To the leading order, geostrophy establishes the lateral slope of the plume interface at the sill. This allows for a transport estimate that is consistent with the model results by evaluating a weir relation at the sill.

Depth dependence of westward-propagating North Atlantic features diagnosed from altimetry and a numerical 1/6° model - A. Lecointre; T. Penduff; P. Cipollini; R. Tailleux; B. Barnier (Ocean Science, 2008)
A 1/6° numerical simulation is used to investigate the vertical structure of westward propagation between 1993 and 2000 in the North Atlantic ocean. The realism of the simulated westward propagating signals, interpreted principally as the signature of first-mode baroclinic Rossby waves (RW), is first assessed by comparing the simulated amplitude and zonal phase speeds of Sea Level Anomalies (SLA) against TOPEX/Poseidon-ERS satellite altimeter data. Then, the (unobserved) subsurface signature of RW phase speeds is investigated from model outputs by means of the Radon Transform which was specifically adapted to focus on first-mode baroclinic RW. The analysis is performed on observed and simulated SLA and along 9 simulated isopycnal displacements spanning the 0–3250 m depth range. Simulated RW phase speeds agree well with their observed counterparts at the surface, although with a slight slow bias. Below the surface, the simulated phase speeds exhibit a systematic deceleration with increasing depth, by a factor that appears to vary geographically. Thus, while the reduction factor is about 15–18% on average at 3250 m over the region considered, it appears to be much weaker (about 5–8%) in the eddy-active Azores Current, where westward propagating structures might be more coherent in the vertical. In the context of linear theories, these results question the often-made normal mode assumption of many WKB-based theories that the phase speed is independent of depth. Alternatively, these results could also suggest that the vertical structure of westward propagating signals may significantly depend on their degree of nonlinearity, with the degree of vertical coherence possibly increasing with the degree of nonlinearity.

Ice-shelf – ocean interactions at Fimbul Ice Shelf, Antarctica from oxygen isotope ratio measurements - M. R. Price; K. J. Heywood; K. W. Nicholls (Ocean Science, 2008)
Melt water from the floating ice shelves at the margins of the southeastern Weddell Sea makes a significant contribution to the fresh water budget of the region. In February 2005 a multi-institution team conducted an oceanographic campaign at Fimbul Ice Shelf on the Greenwich Meridian as part of the Autosub Under Ice programme. This included a mission of the autonomous submarine Autosub 25 km into the cavity beneath Fimbul Ice Shelf, and a number of ship-based hydrographic sections on the continental shelf and adjacent to the ice shelf front. The measurements reveal two significant sources of glacial melt water at Fimbul Ice Shelf: the main cavity under the ice shelf and an ice tongue, Trolltunga, that protrudes from the main ice front and out over the continental slope into deep water. Glacial melt water is concentrated in a 200 m thick Ice Shelf Water (ISW) layer below the base of the ice shelf at 150–200 m, with a maximum glacial melt concentration of up to 1.16%. Some glacial melt is found throughout the water column, and much of this is from sources other than Fimbul Ice Shelf. However, at least 0.2% of the water in the ISW layer cannot be accounted for by other processes and must have been contributed by the ice shelf. Just downstream of Fimbul Ice Shelf we observe locally created ISW mixing out across the continental slope. The ISW formed here is much less dense than that formed in the southwest Weddell Sea, and will ultimately contribute a freshening (and reduction in δ18O) to the upper 100–150 m of the water column in the southeast Weddell Sea.

Variability in the subtropical-tropical cells and its effect on near-surface temperature of the equatorial Pacific: a model study - J. F. Lübbecke; C. W. Böning; A. Biastoch (Ocean Science, 2008)
A set of experiments utilizing different implementations of the global ORCA-LIM model with horizontal resolutions of 2°, 0.5° and 0.25° is used to investigate tropical and extra-tropical influences on equatorial Pacific SST variability at interannual to decadal time scales. The model experiments use a bulk forcing methodology building on the global forcing data set for 1958 to 2000 developed by Large and Yeager (2004) that is based on a blend of atmospheric reanalysis data and satellite products. Whereas representation of the mean structure and transports of the (sub-) tropical Pacific current fields is much improved with the enhanced horizontal resolution, there is only little difference in the simulation of the interannual variability in the equatorial regime between the 0.5° and 0.25° model versions, with both solutions capturing the observed SST variability in the Niño3-region. The question of remotely forced oceanic contributions to the equatorial variability, in particular, the role of low-frequency changes in the transports of the Subtropical Cells (STCs), is addressed by a sequence of perturbation experiments using different combinations of fluxes. The solutions show the near-surface temperature variability to be governed by wind-driven changes in the Equatorial Undercurrent. The relative contributions of equatorial and off-equatorial atmospheric forcing differ between interannual and longer, (multi-) decadal timescales: for the latter there is a significant impact of changes in the equatorward transport of subtropical thermocline water associated with the lower branches of the STCs, related to variations in the off-equatorial trade winds. A conspicuous feature of the STC variability is that the equatorward transports in the interior and along the western boundary partially compensate each other at both decadal and interannual time scales, with the strongest transport extrema occurring during El Niño episodes. The behaviour is rationalized in terms of a wobbling in the poleward extents of the tropical gyres, which is manifested also in a meridional shifting of the bifurcation latitudes of the North and South Equatorial Current systems.

Operational ocean models in the Adriatic Sea: a skill assessment - J. Chiggiato; P. Oddo (Ocean Science, 2008)
In the framework of the Mediterranean Forecasting System (MFS) project, the performance of regional numerical ocean forecasting systems is assessed by means of model-model and model-data comparison. Three different operational systems considered in this study are: the Adriatic REGional Model (AREG); the Adriatic Regional Ocean Modelling System (AdriaROMS) and the Mediterranean Forecasting System General Circulation Model (MFS-GCM). AREG and AdriaROMS are regional implementations (with some dedicated variations) of POM and ROMS, respectively, while MFS-GCM is an OPA based system. The assessment is done through standard scores. In situ and remote sensing data are used to evaluate the system performance. In particular, a set of CTD measurements collected in the whole western Adriatic during January 2006 and one year of satellite derived sea surface temperature measurements (SST) allow to asses a full three-dimensional picture of the operational forecasting systems quality during January 2006 and to draw some preliminary considerations on the temporal fluctuation of scores estimated on surface quantities between summer 2005 and summer 2006. The regional systems share a negative bias in simulated temperature and salinity. Nonetheless, they outperform the MFS-GCM in the shallowest locations. Results on amplitude and phase errors are improved in areas shallower than 50 m, while degraded in deeper locations, where major models deficiencies are related to vertical mixing overestimation. In a basin-wide overview, the two regional models show differences in the local displacement of errors. In addition, in locations where the regional models are mutually correlated, the aggregated mean squared error was found to be smaller, that is a useful outcome of having several operational systems in the same region.

Ecological niches of three teuthophageous odontocetes in the northwestern Mediterranean Sea - E. Praca; A. Gannier (Ocean Science, 2008)
In the northwestern Mediterranean Sea, sperm whales, pilot whales and Risso's dolphins prey exclusively or preferentially on cephalopods. In order to evaluate their competition, we modelled their habitat suitability with the Ecological Niche Factor Analysis (ENFA) and compared their ecological niches using a discriminant analysis. We used a long term (1995–2005) small boat data set, with visual and acoustic (sperm whale) detections. Risso's dolphin had the shallowest and the more spatially restricted principal habitat, mainly located on the upper part of the continental slope (640 m mean depth). With a wider principal habitat, at 1750 m depth in average, the sperm whale used a deeper part of the slope as well as the closest offshore waters. Finally, the pilot whale has the most oceanic habitat (2500 m mean depth) mainly located in the central Ligurian Sea and Provençal basin. Therefore, potential competition for food between these species may be reduced by the differentiation of their habitats.

Operational ocean forecasting in the Eastern Mediterranean: implementation and evaluation - G. Zodiatis; R. Lardner; D. R. Hayes; G. Georgiou; S. Sofianos; N. Skliris; A. Lascaratos (Ocean Science, 2008)
The Cyprus Coastal Ocean Forecasting and Observing System (CYCOFOS) has been producing operational flow forecasts of the northeastern Levantine Basin since 2002 and has been substantially improved in 2005. CYCOFOS uses the POM flow model, and recently, within the frame of the MFSTEP project, the flow model was upgraded to use the hourly SKIRON atmospheric forcing, and its resolution was increased from 2.5 km to 1.8 km. The CYCOFOS model is now nested in the ALERMO regional model from the University of Athens, which is nested within the MFS basin model. The Variational Initialization and FOrcing Platform (VIFOP) has been implemented to reduce the numerical transient processes following initialization. Moreover, a five-day forecast is repeated every day, providing more detailed and more accurate information. Forecast results are posted on the web page http://www.oceanography.ucy.ac.cy/cycofos. The new, daily, high-resolution forecasts agree well with the ALERMO regional model. The agreement is better and results more reasonable when VIFOP is used. Active and slave experiments suggest that a four-week active period produces realistic results with more small-scale features. For runs in September 2004, biases with remote sensing sea surface temperature are less than 0.6°C with similar expressions of the flow field present in both. Remotely-observed coastal upwelling south of Cyprus and advection of cool water from the Rhodes Gyre to the southern shores of Cyprus are also modeled. In situ observed hydrographic data from south of Cyprus are similar to the corresponding forecast fields. Both indicate the relatively fresh subsurface Atlantic Water and a near-surface anticyclone south of Cyprus for August/September of 2004 and September 2005. Plans for further model improvement include assimilation of observed XBT temperature profiles, CTD profiles from drifters and gliders, and CT data from the CYCOFOS ocean observatory.

An Oceanographer's Guide to GOCE and the Geoid - C. W. Hughes; R. J. Bingham (Ocean Science, 2008)
A review is given of the geodetic concepts necessary for oceanographers to make use of satellite gravity data to define the geoid, and to interpret the resulting product. The geoid is defined, with particular attention to subtleties related to the representation of the permanent tide, and the way in which the geoid is represented in ocean models. The usual spherical harmonic description of the gravitational field is described, together with the concepts required to calculate a geoid from the spherical harmonic coefficients. A brief description is given of the measurement system in the GOCE satellite mission, scheduled for launch shortly. Finally, a recipe is given for calculation of the ocean dynamic topography, given a map of sea surface height above a reference ellipsoid, a set of spherical harmonic coefficients for the gravitational field, and defining constants.

A high-resolution free-surface model of the Mediterranean Sea - M. Tonani; N. Pinardi; S. Dobricic; I. Pujol; C. Fratianni (Ocean Science, 2008)
This study describes a new model implementation for the Mediterranean Sea with what is currently the highest vertical resolution over the Mediterranean basin. The resolution is of 1/16°×1/16° in the horizontal and has 72 unevenly spaced vertical levels. This model has been developed in the frame of the EU-MFSTEP project and is the operational forecast model currently used at the basin scale. The model considers an implicit free surface and this characteristic enhances the model's capability to simulate the sea surface height variability and the net transport at the Strait of Gibraltar. In this study we show the calibration/validation experiments performed before and after the model was used for forecasting. The first experiment consists of a six-year simulation forced by a perpetual year forcing, and the other experiment is a simulation from January 1997 to December 2004, forcing the model with 6-h atmospheric forcing fields from ECMWF. The model Sea Level Anomaly has been compared for the first time with satellite SLA and with ARGO data to provide evidence of the quality of the simulation. The results show that this model is capable of reproducing most of the variability of the general circulation in the Mediterranean Sea. However, some basic model inadequacies stand out and should be corrected in the near future.

Influence of numerical schemes on current-topography interactions in 1/4° global ocean simulations - T. Penduff; J. Le Sommer; B. Barnier; A.-M. Treguier; J.-M. Molines; G. Madec (Ocean Science, 2007)
The combined use of partial steps and of an energy-enstrophy conserving momentum advection scheme was shown by Barnier et al. (2006) to yield substantial improvements in the surface solution of the DRAKKAR ¼° global sea-ice/ocean model. The present study extends this investigation below the surface with a special focus on the Atlantic and reveals many improvements there as well: e.g. more realistic path, structure and transports of major currents (Gulf Stream, North Atlantic Current, Confluence region, Zapiola anticyclone), behavior of shedded rings, narrower subsurface boundary currents, stronger mean and eddy flows (MKE and EKE) at depth, beneficial enhancement of cyclonic (anticyclonic) flows around topographic depressions (mountains). Interestingly, adding a no-slip boundary condition to this improved model setup cancels most of these improvements, bringing back the biases diagnosed without the improved momentum advection scheme and partial steps (these biases are typical of other models at comparable or higher resolutions). This shows that current-topography interactions and full-depth eddy-admitting model solutions can be seriously deteriorated by near-bottom sidewall friction, either explicit or inherent to inadequate numerical schemes.

The backward Îto method for the Lagrangian simulation of transport processes with large space variations of the diffusivity - D. Spivakovskaya; A. W. Heemink; E. Deleersnijder (Ocean Science, 2007)
Random walk models are a powerful tool for the investigation of transport processes in turbulent flows. However, standard random walk methods are applicable only when the flow velocities and diffusivity are sufficiently smooth functions. In practice there are some regions where the rapid but continuous change in diffusivity may be represented by a discontinuity. The random walk model based on backward Îto calculus can be used for these problems. This model was proposed by LaBolle et al. (2000). The latter is best suited to the problems under consideration. It is then applied to two test cases with discontinuous diffusivity, highlighting the advantages of this method.

Southern Ocean overturning across streamlines in an eddying simulation of the Antarctic Circumpolar Current - A. M. Treguier; M. H. England; S. R. Rintoul; G. Madec; J. Le Sommer; J.-M. Molines (Ocean Science, 2007)
An eddying global model is used to study the characteristics of the Antarctic Circumpolar Current (ACC) in a streamline-following framework. Previous model-based estimates of the meridional circulation were calculated using zonal averages: this method leads to a counter-intuitive poleward circulation of the less dense waters, and underestimates the eddy effects. We show that on the contrary, the upper ocean circulation across streamlines agrees with the theoretical view: an equatorward mean flow partially cancelled by a poleward eddy mass flux. Two model simulations, in which the buoyancy forcing above the ACC changes from positive to negative, suggest that the relationship between the residual meridional circulation and the surface buoyancy flux is not as straightforward as assumed by the simplest theoretical models: the sign of the residual circulation cannot be inferred from the surface buoyancy forcing only. Among the other processes that likely play a part in setting the meridional circulation, our model results emphasize the complex three-dimensional structure of the ACC (probably not well accounted for in streamline-averaged, two-dimensional models) and the distinct role of temperature and salinity in the definition of the density field. Heat and salt transports by the time-mean flow are important even across time-mean streamlines. Heat and salt are balanced in the ACC, the model drift being small, but the nonlinearity of the equation of state cannot be ignored in the density balance.

Comparing the steric height in the Northern Atlantic with satellite altimetry - V. O. Ivchenko; S. D. Danilov; D. V. Sidorenko; J. Schröter; M. Wenzel; D. L. Aleynik (Ocean Science, 2007)
Anomalies of dynamic height derived from an analysis of Argo profiling buoys data are analysed to assess the relative roles of contributions from temperature and salinity over the North Atlantic for the period of 1999–2004. They are compared with dynamic topography anomalies based on TOPEX/Poseidon and Jason altimetry. It is shown that the halosteric contribution to the anomalies of dynamic height is comparable in magnitude to the thermosteric one over the period analyzed. Taking both salinity and temperature into account improves the agreement between zonally averaged trends in the satellite dynamic topography and dynamic height increasing the correlation between them to 0.73 from 0.63 when only temperature variability is taken into account. The implication of this result is that the salinity contribution cannot be neglected in the North Atlantic and that one cannot rely on estimating the thermosteric part by anomalies in the sea surface dynamic topography derived from the satellite altimetry.

Corrigendum to "Spectral studies of ocean water with space-borne sensor SCIAMACHY using Differential Optical Absorption Spectroscopy (DOAS)" published in Ocean Sci., 3, 429–440, 2007 - M. Vountas; T. Dinter; A. Bracher; J. P. Burrows; B. Sierk (Ocean Science, 2007)

Effects of mesoscale eddies on global ocean distributions of CFC-11, CO 2 , and Δ 14 C - Z. Lachkar; J. C. Orr; J.-C. Dutay; P. Delecluse (Ocean Science, 2007)
Global-scale tracer simulations are typically made at coarse resolution without explicitly modelling eddies. Here we ask what role do eddies play in ocean uptake, storage, and meridional transport of transient tracers. We made global anthropogenic transient-tracer simulations in coarse-resolution (2°cosφ×2°, ORCA2) and eddy-permitting (½°cosφ×½°, ORCA05) versions of the ocean general circulation model OPA9. Our focus is on surface-to-intermediate waters of the southern extratropics where air-sea tracer fluxes, tracer storage, and meridional tracer transport are largest. Eddies have little effect on global and regional bomb Δ14C uptake and storage. Yet for anthropogenic CO2 and CFC-11, refining the horizontal resolution reduced southern extratropical uptake by 25% and 28%, respectively. There is a similar decrease in corresponding inventories, which yields better agreement with observations. With higher resolution, eddies strengthen upper ocean vertical stratification and reduce excessive ventilation of intermediate waters by 20% between 60° S and 40° S. By weakening the residual circulation, i.e., the sum of Eulerian mean flow and the opposed eddy-induced flow, eddies reduce the supply of tracer-impoverished deep waters to the surface near the Antarctic divergence, thus reducing the air-sea tracer flux. Thus in the eddy permitting model, surface waters in that region have more time to equilibrate with the atmosphere before they are transported northward and subducted. As a result, the eddy permitting model's inventories of CFC-11 and anthropogenic CO2 are lower in that region because mixed-layer concentrations of both tracers equilibrate with the atmosphere on relatively short time scales (15 days and 6 months, respectively); conversely, bomb Δ14C's air-sea equilibration time of 6 years is so slow that, even in the eddy permitting model, there is little time for surface concentrations to equilibrate with the atmosphere, i.e., before surface waters are subducted.

Seasonal variability of phytoplankton fluorescence in relation to the Straits of Messina (Sicily) tidal upwelling - F. Azzaro; F. Decembrini; F. Raffa; E. Crisafi (Ocean Science, 2007)
In the Straits of Messina, large gradients of tidal displacements, as well as the topographic constrictions, determine the upwelling of deeper waters in the surface layer. This work describes the seasonal variability in the surface distribution of phytoplankton biomass depending on the upwelling phenomena. Temperature, salinity, nitrates and phytoplankton fluorescence were measured in 1994 and 1995 by continuous underway surface real-time measurements onboard dedicated research boats. Each survey was performed following the dynamic phases of flooding and ebbing tides. Tidal currents are essentially southward during high tide and northward during low tide. During the low water slack tide, large spatial gradients of physical-chemical and biological parameters were found, while at the high water slack tide, a diffused phytoplankton fluorescence was observed only in autumn due to a seasonal thermocline. Salinity, nitrate and chlorophyll-a fluorescence data revealed a significant positive intercorrelation, whereas they were inversely correlated with temperature. Generally, the upwelling distribution was limited to narrow zones during winter, while in summer it was found in the middle of the Straits and in the southern zones. During spring in the southern zone of the Straits, the maximum chlorophyll-a fluorescence was detected (May 1995, 0.32 μg-Chla l−1); in summer, when back and forth tidal movements between the Tyrrhenian and the Ionian seas intensify, decreased values were observed throughout the study area. The data set obtained through continuous and repeatable samplings has allowed the study of different time-space scales in the Straits of Messina, a very strong and dynamic environment. The Straits system could be compared to an "intermittent pump" which, during the different seasons, initially enriches itself and subsequently provides nutrients to the surrounding basins.

Internal tides and energy fluxes over Great Meteor Seamount - T. Gerkema; H. van Haren (Ocean Science, 2007)
Internal-tide energy fluxes are determined halfway over the southern slope of Great Meteor Seamount (Canary Basin), using data from combined CTD/LADCP yoyoing, covering the whole water column. The strongest signal is semi-diurnal and is concentrated in the upper few hundred meters of the water column. An indeterminacy in energy flux profiles is discussed; it is argued that a commonly applied condition used to determine these profiles is in fact invalid over sloping bottoms. However, the vertically integrated flux can be established unambiguously; the observed results are compared with the outcome of a numerical internal-tide generation model. For the semi-diurnal internal tide, the vertically integrated flux found in the model corresponds well to the observed one. The observed diurnal signal appears to be largely of non-tidal origin.

Spectral studies of ocean water with space-borne sensor SCIAMACHY using Differential Optical Absorption Spectroscopy (DOAS) - M. Vountas; T. Dinter; A. Bracher; J. P. Burrows; B. Sierk (Ocean Science, 2007)
Methods enabling the retrieval of oceanic parameter from the space borne instrumentation Scanning Imaging Absorption Spectrometer for Atmospheric ChartographY (SCIAMACHY) using Differential Optical Absorption Spectroscopy (DOAS) are presented. SCIAMACHY onboard ENVISAT measures back scattered solar radiation at a spectral resolution (0.2 to 1.5 nm). The DOAS method was used for the first time to fit modelled Vibrational Raman Scattering (VRS) in liquid water and in situ measured phytoplankton absorption reference spectra to optical depths measured by SCIAMACHY. Spectral structures of VRS and phytoplankton absorption were clearly found in these optical depths. Both fitting approaches lead to consistent results. DOAS fits correlate with estimates of chlorophyll concentrations: low fit factors for VRS retrievals correspond to large chlorophyll concentrations and vice versa; large fit factors for phytoplankton absorption correspond with high chlorophyll concentrations and vice versa. From these results a simple retrieval technique taking advantage of both measurements is shown. First maps of global chlorophyll concentrations were compared to the corresponding MODIS measurements with very promising results. In addition, results from this study will be used to improve atmospheric trace gas DOAS-retrievals from visible wavelengths by including these oceanographic signatures.

Atmospheric forcing of salinity in the overflow of Denmark Strait - J. Holfort; T. Albrecht (Ocean Science, 2007)
The temporal evolution of the characteristics of Denmark Strait Overflow Water (DSOW) is reconstructed using hydrographic data and compared with possible atmospheric forcing mechanisms. It is concluded that the main factor influencing the DSOW characteristics at a time scale of one to several years is the difference in mean sea level pressure across Denmark Strait or, in other words, the wind along Denmark Strait. At these time scales upstream changes in the characteristics of the different water masses involved in the formation of DSOW are only of minor importance. The main process responsible for the observed salinity changes in the DSOW is mixing in Denmark Strait. Triggered by the wind, different water masses contribute with changing amounts to the formation of DSOW, leading to the observed changes in the salinity of DSOW.

On the fast response of the Southern Ocean to changes in the zonal wind - D. J. Webb; B. A. de Cuevas (Ocean Science, 2007)
Model studies of the Southern Ocean, reported here, show that the Antarctic Circumpolar Current responds within two days to changes in the zonal wind stress at the latitudes of Drake Passage. Further investigation shows that the response is primarily barotropic and that, as one might expect, it is controlled by topography. Analysis of the results show that the changes in the barotropic flow are sufficient to transfer the changed surface wind stress to the underlying topography and that during this initial phase baroclinic processes are not involved. The model results also show that the Deacon Cell responds to changes in the wind stress on the same rapid time scale. It is shown that the changes in the Deacon Cell can also be explained by the change in the barotropic velocity field, an increase in the zonal wind stress producing an increased northward flow in shallow regions and southward flow where the ocean is deep. This new explanation is unexpected as previously the Deacon Cell has been thought of as a baroclinic feature of the ocean. The results imply that where baroclinic processes do appear to be involved in either the zonal momentum balance of the Southern Ocean or the formation of the Deacon Cell, they are part of the long term baroclinic response of the ocean's density field to the changes in the barotropic flow.

Assimilation of ocean colour data into a Biogeochemical Flux Model of the Eastern Mediterranean Sea - G. Triantafyllou; G. Korres; I. Hoteit; G. Petihakis; A. C. Banks (Ocean Science, 2007)
An advanced multivariate sequential data assimilation system has been implemented within the framework of the European MFSTEP project to fit a three-dimensional biogeochemical model of the Eastern Mediterranean to satellite chlorophyll data from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS). The physics are described by the Princeton Ocean Model (POM) while the biochemistry of the ecosystem is tackled with the Biogeochemical Flux Model (BFM). The assimilation scheme is based on the Singular Evolutive Extended Kalman (SEEK) filter, in which the error statistics were parameterized by means of a suitable set of Empirical Orthogonal Functions (EOFs). To avoid spurious long-range correlations associated with the limited number of EOFs, the filter covariance matrix was given compact support through a radius of influence around every data point location. Hindcast experiments were performed for one year over 1999 and forced with ECMWF 6 h atmospheric fields. The solution of the assimilation system was evaluated against the assimilated data and the MedAtlas climatology, and by assessing the impact of the assimilation on non-observed biogeochemical processes. It is found that the assimilation of SeaWiFS data improves the overall behavior of the BFM model and efficiently removes long term biases from the model despite some difficulties during the spring bloom period. Results, however, suggest the need of subsurface data to enhance the estimation of the ecosystem variables in the deep layers.

MedArgo: a drifting profiler program in the Mediterranean Sea - P.-M. Poulain; R. Barbanti; J. Font; A. Cruzado; C. Millot; I. Gertman; A. Griffa; A. Molcard; V. Rupolo; S. Le Bras; L. Petit de la Villeon (Ocean Science, 2007)
In the framework of the EU-funded MFSTEP project, autonomous drifting profilers were deployed throughout the Mediterranean Sea to collect temperature and salinity profile data and to measure subsurface currents. The realization of this profiler program in the Mediterranean, referred to as MedArgo, is described and assessed using data collected between June 2004 and December 2006 (including more than 2000 profiles). Recommendations are provided for the permanent future implementation of MedArgo in support of operational oceanography in the Mediterranean Sea. More than twenty drifting profilers were deployed from research vessels and ships-of-opportunity in most areas of the Mediterranean. They were all programmed to execute 5-day cycles with a drift at a parking depth of 350 m and CTD profiles from either 700 or 2000 m up to the surface. They stayed at the sea surface for about 6 h to be localised by, and transmit the data to, the Argos satellite system. The temperature and salinity data obtained with pumped Sea-Bird CTD instruments were processed and made available to the scientific community and to operational users in near-real time using standard Argo protocols, and were assimilated into Mediterranean numerical forecasting models. In general, the cycling and sampling characteristics chosen for the MedArgo profilers were found to be adequate for the Mediterranean. However, it is strongly advised to use GPS and global cellular phone telemetry or the future Argos bi-directional satellite system in order to avoid data compression and losses, for the continuation of the Mediterranean drifting profiler program.

Variability of Antarctic intermediate Water properties in the South Pacific Ocean - M. Tomczak (Ocean Science, 2007)
Argo float time series data are used to study the salinity field at the depth of the salinity minimum produced by Antarctic Intermediate Water (AAIW). It is found that far from showing the smooth erosion of the minimum that would result from diffusive flow, the salinity field is characterized by features of geostrophic turbulence such as fronts, eddies and intrusions. Comparison of the Argo float observations with the climatology of the World Ocean Atlas (WOA) reveals significant differences between the two data sets. Some of the differences may have their origin in problems with the WOA data density in remote regions of the South Pacific, but most are more likely produced by interannual variations of the AAIW salinity field.

Assessment of the importance of the current-wave coupling in the shelf ocean forecasts - G. Jordà; R. Bolaños; M. Espino; A. Sánchez-Arcilla (Ocean Science, 2007)
The effects of wave-current interactions on shelf ocean forecasts is investigated in the framework of the MFSTEP (Mediterranean Forecasting System Project Towards Enviromental Predictions) project. A one way sequential coupling approach is adopted to link the wave model (WAM) to the circulation model (SYMPHONIE). The coupling of waves and currents has been done considering four main processes: wave refraction due to currents, surface wind drag and bottom drag modifications due to waves, and the wave induced mass flux. The coupled modelling system is implemented in the southern Catalan shelf (NW Mediterranean), a region with characteristics similar to most of the Mediterranean shelves. The sensitivity experiments are run in a typical operational configuration. The wave refraction by currents seems to be not very relevant in a microtidal context such as the western Mediterranean. The main effect of waves on current forecasts is through the modification of the wind drag. The Stokes drift also plays a significant role due to its spatial and temporal characteristics. Finally, the enhanced bottom friction is just noticeable in the inner shelf.

Unpredictability of internal M 2 - H. van Haren (Ocean Science, 2007)
Current observations from a shelf sea, continental slopes and the abyssal North-East Atlantic Ocean are all dominated by the semidiurnal lunar (M2) tide. It is shown that motions at M2 vary at usually large barotropic and coherent baroclinic scales, >50 km horizontally and >0.5 H vertically. H represents the waterdepth. Such M2-scales are observed even close to topography, the potential source of baroclinic, "internal" tidal waves. In contrast, incoherent small-scale, ~10 km horizontally and ~0.1 H vertically, baroclinic motions are dominated around f, the local inertial frequency, and/or near 2Ω≈S2, the semidiurnal solar tidal frequency. Ω represents the Earth's rotational vector. This confirms earlier suggestions that small-scale baroclinic M2-motions generally do not exist in the ocean in any predictable manner, except in beams very near,

On the assimilation of ice velocity and concentration data into large-scale sea ice models - V. Dulière; T. Fichefet (Ocean Science, 2007)
Data assimilation into sea ice models designed for climate studies has started about 15 years ago. In most of the studies conducted so far, it is assumed that the improvement brought by the assimilation is straightforward. However, some studies suggest this might not be true. In order to elucidate this question and to find an appropriate way to further assimilate sea ice concentration and velocity observations into a global sea ice-ocean model, we analyze here results from a number of twin experiments (i.e. experiments in which the assimilated data are model outputs) carried out with a simplified model of the Arctic sea ice pack. Our objective is to determine to what degree the assimilation of ice velocity and/or concentration data improves the global performance of the model and, more specifically, reduces the error in the computed ice thickness. A simple optimal interpolation scheme is used, and outputs from a control run and from perturbed experiments without and with data assimilation are thoroughly compared. Our results indicate that, under certain conditions depending on the assimilation weights and the type of model error, the assimilation of ice velocity data enhances the model performance. The assimilation of ice concentration data can also help in improving the model behavior, but it has to be handled with care because of the strong connection between ice concentration and ice thickness. This study is first step towards real data assimilation into NEMO-LIM, a global sea ice-ocean model.