An ankylosaurian dinosaur from the Cenomanian Dunvegan Formation of northeastern British Columbia, Canada

Fragmentary but associated dinosaur bones collected in 1930 from the Pine River of northeastern British Columbia are identified here as originating from an ankylosaur. The specimen represents only the second occurrence of dinosaur skeletal material from the Cenomanian Dunvegan Formation and the first from Dunvegan outcrops in the province of British Columbia. Nodosaurid ankylosaur footprints are common ichnofossils in the formation, but the skeletal material described here is too fragmentary to confidently assign to either a nodosaurid or ankylosaurid ankylosaur. The Cenomanian is a time of major terrestrial faunal transitions in North America, but many localities of this age are located in the southern United States; the discovery of skeletal fossils from the Pine River demonstrates the potential for the Dunvegan Formation to produce terrestrial vertebrate fossils that may provide important new data on this significant transitional period during the Cretaceous.


Introduction
Terrestrial vertebrate fossils have rarely been reported from the Cenomanian-aged Dunvegan Formation (Burns and Vavrek et al., 2014), which crops out in northeastern British Columbia and northwestern Alberta. Here we describe dinosaur bones recovered from outcrops of the formation along the Pine River in northeastern British Columbia (Fig. 1a). The specimen (CMN 59667), a block of sandstone with associated ribs and vertebrae, was collected in 1930 from the Pine River. To our knowledge this is the first dinosaur skeletal fossil discovered in the province, predating by more than 40 years the discovery of Ferrisaurus sustutensis (Arbour and Evans, 2019), first reported by Arbour and Graves (2008).
Although bony fish and sharks have been recorded from the Dunvegan Formation in Alberta (Hay et al., 2007;Cook et al., 2008;Vavrek et al., 2014), the terrestrial vertebrate body fossil record of the formation has thus far been limited to a few ankylosaur ossicles from Alberta (Burns and . Dinosaur bones found in the Tumbler Ridge area, first described as originating from the Dunvegan Formation (Mc-Crea, 2003), were later reinterpreted as originating from the Turonian Kaskapau Formation (McCrea and Buckley, 2004;Rylaarsdam et al., 2006). Therefore, well-preserved dinosaur fossils, particularly from an associated skeleton, would contribute significantly to understanding the fauna of the Dunvegan Formation and preservational potential within this unit.
CMN 59667 is noteworthy for being one of very few Cenomanian-aged dinosaur skeletal fossils from Canada (Brown et al., 2015). Terrestrial Cenomanian assemblages are rare in North America but capture an important time which preserves a distinct and transitional assemblage between characteristic faunas of the Early and latest Late Cretaceous. The nearest comparative material of this age is derived from Dunvegan Formation outcrops in Alberta, the Cedar Mountain Formation of Utah (Avrahami et al., 2018;Tucker Figure 1. CMN 59667 was collected along the Pine River in northeastern British Columbia, approximately 1.6 km downstream of its junction with the Murray River. Sedimentological and palaeontological features in this area suggest these are outcrops of the upper two-thirds of allomember E of the Dunvegan Formation. (a) Location of the collection site for CMN 59667 in British Columbia, Canada, indicated by a red star (inset map shows location in western Canada). (b) Tetrapodosaurus pes natural cast (RBCM P1473) discovered at this site; anterior is to the right. (c) RBCM P1028, sample with Metasequoia leaf and obovate or elliptical angiosperm leaf with entire margin (leaf morphotype 122; Ash et al., 1999). (d) RBCM P1030, palmately lobed angiosperm leaf (leaf morphotypes 140-143; Ash et al., 1999). (e) Cliff-forming outcrops on the eastern bank of the Pine River, 1.6 km from its junction with the Murray River, as viewed from the west bank. (f) Gymnosperm cone in cross section (not collected), observed on the west bank of the Pine River. et al., 2020), the Woodbine (Noto et al., 2012) and Paw Paw formations of Texas, the Frontier Formation of Wyoming (Lull, 1921), the Wayan Formation of Idaho (Krumenacker et al., 2019), and the Dakota Formation of the western US (Eaton, 1960;Eaton, 1993), all of which are either areas that are being actively collected and researched or not well explored or documented. The Dunvegan Formation is the furthest north of the known Cenomanian-aged vertebrate fossilbearing terrestrial units in North America and is therefore of particular interest for testing proposed latitudinal faunal provinciality in the early Late Cretaceous (e.g. Lehman, 2001;Sampson et al., 2010) and understanding its origin.

Materials and methods
CMN 59667 was loaned to the Royal Ontario Museum for preparation and study. Mechanical preparation using air scribes exposed the bones from the surrounding matrix. Parts of several bones were preserved as natural impressions of their external surface in the sandstone block, in which threedimensional fossil bone itself had been weathered away. The easy-to-interpret impressions of the rib shafts were left in place, but a portion of what turned out to be a neural arch was more complex and difficult to interpret as an impression. In order to recover this anatomical data, the impression was used as a natural mold and infilled with grey epoxy putty. Preparation then proceeded to the epoxy cast of the natural mold surface, providing a complete surface view of the right lateral side of the neural arch.
Several of us (Victoria M. Arbour, Derek Larson, and David Evans) undertook a brief visit in May 2019 to outcrops near the junction of the Murray and Pine rivers (the site of origin for the specimen) to attempt to relocate the original site where the fossil had been collected and collect geological and stratigraphic data at the locality (Fig. 1a). River levels were high, preventing extensive exploration of the steep cliff-forming outcrops. 3 Geological setting CMN 59667 was collected by Merton Yarwood Williams in 1930; this was likely during a survey he headed, the Non-Metallic Mineral Investigation on the Pacific Great Eastern Railway Survey of Resources of the Peace River, West Cariboo and West Lillooet Blocks, B.C., during the summers of 1929and 1930(Hives, 1988. Notes associated with CMN 59667 indicate it was collected "1 mile below junction of Murray and Pine Rivers (Peace River district)". Charles M. Sternberg later that year identified the block as containing an anterior dorsal vertebra and rib fragment of an ornithischian resembling Camptosaurus, according to the label associated with the specimen.
The early to middle Cenomanian Dunvegan Formation crops out in this area. The Dunvegan Formation represents a major deltaic complex deposited over the span of about 2 million years (Plint, 2000). Characterized by a succession of alluvial and shallow marine sandstones, siltstones, and shales, it has an interfingering and diachronous relationship with the underlying marine Shaftesbury Formation and overlying marine Kaskapau Formation (Bhattacharya, 1994;Plint, 2000).
We undertook a reconnaissance field trip to the Pine River-Murray River area in May 2019 to document the lithology and fossils of the original collection area. We were able to access the west bank of the Pine River on foot 1.3 km downstream from its junction with the Murray River, very close to the locality described in the original notes associated with the specimen. The east bank of the river exposes a large section of Dunvegan Formation approximately 22 m thick ( Fig. 1e). At the base of the cliff we observed a distinctive orangestained palaeosol overlain by two thin dark-grey coal seams or palaeosols. Above these layers, there were seven packages of fine-grained sediments overlain by a more resistant layer of sandstone. The west side of the river lacks cliff-forming outcrop, but small outcrops are present as well as abundant loose pieces of sandstone.
Although we did not encounter any additional vertebrate bone, we observed current ripple marks, siderite nodules, and abundant ichnofossils, including vertically oriented burrows, shallow surface burrows or feeding traces, and the dinosaur ichnotaxon Tetrapodosaurus (RBCM P1473), made by a nodosaurid ankylosaur (Fig. 1b). Previously, Storer (1975) described a slab with tridactyl dinosaur footprints discovered on the north bank of the Pine River "at the site of the present bridge at East Pine", very close to the reported collection locality of CMN 59667. These were identified as Colombosauripus ungulatus and interpreted as belonging to a small coelurosaur similar to an ornithomimid. Terrestrial vertebrate footprints and trackways are widespread throughout the Dunvegan formation (both geographically and strati-graphically) and include tracks made by turtles, crocodilians, nodosaurs (Tetrapodosaurus), large ornithopods, small theropods (cf. Irenichnites; Columbosauripus), mediumsized theropods (cf. Magnoavipes), and birds; tracks attributed to large theropods have yet to be reported (McCrea et al., 2014). Plant remains in this area included leaf compressions (including the representative samples we collected, RBCM P1028-P1032), carbonized leaves and cones, and wood impressions (Fig. 1c, d, f). These sedimentological and palaeontological features are consistent with the upper twothirds of allomember E of the Dunvegan Formation, dominated by lakes, crevasse deltas, and poorly drained palaeosols (Lumsdon-West and Plint, 2005). This region would have been located at about 60 • N during the Cenomanian (van Hinsbergen et al., 2015).

Material
CMN 59667 consists of the posterior half of a dorsal vertebral centrum, an almost complete dorsal neural arch, and small sections of two dorsal ribs. All of these bones occur closely associated within a single small block of sandstone, and their size and morphology are consistent with being from a single individual, despite the lack of articulation.

Locality and horizon
CMN 59667 was collected from Dunvegan Formation exposures along the Pine River, ∼ 1.6 km north of the junction of the Murray and Pine rivers; precise coordinates are on file at the CMN. CMN 59667 is hosted in a dark-grey, mediumgrained litharenite.

Description and comparisons
Preserved in association in the block are a partial dorsal neural arch and a centrum in close proximity but disarticulated, so it is unclear whether these represent a single vertebra or are parts of two vertebrae (Fig. 2, Table 1). The centrum is taller than wide as preserved on the centrum face. It is constricted at the mid-length, giving it an hourglass-like appearance (Fig. 3b). Although one of the articular faces is damaged, the anterior and posterior ends of the centrum are dorsoventrally offset from each other. The neural arch includes the right transverse process, neural spine, and right postzygapophysis. The neural spine is dorsoventrally low, mediolaterally narrow, and anteroposteriorly long, with an overall square outline ( Fig. 2a-b). The dorsal edge of the neural spine is lightly rugose but not strongly mediolaterally expanded. The transverse process is broken and slightly offset from the neural arch but would have been slightly dorsally directed. The diapophysis is located ventrally on the transverse process towards its distal terminus and is an elongate, triangular groove (Fig. 3c).
Two partial dorsal rib shafts closely associated with the vertebral remains are preserved in the block (Fig. 2c-d). The largest section of rib is preserved in part as an impression in the sandstone block. The preserved portion of the fossil (133.4 mm) is 78 % of the total length of the section, including the impression (17.1 mm). The shaft is anteroposteriorly broad and robust, with an oval cross section that measures 32.2 mm wide and 8.7 mm tall (measured at its distal-most end). The flat lateral and medial surfaces indicate that the preserved section originated from the mid to distal end of the rib. The second rib is represented by an 80.5 mm thin splint of preserved bone from its anterior or posterior margin. The long axis of the second rib measures 160.3 mm long (including its impression), is oriented parallel to the first rib, and is on the same plane, suggestive of articulation.
CMN 59667 is identified here as an ankylosaur based on several aspects of its vertebral morphology. In iguanodontians and ceratopsians, the neural spine of the dorsal vertebrae is relatively dorsoventrally tall and elongate and anteroposteriorly short. The preserved neural spine of CMN 59667, in contrast, is dorsoventrally short but anteroposteriorly long, with a square outline in lateral view (Fig. 3a). This morphology is typical of the dorsals of ankylosaurs (Fig. 4) Fig. 4b), Polacanthus (Blows, 1987), Sauropelta (Ostrom, 1970), and Europelta (Kirkland et al., 2013). Most importantly, the morphology of the transverse process in CMN 59667 is most similar to ankylosaurs compared with other ornithischians. In most ornithischians, the cross section of the distal transverse process is triangular or t-shaped and lacks a diapophyseal sulcus. The grooved morphology in CMN 59667 is consistent with the typical morphology of the transverse process found in ankylosaurids (Fig. 4), such as Akainacephalus (Wiersma and Irmis, 2018), Ankylosaurus (AMNH 5895; Carpenter, 2004), Euoplocephalus (AMNH 5337; Arbour and Currie, 2013), and cf. Talarurus (MPC KID 155), which accepts the rib head; this groove seems to be absent in many nodosaurid ankylosaurs, such as Europelta (Kirkland et al., 2013), Hungarosaurus (Ősi, 2005, Peloroplites (Carpenter et al., 2008), Polacanthus (NHMUK R9293), and Silvisaurus (Carpenter and Kirkland, 1998). We are unaware of this morphology in any other group of ornithischian dinosaurs as well. The incomplete ribs of CMN 59667 are consistent with the morphology seen in other large-bodied ornithischian dinosaurs that are known from middle Cretaceous strata of North America, including ankylosaurs and iguanodontians, but cannot be identified more precisely.

Discussion
Although the Early to mid-Cretaceous dinosaur fossil record of western Canada is currently less well documented (Brown et al., 2015) in comparison to the southwestern United States, recent discoveries reveal the potential for much more information to be gleaned from this time period. Notably, the Early and mid-Cretaceous dinosaur skeletal record of western Canada is almost entirely represented by ankylosaurs, including osteoderms from the Berriasian-aged Pocaterra Creek Member of the Cadomin Formation in the Rocky Mountains of southwest Alberta (Nagesan et al., 2020), the exceptionally preserved Borealopelta from the Aptian Clearwater Formation (Brown et al., 2017), and isolated ossicles from the Dunvegan Formation of Alberta (Burns and .
Ankylosaur skeletal fossils are widespread across North America during the Campanian-Maastrichtian, but they are never common components of the fossil record for a given formation. Even in the well-sampled Dinosaur Park Formation of southern Alberta, ankylosaurs make up less than 15 % of the skeletal record (Brinkman et al., 1998(Brinkman et al., , 2005. Although similar statistical analyses or censuses have yet to be undertaken for many Early and mid-Cretaceous formations in North America, ankylosaurs appear to make up a greater proportion of the ornithischian fauna compared to the Campanian-Maastrichtian in several localities (e.g. Britt et al., 2009). Most Early and mid-Cretaceous Laramidian ankylosaurs were nodosaurids (Carpenter and Kirkland, 1998;Arbour et al., 2016;Brown et al., 2017), with only a single ankylosaurid, Cedarpelta bilbeyhallorum, known from the Cenomanian-aged Mussentuchit Member of the Cedar Mountain Formation (Carpenter et al., , 2008. Ankylosaur footprints are also rare during the Campanian-Maastrichtian of North America, with only a few potential records from the Blackhawk Formation of Utah (McCrea et al., 2001), the Wapiti Formation of northern Alberta (Fanti et al., 2013), and the Chignik Formation of Alaska (Fiorillo et al., 2019). In contrast, the ichnotaxon Tetrapodosaurus, generally considered to have been made by tetradactyl ankylosaurs (McCrea et al., 2001(McCrea et al., , 2014, is widespread in North America during the Early and mid-Cretaceous, with records from the Nanushuk and Chandler formations of Alaska (Mc-Crea et al., 2001;May and Druckenmiller, 2009); the Gates and Dunvegan formations of Alberta (McCrea et al., 2001); the Gorman Creek, Gething, Dunvegan, Gates, and Kaskapau formations of British Columbia (McCrea et al., 2001(McCrea et al., , 2014; the Dakota Group of Colorado (Lockley and Gierlinski, 2014); the Naturita and Cedar Mountain formations of Utah (Lockley et al., 1999(Lockley et al., , 2018; and the Ross River Block in the Yukon (Gangloff et al., 2004). Given the rarity of ankylosaurian footprints in the Campanian-Maastrichtian of North America but their abundance in the Early and mid-Cretaceous (including the Dunvegan Formation), it is notable that the only dinosaur skeletal material thus far described from the Dunvegan Formation, both recorded here and by Burns and Vavrek (2014), is ankylosaurian in nature.
Based on the prevalence of Tetrapodosaurus in the Dunvegan Formation (McCrea et al., 2014), including at the CMN 59667 locality (e.g. RBCM P1473), and the fact that its trackmaker was probably not an ankylosaurid ankylosaur, we might expect ankylosaur skeletal material from the Dunvegan Formation to derive from a nodosaurid ankylosaur. However, the preserved transverse process of CMN 59667 has a diapophyseal sulcus on its ventral side, a feature that seems to be present in ankylosaurids but absent in nodosaurids. The distribution of the diapophyseal sulcus within ankylosaurs requires additional specimens and investigation, and CMN 59667 is too fragmentary to confidently assess its affinities within Ankylosauria. Nevertheless, the presence of both nodosaurids and ankylosaurids further south in North America during the Cenomanian means that an ankylosaurid identification for this specimen cannot be ruled out, despite the prevalence of nodosaurid footprints in the Dunvegan Formation. Burns and  were similarly cautious in their taxonomic appraisal of isolated ankylosaur ossicles from Dunvegan outcrops in Alberta, noting that the morphology and histology of such small elements cannot be as easily mapped to higher-level clades in comparison to larger ankylosaur osteoderms. We note, however, that the ossicles reported by Burns and  are extremely small (∼ 1 mm diameter), a size that is much more consistent with small interstitial osteoderms from ankylosaurid ankylosaurs such as Jinyunpelta (Zheng et al., 2018), Scolosaurus (NHMUK R5161; Nopcsa, 1928), and Zuul (ROM 75860; Arbour and Evans, 2017). Smaller interstitial osteoderms in nodosaurids are sometimes absent because of closer osteoderm packing (e.g. Borealopelta; Brown, 2017;Brown et al., 2017) or are >8 mm in diameter, such as Sauropelta (AMNH 3032) and Panoplosaurus (CMN 2759).
The Cenomanian is a key stratigraphic interval for interpreting the palaeo-biogeographic and evolutionary history of ankylosaurid ankylosaurs. The earliest records of the characteristic tail club occur in Albian to Cenomanian-aged formations in Asia (Arbour and Currie, 2015;Zheng et al., 2018). Ankylosaurids are also first documented in North America during this time period in the Mussentuchit Member of the Cedar Mountain Formation but then disappear from the North American fossil record until the Santonian (Baszio, 1997;Eaton et al., 1999;Parrish, 1999), possibly indicating a response to high sea levels in Laramidia during the early Late Cretaceous (Arbour et al., 2016), although poor sampling in the Turonian-Santonian cannot be ruled out (e.g. Upchurch et al., 2011). In contrast, nodosaurids have a nearly continuous fossil record throughout the entire Cretaceous of North America (Arbour et al., 2016).
CMN 59667 represents the first associated skeletal material from a dinosaur recovered from the Dunvegan Formation. This new, well-preserved specimen therefore shows that this important taphonomic mode (articulated and associated skeletons, modes A1 and A2 of Eberth and Currie, 2005) can occur in the unit and raises the possibility that continued prospecting efforts could yield substantially more complete dinosaur material from this poorly known time interval in the Cretaceous fossil record of Canada. The presence or absence of ankylosaurids in the Dunvegan Formation can only be determined through the discovery of additional, more di-agnostic material, but the morphology of the currently known ankylosaur skeletal remains is tantalizing potential evidence for Cenomanian-aged ankylosaurids in Canada. Similarly, the discovery of unequivocal, diagnostic nodosaurid remains would provide the opportunity to link a trackmaker to the abundant Tetrapodosaurus trackways known from this formation.
Data availability. All specimens discussed within this paper are permanently reposited within museum collections. Precise locality data for CMN 59667 are on file at the Canadian Museum of Nature.
Author contributions. VMA, DL, and DE collected and analyzed the data. VMA wrote the manuscript, and all authors edited and revised the manuscript prior to submission.
Competing interests. The authors declare that they have no conflict of interest.