Consistency of the current global ocean observing systems from an Argo perspective - K. von Schuckmann; J.-B. Sallée; D. Chambers; P.-Y. Le Traon; C. Cabanes; F. Gaillard; S. Speich; M. Hamon (Ocean Science, 2014)
Variations in the world's ocean heat storage and its associated volume changes are a key factor to gauge global warming and to assess the earth's energy and sea level budget. Estimating global ocean heat content (GOHC) and global steric sea level (GSSL) with temperature/salinity data from the Argo network reveals a positive change of 0.5 ± 0.1 W m−2 (applied to the surface area of the ocean) and 0.5 ± 0.1 mm year−1 during the years 2005 to 2012, averaged between 60° S and 60° N and the 10–1500 m depth layer. In this study, we present an intercomparison of three global ocean observing systems: the Argo network, satellite gravimetry from GRACE and satellite altimetry. Their consistency is investigated from an Argo perspective at global and regional scales during the period 2005–2010. Although we can close the recent global ocean sea level budget within uncertainties, sampling inconsistencies need to be corrected for an accurate global budget due to systematic biases in GOHC and GSSL in the Tropical Ocean. Our findings show that the area around the Tropical Asian Archipelago (TAA) is important to closing the global sea level budget on interannual to decadal timescales, pointing out that the steric estimate from Argo is biased low, as the current mapping methods are insufficient to recover the steric signal in the TAA region. Both the large regional variability and the uncertainties in the current observing system prevent us from extracting indirect information regarding deep-ocean changes. This emphasizes the importance of continuing sustained effort in measuring the deep ocean from ship platforms and by beginning a much needed automated deep-Argo network.

Assessment of the structure and variability of Weddell Sea water masses in distinct ocean reanalysis products - T. S. Dotto; R. Kerr; M. M. Mata; M. Azaneu; I. Wainer; E. Fahrbach; G. Rohardt (Ocean Science, 2014)
We assessed and evaluated the performance of five ocean reanalysis products in reproducing essential hydrographic properties and their associated temporal variability for the Weddell Sea, Antarctica. The products used in this assessment were ECMWF ORAS4 (European Centre for Medium-Range Weather Forecasts Ocean Reanalysis System 4), CFSR (Climate Forecast System Reanalysis), MyOcean UR025.4 (University of Reading), ECCO2 (Estimating the Circulation and Climate of the Ocean, Phase II) and SODA (Simple Ocean Data Assimilation). The present study focuses on the Weddell Sea deep layer, which is composed of the following three main water masses: Warm Deep Water (WDW), Weddell Sea Deep Water (WSDW) and Weddell Sea Bottom Water (WSBW). The MyOcean UR025.4 product provided the most accurate representation of the structure and thermohaline properties of the Weddell Sea water masses when compared with observations. All the ocean reanalysis products analyzed exhibited limited capabilities in representing the surface water masses in the Weddell Sea. The CFSR and ECCO2 products were not able to represent deep water masses with a neutral density ≥ 28.40 kg m−3, which was considered the WSBW's upper limit throughout the simulation period. The expected WDW warming was only reproduced by the SODA product, whereas the ECCO2 product was able to represent the trends in the WSDW's hydrographic properties. All the assessed ocean reanalyses were able to represent the decrease in the WSBW's density, except the SODA product in the inner Weddell Sea. Improvements in parameterization may have as much impact on the reanalyses assessed as improvements in horizontal resolution primarily because the Southern Ocean lacks in situ data, and the data that are currently available are summer-biased. The choice of the reanalysis product should be made carefully, taking into account the performance, the parameters of interest, and the type of physical processes to be evaluated.

Extreme winter 2012 in the Adriatic: an example of climatic effect on the BiOS rhythm - M. Gačić; G. Civitarese; V. Kovačević; L. Ursella; M. Bensi; M. Menna; V. Cardin; P.-M. Poulain; S. Cosoli; G. Notarstefano; C. Pizzi (Ocean Science, 2014)
Adriatic and Ionian seas are Mediterranean sub-basins linked through the Bimodal Oscillating System mechanism responsible for decadal reversals of the Ionian basin-wide circulation. Altimetric maps showed that the last cyclonic mode started in 2011 but unexpectedly in 2012 reversed to anticyclonic. We related this "premature" inversion to the extremely strong winter in 2012, which caused the formation of very dense Adriatic waters, flooding Ionian flanks in May and inverting the bottom pressure gradient. Using Lagrangian float measurements, the linear regression between the sea surface height and three isopycnal depths suggests that the southward deep-layer flow coincided with the surface northward geostrophic current and the anticyclonic circulation regime. Density variations at depth in the northwestern Ionian revealed the arrival of Adriatic dense waters in May and maximum density in September. Comparison between the sea level height in the northwestern Ionian and in the basin centre showed that in coincidence with the arrival of the newly formed Adriatic dense waters the sea level was lowered in the northwestern flank, inverting the surface pressure gradient. Toward the end of 2012, the density gradient between the basin flanks and its centre went to zero, coinciding with the weakening of the anticyclonic circulation and eventually with its return to the cyclonic pattern. Thus, the premature and transient reversal of Ionian surface circulation originated from the extremely harsh winter in the Adriatic, resulting in the formation and spreading of highly dense bottom waters. The present study highlights the remarkable sensitiveness of the Adriatic–Ionian BiOS to climatic forcing.

Equilibrator-based measurements of dissolved nitrous oxide in the surface ocean using an integrated cavity output laser absorption spectrometer - I. Grefe; J. Kaiser (Ocean Science, 2014)
Dissolved nitrous oxide (N2O) concentrations are usually determined by gas chromatography (GC). Here we present laboratory tests and initial field measurements using a novel setup comprising a commercially available laser-based analyser for N2O, carbon monoxide and water vapour coupled to a glass-bed equilibrator. This approach is less labour-intensive and provides higher temporal and spatial resolution than the conventional GC technique. The standard deviation of continuous equilibrator or atmospheric air measurements was 0.2 nmol mol−1 (averaged over 5 min). The short-term repeatability for reference gas measurements within 1 h of each other was 0.2 nmol mol−1 or better. Another indicator of the long-term stability of the analyser is the standard deviation of the calibrated N2O mole fraction in marine air, which was between 0.5 and 0.7 nmol mol−1. The equilibrator measurements were compared with purge-and-trap gas chromatography–mass spectrometry (GC-MS) analyses of N2O concentrations in discrete samples from the Southern Ocean and showed agreement to within the 2% measurement uncertainty of the GC-MS method. The equilibrator response time to concentration changes in water was from 142 to 203 s, depending on the headspace flow rate. The system was tested at sea during a north-to-south transect of the Atlantic Ocean. While the subtropical gyres were slightly undersaturated, the equatorial region was a source of nitrous oxide to the atmosphere, confirming previous findings (Forster et al., 2009). The ability to measure at high temporal and spatial resolution revealed submesoscale variability in dissolved N2O concentrations. Mean sea-to-air fluxes in the tropical and subtropical Atlantic ranged between −1.6 and 0.11 μmol m−2 d−1 and confirm that the subtropical Atlantic is not an important source region for N2O to the atmosphere, compared to global average fluxes of 0.6–2.4 μmol m−2 d−1. The system can be easily modified for autonomous operation on voluntary observing ships (VOS). Future work should include an interlaboratory comparison exercise with other methods of dissolved N2O analyses.

The land-ice contribution to 21st-century dynamic sea level rise - T. Howard; J. Ridley; A. K. Pardaens; R. T. W. L. Hurkmans; A. J. Payne; R. H. Giesen; J. A. Lowe; J. L. Bamber; T. L. Edwards; J. Oerlemans (Ocean Science, 2014)
Climate change has the potential to influence global mean sea level through a number of processes including (but not limited to) thermal expansion of the oceans and enhanced land ice melt. In addition to their contribution to global mean sea level change, these two processes (among others) lead to local departures from the global mean sea level change, through a number of mechanisms including the effect on spatial variations in the change of water density and transport, usually termed dynamic sea level changes. In this study, we focus on the component of dynamic sea level change that might be given by additional freshwater inflow to the ocean under scenarios of 21st-century land-based ice melt. We present regional patterns of dynamic sea level change given by a global-coupled atmosphere–ocean climate model forced by spatially and temporally varying projected ice-melt fluxes from three sources: the Antarctic ice sheet, the Greenland Ice Sheet and small glaciers and ice caps. The largest ice melt flux we consider is equivalent to almost 0.7 m of global mean sea level rise over the 21st century. The temporal evolution of the dynamic sea level changes, in the presence of considerable variations in the ice melt flux, is also analysed. We find that the dynamic sea level change associated with the ice melt is small, with the largest changes occurring in the North Atlantic amounting to 3 cm above the global mean rise. Furthermore, the dynamic sea level change associated with the ice melt is similar regardless of whether the simulated ice fluxes are applied to a simulation with fixed CO2 or under a business-as-usual greenhouse gas warming scenario of increasing CO2.

Sources of 21st century regional sea-level rise along the coast of northwest Europe - T. Howard; A. K. Pardaens; J. L. Bamber; J. Ridley; G. Spada; R. T. W. L. Hurkmans; J. A. Lowe; D. Vaughan (Ocean Science, 2014)
Changes in both global and regional mean sea level, and changes in the magnitude of extreme flood heights, are the result of a combination of several distinct contributions most, but not all, of which are associated with climate change. These contributions include effects in the solid earth, gravity field, changes in ocean mass due to ice loss from ice sheets and glaciers, thermal expansion, alterations in ocean circulation driven by climate change and changing freshwater fluxes, and the intensity of storm surges. Due to the diverse range of models required to simulate these systems, the contributions to sea-level change have usually been discussed in isolation rather than in one self-consistent assessment. Focusing on the coastline of northwest Europe, we consider all the processes mentioned above and their relative impact on 21st century regional mean sea levels and the 50-year return flood height. As far as possible our projections of change are derived from process-based models forced by the A1B emissions scenario to provide a self-consistent comparison of the contributions. We address uncertainty by considering both a mid-range and an illustrative high-end combination of the different components. For our mid-range ice loss scenario we find that thermal expansion of seawater is the dominant contributor to change in northwest European sea level by 2100. However, the projected contribution to extreme sea level, due to changes in storminess alone, is in some places significant and comparable to the global mean contribution of thermal expansion. For example, under the A1B emissions scenario, by 2100, change in storminess contributes around 15 cm to the increase in projected height of the 50-year storm surge on the west coast of the Jutland Peninsula, compared with a contribution of around 22 cm due to thermal expansion and a total of 58 cm from all of the contributions we consider. An illustrative combination of our high-end projections suggests increases in the 50-year return level of 86 cm at Sheerness, 95 cm at Roscoff, 106 cm at Esbjerg, and 67cm at Bergen. The notable regional differences between these locations arise primarily from differences in the rates of vertical land movement and changes in storminess.

Modelling seasonal circulation and thermohaline structure of the Caspian Sea - M. Gunduz; E. Özsoy (Ocean Science, 2014)
The wind- and buoyancy-driven seasonal circulation of the Caspian Sea is investigated for a better understanding of its basin-wide and mesoscale dynamics, mixing and transport. The model successfully reproduces the following basic elements of the circulation: the southward-flowing current systems along the eastern and western coasts, the upwelling along the eastern coast, the cyclonic circulation in the Middle Caspian Sea (MCS), especially in winter, and the cyclonic and anticyclonic cells of circulation in the South Caspian Sea (SCS). The observed seasonal variability of sea level and sea surface temperature (SST) is well reproduced. Mesoscale structures, not always evident from hydrographic observations, are identified.

Ventilation of the Mediterranean Sea constrained by multiple transient tracer measurements - T. Stöven; T. Tanhua (Ocean Science, 2014)
Ventilation is the primary pathway for atmosphere–ocean boundary perturbations, such as temperature anomalies, to be relayed to the ocean interior. It is also a conduit for gas exchange between the interface of atmosphere and ocean. Thus it is a mechanism whereby, for instance, the ocean interior is oxygenated and enriched in anthropogenic carbon. The ventilation of the Mediterranean Sea is fast in comparison to the world ocean and has large temporal variability. Here we present transient tracer data from a field campaign in April 2011 that sampled a unique suite of transient tracers (SF6, CFC-12, 3H and 3He) in all major basins of the Mediterranean. We apply the transit time distribution (TTD) model to the data in order to constrain the mean age, the ratio of the advective / diffusive transport and the number of water masses significant for ventilation. We found that the eastern part of the eastern Mediterranean can be reasonably described with a one-dimensional inverse Gaussian TTD (IG-TTD), and thus constrained with two independent tracers. The ventilation of the Ionian Sea and the western Mediterranean can only be constrained by a linear combination of IG-TTDs. We approximate the ventilation with a one-dimensional, two inverse Gaussian TTD (2IG-TTD) for these areas and demonstrate a possibility of constraining a 2IG-TTD from the available transient tracer data. The deep water in the Ionian Sea has a mean age between 120 and 160 years and is therefore substantially older than the mean age of the Levantine Basin deep water (60–80 years). These results are in contrast to those expected by the higher transient tracer concentrations in the Ionian Sea deep water. This is partly due to deep water of Adriatic origin having more diffusive properties in transport and formation (i.e., a high ratio of diffusion over advection), compared to the deep water of Aegean Sea origin that still dominates the deep Levantine Basin deep water after the Eastern Mediterranean Transient (EMT) in the early 1990s. The tracer minimum zone (TMZ) in the intermediate of the Levantine Basin is the oldest water mass with a mean age up to 290 years. We also show that the deep western Mediterranean has contributed approximately 40% of recently ventilated deep water from the Western Mediterranean Transition (WMT) event of the mid-2000s. The deep water has higher transient tracer concentrations than the mid-depth water, but the mean age is similar with values between 180 and 220 years.

Fine-scale features on the sea surface in SAR satellite imagery – Part 2: Numerical modeling - S. Matt; A. Fujimura; A. Soloviev; S. H. Rhee; R. Romeiser (Ocean Science, 2014)
With the advent of the new generation of synthetic aperture radar (SAR) satellites, it has become possible to resolve fine-scale features on the sea surface on the scale of meters. The proper identification of sea surface signatures in SAR imagery can be challenging, since some features may be due to atmospheric distortions (gravity waves, squall lines) or anthropogenic influences (slicks), and may not be related to dynamic processes in the upper ocean. In order to improve our understanding of the nature of fine-scale features on the sea surface and their signature in SAR, we have conducted high-resolution numerical simulations combining a three-dimensional non-hydrostatic computational fluid dynamics model with a radar imaging model. The surface velocity field from the hydrodynamic model is used as input to the radar imaging model. The combined approach reproduces the sea surface signatures in SAR of ship wakes, low-density plumes, and internal waves in a stratified environment. The numerical results are consistent with observations reported in a companion paper on in situ measurements during SAR satellite overpasses. Ocean surface and internal waves are also known to produce a measurable signal in the ocean magnetic field. This paper explores the use of computational fluid dynamics to investigate the magnetic signatures of oceanic processes. This potentially provides a link between SAR signatures of transient ocean dynamics and magnetic field fluctuations in the ocean. We suggest that combining SAR imagery with data from ocean magnetometers may be useful as an additional maritime sensing method. The new approach presented in this work can be extended to other dynamic processes in the upper ocean, including fronts and eddies, and can be a valuable tool for the interpretation of SAR images of the ocean surface.

On the tides and resonances of Hudson Bay and Hudson Strait - D. J. Webb (Ocean Science, 2014)
The resonances of Hudson Bay, Foxe Basin and Hudson Strait are investigated using a linear shallow water numerical model. The region is of particular interest because it is the most important region of the world ocean for dissipating tidal energy. The model shows that the semi-diurnal tides of the region are dominated by four nearby overlapping resonances. It shows that these not only affect Ungava Bay, a region of extreme tidal range, but they also extend far into Foxe Basin and Hudson Bay and appear to be affected by the geometry of those regions. The results also indicate that it is the four resonances acting together which make the region such an important area for dissipating tidal energy.

The effects of global climate change on the cycling and processes of persistent organic pollutants (POPs) in the North Sea - K. O'Driscoll; B. Mayer; J. Su; M. Mathis (Ocean Science, 2014)
The fate and cycling of two selected legacy persistent organic pollutants (POPs), PCB 153 and γ-HCH, in the North Sea in the 21st century have been modelled with combined hydrodynamic and fate and transport ocean models (HAMSOM and FANTOM, respectively). To investigate the impact of climate variability on POPs in the North Sea in the 21st century, future scenario model runs for three 10-year periods to the year 2100 using plausible levels of both in situ concentrations and atmospheric, river and open boundary inputs are performed. This slice mode under a moderate scenario (A1B) is sufficient to provide a basis for further analysis. For the HAMSOM and atmospheric forcing, results of the IPCC A1B (SRES) 21st century scenario are utilized, where surface forcing is provided by the REMO downscaling of the ECHAM5 global atmospheric model, and open boundary conditions are provided by the MPIOM global ocean model. Dry gas deposition and volatilization of γ-HCH increase in the future relative to the present by up to 20% (in the spring and summer months for deposition and in summer for volatilization). In the water column, total mass of γ-HCH and PCB 153 remain fairly steady in all three runs. In sediment, γ-HCH increases in the future runs, relative to the present, while PCB 153 in sediment decreases exponentially in all three runs, but even faster in the future, due to the increased number of storms, increased duration of gale wind conditions and increased water and air temperatures, all of which are the result of climate change. Annual net sinks exceed sources at the ends of all periods. Overall, the model results indicate that the climate change scenarios considered here generally have a negligible influence on the simulated fate and transport of the two POPs in the North Sea, although the increased number and magnitude of storms in the 21st century will result in POP resuspension and ensuing revolatilization events. Trends in emissions from primary and secondary sources will remain the key driver of levels of these contaminants over time.

