Marie Salpin

Johann Schnyder

François Baudin

Guillaume C Suan

Jean-Pierre Suc

Speranta-Maria Popescu

Séverine Fauquette

Lutz Reinhardt

Mark D Schmitz

Loic Labrousse

Evidence for subtropical warmth in the Canadian Arctic (Beaufort-Mackenzie, Northwest Territories, Canada) during the early Eocene

établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

INTRODUCTION

The Paleocene to early Eocene has been recognized as the warmest period of the Cenozoic Era. Paleocene-Eocene global warming began in the late Paleocene and continued into the early Eocene (58-50 Ma;Vandenbergh e et al., 2012), reaching an optimum in the early Eocene (52-50 Ma; [START_REF] Zachos | Trends, rhythms, and aberrations in global climate 65 Ma to present[END_REF][START_REF] Vandenberghe | The Paleogene period[END_REF]. This time interval is punctuated by several short-lived (150-200 ky) hyperthermal events [START_REF] Aubry | Causes and consequences of globally warm climates in the early Paleogene[END_REF][START_REF] Röhl | On the duration of the Paleocene-Eocene thermal maximum[END_REF], the greatest in magnitude being the Paleocene-Eocene Thermal Maximum (PETM, 56 Ma) and the Eocene Thermal Maximum 2 (ETM2, 2 My after the PETM) [START_REF] Zachos | Trends, rhythms, and aberrations in global climate 65 Ma to present[END_REF][START_REF] Brinkhuis | Episodic fresh surface waters in the Eocene Arctic Ocean[END_REF][START_REF] Thomas | Descent into the Icehouse[END_REF][START_REF] Stap | High-resolution deep-sea carbon and oxygen isotope records of Eocene Thermal Maximum 2 and H2[END_REF]Storme et al., 2012a;[START_REF] Garel | Paleohydrological and paleoenvironmental changes recorded in terrestrial sediments of the Paleocene-Eocene boundary (Normandy, France)[END_REF]. During these hyperthermal events, the global average temperature rose by as much as 8 °C above pre-event values [START_REF] Zachos | Trends, rhythms, and aberrations in global climate 65 Ma to present[END_REF][START_REF] Jones | Climate model and proxy data constraints on ocean warming across the Paleocene-Eocene Thermal Maximum[END_REF]. These extreme global warming events are systematically associated with negative carbon isotope excursions (CIE) in both marine and terrestrial carbonate (~-3 ‰) as well as in organic carbon (~6 ‰, e.g., [START_REF] Aubry | The Global Standard Stratotype-section and Point (GSSP) for the base of the Eocene Series in the Dababiya section (Egypt): Episodes[END_REF]Storme et al., 2012b), which is thought to refl ect the massive release of isotopically light carbon dioxide and methane into the ocean and atmosphere carbon pools [START_REF] Kennett | Abrupt deep sea warming, paleoceanographic changes and benthic extinctions at the end of the Paleocene[END_REF][START_REF] Pearson | Atmospheric carbon dioxide concentrations over the past 60 million years[END_REF][START_REF] Zachos | Trends, rhythms, and aberrations in global climate 65 Ma to present[END_REF][START_REF] Zachos | An early Cenozoic perspective on greenhouse warming and carbon-cycle dynamics[END_REF][START_REF] Lourens | Astronomical pacing of late Palaeocene to early Eocene global warming events[END_REF]Schouten et al., 2007;[START_REF] Deconto | Past extreme warming events linked to massive carbon release from thawing permafrost[END_REF]. These late Paleocene and early Eocene hyperthermal events induced major turnovers in the assemblage of planktonic and benthic foraminifera, calcareous nannofossils, and dinoflagellates, with an especially pronounced increase in the relative abundance of the dinoflagellate genus Apectodinium during the PETM [START_REF] Kennett | Abrupt deep sea warming, paleoceanographic changes and benthic extinctions at the end of the Paleocene[END_REF][START_REF] Kaiho | Latest Paleocene benthic foraminiferal extinction and environmental changes at Tawanui[END_REF][START_REF] Thomas | Deep-sea environments on a warm earth: Latest Paleocene-early Eocene[END_REF][START_REF] Crouch | The Apectodinium acme and terrestrial discharge during the Paleocene-Eocene thermal maximum: new palynological, geochemical and calcareous nannoplankton observations at Tawanui, New Zealand[END_REF][START_REF] Zachos | Rapid acidifi cation of the ocean during the Paleocene-Eocene thermal maximum[END_REF]. These biotic changes coincided with major changes in ocean and atmosphere circulations and repeated pulses of lysocline and carbonate compensation depth shoaling [START_REF] Zachos | Rapid acidifi cation of the ocean during the Paleocene-Eocene thermal maximum[END_REF][START_REF] Schouten | The Paleocene-Eocene carbon isotope excursion in higher plant organic matter: Differential fractionation of angiosperms and conifers in the Arctic[END_REF]. Extensive paleobotanic studies have been performed in the Paleogene sediments of the Arctic domain. Many of them were based on plant macroremains with identifications most often at the genus level, establishing that many of the modern botanical genera already existed in Paleogene times (e.g.: [START_REF] Basinger | Early Tertiary vegetation of Arctic Canada and its relevance to paleoclimatic interpretation[END_REF][START_REF] Greenwood | The paleoecology of high-latitude Eocene swamp forests from Axel Heiberg Island[END_REF][START_REF] Jahren | The Arctic forest of the middle Eocene[END_REF][START_REF] Mciver | Early Tertiary floral evolution in the Canadian high Arctic[END_REF][START_REF] Sunderlin | Paleoclimatic and paleoecological implications of the Paleocene-Eocene fossil leaf assemblage[END_REF][START_REF] Eberle | Life at the top of the greenhouse Eocene world-A review of the Eocene flora and vertebrate fauna from Canada's High Arctic[END_REF][START_REF] Eldrett | A seasonality trigger for carbon injection at the Paleocene-Eocene Thermal Maximum[END_REF][START_REF] West | Was the Arctic Eocene 'rainforest' monsoonal? Estimates of seasonal precipitation from early Eocene megafl oras from Ellesmere Island, Nunavut[END_REF][START_REF] Harrington | Arctic plant diversity in the Early Eocene greenhouse[END_REF]and older references). Simultaneously, studies based on pollen grain identification under the transmitted light microscope developed an identification approach mixing the artificial and botanic nomenclatures [START_REF] Ridgway | Climatically induced floristic changes across the Eocene-Oligocene transition in the northern high latitudes[END_REF][START_REF] Kalkreuth | Petrological, palynological and geochemical characteristics of Eureka Sound Group coals (Stenkul Fiord, southern Ellesmere Island, Arctic Canada[END_REF][START_REF] Frederiksen | Palynology of Eocene strata in the Sagavanirktok and Canning formations of the north slope of Alaska[END_REF][START_REF] White | Palynology, age, correlation and paleoclimatology from JAPEX/JNOC/GSC Mallik 2L-38 gas hydrate research well and the significance for gas hydrates: a new approach: Geological Survey of Canada[END_REF][START_REF] Sluijs | Warm and wet conditions in the Arctic region during Eocene Thermal Maximum 2[END_REF]. The result was a taxonomic gap between macroremain and pollen floras. Palynologists have put forth serious efforts in botanical identification aimed at reducing and closing this gap, first using only the transmitted light microscope [START_REF] Eldrett | Increased seasonality through the Eocene to Oligocene transition in northern high latitudes[END_REF][START_REF] Barke | Orbitally forced Azolla blooms and Middle Eocene Arctic hydrology: clues from palynology[END_REF][START_REF] Harrington | Arctic plant diversity in the Early Eocene greenhouse[END_REF][START_REF] Zaporozhets | The Middle and Upper Eocene sections of the Omsk Trough, West Siberian Platform: palynological, stratigraphic, hydrologic and climatic aspects[END_REF], then with the scanning electronic microscope [START_REF] Liu | Fossil Cathaya (Pinaceae) pollen from the Canadian High Arctic[END_REF][START_REF] Grimsson | Middle Eocene palynoflora from a resinite-rich coal bed on Hareøen (Qeqetarsuatsiaq)[END_REF][START_REF] Suan | Subtropical climate conditions and mangrove growth in Arctic Siberia during the early Eocene[END_REF]. The IODP-ACEX wells (2004) were cored on the Lomonosov Ridge in the central Arctic Ocean in order to increase information thanks to a long and continuous record [START_REF] Moran | The Cenozoic palaeo environment of the Arctic Ocean[END_REF][START_REF] Sluijs | Arctic late Paleocene-early Eocene paleoenvironments with special emphasis on the Paleocene-Eocene thermal maximum[END_REF][START_REF] Backman | Expanding the Cenozoic paleoceanographic record in the Central Arctic Ocean: IODP Expedition 302 Synthesis[END_REF]. The results from the ACEX expedition indicate euxinic conditions and high primary productivity in the central Arctic Ocean during the Paleogene and early Eocene [START_REF] Stein | Anoxia and high primary production in the Paleogene central Arctic Ocean: First detailed records from Lomonosov Ridge[END_REF]. Tetraether lipid based proxies (TEX86 and MBT) revealed an increase of air and sea surface temperatures (SST) of up to 8 °C during the PETM [START_REF] Sluijs | Subtropical Arctic Ocean temperatures during the Palaeocene/Eocene thermal maximum[END_REF][START_REF] Sluijs | Warm and wet conditions in the Arctic region during Eocene Thermal Maximum 2[END_REF][START_REF] Weijers | Warm arctic continents during the Palaeocene-Eocene thermal maximum[END_REF] and a subsequent 3-5 °C SST warming during the ETM2 [START_REF] Sluijs | Warm and wet conditions in the Arctic region during Eocene Thermal Maximum 2[END_REF]. On the Siberian margin, floral assemblages indicate a paratropical humid climate [START_REF] Akhmetiev | Paleogene floral assemblages around epicontinental seas and straits in Northern Central Eurasia: proxies for climatic and paleogeographic evolution[END_REF][START_REF] Suan | Subtropical climate conditions and mangrove growth in Arctic Siberia during the early Eocene[END_REF]. Floral and faunal assemblages indicative of exceptionally warmer and more humid conditions during the late Paleocene-early Eocene than today have been also reported in other localities around the Arctic Ocean (Ellesmere Island, e.g., [START_REF] Eberle | Seasonal variability in Arctic temperatures during early Eocene time[END_REF][START_REF] Schubert | A summertime rainy season in the Arctic forests of the Eocene[END_REF]Northeast Yukon, Vavrek et al., 2012;Axel Heiberg Island, e.g., McIver and Basinger, 1999;Northwest Territories (NWT), Ellesmere and Axel Heiberg Islands: [START_REF] Greenwood | How wet was the Arctic Eocene rain forest? Estimates of precipitation from Paleogene Arctic macrofloras[END_REF]Canadian Arctic: Eberle and Greenwood, 2012;Alaska: Sunderlin et al., 2011). Although these studies have provided valuable constraints for climate conditions in the Arctic during the Paleocene-Eocene, they are limited to a few localities (Fig. 1). Therefore, additional studies are required to ascertain the effects of Paleocene-Eocene global environmental changes on western Arctic regions. This study aims to constrain paleoclimatic conditions of the Canadian margin during the Paleocene-Eocene. We focus on the Mackenzie Delta of the Beaufort-Mackenzie Basin in the Arctic Basin with a multiproxy approach comprising organic matter characterization, clay mineralogy, comprehensive botanical palynology, and carbon isotope geochemistry on dispersed organic matter.

GEOLOGICAL SETTING AND SAMPLES

The Mackenzie Delta lies north of the Arctic Circle (68-69°N), on the coastline of the Northwest Territories, Canada (Fig. 1). This fluvial-dominated delta is 200 km long, and 80 km wide, oriented along a NNW axis [START_REF] Dixon | Cretaceous to Pleistocene stratigraphy and paleogeography, northern Yukon and northwestern District of Mackenzie[END_REF], and was located at a paleolatitutde of ~75°N during the Eocene [START_REF] Van Hinsbergen | A paleolatitude calculator for paleoclimate studies[END_REF]. A high sedimentation rate has been recorded in the Beaufort-Mackenzie Basin, with Cretaceous to Cenozoic deposits being 12-16 km thick [START_REF] Dixon | A review of Cretaceous and Tertiary stratigraphy in the northern Yukon and adjacent Northwest Territories: Ottawa, Canada[END_REF]. The Cenozoic strata have been penetrated by many petroleum wells but are rarely exposed in outcrops. The most extensive outcrops are found along the Caribou Hills on the eastern part of the Mackenzie Delta (Fig. 1). Based on outcrop and well data, twelve transgressive-regressive sedimentary sequences have been recognized within Upper Cretaceous to Quaternary sedimentary succession in the Beaufort-Mackenzie Basin, each sequence being bounded by major unconformities [START_REF] Dixon | Geology of the Late Cretaceous to Cenozoic Beaufort-Mackenzie Basin, Canada[END_REF]. The basin infill seems partly controlled by the local tectonic activity, at least until the late Miocene. The major tectonic pulses of deformation in the region spanned from Paleocene to Miocene during the development of the Canadian Cordillera ending to the North in the Richardson Mountains (Fig. 1) and the Brooks Range-British Mountains belt along the Arctic Ocean ending offshore with the Beaufort Fold Belt [START_REF] Harrison | Correlation of Cenozoic sequences of the Canadian Arctic region and Greenland; implications for the tectonic history of northern North America[END_REF][START_REF] Dixon | Geology of the Late Cretaceous to Cenozoic Beaufort-Mackenzie Basin, Canada[END_REF]. The dating of these sedimentary sequences is based on the identification of foraminiferal and palynological assemblages (spore, pollen grains, and dinoflagellate cysts; [START_REF] Mcneil | Foraminiferal zonation and biofacies analysis of Cenozoic strata in the Beaufort-Mackenzie Basin of Arctic Canada[END_REF][START_REF] Dixon | A review of Cretaceous and Tertiary stratigraphy in the northern Yukon and adjacent Northwest Territories: Ottawa, Canada[END_REF]. The Paleogene strata constitute a series of predominantly delta plain deposits made up of unconsolidated sandy to gravel sedimentary packages corresponding to fluvial channel belts, as well as groups of shaly horizons indicative of levee or deltaic plain environment, and occasional thin coal beds of wetland environment [START_REF] Price | Revision of the Tertiary Reindeer Formation in the Caribou Hills, District of Mackenzie[END_REF]. The occurrence of the benthic foraminifera Reticulophragmium sp. at the base of the sequence and the associations of pollen grains, dinoflagellates, and foraminifera above led to the assignment of 720 m out of the total 1200 m of Cenozoic strata exposed at Caribou Hills to the late Paleocene-middle Eocene Aklak and Taglu depositional sequences [START_REF] Parsons | Palynology of Paleogene strata in the Caribou Hills[END_REF]Fig. 2). The rest of the Cenozoic succession is thought to be Oligocene to Pliocene in age (Kugmalitt to Iperk sequences; [START_REF] Parsons | Palynology of Paleogene strata in the Caribou Hills[END_REF]. The boundary between Aklak and Taglu sequences is marked by an angular unconformity, which is well documented in seismic data but difficult to observe on outcrops [START_REF] Dixon | A review of Cretaceous and Tertiary stratigraphy in the northern Yukon and adjacent Northwest Territories: Ottawa, Canada[END_REF]. For this study, two sections have been sampled along the Caribou Hills exposure during the CASE 15 (Circum-Arctic Structural Events) expedition in 2013 by L. Labrousse (ISTeP, France) and L. Reinhardt (BGR, Germany). These two sections are 26 m thick (CaH1) and 50 m thick (CaH2) and belong to the Taglu and/or Aklak sequence. CaH2 is estimated to be stratigraphically ca. 400 m above CaH1. The CaH1 section is mainly composed of silty clays more or less rich in organic matter, which are interrupted by three sandy levels at the base of the section. A bentonite layer occurs at 23 m, which is characterized by its prominent orange color when weathered. It consists mainly of smectite with some traces of goethite and gypsum. The fresh unweathered bentonite of the ~10 cm thick layer has a light beige color. Above this layer, an additional set of ~4 cm thick, orange bentonite layers are preserved in a sequence of brownish siltstones. The horizontally bedded sediments lack signs of reworking. The CaH2 section is composed of a succession of coarse-grained sand, fi ne-sand, and silty clay beds, the latter being sometimes rich in organic matter. Four lignite seams are interbedded in the fi nest deposits. Assuming that the bentonite level measured at 23 m in the CaH1 section is the same as that identified by [START_REF] Price | Revision of the Tertiary Reindeer Formation in the Caribou Hills, District of Mackenzie[END_REF]Fig. 3), CaH1 and CaH2 are respectively equivalent to section 2 and 6 of [START_REF] Parsons | Palynology of Paleogene strata in the Caribou Hills[END_REF]. Moreover, bentonite levels were found in Ellesmere Island and were assumed to result from explosive volcanism during the opening of the Labrador Sea, Baffin Bay, and northern Atlantic [START_REF] Reinhardt | Altered volcanic ashes in Palaeocene and Eocene sediments of the Eureka Sound Group (Ellesmere Island, Nunavut, Arctic Canada)[END_REF]. Facies correspond to migrating channel belts, laterally associated with levee, crevasse splay, and swamp deposits, characteristic of a delta plain environment. A total of 35 and 22 samples were collected along the CaH1 and CaH2 sections, respectively.

METHODS

Clay Mineralogy

Clay mineral (<2μm) associations were identified on all samples by X-ray diffraction (XRD) using a Bruker D2 phaser X-ray powder diffractometer equipped with a LynxEye detector (CuKradiation, = 1.54Å and Ni filter). Oriented mounts were prepared for clay minerals analyses and three X-ray diagrams per sample were performed following the analytical procedure described in [START_REF] Holtzapffel | Les minéraux argileux: préparation, analyse diffractométrique et détermination[END_REF] and [START_REF] Moore | X-ray Diffraction and the Identification and Analysis of Clay Minerals[END_REF]. Semiquantitative proportions of the identified minerals were obtained using the peak area of the basal reflection for each mineral, using MacDiff® 4.2.5 software. Some rock pieces were observed under Scanning Electron Microscopy to determine crystal morphology.

U-Pb Isotopic Dating (ID-TIMS)

Sample CASE15_YU067 (coordinates N 68° 3819.5/ W 134° 0337.4, handheld GPS) from the 10 cm thick bentonite layer in section CaH1 was processed for U-Pb zircon geochronology and measured (ID-TIMS) at the Boise State University Isotope Geology Laboratory, Idaho (USA). The results are shown in Table 1 andFigure 3. Heavy minerals were separated from 0.5 kg aliquots of volcanic tuff via progressive washing under ultrasonication following the method of [START_REF] Hoke | An ultrasonic method for isolating nonclay components from clay-rich material[END_REF], followed by conventional density and magnetic methods. Separate from each sample, the entire heavy mineral was placed in a muffle furnace at 900 °C for 60 h in quartz beakers to anneal minor radiation damage in zircon. Annealing enhances cathodoluminescence (CL) emission [START_REF] Nasdala | Annealing radiation damage and the recovery of cathodoluminescence[END_REF], promotes more reproducible interelement fractionation during laser ablation-inductively coupled plasma mass spectrometry (LA-ICPMS) [START_REF] Allen | Identification and elimination of a matrix-induced systematic error in LA-ICP-MS 206 Pb/ 238 U dating of zircon[END_REF], and prepares the crystals for subsequent chemical abrasion [START_REF] Mattinson | Zircon U-Pb chemical abrasion ("CA-TIMS") method: combined annealing and multi-step partial dissolution analysis for improved precision and accuracy of zircon ages[END_REF]. Following annealing, individual grains were handpicked for further isotopic analysis. U-Pb geochronology methods for isotope dilution thermal ionization mass spectrometry followed those previously published by [START_REF] Davydov | Highprecision U-Pb zircon age calibration of the global Carboniferous time scale and Milankovitch-band cyclicity in the Donets Basin, eastern Ukraine[END_REF] and [START_REF] Schmitz | Quantitative radiometric and biostratigraphic calibration of the global Pennsylvanian-Early Permian time scale[END_REF]. Zircon crystals were subjected to a modified version of the chemical abrasion method of [START_REF] Mattinson | Zircon U-Pb chemical abrasion ("CA-TIMS") method: combined annealing and multi-step partial dissolution analysis for improved precision and accuracy of zircon ages[END_REF], whereby single crystal fragments plucked from grain mounts were individually abraded in a single step with concentrated HF at 180 °C for 12 h. All analyses were undertaken on crystals previously mounted, polished, and imaged by cathodoluminence (CL), and selected on the basis of zoning patterns. U-Pb dates and uncertainties for each analysis were calculated using the algorithms of [START_REF] Schmitz | Derivation of isotope ratios, errors and error correlations for U-Pb geochronology using 205 Pb-235 U-( 233 U)-spiked isotope dilution thermal ionization mass spectrometric data[END_REF] and the U decay constants of [START_REF] Jaffey | Precision measurements of half-lives and specifi c activities of 235 U and 238[END_REF]. Uncertainties are based upon non-systematic analytical errors, including counting statistics, instrumental fractionation, tracer subtraction, and blank subtraction. These error estimates should be considered when comparing our 206 Pb/ 238 U dates with those from other laboratories that used tracer solutions calibrated against the EARTHTIME gravimetric standards. When comparing our dates with those derived from other decay schemes (e.g., 40 Ar/ 39 Ar, 187 Re-187 Os), the uncertainties in tracer calibration (0.03%; [START_REF] Condon | Metrology and traceability of U-Pb isotope dilution geochronology (EARTHTIME Tracer Calibration Part I)[END_REF][START_REF] Mclean | Evaluating uncertainties in the calibration of isotopic reference materials and multi-element isotopic tracers (EARTHTIME Tracer Calibration Part II)[END_REF] and U decay constants (0.108%; [START_REF] Jaffey | Precision measurements of half-lives and specifi c activities of 235 U and 238[END_REF]) should be added to the internal error in quadrature. Quoted errors for calculated weighted means are thus of the form X(Y) [Z], where X is solely analytical uncertainty, Y is the combined analytical and tracer uncertainty, and Z is the combined analytical, tracer, and 238 U decay constant uncertainty.

Rock-Eval Pyrolysis

Rock-Eval pyrolysis was performed on all samples using a Rock-Eval 6 under standard conditions [START_REF] Espitalié | La pyrolyse Rock-Eval et ses applications[END_REF][START_REF] Behar | Rock-Eval 6 technology: performances and developments[END_REF]. The following parameters were determined: the Total Organic Carbon (TOC) content in weight %, reflecting the quantity of pyrolysed carbon and residual organic carbon; the Hydrogen Index (HI) expressed in mg of hydrocarbons by gram of TOC; the Oxygen Index (OI) expressed in mg of carbon dioxide by gram of TOC, constraining kerogen types and their preservation state; and the Tmax expressed in °C, used to estimate the degree of thermogenic maturation of the kerogen (Fig. 4).

Palynofacies

A total of 22 samples were selected for palynofacies observations, i.e., the identification of all particulate organic matter under the microscope in transmitted light [START_REF] Tyson | Sedimentary Organic Matter: Organic Facies and Palynofacies[END_REF]. Samples were crushed into millimetric pieces using a ceramic mortar, oven-dried, and treated with HCl and HF to remove the mineral matrix, and both total and sieved (fraction below 10 micron being removed) slides were prepared with the residues. The observations were carried out with an Axioplan2 Imaging Zeiss optical microscope under transmitted light and under UV excitation (Zeiss HBO 100 Microscope Illuminating System, mercury short-arc lamp). Over 1000 basic surface units per sample were counted to estimate the relative proportions of each organic group. The organic particles were grouped into the following categories: amorphous organic matter (AOM), phytoclast, and palynomorph groups. Based on shape, color, and fluorescence, the AOM can be subdivided into diffuse, granular, and gelified. The diffuse form is mainly yellow in natural light and exhibits fluffy shapes and weak to no fluorescence. The granular form is mainly gray to brown and exhibits irregular aggregated shapes and weak to no fluorescence. It sometimes contains small framboidal pyrite. The diffuse and granular AOM is thought to derive from phytoplankton, algae, and bacteria that accumulated in O2-depleted marine or lacustrine waters [START_REF] Tyson | Sedimentary Organic Matter: Organic Facies and Palynofacies[END_REF], but may also derive from degraded terrestrial organic matter [START_REF] Batten | Palynofacies and paleoenvironmental interpretation[END_REF]. The gelified AOM is orange to brown and lacks internal structures. As gelified AOM is here likely to correspond to degraded gelified phytoclasts [START_REF] Tyson | Sedimentary Organic Matter: Organic Facies and Palynofacies[END_REF], we will further refer to gelified phytoclasts plus gelified AOM as gelified OM. Phytoclasts are plant-derived fragments, including cuticles, woody tissues, and charcoals. In this study the phytoclast group is separated into opaque particles, translucent ligno-cellulosic fragments, and gelified particles, some very dark in color that look like gelified semi-opaque particles.

Palynology

A total of 35 samples were analyzed specifically for palynology. Pollen slide preparation followed a standard gentle method resulting in no damage to palynomorph membranes [START_REF] Cour | Nouvelles techniques de détection des fl ux et retombées polliniques: Etude de la sédimentation des pollens et des spores à la surface des sols[END_REF]. Samples (15-25 g dry weight) were treated with HCl (35%) and HF (70%) to remove carbonates and silicates. Denser particles were separated from the organic residue using ZnCl2 (density = 2.0). Residues were filtered through a 150 μm nylon sieve to eliminate the coarser particles, including organic macroremains. Palynomorphs were further concentrated using a 10 μm nylon sieve. The final residue was then homogenized and mounted onto microscope slides with glycerol. A transmitted light microscope was used for counting (at magnification × 200) and identification (at magnification × 1000) of palynomorphs. Pollen grains are particularly well preserved, a condition that allows a detailed examination of their morphological characters and their identification at the botanical family or genus level supported by comparison of the fossil specimens with those from databases of present-day species (collection of pollen slides, photographs, and iconographic atlases). Information on some critical pollen identifications is provided in the GSA Data Repository Item1. A minimum of 150 pollen grains were counted per sample in addition to those of Pinus. Detailed data are provided in Tables 2 and3. Results are summarized in synthetic pollen diagrams presented in Figure 5B. Botanical identification of pollen grains allows the grouping of the recorded elements with respect to [START_REF] Nix | Environmental determinants of biogeography and evolution in Terra Australis[END_REF] classification of living plants as reported in Tables 2 and3: megatherm (tropical) elements (mean annual temperature (MAT)>24 °C); mega-mesotherm (subtropical) elements (20 °C<MAT<24 °C); Cathaya, a conifer living today in tropical-subtropical regions at mid-to high-elevations but below the Abies-Picea belt; mesotherm (warm-temperate) elements (14 °C<MAT<20 °C); mesomicrotherm (cool-temperate) elements (12 °C<MAT<14 °C); and microtherm (boreal) elements (MAT<12 °C). The other elements are composed of: elements without climatic signification apart from some Cupressaceae (Juniperus-Cupressus-type), which occupy various biotopes (the pollen of which cannot be identified at the genus level); hygrophilous elements (i.e., water plants); and various herbs (Fig. 5B; Tables 2 and3). A pollen ratio is calculated between "thermophilous" (megatherms + megamesotherms + mesotherms) and "non-thermophilous" (Cathaya and Podocarpus-type + meso-microtherms + microtherms) taxa, which is regarded as a representative index of temperature variation (Tables 2 and3). Values of this ratio significantly >1 characterize warm periods [START_REF] Popescu | Repetitive changes in Early Pliocene vegetation revealed by high-resolution pollen analysis: revised cyclostratigraphy of southwestern Romania[END_REF][START_REF] Popescu | Pollen record and integrated high-resolution chronology of the early Pliocene Dacic Basin (southwestern Romania)[END_REF].

Carbon Isotope Analyses

The organic carbon-isotope composition (13Corg) was determined for 42 samples that were powdered using a grinder and decarbonated using 1N HCl for 2 h at 25 °C prior to analysis. Depending on the TOC contents, between 0.2 and 10 mg of fully decarbonated sediment powder was weighted into tin capsules and placed in a PYRO cube® elemental analyzer connected to an Elementar Isoprime® 100 isotope-ratio mass spectrometer in continuous flow. Each analytical run contained four sets of two standards (IAEA CH7, δ 13 C = -32.15 ‰; Caseine, δ 13 C = -27.18 ‰) to monitor analytical precision and accuracy. The carbon isotope results are reported relative to the "Vienna PeeDee belemnite" (VPDB) in delta notation. The precision, based on replicate analysis of standard and samples, was better than 0.13 ‰.

Paleoclimatic Reconstruction Based on Palynology

The climate has been reconstructed, thanks to palynology, for the Caribou Hills deposits using the "Climatic Amplitude Method" in which the climatic requirements of the maximum number of modern taxa are transposed to the fossil data (Fauquette et al., 1998a(Fauquette et al., , 1998b)). This approach relies on the relationship between the relative pollen abundance of each individual taxon and climate. It takes into account not only presence/absence criterion but also pollen percentages to provide more reliable reconstructions. The estimated mean annual temperatures (MAT), coldest month mean temperature (CMMT), warmest month mean temperature (WMMT), mean annual precipitation (MAP), and available moisture concern the low-altitude vegetation because meso-microthermic and microthermic taxa have been excluded from the process to avoid a cold bias linked to transport from higher elevations. The excluded taxa are defined on the basis of their occurrence in modern vegetation zones [START_REF] Walter | Vegetation of the Earth[END_REF], as vegetation types described for studied time period are found today in southeastern China from around 25-30°N and 110-120°E [START_REF] Hou | Vegetation of China with reference to its geographical distribution[END_REF]. In this region of China, the vertical distribution of the vegetation is characterized from the base to the top of the massifs (Taibai Shan in the Tsinling Massif, Shaanxi, Taba Shan, eastern Sichuan, or Hua Ping, Guangxi), by the evergreen broad-leaved forest, mixed evergreen and deciduous (Betula, Acer) broad-leaved forest, Cathaya/Tsuga forest, Picea/Abies forest, and high mountain meadows [START_REF] Wang | The Forests of China with a Survey of Grassland and Desert Vegetation[END_REF][START_REF] Hou | Vegetation of China with reference to its geographical distribution[END_REF]. This modern distribution allows us to exclude the modern highaltitude taxa from the climatic reconstruction to obtain the low to middle-low altitude climate. The most probable climate for a fossil pollen assemblage is estimated as the climatic interval in which the highest number of taxa can exist and a "Most Likely Value" (MLV) corresponding to a weighted mean, according to the size of the climatic intervals of all taxa. As the precision of the information obtained from a taxon's climatic interval is inversely related to the breadth of this interval, the weights are greater for taxa with smaller intervals (for more details, see Fauquette et al., 1998a).

RESULTS

Clay Mineralogy

The clay fraction (<2 μm) mineralogy of samples of CaH1 and CaH2 sections is composed of kaolinite (0-75%, average value 30%), illite (0-70%, average value 30%), smectite (0-100%, average value 30%), and chlorite (0-30%, average value 15%) (Fig. 5A). Important differences in clay mineralogy are obvious between the two sections. Section CaH1 is characterized by high proportions of smectite (10-100%, mean value 50%) when compared to CaH2, which is dominated by kaolinite (20-75%, average value 45%) (Fig. 5A). Illite is the second clay mineral represented in both sections and may correspond to the background physical erosion from the hinterland (Chamley, 1989). The illite crystallinity does not show continuous improvement along the sections, so no diagenetic overprint is highlighted. Chlorite is more abundant in CaH2 than in CaH1, a feature that may reflect local erosional inputs, as chloritic minerals are often thought to result from the physical erosion of catchment areas (Chamley, 1989). Kaolinite is generally poorly crystallized according to Scanning Electron Microscopy (SEM) observations. Only small to moderate amounts of well-crystallized crystals have been observed in some pure sand samples, showing that kaolinite has mainly detrital origins, and only localized, small amounts of kaolinite may be of diagenetic origin. Kaolinite forms under intense hydrolyzing conditions and is considered to reflect warm and wet climatic conditions (Chamley, 1989;[START_REF] Velde | The Origin of Clay Minerals in Soils and Weathered Rocks[END_REF][START_REF] Tramoy | Stepwise palaeoclimate change across the Eocene-Oligocene transition recorded in continental NW Europe by mineralogical assemblages and δ 15 Norg[END_REF]. Smectite forms under colder and/or dryer conditions, in poorly drained environments or under contrasted seasonality (Chamley, 1989;[START_REF] Tramoy | Stepwise palaeoclimate change across the Eocene-Oligocene transition recorded in continental NW Europe by mineralogical assemblages and δ 15 Norg[END_REF]. The (Smectite + Kaolinite) / (Illite + Chlorite) ratio (S+K / I+C) may be used as a relative weathering indicator as smectite and kaolinite are secondary form clay minerals, whereas illite and chlorite are derived from physical erosion of parent rocks. The higher the ratio, the more elevated were the inland hydrolyzing conditions compared to physical erosion. The Smectite/Kaolinite ratio (S/K) may be used as a relative paleoclimatic proxy, with a high S/K ratio indicating dryer and/or colder climate conditions. The (S+K / I+C) ratio is around 1 in most of the CaH1 section and decreases toward the top of the section. It varies from 1.5 to 3 in CaH2, and reaches high values between 20 and 30 m and around 40-45 m (Fig. 5A). The S/K ratio is high in the CaH1 section, mostly above 1 with three really high values (up to 4.9), whereas it is below 0.5 in most of the CaH2 section. In CaH1, the kaolinite content is higher between 19-21 m and 24-26 m, when compared to the rest of the section. Conversely, smectite proportions are higher at the base and between 21 and 24 m when compared to the rest of the section. The bentonite level at the top of the section contains mostly smectite and traces of illite and kaolinite. In CaH2, the kaolinite content is higher than in CaH1 and decreases from the base of the section up to 21 m, where the illite proportion is higher than the kaolinite proportion. From 21 m to 30 m, the proportion of kaolinite increases and then decreases toward the top of the section (Fig. 5A).

U-Pb Isotopic Dating (ID-TIMS)

The heavy mineral separation of sample CASE15_YU067 resulted in an abundant population of sharply facetted, prismatic zircon crystals, with accompanying abundant prismatic apatite. Five grains were selected for CA-TIMS analysis on the basis of external morphology, clarity, and consistency. Chemical abrasion in concentrated HF at 190° for 12 h resulted in only slight dissolution of the zircon crystals. Moderately high uranium concentration of these crystals resulted in modest radiogenic Pb contents and precisely measured isotope ratios. All five crystals yielded concordant and equivalent isotope ratios, with a weighted mean 206 Pb/ 238 U date of 54.304 ± 0.029(0.039)[0.070] Ma (MSWD = 1.18; probability of fit = 0.317), which is interpreted as dating the eruption and deposition of this tuff.

Rock-Eval Pyrolysis

The Rock Eval pyrolysis analyses indicate that TOC values vary from 0 to 56% and the HI values are below 350 mgHC/gTOC. These values reveal a Type III (terrestrial derived) organic matter for most of the samples with minor marine influence for some samples (Fig. 5A). Tmax values are mainly below 435 °C (359-438 °C, with an average of 411 °C), indicative of immature kerogen. This Tmax variation range is not uncommon for outcrop samples and cannot be assigned to differences in maturity of the organic matter but is mostly due to weathering effect on the organic matter. The TOC values from the CaH1 section (0-8.65%, with an average of 2.36%) are lower than TOC values from CaH2 (0-51%, with an average of 12%), reflecting the occurrence of coal seams in CaH2. On average the CaH2 samples have higher HI values (26-322 mgHC/gTOC, with an average of 192 mgHC/gTOC) compared to that of CaH1 (23-158 mgHC/gTOC, with an average of 60 mgHC/gTOC) (Fig. 5B), suggesting that organic matter is probably better preserved in CaH2 than in CaH1. In CaH2, TOC values are higher between 25 and 50 m when compared to the rest of the section (Fig. 5A) due to the occurrence of lignite beds.

Palynofacies

Palynomorphs are observed in minor proportions and are dominated by bisaccate pollen grains. Spores are also present. Very few palynomorphs belonging to freshwater algae (Pediastrum) and dinoflagellate cysts have been observed. The development of Pediastrum algae is generally associated with an eutrophic freshwater environment [START_REF] Tyson | Sedimentary Organic Matter: Organic Facies and Palynofacies[END_REF][START_REF] Garel | Paleohydrological and paleoenvironmental changes recorded in terrestrial sediments of the Paleocene-Eocene boundary (Normandy, France)[END_REF]. The samples from CaH1 are dominated by gelified organic matter, palynomorphs (i.e., pollen grains and/or spores and/or dinoflagellate cysts), and opaque phytoclasts, confirming a dominant terrestrial origin of the organic matter. Diffuse AOM also occurs in minor proportions. CaH2 samples also contain a large amount of gelified organic matter, which, together with phytoclasts, opaque phytoclasts, and palynomorphs, again indicates a terrestrial origin of the organic matter. Moreover, lower proportions of opaque phytoclasts and palynomorphs are recorded when compared to the CaH1 section. Finally, CaH2 samples contain a larger amount of translucent phytoclasts (Fig. 5B). At the base of the CaH1 section, palynofacies are dominated by granular AOM associated with sandy deposits and also characterized by lower TOC values. In contrast, the top of the section is rich in gelified OM correlated to higher TOC values. It also contains important palynomorph proportions. At the base of CaH2, granular AOM dominates the palynofacies in sandy deposits with low TOC. The rest of the section is dominated by gelified OM, phytoclasts associated with higher TOC (Fig. 5B) corresponding to well-preserved OM from lignite formed in swamps.

Pollen Grains and Dinoflagellate Cysts

The pollen flora of CaH1 includes 69 taxa that have been grouped in a synthetic diagram (Fig. 5B) according to the mean annual temperature (MAT) range of their present-day relatives [START_REF] Nix | Environmental determinants of biogeography and evolution in Terra Australis[END_REF]. Mega-mesotherm plants and Cupressaceae are dominant. Mesotherm plants, mesomicrotherm plants, and microtherm plants were also recorded but in comparatively lower proportions. Several specimens of the dinoflagellate cyst Apectodinium spp. have been recorded at the top of the section (24.8 m). The proportion of Pinus is very low along the section and decreases to the top with the exception of one sample at 19.5 m. The "thermophilous/nonthermophilous taxa" ratio, continuously >1, increases toward the top of the section (Table 2). The pollen flora of CaH2 yielded 102 taxa and is also dominated by mega-mesotherm plants.

In contrast with CaH1, megatherm plants including Avicennia (Fig. 6A-6D; Data Repository Plate DR4) and mesotherm plants are abundant (Fig. 5B). Pollen grains of microtherm plants are rare but with larger percentage than in CaH1. At the top of the section (near 41 and 50 m) a large number of individual Azolla microspores, a freshwater fern, have been recorded. The proportion of Pinus is higher than in CaH1 and the thermophilous/non-thermophilous taxa ratio, still >1, shows very high values in the uppermost part of the section (Table 3).

Carbon Isotopes

The δ 13 Corg values measured in CaH1 and CaH2 range from -27.7 to -23.4‰ (Fig. 5B). The average δ 13 Corg values for CaH1 are -26.1‰ and -25.6‰ for CaH2. These values are similar to those of modern C3 plants and to other δ 13 Corg values recorded for Paleocene-Eocene terrestrial organic matter around the world [START_REF] Magioncalda | Paleocene-Eocene carbon isotope excursion in organic carbon and pedogenic carbonate: Direct comparison in a continental stratigraphic section[END_REF][START_REF] Jaramillo | Effects of rapid global warming at the Paleocene-Eocene boundary on neotropical vegetation[END_REF][START_REF] Garel | Paleohydrological and paleoenvironmental changes recorded in terrestrial sediments of the Paleocene-Eocene boundary (Normandy, France)[END_REF]. In CaH1 the values fluctuate around -26‰, with one sample at 20.95 m recording a value of -27.7‰ and two samples at 23 and 24 m showing slightly higher values of -25‰. In CaH2, δ 13 Corg values show an increase from the base to 30 m, interrupted by three excursions to lower values (< -26‰) at 0.2, 3.2, and 40.35 m. One sample near 8 m has a distinctly higher value of -23.4‰.

Palaeoclimate Quantification Based on Palynology

The climatic quantification based on the pollen flora of the CaH1 section shows high temperatures (MAT, CMMT, WMMT) and high precipitation at sea level for this part of the globe. Mean annual temperatures are between 16 and 25 °C, with "Most Likely Values" (MLV) between 18.2 and 19.3 °C, CMMT between 5 and 15.5 °C with MLV between 5.6 and 12 °C, and WMMT between 25 and 28 °C with MLV between 25.8 and 26.4 °C along the section. Mean annual precipitation and available moisture are continuously high, at 1100-1600 mm (with MLV of 1160-1400 mm) and 66-100% (with MLV of 77-87%) (Fig. 5B), respectively. The quantification based on the pollen flora of the CaH2 section shows, from the base of the section to 41.2 m, higher temperatures compared to CaH1, especially for CMMT (Fig. 5B). Mean annual precipitation estimates are similar in the two sections (Fig. 5B). From the base to 41.2 m, MAT are between 16 and 25 °C with MLV between 21.3 and 22 °C (except 2 levels at 18.5 and 19 °C; Fig. 5B), CMMT between 10 and 15.5 °C with MLV of 10.5-14 °C, and WMMT between 25 and 30 °C with MLV of 26.2-27.7 °C along the section. MLV of MAP and available moisture are between 1200 mm and 1360 mm and between 75% and 92%, respectively. From 41.5 m to the top, with the exception of 2 points showing higher temperatures, MAT are between 15 and 25 °C with MLV of 17.8-19.7 °C, CMMT between 5 and 15.5 °C with MLV between 7.6 and 11.6 °C, and WMMT between 25 and 28 °C with MLV between 25.8 and 26.3 °C. Most Likely Values of MAP and available moisture are between 1160 and 1400 mm and between 81% and 92%, respectively.

DISCUSSION

Age of the Sections

Caribou Hills 1 (CaH1) [START_REF] Parsons | Palynology of Paleogene strata in the Caribou Hills[END_REF], in his section 2, equivalent to CaH1, reported a foraminifer taxon (Portatrocharnmina tagluensis), which is characteristic of the Taglu sequence and is indicative of an early to middle Eocene age on the basis of associated dinoflagellate cysts and a rich assemblage of pollen grains and spores as defined by McNeil (1997) on his regional biostratigraphic chart. Thus, based on equivalence with section 2 in [START_REF] Parsons | Palynology of Paleogene strata in the Caribou Hills[END_REF], CaH1 probably corresponds to the lowest part of the Taglu sequence or the upper part of the Aklak sequence and is here considered as early Eocene in age. In our samples, Aquilapollenites pollen grains have been found in the CaH1 succession and could also indicate an early Eocene age [START_REF] Harrison | Correlation of Cenozoic sequences of the Canadian Arctic region and Greenland; implications for the tectonic history of northern North America[END_REF]. This supposition is supported by the bentonite level U-Pb age of 54.304 Ma for the top of the CaH1 section (Figs. 3, 5A, and 5B; Table 1). According to [START_REF] Westerhold | Astronomical calibration of the Ypresian timescale: implications for seafloor spreading rates and the chaotic behavior of the solar system?[END_REF], the onset of the PETM and the ETM-2 hyperthermals have astronomically calibrated ages of 55.930 Ma and 54.050 Ma, respectively. Therefore, the bentonite layer on top of section CaH1 was deposited just before the onset of the ETM-2 hyperthermal. The top of the CaH1 section yielded dinoflagellate cysts of the genus Apectodinium spp. Sporadic occurrences of Apectodinium were also recorded in the partly equivalent section 2 of [START_REF] Parsons | Paleogene fungi from the Caribou Hills[END_REF]. The global Apectodinium acme coinciding with the negative carbon isotope excursions [START_REF] Crouch | The Apectodinium acme and terrestrial discharge during the Paleocene-Eocene thermal maximum: new palynological, geochemical and calcareous nannoplankton observations at Tawanui, New Zealand[END_REF][START_REF] Sluijs | Subtropical Arctic Ocean temperatures during the Palaeocene/Eocene thermal maximum[END_REF]) is a secondary marker of the PETM [START_REF] Aubry | The Global Standard Stratotype-section and Point (GSSP) for the base of the Eocene Series in the Dababiya section (Egypt): Episodes[END_REF]. The occurrence of the species Apectodinium augustum, which is restricted to the CIE interval, is furthermore a useful criterion for identification of the PETM [START_REF] Aubry | The Global Standard Stratotype-section and Point (GSSP) for the base of the Eocene Series in the Dababiya section (Egypt): Episodes[END_REF][START_REF] Sluijs | Arctic late Paleocene-early Eocene paleoenvironments with special emphasis on the Paleocene-Eocene thermal maximum[END_REF]. However, in accordance with McNeil and Parsons (2013), we did not find Apectodinium augustum in the Caribou Hills sections that we studied and the few Apectodinium spp. occurrences found do not account for an acme. In a section located stratigraphically 100 m below CaH1, [START_REF] Ioannides | A preliminary palynological study of the Caribou Hills outcrop section along the Mackenzie River, District of Mackenzie: Current Research, Part A[END_REF] reported Apectodinium parvum, which may co-occur with Apectodinium augustum in Beaufort-Mackenzie wells but ranges considerably higher in the early Eocene [START_REF] Damassa | Correlation of Paleogene dinoflagellate assemblages to standard nannofossil zonation in North Atlantic DSDP sites[END_REF][START_REF] Mcneil | The Paleocene-Eocene thermal maximum in the Arctic Beaufort-Mackenzie Basin-Palynomorphs, carbon isotopes and benthic foraminiferal turnover[END_REF]. The carbon isotope measurements performed on both CaH1 and CaH2 sections can be used to refine this age assignation, if the obtained values are not biased by changing organic matter source or abundance. In both sections, neither linear nor logarithmic correlation is seen between δ 13 Corg and TOC, suggesting that carbon isotope values are not influenced by the amount of remaining organic carbon. To be sure that the correlation coefficients are significant, the autocorrelation of each proxy (clay minerals, organic matter parameters, and carbon isotopes) was tested with the R software (R Core Team, 2013), and shows no autocorrelation in the data (acf function, critic value <0.2). This lack of correlation suggests that the recorded changes in δ 13 Corg values are mostly not attributable to changing organic sources or different stages of organic matter degradation, except for the coarser sand beds (e.g., CaH1 at 20.95 m), where very low δ 13 Corg values possibly reflect enrichment of lignin or other more resistant 13 Cdepleted components [START_REF] Benner | Depletion of 13C in lignin and its implications for stable carbon isotope studies[END_REF][START_REF] Wynn | Rayleigh distillation and the depth profi le of 13C/12C ratios of soil organic carbon from soils of disparate texture in Iron Range National Park, Far North Queensland[END_REF][START_REF] Foreman | Evaluation of paleoclimatic conditions east and west of the southern Canadian Cordillera in the mid-late Paleocene using bulk organic δ13C records: Palaeogeography[END_REF]. The δ 13 Corg values at CaH1 are slightly lower than those recorded at CaH2, which could be associated with the interval showing lower δ 13 C values in lower Eocene marine carbonate materials [START_REF] Vandenberghe | The Paleogene period[END_REF]. This association is supported by the warming trend recorded by clay mineralogy and pollen assemblages at CaH1 toward the top of the section with a slight increase in kaolinite highlighted by several kaolinite peaks, an increase of megamesotherm plant proportions, and an increase in the thermophilous/non-thermophilous taxa ratio. Given the above age constraints, we associate this record with the overall warming trend associated with the Early Eocene Thermal Optimum and proximity to the Eocene Thermal Maximum-2 (ETM-2) hyperthermal (Fig. 7).

Caribou Hills 2 (CaH2)

In the CaH2 section, pollen grains of megatherm plants are abundant, especially the mangrove element Avicennia (Fig. 5B), and thermophilous plants display a high diversity indicative of well-developed forest ecosystems and stable climatic warm conditions. Modern representatives of Avicennia are pollinated by animals and, thus, its pollen is expected to be present in very low quantity in coastal sediments, as suggested by analysis of modern coastal sediments offshore Ivory Coast and Thailand [START_REF] Caratini | Sédimentation palynologique actuelle sur le plateau continental de Côte d'Ivoire[END_REF][START_REF] Somboon | Palynological study of mangrove and marine sediments of the Gulf of Thailand[END_REF]. A depleted mangrove with Avicennia showing pollen percentages similar to those observed in the CaH2 section are recorded in early mid-Miocene coastal marine sediments controlled by river input from the Northwestern Mediterranean region [START_REF] Bessedik | Le couple récif-mangrove à Sant Pau d'Ordal (Vallès-Pénédès, Espagne), témoin du maximum transgressif en Méditerranée nord occidentale[END_REF][START_REF] Jiménez-Moreno | Middle Miocene latitudinal climatic gradient in Western Europe: evidence from pollen records: Palaeogeography[END_REF]. The Miocene context for pollen transport resembles that of the studied Caribou Hills sections: the active pollen transport by a river may explain the relatively high percentages of Avicennia (Table 3). In present-day surface sediments offshore Indonesia, the percentage of Avicennia pollen varies with distance from shoreline, and percentages matching with those of Caribou Hills must be interpreted as indicating a nearby mangrove [START_REF] Engelhart | Mangrove pollen of Indonesia and its suitability as a sea-level indicator[END_REF]. Moreover, microspores of the freefloating freshwater fern Azolla [START_REF] Collinson | The ecology of Cainozoic ferns[END_REF] are abundant at the top of the section. A massive Azolla bloom is recorded in the central Arctic Ocean and in all Nordic Seas in strata of latest early Eocene age (49 Ma;[START_REF] Van Der Burgh | The freshwater fern Azolla (Azollaceae) from Eocene Arctic and Nordic Sea sediments: New species and their stratigraphic distribution[END_REF]chron 22n-21r;[START_REF] Barke | Coeval Eocene blooms of the freshwater fern Azolla in and around Arctic and Nordic seas: Palaeogeography[END_REF], interpreted as reflecting a massive freshwater discharge into the Arctic Ocean [START_REF] Brinkhuis | Episodic fresh surface waters in the Eocene Arctic Ocean[END_REF][START_REF] Collinson | A new species of the freshwater fern Azolla (Azollaceae) from the Eocene Arctic Ocean[END_REF][START_REF] Greenwood | How wet was the Arctic Eocene rain forest? Estimates of precipitation from Paleogene Arctic macrofloras[END_REF][START_REF] Barke | Coeval Eocene blooms of the freshwater fern Azolla in and around Arctic and Nordic seas: Palaeogeography[END_REF] due to elevated precipitations in adjacent drainage areas. The duration of this Azolla bloom is estimated to have been ~800 ky [START_REF] Brinkhuis | Episodic fresh surface waters in the Eocene Arctic Ocean[END_REF] or 1.2 My [START_REF] Barke | Coeval Eocene blooms of the freshwater fern Azolla in and around Arctic and Nordic seas: Palaeogeography[END_REF]. These comparisons suggest that the CaH2 sediments were deposited during the termination of the Early Eocene Climate Optimum and more precisely just before and at the onset of the Azolla event (49 Ma), in agreement with the late early to early middle Eocene age suggested by [START_REF] Parsons | Palynology of Paleogene strata in the Caribou Hills[END_REF] for his section 6, which is partly equivalent to CaH2 (Fig. 2). These dating propositions are also supported by δ 13 Corg data from CaH2, which, similarly to CaH1, show no correlation with TOC, HI, or the type of organic particles described in the palynofacies section. The high frequency, small amplitude changes in δ 13 Corg values recorded at CaH2 could, however, reflect differences in the composition of the organic matter or even annual change in temperature and humidity, which are known to influence δ 13 Corg values of modern C3 plants [START_REF] Jahren | The Arctic forest of the middle Eocene[END_REF]. In this regard, it is noteworthy that the negative δ 13 Corg excursion recorded between 37 and 40 m coincides with a marked increase in the proportion of pollen grains belonging to megatherm plants (Fig. 5B). As several studies have shown that higher precipitation decreases the δ 13 C values of C3 plants [START_REF] Stewart | 13 C natural abundance in plant communities along a rainfall gradient: a biological integrator of water availability[END_REF][START_REF] Diefendorf | Global patterns in leaf 13C discrimination and implications for studies of past and future climate[END_REF][START_REF] Kohn | Carbon isotope compositions of terrestrial C3 plants as indicators of (paleo) ecology and (paleo) climate[END_REF][START_REF] Schubert | A summertime rainy season in the Arctic forests of the Eocene[END_REF], this excursion could be attributed to a transient period of higher precipitation. Nevertheless, significant changes in mean annual precipitation are not supported by our climatic reconstruction or in similar studies for the same interval [START_REF] Eldrett | Increased seasonality through the Eocene to Oligocene transition in northern high latitudes[END_REF][START_REF] Eldrett | A seasonality trigger for carbon injection at the Paleocene-Eocene Thermal Maximum[END_REF] or for the middle to late Paleocene in the Northwest Territories [START_REF] Greenwood | How wet was the Arctic Eocene rain forest? Estimates of precipitation from Paleogene Arctic macrofloras[END_REF]. Notwithstanding this excursion, the CaH2 succession shows average δ 13 Corg values that are significantly higher than in CaH1, and may thus correlate with marine strata showing higher values toward the termination of the EECO (Fig. 7). Altogether, these palynological and geochemical data support a late early Eocene age for the CaH2 section.

Reconstruction of Depositional Environments

The high amounts of smectite, the absence of illite-smectite mixed layers, the stratigraphically invariant illite crystallinity, and the low average Tmax values (<435 °C) recorded in the two studied sections indicate limited burial and thermal overprinting and allow the interpretation of organic matter and clay assemblages in terms of paleoenvironmental conditions. The pseudo Van Krevelen diagram (Fig. 4) and the palynofacies suggest that the organic matter of the two sections mostly corresponds to terrestrial higher plants debris. The numerous plant cell structures observed in the palynofacies suggest a rather well-preserved organic matter and the relative proximity of terrestrial organic sources. Few dinoflagellate cysts have been observed along the sections (stars on Fig. 5B), indicating some marine incursions and relative proximity of the coastline. In CaH1, dinoflagellate cysts are found at 0.65, 18, 20.5, 21.4, 22.4, 23.4, and 24.8 m. In CaH2, only three samples exhibit very few dinoflagellate cysts, at 3.7, 7.6, and 51.2 m, and thus indicate a more proximal environment compared to CaH1. This is in agreement with the higher proportions of small equant opaque phytoclasts (50 μm) and pollen grains in the palynofacies assemblages of CaH1 when compared to those of CaH2. Small equantopaque phytoclasts particles have a high buoyancy [START_REF] Mcarthur | Characterization of deep-marine channel-levee complex architecture with palynofacies: An outcrop example from the Rosario Formation[END_REF] and thus suggest a longer aquatic transport for organic particles at CaH1 when compared to CaH2.

The pollen assemblage of CaH1 is characteristic of littoral forests with freshwater to marine water swamp trees (Glyptostrobus-Taxodium-type) and fern (Osmunda). Freshwater marshes were also present as illustrated by Myrica, Cyperaceae, and the association of Cuphea, Menyanthaceae, Myriophyllum, Oenotheraceae, Restionaceae, and Typha. Very few specimens of brackish marshes species are found, such as many Amaranthaceae species and Tamarix. The hinterland was probably covered by evergreen broadleaf forests (Engelhardia, Amanoa, Distylium, Cupressaceae, Sapotaceae) and with increasing altitude by deciduous broadleaf mesophilous forests (Alnus, Carpinus, Carya, Liquidambar, Populus, Pterocarya, Quercus, Ulmus, Zelkova), both growing in a humid climate. Proximity of a higher elevation massif is indicated by Cathaya, Podocarpus-type, Cedrus, Tsuga, Abies, and Picea. The CaH1 section is characterized by massive silty clay sediments without root traces. Such fine-grained deposits more or less enriched in organic matter are indicative of low-energy conditions that we attribute to outer levees or swamps of deltaic plains. The centimeter-thick laminated sand located at 20 m could correspond to a crevasse splay deposit. In summary, the CaH1 succession corresponds to a proximal, dominantly freshwater littoral swamp environment with intermittent connections with marine areas. Pollen assemblage of CaH2 suggests a freshwater to marine water swamp forest (Glyptostrobus-Taxodium-type), but in a warmer climatic context as illustrated by the presence of megatherm elements (Amanoa, Bombax, Buxus bahamensis-type, Croton, Fabaceae Mimosoideae, Icacinaceae, Resedaceae, and Theaceae). The occurrence of Avicennia pollen on CaH2 is indicative of a mangrove swamp environment, and confirms the presence of this genus in the Eocene Arctic as initially suggested by pollen data of coastal plain sediments from the New Siberia Islands [START_REF] Suan | Subtropical climate conditions and mangrove growth in Arctic Siberia during the early Eocene[END_REF]. Mangrove is a salt tolerant forest environment containing diverse assemblages of trees and shrubs and developing along sheltered tropical to subtropical coasts [START_REF] Mao | Key to mangrove pollen and spores of southern China: an aid to palynological interpretation of Quaternary deposits in the South China Sea[END_REF], where the average annual temperature of seawater is at or above 18 °C [START_REF] Chapman | Mangrove Vegetation: Vaduz[END_REF]. Avicennia is the mangrove genus that can support relatively cooler temperatures compared to other mangrove taxa and may extend into marginally subtropical areas in the Southern Hemisphere [START_REF] Duke | Genetic structure and evolution of species in the mangrove genus Avicennia (Avicenniaceae) in the Indo-West Pacific[END_REF][START_REF] Quisthoudt | Temperature variation among mangrove latitudinal range limits worldwide: Trees[END_REF]. Modern Avicennia mangrove forest develops on soft, exposed mudflats of low to middle intertidal zones inundated by high tides, although they may occasionally develop in the high intertidal zone [START_REF] Duke | Aging Rhizophora seedlings from leaf scar nodes: a technique for studying recruitment and growth in mangrove forests[END_REF]. This taxon indicates a mangrove environment inundated by high tides, as recorded in modern coastal areas of Southeast Asia [START_REF] Giesen | Mangrove guidebook for Southeast Asia: Bangkok[END_REF]. A depleted mangrove with Nypa probably out of its typical latitudinal climate zone was already evidenced at the same time period in Tasmania at a similar paleolatitude [START_REF] Pole | Eocene Nypa from Regatta Point[END_REF][START_REF] Carpenter | Neartropical Early Eocene terrestrial temperatures at the Australo-Antarctic margin[END_REF]. This mangrove swamp would have been located near a fluvial channel belt, as suggested by the occurrence of unconsolidated sand and gravel deposits surrounded by fine sand, silty clay, and coal deposits. These fining upward sequences are characteristic of channel, levee, and backswamp deposits. At the base of the CaH2 section, a sand layer shows parallel laminations and as suggested for a bed of similar aspect at CaH1, is here interpreted as crevasse splay deposit. As for CaH1, the pollen data of CaH2 indicate that the hinterland forest environments contained evergreen broadleaf megathermic and mega-mesothermic taxa (Amanoa, Sapotaceae, Engelhardia, Distylium) and mesophilous taxa (Alnus, Carpinus, Carya, Liquidambar, Populus, Pterocarya, Quercus, Ulmus, Zelkova). Pollen grains of deciduous plants and mesomicrotherm and microtherm taxa are more abundant in CaH2 than in CaH1, which appears surprising given that the CaH2 section also yielded numerous and diversified pollen grains of megatherm and mega-mesotherm elements (Fig. 6B). The higher abundance of these elements could be explained by an uplift phase of the neighboring hinterland heights; this uplift would have either decreased the distance between the main areas of these plants and the studied coastal deposits, or increased the areas favorable to their development. This hypothesis is also supported by higher amounts of chlorite minerals in CaH2, which could be related to an erosion pulse of hinterland areas. Alternatively, the higher abundance of meso-microtherm and microtherm plants in CaH2 could reflect a shift in the fluvial network draining higher altitudinal source areas. Source areas for Beaufort-Mackenzie Basin sediments from the Eocene consisted of the northern Cordillera, including the Ogilvie and Richardson Mountains, the Porcupine River drainage basin, and the northeastern Brooks Range [START_REF] Duk-Rodkin | Tertiary-Quaternary drainage of the pre-glacial Mackenzie Basin: Quaternary International[END_REF]. During the Paleocene and Eocene these mountains underwent major uplift due to the accretion of terranes along the west margin of North America. In summary, the CaH2 succession was in a proximal, freshwater swamp and mangrove environment with a more limited connection to the open sea that could be inundated only during high tides.

Paleoclimatic Conditions

Clay assemblages from the CaH1 section are dominated by smectite, suggesting drier and/or colder and/or seasonally contrasted climate (Chamley, 1989) when compared to CaH2. CaH1 shows kaolinite proportions mostly around 10%, indicating colder and/or dryer paleoclimate compared to the high kaolinite proportions (40% in average) recorded in CaH2. This is in agreement with the lower proportions of megatherm and megamesotherm plants in CaH1 and the climate estimates that indicate lower temperatures for CaH1 than CaH2 (Fig. 6b). Conversely, a warm and humid paleoclimate is suggested by the higher kaolinite content in CaH2, in agreement with high values of the S+K/I+C ratio, low values of S/K, high pollen proportions of megatherm and mega-mesotherm plants (Amanoa, Bombax, Buxus bahamensistype), and high temperature estimates, especially from the base to 41.2 m. Furthermore, at CaH2, the Most Likely Values of MAT (between 21.3 and 22 °C) and MAP (between 1200 and 1360 mm) are compatible with kaolinite formation in soils. In both sections, illite and chlorite are present in high proportions, reflecting a high detrital input and active erosion of the hinterland, consistent with the record of pollen of meso-microtherm and microtherm plants that were most likely sourced from hinterland heights. The two sections also record short-term variations in temperature or precipitation. In CaH1, as stated previously, the higher proportion of kaolinite toward the top of the section is associated with a lower S+K/I+C ratio that could reflect higher physical erosion due to higher precipitation. The higher proportion of mega-mesotherm plants, thermophilous taxa to the top of the section also attest for a warming trend from the base to top of CaH1 section. In CaH2, the proportion of pollen grains of megatherm plants decreases toward the top of the section above 42 m in association with higher proportions of mesotherm and microtherm plants. Above 41.5 m, temperatures become lower than during the first part of the sequence, especially for the coldest month mean temperature, equivalent to CaH1 temperatures (Fig. 5B). This interval shows lower kaolinite proportions and S+K/I+C ratios lower than 1. This interval could therefore be linked to the onset of the cooling phase following the EECO starting near 50 Ma (Fig. 7). Later in this region, SST rose from 10 to 13 °C at the termination of the Azolla event. Along CaH2, high S+K/I+C and S/K ratios at 3.2, 27, and 40.75 m indicate an increase in smectite contents and thus likely point to transient intervals of cooler or dryer climate conditions, in line with the drop in abundance of megatherm taxa in corresponding strata and with the occurrence of Azolla at 40 m. As Azolla is thought to indicate higher precipitation [START_REF] Brinkhuis | Episodic fresh surface waters in the Eocene Arctic Ocean[END_REF], the increase in smectite at 40 m likely reflects a cooler climate. However, Avicennia is still present, and climatic estimates do not show changes. The record of warm and wet climate conditions in the Mackenzie Delta of the Eocene Arctic is consistent with other studies of coeval strata in Ellesmere Island, Axel Heiberg Island, Greenland, Norwegian Sea, and New Siberian Islands (Table 4; Fig. 1), where both macro and micro remains of humid temperate swamp forests and lush rain forests have been also reported [START_REF] Kalkreuth | The petrology, organic geochemistry and palynology of Tertiary age Eureka Sound Group coals[END_REF][START_REF] Eldrett | Increased seasonality through the Eocene to Oligocene transition in northern high latitudes[END_REF][START_REF] Eldrett | A seasonality trigger for carbon injection at the Paleocene-Eocene Thermal Maximum[END_REF][START_REF] Eberle | Seasonal variability in Arctic temperatures during early Eocene time[END_REF][START_REF] Greenwood | How wet was the Arctic Eocene rain forest? Estimates of precipitation from Paleogene Arctic macrofloras[END_REF][START_REF] Harrington | Arctic plant diversity in the Early Eocene greenhouse[END_REF][START_REF] West | Was the Arctic Eocene 'rainforest' monsoonal? Estimates of seasonal precipitation from early Eocene megafl oras from Ellesmere Island, Nunavut[END_REF][START_REF] Suan | Subtropical climate conditions and mangrove growth in Arctic Siberia during the early Eocene[END_REF]. Mild warm-temperate climates with mean annual temperatures (MAT) around 8-15 °C have been estimated for the early and middle Eocene around the Arctic. Warm and wet summers prevailing in the Arctic during the early and early middle Eocene, with the WMMT between 18 and 26 °C and MAP values around 1200 mm, and mild winters (CMMT between 0 and 13 °C; [START_REF] Basinger | Early Tertiary vegetation of Arctic Canada and its relevance to paleoclimatic interpretation[END_REF][START_REF] Greenwood | Eocene continental climates and latitudinal temperature gradients[END_REF][START_REF] Brinkhuis | Episodic fresh surface waters in the Eocene Arctic Ocean[END_REF][START_REF] Sluijs | Warm and wet conditions in the Arctic region during Eocene Thermal Maximum 2[END_REF][START_REF] Greenwood | How wet was the Arctic Eocene rain forest? Estimates of precipitation from Paleogene Arctic macrofloras[END_REF][START_REF] Eberle | Seasonal variability in Arctic temperatures during early Eocene time[END_REF][START_REF] Eberle | First record of Eocene bony fi shes and crocodyliforms from Canada's Western Arctic[END_REF][START_REF] Suan | Subtropical climate conditions and mangrove growth in Arctic Siberia during the early Eocene[END_REF] enabled the growth of a forest assemblage of deciduous and evergreen conifers. Such MAT and humidity conditions are presently found in regions such as southeast North America or southeast China, which have plant associations similar to those recorded in Eocene Arctic forests [START_REF] Kumagai | Growth-ring analysis of Early Tertiary conifer woods from the Canadian High Arctic and its paleoclimatic interpretation: Palaeogeography[END_REF]. Warm climates during the early and latest early Eocene are recorded in many localities around the world, for example in Wyoming, Texas, Canada, England, and high latitude Antarctica [START_REF] Robert | Antarctic subtropical humid episode at the Paleocene-Eocene boundary: Clay-mineral evidence[END_REF][START_REF] Andreasson | Winter and summer temperatures of the early middle Eocene of France from Turritella δ 18 O profiles[END_REF][START_REF] Wilf | Late Paleocene-early Eocene climate changes in southwestern Wyoming: Paleobotanical analysis[END_REF][START_REF] Pross | Persistent near-tropical warmth on the Antarctic continent during the early Eocene epoch[END_REF][START_REF] Smith | Early Eocene plant diversity and dynamics in the Falkland flora[END_REF][START_REF] Jacques | An early-middle Eocene Antarctic summer monsoon: evidence of 'fossil climates[END_REF][START_REF] Gushulak | Paleoclimate and precipitation seasonality of the Early Eocene McAbee megaflora[END_REF], with temperatures up to 10-30 °C higher when compared to those characterizing modern high-and midlatitude areas. These warm temperatures are associated either with high precipitation, especially in high-and mid-latitudes [START_REF] Robert | Paleocene and Eocene kaolinite distribution in the South Atlantic and Southern Ocean: Antarctic climatic and paleoceanographic implications[END_REF][START_REF] Wilf | Late Paleocene-early Eocene climate changes in southwestern Wyoming: Paleobotanical analysis[END_REF] or monsoon-like climates at lower latitudes [START_REF] Shukla | Cool equatorial terrestrial temperatures and the South Asian monsoon in the Early Eocene: evidence from the Gurha Mine, Rajasthan, India[END_REF][START_REF] Bougeois | A high resolution study of trace elements and stable isotopes in oyster shells to estimate Central Asian Middle Eocene seasonality[END_REF][START_REF] West | Was the Arctic Eocene 'rainforest' monsoonal? Estimates of seasonal precipitation from early Eocene megafl oras from Ellesmere Island, Nunavut[END_REF][START_REF] Gushulak | Paleoclimate and precipitation seasonality of the Early Eocene McAbee megaflora[END_REF][START_REF] Spicer | Asian Eocene monsoons as revealed by leaf architectural signatures[END_REF]. All these records provide strong evidence for a reduced latitudinal temperature gradient during the Eocene [START_REF] Zachos | Evolution of early Cenozoic marine temperatures[END_REF][START_REF] Andreasson | Temperature seasonality in the early middle Eocene North Atlantic region: Evidence from stable isotope profiles of marine gastropod shells[END_REF][START_REF] Sluijs | Subtropical Arctic Ocean temperatures during the Palaeocene/Eocene thermal maximum[END_REF] and tend to invalidate recent hypotheses postulating that the Eocene high-latitude warmth is an artifact of paleothermometer miscalibrations or oxygen isotope diagenetic overprint [START_REF] Ho | Flat meridional temperature gradient in the early Eocene in the subsurface rather than surface ocean[END_REF][START_REF] Bernard | Burialinduced oxygen-isotope reequilibration of fossil foraminifera explains ocean paleotemperature paradoxes[END_REF]. Our results thus provide further evidence for climate conditions that are exceptionally warmer and more humid than today in the high latitude regions of the northern hemisphere during the early Eocene, and, in this regard, can be used as critical climatic constraints for future climate modeling work, similarly to the modelling works by [START_REF] Huber | The Early Eocene equable climate problem revisited: Climate of the Past[END_REF] and [START_REF] Huber | Eocene monsoons[END_REF], which used climate condition reconstructions determined from studies on sections. Furthermore, we note that pollen grains of Avicennia are, for the first time, recognized in the Canadian margin of the Arctic Basin. This strengthens the idea that the exceptional Arctic warmth of the Eocene was a phenomenon characterizing both coastal marine and atmospheric environments [START_REF] Suan | Subtropical climate conditions and mangrove growth in Arctic Siberia during the early Eocene[END_REF]. Beyond these climatic implications, the occurrence of Avicennia pollen in the Arctic Canada, through analogy with modern representatives of this genus that only spread through a marine transport of their propagules [START_REF] Quisthoudt | Temperature variation among mangrove latitudinal range limits worldwide: Trees[END_REF], implies a marine connection between the Arctic Basin and the mid-latitude oceans. Avicennia pollen is known today to be transported a very short distance (<2 km) from its source [START_REF] Somboon | Palynological study of mangrove and marine sediments of the Gulf of Thailand[END_REF][START_REF] Behling | Studies on Holocene mangrove ecosystem dynamics of the Bragança Peninsula in north-eastern Pará, Brazil[END_REF], thus indicating the presence of this mangrove around the Arctic Ocean in the early Eocene. According to several authors [START_REF] Ellison | Origins of mangrove ecosystems and the mangrove biodiversity anomaly[END_REF][START_REF] Duke | Global distribution and genetic discontinuities of mangrovesemerging patterns in the evolution of Rhizophora: Trees[END_REF][START_REF] Lo | Phylogeographic pattern of Rhizophora (Rhizophoraceae) reveals the importance of both vicariance and long-distance oceanic dispersal to modern mangrove distribution[END_REF] the mangrove evolved around the Tethys Sea during the Late Cretaceous and spread over long distances through oceanic currents. The only gateways through which this dispersal may have occurred during this time interval were the Greenland-Norway Seaway connecting the Atlantic and Arctic oceans [START_REF] Roberts | Climate sensitivity to Arctic seaway restriction during the early Paleogene[END_REF]Fig. 1) and the Turgai Strait connecting the Tethys and Arctic Oceans (open at least until the Lutetian according to Iakovleva et al., 2001, and until the Oligocene according to [START_REF] Rögl | Mediterranean and Paratethys palaeogeography during the Oligocene and Miocene[END_REF]. We can therefore assume that oceanic currents flowing toward the Arctic Ocean existed through at least one of these gateways at this time. Their influence on the mean salinity of the Arctic Ocean was, however, probably limited, as shown by the development of the Azolla freshwater fern. Moreover, the occurrence of Avicennia on both the Siberian and Canadian margins also implies the existence of oceanic currents within the Arctic Ocean during the Eocene. Although further data are needed to constrain its detailed patterns, such an oceanic circulation could have been similar to the modern Beaufort gyre and Transpolar drift, which has been previously suggested to have developed only after at 46 Ma [START_REF] Backman | Expanding the Cenozoic paleoceanographic record in the Central Arctic Ocean: IODP Expedition 302 Synthesis[END_REF]. The occurrence of Avicennia pollen grains on both sides of the Arctic, as highlighted by our new results, suggests that the onset of such a trans-Arctic oceanic circulation may have occurred as early as during the early Eocene.

CONCLUSIONS

This multiproxy study leads to several conclusions dealing with the climatic variations during the early Eocene in the Canadian margin of the Arctic Ocean: • The two studied sections (CaH1 and CaH2) of the Caribou Hills represent high-latitude records, near 75°N, of early and late early Eocene ages, respectively. The new U-Pb zircon (ID-TIMS) age of 54.304 (±0.029) Ma of the bentonite layer in section CaH1 proves deposition of the sediments just before the Early Eocene Thermal Maximum 2 hyperthermal. CaH1 and CaH2 sediments were deposited under a subtropical climate and warmer conditions during the late early Eocene, respectively, as shown by consistent mineralogical and palynological records. For the two sections, the paleoclimatic reconstruction indicates most likely values of MAT between 18 and 22 °C, CMMT between 6 and 14 °C, WMMT between 26 and 28 °C, MAP between 1200 and 1400 mm, and available moisture between 75 and 95%, compatible with kaolinite formation in soils.

• The two portions of the Caribou Hills outcrop represent two delta plain environments with littoral forests including swamps in a hinterland landscape dominated by evergreen broadleaf forest and with increased altitude by deciduous forest.

• We report the first occurrence of Avicennia pollen in the Canadian Arctic, which indicates that mangrove environments developed along Canadian Arctic coastlines during the early Eocene. Such mangrove environments, recently recognized in Eocene deposits of the New Siberia Islands [START_REF] Suan | Subtropical climate conditions and mangrove growth in Arctic Siberia during the early Eocene[END_REF], imply that warm and wet conditions were widespread along the Arctic coast during the early Eocene. These Avicennia pollen records also imply that oceanic currents connected Siberian and Arctic Canadian margins as early as during the early Eocene.

FIGURESFigure 1 .

 1 FIGURES

Figure 2 .

 2 Figure 2. Composite section of the Caribou Hills outcrop showing the sections studied in the past 30 years by various authors and the interpreted long-term sedimentary sequences (after Parsons, 2000; Price et al., 1980).

Figure 3 .

 3 Figure 3. Results of U-Pb zircon analyses (ID-TIMS) of the bentonite layer (sample CASE15_YU067) in section CaH1. Cathodoluminescence photograph illustrates the zircons from analyzed bentonite layer.

Figure 4 .

 4 Figure4. Rock Eval results of the two Caribou Hills sections. HI-Hydrogen Index in mg HC/g TOC; OI-Oxygen Index in mg CO2/g TOC; S2-mg HC/g rock; Tmax-degree of thermogenic maturation of the kerogen in °C. Type I-lacustrine organic matter. Type II-marine organic matter. Type III-terrestrial organic matter. Filled circles-CaH1; Open circles-CaH2. For the two sections, the organic matter corresponds to a Type III and is immature.

Figure 5 (

 5 Figure 5 (A) From left to right, for the two sections of Caribous Hills: results of Rock Eval pyrolysis, clay mineralogy, and stable carbon isotope geochemistry. (B) From left to right, for the two sections of Caribou Hills: results of palynofacies, stable carbon isotope geochemistry, palynology, and climatic parameters. b-bentonite; p-Pediastrum; A-Avicennia; stars-dinoflagellate cysts; dots-Azolla microspores.

Figure 6 .

 6 Figure 6. Photographs of pollen grains using a scanning electronic microscope (CaH2 section, sample 8). (A-D) Avicennia (Acanthaceae). (A) Equatorial view showing an aperture (scale bar = 2 μm). (B) Detail of the reticulate ornamentation in the mesocolpium (scale bar = 1 μm). (C) Polar view showing the small polar area (scale bar = 2 μm). (D) Detail of the reticulate ornamentation at the pole (scale bar = 1 μm). (E-F) Alfaroa-type (Juglandaceae). (E) Polar view showing the three pores (scale bar = 2 μm). (F) Sculpture made of dense spinules (scale bar = 1 μm). (G-H) Engelhardia (Juglandaceae). (G) Polar view showing the three pores (scale bar = 2 μm). (H) Sculpture made of scattered and uniformly distributed spinules (scale bar = 1 μm). (I-J) Glyptostrobus-Taxodium-type (Glyptostrobus) (Cupressaceae). (I) Pollen showing a small papilla (scale bar = 2 μm). (J) Detail of its composite heteromorphic ornamentation (large verrucae integrated within numerous smaller ones) (scale bar = 1 μm). (K-L) Cathaya (Pinaceae). (K) Proximal view showing the haploxylonoid outline of sacci and corpus (scale bar = 8 μm). (L) Low spinules on the corpus (scale bar = 1 μm). (M-N) Pinus haploxylon-type (Pinaceae). (M) Profi le view showing the haploxylonoid outline of sacci and corpus (scale bar = 8 μm). (N) Verrucae on the corpus (scale bar = 1 μm).

Figure 7 .

 7 Figure 7. Proposed stratigraphic correlation of the two studied sections (green lines) based on bentonite dating, the pollen assemblages, and the climatic clues of the total set of analyses. The Paleocene-Eocene hyperthermal events are indicated in red on the stable isotopes curve. PETM-Paleocene-Eocene Thermal Maximum; ETM2-Eocene Thermal Maximum 2; ETM3-Eocene Thermal Maximum 3; EECO-Early Eocene Climatic Optimum; MECO-Middle Eocene Climatic Optimum. The Azolla event is indicated in black.

ACKNOWLEDGMENTS

This study was funded by TOTAL through the project Groupement Recherche Industrie (GRI) zones péri-arctiques. We acknowledge Florence Savignac, François Fourel, Ingrid Antheaume, Omar Boudouma, and Vincent Grossi for their help with analytical procedures. We wish to thank two anonymous reviewers for their constructive remarks.