The Pan-Arctic Continental Slope: Sharp Gradients of Physical Processes Affect Pelagic and Benthic Ecosystems

Continental slopes – steep regions between the shelf break and abyssal ocean – play key roles in the climatology and ecology of the Arctic Ocean. Here, through review and synthesis, we find that the narrow slope regions contribute to ecosystem functioning disproportionately to the size of the habitat area (∼6% of total Arctic Ocean area). Driven by inflows of sub-Arctic waters and steered by topography, boundary currents transport boreal properties and particle loads from the Atlantic and Pacific Oceans along-slope, thus creating both along and cross-slope connectivity gradients in water mass properties and biomass. Drainage of dense, saline shelf water and material within these, and contributions of river and meltwater also shape the characteristics of the slope domain. These and other properties led us to distinguish upper and lower slope domains; the upper slope (shelf break to ∼800 m) is characterized by stronger currents, warmer sub-surface temperatures, and higher biomass across several trophic levels (especially near inflow areas). In contrast, the lower slope has slower-moving currents, is cooler, and exhibits lower vertical carbon flux and biomass. Distinct zonation of zooplankton, benthic and fish communities result from these differences. Slopes display varying levels of system connectivity: (1) along-slope through property and material transport in boundary currents, (2) cross-slope through upwelling of warm and nutrient rich water and down-welling of dense water and organic rich matter, and (3) vertically through shear and mixing. Slope dynamics also generate separating functions through (1) along-slope and across-slope fronts concentrating biological activity, and (2) vertical gradients in the water column and at the seafloor that maintain distinct physical structure and community turnover. At the upper slope, climatic change is manifested in sea-ice retreat, increased heat and mass transport by sub-Arctic inflows, surface warming, and altered vertical stratification, while the lower slope has yet to display evidence of change. Model projections suggest that ongoing physical changes will enhance primary production at the upper slope, with suspected enhancing effects for consumers. We recommend Pan-Arctic monitoring efforts of slopes given that many signals of climate change appear there first and are then transmitted along the slope domain.

Continental slopes -steep regions between the shelf break and abyssal ocean -play key roles in the climatology and ecology of the Arctic Ocean. Here, through review and synthesis, we find that the narrow slope regions contribute to ecosystem functioning disproportionately to the size of the habitat area (∼6% of total Arctic Ocean area). Driven by inflows of sub-Arctic waters and steered by topography, boundary currents transport boreal properties and particle loads from the Atlantic and Pacific Oceans along-slope, thus creating both along and cross-slope connectivity gradients in water mass properties and biomass. Drainage of dense, saline shelf water and material within these, and contributions of river and meltwater also shape the characteristics of the slope domain. These and other properties led us to distinguish upper and lower slope domains; the upper slope (shelf break to ∼800 m) is characterized by stronger currents, warmer sub-surface temperatures, and higher biomass across several trophic levels (especially near inflow areas). In contrast, the lower slope has slower-moving currents, is cooler, and exhibits lower vertical carbon flux and biomass. Distinct zonation of zooplankton, benthic and fish communities result from these differences. Slopes display varying levels of system connectivity: (1) along-slope through property and material transport in boundary currents, (2) cross-slope through upwelling of warm and nutrient rich water and down-welling of dense water and organic rich matter, and (3) vertically through shear and mixing. Slope dynamics also generate separating functions through (1) along-slope and across-slope fronts concentrating biological activity, and (2) vertical gradients in the water column and at the seafloor that maintain distinct physical structure and community turnover. At the upper slope, climatic change is manifested in sea-ice retreat, increased heat and mass transport by sub-Arctic inflows, surface warming, and altered vertical stratification, while the lower slope has yet to display evidence of

INTRODUCTION: MOTIVATION AND DEFINITIONS
Continental slopes are a ubiquitous tectonic feature of the global ocean, long recognized for their central climatological and ecological roles at the interface of shelves and basins (Springer et al., 1996;Colloca et al., 2004;Azzellino et al., 2008;Bertram et al., 2017). In general, slope system dynamics are associated with strong vertical and cross-slope gradients that contrast with alongslope bands of relatively more uniform conditions (Cacchione et al., 2002). This situation strongly applies to the Arctic Ocean, as was first noted by oceanographer and geographer Gakkel (1957). Using geomorphological and visual observations from aerial reconnaissance, he identified the Arctic circumpolar continental slope and emphasized its dynamic ice regime, along-slope circulation, and enhanced winter heat loss (Gakkel, 1957). Also recognizing the distinct role of slopes, zoologist Uspenskiy (1973) introduced the term the Arctic ring of life, referring roughly to the continental slope region and highlighting its biological richness compared to adjacent shelf or basin areas, specifically noting more frequent observations of higher trophic level predators including polar bears and narwhals. Given this role, Arctic slopes can be considered as their own pan-Arctic contiguous domain (i.e., functional unit, Carmack and Wassmann, 2006) when considering holistic functioning of the Arctic Ocean. This requires, however, an integration of earlier findings into a pan-Arctic perspective, which currently is lacking.
The continental slope is that region starting seaward of the continental shelf break which globally is often marked by the 200 m isobath (but see modification below) and a strong sloping angle (typically > 4 • ) of the seafloor. The continental slope Abbreviations: ACBC, Arctic Circumpolar Boundary Current; AW, Atlantic Water; PW, Pacific Water; PSW, Pacific Summer Water; PWW, Pacific Winter Water. extends to the continental rise of the ocean floor which often is at ∼2000-2500 m or where the angle becomes <∼1 • 25 (Hay, 2016). Combined, the shelf, slope and rise are commonly referred to as the continental margin. As in other regions of the globe, the slopes of the Arctic Ocean are intersected by numerous canyons, troughs, ridges and straits, resulting in complex morphological, oceanographic and biological structures (Jakobsson et al., 2012). Arctic Ocean here refers to the ocean area bounded by Bering Strait on the Pacific side, by landmasses in the Arctic Ocean interior, and by Fram Strait and the western Barents Sea shelf break on the Atlantic side ( Figure 1A). The geological history of the Arctic Ocean has resulted in shelf-break depths that lie at as little as ∼60 m in places off the Siberian shelves to ∼400 m off the Barents Sea (Jakobsson, 2002; Figure 1B). The upper boundary of the continental rise in the Arctic ranges from ∼2000 to 3000 m depending on location (Jakobsson, 2002). For the purpose of this paper we generally focus on the slope depth range of ∼200-2500 m (starting at ∼400 m in the deeper Barents Sea). This area constitutes roughly 6% of the total Arctic Ocean area.
Property gradients in the narrow band above the slope region are expressed in three directions that we use as guiding structure throughout the present paper: (1) the along-slope (azimuthal) direction around the entire Arctic perimeter, and (2) the vertical direction (downwards) from ocean surface to the seafloor, and (3) the cross-slope (radial in the Arctic) direction from the continental shelf break down-slope ( Figure 1C). Beside the three directions, we further distinguish two vertical zones: the upper slope (approximately 200-800 m on the Pacific side and approximately 400-800 on the Atlantic side), and the lower slope (approximately 800-2500 m) ( Figure 1A). The basis for this designation is the dominance of sub-Arctic inflows from the north Atlantic and north Pacific that flow around the basin perimeter above the upper slope (details in section "Physical Oceanography of Arctic Slopes;" Carmack and Wassmann, 2006). The lower boundary of the upper slope is associated with the bottom of the Atlantic layer, defined conventionally -yet somewhat arbitrarily -by the 0 • C isobath. The lower slope then begins with the transition to the Arctic Deep Water (Aksenov et al., 2011;Pnyushkov et al., 2018), extends to the approximate transition to the continental rise, and is associated with greatly declining and less advected biomass Vedenin et al., 2018).
Once we have described bio-physical features along the above gradients and zones, we integrate physical (section "Physical Oceanography of Arctic Slopes") and biological (section "Gradients in Biological Communities at Arctic Slopes") information to conceptualize the ecological functions that Arctic  Figure 2 (with the exception of transect 9) and have been centered so that X = 0 is on the 1000 m isobaths for consistency. Sections are trimmed to show -150 to +250 km from the 1000 m isobath. Gray shows the whole profile while brown is the section that overlaps with temperature transects in Figure 2; black is land. Black horizontal lines mark the slope range focused on in the present paper; gray line marks the lower boundary of the Atlantic Layer (designated as the lower boundary of the "upper slope domain"). (C) The three dimensions of Arctic continental slopes considered in this paper: (1) along-slope (azimuthal) along the basin perimeter, (2) vertical (ocean surface to seafloor), (3) cross-slope (from the shelf break off-shore down slope).
continental slopes add to the Arctic Ocean system in section "Functions of the Slope: System Connectivity vs. Separation." Some of these functions result from the geomorphology of continental slopes in general (Nash et al., 2004;Levin and Dayton, 2009), and others from specific regional hydrographic settings in the Arctic in particular. These functions can act through both connecting and separating processes between adjacent shelves and basins. Finally, in section "Towards a Future Arctic Slope System, " we summarize ongoing climatic changes overlying the Arctic slopes; specifically warming (Polyakov et al., 2020a, this issue) and sea ice loss (Kwok et al., 2013), and then briefly present a modeling exercise to investigate potential associated biological changes.

