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: Halocline water modification and along slope advection at the Laptev Sea continental margin Ocean Science Discussions, 10, 1581-1617, 2013 Author(s): D. Bauch, S. Torres-Valdes, I. Polyakov, A. Novikhin, I. Dmitrenko, J. McKay, and A. Mix A general pattern in water mass distribution and potential shelf-basin exchanges is revealed at the Laptev Sea continental slope based on hydrochemical and stable oxygen isotope data from summers 2005–2009. Despite considerable interannual variations, a frontal system can be inferred between shelf, continental slope and central Eurasian Basin waters in the upper 100 m of the water column along the continental slope. Net sea-ice melt is consistently found at the continental slope; however the sea-ice meltwater signal is independent from the local retreat of the ice cover and appears to be advected from upwind locations. In addition to the along-slope frontal system at the continental shelf break a strong gradient is identified on the Laptev Sea shelf between 122 and 126° E with an eastward increase of riverine and sea-ice related brine water contents. These waters cross the shelf break at ~ 140° E and feed the Low Salinity Halocline Water (LSHW, salinity S 〈 33) in the upper 50 m of the water column. Extremely high silicate concentrations in Laptev Sea bottom waters may lead to speculation on a link to the local silicate maximum found within the salinity range of ~ 33 to 34.5, typical for the Lower Halocline Water (LHW) at the continental slope. But brine signatures and nutrient ratios from the central Laptev Sea differ from those at the continental slope. Thus a significant contribution of Laptev Sea bottom waters to the LHW at the continental slope can be excluded. The silicate maximum within the LHW at the continental slope may be formed locally or at the outer Laptev Sea shelf. Similar to the advection of the sea-ice melt signal along the Laptev Sea continental slope the nutrient signal at 50–70 m water depth within the LHW might also be fed by advection parallel to the slope. Thus, our analyses suggest that advective processes from upwind locations play a significant role in the halocline formation in the northern Laptev Sea.

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

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: The effects of climate change on persistent organic pollutants (POPs) in the North Sea Ocean Science Discussions, 10, 1525-1557, 2013 Author(s): K. O'Driscoll, B. Mayer, J. Su, and M. Mathis The fate and cycling of two selected legacy persistent organic pollutants (POPs), PCB 153 and γ-HCH, in the North Sea in the 21st century have been modelled with combined hydrodynamic and fate and transport ocean models. To investigate the impact of climate variability on POPs in the North Sea in the 21st century, future scenario model runs for three 10 yr periods to the year 2100 using plausible levels of both in situ concentrations and atmospheric, river and open boundary inputs are performed. Since estimates of future concentration levels of POPs in the atmosphere, oceans and rivers are not available, our approach was to reutilise 2005 values in the atmosphere, rivers and at the open ocean boundaries for every year of the simulations. In this way, we attribute differences between the three 10 yr simulations to climate change only. For the HAMSOM and atmospheric forcing, results of the IPCC A1B (SRES) 21st century scenario are utilised, where surface forcing is provided by the REMO downscaling of the ECHAM5 global atmospheric model, and open boundary conditions are provided by the MPIOM global ocean model. Dry gas deposition and volatilisation of γ-HCH increase in the future relative to the present. In the water column, total mass of γ-HCH and PCB 153 remain fairly steady in all three runs. In sediment, γ-HCH increases in the future runs, relative to the present, while PCB 153 in sediment decreases exponentially in all three runs, but even faster in the future, both of which are the result of climate change. Annual net sinks exceed sources at the ends of all periods.

