First steps in reconstructing Early Jurassic sea water temperatures in the Andean Basin of northern Chile based on stable isotope analyses of oyster and brachiopod shells

Matthias Alberti,Franz T.Fürsich and Nils Andersen

Abstract The stable isotope (δ13C, δ18O) composition of a collection of Lower Jurassic brachiopods and oysters from the Andean Basin of northern Chile was analyzed. The results allow the first reconstruction of absolute water temperatures for several ammonite zones in the Lower Jurassic of South America. The temperature record starts with comparatively high values in the Late Sinemurian (average: 27.0°C; Raricostatum Zone). Just before the Sinemurian-Pliensbachian transition, temperatures dropped to an average of 24.3°C. The lowest temperature value in the dataset was recorded for a brachiopod shell of the latest Pliensbachian Spinatum Zone (19.6°C). No data are available for the Early Toarcian, but results for the late Toarcian show again comparatively warm conditions(average: 24.4°C; Thouarsense–Levesquei zones). Even though more material and analyses are necessary to corroborate the recorded temperatures, the present dataset seems to indicate the global nature of the Late Pliensbachian Cooling Event. In contrast, the global warming during the Toarcian Oceanic Anoxic Event has not been recorded due to a lack of Early Toarcian material. The δ13C record of brachiopods and oysters documents a gradual increase in values representing background conditions. Oyster shells were used for high-resolution stable isotope analyses and show seasonal temperature fluctuations over a period of around 3 years in the life time of the bivalves. If explained only by temperatures, the δ18O values point to a minimum estimate for the seasonality in the late Toarcian of slightly more than 3°C.

Keywords: Early Jurassic, Chile,Stable isotopes,Palaeoclimate, Water temperatures, Seasonality

1 Introduction

The climate conditions of the Jurassic world have traditionally been described as warmer than today,relatively stable through time, and with weak latitudinal temperature gradients leading to a lack of polar glaciations (e.g., Chandler et al. 1992; Valdes and Sellwood 1992; Rees et al. 2000; Hallam 2001; Bailey et al. 2003; Sellwood and Valdes 2006). In the last few decades, this simple view has been challenged by a series of articles, suggesting polar ice caps for at least short intervals in the Jurassic (e.g., Price 1999;Dromart et al. 2003; Suan et al. 2008, 2010) and larger temperature fluctuations (e.g., Dera et al. 2011;Korte et al. 2015; Martinez and Dera 2015; Silva and Duarte 2015; Price et al. 2016). Only recently, Alberti et al. (2017, 2019) presented data pointing to latitudinal temperature gradients in the Middle and Late Jurassic as steep as today. Much of our knowledge on the Jurassic climate is based on seawater temperature reconstructions using stable isotope (δ18O) analyses of fossil hard parts. This methodology has been fundamental in improving our understanding of the climate throughout the history of Earth (e.g., Epstein et al.1951; Urey et al. 1951). During the last decades a large amount of data has been assembled by analyzing mainly calcitic fossils such as belemnite rostra and bivalve or brachiopod shells. However, research on temperature reconstructions has concentrated largely on European localities with data from other regions of the Earth being relatively scarce and scattered.Consequently, it seems beneficial to gather more information from non-European sections in order to fully understand global climate conditions in the Jurassic. The present study is part of a series of articles trying to improve our knowledge of the seawater temperature conditions of Gondwanan localities(Alberti et al. 2012a, 2012b, 2013, 2017, 2019).

The Early Jurassic was an epoch which experienced several events with a global impact. Most of these have been identified via geochemical analyses of sedimentary rocks or their fossil content (e.g., Jenkyns 1988; Hesselbo et al. 2000a; Jenkyns et al. 2002; Korte et al. 2009; Bodin et al. 2016) and could be connected to pulses in the volcanic activity of the Karoo-Ferrar large igneous province influencing the atmospheric composition and thus the global climate. The strongest of these events was the Toarcian Oceanic Anoxic Event (TOAE) which is evident in the geological record by a prominent carbon isotope excursion and has been studied thoroughly over the last decades (e.g., Jenkyns 1988; Mettraux et al. 1989;Hesselbo et al. 2000a, 2007; McArthur et al. 2000;Schmid-Röhl et al. 2002; van de Schootbrugge et al.2005a; Metodiev and Koleva-Rekalova 2008; Suan et al.2008; Rodríguez-Tovar and Reolid 2013; Huang and Hesselbo 2014; Ullmann et al. 2014; Ruhl et al. 2016;Bougeault et al. 2017). Following a phase of very low temperatures during the Late Pliensbachian, increased volcanic activity in the Early Toarcian in combination with other factors (such as the release of methane hydrates from continental shelfs; Hesselbo et al. 2000a)seems to have caused perturbations in the carbon cycle,intense global warming, ocean acidification, and ultimately mass extinctions (e.g., Aberhan and Fürsich 1996;Pálfy and Smith 2000; Aberhan and Baumiller 2003;Cecca and Macchioni 2004; Gómez et al. 2008; Dera et al. 2010; Dera and Donnadieu 2012; Danise et al.2013; Huang and Hesselbo 2014; Krencker et al. 2014).While the global nature of many of the Early Jurassic events (particularly the TOAE) has been illustrated by geochemical analyses of rocks outside Europe (e.g., Al-Suwaidi et al. 2010, 2014, 2016; Mazzini et al. 2010; Caruthers et al. 2011; Gröcke et al. 2011; Suan et al. 2011;Yi et al. 2013; Kemp and Izumi 2014; Ros-Franch et al.2019), reconstructions of absolute water temperatures for the Early Jurassic are almost completely restricted to the northwestern Tethys. In the present study, a collection of Sinemurian to Toarcian oyster and brachiopod shells from the Andean Basin of Chile has been analyzed for their stable isotope (δ13C, δ18O) content in a first step to remedy this lack of data.

