Fully coupled generalised hybrid-mixed finite element approximation of two-phase two-component flow in porous media. Part II: numerical scheme and numerical results

We consider the modeling and simulation of compositional two-phase flow in a porous medium, where one phase is allowed to vanish or appear. The modeling of Marchand et al. (in review) leads to a nonlinear system of two conservation equations. Each conservation equation contains several nonlinear diffusion terms, which in general cannot be written as a function of the gradients of the two principal unknowns. Also the diffusion coefficients are not necessarily explicit local functions of them. For the generalised mixed finite elements approximation, Lagrange multipliers associated to each principal unknown are introduced, the sum of the diffusive fluxes of each component is explicitly eliminated and the static condensation leads to a “global” nonlinear system of equations only in the Lagrange multipliers also including complementarity conditions to cope with vanishing or appearing phases. After time discretisation, this system can be solved at each time step using a semi-smooth Newton method. The static condensation involves “local” nonlinear systems of equations associated to each element, solved also by a semismooth Newton method. The algorithm is successfully applied to 1D and 2D examples of water–hydrogen flow involving gas phase appearance and disappearance.     

Topic 3: State of the art of mathematical modelling of agricultural impact on groundwater pollution

The state of the art of mathematical modelling of groundwater pollution by agricultural activities is well advanced at the present time and a large number of models are available for a wide variety of water pollution problems. Also a variety of computational methods are available from very simple ones to the most advanced. In general, there are two ways of obtaining a water quality model: “from below,” when water quality data based on actual observation or measurement are explicitly analyzed by the calculus; and “from above,” when the calculus is available in advance as a theory and the user tries to find its interpretation for a given case and from rather limited data. Classifications may be based on the type of modelled processes, the mathematical approach, or on the purpose of the models. Eleven papers were presented in the Symposium on this topic, with most of the models being of the dynamic type.     

Geochemical studies on Dachang antimony ore deposit in Qinglong, Guizhou Province

The Dachang antimony deposit in Qinglong,Guizhou Province,is strictly controlled by the “Dachang Layer” which is a complex altered rock occurring at unconformity between the Permian Emeishan basalt and the Maokou limestone.Based on the studies of the hanging-and foot-wall rocks,the trace elements and REE contents of the rocks and ores and heavy placer minerals in the basalt,this paper is focused on the relations between these data and the “Dachang Layer”and its hanging- and oot-wall rocks.The author pointed out that the “Dachang Layer” and basalt are the source-beds of antimony;ilmenite and magnetite are the major mineral carriers of antimony.In the processes of halmyrosis and burial metamorphism of the “Dachang Layer” an basalt,antimony was mobilized along with the mobilization of iron and was preliminarily concentrated in the“ Dachang Layer”.     

Plagioclase residence times at two island arc volcanoes (Kameni Islands, Santorini, and Soufriere, St. Vincent) determined by Sr diffusion systematics

The diffusive relaxation of trace element profiles in plagioclase phenocrysts may provide important constraints on magma residence times in crustal magma chambers. Initial trace element profiles in plagioclase phenocrysts are governed by variations in the concentration of a trace element in the melt and by the plagioclase-melt partition coefficient. Trace element diffusion will subsequently act to modify this initial profile and – given enough time – produce a profile that is in equilibrium with the anorthite variations within the crystal. We argue that the trace element partition coefficient D^a/b between two parts a and b of a plagioclase crystal of variable anorthite content is equal to the ratio of their crystal-liquid partition coefficients, and that the equilibrium profile of the crystal can be calculated. The time required to establish diffusive equilibrium is dependent on the wavelength and amplitude of the initial trace element concentration range and on the diffusivity of the trace element in plagioclase. Strontium plagioclase-melt partition coefficients and diffusivities are calculated for a range of magmatic temperatures and plagioclase compositions. A one-dimensional diffusion model is developed that describes the diffusive destruction of oscillatory trace element zoning with time and allows the calculation of upper limits for plagioclase crystal residence times in a magma reservoir. The model is tested using major and trace element concentrations measured along crystal traverses of plagioclase phenocrysts from the Kameni Island dacites, Santorini, and from the 1979 Soufriere andesite, St. Vincent. Three out of eight plagioclase phenocrysts have Sr concentration profiles that are not in diffusive equilibrium. For these, the diffusion model is employed to calculate maximum crystal residence times from incomplete diffusive equilibration of trace element zoning in plagioclase. Maximum crystal residence times range from 100 to 450 years. This is in good agreement with estimates from crystal size distribution and from Ra-Th disequilibrium studies for the Kameni Islands. For Soufriere, however, such short residence times are incompatible with U-Th mineral errorchron data that suggest residence times of >40 ka in a thermally buffered magma reservoir. To reconcile these apparently different ages, we invoke a more complicated magmatic history for Soufriere where an initially buffered magma reservoir is disturbed by magma mixing and suffers limited additional crystal fractionation prior to eruption. Received: 20 October 1998 / Accepted: 22 March 1999     

