Late Quaternary Deglaciation, Glaciomarine Sedimentation and Glacioisostatic Recovery in the Rivière Nastapoka Area, Eastern Hudson Bay, Northern Québec

Lajeunesse, Patrick; Allard, Michel

doi:https://doi.org/10.7202/010331ar

Introduction

Late Quaternary glacial geology of Hudson Bay has been relatively well described from research carried out at sea using marine geophysical methods supported by a few cores (Josenhans et al., 1988; Bilodeau et al., 1990; Josenhans and Zevenhuizen, 1990; Gonthier et al., 1993; Zevenhuizen, 1996; Hill et al., 1999). However, in the terrestrial coastal sectors of eastern Hudson Bay, only a small number of papers (Hillaire-Marcel, 1976; Hillaire-Marcel et al., 1981; Allard and Seguin, 1985; Parent and Paradis, 1994; Parent et al., 1995) have addressed the regional history of ice retreat. Recently, Lajeunesse and Allard (2002, 2003) reported on the presence of a series of emerged ice-proximal glaciomarine fans that were deposited during a short period of stabilization of the ice margin in the Rivière Nastapoka area in association with the first phase of deglaciation of southeastern Hudson Bay. The regional associations of these submarine fans with other glacial landforms, local deglaciation dynamics and the rapid glacioisostatic rebound remain nevertheless poorly described in this sector of Hudson Bay. The influence of basin topography on ice retreat and glaciomarine sediment deposition dynamics, which has been well demonstrated in other regions of the world (e.g., Meier and Post, 1987; Stewart, 1991; Seramur et al., 1997; Ó Cofaigh, 1998), remains unexamined in regards to the deglacial marine transgression in the region (Tyrrell Sea).

In this paper, we describe from on-land evidence the geomorphic, sedimentary, stratigraphic and chronologic record of Late Quaternary ice recession, marine transgression, sedimentation, and glacioisostatic rebound along a stretch of coast located between the northern tip of Lac Guillaume-Delisle and Rivière Richard, eastern Hudson Bay. This integrated approach leads to point out the mechanisms and chronology by which deglaciation and glaciomarine sedimentation occurred in the context of a falling relative sea-level (RSL) due to glacioisostatic recovery.

Study area

The study area is located between 56°-57° N and 76° 40'- 75° 30' W, on the eastern coast of Hudson Bay (Fig. 1). Bedrock geology of most of the study area consists of Archean granitic-gneiss (Stevenson, 1968; Avramtchev, 1982). In the southeast sector of the study area, important outcrops of Proterozoic volcano-sedimentary rocks such as basalt, quartzite and dolomite make up high cuesta ridges that face east and enclose Lac Guillaume-Delisle. The area has four physiographic sectors: 1) the coastal slope, which is a narrow, 4-7 km wide, coastal fringe from the actual seashore to the foot of the hills, following a gently inclined <3° slope; 2) the Nastapoka Hills, a hill range parallel to the coast located between ~4-14 km inland reaching mean elevations of 280-320 m above sea level (asl) and maximum elevations of 420 m asl; 3) an inner basin, i.e. an area of lower and less rugged topography with elevations generally <200 m asl and located between 14-23 km from the seashore; and 4) the Hudsonian Plateau, a peneplain with elevations varying between 200-300 m asl extending from ~23 km inland to well beyond the eastern limit of the study area. The most important river, Rivière Nastapoka, has two waterfalls along its course within the study area. The first falls are located 1 km from the shore and create a drop of >40 m. The second falls are located 7 km inland and are 46 m high. These two falls divide the Rivière Nastapoka into three major segments: 1) a short estuary downstream from the first falls; 2) a second segment between the first and the second falls at 40-80 m asl; and 3) a third one upstream from the second falls at >125 m asl.

Methods

The landforms associated with raised coastlines were first mapped from interpretation of 1:10 000 and 1:20 000 aerial photographs. The marine limit was defined at numerous locations, mainly with the upper limit of the washing zone of ablation till on the tops and sides of rocky hills (Allard and Seguin, 1985); several glaciofluvial deltas are also at the same elevation. Downwards from this maximum level reached by the Tyrrell Sea at deglaciation time, other features indicative of past sea levels include raised boulder fields, many of them having wave-made boulder ridges, sandy beaches and the surface of perched deltas along river courses. Their elevation was measured at numerous places during field traverses. Natural sections of glacial, marine, deltaic, coastal and eolian sediments, exposed along gullies and rivers were described and logged (Fig. 2). Landform and sample elevations were determined with the help of a precision altimeter (0.1 m) using a provincial geodetic benchmark located near base camp as reference. Marine shells, wood fragments and organic paleosoils were collected for ¹⁴C dating. Ages are reported here in conventional δ¹³C normalized ¹⁴C and calibrated (cal.) years BP (Tabl. I; Fig. 3), in contrast to measured ¹⁴C dates in Lajeunesse and Allard (2003). The dates were calibrated to calendar years by using the CALIB 4.3 program (Stuiver and Reimer, 1993). Calibration of shell dates was made using the ΔR regional mean value of 310 years for the southeastern Hudson Bay area (¹⁴C_meas. – contemporaneous surface ocean ¹⁴C age) (Barber et al., 1999). The intercept calendar year belonging to the range with the highest probability was retained as the calibrated date.

Geomorphology

Late-Wisconsinan ice flow

