Biometric characteristics of Orbulina universa (d'Orbigny) were used to differentiate two morphotypes present in sediment trap samples collected from the Cariaco Basin, Venezuela. Specifically, wall thickness and weight-area relationships were used to separate shells into thin (Mthin) and thick (Mthick) morphotypes. Mthick (mean thickness=19-41μm) comprises 75% of the total O. universa in these samples and has morphometric characteristics similar to that of the previously described Type I Caribbean genotype, whereas Mthin (mean thickness=6-22μm) is comparable to the Type III Mediterranean genotype. The flux of Mthick increases during periods of upwelling, whereas Mthin flux shows no systematic relationship with changing hydrographic regimes in the basin. The δ18O and δ13C of Mthick are on average 0.34% higher and 0.38% lower, respectively, than those of Mthin, suggesting that they calcify their final spherical chamber at different depths in the water column and/or differ in their vital effects on shell geochemistry. Additionally, the absolute offset in the stable isotopic compositions of the two morphotypes varies as a function of surface ocean stratification. During periods of upwelling, the δ18O and δ13C absolute offsets between Mthin and Mthick are on average -0.15 and 0.25%, whereas average δ18O and δ13C absolute offsets between these two morphotypes increase to -0.54 and 0.41% during non-upwelling periods. We attribute the higher δ18O and δ13C offsets during non-upwelling periods to a difference in final chamber calcification depth and the combined effects of temperature (δ18O), δ13CDIC (δ13C), irradiance (δ13C) and [CO32-] (δ18O and δ13C). These data provide field evidence that thin and thick morphotypes of O. universa likely experience different environmental conditions during the formation of their final chamber and, therefore, should not be combined in geochemical analyses for reconstructing past surface ocean conditions. Finally, we introduce a new proxy for reconstructing past surface ocean stratification changes via the use of Mthin and Mthick δ18O and δ13C offsets.

Mass extinction at the Cretaceous-Paleogene (K-Pg) boundary coincides with the Chicxulub bolide impact and also falls within the broader time frame of Deccan trap emplacement. Critically, though, empirical evidence as to how either of these factors could have driven observed extinction patterns and carbon cycle perturbations is still lacking. Here, using boron isotopes in foraminifera, we document a geologically rapid surface-ocean pH drop following the Chicxulub impact, supporting impact-induced ocean acidification as a mechanism for ecological collapse in the marine realm. Subsequently, surface water pH rebounded sharply with the extinction of marine calcifiers and the associated imbalance in the global carbon cycle. Our reconstructed water-column pH gradients, combined with Earth system modeling, indicate that a partial ∼50% reduction in global marine primary productivity is sufficient to explain observed marine carbon isotope patterns at the K-Pg, due to the underlying action of the solubility pump. While primary productivity recovered within a few tens of thousands of years, inefficiency in carbon export to the deep sea lasted much longer. This phased recovery scenario reconciles competing hypotheses previously put forward to explain the K-Pg carbon isotope records, and explains both spatially variable patterns of change in marine productivity across the event and a lack of extinction at the deep sea floor. In sum, we provide insights into the drivers of the last mass extinction, the recovery of marine carbon cycling in a postextinction world, and the way in which marine life imprints its isotopic signal onto the geological record.

The cause of the end-Cretaceous mass extinction is vigorously debated, owing to the occurrence of a very large bolide impact and flood basalt volcanism near the boundary. Disentangling their relative importance is complicated by uncertainty regarding kill mechanisms and the relative timing of volcanogenic outgassing, impact, and extinction. We used carbon cycle modeling and paleotemperature records to constrain the timing of volcanogenic outgassing. We found support for major outgassing beginning and ending distinctly before the impact, with only the impact coinciding with mass extinction and biologically amplified carbon cycle change. Our models show that these extinction-related carbon cycle changes would have allowed the ocean to absorb massive amounts of carbon dioxide, thus limiting the global warming otherwise expected from postextinction volcanism.