Previously published anomalous whole-rock stable isotopic values from the Poison Strip Sandstone Member of the Cretaceous Cedar Mountain Formation (CMF) of eastern Utah are of uncertain origins. This study investigated the diagenetic history and the processes responsible for these anomalous data. Accordingly, we integrated photomicroscopic techniques including polarized light microscopy, epifluorescence and cathodoluminescence (CL) imaging, micromilling of stable isotope samples, and fluid-inclusion heating and freezing measurements to this end. The key observations involved the microscopic mapping of calcite cement stratigraphy using CL imaging to permit the analysis of stable isotopes of calcite cements that crystallized during early and late diagenesis. The mapping of calcite cement zones of sufficient submillimeter size to mill out and isolate microgram-sized stable isotope samples enabled this discrimination. Early diagenetic calcite cements have the most positive δ18O values (-10 to -8.5‰ Vienna Pee Dee Belemnite [VPDB]) of all components. The pattern of δ13C and δ18O variation in this early diagenetic cement indicates affinities with early meteoric diagenesis previously documented in published literature on the CFM. The late diagenetic calcite cements yield the most negative δ18O values (-18 to -16‰ VPDB). We interpret the late diagenetic cements to be responsible for the anomalously low whole-rock δ18O values previously reported from the Poison Strip Sandstone Member. Our discoveries of bitumen in late-stage pore fillings and liquid petroleum in the fluid inclusions of late diagenetic calcite cements of the Poison Strip Sandstone Member explain the lower whole-rock organic matter δ13C values and anomalous Δ13C values reported from the unit. Comparatively lower carbonate δ18O and organic δ13C values originally derived from whole-rock analyses of samples from the Poison Strip Sandstone Member resulted from high-temperature basinal diagenesis (hydrothermal circulation and/or petroleum migration), rather than the alternative interpretation of early diagenesis related to a Cretaceous paleoclimatic perturbation. Our results are illustrative of methods to resolve the long-standing geologic problem of discriminating and characterizing products of early vs. late diagenesis in terrigenous clastic sedimentary strata.

The Cretaceous Cedar Mountain Formation (CMF) of eastern Utah is a widely recognized focus of research activity as an archive of the terrestrial record of early Cretaceous global change (Ludvigson et al., 2010, 2015; C. Suarez et al., 2014, 2020; M. Suarez et al., 2020). Studies of the CMF (Berriasian–Albian) on the east side of the San Rafael Swell have provided important δ13C and δ18O chemostrati-graphic data, U–Pb radiometric dates, vertebrate faunas, and palynomorph assemblages. It has also provided the means for interpreting Cretaceous terrestrial paleoclimates, depositional systems, and paleoenvironments in present western North America (e.g., Kirkland et al., 1997, 2016a, 2016b; Ludvigson el al., 2010, 2015; C. Suarez et al., 2014, 2020; Joeckel et al., 2017, 2019; M. Suarez et al., 2017, 2020). Although stable isotopic data from multiple mineral substrates in the CMF have previously been published, this paper is the first to resolve multiple diagenetic effects through the use of calcite cement stratigraphy, thus placing those effects in a temporal context.

Anomalously low whole-rock carbonate δ18O values of -16.0 to -13.5‰relative to the Vienna Pee Dee Belemnite (VPDB) isotope standard from sandstones of the Poison Strip Sandstone Member of the CMF at the Ruby Ranch Road locality (Fig. 1; see Ludvigson et al., 2015, their fig. 10) are markedly different from two distinct subsets of calcite δ18O values previously published from the CMF: (1) calcite cement δ18O values of approximately -8.0‰ VPDB are relatable to early meteoric diagenesis; and (2) vein-calcite δ18O values of -11.0 to -9.0 VPDB likely formed during burial (Ludvigson et al., 2010, 2015). Another important paleoclimate indicator, the Δ13C paleo-pCO2 proxy (carbonate δ13C values - organic δ13C values) yields problematic results in the Poison Strip Sandstone Member (Fig. 1; hereafter referred to as the ‘Poison Strip Member’). The Δ13C values from sedimentary rocks are a critical parameter in transfer functions to calculate ancient pCO2 values. We here consider the two alternative ideas posed by Ludvigson et al. (2015) to explain the origins of the anomalously low whole-rock carbonate δ18O values and organic δ13C values from the Poison Strip Member.

Depositional or Late Diagenetic Controls on Whole-Rock Oxygen Isotope Values?

