Altimeter assimilation in the OCCAM global model: Part I: A twin experiment

An altimeter data assimilation scheme has been tested in the OCCAM (Ocean Circulation and Climate Advanced Modelling) global 1/4°, 36-level model using a twin experiment format. The Cooper and Haines displacement scheme is used. The method works well in most regions and depths. Currents and densities in the top 1000 m generally improve by over 50–70% after 5 months of sea level assimilation every 15 days. Below 1000 m, an error reduction of up to 50% is achieved. The errors remain low during a further 60-day run without assimilation. Diagnostics for the North Atlantic, the Tropical Pacific and the Antarctic Circumpolar Current are shown alongside the global averages.The main problems encountered were in weakly stratified regions of the Antarctic and Arctic seas. A scale selective filter is developed to avoid assimilating scales much larger than the local deformation radius, and this avoids the adverse assimilation effects in the southern oceans. A companion paper uses this scheme to assimilate TOPEX and ERS-1 altimeter maps.     

Validation of a hybrid optimal interpolation and Kalman filter scheme for sea surface temperature assimilation

A hybrid data assimilation scheme designed for operational assimilation of satellite sea surface temperatures (SST) into an ocean model has been developed and validated against in-situ observations. The scheme consists of an optimal interpolation (OI) part and a greatly simplified Kalman filter (KF) part.The OI is performed only in the longitudinal and latitudinal directions. A climatological field is used as a background field for the interpolation. It is constructed by fitting daily averages of satellite SST to the annual mean, annual, and semiannual harmonics in a 20 km by 20 km grid. The background error covariance is approximated by a spatially varying two-dimensional exponential covariance model. The parameters of the covariance model are fitted to the deviations of the satellite data from the background field using data from a full year.The simplified KF uses ocean model forecasts as a background field. It is based on the assumption that it is possible to neglect horizontal SST covariances in the filter and that the typical time scale for vertical mixing in the mixed layer is much shorter than the average time between observations. We therefore assume that the error variance in a column of water is evenly spread out throughout the mixed layer. The result of these simplifications is a computationally very efficient KF.A one year validation of the scheme is performed for year 2001 using an operational eddy resolving ocean model covering the North Sea and the Baltic Sea. It is found that assimilation of sea surface temperature data reduces the model root mean square error from 1.13 °C to 0.70 °C. The hybrid scheme is found to reduce the root mean square error slightly more than the simplified KF without OI to 0.66 °C. The inclusion of spatially varying satellite error variances does not improve the performance of the scheme significantly.     

1/32° real-time global ocean prediction and value-added over 1/16° resolution

A 1/32° global ocean nowcast/forecast system has been developed by the Naval Research Laboratory at the Stennis Space Center. It started running at the Naval Oceanographic Office in near real-time on 1 Nov. 2003 and has been running daily in real-time since 1 Mar. 2005. It became an operational system on 6 March 2006, replacing the existing 1/16° system which ceased operation on 12 March 2006. Both systems use the NRL Layered Ocean Model (NLOM) with assimilation of sea surface height from satellite altimeters and sea surface temperature from multi-channel satellite infrared radiometers. Real-time and archived results are available online at http://www.ocean.nrlssc.navy.mil/global_nlom. The 1/32° system has improvements over the earlier system that can be grouped into two categories: (1) better resolution and representation of dynamical processes and (2) design modifications. The design modifications are the result of accrued knowledge since the development of the earlier 1/16° system. The improved horizontal resolution of the 1/32° system has significant dynamical benefits which increase the ability of the model to accurately nowcast and skillfully forecast. At the finer resolution, current pathways and their transports become more accurate, the sea surface height (SSH) variability increases and becomes more realistic and even the global ocean circulation experiences some changes (including inter-basin exchange). These improvements make the 1/32° system a better dynamical interpolator of assimilated satellite altimeter track data, using a one-day model forecast as the first guess. The result is quantitatively more accurate nowcasts, as is illustrated by several model-data comparisons. Based on comparisons with ocean color imagery in the northwestern Arabian Sea and the Gulf of Oman, the 1/32° system has even demonstrated the ability to map small eddies, 25–75 km in diameter, with 70% reliability and a median eddy center location error of 22.5 km, a surprising and unanticipated result from assimilation of altimeter track data. For all of the eddies (50% small eddies), the reliability was 80% and the median eddy center location error was 29 km. The 1/32° system also exhibits improved forecast skill in relation to the 1/16° system. This is due to (a) a more accurate initial condition for the forecast and (b) better resolution and representation of critical dynamical processes (such as upper ocean – topographic coupling via mesoscale flow instabilities) which allow the model to more accurately evolve these features in time while running in forecast mode (forecast atmospheric forcing for the first 5 days, then gradually reverting toward climatology for the remainder of the 30-day forecast period). At 1/32° resolution, forecast SSH generally compares better with unassimilated observations and the anomaly correlation of the forecast SSH exceeds that from persistence by a larger amount than found in the 1/16° system.     

