A three-dimensional model of the Pyrenees and their foreland basins from geological and gravimetric data

We construct a three-dimensional geological model of the Pyrenees and their foreland basins with the Geomodeller. This model, which accounts for different sources of geological and geophysical informations, covers the whole Pyrenees, from the Atlantic Ocean to the Mediterranean Sea, and from the Iberian range to the Massif Central, down to 70 km depth. We model the geological structure with a stratigraphic column composed of a superposition of layers representing the mantle, lower, middle, and upper crusts. The sedimentary basins are described by two layers which allow us to make the distinction between Mesozoic and Cenozoic sediments, which are characterized by markedly different densities and seismic velocities. Since the Pyrenees result from the convergence between the Iberian and European plates, we ascribe to each plate its own stratigraphic column in order to be able to model the imbrication of Iberian and European crusts along this fossile plate boundary. We also introduce two additional units which describe the orogenic prism and the water column in the Bay of Biscay and in the Mediterranean Sea. The last ingredient is a unit that represents bodies of shallow exhumed and partly serpentinized lithospheric mantle, which are assumed to produce the positive Bouguer gravity anomalies in the North Pyrenean Zone. A first 3D model is built using only the geological information coming from geological maps, drill-holes, and seismic sections. We use the potential field method implemented in Geomodeller to interpolate these geological data. This model is then refined in order to better explain the observed Bouguer anomalies by adding new constraints on the main crustal interfaces. The final model explains the observed Bouguer anomalies with a standard deviation less than 3.4 mGal, and reveals anomalous deep structures beneath the eastern Pyrenees.