The Mesozoic terminated in boreal spring.

The Cretaceous-Palaeogene mass extinction around 66 million years ago was triggered by the Chicxulub asteroid impact on the present-day Yucatán Peninsula1,2. This event caused the highly selective extinction that eliminated about 76% of species3,4, including all non-avian dinosaurs, pterosaurs, ammonites, rudists and most marine reptiles. The timing of the impact and its aftermath have been studied mainly on millennial timescales, leaving the season of the impact unconstrained. Here, by studying fishes that died on the day the Mesozoic era ended, we demonstrate that the impact that caused the Cretaceous-Palaeogene mass extinction took place during boreal spring. Osteohistology together with stable isotope records of exceptionally preserved perichondral and dermal bones in acipenseriform fishes from the Tanis impact-induced seiche deposits5 reveal annual cyclicity across the final years of the Cretaceous period. Annual life cycles, including seasonal timing and duration of reproduction, feeding, hibernation and aestivation, vary strongly across latest Cretaceous biotic clades. We postulate that the timing of the Chicxulub impact in boreal spring and austral autumn was a major influence on selective biotic survival across the Cretaceous-Palaeogene boundary.

Main

The Cretaceous-Palaeogene (K–Pg) mass extinction event affected biodiversity with high but poorly understood taxonomic selectivity. Among archosaurs, for example, all pterosaurs and non-avian dinosaurs succumbed in the K–Pg mass extinction, while crocodilians and birds survived into the Palaeogene period[3,4]. Direct consequences of the impact, including impact glass fallout, large-scale forest fires and tsunamis, are geologically documented more than 3,500 km from the Chicxulub impact crater[5-8]. Although direct effects of the impact devastated a vast geographical area, the global mass extinction probably unfolded during its aftermath, which involved rapid climatic deterioration estimated to have lasted up to several thousands of years[9-11]. Whether seasonal timing of the onset of these marked changes affected the selectivity of the K–Pg extinction could not yet be established owing to the lack of suitable records.

The Tanis event deposit in North Dakota (USA) is an exceptional seiche deposit preserving a rich thanatocoenosis (that is, a mass death assemblage) of latest Cretaceous biota at the top of the Hell Creek Formation. The majority of macrofossils encountered at the Tanis locality represent direct casualties of the K–Pg bolide impact that were buried within the impact-induced seiche deposit[5]. Tens of minutes after the impact, the seiche agitated large volumes of water and soil in the estuary of the Tanis river[5]. As the seiche proceeded upstream, it advected bones, teeth, bivalves, ammonites, benthic foraminifera (Extended Data Fig. 1a–c) and plant matter in the suspended load while impact spherules rained down from the sky[5]. Within the thanatocoenotic accumulation, abundant acipenseriforms—sturgeons and paddlefishes—became oriented along the seiche flow directions and buried alive with numerous impact spherules in their gills[5] (Fig. 1, Extended Data Fig. 2a, b).

Extended Data Fig. 1

Benthic foraminifera and acipenseriforms in the field at the Tanis locality and their respective sampling locations.

a,b,c, and d, undetermined benthic foraminifera recovered from the Tanis deposit. Scale bars 1 mm. e. Sturgeon pectoral fin spines VUA.GG.2017.X-2743M and VUA.GG.2017.X-2744M, preserved in anatomical position, relative to their respective sturgeon carcasses in the Tanis deposit. f, Perichondral sturgeon bone sample VUA.GG.2017.MDX-3. g, perichondral sturgeon bone sample VUA.GG.2017.X-2743M. h, Perichondral sturgeon bone sample VUA.GG.2017.X-2744M. i, Dermal paddlefish bone sample VUA.GG.2017.X-2724. j, VUA.GG.2017.X-2733A in left lateral (left) and right lateral (right) view (sampling locations indicated between lines). The sediment matrix of sample VUA.GG.2017.X-2733A (bottom) also contains specimen VUA.GG.2017.X-2733B (covered). The fragile bone-bearing matrix was stabilised in epoxy resin prior to cutting, which obscured VUA.GG.2017.X-2733B from view. Scale bars 1 cm.

Fig. 1

Reconstruction of a paddlefish with impact spherules in the gill rakers.

a, Three-dimensional rendering of paddlefish FAU.DGS.ND.161.4559.T in left lateral view with the location of a higher-resolution scan (depicted in b) indicated (white outline). b, Three-dimensional rendering of the subopercular and gills in a with trapped impact spherules (yellow). Scale bars, 2 cm. Two-dimensional tomographic data and fully annotated three-dimensional renderings are provided in Extended Data Fig. 2. A three-dimensional animated rendering of FAU.DGS.ND.161.4559.T is provided as Supplementary Video 1.