2 Geological overview

During the Mesozoic, a narrow but long back-arc basin existed along the South American Pacific margin more or less parallel to the modern-day Andes (e.g., Gröschke et al. 1988). This Andean Basin in present-day Chile and Argentina was divided into subbasins by several ridges and separated from the open ocean by volcanic arcs, but connections allowing the exchange of water masses existed most of the time (e.g., Gröschke et al. 1988;Aberhan 1992; Vicente 2006). A marine transgression reached the basin in the Late Triassic and marine conditions prevailed throughout the Jurassic (Aberhan 1992)with sea-level maxima in the Toarcian and Bajocian(Fantasia et al. 2018). The study area was inundated in the Early Sinemurian (von Hillebrandt 1971, 1973). The thickness of the Lower and Middle Jurassic successions(up to the basal Callovian) is more than 1000 m in the west, but only 30–40 m in the east (Aberhan 1992; von Hillebrandt 2002).The sampled Lower Jurassic strata belong to the Montandón Formation (Pérez 1982).

The fossil material used in the current study comes from four sections measured in detail by Aberhan (1992;sections with sample positions are illustrated in chapter 5). Section 2 was measured near Potrerillos (Fig. 1c),while sections 4, 5, and 6 were measured around El Transito approximately 280 km further south (Fig. 1d;compare Aberhan 1992). Their exact palaeogeographic position differs between the different available reconstructions, but was probably somewhere in the subtropical zone between 20°S to 35°S during the studied time interval (based on the Paleolatitude Calculator; van Hinsbergen et al. 2015). The analyzed specimens come from 22 horizons within the four sections and have a Sinemurian to Toarcian age.

As in other regions, the Jurassic biostratigraphy of South America relies heavily on ammonites, but age assignments and correlations based on these index fossils are still continuously improved and have not been finalized for all localities and/or time intervals (compare von Hillebrandt 1987; Riccardi 2008). The fossils used in the current study have been assigned to ammonite zones based on results of previous studies (von Hillebrandt 1973, 1987, 2002; von Hillebrandt and Schmidt-Effing 1981; Aberhan and von Hillebrandt 1996; Pérez et al.2008; Fantasia et al. 2018).

Sinemurian ammonites of Chile were discussed in detail by von Hillebrandt (2002). The lower part of the Jurassic succession near El Transito contains only few ammonites, but von Hillebrandt (2002) described some specimens, which he could assign to the Raricostatum Zone. Even though the Sinemurian at these localities has a thickness of more than 200 m, it seems to belong entirely to this ammonite zone indicating a late Sinemurian age (von Hillebrandt 2002).

Fig.1 a-d Schematic maps illustrating the position of the studied sections in northern Chile(modified after Aberhan 1992)

Section 2 is a classic locality in the lower Quebrada El Asiento near Potrerillos examined or mentioned by numerous previous researchers (e.g., von Hillebrandt and Schmidt-Effing 1981; Pérez 1982; Aberhan 1992; von Hillebrandt 2006; Pérez et al. 2008; Fantasia et al. 2018).Shells used in the present study were collected from the Pliensbachian part of the succession, which was divided into three ammonite zones by von Hillebrandt and Schmidt-Effing (1981): Fig. 3): the Davoei, Margaritatus,and Spinatum zones. This scheme was also adopted by Fantasia et al. (2018), who additionally studied the nannofossil biostratigraphy of the same section. These studies were the foundation to assign Bed 2-1 to the upper Davoei to lower Margaritatus Zone and Bed 2-2 to the Margaritatus Zone (for sections with bed numbers see chapter 5). Based on the stratigraphic correlation by Aberhan(1992),Beds 6-7,6-8,and 6-9 should be slightly younger than Beds 2-1 and 2-2. Consequently, they can be assigned to the Spinatum Zone.

