High-Latitude Depositional Systems, Provenance, and Basinal Setting of the Late Cretaceous Cantwell Basin, Denali National Park and Preserve, Alaska: A Stratigraphic Framework for Paleontological and Paleoclimatic Studies

: The Cantwell Formation of the central Alaska Range provides a robust archive of high-latitude, Late Cretaceous depositional systems and paleo-ﬂoral/faunal assemblages. Our stratigraphic analysis deﬁnes two mappable members. The lower member (1500–2000 m thick) represents vegetated alluvial fan and braided ﬂuvial systems that transition up-section to ﬂuvial–estuarine systems that drained into an inland continental seaway. The upper member (~2000 m thick) represents estuarine–marginal marine and lacustrine systems. Previous paleontological studies demonstrate that the Cantwell basin was populated by various dinosaurs, ﬁshes, bivalves, birds


Introduction
The Cantwell basin, located in the central Alaska Range, contains a rich archive of Late Cretaceous high-latitude depositional systems, paleobiological ecosystems, and tectonics of the northwestern Cordillera of North America. The basin consists of 4 km of siliciclastic strata, the Cantwell Formation, exposed between the Hines Creek and Denali fault systems (Figures 1 and 2) [1,2]. The Cantwell Formation is characterized by thick successions of conglomerate, sandstone, mudstone, and minor coal and limestone interpreted to represent alluvial fan, braided stream, lacustrine, and marginal marine environments [2,3]. Since the initial geologic mapping and sedimentary studies of the Cantwell Formation, abundant trace fossils of theropods, ceratopsians, hadrosaurs, therizinosaurs, fishes, and pterosaurs have been documented, along with high-latitude paleoclimate records e.g., Refs. [4][5][6][7][8][9]. The Cantwell Formation also contains a record of the syndepositional and post-depositional deformational processes that constructed the alpine topography of the central Alaska Range [1,2,10,11]. Despite the importance of the Cantwell basin from a stratigraphic, paleobiologic, and regional tectonic perspective, an integrated stratigraphic and structural framework has not been fully developed.
The primary goal of our study is to establish a detailed stratigraphic framework for the Cantwell basin, integrating new and previously published geologic mapping data to reconstruct depositional systems at a basin scale from deformed strata of the Cantwell Formation. Additionally, we use modern techniques for provenance analysis to develop a more regional perspective on the paleogeographic and tectonic context of the Cantwell basin. With these goals in mind, we present the following components of our analysis: (1) a composite representative section that we correlate throughout the basin; (2) an updated geologic map compilation and three structural cross-sections within the most accessible and continuous exposures of the Cantwell Formation in Denali National Park and Preserve (i.e., the western part of the basin); (3) U-Pb detrital zircon ages (N = 13; n = 4657) with stratigraphic control to add a new perspective to a previous composition-based provenance analysis [2]; and (4) an integrated interpretation of the basin configuration and regional paleogeographic setting of the Cantwell basin.
The results of our study provide an updated basin framework relevant to future studies of Late Cretaceous paleontology, paleoclimate, and tectonics in the northern North American Cordillera. Reconstruction of paleo-environments associated with vertebrate fossil-bearing assemblages has the potential to reveal how high-latitude Cretaceous depositional systems and their related ecosystems responded to changing climate conditions e.g., Ref. [7]. Changes in climate are known to have particularly profound impacts in Arctic environments and these changes are often best preserved in the sedimentary record e.g., Refs. [12,13]. Our analysis, hopefully, will also provide additional geologic content for the over 500,000 annual visitors to Denali National Park and Preserve that view and interact with this scenic landscape.