Fluted till, striae, chattermarks, roches moutonnées and drumlins indicate a general westward (270°) flow that correspond to the regional ice flow reported in Parent et al. (1995). Flutings are located >200 m asl and concentrated in the NE sector of the study area, although some are present north of Rivière du Nord (Fig. 4). They occur in sectors where the Tyrrell Sea waters did not invade the land, and above marine limit in areas where it did (Figs. 4, 5). Bedrock outcrops exposing fresh and intact striae and associated erosion features are relatively rare since they were either eroded by periglacial coastal processes, affected by weathering or buried by Holocene sediments.

Ice grounding and recession positions

Temporary stabilization phases of the ice margin are recorded by ice-contact landforms and sediments. These ice front stillstands occurred in marine and terrestrial environments and led to the construction of frontal, ablation and marine DeGeer moraines as well as ice-contact glaciomarine and glaciofluvial sediment bodies. The distribution of these features is almost exclusively restricted to sectors of relatively higher relief (Fig. 4). Frontal moraines are generally parallel to sub-parallel to the actual seashore and occur on some summits of the Nastapoka Hills, where they are clearly visible above marine limit. Frontal moraines also occur in the bottom of a wide erosion gully in a sandur-delta (S44) and near the downstream portion of an esker in the southeastern part of the area (Fig. 4). Grounding lines and morainal banks occur in a general north-south trend along the Nastapoka Hills and form a distinctive mark between sandy sediments westward and silty-clay sediments eastward.

Ablation till consists of sub-angular to sub-rounded boulders and gravel with a very thin and irregularly distributed blanket of poorly sorted sand. Boulder and gravel often lie directly on glacially moulded bedrock outcrops. Generally, the finer-grained material veneer locally smoothes the landscape by filling out small depressions in bedrock. This moraine was deposited in the ice marginal zone as the ice sheet was retreating. It occurs on most hilltops of the Nastapoka Hills and is almost completely eroded and reworked below the marine washing zone, which forms a clear distinctive mark of marine limit.

DeGeer moraines were observed at 180 m asl and ~15 km north-east of the Rivière Nastapoka mouth. They are oriented parallel to the coast of Hudson Bay. As the marine limit in this sector is ~200-210 m asl, these moraines were formed under ~30 m of water. The mean spacing of 108 m between the ridges (Allard and Seguin, 1985), if taken as the yearly amount of retreat of the ice front, would suggest that it took ~330 yr for the ice front to retreat from the Nastapoka Hills to the inner limit of the marine transgression, some 35 km inland. This rate is less than the 150 m and 217 m/yr measured further south in the James Bay area by Vincent (1977). This slower rate of retreat for the Nastapoka sector can be attributed to more pronounced local topographic variations, marked by deep valleys and hills higher than marine limit (Allard and Seguin, 1985).

Ice-contact glaciomarine sediment accumulations occur as submarine fans, sandur-deltas and deltas. These predominantly sandy bodies were constructed near the ice front during stabilization phases of the eastward retreating ice margin. These deposits form two drift belts: one along the Nastapoka Hills range and the other at the inner limit of marine inundation on the Hudsonian Plateau (Fig. 4). The first belt consists of ice-contact submarine fans that were deposited in a glaciomarine sedimentary environment during a period of high RSL (Lajeunesse and Allard, 2002) and a pause in ice recession (Lajeunesse and Allard, 2003). Some fans become coalescent and form aprons ~5 km wide on the western slopes of the Nastapoka Hills. In many cases, sediments composing these fans are constricted in elongated bedrock basins and become coalescent downslope. The fans start upstream at ~200-220 m asl along the hill range (Lajeunesse and Allard, 2002) and slope down to below the present day sea level underneath Hudson Bay waters (Lavoie et al., 2002). Holocene gully erosion has exposed >20 m of sand, silt, gravel beds containing abundant fossil shells. High energy jets of glacial meltwater carved large potholes of 1-1.5 m of diameter and 2-3 m deep in bedrock side-walls.

The second drift belt consists of ice-contact sandur-deltas located in the eastern part of the study area (Fig. 4). They are located at 195-205 m asl and are associated with both the marine limit and the inland postglacial marine overlap zone. Sediments composing sandur-deltas were deposited in both subaerial and submarine environments and contain topsets, foresets and bottomsets. The subaerial sandurs bear networks of braided channels. The length of the outwash extends on more than 12 km in some cases and are often pitted with kettles. Many sandur-deltas are located at the downstream end of long and sinuous eskers indicating a flow from east to west. A cross-section analysed in S44 located in the downstream sector of a sandur-delta exposes predominantly undisturbed fine sand and silt beds deposited by suspension-settling. Sediment deformation and minor loading structures occur, indicating rapid sedimentation. There are no fossil shells at this site, which can be attributed to high inputs of fresh water in a restrained embayment. The distance between two moraines of the DeGeer type buried under this section indicates a very slow glacial retreat of 15-20 m/yr, i.e. a near stillstand. This second drift belt either marks a pause or a drastic slowdown of glacial retreat as the glacier front became separated from the effect of sea waters. Mass wasting became regulated only by melt without any influence of calving.

Table I

**Radiocarbon and calibrated ages, with 68.3 % (1σ) confidence limits, from fossil shells, driftwood, terrestrial wood, peat and paleosoils**

Marine limit and uplifted shorelines

