Developing a Marine Information System by Integrating Existing Ocean Models Using Object-Oriented Technology

Information systems developed for different applications within the environmental domain have common characteristics, which can potentially be abstracted for sharing and reuse of design and software modules. This article presents an approach to designing for reuse by abstracting commonalities in the design of a Marine Information System (MIS), facilitating data management in a prototype operational monitoring, forecasting, and management system for the North Atlantic and the Nordic Seas. A detailed study of the requirements and data analysis was carried out, and Object-Oriented Technology (OOT) is employed to encapsulate abstractions and to promote reuse of code and design. This article identifies the Object-Oriented Frameworks (OOFW) required to build the MIS. It also provides guidelines to environmental scientists for restructuring legacy software and employing modern programming techniques.     

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.     

Rare-earth abundances in chondritic meteorites

Fifteen chondrites, including eight carbonaceous chondrites, have been analyzed for rare earth element (REE) abundances by isotope dilution. These analyses complement and extend earlier isotope dilution REE determinations in chondrites, performed in other laboratories, so that coverage of major chondrite classes is now complete. An examination of this body of precise and comparable REE data from individual chondrites reveals that only a small proportion of the analyses have flat, unfractionated REE patterns within experimental error. A statistical procedure is used to derive revised chondritic abundances of REE by selection of unfractionated patterns. A number of the remaining analyses show Eu anomalies and fractionated patterns consistent with magmatic fractionation as encountered in the products of planetary differentiation. However, many patterns exhibit features not readily explicable by known magmatic processes; in particular, positive Ce anomalies are often encountered. Abundance anomalies can be quantitatively determined by the use of a least-squares curve fitting procedure. The wide variety of anomalous patterns and the uncertainties in model parameters preclude detailed modeling of the origin of anomalies, but it is probable that at least some arise from fractional condensation in the solar nebula, as has been demonstrated for Allende inclusions. Elemental abundance anomalies are found in all major chondrite classes. If these anomalies are ignored, the range and nature of variation within chondrite classes are consistent with a parent body model, in which solid-liquid or solid-solid equilibria induce variations from an unfractionated bulk composition. Absolute abundances in the H, L and LL parent bodies are almost twice those of the E parent body.The persistence of anomalies in chondritic materials relatively removed from direct condensational processes implies that anomalous components are resistant to equilibration or were introduced at a late stage of chondrite formation. Large scale segregation of gas and condensate is also implied, and raises the possibility of bulk variations in REE abundances between planetary bodies.     

Experimental silicate mineral/melt partition coefficients for beryllium and the crustal Be cycle from migmatite to pegmatite

Partition coefficients (D_Be^mineral/melt) for beryllium between hydrous granitic melt and alkali feldspars, plagioclase feldspars, quartz, dark mica, and white mica were determined by experiment at 200 MPa H₂O as a function of temperature (650-900°C), activity of Be in melt (trace levels to beryl saturation), bulk composition, and thermal run direction. At trace levels, Be is compatible in plagioclase of An₃₁ (1.84 at 700°C) and muscovite (1.35 at 700°C) but incompatible in biotite (0.39-0.54 from 650-800°C), alkali feldspar (0.38-0.19 from 680-850°C), quartz (0.24 at 800°C), and albite (0.10 at 750°C). The partition coefficients are different at saturation of the melt in beryl: lower in the case of plagioclase of An₃₁ (0.89 at 700°C), muscovite (0.87 at 700°C), biotite (0.18-0.08 from 675-800°C), alkali feldspar (0.18-0.14 from 680-700°C), and quartz (0.17-0.08 from 750-800°C), but higher in the case of albite (0.37 at 750°C).With other data sources, these new partition coefficients were utilized to track, first, the distribution of Be between aluminous quartzofeldspathic source rocks and their anatectic melts, and second, the dispersion or concentration of Be in melt through igneous crystal fractionation of different magma types (e.g., S-type, I-type) up to beryl-saturated granitic pegmatites and, finally, into their hydrothermal aureoles. Among the rock-forming minerals, cordierite, calcic oligoclase, and muscovite (in this order) control the fate of Be because of the compatibility of Be in these phases. In general, beryl-bearing pegmatites can arise only after extended crystal fractionation of large magma batches (to F, fraction of melt remaining, ≤0.05); granitic magmas that originate from cordierite-bearing protoliths or that contain large modal quantities of calcic oligoclase will not achieve beryl saturation at any point in their evolution.     

Experimental partitioning of Be,Cs, and other trace elements between cordierite and felsic melt,and the chemical signature of S-type granite

