# Integrated Paleoenvironmental Reconstruction and Taphonomy of a Unique Upper Cretaceous Vertebrate-Bearing Locality (Velaux, Southeastern France)

Aude Cincotta, Johan Yans, Pascal Godefroit, Géraldine Garcia, Jean Dejax, Mouloud Benammi, Sauveur Amico, Xavier Valentin

► To cite this version:

Aude Cincotta, Johan Yans, Pascal Godefroit, Géraldine Garcia, Jean Dejax, et al. Integrated Paleoenvironmental Reconstruction and Taphonomy of a Unique Upper Cretaceous Vertebrate-Bearing Locality (Velaux, Southeastern France). PLoS ONE, 2015, 10 (8), pp.e0134231. 10.1371/journal.pone.0134231. hal-01232372

HAL Id: hal-01232372

https://hal.sorbonne-universite.fr/hal-01232372

Submitted on 23 Nov 2015

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Integrated Paleoenvironmental Reconstruction and Taphonomy of a Unique Upper Cretaceous Vertebrate-Bearing Locality (Velaux, Southeastern France)

Aude Cincotta¹,²*, Johan Yans¹, Pascal Godefroit², Géraldine Garcia³, Jean Dejax⁴, Mouloud Benammi³, Sauveur Amico⁵, Xavier Valentin³,⁶

1 Department of Geology, NaGRIDD, University of Namur, Namur, Belgium, 2 Directorate 'Earth and History of Life', Royal Belgian Institute of Natural Sciences, Brussels, Belgium, 3 Institut de Paléoprimatologie, Paléontologie Humaine: Evolution et Paléoenvironnements UMR 7262, University of Poitiers TSA 51106, Poitiers, France, 4 Centre de Recherches sur la Paléobiodiversité et les Paléoenvironnements (CR2P, UMR 7207), Sorbonne Universités- MNHN, CNRS, UPMC Paris 6. Muséum national d'Histoire naturelle, Paris. France, 5 Hôtel du Département, Direction de l'Environnement, Conseil Général des Bouches-du-Rhône, Marseille, France, 6 Palaios association, Valdivienne, France

* aude.cincotta@unamur.be

<figure><img src="image_2.png" /></figure>

CrossMark
click for updates

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Cincotta A, Yans J, Godefroit P, Garcia G, Dejax J, Benammi M, et al. (2015) Integrated Paleoenvironmental Reconstruction and Taphonomy of a Unique Upper Cretaceous Vertebrate-Bearing Locality (Velaux, Southeastern France). PLoS ONE 10(8): e0134231. doi:10.1371/journal.pone.0134231

**Editor:** Peter Dodson, University of Pennsylvania, UNITED STATES

**Received:** January 30, 2015

**Accepted:** July 7, 2015

**Published:** August 19, 2015

**Copyright:** © 2015 Cincotta et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper.

**Funding:** AC disposes of a FRIA grant provided by the 'Fond National pour la Recherche Scientifique' (FNRS). The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

## Abstract

The Velaux-La Bastide Neuve fossil-bearing site (Bouches-du-Rhône, France) has yielded a diverse vertebrate assemblage dominated by dinosaurs, including the titanosaur *Atsinga-nosaurus velauciensis*. We here provide a complete inventory of vertebrate fossils collected during two large-scale field campaigns. Numerous crocodilian teeth occur together with complete skulls. Pterosaur, hybodont shark and fish elements are also represented but uncommon. Magnetostratigraphic analyses associated with biostratigraphic data from dinosaur eggshell and charophytes suggest a Late Campanian age for the locality. Lithologic and taphonomic studies, associated with microfacies and palynofacies analyses, indicate a fluvial setting of moderate energy with broad floodplain. Palynomorphs are quite rare; only three taxa of pollen grains occur: a bisaccate taxon, a second form probably belonging to the Normapolles complex, and another tricolporate taxon. Despite the good state of preservation, these taxa are generally difficult to identify, since they are scarce and have a very minute size. Most of the vertebrate remains are well preserved and suggest transport of the carcasses over short distances before accumulation in channel and overbank facies, together with reworked Aptian grains of glauconite, followed by a rapid burial. The bones accumulated in three thin layers that differ by their depositional modes and their taphonomic histories. Numerous calcareous and iron oxides-rich paleosols developed on the floodplain, suggesting an alternating dry and humid climate in the region during the Late Campanian.
**Introduction**

Late Cretaceous continental deposits are widely exposed in southern France and have yielded numerous and diverse vertebrate remains [1] especially in the Aix-en-Provence Basin (e.g. [2], [3], and [4]). However, few studies of both sedimentology (including lithofacies, microfacies, and palynofacies) and vertebrate taphonomy have been conducted at these localities (e.g. [5]). The Velaux-La Bastide Neuve site (Fig 1) was discovered in 1992 by one of us (X.V.) and hundreds of vertebrate remains have been collected during an initial survey (2002) and two large-scale field campaigns (2009 and 2012). The vertebrate assemblage represents a highly diverse fauna including chelonian, crocodilian, dinosaur and pterosaurs. Among the dinosaurs, a new titanosaur genus, *Atsinganosaurus velauciensis*, was described based on partially articulated skeletons [6]. The present paper aims to describe the sedimentological and taphonomic context of the Velaux locality and to reconstruct its paleoenvironment, using lithofacies, microfacies, palynofacies associated with taphonomy and fossil descriptions, and also provide a complete inventory of the vertebrate taxa found in this area.

**Geological Setting and Stratigraphy**

The Velaux-La Bastide Neuve fossil-bearing site is located in the western part of the Aix-en-Provence Basin (Bouches-du-Rhône department, southeastern France), an east-west oriented syncline of about 400 km² [7]. The basin fill is composed of fluvio-lacustrine deposits dating from Santonian to Lutetian age [8] that were deposited following an episode of epeirogeny at the end of the Santonian [9]. The age of the site is based on continental biostratigraphic data such as charophytes and dinosaur eggshell. The continental layers exposed in Velaux-La Bastide Neuve site were previously attributed to the “Begudian” local stage (Fig 1), correlated to the Late Campanian [10]. Two charophyte biozones have been recognized at Velaux-La Bastide Neuve. The *Peckichara pectinata* biozone is correlated to the middle-early Late Campanian,

<figure><img src="image_4.png" /><figcaption>**Fig 1.** Geographical and Geological maps of Velaux-La Bastide Neuve (VBN) area. VBN is located in southeastern France, between Marseille and Aix-en-Provence. Our fossil-bearing locality is indicated by the bone, between two major roads (D113 and A7). The site belongs to the "Begudian" local stage.</figcaption></figure>