Description: Observations of phytoplankton spring bloom onset triggered by a density front in NW Mediterranean Ocean Science Discussions, 10, 1559-1580, 2013 Author(s): A. Olita, S. Sparnocchia, S. Cusí, L. Fazioli, R. Sorgente, J. Tintoré, and A. Ribotti Phytoplankton bloom in NW Mediterranan sea is a seasonal event that mainly occurrs in a limited area (Gulf of Lyon and Provençal basin) where this phenomenon is promoted by a cyclonic circulation, strong wind-driven mixing and subsequent spring restratification. At the southern boundary of this area a density front (North Balearic Front) separating denser waters from the lighter Modified Atlantic Waters reservoir at south is suspected to trigger weaker and earlier (late winter) blooms by (a) enhanced pumping of nutrients into the euphotic layer and (b) promoting an early restratification of the water column (by frontal instabilities). A multisensor glider round trip, equipped with CTD and fluorimeter, crossing the frontal area in February–March 2013, allowed to observe the bloom triggering after the decrease of intense wind-driven turbulent convection and mixing. Satellite imagery supports and confirms in-situ observations. It was shown that frontal activity has a relevant role in the promotion and acceleration of the dynamical restratification, with a consequent biological response in terms of primary production. Restratification is necessary preconditioning factor for bloom triggering in frontal area, net of other involved mechanism promoting the bloom as the enhanced biological pump. So, like for high-latitude fronts (Taylor and Ferrari, 2011a), also for this mid-latitude oligotrophic region front seems to promote new production by dynamically enahnced restratification inhibiting mixing. Finally, we argued that Sverdrup's Critical Depth criterion seems to apply in the northern well-mixed area, where the zeroing of heat fluxes (and related turbulent convection) does not correspond to a prompt onset of the bloom (which appeared 1 month later).

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

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: Global representation of tropical cyclone-induced ocean thermal changes using Argo data – Part 2: Estimating air–sea heat fluxes and ocean heat content changes Ocean Science Discussions, 11, 2907-2937, 2014 Author(s): L. Cheng, J. Zhu, and R. L. Sriver We use Argo temperature data to examine changes in ocean heat content (OHC) and air–sea heat fluxes induced by tropical cyclones (TC)s on a global scale. A footprint technique that analyzes the vertical structure of cross-track thermal responses along all storm tracks during the period 2004–2012 is utilized (see part I). We find that TCs are responsible for 1.87 PW (11.05 W m −2 when averaging over the global ocean basin) of heat transfer annually from the global ocean to the atmosphere during storm passage (0–3 days) on a global scale. Of this total, 1.05 ± 0.20 PW (4.80 ± 0.85 W m −2 ) is caused by Tropical storms/Tropical depressions (TS/TD) and 0.82 ± 0.21 PW (6.25 ± 1.5 W m −2 ) is caused by hurricanes. Our findings indicate that ocean heat loss by TCs may be a substantial missing piece of the global ocean heat budget. Net changes in OHC after storm passage is estimated by analyzing the temperature anomalies during wake recovery following storm events (4–20 days after storm passage) relative to pre-storm conditions. Results indicate the global ocean experiences a 0.75 ± 0.25 PW (5.98 ± 2.1W m −2 ) net heat gain annually for hurricanes. In contrast, under TS/TD conditions, ocean experiences 0.41 ± 0.21 PW (1.90 ± 0.96 W m −2 ) net ocean heat loss, suggesting the overall oceanic thermal response is particularly sensitive to the intensity of the event. The net ocean heat uptake caused by all storms is 0.34 PW.

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: Dynamics of turbulent western boundary currents at low latitude in a shallow water model Ocean Science Discussions, 11, 2461-2493, 2014 Author(s): C. Q. C. Akuetevi and A. Wirth The dynamics of low latitude turbulent western boundary currents crossing the equator is considered using numerical results from integrations of a reduced gravity shallow-water model. For viscosity values of 1000 m 2 s −1 and more, the boundary layer dynamics compares well to the analytical Munk-layer solution. When the viscosity is reduced, the boundary layer becomes turbulent and coherent structures in form of anticyclonic eddies, bursts (violent detachments of the viscous sub-layer) and dipoles appear. Three distinct boundary layers emerge, the viscous sub-layer, the advective boundary layer and the extended boundary layer. The first is characterized by a dominant vorticity balance between the viscous transport and the advective transport of vorticity. The second by a balance between the advection of planetary vorticity and the advective transport of relative vorticity. The extended boundary layer is the area to which turbulent motion from the boundary extends. The scaling of the three boundary layer thicknesses with viscosity is evaluated. Characteristic scales of the dynamics and dissipation are determined. A pragmatic approach to determine the eddy viscosity diagnostically for coarse resolution numerical models is proposed.

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: Calculating the water and heat balances of the Eastern Mediterranean basin using ocean modelling and available meteorological, hydrological, and ocean data Ocean Science Discussions, 8, 1301-1338, 2011 Author(s): M. Shaltout and A. Omstedt This paper analyses the Eastern Mediterranean water and heat balances over a 52-yr period. The modelling uses a process-oriented approach resolving the one-dimensional equations of momentum, heat, and salt conservation, with turbulence modelled using a two-equation model. The exchange through the Sicily Channel connecting the Eastern and Western basins is calculated from satellite altimeter data. The results illustrates that calculated surface temperature and salinity follow the reanalysed data well and with biases of −0.4 °C and −0.004, respectively. Monthly and yearly temperature and salinity cycles are also satisfactory simulated. Reanalysed data and calculated water mass structure and heat balance components are in good agreement, indicating that the air-sea interaction and the turbulent mixing are realistically simulated. The study illustrates that the water balance in the Eastern Mediterranean basin is controlled by the difference between inflows/outflows through the Sicily Channel and by the net precipitation rates. The heat balance is controlled by the heat loss from the water surface, sun radiation into the sea, and heat flow through the Sicily Channel, the first two displaying both climate trends. An annual net heat loss of approximately 9 W m −2 was balanced by net heat in flow through the Sicily Channel.

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

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: Tracer distribution in the Pacific Ocean following a release off Japan – what does an oceanic general circulation model tell us? Ocean Science Discussions, 8, 1441-1466, 2011 Author(s): H. Dietze and I. Kriest In the aftermath of an earthquake and tsunami on 11 March 2011 considerable amounts of radioactive materials were accidentally released into the sea off Fukushima-Daiichi, Japan. This study uses a three-dimensional eddy-resolving oceanic general circulation model to explore potential pathways of a tracer, similar to 137 Cs, from the coast to the open ocean. Results indicate that enhanced concentrations meet a receding spring bloom offshore and that the area of enhanced concentrations offshore is strongly determined by surface mixed layer dynamics. However, huge uncertainties remain. Among them are the realism of the simulated cross-shelf transport and apparently inconsistent estimates of the particle reactivity of 137 Cs which are discussed in a brief literature review. We argue that a comprehensive set of 137 Cs measurements, including sites offshore, could be a unique opportunity to both evaluate and advance the evaluation of oceanic general circulation models.

Description: Usefulness of high resolution coastal models for operational oil spill forecast: the Full City accident Ocean Science Discussions, 8, 1467-1504, 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 decision support system (DeSS) 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 is 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 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 including an analysis based on a higher resolution model (1.5 km resolution) for the area the model system show results that compare well with observations. The study also shows that an ensemble using three different models is useful when predicting/analyzing oil spill in coastal areas.

