ALBERT

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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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.

Description: Technical Note: Mean sea level variation in the Singapore Strait from long-term tide data Ocean Science Discussions, 9, 2255-2271, 2012 Author(s): P. Tkalich, M. T. Babu, and P. Vethamony Winds over the South China Sea (SCS) are primarily responsible for the observed variability in sea level anomalies (SLAs) in the Singapore Strait (SS). The present study focuses on remote forcing contributing to local mean sea level changes in the SS in seasonal and inter-annual scales, and relating the long term mean sea level variation to El Niño/ENSO. As Tanjong Pagar (TP) tide station in the SS has nearly 23.5 yr (1984–2007) of time series data with less data gaps, this data was subject to harmonic and sea level analyses. The mean sea level changes suggest that the fluctuations are quasi-periodic. Rising and falling of sea level is noticed atleast 7 times in a period of 15 yr, with 3 distinct sharp falls (1984–1987, 1989–1992 and 1995–1996) and 4 sharp rises (1987–1988, 1992–1993, 1994–1995 and 1997–1999). These sea level falls are related to El Niño events. When we segregated the results into 2 time spans, we find that from 1984 to 1999 the sea level was on the rising trend in spite of sharp falls, and from 1999 to 2007 on gradual falling trend. More or less similar trend was observed by other researchers for the SCS with altimetry data. During the El Niño periods of 1987 and 1992, the inter-annual MSL variability is the highest, of the order of 7 cm. In one of the events, sea level recovered from a fall of 60 mm (in 1987) to a rise of 40 mm (in 1988). During 1992 to 1999, sea level was continuously on rising trend (from −50 mm to +60 mm), except in one year (1995–1996). The analysis shows a MSL rise rate of 15.7 mm yr −1 , which is very closer to MSL in the SCS. The average rate of sea level rise around Singapore as shown by the Tanjong Pagar tidal station is 1.6 mm yr −1 , and this matches with the global sea level rise.

Description: Impact of the sea surface temperature forcing on hindcasts of Madden-Julian Oscillation events using the ECMWF model Ocean Science Discussions, 9, 2535-2559, 2012 Author(s): E. de Boisséson, M. A. Balmaseda, F. Vitart, and K. Mogensen This paper explores the sensitivity of the prediction of Madden Julian Oscillation (MJO) events to different aspects of the sea surface temperature (SST) in the European Centre for Medium-range Weather Forecasts (ECMWF) model. The impact of temporal resolution of SST on the MJO is first evaluated via a set of monthly hindcast experiments. The experiments are conducted with an atmosphere forced by persisted SST anomalies, monthly and weekly SSTs. Skill scores are clearly degraded when weekly SSTs are replaced by monthly values or persisted anomalies. The new high resolution OSTIA SST daily reanalysis would in principle allow to establish the impact of daily versus weekly SST values on the MJO prediction. It is found however that OSTIA SSTs provide lower skill scores than the weekly product. Further experiments show that this loss of skill cannot be attributed to either the mean state or the daily frequency of OSTIA SSTs. Additional diagnostics show that the phase relationship between OSTIA SSTs and tropical convection is not optimal with repspect to observations. Such result suggests that capturing the correct SST-convection phase relationship is a major challenge for the MJO predictions.

Description: Surface signature of Mediterranean water eddies in the North-East Atlantic: effect of the upper ocean stratification Ocean Science Discussions, 9, 2457-2491, 2012 Author(s): I. Bashmachnikov and X. Carton Meddies, intra-thermocline eddies of Mediterranean water, are often visible at the sea surface as positive sea-level anomalies. Here we study the surface signature of several meddies tracked with RAFOS floats and AVISO altimetry. Then, theoretical estimates of the surface signature of a meddy are derived, based on geostrophy and potential vorticity balance. The intensity of the surface signature is proportional to the meddy core radius and to the Coriolis parameter, and inversely proportional to the core depth and buoyancy frequency. This indicates that surface signature of a meddy may be strongly reduced by the upper ocean stratification. Estimates suggest that the southernmost limit for detection in altimetry of small meddies (with radii on the order of 15 km) should lie in the northern subtropics, while large meddies (with radii of 25–30 km) could be detected as far south as the northern tropics. During the initial period of meddy acceleration after meddy formation or a stagnation stage, a cyclonic signal also is generated at the sea-surface, but mostly the anticyclonic surface signal follows the meddy.

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

Description: Measurement of turbulence in the oceanic mixed layer using Synthetic Aperture Radar (SAR) Ocean Science Discussions, 9, 2851-2883, 2012 Author(s): S. G. George and A. R. L. Tatnall Turbulence in the surface layer of the ocean contributes to the transfer of heat, gas and momentum across the air-sea boundary. As such, study of turbulence in the ocean surface layer is becoming increasingly important for understanding its effects on climate change. Direct Numerical Simulation (DNS) techniques were implemented to examine the interaction of small-scale wake turbulence in the upper ocean layer with incident electromagnetic radar waves. Hydrodynamic-electromagnetic wave interaction models were invoked to demonstrate the ability of Synthetic Aperture Radar (SAR) to observe and characterise surface turbulent wake flows. A range of simulated radar images are presented for a turbulent surface current field behind a moving surface vessel, and compared with the surface flow fields to investigate the impact of turbulent currents on simulated radar backscatter. This has yielded insights into the feasibility of resolving small-scale turbulence with remote-sensing radar and highlights the potential for extracting details of the flow structure and characteristics of turbulence using SAR.

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