Along-Slope: Boundary Currents
The Arctic continental slopes provide a "handrail" for incoming sub-Arctic waters, meaning the slope topography steers these waters counterclockwise (cyclonically) around the basin perimeter. The most prominent pan-Arctic circulation feature transporting inflowing Atlantic Water (AW) is the Arctic Circumpolar Boundary Current (ACBC), which funnels Atlantic-origin water as a narrow, contiguous stream that encircles all Arctic Basins (Aagaard, 1989;Rudels et al., 1994;Aksenov et al., 2011). Similar dynamics hold for incoming waters of Pacific origin, but these are more complicated and largely retained within the Amerasian Basin (McLaughlin et al., 1996;Shimada et al., 2006;Pickart et al., 2009Pickart et al., , 2013aAksenov et al., 2016). Such boundary currents are common to continental slopes globally (Huthnance, 1981) but are especially strong in the Arctic owing to the joint effects of the increase in Coriolis force with latitude and the exposure to ice-cover. The ACBC initially forms from the entry of warm and saline AW through Fram Strait (as the Fram Strait Branch) and then flows along the northern Barents Sea slope. A second branch of AW crosses the Barents Sea (the Barents Sea Branch), then enters the Arctic Ocean, mainly through St. Anna Trough, where the two branches remerge in the northern Kara Sea slope (Hanzlick and Aagaard, 1980). These are then joined by the drainage of Siberian Shelf waters before continuing around the Arctic Ocean perimeter (Aksenov et al., 2011). This circulation pattern extends through the full depth of the pan-Arctic continental slope. Current velocities, however, are strongest at the upper slope, especially along the Barents Sea slope (peak flow rates of >20 cm/s), and noticeably weaker in deeper water and along the basin perimeter on the Amerasian side (Aksenov et al., 2011;Pnyushkov et al., 2015;Menze et al., 2019). In the surface layer, the subsurface ACBC can be opposed by wind-driven circulation, such as the clockwise-flowing Beaufort Gyre (Proshutinsky et al., 2009).
To demonstrate the progressive evolution of AW properties while flowing counterclockwise along the upper pan-Arctic slope we here use a near-synoptic view of the Arctic shelfto-slope structure in summer, 2015, in temperature sections (Figure 2), and temperature/salinity plots (Figures 3A,B), and vertical salinity profiles ( Figure 3C). Initially, AW outcrops at the surface in Fram Strait and north of Svalbard; typically, in winter Renner et al., 2018). When sea ice melts, either due to contact with the inflowing warm AW, or due to solar input during spring and summer, a cold and fresh (and less dense) surface layer develops (Untersteiner, 1988;Rudels et al., 2013Rudels et al., , 2014. When this happens in summer, there is substantial solar heating of the shallow, seasonal surface layer, making it even warmer than the AW below (e.g., north of the Svalbard slope) [Figures 2(1), 3A(1); note that an expanded temperature and salinity range would reveal warmer and fresher varieties of surface water north of the Svalbard slope]. Slightly further east the core of the AW subducts and becomes capped by the cold and fresh near-surface layer [Figures 2(2), 3A(2)], and subsequently cools and deepens during its translation along the slopes of the Barents and Kara Seas [Figures 2(3), 3A(3),B]. Just upstream of 95 • E at St. Anna Trough, the Barents Sea-modified branch of AW debouches the continental slope region and interacts with the Fram Strait Branch. As a result, the AW core is further cooled, with maximum temperatures decreasing by at least 1 • C, while also gradually freshening [Figures 2(4,5), 3A(4,5)]. Further east in the Laptev Sea in particular, cold and saline waters created in strong polynyas that form in winter and spring (Bareiss and Görgen, 2005) cascade down the continental slope and contribute to the halocline complex (Aagaard et al., 1981;Martin and Cavalieri, 1989;Ivanov and Golovin, 2007;Walsh et al., 2007). Numerous canyons and passages in the Severnaya Zemlya region (100-110 • E) (Shokalsky Strait, Vilkitsky Strait) guide the export of these cold and dense waters toward the slopes (Janout et al., , 2017. The injection of these waters then results in further cooling of the AW as it continues eastward along the Siberian shelves [Figures 2(6), 3A(6)]. At the same time, large rivers dilute the Arctic Ocean surface and strengthen the fresh cap over the AW as it propagates eastward along-slope. Note that prior to leaving the Nansen Basin the salinity gradient below the core (i.e., T max ) of AW is negative, thus permitting salt-fingering, while beyond 135 • E it is positive. Upon reaching the Chukchi Sea the joint effects of the Pacific Water (PW) inflow and the clockwise wind field over the Beaufort Sea result in a further deepening of the AW core [Figures 2(7,8), 3A(7)]. Beyond this point the AW begins its exit into the Greenland Sea via the western Fram Strait and the East Greenland Current (Rudels et al., 2012;Håvik et al., 2017). During its propagation along the Arctic slopes, the AW thus decreases from a >3 • C warm, several hundred-meterthick layer near the inflow region to a thinning, <1 • C-layer that is centered around a depth of 400 m along the slope of the Canada Basin (Figures 4A,B). Also, along this transit the salinity of the core layer, referenced to the surface, freshens substantially from about 35.0 to 34.8 (Figure 3B), the density decreases slightly from about 27.94 to 27.92, while stratification (i.e., Brunt-Vaisala buoyancy frequency) of the overlying water column increases by a factor of 3-10 (cf. Polyakov et al., 2018).
Boundary currents also form along the Pacific Arctic slopes because of inflowing PW through Bering Strait, but with different flow patterns than those of the AW ( Figure 4B). Incoming PW is strongly modified by seasonal processes while crossing the broad Bering and Chukchi Sea shelves, so that it arrives at the shelf slope boundary as either Pacific Summer Water (PSW; relatively warm and fresh) or Pacific Winter Water (PWW; colder Frontiers in Marine Science | www.frontiersin.org FIGURE 2 | Pan-Arctic potential temperature transects of the Arctic shelf-to-slope structure in summer of 2015 showing Atlantic Water (AW) position over the slopes; transects are the same as in Figure 1. Warm AW near the surface in Fram Strait / north of Svalbard (1, 2) becomes capped by a colder and fresher surface, and gradually cools and deepens along the slopes of the Barents,Kara,and Laptev Seas (3,4,5). Canyons and passages bring cold and dense shelf water to the slopes, further cooling the AW on the E Siberian shelf (6). Pacific Water inflow and clockwise wind fields over the Beaufort Sea result in further deepening of the AW core (7, 8) before beginning to move towards exiting via western Fram Strait and the East Greenland current. We lack oceanographic data from transect (9) in Figure 1. and more saline) (Weingartner et al., 2005;Pickart et al., 2016). Much of PSW exits via Barrow Canyon, while smaller fractions exit through Herald Canyon and Central Channel (Weingartner et al., 2017). Some fraction turns eastward as a subsurface flow at ∼40-80 m along the shelf-slope boundary, referred to either as the Beaufort Undercurrent (Aagaard, 1984) or shelf-break jet (Pickart et al., 2005; Figure 4B). This thermohaline feature differs from the highly transient, wind-driven jet that forms at the Pacific Arctic slopes during upwelling and downwelling events . The remaining fraction of PSW exiting  (A) Temperature/Salinity diagrams for selected slope transects 1-8 shown in Figure 2 for temperature range -2-4 • C and salinity range 32.0-35.2. Individual stations are color-coded by longitude and refer to transect locations in Figure 2. Numbers in circles (from 1 to 11), in contrast, refer to the following noteworthy features: (1) crossing the Fram Strait Branch (FSB) northwest of Svalbard prior to encountering ice and Arctic surface water, (2) flowing north of Svalbard subsequent to mixing with and subduction of the core of AW below Artic surface water, (3) the core of AW as it continues to cool and freshen along the Barents Sea slope, (4) water from the Barents Sea Branch (BSB) entering via the St. Anna Trough and converging with the FSB, (5) the warm core of the FSB is cooled, freshened and displaced offshore by the BSB, (6) continued cooling and freshening of the AW core along the Siberian shelves, (7) after crossing the Lomonosov Ridge into the Amerasian Basin, (8) the Atlantic-Pacific halocline front where the gray dashed line indicates the mixing line leading to bottom water in the Nansen basin, (9) the T/S correlation curve for modified AW in the BSB prior upon entering the slope domain through St. Anna Trough, and 10) the maximum salinity reached by Cold Halocline Water. (B) T/S for salinity range 34.2-35.2. Orange: T/S from 95 • E. The red dashed line indicates a straight line tangent to σ 0 = 27.94 kg m −3 ; the point at which this line crosses the freezing temperature yields the maximum salinity that surface waters can reach without triggering the caballing instability (34.53) (cf. Fofonof, 1956). (C) Vertical profiles of salinity for slope stations at transects 1-8 from Figure 2, showing progressive freshening. Stars show the depth of the AW core (T max ) at each station location; note the progressive deepening of the AW core layer as the boundary current progresses around the Arctic basin.
through Barrow Canyon is drawn offshore into the Beaufort Gyre  or ejected off-slope by mesoscale eddies (cf. D 'Asaro, 1988;Mathis et al., 2007). In addition to PW joining the shelf-break jet, recent studies indicate that a portion of the PWW turns westward over the slope, forming a boundary current that appears to follow the clockwise circulation pattern of the Beaufort Gyre (Corlett and Pickard, 2017;Spall et al., 2018;Li et al., 2019). Long-term observations in the Canada Basin showed an increase in PWW volume by 18% from 2002 to 2016 (Zhong et al., 2019), thus underlining the increasing role of PWW as a sub-surface freshwater source to the Beaufort Gyre.
In contrast to the dynamic shelf break and upper slope described thus far, the lower slope lacks large gradients and is comparatively quiescent. Below the core of AW and to the depths of the Lomonosov and Alpha-Mendeleev ridges, potential temperature decreases and salinity increases slightly to the seabed along the lower slope (Timmermans et al., 2003;Björk and Winsor, 2006). The lower slope has current velocities that are generally below 2 cm/s and is dominated by smaller-scale processes (Bluhm et al., 2015 and references therein). Visual evidences of abundant and persistent animal tracks despite low faunal densities support the notion of extremely low flow rates (Zhulay et al., 2019). Overall, little is actually known about the lower slope environment and more work is required to characterize it.

Vertical and Cross-Slope Structure and Processes
The water column above the slopes is generally salt-stratified, with relatively fresh surface water and a halocline above the Atlantic Layer ( Figure 3C). Specifically, the combination of slope, and lower circle denoting eastward directed AW. The Eurasian view is from the west to the east; the Amerasian view is from the east to the west. Dashed lines denote velocities, while solid lines depict temperature; depth is shown as a square root scale. (C) Upwelling circulation across a shelf with a shallow shelf break, typical of interior shelves. Here the wind is directed into the panel such that the offshore Ekman transport (pink arrow) is compensated by an onshore flow across the shelf/slope boundary (blue arrow); in this case drawing cooler and presumably nutrient-rich waters onto the shelf. Also shown is a transient shelf/slope break jet, or bottom-boundary layer, also directed into the panel. (D) Brine release during sea ice formation and brine-driven circulation across a shelf and down slope. Q is heat flux (Q) and stars indicate frazil ice formation. Inset profiles are shown for density (ρ) and velocity (v). Circles with dots denote flow out of the panel while circles with crosses denote flow into the panel.
freshwater inputs (seasonal sea ice melt, river water, glacial melt water), surface warming and wind-driven mixing result in a seasonal mixed layer (approximately 40-60 m thick). Underneath it, the halocline complex with (Amerasian slopes) or without PW (Eurasian slopes) (McLaughlin et al., 1996;Rudels et al., 2004) forms a boundary that inhibits upward mixing of nutrients , but see section "Separation Generates and Maintains Structure." Multiple steps in the halocline reflect the different sources that supply water to the halocline, e.g., sea ice melt, river discharge, net precipitation, Pacific inflows and Siberian shelf drainage (Brown et al., 2020). The Atlantic Layer below (to ∼800 m) is warmer and saltier than the underlying Arctic Ocean Deep Water at the lower slope (Aagaard et al., 1985;Rudels et al., 2012). Variations in the vertical stratification resulting from this layering are primarily a result of the interplay of the warm and salty AW, sea ice melt, incoming rivers and the PW inflow, leading to distinct regimes in vertical stratification along the pan-Arctic slope. These can be tied directly to the shelf typology as discussed by Carmack and Wassmann (2006) and Bluhm et al. (2015): At the slope of the Atlantic inflow shelf region north of Svalbard/Barents Sea where AW seasonally extends fully to the surface, a halocline is initially lacking and is first established by the interaction of AW with sea ice melt water (Figure 3C). At the slopes of the interior shelves subsequent to the St. Anna Trough Siberian river water reaches the slopes, and brine-enriched Siberian shelf waters flowing offshelf form additional halocline waters. At the slopes of the Pacific inflow shelf (the Chukchi Sea) and into the Pacific Interior shelf (Beaufort Sea) the combined effects of the PW inflow and the Beaufort Gyre circulation greatly increase vertical stratification. And finally, at the slopes of the outflow shelves (the Canadian Arctic Archipelago and the northern Greenland slopes) vertical stratification is thought to remain largely unchanged until waters exit along East Greenland. Increasing glacial runoff to the East Greenland shelf, however, may in the future not only increase the local stratification and is likely to speed up the coastal circulation for instance in the wind-and buoyancy-forced East Greenland Coastal Current (Sutherland and Pickart, 2008), but also already has the potential to reach the outer shelf and boundary currents (Hendry et al., 2019). Overall, the outflow shelves are strongly impacted by liquid and solid freshwater export from the Arctic Ocean (Michel et al., 2015), and will likely continue to receive additional freshwater loads from glacial runoff with not yet understood impacts on the slope system.
The above described progressive cooling and freshening in the along-slope direction from Eurasian to Amerasian slopes involves numerous vertical and cross-slope displacements of water mass layering and frontal zone structures. For example, vertical displacements, forced by winds, tides, eddies and seasonal heating and cooling all give rise to significant (>1 • C) temperature fluctuations at any given location on the upper slope (Baumann et al., 2018). Properties of the upper slope domain are especially impacted by cross-slope mechanisms that act to break the dynamic barrier between shelves and the Arctic Ocean's interior, such as upwelling and brine drainage (see section "Connectivity") ( Figures 4A,B). Cross-slope exchange is especially amplified at sudden topographic junctions (straits, canyons, relaxations in the slope angle, etc.). One example is where Barrow Canyon intersects the Beaufort/Chukchi Sea continental slope at nearly a right angle, and where on average 25-50% of northward flowing PW exits the Chukchi shelf (Itoh et al., 2013;Weingartner et al., 2017). Downcanyon flow can only partially adjust to the sharply turning isobaths when the canyon abuts the slope, thus leading to an overshoot and injection of shelf waters onto the upper slope or into the basin. Consequently, regions of cross-slope flow often are associated with the formation of energetic eddies that carry mass, heat, and biogeochemical constituents from the shelf regions to the slope and the basin interior Mathis et al., 2007;Hattermann et al., 2016;Våge et al., 2016;Pnyushkov et al., 2018). Another example occurs in the Mackenzie Trough off the Canadian Beaufort shelf, where upwelling-favorable winds draw Atlantic waters as much as 400 m or more up canyon, which subsequently collapse back into the basin, generating eastward propagating internal Kelvin waves along the slope (Carmack and Kulikov, 1998). Slope angle also appears to govern contrasts in shelf-basin exchange; for example, the basin seaward of the steep Laptev Sea slope shows comparatively little signal of Laptev Sea shelf waters egressing seaward, suggesting that shelf-basin exchange here is constrained. In contrast, further east along the wide East Siberian slope, enhanced transport of shelf waters into the basin is found (Anderson et al., 2017). These exchange windows are particularly relevant near the large rivers, as freshwater is an important contribution for the water column stratification of the Arctic Ocean (Carmack et al., 2016), as discussed further in section "Functions of the Slope: System Connectivity vs. Separation."

GRADIENTS IN BIOLOGICAL COMMUNITIES AT ARCTIC SLOPES
Gradients in standing stock, community composition, and production are prominent at slope areas, both vertically in the water column and across-slope for both pelagic and seafloor communities. Spatial patterns in primary production levels reflect oceanographic patterns and processes described in section "Physical Oceanography of Arctic Slopes, " and contribute to driving spatial patterns in standing stocks across trophic levels [of which we here consider zooplankton and benthos (section "Gradients in Primary Production and Lower Trophic Level Biomass"), and fishes, seabirds and marine mammals (section "Higher Trophic Level Biomass and Distribution") and their community composition (section "Community Structure and Biodiversity Trends")]. Geographical variation in each trophic level is discussed beginning at the Atlantic gateway slope and proceeding eastward around the basin perimeter, as data coverage allows.

Gradients in Primary Production and Lower Trophic Level Biomass
The hydrography over the slopes described above, combined with the steep depth gradient, provide the backdrop for biological gradients along and across the Arctic slopes. In the alongslope direction, levels of primary production vary as a result of gradients in nutrient concentrations, supply and sources. The AW inflow is the primary nutrient source for the Eurasian Arctic slope  and these advective inputs fuel gross primary production levels above the slope northwest of Svalbard of 0.1-0.9 mg C m −2 d −1 (May and August, Svensen et al., 2019), with annual estimates of ∼70-100 g C m −2 year −1 for that area (Slagstad et al., 2011;Frey et al., 2018). In the slope areas of the Siberian and North-American interior shelves, the large rivers add only limited amounts of nutrients (specifically nitrogen and silicic acid, Codispoti et al., 2013), but introduce substantial loads of CDOM and inorganic sediments which can degrade the light regime and limit primary production (Popova et al., 2012). As a result of this and increasing distance from sub-Arctic inputs along the basin perimeter, primary production levels at the slope in these interior shelf regions are less than ∼30 g C m −2 year −1 (Slagstad et al., 2011;Matrai et al., 2013). By comparison, near the Pacific inflow, nutrient-rich PW enhances production to levels of ∼170 g C m −2 year −1 on the Chukchi slope (Hill et al., 2018). Then, production drops off sharply in the Beaufort Sea (∼90 g C m −2 year −1 , Hill et al., 2018) and over the slopes towards the Amerasian Basin (∼15-30 g C m −2 year −1 , Carmack et al., 2004;Frey et al., 2018). Primary production levels along the slopes of the outflow shelves are comparatively poorly quantified, but tend to be lower than in inflow shelves (∼60 g C m −2 year −1 and dropping off beyond the slopes; Matrai et al., 2013;Michel et al., 2015;Frey et al., 2018). We note, however, that changes in primary productivity are ongoing and predicted (see section "Towards a Future Arctic Slope System") and unusual peaks in algal biomass were for example recently observed along the Greenland Sea slope (Frey et al., 2019).
Enhanced abundance and biomass of zooplankton, benthos, and upper trophic level consumer communities are apparent over the inflow slopes and to some degree beyond. For zooplankton, this enhancement is evident in cross-slope transects of integrated mesozooplankton biomass in both the Eurasian and Amerasian sectors of the slope ( Figure 5A). Vertically integrated biomass levels from surface to near-bottom increase from 1 to 7 g dry weigh (DW) m −2 near the shelf-break to 5-15 g DW m −2 over the slope ( Figure 5A). Vertically, zooplankton concentrations are highest in the 50-200 m layer (Kosobokova and Hirche, 2009), suggesting the higher values over the slope are not merely a bias of integration depth. In the along-slope direction biomass is highest near sub-Arctic inflows, in particular close to the core of the Atlantic inflow (Kosobokova and Hirche, 2009;Kosobokova, 2012;Basedow et al., 2018). This biomass dense area (peak measured at 24 g DW m −2 NE of Severnaya Zemlya; Kosobokova and Hirche, 2009) is dominated by the expatriate copepod Calanus finmarchicus which, after its injection north of Svalbard, is transported alongslope far into the slope areas of the interior Siberian shelves (Kosobokova and Hirche, 2009;Kosobokova, 2012;Wassmann et al., 2015). Similarly, zooplankton biomass in surface layers in the western part of the Beaufort Sea slope is enriched by Pacific zooplankton inputs at the Pacific inflow (Berline et al., 2008;Smoot and Hopcroft, 2017). Here, copepods in the genus Neocalanus in particular contribute advected biomass that is highest in the upper portion (200-500 m) of the Atlantic Layer (Smoot and Hopcroft, 2017) where biomass levels are comparable to the slope in the Atlantic inflow ( Figure 5A). Within the western Beaufort Sea, biomass of these Pacific species declines moving eastward.
Benthic and demersal fish biomass generally declines with depth across the slopes towards abyssal plains. This trend applies to all size fractions including macroinfauna (≥0.5 or 1 mm, inside sediment; Figure 5B) (Wlodarska-Kowalczuk et al., 2004;Bluhm et al., 2011;Grebmeier, 2012;Grebmeier et al., 2015;Vedenin et al., 2018), mega-epifauna (benthos in trawls or on photographs) (Ravelo et al., 2020, Jørgensen et al., unpubl. data), and demersal fish (Majewski et al., 2017;Norcross et al., 2017), with a less pronounced decrease for the smaller meiofauna (≥32 µm -0.5 or 1 mm; Vanreusel et al., 2000;Wei et al., 2010). The underlying reason for benthic biomass declines with depth is primarily the diminishing vertical flux of particulate organic matter (i.e., food particles) with increasing depth (Wiedmann et al., 2020, and references therein). Macrofaunal biomass drops off from peaks of 10-20 g C m −2 on inflow shelves (Grebmeier, 2012), and <1 g C m −2 on the interior Laptev Sea shelf (Vedenin et al., 2018), to an average ∼0.5 g C m −2 at the upper slope and <0.2 g C m −2 at the lower slope (Wlodarska-Kowalczuk et al., 2004;Bluhm et al., 2005Bluhm et al., , 2011Grebmeier et al., 2006;Nelson et al., 2014;Vedenin et al., 2018; estimated from replicate van Veen grab or box core samples at discrete sampling depths). Gross estimates of wet weight (WW) biomass of epibenthic and demersal fish at the Amerasian slope decrease from 2 to 20 kg WW 1000 m −2 at the outer shelf to generally <5 kg WW 1000 m −2 at 1000 m, with fish biomass contributing generally <15% to total biomass (Norcross et al., 2017;Ravelo et al., 2020; estimated from trawl hauls at discrete depths). In slope areas near advective inflows from the Atlantic and Pacific, however, the down-slope profile in biomass exhibits a peak at the depth of the shelf break where benthic biomass can be enhanced ( Figure 5B). In the PW inflow specifically, macrobenthic biomass is clearly elevated at the head of Barrow Canyon (Grebmeier, 2012), and a tongue of elevated biomass tracks the outflowing nutrient and particle-rich PW along the upper slope of the western Beaufort Sea to roughly 150 • W (Bilyard and Carey, 1979;Dunton et al., 2005;Ravelo et al., 2015). The active fishery on the deep-water shrimp Pandalus borealis at the upper slope north of Svalbard (Misund et al., 2016;Haug et al., 2017) gives testimony of enhanced biomass at the Atlantic inflow slope area as well, although exactly comparable numbers are yet unpublished.

Higher Trophic Level Biomass and Distribution
For the most abundant Arctic fish, Boreogadus saida (Polar cod in European and Arctic cod in North American nomenclature), biomass of certain age classes also appears to be elevated along the upper Arctic slopes. This pattern is, however, so far only confirmed for the Amerasian slopes (Crawford et al., 2012;Geoffroy et al., 2015;Majewski et al., 2016;Logerwell et al., 2018) where biomass levels are higher in the western Beaufort than farther east; it remains unclear if this "cod belt" is a pan-Arctic phenomenon. In the along-slope direction, distribution patterns of (pelagic) B. saida larvae and juveniles were highest over slope areas in the central Kara, the Laptev and East Siberian seas as well as seasonally in the eastern Beaufort Sea (Suzuki et al., 2015;Mishin et al., 2018), while they were more abundant at the coast in the Barents and southwestern Kara Seas. Peaks in older B. saida biomass are documented for the mouths of Barrow and Mackenzie canyons, but with occurrence generally all along the slope of the Beaufort Sea (Crawford et al., 2012;Logerwell et al., 2018). Evidence for B. saida concentrations west of the Canadian Arctic Archipelago into McClintock Channel could be suspected based on observations of belugas echolocating for prey in the AW layer in that region (Carmack pers. obs.). In the vertical direction, the band of enhanced polar cod biomass in the Beaufort Sea is located where the relatively warm PSW and upper Atlantic Layer encroach the bottom along the continental slope (Parker-Stetter et al., 2011;Crawford et al., 2012;Majewski et al., 2016). This peak, reaching estimates of >150,000 individuals hectare −1 (August; Parker-Stetter et al., 2011), is in part driven by the vertical segregation of age classes with large, adult fish aggregating in those layers, albeit with some seasonal variations (Geoffroy et al., 2015). Peak densities and/or biomass coincided with waters above 0 • C (Crawford et al., 2012;Majewski et al., 2016), suggesting a possible role of this comparatively warm upper slope habitat from a thermal optimum perspective (Drost et al., 2014), in conjunction with zooplankton prey concentrations. On the Eurasian inflow slope, no comparable concentrations of B. saida have been reported, but several observations perhaps indirectly point to concentrations of B. saida abundance over Eurasian slopes: First, a mesopelagic layer encountering the seafloor at upper slope depths contains large zooplankton and a variety of fish species (250-600 m, Knutsen et al., 2017) includingand seasonally dominated by -B. saida (Geoffroy et al., 2019). And second, young Greenland halibut follow the continental slope and concentrate in the warmer (-0.6 to 1.2 • C) water layer between ∼400 and 780 m in St Anna and Voronin Troughs, where they actively feed on polar cod as their dominant prey (Dolgov and Benzik, 2017).
Continental slopes and shelf breaks with their associated frontal zones attract high numbers of foraging seabirds in many high latitude oceanic regions (e.g., Ainley and Jacobs, 1981;Baird and Mormede, 2014). Examples from temperate and sub-Arctic regions include: little auks in the NE Atlantic (Follestad, 1990), Cassin's auklet Ptychoramphus aleuticus in NE Pacific (Bertram et al., 2017;Wilkinson et al., 2018), and northern fulmars Fulmarus glacialis and fork-tailed storm-petrels Oceanodroma furcata in the Bering Sea (Schneider, 1982). While in the Arctic the role of slopes for higher trophic levels has not yet been summarized, we find that seabird accumulations at shelf break and continental slope area have in fact been documented for several species at both the Eurasian and Amerasian slopes. This is particularly the case for the small yet enormously numerous planktivorous Arctic endemic little auk, Alle. Little auks in Svalbard preferably fed far from the colony north of Svalbard in the marginal ice zone; with highest bird densities found over the slope area (concluded from data by Jakubas et al., 2017). Further east in the Kara and Laptev Seas, breeding distribution of the little auk also showed that nesting colonies were disproportionately located on coasts proximal to the shelf break (concluded from data by Gavrilo et al., 2011). This observation was consistent with Amélineau et al. (2016) who clearly demonstrated that slope waters off East Greenland were selected as profitable foraging area by the little auk irrespective of sea ice conditions, although these birds had to commute some 100 km from their breeding colonies to these preferred feeding grounds at the slope. These little auk distributions are influenced by the availability of large, lipid-rich copepod prey (Jakubas et al., 2017 and references therein), which we show can be concentrated in these slope regions. Attraction to shelf break areas was also documented for Ross's gull Rodostethia rosea by Hjort et al. (1997), who -during their Arctic Ocean crossing -found maximum concentrations of this species in the area where the shelf break is cut by St. Anna Trough. In the same area, Nansen already found these gulls to be common and locally numerous in the summer of 1895 (Collett and Nansen, 1900). For ivory gulls, the postbreeding staging areas are restricted to the marginal ice zone (Gilg et al., 2010), yet these staging areas also largely coincide with the continental slope from north of Svalbard to the east of the Severnaya Zemlya Archipelago (data from Gilg et al., 2010;Gavrilo et al., 2011). Examples of seabird concentrations at the slope of the Pacific inflow have also been documented. Total seabird abundance of both sea surface feeding bird species (kittiwakes Rissa tridactyla, short-tailed shearwaters Ardenna tenuirostris and Arctic terns Sterna paradisaea) and diving species (thick-billed murre Uria lomvia, horned puffin Fratercula corniculata and auklets Aethia spp.) peaked at the mouth of Barrow Canyon and the adjacent upper slope of the Beaufort Sea (Kuletz et al., 2015Moore and Kuletz, 2019). We pose that for some of these seabird species (e.g., the little auk) the shelf break might be potentially more important than sea ice distribution. This hypothesis draws attention to the slope area in the context of Arctic climate change, since topographically controlled and biologically important oceanographic features could be resilient even under conditions of a shrinking sea ice cover. To test this hypothesis, further studies of birds distribution patterns in slope areas are warranted.
Patterns in marine mammal distributions are largely consistent with the described patterns in bird observations at Arctic continental slopes. In the area of western Fram Strait and north off Svalbard bowhead whale and narwhal sightings were also associated with slope regions (Storrie et al., 2018;de Boer et al., 2019). Seasonally occurring whales such as blue whales were also seen at the shelf break (Falk-Petersen et al., 2015;Vacquié-Garcia et al., 2017) argue that shelf-break upwelling in that area (particular the part called Whaler's Bay) and associated dense occurrence of Calanus spp. copepods may have sustained the once high abundance of bowhead whales in that area. Similarly, extensive surveys for marine mammals along Amerasian Arctic slopes revealed several marine mammal hotspots associated with the continental slope of the Beaufort Sea and especially with canyon areas and strong oceanographic fronts (Moore et al., 2010;Okkonen et al., 2011;Moore and Kuletz, 2019). Specifically, hotspots for bowhead whales Balaena mysticetus, gray whales and beluga whales were located at the mouth of Barrow Canyon at the Chukchi Sea slope, or along the slope of the Beaufort Sea between Barrow Canyon and Mackenzie Canyon. For bowhead whales in particular, a shelfbreak front in the Barrow Canyon area, when present under certain conditions, is thought to enhance feeding conditions on zooplankton for those whales (Okkonen et al., 2011;Citta et al., 2015). Occurrence of beluga whale aggregations in the Beaufort slope region (Hauser et al., 2015;Stafford et al., 2018) has been linked to prey concentrations of B. saida and zooplankton also associated with local advection and upwelling processes . All of these findings suggest that slope and canyon bathymetry is a functionally important feature for many higher trophic level species.

Community Structure and Biodiversity Trends
As with biomass, the composition of biological communities (i.e., the members of all populations of species in a given area) vary strongly in their composition in the vertical and cross-slope dimensions, and to a lesser extent along-slope, both in the global ocean (Buhl-Mortensen et al., 2012) and in the Arctic. In all biological compartments, slope communities differ markedly from those on the shelves. Around the Arctic Basin perimeter, a vertical zonation of clearly distinct communities is identifiable across trophic levels as water depth increases from the shelf break to basin depths, from zooplankton and pelagic fishes in the water column to benthos and demersal fishes at the seafloor. For zooplankton communities, this vertical structure is evident near the slope but then extends horizontally throughout the basins, with epi-, meso-, and bathypelagic communities being correlated with vertical water mass structure ( Figure 6A; Kosobokova et al., 2011;Smoot and Hopcroft, 2017). In surface waters and into mesopelagic layers, the expatriate fraction of the community can be used to track the strength of influence of Atlantic (e.g., Oithona atlantica, Calanus finmarchicus, Meganyctiphanes norvegica) and Pacific (e.g., Neocalanus spp., Metrida pacifica) water masses as they spread along the slope perimeter (Hirche and Kosobokova, 2007;Kosobokova and Hirche, 2009;Nelson et al., 2009Nelson et al., , 2014Kosobokova and Hopcroft, 2010;Wassmann et al., 2015).
Benthic invertebrate communities also shift in community structure vertically (Figure 6B), with vertical shifts inherently tied to cross-slope gradients at the seafloor. As with zooplankton, taxonomic shifts often occur to other species or families within the same class or phylum, rather than to entirely different organisms at the phylum or class level. Down-slope community shifts in infaunal macrobenthos below the halocline, for example, are documented through changes in dominant polychaete and bivalve species on both Amerasian slopes (Chukchi and western Beaufort Sea: Bilyard and Carey, 1979;Pirtle-Levy, 2006;Grebmeier, 2012;Nelson et al., 2014) (Jørgensen et al., 2020) reflective of the highflow environment of the boundary current in this area (see section "Physical Oceanography of Arctic Slopes"), and different suspension feeders extend to the mid-slope (800-1500 m). Suspension-feeders are also more prominent in the Barrow Canyon outflow than on adjacent slope regions coinciding with the high velocity of transport of particle-rich water required to sustain them (Grebmeier, 2012;Pisareva et al., 2015). Outside the immediate inflow regions, slope regions tend to have high proportions of deposit-feeding taxa Bell et al., 2016) reflecting fine sediment organic carbon deposition on slopes from shelves. On the lower slope to the rise, various smaller sponges, sea anemones, crustaceans, and sea cucumbers are characteristic epifaunal taxa at both Eurasian and Amerasian slopes. Glacial dropstones reach slopes and adjacent basins alike and house higher diversity than the otherwise dominating soft sediments Meyer et al., 2016;Zhulay et al., 2019).
Communities of demersal fish on the slopes also shift in taxonomic composition from the adjacent shelf areas on both Eurasian and Amerasian slopes with differences also between upper and lower slope communities. Characteristic for the upper slope are high contributions of polar cod, the snail fish Liparis tunicatus and various eelpout species on the Pacific inflow slope (Beaufort Sea slope, Rand and Logerwell, 2011;Majewski et al., 2017;Norcross et al., 2017). At the Atlantic inflow, larger fish species are characteristic of the upper slope, such as Greenland halibut Reinhardtius hippoglossoides, Atlantic cod Gadhus morhua and the redfish Sebastes mentella (Bergstad et al., 2018;Jørgensen et al., unpubl. data); in warm years these species Non-metric multidimensional scaling plots show community shifts for shelf (light blue), upper slope (mid blue), and (partial) lower slope. For meso-zooplankton each circle represents one discrete depth layer of a multinet haul overall covering surface to near-bottom depth; stippled lines denote three sampling expeditions for the most part representing the three indicated slope areas. For benthic epifauna, each circle represents a single demersal trawl haul for a shallower Amerasian shelf-slope transition (upper panel; beam trawl) and a deeper Eurasian shelf-slope transition (lower panel; Campelen shrimp trawl); trawl depth was limited to ca. 1000 m.
even reach as far as the Kara Sea slope (Dolgov, 2013). On the mid to lower slope Arctic skate Amblyraja hyperborea, and the snail fish Liparis fabricii join Greenland halibut and the eelpouts, in both Atlantic and Pacific inflow areas (Beaufort  Dolgov, 2013). The major transition in fish species composition at the Atlantic inflow area north of Svalbard occurs between the slope (>∼500 m) and the upper slope/shelf break assemblages (<500 m) with more Arctic species in the deeper and more boreal species in the shallower community (Bergstad et al., 2018). A similar transition is equally distinct on the Amerasian slopes and interior Kara Sea slope, though it occurs shallower (matching the shallower shelf break ∼200 m), with typical Arctic shelf taxa such as Myoxocephalus, Artediellus, and other sculpins, alligatorfish Aspidophoroides (formerly Ulcina) olrikki and related species disappearing at upper slope depths Dolgov, 2013;Majewski et al., 2017;Norcross et al., 2017). Several fish species of the slopes extend into the basin (Stein et al., 2005;Zhulay et al., 2019), though few surveys have sampled beyond 1500 m due to logistical and time constraints. A difference between Amerasian and Eurasian slopes is that commercial densities of certain fish and shrimps reach the Eurasian inflow slope north of Svalbard. Finally, biodiversity and biogeography patterns across some ecological groups also show distinct gradients from the shelf down the slope, which may have implications for potential dispersal trajectories of new and potentially invasive species entering the Arctic Ocean. Zooplankton biodiversity peaks beyond the shelf break and in vertical layers between 200 and 2000 m over slopes (and in the basins) Kosobokova, 2012;Smoot and Hopcroft, 2017;Ershova and Kosobokova, 2019). The lower slopes stand out in that endemic zooplankton species begin to be more numerous than on the shelves (Markhaseva, 1998;Markhaseva and Kosobokova, 1998;Andronov and Kosobokova, 2011;Kosobokova et al., 2011). As discussed above, expatriate zooplankton initially follow boundary currents, but then decline with increasing distance from the inflow source areas. For benthos, diversity may also be higher at the upper-and mid-slope than on adjacent shelves and basins (Wlodarska-Kowalczuk et al., 2004;Vedenin et al., 2018), and the fraction of endemic species increases at the lower slope (Bluhm et al., 2011;Mironov et al., 2013). Similarly, fish species richness is higher at the upper slope than adjacent shelf and slower slope in the Beaufort Sea (Majewski et al., 2017) and northern Kara Sea (Dolgov and Benzik, 2017). Biogeographically, the Arctic continental slope seafloor is primarily the domain of species of Atlantic-Arctic affinity, cosmopolitan deep-sea and -at the upper slope -amphiboreal species, while Pacific-Arctic species are largely absent in and below the Atlantic Layer, even on the Amerasian side (benthos: Bilyard and Carey, 1979;Krylova et al., 2013;Zhulay et al., 2019;Ravelo et al., 2020;fish: Dolgov, 2013;Mecklenburg et al., 2016).