One idea is that the whole-rock carbonate δ18O data from the Poison Strip Member result from authigenic processes during an Early Cretaceous paleoclimatic perturbation, such as a deglaciation in the nearby Sevier orogenic highlands (see Ludvigson et al., 2015, their fig. 3). The anomalously low δ18O values in the Poison Strip Member are similar to those reported from well-preserved Early and Late Cretaceous unionid bivalves that inhabited ancient rivers presumed to have carried some ancient alpine snowmelt (Glancy et al., 1993; Dettman and Lohmann, 2000). An alternative idea is that that the anomalous δ18O values result from much later, higher-temperature basinal diagenetic processes under burial conditions, directly related to either late-stage hydrothermal circulation and/or petroleum migration. This concept relates to the well-documented observation that burial diagenetic processes yield much lower carbonate δ18O values than those of earlier authigenic carbonates (Choquette and James, 1987; Algeo et al., 1992; Hasiuk et al., 2016).

Depositional or Late Diagenetic Controls on Whole-Rock Organic Carbon Isotope Values?

The lower organic matter δ13C values in the coarse-grained Poison Strip Member could possibly be related to the particulate organic matter inputs from C3 terrestrial vegetation of alpine paleoenvironments from higher elevations in the drainage basin. This organic matter could have differed from the more water-stressed lowland paleofloras characteristic of the remainder of the CMF. Lower δ13C values in C3 plants are associated with lower temperatures and higher relative humidities, both of which are compatible with alpine environments (Edwards et al., 2000). Water-stressed lowland paleofloras of the CMF (Ludvigson et al., 2015; Joeckel et al., 2017) could possibly have included crassulacean acid metabolism (CAM) plants with higher δ13C values than C3 plants, considering speculations about the appearance of CAM during the Mesozoic (Ehleringer and Monson, 1993). The δ13C values of authigenic carbonates throughout the CMF are interpreted to have been controlled by changes in paleoatmospheric δ13C values (Ludvigson et al., 2010, 2015; M. Suarez et al., 2020). Therefore, changes in the type of organic matter and the lower organic δ13C values in the Poison Strip Member would give rise to the anomalous Δ13C values in the unit. An alternative idea is that organic matter with lower δ13C values might have been introduced into the Poison Strip Member during later diagenetic processes under burial conditions. For example, migrated petroleum could alter the original depositional whole-rock Δ13C values in the Poison Strip Member.

The scientific results presented here address the lingering questions posed by the published whole-rock chemostratigraphic profile of the CMF at the Ruby Ranch Road section. An important goal of presenting these research methods and results is to provide an example of the framework for deconvolving early diagenesis from later alteration, and differentiating these products and processes in ancient terrestrial stratigraphic successions.

The Cedar Mountain Formation (CMF) is a distinct interval of continental sedimentary strata divided into six members in eastern Utah and westernmost Colorado (Kirkland et al., 2016b). The Poison Strip Member is underlain by the Yellow Cat Member—the basal member of the CMF that only occurs between the Green and Colorado rivers in Utah (Kirkland et al., 2016b; Joeckel et al., 2017, 2019). The Poison Strip Member is overlain by the more widespread Ruby Ranch Member across a large part of the Colorado Plateau (Kirkland et al., 2016b). The Poison Strip Member is a 2- to 15-m-thick stratigraphic package of very fine- to very coarse-grained sandstone and matrix-supported granule to pebble conglomerates with framework grain populations dominated by quartz and having minor feldspar and lithic fragments (Stikes, 2007). The Poison Strip Member is the amalgamated channel-fill and bar deposits of major fluvial systems that flowed northeastward to eastward from the Sevier thrust front and Mogollon Highlands across the Sevier foreland (Jordan 1981; Lawton, 1994; Lawton et al., 2010; Currie, 2002; DeCelles and Coogan, 2006; Stikes, 2007; Dickinson and Gehrels, 2008; Kirkland et al., 2016b). Lacustrine and stromatolitic limestones may be interbedded with, or assume a laterally equivalent position as the fluvial units (Kirkland et al., 2016b). In contrast, both the underlying Yellow Cat and overlying Ruby Ranch Members are dominated by drab to variegated mudrock, which is generally interpreted as floodplain and lacustrine sediments (Kirkland et al., 2016b). The Ruby Ranch Member generally ranges from 15 to 30 m in thickness, although it thickens locally and to the west across the San Rafael Swell, and it is dominated by variegated mudstone and shale with prominent strata of nodular carbonates or “calcretes,” for which multiple lines of evidence suggest a palustrine origin (Ludvigson et al., 2010, 2015).

