The Roussillon Basin is a non-silled Miocene sedimentary basin filling a late Oligocene-early Miocene graben. The basin was intensively impacted by the 2 Messinian fluvial erosion, as evidenced in exposed sections, in seismic profiles and in deep boreholes drilled for hydrocarbon exploration. As the basin was open to the Mediterranean Sea, the huge drop in sea level at the peak of the Messinian Salinity Crisis is clearly recorded, along with the subsequent sudden marine reflooding and the resulting prograding sedimentary filling, particularly in Gilbert-type fan deltas. Here, the Messinian Erosional Surface (MES) is accurately mapped in a high-resolution document, which corrects the confusion resulting from the set of 1:50,000 scale regional maps. Aim of the 3D reconstruction of the MES is to modernize geological mapping, a crucial challenge for Mediterranean and peripheral areas. Thanks to a reliable chronostratigraphy provided by planktonic foraminifers, calcareous nannofossils, micro-and macro-mammal remains, paleomagnetism and a 10 Be cosmogenic nuclide-derived study, our reconstruction is one of the most comprehensive models of changes in sea level from 6 to 3 Ma. After the marine reflooding of the Mediterranean Basin at 5.46 Ma, the fluctuations in sea level recorded in the Roussillon Basin were forced by global changes. Following reflooding, the Prades large olistostrome collapsed prepared by the previous exhumation along the Canigou fault. The olistostrome is a good example of a local accident resulting from Messinian events. The exceptional changes in sea level at the peak of the Messinian Salinity Crisis deeply marked the Roussillon Basin, momentarily overprinting the Pyrenean orogenesis.
A thermokarst is a collapse feature resulting from the thawing of ice-rich permafrost or of massive ice of various origins. Little attention has been paid to the sedimentary fabric resulting from this type of collapse, except for glaciotectonic features. In western Europe, two palaeo-forms are commonly studied: lithalsas and ice-wedge casts. Collapsed pingos are much rarer. Very few papers have compiled present-day and fossil data. Here, field data collected from quarries in the eastern Paris Basin were analysed, providing useful records of thermokarst collapses in alluvial calcareous silts, sands, and gravels. These forms have a circular shape when viewed on satellite images. Permafrost is attested regionally by the recurrent occurrence of meter-sized pattern grounds at the surface of the chalk and of ice-wedge casts. Traces of segregation and reticulate ice are common. These features are primarily connected to a major interstadial, c. 150 ka BP, orbitally forced and commonly associated with a major glacial retreat. They occur both in drained and waterlogged situations, resulting in a specific pattern of deformation. They are controlled by the brittle and plastic behaviour of sediments and resemble passive glaciotectonism. Normal and reverse faults, with the offset decreasing downward, are common, and those with local shear are reported. Lithalsas, seasonal frost blisters, spring frost blisters and perhaps pingos seem to have formed. Most of these deformations correspond to thermokarst sinkholes bordered by gravitational collapse faults. The offset of these faults increases towards the surface, and the faults have been recurrently confused with neotectonism triggered by palaeo-earthquakes. However, there are no faults beneath the observed deformation features, and the region lacks recorded seismic activity over the last century. Our data may be helpful in interpreting similar structures elsewhere.Brigitte Van Vliet-Lano€ e (brigitte.vanvlietlanoe@univ-brest.fr), Domaines Oc eaniques Laboratory, UMR 6538
Abstract. This study presents a novel workflow to model the internal heterogeneity of complex aquifers using the multiple-point statistics algorithm DeeSse. We illustrate the applicability of this workflow on the Roussillon's aquifer in the region of Perpignan (southern France). This work is part of a project aiming at assessing the groundwater dynamics of this Mediterranean aquifer in the context of a growing population, climate change, and increasing pressure on the freshwater resources. We focus here on the geological heterogeneity of the Continental Pliocene layer because it is expected to influence possible saltwater intrusion process and its corresponding uncertainty quantification. The main aim of the paper is therefore to describe the procedure that is used to model the aquifer heterogeneity with a relatively small number of direct geological observations and a well defined geological concept. When few direct observations are available, the traditional geostatistical approaches cannot be applied easily because variogram inference is difficult. On the opposite, multiple-point statistics simulations can rely on a conceptual geological model. Here, the conceptual model consists not only of a training image displaying the spatial organization of the main sedimentological elements in space, but also in a set of additional information such as general trends and paleo orientations of the sedimentological features. The direct sampling algorithm DeeSse can be used in this context to model the expected heterogeneity. The workflow involves creating 2D non-stationary training images (TI) coupled during simulation with auxiliary information and controlled by hard conditioning data obtained from interpreted electrofacies. To control the non-stationarity, a 3D trend map is obtained by solving numerically the diffusivity equation as a proxy to describe the spatial evolution of the sedimentary patterns, from the source of the sediments to the outlet of the system. A 3D continuous rotation map is estimated from paleo orientations of the fluvial system. Both trend and orientation maps are derived from geological insights gathered from outcrops and general knowledge of processes occurring in these types of sedimentary environments. Finally, the 3D model is obtained by stacking 2D simulations following the paleo-topography of the aquifer. The vertical facies transition between two 2D simulations is controlled by both the hard conditioning data set and by simulating conditional data points from one simulation to another. This process allows to bypass the creation of a 3D training image while preserving the vertical continuity of the sedimentary objects.
Abstract. This study introduces a novel workflow to model the heterogeneity of complex aquifers using the multiple-point statistics algorithm DeeSse. We illustrate the approach by modeling the Continental Pliocene layer of the Roussillon aquifer in the region of Perpignan (southern France). When few direct observations are available, statistical inference from field data is difficult if not impossible and traditional geostatistical approaches cannot be applied directly. By contrast, multiple-point statistics simulations can rely on one or several alternative conceptual geological models provided using training images (TIs). But since the spatial arrangement of geological structures is often non-stationary and complex, there is a need for methods that allow to describe and account for the non-stationarity in a simple but efficient manner. The main aim of this paper is therefore to propose a workflow, based on the direct sampling algorithm DeeSse, for these situations. The conceptual model is provided by the geologist as a 2D non-stationary training image in map view displaying the possible organization of the geological structures and their spatial evolution. To control the non-stationarity, a 3D trend map is obtained by solving numerically the diffusivity equation as a proxy to describe the spatial evolution of the sedimentary patterns, from the sources of the sediments to the outlet of the system. A 3D continuous rotation map is estimated from inferred paleo-orientations of the fluvial system. Both trend and orientation maps are derived from geological insights gathered from outcrops and general knowledge of processes occurring in these types of sedimentary environments. Finally, the 3D model is obtained by stacking 2D simulations following the paleo-topography of the aquifer. The vertical facies transition between successive 2D simulations is controlled partly by the borehole data used for conditioning and by a sampling strategy. This strategy accounts for vertical probability of transitions, which are derived from the borehole observations, and works by simulating a set of conditional data points from one layer to the next. This process allows us to bypass the creation of a 3D training image, which may be cumbersome, while honoring the observed vertical continuity.
<p>More than 80 million m3 per year are pumped into the Roussillon plain coastal aquifer, covering 850 km&#178; and located between the Pyrenean massif to the west and the Mediterranean Sea to the east, south of France. This is a multilayer aquifer of more than 350 m thick, made up of sandy layers embedded in low-permeability clayey material from the Pliocene and topped by alluvial formations from the Quaternary. Its groundwater resource is primarily used for the supply of drinking water, but also contributes to the irrigation of some 13,000 hectares.</p><p>For more than 40 years, this aquifer has been undergoing a general decline in its piezometric level due to pumping and water demand is expected to increase (growing irrigation areas and climatic demand). Moreover, given its flat topography, the Roussillon plain is likely to suffer sea water intrusions and marine submersion, due to the sea level rise, which could reach 1 m by 2100.</p><p>This context shaped the Dem'Eaux Roussillon project, which brought together nearly ten partners from the Occitanie region (research units, consultancies and local authorities). Its objective was to characterise the behaviour of the groundwater resource in this aquifer, in order to be able to project its future situation, in the context of climate change, rising sea levels (risk of saline intrusion) and changes in water use. A detailed characterisation of the geological reservoir highlighted the need to consider the offshore extension of this coastal aquifer. The analysis of the piezometric evolution at the scale of the Roussillon plain over the last 50 years allowed the spatialized characterization of the hydrodynamic parameters and the understanding of the vertical drainage processes that control the hydraulic equilibrium between the Quaternary and the Pliocene water tables. Two high-resolution hydro-geophysical observatories have been set up to quantify these processes and improve understanding of saline intrusions processes. Finally, a conceptual model presenting the main features of the main processes controlling the groundwater evolution and the sea water intrusion risk was obtained ready to launch a numerical modelling work.</p>
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.