Advantages of the analytic element method for the solution of groundwater management problems

In the simulation‐optimization approach, a coupled optimization and groundwater flow/transport model is used to solve groundwater management problems. The efficiency of the numerical method, which is used to simulate the groundwater flow, is one the major reason to obtain the best solution for a management problem. This study was carried out to examine the advantages of the analytic element method (AEM) in the simulation‐optimization approach, for the solution of groundwater management problems. For this study, the AEM and finite difference method (FDM) based flow models were developed and coupled with the particle swarm optimization (PSO)‐based optimization model. Furthermore, the AEM‐PSO and FDM‐PSO models developed were applied in hypothetical as well as real field conditions to address groundwater management problems and the results were compared. For the real field situation, the models developed were applied to the Dore River basin in France to minimize the installation and operational cost of new pumping wells taking the location and discharge of the pumping wells as decision variables. The constraints of the problem were identified with the help of stakeholders and water authority officials. The AEM flow model was developed to facilitate the management model particularly when at each iteration, the optimization model calls for a simulation model to calculate the values of groundwater heads. The results show that, at some points, the AEM‐PSO model is efficient in identifying the optimal location of wells and consequently results in optimal costs, sometimes difficult when using the FDM. Copyright © 2011 John Wiley & Sons, Ltd.     

Stratigraphic simulations of the shelf of the Gulf of Lions: testing subsidence rates and sea‐level curves during the Pliocene and Quaternary

Determining the relative importance of factors interacting to control stratigraphic organization is a key issue in sedimentology. The Pliocene‐Quaternary chronostratigraphy on the Gulf of Lions platform is still poorly constrained, giving rise to different interpretations of the evolution of its subsidence through time. This paper examines the Pliocene‐Quaternary sedimentary filling of the Gulf of Lion's shelf with Dionisos, a numerical stratigraphic model. Our results show that a constant subsidence rate accurately reproduces the observed geometries, whereas a varying subsidence rate reproduces them only if the acceleration of subsidence is limited. At this time‐scale, a third‐order eustatic curve is also reappraised: a higher resolution curve (built using δ¹⁸O measurements) gives a more realistic restitution of our stratigraphic markers. Finally, the constant subsidence rate and sediment fluxes implied in these modellings are discussed relative to climate and local factors of sedimentation.     

Estimating water fluxes in the critical zone using water stable isotope approaches in the Groundnut and Ferlo basins of Senegal

Sustainable water management in semi-arid agriculture practices requires quantitative knowledge of water fluxes within the soil-vegetation-atmosphere system. Therefore, we used stable-isotope approaches to evaluate evaporation (E_a), transpiration (T_a), and groundwater recharge (R) at sites in Senegal's Groundnut basin and Ferlo Valley pasture region during the pre-monsoon, monsoon, and post-monsoon seasons of 2021. The approaches were based upon (i) the isothermal evaporation model (for quantifying E_a); (ii) water and isotope mass balances (to partition E_a and T_a for groundnut and pasture); and (iii) the piston displacement method (for estimating R). E_a losses derived from the isothermal evaporation model corresponded primarily to Stage II evaporation, and ranged from 0.02 to 0.09 mm d⁻¹ in the Groundnut basin, versus 0.02–0.11 mm d⁻¹ in Ferlo. At the groundnut site, E_a rates ranged from 0.01 to 0.69 mm d⁻¹; T_a was in the range 0.55–2.29 mm d⁻¹; and the T_a/ET_a ratio was 74%–90%. At the pasture site, the ranges were 0.02–0.39 mm d⁻¹ for E_a; 0.9–1.69 mm d⁻¹ for T_a; and 62–90% for T_a/ET_a. The ET_a value derived for the groundnut site via the isotope approach was similar to those from eddy covariance measurements, and also to the results from the previous validated HYDRUS-1D model. However, the HYDRUS-1D model gave a lower T_a/ET_a ratio (23.2%). The computed groundwater recharge for the groundnut site amounted to less than 2% of the local annual precipitation. Recommendations are made regarding protocols for preventing changes to isotopic compositions of water in samples that are collected in remote arid regions, but must be analysed days later. The article ends with suggestions for studies to follow up on evidence that local aquifers are being recharged via preferential pathways.     

