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.     

探讨船舶触损海上养殖区事故的调查处理

海上养殖业也处于高速发展期,并成为部分地区重要的区域支柱产业,但随意建立或扩大养殖区域对船舶正常的海上航行安全构成严重的威胁,海上交通安全与海洋养殖的矛盾已日趋突出。本文介绍了海上养殖区的设立应具备的必要条件,并客观的提出了船舶触损海上养殖区事故调处的法律法规理解和处理方式。     

A 3D finite-element model of the Adriatic tides

A 3D finite-element numerical model is applied to the Adriatic Sea to simulate its tidal motions. This fully nonlinear model includes a free surface, very realistic topography, and an advanced turbulence closure. Comparison with available tidal elevations at coastal stations and with tidal ellipses at a few locations in the open sea demonstrates that the model simulations are highly accurate. The results are then used to determine the 3D distribution of the tidal residual currents.     

论上海北外滩航运服务功能的拓展

上海国际航运中心的发展和黄浦江的开发为北外滩航运服务功能的拓展提供了契机。本文分析了北外滩航运服务功能历史发展概况和现状 ,阐述了北外滩航运服务功能开发的机遇和挑战、优势和劣势 ,最后提出了北外滩航运服务功能发展的要求和建议。     

Matching biological traits to environmental conditions in marine benthic ecosystems

The effects of variability in environmental conditions on species composition in benthic ecosystems are well established, but relatively little is known about how environmental variability relates to ecosystem functioning. Benthic invertebrate assemblages are heavily involved in the maintenance of ecological processes and investigation of the biological characteristics (traits) expressed in these assemblages can provide information about some aspects of functioning. The aim of this study was to establish and explore relationships between environmental variability and biological traits expressed in megafauna assemblages in two UK regions. Patterns of trait composition were matched to environmental conditions and subsets of variables best describing these patterns determined. The nature of the relationships were subsequently examined at two separate scales, both between and within the regions studied. Over the whole area, some traits related to size, longevity, reproduction, mobility, flexibility, feeding method, sociability and living habit were negatively correlated with salinity, sea surface temperature, annual temperature range and the level of fishing effort, and positively associated with fish taxon richness and shell content of the substratum. Between the two regions, reductions in temperature range and shell content were associated with infrequent relative occurrences of short-lived, moderately mobile, flexible, solitary, opportunistic, permanent-burrow dwelling fauna and those exhibiting reproductive strategies based on benthic development. Relationships between some traits and environmental conditions diverged within the two regions, with increases in fishing effort and shell content of the substratum being associated with low frequencies of occurrence of moderately mobile and moderately to highly flexible fauna within one region, but high frequencies in the other. These changes in trait composition have implications for ecosystem processes, with, for example, reductions in permanent-burrow dwellers within one of the regions potentially compromising the ability of the assemblages to process and store chemicals and waste products. However, the connections between environmental conditions and trait composition are complex and incorporate many factors. Experimental investigations will be necessary to determine the extent and consequences of these important relationships.     

Study of the seasonal cycle of the biogeochemical processes in the Ligurian Sea using a 1D interdisciplinary model

A one-dimensional coupled physical–biogeochemical model has been built to study the pelagic food web of the Ligurian Sea (NW Mediterranean Sea). The physical model is the turbulent closure model (version 1D) developed at the GeoHydrodynamics and Environmental Laboratory (GHER) of the University of Liège. The ecosystem model contains 19 state variables describing the carbon and nitrogen cycles of the pelagic food web. Phytoplankton and zooplankton are both divided in three size-based compartments and the model includes an explicit representation of the microbial loop including bacteria, dissolved organic matter, nano-, and microzooplankton. The internal carbon/nitrogen ratio is assumed variable for phytoplankton and detritus, and constant for zooplankton and bacteria. Silicate is considered as a potential limiting nutrient of phytoplankton's growth. The aggregation model described by Kriest and Evans in (Proc. Ind. Acad. Sci., Earth Planet. Sci. 109 (4) (2000) 453) is used to evaluate the sinking rate of particulate detritus. The model is forced at the air–sea interface by meteorological data coming from the “Côte d'Azur” Meteorological Buoy. The dynamics of atmospheric fluxes in the Mediterranean Sea (DYFAMED) time-series data obtained during the year 2000 are used to calibrate and validate the biological model. The comparison of model results within in situ DYFAMED data shows that although some processes are not represented by the model, such as horizontal and vertical advections, model results are overall in agreement with observations and differences observed can be explained with environmental conditions.     

