Search

Ocean acidification (OA) during the Paleocene-Eocene thermal maximum (PETM) likely caused a biocalcification crisis. The calcium isotope composition (δ44/40Ca) of primary carbonate producers may be sensitive to OA. To test this hypothesis, we constructed the first high-resolution, high-precision planktic foraminiferal δ44/40Ca records before and across the PETM. The records employ specimens of Morozovella spp. collected from Ocean Drilling Program Sites 1209 (Shatsky Rise, Pacific Ocean) and 1263 (Walvis Ridge, Atlantic Ocean). At Site 1209, δ44/40Ca values start at −1.33‰ during the Upper Paleocene and increase to a peak of −1.15‰ immediately before the negative carbon isotope excursion (CIE) that marks the PETM onset. Values remain elevated through the PETM interval and decrease into the earliest Eocene. A shorter-term record for Site 1263 shows a similar trend, although δ44/40Ca values are on average 0.22‰ lower and decrease shortly after the CIE onset. The trends support neither diagenetic overprinting, authigenic carbonate additions, nor changes in the δ44/40Ca value of seawater. Rather, they are consistent with a kinetic isotope effect, whereby calcite δ44/40Ca values inversely correlate with precipitation rate. Geologically rapid Ca isotope shifts appear to reflect the response of Morozovella to globally forced changes in the local carbonate geochemistry of seawater. All data combined suggest that the PETM-OA event occurred near the peak of a gradual reduction in seawater carbonate ion concentrations during a time of elevated atmospheric pCO2, potentially driven by North Atlantic igneous province emplacement.

Large igneous province (LIP) eruptions are hypothesized to trigger biocalcification crises. The Aptian nannoconid crisis, which correlates with emplacement of the Ontong Java Plateau and Oceanic Anoxic Event 1a (OAE 1a, ca. 120 Ma), represents one such example. The Ca isotope (δ44/40Ca) system offers potential to detect biocalcification fluctuations in the rock record because Ca isotope fractionation is sensitive to precipitation rate. However, other primary and secondary processes, such as input-output flux perturbations and early diagenesis, can produce similar signals. Here, we exploit emergent properties of the stable Sr isotope (δ88/86Sr) system to resolve the origin of δ44/40Ca variability during OAE 1a. This study reports high-precision thermal ionization mass spectrometry (TIMS) δ44/40Ca, δ88/86Sr, and 87Sr/86Sr records for Hole 866A of Ocean Drilling Program Leg 143 drilled in Resolution Guyot, mid-Pacific Ocean. The samples span ∼27 m.y. from the Barremian (ca. 127 Ma) to the Albian (ca. 100 Ma). The δ44/40Ca and δ88/86Sr secular trends differ from the 87Sr/86Sr record but mimic each other. δ44/40Ca and [Sr], as well as δ44/40Ca and δ88/86Sr, strongly correlate and yield slopes predicted for kinetic control, which demonstrates that variable mass-dependent fractionation rather than end-member mixing dominated the isotopic relationship between carbonates and seawater. Positive δ44/40Ca and δ88/86Sr shifts that begin before OAE 1a and peak within the interval are consistent with reduced precipitation rates. All results combined point to a cascade of effects on rate-dependent Ca and Sr isotope fractionation, which derive from the dynamic interplay between LIP eruptions and biocalcification feedbacks.

