A roadmap for a dedicated Earth Science Grid platform

Due to its intensive data processing and highly distributed organization, the multidisciplinary Earth Science applications community is uniquely positioned for the uptake and exploitation of Grid technologies. Currently Enabling Grids for E-sciencE, and other large Grid infrastructures are already deployed and capable of operational services. So far however, the adoption and exploitation of Grid technology throughout the Earth Science community has been slower than expected. The Dissemination and Exploitation of GRids in Earth sciencE project, proposed by the European Commission to assist and accelerate this process in a number of different ways, had between its main goals the creation of a roadmap towards Earth Science Grid platform. This paper presents the resulting roadmap.     

An adaptive finite element model for wellbore breakout simulation

A finite element formulation is proposed and implemented for analysing the stability of excavated wells using the DiMaggio-Sandler constitutive elastoplastic model with a typical carbonate reservoir configuration. The quality of the finite element approximation is ensured by applying smooth curved elements adapted to the wellbore geometry, and h − p adaptive finite element meshes in the plastic zone. General purpose procedures are defined to transfer the elastoplastic deformation history to newly created integration points. A breakout damage criterion is proposed based on the second invariant of the deviatoric plastic deformation tensor. This damage criterion is used to apply a mesh movement algorithm to represent material collapse. The automatic successive application of the breakout damage criterion results in elliptical realistically looking geometries obtained in experiments reported in the literature.     

Surface chemistry and reactivity of biogenic silica

The surface chemistry of cultured diatoms was compared to that of biosiliceous material in Southern Ocean sediments, using potentiometric titrations and aluminum adsorption experiments. Aerosil 200, a well-studied synthetic amorphous silica, served as reference solid. Surface charge development and aluminum adsorption on cultured diatom shells were comparable to Aerosil. The surface chemical properties of biosiliceous material buried to depths of 15-25 cm in Southern Ocean sediments, however, deviated markedly from those of the cultured diatoms. In pH range 4-8.5, the surface charge density was systematically lower for biogenic silica from the sediments. In addition, the aluminum adsorption edge on the biosiliceous sediments was shifted to higher pHs by about 0.4 units. The results indicate that ionizable surface silanol groups on diagenetically altered diatom shells are less abundant and, possibly, less acidic than on freshly cultured diatoms. The observed differences in surface chemical structure are consistent with the progressive loss of reactivity, or aging, of biogenic silica which promotes its preservation in sediments.     

Modeling of in-situ crystallization processes in the Permian mafic layered intrusion of Mont Collon (Dent Blanche nappe, western Alps)