Evaluation of MERIS products from Baltic Sea coastal waters rich in CDOM - J. M. Beltrán-Abaunza; S. Kratzer; C. Brockmann (Ocean Science, 2014)
In this study, retrievals of the medium resolution imaging spectrometer (MERIS) reflectances and water quality products using four different coastal processing algorithms freely available are assessed by comparison against sea-truthing data. The study is based on a pair-wise comparison using processor-dependent quality flags for the retrieval of valid common macro-pixels. This assessment is required in order to ensure the reliability of monitoring systems based on MERIS data, such as the Swedish coastal and lake monitoring system (http://vattenkvalitet.se). The results show that the pre-processing with the Improved Contrast between Ocean and Land (ICOL) processor, correcting for adjacency effects, improves the retrieval of spectral reflectance for all processors. Therefore, it is recommended that the ICOL processor should be applied when Baltic coastal waters are investigated. Chlorophyll was retrieved best using the FUB (Free University of Berlin) processing algorithm, although overestimations in the range 18–26.5%, dependent on the compared pairs, were obtained. At low chlorophyll concentrations (

Numerical modelling of physical processes governing larval transport in the southern North Sea - M. C. H. Tiessen; L. Fernard; T. Gerkema; J. van der Molen; P. Ruardij; H. W. van der Veer (Ocean Science, 2014)
A three-dimensional hydrodynamic model (GETM) was coupled with a particle tracking routine (GITM) to study the inter-annual variability in transport paths of particles in the North Sea and English Channel. For validation, a comparison with observed drifter trajectories is also presented here. This research investigated to what extent variability in the hydrodynamic conditions alone (reflecting passive particle transport) contributed to inter-annual variability in the transport of eggs and larvae. In this idealised study, no a priori selection of specific spawning grounds or periods was made and no active behaviour (vertical migration) or mortality was included. In this study, egg and larval development towards coastal nursery areas was based solely on sea water temperature, while settlement areas were defined by a threshold water depth. Results showed strong inter-annual variability in drift direction and distance, caused by a combination of wind speed and direction. Strong inter-annual variability was observed both in absolute amount of settlement in several coastal areas, and in the relative importance of the different areas. The effects of wind and temperature variability are minor for settlement along the western shores of the North Sea and in the English Channel, but have a very significant impact on settlement along the eastern shores of the North Sea. Years with strong south-westerly winds across the Dover Straight resulted in higher settlement figures along its eastern shores of the North Sea (standard deviation 37% of the mean annual settlement value). Settlement in the western Dutch Wadden Sea did not only show inter-annual variability, but patterns were also variable within each year and revealed seasonal changes in the origin of particles: during winter, stronger currents along with colder temperatures generally result in particles originating from further away.

Interannual correlations between sea surface temperature and concentration of chlorophyll pigment off Punta Eugenia, Baja California, during different remote forcing conditions - H. Herrera-Cervantes; S. E. Lluch-Cota; D. B. Lluch-Cota; G. Gutiérrez-de-Velasco (Ocean Science, 2014)
Interannual correlation between satellite-derived sea surface temperature (SST) and surface chlorophyll a (Chl a) are examined in the coastal upwelling zone off Punta Eugenia on the west coast of the Baja California Peninsula, an area than has been identified as having intense biological productivity and oceanographic transition between midlatitude and tropical ocean conditions. We used empirical orthogonal functions (EOF) analysis separately and jointly on the two fields from 1997 through 2007, a time period dominated by different remote forcing: ENSO (El Niño–Southern Oscillation) conditions (weak, moderate and strong) and the largest intrusion of subarctic water reported in the last 50 years. Coastal upwelling index anomalies (CUI) and the multivariate ENSO index (MEI) were used to identify the influence of local (wind stress) and remote (ENSO) forcing over the interannual variability of both variables. The spatial pattern of the individual EOF1 analysis showed the greater variability of SST and Chl a offshore, their corresponding amplitude time series presented the highest peaks during the strong 1997–2000 El Niño–La Niña cycles and during the 2002–2004 period associated to the intrusion of subarctic water. The MEI is well correlated with the individual SST principal component (R ≈ 0.67, P 0.4) mainly regulated by ENSO cycles. This was spatially revealed when we calculated the homogeneous correlations for the 1997–1999 El Niño–La Niña period and during the 2002–2004 period, the intrusion of subarctic water period. Both, SST and Chl a showed higher coupling and two distinct physical–biological responses: on average ENSO influence was observed clearly along the coast mostly in SST, while the subarctic water influence, observed offshore and in Bahía Vizcaíno, mostly in Chl a. We found coastal chlorophyll blooms off Punta Eugenia during the 2002–2003 period, an enrichment pattern similar to that observed off the coast of Oregon. These chlorophyll blooms are likely linked to high wind stress anomalies during 2002, mainly at high latitudes. This observation may provide an explanation of why Punta Eugenia is one of the most important biological action centers on the Pacific coast.

Adapting to life: ocean biogeochemical modelling and adaptive remeshing - J. Hill; E. E. Popova; D. A. Ham; M. D. Piggott; M. Srokosz (Ocean Science, 2014)
An outstanding problem in biogeochemical modelling of the ocean is that many of the key processes occur intermittently at small scales, such as the sub-mesoscale, that are not well represented in global ocean models. This is partly due to their failure to resolve sub-mesoscale phenomena, which play a significant role in vertical nutrient supply. Simply increasing the resolution of the models may be an inefficient computational solution to this problem. An approach based on recent advances in adaptive mesh computational techniques may offer an alternative. Here the first steps in such an approach are described, using the example of a simple vertical column (quasi-1-D) ocean biogeochemical model. We present a novel method of simulating ocean biogeochemical behaviour on a vertically adaptive computational mesh, where the mesh changes in response to the biogeochemical and physical state of the system throughout the simulation. We show that the model reproduces the general physical and biological behaviour at three ocean stations (India, Papa and Bermuda) as compared to a high-resolution fixed mesh simulation and to observations. The use of an adaptive mesh does not increase the computational error, but reduces the number of mesh elements by a factor of 2–3. Unlike previous work the adaptivity metric used is flexible and we show that capturing the physical behaviour of the model is paramount to achieving a reasonable solution. Adding biological quantities to the adaptivity metric further refines the solution. We then show the potential of this method in two case studies where we change the adaptivity metric used to determine the varying mesh sizes in order to capture the dynamics of chlorophyll at Bermuda and sinking detritus at Papa. We therefore demonstrate that adaptive meshes may provide a suitable numerical technique for simulating seasonal or transient biogeochemical behaviour at high vertical resolution whilst minimising the number of elements in the mesh. More work is required to move this to fully 3-D simulations.

Physical forcing and physical/biochemical variability of the Mediterranean Sea: a review of unresolved issues and directions for future research - P. Malanotte-Rizzoli; V. Artale; G. L. Borzelli-Eusebi; S. Brenner; A. Crise; M. Gacic; N. Kress; S. Marullo; M. Ribera d'Alcalà; S. Sofianos; T. Tanhua; A. Theocharis; et al. (Ocean Science, 2014)
This paper is the outcome of a workshop held in Rome in November 2011 on the occasion of the 25th anniversary of the POEM (Physical Oceanography of the Eastern Mediterranean) program. In the workshop discussions, a number of unresolved issues were identified for the physical and biogeochemical properties of the Mediterranean Sea as a whole, i.e., comprising the Western and Eastern sub-basins. Over the successive two years, the related ideas were discussed among the group of scientists who participated in the workshop and who have contributed to the writing of this paper. Three major topics were identified, each of them being the object of a section divided into a number of different sub-sections, each addressing a specific physical, chemical or biological issue: 1. Assessment of basin-wide physical/biochemical properties, of their variability and interactions. 2. Relative importance of external forcing functions (wind stress, heat/moisture fluxes, forcing through straits) vs. internal variability. 3. Shelf/deep sea interactions and exchanges of physical/biogeochemical properties and how they affect the sub-basin circulation and property distribution. Furthermore, a number of unresolved scientific/methodological issues were also identified and are reported in each sub-section after a short discussion of the present knowledge. They represent the collegial consensus of the scientists contributing to the paper. Naturally, the unresolved issues presented here constitute the choice of the authors and therefore they may not be exhaustive and/or complete. The overall goal is to stimulate a broader interdisciplinary discussion among the scientists of the Mediterranean oceanographic community, leading to enhanced collaborative efforts and exciting future discoveries.

Simulated melt rates for the Totten and Dalton ice shelves - D. E. Gwyther; B. K. Galton-Fenzi; J. R. Hunter; J. L. Roberts (Ocean Science, 2014)
The Totten Glacier is rapidly losing mass. It has been suggested that this mass loss is driven by changes in oceanic forcing; however, the details of the ice–ocean interaction are unknown. Here we present results from an ice shelf–ocean model of the region that includes the Totten, Dalton and Moscow University ice shelves, based on the Regional Oceanic Modeling System for the period 1992–2007. Simulated area-averaged basal melt rates (net basal mass loss) for the Totten and Dalton ice shelves are 9.1 m ice yr−1 (44.5 Gt ice yr−1) and 10.1 m ice yr−1 (46.6 Gt ice yr−1), respectively. The melting of the ice shelves varies strongly on seasonal and interannual timescales. Basal melting (mass loss) from the Totten ice shelf spans a range of 5.7 m ice yr−1 (28 Gt ice yr−1) on interannual timescales and 3.4 m ice yr−1 (17 Gt ice yr−1) on seasonal timescales. This study links basal melt of the Totten and Dalton ice shelves to warm water intrusions across the continental shelf break and atmosphere–ocean heat exchange. Totten ice shelf melting is high when the nearby Dalton polynya interannual strength is below average, and vice versa. Melting of the Dalton ice shelf is primarily controlled by the strength of warm water intrusions across the Dalton rise and into the ice shelf cavity. During periods of strong westward coastal current flow, Dalton melt water flows directly under the Totten ice shelf further reducing melting. This is the first such modelling study of this region to provide a valuable framework for directing future observational and modelling efforts.

First laboratory study of air–sea gas exchange at hurricane wind speeds - K. E. Krall; B. Jähne (Ocean Science, 2014)
In a pilot study conducted in October and November 2011, air–sea gas transfer velocities of the two sparingly soluble trace gases hexafluorobenzene and 1,4-difluorobenzene were measured in the unique high-speed wind-wave tank at Kyoto University, Japan. This air–sea interaction facility is capable of producing hurricane strength wind speeds of up to u10 =67 m s−1. This constitutes the first lab study of gas transfer at such high wind speeds. The measured transfer velocities k600 spanned two orders of magnitude, lying between 11 cm h−1 and 1180 cm h−1 with the latter being the highest ever measured wind-induced gas transfer velocity. The measured gas transfer velocities are in agreement with the only available data set at hurricane wind speeds (McNeil and D'Asaro, 2007). The disproportionately large increase of the transfer velocities found at highest wind speeds indicates a new regime of air–sea gas transfer, which is characterized by strong wave breaking, enhanced turbulence and bubble cloud entrainment.

Meridional transport of salt in the global ocean from an eddy-resolving model - A. M. Treguier; J. Deshayes; J. Le Sommer; C. Lique; G. Madec; T. Penduff; J.-M. Molines; B. Barnier; R. Bourdalle-Badie; C. Talandier (Ocean Science, 2014)
The meridional transport of salt is computed in a global eddy-resolving numerical model (1/12° resolution) in order to improve our understanding of the ocean salinity budget. A methodology is proposed that allows a global analysis of the salinity balance in relation to surface water fluxes, without defining a "freshwater anomaly" based on an arbitrary reference salinity. The method consists of a decomposition of the meridional transport into (i) the transport by the time–longitude–depth mean velocity, (ii) time–mean velocity recirculations and (iii) transient eddy perturbations. Water is added (rainfall and rivers) or removed (evaporation) at the ocean surface at different latitudes, which creates convergences and divergences of mass transport with maximum and minimum values close to ±1 Sv. The resulting meridional velocity effects a net transport of salt at each latitude (±30 Sv PSU), which is balanced by the time–mean recirculations and by the net effect of eddy salinity–velocity correlations. This balance ensures that the total meridional transport of salt is close to zero, a necessary condition for maintaining a quasi-stationary salinity distribution. Our model confirms that the eddy salt transport cannot be neglected: it is comparable to the transport by the time–mean recirculation (up to 15 Sv PSU) at the poleward and equatorial boundaries of the subtropical gyres. Two different mechanisms are found: eddy contributions are localized in intense currents such as the Kuroshio at the poleward boundary of the subtropical gyres, while they are distributed across the basins at the equatorward boundaries. Closer to the Equator, salinity–velocity correlations are mainly due to the seasonal cycle and large-scale perturbations such as tropical instability waves.

The role of subpolar deep water formation and Nordic Seas overflows in simulated multidecadal variability of the Atlantic meridional overturning circulation - K. Lohmann; J. H. Jungclaus; D. Matei; J. Mignot; M. Menary; H. R. Langehaug; J. Ba; Y. Gao; O. H. Otterå; W. Park; S. Lorenz (Ocean Science, 2014)
We investigate the respective role of variations in subpolar deep water formation and Nordic Seas overflows for the decadal to multidecadal variability of the Atlantic meridional overturning circulation (AMOC). This is partly done by analysing long (order of 1000 years) control simulations with five coupled climate models. For all models, the maximum influence of variations in subpolar deep water formation is found at about 45° N, while the maximum influence of variations in Nordic Seas overflows is rather found at 55 to 60° N. Regarding the two overflow branches, the influence of variations in the Denmark Strait overflow is, for all models, substantially larger than that of variations in the overflow across the Iceland–Scotland Ridge. The latter might, however, be underestimated, as the models in general do not realistically simulate the flow path of the Iceland–Scotland overflow water south of the Iceland–Scotland Ridge. The influence of variations in subpolar deep water formation is, on multimodel average, larger than that of variations in the Denmark Strait overflow. This is true both at 45° N, where the maximum standard deviation of decadal to multidecadal AMOC variability is located for all but one model, and at the more classical latitude of 30° N. At 30° N, variations in subpolar deep water formation and Denmark Strait overflow explain, on multimodel average, about half and one-third respectively of the decadal to multidecadal AMOC variance. Apart from analysing multimodel control simulations, we have performed sensitivity experiments with one of the models, in which we suppress the variability of either subpolar deep water formation or Nordic Seas overflows. The sensitivity experiments indicate that variations in subpolar deep water formation and Nordic Seas overflows are not completely independent. We further conclude from these experiments that the decadal to multidecadal AMOC variability north of about 50° N is mainly related to variations in Nordic Seas overflows. At 45° N and south of this latitude, variations in both subpolar deep water formation and Nordic Seas overflows contribute to the AMOC variability, with neither of the processes being very dominant compared to the other.

Antarctic circumpolar transport and the southern mode: a model investigation of interannual to decadal timescales - C. W. Hughes; Joanne Williams; A. C. Coward; B. A. de Cuevas (Ocean Science, 2014)
It is well-established that, at periods shorter than a year, variations in Antarctic circumpolar transport are reflected in a barotropic mode, known as the southern mode, in which sea level and bottom pressure varies coherently around Antarctica. Here, we use two multidecadal ocean model runs to investigate the behaviour of the southern mode at timescales on which density changes become important, leading to a baroclinic component to the adjustment. We find that the concept of a southern mode in bottom pressure remains valid, and remains a direct measure of the circumpolar transport, with changes at the northern boundary playing only a small role even on decadal timescales. However, at periods longer than about 5 years, density changes start to play a role, leading to a surface intensification of the vertical profile of the transport. We also find that barotropic currents on the continental slope account for a significant fraction of the variability, and produce surface intensification in the meridional-integral flow. Circumpolar sea level and transport are related at all investigated timescales. However, the role of density variations results in a ratio of sea level change to transport which becomes larger at longer timescales. This means that any long-term transport monitoring strategy based on present measurement systems must involve multiplying the observed quantity by a factor which depends on frequency.

Intrinsic variability of the Antarctic Circumpolar Current system: low- and high-frequency fluctuations of the Argentine Basin flow - G. Sgubin; S. Pierini; H. A. Dijkstra (Ocean Science, 2014)
In this paper, the variability of the Antarctic Circumpolar Current system produced by purely intrinsic nonlinear oceanic mechanisms is studied through a sigma-coordinate ocean model, implemented in a large portion of the Southern Ocean at an eddy-permitting resolution under steady surface heat and momentum fluxes. The mean transport through the Drake Passage and the structure of the main Antarctic Circumpolar Current fronts are well reproduced by the model. Intrinsic variability is found to be particularly intense in the Subantarctic Front and in the Argentine Basin, on which further analysis is focused. The low-frequency variability at interannual timescales is related to bimodal behavior of the Zapiola Anticyclone, with transitions between a strong and collapsed anticyclonic circulation in substantial agreement with altimeter observations. Variability on smaller timescales shows clear evidence of topographic Rossby-wave propagation along the eastern and southern flanks of the Zapiola Rise and of mesoscale eddies, also in agreement with altimeter observations. The analysis of the relationship between the low- and high-frequency variability suggests possible mechanisms of mutual interaction.

Secchi depth in the Oslofjord–Skagerrak area: theory, experiments and relationships to other quantities - E. Aas; J. Høkedal; K. Sørensen (Ocean Science, 2014)
The Secchi depth and its relationships to other properties of the sea water in the Oslofjord–Skagerrak area have been investigated. White and black disks of different sizes have been applied, and the Secchi depth has been observed with the naked eye, through colour filters and with a water telescope. Spectral luminances and illuminances have been calculated from recordings of radiance and irradiance, and attenuation coefficients have been determined. A theoretical expression for the Secchi depth based on luminances has been tested against field observations, and it is found that the field results for the product of Secchi depth and attenuation coefficients are on average only 4% less than the predicted value for the white disk. For the Secchi depths observed through colour filters or for the black disk, the average field results are more than 30% smaller than the theoretical estimates. The reduction in the disk diameter from 30 to 10 cm should theoretically reduce the Secchi depths by 13–22%, while the field observations show an average reduction of 10–20%. Similarly we find from theory that the removal of sun glitter should increase the Secchi depth by 12%, while the observed increase is 14% on average for the white disk. Our overall conclusion is that the theoretical expression works well for the white disk, but less so for the colour filter observations and the black disk. Statistical relationships between Secchi depths and attenuation coefficients have been determined, and it is found that the root-mean-square errors relative to the mean value are smaller for the beam attenuation coefficients (12–24%, white disk) than for the vertical attenuation coefficients (16–65%, white disk). The depth of the 1% level of surface quanta irradiance (PAR) can be estimated with a relative root-mean-square error of 23% from observations of the white Secchi depth. Similar estimates of chlorophyll a and total suspended material will have rms errors in the range 40–90%. Our conclusion becomes that the Secchi depth observation is a very useful tool for checking the value and order of magnitude of other related quantities in the Oslofjord–Skagerrak area.

Turbulent dispersion properties from a model simulation of the western Mediterranean Sea - H. Nefzi; D. Elhmaidi; X. Carton (Ocean Science, 2014)
Using a high-resolution primitive equation model of the western Mediterranean Sea, we analyzed the dispersion properties of a set of homogeneously distributed, passive particle pairs. These particles were initially separated by different distances D0 (D0 = 5.55, 11.1 and 16.65 km), and were seeded in the model at initial depths of 44 and 500 m. This realistic ocean model, which reproduces the main features of the regional circulation, puts into evidence the three well-known regimes of relative dispersion. The first regime due to the chaotic advection at small scales lasts only a few days (3 days at 44 m depth, a duration comparable with the integral timescale), and the relative dispersion is then exponential. In the second regime, extending from 3 to 20 days, the relative dispersion has a power law tα where α tends to 3 as D0 becomes small. In the third regime, a linear growth of the relative dispersion is observed starting from the twentieth day. For the relative diffusivity, the D2 growth is followed by the Richardson regime D4/3. At large scales, where particle velocities are decorrelated, the relative diffusivity is constant. At 500 m depth, the integral timescale increases (> 4 days) and the intermediate regime becomes narrower than that at 44 m depth due to the weaker effect of vortices (this effect decreases with depth). The turbulent properties become less intermittent and more homogeneous and the Richardson law takes place.

Geostrophic currents and kinetic energies in the Black Sea estimated from merged drifter and satellite altimetry data - M. Menna; P.-M. Poulain (Ocean Science, 2014)
Drifter measurements and satellite altimetry data are merged to reconstruct the surface geostrophic circulation of the Black Sea in the period 1999–2009. This combined data set is used to estimate pseudo-Eulerian velocity statistics for different time periods. Seasonal and interannual variability of currents and kinetic energy fields are described with particular attention to the mesoscale and sub-basin coastal eddies. The mean currents are generally stronger in winter and enhanced speeds are observed in the period 2002–2006. The most intense activity of sub-basin Batumi Eddy occurs in summer with greater speeds and dimensions in 2006 and 2008. The sub-basin Sevastopol Eddy is generated in spring from a meander of the Rim Current. Mesoscale eddies located along the Anatolia, Caucasus and Crimea coasts are permanent, quasi-permanent or intermittent features and can interact and merge with each other, showing high values of kinetic energy.

Halocline water modification and along-slope advection at the Laptev Sea continental margin - D. Bauch; S. Torres-Valdes; I. Polyakov; A. Novikhin; I. Dmitrenko; J. McKay; A. Mix (Ocean Science, 2014)
A general pattern in water mass distribution and potential shelf–basin exchange is revealed at the Laptev Sea continental slope based on hydrochemical and stable oxygen isotope data from the summers 2005–2009. Despite considerable interannual variations, a frontal system can be inferred between shelf, continental slope and central Eurasian Basin waters in the upper 100 m of the water column along the continental slope. Net sea-ice melt is consistently found at the continental slope. However, the sea-ice meltwater signal is independent from the local retreat of the ice cover and appears to be advected from upwind locations. In addition to the along-slope frontal system at the continental shelf break, a strong gradient is identified on the Laptev Sea shelf between 122° E and 126° E with an eastward increase of riverine and sea-ice related brine water contents. These waters cross the shelf break at ~ 140° E and feed the low-salinity halocline water (LSHW, salinity S

Decadal variability of heat content in the South China Sea inferred from observation data and an ocean data assimilation product - Wei Song; Jian Lan; Qinyan Liu; Dandan Sui; Lili Zeng; Dongxiao Wang (Ocean Science, 2014)
Using an observation data set of temperature and the Simple Ocean Data Assimilation (SODA), the decadal variability of upper ocean heat content (0–400 m; hereafter, OHC) in the South China Sea (SCS) was investigated for the period from 1958 to 2007. Decadal variability was identified as the dominant mode of upper OHC besides the seasonal cycle. According to decreasing or increasing OHC, four periods were chosen to discuss detailed processes behind OHC variability in the SCS; the four periods are 1958–1968, 1969–1981, 1982–1992, and 1993–2003. Results show that advection was the major factor for decreasing (increasing) OHC during 1958–1968 (1968–1981). During 1982–1992 and 1993–2003, the net surface heat flux was the main contributor to the variability of OHC besides the advection. The OHC, advection and net surface heat flux had significant positive trends during 1992–2003. The spatial characteristics of OHC variability and heat budget in the Luzon Strait, west of Luzon Island, and in the Xisha warm eddy region are also discussed.

The instability of diffusive convection and its implication for the thermohaline staircases in the deep Arctic Ocean - S.-Q. Zhou; L. Qu; Y.-Z. Lu; X.-L. Song (Ocean Science, 2014)
In the present study, the classical description of diffusive convection is updated to interpret the instability of diffusive interfaces and the dynamical evolution of the bottom layer in the deep Arctic Ocean. In the new consideration of convective instability, both the background salinity stratification and rotation are involved. The critical Rayleigh number of diffusive convection is found to vary from 103 to 1011 in the deep Arctic Ocean as well as in other oceans and lakes. In such a wide range of conditions, the interface-induced thermal Rayleigh number is shown to be consistent with the critical Rayleigh number of diffusive convection. In most regions, background salinity stratification is found to be the main hindrance to the occurrence of convecting layers. With the new parameterization, it is predicted that the maximum thickness of the bottom layer is 1051 m in the deep Arctic Ocean, which is close to the observed value of 929 m. The evolution time of the bottom layer is predicted to be ~ 100 yr, which is on the same order as that based on 14C isolation age estimation.

Fate of colloids during estuarine mixing in the Arctic - O. S. Pokrovsky; L. S. Shirokova; J. Viers; V. V. Gordeev; V. P. Shevchenko; A. V. Chupakov; T. Y. Vorobieva; F. Candaudap; C. Causserand; A. Lanzanova; C. Zouiten (Ocean Science, 2014)
The estuarine behavior of organic carbon (OC) and trace elements (TE) was studied for the largest European sub-Arctic river, which is the Severnaya Dvina; this river has a deltaic estuary covered in ice during several hydrological seasons: summer (July 2010, 2012) and winter (March 2009) baseflow, and the November–December 2011 ice-free period. Colloidal forms of OC and TE were assessed for three pore size cutoffs (1, 10, and 50 kDa) using an in situ dialysis procedure. Conventionally dissolved (

Temporal variations of zooplankton biomass in the Ligurian Sea inferred from long time series of ADCP data - R. Bozzano; E. Fanelli; S. Pensieri; P. Picco; M. E. Schiano (Ocean Science, 2014)
Three years of 300 kHz acoustic doppler current profiler data collected in the central Ligurian Sea are analysed to investigate the variability of the zooplankton biomass and the diel vertical migration in the upper thermocline. After a pre-processing phase aimed at avoiding the slant range attenuation, hourly volume backscattering strength time series are obtained. Despite the lack of concurrent net samples collection, different migration patterns are identified and their temporal variability examined by means of time–frequency analysis. The effect of changes in the environmental condition is also investigated. The highest zooplankton biomasses are observed in April–May just after the peak of surface primary production in March–April. The main migration pattern found here points to a "nocturnal" migration, with zooplankton organisms occurring deeper in the water column during the day and shallower at night. Also, twilight migration is highlighted during this study. The largest migrations are recorded in November–December, corresponding to lowest backscattering strength values and they are likely attributable to larger and more active organisms (i.e. euphausiids and mesopelagic fish). The results suggest further applications of the available historical acoustic doppler current profiler time series.

The CO 2 system in the Mediterranean Sea: a basin wide perspective - M. Álvarez; H. Sanleón-Bartolomé; T. Tanhua; L. Mintrop; A. Luchetta; C. Cantoni; K. Schroeder; G. Civitarese (Ocean Science, 2014)
The Mediterranean Sea (MedSea) is considered a "laboratory basin" being an ocean in miniature, suffering dramatic changes in its oceanographic and biogeochemical conditions derived from natural and anthropogenic forces. Moreover, the MedSea is prone to absorb and store anthropogenic carbon due to the particular CO2 chemistry and the active overturning circulation. Despite this, water column CO2 measurements covering the whole basin are scarce. This work aims to be a base-line for future studies about the CO2 system space-time variability in the MedSea combining historic and modern CO2 cruises in the whole area. Here we provide an extensive vertical and longitudinal description of the CO2 system variables (total alkalinity – TA, dissolved inorganic carbon – DIC and pH) along an East-West transect and across the Sardinia-Sicily passage in the MedSea from two oceanographic cruises conducted in 2011 measuring CO2 variables in a coordinated fashion, the RV Meteor M84/3 and the RV Urania EuroFleets 11, respectively. In this sense, we provide full-depth and length CO2 distributions across the MedSea, and property-property plots showing in each sub-basin post-Eastern Mediterranean Transient (EMT) situation with regard to TA, DIC and pH. The over-determined CO2 system in 2011 allowed performing the first internal consistency analysis for the particularly warm, high salinity and alkalinity MedSea waters. The CO2 constants by Mehrbach et al. (1973) refitted by Dickson and Millero (1987) are recommended. The sensitivity of the CO2 system to the atmospheric CO2 increase, DIC and/or TA changes is evaluated by means of the Revelle and buffer factors.

Tidal variability of the motion in the Strait of Otranto - L. Ursella; V. Kovačević; M. Gačić (Ocean Science, 2014)
Various current data, collected in the Strait of Otranto during the period 1994–2007, have been analysed with the aim of describing the characteristics of the tidal motions and their contribution to the total flow variance. The principal tidal constituents in the area were the semi-diurnal (M2) and the diurnal (K1), with the latter one predominant. The total flow was, in general, more energetic along the flanks than in the middle of the strait. Specifically, it was most energetic over the western shelf and in the upper layer along the eastern flank. In spite of the generally low velocities (a few cm s−1) of the principal tidal constituents, the tidal variance has a pattern similar to that of the total flow variance, that is, it was large over the western shelf and low in the middle. The proportion of non-tidal (comprising the inertial and sub-inertial low-frequency bands) to tidal flow variances was quite variable in both time and space. The low-frequency motions dominated over the tidal and inertial ones in the eastern portion of the strait during the major part of the year, particularly in the upper and intermediate layers. In the deep, near-bottom layer the variance was evenly distributed between the low frequency, diurnal and semi-diurnal bands. An exception was observed near the western shelf break during the summer season when the contribution of the tidal signal to the total variance reached 77%. This high contribution was mainly due to the intensification of the diurnal signal at that location at both upper and bottom current records (velocities of about 10 cm s−1). Local wind and sea level data were analysed and compared with the flow to find the possible origin of this diurnal intensification. Having excluded the sea-breeze impact on the intensification of the diurnal tidal signal, the most likely cause remains the generation of the topographically trapped internal waves and the diurnal resonance in the tidal response. These waves were sometimes generated by the barotropic tidal signal in the presence of summer stratification and the strong bottom slope. This phenomenon may stimulate diapycnal mixing during the stratified season and enhance ventilation of the near-bottom layers.

Enhancing the accuracy of automatic eddy detection and the capability of recognizing the multi-core structures from maps of sea level anomaly - J. Yi; Y. Du; Z. He; C. Zhou (Ocean Science, 2014)
Automated methods are important for automatically detecting mesoscale eddies in large volumes of altimeter data. While many algorithms have been proposed in the past, this paper presents a new method, called hybrid detection (HD), to enhance the eddy detection accuracy and the capability of recognizing eddy multi-core structures from maps of sea level anomaly (SLA). The HD method has integrated the criteria of the Okubo–Weiss (OW) method and the sea surface height-based (SSH-based) method, two commonly used eddy detection algorithms. Evaluation of the detection accuracy shows that the successful detection rate of HD is ~ 96.6% and the excessive detection rate is ~ 14.2%, which outperforms the OW and those methods using SLA extrema to identify eddies. The capability of recognizing multi-core structures and its significance in tracking eddy splitting or merging events have been illustrated by comparing with the detection results of different algorithms and observations in previous literature.

Observed decline of the Atlantic meridional overturning circulation 2004–2012 - D. A. Smeed; G. D. McCarthy; S. A. Cunningham; E. Frajka-Williams; D. Rayner; W. E. Johns; C. S. Meinen; M. O. Baringer; B. I. Moat; A. Duchez; H. L. Bryden (Ocean Science, 2014)
The Atlantic meridional overturning circulation (AMOC) has been observed continuously at 26° N since April 2004. The AMOC and its component parts are monitored by combining a transatlantic array of moored instruments with submarine-cable-based measurements of the Gulf Stream and satellite derived Ekman transport. The time series has recently been extended to October 2012 and the results show a downward trend since 2004. From April 2008 to March 2012, the AMOC was an average of 2.7 Sv (1 Sv = 106 m3 s−1) weaker than in the first four years of observation (95% confidence that the reduction is 0.3 Sv or more). Ekman transport reduced by about 0.2 Sv and the Gulf Stream by 0.5 Sv but most of the change (2.0 Sv) is due to the mid-ocean geostrophic flow. The change of the mid-ocean geostrophic flow represents a strengthening of the southward flow above the thermocline. The increased southward flow of warm waters is balanced by a decrease in the southward flow of lower North Atlantic deep water below 3000 m. The transport of lower North Atlantic deep water slowed by 7% per year (95% confidence that the rate of slowing is greater than 2.5% per year).

A parameter model of gas exchange for the seasonal sea ice zone - B. Loose; W. R. McGillis; D. Perovich; C. J. Zappa; P. Schlosser (Ocean Science, 2014)
Carbon budgets for the polar oceans require better constraint on air–sea gas exchange in the sea ice zone (SIZ). Here, we utilize advances in the theory of turbulence, mixing and air–sea flux in the ice–ocean boundary layer (IOBL) to formulate a simple model for gas exchange when the surface ocean is partially covered by sea ice. The gas transfer velocity (k) is related to shear-driven and convection-driven turbulence in the aqueous mass boundary layer, and to the mean-squared wave slope at the air–sea interface. We use the model to estimate k along the drift track of ice-tethered profilers (ITPs) in the Arctic. Individual estimates of daily-averaged k from ITP drifts ranged between 1.1 and 22 m d−1, and the fraction of open water (f) ranged from 0 to 0.83. Converted to area-weighted effective transfer velocities (keff), the minimum value of keff was 10−55 m d−1 near f = 0 with values exceeding keff = 5 m d−1 at f = 0.4. The model indicates that effects from shear and convection in the sea ice zone contribute an additional 40% to the magnitude of keff, beyond what would be predicted from an estimate of keff based solely upon a wind speed parameterization. Although the ultimate scaling relationship for gas exchange in the sea ice zone will require validation in laboratory and field studies, the basic parameter model described here demonstrates that it is feasible to formulate estimates of k based upon properties of the IOBL using data sources that presently exist.

Changes in ventilation of the Mediterranean Sea during the past 25 year - A. Schneider; T. Tanhua; W. Roether; R. Steinfeldt (Ocean Science, 2014)
Significant changes in the overturning circulation of the Mediterranean Sea has been observed during the last few decades, the most prominent phenomena being the Eastern Mediterranean Transient (EMT) in the early 1990s and the Western Mediterranean Transition (WMT) during the mid-2000s. During both of these events unusually large amounts of deep water were formed, and in the case of the EMT, the deep water formation area shifted from the Adriatic to the Aegean Sea. Here we synthesize a unique collection of transient tracer (CFC-12, SF6 and tritium) data from nine cruises conducted between 1987 and 2011 and use these data to determine temporal variability of Mediterranean ventilation. We also discuss biases and technical problems with transient tracer-based ages arising from their different input histories over time; particularly in the case of time-dependent ventilation. We observe a period of low ventilation in the deep eastern (Levantine) basin after it was ventilated by the EMT so that the age of the deep water is increasing with time. In the Ionian Sea, on the other hand, we see evidence of increased ventilation after year 2001, indicating the restarted deep water formation in the Adriatic Sea. This is also reflected in the increasing age of the Cretan Sea deep water and decreasing age of Adriatic Sea deep water since the end of the 1980s. In the western Mediterranean deep basin we see the massive input of recently ventilated waters during the WMT. This signal is not yet apparent in the Tyrrhenian Sea, where the ventilation seems to be fairly constant since the EMT. Also the western Alboran Sea does not show any temporal trends in ventilation.

A practical scheme to introduce explicit tidal forcing into an OGCM - K. Sakamoto; H. Tsujino; H. Nakano; M. Hirabara; G. Yamanaka (Ocean Science, 2013)
A practical scheme is proposed to explicitly introduce tides into ocean general circulation models (OGCM). In this scheme, barotropic linear response to the tidal forcing is calculated by the time differential equations modified for ocean tides, instead of the original barotropic equations of an OGCM. This allows for the usage of various parameterizations specified for tides, such as the self-attraction/loading (SAL) effect and energy dissipation due to internal tides, without unintentional violation of the original dynamical balances in an OGCM. Meanwhile, secondary nonlinear effects of tides, e.g., excitation of internal tides and advection by tidal currents, are fully represented within the framework of the original OGCM equations. That is, this scheme drives the OGCM by the barotropic linear tidal currents which are predicted progressively by a tuned tide model, instead of the equilibrium tide potential, without large additional numerical costs. We incorporated this scheme into Meteorological Research Institute Community Ocean Model and executed test experiments with a low-resolution global model. The results showed that the model can simulate both the non-tidal circulations and the tidal motion simultaneously. Owing to the usage of tidal parameterizations such as a SAL term, a root-mean-squared error in the tidal heights is found to be as small as 10.0 cm, which is comparable to that of elaborately tuned tide models. In addition, analysis of the speed and energy of the barotropic tidal currents is found to be consistent with that of past tide studies. The model also showed active excitement of internal tides and tidal mixing. In the future, the impacts of internal tides and tidal mixing should be examined using a model with a finer resolution, since explicit and precise introduction of tides into an OGCM is a significant step toward the improvement of ocean models.

A new method for continuous measurements of oceanic and atmospheric N 2 O, CO and CO 2 : performance of off-axis integrated cavity output spectroscopy (OA-ICOS) coupled to non-dispersive infrared detection (NDIR) - D. L. Arévalo-Martínez; M. Beyer; M. Krumbholz; I. Piller; A. Kock; T. Steinhoff; A. Körtzinger; H. W. Bange (Ocean Science, 2013)
A new system for continuous, highly resolved oceanic and atmospheric measurements of N2O, CO and CO2 is described. The system is based upon off-axis integrated cavity output spectroscopy (OA-ICOS) and a non-dispersive infrared analyzer (NDIR), both coupled to a Weiss-type equilibrator. Performance of the combined setup was evaluated by testing its precision, accuracy, long-term stability, linearity and response time. Furthermore, the setup was tested during two oceanographic campaigns in the equatorial Atlantic Ocean in order to explore its potential for autonomous deployment onboard voluntary observing ships (VOS). Improved equilibrator response times for N2O (2.5 min) and CO (45 min) were achieved in comparison to response times from similar chamber designs used by previous studies. High stability of the OA-ICOS analyzer was demonstrated by low optimal integration times of 2 and 4 min for N2O and CO respectively, as well as detection limits of

Atlantic meridional ocean heat transport at 26° N: impact on subtropical ocean heat content variability - M. Sonnewald; J. J.-M. Hirschi; R. Marsh; E. L. McDonagh; B. A. King (Ocean Science, 2013)
Local climate is significantly affected by changes in the oceanic heat content on a range of timescales. This variability is driven by heat fluxes from both the atmosphere and the ocean. In the Atlantic the meridional overturning circulation is the main contributor to the oceanic meridional heat transport for latitudes south of about 50° N. The RAPID project has been successfully monitoring the Atlantic meridional overturning at 26° N since 2004. This study demonstrates how these data can be used to estimate the variability of the basin-wide ocean heat content in the upper 800 m between 26° and 36° N. Traditionally the atmosphere is seen to dominate the ocean heat content variability. However, previous studies have looked at smaller areas in the Gulf Stream region, finding that the ocean dominates deseasoned fluctuations of ocean heat content, while studies of the whole North Atlantic region suggest that the atmosphere may be dominant. In our study we use a box model to investigate fluctuations of the ocean heat content in the subtropical North Atlantic between 26° and 36° N. The box model approach is validated using 19 yr of high-resolution general circulation model (GCM) data. We find that in both the GCM- and RAPID-based data the ocean heat transport dominates the deseasoned heat content variability, while the atmosphere's impact on the ocean heat content evolution stabilizes after 6 months. We demonstrate that the utility of the RAPID data goes beyond monitoring the overturning circulation at 26° N, and that it can be used to better understand the causes of ocean heat content variability in the North Atlantic. We illustrate this for a recent decrease in ocean heat content which was observed in the North Atlantic in 2009 and 2010. Our results suggest that most of this ocean heat content reduction can be explained by a reduction of the meridional ocean heat transport during this period.

Time and space variability of freshwater content, heat content and seasonal ice melt in the Arctic Ocean from 1991 to 2011 - M. Korhonen; B. Rudels; M. Marnela; A. Wisotzki; J. Zhao (Ocean Science, 2013)
Changes in the hydrography of the Arctic Ocean have recently been reported. The upper ocean has been freshening and pulses of warm Atlantic Water have been observed to spread into the Arctic Ocean. Although these changes have been intensively studied, salinity and temperature variations have less frequently been considered together. Here hydrographic observations, obtained by icebreaker expeditions conducted between 1991 and 2011, are analyzed and discussed. Five different water masses in the upper 1000 m of the water column are examined in five sub-basins of the Arctic Ocean. This allows for studying the variations of the distributions of the freshwater and heat contents in the Arctic Ocean not only in time but also laterally and vertically. In addition, the seasonal ice melt contribution is separated from the permanent, winter, freshwater content of the Polar Mixed Layer. Because the positions of the icebreaker stations vary between the years, the icebreaker observations are at each specific point in space and time compared with the Polar Science Center Hydrographic Climatology to separate the effects of space and time variability on the observations. The hydrographic melt water estimate is discussed and compared with the potential ice melt induced by atmospheric heat input estimated from the ERA–Interim and NCEP/NCAR reanalyses. After a period of increased salinity in the upper ocean during the 1990s, both the Polar Mixed Layer and the upper halocline have been freshening. The increase in freshwater content in the Polar Mixed Layer is primarily driven by a decrease in salinity, not by changes in Polar Mixed Layer depth, whereas the freshwater is accumulating in the upper halocline mainly through the increasing thickness of the halocline. This is especially evident in the Northern Canada Basin, where the most substantial freshening is observed. The warming, and to some extent also the increase in salinity, of the Atlantic Water during the early 1990s extended from the Nansen Basin into the Amundsen and Makarov basins, while the warm and saline inflows occurring during the 2000s appear to be confined to the Nansen Basin, suggesting that the warm and saline inflow through Fram Strait largely recirculates in the Nansen Basin.

Microstructure measurements and estimates of entrainment in the Denmark Strait overflow plume - V. Paka; V. Zhurbas; B. Rudels; D. Quadfasel; A. Korzh; D. Delisi (Ocean Science, 2013)
To examine processes controlling the entrainment of ambient water into the Denmark Strait overflow (DSO) plume/gravity current, measurements of turbulent dissipation rate were carried out by a quasi-free-falling (tethered) microstructure profiler (MSP). The MSP was specifically designed to collect data on dissipation-scale turbulence and fine thermohaline stratification in an ocean layer located as deep as 3500 m. The task was to perform microstructure measurements in the DSO plume in the lower 300 m depth interval including the bottom mixed layer and the interfacial layer below the non-turbulent ambient water. The MSP was attached to a Rosette water sampler rack equipped with a SeaBird CTDO and an RD Instruments lowered acoustic Doppler current profiler (LADCP). At a chosen depth, the MSP was remotely released from the rack to perform measurements in a quasi-free-falling mode. Using the measured vertical profiles of dissipation, the entrainment rate as well as the bottom and interfacial stresses in the DSO plume were estimated at a location 200 km downstream of the sill at depths up to 1771 m. Dissipation-derived estimates of entrainment were found to be much smaller than bulk estimates of entrainment calculated from the downstream change of the mean properties in the plume, suggesting the lateral stirring due to mesoscale eddies rather than diapycnal mixing as the main contributor to entrainment. Dissipation-derived bottom stress estimates are argued to be roughly one third the magnitude of those derived from log velocity profiles. In the interfacial layer, the Ozmidov scale calculated from turbulence dissipation rate and buoyancy frequency was found to be linearly proportional to the overturning scale extracted from conventional CTD data (the Thorpe scale), with a proportionality constant of 0.76, and a correlation coefficient of 0.77.

An optical model for deriving the spectral particulate backscattering coefficients in oceanic waters - S. P. Tiwari; P. Shanmugam (Ocean Science, 2013)
An optical model is developed based on the diffuse attenuation coefficient (Kd) to estimate particulate backscattering coefficients bbp(λ) in oceanic waters. A large in situ data set is used to establish robust relationships between bbp(530) and bbp(555) and Kd(490) using an efficient nonlinear least-square method which uses the trust region algorithm with Bisquare weights scheme to adjust the coefficients. These relationships are obtained with good correlation coefficients (R2 = 0.786 and 0.790), low root mean square error (RMSE = 0.00076 and 0.00072) and 95% confidence bounds. The new model is tested with three independent data sets: the NOMAD SeaWiFS Match ups, OOXIX IOP algorithm workshop evaluation data set (Version 2.0w APLHA), and IOCCG simulated data set. Results show that the new model makes good retrievals of bbp at all key wavelengths (from 412–683 nm), with statistically significant improvements over other inversion models. Thus, the new model has the potential to improve our present knowledge of particulate matter and their optical variability in oceanic waters.

Near-surface diurnal warming simulations: validation with high resolution profile measurements - B. Scanlon; G. A. Wick; B. Ward (Ocean Science, 2013)
Sea surface temperature (SST) is an important property for governing the exchange of energy between the ocean and the atmosphere. Common in situ methods of measuring SST often require a cool-skin and warm-layer adjustment in the presence of diurnal warming effects. A critical requirement for an ocean submodel is that it can simulate the change in SST over diurnal, seasonal and annual cycles. In this paper we use high-resolution near-surface profiles of SST to validate simulated near-surface temperature profiles from a modified version of the Kantha and Clayson 1-D mixed-layer model. Additional model enhancements such as the incorporation of a more recent parameterization of turbulence generated by wave breaking and a recent solar absorption model are also validated. The model simulations show a strong variability in highly stratified conditions, with different models providing the best results depending on the specific criteria and conditions. In general, the models with enhanced wave breaking effects provided underestimated temperature profiles while the more coarse baseline and blended approaches produced the most accurate SST estimates.

Frontal structures in the West Spitsbergen Current margins - W. Walczowski (Ocean Science, 2013)
The structures of the hydrographic fronts separating the Atlantic-origin waters from ambient waters in the northern Nordic Seas are discussed. Flows of the western and eastern branches of the West Spitsbergen Current create the Atlantic domain borders and maintain these fronts. This work is based on previous research and on investigations carried out in the project DAMOCLES (Developing Arctic Modelling and Observational Capabilities for Long-term Environmental Studies). Most of the observational data were collected during the R/V Oceania cruises. The main focus of the paper is the western border of the Atlantic domain – the Arctic Front, alongfrontal and transfrontal transports, and the front instability and variability. The alongfrontal baroclinic jet streams were described as a significant source of the Atlantic Water and heat in the Nordic Seas. The baroclinic instability and advection of baroclinic eddies which occurs due to this instability were found to be the main transfrontal transport processes. Most of the Atlantic Water transported by the western branch recirculates west and southward. The eastern branch of the West Spitsbergen Current provides most of the Atlantic Water entering the Arctic Ocean. Both processes are very important for the Arctic and global thermohaline circulation.

The circulation of Icelandic waters – a modelling study - K. Logemann; J. Ólafsson; Á. Snorrason; H. Valdimarsson; G. Marteinsdóttir (Ocean Science, 2013)
The three-dimensional flow, temperature and salinity fields of the North Atlantic, including the Arctic Ocean, covering the time period 1992 to 2006 are simulated with the numerical ocean model CODE. The simulation reveals several new insights and previously unknown structures which help us to clarify open questions on the regional oceanography of Icelandic waters. These relate to the structure and geographical distribution of the coastal current, the primary forcing of the North Icelandic Irminger Current (NIIC) and the path of the Atlantic Water south-east of Iceland. The model's adaptively refined computational mesh has a maximum resolution of 1 km horizontal and 2.5 m vertical in Icelandic waters. CTD profiles from this region and the river discharge of 46 Icelandic watersheds, computed by the hydrological model WaSiM, are assimilated into the simulation. The model realistically reproduces the established elements of the circulation around Iceland. However, analysis of the simulated mean flow field also provides further insights. It suggests a distinct freshwater-induced coastal current that only exists along the south-west and west coasts, which is accompanied by a counter-directed undercurrent. The simulated transport of Atlantic Water over the Icelandic shelf takes place in a symmetrical system of two currents, with the established NIIC over the north-western and northern shelf, and a hitherto unnamed current over the southern and south-eastern shelf, which is simulated to be an upstream precursor of the Faroe Current (FC). Both currents are driven by barotropic pressure gradients induced by a sea level slope across the Greenland–Scotland Ridge. The recently discovered North Icelandic Jet (NIJ) also features in the model predictions and is found to be forced by the baroclinic pressure field of the Arctic Front, to originate east of the Kolbeinsey Ridge and to have a volume transport of around 1.5 Sv within northern Denmark Strait. The simulated multi-annual mean Atlantic Water transport of the NIIC increased by 85% during 1992 to 2006, whereas the corresponding NIJ transport decreased by 27%. Based on our model results we propose a new and further differentiated circulation scheme of Icelandic waters whose details may inspire future observational oceanography studies.

Estuarine circulation reversals and related rapid changes in winter near-bottom oxygen conditions in the Gulf of Finland, Baltic Sea - T. Liblik; J. Laanemets; U. Raudsepp; J. Elken; I. Suhhova (Ocean Science, 2013)
The reversal of estuarine circulation caused by southwesterly wind forcing may lead to vanishing of stratification and subsequently to oxygenation of deep layers during the winter in the Gulf of Finland. Six conductivity, temperature, depth (CTD)+oxygen transects (130 km long, 10 stations) were conducted along the thalweg from the western boundary to the central gulf (21 December 2011–8 May 2012). Two bottom-mounted ADCP were installed, one near the western border and the second in the central gulf. A CTD with a dissolved oxygen sensor was deployed close to the western ADCP. Periods of typical estuarine circulation were characterized by strong stratification, high salinity, hypoxic conditions and inflow to the gulf in the near-bottom layer. Two circulation reversals were observed: one in December–January and one in February. The first reversal event was well developed; it caused the disappearance of the stratification and an increase in the oxygen concentration from hypoxic values to 270 μmol L−1 (to 6 mL L−1) throughout the water column along the thalweg and lasted approximately 1.5 months. Shifts from estuarine circulation to reversed circulation and vice versa were both associated with strong longitudinal (east–west) gulf currents (up to 40 cm s−1) in the deep layer. The change from oxygenated to hypoxic conditions in the western near-entrance area of the gulf occurred very rapidly, within less than a day, due to the intrusion of the hypoxic salt wedge from the NE Baltic Proper. In the eastern part of the gulf, good oxygen conditions caused by reversals remained for a few months.

From satellite altimetry to Argo and operational oceanography: three revolutions in oceanography - P. Y. Le Traon (Ocean Science, 2013)
The launch of the French/US mission Topex/Poseidon (T/P) (CNES/NASA) in August 1992 was the start of a revolution in oceanography. For the first time, a very precise altimeter system optimized for large-scale sea level and ocean circulation observations was flying. T/P alone could not observe the mesoscale circulation. In the 1990s, the ESA satellites ERS-1/2 were flying simultaneously with T/P. Together with my CLS colleagues, we demonstrated that we could use T/P as a reference mission for ERS-1/2 and bring the ERS-1/2 data to an accuracy level comparable to T/P. Near-real-time high-resolution global sea level anomaly maps were then derived. These maps have been operationally produced as part of the SSALTO/DUACS system for the last 15 yr. They are now widely used by the oceanographic community and have contributed to a much better understanding and recognition of the role and importance of mesoscale dynamics. Altimetry needs to be complemented with global in situ observations. At the end of the 90s, a major international initiative was launched to develop Argo, the global array of profiling floats. This has been an outstanding success. Argo floats now provide the most important in situ observations to monitor and understand the role of the ocean on the earth climate and for operational oceanography. This is a second revolution in oceanography. The unique capability of satellite altimetry to observe the global ocean in near-real-time at high resolution and the development of Argo were essential for the development of global operational oceanography, the third revolution in oceanography. The Global Ocean Data Assimilation Experiment (GODAE) was instrumental in the development of the required capabilities. This paper provides an historical perspective on the development of these three revolutions in oceanography which are very much interlinked. This is not an exhaustive review and I will mainly focus on the contributions we made together with many colleagues and friends.

Tidally induced lateral dispersion of the Storfjorden overflow plume - F. Wobus; G. I. Shapiro; J. M. Huthnance; M. A. M. Maqueda; Y. Aksenov (Ocean Science, 2013)
We investigate the flow of brine-enriched shelf water from Storfjorden (Svalbard) into Fram Strait and onto the western Svalbard Shelf using a regional set-up of NEMO-SHELF, a 3-D numerical ocean circulation model. The model is set up with realistic bathymetry, atmospheric forcing, open boundary conditions and tides. The model has 3 km horizontal resolution and 50 vertical levels in the sh-coordinate system which is specially designed to resolve bottom boundary layer processes. In a series of modelling experiments we focus on the influence of tides on the propagation of the dense water plume by comparing results from tidal and non-tidal model runs. Comparisons of non-tidal to tidal simulations reveal a hotspot of tidally induced horizontal diffusion leading to the lateral dispersion of the plume at the southernmost headland of Spitsbergen which is in close proximity to the plume path. As a result the lighter fractions in the diluted upper layer of the plume are drawn into the shallow coastal current that carries Storfjorden water onto the western Svalbard Shelf, while the dense bottom layer continues to sink down the slope. This bifurcation of the plume into a diluted shelf branch and a dense downslope branch is enhanced by tidally induced shear dispersion at the headland. Tidal effects at the headland are shown to cause a net reduction in the downslope flux of Storfjorden water into the deep Fram Strait. This finding contrasts previous results from observations of a dense plume on a different shelf without abrupt topography.

Optimal adjustment of the atmospheric forcing parameters of ocean models using sea surface temperature data assimilation - M. Meinvielle; J.-M. Brankart; P. Brasseur; B. Barnier; R. Dussin; J. Verron (Ocean Science, 2013)
In ocean general circulation models, near-surface atmospheric variables used to specify the atmospheric boundary condition remain one of the main sources of error. The objective of this research is to constrain the surface forcing function of an ocean model by sea surface temperature (SST) data assimilation. For that purpose, a set of corrections for ERAinterim (hereafter ERAi) reanalysis data is estimated for the period of 1989–2007, using a sequential assimilation method, with ensemble experiments to evaluate the impact of uncertain atmospheric forcing on the ocean state. The control vector of the assimilation method is extended to atmospheric variables to obtain monthly mean parameter corrections by assimilating monthly SST and sea surface salinity (SSS) climatological data in a low resolution global configuration of the NEMO model. In this context, the careful determination of the prior probability distribution of the parameters is an important matter. This paper demonstrates the importance of isolating the impact of forcing errors in the model to perform relevant ensemble experiments. The results obtained for every month of the period between 1989 and 2007 show that the estimated parameters produce the same kind of impact on the SST as the analysis itself. The objective is then to evaluate the long-term time series of the forcing parameters focusing on trends and mean error corrections of air–sea fluxes. Our corrections tend to equilibrate the net heat-flux balance at the global scale (highly positive in ERAi database), and to remove the potentially unrealistic negative trend (leading to ocean cooling) in the ERAi net heat flux over the whole time period. More specifically in the intertropical band, we reduce the warm bias of ERAi data by mostly modifying the latent heat flux by wind speed intensification. Consistently, when used to force the model, the corrected parameters lead to a better agreement between the mean SST produced by the model and mean SST observations over the period of 1989–2007 in the intertropical band.

Mapping flow distortion on oceanographic platforms using computational fluid dynamics - N. O'Sullivan; S. Landwehr; B. Ward (Ocean Science, 2013)
Wind speed measurements over the ocean on ships or buoys are affected by flow distortion from the platform and by the anemometer itself. This can lead to errors in direct measurements and the derived parametrisations. Here we computational fluid dynamics (CFD) to simulate the errors in wind speed measurements caused by flow distortion on the RV Celtic Explorer. Numerical measurements were obtained from the finite-volume CFD code OpenFOAM, which was used to simulate the velocity fields. This was done over a range of orientations in the test domain from −60 to +60° in increments of 10°. The simulation was also set up for a range of velocities, ranging from 5 to 25 m s−1 in increments of 0.5 m s−1. The numerical analysis showed close agreement to experimental measurements.

The transient distributions of nuclear weapon-generated tritium and its decay product 3 He in the Mediterranean Sea, 1952–2011, and their oceanographic potential - W. Roether; P. Jean-Baptiste; E. Fourré; J. Sültenfuß (Ocean Science, 2013)
We present a comprehensive account of tritium and 3He in the Mediterranean Sea since the appearance of the tritium generated by the atmospheric nuclear-weapon testing in the 1950s and early 1960s, based on essentially all available observations. Tritium in surface waters rose to 20–30 TU in 1964 (TU = 1018 × [3H]/H]), a factor of about 100 above the natural level, and thereafter declined 30-fold up to 2011. The decline was largely due to radioactive tritium decay, which produced significant amounts of its stable daughter 3He. We present the scheme by which we separate the tritiugenic part of 3He and the part due to release from the sea floor (terrigenic part). We show that the tritiugenic component can be quantified throughout the Mediterranean waters, typically to a ± 0.15 TU equivalent, mostly because the terrigenic part is low in 3He. This fact makes the Mediterranean unique in offering a potential for the use of tritiugenic 3He as a tracer. The transient distributions of the two tracers are illustrated by a number of sections spanning the entire sea and relevant features of their distributions are noted. By 2011, the 3He concentrations in the top few hundred metres had become low, in response to the decreasing tritium concentrations combined with a flushing out by the general westward drift of these waters. Tritium-3He ages in Levantine Intermediate Water (LIW) were obtained repeated in time at different locations, defining transit times from the LIW source region east of Rhodes. The ages show an upward trend with the time elapsed since the surface-water tritium maximum, which arises because the repeated observations represent increasingly slower moving parts of the full transit time spectrum of LIW. The transit time dispersion revealed by this new application of tritium-3He dating is considerable. We find mean transit times of 12 ± 2 yr up to the Strait of Sicily, 18 ± 3 yr up to the Tyrrhenian Sea, and 22 ± 4 yr up into the Western Mediterranean. Furthermore, we present full Eastern Mediterranean sections of terrigenic 3He and tritium-3He age in 1987, the latter one similarly showing an effect of the transit time dispersion. We conclude that the available tritium and 3He data, particularly if combined with other tracer data, are useful for constraining the subsurface circulation and mixing of the Mediterranean Sea.

Transport of Antarctic bottom water through the Kane Gap, tropical NE Atlantic Ocean - E. G. Morozov; R. Y. Tarakanov; H. van Haren (Ocean Science, 2013)
We study low-frequency properties of the Antarctic Bottom Water (AABW) flow through the Kane Gap (9° N) in the Atlantic Ocean. The measurements in the Kane Gap include five visits with CTD (Conductivity-Temperature-Depth) sections in 2009–2012 and a year-long record of currents on a mooring using three AquaDopp current meters. We found an alternating regime of flow, which changes direction several times during a year. The seasonal signal seems to dominate. The maximum daily average values of southerly velocities reach 0.20 m s−1, while the greatest north-northwesterly velocity is as high as 0.15 m s−1. The velocity and transport at the bottom are aligned along the slope of a local hill near the southwestern side of the gap. The distribution of velocity directions at the upper boundary of AABW is wider. The transport of AABW (Θ

Chaotic variability of the meridional overturning circulation on subannual to interannual timescales - J. J.-M. Hirschi; A. T. Blaker; B. Sinha; A. Coward; B. de Cuevas; S. Alderson; G. Madec (Ocean Science, 2013)
Observations and numerical simulations have shown that the meridional overturning circulation (MOC) exhibits substantial variability on sub- to interannual timescales. This variability is not fully understood. In particular it is not known what fraction of the MOC variability is caused by processes such as mesoscale ocean eddies and waves which are ubiquitous in the ocean. Here we analyse twin experiments performed with a global ocean model at eddying (1/4°) and non-eddying (1°) resolutions. The twin experiments are forced with the same surface fluxes for the 1958 to 2001 period but start from different initial conditions. Our results show that on subannual to interannual timescales a large fraction of MOC variability directly reflects variability in the surface forcing. Nevertheless, in the eddy-permitting case there is an initial-condition-dependent MOC variability (hereinafter referred to as "chaotic" variability) of several Sv (1Sv = 106 m3 s−1) in the Atlantic and the Indo-Pacific. In the Atlantic the chaotic MOC variability represents up to 30% of the total variability at the depths where the maximum MOC occurs. In comparison the chaotic MOC variability is only 5–10% in the non-eddying case. The surface forcing being almost identical in the twin experiments suggests that mesoscale ocean eddies are the most likely cause for the increased chaotic MOC variability in the eddying case. The exact formation time of eddies is determined by the initial conditions which are different in the two model passes, and as a consequence the mesoscale eddy field is decorrelated in the twin experiments. In regions where eddy activity is high in the eddy-permitting model, the correlation of sea surface height variability in the twin runs is close to zero. In the non-eddying case in contrast, we find high correlations (0.9 or higher) over most regions. Looking at the sub- and interannual MOC components separately reveals that most of the chaotic MOC variability is found on subannual timescales for the eddy-permitting model. On interannual timescales the amplitude of the chaotic MOC variability is much smaller and the amplitudes are comparable for both the eddy-permitting and non-eddy-permitting model resolutions. Whereas the chaotic MOC variability on interannual timescales only accounts for a small fraction of the total chaotic MOC variability in the eddy-permitting case, it is the main contributor to the chaotic variability in the non-eddying case away from the Equator.

The Mediterranean Sea system: a review and an introduction to the special issue - T. Tanhua; D. Hainbucher; K. Schroeder; V. Cardin; M. Álvarez; G. Civitarese (Ocean Science, 2013)
The Mediterranean Sea is a semi-enclosed sea characterized by high salinities, temperatures and densities. The net evaporation exceeds the precipitation, driving an anti-estuarine circulation through the Strait of Gibraltar, contributing to very low nutrient concentrations. The Mediterranean Sea has an active overturning circulation, one shallow cell that communicates directly with the Atlantic Ocean, and two deep overturning cells, one in each of the two main basins. It is surrounded by populated areas and is thus sensitive to anthropogenic forcing. Several dramatic changes in the oceanographic and biogeochemical conditions have been observed during the past several decades, emphasizing the need to better monitor and understand the changing conditions and their drivers. During 2011 three oceanographic cruises were conducted in a coordinated fashion in order to produce baseline data of important physical and biogeochemical parameters that can be compared to historic data and be used as reference for future observational campaigns. In this article we provide information on the Mediterranean Sea oceanographic situation, and present a short review that will serve as background information for the special issue in Ocean Science on "Physical, chemical and biological oceanography of the Mediterranean Sea". An important contribution of this article is the set of figures showing the large-scale distributions of physical and chemical properties along the full length of the Mediterranean Sea.

Corrigendum to "NEMO on the shelf: assessment of the Iberia–Biscay–Ireland configuration" published in Ocean Sci., 9, 745–771, 2013 - C. Maraldi; J. Chanut; B. Levier; N. Ayoub; P. De Mey; G. Reffray; F. Lyard; S. Cailleau; M. Drévillon; E. A. Fanjul; M. G. Sotillo; P. Marsaleix (Ocean Science, 2013)

Impact of the Indonesian Throughflow on Agulhas leakage - D. Le Bars; H. A. Dijkstra; W. P. M. De Ruijter (Ocean Science, 2013)
Using ocean models of different complexity we show that opening the Indonesian Passage between the Pacific and the Indian oceans increases the input of Indian Ocean water into the South Atlantic via the Agulhas leakage. In a strongly eddying global ocean model this response results from an increased Agulhas Current transport and a constant proportion of Agulhas retroflection south of Africa. The leakage increases through an increased frequency of ring shedding events. In an idealized two-layer and flat-bottom eddy resolving model, the proportion of the Agulhas Current transport that retroflects is (for a wide range of wind stress forcing) not affected by an opening of the Indonesian Passage. Using a comparison with a linear model and previous work on the retroflection problem, the result is explained as a balance between two mechanisms: decrease retroflection due to large-scale momentum balance and increase due to local barotropic/baroclinic instabilities.

NEMO on the shelf: assessment of the Iberia–Biscay–Ireland configuration - C. Maraldi; J. Chanut; B. Levier; N. Ayoub; P. De Mey; G. Reffray; F. Lyard; S. Cailleau; M. Drévillon; E. A. Fanjul; M. G. Sotillo; P. Marsaleix (Ocean Science, 2013)
This work describes the design and validation of a high-resolution (1/36°) ocean forecasting model over the "Iberian–Biscay–Irish" (IBI) area. The system has been set-up using the NEMO model (Nucleus for European Modelling of the Ocean). New developments have been incorporated in NEMO to make it suitable to open- as well as coastal-ocean modelling. In this paper, we pursue three main objectives: (1) to give an overview of the model configuration used for the simulations; (2) to give a broad-brush account of one particular aspect of this work, namely consistency verification; this type of validation is conducted upstream of the implementation of the system before it is used for production and routinely validated; it is meant to guide model development in identifying gross deficiencies in the modelling of several key physical processes; and (3) to show that such a regional modelling system has potential as a complement to patchy observations (an integrated approach) to give information on non-observed physical quantities and to provide links between observations by identifying broader-scale patterns and processes. We concentrate on the year 2008. We first provide domain-wide consistency verification results in terms of barotropic tides, transports, sea surface temperature and stratification. We then focus on two dynamical subregions: the Celtic shelves and the Bay of Biscay slope and deep regions. The model–data consistency is checked for variables and processes such as tidal currents, tidal fronts, internal tides and residual elevation. We also examine the representation in the model of a seasonal pattern of the Bay of Biscay circulation: the warm extension of the Iberian Poleward Current along the northern Spanish coast (Navidad event) in the winter of 2007–2008.

On the shelf resonances of the English Channel and Irish Sea - D. J. Webb (Ocean Science, 2013)
The resonances of the English Channel and Irish Sea are investigated using the methods of Webb (2012a) together with an Arakawa C-grid model of the region under study. In the semi-diurnal tidal band, the high tides of the Bristol Channel and Gulf of St. Malo are shown to be due to two shelf resonances which strongly couple the two regions. In the diurnal band, the response is complicated by the presence of continental shelf waves.

Intercomparison of the Charnock and COARE bulk wind stress formulations for coastal ocean modelling - J. M. Brown; L. O. Amoudry; F. M. Mercier; A. J. Souza (Ocean Science, 2013)
The accurate parameterisation of momentum and heat transfer across the air–sea interface is vital for realistic simulation of the atmosphere–ocean system. In most modelling applications accurate representation of the wind stress is required to numerically reproduce surge, coastal ocean circulation, surface waves, turbulence and mixing. Different formulations can be implemented and impact the accuracy of the instantaneous and long-term residual circulation, the surface mixed layer, and the generation of wave-surge conditions. This, in turn, affects predictions of storm impact, sediment pathways, and coastal resilience to climate change. The specific numerical formulation needs careful selection to ensure the accuracy of the simulation. Two wind stress parameterisations widely used in the ocean circulation and the storm surge communities respectively are studied with focus on an application to the NW region of the UK. Model–observation validation is performed at two nearshore and one estuarine ADCP (acoustic Doppler current profiler) stations in Liverpool Bay, a hypertidal region of freshwater influence (ROFI) with vast intertidal areas. The period of study covers both calm and extreme conditions to test the robustness of the 10 m wind stress component of the Coupled Ocean–Atmosphere Response Experiment (COARE) bulk formulae and the standard Charnock relation. In this coastal application a realistic barotropic–baroclinic simulation of the circulation and surge elevation is set-up, demonstrating greater accuracy occurs when using the Charnock relation, with a constant Charnock coefficient of 0.0185, for surface wind stress during this one month period.

Transit and residence times in the Adriatic Sea surface as derived from drifter data and Lagrangian numerical simulations - P.-M. Poulain; S. Hariri (Ocean Science, 2013)
Statistics of transit and residence times in the Adriatic Sea surface, a semi-enclosed basin of the Mediterranean, are estimated from drifter data and Lagrangian numerical simulations. The results obtained from the drifters are generally underestimated given their short operating lifetimes (half life of ∼40 days) compared to the transit and residence times. This bias can be removed by considering a large amount of numerical particles whose trajectories are integrated over a long time (750 days) with a statistical advection–dispersion model of the Adriatic surface circulation. Numerical particles indicate that the maximum transit time to exit the basin is about 216–260 days for particles released near the northern tip of the Adriatic, and that a particle entering on the eastern Otranto Channel will typically exit on the other side of the channel after 170–185 days. A duration of 150–168 days is estimated as the residence time in the Adriatic Basin.

Comparing historical and modern methods of sea surface temperature measurement – Part 2: Field comparison in the central tropical Pacific - J. B. R. Matthews; J. B. Matthews (Ocean Science, 2013)
Discrepancies between historical sea surface temperature (SST) datasets have been partly ascribed to use of different adjustments to account for variable measurement methods. Until recently, adjustments had only been applied to bucket temperatures from the late 19th and early 20th centuries, with the aim of correcting their supposed coolness relative to engine cooling water intake temperatures. In the UK Met Office Hadley Centre SST 3 dataset (HadSST3), adjustments have been applied over its full duration to observations from buckets, buoys and engine intakes. Here we investigate uncertainties in the accuracy of such adjustments by direct field comparison of historical and modern methods of shipboard SST measurement. We compare wood, canvas and rubber bucket temperatures to 3 m seawater intake temperature along a central tropical Pacific transect conducted in May and June 2008. We find no average difference between the temperatures obtained with the different bucket types in our short measurement period (∼1 min). Previous field, lab and model experiments have found sizeable temperature change of seawater samples in buckets of smaller volume under longer exposure times. We do, however, report the presence of strong near-surface temperature gradients day and night, indicating that intake and bucket measurements cannot be assumed equivalent in this region. We thus suggest bucket and buoy measurements be considered distinct from intake measurements due to differences in sampling depth. As such, we argue for exclusion of intake temperatures from historical SST datasets and suggest this would likely reduce the need for poorly field-tested bucket adjustments. We also call for improvement in the general quality of intake temperatures from Voluntary Observing Ships. Using a physical model we demonstrate that warming of intake seawater by hot engine room air is an unlikely cause of overly warm intake temperatures. We suggest that reliable correction for such warm errors is not possible since they are largely of unknown origin and can be offset by real near-surface temperature gradients.

Comparing historical and modern methods of sea surface temperature measurement – Part 1: Review of methods, field comparisons and dataset adjustments - J. B. R. Matthews (Ocean Science, 2013)
Sea surface temperature (SST) has been obtained from a variety of different platforms, instruments and depths over the past 150 yr. Modern-day platforms include ships, moored and drifting buoys and satellites. Shipboard methods include temperature measurement of seawater sampled by bucket and flowing through engine cooling water intakes. Here I review SST measurement methods, studies analysing shipboard methods by field or lab experiment and adjustments applied to historical SST datasets to account for variable methods. In general, bucket temperatures have been found to average a few tenths of a °C cooler than simultaneous engine intake temperatures. Field and lab experiments demonstrate that cooling of bucket samples prior to measurement provides a plausible explanation for negative average bucket-intake differences. These can also be credibly attributed to systematic errors in intake temperatures, which have been found to average overly-warm by >0.5 °C on some vessels. However, the precise origin of non-zero average bucket-intake differences reported in field studies is often unclear, given that additional temperatures to those from the buckets and intakes have rarely been obtained. Supplementary accurate in situ temperatures are required to reveal individual errors in bucket and intake temperatures, and the role of near-surface temperature gradients. There is a need for further field experiments of the type reported in Part 2 to address this and other limitations of previous studies.

Corrigendum to "Biogeography of planktonic bacterial communities across the whole Mediterranean Sea" published in Ocean Sci., 9, 585–595, 2013 - F. Mapelli; M. M. Varela; M. Barbato; R. Alvariño; V. Hernando-Morales; M. Fusi; M. Álvarez; G. Merlino; D. Daffonchio; E. Teira; S. Borin (Ocean Science, 2013)

Daily scale wintertime sea surface temperature and IPC-Navidad variability in the southern Bay of Biscay from 1981 to 2010 - G. Esnaola; J. Sáenz; E. Zorita; A. Fontán; V. Valencia; P. Lazure (Ocean Science, 2013)
The combination of remotely sensed gappy Sea surface temperature (SST) images with the missing data filling DINEOF (data interpolating empirical orthogonal functions) technique, followed by a principal component analysis of the reconstructed data, has been used to identify the time evolution and the daily scale variability of the wintertime surface signal of the Iberian Poleward Current (IPC), or Navidad, during the 1981–2010 period. An exhaustive comparison with the existing bibliography, and the vertical temperature and salinity profiles related to its extremes over the Bay of Biscay area, show that the obtained time series accurately reflect the IPC-Navidad variability. Once a time series for the evolution of the SST signal of the current over the last decades is well established, this time series is used to propose a physical mechanism in relation to the variability of the IPC-Navidad, involving both atmospheric and oceanic variables. According to the proposed mechanism, an atmospheric circulation anomaly observed in both the 500 hPa and the surface levels generates atmospheric surface level pressure, wind-stress and heat-flux anomalies. In turn, those surface level atmospheric anomalies induce mutually coherent SST and sea level anomalies over the North Atlantic area, and locally, in the Bay of Biscay area. These anomalies, both locally over the Bay of Biscay area and over the North Atlantic, are in agreement with several mechanisms that have separately been related to the variability of the IPC-Navidad, i.e. the south-westerly winds, the joint effect of baroclinicity and relief (JEBAR) effect, the topographic β effect and a weakened North Atlantic gyre.

Combining in situ measurements and altimetry to estimate volume, heat and salt transport variability through the Faroe–Shetland Channel - B. Berx; B. Hansen; S. Østerhus; K. M. Larsen; T. Sherwin; K. Jochumsen (Ocean Science, 2013)
From 1994 to 2011, instruments measuring ocean currents (Acoustic Doppler Current Profilers; ADCPs) have been moored on a section crossing the Faroe–Shetland Channel. Together with CTD (Conductivity Temperature Depth) measurements from regular research vessel occupations, they describe the flow field and water mass structure in the channel. Here, we use these data to calculate the average volume transport and properties of the flow of warm water through the channel from the Atlantic towards the Arctic, termed the Atlantic inflow. We find the average volume transport of this flow to be 2.7 ± 0.5 Sv (1 Sv = 106 m3 s–1) between the shelf edge on the Faroe side and the 150 m isobath on the Shetland side. The average heat transport (relative to 0 °C) was estimated to be 107 ± 21 TW (1 TW = 1012 W) and the average salt import to be 98 ± 20 × 106 kg s−1. Transport values for individual months, based on the ADCP data, include a large level of variability, but can be used to calibrate sea level height data from satellite altimetry. In this way, a time series of volume transport has been generated back to the beginning of satellite altimetry in December 1992. The Atlantic inflow has a seasonal variation in volume transport that peaks around the turn of the year and has an amplitude of 0.7 Sv. The Atlantic inflow has become warmer and more saline since 1994, but no equivalent trend in volume transport was observed.

Image of a subsurface current core in the southern South China Sea - Q. S. Tang; D. X. Wang; J. B. Li; P. Yan; J. Li (Ocean Science, 2013)
A legacy seismic transect acquired on 30 and 31 May 2009 in the southern South China Sea (SCS) was reprocessed to reveal the thermohaline structure of the water column. In the study region, a mesoscale subsurface lens with extraordinary features was detected at 113.5° E, 11.5° N. It is centred at 450 m depth, occupies both the subsurface and intermediate water from 250 to 600 m, and has an intersection diameter of around 60 km. The simulated results from Hybrid Coordinate Ocean Model reveal an eddy-induced subsurface current running southwestward along the deep basin edge and suggest that the imaged lens is a snapshot of the subsurface current core rather than a subsurface eddy.

A comparison between gradient descent and stochastic approaches for parameter optimization of a sea ice model - H. Sumata; F. Kauker; R. Gerdes; C. Köberle; M. Karcher (Ocean Science, 2013)
Two types of optimization methods were applied to a parameter optimization problem in a coupled ocean–sea ice model of the Arctic, and applicability and efficiency of the respective methods were examined. One optimization utilizes a finite difference (FD) method based on a traditional gradient descent approach, while the other adopts a micro-genetic algorithm (μGA) as an example of a stochastic approach. The optimizations were performed by minimizing a cost function composed of model–data misfit of ice concentration, ice drift velocity and ice thickness. A series of optimizations were conducted that differ in the model formulation ("smoothed code" versus standard code) with respect to the FD method and in the population size and number of possibilities with respect to the μGA method. The FD method fails to estimate optimal parameters due to the ill-shaped nature of the cost function caused by the strong non-linearity of the system, whereas the genetic algorithms can effectively estimate near optimal parameters. The results of the study indicate that the sophisticated stochastic approach (μGA) is of practical use for parameter optimization of a coupled ocean–sea ice model with a medium-sized horizontal resolution of 50 km × 50 km as used in this study.

Wave-turbulence scaling in the ocean mixed layer - G. Sutherland; B. Ward; K. H. Christensen (Ocean Science, 2013)
Microstructure measurements were collected using an autonomous freely rising profiler under a variety of different atmospheric forcing and sea states in the open ocean. Here, profiles of turbulent kinetic energy dissipation rate, ε, are compared with various proposed scalings. In the oceanic boundary layer, the depth dependence of ε was found to be largely consistent with that expected for a shear-driven wall layer. This is in contrast with many recent studies which suggest higher rates of turbulent kinetic energy dissipation in the near surface of the ocean. However, some dissipation profiles appeared to scale with the sum of the wind and swell generated Stokes shear with this scaling extending beyond the mixed layer depth. Integrating ε in the mixed layer yielded results that 1% of the wind power referenced to 10 m is being dissipated here.

Biogeography of planktonic bacterial communities across the whole Mediterranean Sea - F. Mapelli; M. M. Varela; M. Barbato; R. Alvariño; M. Fusi; M. Álvarez; G. Merlino; D. Daffonchio; S. Borin (Ocean Science, 2013)
Bacterial population distribution in the Mediterranean Sea has been mainly studied by considering small geographical areas or specific phylogenetic groups. The present study is a molecular microbiology investigation aimed to identify the environmental factors driving total bacterioplankton community composition of seawater samples collected along a transect covering the whole Mediterranean Sea. We performed automated ribosomal intergenic sequence analysis (ARISA) and microscope evaluation of prokaryotic abundance of seawater sampled across both vertical profiles and longitudinal transects in the whole basin. Prokaryotic abundance decreased with depth at all the stations and presented similar values in epi-, meso- and bathypelagic layers across the whole Mediterranean Sea. However, peculiar bacterial assemblages were selected along a longitudinal transect in the epipelagic layers of the eastern and western sub-basins. Vertical differences of the bacterial communities were observed only when considering the epi- and bathypelagic waters, while the study of the structure of bacterial communities at a finer scale across the water column displayed higher variability at the intermediate layers. Nonetheless, different physico-chemical factors were significantly related to the distribution of bacterial populations, both according to geographic position and down the water column in the whole Mediterranean Sea. These results demonstrated that bacterial assemblages are putatively correlated to different water masses of the complex hydrographical systems of the eastern and western Mediterranean sub-basins.

Investigation of saline water intrusions into the Curonian Lagoon (Lithuania) and two-layer flow in the Klaipėda Strait using finite element hydrodynamic model - P. Zemlys; C. Ferrarin; G. Umgiesser; S. Gulbinskas; D. Bellafiore (Ocean Science, 2013)
This work is focused on the application of a modelling system to simulate 3-D interaction between the Curonian Lagoon and the Baltic Sea coastal waters and to reflect spatiotemporal dynamics of marine waters in the Curonian Lagoon. The model system is based on the finite element programme package SHYFEM which can be used to resolve the hydrodynamic equations in lagoons, coastal seas, estuaries and lakes. The results of a one year (2009) 3-D model simulation with real weather and hydrological forcing show that the saline water intrusions from the sea through Klaipėda Strait are gradually decreasing with distance from the sea and become negligible (average annual salinity about 0.57‰) at a distance of about 20 km to the south of Kiaulės Nugara island. Analyses of the simulation results also show this area to be highly heterogeneous according to the vertical salinity distribution. While in the deeper Klaipėda Strait (harbour waterway) differences in average salinity between near bottom and surface layers varies in the range 2–2.5‰, in the rest of the Curonian Lagoon it is less than 0.5‰. The exchange flow showed vertical structure, but was horizontally uniform with the presence of a two-directional flow that from time to time changes to either saline water one-directional flow to the Curonian Lagoon or fresh water one-directional flow to the sea. Two-directional flow duration decreases with a distance from sea entrance in Klaipėda Strait from around 180 days yr−1 close to the sea entrance to 50 days yr−1 just behind Kiaulės Nugara island. One-directional outflow duration is increasing with a distance from the sea entrance from 100 to 225 days yr−1. One-directional inflow duration occurs in the range of 70–100 days yr−1. The analysis of the ratio of buoyancy layer thickness to water depth (hb/H) and the Wedderburn number identified the main importance of wind action on the flow structure. Strong winds from the North and NW determine a barotropic inflow which is mostly responsible for the salt water intrusion into the Curonian Lagoon. Absence of wind or cross-strait wind regimes allows the maintenance of a two-layer flow typical of estuarine dynamics.

Exceptional dense water formation on the Adriatic shelf in the winter of 2012 - H. Mihanović; I. Vilibić; S. Carniel; M. Tudor; A. Russo; A. Bergamasco; N. Bubić; Z. Ljubešić; D. Viličić; A. Boldrin; V. Malačič; M. Celio; et al. (Ocean Science, 2013)
In this paper we document dense water formation throughout the Adriatic shelf and coastal area in January/February 2012, resulting in record-breaking densities observed during and after the event. The unprecedented dense water generation was preconditioned by a dry and warm year which resulted in a significant reduction of coastal freshwaters, superimposed on a long-term basin-wide salinity increase. The final event that triggered the dense water formation was an extended period of cold weather with strong and severe winds. Record-breaking potential density anomalies (above 30 kg m−3) were measured at several formation sites. Accumulated surface net heat and water losses in some coastal regions exceeded 1.5 GJ m−2 and 250 kg m−2 over 21 days, respectively. Excessiveness, importance of shelf-type dense water formation and effects on the thermohaline circulation and deep aquatic systems are discussed.

Ocean state indicators from MyOcean altimeter products - L. Bessières; M. H. Rio; C. Dufau; C. Boone; M. I. Pujol (Ocean Science, 2013)
The European MyOcean project (http://www.myocean.eu.org) provides observations of the ocean dynamic topography from altimeter measurements. Three specific indicators have been developed, based on altimeter data only, in order to monitor the ocean state. The first ocean indicator observes the positive and negative phases of the ENSO events in the tropical Pacific, the El Niño/La Niña events, since 1992. The second ocean indicator tracks the contracted or extended state of the Kuroshio Extension. The last ocean indicator is dedicated to the Ionian Basin in the Mediterranean Sea and permits separation of "zonal-cyclonic" state (1998–2005 and since 2011 up to now) from the "anticyclonic" state (1993–1996) usually discussed in the literature. In addition, it allows identifying a third state in which both the anticyclonic circulation around the northern part of the basin and the strong zonal Mid-Ionian Jet co-exist (2008–2010).

ENSO components of the Atlantic multidecadal oscillation and their relation to North Atlantic interannual coastal sea level anomalies - J. Park; G. Dusek (Ocean Science, 2013)
The El Niño Southern Oscillation (ENSO) and the Atlantic Multidecadal Oscillation (AMO) are known to influence coastal water levels along the East Coast of the United States. By identifying empirical orthogonal functions (EOFs), which coherently contribute from the Multivariate ENSO Index (MEI) to the AMO index (AMOI), we characterize both the expression of ENSO in the unsmoothed AMOI, and coherent relationships between these indices and interannual sea level anomalies at six stations in the Gulf of Mexico and western North Atlantic. Within the ENSO band (2–7 yr periods) the total contribution of MEI to unsmoothed AMOI variability is 79%. Cross correlation suggests that the MEI leads expression of the ENSO signature in the AMOI by six months, consistent with the mechanism of an atmospheric bridge. Within the ENSO band, essentially all of the coupling between the unsmoothed AMOI and sea level anomalies is the result of ENSO expression in the AMOI. At longer periods we find decadal components of sea level anomalies linked to the AMOI at three southern stations (Key West, Pensacola, Charleston), but not at the northern stations (Baltimore, Boston, Portland), with values of coherence ranging from 20 to 50%. The coherence of MEI to coastal sea level anomalies has a different structure and is generally weaker than that of the ENSO expressed AMOI influence, suggesting distinct physical mechanisms are influencing sea level anomalies due to a direct ENSO teleconnection when compared to teleconnections based on ENSO expression in the AMOI. It is expected that applying this analysis to extremes of sea level anomalies will reveal additional influences.

Assessment of MERIS ocean color data products for European seas - G. Zibordi; F. Mélin; J.-F. Berthon; E. Canuti (Ocean Science, 2013)
The accuracy of marine data products from the Medium Resolution Imaging Spectrometer (MERIS) operated on board the Envisat platform is investigated with the aid of in situ geographically distributed measurements from different European seas. The assessment focuses on standard products from the 2012 data update commonly identified as 3rd Reprocessing. Results indicate atmospherically corrected data affected by a negative bias of several tens percent at the 413 nm center wavelength, significantly decreasing to a few percent at 560 nm and increasing again at 665 nm. Such an underestimate at the blue center wavelengths leads to an average overestimate of the algal-1 MERIS pigment index largely exceeding 100% for the considered European seas. A comparable overestimate is also observed for the algal-2 pigment index independently determined from top-of-atmosphere radiance through the application of neural networks.

Recirculation in the Fram Strait and transports of water in and north of the Fram Strait derived from CTD data - M. Marnela; B. Rudels; M.-N. Houssais; A. Beszczynska-Möller; P. B. Eriksson (Ocean Science, 2013)
The volume, heat and freshwater transports in the Fram Strait are estimated from geostrophic computations based on summer hydrographic data from 1984, 1997, 2002 and 2004. In these years, in addition to the usually sampled section along 79° N, a section between Greenland and Svalbard was sampled further north. Quasi-closed boxes bounded by the two sections and Greenland and Svalbard can then be formed. Applying conservation constraints on these boxes provides barotropic reference velocities. The net volume flux is southward and varies between 2 and 4 Sv. The recirculation of Atlantic water is about 2 Sv. Heat is lost to the atmosphere and the heat loss from the area between the sections averaged over the four years is about 10 TW. The net heat (temperature) transport is 20 TW northward into the Arctic Ocean, with large interannual differences. The mean net freshwater added between the sections is 40 mSv and the mean freshwater transport southward across 79° N is less than 60 mSv, indicating that most of the liquid freshwater leaving the Arctic Ocean through Fram Strait in summer is derived from sea ice melt in the northern vicinity of the strait. In 1997, 2001 and 2003 meridional sections along 0° longitude were sampled and in 2003 two smaller boxes can be formed, and the recirculation of Atlantic water in the strait is estimated by geostrophic computations and continuity constraints. The recirculation is weaker close to 80° N than close to 78° N, indicating that the recirculation is mainly confined to the south of 80° N. This is supported by the observations in 1997 and 2001, when only the northern part of the meridional section, from 79° N to 80° N, can be computed with the constraints applied. The recirculation is found strongest close to 79° N.

Springtime contribution of dinitrogen fixation to primary production across the Mediterranean Sea - E. Rahav; B. Herut; A. Levi; M. R. Mulholland; I. Berman-Frank (Ocean Science, 2013)
Dinitrogen (N2) fixation rates were measured during early spring across the different provinces of Mediterranean Sea surface waters. N2 fixation rates, measured using 15N2 enriched seawater, were lowest in the eastern basin and increased westward with a maximum at the Strait of Gibraltar (0.10 to 2.35 nmol N L−1 d−1, respectively). These rates were 3–7 fold higher than N2 fixation rates measured previously in the Mediterranean Sea during summertime and we estimated that methodological differences alone did not account for the seasonal changes we observed. Higher contribution of N2 fixation to primary production (4–8%) was measured in the western basin compared to the eastern basin (∼2%). Our data indicates that these differences between basins may be attributed to changes in N2-fixing planktonic communities and that heterotrophic diazotrophy may play a significant role in the eastern Mediterranean while autotrophic diazotrophy has a more dominant role in the western basin.

MERIS-based ocean colour classification with the discrete Forel–Ule scale - M. R. Wernand; A. Hommersom; H. J. van der Woerd (Ocean Science, 2013)
Multispectral information from satellite borne ocean colour sensors is at present used to characterize natural waters via the retrieval of concentrations of the three dominant optical constituents; pigments of phytoplankton, non-algal particles and coloured dissolved organic matter. A limitation of this approach is that accurate retrieval of these constituents requires detailed local knowledge of the specific absorption and scattering properties. In addition, the retrieval algorithms generally use only a limited part of the collected spectral information. In this paper we present an additional new algorithm that has the merit of using the full spectral information in the visible domain to characterize natural waters in a simple and globally valid way. This Forel–Ule MERIS (FUME) algorithm converts the normalized multiband reflectance information into a discrete set of numbers using uniform colourimetric functions. The Forel–Ule (FU) scale is a sea colour comparator scale that has been developed to cover all possible natural sea colours, ranging from indigo blue (the open ocean) to brownish-green (coastal water) and even brown (humic-acid dominated) waters. Data using this scale have been collected since the late nineteenth century, and therefore, this algorithm creates the possibility to compare historic ocean colour data with present-day satellite ocean colour observations. The FUME algorithm was tested by transforming a number of MERIS satellite images into Forel–Ule colour index images and comparing in situ observed FU numbers with FU numbers modelled from in situ radiometer measurements. Similar patterns and FU numbers were observed when comparing MERIS ocean colour distribution maps with ground truth Forel–Ule observations. The FU numbers modelled from in situ radiometer measurements showed a good correlation with observed FU numbers (R2 = 0.81 when full spectra are used and R2 = 0.71 when MERIS bands are used).

A new 3-D modelling method to extract subtransect dimensions from underwater videos - L. Fillinger; T. Funke (Ocean Science, 2013)
Underwater video transects have become a common tool for quantitative analysis of the seafloor. However a major difficulty remains in the accurate determination of the area surveyed as underwater navigation can be unreliable and image scaling does not always compensate for distortions due to perspective and topography. Depending on the camera set-up and available instruments, different methods of surface measurement are applied, which make it difficult to compare data obtained by different vehicles. 3-D modelling of the seafloor based on 2-D video data and a reference scale can be used to compute subtransect dimensions. Focussing on the length of the subtransect, the data obtained from 3-D models created with the software PhotoModeler Scanner are compared with those determined from underwater acoustic positioning (ultra short baseline, USBL) and bottom tracking (Doppler velocity log, DVL). 3-D model building and scaling was successfully conducted on all three tested set-ups and the distortion of the reference scales due to substrate roughness was identified as the main source of imprecision. Acoustic positioning was generally inaccurate and bottom tracking unreliable on rough terrain. Subtransect lengths assessed with PhotoModeler were on average 20% longer than those derived from acoustic positioning due to the higher spatial resolution and the inclusion of slope. On a high relief wall bottom tracking and 3-D modelling yielded similar results. At present, 3-D modelling is the most powerful, albeit the most time-consuming, method for accurate determination of video subtransect dimensions.

Large-scale temperature and salinity changes in the upper Canadian Basin of the Arctic Ocean at a time of a drastic Arctic Oscillation inversion - P. Bourgain; J. C. Gascard; J. Shi; J. Zhao (Ocean Science, 2013)
Between 2008 and 2010, the Arctic Oscillation index over Arctic regions shifted from positive values corresponding to more cyclonic conditions prevailing during the 4th International Polar Year (IPY) period (2007–2008) to extremely negative values corresponding to strong anticyclonic conditions in 2010. In this context, we investigated the recent large-scale evolution of the upper western Arctic Ocean, based on temperature and salinity summertime observations collected during icebreaker campaigns and from ice-tethered profilers (ITPs) drifting across the region in 2008 and 2010. Particularly, we focused on (1) the freshwater content which was extensively studied during previous years, (2) the near-surface temperature maximum due to incoming solar radiation, and (3) the water masses advected from the Pacific Ocean into the Arctic Ocean. The observations revealed a freshwater content change in the Canadian Basin during this time period. South of 80° N, the freshwater content increased, while north of 80° N, less freshening occurred in 2010 compared to 2008. This was more likely due to the strong anticyclonicity characteristic of a low AO index mode that enhanced both a wind-generated Ekman pumping in the Beaufort Gyre and a possible diversion of the Siberian River runoff toward the Eurasian Basin at the same time. The near-surface temperature maximum due to incoming solar radiation was almost 1 °C colder in the southern Canada Basin (south of 75° N) in 2010 compared to 2008, which contrasted with the positive trend observed during previous years. This was more likely due to higher summer sea ice concentration in 2010 compared to 2008 in that region, and surface albedo feedback reflecting more sun radiation back in space. The Pacific water (PaW) was also subjected to strong spatial and temporal variability between 2008 and 2010. In the Canada Basin, both summer and winter PaW signatures were stronger between 75° N and 80° N. This was more likely due to a strong recirculation within the Beaufort Gyre. In contrast, south of 75° N, the cooling and warming of the summer and winter PaW, respectively, suggest that either the PaW was less present in 2010 than in 2008 in this region, and/or the PaW was older in 2010 than in 2008. In addition, in the vicinity of the Chukchi Sea, both summer and winter PaW were significantly warmer in 2010 than in 2008, as a consequence of a general warming trend of the PaW entering in the deep Arctic Ocean as of 2008.

From the chlorophyll a in the surface layer to its vertical profile: a Greenland Sea relationship for satellite applications - A. Cherkasheva; E.-M. Nöthig; E. Bauerfeind; C. Melsheimer; A. Bracher (Ocean Science, 2013)
Current estimates of global marine primary production range over a factor of two. Improving these estimates requires an accurate knowledge of the chlorophyll vertical profiles, since they are the basis for most primary production models. At high latitudes, the uncertainty in primary production estimates is larger than globally, because here phytoplankton absorption shows specific characteristics due to the low-light adaptation, and in situ data and ocean colour observations are scarce. To date, studies describing the typical chlorophyll profile based on the chlorophyll in the surface layer have not included the Arctic region, or, if it was included, the dependence of the profile shape on surface concentration was neglected. The goal of our study was to derive and describe the typical Greenland Sea chlorophyll profiles, categorized according to the chlorophyll concentration in the surface layer and further monthly resolved profiles. The Greenland Sea was chosen because it is known to be one of the most productive regions of the Arctic and is among the regions in the Arctic where most chlorophyll field data are available. Our database contained 1199 chlorophyll profiles from R/Vs Polarstern and Maria S. Merian cruises combined with data from the ARCSS-PP database (Arctic primary production in situ database) for the years 1957–2010. The profiles were categorized according to their mean concentration in the surface layer, and then monthly median profiles within each category were calculated. The category with the surface layer chlorophyll (CHL) exceeding 0.7 mg C m−3 showed values gradually decreasing from April to August. A similar seasonal pattern was observed when monthly profiles were averaged over all the surface CHL concentrations. The maxima of all chlorophyll profiles moved from the greater depths to the surface from spring to late summer respectively. The profiles with the smallest surface values always showed a subsurface chlorophyll maximum with its median magnitude reaching up to three times the surface concentration. While the variability of the Greenland Sea season in April, May and June followed the global non-monthly resolved relationship of the chlorophyll profile to surface chlorophyll concentrations described by the model of Morel and Berthon (1989), it deviated significantly from the model in the other months (July–September), when the maxima of the chlorophyll are at quite different depths. The Greenland Sea dimensionless monthly median profiles intersected roughly at one common depth within each category. By applying a Gaussian fit with 0.1 mg C m−3 surface chlorophyll steps to the median monthly resolved chlorophyll profiles of the defined categories, mathematical approximations were determined. They generally reproduce the magnitude and position of the CHL maximum, resulting in an average 4% underestimation in Ctot (and 2% in rough primary production estimates) when compared to in situ estimates. These mathematical approximations can be used as the input to the satellite-based primary production models that estimate primary production in the Arctic regions.

Seasonality of intermediate waters hydrography west of the Iberian Peninsula from an 8 yr semiannual time series of an oceanographic section - E. Prieto; C. González-Pola; A. Lav& iacute; n; R. F. Sánchez; M. Ruiz-Villarreal (Ocean Science, 2013)
Seasonality of hydrographical properties at depth in the western Iberian margin (eastern North Atlantic) is analysed from a 2003–2010 time series of a semiannual oceanographic section extending ∼200 nm off Cape Finisterre (43° N). All water masses down to the permanent thermocline (2000 dbar) show a consistent seasonal signature in their thermohaline properties and there is a notable asymmetry between the slope region and the outer ocean (in the surroundings of the Galicia Bank). There is overall cooling and freshening of eastern North Atlantic central waters in summertime, which is larger and deeper-reaching on the slope. In summertime, Mediterranean Water (MW) gets tightly attached against the slope and is uplifted, reinforcing its thermohaline signature and diminishing its presence at the outer ocean. In wintertime the situation reverses, MW seems to detach from the slope and spreads out to the open ocean, even being observed a secondary branch around the Galicia Bank. Thermohaline seasonality at depth shows values up to 0.4 °C and 0.08 in salinity at the lower MW, of the order of 20% of the overall interannual variability observed during the whole period. Decomposition of thermohaline changes at isobaric levels to changes along isoneutral surfaces and changes due to vertical displacements help analyse the physical processes behind the observed seasonality in terms of (1) the large-scale seasonality of the subtropical gyre in response to the seasonal migration of the subtropical high pressure system and subsequent anomalies in Ekman transport and wind stress curl, (2) the continental slope dynamics, characterized by summer upwelling, winter development of the Iberian Poleward Current and Mediterranean water spreading, and (3) the possible influence of seasonal changes of water mass properties at their formation sources.

Variability in the air–sea interaction patterns and timescales within the south-eastern Bay of Biscay, as observed by HF radar data - A. Fontán; G. Esnaola; J. Sáenz; M. González (Ocean Science, 2013)
Two high-frequency (HF) radar stations were installed on the coast of the south-eastern Bay of Biscay in 2009, providing high spatial and temporal resolution and large spatial coverage of currents in the area for the first time. This has made it possible to quantitatively assess the air–sea interaction patterns and timescales for the period 2009–2010. The analysis was conducted using the Barnett–Preisendorfer approach to canonical correlation analysis (CCA) of reanalysis surface winds and HF radar-derived surface currents. The CCA yields two canonical patterns: the first wind–current interaction pattern corresponds to the classical Ekman drift at the sea surface, whilst the second describes an anticyclonic/cyclonic surface circulation. The results obtained demonstrate that local winds play an important role in driving the upper water circulation. The wind–current interaction timescales are mainly related to diurnal breezes and synoptic variability. In particular, the breezes force diurnal currents in waters of the continental shelf and slope of the south-eastern Bay. It is concluded that the breezes may force diurnal currents over considerably wider areas than that covered by the HF radar, considering that the northern and southern continental shelves of the Bay exhibit stronger diurnal than annual wind amplitudes.

On the use of the Stokes number to explain frictional tidal dynamics and water column structure in shelf seas - A. J. Souza (Ocean Science, 2013)
In recent years coastal oceanographers have suggested using the "Strouhal" number or its inverse, the "Stokes" number, to describe the effect of bottom boundary layer turbulence on the vertical structure of both density and currents. These are defined as the ratios of the frictional depth (δ) to the water column depth (h) or vice versa. Although many researchers have mentioned that the effects of the earth's rotation should be important, they have tended to omit it. Rotation may have an important influence on tidal currents, as the frictional depth from a fully cyclonic to a fully anticyclonic tidal ellipse can vary by up to an order of magnitude at mid latitudes. The Stokes number might appear smaller for cyclonic current ellipses (larger for anticyclonic) than it is without rotation, resulting in frictional effects being underestimated (overestimated). Here, a way to calculate a Stokes number is proposed, in which the effect of the earth's rotation is taken into account. The standard Stokes and the rotational Stokes numbers are used as predictors for the position of the tidal mixing fronts in the Irish Sea. Results show that use of the rotational number improves the predictions of fronts in shallow cyclonic areas of the eastern Irish Sea. This suggests that the effect of rotation on the water column structure will be more important in shallow shelf seas and estuaries with strong rotational currents.

The effect of various vertical discretization schemes and horizontal diffusion parameterization on the performance of a 3-D ocean model: the Black Sea case study - G. Shapiro; M. Luneva; J. Pickering; D. Storkey (Ocean Science, 2013)
Results of a sensitivity study are presented from various configurations of the NEMO ocean model in the Black Sea. The standard choices of vertical discretization, viz. z levels, s coordinates and enveloped s coordinates, all show their limitations in the areas of complex topography. Two new hybrid vertical coordinate schemes are presented: the "s-on-top-of-z" and its enveloped version. The hybrid grids use s coordinates or enveloped s coordinates in the upper layer, from the sea surface to the depth of the shelf break, and z-coordinates are set below this level. The study is carried out for a number of idealised and real world settings. The hybrid schemes help reduce errors generated by the standard schemes in the areas of steep topography. Results of sensitivity tests with various horizontal diffusion formulations are used to identify the optimum value of Smagorinsky diffusivity coefficient to best represent the mesoscale activity.

Numerical tools to estimate the flux of a gas across the air–water interface and assess the heterogeneity of its forcing functions - V. M. N. C. S. Vieira; F. Martins; J. Silva; R. Santos (Ocean Science, 2013)
A numerical tool was developed for the estimation of gas fluxes across the air–water interface. The primary objective is to use it to estimate CO2 fluxes. Nevertheless application to other gases is easily accomplished by changing the values of the parameters related to the physical properties of the gases. A user-friendly software was developed allowing to build upon a standard kernel a custom-made gas flux model with the preferred parameterizations. These include single or double layer models; several numerical schemes for the effects of wind in the air-side and water-side transfer velocities; the effects of atmospheric stability, surface roughness and turbulence from current drag with the bottom; and the effects on solubility of water temperature, salinity, air temperature and pressure. An analysis was also developed which decomposes the difference between the fluxes in a reference situation and in alternative situations into its several forcing functions. This analysis relies on the Taylor expansion of the gas flux model, requiring the numerical estimation of partial derivatives by a multivariate version of the collocation polynomial. Both the flux model and the difference decomposition analysis were tested with data taken from surveys done in the lagoon system of Ria Formosa, south Portugal, in which the CO2 fluxes were estimated using the infrared gas analyzer (IRGA) and floating chamber method, whereas the CO2 concentrations were estimated using the IRGA and degasification chamber. Observations and estimations show a remarkable fit.

Influence of winds on temporally varying short and long period gravity waves in the near shore regions of the eastern Arabian Sea - J. Glejin; V. Sanil Kumar; T. M. Balakrishnan Nair; J. Singh (Ocean Science, 2013)
Wave data collected off Ratnagiri, west coast of India, during 1 May 2010 to 30 April 2012 are used in this study. Seasonal and annual variations in wave data controlled by the local wind system such as sea breeze and land breeze, and remote wind generated long period waves are also studied. The role of sea breeze on the sea state during pre- and postmonsoon seasons is studied and it is found that the maximum wave height is observed at 15:00 UTC during the premonsoon season, with an estimated difference in time lag of 1–2 h in maximum wave height between premonsoon and postmonsoon seasons. Observed waves are classified in to (i) short waves (Tp 13 s) based on peak period (Tp) and the percentages of occurrence of each category are estimated. Long period waves are observed mainly during the pre- and the postmonsoon seasons. During the southwest monsoon period, the waves with period > 13 s are a minimum. An event during 2011 is identified as swells propagated from the Southern Ocean with an estimated travelling time of 5–6 days. The swells reaching the Arabian Sea from the south Indian Ocean and Southern Ocean, due to storms during the pre- and postmonsoon periods, modify the near surface winds due to higher phase wave celerity than the wind speed. Estimation of inverse wave age using large-scale winds such as NCEP (National Centers for Environmental Prediction) reflects the presence of cyclonic activity during pre- and postmonsoon seasons but not the effect of the local sea breeze/land breeze wind system.

Sea wave modeling with X-band COSMO-SkyMed © SAR-derived wind field forcing and applications in coastal vulnerability assessment - G. Benassai; A. Montuori; M. Migliaccio; F. Nunziata (Ocean Science, 2013)
In this paper, X-band COSMO-SkyMed© synthetic aperture radar (SAR) wind field data are first used to force coastal wind wave modeling for both sea wave numerical simulation and coastal vulnerability assessment purposes. The SAR-based wind field retrieval is accomplished by resolving the SAR-based wind speed and wind direction retrieval problems independently. The sea surface wind speed is retrieved through the azimuth cut-off procedure, and the sea surface wind direction is determined by the multi-resolution analysis of the discrete wavelet transform. The wind wave modeling is based on the third-generation Simulating WAves Nearshore (SWAN) model, which is used for sea wave state estimation in coastal and inland regions. The coastal vulnerability assessment is provided by means of a key parameter, known as impact index, which evaluates the coastal risk due to the inundation of the inshore land. Experiments consist of SWAN numerical simulations run with respect to some relevant wave storms recorded in the southern Tyrrhenian Sea on 2010, with applications in coastal vulnerability assessment along the Sele coastal plain. Experimental results show the benefits of blended wind field products, provided by European Centre for Medium Weather Forecast (ECMWF) model winds and SAR-based wind field estimations, for both wind wave modeling and coastal vulnerability assessment purposes.

Long-term monitoring programme of the hydrological variability in the Mediterranean Sea: a first overview of the HYDROCHANGES network - K. Schroeder; C. Millot; L. Bengara; S. Ben Ismail; M. Bensi; M. Borghini; G. Budillon; V. Cardin; L. Coppola; C. Curtil; A. Drago; B. El Moumni; et al. (Ocean Science, 2013)
The long-term monitoring of basic hydrological parameters (temperature and salinity), collected as time series with adequate temporal resolution (i.e. with a sampling interval allowing the resolution of all important timescales) in key places of the Mediterranean Sea (straits and channels, zones of dense water formation, deep parts of the basins), constitute a priority in the context of global changes. This led CIESM (The Mediterranean Science Commission) to support, since 2002, the HYDROCHANGES programme (http//www.ciesm.org/marine/programs/hydrochanges.htm), a network of autonomous conductivity, temperature, and depth (CTD) sensors, deployed on mainly short and easily manageable subsurface moorings, within the core of a certain water mass. The HYDROCHANGES strategy is twofold and develops on different scales. To get information about long-term changes of hydrological characteristics, long time series are needed. But before these series are long enough they allow the detection of links between them at shorter timescales that may provide extremely valuable information about the functioning of the Mediterranean Sea. The aim of this paper is to present the history of the programme and the current set-up of the network (monitored sites, involved groups) as well as to provide for the first time an overview of all the time series collected under the HYDROCHANGES umbrella, discussing the results obtained thanks to the programme.

Sea level trend and variability in the Singapore Strait - P. Tkalich; P. Vethamony; Q.-H. Luu; M. T. Babu (Ocean Science, 2013)
Sea level in the Singapore Strait (SS) exhibits response to various scale phenomena, from local to global. Longest tide gauge records in SS are analysed to derive local sea level trend and annual, inter-annual and multi-decadal sea level variability, which then are attributed to regional and global phenomena. Annual data gaps are reconstructed using functions correlating sea level variability with ENSO. At annual scale, sea level anomalies in SS are (quasi-periodic) monsoon-driven, of the order of ±20 cm, the highest during northeast monsoon and the lowest during southwest monsoon. Interannual regional sea level drops are associated with El Niño events, while the rises are correlated with La Niña episodes; both variations are in the range of ±5 cm. At multi-decadal scale, annual measured sea levels in SS are varying with global mean sea level, rising at the rate 1.2–1.7 mm yr−1 for 1975–2009, 1.8–2.3 mm yr−1 for 1984–2009 and 1.9–4.6 mm yr−1 for 1993–2009. When SS rates are compared with the global trends (2.0, 2.4 and 2.8 mm yr−1, respectively) derived from tide gauge measurements for the same periods, they are smaller in the earlier era and considerably larger in the recent one. Taking into account the first estimate of land subsidence rate, 1–1.5 mm yr−1 in Singapore, the recent trend of absolute sea level rise in SS follows regional tendency.

Co-existence of wind seas and swells along the west coast of India during non-monsoon season - R. Rashmi; V. M. Aboobacker; P. Vethamony; M. P. John (Ocean Science, 2013)
An attempt has been made to understand the co-existence of wind seas and swells along the west coast of India during non-monsoon season. Wave data were collected in different years during non-monsoon season (off Goa during May 2005, off Ratnagiri during January–February 2008 and off Dwarka during December 2007–January 2008), which is fairly a calm weather season along these regions. Diurnal variation in wave parameters is noticeable along the central west coast of India (off Goa and Ratnagiri), which is due to the interaction of multidirectional waves (both wind seas and swells) of varying magnitudes and frequencies. Swells are predominantly mature (91%) and old (88%) during late pre-monsoon and post-monsoon seasons, respectively. Sea Swell Energy Ratio quantifies wind sea, swell and mixed seas prevailing in these regions during non-monsoon season. Intermodal distance (ID) between the energy peaks is moderately separated during non-monsoon season, whereas, during the shamal events, energy peaks are very close to each other (ID ∼ 0). However, pure wind seas (ID ∼ 1) are weakly present and found to co-exist with the swells almost all the time during non-monsoon season. Wind sea growth has been found while the swell propagates opposite to the direction of the wind and wind sea. Wind seas have minimum angular spreads in multimodal state. Under low winds, the interaction between wind sea and swell dominates and thereby the multimodal state reduces to unimodal state. The fetch available for the evolution of the wind sea spectrum has been estimated, and it is found to be less than 150 km. For the fetch limited condition, a non-dimensional empirical relation has been derived relating the significant wind sea height in terms of wind speed and peak wind sea period, and this relation fits for the west coast of India.

Towards an integrated forecasting system for fisheries on habitat-bound stocks - A. Christensen; M. Butenschön; Z. Gürkan; I. J. Allen (Ocean Science, 2013)
First results of a coupled modelling and forecasting system for fisheries on habitat-bound stocks are being presented. The system consists currently of three mathematically, fundamentally different model subsystems coupled offline: POLCOMS providing the physical environment implemented in the domain of the north-west European shelf, the SPAM model which describes sandeel stocks in the North Sea, and the third component, the SLAM model, which connects POLCOMS and SPAM by computing the physical–biological interaction. Our major experience by the coupling model subsystems is that well-defined and generic model interfaces are very important for a successful and extendable coupled model framework. The integrated approach, simulating ecosystem dynamics from physics to fish, allows for analysis of the pathways in the ecosystem to investigate the propagation of changes in the ocean climate and to quantify the impacts on the higher trophic level, in this case the sandeel population, demonstrated here on the basis of hindcast data. The coupled forecasting system is tested for some typical scientific questions appearing in spatial fish stock management and marine spatial planning, including determination of local and basin-scale maximum sustainable yield, stock connectivity and source/sink structure. Our presented simulations indicate that sandeel stocks are currently exploited close to the maximum sustainable yield, even though periodic overfishing seems to have occurred, but large uncertainty is associated with determining stock maximum sustainable yield due to stock inherent dynamics and climatic variability. Our statistical ensemble simulations indicates that the predictive horizon set by climate interannual variability is 2–6 yr, after which only an asymptotic probability distribution of stock properties, like biomass, are predictable.

Manifestation of two meddies in altimetry and sea-surface temperature - I. Bashmachnikov; D. Boutov; J. Dias (Ocean Science, 2013)
Two meddies were identified in the Iberian Basin using shipboard ADCP (Meddy 1) and Argo float (Meddy 2) in contrasting background conditions. Meddy 1 was observed while interacting with the Azores Current (AzC), while Meddy 2 was observed in a much calmer dynamical background, north from the AzC jet. In both cases the meddies formed a clear anticyclonic surface signal, detectable in altimetry as well as in sea-surface temperature (SST). Analysis of the in situ observations of the dynamic signal over Meddy 1 showed that the signal, generated by the moving meddy, dominated the AzC dynamics at least up to the base of the seasonal thermocline even at the late stages of its interaction with the jet. The centre of rotation of the surface signal was shifted south-westward from the axis of the meddy by about 18 km, and its dynamic radius was 2 times bigger than that of the meddy. In the centre of the anticyclonic surface signals of both meddies, SST was colder than that of the surrounding water, in contrast to warm SST anomalies in the cores of surface anticyclones generated by meandering surface currents. The latter difference gives ground for identification of meddies (as well as other sub-surface anticyclones) in comparatively dynamically calm regions using coupled altimetry–SST remote sensing data. An identification of Meddy 1 prior to the shipboard ADCP measurements was the first successful experience. At the same time, SST anomalies over the meddies were rather weak, often unstable and statistically significant only over periods of months.

Arctic rapid sea ice loss events in regional coupled climate scenario experiments - R. Döscher; T. Koenigk (Ocean Science, 2013)
Rapid sea ice loss events (RILEs) in a mini-ensemble of regional Arctic coupled climate model scenario experiments are analyzed. Mechanisms of sudden ice loss are strongly related to atmospheric circulation conditions and preconditioning by sea ice thinning during the seasons and years before the event. Clustering of events in time suggests a strong control by large-scale atmospheric circulation. Anomalous atmospheric circulation is providing warm air anomalies of up to 5 K and is forcing ice flow, affecting winter ice growth. Even without a seasonal preconditioning during winter, ice drop events can be initiated by anomalous inflow of warm air during summer. It is shown that RILEs can be generated based on atmospheric circulation changes as a major driving force without major competing mechanisms, other than occasional longwave effects during spring and summer. Other anomalous seasonal radiative forcing or short-lived forcers (e.g., soot) play minor roles or no role at all in our model. RILEs initiated by ocean forcing do not occur in the model, although cannot be ruled out due to model limitations. Mechanisms found are qualitatively in line with observations of the 2007 RILE.

Global surface-ocean p CO 2 and sea–air CO 2 flux variability from an observation-driven ocean mixed-layer scheme - C. Rödenbeck; R. F. Keeling; D. C. E. Bakker; N. Metzl; A. Olsen; C. Sabine; M. Heimann (Ocean Science, 2013)
A temporally and spatially resolved estimate of the global surface-ocean CO2 partial pressure field and the sea–air CO2 flux is presented, obtained by fitting a simple data-driven diagnostic model of ocean mixed-layer biogeochemistry to surface-ocean CO2 partial pressure data from the SOCAT v1.5 database. Results include seasonal, interannual, and short-term (daily) variations. In most regions, estimated seasonality is well constrained from the data, and compares well to the widely used monthly climatology by Takahashi et al. (2009). Comparison to independent data tentatively supports the slightly higher seasonal variations in our estimates in some areas. We also fitted the diagnostic model to atmospheric CO2 data. The results of this are less robust, but in those areas where atmospheric signals are not strongly influenced by land flux variability, their seasonality is nevertheless consistent with the results based on surface-ocean data. From a comparison with an independent seasonal climatology of surface-ocean nutrient concentration, the diagnostic model is shown to capture relevant surface-ocean biogeochemical processes reasonably well. Estimated interannual variations will be presented and discussed in a companion paper.

Sea surface freshening inferred from SMOS and ARGO salinity: impact of rain - J. Boutin; N. Martin; G. Reverdin; X. Yin; F. Gaillard (Ocean Science, 2013)
The sea surface salinity (SSS) measured from space by the Soil Moisture and Ocean Salinity (SMOS) mission has recently been revisited by the European Space Agency first campaign reprocessing. We show that, with respect to the previous version, biases close to land and ice greatly decrease. The accuracy of SMOS SSS averaged over 10 days, 100 × 100 km2 in the open ocean and estimated by comparison to ARGO (Array for Real-Time Geostrophic Oceanography) SSS is on the order of 0.3–0.4 in tropical and subtropical regions and 0.5 in a cold region. The averaged negative SSS bias (−0.1) observed in the tropical Pacific Ocean between 5° N and 15° N, relatively to other regions, is suppressed when SMOS observations concomitant with rain events, as detected from SSM/Is (Special Sensor Microwave Imager) rain rates, are removed from the SMOS–ARGO comparisons. The SMOS freshening is linearly correlated to SSM/Is rain rate with a slope estimated to −0.14 mm−1 h, after correction for rain atmospheric contribution. This tendency is the signature of the temporal SSS variability between the time of SMOS and ARGO measurements linked to rain variability and of the vertical salinity stratification between the first centimeter of the sea surface layer sampled by SMOS and the 5 m depth sampled by ARGO. However, given that the whole set of collocations includes situations with ARGO measurements concomitant with rain events collocated with SMOS measurements under no rain, the mean −0.1 bias and the negative skewness of the statistical distribution of SMOS minus ARGO SSS difference are very likely the mean signature of the vertical salinity stratification. In the future, the analysis of ongoing in situ salinity measurements in the top 50 cm of the sea surface and of Aquarius satellite SSS are expected to provide complementary information about the sea surface salinity stratification.

A clustering analysis of eddies' spatial distribution in the South China Sea - J. Yi; Y. Du; X. Wang; Z. He; C. Zhou (Ocean Science, 2013)
Spatial variation is important for studying the mesoscale eddies in the South China Sea (SCS). To investigate such spatial variations, this study made a clustering analysis on eddies' distribution using the K-means approach. Results showed that clustering tendency of anticyclonic eddies (AEs) and cyclonic eddies (CEs) were weak but not random, and the number of clusters were proved greater than four. Finer clustering results showed 10 regions where AEs densely populated and 6 regions for CEs in the SCS. Previous studies confirmed these partitions and possible generation mechanisms were related. Comparisons between AEs and CEs revealed that patterns of AE are relatively more aggregated than those of CE, and specific distinctions were summarized: (1) to the southwest of Luzon Island, AEs and CEs are generated spatially apart; AEs are likely located north of 14° N and closer to shore, while CEs are to the south and further offshore. (2) The central SCS and Nansha Trough are mostly dominated by AEs. (3) Along 112° E, clusters of AEs and CEs are located sequentially apart, and the pairs off Vietnam represent the dipole structures. (4) To the southwest of the Dongsha Islands, AEs are concentrated to the east of CEs. Overlaps of AEs and CEs in the northeastern and southern SCS were further examined considering seasonal variations. The northeastern overlap represented near-concentric distributions while the southern one was a mixed effect of seasonal variations, complex circulations and topography influences.

Observations of water masses and circulation with focus on the Eurasian Basin of the Arctic Ocean from the 1990s to the late 2000s - B. Rudels; U. Schauer; G. Björk; M. Korhonen; S. Pisarev; B. Rabe; A. Wisotzki (Ocean Science, 2013)
The circulation and water mass properties in the Eurasian Basin are discussed based on a review of previous research and an examination of observations made in recent years within, or parallel to, DAMOCLES (Developing Arctic Modeling and Observational Capabilities for Long-term Environmental Studies). The discussion is strongly biased towards observations made from icebreakers and particularly from the cruise with R/V Polarstern 2007 during the International Polar Year (IPY). Focus is on the Barents Sea inflow branch and its mixing with the Fram Strait inflow branch. It is proposed that the Barents Sea branch contributes not just intermediate water but also most of the water to the Atlantic layer in the Amundsen Basin and also in the Makarov and Canada basins. Only occasionally would high temperature pulses originating from the Fram Strait branch penetrate along the Laptev Sea slope across the Gakkel Ridge into the Amundsen Basin. Interactions between the Barents Sea and the Fram Strait branches lead to formation of intrusive layers, in the Atlantic layer and in the intermediate waters. The intrusion characteristics found downstream, north of the Laptev Sea are similar to those observed in the northern Nansen Basin and over the Gakkel Ridge, suggesting a flow from the Laptev Sea towards Fram Strait. The formation mechanisms for the intrusions at the continental slope, or in the interior of the basins if they are reformed there, have not been identified. The temperature of the deep water of the Eurasian Basin has increased in the last 10 yr rather more than expected from geothermal heating. That geothermal heating does influence the deep water column was obvious from 2007 Polarstern observations made close to a hydrothermal vent in the Gakkel Ridge, where the temperature minimum usually found above the 600–800 m thick homogenous bottom layer was absent. However, heat entrained from the Atlantic water into descending, saline boundary plumes may also contribute to the warming of the deeper layers.

Near-surface measurements of sea spray aerosol production over whitecaps in the open ocean - S. J. Norris; I. M. Brooks; B. I. Moat; M. J. Yelland; G. de Leeuw; R. W. Pascal; B. Brooks (Ocean Science, 2013)
Simultaneous measurements of near-surface aerosol (0.12

X-band COSMO-SkyMed wind field retrieval, with application to coastal circulation modeling - A. Montuori; P. de Ruggiero; M. Migliaccio; S. Pierini; G. Spezie (Ocean Science, 2013)
In this paper, X-band COSMO-SkyMed© synthetic aperture radar (SAR) wind field retrieval is investigated, and the obtained data are used to force a coastal ocean circulation model. The SAR data set consists of 60 X-band Level 1B Multi-Look Ground Detected ScanSAR Huge Region COSMO-SkyMed© SAR data, gathered in the southern Tyrrhenian Sea during the summer and winter seasons of 2010. The SAR-based wind vector field estimation is accomplished by resolving both the SAR-based wind speed and wind direction retrieval problems independently. The sea surface wind speed is retrieved by means of a SAR wind speed algorithm based on the azimuth cut-off procedure, while the sea surface wind direction is provided by means of a SAR wind direction algorithm based on the discrete wavelet transform multi-resolution analysis. The obtained wind fields are compared with ground truth data provided by both ASCAT scatterometer and ECMWF model wind fields. SAR-derived wind vector fields and ECMWF model wind data are used to construct a blended wind product regularly sampled in both space and time, which is then used to force a coastal circulation model of a southern Tyrrhenian coastal area to simulate wind-driven circulation processes. The modeling results show that X-band COSMO-SkyMed© SAR data can be valuable in providing effective wind fields for coastal circulation modeling.

The coherence of small island sea level with the wider ocean: a model study - Joanne Williams; Chris W. Hughes (Ocean Science, 2013)
Studies comparing tide gauge measurements with sea level from nearby satellite altimetry have shown good agreement for some islands and poor agreement for others, though no explanation has been offered. Using the 1/12° OCCAM ocean model, we investigate the relationship between sea level at small, open-ocean islands and offshore sea level. For every such island or seamount in the model, we compare the shallow-water sea level with the steric and bottom pressure variability in a neighbouring ring of deep water. We find a latitude-dependent range of frequencies for which off-shore sea level is poorly correlated with island sea level. This poor coherence occurs in a spectral region for which steric signals dominate, but are unable to propagate as baroclinic Rossby waves. This mode of decoupling does not arise because of islands bathymetry, as the same decoupling is seen between deep ocean points and surrounding rings.

Liquid export of Arctic freshwater components through the Fram Strait 1998–2011 - B. Rabe; P. A. Dodd; E. Hansen; E. Falck; U. Schauer; A. Mackensen; A. Beszczynska-Möller; G. Kattner; E. J. Rohling; K. Cox (Ocean Science, 2013)
We estimated the magnitude and composition of southward liquid freshwater transports in the East Greenland Current near 79° N in the Western Fram Strait between 1998 and 2011. Previous studies have found this region to be an important pathway for liquid freshwater export from the Arctic Ocean to the Nordic Seas and the North Atlantic subpolar gyre. Our transport estimates are based on six hydrographic surveys between June and September and concurrent data from moored current meters. We combined concentrations of liquid freshwater, meteoric water (river water and precipitation), sea ice melt and brine from sea ice formation, and Pacific Water, presented in Dodd et al. (2012), with volume transport estimates from an inverse model. The average of the monthly snapshots of southward liquid freshwater transports between 10.6° W and 4° E is 100 ± 23 mSv (3160 ± 730 km3 yr−1), relative to a salinity of 34.9. This liquid freshwater transport consists of about 130% water from rivers and precipitation (meteoric water), 30% freshwater from the Pacific, and −60% (freshwater deficit) due to a mixture of sea ice melt and brine from sea ice formation. Pacific Water transports showed the highest variation in time, effectively vanishing in some of the surveys. Comparison of our results to the literature indicates that this was due to atmospherically driven variability in the advection of Pacific Water along different pathways through the Arctic Ocean. Variations in most liquid freshwater component transports appear to have been most strongly influenced by changes in the advection of these water masses to the Fram Strait. However, the local dynamics represented by the volume transports influenced the liquid freshwater component transports in individual years, in particular those of sea ice melt and brine from sea ice formation. Our results show a similar ratio of the transports of meteoric water and net sea ice melt as previous studies. However, we observed a significant increase in this ratio between the surveys in 1998 and in 2009. This can be attributed to higher concentrations of sea ice melt in 2009 that may have been due to enhanced advection of freshwater from the Beaufort Gyre to the Fram Strait. Known trends and variability in the Arctic liquid freshwater inflow from rivers are not likely to have had a significant influence on the variation of liquid freshwater component transports between our surveys. On the other hand, known freshwater inflow variability from the Pacific could have caused some of the variation we observed in the Fram Strait. The apparent absence of a trend in southward liquid freshwater transports through the Fram Strait and recent evidence of an increase in liquid freshwater storage in the Arctic Ocean raise the question: how fast will the accumulated liquid freshwater be exported from the Arctic Ocean to the deep water formation regions in the North Atlantic and will an increased export occur through the Fram Strait.

Salinity in the Sicily Channel corroborates the role of the Adriatic–Ionian Bimodal Oscillating System (BiOS) in shaping the decadal variability of the Mediterranean overturning circulation - M. Gačić; K. Schroeder; G. Civitarese; S. Cosoli; A. Vetrano; G. L. Eusebi Borzelli (Ocean Science, 2013)
Previous studies have demonstrated that the salinity in the Levantine basin depends on the intensity of the Atlantic water (AW) inflow. Moreover, its spreading eastward (to the Levantine basin) or northward (to the Ionian Sea) is determined by the Ionian circulation pattern, i.e. by the Adriatic–Ionian Bimodal Oscillating System (BiOS) mechanism. The aim of this paper is to relate salinity variations in the Levantine basin to the salt content variability in the core of the Levantine Intermediate Water (LIW) passing through the Sicily Channel (SC) and its possible impact on the Western Mediterranean Transition – WMT (i.e. the sudden salinity and temperature increase in the deep layer of the Algero-Provençal subbasin occurring since 2004). From the historical data set MEDAR/MEDATLAS in the Levantine and northern Ionian, we present evidence of decadal occurrences of extreme salinities associated with the varying influx of AW over the last 60 yr. Furthermore, we show that the salinity variations in the two subbasins are out of phase. High-salinity episodes in the Levantine are a pre-conditioning for the potential occurrence of the events like the Eastern Mediterranean Transient (EMT). Cross-correlation between the salinity time series in the Levantine basin and in the SC suggests that the travel time of the LIW is between 10 and 13 yr. Comparing the timing of the salinity increase associated with the WMT and the salinity in the LIW core in the SC, we estimate that the total time interval needed for the signal propagating from the Levantine to reach the deep mixed layers of the Algero-Provençal subbasin is about 25 yr. We also showed that the extra salt input from the eastern Mediterranean contribute up to about 60% to the salt content increase in the bottom layer of the western Mediterranean.

Evaluation of global monitoring and forecasting systems at Mercator Océan - J.-M. Lellouche; O. Le Galloudec; M. Drévillon; C. Régnier; E. Greiner; G. Garric; N. Ferry; C. Desportes; C.-E. Testut; C. Bricaud; R. Bourdallé-Badie; B. Tranchant; et al. (Ocean Science, 2013)
Since December 2010, the MyOcean global analysis and forecasting system has consisted of the Mercator Océan NEMO global 1/4° configuration with a 1/12° nested model over the Atlantic and the Mediterranean. The open boundary data for the nested configuration come from the global 1/4° configuration at 20° S and 80° N. The data are assimilated by means of a reduced-order Kalman filter with a 3-D multivariate modal decomposition of the forecast error. It includes an adaptive-error estimate and a localization algorithm. A 3-D-Var scheme provides a correction for the slowly evolving large-scale biases in temperature and salinity. Altimeter data, satellite sea surface temperature and in situ temperature and salinity vertical profiles are jointly assimilated to estimate the initial conditions for numerical ocean forecasting. In addition to the quality control performed by data producers, the system carries out a proper quality control on temperature and salinity vertical profiles in order to minimise the risk of erroneous observed profiles being assimilated in the model. This paper describes the recent systems used by Mercator Océan and the validation procedure applied to current MyOcean systems as well as systems under development. The paper shows how refinements or adjustments to the system during the validation procedure affect its quality. Additionally, we show that quality checks (in situ, drifters) and data sources (satellite sea surface temperature) have as great an impact as the system design (model physics and assimilation parameters). The results of the scientific assessment are illustrated with diagnostics over the year 2010 mainly, assorted with time series over the 2007–2011 period. The validation procedure demonstrates the accuracy of MyOcean global products, whose quality is stable over time. All monitoring systems are close to altimetric observations with a forecast RMS difference of 7 cm. The update of the mean dynamic topography corrects local biases in the Indonesian Throughflow and in the western tropical Pacific. This improves also the subsurface currents at the Equator. The global systems give an accurate description of water masses almost everywhere. Between 0 and 500 m, departures from in situ observations rarely exceed 1 °C and 0.2 psu. The assimilation of an improved sea surface temperature product aims to better represent the sea ice concentration and the sea ice edge. The systems under development are still suffering from a drift which can only be detected by means of a 5-yr hindcast, preventing us from upgrading them in real time. This emphasizes the need to pursue research while building future systems for MyOcean2 forecasting.

Toward a multivariate reanalysis of the North Atlantic Ocean biogeochemistry during 1998–2006 based on the assimilation of SeaWiFS chlorophyll data - C. Fontana; P. Brasseur; J.-M. Brankart (Ocean Science, 2013)
Today, the routine assimilation of satellite data into operational models of ocean circulation is mature enough to enable the production of global reanalyses describing the ocean circulation variability during the past decades. The expansion of the "reanalysis" concept from ocean physics to biogeochemistry is a timely challenge that motivates the present study. The objective of this paper is to investigate the potential benefits of assimilating satellite-estimated chlorophyll data into a basin-scale three-dimensional coupled physical–biogeochemical model of the North Atlantic. The aim is on the one hand to improve forecasts of ocean biogeochemical properties and on the other hand to define a methodology for producing data-driven climatologies based on coupled physical–biogeochemical modeling. A simplified variant of the Kalman filter is used to assimilate ocean color data during a 9-year period. In this frame, two experiments are carried out, with and without anamorphic transformations of the state vector variables. Data assimilation efficiency is assessed with respect to the assimilated data set, nitrate of the World Ocean Atlas database and a derived climatology. Along the simulation period, the non-linear assimilation scheme clearly improves the surface analysis and forecast chlorophyll concentrations, especially in the North Atlantic bloom region. Nitrate concentration forecasts are also improved thanks to the assimilation of ocean color data while this improvement is limited to the upper layer of the water column, in agreement with recent related literature. This feature is explained by the weak correlation taken into account by the assimilation between surface phytoplankton and nitrate concentrations deeper than 50 meters. The assessment of the non-linear assimilation experiments indicates that the proposed methodology provides the skeleton of an assimilative system suitable for reanalyzing the ocean biogeochemistry based on ocean color data.

Assimilation of sea-ice concentration in a global climate model – physical and statistical aspects - S. Tietsche; D. Notz; J. H. Jungclaus; J. Marotzke (Ocean Science, 2013)
We investigate the initialisation of Northern Hemisphere sea ice in the global climate model ECHAM5/MPI-OM by assimilating sea-ice concentration data. The analysis updates for concentration are given by Newtonian relaxation, and we discuss different ways of specifying the analysis updates for mean thickness. Because the conservation of mean ice thickness or actual ice thickness in the analysis updates leads to poor assimilation performance, we introduce a proportional dependence between concentration and mean thickness analysis updates. Assimilation with these proportional mean-thickness analysis updates leads to good assimilation performance for sea-ice concentration and thickness, both in identical-twin experiments and when assimilating sea-ice observations. The simulation of other Arctic surface fields in the coupled model is, however, not significantly improved by the assimilation. To understand the physical aspects of assimilation errors, we construct a simple prognostic model of the sea-ice thermodynamics, and analyse its response to the assimilation. We find that an adjustment of mean ice thickness in the analysis update is essential to arrive at plausible state estimates. To understand the statistical aspects of assimilation errors, we study the model background error covariance between ice concentration and ice thickness. We find that the spatial structure of covariances is best represented by the proportional mean-thickness analysis updates. Both physical and statistical evidence supports the experimental finding that assimilation with proportional mean-thickness updates outperforms the other two methods considered. The method described here is very simple to implement, and gives results that are sufficiently good to be used for initialising sea ice in a global climate model for seasonal to decadal predictions.

The CORA dataset: validation and diagnostics of in-situ ocean temperature and salinity measurements - C. Cabanes; A. Grouazel; K. von Schuckmann; M. Hamon; V. Turpin; C. Coatanoan; F. Paris; S. Guinehut; C. Boone; N. Ferry; C. de Boyer Montégut; T. Carval; et al. (Ocean Science, 2013)
The French program Coriolis, as part of the French operational oceanographic system, produces the COriolis dataset for Re-Analysis (CORA) on a yearly basis. This dataset contains in-situ temperature and salinity profiles from different data types. The latest release CORA3 covers the period 1990 to 2010. Several tests have been developed to ensure a homogeneous quality control of the dataset and to meet the requirements of the physical ocean reanalysis activities (assimilation and validation). Improved tests include some simple tests based on comparison with climatology and a model background check based on a global ocean reanalysis. Visual quality control is performed on all suspicious temperature and salinity profiles identified by the tests, and quality flags are modified in the dataset if necessary. In addition, improved diagnostic tools have been developed – including global ocean indicators – which give information on the quality of the CORA3 dataset and its potential applications. CORA3 is available on request through the MyOcean Service Desk (http://www.myocean.eu/).

A global algorithm for estimating Absolute Salinity - T. J. McDougall; D. R. Jackett; F. J. Millero; R. Pawlowicz; P. M. Barker (Ocean Science, 2012)
The International Thermodynamic Equation of Seawater – 2010 has defined the thermodynamic properties of seawater in terms of a new salinity variable, Absolute Salinity, which takes into account the spatial variation of the composition of seawater. Absolute Salinity more accurately reflects the effects of the dissolved material in seawater on the thermodynamic properties (particularly density) than does Practical Salinity. When a seawater sample has standard composition (i.e. the ratios of the constituents of sea salt are the same as those of surface water of the North Atlantic), Practical Salinity can be used to accurately evaluate the thermodynamic properties of seawater. When seawater is not of standard composition, Practical Salinity alone is not sufficient and the Absolute Salinity Anomaly needs to be estimated; this anomaly is as large as 0.025 g kg−1 in the northernmost North Pacific. Here we provide an algorithm for estimating Absolute Salinity Anomaly for any location (x, y, p) in the world ocean. To develop this algorithm, we used the Absolute Salinity Anomaly that is found by comparing the density calculated from Practical Salinity to the density measured in the laboratory. These estimates of Absolute Salinity Anomaly however are limited to the number of available observations (namely 811). In order to provide a practical method that can be used at any location in the world ocean, we take advantage of approximate relationships between Absolute Salinity Anomaly and silicate concentrations (which are available globally).

An Overview of the Synoptic Antarctic Shelf-Slope Interactions (SASSI) project for the International Polar Year - K. J. Heywood; R. Muench; G. Williams (Ocean Science, 2012)

An Overview of the Synoptic Antarctic Shelf-Slope Interactions (SASSI) project for the International Polar Year - K. J. Heywood; R. Muench; G. Williams (Ocean Science, 2012)

Predictions for oil slicks detected from satellite images using MyOcean forecasting data - G. Zodiatis; R. Lardner; D. Solovyov; X. Panayidou; M. De Dominicis (Ocean Science, 2012)
The pan-European capacity for the Ocean Monitoring and Forecasting (MyOcean) Marine Core Service, implementing the Global Monitoring for Environment and Security (GMES) objectives, targets the provision of ocean state observations from various platforms and analysis and forecasting products to assist, among other downscaling activities, the needs of the operational response to marine safety, particularly concerning oil spills. The MEDSLIK oil spill and trajectory prediction system makes use of the MyOcean regional and Cyprus Coastal Ocean Forecasting and Observing System (CYCOFOS) downscaled forecasting products for operational application in the Mediterranean and pre-operational use in the Black Sea. Advanced Synthetic Aperture Radar (ASAR) satellite remote sensing images from the European Space Agency (ESA) and European Maritime Safety Agency CleanSeaNet (EMSA-CSN) provide the means for routine monitoring of the southern European seas for the detection of illegal oil discharges. MEDSLIK offers various ways, to be described in this paper, of coupling the MyOcean forecasting data with ASAR images to provide both forecasts and hindcasts for such remotely observed oil slicks. The main concern will be the drift of the oil slick and also, in the case of the forecast mode, its diffusive spreading, although some attempt is also made to estimate the changes in the state of the oil. The successful link of the satellite-detected oil slicks with their operational predictions using the MyOcean products contributes to the operational response chain and the strengthening of maritime safety for accidental or illegal spills, in implementation of the Mediterranean Decision Support System for Marine Safety (MEDESS-4MS) regarding oil spills.

Seawater capacitance – a promising proxy for mapping and characterizing drifting hydrocarbon plumes in the deep ocean - J. C. Wynn; J. A. Fleming (Ocean Science, 2012)
Hydrocarbons released into the deep ocean are an inevitable consequence of natural seep, seafloor drilling, and leaking wellhead-to-collection-point pipelines. The Macondo 252 (Deepwater Horizon) well blowout of 2010 was even larger than the Ixtoc event in the Gulf of Campeche in 1979. History suggests it will not be the last accidental release, as deepwater drilling expands to meet an ever-growing demand. For those who must respond to this kind of disaster, the first line of action should be to know what is going on. This includes knowing where an oil plume is at any given time, where and how fast it is moving, and how it is evolving or degrading. We have experimented in the laboratory with induced polarization as a method to track hydrocarbons in the seawater column and find that finely dispersed oil in seawater gives rise to a large distributed capacitance. From previous sea trials, we infer this could potentially be used to both map and characterize oil plumes, down to a ratio of less than 0.001 oil-to-seawater, drifting and evolving in the deep ocean. A side benefit demonstrated in some earlier sea trials is that this same approach in modified form can also map certain heavy placer minerals, as well as communication cables, pipelines, and wrecks buried beneath the seafloor.

Imbalance of energy and momentum source terms of the sea wave transfer equation for fully developed seas - G. V. Caudal (Ocean Science, 2012)
In the concept of full development, the sea wave spectrum is regarded as a nearly stationary solution of the wave transfer equation, where source and sink terms should be in balance with respect to both energy and momentum. Using a two-dimensional empirical sea wave spectral model at full development, this paper performs an assessment of the compatibility of the energy and momentum budgets of sea waves over the whole spectral range. Among the various combinations of model functions for wave breaking and wind source terms tested, not one is found to fulfill simultaneously the energy and momentum balance of the transfer equation. Based on experimental and theoretical grounds, wave breaking is known to contribute to frequency downshift of a narrow-banded wave spectrum when the modulational instability is combined with wave breaking. On those grounds, it is assumed that, in addition to dissipation, wave breaking produces a spectral energy flux directed toward low wavenumbers. I show that it is then possible to remove the energy and momentum budget inconsistency, and correspondingly the required strength of this spectral flux is estimated. Introducing such a downward spectral flux permits fulfilling both energy and momentum balance conditions. Meanwhile, the consistency between the transfer equation and empirical spectra, estimated by means of a cost function K, is either improved or slightly reduced, depending upon the wave breaking and wind source terms chosen. Other tests are performed in which it is further assumed that wave breaking would also be associated with azimuthal diffusion of the spectral energy. This would correspondingly reduce the required downward spectral flux by a factor of up to 5, although it would not be able to remove it entirely.

Impact of the sea surface temperature forcing on hindcasts of Madden-Julian Oscillation events using the ECMWF model - E. de Boisséson; M. A. Balmaseda; F. Vitart; K. Mogensen (Ocean Science, 2012)
This paper explores the sensitivity of hindcasts of the Madden Julian Oscillation (MJO) to the use of different sea surface temperture (SST) products as lower boundary conditions in the European Centre for Medium-range Weather Forecasts (ECMWF) atmospheric model. Three sets of monthly hindcast experiments are conducted, starting from initial conditions from the ERA interim reanalysis. First, as a reference, the atmosphere is forced by the SST used to produce ERA interim. In the second and third experiments, the SST is switched to the OSTIA (Operational Sea Surface Temperature and Sea-Ice Analysis) and the AVHRR-only (Advanced Very High Resolution Radiometer) reanalyses, respectively. Tests on the temporal resolution of the SST show that monthly fields are not optimal, while weekly and daily resolutions provide similar MJO scores. When using either OSTIA or AVHRR, the propagation of the MJO is degraded and the resulting scores are lower than in the reference experiment. Further experiments show that this loss of skill cannot be attributed to either the difference in mean state or temporal variability between the SST products. Additional diagnostics show that the phase relationship between either OSTIA or AVHRR SST and the MJO convection is distorted with respect to satellite observations and the ERA interim reanalysis. This distortion is expected to impact the MJO hindcasts, leading to a relative loss of forecast skill. A realistic representation of ocean–atmosphere interactions is thus needed for MJO hindcasts, but not all SST products – though accurate for other purposes – fulfill this requirement.

Improvement to the PhytoDOAS method for identification of coccolithophores using hyper-spectral satellite data - A. Sadeghi; T. Dinter; M. Vountas; B. B. Taylor; M. Altenburg-Soppa; I. Peeken; A. Bracher (Ocean Science, 2012)
The goal of this study was to improve PhytoDOAS, which is a new retrieval method for quantitative identification of major phytoplankton functional types (PFTs) using hyper-spectral satellite data. PhytoDOAS is an extension of the Differential Optical Absorption Spectroscopy (DOAS, a method for detection of atmospheric trace gases), developed for remote identification of oceanic phytoplankton groups. Thus far, PhytoDOAS has been successfully exploited to identify cyanobacteria and diatoms over the global ocean from SCIAMACHY (SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY) hyper-spectral data. This study aimed to improve PhytoDOAS for remote identification of coccolithophores, another functional group of phytoplankton. The main challenge for retrieving more PFTs by PhytoDOAS is to overcome the correlation effects between different PFT absorption spectra. Different PFTs are composed of different types and amounts of pigments, but also have pigments in common, e.g. chl a, causing correlation effects in the usual performance of the PhytoDOAS retrieval. Two ideas have been implemented to improve PhytoDOAS for the PFT retrieval of more phytoplankton groups. Firstly, using the fourth-derivative spectroscopy, the peak positions of the main pigment components in each absorption spectrum have been derived. After comparing the corresponding results of major PFTs, the optimized fit-window for the PhytoDOAS retrieval of each PFT was determined. Secondly, based on the results from derivative spectroscopy, a simultaneous fit of PhytoDOAS has been proposed and tested for a selected set of PFTs (coccolithophores, diatoms and dinoflagellates) within an optimized fit-window, proven by spectral orthogonality tests. The method was then applied to the processing of SCIAMACHY data over the year 2005. Comparisons of the PhytoDOAS coccolithophore retrievals in 2005 with other coccolithophore-related data showed similar patterns in their seasonal distributions, especially in the North Atlantic and the Arctic Sea. The seasonal patterns of the PhytoDOAS coccolithophores indicated very good agreement with the coccolithophore modeled data from the NASA Ocean Biochemical Model (NOBM), as well as with the global distributions of particulate inorganic carbon (PIC), provided by MODIS (MODerate resolution Imaging Spectroradiometer)-Aqua level-3 products. Moreover, regarding the fact that coccolithophores belong to the group of haptophytes, the PhytoDOAS seasonal coccolithophores showed good agreement with the global distribution of haptophytes, derived from synoptic pigment relationships applied to SeaWiFS chl a. As a case study, the simultaneous mode of PhytoDOAS has been applied to SCIAMACHY data for detecting a coccolithophore bloom which was consistent with the MODIS RGB image and the MODIS PIC map of the bloom, indicating the functionality of the method also in short-term retrievals.

Deep Western Boundary Current transport variability in the South Atlantic: preliminary results from a pilot array at 34.5° S - C. S. Meinen; A. R. Piola; R. C. Perez; S. L. Garzoli (Ocean Science, 2012)
The first direct estimates of the temporal variability of the absolute transport in the Deep Western Boundary Current (DWBC) at 34.5° S in the South Atlantic Ocean are obtained using just under one year of data from a line of four pressure-equipped inverted echo sounders. Hydrographic sections collected in 2009 and 2010 confirm, based on neutral density, temperature, salinity, and oxygen values, the presence of the DWBC, one of the main deep pathways of the Meridional Overturning Circulation. Both data sets indicate that the DWBC reconstitutes itself after breaking into eddies in the western sub-tropical Atlantic near 8° S. The amplitude and spectral character of the DWBC transport variability are comparable with those observed in the North Atlantic, where longer records exist, with the DWBC at 34.5° S exhibiting a transport standard deviation of 25 Sv and variations of ∼ 40 Sv occurring within periods as short as a few days. There is little indication of an annual cycle in the DWBC transports, although the observational records are too short to be definitive. A Monte Carlo-style analysis using 27 yr of model output from the same location as the observations indicates that about 48–60 months of data will be required to fully assess the deep transport variability. The model suggests the presence of an annual cycle in DWBC transport, however its statistical significance with even 27 yr of model output is low, suggesting that seasonal variations in the model are weak.

Effect of variable winds on current structure and Reynolds stresses in a tidal flow: analysis of experimental data in the eastern English Channel - K. A. Korotenko; A. V. Sentchev; F. G. Schmitt (Ocean Science, 2012)
Wind and wave effects on tidal current structure and turbulence throughout the water column are examined using an upward-looking acoustic Doppler current profiler (ADCP). The instrument has been deployed on the seafloor of 18-m mean depth, off the north-eastern French coast in the eastern English Channel, over 12 tidal cycles, and covered the period of the transition from mean spring to neap tide, and forcing regimes varied from calm to moderate storm conditions. During storms, we observed gusty winds with magnitudes reaching 15 m s−1 and wave heights reaching up to 1.3 m. Analysis of velocity spectra revealed a noticeable contribution of wind-induced waves to spectral structure of velocity fluctuations within the subsurface layer. Near the surface, stormy winds and waves produced a significant intensification of velocity fluctuations, particularly when the sustained wind blew against the ebb tide flow. As during wavy periods, the variance-derived Reynolds stress estimates might include a wave-induced contamination, we applied the Variance Fit method to obtain unbiased stresses and other turbulent quantities. Over calm periods, the turbulent quantities usually decreased with height above the seabed. The stresses were found to vary regularly with the predominantly semidiurnal tidal flow. The along-shore stress being generally greater during the flood flow (~2.7 Pa) than during the ebb flow (~−0.6 Pa). The turbulent kinetic energy production rate, P, and eddy viscosity, Az, followed a nearly regular cycle with close to a quarter-diurnal period. As for the stresses, near the seabed, we found the maximum values of estimated quantities of P and Az to be 0.1 Wm−3 and 0.5 m2 s−1, respectively, during the flood flow. Over the storm periods, we found the highest unbiased stress values (~−2.6 Pa) during ebb when tidal currents were opposite to the southwesterly winds while, during the flood, the surface stresses slightly exceeded those estimated for a calm period. A comparison of obtained results gives a good agreement with those of other researchers working on direct measurements of turbulence in tidal flows.

Marine atmospheric boundary layer over some Southern Ocean fronts during the IPY BGH 2008 cruise - C. Messager; S. Speich; E. Key (Ocean Science, 2012)
A set of meteorological instruments was added to an oceanographic cruise crossing the Southern Ocean from Cape Town to 57°33' S during the summer of 2008. The Cape Cauldron, the Subtropical, Subantarctic, Polar and southern Antarctic Circumpolar current fronts were successively crossed. The recorded data permitted to derive the exchange of momentum, heat and water vapour at the ocean-atmosphere interface. A set of 38 radiosonde releases complemented the dataset. The marine atmospheric boundary layer characteristics and air-sea interaction when the ship crossed the fronts and eddies are discussed. The specific role of the atmospheric synoptic systems advection on the air-sea interaction over these regions is highlighted. Additionally, the Subantarctic front mesoscale variability induced an anticyclonic eddy considered as part of the Subantarctic front. The specific influence of this Agulhas ring on the aloft atmosphere is also presented.

Validation of the NEMO-ERSEM operational ecosystem model for the North West European Continental Shelf - K. P. Edwards; R. Barciela; M. Butenschön (Ocean Science, 2012)
This paper details updates to the Met Office's operational coupled hydrodynamic-ecosystem model from the 7 km Medium-Resolution Continental Shelf – POLCOMS-ERSEM (MRCS-PE) system (Siddorn et al., 2007) to the 7 km Atlantic Margin Model NEMO-ERSEM (AMM7-NE) system. We also provide a validation of the ecosystem component of the new operational system. Comparisons have been made between the model variables and available in situ, satellite and climatological data. The AMM7-NE system has also been benchmarked against the MRCS-PE system. The transition to the new AMM7-NE system was successful and it has been running operationally since March 2012 and has been providing products through MyOcean (http://www.myocean.eu.org) since that time. The results presented herein show the AMM7-NE system performs better than the MRCS-PE system with the most improvement in the model nutrient fields. The problem of nutrient accumulation in the MRCS-PE system appears to be solved in the new AMM7-NE system with nutrient fields improved throughout the domain as discussed in Sect. 4. Improvements in model chlorophyll are also seen but are more modest.

The link between the Barents Sea and ENSO events simulated by NEMO model - V. N. Stepanov; H. Zuo; K. Haines (Ocean Science, 2012)
An analysis of observational data in the Barents Sea along a meridian at 33°30' E between 70°30' and 72°30' N has reported a negative correlation between El Niño/La Niña Southern Oscillation (ENSO) events and water temperature in the top 200 m: the temperature drops about 0.5 °C during warm ENSO events while during cold ENSO events the top 200 m layer of the Barents Sea is warmer. Results from 1 and 1/4-degree global NEMO models show a similar response for the whole Barents Sea. During the strong warm ENSO event in 1997–1998 an anomalous anticyclonic atmospheric circulation over the Barents Sea enhances heat loses, as well as substantially influencing the Barents Sea inflow from the North Atlantic, via changes in ocean currents. Under normal conditions along the Scandinavian peninsula there is a warm current entering the Barents Sea from the North Atlantic, however after the 1997–1998 event this current is weakened. During 1997–1998 the model annual mean temperature in the Barents Sea is decreased by about 0.8 °C, also resulting in a higher sea ice volume. In contrast during the cold ENSO events in 1999–2000 and 2007–2008, the model shows a lower sea ice volume, and higher annual mean temperatures in the upper layer of the Barents Sea of about 0.7 °C. An analysis of model data shows that the strength of the Atlantic inflow in the Barents Sea is the main cause of heat content variability, and is forced by changing pressure and winds in the North Atlantic. However, surface heat-exchange with the atmosphere provides the means by which the Barents sea heat budget relaxes to normal in the subsequent year after the ENSO events.

Arctic surface temperatures from Metop AVHRR compared to in situ ocean and land data - G. Dybkjær; R. Tonboe; J. L. Høyer (Ocean Science, 2012)
The ice surface temperature (IST) is an important boundary condition for both atmospheric and ocean and sea ice models and for coupled systems. An operational ice surface temperature product using satellite Metop AVHRR infra-red data was developed for MyOcean. The IST can be mapped in clear sky regions using a split window algorithm specially tuned for sea ice. Clear sky conditions prevail during spring in the Arctic, while persistent cloud cover limits data coverage during summer. The cloud covered regions are detected using the EUMETSAT cloud mask. The Metop IST compares to 2 m temperature at the Greenland ice cap Summit within STD error of 3.14 °C and to Arctic drifting buoy temperature data within STD error of 3.69 °C. A case study reveals that the in situ radiometer data versus satellite IST STD error can be much lower (0.73 °C) and that the different in situ measurements complicate the validation. Differences and variability between Metop IST and in situ data are analysed and discussed. An inter-comparison of Metop IST, numerical weather prediction temperatures and in situ observation indicates large biases between the different quantities. Because of the scarcity of conventional surface temperature or surface air temperature data in the Arctic, the satellite IST data with its relatively good coverage can potentially add valuable information to model analysis for the Arctic atmosphere.

Microstructure observations during the spring 2011 STRATIPHYT-II cruise in the northeast Atlantic - E. Jurado; H. A. Dijkstra; H. J. van der Woerd (Ocean Science, 2012)
Small-scale temperature and conductivity variations have been measured in the upper 100 m of the northeast Atlantic during the STRATIPHYT-II cruise (Las Palmas–Reykjavik, 6 April–3 May 2011). The measurements were done at midday and comprised 2 to 15 vertical profiles at each station. The derived turbulent quantities show a transition between weakly-stratified (mixed layer depth, MLD, 100), which was centered at about 48° N. The temperature eddy diffusivities, KT, range from 10−5 to 100 m2 s−1 in the weakly-stratified stations, and range from 3 × 10−4 to 2 × 100 m2 s−1 in the well-mixed stations. The turbulent kinetic energy dissipation rates, ε, range from 3 × 10−8 to 2 × 10−6 m2 s−3 south of the transition zone, and from 10−7 to 10−5 m2 s−3 north of the transition zone. The station-averaged KT values throughout the mixed layer increase exponentially with the wind speed. The station-averaged ε values throughout the mixed layer scale with the wind stress similarity variable with a scaling factor of about 1.8 in the wind-dominated stations (ε ≈ 1.8 u☆3/(−κz)). The values of KT and ε are on average 10 times higher compared to the values measured at the same stations in July 2009. The results presented here constitute a unique data set giving large spatial coverage of upper ocean spring turbulence quantities.

Surface signature of Mediterranean water eddies in the Northeastern Atlantic: effect of the upper ocean stratification - I. Bashmachnikov; X. Carton (Ocean Science, 2012)
Meddies, intra-thermocline eddies of Mediterranean water, can often be detected at the sea surface as positive sea-level anomalies. Here we study the surface signature of several meddies tracked with RAFOS floats and AVISO altimetry. While pushing its way through the water column, a meddy raises isopycnals above. As a consequence of potential vorticity conservation, negative relative vorticity is generated in the upper layer. During the initial period of meddy acceleration after meddy formation or after a stagnation stage, a cyclonic signal is also generated at the sea-surface, but mostly the anticyclonic surface signal follows the meddy. Based on geostrophy and potential vorticity balance, we present theoretical estimates of the intensity of the surface signature. It appears to be proportional to the meddy core radius and to the Coriolis parameter, and inversely proportional to the core depth and buoyancy frequency. This indicates that surface signature of a meddy may be strongly reduced by the upper ocean stratification. Using climatic distribution of the stratification intensity, we claim that the southernmost limit for detection in altimetry of small meddies (with radii on the order of 10–15 km) should lie in the subtropics (35–45° N), while large meddies (with radii of 25–30 km) could be detected as far south as the northern tropics (25–35° N). Those results agree with observations.

Characteristics of the seasonal cycle of surface layer salinity in the global ocean - F. M. Bingham; G. R. Foltz; M. J. McPhaden (Ocean Science, 2012)
The seasonal variability of surface layer salinity (SLS), evaporation (E), precipitation (P), E-P, advection and vertical entrainment over the global ocean is examined using in situ salinity data, the National Centers for Environmental Prediction's Climate System Forecast Reanalysis and a number of other ancillary data. Seasonal amplitudes and phases are calculated using harmonic analysis and presented in all areas of the open ocean between 60° S and 60° N. Areas with large amplitude SLS seasonal variations include: the intertropical convergence zone (ITCZ) in the Atlantic, Pacific and Indian Oceans; western marginal seas of the Pacific; and the Arabian Sea. The median amplitude in areas that have statistically significant seasonal cycles of SLS is 0.19. Between about 60° S and 60° N, 37% of the ocean surface has a statistically significant seasonal cycle of SLS and 75% has a seasonal cycle of E-P. Phases of SLS have a bimodal distribution, with most areas in the Northern Hemisphere peaking in SLS in March/April and in the Southern Hemisphere in September/October. The seasonal cycle is also estimated for surface freshwater forcing using a mixed-layer depth climatology. With the exception of areas near the western boundaries of the North Atlantic and North Pacific, seasonal variability is dominated by precipitation. Surface freshwater forcing also has a bimodal distribution, with peaks in January and July, 1–2 months before the peaks of SLS. Seasonal amplitudes and phases calculated for horizontal advection show it to be important in the tropical oceans. Vertical entrainment, estimated from mixed-layer heaving, is largest in mid and high latitudes, with a seasonal cycle that peaks in late winter. The amplitudes and phases of SLS and surface fluxes compare well in a qualitative sense, suggesting that much of the variability in SLS is due to E-P. However, the amplitudes of SLS are somewhat different than would be expected and the peak of SLS comes typically about one month earlier than expected. The differences of the amplitudes of the two quantities is largest in such areas as the Amazon River plume, the Arabian Sea, the ITCZ and the eastern equatorial Pacific and Atlantic.

Modelling temperature and salinity in Liverpool Bay and the Irish Sea: sensitivity to model type and surface forcing - C. K. O'Neill; J. A. Polton; J. T. Holt; E. J. O'Dea (Ocean Science, 2012)
Three shelf sea models are compared against observed surface temperature and salinity in Liverpool Bay and the Irish Sea: a 7 km NEMO (Nucleus for European Modelling of the Ocean) model, and 12 km and 1.8 km POLCOMS (Proudman Oceanographic Laboratory Coastal Ocean Modelling System) models. Each model is run with two different surface forcing datasets of different resolutions. Comparisons with a variety of observations from the Liverpool Bay Coastal Observatory show that increasing the surface forcing resolution improves the modelled surface temperature in all the models, in particular reducing the summer warm bias and winter cool bias. The response of surface salinity is more varied with improvements in some areas and deterioration in others. The 7 km NEMO model performs as well as the 1.8 km POLCOMS model when measured by overall skill scores, although the sources of error in the models are different. NEMO is too weakly stratified in Liverpool Bay, whereas POLCOMS is too strongly stratified. The horizontal salinity gradient, which is too strong in POLCOMS, is better reproduced by NEMO which uses a more diffusive horizontal advection scheme. This leads to improved semi-diurnal variability in salinity in NEMO at a mooring site located in the Liverpool Bay ROFI (region of freshwater influence) area.

Towards high resolution mapping of 3-D mesoscale dynamics from observations - B. Buongiorno Nardelli; S. Guinehut; A. Pascual; Y. Drillet; S. Ruiz; S. Mulet (Ocean Science, 2012)
The MyOcean R&D project MESCLA (MEsoSCaLe dynamical Analysis through combined model, satellite and in situ data) was devoted to the high resolution 3-D retrieval of tracer and velocity fields in the oceans, based on the combination of in situ and satellite observations and quasi-geostrophic dynamical models. The retrieval techniques were also tested and compared with the output of a primitive equation model, with particular attention to the accuracy of the vertical velocity field as estimated through the Q vector formulation of the omega equation. The project focused on a test case, covering the region where the Gulf Stream separates from the US East Coast. This work demonstrated that innovative methods for the high resolution mapping of 3-D mesoscale dynamics from observations can be used to build the next generations of operational observation-based products.

The Mediterranean Ocean Colour Observing System – system development and product validation - G. Volpe; S. Colella; V. Forneris; C. Tronconi; R. Santoleri (Ocean Science, 2012)
This paper presents the Mediterranean Ocean Colour Observing System in the framework of the growing demand of near real-time data emerging within the operational oceanography international context. The main issues related to the satellite operational oceanography are tied to the following: (1) the near real-time ability to track data flow uncertainty sources; (2) in case of failure, to provide backup solutions to end-users; and (3) to scientifically assess the product quality. We describe the major scientific and technological steps made to develop, maintain and improve the operational system and its products. A method for assessing the near real-time product quality is developed and its limitation discussed. Main results are concerned with the degradation, starting from mid-2010, of the MODIS Aqua channel at 443 nm with its successive recovery thanks to the new calibration scheme implemented in the recently released SeaDAS version 6.4. The product validation analysis highlights that SeaWiFS chlorophyll product over the Mediterranean Sea is the best performing in comparison with those of MODIS and MERIS. Despite their general good agreement with in situ observations, MODIS- and MERIS-derived chlorophyll present a slight and systematic underestimation of the in situ counter part. The most relevant implications induced by these results are discussed from an operational point of view.

Evaluation of Release-05 GRACE time-variable gravity coefficients over the ocean - D. P. Chambers; J. A. Bonin (Ocean Science, 2012)
The latest release of GRACE (Gravity Recovery and Climate Experiment) gravity field coefficients (Release-05, or RL05) are evaluated for ocean applications. Data have been processed using the current methodology for Release-04 (RL04) coefficients, and have been compared to output from two different ocean models. Results indicate that RL05 data from the three Science Data Centers – the Center for Space Research (CSR), GeoForschungsZentrum (GFZ), and Jet Propulsion Laboratory (JPL) – are more consistent among themselves than the previous RL04 data. Moreover, the variance of residuals with the output of an ocean model is 50–60% lower for RL05 data than for RL04 data. A more optimized destriping algorithm is also tested, which improves the results slightly. By comparing the GRACE maps with two different ocean models, we can better estimate the uncertainty in the RL05 maps. We find the standard error to be about 1 cm (equivalent water thickness) in the low- and mid-latitudes, and between 1.5 and 2 cm in the polar and subpolar oceans, which is comparable to estimated uncertainty for the output from the ocean models.

High resolution 3-D temperature and salinity fields derived from in situ and satellite observations - S. Guinehut; A.-L. Dhomps; G. Larnicol; P.-Y. Le Traon (Ocean Science, 2012)
This paper describes an observation-based approach that efficiently combines the main components of the global ocean observing system using statistical methods. Accurate but sparse in situ temperature and salinity profiles (mainly from Argo for the last 10 yr) are merged with the lower accuracy but high-resolution synthetic data derived from satellite altimeter and sea surface temperature observations to provide global 3-D temperature and salinity fields at high temporal and spatial resolution. The first step of the method consists in deriving synthetic temperature fields from altimeter and sea surface temperature observations, and salinity fields from altimeter observations, through multiple/simple linear regression methods. The second step of the method consists in combining the synthetic fields with in situ temperature and salinity profiles using an optimal interpolation method. Results show the revolutionary nature of the Argo observing system. Argo observations now allow a global description of the statistical relationships that exist between surface and subsurface fields needed for step 1 of the method, and can constrain the large-scale temperature and mainly salinity fields during step 2 of the method. Compared to the use of climatological estimates, results indicate that up to 50% of the variance of the temperature fields can be reconstructed from altimeter and sea surface temperature observations and a statistical method. For salinity, only about 20 to 30% of the signal can be reconstructed from altimeter observations, making the in situ observing system essential for salinity estimates. The in situ observations (step 2 of the method) further reduce the differences between the gridded products and the observations by up to 20% for the temperature field in the mixed layer, and the main contribution is for salinity and the near surface layer with an improvement up to 30%. Compared to estimates derived using in situ observations only, the merged fields provide a better reconstruction of the high resolution temperature and salinity fields. This also holds for the large-scale and low-frequency fields thanks to a better reduction of the aliasing due to the mesoscale variability. Contribution of the merged fields is then illustrated to describe qualitatively the temperature variability patterns for the period from 1993 to 2009.

A 20-year reanalysis experiment in the Baltic Sea using three-dimensional variational (3DVAR) method - W. Fu; J. She; M. Dobrynin (Ocean Science, 2012)
A 20-year retrospective reanalysis of the ocean state in the Baltic Sea is constructed by assimilating available historical temperature and salinity profiles into an operational numerical model with three-dimensional variational (3DVAR) method. To determine the accuracy of the reanalysis, the authors present a series of comparisons to independent observations on a monthly mean basis. In the reanalysis, temperature (T) and salinity (S) fit better with independent measurements than the free run at different depths. Overall, the mean biases of temperature and salinity for the 20 year period are reduced by 0.32 °C and 0.34 psu, respectively. Similarly, the mean root mean square error (RMSE) is decreased by 0.35 °C for temperature and 0.3 psu for salinity compared to the free run. The modeled sea surface temperature, which is mainly controlled by the weather forcing, shows the least improvements due to sparse in situ observations. Deep layers, on the other hand, witness significant and stable model error improvements. In particular, the salinity related to saline water intrusions into the Baltic Proper is largely improved in the reanalysis. The major inflow events such as in 1993 and 2003 are captured more accurately as the model salinity in the bottom layer is increased by 2–3 psu. Compared to independent sea level at 14 tide gauge stations, the correlation between model and observation is increased by 2%–5%, while the RMSE is generally reduced by 10 cm. It is found that the reduction of RMSE comes mainly from the reduction of mean bias. In addition, the changes in density induced by the assimilation of T/S contribute little to the barotropic transport in the shallow Danish Transition zone. The mixed layer depth exhibits strong seasonal variations in the Baltic Sea. The basin-averaged value is about 10 m in summer and 30 m in winter. By comparison, the assimilation induces a change of 20 m to the mixed layer depth in deep waters and wintertime, whereas small changes of about 2 m occur in summer and shallow waters. It is related to the strong heating in summer and the dominant role of the surface forcing in shallow water, which largely offset the effect of the assimilation.

High frequency fluctuations in the heat content of an ocean general circulation model - A. M. Huerta-Casas; D. J. Webb (Ocean Science, 2012)
The transport and storage of heat by the ocean is of crucial importance because of its effect on ocean dynamics and its impact on the atmosphere, climate and climate change. Unfortunately, limits to the amount of data that can be collected and stored mean that many experimental and modelling studies of the heat budget have to make use of mean datasets where the effects of short term fluctuations are lost. In this paper we investigate the magnitude of the resulting errors by making use of data from OCCAM, a high resolution global ocean model. The model carries out a proper heat balance every time step so any imbalances that are found in the analysis must result from the use of mean fields. The study concentrates on two areas of the ocean affecting the El Nino. The first is the region of tropical instability waves north of the Equator. The second is in the upwelling region along the Equator. It is shown that in both cases, processes with a period of less than five days can have a significant impact on the heat budget. Thus, analyses using data averaged over five days or more are likely to have significant errors. It is also shown that if a series of instantaneous values is available, reasonable estimates can be made of the size of the errors. In model studies, such values are available in the form of the datasets used to restart the model. In experimental studies they may be in the form of individual unaveraged observations.

Obstacles and benefits of the implementation of a reduced-rank smoother with a high resolution model of the tropical Atlantic Ocean - N. Freychet; E. Cosme; P. Brasseur; J.-M. Brankart; E. Kpemlie (Ocean Science, 2012)
Most of oceanographic operational centers use three-dimensional data assimilation schemes to produce reanalyses. We investigate here the benefits of a smoother, i.e. a four-dimensional formulation of statistical assimilation. A square-root sequential smoother is implemented with a tropical Atlantic Ocean circulation model. A simple twin experiment is performed to investigate its benefits, compared to its corresponding filter. Despite model's non-linearities and the various approximations used for its implementation, the smoother leads to a better estimation of the ocean state, both on statistical (i.e. mean error level) and dynamical points of view, as expected from linear theory. Smoothed states are more in phase with the dynamics of the reference state, an aspect that is nicely illustrated with the chaotic dynamics of the North Brazil Current rings. We also show that the smoother efficiency is strongly related to the filter configuration. One of the main obstacles to implement the smoother is then to accurately estimate the error covariances of the filter. Considering this, benefits of the smoother are also investigated with a configuration close to situations that can be managed by operational center systems, where covariances matrices are fixed (optimal interpolation). We define here a simplified smoother scheme, called half-fixed basis smoother, that could be implemented with current reanalysis schemes. Its main assumption is to neglect the propagation of the error covariances matrix, what leads to strongly reduce the cost of assimilation. Results illustrate the ability of this smoother to provide a solution more consistent with the dynamics, compared to the filter. The smoother is also able to produce analyses independently of the observation frequency, so the smoothed solution appears more continuous in time, especially in case of a low frenquency observation network.

Assimilation of SLA along track observations in the Mediterranean with an oceanographic model forced by atmospheric pressure - S. Dobricic; C. Dufau; P. Oddo; N. Pinardi; I. Pujol; M.-H. Rio (Ocean Science, 2012)
A large number of SLA observations at a high along track horizontal resolution are an important ingredient of the data assimilation in the Mediterranean Forecasting System (MFS). Recently, new higher-frequency SLA products have become available, and the atmospheric pressure forcing has been implemented in the numerical model used in the MFS data assimilation system. In a set of numerical experiments, we show that, in order to obtain the most accurate analyses, the ocean model should include the atmospheric pressure forcing and the observations should contain the atmospheric pressure signal. When the model is not forced by the atmospheric pressure, the high-frequency filtering of SLA observations, however, improves the quality of the SLA analyses. It is further shown by comparing the power density spectra of the model fields and observations that the model is able to extract the correct information from noisy observations even without their filtering during the pre-processing.

Fate of river Tiber discharge investigated through numerical simulation and satellite monitoring - R. Inghilesi; L. Ottolenghi; A. Orasi; C. Pizzi; F. Bignami; R. Santoleri (Ocean Science, 2012)
The aim of this study was to determine the dispersion of passive pollutants associated with the Tiber discharge into the Tyrrhenian Sea using numerical marine dispersion models and satellite data. Numerical results obtained in the simulation of realistic discharge episodes were compared with the corresponding evolution of the spatial distributions of MODIS diffuse light attenuation coefficient at 490 nm (K490), and the results were discussed with reference to the local climate and the seasonal sub-regional circulation regime. The numerical model used for the simulation of the sub-tidal circulation was a Mediterranean sub-regional scale implementation of the Princeton Ocean Model (POM), nested in the large-scale Mediterranean Forecasting System. The nesting method enabled the model to be applied to almost every area in the Mediterranean Sea and also to be used in seasons for which imposing climatological boundary conditions would have been questionable. Dynamical effects on coastal circulation and on water density due to the Tiber discharge were additionally accounted for in the oceanographic model by implementing the river estuary as a point source of a buoyant jet. A Lagrangian particle dispersion model fed with the POM current fields was then run in order to reproduce the effect of the turbulent transport of passive tracers mixed in the plume with the coastal flow. Two significant episodes of river discharge in both winter and summer conditions were discussed in this paper. It was found that the winter regime was characterized by the presence of a strong coastal jet flowing with the ambient current. In summer the prevailing wind regime induced coastal downwelling conditions, which tended to confine the riverine waters close to the shore. In such conditions sudden wind reversals due to local weather perturbations, causing moderate local upwelling, proved to be the only effective way to disperse the tracers offshore, moving the plume from the coast and detaching large pools of freshwater.

Assimilating GlobColour ocean colour data into a pre-operational physical-biogeochemical model - D. A. Ford; K. P. Edwards; D. Lea; R. M. Barciela; M. J. Martin; J. Demaria (Ocean Science, 2012)
As part of the GlobColour project, daily chlorophyll a observations, derived using remotely sensed ocean colour data from the MERIS, MODIS and SeaWiFS sensors, are produced. The ability of these products to be assimilated into a pre-operational global coupled physical-biogeochemical model has been tested, on both a hindcast and near-real-time basis, and the impact on the system assessed. The assimilation was found to immediately and considerably improve the bias, root mean square error and correlation of modelled surface chlorophyll concentration compared to the GlobColour observations, an improvement which was sustained throughout the year and in every ocean basin. Errors against independent in situ chlorophyll observations were also reduced, both at and beneath the ocean surface. However, the model fit to in situ observations was not consistently better than that of climatology, due to errors in the underlying model. The assimilation scheme used is multivariate, updating all biogeochemical model state variables at all depths. The other variables were not degraded by the assimilation, with annual mean surface fields of nutrients, alkalinity and carbon variables remaining of similar quality compared to climatology. There was evidence of improved representation of zooplankton concentration, and reduced errors were seen against in situ observations of nitrate and pCO2, but too few observations were available to conclude about global model skill. The near-real-time GlobColour products were found to be sufficiently reliable for operational purposes, and of benefit to both operational-style systems and reanalyses.

On the shelf resonances of the Gulf of Carpentaria and the Arafura Sea - D. J. Webb (Ocean Science, 2012)
A numerical model is used to investigate the resonances of the Gulf of Carpentaria and the Arafura Sea, and the additional insights that come from extending the analysis into the complex angular velocity plane. When the model is forced at the shelf edge with physically realistic real values of the angular velocity, the response functions at points within the region show maxima and other behaviour which imply that resonances are involved but provide little additional information. The study is then extended to complex angular velocities, and the results then show a clear pattern of gravity wave and Rossby wave like resonances. The properties of the resonances are investigated and used to reinterpret the response at real values of angular velocity. It is found that in some regions the response is dominated by modes trapped between the shelf edge and the coast or between opposing coastlines. In other regions the resonances show cooperative behaviour, possibly indicating the importance of other physical processes.

An operational model for the West Iberian coast: products and services - M. Mateus; G. Riflet; P. Chambel; L. Fernandes; R. Fernandes; M. Juliano; F. Campuzano; H. de Pablo; R. Neves (Ocean Science, 2012)
This paper presents the structure and application of a regional scale operational modelling tool for the West Iberian coast, and discusses its potential for products and services for both scientific and coastal management activities. The forecasting suite includes nested hydrodynamic models forced with up-to-date meteorological forecast data and large-scale model results. The present status of the system and its recent upgrades are reviewed, offering a general description of the main components of the system: the forcing data, the circulation model, the model outputs and the validation methodology of model results. Seasonal differences in temperature, salinity and current velocity fields are illustrated and show satisfactory reproduction of the top and deep layer thermodynamics. The system provides boundary forcing for a number of local-scale model applications via downscaling of the solution and enables potential products and services from which civil society will benefit.

Net primary productivity, upwelling and coastal currents in the Gulf of Ulloa, Baja California, México - E. González-Rodríguez; A. Trasviña-Castro; G. Gaxiola-Castro; L. Zamudio; R. Cervantes-Duarte (Ocean Science, 2012)
The Gulf of Ulloa, a highly productive area off the western coast of the Baja California Peninsula, is examined for five successive years (2003–2007) by using satellite data and seasonal net primary productivity (NPP) estimates obtained from a vertical generalised production model. The results identify that northwestern winds blow parallel to the coast throughout the year. However, highest NPP occurs from March to June. During this period, an equatorward coastal current transports water from neighbouring upwelling areas to the northern Gulf of Ulloa and in combination with local upwelling, which injects nutrients into the euphotic zone, produce the observed increase in NPP. The opposite situation occurs in late summer when a warm poleward current of tropical characteristics arrives and inhibits the productivity in the whole region and generates the yearly lowest NPP levels. Our findings reveal the importance of lateral advection in the modulation of the primary productivity in this subtropical upwelling region.

Assessment of a physical-biogeochemical coupled model system for operational service in the Baltic Sea - Z. Wan; J. She; M. Maar; L. Jonasson; J. Baasch-Larsen (Ocean Science, 2012)
Thanks to the abundant observation data, we are able to deploy the traditional point-to-point comparison and statistical measures in combination with a comprehensive model validation scheme to assess the skills of the biogeochemical model ERGOM in providing an operational service for the Baltic Sea. The model assessment concludes that the operational products can resolve the main observed seasonal features for phytoplankton biomass, dissolved inorganic nitrogen, dissolved inorganic phosphorus and dissolved oxygen in euphotic layers as well as their vertical profiles. This assessment reflects that the model errors of the operational system at the current stage are mainly caused by insufficient light penetration, excessive organic particle export downward, insufficient regional adaptation and some from improper initialization. This study highlights the importance of applying multiple schemes in order to assess model skills rigidly and identify main causes for major model errors.

Three-dimensional modelling of wave-induced current from the surf zone to the inner shelf - H. Michaud; P. Marsaleix; Y. Leredde; C. Estournel; F. Bourrin; F. Lyard; C. Mayet; F. Ardhuin (Ocean Science, 2012)
We develop and implement a new method to take into account the impact of waves into the 3-D circulation model SYMPHONIE (Marsaleix et al., 2008, 2009a) following the simplified equations of Bennis et al. (2011) which use glm2z-RANS theory (Ardhuin et al., 2008c). These adiabatic equations are completed by additional parameterizations of wave breaking, bottom friction and wave-enhanced vertical mixing, making the forcing valid from the surf zone through to the open ocean. The wave forcing is performed by wave generation and propagation models WAVEWATCH III® (Tolman, 2008, 2009; Ardhuin et al., 2010) and SWAN (Booij et al., 1999). The model is tested and compared with other models for a plane beach test case, previously tested by Haas and Warner (2009)and Uchiyama et al. (2010). A comparison is also made with the laboratory measurements of Haller et al. (2002) of a barred beach with channels. Results fit with previous simulations performed by other models and with available observational data. Finally, a realistic case is simulated with energetic waves travelling over a coast of the Gulf of Lion (in the northwest of the Mediterranean Sea) for which currents are available at different depths as well as an accurate bathymetric database of the 0–10 m depth range. A grid nesting approach is used to account for the different forcings acting at different spatial scales. The simulation coupling the effects of waves and currents is successful to reproduce the powerful northward littoral drift in the 0–15 m depth zone. More precisely, two distinct cases are identified: When waves have a normal angle of incidence with the coast, they are responsible for complex circulation cells and rip currents in the surf zone, and when they travel obliquely, they generate a northward littoral drift. These features are more complicated than in the test cases, due to the complex bathymetry and the consideration of wind and non-stationary processes. Wave impacts in the inner shelf are less visible since wind and regional circulation seem to be the predominant forcings. Besides, a discrepancy between model and observations is noted at that scale, possibly linked to an underestimation of the wind stress. This three-dimensional method allows a good representation of vertical current profiles and permits the calculation of the shear stress associated with waves and currents. Future work will focus on the combination with a sediment transport model.

TOPAZ4: an ocean-sea ice data assimilation system for the North Atlantic and Arctic - P. Sakov; F. Counillon; L. Bertino; K. A. Lisæter; P. R. Oke; A. Korablev (Ocean Science, 2012)
We present a detailed description of TOPAZ4, the latest version of TOPAZ – a coupled ocean-sea ice data assimilation system for the North Atlantic Ocean and Arctic. It is the only operational, large-scale ocean data assimilation system that uses the ensemble Kalman filter. This means that TOPAZ features a time-evolving, state-dependent estimate of the state error covariance. Based on results from the pilot MyOcean reanalysis for 2003–2008, we demonstrate that TOPAZ4 produces a realistic estimate of the ocean circulation in the North Atlantic and the sea-ice variability in the Arctic. We find that the ensemble spread for temperature and sea-level remains fairly constant throughout the reanalysis demonstrating that the data assimilation system is robust to ensemble collapse. Moreover, the ensemble spread for ice concentration is well correlated with the actual errors. This indicates that the ensemble statistics provide reliable state-dependent error estimates – a feature that is unique to ensemble-based data assimilation systems. We demonstrate that the quality of the reanalysis changes when different sea surface temperature products are assimilated, or when in-situ profiles below the ice in the Arctic Ocean are assimilated. We find that data assimilation improves the match to independent observations compared to a free model. Improvements are particularly noticeable for ice thickness, salinity in the Arctic, and temperature in the Fram Strait, but not for transport estimates or underwater temperature. At the same time, the pilot reanalysis has revealed several flaws in the system that have degraded its performance. Finally, we show that a simple bias estimation scheme can effectively detect the seasonal or constant bias in temperature and sea-level.

The analysis of large-scale turbulence characteristics in the Indonesian seas derived from a regional model based on the Princeton Ocean Model - K. O'Driscoll; V. Kamenkovich (Ocean Science, 2012)
Turbulence characteristics in the Indonesian seas on the horizontal scale of order of 100 km were calculated with a regional model of the Indonesian seas circulation in the area based on the Princeton Ocean Model (POM). As is well known, the POM incorporates the Mellor–Yamada turbulence closure scheme. The calculated characteristics are: twice the turbulence kinetic energy per unit mass, q2; the turbulence master scale, ℓ; mixing coefficients of momentum, KM; and temperature and salinity, KH; etc. The analyzed turbulence has been generated essentially by the shear of large-scale ocean currents and by the large-scale wind turbulence. We focused on the analysis of turbulence around important topographic features, such as the Lifamatola Sill, the North Sangihe Ridge, the Dewakang Sill, and the North and South Halmahera Sea Sills. In general, the structure of turbulence characteristics in these regions turned out to be similar. For this reason, we have carried out a detailed analysis of the Lifamatola Sill region because dynamically this region is very important and some estimates of mixing coefficients in this area are available. Briefly, the main results are as follows. The distribution of q2 is quite adequately reproduced by the model. To the north of the Lifamatola Sill (in the Maluku Sea) and to the south of the Sill (in the Seram Sea), large values of q2 occur in the deep layer extending several hundred meters above the bottom. The observed increase of q2 near the very bottom is probably due to the increase of velocity shear and the corresponding shear production of q2 very close to the bottom. The turbulence master scale, ℓ, was found to be constant in the main depth of the ocean, while ℓ rapidly decreases close to the bottom, as one would expect. However, in deep profiles away from the sill, the effect of topography results in the ℓ structure being unreasonably complicated as one moves towards the bottom. Values of 15 to 20 × 10−4 m2 s−1 were obtained for KM and KH in deep water in the vicinity of the Lifamatola Sill. These estimates agree well with basin-scale averaged values of 13.3 × 10−4 m2 s−1 found diagnostically for KH in the deep Banda and Seram Seas (Gordon et al., 2003) and a value of 9.0 × 10−4 m2 s−1 found diagnostically for KH for the deep Banda Sea system (van Aken et al., 1988). The somewhat higher simulated values can be explained by the presence of steep topography around the sill.

Variability of synoptic-scale quasi-stationary thermohaline stratification patterns in the Gulf of Finland in summer 2009 - T. Liblik; U. Lips (Ocean Science, 2012)
We present and analyze high-resolution observational data of thermohaline structure and currents acquired in the Gulf of Finland (Baltic Sea), using an autonomous buoy profiler and bottom-mounted acoustic Doppler current profiler during July–August 2009. Vertical profiles of temperature and salinity were measured in the upper 50-m layer with a 3 h time resolution, and vertical profiles of current velocity and direction were recorded with a 10 min time resolution. Although large temporal variations of vertical temperature and salinity distributions were revealed, it was possible to define several periods with quasi-stationary vertical thermohaline structure. These quasi-stationary stratification patterns persisted for 4–15 days and were dominated by certain physical processes: upwelling, relaxation of upwelling, estuarine circulation and its wind-induced reversal, and downwelling. Vertical profiles of current velocities supported the concept of synoptic-scale, quasi-stationary periods of hydrophysical fields, characterized by distinct layered flow structures and current oscillations. To estimate the contribution of different processes to the changes in stratification, a simple conceptual model was developed. The model accounts for heat flux through the sea surface, wind mixing, wind-induced transport (parallel to the horizontal salinity gradient) in the upper layer, and estuarine circulation. It reproduced observed changes in vertical stratification reasonably well. The largest discrepancies between observations and model results were found when water motions across the Gulf and associated vertical displacements of isopycnals (upwelling or downwelling) were dominant processes.

Estimation of positive sum-to-one constrained zooplankton grazing preferences with the DEnKF: a twin experiment - E. Simon; A. Samuelsen; L. Bertino; D. Dumont (Ocean Science, 2012)
We consider the estimation of the grazing preferences parameters of zooplankton in ocean ecosystem models with ensemble-based Kalman filters. These parameters are introduced to model the relative diet composition of zooplankton that consists of phytoplankton, small size-classes of zooplankton and detritus. They are positive values and their sum is equal to one. However, the sum-to-one constraint cannot be guaranteed by ensemble-based Kalman filters when parameters are bounded. Therefore, a reformulation of the parameterization is proposed. We investigate two types of variable transformations for the estimation of positive sum-to-one constrained parameters that lead to the estimation of a new set of parameters with normal or bounded distributions. These transformations are illustrated and discussed with twin experiments performed with the 1-D coupled model GOTM-NORWECOM with Gaussian anamorphosis extensions of the deterministic ensemble Kalman filter (DEnKF).

In situ determination of the remote sensing reflectance: an inter-comparison - G. Zibordi; K. Ruddick; I. Ansko; G. Moore; S. Kratzer; J. Icely; A. Reinart (Ocean Science, 2012)
Inter-comparison of data products from simultaneous measurements performed with independent systems and methods is a viable approach to assess the consistency of data and additionally to investigate uncertainties. Within such a context the inter-comparison called Assessment of In Situ Radiometric Capabilities for Coastal Water Remote Sensing Applications (ARC) was carried out at the Acqua Alta Oceanographic Tower in the northern Adriatic Sea to explore the accuracy of in situ data products from various in- and above-water optical systems and methods. Measurements were performed under almost ideal conditions, including a stable deployment platform, clear sky, relatively low sun zenith angles and moderately low sea state. Additionally, all optical sensors involved in the experiment were inter-calibrated through absolute radiometric calibration performed with the same standards and methods. Inter-compared data products include spectral water-leaving radiance Lw (λ), above-water downward irradiance Ed(0+,λ) and remote sensing reflectance Rrs(λ). Data products from the various measurement systems/methods were directly compared to those from a single reference system/method. Results for Rrs(λ) indicate spectrally averaged values of relative differences comprised between −1 and +6%, while spectrally averaged values of absolute differences vary from approximately 6% for the above-water systems/methods to 9% for buoy-based systems/methods. The agreement between Rrs(λ) spectral relative differences and estimates of combined uncertainties of the inter-compared systems/methods is noteworthy.

Validation of FOAM near-surface ocean current forecasts using Lagrangian drifting buoys - E. W. Blockley; M. J. Martin; P. Hyder (Ocean Science, 2012)
In this study, the quality of near-surface current forecasts from the FOAM ocean forecasting system is assessed using the trajectories of Lagrangian drifting buoys. A method is presented for deriving pseudo-Eulerian estimates of ocean currents from the positions of Surface Velocity Program drifters and the resulting data are compared to velocities observed by the global tropical moored buoy array. A quantitative analysis of the global FOAM velocities is performed for the period 2007 and 2008 using currents derived from over 3000 unique drifters (providing an average of 650 velocity observations per day). A potential bias is identified in the Southern Ocean which appears to be caused by wind-slip in the drifter dataset as a result of drogue loss. The drifter-derived currents are also used to show how the data assimilation scheme and a recent system upgrade impact upon the quality of FOAM current forecasts.

The 2011 marine heat wave in Cockburn Sound, southwest Australia - T. H. Rose; D. A. Smale; G. Botting (Ocean Science, 2012)
Over 2000 km of Western Australian coastline experienced a significant marine heat wave in February and March 2011. Seawater temperature anomalies of +2–4 °C were recorded at a number of locations, and satellite-derived SSTs (sea surface temperatures) were the highest on record. Here, we present seawater temperatures from southwestern Australia and describe, in detail, the marine climatology of Cockburn Sound, a large, multiple-use coastal embayment. We compared temperature and dissolved oxygen levels in 2011 with data from routine monitoring conducted from 2002–2010. A significant warming event, 2–4 °C in magnitude, persisted for > 8 weeks, and seawater temperatures at 10 to 20 m depth were significantly higher than those recorded in the previous 9 yr. Dissolved oxygen levels were depressed at most monitoring sites, being ~ 2 mg l−1 lower than usual in early March 2011. Ecological responses to short-term extreme events are poorly understood, but evidence from elsewhere along the Western Australian coastline suggests that the heat wave was associated with high rates of coral bleaching; fish, invertebrate and macroalgae mortalities; and algal blooms. However, there is a paucity of historical information on ecologically-sensitive habitats and taxa in Cockburn Sound, so that formal examinations of biological responses to the heat wave were not possible. The 2011 heat wave provided insights into conditions that may become more prevalent in Cockburn Sound, and elsewhere, if the intensity and frequency of short-term extreme events increases as predicted.

Propagation and dissipation of internal tides in the Oslofjord - A. Staalstrøm; E. Aas; B. Liljebladh (Ocean Science, 2012)
Observations of velocity, pressure, temperature and salinity in the inner Oslofjord have been analysed to provide new information about the relationships between internal tides generated by tidal currents across the Drøbak Sill and dissipation and diffusivity in the fjord. The most energetic vertical displacement of density surfaces inside the sill is associated with the first internal mode that has maximum amplitude around sill depth. The amplitude of the vertical displacement around sill depth correlates with the amplitude of the surface elevation, and, at a distance of 1 km inside the sill, the ratio between the amplitudes is 38, decreasing to 11 at a distance of 10 km. The greatest vertical displacements inside the sill, however, are found at 40 m depth. These latter internal waves are not associated with a first-mode internal tide, but are rather associated with higher internal modes controlled by stratification. The energy flux of the internal wave propagating from the Drøbak Sill into the inner fjord on the east side of the Håøya Island is estimated to vary in the range 155–430 kW. This is the same order of magnitude as the estimated barotropic energy loss over the Drøbak Sill (250 kW), but only 4–10% of the total barotropic flux. Approximately 40–70% of the internal energy flux is lost within a distance of 10 km from the sill. The mean diffusivity below 90 m depth in this area (~20 cm2 s−1) is more than four times higher than in the rest of the fjord (~5 cm2 s−1 or less).

Atlantic transport variability at 25° N in six hydrographic sections - C. P. Atkinson; H. L. Bryden; S. A. Cunningham; B. A. King (Ocean Science, 2012)
In January and February 2010, a sixth transatlantic hydrographic section was completed across 25° N, extending the hydrographic record at this latitude to over half a century. In combination with continuous transport measurements made since 2004 at 26.5° N by the Rapid-WATCH project, we reassess transport variability in the 25° N hydrographic record. Past studies of transport variability at this latitude have assumed transport estimates from each hydrographic section to represent annual average conditions. In this study the uncertainty in this assumption is assessed through use of Rapid-WATCH observations to quantify sub-seasonal and seasonal transport variability. Whilst in the upper-ocean no significant interannual or decadal transport variability are identified in the hydrographic record, in the deep ocean transport variability in both depth and potential temperature classes suggests some interannual or decadal variability may have occurred. This is particularly striking in the lower North Atlantic Deep Water where southward transports prior to 1998 were greater than recent transports by several Sverdrups. Whilst a cooling and freshening of Denmark Straits Overflow Water has occurred which is coincident with these transport changes, these water mass changes appear to be density compensated. Transport changes are the result of changing velocity shear in the vicinity of the Deep Western Boundary Current.

Impact of SLA assimilation in the Sicily Channel Regional Model: model skills and mesoscale features - A. Olita; S. Dobricic; A. Ribotti; L. Fazioli; A. Cucco; C. Dufau; R. Sorgente (Ocean Science, 2012)
The impact of the assimilation of MyOcean sea level anomalies along-track data on the analyses of the Sicily Channel Regional Model was studied. The numerical model has a resolution of 1/32° degrees and is capable to reproduce mesoscale and sub-mesoscale features. The impact of the SLA assimilation is studied by comparing a simulation (SIM, which does not assimilate data) with an analysis (AN) assimilating SLA along-track multi-mission data produced in the framework of MyOcean project. The quality of the analysis was evaluated by computing RMSE of the misfits between analysis background and observations (sea level) before assimilation. A qualitative evaluation of the ability of the analyses to reproduce mesoscale structures is accomplished by comparing model results with ocean colour and SST satellite data, able to detect such features on the ocean surface. CTD profiles allowed to evaluate the impact of the SLA assimilation along the water column. We found a significant improvement for AN solution in terms of SLA RMSE with respect to SIM (the averaged RMSE of AN SLA misfits over 2 years is about 0.5 cm smaller than SIM). Comparison with CTD data shows a questionable improvement produced by the assimilation process in terms of vertical features: AN is better in temperature while for salinity it gets worse than SIM at the surface. This suggests that a better a-priori description of the vertical error covariances would be desirable. The qualitative comparison of simulation and analyses with synoptic satellite independent data proves that SLA assimilation allows to correctly reproduce some dynamical features (above all the circulation in the Ionian portion of the domain) and mesoscale structures otherwise misplaced or neglected by SIM. Such mesoscale changes also infer that the eddy momentum fluxes (i.e. Reynolds stresses) show major changes in the Ionian area. Changes in Reynolds stresses reflect a different pumping of eastward momentum from the eddy to the mean flow, in turn influencing transports through the channel.

Operational SAR-based sea ice drift monitoring over the Baltic Sea - J. Karvonen (Ocean Science, 2012)
An algorithm for computing ice drift from pairs of synthetic aperture radar (SAR) images covering a common area has been developed at FMI. The algorithm has been developed based on the C-band SAR data over the Baltic Sea. It is based on phase correlation in two scales (coarse and fine) with some additional constraints. The algorithm has been running operationally in the Baltic Sea from the beginning of 2011, using Radarsat-1 ScanSAR wide mode and Envisat ASAR wide swath mode data. The resulting ice drift fields are publicly available as part of the MyOcean EC project. The SAR-based ice drift vectors have been compared to the drift vectors from drifter buoys in the Baltic Sea during the first operational season, and also these validation results are shown in this paper. Also some navigationally useful sea ice quantities, which can be derived from ice drift vector fields, are presented.

On the influence of wind and waves on underwater irradiance fluctuations - M. Hieronymi; A. Macke (Ocean Science, 2012)
The influence of various wind and wave conditions on the variability of downwelling irradiance Ed (490 nm) in water is subject of this study. The work is based on a two-dimensional Monte Carlo radiative transfer model with high spatial resolution. The model assumes conditions that are ideal for wave focusing, thus simulation results reveal the upper limit for light fluctuations. Local wind primarily determines the steepness of capillary-gravity waves which in turn dominate the irradiance variability near the surface. Down to 3 m depth, maximum irradiance peaks that exceed the mean irradiance Ed by a factor of more than 7 can be observed at low wind speeds up to 5 m s−1. The strength of irradiance fluctuations can be even amplified under the influence of higher ultra-gravity waves; thereby peaks can exceed 11 Ed. Sea states influence the light field much deeper; gravity waves can cause considerable irradiance variability even at 100 m depth. The simulation results show that under realistic conditions 50% radiative enhancements compared to the mean can still occur at 30 m depth. At greater depths, the underwater light variability depends on the wave steepness of the characteristic wave of a sea state; steeper waves cause stronger light fluctuations.

Transport of warm Upper Circumpolar Deep Water onto the western Antarctic Peninsula continental shelf - D. G. Martinson; D. C. McKee (Ocean Science, 2012)
Five thermistor moorings were placed on the continental shelf of the western Antarctic Peninsula (between 2007 and 2010) in an effort to identify the mechanism(s) responsible for delivering warm Upper Circumpolar Deep Water (UCDW) onto the broad continental shelf from the Antarctic Circumpolar Current (ACC) flowing over the adjacent continental slope. Historically, four mechanisms have been suggested: (1) eddies shed from the ACC, (2) flow into the cross-shelf-cutting canyons with overflow onto the nominal shelf, (3) general upwelling, and (4) episodic advective diversions of the ACC onto the shelf. The mooring array showed that for the years of deployment, the dominant mechanism is eddies; upwelling may also contribute but to an unknown extent. Mechanism 2 played no role, though the canyons have been shown previously to channel UCDW across the shelf into Marguerite Bay. Mechanism 4 played no role independently, though eddies may be advected within a greater intrusion of the background flow.

Support to oil spill emergencies in the Bonifacio Strait, western Mediterranean - A. Cucco; A. Ribotti; A. Olita; L. Fazioli; B. Sorgente; M. Sinerchia; A. Satta; A. Perilli; M. Borghini; K. Schroeder; R. Sorgente (Ocean Science, 2012)
An innovative forecasting system of the coastal marine circulation has been implemented in the Bonifacio Strait area, between Corsica and Sardinia, using a numerical approach to facilitate the rapid planning and coordination of remedial actions for oil spill emergencies at sea by local authorities. Downscaling and nesting techniques from regional to coastal scale and a 3-D hydrodynamic numerical model, coupled with a wind wave model, are the core of the integrated Bonifacio Strait system. Such a system is capable of predicting operationally the dispersion of hydrocarbon spills in the area, both in forward and backward mode, through an easy-to-use graphical user interface. A set of applications are described and discussed including both operational applications aimed at providing rapid responses to local oil spill emergences and managing applications aimed at mitigating the risk of oil spill impacts on the coast.

Influence of Ross Sea Bottom Water changes on the warming and freshening of the Antarctic Bottom Water in the Australian-Antarctic Basin - K. Shimada; S. Aoki; K. I. Ohshima; S. R. Rintoul (Ocean Science, 2012)
Changes to the properties of Antarctic Bottom Water in the Australian-Antarctic Basin (AA-AABW) between the 1990s and 2000s are documented using data from the WOCE Hydrographic Program (WHP) and repeated hydrographic surveys. Strong cooling and freshening are observed on isopycnal layers denser than γn = 28.30 kg m−3. Changes in the average salinity and potential temperature below this isopycnal correspond to a basin-wide warming of 1300 ± 200 GW and freshening of 24 ± 3 Gt year−1. Recent changes to dense shelf water in the source regions in the Ross Sea and George V Land can explain the freshening of AA-AABW but not its extensive warming. An alternative mechanism for this warming is a decrease in the supply of AABW from the Ross Sea (RSBW). Hydrographic profiles between the western Ross Sea and George V Land (171–158° E) were analyzed with a simple advective-diffusive model to assess the causes of the observed changes. The model suggests that the warming of RSBW observed between the 1970s and 2000s can be explained by a 21 ± 23% reduction in RSBW transport and the enhancement of the vertical diffusion of heat resulting from a 30 ± 7% weakening of the abyssal stratification. The documented freshening of Ross Sea dense shelf water leads to a reduction in both salinity and density stratification. Therefore the direct freshening of RSBW at its source also produces an indirect warming of the RSBW. A simple box model suggests that the changes in RSBW properties and volume transport (a decrease of 6.7% is assumed between the year 1995 and 2005) can explain 51 ± 6% of the warming and 84 ± 10% of the freshening observed in AA-AABW.

On the chemical dynamics of extracellular polysaccharides in the high Arctic surface microlayer - Q. Gao; C. Leck; C. Rauschenberg; P. A. Matrai (Ocean Science, 2012)
The surface microlayer (SML) represents a unique system of which the physicochemical characteristics may differ from those of the underlying subsurface seawater (SSW). Within the Arctic pack ice area, the SML has been characterized as enriched in small colloids of biological origin, resulting from extracellular polymeric secretions (EPS). During the Arctic Summer Cloud Ocean Study (ASCOS) in August 2008, particulate organic matter (POM, with size range > 0.22 μm) and dissolved organic matter (DOM, 5 kDa) and 64.4 ± 14.5 nM in POM. The enrichment of polysaccharides in the SML appeared to be a common feature, with EFs ranging from 1.7 to 7.0 for particulate polysaccharides and 3.5 to 12.1 for polysaccharides in the HMW DOM fraction. A calculated monosaccharide yield suggests that polymers in the HMW DOM fraction were scavenged, without substantial degradation, into the SML. Bubble scavenging experiments showed that newly aggregated particles could be formed abiotically by coagulation of low molecular weight nanometer-sized gels. Aerosol particles, artificially generated by bubbling experiments, were enriched in polysaccharides by factors of 22–70, relative to the source seawater. We propose that bubble scavenging of surface-active polysaccharides could be one of the possible mechanisms for the enrichment of polysaccharides in the high Arctic open lead SML.

Particle aggregation at the edges of anticyclonic eddies and implications for distribution of biomass - A. Samuelsen; S. S. Hjøllo; J. A. Johannessen; R. Patel (Ocean Science, 2012)
Acoustic measurements show that the biomass of zooplankton and mesopelagic fish is redistributed by mesoscale variability and that the signal extends over several hundred meters depth. The mechanisms governing this distribution are not well understood, but influences from both physical (i.e. redistribution) and biological processes (i.e. nutrient transport, primary production, active swimming, etc.) are likely. This study examines how hydrodynamic conditions and basic vertical swimming behavior act to distribute biomass in an anticyclonic eddy. Using an eddy-resolving 2.3 km-resolution physical ocean model as forcing for a particle-tracking module, particles representing passively floating organisms and organisms with vertical swimming behavior are released within an eddy and monitored for 20 to 30 days. The role of hydrodynamic conditions on the distribution of biomass is discussed in relation to the acoustic measurements. Particles released close to the surface tend, in agreement with the observations, to accumulate around the edge of the eddy, whereas particles released at depth gradually become distributed along the isopycnals. After a month they are displaced several hundreds meters in the vertical with the deepest particles found close to the eddy center and the shallowest close to the edge. There is no evidence of aggregation of particles along the eddy rim in the last simulation. The model results points towards a physical mechanism for aggregation at the surface, however biological processes cannot be ruled out using the current modeling tool.

On the outflow of dense water from the Weddell and Ross Seas in OCCAM model - R. Kerr; K. J. Heywood; M. M. Mata; C. A. E. Garcia (Ocean Science, 2012)
We describe the seasonal and interannual variability of volume transports in the Weddell and Ross Seas using the 1/12° 20-yr simulation of the OCCAM global ocean general circulation model. The average simulated full-depth cumulative volume transports were 28.5 ± 2.9 Sv (1 Sv ≡ 106 m3 s−1) and 13.4 ± 5.2 Sv, across the main export regions of the Weddell and Ross Seas, respectively. The values of mean outflow of Antarctic Bottom Water (AABW) (defined by neutral density γn ≥ 28.27 kg m−3) from the Weddell and Ross Seas of 10.6 ± 3.1 Sv and 0.5 ± 0.7 Sv, respectively, agree with the range reported in historical observational studies. The export of Weddell Sea dense water in OCCAM is primarily determined by the strength of the Weddell Gyre. Variability in AABW export is predominantly at periods of ~1 yr and 2–4 yr.

A vertical-mode decomposition to investigate low-frequency internal motion across the Atlantic at 26° N - Z. B. Szuts; J. R. Blundell; M. P. Chidichimo; J. Marotzke (Ocean Science, 2012)
Hydrographic data from full-depth moorings maintained by the Rapid/\\-MOCHA project and spanning the Atlantic at 26° N are decomposed into vertical modes in order to give a dynamical framework for interpreting the observed fluctuations. Vertical modes at each mooring are fit to pressure perturbations using a Gauss-Markov inversion. Away from boundaries, the vertical structure is almost entirely described by the first baroclinic mode, as confirmed by high correlation between the original signal and reconstructions using only the first baroclinic mode. These first baroclinic motions are also highly coherent with altimetric sea surface height (SSH). Within a Rossby radius (45 km) of the western and eastern boundaries, however, the decomposition contains significant variance at higher modes, and there is a corresponding decrease in the agreement between SSH and either the original signal or the first baroclinic mode reconstruction. Compared to the full transport signal, transport fluctuations described by the first baroclinic mode represent

Transports and budgets in a 1/4 ° global ocean reanalysis 1989–2010 - K. Haines; M. Valdivieso; H. Zuo; V. N. Stepanov (Ocean Science, 2012)
Large-scale ocean transports of heat and freshwater have not been well monitored, and yet the regional budgets of these quantities are important to understanding the role of the oceans in climate and climate change. In contrast, atmospheric heat and freshwater transports are commonly assessed from atmospheric reanalysis products, despite the presence of non-conserving data assimilation based on the wealth of distributed atmospheric observations as constraints. The ability to carry out ocean reanalyses globally at eddy-permitting resolutions of 1/4 ° or better, along with new global ocean observation programs, now makes a similar approach viable for the ocean. In this paper we examine the budgets and transports within a global high resolution ocean model constrained by ocean data assimilation, and compare them with independent oceanic and atmospheric estimates.

137 Cs off Fukushima Dai-ichi, Japan – model based estimates of dilution and fate - H. Dietze; I. Kriest (Ocean Science, 2012)
In the aftermath of an earthquake and tsunami on 11 March 2011 radioactive 137Cs was discharged from a damaged nuclear power plant to the sea off Fukushima Dai-ichi, Japan. Here we explore its dilution and fate with a state-of-the-art global ocean general circulation model, which is eddy-resolving in the region of interest. We find apparent consistency between our simulated circulation, estimates of 137Cs discharged ranging from 0.94 p Bq (Japanese Government, 2011) to 3.5 ± 0.7 p Bq (Tsumune et al., 2012), and measurements by Japanese authorities and the power plant operator. In contrast, our simulations are apparently inconsistent with the high 27 ± 15 p Bq discharge estimate of Bailly du Bois et al. (2012). Expressed in terms of a diffusivity we diagnose, from our simulations, an initial dilution on the shelf of 60 to 100 m2 s−1. The cross-shelf diffusivity is at 500 ± 300 m2 s−1 significantly higher and variable in time as indicated by its uncertainty. Expressed as an effective residence time of surface water on the shelf, the latter estimate transfers to 43 ± 16 days. As regards the fate of 137Cs, our simulations suggest that activities up to 4 mBq l−1 prevail in the Kuroshio-Oyashio Interfrontal Zone one year after the accident. This allows for low but detectable 0.1 to 0.3 m Bq l−1 entering the North Pacific Intermediate Water before the 137Cs signal is flushed away. The latter estimates concern the direct release to the sea only.

Laminar and weakly turbulent oceanic gravity currents performing inertial oscillations - A. Wirth (Ocean Science, 2012)
The small scale dynamics of a weakly turbulent oceanic gravity current is determined. The gravity current considered is initially at rest and adjusts by performing inertial oscillations to a geostrophic mean flow. The dynamics is explored with a hierarchy of mathematical models. The most involved are the fully 3-D Navier-Stokes equations subject to the Boussinesq approximation. A 1-D and 0-D mathematical model of the same gravity current dynamics are systematically derived. Using this hierarchy and the numerical solutions of the mathematical models, the turbulent dynamics at the bottom and the interface is explored and their interaction investigated. Three different regimes of the small scale dynamics of the gravity current are identified, they are characterised by laminar flow, coherent roll vortices and turbulent dynamics with coherent streaks and bursts. The problem of the rectification of the turbulent fluxes, that is, how to average out the fluctuations and calculate their average influence on the flow, is considered. It is shown that two different regimes of friction are superposed, an Ekman friction applies to the average geostrophic flow and a linear friction, not influenced by rotation, to the inertial oscillations. The combination of the two makes the bulk friction non-local in time for the 0-D model. The implications of the results for parametrisations of the Ekman dynamics and the small scale turbulent fluxes in the planetary boundary layer are discussed.

Wave climate in the Arkona Basin, the Baltic Sea - T. Soomere; R. Weisse; A. Behrens (Ocean Science, 2012)
The basic features of the wave climate in the Southwestern Baltic Sea (such as the average and typical wave conditions, frequency of occurrence of different wave parameters, variations in wave heights from weekly to decadal scales) are established based on waverider measurements at the Darss Sill in 1991–2010. The measured climate is compared with two numerical simulations with the WAM wave model driven by downscaled reanalysis of wind fields for 1958–2002 and by adjusted geostrophic winds for 1970–2007. The wave climate in this region is typical for semi-enclosed basins of the Baltic Sea. The maximum wave heights are about half of those in the Baltic Proper. The maximum recorded significant wave height HS =4.46 m occurred on 3 November 1995. The wave height exhibits no long-term trend but reveals modest interannual (about 12 % of the long-term mean of 0.76 m) and substantial seasonal variation. The wave periods are mostly concentrated in a narrow range of 2.6–4 s. Their distribution is almost constant over decades. The role of remote swell is very small.

Arctic Ocean circulation and variability – advection and external forcing encounter constraints and local processes - B. Rudels (Ocean Science, 2012)
The first hydrographic data from the Arctic Ocean, the section from the Laptev Sea to the passage between Greenland and Svalbard obtained by Nansen on his drift with Fram 1893–1896, aptly illustrate the main features of Arctic Ocean oceanography and indicate possible processes active in transforming the water masses in the Arctic Ocean. Many, perhaps most, processes were identified already by Nansen, who put his mark on almost all subsequent research in the Arctic. Here we shall revisit some key questions and follow how our understanding has evolved from the early 20th century to present. What questions, if any, can now be regarded as solved and which remain still open? Five different but connected topics will be discussed: (1) The low salinity surface layer and the storage and export of freshwater. (2) The vertical heat transfer from the Atlantic water to sea ice and to the atmosphere. (3) The circulation and mixing of the two Atlantic inflow branches. (4) The formation and circulation of deep and bottom waters in the Arctic Ocean. (5) The exchanges through Fram Strait. Foci will be on the potential effects of increased freshwater input and reduced sea ice export on the freshwater storage and residence time in the Arctic Ocean, on the deep waters of the Makarov Basin, and on the circulation and relative importance of the two inflows, over the Barents Sea and through Fram Strait, for the distribution of heat in the intermediate layers of the Arctic Ocean.

Variational assimilation of Lagrangian trajectories in the Mediterranean ocean Forecasting System - J. A. U. Nilsson; S. Dobricic; N. Pinardi; P.-M. Poulain; D. Pettenuzzo (Ocean Science, 2012)
A novel method for three-dimensional variational assimilation of Lagrangian data with a primitive-equation ocean model is proposed. The assimilation scheme was implemented in the Mediterranean ocean Forecasting System and evaluated for a 4-month period. Four experiments were designed to assess the impact of trajectory assimilation on the model output, i.e. the sea-surface height, velocity, temperature and salinity fields. It was found from the drifter and Argo trajectory assimilation experiment that the forecast skill of surface-drifter trajectories improved by 15 %, that of intermediate-depth float trajectories by 20 %, and moreover, that the forecasted sea-surface height fields improved locally by 5 % compared to satellite data, while the quality of the temperature and salinity fields remained at previous levels. In conclusion, the addition of Lagrangian trajectory assimilation proved to reduce the uncertainties in the model fields, thus yielding a higher accuracy of the ocean forecasts.

Mesoscale variability of water masses in the Arabian Sea as revealed by ARGO floats - X. Carton; P. L'Hegaret; R. Baraille (Ocean Science, 2012)
By analysing ARGO float data over the last four years, a few aspects of the mesoscale variability of water masses in the Arabian Sea are described. The Red Sea Outflow Water (RSOW) is concentrated in the Southwestern Gulf of Aden, in particular when a cyclonic gyre predominates in this region. Salinities of 36.5 and temperatures of 16 °C are found in this area at depths between 600 and 1000 m. RSOW is more dilute in the eastern part of the Gulf, where intense and relatively barotropic gyres mix it with Indian ocean Central Water. RSOW is also detected along the northeastern coast of Socotra, and fragments of RSOW are found between one and three degrees of latitude north of this island. In the whole Gulf of Aden, the correlation between the deep motions of the floats and the sea-level anomaly measured by altimetry is strong, at regional scale. The finer scale details of the float trajectories are not sampled by altimetry and are often related to the anomalous water masses that the floats encounter. The Persian Gulf Water (PGW) is found in the float profiles near Ras ash Sharbatat (near 57° E, 18° N), again with 36.5 in salinity and about 18–19 °C in temperature. These observations were achieved in winter when the southwestward monsoon currents can advect PGW along the South Arabian coast. Fragments of PGW were also observed in the Arabian Sea between 18 and 20° N and 63 and 65° E in summer, showing that this water mass can escape the Gulf of Oman southeastward, during that season. Kinetic energy distributions of floats with respect to distance or angle share common features between the two regions (Gulf of Aden and Arabian Sea), in particular peaks at 30, 50 and 150 km scales and along the axis of monsoon currents. Hydrological measurements by floats are also influenced by the seasonal variations of PGW and RSOW in these regions.

ENSURF: multi-model sea level forecast – implementation and validation results for the IBIROOS and Western Mediterranean regions - B. Pérez; R. Brouwer; J. Beckers; D. Paradis; C. Balseiro; K. Lyons; M. Cure; M. G. Sotillo; B. Hackett; M. Verlaan; E. A. Fanjul (Ocean Science, 2012)
ENSURF (Ensemble SURge Forecast) is a multi-model application for sea level forecast that makes use of several storm surge or circulation models and near-real time tide gauge data in the region, with the following main goals: 1. providing easy access to existing forecasts, as well as to its performance and model validation, by means of an adequate visualization tool; 2. generation of better forecasts of sea level, including confidence intervals, by means of the Bayesian Model Average technique (BMA). The Bayesian Model Average technique generates an overall forecast probability density function (PDF) by making a weighted average of the individual forecasts PDF's; the weights represent the Bayesian likelihood that a model will give the correct forecast and are continuously updated based on the performance of the models during a recent training period. This implies the technique needs the availability of sea level data from tide gauges in near-real time. The system was implemented for the European Atlantic facade (IBIROOS region) and Western Mediterranean coast based on the MATROOS visualization tool developed by Deltares. Results of validation of the different models and BMA implementation for the main harbours are presented for these regions where this kind of activity is performed for the first time. The system is currently operational at Puertos del Estado and has proved to be useful in the detection of calibration problems in some of the circulation models, in the identification of the systematic differences between baroclinic and barotropic models for sea level forecasts and to demonstrate the feasibility of providing an overall probabilistic forecast, based on the BMA method.

An empirical model for the statistics of sea surface diurnal warming - M. J. Filipiak; C. J. Merchant; H. Kettle; P. Le Borgne (Ocean Science, 2012)
A statistical model is derived relating the diurnal variation of sea surface temperature (SST) to the net surface heat flux and surface wind speed from a numerical weather prediction (NWP) model. The model is derived using fluxes and winds from the European Centre for Medium-Range Weather Forecasting (ECMWF) NWP model and SSTs from the Spinning Enhanced Visible and Infrared Imager (SEVIRI). In the model, diurnal warming has a linear dependence on the net surface heat flux integrated since (approximately) dawn and an inverse quadratic dependence on the maximum of the surface wind speed in the same period. The model coefficients are found by matching, for a given integrated heat flux, the frequency distributions of the maximum wind speed and the observed warming. Diurnal cooling, where it occurs, is modelled as proportional to the integrated heat flux divided by the heat capacity of the seasonal mixed layer. The model reproduces the statistics (mean, standard deviation, and 95-percentile) of the diurnal variation of SST seen by SEVIRI and reproduces the geographical pattern of mean warming seen by the Advanced Microwave Scanning Radiometer (AMSR-E). We use the functional dependencies in the statistical model to test the behaviour of two physical model of diurnal warming that display contrasting systematic errors.

Black Sea coastal forecasting system - A. I. Kubryakov; G. K. Korotaev; V. L. Dorofeev; Y. B. Ratner; A. Palazov; N. Valchev; V. Malciu; R. Mateescu; T. Oguz (Ocean Science, 2012)
The Black Sea coastal nowcasting and forecasting system was built within the framework of EU FP6 ECOOP (European COastalshelf sea OPerational observing and forecasting system) project for five regions: the south-western basin along the coasts of Bulgaria and Turkey, the north-western shelf along the Romanian and Ukrainian coasts, coastal zone around of the Crimea peninsula, the north-eastern Russian coastal zone and the coastal zone of Georgia. The system operates in the real-time mode during the ECOOP project and afterwards. The forecasts include temperature, salinity and current velocity fields. Ecosystem model operates in the off-line mode near the Crimea coast.

Technical Note: Detection of gas bubble leakage via correlation of water column multibeam images - J. Schneider von Deimling; C. Papenberg (Ocean Science, 2012)
Hydroacoustic detection of natural gas release from the seafloor has been conducted in the past by using singlebeam echosounders. In contrast, modern multibeam swath mapping systems allow much wider coverage, higher resolution, and offer 3-D spatial correlation. Up to the present, the extremely high data rate hampers water column backscatter investigations and more sophisticated visualization and processing techniques are needed. Here, we present water column backscatter data acquired with a 50 kHz prototype multibeam system over a period of 75 seconds. Display types are of swath-images as well as of a "re-sorted" singlebeam presentation. Thus, individual and/or groups of gas bubbles rising from the 24 m deep seafloor clearly emerge in the acoustic images, making it possible to estimate rise velocities. A sophisticated processing scheme is introduced to identify those rising gas bubbles in the hydroacoustic data. We apply a cross-correlation technique adapted from particle imaging velocimetry (PIV) to the acoustic backscatter images. Temporal and spatial drift patterns of the bubbles are assessed and are shown to match very well to measured and theoretical rise patterns. The application of this processing to our field data gives clear results with respect to unambiguous bubble detection and remote bubble rise velocimetry. The method can identify and exclude the main source of misinterpretations, i.e. fish-mediated echoes. Although image-based cross-correlation techniques are well known in the field of fluid mechanics for high resolution and non-inversive current flow field analysis, we present the first application of this technique as an acoustic bubble detector.

Preface "An historical perspective on the development of the Thermodynamic Equation of Seawater – 2010" - R. Pawlowicz; T. McDougall; R. Feistel; R. Tailleux (Ocean Science, 2012)

Towards a regional ocean forecasting system for the IBI (Iberia-Biscay-Ireland area): developments and improvements within the ECOOP project framework - S. Cailleau; J. Chanut; J.-M. Lellouche; B. Levier; C. Maraldi; G. Reffray; M. G. Sotillo (Ocean Science, 2012)
The regional ocean operational system remains a key element in downscaling from large scale (global or basin scale) systems to coastal ones. It enables the transition between systems in which the resolution and the resolved physics are quite different. Indeed, coastal applications need a system to predict local high frequency events (inferior to the day) such as storm surges, while deep sea applications need a system to predict large scale lower frequency ocean features. In the framework of the ECOOP project, a regional system for the Iberia-Biscay-Ireland area has been upgraded from an existing V0 version to a V2. This paper focuses on the improvements from the V1 system, for which the physics are close to a large scale basin system, to the V2 for which the physics are more adapted to shelf and coastal issues. Strong developments such as higher regional physics resolution in the NEMO Ocean General Circulation Model for tides, non linear free surface and adapted vertical mixing schemes among others have been implemented in the V2 version. Thus, regional thermal fronts due to tidal mixing now appear in the latest version solution and are quite well positioned. Moreover, simulation of the stratification in shelf areas is also improved in the V2.

Towards an improved description of ocean uncertainties: effect of local anamorphic transformations on spatial correlations - J.-M. Brankart; C.-E. Testut; D. Béal; M. Doron; C. Fontana; M. Meinvielle; P. Brasseur; J. Verron (Ocean Science, 2012)
The objective of this paper is to investigate if the description of ocean uncertainties can be significantly improved by applying a local anamorphic transformation to each model variable, and by making the assumption of joint Gaussianity for the transformed variables, rather than for the original variables. For that purpose, it is first argued that a significant improvement can already be obtained by deriving the local transformations from a simple histogram description of the marginal distributions. Two distinctive advantages of this solution for large size applications are the conciseness and the numerical efficiency of the description. Second, various oceanographic examples are used to evaluate the effect of the resulting piecewise linear local anamorphic transformations on the spatial correlation structure. These examples include (i) stochastic ensemble descriptions of the effect of atmospheric uncertainties on the ocean mixed layer, and of wind uncertainties or parameter uncertainties on the ecosystem, and (ii) non-stochastic ensemble descriptions of forecast uncertainties in current sea ice and ecosystem pre-operational developments. The results indicate that (i) the transformation is accurate enough to faithfully preserve the correlation structure if the joint distribution is already close to Gaussian, and (ii) the transformation has the general tendency of increasing the correlation radius as soon as the spatial dependence between random variables becomes nonlinear, with the important consequence of reducing the number of degrees of freedom in the uncertainties, and thus increasing the benefit that can be expected from a given observation network.

Modeling of wave-induced irradiance variability in the upper ocean mixed layer - M. Hieronymi; A. Macke; O. Zielinski (Ocean Science, 2012)
A Monte Carlo based radiative transfer model has been developed for calculating the availability of solar radiation within the top 100 m of the ocean. The model is optimized for simulations of spatial high resolution downwelling irradiance Ed fluctuations that arise from the lensing effect of waves at the water surface. In a first step the accuracy of simulation results has been verified by measurements of the oceanic underwater light field and through intercomparison with an established radiative transfer model. Secondly the potential depth-impact of nonlinear shaped single waves, from capillary to swell waves, is assessed by considering the most favorable conditions for light focusing, i.e. monochromatic light at 490 nm, very clear oceanic water with a low chlorophyll a content of 0.1 mg m−3 and high sun elevation. Finally light fields below irregular wave profiles accounting for realistic sea states were simulated. Our simulation results suggest that under open ocean conditions light flashes with 50% irradiance enhancements can appear down to 35 m depth, and light variability in the range of ±10% compared to the mean Ed is still possible in 100 m depth.

A multi-decadal meridional displacement of the Subpolar Front in the Newfoundland Basin - I. Núñez-Riboni; M. Bersch; H. Haak; J. H. Jungclaus; K. Lohmann (Ocean Science, 2012)
Observations since the 1950s show a multi-decadal cycle of a meridional displacement of the Subpolar Front (SPF) in the Newfoundland Basin (NFB) in the North Atlantic. The SPF displacement is associated with corresponding variations in the path of the North Atlantic Current. We use the ocean general circulation model MPIOM with enhanced horizontal and vertical resolutions and forced with NCEP/NCAR reanalysis data to study the relation of the SPF displacement to atmospheric forcing, intensities of the subpolar gyre (SPG) and Meridional Overturning Circulation (MOC), and Labrador Sea Water (LSW) volume. The simulations indicate that the SPF displacement is associated with a circulation anomaly between the SPG and the subtropical gyre (STG), an inter-gyre gyre with a multi-decadal time scale. A sensitivity experiment indicates that both wind stress curl (WSC) and heat fluxes (which match LSW changes) contribute to the circulation anomalies in the frontal region and to the SPF displacement. An anticyclonic inter-gyre gyre is related to negative WSC and LSW anomalies and to a SPF north of its climatological position, indicating an expanding STG. A cyclonic inter-gyre gyre is related to positive WSC and LSW anomalies and a SPF south of its climatological position, indicating an expanding SPG. Therefore, the mean latitudinal position of the SPF in the NFB (a "SPF index") could be an indicator of the amount of LSW in the inter-gyre region. Spreading of LSW anomalies intensifies the MOC, suggesting our SPF index as predictor of the MOC intensity at multi-decadal time scales. The meridional displacement of the SPF has a pronounced influence on the meridional heat transport, both on its gyre and overturning components.

Sea surface temperature anomalies, seasonal cycle and trend regimes in the Eastern Pacific coast - A. Ramos-Rodríguez; D. B. Lluch-Cota; S. E. Lluch-Cota; A. Trasviña-Castro (Ocean Science, 2012)
We used the extended reconstruction of sea surface temperature (ERSST) to analyze the variation of surface temperature and the seasonal cycle along the coast of the eastern Pacific (60° N–60° S, 61 pixels alongshore) from 1950 to 2010 (732 months). First, we analyzed the monthly anomalies and looked for a relationship of such anomalies with total solar irradiance (TSI) and then the Regime Shift Detector (RSD) was applied to detect possible temperature regimes in the series. Afterwards, we calculated a yearly temperature range per pixel (amplitude of seasonal cycle) and through the subtraction of a latitudinal theoretical curve of temperature based on solar irradiance, the residuals of the seasonal cycle were obtained. The results showed an almost complete spatial synchrony and dominance of negative anomalies from 1950 to mid-late 1970's, with a switch to near-zero and positive anomalies that lasted up to late 1990's when a new shift to negative values was detected. Such a shift lasted until the early 2000's when positive anomalies appeared again but there was a change to negative anomalies in the late 2000's. These results were supported by the RSD. The TSI variability shows a clear relationship with that of sea surface temperature anomalies and with the regime changes. This is probably due to a difference in the amount of energy received from the sun. Comparing the "cool regime" versus the "warm regime", the second one received 0.39% more energy (approximately 3 × 108 J m−2) from the sun. Seasonal cycles show larger ranges at northern latitudes (>40° N), northern tropical-temperate transition zone (20°–26° N) and in the tropical-equatorial band (0°–30° S). The smallest ranges occur at 0°–16° N and 50°–60° S. The residuals (seasonal minus the theoretical curve) indicated a clear modulation due to advection by ocean currents.

Impact of combining GRACE and GOCE gravity data on ocean circulation estimates - T. Janjić; J. Schröter; R. Savcenko; W. Bosch; A. Albertella; R. Rummel; O. Klatt (Ocean Science, 2012)
With the focus on the Southern Ocean circulation, results of assimilation of multi-mission-altimeter data and the GRACE/GOCE gravity data into the finite element ocean model (FEOM) are investigated. We use the geodetic method to obtain the dynamical ocean topography (DOT). This method combines the multi-mission-altimeter sea surface height and the GRACE/GOCE gravity field. Using the profile approach, the spectral consistency of both fields is achieved by filtering the sea surface height and the geoid. By combining the GRACE and GOCE data, a considerably shorter filter length can be used, which results in more DOT details. We show that this increase in resolution of measured DOT carries onto the results of data assimilation for the surface data. By assimilating only absolute dynamical topography data using the ensemble Kalman filter, we were able to improve modeled fields. Results are closer to observations which were not used for assimilation and lie outside the area covered by altimetry in the Southern Ocean (e.g. temperature of surface drifters or deep temperatures in the Weddell Sea area at 800 m depth derived from Argo composite.)

A simple and self-consistent geostrophic-force-balance model of the thermohaline circulation with boundary mixing - J. Callies; J. Marotzke (Ocean Science, 2012)
A simple model of the thermohaline circulation (THC) is formulated, with the objective to represent explicitly the geostrophic force balance of the basinwide THC. The model comprises advective-diffusive density balances in two meridional-vertical planes located at the eastern and the western walls of a hemispheric sector basin. Boundary mixing constrains vertical motion to lateral boundary layers along these walls. Interior, along-boundary, and zonally integrated meridional flows are in thermal-wind balance. Rossby waves and the absence of interior mixing render isopycnals zonally flat except near the western boundary, constraining meridional flow to the western boundary layer. The model is forced by a prescribed meridional surface density profile. This two-plane model reproduces both steady-state density and steady-state THC structures of a primitive-equation model. The solution shows narrow deep sinking at the eastern high latitudes, distributed upwelling at both boundaries, and a western boundary current with poleward surface and equatorward deep flow. The overturning strength has a 2/3-power-law dependence on vertical diffusivity and a 1/3-power-law dependence on the imposed meridional surface density difference. Convective mixing plays an essential role in the two-plane model, ensuring that deep sinking is located at high latitudes. This role of convective mixing is consistent with that in three-dimensional models and marks a sharp contrast with previous two-dimensional models. Overall, the two-plane model reproduces crucial features of the THC as simulated in simple-geometry three-dimensional models. At the same time, the model self-consistently makes quantitative a conceptual picture of the three-dimensional THC that hitherto has been expressed either purely qualitatively or not self-consistently.

The impacts of physical processes on oxygen variations in the North Sea-Baltic Sea transition zone - L. Jonasson; J. L. S. Hansen; Z. Wan; J. She (Ocean Science, 2012)
The bottom water of the North Sea–Baltic Sea transition zone suffers from seasonal hypoxia, usually during late summer and autumn. These hypoxic events are critical for the benthic ecosystems and the concentration of dissolved oxygen is an important measure of the water quality. However, to model the subsurface dissolved oxygen is a major challenge, especially in estuaries and coastal regions. In this study a simple oxygen consumption model is coupled to a 3-D hydrodynamical model in order to analyse oxygen variations in the transition zone. The benthic and pelagic consumption of oxygen is modelled as a function of water temperature and oxygen concentration. A quantitative assessment of the model demonstrates that the model is able to resolve both seasonal and interannual variations in dissolved oxygen. Results from several experimental simulations highlight the importance of physical processes in the regulation of dissolved oxygen. Advective oxygen transport and wind induced mixing are two key processes that control the extent of hypoxia in the transition zone.

The vertical structure of oceanic Rossby waves: a comparison of high-resolution model data to theoretical vertical structures - F. K. Hunt; R. Tailleux; J. J.-M. Hirschi (Ocean Science, 2012)
Tests of the new Rossby wave theories that have been developed over the past decade to account for discrepancies between theoretical wave speeds and those observed by satellite altimeters have focused primarily on the surface signature of such waves. It appears, however, that the surface signature of the waves acts only as a rather weak constraint, and that information on the vertical structure of the waves is required to better discriminate between competing theories. Due to the lack of 3-D observations, this paper uses high-resolution model data to construct realistic vertical structures of Rossby waves and compares these to structures predicted by theory. The meridional velocity of a section at 24° S in the Atlantic Ocean is pre-processed using the Radon transform to select the dominant westward signal. Normalized profiles are then constructed using three complementary methods based respectively on: (1) averaging vertical profiles of velocity, (2) diagnosing the amplitude of the Radon transform of the westward propagating signal at different depths, and (3) EOF analysis. These profiles are compared to profiles calculated using four different Rossby wave theories: standard linear theory (SLT), SLT plus mean flow, SLT plus topographic effects, and theory including mean flow and topographic effects. Our results support the classical theoretical assumption that westward propagating signals have a well-defined vertical modal structure associated with a phase speed independent of depth, in contrast with the conclusions of a recent study using the same model but for different locations in the North Atlantic. The model structures are in general surface intensified, with a sign reversal at depth in some regions, notably occurring at shallower depths in the East Atlantic. SLT provides a good fit to the model structures in the top 300 m, but grossly overestimates the sign reversal at depth. The addition of mean flow slightly improves the latter issue, but is too surface intensified. SLT plus topography rectifies the overestimation of the sign reversal, but overestimates the amplitude of the structure for much of the layer above the sign reversal. Combining the effects of mean flow and topography provided the best fit for the mean model profiles, although small errors at the surface and mid-depths are carried over from the individual effects of mean flow and topography respectively. Across the section the best fitting theory varies between SLT plus topography and topography with mean flow, with, in general, SLT plus topography performing better in the east where the sign reversal is less pronounced. None of the theories could accurately reproduce the deeper sign reversals in the west. All theories performed badly at the boundaries. The generalization of this method to other latitudes, oceans, models and baroclinic modes would provide greater insight into the variability in the ocean, while better observational data would allow verification of the model findings.

Mapping turbidity layers using seismic oceanography methods - E. A. Vsemirnova; R. W. Hobbs; P. Hosegood (Ocean Science, 2012)
Using a combination of seismic oceanographic and physical oceanographic data acquired across the Faroe-Shetland Channel we present evidence of a turbidity layer that transports suspended sediment along the western boundary of the Channel. We focus on reflections observed on seismic data close to the sea-bed on the Faroese side of the Channel below 900 m. Forward modelling based on independent physical oceanographic data show that thermohaline structure does not explain these near sea-bed reflections but they are consistent with optical backscatter data, dry matter concentrations from water samples and from seabed sediment traps. Hence we conclude that an impedance contrast in water column caused by turbidity layers is strong enough to be seen in seismic sections and this provides a new way to visualise this type of current and its lateral structure. By inverting the seismic data we estimate a sediment concentration in the turbidity layers, present at the time of the survey, of 45 ± 25 mg l−1. We believe this is the first direct observation of a turbidity current using Seismic Oceanography.

Using dissolved oxygen concentrations to determine mixed layer depths in the Bellingshausen Sea - K. Castro-Morales; J. Kaiser (Ocean Science, 2012)
Concentrations of oxygen (O2) and other dissolved gases in the oceanic mixed layer are often used to calculate air-sea gas exchange fluxes. The mixed layer depth (zmix) may be defined using criteria based on temperature or density differences to a reference depth near the ocean surface. However, temperature criteria fail in regions with strong haloclines such as the Southern Ocean where heat, freshwater and momentum fluxes interact to establish mixed layers. Moreover, the time scales of air-sea exchange differ for gases and heat, so that zmix defined using oxygen may be different than zmix defined using temperature or density. Here, we propose to define an O2-based mixed layer depth, zmix(O2), as the depth where the relative difference between the O2 concentration and a reference value at a depth equivalent to 10 dbar equals 0.5 %. This definition was established by analysis of O2 profiles from the Bellingshausen Sea (west of the Antarctic Peninsula) and corroborated by visual inspection. Comparisons of zmix(O2) with zmix based on potential temperature differences, i.e., zmix(0.2 °C) and zmix(0.5 °C), and potential density differences, i.e., zmix(0.03 kg m−3) and zmix(0.125 kg m−3), showed that zmix(O2) closely follows zmix(0.03 kg m−3). Further comparisons with published zmix climatologies and zmix derived from World Ocean Atlas 2005 data were also performed. To establish zmix for use with biological production estimates in the absence of O2 profiles, we suggest using zmix(0.03 kg m−3), which is also the basis for the climatology by de Boyer Montégut et al. (2004).

Extraction of spatial-temporal rules from mesoscale eddies in the South China Sea based on rough set theory - Y. Du; X. Fan; Z. He; F. Su; C. Zhou; H. Mao; D. Wang (Ocean Science, 2011)
In this paper, a rough set theory is introduced to represent spatial-temporal relationships and extract the corresponding rules from typical mesoscale-eddy states in the South China Sea (SCS). Three decision attributes are adopted in this study, which make the approach flexible in retrieving spatial-temporal rules with different features. The results demonstrate that this approach is effective, and therefore provides a powerful approach to forecasts in the future studies. Spatial-temporal rules in the SCS indicate that warm eddies following the rules are generally in the southeastern and central SCS around 2000 m isobaths in winter. Their intensity and vorticity are weaker than those of cold eddies. They usually move a shorter distance. By contrast, cold eddies are in 2000 m and deeper regions of the southwestern and northeastern SCS in spring and fall. Their intensity and vorticity are strong. Usually they move a long distance. In winter, a few rules are followed by cold eddies in the northern tip of the basin and southwest of Taiwan Island rather than warm eddies, indicating cold eddies may be well-regulated in the region. Several warm-eddy rules are achieved west of Luzon Island, indicating warm eddies may be well-regulated in the region as well. Otherwise, warm and cold eddies are distributed not only in the jet flow off southern Vietnam induced by intraseasonal wind stress in summer-fall, but also in the northern shallow water, which should be a focus of a future study.

Wind forcing of salinity anomalies in the Denmark Strait overflow - S. Hall; S. R. Dye; K. J. Heywood; M. R. Wadley (Ocean Science, 2011)
The overflow of dense water from the Nordic Seas to the North Atlantic through Denmark Strait is an important part of the global thermohaline circulation. The salinity of the overflow plume has been measured by an array of current meters across the continental slope off the coast of Angmagssalik, southeast Greenland since September 1998. During 2004 the salinity of the overflow plume changed dramatically; the entire width of the array (70 km) freshened between January 2004 and July 2004, with a significant negative salinity anomaly of about 0.06 in May. The event in May represents a fresh anomaly of over 3 standard deviations from the mean since recording began in 1998. The OCCAM 1/12° Ocean General Circulation Model not only reproduces the 2004 freshening event (r=0.96, p

Usefulness of high resolution coastal models for operational oil spill forecast: the "Full City" accident - G. Broström; A. Carrasco; L. R. Hole; S. Dick; F. Janssen; J. Mattsson; S. Berger (Ocean Science, 2011)
Oil spill modeling is considered to be an important part of a decision support system (DeSS) for oil spill combatment and is useful for remedial action in case of accidents, as well as for designing the environmental monitoring system that is frequently set up after major accidents. Many accidents take place in coastal areas, implying that low resolution basin scale ocean models are of limited use for predicting the trajectories of an oil spill. In this study, we target the oil spill in connection with the "Full City" accident on the Norwegian south coast and compare operational simulations from three different oil spill models for the area. The result of the analysis is that all models do a satisfactory job. The "standard" operational model for the area is shown to have severe flaws, but by applying ocean forcing data of higher resolution (1.5 km resolution), the model system shows results that compare well with observations. The study also shows that an ensemble of results from the three different models is useful when predicting/analyzing oil spill in coastal areas.

Operational forecast of hydrophysical fields in the Georgian Black Sea coastal zone within the ECOOP - A. A. Kordzadze; D. I. Demetrashvili (Ocean Science, 2011)
One of the parts of the Black Sea Nowcasting/Forecasting System is the regional forecasting system for the easternmost part of the Black Sea (including the Georgian water area), which has been developed within the context of the EU International projects ARENA and ECOOP. A core of the regional system is a high-resolution baroclinic regional model of the Black Sea dynamics developed at M. Nodia Institute of Geophysics (RM-IG). This model is nested in the basin-scale model of Marine Hydrophysical Institute (MHI, Sevastopol/Ukraine). The regional area is limited to the Caucasian and Turkish coastal lines and the western liquid boundary coinciding with the meridian 39.36° E. Since June 2010 we have regularly been computing 3 days' forecasts of current, temperature and salinity for the easternmost part of the Black Sea with 1 km spacing. In this study the results of two forecasts are presented. The first forecast corresponds to summer season and covers the prognostic interval from 00:00 h, 6 August to 00:00 h, 9 August 2010. The second one corresponds to autumn season and covers the prognostic interval from 00:00 h, 26 October to 00:00 h, 29 October 2010. Data needed for the forecasts – the initial and prognostic hydrophysical fields on the open boundary, also 2-D prognostic meteorological fields at the sea surface – wind stress, heat fluxes, evaporation and precipitation rates for our regional area are being placed on the MHI server every day and we are available to use these data operatively. Prognostic hydrophysical fields are results of forecast by the basin-scale model of MHI and 2-D meteorological boundary fields represent the results of forecast by regional atmospheric model ALADIN. All these fields are given on the grid of basin-scale model with 5 km spacing and with one-hour time step frequency for the integration period. The analysis of predicted fields shows that to use the model with high resolution is very important factor for identification of nearshore eddies of small sizes. It should be noted the very different character of regional circulation in summer and autumn seasons in the easternmost part of the Black Sea.

How well can we derive Global Ocean Indicators from Argo data? - K. von Schuckmann; P.-Y. Le Traon (Ocean Science, 2011)
Argo deployments began in the year 2000 and by November 2007, the array reached its initial goal of 3000 floats operating worldwide. In this study, Argo temperature and salinity measurements during the period 2005 to 2010 are used to estimate Global Ocean Indicators (GOIs) such as global ocean heat content (GOHC), global ocean freshwater content (GOFC) and global steric sea level (GSSL). We developed a method based on a simple box averaging scheme using a weighted mean. Uncertainties due to data processing methods and choice of climatology are estimated. This method is easy to implement and run and can be used to set up a routine monitoring of the global ocean. Over the six year time period, trends of GOHC and GSSL are 0.54 ± 0.1 W m−2 and 0.75 ± 0.15 mm yr−1, respectively. The trend of GOFC is barely significant. Results show that there is significant interannual variability at global scale, especially for GOFC. Annual mean GOIs from the today's Argo sampling can be derived with an accuracy of ±0.11 cm for GSSL, ±0.22 × 108 J m−2 for GOHC, and ±700 km3 for GOFC. Long-term trends (15 yr) of GOIs based on the complete Argo sampling for the upper 1500 m depth can be estimated with an accuracy of ±0.04 mm yr−1 for GSSL, ±0.02 W m−2 for GOHC and ±20 km3 yr−1 for GOFC – under the assumption that no systematic errors remain in the observing system.

A pre-operational three Dimensional variational data assimilation system in the North/Baltic Sea - S. Y. Zhuang; W. W. Fu; J. She (Ocean Science, 2011)
This paper describes the implementation and evaluation of a pre-operational three dimensional variational (3DVAR) data assimilation system for the North/Baltic Sea. Univariate analysis for both temperature and salinity is applied in a 3DVAR scheme in which the horizontal component of the background error covariance is modeled by an isotropic recursive filter (IRF) and the vertical component is represented by dominant Empirical Orthogonal Functions (EOFs). Observations of temperature and salinity (T/S) profiles in the North/Baltic Sea are assimilated in the year of 2005. Effect of the 3DVAR scheme is assessed by a comparison between data assimilation run and control run. The statistical analysis indicates that the model simulation is significantly improved with the 3DVAR scheme. On average, the root mean square errors (RMSE) of temperature and salinity are reduced by 0.2 °C and 0.25 psu in the North/Baltic Sea. In addition, the bias of temperature and salinity is also decreased by 0.1 °C and 0.2 psu, respectively. Starting from an analyzed initial state, one month simulation without assimilation is carried out with the aim of examining the persistence of the initial impact. It is shown that the assimilated initial state can impact the model simulation for nearly two weeks. The influence on salinity is more pronounced than temperature.

An empirical stochastic model of sea-surface temperatures and surface winds over the Southern Ocean - S. Kravtsov; D. Kondrashov; I. Kamenkovich; M. Ghil (Ocean Science, 2011)
This study employs NASA's recent satellite measurements of sea-surface temperatures (SSTs) and sea-level winds (SLWs) with missing data filled-in by Singular Spectrum Analysis (SSA), to construct empirical models that capture both intrinsic and SST-dependent aspects of SLW variability. The model construction methodology uses a number of algorithmic innovations that are essential in providing stable estimates of the model's propagator. The best model tested herein is able to faithfully represent the time scales and spatial patterns of anomalies associated with a number of distinct processes. These processes range from the daily synoptic variability to interannual signals presumably associated with oceanic or coupled dynamics. Comparing the simulations of an SLW model forced by the observed SST anomalies with the simulations of an SLW-only model provides preliminary evidence for the ocean driving the atmosphere in the Southern Ocean region.

Mean Dynamic Topography of the Black Sea, computed from altimetry, drifter measurements and hydrology data - A. A. Kubryakov; S. V. Stanichny (Ocean Science, 2011)
Mean Dynamic Topography (MDT) is a crucial parameter for estimating dynamic topography, and, therefore, geostrophic circulation from satellite altimetry measurements. In this work we use drifting buoy measurements, hydrographic profiles and along-track Sea Level Anomalies (SLA) to reconstruct MDT of the Black Sea by the "synthetic" method. Obtained MDT shows a lot of mesoscale features, which are not present in previous MDT fields of the Black Sea, mostly based on climatic data. Moreover, gradients of sea level in the synthetic MDT are significantly higher compared to other fields, which is evidence of more intense currents in the basin. Validation of determined MDT field with independent dynamic heights and drifter buoy velocities shows good quantitative and qualitative coincidence over all Black Sea basin and improvements compare to previous fields. New Black Sea MDT will improve quality of altimetry-derived geostrophic velocities and lead to better understanding of the spatial and temporal features of the upper layer dynamics.

A computational method for determining XBT depths - J. Stark; J. Gorman; M. Hennessey; F. Reseghetti; J. Willis; J. Lyman; J. Abraham; M. Borghini (Ocean Science, 2011)
A new technique for determining the depth of expendable bathythermographs (XBTs) is developed. This new method uses a forward-stepping calculation which incorporates all of the forces on the XBT devices during their descent. Of particular note are drag forces which are calculated using a new drag coefficient expression. That expression, obtained entirely from computational fluid dynamic modeling, accounts for local variations in the ocean environment. Consequently, the method allows for accurate determination of depths for any local temperature environment. The results, which are entirely based on numerical simulation, are compared with the experiments of LM Sippican T-5 XBT probes. It is found that the calculated depths differ by less than 3% from depth estimates using the standard fall-rate equation (FRE). Furthermore, the differences decrease with depth. The computational model allows an investigation of the fluid flow patterns along the outer surface of the probe as well as in the interior channel. The simulations take account of complex flow phenomena such as laminar-turbulent transition and flow separation.

Annual cycles of chlorophyll- a , non-algal suspended particulate matter, and turbidity observed from space and in-situ in coastal waters - F. Gohin (Ocean Science, 2011)
Sea surface temperature, chlorophyll, and turbidity are three variables of the coastal environment commonly measured by monitoring networks. The observation networks are often based on coastal stations, which do not provide a sufficient coverage to validate the model outputs or to be used in assimilation over the continental shelf. Conversely, the products derived from satellite reflectance generally show a decreasing quality shoreward, and an assessment of the limitation of these data is required. The annual cycle, mean, and percentile 90 of the chlorophyll concentration derived from MERIS/ESA and MODIS/NASA data processed with a dedicated algorithm have been compared to in-situ observations at twenty-six selected stations from the Mediterranean Sea to the North Sea. Keeping in mind the validation, the forcing, or the assimilation in hydrological, sediment-transport, or ecological models, the non-algal Suspended Particulate Matter (SPM) is also a parameter which is expected from the satellite imagery. However, the monitoring networks measure essentially the turbidity and a consistency between chlorophyll, representative of the phytoplankton biomass, non-algal SPM, and turbidity is required. In this study, we derive the satellite turbidity from chlorophyll and non-algal SPM with a common formula applied to in-situ or satellite observations. The distribution of the satellite-derived turbidity exhibits the same main statistical characteristics as those measured in-situ, which satisfies the first condition to monitor the long-term changes or the large-scale spatial variation over the continental shelf and along the shore. For the first time, climatologies of turbidity, so useful for mapping the environment of the benthic habitats, are proposed from space on areas as different as the southern North Sea or the western Mediterranean Sea, with validation at coastal stations.

N/P ratio of nutrient uptake in the Baltic Sea - Z. Wan; L. Jonasson; H. Bi (Ocean Science, 2011)
The N/P ratio of nutrient uptake, the change of dissolved inorganic nitrogen (DIN) relative to the change of dissolved inorganic phosphorus (DIP), is a key parameter for many ecological models. In the Baltic Sea ecosystem, the N/P ratio of nutrient uptake varies among different basins and different seasons. The N/P ratio of nutrient alteration, i.e., the ratio of DIN to DIP altered before and after spring blooms, is not the same as the N/P ratio of nutrient uptake, but the former can be regarded as an indicator of the latter in the Baltic Sea. Based on the observed N/P ratio of nutrient alteration, we hypothesize a non-Redfield N/P ratio of nutrient uptake. The 3-D-ecosystem model ERGOM coupled with the circulation model DMI-BSHcmod was used to test this hypothesis. When the Redfield ratio was used in the model, the DIP surplus after spring blooms was too high and resulted in excessive growth of cyanobacteria and too much nitrogen fixation. When the non-Redfield ratio was used in the model, these problems tended to disappear. In summary, we show that: (1) the Redfield N/P ratio of nutrient uptake in the Baltic Sea tends to be too high; (2) a N/P ratio of 10:1 appears to work better than the Redfield value; and (3) the N/P ratio of nutrient uptake in the Baltic Proper during spring blooms is around 6:1. As the model limitation using one identical value for two N/P ratios for nutrient uptake and remineralization, the quantitative conclusions are only convincing as a model parameter even though it obviously improves model predictions. Whether this model parameter is consistent with the biological nutrient uptake is worth being further verified with some laboratory investigations or simulations using a more sophisticated model with independent N/P ratios for nutrient uptake and remineralization.

ENSO-correlated fluctuations in ocean bottom pressure and wind-stress curl in the North Pacific - D. P. Chambers (Ocean Science, 2011)
We examine the output of an ocean model forced by ECMWF winds to study the theoretical relationship between wind-induced changes in ocean bottom pressure in the North Pacific between 1992 until 2010 and ENSO. Our analysis indicates that while there are significant fluctuations correlated with some El Niño and La Niña events, the correlation is still relatively low. Moreover, the ENSO-correlated variability explains only 50 % of the non-seasonal, low-frequency variance. There are significant residual fluctuations in both wind-stress curl and ocean bottom pressure in the region with periods of 4-years and longer. One such fluctuation began in late 2002 and has been observed by the Gravity Recovery and Climate Experiment (GRACE). Even after accounting for possible ENSO-correlated variations, there is a significant trend in ocean bottom pressure in the region, equivalent to 0.7 ± 0.3 cm yr−1 of sea level from January 2003 until December 2008, which is confirmed with steric-corrected altimetry. Although this low-frequency fluctuation does not appear in the ocean model, we show that ECMWF winds have a significantly reduced trend that is inconsistent with satellite observations over the same time period, and so it appears that the difference is due to a forcing error in the model and not an intrinsic error.

The Aegean sea marine security decision support system - L. Perivoliotis; G. Krokos; K. Nittis; G. Korres (Ocean Science, 2011)
As part of the integrated ECOOP (European Coastal Sea Operational observing and Forecasting System) project, HCMR upgraded the already existing standalone Oil Spill Forecasting System for the Aegean Sea, initially developed for the Greek Operational Oceanography System (POSEIDON), into an active element of the European Decision Support System (EuroDeSS). The system is accessible through a user friendly web interface where the case scenarios can be fed into the oil spill drift model component, while the synthetic output contains detailed information about the distribution of oil spill particles and the oil spill budget and it is provided both in text based ECOOP common output format and as a series of sequential graphics. The main development steps that were necessary for this transition were the modification of the forcing input data module in order to allow the import of other system products which are usually provided in standard formats such as NetCDF and the transformation of the model's calculation routines to allow use of current, density and diffusivities data in z instead of sigma coordinates. During the implementation of the Aegean DeSS, the system was used in operational mode in order to support the Greek marine authorities in handling a real accident that took place in North Aegean area. Furthermore, the introduction of common input and output files by all the partners of EuroDeSS extended the system's interoperability thus facilitating data exchanges and comparison experiments.

Eddy characteristics in the northern South China Sea as inferred from Lagrangian drifter data (Ocean Science, 2011)
Cyclonic and anticyclonic eddies from large scale to submesoscale in the northern South China Sea (NSCS) have been statistically characterized based on the satellite-tracked Lagrangian drifters using our developed geometric eddy identification method. There are in total 2208 eddies identified, 70% of which are anticyclonic eddies. If the submesoscale eddies are eliminated, the other eddies in the NSCS will show a 1.2:1 ratio of the number of anticyclones (210) to the number of cyclones (171). The spatial distribution of the eddies is regional: in southwest of Taiwan, the number of anticyclones dominates the number of cyclones, and most of them are the submesoscale anticyclones with small radii; in contrast, the large and medium cyclonic eddies are a little more than the same scale anticyclonic eddies in northwest of Luzon. The temporal distribution of eddy number in the NSCS has a close relation with the Asian monsoon. The number of the large and medium eddies peaks during the winter monsoon, while the submesoscale eddies are apt to generate in the summer monsoon. The spatial and temporal patterns have a good agreement with the results of the sea surface height anomaly (SSHA). The maximum and mean tangential velocities of anticyclones (cyclones) are 40 (30) cm s−1 and 25 (15) cm s−1, respectively. The calculated normalized vorticities from drifters suggest that although the mesoscale eddies may be considered in geostrophic balance, ageostrophic dynamics and centrifugal effects may play an important role for the growth and decay of the mesoscale cores.

About uncertainties in practical salinity calculations - M. Le Menn (Ocean Science, 2011)
In the current state of the art, salinity is a quantity computed from conductivity ratio measurements, with temperature and pressure known at the time of the measurement, and using the Practical Salinity Scale algorithm of 1978 (PSS-78). This calculation gives practical salinity values S. The uncertainty expected in PSS-78 values is ±0.002, but no details have ever been given on the method used to work out this uncertainty, and the error sources to include in this calculation. Following a guide published by the Bureau International des Poids et Mesures (BIPM), using two independent methods, this paper assesses the uncertainties of salinity values obtained from a laboratory salinometer and Conductivity-Temperature-Depth (CTD) measurements after laboratory calibration of a conductivity cell. The results show that the part due to the PSS-78 relations fits is sometimes as significant as the instrument's. This is particularly the case with CTD measurements where correlations between variables contribute mainly to decreasing the uncertainty of S, even when expanded uncertainties of conductivity cell calibrations are for the most part in the order of 0.002 mS cm−1. The relations given here, and obtained with the normalized GUM method, allow a real analysis of the uncertainties' sources and they can be used in a more general way, with instruments having different specifications.

Development of Black Sea nowcasting and forecasting system - G. K. Korotaev; T. Oguz; V. L. Dorofeyev; S. G. Demyshev; A. I. Kubryakov; Yu. B. Ratner (Ocean Science, 2011)
The paper presents the development of the Black Sea community nowcasting and forecasting system under the Black Sea GOOS initiative and the EU framework projects ARENA, ASCABOS and ECOOP. One of the objectives of the Black Sea Global Ocean Observing System project is a promotion of the nowcasting and forecasting system of the Black Sea, in order to implement the operational oceanography in the Black Sea region. The first phase in the realization of this goal was the development of the pilot nowcasting and forecasting system of the Black Sea circulation in the framework of project ARENA funded by the EU. The ARENA project included the implementation of advanced modeling and data assimilation tools for near real time prediction. Further progress in development of the Black Sea nowcasting and forecasting system was made in the frame of ASCABOS project, which was targeted at strengthening the communication system, ensuring flexible and operative infrastructure for data and information exchange between the Black Sea partners and end-users. The improvement of the system was made in the framework of the ECOOP project. As a result it was transformed into a real-time mode operational nowcasting and forecasting system. The paper provides the general description of the main parts of the system: circulation and ecosystem models, data assimilation approaches, the system architecture as well as their qualitative and quantitative calibrations.

An eddy resolving tidal-driven model of the South China Sea assimilating along-track SLA data using the EnOI - J. Xie; F. Counillon; J. Zhu; L. Bertino (Ocean Science, 2011)
The upper ocean circulation in the South China Sea (SCS) is driven by the Asian monsoon, the Kuroshio intrusion through the Luzon Strait, strong tidal currents, and a complex topography. Here, we demonstrate the benefit of assimilating along-track altimeter data into a nested configuration of the HYbrid Coordinate Ocean Model that includes tides. Including tides in models is important because they interact with the main circulation. However, assimilation of altimetry data into a model including tides is challenging because tides and mesoscale features contribute to the elevation of ocean surface at different time scales and require different corrections. To address this issue, tides are filtered out of the model output and only the mesoscale variability is corrected with a computationally cheap data assimilation method: the Ensemble Optimal Interpolation (EnOI). This method uses a running selection of members to handle the seasonal variability and assimilates the track data asynchronously. The data assimilative system is tested for the period 1994–1995, during which time a large number of validation data are available. Data assimilation reduces the Root Mean Square Error of Sea Level Anomalies from 9.3 to 6.9 cm and improves the representation of the mesoscale features. With respect to the vertical temperature profiles, the data assimilation scheme reduces the errors quantitatively with an improvement at intermediate depth and deterioration at deeper depth. The comparison to surface drifters shows an improvement of surface current by approximately −9% in the Northern SCS and east of Vietnam. Results are improved compared to an assimilative system that does not include tides and a system that does not consider asynchronous assimilation.

Spectrophotometric high-precision seawater pH determination for use in underway measuring systems - S. Aßmann; C. Frank; A. Körtzinger (Ocean Science, 2011)
Autonomous sensors are required for a comprehensive documentation of the changes in the marine carbon system and thus to differentiate between its natural variability and anthropogenic impacts. Spectrophotometric determination of pH – a key variable of the seawater carbon system – is particularly suited to achieve precise and drift-free measurements. However, available spectrophotometric instruments are not suitable for integration into automated measurement systems (e.g. FerryBox) since they do not meet the major requirements of reliability, stability, robustness and moderate cost. Here we report on the development and testing of a~new indicator-based pH sensor that meets all of these requirements. This sensor can withstand the rough conditions during long-term deployments on ships of opportunity and is applicable to the open ocean as well as to coastal waters with a complex matrix and highly variable conditions. The sensor uses a high resolution CCD spectrometer as detector connected via optical fibers to a custom-made cuvette designed to reduce the impact of air bubbles. The sample temperature can be precisely adjusted (25 °C ± 0.006 °C) using computer-controlled power supplies and Peltier elements thus avoiding the widely used water bath. The overall setup achieves a measurement frequency of 1 min−1 with a precision of ±0.0007 pH units, an average offset of +0.0005 pH units to a reference system, and an offset of +0.0081 pH units to a certified standard buffer. Application of this sensor allows monitoring of seawater pH in autonomous underway systems, providing a key variable for characterization and understanding of the marine carbon system.

Seasonal and inter-annual temperature variability in the bottom waters over the western Black Sea shelf - G. I. Shapiro; F. Wobus; D. L. Aleynik (Ocean Science, 2011)
Long-term changes in the state of the Bottom Shelf Water (BSW) on the Western shelf of the Black Sea are assessed using analysis of intra-seasonal and inter-annual temperature variations. For the purpose of this study the BSW is defined as such shelf water mass between the seabed and the upper mixed layer (bounded by the σθ = 14.2 isopycnal) which has limited ability to mix vertically with oxygen-rich surface waters during the warm season due to formation of a seasonal pycnocline. A long-term time series of temperature anomalies in the BSW is constructed from in-situ observations taken over the 2nd half of the 20th century. The BSW is shown to occupy nearly half of the shelf area during the summer stratification period (May–November).The results reveal a warm phase in the 1960s/70s, followed by a cold phase between 1985 and 1995 and a further warming after 1995. The transition between the warm and cold periods coincides with a regime shift in the Black Sea ecosystem. While it was confirmed that the memory of winter convection is well preserved over the following months in the deep sea, the signal of winter cooling in the BSW significantly reduces during the warm season. The potential of the BSW to ventilate horizontally during the warm season with the deep-sea waters is assessed using isopycnic analysis of temperature variations. It is shown that temperature in the BSW is stronger correlated with the temperature of Cold Intermediate Waters (CIW) in the deep sea than with the severity of the previous winters, thus indicating that the isopycnal exchanges with the deep sea are more important for inter-annual/inter-decadal variability of the BSW on the western Black Sea shelf than effects of winter convection on the shelf itself.