doi:10.1371/journal.pone.0134231.g001
and the *Peckichara cancellata* biozone to the Late Campanian [5]. In addition, dinosaur egg-
shell from this locality belong to the *Megaloolithus aureliensis* biozone, correlated to the Late
Campanian [10].

**Material and Methods**

During the two field campaigns a surface of 375 m² to a depth of 1.2 m was excavated, which resulted in 100 m³ of overburden and matrix. During the first campaign (2009), several small areas separated from the main fossil-bearing layer were worked. In 2012, two large sections were excavated to complete the previous sampling. A total of 308 fossil specimens were inventoried during the two field campaigns- they are housed at the Moulin seigneurial Museum and at the Henri-Ricard Archaeological repository, both located in Velaux.

Rock samples for sedimentological analyses were collected from each layer of the sedimentological section (Fig 2). Sampling for several types of analysis was performed at the same time; 18 samples were collected for magnetostratigraphic analyses. Forty one samples were used to study the lithofacies; among these, eight were used for palynological processing, and thirteen were used to prepare thin sections for microfacies analysis. Palynological samples were prepared using standard methods [11]. Two samples (VBN-18 and VBN-20A) collected from sandy horizons were used for isotopic dating using standard K-Ar methods [12].

Magnetostratigraphy

Oriented block samples were collected from 27 stratigraphic levels in the 16.3 m thick section.
The paleomagnetic analyses were carried out at the Laboratory of the iPHEP (University of
Poitiers). The intensity and direction of remanent magnetization were measured with a JR6
spinner magnetometer using four position standard specimen holders.

Taphonomy

The position of each bone collected during the 2009 and 2012 field seasons was mapped (Fig 3). Data were also collected as to the nature, dimensions and orientations in the case of elongate bones. The fossil-bearing area covers a surface of about 140 m². Features of the bones (weathering stage, fractures) were studied in laboratory (University of Poitiers and IRSNB), in order to document the taphonomy of the locality. Stereonet free program was used for drawing bidirectional rose diagrams and the nonparametric Rayleigh's test was performed for testing uniformity of the data.

Results

Magnetostratigraphy

Concerning the magnetostratigraphy, only one paleomagnetic component could be recognized in the studied samples (Fig 4a and 4b). The greater part of remanent magnetization was removed at 60 mT, which may indicate that magnetite is the main carrier. The average directions were determined after tilt correction (Fig 4c): the mean direction of the characteristic components is (I = 56.7, D = 4.42, a = 5.4, kappa = 42; n = 18). The mean direction is different from the present day magnetic field at the sampling site suggesting a primary remanence of magnetization. The rocks exhibit a single normal polarity. According to the biochronological data, the correlation with the Geomagnetic Polarity Time Scale (GPTS) could be done with the normal chron of the chron C32 (C32n.1n, C32n.2n or C32r.1n). Accordingly, an age of between 70.9 and 73.3 Ma, or an adjusted age of between 71.6 and 74.0 Ma [12], can be assigned to the Velaux-La Bastide Neuve section.
<figure><img src="image_5.png" /><figcaption>**Fig 2. Sedimentological succession of Velaux-La Bastide Neuve locality.** Three vertebrate-bearing layers are present: in the sequences 1, 2 and 3. Numbers on the left side correspond to samplings; blue numbers are palynofacies samples; red numbers are microfacies samples and black ones are lithofacies samples. The three tables on the right discriminate the different facies for these samples. Location of magnetostratigraphic samples are indicated by arrows, and the letters MA.</figcaption></figure>

doi:10.1371/journal.pone.0134231.g002
<figure><img src="image_6.png" /><figcaption>**Fig 3.** Map of the fossil-bearing site at Velaux-La Bastide Neuve. Two zones were explored, during the two field campaigns (2009 and 2012) on a total surface of about 140 m². The three taphonomic modes (TM) are represented by three colors—TM1 in green, TM2 in black, and TM3 in red—, which represent the three successive deposits.</figcaption></figure>

doi:10.1371/journal.pone.0134231.g003

## Section, Sequences, and Facies

The studied section has a thickness of 16.3 m from the basal level located below the lowermost bone-bed to the top of the upper limestone bed. The section can be divided into nine depositional sequences (Fig 2). The lower contact of each sequence is erosional and each sequence fines upward. Four sequences are characterized by the same facies succession: sandstones at the base followed by calcareous siltstones. Only the uppermost sequence contains limestones, indicating a different sub-environment.

Oxidized glauconite grains were found in sandstones from sequence 3 (samples VBN-18 and VBN-20A) and were dated using K-Ar analysis. Results provide ages of 122.2 ± 3.2 Ma and 123.6 ± 3.4 Ma (Table 1), corresponding to the early Aptian stage [12], [13].

The entire succession comprises five dominant lithofacies: sandstones (Fig 5a and 5b), siltstones (Fig 5c and 5d), limestones (Fig 5e and 5f), mudstones (Fig 5f) and lignite layers. The sedimentological section is dominated by variegated siltstones (52% of the section thickness) and sandstones (34%; Fig 2). There are five microfacies: bioclastic micro-conglomerate (Fig 6a),

<figure><img src="image_7.png" /><figcaption>**Fig 4.** Zijderveld diagrams. (a) and (b) representative Zijderveld diagrams of the samples (c) equal-area stereographic projection of characteristic direction.</figcaption></figure>

doi:10.1371/journal.pone.0134231.g004
<table><caption>Table 1. Radiometric results on glauconite grains belonging to two sandstone samples.</caption><thead><tr><th>Sample</th><th>K (%)</th><th><sup>40</sup>Ar rad (nl/g)</th><th><sup>40</sup>Ar air (%)</th><th>Age (Ma)</th></tr></thead><tbody><tr><td>VBN-18</td><td>1.68</td><td>8.073</td><td>37.5</td><td>122.2 ± 3.2</td></tr><tr><td>VBN-20A</td><td>1.69</td><td>8.265</td><td>48.3</td><td>123.6 ± 3.4</td></tr></tbody></table>

doi:10.1371/journal.pone.0134231.t001

bioclastic sandstone (Fig 6b), bioclastic siltstone (Fig 6c and 6d), micritic wackestone (Fig 6e) and azoic oxidized siltstone (Fig 6f). Palynological analysis allow discrimination of three palynofacies: a high diversity palynofacies (Fig 7a and 7b), a charcoal-rich palynofacies (Fig 7c and 7d) and a phytoclast-rich palynofacies (Fig 7e and 7f). Lithofacies are associated with microfacies and palynofacies to specify the nature of sub-environments. Nonetheless, some of the same microfacies and palynofacies are associated with different, but related, sub-environments.