Redox environment effect on redox sensitive elements in surface sediments of the Changjiang Estuary hypoxia zone

The grain size and element (including redox sensitive elements and terrigenous elements) concentration of surface sediments from the Changjiang Estuary hypoxia zone and its adjacent sea area were measured in this research. Based on the obtained data, the hypoxic environment’s influence on the distribution of elements in surface sediments was further studied. We believe that the “redox environment effect” greatly influences the distribution of the RSE, which reveals the “patchy enrichment pattern” offshore in the hypoxia zone, while the distribution of the terrigenous elements which shows the “stripped enrichment pattern” near shore is mainly affected by “granularity effects”. Due to the existence of the hypoxia zone of the Changjiang Estuary, the distribution of the RSE such as Mo, Cd and V in the study area exhibits the characteristics of “redox environment effects”. __________ Translated from Marine Geology & Quaternary Geology, 2007, 27(3):1–8 [译自: 海洋地质与第四纪地质]     

Chemodynamics of trace metal fractions in surface sediments of the Pandoh Lake,Lesser Himalaya,India

The seasonal variation in the trace metals’ concentrations (Cd, Co, Cu, Fe, Mn, Ni, Pb, and Zn) were investigated in surface sediments of the Pandoh Lake. The horizontal distribution of TC, TN, and TP reflects spatial and temporal differences in sedimentary organic production. The chemical sequential extraction of heavy metals was carried out by seven-step fractionation scheme (Leleyter and Probst in Int J Environ Chem 73:109–128, 1999). The significant concentrations of Ni and Cd were associated with “water soluble (Eua)” fraction in the monsoon and winter, respectively, while “exchangeable (Exch)” and “carbonate-bound (Carb)” fractions for Ni and Cd were abundant in winter and summer. The Cd, Cu, and Pb associated with “Exch” fraction in the summer season support their availability on exchange sites due to oxidized nature of surface sediments. Enrichment of Co, Fe, Mn, and Zn in “AFeO” fraction showed poor bioavailability, while Cd, Cu, and Mn in the monsoon, Co in the winter and summer, and Zn in the winter season showed significant “organically bound (Org)” fraction. The ANOVA was significant for chemical fractions of trace elements except “Carb” fraction of Pb and Zn and “CFeO” fraction of Pb. Factor analysis revealed that the “Eua”, “Exch”, and “Carb” fractions together control the metal enrichment of “MnO”, “AFeO”, and “CFeO” fractions in the summer season.     

Age constraints on faulting and fault reactivation: a multi-chronological approach

