A singular evolutive extended Kalman filter to assimilate ocean color data in a coupled physical–biochemical model of the North Atlantic ocean

Within the European DIADEM project, a data assimilation system for coupled ocean circulation and marine ecosystem models has been implemented for the North Atlantic and the Nordic Seas. One objective of this project is to demonstrate the relevance of sophisticated methods to assimilate satellite data such as altimetry, surface temperature and ocean color, into realistic ocean models. In this paper, the singular evolutive extended Kalman (SEEK) filter, which is an advanced assimilation scheme where three-dimensional, multivariate error statistics are taken into account, is used to assimilate ocean color data into the biological component of the coupled system. The marine ecosystem model, derived from the FDM model [J. Mar. Res. 48 (1990) 591], includes 11 nitrogen and carbon compartments and describes the synthesis of organic matter in the euphotic zone, its consumption by animals of upper trophic levels, and the recycling of detritic material in the deep ocean. The circulation model coupled to the ecosystem is the Miami isopycnic coordinate ocean model (MICOM), which covers the Atlantic and the Arctic Oceans with an enhanced resolution in the North Atlantic basin. The model is forced with realistic ECMWF ocean/atmosphere fluxes, which permits to resolve the seasonal variability of the circulation and mixed layer properties. In the twin assimimation experiments reported here, the predictions of the coupled model are corrected every 10 days using pseudo-measurements of surface phytoplankton as a substitute to chlorophyll concentrations measured from space. The diagnostics of these experiments indicate that the assimilation is feasible with a reduced-order Kalman filter of small rank (of order 10) as long as a sufficiently good identification of the error structure is available. In addition, the control of non-observed quantities such as zooplankton and nitrate concentrations is made possible, owing to the multivariate nature of the analysis scheme. However, a too severe truncation of the error sub-space downgrades the propagation of surface information below the mixed layer. The reduction of the actual state vector to the surface layers is therefore investigated to improve the estimation process in the perspective of sea-viewing wide field-of-view sensor (SeaWiFS) data assimilation experiments.     

The buckling of plastic oblate hemi-ellipsoidal dome shells under external hydrostatic pressure

The paper presents a theoretical and an experimental investigation into the buckling of seven oblate hemi-ellipsoidal dome shells under external hydrostatic pressure. Four of the shells were made in glass reinforced plastic and three were made from a thermosetting plastic called solid urethane plastic. All the vessels were tested to destruction. The theoretical study was made with the aid of a non-linear finite element solution, where both geometrical and material non-linearity were allowed for. Good agreement was found between experiment and theory for all the vessels. The very oblate domes failed axisymmetrically. Theoretical convergence was good for the more oblate domes but it was not as good as for the less oblate domes. This may have been because the less oblate domes did not fail in a classical axisymmetric manner as was expected. This work is of much importance in ocean engineering.     

Parameter identification of a compliant nonlinear SDOF system in random ocean waves by reverse MISO method

The determination of the drag and inertia coefficients, which enter into the wave force model given by Morison's equation, is particularly uncertain and difficult when a linear spectral model is used for ocean waves, and the structure is compliant and has nonlinear dynamic response. In this paper, a nonlinear System Identification method, called Reverse Multiple Inputs–Single Output (R–MISO) is applied to identify the hydrodynamic coefficients as well as the nonlinear stiffness parameter for a compliant single-degree-of-freedom system. Four different types of problems have been identified for use in various situations and the R–MISO has been applied to all of them. One of the problems requires iterative solution strategy to identify the parameters. The method has been found to be efficient in predicting the parameters with reasonable accuracy and has the potential for use in the laboratory experiments on compliant nonlinear offshore systems.     

Impact of Satellite Altimetry on Simulations of Sea Level Variability by an Indian Ocean Model

This article describes the impact of satellite altimeter data on the simulations of sea level variability (SLV) by a nonlinear reduced gravity model of the entire Indian Ocean. The model has been forced by 6-hourly analyzed wind stress data containing SSM/I observations and has been able to produce realistic circulation features. However, SLV values observed by Topex/Poseidon altimeter do not fit these simulations because of imperfect initial data. Hence an attempt has been made to initialize the model using altimeter data. The initialized model-generated SLVvalues have been compared with SLV derived by altimeter for monsoon as well as nonmonsoon months of 1996. Experimental runs have been performed for 10 days, 20 days, and one month. It has been found that the initialized model results on the final day of these experiments are in very good agreement with altimeter data of the same day. It is thus possible, in principle, to hindcast and forecast sea level variations in the time scale of 10 days to one month with the availability of good quality wind data for forcing the model and altimeter observations of sea level for initializing it.     

