Stratigraphic architecture of the Belly River Group (Campanian, Cretaceous) in the plains of southern Alberta: Revisions and updates to an existing model and implications for correlating dinosaur-rich strata

Abstract

The Upper Cretaceous (Campanian) Belly River Group (BRG) of southern Alberta has a complex internal stratigraphic architecture derived from differential geometries of its component formations that resulted from regionalized tectonic influences and shifting source areas. A full understanding of BRG architecture has been compromised heretofore by a limited understanding of subsurface data in southwestern- and southeastern-most Alberta. In this study outcrop exposures throughout southern Alberta are tied to reference well logs and subsurface cross-sections allowing a more precise understanding of BRG architecture and how it relates to well-known vertebrate fossil producing areas. Modifications to an existing stratigraphic model of the BRG show that the Oldman and the Dinosaur Park formations have reciprocal north-to-south wedge-shaped geometries and a diachronous contact that become prominently expressed south of Twp 12. The updated model also demonstrates that the Oldman Formation thickens stratigraphically up-section to the south, and that the Foremost-Oldman contact is, essentially, a datum across much of southern Alberta. Identification of the Oldman Formation in the subsurface remains based on its relatively high gamma-ray response in mudstone successions, but it is also recognized that many of its sandstones exhibit relatively low gamma-ray responses like those in underlying and overlying formations. Nomenclature and subdivisions of the Oldman Formation are revised to accommodate this updated understanding, and modifications are also made to the definition of the Judith River-Belly River discontinuity, a newly recognized surface that marks the onset of accommodation and eustatic rise in sea-level in the northern Western Interior Basin at ~76.3 Ma.

Introduction

The complex internal stratigraphic architecture of the Belly River Group (BRG) clastic wedge in the plains of southern Alberta derives from differential geometries of its component formations across its regional expression [1–9]. Until now, aspects of this complex architecture have been documented within limited stratigraphic intervals and geographic areas, creating a piece-meal understanding of the wedge, and leaving a “big picture” understanding of the wedge and its formations largely unrealized. This is unfortunate because the BRG is a significant source of Campanian-age oil, gas, and coal resources as well as dinosaur and other vertebrate fossils in the Western Interior Basin (WIB) of North America [6, 9–12]. Accordingly, accurate litho- and chronostratigraphic correlation within the BRG at numerous widely separated localities across the southern Alberta plains and beyond Alberta depends on a clear understanding of how the group’s formational geometries relate to one another as well as to the units that bound the wedge (Pakowki/Lea Park and Bearpaw formations). Furthermore, advances continue to be made in our understanding of the litho- and chronostratigraphy of adjacent correlative units in Saskatchewan and Montana [e.g., 9, 13, 14] making it increasingly important that the BRG be understood in terms of both its parts and the whole to facilitate accurate correlation and meaningful interpretations of its depositional history and preservation within the northern WIB.

Reconstruction of Belly River Group architecture across southern Alberta is now possible due to a combination of location-specific studies (outcrop and subsurface) that have been conducted during the past 80 years, and the relatively recent availability of a rich assemblage of geophysical wireline logs that include data from the BRG strata that occur at very shallow depths and behind casing [8]. This report utilizes 105 well logs (Fig 1; Table 1; S1–S4 Figs) selected during a review of thousands to update and refine recognition of formational contacts in the subsurface, and to document the extent and geometry of lithostratigraphic units and markers within the BRG. Fourteen reference well logs are tied to key outcrop sections (and fossiliferous locations) in southern Alberta and Saskatchewan (Fig 1) thus ensuring accurate calibration and correlation of stratigraphic cross-sections throughout the field area. Emphasis is placed on southeastern Alberta where very little subsurface work has been previously conducted and fossiliferous exposures of the BRG are abundant. Throughout this report, the following abbreviations are used for conciseness: BFm, Bearpaw Formation; BRG, Belly River Group; DPFm, Dinosaur Park Formation; FFm, Foremost Formation; Fm, Formation; IHS, inclined heterolithic strata; km, kilometers; m, meters; MRFm, Milk River Formation; MRS/MRs, Milk River shoulder; OFm, Oldman Formation; PFm, Pakowki Formation; Twp, Township; Rg, Range; S, Section; W, west of.

[Figure omitted. See PDF.]

Light green indicates Belly River Group exposures in southern Alberta and southwest Saskatchewan. Inset shows field area in western Canada. Township and range numbers along margins. Cross-sections (A–A’, B–B’, C–C’) indicated by solid lines. Well logs indicated by dots and numbers. Reference logs indicated by yellow dots and bold numbers. Outcrop locations indicated by red stars. Comrey sandstone zone study area and well log locations indicated southeast of Medicine Hat. Depositional edge of Dinosaur Park Fm indicated by dashed line southeast of Lethbridge. Abbreviations: DC, Devil’s Coulee; DPP, Dinosaur Provincial Park; H/RH, Hays/Rolling Hills; OF, Onefour; MB, Manyberries; ML, Muddy Lake; NF, Nesmo Farm; RN, Rapid Narrows; RC, Redcliff; SF, Suffield Pumping Station; ss, sandstone.

[Figure omitted. See PDF.]

Geologic context

The Belly River Group is an eastward-thinning paralic-to-nonmarine clastic wedge (Fig 2A) that was deposited as part of a 3^rd order regressive-transgressive cycle in the Western Interior Basin (WIB) from ~81 Ma to ~74 Ma during the mid-to-late Campanian [15]. As redefined by Jerzykiewicz and Norris [2] and subsequently confirmed by numerous researchers [3–7], the BRG is bounded below and above by marine shales of the Pakowki/Lea Park (Claggett equivalent in the USA) and Bearpaw formations, respectively. The geographic limits of the BRG have been described previously [3–8] and are briefly reviewed here. The BRG can be recognized as far north as Edmonton and as far west as the deformed belt where Bearpaw shales pinch out. North of Edmonton, the BRG is equivalent to the lower Wapiti Formation. West of Calgary and into the deformed belt, the BRG is equivalent to the lower Brazeau Group. In the southwest foothills deformed area—beyond the range of reliable subsurface data and thus the limits of our study—the BRG has been divided into the Connelly Creek, Lundbreck, and Drywood Creek formations [2]. To the east in central Saskatchewan, the BRG thins stratigraphically and eventually pinches out in the marine shales of the Lea Park Formation. South of the international border, Belly River Group strata are equivalent to the Two Medicine and Judith River formations in northwestern and northcentral Montana, respectively [9].

[Figure omitted. See PDF.]

A, Stratigraphic chart for the Belly River Group. Note the complex geometries of the Foremost, Oldman, and Dinosaur Park formations, and the step-like nature of the Dinosaur Park-Oldman formational contact in southeastern Alberta. B, Stratigraphic model for the Belly River Group showing shifting sediment supply directions and intensities through time. Abbreviations: Ma, Mega-annum; nw, northwest; se, southeast.

Throughout most of the plains, and in ascending order, the BRG consists of the Foremost, Oldman, and Dinosaur Park formations (Fig 2A) [1, 3–8]. The three formations record significant adjustments in provenance, sediment supply and distributional directions over millions of years (Fig 2B), with the Oldman Fm representing a wedge of sediment sourced largely from volcanic terranes in northwestern Montana, and the Foremost and Dinosaur Park formations sourced farther north from centers in the Cordillera of central British Columbia [1, 6, 7, 16].

The Foremost Fm is a coaly, marine-to-paralic unit that prograded eastward and marked the retreat from Alberta of marine environments (Pakowki/Lea Park Fm) from ~80–78 Ma [8]. The Oldman Fm is a predominantly nonmarine alluvial unit that records the maximum retreat of the WIS from Alberta and western Saskatchewan at ~78–77 Ma [1, 7, 11, 12]. It thickens considerably in southwestern Alberta where it replaces the Dinosaur Park Fm and becomes younger up-section. To the north, it is eventually replaced by the Foremost Fm. As the Oldman Fm loses distinction to the north (~Twp 40), the Foremost and Dinosaur Park formations amalgamate, but are difficult to differentiate. The Dinosaur Park Fm is a marine-paralic and nonmarine alluvial unit that records the transgression of the Bearpaw Sea across western Saskatchewan and southern Alberta from ~77–74 Ma [1, 7, 9, 11, 12, 14]. It pinches out stratigraphically to the southwest in the Lethbridge area and thins toward the southeast (Manyberries area) where it extends into northern Montana. There, it may be the stratigraphic equivalent of the Coal Ridge Member of the Judith River Fm [1, 14].