**Facies 1.** Facies 1 corresponds to predominantly sandstone deposits that in all cases overlie an erosional base (Fig 5a). The facies is composed of yellow to grey sandstones and displays various grain sizes, from conglomeratic to coarse-grained sandstones. Bed thickness varies from 5 cm to 80 cm. Sandstones exhibit small-scale planar cross-stratification and horizontal lamination. At the base of sequence 7 (Fig 2), a thin clast-supported intraformational conglomerate is present. Sandstones are in some cases interstratified with thin layers of lignite. Sandstone abundance decreases from the base of the section to the top. Thin black lignite laminae are interstratified within sandy beds, especially at the top of sequences 2 and 3 (Fig 2). Fossil trunk wood traces are also found in these layers. Facies 1 comprises two distinct microfacies. The first one is the coarsest, represented by a micro-conglomerate composed of subangular to subrounded quartz and feldspar grains. It contains micritic grains as algal balls, glaebules, peloids, calcareous intraclasts (Fig 6a), fragments of mollusk shells (Fig 6b), “algae” and foraminifera. The second microfacies is a bioclastic sandstone that contains the same skeletal elements but is matrix-supported. The palynological assemblage comprises abundant phytoclasts (rounded charcoal grains and various lignitic debris), bisaccate pollen grains, fungal spores (Fig 7a) and fungal hyphae (Fig 8a).

**Facies 2.** The deposits of facies 2 are characterized by medium to coarse-grained sandstones (Fig 5b). The sandstones are in some cases variegated, contain mollusk shells, calcareous nodules and, locally, orange mud clasts. Sediments are mixed with reworked debris, including marine invertebrate fossils (foraminifera and “red algae”). Feldspar grains, foraminifera, *Microcodium* (calcite prisms induced by biogenic processes [14], [15]), freshwater fossils as Characeae gyrogonites, algal balls, mud coated grains and stromatolites are also present in these deposits. Skeletal grains are in some cases highly micritized and not easily identifiable (Fig 6c). The palynofacies contains the same palynomorphs as facies 1; including phytoclasts, abundant charcoal and bisaccate pollen grains (Fig 8b), fungal hyphae and fungal spores.

**Facies 3.** The third facies consists of variegated and fine-grained siltstones. The siltstones are bioclastic—and contain mollusk shells, foraminifera, *Microcodium*, algal balls, mud-coated grains and stromatolites (Fig 6d). Skeletal grains are highly micritized. The palynological material constitutes mostly phytoclasts, mainly with charcoal and lignitic debris (Fig 7e); but fungal hyphae are also present. Pollen grains (Fig 7f) are rare. Charcoal grains are less abundant than in the facies 1. Tracheid fragments are bigger and better preserved than charcoal.

**Facies 4.** Sediments of facies 4 are variegated, but mostly reddish and fine-grained siltstones. No sedimentary structures are present in these siltstones, but pedogenic features, such as carbonate nodules are well developed. Skeletal elements are also lacking (Fig 6f). There are only rare reworked glauconite grains, and mineral grains are coated by iron oxides. The corresponding palynofacies contains nearly the same organic particles as facies 1 and 2, although
<figure><img src="image_8.png" /><figcaption>**Fig 5. Lithofacies.** (a) small-scale cross-laminated sandstones (seq. 1) (b) beige fine-grained sandstones showing no sedimentary features (seq. 3) (c) variegated siltstones overlying grey sandstones (seq. 6) (d) blue-grey siltstones overlaying variegated siltstones (seq. 7) (e) both first beds from the massive limestone of the sequence 8 (f) thin muddy layers below the lower limestone bed, and the undulated limestone layer at the base of the lower bed (seq. 8).</figcaption></figure>

doi:10.1371/journal.pone.0134231.g005
<figure><img src="image_9.png" /><figcaption>**Fig 6. Microfacies.** (a) MF1—sandy conglomerate including mollusk valves and micritic grains (b) MF1—carbonated intraclast including *Microcodium* (c) MF2—fine-grained sandstone with algal ball and mollusk shells (d) MF3—coarse-grained siltstone including micritic grains and mollusk shells (e) MF4—micritic limestone including *Microcodium* and shells (f) MF5—fine-grained siltstone coated with iron oxides. Scale bars = 500 µm.</figcaption></figure>

doi:10.1371/journal.pone.0134231.g006
<figure><img src="image_10.png" /><figcaption>**Fig 7.** Palynofacies preparations from the Velaux-La Bastide Neuve sedimentological succession. (a) PF1- high diversity section with lignitic debris, charcoal and abundant fungal hyphae (scale bar = 100 µm) (b) PF1- fungal hyphae showing septae, lignitic debris, charcoal and amorphous organic matter (scale bar = 50 µm) (c) PF2- numerous rectangular and rounded charcoal grains, and lignitic debris (scale bar = 50 µm) (d) PF2- A.O.M., charcoal and pollen grain (scale bar = 50 µm) (e) PF3- charcoal and tracheid fragments (f) PF3- fungal spore, diverse plant debris (e.g. tracheid fragments), charcoal and A.O.M. (scale bar = 50 µm).</figcaption></figure>

doi:10.1371/journal.pone.0134231.g007
<figure><img src="image_11.png" /><figcaption>**Fig 8.** Fungal hyphae, spores and palynomorphs. (a) fungal spore with thin wall and folds and connected to an hypha (*Paleomycites*) [12-VBN-05; sequence 1] (b) bisaccate pollen grain among coalified plant debris and hyphae (arrow) [12-VBN-05; sequence 1] (c) prolate tricolporate pollen grain in meridian view (d) pollen grain probably belonging to the Normapolles complex (polar view) (e) monoaperturate pollen grain, probably monosulcate pollen grain (f) semitectate, reticulate, probably monosulcate pollen grain (g) and (h) monosulcate pollen grain; same taxon in median (g) and high (h) focus.</figcaption></figure>

doi:10.1371/journal.pone.0134231.g008
fungal hyphae (Fig 7c) are absent. Charcoal grains are very abundant and lignitic debris is less common although less degraded than in facies 3. A few bisaccate (Fig 8b), tricolporate (Fig 8c) and monosulcate (Fig 8e) pollen grains, and also some Normapolles (Fig 8d) occur.