Marine limit is at 248 m asl near the Rivière Nastapoka mouth, ~220-227 m asl north of Rivière Nastapoka and 230 m asl at Rivière Devaux (Fig. 5). East of the Nastapoka Hills, the plane of the marine limit slopes down to 204 m asl at the maximum inland reach of the Tyrrell Sea where it forms an irregular paleo-coastline pattern (Allard and Seguin, 1985). Marine invasion went as far as 56 km inland from the present day coastline at Lac Kakiattualuk, 48 km along Rivière Nastapoka, 23 km along Rivière Devaux and 67 km in the southern portion of the study area in the depression of Lac Guillaume-Delisle.

Postglacial emergence is recorded by boulder ridges and pavements, raised beaches, perched deltas and terraces as well as coastal gelifraction of bedrock. Boulder ridges and pavements occur at any elevation below marine limit. Between 248 and 200 m asl, where bedrock has absolutely no cover of fines, they represent, together with gelifracted outcrops, the main geomorphic record of land emergence. Boulder ridges generally form series of linear or lightly curved ridges on sidewalls of small valleys. Boulder pavements occur on flat or gently sloping surfaces. Coastal periglacial erosion features (i.e. landforms made by frost-riving and ice-push processes on the shore) can be found up to the very marine limit. Below this limit, the smooth glacially moulded surface of bedrock was intensely disrupted by frost action in a wet coastal environment into angular clasts which were then carried away by drift ice and storm waves (Allard and Tremblay, 1983b). A large number of these clasts are present in boulder barricades and pavements and are mixed with subrounded allochtonous clasts derived from the washing out of glacial drift (Dionne, 1978). Coastal periglacial processes are still active along the present day shoreline. Series of beach ridges occur below 200 m asl, where they often form almost regular staircases in sand deposits. Vertical spacing between these ridges ranges from 2 to 7 m, as in the Lac Guillaume-Delisle region, just south of the study area (Hillaire-Marcel, 1976). These ridges are common features on the surface of the emerged deposits of the western slopes of the Nastapoka Hills. At least 60 beach crests occur between 0 and 100 m asl on a distance of 2.5 km near the Rivière Nastapoka mouth.

Sedimentology and stratigraphy

Raised late quaternary deposits

With the exception of till and moraines, Late Quaternary sediments of the study area record depositional and reworking processes that took place during the initial phase of marine transgression and the subsequent falling RSL. Sediments exposed in sections were deposited in ice-proximal (De1) and ice-distal (De2) glaciomarine, paraglacial (Pa) and emergence (Em) depositional environments (Tabl. II; Fig. 6). De1 deposits are associated with the two stillstands of the ice front (Nastapoka Hills glaciomarine fan complex and sandur-delta belt), De2 sediments are associated with the retreat between these two ice front positions. Pa sediments are associated with fluvial erosion and sedimentation as rivers made new courses by down-cutting and re-depositing existing sediments. A coarsening upward sequence representing the Holocene falling RSL is observed in most sections.

Deglacial deposits: ice-proximal sediments (De1)

Ice-proximal sediments (De1) are dominated by sand (Fig. 6). These deposits are concentrated on the coastal slope between 80 and 220 m asl, where they form thick and extensive submarine fans (Lajeunesse and Allard, 2002), and in the eastern part of the study area between ~185-220 m asl, where they form deltas and sandur-deltas. The largest deposits occur near the Rivière Nastapoka mouth and along the river in the eastern end of the studied sector. Smaller and restricted accumulations of De1 sediments occur in the inner basin and are associated with ice-contact deltas. In the Nastapoka Hills, De1 deposits extend from 0 to 4 km downstream from the former ice front. Exposures of De1 are up to 23 m high in sections S32, S39-S40 and S42. De1 are exposed in S1a, S6-S9, S12-S19, S21-S35, S38-S40 and S42-S44 (Fig. 2). De1 is divided in members De1a and De1b, which correspond respectively to FA I and FA II in Lajeunesse and Allard (2002). De1a sediments consist of interbedded and interlaminated sand, silt and silty clay and De1b consist of massive diamicts. Sand layers (De1a) located nearer to the former ice front are coarse and massive and interfinger with diamicts (De1b) while those located at relative distances from it are fine and stratified. They are interlayered with silty sand and silt and silty clay massive and laminated beds. General bed inclination result from westward flowing paleo-currents. Deposition of these facies was dominated by suspension-settling and submarine gravity flows. Ice rafted debris are also incorporated, but in a lesser proportion. Sediment mass-movement processes were very active in the Nastapoka Hills submarine fans, where steep slopes contributed largely to sediment instability. Dewatering structures are present in many exposures and suggest rapid sedimentation. Sediment deformation and semi-consolidation caused by glaciotectonism occur in S9. Exposures of ice-proximal deposits in S8, S9 and S43 are located near former grounding lines and S32, S39 and S44 are close to frontal moraines.

Table II

**Depositional processes and environments associated to the emerged marine deposits in the Rivière Nastapoka area**

Deglacial deposits: ice-distal sediments (De2)

