1. Introduction
There are no well-constrained radiometric dates, closely tied to ammonite biostratigraphy, currently available for the whole of the Upper Jurassic (Gradstein et al. Reference Gradstein, Ogg, Schmitz and Ogg2012). Some Upper Jurassic Ar–Ar dates are integrated as secondary guides into the GTS2012: (1) a suite of dates from the almost totally non-marine Morrison Formation in the USA (Gradstein, Ogg & Smith, Reference Gradstein, Ogg and Smith2004; Ogg, Ogg & Gradstein, Reference Ogg, Ogg and Gradstein2008); (2) dates from Oxfordian tuffs intercalated with terrestrial sediments in China (Chang et al. Reference Chang, Zhang, Renne and Fang2009); and (3) dates from ocean-floor basalt veins in the Pacific (Gradstein et al. Reference Gradstein, Ogg, Schmitz and Ogg2012). A single Re–Os date is available from ammonite-bearing marine sedimentary successions in the Lower Kimmeridgian (Selby, Reference Selby2007). As a consequence, the Late Jurassic Time Scale derives mainly from the Pacific seafloor-spreading numerical model of the M-sequence magnetic polarity pattern and from limited recent cyclostratigraphic studies (Ogg & Smith, Reference Ogg, Smith, Gradstein, Ogg and Smith2004; Ogg et al. Reference Ogg, Coe, Przybylski and Wright2010; Gradstein et al. Reference Gradstein, Ogg, Schmitz and Ogg2012). Magnetostratigraphy can be calibrated with ammonite assemblage biochronology, which is mainly defined in northwestern European domains (Cariou & Hantzpergue, Reference Cariou and Hantzpergue1997; Morton, Reference Morton2006). However, provincialism in Boreal, sub-Boreal, sub-Mediterranean and Tethyan domains prevents unequivocal zonation correlation, especially for certain intervals, and hence introduces a temporal bias in the magnetostratigraphic model. Despite recent progress in reducing this bias (Ogg et al. Reference Ogg, Coe, Przybylski and Wright2010; Przybylski et al. Reference Przybylski, Ogg, Wierzbowski, Coe, Hounslow, Wright, Atrops and Settles2010; Gradstein et al. Reference Gradstein, Ogg, Schmitz and Ogg2012), the scarcity of interbedded volcanic units in ammonite-bearing marine successions hinders the accurate numerical calibration of the Late Jurassic Time Scale, even with the progress made in the GTS2012, including improved numerical ages for stage boundaries, obtained by selecting only single-zircon U–Pb ages, recalculating 40Ar–39Ar dates and more precise magnetostratigraphy and cyclostratigraphy.
Therefore, to obtain radiometrically calibrated tie-points for the Late Jurassic, biostratigraphically constrained volcanic ash layers in Tethyan basins have been studied (Pellenard et al. Reference Pellenard, Deconinck, Huff, Thierry, Marchand, Fortwengler and Trouiller2003; Pellenard & Deconinck, Reference Pellenard and Deconinck2006). Here, we focus on eight volcanic ash layers, weathered into bentonites, sampled in pelagic cherty limestones from the Altopiano di Asiago (Trento Plateau domain, northeastern Italy; Bernoulli & Peters, Reference Bernoulli and Peters1970; Martire, Reference Martire1996). We present a new 40Ar–39Ar radiometric date from one of these bentonites, providing the first radiometric tie-point from biostratigraphically well-constrained sedimentary strata for the Middle Oxfordian and discuss these volcanic events and their potential sources.
2. Material and method
Six bentonite layers were identified by their field characteristics, mineralogy and geochemical features at the Serrada section and a further five, 28 km away, at the Echar and Kaberlaba sections, in the Altopiano di Asiago (Trentino Alto Adige and Veneto regions, Italy; Fig. 1a). Weathering of volcanic ashes into clays produced bentonite deposits during the early stages of diagenesis at the sediment/seawater interface. In the Rosso Ammonitico Veronese (RAV), bentonites appear as continuous, centimetre-thick red or white plastic clay-rich horizons, interbedded with limestones and cherts (Figs 1b, 3a). The RAV is an Upper Bajocian to Tithonian pelagic limestone succession, which can be divided into three units (Figs 1b, 2a; Sarti, Reference Sarti1985; Martire, Reference Martire1992; Martire et al. Reference Martire, Clari, Lozar and Pavia2006). The lower unit (Rosso Ammonitico Inferiore: RAI) and the upper unit (Rosso Ammonitico Superiore: RAS) are composed of massive nodular limestones, while the Rosso Ammonitico Middle unit (RAM), containing all the bentonite layers, consists of thin, evenly bedded, non-nodular, chert-rich limestones. The RAM unit reaches a maximum thickness of 10 m, although it occasionally thins out and disappears (Martire, Reference Martire1996; Fig. 2a).
(a) Locality map of the sections sampled in the Altopiano di Asiago region: Kaberlaba (45°50′27.38″N, 11°29′47.56″E), Echar (45°51′21.24″N, 11°34′44.47″E) and Serrada (45°53′16.10″N, 11°09′12.40″E); (b) view of the three members of the Rosso Ammonitico Veronese (RAV) in the Kaberlaba quarry. Scale bar corresponds to 5.2 m, the thickness of the middle unit (RAM). The RAM unit contains typical chert layers and several bentonite layers. The position of the thickest bentonite AB4 is indicated (white arrow).
(a) Detailed logs of the three sections measured, showing correlations between bentonite layers and the 40Ar–39Ar dated AB4 bentonite (Kaberlaba), attributed to the G. transversarium Biozone. (b) PCA based on Al2O3-normalized Hf, Ga, Th, Ta, La, Zr and Ti concentrations. Circles correspond to proposed correlations between bentonite layers. E1 to E8 number the volcanic events.
(a) Photograph of AB4 and AB5 bentonites intercalated between nodular limestones of the RAM unit (Middle Oxfordian) and overlain by the unconformity between the RAM and RAS units. (b) Apparent age spectra for AB4, showing a well-defined plateau age of 156.1 ± 0.89 Ma (2σ external). (c) MORB-normalized multi-element diagram for bentonite layers, pyroclastic deposits from the Vardar domain and comparative patterns for standard rocks from various geodynamic contexts.
Mineralogical (X-ray diffraction, Biogéosciences, Dijon, France) and elemental analyses (inductively coupled plasma-optical emission spectrometry (ICP-OES) and inductively coupled plasma-mass spectrometry (ICP-MS), CRPG Nancy, France) were performed on all powdered samples to confirm their volcanic nature (online Supplementary Material Table S1 available at http://journals.cambridge.org/geo). Principal Component Analysis (PCA) was used to evaluate the number of volcanic events. Prior to the correlation matrix-based PCA, trace element concentrations were re-expressed assuming an initial volcanic concentration of 15% Al2O3 (Spears et al. Reference Spears, Kanaris-Sotiriou, Riley and Krause1999; Pellenard et al. Reference Pellenard, Deconinck, Huff, Thierry, Marchand, Fortwengler and Trouiller2003). This procedure reduces variability in lithophile element concentration, which could be owing to post-depositional diagenetic processes, such as dilution by authigenic phases, or concentration by dissolution of less stable minerals.
The 40Ar–39Ar dating (OSIRIS reactor CEA Saclay, France) was performed by step-heating about 30 small (< 100 μm) transparent sanidines, carefully hand-picked under a binocular microscope after several treatments from the Kaberlaba section AB4 bentonite (original sample weight 2 kg, see online Supplementary Material available at http://journals.cambridge.org/geo for details). Each Ar isotope measurement consists of 20 cycles by peak switching between the different Ar isotopes. The J value was determined using three single ACs (Alder Creek sanidine) grains taken from the same hole as the sample. Recently, Renne et al. (Reference Renne, Mundil, Balco, Min and Mudwig2010, Reference Renne, Mundil, Balco, Min and Mudwig2011) published an optimization model for estimating the partial decay constants of 40K and 40Ar*/40K ratio of FCs (Fish Canyon sanidine). This calibration reduces systematic uncertainties in the 40Ar–39Ar system from c. 2.5% (Steiger & Jäger, Reference Steiger and Jäger1977) to 0.27%. The optimization model yields an age for ACs of 1.2056 Ma, equivalent to FCs of 28.294 Ma, that overlaps at the 2σ confidence level the astronomically tuned ACs and FCs ages reported by Kuiper et al. (Reference Kuiper, Deino, Hilgen, Krijgsman, Renne and Wijbrans2008). The optimization model of Renne et al. (Reference Renne, Mundil, Balco, Min and Mudwig2010, Reference Renne, Mundil, Balco, Min and Mudwig2011) used pairs of 238U–206Pb and 40Ar–39Ar data as inputs. Therefore, 40Ar–39Ar ages calibrated with this optimization model could be directly compared to U–Pb. The corresponding J value (0.0006846 ± 0.00000137, 2σ) was calculated using the Renne et al. (Reference Renne, Mundil, Balco, Min and Mudwig2011) calibration of ACs. The J uncertainty corresponds to the standard deviation of the weighted mean of three ACs single grains (see Nomade et al. Reference Nomade, Gauthier, Guillou and Pastre2010, Reference Nomade, Muttoni, Guillou, Robin and Scardia2011 and online Supplementary Material available at http://journals.cambridge.org/geo for detailed methodology).
3. Biostratigraphy and correlation of ash layers
At Kaberlaba, calcareous nannofossil assemblages indicate a Late Callovian age for the base of the RAM unit, while the ammonite assemblage composed of Gregoryceras fouquei, Passendorferia (Enayites) birmensdorfensis, Passendorferia cf. ziegleri, Perisphinctes (Otosphinctes) nectobrigensis, Perisphinctes (Dichotomosphinctes) aff. elisabethae, Sequeirosia (Gemmellarites) aff. trichoplocus and Subdiscosphinctes richei, which is characteristic of the Gregoryceras transversarium Biozone, indicates a Middle Oxfordian age for the top of the RAM unit (Clari, Martire & Pavia, Reference Clari, Martire, Pavia, Pallini, Cecca, Cresta and Santantonio1990; Martire, Reference Martire1992, Reference Martire1996; Martire et al. Reference Martire, Clari, Lozar and Pavia2006; see Fig. S1 in online Supplementary Material available at http://journals.cambridge.org/geo for photographs of typical ammonites of the G. transversarium Biozone). All these ammonite taxa come from the bed between bentonites AB3 and AB4 at Kaberlaba, where preservation is better than in the rest of the section. They are all exclusive to the G. transversarium Biozone, except for G. fouquei, which spans both the G. transversarium Biozone and the Perisphinctes (Dichotomoceras) bifurcatus Biozone. The overlying RAS unit contains ammonites such as Orthosphinctes (Ardescia) gr. inconditus, Crussoliceras aceroides and Idoceras (Lessiniceras) sp., which are characteristic of the Taramelliceras strombecki and Presimoceras herbichi biozones of the Lower Kimmeridgian (Sarti, Reference Sarti1993; Clari, Martire & Pavia, Reference Clari, Martire, Pavia, Pallini, Cecca, Cresta and Santantonio1990; Martire, Reference Martire1992, Reference Martire1996). Therefore, at Kaberlaba, there is a major hiatus (four ammonite biozones) between the upper part of the Middle Oxfordian and the lowermost part of the Lower Kimmeridgian (Fig. 2a). However, the RAM unit of the Echar section provides a good biostratigraphic framework for bentonites AB4 and AB5, as here the overlying sediments are well dated, with no hiatus. The RAM unit at Echar contains the same five bentonites and is overlain by three stromatolitic beds, the first of which belongs to the G. transversarium Biozone, with the same taxa as Kaberlaba. The second stromatolitic bed is dated to the Lower Kimmeridgian (Sowerbyceras silenum Biozone), on the basis of the following assemblage: Taramelliceras cf. rigidum, Idoceras (Lessiniceras) cf. raschii, Lithacosphinctes cf. stromeri, Mesosimoceras evolutum and Euaspidoceras (Epaspidoceras) sp. The third stromatolitic bed belongs to the P. herbichi Biozone (Lower Kimmeridgian, Fig. 2a). In the Serrada section, the RAM unit extends from the Upper Callovian to the Middle Oxfordian (G. transversarium Biozone). As all bentonites sampled were from the RAM unit, they therefore date from the Upper Callovian to the Middle Oxfordian. As ammonites diagnostic of the G. transversarium Biozone were found just below AB4 at Kaberlaba and just above AB5 at Echar, the two uppermost bentonite beds in these sections (AB4 and AB5), easily recognizable because of their thickness and vivid red colour (Fig. 3a), are attributed to the G. transversarium Biozone (Fig. 2a).
The bentonites studied, which correspond to pure smectite horizons containing occasional volcanic crystals (e.g. sanidine, quartz, biotite), are marked by positive anomalies in Th, Ta, Hf and Ga, characteristic of bentonite deposits (Spears et al. Reference Spears, Kanaris-Sotiriou, Riley and Krause1999; Pellenard et al. Reference Pellenard, Deconinck, Huff, Thierry, Marchand, Fortwengler and Trouiller2003; Table S1 in online Supplementary Material available at http://journals.cambridge.org/geo). PCA was used to examine possible similarities between ash layers in the Serrada and Kaberlaba sections, 28 km apart, in order to correlate the bentonites and to evaluate the number of volcanic events and their preservation in the Trento Plateau domain. The most typically immobile, volcanogenic elements were selected for this analysis: Hf, Ga, Th, Ta, La, Zr and Ti (Fig. 2b). In the F2 v. F1 diagram (Fig. 2b), representing more than 80% of the total variance, four groups consistent with the stratigraphy can be clearly identified: (i) AB2, AB3, (ii) SB1, SB2, AB1, (iii) SB3, AB4, and (iv) SB4, SB5, SB6, AB5. The first Kaberlaba level, AB1, corresponds either to the first Serrada level SB1 or possibly to SB2. Samples AB2 and AB3 (Kaberlaba) have no equivalents in the Serrada section, indicating that these events were not systematically preserved. Sample AB4, a thick red bentonite from Kaberlaba, is geochemically similar to SB3, the thickest bentonite from Serrada. Sample AB5 from Kaberlaba probably corresponds to SB4, perhaps to SB5 or SB6. At least eight individual volcanic events are therefore identified using PCA (Fig. 2b), with correlations over a large geographic area, coherent with the biostratigraphic framework.
4. 40Ar–39Ar results
We used the laser-fusion step-heating 40Ar–39Ar method to date level AB4, which contains the highest abundance of well-preserved sanidines and which is also biostratigraphically the most precisely constrained. The apparent age spectrum obtained for the AB4 sanidines is 100% concordant (Fig 3b, details in online Supplementary Material Tables S2 and S3 available at http://journals.cambridge.org/geo): all steps yield indistinguishable ages, with a well-defined plateau age of 156.1 ± 0.89 Ma (2σ full uncertainty propagation). As the inverse isochron displays low scatter because of its highly radiogenic content, it was not used, given the imprecise initial atmospheric40Ar/36Ar ratio obtained. The plateau age we obtain can be directly compared to U–Pb ages available for the Jurassic Time Scale (GTS2004, GTS2012 and Pálfy, Reference Pálfy2008). We, however, present alternative calculations (Table 1) using several standards (Nomade et al. Reference Nomade, Renne, Vogel, Deino, Sharp, Becker, Jaouni and Mundil2005; Kuiper et al. Reference Kuiper, Deino, Hilgen, Krijgsman, Renne and Wijbrans2008) and total 40K decay constants (Steiger & Jäger, Reference Steiger and Jäger1977; Renne et al. Reference Renne, Mundil, Balco, Min and Mudwig2011) since the Mesozoic GSSP time scale (GTS2004 and GTS2012) was based on many 40Ar–39Ar ages, using different 40K constant and various standards. The full uncertainty propagation of the Steiger & Jäger (Reference Steiger and Jäger1977) 40K total decay constant (c. 2.5% at 2σ) results in an AB4 error of about 4.0 Ma, while the Min et al. (Reference Min, Mundil, Renne and Ludwig2000) 40K decay constant, proposed by Kuiper et al. (Reference Kuiper, Deino, Hilgen, Krijgsman, Renne and Wijbrans2008), could not be retained because of its high degree of uncertainty of 3.9% at 2σ, compared to the 0.57% from Renne et al. (Reference Renne, Mundil, Balco, Min and Mudwig2010, Reference Renne, Mundil, Balco, Min and Mudwig2011) that has been adopted in this study.
Calculated ages and corresponding uncertainties using various K total decay constants
*The uncertainty reported is the full propagated uncertainty. (1) – Nomade et al. (Reference Nomade, Renne, Vogel, Deino, Sharp, Becker, Jaouni and Mundil2003); (2) – Kuiper et al. (Reference Kuiper, Deino, Hilgen, Krijgsman, Renne and Wijbrans2008); (3) – Renne et al. (Reference Renne, Mundil, Balco, Min and Mudwig2011). ACs – Alder Creek sanidine; FCs – Fish Creek sanidine.
5. Nature and source of volcanic events
The bentonite profile in the MORB-normalized multi-element plot clearly shows that the initial ash layers result from an evolved calc-alkaline magma (Fig. 3c and online Supplementary Material Fig. S2 available at http://journals.cambridge.org/geo). The characteristic Nb depletion and the Hf–Th–Ta diagram are typical of subduction-related arc materials, while the Zr/TiO2 v. Nb/Y diagram indicates mainly andesite to rhyodacite products (online Supplementary Material Fig. S2 available at http://journals.cambridge.org/geo). As no lavas or thick pyroclastic deposits have been identified in or nearby the Trento domain within the Upper Jurassic (Bernoulli & Peters, Reference Bernoulli and Peters1970; Pellenard et al. Reference Pellenard, Deconinck, Huff, Thierry, Marchand, Fortwengler and Trouiller2003), sources must be distant. In addition, fine-grained ashes emitted by highly explosive eruptions are known to be distributed over long distances (> 1000 km). This hypothesis is supported by (i) the correlation indicated by the PCA of several events with similar features (e.g. thickness), over a large area in the Venetian Pre-Alps, and (ii) the size (50–100 μm) of the preserved pyroclastic minerals (i.e. sanidine and quartz). Emissions of tholeiitic basalts, andesites and pyroclastites are reported for the Middle–Late Jurassic from the island-arc magmatism in the eastern Rhodope–Thrace region in Bulgaria and Greece (Bonev & Stampfli, Reference Bonev and Stampfli2008). This volcanism was associated with the southward subduction of the Meliata–Maliac Ocean under the supra-subduction back-arc Vardar Ocean/island-arc system (Bonev & Stampfli, Reference Bonev and Stampfli2008). The Vardar geodynamic context undoubtedly produced huge eruptions and subsequent widespread ashes. The age of the Vardar subduction, ranging from the Early Jurassic incipient proto-arc to the Middle–Late Jurassic arc–back-arc spreading, is coherent with the biostratigraphic age of the bentonites studied here, whose geochemical fingerprint is similar to that of the Vardar pyroclastics (Fig. 3c). This evidence supports Vardar island-arc volcanism as the probable source of the ash layers found in the Venetian Pre-Alps.
6. A new tie-point for the Late Jurassic Time Scale
There are few biostratigraphically well-constrained radiometric tie-points for the Middle–Late Jurassic. For the Middle Jurassic, the only available U–Pb ages are from (i) British Columbia bentonites, ascribed to the early Late Bathonian (Pálfy, Reference Pálfy2008), and (ii) an ash layer (164.6 ± 0.2 Ma) in the Neuquén province (Argentina), at the Bathonian–Callovian boundary (Kamo & Riccardi, Reference Kamo and Riccardi2009). There are no biostratigraphically well-constrained radiometric ages for the Oxfordian–Tithonian interval, while only a few 40Ar–39Ar dates from oceanic basalts are retained in the current GTS2012: (i) 159.86 ± 3.33 (2σ) Ma and 161.17 ± 0.74 (2σ) Ma from Pacific tholeiitic basalts (site 801) assigned to the Oxfordian, based on radiolarian calibration, (ii) a revised 156.3 ± 3.4 (2σ) Ma reported for the Hawaiian basalt seafloor (site 765), correlated to the base of the Kimmeridgian (P. baylei ammonite zone) using the M26r magnetochron (Gradstein et al. Reference Gradstein, Ogg, Schmitz and Ogg2012, appendix 2, p. 1045), and (iii) an earliest Berriasian 40Ar–39Ar date of 145.5 ± 0.8 (2σ) Ma from an oceanic basalt sill in the Pacific Ocean (Mahoney et al. Reference Mahoney, Duncan, Tejada, Sager and Bralower2005). Robust 40Ar–39Ar ages of 160.7 ± 0.4 (2σ) Ma and 158.7 ± 0.6 (2σ) Ma have recently been obtained from two tuffs in the Lanqi Formation in northeastern China, but the terrestrial fossils do not allow the attribution of a more precise stratigraphy than a Late Jurassic age (Chang et al. Reference Chang, Zhang, Renne and Fang2009). The only biostratigraphically well-constrained age, documented by Selby (Reference Selby2007) on an organic-rich mudrock deposit from the Isle of Skye (Scotland), yields a Re–Os age of 154.1 ± 2.2 Ma (2σ) in the Lower Kimmeridgian, just above the proposed Oxfordian/Kimmeridgian GSSP.
As a consequence, Middle–Upper Jurassic biozone duration and stage boundary ages are mainly estimated from secondary radiometric guides, indirect methods and mathematical interpolations. These approaches combine a magnetostratigraphic age model based on the cycle-scaling of the M-sequence spreading rate model correlated to the magnetostratigraphy of outcrops (Ogg et al. Reference Ogg, Coe, Przybylski and Wright2010; Przybylski et al. Reference Przybylski, Ogg, Wierzbowski, Coe, Hounslow, Wright, Atrops and Settles2010; Gradstein et al. Reference Gradstein, Ogg, Schmitz and Ogg2012) and cycle-derived durations of ammonite zones from cyclostratigraphy (Boulila et al. Reference Boulila, Hinnov, Huret, Collin, Galbrun, Fortwengler, Marchand and Thierry2008, Reference Boulila, de Rafélis, Hinnov, Gardin, Galbrun and Collin2010; Ogg, Ogg & Gradstein, Reference Ogg, Ogg and Gradstein2008; Huang, Hesselbo & Hinnov, Reference Huang, Hesselbo and Hinnov2010; Gradstein et al. Reference Gradstein, Ogg, Schmitz and Ogg2012). Cyclostratigraphy from SE France has considerably modified ammonite biozone durations. Using a condensed section in Britain, the entire Oxfordian stage had previously been fixed at 0.6 Ma, in the GTS2004. New data from cyclostratigraphy suggest that the Oxfordian spanned 6.0 Ma, with 2 Ma attributed to the Quenstedtoceras mariae Zone alone (Boulila et al. Reference Boulila, Hinnov, Huret, Collin, Galbrun, Fortwengler, Marchand and Thierry2008; Gradstein et al. Reference Gradstein, Ogg, Schmitz and Ogg2012). The age of the Oxfordian/Kimmeridgian boundary is now set at 157.3 ± 1.0 Ma in the GTS2012, whereas it was 155.6 ± 4.0 Ma in the GTS2004 and GTS2008 (Gradstein, Ogg & Smith, Reference Gradstein, Ogg and Smith2004; Ogg, Ogg & Gradstein, Reference Ogg, Ogg and Gradstein2008). In this study, the 40Ar–39Ar age of 156.1 ± 0.89 Ma (2σ full uncertainty propagation), attributed to the G. transversarium Biozone (Middle Oxfordian), is consistent with the existing Re–Os age and the 40Ar–39Ar ages retained as secondary guides in the GTS2012. Nevertheless, it falls outside of the current base and top limits of the G. transversarium Biozone proposed, respectively, at 160.09 ± 1.0 Ma (2σ) and 159.44 ± 1.0 Ma (2σ), both interpolated from Oxfordian stage boundaries (Gradstein et al. Reference Gradstein, Ogg, Schmitz and Ogg2012). The age proposed here remains compatible with the Oxfordian boundaries (163.5 ± 1.1 Ma and 157.3 ± 1.0 Ma) proposed by the GTS2012 if maximum uncertainties are taken into account. However, there is a better fit with the previous Oxfordian base (161.2 ± 4.0 Ma) and top (155.6 ± 4.0 Ma) from the GTS2004 and GTS2008, where the proposed boundaries were around 2 Ma younger.
The age proposed here, well constrained within the standard Jurassic biostratigraphic zonation (Cariou & Hantzpergue, Reference Cariou and Hantzpergue1997), provides the first accurate and reliable numerical age currently available for the Late Jurassic Time Scale. This precise new tie-point can be used to anchor floating cyclostratigraphy and magnetostratigraphy, thus contributing to the improvement of seafloor-spreading models and, above all, will aid in the calibration of the Late Jurassic Time Scale.
Acknowledgements
This work was supported by the Centre National de la Recherche Scientifique and the Commissariat à l'Energie Atomique et aux énergies alternatives. The authors thank Carmela Chateau-Smith for English proof-reading. We are grateful to J. F. Deconinck and D. Bernoulli for support and discussions about bentonites and G. Pavia and G. Meléndez for help with ammonite biostratigraphy. We particularly thank the reviewers A. L. Coe and F. Jourdan whose constructive criticisms and detailed remarks have greatly contributed to improve the manuscript.