Discussions of the complex nomenclatural history of the Belly River Group and its lithostratigraphic subdivisions in Canada are available in a variety of sources [1, 3, 7–9, 17, 18] and are not presented or reviewed here.

Methods

Raster well-log access and use was graciously provided by Divestco Geoscience Ltd. (Calgary, Alberta) through their “Energisite” website and platform. All well logs employed herein are listed in Table 1, and are presented here in annotated form in Figs 3–8, 12–14, and S1–S4 Figs. Well log data in Table 1 and illustrated well logs include formational picks presented as depth-from-surface measurements in meters. Gamma-ray and “porosity” logs (density, neutron, and sonic) are required to accurately differentiate the Foremost, Oldman, and Dinosaur Park formations in the Belly River Group [1, 4–6, 8, 9], and gamma-density (or neutron) pairs were mostly utilized here. Well log locations are reported in Table 1 and log headings using the standard legal land location designation for Alberta (Legal Subdivision-Section-Township-Range). Unless otherwise indicated (e.g., wells and outcrop sections in southwestern Saskatchewan), locations occur in Alberta west of the 4^th Meridian.

[Figure omitted. See PDF.]

Inset shows location (red dot) in southern Alberta relative to other reference logs (circles). Gamma (red) and neutron-porosity log patterns are typical of the Belly River Group north of Twp 19 in southern Alberta. Note relative thinness of the Oldman Formation compared to more southerly logs (Figs 4 and 5). Abbreviations: LCZ, Lethbridge coal zone; HSS, Herronton sandstone zone; CZ, coal zone; TCZ, Taber coal zone; MRS, Milk River shoulder. Formational abbreviations as in text. Abbreviations in inset indicate locations of Calgary, Lethbridge, Medicine Hat, and Red Deer.

[Figure omitted. See PDF.]

Note extreme thickness of the Oldman Fm in this area. Also note the measured section of fossiliferous strata from Devil’s Coulee, west of Warner (see Fig 1). Measured section was stratigraphically correlated to well log #68 using field observations of DAE. Abbreviations: B, datable bentonite; c, claystone; cs, coarse sandstone; fs, fine sandstone; masl, meters above sea level; ms, medium sandstone; PFm, Pakowki Formation; si, siltstone.

[Figure omitted. See PDF.]

Note the thickness of the Oldman Fm. Also included are measured sections of fossiliferous strata from Onefour, Nesmo Farm, and South Manyberries Creek (southeast of Manyberries). All sections are stratigraphically correlated to one another and well logs using local marker beds and formational contacts observed by DAE. Abbreviations: CSS, Comrey Sandstone; LCZ, Lethbridge coal zone; MCZ, McKay coal zone; NT, Nesmo tuff.

[Figure omitted. See PDF.]

Also included are measured sections of fossiliferous strata from the Bow and South Saskatchewan rivers near Medicine Hat. All sections are correlated to one another and well logs using local marker beds, formational contacts, and formation thicknesses observed by DAE as well as Hamblin [4]. Note the Herronton sandstone at the base of the Oldman in the measured sections.

[Figure omitted. See PDF.]

Also included is a measured section of fossiliferous strata from the Red Deer River at Dinosaur Provincial Park. All sections are correlated to one another and well logs using local marker beds, formational contacts, and formation thicknesses [1, 7, 11, 12]. The five U-Pb radiometric ages from the Oldman, Dinosaur Park and Bearpaw formations at Dinosaur Provincial Park are from Ramezani et al. [11] and Eberth et al. [12]. Also note the sandy lower one-half and muddy upper one-half of the Dinosaur Park Fm.

[Figure omitted. See PDF.]

Also included is a measured section from Muddy Lake, Saskatchewan [20]. Note the westward stratigraphic thickening of the Foremost and Dinosaur Park formations relative to the Oldman Fm.

Most well logs used here employ the Milk River “shoulder” as a datum for stratigraphic correlation and the compilation of cross-sections. The shoulder reflects a widespread unconformity and a sharp up-section transition from coarser (mudstones and sandstones) to finer sediments (shales) at the base of the Pakowki Fm (in southern Alberta) and the top of the Alderson Member within the Lea Park Fm at more distal locations to the north and east [4, 19]. The MRs can be traced to the southeast into Montana [9] but loses distinction in southwestern Alberta, near the international border, due to thinning and stratigraphic pinch-out of the marine Pakowki Fm. In that area, the top of the Milk River Fm is used as a datum.

Fourteen reference logs (Figs 3–8) were utilized here based on their clarity and effectiveness in characterizing typical well-log signatures for the stratigraphic units discussed here in different geographic regions of southern Alberta. Three subsurface cross-sections were compiled covering the full range of expression of the BRG in southern Alberta (Fig 1). A one-well-per-township spacing was assessed as sufficient to illustrate the formation-scale stratigraphic patterns described here. In areas north of Twp 30, where the Oldman Fm thins to less than 20 m and Belly River stratigraphy becomes more monotonous, well-log spacing was increased to 2–5 townships or more. In contrast, a one-well-per-section spacing was employed to test hypotheses regarding the continuity of the “Comrey Sandstone” and “Comrey sandstone zone” in southeastern Alberta (Fig 1). Subsurface data from the deformed belt, which is characterized by structurally complex strata with repeated stratigraphic intervals and a different stratigraphic nomenclature [2], are not included in this study.

Ten measured outcrop locations are indicated in Fig 1 and listed in Table 2. They include vertebrate-fossil-rich strata reported in the literature and typically include at least one formational contact to facilitate correlation with geophysical data from nearby wells. Sections were measured using a 1.5 m Jacob’s staff with a top-mounted Abney level. As presented here, section graphics depict grain size and basic sedimentary structures. Colors in the graphics reflect increasing amounts of organic content in fine grained deposits (yellows-to-browns), iron-staining in sandstones (red-to-purple), and the presence of airfall bentonites (green).

[Figure omitted. See PDF.]

Characterizing the Belly River Group in well logs

The lower and upper boundaries of the BRG are marked by interfingering and thus diachronous transitions between marine and paralic facies that make precise placement of formational contacts somewhat arbitrary in any given location or well [17, 18, 21, 22]. For consistency and convenience in this report, the lower contact of the Foremost Fm on the Pakowki/Lea Park Fm marine shales is placed at the base of the first prominent upward-coarsening sandstone that also marks the base of a succession of interbedded shales, coals, siltstones, and sandstones that is ultimately overlain by the non-marine Oldman Fm (e.g., Fig 3). In this context, marine shales may be included within the formation but do not dominate. Exceptions to this pattern occur in southwestern Alberta near the international border. There, the Pakowki Fm thins and carbonaceous, paralic strata of the lower Foremost Fm approach or amalgamate with sandy paralic strata of the Milk River Fm (e.g., Fig 4). Likewise, the upper boundary of the BRG, which consists of the contact of Bearpaw Fm (marine shales) on the Dinosaur Park Fm (paralic sediments) is placed above the highest sandstones, siltstones, coals, and carbonaceous shales of the Dinosaur Park Fm, at a position where shales begin to dominate the overlying succession (Fig 3). Exceptions to this pattern occur west of Calgary where marine shales of the Bearpaw Fm pinch out and the Belly River Group transitions laterally into the lower portion of the Campanian-Maastrichtian age non-marine Brazeau Group.