Toarcian oyster shells and one brachiopod were available from sections 4, 5, and 6 near El Transito. Beds 4-1,5-7, and 6-10 correspond to a conspicuous oolitic ironstone to bioclastic limestone unit almost directly above a thick package of sandstones. This stratigraphic interval has been assigned to the Upper Toarcian Copiapense and Tenuicostatum zones (not to be confused with the European Lower Toarcian Tenuicostatum Zone;compare von Hillebrandt 1973, 1987; Aberhan and von Hillebrandt 1996), which are equivalent to the Thouarsense and lowermost Levesquei zones of Europe (e.g.,Pérez et al. 2008). The overlying marly limestones with Beds 4-2 and 4-3 have been assigned to the Lotharingica and Fluitans zones (compare Aberhan and von Hillebrandt 1996), which are equivalent to the middle and upper Levesquei Zone of Europe (e.g., Pérez et al. 2008).

3 Material and methods

Fifty-nine shells of oysters and brachiopods from the Lower Jurassic strata of northern Chile were used in the current study. The fossil material was originally collected by Martin Aberhan for a palaeoecological study (Aberhan 1992, 1993) and is now stored in the collections of the Bayerische Staatssammlung für Paläontologie und Geologie in Munich, Germany. The diverse benthic fauna has been identified and described in detail by Aberhan (1992) and contains a number of species with an original calcitic shell enabling stable isotope analyses.Taxa used in the present study include oysters of the genera Gryphaea and Actinostreon as well as rhynchonellid brachiopods of the genera Gibbirhynchia, Tetrarhynchia, Quadratirhynchia, Rudirhynchia, and Rhynchonelloidea. One specimen of the spiriferid brachiopod Spiriferina chilensis was analyzed as well. Illustrations of these taxa can be found in taxonomic studies(e.g., Manceñido 1981; Manceñido and Dagys 1992;Aberhan 1994; Baeza-Carratalá 2013).

All specimens were checked for diagenetic alteration before analyzing their stable isotope (δ13C, δ18O) composition. Sections of 5–10 mm thickness were cut through the shells and their surface ground (Fig. 2a, b).All specimens were then studied with a cold cathodoluminescence microscope at the GeoZentrum Nordbayern of the Friedrich-Alexander-Universität Erlangen-Nürnberg, Germany. Non-luminescent areas were selected for sampling, while specimens with a very thin or luminescent shell were excluded from the collection. Few carbonate samples of presumably altered shells, sediment, and sparitic cement filling shells were measured for comparison purposes. Furthermore, three large oyster shells of the genus Gryphaea were selected for high-resolution stable isotope analyses (Fig. 2c–e).Results of these three shells were used for seasonality reconstructions and listed as well as interpreted separately. Carbonate samples were taken with a hand-held dental drill in the case of thick shells or with a computer-controlled micromill at the GeoZentrum Nordbayern. In summary, the isotopic composition of 105 brachiopod and oyster samples was analyzed using a carbonate preparation device (Kiel IV) connected with a ThermoScientific MAT 253 mass spectrometer at the Leibniz Laboratory for Radiometric Dating and Stable Isotope Research at the Christian-Albrechts-Universität zu Kiel, Germany. The carbonate samples were reacted within the preparation device with 100% orthophosphoric acid at 75°C and the evolved CO2gas was analyzed using the mass spectrometer. On daily routine,different laboratory internal carbonate standards and two international carbonate standards (NBS-19; IAEA-603) were analyzed to control the precision of measured δ13C and δ18O values. All values are reported in per mil relative to the Vienna Pee Dee Belemnite(VPDB) scale using NBS-19. Generally, one sample was collected from each fossil shell and analyzed, but in the case of the three oysters used for high-resolution analyses between 16 to 21 samples were collected. Ten brachiopod shells, which were believed to be well preserved, did not yield enough carbonate powder for stable isotope analyses, because their shells were too thin. For the resulting δ18O values of well-preserved shells, palaeotemperatures were calculated using the equation given by Anderson and Arthur (1983) with a δ18O value of -1‰ VSMOW for seawater during shell precipitation (as suggested for an ice-free Jurassic world; Shackleton and Kennett 1975).

4 Preservation of the fossil material

Calcareous shells of marine organisms can store information on sea water temperatures in Earth’s history as shown by a large number of geochemical studies (e.g.,Epstein et al. 1951; Urey et al. 1951). However, the chemical composition of shells can change after the death of the animal, for example, when pore fluids enter the shells along fractures. It is therefore very important in any such study to carefully examine the preservation of the fossil material used. Consequently, cathodoluminescence microscopy has become a standard procedure for stable isotope research on fossil shells (e.g., Wierzbowski 2002, 2004; Wierzbowski and Joachimski 2007;Ullmann and Korte 2015; Arabas 2016). Most shells in the present study revealed a good preservation with largely dark and non-luminescent shells (Fig. 3a). Only some shells showed a red to orange luminescence thought to be caused by Mn2+which entered the calcite after the burial of the shell (e.g., Fürsich et al. 2005;Wierzbowski et al. 2009; also compare Barbin 2013). In the present cases, these alterations were strongest along fractures,which might have served as pathways for fluids entering the shells (Fig. 3b). Only specimens with large non-luminescent areas, believed to be well-preserved,were used for oxygen isotope analyses and temperature reconstructions, while some altered shells were analyzed for comparison purposes. Cathodoluminescence microscopy led to an exclusion of 15 specimens from the total collection of 59 shells.