Influence of the heavy fuel spill from the Prestige tanker wreckage in the overlying seawater column levels of copper, nickel and vanadium (NE Atlantic ocean)

The water column above the Prestige wreckage was sampled during two consecutive campaigns: Prestinaut (December 2002) two weeks after the tanker sunk and HidroPrestige0303 (March 2003) one month after the sealing of the main fuel leaks. Samples of the original cargo fuel and the emulsified fuel in the surface of the ocean were also collected. Analysis of the fuel indicated the release of 135 kg of Cu, 1700 kg of Ni and 5300 kg of V from the original fuel to the water column, remaining 35 kg of Cu, 3100 kg of Ni and 13,800 kg of V in the emulsified fuel. The metal partitioning between the water column and the emulsioned floating fuel, Cu > Ni ~ V, are in accordance with the stability index for the metal–nitrogen bond in metalloporphyrins. This release had an impact on dissolved trace metal concentrations in the water column. An increase on dissolved copper (2.8–4.7 nM) and nickel (2.2–8.0 nM) with respect to natural values (1–3 nM for Cu and 1.6–5 nM for Ni) was observed. Values for vanadium (28–35 nM) were in the range of pristine North Atlantic waters (30–36 nM). This contamination was especially observed in the upper water column (0–50 m), associated with the mixing of seawater with the fuel moving upwards, and in deep waters, where the residence time of fuel is higher. Future research in this field should focus on the environmental variables and the processes that control the release of contaminants from fuels for a better assessment of the contamination in oil-spill events.     

Parameterization of air–sea gas fluxes at extreme wind speeds