Paleoclimate during the Cenomanian-Turonian Oceanic Anoxic Event 2 (OAE 2, 94.5–93.9 Ma) was characterized by elevated atmospheric CO2 concentrations and peak global temperatures. In this study, we employ δ66Zn measured on samples spanning OAE 2 in an expanded hemipelagic section in Tibet to trace changes in the major fluxes that influence biogeochemical cycles. The prominent feature of the δ66Zn record in the studied section is a continuous decrease from ∼1‰ at the onset of OAE 2 to a minimum of ∼0.2‰ within the Plenus Cold Event (ca. 94.3 Ma), followed by a stepwise recovery through the upper part of OAE 2. The negative shift in δ66Zn corresponds with higher terrigenous inputs, as revealed by previously published detrital index and TOC/TN (total organic carbon to total nitrogen) ratio records, and covaries with a notable decreasing trend recorded in compiled pCO2 data of different basins. We propose that influx of isotopically light Zn from weathered volcanic rocks associated with submarine large igneous provinces and/or (sub)tropic Indian continental volcanics is likely responsible for the δ66Zn decrease. We infer that the recovery of δ66Zn was caused by continued high primary production and an inevitable decline in the flux of light Zn as volcanic terrains were progressively weathered. The ultimate cessation of OAE 2 may have been a consequence of the same effect, with the nutrient supply from weathering reaching a minimum threshold to maintain productivity-anoxia feedback.

Carbon dioxide release during Deccan Traps volcanism and the Chicxulub impact likely contributed to the Cretaceous-Paleogene (K-Pg) mass extinction; however, the intensity and duration of CO2 input differed between the two events. Large and rapid addition of CO2 to seawater causes transient decreases in pH, [CO32–], and carbonate mineral saturation states. Compensating mechanisms, such as dissolution of seafloor sediment, reduced biomineralization, and silicate weathering, mitigate these effects by increasing the same parameters. The calcium isotope ratios (δ44/40Ca) of seawater and marine carbonates are hypothesized to respond to these perturbations through weathering/carbonate deposition flux imbalances and/or changes in fractionation between carbonate minerals and seawater. We used a high-precision thermal ionization mass spectrometry method to measure δ44/40Ca values of aragonitic bivalve and gastropod mollusk shells from the K-Pg interval of the López de Bertodano Formation on Seymour Island, Antarctica. Well-preserved shells spanning the late Maastrichtian (ca. 67 Ma) to early Danian (ca. 65.5 Ma) have δ44/40Ca values ranging from −1.89‰ to −1.57‰ (seawater [sw]). Shifts in δ44/40Ca inversely correlate with sedimentological indicators of saturation state. A negative excursion begins before and continues across the K-Pg boundary. According to a simple mass-balance model, neither input/output flux imbalances nor change in the globally integrated bulk carbonate fractionation factor can produce variations in seawater δ44/40Ca sufficient to explain the measured trends. The data are consistent with a dynamic molluscan Ca isotope fractionation factor sensitive to the carbonate geochemistry of seawater. The K-Pg extinction appears to have occurred during a period of carbonate saturation state variability caused by Deccan volcanism.

Cretaceous oceanic anoxic event 2 (OAE2) is thought to have been contemporary with extensive volcanism and the release of large quantities of volcanic CO2 capable of triggering marine anoxia through a series of biogeochemical feedbacks. High-resolution reconstructions of atmospheric CO2 concentrations across the initiation of OAE2 suggest that there were also two distinct pulses of CO2 drawdown coeval with increased organic carbon burial. These fluctuations in CO2 likely led to significant climatic changes, including fluctuations in temperatures and the hydrological cycle. Paleofire proxy records suggest that wildfire was a common occurrence throughout the Cretaceous Period, likely fueled by the estimated high atmospheric O2 concentrations at this time. However, over geological time scales, the likelihood and behavior of fire are also controlled by other factors such as climate, implying that CO2-driven climate changes should also be observable in the fossil charcoal record. We tested this hypothesis and present a high-resolution study of fire history through the use of fossil charcoal abundances across the OAE2 onset, and we compared our records to the estimated CO2 fluctuations published from the same study sites. Our study illustrates that inferred wildfire activity appears to relate to changes in CO2 occurring across the onset of OAE2, where periods of CO2 drawdown may have enabled an increase in fire activity through suppression of the hydrological cycle. Our study provides further insight into the relationships between rapid changes in the carbon cycle, climate, and wildfire activity, illustrating that CO2 and climate changes related to inferred wildfire activity can be detected despite the estimated high Cretaceous atmospheric O2 concentrations.