This study assesses three samples from the Poison Strip Member in great detail. Two of the specimens come from an outcrop, which is the Ruby Ranch Road (RRR) locality of Ludvigson et al. (2010, 2015). These samples (RRR-12 and RRR-15) are both medium-grained sandstones. The third sample is a pebble conglomerate from a Utah Geological Survey (UGS) core drilled at Green River, Utah. This core sample is UGS-GRN-184 at the 178-foot-depth level (Figs. 23). Several serial cuts were made from this rock sample, and one of them, identified herein as UGS-GRN-184a, was used for stable isotope microsampling. All three rock samples used in this study came from a large research collection directly associated with earlier publications (Ludvigson et al., 2010, 2015). They were selected to address specific questions posed in these prior publications.

Petrographic, diagenetic, stable isotopic, and fluid-inclusion analyses of calcite cements in these samples discriminate the early and late temporal aspects of isotopic variability and crystallization temperatures by establishing the “cement stratigraphy” (Meyers and Lohmann, 1985; Meyers, 1991). Precise petrographic mapping and isotopic microsampling, which avoided inadvertent cross-contamination of cement types, were prerequisites for these analyses.


Rock samples were slabbed and vacuum-impregnated with blue-dyed epoxy resins. Epoxy-impregnated slabs were then cut and micropolished. Counterpart micropolished thin (0.03 mm) and thick (1 mm) sections were used for microsampling, following the techniques of Meyers and Lohmann (1985). Petrographic assessments of thin sections were carried out using an Olympus BX51 petrographic microscope with a digital 10.5 megapixel (Mpx) Olympus SC100 color camera. Images were collected under paired plane-polarized and cross-polarized illumination. Epifluorescence images were collected using an Olympus BX53M microscope with a digital 18 Mpx Olympus SC180 color camera. A cathodoluminescence (CL) imaging system allowed for the discrimination of cement types and zonation. A Relion Industries RELIOTRON III cold-cathode chamber—with operating conditions consisting of a rarified helium atmosphere at 50 millitorr, accelerating voltage of 10 kV, and beam current of 0.5 mA—was used for imaging. Macroscale imaging through the 50 mm top window of the chamber was carried out using a 16 Mpx Canon EOS SL1 digital single-lens reflex camera with a macro lens suspended over the CL chamber. Microscale imaging was carried out with the CL chamber mounted on a modified Olympus BX41 compound microscope using a Peltiercooled 17 Mpx Olympus DP73 color camera. Imaging modes included transmitted light (TL), reflected light (RL), and CL. Macroscale imaging was used to quickly scan 2 × 3 inch and 1 × 2 inch thin and thick sections for the identification of microscopic domains of special interest. Paired microscopic RL and CL images were also collected to map microsampling targets from micropolished thick sections at sub-millimeter scales. All petrographic, fluorescence, and CL imaging was performed at the Kansas Geological Survey.

Intragranular volumes (IGVs) of sandstone samples RRR-12 and RRR-15 were calculated from CL images of luminescent calcite cements to estimate approximate burial depths of cementation using ImageJ, a Java-based image processing program developed at the National Institutes of Health and the University of Wisconsin’s Laboratory for Optical and Computational Instrumentation.

Carbonate Microsampling and Stable Isotope Analyses

Carbonate microsampling for stable isotopes was guided by: (1) digital maps of micropolished thick sections acquired with a flatbed scanner; and (2) microscopic RL and CL images that were registered in graphics software to overlap perfectly. A high-precision dental drill with a 0.3 mm-diameter carbide bit was used to extract sample powders (0.05–0.1 mg/sample) from the polished thick sections (Fig. 3). Sample sites were selected and mapped on the matching CL and RL images of polished thick sections in advance to repetitively mill out the calcite components in sub-millimeter space. Thus, we were able to characterize the carbon and oxygen isotopic heterogeneity within each sequential cement zone. Sample powders were roasted in vacuum for 1 hour at 200°C to remove volatile contaminants prior to analysis. The roasted powders were then loaded into single reaction vessels in an automated Kiel III acid-dosing device and reacted with anhydrous phosphoric acid to produce a cryogenically-purified CO2 gas that was carried to a dual-inlet Thermo Finnigan MAT 253 stable isotope mass spectrometer (University of Kansas W.M. Keck-National Science Foundation Paleoenvironmental and Environmental Stable Isotope Laboratory). The subsequent isotopic results were calibrated using the Vienna Pee Dee Belemnite (VPDB) carbonate isotope standard (as in Ludvigson et al., 2015). The uncertainties of the δ13C and δ18O measurements are better than ± 0.1‰.

Fluid-Inclusion Analysis

