Error evolution in the dynamics of an ocean general circulation model

The problem of error propagation is considered for spatially uncorrelated errors of the barotropic stream function in an oceanic general circulation model (OGCM). Such errors typically occur when altimetric data from satellites are assimilated into ocean models. It is shown that the error decays at first due to the dissipation of the smallest scales in the error field. The error then grows exponentially before it saturates at the value corresponding to the difference between independent realizations. A simple analytic formula for the error behavior is derived; it matches the numerical results documented for the present primitive-equation ocean model, and other models in the literature.     

The crystal and magnetic structure of Li-aegirine LiFe3+Si2O6: a temperature-dependent study

Single crystals of Li-aegirine LiFe³⁺Si₂O₆ were synthesized at 1573?K and 3?GPa, and a polycrystalline sample suitable for neutron diffraction was produced by ceramic sintering at 1223?K. LiFe³⁺Si₂O₆ is monoclinic, space group C2/c, a=9.6641(2)?Å, b= 8.6612(3)?Å, c=5.2924(2)?Å, β=110.12(1)^° at 300?K as refined from powder neutron data. At 229?K Li-aegirine undergoes a phase transition from C2/c to P2₁ /c. This is indicated by strong discontinuities in the temperature variation of the lattice parameters, especially for the monoclinic angle β and by the appearance of Bragg reflections (hkl) with h+k≠2n. In the low-temperature form two non-equivalent Si-sites with 〈SiA–O〉=1.622?Å and 〈SiB–O〉=1.624?Å at 100?K are present. The bridging angles of the SiO₄ tetrahedra O3–O3–O3 are 192.55(8)° and 160.02(9)° at 100?K in the two independent tetrahedral chains in space group P2₁ /c, whereas it is 180.83(9)° at 300?K in the high-temperature C2/c phase, i.e. the chains are nearly fully expanded. Upon the phase transition the Li-coordination changes from six to five. At 100?K four Li–O bond lengths lie within 2.072(4)–2.172(3)?Å, the fifth Li–O bond length is 2.356(4)?Å, whereas the Li–O3?A bond lengths amount to 2.796(4)?Å. From ⁵⁷Fe Mössbauer spectroscopic measurements between 80 and 500?K the structural phase transition is characterized by a small discontinuity of the quadrupole splitting. Temperature-dependent neutron powder diffraction experiments show first occurrence of magnetic reflections at 16.5?K in good agreement with the point of inflection in the temperature-dependent magnetization of LiFe³⁺Si₂O₆. Distinct preordering phenomena can be observed up to 35?K. At the magnetic phase transition the unit cell parameters exhibit a pronounced magneto-striction of the lattice. Below T _N Li-aegirine shows a collinear antiferromagnetic structure. From our neutron powder diffraction experiments we extract a collinear antiferromagnetic spin arrangement within the a–c plane.     

Compositional variation in Hg-bearing sphalerite from the polymetallic Eskay Creek deposit, British Columbia, Canada

The Eskay Creek, British Columbia, Canada, is a polymetallic, gold- and silver-rich, volcanic-hosted, massive sulfide deposit. The ore in the deposit is divided into subzones distinguished by mineralogy, texture, grade and metallurgical characteristics. This study presents the results of a mineralogical examination of three composite field samples, with emphasis on the chemistry of sphalerite. Sphalerite is associated with variable amounts of Hg-tetrahedrite and cinnabar, and an array of sulfides, sulphosalts and non-opaque minerals. Electron micro probe analyses of sphalerite in the three composite samples reveal wide variations in compositions. The Hg content in sphalerite in the three samples varies between 0.08 and 16.35 wt%, whereas the Fe content ranges from 0.33 to 2.29 wt%. The chemical formula of the sphalerite shows the compositional range (Zn_0.89–0.98Hg_0.01–0.09.Fe_0.005–0.02)S. Sphalerite exhibits an almost perfect substitution of Hg and Zn, as shown by the negative covariance between them. Sphalerite with the highest Hg contents tends to have the lowest Fe concentrations. The highest Hg contents in sphalerite are recorded in the samples with the highest bulk Hg concentrations and with the highest cinnabar contents.