Air–sea flux measurements of O₂ and N₂ obtained during Hurricane Frances in September 2004 [D'Asaro, E. A. and McNeil, C. L., 2006. Measurements of air–sea gas exchange at extreme wind speeds. Journal Marine Systems, this edition.] using air-deployed neutrally buoyant floats reveal the first evidence of a new regime of air–sea gas transfer occurring at wind speeds in excess of 35 m s^− 1. In this regime, plumes of bubbles 1 mm and smaller in size are transported down from near the surface of the ocean to greater depths by vertical turbulent currents with speeds up to 20−30 cm s^− 1. These bubble plumes mostly dissolve before reaching a depth of approximately 20 m as a result of hydrostatic compression. Injection of air into the ocean by this mechanism results in the invasion of gases in proportion to their tropospheric molar gas ratios, and further supersaturation of less soluble gases. A new formulation for air–sea fluxes of weakly soluble gases as a function of wind speed is proposed to extend existing formulations [Woolf, D.K, 1997. Bubbles and their role in gas exchange. In: Liss, P.S., and Duce, R.A., (Eds.), The Sea Surface and Global Change. Cambridge University Press, Cambridge, UK, pp. 173–205.] to span the entire natural range of wind speeds over the open ocean, which includes hurricanes. The new formulation has separate contributions to air–sea gas flux from: 1) non-supersaturating near-surface equilibration processes, which include direct transfer associated with the air–sea interface and ventilation associated with surface wave breaking; 2) partial dissolution of bubbles smaller than 1 mm that mix into the ocean via turbulence; and 3) complete dissolution of bubbles of up to 1 mm in size via subduction of bubble plumes. The model can be simplified by combining “surface equilibration” terms that allow exchange of gases into and out of the ocean, and “gas injection” terms that only allow gas to enter the ocean. The model was tested against the Hurricane Frances data set. Although all the model parameters cannot be determined uniquely, some features are clear. The fluxes due to the surface equilibration terms, estimated both from data and from model inversions, increase rapidly at high wind speed but are still far below those predicted using the cubic parameterization of Wanninkhof and McGillis [Wannikhof, R. and McGillis, W.R., 1999. A cubic relationship between air–sea CO₂ exchange and wind speed. Geophysical Research Letters, 26:1889–1892.] at high wind speed. The fluxes due to gas injection terms increase with wind speed even more rapidly, causing bubble injection to dominate at the highest wind speeds.     

Microscale variability in the distributions of large fluorescent particles observed in situ with a planar laser imaging fluorometer

It has been known for decades that particle-size and biomass spectra show regular patterns in the ocean, and that these patterns often show systematic variations with other properties such as total biomass, nutrient concentration, season, and distance (both vertical and horizontal). The recent finding of the ubiquitous nature of layers of phytoplankton < 1 m thick prompted us to explore the fine- and microscale vertical variations of size- and fluorescence-abundance spectra in the ocean. Using a two-dimensional planar laser imaging system mounted on a free-falling platform, we quantified the properties of large fluorescent particles ( 20 μm–2 cm) through the water column, obtaining images every 10–30 cm. These images showed systematic relationships of the spectral properties to total chlorophyll: increased proportions of the smallest particles at high chlorophyll concentrations, and a lengthening of the spectral size range at high total chlorophyll concentrations (more large particles at high chlorophyll concentrations). Further, we observed significant variations of the spectral properties over scales of 1 m and less, and recorded the frequent occurrence of unusual layers of large particles. Our new instrument, which is sensitive to thin layers of enhanced phytoplankton biomass, shows the planktonic community to be highly structured vertically on scales of 1–2 m, particularly within the DCM.     

Comparison between a reanalyzed product by 3-dimensional variational assimilation technique and observations in the Ulleung Basin of the East/Japan Sea

Reanalyzed products from a MOM3-based East Sea Regional Ocean Model with a 3-dimentional variational data assimilation module (DA-ESROM), have been compared with the observed hydrographic and current datasets in the Ulleung Basin (UB) of the East/Japan Sea (EJS). Satellite-borne sea surface temperature and sea surface height data, and in-situ temperature profiles have been assimilated into the DA-ESROM. The performance of the DA-ESROM appears to be efficient enough to be used in an operational ocean forecast system.Comparing with the results from Mitchell et al. [Mitchell, D. A., Watts, D. R., Wimbush, M., Teague, W.J., Tracey, K. L., Book, J. W., Chang, K.-I., Suk, M.-S., Yoon, J.-H., 2005a. Upper circulation patterns in the Ulleung Basin. Deep-Sea Res. II, 52, 1617-1638.], the DA-ESROM fairly well simulates the high variability of the Ulleung Warm Eddy and Dok Cold Eddy as well as the branching of the Tsushima Warm Current in the UB. The overall root-mean-square error between 100 m temperature field reproduced by the DA-ESROM and the observed 100-dbar temperature field is 2.1 °C, and the spatially averaged grid-to-grid correlation between the two temperature fields is high with a mean value of 0.79 for the inter-comparison period.The DA-ESROM reproduces the development of strong southward North Korean Cold Current (NKCC) in summer consistent with the observational results, which is thought to be an improvement of the previous numerical models in the EJS. The reanalyzed products show that the NKCC is about 35 km wide, and flows southward along the Korean coast from spring to summer with maximum monthly mean volume transport of about 0.8 Sv in August–September.     

A mesoscale index to describe the regional circulation around the Balearic Islands

Historical oceanographic surveys carried out around the Balearic Islands (western Mediterranean) suggest two different scenarios for the regional ocean circulation. In one scenario, occurring during cold winters, cool water is formed at intermediate layers (100–300 m) in the Gulf of Lions. This Western Mediterranean Intermediate Water (WIW) usually moves southward reaching the Balearic Channels, deflecting the warmer Levantine Intermediate Water (LIW) coming from the Eastern Mediterranean, and even blocking the Ibiza Channel. On the other hand, during mild winters, less WIW is formed and then LIW flows through the channels, appearing at their characteristic depths. The oceanographic surveys around the Balearic Islands (1985–2004) have provided a qualitative index, indicating the presence or not of WIW in the Ibiza Channel, based on the analyses of θS diagrams. A quantitative index based on mean water temperature between 100 and 300 m depth in the channels may also be defined. Both indexes, the qualitative and the quantitative, give consistent information on WIW presence for the period 1985–2004, however, both are short in time and have gaps in the series. In order to obtain a longer and continuous index for WIW presence and then for regional circulation, air–sea heat fluxes at the Gulf of Lions during winter months were obtained from the meteorological NCEP/NCAR reanalysis dataset and compared with other meteorological data such as surface air temperature. The standardized air temperature anomalies at 1000 hPa in the Gulf of Lions during winter (December–March) has been shown to be the simplest and best indicator of absence/presence of WIW in the Balearic Islands channels in late spring. Values above 1.0 of the standardized temperature anomaly would indicate absence of WIW in the Ibiza Channel. The high correlation obtained with available in-situ oceanographic data allows the use of this index as an indicator of presence of WIW and then of different regional circulation scenarios backwards in time and in those years for which the oceanographic data are missing or scarce.     

Turbulent mixing and internal tides in Gaoping (Kaoping) Submarine Canyon, Taiwan

Turbulent overturning on scales greater than 10 m is observed near the bottom and in mid-depth layers within the Gaoping (formerly spelled Kaoping) Submarine Canyon (KPSC) in southern Taiwan. Bursts of strong turbulence coexist with bursts of strong sediment concentrations in mid-depth layers. The turbulence kinetic energy dissipation rate in some turbulence bursts exceeds 10^− 4 W kg^− 1, and the eddy diffusivity exceeds 10^− 1 m² s^− 1. Within the canyon, the depth averaged turbulence kinetic energy dissipation rate is ~ 7 × 10^− 6 W kg^− 1, and the depth averaged eddy diffusivity is ~ 10^− 2 m² s^− 1. These are more than two orders of magnitude greater than typical values in the open ocean, and are much larger than those found in the Monterey Canyon where the strong turbulent mixing has also been. The interaction of tidal currents with the complex topography in Gaoping Submarine Canyon is presumably responsible for the observed turbulent overturning via shear instability and the breaking of internal tides and internal waves at critical frequencies. Strong 1st-mode internal tides exist in KPSC. The depth averaged internal tidal energy near the canyon mouth is ~ 0.17 m² s^− 2. The depth integrated internal tidal energy flux at the mouth of the canyon is ~ 14 kW m^− 1, propagating along the axis of the canyon toward the canyon head. The internal tidal energy flux in the canyon is 3–7 times greater than that found in Monterey Canyon, presumably due to the more than 10 times larger barotropic tide in the canyon. Simple energy budget calculations conclude that internal tides alone may provide energy sufficient to explain the turbulent mixing estimated within the canyon. Further experiments are needed in order to quantify the seasonal and geographical distributions of internal tides in Gaoping Submarine Canyon and their effects on the sediment flux in the canyon.     

Recovery of North-East Atlantic temperature fields from profiling floats: Determination of the optimal float number from sampling and instrumental error analysis

Argo is an international project that is deploying an array of temperature and salinity profiling floats over the global ocean. Here we use the error formulation derived from Optimal Statistical Interpolation to estimate statistical errors associated with the recovery of the temperature field in the North-East Atlantic ocean. Results indicate that with the present distribution of floats (119 in the considered domain), scales of wavelength larger than 500 km can be recovered with a relative uncertainty (rms error relative to the standard deviation of the field) of about 7% at 50 m, 8% at 200 m and 10% at 1000 m. This corresponds to mean absolute errors of 0.111 °C at 50 m, 0.104 °C at 200 m and 0.073 °C at 1000 m.The splitting of total errors into instrumental and sampling contributions reveals that, in the present scenario, errors are more due to the small number of floats than to instrumental errors, especially at upper levels. For scales larger than 500 km this will hold true until 200–250 floats are deployed (less than 200 for deep levels). In such a simulated scenario, the number of observations and the technology become approximately equally limiting factors for the accuracy of the temperature field mapping, with total relative errors of less than 2% at upper levels and about 3% at 1000 m.     

Bathymetry impacts on water exchange modelling through the Danish Straits

Deep and narrow channels in Danish Straits are one of the governing factors for the Baltic–North Sea water and salt exchange. The channels have a depth up to 50 m and a horizontal scale of a few hundred meters. The typical horizontal resolution used in current operational three dimensional Baltic–North Sea models is 1 nautical mile (nm) which can not well resolve these deep channels. In this paper, an alternative method is used to generate the 1 nm resolution bathymetry so that the deep channel is well resolved and at the same time the total water volume is roughly conserved. The impact of the new bathymetry on modelling water and salt transports as well as temperature and salinity structure is assessed by comparing a 3-year model run with the adjusted bathymetry and a control run with the averaged bathymetry. Volume and salt transports through the Great Belt are examined in the two runs. The results show that the model ocean is dominated by a typical two-layer transport (i.e., upper brackish Baltic outflow and lower saltier inflow), and the new bathymetry significantly enhances the two-layer transport. The lighter Baltic outflow is increased by 18% in the upper 10 m and saltier deep inflow is increased by 300% (in comparison with the old bathymetry) below 10 m. The total net transport into the Baltic Sea is increased by 13%. The temperature and salinity structure is also significantly influenced by the bathymetry, especially during inflow events. The stratification is strengthened and the bottom salinity is increased in Danish Straits and adjacent waters. The bathymetry impact is found significant through the entire 3-year model run period, and the signal is propagated to a large area covering the Stopple Channel. Comparison with observations show that such changes are positive improvements to the models. The results suggest that the deep channels in the Danish Straits have to be carefully resolved in order to correctly simulate the Baltic–North Sea water exchange.     

A Mean Dynamic Topography of the Mediterranean Sea computed from altimetric data, in-situ measurements and a general circulation model

In the Mediterranean Sea, where the mean circulation is largely unknown and characterized by smaller scales and less intensity than in the open ocean, the interpretation of altimetric Sea Level Anomalies (SLA) is rather difficult. In the context of operational systems such as MFS (Mediterranean Forecasting System) or MERCATOR, that assimilate the altimetric information, the estimation of a realistic Mean Dynamic Topography (MDT) consistent with altimetric SLA to be used to reconstruct absolute sea level is a crucial issue. A method is developed here to estimate the required MDT combining oceanic observations as altimetric and in-situ measurements and outputs from an ocean general circulation model (OGCM).In a first step, the average over the 1993–1999 period of dynamic topography outputs from MFS OGCM provides a first guess for the computation of the MDT. Then, in a second step, drifting buoy velocities and altimetric data are combined using a synthetic method to obtain local estimates of the mean geostrophic circulation which are then used to improve the first guess through an inverse technique and map the MDT field (hereafter the Synthetic Mean Dynamic Topography or SMDT) on a 1/8° resolution grid.Many interesting current patterns and cyclonic/anticyclonic structures are visible on the SMDT obtained. The main Mediterranean coastal currents are well marked (as the Algerian Current or the Liguro–Provenço–Catalan Current). East of the Sicily channel, the Atlantic Ionian Stream divides into several main branches crossing the Ionian Sea at various latitudes before joining at 19°E into a unique Mid-Mediterranean Jet. Also, strong signatures of the main Mediterranean eddies are obtained (as for instance the Alboran gyre, the Pelops, Ierapetra, Mersa-Matruh or Shikmona anticyclones and the Cretan, Rhodes or West Cyprius cyclones). Independent in-situ measurements from Sea Campaigns NORBAL in the North Balearic Sea and the North Tyrrhenian Sea and SYMPLEX in the Sicily channel are used to validate locally the SMDT: deduced absolute altimetric dynamic topography compares well with in-situ observations. Finally, the SMDT is used to compute absolute altimetric maps in the Alboran Sea and the Algerian Current. The use of absolute altimetric signal allows to accurately follow the formation and propagation of cyclonic and anticyclonic eddies in both areas.     

Effect of wave boundary layer on sea-to-air dimethylsulfide transfer velocity during typhoon passage

A full-spectral third-generation ocean wind–wave model (Wavewatch-III) implemented in the South China Sea is used to investigate the effects of the wave boundary layer on the drag coefficient and the sea-to-air transfer velocity of dimethylsulfide (DMS) during passage of Typhoon Wukong (September 5–11, 2000) with a maximum sustained wind speed of 38 m s^− 1. The model is driven by the reanalyzed surface winds (1° × 1°, four times daily) from the National Centers for Environmental Prediction. It is found that the wave boundary layer evidently enhances (16.5%) the drag coefficient (in turn increases the momentum flux across the air–sea interface), and reduces (13.1%) the sea-to-air DMS transfer velocity (in turn decreases the sea-to-air DMS flux). This indicates the possibility of important roles of wave boundary layer in atmospheric DMS contents and global climate system.     

Multivariate assimilation in MERCATOR project: New statistical parameters from forecast error estimation

The new operational prototype of Mercator (french Global Ocean Data Assimilation Experiment contribution) is composed of a North Atlantic primitive equation ocean model OPA (Ocean Parallel Algorithm between 20°S and 70°N, [Madec, G., P. Delecluse, M. Imbard and C. Lévy (1998). OPA8.1 ocean general circulation model reference manuel. Notes du pôle de modélisation IPSL. n°11: 91p]) and of a multivariate and multidata assimilation scheme [De Mey, P. and M. Benkiran (2002). “A multivariate reduced-order optimal interpolation method and its application in Mediterranean basin-scale circulation.” Ocean Forecasting : Conceptual basis and application, Pinardi, N., Springer Verlag.] This system has already given some significant improvements from previous Mercator configurations (M. Benkiran, personal communication). However some biases on ocean state still remain in the tropics where the reduced-order optimal interpolation scheme is suspected to be ill-parameted in the model forecast error. Indeed the guess error covariance matrix is decomposed into an error variance value and a spatio-temporal correlation function which are assumed to have some “good” properties (spatial homogeneity of the correlation function, constant ratio between signal and error variance). This study shows how we can use ensemble methods to validate these assumptions. We can see that the correlation function can reach negative values locally, mostly in regions of high variability contradictory with the homogeneous hypothesis. The reduced space used in the operational configuration is based on the signal seasonal Empirical Orthogonal Functions (EOFs). An empirical relationship between signal and error variance has been set and the correlation function is the same on every dimension of the reduced space. By projection of the estimated guess error variance onto the reduced space, we find a repartition of this quantity quite different to what was set in the system. The error statistics is found to be inhomogeneous compared to hypothesis made in the assimilation scheme. These two new parameters tested separately in the assimilation scheme gives significant improvements of the forecast and analysis results. This is particularly obvious in the tropics. But relationship between signal and error statistics (as assumed in the optimal interpolation) is found to be complex.     

Inter-comparing five forecast operational systems in the North Atlantic and Mediterranean basins: The MERSEA-strand1 methodology

The methodology to achieve a real time inter-comparison of five state-of-the-art operational forecast systems for the North Atlantic and Mediterranean basins is presented. All systems provide analysis and near real-time prediction of the three-dimensional ocean through Opendap servers. A standard set of diagnostics called metrics, is described. Definition and examples of metrics are given. An inter-comparison of the five systems is conducted over a 1 year period using those metrics. It is shown that the methodology developed allows a successful inter-comparison. It has been adopted by the GODAE community. It is also shown that the systems are consistent with the current knowledge of the ocean circulation and climatologies. Systems are deficient in the representation of specific water masses characteristics as Mode waters. Data assimilation of vertical profiles of temperature and salinity solve such deficiencies. Metrics also allow a monitoring of the system's North Atlantic overturning stream function and will allow detecting any changes in the coming year system's thermohaline circulation.     

Assimilation of SeaWiFS ocean chlorophyll data into a three-dimensional global ocean model

Sea-viewing Wide Field-of-view Sensor (SeaWiFS) chlorophyll data were assimilated with an established three-dimensional global ocean model. The assimilation improved estimates of chlorophyll relative to a free-run (no assimilation) model. Compared to SeaWiFS, annual bias of the assimilation model was 5.5%, with an uncertainty of 10.1%. The free-run model had a bias of 21.0% and an uncertainty of 65.3%. In situ data were compared to the assimilation model over a 6-year time period from 1998 through 2003, indicating a bias of 0.1%, and an uncertainty of 33.4% for daily coincident, co-located data. SeaWiFS bias was slightly higher at − 1.3% and nearly identical uncertainty at 32.7%. The free-run bias and uncertainty at − 1.4% and 61.8%, respectively, indicated how much the assimilation improved model results. Annual primary production estimates for the 1998–2003 period produced a nearly 50% improvement by the assimilation model over the free-run model as compared to a widely used algorithm using SeaWiFS chlorophyll data. These results suggest the potential of assimilation of satellite ocean chlorophyll data for improving model results.     

Estimation of the surface and mid-depth currents from Argo floats in the Pacific and error analysis

With the large deployment, the Array for Real-time Geostrophic Oceanography program has great potential for measuring the ocean currents both on the surface and at mid-depth. However the positioning error of fixes in a trajectory varies from 150 m to 1000 m, and thus created difficulty for accurate estimations of the surface and mid-depth currents. Also the reliability of the estimated surface and mid-depth currents requires accurate error estimations.In this study a new sequential method of Argo float surface trajectory tracking and extrapolating is proposed based on Kalman Filter (KF), under the presumption that a surface trajectory of Argo float is dominated by a constant current plus inertial oscillation. This trajectory tracking and extrapolating method is able to reduce the positioning uncertainties of Argo surface trajectories and provides error estimations. When this method was applied to extrapolate the positions when float resurfacing and descending, the estimation error of the mid-depth currents can be reduced. Utilizing this method in the Pacific, surface and mid-depth currents were estimated from surface trajectories of Argo floats from 2001 to 2004, along with their detailed error estimations. The average error for surface currents is about 4.4 cm s^− 1 which is equivalent to the accuracy order (5 cm s^− 1) of the Surface Velocity Program drifters. The estimation error of the mid-depth currents at 1000 db is reduced to about 0.21 cm s^− 1 without considering the effect of vertical shear.This study shows that the surface trajectory from Argo float provides a new means to measure surface circulations in the global ocean at real time, and that the estimated mid-depth current could be one of the important sources to improve the understanding for ocean dynamic.     

Silicon dynamics in the Oder estuary, Baltic Sea

Studies on dissolved silicate (DSi) and biogenic silica (BSi) dynamics were carried out in the Oder estuary, Baltic Sea in 2000–2005. The Oder estuary proved to be an important component of the Oder River–Baltic Sea continuum where very intensive seasonal DSi uptake during spring and autumn, but also BSi regeneration during summer take place. Owing to the regeneration process annual DSi patterns in the river and the estuary distinctly differed; the annual patterns of DSi in the estuary showed two maxima and two minima in contrast to one maximum- and one minimum-pattern in the Oder River. DSi concentrations in the river and in the estuary were highest in winter (200–250 μmol dm^− 3) and lowest (often less than 1 μmol dm^− 3) in spring, concomitant with diatom growth; such low values are known to be limiting for new diatom growth. Secondary DSi summer peaks at the estuary exit exceeded 100 μmol dm^− 3, and these maxima were followed by autumn minima coinciding with the autumn diatom bloom. Seasonal peaks in BSi concentrations (ca. 100 μmol dm^− 3) occurred during the spring diatom bloom in the Oder River. Mass balance calculations of DSi and BSi showed that DSi + BSi import to the estuary over a two year period was 103.2 kt and that can be compared with the DSi export of 98.5 kt. The difference between these numbers gives room for ca. 2.5 kt BSi to be annually exported to the Baltic Sea. Sediment cores studies point to BSi annual accumulation on the level of 2.5 kt BSi. BSi import to the estuary is on the level of ca. 10.5 kt, thus ca. 5 kt of BSi is annually converted into the DSi, increasing the pool of DSi that leaves the system. BSi concentrations being ca. 2 times higher at the estuary entrance than at its exit remain in a good agreement with the DSi and BSi budgeting presented in the paper.     

Hydrographic structure and zooplankton abundance and diversity off Paita, northern Peru (1994 to 2004) — ENSO effects, trends and changes

The objective of the present study was to verify possible spatial, seasonal, and inter-annual changes in the zooplankton off Paita (northern Peru), an upwelling area located closely to the limits of cold Humboldt Current and warm Equatorial Surface Waters. Zooplankton was sampled at subsurface on 53 occasions from August 1994 to December 2004 at four stations located 2 to 30 km offshore with a WP-2 net (300 µm). Extremely high surface water temperatures combined with low salinities were observed during the 1997/98 El Niño up to 29.0 °C) and in April 2002 (up to 25.0 °C). Temperatures more than 2 °C above monthly average were also observed in October 1994, in April 2000, and in November 2004. Significant trends were observed for oxygen concentration (increase) and several horizontal and vertical gradients. Among the copepods (72% of all individuals), the most abundant species were Paracalanus parvus (28%), Acartia tonsa (26%), and Calanus sp. (10%). The strong 1997–98 El Niño (EN) event led to drastic changes in species composition that were reversed during the 1998–99 La Niña (LN) event. Community parameters such as total abundance, diversity, equitability and species richness displayed marked variations associated with the 1997–98 EN and long-term trends. Long-term trends were significant for several vertical and horizontal temperature and oxygen gradients, indicating an increase in upwelling intensity at the shelf during the study period. 10-year-trends were also significant for total zooplankton abundance (increase) and community evenness (J, decline). Our data confirmed the importance of the weak EN in 2002/03 for the study region. Within the trend of increasing zooplankton abundance, a sharp step or shift was observed from 1999 to 2000. When using sequential t-tests to detect shifts in (x + 1) transformed abundance data, a significant rupture was found between the last sampling in 1999 and the first sampling in 2000. Also, a substantial decrease in diel variability occurred after 1999, probably due to changes in vertical migration patterns. The considerable increase in zooplankton abundance over the study period, the ENSO effect, and the 1999–2000 transition are discussed with regard to synchronicity with other zooplankton time series. The present study contributes with the first evidence from an important area located in the Humboldt Current for synchronous trends and changes that were previously observed elsewhere in the Pacific. Our results demonstrated the importance of long-term zooplankton monitoring studies in upwelling areas, and confirms the idea of dramatic changes in pelagic ecosystem structure occurring in the East Pacific.     

