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.     

A pronounced dry event recorded around 4.2 ka in brine sediments from the northern Red Sea

Partly laminated sediments were sampled from the brine-filled, anoxic Shaban Deep basin in the northern Red Sea. At about 4200 cal yr BP more than two millennia of anoxic sedimentation is replaced by a sub-oxic facies strongly suggesting the episodic absence of the brine. At the same time stable oxygen isotopes from surface dwelling foraminifera show a sharp increase (within less than 100 yr) pointing to a strong positive salinity anomaly at the sea surface. This major evaporation event significantly enhanced the renewal of deep water and the subsequent ventilation of the small Shaban Deep basin. The timing and strength of the reconstructed environmental changes around 4200 cal yr BP suggest that this event is the regional expression of a major drought event, which is widely observed in the neighboring regions, and which strongly affected Middle East agricultural civilizations.     

Effects of fracture lineaments and in-situ rock stresses on groundwater flow in hard rocks: a case study from Sunnfjord, western Norway

Systematic field mapping of fracture lineaments observed on aerial photographs shows that almost all of these structures are positively correlated with zones of high macroscopic and mesoscopic fracture frequencies compared with the surroundings. The lineaments are subdivided into zones with different characteristics: (1) a central zone with fault rocks, high fracture frequency and connectivity but commonly with mineral sealed fractures, and (2) a damage zone divided into a proximal zone with a high fracture frequency of lineament parallel, non-mineralized and interconnected fractures, grading into a distal zone with lower fracture frequencies and which is transitional to the surrounding areas with general background fracturing. To examine the possible relations between lineament architecture and in-situ rock stress on groundwater flow, the geological fieldwork was followed up by in-situ stress measurements and test boreholes at selected sites. Geophysical well logging added valuable information about fracture distribution and fracture flow at depths. Based on the studies of in-situ stresses as well as the lineaments and associated fracture systems presented above, two working hypotheses for groundwater flow were formulated: (i) In areas with a general background fracturing and in the distal zone of lineaments, groundwater flow will mainly occur along fractures parallel with the largest in-situ rock stress, unless fractures are critically loaded or reactivated as shear fractures at angles around 30° to σ_H; (ii) In the influence area of lineaments, the largest potential for groundwater abstraction is in the proximal zone, where there is a high fracture frequency and connectivity with negligible fracture fillings. The testing of the two hypotheses does not give a clear and unequivocal answer in support of the two assumptions about groundwater flow in the study area. But most of the observed data are in agreement with the predictions from the models, and can be explained by the action of the present stress field on pre-existing fractures.     

The BEAM-project: prediction and assessment of mixture toxicities in the aquatic environment

Freshwater and marine ecosystems are exposed to various multi-component mixtures of pollutants. Nevertheless, most ecotoxicological research and chemicals regulation focus on hazard and exposure assessment of individual substances only, the problem of chemical mixtures in the environment is ignored to a large extent. In contrast, the assessment of combination effects has a long tradition in pharmacology, where mixtures of chemicals are specifically designed to develop new products, e.g. human and veterinary drugs or agricultural and non-agricultural pesticides. In this area, two concepts are frequently used and are thought to describe fundamental relationships between single substance and mixture effects: Independent Action (Response Addition) and Concentration Addition. The question, to what extent these concepts may also be applied in an ecotoxicological and regulatory context may be considered a research topic of major importance, as the concepts would allow to make use of already existing single substance toxicity data for the predictive assessment of mixture toxicities. Two critical knowledge gaps are identified: (a) There is a lack of environmental realism, as a huge part of our current knowledge about the applicability of the concepts is restricted to artificial situations with respect to mixture composition or biological effect assessment. (b) The knowledge on what exactly is needed for using the concepts as tools for the predictive mixture toxicity assessment is insufficient. Both gaps seriously hamper the necessary, scientifically sound consideration of mixture toxicities in a regulatory context.In this paper, the two concepts will be briefly introduced, the necessity of considering the toxicities of chemical mixtures in the environment will be demonstrated and the applicability of Independent Action and Concentration Addition as tools for the prediction and assessment of mixture toxicities will be discussed. An overview of the specific aims and approaches of the BEAM project to fill in the identified knowledge gaps is given and first results are outlined.     

Description and evaluation of the bergen climate model: ARPEGE coupled with MICOM

