Construction of three-phase data to model multiphase flow in porous media: Comparing an optimization approach to the finite element approach

Multiphase flow modelling is a major issue in the assessment of groundwater pollution. Three-phase flows are commonly governed by mathematical models that associate a pressure equation with two saturation equations. These equations involve a number of secondary variables that reflect the fluid behaviour in a porous medium. To improve the computational efficiency of multiphase flow simulators, several simplified reformulations of three-phase flow equations have been proposed. However, they require the construction of new secondary variables adapted to the reformulated flow equations. In this article, two different approaches are compared to quantify these variables. A numerical example is given for a typical fine sand.     

Präparationstechnik für Aufbereitungsprodukte zur quantitativen erzmikroskopischen Untersuchung

Zusammenfassung Mikroskopische Anschliffe feinkörniger Aufbereitungsprodukte können nur dann quantitativ ausgewertet werden, wenn ihre Schliffqualität so hoch ist, daß auch feinste Partien sicher zu bestimmen sind. Daneben muß gewährleistet sein, daß die an der Schliffoberfläche erfaßten Körner die tatsächliche Verteilung der Mineralien der Probe repräsentieren. Diese Forderung kann nur durch zweckentsprechende Präparation, die dem eigentlichen Schleif- und Polierprozeß vorausgeht, erfüllt werden. Als günstigste Lösung erscheint es, das Korngemisch exakt als Einkornschicht zu präparieren und anschließend wie Stückproben anzuschleifen. Die Schliffe können dann nach den gleichen statistischen Gesetzmäßigkeiten ausgewertet werden, wie sie bei praktisch zweidimensionalen Dünnschliffen oder bei Streupräparaten nichtopaker Korngemische vorliegen.

---

Quantitative examination of fine-grained ore dressing products depends on good quality of the sections and on whether very small particles and areas of grains can be determined. Besides, it must absolutely be guaranteed that the grains at the section surface represent the actual distribution of the minerals in the sample. This condition can only be realized by suitable grain mounting prior to grinding and polishing. The best method for this purpose appears to be the preparation of a mono-grain layer. The polished sections can then be examined according to the same statistical laws as two-dimensional thin sections or immersion-liquid slides of non-opaque mixtures of grains.

---

---

Mitteilung Nr. 141 aus dem Forschungsinstitut für Aufbereitung der Deutschen Akademie der Wissenschaften zu Berlin, Freiberg (Sachsen). Direktor: Prof. Dr.-Ing. habil. H. Kirchberg.     

Holocene hydrological variability of Lake Ladoga,northwest Russia,as inferred from diatom oxygen isotopes

This article presents a new comprehensive assessment of the Holocene hydrological variability of Lake Ladoga, northwest Russia. The reconstruction is based on oxygen isotopes of lacustrine diatom silica (δ¹⁸O_diatom) preserved in sediment core Co 1309, and is complemented by a diatom assemblage analysis and a survey of modern isotope hydrology. The data indicate that Lake Ladoga has existed as a freshwater reservoir since at least 10.8 cal. ka BP. The δ¹⁸O_diatom values range from +29.8 to +35.0‰, and relatively higher δ¹⁸O_diatom values around +34.7‰ between c. 7.1 and 5.7 cal. ka BP are considered to reflect the Holocene Thermal Maximum. A continuous depletion in δ¹⁸O_diatom since c. 6.1 cal. ka BP accelerates after c. 4 cal. ka BP, indicating Middle to Late Holocene cooling that culminates during the interval 0.8–0.2 cal. ka BP, corresponding to the Little Ice Age. Lake‐level rises result in lower δ¹⁸O_diatom values, whereas lower lake levels cause higher δ¹⁸O_diatom values. The diatom isotope record gives an indication for a rather early opening of the Neva River outflow at c. 4.4–4.0 cal. ka BP. Generally, overall high δ¹⁸O_diatom values around +33.5‰ characterize a persistent evaporative lake system throughout the Holocene. As the Lake Ladoga δ¹⁸O_diatom record is roughly in line with the 60°N summer insolation, a linkage to broader‐scale climate change is likely.     

Redistribution of HFSE elements during rutile replacement by titanite

