Identification of Groundwater Parameters Using an Adaptative Multiscale Method

The identification of groundwater parameters in heterogeneous systems is a major challenge in groundwater modeling. Flexible parameterization methods are needed to assess the complexity of the spatial distributions of these parameters in real aquifers. In this article, we introduce an adaptative parameterization to identify the distribution of hydraulic conductivity within the large‐scale (4400 km²) Upper Rhine aquifer. The method is based on adaptative multiscale triangulation (AMT) coupled with an inverse problem procedure that identifies the parameters' distributions by reducing the error between measured and simulated heads. The AMT method has the advantage of combining both zonation and interpolation approaches. The AMT method uses area‐based interpolation rather than an interpolation based on stochastic features. The method is applied to a standard 2D groundwater model that takes into account the interactions between the aquifer and surface water bodies, groundwater recharge, and pumping wells. The simulation period covers 204 months, from January 1986 to December 2002. Recordings at 109 piezometers are used for model calibration. The simulated heads are globally quite accurate and reproduce the main dynamics of the system. The local hydraulic conductivities resulting from the AMT method agree qualitatively with existing local experimental observations across the Rhine aquifer.     

Generation mechanisms for mesoscale eddies in the Gulf of Lions: radar observation and modeling

Coastal mesoscale eddies were evidenced during a high-frequency radar campaign in the Gulf of Lions (GoL), northwestern Mediterranean Sea, from June 2005 to January 2007. These anticyclonic eddies are characterized by repeated and intermittent occurrences as well as variable lifetime. This paper aims at studying the link between these new surface observations with similar structures suggested at depth by traditional acoustic Doppler current profiler measurements and investigates the eddy generation and driving mechanisms by means of an academic numerical study. The influence of the wind forcing on the GoL circulation and the eddy generation is analyzed, using a number of idealized configurations in order to investigate the interaction with river discharge, buoyancy, and bathymetric effects. The wind forcing is shown to be crucial for two different generation mechanisms: A strong northerly offshore wind (Mistral) generates a vortex column due to the bathymetric constraint of a geostrophic barotropic current, which can surface after the wind relaxes; a southerly onshore wind generates a freshwater bulge from the Rhône river discharge, which detaches from the coast and forms a well-defined surface anticyclonic eddy based on buoyancy gradients. These structures are expected to have important consequences in terms of dispersion or retention of biogeochemical material at local scales.     

Etude des Déformations et Inclinaisons du Sous-sol: Développements Instrumentaux et Applications dans la Région Sismique d’Arette,France

We present a modified version of a tilt and strain meter conceived to have great sensibility and good stability. These instruments were deployed in a cave at the seismic Pyrenean region of Arette (France). We set in evidence several types of phenomena, of which deformations brought about by water recharge filling up fissures on the massif are of importance. Owing to the redundance in the data, the study of the co-seismic deformations allows one to separate instrumental displacements and local readjustments from regional motions.     

A preliminary investigation of strong-motion data from the French Antilles

Strong-motion networks have been operating in the Caribbean region since the 1970s, however, until the mid-1990s only a few analogue stations were operational and the quantity of data recorded was very low. Since the mid-1990s, digital accelerometric networks have been established on islands within the region. At present there are thought to be about 160 stations operating in this region with a handful on Cuba, 65 on the French Antilles (mainly Guadeloupe and Martinique), eight on Jamaica, 78 on Puerto Rico (plus others on adjacent islands) and four on Trinidad.After briefly summarising the available data from the Caribbean islands, this article is mainly concerned with analysing the data that has been recorded by the networks operating on the French Antilles in terms of their distribution with respect to magnitude, source-to-site distance, focal depth and event type; site effects at certain stations; and also with respect to their predictability by ground motion estimation equations developed using data from different regions of the world. More than 300 good quality triaxial acceleration time-histories have been recorded on Guadeloupe and Martinique at a large number of stations from earthquakes with magnitudes larger than 4.8, however, most of the records are from considerable source-to-site distances. From the data available it is found that many of the commonly-used ground motion estimation equations for shallow crustal earthquakes poorly estimate the observed ground motions on the two islands; ground motions on Guadeloupe and Martinique have smaller amplitudes and are more variable than expected. This difference could be due to regional dependence of ground motions because of, for example, differing tectonics or crustal structures or because the ground motions so far recorded are, in general, from smaller earthquakes and greater distances than the range of applicability of the investigated equations.     