Fig. 2 Preparation of the fossil material (modified after Alberti et al. 2012a, 2012b, 2017, 2019). a, b Sections of 5–10 mm thickness were cut from each shell and studied with a cathodoluminescence microscope to select sample areas; c–e Three oyster specimens (CH18-006, CH18-017, and CH18-030)were used for high-resolution stable isotope analyses

Another possible indicator for alteration is a correlation between the δ18O and δ13C values of the analyzed shells. Sedimentary rocks are commonly characterized by lower δ18O and δ13C values than fossil shells, especially after diagenetic alteration or influence by meteoric waters (e.g., Hodgson 1966; Hudson 1977; Nelson and Smith 1996). Similarly, stable isotope ratios of fossil shells shift towards lower values with increasing alteration leading to a correlation between both values in the dataset. Such a correlation should therefore be treated with caution,although in some cases it can be primary in origin,for example, because of an increase in primary productivity with higher water temperatures.

Fig. 3 Cathodoluminescence microscopy. a Well-preserved shell of Rudirhynchia aff. rudis from the Upper Pliensbachian (CH18-001); b Poorly preserved shell of Gryphaea cf.dumortieri from the Upper Pliensbachian(CH18-050)

Fig. 4 a δ18O versus δ13C values of well-preserved oysters and brachiopods compared to samples of poorly preserved oysters, sediment, and cement; b–d δ18O versus δ13C values of the three oysters used for high-resolution analyses. The δ18O and δ13C values of specimen CH18-006 are strongly correlated and believed to show the influence of diagenetic alteration. In contrast, the results of CH18-017 and CH18-030 show no correlation and are believed to represent original compositions

The present dataset of stable isotope (δ18O, δ13C) values of well-preserved oyster and brachiopod shells as well as of altered oysters, sediment, and cement samples is illustrated in Fig. 4a. While the δ13C values of the individual samples are comparatively similar, the δ18O values range widely. In particular, the samples of poorly preserved oysters and cement show very low δ18O values (less than -8‰). The Spearman correlation coefficient indicates no significant correlation between δ18O and δ13C values of the wellpreserved brachiopods (rs=-0.08; p=0.74). The values of seemingly well-preserved oyster shells are moderately correlated (rs=0.56; p=0.05). This correlation is mainly due to specimen CH18–006 with a comparatively low δ18O value.If this sample is excluded from the dataset, the correlation is statistically not significant with a p-value above 5%.

Three oysters were selected for high-resolution stable isotope analyses (Fig. 2c–e). The data of specimen CH18–006 show a strong correlation between δ18O and δ13C values(Fig. 4b; rs=0.89; p ＜0.05), which is believed to be due to the influence of diagenetic alteration. Results of this specimen were therefore not used for temperature reconstructions. Stable isotope data of the other two specimens CH18-017 and CH18-030 show no significant correlation (Fig. 4c, d) and are therefore believed to represent original conditions.

In summary, cathodoluminescence microscopy and the evaluation of a potential correlation between δ18O and δ13C values indicate a good preservation of the shells used for temperature reconstructions (except for specimen CH18-006).

5 Results of the stable isotope analyses

The results of the stable isotope (δ13C, δ18O) analyses of well-preserved oysters and brachiopods are listed in Table 1. In addition, Table 1 lists the data from poorly preserved oysters, sediment and cement samples. The data obtained from high-resolution analyses of the three oyster shells can be found in Table 2.

Figure 5 illustrates the stratigraphic position of the sampled shells and shows the δ13C and δ18O values of the well-preserved oysters and brachiopods. While the data for the Sinemurian and Pliensbachian are retrievedfrom brachiopods, the Toarcian is mostly represented by oyster shells. The δ13C values show an overall slight increase from the Sinemurian to the Toarcian.They fluctuate around an average of 1.37‰ in the Upper Sinemurian(range:0.72‰to 1.97‰).Although the data for the Pliensbachian is limited,δ13C values seem to be higher around an average of 1.68‰(range:1.48‰to 2.01‰).Apart from one outlier (specimen CH18-046: -0.05‰), the highest δ13C values have been recorded in the upper Toarcian with an average of 2.00‰(range:1.76‰to 2.59‰).

Table 1 Results of the stable isotope (δ13C, δ18O) analyses of oysters and brachiopods from the Lower Jurassic of northern Chile. In addition to well-preserved shells, some luminescent and altered oyster shells as well as sediment and cement samples were analyzed. Please note that the high-resolution analysis of specimen CH18-006 indicated a diagenetic alteration of this shell, which has consequently not been used for further interpretations

