1. Introduction
The age of the Chicxulub impact has been a contentious issue ever since the discovery of the impact structure on the Yucatan peninsula and its initial tentative link to the Cretaceous–Tertiary boundary (KTB; also called K–P or K–Pg) mass extinction, based mainly on the stratigraphic position of impact spherules relative to the KTB in Haiti and NE Mexico (e.g. Maurrasse & Sen, Reference Maurrasse and Sen1991; Smit et al. Reference Smit, Montanari, Swinburne, Alvarez, Hildebrand, Margolis, Claeys, Lorie and Asaro1992; Stinnesbeck et al. Reference Stinnesbeck, Barbarin, Keller, Lopez-Oliva, Pivnik, Lyons, Officer, Adatte, Graup, Rocchia and Robin1993; Smit, Reference Smit1999). Over the past 20 years, data from over 80 localities throughout the Caribbean, Central and North America (Fig. 1) have provided variable results with impact spherule layers in latest Maastrichtian and early Danian sediments (Keller et al. Reference Keller, Stinnesbeck, Adatte, Holland, Stueben, Harting, De Leon and De La Cruz2003a , Reference Keller, Stinnesbeck, Adatte and Stueben a , Reference Keller, Adatte, Berner, Harting, Baum, Prauss, Tantawy and Stueben2007, Reference Keller, Adatte, Berner, Pardo and Lopez-Oliva2009, Reference Keller, Keller and Adatte2011). To understand the age and nature of spherule deposition we must comprehend this regional stratigraphic variation and consider all evidence to determine what really happened over 65 Ma ago.
The strongest age claims for a KTB age for the Chicxulub impact are based on a select few deep-sea localities where a single thin impact spherule layer is reported directly between the Cretaceous and Paleocene. These localities include the Bass River Core of New Jersey, NW Atlantic Site 1049, Caribbean Site 1001 and SW Atlantic Site 1259 (Olsson et al. Reference Olsson, Miller, Browning, Habib and Sugarmann1997; Sigurdsson et al. Reference Sigurdsson, Leckie, Acton, Miller and Maddox1997; Norris, Huber & Self-Trail Reference Norris, Huber and Self-Trail1999; Norris et al. Reference Norris, Firth, Blusztajn and Ravizza2000; Klaus et al. Reference Klaus, Norris, Kroon and Smit2000; Martinez-Ruiz et al. Reference Martinez-Ruiz, Ortega-Huertas, Palomo-Delgado, Smit, Kroon, Norris and Klaus2001; Huber, MacLeod & Norris, Reference Huber, MacLeod, Norris, Koeberl and MacLeod2002; MacLeod et al. Reference MacLeod, Whitney, Huber and Koeberl2007). This single spherule layer is commonly cited as proof that the Chicxulub impact caused the KTB mass extinction. In this scenario, the presence of multiple spherule layers in most other localities is attributed to an impact generated tsunami that caused gravity flows, slumps or liquefaction sucking spherules from the KTB into Maastrichtian sediments (Schulte et al. Reference Schulte, Alegret and Arenillas2010).
The simplicity of a thin spherule layer separating Maastrichtian and Danian sediments is appealing as an age indicator but difficult to reconcile with regional data. Are these sections complete across the KTB as claimed? If not, is there a pattern of spherule erosion and redeposition? To answer these questions a high-resolution quantitative planktic foraminiferal biostratigraphy is required, which does not exist for these localities. This type of data can determine not only the presence or absence of biozones, but also whether sedimentation is complete in each biozone based on the abrupt onset or termination of species populations (e.g. the sudden appearance or disappearance of species that dominate the assemblage). For these reasons we examined the KTB sequences at Bass River, New Jersey, and ODP Sites 1049A, 1049C, 1050C, 999B, 1001B and 1259B (Fig. 1) based on: (1) high-resolution planktic foraminiferal biostratigraphy; (2) quantitative species abundance data as indicators of biozone truncation and hiatuses; (3) completeness of the KT boundary transition; (4) stable isotopes; and (5) stratigraphic position of the spherule layer relative to the KTB. The KTB was identified based on the five criteria of the El Kef stratotype: KTB characteristic clay layer, δ13C shift, iridium anomaly, mass extinction and evolution of Danian species (Cowie, Ziegler & Remane, Reference Cowie, Ziegler and Remane1989; Keller, Li & MacLeod, Reference Keller, Li and MacLeod1995; Remane et al. Reference Remane, Keller, Hardenbol and Ben Haj Ali1999).
2. Chicxulub impact controversy – background
The ongoing KT impact controversy is now over three decades old and was recently summarized from very different viewpoints in Schulte et al. (Reference Schulte, Alegret and Arenillas2010) and Keller (Reference Keller, Keller and Adatte2011a ). At the core of this controversy is the age of the Chicxulub impact and the cause of the KTB mass extinction. The most critical evidence comes from NE Mexico, Texas and the Chicxulub impact crater core Yaxcopoil-1. In NE Mexico and Texas impact spherule layers are at the base of a sandstone complex that infills submarine channels along the slope of the Gulf of Mexico (Sohl et al. Reference Sohl, Martinez, Salmerón-Ureña, Soto-Jaramillo and Salvador1991) and at the margin of the Western Interior Seaway (Keller & Adatte, Reference Keller and Adatte2011), respectively. The controversy revolves over the origin of this sandstone complex. Was it deposited as the result of an impact-generated tsunami at the KTB (Smit et al. Reference Smit, Roep, Alvarez, Montanari, Claeys, Grajales-Nishimura and Bermudez1996; Schulte et al. Reference Schulte, Alegret and Arenillas2010), or was deposition long term and associated with erosion, transport and redeposition during the sea-level fall near the end of the Maastrichtian (Keller et al. Reference Keller, Lopez-Oliva, Stinnesbeck and Adatte1997, Reference Keller, Adatte, Berner, Pardo and Lopez-Oliva2009; Adatte, Stinnesbeck & Keller, Reference Adatte, Stinnesbeck, Keller, Ryder, Fastovsky and Gartner1996; Adatte, Keller & Stinnesbeck, Reference Adatte, Keller and Stinnesbeck2002).
The sandstone complex is deposited in submarine channels at depths of 500–1000 m in NE Mexico and in shallow inner neritic environments (<50 m depth) in Texas. Three units are generally present: Unit 1 at the base consists of one to three layers of spherule-rich, sandy and glauconitic sediments enriched in reworked shallow water debris (e.g. plant remains, shallow water foraminifera; Keller et al. Reference Keller, Lopez-Oliva, Stinnesbeck and Adatte1997; Alegret, Molina & Thomas, Reference Alegret, Molina and Thomas2001; Alegret et al. Reference Alegret, Arenillas, Arz, Liesa, Melendez, Molina, Soria and Thomas2002). Thin sandstone or limestone layers separate the spherule layers (Keller et al. Reference Keller, Adatte, Berner, Harting, Baum, Prauss, Tantawy and Stueben2007, Reference Keller, Adatte, Berner, Pardo and Lopez-Oliva2009; Adatte, Keller & Baum, Reference Adatte, Keller, Baum, Keller and Adatte2011). Unit 2 in NE Mexico consists of a thick sandstone with two to three thin volcanic-rich (zeolite) layers (Adatte, Stinnesbeck & Keller, Reference Adatte, Stinnesbeck, Keller, Ryder, Fastovsky and Gartner1996). Unit 3 at the top consists of alternating sandy shale and marl layers with abundant burrows (Ekdale & Stinnesbeck, Reference Ekdale and Stinnesbeck1998). The KTB is above this unit.
The number of spherule layers in the sandstone complex, or in the Danian above and Maastrichtian below, their mineralogical and fossil contents and nature of deposition are variable. For example, at El Penon and El Mimbral in NE Mexico two spherule layers at the base of the sandstone complex are separated by a 15–20 cm thick limestone bed with occasional large burrows infilled with spherules, in other localities this limestone bed is truncated by erosion (e.g. La Lajilla, Fig. 1; Lopez-Oliva & Keller, Reference Lopez-Oliva and Keller1996; Adatte, Stinnesbeck & Keller, Reference Adatte, Stinnesbeck, Keller, Ryder, Fastovsky and Gartner1996; Keller et al. Reference Keller, Stinnesbeck, Adatte and Stueben2003b ; Stueben et al. Reference Stueben, Kramar, Harting, Stinnesbeck and Keller2005). At El Penon an additional spherule layer is interbedded in upper Maastrichtian marls up to 4–5 m below the sandstone complex (Keller et al. Reference Keller, Adatte, Berner, Pardo and Lopez-Oliva2009), as also observed at Loma Cerca (Schulte et al. Reference Schulte, Stinnesbeck, Stueben, Kramar, Berner, Keller and Adatte2003). In the Brazos sections of Texas two to three spherule layers are separated by thin sandstones at the base of the sandstone complex and up to 1 m of latest Maastrichtian sediments are between the KTB and the top of the sandstone complex (Gale, Reference Gale2006; Schulte et al. Reference Schulte, Speijer, Mai and Kontny2006; Prauss, Reference Prauss2009; Keller et al. Reference Keller, Adatte, Berner, Harting, Baum, Prauss, Tantawy and Stueben2007, Reference Keller, Abramovich, Adatte, Berner, Keller and Adatte2011; Adatte, Keller & Baum, Reference Abramovich, Keller, Berner, Cymbalista, Rak, Keller and Adatte2011). In the Chicxulub impact crater core Yaxcopoil-1 a 50 cm thick limestone separates the impact breccia from the KTB (Arz, Alegret & Arenillas, Reference Arz, Alegret and Arenillas2004; Keller, Adatte & Stinnesbeck, Reference Keller, Adatte and Stinnesbeck2004a ; Keller et al. Reference Keller, Adatte, Stinnesbeck, Rebolledo-Vieyra, Urrutia Fucugauchi, Kramar and Stueben2004b ; Smit, van der Gaast & Lustenhouwer, Reference Smit, van der Gaast and Lustenhouwer2004).
In central and southern Mexico (e.g. Tlaxcalantongo/La Ceiba, Coxquihui, Bochil, Guayal) impact spherules are found at or just below the KTB (Arz et al. Reference Arz, Arenillas, Soria, Alegret, Grajales-Nishimura, Liesa, Meléndez, Molina and Rosales2001; Arenillas et al. Reference Arenillas, Arz, Grajales-Nishimura, Murillo-Muneton, Alvarez, Camargo-Zanguera, Molina and Rosales-Dominguez2006) as well as in the early Danian zone P1a (Stinnesbeck et al. Reference Stinnesbeck, Keller, Schulte, Stueben, Berner, Kramar and Lopez-Oliva2002; Keller et al. Reference Keller, Stinnesbeck, Adatte and Stueben2003b ; Keller, Reference Keller, Evans, Horton and King2008; Fig. 1). In Guatemala and Belize, all impact spherules are in early Danian zone P1a sediments (Keller et al. Reference Keller, Stinnesbeck, Adatte, Holland, Stueben, Harting, De Leon and De La Cruz2003a ). In Cuba impact spherules are reported at the top of a conglomerate (Alegret et al. Reference Alegret, Arenillas, Arz, Diaz, Grajales-Nishimura, Melendez, Molina, Rojas and Soria2005), as well as a 2 cm thick layer in the overlying early Danian zone P1a (Keller, Reference Keller, Evans, Horton and King2008). Four sections analyzed in Haiti also revealed impact spherules in early Danian zone P1a sediments (Keller et al. Reference Keller, Adatte, Stinnesbeck, Stueben and Berner2001, Reference Keller, Stinnesbeck, Adatte, Holland, Stueben, Harting, De Leon and De La Cruz2003b ). From the North Atlantic (Bass River core, ODP Sites 1049, 1001B and 1259) a thin spherule layer is reported at the KTB (Olsson et al. Reference Olsson, Miller, Browning, Habib and Sugarmann1997; Sigurdsson et al. Reference Sigurdsson, Leckie, Acton, Miller and Maddox1997; Norris et al. Reference Norris, Kroon, Klaus, Baez and Kapitan-White1998; Norris, Huber & Self-Trail, Reference Norris, Huber and Self-Trail1999; MacLeod et al. Reference MacLeod, Whitney, Huber and Koeberl2007).
2.a. Nature of impact spherule deposition
The earliest and most commonly advocated scenario for the sandstone complex and impact spherules as impact generated tsunami deposits precisely at the KTB is based on its proximity to the KTB mass extinction and the KTB iridium anomaly, and the belief that no other major event occurred or could have caused the mass extinction (e.g. Bourgeois et al. Reference Bourgeois, Hansen, Wiberg and Kauffman1988; Smit et al. Reference Smit, Montanari, Swinburne, Alvarez, Hildebrand, Margolis, Claeys, Lorie and Asaro1992, Reference Smit, Roep, Alvarez, Montanari, Claeys, Grajales-Nishimura and Bermudez1996; Alegret et al. Reference Alegret, Arenillas, Arz, Liesa, Melendez, Molina, Soria and Thomas2002, Reference Alegret, Arenillas, Arz, Diaz, Grajales-Nishimura, Melendez, Molina, Rojas and Soria2005; Arenillas et al. Reference Arenillas, Arz, Grajales-Nishimura, Murillo-Muneton, Alvarez, Camargo-Zanguera, Molina and Rosales-Dominguez2006; Schulte et al. Reference Schulte, Stinnesbeck, Stueben, Kramar, Berner, Keller and Adatte2003, Reference Schulte, Speijer, Mai and Kontny2006, Reference Schulte, Alegret and Arenillas2010). In this scenario, the KTB clay and global iridium anomaly is attributed to the Chicxulub impact, whereas pyrite, glauconite, and sanidine spherules commonly present in the KTB clay are interpreted as altered Chicxulub impact spherules (Smit, Reference Smit1999), although there is no evidence for these claims. Various sedimentological, mineralogical, stratigraphical and palaeontological data argue against this KTB impact-tsunami scenario. For example, in NE Mexico the multiple spherule layers separated by limestone, containing rare burrows infilled with spherules, could not have been deposited within hours to days (Keller et al. Reference Keller, Lopez-Oliva, Stinnesbeck and Adatte1997, Reference Keller, Adatte, Berner, Pardo and Lopez-Oliva2009). The shallow water debris and benthic foraminifera in the spherule layers of deep-water deposits (Alegret, Molina & Thomas, Reference Alegret, Molina and Thomas2001; Alegret et al. Reference Alegret, Arenillas, Arz, Liesa, Melendez, Molina, Soria and Thomas2002; Alegret & Thomas, Reference Alegret and Thomas2004) require prior spherule deposition in shallow water followed by erosion and transport. The two to three correlative zeolite layers in the middle unit of the sandstone complex require depositional events separated by time (Adatte, Stinnesbeck & Keller, Reference Adatte, Stinnesbeck, Keller, Ryder, Fastovsky and Gartner1996). The burrowed intervals in the alternating sandy shale and marl layers of the upper part of the sandstone complex indicate colonization of the seafloor during deposition (Ekdale & Stinnesbeck, Reference Ekdale and Stinnesbeck1998). Similarly, the truncated burrows across the hummocky sandstones in the Brazos sections of Texas and in the overlying layers mark multiple recolonization of the seafloor (Gale, Reference Gale2006; Adatte, Keller & Baum, Reference Adatte, Keller, Baum, Keller and Adatte2011). At the base of the Brazos sandstone complex the presence of lithified clasts with impact spherules, including some in desiccation cracks, marks a history of prior spherule deposition followed by lithification, emersion, erosion, transport and redeposition (Adatte, Keller & Baum, Reference Adatte, Keller, Baum, Keller and Adatte2011; Keller et al. Reference Keller, Adatte, Berner, Harting, Baum, Prauss, Tantawy and Stueben2007, Reference Keller, Abramovich, Adatte, Berner, Keller and Adatte2011).
All these characteristics of the sandstone complex suggest long-term deposition interrupted by gravity flows and channelized downslope transport possibly related to climate and sea level changes. Schulte et al. (Reference Schulte, Alegret and Arenillas2010, p. 1215) deny the presence of most of these well-documented features and maintain rapid deposition by tsunami and gravity flows. Incongruously they claim that ‘the presence of shallow water foraminifera in the clastic unit contradicts a long-term depositional sequence; if in situ, their presence requires unrealistically rapid relative sea-level changes of >500 m.’ However, there never was a claim of in situ deposition of shallow benthic foraminifera in bathyal deposits, but rather downslope transport (Alegret, Molina & Thomas, Reference Alegret, Molina and Thomas2001; Alegret et al. Reference Alegret, Arenillas, Arz, Liesa, Melendez, Molina, Soria and Thomas2002; Keller et al. Reference Keller, Stinnesbeck, Adatte and Stueben2003b , Reference Keller, Adatte, Berner, Harting, Baum, Prauss, Tantawy and Stueben2007, Reference Keller, Adatte, Baum and Berner2008). In contrast, to NE Mexico and Texas, impact spherules are generally found in early Danian sediments (zone P1a) of central and southern Mexico, Guatemala, Belize, Haiti and Cuba. This can be easily explained by reworking and redeposition of spherule deposits by intensified Gulf Stream circulation and sea level changes (e.g. Keller et al. Reference Keller, Lyons, MacLeod and Officer1993, Reference Keller, Stinnesbeck, Adatte, Holland, Stueben, Harting, De Leon and De La Cruz2003a ; Watkins & Self-Trail, Reference Watkins and Self-Trail2005; Keller, Reference Keller, Evans, Horton and King2008). These regions are rarely mentioned in the Chicxulub controversy.
2.b. Age of Chicxulub impact
The most critical part of the Chicxulub controversy is the question of the age of this impact. Based on the original five KTB-identifying criteria (mass extinction, evolution of the first Danian species, KTB clay and red layer, iridium anomaly, and δ13C shift; Cowie, Ziegler & Remane, Reference Cowie, Ziegler and Remane1989; Keller, Li & MacLeod, Reference Keller, Li and MacLeod1995; Remane et al. Reference Remane, Keller, Hardenbol and Ben Haj Ali1999), no KTB age can be assigned to the Chicxulub impact. A strong belief that the Chicxulub impact is the cause for the KTB catastrophe led to the proposal to use the mass extinction and impact ejecta as the sole KTB markers (Molina et al. Reference Molina, Alegret, Arenillas, Arz, Gallala, Hardenbol, Von Salis, Steurbaut, Vandenberghe and Zaghbib-Turki2006). This re-definition was embraced by KTB impact proponents who then claimed that the sandstone complex in NE Mexico and Texas and the clastic complex of southern Mexico and Cuba are equivalent to the KTB red layer and iridium anomaly at the El Kef stratotype in Tunisia (e.g. Alegret et al. Reference Alegret, Arenillas, Arz, Diaz, Grajales-Nishimura, Melendez, Molina, Rojas and Soria2005; Arenillas et al. Reference Arenillas, Arz, Grajales-Nishimura, Murillo-Muneton, Alvarez, Camargo-Zanguera, Molina and Rosales-Dominguez2006; Schulte et al. Reference Schulte, Speijer, Mai and Kontny2006, Reference Schulte, Speijer, Brinkhuis, Kontny, Caleys, Galeotti and Smit2008, Reference Schulte, Alegret and Arenillas2010; see review in Keller, Reference Keller, Evans, Horton and King2008, Reference Keller, Keller and Adatte2011b ). Rather than solving the age problem, defining Chicxulub as KTB age has led to circular reasoning; one cannot evaluate the age of the Chicxulub impact by defining it.
A wealth of biostratigraphical data contradicts a KTB age for the Chicxulub impact. For example, in the Brazos sections the 0.30 m to 1.0 m burrowed claystone between the sandstone complex and the true KTB contains latest Maastrichtian palynomorphs and planktic foraminifera that indicate deposition prior to the KTB mass extinction (Prauss, Reference Prauss2009; Keller et al. Reference Keller, Adatte, Berner, Harting, Baum, Prauss, Tantawy and Stueben2007, Reference Keller, Adatte, Berner, Pardo and Lopez-Oliva2009, Reference Keller, Abramovich, Adatte, Berner, Keller and Adatte2011). Schulte et al. (Reference Schulte, Speijer, Mai and Kontny2006, Reference Schulte, Speijer, Brinkhuis, Kontny, Caleys, Galeotti and Smit2008) reported this interval (1.6 m) as barren and attributed it to the earliest Danian on the basis that the sandstone complex with impact spherules is by definition equivalent to the KTB at El Kef (see Keller, Reference Keller, Keller and Adatte2011b ). In addition, the lithified clasts with impact spherules at the base of the sandstone complex reveal a history of prior deposition, lithification, erosion, transport and redeposition all prior to the KTB (Keller et al. Reference Keller, Adatte, Berner, Harting, Baum, Prauss, Tantawy and Stueben2007; Adatte et al. Reference Adatte, Keller, Baum, Keller and Adatte2011).
A major problem for a Chicxulub KTB age is the up to 2 m thick spherule layer in undisturbed bedded marls of upper Maastrichtian zone CF1 about 4 m below the sandstone complex at El Penon, which can be traced over 90 m before the outcrop dips below the surface (Keller et al. Reference Keller, Adatte, Berner, Pardo and Lopez-Oliva2009). Schulte et al. (Reference Schulte, Alegret and Arenillas2010, p. 1215) interpret this and other spherule layers below the sandstone complex as ‘lenslike ejecta deposits in Mexico generated by impact-related liquefaction and slumping, consistent with the single very-high-energy Chucxulub impact.’ They do not explain the process by which spherules could have been injected as a layer at the same stratigraphic interval over 90 m and 4 m below their presumed origin, or the mechanism of liquefaction that could suck spherules 4 m into the underlying sediments and spread them out over 90 m, or why shallow water debris is absent and spherules are cemented by calcite or consist of welded melt rock glass, which is inconsistent with the debris-rich spherule layers at the base of the sandstone complex from which they are supposed to have derived. A more parsimonious explanation is primary deposition and pre-KTB age for the stratigraphically older layer, but erosion and downslope transport of the spherule layers in the sandstone complex at a much later time during the sea-level fall near the end of the Maastrichtian.
Similarly problematic for a Chicxulub KTB age interpretation is the 50 cm thick limestone between the impact breccia and the true KTB in the impact crater core Yaxcopoil-1. Keller et al. (Reference Keller, Adatte, Stinnesbeck, Rebolledo-Vieyra, Urrutia Fucugauchi, Kramar and Stueben2004a , b) demonstrated that this limestone is laminated, contains five glauconitic clay layers with burrows, three 1 cm layers of oblique bedding near the base (called cross bedding by Schulte et al. Reference Schulte, Alegret and Arenillas2010), late Maastrichtian (zone CF1) planktic foraminifera (see also Arz, Alegret & Arenillas, Reference Arz, Alegret and Arenillas2004), stable isotope signals, and palaeomagnetic C29R signals. All of these features indicate long-term deposition near the end of the Maastrichtian and well after the Chicxulub impact. Smit, van der Gaast & Lustenhouwer, (Reference Smit, van der Gaast and Lustenhouwer2004) interpreted this limestone as backwash and crater infill. Schulte et al. (Reference Schulte, Alegret and Arenillas2010) described the limestone layer as dolomitic sandstone, in part cross-bedded and containing ejecta clasts indicative of erosion and reworking. They ignored the glauconite layers, burrowing, stable isotopes, the latest Maastrichtian index species of zone CF1 and C29R palaeomagnetic data.
3. Methods
Cores across the KTB interval at ODP Sites 999B, 1001B, 1049A, 1049C, 1050C and 1259B were sampled at 1 cm, 5 cm and 10 cm intervals. Samples of about 1 cm3 were taken from central portions of the cores to avoid down-core contamination. For palaeontological analyses, samples were processed and washed through 36 μm and 63 μm screens following the procedure described in Keller et al. (Reference Keller, Li and MacLeod1995).
Quantitative palaeontological analysis was performed on aliquots of 300–400 specimens in the > 63 μm size fraction. The 36–63 μm size fraction was examined for small species and counts were done for intervals where Danian species are more abundant in this size fraction. All specimens were identified and mounted on microslides for a permanent record. Preservation of foraminifera at Blake Nose ODP Sites 1049–1050 is excellent, though some calcite overgrowth and dissolution effects are apparent in sugary textures in the spherule layer. At Bass River and Site 1259B moderate dissolution has resulted in holes in the chambers of some Danian species. Strong dissolution and poor preservation at Sites 999B and 1001B prevented quantitative analysis and the species census data was obtained from both washed residues and thin sections.
Stable isotope analysis of samples from ODP Sites 1049A, 1049C, 1050C and 1259B was performed on the bulk fine fraction (<63 μm). The abundance of reworked Cretaceous species in the early Danian of Sites 1049A and 1049C necessitated analysis of the Danian planktic species Parvularugoglobigerina eugubina and Chiloguembelina morsei to obtain a reliable Danian record. Analyses were performed at the Institute of Mineralogy and Geochemistry, Lausanne University, using a Thermo Fisher Scientific (Bremen, Germany) GasBench II preparation device interfaced with a Thermo Fisher Scientific Delta Plus XL continuous flow isotope ratio mass spectrometer (IRMS) (Revesz, Landwehr & Keybl, Reference Revesz, Landwehr and Keybl2001). Stable carbon and oxygen isotope ratios are reported in the delta notation as the per mil (‰) deviation relative to the Vienna Pee Dee belemnite standard (VPDB). Analytical uncertainty (2σ) monitored by replicate analyses of the international calcite standard NBS-19 and the laboratory standards Carrara Marble and Binn Dolomite is no greater than ±0.05‰ for δ13C and ±0.1‰ for δ18O.
4. KT boundary definition and placement
The controversy over the age of the Chicxulub impact is in no small measure the result of redefining the KTB as the impact ejecta layer based on the assumption that this impact caused the end-Cretaceous mass extinction. Unfortunately, this has led to confusion and circular reasoning. For example, Gradstein, Ogg & Smith, (Reference Gradstein, Ogg and Smith2004) proposed to define the KTB based on just two criteria ‘the Ir anomaly associated with a major mass extinction.’ This has proved impractical because small (~ 1 ppb) Iridium anomalies may be remobilized and concentrated at redox boundaries thus leading to erroneous KTB placement (e.g. Graup & Spettel, Reference Graup and Spettel1989; Kramar et al. Reference Kramar, Stueben, Berner, Stinnesbeck, Philipp and Keller2001; Colodner et al. Reference Colodner, Boyle, Edmond and Thomson1992; Lee, Wasserburg & Kyte, Reference Lee, Wasserburg and Kyte2003; Miller et al. Reference Miller, Sherrell, Browning, Field, Gallagher, Olsson, Sugarman, Tuorto and Wahyudi2010; Gertsch et al. Reference Gertsch, Keller, Adatte, Berner, Keller and Adatte2011). Partly as a solution, Molina et al. (Reference Molina, Alegret, Arenillas, Arz, Gallala, Hardenbol, Von Salis, Steurbaut, Vandenberghe and Zaghbib-Turki2006) proposed that any impact ejecta define the KTB (e.g. spherules, breccia, iridium anomaly, shocked quartz, spinels etc.). However, defining any presumed impact related material as KTB age leads to circular reasoning and rules out evaluating the true age of the Chicxulub impact. This has led to placing the KTB based on the select occurrence of impact spherules at the base of clastic deposits, at the base of breccia deposits of dubious origin in the Maastrichtian, and at spherule layers reworked into Danian sediments (e.g. Smit, Reference Smit1999; Norris, Huber & Self-Trail, Reference Norris, Huber and Self-Trail1999; Arenillas et al. Reference Arenillas, Arz, Grajales-Nishimura, Murillo-Muneton, Alvarez, Camargo-Zanguera, Molina and Rosales-Dominguez2006; Schulte et al. Reference Schulte, Speijer, Mai and Kontny2006, Reference Schulte, Alegret and Arenillas2010).
These pitfalls can be avoided by returning to the unique defining and supporting KTB criteria of the GSSP, which consistently provide the most reliable KTB identification globally (Keller, Reference Keller, Keller and Adatte2011b ). The KTB is globally defined based on five criteria developed at the El Kef stratotype section and point (GSSP) and Elles parastratotype section in Tunisia (Cowie, Zielgler & Remane, Reference Cowie, Ziegler and Remane1989; Keller, Li & MacLeod, Reference Huber, MacLeod and Turl995; Keller et al. Reference Keller, Adatte, Stinnesbeck, Luciani, Karoui and Zaghbib-Turki2002; Remane et al. Reference Remane, Keller, Hardenbol and Ben Haj Ali1999). Two of the five are unique defining criteria: the mass extinction of 66% of planktic foraminiferal taxa and a few centimetres above the extinction horizon the first appearance (FA) of five Danian species (Parvularugoglobigerina extensa, Woodringina hornerstownensis, W. claytonensis, Globoconusa daubjergensis, Eoglobigerina eobulloides). Three are supporting criteria, which are not unique in the geological record but characterize the KTB: the clay layer, iridium anomaly and negative shift in δ13C values. All five criteria have been tested and successfully applied to hundreds of KTB sequences worldwide and their presence is powerful confirmation of environmental changes across this boundary event. These KTB criteria are applied in this study. However, the iridium anomaly in all sections is small, located in early Danian sediments and apparently represents reworking and remobilization of iridium (Kramar et al. Reference Kramar, Stueben, Berner, Stinnesbeck, Philipp and Keller2001; Gertsch et al. Reference Gertsch, Keller, Adatte, Berner, Keller and Adatte2011).
5. Lithology and biostratigraphy
Age control is based on high-resolution quantitative planktic foraminiferal biostratigraphy, which provides the best age resolution for palaeomagnetic chron C29r spanning the KTB (Fig. 2). C29r below the KTB corresponds to planktic foraminiferal zones CF1 and CF2. The calcareous nannofossil zone Micula prinsii spans from the top of CF3 to the KTB. Above the KTB, C29r corresponds to Danian zones P0, P1a with subzones P1a(1) and P1a(2), which is encompassed by the lower part of nannofossil zone CP1a (or NP1a). Previous studies of the sections analyzed in this report have used standard nannofossil and/or planktic foraminiferal biostratigraphy without quantitative data, which is not well suited for evaluating the completeness of the KTB transition.
5.a. Bass River core, New Jersey
At Bass River, New Jersey (39°36′42″N, 74°26′12″W) the KTB transition was recovered in a relatively shallow water environment (palaeodepth ~100 m; Miller et al. Reference Miller, Sugarman, Browning, Olsson, Pekar, Reilly, Cramer, Aubry, Lawrence, Curran, Stewart, Metzger, Uptegrove, Bukry, Burckle, Wright, Feigenson, Brenner and Dalton1998). Sediments consist of calcareous glauconitic silty shale with bivalve shells that disconformably underlie a 6 cm thick glauconitic impact spherule layer with clasts (1–6 cm, Fig. 3). A mm thin red-brown limonite layer separates the spherule layer from the overlying 6 cm thick glauconitic silty-sandy shale, which is interrupted by an 8–10 cm thick dense layer of phosphate clasts and Steinkerns with an erosion surface at the base. About 10 cm above this layer is a 5 cm thick burrowed interval truncated by a 16 cm core gap. The glauconitic silty shale above the core gap contains occasional phosphate clasts.
Planktic foraminifera from the sediments below the 6 cm thick spherule layer are relatively well preserved, though small, thin walled and fragile. A typical low diversity middle shelf planktic foraminiferal assemblage marks zone CF3 as indicated by the presence of Gansserina gansseri, the index species that marks the CF3–CF2 boundary (Fig. 3). This indicates erosion of the CF2–CF1 interval that spans the top 450 ka of the Maastrichtian based on the time scale of Cande & Kent (Reference Cande and Kent1995), or 280 ka based on Gradstein, Ogg & Smith, (Reference Gradstein, Ogg and Smith2004) (Fig. 2). Olsson et al. (Reference Olsson, Miller, Browning, Habib and Sugarmann1997) reported Micula prinsii in the 8 cm below the spherule layer, which is consistent with the overlap of this species in the upper part of zone CF3.
The spherule layer contains rare Cretaceous planktic foraminifera (e.g. Planoglobulina brazoensis, P. carseyae, Rugoglobigerina rugosa, R. macrocephala; Fig. 3). The 6 cm immediately above the spherule layer contains common small Cretaceous species (KT survivors or reworked) and a diverse (11 species) early Danian assemblage including Parvularugoglobigerina eugubina, Parasubbotina pseudobulloides and Subbotina triloculinoides, which is characteristic of subzone P1a(2) (upper part of P1a zone). Thus the total interval missing at the KTB unconformity includes zones CF1, CF2, P0, P1a(1) and most of P1a(2), and spans about 620 ka; Fig. 2). Just 6 cm above the spherule layer this Danian P1a(2) assemblage and Cretaceous species terminate at the phosphate clast layer, which marks another hiatus between P1a(2) and P1b (Figs 2, 3). Above this hiatus the zone P1b assemblage is dominated by Chiloguembelina morsei, Praemurica taurica and Parasubbotina pseudobulloides.
5.2. North Atlantic ODP Sites 1049A, 1049C and 1050C
Ocean Drilling Program (ODP) Leg 171 drilled three holes at Site 1049 (A, B, C) on the eastern margin of Blake Nose (30°08.5436′N, 76°06.7312′W) at a water depth of 2671 m (Norris et al. Reference Norris, Kroon, Klaus, Baez and Kapitan-White1998; Norris, Huber & Self-Trail, 1999). In each hole, located just 10 m apart from the next, the KTB transition was recovered along with a 10–15 cm thick impact spherule-rich layer. Site 1050C was drilled at a water depth of 2300 m, 330 m upslope from Site 1049 and at a distance of 10 km (30°05.9953′N, 76°14.0997′W; Fig. 1). The KTB transition was also recovered but no impact spherules. We analyzed two holes at Site 1049 (A, C) and Site 1050C in order to insure replicable results and exclude artifacts of sampling or processing. All three sites analyzed show very similar lithologies, except for the absence of a spherule layer at Site 1050C (Fig. 4).
Upper Maastrichtian sediments at Sites 1049A and 1049C consist of disturbed, bioturbated, mottled light and darker grey nannofossil and foraminiferal ooze, which underlie the 10–15 cm thick spherule layer (Fig. 4). The spherule layer has a carbonate ooze matrix, irregular bedding, clasts of limestone, chalk and dolomite and is bounded by irregular surfaces. At Site 1049A a 2 mm thin red-brown limonite layer and 8 cm thick white ooze overlie the spherule layer followed by grey carbonate ooze. At Site 1049C, the limonite layer is also present but discontinuous and white ooze is present in the middle of the spherule layer. Above the limonite layer the 8 cm dark and light grey ooze is strongly mottled and bioturbated, followed by light grey ooze. At Site 1050C the KTB transition consists of dark and light grey burrowed ooze with some large horizontal burrows infilled with grey ooze that contains rare Danian species below the KTB (Fig. 4). An erosional surface and darker grey layer is correlative with the spherule layer.
Upper Maastrichtian faunal assemblages at all three sites are diverse, abundant and well preserved but lack Plummerita hantkeninoides, the zone CF1 index species for the topmost Maastrichtian (Figs 5–7). However, the presence of Gansserina gansseri (and G. wiedenmayeri in Site 1050C) and Pseudoguembelina hariaensis mark this interval as zone CF3 (Abramovich et al. Reference Abramovich, Yovel‐Corem, Almogi‐Labin and Benjamini2010). Micula prinsii was reported from a short interval below the spherule layer (Norris, Huber & Self-Trail, Reference Norris, Huber and Self-Trail1999), which is consistent with the upper part of zone CF3 (Fig. 2). This indicates a major hiatus that spans Maastrichtian zones CF1 and CF2 (~ 280 ka, time scale of Gradstein, Ogg & Smith, Reference Gradstein, Ogg and Smith2004; Fig. 2).
5.b.1. Impact spherule layer with early Danian species
The spherule-rich layer above the KTB unconformity at Sites 1049A and 1049C is 15 cm and 10 cm thick, respectively (Fig. 4). Reworked Cretaceous species are common (Figs 5, 6) along with shallow and deep-water benthic foraminifera that indicate reworking and transport of shallow water sediments into the basin (Alegret & Thomas, Reference Alegret and Thomas2004). We also observed a dominance of large and thick-shelled species with variable preservation, including discoloration, abrasion and broken shells, and low abundance of small species, all of which indicate reworking and transport. On a quantitative basis 83–95% (>63 μm) of the Maastrichtian species are reworked in the spherule layer and 5–17% are small Danian species.
The relative abundance of Danian species was calculated based on the 63–100 μm size fraction from the spherule layer (Figs 5, 6). The presence of diverse assemblages (12 species), including P. eugubina and P. longiapertura, but absence of P. pseudobulloides and/or S. triloculinoides indicates a subzone P1a(1) assemblage in the spherule layer with the lower part of this subzone and zone P0 missing. A similar Danian assemblage is present in the 8 cm above the spherule layer at Site 1049A and in the 20 cm above the spherule layer at Site 1049C, whereas the abundance of reworked Cretaceous species rapidly decreases to <10% (Figs 5, 6). The first appearances P. pseudobulloides, S. triloculinoides and the sharp increase in Chiloguembelina morsei and Guembelitria species 8 cm and 20 cm above the spherule layer at Sites 1049A and 1049C marks a hiatus between P1a(1) and P1a(2). At both sites P1a(2) is dominated by two opportunistic species, C. morsei and G. cretacea, that indicate increased environmental stress. Reworked Cretaceous species are still present but rare (<2%).
At Site 1050C the comparable P1a(1) and P1a(2) assemblages spans the 60 cm above the KTB unconformity. Reworked species are very rare (Fig. 7). The basal Danian sample is dominated by P. eugubina (50%), which indicates that zone P0 and the lower part of P1a(1) is missing along with zones CF1 and CF2 (~400 ka). Subzone P1a(1) is dominated by P. eugubina (55–64%) with common Eoglobigerina edita and variable presence of seven other Danian species (Fig. 7). The P1a(1)–P1a(2) boundary is marked by a sharp decrease in P. eugubina, concurrent increase in Guembelitria spp. and Chiloguembelina morsei, first appearances of P. pseudobulloides and P. taurica. For the next 30 cm P. eugubina and E. edita dominate along with C. morsey. At the top of this interval E. edita becomes rare and P. eugubina disappears, which marks the P1a–P1b boundary. At this interval the sudden disappearance of P. eugubina, decrease in E. edita, sudden increase (to 35–40%) in Guembelitria spp. and coincident drop in δ13C values by nearly 1‰, suggests another hiatus. During zone P1b Guembelitria spp. and C. morsei dominate with little change.
Comparison of Sites 1049A, 1049C and 1050C reveals both similarities in faunal assemblages and differences in relative abundances of species populations. The latter are primarily due to variable erosion at the KTB and P1a(1)–P1a(2) hiatuses. In zone P1b all sections show the same faunal assemblages dominated by Guembelitria sp. and Chiloguembelina morsei. Carbon isotope signals also reflect great similarities among these sequences and each hiatus is marked by abrupt shifts in δ13C values (see Section 6).
5.c. Caribbean ODP Sites 999B and 1001B
ODP Sites 1001A and B were drilled 30 m apart on the Hess Escarpment (lower Nicaragua Rise, 15°45.427′N, 74°54.627′W) at depth of 3259 m (Fig. 1). The KTB transition was recovered at Site 1001B, core 18R-5 (Sigurdsson et al. Reference Sigurdsson, Leckie, Acton, Miller and Maddox1997). Maastrichtian sediments consist of light grey mottled, burrowed limestone, which in the shipboard photo is speckled with dark flakes from the overlying claystone, but not in the archive core photographed in 2006 (Fig. 8). This limestone contains common Heterohelix globulosa and Ruguglobigerina rugosa but all other species are rare. This low diversity assemblage is partly due to poor preservation and thin-section identification. However, the rarity of globotruncanids, which are more dissolution and break-resistant, and absence of upper Maastrichtian index species indicates a lower upper Maastrichtian zone CF5 age in agreement with the age previously assigned as base Abathomphalus mayaroensis (Sigurdsson et al. Reference Sigurdsson, Leckie, Acton, Miller and Maddox1997). This indicates a major KTB hiatus with all or most of the upper Maastrichtian missing (zones CF1–CF4; Fig. 2).
Above the unconformity is a 1.5–2.0 cm thick green-grey claystone with impact spherules altered to smectite and containing a lower Danian subzone P1a(1) assemblage (Fig. 8). The overlying 1.5–2.0 cm thick dark green-grey claystone contains abundant clay-altered spherule debris but no foraminifera. No sample was available for this study from the overlying 1.5 cm thick brown-grey claystone. Above this claystone at Site 1001B is a 2 cm thick white-grey chalk, which is not observed at Hole 1001A but well developed (~10 cm) at Site 999B, where it overlies a white limestone (Fig. 8). The 3–5 cm thick overlying claystone at holes 1001A and B contains a more diverse subzone P1a(1) assemblage of nine species (Fig. 8).
Sigurdsson et al. (Reference Sigurdsson, Leckie, Acton, Miller and Maddox1997) reported the first appearance of lower Danian species in this white-grey chalk of Hole 1001B. In the overlying claystone the abrupt disappearance of the early Danian species (P. eugubina, E. edita, W. claytonensis and W. hornerstownensis) and simultaneous first appearance of S. triloculinoides. P. pseudobulloides and S. trivialis (core 18R-5, 14 cm) marks a hiatus between P1a(1) and P1b with subzone P1a(2) missing. This data suggests a similar KTB and early Danian erosion pattern to Bass River and Blake Nose with impact spherules reworked into early Danian subzone P1a(1). However, erosion of Maastrichtian sediments at the KTB unconformity is more extensive spanning zones CF1–CF4 (~3 Ma; Figs 2, 8).
ODP Site 999B was drilled on Kogi Rise, Colombian Basin (12°44.639′N, 78°44.360′W) at 2828 m depth (Fig. 1). Recovery across the KTB was incomplete. Sigurdsson et al. (Reference Sigurdsson, Leckie, Acton, Miller and Maddox1997) assembled a composite section based on three segments (59R-3, 59R-CC and the top of core 60R-1) and placed the KTB at the base of a white limestone (60R-1, 1 cm), correlative with the white-grey chalk at Site 1001B (Fig. 8). No spherules are present and it was suggested they may have been lost in the drilling process.
Samples were obtained from the KTB transition. Preservation of foraminifera is poor due to dissolution, and the biostratigraphy was determined from species in washed residues and from thin sections. The Maastrichtian limestone at Site 999B (up to 60R-1, 15–18 cm) contains the same impoverished assemblage indicative of zone CF5 as at Site 1001B, which is indicative of the lower upper Maastrichtian (Fig. 2). Sigurdsson et al. (Reference Sigurdsson, Leckie, Acton, Miller and Maddox1997, p. 154) suggested an upper Maastrichtian age of CC26 (Nephrolithus frequens) or underlying CC25 zone (Arkhangelskiella cymbiformis). They concluded: ‘it is not possible to determine whether the K/T boundary at Site 999 is complete from a nannofossil biostratigraphic viewpoint.’
Between the grey Maastrichtian limestone and overlying white limestone is a 1–2 cm thick grey claystone with rare Maastrichtian species but abundant Guembelitria cretacea and few W. hornerstownensis, G. daubjergensis and P. eugubina, which indicate deposition in the lower Danian subzone P1a(1). The overlying 10 cm thick white limestone also contains a lower Danian subzone P1a(1) assemblage, similar to the thin grey-white chalk at Site 1001B (Fig. 8). Based on this data we place the KTB unconformity between the Maastrichtian limestone and grey claystone at the first Danian assemblage.
The white-grey limestone (59R-CC, 1–15 cm) above the white limestone contains rare and poorly preserved Danian species. Above this layer is a darker grey clay-rich marl that contains a diverse lower Danian subzone P1a(2) assemblage (G. cretacea, Eoglobigerina edita, P. eugubina, Woodringina hornerstownensis, W. claytonensis, Globoconusa daubjergensis, Praemurica taurica, S. triloculinoides, P. pseudobulloides, G. compressa), which is missing at Site 1001B (Fig. 8). However, the overlying zone P1b assemblage (59R-3, 74–82 cm) is present in both sites. These data indicate a more expanded sediment record across the KTB at Site 999B. However, because this record was pieced together from three different core segments the true sediment record remains uncertain.
5.4. Demerara Rise ODP Site 1259B
ODP Site 1259 was drilled on Demerara Rise (9°18.048′N, 54°11.945′W) at a water depth of 2354 m (Fig. 1). A 2 cm thick impact spherule layer was recovered at Sites 1259B and 1259C (as well as Sites 1258A, 1260A, Shipboard Scientific Party, Reference Erbacher, Mosher and Malone2004; MacLeod et al. Reference MacLeod, Whitney, Huber and Koeberl2007; Schulte et al. Reference Schulte, Deutsch, Salge, Berndt, Kontny, MacLeod, Neuser and Krumm2009). Most KTB studies have focused on Sites 1259B and 1259C. We obtained samples from Site 1259B where the spherule layer is in core 13-1, 47–49 cm; 31 samples were analyzed for this report. Sediments below the spherule layer consist of mottled, bioturbated marly chalk with a relatively sharp upper contact. The spherule layer is graded and thin white chalk layers are present at the base and top, similar to Site 1001B. Danian sediments consist of strongly mottled, bioturbated light and darker grey clayey chalk.
Preservation of planktic foraminifera is good for the Danian and moderate in the Maastrichtian with dissolution effects (thin walls, holes in test chambers). The marly chalk below the spherule layer contains a diverse late Maastrichtian planktic foraminiferal assemblage, including rare Plummerita hantkeninoides, the index species for zone CF1 (Fig. 9). Cretaceous assemblages in the >63 μm size fraction are dominated by small biserial species (Pseudoguembelina costulata, Heterohelix planata, H. navarroensis), whereas Guembelitria cretacea is common in the smaller size fraction (36–63 mm). Rare Danian species were obtained from the 10 cm thick bioturbated light and dark grey clay below the KTB unconformity and spherule layer; these specimens are likely due to burrowing from the overlying Danian.
One sample was obtained from the spherule layer. Upper Maastrichtian planktic foraminifera are common in the >63 μm size fraction, and small Danian species are common in the 36–63 μm size fraction where Guembelitria species dominate (83–96%) in the spherule layer and up to 8 cm above it (Fig. 9). The presence of P. eugubina, G. daubjergensis, W. hornerstownensis, W. claytonensis and E. eobulloides in the spherule layer and the 8 cm above indicate subzone P1a(1). The first P. pseudobulloides and Subbotina triloculinoides appear in sample 13-1, at 37–39 cm, and mark the P1a(1)–P1a(2) subzone boundary just 8 cm above the spherule layer. A short hiatus is present at this level as indicated by the abrupt decrease in Guembelitria spp. and correlative abrupt recovery in δ13C values. Above this hiatus the rapid increase in Globoconusa daubjergensis followed by C. morsei and the extinction of P. eugubina 37 cm above the spherule layer marks the P1a(2)–P1b boundary. A short hiatus may also be present at this level as indicated by the abrupt assemblage change to dominant Guembelitria spp. (85–90%) and abrupt decrease in C. morsei.
These data suggest a relatively complete latest Maastrichtian (zone CF1 index species present) but a short KTB hiatus with most of the Danian subzone P1a(1) and underlying zone P0 missing. In addition, a short hiatuses is present at the P1a(1)–P1a(2) boundary and probably P1a(2)–P1b boundary. Similar early Danian erosion and redeposition patterns are observed at Sites 999B, 1001B, 1049A, 1049C, 1050C and Bass River, although the extent of erosion varies at each location. This is indicated by the variable species abundance patterns between sites although in the same biozones P1a(1), P1a(2) and P1b and presence of similar δ13C signals. For example, at Site 1259B P. eugubina is relatively rare, subzone P1a(1) is only 8 cm thick, dominated by Guembelitria spp., but with a diverse (14 species) assemblage that indicates an upper P1a(1) interval not observed in the other sections. Subzone P1a(2) is dominated by C. morsei, similar to all other sites, but fewer Guembelitria spp. Zone P1b differs in its great abundance of Guembelitria (~90%) and low abundance of all other species. The difference here is likely due to the core gap of 1.33 m where a critical part of the faunal transition is missing (Fig. 9).
6. Stable isotope stratigraphy
One of the characteristic KTB markers is a 2–3‰ negative δ13C shift at the base of the Danian, which reflects the productivity crash after the mass extinction. Similarly, the δ13C variations in the early Danian mark the delayed marine recovery and setbacks in the aftermath of the mass extinction. Therefore, carbon isotope stratigraphy, coupled with planktic foraminiferal biostratigraphy, is an excellent and independent method for determining the completeness of the sediment record. No such high-resolution data is published to date for any of the sites analyzed for this study. For Blake Nose (Sites 1049A, 1049B and 1050) and Demerara Site 1259B we analyzed the bulk rock fine fraction (36–63 μm), which consists primarily of nannoplankton and small planktic foraminifera and is generally a good indicator of primary productivity (Zachos, Arthur & Dean, Reference Zachos, Arthur and Dean1989; Schrag, DePaolo & Richter, Reference Schrag, DePaolo and Richter1995). For the Maastrichtian, δ13C values vary between 1.8 and 2.2‰, except for Site 1050C which is slightly heavier (2.2–2.4‰; Fig. 10). These values are comparable to those observed in open marine environments, such as South Atlantic Site 525A (Li & Keller, Reference Li and Keller1998a ). No significant changes in δ13C are apparent for the intervals of CF3 and CF1 analyzed below the KTB unconformity.
At Sites 1050C and 1259B, the KTB unconformity is marked by a sudden drop in δ13C values of 1.2‰ and 2.0‰, respectively. At Sites 1049A and 1049C the bulk fine fraction yielded a gradual decrease in δ13C values well into the early Danian because of the abundance of reworked Cretaceous species (Figs 5, 6). We obtained an uncompromised lower Danian record from monospecific samples of the Danian species P. eugubina and C. morsei. In some samples both species were analyzed and the interspecies variation is less than 0.3‰. Preservation of Danian foraminiferal tests is very good and shell calcite is not recrystallized, but there is minor calcite overgrowth on shell surfaces, which may account for some of the interspecies variation. About 200–300 specimens were picked per sample for stable isotope analysis. In the lower and middle part of the spherule layer only about 100–150 specimens per species could be obtained in the available sample (in addition to the quantitative counts that are archived) and this quantity turned out to be insufficient for reliable isotope signals. Nevertheless, in the lower Danian subzone P1a(1), including the top of the spherule layer, the δ13C values of Danian species reached lows of −0.5‰ suggesting a drop of about 2‰ across the KTB, similar to Site 1259B (Fig. 10). These values are characteristic for the lower Danian above the impact spherule layer in Haiti and the equivalent P1a(1) interval in Tunisia (Keller & Lindinger, Reference Keller and Lindinger1989; Keller et al. Reference Keller, Adatte, Stinnesbeck, Stueben and Berner2001; Stueben et al. Reference Stueben, Kramar, Berner, Eckhardt, Stinnesbeck, Keller, Adatte and Heide2002).
Major variations in δ13C signals through the early Danian can be correlated despite the frequent hiatuses and variable erosion identified from planktic foraminiferal biostratigraphy and species abundance patterns. The lowest δ13C values are known from the early Danian in the aftermath of the mass extinction. In all sites analyzed, at least part of this interval is recorded in fine fraction and planktic foraminifera of subzone P1a(1) that includes the spherule layer (Fig. 10). An abrupt positive shift of about 1‰ coincides with the P1a(1)–P1a(2) hiatus. In P1a(2), δ13C signals remain relatively high but do not reach pre-KTB values. The P1a(2) record is fragmented in the sites analyzed due to variable erosion similar to P1a(1). At Site 1050C a negative δ13C shift at the P1a(2)–P1b boundary marks the hiatus. This shift was not detected at Sites 1049A and 1259B because of the 0.80 m and 1.33 m core gaps, respectively.
7. Discussion
7.a. KTB unconformity
A thin impact spherule layer between Maastrichtian and Paleocene sediments in some KTB sequences from the North Atlantic and Caribbean, is commonly cited as proof that the Chicxulub impact is KTB in age and the sedimentation record complete (e.g. Olsson et al. Reference Olsson, Miller, Browning, Habib and Sugarmann1997; Sigurdsson et al. Reference Sigurdsson, Leckie, Acton, Miller and Maddox1997; Norris et al. Reference Norris, Kroon, Klaus, Baez and Kapitan-White1998; Norris, Huber & Self-Trail, Reference Norris, Huber and Self-Trail1999; Martinez-Ruiz et al. Reference Martinez-Ruiz, Ortega-Huertas, Palomo-Delgado, Smit, Kroon, Norris and Klaus2001; MacLeod et al. Reference MacLeod, Whitney, Huber and Koeberl2007; Schulte et al. Reference Schulte, Alegret and Arenillas2010). Results based on high-resolution biostratigraphy, quantitative faunal analyses and stable isotopes show that far from being complete all seven sections analyzed show a major unconformity across the KTB with variable erosion of Maastrichtian sediments.
The maximum extent of erosion can be estimated based on Paleocene and Maastrichtian sediments above and below the unconformity, respectively. At Bass River, New Jersey, the KTB unconformity spans from the early Danian subzone P1a(2) through to the upper Maastrichtian zones CF1 and CF2 (~620 ka; Figs 2, 3), in North Atlantic Sites 1049A, 1049C and 1050C from P1a(1) through to CF1 and CF2 (~400 ka; Figs 5–7), and in the Caribbean Sites 999B and 1001B from P1a(1) through to zones CF1, CF2, CF3 and CF4 (~3 Ma; Fig. 8). Only at Demerara Rise Site 1259B is erosion or non-deposition more limited, spanning most of subzone P1a(1) and zone P0 (~150–200 ka; Fig. 9). The nearby Site 1258A provides some clues to possible erosion at the top of the Maastrichtian. Just below the KTB, Thibault & Gardin (Reference Thibault and Gardin2006) recorded the abrupt end of the nannofossils Micula murus acme that marks the CF1 warm event. In complete sequences, this global warm event ends with the onset of cooling about 50–100 ka prior to the KTB (Li & Keller, Reference Li and Keller1998a ; Abramovich & Keller, Reference Abramovich and Keller2003; Abramovich et al. Reference Abramovich, Keller, Berner, Cymbalista, Rak, Keller and Adatte2011), which suggests erosion of the latest Maastrichtian at Site 1258A and likely at Site 1259B.
We can estimate the missing interval above the KTB based on Danian sediments present within and above the spherule layer. In a complete sequence, biozones P0 + P1a(1) + P1a(2) span C29R above the KTB, or 380 ka based on Gradstein, Ogg & Smith, (Reference Gradstein, Ogg and Smith2004), and half of this is P0 + P1a(1) (Fig. 2). In all sections examined, zone P0 is absent and P1a(1) is very reduced (missing at Bass River) with abrupt changes in species abundances and the δ13C shift reflecting a hiatus (Figs 5–9). Similarly, subzone P1a(2) is reduced with abrupt changes in δ13C and species abundances marking a hiatus at the P1a(2)–P1b boundary (Figs 2–9).
Early Danian P1a hiatuses have been observed in many sections worldwide (e.g. MacLeod & Keller, Reference MacLeod and Keller1991a , Reference MacLeod and Keller b ; Keller et al. Reference Keller, Adatte, Stinnesbeck, Luciani, Karoui and Zaghbib-Turki2002, Reference Keller, Stinnesbeck, Adatte, Holland, Stueben, Harting, De Leon and De La Cruz2003a , b; Keller, Reference Keller, Evans, Horton and King2008) and may reflect global climate and sea level changes (Adatte, Keller & Stinnesbeck, Reference Adatte, Keller and Stinnesbeck2002). However, North Atlantic and Caribbean sites differ in the extensive erosion and/or non-deposition in the early Danian (P0–P1a–P1b) and spanning well into the upper Maastrichtian (CF1–CF2; Fig. 2). We attribute this primarily to tectonic activity in the Caribbean, climate and sea level changes, and intensified current circulation (e.g. Gulf Stream, Keller et al. Reference Keller, Lyons, MacLeod and Officer1993; Watkins & Self-Trail, Reference Watkins and Self-Trail2005). The Chicxulub impact may have contributed mainly via mass wasting of Maastrichtian sediments as suggested by Norris et al. (Reference Norris, Firth, Blusztajn and Ravizza2000) and Klaus et al. (Reference Klaus, Norris, Kroon and Smit2000), although these sections do not support this scenario.
7.b. Age and origin of impact spherules
A KTB age for the Chicxulub impact was assumed based on spherules between Maastrichtian and Danian sediments at the sites investigated (Fig. 11; e.g. Olsson et al. Reference Olsson, Miller, Browning, Habib and Sugarmann1997; Sigurdsson et al. Reference Sigurdsson, Leckie, Acton, Miller and Maddox1997; Norris, Huber & Self-Trail, Reference Norris, Huber and Self-Trail1999; MacLeod et al. Reference MacLeod, Whitney, Huber and Koeberl2007). There are three key factors that argue against this interpretation: (1) there is a major KTB unconformity in all sequences, except for Site 1259B; (2) the spherule-rich sediments contain common reworked Maastrichtian species and shallow water benthic foraminifera (Alegret & Thomas, Reference Alegret and Thomas2004), which argues for erosion and redeposition, and (3) the spherule layer contains a diverse early Danian P1a(1) species assemblage, which argues for deposition during the early Danian, well after the mass extinction and evolution of the first Danian species. Earlier studies reported the first Danian species above the spherule layer and thereby justified placing the KTB at its base. It is possible that the early Danian assemblages were not observed because of their small size.
One might argue that the Danian species in the spherule layer are an artifact of bioturbation or leaking of sediments into the more porous spherule layer. Some bioturbation is clearly evident in the lithology of Sites 1049A, 1049C, 1050C (Fig. 4). Similarly at Site 1259B, rare Danian species were observed in burrows below the KTB unconformity (Fig. 9), contrary to Schulte et al. (Reference Schulte, Deutsch, Salge, Berndt, Kontny, MacLeod, Neuser and Krumm2009, p. 1185) who noted that ‘burrows or other signs of bioturbation are totally lacking.’ Although bioturbation may account for some Danian specimens in the spherule layer, this could not account for the diverse assemblage with hundreds of specimens and similar species abundances and diversity as in assemblages above the spherule layer.
The presence of these Danian assemblages in the spherule layer above the unconformity indicate that sedimentation resumed in the upper subzone P1a(1) about 150–200 ka after the KTB mass extinction, except for Bass River where the age is P1a(2). This indicates that the spherules are eroded, transported and redeposited long after the Chicxulub impact. Because the normal early Danian evolutionary pattern is present, except for the missing basal Danian, we can assume that the Danian fauna represents in situ deposition but with added downslope transport of spherule-rich reworked debris and Maastrichtian foraminifera. The absence of spherules in Site 1050C can be explained by its distant location from the submarine channel that transported debris to Sites 1049. Similar reworking of shallow water sediments and spherules, transport and redeposition into deeper waters is also observed in NE Mexico sections and Site 1259 (Alegret, Molina & Thomas, 2001; Alegret et al. Reference Alegret, Arenillas, Arz, Liesa, Melendez, Molina, Soria and Thomas2002; Alegret & Thomas, Reference Alegret and Thomas2004). Based on the results from this study, the spherules in lower Danian sediments must be reworked from an older primary deposit. Reworked impact spherules have previously been observed in Danian sediments from Haiti, Cuba, Belize, Guatemala and southern Mexico (Keller et al. Reference Keller, Adatte, Stinnesbeck, Stueben and Berner2001, Reference Keller, Stinnesbeck, Adatte, Holland, Stueben, Harting, De Leon and De La Cruz2003a ; Alegret et al. Reference Alegret, Arenillas, Arz, Diaz, Grajales-Nishimura, Melendez, Molina, Rojas and Soria2005; Arenillas et al. Reference Arenillas, Arz, Grajales-Nishimura, Murillo-Muneton, Alvarez, Camargo-Zanguera, Molina and Rosales-Dominguez2006; Keller, Reference Keller, Evans, Horton and King2008).
The strongest case for primary spherule deposition precisely at the KTB, and hence precise KTB age for the Chicxulub impact, was made for Demerara Site 1259B by Schulte et al. (Reference Schulte, Deutsch, Salge, Berndt, Kontny, MacLeod, Neuser and Krumm2009) based on geochemical and mineralogical analyses. Based on their results (p. 1185) they divided the 1.5–2 cm thick spherule layer into two distinct parts: the upper 0.5–0.7 mm with (1) ‘locally enhanced Si, Al, and K contents due to the presence of silicic detritus (i.e. quartz, feldspars, mica), and (2) local Ca- and Mg-enrichment, indicating the occurrence of calcite and dolomite clasts.’ Moreover, they observed a strong negative Nd anomaly, which indicates uptake from contemporaneous ocean water during alteration, which in turn implies a long term residence time for the spherules in seawater. Schulte et al. (Reference Schulte, Deutsch, Salge, Berndt, Kontny, MacLeod, Neuser and Krumm2009, p. 1200) concluded that the Nd isotope systematics ‘is best explained by uptake of neodymium from the ocean during the alteration process; this Nd directly reflects fluviatile input from the rather old crustal material in the South American hinterland.’ Rather than primary deposition at the time of the Chicxulub impact as argued by Schulte et al. (Reference Schulte, Deutsch, Salge, Berndt, Kontny, MacLeod, Neuser and Krumm2009), these data indicate redeposition of spherules and other detritus that originally accumulated somewhere close to the Guyana craton. This interpretation is consistent with the biostratigraphic results of this study.
7.c. KTB unconformity and the Chicxulub impact
One might argue that the KTB unconformity is unusual and must have an extraordinary cause, such as the Chicxulub impact (e.g. Bralower, Paull & Leckie, Reference Bralower, Paull and Leckie1998; Arenillas et al. Reference Arenillas, Arz, Grajales-Nishimura, Murillo-Muneton, Alvarez, Camargo-Zanguera, Molina and Rosales-Dominguez2006; Schulte et al. Reference Schulte, Alegret and Arenillas2010). However, this unconformity is not unique in the NW–SW Atlantic, Caribbean and Gulf of Mexico, where major unconformities are the norm throughout the Cretaceous (Keller et al. Reference Keller, Lyons, MacLeod and Officer1993), though most workers have concentrated on the KTB (Bralower, Paull & Leckie, Reference Bralower, Paull and Leckie1998; Watkins & Self-Trail, Reference Watkins and Self-Trail2005). Shorter pulses of erosion in the early Danian, as documented in this study, are also widespread (MacLeod & Keller, Reference MacLeod and Keller1991a , Reference MacLeod and Keller b ; Keller et al. Reference Keller, Stinnesbeck, Adatte and Stueben2003b ; this study). Moreover, impact spherules and frequently small iridium anomalies are present above the KTB unconformity in early Danian P1a sediments in Cuba and southern Mexico (Alegret et al. Reference Alegret, Arenillas, Arz, Diaz, Grajales-Nishimura, Melendez, Molina, Rojas and Soria2005; Arenillas et al. Reference Arenillas, Arz, Grajales-Nishimura, Murillo-Muneton, Alvarez, Camargo-Zanguera, Molina and Rosales-Dominguez2006), in Haiti, Belize, and Guatemala (Keller et al. Reference Keller, Adatte, Stinnesbeck, Stueben and Berner2001, Reference Keller, Stinnesbeck, Adatte, Holland, Stueben, Harting, De Leon and De La Cruz2003a ; Keller, Reference Keller, Evans, Horton and King2008) and in the Caribbean and N. Atlantic (this study). Impact spherules are abundant and restricted to the latest Maastrichtian zone CF1 (~130–150 ka below KTB) in NE Mexico and Texas where KTB sequences are more complete. Their absence below the KTB in the Atlantic, Caribbean, Gulf of Mexico, southern Mexico, Haiti, Belize and Guatemala is likely due to erosion linked to intensified current circulation and the impact (Keller et al. Reference Keller, Lyons, MacLeod and Officer1993; Watkins and Self-Trail, Reference Watkins and Self-Trail2005).
Bralower, Paull & Leckie, (Reference Bralower, Paull and Leckie1998) attributed any reworked microfossils, impact derived materials, lithic fragments and unconformities in the Caribbean and Gulf of Mexico to the Chicxulub impact and proposed that such deposits (labeled impact cocktail) are a reliable way to identify the KTB. All of the sections discussed in their paper show faunal and hiatus patterns similar to the Atlantic and Caribbean sections of this study where impact spherules are reworked in the early Danian and provide no evidence for a KTB age for the Chicxulub impact. Thus, caution must be exercised in interpreting iridium anomalies and impact spherules as KTB age markers for the Chicxulub impact, limestone conglomerates and breccias as related to the Chicxulub impact, and unconformities and hiatuses as caused by this impact. Sediment depositional patterns, impact spherule layers, minor iridium anomalies, short hiatuses and major unconformities through N. Atlantic, Caribbean, Gulf of Mexico and Central America are far too variable and complex for a simple one-size-fits-all interpretation.
8. Conclusions
Cretaceous–Tertiary sequences from Bass River (New Jersey), Blake Nose to Demerara Rise have been hailed as ultimate proof that the Chicxulub impact hit Yucatan at KTB time and caused the mass extinction. However, our detailed high-resolution biostratigraphic, faunal and stable isotope investigation reveals a fragmented sedimentary record for the KTB transition with a major unconformity that spans from the lower Danian across the KTB and well into the upper Maastrichtian, and impact spherules reworked and redeposited in early Danian Zone P1a (Figs 2–7). Similar early Danian records with multiple impact spherule layers and minor iridium anomalies have been published from the NW Atlantic, Haiti, Belize, Guatemala and southern Mexico, and Texas (Kramar et al. Reference Kramar, Stueben, Berner, Stinnesbeck, Philipp and Keller2001; Stueben et al. Reference Stueben, Kramar, Berner, Eckhardt, Stinnesbeck, Keller, Adatte and Heide2002; Keller, Reference Keller, Evans, Horton and King2008; Miller et al. Reference Miller, Sherrell, Browning, Field, Gallagher, Olsson, Sugarman, Tuorto and Wahyudi2010; Gertsch et al. Reference Gertsch, Keller, Adatte, Berner, Keller and Adatte2011). Therefore, sequences in these regions have very incomplete sedimentary records with reworked impact spherule layers that provide no insight into the age of the Chicxulub impact.
In contrast, KTB sequences in NE Mexico and Texas are more complete through to the upper Maastrichtian to early Danian and contain the stratigraphically oldest impact spherule deposits in the lower part of the latest Maastrichtian zone CF1, about 130–150 ka below the KTB. No significant species extinctions are associated with this impact spherule horizon (Keller et al. Reference Keller, Adatte, Berner, Pardo and Lopez-Oliva2009, Reference Keller, Abramovich, Adatte, Berner, Keller and Adatte2011). These results confirm our previous reports that the Chicxulub impact predates the KTB and cannot be the cause of the mass extinction, although it could have contributed through long term environmental changes. The real cause for the KTB mass extinction may have been India's Deccan Traps volcanism, where the main Deccan phase (80%) erupted during C29r over a short interval, coincident with rapid extinctions in planktic foraminifera, and ended at the KTB mass extinction (Keller et al. Reference Keller, Adatte, Bhowmick, Upadhyay, Dave, Reddy and Jaiprakash2012).
The following biostratigraphic results and age determination of missing intervals are based on this study:
(1) Bass River core, NJ: The KTB unconformity spans from the early Danian subzone P1a(2) to P1a(1), P0, CF1 and CF2 for a total of about 620 ka (Figs 2, 3, 11). The 6 cm thick impact spherule layer consists of eroded and redeposited debris.
(2) Blake Nose sites 1049A, 1049C, 1050C: The KTB unconformity spans from the early Danian subzone P1a(1) to upper Maastrichtian zones CF1 and CF2 for a total of about 400 ka (Figs 5, 6, 11). The 10 cm and 15 cm thick spherule layer was eroded and redeposited in the early Danian. The same unconformity is also present in Site 1050C where no spherules are present (Fig. 7).
(3) Caribbean sites 999B and 1001B: The KTB unconformity at Sites 999B and 1001B spans from Danian subzone P1a(1) to upper Maastrichtian zones CF1, CF2, CF3, CF4 for a total of about 3 Ma (Figs 2, 8). Impact spherules at Site 1001B are eroded and redeposited in the early Danian (no spherules in site 999B).
(4) Demerara Site 1259B: A KTB hiatus spans most of subzone P1a(1) and zone P0 (total of 150–200 ka). Spherules appear to be reworked into early Danian sediments.
(5) Short hiatuses also occur in the early Danian at the P1a(1)–P1a(2) and Pla(2)–P1b boundaries in all sections examined and are likely linked to climate and sea level fluctuations. Similar hiatuses have been identified in many early Danian sequences worldwide.
(6) Abrupt negative and positive δ13C shifts observed in Sites 1049A, 1049C, 1050C, 1259B support the KTB and early Danian hiatuses and reflect low primary productivity and delayed post-KTB recovery of the ecosystem into zone P1b.
(7) The widespread erosion and sedimentation pattern through the NW–SW Atlantic and Caribbean and the reworked impact spherules in lower Danian sediments rule this area out as indicator of the age of the Chicxulub impact.
Acknowledgements
We are grateful for the helpful comments and suggestions of the reviewers. Samples for this research were obtained from the Integrated Ocean Drilling Program (IODP), Bremen Core Repository (BCR), Germany, during two sampling visits in 2006 and 2011. We are grateful to Walter Hale, Repository Superintendent for his kind assistance during sampling. The material for this study is based upon work supported by the US National Science Foundation through the Continental Dynamics Program and Sedimentary Geology and Paleobiology Program under NSF Grants EAR-0207407 and EAR-0447171.