A new coupled atmosphere–ocean–sea ice model has been developed, named the Bergen Climate Model (BCM). It consists of the atmospheric model ARPEGE/IFS, together with a global version of the ocean model MICOM including a dynamic–thermodynamic sea ice model. The coupling between the two models uses the OASIS software package. The new model concept is described, and results from a 300-year control integration is evaluated against observational data. In BCM, both the atmosphere and the ocean components use grids which can be irregular and have non-matching coastlines. Much effort has been put into the development of optimal interpolation schemes between the models, in particular the non-trivial problem of flux conservation in the coastal areas. A flux adjustment technique has been applied to the heat and fresh-water fluxes. There is, however, a weak drift in global mean sea-surface temperature (SST) and sea-surface salinity (SSS) of respectively 0.1 °C and 0.02 psu per century. The model gives a realistic simulation of the radiation balance at the top-of-the-atmosphere, and the net surface fluxes of longwave, shortwave, and turbulent heat fluxes are within observed values. Both global and total zonal means of cloud cover and precipitation are fairly close to observations, and errors are mainly related to the strength and positioning of the Hadley cell. The mean sea-level pressure (SLP) is well simulated, and both the mean state and the interannual standard deviation show realistic features. The SST field is several degrees too cold in the equatorial upwelling area in the Pacific, and about 1 °C too warm along the eastern margins of the oceans, and in the polar regions. The deviation from Levitus salinity is typically 0.1 psu – 0.4 psu, with a tendency for positive anomalies in the Northern Hemisphere, and negative in the Southern Hemisphere. The sea-ice distribution is realistic, but with too thin ice in the Arctic Ocean and too small ice coverage in the Southern Ocean. These model deficiencies have a strong influence on the surface air temperatures in these regions. Horizontal oceanic mass transports are in the lower range of those observed. The strength of the meridional overturning in the Atlantic is 18 Sv. An analysis of the large-scale variability in the model climate reveals realistic El Niño – Southern Oscillation (ENSO) and North Atlantic–Arctic Oscillation (NAO/AO) characteristics in the SLP and surface temperatures, including spatial patterns, frequencies, and strength. While the NAO/AO spectrum is white in SLP and red in temperature, the ENSO spectrum shows an energy maximum near 3 years.     

Progradational and backstepping shoreface deposits in the Ladinian to Early Norian Snadd Formation of the Barents Sea

Shoreface sandstone deposits within the Early Carnian part of the Snadd Formation of the Norwegian Barents Sea can be traced for hundreds of kilometres in the depositional strike direction and for tens of kilometres in the depositional‐dip direction. This study uses three‐dimensional seismic attribute mapping and two‐dimensional regional seismic profiles to visualize the seismic facies of these shoreface deposits and to map their internal stratigraphic architecture at a regional scale. The shoreface deposits are generally elongate but show variable width from north‐east to south‐west, which corresponds to a sediment source in the northern part of the basin and a southward decrease in longshore sediment transport. The Snadd Formation presents an example of how large‐scale progradational shoreface deposits develop. The linear nature of its shoreface deposits contrasts with more irregular, cuspate wave‐dominated deltaic shorelines that contain river outlets, and instead implies longshore drift as the main sediment source. In map view, discrete sets of linear features bounded by truncation surfaces scale directly to beach ridge sets in modern counterparts. The shoreface deposits studied here are characteristic in terms of scale and basin‐wide continuity, and offer insight into the contrast between shallow marine deposition under stable Triassic Greenhouse and fluctuating Holocene Icehouse climates. Findings presented herein are also important for hydrocarbon exploration in the Barents Sea, because they describe a hitherto poorly understood reservoir play in the Triassic interval, wherein the most prominent reservoir plays have so far been considered to be found in channelized deposits in net‐progradational delta‐plain strata that form the topsets to shelf‐edge clinoforms. The documented presence of widespread wave‐dominated shoreface deposits also has implications for how the relative importance of different sedimentary processes is considered within the basin during this period.     

Mercury monitoring: Mercury stability in bird feathers

The influence of ultraviolet light, heating, freezing and weathering on the mercury concentration in the primary feathers from Guillemot and Black Guillemot has been examined. Even within 8 months of exposure variation in mercury concentration due to either loss of mercury or weight loss of the feathers has been found to be less than 10% relative.     

Comment on “Is there a paleolimnological explanation for ‘walking on water’ in the Sea of Galilee?”

In a recent paper, Nof et al. (J Paleolimnol 35:417–439, 2006) suggest a physical mechanism which could account for the formation of ice on Lake Kinneret (Sea of Galilee) in northern Israel. Based on the sea surface temperature record of sediment cores from the Mediterranean Sea the authors argue that centennial-scale cold events had the potential to trigger local ‘springs ice’ formation on the lake in the past. Here, we demonstrate that a closer inspection of the paleoceanographic record in combination with correlation and regression analyses of meteorological data provides no evidence for such cold events in the lake region during the last 10,000 years. Thus, the formation of ‘springs ice’ on Lake Kinneret was unlikely at least since the beginning of the Neolithic.