Table 1 Results of the stable isotope(δ13C,δ18O) analyses of oysters and brachiopods from the Lower Jurassic of northern Chile. In addition to well-preserved shells, some luminescent and altered oyster shells as well as sediment and cement samples were analyzed. Please note that the high-resolution analysis of specimen CH18-006 indicated a diagenetic alteration of this shell,which has consequently not been used for further interpretations(Continued)

The δ18O values show some fluctuations through time.In the late Sinemurian,the values range around an average of -3.22‰ (range: -4.09‰ to -2.31‰) with the highest values close to the Sinemurian–Pliensbachian boundary(Fig. 5). In the Pliensbachian, there are two outliers with comparatively low values (CH18-032: -4.88‰; CH18-034:-3.86‰). These two shells (which show δ13C values similar to the other shells) are from section 2 near Potrerillos around 280 km north of all other analyzed shells (Fig. 1).The other two Pliensbachian shells from a higher stratigraphic interval near the boundary to the Toarcian show comparatively high δ18O values around an average of-2.46‰ (CH18-001: -3.07‰; CH18-002: -1.85‰). In the late Toarcian, the data fluctuate around an average of-2.92‰(range:-3.69‰to-2.46‰).This interval also contains specimen CH18-006, which has been excluded from interpretation, as it has shown signs of diagenetic alteration during high-resolution analysis(see above).

In summary, the δ13C and δ18O values show differing trends through time. While the δ13C values show an overall steady increase through time, the changes in the δ18O values are somewhat more complex. Over most of the late Sinemurian, the δ18O values are below -3‰. Just below the Sinemurian–Pliensbachian boundary, the average δ18O value rises above -3‰.The highest δ18O value in the entire dataset is reached at the end of the Pliensbachian. The δ18O values in the late Toarcian are again lower but not as low as in the late Sinemurian (Fig. 5).

The results of the three oysters used for highresolution stable isotope analyses are shown in Fig. 6. As described already above, the δ13C and δ18O values of specimen CH18-006 are correlated strongly (Fig. 6a).Throughout the shell, the δ18O values decrease to lower values with a few fluctuations. This pattern is almost perfectly mirrored in the δ13C values. These parallel trends are believed to be caused by diagenetic alteration and specimen CH18-006 was excluded from further interpretation. In contrast, the δ13C and δ18O values of the oysters CH18-017 and CH18-030 indicate an original stable isotope composition. Both shells show cyclic patterns in their δ18O values, while the δ13C values show very little and unrelated fluctuations (Fig. 6b, c). The δ18O values of specimen CH18-017 vary around an average of -2.82‰ (range: -3.18‰ to -2.47‰). The data show two minima and three maxima. The δ18O values of specimen CH18-030 show a somewhat similar pattern,but are slightly higher with an average of -2.28‰ (range:-2.56‰ to -1.88‰). The cycles are less clear compared to CH18-017, but three minima and two maxima can be discerned (most easily by using a two-point running average).

6 Discussion

6.1 Palaeoenvironment

Aberhan (1992) assigned three different depositional environments to the Lower Jurassic strata of the study area: (1) a coastal zone dominated by siliciclasticsediments with sedimentary structures pointing to a setting above or just below the fair-weather wave-base(mainly sand- or siltstones), (2) an intermediate zone dominated by carbonates, but still showing a considerable clay and silt content (mainly bioclastic wacke- to packstones), and (3) a distal zone with low sedimentation rates characterized by very fine-grained terrigenous siliciclastics and carbonate mud (mainly mud- to wackestones and marl). While zones 1 and 2 have comparatively high faunal diversities, zone 3 is characterized by very low diversities (possibly caused by low oxygen availability and soupy substrates;Aberhan 1992).

Table 2 Results of the high-resolution stable isotope (δ13C,δ18O) analyses of three oysters with a late Toarcian age

Table 2 Results of the high-resolution stable isotope (δ13C,δ18O) analyses of three oysters with a late Toarcian age(Continued)

Shells used in the present study are from the deeper part of zone 1 and zone 2 (compare Fig. 5). Some shells come from siltstones or mixed carbonate–siliciclastics(partly with ferruginous ooids), which represent somewhat deeper water around or just below the fair-weather wave-base (zone 1 of Aberhan 1992). Most fossils come from carbonates (zone 2 of Aberhan 1992), which are occasionally marly or are developed as shell beds possibly formed by storm action (the repeated occurrence of higher energy conditions such as storms has also been suggested by Fantasia et al. 2018). No direct correlation between the results of the stable isotope analyses and the facies types (= water depth) can be seen (Fig. 5) and it can be assumed that the analyzed shells were formed above the thermocline and therefore represent seasurface water conditions.

Fig. 5 Lithologs (modified after Aberhan 1992) with the position of the individual samples and the results of the stable isotope (δ13C, δ18O)analyses of well-preserved oysters and brachiopods. Note the inverted scale for δ18O values

