On mountain wave drag over complex terrain

Summary Mountain wave drag is calculated for rotating, stratified, nonhydrostatic Boussinesq flow over a mountain ridge using linear theory for a variety of mountain profiles representing complex/irregular terrain. The inclusion of a sinusoidal corrugation to the familiar witch-of-Agnesi profile creates a stegosaurus profile. The associated drag is greatly enhanced for mesoscale mountains when the corrugation wave-number matches that for the dominant inertia-gravity wave contribution to the cross-mountain surface pressure gradient. Similarly, increasing the jaggedness (by decreasing the exponentb) increases the drag for mesoscale mountains whose topographic spectral intensity,M(k), has the form of a power law:M(k)=mk ^–b wherek is the zonal wavenumber.Spectral analysis of one-kilometer resolution topographic data for the Appalachian Mountains suggests that a power law profile withb=1.7 accurately represents the topographic spectral intensity and that it yields good estimates of the drag.The application of these results to the parameterization of mountain wave drag in general circulation models is discussed.With 7 Figures     

Transfer and early diagenesis of biogenic silica oxygen isotope signals during settling and sedimentation of diatoms in a temperate freshwater lake (Lake Holzmaar, Germany)

We have investigated the transfer of oxygen isotope signals of diatomaceous silica (δ¹⁸O_diatom) from the epilimnion (0-7 m) through the hypolimnion to the lake bottom (∼20 m) in freshwater Lake Holzmaar, Germany. Sediment-traps were deployed in 2001 at depths of 7 and 16 m to harvest fresh diatoms every 28 days. The 7 m trap collected diatoms from the epilimnion being the main zone of primary production, while the 16 m trap collected material already settled through the hypolimnion. Also a bottom sediment sample was taken containing diatom frustules from approximately the last 25 years. The δ¹⁸O_diatom values of the 7 m trap varied from 29.4‰ in spring/autumn to 26.2‰ in summer according to the temperature dependence of oxygen isotope fractionation and represent the initial isotope signal in this study. Remarkably, despite the short settling distance δ¹⁸O_diatom values of the 7 and the 16 m trap were identical only during spring and autumn seasons while from April to September δ¹⁸O_diatom values of the 16 m trap were roughly ∼1.5‰ enriched in ¹⁸O compared to those of the 7 m trap. Isotopic exchange with the isotopically lighter water of the hypolimnion would shift the δ¹⁸O_diatom value to lower values during settling from 7 to 16 m excluding this process as a cause for the deviation. Dissolution of opal during settling with intact organic coatings of the diatom cells and near neutral pH of the water should only cause a minor enrichment of the 16 m values. Nevertheless, opal from the bottom sediment was found to be 2.5‰ enriched in ¹⁸O compared to the weighted average of the opal from the 7 m trap. Thus, resuspension of bottom material must have contributed to the intermediate δ¹⁸O_diatom signal of the 16 m trap during summer. Dissolution experiments allowed further investigation of the cause for the remarkably enriched δ¹⁸O_diatom value of the bottom sediment. Experiments with different fresh diatomaceous materials show an increase of opaline ¹⁸O at high pH values which is remarkably reduced when organic coatings of the cells still exist or at near neutral pH. In contrast, high pH conditions do not affect the δ¹⁸O_diatom values of sub-fossil and even fossil opal. IR analyses show that the ¹⁸O enrichment of the sedimentary silica is associated with a decrease in Si-OH groups and the formation of Si-O-Si linkages. This indicates a silica dehydroxylation process as cause for the isotopic enrichment of the bottom sediment. Silica dissolution and dehydroxylation clearly induce a maturation process of the diatom oxygen isotope signal presumably following an exponential behaviour with a rapid initial phase of signal alteration. The dynamics of this process is of particular importance for the quantitative interpretation of sedimentary δ¹⁸O_diatom values in terms of palaeothermometry.     

Geochemistry and Rb–Sr geochronology of a ductile shear zone in the Orlica-Śnieżnik dome (West Sudetes, Poland)