De2 sediments consist of interbedded and interlaminated silty clay and sand (Fig. 7a) and are, in some cases, similar to FA II silty clay. They form the lowermost unit in exposed sections. Silty clay is the predominant size fraction. Extensive zones of De2 occur in the inner basin below 200 m asl and are generally restricted to areas located along Nastapoka, Sheldrake and du Nord rivers. These deposits can be as much as 85 m thick in deep bedrock basins (Allard and Seguin, 1985). The considerable thickness of this unit in valleys of the shield necessarily result from an abundant and rapid glaciomarine sedimentation, considering the short period of only 300-400 yr for the ice to retreat and the rapid emergence rate at deglaciation time. Such silty clay deposits are one of the most distinctive constituent of meltwater derived sediments in temperate and subpolar glaciomarine environments (Cowan and Powell, 1990). De2 deposits, exposed in S1a, S20 and S45, are in some cases deformed, showing wavy parallel and chaotic bedding and normal faults. Outsized clasts and fossil shells are sparsely distributed in this unit. Some layers of De2 are composed of blue-grey silty clay with a strong sulphur odour (Fig. 7b). These beds are varved and contain black speckles. The blue-grey color of this silty clay indicate a direct glacial input (Powell, 1983). Thin black layers are monosulphide layers (Elverhøi et al., 1980, 1983), which are typical in glaciomarine high accumulation rate conditions (Görlich, 1986). These layers were preserved by rapid sedimentation that limited the effect and rate of their transformation into pyrite and their metabolisation to organic matter (Berner, 1984; Stevens, 1990). They correspond to annual organic varves probably formed during a time-lag between the important spring diatom bloom and clastic summer sedimentation (Elverhøi et al., 1980). Textural and structural characteristics indicate that De2 sediments were transported by turbid meltwater plumes emanating from the glacial margin (hypopycnal flow) and were deposited by suspension-settling (Elverhøi et al., 1980, 1983; Powell, 1981; Görlich, 1986; Cowan and Powell, 1990). The presence of thin sand layers within silty clay beds may be associated with turbulence caused by shifting of the meltwater jet (McCabe et al., 1993) or by submarine gravity flows. In S1a, sand beds associated with De1 sediments occur above older ice-distal De2 laminated silty clay. This can be explained by a remobilization of the original ice-proximal sand deposition and by their progradation on silty clay.

Paraglacial deposits (Pa)

Paraglacial sediments (Pa) are coarse-grained fluvial deposits transported to the marine environment under conditions of abnormally high transport rates associated with the rapid ablation of ice sheets in the hinterland (Church and Ryder, 1972; Syvitski et al., 1987, Syvitski, 1991). During this period, sediment yields were high due to fluvial-related erosion and transportation of large volumes of pre-deposited glacigenic and marine sediments (Syvitski, 1991). In the study area, Pa sediments were deposited as the glacial margin was retreating eastward towards central Québec-Labrador, inland from the marine limit. Pa deposits form extensive zones of thick sand bodies on the lower western slopes of the Nastapoka Hills. They are concentrated near the Nastapoka, Sheldrake and Umiujaq river mouths, where they are exposed between 0 and 100 m asl. They are considerably more extensive near the Rivière Nastapoka mouth. Pa deposits consist of well sorted and bedded sand bodies associated with fan deltas and comprise from base to top three members: 1) interbedded and interlaminated silty clay and sand (Pa1); 2) interbedded sand, silt and silty clay (Pa2); and 3) laminated and massive well sorted fine to medium sand beds (Pa3).

Pa1 sediments form sequences of massive and laminated silty clay irregularly alternating with laminated and massive fine sand, silty clay being the predominant member. They are well exposed and more concentrated near the Nastapoka (S2, S3) and Umiujaq (S41) river mouths. These sediments overlie ice-distal beds with unclear contacts. They are easily distinguishable from ice-distal sediments by their much higher sand composition and concentrations of sand laminations. This unit is located at relatively lower altitudes (<53 m asl) and underlies Pa2 rhythmic beds with a gradation contact. However, in S41, they are found in sharp contact with Pa3. Contacts between Pa silty clay and sand beds (Fig. 7c) are sharp, erosive or gradational. Pa1 were deposited during periods of predominant suspension settling and bedload. Variations in overall relative bed thickness of fines and sand is attributed to changes of the sediment source. Increasing thickness of sand beds in the upper parts of this unit corresponds to increasing proximity of the sediment source and higher energy due to glacioisostatic uplift. The muddy texture and the position of Pa1 below Pa2 rhythmic beds suggest that they are bottomsets. Terrestrial organic matter and wood fragments occur in the deposits.

Pa2 sediments are composed of rhythmically interbedded well sorted to poorly sorted fine sand and poorly sorted silty sand and/or poorly sorted silty clay (Fig. 7d). They are restricted to relatively low elevations (<70 m asl) and are mainly concentrated near the Nastapoka (S2-S5), Devaux (S37) and Umiujaq (S41) river mouths. They overlie Pa1 beds with sharp or gradation upper contact. In areas unaffected by erosion, Pa2 beds always underlie Pa3 beds. They generally present a slight westward inclination. The rhythmites resemble closely to varves. Sand beds become progressively thicker and more numerous than the silty sand and silty clay beds in the upper parts of this unit. These deposits contain no fossil shells and are usually not bioturbated, although ichnofossils were observed in some exposures. Processes responsible for deposition of Pa2 are suspension, bedload and recurrent small-scale gravity flow. A coarsening-upward and gradation upper and lower contacts indicate progradation due to glacioisostatic uplift. Rhythmic bedding of fines with sand found in these beds are typical of a prodeltaic environment and are interpreted here as transitional deposits between foresets and bottomsets. The generally rare bioturbations suggests rapid sedimentation.

Figure 6