Apart from the water temperature, the stable isotope(δ18O) composition of shells in marine environments can be influenced by salinity fluctuations. A deviation from the normal marine salinity due to freshwater influx(rain fall or river discharge) or strong evaporation (e.g.,in an enclosed lagoon) can considerably alter the oxygen isotope composition of the ocean water and consequently of the shells. As described above, the fossil material used in the present study represents a shallow water depth slightly above the fair-weather wave-base or just below, but still above the thermocline and within reach of occasional stronger storms.No fossils have been used from the cross-bedded sandstones which formed in a very shallow setting.Freshwater influx via stronger rain fall affecting the uppermost water layers can therefore be excluded. Strong deviations from normal marine salinity are also unlikely considering that the fossils collected in the stratigraphic interval studied include stenohaline forms such as corals, brachiopods, crinoids,and ammonites (Fig. 5; Aberhan 1992; compare also Fantasia et al. 2018).

6.2 Carbon isotope record

The δ13C data (Fig. 5) indicate a trend towards higher values, however, it should be noted that the Sinemurian and Pliensbachian samples are represented by brachiopods,while the Toarcian data come mainly from oysters. The δ13C values of oysters are consistently higher than those of the brachiopods (Fig. 4a), even though the δ18O values are comparable.Vital effects causing higher δ13C values of oysters might therefore explain the recorded “trend” from the Sinemurian to the Toarcian at least in parts. This is in broad agreement with the dataset of Korte et al. (2015) for the Early Jurassic, which shows mostly higher δ13C values in bivalves than in brachiopods of the same age. Other studies, however, show that co-occurring bivalves and brachiopods have very similar δ13C values (e.g., Alberti et al.2017,2019).Further data of different faunal groups is therefore necessary to judge whether there actually is a trend in the δ13C values through time or whether the difference is merely caused by vital effects.Apart from this problem, the absolute δ13C values are comparable to those of shells from Europe(e.g.,Korte et al.2015).

Fig. 6 Results of the high-resolution stable isotope (δ13C, δ18O)analyses of three oysters with a late Toarcian age. While the data of specimen CH18-006 (a) shows evidence of diagenetic alteration, the values of CH18-017 (b) and CH18-030 (c) are assumed to reflect original conditions representing seasonal temperature fluctuations.In addition to the overall average for each dataset, a two-point running average is shown.Note the inverted scale for δ18O values

Lower Jurassic rocks include a number of prominent carbon isotope excursions including one during the Toarcian Oceanic Anoxic Event (e.g., Jenkyns 1988; Hesselbo et al. 2000a, 2007; Pálfy and Smith 2000; van de Schootbrugge et al. 2005a; Gómez et al. 2008; Suan et al.2008, 2011; Littler et al. 2010; Korte and Hesselbo 2011;Silva et al. 2011; Caswell and Coe 2012; Riding et al.2013; Kemp and Izumi 2014; Krencker et al. 2014; Ait-Itto et al. 2017; Arabas et al. 2017; Bougeault et al. 2017;Menini et al. 2019; Mercuzot et al. 2019). The TOAE has also been reported from Argentina (Al-Suwaidi et al.2010, 2014, 2016; Mazzini et al. 2010) and Chile (Fantasia et al. 2018). Fantasia et al. (2018) reported several negative carbon isotope excursions based on analyses of whole-rock samples (i.e. within the Margaritatus Zone,at the Pliensbachian–Toarcian boundary, and during the TOAE). These excursions are not visible in the present δ13C values of fossil shells (Fig. 5), because they were of relatively short duration and/or the stratigraphic intervals where they are most pronounced are not represented by fossils (compare Boulila and Hinnov 2017).The variation in the δ13C values of fossil shells is relatively small compared to that recorded by Fantasia et al.(2018) from whole-rock analyses. The only exception is specimen CH18-046 from the upper Toarcian, which shows a much lower δ13C value than all other shells analyzed. However, this shell has a considerably younger age than the previously reported negative carbon isotope excursions, and since other shells from the same horizon show higher values, it might be less reliable. It can be therefore presumed that the presented dataset represents background conditions prevalent over most of the Early Jurassic and not exceptional events.

6.3 Water temperatures

Considering that the shells formed above the thermocline and in normal marine salinity conditions, the recorded δ18O values can be translated into sea-surface watertemperatures. In accordance with the majority of other studies on Jurassic stable isotopes, the equation of Anderson and Arthur(1983)with a δ18O value of-1‰VSMOW for seawater during shell precipitation (as suggested for an ice-free Jurassic world; Shackleton and Kennett 1975) was used. The reconstructed water temperatures show several fluctuations through time (Fig. 5). During most of the late Sinemurian, temperatures were relatively high (average:27.0°C). Just before the Sinemurian–Pliensbachian transition,temperatures dropped to an average of 24.3°C.At the Quebrada El Asiento(section 2)near Potrerillos,the results point to very warm conditions during the Pliensbachian(average: 31.4°C). Further south near El Transito (sections 4, 5, 6) temperatures were much lower (average: 22.3°C).The lowest temperature in the entire dataset was recorded for the latest Pliensbachian (specimen CH18-002: 19.6°C).No data is available for the early Toarcian,but results indicate again higher temperatures for the late Toarcian (average:24.4°C)comparable to those of the late Sinemurian.

