Comparison of Galerkin-type discretization techniques for two-phase flow in heterogeneous porous media

A comparison of Standard Galerkin, Petrov-Galerkin, and Fully-Upwind Galerkin methods for the simulation of two-phase flow in heterogeneous porous media is presented. On the basis of the coupled pressure-saturation equations, a generalized formulation for all three finite element methods is derived and analysed. For flow in homogeneous media, the Petrov-Galerkin method gives excellent results. But this method fails miserably for problems with heterogeneous media. This is because it is not able to capture correctly processes that take place at interfaces when, for instance, the capillary pressure-saturation relationship after Brooks and Corey is assumed. The Fully-Upwind Galerkin method is superior to the Petrov-Galerkin approach because it is able to give correct results for flow in homogeneous and heterogeneous media for the two models of van Genuchten and Brooks-Corey. The widely used formulation which is correct for the homogeneous case cannot be used for heterogeneous media. Instead the straightforward approach of gradp_c in combination with a chord-slope technique must be utilized.     

Node-centered finite volume discretizations for the numerical simulation of multiphase flow in heterogeneous porous media

Three node-centered finite volume discretizations for multiphase porous media flow are presented and compared. By combination of these methods two additional discretization methods are generated. The ability of these schemes to describe flows at textural interfaces of different geologic formations is investigated. It was found that models with nonzero-entry pressures for the capillary pressure-saturation relationship in conjunction with the Box discretization may give rise to spurious oscillations for flows around low permeable lenses. Furthermore, the applicability and sensitivity of the discretization methods with regard to the used computational grids is discussed. The schemes are used for the numerical study of two-phase flow in porous media with zones of different material properties. This revised version was published online in July 2006 with corrections to the Cover Date.     

Methane exchange between coal-bearing basins and the atmosphere: the Ruhr Basin and the Lower Rhine Embayment, Germany

A precise knowledge of methane exchange processes is required to fully understand the recent rise of atmospheric methane concentration. Three of these processes take place at the lithosphere/atmosphere boundary: bacterial consumption of methane and emission of bacterial or thermogenic methane. This study was initiated to quantify these processes on a regional scale in the Ruhr Basin and the Lower Rhine Embayment. Since these areas are subject to bituminous coal and lignite mining, natural and anthropogenically-induced methane exchange processes could be studied. The methane emission and consumption rates and their carbon isotope signal were measured at the lithosphere/atmosphere boundary using flux chambers. On most of the soils studied, methane consumption by bacteria was identified. Thermogenic methane was released only at some of the natural faults examined. In active and abandoned bituminous coal mining areas methane emissions were restricted to small areas, where high emission rates were measured. The carbon isotope composition of methane at natural faults and in mining subsidence troughs was typical of thermogenic methane (−45 to −32 ‰ δ¹³C). Methane exchange balancing revealed that natural methane emissions from these two basins represent no source of atmospheric importance. However, methane release by upcast mining shafts dominates the methane exchange processes and is by about two orders of magnitude greater than methane consumption by bacterial oxidation in the soils.     

Assimilation and diffusion during xenolith-magma interaction: a case study of the Variscan Karkonosze Granite, Bohemian Massif

Magmatic enclaves from the Rudolfov quarry near Liberec (Czech Republic) are interpreted to represent remnants of lamprophyric melt that intruded the Karkonosze granite at a stage at which the granite was not fully solidified. Based on the observation that larger enclaves are generally more circular than the smaller ones, we conclude that bigger blobs of mafic magma became more spherical during flow in the gravity field (sink or float). This flow is also interpreted to be responsible for the incorporation of mineral grains into the enclaves and may have facilitated the assimilation of granitic melt. Linear mixing trends on Harker diagrams for most network-forming and mainly slow-diffusing or fluid-immobile elements indicate such an assimilation process between granite and lamprophyre. In contrast, all fast-diffusing or fluid-mobile elements display scattered element distributions, implying that chemical diffusion also played a role. Pressure and temperature for this late magmatic stage are estimated at around 1 kbar and 500°C. These results suggest that two processes modified the composition of the enclaves in the Karkonosze granite: (1) assimilation (mechanical mixing) of granitic melt during the injection of the lamprophyric melt and the subsequent flow of the forming enclaves in the gravity field (responsible for the linear mixing trends) and (2) diffusion-controlled redistribution of elements between both solidifying rock types at the magmatic stage and/or due to late-stage magmatic fluids (responsible for the scattering and deviation from the linear mixing trends).     

