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.

Wave-driven wind jets in the marine atmospheric boundary layer

The interaction between ocean surface waves and the overlying wind leads to a transfer of momentum across the air–sea interface. Atmospheric and oceanic models typically allow for momentum transfer to be directed only downward, from the atmosphere to the ocean. Recent observations have suggested that momentum can also be transferred upward when long wavelength waves, characteristic of remotely generated swell, propagate faster than the wind speed. The effect of upward momentum transfer on the marine atmospheric boundary layer is investigated here using idealized models that solve the momentum budget above the ocean surface. A variant of the classical Ekman model that accounts for the wave-induced stress demonstrates that, although the momentum flux due to the waves penetrates only a small fraction of the depth of the boundary layer, the wind profile is profoundly changed through its whole depth. When the upward momentum transfer from surface waves sufficiently exceeds the downward turbulent momentum flux, then the near-surface wind accelerates, resulting in a low-level wave-driven wind jet. This increases the Coriolis force in the boundary layer, and so the wind turns in the opposite direction to the classical Ekman layer. Calculations of the wave-induced stress due to a wave spectrum representative of fast-moving swell demonstrate upward momentum transfer that is dominated by contributions from waves in the vicinity of the peak in the swell spectrum. This is in contrast to wind-driven waves whose wave-induced stress is dominated by very short wavelength waves. Hence the role of swell can be characterized by the inverse wave age based on the wave phase speed corresponding to the peak in the spectrum. For a spectrum of waves, the total momentum flux is found to reverse sign and become upward, from waves to wind, when the inverse wave age drops below the range 0.15–0.2, which agrees reasonably well with previously published oceanic observations.

Enhancement of the aerosol direct radiative effect by semi-volatile aerosol components: airborne measurements in North-Western Europe

A case study of atmospheric aerosol measurements exploring the impact of the vertical distribution of aerosol chemical composition upon the radiative budget in North-Western Europe is presented. Sub-micron aerosol chemical composition was measured by an Aerodyne Aerosol Mass Spectrometer (AMS) on both an airborne platform and a ground-based site at Cabauw in the Netherlands. The examined period in May 2008 was characterised by enhanced pollution loadings in North-Western Europe and was dominated by ammonium nitrate and Organic Matter (OM). Both ammonium nitrate and OM were observed to increase with altitude in the atmospheric boundary layer. This is primarily attributed to partitioning of semi-volatile gas phase species to the particle phase at reduced temperature and enhanced relative humidity. Increased ammonium nitrate concentrations in particular were found to strongly increase the ambient scattering potential of the aerosol burden, which was a consequence of the large amount of associated water as well as the enhanced mass. During particularly polluted conditions, increases in aerosol optical depth of 50–100% were estimated to occur due to the observed increase in secondary aerosol mass and associated water uptake. Furthermore, the single scattering albedo was also shown to increase with height in the boundary layer. These enhancements combined to increase the negative direct aerosol radiative forcing by close to a factor of two at the median percentile level. Such increases have major ramifications for regional climate predictions as semi-volatile components are often not included in aerosol models. The results presented here provide an ideal opportunity to test regional and global representations of both the aerosol vertical distribution and subsequent impacts in North-Western Europe. North-Western Europe can be viewed as an analogue for the possible future air quality over other polluted regions of the Northern Hemisphere, where substantial reductions in sulphur dioxide emissions have yet to occur. Anticipated reductions in sulphur dioxide in polluted regions will result in an increase in the availability of ammonia to form ammonium nitrate as opposed to ammonium sulphate. This will be most important where intensive agricultural practises occur. Our observations over North-Western Europe, a region where sulphur dioxide emissions have already been reduced, indicate that failure to include the semi-volatile behaviour of ammonium nitrate will result in significant errors in predicted aerosol direct radiative forcing. Such errors will be particularly significant on regional scales.

Variation of the mixing state of Saharan dust particles with atmospheric transport

Mineral dust is an important aerosol species in the Earth’s atmosphere and has a major source within North Africa, of which the Sahara forms the major part. Aerosol Time of Flight Mass Spectrometry (ATOFMS) is first used to determine the mixing state of dust particles collected from the land surface in the Saharan region, showing low abundance of species such as nitrate and sulphate internally mixed with the dust mineral matrix. These data are then compared with the ATOFMS single particle mass spectra of Saharan dust particles detected in the marine atmosphere in the vicinity of the Cape Verde islands, which are further compared with those from particles with longer atmospheric residence sampled at a coastal station at Mace Head, Ireland. Saharan dust particles collected near the Cape Verde Islands showed increased internally mixed nitrate but no sulphate, whilst Saharan dust particles collected on the coast of Ireland showed a very high degree of internally mixed secondary species including nitrate, sulphate and methanesulphonate. This uptake of secondary species will change the pH and hygroscopic properties of the aerosol dust and thus can influence the budgets of other reactive gases, as well as influencing the radiative properties of the particles and the availability of metals for dissolution.