Extended Data Fig. 2

PPC- SRμCT data of FAU.DGS.ND.161.4559.T, a partial paddlefish from the Tanis locality.

a, Orthogonal virtual thin sections (100 μm thick, average-value projections) obtained in front, top, and right view. b, Impact spherules in virtual thin sections of a, indicated with yellow circles. Scale bars a 1 mm. c, Three-dimensional rendering (in left lateral view) with virtual cross sections of d (blue), e (green), and f (red) indicated. d, Coronal virtual slice. e, Sagittal virtual slice. f, Axial virtual slice, brain-enveloping tissues indicated with red arrows. g, Three-dimensional rendering in right lateral view with anatomical labels. h, Three-dimensional rendering in left lateral view with anatomical labels.

During the Maastrichtian (that is, the last age of the Cretaceous), the climate of present-day North Dakota involved four seasons that were documented in tree-ring records recovered from other Upper Cretaceous sites in the Hell Creek Formation[12,13]. Tanis was located at approximately 50° N during the latest Cretaceous and experienced distinct seasonality in rainfall and temperature[14]. Regional air temperatures were reconstructed to range from 4–6 °C in winter up to an average of about 19 °C in summer[13,14]. To uncover the season of the K–Pg bolide impact, we analysed osteohistological records of acipenseriform bone apposition in three paddlefish dentaries and three sturgeon pectoral fin spines that were excavated at the Tanis site in 2017 (Extended Data Fig. 1d–j). These skeletal elements preserve unaltered growth records from embryonic development up to death, making them highly suitable for life history reconstructions[15,16].

Growth records of end-Cretaceous fishes

To trace appositional growth and pinpoint the season in which bone apposition terminated, we first assessed the preservation of bone growth patterns across the studied specimens. We prepared dermal bone slices of six acipenseriform specimens as microscopic slides and subjected these to osteohistological assessment, during which lines of arrested growth (LAGs) were easily recognized (Fig. 2). To corroborate the annual nature of the LAGs using virtual high-resolution osteohistology[17,18], three-dimensional (3D) volumes were produced with propagation phase-contrast synchrotron radiation micro-computed tomography[19] on beamline BM05 of the European Synchrotron Radiation Facility, France. The 3D nature of the synchrotron data enables optimal projection of the bone deposition pattern across multiple cross-sectional planes and resolved the exact relationship between seasonality and cyclical bone apposition in superb detail[20]. In addition, virtual osteohistology allowed us to visualize the seasonal fluctuations of osteocyte lacunar density and volume, which are poorly expressed in the physical 2D thin sections[18] (Fig. 3c, d) . The osteohistological data (Figs. 2, 3, Extended Data Figs. 3–6) were complemented with an incremental carbon isotope record extracted from one of the paddlefish dentaries (VUA.GG.2017.X-2724).

Fig. 2

Osteohistological thin sections of five acipenseriform fishes.

a–e, Thin sections in transmitted light of VUA.GG.2017.MDX-3 (a), VUA.GG.2017.X-2743M (b), VUA.GG.2017.X-2744M (c), VUA.GG.2017.X-2733A (d) and VUA.GG.2017.X-2733B (e), showing congruent pacing of bone apposition during the final years of life, terminating in spring. Red arrows indicate LAGs. Scale bars, 0.5 mm.

Fig. 3

Carbon isotope record alongside the incremental growth profiles across the dentary of paddlefish VUA.GG.2017.X-2724.

a, 𝛿13C record expressed as ‰ on the Vienna Pee Dee Belemnite (VPDB) reference scale. The colour gradient highlights the theoretical range between maximum values during seasonal (summer) trophic increase of 13C (yellow) and minimum values during trophic decrease of 13C (winter) (blue). b, Virtual thick section (average-value projection with 0.1-mm depth) showing growth zones during the favourable growth seasons and annuli and LAGs outside the favourable growth seasons. c, Cell density map[51] of a virtual thick section (minimum-value projection with 0.2-mm depth) showing fluctuating osteocyte lacunar densities and sizes, with higher densities and largest sizes recorded during the favourable growth seasons (orange) and lower densities and smaller sizes outside the favourable growth seasons (purple). A comparative image of a larger section of bone with scale is provided in Extended Data Fig. 6. d, Microscopic thin section in transmitted light showing LAGs (red arrows) and a single growth mark indicated (bracket) spanning the distance between two subsequent LAGs and including a zone and an annulus (Extended Data Fig. 10b). Scanning data visualized in b and c were obtained approximately 10-mm distal from the physically sectioned thin slice of d, which itself was located directly proximal to the thick section sampled for a. Scale bars, 1 mm. Corresponding osteohistological data of the other five sampled acipenseriform fishes are presented in Extended Data Figs. 3–5.