Proximal marine strata of the North American Western Interior Basin (WIB) preserve a rich record of biotic turnover during Oceanic Anoxic Event 2 (OAE2; ca. 94 Ma), a pronounced Late Cretaceous carbon cycle perturbation interpreted to reflect global warming, widespread hypoxia, and possible ocean acidification. To develop a more robust synthesis of paleobiologic and geochemical data sets spanning this Earth-life transition, we drilled the 131 m Smoky Hollow #1 Core (SH#1), on the Kaiparowits Plateau of southern Utah, USA, recovering the Cenomanian–Turonian Boundary (CTB) interval in the Tropic Shale Formation. A 17.5 m positive excursion in high-resolution bulk carbon isotope chemostratigraphy (δ13Corg) of SH#1 characterizes the most expanded OAE2 record recovered from the mid-latitudes of the WIB.Depleted values in a paired carbonate carbon isotope (δ13Ccarb) chemostratigraphy cyclically punctuate the OAE2 excursion. These depletions correspond to intervals in the core with a higher degree of carbonate diagenesis and correlate well to an existing sequence stratigraphic framework of flooding surfaces in the shoreface facies of the Markagunt Plateau (∼100 km west). We detect statistically significant evidence for astronomical cycles in the δ13Ccarb data set, imparted by diagenesis at flooding surfaces, and develop a floating astronomical time scale (ATS) for the study interval. Stable eccentricity cycles (405 k.y.) align with stratigraphic sequences and associated trends in sedimentation rate, and short eccentricity cycles (∼100 k.y.) pace nested parasequences. These results confirm an astronomical signal and, therefore, climatic forcing of relative sea level during OAE2 in the WIB. Furthermore, cross-basin correlation of the ATS and expanded δ13C chemostratigraphy of SH#1 suggests that these transgressive-regressive parasequences modulated siliciclastic sediment delivery in the seaway and contributed to deposition of prominent rhythmically bedded CTB units across the WIB, including the Bridge Creek Limestone. The presented approach to analysis of these proximal offshore siliciclastic facies links early diagenetic influences on chemostratigraphy to astronomically modulated sequence stratigraphic horizons, and helps to resolve rates of paleobiologic and paleoenvironmental change during a significant Mesozoic carbon cycle perturbation.

Fluctuations in depositional conditions during the onset of severe climate events in Earth history predispose stratigraphic archives to hiatuses, often hindering complete reconstructions of paleoclimate events and their triggers. Several studies have proposed that a hiatus of unknown duration exists at the base of Oceanic Anoxic Event 2 (OAE2) in the North American Western Interior Basin at the base Turonian global boundary stratotype section and point (GSSP) in Pueblo, Colorado, which potentially influences integrated radioisotopic, biostratigraphic, and astrochronologic age models of the Cenomanian-Turonian boundary interval. To quantify the duration of this hiatus, refine the chronology of OAE2, and assess marine geochemical perturbations associated with the onset of the event, we present new 40Ar/39Ar dates from regional bentonites along with a new proximal-distal chemostratigraphic transect of the epeiric Western Interior Basin (WIB), including initial osmium isotope (Osi) and stable carbon isotope (δ13C) data. The new 40Ar/39Ar age determinations confirm and further constrain previous estimates of Cenomanian-Turonian boundary timing. Further, the regional chemostratigraphic synthesis demonstrates the conformity of the OAE2 successions correlated to Pueblo, shows that the duration of the lag between the onset of the Osi and δ13C excursions is ∼60 k.y., and thus constrains the magnitude of the pre-OAE2 hiatus in Pueblo to less than this value. The new astronomically tuned, conformable Osi record across the onset of OAE2 captures a geologically rapid onset of large igneous province volcanism, consistent with other records, such that the addition of CO2 to the ocean-atmosphere system may have driven changes in marine carbonate chemistry. Additionally, the refined chronostratigraphy of OAE2 and the Cenomanian-Turonian boundary in the central WIB improves correlation with other records, such as those in the Eagle Ford Group, Texas. The correlations highlight that discrepancies among OAE2 age models from globally distributed sections commonly stem from differing definitions of the event and uncertainties associated with astronomical tuning, in addition to stratigraphic preservation.