FUNCTIONS OF THE SLOPE: SYSTEM CONNECTIVITY VS. SEPARATION
In this section we show that the geomorphological setting and physical processes reviewed in section "Physical Oceanography of Arctic Slopes" define the ecological functions that govern biological patterns described in section "Gradients in Biological Communities at Arctic Slopes." These functions are a consequence of the slope's role as both a connector and as a separator, in each of the three dimensions, along-slope, cross-slope and vertically (cf. Figure 1C).

Along-Slope: Boundary Currents as Sub-Arctic Messenger
Along-slope connectivity of energy, material properties and organism transport around the basin perimeter and partly into the basin interior is primarily mediated through the ACBC, and to a lesser extent through the various and highly variable Pacific inflow boundary currents. This "plumbing system" effectively enhances along-slope transport, a key functional feature affecting the Arctic Ocean's sea ice, climate and ecosystem structure at and far beyond the slope region . Given the way this flow is dynamically constrained to follow isobaths, it follows that the steeper the slope, the more tightly constrained and intense the transport (e.g., Polyakov, 2001;Aksenov et al., 2011). For example, the Laptev Sea is steeper, and the boundary current is hence more constrained than at, for example, the East Siberian Sea slope.
For ecosystem productivity, the boundary current system is a very effective conduit of sub-Arctic nutrients, detrital particles, and living phytoplankton, zooplankton and fish larvae (Nelson et al., 2009;Torres-Valdés et al., 2013;Gjøsaeter et al., 2017;Knutsen et al., 2017;Hop et al., 2019). The species composition within this flow changes along-slope, with numbers of Atlantic/Pacific taxa decreasing and proportions of Arctic ones increasing. Given that the copepod Calanus finmarchicus is incapable of reproduction during its transit along-slope (Hirche and Kosobokova, 2007) biomass gradually diminishes through consumption and natural mortality . At the slope north of Svalbard, for example, its proportion is ∼40% while it drops to 5% or less at the New Siberian Islands (Kosobokova, 2012); low numbers of C. finmarchicus reach as far as the East Siberian Sea slope . Estimates of the amount of biomass injected into the slopes at the inflows are on the order of 1.75 million tons of zooplankton DW year −1 (∼0.8 million tons C year −1 ) in the Pacific inflow (Springer et al., 1996;Wassmann et al., 2015). Estimates for the Atlantic inflow were even higher than that (0.5 million tons C year −1 for C. finmarchicus alone, Basedow et al., 2018, or 0.9 million tons C year −1 of C. finmarchicus and C. glacialis combined, from model estimates, Wassmann et al., 2015). It is the persistent through-flow of these zooplankters that ensures plentiful feeding opportunities for their predators (Wassmann et al., 2019).
Several factors contribute to variability and loss of material transport along the boundary current path and these pertain, one, to the water transport, and two, to the transport of its material load. Loss of water and biomass from the boundary current occurs at intersections with ocean ridges (e.g., the Lomonosov, Alpha-Mendeleev, and Northwind Ridges), through flow instabilities that eject water into the basin interior (e.g., eddies, intrusion), and through relaxed slope steepness (releasing the tightness of the flow). The combined effects of these contribute to the downstream weakening of the boundary current. For example, the intersection of the Lomonosov Ridge with the slope and reduced slope angle towards the Amerasian Basin serve to bifurcate the boundary current and send a fraction across the basin as part of the Trans-Polar Drift (Rudels et al., 1994;McLaughlin et al., 1996). Likewise, the complex topography of the Chukchi Borderland acts to bifurcate the boundary current into two branches, one flowing north and the other through a gap south of the Northwind Ridge. Similarly, the PW inflow and its material loads may be weakened or even reversed by the clockwise (i.e., counter-opposing) wind-driven Beaufort Gyre and plumes of slope-constrained PW may be lost to the Amerasian Basin interior (e.g., Shimada et al., 2006). The magnitude of pelagic biomass transport within the boundary current on in the Atlantic inflow is seasonally variable which is related to seasonally changing vertical distribution of these zooplankton rather than variations in water transport (Basedow et al., 2018). Model simulations suggest that variability is also possibly related to variation in upstream source regions of the zooplankton (Wassmann et al., 2019). In the Pacific inflow changes in water temperature associated with sea ice retreat correlate with zooplankton biomass in the boundary current at the Chukchi Sea slope, and into the Canada Basin (Itoh et al., unpubl. data). Grazing on the algal bloom, natural mortality, consumption of copepods by fish, seabirds and mammals and other trophic relationships all also contribute to reduced biomass transport in the boundary current as the distance from the inflow decreases.

