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.     

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.     

The effect of long-term low-temperature exposure on apatite fission track stability: A natural annealing experiment in the deep ocean

Since studies on deep-sea cores were carried out in the early 1990s it has been known that ambient temperature may have a marked affect on apatite fission track annealing. Due to sluggish annealing kinetics, this effect cannot be quantified by laboratory annealing experiments. The unknown amount of low-temperature annealing remains one of the main uncertainties for extracting thermal histories from fission track data, particularly for samples which experienced slow cooling in shallow crustal levels. To further elucidate these uncertainties, we studied volcanogenic sediments from five deep-sea drill cores, that were exposed to maximum temperatures between ∼10° and 70 °C over geological time scales of ∼15-120 Ma. Mean track lengths (MTL) and etch pit diameters (Dpar) of all samples were measured, and the chemical composition of each grain analysed for age and track length measurements was determined by electron microprobe analysis. Thermal histories of the sampled sites were independently reconstructed, based on vitrinite reflectance measurements and/or 1D numerical modelling. These reconstructions were used to test the most widely used annealing models for their ability to predict low-temperature annealing. Our results show that long-term exposure to temperatures below the temperature range of the nominal apatite fission track partial annealing zone results in track shortening ranging between 4 and 11%. Both chlorine content and Dpar values explain the downhole annealing patterns equally well. Low chlorine apatite from one drill core revealed a systematic relation between Si-content and Dpar value. The question whether Si-substitution in apatite has direct and systematic effects on annealing properties however, cannot be addressed by our data. For samples, which remained at temperatures <30 °C, and which are low in chlorine, the Laslett et al. [Laslett G., Green P., Duddy I. and Gleadow A. (1987) Thermal annealing of fission tracks in apatite. Chem. Geol. 65, 1-13] annealing model predicts MTL up to 0.6 μm longer than those actually measured, whereas for apatites with intermediate to high chlorine content, which experienced temperatures >30 °C, the predictions of the Laslett et al. (1987) model agree with the measured MTL data within error levels. With few exceptions, predictions by the Ketcham et al. [Ketcham R., Donelick R. and Carlson W. (1999) Variability of apatite fission-track annealing kinetics. III: Extrapolation to geological time scales. Am. Mineral. 84/9, 1235-1255] annealing model are consistent with the measured data for samples which remained at temperatures below ∼30 °C. For samples which experienced maximum temperatures between ∼30 and 70 °C, and which are rich in chlorine, the Ketcham et al. (1999) model overestimates track annealing.     

Facies architecture of an isolated carbonate platform: tracing the cycles of the Latemàr (Middle Triassic, northern Italy)

The 720-m-thick succession of the Middle Triassic Latemàr Massif (Dolomites, Italy) was used to reconstruct the lagoonal facies architecture of a small atoll-like carbonate platform. Facies analysis of the lagoonal sediments yields a bathymetric interpretation of the lateral facies variations, which reflect a syndepositional palaeorelief. Based on tracing of lagoonal flooding surfaces, the metre-scale shallowing-upward cycles are interpreted to be of allocyclic origin. Short-term sea-level changes led to subaerial exposure of wide parts of the marginal zone, resulting in the development of a tepee belt of varying width. Occasional emergence of the entire lagoon produced lagoon-wide decimetre-thick red exposure horizons. The supratidal tepee belt in the backreef area represented the zone of maximum elevation, which circumscribed the sub- to peritidal lagoonal interior during most of the platform's development. This tepee rim, the subtidal reef and a sub- to peritidal transition zone in between stabilized the platform margin. The asymmetric width of facies belts within individual metre-scale cycles was caused by redistribution processes that reflect palaeowinds and storm paths from the present-day south and west. The overall succession shows stratigraphic changes on a scale of tens of metres from a basal subtidal unit, overlain by three tepee-rich intervals, separated by tepee-poor units composed of subtidal to peritidal facies. This stacking pattern reflects two third-order sequences during the late Anisian to early middle Ladinian.     

Sublimation pressure and sublimation enthalpy of H2O ice Ih between 0 and 273.16 K

The sublimation of water in the region of extremely low pressures and temperatures meets increasing interest for the exploration of icy cosmic bodies. At temperatures below 130 K, the shape of the sublimation curve of H₂O is not known; neither experimental data nor theoretical treatments exist for this region. Based upon theoretical upper and lower bounds for the heat capacity of water vapor in this range, a narrow region in the pressure-temperature diagram is identified which must necessarily enclose the sublimation curve down to virtually any lower pressures and temperatures. Within this region, an approximate sublimation curve is computed from the 2006 Gibbs potential of ice Ih, using available heat capacity data points of water vapor between 10 and 130 K. The theoretical zero-point limiting law of the sublimation pressure is derived. Valid between 20 and 273.16 K, correlation equations for the sublimation enthalpy and the sublimation pressure are fitted to the computed data. All quantities are expressed in the 1990 temperature scale ITS-90. Under cosmic conditions, our results suggest that the sublimation of ice is unlikely below 50 K and impossible below 23 K.     