Early Cretaceous (145–100 Ma) rocks record a ∼5‰ negative shift in the sulfur isotope composition of marine sulfate, the largest shift observed over the past 130 m.y. Two hypotheses have been proposed to explain this shift: (1) massive evaporite deposition associated with rifting during opening of the South Atlantic, and (2) increased inputs of volcanically derived sulfur due to eruption of large igneous provinces. Each process produces a very different impact on marine sulfate concentrations, which in turn affects several biogeochemical phenomena that regulate the global carbon cycle and climate. Here we present sulfur isotope data from Resolution Guyot, Mid-Pacific Mountains (North Pacific Ocean), that track sympathetically with strontium isotope records through the ∼5‰ negative sulfur isotope shift. We employ a linked sulfur-strontium isotope mass-balance model to identify the mechanisms driving the sulfur isotope evolution of the Cretaceous ocean. The model only reproduces the coupled negative sulfur and strontium isotope shifts when both hydrothermal and weathering fluxes increase. Our results indicate that marine sulfate concentrations increased significantly during the negative sulfur isotope shift and that enhanced hydrothermal and weathering input fluxes to the ocean played a dominant role in regulating the marine sulfur cycle and CO 2 exchange in the atmosphere-ocean system during this interval of rapid biogeochemical change.

This study revises and improves the chronostratigraphic framework for late Turonian through early Campanian time based on work in the Western Interior U.S. and introduces new methods to better quantify uncertainties associated with the development of such time scales. Building on the unique attributes of the Western Interior Basin, which contains abundant volcanic ash beds and rhythmic strata interpreted to record orbital cycles, we integrate new radioisotopic data of improved accuracy with a recently published astrochronologic framework for the Niobrara Formation. New 40Ar/39Ar laser fusion ages corresponding to eight different ammonite biozones are determined by analysis of legacy samples, as well as newly collected material. These results are complemented by new U-Pb (zircon) chemical abrasion–isotope dilution–thermal ionization mass spectrometry ages from four biozones in the study interval. When combined with published radioisotopic data from the Cenomanian-Turonian boundary, paired 206Pb/238U and 40Ar/39Ar ages spanning Cenomanian to Campanian time support an astronomically calibrated Fish Canyon sanidine standard age of 28.201 Ma. Stage boundary ages are estimated via integration of new radioisotopic data with the floating astrochronology for the Niobrara Formation. The ages are determined by anchoring the long eccentricity bandpass from spectral analysis of the Niobrara Formation to radioisotopic ages with the lowest uncertainty proximal to the boundary, and adding or subtracting time by parsing the 405 k.y. cycles. The new stage boundary age determinations are: 89.75 ± 0.38 Ma for the Turonian-Coniacian, 86.49 ± 0.44 Ma for the Coniacian-Santonian, and 84.19 ± 0.38 Ma for the Santonian-Campanian boundary. The 2σ uncertainties for these estimates include systematic contributions from the radioisotopic measurements, astrochronologic methods, and geologic uncertainties (related to stratigraphic correlation and the presence of hiatuses). The latter geologic uncertainties have not been directly addressed in prior time scale studies and their determination was made possible by critical biostratigraphic observations. Each methodological approach employed in this study—new radioisotopic analysis, stratigraphic correlation, astrochronology, and ammonite and inoceramid biostratigraphy—was critical for achieving the final result.