The compositional variations of sphalerite are important because they can be used in mapping ore forming fluids and indicate possible temporal variations. Second, determination of the compositional variation of the sphalerite in the mine has metallurgical implications because the mineral is an important Zn source. The mineralogical data indicate that non-physical processes (e.g. pyrometallurgy) must be used to separate Hg from Zn concentrates, with direct environmental implications, that is, release of metals, such as Hg, into the environment during mining and processing.     

Origin of ore fluids in the Muruntau gold system: Constraints from noble gas, carbon isotope and halogen data

Hydrothermal vein minerals directly associated with native gold mineralization in the Muruntau vein system (Uzbekistan) have been studied for noble gas, carbon isotope and halogen chemistry of the trapped ore-related fluids. Helium trapped in early arsenopyrite 1, which has preserved the original fluid signature better than associated scheelite and quartz, indicates a small input from a mantle source (?5% of total He). However, the overwhelming majority of the He in the fluid (∼95%) is from crustal sources. The noble gases Ne, Kr and Xe in the sample fluids are dominated by gases of atmospheric origin. The carbon isotope (δ¹³C: −2.1‰ to −5.3‰) and halogen characteristics of the fluids (log Br/Cl: −2.64 to −3.23) lend further support for the activity of juvenile fluids during the main ore stage. The high proportion of crustal components in the ore-forming fluids may be explained by intense fluid-rock interaction and is also supported by previous Nd and Sr isotope studies. The involvement of a juvenile fluid component during the main stage of hydrothermal activity at Muruntau (∼275 Ma) can be linked to the emplacement of lamprophyric dikes at Muruntau, due to apparently overlapping ages for high-temperature alteration, related ore vein formation and intrusion of the dikes. The input of mantle-derived fluids, possibly related to the Hercynian collisional event in the western Tien Shan, stimulated intense fluid-rock interaction in the crust. In this context, the mantle-derived fluid should be considered as one possible carrier of metals. Significant amounts of external meteoric fluids circulating in fracture systems are interpreted to have modified the noble gas signature of fluid in quartz, mostly during late, low temperature fluid circulation.     

Elastic and visco-elastic stress triggering in the South Iceland Seismic Zone due to large earthquakes since 1706

Damaging earthquakes in the South Iceland Seismic Zone (SISZ) occur fairly regularly and often as a series of events with a few days only between individual events. Tolerably reliable information on epicentre locations and mechanisms are available for 13 M ≥ 6 events between 1706 and 2000. For these events, we computed the co- and post-seismic stress fields, hereby approximating the SISZ by a mixed elastic/visco-elastic layered half-space. The horizontal shear stress and the Coulomb stress changes were analysed to detect possible trigger mechanisms, which may aid future earthquake mitigation efforts. We tested several criteria but must conclude that the start of an earthquake series in the SISZ cannot be explained by triggering through previous events. Inside an individual series, however, one may infer triggering. Our results are in contradiction with the findings in other regions of the world. The reason might be related to the fact that the SISZ is not a mature fault zone, in which old faults are re-activated if a certain stress level threshold is passed. In addition, uncertainties in the model parameters as well as the neglect of horizontal variations in the model and of possible stress transfer due to volcanic activity further complicate the evaluation of our results and need to be taken into account in future studies.     

Preliminary reconstruction of sediment-source linkages for the past 6000 yr at the Petit Lac d'Annecy, France, based on mineral magnetic data

An 8 m core from the central plain of the Petit Lac d'Annecy, France, two floodplain cores, river bedload sediments and several hundred soil samples from the catchment have been studied using magnetic techniques. The soils, mainly developed on limestones and local glacial tills, show widespread magnetic enhancement with higher ferrimagnetic concentrations and contents of SP grains than found in the lake sediments. Some soils show significant concentrations of canted antiferromagnetic minerals (mainly haematite). Using magnetic quotient parameters the surface soils are classified into four mineralogical types. The lake and floodplain sediment properties over the past 6000 yrs can largely be explained by the erosion and deposition of these sources, with a smaller superimposed biogenic (magnetosomes) signal. Derived sediment-source linkages allow the construction of several hypotheses about geomorphological changes in the catchment system: (i) the long-term erosion of high altitude unweathered substrates has gradually increased towards the present day; (ii) the erosion of high altitude soils has increased within the last 1000 yrs, possibly during the period of the 'Little Ice Age'; (iii) shifts towards an increased erosion of surface lowland soil occurred ~2000 and 1000 yrs ago and may be linked to an accelerated accretion of floodplain overbank deposits; (iv) there has been a significant storage of surface soil within floodplains, which leads to an underestimation of the importance of soil erosion in the lake sediment records; (v) the sediment transported by high magnitude, low frequency flood events has shifted in source from high altitude soils before ~1000 cal. yr BP to lowland and mid-altitude free draining soils after ~1000 cal. yr BP.     