Description: Impact of combining GRACE and GOCE gravity data on ocean circulation estimates Ocean Science Discussions, 8, 1535-1573, 2011 Author(s): T. Janjić, J. Schröter, R. Savcenko, W. Bosch, A. Albertella, R. Rummel, and O. Klatt In this work we examine the impact of assimilation of multi-mission-altimeter data and the GRACE/GOCE gravity fields into the finite element ocean model (FEOM), with the focus on the Southern Ocean circulation. In order to do so, we use the geodetic approach for obtaining the dynamical ocean topography (DOT), that combines the multi-mission-altimeter data and the GRACE/GOCE gravity fields, and requires that both fields be spectrally consistent. The spectral consistency is achieved by filtering of the sea surface height and the geoid using profile approach. Combining the GRACE and GOCE data, a considerably shorter filter length resolving more DOT details can be used. In order to specify the spectrally consistent geodetic DOT we applied the Jekeli-Wahr filter corresponding to 241 km, 121 km, 97 km and 81 km halfwidths for the GRACE/GOCE based gravity field model GOCO01S and to the sea surface. More realistic features of the ocean assimilation were obtained in the Weddel gyre area due to increased resolution of the data fields, particularly for temperature field at the 800 m depth compared to Argo data.

Description: Wind forcing of salinity anomalies in the Denmark Strait overflow Ocean Science Discussions, 8, 1403-1440, 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, with the entire width of the array (70 km) freshening 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. We show that the OCCAM 1/12° Ocean General Circulation Model not only reproduces the 2004 freshening event ( r =0.96, p

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: Mesoscale variability of water masses in the Arabian Sea as revealed by ARGO floats Ocean Science Discussions, 8, 1369-1402, 2011 Author(s): X. Carton and P. L'Hegaret By analysing ARGO float data over the last four years, some aspects of the mesoscale variability of water masses in the Arabian Sea are described. The Red Sea Water outflow is strong in the Southwestern Gulf of Aden, in particular when a cyclonic gyre predominates in this region. Salinities of 36.5 and temperatures of 16 °C are found there between 600 and 1000 m depths. The Red Sea Water is more dilute in the eastern part of the Gulf, and fragments of this water mass can be advected offshore across the gulf or towards its northern coast by the regional gyres. The Red Sea Water outflow is also detected along the northeastern coast of Socotra, and fragments of RSW are found between one and three degrees of latitude north of this island. In the whole Gulf of Aden, the correlation between the deep motions of the floats and the SSH measured by altimetry is strong, at regional scale. The finer scale details of the float trajectories are more often related to the anomalous water masses that they encounter. The Persian Gulf Water (PGW) is found in the float profiles near Ras ash Sharbatat (near 57° E, 18° N), again with 36.5 in salinity and about 18–19 °C in temperature. These observations were achieved in winter when the southwestward monsoon currents can advect PGW along the South Arabian coast. Fragments of PGW are found in the Arabian Sea between 18 and 20° N and 63 and 65° E, showing that this water mass can escape the Gulf of Oman southeastward, in particular during summer.

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

Description: N/P ratio of nutrient uptake in the Baltic Sea Ocean Science Discussions, 8, 1233-1259, 2011 Author(s): Z. Wan, L. Jonasson, and H. Bi The N/P ratio of nutrient uptake, i.e., the ratio of dissolved inorganic nitrogen (DIN) to dissolved inorganic phosphorus (DIP) taken by primary producers, varies in different basins and in different seasons in the Baltic Sea. The N/P ratio of nutrient alteration fore and after spring blooms is not same as the N/P ratio of nutrient uptake, but the former can be regarded as an indicator for the later 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 3D-ecosystem model ERGOM coupled with the circulation model DMI-BSHcmod was used to test the hypothesis. When the Redfield ratio was used in the model, the DIP surplus after spring blooms was too high and resulted in the overly growth of cyanobacteria and too much nitrogen fixation. When the non-Redfield ratio was used in the model, the corresponding problem 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 lower N/P ratio 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.

Description: Spectrophotometric high-precision seawater pH determination for use in underway measuring systems Ocean Science Discussions, 8, 1339-1367, 2011 Author(s): S. Aßmann, C. Frank, and A. Kötzinger 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 opportunities and is applicable on the open ocean as well as in coastal waters with complex background 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 and an average offset of +0.0018 pH units to a pH reference during shipboard operation. Application of this sensor allows monitoring of seawater pH in autonomous underway systems, providing a key variable for characterization and understanding the marine carbon system.

Description: A pre-operational 3-D variational data assimilation system in the North/Baltic Sea Ocean Science Discussions, 8, 1131-1160, 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. The univariate analysis for 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) of the background error. Observations of temperature and salinity ( T / S ) profiles in the North/Baltic Sea are assimilated in the year of 2005. Effectiveness of the data assimilation scheme is assessed by comparison with the control run that no assimilation is done. The statistical analysis indicates that the model simulation is significantly improved with the 3DVAR scheme. On average, the root mean square error (RMSE) of temperature and salinity is 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: 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: Numerical simulation and decomposition of kinetic energies in the Central Mediterranean Sea: insight on mesoscale circulation and energy conversion Ocean Science Discussions, 8, 1161-1214, 2011 Author(s): R. Sorgente, A. Olita, P. Oddo, L. Fazioli, and A. Ribotti The spatial and temporal variability of eddy and mean kinetic energy of the Central Mediterranean Sea has been investigated, from January 2008 to December 2010, by mean of a numerical simulation mainly to quantify the mesoscale dynamics and their relationships with physical forcing. In order to understand the energy redistribution processes, the baroclinic energy conversion has been analysed, suggesting hypotheses about the drivers of the mesoscale activity in this area. The ocean model used is based on the Princeton Ocean Model implemented at 1/32° horizontal resolution. Surface momentum and buoyancy fluxes are interactively computed by mean of standard bulk formulae using predicted model Sea Surface Temperature and atmospheric variables provided by the European Centre for Medium Range Weather Forecast operational analyses. At its lateral boundaries the model is one-way nested within the Mediterranean Forecasting System operational products. The model domain has been subdivided in four sub-regions: Sardinia channel and southern Tyrrhenian Sea, Sicily channel, eastern Tunisian shelf and Libyan Sea. Temporal evolution of eddy and mean kinetic energy has been analysed, on each of the four sub-regions composing the model domain, showing different behaviours. On annual scales and within the first 5 m depth, the eddy kinetic energy represents approximately the 60 % of the total kinetic energy over the whole domain, confirming the strong mesoscale nature of the surface current flows in this area. The analyses show that the model well reproduces the path and the temporal behaviour of the main known sub-basin circulation features. New mesoscale structures have been also identified, from numerical results and direct observations, for the first time as the Pantelleria Vortex and the Medina Gyre. The classical the kinetic energy decomposition (eddy and mean) allowed to depict and to quantify the stable and fluctuating parts of the circulation in the region, and to differentiate the four sub-regions as function of relative and absolute strength of the mesoscale activity. Furthermore the Baroclinic Energy Conversion term shows that in the Sardinia Channel the mesoscale activity, due to baroclinic instabilities, is significantly larger than in the other sub-regions, while a negative sign of the energy conversion, meaning a transfer of energy from the Eddy Kinetic Energy to the Eddy Available Potential Energy, has been recorded only for the surface layers of the Sicily Channel during summer.