Description: Arctic rapid sea ice loss events in regional coupled climate scenario experiments Ocean Science Discussions, 9, 2327-2373, 2012 Author(s): R. Döscher and T. Koenigk Rapid sea ice loss events (RILEs) in a mini-ensemble of regional Arctic coupled climate model scenario experiments are analyzed. Mechanisms of sudden ice loss are strongly related to atmospheric circulation conditions and preconditioning by sea ice thinning during the seasons and years before the event. Clustering of events in time suggests a strong control by large scale atmospheric circulation. Anomalous atmospheric circulation is forcing ice flow and providing warm air affecting winter ice growth. Even without a seasonal preconditioning during winter, ice drop events can be initiated by anomalous inflow of warm air from the Atlantic sector during summer. It is shown that RILE events can be generated solely based on atmospheric circulation changes without possible competing mechanisms, such as anomalous seasonal radiative forcing or short-lived forcers (e.g. soot). Such forces do merely play minor roles or no role at all in our model. Mechanisms found are qualitatively in line with observations of the 2007 RILE.

Description: Evaluation of Release-05 GRACE time-variable gravity coefficients over the Ocean Ocean Science Discussions, 9, 2187-2214, 2012 Author(s): D. P. Chambers and J. A. Bonin The latest release of GRACE 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 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 sub-polar oceans, which is comparable to estimated uncertainty for the output from the ocean models.

Description: Sea-air CO 2 flux estimated from SOCAT surface-ocean CO 2 partial pressure data and atmospheric CO 2 mixing ratio data Ocean Science Discussions, 9, 2273-2326, 2012 Author(s): C. Rödenbeck, R. F. Keeling, D. C. E. Bakker, N. Metzl, A. Olsen, C. Sabine, and M. Heimann Surface-ocean CO 2 partial pressure data have been assimilated into a simple diagnostic model of surface-ocean biogeochemistry to estimate the spatio-temporal CO 2 partial pressure field and ultimately the sea-air CO 2 fluxes. Results compare well with the widely used monthly climatology by Takahashi et al. (2009) but also contain some short-term and interannual variations. Fitting the same model to atmospheric CO 2 data yields less robust but consistent estimates, confirming that using the partial pressure based estimates as ocean prior in atmospheric CO 2 inversions may improve land CO 2 flux estimates. Estimated seasonality of ocean-internal carbon sources and sinks is discussed in the light of observed nutrient variations.

