Generation of 3D Spatially Variable Anisotropy for Groundwater Flow Simulations

Sedimentary units generally present anisotropy in their hydraulic properties, with higher hydraulic conductivity along bedding planes, rather than perpendicular to them. This common property leads to a modeling challenge if the sedimentary structure is folded. In this paper, we show that the gradient of the geological potential used by implicit geological modeling techniques can be used to compute full hydraulic conductivity tensors varying in space according to the geological orientation. For that purpose, the gradient of the potential, a vector normal to the bedding, is used to construct a rotation matrix that allows the estimation of the 3D hydraulic conductivity tensor in a single matrix operation. A synthetic 2D cross section example is used to illustrate the method and show that flow simulations performed in such a folded environment are highly influenced by this rotating anisotropy. When using the proposed method, the streamlines follow very closely the folded formation. This is not the case with an isotropic model.     

New Methods to Estimate 2D Water Level Distributions of Dynamic Rivers

River restoration measures are becoming increasingly popular and are leading to dynamic river bed morphologies that in turn result in complex water level distributions in a river. Disconnected river branches, nonlinear longitudinal water level profiles and morphologically induced lateral water level gradients can evolve rapidly. The modeling of such river‐groundwater systems is of high practical relevance in order to assess the impact of restoration measures on the exchange flux between a river and groundwater or on the residence times between a river and a pumping well. However, the model input includes a proper definition of the river boundary condition, which requires a detailed spatial and temporal river water level distribution. In this study, we present two new methods to estimate river water level distributions that are based directly on measured data. Comparing generated time series of water levels with those obtained by a hydraulic model as a reference, the new methods proved to offer an accurate and faster alternative with a simpler implementation.     

Sedimentary archives of the French Atlantic coast (inner Bay of Vilaine,south Brittany): Depositional history and late Holocene climatic and environmental signals

The late Holocene is of particular interest to our understanding of the evolution of coastal sedimentary systems because this period encompasses warmer and cooler periods, and rising sea level in northern Europe. Based on an approach combining AMS ¹⁴C, sedimentological and rock magnetic analyses on sediment cores complemented with seismic data collected in the macrotidal Bay of Vilaine (south Brittany), we document the depositional history of the inner bay coeval to the mid- to late-Holocene transgression in south Brittany. Correlation between sedimentary archives revealed the main sedimentary infilling phases during the last 6000 years. Four units (U1–U4) are recognized in the coastal sediment wedge of the system, corresponding to the stepwise marine invasion of the bay. We show that (1) marine inundation, due to the steep morphology of the bedrock, is diachronous between distal and proximal records. A time lag of ～1000 years is inferred over a distance of less than 5 km; (2) in the outer areas, the sedimentation has been condensed since 3000 years; (3) proximal estuarine archives offer the best record of sedimentary processes covering the last 2000 years, including the Medieval Warm Period (MWP).Correlations in proximal records in the Bay of Vilaine assess the connection between coastal sedimentary dynamics, climatic conditions and anthropogenic activities during the MWP. We match the preservation of clay deposits to increased river-borne suspended matter transported to the estuary probably as a result of accelerated land-use development (higher soil erosion) in the catchment area between ca. 880 and 1050 AD. Because the preservation of estuarine sedimentary successions is favoured when coastal wave sediment reworking is minimal, it is proposed that the prevailing climatic regime in south Brittany during the MWP likely resembled to that of the preferred negative phase of the North Atlantic Oscillation (NAO). Our data are fairly consistent with other late Holocene records from northern Europe including the Atlantic seaboard. However, they outline the difficulty in interpreting climatic and anthropogenic signatures in coastal sedimentary records where high-resolution chronologies required to unravel their respective influences are still missing.     

Calderas,landslides and paleo-canyons on Piton de la Fournaise volcano (La Réunion Island,Indian Ocean)

Based upon a re-interpretation of previous data and a new field campaign, a structural evolution is proposed for the early history of Piton de la Fournaise volcano from 500,000 to 50,000 years. Conceptually, it is shown that the formation of a caldera in which lava flows are contained inside the caldera depression, gives time for erosion to excavate deep canyons on the external slopes of the volcano, for example, the Rivière des Remparts, the Rivière Langevin and the Rivière de l'Est canyons on Piton de la Fournaise volcano. These canyons are infilled when lavas, filling the caldera and overflowing its rim, are able again to flow on the external slopes of the volcano. In the past, this excavating/infilling process has occurred twice following the formation of the Rivière des Remparts and Morne Langevin calderas. The formation of the third caldera, the Plaine des Sables caldera, was followed by the excavation of the current canyons. In addition to this process, two large landslides have been documented in the field. The first, which happened about 300,000 years ago, is apparently the first episode of the break up of Piton de la Fournaise volcano, predating the formation of the four large calderas. The second landslide, which occurred 150,000 years ago and is considered to be less extensive, has carried away the entire southern flank of the Rivière des Remparts caldera.     