Representative logs of cross-sections through ice-proximal (De1), ice-distal (De2), paraglacial (Pa) and emergence (Em) lithofacies associations, Rivière Nastapoka.

Coupes stratigraphiques représentatives dans des sédiments glaciaires proximaux (De1), distaux (De2), paraglaciaires (Pa) et d’émersion (Em), dans la rivière Nastapoka.

-> Voir la liste des figures

Pa3 are composed of laminated and massive well sorted fine to medium sand (Fig. 7e). They occur at elevations lower than 100 m asl near the Nastapoka (S3-S5), Devaux (S36, S37) and Umiujaq (S41) river mouths. These beds form the main body of a very large raised delta 63 m asl near the Rivière Nastapoka mouth. Pa3 reaches maximum thickness of 16 m in S5. Coarser sand and gravel beds occur in this unit and usually tend to be poorly sorted. Laminated beds are generally inclined (25°) and oriented WNW. Inclined (25°) and WNW oriented <20 cm thick rhythmically interbedded sand and silt of Pa2 occur with sharp contacts within these deposits. Pa3 beds overlie Pa2 beds with gradation or sharp contacts. Upper contacts of Pa3 with emergence beds are erosive or scoured. The high-angle dipping, prograding sand and occasional gravel beds of Pa3 sediments are characteristics of foresets. Paleocurrent data also suggest a sediment source in the coastal hills flowing towards the Hudson Bay coastline, which is similar to modern river flow.

Emergence deposits (Em)

Emergence deposits were laid in place in a wave affected environment when the glaciomarine and paraglacial sediments were exposed to coastal processes during emergence. Em consists of interbedded sand and gravel (Fig. 7f). They are exposed in almost all the studied sections where they form the uppermost unit. They tend to become thicker at lower elevations. Em are composed of well to poorly sorted massive to crudely bedded gravel, granules and coarse sand alternating with well to poorly sorted laminated and massive coarse to fine sand. They are generally 0.75-2 m thick but can vary from a few cm to 5 m thick. Lower contacts are erosive and scoured and often form lag deposits. Presence of cross-bedding, plane-parallel lamination and gravel imbrications indicate bottom-current processes. Discontinuous and wavy sand laminations, frequent heavy mineral layers and very good sorting are typical of beach deposits. Swash-and-backswash and ice rafting processes disturbed, eroded and re-transported previous laid deposits, thus affecting their original character. Presence of under-consolidated very fine sand and diffuse laminations in the upper parts of this unit is attributed to some eolian reworking that most likely started in the beach environment. This process formed dunes on the surface of raised sandy accumulations. Lulls in eolian activity lead to development of soils. Paleosoils were preserved by the return of eolian activity (Filion et al., 1991; Saint-Laurent and Filion, 1992) and provided minimum ages for sediment emergence.

Figure 7

-> Voir la liste des figures

Chronology

A great part of the radiocarbon dated material was recovered in exposures located near the Rivière Nastapoka mouth (Tabl. I; Fig. 3). The oldest radiocarbon dates are concentrated along the Nastapoka Hills and provide minimum ages for deglaciation. Many of these dates were recovered in De1a, in many cases associated with submarine gravity flows lithofacies. An age of 8580 cal. yr BP (UL-2346; S47) was obtained on paired valves Macoma balthica at 4 m asl in De2b sediments. This sample, recovered in ice-distal glaciomarine sediments deposited during the Nastapoka Hills stillstand, yielded the oldest radiocarbon date of the study area. However, since ΔR varies greatly in eastern Hudson Bay, reaching in some cases exceptionally high values of nearly 600 years (Barber et al., 1999), it is possible that sample UL-2346 could overestimate by a few hundred years the real age of the Nastapoka stillstand. In order to present a more cautious interpretation of deglacial events in the sector, this sample is disregarded in this study, even if its age is still enclosed within the extremes of the 1σ age range for the sudden drainage event of Glacial Lake Ojibway of 8160-8740 cal. yr BP (Barber et al., 1999). A younger age of 8300 cal. yr BP was provided by in situ Mya truncata, Hiatella arctica and Macoma balthica recovered 1 km downstream from the upper terrace of an ice-contact submarine fan at 213 m asl (UL-1932; S29; 145 m asl). Although the age of this sample could also slightly overestimate the timing of deglaciation due to a possible abnormally high ΔR value, it is contemporaneous with the 8280 cal. yr BP age for the drainage of Glacial Lake Ojibway reported in Barber et al. (1999) for southeastern Hudson Bay. UL-1932 is thus used to date deglaciation and as a minimal age for the 248 m marine limit in the Rivière Nastapoka area. This age also supports the model of a contemporaneous ice-stillstand in southeastern Hudson Bay and eastern James Bay reported by Lajeunesse and Allard (2003). In the same section, 9 m above UL-1932, a younger age of 7840 cal. yr BP was provided by a sample of in situ Mya truncata and Macoma balthica (TO-7411; S29; 154 m). A slightly older age of 7920 cal. yr BP (UL-1951; S19) was yielded by valves of Hiatella arctica at 108 m asl in a massive silty sand bed of De1a associated with another fan exposed along Stream A. An age of 8180 cal. yr BP was yielded by a driftwood sample recovered in Em sediments at 5 m asl along the shoreline (CAMS-37437; S10). This sample was recovered in modern sediments and in association with a sample of Chlamys islandicus that gave an age of 4630 cal. yr BP (CAMS-37434; S10). A few hundred metres south of Rivière Nastapoka, in the hill range area, valves and fragments of Hiatella arctica found on the surface of mudboils (De2a) gave an age of 7750 cal. yr BP (UL-1740; S45; 115 m). This sample recovered in ice-distal deposits slightly east of the hills indicate that the Nastapoka Hills were already deglaciated by 7.8 ka cal. BP.