Description: Long-term monitoring programme of the hydrological variability in the Mediterranean Sea: a first overview of the HYDROCHANGES network Ocean Science Discussions, 9, 1741-1812, 2012 Author(s): K. Schroeder, C. Millot, L. Bengara, S. Ben Ismail, M. Bensi, M. Borghini, G. Budillon, V. Cardin, L. Coppola, C. Curtil, A. Drago, B. El Moumni, J. Font, J. L. Fuda, J. García-Lafuente, G. P. Gasparini, H. Kontoyiannis, D. Lefevre, P. Puig, P. Raimbault, G. Rougier, J. Salat, C. Sammari, J. C. Sánchez Garrido, A. Sanchez-Roman, S. Sparnocchia, C. Tamburini, I. Taupier-Letage, A. Theocharis, M. Vargas-Yáñez, and A. Vetrano The long-term monitoring of basic hydrological parameters (temperature and salinity), collected as time series with adequate temporal resolution (i.e. with a sampling interval allowing the resolution of all important time scales) in key-places of the Mediterranean Sea (straits and channels, zones of dense water formation, deep parts of the basins), constitute a priority in the context of global changes. This led CIESM (The Mediterranean Science Commission) to support, since 2002, the HYDROCHANGES programme ( http://www.ciesm.org/marine/programs/hydrochanges.htm ), a network of autonomous conductivity-temperature-depth (CTD) sensors, deployed on mainly short and easily manageable sub-surface moorings. The HYDROCHANGES strategy is twofold and develops on different scales. To get information about long-term changes of hydrological characteristics, long time series are needed. But before these series are long enough they allow detecting links between them at shorter time scales that may provide extremely valuable information about the functioning of the Mediterranean Sea. Here we present the history of the programme and the current setup of the network (monitored sites, involved groups) and provide an overview of all the available time series, discussing some of the results obtained thanks to the programme.

Description: Toward a multivariate reanalysis of the North Atlantic ocean biogeochemistry during 1998–2006 based on the assimilation of SeaWiFS chlorophyll data Ocean Science Discussions, 9, 1887-1931, 2012 Author(s): C. Fontana, P. Brasseur, and J.-M. Brankart Today, the routine assimilation of satellite data into operational models of the 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 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 modelling. A simplified variant of the Kalman filter is used to assimilate ocean color data during a 9 year-long period. In this frame, two experiences 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, the nitrate World Ocean Atlas database and a derived climatology. Along the simulation period, the non-linear assimilation scheme clearly improves the surface chlorophyll concentrations analysis and forecast, 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 litterature. This feature is explained by the weak correlation taken into account by the assimilation between surface phytoplankton and nitrate concentration deeper than 50 m. The assessement of the non-linear assimilation experiments indicates that the proposed methodology provides the skeleton of an assimilative system suitable for reanalysing the ocean biogeochemistry based on ocean color data.

Description: Temporal and spatial distribution of the meiobenthic community in Daya Bay, South China Sea Ocean Science Discussions, 9, 1853-1885, 2012 Author(s): L. Tang, H. X. Li, and Y. Yan Spatial and temporal biodiversity patterns of the meiobenthos were studied for the first time in Daya Bay, which is a tropical semi-enclosed basin located in the South China Sea. The abundance, biomass, and composition of the meiobenthos and the basic environmental factors in the bay were investigated. The following 19 taxonomic groups were represented in the meiofauna: Nematoda, Copepoda, Polychaeta, Oligochaeta, Kinorhyncha, Gastrotricha, Ostracoda, Bivalvia, Turbellaria, Nemertinea, Sipuncula, Hydroida, Amphipoda, Cumacea, Halacaroidea, Priapulida, Echinodermata, Tanaidacea, and Rotifera. Total abundance and biomass of the meiobenthos showed great spatial and temporal variation, with mean values of 993.57 ± 455.36 ind cm −2 and 690.51 ± 210.64 μg 10 cm −2 , respectively. Nematodes constituted 95.60 % of the total abundance and thus had the greatest effect on meiofauna quantity and distribution, followed by copepods (1.55 %) and polychaetes (1.39 %). Meiobenthos abundance was significantly negatively correlated with water depth at stations ( r =−0.747, P

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

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

Description: A 20-yr reanalysis Experiment in the Baltic Sea Using three Dimensional Variational (3DVAR) method Ocean Science Discussions, 9, 1933-1971, 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 using three dimensional variational (3DVAR) data assimilation combining an operational numerical model with available historical temperature ( T ) and salinity ( S ) profiles. To determine the accuracy of the reanalysis, the authors present a series of comparisons with independent observations on a monthly mean basis. The performance of the assimilation in deep/shallow waters is investigated. With assimilation, temperature and salinity in the reanalysis fit better than the free run with independent measurements at different depths. Overall, the mean biases of temperature and salinity are reduced by 0.32 °C and 0.34 psu, respectively. Similarly, the mean root mean square error (RMSE) of the reanalysis is decreased by 0.35 °C and 0.3 psu compared to the free run. In space, the model error is inhomogeneous and strongly steered by the model error dynamics. Seasonally varying error of the modeled sea surface temperature is mainly controlled by the weather forcing, and shows the least improvements due to sparse 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 in the reanalysis as the model salinity in the bottom layer is increased by 2–3 psu. Sea level is also improved due to an improved density field. The correlation between model and observation is increased by 2 %–5 %, and the RMSE is generally reduced by 10 cm in the reanalysis compared to the free run. The reduction of RMSE is mainly due to the reduction of mean bias. Assimilation of T/S contributes 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. In addition, assimilation of T/S profiles results in changes of about 20 m for the mixed layer depth in the Baltic Proper region in winter. Comparisons of mixed layer depth show that the assimilation induces more changes in deep water of winter time whereas the mixed layer depth is changed only about 2 m in summer time and shallow waters. One reason could be that the effect of the assimilation is counterbalanced by the effect of heating in summer and the dominant role of the surface forcing in shallow water. The significant impact in deep waters suggests that the T/S assimilation mainly adjusts the baroclinic transport by redistributing the density field.

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 Discussions, 9, 2001-2038, 2012 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 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 Platforms (ITP) 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 and Atlantic Oceans into the deep 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 diversion of the Siberian rivers 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 waters were also subjected to strong spatial and temporal variability between 2008 and 2010. In the Canada basin, both Summer and Winter Pacific waters influence increased 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 PaW influence decreased indicative of the fact that they were not responsible for the freshening already mentioned, due to other sources. In addition, in the vicinity of the Chukchi Sea, both Summer and Winter Pacific waters were significantly warmer in 2010 than in 2008 as a consequence of a general warming trend of the Pacific waters entering in the deep Arctic Ocean since 2008. Finally, the warm Atlantic water remained relatively stable between 2008 and 2010 in the Canadian basin despite strong atmospheric shift, probably because of large time lag response. Atlantic water variability resulting from the presence of a warm "pulse-like" event in this region since 2005 was still noticeable even if a cooling effect was observed at a rate of 0.015 °C yr −1 between 2008 and 2010 in that region.

Description: Modeling long-term changes of the Black Sea ecosystem characteristics Ocean Science Discussions, 9, 2039-2080, 2012 Author(s): V. L. Dorofeyev, T. Oguz, L. I. Sukhikh, V. V. Knysh, A. I. Kubryakov, and G. K. Korotaev A three dimensional coupled physical-biological model is provided for the Black Sea to investigate its long-term changes under the synergistic impacts of eutrophication, climatic changes and population outbreak of the gelatinous invader Mnemiopsis leidyi . The model circulation field is simulated using the high frequency ERA40 atmospheric forcing as well as assimilation of the available hydrographic and altimeter sea level anomaly data for the 30 yr period of 1971–2001. The circulation dynamics are shown to resolve well the different temporal and spatial scales from mesoscale to sub-basin scale and from seasonal peaks to decadal scale trend-like changes. The biogeochemical model includes the main vertical biological and chemical interactions and processes up to the anoxic interface zone. Its food web structure is represented by two phytoplankton and zooplankton size groups, bacterioplankton, gelatinous carnivores Mnemiopsis and Aurelia , opportunistic species Noctiluca scientillans . The nitrogen cycling is accommodated by the particulate and dissolved organic nitrogen compartments and the dissolved inorganic nitrogen in the forms of ammonium, nitrite and nitrate. The ecosystem model is able to simulate successfully main observed features and trends of the intense eutrophication phase (from the early 1970s to the early 1990s), but points to its modification to simulate better the ecosystem conditions of the post-eutrophication phase.

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

Description: Wave climatology in the Arkona Basin, the Baltic Sea Ocean Science Discussions, 8, 2237-2270, 2011 Author(s): T. Soomere, R. Weisse, and A. Behrens The basic features of the wave climate in the South-Eastern Baltic Sea are studied based on available long-term measurements and simulations. The analysis of average, typical and extreme wave conditions, frequency of occurrence of different wave parameters, variations in wave heights from weekly to decadal scales, etc., is performed based on waverider measurements at the Darss Sill since 1991. The measured climatology is compared against numerical simulations with the WAM wave model driven by downscaled reanalysis of wind fields for 1958–2002 and by adjusted geostrophic winds for 1970–2007. The wave climate in this region is typical for semi-enclosed basins of the Baltic Sea. The maximum wave heights are about half of those in the Baltic Proper. The overall reliably recorded maximum significant wave height H S =4.46 m occurred during a severe S-SW storm in 1993 when the 10-min average wind speed reached 28 m s −1 . The long-term average significant wave height (0.75 m) shows modest interannual (about 12 % of the long-term mean) and substantial seasonal variation. The wave periods are mostly concentrated in a narrow range of 2.5–4 s and their distribution is almost constant over decades. The role of remote swell is very small. The annual wave properties show large interannual variability but no long-term trends in average and extreme wave heights can be observed.

Description: Improvements to the PhytoDOAS method for identification of major phytoplankton groups using hyper-spectral satellite data Ocean Science Discussions, 8, 2271-2311, 2011 Author(s): A. Sadeghi, T. Dinter, M. Vountas, B. Taylor, I. Peeken, and A. Bracher The goal of this study was to improve PhytoDOAS, which is a new retrieval method for quantitative identification of major Phytoplankton Functional Types (PFTs) using hyper-spectral satellite data. PhytoDOAS is an extension of the Differential Optical Absorption Spectroscopy (DOAS, a method for detection of atmospheric trace gases), developed for remote identification of oceanic phytoplankton groups. Thus far, PhytoDOAS has been successfully exploited to identify cyanobacteria and diatoms over the global ocean from SCIAMACHY (SCanning Imaging Absorption spectroMeter for Atmospheric CartograpHY) hyper-spectral data. The main challenge for retrieving more PFTs by PhytoDOAS is to overcome the correlation effects between different PFTs' absorption spectra. Different PFTs are composed of different types and amounts of pigments, but also have pigments in common, e.g., chl- a , causing correlation effects in the usual performance of the PhytoDOAS retrieval. Two ideas have been implemented to improve PhytoDOAS for the PFT retrieval of more phytoplankton groups. Firstly, using the fourth-derivative spectroscopy, the peak positions of the main pigment components in each absorption spectrum have been derived. After comparing the corresponding results of major PFTs, the optimized fit-window for the PhytoDOAS retrieval of each PFT was determined. Secondly, based on the results from derivative spectroscopy, simultaneous fit of PhytoDOAS has been proposed and tested for a selected set of PFTs (coccolithophores, diatoms and dinoflagllates) within an optimized fit-window. The method was then applied to the processing of SCIAMACHY data over the year 2005. Comparisons of the PhytoDOAS PFT retrievals in 2005 with the modeled PFT data from the NASA Ocean Biochemical Model (NOBM) showed similar patterns in their seasonal distributions for diatoms and coccolithophores, especially in the northern parts of the global ocean. The seasonal patterns of the PhytoDOAS coccolithophores indicated very good agreement with the global distributions of Particulate Inorganic Carbon (PIC) provided by MODIS (MODerate resolution Imaging Spectroradiometer)-Aqua level-3 products. Since PIC is known as a proxy for the abundance of coccolithophores (in open ocean), the latter agreement indicates the basic functionality of the method in retrieving coccolithophores. Moreover, as a case study, the simultaneous mode of PhytoDOAS has been applied to SCIAMACHY data for detecting a coccolithophore bloom around New Zealand (reported by NASA from MODIS imagery in December 2009); the result was quite consistent with the MODIS RGB image and the MODIS PIC map of the bloom, indicating the functionality of the method in short-term retrievals.

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

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

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

Description: Towards a regional ocean forecasting system for the IBI (Iberia-Biscay-Ireland area): developments and improvements within the ECOOP project framework Ocean Science Discussions, 8, 1937-1977, 2011 Author(s): S. Cailleau, J. Chanut, J.-M. Lellouche, B. Levier, C. Maraldi, G. Reffray, and M. Garcia Sotillo ECOOP project allowed the improvement of the regional and coastal operational forecasting systems for the different European Seas. In the Iberia-Biscay-Ireland area (IBI) a regional system has been developed and improved for the project in order to provide IBI partners the best initial and boundary conditions to their embedded coastal systems. End users could also get access to the regional hindcasts and forecasts through the ECOOP website. This system has been upgraded as follow: a first existing version V0, a second base-line version V1 ready for the ECOOP Target Operational Period and a third one V2 which consists in a new generation regional system. This paper especially pays attention to the improvements from the V1 system, whose physics are close to a large scale basin system, to the V2 one which physics are more adapted to shelf and coastal issues. Strong developments which allow further regional physics resolution in the NEMO OGCM such as tide, non linear free surface and adapted vertical mixing scheme among other have been carried out from V1 to V2 versions for the project. Thus, regional thermal fronts due to tidal mixing appear in V2 solution and are quite well placed. Moreover, simulation of the stratification in shelf areas is also improved in V2.

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

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

Description: Influence of Ross Sea Bottom Water changes on the warming and freshening of the Antarctic Bottom Water in the Australian-Antarctic Basin Ocean Science Discussions, 8, 2197-2235, 2011 Author(s): K. Shimada, S. Aoki, K. I. Ohshima, and S. R. Rintoul The WOCE Hydrographic Program (WHP) and repeated hydrographic data were used to document overall property changes of the Antarctic Bottom Water (AABW) in the Australian-Antarctic Basin between the 1990s and 2000s. Strong cooling and freshening is observed on isopycnals for layers denser than γ n = 28.30. Changes in average salinity and potential temperature below this isopycnal correspond to basin-wide warming of 1300 ± 200 TW and freshening of 24 ± 3 Gt yr −1 . While freshening can be explained by freshening of major source waters, i.e., the High Salinity Shelf Water (HSSW) of the Ross Sea and the dense shelf water formed in the Adélie and George V Land (AGVL) region, extensive warming of the AABW cannot be explained by warming of the source waters. A possible cause of warming of the AABW is a decrease in supply of the Ross Sea Bottom Water (RSBW). Hydrographic profiles between the Drygalski Trough of the Western Ross Sea and 150° E were analyzed in the context of a simple advective-diffusive model to assess the causes of the observed changes. The RSBW has also warmed by a larger amount than its source water (the HSSW). The model suggests that the warming of the RSBW observed between the 1970s and 2000s can be explained by a 21 ± 23% reduction in transport of the RSBW and an enhancement of the vertical diffusion of heat as a result of a 30 ± 7% weakening of the abyssal stratification. Freshening of the HSSW reduced the salinity and density stratification between the bottom water layer and overlying ambient water. Hence, freshening of the HSSW both directly freshened and indirectly warmed the RSBW by enhancing the vertical mixing. A simple box model suggest that changes in property and volume transport (decrease of 6.7% is assumed between the year 1995 and 2005) of the RSBW can explain 51 ± 6% of the warming and 84 ± 10% of the freshening observed in the AABW. These facts demonstrate that changes in both property and volume transport of the RSBW have contributed to the warming and freshening of the AABW in the Australian-Antarctic Basin.

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

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

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

Description: Influence of climate parameters on long-term variations of the distribution of phytoplankton biomass and nutrient concentration in the Baltic Sea simulated by a 3-D model Ocean Science Discussions, 8, 533-564, 2011 Author(s): L. Dzierzbicka-Głowacka, J. Piskozub, J. Jakacki, M. Janecki, and A. Nowicki Influence of long term trends in the climate parameters (average temperature, wind speed and solar irradiance) on phytoplankton and nutrient in the Baltic Sea was studied with an integrated three-dimensional coupled sea-ice ecological model. Simple ecosystem has been added to the sea-ice model and it has been used to estimate variability of the phytoplankton and nutrient during long term changes of the main atmospheric forces. Several numerical experiments were conducted to test the sensitivity of the model to changes of the main physical parameters such as temperature, wind speed, solar and thermal radiation (in different configurations). Influences of variability of these parameters on phytoplankton and nutrient (total inorganic nitrogen) is presented and discussion on relevance to expected future climate change is provided.

Description: A Study of the hydrographic conditions in the Adriatic Sea from numerical modelling and direct observations (2000–2008) Ocean Science Discussions, 8, 565-611, 2011 Author(s): P. Oddo and A. Guarnieri The inter-annual variability of Adriatic Sea hydrographic characteristics is investigated by means of numerical simulation and direct observation. The period investigated runs from the beginning of 2000 to the end of 2008. The model used to carry out the simulation is derived from the primitive equation component of the Adriatic Forecasting System (AFS). The model is based on the Princeton Ocean Model (POM) adapted in order to reproduce the features of the Adriatic. Both numerical findings and observations agree in depicting a strong inter-annual variability in the entire Adriatic Sea and its sub-basins. The dense water formation process has been found to be intermittent. In addition to inter-annual variability, a long-scale signal has been observed in the salinity content of the basin as a consequence of a prolonged period of reduced Po river runoff and high evaporation rates. As a result, the temperature and salinity of the northern Adriatic dense water vary considerably between the beginning and the end of the period investigated.

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

Description: A study on distribution of chlorophyll- a in the coastal waters of Anzali Port, south Caspian Sea Ocean Science Discussions, 8, 435-451, 2011 Author(s): S. Jamshidi and N. Bin Abu Bakar Phytoplankton as chlorophyll-containing organisms is the first step of production in most marine processes and food chains. Nutrient enhancement in the seawater due to the discharge of agricultural, industrial, and urban wastes threatens the Caspian Sea environment. Increasing concentrations of chlorophyll- a in seawater, in reaction to the elevation of nutrient supply can have severely damaging effects on the marine environment of the Caspian. In this research, seasonal variability of the chlorophyll- a concentrations in the western part of the southern coastal waters of the Caspian Sea near Iranian coast was examined using field observations. The data showed that the most chlorophyll- a was found below the sea surface. The thermal stratification in water column and outflow of the Anzali Lagoon affect the chlorophyll- a concentrations in the region. Concentrations of chlorophyll- a were recorded in midsummer in a range of 0.2–3.4 mg m −3 .

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

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

Description: Empirical correction of XBT fall rate and its impact on heat content analysis Ocean Science Discussions, 8, 291-320, 2011 Author(s): M. Hamon, P. Y. Le Traon, and G. Reverdin We used a collocation method between XBT and CTD/OSD (Ocean Station Data including bottle cast and low resolution CTD) from WOD05 (1°×2°×15 days) to statistically correct the XBT fall rate. An analysis of the annual median bias on depth showed that it is necessary to apply a thermal correction linked to probe calibration error, a second order correction on the depth as well as a depth offset representing measurement errors during XBT deployment. We had to separate data in several categories: shallow and deep XBT and deployment sea temperatures (below or above 10 °C). We also processed separately XBT measurements close to Japan between 1968 and 1985 due to large regional biases. Once the corrections have been applied, the analysis of heat content signal is derived from corrected XBT. From this analysis, we confirm that the maximum heat content in the top 700 m found during the 70's in early papers can be explained by the XBT biases. In addition, a trend of 0.32.10 22 J/year is observed between the period 1970 and 2008.

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

Description: Seasonal and inter-annual temperature variability in the bottom waters over the Black Sea shelf Ocean Science Discussions, 8, 321-352, 2011 Author(s): G. I. Shapiro, F. Wobus, and D. L. Aleynik Long-term changes in the state of the Bottom Shelf Water (BSW) on the Western shelf of the Black Sea are assessed using analysis of intra- and inter-annual variations of temperature as well as their relations to physical parameters of both shelf and deep-sea waters. First, large data sets of in-situ observations over the 20th century are compiled into high-resolution monthly climatology at different depth levels. Then, the temperature anomalies from the climatic mean are calculated and aggregated into spatial compartments and seasonal bins to reveal temporal evolution of the BSW. For the purpose of this study the BSW is defined as such shelf water body between the seabed and the upper mixed layer (bounded by the σ θ = 14.2 isopycnal) which has limited ability to mix vertically with oxygen-rich surface waters during the warm season (May–November) due to the formation of a seasonal pycnocline. The effects of atmospheric processes at the surface on the BSW are hence suppressed as well as the action of the "biological pump". The vertical extent of the near- bottom waters is determined based on energy considerations and the structure of the seasonal pycnocline, whilst the horizontal extent is controlled by the shelf break, where strong along-slope currents hinder exchanges with the deep sea. The BSW is shown to occupy nearly half of the area of the shelf during the summer stratification period. The potential of the BSW to ventilate horizontally during the warm season with the deep-sea waters is assessed using isopycnic analysis of temperature variations. A long-term time series of temperature anomalies in the BSW is constructed from observations during the May–November period for the 2nd half of the 20th century. The results reveal a warm phase in the 1960s/70s, followed by cooling of the BSW during 1980–2001. The transition between the warm and cold periods coincides with a regime shift in the Black Sea ecosystem. While it was confirmed that the memory of winter convection events is well preserved over the following months in the deep sea, the signal of winter cooling in the Bottom Shelf Waters significantly reduces during the warm season. The time series of temperature in the BSW is highly correlated with the temperature of Cold Intermediate Waters in the deep sea thus indicating that the isopycnal exchanges with the deep sea are more important for inter-annual/inter-decadal variability of the BSW on the Western Black Sea shelf than winter convection on the shelf itself.

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

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

Description: Technical Note: Detection of gas bubble leakage via correlation of water column multibeam images Ocean Science Discussions, 8, 1757-1775, 2011 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. However, up to the present, the extremely high data rate hampers water column backscatter investigations. More sophisticated visualization and processing techniques for water column backscatter analysis are still under development. We here present such water column backscattering data gathered with a 50 kHz prototype multibeam system. Water column backscattering data is presented in videoframes grabbed over 75 s and a "re-sorted" singlebeam presentation. Thus individual gas bubbles rising from the 24 m deep seafloor clearly emerge in the acoustic images and rise velocities can be determined. A sophisticated processing scheme is introduced to identify those rising gas bubbles in the hydroacoustic data. It applies a cross-correlation technique similar to that used in Particle Imaging Velocimetry (PIV) to the acoustic backscatter images. Tempo-spatial drift patterns of the bubbles are assessed and match very well measured and theoretical rise patterns. The application of this processing scheme to our field data gives impressive results with respect to unambiguous bubble detection and remote bubble rise velocimetry. The method can identify and exclude the main driver for misinterpretations, i.e. fish-mediated echoes. Even though image-based cross-correlation techniques are well known in the field of fluid mechanics for high resolution and non-inversive current flow field analysis, this technique was never applied in the proposed sense for an acoustic bubble detector.

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

Description: Sea surface temperature anomalies, seasonal cycle and trend regimes in the eastern Pacific coast Ocean Science Discussions, 8, 1215-1232, 2011 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 radiation (TSI) and then, the regime shift detector (RSD) was applied to detect possible temperature regimes in the series. Posterior to this, 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 in the late 1990's when a shift to negative values is detected. Such a shift lasted until the early 2000's when positive anomalies appear again but there is 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 in sea surface temperature anomalies and with the regime changes. This would be due to a difference in the amount of energy received from the sun. Comparing two consecutive periods, 1952–1975 with 1977–1999, the second received 0.39 % more energy (approximately 3 × 10 8 J m −2 ) from the sun. Seasonal cycles show larger range at northern latitudes (〉40° N), northern tropical-template transition zone (20°–26° N) and in the tropical-equatorial band (0°–30° S). The smaller ranges occur at 0°–16° N and 50°–60° S. The residuals (seasonal minus the theoretical curve) indicate a clear modulation due to advection by ocean currents.

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

Description: Modelling the variability of the Antarctic Slope Current Ocean Science Discussions, 8, 1-38, 2011 Author(s): P. Mathiot, H. Goosse, T. Fichefet, B. Barnier, and H. Gallée One of the main features of the oceanic circulation along Antarctica is the Antarctic Slope Current (ASC). This circumpolar current flows westward and allows communication between the three major basins around Antarctica. The ASC is not very well known due to difficult access and the presence of sea ice during several months, allowing in situ study only during summertime. Moreover, only few numerical studies of this current have been carried out. Here, we investigate the sensitivity of this current to two different atmospheric forcing sets and to four different resolutions in a coupled ocean-sea ice model (NEMO-LIM). Two sets of simulation are conducted. For the first set, global model configurations are run at coarse (2°) to eddy permitting resolutions (0.25°) with the same atmospheric forcing. For the second set, simulations with two different atmospheric forcing sets are performed with a regional circumpolar configuration (south of 30° S) at 0.5° resolution. The first atmospheric forcing set is based on ERA40 reanalysis and CORE data, while the second one is based on a downscaling of the reanalysis ERA40 by the MAR regional atmospheric model. Sensitivity experiments to resolution show that a minimum model resolution of 0.5° is needed to capture the dynamics of the ASC in term of transport and recirculation. Sensitivity of the ASC to atmospheric forcing fields shows that the wind speed along the Antarctic coast strongly controls the transport and the seasonal cycle of the ASC. An increase of the Easterlies by about 30% leads to an increase of the mean transport of ASC by about 40%. Similar effects are obtained on the seasonal cycle: using a forcing fields with a stronger amplitude of the seasonal cycle leads to double the amplitude of the seasonal cycle of the ASC. To confirm the importance of the wind speed, a simulation, where the seasonal cycle of the wind speed is removed, is carried out. This simulation shows a decrease by more than 50% of the amplitude of the seasonal cycle without changing the mean value of ASC transport.

Description: New constraints on the Eastern Mediterranean δ 18 O: δ D relationship Ocean Science Discussions, 8, 39-53, 2011 Author(s): K. A. Cox, E. J. Rohling, G. A. Schmidt, R. Schiebel, S. Bacon, D. A. Winter, M. Bolshaw, and H. J. Spero Previous work on oxygen and hydrogen isotope data from Eastern Mediterranean water samples has defined a mixing relationship in this region that is different from the world surface ocean. This prompted speculations about the hydrological processes in the Mediterranean region. We present new δ 18 O and δ D data from the Eastern Mediterranean region and the East Greenland Current system, spanning a wide salinity range. These data define δ 18 O: δ D relationships for both regions that are consistent with the world surface ocean δ 18 O: δ D relationship, despite the highly evaporative conditions that prevail in the Mediterranean region. These new geochemical data have suggested that the world surface ocean &delta 18 O: δ D relationship holds throughout almost the entire global salinity range.

Description: Tidal generation of large sub-mesoscale eddy dipoles Ocean Science Discussions, 8, 723-760, 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 Gwai (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 to the interior of the subpolar gyre.

Description: ENSURF: multi-model sea level forecast – implementation and validation results for the IBIROOS and Western Mediterranean regions Ocean Science Discussions, 8, 761-800, 2011 Author(s): B. Pérez, R. Brower, J. Beckers, D. Paradis, C. Balseiro, K. Lyons, M. Cure, M. G. Sotillo, B. Hacket, M. Verlaan, and E. Alvarez Fanjul ENSURF (Ensemble SURge Forecast) is a multi-model application for sea level forecast that makes use of existing storm surge or circulation models today operational in Europe, as well as near-real time tide gauge data in the region, with the following main goals: – providing an easy access to existing forecasts, as well as to its performance and model validation, by means of an adequate visualization tool – generation of better forecasts of sea level, including confidence intervals, by means of the Bayesian Model Average Technique (BMA) The system was developed and implemented within ECOOP (C.No. 036355) European Project for the NOOS and the IBIROOS regions, based on MATROOS visualization tool developed by Deltares. Both systems are today operational at Deltares and Puertos del Estado respectively. The Bayesian Modelling Average technique generates an overall forecast probability density function (PDF) by making a weighted average of the individual forecasts PDF's; the weights represent the probability that a model will give the correct forecast PDF and are determined and updated operationally based on the performance of the models during a recent training period. This implies the technique needs the availability of sea level data from tide gauges in near-real time. Results of validation of the different models and BMA implementation for the main harbours will be presented for the IBIROOS and Western Mediterranean regions, where this kind of activity is performed for the first time. The work has proved to be useful to detect problems in some of the circulation models not previously well calibrated with sea level data, to identify the differences on baroclinic and barotropic models for sea level applications and to confirm the general improvement of the BMA forecasts.

Description: Joint use of satellite and in-situ data for coastal monitoring Ocean Science Discussions, 8, 955-998, 2011 Author(s): F. Gohin Sea surface Temperature, Chlorophyll and turbidity are three variables of the coastal environment commonly measured by monitoring networks. The observation networks are often based on coastal stations which do not provide a sufficient coverage to val-idate the model outputs or to be used in assimilation over the continental shelf. Conversely, the products derived from satellite reflectance show generally a decreasing quality shoreward and an accurate assessment of these data is required. In this text, we show that the satellite-derived chlorophyll products, obtained through a dedicated coastal algorithm, fulfil the first requirement of a monitoring system: the ability to represent correctly the mean annual cycle. The annual cycle, mean and percentile 90 of the chlorophyll concentration, derived from MERIS/ESA and MODIS/NASA, have been compared to in-situ observations at twenty six selected stations from the Mediterranean Sea to the North-Sea. Keeping in mind the validation, the forcing or the assimilation in hydrological, sediment-transport or ecological models, the non-algal Suspended Particulate Matter (SPM) is also a parameter which is expected from the satellite imagery. However, the monitoring networks measure essentially the turbidity and a consistency between chlorophyll, representative of the phytoplankton biomass, non-algal SPM, and turbidity is required. In this study, we derive the satellite turbidity from chlorophyll and non-algal SPM with a common formula applied to in-situ or satellite observations. The distribution of the satellite-derived turbidity shows the same main statistical characteristics that measured in-situ; which satisfies our first condition to monitor the long-term changes or the large-scale spatial variation over the continental shelf and along the shore. For the first time, maps of turbidity, so useful for the surveillance of the benthic habitats, are proposed operationally from space on areas as different as the Southern North-Sea or the Western Mediterranean Sea, with validation at coastal stations.