Movement within the Earth’s upper crust is commonly accommodated by faults or shear zones, ranging in scale from micro-displacements to regional tectonic lineaments. Since faults are active on different time scales and can be repeatedly reactivated, their displacement chronology is difficult to reconstruct. This study represents a multi-geochronological approach to unravel the evolution of an intracontinental fault zone locality along the Danube Fault, central Europe. At the investigated fault locality, ancient motion has produced a cataclastic deformation zone in which the cataclastic material was subjected to hydrothermal alteration and K-feldspar was almost completely replaced by illite and other phyllosilicates. Five different geochronological techniques (zircon Pb-evaporation, K–Ar and Rb–Sr illite, apatite fission track and fluorite (U-Th)/He) have been applied to explore the temporal fault activity. The upper time limit for initiation of faulting is constrained by the crystallization age of the primary rock type (known as “Kristallgranit”) at 325 ± 7 Ma, whereas the K–Ar and Rb–Sr ages of two illite fractions <2 μm (266–255 Ma) are interpreted to date fluid infiltration events during the final stage of the cataclastic deformation period. During this time, the “Kristallgranit” was already at or near the Earth’s surface as indicated by the sedimentary record and thermal modelling results of apatite fission track data. (U–Th)/He thermochronology of two single fluorite grains from a fluorite–quartz vein within the fault zone yield Cretaceous ages that clearly postdate their Late-Variscan mineralization age. We propose that later reactivation of the fault caused loss of helium in the fluorites. This assertion is supported by geological evidence, i.e. offsets of Jurassic and Cretaceous sediments along the fault and apatite fission track thermal modelling results are consistent with the prevalence of elevated temperatures (50–80°C) in the fault zone during the Cretaceous.     

Landslide Hazard Zonation using Remote Sensing and GIS: a case study of Dikrong river basin, Arunachal Pradesh, India

Landslides are among the most costly and damaging natural hazards in mountainous regions, triggered mainly under the influence of earthquakes and/or rainfall. In the present study, Landslide Hazard Zonation (LHZ) of Dikrong river basin of Arunachal Pradesh was carried out using Remote Sensing and Geographic Information System (GIS). Various thematic layers namely slope, photo-lineament buffer, thrust buffer, relative relief map, geology and land use / land cover map were generated using remote sensing data and GIS. The weighting-rating system based on the relative importance of various causative factors as derived from remotely sensed data and other thematic maps were used for the LHZ. The different classes of thematic layers were assigned the corresponding rating value as attribute information in the GIS and an “attribute map” was generated for each data layer. Each class within a thematic layer was assigned an ordinal rating from 0 to 9. Summation of these attribute maps were then multiplied by the corresponding weights to yield the Landslide Hazard Index (LHI) for each cell. Using trial and error method the weight-rating values have been re-adjusted. The LHI threshold values used were: 142, 165, 189 and 216. A LHZ map was prepared showing the five zones, namely “very low hazard”, “low hazard”, “moderate hazard”, “high hazard” and “very high hazard” by using the “slicing” operation.     

Modelling the structure of the metastable phases in the reaction sequence kaolinite-mullite by X-ray scattering experiments

The crystal structure of two kaolinites with a different degree of planar disorder and their high temperature (HT) decomposition products have been investigated by means of symmetrical transmission powder X-ray diffractometry and the subsequent calculation of the radial distribution function (RDF). An ordered and a disordered kaolinite (international standard Kga-1 and Kga-2) were selected along with their HT reaction products at 600°C (metakaolinite) and at 900°C (mullite precursor). The calcinations and quenches have been carried out “ex-situ” prior to the data collection. The procedure for collection and reduction of the data and their interpretation is outlined. The structural characterization of the HT phase transformations is obtained by qualitative analysis of the RDF curves and is related to the extent of planar disorder in the original material. A model for the structure of metakaolinite is proposed on the basis of the present experiments, and of a review of literature data. According to this model, the structure of metakaolinite displays an almost uncollapsed structure with a short-range order in which Al is mainly in four-fold coordination. An “ab initio” modelling of the RDF curves provides a quantitative evaluation of the proposed structural model. Differences in the HT products from different kaolinites are directly correlated with the density of stacking faults in the starting kaolinites. Moreover the structure of the precursor of mullite is described and discussed with respect to the literature data. The nature and role of γ-alumina is discussed as well. Finally the DTA exothermic peak of kaolinite is interpreted in the light of the proposed general model. Received: 18 March 1997 / Revised, accepted: 26 October 1997     

Dualistic distribution coefficients of elements in the system mineral-hydrothermal solution. I. Gold accumulation in pyrite