**Facies 5.** The fifth facies is only present in sequence 8 (Fig 2) and consists of resistant limestone beds interstratified with thin greyish mudstones. The limestones reach 1.2 m in thickness, and include very fine-grained beds containing mollusk shells and ostracod valves. Undulatory laminae are present at the base of the lowermost bed (Fig 5f). Upper beds are massive without any sedimentary structure. Very few mudstone layers occur in this facies; where present, they are horizontally laminated. Limestones are overlain by reddish siltstones presenting color motting. The microfacies is a fine micritic limestone, containing abundant *Microcodium* (Fig 6e) but also charophyte gyrogonites and some ostracods. *Microcodium* are in some cases wrapped in thin mud layers. No detrital grains are present in this facies, but only a fine micritic mud with a wackestone texture including *Microcodium* features and shells dispersed in a micritic matrix.

## Taphonomy

**Taxonomic and skeletal composition.** Three hundred and eight vertebrate fossils—including bones, plates, osteoderms, and teeth—belonging to several taxa were inventoried at Velaux during the two field-campaigns. Dinosaur remains dominate the vertebrate fossils (Fig 9—above), representing over 38% of the total assemblage collected during the 2009 and 2012 campaigns. Crocodilian elements are particularly abundant (16%), and mainly represented by shed teeth (Fig 9—below right). Numerous crocodilian (non-shed) teeth occur together with incomplete juvenile skulls and one complete adult skull (MMS/VBN-12-10A). Chelonian carapace and plastron parts are abundant (64 specimens collected so far), representing 22% of the assemblage. Pterosaur, hybodont shark and fish elements are also represented in the assemblage although they are rather rare (17, 10 and 4 elements, respectively). Among the identifiable dinosaur remains, ornithopod (39%) elements are the most abundant. Various ankylosaur (24.5%) elements were found including osteoderms, long bones and one scapula. Titanosaur elements (22.4%) comprise vertebrae and long bones. Theropods (14%) are mainly known

<figure><img src="image_12.png" /><figcaption>**Fig 9. Vertebrate remains distribution and bone nature.** Dinosaurs and chelonians dominate the fossil assemblage. Elements as teeth or plates are particularly abundant. Data were taken during two field campaigns, in 2009 and 2012.</figcaption></figure>

doi:10.1371/journal.pone.0134231.g009
from isolated teeth, with few long bones and vertebrae. Small fossils are numerous, dominated
by teeth and chelonian shell parts. Fragile bones such as ribs, chevrons or cranial elements are
under-represented (Fig 9—below left).

**Accumulation modes.** The vertebrate fossils of the Velaux-La Bastide Neuve site occur in three thin layers with different accumulation modes reflecting different depositional histories. Investigations were concentrated on the three vertebrate-bearing layers (Fig 3), located in the lower part of the sedimentological section in sequences 1, 2 and 3 (Fig 2). Bones are particularly abundant in the second fossiliferous layer, which consists of fine to coarse-grained sandstone.

Isolated and disarticulated elements belonging to a single taxon were found in variegated coarse-grained siltstones from sequence 1, and were interpreted as an overbank facies. Two long (>2m) fragile but very well preserved parallel bones were found in 2009 in these variegated siltstones, in association with skin impressions. These elements can be identified as sauropod cervical ribs (Koen Stein, pers. comm.). The lowermost fossiliferous layer cannot be regarded as a bone bed, in the current state of excavation, as all the elements possibly belong to a single individual [16].

The second fossiliferous layer is represented by a multitaxic accumulation of elements included in coarse-grained conglomeratic sandstones from sequence 2. The following elements are represented in this multitaxic bone bed; ornithopod bones and teeth, bones and teeth from the titanosaur *Atsinganosaurus velauciensis* [6], theropod teeth, crocodilian teeth, isolated bones, skulls (Fig 10a) and one osteoderm, chelonians carapace and plastron elements and long bones belonging to two different taxa—*Solemys* (Fig 10b) and *Polysternon*—, hybodont shark teeth, fish teeth and pterosaur bones. Most of bones are disarticulated, except partially articulated titanosaur skeletons [6]. Elements, presenting various sizes and shapes, are mixed together and are often entangled in the sandy matrix. Long bones in this layer tend to be aligned along a NW-SE axis (Fig 11). Predation marks are extremely rare, they are even absent on titanosaur bones, although crocodilian shed teeth are rather abundant. A single deep tooth mark is present near the ocular fenestrae of a complete adult crocodilian skull discovered in 2012 in this bone bed (Fig 10a). Trample marks have not been identified in this material; these are characterized by shallow, sub-parallel scratch striae on the bone surfaces [17], [18].

A third taphonomic mode consists of disarticulated bones (Fig 10d and 10e), plates, and abundant teeth found in sequence 3. The bones are more dispersed than in the other modes and belong to several taxa. Ankylosaur long bones and osteoderms are only present in this silty sandstone layer. Theropod teeth, ornithopod teeth, crocodilian osteoderms and teeth, and chelonian plates belonging to both *Solemys* and *Polysternon* (Fig 10c) occur. Neither predation marks nor tramples marks are present. The chelonian remains are very well preserved and usually less degraded than dinosaur bones. Long bones (mostly belonging to the ankylosaur) are mostly oriented NE-SW but other elements present orientations that are quite varied (Fig 11).

Discussion

Paleoenvironment: integrated litho-, micro-, and palynofacies

The observed lithologies, microfacies and palynomorphs documented in the nine depositional sequences exposed at Velaux-La Bastide Neuve indicate several alluvial sub-environments, ranging from channel lag to floodplain.

**Facies 1.** Conglomeratic beds (e.g. sequence 2) correspond to the basal deposits in a stream channel, whereas all finer sandstones are interpreted as channel levee deposits or crevasse splay deposits. The presence of small-scale planar cross-stratification indicates low-energy sedimentation [19], [20], whereas horizontal bedding corresponds to upper flow regime conditions [19], [21]. Mud clasts could indicate reworking of channel banks within the river [22] or basal
<figure><img src="image_13.png" /><figcaption>**Fig 10.** Plate illustrating the diversity of vertebrate remains collected from the two bone beds of Velaux-La Bastide Neuve. (a) complete crocodilian skull (MMS/VBN-12-10A) from an adult specimen that shows a deep tooth pit near the ocular fenestra (square) (b) ornamented chelonian plate attributed to the genus *Solemys* (c) non-ornamented chelonian plate attributed to the genus *Polysternon* (d) indetermined dinosaur phalanx (e) dinosaur vertebra.</figcaption></figure>

doi:10.1371/journal.pone.0134231.g010