The income requirements of marine protected areas

Given the growing impact of human activities on the sea, managers are increasingly turning to marine protected areas (MPAs) to protect marine habitats and species. Many MPAs have been unsuccessful, however, and lack of income has been identified as a primary reason for failure. In this study, data from a global survey of 79 MPAs in 36 countries were analysed and attempts made to construct predictive models to determine the income requirements of any given MPA. Statistical tests were used to uncover possible patterns and relationships in the data, with two basic approaches. In the first of these, an attempt was made to build an explanatory "bottom-up" model of the cost structures that might be required to pursue various management activities. This proved difficult in practice owing to the very broad range of applicable data, spanning many orders of magnitude. In the second approach, a "top-down" regression model was constructed using logarithms of the base data, in order to address the breadth of the data ranges. This approach suggested that MPA size and visitor numbers together explained 46% of the minimum income requirements (P < 0.001), with area being the slightly more influential factor. The significance of area to income requirements was of little surprise, given its profile in the literature. However, the relationship between visitors and income requirements might go some way to explaining why northern hemisphere MPAs with apparently high incomes still claim to be under-funded. The relationship between running costs and visitor numbers has important implications not only in determining a realistic level of funding for MPAs, but also in assessing from where funding might be obtained. Since a substantial proportion of the income of many MPAs appears to be utilized for amenity purposes, a case may be made for funds to be provided from the typically better resourced government social and educational budgets as well as environmental budgets. Similarly visitor fees, already an important source of funding for some MPAs, might have a broader role to play in how MPAs are financed in the future. (C) 2007 Elsevier Ltd. All rights reserved.

Effects of forage source and soybean and marine algal oil supplementation on milk fatty acid composition in ewes

The objective of the present studies was to determine effects of basal dietary forage source on the response of milk fatty acid composition to an oil supplement based (2:1, respectively, w/w) on soybean oil and marine algae biomass oil high in cis-9, cis-12 C18:2n − 3 and C22:6n − 3, respectively. In Study 1, Hampshire × Dorset ewes (48) were randomly assigned to one of four treatments and 12 pens in a completely randomized design blocked on the basis of lambing date and number of lambs suckled. Control rations (60:40 forage:concentrate, dry matter (DM) basis) based on alfalfa pellets (AP) or corn silage (CS) were fed from lambing. Beginning at 22 days postpartum, three pens of ewes fed AP and three pens of ewes fed CS were supplemented with oil (30 g/kg of ration DM) in place of corn meal. Average ewe DM intake (DMI) and average daily gain (ADG) were measured weekly. Milk yield and composition were measured at 42 days postpartum. DMI was lower (P<0.02) for CS and for oil, but milk yield was not affected by forage source or oil supplementation. Milk fat content was higher for oil (P<0.10) and milk protein content was higher for AP (P<0.04). Total CLA concentration (g/100 g fatty acids) increased (P<0.01) with CS and oil, and the response to oil was greater for AP (P<0.04). Similarly, total trans-C18:1 and C22:6ω−3 concentrations were higher for CS and oil, but the response to oil was greater for CS (P<0.06 and P<0.01, respectively). In Study 2, the experiment was repeated using alfalfa haylage (AH) instead of AP. The DMI decreased (P<0.05) with oil feeding, but was not affected by forage source. Milk yield was decreased by feeding oil with AH, but not by feeding oil with CS (P<0.03). Milk fat content tended to be increased by feeding oil with AH, but tended to be decreased by feeding oil with CS (P<0.08). Total CLA concentration was increased (P<0.01) for AH versus CS and by oil, and the response to oil supplementation was greater for AH (P<0.01). In contrast, total trans-C18:1 concentration was higher for CS versus AH, with a greater response to oil for CS (P<0.05). Feeding marine oil increased the C22:6ω−3 (P<0.01) concentration, and the response was greater for AH (P<0.04). To further characterize the response of milk fat composition to dietary oil in ewes, a third study used six pens of three ewes each assigned to either the control CS diet used for Study 2 or the same diet supplemented with 45 g/kg (DM basis) of the oil mixture. Feeding oil had no effect on DMI, milk yield or milk fat concentration, but again increased (P<0.001) total trans-C18:1 and C22:6ω−3 concentrations and numerically increased (114%) total CLA concentration. Milk fatty acid composition responses to supplemental vegetable and marine oils were affected by forage source. Milk trans-C18:1 concentration was higher when CS was fed in Studies 1 and 2, but the effect of forage species on CLA concentration differed between studies, which may reflect differences in diet PUFA content and consumption, as well as amounts of dietary starch and fiber consumed. Despite large increases in trans-C18:1 concentration, milk fat content was not decreased by feeding unsaturated oils to ewes, even at diet levels of 45 g/kg of ration DM.