The use of trace elements (TE) as geochemical indicators is complicated by the dualism of their distribution coefficients D due to the additional (i.e., above the concentrations of an isomorphic component) incorporation of elements at structural defects of various nature (including the surface of the crystal). A pressing problem in this situation is to determine the true D values that pertain to the structural component of an admixture D ^str and evaluate effects of other modes of TE occurrence. Only upon distinguishing D ^str in the bulk coefficient D ^bulk it is possible to evaluate the ore potential of fluid in terms of certain TE from the composition of a mineral containing the TE. Pyrite synthesized in solutions of variable pH at 450°C and 1 kbar (100 MPa) at fluid portions sampled in a trap is utilized to demonstrate the role of a surface nonautonomous phase (NP) in the incorporation of gold in this mineral. The distribution coefficient of gold between pyrite and hydrothermal solution is 0.14 for “pure” pyrite and 0.05 for As-bearing pyrite (containing 0.02–0.05 wt % As), and these coefficients for NP are 310 and 170, respectively. This increases the D ^bulk for evenly distributed (“invisible”) gold by factors of four and nine. In contrast to the results of earlier studies conducted at room temperature and pressure or parameters close to them, our data demonstrate that the accumulation of “invisible” Au in pyrite is controlled not only by reducing adsorption with the development of Au(0) particles and films but also by Au incorporation in NP developing in the surface layer of the crystal approximately 500 nm thick as chemically bound Au [most likely as Au(I)]. The possible reason for the high absorption capacity of NP is the defect (pyrrhotite-like) structure, which is not saturated with bonds of excess S and sulfoxi onions.     

Preservation of oxygen isotope compositions in granulites from Northwestern Canada and Enderby Land,Antarctica: implications for high-temperature isotopic thermometry

 Preservation of high-temperature mineral isotopic compositions is necessary for successful high-temperature isotopic thermometry. Other requirements include large fractionations between constituent minerals, well-calibrated equilibria, carefully designed sampling strategies and data handling techniques that quantitatively account for retrograde exchange. Here, we apply isotopic thermometry and data handling techniques to calculate and contrast mineral-pair apparent temperature data and observed closure temperature data (T _c-observed) (cf. Farquhar et al. 1993) for the very high temperature (>900°C), dry granulites of the Taltson Magmatic Zone of Northwestern Canada and the Napier Complex of Enderby Land, Antarctica. The isotopic compositions of garnet grains from both terrains reflect high temperature conditions (>950°C) and point to this mineral as an excellent candidate for isotopic thermometry. The isotopic compositions of quartz, pyroxene, ilmenite and magnetite indicate that they equilibrated to lower temperature conditions (<900°C) due to faster rates of oxygen diffusion in these minerals, possibly enhanced by exsolution and ductile deformation, compared with garnet. Our temperature data for garnet and pyroxene are ≈200°C higher than is possible to explain by existing “wet” diffusion data, but is consistent with “dry” diffusion data, suggesting that the extremely dry nature of these rocks may have played a significant role in the preservation of high-temperature isotopic compositions. Both quartz and magnetite exhibit subgrain features, indicative of ductile deformation. Quartz-magnetite temperatures from the Napier complex are similar to those inferred for a late (D3) deformation and are lower than those predicted by “dry” diffusion data. We infer that the quartz-magnetite isotopic fractionations reflect deformation-enhanced exchange that accompanied D3. Garnet in these same samples did not undergo ductile deformation and did not exchange oxygen with coexisiting phases during cooling. This may reflect strain partitioning between less easily deformed, low abundance garnet and more easily deformed matrix quartz and magnetite. The resistance of garnet to ductile deformation in these rocks is a second reason why garnet is suitable for isotopic thermometry. Received: 6 February 1996 / Accepted: 25 April 1996     

Mantle Pb paradoxes: the sulfide solution