Fossil shells recovered in highly deformed sediments provide maximum and minimum ages for submarine gravity flow events. Fragments of Hiatella arctica from De1a disturbed laminated sand at 110 m asl exposed along a stream near Rivière Nastapoka gave an age of 7700 cal. yr BP (UL-1728; S7). Fragments of Mytilus edulis and Hiatella arctica found at 132 m asl in De1 silty clay exposed along Stream B yielded an age of 7500 cal. yr BP (UL-1722; S27). In similar deposits of the same section, in situ Hiatella arctica and Mytilus edulis gave an age of 7360 cal. yr BP (UL-1718; 112 m). Near Rivière Devaux, fragments of Hiatella arctica and Mytilus edulis found in an exposure perched along a steep valley wall in De1a at 146 m asl gave an age of 7130 cal. yr BP (Beta-121551; S40). In “Vallée-des-Trois”, near Lac Guillaume-Delisle, fragmented Hiatella arctica and Mytilus edulis at 100 m asl in an exposure of De1b diamict bed deposited by submarine mass movement dated 7020 cal. yr BP (Beta-125152; S43). Younger radiocarbon dates in similar deposits were provided by paired valves of Mya truncata at 109 m asl that dated 6530 cal. yr BP (TO-7409; S13) and paired valves of Mytilus edulis and Hiatella arctica at 103 m asl that dated 6570 cal. yr BP (UL-1710; S17). Nearby S17, valves of Hiatella arctica in highly disturbed De1 laminated sand at 97 m asl dated 6660 cal. yr BP (UL-1540; S23). The same age was provided by a sample collected by Allard and Seguin (1985) at 118 m asl in similar deposits near S23 (UQ-540). Fragments of remobilized Hiatella arctica and Mytilus edulis in a silty sand layer at 84 m asl gave an age of 6730 cal. yr BP (UL-1952; S38). This sample provides a maximum age for a submarine landslide that occurred in a high RSL. The radiocarbon database from De1 disturbed beds suggests that submarine mass-movements were still active in the submarine fans about 6.7 ka cal. BP, some 1500 years after their deposition in shallowing marine conditions.

Submerged driftwood collected at 51 m in Pa1 deltaic sediments near the Rivière Nastapoka mouth dated 5760 cal. yr BP (UL-1556; S3). The upper terrace of the delta is located at ~100 m asl. UL-1556 therefore provides a maximum age of 5.8 ka cal. BP for the formation of this delta terrace. Near the Rivière Devaux mouth, fragments and valves of mixed shells at 15 m asl in Pa3 gravel dated 4340 cal. yr BP (UL-1919; S36) and at 13 m asl in similar sediments dated 3870 cal. yr BP (UL-1917; S36). The stratigraphic position of UL-1917 under 10 m of sediments implies that an important sedimentation was still prevailing in relatively shallow water offshore after 3.9 ka cal. BP. An age of 3550 cal. yr BP (UL-1546; S1a) was yielded by undifferentiated shells fragments in Em alluvial bedded sand at 12 m asl. Paired valves of Hiatella arctica recovered at a lower elevation east of S43 in Pa2 sediments dated 3040 cal. yr BP (UL-1371).

Organic samples of buried paleosoils were collected at high elevations near the Nastapoka Hills. Paleosoil samples UL-1904, at 98 m asl, and UL-1878, at 136 m asl, gave ages of 5420 cal. yr BP (S14) and 4750 cal. yr BP (S46), respectively. These samples provide minimum ages for emergence of their enclosing sediments. Paleosoil samples UL-1903 (1960 cal. yr BP; S24; 150 m asl) and UL-1949 (730 cal. yr BP; S9; 133 m asl) developed long after emergence and are not representative of land uplift. However, a younger paleosoil sampled in the trough of an ice-wedge (polygon) at 12 m asl near the shore gave an age of 310 cal. yr BP (UL-1858; S1b) and provides a minimum age of emergence of the coastal beach gravel (Em) in which it developed. Two driftwood samples collected at 4 m asl under ~40-60 cm of beach sand (Em) in the vicinity of S1b yielded an age of 510 cal. yr BP (UL-1536; UL-1528), providing a minimal age for emergence of the deposits. Basal peat and wood samples recovered at 170 m asl over De2 sediments east of the Nastapoka Hills south of Rivière Nastapoka dated 6640 cal. yr BP (I-13130 and I-13131) (Allard and Seguin, 1985). Another sample collected by the same authors at Lac Sheldrake at 180 m asl dated 6910 cal. yr BP (I-13406), indicating that the underlying sediments were emerged by 6.9 ka cal. BP.

Discussion

Land uplift and emergence curve

Fossil shells and driftwood buried under marine or coastal sediments provide maximum ages for the emergence of their overlying material whereas paleosoils, wood and peat buried by eolian processes or organic accumulation provide minimum ages for emergence. Hence, radiocarbon ages on shells and driftwood are located below the RSL curve while those from paleosoils, terrestrial wood and peat are above it. The assemblages of fossil shells recovered are not typical of the shore zone and the emergence facies are not fossiliferous; therefore no single raised beach or shoreline could be dated, such as what was done on the Manitounuk Islands, south of the study area (Allard and Tremblay, 1983a). The available data allows only a constrained RSL curve to be deducted (Fig. 8).