Two doubly polished thick sections (UGS-GRN-184a and UGS-GRN-184b) were prepared from serial cuts of cored rock sample UGS-GRN-184 to enlarge the surface area available for inspection. Our petrographic observations identified both one-phase (liquid) and two-phase (liquid and vapor) fluid inclusions within calcite cements in doubly polished thick sections of the core sample. Consequently, we undertook fluid-inclusion microthermometry measurements using a Linkam THMSG600 stage and a Fluid Flow gas flow heating and cooling system. This system is mounted on an Olympus BX53M infrared microscope that is equipped with 4X, 10X, 40X, and 100X long focal distance objective lenses and with Link Interface Software. Fluid-inclusion micrographs were collected using a 2.7 Mpx FLIR Grasshopper® 3 USB 3.0 digital camera. Samples were examined to ensure full characterization of petrographic associations of fluid inclusions and to establish inclusion origins at the level of fluid-inclusion assemblages (Goldstein and Reynolds, 1994). Salinities are calculated from Tm measurements using equations from Bodnar (1992). Petroleum inclusions were identified using UV epifluorescence and the °API gravity of the oil-bearing inclusions as outlined by Goldstein and Reynolds (1994).


Core UGS-GRN-184 includes slightly more than 42 m (138 ft) of the Cedar Mountain Formation (CMF) (Fig. 3). There are no chronostratigraphic data from this core, but the distinct characteristics of three informal stratigraphic intervals within the core (Fig. 3: Intervals I, II, and III) provide a basis for correlations to regional outcrops. Interval I in the core (Fig. 3) is cross-stratified, very fine to medium sandstone and matrix-supported granule to pebble conglomerates. This interval directly overlies the Brushy Basin Member of the Jurassic Morrison Formation. The Yellow Cat Member—which overlies the Morrison Formation to the east of the city of Green River near Arches National Park, and which is overlain by the Poison Strip Member—is absent where the UGS-GRN-184 core was drilled (Kirkland et al., 2016b). Thus, both the stratigraphic position and the lithologic characteristics of Interval I lead us to correlate it as the Poison Strip Member, which is dominated by sandstone and pebbly sandstone (Kirkland, 2016b). Likewise, the variegated mudrock and the seven or eight interspersed nodular carbonate units (“calcretes”) comprising the overlying Interval II (Fig. 3) are wholly compatible with the characteristics of the Ruby Ranch Member in outcrops across the enclosing area (Kirkland et al., 2016b). Therefore, we correlate Interval II as a part of the Ruby Ranch Member. Interval III—another sandstone-dominated interval—is also likely to be a part of the Ruby Ranch Member, which contains sandstone bodies in its outcrop belt (Kirkland et al., 1997, 2016b).


Samples RRR-12 and RRR-15 are medium-grained sandstone with framework-grain populations dominated by quartz and chert. The optical continuity of 3- to 5-mm-diameter calcite cement crystals in cross-polarized light engulfs the 1- to 2-mm-diameter framework grains, resulting in a poikilotopic fabric. An average intergranular volume of 26.5% was calculated for the sandstone samples of RRR-12 and RRR-15 based on the area filled by the CL luminescent poikilotopic cement in thin-section view (Fig. 4). In addition to the poikilotopic calcite cement, another earlier zone of fine-grained isopachous cement crystals less than 10 microns in width coat the framework grains in samples RRR-12 and RRR-15 (Robertson, 2019). Samples from the Ruby Ranch Road section also contain late-phase, pore-filling bitumen that has undergone shrinkage on rock exposure (Fig. 5).

The poikilotopic cements in RRR-12 and RRR-15 are uniformly luminescent, and lack any variability in CL brightness or color (Figs. 6A–B). Blue-luminescing quartz framework grains, likely of high-temperature origin (Boggs and Krinsley, 2006), as well as a purple-violet luminescent microquartz that overprints the poikilotopic calcite cement, exist in both of the Ruby Ranch Road samples (Figs. 6A–B).

Sample UGS-GRN-184a is a pebble conglomerate near the base of the Poison Strip Member. It has large, sheltered, cm-scale intergranular voids infilled by calcite cements with CL zones of sufficient size (i.e., sampling targetzone diameters that exceed the diameter of the microsampling drill bit) to permit the mapping and microsampling of the full paragenetic sequence of cementation using a 0.3-mm dental burr (Figs. 7A–B). Cement development follows a paragenetic sequence of: (1) calcite cement Zone 1, which is nonluminescent dogtooth isopachous spar 0.5mm to 1 mm in width coating framework grains; (2) nonluminescent chalcedony cement, which replaces calcite cement Zone 1; and (3) calcite cement Zone 2, which is brightly luminescent calcite with equant crystal sizes ranging from 0.1 to 2 mm in size and infills any remaining pore space (Figs. 6C–D).The nonluminescent chalcedony cement preserves ghost structures of dogtooth crystal growth surfaces of calcite cement Zone 1 (Robertson, 2019). Calcite cement Zone 1 is present in all of our samples, but was only developed at a scale of sufficient size to be microsampled in UGS-GRN-184a. This paper seeks to understand the relationship between calcite cement zones 1 and 2 in sample UGS-GRN-184a, and the poikilotopic calcite cements in samples RRR-12 and RRR-15.