Pore water phosphate and ammonia below the permanent halocline in the south-eastern Baltic Sea and their benthic fluxes under anoxic conditions

In this paper the results of a study on the distribution of pore water phosphates and ammonia, and their fluxes under anoxic condition in a deep (> 70 m) accumulation-type bottom of the south-eastern Baltic Sea, namely in the Gdańsk Deep and the adjacent areas, are presented. All measurements were taken during the growth period, i.e. in September 2000, April 2001 and June 2002. Benthic phosphate and ammonia fluxes were estimated using Fick's First Law. Phosphate and ammonia concentrations ranged from 7.5 to 266.3 μmol dm^− 3 and from 53.6 to 1248.3 μmol dm^− 3, respectively. The values recorded in the central part of the Gdańsk Deep were lower than those found both on its slopes and on the SW slope of the Gotland Deep. The lowest phosphate contents were typical of the Oblique Sill which separates the Gdańsk and Gotland Deeps.In 1993–2002, as a result of anoxia the sediments in the Gdańsk Deep released about 5.1 × 10³ t P and 22.8 × 10³ t N. These loads supplied on average 1.5% and 0.9% of phytoplankton's demand for P and N, respectively. In comparison to the total external load of nutrients discharged to the Gulf of Gdańsk (i.e. 8.79 × 10³ t year^− 1 P_tot and 130.79 × 10³ t year^− 1 N_tot; [Witek, Z., Humborg, Ch., Savchuk, O., Grelowski, A. and Łysiak-Pastuszak, E., 2003. Nitrogen and phosphorus budgets of the Gulf of Gdańsk (Baltic Sea). Est. Coast. Shelf Sci., 57:239–248.]), the return flux of P and N from the anoxic sediments to the water column in the Gdańsk Deep was a minor source of these elements.     

Large-scale physical controls on phytoplankton growth in the Irminger Sea, Part II: Model study of the physical and meteorological preconditioning

The upper water column in the Irminger Sea is characterized by cold fresh arctic and subarctic waters and warm saline North Atlantic waters. In this study the local physical and meteorological preconditioning of the phytoplankton development over an annual cycle in the upper water column in four physical zones of the Irminger Sea is investigated. Data from four cruises of the UK's Marine Productivity programme are combined with results from a coupled biological–physical nitrogen–phytoplankton–zooplankton–detritus model run using realistic forcing. The observations and model predictions are compared and analyzed to identify the key parameters and processes which determine the observed heterogeneity in biological production in the Irminger Sea. The simulations show differences in the onset of the bloom, in the time of the occurrence of the maximum phytoplankton biomass and in the length of the bloom between the zones. The longest phytoplankton bloom of 90 days duration was predicted for the East Greenland Current of Atlantic origin zone. In contrast, for the Central Irminger Sea zone a phytoplankton bloom with a start at the beginning of May and the shortest duration of only 70 days was simulated. The latest onset of the phytoplankton bloom in mid May and the latest occurrence of the maximum biomass (end of July) were predicted for the Northern Irminger Current zone. Here the bloom lasted for 80 days. In contrast the phytoplankton bloom in the Southern Irminger Current zone started at the same time as in Central Irminger Sea, but peaked end of June and lasted for 80 days. For all four zones relatively low daily (0.3–0.5 g C m^− 2d^− 1) and annual primary production was simulated, ranging between 35.6 g C m^− 2y^− 1 in the East Greenland Current of Atlantic origin zone and 45.6 g C m^− 2y^− 1 in the Northern Irminger Current zone. The model successfully simulated the observed regional and spatial differences in terms of the maximum depth of winter mixing, the onset of stratification and the development of the seasonal thermocline, and the differences in biological characteristics between the zones. The initial properties of the water column and the seasonal cycle of physical and meteorological forcing in each of the zones are responsible for the observed differences during the Marine Productivity cruises. The timing of the transition from mixing to stratification regime, and the different prevailing light levels in each zone are identified as the crucial processes/parameters for the understanding of the dynamics of the pelagic ecosystem in the Irminger Sea.     