Amphibolite-facies orthogneisses of the Orlica-vnie™nik dome in the West Sudetes (Poland) show a local continuous transition from weakly deformed augen gneisses to finely laminated mylonites. Field evidence indicates that ductile shearing developed pre- or syntectonically to a migmatization event. Bulk-rock compositions of variably deformed samples yield no indications for deformation- and/or fluid-enhanced element mobility and redistribution. ⁸⁷Rb-⁸⁶Sr geochronology (biotite, phengite, whole rock) places time constraints on the deformation process and the post-orogenic cooling history. Phengite- and biotite-whole-rock pairs yield Rb-Sr ages of 340 to 334 Ma and 335 to 294 Ma, respectively, independent of the degree of deformation. The weighted mean of phengite-whole-rock pairs indicates an age of 337.4DŽ.3 Ma. Combining most of the biotite-whole-rock data yields a weighted mean age of 328.6dž.4 Ma. Because of their different closure temperatures for the Rb-Sr system, these differences are interpreted to date cooling after a thermal event. Direct dating of the deformation is not possible, but the cooling history record defines a minimum age for the development of ductile shearing and the last migmatization event. These time constraints provide evidence for the initiation of crustal collapse during or immediately following peak metamorphic conditions. The results of this study further document the importance of Variscan metamorphism in the Orlica-vnie™nik dome.     

3D shear-wave velocity structure of the Eifel plume, Germany

The Quaternary Eifel volcanic fields, situated on the Rhenish Massif in Germany, are the focus of a major interdisciplinary project. The aim is a detailed study of the crustal and mantle structure of the intraplate volcanic fields and their deep origin. Recent results from a teleseismic P-wave tomography study reveal a deep low-velocity structure which we infer to be a plume in the upper mantle underneath the volcanic area [J.R.R. Ritter et al., Earth Planet. Sci. Lett. 186 (2001) 7-14]. Here we present a travel-time investigation of 5038 teleseismic shear-wave arrivals in the same region. First, the transverse (T) and radial (R) component travel-time residuals are treated separately to identify possible effects of seismic anisotropy. A comparison of 2044 T- and 2994 R-component residuals demonstrates that anisotropy does not cause any first-order travel-time effects. The data sets reveal a deep-seated low-velocity anomaly beneath the volcanic region, causing a delay for teleseismic shear waves of about 3 s. Using 3773 combined R- and T-component residuals, an isotropic non-linear inversion is calculated. The tomographic images reveal a prominent S-wave velocity reduction in the upper mantle underneath the Eifel region. The anomaly extends down to at least 400 km depth. The velocity contrast to the surrounding mantle is depth-dependent (from −5% at 31-100 km depth to at least −1% at 400 km depth). At about 170-240 km depth the anomaly is nearly absent. The resolution of the data is sufficient to recover the described features, however the anomaly in the lower asthenosphere is underestimated due to smearing and damping. The main anomaly is similar to the P-wave model except the latter lacks the ‘hole’ near 200 km depth, and both are consistent with an upper mantle plume structure. For plausible anhydrous plume material in the uppermost 100 km of the mantle, an excess temperature as great as 200-300 K is estimated from the seismic anomaly. However, 1% partial melt reduces the required temperature anomaly to about 100 K. The temperature anomaly associated with the deeper part of the plume (250 to about 450 km depth) is at least 70 K. However, this estimate is quite uncertain, because the amplitude of the shear-wave anomaly may be larger than the modelled one. Another possibility is water in the upwelling material. The gap at 170-240 km depth could arise from an increase of the shear modulus caused by dehydration processes which would not affect P-wave velocities as much. An interaction of temperature and compositional variations, including melt and possibly water, makes it difficult to differentiate quantitatively between the causes of the deep-seated low-velocity anomaly.     

Experimental determination of the solubility of Au in silicate melts