The data retrieved by high-resolution stable isotope(δ18O) analyses of the two oysters CH18-017 and CH18-030 can be translated into seasonal temperature fluctuations (Fig. 6b, c). Both shells show similar conditions in the late Toarcian reflecting temperature cycles during a period of around 3 years in the life time of the oysters. If explained only by temperature, the δ18O values of specimen CH18-017 record an average temperature of 23.9°C with a minimum of 22.4°C and a maximum of 25.6°C. The δ18O values of specimen CH18-030 record a slightly lower average temperature of 21.5°C with a minimum of 19.7°C and a maximum of 22.8°C. Since the oyster shells are not very large (Fig. 2d, e), the temporal resolution of samples is not very high for both shells with a year being represented by less than ten samples. Therefore true temperature maxima and minima occurring over the life time of the oysters are not recorded. Each value instead represents an average temperature of a period of one to 2 months. The reconstructed seasonality of slightly more than 3°C for both shells is therefore only a minimum estimate.

6.4 Comparisons with previously published data

Fig. 7 The Early Jurassic δ18O data from a Chile (for legend see Fig. 5) compared to the large datasets of b Martinez and Dera (2015) and c Korte et al. (2015) from the northwestern Tethys for the same time interval. The dataset of Martinez and Dera (2015) consists of 1662 individual belemnite samples (from Hesselbo et al. 2000b; McArthur et al. 2000; Jenkyns et al. 2002; Bailey et al. 2003; Gómez et al. 2008; Metodiev and Koleva-Rekalova 2008; Dera et al. 2009b, 2011; Price 2010; Korte and Hesselbo 2011; Armendáriz et al. 2012; Li et al. 2012; Harazim et al. 2013;Riding et al. 2013). The dataset of Korte et al. (2015) consists of 823 individual belemnite samples in addition to the results of 140 bivalves and 46 brachiopods (integrating their new data with that of McArthur et al. 2000; Cresta et al. 2001; Jenkyns et al. 2002; Suan et al. 2008; Price 2010;Korte and Hesselbo 2011; Li et al.2012;Ullmann et al.2014).Note the inverted scale for δ18O values

So far, reconstructions of absolute seawater temperatures for the Early Jurassic using the stable isotope(δ18O) composition of fossil hard parts are almost completely restricted to European and neighboring localities(e.g., Sælen et al. 1996; Hesselbo et al. 2000b; McArthur et al. 2000, 2007; Bailey et al. 2003; Nori and Lathuilière 2003; Rosales et al. 2004a, 2004b, 2018; van de Schootbrugge et al. 2005b; Gómez et al. 2008, 2016; Metodiev and Koleva-Rekalova 2008; Suan et al. 2008, 2010; Dera et al. 2009a, 2009b, 2011; Gómez and Arias 2010; Price 2010; Korte and Hesselbo 2011; Armendáriz et al. 2012;Li et al. 2012; Harazim et al. 2013; Teichert and Luppold 2013; Krencker et al. 2014, 2015; Metodiev et al. 2014;Korte et al. 2015; Martinez and Dera 2015; Price et al.2016; Arabas et al. 2017). In these regions of the northwestern Tethys, previously published results show decreasing temperatures from the late Sinemurian throughout the Pliensbachian with a temperature minimum in the late Pliensbachian and around the Pliensbachian–Toarcian boundary. Temperatures then strongly increased in the early Toarcian with the development of the Toarcian Oceanic Anoxic Event, but remained relatively high throughout the Toarcian until they again decreased strongly into the Aalenian (Korte et al. 2015).Even though most authors agree that temperatures increased dramatically throughout the early Toarcian,some authors have questioned the reconstructed amplitude of this change and attributed at least part of it to a parallel change in the δ18O values of the seawater (e.g.,by the disappearance of polar ice shields or an increased freshwater runoff in the Toarcian; Sælen et al. 1996;Rosales et al. 2004b; Suan et al. 2008; Dera et al. 2009b,2011; Dera and Donnadieu 2012; Korte et al. 2015).

Temperature reconstructions for regions outside Europe especially those based on stable isotope (δ18O) analyses are very scarce(Volkheimer et al.2008).Bowen(1963)analyzed eight poorly dated belemnites from the Lower Jurassic of Argentina. After dismissing one specimen due to unlikely analytical results,his δ18O results vary between-0.1‰and-2.9‰(equivalent to temperatures of 12.4°C to 24.3°C using a δ18O for seawater of-1‰and the equation of Anderson and Arthur 1983). Interestingly, Bowen(1963)reconstructed a seasonality of 3.3°C to 3.7°C based on his belemnites,which is quite comparable to the results of the present study. However, since his data is only known to be Lower Jurassic,it can hardly be compared to the current dataset considering the likelihood of prominent climate changes during this time interval.