Dynamics of suprabenthos-zooplankton communities around the Balearic Islands (western Mediterranean): Influence of environmental variables and effects on the biological cycle of Aristeus antennatus

Dynamics of suprabenthos and zooplankton were analyzed in two areas located in the NW (off Sóller harbour) and S (off Cabrera Archipelago) of Mallorca (Balearic Islands, western Mediterranean) at depths ranging between 135–780 m. Four stations situated respectively at 150 m (shelf-slope break), and at bathyal depths of 350, 650 and 750 m were sampled at bi-monthly intervals during six cruises performed between August 2003 and June 2004. Suprabenthos showed maximum biomass in both areas from late spring to summer (April to August), while minimum biomass was found in autumn (September–November). Though variable, temporal dynamics of zooplankton showed peaks of biomass in late winter and summer (February and June), while minimals occurred in autumn (August–September) and, at bathyal depths, in April. Suprabenthos (abundance; MDS analyses) showed a sample aggregation as a function of depth (3 groups corresponding to the shelf-slope break, upper slope — over 350 m; and the middle, deeper part of the slope — over 650–750 m), without any separation of hauls by season. By contrast, zooplankton samples were separated by season and not by depth. There was evidence of three seasonal groups corresponding to summer (June 2004–August 2003), autumn–winter (September and November 2003, February 2004), and spring (April 2004), being especially well established off Sóller. In general, suprabenthos was significantly correlated with the sediment variables (e.g. total organic matter content (% OM), potential REDOX), whereas zooplankton was almost exclusively dependent on Chl a at the surface, which suggests two different food sources for suprabenthos and zooplankton. The increase of suprabenthos abundance in April–June was paralleled by a sharp increase (ca. 2.8 times) in the %OM on sediment during the same period, coupled ca. 1–2 months of delay with the peak of surface Chl a recorded in February–March (from satellite imagery data). Suprabenthos biomass was also correlated with salinity close to the bottom, suggesting a link between suprabenthos abundance and changes in the oceanographic condition of water masses close to the bottom. It is suggested that a higher suprabenthos biomass recorded off Sóller in comparison to that off Cabrera in June could, in turn, be related to a seasonal inflow of Levantine Intermediate Water (LIW) in April–June in this area at mid bathyal depths (350–650 m). This trend would be based on: 1) it was evident only at mid-slope depths between 350–750 m, coinciding with the LIW distribution, and 2) it was not recorded among zooplankton (collected throughout the water column). The possible effect of the fluctuations of suprabenthos and zooplankton on higher trophic levels has been explored studying the diet and food consumption rates of the red shrimp Aristeus antennatus, as indicator species by its dominance in bathyal communities. A. antennatus increased its food consumption from February to April–June 2004 off Sóller, which in the case of large (CL > 40 mm) specimens was found in both areas. In addition, there was a shift of diet from winter to spring–early summer. In this last period, A. antennatus preyed upon euphausiids and mesopelagic decapods and fish, while benthos (e.g. polychaetes and bivalves) decreased in the diet. This indicates an increase in the food consumption and probably in the caloric content of the diet in pre-spawning females in April–June 2004, which is synchronized with the period when gonad development begins in A. antennatus females (May–June). Anyway, macrozooplankton, and not suprabenthos, is crucial as a high energetic food source in the coupling between food intake and reproduction in the red shrimp.