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.     

Characterization of a major refractory component of marine dissolved organic matter

Refractory carboxyl-rich alicyclic molecules (CRAM) are characterized in marine dissolved organic matter (DOM) using nuclear magnetic resonance spectroscopy and ultrahigh resolution mass spectrometry. CRAM are distributed throughout the water column and are the most abundant components of deep ocean DOM ever characterized. CRAM are comprised of a complex mixture of carboxylated and fused alicyclic structures with a carboxyl-C:aliphatic-C ratio of 1:2 to 1:7. CRAM are expected to constitute a strong ligand for metal binding, and multiple coordination across cations could promote aggregation and marine gel formation thereby affecting CRAM reactivity and the bioavailability of nutrients and trace metals. It appears CRAM are ultimately derived from biomolecules with structural similarities to sterols and hopanoids. The occurrence of CRAM in freshwater and terrestrial environments seems likely, considering the global distribution of biomolecules and the similarities of biogeochemical processes among environments.     

Benchmarking spatial joins a la carte

Spatial joins are join operations that involve spatial data types and operators. Spatial access methods are often used to speed up the computation of spatial joins. This paper addresses the issue of benchmarking spatial join operations. For this purpose, we first present a WWW-based benchmark generator to produce sets of rectangles. Using a Web browser, experimenters can specify the number of rectangles in a sample, as well as the statistical distributions of their sizes, shapes, and locations. Second, using the generator and a well-defined set of statistical models we define several tests to compare the performance of three spatial join algorithms: nested loop, scan-and-index, and synchronized tree traversal. We also added two real-life data sets from the Sequoia 2000 storage benchmark. Our results show that the relative performance of the different techniques mainly depends on the selectivity factor of the join predicate. All of the statistical models and algorithms are available on the Web, which allows for easy verification and modification of our experiments.     

Metamorphism at the northwest contact of the Stanhope Pluton,Quebec Appalachians: Mineral equilibria in interbedded pelite and calc-schist

The interbedding of pelite and calc-schist in part of the contact aureole of a Devonian (biotite K-Ar age 346±4 Ma) granite pluton that straddles the Quebec-Vermont border offers an opportunity to compare metamorphic conditions prevailing in both rock types, and to test for internal consistency among several different geothermometers and geobarometers.Microprobe analyses and recent thermodynamic data are used with simple activity-composition models to estimate P-T-X _fluid conditions near the sillimanite zone. Calculations yield the following results: P=2 to 3 kb, and T =400° to 600° C in the inner one kilometre of the aureole; was around 0.8. Estimates are based on the calc-silicate isobaric invariant assemblages tremolite-K-feldspar-plagioclase-clinozoisite-phlogopite-calcite-quartz and diopside-tremolite-plagiolase-clinozoisite-calcite-quartz, and the paragonite-andalusite-sillimanite-albite-quartz equilibrium. The solid-solid reaction phlogopite+2 diopside + 4 quartz=tremolite+K-feldspar in the calc-schist (combined with the andalusite-sillimanite equilibrium), and phase relations in the granite yield apparently inconsistent results.The implied 6 to 10 km of cover at the time of intrusion may have been provided by subsequently eroded thrust sheets.     

The sandstone-hosted stratiform copper mineralization at Mwitapile and its relation to the mineralization at Lufukwe, Lufilian foreland, Democratic Republic of Congo

The Lufilian foreland is a triangular-shaped area located in the SE of the Democratic Republic of Congo and to the NE of the Lufilian arc, which hosts the well-known Central African Copperbelt. The Lufilian foreland recently became an interesting area with several vein-type (e.g., Dikulushi) and stratiform (e.g., Lufukwe and Mwitapile) copper occurrences. The Lufilian foreland stratiform Cu mineralization is, to date, observed in sandstone rock units belonging to the Nguba and Kundelungu Groups (Katanga Supergroup).The Mwitapile sandstone-hosted stratiform Cu prospect is located in the north eastern part of the Lufilian foreland. The host rock for the Cu mineralization is the Sonta Sandstone of the Ngule Subgroup (Kundelungu Group). A combined remote sensing, petrographic and fluid inclusion microthermometric analysis was performed at Mwitapile and compared with similar analysis previously carried out at Lufukwe to present a metallogenic model for the Mwitapile- and Lufukwe-type stratiform copper deposits. Interpretation of ETM+ satellite images for the Mwitapile prospect and the surrounding areas indicate the absence of NE–SW or ENE–WSW faults, similar to those observed controlling the mineralization at Lufukwe. Faults with these orientations are, however, present to the NW, W, SW and E of the Mwitapile prospect. At Mwitapile, the Sonta Sandstone host rock is intensely compacted, arkosic to calcareous with high silica cementation (first generation of authigenic quartz overgrowths). In the Sonta Sandstone, feldspar and calcite are present in disseminated, banded and nodular forms. Intense dissolution of these minerals caused the presence of disseminated rectangular, pipe-like and nodular dissolution cavities. Sulfide mineralization is mainly concentrated in these cavities. The hypogene sulfide minerals consist of two generations of pyrite, chalcopyrite, bornite and chalcocite, separated by a second generation of authigenic quartz overgrowth. The hypogene sulfide minerals are replaced by supergene digenite and covellite. Fluid inclusion microthermometry on the first authigenic quartz phase indicates silica precipitation from an H₂O–NaCl–CaCl₂ fluid with a minimum temperature between 111 and 182 °C and a salinity between 22.0 and 25.5 wt.% CaCl₂ equiv. Microthermometry on the second authigenic quartz overgrowths and in secondary trails related to the mineralization indicate that the mineralizing fluid is characterized by variable temperatures (Th = 120 to 280 °C) and salinities (2.4 to 19.8 wt.% NaCl equiv.) and by a general trend of increasing temperatures with increasing salinities.Comparison between Mwitapile and Lufukwe indicates that the stratiform Cu mineralization in the two deposits is controlled by similar sedimentary, diagenetic and structural factors and likely formed from a similar mineralizing fluid. A post-orogenic timing is proposed for the mineralization in both deposits. The main mineralization controlling factors are grain size, clay and pyrobitumen content, the amount and degree of feldspar and/or calcite dissolution and the presence of NE–SW to ENE–WSW faults. The data support a post-orogenic fluid-mixing model for the Mwitapile- and Lufukwe-type sandstone-hosted stratiform Cu deposits, in which the mineralization is related to the mixing between a Cu-rich hydrothermal fluid, with a temperature up to 280 °C and a maximum salinity of 19.8 wt.% NaCl equiv., with a colder low salinity reducing fluid present in the sandstone host rock. The mineralizing fluid likely migrated upwards to the sandstone source rocks along NE–SW to ENE–WSW orientated faults. At Lufukwe, the highest copper grades at surface outcrops and boreholes were found along and near to these faults. At Mwitapile, where such faults are 2 to 3 km away, the Cu grades are much lower than at Lufukwe. Copper precipitation was possibly promoted by reduction from pre-existing hydrocarbons and non-copper sulfides and by the decrease in fluid salinity and temperature during mixing. Based on this research, new Cu prospects were proposed at Lufukwe and Mwitapile and a set of recommendations for further Cu exploration in the Lufilian foreland is presented.     

