Automatic correction and simplification of geological maps and cross-sections for numerical simulations

Incorporating prior geological knowledge in geophysical process models often meets practical meshing challenges and raises the question of how much detail is to be included in the geometric model. We introduce a strategy to automatically repair and simplify geological maps, geological cross-sections and the associated meshes while preserving elementary consistency rules. To identify features breaking validity and/or the thin features potentially problematic when generating a mesh, we associate an exclusion zone with each model feature (horizon, fault). When these zones overlap, both the connectivity and the geometry of the geological layers are automatically modified. The output model enforces specific practical quality criteria on the model topology and geometry that facilitates the generation of a mesh with lower bounds on minimum angles and minimum local entity sizes. Our strategy is demonstrated on an invalid geological cross-section from a real-case study in the Lorraine coal basin. We further explore the impacts of the model modifications on wave propagation simulation. We show that the differences on the seismograms due to model simplifications are relatively small if the magnitude of simplifications is adapted to the physical problem parameters.     

H2O Content Measurement in Phengite by Secondary Ion Mass Spectrometry: A New Set of Reference Materials

Five new natural white mica reference materials (RMs) were developed for in situ H₂O content analyses by secondary ion mass spectrometry at the SwissSIMS laboratory of Lausanne University, Switzerland. The white mica reference materials cover a large part of the natural muscovite–phengite compositional range and are therefore suitable as reference materials for the analysis of natural rocks as well as individual minerals. The independent H₂O content of the reference materials UNIL_WM1 to UNIL_WM5 was obtained by thermal conversion elemental analyser and corresponds to 4.35 ± 0.02, 4.33 ± 0.03, 4.30 ± 0.07, 4.50 ± 0.02 and 4.42 ± 0.11 (% m/m, ± 1s), respectively. SIMS determinations of H₂O content revealed a matrix effect correlated to the FeO content of white mica. The compositional range in FeO of the reference materials that were calibrated for H₂O determination is from 1.13% to 3.67% m/m. No crystallographic orientation dependency was observed at the level of homogeneity of these reference materials. An analytical precision of 0.02% to 0.08% m/m (1SE) is expected for the final uncertainty on measurements of unknown white micas in natural samples.     

The Planeterrella,a pedagogic experiment in planetology and plasma physics

We present here a plasma physics experiment which makes it possible to simulate, in a naive yet useful way, the formation of polar lights. It involves shooting electrons at a magnetized sphere placed in a vacuum chamber. This experiment, inspired by K. Birkeland’s Terrella, built at the turn of 19th century, allows the visualization of very many geophysical and astrophysical situations. Although delicate, it is feasible at undergraduate level.     

Reconstruction of Channelized Systems Through a Conditioned Reverse Migration Method

Geological heterogeneities directly control underground flow. In channelized sedimentary environments, their determination is often underconstrained: it may be possible to observe the most recent channel path and the abandoned meanders on seismic or satellite images, but smaller-scale structures are generally below image resolution. In this paper, reconstruction of channelized systems is proposed with a stochastic inverse simulation reproducing the reverse migration of the system. Maps of the recent trajectories of the Mississippi river were studied to define appropriate relationships between simulation parameters. Measurements of curvature and migration vectors showed (i) no significant correlation between curvature and migration offset and (ii) correlation trends of downstream and lateral migration offsets versus the curvature at half-meander scale. The proposed reverse migration method uses these trends to build possible paleo-trajectories of the river starting from the last stage of the sequence observed from present-day (satellite or seismic) data. As abandoned meanders provide clues about the paleo-locations of the river, they are integrated time step by time step during the reverse simulation process. We applied the method to a satellite image of a fluvial system. Each of the different resulting geometries of the system honored most of the available observations and presented meandering patterns similar to the observed ones.     

The importance of chemical data for the study of interstellar shocks

Shock waves perturb the chemical state of the interstellar gas. We consider the effects of C-shocks on the composition of molecular clouds, for a range of values of the pre-shock gas density and magnetic induction. The time required to re-establish equilibrium in the post-shock gas depends on the initial conditions and can become very large. The significance of the two known chemical phases of dark clouds and of bistability is considered in this context.     