scour during channel incision. Lignite layers are interpreted as organic matter accumulations introduced by bank erosion into the river channel, collected by the river stream and then deposited along with bedload sand during lower flow regime conditions [22]. Glauconite grains dated from 122.2 ± 3.2 Ma and 123.6 ± 3.4 Ma are interpreted as reworked grains probably originating from Aptian marine limestones, widespread in this area (e.g. [23]). Foraminifera observed in thin sections are also probably reworked from older carbonate rocks. *Microcodium* [15] may also be reworked. Broken and rounded grains occur in the upper flow regime sediments, resulting from deposition in the channel or in a crevasse play. Phytoclasts and in
<figure><img src="image_14.png" /><figcaption>**Fig 11.** Rose diagrams showing the long-axis orientation for bones collected at Velaux-La Bastide Neuve. (a) represents the south-east/north-west trend for long bones from the second taphonomic layer (BB1; n = 23); and (b) shows the orientation pattern for long bones of the third taphonomic layer (BB2; n = 12). Statistical analyses show that the null hypothesis of uniformly distributed data can be rejected for (a) but cannot be rejected for (b); Rayleigh's test, (b) z = 2.969, p < 0.05; (b) z = 0.661, p > 0.05.</figcaption></figure>

doi:10.1371/journal.pone.0134231.g011

particular lignitic debris reflect the presence of vascular plants in the environment, either from riparian forests (if found in channel deposits) or other remote forests. The roundness of most charcoal indicates a prolonged transportation of these particles in the river channel. Fungal hyphae indicate recycling of wood debris from surrounding soils [24].

**Facies 2.** Facies 2 is interpreted as levee deposits. Because feldspar crystals are unstable and quickly weathered, their presence in the sandstones reflects rapid transport and deposition from the sediment source (e.g. [25]). The micritization of many skeletal elements reveals shallow water conditions [26], [27]. The rounded shape of charcoal still indicates a prolonged transportation from the plant source. These charcoal grains potentially reflect different stages of organic matter maturity. Indeed, they are common features in palynological assemblages, corresponding to charring states of diverse organic debris; they may be a result of wild fires. Vegetation in levee facies was likely abundant, suggesting rare inundation episodes [28].

**Facies 3.** The presence of color mottling, in addition to the absence of stratification, in the fine-grained siltstones of facies 3 suggests soil formation, resulting from periods of subaerial exposure (e.g. [29], [30]). The paleosol profiles are particularly rich in carbonates (Table 2). Broken mollusk shells and other skeletal grains indicate that they were transported by relatively high currents before deposition within fine sediments. Moreover, mud coated grains reflect rather shallow water conditions [26], [27]. Charcoal grains are less abundant than in the previous facies. The large size of tracheid fragments suggests better preservational conditions, and/or the proximity of source plants to the depositional site, and/or low hydrodynamic conditions. The presence of both broken skeletal elements and well preserved phytoclasts in fine-grained

<table><caption>**Table 2.** Carbonate concentration of paleosol samples. The table displays the CaCO₃ content of calcareous siltstones. All samples are particularly rich in carbonates. See Fig 2 for correlation with these samples.</caption><thead><tr><th>Sample</th><th>VBN-01</th><th>VBN-02</th><th>VBN-22</th><th>VBN-23</th><th>VBN-25</th><th>VBN-27</th><th>VBN-28</th><th>VBN-31</th><th>VBN-32</th><th>VBN-33</th><th>VBN-41</th></tr></thead><tbody><tr><td>CaCO<sub>3</sub> (%)</td><td>23.8</td><td>24.3</td><td>48.5</td><td>93.4</td><td>38</td><td>30.9</td><td>34.1</td><td>50</td><td>32.7</td><td>29.7</td><td>38.1</td></tr></tbody></table>

doi:10.1371/journal.pone.0134231.t002
sediments indicate an intermediate sub-environment, intermediate between the active channel and surrounding floodplain. Therefore, this facies is tentatively interpreted as a proximal floodplain environment.

**Facies 4.** Facies 4 is the finest grained, consisting predominantly of siltstones and mudstones. It is interpreted as a low-energy alluvial plain facies rather far from the active channel, according to the size of the grains, the oxide-coated grains, the quantity of mud and the lack of bioclasts. Deposits are variegated, probably indicating local changes from oxidizing to reducing conditions [29]. The silty sediments were exposed subaerially and subject to soil-forming processes leading to paleosol formation in a well-drained floodplain as it is shown by mottled sediments, and calcareous nodules (e.g. [31]). The mudstones are horizontally laminated, and are interpreted as a result of settling of suspended sediment in overbank areas during flood events [32]. The low amount of lignitic particles in palynological assemblages, although they are well preserved, suggests long-range transportation from the source.

**Facies 5.** Undulatory laminae at the base of the lowermost limestone bed (Fig 5f) suggest deposition in a shallow lake with bottom currents [19]. Upper beds are massive without sedimentological structure. Because of its very fine texture and the presence of ostracods, this calcareous facies may be interpreted as freshwater limestone deposits within an alluvial setting [33]. Lake formation could have resulted from avulsion when channels were abandoned [18]. During periods of high discharge, flows can spread out onto the floodplain and then fill abandoned channel to form a lake [34]. On the other hand, such a lake could have been formed independently from the fluvial system. Absence of any form of sedimentary structure in the upper carbonate beds probably represents low-energy lake sedimentation with high biogenic productivity [19], [31]. These micritic limestones correspond to autochthonous deposition of carbonate mud. The microfacies reveal very fine matrix and the presence of freshwater “algae”, mollusk shells and ostracod valves [35]. The fine-grained sediments reflect a very calm depositional environment without channel sediment contribution. Moreover, the presence of fossils and the homogeneous texture of the lacustrine wackestone indicate a shallow low-energy depositional context [36], [37], as well as the decrease in the detrital input.

## Taphonomy

The vertebrate assemblage collected in 2009 and 2012 includes two chelonians, *Solemys* and *Polysternon* [38], which are freshwater taxa that lived in fluvial environments [39], [40]. Chelonian specimens are abundant within the different taphonomic modes and are usually well preserved (Fig 9), suggesting a local source for the turtle carcasses. The crocodilian skull and a mandible probably belonging to a basal Alligatoroïdae are similar to a specimen from this crocodilian family [41] previously found in Romanian pond deposits, suggesting that they lived in lacustrine habitats [40]. The excellent preservation of several crocodilian specimens in fluvial deposits from the Velaux-La Bastide Neuve locality also indicates short transport followed by a rapid burial at the local site and suggests that this crocodilian could also live in fluvial environments. The presence of teeth of the hybodont shark *Meristonoides* (Gilles Cuny, pers. com.) is rather unusual. However, hybodont sharks are abundant in freshwater environments during the Cretaceous, although it is difficult to determine precisely if they lived only or occasionally in freshwater environments [42].