Extended Data Fig. 3

Osteohistology of acipenseriforms from the Tanis locality.

a, Thin section of paddlefish dentary VUA.GG.2017.X-2724 under transmitted light. b, Detail of VUA.GG.2017.X-2724 thin section (white box in a), scale bar 100 μm. c, Detail of VUA.GG.2017.X-2724 thin section (white box in b), scale bar 100 μm. d, Thin section of sturgeon pectoral fin spine VUA.GG.2017.MDX-3 under transmitted light. e, Detail of VUA.GG.2017.MDX-3 thin section (white box in d), scale bar 100 μm. f, Detail of VUA.GG.2017.MDX-3 thin section (white box in e), scale bar 100 μm. g, Thin section of paddlefish dentary VUA.GG.2017.X-2733A under transmitted light. h, Detail of VUA.GG.2017.X-2733A thin section (white box in g) with red arrows indicating Lines of Arrested Growth (LAGs), scale bar 100 μm. i, Detail of VUA.GG.2017.X-2724 thin section (white box in h), scale bar 100 μm.

Extended Data Fig. 6

Osteocyte lacunar density and volume distribution in paddlefish dentary VUA.GG.2017.X-2724 revealed by PPC-SRμCT.

a, Virtual thin section (100 μm thick, average-value projection) with osteohistological features indicated. b, Segmented osteocyte lacunar distribution (black dots). c, Osteocyte lacunar density ma 54= Sanchez et al., 2013 p54 with gradient scale. d, Virtual thick section (100 μm thick, minimum-value projection). e, Three-dimensional osteocyte lacunar distribution at d (depth circa 1230 μm), colour-coded by volume. f, Oscillating osteocyte lacunar volumes towards periosteal margin (10-point moving average) in dashed box in e. Successive annual growth climaxes (green arrows) and growth cessations (LAGs; red arrows) indicated. Typical annual maximum osteocyte lacunar volumes and highest osteocyte lacunar densities were not yet achieved in the year of death, indicating that growth ceased prior to the annual growth climax projected to occur during summer. Because the data resolution (voxel size of 4.35 μm: appropriate for assessing growth marks and osteocyte lacunar distributions) is sixfold lower than that demonstrated in detailed osteocyte lacunar volume reconstructions[18] (voxel size of 0.7 μm), these results should be considered qualitatively.

Extended Data Fig. 8

Elemental distribution maps of acipenseriform elements from the Tanis locality obtained with micro-X-ray fluorescence.

a, Ca, P, and Mn distribution in paddlefish dentary VUA.GG.2017.X-2724. b, Ca, P, and Mn distribution in sturgeon pectoral fin spine VUA.GG.2017.MDX-3. c, Ca, P, and Mn distribution in paddlefish dentaries VUA.GG.2017.X-2733A, VUA.GG.2017.X-2733B, and the surrounding sediment matrix. d, Ca, P, and Mn distribution in sturgeon pectoral fin spine VUA.GG.2017.X-2743M. e, Ca, P, and Mn distribution in sturgeon pectoral fin spine VUA.GG.2017.X-2744M. f, K, Si, and Fe distribution in paddlefish dentary VUA.GG.2017.X-2724. g, K, Si, and Fe distribution in sturgeon pectoral fin spine VUA.GG.2017.MDX-3. h, K, Si, and Fe distribution in paddlefish dentaries VUA.GG.2017.X-2733A, VUA.GG.2017.X-2733B and the surrounding sediment matrix. i, K, Si, and Fe distribution in sturgeon pectoral fin spine VUA.GG.2017.X-2743M. j, K, Si, and Fe distribution in sturgeon pectoral fin spine VUA.GG.2017.X-2744M.

Extended Data Fig. 10

δ13Csc and δ18Osc data of structurally-bound carbonate across the growth record in paddlefish dentary VUA.GG.2017.X-2724 and osteohistological diagram of paddlefish dentary VUA.GG.2017.X-2724.

a, Incremental records of the carbon (δ13Csc; left) and oxygen (δ18Osc; right) isotopic composition of structurally-bound carbonate in the apatitic matrix of VUA.GG.2017.X-2724 with uncertainty limits (1 s.d.). The seasonal cyclicity of δ13Csc that reflects alternations in seasonal food availability is, to a lesser extent, also expressed in the δ18Osc record. b, Schematic diagram of paddlefish dentary VUA.GG.2017.X-2724 with osteohistological features indicated. In the acipenseriform dermal bones examined in this study, annual growth cyclicity can be traced through Growth Marks (GMs). A GM spans a single growth cycle that typically lasts one year and can be divided into a zone, an annulus, and a LAG[20,52]. The zone is deposited during a period of relative rapid growth in the active or favourable growth season[20]. The annulus is subsequently formed when growth slows down towards the end of the growth season[20]. Finally, a LAG forms when growth periodically ceases until the next growth season starts and a new zone is deposited[20].

