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.     

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.     

Experimental growth of åkermanite reaction rims between wollastonite and monticellite: evidence for volume diffusion control

Growth rates of monomineralic, polycrystalline åkermanite (Ca₂MgSi₂O₇) rims produced by solid-state reactions between monticellite (CaMgSiO₄) and wollastonite (CaSiO₃) single crystals were determined at 0.5 GPa dry argon pressure, 1,000–1,200°C and 5 min to 60 h, using an internally heated pressure vessel. Inert Pt-markers, initially placed at the monticellite–wollastonite interface, indicate symmetrical growth into both directions. This and mass balance considerations demonstrate that rim growth is controlled by transport of MgO. At 1,200°C and run durations between 5 min and 60 h, rim growth follows a parabolic rate law with rim widths ranging from 0.4 to 16.3 μm indicating diffusion-controlled rim growth. The effective bulk diffusion coefficient \( D_{\text{eff,MgO}}^{\text{Ak}} \) is calculated to 10^?15.8±0.1 m² s^?1. Between 1,000°C and 1,200°C, the effective bulk diffusion coefficient follows an Arrhenius law with E _a = 204 ± 18 kJ/mol and D ₀ = 10^?8.6±1.6 m² s^?1. Åkermanite grains display a palisade texture with elongation perpendicular to the reaction interface. At 1,200°C, average grain widths measured normal to elongation, increase with the square root of time and range from 0.4 to 5.4 μm leading to a successive decrease in the grain boundary area fraction, which, however, does not affect \( D_{\text{eff,MgO}}^{\text{Ak}} \) to a detectible extent. This implies that grain boundary diffusion only accounts for a minor fraction of the overall chemical mass transfer, and rim growth is essentially controlled by volume diffusion. This is corroborated by the agreement between our estimates of the effective MgO bulk diffusion coefficient and experimentally determined volume diffusion data for Mg and O in åkermanite from the literature. There is sharp contrast to the MgO–SiO₂ binary system, where grain boundary diffusion controls rim growth.     

It starts at home? Climate policies targeting household consumption and behavioral decisions are key to low-carbon futures

Through their consumption behavior, households are responsible for 72% of global greenhouse gas emissions. Thus, they are key actors in reaching the 1.5 °C goal under the Paris Agreement. However, the possible contribution and position of households in climate policies is neither well understood, nor do households receive sufficiently high priority in current climate policy strategies. This paper investigates how behavioral change can achieve a substantial reduction in greenhouse gas emissions in European high-income countries. It uses theoretical thinking and some core results from the HOPE research project, which investigated household preferences for reducing emissions in four European cities in France, Germany, Norway and Sweden. The paper makes five major points: First, car and plane mobility, meat and dairy consumption, as well as heating are the most dominant components of household footprints. Second, household living situations (demographics, size of home) greatly influence the household potential to reduce their footprint, even more than country or city location. Third, household decisions can be sequential and temporally dynamic, shifting through different phases such as childhood, adulthood, and illness. Fourth, short term voluntary efforts will not be sufficient by themselves to reach the drastic reductions needed to achieve the 1.5 °C goal; instead, households need a regulatory framework supporting their behavioral changes. Fifth, there is a mismatch between the roles and responsibilities conveyed by current climate policies and household perceptions of responsibility. We then conclude with further recommendations for research and policy.     

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.     

A diverse Pleistocene marsupial trackway assemblage from the Victorian Volcanic Plains,Australia

A diverse assemblage of late Pleistocene marsupial trackways on a lake bed in south-western Victoria provides the first information relating to the gaits and morphology of several megafaunal species, and represents the most speciose and best preserved megafaunal footprint site in Australia. The 60–110 ka volcaniclastic lacustrine sedimentary rocks preserve trackways of the diprotodontid Diprotodon optatum, a macropodid (probably Protemnodon sp.) and a large vombatid (perhaps Ramsayia magna or ‘Phascolomys’ medius) and possible prints of the marsupial lion, Thylacoleo carnifex. The footprints were imprinted within a short time period, demonstrating the association of the taxa present, rather than the time-averaged accumulations usually observed in skeletal fossil deposits. Individual manus and pes prints are distinguishable in some trackways, and in many cases some digital pad morphology is also present. Several parameters traditionally used to differentiate ichnotaxa, including trackway gauge and the degree of print in-turning relative to the midline, are shown to be subject to significant intraspecific variation in marsupials. Sexual dimorphism in the trackway proportions of Diprotodon, and its potential for occurrence in all large bodied, quadrupedal marsupials, is identified here for the first time.     

The Glarus thrust: excursion guide and report of a field trip of the Swiss Tectonic Studies Group (Swiss Geological Society, 14.–16. 09. 2006)

This excursion guide results form a field trip to the Glarus nappe complex organized by the Swiss Tectonic Studies Group in 2006. The aim of the excursion was to discuss old and recent concepts related to the evolution of the Glarus thrust. The major aspects were (i) the interplay between deformation, fluid flow and geochemical alteration, (ii) episodic versus continuous deformation and fluid flow, and (iii) the link between large-scale structures, microstructures, and geochemical aspects. Despite 150 years of research in the Glarus nappe complex and the new results discussed during the excursion, there exist controversies that still are unsolved.     

The challenge of dating early pleistocene fossil teeth by the combined uranium series–electron spin resonance method: the Venta Micena palaeontological site (Orce,Spain)