High-resolution sea ice in long-term global ocean GCM integrations

The resolution of the sea-ice component of a coarse-resolution global ocean general circulation model (GCM) has been enhanced to about 22 km in the Southern Ocean. The ocean GCM is designed for long-term integrations suitable for investigations of the deep-ocean equilibrium response to changes in southern hemisphere high-latitude processes. The space and time scales of the high-resolution sea-ice component are commensurate with those of the resolution of satellite passive-microwave sea-ice data. This provides the opportunity for a rigorous evaluation of simulated sea-ice characteristics. It is found that the satellite-derived continuous high ice concentration of the interior winter ice pack can only be captured when vertical oceanic mixing is modified in a way that less local, intermittent convection occurs. Furthermore, the width and the variability of the coastal polynyas around the Antarctic continent and its ice shelves are best captured when some form of ice-shelf melting is accounted for. The width of the wintertime ice edge is reasonably reproduced, while its variability remains underestimated, closely following the coarse-grid pattern of the ocean model due to its high dependence on ocean temperature. Additional variability besides daily winds, e.g. in form of idealized tidal currents, improves the temporal and spatial ice-edge variability, while leads in the interior ice pack become more abundant, more in line with the fine-scale satellite-derived texture. The coast- or ice-shelf line is described on the fine grid based on satellite passive-microwave data. This method requires parts of a coarse coastal ocean grid cell to be covered by an inert layer of “fast ice” or “ice shelf”. Reasonable long-term global deep-ocean properties can only be achieved when these areas are not inert, i.e. are exposed to heat flux and ice growth, or when the vertical mixing parameterization allows for excessive open-ocean convection. The model area exposed to cold high-latitude atmospheric conditions thus being most decisive for a realistic representation of the long-term deep-ocean properties, suggests that high-latitude coastlines are definitely in need of being represented at high resolution, including ice sheets and their effects on the heat and freshwater flux for the ocean.     

Evolution of a Holsteinian (MIS 11c) palaeolake based on a 12‐ka‐long diatom record from Dethlingen (northern Germany)

To provide insights into the long‐term evolution of aquatic ecosystems without human interference, we here evaluate a decadal‐ to centennial‐scale‐resolution diatom record spanning about 12 ka of the Holsteinian interglacial (Marine Isotope Stage 11c). Using a partially varved sediment core from the Dethlingen palaeolake (northern Germany), which has previously been studied for palynological and microfacies signals, we document the co‐evolution of the aquatic and surrounding terrestrial environment. The diatom record is dominated by the genera Stephanodiscus, Aulacoseira, Ulnaria and Fragilaria. Based on the diatom assemblages and physical sediment properties, the evolution of the Dethlingen palaeolake can be subdivided into three major phases. During the oldest phase (lasting ～1900 varve years), the lake was ～10–15 m deep and characterized by anoxic bottom‐water conditions and a high nutrient content. The following ～5600 years exhibited water depths >20 m, maximum diatom and Pediastrum productivity, and a peak in allochtonous nutrient input. During this phase, water‐column mixing became more vigorous, resulting in a breakdown of anoxia. The youngest lake phase (～4000–5000 years) was characterized by decreasing water depth, turbulent water conditions and decreased nutrient loading. Based on our palaeolimnological data, we conclude that the evolution of the Dethlingen palaeolake during the Holsteinian interglacial responded closely to (i) changes within the catchment area (as documented by vegetation and sedimentation) related to the transition from closed forests growing on nutrient‐rich soils (mesocratic forest phase) to open forests developing on poor soils (oligocratic forest phase), and (ii) short‐term climate variability as reflected in centennial‐scale climate perturbations.