Extended Data Fig. 5

PPC-SRμCT virtual osteohistology of acipenseriforms from the Tanis locality at 4.53 μm voxel size.

a, Virtual thin section (100 μm thick, average-value projection; avp 100) of paddlefish dentary VUA.GG.2017.X-2724. b, Virtual thin section (avp 100) of sturgeon pectoral fin spine VUA.GG.2017.MDX-3. c, Virtual thin section (avp 100) of sturgeon pectoral fin spine VUA.GG.2017.X-2743M. d, Virtual thin section (avp 100) of sturgeon pectoral fin spine VUA.GG.2017.X-2744M. e, Thin section, virtual thin section, and virtual thick section of sturgeon pectoral fin spine VUA.GG.2017.MDX-3. e,I. Microscopic thin section under transmitted light showing LAGs. e,II. Virtual thin section (avp 100) showing LAGs. e,III. Virtual thick section (minimum-value projection with 200 μm depth (mvp 200)) showing oscillating osteocyte lacunar densities and sizes, with highest densities and largest sizes recorded during the favourable growth seasons and lowest densities and smallest sizes outside the favourable growth seasons. f. Thin section, virtual thin section, and virtual thick section of sturgeon pectoral fin spine VUA.GG.2017.X-2743M. f,I. Microscopic thin section under transmitted light showing LAGs. f,II. Virtual thin section (avp 100) showing LAGs. f,III. Virtual thick section (mvp 200) showing oscillating osteocyte lacunar densities and sizes, with highest densities and largest sizes recorded during the favourable growth seasons and lowest densities and smallest sizes outside the favourable growth seasons. g, Thin section, virtual thin section, and virtual thick section of sturgeon pectoral fin spine VUA.GG.2017.X-2744M.g,I. Microscopic thin section under transmitted light showing LAGs. g,II. Virtual thin section (avp 100) showing LAGs. g,III. Virtual thick section (mvp 200) showing oscillating osteocyte lacunar densities and sizes, with highest densities and largest sizes recorded during the favourable growth seasons and lowest densities and smallest sizes outside the favourable growth seasons. Scanning data visualised in e,II,III, f,II,III, and g,II,III, were obtained approximately 10 mm distal to the physical thin sections in e,I, f,I, and g,I, respectively. Red arrows indicate LAGs, green arrows indicate oncentrations of larger osteocyte lacunae. Scale bars in e, f, and g, 1 mm.

The tomographic data show that impact spherules associated with the paddlefish skeleton are present exclusively in its gill rakers[5] and are absent elsewhere in the preserved specimen (Fig. 1). The absence of impact spherules outside the gill rakers demonstrates that spherules were filtered out of the surrounding waters but had not yet proceeded into the oral cavity or further down the digestive tract, and had not impacted the fish carcases during perimortem exposure. Impact spherule accumulation in the gill rakers and the arrival of the seiche waves must therefore have occurred simultaneously[5], which implies that the acipenseriforms were alive and foraging during the bolide impact and the last minutes of the Cretaceous.

Well-conserved bone growth archives

The degree of preservation of the sampled acipenseriform bones was assessed using micro-X-ray fluorescence (Methods, Extended Data Figs. 7–9), which would reveal potential taphonomic elemental exchange that may have affected the primary stable isotope composition. The micro-X-ray fluorescence maps show that Fe and Mn oxides are present in the bone vascular canals and surrounding sediments (Extended Data Fig. 8), but have not invaded the bone apatite (Ca5(PO4,CO3)3(OH,F,Cl)). Detrital components, characterized by high concentrations of K and Si, remain restricted to the sediment matrix (Extended Data Fig. 8f–j). The bone apatite conserves a highly homogeneous distribution of P and Ca (Extended Data Fig. 9), which corroborates the unaltered preservation of these apatitic tissues. Skeletal remains of the paddlefishes and sturgeons thus experienced negligible diagenetic alteration, probably as a consequence of rapid burial and possibly aided by early Mn and Fe oxide seam formation[21,22]. The exquisite 3D preservation of delicate structures, including non-ossified tissues that originally enveloped the brain (Extended Data Fig. 2c–f), further demonstrates the excellent preservation of the fossils and absence of taphonomic reorganization[23].

Extended Data Fig. 7

Micro-X-ray fluorescence spectra of acipenseriform elements from the Tanis locality.