Tertiary tectono-metamorphic evolution of the European margin during Alpine collison: example of the Leventina Nappe (Central Alps,Switzerland)

The Leventina Nappe represents one of the lowermost exposed units in the Alpine nappe stack and corresponds to a slice of the European margin that was entrained into the Alpine continental accretionary prism during the Tertiary tectonic event. This study yields details regarding the tectonic and metamorphic history of the Leventina Nappe, through detailed analysis of structures and shear zone patterns, and the examination of the Si-content of white mica along a north-south profile. The Leventina Nappe underwent three phases of ductile deformation. Foliation S1 is mostly sub-parallel to the regionally dominant structural fabric (the S2 foliation). S2 foliation is penetratively developed in the structurally higher portions of the Leventina Nappe toward the Simano Nappe, while it is only weakly developed in the core of the Leventina Nappe. A 50 to 200 m wide mylonite zone, with a D2 top-to-NW sense of shear marks the boundary to the Simano Nappe. Throughout the Leventina Nappe only small-scale D2 shear bands (mm to cm wide) are observed, showing a top-to-NW sense of shear. Deformation phase D3 locally generated a vertical axial plane foliation (S3) associated with the large-scale D3 Leventina antiform.Microtextural evidence and phengite geobarometry were used to constrain the temperature and pressure conditions of equilibration of the Leventina Gneisses. Highest Si (pfu) values are preserved in the core of phengitic micas and reflect pressure and temperature conditions of around 8 kbar at 550 °C and 10 kbar at 650 °C in the northern and southern parts of the Leventina Nappe, respectively. Lower Si (pfu) values from the rims of white micas correspond to a metamorphic pressure of ca. 5 kbar during the exhumation of the unit. These metamorphic conditions are related to the underthrusting of the thinned European margin into the continental accretionary prism during late Eocene time. These new data allow us to propose a kinematic model for the Leventina Nappe during the Tertiary Alpine tectonics.     

Petrographical and geochemical characterization of Comblanchien limestone (Bourgogne,France): A fingerprint of the building stone provenance

The Comblanchien limestone is a well-known building stone that deserves a quantitative characterization in order to be able to certify its origin. We compare petrographical and geochemical characteristics of limestones commercialized under the name of Comblanchien to those of one of their competitors (VATC, a Portuguese coated grain-bearing limestone). From macroscopic and microscopic observations, the distinction is somehow difficult but the coated grain-bearing Comblanchien samples have δ¹³C values (from +1.75 to +2.25‰) different from those of the VATC sample (negative δ¹³C value). As an example of application of these results, three samples from a building site near Paris are also studied in order to verify their provenance. We conclude that they do not come from quarries from the Comblanchien area, but likely from Portugal. The study shows that a detailed analysis of building limestones may provide a robust tool for discriminating between currently mined stones of different provenances, and thus open the possibility to attach to a given limestone a fingerprint, an identity card, and if required a commercial label.     

Progressive degradation of rock properties and time-dependent behavior of deep tunnels

Time-dependent response of deep tunnels is studied considering the progressive degradation of the mechanical properties of the rock mass. The constitutive model is based on a rock-aging law for the uniaxial strength of the rock and for the Young’s modulus. A semi-analytical solution is developed for the stresses and displacements around a deep circular tunnel taking into account the face advance. The evolution of the plastic and damage zones over time is determined. Numerical examples are presented for the case of Saint-Martin-La-Porte access adit in France of the Lyon–Turin Base Tunnel. The computed results which are compared with the field data in terms of the convergence of tunnel wall and of the displacements inside the rock mass monitored by multi-point extensometers show the efficiency of the approach to simulate the time-dependent deformation of a tunnel excavated in squeezing ground. Simple relationships are proposed to evaluate the parameters of the constitutive model directly from those of the empirical convergence law presented in previous work.     

Trace element geochemistry by laser ablation ICP-MS of micas associated with Ta mineralization in the Tanco pegmatite,Manitoba, Canada