Eruptive history of the Colima volcanic complex (Mexico)

The evolution of the Colima volcanic complex can be divided into successive periods characterized by different dynamic and magmatic processes: emission of andesitic to dacitic lava flows, acid-ash and pumice-flow deposits, fallback nuées ardentes leading to pyroclastic flows with heterogeneous magma, plinian air-fall deposits, scoriae cones of alkaline and calc-alkaline nature. Four caldera-forming events, resulting either from major ignimbrite outbursts or Mount St. Helens-type eruptions, separate the main stages of development of the complex from the building of an ancient shield volcano (25 × 30 km wide) up to two summit cones, Nevado and Fuego.The oldest caldera, C1 (7–8 km wide), related to the pouring out of dacitic ash flows, marks the transition between two periods of activity in the primitive edifice called Nevado I: the first one, which is at least 0.6 m.y. old, was mainly andesitic and effusive, whereas the second one was characterized by extrusion of domes and related pyroclastic products. A small summit caldera, C2 (3–3.5 km wide), ended the evolution of Nevado I.Two modern volcanoes then began to grow. The building of the Nevado II started about 200,000 y. ago. It settled into the C2 caldera and partially overflowed it. The other volcano, here called Paleofuego, was progressively built on the southern side of the former Nevado I. Some of its flows are 50,000 y. old, but the age of its first outbursts is not known. However, it is younger than Nevado II. These two modern volcanoes had similar evolutions. Each of them was affected by a huge Mount St. Helens-type (or Bezymianny-type) event, 10,000 y. ago for the Paleofuego, and hardly older for the Nevado II. The landslides were responsible for two horseshoe-shaped avalanche calderas, C3 (Nevado) and C4 (Paleofuego), each 4–5 km wide, opening towards the east and the south. In both cases, the activity following these events was highly explosive and produced thick air-fall deposits around the summit craters.The Nevado III, formed by thick andesitic flows, is located close to the southwestern rim of the C3 caldera. It was a small and short-lived cone. Volcan de Fuego, located at the center of the C4 caldera, is nearly 1500 m high. Its activity is characterized by an alternation of long stages of growth by flows and short destructive episodes related to violent outbursts producing pyroclastic flows with heterogeneous magma and plinian air falls.The evolution of the primitive volcano followed a similar pattern leading to formation of C1 and then C2. The analogy between the evolutions of the two modern volcanoes (Nevado II–III; Paleofuego-Fuego) is described. Their vicinity and their contemporaneous growth pose the problem of the existence of a single reservoir, or two independent magmatic chambers, after the evolution of a common structure represented by the primitive volcano.     

Towards Improved Remedial Outcomes in Categorical Aquifers with an Iterative Ensemble Smoother

Categorical parameter distributions consisting of geologic facies with distinct properties, for example, high-permeability channels embedded in a low-permeability matrix, are common at contaminated sites. At these sites, low-permeability facies store solute mass, acting as secondary sources to higher-permeability facies, sustaining concentrations for decades while increasing risk and cleanup costs. Parameter estimation is difficult in such systems because the discontinuities in the parameter space hinder the inverse problem. This paper presents a novel approach based on Traveling Pilot Points (TRIPS) and an iterative ensemble smoother (IES) to solve the categorical inverse problem. Groundwater flow and solute transport in a hypothetical aquifer with a categorical parameter distribution are simulated using MODFLOW 6. Heads and concentrations are recorded at multiple monitoring locations. IES is used to generate posterior ensembles assuming a TRIPS prior and an approximate multi-Gaussian prior. The ensembles are used to predict solute concentrations and mass into the future. The evaluation also includes an assessment of how the number of measurements and the choice of the geological prior determine the characteristics of the posterior ensemble and the resulting predictions. The results indicate that IES was able to efficiently sample the posterior distribution and showed that even with an approximate geological prior, a high degree of parameterization and history matching could lead to parameter ensembles that can be useful for making certain types of predictions (heads, concentrations). However, the approximate geological prior was insufficient for predicting mass. The analysis demonstrates how decision-makers can quantify uncertainty and make informed decisions with an ensemble-based approach.     

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.     