The Grand Brûlé exploration drilling: New data on the deep framework of the Piton de la Fournaise volcano. Part 2: Secondary minerals

The Grand Brûlé borehole intersects a thick pile of basaltic lavas, down to 1010 m, and a basic-ultrabasic intrusive complex, from 1010 to 3003.50 m.The lavas are, in general, unaltered except in two fractured zones, where hydrothermal fluids circulated at temperatures not exceeding 350 ° C. The main secondary minerals there are pyroxene, feldspar, epidote, actinolite and chlorite.The entire thickness of the intrusive body intersected contains secondary minerals representing three stages of cooling:

1. 1. A late magmatic episode (600–900°C) characterized by biotite, kaersutite, edenite and pargasite.

2. 2. A hydrothermal episode (T<- 350°C) characterized by epidote, albite, biotite, actinolite and chlorite.

3. 3. A phase of serpentinization (T≤ 350°C).

It is very likely that the two later events occurred simultaneously, with physical and chemical interference.     

A 3-D model for the Antarctic ice sheet: a sensitivity study on the glacial-interglacial contrast

On the longer climatic time scales, changes in the elevation and extent of the Antarctic ice sheet have an important role in modulating global atmospheric and oceanographic processes, and contribute significantly to world-wide sea levels. In this paper, a 3-D time-dependent thermomechanical model for the entire ice sheet is presented, that is subsequently used to examine the effects of glacial-interglacial shifts in environmental boundary conditions on its geometry. The model takes into account a coupled ice shelf, grounding-line dynamics, basal sliding and isostatic bed adjustment and considers the fully coupled velocity and temperature fields. Ice flow is calculated on a fine mesh (40 km horizontal grid size and 10 layers in the vertical) for grounded and floating ice and a stress transition zone in between at the grounding line, where all stress components contribute in the effective stress in the flow law. There is free interaction between ice sheet and ice shelf, so that the entire geometry is internally generated. A simulation of the present ice sheet reveals that the model is able to yield realistic results. A series of sensitivity experiments are then performed, in which lower temperatures, reduced accumulation rates and lower global sea level stands are imposed, either singly or in combination. By comparing results of pairs of experiments, the effects of each of these environmental changes can be determined. In agreement with glacial-geological evidence, we found that the most pronounced changes show up in the West Antarctic ice sheet configuration. They appear to be essentially controlled by variations in eustatic sea level, whereas typical glacial-interglacial changes in temperature and ice deposition rates tend to balance one another. These findings support the hypothesis that the Antarctic ice sheet basically follows glacial episodes in the northern hemisphere by means of sea-level teleconnections. Grounding occurs more readily in the Weddell sea than in the Ross sea and long time scales appear to be involved: it may take up to 30–40000 years for these continental shelf areas to become completely grounded after an initial stepwise perturbation in boundary conditions. According to these reconstructions, a steady state Antarctic ice sheet may contribute some 16 m to global sea level lowering at maximum glaciation.     

Isotopic study of the Manaslu granite (Himalaya,Nepal): inferences on the age and source of Himalayan leucogranites

A detailed isotopic study of the Manaslu leucogranite was carried out. A U-Pb age of 25 Ma and a whole rock Rb-Sr age isochron of 18 Ma were obtained, suggesting that the magmatic activity lasted at least 7 Ma. Initial Sr isotopic ratios are very high (0.740 to 0.760) and initial Nd isotopic ratios are low ( _Nd ⁱⁿ : –13 to –16), and they show the existence of large isotopic variations even at the metre scale. These are not the result of perturbations by fluids but rather they reflect the initial isotopic heterogeneity of the source material which has not been obliterated by magmatic processes (e.g. fusion, mixing by convection). These results also support the crustal origin of this leucogranite. The Tibetan slab paragneisses, whose Sr and Nd isotopic ratios are very similar to those of the granite at an age of 20 Ma, are the most probable parental material. Nd model ages for both the leucogranite and the gneisses are in the range 1.5–2 Ga. A model of formation of the Manaslu granite by coalescence of different batches of magma is in agreement with the present data.