Predation marks are scarce: only one supposed predation mark is observed among the total assemblage. Tooth marks on dinosaur bones are usually attributed to scavenging and prey carcass utilization [43], [16], and [21]. If so the scarcity of such marks on the bones of Velaux-La Bastide Neuve assemblage might reflect a low incidence of scavenging at the locality. Alternatively it may also be hypothesized that prey were particularly numerous and that scavengers
only ate the fleshy part of carcasses [44] or preferred some portions of the skeletal remains [45]. Tooth-marked bones are never frequent in dinosaur localities and that theropod dinosaurs did not routinely bite bones during prey carcass utilization [21].

*Post-mortem* fractures are frequent in the fossil assemblage, but bone edges usually remain sharp, indicating that the transport that affected the dead carcasses was rather short [43]. However, numerous parallel transverse fractures (perpendicular to bone fibers) represent *post-mortem* fracturing that affects bones during burial compaction [46], while longitudinal breakages represent weathering, a pre-diagenetic process [16]. Oblique and in rungs breakages are also common and consist of pre-fossilization fracturing, revealing fresh bone damage [44].

Perthotaxic features (bone modification processes active on the land surface [47]) are not found on the bones inspected here. Weathering of bone surfaces is weak in the assemblage. This indicates that the bones were not exposed subaerially for any significant time before or after reworking and that they were quickly buried. The chelonian remains are very well preserved and usually less degraded than dinosaur bones, suggesting that their transport was of shorter duration and that they lived directly in (or very close to) the area.

Most of vertebrate remains are very well preserved in the three fossil-bearing layers. The exceptional preservation of two parallel sauropod cervical ribs, associated with skin impressions in the lowermost layer, reflect rapid burial [48] in fine-grained sediments after short transport in the alluvial plain. Associated elements (partially articulated titanosaur skeletons) from the lowermost bone bed are not complete suggesting that the rest of the skeleton was dissociated earlier.

The taxonomic diversity of the two bone beds suggests distinct accumulations of vertebrate remains belonging to individuals that died in different places and/or at different times. Their remains were then transported by river streams from their different habitats, deposited and finally gathered in the same area in a channel.

It appears that the lowermost bone bed presents a non-random orientation pattern with a major south-east trend for the elongate elements (Fig 11). In this sandy layer elements are heterogeneous with various sizes and shape ranges. This orientation pattern with heterogeneous elements results probably from larger elements acting as obstacles for the other, smaller, elements [49]. No clear preferential orientation in the uppermost layer probably results from low-energy depositional conditions. In this case, currents were probably too weak to orient all the elements in the same direction. The apparently random orientation in this layer could also originated in the abundance of elements or by their sudden deposition. It has been demonstrated that bones transported in shallow water adopt a parallel orientation to the flow, suggesting a NE-SW direction for the currents that transported these bones [49].

## Sedimentology and taphonomy

Sedimentology and taphonomy, together, permit a reconstruction of the paleoenvironments in Velaux-La Bastide Neuve locality. Interpretations resulting from both sedimentological and taphonomic studies, indicate that the bones were deposited in a fluvial environment. The size of the grains, related to downstream decrease in discharge, is a good indicator for assessing the distance between the depositional environments from the main stream channel. Coarser sediments (conglomeratic and sandy facies) deposited close to the channel, were accompanied by heavier vertebrate elements. Finer sediments were deposited further in the alluvial plain, with lighter elements. Sandy fossil-bearing deposits are interpreted as either channel or crevasse splay facies, while silty deposits are interpreted as overbank facies.

Organic particles are typically found in the finer sediments because of their low densities. Thus, organic matter reflects a low energy depositional environment far from the river stream.
Numerous charcoal fragments might indicate the presence of wildfires in the paleoenvironment. Moreover wildfires might intensify landslides and thus influence sedimentation [50]. On the other hand, charcoal might also be the result of plant debris maturation. Charcoal is common in palynological assemblages. Palynofacies therefore confirm the results of the sedimentological studies. Preservation of charcoal, lignitic debris and other organic particles probably reflect flooding events that concentrated phytoclasts. The scarcity of pollen grains in the sediment is more likely the result of the distance from the vegetation sources than their non-preservation.

The vertebrate assemblage discovered at Velaux-La Bastide Neuve clearly represents a mixture of aquatic and terrestrial animals that lived in different habitats. Transported vertebrate carcasses appear to be well preserved both in channel-fills and overbank deposits, due to rather short transport in a low-energy fluvial system. Chelonians and crocodilians were potentially (para) "autochthonous"—their carcasses experienced very little transport inside their life environment-, whereas dinosaurs and sharks were "allochthonous": the latter did not live directly inside the fluvial realm, but in its vicinity and their remains were transported.

During flood events, the river water could overflow channel banks and spread out over the alluvial plain. The finest deposits result from settling of material further from the active channel and accumulated in overbank areas. The river system is interpreted to have had a broad floodplain subjects to periods of soil development. Major vertebrate taxa (dinosaurs, turtles, pterosaurs and crocodiles) lived in this fluvial environment. Palynological analysis confirm that this environment was surrounded by flowering plants and conifers. Repeatedly rainfall could induce flooding events, which resulted in transport of vertebrate remains and organic matter within the system. Different sub-environments then succeeded in the same area and the sedimentological facies in Velaux record these fluvial deposit successions.

## Conclusion

The sedimentological record observed in Velaux-La Bastide Neuve site represents a vertical succession of five facies that are interpreted as channel and floodplain deposits in the medial part of a fluvial depositional system; moreover, the section shows evidence for both river and lacustrine environments. There are nine sedimentological sequences dominated by overbank deposits representing successive flooding episodes. Each sequence is characterized by erosive and depositional episodes. Relatively fine-grained sediments and amalgamation of channel-fills suggest a river system surrounded by a floodplain. Lack of floodplain deposits overlying the sandstones in sequences 2 and 3 suggest several episodes in the river channel down cutting within this zone.

Palynological study corroborates the environmental model proposed, confirming relatively low hydrodynamic depositional conditions and dry climate alternating with wet periods. Vascular plants were present in the paleoenvironment. Palynological assemblages show that angiosperms were also represented in the environment. A few triaperturate pollen grains possibly belonging to Normapolles complex indicate that the vegetation from the Velaux region was probably composed of both gymnosperms (Coniferales) and angiosperms. Transported vertebrate carcasses appear to be well preserved both in channel-fills and overbank deposits, due to rather short transport in a low-energy fluvial system.

## Acknowledgments