We develop an intercalibrated astrochronologic and radioisotopic time scale for the Cenomanian-Turonian boundary (CTB) interval near the Global Stratotype Section and Point in Colorado, USA, where orbitally influenced rhythmic strata host bentonites that contain sanidine and zircon suitable for 40Ar/39Ar and U-Pb dating. Paired 40Ar/39Ar and U-Pb ages are determined from four bentonites that span the Vascoceras diartianum to Pseudaspidoceras flexuosum ammonite biozones, utilizing both newly collected material and legacy sanidine samples of J. Obradovich. Comparison of the 40Ar/39Ar and U-Pb results underscores the strengths and limitations of each system, and supports an astronomically calibrated Fish Canyon sanidine standard age of 28.201 Ma. The radioisotopic data and published astrochronology are employed to develop a new CTB time scale, using two statistical approaches: (1) a simple integration that yields a CTB age of 93.89 ± 0.14 Ma (2σ; total radioisotopic uncertainty), and (2) a Bayesian intercalibration that explicitly accounts for orbital time scale uncertainty, and yields a CTB age of 93.90 ± 0.15 Ma (95% credible interval; total radioisotopic and orbital time scale uncertainty). Both approaches firmly anchor the floating orbital time scale, and the Bayesian technique yields astronomically recalibrated radioisotopic ages for individual bentonites, with analytical uncertainties at the permil level of resolution, and total uncertainties below 2‰. Using our new results, the duration between the Cenomanian-Turonian and the Cretaceous-Paleogene boundaries is 27.94 ± 0.16 Ma, with an uncertainty of less than one-half of a long eccentricity cycle.

Molybdenum (Mo) isotopes in ancient sediments are promising recorders of global ocean paleoredox conditions. Organic-rich black shales can be used to reconstruct ancient ocean Mo isotope compositions if these sediments record the isotopic composition of contemporaneous seawater. Comparison of δ 98/95 Mo in two Devonian shale sequences of similar age, the New York Oatka Creek and Geneseo Formations, reveals that this assumption cannot be applied to all organic-rich shales. Although both sequences contain laminated intervals, elevated organic carbon, and enrichments of redox-sensitive metals, the mean δ 98/95 Mo differs systematically between the formations by ~0.59‰. Independent paleoredox indicators reveal that portions of the Oatka Creek Formation were deposited under pervasively euxinic (anoxic and sulfidic) conditions, whereas conditions during deposition of the Geneseo Formation were intermittently euxinic to suboxic (oxygen deficient but not sulfidic in the water column). We infer that reconstruction of ancient ocean δ 98/95 Mo from organic-rich shales requires independent verification of persistent local euxinia. With these considerations in mind, our data point to δ 98/95 Mo in the Devonian oceans ~0.6‰ lighter than in today9s oceans, consistent with expanded anoxia.

The Late Devonian mass extinction was unusually protracted and ecologically selective, with preferential diversity losses among reefbuilding organisms and tropical, shallow-water faunas in general. We have investigated the link between the extinction’s unique characteristics and changes in biogeochemical cycling through analyses of the δ 13 C and C:N:P atomic ratios of organic matter buried across the Kellwasser Horizons in western New York State. Each horizon is characterized by (1) a long-term, + 4‰–5‰ excursion in δ 13 C, ~3‰ of which occurs within the horizon, and (2) a dramatic increase in the burial ratios of C:N:P, from values of ~100:15:1 to an average of ~5000:170:1. On the basis of these results, we propose that (1) increased efficiency of biolimiting nutrient recycling, resulting from cyclic water column stratification and mixing, promoted eutrophication during Kellwasser deposition in New York, and (2) the isotope excursions represent the composite effect of long-term, global organic C burial, and local changes in photosynthetic C isotope fractionation related to nutrient availability. This eutrophication model forges a mechanistic link between proposed Late Devonian climatic cooling and the selective demise of taxa likely to have been narrowly adapted to oligotrophic conditions.

