Numerical modeling of two-phase hysteresis combined with an interface condition for heterogeneous porous media

This paper presents a numerical implementation of two-phase capillary hysteresis and its combination with a capillary interface condition for the treatment of heterogeneities. The hysteresis concepts chosen in this work are first implemented in a node-centered FV discretization scheme and subsequently combined with the interface condition that predicts sharp saturation discontinuities at material interfaces, based on a pressure equilibrium concept. This approach allows for the approximation of history-dependent, and at the same time discontinuous, saturations at material interfaces. The resulting model provides a well-defined evolution of the hysteretic capillary pressure–saturation relationships at material interfaces that is independent of the grid spacing. As demonstrated with a simple 1-D example, this concept therefore offers the advantage that the solution of a two-phase flow problem involving hysteresis does not relate to the grid resolution at the material interfaces.     

Garnet reaction rims from the breakdown of Staurolite in polymetamorphic micashists from the Rappold complex, Austroalpine basement, Eastern Alps

In polymetamorphic pelites of the Rappold complex in the Wölz crystalline basement (Eastern Alps) reaction rim garnets at staurolite-quartz interfaces (type I) and single grain garnets along previous staurolite-white mica interfaces (type II) were formed. The garnet reaction rims were formed during the Cretaceous amphibolite facies metamorphic overprint of the pre-existing mineral assemblages comprising garnet, staurolite, and kyanite from an amphibolite facies metamorphic event probably of Variscian age. The newly formed garnet may take the form of reaction rims along the margins of large pre-existing staurolite blasts. The initial growth increments of garnet have low grossular content, and reaction rim growth was controlled by the transfer of Fe, Mg and Mn components from the staurolite-garnet interface to the quartz-garnet interface. Later garnet growth increments have relatively high grossular content due to consumption of matrix plagioclase, which was destabilized by successive pressure increase. The grossular content of newly formed garnet shows systematic increase towards sites where plagioclase breaks down indicating that transport of calcium through the matrix was sluggish. On the basis of reaction microstructures it is demonstrated that the mineral assemblage garnet?+?kyanite?+?biotite?+?paragonite was formed at the conditions of eo-alpine amphibolite facies overprint while staurolite and plagioclase broke down successively with increasing pressure.     

Compositional zoning of garnet porphyroblasts from the polymetamorphic Wölz Complex, Eastern Alps

We employ garnet isopleth thermobarometry to derive the P–T conditions of Permian and Cretaceous metamorphism in the Wölz crystalline Complex of the Eastern Alps. The successive growth increments of two distinct growth zones of the garnet porphyroblasts from the Wölz Complex indicate garnet growth in the temperature interval of 540°C to 560°C at pressures of 400 to 500 MPa during the Permian and temperatures ranging from 550°C to 570°C at pressures in the range of 700 to 800 MPa during the Cretaceous Eo-Alpine event. Based on diffusion modelling of secondary compositional zoning within the outermost portion of the first garnet growth zone constraints on the timing of the Permian and the Eo-Alpine metamorphic events are derived. We infer that the rocks remained in a temperature interval between 570°C and 610°C over about 10 to 20 Ma during the Permian, whereas the high temperature stage of the Eo-Alpine event only lasted for about 0.2 Ma. Although peak metamorphic temperatures never exceeded 620°C, the prolonged thermal annealing during the Permian produced several 100 µm wide alteration halos in the garnet porphyroblasts and partially erased their thermobarometric memory. Short diffusion profiles which evolved around late stage cracks within the first garnet growth zone constrain the crack formation to have occurred during cooling below about 450°C after the Eo-Alpine event.     

Late Glacial ice advances in southeast Tibet

The sensitivity of Tibetan glacial systems to North Atlantic climate forcing is a major issue in palaeoclimatology. In this study, we present surface exposure ages of erratic boulders from a valley system in the Hengduan Mountains, southeastern Tibet, showing evidence of an ice advance during Heinrich event 1. Cosmogenic nuclide analyses (¹⁰Be and ²¹Ne) revealed consistent exposure ages, indicating no major periods of burial or pre-exposure. Erosion-corrected (3 mm/ka) ¹⁰Be exposure ages range from 13.4 to 16.3 ka. This is in agreement with recalculated exposure ages from the same valley system by [Tschudi, S., Schäfer, J.M., Zhao, Z., Wu, X., Ivy-Ochs, S., Kubik, P.W., Schlüchter, C., 2003. Glacial advances in Tibet during the Younger Dryas? Evidence from cosmogenic ¹⁰Be, ²⁶Al, and ²¹Ne. Journal of Asian Earth Sciences 22, 301–306.]. Thus this indicates that local glaciers advanced in the investigated area as a response to Heinrich event 1 cooling and that periglacial surface adjustments during the Younger Dryas overprinted the glacial morphology, leading to deceptively young exposure ages of certain erratic boulders.     