Titanite growth at the expense of rutile during retrograde hydration of eclogite into amphibolite is a common phenomenon. We investigated an amphibolite sample from the Tromsø eclogite facies terrain in Northern Norway to gain insight into the trace element distribution between rutile and titanite during incomplete resorption of the rutile by titanite. Patchy compositional zoning of Al, Ti, and F in titanite relates to the presence of a fluid with variable Ti/Al and/or F during its growth. Laser ablation ICP–MS and electron microprobe data for high field strength elements (HFSE: Nb, Zr, Ta, and Hf) of rutile resorbed by titanite indicate a pronounced enrichment of these elements in the rim of a large single rutile crystal (~8 mm) and a systematic decrease towards uniform HFSE contents in the large core. HFSE contents of smaller rutile grains (~0.5 mm) and rutile inclusions (<100 μm) in the titanite overgrowth are similar or higher than in the rims of large rutile crystals. Element profiles from the rim inward demonstrate that HFSE enrichment in rutile is controlled by diffusion. HFSE ratios in diffusion-altered rutile show systematic variations compared with the uniform core composition of the large rutile. Modelling of Zr and Nb diffusion in rutile indicates that diffusion coefficients in rutile in fluid-dominated natural systems must be considerably higher than those determined experimentally at 1 bar in dry systems. Variations of HFSE contents in the newly formed titanite show no systematic spatial distribution. HFSE ratios in titanite and the rims of rutile are different, indicating different solid/fluid distribution coefficients in these minerals. Element fractionation by diffusion into the relict rutile and during fluid-mediated growth of new titanite could substantially change the HFSE budget of these minerals and could affect their use for geochemical tracing and other applications, such as Zr-based geothermobarometry.     

Mutual validation of GNSS height measurements and high-precision geometric-astronomical leveling

The method of geometric-astronomical leveling is presented as a suited technique for the validation of GNSS (Global Navigation Satellite System) heights. In geometric-astronomical leveling, the ellipsoidal height differences are obtained by combining conventional spirit leveling and astronomical leveling. Astronomical leveling with recently developed digital zenith camera systems is capable of providing the geometry of equipotential surfaces of the gravity field accurate to a few 0.1 mm per km. This is comparable to the accuracy of spirit leveling. Consequently, geometric-astronomical leveling yields accurate ellipsoidal height differences that may serve as an independent check on GNSS height measurements at local scales. A test was performed in a local geodetic network near Hanover. GPS observations were simultaneously carried out at five stations over a time span of 48 h and processed considering state-of-the-art techniques and sophisticated new approaches to reduce station-dependent errors. The comparison of GPS height differences with those from geometric-astronomical leveling shows a promising agreement of some millimeters. The experiment indicates the currently achievable accuracy level of GPS height measurements and demonstrates the practical applicability of the proposed approach for the validation of GNSS height measurements as well as the evaluation of GNSS height processing strategies.     

Impact of loading displacements on SLR-derived parameters and on the consistency between GNSS and SLR results

Displacements of the Earth’s surface caused by tidal and non-tidal loading forces are relevant in high-precision space geodesy. Some of the corrections are recommended by the international scientific community to be applied at the observation level, e.g., ocean tidal loading (OTL) and atmospheric tidal loading (ATL). Non-tidal displacement corrections are in general recommended not to be applied in the products of the International Earth Rotation and Reference Systems Service, in particular atmospheric non-tidal loading (ANTL), oceanic and hydrological non-tidal corrections. We assess and compare the impact of OTL, ATL and ANTL on SLR-derived parameters by reprocessing 12 years of SLR data considering and ignoring individual corrections. We show that loading displacements have an influence not only on station long-term stability, but also on geocenter coordinates, Earth Rotation Parameters, and satellite orbits. Applying the loading corrections reduces the amplitudes of annual signals in the time series of geocenter and station coordinates. The general improvement of the SLR station 3D coordinate repeatability when applying OTL, ATL and ANTL corrections are 19.5 %, 0.2 % and 3.3 % respectively, w.r.t. the solutions without loading corrections. ANTL corrections play a crucial role in the combination of optical (SLR) and microwave (GNSS, VLBI, DORIS) space geodetic observation techniques, because of the so-called Blue-Sky effect: SLR measurements can be carried out only under cloudless sky conditions—typically during high air pressure conditions, when the Earth’s crust is deformed, whereas microwave observations are weather-independent. Thus, applying the loading corrections at the observation level improves SLR-derived products as well as the consistency with microwave-based results. We assess the Blue-Sky effect on SLR stations and the consistency improvement between GNSS and SLR solutions when ANTL corrections are included. The omission of ANTL corrections may lead to inconsistencies between SLR and GNSS solutions of up to 2.5 mm for inland stations. As a result, the estimated GNSS–SLR coordinate differences correspond better to the local ties at the co-located stations when applying ANTL corrections.     