The palaeontological site of Venta Micena (Orce, Andalusia, Spain) lies in the eastern sector of the Guadix–Baza basin, one of the best documented areas in Europe for Plio‐Pleistocene biostratigraphy. The combination of biochronological and palaeomagnetic results, combined with the radiometric data obtained for Atapuerca Sima del Elefante, indicated that the Venta Micena stratum was formed between the Jaramillo and Olduvai palaeomagnetic events, most likely between 1.22 and 1.77 Ma. Five fossil teeth from two outcrops (sites A and B) were selected to assess the potential of combined uranium series–electron spin resonance (US‐ESR) dating of Early Pleistocene sites. Although the US‐ESR results of the first outcrop showed a large scatter between the three teeth, the mean age of 1.37 ± 0.24 Ma can be considered a reasonable age estimate for Venta Micena. The mean ESR age of 0.62 ± 0.03 Ma obtained for site B seems to be a severe underestimation when compared with the independent age control. This underestimation is attributed to a relative recent U‐mobilization event that led to some U‐leaching. The results show that any ESR age calculations of old samples are extremely sensitive to variations in the measured ²³⁰Th/²³⁴U ratios in dental tissues. Although the results demonstrate that ESR can in principle be applied to Early Pleistocene sites, they also reveal the complexity of dating such old teeth. It is necessary to continue research in several directions, such as study of the behaviour of ESR signals in old teeth and understanding recent U‐mobilization processes, to improve the reliability of the combined US‐ESR dating method applied to Early Pleistocene times, a period for which the number of available numerical dating techniques is very limited. Copyright © 2011 John Wiley & Sons, Ltd.     

Perthite microstructure in magmatic alkali feldspar with oscillatory zoning; Weinsberg Granite, Upper Austria

Feldspar megacrysts from the Weinsberg granite (Moldanubian Zone) show oscillatory zoning with respect to the albite- and orthoclase components. All growth zones show perthitic exsolutions which take the form of bleb- and lens shaped albite-rich precipitates in an orthoclase-rich host. The average sizes and shapes of the precipitates show systematic variation with the integrated bulk compositions of the respective growth zones. The precipitates are abundant and relatively small in growth zones with intermediate bulk composition (Or50Ab41An09 - Or80Ab18An02), and they are less abundant and larger in more orthoclase-rich zones (Or88Ab11An01). Small precipitates have a relatively high aspect ratio, whereas the large precipitates in the potassium-rich zones are more spherical. The relation between microstructure and integrated bulk composition suggests that exsolution and subsequent growth and coarsening occurred by different mechanisms in the respective growth zones. Numerical modeling shows that rapid growth of precipitates over extended periods of time and attainment of relatively large final size is favored, if only few nuclei are formed in an oversaturated host. In contrast, precipitates can grow rapidly only over limited time intervals and remain relatively small, if abundant nuclei are present. During cooling of the oscillatorily zoned alkali-feldspar, exsolution started at relatively high temperatures in growth zones of intermediate integrated bulk composition as compared to exsolution in the more orthoclase-rich growth zones. Irrespective of whether exsolution occurred by spinodal decomposition or by nucleation at relatively high temperatures in the growth zones of intermediate integrated bulk composition, it produced abundant nuclei and resulted in relatively small precipitates. In contrast, comparatively few nuclei were formed in the orthoclase-rich growth zones resulting in large precipitates. The Na/K partitioning between precipitates and the host is independent of the integrated bulk composition of the respective growth zone reflecting re-equilibration during cooling down to relatively low temperatures (<?400°C). The shape of the precipitates probably has evolved from an initially lamellar or spindle-like geometry with high aspect ratio to more isometric, spheroidal shapes during precipitate growth and coarsening. Host/precipitate interfaces served as fluid pathways during late stage deuteric alteration.     

The competing models for the origin and internal evolution of granitic pegmatites in the light of melt and fluid inclusion research

In this paper we discuss the main petrogenetic models for granitic pegmatites and how these models have evolved over time. We suggest that the present state of knowledge requires that some aspects of these models to be modified, or absorbed into newer ones. Pegmatite formation and internal evolution have long supposed the need for highly water- and flux-enriched magmas to explain the differences between pegmatites and other intrusives of similar major element composition. Compositions and textural characteristics of fluid and melt inclusions in pegmatite minerals provide strong evidence for such magmas. Furthermore, we show that melt inclusion research has increased the number of potential flux components, which may include H₂O, OH^?, CO₂, HCO ₃ ^? , CO ₃ ^2? , SO ₄ ^2? , PO ₄ ^3? , H₃BO₃, F , and Cl, as well as the elements Li, Na, K, Rb, Cs, and Be, herein described as melt structure modifiers. In this paper we emphasize that the combined effect which these components have on the properties of pegmatite melts is difficult to deduce from experimental studies using only a limited number of these components. The combination and the amount of the different magmatic species, together with differences in the source region, and variations in pressure and temperature cause the great diversity of the pegmatites observed. Some volatile species, such as CO ₃ ^2? and alkalis, have the capacity to increase the solubility of H₂O in silicate melt to an extraordinary degree, to the extent that melt-melt-fluid immiscibility becomes inevitable. It is our view that the formation of pegmatites is connected with the complex interplay of many factors.