himu-em: The French Polynesian connection

himu, em i andem ii are three of the main geochemical mantle components that give rise to oceanic island basalts [1]. They represent the end members that produce the extreme isotopic compositions measured on intraplate volcanics. In French Polynesia, all three mantle components are represented in volcanic rocks. The characteristichimu signature is found in Tubuai, Mangaia and Rurutu,em i is present in the source of Rarotonga and Pitcairn volcanics andem ii dominates the composition of most Society Islands. Intermediate values between the three end members are found on most islands.We suggest that the three components are not independent but are physically related in the mantle. Thehimu component is thought to be recycled oceanic crust that lost part of its Pb through hydrothermal processes prior to and during subduction.em i andem ii are believed to acquire their isotopic and trace element characteristics through entrainment of sediments that were subducted together with the oceanic crust.The trace element pattern and the isotopic composition ofhimu lavas can be quantitatively modelled using a mixture of 25% old recycledmorb crust and 75% mantle peridotite. The extreme Pb composition is modelled assuming that Pb was lost from oceanic crust when hydrothermal alteration at the ridge leached Pb from the basalt to redeposit it as sulphides on top of and throughout the crust, followed by preferential dissolution of sulphides during dehydration in the subduction zone. These processes led to a drastic increase of theU/Pb ratio of the subducted material which evolved over 2 Ga to very radiogenic Pb isotopic compositions. Pb isotopic compositions similar to those ofem i andem ii are modelled assuming that sediments with average crustal Pb isotopic compositions were subducted and recycled into the mantle together with the underlyingmorb oceanic crust. Pelagic sediments (μ 5 andκ 6) account for the Pb isotopic composition ofem i whereas terrigenous sediments (μ 10 andκ 4.5) evolve towards theem ii end member. A few percent of sediment in the recycled crust-sediment mixture will destroy the characteristic Pb isotopic signature of thehimu component. This, together with the low probability of isolating oceanic crust in the mantle for 2 Ga, explains why the extremehimu composition, as seen on Tubuai and St Helena, is sampled so rarely by oceanic volcanism.     

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.     

Early contamination of European flounder (Platichthys flesus) by PCDD/Fs and dioxin-like PCBs in European waters

Contamination levels and profiles of 7 polychlorinated-p-dioxins, 10 polychlorinated furans (PCDD/Fs) and 12 dioxin-like polychlorinated biphenyls (dl-PCBs) were investigated in juvenile European flounder (Platichthys flesus) captured in different nursery areas in the northeastern Atlantic coast across its geographical distribution range. The toxic equivalent concentrations (WHO-TEQ_fish) were also determined in order to evaluate which P. flesus population was more exposed to dioxin-like toxicity. Juveniles caught in the Sørfjord (Norway) showed the lowest WHO-TEQ_fish concentration (0.052 pg WHO-TEQ_fish g⁻¹ wet weight) whereas the highest value was observed in fish from the Wadden Sea (The Netherlands; 0.291 pg WHO-TEQ_fish g⁻¹ ww), mainly due to the greater contribution of 2,3,7,8-tetrachlorodibenzo-p-dioxin, the most toxic congener. Nonetheless, when comparing the results with existent tissue residue-based toxicity benchmarks, no adverse effects resulting from PCDD/Fs and dl-PCBs are expected to occur in flounder from the studied systems.