Distribution, abundance and life history of Mysis relicta (Lovén) in the Feldberg Lake District, Germany

Distribution, abundance and life history characteristics of Mysis relicta were studied in the Feldberg Lake District (Lake Breiter Luzin, Lake Schmaler Luzin, Lake Zansen) located in northeastern Germany. Between July 2001 and November 2002 mysids were collected by vertical net hauls. In order to determine the impact of the current trophic conditions on the distribution of mysids in these lakes, oxygen concentration, total phosphorus, chlorophyll a and water transparency were also measured. All investigated lakes are mesotrophic at present. Lake Breiter Luzin exhibited great seasonal and spatial variations in mysid abundance. Density of adults and juveniles had a mean of 44.9 ± 57.1 and 68.7 ± 99.6 m⁻², respectively. Highest abundance of adults was 110.4 ± 76.5 m⁻² in summer, lowest abundances of 2.0 ± 4.0 m⁻² occurred in spring. For juveniles, highest density of 218.4 ± 174.6 m⁻² was detected in summer and lowest of 0.8 ± 1.8 m⁻² in winter. No mysids were caught in any of the daytime hauls, but they were widely distributed throughout the water column at night. Size frequency distribution of mysids suggested that reproduction occurred year-round, the most consistent influx of juveniles occurred in early summer and a smaller second cohort in autumn. Highest mysid abundance was 189.2 ± 318.6 adults and 127.0 ± 66.3 juveniles m⁻² in Lake Schmaler Luzin, and 59.6 ± 5.6 adults and 79.4 ± 11.2 juveniles m⁻² in Lake Zansen. There were great spatial differences in abundance in both lakes.     

Laboratory measurements of the association rate coefficients of NO+, O+2, N+, and N+2 ions with N2 and CO2 at temperatures between 100 k and 400 k

Rate coefficients for the association reactions of NO⁺ ions with N₂ and CO₂, O₂⁺ with N₂, and N⁺ and N₂⁺ with N₂ have been determined as a function of gas temperature in a laboratory experiment employing a variable-temperature drift-tube apparatus. The measured rate coefficients were fitted to power laws of the form $\text{k = C}\text{(}\text{T}\text{300}\text{)}{}^{\mathrm{x}}$ where the exponents x ranged from 2.2 to 4.3. The strong temperature dependence observed in the case of the reaction of NO⁺ with N₂ (x = 4.3) supports the thesis by Arnold et al. (1979) that the temperature variability of D-region ion densities is a result of this reaction step in the ion clustering sequence.     

Isotopic and elemental fractionation of solar wind implanted in the Genesis concentrator target characterized and quantified by noble gases

Abstract– We report concentrations and isotopic compositions of He, Ne, and Ar measured with high spatial resolution along a radial traverse of a silicon carbide (SiC) quadrant of the Genesis mission concentrator target. The Ne isotopic composition maps instrumental fractionation as a function of radial position in the target: the maximum observed isotopic fractionation is approximately 33‰ per mass unit between the center and periphery. The Ne fluence is enhanced by a factor of 43 at the target center and decreases to 5.5 times at the periphery relative to the bulk solar wind fluence. Neon isotopic profiles measured along all four arms of the “gold cross” mount which held the quadrants in the concentrator target demonstrate that the concentrator target was symmetrically irradiated during operation as designed. We used implantation experiments of Ne into SiC and gold to quantify backscatter loss and isotopic fractionation and compared measurements with numerical simulations from the code “stopping and range of ions in matter.” The ²⁰Ne fluence curve as a function of radial distance on the target may be used to construct concentration factors relative to bulk solar wind for accurate corrections for solar wind fluences of other light elements to be measured in the concentrator target. The Ne isotopic composition as a function of the radial distance in the SiC quadrant provides a correction for the instrumental mass‐dependent isotopic fractionation by the concentrator and can be used to correct measured solar wind oxygen and nitrogen isotopic compositions to obtain bulk solar wind isotopic compositions.