The Mont Collon mafic complex is one of the best preserved examples of the Early Permian magmatism in the Central Alps, related to the intra-continental collapse of the Variscan belt. It mostly consists (> 95 vol.%) of ol + hy-normative plagioclase-wehrlites, olivine- and cpx-gabbros with cumulitic structures, crosscut by acid dikes. Pegmatitic gabbros, troctolites and anorthosites outcrop locally. A well-preserved cumulative sequence is exposed in the Dents de Bertol area (center of intrusion). P–T calculations indicate that this layered magma chamber emplaced at mid-crustal levels at about 0.5 GPa and 1100 °C. The Mont Collon cumulitic rocks record little magmatic differentiation, as illustrated by the restricted range of clinopyroxene mg-number (Mg#_cpx = 83–89). Whole-rock incompatible trace-element contents (e.g. Nb, Zr, Ba) vary largely and without correlation with major-element composition. These features are characteristic of an in-situ crystallization process with variable amounts of interstitial liquid L trapped between the cumulus mineral phases. LA-ICPMS measurements show that trace-element distribution in the latter is homogeneous, pointing to subsolidus re-equilibration between crystals and interstitial melts. A quantitative modeling based on Langmuir's in-situ crystallization equation successfully duplicated the REE concentrations in cumulitic minerals of all rock facies of the intrusion. The calculated amounts of interstitial liquid L vary between 0 and 35% for degrees of differentiation F of 0 to 20%, relative to the least evolved facies of the intrusion. L values are well correlated with the modal proportions of interstitial amphibole and whole-rock incompatible trace-element concentrations (e.g. Zr, Nb) of the tested samples. However, the in-situ crystallization model reaches its limitations with rock containing high modal content of REE-bearing minerals (i.e. zircon), such as pegmatitic gabbros. Dikes of anorthositic composition, locally crosscutting the layered lithologies, evidence that the Mont Collon rocks evolved in open system with mixing of intercumulus liquids of different origins and possibly contrasting compositions. The proposed model is not able to resolve these complex open systems, but migrating liquids could be partly responsible for the observed dispersion of points in some correlation diagrams. Absence of significant differentiation with recurrent lithologies in the cumulitic pile of Dents de Bertol points to an efficiently convective magma chamber, with possible periodic replenishment.     

Accuracy and efficiency of time integration methods for 1D diffusive wave equation

The diffusive wave approximation of the Saint-Venant equations is commonly used in hydrological models to describe surface flow processes. Numerous numerical approaches can be used to solve this highly nonlinear equation. Nonlinear time integration schemes—also called methods of lines (MOL)—were proven very efficient to solve other nonlinear problems in geosciences but were never considered to deal with surface flow modeling with the diffusive wave equation. In this paper, we study the relative performance of different time and space integration schemes by comparing the results obtained with classical approaches and with nonlinear time integration approaches. The results show that (i) the integration method with a higher order in space shows high accuracy regarding an integrated indicator such as the global mass balance error but is less accurate regarding local indicators, and (ii) nonlinear time integration techniques perform better than classical ones. Overall, it seems that integration techniques combining nonlinear time integration and a low spatial order need to be considered when developing hydrological modeling tools owing to their simplicity of implementation and very good performance.     

Lateglacial palaeoenvironments and palaeoclimates from Conty and Houdancourt, northern France, reconstructed from Beetle remains

Two Lateglacial insect sequences are described (a) from Conty, in the Selle valley and (b) from Houdancourt in the Oise valley, both in northern France. These investigations are part of a multidisciplinary investigation involving archaeology, sedimentology, geomorphology, geochronology and palaeontology (pollen, macroscopic plant remains, vertebrates, molluscs and insect fossils). The sequences of insect assemblages date from the latter part of the Bølling to the end of the Allerød periods. Environmental analysis of these faunas shows that rivers with riffles and pools meandered across flood plains. The river was extensively fringed with reedy vegetation. The only trees growing close to the river were of Salix and/or Populus. All the insect assemblages indicate that the thermal climates during the Bølling and Allerød periods were similar to one another and during both periods were very nearly as warm as that of the present day. No insect fossils were recovered from the sediments attributed to the Older Dryas interval though other evidence from these sites suggests that this event was decidedly colder than those immediately preceding and succeeding it. Comparisons are made between Lateglacial climatic patterns in northern France with those elsewhere in Europe.     

Seismic vulnerability assessment for Montreal

In Canada, Montreal is the second city with the highest seismic risk. This is due to its relatively high seismic hazard, old infrastructures and high population density. The region is characterised by moderate seismic activity with no recent record of a major earthquake. The lack of historical strong ground motion records for the region contributes to large uncertainties in the estimation of hazards. Among the sources of uncertainty, the attenuation function is the main contributor and its effect on estimates of risks is investigated. Epistemic uncertainty was considered by obtaining damage estimates for three attenuation functions that were developed for Eastern North America. The results indicate that loss estimates are highly sensitive to the choice of the attenuation function and suggest that epistemic uncertainty should be considered both for the definition of the hazard function and in loss estimation methodologies. Seismic loss estimates are performed for a 2% in 50 years seismic threat, which corresponds to the design level earthquake in the national building code of Canada, using HAZUS-MH4 for the Montreal region over 522 census tracts. The study estimated that for the average scenario roughly 5% of the building stock would be damaged with direct economic losses evaluated at 1.4 billion dollars for such a scenario. The maximum number of casualties would result in approximately 500 people being injured or dead at a calculated time of occurrence of 2?pm.     

Modeling the impact of microbial activity on redox dynamics in porous media

The present study investigates the interaction between microbial growth and activity and the redox dynamics in natural porous media. The impact the transport regime has on this interaction is also addressed. Expressions for microbial growth are incorporated into a geochemical reaction network linking redox reaction rates to the activity of the microorganisms. A flexible simulation environment, the Biogeochemical Reaction Network Simulator (BRNS) is used for this purpose. Two reactive transport applications relevant to fields of contaminant hydrology and early diagenesis are simulated with the BRNS. Model results are evaluated based on a comparison with comprehensive datasets on the biodegradation of lactate in a sand column experiment and on the distribution of redox-sensitive chemical species in marine sediments of the Skagerrak, Denmark. It is shown that, despite quite different transport regimes, the geomicrobiological model performs equally well in the reproduction of measured chemical species distribution for both applications. This result emphasizes the broad applicability of the proposed approach. Our simulations support that the competitive behavior between various microbial groups is a process controlling the development of redox stratified environments. Furthermore, it is also shown that the transport regime is a key controlling factor for the degree of spatial correlation between microbial biomass distributions and redox reaction rates. Although all our simulations yield a pronounced stratification of the redox processes in the system, the biomass distribution is related to the associated reaction rates only in case of the advection controlled column experiment. In the early diagenetic application, mixing due to bioturbation is the dominant transport process for particulate matter, hence leading to fairly homogeneous distribution of bacterial biomasses which are unrelated to the spatial distribution of redox reaction rates. This homogeneous biomass distribution combined with the 1G carbon degradation model approach might explain why the steady state concentration profiles in such systems can be reproduced by diagenetic models without explicit representation of microbial growth.