There is growing evidence that the budget of Pb in mantle peridotites is largely contained in sulfide, and that Pb partitions strongly into sulfide relative to silicate melt. In addition, there is evidence to suggest that diffusion rates of Pb in sulfide (solid or melt) are very fast. Given the possibility that sulfide melt “wets” sub-solidus mantle silicates, and has very low viscosity, the implications for Pb behavior during mantle melting are profound. There is only sparse experimental data relating to Pb partitioning between sulfide and silicate, and no data on Pb diffusion rates in sulfides. A full understanding of Pb behavior in sulfide may hold the key to several long-standing and important Pb paradoxes and enigmas. The classical Pb isotope paradox arises from the fact that all known mantle reservoirs lie to the right of the Geochron, with no consensus as to the identity of the “balancing” reservoir. We propose that long-term segregation of sulfide (containing Pb) to the core may resolve this paradox. Another Pb paradox arises from the fact that the Ce/Pb ratio of both OIB and MORB is greater than bulk earth, and constant at a value of 25. The constancy of this “canonical ratio” implies similar partition coefficients for Ce and Pb during magmatic processes (Hofmann et al. in Earth Planet Sci Lett 79:33–45, 1986), whereas most experimental studies show that Pb is more incompatible in silicates than Ce. Retention of Pb in residual mantle sulfide during melting has the potential to bring the bulk partitioning of Ce into equality with Pb if the sulfide melt/silicate melt partition coefficient for Pb has a value of ∼ 14. Modeling shows that the Ce/Pb (or Nd/Pb) of such melts will still accurately reflect that of the source, thus enforcing the paradox that OIB and MORB mantles have markedly higher Ce/Pb (and Nd/Pb) than the bulk silicate earth. This implies large deficiencies of Pb in the mantle sources for these basalts. Sulfide may play other important roles during magmagenesis: (1) advective/diffusive sulfide networks may form potent metasomatic agents (in both introducing and obliterating Pb isotopic heterogeneities in the mantle); (2) silicate melt networks may easily exchange Pb with ambient mantle sulfides (by diffusion or assimilation), thus “sampling” Pb in isotopically heterogeneous mantle domains differently from the silicate-controlled isotope tracer systems (Sr, Nd, Hf), with an apparent “de-coupling” of these systems.     

Groundwater acidification in Triassic sandstones: prediction with MAGIC modelling

 Acidification of groundwater lags behind acid deposition due to the relatively long water residence time in conjunction with various buffering processes in the soil zone and deeper aquifer (chemical weathering, cation exchange, sulfate sorption, and N uptake by the biomass). Extensive field data from eight forested catchments in the Bunter Sandstone of the Black Forest, including results from water budget studies and hydrochemical analysis of stream and spring waters, were used to simulate the future evolution of ground-water acidification with the MAGIC model. The present acid deposition exceeds the “critical load” (here meaning buffering due to chemical weathering and protonation of organic acids) in six of eight catchments. Two catchments are well buffered because they contain carbonate-bearing layers in the Upper Bunter sandstone. Transient buffering (i.e., cation exchange, N uptake, the sulfate sorption) thus far prevents worse acidification, but this effect will decline in the future. For one of the poorly buffered catchments (Seebach), a two-layer simulation was carried out, based on extensive data from 10 years of measurements. Validation of the long-term simulations by hydrochemical and soil data was hampered by strong annual variations but generally supported by paleolimnological studies. In the future, reductions in the S deposition by 20% and the N deposition by 10% up to the year 2030 are assumed as the most probable scenario. N uptake through soil and vegetation will come to an end as suggested by decreasing C/N ratios of the organic matter. This process is arbitrarily included in the simulations. In the periglacial soil layer, acidification will decrease until the year 2030 and then approach a steady-state condition. In the fractured aquifer, acidification will also proceed at a decreasing rate; however, sulfate desorption up to the year 2130, the end of simulated period, will prevent earlier remediation. Despite a significant reduction in S deposition since the mid-1980s, further efforts are necessary to reduce the emission of acidifying substances. Liming in the recharge area is partially effective to ameliorate “shallow” groundwater but largely fails to ameliorate “deeper” groundwater in the sandstone aquifer. Received: 30 July 1996/Accepted: 23 January 1997     

Survival of the thinnest: rediscovery of Bauer’s (1898) ichthyosaur tooth sections from Upper Jurassic lithographic limestone quarries,south Germany

