A fast flux tube-based method for solute-transport simulation

A new method to calculate the transport of dissolved species in aquifers is presented. This approach is an extension of the stream tubes which are used for flow computation. The flux tubes defined here are conservative for solutes, but not for water mass. The flux tubes are first defined in a general domain and then calculated in a two-dimensional uniform flow field. The tubes?? computation is based on a parametric solution. The method is extended further in order to deal with heterogeneous media. A particle-tracking algorithm is used where the deviation of the flux-tube boundaries due to dispersion is accounted for. The approximate solution obtained by this approach is compared to classical numerical solutions given by a finite difference code (RT3D) and a finite element code (FEFLOW). This comparison was performed for several test cases with increasing complexity. The differences between the flux-tube approach and the other methods always remain small, even regarding mass conservation. The major advantage of the flux-tube approach is the ability to reach a solution quickly, as the method is hundreds to thousands of times faster than classical finite difference or finite element models.     

Effective velocities in fractured media: a numerical study using the rotated staggered finite‐difference grid

The modelling of elastic waves in fractured media with an explicit finite‐difference scheme causes instability problems on a staggered grid when the medium possesses high‐contrast discontinuities (strong heterogeneities). For the present study we apply the rotated staggered grid. Using this modified grid it is possible to simulate the propagation of elastic waves in a 2D or 3D medium containing cracks, pores or free surfaces without hard‐coded boundary conditions. Therefore it allows an efficient and precise numerical study of effective velocities in fractured structures. We model the propagation of plane waves through a set of different, randomly cracked media. In these numerical experiments we vary the wavelength of the plane waves, the crack porosity and the crack density. The synthetic results are compared with several static theories that predict the effective P‐ and S‐wave velocities in fractured materials in the long wavelength limit. For randomly distributed and randomly orientated, rectilinear, non‐intersecting, thin, dry cracks, the numerical simulations of velocities of P‐, SV‐ and SH‐waves are in excellent agreement with the results of the modified (or differential) self‐consistent theory. On the other hand for intersecting cracks, the critical crack‐density (porosity) concept must be taken into account. To describe the wave velocities in media with intersecting cracks, we propose introducing the critical crack‐density concept into the modified self‐consistent theory. Numerical simulations show that this new formulation predicts effective elastic properties accurately for such a case.     

A~200 ka pollen record from Okinawa Trough:Paleoenvironment reconstruction of glacial-interglacial cycles

Pollen analysis was carried out on the Core MD982194 of 29.78 m retrieved from the Middle Okinawa Trough which was dated as old as~200 ka BP.The results revealed that pollen assemblages mainly presented an alternation of coniferous and herb pollen.The coniferous saccate pollen,principally Pinus and Tsuga,predominated in most parts of the core,especially highlighted in the interstadial stages including MIS 1,3,5 and MIS7,whereas the herb pollen significantly increased in the glacial periods.Thus the pollen flora and their percentages showed the sensitive changes under the influence of ice volume during the glacial and interglacial periods.Our record from this core has first documented that the percentage of Cyperaceae was extremely high in the glacial stages with a notable increase in Artemisia,Gramineae,Asteraceae,Chenopodiaceae,and freshwater algae,which can be used as a proxy for sea-level change at the study site because of their close negative correlation of the orbital-scale changes in sea level.The distance between the continental coastline and the Okinawa Trough has deeply shortened due to the sea-level drop in the LGM.As a result,the sediment materials from Yangtze River were extensively deposited on the flat,exposed continental shelf owing to the rapid decline of river flow speed,leading to that pollen grains from Okinawa Trough are derived mainly from the flat coastal vegetation of exposed continental shelf at glacial stages.Changes of pollen assemblage were consistent with the variation of temperature and humidity,which showed that the percentage of arboreal pollen was highly augmented at MIS 7,5 and MIS 1,corresponding to the strengthening of the East Asian summer monsoon and increasing of rainfall.Moreover,the peak of Pinus percentage in MIS 5.3,5.1 and MIS 3.3 may be closely linked by orbital and sub-orbital cycles of solar radiation and monsoonal variability.The present study of core MD981294 implied lower temperature and precipitation during the lowest sea-level stage(LGM),and more visibly testified that the vegetation of the flat plain on the exposed continental shelf was dominated by intrazonal communities such as halophyte grasslands and freshwater wetlands instead of zonal steppe or semi-arid desert.All above evidence demonstrated that the fundamental changes of pollen assemblage and their origins in Okinawa Trough since~200 ka BP were affected by combine factors including the coastline position and climate fluctuation.Moreover,the substantial shortening of distance between shoreline and the Okinawa Trough driven by orbital insolation cycles was clearly indicated by the pollen spectra,whereas the source-area climate signal of the pollen record was largely weakened.     

Cementation of kerogen-rich marls by alkaline fluids released during weathering of thermally metamorphosed marly sediments. Part I: Isotopic (C,O) study of the Khushaym Matruk natural analogue (central Jordan)