The trace-element signature that cordierite (Crd) imparts to silicic magmas was evaluated by experiment using metapelite mineral mixtures to produce cordierite-bearing peraluminous granitic melts at 200 MPa (P_H2O), from 700 to 850 °C. Most elemental partition coefficients vary with T. Beryllium is strongly compatible, with D_Be^Crd/melt values decreasing linearly from 202.0 to 6.7 as T rises from 700 to 850 °C. Manganese is compatible (D_Mn^Crd/melt=7.67 to 1.92 over the same range of T), and shows similar values to those reported for biotite in silicic melts. Incompatible components include Li, Rb, B, F and P, although Cs is nearly compatible in cordierite, especially at higher T (D_Cs^Crd/melt=~0.19 to 0.60) where the large alkalis are better accommodated structurally. Cordierite appears to be the most effective crystalline reservoir of Be and Cs in metapelites and their anatectic melts. Natural data support the hypothesis that Crd, when present in granitic melts, sequesters Be, Cs and, in the absence of garnet, Mn. S-type granitic rocks containing Crd show consistently low Be contents (mean=0.8 ppm Be with an average range of <1 to 1.20) whereas Crd-free granites (e.g., containing accessory garnet) exhibit distinctly higher Be contents (mean=6 ppm Be with an average range of 3 to 12). These values increase further in evolved facies (mean=69 ppm Be with a an average range of 11 to 145) which commonly give rise to beryl-bearing pegmatites. Whole-rock signatures of Be discriminate source environments of silicic magmas at a resolution equal to the boundaries of the cordierite stability field - e.g., at the P-T-X conditions where cordierite gives rise to garnet+aluminum silicate. Cordierite-bearing granitic rocks contain low Cs contents (mean=1.8 ppm Cs) compared to the Crd-free equivalents (mean=18 ppm Cs). Mn contents also correlate with the presence (mean=0.01 wt% MnO) or absence of Crd (mean=0.09 wt% MnO). Depending on its contribution to anatexis, cordierite may either give or take S-type chemical character from granitic liquids, resulting in a distinctive Crd-associated group of S-type elements. This signature is different from that of micas (high Li, F and, to a lesser degree, Be and Mn). Whole-rock compositions of granites, coupled with notable absences of beryl in their associated pegmatites, indicate that a sizable population of S-type granites originated from Crd-bearing sources. The normative Crd component of silicic peraluminous melts is Д wt% to 850 °C. Higher modal contents of cordierite reflect either restite entrainment or peritectic reactions which produce Crd after magma ascent to shallow depths. The distinctive trace-element signature of cordierite now provides improved resolution of the source mineralogy for S-type magmas.     

Assimilation of SLA and SST data into an OGCM for the Indian Ocean

 Remotely sensed observations of sea-level anomaly and sea-surface temperature have been assimilated into an implementation of the Miami Isopycnic Coordinate Ocean Model (MICOM) for the Indian Ocean using the Ensemble Kalman Filter (EnKF). The system has been applied in a hindcast validation experiment to examine the properties of the assimilation scheme when used with a full ocean general circulation model and real observations. This work is considered as a first step towards an operational ocean monitoring and forecasting system for the Indian Ocean. The assimilation of real data has demonstrated that the sequential EnKF can efficiently control the model evolution in time. The use of data assimilation requires a significant amount of additional processing and computational resources. However, we have tried to justify the cost of using a sophisticated assimilation scheme by demonstrating strong regional and temporal dependencies of the covariance statistics, which include highly anisotropic and flow-dependent correlation functions. In particular, we observed a marked difference between error statistics in the equatorial region and at off-equatorial latitudes. We have also demonstrated how the assimilation of SLA and SST improves the model fields with respect to real observations. Independent in situ temperature profiles have been used to examine the impact of assimilating the remotely sensed observations. These intercomparisons have shown that the model temperature and salinity fields better resemble in situ observations in the assimilation experiment than in a model free-run case. On the other hand, it is also expected that assimilation of in situ profiles is needed to properly control the deep ocean circulation. Received: 8 January 2002 / Accepted: 8 April 2002     

A hybrid coordinate ocean model for shelf sea simulation

The general circulation in the North Sea and Skagerrak is simulated using the hybrid coordinate ocean model (HYCOM). Although HYCOM was originally developed for simulations of the open ocean, it has a design which should make it applicable also for coastal and shallow shelf seas. Thus, the objective of this study has been to examine the skills of the present version of HYCOM in a coastal shelf application, and to identify the areas where HYCOM needs to be further developed. To demonstrate the capability of the vertical coordinate in HYCOM, three experiments with different configurations of the vertical coordinate were carried out. In general, the results from these experiments compares quite well with in situ and satellite data, and the water masses and the general circulation in the North Sea and Skagerrak is reproduced in the simulations. Differences between the three experiments are small compared to other errors, which are related to a combined effect of model setup and properties of the vertical mixing scheme. Hence, it is difficult to quantify which vertical coordinate configuration works best for the coastal region. It is concluded that HYCOM can be used for simulations of coastal and shelf seas, and further suggestions for improving the model results are given. Since HYCOM also works well in open ocean and basin scale simulations, it may allow for a realistic modelling of the transition region between the open ocean and coastal shelf seas.     