Paleontology and Paleoclimate
Ongoing interest in the Cantwell Formation has stemmed from recent discoveries of abundant bird and dinosaur trackways, including theropods, ceratopsians, pterosaurs, therizinosaurs, and hadrosaurs e.g., Refs. [4,6,[14][15][16][17][18]. Various bird trackways, fishes, and crayfish burrows have also been documented [4,5]. The most prolific dinosaur and bird tracksites are associated with fluvial and overbank deposits, primarily in the middle-upper part of the Cantwell Formation e.g., Ref. [6]. These strata are also rich in leaf fossils, wood fragments, and other organic matter. Stable isotopes from these materials have been used to determine annual temperature ranges and serve as a proxy for paleo-precipitation rates. These data yield a mean annual temperature of 7.4 ± 2.4 • C, with an annual temperature range from −2.3 ± 3.8 • C to 17.1 ± 3.2 • C (mean temperatures for the coldest and warmest months, respectively) [9].
Along with the Cantwell basin, vertebrate fossil assemblages, including several dinosaur species, have been reported in Mid-Upper Cretaceous strata throughout Alaska. In southern Alaska, dinosaur trackways have been identified in forearc basin strata of the Matanuska basin [19] and the Chignik Formation of the Alaska Peninsula [20]. Trackways of theropods and ornithopods have been identified in proximal foreland basin strata of the Beaver Lake formation in the Nutzotin Mountains of Wrangell-St. Elias National Park and Preserve (Figure 1) [21]. In interior Alaska, one dinosaur track has been identified in fluvial strata exposed in the Charley River thrust belt within Yukon-Charley Rivers National Preserve [22]. Finally, the Prince Creek Formation exposed on the North Slope of Alaska has yielded a rich assemblage of dinosaur and bird trackways e.g., Refs. [23][24][25].  Figure 2 shows the location of the study area; this area is enlarged in Figure 2.  . Box labeled Figure 2 shows the location of the study area; this area is enlarged in Figure 2.
Workers have also begun to use various proxies archived in the sedimentary record to search for regional links between high-latitude, Late Cretaceous Arctic tracksites and broader climate patterns e.g., Refs. [7,26]. In addition to the Cantwell Formation, terrestrial assemblages that have been the focus of detailed sedimentological and paleoclimate studies include the Prince Creek Formation (North Slope of Alaska) and the Chignik Formation (southern Alaska) e.g., Refs. [5,[7][8][9]27]. Results of these studies suggest that differences in climate over~20 • of paleo-latitude likely served as a key driver of differences observed in large-bodied herbivore populations [7]. Initial studies of paleoclimate have also begun to demonstrate links to more global climate patterns during the late Campanian-early Maastrichtian. For example, δ 13 C isotope analysis of organic material suggests that a mid-Maastrichtian warming event recognized globally is also recorded in the Cantwell Formation [8].

Existing Structural and Stratigraphic Framework 2.2.1. Age Control
The age of the Cantwell Formation was originally determined using plant fossils, which resulted in conflicting age determinations ranging from Cretaceous to Paleocene [28][29][30][31]. For several decades, the accepted age for the Cantwell Formation was Paleocene based on paleobotanical data [30]. More refined age constraints now include palynological data with stratigraphic control [2] and 14 bentonites dated using U-Pb Laser Ablation Inductively Coupled Plasma Mass Spectrometry (LA-ICP-MS) [3,8,11,18]. Ridgway et al. [2] reported results of palynologic analyses from 135 mudstone samples that yielded an age range of late Campanian-early Maastrichtian for the Cantwell Formation. Recently dated bentonites yielded ages between 88.0 ± 0.4 Ma and 67.8 ± 0.2 Ma [3,8,11,18]. Seven locations for which published bentonites are reported are within our general study area, though localities with bentonite ages remain difficult to correlate with measured stratigraphic sections from both this study and previously published studies due to the structural complexity e.g., Ref. [11].

Previous Stratigraphic Analyses of the Cantwell Basin
Initial geologic mapping and stratigraphic studies of the Cantwell basin were part of regional analyses primarily focused on the area east of the Nenana River ( Figure 2) [1,31,32]. The first workers to classify the Cantwell basin interpreted the basin as a Paleogene, strikeslip pull-apart basin related to the Hines Creek and Denali fault systems e.g., Ref. [1]. Ridgway et al. [2] later re-interpreted the Cantwell basin as a Late Cretaceous thrust-top basin related to the final stages of Mesozoic collisional tectonics. The reconstruction of depositional environments from isolated exposures both east and west of the Nenana River led to the development of a basin model by Ridgway et al. [2], in which alluvial fan systems merged into axial braided fluvial systems that drained the basin to the east. Extensive lacustrine deposits were thought to be most prevalent along the southern basin margin, consistent with an asymmetric, southward-dipping basin floor. Marginal marine influence near the top of the Cantwell Formation is suggested by the presence of marginal marine dinoflagellates, and oncolitic and algal-laminated limestone beds at Double Mountain (DM1 in Figure 2) [2]. It is important to note that these early studies included limited data from exposures of the Cantwell Formation located east of the Nenana River ( Figure 2). The exposures in that region are included in regional mapping studies [33,34], but have not been integrated into detailed stratigraphic analyses of the Cantwell basin, including this study. Finally, the most recent stratigraphic studies have focused on the Tattler Creek/Sable Mountain area (Figure 2), where robust paleobotanical and paleontological datasets have been collected with stratigraphic control and integrated with interpretations of alluvialfluvial depositional environments e.g., Refs. [6,9]. Preserve. The stratigraphic chart shows the Middle-Jurassic-Cenozoic stratigraphy of the study area (modified from Trop et al. [11]). The areal extent of this map is shown by the dashed box in Figure 1. The locations of measured sections referenced in Figures 3 and 6 are labelled using the abbreviations shown in the key. Geologic mapping data are a compilation of new data described in the text and the data of Reed [35], Gilbert and Redman [36], Csejtey et al. [33,34], and Cole et al. [10]. The Cantwell Formation is subdivided into the Polychrome member (Kcp-lower member) and the Cabin Creek-Double Mountain member (Kcc-upper member).

Tectonic Setting of the Central Alaska Range
The Cantwell basin formed in a region of polyphase deformation dominated by collisional tectonics from the Late Jurassic through the Cretaceous. This region has been termed the Alaska Range Suture Zone (ARSZ), with the present-day surface expression defined as the region between the Talkeetna fault to the south and the Hines Creek fault to the north (Figure 1) e.g., Refs. [11,[37][38][39][40][41]. The ARSZ has been studied extensively and is divided into three components: (1) the ancestral North American margin to the north, which served as the Mesozoic collisional backstop (Yukon composite terrane in Figure 1); (2) the oceanic Wrangellia composite terrane (WCT) to the south (Figure 1), which is thought to have collided obliquely from south to north along the continental margin [38,42,43]; and (3) deformed and metamorphosed sedimentary basin assemblages that represent the Jurassic-Cretaceous ocean basin between the ancestral North American margin and the WCT that was subsequently segmented and eventually fully closed in the collisional zone (e.g., Kahiltna basin) [43,44]. The part of the Cantwell basin exposed in our study area also overlaps a slightly older, Late Jurassic-Early Cretaceous accretionary zone associated with the interpreted collision of the allochthonous Farewell terrane with the ancestral North American margin [45]. Much of the Cantwell Formation in our study area was deposited unconformably atop Upper Triassic basalt interpreted as part of the Farewell terrane (unit TRbd in Figure 2) e.g., Refs. [35,36,45], whereas along the northern basin margin, the Cantwell Formation overlies deformed Mesozoic marine strata interpreted to be associated with the ancestral North American margin (unit TRcs in Figure 2) [45].

Stratigraphic Nomenclature
Early studies of the Cantwell Formation typically used the terminology "upper" and "lower" Cantwell Formation to refer to the Upper Paleocene-Lower Eocene Teklanika volcanic strata and the Upper Cretaceous siliciclastic Cantwell Formation, respectively e.g., Refs. [2,30,34]. More recently, the long-standing name Teklanika Formation [46] has been used for the volcanic sequence, while "Cantwell Formation" refers only to the sedimentary strata e.g., Ref. [11]. We conform to this terminology as separate formations are warranted by distinct lithologies, a depositional hiatus, and an angular unconformity between the Cantwell Formation and the Teklanika Formation volcanic strata in some parts of the basin. For the purposes of our stratigraphic analysis, we further divide the Cantwell Formation into two informal members. We use the term Polychrome member to refer to a lower member and the term Cabin Creek-Double Mountain member to refer to an upper member. The sedimentological basis for these members is described in the Results section.

Stratigraphic Analysis
Our stratigraphic analysis was based on the construction of a composite representative section for the Cantwell Formation ( Figure 3). We identified several exposures between the Toklat River and the East Fork of the Toklat River that comprise common lithofacies of the Cantwell Formation ( Figure 2). Due to the structural complexity of the basin, this composite representative section does not perfectly record the stratigraphic thickness of the Cantwell Formation. We accounted for these challenges by completing our own mapping in select parts of the study area to develop as complete a stratigraphic framework as possible, particularly where measured sections were separated by faults. Stratigraphic data were collected using the standard methods of measuring a section with a Jacob's staff. Data collected from stratigraphic sections included bedding thickness, bedding geometry, grain size, sedimentary structures, color, compositional clast counts of conglomerate, maximum particle size measurements, and the documentation of fossils. The latitude and longitude of the base and top of measured sections that comprised our composite representative section are available in Supplementary File S1.

Geologic Mapping
Geologic mapping was conducted over two field seasons as part of a USGS EDMAP cooperative project. The bulk of our mapping was conducted along the Toklat River corridor and in the region between the Toklat River and the East Fork of the Toklat River north of the Denali Park Road (Figure 2). In this paper, we present results of our mapping compiled with published mapping data of Reed [35], Gilbert and Redman [36], Csejtey et al. [34], and Cole et al. [10], and previously compiled data of Wilson et al. [48]. The data compilation was completed using ESRI ArcMap version 10.6.2. Our compilation also merged multiple generations of mapping along the Healy/Mt. McKinley quadrangle boundary, which lies just east of the Toklat River.

U-Pb Detrital Zircon Geochronology
All detrital zircon samples were crushed, separated, and mounted at the Arizona LaserChron Center via standardized methods of physical and chemical separation [49][50][51][52]. SEM-generated cathodoluminescence (CL) images were used to select grains without complex crystallographic domains, and back-scattered electron (BSE) images were used to avoid grains that were not zircon. Grains were selected randomly without intentional bias toward size or shape, and fractures were avoided to minimize the introduction of common Pb. All U-Pb analyses were conducted using laser ablation inductively coupled mass spectrometry (LA-ICPMS). Some analyses were conducted using a Thermo Element2 single collector ICPMS and some samples were analyzed in 2019 and 2020 using a Nu Plasma multicollector ICPMS [51]. Analyses conducted on the Nu Plasma multicollector ICPMS utilized rapid acquisition (3s/analysis) protocols as outlined by Sundell et al. [53]. A spot diameter of 20-30 µm was used for all analyses. Data reduction for analyses conducted on the Element2 was completed using "AgeCalc" protocols; in-house protocols developed by the ALC. Data reduction for analyses conducted on the Nu Plasma Multicollector was completed using AgeCalcML [53]. Both reductions include a common Pb correction (cf. [54]). Probability density plots of detrital zircon data were generated using Isoplot 4.15 [55], and age peaks for detrital zircon spectra were calculated using Age Pick software [49]. Sample locations are available in Supplementary File S1, probability density plots by sample are available in Supplementary File S2, and U-Pb data are available in Supplementary File S3.
Analysis of detrital zircon data in our study was utilized primarily for provenance interpretations. We do not report maximum depositional ages for our detrital zircon geochronology dataset since, as has been recognized by previous workers, Pb loss is likely prevalent in zircons from the Cantwell Formation [11]. Anomalously young ages that comprise young age tails are present in published U-Pb zircon data from the Cantwell basin [11], as well as our detrital zircon data. Given that these anomalously young ages consistently fall within the time interval of peak magmatism associated with the Teklanika Formation at~57 Ma [10], it is likely that leaching of Pb occurred in the presence of hydrothermal fluids at this time [56].

Stratigraphic Framework
We divide the Cantwell Formation into two informal stratigraphic members: a lower member defined as the Polychrome member, and an upper member defined as the Cabin Creek-Double Mountain member. Figure 3 shows detailed sedimentological data for our composite representative section; Figures 4 and 5 illustrate the distinctive sedimentological features of each member; and Figure 6 shows lithostratigraphic correlations by member between measured sections that comprise our composite section and other published measured sections.      . Stratigraphic cross-section correlating measured sections from this study with previously published measured sections (e.g., Ridgway et al. [2]; Trop [47]). The gradational contact between the Polychrome and Cabin Creek-Double Mountain members is used as a datum. The Cabin Divide fault (CDF) truncates stratigraphy along the northern basin margin and section MS1 is structurally isolated and in unconformable contact with underlying Upper Triassic submarine strata (TRcs).

Representative Locations and Measured Sections
The lower part of the Polychrome member is best exposed on the southern and eastern flanks of Polychrome Mountain, and the upper part is best exposed in the Tattler Creek/Sable Pass area and along the East Fork of the Toklat River ( Figure 2). The Polychrome member is represented by our measured sections PC1L, RC1, and EF2, and also PC1 of Trop [47] (Figure 3). The base of section PC1L corresponds to an unconformable depositional contact between the Cantwell Formation and Upper Triassic basalt (TRbd in Figure 4A).

Sedimentology
The Polychrome member fines gradationally upward and is characterized in the lowermost 180 m (Section PC1L in Figure 3) by clast-supported, cobble-boulder conglomerate with a sandy matrix and interbedded massive, coarse-grained sandstone. Conglomerate clasts are well-rounded but poorly organized, lacking clear imbrication. Individual conglomerate beds are several meters thick, laterally continuous, and sometimes include localized lenses of coarse-grained sandstone. The average maximum particle size is 11.7 cm, though some clasts recorded are as large as 28 cm. Coarse-very-coarse-grained, massive sandstone intervals range in thickness from 1 to 4 m and comprise less than 10% of the total stratigraphic thickness in section PC1L ( Figure 3).
The middle~1200 m of the Polychrome member is characterized by interbedded pebble-cobble conglomerate, medium-grained sandstone, and minor siltstone ( Figure 3). This part of the stratigraphy generally continues to fine upward, though some cyclicity exists between conglomerate-dominated and sandstone-dominated packages within a stratigraphic section of several hundred meters. In conglomerate-dominated successions such as in much of section PC1 (Figure 3), conglomerate beds range from one to several meters thick and are laterally continuous with tabular geometries ( Figure 4B). Conglomerate beds are clast-supported and well organized, and clast imbrication is common. The average maximum particle size is 6 cm (Section PC1 of Trop, [47] and Ridgway et al. [2]). Sheet conglomerates are interbedded with lenses of medium-coarse-grained, trough crossstratified sandstone, massive sandstone, and minor siltstone. Multiple upright tree casts were observed where pebble-cobble conglomerate is concentrated around the base of the fossilized trees, and interbedded pebble conglomerate, sandstone, and siltstone encases the upper part of the trunks ( Figure 4C).
Sandstone-dominated intervals, best represented by our measured section RC1, are characterized by medium-coarse-grained sandstone beds with abundant internal scouring and planar and trough cross-stratification ( Figure 3). Granule-pebble lags are common at the base of scour surfaces and transported fossilized plant fragments are abundant. Shallow, lenticular bodies fine upward from granule-pebble conglomerate at the base to massive, medium-grained sandstone at the top. The average maximum particle size of conglomerate beds is 3.5 cm, and the thickness of these bodies rarely exceeds three meters. Interbedded massive and ripple cross-stratified siltstone and fine-grained sandstone beds display abundant mottling and other pedogenic structures, and thin, organic-rich mudstone beds are sometimes present. Fossils most commonly include leaf imprints in fine-grained beds, though a relatively small, likely juvenile theropod trackway was found on the base of one coarse-grained sandstone bed ( Figure 4E).
Measured section EF2, along the west bank of the East Fork of the Toklat River, documents lithofacies diagnostic of the uppermost Polychrome member (Figure 3). The maximum thickness of this part of the stratigraphy is uncertain, but our observations at section EF2 and along Tattler Creek suggest that at least 450 m of similar strata are present. Section EF2 is characterized by three distinctive lithofacies separated by sharp contacts ( Figure 4F-H): (1) Medium-coarse-grained, massive and trough cross-stratified sandstone with abundant internal scouring (Sm in Figure 4F). Sandstone is internally organized in broadly lenticular bodies several meters thick with individual bedding thicknesses of 10-30 cm and granule conglomerate lags associated with trough cross-sets within 20-50 cm thick beds. (2) Nodular, very thinly bedded carbonaceous siltstone and shale with lignite breaks. This lithofacies comprises intervals approximately 2-5 m thick ( Figure 4F-H).
(3) Fine-grained, ripple cross-stratified sandstone with abundant plant debris and mud drapes in tabular beds ( Figure 4G). Individual beds are less than 10 cm thick and comprise intervals that are 1-3 m thick ( Figure 4F-H). The transition between siltstone and sandstone is typically gradational ( Figure 4G).

Representative Locations and Measured Sections
The Cabin Creek-Double Mountain member is well exposed in the Cabin Divide area, located immediately south of the Wyoming Hills (Figure 2), and at Double Mountain (DM1 in Figure 2). This mudstone-rich member is highly folded in the Cabin Creek area but is much less deformed in the Double Mountain area. Measured sections CAB01 (this study) and DM1 (Ridgway et al. [2]) serve as well-exposed examples of characteristic lithofacies, but do not record the total stratigraphic thickness of the member due to deformation and poor exposure in the more mudstone-rich sections.

Sedimentology
The Cabin Creek-Double Mountain member is characterized by thick sections of blackgray mudstone interbedded with sandstone packages of varying thickness ( Figure 5A). All beds are laterally continuous and easily traced for hundreds of meters ( Figure 5A). Sandstone beds occur as both 1-3 m thick packages interbedded within predominantly mudstone intervals, and as 5-7 m thick packages of more amalgamated sandstone with minor interbedded mudstone ( Figure 5A,B). The thicker, more amalgamated sandstone packages consist internally of 10-30 cm thick beds of massive and ripple cross-stratified, medium-grained sandstone interbedded with massive and laminated mudstone intervals averaging~0.5-1 m thick ( Figure 5B). In the 1-3 m thick sandstone packages interbedded within intervals dominated by mudstone, bedding thickness diminishes from base to top, and the ratio of sand to shale decreases as shale breaks become more common ( Figure 5C). In some sandstone packages, low-angle lateral accretion surfaces (see white triangles in Figure 5C) and minor internal scouring is observed. Broad, shallow erosional scours are most common in the thick, amalgamated sandstone packages.
Sandstone packages in the Cabin Creek-Double Mountain member consist of individual beds that are fine-medium-grained and massive at the base, and grade upward into wavy, flaser, and less commonly lenticular bedding ( Figure 5D). In some cases, sandstone beds are bioturbated ( Figure 5D) and contain soft-sediment deformational structures. The tops of most sandstone beds are dominated by symmetric ripples (Figure 5D,E); in rare cases, interference ripples are also preserved ( Figure 5G). Fossilized wood fragments are commonly concentrated at the base of sandstone-dominated packages with minor coalified fragments interspersed ( Figure 5H). Well-preserved leaf imprints are common in mudstone beds of any thickness.
Cyclicity between sandstone-and mudstone-dominated intervals is observed over the thickest exposed sections of the Cabin Creek-Double Mountain member. In Figure 5I, these cycles are bounded by~5-7 m thick, amalgamated sandstone packages that stand out as the most resistant features of the exposure. The intervening thinner-bedded intervals have a gross thickness of~15-25 m, are dominated by mudstone, and contain within them 1-3 m thick sandstone packages such as the package shown in Figure 5C. Other distinctive lithofacies in the Cabin Creek-Double Mountain member are thin (<30 cm) limestone beds. Two types of limestone have been documented. The first type is a laminated limestone ( Figure 5J). Individual laminae are crenulated and wispy with an individual thickness of less than 5 cm. The second type of limestone contains oncolites at the base and laminated limestone at the top of the bed ( Figure 5K).

Correlation of Published Measured Sections
Here, we use the stratigraphic framework derived from our composite representative section of the Cantwell Formation to establish lithofacies correlations with previously published stratigraphic sections ( Figure 6) (Trop [47]; Ridgway et al. [2]). We use the gradational contact between the Polychrome and Cabin Creek-Double Mountain members as a datum for comparison between sections. Note that the contact between the overlying Teklanika Formation and the Cantwell Formation is not a consistent datum for correlations as this contact is unconformable ( Figure 6).

Correlation of the Cabin Creek-Double Mountain Member
We correlate the upper~1000-1200 m of the previously published measured section PA1/PA2 at Panorama Mountain with the Cabin Creek-Double Mountain member ( Figure 6). The upper part of section PA1/PA2 displays a prominent shift to 10-40 m thick sandstone packages interbedded with 10-50 m thick mudstone packages (FA3 and FA5 of Ridgway et al. [2]). In contrast to sections CAB01 and DM1, sandstone beds are commonly trough cross-stratified and mudstone packages are finely laminated, described by Ridgway et al. [2] as 'varve-like'. Despite these differences, we correlate this part of the stratigraphy at Panorama Mountain with the Cabin Creek-Double Mountain member because of the presence of thick mudstone packages and a clear lithofacies trend throughout section PA1/PA2 that we interpret to represent the gradational transition between members. We interpret conglomerate-dominated intervals in the lower Polychrome member to represent deposition in fluvial-dominated or 'wet' alluvial fan systems [57,58]. Tabular, clast-supported conglomerate beds ( Figure 4B) suggest deposition via sheet flood processes, where sediment is dispersed through unconfined flow from the mouths of proximal channels on alluvial fans e.g., Refs. [59,60]. Large average clast size (Figure 3) suggests an alluvial surface that was altered periodically by high-velocity flood events. Rapid sedimentation events on vegetated fan surfaces are supported by the documentation of upright tree fossils ( Figure 4C), implying the deposition of several meters of sediment relatively rapidly e.g., Ref. [61]. While conglomerate beds deposited via sheet floods are commonly imbricated [59], we interpret the general lack of both organization and imbrication of conglomerate in the lowest 180 m (PC1L in Figure 3) to suggest a more proximal fan position, where coarse sediment was transported in flashy, high-energy events (PC1L in Figure 7A) e.g., Ref. [62]. This interpretation is further supported by boulder-size clasts in the Polychrome member at PC1L.

Middle Polychrome Member-Braided Fluvial System
We interpret sandstone-dominated intervals characterized by trough cross-stratified sandstone and lenticular granule-pebble conglomerate beds with horizontal stratification to represent axial, braided-fluvial systems in the Cantwell basin (RC1 in Figure 7A). Conglomeratic lenticular bodies represent traction-based transport in active channels and on longitudinal barforms, whereas planar and trough cross-stratified sandstone represents the deposition of finer sediment on transverse barforms [63]. Thin-bedded, massive and horizontally stratified sandstone and interbedded siltstone are interpreted to represent overbank and floodplain deposits. Abundant leaf fossils suggest that subaerially exposed parts of bars and overbank environments were vegetated. Paleoflow indicators from the northern and southern basin margins [2,47] suggest that alluvial fans represented by the lower-middle Polychrome member were sourced from high topography on both sides of the basin, and that these systems merged into axial fluvial systems that drained the basin to the east, consistent with the depositional model of Ridgway et al. [2].

Upper Polychrome Member-Coastal Fluvial System
We interpret the upper Polychrome member to represent a transition to a coastal fluvial system with a vegetated floodplain ( Figure 7B). Broadly channelized, mediumcoarse-grained sandstone bodies suggest deposition within well-established channels with the tractive transport of sediment ( Figure 4F). Tabular, thin-bedded ripple crossstratified sandstone beds ( Figure 4G,H) are interpreted to represent overbank deposits proximal to channels. Finally, thick carbonaceous siltstone and shale intervals ( Figure 4F,H) are interpreted to represent a marshy overbank environment e.g., Ref. [64]. Thin lignite beds less than 5 cm thick indicate times of accumulation of organic matter atop a poorly developed soil profile and are often reworked as transported peat mats e.g., Ref. [65].

Cabin Creek-Double Mountain Member
We interpret the Cabin Creek-Double Mountain member to represent a tidally influenced marginal marine depositional system that included estuarine and coastal lacustrine environments ( Figure 7B,C). Note that our interpretation of an estuarine environment does not necessarily invoke a relationship with an incised paleovalley and related embayment, rather the use of the term is intended to refer generally to a transgressive estuarine system e.g., Ref. [66]. The Cabin Creek-Double Mountain member is characterized by heterolithic intervals dominated by wavy and flaser bedding, mud drapes on internal beds within sandstone packages, and symmetric ripples. The flaser-bedded sandstone, defined by discontinuous mud drapes in ripple cross-stratified sandstone, and the wavy-bedded sandstone ( Figure 5D), defined by continuous mud drapes in ripple cross-stratified sandstone, both reflect abrupt changes in flow velocity that allow the deposition of markedly different grain sizes. These types of sedimentary structures are common in modern estuarine depositional systems within intertidal zones where ripple cross-stratified sandstone is deposited during the maximum flow velocities associated with dominant tidal periods, whereas silt and mud settle out of suspension during slack-water periods [67,68]. These diagnostic sedimentary structures are more common at Double Mountain than Cabin Creek, which we interpret to reflect a more proximal part of the estuarine system preserved at Cabin Creek relative to Double Mountain. In many estuarine systems, tidal currents are often strongest in medial positions where water depths remain shallow, but the influence of fluvial processes and direct sediment input is minimal [66]. Symmetric ripples are abundant in all sandstone packages ( Figure 5E,F) and are interpreted as a product of bi-directional tidal currents re-working sand deposited in distributary channels and as part of tidal bars. Interference ripples ( Figure 5G) are far less common and most likely formed in specific locations where subordinate currents disrupted the dominant bi-directional current, overprinting the primary symmetric ripples that are ubiquitous throughout the section e.g., Ref. [69].
The 5-7 m thick amalgamated sandstone packages consist of mostly massive and ripple cross-stratified sandstone beds with broad internal scours and are interpreted to represent primary distributary channels within the estuarine system. When active, these distributary channels were where most of the sand in the system was deposited and frequently reworked e.g., Ref. [66]. The cyclicity that is highlighted by the occurrence of an amalgamated sandstone package every~25 m ( Figure 5I) is interpreted as the product of avulsion processes where the depositional axis migrated laterally across the estuary. Low-angle lateral accretion surfaces were observed in sandstone packages of varying thickness and are interpreted as evidence of laterally migrating point bars formed within these constantly evolving distributary channels [66]. The 1-3 m thick sandstone packages within mudstone-dominated intervals are interpreted to represent off-axis parts of the system where shallow, short-lived distributary channels formed when sand was routed from the main depositional axis, and mudstone was deposited at times when there was an absence of sand input and corresponding lower energy.
At Cabin Creek, there is an absence of obvious trace fossils, but burrows are common at the Double Mountain section. It is possible that these features were not well preserved in the Cabin Creek area due to the extensive deformation that these strata have experienced in the footwall of the Cabin Creek fault (see Section 6). Alternatively, the absence of burrowing organisms may reflect particularly harsh conditions in the depositional environment during the deposition of the Cabin Creek section. The fluvial-tidal transition is often a particularly challenging environment for biological life to thrive due to frequent re-working of the sediment substrate and variable salinity, including seasonal variations e.g., Ref. [66]. The more common bioturbation in the Double Mountain section may be a function of more stable salinity levels in what we interpret to have been a more medial part of the estuarine system. This interpretation would be consistent with the documentation of several species of dinoflagellates, suggestive of more brackish marine conditions at Double Mountain [2]. The thin limestone beds best developed at Double Mountain are interpreted as algalinfluenced limestones. The laminated limestone beds with crenulated lamina ( Figure 5J) are interpreted to represent amalgamated algal mats e.g., Refs. [70][71][72]. The oncolitic limestones ( Figure 5K) represent deposition when the algal material accumulated around a nucleus that was periodically transported in the intertidal zone e.g., Ref. [73].
Like modern estuarine-dominated coastlines, we interpret the Late Cretaceous estuarine systems of the Cantwell Formation to have been laterally adjacent to lacustrine systems. Sedimentologic data for lacustrine depositional systems is best documented from the previously published measured section at Panorama Mountain ( Figure 6) [2]. At this location, interbedded with the heterolithic intervals are 10-50 m thick packages of laminated mudstone with millimeter-scale laminations rich in plant fragments. These types of sedimentary structures reflect an exceptionally calm environment within which suspension deposition dominated for long periods of time ( Figure 7C). Trough cross-stratified sandstone bodies associated with the laminated mudstone packages likely represent times when coarse sediment prograded into lacustrine environments from the outlets of fluvial channels, possibly within small, prograding deltas e.g., Ref. [74].
The basin-scale shift from alluvial-fluvial systems in the Polychrome member to marginal marine and lacustrine systems in the Cabin Creek-Double Mountain member represents a long-lived transgression. Ridgway et al. [2] interpreted marginal marine influence in the uppermost Cantwell Formation to be part of the Bearpaw transgression. The Bearpaw transgression has been documented in several locations in North America as a diachronous transgression throughout the North American Cordillera and other parts of Late Cretaceous Euramerica e.g., Refs. [2,[75][76][77]. Overall, our depositional model remains consistent with this general interpretation, though more robust chronostratigraphic control would be required to relate the specific timing of transgression in the Cantwell basin to other stratigraphic sections in the Cordilleran foreland basin system where the Bearpaw transgression has been better documented. An alternative, more local driver of the transgression recorded by the Cantwell Formation may be increased subsidence resulting from loading of the crust by syndepositional thrust systems that bordered the basin.

Stratigraphic Context of Trace Fossils in the Cantwell Formation
Our depositional model also provides stratigraphic context to important paleontological datasets, especially bird and dinosaur trackways in Denali National Park and Preserve. The most prolific dinosaur tracksites described in the literature are associated with fluvial strata of the middle-upper Polychrome member, most notably in the Tattler/Igloo Creek region (Figures 2 and 7A,B). This part of the basin has yielded numerous theropods and hadrosaur trackways, the first pterosaur manus trackway found in Alaska, and the northernmost occurrence of a therizinosaur e.g., Refs. [14,16,17,78]. Also of note are numerous documented trackways and feeding traces of wading birds (Aquatilavipes swiboldae, Ignotornis mcconnelli, Magnoavipes denaliensis, Gruipeda vegrandiunus, and Uhangrichnus chuni), which are most common in strata of the uppermost Polychrome member and the lower Cabin Creek-Double Mountain member ( Figure 7B) [4,18,78]. This record of wading birds is important because it provides insight into the paleobiology of the coastal fluvial and estuarine/lacustrine depositional systems that formed through transgression in the Cantwell basin. Similar wading birds in modern estuarine and tidal flat environments often take advantage of low tide to feed in the intertidal zone. Trackways in these environments have especially high preservation potential as new sediment is deposited rapidly by the next high tide e.g., Ref. [79].

Geologic Mapping and Structural Observations
Our mapping and compilation work was focused primarily on the Toklat River corridor and the northern margin of the Cantwell basin between the Toklat River and the East Fork of the Toklat River ( Figure 2). The results of this compilation and our strike and dip measurements along the Toklat River corridor are shown in Figure 8A. The inset map shown in Figure 8B is a detailed map of the northern margin of the study area. Numbered locations in both maps are referenced in the text to indicate the locations where key structural observations were made in the field.

Thrust System along the Northern Basin Margin
The Cantwell Formation along the northern basin margin is deformed within a previously undocumented, south-verging thrust system (Figure 8). The Cabin Divide fault (CDF) marks the northern margin of the continuous exposures of the Cantwell Formation within this system and juxtaposes Jurassic-Cretaceous metasedimentary strata (unit KJf) against Upper Cretaceous strata of the Cantwell Formation (Location 1, Figures 8A,B and  9). The dip of the CDF is difficult to constrain since the fault is mostly covered in a valley. Thus, we did not alter the trace of the fault as shown in the original linework of Gilbert and Redman [36], which depicts the fault as a steeply south-dipping normal or oblique-slip fault. However, based on the following structural characteristics of the Cantwell Formation just to the south, it is apparent that the CDF must have accommodated south-verging, contractional deformation, at least during the deformation of the Cantwell Formation. Open folds in the Polychrome member south of the CDF verge to the south-southeast with approximately east-west trending axial traces (e.g., Location 2 in Figure 8B; folds also shown in a field photo in Figure 8C). Reverse faults associated with these folds dip approximately 40-45 • to the north and typically break through anticlines slightly fore of the fold hinge. Based on this observation and the fold geometries, we interpret the folds that deform the Polychrome member as fault propagation folds. Upper Paleocene-Lower Eocene volcanic strata of the Teklanika Formation unconformably overlie folded strata of the Cantwell Formation within the northern margin thrust belt and are most extensively preserved in synclines ( Figure 8C). Folds in the Teklanika Formation are lower amplitude and disharmonic relative to folds in the Cantwell Formation (Location 3 in Figures 8B,C  and 9B).
In the shale-rich Cabin Creek-Double Mountain member, fold geometries range from open to tight with angular hinges, including chevron-style folds. Local thickness changes are common in deformed shale beds, particularly near and within fold hinges. The axial planes of folds in this part of the Cantwell Formation are typically subvertical and trend roughly east-west. While the difference in deformational style observed between the shalerich and sandstone/conglomerate-rich strata in the Cantwell Formation may be attributable to lithology, it is also possible that folds in the upper part of the section record a component of detachment folding e.g., Ref. [80], though we cannot fully assess this option based on our data.    Figure 8C is marked on this cross-section. (C) Geologic cross-section from C-C . The frontal reverse fault of the south-verging system places the Polychrome member (Kcp) over the Teklanika Formation (Pgtv) (Location 7). The Park Road fault (PRF) passes through Sable Pass with Neogene strata exposed south of the fault. The location of the field photographs shown in Figure 8D,E are marked on this cross-section.

Thrust System along the Southern Basin Margin
The southern basin margin is defined by a north-verging thrust system involving Ordovician-Triassic marine strata for tens of kilometers along strike (units DOs, Dl, TRPs, and TRbd, Figures 2 and 8A). North-verging thrust sheets are particularly well exposed along the upper Toklat River corridor ( Figure 8A). A key structure associated with this system is a north-south-trending tear fault along the east bank of the upper Toklat River (Location 4 in Figure 8A). The leading reverse fault along the northern edge of the thrust sheet bounded by the tear fault places Upper Triassic basalt (TRbd) over the Polychrome member (Kcp) (Location 5 in Figures 8A and 9A). The tear fault itself truncates a depositional contact between the Teklanika Formation (Upper Paleocene-Lower Eocene) and unit TRbd at location 4 ( Figure 8A). This relationship contrasts with location 6, just 5 km to the east, where the Cantwell Formation is in fault contact with unit TRbd ( Figure 8A).

Interpretation of Deformational Timing
Syn-and post-depositional contraction influenced alluvial systems along the basin margins and resulted in the deformation of the Cantwell Formation during and soon after deposition. During the deposition of the Cantwell Formation, high topography within the active, basin-bounding thrust belts provided proximal sources of sediment such as quartz and schist clasts along the northern basin margin and diagnostic limestone clasts along the southern basin margin. Basinward propagation of the thrust belts early in the deposition of the Cantwell Formation is evident from the presence of eroded, lithic-rich coarse-grained sandstone of the lowermost Polychrome member re-deposited as clasts in conglomerate beds of the upper-lower to middle Polychrome member, cf. [81]. Ridgway et al. [2] also documented growth unconformities along the southern basin margin near Panorama Mountain (PA1/2 in Figure 2), which clearly record this process of the syn-kinematic deposition of alluvial fan strata along the active margins of the basin. Additional evidence of the basinward propagation of the northern margin thrust belt is the structural isolation of part of the lower Polychrome member at Mt. Sheldon ( Figure 8A). In this location, the basal cobble-boulder conglomerate of the Polychrome member unconformably overlies Upper Triassic marine strata (TRcs in Figures 2 and 8A) and is folded in the core of an open syncline. This exposure is isolated to the north of the CDF, though the timing of its isolation is unclear as younger strata of the Cantwell Formation may have been deposited and subsequently eroded as the northern basin margin was uplifted.
The volcanic eruptions that produced the Teklanika Formation had commenced by the Late Paleocene [10,46] and coincided with the continued shortening of the Cantwell basin. The contact between the Cantwell Formation and Teklanika Formation is locally an angular unconformity, and the Teklanika Formation is deformed by lower amplitude, more open folds relative to the Cantwell Formation. The most complete sections of the volcanic strata are preserved in synclines in the Cantwell Formation. We interpret these relationships to reflect the syn-deformational deposition of the Teklanika Formation following significant unroofing of parts of the Cantwell Formation, such as at location 3 where the Teklanika Formation is deposited over part of the middle-lower Polychrome Member ( Figure 8B,C). Further evidence for contractional deformation following deposition of the Teklanika Formation includes two key field relationships: (1) at location 7, an east-west-trending, south-verging reverse fault at the frontal edge of the northern thrust belt places the Cantwell Formation over the Teklanika Formation ( Figures 8B,E and 9C); and (2) at location 4, the tear fault along the upper Toklat River corridor truncates the contact between the Teklanika Formation and underlying Triassic basalt ( Figure 8A). Though not the focus of this study, evidence for the post-Eocene overprinting of contractional structures in the Cantwell basin is abundant. One of the better examples is the Park Road fault (PRF), first documented and informally named by Gilbert and Redman [36]. This fault roughly coincides with the location of the Park Road through the middle third of the road, truncating the frontal edge of the south-verging northern margin thrust belt ( Figure 8A). This structure appears to behave like well-described strike-slip faults that often change along strike from extension to strike-slip to contraction over relatively short distances e.g., Ref. [82].

U-Pb Detrital Zircon Geochronology Results
Analysis of detrital zircon populations for 13 samples from the Cantwell Formation allows for the reconstruction of regional sediment pathways into the basin and provides new insights into the paleogeographic setting of the Cantwell basin. The stratigraphic locations of the detrital zircon samples are shown in Figures 3 and 6. Since there is very little relative variation in the age spectra of samples from each respective member, we present composite probability density plots (PDPs) for each member in Figure 10 (Figure 10). Triassic-Jurassic grains comprise the most abundant Phanerozoic population with ages between ca. 240 and 150 Ma and account for 25% of all grains. Early and Late Cretaceous populations account for 4% and 8%, respectively (Figure 10). Paleozoic grains account for 15% of all grains and are dominated by a population  (Figure 10). Early and Late Cretaceous grains are more abundant relative to the Polychrome member, comprising 8% and 29%, respectively ( Figure 10). A total of 30% of grains yield Triassic-Jurassic ages, and 13% yield Paleozoic ages ( Figure 10). The Cabin Creek-Double Mountain member includes a smaller proportion of Precambrian grains (20%) relative to the Polychrome member with similar common populations. Our raw data include ten ages slightly younger than Late Cretaceous, which we attribute to Pb loss based on the correlation with the timing of peak magmatism at 57 Ma [10], along with no existing age control to suggest that any part of the Cantwell Formation was deposited after the Maastrichtian.

Detrital Zircon Provenance Interpretation
Comparison of detrital zircon ages from the Cantwell Formation with published ages and locations for potential primary (igneous) and secondary (meta)sedimentary sources yields the following insights: (1) The youngest population (ca. 77-66 Ma) overlaps with the deposition of the Cantwell Formation and was derived from a coeval magmatic arc system located south of the basin. (2) Cretaceous grains with ages that precede the deposition of the Cantwell basin were derived from older exhumed Cretaceous plutons exposed throughout south-central Alaska, primarily south of the Cantwell basin. (3) Sediment was recycled from Paleozoic and Mesozoic metasedimentary strata that were locally uplifted in syndepositional thrust systems along the southern and northern basin margins. (4) Sediment was recycled from Jurassic-Cretaceous sedimentary sources that were exhumed in a Cretaceous suture zone presently exposed south of the basin. (5) Precambrian grains were recycled from metamorphic rocks exposed north of the Cantwell basin that represent the ancestral North American margin.
For brevity, we have summarized our provenance interpretations in Table 1, along with key references.  [40], and references therein * Arc sources interpreted in the Cantwell Formation as recycled from the Kahiltna assemblage Note: Colors for primary sources correspond to colored bands in Figure 11 which represent these sources We also refer the reader to detailed compilations of regional detrital zircon sources documented for strata exposed in the central Alaska Range by Hampton et al. [42,83], Brennan and Ridgway [120], Lease et al. [121], and Romero et al. [40], and in the western Alaska Range by Box et al. [43].

Paleogeographic Insights from Provenance Interpretations
Late Cretaceous detrital zircon ages between ca. 75 and 66 Ma overlap with the Campanian-Maastrichtian biostratigraphic age ranges based on palynology of Ridgway et al. [2], and with U-Pb ages of bentonites interbedded within the Cantwell Formation e.g., Refs. [3,8,11,18]. This overlap suggests that there was an active silicic volcanic source supplying sediment to the basin during the deposition of the Cantwell Formation ( Figure 11). The coeval magmatic arc system was located south of the Cantwell basin (present day coordinates) and is represented by plutons exposed throughout the Alaska Range (Figures 1  and 12; Table 1). The youngest part of the arc consists of plutons mostly less than 75 Ma and is well exposed in the western Alaska Range, south of the Denali fault ( Figure 1) [83]. Restoration of~480 km of dextral displacement on the Denali fault [41,[122][123][124] places these plutons proximal to the southern margin of the Cantwell basin. While grains derived from the active arc are present in all samples from the Cantwell Formation (Figure 9), we observe a marked increase in the proportion of Late Cretaceous grains in the Cabin Creek-Double Mountain member of the Cantwell Formation. We interpret the increase in Late Cretaceous grains to record an uptick in magmatism, consistent with magmatic ages from the arc (Figure 11) [83]. North of the Cantwell basin, parts of the Yukon composite terrane (ancestral North American margin) comprised topographic highs during the Late Cretaceous, including a regional anticlinorium proximal to the northern margin of the basin ( Figure 12) [2,125]. We interpret significant sediment contribution from this region to the Cantwell basin. For example, all detrital zircon samples from the Cantwell Formation, regardless of geographic or stratigraphic location, yield a significant population between ca. 220 and 180 Ma (Figures 10 and 11). We interpret these grains to reflect plutonic sources located northeast of the Cantwell basin, the most proximal of which is the Taylor Mountain batholith (TMB) ( Table 1). The TMB intruded into the Yukon composite terrane, and published magmatic ages range from 216 to 181 Ma. This age range accounts for 23% of grains in the Cantwell Formation. Evidence for much more proximal sources of sediment from the YCT includes the presence of abundant quartz and schist clasts in conglomerate of the lower Polychrome member along the northern basin margin. These conglomerate beds are also associated with southward paleoflow indicators [2,81].  [83] and references therein). Latest Cretaceous detrital zircon ages overlap with a period of increasing magmatic intensity in the arc, which was located south of the Cantwell basin during the time of deposition. Colored bands delineate key time intervals that correspond with interpreted primary sources for detrital zircons in the Cantwell Formation. The gray band corresponds to the timing of regionally extensive magmatism represented in the study area by the Teklanika Formation. Abbreviated cardinal directions in parentheses indicate the paleo-position of each source relative to the Cantwell basin using modern coordinates. WAR-Western Alaska Range.
North of the Cantwell basin, parts of the Yukon composite terrane (ancestral North American margin) comprised topographic highs during the Late Cretaceous, including a regional anticlinorium proximal to the northern margin of the basin ( Figure 12) [2,125]. We interpret significant sediment contribution from this region to the Cantwell basin. For example, all detrital zircon samples from the Cantwell Formation, regardless of geographic or stratigraphic location, yield a significant population between ca. 220 and 180 Ma (Figures 10 and 11). We interpret these grains to reflect plutonic sources located northeast of the Cantwell basin, the most proximal of which is the Taylor Mountain batholith (TMB) ( Table 1). The TMB intruded into the Yukon composite terrane, and published magmatic ages range from 216 to 181 Ma. This age range accounts for 23% of grains in the Cantwell Formation. Evidence for much more proximal sources of sediment from the YCT includes the presence of abundant quartz and schist clasts in conglomerate of the lower Polychrome member along the northern basin margin. These conglomerate beds are also associated with southward paleoflow indicators [2,81].  Figures 11 and 12). The Kahiltna assemblage is exposed throughout the central and western Alaska Range, as well as the northern Talkeetna Mountains (Figure 1). Since most Early Cretaceous and older grains in the Cantwell Formation overlap with dominant populations in the Kahiltna Assemblage, it is difficult to accurately interpret its significance as a secondary sediment source versus both units being sourced from the same primary sources. However, there are minor, yet unique populations present in the Cantwell Formation that were almost certainly recycled from the Kahiltna assemblage. Specifically, the Talkeetna, Chisana, and Chitina arcs that were built on the Wrangellia composite terrane are the most likely sources for grains between~175 and 135 Ma in the Cantwell Formation (Table 1). From a paleogeographic perspective, it is most easily reconciled to have these specific zircons recycled from the more proximal Kahiltna assemblage rather than the more distal exposures of the plutons representing these three arc systems. It is important to note that reconstruction of~480 km of dextral displacement on the Denali fault system places the Kahiltna assemblage now exposed in the western Alaska Range directly south of the Cantwell basin.

Recycling Driven by Syndepositional Exhumation of the Cantwell Basin Margins
Minor detrital zircon populations between ca. 700 and 400 Ma (n = 112) and 1000 and 900 Ma (n = 20), observed primarily in the Polychrome member ( Figure 10), are unique to metasedimentary strata exposed in the syndepositional thrust sheets that border the Cantwell basin (Table 1; Figures 2, 8 and 9). Specific map units include Devonian-Ordovician and Permian-Middle Triassic marine strata (DOs and TRPs in Figure 2), Upper Triassic marine strata (TRcs in Figure 2), and Upper Jurassic-Lower Cretaceous marine strata (KJf in Figure 2) e.g., Refs. [33,34,127]. U-Pb detrital zircon ages from these strata include late Mesoproterozoic, Neoproterozoic, and Silurian-Ordovician populations that are extremely uncommon in rocks of the YCT or Wrangellia composite terrane [45]. Note that while unit KJf in the study area has been interpreted by regional mapping studies as part of the Kahiltna assemblage [34], recent U-Pb detrital zircon analyses from unit KJf along the northern basin margin yielded age spectra that are distinct from strata of the Kahiltna assemblage exposed south of the Denali fault [45]. Direct sediment input from basin-bounding thrust systems during the deposition of the Cantwell basin is also consistent with conglomerate clast compositions in alluvial fan strata along both margins documented in previous studies [2,81]. For example, limestone clasts common in conglomerate along the southern margin of the Cantwell basin can be directly linked to Upper Devonian sources (Dl in Figures 2 and 8A) in the southern thrust system based on conodont assemblages [2].

Regional Sediment Recycling from the Yukon Composite Terrane
Precambrian populations comprise significant components of the age spectra for the Polychrome member (48%) and the Cabin Creek-Double Mountain member (20%) (Figure 10). We interpret the more common populations of 1100-1000 Ma, 2000-1800 Ma, and 2800-2600 Ma to have been derived from the extensively exposed metasedimentary rocks of the Yukon composite terrane located north of the Cantwell basin ( Figure 1). Zircons of these ages have been well documented from primary and secondary sources across Laurentia (Table 1).

Regional Paleogeography of Mid-to Late Cretaceous Foreland Basin Depositional Systems: Next Steps
We interpret the Cantwell basin to have formed in a proximal part of a now-dismembered, regionally extensive Cretaceous foreland basin system that was bordered to the south by an active magmatic arc system ( Figure 12). During the deposition of the Cantwell Formation, defining components of the southern Alaska convergent margin located south of the paired foreland basin and arc system included marine and terrestrial depositional systems of the Matanuska forearc basin and the Chugach accretionary complex ( Figure 12) [83,[128][129][130][131][132]. Inboard of the magmatic arc, Mid-Late Cretaceous foreland basin depositional systems underwent a general transition from marine to non-marine depositional environments. The proximal terrestrial depocenters appear to have been spatially isolated and concentrated near the hinterland and magmatic arc ( Figure 12). These depocenters include the Beaver Lake formation in the eastern Alaska Range (mid-Cretaceous) [126], the Caribou Pass formation in the northern Talkeetna Mountains (mid-Cretaceous) [103], and possibly the Panorama Peak formation in the central Alaska Range (Upper Cretaceous) [133]. It is important to recognize that all Cretaceous terrestrial strata located south of the Denali fault system have been displaced along the fault away from their original depositional location.
An important unresolved issue is the stratigraphic and temporal connections between the Cantwell Formation and other Cretaceous foreland basin strata that cover a large part of interior Alaska. Eight to twelve kilometers of Cretaceous strata of the Kuskokwim Group are exposed in central and southwestern Alaska [134][135][136]. Most of the marine strata of the Kuskokwim Group range in age from Albian to Campanian [134]. A tuff in the Kuskokwim Group with an age of 77 Ma reported by Miller and Bundtzen [137], however, makes the upper part of this section temporally equivalent with the deposition of the Cantwell Formation. Roughly 200 km north of the Cantwell basin are exposures of the poorly dated and highly deformed Jurassic-Cretaceous strata of the Manley basin e.g., Refs. [138,139]. Northeast of the Cantwell basin, extensive Cretaceous strata of the Kandik basin have been incorporated into the Charley River thrust belt [140,141]. We tentatively interpret the Kuskokwim, Manley, and Kandik basins to represent more distal, largely marine parts of the Cretaceous foreland basin system that were connected to the Cantwell basin by the Maastrichtian but have subsequently been dismembered by Cenozoic strike-slip faults and regional shortening e.g., Refs. [110,135,136].
Future studies are required to better define the nature of the sedimentary relationships between the Cantwell basin and other Cretaceous strata in the interior of Alaska. Results of these studies will be important for refining the regional extent and paleogeographic configuration of the Cretaceous foreland basin system and interior seaways. Further integration of sedimentological and paleontological data from southern, interior, and northern Alaska should result in an improved understanding of high-latitude ecosystems and the evolution of their vertebrate fauna. To resolve these connections between northern and southern Alaska would represent a significant step forward in addressing questions such as how temperature and precipitation trends over the represented range of paleolatitudes impacted the evolution and seasonal or longer-term migration of large-bodied herbivore populations e.g., Ref. [7]. Ultimately, the integration of Cretaceous high-latitude depositional systems, paleo flora and fauna, and paleoclimate in Alaska has the potential to bear on paleogeographic connections between northwestern North America and Eurasia, but much more work is needed.

Conclusions
(1) A thorough geologic and paleogeographic characterization of the Cantwell Formation is important because it provides a stratigraphic framework for recent and future paleontological and paleoclimate investigations of a Late Cretaceous high-latitude ecosystem. Our stratigraphic data allow the Cantwell Formation to be divided into two mappable members. The lower member, the Polychrome member (minimum thickness:~1500 m), consists of wet alluvial fan, braided stream, and coastal fluvial strata. The upper member, the Cabin Creek-Double Mountain member (minimum thickness:~1200 m), was deposited in estuarine and coastal lacustrine depositional environments. The coastal fluvial and tidal environments are associated with the remarkable preservation of vertebrate trackways such as theropods, hadrosaurs, and wading birds.
(2) The configuration of the Cantwell basin was heavily influenced by syndepositional deformation along the northern and southern basin margins. In our study area, the basin margins are defined by opposing thrust belts, including a previously undocumented south-verging thrust system along the northern basin margin. The Cantwell Formation was deformed within these basinward propagating thrust systems from the Late Cretaceous through the Paleocene. Upper Paleocene-lower Eocene volcanic strata of the Teklanika formation were deposited unconformably atop the Cantwell Formation, and contractional deformation of both formations continued through at least the early Eocene. (3) U-Pb detrital zircon data from the Cantwell Formation (n = 4657) indicate that sediment was sourced from a coeval Late Cretaceous magmatic arc, as well as older exhumed Cretaceous plutonic sources. Plutons of the Late Cretaceous arc system that bordered the southern margin of the basin during the deposition of the Cantwell Formation are now exposed~400 km to the southwest in the western Alaska Range due to dextral offset along the Denali fault system. Sediment sources north of the Cantwell basin included Upper Triassic-Lower Jurassic plutonic rocks of the Taylor Mountain batholith, as well as metamorphic rocks of the ancestral North American margin that provided recycled Proterozoic and Archean detrital zircons. Other secondary (recycled sedimentary) sources of sediment include Paleozoic and Mesozoic metasedimentary strata exposed in syndepositional thrust sheets along the margins of the basin and the Upper Jurassic-Cretaceous Kahiltna assemblage. (4) The Cantwell basin is interpreted to have formed in the proximal part of a nowdismembered foreland basin system and drained to the northeast into a Late Cretaceous interior seaway. Future studies may seek to better define the stratigraphic relationships between the interior seaway, likely represented by relatively unstudied Cretaceous strata of interior Alaska, and terrestrial and estuarine depositional systems of south-central and northern Alaska. These regional datasets would serve to better integrate disparate records of the Cretaceous high-latitude paleo-landscape of Alaska and would contribute to our understanding of relationships between depositional systems, faunal populations, and basinal settings. These datasets may also be critical for understanding paleogeographic and faunal connections between northwestern North American and Eurasia.