Stable Isotope Analyses

A total of 18 carbon and oxygen isotope microsamples were taken from the poikilotopic calcite cements in the micropolished slabs of RRR-12 and RRR-15.Tenmicrosamples came from calcite cement Zone 1 in the micropolished slab of rock sample UGS-GRN-184a. Seventeen additional microsamples of calcite cement Zone 2 were extracted from the same micropolished slab. Rock samples RRR-12 and RRR-15 yielded similar isotopic results. RRR-12 yielded δ18O values ranging from -18 to -15‰ VPDB, whereas RRR-15 yielded a slightly narrower range of -17 to -16‰ VPDB. The δ13C values ranged from -8 to -7‰ VPDB in RRR-12, and -7 to -6 in RRR-15 (Fig. 8). Calcite cements in both sandstone samples yielded relatively invariant carbon and oxygen isotopic compositions in comparison to the total spread of isotopic values (Fig. 8). These values are slightly lower than the oxygen isotopic range of -16 to -13.5‰ VPDB recorded by earlier studies of these same samples using whole-rock analysis (Ludvigson et al., 2015). Calcite cement zones 1 and 2 in the micropolished slab of rock sample UGS-GRN-184a yielded overlapping isotopic values. The initial nonluminescent dogtooth sparry calcite cement of Zone 1 yielded δ18O values ranging from -10.0 to -8.5‰ VPDB, and δ13C values ranging from -5.2 to -4.7‰ VPDB (Fig. 8). The pore-filling brightly luminescent sparry calcite of Zone 2 exhibited a wider range of isotopic compositions—δ18O values ranged from -17.0 to -8.5‰ VPDB, and δ13C values ranged from -6.5 to -5.0‰ VPDB (Fig. 8).

Fluid-Inclusion Analysis

Petrographic examination of calcite cements in the two doubly polished thick sections of the cored rock sample UGS-GRN-184 (Figs. 23) identified fluid inclusions that were the subject of further study. Calcite cements from two serial sections in the core contain one phase (liquid) and two-phase (liquid and vapor) fluid inclusions (Fig. 9). Homogenization temperature (Th) values range from 86°C to 115°C, and Tm values range from -0.2°C to -0.8°C (Table 1). Petroleum inclusions also were observed in calcite cements with °API gravity ranges estimated from 40 to 45 based on the color of UV fluorescence (Goldstein and Reynolds, 1994) (Figs. 10A–B). Several positive Tm values represent evidence of superheated ice under extreme negative pressure (Adams and Gibson, 1930; Table 1).

Cement Paragenesis

The paragenetic sequence of cementation in the Poison Strip Member is readily apparent from petrographic observations of the cement stratigraphy in micropolished thin section and slabs of rock sample UGS-GRN-184a (Fig. 11). The nonluminescent dogtooth calcite of calcite cement Zone 1 is the first phase of cementation, and it coated the largest framework grains in rock sample UGS-GRN-184a. Zone 1 provides the earliest record of diagenetic processes (Figs. 67) in the Poison Strip Member. Nonluminescent chalcedony cement precipitated after calcite cement Zone 1 (Fig. 7). This chalcedony cement replaced Zone 1 calcite cement, and it preserves ghost crystal growth structures of that predecessor cement. Brightly luminescent Zone 2 cement fills most of the remaining intergranular volume in the conglomerate, post-dating the earlier calcite cement Zone 1. The brightly luminescent Zone 2 calcite cement in micropolished thin sections and slabs of rock sample UGS-GRN-184a partially overprints the nonluminescent chalcedony and calcite cement Zone 1 calcite as a later pore-filling phase, providing further evidence for its relatively younger age and its placement in the paragenetic sequence (Fig. 7).

The poikilotopic calcite cements in the RRR-12 and RRR-15 samples fill intergranular volumes that range from 23.5% to 26.5%. Intergranular volumes in medium-grained sandstones in that range are characteristic of cementation of rigid-grained sandstones after pore space has collapsed under lithostatic loading at deep burial conditions of 2 to 3 km (Paxton et al., 2002, their fig. 6).