a, Elemental spectrum of paddlefish dentary VUA.GG.2017.X-2724. b, Elemental spectrum of sturgeon pectoral fin spine VUA.GG.2017.MDX-3. c, Elemental spectrum of paddlefish dentary VUA.GG.2017.X-2733A, VUA.GG.2017.X-2733B, and surrounding matrix. d, Elemental spectrum of sturgeon pectoral fin spine VUA.GG.2017.X-2743M. e, Elemental spectrum of sturgeon pectoral fin spine VUA.GG.2017.X-2744M.

Extended Data Fig. 9

Elemental heat maps of acipenseriform elements from the Tanis locality, obtained with micro-X-ray fluorescence, showing homogenous distributions of Ca and P.

a, Ca heat map of paddlefish dentary X-2723. b, Ca heat map of sturgeon pectoral fin spine VUA.GG.2017.MDX-3. c, Ca heat map of paddlefish dentaries VUA.GG.2017.X-2733A, VUA.GG.2017.X-2733B, and the surrounding sediment matrix. d, Ca heat map of sturgeon pectoral fin spine VUA.GG.2017.X-2743M. e, Ca heat map of sturgeon pectoral fin spine VUA.GG.2017.X-2744M. f, P heat map of paddlefish dentary VUA.GG.2017.X-2724. g, P heat map of sturgeon pectoral fin spine VUA.GG.2017.MDX-3. h, P heat map of paddlefish dentaries VUA.GG.2017.X-2733A, VUA.GG.2017.X-2733B, and the surrounding sediment matrix. i, P heat map of sturgeon pectoral fin spine VUA.GG.2017.X-2743M. j, P heat map of sturgeon pectoral fin spine VUA.GG.2017.X-2744M. Red and blue indicate higher,respectively lower abundance of Ca (a, b, c, d, and e) and P (f, g, h, i, and j).

Consistent records of a spring death

Paddlefish dentaries form through perichondral ossification around the Meckel’s cartilage[24]. Sturgeon pectoral fin spines consist of dermal bone—an intramembranous skeletal tissue that forms in the mesenchyme (mesodermal embryonic tissue)[25]. Unlike endochondral bone, perichondral and dermal bone do not originate through mineralization of cartilaginous precursors[26-28] but grow exclusively through incremental bone matrix apposition by secretion of a row of osteoblasts[24,26-28]. The thickness of one annual growth mark cumulatively spans a thick (favourable) growth zone, a thinner (slowly deposited) annulus and, ultimately, a LAG[20]. Our microscopic and virtual osteohistological data consistently show that the six fishes perished (that is, stopped growing) while forming a growth zone shortly after a LAG was deposited (Figs. 2, 3, Extended Data Figs. 3–6), which coincides with an early stage of the favourable growth season[20]. The outermost cortices of all six acipenseriform individuals studied here also exhibit increasing osteocyte lacunar densities and sizes towards their periosteal surfaces (Fig. 3c, Extended Data Figs. 5, 6). In all specimens, this density remained lower than the highest densities and average sizes recorded in previous years (Fig. 3c, Extended Data Figs. 3–6, 10b). As osteocyte lacunar density and size patterns were consistently cyclical across the preceding years during which they peaked at the climaxes of the growth seasons, the last recorded growth season had thus not yet climaxed at the time of death (Figs. 2, 3, Extended Data Figs. 3–6, 10b).

The inferred annual growth cycles are independently corroborated by a stable carbon isotope (𝛿13Csc) archive that recorded several years of seasonal dietary fluctuations in growing bone. Paddlefish VUA.GG.2017.X-2724 also yielded, in addition to this 𝛿13Csc archive, an oxygen isotope (𝛿18Osc) record across the final six years of its life (Supplementary Data Table 1, Extended Data Fig. 10a, Methods). The low and constant 𝛿18Osc values in VUA.GG.2017.X-2724 reflect exclusive inhabitation of freshwater environments by the paddlefishes. This implies that their osteohistological records must have captured seasonal variability rather than, for example, migration between saline and freshwater habitats. Although modern sturgeons are known to have anadromous lifestyles[29,30], this remains to be confirmed for the fossil sturgeons at Tanis, as isotopic data from sturgeon pectoral fin spines could not be secured (Methods, ‘Micromill’). Notably, the osteohistological records of all our sturgeons and paddlefishes converge on the same annual growth phase, despite their potential different lifestyles.

Like their modern-day relatives, the latest Maastrichtian paddlefishes of Tanis were filter feeders that presumably consumed copepods and other zooplankton[29-31]. These fishes probably experienced an annual feeding pattern, determined by fluctuating food availability, that peaked between spring and autumn[31]. During maximum productivity, ingested zooplankton enriches the growing skeleton of filter-feeding fishes with 13C relative to 12C[32,33]. Thus, the cyclically elevated 13C/12C ratios in paddlefish VUA.GG.2017.X-2724 (Fig. 3a) reflect distinct episodes of high food availability and consumption. Carbon isotope records across the growth record of Paddlefish VUA.GG.2017.X-2724 indicate that peak annual growth rate was not yet attained and the feeding season had thus not yet climaxed—corroborating a boreal spring death.