Description: Modelling of the anthropogenic tritium transient and its decay product helium-3 in the Mediterranean Sea using a high-resolution regional model Ocean Science Discussions, 11, 2691-2732, 2014 Author(s): M. Ayache, J.-C. Dutay, P. Jean-Baptiste, K. Beranger, T. Arsouze, J. Beuvier, J. Palmieri, B. Le-vu, and W. Roether This numerical study provides the first simulation of the anthropogenic tritium invasion and its decay product helium-3 ( 3 He) in the Mediterranean Sea. The simulation covers the entire tritium ( 3 H) transient generated by the atmospheric nuclear-weapon tests performed in the 1950s and early 1960s and run till 2011. Tritium, helium-3 and their derived age estimates are particularly suitable for studying intermediate and deep-water ventilation and spreading of water masses at intermediate/deep levels. The simulation is made using a high resolution regional model NEMO-MED12 forced at the surface with prescribed tritium evolution derived from observations. The simulation is compared to measurements of tritium and helium-3 performed along large-scale transects in the Mediterranean Sea during the last few decades on cruises of Meteor M5/6, M31/1, M44/4, M51/2, M84/3, and Poseidon 234. The results show that the input function used for the tritium, generates a realistic distribution of the main hydrographic features of the Mediterranean Sea circulation. In the eastern basin, the results highlight the weak formation of Adriatic Deep Water in the model, which explains its weak contribution to the Eastern Mediterranean Deep Water in the Ionian sub-basin. It produces a realistic representation of the Eastern Mediterranean Transient signal, simulating a deep-water formation in the Aegean sub-basin at the beginning of the 1993, with a realistic timing of deep-water renewal in the eastern basin. In the western basin, the unusual intense deep convection event of winter 2005 in the Gulf of Lions during the Western Mediterranean Transition is simulated. However the spreading of the recently ventilated deep water toward the South is too weak. The ventilation and spreading of the Levantine Intermediate Water from the eastern basin toward the western basin is simulated with realistic tracer-age distribution compared to observation-based estimates.

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: The Aegean Sea marine security decision support system Ocean Science Discussions, 8, 1025-1053, 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 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: How well can we derive Global Ocean Indicators from Argo data? Ocean Science Discussions, 8, 999-1024, 2011 Author(s): K. von Schuckmann and P.-Y. Le Traon Argo deployments began in the year 2000 and by November 2007 the array was 100 % complete, covering the global ocean from the surface down to 2000 m depth. In this study, Argo temperature and salinity measurements during the period 2005 to 2010 are used to develop a revised estimation of Global Ocean Indicators (GOIs) such as heat content variability, freshwater content and steric height. These revised indices are based on a simple box averaging scheme using a weighted mean. They include a proper estimation of the errors due to data handling methods and climatology uncertainties. A global ocean heat content change (OHC) trend of 0.55 ± 0.1 W m −2 is estimated over the time period 2005–2010. Similarly, a global steric sea level (GSSL) rise of 0.69 ± 0.14 mm yr −1 is observed. The global ocean freshwater content (OFC) trend is barely significant. Results show that there is significant interannual variability at global scale, especially for global OFC. Annual mean GOIs from the today's Argo samling can be derived with an accuracy of ±0.10 cm for GSSL, ±0.21 × 10 8 J m −2 for global OHC, and ±700 km 3 for global OFC. Long-term trends (15 yr) of GOIs based on the complete Argo sampling (10–1500 m depth) can be performed with an accuracy of about ±0.03 mm yr −1 for steric rise, ±0.02 W m −2 for ocean warming and ±20 km 3 yr −1 for global OFC trends – under the assumption that no systematic errors remain in the observing system.

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