Genotypic sex determination enabled adaptive radiations of extinct marine reptiles

Adaptive radiations often follow the evolution of key traits, such as the origin of the amniotic egg and the subsequent radiation of terrestrial vertebrates. The mechanism by which a species determines the sex of its offspring has been linked to critical ecological and life-history traits(1-3) but not to major adaptive radiations, in part because sex-determining mechanisms do not fossilize. Here we establish a previously unknown coevolutionary relationship in 94 amniote species between sex-determining mechanism and whether a species bears live young or lays eggs. We use that relationship to predict the sex-determining mechanism in three independent lineages of extinct Mesozoic marine reptiles (mosasaurs, sauropterygians and ichthyosaurs), each of which is known from fossils to have evolved live birth(4-7). Our results indicate that each lineage evolved genotypic sex determination before acquiring live birth. This enabled their pelagic radiations, where the relatively stable temperatures of the open ocean constrain temperature-dependent sex determination in amniote species. Freed from the need to move and nest on land(4,5,8), extreme physical adaptations to a pelagic lifestyle evolved in each group, such as the fluked tails, dorsal fins and wing-shaped limbs of ichthyosaurs. With the inclusion of ichthyosaurs, mosasaurs and sauropterygians, genotypic sex determination is present in all known fully pelagic amniote groups (sea snakes, sirenians and cetaceans), suggesting that this mode of sex determination and the subsequent evolution of live birth are key traits required for marine adaptive radiations in amniote lineages.

Low genetic diversity in a marine nature reserve: re-evaluating diversity criteria in reserve design

Little consideration has been given to the genetic composition of populations associated with marine reserves, as reserve designation is generally to protect specific species, communities or habitats. Nevertheless, it is important to conserve genetic diversity since it provides the raw material for the maintenance of species diversity over longer, evolutionary time-scales and may also confer the basis for adaptation to environmental change. Many current marine reserves are small in size and isolated to some degree (e.g. sea loughs and offshore islands). While such features enable easier management, they may have important implications for the genetic structure of protected populations, the ability of populations to recover from local catastrophes and the potential for marine reserves to act as sources of propagules for surrounding areas. Here, we present a case study demonstrating genetic differentiation, isolation, inbreeding and reduced genetic diversity in populations of the dogwhelk Nucella lapillus in Lough Hyne Marine Nature Reserve (an isolated sea lough in southern Ireland), compared with populations on the local adjacent open coast and populations in England, Wales and France. Our study demonstrates that this sea lough is isolated from open coast populations, and highlights that there may be long-term genetic consequences of selecting reserves on the basis of isolation and ease of protection.

Computer modelling and estimation of recruitment patterns of non-branching coral colonies at three sites in Wakatobi Marine Park, S.E. Sulawesi, Indonesia; implications for coral reef conservation