Implications for selective K–Pg survival

The Chicxulub bolide impact caused a global heat pulse that ignited widespread wildfires[9,34]. After this heat wave, the last boreal spring of the Mesozoic transitioned to a global impact winter[10]. Although a June timing for the K–Pg impact has been suggested on the basis of palaeobotanical indications for anomalous freezing in this region (Wyoming, USA)[35], the palaeobotanical identities, taphonomic inferences and stratigraphic assumptions underlying that conclusion have since all been refuted[36-39]. Moreover, post-impact cooling happened in the first months to decades following the K–Pg impact[10], which renders proxies registering post-impact freezing conditions asynchronous with the impact event itself.

A suite of impact-induced phenomena contributed to the K–Pg extinction on differing timescales[40,41]. In the days to months following the impact, its instantaneous effects, such as intense infrared radiation caused by ejecta reentry[34], resulting wildfires[9,34] and the spread of sulfurous aerosols leading to acid precipitation[42] must have predominantly afflicted the exposed continental environments. Although negotiating these hostile conditions would not have guaranteed survival, an early clade-wide eradication would always have meant immediate extinction[41].

The seasonal timing of the catastrophic end-Cretaceous bolide impact places the event at a particularly sensitive stage for biological life cycles in the Northern Hemisphere. In many taxa, annual reproduction and growth take place during spring. Species with longer incubation times, such as non-avian reptiles, including pterosaurs and most dinosaurs, were arguably more vulnerable to sudden environmental perturbations than other groups[43] (for example, birds). Southern Hemisphere ecosystems, which were struck during austral autumn, appear to have recovered up to twice as fast as Northern Hemisphere communities[44], consistent with a seasonal effect on biotic recovery.

Subterranean sheltering conceivably contributed to the cynodont survival of the Permo-Triassic (PT) crisis[45]. Similarly, large-scale wildfires raging across the Southern Hemisphere[9,34,41] may have been evaded by hibernating mammals that were already sheltered in burrows[34,41] in anticipation of austral winter. Additional modes of seasonal dormancy, torpor and/or aestivation, which are nowadays practised by various mammals[46,47] as well as certain amphibians, birds and crocodilians[48], could have facilitated further underground survival. In the aftermath of the K–Pg event, ecological networks collapsed from the bottom up. Floral necrosis[9] and extinction immediately affected species dependent on primary producers, while some animals capable of exploiting alternative resources—for example, certain birds and mammals[49,50]—persisted.

Conclusions

Seasonal timing of the Chicxulub impact in boreal spring and austral autumn will aid in further calibrating evolutionary models exploring the selectivity of the K–Pg extinction and the asymmetry in extinction and recovery patterns between the two hemispheres. Decoupling short- and long-term effects of the bolide impact on the K–Pg mass extinction will also aid in identifying extinction risks and modes of ecological deterioration caused by the forthcoming global climate change. The uniquely constrained Tanis site[5] offers valuable proxies for reconstructing the environmental, climatological and biological conditions that prevailed locally when the Mesozoic ended.

Methods

Fieldwork

Excavation at the Tanis locality in south-western North Dakota took place between 10 August and 20 August 2017. Sections of dentaries of paddlefishes and pectoral fin spines of sturgeons were collected in the field for histological study.

Thin sectioning

Four out of the six samples were excavated from the sediment matrix. These included all sturgeon pectoral fin spines (VUA.GG.2017.X-2743M, VUA.GG.2017.X-2744M, and VUA.GG.2017.MDX-3) and one of the paddlefish dentaries (VUA.GG.2017.X-2724). Paddlefish dentaries VUA.GG.2017.X-2733A and VUA.GG.2017.X-2733B, belong to two individuals, were uncovered aligned to each other and fractured upon discovery. To avoid further damage, the samples were embedded in epoxy resin prior to thin sectioning. All specimens were cut with a diamond saw and polished to obtain microscopic thin sections (about 50-μm thick) and thick sections for micromilling (about 200-μm thick). See Extended Data Fig. 1e–j for images of the specimens and the sampling locations.

Osteohistological analysis

In the acipenseriform dermal bones examined in this study, annual growth cyclicity can be traced through growth marks (GMs).

A GM spans a single growth cycle that typically lasts one year and can be divided into a zone, an annulus, and a LAG[20,52]. The zone is deposited during a period of relative rapid growth in the active or favourable growth season[20]. The annulus is subsequently formed when growth slows down towards the end of the growth season[20]. Finally, a LAG forms when growth periodically ceases until the next growth season starts and a new zone is deposited[20].