Summary We report here new data on the solubility of Au in silicate melts of anorthite-diopside eutectic composition at a wide range of oxygen fugacities, from pure oxygen to 10^–8 atm, and at a temperature range of 1300 °C to 1480 °C. Because experiments were done with metal loops at temperatures above the Au-melting temperature, PdAu-metal-alloys had to be used. Pd-solubility data derived from the same set of experiments agree with earlier data obtained from experiments with pure Pd-metal (Borisov et al., 1994a). The results of the present experiments show that Pd-solubilities are by a factor of 2 to 6 higher than Au-solubilities. Both, Au and Pd solubilities decrease with decreasing oxygen fugacity. At oxygen fugacities below the iron-wiistite buffer (IW) Au solubility increases with decreasing fO₂ probably reflecting formation of Au-silicides at such reducing conditions. Compared to Pd, Au has higher activity coefficients in Fe-metal and lower solubility in silicate melts. This leads to similar metal-silicate partition coefficients for both elements. At a temperature of 1350 °C and an oxygen fugacity corresponding to IW-2 D^Au (met/sil) is about 2.5 · 10⁷ and D^Pd (met/sil) about 1.6 · 10⁷. Thus similar behavior is expected during metal separation in planetary bodies including core formation in the Earth. The metal/silicate partition coefficient of Ir is, however, by several orders of magnitudes higher (Borisov and Palme, 1995a). Equilibration with chondritic metal will therefore lead to grossly non-chondritic Pd/Ir or Au/Ir ratios in coexisting silicate phases. Chondritic ratios are thus indicative of the presence of unfractionated meteoritic components. Samples from the upper mantle of the Earth, for example, reflect the admixture of a late unfractionated (chondritic) veneer (e.g.,Kimura et al., 1974;Jagoutz et al., 1979).Solubilities of Pd and Au in silicate melts are much higher than the contents in terrestrial basalts implying that the abundances of these two elements are not buffered by residual PGE- and Au-containing alloys. The most likely process for fractionating PGEs in terrestrial magmas are mineral-melt (e.g., olivine/melt) equilibria.

---

Experimentelle Bestimmung der Löslichkeit von Au in Silikatschmelzen

Zusammenfassung In der vorliegenden Arbeit wird über die Ergebnisse der Bestimmung der Löslichkeit von Au in Silikatschmelzen mit der Zusammensetzung des Anorthit-Diopsid Eutektikums berichtet. Die Versuche wurden mittels Metallschlaufe über einen weiten Sauerstoffpartialdruckbereich, von reinem Sauerstoff bis zu 10^–8 atm und in einem Temperaturbereich von 1300 °C bis 1480 °C, durchgeführt. Da diese Temperaturen jedoch den Au-Schmelzpunkt überschreiten, wurde mit AuPd-Legierungen gearbeitet. Die Ergebnisse der dadurch zusätzlich erhaltenen Pd-Versuche stimmen mit früher bestimmten, mit reinen Pd-Schlaufen durchgeführten Pd-Löslichkkeiten überein (Borisov et al., 1994a). Die auf reine Metalle zurückgerechneten Löslichkeiten von Pd sind um einen Faktor 2 bis 6 mal höher als die entsprechenden Au-Löslichkeiten. Die Löslichkeiten beider Metalle nehmen mit abnehmendem Sauerstoffpartialdruck ab. Unter noch stärker reduzierenden Bedingungen (Eisen-Wüstit Gleichgewicht) nimmt die Löslichkeit von Au jedoch zu. Dies könnte auf die Bildung von Au-Siliziden zurückzuführen sein.Im Vergleich zu Pd sind die Aktivitätskoeffizienten von Au in metallischem Eisen höher, die Löslichkeiten in Silikatschmelzen jedoch niedriger. Das führt zu ähnlichen Metall/Silikat Verteilungskoeffizienten von Au und Pd. Bei einer Temperatur von 1350 °C und einer Sauerstoffugazität von IW-2 ergeben sich für D^Au (met/sil) 2.5 · 10⁷ und für D^Pd (met/sil) 1.6* 10⁷. Der Metall/Silikat-Verteilungskoeffizient von Ir ist jedoch unter den gleichen Bedingungen um mehrere Größenordnungen höher (Borisov andPalme, 1995a). Ein chondritisches Pd/Ir- oder Au/Ir-Verhältnis kann also auf die Anwesenheit einer unfraktionierten chondritischen Komponente zurückgeführt werden. Dies gilt beispielsweise für Proben aus dem oberen Erdmantel. Hier handelt es sich vermutlich um Zumischung einer späten chondritischen Akkretionskomponente, die sich nicht mehr mit einer metallischen Phase (Kern) ins Gleichgewicht gesetzt hat (z.B.Kimura et al., 1974,Jagoutz et al., 1979).Die Löslichkeiten von Pd und Au in Silikatschmelzen sind wesentlich höher als ihre Gehalte in basaltischen und komatiitischen Laven. Dies bedeutet, daß Au und Pd in Schmelzen aus dem Erdmantel nicht durch residuale Au- und/oder Pd-haltige Metall phasen bestimmt sind. Gleichgewichte zwischen Schmelze und Mineralen (z.B. Olivin) sind die wahrscheinlichsten Fraktionierungsmechanismen für Platingruppenelemente in terrestrischen Magmen.

---

---

With 5 Figures