Dynamics of coastline changes in Mexico

Journal of Geographical Sciences - The goal of our work was to locate and quantify changes that occurred in 66% of the Mexican coastline, based on four land cover maps generated by the Mexican...     

Interpretation of shallow refraction seismic data by reflection/refraction tomography

A new algorithm for tomographic inversion of traveltimes of reflected and refracted seismic waves is developed. The inversion gives interface configurations and velocity distributions in layers. The important features of the algorithm are: (a) the inclusion of shot time delays in the list of unknown parameters; (b) the regularization is applied in such a way that the most probable model is characterized by the similarity of neighbouring interfaces. As the problem under consideration is non-linear, several iterations are necessary in order to obtain the final model. In the case of a very inexact initial model, a 'layer-by-layer' inversion strategy is recommended as a first inversion step. The inversion program is supplied with a user interface, thanks to which one can: (a) pick interactively and identify seismic traveltimes; (b) build and edit depth/velocity models; and (c) display calculated traveltime curves and compare them with picked traveltimes as well as with the original seismic sections. The efficiency of the inversion software developed is illustrated by a numerical example and a field example in which shallow seismic data are considered. Application to wide-aperture reflection/refraction profiling (WARRP) data is also possible.     

Thermal evolution of an extensional detachment as constrained by organic metamorphic data and thermal modeling: Graz Paleozoic Nappe Complex (Eastern Alps)

Following Early Cretaceous nappe stacking, the Eastern Alps were affected by late-orogenic extension during the Late Cretaceous. In the eastern segment of this range, a Late Cretaceous detachment separates a very low- to low-grade metamorphic cover (Graz Paleozoic Nappe Complex, GPNC) above a low- to high-grade metamorphic basement. Synchronously, the Kainach Gosau Basin (KGB) collapsed and subsided on top of the section.Metamorphism of organic material within this section has been investigated using vitrinite reflectance data and Raman spectra of extracted carbonaceous material. In the southern part of the GPNC, vitrinite reflectance indicates a decrease in organic maturity towards the stratigraphic youngest unit. The remaining part of the GPNC is characterized by an aureole of elevated vitrinite reflectance values and Raman R2 ratios that parallels the margins of the GPNC. Vitrinite reflectance in the KGB shows a steep coalification gradient and increases significantly towards the western basin margin. The observed stratigraphic trend in the southern GPNC is a result of deep Paleozoic to Early Cretaceous burial. This maturity pattern was overprinted along the margins by advective heat and convective fluids during Late Cretaceous to Paleogene exhumation of basement rocks.During shearing, the fault zone was heated up to ca. 500 °C. This overprint is explained by a two-dimensional thermal model with a ramp-flat fault geometry and a slip rate of 1 to 1.5 cm/year during 5 Ma fault movement. The collapse basin above the detachment subsided in a thermal regime which was characterized by relaxing isotherms.     

Meteoritic material at four Canadian impact craters

Eleven impact melt and 6 basement rock samples from 4 craters were analyzed by neutron activation for Au, Co, Cr, Fe, Ge, Ir, Ni, Os, Pd, Re and Se. Wanapitei Lake, Ontario: the impact melts show uniform enrichments corresponding to 1–2% C1-chondrite material. Interelement ratios ($\text{Co}\text{Cr}$, $\text{Ni}\text{Cr}$, $\text{Ni}\text{Ir}$) suggest that the impacting body was a Cl-, C2-, or LL-chondrite. Nicholson Lake, North West Territory: Ni, Cr and Co are distinctly more enriched than Ir and Au which tentatively suggests an olivine-rich achondrite (nakhlite or ureilite). Gow Lake, Saskatchewan and Mistastin, Labrador: small enrichments in Ir and Ni; both the low $\text{Ir}\text{Ni}$ ratios and low Cr content suggest iron meteorites, but the signals are too weak for conclusive identification.A tentative comparison of meteoritic signatures at 10 large, ≥4km craters and their presumed celestial counterparts (13 Apollo and Amor asteroids) shows more irons and achondrites among known projectile types, and a preponderance of S-type objects, having no known meteoritic equivalent, among asteroids. It is not yet clear that these differences are significant, in view of the tentative nature of the crater identifications (achondrites in particular), and the limited statistics.