Formation and chemical evolution of magnesium chloride brines by evaporite dissolution processes—Implications for evaporite geochemistry

The evolution of magnesium chloride brines with high bromide contents via a multistage reaction and dissolution process has been studied in brine seeps of a German potash mine. The observed chemical trends and phase equilibria can be modeled and interpreted in terms of a NaCl solution (cap rock brine) infiltrating into a potash zone characterized by the metamorphic mineral assemblage kieserite + sylvite + halite + anhydrite. Establishment of a persistent, stable equilibrium assemblage and constant fluid composition in the invariant point IP1 of the six component (Na-K-Mg-Ca-Cl-SO₄-H₂O) system of oceanic salts is prevented by the perpetually renewed input of NaCl-brine and by the intermittent exposure of incompatible kieserite. Instead, the solutions develop towards the metastable invariant point IP1(gy), with the mineral assemblage carnallite + polyhalite + sylvite + halite + gypsum, where gypsum takes the place of anhydrite (stage I). The temporary exposure of kieserite and the ensuing formation of polyhalite effectively buffer the solutions along the metastable polyhalite phase boundary during stages II and III. Eventually, in stage IV, polyhalite becomes depleted and admixture of more NaCl brine leads to low sulfate solution compositions, which are now only constrained by carnallite + sylvite + halite, and the once hexary system degenerates to a quaternary one (Na-K-Mg-Cl-H₂O) in point E. Bromide in brines shows equilibrium partitioning with respect to the wall rock minerals. The pattern of evolving brine compositions may serve as a model for similar brine occurrences, which in some cases may have been misinterpreted as remains of fossil, highly concentrated and chemically modified seawater. Similar magnesium chloride brines of salt lakes (e.g., Dead Sea, Dabusun Lake) show subtle differences and are constrained by fewer mineral equilibria (more degrees of freedom), and their low sulfate contents are due to gypsum precipitation, driven by calcium chloride input from dolomitization reactions. Finally, the observed reaction sequence is generalized, and a model for the formation of magnesium sulfate depleted, chloride-type potash salts and bischofite deposits by leaching of sulfate-type evaporites is proposed.     

Instruments, Detectors and the Future of Astronomy with Large Ground Based Telescopes

Results of a survey of instrumentation and detector systems, either currently deployed or planned for use at telescopes larger than 3.5 m, in ground based observatories world-wide, are presented. This survey revealed a number of instrumentation design trends at optical, near, and mid-infrared wavelengths. Some of the most prominent trends include the development of vastly larger optical detector systems (> 10⁹ pixels) than anything built to date, and the frequent use of mosaics of near-infrared detectors – something that was quite rare only a decade ago in astronomy. Some future science applications for detectors are then explored, in an attempt to build a bridge between current detectors and what will be needed to support the research ambitions of astronomers in the future.     

Results from a Multi-disciplinary Sedimentary Pilot Study of Tectonic Lake Iznik (NW Turkey) – Geochemistry and Paleolimnology of the Recent Past