Cross-Slope "Leaks" Facilitate Connectivity
Earlier we outlined that a number of cross-shelf processes manage to break through the barriers of the along-slope fronts. One of these is shelf-break upwelling, a globally common phenomenon (Kämpf and Chapman, 2016) where surface waters above the shelf are driven offshore by upwelling favorable winds -in the Arctic generally easterlies -or by ice drift, to be replaced by deeper waters overlying the slope that are drawn onshore (Figure 4C). In the Arctic setting, halocline and upper AW is moved upwards cross-slope onto the shelves, facilitated by the rapidly changing water depth above the slope during situations of upwelling-favorable wind (Carmack and Chapman, 2003;Cottier et al., 2007;Pickart et al., 2009Pickart et al., , 2013aOkkonen et al., 2009;Williams and Carmack, 2015;. Convergence or divergence of isobaths along the slope can also result in locally intensified upwelling through acceleration or deceleration of isobath-following flows . In particular, submarine canyoncutting slopes result in very sharp changes in the direction (curvature) of isobaths, and the dynamics of isobath-following flows and thus can intensify upwelling and downwelling across the shelf-slope boundary (Carmack and Kulikov, 1998;Williams et al., 2006;. Deeper shelves like the Barents Sea, in contrast, are less likely to experience upwelling . The topic has been a key research focus around the Arctic due to strong implications for heat budgets, nutrient transport and primary production (cf. Tremblay et al., 2014). Phytoplankton production is enhanced by upwelled nutrients in the Beaufort Sea (Tremblay et al., 2011(Tremblay et al., , 2014 and such blooms are thought to boost grazing zooplankton biomass (Tremblay et al., 2011;Falk-Petersen et al., 2015) and benthic carbon demand (Tremblay et al., 2011). Falk-Petersen et al. (2015 hypothesize that upwelling with subsequent zooplankton biomass increase could be partly related to high bowhead whale numbers during the whaling period on the northern Svalbard shelf break. The current literature suggests that the relative roles of advected vs. upwelled nutrients and particulate constituents shift along the slopes from a large role of the former at inflow slopes to a larger role of the latter at interior slopes. Downwelling, in contrast, occurs when westerly winds over open water or drifting ice force surface waters above the shelf and slope onshore, thus drawing sub-surface shelf waters from offshore over the slope and into the basin (Sverdrup et al., 1942). Functions associated with downwelling are less well understood and draw less attention than those of upwelling. If sustained over a summer, ice-free period, downwelling may precondition winter conditions by forcing low salinity and nutrient poor waters onshore and constraining the seaward spreading of incoming river waters. If sustained over a sufficiently broad region, downwelling may thereby help the Riverine Coastal Domain (cf. Carmack et al., 2015) to remain intact along the coastline; thus, we hypothesize, supporting along-shore transport of material properties and coastal organisms including meroplankton and fish (cf. Griffiths et al., 1983;Craig, 1984;Fechhelm et al., 2007). Similar to upwelling, we suspect this process may be more prominent along narrow and shallow shelves.
In addition to wind-driven upwelling and downwelling, density-driven flows resulting from brine drainage during sea ice formation also provide a connection pathway between the shelves and the slope domain (Aagaard et al., 1981;Melling and Lewis, 1982;Schauer et al., 1997; Figure 4D). Not only do these waters provide important contributions to the halocline complex and potentially to the ventilation of the deep Arctic Ocean (Aagaard et al., 1985;Aagaard and Carmack, 1994;Rudels et al., 1994), but they also drive contour-following currents along the slope with the descending plume turning to the right by the Coriolis force . Dense water formation is particularly prominent in polynyas, sea ice formation factories that occur in all coastal regions around the Arctic (Ito et al., 2015) where heat flux and frazil ice formation are great (Figure 4D), and in shallow water where less dilution takes place. Densification of surface waters and subsequent sinking and off-shelf flow outside polynyas are episodic and difficult to observe, but also occur. Dense water production hotspots include the northeastern Kara Sea and northwestern Laptev Sea, where modeling studies (Ivanov and Golovin, 2007) discussed cascading of dense shelf waters down the continental slope to a depth of ∼400 m. The depths to which these winter-formed dense waters can reach depend on plume entrainment rates, upper ocean salinity prior to freeze-up (i.e., preconditioning) and ice production rates, which are both quite high in that region. However, with their close proximity to the large Arctic rivers, surface waters are becoming increasingly fresher (Janout et al., 2020, this issue), and thus even strong ice formation rates may not sufficiently densify these shelf waters to flow far down the slope. Densification and drainage are especially complex near the numerous topographic features (submarine valleys, troughs, canyons) that are too deep to locally form dense water by brine convection, but which can collect dense water and serve as drainage pathways (Janout et al., 2017).
Physical processes facilitating downwelling (dense-water formation and wind) can transport organic and inorganic material to depth, thus enhancing the biological carbon pump in the shelf-break region. We suggest that fixed carbon and organic detritus, both in surface waters and in advected waters associated with the ACBC, can in fact exhibit enhanced export from surface and intermediate waters to the deep Arctic basin, where it can be sequestered for the long-term. Removal of fixed carbon, as well as dissolved CO 2 by these physical processes should increase the air-sea gradient in inorganic carbon and lead to greater uptake of CO 2 from the atmosphere. The relative importance of these processes and the inverse, where upwelling brings inorganic carbon and perhaps even organic material to the surface, has not been addressed on a basin-wide scale.
Organic and inorganic constituents are also transported cross-slope with transport pathways associated with turbidity plumes, eddies interacting with the seafloor, and internal waves (Fohrmann et al., 2001). This cross-slope transport is -on Arctic and global slopes -reflected in enhanced content of total organic carbon in the sediment associated with a higher percentage of silt and clay sediments as fine particles settle out on the slope regions (Grebmeier et al., 2006;Mathis et al., 2014). Depending on location around the basin perimeter, composition of organic matter delivered to slopes has differing proportions of marine vs. terrestrial sources, with higher terrestrial fractions in areas where river influence reaches the slope (Fahl and Stein, 1997), i.e., generally in interior slopes, thus affecting food quality for benthic consumers (Bell et al., 2016). For example, water from the Mackenzie River, the Arctic river with the highest load of terrestrial organic matter (Rachold et al., 2004), provides a carpet of presumed refractory material that reaches the Beaufort Sea slope (Magen et al., 2010;Bell et al., 2016). At inflow (upper) slopes, in contrast, marine-derived carbon has a larger role (Magen et al., 2010;Divine et al., 2015). That this material is efficiently incorporated into the food web is, for example, also indicated in enhanced biological oxygen consumption in the sediments of the upper slope in and near Barrow Canyon and reaching down-slope with even a small increase at depths just past 2000 m (Moran et al., 2005;Grebmeier et al., 2006;Grebmeier and Barry, 2007, their Figure 6). In addition to this carbon supply, internal waves can mobilize sediments (Fahrbach and Meincke, 1982), thereby contributing to structuring seafloor communities since they respond to grain size composition (Ravelo et al., 2015).
We note that not all forcing mechanisms have uni-directional outcomes, for example upwelling and downwelling winds act primarily in the cross-slope direction, but will also have an along-slope function by generating a transient, along-slope jet in the direction of the wind along the upper slope ; Figure 4C). This is because over the upper slope there is less water for the wind to accelerate than in deeper waters farther offshore, and so the acceleration is greater there, leading to faster flow (the transient jet). This effect is more pronounced for shallow shelf-breaks where there is less water to accelerate over the upper slope . Functionally, the formation of transient, shelf-break jets associated with upwelling and downwelling may transport biota, and winnow and re-suspend sediments along the upper slope (cf. Williams and Carmack, 2015), again with effects on benthic community composition given their substrate affinity.

Vertical Connectivity at the Slope by Differential Weakening of Vertical Stratification
In general, halocline stratification throughout the Arctic Basins strongly constrains the vertical exchange of heat, nutrients, oxygen, biogeochemical tracers and other properties (McLaughlin and Carmack, 2010;Nishino et al., 2019;Randelhoff et al., 2020). At the slopes, however, stronger currents infer greater vertical shear and thus stronger mixing between AW and the overlying halocline compared to adjacent basins, and this results in a relative weakening of vertical gradients of temperature and salinity at the slope (Polyakov et al., 2020). The rapidly decreasing intensity of horizontal currents with depth, i.e., strong vertical shear, is a robust feature of the pan-Arctic along-slope boundary current; examples include: northeast of Svalbard (Ivanov et al., 2009), St. Anna Trough (Dmitrenko et al., 2015, the eastern Eurasian Basin slope , and the Alaskan Arctic slope (Li et al., 2019). Moreover, baroclinic tides are amplified over steep topography (e.g., Rippeth et al., 2015;Luneva et al., 2015), and thus also contribute to strong vertical shear. Stronger vertical shear of mean and tidal currents over the Arctic slopes is one important contributor to enhanced mixing, for which the characteristic signature is reduced vertical stratification within the halocline and AW layers . Weakening of oceanic stratification in the halocline thus allows greater vertical transfer of AW heat and other scalar properties. These patterns of weakened temperature and salinity gradients represent proxies for exchange of nutrients and other biogeochemical variables through the halocline. For example, the progression of the boundary current eastward displays a cooling of the AW (section "Along-Slope: Boundary Currents") and a concurrent warming of the halocline, implying heat transfer upwards (Polyakov et al., 2010). We suggest that this process is accompanied by vertical flux of nutrients and other biogeochemical properties, and hence helps enhance levels of new production above the slope regions. These fluxes are not trivial to estimate given that uptake by phytoplankton, nitrification and horizontal advection may seasonally contribute to overall fluxes (Randelhoff et al., 2015). These authors argue, however, that the vertical nutrient flux component dominates during fall and winter when light levels limit phytoplankton growth, but that the relative roles of vertical and advective fluxes would vary along the Arctic slopes. They further suggest that this vertical nitrate replenishment in winter could potentially support increased production with decreasing sea ice in the inflow slopes in the future (also see section "Towards a Future Arctic Slope System").

Fronts Provide Biological Concentration Mechanism Along and Cross-Slope
Oceanic frontal zones are natural water mass boundaries often expressed by steeply sloping isolines of temperature and salinity; fronts are typically maintained by either divergent or convergent water movement. The Arctic fronts of interest here are (1) alongslope fronts created by the density structure associated with the ACBC, (2) cross-slope fronts occurring at sites of subduction of incoming sub-Arctic waters, or of lateral injection of subsurface waters at canyons (St. Anna Trough and Barrow Canyon), and (3) front associated with shelf-break jets, both surface and near-bottom. Due to the sustained intensity of the boundary current system these frontal zones are maintained and can act as an effective barrier between the abyssal ocean and shelves (Tverberg and Nøst, 2009). For example, analysis of 2013-2015 cross-slope mooring observations in the eastern European Basin revealed that currents, even in the very surface layer, are mostly aligned with the underlying topography throughout all averaging periods, regardless of wind direction (Baumann et al., 2018). Such fronts are found around the pan-Arctic slopes and appear to facilitate concentration of pelagic prey, as is also well-described for fronts globally (Pakhomov and McQuald, 1996;Genin, 2004;Trudnowska et al., 2016). As a result, horizontal property gradients establish hydrographic fronts that are observed along the AW pathway following the continental slope of the Eurasian Basin. The front at the eastern flank of the St. Anna Trough extends vertically throughout the entire water column with a strong horizontal density gradient  and is associated with bird concentrations (section "Gradients in Biological Communities at Arctic Slopes"). Similarly, Bauch et al. (2014) observed a front at the continental slope of the Laptev Sea separating shelf and basin waters, and hypothesized that the front is maintained by the ACBC. Also, a strong front near the shelf-break at the southern edge of Barrow Canyon is maintained under certain conditions and concentrates zooplankton which in turn attracts feeding groups of bowhead whales (Okkonen et al., 2011). Such zooplankton concentrations, explained through sharp density gradients (Prairie et al., 2012;Trudnowska et al., 2016), in turn explain the widespread aggregation of certain planktivorus seabirds and marine mammals (section "Gradients in Biological Communities at Arctic Slopes") at the Arctic inflow slopes.

Light, Hydrography, and Depth Gradients Maintain Vertical Domains
While vertical connectivity is enhanced over slopes compared to the basins (section "Vertical Connectivity at the Slope by Differential Weakening of Vertical Stratification") we briefly stress here that for the food web at slopes, vertical gradients in light, food supply and hydrographic properties remain strong, maintaining strong functional gradients in the vertical. The combined effects of primary productivity constrained to upper water, consumer biomass peaks in form of a 'lipid belt' of enhanced zooplankton and -regionally -polar cod, strong vertical flux attenuation prevalent anywhere in the global ocean, and the strong belts of boundary currents, the upper slope is clearly separated from the lower slope. Within it, the euphotic zone is distinct from the basins through enhanced mixing as shown in "Vertical Connectivity at the Slope by Differential Weakening of Vertical Stratification." Clearly separated, the lower slope is generally more characteristic of global deep-sea conditions in terms of low food supply, a community structure of more detrital and predatory taxa, and hydrographic stability. The described down-slope processes and possibly enhanced carbon pump, however, set this zone apart from the adjacent basin conditions.

TOWARDS A FUTURE ARCTIC SLOPE SYSTEM
In summary, we have shown that the Arctic continental slope is a distinct and dynamically active domain of the Arctic Ocean, manifest as narrow, horizontal bands with specific combinations of physical, chemical and biological properties that encircle the basin perimeter. Though spatially occupying only a small fraction (∼6%) of the Arctic Ocean surface area, this domain exerts a disproportionately large influence on the overall functioning of Arctic Ocean ecosystems. Along-slope advective inputs from sub-Arctic seas create and maintain gradients, especially at the upper slope, in water mass and biotic properties moving along the basin perimeter in boundary currents that facilitate connectivity. The associated near-continuous carbon injections moving along-slope support multiple trophic levels in the water column , but also enhance food supply to the underlying seabed. The lower slope, in contrast, is a quiescent, low energy environment, less susceptible to climate forcing, and with comparatively low biological stocks and activity. Connectivity, both cross-slope and vertical, is generated by processes cutting through the boundary currents such as shelf break upwelling/ downwelling and brine drainage, particularly through canyons (Figures 4C,D). Transport of organic and inorganic matter is associated with these processes, providing a mechanism to enhance biological productivity subsequently attracting multiple trophic levels. Vertical connectivity also results from generally weaker stratification at the slopes than in the adjacent basins through shear and turbulent mixing, enhancing upward nutrient fluxes. At the same time, however, salt-stratification, light attenuation (determining the euphotic zone) and water depth (indirectly affecting vertical flux of organic matter) maintain distinct hydrographic layers, biological communities and production regimes. Water mass boundaries at fronts along-slope and cross-slope near canyons tend to be places where biological concentrations occur, often evident in predator-prey associations.
These factors contribute to the upper Arctic slope domain exhibiting signals of a changing climate rapidly and acutely because the upstream signals of change originating in sub-Arctic source waters and inflow shelves are rapidly propagated by the ACBC (Polyakov et al., 2020a, this issue). For example, the Arctic slope has warmed, as shown by an increase in the AW core temperature since 1980 (Figure 7). This trend, however, is not spatially uniform, as warming rates in the western Eurasian Basin (0.03 • C per decade) exceed those in the Amerasian Basin (0.01 • C per decade) by a factor of 2-3, although along this track there has been little change of water transports over the past two decades -2018Pnyushkov, Polyakov et al., pers. com.). Freshening is also visible in the AW core over the slope regions of the Amerasian Basin while AW salinity over the Eurasian Basin slopes shows a uniformly increasing trend (Figure 7), consistent with the ongoing Atlantification of that part of the Arctic Ocean (e.g., Polyakov et al., 2017Polyakov et al., , 2018. Polyakov et al. (2020a, this issue) also documented that the halocline in both the Eurasian and Amerasian basins has shown a clear warming tendency, while salinity trends differ, increasing in the Eurasian Basin whereas the Amerasian Basin halocline has become increasingly fresher. Thus, the two basins and their slopes are structurally responding in opposite directions, reflected in weakening halocline stratification over the Eurasian Basin slope and strengthening over the Amerasian Basin slope areas. As a result, more favorable conditions are now established for higher biological productivity at the Eurasian Basin's margin while at the same time, conditions at the upper slope in the Amerasian Basin experiences increased constraints to vertical mixing and the flux of nutrients to the surface layer. The increased influx of warmer PW through Bering Strait since the early 1990s (Woodgate et al., 2006;Woodgate, 2018) implies enhanced along-slope PW transports of nutrients, pelagic biomass and propagules of benthic biota in the Amerasian Basin. There is, however, no direct observational evidence for this intensification above the slope. Benthic indicator species are, in fact, now actually distributed less far to the east than previously (Ravelo et al., 2015) which these authors suggested could be related to prolonged reversals of water flow from east to west reported by von Appen and Pickart (2012). Instead, sea-ice retreat in summer and stronger air-sea coupling may be detaching a greater fraction of this water from the shelfbreak current into the basin interior .
Sea ice decline, in particular, has acted to accelerate the previously existing processes of shelf-slope exchange, including shelf-break upwelling and downwelling (Carmack and Chapman, 2003), cross-slope flow of brine-enriched dense water (Janout et al., 2017), and downslope transport of organic matter (Fahl and Stein, 1997). Now, the seasonal ice edge retreats well FIGURE 7 | Increase of AW core temperature (defined by potential temperature maximum) and associated salinity from 1980 to 2015 over the Arctic slope from 1980 to 2015. Locations are shown by red transects. Solid blue lines connect annual measurements with no gaps in between whereas dash-dotted lines are used to fill gaps. Light blue dots show original data taken from CTD profiles. Red lines show linear trends, their values complemented by statistical significance at 95% are also shown. beyond the slopes into the basin in many regions (Figure 8), whereas only a few decades ago, seasonally ice-free waters were limited to near-shore regions along the periphery (Stroeve et al., 2012;Kwok et al., 2019). The decline in ice cover exposes waters and drifting ice above the slope to enhanced wind and radiative forcing, and thus both upwelling and downwelling favorable winds can now drive increased cross-slope exchange. This trend of increased exchange has demonstrated the potential to increase the flux of nutrient-rich, sub-Arctic Atlantic and Pacific waters onto the shelf, and potentially increase primary production (Tremblay et al., 2011;Spall et al., 2014).
Multiple processes were mentioned that may increase primary production over the slopes. To evaluate such potential changes on the pan-Arctic scale, we used the regional coupled physical and biological SINMOD system, run with atmospheric forcing data from the global MPI-ESM model system (Notz et al., 2013). Except for this forcing, the configuration is otherwise similar as presented in Slagstad et al. (2015). To eliminate interannual variability, the change is calculated by subtracting the decadal mean annual primary production over the period from 2090 to 2099 from the mean annual primary production over the period from 2006 to 2015. Simulated yearly primary production in the Arctic is particularly sensitive to how well the sea ice is represented by the ESM (Slagstad et al., 2011 and therefore the future projection of new production depends on future projection of summer sea ice cover. An ensemble simulation of the IPCC RCP 8.5 scenario that predicts a sharp decline in ice cover in September after 2050s (Notz et al., 2013) alters primary production (Figure 9). The magnitude of the change in production will depend on the forcing scenario, but qualitatively the results are consistent whether we run with forcing from other FIGURE 9 | Modeled change in annual primary production by the end of the century compared to the present climate based on SINMOD results for the RCP8.5 scenario. Predictions show peak increases (oranges and reds) partially coinciding with the continental slopes along the Atlantic inflow and Siberian slopes and enhanced production (yellow) over essentially the entire slope band.
EMSs or do more idealized blue Arctic scenarios . The results can be summarized in the Eurasian Arctic as an increase in new primary production in the northern Barents Sea and along the inflow path of AW north of the Barents Sea and further along the western Eurasian slope. Along the western Eurasian slope there is a projected increase of 20-40 g C m −2 year −1 from today's values of 70-100 g C m −2 year −1 with peaks in increase projected for the eastern Barents and western Kara seas. The Amerasian slope values are projected to increase by 10-20 g C m −2 year −2 from today's values of generally <30 g C m −2 year −1 , with highest projected increase for the Laptev and East Siberian and lowest for parts of the Beaufort Sea slopes. With the projected change in ice conditions, the productive season will be longer though phytoplankton growth may eventually become nutrient limited. The strong stratification in parts of the Arctic constrains turbulence levels (Randelhoff et al., 2015), and future increases in new production will, therefore, depend on an increase in vertical fluxes of nutrients. SINMOD projects this to occur mainly in the northern Barents Sea and along the Siberian slopes due to weakening of the halocline and deeper mixing with underlying AW above the slope, a trend that is already taking place in these regions (Polyakov et al., 2017;Lind et al., 2018). These projected patterns are rather consistent with the patterns in documented increases in primary production and algal biomass in the last decade (Lewis et al., 2020), where these authors in fact suggest that regional influx of new nutrients played a role. Simulations also show a high projected increase in new production in the Chukchi Sea related to changes in ice conditions, a trend that again is in agreement with satellite based observations (Arrigo and van Dijken, 2015;Frey et al., 2019). Whether and to what extent this new production would reach the slopes and how that might relate to the increased volume transport of nutrient-rich water through Bering Strait from 0.7 to 1.2 Sv since the 1990s (Woodgate, 2018) warrants a thorough analysis.
Microbial and consumer communities above continental slopes are presented with changed proportions of carbon sources as well as temperature regimes. A new mix of carbon sources can be expected given the combination of the predicted increase in primary production, decrease in ice cover (and hence likely ice-algal production), potentially changed cross-slope transport of organic carbon delivered with sea ice, river run-off and permafrost (Holmes et al., 2002(Holmes et al., , 2012Krumpen et al., 2019), and perhaps even macroalgal carbon reaching some slopes from nearby island groups with rocky shores. Overall, consumer responses are difficult to predict. Some changes would support increases in consumer production and biomass: increased pelagic primary production, upwelling, higher temperatures and associated enhanced survival and subsequent reproduction of advected zooplankton. In contrast, other changes might reduce production capacity, e.g., smothering by down-slope inorganic particles, increased metabolic demands related to temperature increase, and reduced food quality with increased fractions of terrestrial matter. In addition, the consumer communities at slopes themselves are beginning to change, at least near the sub-Arctic inflows. Penetration of sub-Arctic Atlantic communities farther into the Arctic along-slope is documented for inflow shelves (e.g., Fossheim et al., 2015), and is now emerging also for at least individual surveyed species at the matching slopes. For example, recent studies found Greenland halibut Reinhardtius hippoglossoides to be distributed much further eastwards from the northern Barents Sea to the Kara Sea and Laptev Sea (Borkin et al., 2008;Sentyabov and Smirnov, 2010).
In conclusion, the distinct functions of the Arctic slope domain and the ongoing changes in this domain should be incorporated into future Arctic conceptual models and research planning. Based on the role of the slope in climate changes, both as an early warning network and as pathway of upstream change to the entire Arctic marine system, we recommend an interdisciplinary, international slope study be conducted in order to synoptically characterize connecting and separating processes over slopes. The developing international Synoptic Arctic Survey (SAS) for example, composed of regional shelfto-basin transects for multidisciplinary studies into the Arctic Basin and recommended interdecadal follow-up along timeseries lines (Paasche et al., 2019) can yield valuable pan-Arctic slope information on the status and change of the Arctic marine ecosystem also along slopes. Given the growing appreciation for the role of the numerous troughs and canyons, vertical and lateral physical and biological exchange and mixing mechanisms of AW with peripheral waters are crucial study subjects to understand the future of the Arctic ecosystem. Also, a suite of slope transects should monitor the potential arrival of invasive species or those expanding their current distribution range as is already seen along the Atlantic inflow slope.

AUTHOR CONTRIBUTIONS
BB and EC conceived the idea and drafted the initial draft of the manuscript. BB led the writing of sections "Introduction: Motivation and Definitions" and "Gradients in Biological Communities at Arctic Slopes." MJ led section "Physical Oceanography of Arctic Slopes." BB, EC, and IP led section "Functions of the Slope: System connectivity vs. Separation." BB, IE, and EC led section "Towards a Future Arctic slope System." MJ made Figures 1A, 2