The RSL curve is characterized by an inflection between 5 and 6 ka cal. BP. This inflection indicates that the emergence occurred in an abruptly decelerating mode and divides the curve into three phases: 1) before 6 ka cal. BP; 2) between 6 and 5 ka cal. BP; and 3) after 5 ka cal. BP. Most fossils shells reported in this study were recovered in sediments deposited or remobilized by submarine mass movements before 6.5 ka cal. BP. A very important part of the 248 m total Holocene emergence took place in the first phase of glacioisostatic rebound.

After reaching a maximum elevation of 248 m asl about 8.3 ka cal. BP in the Nastapoka Hills, the Tyrrell Sea regressed rapidly at a rate of ~7 m/100 yr during deglaciation. This initial uplift rate is similar to the previously reported maximum rate of 6 m/100 yr on Ottawa Islands by Andrews and Falconer (1969). This emergence was accompanied by extensive and numerous submarine gravity flows on the steep slopes of ice-contact submarine fans located on the structural slope west of the Nastapoka Hills. These events cannot be directly attributed to postglacial seismicity linked to glacioisostatic recovery since they occurred in a setting prone to mass-movements (Lajeunesse and Allard, 2002). Slope instabilities in soft sediments, possibly triggered by bottom wave disturbances in shallowing waters might have been a key factor for submarine landslides during this phase. Following this period of homogeneous uplift rates (upper segment of the RSL curve), emergence decelerated to reach rates of ~4 m/100 yr between 6 and 5 ka cal. BP (inflection segment of the curve). This phase was characterized by a period of vegetation development on emerged sediments, in accordance with the timing of forest-tundra colonisation inland (Payette, 1983). Since 5 ka cal. BP, ~85 m of uplift has occurred at rates of ~1.7 m/100 yr (lower segment of the curve). Our data suggest an emergence rate slightly slowing down to ~1.6 m/100 yr over the last 1000 years, which is higher that the 1.3-1.5 cm/yr at Manitounouk Islands (Allard and Tremblay, 1983a; Bégin et al., 1993), 1.4 cm/yr near Inukjuak (Gray et al., 1993) and 1.1 cm/yr at Lac Guillaume-Delisle (Hillaire-Marcel, 1976). However, more dates are needed on this segment of the curve in order to improve the resolution of the recent uplift rate.

Land emergence led to extensive incision by rivers and streams through sand and silty clay deposits, which remobilized and carried downstream important volumes of sediments. RSL fall exposed marine sediments deposited during the deglaciation to the action of coastal processes. Exposed bedrock outcrops were disrupted by gelifraction; ensuing large rock fragments were redistributed by drift ice and added to glacial erratics to form boulder pavements and ridges. Accumulations of sand deposits were reworked by swash-and-backswash and eolian processes to form successive beach ridges. The upper portion of emergence deposits was reworked by eolian processes that formed dunes and deflation surfaces. During the late phases of emergence, the sedimentary regime of the Rivière Nastapoka was strongly affected by the formation of two important falls near the hill range. These bedrock sills formed inland basins that restrained coarse material from reaching the seashore and reduced greatly sand transport to the Rivière Nastapoka mouth (Lavoie, 2000; Lavoie et al., 2002). Ice-distal glaciomarine silty clay (De2), which was previously buried under deltaic deposits, will probably be the predominant sediment emerging in the future.

Considering the very important postglacial emergence of the area (~250 m) and the very high rates of regression still prevailing today (~8.3 ka cal. BP following deglaciation), it is questionable if this uplift can be attributed totally to glacioisostasy. In fact, an abnormally cold and dense downwelling convective region in the earth’s mantle located below the North American craton (Mitrovica, 1997; Simons and Hager, 1997) may have partially acted to depress the earth’s crust and/or contributed to create a sensitive zone for isostatic depression during the last glaciation.

Deglaciation and sedimentation

Geomorphic, sedimentary and radiometric data indicate that deglaciation of the Rivière Nastapoka area was marked by four phases. During Phase 1, an ice stillstand occurred at 8.3 ka cal. BP along the Nastapoka Hills after a period of very rapid down-wasting and collapse over Hudson Bay (Lajeunesse and Allard, 2003). Large ice-contact glaciomarine fans were laid from the hills on the sea bottom, where now extends the coastal slope. At that time, theTyrrell Sea was at its maximum elevation and the emergence rate was extremely rapid. Sedimentation was dominated by suspension-settling from turbid glacial meltwater outflows, submarine gravity flows and ice rafting (Lajeunesse and Allard, 2002).

Phase 2 was marked by the unpinning of the ice front from the Nastapoka Hills and its rapid eastward retreat in the low elevation inner basin to the easternmost reaches of the Tyrrell Sea, about 35 km inland. This distance, together with the rate of retreat of 108 m/yr suggested by the DeGeer moraines, indicate that this phase lasted for 330 years. At the end of this phase, RSL had dropped to about 200-190 m asl. The high RSL in this deep basin exerted an important control on the glacial margin by promoting high rates of calving and fast ice flow and thus favouring ice sheet collapse and drawdown (e.g., Hugues, 1987; McCabe et al., 1993). Sedimentation was dominated by glacial meltwater outflows that deposited a thick silty clay drape from suspension-settling in submerged fjord-like valleys. During this period, the ice-contact glaciomarine fans of the Nastapoka Hills complex were intensively affected by submarine landslides. Timing control of this phase from radiocarbon dates is needed. However, considering that the end of Phase 2 sedimentation occurred in a RSL of ~200-190 m, it can be deducted from the uplift curve (Fig. 8) that this phase occurred prior to ~7.5 ka cal. BP.