We have studied growth and estimated recruitment of massive coral colonies at three sites, Kaledupa, Hoga and Sampela, separated by about 1.5 km in the Wakatobi Marine National Park, S.E. Sulawesi, Indonesia. There was significantly higher species richness (P<0.05), coral cover (P<0.05) and rugosity (P<0.01) at Kaledupa than at Sampela. A model for coral reef growth has been developed based on a rational polynomial function, where dx/dt is an index of coral growth with time; W is the variable (for example, coral weight, coral length or coral area), up to the power of n in the numerator and m in the denominator; a1……an and b1…bm are constants. The values for n and m represent the degree of the polynomial, and can relate to the morphology of the coral. The model was used to simulate typical coral growth curves, and tested using published data obtained by weighing coral colonies underwater in reefs on the south-west coast of Curaçao [‘Neth. J. Sea Res. 10 (1976) 285’]. The model proved an accurate fit to the data, and parameters were obtained for a number of coral species. Surface area data was obtained on over 1200 massive corals at three different sites in the Wakatobi Marine National Park, S.E. Sulawesi, Indonesia. The year of an individual's recruitment was calculated from knowledge of the growth rate modified by application of the rational polynomial model. The estimated pattern of recruitment was variable, with little numbers of massive corals settling and growing before 1950 at the heavily used site, Sampela, relative to the reef site with little or no human use, Kaledupa, and the intermediate site, Hoga. There was a significantly greater sedimentation rate at Sampela than at either Kaledupa (P<0.0001) or Hoga (P<0.0005). The relative mean abundance of fish families present at the reef crests at the three sites, determined using digital video photography, did not correlate with sedimentation rates, underwater visibility or lack of large non-branching coral colonies. Radial growth rates of three genera of non-branching corals were significantly lower at Sampela than at Kaledupa or at Hoga, and there was a high correlation (r=0.89) between radial growth rates and underwater visibility. Porites spp. was the most abundant coral over all the sites and at all depths followed by Favites (P<0.04) and Favia spp. (P<0.03). Colony ages of Porites corals were significantly lower at the 5 m reef flat on the Sampela reef than at the same depth on both other reefs (P<0.005). At Sampela, only 2.8% of corals on the 5 m reef crest are of a size to have survived from before 1950. The Scleractinian coral community of Sampela is severely impacted by depositing sediments which can lead to the suffocation of corals, whilst also decreasing light penetration resulting in decreased growth and calcification rates. The net loss of material from Sampela, if not checked, could result in the loss of this protective barrier which would be to the detriment of the sublittoral sand flats and hence the Sampela village.

Simulating the formation of secondary organic aerosol from the photooxidation of aromatic hydrocarbons

The formation and composition of secondary organic aerosol (SOA) from the photooxidation of benzene, p-xylene, and 1,3,5-trimethylbenzene has been simulated using the Master Chemical Mechanism version 3.1 (MCM v3.1) coupled to a representation of the transfer of organic material from the gas to particle phase. The combined mechanism was tested against data obtained from a series of experiments conducted at the European Photoreactor (EUPHORE) outdoor smog chamber in Valencia, Spain. Simulated aerosol mass concentrations compared reasonably well with the measured SOA data only after absorptive partitioning coefficients were increased by a factor of between 5 and 30. The requirement of such scaling was interpreted in terms of the occurrence of unaccounted-for association reactions in the condensed organic phase leading to the production of relatively more nonvolatile species. Comparisons were made between the relative aerosol forming efficiencies of benzene, toluene, p-xylene, and 1,3,5-trimethylbenzene, and differences in the OH-initiated degradation mechanisms of these aromatic hydrocarbons. A strong, nonlinear relationship was observed between measured (reference) yields of SOA and (proportional) yields of unsaturated dicarbonyl aldehyde species resulting from ring-fragmenting pathways. This observation, and the results of the simulations, is strongly suggestive of the involvement of reactive aldehyde species in association reactions occurring in the aerosol phase, thus promoting SOA formation and growth. The effect of NO, concentrations on SOA formation efficiencies (and formation mechanisms) is discussed.

Simulating the detailed chemical composition of secondary organic aerosol formed on a regional scale during the TORCH 2003 campaign in the southern UK

Following on from the companion study (Johnson et al., 2006), a photochemical trajectory model (PTM) has been used to simulate the chemical composition of organic aerosol for selected events during the 2003 TORCH (Tropospheric Organic Chemistry Experiment) field campaign. The PTM incorporates the speciated emissions of 124 nonmethane anthropogenic volatile organic compounds (VOC) and three representative biogenic VOC, a highly-detailed representation of the atmospheric degradation of these VOC, the emission of primary organic aerosol (POA) material and the formation of secondary organic aerosol (SOA) material. SOA formation was represented by the transfer of semi and non-volatile oxidation products from the gas-phase to a condensed organic aerosol-phase, according to estimated thermodynamic equilibrium phase-partitioning characteristics for around 2000 reaction products. After significantly scaling all phase-partitioning coefficients, and assuming a persistent background organic aerosol (both required in order to match the observed organic aerosol loadings), the detailed chemical composition of the simulated SOA has been investigated in terms of intermediate oxygenated species in the Master Chemical Mechanism, version 3.1 ( MCM v3.1). For the various case studies considered, 90% of the simulated SOA mass comprises between ca. 70 and 100 multifunctional oxygenated species derived, in varying amounts, from the photooxidation of VOC of anthropogenic and biogenic origin. The anthropogenic contribution is dominated by aromatic hydrocarbons and the biogenic contribution by alpha-and beta-pinene (which also constitute surrogates for other emitted monoterpene species). Sensitivity in the simulated mass of SOA to changes in the emission rates of anthropogenic and biogenic VOC has also been investigated for 11 case study events, and the results have been compared to the detailed chemical composition data. The role of accretion chemistry in SOA formation, and its implications for the results of the present investigation, is discussed.