In the Tanco pegmatite, one of the world’s major Ta deposits, tantalum mineralization shows a complexity that reflects the complex petrogenesis of its host pegmatite. Micas are common in most of the pegmatite units and are intimately associated with the successive stages of Ta mineralization, from the wall zone to the central zones where micaceous replacement is pervasive. Different generations of micas, both primary and secondary, associated with Ta oxides, were selected for electron microprobe and laser ablation ICP-MS investigation. Their chemical trends are used to constrain the magmatic versus hydrothermal processes that played a role in their crystallization and their associated Ta mineralization. Micas range from dioctahedral muscovite to trioctahedral lepidolite through Al↔Li substitution. Unexpectedly, the most evolved compositions (low K/Rb ratios and high Li contents) occur in the wall zone; they are interpreted to reflect nonequilibrium crystallization from an undercooled melt, with or without boundary layer effects. In the central zones, the fine-grained mica–quartz assemblage hosts some coarser-grained Li-muscovite, which has the highest Ta contents (up to 400 ppm). These Li–F–a-rich micas are interpreted to originate from a magmatic metasomatic event, which was also at the origin of the MQM-style Ta mineralization at Tanco. However, the Li–Ta-poor, muscovite end-member compositions of fine-grained alteration micas suggest crystallization from an aqueous fluid, during a metasomatic (hydrothermal) event involving late pegmatitic fluids. The low Ta concentrations (around 50 ppm) of this fine-grained muscovite suggest that this fluid transported at least small amounts of Ta.     

单个裂隙中油运移实验及特征分析

利用狭窄平行板裂隙模型来模拟和分析单个裂隙内油运移的过程。实验中采用表面起伏的玻璃板制作平行板裂隙模型,通过改变水中蔗糖的浓度控制两相流体间的密度差,来分析不同密度差和注入速率条件下油的运移过程和特征。观察发现,油在饱含水的狭窄裂隙空间内的运移特征与在孔隙介质中的运移有相似之处:运移的路径只占通道的一部分,都具有不规则的分形特征,也可形成活塞式、指进式和路径式等模式。而这些模式的形成也同样受各种因素的影响。利用孔隙介质中油运移实验结果总结出来的、以Bo数和Ca数为纵横坐标运移相图完全可以用来表征单个裂隙中油的运移模式。     

The loess sequence of Dolní Věstonice,Czech Republic: A new OSL‐based chronology of the Last Climatic Cycle

The Dolní Věstonice loess section in the Czech Republic is well known for its high‐resolution loess–palaeosol sequence of the last interglacial–glacial climatic cycle (Upper Pleistocene). The loess section is situated in a climatic transition zone between oceanic and continental climates and is therefore of great value in reconstructing past regional climate conditions and their interaction with climate systems, in particular that of the North Atlantic. Based on a combination of optically stimulated luminescence (OSL) ages, stratigraphic field observations and magnetic susceptibility measurements, a chrono‐climatic interpretation of the Dolní Věstonice loess section is presented. To establish a reliable Upper Pleistocene chronology, a quartz OSL approach was applied for equivalent dose (D_e) determination. Monomineralic quartz extracts of three distinct grain sizes, fine (4–11 μm), middle (38–63 μm) and coarse (90–200 μm), were used and compared. Within error limits, the calculated OSL ages are the same for the different grain sizes, and the OSL ages are in stratigraphic order. The established OSL chronology is in agreement with a Weichselian litho‐ and pedostratigraphy. The Dolní Věstonice loess section is characterized by four pedosedimentary subsequences. At the base of the profile, subsequence I is characterized by a distinct Early Glacial soil complex, OSL‐dated to c. 110 to 70 ka, representing one of the most complete records of environmental change in the European loess belt. Subsequence II is allocated to the Lower Pleniglacial and is characterized by laminated sandy loess. Middle Pleniglacial subsequence III is represented by a brown soil complex, and is followed by the uppermost subsequence IV, characterized by a thick body of laminated sandy loess, indicating strong wind activity and a high sedimentation rate of more than ～1 mm a^?1 during the Upper Pleniglacial. According to the OSL chronology, as well as to the sedimentological and palaeopedological investigations, it is likely that the sequence at Dolní Věstonice has recorded most of the climatic events expressed in the NGRIP δ¹⁸O reference record between 110 and 70 ka.