On the electrical conductivity of air inside buildings

Summary Simultaneous measurements of the electrical conductivity of air inside and outside of buildings have been performed with continuously recording electronic instruments of high sensitivity and short time constant. It appears that the electrical conductivity, a measure of the natural ionization of air, is one atmospheric-electric parameter that can manifest itself inside occupied buildings without appreciable difference in physical character from outside values.The probable importance of these experimental results for the field of bioclimatology is pointed out and the possibility is discussed that the electrical conductivity of air might be one of the searched—for parameters in associating physiological phenomena with weather changes.

---

Zusammenfassung Gleichzeitige Messungen der luftelektrischen Leitfähigkeit im Freien und innerhalb von Gebäuden wurden mit kontinuierlich registrierenden Röhreninstrumenten von großer Empfindlichkeit und kurzer Zeitkonstante durchgeführt. Es wurde gefunden, daß die Leitfähigkeit als Maß der natürlichen Luftionisierung eine luftelektrische Größe darstellt, die ohne wesentliche physikalische Beeinflussung in bewohnte Gebäude einzudringen vermag.Die mögliche Bedeutung dieser experimentellen Ergebnisse für das Gebiet der Bioklimatologie wird angedeutet und die Wahrscheinlichkeit diskutiert, daß die luftelektrische Leitfähigkeit einen der noch unbestimmten Faktoren des Zusammenhanges physiologischer Erscheinungen mit Wetterveränderungen darstellen könnte.

---

Résumé Des instruments électroniques à enregistrement continu de haute fidelité et ayant une constante de courte durée ont permis de prendre des mesures simultanées de la conductivité électrique de l'air à l'intérieur et à l'extérieur de bâtiments. Il semble que la conductivité électrique, mesure de l'ionisation naturelle de l'air, soit un paramètre de l'électricité atmosphérique qui peut se manifester à l'intérieur de bâtiments habités sans différence appréciable de caractéristiques physiques avec celles des valeurs extérieures.L'importance probable de ces résultats expérimentaux dans le champ de la bioclimatologie est indiquée, et la possibilité que la conductivité atmosphérique puisse être un des paramètres recherchés dans l'association des phénomènes physiologiques avec des changements du temps est étudiée.

---

---

With 1 Figure.

---

Publication No. 25 of the Institute of Geophysics, University of California at Los Angeles. The research reported in this paper has been sponsored in part by the Geophysics Research Directorate of the Air Force Cambridge Research Center, Air Research and Development Command, under Contract No. AF 19(122)-254.     

Complete recycling of a magmatic arc: evidence from chemical and isotopic composition of Quaternary trench sediments in Chile (36°–40°S)

Marine Quaternary trench and slope sediments were sampled along the margin of the Southern Andes, Chile between 36° and 40°S. Major and trace element contents indicate only minor influence of weathering and transport fractionation. The whole rock composition of the sediments is similar to the average rock of the Cretaceous to Holocene magmatic arc of this section of the southern volcanic zone. Sr, Nd, and Pb isotope composition of the sediments also resembles closely the average composition of the magmatic arc. The contribution of compositionally distinct Palaeozoic crust, which makes up most of the volume of the forearc, is ~0–20% crustal Sr, Nd, and Pb according to the isotope record of the trench and slope sediments. Input of sediments from the continent into the subduction system was dominated by detritus from the magmatic arc at least for the last 20 My on the basis of the Oligocene to Holocene exhumation history of the margin.     

Distribution of chlorophyll a and diatom in the sea ice of the Weddell Sea, Antarctica

Diatoms are major primary producers of microbial biomass in the Antarctica. They are found in the water and sea ice. The distribution, abundance of the ice diatoms and their relation to the environmental factors inside and outside the ice have been studied for its special role in the Antarctic Ocean ecology. In this paper we describe the abundance, distribution and composition of diatom assemblages in