During the late Early Cretaceous, the shallow-water domains of the western Tethys are characterized by the widespread deposition of Urgonian-type carbonates rich in rudists, corals and other oligotrophic, shallow-marine organisms. In the Helvetic Alps, the Urgonian occurrences have been dated by ammonite biostratigraphy as Late Barremian and Early Aptian. For the more proximal occurrences in the western Swiss Jura, a recent age model based on bio-, chemo- and sequence stratigraphy has been proposed, which allows for an improved correlation with the Helvetic counterparts. In order to corroborate the recently proposed age model for the Jura, a set of well-preserved rhynchonellids collected from five different lithostratigraphical formations and members (“Marnes bleues d’Hauterive”, “Marnes d’Uttins”, basal marly layers within the “Urgonien Jaune”, “Marnes de la Russille”, “Urgonien Blanc”) has been analysed for its strontium–isotope ratios (⁸⁷Sr/⁸⁶Sr). In addition, K–Ar dating was performed on well-preserved glauconite grains from two different levels (“Marnes d’Uttins” and a basal layer within the “Urgonien Jaune”). The correlation of the Sr–isotope data set with a belemnite-based, ammonite-calibrated reference curve provides an age model which is coherent with recently published ages based on calcareous nannofossil biostratigraphy and the correlation of trends in chemo- and sequence stratigraphy. K–Ar dating on well-preserved glauconite grains from the “Marnes d’Uttins” and lowermost part of the “Urgonien Jaune” delivered ages of 127.5 ± 2.3 and 130.7 ± 2.6 Ma, respectively. Whereas the age of the glauconitic level near the base of the “Urgonien Jaune” is chronostratigraphically meaningful, the K–Ar age of the “Marnes d’Uttins” appears too young relative to the presently used time scale. This may be related to rejuvenation of the K–Ar chronometer due to post-depositional Ar loss, most likely during hardground formation. The ages obtained here confirm the Late Barremian age for the onset of the Urgonian platform, an age which is conform with ages obtained in the Helvetic Alps and elsewhere along the northern Tethyan margin.     

On the local multiscale determination of the Earth’s disturbing potential from discrete deflections of the vertical

As a first approximation, the Earth is a sphere; as a second approximation, it may be considered an ellipsoid of revolution. The deviations of the actual Earth’s gravity field from the ellipsoidal “normal” field are so small that they can be understood to be linear. The splitting of the Earth’s gravity field into a “normal” and a remaining small “disturbing” field considerably simplifies the problem of its determination. Under the assumption of an ellipsoidal Earth model, high observational accuracy is achievable only if the deviation (deflection of the vertical) of the physical plumb line, to which measurements refer, from the ellipsoidal normal is not ignored. Hence, the determination of the disturbing potential from known deflections of the vertical is a central problem of physical geodesy. In this paper, we propose a new, well-promising method for modelling the disturbing potential locally from the deflections of the vertical. Essential tools are integral formulae on the sphere based on Green’s function with respect to the Beltrami operator. The determination of the disturbing potential from deflections of the vertical is formulated as a multiscale procedure involving scale-dependent regularized versions of the surface gradient of the Green function. The modelling process is based on a multiscale framework by use of locally supported surface curl-free vector wavelets.      

Garnet and diopside-bearing diamondites (framesites)

Summary ?Rocks consisting almost entirely of diamonds (diamondites) that contain minor amounts of silicates were analyzed for trace element abundances in the silicates by Laser Ablation ICP Mass Spectrometry for the first time. Diamondites, previously described as polycrystalline diamond “aggregates” and “framesite”, extend the range of monomineralic rocks known from the Earth’s upper mantle. Our samples are intergrowths of diamonds with abundant open cavities and some interstitial silicates. The most common silicate is pyrope which occurs in two different colors (and chemical compositions): orange and lilac similar to garnet inclusions in diamonds and garnets known from upper mantle eclogites and garnet peridotites, respectively. In our sample, the “peridotitic” garnet is accompanied by Cr-rich diopside whereas the “eclogitic” garnet is unaccompanied. Trace element abundances suggest that both types of garnet formed from upper mantle fluids of similar origin which were rich in a carbonatitic component. The diamondites likely formed from the same fluids. Diamonds precipitated first and – in smaller amounts – contemporaneously with the silicates. Major upper mantle minerals like olivine, orthopyroxene and omphacite are missing, possibly indicating that these minerals behaved as refractory phases and were not mobilized by fluids. The chemical composition of “eclogite” and “peridotite” garnets differ in Cr and high field strength elements contents but not in the moderately compatible elements. They also have the same low Fe/Mg ratio which indicates a peridotitic source for the fluids. The compositional difference in minor and trace elements appears to be the result of different fluid processing rather than of a different source, i.e., peridotite or eclogite.