During the formation of a growth zone, the density and volumes of osteocyte lacunae (OL; subcircular dark features in Extended Data Fig. 10a) initially increase when growth accelerates. Subsequently, towards and into the annulus, OL density and volume decrease as growth slows down[18]. Because a LAG coincides with a temporary arrest of local osteogenesis, it is only expressed when deposition of a new growth zone has commenced. All six studied specimens show a LAG relatively close to the outermost partial growth zone.

In fossil bone, LAGs often appear as sharply defined dark lines[53] that typically constitute a poorly coherent interface between adjacent bone layers, thus facilitating (local) delamination between adjacent cortical layers[53]. During fossilization, percolation products can accumulate in these gaps and thereby (locally) accentuate the LAGs[51,53] (figure 31.3G of ref. [52]). Based on this well-understood expression of LAGs (that we recognize from our own experience as well; S.S. personal observation), we have consistently identified the LAGs as locally stained dark lines that may be associated with circumferentially propagated cracked surfaces which are oriented parallel to the periosteal deposits.

Besides cyclical seasonal factors that synchronize GM accretion, stress may induce additional diapause stages that result in supplementary marks within a single year[54]. Cessation of growth for the duration of several weeks can provoke the formation of a LAG[54]. However, such non-cyclical marks “tend to be haphazard rather than regular (that is, they do not reflect a particular spacing or rhythm)” and do not encircle the cortex of the skeletal element but “tend to be locally confined to an arc”[55].

As the studied bones yield only regularly spaced GMs along their complete circumference, we confidently identify the preserved GMs as annual cycles. Moreover, the fluctuating quantified density and volumes of osteocyte lacunae (Extended Data Fig. 6d–f) and the carbon isotopic record (Fig. 3a, Extended Data Fig. 10a) across the final seven years of growth of VUA.GG.2017.X-2724 are exclusively consistent with the identification of annual LAGs in corresponding physical thin sections. In all studied specimens, bone growth terminated during the process of zonal bone growth.

Micro-X-ray fluorescence

Fragments of the paddlefish and sturgeon samples that remained after thin sectioning were analysed with microX-ray fluorescence. High-resolution elemental mapping was conducted using a Bruker M4 Tornado 2D spectrometer at 50 kV and 600 μA, without a filter, and at an acquisition rate of 20 μm per 5 ms at the Vrije Universiteit Brussel.

Micromill

The growth increments were sampled in the thick sections (about 200-μm thick) at the highest possible accuracy using a Micromill (Merkantek). Drill transects were assigned in the accompanying software and after each individual sample was collected, the drill bit was cleaned with ethanol. Not all thick sections were suitable for micromilling. The lobed anatomy of the sturgeon fin spines (VUA.GG.2017.X-2743M and VUA.GG.2017.X-2744M) proved too complex to reliably sample single growth increments with the micromill. Paddlefish dentaries VUA.GG.2017.X-2733A and VUA.GG.2017.X-2733B only exposed a few growth lines that were too narrow to sample with the micromill. Sturgeon pectoral fin spine VUA.GG.2017.MDX-3 and paddlefish dentary VUA.GG.2017.X-2724 were sampled up to the outermost growth increment.

Stable isotope analysis

Micromilled hydroxyapatite samples of specimen VUA.GG.2017.X-2724 weighing about 50 μg were placed in Exetainer vials (Labco) and flushed with purified helium gas. For reference, the analysed amounts of structural carbonate are equivalent to anout 5 μg of CaCO3. Orthophosphoric acid was subsequently added and allowed to react for 24 h at 45 °C. VUA.GG.2017.MDX-3 was routinely analysed with a Thermo Finnigan Deltaplus mass spectrometer connected to a Thermo Finnigan GasBench II at the Earth Sciences Stable Isotope Laboratory (Vrije Universiteit, Amsterdam). However, the amount of CO2 generated was found to be too small to permit reliable isotopic determinations. To alleviate this, the GasBench was provisionally interfaced with a cold trap in which the CO2 was frozen with liquid nitrogen during a 2 min period. After trapping for 2 min, an accurate single-pulse measurement was performed for each of the apatitic samples and standards. Each isotopic sample determination was preceded by six pulses of monitoring CO2 with a calibrated isotopic composition to assure stable conditions of the mass spectrometer. The isotopic measurements of the weighted micromilled samples were bracketed by the analyses of the inter-laboratorial apatite standard (Ag-Lox) to account for the linearity effect[56]. After corrections, the uncertainties for 𝛿13C and 𝛿18O of the Ag-Lox (n = 4) were 0.16 ‰ and 0.39 ‰ (1 s.d.) respectively. Although the amount of extracted and analysed structural carbonate remains insufficient for optimal isotopic determination, the relatively large recovered 𝛿13C variability still yields a meaningful record across the appositional bone archive. The 𝛿18O values of structural carbonate, unlike those of phosphate (PO4)[57], do not offer a sensitive palaeo-environmental proxy for accurate seasonal temperature reconstructions[58]. However, the relatively constant 𝛿18O values of structural carbonate precludes large 𝛿18O changes in ambient water, such as shifts between freshwater and saline environments.