Assessment of wind quality for oceanographic modelling in semi-enclosed basins

The quality of surface winds derived from four meteorological models is assessed in the semi-enclosed Adriatic Sea over a 2-month period: a global hydrostatic model ECMWF T511 (40 km resolution), a hydrostatic limited area model LAMBO (20 km), and two non-hydrostatic limited area models: LAMI (7 km) and COAMPS™ (4 km). These wind models are used to drive a 2 km resolution wave model (SWAN) of the Adriatic, and wind and wave results are compared with observations at the ISMAR oceanographic tower off Venice. Waves are also compared at buoy locations near Ancona and Ortona. Consistently with earlier studies, the ECMWF fields underestimate the wind magnitude and do not reproduce the known spatial structure of strong wind events. The results show that the higher-resolution, limited area models LAMI and COAMPS exhibit better amplitude response than the coarser ECMWF: there is a 3- to 4-fold reduction of the wind underestimation at the platform (from 36% to 8–11%). The wave response is also improved with LAMI and COAMPS: there is a 2-fold reduction in the underestimation of wave heights at the platform. These non-hydrostatic models also produce wind fields with more realistic small-scale, spatial structure during strong wind events. The temporal correlation between observed and modelled wind, however, is highest with the global ECMWF model due to the fact that large-scale features can be predicted deterministically, whereas small-scale features can only be predicted stochastically. Models with less small-scale structure have better correlation because they have less “noise.” This explanation is supported by increased correlation between modelled and observed waves, the waves representing a smoothing of the wind over fetch and duration. Although there is room for improvement, the high-resolution, non-hydrostatic models (LAMI and COAMPS) offer significant advantages for driving oceanographic simulations in semi-enclosed basins such as the Adriatic Sea.     

The Mercator global ocean operational analysis system: Assessment and validation of an 11-year reanalysis

This paper presents Prototype Système 2 Global (PSY2G), the first Mercator global Ocean General Circulation Model (OGCM) to assimilate along-track sea level anomaly (SLA) satellite data. Based on a coarse resolution ocean model, this system was developed mainly for climatic purposes and will provide, for example, initial oceanic states for coupled ocean-atmosphere seasonal predictions. It has been operational since 3 September 2003 and produces an analysis and a two-week forecast for the global ocean every week. The PSY2G system uses an incremental assimilation scheme based on the Cooper and Haines [Cooper, M., Haines, K., 1996. Data assimilation with water property conservation. J. Geophys. Res., 101, 1059-1077.] lifting–lowering of isopycnals. The SLA increment is obtained using an optimal interpolation method then the correction is partitioned into baroclinic and barotropic contributions. The baroclinic ocean state correction consists of temperature, salinity and geostrophic velocity increments and the barotropic correction is a barotropic velocity increment. A reanalysis (1993–2003) was carried out that enabled the PSY2G system to perform its first operational cycle. All available SLA data sets (TOPEX/Poséïdon, ERS2, Geosat-Follow-On, Jason1 and Envisat) were assimilated for the 1993–2003 period. The major objective of this study is to assess the reanalysis from both an assimilation and a thermodynamic point of view in order to evaluate its realism, especially in the tropical band which is a key region for climatic studies. Although the system is also able to deliver forecasts, we have mainly focused on analysis. These results are useful because they give an a priori estimation of the qualities and capabilities of the operational ocean analysis system that has been implemented. In particular, the reanalysis identifies some regional biases in sea level variability such as near the Antarctic Circumpolar Current, in the eastern Equatorial Pacific and in the Norwegian Sea (generally less than 1 cm) with a small seasonal cycle. This is attributed to changes in mean circulation and vertical stratification caused by the assimilation methodology. But the model's low resolution, inaccurate physical parameterisations (especially for ocean–ice interactions) and surface atmospheric forcing also contribute to the occurrence of the SLA biases. A detailed analysis of the thermohaline structure of the ocean reveals that the isopycnal lifting–lowering tends to diffuse vertically the main thermocline. The impact on temperature is that the surface layer (0–200 m) becomes cooler whereas in deeper waters (from 500 to 1500 m), the ocean becomes slightly warmer. This is particularly true in the tropics, between 30°N and 30°S. However it can be demonstrated that the assimilation improves the variability in both surface currents and sub-surface temperature in the Equatorial Pacific Ocean.