ALBERT

Description: The Mediterranean Sea system: a review and an introduction to the special issue Ocean Science, 9, 789-803, 2013 Author(s): T. Tanhua, D. Hainbucher, K. Schroeder, V. Cardin, M. Álvarez, and G. Civitarese 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.

Description: Corrigendum to "NEMO on the shelf: assessment of the Iberia–Biscay–Ireland configuration" published in Ocean Sci., 9, 745–771, 2013 Ocean Science, 9, 787-787, 2013 Author(s): 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, and the Mercator Research and Development Team No abstract available.

Description: ENSO components of the Atlantic multidecadal oscillation and their relation to North Atlantic interannual coastal sea level anomalies Ocean Science, 9, 535-543, 2013 Author(s): J. Park and G. Dusek 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.

Description: From the chlorophyll a in the surface layer to its vertical profile: a Greenland Sea relationship for satellite applications Ocean Science, 9, 431-445, 2013 Author(s): A. Cherkasheva, E.-M. Nöthig, E. Bauerfeind, C. Melsheimer, and A. Bracher 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 C tot (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.

Description: Large-scale temperature and salinity changes in the upper Canadian Basin of the Arctic Ocean at a time of a drastic Arctic Oscillation inversion Ocean Science, 9, 447-460, 2013 Author(s): P. Bourgain, J. C. Gascard, J. Shi, and J. Zhao 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.

Description: Seasonality of intermediate waters hydrography west of the Iberian Peninsula from an 8 yr semiannual time series of an oceanographic section Ocean Science, 9, 411-429, 2013 Author(s): E. Prieto, C. González-Pola, A. Lavín, R. F. Sánchez, and M. Ruiz-Villarreal 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.

Description: Variability in the air–sea interaction patterns and timescales within the south-eastern Bay of Biscay, as observed by HF radar data Ocean Science, 9, 399-410, 2013 Author(s): A. Fontán, G. Esnaola, J. Sáenz, and M. González 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.

Description: Chaotic variability of the meridional overturning circulation on subannual to interannual timescales Ocean Science, 9, 805-823, 2013 Author(s): J. J.-M. Hirschi, A. T. Blaker, B. Sinha, A. Coward, B. de Cuevas, S. Alderson, and G. Madec 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 = 10 6 m 3 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.

Description: Transport of Antarctic bottom water through the Kane Gap, tropical NE Atlantic Ocean Ocean Science, 9, 825-835, 2013 Author(s): E. G. Morozov, R. Y. Tarakanov, and H. van Haren 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 (Θ 〈 1.9 °C) based on the mooring and LADCP (Lowered Acoustic Doppler Current Profiler) data varies approximately within ±0.35 Sv in the northern and southern directions. The annual mean AABW transport through the Kane Gap is almost zero.

Description: Changes in extreme regional sea surface height due to an abrupt weakening of the Atlantic meridional overturning circulation Ocean Science, 10, 881-891, 2014 Author(s): S.-E. Brunnabend, H. A. Dijkstra, M. A. Kliphuis, B. van Werkhoven, H. E. Bal, F. Seinstra, J. Maassen, and M. van Meersbergen As an extreme scenario of dynamical sea level changes, regional sea surface height (SSH) changes that occur in the North Atlantic due to an abrupt weakening of the Atlantic meridional overturning circulation (AMOC) are simulated. Two versions of the same ocean-only model are used to study the effect of ocean model resolution on these SSH changes: a high-resolution (HR) strongly eddying version and a low-resolution (LR) version in which the effect of eddies is parameterised. The weakening of the AMOC is induced in both model versions by applying strong freshwater perturbations around Greenland. A rapid decrease of the AMOC in the HR version induces much shorter return times of several specific regional and coastal extremes in North Atlantic SSH than in the LR version. This effect is caused by a change in main eddy pathways associated with a change in separation latitude of the Gulf Stream.

Description: A century of sea level data and the UK's 2013/14 storm surges: an assessment of extremes and clustering using the Newlyn tide gauge record Ocean Science, 10, 1031-1045, 2014 Author(s): M. P. Wadey, I. D. Haigh, and J. M. Brown For the UK's longest and most complete sea level record (Newlyn), we assess extreme high waters and their temporal clustering; prompted by the 2013/2014 winter of storms and flooding. These are set into context against this almost 100-year record. We define annual periods for which storm activity and high sea levels can be compared on a year-by-year basis. Amongst the storms and high tides which affected Newlyn, the recent winter produced the largest recorded high water level (3 February 2014) and five other high water events above a 1 in 1-year return period. The large magnitude of tide and mean sea level, and the close inter-event spacings (of large return period high waters), suggests that the 2013/2014 extreme high water level "season" can be considered the most extreme on record. However, storm and sea level events may be classified in different ways. For example, in the context of sea level rise (which we calculate linearly as 1.81 ± 0.1 mm yr −1 from records between 1915 to 2014), a lower probability combination of surge and tide occurred on 29 January 1948, whilst the 1995/1996 storm surge season saw the most high waters of ≥ the 1 in 1-year return period. We provide a basic categorisation of the four types of extreme high water level cluster, ranging from consecutive tidal cycles to multiple years. The assessment is extended to other UK sites (with shorter sea level records and different tide-surge characteristics), which suggests 2013/2014 was particularly unusual. Further work will assess clustering mechanisms and flood system "memory".

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

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

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

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

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

Description: Influence of frontal cyclone evolution on the 2009 (Ekman) and 2010 (Franklin) Loop Current eddy detachment events Ocean Science, 10, 947-965, 2014 Author(s): Y. S. Androulidakis, V. H. Kourafalou, and M. Le Hénaff The anticyclonic Loop Current Eddy (LCE) shedding events are strongly associated with the evolution of Loop Current Frontal Eddies (LCFEs) over the eastern Gulf of Mexico (GoM). A numerical simulation, in tandem with in situ measurements and satellite data, was used to investigate the Loop Current (LC) evolution and the surrounding LCFE formation, structure, growth and migration during the Eddy Ekman and Eddy Franklin shedding events in the summers of 2009 and 2010, respectively. During both events, northern GoM LCFEs appeared vertically coherent to at least 1500 m in temperature observations. They propagated towards the base of the LC, where, together with the migration of Campeche Bank (southwest GoM shelf) eddies from south of the LC, contributed to its "necking-down". Growth of Campeche Bank LCFEs involved in Eddy Franklin was partially attributed to Campeche Bank waters following upwelling events. Slope processes associated with such upwelling included offshore exports of high positive potential vorticity that may trigger cyclone formation and growth. The advection and growth of LCFEs, originating from the northern and southern GoM, and their interaction with the LC over the LCE detachment area favor shedding conditions and may contribute to the final separation of the LCE.

Description: The Rossby radius in the Arctic Ocean Ocean Science, 10, 967-975, 2014 Author(s): A. J. G. Nurser and S. Bacon The first (and second) baroclinic deformation (or Rossby) radii are presented north of ~60° N, focusing on deep basins and shelf seas in the high Arctic Ocean, the Nordic seas, Baffin Bay, Hudson Bay and the Canadian Arctic Archipelago, derived from climatological ocean data. In the high Arctic Ocean, the first Rossby radius increases from ~5 km in the Nansen Basin to ~15 km in the central Canadian Basin. In the shelf seas and elsewhere, values are low (1–7 km), reflecting weak density stratification, shallow water, or both. Seasonality strongly impacts the Rossby radius only in shallow seas, where winter homogenization of the water column can reduce it to below 1 km. Greater detail is seen in the output from an ice–ocean general circulation model, of higher resolution than the climatology. To assess the impact of secular variability, 10 years (2003–2012) of hydrographic stations along 150° W in the Beaufort Gyre are also analysed. The first-mode Rossby radius increases over this period by ~20%. Finally, we review the observed scales of Arctic Ocean eddies.

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

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

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

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

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

Description: Spectrophotometric high-precision seawater pH determination for use in underway measuring systems Ocean Science, 7, 597-607, 2011 Author(s): S. Aßmann, C. Frank, and A. Körtzinger 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.

Description: An eddy resolving tidal-driven model of the South China Sea assimilating along-track SLA data using the EnOI Ocean Science, 7, 609-627, 2011 Author(s): J. Xie, F. Counillon, J. Zhu, and L. Bertino 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.

Description: A computational method for determining XBT depths Ocean Science, 7, 733-743, 2011 Author(s): J. Stark, J. Gorman, M. Hennessey, F. Reseghetti, J. Willis, J. Lyman, J. Abraham, and M. Borghini 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.

Description: Mean Dynamic Topography of the Black Sea, computed from altimetry, drifter measurements and hydrology data Ocean Science, 7, 745-753, 2011 Author(s): A. A. Kubryakov and S. V. Stanichny 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.

Description: Towards a regional ocean forecasting system for the IBI (Iberia-Biscay-Ireland area): developments and improvements within the ECOOP project framework Ocean Science, 8, 143-159, 2012 Author(s): S. Cailleau, J. Chanut, J.-M. Lellouche, B. Levier, C. Maraldi, G. Reffray, and M. G. Sotillo 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.

Description: Technical Note: Detection of gas bubble leakage via correlation of water column multibeam images Ocean Science, 8, 175-181, 2012 Author(s): J. Schneider von Deimling and C. Papenberg 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.

Description: Black Sea coastal forecasting system Ocean Science, 8, 183-196, 2012 Author(s): A. I. Kubryakov, G. K. Korotaev, V. L. Dorofeev, Y. B. Ratner, A. Palazov, N. Valchev, V. Malciu, R. Matescu, and T. Oguz 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.

Description: Sea surface temperature anomalies, seasonal cycle and trend regimes in the Eastern Pacific coast Ocean Science, 8, 81-90, 2012 Author(s): A. Ramos-Rodríguez, D. B. Lluch-Cota, S. E. Lluch-Cota, and A. Trasviña-Castro 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 × 10 8 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.

Description: Assimilation of sea-ice concentration in a global climate model – physical and statistical aspects Ocean Science, 9, 19-36, 2013 Author(s): S. Tietsche, D. Notz, J. H. Jungclaus, and J. Marotzke 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.

Description: Toward a multivariate reanalysis of the North Atlantic Ocean biogeochemistry during 1998–2006 based on the assimilation of SeaWiFS chlorophyll data Ocean Science, 9, 37-56, 2013 Author(s): C. Fontana, P. Brasseur, and J.-M. Brankart 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.

Description: A clustering analysis of eddies' spatial distribution in the South China Sea Ocean Science, 9, 171-182, 2013 Author(s): J. Yi, Y. Du, X. Wang, Z. He, and C. Zhou 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.

Description: Sea surface freshening inferred from SMOS and ARGO salinity: impact of rain Ocean Science, 9, 183-192, 2013 Author(s): J. Boutin, N. Martin, G. Reverdin, X. Yin, and F. Gaillard 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 km 2 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.

Description: Microstructure observations during the spring 2011 STRATIPHYT-II cruise in the northeast Atlantic Ocean Science, 8, 945-957, 2012 Author(s): E. Jurado, H. A. Dijkstra, and H. J. van der Woerd 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, K T , range from 10 −5 to 10 0 m 2 s −1 in the weakly-stratified stations, and range from 3 × 10 −4 to 2 × 10 0 m 2 s −1 in the well-mixed stations. The turbulent kinetic energy dissipation rates, ε, range from 3 × 10 −8 to 2 × 10 −6 m 2 s −3 south of the transition zone, and from 10 −7 to 10 −5 m 2 s −3 north of the transition zone. The station-averaged K T 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 K T 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.

Description: Impact of the sea surface temperature forcing on hindcasts of Madden-Julian Oscillation events using the ECMWF model Ocean Science, 8, 1071-1084, 2012 Author(s): E. de Boisséson, M. A. Balmaseda, F. Vitart, and K. Mogensen 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.

Description: Imbalance of energy and momentum source terms of the sea wave transfer equation for fully developed seas Ocean Science, 8, 1085-1098, 2012 Author(s): G. V. Caudal 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.

Description: High frequency fluctuations in the heat content of an ocean general circulation model Ocean Science, 8, 813-825, 2012 Author(s): A. M. Huerta-Casas and D. J. Webb 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.

Description: A 20-year reanalysis experiment in the Baltic Sea using three-dimensional variational (3DVAR) method Ocean Science, 8, 827-844, 2012 Author(s): W. Fu, J. She, and M. Dobrynin 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.

Description: Evaluation of Release-05 GRACE time-variable gravity coefficients over the ocean Ocean Science, 8, 859-868, 2012 Author(s): D. P. Chambers and J. A. Bonin 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.

Description: The Mediterranean Ocean Colour Observing System – system development and product validation Ocean Science, 8, 869-883, 2012 Author(s): G. Volpe, S. Colella, V. Forneris, C. Tronconi, and R. Santoleri 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.

Description: The coherence of small island sea level with the wider ocean: a model study Ocean Science, 9, 111-119, 2013 Author(s): Joanne Williams and Chris W. Hughes 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.

Description: X-band COSMO-SkyMed wind field retrieval, with application to coastal circulation modeling Ocean Science, 9, 121-132, 2013 Author(s): A. Montuori, P. de Ruggiero, M. Migliaccio, S. Pierini, and G. Spezie 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.

Description: Net primary productivity, upwelling and coastal currents in the Gulf of Ulloa, Baja California, México Ocean Science, 8, 703-711, 2012 Author(s): E. González-Rodríguez, A. Trasviña-Castro, G. Gaxiola-Castro, L. Zamudio, and R. Cervantes-Duarte 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.

Description: Assessment of a physical-biogeochemical coupled model system for operational service in the Baltic Sea Ocean Science, 8, 683-701, 2012 Author(s): Z. Wan, J. She, M. Maar, L. Jonasson, and J. Baasch-Larsen 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.

Description: An operational model for the West Iberian coast: products and services Ocean Science, 8, 713-732, 2012 Author(s): M. Mateus, G. Riflet, P. Chambel, L. Fernandes, R. Fernandes, M. Juliano, F. Campuzano, H. de Pablo, and R. Neves 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.

Description: On the shelf resonances of the Gulf of Carpentaria and the Arafura Sea Ocean Science, 8, 733-750, 2012 Author(s): D. J. Webb 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.

Description: The analysis of large-scale turbulence characteristics in the Indonesian seas derived from a regional model based on the Princeton Ocean Model Ocean Science, 8, 615-631, 2012 Author(s): K. O'Driscoll and V. Kamenkovich 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, q 2 ; the turbulence master scale, ℓ; mixing coefficients of momentum, K M ; and temperature and salinity, K H ; 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 q 2 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 q 2 occur in the deep layer extending several hundred meters above the bottom. The observed increase of q 2 near the very bottom is probably due to the increase of velocity shear and the corresponding shear production of q 2 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 m 2 s −1 were obtained for K M and K H in deep water in the vicinity of the Lifamatola Sill. These estimates agree well with basin-scale averaged values of 13.3 × 10 −4 m 2 s −1 found diagnostically for K H in the deep Banda and Seram Seas (Gordon et al., 2003) and a value of 9.0 × 10 −4 m 2 s −1 found diagnostically for K H 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.

Description: TOPAZ4: an ocean-sea ice data assimilation system for the North Atlantic and Arctic Ocean Science, 8, 633-656, 2012 Author(s): P. Sakov, F. Counillon, L. Bertino, K. A. Lisæter, P. R. Oke, and A. Korablev 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.

Description: Atlantic transport variability at 25° N in six hydrographic sections Ocean Science, 8, 497-523, 2012 Author(s): C. P. Atkinson, H. L. Bryden, S. A. Cunningham, and B. A. King 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.

Description: Propagation and dissipation of internal tides in the Oslofjord Ocean Science, 8, 525-543, 2012 Author(s): A. Staalstrøm, E. Aas, and B. Liljebladh 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 cm 2 s −1 ) is more than four times higher than in the rest of the fjord (~5 cm 2 s −1 or less).