Simulating regional scale secondary organic aerosol formation during the TORCH 2003 campaign in the southern UK

A photochemical trajectory model has been used to simulate the chemical evolution of air masses arriving at the TORCH field campaign site in the southern UK during late July and August 2003, a period which included a widespread and prolonged photochemical pollution episode. The model incorporates speciated emissions of 124 nonmethane anthropogenic VOC and three representative biogenic VOC, coupled with a comprehensive description of the chemistry of their degradation. A representation of the gas/aerosol absorptive partitioning of ca. 2000 oxygenated organic species generated in the Master Chemical Mechanism (MCM v3.1) has been implemented, allowing simulation of the contribution to organic aerosol (OA) made by semi- and non-volatile products of VOC oxidation; emissions of primary organic aerosol (POA) and elemental carbon (EC) are also represented. Simulations of total OA mass concentrations in nine case study events (optimised by comparison with observed hourly-mean mass loadings derived from aerosol mass spectrometry measurements) imply that the OA can be ascribed to three general sources: (i) POA emissions; (ii) a '' ubiquitous '' background concentration of 0.7 mu g m(-3); and (iii) gas-to-aerosol transfer of lower volatility products of VOC oxidation generated by the regional scale processing of emitted VOC, but with all partitioning coefficients increased by a species-independent factor of 500. The requirement to scale the partitioning coefficients, and the implied background concentration, are both indicative of the occurrence of chemical processes within the aerosol which allow the oxidised organic species to react by association and/or accretion reactions which generate even lower volatility products, leading to a persistent, non-volatile secondary organic aerosol (SOA). The contribution of secondary organic material to the simulated OA results in significant elevations in the simulated ratio of organic carbon (OC) to EC, compared with the ratio of 1.1 assigned to the emitted components. For the selected case study events, [OC]/[EC] is calculated to lie in the range 2.7-9.8, values which are comparable with the high end of the range reported in the literature.

Oxidation of biogenic and water-soluble compounds in aqueous and organic aerosol droplets by ozone: a kinetic and product analysis approach using laser Raman tweezers

The results of an experimental study into the oxidative degradation of proxies for atmospheric aerosol are presented. We demonstrate that the laser Raman tweezers method can be used successfully to obtain uptake coeffcients for gaseous oxidants on individual aqueous and organic droplets, whilst the size and composition of the droplets is simultaneously followed. A laser tweezers system was used to trap individual droplets containing an unsaturated organic compound in either an aqueous or organic ( alkane) solvent. The droplet was exposed to gas- phase ozone and the reaction kinetics and products followed using Raman spectroscopy. The reactions of three different organic compounds with ozone were studied: fumarate anions, benzoate anions and alpha pinene. The fumarate and benzoate anions in aqueous solution were used to represent components of humic- like substances, HULIS; a alpha- pinene in an alkane solvent was studied as a proxy for biogenic aerosol. The kinetic analysis shows that for these systems the diffusive transport and mass accommodation of ozone is relatively fast, and that liquid- phase di. ffusion and reaction are the rate determining steps. Uptake coe. ffcients, g, were found to be ( 1.1 +/- 0.7) x 10(-5), ( 1.5 +/- 0.7) x 10 (-5) and ( 3.0 - 7.5) x 10 (-3) for the reactions of ozone with the fumarate, benzoate and a- pinene containing droplets, respectively. Liquid- phase bimolecular rate coe. cients for reactions of dissolved ozone molecules with fumarate, benzoate and a- pinene were also obtained: k(fumarate) = ( 2.7 +/- 2) x 10 (5), k(benzoate) = ( 3.5 +/- 3) x 10 (5) and k(alpha-pinene) = ( 1-3) x 10(7) dm(3) mol (-1) s (- 1). The droplet size was found to remain stable over the course of the oxidation process for the HULIS- proxies and for the oxidation of a- pinene in pentadecane. The study of the alpha- pinene/ ozone system is the first using organic seed particles to show that the hygroscopicity of the particle does not increase dramatically over the course of the oxidation. No products were detected by Raman spectroscopy for the reaction of benzoate ions with ozone. One product peak, consistent with aqueous carbonate anions, was observed when following the oxidation of fumarate ions by ozone. Product peaks observed in the reaction of ozone with alpha- pinene suggest the formation of new species containing carbonyl groups.