Upper Cretaceous limestone-shale couplets developed within the late transgressive stage of the Greenhorn cyclothem may be correlated from carbonate-dominated (basinal) sequences in central Kansas and Colorado westward to clastic cycles in southern Utah. Six such basinal couplets have been traced to corresponding upward-coarsening progradational cycles developed on the western margin of the Western Interior basin. In the central basin in Colorado and Kansas, these sedimentary cycles are represented by limestone-shale and marlstone-shale couplets ∼0.5-1.0 m in thickness. More calcareous parts of these couplets may be correlated westward into condensed, fossiliferous concretion and shell beds in proximal offshore lithofacies of Arizona and Utah. These concretion and shell beds are physically traceable farther landward (westward) into bioturbated, fossil-rich, transgressive lag deposits that bound 10- to 20-m-thick coarsening-upward progradational strand-plain deposits (parasequences) in southwestern Utah. Thus, the progra-dational phase of parasequence deposition correlates with accumulation of clay-rich sediment in the central basin, and the transgressive phase is characterized by reduced terrigenous input and deposition of carbonate-rich sediment. We consider Milankovitch-style orbital forcing of climate and tectonically induced fluctuations in rates of foredeep basin subsidence as possible forcing mechanisms for these basinwide events. Based on the widespread distribution of the limestone-shale couplets, as well as on estimated sedimentation rates and geochronology, it has been widely speculated that these carbonate cycles reflect Milankovitch cycles with periodicities on the order of 20 k.y. to 100 k.y. If so, then stratigraphic data suggest that orbital forcing of climate affected eustasy and/or sediment input and biogenic production in the Western Interior basin. Alternatively, thrusting events in the Sevier orogenic belt may have produced episodic changes in the rates of foredeep basin subsidence and consequent changes in base level, which could have controlled the deposition of the Greenhorn parasequences and carbonate cycles. In either case, correlation of these units demonstrates a consistent basinwide sedimentary response to high-frequency base level or sediment input changes in the Western Interior epicontinental basin.

A combination of biostratigraphic markers (ammonites, inoceramid bivalves) and carbon isotope excursions is employed to establish a high-resolution correlation between the middle to late Cenomanian successions of the Western Interior Basin (USA) and the Anglo-Paris Basin (southern UK). Sequences identified from sedimentologic criteria in the Pueblo succession and elsewhere in the Western Interior Basin are shown to coincide precisely with globally recognized sea-level events and were therefore under eustatic control. This evidence refutes arguments that Cenomanian sequences in the Western Interior Basin were formed by local tectonic events. The interaction of longer-term tectonic movements and more rapid eustatic change may have simply enhanced the amount of erosion associated with sequence boundaries. A crossplot of radiometric ages derived from North American bentonites against an orbitally tuned time scale developed in the Anglo-Paris Basin provides support for the argument that the sequences were controlled by the 405-k.y.-long eccentricity cycle. © 2008 The Geological Society of America.

Previous time scales for the Cenomanian-Turonian boundary (CTB) interval containing Oceanic Anoxic Event II (OAE II) vary by a factor of three. In this paper we present a new orbital time scale for the CTB stratotype established independently of radiometric, biostratigraphic, or geochemical data sets, update revisions of CTB biostratigraphic zonation, and provide a new detailed carbon isotopic record for the CTB study interval. The orbital time scale allows an independent assessment of basal biozone ages relative to the new CTB date of 93.55 Ma (GTS04). The d 13 Corg data document the abrupt onset of OAE II, significant variability in d 13 Corg values, and values enriched to almost 222‰. These new data underscore the difficulty in defining OAE II termination. Using the new isotope curve and time scale, estimates of OAE II duration can be determined and exported to other sites based on integration of well-established chemostratigraphic and biostratigraphic datums. The new data will allow more accurate calculations of biogeochemical and paleobiologic rates across the CTB.