Description: Operational SAR-based sea ice drift monitoring over the Baltic Sea Ocean Science, 8, 473-483, 2012 Author(s): J. Karvonen 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.

Description: Impact of SLA assimilation in the Sicily Channel Regional Model: model skills and mesoscale features Ocean Science, 8, 485-496, 2012 Author(s): A. Olita, S. Dobricic, A. Ribotti, L. Fazioli, A. Cucco, C. Dufau, and R. Sorgente 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.

Description: Three-dimensional modelling of wave-induced current from the surf zone to the inner shelf Ocean Science, 8, 657-681, 2012 Author(s): H. Michaud, P. Marsaleix, Y. Leredde, C. Estournel, F. Bourrin, F. Lyard, C. Mayet, and F. Ardhuin 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.

Description: Particle aggregation at the edges of anticyclonic eddies and implications for distribution of biomass Ocean Science, 8, 389-400, 2012 Author(s): A. Samuelsen, S. S. Hjøllo, J. A. Johannessen, and R. Patel 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.

Description: Influence of Ross Sea Bottom Water changes on the warming and freshening of the Antarctic Bottom Water in the Australian-Antarctic Basin Ocean Science, 8, 419-432, 2012 Author(s): K. Shimada, S. Aoki, K. I. Ohshima, and S. R. Rintoul 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.