The Khushaym Matruk site in central Jordan may represent a natural analogue depicting the interaction of alkaline solutions with a clayey sedimentary formation or with clay-rich confining barriers at the interface with concrete structures in waste disposal sites. In this locality, past spontaneous combustion of organic matter in a clayey biomicritic formation produced a ca. 60 m-thick layer of cement-marble containing some of the high-temperature phases usually found in industrial cements (e.g., spurrite, brucite, and Ca-aluminate). A vertical cross-section of the underlying sediments was used in order to study the interaction between cement-marbles and neighbouring clayey limestones under weathering conditions. A thermodynamic approach of the alteration parageneses (calcite–jennite–afwillite–brucite and CSH phases) in the cement-marbles constrains the interacting solutions to have had pH-values between 10.5 and 12. Over 3 m, the sediments located beneath the metamorphic unit were compacted and underwent carbonation. They display large C and O isotopic variations with respect to “pristine” sediments from the bottom of the section. Low δ¹³C-values down to −31.4‰/PDB show the contribution of CO₂ derived from the oxidization of organic matter and from the atmosphere to the intense carbonation process affecting that particular sedimentary level. The size of the C isotopic anomalies, their geometrical extent and their coincidence with the variations of other markers like the Zn content, the structure of organic matter, the mineralogical composition, all argue that the carbonation process was induced by the percolation of high pH solutions which derived from the alteration of cement-marbles. The temperature of the carbonation process remains conjectural and some post-formation O isotopic reequilibration likely affected the newly-formed carbonate. Carbonation induced a considerable porosity reduction, both in fractures and matrixes. The Khushaym Matruk site may have some bearing to the early life of a repository site, when water saturation of the geological formations hosting the concrete structures is incomplete, enabling simultaneous diffusion of alkaline waters and gaseous CO₂ in the near field.     

Stratospheric ozone measurements by solar ultraviolet absorption

Vertical distributions of ozone in the stratosphere have been determined by means of u.v. differential absorption technique, using the Sun as a source. The filter radiometer and the data processing are described. Results obtained between altitudes of 22 and 32 km during the Intercomparison Ozone Campaign held in France in June 1981 are presented and discussed.     

Predictive Geological Mapping Using Closed-Form Non-stationary Covariance Functions with Locally Varying Anisotropy: Case Study at El Teniente Mine (Chile)

This paper is concerned with the problem of predicting the surface elevation of the Braden breccia pipe at the El Teniente mine in Chile. This mine is one of the world’s largest and most complex porphyry-copper ore systems. As the pipe surface constitutes the limit of the deposit and the mining operation, predicting it accurately is important. The problem is tackled by applying a geostatistical approach based on closed-form non-stationary covariance functions with locally varying anisotropy. This approach relies on the mild assumption of local stationarity and involves a kernel-based experimental local variogram a weighted local least-squares method for the inference of local covariance parameters and a kernel smoothing technique for knitting the local covariance parameters together for kriging purpose. According to the results, this non-stationary geostatistical method outperforms the traditional stationary geostatistical method in terms of prediction and prediction uncertainty accuracies.     

Post-collisional transition from calc-alkaline to alkaline volcanism during the Neogene in Oranie (Algeria): magmatic expression of a slab breakoff

During the Neogene, a magmatic change from calc-alkaline to alkaline types occurred in all the regions surrounding the western Mediterranean. This change has been studied in Oranie (western Algeria). In this area, potassic to shoshonitic calc-alkaline andesites (with La/Nb ratios in the range 4–6) were mainly erupted between 12 and 9 Ma. They were followed (between 10 and 7 Ma) by basalts displaying geochemical features which are transitional between calc-alkaline and alkaline lavas (La/Nb=1–1.7). After a ca. 3-Ma quiescence period, volcanic activity resumed, with the eruption of OIB-type alkaline basalts (La/Nb=0.5–0.6), from 4 to 0.8 Ma. A combined geochemical approach, using incompatible elements and Sr, Nd and O isotopes, allows us to conclude that the transitional basalts derived from the melting of a heterogeneous mantle source, at the boundary between lithosphere and asthenosphere. We propose that melting of a previously subduction-modified lithospheric mantle occurred between 12 and 10 Ma, in response to the upwelling of hot asthenosphere flowing up into an opening gap above a detached sinking slab. As a result, calc-alkaline magmas were formed. From 10 to 7 Ma, the transitional basalts were generated through melting of the boundary mantle zone between the lithosphere and the upwelling asthenosphere. During that stage, the contribution of the lithospheric source was still predominant. Then, as sinking of the oceanic slab progressed, the increasing uprise of the asthenosphere led to the formation and emplacement (from 4 to 0.8 Ma) of typical within-plate alkaline basalts derived from a plume-modified asthenospheric mantle.     

Interactions climat-eustatisme-tectonique. Les enseignements et perspectives du Crétacé supérieur (Cénomanien-Coniacien)

Résumé

Une analyse de séries sédimentaires d’âge Crétacé supérieur dans le Sud-Est de la France a été entreprise et intégrée dans un cadre de stratigraphie séquentielle. Les limites Cénomanien–Turonien et la limite Turonien–Coniacien sont marquées par des variations rapides et de fortes amplitudes du niveau marin relatif. La comparaison effectuée avec d’autres bassins mondiaux nous amène à envisager le problème dans un cadre global. On constate : (1) le synchronisme d’événements à haute fréquence dans des contextes géodynamiques différents. (2) les effets de la superposition de cycles eustatiques hiérarchiquement différents, c’est-à-dire la superposition d’oscillations à haute fréquence sur une composante de 3^e ordre. Ces observations sont confrontées à deux hypothèses : le glacio-eustatisme et la tectonique à haute fréquence. © Elsevier, Paris.