A limnogeological reconnaissance study was carried out on Lake Iznik, located in the southeast of the Marmara region of Turkey, involving a seismic survey and collection of short sediment cores. This lake is located on the middle branch of the North Anatolian Fault Zone (NAFZ), a transform plate boundary between the Eurasian and Anatolian Plates. It is, therefore, tectonically active and offers an opportunity to investigate the interplay of sedimentary and seismo-tectonic processes, as well as climate change and human impact in the region. Short cores of the three sub-basins, maximum length of 35.5 cm, recovered non-laminated, blackish clays and silts with varying amounts of biogenic and minerogenic (allochthonous, autochthonous) material, which documented almost the last 80 years of deposition and environmental history. High sedimentation rates in the deeper core sections are accompanied by changes in land use (conversion of woodland to farmland) in the northern areas of Lake Iznik, which caused the deposition of more weathered material (high K/Na ratios) and higher contents of Mn in the lake. A tendency towards eutrophic conditions within the last 20 years is indicated by high nutrient content (N, TOC, P), decreasing C/N-ratios, and characteristic diatom and cladoceran associations. Also increased pollution is revealed by higher Pb, Cu, and Zn contents and increased supply of human and animal faeces (high coprostanol content) during the last two decades. But simultaneous lower sedimentation rates towards the core tops complicate the reconstruction of recent and past eutrophication and pollution states of Lake Iznik. This requires an extension of the pilot study and deeper sediment cores, to recover non-anthropogenic influenced sediment levels.     

Volume diffusion of Ytterbium in YAG: thin-film experiments and combined TEM–RBS analysis

In this study, we address volume diffusion of ytterbium in yttrium aluminum garnet (YAG) using thin-film single crystal diffusion couples. We employ analytical transmission electron microscopy (ATEM) as a tool for combined microstructural and microchemical analysis and compare the results to Rutherford backscattering (RBS) analysis. Given the high spatial resolution of the method, we focus on microstructural changes of the thin-film diffusant source during the diffusion anneal. We evaluate the potential influence of the associated changes in its transport properties on the evolution of concentration profiles in the single crystal substrate. This approach allows us to test the reliability of determination of volume diffusion coefficients from thin-film diffusion experiments. We found that for the chosen experimental setting, the influence of thin-film re-crystallization is small when compared with the experimental uncertainty and good estimates for the volume diffusion coefficients of Yb in YAG can be obtained using standard assumptions. Both Yb-concentration profiles analyzed with ATEM and with RBS give similar results. At 1,450°C and 1 bar, we infer log D _Yb (m²/s) values of −19.37 ± 0.07 (TEM) and −19.84 ± 0.02 (RBS). Although the change in thin-film transport properties associated with successive crystallization during the diffusion anneal does not play a major role for our experimental setup, this effect cannot generally be ignored.     

The first galaxies: assembly, cooling and the onset of turbulence

We investigate the properties of the first galaxies at   z ≳ 10  with highly resolved numerical simulations, starting from cosmological initial conditions and taking into account all relevant primordial chemistry and cooling. A first galaxy is characterized by the onset of atomic hydrogen cooling, once the virial temperature exceeds  ≃10⁴ K  , and its ability to retain photoheated gas. We follow the complex accretion and star formation history of a  ≃5 × 10⁷ M_⊙  system by means of a detailed merger tree and derive an upper limit on the number of Population III (Pop III) stars formed prior to its assembly. We investigate the thermal and chemical evolution of infalling gas and find that partial ionization at temperatures  ≳10⁴ K  catalyses the formation of  H₂  and hydrogen deuteride, allowing the gas to cool to the temperature of the cosmic microwave background. Depending on the strength of radiative and chemical feedback, primordial star formation might be dominated by intermediate-mass Pop III stars formed during the assembly of the first galaxies. Accretion on to the nascent galaxy begins with hot accretion, where gas is accreted directly from the intergalactic medium and shock heated to the virial temperature, but is quickly accompanied by a phase of cold accretion, where the gas cools in filaments before flowing into the parent halo with high velocities. The latter drives supersonic turbulence at the centre of the galaxy and could lead to very efficient chemical mixing. The onset of turbulence in the first galaxies thus likely marks the transition to Pop II star formation.     

Freshwater snail communities and lake classification. An example from the Åland Islands, southwestern Finland

The freshwater snail fauna in 51 lakes on the Åland Islands was investigated. By means of cluster analysis it appears that lakes may be divided into four groups on the basis of snail species composition. The lakes of the four groups show a spatial distribution which is related to topographical differences, terrestrial vegetation and land use. Snails may be divided on the basis of water hardness into hard water species and calcium indifferent species or into demanding and modest/indifferent species when clustering is based on presence/absence.