Sediment core ARC4-BN05 collected from the Canada Basin, Arctic
Ocean, covers the late to middle Quaternary (Marine Isotope Stage – MIS – 1–15,
ca. 0.5–0.6 Ma) as estimated by correlation to earlier proposed Arctic
Ocean stratigraphies and AMS
The advances and decays of continental ice sheets play a significant role in the alteration of global climatic system, such as changing atmospheric circulations, creating large-area albedo anomalies and regulating the global sea level fluctuations (Clark et al., 1990). Reconstruction of the history of ice sheets is therefore important not only for a better understanding of feedbacks of the future climate change and its impact on regional climates but also for getting insights into the mechanisms of abrupt climate change.
Background map showing the location of core ARC4-BN05,the main Arctic rivers and the two major surface current systems: Beaufort Gyre (BG), Transpolar Drift (TPD) and Bering Strait (BS). Schematic geological map shows the distribution and prevailing lithology of the main terrains adjacent to the Arctic Ocean (Fagel et al., 2014).
Studies of Pleistocene glaciations around the Arctic Ocean dealt mostly with the late Quaternary history of the Eurasian Ice Sheet during Marine Isotope Stages (MIS) 1–6 (e.g., Svendsen et al.,2004; Larsen et al., 2006) or the Laurentide Ice Sheet (LIS) with a special attention to the Last Glacial Maximum (LGM) (e.g., Dyke et al., 2002; England et al., 2009). In addition to terrestrial data, studies of sediment cores from the Arctic Ocean are critical for comprehending the history of glacial advances and retreats (e.g., Polyak et al., 2004, 2009; Spielhagen et al., 2004; Stein et al., 2012; Kaparulina et al., 2016). However, the long-term history of circum-Arctic glaciations is still poorly understood, especially with respect to the western Arctic including the North America and East Siberia. A major impact of the North American ice sheets on circulation and depositional environments in the Arctic Ocean is indicated by various marine and terrestrial data (e.g., Phillips and Grantz, 2001; Stokes et al., 2005), whereas the East Siberian Ice Sheet (ESIS) remained largely hypothetical until recently. While terrestrial data are limited and remain to be better investigated (Grosswald, 1989; Basilyan et al., 2010; Ivanova, 2012), seafloor mapping data now provide ample evidence for the existence of considerable ice masses on the East Siberian margin (Niessen et al., 2013; Dove et al., 2014; Jakobsson et al., 2014, 2016), but the timing and extent of these glaciations is virtually unknown. Marine sedimentary records from the Arctic Ocean adjacent to the East Siberian margin could add valuable information to this intriguing paleoglaciological problem.
In this paper, we present a multiproxy study of glacial–interglacial changes during the late to middle Pleistocene based on sediment core ARC4-BN05 from the Canada Basin north of the Chukchi Plateau and east of the Mendeleev Ridge (Fig. 1). This location can be affected by the two main Arctic Ocean circulation systems, the Beaufort Gyre and the Transpolar Drift, which carry sea ice, icebergs, and sediment discharge from North America and Siberia, respectively. As this circulation along with sedimentary environments and sources varied greatly during the Pleistocene climate cycles, resulting variations in sediment delivery and deposition make for a valuable paleoclimatic record for the western Arctic. Biogenic proxies (such as foraminifers) have uneven and overall limited distribution in Arctic Ocean sediments, while the terrigenous component provides a more consistent material for paleoceanographic studies (e.g. Stein, 2008; Polyak et al., 2009). As sediments in the Arctic Ocean are primarily transported by sea ice and/or icebergs during glacial events, sediment composition yields important information not only on the provenance and transport pathways but also on the attendant glacial and paleoclimatic history (e.g. Spielhagen et al., 1997; Vogt et al., 2001; Knies et al., 2001). By using clay and bulk mineralogy, along with grain size and the content of major elements Ca and Mn, we reconstruct depositional environments and sediment provenance to provide clues to the history of western Arctic ice sheets and their interaction with the Arctic Ocean.
The Arctic Ocean is surrounded by land masses composed of an assortment of lithologies and situated in a variety of climatic, tectonic, and physiographic settings. Figure 1 depicts a schematic geological map showing the main terrains and associated lithologies (Fagel et al., 2014). The West Siberian Basin and East Siberian platform of the Eurasian continent are mainly composed of terrigenous sediment (Fagel et al., 2014). The Siberian (Putorana) traps constitute one of the largest flood basalts in the world (Sharma et al., 1992). The western Okhotsk–Chukotsk volcanic belt contains acidic to intermediate rocks, whereas intermediate to basic rocks are more characteristic of the eastern side (Viscosi-Shirley et al., 2003). The Kara Plate and the Taymyr foldbelt, as well as the Ural and Novaya Zemlya foldbelt, are mainly composed of intrusive and metamorphic rocks (Fagel et al., 2014).
The geology of outcropping terraines of Alaska mainly includes Canadian–Alaskan Cordillera, Brooks Range, and part of the North American platform containing mostly intrusive, metamorphic, and some clastic rocks (Fagel et al., 2014). The outcrops of the Canadian Arctic Archipelago are mainly composed of carbonate and clastic rocks (Phillips and Grantz, 2001; Fagel et al., 2014), whereas intrusive and clastic rocks are mostly characteristic for Greenland (Fagel et al., 2014).
Dissolved and suspended matter is transported to the Arctic Ocean by voluminous rivers, with the Lena and Mackenzie rivers being the largest on the Siberian and North American side, respectively, both directly affecting the western Arctic Ocean. The transported material is further distributed across the Arctic Ocean in water and/or ice by currents. The two main surface, wind-driven circulation systems are the clockwise Beaufort Gyre (BG) in the western Arctic and the Transpolar Drift (TPD), which carries water and ice from the Siberian margin to the Norwegian–Greenland Sea (e.g., Rudels, 2009). The strength and trajectories of these current systems may vary depending on changes in atmospheric pressure fields known as the Arctic Oscillation (Rigor et al., 2002).
Sedimentation in the Arctic Ocean is strongly controlled by sea ice that acts as sediment carrier but can also suppress sediment deposition under thick and persistent ice cover (Darby et al., 2006; Polyak et al., 2009). During glacial/deglacial events, multiple icebergs discharged into the Arctic Ocean from the termini of marine-based ice sheets and strongly affected sediment dispersal and deposition (e.g., Spielhagen et al., 2004; Polyak et al., 2009). Fine-grained sediments can also be transported by subsurface and deep-water currents, such as the Atlantic water (Winkler et al., 2002), but their role in the overall Arctic Ocean sedimentation is not well understood.
Gravity core ARC4-BN05 (referred to hereafter as BN05) was collected from the
Canada Basin in the vicinity of the Mendeleev Ridge (80
Index map showing the location of core ARC4-BN05 (yellow circle) and other cores from previous studies mentioned in this paper (red circles). LR, MR, AR, and NR are Lomonosov, Mendeleev, Alpha, and Northwind ridges, respectively; NGS is Norwegian–Greenland Sea. White lines show maximal Pleistocene limits reconstructed for Greenland, Laurentide, Eurasian, and East Siberian ice sheets (GIS, LIS, EAIS and ESIS; England et al., 2009; Svendsen et al., 2004; Niessen et al., 2013). Proposed flow lines for grounded ice sheets and ice shelves (red and white arrows, respectively) are after Niessen et al. (2013).
Lithostratigraphy and major proxies in core BN05: core
photograph with brown layer indices, lightness, Ca and Mn content (bulk
XRF – gray line, ICP-OES – black line), paleomagnetic inclination, magnetic susceptibility, planktic
foraminiferal abundance, and AMS
For age constraint within the radiocarbon range, accelerator mass
spectrometry
For grain-size analysis,
Coarse sediment
Concentrations of major elements, such as Ca and Mn, were determined on point samples by ICP-OES (iCAP6300) at the First Institute of Oceanography, SOA, China, following the standard procedures. For a more detailed downcore distribution, relative elemental abundances were obtained at 1 cm resolution using the Itrax XRF core scanner at the Polar Research Institute of China, set at 20 s count times, 10 kV X-ray voltage and 20 mA X-ray current. A good match of the ICP-OES and Itrax XRF data (Fig. 3) verifies the consistency of results. To account for the dilution effects on the background sedimentation, such as by coarse debris and biogenic processes, element contents were normalized to Al (e.g., März et al., 2011).
Color reflectance was measured using a hand-held Minolta CM-2002
spectrophotometer at 1 cm intervals. Only the grayscale lightness index (
A total of sixty 2 cm thick samples were collected at 4 cm intervals for
paleomagnetic measurements performed at the Paleomagnetism and Geochronology
Laboratory of the Institute of Geology and Geophysics, Chinese Academy of
Science. Magnetic susceptibility was measured using the KLY-4s Kappabridge
instrument. Subsequently, stepwise alternating-field (AF) demagnetization of
natural remanent magnetization (NRM) was conducted using the 2-G Enterprises
model 760-R cryogenic magnetometer (2G760) installed in a magnetically
shielded (
For bulk sediment mineralogy
Minerals actively sought in diffraction data analysis.
Samples for clay minerals determination (
To enhance the identification of potential contributions from various
sediment sources, and thus the interpretation of downcore proxy
distributions, principal component analysis (PCA) was performed in MATLAB
(MathWorks, 2014). To account for proxies potentially indicative of sediment
provenance and depositional processes and environments, PCA included all
analyzed mineralogical proxies along with main grain-size groups (clay,
silt, fine to medium sand (63–250
As common for sediment cores from the Arctic Ocean (e.g., Jakobsson et al., 2000; Polyak et al., 2004, 2009; Spielhagen et al., 2004; Stein et al., 2010a, b), core ARC4-BN05 displays distinct cycles in sediment color and composition expressed in interlamination of dark-brownish and lighter-colored grayish muds (silty clays, clay silts and sandy silt), with coarser dropstones occurring in several layers. The color cyclicity is approximated by changes in sediment lightness that largely mirrors the content of Mn (Fig. 3), consistent with other studies from the Arctic Ocean (e.g., Jakobsson et al., 2000; Polyak et al., 2004; Löwemark et al., 2008; Adler et al., 2009). We identify 18 distinctly brown units, from B1 to B18, characterized by elevated content of Mn (Fig. 3). Another prominent lithostratigraphic feature in the western Arctic Ocean, widely used for core correlation, is pink-white to whitish layers (PW) rich in detrital carbonates (e.g., Clark et al., 1980; Polyak et al., 2009; Stein et al., 2010a, b). We identify three major PW layers expressed both visually and in high Ca content (Fig. 3). Lower Ca peaks occur throughout the record without being clearly expressed in the core macroscopic appearance.
AMS
Foraminiferal abundances are generally high (mostly
The measured AMS
While detailed paleomagnetic investigation is not an objective of this
paper, we utilize the inclination data for an independent stratigraphic
constraint in line with earlier studies (e.g., Jakobsson et al., 2000;
Spielhagen et al., 2004; Polyak et al., 2009). Paleomagnetic inclination in
core ARC4-BN05 shows mostly positive values oscillating around
Other paleomagnetic parameters, such as magnetic susceptibility (MS), can provide additional correlation means (e.g., Sellén et al., 2010). Two prominent peaks in MS occur in the intervals between units B7/B8 and B10/B11 (Fig. 3).
Based on the results of grain-size analysis, sediment in core BN05 can be
generally classified as sandy, poorly sorted mud (e.g., Blott and Pye,
2012). Overall, silt and clay predominate grain-size composition
(33–60 and 23–61 %, respectively), but coarser particles also make a considerable
contribution, with up to
Relative weight contents of major clay mineral groups in
the clay fraction (
Grain-size distribution is mostly polymodal with three distinct major modes
centered at
Several core intervals contain large rock fragments
The clay assemblage in samples from core ARC4-BN05 mainly consists of illite, chlorite, kaolinite and smectite (Fig. 5). The illite group is overall the major constituent of the clay mineral fraction, ranging between 43 and 73 %. Its downcore distribution pattern is opposite to that of the three other major clay-mineral groups – kaolinite, chlorite, and smectite (mostly present in very low contents). These three groups largely co-vary except for some lithostratigraphic intervals, such as PW layers. Elevated content of these clay minerals is characteristic for grayish sedimentary units.
The bulk mineral assemblage in core ARC4-BN05 mainly consists of quartz, K-feldspar, plagioclase, calcite, dolomite and pyroxene (Fig. 5). Quartz is generally the most abundant mineral, ranging between 20 and 51 % and typically peaking in grayish sediment units. K-feldspar, plagioclase and pyroxene (mainly augitic) mostly co-vary, with peaks in gray units in the upper part of the core, but more in brown units in the lower part starting from unit B10. Calcite has a high content in brown units of the upper part and much lower values below unit B9. Dolomite distribution shows distinct peaks reaching up to 53 %, with the highest peaks occurring in or adjacent to the PW layers. Similar to other minerals, the pattern of dolomite distribution changes around unit B10, with maxima in thick gray units below and in thin interlayers within brown units above this stratigraphic level.
Stratigraphic correlation of core BN05 with PS72/392-5 (Stein et al., 2010a) based on sediment lightness, magnetic susceptibility, calcite and dolomite content. See Fig. 3 for other stratigraphic proxies and lithostratigraphy explanation. The vertical magenta bar indicates the position of the foraminiferal peak in B14–15.
Loading scores for variables used in the PCA.
Scores
The first five principal components identified by PCA with a Varimax rotation account for 77 % of the total variance, with relatively evenly distributed communalities (Table 3). This pattern presumably reflects a complexity of multi-proxy variables characterizing sedimentary environments and provenance, as well as their strong variability occurring over multiple climatic cycles. To further test the PCA performance, we have also run a factor analysis with the maximum likelihood extraction, which produced similar factor loadings and variance explained, thus indicating the robustness of the results.
As no single existing chronostratigraphic method can comprehensively
constrain the age of the Arctic Ocean Pleistocene sediments, the age model
for core ARC4-BN05 was developed by correlating multiple proxies (such as
paleomagnetic, foraminiferal, and lithological; see Figs. 3 and 5), combined
with
The two
Two
An abrupt increase in sediment age between closely spaced B1 and B2 in core ARC4-BN05 suggests a very condensed section or a hiatus between MIS 1 and MIS 3. This age distribution is common for the western Arctic Ocean and has been attributed to very low to no sedimentation due to a very solid sea-ice cover or an ice shelf during the LGM in MIS 2 (e.g. Polyak et al., 2009; Wang et al., 2013).
Below the range of
According to this approach, we identify foraminifera- and Mn-rich brown units B3–B7 and B8–B10 as warm substages of MIS 5 and 7, respectively (Figs. 3 and 6). This age assignment is corroborated by the prominent detrital carbonate peaks PW 2 and 1 near the bottom of MIS 5 and 7, respectively. Furthermore, the principal drop in paleomagnetic inclination in core ARC4-BN05 occurs in the lower part of MIS 7, consistent with many cores from the Arctic Ocean (e.g., Jakobsson et al., 2000; Spielhagen et al., 2004; Adler et al., 2009; Polyak et al., 2009). A solidly grayish, foraminifera- and Mn-poor unit separating brown units B2 and B3 is accordingly considered as related to glacial MIS 4, and a similar unit between B7 and B8 related to MIS 6. It is possible, however, that most of the fine-grained, grayish sediment was deposited during deglaciations following the actual glacial intervals, which may have been very compressed, similar to the LGM.
Stratigraphy below MIS 7 has been less investigated in prior studies, and is
more difficult to address due to often less distinct units and scarce to
absent foraminifers, probably resulting from stronger dissolution (e.g.,
Lazar and Polyak, 2016). Therefore the age model for the lower part of the
core is more tentative. Nevertheless, a prominent oldest foraminiferal peak
in units B14–B15 (Fig. 3) allows us to identify these units as MIS 11 by
comparison with other microfaunal records reported from the western Arctic
Ocean (e.g., Cronin et al., 2013; Polyak et al., 2013). While individual
species have not been counted in ARC4-BN05, predominant planktic
foraminifers in this peak are identifiable as
Distribution of various terrigenous components in Arctic sediment records
carries information on sediment sources and depositional environments, and
thus paleocirculation and changes in paleoclimatic conditions, such as
connection to other oceans and build-up/disintegration of ice sheets (e.g.,
Bischof and Darby, 1997; Krylov et al., 2008; Polyak et al., 2009; Stein et
al., 2010a, b; Yurco et al., 2010; Fagel et al., 2014). We utilize the data
on clay and bulk minerals along with the grain size and total Ca and Mn
distribution in core ARC4-BN05 to reconstruct changes in glacial conditions
and circulation in the western Arctic Ocean during several glacial cycles
extending to estimated ca. 0.5–0.6 Ma. In this work we capitalize on earlier
studies on the distribution of bulk and/or clay minerals in surface and
downcore Arctic Ocean sediments (e.g., Vogt, 1997; Stein, 2008; Krylov et
al., 2014; Zou, 2016), corroborated by more targeted provenance proxies,
such as radiogenic isotopes (Fagel et al., 2014; Bazhenova et al., 2017),
heavy minerals (Stein, 2008; Kaparulina et al., 2016), composition of coarse
debris (Bischof et al., 1996; Wang et al., 2013), and iron oxide grains
(e.g., Bischof and Darby, 1997; Darby et al., 2002). To optimize the PCA
results for clarifying relationships between various sedimentary proxies, we
plotted the leading PC loading scores as biplots in the PC 1–2 and PC 3–4
space (Fig. 7a). These plots help to identify several sedimentary
variable groups with high loadings (
To gain insight into stratigraphic changes in sedimentary environments and provenance, we plotted the distribution of the identified variable groups 1–4 using the combined downcore scores of PC 1–2 and PC 3–4 (Fig. 7b). A combination of the PC group composition and downcore variability provides useful guidance for interpreting major sedimentary controls and their stratigraphic evolution.
A variable, mostly multimodal distribution of grain size in core BN05
indicates multiple controls on sediment delivery and/or deposition. The
prevailing mode 1 at
Mode 2, centered at 7–7.5
Additionally, modes 1 and 2 make up a bimodal distribution in the lowermost
part of the core – mostly in estimated MIS 13/15 and near the bottom of
MIS 11. The predominant stratigraphic position in brown units makes
the glacigenic origin of this sediment unlikely. We hypothesize that this grain-size
pattern reflects a combination of “normal” interglacial environments with
winnowed silts deposited by downwelling of shelf waters enriched in dense
brines. Although no observational evidence exists for such waters
penetrating deeper than the halocline (
Coarse sediment, up to dropstones of several centimeters large, is a consistent feature in core BN05. In the apparent absence of strong current control on sedimentation, except for some shelf areas, and a pervasive presence of floating ice, coarse sediment in the Arctic Ocean is typically attributed to ice rafting, including sea ice and icebergs (e.g., Stein, 2008; Polyak et al., 2010, and references therein). Sedimentological studies in areas of sea-ice formation or melting and in ice itself indicate that sediment carried by sea ice in the Arctic Ocean is predominated by silt and clay, while coarser fractions are of minor importance (Clark and Hanson, 1983; Nürnberg et al., 1994; Hebbeln, 2000; Darby, 2003; Dethleff, 2005; Darby et al., 2009). Some studies suggest a higher content of sand in ice formed at the sea floor (anchor ice) (Darby et al., 2011), but the contribution of this source still needs to be evaluated. Furthermore, the role of sea ice on sedimentation in the Arctic Ocean is not clear for glacial intervals, when most of the sediment entrainment areas were exposed or covered by ice sheets. In contrast, in iceberg-rafted sediment, deposited mostly in glacial/deglacial environments, the content of large size fractions, from sand to boulders, is typically high, in excess of 10–20 % (Clark and Hanson, 1983; Dowdeswell et al., 1993; Andrews, 2000). Thus, elevated content of coarse sediment can be regarded as a good indicator of intense iceberg rafting. Such events are not probable during full interglacials, exemplified by modern conditions, but most likely occurred at times of instability and disintegration of ice sheets that extended to the Arctic Ocean in the past (e.g., Spielhagen et al., 2004; Stokes et al., 2005; Polyak et al., 2009).
In core BN05, coarse fractions (from coarse silt to sand) measured at
different sizes show very similar distribution patterns (Fig. 4a),
indicating the same predominant delivery mechanism, that is, iceberg
rafting. This pattern is reflected in a good correlation of fine to medium
sand (63–250
A common occurrence (separate or combined) of two coarse grain modes, around
85–90 and 400–450
One of the most robust sedimentary variable groups is distinctly
characterized by high loadings of dolomite along with total Ca content
(group 4: Fig. 7a, Table 3). Dolomite has been proposed as the main
contributor of Ca in sediment cores from the western Arctic Ocean, with an
especially high content in multiple coarse-grain peaks of detrital
carbonates (Bischof et al., 1996; Phillips and Grantz, 2001; Polyak et al.,
2009; Stein et al., 2010a, b). A high correlation (
The main western Arctic source for dolomite is the extensive, carbonate-rich Paleozoic terrane in northern Canada (North American Platform; Fig. 1; Okulitch, 1991; Harrison et al., 2008). During the Pleistocene this terrane was repeatedly impacted by the LIS with a subsequent transport of eroded material into the western Arctic Ocean (e.g., Stokes et al., 2005; England et al., 2009). The distribution of dolomite in Arctic sediment cores is thus a robust indicator of the North American provenance and can be used for reconstructing the history of the LIS sedimentary inputs.
Consistent with other cores from the western Arctic Ocean, overall high dolomite content in core ARC4-BN05 has major peaks corresponding to visually identifiable PW/W layers enriched in coarse debris (Fig. 5). As has been suggested in earlier studies (e.g., Stokes et al., 2005; Polyak et al., 2009), we infer that the dolomite peaks are related to pulses of massive iceberg discharge from the LIS during the periods of its destabilization and disintegration. Furthermore, radiogenic isotope studies demonstrate that fine sediment in the dolomitic peaks also has North American provenance (Fagel et al., 2014; Bazhenova et al., 2017). These results indicate that dolomite may have been transported not only by icebergs but also in meltwater plumes coming during deglaciations from the Canadian Archipelago or the Mackenzie River.
As noted above, a change in the stratigraphic pattern of dolomite distribution occurs around unit B10 estimated to correspond to the lower part of MIS 7 (Fig. 6). In older sediments dolomite maxima co-occur with glacial (predominantly gray) intervals, whereas in the younger stratigraphy dolomite peaks in brown sediment or grayish interlayers within brown units (MIS 3, 5, and 7), presumably corresponding to transitional paleoclimatic environments, such as interstadials or stadials within complex interglacial stages.
Other potential mineral indicators related to the North American provenance
are quartz/feldspar and K-feldspar/plagioclase ratios as exemplified by the
BN-05 PCA results (group 4: Fig. 7a), consistent with earlier studies (e.g.,
Vogt, 1997; Zou, 2016; Kobayashi et al., 2016). High Qz
Mineral proxies potentially linked to Siberian provenance make two distinct groups, as reflected in the PCA results (groups 2 and 3: Fig. 7a, Table 3). Group 3 comprises primarily pyroxene, feldspar, and plagioclase, and strongly anticorrelates with the North American proxies, primarily dolomite. The downcore distribution pattern of this group changes from the affinity to interglacials in the lower part of the record to peaks in glacial/deglacial intervals related to MIS 4 and 6 (Fig. 7b). The major source for pyroxene in the Arctic Ocean is the Siberian trap basaltic province that drains to the Kara Sea and western Laptev Sea (Fig. 1; Washner et al., 1999; Schoster et al., 2000; Krylov et al., 2008). However, basaltic rocks related to the Okhotsk–Chukotka province (Fig. 1) may have also provided a significant source of pyroxenes, as exemplified in surface sediments by a relative pyroxene enrichment in the Chukchi Basin on the background of overall low values in the western Arctic Ocean (Dong et al., 2014). Distributions of feldspar and plagioclase at the Siberian margin show elevated contents occurring both in the western Laptev Sea and the East Siberian Sea (Zou, 2016).
Based on a considerable affinity of the pyroxene-feldspar group to brown units and a lack of correlation with coarse sediment fractions, we infer that it is primarily related to sea-ice transport during interglacial/deglacial intervals, with sources potentially including the East Siberian margin and more westerly areas. The difference in both the sources and delivery processes from the LIS proxies may explain an especially strong opposition of these groups. Multiple studies suggest that sea ice from the Kara and Laptev seas may transport sediments to the Canada Basin under favorable atmospheric conditions, such as the positive phase of the Arctic Oscillation (Behrends, 1999; Darby, 2003; Darby et al., 2004, 2012; Yurco et al., 2010), although it remains to be investigated to what extent this circulation pattern could have provided a significant sediment source for the western Arctic Ocean in the Pleistocene.
Another leading sedimentary variable group comprises primarily clay minerals
smectite, kaolinite, and chlorite, and shows affinity to coarse sediment,
especially consistently to fine sand (63–250
We infer that sediment with a concerted enrichment in smectite, kaolinite, and chlorite clay minerals associated with coarse fractions was transported to the Canada Basin primarily in relation to the existence of large ice sheets in northern East Siberia during glacial periods. Radiogenic isotope signature in upper Quaternary records from the Mendeleev Ridge also indicates that the Okhotsk–Chukotka volcanic rocks provided one of the principal end members, especially during MIS 4 and 6 (Fagel et al., 2014; Bazhenova et al., 2017). This sediment had to be transported into the Arctic Ocean directly from the East Siberian/Chukchi margin as the alternative pathway via the Bering Sea only operated at high interglacial sea levels, when the Bering Strait was open for throughflow (e.g., Keigwin et al., 2006; Ortiz et al., 2009). Considering an affinity of the kaolinite-smectite-chlorite group with sediments coarser than clays, corresponding to grain-size modes 2 and 3, their distribution across the basin was likely related to iceberg rafting and glacial underflows, as discussed above in Sect. 5.2.1. A relatively fast and direct delivery mechanism by debris flows and ensuing turbidites may explain a good preservation of fragile clay minerals, normally not resistant to physical erosion.
Some early paleoglaciological studies proposed the existence of a thick Pleistocene ice sheet centered over the East Siberian shelf (Hughes et al., 1977; Grosswald and Hughes, 2002). The inference of former ice sheets/shelves in this region is now corroborated by multibeam bathymetry and sub-bottom data revealing multiple glacigenic features on the top and slopes of the Chukchi and East Siberian margin (Polyak et al., 2001, 2007; Jakobsson et al., 2008, 2014, 2016; Niessen et al., 2013; Dove et al., 2014). ESIS has also been reproduced by numerical paleoclimatic modeling for a large Pleistocene glaciation exemplified by MIS 6 (Colleoni et al., 2016). Sedimentary proxies indicative of the Okhotsk–Chukotka provenance in cores from the Canada Basin may provide an additional tool for reconstructing the ESIS history.
Data points from brown units make up a distinct sedimentary variable group with Mn, foraminiferal numbers, calcite, and fine sediment as lead variables (group 1: Fig. 7a; Table 3). This composition is consistent with the modern-type Arctic Ocean environments characterized by predominant controls of sediment deposition by sea ice, considerable biological activity in summer, and high sea levels. The last of these is important for providing supply of Mn from the surrounding shelves (März et al., 2011; Löwemark et al., 2014). The same condition may also control biological production, and thus foraminiferal numbers, via export of nutrients from the marginal seas (e.g., Xiao et al., 2014), although interaction of this factor with sea-ice conditions still needs to be clarified. We note that the absence (dissolution) of foraminiferal tests in brown units corresponding to MIS 9 and below MIS 11 likely weakens their relationship to other interglacial proxies. Nevertheless, the foraminiferal variable shows a consistent proximity to Mn, clay, and calcite in the PCA results (Fig. 7a).
The mineral having the closest distribution to the main constituents of PC
group 1 is illite, consistent with a predominant occurrence in brown,
interglacial/major interstadial units (Figs. 5 and 7a). Illite is atypical
high-latitude clay mineral, mainly supplied by physical weathering of
metasedimentary and plutonic rocks (Chamley, 1989; Junttila, 2007). High
illite concentrations in surficial Arctic Ocean sediments have been found in
many areas including the Alaska margin and adjacent Canada Basin (Dong et
al., 2014; Kobayashi et al., 2016), East Siberian Sea and the adjacent part
of the Laptev Sea (Wahsner et al., 1999; Kalinenko, 2001; Viscosi-Shirley
et al., 2003; Dethleff, 2005; Zou, 2016), and northern Greenland and
Svalbard regions (Stein et al., 1994). In core ARC4-BN05 illite has
consistently high values in generally fine-grained brown units (Fig. 5),
although peak values may not exactly coincide with those of Mn or
foraminiferal numbers. In addition, illite shows a prominent peak in a very
fine-grained interval at
High contents of calcite in core ARC4-BN05 mostly co-occur with high numbers of foraminifers (Fig. 7a; Table 3), indicating that calcite in these sediments is to a large extent biogenic, consistent with earlier results from the study area (Stein et al., 2010a). Nevertheless, in the lower part of the record, where calcareous fossils are mostly not preserved, calcite shows a considerable affinity to dolomite, which corroborates a mixed, biogenic and detrital nature of calcite in Arctic Ocean sediments (e.g., Vogt, 1997).
The stratigraphically changing pattern of sediment delivery and deposition, including cyclic glacial–interglacial fluctuations and longer-term changes, indicates complex interactions of climatic and oceanographic factors controlling depositional environments in both glacial and interglacial intervals. A long-term trend in interglacial environments is indicated by a shift from predominantly Siberian to more North American provenance, especially strong in MIS 5 and 1, and increasingly high scores of interglacial proxies (group 1), with a threshold around the bottom of MIS 7 (Fig. 7b). Glacial environments show an apparently more complex provenance change, with Siberian sources predominating MIS 4 and 6, and Laurentide provenance controlling MIS 8 and 10 (Fig. 7b). Earlier glaciations, exemplified by a prominent MIS 12 unit, have a mixed signature of high smectite and dolomite contents, likely reflecting a combination of East Siberian and LIS inputs. In addition, an interglacial-type signature (group 1) characterizes some intervals in MIS 4 and 6 as well as intermittent (stadial) intra-MIS 3, 5, and 7 events. We note that MIS 2 is not represented in these data due to its very compressed nature.
The identified changes in sedimentary environments and provenance can be explained by several types of controls, including configuration of ice sheets against sea level and climatic conditions, sediment delivery mechanisms, and circulation. Ice sheet sites and geometry at specific time intervals dictate the timing and location of major sediment discharge events into the Arctic Ocean. Transportation mechanisms, such as by icebergs, debris flows, or suspension plumes, further control sediment delivery to specific sites. Finally, oceanic circulation affects the distribution of sediment across the oceanic basins. This may include surface circulation driving sea ice, icebergs, and surface plumes; deep circulation affecting turbidite/contourite pathways; and downwelling of sediment-laden dense waters.
We infer that sedimentary variations observed in core BN05 and correlative records from the western Arctic Ocean can be explained by the evolution of surrounding ice sheets and associated changes in oceanic conditions, such as circulation, sea ice, and biota. It has been known from early studies (e.g., Clark et al., 1980; Winter et al., 1997) that glacial, notably LIS, impact on the western Arctic Ocean has been steadily increasing over the time span covered by sediment cores from this region. A recent investigation utilizing a more up-to-date stratigraphic paradigm estimated the timing of a step increase in LIS inputs as ca. 0.8 Ma (Polyak et al., 2013), consistent with the onset of major glaciations in the Northern Hemisphere (Head and Gibbard, 2015). Core BN05 provides a record of sediment deposition in the Canada Basin, and thus glacial inputs into the western Arctic Ocean during most of the time interval to follow.
Considering the overall gradual growth of Pleistocene Arctic ice sheets, we infer that the shift from Siberian to North American sources between MIS 12 and 10 was primarily related to the expansion of the LIS, especially the northwestern Keewatin sector that discharges into the western Arctic Ocean. However, its further growth may have had an opposite effect due to a more massive ice sheet that required warmer climatic conditions and/or higher sea levels to destabilize it. Based on data for the last glacial cycle, the Keewatin sector of the LIS rested mostly on relatively elevated terrane of the Canadian Archipelago and adjacent mainland, fringed by a narrow continental shelf and dissected by numerous channels providing conduits for ice streams and evacuation of icebergs at rising sea levels (Stokes et al., 2005, 2009; England et al., 2009; Margold et al., 2015). The latter events are illustrated in BN05 data by intra-MIS 5 stadials with a consistent LIS signature (group 4: Fig. 7b). Especially high LIS scores characterize PW layers 2 and 3 attributed to MIS 5d and late MIS 3, respectively. A similar, LIS-dominated pattern likely represents the last deglaciation as indicated by a number of provenance studies (e.g., Stokes et al., 2005; Jang et al., 2013; Bazhenova et al., 2017).
In comparison to the LIS, a presumably much smaller ESIS, formed on a broad and overall flat East Siberian/Chukchi margin (Niessen et al., 2013; Dove et al., 2014; Colleoni et al., 2016), had to be responsive to sea-level changes even at low levels. It may be possible that the ESIS also increased in size by MIS 6, known as a time of a dramatic increase in glacial inputs from the Barents–Kara Ice Sheet into the eastern Arctic Ocean (e.g., O'Regan et al., 2008). A synchronous MIS 6 expansion of both North American and Siberian ice sheets and related ice shelves might explain the deep-keel glacial erosion of the Lomonosov Ridge at modern water depths exceeding 1000 m (Jakobsson et al., 2016, and references therein).
A concurrent interpretation can be proposed with a focus on sediment transportation processes as deposits of some glacial intervals, notably MIS 12 and parts of MIS 4 and 6, are associated with grain-size mode 2, potentially indicating glacial debris flow/turbidite emplacement. Large debris flows entering the Chukchi Basin and continuing as turbidites into Canada Basin, as exemplified by sub-bottom sonar profiles (Niessen et al., 2013; Dove et al., 2014), may have overprinted deposition from icebergs. We note that deposits of MIS 4 and 6 also contain intervals where Siberian provenance is combined with interglacial positive scores (group 1: Fig. 7b) due to their fine-grained composition along with high illite content. These sediments likely represent deposition from suspension plumes, potentially marking especially strong deglacial meltwater discharge. A prominent fine-grained, finely laminated interval within MIS 4 deglaciation (possibly extending into MIS 3) has been reported from multiple cores across the Chukchi Basin–Mendeleev Ridge area (Adler et al., 2009; Matthiessen et al., 2010; Wang et al., 2013; Bazhenova et al., 2017).
Under modern conditions the BN05 site is mostly controlled by the Beaufort Gyre current circulation system, although it can also be affected by the Transpolar Drift during strong shifts in the Arctic Oscillation (Rigor et al., 2002). This setting probably applies to the Holocene and comparable interglacial conditions (Darby and Bischof, 2004). Some authors have suggested that during glacial periods the surface circulation that controls pathways of iceberg and sea-ice drift may have been considerably different from the modern pattern, with both North American and Siberian sources shortcutting the Arctic Ocean towards the Fram Strait (Bischof and Darby, 1997; Stärz et al., 2012). These changes would have potentially affected the study area, possibly making it more exposed to the Siberian provenance than under present conditions. However, the existing reconstructions based on very limited records with only crude stratigraphic controls, need to be elaborated by spatially and stratigraphically more representative data constraining past circulation changes. In particular, glacial maxima may be elusive, especially in the western Arctic Ocean, due to extremely low sedimentation rates or a hiatus, as exemplified by the LGM (Polyak et al., 2009; Poirier et al., 2012).
Schematic reconstruction of glacial environments in the
western Arctic Ocean and factors controlling sedimentation at the BN05 site
(yellow circle): surface circulation (red and green arrows),
glacioturbidites (orange filled arrow), and relative ice-sheet size (red and
green crosses). See Fig. 1 for modern circulation.
An overall integration of potential controls on sediment deposition in the study area during major identified types of glacial environments are illustrated in Fig. 8. More studies are needed to discriminate between different controls, including proxy records providing higher resolution for target intervals as well as modeling experiments to test spatial and stratigraphic variability in such factors as iceberg and meltwater discharge and their ensuing distribution pathways.
The long-term trend in interglacial environments reflected in a shift from negative to increasingly positive scores of interglacial proxies (group 1: Fig. 7b), with a threshold around the bottom of MIS 7, can be partially explained by the absence of calcareous foraminifers in the lower part of the record. However, even MIS 11, which has abundant foraminifers, has low interglacial scores, suggesting more controls. One possibility is that this trend was related to the evolution of circum-Arctic ice sheets that would have inevitably incurred changes in oceanic conditions, such as circulation and sea ice. An expansion of perennial sea ice in the western Arctic Ocean near the MIS 7 bottom has been proposed based on foraminiferal assemblages (Polyak et al., 2013; Lazar and Polyak, 2016). This step change has been tentatively attributed to the LIS growth that may have affected sea-ice conditions via increased albedo and/or higher meltwater inputs. This inference is consistent with a coeval change from mostly Siberian (group 3) to North American (group 4) provenance during interglacials in BN05 (Fig. 7b). In addition to a more lingering LIS during interstadials/interglacials, this shift in provenance could be related to a strengthening of the Beaufort Gyre as more sea ice filled the western Arctic Ocean.
More limited sea-ice cover in the older part of the middle Pleistocene could have also enhanced the production of dense brines at the Siberian margin, resulting in a deeper convection and cascading of shelf sediments to the deep basin. This scenario would explain an unusual grain-size composition of sediments in the older interglacials combining mode 2, indicative of winnowed silt, with a typical interglacial fine-grained mode 1.
Sediment core ARC4-BN05 was collected from the Canada Basin in the vicinity
of the Chukchi Plateau and the Mendeleev Ridge, Arctic Ocean, on the fourth
Chinese National Arctic Research Expedition (CHINARE-IV). Based on
correlation to earlier proposed Arctic Ocean stratigraphies (e.g., Adler et
al., 2009; Stein et al., 2010a; Polyak et al., 2013) and AMS
Glacially derived sediment can be discriminated between the North American and Siberian provenance by their mineralogical and textural signature. In particular, peaks of dolomite debris, including large dropstones, track the Laurentide Ice Sheet (LIS) discharge events, while the East Siberian Ice Sheet (ESIS) inputs are inferred from combined peaks of smectite, kaolinite, and chlorite associated with coarse sediment. Siberian provenance is also identified from high content of pyroxene, feldspar, and plagioclase, unrelated to coarse sediment. This sedimentary signature is interpreted to indicate sea-ice transport from the Siberian margin during interglacial/deglacial intervals. Full interglacial environments are characterized by overall fine grain size, high content of Mn (and resulting dark-brown sediment color), and elevated contents of calcite and chlorite. Foraminiferal tests are abundant in interglacial units in the upper part of the record (MIS 1–7) and estimated MIS 11 but have very low numbers in other interglacials older than MIS 7, apparently due to dissolution.
In addition to glacial–interglacial cyclicity, the investigated record indicates variable impacts of LIS vs. ESIS on sediment inputs at different glacial events, along with a long-term change in middle to late Quaternary sedimentary environments. Based on the age model employed, major LIS inputs to the study area occurred during MIS 3, intra-MIS 5 and 7 events, MIS 8, and MIS 10, while ESIS signature is characteristic for MIS 4, MIS 6 and MIS 12. These differences may be related to ice-sheet configurations at different sea levels, sediment delivery mechanisms (iceberg rafting, suspension plumes, and debris flows), and surface circulation. A long-term shift in the pattern of sediment inputs shows an apparent step change near the estimated MIS 7/8 boundary (ca. 0.25 Ma), consistent with more sea-ice growth in the Arctic Ocean inferred from benthic foraminiferal assemblages (Lazar and Polyak, 2016). This development of Arctic Ocean paleoenvironments possibly indicates an overall glacial expansion at the western Arctic margins, especially in North America. Such expansion may have affected not only glacial but also interglacial conditions via increased albedo and/or higher meltwater inputs, as well as a strengthening of the Beaufort Gyre circulation as more sea ice filled the western Arctic Ocean.
No data sets were used in this article.
The authors declare that they have no conflict of interest.
We are grateful to the team of the 4th Chinese Arctic Research
Expedition for their assistance with sample collection. Special thanks to
Shijuan Yan for help with sampling and to Quanshu Yan for help in paper
editing. This work was jointly supported by the Research Foundation of the
First Institute of Oceanography, State Oceanic Administration of China
(no. 2013G07, 2014G30), the Chinese Polar Environment Comprehensive
Investigation & Assessment Programmes (no. CHINARE 2017-03-02), and the National Natural
Science Foundation of China (no. 41306205, 41676053, 40176136). Leonid Polyak's
participation was supported by the US National Science Foundation award
ARC