Pore-filling bitumen (Fig. 5) is an irrefutable indicator that petroleum migrated through the Poison Strip Member during a late stage of diagenesis. This petroleum migration filled pore space remaining after the preceding stages of cementation, and terminated the precipitation of poikilotopic calcite cements in RRR-12 and RRR-15. We hypothesize that pore space that now exists between any remaining bitumen and the surrounding calcite cement resulted from volume reduction as organic volatiles were lost with exposure at Earth’s surface (Fig. 5). The origin and age of the bright purple-violet luminescent microquartz cements in RRR-12 and RRR-15 remain unclear. These silica cements overprint the brightly luminescent poikilotopic calcite cements, placing them at later stage in the paragenetic sequence (Figs. 6A–B). Furthermore, CL properties of the bright violet-colored silica cement are unlike those of the nonluminescent chalcedony cements observed in thin sections from rock sample UGS-GRN-184a. Therefore, the two silica cements from two different locations probably are not genetically related.

Diagenetic Analysis of Strata from the Early Aptian Ap10 Interval

The early Aptian Ap10 interval was defined by Herrle et al. (2004) from high-resolution chemostratigraphic profiles of Cretaceous C-isotope stratigraphy. It is characterized by a pronounced minimum in δ13C values occurring between higher δ13C values of the early Aptian Ap7 and the late Aptian Ap15 intervals of Herrle et al. (2004). Ludvigson et al. (2010, 2015) correlated chemostratigraphic profiles of the Ap10 interval between the Ruby Ranch Road (RRR) and the Price River (PRC) stratigraphic sections (Fig. 2). From the Price River stratigraphic section (Ludvigson et al., 2010, 2015), rock sample PRC-23 from the Ap10 interval correlates with RRR-12 and RRR-15 from the Ruby Ranch Road stratigraphic section. Diagenetic evaluations of carbon-oxygen isotope data from the RRR rock samples, along with the correlative interval in the UGS-GRN-184 core, provide the basis for interpreting the early-to-late diagenetic evolution of strata from the Ap10 interval in the Cretaceous foreland basin. The special utility of this cored conglomerate sample at Green River, Utah, is that the large sheltered intergranular voids filled by calcite cements preserve CL-mapped zones of sufficient thickness to define and microsample the cement stratigraphy (Meyers, 1991). This approach permits the interpretation of the mineral paragenesis and diagenetic history of the Ap10 interval (Fig. 7).

Ludvigson et al. (2010) presented a diagenetic evaluation of rock sample PRC-23 from the chemostratigraphic profile of the Price River stratigraphic section, contrasting results from both the authigenic micritic or microspar components as well as later diagenetic vein or intergranular spar components. By combining these data (Ludvigson et al., 2010) with results from stratigraphically equivalent rock samples from the Poison Strip Member at Ruby Ranch Road (RRR-12 and RRR-15) and the UGS-GRN-184 core (Fig. 8), we can clearly define carbon-oxygen isotope values that are characteristic of early meteoric diagenesis, enabling the discrimination of products of early from later burial diagenesis in the CMF.

Diagenetic Interpretations of Carbonate Component Stable Isotope Data

Interpretations of the carbon and oxygen isotope data from calcite cements in the Poison Strip Member are based on recognition of distinctive patterns in carbon and oxygen isotope space that are characteristic of individual carbonate components, and on the well-documented scientific literature on interpretations of carbonate diagenesis. Burial cements typically have highly variable δ18O values and comparatively invariant δ13C values (Choquette and James, 1987). This covariant isotope pattern is primarily the result of increasing temperatures of precipitating waters under burial conditions, and the retrograde solubility of calcite (Choquette and James, 1987; Algeo et al., 1992; Hasiuk et al., 2016). A markedly different covariant spatial pattern is characteristic of early meteoric diagenesis, where a vertical trend of highly variable δ13C values with relatively invariant δ18O values, known as the “meteoric calcite line” (Choquette and James, 1987; Lohmann, 1988; Hasiuk et al., 2016). The carbon and oxygen isotope pattern observed in burial diagenetic environments was verified by Hasiuk et al. (2016) in a study that examined 28 subsurface petroleum reservoirs. The carbonate isotopic data taken from the Poison Strip Member clearly displays highly variable δ18O values with a covariant carbon and oxygen isotope linear trend with positive slope, matching the signature typically associated with carbonates influenced by burial diagenetic environments (Fig. 8).