Many ancient rhythmic hemipelagic sequences have been interpreted to record orbital variations, but the exact nature of the climatic and depositional transfer functions responsible for this link remains poorly understood. Two-dimensional numerical simulations were used to explore selected aspects of orbital signal distortion in linked siliciclastic and hemipelagic systems. The models suggest that transfer of multiorder (e.g., 20, 100, and 400 k.y.) oscillations in relative sea level into the hemipelagic record produces an inherent amplitude distortion of the shorter-period (e.g., 20 k.y.) cycle. This distortion gives rise to amplitude modulation (AM), which is qualitatively similar to AM of orbitally driven changes in insolation (e.g., eccentricity modulation of precession-driven cycles). However, unlike the orbitally driven AM, synthesized AM is distinctly phase shifted relative to the stratigraphic record of the long-period (e.g., 100 k.y., 400 k.y.) cycle as a result of sealevel‐driven changes in the storage capacity of nearshore through alluvial parts of the source siliciclastic system. Hence, multiorder changes in sea level can leave a distinct AM signature in dilution-affected hemipelagic records, thus making hemipelagic rhythms due to eccentricity-forced sea-level changes distinguishable from other types of orbitally driven hemipelagic cyclicity.

The first systematic test for a predictive relationship between organic carbon content and stratigraphic hierarchy in a deep-water slope to basin-floor deposit was performed. The studied section includes the Pipeline Shale, the Brushy Canyon Formation, and the lower part of the Cherry Canyon Formation of the Delaware Mountain Group, West Texas. This interval represents one large-scale, 3rd-order genetic sequence within which 4th- and 5th-order stratigraphic cycles are recognized. Samples of fine-grained facies throughout the section were collected from outcrop and analyzed for organic carbon content and hydrogen index. Degree of pyritization was also determined for a subset of the samples. The results indicate that organic enrichment is closely correlated to the stratigraphic hierarchy at the 3rd-, 4th-, and 5th-order levels. The data suggest that quantity and quality of preserved organic matter are controlled by changes in bulk sedimentation rate (dilution vs. condensation), which affect organic matter inputs to the sediment, as well as the balance between (1) burial and preservation of organic matter and (2) its degradation on the sea floor during times of sediment starvation.

High-resolution sedimentologic, paleontologic, and stratigraphic data were collected for three Cenomanian-Turonian skeletal limestone packages in the Western Interior basin. Results of the study have significant implications for sequence stratigraphy in fine-grained basinal facies. (1) Skeletal limestones represent calcareous tempestites deposited at or near storm wave base, and their stratigraphic distribution provides a constraint on paleobathymetry. (2) Skeletal limestones are widely traceable and can be correlated to coeval progradational/retrogradational shoreface deposits, placing them in a relative sea-level context. (3) Skeletal limestone origin can be described by a comprehensive depositional model. The model incorporates effects of winnowing/erosion on relative fall, and sediment starvation/condensation on subsequent rise, relating both to the absence of an effective depositional gradient. (4) Application of the model suggests a reinterpretation of relative sea-level history for the Cenomanian-Turonian in the Western Interior that is better constrained than previous estimates.

To understand the patterns of lithofacies, marine faunas, organic-carbon enrichment, isotopes, and trace elements deposited in the early Turonian Western Interior seaway, we conducted circulation experiments using a three-dimensional, turbulent flow, coastal ocean model driven by GENESIS, a climate model developed at the National Center for Atmospheric Research (NCAR). Circulation and chemical evolution of the seaway waters are computed under the following initial and boundary conditions: (1) paleobathymetry according to a new interpretation of the lithostratigraphy and biostratigraphy; (2) temperatures and salinities of the Boreal and Tethys oceans and adjacent drainage basins based on isotopic data, atmospheric temperatures, and precipitation-evaporation magnitudes computed by GENESIS; and (3) mean annual wind stresses over the seaway computed by GENESIS. Results show that the seaway exported freshened water much like Hudson Bay today. Runoff from eastern drainages exited the seaway as a northern coastal jet; runoff from western drainages exited as a southern coastal jet. Both jets simultaneously drew in surface Tethyan and Boreal waters, creating a strong counterclockwise gyre occupying the entire north-south extent of the seaway. The curious stratal and faunal variations of the early Turonian deposits arise from this gyre and its associated water masses.