---

Zusammenfassung Granat- und Diopsid-führende Diamantite (Framesite) Silikat-führende Gesteine aus Diamanten (Diamantite) wurden erstmals mittels Laser Ablation ICP Massenspektrometrie auf ihre Gehaltean Spurenelementen analysiert. Diese ursprünglich als polykristalline Diamant-“Aggregate” und Framesite beschriebenen Diamantite werden den monomineralischen Gesteinen des Erdmantels zugeordnet. Unsere Proben sind Verwachsungen von Diamanten mit verbreitet offenen Hohlr?umen und wenig interstitialen Silikaten. Das verbreitetste Silikat ist Pyrop, welcher zwei verschiedene Farben und chemische Zusammensetzungen hat: orange und fliederfarben, ?hnlich den Granat-Einschlüssen in Diamanten und den Granaten aus Erdmantel-Eklogiten und Peridotiten. In unserer Probe ist der “peridotitische” Granat von Cr-Diopsid begleitet, der “eklogitische” Granat ist unbegleitet. Die Spurenelement-H?ufigkeiten in beiden Granaten machen es wahrscheinlich, da? diese von Erdmantel-Fluiden gebildet wurden, die reich an einerkarbonatitischen Komponente waren und aus ?hnlichen Quellen stammten. Die Diamanten wurden wahrscheinlich von denselben Fluiden gebildet. Diamanten wurden zuerst ausgef?llt und wuchsen – in geringerem Ausma?– auch gemeinsam mit den Silikaten. Hauptminerale des Erdmantels wie Olivin, Orthopyroxen und Omphazit fehlen. M?glicherweise verhielten sich diese Minerale refrakt?r und wurden von den Fluiden nicht mobilisiert. Die chemische Zusammensetzung von “eklogitischem” und “peridotitischem” Granat unterscheidet sich in den Gehaltenan Cr und den Elementen, welche lonen hoher Feldst?rke bilden, aber nicht in den moderat kompatiblen Elementen. Die beiden Granate haben auch das gleiche niedrige Fe/Mg-Verh?ltnis, welches auf eine peridotitische Quelle für die Fluide hinweist. Die Unterschiede in der H?ufigkeit von Neben- und Spurenelementen in den Granaten k?nnte auf verschiedene Entwicklung der Fluide zurückzuführen sein und wahrscheinlich nicht auf eine Herkunft der Fluide aus verschiedenen Quellen, wie Peridotiten oder Eklogiten.

---

---

Received July 10, 1999;/revised version accepted November 8, 1999     

Thermal instability of the fluid column in a borehole: application to the Yaxcopoil hole (Mexico)

Observational evidence proved that even when a borehole is in “fully” stabilized conditions, temperature data may exhibit certain unrest resembling irregular oscillations in the order of hundredths or (in the extreme case) even tenths of degree. Temperature was monitored in complicated hydrogeological conditions in the Yaxcopoil-1 hole (Chicxulub impact structure, Mexico). Two experiments are reported: (a) 20-day monitoring when a logger was located in the center of the high temperature gradient anomaly produced by the cold wave slowly propagating downwards and (b) simultaneous three-loggers 18-day monitoring with loggers located above, in and below the anomaly. All observed temperature–time series displayed intermittent oscillations of temperature with sharp gradients and large fluctuations over all observed time scales. While the “upper” and “lower” records revealed quasi-periodic temperature variations, the “central” record shows fast temperature oscillations with strong up-and-down reversals, all with amplitudes up to a few tenths of degree. The observed temperature–time series were processed by recurrence and recurrence interval quantification as well as by spectral analyses. It is shown that fluid in a borehole, subject to thermal gradient, is stable, as far as the gradient remains below a certain critical value. At higher Rayleigh numbers, the periodic character of oscillations typical for “quiescent” regime is superseded by stochastic features. This “oscillatory” convection occurs due to instability of the horizontal boundary layers. In the specific case of the Yaxcopoil hole, the time series above and below the cold wave (characterized by relatively lower temperature gradients between 20 and 50 mK/m) contain a clear low frequency component produced by tidal forcing. This component dominates over the high frequency domain (periods from 10–15 to 1 min), which exhibit a scaling behavior. This pattern conspicuously changes in the center part of the cold wave, where the local temperature gradient exceeds 200 mK/m and where tidal forcing composes only ~3% of the signal.     