Description: Quality control of automated hyperspectral remote sensing measurements from a seaborne platform Ocean Science Discussions, 8, 613-638, 2011 Author(s): S. P. Garaba, M. R. Wernand, and O. Zielinski In this study four data quality flags are presented for automated and unmanned above-water hyperspectral optical measurements collected underway in the North Sea, The Minch, Irish Sea and Celtic Sea in April/May 2009. Coincident to these optical measurements a DualDome D12 (Mobotix, Germany) camera system was used to capture sea surface and sky images. The first three flags are based on meteorological conditions, to select erroneous incoming solar irradiance ( E S ) taken during dusk, dawn, before significant incoming solar radiation could be detected or under rainfall. Furthermore, the relative azimuthal angle of the optical sensors to the sun is used to identify possible sunglint free sea surface zones. A total of 629 spectra remained after applying the meteorological masks (first three flags). Based on this dataset, a fourth flag for sunglint was generated by analysing and evaluating water leaving radiance ( L W ) and remote sensing reflectance ( R RS ) spectral behaviour in the presence and absence of sunglint salient in the simultaneously available sea surface images. Spectra conditions satisfying "mean L W (700–950 nm) 〈 2 mW m −2 nm −1 Sr −1 " or alternatively "minimum R RS (700–950 nm) 〈 0.010 Sr −1 ", mask the most measurements affected by sunglint, providing efficient flagging of sunglint in automated quality control. It is confirmed that valid optical measurements can be performed 0° ≤ Φ ≤ 360° although 90° ≤ Φ ≤ 135° is recommended.

Description: Numerical modelling of POC yearly dynamics in the southern Baltic under variable scenarios of nutrients, light and temperature Ocean Science Discussions, 8, 675-700, 2011 Author(s): L. Dzierzbicka-Glowacka, K. Kuliński, A. Maciejewska, J. Jakacki, and J. Pempkowiak This paper presents various scenarios of the particulate organic carbon (POC) in the southern Baltic Sea. The study is based on a one-dimensional Particulate Organic Carbon model (1-D POC). Mathematically, the pelagic variables of 1-D POC model are described by a second-order partial differential equations of the diffusion type with biogeochemical sources and sinks. The POC concentration is determined as the sum of phytoplankton, zooplankton and dead organic matter (detritus) concentrations. The temporal changes in the phytoplankton biomass are caused by primary production, mortality, grazing by zooplankton and sinking. The zooplankton biomass is affected by ingestion, excretion, faecal production, mortality, and carnivorous grazing. The changes in the pelagic detritus concentration are determined by input of: dead phytoplankton and zooplankton, natural mortality of predators, faecal pellets, and sinks: sedimentation, zooplankton grazing and biochemical decomposition. The 1-D POC model was used to simulate temporal dynamics of POC in the southern Baltic Sea (Gdansk Deep, Bornholm Deep and Gotland Deep) under scenarios characterized by different temperature, nutrients and light. Daily, monthly, seasonal and annual variabilities of POC in the upper water layer are presented for the different scenarios. The starting-point of the numerical simulations was assumed as average values of the investigated pelagic variables for 1965–1998 period. Two- to three-fold increases of POC concentrations in late spring were revealed as well as the shift towards postponed maximum POC concentration. It is speculated that, due to POC increase, oxygenation of under-halocline water layer will decrease, while supply of food to organisms from higher trophic level should increase.

Description: A computational method for determining XBT depths Ocean Science Discussions, 8, 1777-1802, 2011 Author(s): J. Stark, J. Gorman, M. Hennessey, F. Reseghetti, J. Willis, J. Lyman, and J. Abraham A new technique for determining the depth of expendable bathythermographs (XBTs) is developed. This new method combines 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 an experimental descent of an LM-Sippican T-5 XBT. It is found that the calculated depths differ by less than 3 % from depth estimates using the industry standard FRE. Furthermore, the differences decrease with depth. The computational model allows an investigation of the fluid 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: 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: A vertical-mode decomposition to investigate low-frequency internal motion across the Atlantic at 26° N Ocean Science Discussions, 8, 2047-2100, 2011 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 that span 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). On both 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. At the boundaries, the transport fluctuations described by the first baroclinic mode represent less than 10% of the variance of the full transport signal. At the eastern boundary, a linear combination of many baroclinic modes is required to explain the observed vertical density profile of the seasonal cycle, a result that is consistent with the oceanic response to wind-forcing not propagating far from the eastern boundary.

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

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: Assessment of a physical-biogeochemical coupled model system for operational service in the Baltic Sea Ocean Science Discussions, 9, 835-876, 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 of 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: Numerical tools to estimate the flux of a gas across the air-water interface and assess the heterogeny of its forcing functions Ocean Science Discussions, 9, 909-975, 2012 Author(s): V. M. N. de C. da S. Vieira A numerical tool was developed for the estimation of gas fluxes across the air water interface. The primary objective is to use it to estimate CO 2 fluxes. Nevertheless application to other gases is easily accomplished by changing the values of the parameters related to the physical properties of the gases. A user friendly software was developed allowing to build upon a standard kernel a custom made gas flux model with the preferred parametrizations. These include single or double layer models; several numerical schemes for the effects of wind in the air-side and water-side transfer velocities; the effect of turbulence from current drag with the bottom; and the effects on solubility of water temperature, salinity, air temperature and pressure. It was also developed an analysis which decomposes the difference between the fluxes in a reference situation and in alternative situations into its several forcing functions. This analysis relies on the Taylor expansion of the gas flux model, requiring the numerical estimation of partial derivatives by a multivariate version of the collocation polynomial. Both the flux model and the difference decomposition analysis were tested with data taken from surveys done in the lagoonary system of Ria Formosa, south Portugal, in which the CO 2 fluxes were estimated using the IRGA and floating chamber method whereas the CO 2 concentrations were estimated using the IRGA and degasification chamber. Observations and estimations show a remarkable fit.

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

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: Towards high resolution mapping of 3-D mesoscale dynamics from observations: preliminary comparison of retrieval techniques and models within MESCLA project Ocean Science Discussions, 9, 1045-1083, 2012 Author(s): B. Buongiorno Nardelli, S. Guinehut, A. Pascual, Y. Drillet, S. Ruiz, and S. Mulet Within the MyOcean R&D project MESCLA (MEsoSCale dynamical Analysis through combined model, satellite and in situ data), different estimates of the vertical velocities derived from observations have been compared. Two main approaches have been considered, one based on the retrieval of 3-D fields from the observations alone and one based on the analyses provided by MyOcean MERCATOR models. The motivation for this double approach is that, while data assimilation in numerical models is crucial to obtain more accurate analyses and forecasts, its results might be significantly influenced by specific model configurations (e.g. forcing, parameterization of smaller scale processes and spatial resolution). On the other hand, the purely observation-based approach is limited by the underlying assumptions of simplified dynamical models and by the relatively low resolution of present products. MESCLA tested innovative methods for the high resolution mapping of 3-D mesoscale dynamics from observations, developing new products that might be used to gradually build the next generations of operational observation-based products.