Description: Support to oil spill emergencies in the Bonifacio Strait, western Mediterranean Ocean Science, 8, 443-454, 2012 Author(s): A. Cucco, A. Ribotti, A. Olita, L. Fazioli, B. Sorgente, M. Sinerchia, A. Satta, A. Perilli, M. Borghini, K. Schroeder, and R. Sorgente 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.

Description: On the chemical dynamics of extracellular polysaccharides in the high Arctic surface microlayer Ocean Science, 8, 401-418, 2012 Author(s): Q. Gao, C. Leck, C. Rauschenberg, and P. A. Matrai 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, 〈 0.22 μm, obtained after filtration) samples were collected and chemically characterized from the SML and the corresponding SSW at an open lead centered at 87.5° N and 5° E. Total organic carbon was persistently enriched in the SML with a mean enrichment factor (EF) of 1.45 ± 0.41, whereas sporadic depletions of dissolved carbohydrates and amino acids were observed. Monosaccharide compositional analysis reveals that EPS in the Arctic lead was formed mainly of distinctive heteropolysaccharides, enriched in xylose, fucose and glucose. The mean concentrations of total hydrolysable neutral sugars in SSW were 94.9 ± 37.5 nM in high molecular weight (HMW) 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.

Description: On the influence of wind and waves on underwater irradiance fluctuations Ocean Science, 8, 455-471, 2012 Author(s): M. Hieronymi and A. Macke The influence of various wind and wave conditions on the variability of downwelling irradiance E d (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 E d 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 E d . 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.

Description: Transport of warm Upper Circumpolar Deep Water onto the western Antarctic Peninsula continental shelf Ocean Science, 8, 433-442, 2012 Author(s): D. G. Martinson and D. C. McKee 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.

Description: On the outflow of dense water from the Weddell and Ross Seas in OCCAM model Ocean Science, 8, 369-388, 2012 Author(s): R. Kerr, K. J. Heywood, M. M. Mata, and C. A. E. Garcia 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 ≡ 10 6 m 3 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.

Description: Transports and budgets in a 1/4 ° global ocean reanalysis 1989–2010 Ocean Science, 8, 333-344, 2012 Author(s): K. Haines, M. Valdivieso, H. Zuo, and V. N. Stepanov 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.

Description: 137 Cs off Fukushima Dai-ichi, Japan – model based estimates of dilution and fate Ocean Science, 8, 319-332, 2012 Author(s): H. Dietze and I. Kriest In the aftermath of an earthquake and tsunami on 11 March 2011 radioactive 137 Cs 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 137 Cs 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 m 2 s −1 . The cross-shelf diffusivity is at 500 ± 300 m 2 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 137 Cs, 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 137 Cs signal is flushed away. The latter estimates concern the direct release to the sea only.

Description: A vertical-mode decomposition to investigate low-frequency internal motion across the Atlantic at 26° N Ocean Science, 8, 345-367, 2012 Author(s): Z. B. Szuts, J. R. Blundell, M. P. Chidichimo, and J. Marotzke 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

Description: Wave climate in the Arkona Basin, the Baltic Sea Ocean Science, 8, 287-300, 2012 Author(s): T. Soomere, R. Weisse, and A. Behrens 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 H S =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.

Description: Laminar and weakly turbulent oceanic gravity currents performing inertial oscillations Ocean Science, 8, 301-317, 2012 Author(s): A. Wirth 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.

Description: How well can we derive Global Ocean Indicators from Argo data? Ocean Science, 7, 783-791, 2011 Author(s): K. von Schuckmann and P.-Y. Le Traon 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 × 10 8 J m −2 for GOHC, and ±700 km 3 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 km 3 yr −1 for GOFC – under the assumption that no systematic errors remain in the observing system.

Description: Wind forcing of salinity anomalies in the Denmark Strait overflow Ocean Science, 7, 821-834, 2011 Author(s): S. Hall, S. R. Dye, K. J. Heywood, and M. R. Wadley 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

Description: A pre-operational three Dimensional variational data assimilation system in the North/Baltic Sea Ocean Science, 7, 771-781, 2011 Author(s): S. Y. Zhuang, W. W. Fu, and J. She 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.

Description: About uncertainties in practical salinity calculations Ocean Science, 7, 651-659, 2011 Author(s): M. Le Menn 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.

Description: Eddy characteristics in the northern South China Sea as inferred from Lagrangian drifter data Ocean Science, 7, 661-669, 2011 Author(s): Jiaxun Li, Ren Zhang, and Baogang Jin 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.

Description: The Aegean sea marine security decision support system Ocean Science, 7, 671-683, 2011 Author(s): L. Perivoliotis, G. Krokos, K. Nittis, and G. Korres 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.

Description: N/P ratio of nutrient uptake in the Baltic Sea Ocean Science, 7, 693-704, 2011 Author(s): Z. Wan, L. Jonasson, and H. Bi 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.

Description: Annual cycles of chlorophyll- a , non-algal suspended particulate matter, and turbidity observed from space and in-situ in coastal waters Ocean Science, 7, 705-732, 2011 Author(s): F. Gohin 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.

Description: An empirical stochastic model of sea-surface temperatures and surface winds over the Southern Ocean Ocean Science, 7, 755-770, 2011 Author(s): S. Kravtsov, D. Kondrashov, I. Kamenkovich, and M. Ghil 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.

Description: Operational forecast of hydrophysical fields in the Georgian Black Sea coastal zone within the ECOOP Ocean Science, 7, 793-803, 2011 Author(s): A. A. Kordzadze and D. I. Demetrashvili 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.

Description: Usefulness of high resolution coastal models for operational oil spill forecast: the "Full City" accident Ocean Science, 7, 805-820, 2011 Author(s): G. Broström, A. Carrasco, L. R. Hole, S. Dick, F. Janssen, J. Mattsson, and S. Berger 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.

Description: Modeling of wave-induced irradiance variability in the upper ocean mixed layer Ocean Science, 8, 103-120, 2012 Author(s): M. Hieronymi, A. Macke, and O. Zielinski 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 E d 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 E d is still possible in 100 m depth.

Description: Towards an improved description of ocean uncertainties: effect of local anamorphic transformations on spatial correlations Ocean Science, 8, 121-142, 2012 Author(s): J.-M. Brankart, C.-E. Testut, D. Béal, M. Doron, C. Fontana, M. Meinvielle, P. Brasseur, and J. Verron 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.

Description: ENSO-correlated fluctuations in ocean bottom pressure and wind-stress curl in the North Pacific Ocean Science, 7, 685-692, 2011 Author(s): D. P. Chambers 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.

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

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

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

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

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

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

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

Description: Preface "Deep Ocean Exchange with the Shelf (DOES)" Ocean Science, 7, 101-109, 2011 Author(s): J. Johnson and P. Chapman No abstract available.

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

Description: Phytoplankton distribution and nitrogen dynamics in the southwest indian subtropical gyre and Southern Ocean waters Ocean Science, 7, 113-127, 2011 Author(s): S. J. Thomalla, H. N. Waldron, M. I. Lucas, J. F. Read, I. J. Ansorge, and E. Pakhomov During the 1999 Marion Island Oceanographic Survey (MIOS 4) in late austral summer, a northbound and reciprocal southbound transect were taken along the Southwest Indian and Madagascar Ridge, between the Prince Edward Islands and 31° S. The sections crossed a number of major fronts and smaller mesoscale features and covered a wide productivity spectrum from subtropical to subantarctic waters. Associated with the physical survey were measurements of size fractionated chlorophyll, nutrients and nitrogen (NO 3 , NH 4 and urea) uptake rates. Subtropical waters were characterised by low chlorophyll concentrations (max = 0.27.3 mg m −3 dominated by pico-phytoplankton cells (〉 81%) and very low f-ratios ( 〈 0.1), indicative of productivity based almost entirely on recycled ammonium and urea. Micro-phytoplankton growth was limited by the availability of NO 3 ( 〈 0.5 mmol m −3 and Si(OH) 4 ( 〈 1.5 mmol m −3 through strong vertical stratification preventing the upward flux of nutrients into the euphotic zone. Biomass accumulation of small cells was likely controlled by micro-zooplankton grazing. In subantarctic waters, total chlorophyll concentrations increased (max = 0.74 mg m −3 relative to the subtropical waters and larger cells became more prevalent, however smaller phytoplankton cells and low f-ratios ( 〈 0.14) still dominated, despite sufficient NO 3 availability. The results from this study favour Si(OH) 4 limitation, light-limited deep mixing and likely Fe deficiency as the dominant mechanisms controlling significant new production by micro-phytoplankton. The percentage of micro-phytoplankton cells and rates of new production did however increase at oceanic frontal regions (58.6% and 11.22%, respectively), and in the region of the Prince Edward archipelago (61.4% and 14.16%, respectively). Here, water column stabilization and local Fe-enrichment are thought to stimulate phytoplankton growth rates. Open ocean regions such as these provide important areas for local but significant particulate organic carbon export and biological CO 2 draw-down in an overall high nutrient low chlorophyll Southern Ocean.

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

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

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

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

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

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

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

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

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