Whereas picking the lower and upper bounding contacts of the BRG in the subsurface is a relatively straightforward process, picking the FFm-OFm and OFm-DPFm contacts requires greater consideration. Dowling’s [23] two-fold division of the BRG (Foremost and Pale beds) was further refined by Russell and Landes [17] who elevated those units to formational rank and renamed the Pale beds the Oldman Fm (based on “type” exposures along the Oldman River between Lethbridge and the downstream confluence with the Bow River). Since that time, the Foremost-Oldman formational contact in the southern Alberta plains has been placed in the middle of the BRG, at the top of the upper-most coal/carbonaceous shale succession in the Foremost Fm and/or at the base of the first fluvial sandstone in the overlying silty-sandy, non-coaly alluvial succession (Oldman Fm; e.g., [1, 4–6, 21]).

McLean [18] emphasized electric logs in his subsurface study of the Belly River Group (his Judith River Fm) but could not differentiate with any consistency the FFm-OFm formational contact in either outcrop or subsurface as described by previous authors. In contrast, Eberth and Hamblin [1] and subsequent workers [3–9] have emphasized gamma-ray logs in their study of the BRG and have recognized a consistent and significantly high gamma-ray response through muddy portions of the Oldman Fm (Fig 3; also visible in the subsurface cross-sections of MacDonald et al. [21]. Eberth and Hamblin [1] studied the region surrounding Dinosaur Provincial Park (Twp 20-Rg 12) in detail, formally recognized the uppermost ~80 m of the Oldman Fm as a new unit, the Dinosaur Park Fm, and used the relatively high gamma-ray response in the remaining Oldman Fm as a means to identify and differentiate that formation in the subsurface and to identify and map the OFm-DPFm contact, which they interpreted as a regional discontinuity (ODPD). Specifically, they recognized a high gamma-ray log signature in the Oldman Fm followed by an abrupt leftward gamma-ray deflection (drop) upward across the base of the Dinosaur Park Fm as a unique feature that defines the OFm-DPFm contact and discontinuity in well logs from that area. However, they did not fully assess the utility of the high gamma-ray response as a diagnostic feature of the Oldman Fm elsewhere in southern Alberta. They also did not trace the OFm-DPFm contact and discontinuity in the subsurface to the southeast or southwest corners of the province, where the Oldman Fm had been observed to thicken considerably in outcrop [1, 17].

Hamblin [4, 5] and Hamblin and Lee [6] also recognized the Dinosaur Park Fm, and redefined the Oldman Fm in the subsurface as a two-fold unit consisting of a lower interval of sandstone(s) (their “Comrey Member sandstone”) overlying the Taber Coal zone (and its equivalents), and an upper interval of mudstones with a relatively high gamma-ray response (their “upper siltstone member”). Hamblin [4] presented a gamma-density log pair from a well near Gleichen (see our Fig 3 as an example of a log pair from that area) as characteristic of the Oldman Fm’s subsurface features across southern Alberta and north to Twp 45 and utilized gamma-ray-porosity log pairs to identify and trace the formation in the subsurface. Hamblin [4] showed the Oldman Fm having a thickness of about 50 m or less north of Twp 22 but provided no subsurface documentation south of Twp 19 (his Figs 3–6). Subsequent studies have relied on the data provided by Eberth and Hamblin [1] and Hamblin [4] to define the FFm-OFm and OFm-DPFm contacts [e.g., 7–9].

Whereas the formational picks and interpretations of thickness for the Oldman Fm as presented by Eberth and Hamblin [1] and Hamblin [4] are accepted here for areas of Alberta north of Twp 19, it is clear based on data gathered in this report that the subsurface characterization of the Oldman Fm in southernmost Alberta (south of Twp 19) as described by those authors remains limited and requires modification. Specifically, previous studies neither describe nor depict the subsurface (well log) features that characterize the greatly thickened Oldman Fm in southeastern and southwestern Alberta. The absence of an adequate subsurface characterization for the Oldman Fm in those areas has important implications for delineating the geometries and correlations of the formations that make up the BRG across southern Alberta and beyond the limits of the province.

Characterizing the Oldman Formation in the subsurface

Although the Oldman Fm is easily differentiated in outcrop from overlying and underlying formations by a large variety of sedimentological data [e.g., 1, 17], the unit has been difficult to distinguish from other formations using geophysical data alone [18]. An examination of well logs from southern Alberta in this study, specifically the southeastern and southwestern corners of the province, confirms that the most consistent feature that characterizes the Oldman Fm in the subsurface—and thus differentiates the three formations of the Belly River Group—is a relatively high gamma-ray response in Oldman Fm mudstones (an increase of ~10–30 API) compared to mudstones in the immediately underlying portions of the Foremost and overlying Dinosaur Park formations. Although relatively high gamma-ray responses may occur within single mudstone beds in the non-marine portions of the Foremost and Dinosaur Park formations, multi-meter successions of mudstones with similarly high gamma-ray responses remain absent. However, as the Foremost and Dinosaur Park formations transition down- and up-section, respectively, into marine shales of the Pakowki and Bearpaw formations, extensive intervals with high gamma-ray responses are re-established (Fig 3).

Well log data from six reference wells (Figs 3–5), many of which are tied to nearby exposures, show that as the Oldman Fm thickens toward the international border in both the southwest and southeast regions of the province, a high gamma-ray response is maintained throughout the formation. In each well log there are multiple intervals of mudstones with gamma-ray responses that are higher overall than those in the Foremost and Dinosaur Park formations. However, this pattern is made complex by the presence of multiple sandstone bodies throughout the Oldman Fm, many of which that show low gamma-ray responses like those in the Foremost and Dinosaur Park formations (Figs 4 and 5). These observations underscore that the relatively higher gamma-ray response that characterizes the Oldman Fm is mostly limited to the formation’s mudstones, and that gamma-ray comparisons between the Foremost, Oldman and Dinosaur Park formations should, accordingly, focus on their mudstone units. Unfortunately, this point was only briefly touched on by Eberth and Hamblin [1, p. 183], and has subsequently been overlooked by subsequent researchers who have based their understanding of the Oldman-Dinosaur Park formational contact in the subsurface on well logs from north of Twp 19. Those logs show a clear dichotomy between a single interval of high-gamma mudstones in the Oldman Fm overlain by a prominent low gamma-ray response in sandstones that characterize the base of the Dinosaur Park Fm (Fig 3).

Whereas the Oldman Fm in subsurface can be characterized by the higher gamma-ray response of its mudstones (and some of its sandstones) relative to those in the Foremost and Dinosaur Park formations, its basal contact on the Foremost Fm is typically marked by an upward-fining (fluvial) sandstone(s) that marks the transition from carbonaceous mudstones, coals, and/or shallow-marine sandstones of the Foremost Fm to the high gamma-ray-response alluvial mudstones of the Oldman Fm [4] (Figs 3, 5–7). Similarly, following the convention first proposed by Eberth and Hamblin [1, p. 183], the Oldman-Dinosaur Park formational contact is placed at the base of the first sandstone that fines upward into one or more mudstone successions in the Dinosaur Park Fm with gamma-ray responses lower than those in the underlying Oldman Fm. As described by Hamblin [4], the contact with the overlying Dinosaur Park Fm is easily picked north of Twp 19 where the Oldman Fm is quite thin (~50 m thick, Fig 3). In the southeastern and southwestern corners of the province the Oldman Fm thickens considerably and placing the Oldman-Dinosaur Park formational contact requires consideration of the gamma-ray response within multiple mudstone successions between the Taber and Lethbridge coal zones.

Figs 6–8 show additional reference well logs from the southeast, central and northern regions of our field area that characterize variations in the gamma-ray response, thickness, and geometry of the three formations that make up the BRG. Measured section data from nearby outcrops that exhibit formational contacts are included and allow formational picks in well logs to be ‘ground-truthed’. Specifically, Fig 6 includes two well logs in southcentral and southeast Alberta that can be directly correlated to outcrops on the Bow and South Saskatchewan rivers that exhibit the FFm-OFm and OFm-DPFm contacts. Fig 7 includes a well log from Dinosaur Provincial Park that can be correlated to the type outcrop section of the Dinosaur Park Fm near Iddesleigh (along the Red Deer River) and that exhibits the OFm-DPFm and DPFm-BFm contacts. Lastly, Fig 8 includes two well logs from the northeastern portion of our field area that can be correlated with fossiliferous exposures of the Oldman and Dinosaur Park formations at Muddy Lake, Saskatchewan (Twp 39, Rg 23W3), and that exhibit the FFm-OFm and OFm-DPFm contacts.

Regional cross-sections

Regional cross-sections were assembled using 85 gamma-porosity log pairs (Table 1, S1–S3 Figs). These delineate Belly River Group formational geometries across southern Alberta (Figs 9–11). Although each cross-section is based on data from tens of well logs, individual logs are not shown in the cross-sections due to image-resolution constraints. Instead, the cross-sections exhibit color-coded formational picks and thicknesses at each well location, which are projected between locations.

[Figure omitted. See PDF.]

Cross-section A-A’ is based on data from 35 wells logs (Fig 1, Table 1, S1 Fig). Logs are not illustrated in the cross-section due to image-size constraints. Instead, the cross-section exhibits color-coded formational thicknesses within and between adjacent wells. Numbers relate to well locations shown in Fig 1. Numbers in red indicate reference logs shown in detail in Figs 5–8. Note the southward thickening (mostly up-section) of the Oldman Formation and the commensurate thinning of the Dinosaur Park Fm in southeastern-most, as originally hypothesized by Eberth and Hamblin [1]. The Foremost-Oldman formational contact rises slightly to the north due to stratigraphic thinning of the Oldman Fm. Also note the occurrence of a well-developed two-fold stratigraphic division of the Oldman Fm in the southern-most corner of the field area (Twp 01–05). This pattern of a lower sandstone dominated interval overlain by a mudstone dominated interval compares closely to that described by Rogers et al. [9, 14] for the Judith River Formation. The contact between the two intervals in southeast Alberta thus likely marks the mid-Judith discontinuity (MJD) described by those authors.

[Figure omitted. See PDF.]

Cross-section B-B’ is based on data from 28 well logs (Fig 1, Table 1, S2 Fig) and shows the mostly tabular geometry of the Oldman Fm that is maintained W–E (proximal-distal) at Twp 21. In contrast, both the Foremost and Dinosaur Park formations exhibit stratigraphic thickening to the west (toward source). Thus, the mostly tabular W-E geometry of the Oldman Fm confirms its separate and more southwesterly source area. Note also that west of Rg 28 (Calgary) the Bearpaw Fm pinches out and the Belly River Group is equivalent to the lower portion of the much thicker non-marine Brazeau Group.

[Figure omitted. See PDF.]

Cross-section C-C’ is based on 24 well logs (Fig 1, Table 1, S3 Fig) and shows the extreme thickening of the Oldman Fm to the southwest, which underscores the presence of its unique, southwesterly source area. Note the southward termination (depositional edge) of the Dinosaur Park Fm near Lethbridge (Fig 1). Southward across the international border in this region, the Belly River Group correlates with the Two Medicine Formation.

Cross-section A–A’ (Figs 1 and 9) runs ~400 km north-south along the Alberta-Saskatchewan border from Provost to the international border. Based on 35 well logs (Fig 1, Table 1, S1 Fig), it provides a strike section through the more distal portion of the Belly River Group clastic wedge at the eastern edge of our field area. From north to south the Oldman Fm mostly thickens stratigraphically up-section. At its most northern extent (Twp 40) the Oldman Fm is essentially a 10-m-thick marker bed that continues to thin and become unrecognizable farther north. At Twp 36, it is 25–35 m thick and maintains that thickness south to Twp 23 where it thickens sharply to more than 50 m. It maintains a thickness of 50–70 m farther south to Twp 11, where it thickens sharply to more than 80 m. From Twp 11 south to the international border the Oldman continues to thicken gradually, reaching a maximum of 126 m at Twp 02. From Twp 05 to the international border there is evidence for a two-fold subdivision of the OFm into a lower sandstone dominated interval and an upper mudstone dominated interval.

In contrast to the Oldman Fm, both the Foremost and Dinosaur Park formations in cross-section A–A’ exhibit minor stratigraphic variations in thickness. The Foremost Fm exhibits a total stratigraphic thickness of ~70–80 m at Twp 40, thickens to ~100 m at Twp 22, and then thins gradually to 50–60 m south of Twp 04, near the international border. The Dinosaur Park Fm exhibits a more consistent southward-thinning trend with a stratigraphic thickness of 50–70 m from Twp 40 south to Twp 07, and then thinning to 30–40 m from there to the international border. The FFm-OFm contact remains essentially flat relative to the MRs throughout its extent in cross-section A–A’, with a height above the MRs that ranges from 190–225 m. Its highest stratigraphic occurrence above the MRs is limited to the north (Twps 30–40) where it begins to pinch out, suggesting that the Foremost Fm thickens slightly up-section in that direction.

Cross-section B–B’ (Figs 1 and 10) runs ~230 km west-east in Twp 21 from Calgary to the Alberta-Saskatchewan border near Empress. It is based on 28 well logs (Fig 1, Table 1, S2 Fig), reveals a dip section through the Belly River clastic wedge in southern Alberta, and exhibits pronounced eastward-stratigraphic-thinning of both the Foremost and Dinosaur Park formations as the BRG wedge extends farther into the basin. Similarly, the FFm-OFm contact drops gradually from east to west relative to the MRs reflecting the overall proximal-distal thinning of the entire wedge. Conversely, the Oldman Fm maintains a relatively constant thickness (40–60 m) and shows a slight stratigraphic thickening (up to 70 m) only in its most eastern locations (Rg 06–01).

Cross-section C–C’ (Figs 1 and 11) extends for 340 km north-south from Big Valley to near the international border in southwestern Alberta. It is based on 24 well logs (Fig 1, Table 1, S3 Fig) and exhibits a strike section through a more proximal portion of the Belly River Group clastic wedge at the western edge of the field area. At Twp 35 and locations farther north, the Oldman Fm is a four-meter-thick marker bed that is difficult to identify in well logs; it eventually disappears north of Twp 45, marking the point where the Dinosaur Park and Foremost formations amalgamate and become indistinguishable. South of Twp 35 to the international border, the Oldman Fm thickens stratigraphically and mostly up-section. At Twp 10 it thickens sharply to ~160 m. At Twp 8 it thickens to more than 230 m, completely replaces the Dinosaur Park Fm, and is overlain by the Bearpaw Fm. From there south to the border it maintains a thickness of ~200 m whereas the Foremost Fm continues to expand down-section, eventually amalgamating with the underlying Milk River Fm at Twp 01. Like cross-section A–A’ the FFm-OFm contact remains flat relative to the MRs with an average height of ~200 m above the MRs. However, the contact can exhibit a large amount of variation in stratigraphic position south of Twp 08. The Dinosaur Park Fm pinches out sharply at Lethbridge (Twp 08–09), an observation previously made by Hamblin [5].

Discussion

Updating the Belly River Group stratigraphic model

Data presented here support and add detail to the longstanding Belly River stratigraphic model [1, 4–7]. The existing model consists of three distinct stages in sediment supply and distribution, and shoreline orientation through the Belly River Group (Fig 2B). Each stage relates exclusively to one of the three formations that make up the BRG. The Foremost Fm reflects the earliest stage and is characterized by widespread western-to-northwestern sediment sources in the British Columbia cordillera supplying extensive north-south oriented shorelines that prograded cyclically to the east and ultimately formed stacked successions of upward-coarsening parasequences across southern Alberta [3, 6, 24]. The Oldman Fm records a second stage of stratigraphic evolution and is a mostly non-coaly nonmarine unit sourced primarily from northwestern Montana and from limited regions of southwestern Alberta and southeastern British Columbia. Oldman Fm sediments prograded rapidly to the east and northeast across peat-forming wetlands of the uppermost Foremost Fm. Oldman Fm sediments ultimately prograded into central Saskatchewan where they recorded the maximum extent of the BRG wedge, as well as the maximum regression of the Pakowki/Lea Park sea [1, 4]. The third stage in the evolution of the BRG wedge records a variety of lithologic, mineralogic, sedimentary, geophysical, paleontological, and taphonomic differences between the Oldman and overlying Dinosaur Park formations, reflecting a significant reorganization in sediment supply and distribution across southern Alberta during the onset of transgression of the Bearpaw Sea (Fig 2B). Eberth and Hamblin [1] and others [4–7] interpreted the reorganization as due to the onset of tectonic rebound and unroofing in the proximal foredeep of the central cordillera (central British Columbia), an event that supplied a large volume of sediments and initially overwhelmed accommodation created during eustatic sea-level rise (Bearpaw transgression) and tectonic adjustments in southern Alberta. They concluded that Dinosaur Park Fm sediments expanded toward the southeast during transgression, gradually overstepping Oldman Fm wedge and expanding into northern-most Montana.

Data presented here support the existing three-stage tectono-stratigraphic model for the Belly River Group, but more clearly document the unique geometries and well-log signatures of the Oldman and Dinosaur Park formations that are important elements of the model. Specifically, the dramatic southwestward thickening of the Oldman Fm and overall higher gamma-ray response of the formation’s fine-grained sediments confirms the presence of a geographically unique source area in northwestern Montana and southwestern Alberta. The data also document the development of a widespread and broadly isochronous horizon of carbonaceous shales and coals at the top of the Foremost Fm (Taber coal zone), upon which the Oldman Fm alluvial sediments were deposited. Regarding the Dinosaur Park Fm, the depositional edge of the formation can now be projected as a nnw–sse trending feature across southern Alberta (Fig 1). The southward thinning of the Dinosaur Park Fm is also more clearly documented relative to the thickening of the Oldman Fm, and thus further confirms that these two units were sourced from separate regions of the cordillera throughout their histories.

Revising Oldman Formation lithostratigraphy and nomenclature

The data presented here allow us to clarify, update, and revise aspects of Oldman Fm stratigraphic nomenclature that have evolved during the past quarter-century. Eberth and Hamblin [1] modified the original definition of Oldman Fm [17] when they assigned the uppermost portion of the original OFm in the region surrounding Dinosaur Provincial Park to a new stratigraphic unit, the Dinosaur Park Fm, based on a variety of sedimentological, geophysical and paleontological criteria. Hamblin [4] then subdivided the Oldman Fm. He recognized two informal members making up the formation north of Twp 19: the lower “Comrey Member sandstone” and the “upper siltstone member”. In his Figs 4, 6, 9 and 11, he consistently shows the base of his “Comrey Member sandstone” as more or less coincident with the base of the Oldman Fm on, or up to 2 m above the Taber coal zone, and 210–220 m above the Milk River shoulder.

However, outcrop data from southeastern Alberta along the Milk River canyon area (Fig 5) show that the base of the type Comrey Sandstone—originally named for the abandoned Comrey settlement near Manyberries (Twp 2, Rg 6; [17, 25])—lies ~60 m or more above the Taber coal zone in the middle of the Oldman Fm, and ~260–265 m above the Milk River shoulder datum. Using these data, Eberth [7] applied the term “Comrey sandstone zone” to mineralogically mature sandstones in the middle of the Oldman Fm throughout southeastern Alberta, and proposed them to be stratigraphically equivalent to the type Comrey Sandstone exposed along the Milk River canyon, just south of Comrey.

Eberth [7] recognized the same sandstone zone recognized by Hamblin (1997a) on or just above the Taber coal zone throughout southern Alberta in logs and outcrop but referred it to the “Herronton sandstone zone” after the hamlet of Herronton, which lends its name to the productive, shallow-gas-field that produces mostly from this sandstone horizon southeast of Calgary. Eberth [7] assigned the Herronton sandstone zone to the top of the Foremost Fm due to its distinctive mineralogy (immature and bentonitic fine-to-medium-grained sandstones) and sedimentology (locally incised multistoried high-sinuosity paleochannels characterized by clean sandstone and inclined heterolithic stratification [IHS]). Data presented here from numerous locations along the Bow and South Saskatchewan Rivers and the Milk River canyon area of southeastern Alberta (Figs 5 and 6) indicate the presence of a well-developed multimeter thick sandstone zone, often consisting of IHS paleochannel deposits, on or just above the Taber coal zone, and thus equivalent to the Herronton sandstone zone of Eberth [7].

A comparison of the interpretations of both Hamblin [4, Fig 14] and Eberth [7, Fig 3.9], as well as a review of their sedimentological data, show that the “Comrey Member sandstone” of Hamblin [4] and “Herronton sandstone zone” of Eberth [7] refer to the same sandy gas-producing stratigraphic interval that rests on, or slightly above, the Taber coal zone. For reasons of nomenclatural clarity, and to maintain consistent reference to and accurate stratigraphic placement of the type Comrey Sandstone in southeastern Alberta as originally presented by Troke [25], we propose that the sandy, IHS-rich interval that is positioned at the base of the Oldman Fm on or just above the Taber coal zone be henceforth referred to as the “Herronton sandstone zone” as proposed by Eberth [7] and in line with industry nomenclatural practice. Although future work may reveal that the Herronton sandstones are mineralogically unlike those higher in the Oldman Fm (as suggested by Eberth [7] and supported by the data of Eberth and Hamblin [1, their Fig 6], it is proposed here that the Herronton sandstone zone be included, for the time being, at the base of the Oldman Fm. This follows lithostratigraphic and historical precedents set by most researchers who have sought to distinguish the Foremost and Oldman (Pale Beds) formations based on paleoenvironmental criteria (Foremost paralic facies versus Oldman alluvial facies). More significantly, assigning the Herronton sandstone zone to the base of the Oldman Fm continues to allow for meaningful comparisons with the large amount of previously published and archived Oldman Fm subsurface data (e.g., picks, geometries, and isopach-calculations) that are based on the inclusion of the Herronton sandstone zone in that formation [3–6, 8].

The present study also demonstrates that the Oldman Fm thickens substantially to the southeast and southwest, stratigraphically extending the Oldman Fm (and its multiple mudstone and sandstone successions) up-section, well above the “upper siltstone member” as defined and depicted by Hamblin [4]. By Twp 13 (Fig 6) and farther south (Fig 5), the stratigraphically thickened Oldman Fm exhibits multiple stratigraphic intervals dominated by mudstone and sandstones, and no single interval can be deemed equivalent to the “upper siltstone member”. Accordingly, because that term ceases to have any lithostratigraphic meaning south of Twp 16 and is deemed here to be misleading—the “upper siltstone” is equivalent to the lowest portions of the Oldman Fm in southern-most Alberta—it is recommended that the term be dropped altogether in favor of the simpler and more inclusive “Oldman Formation”.

Is there a “Comrey sandstone zone” in southern Alberta?

The Comrey Sandstone was defined and described by Troke [25] in a study of its type exposures in the walls of the Milk River canyon in southeastern Alberta (Twp 02, Rg 06). It occurs stratigraphically in the middle of the Oldman Fm, ~60 m above the base of the formation (Fig 5) and has been successfully utilized by vertebrate paleontologists as a local marker with which to divide the Oldman Fm into lower and upper intervals in the Manyberries region of southeastern Alberta [26, 27]. Similarly, Eberth [7] proposed that the Comrey Sandstone may be part of a laterally continuous and stratigraphically identifiable sandstone zone across southeastern Alberta based on informal observations of outcrop along the South Saskatchewan River, north of Medicine Hat. However, the existence of such a continuous sandstone zone was not tested with subsurface data, and thus remains unproven.

Well logs from cross-section A–A’ (S1 Fig) document the presence of numerous well-developed fluvial sandstones throughout the Oldman Fm, including those with stratigraphic positions that are approximately equivalent to the type Comrey Sandstone (~260–265 m above the MRs). However, because of the ~10 km distance between wells in that cross-section, it is difficult to determine if sandstone bodies are laterally continuous between adjacent wells. Accordingly, a higher resolution study was undertaken in a smaller area along the A–A’ transect (Twp 08–10, Rg 1–4) using data from 22 wells (Fig 1, Table 1, S4 Fig) with one-well-per-section spacing (Figs 12 and 13). This area was chosen because of the large number of available wells that include complete sections of the Oldman Formation, and because the location of the study area is midway between the Milk River canyon area and Oldman Fm exposures along the South Saskatchewan River north of Medicine Hat—the areas that were the basis for the original hypothesis of Eberth [7]. These high-resolution cross-sections (Figs 12 and 13) reveal that although sandstone bodies are abundant throughout the Oldman Fm, there is neither a distinct (e.g., uniquely thick) nor a single, laterally continuous sandstone horizon at the expected stratigraphic position of the type Comrey Sandstone (or anywhere else in the Oldman Fm). The maximum lateral continuity of sandstone bodies in east-west or north-south directions is 16 km (Figs 12 and 13). Based on the absence of unique or laterally continuous sandstones, and the absence of evidence for lateral continuity of the type Comrey Sandstone itself beyond its outcrops south of Manyberries, the hypothesis of a discrete and traceable “Comrey sandstone zone” throughout southern Alberta that is laterally continuous with the type Comrey Sandstone is rejected here.

[Figure omitted. See PDF.]

Based on logs from 10 wells (Fig 1, Table 1, S4 Fig), the cross-section documents the presence of numerous alluvial sandstone bodies (red bars in well logs) throughout the Oldman Fm. Maximum e–w distance of ss-body correlation in this area is ~16 km (red lines), making regional correlation of the Comrey Sandstone or a “Comrey sandstone zone” unlikely.

[Figure omitted. See PDF.]

Based on logs from 12 wells (Fig 1, Table 1, S4 Fig), this cross-section documents numerous alluvial sandstone bodies (red bars in well logs) throughout the Oldman Fm. Similar to the pattern shown in Fig 12, sandstone bodies show a maximum lateral extent of ~10 km, making regional correlation of the Comrey Sandstone or identification of a “Comrey sandstone zone” unlikely.

Clarifying and refining the Judith River-Belly River “problem”

The data presented here can be used to clarify and refine results of a recent study of the Judith River-Belly River clastic wedge by Rogers et al. [9]. That study seeks to resolve lithostratigraphic correlations and nomenclatural differences between the Judith River Formation (USA) and Belly River Group (Canada) and proposes the physical continuity and isochroneity of two previously recognized discontinuities (the mid-Judith discontinuity in Montana [MJD] and the Oldman-Dinosaur Park discontinuity in Alberta [ODPD]). It establishes the single term, Judith River-Belly River discontinuity (JRBRD), to refer to the widespread north-south expression of the combined surfaces. The study further emphasizes the consistent stratigraphic thickness of ~250–265 m between the MRs datum and the JRBRD throughout the field area and, using the U-Pb geochronology of Ramezani et al. [11], assigns an age of ~76.3 Ma to the JRBRD. Lastly, the study builds on the work of Rogers et al. [14] and interprets the combined discontinuities as reflecting a widespread pulse of accommodation resulting from tectonically induced subsidence and the onset of transgression of the Bearpaw Sea at ~76.3 Ma.

Data and interpretations presented here differ somewhat from those of Rogers et al. [9] in presenting a more nuanced picture of the two discontinuities and their geographic relationships. Whereas the Oldman-Dinosaur Park formational contact and discontinuity (ODPD) conforms to the pattern described by Rogers et al. [9] throughout the central and northern portions of our field area (north of Twp 12) the ODPD is neither continuous nor isochronous with the mid-Judith discontinuity (MJD) in the southeastern and southwestern corners of Alberta south of Twp 12 (Figs 9, 11, 14). Specifically, data presented here demonstrate that at Twp 12 and locations farther south to the international border in southeastern Alberta, the ODPD climbs significantly higher stratigraphically and, ultimately, becomes established ~320 m above the Milk River shoulder, ~50 m higher than the MJD (Figs 9 and 14). Similarly, in southwestern-most Alberta at Lethbridge (Twp 10) the ODPD rises sharply higher stratigraphically and becomes established ~357 m above the Milk River shoulder (Fig 11). South of Lethbridge at Twp 08 the Dinosaur Park Fm pinches out altogether and the ODPD ceases to exist (Fig 11; cf. Hamblin [5]).

[Figure omitted. See PDF.]

The model is based on data and correlations presented here from cross-section A–A’ (Fig 1, Table 1, S1 Fig) and differs from that of Rogers et al. [9] in that the ODPD and MJD are litho- and chronostratigraphically separate in southeastern-most Alberta due to up-section thickening of the Oldman Formation and thinning of the Dinosaur Park Fm. More specifically, in southeastern-most Alberta, the ODPD is positioned 320 or more meters above the MRs, whereas the MJD is positioned ~270 m above the MRs. North of Twp 12 in southeastern Alberta, however, the ODPD maintains a consistent position at about 240–260 above the MRs, thereby allowing the recognition of the combined Judith River-Belly River discontinuity (JRBRD) as originally proposed by Rogers et al. [9]. Well log numbers relate to positions shown in Fig 1.

Differences in interpretation between Rogers et al. [9] and the present study are understandable. First, Rogers et al. [9] did not include locations farther west than Rg 12 in their study. Accordingly, they did not integrate the obvious and extreme southwestward thinning and pinch-out of the Dinosaur Park Fm near Lethbridge in their interpretations. Secondly, Rogers et al. [9] failed to integrate into their study the southeastern thickening of the Oldman Fm to the international border as originally described by Eberth and Hamblin [1] and later underscored by isopach (sediment dispersal) data [4–6, 8]. Accordingly, south of Twp 12 in southeastern Alberta and southwestern Saskatchewan, Rogers et al. [9] place many of their well log picks for the ODPD too low in section (their Figs 7, 8). Their misplaced picks likely also resulted from confusion as to how to identify the Oldman-Dinosaur Park formational contact (and thus the ODPD) in gamma-ray logs. As discussed above, Eberth and Hamblin [1] limited their application of the characteristic ODPD gamma-ray log signature to areas in and around Dinosaur Provincial Park (Twp 20, Rg 12) where the ODPD is commonly identified as occurring at the base of a sandstone with low gamma-ray response and sharply in contact with a limited stratigraphic interval of mudstones with high gamma-ray peaks (e.g., Figs 3 and 7). Because Eberth and Hamblin [1] did not emphasize the importance of assessing the gamma-ray responses in the multiple mudstone successions that occur between the Taber and Lethbridge coal zones as the Oldman Formation thickens to the south in Alberta, correct placement of the ODPD has remained problematic. Rogers et al. [9] followed the practice developed by Eberth and Hamblin [1] for identifying the ODPD in subsurface in and around Dinosaur Provincial Park and placed their ODPD picks throughout southern Alberta at the base of the first prominent sandstone that they recognized above the first prominent occurrence of a high gamma-ray response. In so doing, they failed to incorporate the complexity required to accurately identify the ODPD in a significantly thickened Oldman Formation in southeastern Alberta. Lastly, Rogers et al. [9] did not employ measured section outcrop data with which they could calibrate and test their interpretations. Accordingly, whereas they correctly correlated U-Pb ages of the MJD and ODPD in central Montana and at Dinosaur Provincial Park, respectively, and correctly linked two different stratigraphic expressions for the widespread tectonic event(s) that initiated distal foredeep subsidence and influenced the eustatic rise of sea level and marine transgression throughout the northern portion of the WIB, they failed to fully address the changes in BRG stratigraphic architecture just north of the international border.

It is proposed here that the results of Rogers et al. [9] need only be modified slightly by recognizing that the ODPD climbs in section south of Twp 12 (in southeastern Alberta and southwestern Saskatchewan) and that, from Twp 12 to the international border, the two discontinuities (MJD, ODPD) remain stratigraphically separate (Fig 14). In this context, Rogers et al.’s definition of the Judith River Belly River discontinuity (JRBRD) should be revised as representing a widespread isochronous (~76.3 Ma) surface that includes the entire MJD (north-central Montana, southeast Alberta, and southwest Saskatchewan), and that portion of the ODPD that lies north of Twp 12 in southern Alberta and southwestern Saskatchewan. This revised definition necessarily calls attention to a ~60–80 km north-south gap or geographic zone between Twp 14 and Twp 05 in southeastern Alberta (Fig 14) where there exists little evidence for either the MJD or ODPD 250–265 m above the MRs, and thus, no direct evidence for the JRBRD. Here, the geographic gap is interpreted as reflecting a complex interaction between two distinct and regionalized forms of tectonic activity that occurred during the evolution of the combined Judith River-Belly River wedge: (1) the development of a subsidence-induced expansion surface in the Judith River Fm (mid-Judith discontinuity, MJD) of north-central Montana and (2) the tectonically induced rebound of the proximal foredeep in central British Columbia that resulted in the advance of the Dinosaur Park Fm clastic wedge across southern Alberta.

Implications for vertebrate fossil occurrences

The complex internal stratigraphic architecture of the Belly River Group described here has implications for stratigraphic assignment of vertebrate fossils, especially in the Oldman and Dinosaur Park formations. For example, the reciprocal wedged-shaped geometries of the Oldman and Dinosaur Park formations across southern Alberta suggest that some fossil locations that occur in the lower one-half of the Dinosaur Park Fm at northern locations (e.g., Hilda, Sandy Point, Jenner, and Dinosaur Provincial Park) likely are time equivalent with fossil locations that occur in the upper one half of the Oldman Fm at south-eastern locations (e.g., Manyberries area). This possibility has been mentioned or discussed previously [26–28] based on taxonomic comparisons, and the hypothesis can now be tested more rigorously using high-resolution U-Pb CA-ID-TIMS dating of bentonites from each location [9, 11, 12].

Additionally, based on the relatively consistent stratigraphic distance between the Milk River shoulder and the FFm-OFm contact documented here (especially in eastern Alberta) fossil occurrences known to occur stratigraphically close to the FFm-OFm contact [29, 30] may prove to be close in age. This suggests broadly isochronous ages for other more widely dispersed vertebrate fossil occurrences from Milk River Canyon, Redcliff, Suffield, Hays/Rolling Hills, Jenner, Sandy Point, and Rapid Narrows, many of which remain unpublished.

Summary and conclusions

This comprehensive subsurface study of the Belly River Group extends across all southern Alberta and thereby allows previously proposed interpretations of stratigraphic architecture based on smaller geographic areas and more limited stratigraphic intervals to be tested, refined, modified, or rejected. Principal among these is the geometry of the Oldman Fm relative to both the underlying Foremost Fm and the overlying Dinosaur Park Fm. Subsurface data presented here confirms thickening of the Oldman Fm to the south and principally to the southwest into Montana, and documents stratigraphic thickening of the Oldman Fm primarily up-section at the expense of the overlying Dinosaur Park Fm. It also underscores the increasingly younger, diachronous nature of the Oldman-Dinosaur Park formational contact to the south, and the widespread possibly isochronous nature of the Foremost-Oldman formational contact across southern Alberta. Lastly, as the Dinosaur Park Fm thins to the south and pinches out to the southwest on the Oldman Fm, the depositional edge of the formation is resolvable as a wnw–ese trending line, and low gamma/high porosity signatures characteristic of the Lethbridge coal zone continue south into the uppermost Oldman Fm.

Up-section stratigraphic thickening of the Oldman Fm requires that our understanding of the formation in the subsurface be modified and broadened in several ways. First, although the higher gamma-ray response of mudstone successions through the formation remains a reliable means for identifying the formation in well logs and distinguishing it from the Foremost and Dinosaur Park formations it is also clear that there exist many sandstones with low gamma-ray responses within the Oldman Fm that are equivalent in response to those in the underlying and overlying formations. Thus, picking the top or the bottom of the formation requires careful consideration of the entire gamma-ray response pattern between the Taber and Lethbridge coal zones of the Belly River Group, an approach that has yet been applied consistently by students of BRG stratigraphy. Up-section thickening of the Oldman Fm also forces a reconsideration of two terms (“Comrey Member sandstone” and “upper siltstone member)” that were previously applied in reference to a two-fold subdivision of the Oldman Fm north of Twp 19. Recognition that the type Comrey Sandstone occurs in the middle of the Oldman Formation in southeastern Alberta suggests that previous use of that term for the sandstone zone that forms the base of the Oldman Fm throughout large portions of southern Alberta is incorrect and should be replaced by the term Herronton sandstone zone. The Herronton sandstone zone is also recognized at the base of the Oldman Fm above the Taber coal zone in local outcrops along the Red Deer, Bow, South Saskatchewan, and Milk rivers of southern Alberta. The “upper siltstone member” was originally proposed for a stratigraphically limited interval (typically less than 50 m thick). However, in both southeastern and southwestern Alberta, where the formation thickens to more than 100 m or more, the term is neither stratigraphically accurate nor useful, and thus should be avoided.

Rogers et al.’s [9] correlation of discontinuities from both the Judith River Formation and the Belly River Group across the international border also requires reconsideration given the improved understanding of Belly River Group stratigraphic architecture provided here. It is clear now that different and regionalized tectonic influences on stratigraphic architecture were at work in the Judith River Formation and the Belly River Group during the onset of the Bearpaw transgression. In particular, the development of a subsidence-induced expansion surface in the Judith River Fm (mid-Judith discontinuity, MJD) versus a tectonically induced rebound of the proximal foredeep that resulted in a reorganization of sediment distribution systems across southern Alberta (the Oldman-Dinosaur Park discontinuity, ODPD). New U-Pb geochronologic data indicate that the initiation of these regionalized events was closely related in time, likely tectonically related, and thus correlable [11, 12, 14]. However, data presented here also show that the ODPD is locally diachronous and becomes significantly younger south of Twp 12 as the Dinosaur Park Fm overrides the Oldman Fm in a step-wise fashion. This pattern indicates that Judith River Belly River discontinuity—the term referring to the combined MJD and ODPD—consists of the MJD but only the oldest portions of the ODPD (north of Twp 12).

These data and modifications to our understanding of Belly River stratigraphic architecture will likely improve correlation of fossiliferous exposures in southern Alberta as well as correlations of the Belly River Group with the Judith River and Two Medicine formations in northern Montana.

Supporting information

S1 Fig. Annotated well logs (n = 35) used in cross-section A–A’.

https://doi.org/10.1371/journal.pone.0292318.s001

S2 Fig. Annotated well logs (n = 28) used in cross-section B–B’.

https://doi.org/10.1371/journal.pone.0292318.s002

S3 Fig. Annotated well logs (n = 24) used in cross-section C–C’.

https://doi.org/10.1371/journal.pone.0292318.s003

S4 Fig. Well logs (n = 22) used in high-resolution cross-sections in Figs 12 and 13.

https://doi.org/10.1371/journal.pone.0292318.s004

Acknowledgments

Divestco Geoscience Ltd. (Calgary) and staff are thanked for providing access to their raster and digital well log library, and for assistance in use of their digital platform (Energisite). Specifically, I thank Tania Uszacki (Product Manager) and Veronica Holmes for support that made this project possible. Ray Rogers was instrumental in revitalizing my interest in this project and incentivising the completion of research initiated in 1986. I thank him for the many interesting and stimulating discussions concerning discontinuities and upper Belly River stratigraphy discussed here. I thank Federico Fanti and David Evans for graciously sharing their Belly River outcrop data and opinions from the Milk River canyon area of southeastern Alberta. Field assistance was provided by Tyrrell Museum technical staff over many decades and additional support was supplied by the Tyrrell Museum Cooperating Society. Special thanks to Marty who shared in the day-to-day advances and setbacks inherent in developing a data-heavy manuscript.

Citation: Eberth DA (2024) Stratigraphic architecture of the Belly River Group (Campanian, Cretaceous) in the plains of southern Alberta: Revisions and updates to an existing model and implications for correlating dinosaur-rich strata. PLoS ONE 19(1): e0292318. https://doi.org/10.1371/journal.pone.0292318

About the Authors:

David A. Eberth

Roles: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

E-mail: [email protected]

Affiliation: Royal Tyrrell Museum of Palaeontology, Drumheller, Alberta, Canada

ORICD: https://orcid.org/0009-0009-4671-7078

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

References

1. Eberth DA, Hamblin AP. Tectonic, stratigraphic, and sedimentologic significance of a regional discontinuity in the upper Judith River Group (Belly River Wedge) of southern Alberta, Saskatchewan, and northern Montana. Canadian Journal of Earth Sciences. 1993; 30: 174–200.

2. Jerzykiewicz T, Norris DK. Stratigraphy, structure and syntectonic sedimentation of the Campanian ’Belly River’ clastic wedge in the southern Canadian Cordillera. Cretaceous Research. 1994; 15: 367–399.

3. Hamblin AP, Abrahamson BW. Stratigraphic architecture of "Basal Belly River" cycles, Foremost Formation Belly River Group, subsurface of southern Alberta and southwestern Saskatchewan. Bulletin of Canadian Petroleum Geology. 1996; 44: 654–673.

4. Hamblin AP. Stratigraphic architecture of the Oldman Formation, Belly River Group, surface and subsurface of southern Alberta. Bulletin of Canadian Petroleum Geology. 1997; 45: 155–177.

5. Hamblin AP. Regional distribution and dispersal of the Dinosaur Park Formation, Belly River Group, surface and subsurface of southern Alberta. Bulletin of Canadian Petroleum Geology. 1997; 45: 377–399.

6. Hamblin AP, Lee PJ. Uppermost Cretaceous, Post-Colorado Group gas resources of the Western Canada Sedimentary Basin, Interior Plains. Geological Survey of Canada Bulletin. 1997; 518, 88p.

7. Eberth DA. The Geology. In Dinosaur Provincial Park: A Spectacular Ancient Ecosystem Revealed. Edited by PJ Currie and EB Koppelhus. Indiana University Press, Bloomington, 2005; pp. 54–82.

8. Glombick P. Subsurface stratigraphic picks for the top of the Oldman Formation (base of Dinosaur Park Formation), Alberta Plains: Energy Resources Conservation Report/Alberta Geological Survey, Open File Report. 2011; 13, 34 p.

9. Rogers RR, Eberth DA, Ramezani J. The “Judith River-Belly River Problem” revisited: New perspectives on the correlation of Campanian dinosaur-bearing strata in Montana and Alberta based on a revised stratigraphic model updated with CA-ID-TIMS U-Pb geochronology. Geological Society of America Bulletin. 2023.

10. Brown CM, Ryan MJ, Evans DC. A census of Canadian dinosaurs: more than a century of discovery. In All Animals Are Interesting: A Festschrift in Honour of Anthony P. Russell. Edited by ORP Bininda-Edmonds, Powell GL, Jamniczky HA, Bauer AM, Theodor J. BIS-Verlag, Oldenburg, Germany. 2015; pp. 151–209.

11. Ramezani J, Beveridge TL, Rogers RR, Eberth DA, Roberts EM. Calibrating the zenith of dinosaur diversity: high-resolution Campanian chronostratigraphy of the Western Interior Basin based on CA-ID-TIMS U-Pb zircon geochronology: Scientific Reports 2022; 12:16026.

12. Eberth DA, Evans DC, Ramezani J, Kamo SL, Brown CM, Currie PJ, et al. Calibrating geologic strata, dinosaurs, and other fossils at Dinosaur Provincial Park (Alberta, Canada) using a new CA-ID-TIMS U–Pb geochronology. Canadian Journal of Earth Sciences. 2023;

13. Gilbert MM. Sedimentology, ichnology, sequence stratigraphy, and vertebrate paleontology of the Belly River Group, southwestern Saskatchewan, Canada [Ph.D. thesis]: Saskatoon, University of Saskatchewan. 2019; 398 p.

14. Rogers RR, Kidwell SM, Deino A, Mitchell JP, Nelson K. Age, correlation, and lithostratigraphic revision of the Upper Cretaceous (Campanian) Judith River Formation in its type area (north-central Montana), with a comparison of low- and high-accommodation alluvial records. Journal of Geology. 2016; 124: 99–135,

15. Kauffman EG, Caldwell WGE. The Western Interior Basin in space and time. In Evolution of the Western Interior Basin. Edited by EG Kauffman and WGE Caldwell. Geological Association of Canada, Special Paper. 1993; 39: 1–30.

16. Williams GD, and Burk CF. Upper Cretaceous. In: Geological History of Western Canada. Edited by RG McCrossan and RP Glaister. Alberta Society of Petroleum Geologists. 1964; p. 169–189.

17. Russell LS, Landes RW. Geology of the southern Alberta plains. Geological Survey of Canada, Memoir. 1940; 221, 223 p.

18. McLean JR. Stratigraphy of the Upper Cretaceous Judith River Formation in the Canadian Great Plains. Saskatchewan Research Council Geology Division Report. 1971; 11, 97 p.

19. Glombick P, Mumpy AJ. Subsurface stratigraphic picks for the top of the Milk River ’shoulder’, Alberta Plains: Including tops for the Milk River Formation and Alderson Member of the Lea Park Formation. AER/AGS Open File Report. 2013; 17.

20. Eberth DA, Braman DR, Tokaryk TT. Stratigraphy, sedimentology and vertebrate paleontology of the Judith River Formation (Campanian) near Muddy Lake, west-central Saskatchewan: Bulletin of Canadian Petroleum Geology. 1990; 38: 387–406.

21. Macdonald DE, Ross TC, McCabe PJ, Bosman A. An evaluation of the coal resources of the Belly River Group to a depth of 400 m in the Alberta plains: Alberta Research Council, Open File Report. 1987; 8, 118 p.

22. Dawson FM, Evans CG, Marsh R, Richardson R. Chapter 24: Uppermost Cretaceous and Tertiary strata of the Western Canada Sedimentary Basin. In Geological Atlas of the Western Canada Sedimentary Basin, Compiled by GD Mossop, I Shetsen. Canadian Society of Petroleum Geologists and Alberta Research Council, 1994. https://ags.aer.ca/atlas-the-western-canada-sedimentary-basin/chapter-24-uppermost-cretaceous-and-tertiary-strata. Accessed on line April 3, 2023.

23. Dowling DB. The southern plains of Alberta: Geological Survey of Canada, Memoir 1917; 93, 200 p.

24. Gordon J. Stratigraphy< and sedimentology of the Foremost Formation in southeastern Alberta and southwestern Saskatchewan. M.Sc. thesis. University of Regina, Regina. 2000; 110 p.

25. Troke CG. Sedimentology, stratigraphy and calcretes of the Comrey Member, Oldman Formation (Campanian), Milk River Canyon, southeastern Alberta: MSc. Thesis, University of Calgary, Calgary. 1993; 146 p.

26. Peng J-h. Palaeoecology of vertebrate assemblages from the Upper Cretaceous Judith River Group (Campanian) of southeastern Alberta, Canada. Ph.D. Thesis, University of Calgary, Calgary. 1998; 316 p.

27. Chiba K, Ryan MJ, Eberth DA, Braman D, Kobayashi Y, Evans DC. Taphonomy of a Centrosaurus apertus (Dinosauria: Ceratopsidae) bonebed from the Oldman Formation of Alberta, Canada. PALAIOS. 2015; 30(9): 655–667.

28. Ryan MJ, Russell AP. A new centrosaurine ceratopsid from the Oldman Formation of Alberta and its implications for centrosaurine taxonomy and systematics. Canadian Journal of Earth Sciences. 2005; 42: 1369–1387.

29. McKellar RC, Wolfe AP. Canadian amber. In Biodiversity of Fossils in Amber from the Major World Deposits. Edited by D Penney. Siri Scientific Press, Rochdale, UK. 2010; pp. 96–113.

30. Cullen T, Fanti F, Capobianco C, Ryan MJ, Evans DC. A vertebrate microsite from a marine-terrestrial transition in the Foremost Formation (Campanian) of Alberta, Canada, and the use of faunal assemblage data as a paleoenvironmental indicator. Palaeogeography, Palaeoclimatology, Palaeoecology. 2016; 444: 101–114.