Bed sediment-associated trace metals in an urban stream, Oahu, Hawaii

 Of the 117 stream and lake systems sampled nationwide, fish from Manoa Stream on Oahu, Hawaii, have consistently shown the highest Pb concentrations. Therefore a detailed study was conducted to examine total metal contents in bed sediments from a 5.8-km stretch of Manoa Stream. A total of 123 samples (<63 μm) were examined for 18 elements and 14 samples for 21 elements. Selected samples were also examined using different leach solutions to examine metal phase associations. All trace metal data, computations of enrichment ratios and the modified index of geoaccumulation point to mineralogical control for Cr and Ni; minor anthropogenic contamination for Ba, Cd, Cu, Hg and Zn; and a very strong contamination signal for Pb. Maximum Pb contents (up to 1080 mg kg⁻¹) were associated with anthropogenic material dumping in minor tributaries, storm sewer sediments and sediments in the “lower” section of the basin. Proportionally Pb had the highest non-residual component of elements examined; dominantly in the reducible phase associated with Mn and amorphous Fe oxyhydroxides. The contamination signal was typically lowest in the “undisturbed” headwater reach of the basin (above 5.1 km) with significant increases throughout the “residential” and “commercial-institutional” zones of the mid-basin. The spatial pattern of bed sediment contamination and evidence from storm sewer-outlet sediments strongly indicates that Pb, and to a lesser degree some other metals, is still being transported to the stream and the primary agent is soil erosion and transport of metals sorbed to sediments. The primary source of sediment-associated metals is considered to be the automobile, though other minor sources can not be ruled out. Received: 3 November 1998 · Accepted: 26 January 1999     

Calculation of the configurational entropy of Al, Si in layer silicates using the cluster variation method

Entropies of Al-Si in layer silicates have been calculated using a series of CVM approximations for the honeycomb lattice. The parameters of the models have been constrained by ²⁹Si NMR data. The results of low order approximations such as “pair” and “star” have been rejected because of their low accuracy at high Al/(Al+Si) ratios. Reasonably accurate results have been achieved with the help of the “hexagon” and “star-hexagon” approximations. Received: 31 October 1997 / Revised, accepted: 17 April 1998     

Some structural features of the convective-velocity field in the solar photosphere

An algorithm for measuring horizontal photospheric velocities previously employed to process aerospace images is adapted for problems in solar physics and realized in a computational code. It differs from the standard procedure of local correlation tracking in a special choice of trial areas (“targets”), whose displacements are determined bymaximizing the correlation between the original and various shifted positions of the target. Specifically, an area is chosen as a target in a certain neighborhood of each node of a predefined grid if either the contrast or the entropy of the brightness distribution reaches its maximum in this area. The horizontal velocities obtained are then interpolated to the positions of imaginary “corks” using the Delaunay triangulation and affine transformations specified by the deformation of the obtained triangles at the time step considered. The motion of the corks is represented by their trajectories. A superposition of flows on different scales, from mesogranular to supergranular, can clearly be seen. “Large mesogranules” with sizes of order 15 Mm are revealed. In many cases, these are stellate in shape. Areas of strong convergence of the horizontal flows are detected; this convergence is sometimes accompanied by swirling. Evidence is found for the possible coexistence of convection cells with different circulation directions, so-called l-type and g-type cells.