Phase 3 was characterized by a second temporary pause or slowdown in ice retreat at the eastern margin of the Tyrrell Sea invasion, where elevations are slightly higher than in the inner basin. This phase began when RSL was still ~200-190 m asl, i.e. prior to ~7.5 ka cal. BP. The topography of this part of the study area and the rapid RSL fall caused this stabilisation by allowing the ice margin to become terrestrial, thus considerably reducing its ablation rate in shallower marine waters at the terminus. This phase lead to the deposition of a second regional drift belt consisting of sandur-deltas in association with important esker systems characterized by the same direction of flow than major rivers. Deposition of this sediment belt corresponds to either a complete stillstand of the ice margin or a change in sedimentation regime when the ice margin changed from a tidewater to a terrestrial environment. Sedimentation was dominated by glaciofluvial processes on the surface of sandurs and by turbid plumes in their submerged zones that deposited ice-proximal sand and silt. These plumes transported fines westward from these glaciofluvial systems and deposited silty clay along the fjord-like bedrock basins. According to shells recovered in De2 east of the Nastapoka Hills, ice-distal sedimentation continued until at least 7240 cal. yr BP (UQ-545). There is no evidence in support of a complete belt of ice recession for a set period on the Hudsonian Plateau. However, some far inland sectors were ice free at 6.6 ka cal. BP since organic remains had by then begun to accumulate in alluvial deposits over clay in the Rivière Sheldrake basin (Allard and Seguin, 1985). This phase was also early characterized by the development of a discontinuous permafrost (Allard and Seguin, 1987) and periglacial landforms (Lajeunesse, 2000).

Phase 4 was associated with the complete retreat of the ice margin from the area and its terrestrial ablation inland. Although the ice terminus was distant during this phase, sedimentation dynamics were still controlled by meltwaters. Sedimentation was dominated by fluvial processes that transported large amounts of glacigenic material along major river systems, which is typical of paraglacial sedimentary environments (Church and Ryder, 1972; Syvitski, 1993). During this period, sediment delivery was essentially controlled by river downcutting through pre-existing emerged deposits. Within the inner basin and where silty clay deposits are present, elevations of the third segment of the Rivière Nastapoka upstream from the second falls are between 125 and 160 m asl. This implies that the RSL was probably <160 m asl, considering that paraglacial deposits do not occur along the upper sections of Rivière Nastapoka. At the mouth of the paleo-Rivière Nastapoka, a thick and extensive fan delta developed <100 m asl. The RSL fall also lead to the incision, erosion and remobilization of ice-proximal and ice-distal glaciomarine and glacial sediments located along the fluvial system. These sediments were transported and deposited near the actual Rivière Nastapoka mouth as well as near other less important river mouths when the ice margin was tens and probably hundreds of kilometres distant. According to wood fragments (UL-1556) recovered in a fossil rich Pa1 bed, paraglacial sedimentation had already begun in some valleys by 5760 cal. yr BP. The fact that this rich fossil bed is overlain by more than 19 m of deposits points out abnormally high sedimentation rates that prevailed during this period, which is not representative of typical postglacial sedimentation conditions in Hudson Bay (Josenhans et al., 1988). This considerable sedimentation is attributed to the disintegration of Québec-Labrador ice and sediment reworking. Deposits similar to those reported here have been observed near the mouth of Grande Rivière de la Baleine (Gonthier et al., 1993).

Conclusion

This study integrating geomorphological, sedimentary, stratigraphic and radiometric data presents new information on the style and timing of deglaciation and glaciomarine deposition, and provides an additional uplift curve for eastern Hudson Bay. In the study area, deglaciation and sedimentation were controlled by topography and RSL fall and were marked by four phases: 1) a stillstand of the Québec-Labrador ice margin at 8.3 ka cal. BP along the Nastapoka Hills that deposited a belt of ice-contact glaciomarine sediments. During this phase, RSL was high and was falling very rapidly; 2) the unpinning of the margin from the hills followed by a retreat phase during a still rapidly falling RSL that deposited a thick drape of fine-grained glaciomarine sediments. DeGeer moraines suggest that this retreat took place in about 330 years; 3) a second stillstand along the margin of the higher relief of the Hudsonian Plateau which occurred when RSL had fallen to ~2 ka cal. BP. It deposited ice-proximal sandy sediments in sandur-deltas as well as ice-distal silty-clay in submerged valleys; and 4) a complete retreat of the Québec-Labrador ice from the study area at its disintegration further inland. This event generated meltwaters that eroded and transported important volumes of sediments along the coast in a decelerating fall of RSL. Paraglacial sedimentation along the coast of Hudson Bay was already underway at 5.8 ka cal BP. Further investigations are needed in the area in order to evaluate the timing and duration of the second stillstand or retreat slowdown phase along the Hudsonian Plateau (Phase 3) and the rates of retreat toward central Québec-Labrador (Phase 4).

Abstract

Résumé

Resumen

Introduction

Study area

Methods