Data assimilation for marine monitoring and prediction: The MERCATOR operational assimilation systems and the MERSEA developments

During the past 15 years, a number of initiatives have been undertaken at national level to develop ocean forecasting systems operating at regional and/or global scales. The co-ordination between these efforts has been organized internationally through the Global Ocean Data Assimilation Experiment (GODAE). The French MERCATOR project is one of the leading participants in GODAE. The MERCATOR systems routinely assimilate a variety of observations such as multi-satellite altimeter data, sea-surface temperature and in situ temperature and salinity profiles, focusing on high-resolution scales of the ocean dynamics. The assimilation strategy in MERCATOR is based on a hierarchy of methods of increasing sophistication including optimal interpolation, Kalman filtering and variational methods, which are progressively deployed through the Syst`eme d’Assimilation MERCATOR (SAM) series. SAM-1 is based on a reduced-order optimal interpolation which can be operated using ‘altimetry-only’ or ‘multi-data’ set-ups; it relies on the concept of separability, assuming that the correlations can be separated into a product of horizontal and vertical contributions. The second release, SAM-2, is being developed to include new features from the singular evolutive extended Kalman (SEEK) filter, such as three-dimensional, multivariate error modes and adaptivity schemes. The third one, SAM-3, considers variational methods such as the incremental four-dimensional variational algorithm. Most operational forecasting systems evaluated during GODAE are based on least-squares statistical estimation assuming Gaussian errors. In the framework of the EU MERSEA (Marine EnviRonment and Security for the European Area) project, research is being conducted to prepare the next-generation operational ocean monitoring and forecasting systems. The research effort will explore nonlinear assimilation formulations to overcome limitations of the current systems. This paper provides an overview of the developments conducted in MERSEA with the SEEK filter, the Ensemble Kalman filter and the sequential importance re-sampling filter.

Conservative upwind difference schemes for the Euler equations for real gas flows in a duct

In a recent paper [P. Glaister, Conservative upwind difference schemes for compressible flows in a Duct, Comput. Math. Appl. 56 (2008) 1787–1796] numerical schemes based on a conservative linearisation are presented for the Euler equations governing compressible flows of an ideal gas in a duct of variable cross-section, and in [P. Glaister, Conservative upwind difference schemes for compressible flows of a real gas, Comput. Math. Appl. 48 (2004) 469–480] schemes based on this philosophy are presented for real gas flows with slab symmetry. In this paper we seek to extend these ideas to encompass compressible flows of real gases in a duct. This will incorporate the handling of additional terms arising out of the variable geometry and the non-ideal nature of the gas.

Mapping and controlling molecular processes: The application of control theory and system identification algorithms to ASOPS and other fast pulse spectroscopies

Asynchronous Optical Sampling (ASOPS) [1,2] and frequency comb spectrometry [3] based on dual Ti:saphire resonators operated in a master/slave mode have the potential to improve signal to noise ratio in THz transient and IR sperctrometry. The multimode Brownian oscillator time-domain response function described by state-space models is a mathematically robust framework that can be used to describe the dispersive phenomena governed by Lorentzian, Debye and Drude responses. In addition, the optical properties of an arbitrary medium can be expressed as a linear combination of simple multimode Brownian oscillator functions. The suitability of a range of signal processing schemes adopted from the Systems Identification and Control Theory community for further processing the recorded THz transients in the time or frequency domain will be outlined [4,5]. Since a femtosecond duration pulse is capable of persistent excitation of the medium within which it propagates, such approach is perfectly justifiable. Several de-noising routines based on system identification will be shown. Furthermore, specifically developed apodization structures will be discussed. These are necessary because due to dispersion issues, the time-domain background and sample interferograms are non-symmetrical [6-8]. These procedures can lead to a more precise estimation of the complex insertion loss function. The algorithms are applicable to femtosecond spectroscopies across the EM spectrum. Finally, a methodology for femtosecond pulse shaping using genetic algorithms aiming to map and control molecular relaxation processes will be mentioned.

Evaluating the efficiency of sampling methods in assessing soil macrofauna communities in arable systems

The soil fauna is often a neglected group in many large-scale studies of farmland biodiversity due to difficulties in extracting organisms efficiently from the soil. This study assesses the relative efficiency of the simple and cheap sampling method of handsorting against Berlese-Tullgren funnel and Winkler apparatus extraction. Soil cores were taken from grassy arable field margins and wheat fields in Cambridgeshire, UK, and the efficiencies of the three methods in assessing the abundances and species densities of soil macroinver-tebrates were compared. Handsorting in most cases was as efficient at extracting the majority of the soil macrofauna as the Berlese-Tullgren funnel and Winkler bag methods, although it underestimated the species densities of the woodlice and adult beetles. There were no obvious biases among the three methods for the particular vegetation types sampled and no significant differences in the size distributions of the earthworms and beetles. Proportionally fewer damaged earthworms were recorded in larger (25 x 25 cm) soil cores when compared with smaller ones (15 x 15 cm). Handsorting has many benefits, including targeted extraction, minimum disturbance to the habitat and shorter sampling periods and may be the most appropriate method for studies of farmland biodiversity when a high number of soil cores need to be sampled. (C) 2008 Elsevier Masson SAS. All rights reserved.

The impact of two arable field margin management schemes on litter decomposition

A primary objective of agri-environment schemes is the conservation of biodiversity; in addition to increasing the value of farmland for wildlife, these schemes also aim to restore natural ecosystem functioning. The management of scheme options can influence their value for delivering ecosystem services by modifying the composition of floral and faunal communities. This study examines the impact of an agri-environment scheme prescription on ecosystem functioning by testing the hypothesis that vegetation management influences decomposition rates in grassy arable field margins. The effects of two vegetation management practices in arable field margins - cutting and soil disturbance (scarification) - on litter decomposition were compared using a litterbag experimental approach in early April 2006. Bags had either small mesh designed to restrict access to soil macrofauna, or large mesh that would allow macrofauna to enter. Bags were positioned on the soil surface or inserted into the soil in cut and scarified margins, retrieved after 44, 103 and 250 days and the amount of litter mass remaining was calculated. Litter loss from the litterbags with large mesh was greater than from the small mesh bags, providing evidence that soil macrofauna accelerate rates of litter decomposition. In the large mesh bags, the proportion of litter remaining in bags above and belowground in the cut plots was similar, while in the scarified plots, there was significantly more litter left in the aboveground bags than in the belowground bags. This loss of balance between decomposition rates above and belowground in scarified margins may have implications for the development and maintenance of grassy arable field margins by influencing nutrient availability for plant communities. (C) 2008 Elsevier B.V. All rights reserved.

Measuring bee diversity in different European habitats and biogeographical regions

Bee pollinators are currently recorded with many different sampling methods. However, the relative performances of these methods have not been systematically evaluated and compared. In response to the strong need to record ongoing shifts in pollinator diversity and abundance, global and regional pollinator initiatives must adopt standardized sampling protocols when developing large-scale and long-term monitoring schemes. We systematically evaluated the performance of six sampling methods (observation plots, pan traps, standardized and variable transect walks, trap nests with reed internodes or paper tubes) that are commonly used across a wide range of geographical regions in Europe and in two habitat types (agricultural and seminatural). We focused on bees since they represent the most important pollinator group worldwide. Several characteristics of the methods were considered in order to evaluate their performance in assessing bee diversity: sample coverage, observed species richness, species richness estimators, collector biases (identified by subunit-based rarefaction curves), species composition of the samples, and the indication of overall bee species richness (estimated from combined total samples). The most efficient method in all geographical regions, in both the agricultural and seminatural habitats, was the pan trap method. It had the highest sample coverage, collected the highest number of species, showed negligible collector bias, detected similar species as the transect methods, and was the best indicator of overall bee species richness. The transect methods were also relatively efficient, but they had a significant collector bias. The observation plots showed poor performance. As trap nests are restricted to cavity-nesting bee species, they had a naturally low sample coverage. However, both trap nest types detected additional species that were not recorded by any of the other methods. For large-scale and long-term monitoring schemes with surveyors with different experience levels, we recommend pan traps as the most efficient, unbiased, and cost-effective method for sampling bee diversity. Trap nests with reed internodes could be used as a complementary sampling method to maximize the numbers of collected species. Transect walks are the principal method for detailed studies focusing on plant-pollinator associations. Moreover, they can be used in monitoring schemes after training the surveyors to standardize their collection skills.

A jackknife variance estimator for unequal probability sampling

The jackknife method is often used for variance estimation in sample surveys but has only been developed for a limited class of sampling designs.We propose a jackknife variance estimator which is defined for any without-replacement unequal probability sampling design. We demonstrate design consistency of this estimator for a broad class of point estimators. A Monte Carlo study shows how the proposed estimator may improve on existing estimators.

Variance estimation with highly stratified sampling designs with unequal probabilities

It is common practice to design a survey with a large number of strata. However, in this case the usual techniques for variance estimation can be inaccurate. This paper proposes a variance estimator for estimators of totals. The method proposed can be implemented with standard statistical packages without any specific programming, as it involves simple techniques of estimation, such as regression fitting.

Variance estimation for systematic sampling from deliberately ordered populations

The systematic sampling (SYS) design (Madow and Madow, 1944) is widely used by statistical offices due to its simplicity and efficiency (e.g., Iachan, 1982). But it suffers from a serious defect, namely, that it is impossible to unbiasedly estimate the sampling variance (Iachan, 1982) and usual variance estimators (Yates and Grundy, 1953) are inadequate and can overestimate the variance significantly (Särndal et al., 1992). We propose a novel variance estimator which is less biased and that can be implemented with any given population order. We will justify this estimator theoretically and with a Monte Carlo simulation study.

Variance estimation with Chao's sampling scheme

We show that the Hájek (Ann. Math Statist. (1964) 1491) variance estimator can be used to estimate the variance of the Horvitz–Thompson estimator when the Chao sampling scheme (Chao, Biometrika 69 (1982) 653) is implemented. This estimator is simple and can be implemented with any statistical packages. We consider a numerical and an analytic method to show that this estimator can be used. A series of simulations supports our findings.