Cretaceous length of day perturbation by mantle avalanche

These last 10 years, numerical models of mantle convection have emphasized the role of the 670 km endothermic phase change in generating avalanches that trigger catastrophic mass transfers between upper and lower mantle. On the other hand, scientists have emphasized the concomitance of large-scale worldwide geophysical and tectonic events, which could find their deep thermal roots in the huge mass transfers induced by the avalanches. In particular, the paleontological records show two periods of length of day (l.o.d.) shortening between 420 and 360, and 200 and 80 Myr BP. This last event is synchronous with a strong true polar wander and a global warming of the upper mantle. In order to study the potential effects of the avalanche on the main component of the Earth’s rotation, the Liouville equation has been solved and the l.o.d. evolution has been calculated from the perturbations of the inertia tensor. The results show that the inertia tensor of the Earth’s is mainly sensitive to the global transfers through the 670 km discontinuity. The l.o.d. perturbations will be synchronous with the global thermal effects of the avalanche. These theoretical results allow proposing a self-consistent physical mechanism to explain periods of the Earth’s rotation acceleration. Within this context, the l.o.d. shortening during the Cenozoic and Cretaceous brings one more clue to the possible participation of a mantle avalanche in generating the concomitant large scale events which have occurred during this very particular period of the Earth’s history.     

Global thermal and dynamical perturbations due to Cretaceous mantle avalanche

The numerical models of mantle convection agree to depict avalanches behaviour according to the level of endothermicity of the spinel → perovskite phase change. Their potential effects on the global thermal and dynamical states of the mantle have been computed thanks to a numerical code, which takes into account both the 400-km exothermic and the 660-km endothermic phase changes. The cycle followed by the avalanches is: local layering, destabilization of the 660-km thermal layer, travelling and spreading on the core, and reappearing of the local layering. Therefore, mantle convection is characterized by quiet periods of partial layering embedded in catastrophic events. During the avalanche, the amplitude of the surface velocity is multiplied by two, which would imply an enhanced plate tectonic and ridge activities. The global thermal effects of the avalanche are compatible with a high mantle temperature and an acceleration of Earth's rotation during the Cretaceous. They also offer a coherent explanation to locate the origin of mantle plumes both within the CMB and just below the transition zone.     

Coupling mechanisms in double sandbar systems. Part 2: impact on alongshore variability of inner‐bar rip channels

Double sandbar systems are common morphological features along sandy, wave‐dominated, micro‐ to meso‐tidal coastlines. In the companion paper, we demonstrated how various alongshore inner‐bar rip‐channel patterns can develop through morphological coupling to an alongshore‐variable outer bar. The simulated coupling patterns are, however, scarcely observed in the field. Instead, inner‐bar rip channels more often possess remarkably smaller and more variable alongshore length scales, suggesting that coupling mechanisms do not play a substantial role in the overall double‐sandbar dynamics. Here we use a numerical model to show that the relative importance of self‐organization and morphological coupling changes in favour of the latter with an increase in waterdepth variability along the outer‐bar crest. Furthermore, we find that the typical alongshore variability in inner‐bar rip‐channel scale is indicative of a mixture of self‐organization and morphological coupling rather than self‐organization alone. Morphological coupling may thus be more important to understanding and predicting the evolution of inner‐bar rip channels than previously envisaged. Copyright © 2010 John Wiley and Sons, Ltd.