Description: On the Shelf Resonances of the Gulf of Carpentaria and the Arafura Sea Ocean Science Discussions, 9, 443-497, 2012 Author(s): D. J. Webb A numerical model is used to investigate the resonances of the Gulf of Carpentaria and the Arafura Sea. The model is forced at the shelf edge, first with physically realistic real values of angular velocity. The response functions at points within the region show maxima and other behaviour which imply that resonances are involved but it is difficult to be more specific. 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 responses 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: 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: NEMO on the shelf: assessment of the Iberia-Biscay-Ireland configuration Ocean Science Discussions, 10, 83-151, 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, and P. Marsaleix The Iberia-Biscay-Ireland (IBI) system serves one of the 7 MyOcean "Monitoring and Forecasting Centres". A high resolution simulation covering the IBI region is set-up over July 2007–February 2009. The NEMO (Nucleus for European Modelling of the Ocean) model is used with a 1/36° horizontal resolution and 50 z-levels in the vertical. New developments have been incorporated in NEMO to make it suitable to open- as well as coastal-ocean modelling. In this paper, we pursue three main objectives: (1) give an overview of the model configuration used for the simulations; (2) give a broad-brush account of one particular aspect of this work, namely consistency verification; this type of validation is conducted upstream of the implementation of the system before it is used for production and routinely validated; it is meant to guide model development in identifying gross deficiencies in the modelling of several key physical processes; (3) show that such a regional modelling system has potential as a complement to patchy observations (an integrated approach) to give information on non-observed physical quantities and to provide links between observations by identifying broader-scale patterns and processes. We concentrate on the year 2008. We first provide domain-wide consistency verification results in terms of barotropic tides, transports, sea surface temperature and stratification. We then focus on two dynamical sub-regions: the Celtic shelves and the Bay of Biscay slope and deep regions. The model-data consistency is checked for variables and processes such as tidal currents, tidal fronts, internal tides, residual elevation. We also examine the representation in the model of a seasonal pattern of the Bay of Biscay circulation: the warm extension of the Iberian Poleward Current along the northern Spanish coast (Navidad event) in winter 2007–2008.

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: Impact of the Indonesian throughflow on Agulhas leakage Ocean Science Discussions, 10, 353-391, 2013 Author(s): D. Le Bars, H. A. Dijkstra, and W. P. M. De Ruijter Using ocean models of different complexity we show that opening the Indonesian Passage between the Pacific and the Indian Ocean increases the input of Indian Ocean water into the South Atlantic via the Agulhas leakage. In a strongly eddying global ocean model this response results from an increased Agulhas Current transport and a constant proportion of Agulhas retroflection south of Africa. The leakage increases through an increased frequency of ring shedding events. In an idealized two-layer and flat-bottom eddy resolving model, the proportion of the Agulhas Current transport that retroflects is (for a wide range of wind stress forcing) not affected by an opening of the Indonesian Passage. A linear ocean model is not able to explain this behavior which reveals the importance of mixed barotropic/baroclinic instabilities in controlling the Agulhas leakage.

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

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: On the use of the Strouhal/Stokes number to explain the dynamics and water column structure on shelf seas Ocean Science Discussions, 9, 3723-3738, 2012 Author(s): A. J. Souza In recent years coastal oceanographers have suggested the use of the "Strouhal" number or it's inverse the "Stokes" number, which have been defined as the ratios of the frictional depth (δ) to the water column depth ( h ) or vice versa, to describe the effect of bottom boundary layer turbulence on the vertical structure of both density and currents. Although they have mention that the effects of rotation should be important, they have tended to omit it. This omission may be important when talking about tidal currents as the frictional depth from a fully cyclonic to a fully anticyclonic tidal ellipse can vary up to an order of magnitude in the mid latitudes; so that the stokes number might appear smaller (larger) than it is resulting in frictional effects being underestimated (overestimated). Here a way to calculate a Stokes number, in which the effect of the Earth's rotation is taken into account, is suggested. Then the standard Stokes and the rotational Stokes numbers are used as predictors for the position of the tidal mixing fronts in the Irish Sea. Results show that the rotational number improves prediction of the front in shallow cyclonic areas of the eastern Irish Sea. This suggest that the effect of rotation on the water column structure will be more important in shallow shelf seas and estuaries with strong rotational currents.

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: Chaotic variability of the meridional overturning circulation on subannual to interannual timescales Ocean Science Discussions, 9, 3191-3238, 2012 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 internal 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 (1 Sv = 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 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 accordance with 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 despite the amplitude of the chaotic variability being larger on subannual than on interannual timescales, the ratio of chaotic to total MOC variability is larger on interannual than on subannual timescales.

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: Co-existence of wind seas and swells along the west coast of India during non-monsoon season Ocean Science Discussions, 9, 3097-3125, 2012 Author(s): R. Rashmi, V. M. Aboobacker, P. Vethamony, and M. P. John Wave data collected along the west coast of India (off Goa, Ratnagiri and Dwarka) during non-monsoon season have been analysed to study the co-existence of wind seas and swells. Diurnal variation in wind and wave parameters is noticeable along the central west coast of India (off Goa and Ratnagiri), and this is not present along the northwest coast of India (off Dwarka). Swells are predominantly mature (91%) and old (88%) during late pre-monsoon and post-monsoon seasons, respectively. Sea Swell Energy Ratio quantifies wind sea, swell and mixed seas prevailing in the regions during non-monsoon season. Intermodal Distance (ID) between the energy peaks is moderately separated during non-monsoon season, whereas, during the shamal events, energy peaks are very close to each other (ID ~ 0). However, pure wind seas (ID ~ 1) are found to co-exist with the swells during non-monsoon season. Wind seas are growing, when wind and wind seas are opposite to swell direction. Wind seas have minimum angular spreads in multimodal state. Under low winds, the interaction between wind sea and swell dominates and thereby the multimodal state reduces to unimodal state. The fetch available for the evolution of the wind sea spectrum has been estimated, and it is found to be less than 150 km. For the fetch limited condition, a non-dimensional empirical relation has been derived relating the significant wind sea height in terms of wind speed and peak wind sea period, and this relation fits for the west coast of India.

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: Assimilation of SLA along track observations in the Mediterranean with an oceanographic model forced by atmospheric pressure Ocean Science Discussions, 9, 1577-1598, 2012 Author(s): S. Dobricic, C. Dufau, P. Oddo, N. Pinardi, I. Pujol, and M.-H. Rio A large number of SLA observations at a high along track horizontal resolution are an important ingredient of the data assimilation in the Mediterranean Forecasting System (MFS). Recently new higher frequency SLA products have become available, and the atmospheric pressure forcing has been implemented in the numerical model used in the MFS data assimilation system. In a set of numerical experiments we show that in order to obtain the most accurate analyses the ocean model should include the atmospheric pressure forcing and the observations should contain the atmospheric pressure signal. When the model is not forced by the atmospheric pressure the high frequency filtering of SLA observations, however, improves the quality of the analyses. It is further shown that MFS analyses, produced by an assimilation system given by the numerical model and the high frequency SLA observations, have a correct power spectrum at high wave numbers and they filter efficiently the SLA assimilated observations which, on the other hand, are contaminated by high wavenumber noise.

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

Description: An operational model for the West Iberian coast: products and services Ocean Science Discussions, 9, 1651-1689, 2012 Author(s): M. Mateus, G. Riflet, P. Chambel, L. Fernandes, R. Fernandes, M. Juliano, F. Campuzano, H. de Pablo, and R. Neves The paper presents the structure and application of a regional scale operational modelling tool for the West Iberian coast. The forecasting suite includes nested hydrodynamic models forced with up-to-date meteorological forecast data and large scale model results as lateral boundary conditions. The present status of the system and its recent upgrades are reviewed, offering a general description of the main components of the system: circulation model, qualitative and quantitative validation methodology and type of results. Seasonal differences in temperature, salinity and current velocity fields are illustrated with model results, and the validation shows a satisfactory reproduction of the top and deep layers thermodynamics. The system provides boundary forcing for a number of local scale model applications via downscaling of the solution, and its potential for products and services for both scientific and coastal management activities is discussed here.

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: Investigation of saline water intrusions into the Curonian Lagoon (Lithuania) and two-layer flow in the Klaipėda Strait using finite element hydrodynamic model Ocean Science Discussions, 10, 321-352, 2013 Author(s): P. Zemlys, C. Ferrarin, G. Umgiesser, S. Gulbinskas, and D. Bellafiore This work is focused on the application of a modelling system to simulate 3-D interaction between the Curonian Lagoon and the Baltic Sea coastal waters and to reflect spatio-temporal dynamics of marine waters in the Curonian Lagoon. The model system is based on the finite element program package SHYFEM which can be used to resolve the hydrodynamic equations in lagoons, coastal seas, estuaries and lakes. The results of a one year 3-D model simulation with real weather and hydrological forcing show that the saline water intrusions from the sea through Klaipėda Strait are gradually decreasing with distance from the sea and become negligible (average annual salinity about 0.5 ‰) at a~distance of about 20 km to the south of Kiaulės Nugara island. Analyses of the simulation results also show this area being highly heterogeneous according to the vertical salinity distribution. While in the deeper Klaipėda Strait (harbour waterway) differences in average salinity between near bottom and surface layers varies in the range 2–2.5 ‰, in the rest of the Curonian Lagoon it is less than 0.1 ‰. Analyses of the simulation results confirmed the presence of a two-directional flow that from time to time changes to either saline water one-directional flow to the Curonian Lagoon or fresh water one-directional flow to the sea. Two-directional flow duration decreases with a distance from sea entrance in Klaipėda Strait from around 180 days yr −1 close to the sea entrance to 50 days yr −1 just behind Kiaulės Nugara island. One-directional outflow duration is increasing with a distance from the sea entrance from 100 to 225 days yr −1 . One-directional inflow duration occurs in the range 85–100 days yr −1 . The analysis of the ratio of buoyancy layer thickness to water depth ( h b /H) and the Wedderburn number showed three main flow regimes in the strait, identifying the main importance of wind action in the along-strait direction. Absence of wind or cross-strait wind regimes allow the maintenance of an two-layer flow typical of estuarine dynamics.

Description: Variability in the air–sea interaction patterns and time-scales within the Southeastern Bay of Biscay, as observed by HF radar data Ocean Science Discussions, 9, 2793-2815, 2012 Author(s): A. Fontán, G. Esnaola, J. Sáenz, and M. González Two high frequency (HF) radar stations were installed on the Southeastern Bay of Biscay in 2009, providing high spatial and temporal resolution and large spatial coverage currents for the first time in the area. This has enabled to determine quantitatively the air–sea interaction patterns and time-scales for the period 2009–2010. The analysis was conducted by using the Barnett-Preisendorfer approach to canonical correlation analysis (CCA) of reanalysis surface winds and HF radar-derived currents. The results reveal that the CCA yields two canonical patterns. The first wind-current interaction pattern corresponds to the classical Ekman drift at sea surface, whilst the second describes an anticyclonic/cyclonic surface circulation. The results obtained demonstrate that the local winds play an important role in driving the upper water circulation. The wind-current interaction time-scales are mainly related to diurnal breezes and synoptic variability. In particular, the breezes force diurnal currents in the continental shelf and slope of the Southeastern 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: 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: MERIS-based ocean colour classification with the discrete Forel–Ule scale Ocean Science Discussions, 9, 2817-2849, 2012 Author(s): M. R. Wernand, A. Hommersom, and H. J. van der Woerd Multispectral information from satellite borne ocean colour sensors is at present used to characterize natural waters via the retrieval of concentration of the three dominant optical constituents; pigments of phytoplankton, non-algal particles and coloured dissolved organic matter. A limitation of this approach is that accurate retrieval of these constituents requires detailed local knowledge of the specific absorption and scattering properties. In addition, the retrieval algorithms generally use only a limited part of the collected spectral information. In this paper we present an additional new algorithm that has the merit to use the full spectral information in the visible domain to characterize natural waters in a simple and globally valid way. This Forel–Ule MERIS (FUME) algorithm converts the normalized multi-band reflectance information into a discrete set of numbers using uniform colourimetric functions. The Forel–Ule scale is a sea colour comparator scale that has been developed to cover all possible natural sea colours, ranging from indigo blue (the open ocean) to brownish-green (coastal water) and even brown (humic-acid dominated) waters. Data using this scale have been collected since the late nineteenth century, and therefore, this algorithm creates the possibility to compare historic ocean colour data with present-day satellite ocean colour observations. The FUME algorithm was tested by transforming a number of MERIS satellite images into Forel–Ule colour index images and comparing in situ observed FU numbers with FU numbers modelled from in situ radiometer measurements.

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: Assimilation of sea-ice concentration in a global climate model – physical and statistical aspects Ocean Science Discussions, 9, 2403-2455, 2012 Author(s): S. Tietsche, D. Notz, J. H. Jungclaus, and J. Marotzke We investigate the initialization 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 significantly reduces assimilation error both in identical-twin experiments and when assimilating sea-ice observations, reducing the concentration error by a factor of four to six, and the thickness error by a factor of two. 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 the strong dependence of thermodynamic ice growth on ice concentration necessitates an adjustment of mean ice thickness in the analysis update. 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 proportional mean-thickness updates are superior to the other two methods considered and enable us to assimilate sea ice in a global climate model using simple Newtonian relaxation.

Description: Sea level variability in the Arctic Ocean observed by satellite altimetry Ocean Science Discussions, 9, 2375-2401, 2012 Author(s): P. Prandi, M. Ablain, A. Cazenave, and N. Picot We investigate sea level variability in the Arctic Ocean from observations. Variability estimates are derived both at the basin scale and on smaller local spatial scales. The periods of the signals studied vary from high frequency (intra-annual) to long term trends. We also investigate the mechanisms responsible for the observed variability. Different data types are used, the main one being a recent reprocessing of satellite altimetry data in the Arctic Ocean. Satellite altimetry data is compared to tide gauges measurements, steric sea level derived from temperature and salinity fields and GRACE ocean mass estimates. We establish a consistent regional sea level budget over the GRACE availability era (2003–2009) showing that the sea level drop observed by altimetry over this period is driven by ocean mass loss rather than steric effects. The comparison of altimetry and tide gauges time series show that the two techniques are in good agreement regarding sea level trends. Coastal areas of high variability in the altimetry record are also consistent with tide gauges records. An EOF analysis of September mean altimetry fields allows identifying two regions of wind driven variability in the Arctic Ocean: the Beaufort Gyre region and the coastal European and Russian Arctic. Such patterns are related to atmospheric regimes through the Arctic Oscillation and Dipole Anomaly.