Structural Setting of the Neoarchean Terrains in the Commonwealth Bay Area (143-145°E), Terre Adélie Craton, East Antarctica

Geological maps of East Commonwealth Bay Unit (ECB), (Terre Adélie and Georges V Land, Antarctica) are presented with a summary of the main structural and metamorphic data for the region. The ECB unit was developed during Neoarchean_Paleoproterozoic event (at 2.5-2.42 Ga), with (i) granulite metamorphism at 9±1.5 kbar and 800±50°C in the lower crust section and amphibolite metamorphism (P=5 kbar, T=750°C) at the upper crustal levels; (ii) the lower crustal granulites were uplifted, and suffered local partial melting and retrogression to the amphibolite facies at 550±50°C_5 kbar. Granulites were extruded in the core of a crustal-scale anticlinal fold, but retrogressed only on the rims of the anticline. Crustal-scale folding, along with other structural features resulted from intense NE-SW shortening that prevailed during the Neoarchean orogenic cycle. Strike-slip and extensional motions were only minor components in that process; (iii) top-to-the-East thrusting and nappe piling had (at least locally) occurred under lower amphibolite to greenschist facies conditions. Finally, it seems that (iv) the Paleoproterozoic 1.7 Ga structural imprint may have only affected the rims of the Archean units. The tectonic context observed in the 1.7 Ga Cape Hunter phyllites features mainly an E-W shortening component and vertical extrusion. The eastern (Mertz) and western (Port Martin) parts of the Archean block were reactivated by localized dextral shearing.     

Quaternary faulting in the Upper Rhine Graben revealed by high-resolution multi-channel reflection seismic

A high-resolution multichannel seismic reflection river profiling campaign was completed in July 2002 in the southern Upper Rhine graben (URG), along the River Rhine. Preliminary results show apparent Quaternary vertical slip rates, on intra-graben faults that are relatively slow, of the order of a few thousandths to a few hundredths of mm/yr. Moreover, kinematical data from the Ludwigshafen area show decreasing vertical slip rates since the Middle Pleistocene and/or a migration of tectonic activity. While still preliminary, these data show inhomogeneous and relatively slow tectonic activity in the URG that could probably not alone have shaped the Quaternary graben morphology. To cite this article: G. Bertrand et al., C. R. Geoscience 338 (2006).     

Aquifer-yield continuum as a guide and typology for science-based groundwater management

Groundwater availability is at the core of hydrogeology as a discipline and, simultaneously, the concept is the source of ambiguity for management and policy. Aquifer yield has undergone multiple definitions resulting in a range of scientific methods to calculate and model availability reflecting the complexity of combined scientific, management, policy, and stakeholder processes. The concept of an aquifer-yield continuum provides an approach to classify groundwater yields along a spectrum, from non-use through permissive sustained, sustainable, maximum sustained, safe, permissive mining to maximum mining yields, that builds on existing literature. Additionally, the aquifer-yield continuum provides a systems view of groundwater availability to integrate physical and social aspects in assessing management options across aquifer settings. Operational yield describes the candidate solutions for operational or technical implementation of policy, often relating to a consensus yield that incorporates human dimensions through participatory or adaptive governance processes. The concepts of operational and consensus yield address both the social and the technical nature of science-based groundwater management and governance.     

Finite Difference Implicit Structural Modeling of Geological Structures

We introduce a new method for implicit structural modeling. The main developments in this paper are the new regularization operators we propose by extending inherent properties of the classic one-dimensional discrete second derivative operator to higher dimensions. The proposed regularization operators discretize naturally on the Cartesian grid using finite differences, owing to the highly symmetric nature of the Cartesian grid. Furthermore, the proposed regularization operators do not require any special treatment on boundary nodes, and their generalization to higher dimensions is straightforward. As a result, the proposed method has the advantage of being simple to implement. Numerical examples show that the proposed method is robust and numerically efficient.