Calcite cement Zone 1—the first phase of pore-filling cement and the earliest record of diagenesis—contains the highest δ18O values (-10.0 to -8.5‰ VPDB). The pattern of carbon and oxygen isotopic variation in Zone 1 calcite converges with the meteoric calcite line documented from correlative rock sample PRC-23 (Ludvigson et al., 2010), suggesting origin from early meteoric diagenesis during initial lithification of the conglomerate. This range of δ18O values is also similar to the δ18O values identified from burial diagenetic cements and those from authigenic components of rock sample PRC-23 from the Ap10 interval at the Price River stratigraphic section (Ludvigson et al., 2010, 2015). Petrographic observations indicate that the brightly luminescent calcite cement of Zone 2 partly replaces the nonluminescent Zone 1 cement, and while this intimate relationship provides cogent evidence for a paragenetic relationship, it also poses a challenge for the isolation of cement zones during microsampling. Thus, our ambiguous δ18O values from calcite cement Zone 1 may possibly reflect inadvertent solid-phase mixing with material from calcite cement Zone 2 during sample extraction (Fig. 8). Although it remains a possibility that Zone 1 is solely an early meteoric diagenetic cement, the overall pattern of the data in carbon and oxygen isotope space could also suggest that Zone 1 may preserve a record of the earliest stages of burial-diagenetic cementation.

Calcite cement Zone 2 yielded δ18O values ranging from -17.0 to -8.5‰ VPDB. These values plot in the gap between the higher Zone 1 δ18O values and the much lower δ18O values observed in the poikilotopic calcite cements samples from RRR-12 and RRR-15 (Fig. 8). The comparatively large oxygen isotopic variation in this cement zone suggests precipitation over a longer time span with changing subsurface temperatures relative to the other calcite cements sampled from the Cretaceous Ap10 interval (Hasiuk et al., 2016).

Poikilotopic calcite cements in samples RRR-12 and RRR-15 yielded the lowest δ18O values of any components sampled from the CMF (Fig. 8). Therefore, we consider them to be responsible for the anomalously low δ18O values originally generated by whole-rock isotopic analysis of the Poison Strip Member (Fig. 1). The stable isotopic data obtained from microsamples of rock sample UGS-GRN-184a clarified the relative timing and diagenetic environments of the poikilotopic cements microsampled from slabs of the Ruby Ranch Road rock samples. The poikilotopic calcite cements from the Ruby Ranch Road section and the Zone 2 cements from micropolished slabs of the UGS-GRN-184a rock sample have overlapping stable isotope compositions and similar crystal morphologies and CL luminescence characteristics. They represent the same stage of burial diagenesis. Poikilotopic calcite cements represent a later stage of pore-filling in the paragenetic sequence of the Poison Strip Member, the only later stages being bitumen pore fills and violet-luminescent microquartz cement (Fig. 11).

Burial Depth and Temperature

The average intergranular volume of 26.5% filled by the poikilotopic cements of samples RRR-12 and RRR-15 (Fig. 6) indicates burial depths ranging from 2 to 3 km (Paxton et al., 2002). This estimate is logically consistent with the thickness of overlying Cretaceous and Cenozoic strata eroded from the adjacent Book Cliffs (Fisher et al., 1960). Applying geothermal gradients ranging between 25° and 45°C per km of depth (Blackett and Wakefield, 2004) to this proposed depth range estimates burial temperatures of 75° to 155°C during the precipitation of the poikilotopic calcites cements. This estimate overlaps with the fluid-inclusion filling temperatures ranging from 87°C to 115°C from the calcite cements in rock sample UGS-GRN-184a.

In order to contrast early depositional from later burial temperatures, the temperatures of early meteoric diagenesis during the deposition of the Poison Strip Member can be estimated through the method described by Suarez et al. (2009). This method employs a second-order polynomial regression (Ufnar et al., 2002; Suarez et al., 2009), based on a Cretaceous latitudinal temperature gradient (Spicer and Corfield, 1992) that was calculated on the basis of fossil-leaf physiognomy (Wolfe and Upchurch, 1987). The equation of Suarez et al. (2009) is:


where t is temperature (°C), and l is paleolatitude. Using a paleolatitude of 34°N specified for the CMF, the Suarez et al. (2009) equation results in a calculated mean annual depositional surface temperature of 23.4°C. This mean annual temperature is also an estimate of the shallow subsurface temperatures of the meteoric waters from which calcite crystallization took place during early diagenesis, because shallow groundwater usually takes on the mean annual temperature of a location (Lohmann, 1988).

Influence of Petroleum Migration