A geochemical study of island-arc and back-arc tholeiites from the Scotia Sea

⁸⁷Sr/⁸⁶Sr and¹⁴³Nd/¹⁴⁴Nd ratios, REE and selected minor and trace elements are presented and compared for present-day volcanic rocks in the Scotia Sea.Tholeiitic basalts from the South Sandwich Islands show widely ranging contents of some lithophile elements, e.g. K₂O (0.09–0.55%) and Rb (1.55–14.2 ppm), but fairly constant Na₂O and Sr. Total REE contents range from about 4–20 times chondritic abundances with significant light-REE depletion and both positive and negative Eu anomalies. The variations in minor and trace element abundances are consistent with low-pressure fractional crystallization of plagioclase and clinopyroxene but only minor amounts of olivine. The⁸⁷Sr/⁸⁶Sr and¹⁴³Nd/¹⁴⁴Nd ratios of the parental magmas are thought be 0.7038–0.7039 and 0.51301–0.51314 respectively, and indicate derivation of at least some⁸⁷Sr from subducted ocean crust.The back-arc tholeiites in the Scotia Sea have lower⁸⁷Sr/⁸⁶Sr ratios (0.7028–0.7033), similar¹⁴³Nd/¹⁴⁴Nd ratios (0.51305) and are variably light-REE-enriched(Ce_N/Yb_N= 1.0–1.6). Total REE contents are comparable to those of the South Sandwich Islands tholeiites.     

Timing and thermochemical constraints on multi-element mineralisation at the Nori/RA Cu–Mo–U prospect, Great Bear magmatic zone, Northwest Territories, Canada

The timing of Cu–Mo–U mineralisation at the Nori/RA prospect in the Paleoproterozoic Great Bear magmatic zone has been investigated using Re–Os molybdenite and ⁴⁰Ar–³⁹Ar biotite geochronology. The Re–Os molybdenite ages presented are the first robust sulphide mineralisation ages derived from the Great Bear magmatic zone. Cu–Mo–U mineralisation is hosted in early to syn-deformational hydrothermal veins consisting of quartz and K-feldspar or more commonly tourmaline-biotite-quartz-K-feldspar, with associated wall-rock alteration assemblages being predominantly biotite. Sulphide and oxide minerals consist of chalcopyrite, molybdenite and uraninite with lesser pyrite and magnetite. Elevated light rare earth elements and tungsten concentrations associated with the Cu–Mo–U mineralisation have also been reported at the prospect by previous workers. Molybdenite and uraninite occur intimately in dravitic tourmaline growth zones and at grain margins, attesting to their syngenetic nature (with respect to hydrothermal veining). Two molybdenite separates yield Re–Os model ages of 1,874.4 ± 8.7 (2σ) and 1,872.4 ± 8.8 Ma (2σ) with a weighted average model age of 1,873.4 ± 6.1 Ma (2σ). Laser step heating of biotite from the marginal alteration of the wall-rock adjacent to the veins yields a ⁴⁰Ar–³⁹Ar maximum cooling age of 1,875 ± 8 Ma (MSWD = 3.8; 2σ), indistinguishable from the Re–Os molybdenite model age and a previously dated ‘syn-tectonic’ aplitic dyke in the region. Dravitic tourmaline hosts abundant primary liquid–vapour–solid-bearing fluid inclusions. Analytical results indicate liquid–vapour homogenisation at >260°C constraining the minimum temperature of mineralisation. The solids, which are possibly trapped, did not homogenise with the liquid–vapour by 400°C. Salinities in the inclusions are variable. Raman spectra identify that at least some of the solids are calcite and anhydrite. Raman spectra also confirm the vapour phases contain some CO₂; whereas clathrates or CH₄ was not observed or detected. Quartz grains only host secondary fluid inclusions, which fluoresce under ultraviolet light, indicating trapped hydrocarbons. We speculate that these resulted from Phanerozoic fluid circulation through the Proterozoic basement. The collective interpretation of the age, hydrothermal character and associated metals, high temperature and variable salinity suggests that the Nori/RA Cu–Mo–U mineralisation can be linked with the earliest stages of plutonism in the Great Bear magmatic zone. From a regional perspective, the mineralisation may pre-date the extensive multi-element mineralisation now recognised as part of the iron oxide copper–gold (IOCG) spectrum of deposits. As IOCG provinces generally contain a variety of mineralisation styles, we interpret this as the earliest phase of the extensive mineralising system.