Directional spread parameter at intermediate water depth

The characteristics of directional spread parameters at intermediate water depth are investigated based on a cosine power ‘2s' directional spreading model. This is based on wave measurements carried out using a Datawell directional waverider buoy in 23 m water depth. An empirical equation for the frequency dependent directional spreading parameter is presented. Directional spreading function estimated based on the Maximum Entropy Method is compared with those obtained using a cosine power ‘2s' parameter model. A set of empirical equations relating the directional spreading parameter corresponding to the peak of wave spectrum to other wave parameters like significant wave height and period are obtained. It shows that the wave directional spreading at peak wave frequency can be related to the non-linearity parameter, which allows estimation of directional spreading without reference to wind information.     

Demersal Continental Shelf and Upper Slope Cephalopod Assemblages from the Balearic Sea (North-Western Mediterranean). Biological Aspects of Some Deep-Sea Species

The analysis of 79 hauls performed by commercial bottom trawlers from 50 to 800 m depth in the Balearic Sea (north-western Mediterranean) from June 1995 to September 1996 yielded a total of 30 cephalopod species belonging to 12 families. Cluster analysis of these data gave as a result two main groups 50–200 m and 200–800 m each subdivided into two other groups (50–100 vs. 100–200 m and 200–600 vs. 600–800 m). These results suggested the existence of two assemblages that could be associated to the continental shelf (50–100 m) and the upper slope (600–800 m) separated by a wide transitional zone (100–600 m) representing a region of overlapping shelf and slope faunas (ecotone). The faunistic bathymetric gradient showed a continuous substitution of species with depth rather than discrete assemblages separated by distinct boundaries. The more coastal species such as Eledone moschata, Loligo vulgarisSepia officinalis and Octopus vulgaris were found on the continental shelf; on the transitional zone, apart from species characteristic of this zone (Illex coindetii, Sepietta oweniana, Rossia macrosoma, Scaeurgus unicirrhus and Pteroctopus tetracirrhus), we also observed species from both the continental shelf and slope. The upper slope was characterized by typical species of deeper waters, such asBathypolypus sponsalis , Histioteuthis reversa, H. bonnellii, Ancistroteuthis lichtensteinii andOnychoteuthis banksii . The octopod O. vulgaris was the dominant species on the continental shelf and upper transitional zone, being substituted by T. sagittatus on the lower transitional zone and upper slope. Mean biomass decreased abruptly from the continental shelf to the transitional zone and from there to the upper slope. Mean species richness and species diversity were higher in the transitional zone than in the continental shelf and upper slope. Finally, some biological aspects of the more abundant deep-sea cephalopod species are studied: Bathypolypus sponsalis, Octopus salutii,Pteroctopus tetracirrhus , Histioteuthis reversa and H. bonnellii.     

Spatial and temporal variability of crustal magnetization of a slowly spreading ridge: Mid-Atlantic Ridge (20°–24° N)

The results of the two- and three-dimensional magnetic inversions performed on data located between 20°–24° N on the Mid-Atlantic Ridge indicate the crustal magnetization has decayed exponentially for the last 10 Ma, and that this decay has been fairly symmetric about the ridge axis. After removal of the mean temporal decay, the residual field is characterized by more positive magnetizations at the second-order discontinuities, regardless of initial magnetization direction. A model that involves the preferential emplacement of serpentinized lithologies near the discontinuities is proposed to explain this correlation. The temporal detrending method also indicates that several ridge-parallel depressions located on the flanks of the ridge axis are regions of more positive magnetizations. These bathymetric depressions may mark the locations of detachment faulting that occurred during amagmatic periods of extension. The general symmetry of the crustal magnetization about the ridge axis does not support the occurrence of continuous detachment faulting proposed to correspond to the inner and outer corners of ridge axis discontinuities.