Pore-filling bitumen (Fig. 5), a late paragenetic stage (Fig. 11), indicates basinal petroleum migration. This observation is unsurprising, because the Cedar Mountain and correlative Burro Canyon Formations are important oil and natural gas reservoirs in the Paradox Basin (Whidden et al., 2014, their fig. 13). This study presents the first published petrographic evidence for the migration of petroleum and other fluids through the Poison Strip Member. It is possible that petroleum migration occurred after or coincided with the precipitation of poikilotopic calcite cements. We originally estimated burial temperatures of 75° to 155°C, and measured fluid-inclusion filling temperatures ranging from 86° to 115°C. These temperature ranges would place the Poison Strip Member within the oil generation window during the growth of the poikilotopic calcite cement (Eremenko and Tverdova, 1980). Therefore, it is plausible that the fluids associated with petroleum migration influenced the precipitation of the poikilotopic calcite cements, and that they were also responsible for the low δ18O values derived from the unit. The burial diagenetic pattern of carbonates with widely ranging δ18O values that are arrayed along linear trends with positive slopes in carbon-oxygen isotope space, as documented in this work, is a characteristic of ancient carbonates influenced by petroleum migration (Hasiuk et al., 2016; Ross et al., 2020).

The timing and origin of the purple-violet CL luminescent silica cements (Figs. 6A–B) are unknown. The color of these cements, however, is compatible with an origin as quartz polymorphs that crystallized from hydrothermal solutions (Götze and Zimmerle, 2000; Götze et al., 2001; Witke et al., 2004; Boggs and Krinsley, 2006). Thus, we speculate that purple-violet luminescent silica cements crystallized from hydrothermal fluids emanating from the La Sal Mountains, Oligocene laccolith remnants some 40 km distant that are part of a larger regional trend of igneous activity and that exhibit common evidence for hydrothermal alteration, including the crystallization of silica (Hunt, 1958; Nelson et al., 1992; Ross, 1998).

Fluid Inclusions

A fluid-inclusion study of the calcite cement in this investigation indicates Th values ranging from 86° to 115°C, compatible with but below the maximum burial temperatures estimated from intergranular volumes preserved in rock samples RRR-12 and RRR-15. The distribution of fluid inclusion results along with carbon and oxygen isotope data indicate that these fluids (secondary and pseudo-secondary inclusions) are associated with migration of inter-strata water and petroleum during later stages of diagenesis. These oil-bearing inclusions with light blue color in epifluorescence light indicate intense ultraviolet with higher °API gravity. These oil-filled inclusions seem to be filled with immature to mature petroleum that migrated in the later stage of diagenesis. This migration is most likely associated with the petroleum system in the Paradox Basin (Whidden et al., 2014).

These data further support the idea that the anomalous lower carbonate δ18O values associated with the Poison Strip Member of the CMF are not associated with low-temperature meteoric calcite line patterns in carbon-oxygen isotope space, but rather with an isotopic pattern associated with higher temperatures during burial diagenesis. Accordingly, the low whole-rock δ18O values are not products of early meteoric diagenesis, and the idea of an authigenic record of a Cretaceous paleoclimatic event in the Poison Strip Member is not supported. Likewise, the anomalous whole-rock Δ13C values reported from the Poison Strip Member are not the result of any possible entrainment of sedimentary organic carbon particles from Cretaceous alpine paleofloras, but rather the presence of pore-filling bitumen that resulted from petroleum migration in the unit under burial conditions.

The low δ18O values of whole-rock carbonate cements in the Poison Strip Sandstone Member, atypical of the Cedar Mountain Formation (CMF), are not representative of authigenic processes in the CMF, but rather are the result of deeper burial diagenesis. The precipitation of late-stage poikilotopic calcite cements in the Poison Strip Member was likely associated with petroleum migration in the Paradox Basin. Nearby Oligocene laccolith emplacement in the La Sal Mountains might have also played an additional role by providing hydrothermal fluids for late silica cementation. Additional geochemical studies are needed to further explore the burial and thermal evolution of the petroleum system in the Paradox Basin. Our results clearly show that sandstone of the Poison Strip Member served as preferential flow paths for higher-temperature fluids, processes that resulted in anomalous whole-rock isotopic values in the CMF.

Co-author James Kirkland thanks the late Richard Dayvault (U.S. Department of Energy) for securing the UGS-GRN-184 core for the Utah Geological Survey Core Repository. We thank Bruce Barnett at the University of Kansas W.M. Keck-National Science Foundation Paleoenvironmental and Environmental Stable Isotope Laboratory for supervision of the stable isotope analyses carried out for this work. We thank our Kansas Geological Survey colleague Franciszek ‘Franek’ Hasiuk for technical review of an earlier draft, and Utah Geological Survey colleagues Donald DeBlieux, Grant Willis, and Stephanie Carney for technical reviews. We thank peer reviewers B. Ronald Frost, Laura Crossey, Sally Potter-McIntyre, and Minyu Yang for helpful suggestions that improved our presentation. Appreciation is extended to Rocky Mountain Geology Managing Editor Brendon Orr for editorial assistance.

Scientific Editor: B. Ronald Frost

Gold Open Access: This article is published under the terms of the CC-BY 3.0 license.