This work was realized with the cooperation of the Conseil Général des Bouches-du-Rhône (N. Mouly, M. Bourrelly and G. Michel) and the Velaux city (J.P. Maggi, L. Melih with its heritage, culture and technical services: M. Calvier, C. Cauhapé and S. Chauvet).
The Bouches-du-Rhône Departmental Council (CG 13) provided the official authorizations (Licence: MAPADGAC23112010-1 and MAPADGAC16012014-1-AAPC) to the PALAIOS Association for excavating the Velaux-La Bastide Neuve site (GPS coordinates: 43° 32' 6.58" N, 5° 12' 42.75" E).

We thank the numerous field volunteers and N. Bardet (MNHN Paris), L. Cavin (MHN Geneve), G. Cuny (Univ. Lyon 1), M. Feist (Univ. Montpellier II), X. Pereda Suberbiola (Univ. del Pais Vasco), C. Verati (Univ. Nice Sophia Antipolis), L. Villier (Univ. Paris 6) and F. Fournier (Univ. of Provence), J.F. Babinot, G. Tronchetti, J. Philip (Univ. Aix-Marseille) for the scientific discussions. We are most grateful P. Lauters, S. Berton and T. Hubin for technical support.

## Author Contributions

Conceived and designed the experiments: AC JY MB. Performed the experiments: AC JD MB. Analyzed the data: AC JD. Contributed reagents/materials/analysis tools: XV GG JY. Wrote the paper: AC JY PG GG XV JD. Gave authorizations to make excavations: SA.

## References

1. Allain R, Pereda-Suberbiola X (2003) Dinosaurs of France. Palevol 2: 27-44.

2. Valentin X, Godefroit P, Tabuce R, Vianey-Liaud M, Wenhao W, Garcia G (2012) First Late Maastrichtian (Late Cretaceous) vertebrate assemblage from Provence (Vitrolles-la-Plaine, southern France). In: Godefroit P., editor. Bernissart dinosaurs and Early Cretaceous terrestrial ecosystems. pp. 582-597.

3. Vila B, Galobart A, Canudo JL, Le Loeuff J, Dinarès-Turell J, Riera V, et al. (2012) The diversity of sauropod dinosaurs and their first taxonomic succession from the latest Cretaceous of southwestern Europe: Clues to demise and extinction. Palaeogeography, Palaeoclimatology, Palaeoecology 350-352: 19-38.

4. Tabuce R, Tortosa T, Vianey-Liaud M, Garcia G, Lebrun R, Godefroit P, et al. (2013) New eutherian mammals from the Late Cretaceous of Aix-en-Provence Basin, south-eastern France. Zoological Journal of the Linnean Society 169: 653-672.

5. Thouand E (2004) Biodiversité reptilienne et paléoenvironnement d'un gisement campanien du Bassin d'Aix-en-Provence: La Bastide Neuve (Velaux), M.Sc. Thesis, Université de Poitiers.

6. Garcia G, Amico S, Fournier F, Thouand E, Valentin X (2010) A new titanosaur genus (Dinosauria, Sauropoda) from the Late Cetaceous of southern France and its paleobiogeographic implications. Bulletin de la Société Géologique de France 181(3): 269-277.

7. Westphal M, Durand JP (1989) An Upper Cretaceous paleomagnetic pole for stable Europe from Aix-en-Provence (France). Earth and Planetary Science Letters 94: 143-150.

8. Cojan I, Renard M, Emmanuel L (2003) Palaeoenvironmental reconstruction of dinosaur nesting sites based on a geochemical approach to eggshells and associated palaeosols (Maastrichtian, Provence Basin, France). Palaeogeography, Palaeoclimatology, Palaeoecology 191: 111-138.

9. Plaziat J (1981) Late Cretaceous to Late Eocene Palaeogeographic evolution of southwest Europe. Palaeogeography, Palaeoclimatology, Palaeoecology 36: 263-280.

10. Garcia G, Vianey-Liaud M (2001) Dinosaur eggshells as biochronological markers in Upper Cretaceous continental deposits. Palaeogeography, Palaeoclimatology, Palaeoecology 169: 153-164.

11. Schnyder J, Dejax J, Keppers E, Nguyen-Tu T, Spagna P, Bouliba S, et al. (2009) An Early Cretaceous lacustrine record: organic matter and organic carbon isotopes at Bernissart (Mons Basin, Belgium). Palaeogeography, Palaeoclimatology, Palaeoecology 281: 79-91.

12. Gradstein F, Ogg J, Schmitz M, Ogg G (2012) International Stratigraphic Chart. International Commission on Stratigraphy. In: The Geologic Time Scale. p. 435.

13. Steiger RH, Jäger E (1977) Subcommission on geochronology: convention on the use of decay constants in geo- and cosmochronology. Earth and Planetary Science Letters 36: 359-362.

14. Murru M, Ferrara C, Da Pelo S, Ibba A (2003) The Palaeocene-Middle Eocene deposits of Sardinia (Italy) and their palaeoclimatic significance. Comptes Rendus Geoscience 335: 227-238.

15. Kabanov P, Anadon P, Krumbein E (2008) *Microcodium*: an extensive review and a proposed non-rhizogenic biologically induced origin for its formation. Sedimentary Geology 205: 81-93.
16. Behrensmeyer AK (2007) Bonebeds through Time. In: Bonebeds: genesis, analysis, and paleobiological significance. Rogers RR, Eberth DA, Fiorillo AR, eds. The University of Chicago Press, Chicago.

17. Fiorillo A (1991b) Prey bone utilization by predatory dinosaurs. Palaeogeography, Palaeoclimatology, Palaeoecology 88: 157–166.

18. Blasco R, Rosell J, Peris JF, Caceres I, Vergès JM (2008) A new element of trampling: an experimental application on the Level XII faunal record of Bolomar Cave (Valencia, Spain). Journal of Archaeological Science 35: 1605–1618.

19. Stewart D (1981) A meander-belt sandstone of the Lower Cretaceous of Southern England. Sedimentology 28: 1–20.

20. Smoot J (1991) Sedimentary facies and depositional environments of early Mesozoic Newark Supergroup basins, eastern North America. Palaeogeography, Palaeoclimatology, Palaeoecology 84: 369–423.

21. Viseras C, Soria JM, Duràn JJ, Pla S, Garrido G, Garcia-Garcia F, et al. (2006) A large-mammal site in a meandering fluvial context (Fonelas P-1, Late Pliocene, Guadix Basin, Spain). Sedimentological keys for its paleoenvironmental reconstruction. Palaeogeography, Palaeoclimatology, Palaeoecology 242: 139–168.

22. Fiorillo A (1991a) Taphonomy and depositional setting of Careless Creek Quarry (Judith River Formation), Wheatland County, Montana, U.S.A. Palaeogeography, Palaeoclimatology, Palaeoecology 81: 281–311.

23. Guyonnet-Benaize C, Lamarche J, Masse JP, Villeneuve M, Viseur S (2010) 3D structural modeling of small-deformations in poly-phase faults pattern. Application to the Mid-Cretaceous Durance uplift, Provençance (SE France). Journal of Geodynamics 50: 81–93.

24. Scott A, Cripps JA, Collinson ME, Nichols GJ (2000) The taphonomy of charcoal following a recent heathland fire and some implications for the interpretation of fossil charcoal deposits. Palaeogeography, Palaeoclimatology, Palaeoecology 164: 1–31.

25. do Nascimento DR, Sawakuchi AO, Guedes CCF, Giannini PCF, Grohmann CH, Ferreira MP (2015) Provenance of sands from the confluence of the Amazon and Madeira rivers based on detrital heavy minerals and luminescence of quartz and feldspar. Sedimentary Geology 316: 1–12.

26. Srivastava P (2006) Meso-Neoproterozoic coated grains and palaeoecology of associated microfossils: The Deoban Limestone, Lesser Himalaya, India. Palaeogeography, Palaeoclimatology, Palaeoecology 239: 241–252.

27. Védrine S, Strasser A, Hug W (2007) Oncoid growth and distribution controlled by sea-level fluctuations and climate (Late Oxfordian, Swiss Jura Mountains). Facies 53: 535–552.

28. Brierley GJ, Ferguson RJ, Woolfe KJ (1997) What is a fluvial levee? Sedimentary Geology 114: 1–9.

29. Wright VP (1992) Paleosol recognition: a guide to early diagenesis in terrestrial settings. Developments in Sedimentology 47: 591–619.

30. Pimentel NL, Wright VP, Azevedo TM (1996) Distinguishing early groundwater alteration effects from pedogenesis in ancient alluvial basins: examples from the Paleogene of southern Portugal. Sedimentary Geology 105: 1–10.

31. Csiki Z, Grigirescu D, Codrea V, Therrien F (2010) Taphonomic modes in the Maastrichtian continental deposits of the Hateg Basin, Romania-Palaeoecological and palaeobiological inferences. Palaeogeography, Palaeoclimatology, Palaeoecology 293: 375–390.

32. Gosh P, Sarkar S, Maulik P (2006) Sedimentology of a muddy alluvial deposit: Triassic Denwa Formation, India. Sedimentary Geology 191: 3–36.

33. Armenteros I, Daley B, Garcia E (1997) Lacustrine and palustrine facies in the Bembridge Limestone (late Eocene, Hampshire Basin) of the Isle of Wight, southern England. Palaeogeography, Palaeoclimatology, Palaeoecology 128: 111–132.

34. Nichols G, Fisher J (2007) Processes, facies and architecture of fluvial distributary system deposits. Sedimentary Geology 195: 75–90.

35. Turkmen I, Aksoy E, Tasgin C (2007) Alluvial and lacustrine facies in an extensional basin: The Mio-cene of Malatya basin, eastern Turkey. Journal of Asian Earth Sciences 30: 181–198.

36. Tofalo O, Pazos P (2010) Paleoclimatic implications (Late Cretaceous- Paleogene) from micromorphology of calcretes, palustrine limestones and silcretes, southern Paraná Basin, Uruguay. Journal of South American Earth Sciences 29: 665–675.

37. Wright V, Alonso Zarza A, Sanz M, Calvo J (1997) Diagenesis of Late Miocene micritic lacustrine carbonates, Madrid Basin, Spain. Sedimentary Geology 114: 81–95.

38. Lapparent de Broin F, Murelaga X (1996) Une nouvelle faune de chéloniens dans le Crétacé supérieur européen. Comptes Rendus de l'Académie des Sciences Paris 323: 729–735.
39. Marmi J, Vila B, Galobart A (2009) *Solemys* (Chelonii, Solemydidae) remains from the Maastrichtian of pyrenees: evidence for a semi-aquatic lifestyle. Cretaceous Research 30: 1307–1312.

40. Marmi J, Lujàn AH, Riera V, Gaete R, Oms O, Galobart A (2012) The youngest species of *Polysternon*: a new bothremydid turtle from the uppermost Maastrichtian of the Southern Pyrenees. Cretaceous Research 35: 133–142.

41. Delfino M, Codrea V, Folie A, Dica P, Godefroit P, Smith T (2008) A complete skull of *Allodaposuchus precedens* Nopsca, 1928 (Eusuchia) and a reassessment of the morphology of the taxon based on Romanian remains. Journal of Vertebrate Paleontology 28(1): 111–122.

42. Cuny G (2012) Freshwater hybodont shark in Early Cretaceous ecosystems: a review. In: Bernissart Dinosaurs and early Cretaceous terrestrial ecosystems. Bloomington: Indiana University Press edit. p. 629.

43. Fiorillo A (1988) Taphonomy of Hazard Homestead Quarry (Ogallala Group), Hitchcock County, Nebraska. University of Wyoming Contributions to Geology 26: 57–98.

44. Eberth D, Getty M (2005) Ceratopsians bonebeds: Occurence, origins, and significance. In: Dinosaur Provincial Park, a Spectacular Ancient Ecosystem Revealed. Indiana University Press. pp. 501–506.

45. Faith J, Behrensmeyer A (2006) Changing patterns of carnivore modification in a landscape bone assemblage, Amboseli Park, Kenya. Journal of Archaeological Science 33: 1718–1733.

46. Kos A (2003) Characterization of post-depositional taphonomic processes in the accumulation of mammals in a pitfall cave deposit from southeastern Australia. Journal of Archaeological Science 30: 781–796.

47. Clark J, Beerbower J, Kietzke K (1967) Oligocene sedimentation, stratigraphy, paleoecology, and paleoclimatology in the Big Badlands of South Dakota. Fieldiana Geology 5: 1–158.

48. Eberth D, Currie P (2005) Vertebrate taphonomy and taphonomic modes. In: Currie, editor. Dinosaur provincial park, a spectacular ancient ecosystem revealed. Bloomington: Indiana University press, pp. 453–477.

49. Domínguez-Rodrigo M, Uribelarrea D, Santonja M, Bunn HT, García-Pérez A, Pérez-González A, et al. (2014) Autochthonous anisotropy of archaeological materials by the action of water: experimental and archaeological reassessment of the orientation patterns at the Olduvai sites. Journal of Archaeological Science 41: 44–68.

50. Brown S, Scott A, Glasspool I, Collinson M (2012) Cretaceous wildfires and their impact on the Earth System. Cretaceous research 36: 162–190.