The current study is the first to present Early Jurassic absolute water temperatures from South America based on material which could be assigned to ammonite zones.Figure 7 interprets the presented data from Chile and compares it to the datasets compiled by Martinez and Dera (2015) and Korte et al. (2015) for the northwestern Tethys. Martinez and Dera (2015) and Korte et al.(2015) compiled high-resolution stable isotope (δ18O)records from European localities which show almost identical patterns. The most striking feature is the Late Pliensbachian Cooling Event (LPCE) and the early Toarcian Oceanic Anoxic Event (TOAE). Even though the South American dataset is limited so far, it seems that it broadly follows the temperature pattern already described from Europe. The water temperatures show a general decrease from the late Sinemurian until a minimum in the latest Pliensbachian Spinatum Zone,followed again by a warming until the late Toarcian.Unfortunately, no material from the lower Toarcian was available in the collection to document the possible presence of the early Toarcian warming. The high temperatures at Potrerillos during the Pliensbachian might indicate a special situation in this part of the basin (possibly the influence of freshwater), but more data is needed in order to come to a final conclusion.

On average, the reconstructed absolute water temperatures in Chile are slightly higher than those of the northwestern Tethys (Fig. 7). This can be attributed to a higher palaeolatitude of the European localities in the northern hemisphere compared to Chile in the southern hemisphere. At the same time, it has to be noted that the majority of data from Europe is based on belemnites,which are known to commonly record lower temperatures than co-occurring brachiopods and bivalves (e.g.,Dera et al. 2011; Alberti et al. 2012a, 2019). The cause of this difference might include different life habits between the faunal groups, but is still strongly debated in literature (Fürsich et al. 2005; Mutterlose et al. 2010; Alberti et al. 2012a, 2012b; Price et al. 2015; Dera et al. 2016;Hoffmann et al. 2016).

7 Conclusions

The dataset from the Andean Basin in northern Chile is the first stable isotope record of well-dated brachiopod and oyster shells from the Lower Jurassic of South America. The recorded δ18O values point to warm conditions in the late Sinemurian (average: 27.0°C), a slight cooling just before the Sinemurian–Pliensbachian transition (average: 24.3°C), lowest temperatures in the latest Pliensbachian (19.6°C), and again warmer temperatures in the late Toarcian(average:24.4°C).Thereby,the dataset supports the global nature of the Late Pliensbachian Cooling Event. More material is necessary to detect the presumed global warming during the early Toarcian Oceanic Anoxic Event, which was not represented by material in the collection analyzed. The recorded δ13C values represent background conditions prevalent over most of the Early Jurassic instead of short-term intervals of carbon-cycle perturbations. High-resolution stable isotope (δ18O) analyses of oyster shells point to a seasonality in northern Chile of more than 3°C in the late Toarcian.

Abbreviations

cl: Clay; g: Gravel; LPCE: Late pliensbachian cooling event; s: Sand;si: Silt;TOAE: Toarcian oceanic anoxic event; VPDB: Vienna pee dee belemnite;VSMOW: Vienna standard mean ocean water

Acknowledgement

The authors thank Winfried Werner for enabling access to the collections of the Bayerische Staatssammlung für Paläontologie und Geologie in Munich,Germany. Birgit Leipner-Mata and Michael Joachimski allowed the use of laboratories at the University of Erlangen-Nürnberg, Germany. Horacio Parent(Rosario, Argentina)and Axel von Hillebrandt (Berlin,Germany) provided copies of scientific articles. The authors thank the Editor-in-Chief Zengzhao Feng and three anonymous reviewers for their support and constructive comments.

Authors’ contributions

MA designed research, performed research, acquired data, and analyzed data. NA performed the stable isotope analyses. MA, FTF, and NA drafted the manuscript. All authors read and approved the final manuscript.

Funding

The study was supported financially by the German Research Foundation(DFG; project AL 1740/3–1).

Availability of data and materials

All data generated and analyzed during this study are included in this published article or in the cited references. The studied material is stored in the collections of the Bayerische Staatssammlung für Paläontologie und Geologie in Munich, Germany.

Competing interests

The authors declare that they have no competing interests.

Author details

1Institut für Geowissenschaften, Christian-Albrechts-Universität zu Kiel,Ludewig-Meyn-Straße 10, 24118 Kiel, Germany.2GeoZentrum Nordbayern,Fachgruppe PaläoUmwelt,Friedrich-Alexander-Universität Erlangen-Nürnberg,Loewenichstraße 28, 91054 Erlangen, Germany.3Leibniz Laboratory for Radiometric Dating and Stable Isotope Research,Christian-Albrechts-Universität zu Kiel, Max-Eyth-Straße 11, 24118 Kiel,Germany.

Received: 12 July 2019 Accepted: 24 October 2019