Propagation phase-contrast synchrotron radiation micro-computed tomography

Paddlefish specimen FAU.DGS.ND.161.4559.T lacks the paddle-shaped rostrum and all aspects caudal to the pectoral girdle. FAU.DGS.ND.161.4559.T was provided by the Palm Beach Museum of Natural History. Data acquisition took place in May 2018 on Beamline BM05 of the European Synchrotron Radiation Facility, Grenoble, France[59]. The complete specimen was scanned at an average energy of 132 keV using the white beam of BM05 filtered with 0.4 mm of Mo and 9 mm of Cu. The detector was composed of a 2-mm-thick LuAG:Ce scintillator optically coupled to a PCO edge 4.2 CLHS sCMOS camera. The resulting voxel size was 43.5 µm. To obtain sufficient propagation phasecontrast, the distance between the sample and the detector was set at 5 m. A total of 205 scans, each consisting of 5,000 projections taken at 7-ms intervals, were performed with a vertical displacement of 1.4 mm at a vertical field of view of 2.8 mm to ensure a double scan of the complete samples. Scans were performed in half-acquisition mode to enlarge the lateral field of view. The volume was reconstructed using a single-distance phase retrieval algorithm coupled with filtered back projection as implemented in the ESRF software PyHST2. Vertical concatenation, 16-bit conversion, and ring artefact corrections were performed using MATLAB scripts developed in-house. The gill region and impact spherules were subsequently scanned at a voxel size of 13.67 μm (filters: 0.4 mm of Mo and 6 mm of Cu, scintillator: LuAG:Ce, 500-μm thick, detected energy: 166 keV, propagation distance: 2.5 m). The samples were scanned in half-acquisition mode in two columns of 77 scans, each consisting of 4,998 projections with exposure times of 0.05 s, that were laterally concatenated after reconstruction. Finally, sample (VUA.GG.2017.X-2724) from the paddlefish dentaries and (VUA.GG.2017.MDX-3, VUA.GG.2017.X-2743M and VUA.GG.2017.X-2744M) of the sturgeon pectoral fin spines were scanned at 4.35 µm voxel size for osteohistological analysis[60] (filters: 3.5 mm of Al plus 11 bars Al with a diameter of 5 mm, scintillator: LuAG:Ce scintillator, 500-µm thick, detected energy: 92 keV, propagation distance: 1.5 m). The samples were scanned in half-acquisition mode in one single column of 22 scans, each consisting of 4,998 projections with exposure times of 60 ms.

Digital 3D extraction of the bones and impact spherules was performed in VGStudio MAX 3.2 (Volume Graphics). VGStudio MAX 3.2 furthermore enabled the creation of virtual thick sections of the osteohistological samples through the ‘thick slab-mode’, which captures the maximum, average, or minimum, grey-level values along the desired field depth. Virtual thick sections were obtained from the average grey-level values at a thickness of 100 µm following optimal 3D alignment of the annuli and LAGs. Additional virtual thick sections were created from the minimum grey-level values at a thickness of 200 µm to best resolve the sizes and distributions of osteocyte lacunae. A coloured map of the density of the osteocyte lacunar distribution was created with a Gaussian filter[51]. Finally, we visualized the annual cyclicity of osteocyte lacunar volumes[18] in paddlefish dentary VUA.GG.2017.X-2724. As the resolution of our data (voxel size of 4.35 μm; appropriate for assessing GMs and osteocyte lacunar distributions) is sixfold lower than that used for earlier osteocyte lacunar volumetric quantification in fish bones[18] (voxel size of 0.7 μm), our result should be considered with appropriate care. Closely spaced (large) osteocyte lacunae may occasionally be conjoined and additional phenomena in the broad size range of osteocyte lacunae may be incidentally included in the visualized distribution. Moreover, in tomographic data, osteocyte lacunae are delimited by slight colour gradients (rather than discrete lines) that scale with voxel size. Because the outermost feature fringe contributes disproportionally to recovered volumes, these values are somewhat skewed relative to the original osteocyte lacunar volumes, which likely produces exaggerated volume values. Therefore, although all rendered features were extracted with a single thresholding operation and relative patterns are conservatively retained, absolute volume values are best considered in a comparative context.