Subdivision of the Cretaceous system according to the ICS, as of 2017.[1]
The Cenomanian-Turonian boundary event, also known as the Cenomanian-Turonian extinction, Cenomanian-Turonian oceanic anoxic event (OAE 2), and referred to also as the Bonarelli event,[2] was one of two anoxicextinction events in the Cretaceous period. (The other being the earlier Selli event, or OAE 1a, in the Aptian.[3]) The Cenomanian-Turonian oceanic anoxic event is considered to be the most recent truly global oceanic anoxic event in Earth's geologic history.[4] Selby et al. in 2009 concluded the OAE 2 occurred approximately 91.5 ± 8.6 Ma,[5] though estimates published by Leckie et al. (2002) are given as 93–94 Ma.[6] The Cenomanian-Turonian boundary has been refined in 2012 to 93.9 ± 0.15 Ma.[7] There was a large carbon cycle disturbance during this time period,[8][9] signified by a large positive carbon isotope excursion.[10][11][12] However, apart from the carbon cycle disturbance, there were also large disturbances in the ocean's nitrogen,[13]oxygen,[14]phosphorus,[15][16][17]sulphur,[18] and iron cycles.[19]
The Cenomanian and Turonian stages were first noted by D'Orbigny between 1843 and 1852. The global type section for this boundary is located in the Bridge Creek Limestone Member of the Greenhorn Formation near Pueblo, Colorado, which are bedded with the Milankovitch orbital signature. Here, a positive carbon-isotope event is clearly shown, although none of the characteristic, organic-rich black shale is present. It has been estimated that the isotope shift lasted approximately 850,000 years longer than the black shale event, which may be the cause of this anomaly in the Colorado type section.[20] A significantly expanded OAE2 interval from southern Tibet documents a complete, more detailed, and finer-scale structures of the positive carbon isotope excursion that contains multiple shorter-term carbon isotope stages amounting to a total duration of 820 ±25 ka.[21]
The level is also known as the Bonarelli event because of 1-to-2-metre (3 ft 3 in to 6 ft 7 in) layer of thick, black shale that marks the boundary and was first studied by Guido Bonarelli (it) in 1891.[22] It is characterized by interbedded black shales, chert and radiolarian sands and is estimated to span a 400,000-year interval. Planktonic foraminifera do not exist in this Bonarelli level, and the presence of radiolarians in this section indicates relatively high productivity and an availability of nutrients.[23] In the Western Interior Seaway, the Cenomanian-Turonian boundary event is associated with the Benthonic Zone, characterised by a higher density of benthic foraminifera relative to planktonic foraminifera, although the timing of the appearance of the Benthonic Zone is not uniformly synchronous with the onset of the oceanic anoxic event and is thus cannot be used to consistently demarcate its beginning.[24]
Timeline
The total duration of OAE2 has been estimated at 0.82 ± 0.025 Myr[21] or 0.71 ± 0.17 Myr.[25] At high latitudes, the event lasted for a shorter time: only ~600 kyr.[26] Biodiversity patterns of planktic foraminifera indicate that the Cenomanian-Turonian extinction occurred in five phases. Phase I, which took place from 313,000 to 55,000 years before the onset of the anoxic event, witnessed a stratified water column and high planktonic foraminiferal diversity, suggesting a stable marine environment. Phase II, characterised by significant environmental perturbations, lasted from 55,000 years before OAE2 until its onset and witnessed a decline in rotaliporids and heterohelicids, a zenith of schackoinids and hedbergellids, a 'large form eclipse' during which foraminifera exceeding 150 microns disappeared, and the start of a trend of dwarfism among many foraminifera. This phase also saw an enhanced oxygen minimum zone and increased productivity in surface waters. Phase III lasted for 100,000 to 900,000 years and was coincident with the Bonarelli Level's deposition and exhibited extensive proliferation of radiolarians, indicative of extremely eutrophic conditions. Phase IV lasted for around 35,000 years and was most notable for the increase in the abundance of hedbergellids and schackoinids, being extremely similar to Phase II, with the main difference being that rotaliporids were absent from Phase IV. Phase V was a recovery interval lasting 118,000 years and marked the end of the 'large form eclipse' that began in Phase II; heterohelicids and hedbergellids remained in abundance during this phase, pointing to continued environmental disturbance during this phase.[27]
Causes
Climate change
Earth pronouncedly warmed just before the beginning of OAE2.[28] The Cenomanian-Turonian interval represents one of the hottest intervals of the entire Phanerozoic eon,[29] and it boasted the highest carbon dioxide concentrations of the Cretaceous period.[30] Even before OAE2, during the late Cenomanian, tropical sea surface temperatures (SSTs) were very warm, about 27-29 °C.[31] The onset of OAE2 was concurrent with a 4-5 °C rise in shelf sea temperatures.[32] Mean tropical SSTs during OAE2 have been conservatively estimated to have been at least 30 °C, but may have reached as high as 36 °C.[33] Minimum SSTs in mid-latitude oceans were >20 °C.[34] This exceptional warmth persisted until the Turonian-Coniacian boundary.[35]
One possible cause of this hothouse was sub-oceanic volcanism. During the middle of the Cretaceous period, the rate of crustal production reached a peak, which may have been related to the rifting of the newly formed Atlantic Ocean.[36] It was also caused by the widespread melting of hot mantle plumes under the ocean crust, at the base of the lithosphere, which may have resulted in the thickening of the oceanic crust in the Pacific and Indian Oceans. The resulting volcanism would have sent large quantities of carbon dioxide into the atmosphere, leading to an increase in global temperatures. Several independent events related to large igneous provinces (LIPs) occurred around the time of OAE2. A multitude of LIPs were active during OAE2: the Madagascar,[37][38]Caribbean,[39][40][41] Gorgona,[42]Ontong Java,[37] and High Arctic LIPs.[43][44][45] The abundance of LIPs at this time reflects a major overturning in mantle convection.[46] Trace metals such as chromium (Cr), scandium (Sc), copper (Cu) and cobalt (Co) have been found at the Cenomanian-Turonian boundary, which suggests that an LIP could have been one of the main basic causes involved in the contribution of the event.[47] The timing of the peak in trace metal concentration coincides with the middle of the anoxic event, suggesting that the effects of the LIPs may have occurred during the event, but may not have initiated the event. Other studies linked the lead (Pb) isotopes of OAE-2 to the Caribbean-Colombian and the Madagascar LIPs.[48] An osmium isotope excursion coeval with OAE2 strongly suggests submarine volcanism as its cause;[49] in the Pacific, an unradiogenic osmium spike began about 350 kyr before the onset of OAE2 and terminated around 240 kyr after OAE2's beginning;[50] the osmium isotope data from a highly expanded OAE2 interval in southern Tibet show multiple osmium excursions with the most pronounced one lagging the onset of OAE2 by ≈50 kyr that was probably related to the ocean connectivity change at ~94.5 Ma.[51] Osmium data also reveal that three distinct pulses of intense volcanism occurred ~60, ~270, and ~400 kyr after OAE2's onset, prolonging it.[52] Positive neodymium isotope excursions provide additional indications of pervasive volcanism as a cause of OAE2.[53] Enrichments in zinc further bolster and reinforce the existence of extensive hydrothermal volcanism,[54] as do extreme negative δ53Cr excursions.[55] The absence of geographically widespread mercury (Hg) anomalies resulting from OAE2 has been suggested to be because of the limited dispersal range of this heavy metal by submarine volcanism.[56] A modeling study performed in 2011 confirmed that it is possible that a LIP may have initiated the event, as the model revealed that the peak amount of carbon dioxide degassing from volcanic LIP degassing could have resulted in more than 90 percent global deep-ocean anoxia.[57]
Later on, when anoxia became widespread, the production of nitrous oxide, a greenhouse gas about 1,000 times more potent than carbon dioxide, drastically increased because of elevated nitrification and denitrification rates. This powerful positive feedback mechanism is what may have enabled extremely hot temperatures to persist in spite of the supercharged organic carbon burial associated with anoxic events.[58]
Plenus Cool Event
Large-scale organic carbon burial acted as a negative feedback loop that partially mitigated the warming effects of volcanic discharge of carbon dioxide, resulting in the Plenus Cool Event during the Metoicoceras geslinianum European ammonite biozone.[59] Global average temperatures fell to around 4 °C lower than they were pre-OAE2.[31] Equatorial SSTs dropped by 2.5–5.5 °C.[60] This cooling event was insufficient at completely stopping the rise in global temperatures. This negative feedback was ultimately overridden, as global temperatures continued to shoot up in sync with continued volcanic release of carbon dioxide following the Plenus Cool Event,[59] although this theory has been criticised and the warming after the Plenus Cool Event attributed to decreased silicate weathering instead.[61]
Ocean acidification
Within the oceans, the emission of SO2, H2S, CO2, and halogens would have increased the acidity of the water, causing the dissolution of carbonate, and a further release of carbon dioxide. Evidence of ocean acidification can be gleaned from δ44/40Ca increases coeval with the extinction event,[62][63][64] as well as coccolith malformation and dwarfism.[65] Ocean acidification was exacerbated by a positive feedback loop of increased heterotrophic respiration in highly biologically productive waters, elevating seawater concentrations of carbon dioxide and further decreasing pH.[66]
Anoxia and euxinia
When the volcanic activity declined, this run-away greenhouse effect would have likely been put into reverse. The increased CO2 content of the oceans could have increased organic productivity in the ocean surface waters. The consumption of this newly abundant organic life by aerobic bacteria would produce anoxia and mass extinction.[67] An acceleration of the hydrological cycle induced by warmer global temperatures drove greater fluxes of nutrient runoff into the oceans, fuelling primary productivity.[68][69] The global environmental disturbance that resulted in these conditions increased atmospheric and oceanic temperatures. Extreme hothouse conditions encouraged ocean stratification.[70] Boundary sediments show an enrichment of trace elements, and contain elevated δ13C values.[71][72] The positive δ13C excursion found at the Cenomanian-Turonian boundary is one of the main carbon isotope events of the Mesozoic. It represents one of the largest disturbances in the global carbon cycle from the past 110 million years. This δ13C excursion indicates a significant increase in the burial rate of organic carbon, indicating the widespread deposition and preservation of organic carbon-rich sediments and that the ocean was depleted of oxygen at the time.[73][74][75] Depletion of manganese in sediments corresponding to OAE2 provides additional strong evidence of severe bottom water oxygen depletion.[54] An increase in the abundance of the planktonic foraminifer Heterohelix provides further evidence still of anoxia.[76][52] The resulting elevated levels of carbon burial would account for the black shale deposition in the ocean basins.[71] The proto-North Atlantic in particular was a hotbed of carbon burial during OAE2 as it was in later, less severe anoxic events.[77]
Sulphate reduction increased during OAE2,[19] causing euxinia, a type of anoxia defined by sulphate reduction and hydrogen sulphide production, to occur during OAE2, as revealed by negative δ53Cr excursions,[78] positive δ98Mo excursions,[79] a low seawater molybdenum inventory,[80] and molecular biomarkers of green sulfur bacteria.[81][82][83]
OAE2 began on the southern margins of the proto-North Atlantic, from where anoxia spread across the rest of the proto-North Atlantic and then into the Western Interior Seaway (WIS) and the epicontinental seas of the Western Tethys.[84] Anoxic waters spread rapidly throughout the WIS due to marine transgression and a powerful cyclonic circulation resulting from an imbalance between precipitation in the north and evaporation in the south.[85] Anoxia was especially intense in the eastern North Sea, evidenced by its very positive δ13C values.[86] Thanks to persistent upwelling, some marine regions, such as the South Atlantic, were able to remain partially oxygenated at least intermittently.[87] Indeed, redox states of oceans vary geographically, bathymetrically and temporally during OAE2.[88]
Milankovitch cycles
It has been hypothesised that the Cenomanian-Turonian boundary event occurred during a period of very low variability in Earth's insolation, which has been theorised to be the result of coincident nodes in all orbital parameters. Barring chaotic perturbations in Earth's and Mars' orbits, the simultaneous occurrence of nodes of orbital eccentricity, axial precession, and obliquity on Earth occurs approximately every 2.45 million years.[89] Numerous other oceanic anoxic events occurred throughout the extremely warm greenhouse conditions of the Middle Cretaceous,[90] and it has been suggested that these Middle Cretaceous ocean anoxic events occurred cyclically in accordance with orbital cycle patterns.[89] The mid-Cenomanian Event (MCE), which occurred in the Rotalipora cushmani planktonic foraminifer biozone, has been argued to be another example supporting this hypothesis of regular oceanic anoxic events governed by Milankovitch cycles.[90] The MCE took place approximately 2.4 million years before the Cenomanian-Turonian oceanic anoxic event, roughly at the time when an anoxic event would be expected to occur given such a cycle.[89] Geochemical evidence from a sediment core in the Tarfaya Basin is indicative of the main positive carbon isotope excursion occurring during a prolonged eccentricity minimum. Carbon isotope shifts smaller in scale observed in this core likely reflected variability in obliquity.[91] Ocean Drilling Program Site 1138 in the Kerguelen Plateau yields evidence of a 20,000 to 70,000 year periodicity in changes in sedimentation, suggesting that either obliquity or precession governed the large-scale burial of organic carbon.[92] Within the OAE2 positive δ13C excursion, short eccentricity scale carbon isotope variability is documented in a significantly expanded OAE2 interval from southern Tibet;[21] periodic negative δ13C excursions paced by the short eccentricity cycle are easily detectable in southwestern Utah too.[93]
Enhanced phosphorus recycling
The phosphorus retention ability of seafloor sediments declined during OAE2,[15][94] revealed by a decline in reactive phosphorus species within OAE2 sediments.[95] The mineralisation of seafloor phosphorus into apatite was inhibited by the significantly lower pH of seawater and much warmer temperatures during the Cenomanian and Turonian compared to the present day, which meant that significantly more phosphorus was recycled back into ocean water after being deposited on the sea floor during this time. This would have intensified a positive feedback loop in which phosphorus is recycled faster into anoxic seawater compared to oxygen-rich water, which in turn fertilises the water, causes increased eutrophication, and further depletes the seawater of oxygen.[16] The influx of volcanically erupted and chemically weathered sulphate into the ocean also inhibited phosphorus burial by increasing hydrogen sulphide production,[96] which hinders the burial of phosphorus through sorption to iron oxyhydroxide phases.[18] OAE2 may have occurred during a peak in a 5-6 Myr cycle governing phosphorus availability; at this and other peaks in this oscillation, an increase in chemical weathering would have increased the marine phosphorus inventory and sparked a positive feedback loop of increasing productivity, anoxia, and phosphorus recycling that was only ended by a negative feedback of increased atmospheric oxygenation and wildfire activity that decreased chemical weathering, a feedback which operated on a much longer timescale.[17] Enhanced phosphorus recycling would have resulted in an abundance of nitrogen fixing bacteria, increasing the availability of yet another limiting nutrient and supercharging primary productivity through nitrogen fixation.[97] The ratio of bioavailable nitrogen to bioavailable phosphorus, which is 16:1 in the present, fell precipitously as the ocean transitioned from being oxic and nitrate-dominated to anoxic and ammonium-dominated.[58] A potent feedback loop of nitrogen fixation, productivity, deoxygenation, nitrogen removal, and phosphorus recycling was created.[13] Bacterial hopanoids indicate populations of nitrogen fixing cyanobacteria were high during OAE2, providing a rich supply of nitrates and nitrites.[98] Negative δ15N values reveal the dominance of ammonium through regenerative nutrient loops in the proto-North Atlantic.[99]
Decreased sulphide oxidation
In the present day, sulphidic waters are generally prevented from spreading throughout the water column by the oxidation of sulphide with nitrate. However, during OAE2, the inventory of seawater nitrate was lower, meaning that chemolithoautotrophic oxidation of sulphides with nitrates was inefficient at preventing the spread of euxinia.[100]
Sea level rise
A marine transgression in the latest Cenomanian resulted in an increase in average water depth, causing seawater to become less eutrophic in shallow, epicontinental seas. Turnovers in marine biota in such epicontinental seas have been suggested to be driven more so by changes in water depth rather than anoxia.[101] Sea level rise also contributed to anoxia by transporting terrestrial plant matter from inundated lands seaward, providing an abundant source of sustenance for eutrophicating microorganisms.[102]
Geological effects
Phosphate deposition
A phosphogenic event occurred in the Bohemian Cretaceous Basin during the peak of oceanic anoxia. Phosphorus liberation in the pore water environment, several centimetres below the interface between seafloor sediments and the water column, enabled the precipitation of phosphate through biological mediation by microorganisms.[103]
Increase in weathering
Strontium and calcium isotope ratios both indicate that silicate weathering increased over the course of OAE2. Because of its effectiveness as a carbon sink on geologic timescales, the uptick in sequestration of carbon dioxide by the lithosphere may have helped to stabilise global temperatures after global temperatures soared.[104] Particularly so at high latitudes, where the increase in weatherability was very pronounced.[105]
Biotic effects
Changes in oceanic biodiversity and its implications
The event brought about the extinction of the pliosaurs, and most ichthyosaurs. Coracoids of Maastrichtian age were once interpreted by some authors as belonging to ichthyosaurs, but these have since been interpreted as plesiosaur elements instead.[106]
Although the cause is still uncertain, the result starved the Earth's oceans of oxygen for nearly half a million years, causing the extinction of approximately 27 percent of marine invertebrates, including certain planktic and benthic foraminifera, mollusks, bivalves, dinoflagellates and calcareous nannofossils.[67] Planktonic foraminifera that dwelt in deeper waters were especially hard hit by the expansion of oxygen minimum zones.[107] In Whadi El Ghaib, a site in Sinai, Egypt, the foraminiferal community during OAE2 was low in diversity and dominated by taxa that were extremely tolerant of low salinity, anoxic water.[108] In the southeastern Indian Ocean, off the coast of Australia, the planktonic foraminifer Microhedbergella was highly abundant,[109] while Heterohelix thrived in reducing waters in the South Atlantic.[76][52] The alterations in diversity of various marine invertebrate species such as calcareous nannofossils are reflective and characteristic of oligotrophy and ocean warmth in an environment with short spikes of productivity followed by long periods of low fertility.[110] A study performed in the Cenomanian-Turonian boundary of Wunstorf, Germany, reveal the uncharacteristic dominance of a calcareous nannofossil species, Watznaueria, present during the event. Unlike the Biscutum species, which prefer mesotrophic conditions and were generally the dominant species before and after the C/T boundary event; Watznaueria species prefer warm, oligotrophic conditions.[111] In the Ohaba-Ponor section in Romania, the presence of Watznaueria barnesae indicates warm conditions, while the abundances of Biscutum constans, Zeugrhabdotus erectus, and Eprolithus floralis peak during cool intervals.[110] Sites in Colorado, England , France , and Sicily show an inverse relationship between atmospheric carbon dioxide levels and the size of calcareous nannoplankton.[112]Radiolarians also suffered heavy losses in OAE2, one of their highest diversity losses in the Cretaceous.[113]Rudist bivalves suffered high extinction rates combined with low origination rates during OAE2.[114]
The diversity of trace fossils sharply plummeted during the beginning of the Cenomanian-Turonian boundary event. The recovery interval after the anoxic event's conclusion features an abundance of Planolites and is characterised overall by a high degree of bioturbation.[115]
At the time, there were also peak abundances of the green algal groups Botryococcus and prasinophytes, coincident with pelagic sedimentation. The abundances of these algal groups are strongly related to the increase of both the oxygen deficiency in the water column and the total content of organic carbon. The evidence from these algal groups suggest that there were episodes of halocline stratification of the water column during the time. A species of freshwater dinocyst—Bosedinia—was also found in the rocks dated to the time and these suggest that the oceans had reduced salinity.[116][117]
Changes in terrestrial biodiversity
No major change in terrestrial ecosystems is known to have been synchronous with the marine transgression associated with OAE2, although the loss of freshwater floodplain habitat has been speculated to have possibly resulted in the demise of some freshwater taxa. In fossiliferous rocks in southwestern Utah, a local extirpation of some metatherians and brackish water vertebrates is associated with the later marine regression following OAE2 in the Turonian. Whatever the nature and magnitude of terrestrial extinctions at or near the Cenomanian-Turonian boundary was, it was most likely caused mainly by other factors than eustatic sea level fluctuations.[118] The effect of the ecological crisis on terrestrial plants has been concluded to have been inconsequential, in contrast to extinction events driven by terrestrial large igneous provinces.[119]
↑G. Bonarelli, Il territorio di Gubbio - Notizie geologiche, Roma 1891
↑G. Parisi, F. Piergiovanni e M. Marcucci, Il livello Bonarelli nell'area umbro-marchigiana, in Stratigrafia del Mesozoico e Cenozoico nell'area Umbro-Marchigiana, Roma, 1989
↑Bryant, Raquel; Belanger, Christina L. (19 January 2023). "Spatial heterogeneity in benthic foraminiferal assemblages tracks regional impacts of paleoenvironmental change across Cretaceous OAE2". Paleobiology49 (3): 431–453. doi:10.1017/pab.2022.47.
↑ 31.031.1Forster, Astrid; Schouten, Stephan; Moriya, Kazuyoshi; Wilson, Paul A.; Sinninghe Damsté, Jaap S. (14 March 2007). "Tropical warming and intermittent cooling during the Cenomanian/Turonian oceanic anoxic event 2: Sea surface temperature records from the equatorial Atlantic". Paleoceanography and Paleoclimatology22 (1): 1–14. doi:10.1029/2006PA001349. Bibcode: 2007PalOc..22.1219F.
↑O'Brien, Charlotte L.; Robinson, Stuart A.; Pancost, Richard D.; Sinninghe Damsté, Jaap S.; Schouten, Stefan; Lunt, Daniel J.; Alsenz, Heiko; Bornemann, André et al. (September 2017). "Cretaceous sea-surface temperature evolution: Constraints from TEX86 and planktonic foraminiferal oxygen isotopes". Earth-Science Reviews172: 224–247. doi:10.1016/j.earscirev.2017.07.012. Bibcode: 2017ESRv..172..224O.
↑ 37.037.1Scaife, J. D.; Ruhl, Micha; Dickson, A. J.; Mather, Tamsin A.; Jenkyns, Hugh C.; Percival, L. M. E.; Hesselbo, Stephen P.; Cartwright, J. et al. (1 November 2017). "Sedimentary Mercury Enrichments as a Marker for Submarine Large Igneous Province Volcanism? Evidence From the Mid-Cenomanian Event and Oceanic Anoxic Event 2 (Late Cretaceous)". Geochemistry, Geophysics, Geosystems18 (12): 4253–4275. doi:10.1002/2017GC007153. Bibcode: 2017GGG....18.4253S.
↑Du Vivier, Alice D. C.; Selby, David; Sageman, Bradley B.; Jarvis, Ian; Gröcke, Darren R.; Voigt, Silke (1 March 2014). "Marine 187Os/188Os isotope stratigraphy reveals the interaction of volcanism and ocean circulation during Oceanic Anoxic Event 2". Earth and Planetary Science Letters389: 23–33. doi:10.1016/j.epsl.2013.12.024. ISSN0012-821X.
↑Naber, T.V.; Grasby, S.E.; Cuthbertson, J.P.; Rayner, N.; Tegner, C. (16 December 2020). "New constraints on the age, geochemistry, and environmental impact of High Arctic Large Igneous Province magmatism: Tracing the extension of the Alpha Ridge onto Ellesmere Island, Canada" (in en). Geological Society of America Bulletin133 (7–8): 1695–1711. doi:10.1130/B35792.1. ISSN0016-7606.
↑Du Vivier, Alice D. C.; Jacobson, Andrew D.; Lehn, Gregory O.; Selby, David; Hurtgen, Matthew T.; Sageman, Bradley B. (15 April 2015). "Ca isotope stratigraphy across the Cenomanian–Turonian OAE 2: Links between volcanism, seawater geochemistry, and the carbonate fractionation factor". Earth and Planetary Science Letters416: 121–131. doi:10.1016/j.epsl.2015.02.001. Bibcode: 2015E&PSL.416..121D.
↑Schlanger, S. O.; Arthur, M. A.; Jenkyns, Hugh C.; Scholle, P. A. (1987). "The Cenomanian-Turonian Oceanic Anoxic Event, I. Stratigraphy and distribution of organic carbon-rich beds and the marine δ 13 C excursion". Geological Society, London, Special Publications26 (1): 371–399. doi:10.1144/GSL.SP.1987.026.01.24. ISSN0305-8719. Bibcode: 1987GSLSP..26..371S.
↑Poulton, Simon W.; Henkel, Susann; März, Christian; Urquhart, Hannah; Flögel, Sascha; Kasten, Sabine; Sinninghe Damsté, Jaap S.; Wagner, Thomas (1 November 2015). "A continental-weathering control on orbitally driven redox-nutrient cycling during Cretaceous Oceanic Anoxic Event 2". Geology43 (11): 963–966. doi:10.1130/G36837.1. Bibcode: 2015Geo....43..963P.
↑Monteiro, F. M.; Pancost, Richard D.; Ridgwell, Andy; Donnadieu, Yannick (15 December 2012). "Nutrients as the dominant control on the spread of anoxia and euxinia across the Cenomanian-Turonian oceanic anoxic event (OAE2): Model-data comparison". Paleoceanography and Paleoclimatology27 (4): 1–17. doi:10.1029/2012PA002351. Bibcode: 2012PalOc..27.4209M.
↑Fonseca, Carolina; Mendonça Filho, João Graciano; Lézin, Carine; de Oliveira, António Donizeti; Duarte, Luís V. (December 2019). "Organic matter deposition and paleoenvironmental implications across the Cenomanian-Turonian boundary of the Subalpine Basin (SE France): Local and global controls". International Journal of Coal Geology218: 103364. doi:10.1016/j.coal.2019.103364.