ORIGINAL RESEARCH article

Front. Mar. Sci., 19 November 2020

Sec. Global Change and the Future Ocean

Volume 7 - 2020 | https://doi.org/10.3389/fmars.2020.548931

Extraordinary Carbon Fluxes on the Shallow Pacific Arctic Shelf During a Remarkably Warm and Low Sea Ice Period

  • 1. College of Fisheries and Ocean Sciences, University of Alaska Fairbanks, Fairbanks, AK, United States

  • 2. Amundsen Science, Université Laval, Québec, QC, Canada

  • 3. Institute of Marine Sciences, University of Alaska Fairbanks, Fairbanks, AK, United States

Abstract

The shallow Pacific Arctic shelf has historically acted as an effective carbon sink, characterized by tight benthic pelagic coupling. However, the strength of the biological carbon pump in the Arctic has been predicted to weaken with climate change due to increased duration of the open-water period for primary production, enhanced nutrient limitation, and increased pelagic heterotrophy. In order to gain insights into how the biological carbon pump is functioning under the recent conditions of extreme warming and sea ice loss on the Pacific Arctic shelf, we measured sinking particulate organic carbon (POC) fluxes with drifting and moored sediment traps, as well as rates of primary production and particle-associated microbial respiration during June 2018. In Bering Shelf/Anadyr Water masses, sinking POC fluxes ranged from 0.8 to 2.3 g C m–2 day–1, making them among the highest fluxes ever documented in the global oceans. Furthermore, high export ratios averaging 82% and low rates of particle-associated microbial respiration also indicated negligible recycling of sinking POC in the water column. These results highlight the extraordinary strength of the biological carbon pump on the Pacific Arctic shelf during an unusually warm and low-sea ice year. While additional measurements and time are needed to confirm the ultimate trajectory of these fluxes in response to ongoing climate change, these results do not support the prevailing hypothesis that the strength of the biological carbon pump in the Pacific Arctic will weaken under these conditions.

Introduction

Arctic marine systems are currently undergoing rapid and profound changes due to the effects of climate change, including reduced sea ice extent, earlier sea ice retreat, protracted ice-free seasons, warming air and ocean temperatures, and shifts in currents and water column stratification (; ). These environmental changes have recently accelerated on the Pacific Arctic’s Bering and Chukchi Sea shelves (; ; ). During 2017 and 2018, bottom water temperatures in the Bering Sea were 3°C higher than the 2005–2016 baseline (), and the four lowest maximum sea ice extents since 1979 in the Bering and Chukchi Seas have occurred after 2015 (). Cascading impacts on the regional ecosystems, biogeochemical cycles, climate, and human communities on and around the Pacific Arctic shelf are expected, although the nature and magnitude of these impacts remain largely speculative (; ; ; ).

The shallow Pacific Arctic shelf, averaging 50 m depth, has historically acted as a strong sink of carbon (; ). Water movement on these shelves is generally northward carrying different water masses of Pacific origin into the Arctic (; ), with a significant seasonal modulation () (see Supplementary Figure 1 for visualization of currents). This region is also characterized by a strong biological carbon pump having pelagic primary productivity (; ), sedimentation (), and benthic productivity (; ) rates that are all amongst the highest measured in any marine system. When light is sufficient, pelagic primary productivity may reach up to 16 g C m–2 day–1 and 470–840 g C m–2 year–1 (; ) due to elevated nutrient concentrations (5–20 μM) () advected into the region with the Anadyr current from deep Pacific upwelling (; ). The spring phytoplankton bloom is typically dominated by large, rapidly sinking sympagic or pelagic diatoms (; , ) that contributed to an annual particulate organic carbon (POC) flux of up to 145 g C m–2 year–1 near Hanna Shoal () and likely facilitate substantial carbon burial in sediments. Similar to other Arctic shelves (), total organic carbon on the Pacific Arctic shelf averages 1% and reaches up to 2% in surface layers (; ; ). This active biological carbon pump supports large populations of benthic-feeding pelagic seabirds and marine mammals (; ), many of which are important to Indigenous communities that rely on subsistence hunting ().

The strength of the biological carbon pump in the Arctic is predicted to weaken with climate change due to increased duration of the open-water period for primary production and enhanced nutrient limitation (; ; ). Warmer waters could increase metabolic rates of pelagic grazers and heterotrophic bacteria and potentially favor smaller phytoplankton and faster-growing grazers that more rapidly recycle organic matter within the water column (; ). Additionally, increased frequency of storms could increase mixing and efflux of carbon dioxide (; ). If these predictions prove accurate, such mechanisms could accelerate feedback processes on the services supported by the biological carbon pump. While the effects of changing ice conditions and warming water on production on the Pacific Arctic shelf have been investigated for many years (; ; ), few studies have directly measured pathways within the biological carbon pump (Fukuchi et al.,1993; ). As part of the Arctic Shelf Growth, Advection, Respiration, and Deposition rate experiments (ASGARD) project, we measured the strength and efficiency of the biological carbon pump by directly quantifying and comparing rates of primary productivity, sinking POC flux, and microbial respiration associated with trap-collected sinking particulate matter during June of 2018 on the Pacific Arctic Shelf. This study occurred after the winter with the lowest maximum sea ice extent on record and in a prolonged time of abnormally warm water (; ; ). These unique environmental conditions, described in detail in , provided an opportunity to test the prevailing hypothesis that the biological carbon pump will decrease in strength with climate change (; ; ).

Materials and Methods

Hydrography Sampling

All cruise operations were performed on the R/V Sikuliaq. The CTD unit consisted of a Seabird SBE16plus unit coupled with WetLabs fluorometer and transmissometer. A Satlantic SUNA V2 instrument was also mounted to the rosette to measure nitrate. To characterize the water mass at each station, surface salinity and temperature data were retrieved from the CTD profiles at each station, then plotted on a temperature–salinity (T/S) diagram.

Drifting Sediment Trap Sampling

A standard Lagrangian-type surface-tethered drifting sediment trap (KC Denmark model number 28.200) was used to collect sinking particles () at seven locations (Figure 1). Two of the four tubes contained a removable clear-bottomed cup filled with 250 mL of viscous polyacrylamide gel. The cups were fitted with a thin sloping ramp to funnel all sinking particles into the gel within the cup and prevent particles from settling between the inside of the tube and the outside of the cup. All four tubes were filled with chilled (0°C) filtered seawater (0.3 μm) collected in Niskin bottles from the same depth and station at which the drifting sediment trap was deployed. The remaining two tubes collected sinking particles in bulk, maintaining in situ chemistry as much as possible. The trap array was deployed at 30 m below the surface at each station, estimated to correspond to the bottom of the euphotic zone, for 3–12 h depending on the timing of other cruise operations (Table 1). We used the same depth for consistency and to reduce issues of resuspension by sampling too close to the seafloor. The drifting sediment trap was fitted with an ARGOS beacon and a go-Tele GPS tracker unit to track its real-time location.

FIGURE 1

), where the 2018 June median ice extent (maroon) is compared with the 1981–2010 median June baseline (purple) and the record breaking 2018 maximum ice extent (olive) is compared with the 1981–2010 median winter maximum (teal).

TABLE 1

Station nameASGARD station nameBottom depth (m)Latitude deploy (degrees decimal minute)Longitude deploy (degrees decimal minute)Date and time of deployment (M/DD HH:MM UTC)Date and time of recovery (M/DD HH:MM UTC)Total time of deployment
ACBE14163° 18.1′−168° 27.0′6/07 15:056/08 2:4011h 35m
BDBO2.45064° 58.6′−169° 52.8′6/11 10:596/11 17:296h 40m
CDBO3.85067° 40.4′−168° 50.1′6/14 23:106/15 4:525h 42m
DIL44267° 28.3′−166° 12.5′6/13 11:576/13 20:448h 47m
EDBO3.34968° 11.1′−167° 18.6′6/15 19:306/15 22:553h 25m
FCL35169° 2.1′−168° 49.4′6/16 19:416/17 0:104h 29m
GCL14668° 57.3′−166° 53.8′6/17 21:186/18 2:205h 2m

Location and duration of drifting sediment trap deployments 30 m below the sea surface in 2018.

Flux Rate Measurements

Sinking particles collected in the drifting sediment trap were used to determine POC fluxes. Once the trap was recovered, the following steps were performed as quickly as possible in an environmental chamber that fluctuated in temperature from 3 to 5°C in order to maintain as close to in situ conditions for particle-associated microbes as possible. Overlying water was siphoned using a vacuum pump down to a boundary layer above the settled particles at the bottom of all four tubes. In the two bulk particle collection tubes, the material that remained in the tubes after siphoning was quantitatively split into four subsamples using a Folsom plankton splitter. Three subsamples were used for triplicate analytical flux measurements. These subsamples were filtered onto pre-combusted 25-mm Whatman GF/F filters and placed in a dehydrator at 60°C for 12 h. Once dried, the filters were sealed in Petri dishes until further analysis (See section “Elemental POC/PN/δ13C/δ15N Analysis” for details). Particulate organic carbon values were converted to daily fluxes depending on the deployment period and the collecting area (g C m–2 day–1).

Respiration Rate Measurements

One subsample from each of the two drifting sediment trap tubes was used to estimate particle-associated microbial respiration rates. This material was homogenized by swirling the container and pipetted with a wide-bore pipette into eight replicate 2-mL glass vials (Batch PSTS-1721-01) fitted with Pre-Sens Oxygen Optode Sensor Spots (Regensburg, Germany) per drifting sediment trap tube, totaling 16 experimental samples. Filtered seawater controls were obtained from a Niskin bottle closed at 30 m depth during a CTD cast upon recovery of the drifting sediment trap. Water for the control samples was filtered (0.3 μm) to remove particles and particle-associated microbes. The filtrate was pipetted into eight replicate vials that were identical to the experimental vials. All 24 vials were checked for air bubbles, and vials were then placed inside a sealable, clear, plastic water-bath and placed on top of a PreSens SDR SensorDish Reader The water-bath was located inside a dimly lit cold room that varied from 3–5°C and connected to a Fisherbrand Isotemp 500LCSU Circulator, now referred to as a chiller, which maintained temperature at precisely 4.0°C during the incubation. The concentration of oxygen in each vial and temperature in the incubation chamber were recorded every 30 s for the duration of the incubation using PreSens – Sensor Dish Reader Version 4 Software. Incubations lasted for between 3 and 12 h. After the incubations, the remaining material from each vial was filtered onto individual GF/F filters and treated the same as the flux measurement samples.

A few modifications were made to the methods used for measuring particle-associated microbial respiration rates during the course of the study in order to try to improve the accuracy of our measurements. During the first incubation at station A, a low-oxygen micro-environment formed around the sensor spot, located at the bottom of the vial. For all subsequent incubations, the entire incubation chamber was repeatedly inverted for 5 s every 3 min to mix the sample. Additionally, we noticed the concentration of oxygen increased over time in a few of the experimental samples (i.e., at stations A, C, and E), suggesting photosynthetic activity. During the last two incubations at stations F and G, a black cloth was used to cover the incubation chamber in order to prevent any light from reaching the samples, theoretically preventing light reactions associated with photosynthesis. However, it should be noted that dark reactions associated with photosynthesis can continue for several hours after the removal of light in cold water.

The data recorded by the PreSens software were downloaded and analyzed using the following steps in Matlab 2017a computing software. Data collected before the incubation temperature stabilized were trimmed so that only the time during which the incubation temperature remained stable was analyzed. The first 188 min of data after temperature stabilization was used to determine respiration rates. Linear regression analysis was performed on the oxygen concentration data from each vial. The average and standard deviation of the eight replicate control slopes (rcontrol) was taken, and for each of the two experimental samples. Then the average slopes for each of the experimental incubations were averaged together and the error was propagated (rexp). We calculated the carbon-specific particle-associated microbial respiration rate (R_PAM) using a 117:170 organic carbon to oxygen molar respiratory quotient (VOC:O2), assuming a one to one relationship with organic carbon degradation and carbon dioxide production (), an incubation volume (vol) of 2 mL, and the final concentration of POC at the end of the incubation ([POC]) (Eq. 1).

The average R_PAM ± 1 standard deviation was compared with those from other studies (; ; ; ; ; ) at different locations.

Sinking Particle Visualization

Collecting particles in polyacrylamide gel kept sinking particles intact and allowed for particle imaging and identification (; ; ). The contents of the cups were photographed within 6 h of sediment trap recovery using a 42.4 MP digital camera equipped with a 90 mm macro-lens and a flash unit. A length to pixel relationship was determined for each image in Adobe Photoshop CS6. These samples were used to qualitatively determine sinking particle type.

Primary Productivity Rate Measurements

13C-15N dual-isotope tracer technique was used to measure integrated rates of primary productivity at the seven stations following a standard protocol (; ). For primary productivity rate measurements, water was collected at six depths corresponding to the 100, 50, 30, 12, 5, and 1% light levels. The 1% light level is estimated to be the minimum amount of light necessary for photosynthesis to occur, i.e., the bottom of the euphotic zone. The incubations lasted between 4 and 7 h, and measurements were extrapolated to daily production by adjusting to total daylight for each incubation site. At the end of the experiment, contents in each incubation bottle were filtered onto pre-combusted GF/F filters and frozen at −80°C until further analysis. These depth-specific rates were then integrated over the entire depth of the euphotic zone to determine total water column primary productivity rates in units of g C m–2 day–1.

Elemental POC/PN/δ13C/δ15N Analysis

All dried or frozen GF/F filters were processed in the Alaska Stable Isotope Facility at the University of Alaska Fairbanks’s Water and Environmental Research Center. Filters were acidified with 10% hydrochloric acid for 24 h to remove particulate inorganic carbon (PIC). The filters were pelletized in tin cups. Stable isotope data were obtained using continuous-flow isotope ratio mass spectrometry. Stable isotope ratios were reported in δ notation as parts per thousand (‰) deviation from the international standards VPDB (carbon) and air (nitrogen). Typically, instrument precision was <0.2‰.

Moored Sediment Trap Sampling

Two 24-cup Hydro-Bios sediment traps were moored at stations B south of Bering Strait (trap depth 37 m, bottom depth 49 m) and C north of Bering Strait (trap depth 35 m, bottom depth 50 m) from June 2017 to June 2018 (see Table 2 and Figure 1 for location of traps). Sample cups were filled with a hyper-saline (salinity 38) 5% formalin solution in filtered seawater to preserve samples during and after deployment (). The carousel rotated at pre-programmed intervals ranging from seven to 40 days.

TABLE 2

Station nameTrap nameSampling periodTrap depth (m)Bottom depth (m)Latitude (degrees decimal minute)Longitude (degrees decimal minute)
BN4June 26, 2017 to June 08, 2018374964° 55. 7′−169° 55.1′
CN6June 17, 2017 to June 08, 2018355067° 40.2′−168° 44.7′

Location and duration of moored sediment traps.

As the trap at station B was recovered before the end of its rotation, the material in the last open cup was excluded from analysis. Subsamples from each cup were filtered onto pre-combusted (500°C overnight) GF/F filters (0.7 μm), exposed to 1 N hydrochloric acid overnight for removal of inorganic carbon, and dried at 60°C overnight before encapsulation for POC measurements (). Particulate organic carbon measurements were conducted on a Perkin Elmer CHNS 2400 Series II elemental analyzer. Particulate organic carbon measurements were converted to daily flux rates depending on the open cup duration of each sample.

Instruments measuring physical and biological parameters in tandem with sinking POC flux were deployed on the moorings at stations B and C. A 300 kHz RDI workhorse ADCP measured bottom current velocity ∼5 m off the seafloor at each site. Lowpass-filtered bottom current velocity were plotted overlaid with a 12-h smoothing. A Seabird SBE16plus unit coupled with a Wetlabs fluorometer measured temperature, salinity and fluorescence at 27 m at station B and 25 m at station C.

Remote Sensing

Daily sea ice concentrations were retrieved from the National Snow and Ice Data Center satellite records for the deployment period at both mooring sites (). Wind velocity was obtained from modeled wind reanalysis for the deployment period at both mooring sites (European Centre for Medium-Range Weather Forecasts, 2019). Twenty-four hour smoothing was performed on wind data.

Results

Environmental Conditions

Sea surface temperatures ranged from 1 to 10°C during the ASGARD expedition in June 2018, with the warmest water temperature above 8°C observed south of Nome and west of Norton Sound. These warm waters were characteristically fresher, with salinities ranging from 30 to 30.5, consistent with Alaska Coastal Water (ACW) characteristics and a shift to wind direction from the south (see Supplementary Figure 1 for a map depicting regional currents in the study area). Wind speed ranged from 0 to 36 kt during the course of the cruise as measured by a vessel-mounted anemometer corrected for ship motion. Throughout the water column, temperatures ranged from −0.5 to 10°C and salinities ranged from 30 to 33.5. Chlorophyll fluorescence ranged from 0 to 15 μg L–1 while nitrate concentration ranged from 0 to 20 μmol L–1 and was highest along the northernmost Cape Lisburne line. SUNA-derived nitrate concentration measurements were confirmed with bottle-derived nitrate concentration measurements. Sea ice was absent during the cruise. A salinity of 31 delimited stations A, B, and C into the Bering Shelf/Anadyr Waters (BSAW) (salinity > 31) and stations D, E, F, and G into ACW waters (salinity < 31) (Figure 2).

FIGURE 2

Drifting Sediment Trap Flux Rate Measurements

The depth of the water column varied little at our seven stations, ranging from 41 to 51 m and averaging 47 m, 17 m deeper than the drifting sediment traps sampling depths at 30 m (1% Photosynthetically Active Radiation (PAR)). The actual 1% PAR varied from between 16 and 30 m, except for station F having a 1% PAR reaching 38 m (Table 3). The overall average depth of the euphotic zone was 26 m, comparable to POC export measurements at 30 m. Sinking POC fluxes were high (up to 2.2 g C m–2 day–1) but spatially variable. Bering Shelf/Anadyr Waters were associated with higher fluxes (1.2–2.2 g C m–2 day–1) at stations A, B, and C, while lower fluxes (0.2–0.5 g C m–2 day–1) were characteristic of ACW at stations D, E, F, and G (Figure 2). For stations in the BSAW, sinking particles consisted mostly of aggregated diatoms and viable diatom cells while the ACW stations contained more diverse particles including fecal pellets, zooplankton, and diatom cells (Figure 3).

TABLE 3

Station NameWater MassBottom Depth (m)Euphotic Zone Depth (m)POC flux (g C m–2 day–1)Primary Productivity (g C m–2 day–1)Export RatioPN (g N m–2 day–1)δ13-C of sinking material (‰)δ15-N of sinking material (‰)
ABSAW41162.20 ± 0.194.240.52 ± 0.050.40 ± 0.04−21.08 ± 0.107.51 ± 0.49
BBSAW50241.18 ± 0.100.871.36 ± 0.120.21 ± 0.02−20.27 ± 0.066.23 ± 0.50
CBSAW50241.39 ± 0.072.150.65 ± 0.050.21 ± 0.01−19.94 ± 0.926.71 ± 0.29
DACW42260.48 ± 0.03**0.481.00 ± 0.070.47 ± 0.03**−19.55 ± 0.04**7.20 ± 1.43**
EACW49300.34 ± 0.030.740.45 ± 0.040.06 ± 0.01−21.64 ± 0.437.37 ± 0.17
FACW5138*0.34 ± 0.050.331.02 ± 0.140.05 ± 0.01−22.30 ± 0.155.63 ± 0.28
GACW46240.17 ± 0.000.230.75 ± 0.060.03 ± 0.00−24.03 ± 0.538.50 ± 1.12

Particulate organic carbon (POC) flux, primary productivity, and export ratios at seven stations in the Bering and Chukchi seas at drifting sediment trap sites.

Drifting sediment trap sampled at 30 m depth. BSAW, Bering Shelf/Anadyr Waters; ACW, Alaska Coastal Waters. *Euphotic zone depth deeper than 30 m, the depth of the drifting sediment trap deployment. **Outlier of 3 standard deviations higher than other 3 replicates removed from analysis Pacific Arctic’s extraordinary carbon fluxes.

FIGURE 3

The particulate nitrogen (PN) flux ranged from 0.03 to 0.47 g N m–2 day–1. Both the highest and lowest PN flux were measured in the ACW (Table 3). PN flux was slightly, though not significantly, higher in the BSAW (0.27 g N m–2 day–1) than in the ACW (0.15 g N m–2 day–1). Delta 13C ranged from −24.03 to −19.94‰ with no clear distinction in δ13C values between the ACW and BSAW. Both the least negative and most negative δ13C values were located at stations in the ACW. A similar pattern was true for δ15N values with values ranging from 5.63 to 8.5‰. There was no spatial pattern in δ15N values and both the highest and lowest δ15N values were found in stations in the ACW.

Respiration Rate Measurements

Overall carbon specific particle-associated microbial respiration ranged from −13.7 to 12.8% day–1 (Table 4). Negative carbon specific particle-associated microbial respiration indicates net respiration of carbon while positive values indicate net production of carbon. Five of the seven stations had carbon specific particle-associated microbial respiration that were indistinguishable from zero where carbon specific particle-associated microbial respiration could not be distinguished from free living microbial respiration. The two respiration measurements that were distinguishable from zero were from stations B (−13.7 ± 10.5% day–1) and C (12.8 ± 6.7% day–1), both of which are in BSAW. Therefore, all carbon specific particle-associated microbial respiration from stations in ACW were indistinguishable from zero. Average carbon specific particle-associated microbial respiration from BSAW stations (3.1 ± 14.6% day–1, average ± 1 standard deviation), ACW stations (1.7 ± 2.9% day–1), as well as all stations combined (2.3 ± 8.7% day–1) were indistinguishable from zero. Additionally, there was no statistical difference in carbon specific particle-associated microbial respiration between the ACW and BSAW.

TABLE 4

Station NameWater Massr_exp (μmol O2 m–1 L–1)r_control (μmol O2 m–1 L–1)r_PAM (μmol O2 m–1 L–1)R_PAM (% day–1)
ABSAW−0.215 ± 0.120−0.076 ± 0.086−0.139 ± 0.14710.3 ± 24.0
BBSAW−0.073 ± 0.070 0.017 ± 0.018−0.090 ± 0.072−13.7 ± 10.5
CBSAW0.087 ± 0.044−0.001 ± 0.0090.088 ± 0.04512.8 ± 6.7
DACW0.048 ± 0.0240.030 ± 0.0120.018 ± 0.027−0.9 ± 5.7
EACW0.019 ± 0.032−0.001 ± 0.0100.020 ± 0.0334.3 ± 8.0
FACW0.007 ± 0.041−0.009 ± 0.0080.016 ± 0.0424.2 ± 9.5
GACW−0.008 ± 0.032−0.003 ± 0.008−0.005 ± 0.033−0.8 ± 7.5

Particle-associated microbial respiration rates and carbon specific rates for sinking material.

Significant rates in bold and italicized. A total of 188 min were analyzed treating all 16 experimental sample as replicates and the final carbon value used to calculate the per carbon rate. r_exp, experimental respiration rate ± 1 SD; r_control, control respiration rate ± 1 SD; r_PAM, particle-associated microbial respiration rate ± 1 SD, R_PAM; carbon specific particle-associated microbial respiration rate ± 1 SD, BSAW; Bering Shelf/Anadyr Waters; ACW, Alaska Coastal Waters.

Primary Productivity Rate Measurements

Primary productivity rates were spatially variable with an overall range of 0.23–4.24 g C m–2 day–1. Station A had the highest rate of primary productivity and had the highest fluorescence signal of the seven stations (full depth CTD profiles are shown in Supplementary Figure 2). Here, a slight chlorophyll-a maximum of 10 mg m–3 fluorescence was present at 5 m depth and the water column was well-mixed. More pronounced, though lower chlorophyll-a maximums occurred at stations C (30 m) and F (38 m). Station F had the deepest 1% PAR depth of 38 m, co-occurring with the chlorophyll-a maximum, but had the second lowest rate of primary productivity of 0.34 g C m–2 day–1. Overall, primary productivity was higher in the BSAW (0.87–4.24 g C m–2 day–1) than in the ACW (0.23–0.74 g C m–2 day–1). All three stations in the BSAW had less stratified water columns, while all four stations in the ACW had more stratified water columns.

Export Ratios

Export ratios are a metric that characterizes the efficiency of the biological carbon pump, calculated using Eq. 2.

Higher export ratios indicate a more efficient biological carbon pump and lower ones less efficient. Export ratios ranged from 0.45 to 1.36 with the lowest export ratio observed at station E and the highest at station B (Table 3). There was no significant difference in export ratios between water masses (BSAW: 0.84 ± 0.45 and ACW: 0.81 ± 0.27 mean ± 1 SD). Three stations (B, D, and F) had export ratios at or above 1 and the overall study average export ratio was 0.82 ± 0.32 (mean ± 1 SD).

Moored Sediment Trap Time Series Flux Rate Measurements

Moored sediment trap-derived POC fluxes provide independent measures to compare with the fluxes observed with drifting sediment trap sampling. Drifting sediment trap sampling at stations B and C took place three and six days following the end of the moored sediment trap sampling, respectively. Particulate organic carbon fluxes of 0.8 and 1.2 g C m–2 day–1 at station B and 2.3 and 1.4 g C m–2 day–1 at station C were obtained with the moored and drifting sediment traps, respectively (Table 3 and Figure 4).

FIGURE 4

; , ; , ; ; , ; ) and black markers global context (). *No corresponding productivity data with (), , , and ().

At station B, POC fluxes were generally low (<0.25 g C m–2 day–1) from June through October 2017 with brief periods of elevated POC fluxes (0.6–1.3 g C m–2 day–1) occurring around the same time as peaks in fluorescence (Figure 5). Particulate organic carbon fluxes increased along with wind speed during November and December 2017 (0.9–1.4 g C m–2 day–1), decreased during January 2018, and remained relatively low (<0.65 g C m–2 day–1) when sea ice was present from January through late April 2018. Fluorescence remained low from mid-October 2017 through April 2018. Sea ice melted at the end of April 2018 and the highest POC fluxes were observed about one month later during late May 2018 (1.5 g C m–2 day–1).

FIGURE 5

). Wind velocity (black line) taken from modeled wind reanalysis (European Centre for Medium-Range Weather Forecasts, 2019). Twenty-four hour smoothing is shown with the thick black line. Lowpass-filtered current velocity (magenta line) taken from AD, with 12 h smoothing shown with the thick magenta line, temperature (red line), salinity (blue line), fluorescence (green line) taken from seabird SBE16plus unit coupled with a WetLabs fluorometer on each of the moorings.

At station C, high POC fluxes were recorded from June to mid-July 2017 (1.2–1.7 g C m–2 day–1), followed by a period of low POC fluxes from mid-July through mid-October 2017 (<0.5 g C m–2 day–1) (Figure 5). Spikes in fluorescence occurred sporadically from June until early October 2017. Particulate organic carbon fluxes increased starting in mid-October 2017 and were elevated throughout November and December 2017 (0.9–1.2 g C m–2 day–1). Particulate organic carbon flux dramatically decreased when sea ice formed during January 2018. A period of low POC fluxes was observed between January and May 2018 (<0.5 g C m–2 day–1), while sea ice was consistently present. The highest POC fluxes were measured at the beginning of June 2018 (2.3 g C m–2 day–1), at the same time as the highest peaks of fluorescence soon after sea ice retreated from this station.

The moored sediment trap time series indicates that POC flux had recently peaked at station B before drifting sediment trap sampling took place during early June 2018, but was likely at or near the period of peak annual flux at station C (Figure 5). The composition of the material collected in the moored sediment trap samples during June 2018 indicated the occurrence of a pelagic phytoplankton bloom at station B, reflected by the export of the exclusively pelagic centric diatoms Chaetoceros spp. and Thalassiosira spp. that usually dominate spring blooms on Arctic shelves (; ). Diatom fluxes, composed of several pennate and centric diatom groups, were nearly three times higher at station C than at station B, reflecting a large diatom bloom.

Discussion

The overall objectives of this study were to characterize the strength and efficiency of the biological carbon pump on the Pacific Arctic shelf during a warm, low-ice year in order to shed light on potential current and future changes in carbon cycling in this region. We addressed this by considering three major aspects of relevance: primary productivity, sinking POC flux, and particle-associated microbial respiration during June of 2018.

Regional Spatial Trends

The largest distinction of regional spatial trends occurred between the two water masses present in this region (ACW and BSAW). Consistent with previous studies, the ACW was warmer and fresher with lower nutrients than the BSAW during June 2018 (; ). Both regions were warmer than normal for this time of year and had experienced much less sea ice than normal (; ).

We measured consistently lower POC flux rates in ACW than BSAW. These results support the previously untested hypothesis that POC fluxes would be higher in the BSAW compared to ACW (). This distinction in primary productivity and POC flux between the ACW and BSAW could partially be attributed to differences in stratification between these regions; we found a less stratified water column in the BSAW, which could contribute to nutrient input to the surface and allow for higher primary productivity rates. Well-mixed water could also help facilitate POC mixing out of the euphotic zone, however it is more likely that this material could be brought back up to the surface through the same mechanism. The differences in primary productivity and POC flux cannot be fully explained by differences in water column stratification. The two regions were characterized by different types of sinking particles: more uniform aggregated diatoms and viable diatoms were found in the BSAW, while more processed material like fecal pellets and zooplankton were found in the ACW in addition to living diatoms. This distinction suggests more processing of POC by zooplankton or heterotrophic bacteria in the ACW. However, there was no difference in export ratio or particle-associated microbial respiration between the two water masses, indicating that POC fluxes are mostly regulated by primary production rates rather than heterotrophic processing. The spatially uniform particle associated microbial respiration rates we measured do not support the postulation that the ACW would have higher particle-associated microbial respiration rates (; ). Nonetheless, the higher primary productivity and POC fluxes in the BSAW demonstrate a stronger biological carbon pump in the BSAW region.

PN, δ13C, and δ15N values were not significantly different between the BSAW and the ACW. Higher PN in sinking material are associated with more nutritious food for the benthos (). Less negative δ13C values tend to indicate a larger influence of ice algae or a marine signature while more negative δ13C values tend to indicate a more coastal or terrigenous signature (). Larger δ15N values were associated with material that is higher on the food chain (i.e., secondary and tertiary producers), while lower δ15N values are associated with material lower on the food chain (i.e., primary producers) ().

The ACW had lower primary productivity than the BSAW. This regional pattern has been well described previously (; ) and is attributed to the lower nutrient concentrations in ACW compared with BSAW (). Primary productivity values in both water masses fell within a typical range for these regions (; ; ). However, a primary productivity rate of 16 g C m–2 day–1 has been previously observed in the Pacific Arctic Shelf (; ), much higher than what we measured and what is typically measured. If these higher production rates were associated with export ratios similar to what we observed here, then the associated fluxes would be even more remarkable than the values we observed with the DSTs and MSTs during this study.

Stations with higher rates of primary productivity tended to have higher rates of POC flux. However, there was not a perfect relationship between primary productivity and POC flux, which caused some variations in the export ratios. The stations with export ratios over 1 and the high average export ratio indicate an extremely efficient biological carbon pump or temporal or spatial decoupling between primary production and flux.

While the dominant regional patterns were associated with water masses, we also expected some patterns falling along a latitudinal gradient. It is difficult to separate the signal of water mass from latitude because most of the stations that were classified as BSAW were located south of the stations classified as ACW. We found higher daily primary productivity rates and POC fluxes at the southern stations in the BSAW than at the northern stations in the ACW. The annual POC flux was higher at station C (215 g C m–2 year–1) than at station B (204 g C m–2 year–1), indicating an increase in POC flux with latitude. However, an annual POC flux lower than these (145 g C m–2 year–1) was measured at about 200 miles north of our study area (). Drifting sediment trap sampling at more stations in the southern portion of the ACW and moored sediment trap sampling in the ACW are needed to better tease apart the differences associated with latitude and water mass in this region.

The annual pattern of POC flux shows some latitudinal distinction between the more northern station C and the more southern station B. The peak annual flux occurred during early June 2018 at station C, while the peak annual flux occurred a couple of weeks earlier at station B (late May 2018). The spring peak flux was higher at station C (2.3 g C m–2 day–1) compared to station B (1.5 g C m–2 day–1). Increased POC flux measurements occurred in the absence of peaks in fluorescence at both stations from November 2017 to January 2018 strongly suggesting episodic resuspension events during fall. Particulate organic carbon fluxes decreased in the presence of sea ice, reducing wind mixing and resuspension. This is particularly evident at station C under higher sea ice concentrations. It is likely that these fall high flux events do not represent increased net flux, as they are likely partly the result of material that previously fell to the seafloor being resuspended and collected in the moored sediment trap again.

Role of Heterotrophy in the Water Column

Bacterial production largely controls how much exported POC reaches the seafloor and might increase in Arctic waters under more acidic, warmer, and lower-ice conditions (; ), which could result in higher pelagic community carbon demand (). One factor that has received a lot of thought for the Pacific Arctic shelf region is how a reduction in ice algae relative to to pelagic phytoplankton as primary producers might impact the benthic-pelagic coupling, with the prediction that smaller pelagic phytoplankton will have slower sinking rates and will therefore be more likely to be consumed by zooplankton or bacteria in the water column, resulting in less material reaching the seafloor (; ; ). A slower particle sinking rate will give more time for bacteria to both colonize and degrade sinking particles. In our study, we mostly collected pelagic diatoms in the drifting sediment traps rather than species associated with ice. However, our direct measurements of microbial respiration rates associated with sinking particles were mostly indistinguishable from zero. This is not unprecedented in high latitude regions ().

Conducting a comparison of measured particle-associated microbial respiration rates from around the globe, we found that particle-associated microbial respiration generally decreases with increasing latitude (Figure 6). In our study, the rate of particle-associated microbial respiration was 2% day–1 on average, with a 95% confidence interval ranging from −24.2 to 34.3% day–1. Given the shallow nature of the Pacific Arctic shelf (20 m average distance from base of euphotic zone to seafloor) and rapid sinking velocity of material caught in the traps (greater than 100 m day–1), even under the fastest respiration rate (−24.2% day–1) we calculated that less than 5% of the exported organic carbon would be remineralized within the water column before being deposited on the seafloor. Considering our conservative estimates, the true consumption is likely much smaller than this value. We conclude that particle-associated microbial respiration does not play a large role in recycling POC below the euphotic zone in this region, implying that most of the material that is exported from the euphotic zone will likely reach the shallow seafloor.

FIGURE 6

2, 3, 4, 5, 6, 7, and 8.

One mitigating factor in how much organic matter is deposited on the seafloor is the role zooplankton and free-living heterotrophic microbes play in consuming organic matter in the water column. Historically, zooplankton have not consumed large proportions of organic matter in the water column (; ; ; ; ), but it is possible they may play a larger role in the future. One study from just north of our study area found that primary production rates were similar to free-living community microbial respiration rates during the summer, indicating a large proportion of primary production could be consumed by free-living microbes (). High export ratios in the present study point to zooplankton and free-living heterotrophic microbes playing a small role in consuming organic matter within the euphotic zone. It is possible that zooplankton and bacteria may play a larger role in consuming POC later in the summer. We suggest measuring export ratios and particle associated microbial respiration rates in August on the Pacific Arctic shelf to answer this remaining question.

Comparison of Drifting and Moored Sediment Trap POC Fluxes

The POC flux measurements measured with the drifting and moored traps at stations B and C, while of similar magnitude, were not the same. Many factors potentially caused variations between POC flux values obtained with drifting and moored traps. One reason is the different sampling times as POC fluxes may have changed on time scales much shorter than three or six days. In addition, the moored sediment trap measured flux over eight days, while the drifting sediment traps measured flux for six and a half hours and station B and five and a half hours at station C. If there is a diurnal cycle in flux regulated by zooplankton or phytoplankton it can be captured in the drifting sediment trap sampling and masked in the moored sediment trap sampling. Our study took place on a shallow Arctic shelf over the summer solstice. It is unlikely there was a diurnal cycle of primary production due the nearly 24 h of sunlight that were present. Additionally, due to the shallowness of the shelf, zooplankton in this region are not know to exhibit diel vertical migration (; ). It is unlikely that any differences in POC flux are a result of any changes in flux as a diurnal cycle, and, if present, are likely due to changes in the rate of primary production controlled by variable cloud cover and nutrient availability. Finally, the locations of the moored sediment trap and drifting sediment trap sampling did not perfectly overlap, although the drifting sediment traps were deployed within half a mile of the moored sediment traps.

Even if the sampling of these two traps perfectly overlapped in time and space, it is unlikely that they would produce the same POC flux values. One reason is because the moored sediment traps sample with a Eulerian framework, being moored in one location sampling various water masses as they flow, while drifting sediment traps sample with a Lagrangian approach, staying with one parcel of water and sampling it continuously as it moves with the currents. Another reason is that each of these trap designs have their own individual biases. In high current environments moored sediment traps can tilt to the side, affecting the collection of sinking particles. However, no tilt occurred at stations B and C based on CTD data. In contrast, drifting sediment traps may reduce the vertical shear in high current environments by floating freely within the water column. Additionally, we minimized other hydrodynamic concerns by using a bungee to dampen surface motion, tubes with a high aspect ratio, and bottom weighted tubes to keep them upright (; , ). However, sinking POC flux may be incorrectly measured with drifting sediment traps because the tubes are open during deployment and recovery, contrary to the moored sediment traps (). This source of error is minimized by deploying the drifting sediment traps at a shallow depth (i.e., 30 m).

With these sources of error, it is helpful to have two independent measurements of POC flux using different methods. The overall range of POC flux values in June 2018 was 1.48 to 2.29 g C m–2 day–1 with the moored sediment traps and 0.17 to 2.20 g C m–2 day–1 with the drifting sediment traps. Comparable maximum flux magnitudes from these two different methods minimize the concerns of collection biases common with sediment traps and provide some supporting evidence of the validity of POC fluxes of this magnitude.

Comparing POC Flux, Primary Productivity, and Export Ratios

Primary productivity rates, sinking flux, and export ratios were compared with previous measurements from the same study area (; ; , ; ), the broader Bering and Chukchi shelf system (; ; ; ), the greater Arctic area (Supplementary Table 1), and from a global compilation (; Figure 4, and Supplementary Figure 2). The upper range of our POC flux measurements (2.2 g C m–2 day–1 from drifting sediment trap and 2.3 g C m–2 day–1 from moored sediment trap) was unprecedentedly high compared to other measurements in this region and among the highest recorded in the surrounding areas, the broader Arctic, and globally.

Five previous studies report particulate flux estimates for the Bering and Chukchi shelves (; ; , ; ), two based on sediment trap measurements. A moored sediment trap deployed (36 m water depth, 49 m bottom depth) from late June to late September 1988, about 500 miles south of Bering Strait, measured flux ranging from 253 to 654 mg C m–2 day–1 (). More recently, a moored sediment trap deployed (37 m water depth, 45 m bottom depth) from August 2015 to July 2016, about 200 miles north of our most northern stations on the Chukchi shelf, measured POC fluxes ranging from 72 to 1184 mg C m–2 day–1 (). An estimate of 456 mg C m–2 day–1 (36 m water depth, 49 m bottom depth) was calculated using the 234Th/238U disequilibrium method on the Chukchi Sea shelf in August 1994 (). measured a POC flux value of 243.8 mg C m–2 day–1 (40 m water depth, 50 m bottom depth) using the 234Th/238U disequilibrium method sometime between July and September on the Chukchi shelf within the bounds of our study area. Finally, during a study on the Chukchi shelf from July to September, a POC flux measurement of 951.1 mg C m–2 day–1 (30 m water depth, 40 m bottom depth) was made using the 234Th/238U disequilibrium method within the bounds of and to the north of our study area (). We selected the peak annual POC flux value from and the flux measurement from and , and plotted them directly on the y-axis indicating no known corresponding primary productivity value (Figure 4). The spatial extent of these studies are shown along with that of this study in Supplementary Figure 3. It should be noted that the 234Th/238U disequilibrium method for calculating sinking POC flux has its own biases, especially in areas with non-steady state flux events and advection and dispersion processes (), such as on a shallow Arctic shelf. The measurements of POC flux we made in this study were the same or higher than previous measurements made in this region.

We expanded our region of comparison to include the Bering and Chukchi shelf breaks. Particulate organic carbon flux has been estimated more frequently on the Bering and Chukchi shelf breaks, with a maximum flux value at the base of the euphotic zone of 1.381 g C m–2 day–1 reported slightly south of our study area on the Bering Sea shelf break (40 m water depth, >125 m bottom depth) in July 2010 (Figure 6; ). Additionally, our average regional flux for the BSAW, 1.59 ± 0.54 g C m–2 day–1 (mean ± 1 SD) is much higher than previous average regional flux estimates from the shelf breaks just north and south of this region, which range from 34 mg C m–2 day–1 (50 m water depth, bottom depth average 1275 m, May and June sampling period) to 376 mg C m–2 day–1 (50 m water depth, bottom depth average 838 m, May and June sampling period) (, ; ; ; ). Although individual measurements of export ratios approaching 1 are somewhat common (; ), our average export ratio of 0.82 ± 0.32 (mean ± 1 SD) is very high. These observations illustrate the exceptional efficiency and strength of the biological carbon pump in the shallow Pacific Arctic shelf and shelf breaks.

Sinking POC flux and primary productivity values were also compared with a global review of POC flux measurements obtained using the 234Th/238U disequilibrium method (). We plotted all values with both POC flux and primary productivity rate measurements from , along with previous measurements from the broader Pacific shelf system, and our specific study area with the values we measured in this study from the drifting sediment traps and the last values of POC flux from the moored sediment traps (Figure 4). When available, sinking POC flux measurements at the base of the euphotic zone were selected, a metric shown to be comparable at sites with different bottom and euphotic zone depths (). From this analysis, it is evident that the primary productivity rates from this study mostly fall within the upper range of what has been measured in this region before. Flux measurements at the BSAW stations (0.8–2.3 g C m–2 day–1) are very high compared to what has been measured previously (0–1.4 g C m–2 day–1).

In addition to the global review by , we compiled POC fluxes, primary productivity rates, and export ratios from other high latitude studies. The results of this review can be found in Supplementary Table 1. We categorized the studies by region, including the Baffin Bay, Baltic Sea, Barents Sea, Beaufort Sea, Bering Sea, Canadian Archipelago, Chukchi Sea, Fram Strait, Greenland Sea, Hudson Bay, Kara Sea, Labrador Sea, Laptev Sea, North Atlantic, North Sea, White Sea, as well as the high Arctic. We considered studies that measured sinking POC flux rates using drifting, moored, or neutrally buoyant sediment traps, marine snow catchers, or 234Th/238U disequilibrium. When possible, we selected samples from as close to the euphotic zone as possible. Primary productivity rates ranged from 0 to 2.6 g C m–2 day–1 and export ratios ranged from 0.03 to 1.67. These ranges are consistent with the primary productivity and export ratios measured in this study. Out of 79 studies, only nine measured rates of POC flux greater than 1 g C m–2 day–1 that were measured in the Baffin Bay (), the Barents Sea (; ; ; ), the Bering Sea (), the Beaufort Sea (), the Chukchi Sea (), and the North Atlantic (). The highest POC flux measurement of 2.5 g C m–2 day–1 was measured in the Beaufort Sea near the Makenzie River drainage in June at 50 m depth and 230 m water depth (). The second highest POC flux measurement of 1.5 g C m–2 day–1 was measured in the Barents Sea in May at 30 m depth and 239 m bottom depth (). We also selected some studies of POC flux from known high productivity and/or high flux regions from around the world (Supplementary Table 2). These POC flux values ranged from 1 to 620 mg C m–2 day–1. Two POC flux values recorded in the present study were among the highest ever recorded at the base of the euphotic zone (2.20 g C m–2 day–1 at station A with a drifting sediment trap and 2.29 g C m–2 day–1 at station C with a moored sediment trap). With individual flux values from different methods at different stations ranking among the highest ever recorded, it is clear that the Pacific Arctic shelf exported a massive amount of organic carbon out of the euphotic zone, even in an anomalously warm year with low sea ice.

Limitations and Implications

Even with strong efforts in place to study the processes on the Pacific Arctic Shelf with the ASGARD program, there is still a lack of available data. While POC flux and water column oceanographic measurements are being obtained more frequently with moored sediment traps, temporally overlapping primary productivity rate measurements are not often available, notably later in the summer or earlier in the spring, when production is highest. Previous primary productivity measurements obtained later in the summer (; ; ) may no longer be representative of current conditions. Additionally, particle-associated microbial respiration rates are likely variable throughout the spring, summer and fall and therefore cannot be extrapolated beyond spring. Therefore, there are still many unknowns regarding how the Pacific Arctic is responding and will respond to climate change. It is possible that only after many years of consistently warm and low-ice conditions changes in the strength and efficiency of the biological carbon pump on the Pacific Arctic shelf will become apparent.

It has been hypothesized that the strength of the biological carbon pump in the Arctic may weaken with climate change due to increased duration of the open-water period for primary production and enhanced nutrient limitation (; ; ). Warmer waters have been predicted to increase metabolic rates of pelagic grazers and heterotrophic bacteria and potentially favor smaller phytoplankton and faster-growing grazers that more rapidly recycle organic matter within the water column (; ). We postulate that high nutrient concentrations, the shallow nature of the Pacific Arctic shelf, and the large-celled, fast-sinking phytoplankton that dominate pelagic productivity create conditions unique to this Arctic shelf (; ). Nutrients are unlikely to become more limited, especially in the BSAW, because of the consistent influx of the Anadyr Current, which is nutrient replete from Pacific upwelling. High nutrient concentrations favor large-cell phytoplankton (). Therefore, it is possible the Pacific Arctic shelf will not experience as dramatic of a shift from large cells to small cells with warming conditions when compared to other Arctic shelves. Finally, the Pacific Arctic shelf is shallower than most other Arctic shelves (averaging only 50 m deep). It will never take very long for sinking material leaving the euphotic zone to reach the seafloor, as it only has to sink about 20 m. Even if there is some increased heterotrophy in zooplankton or bacteria or decrease in cell size, the shallow nature of the shelf will allow for a higher proportion of organic matter to reach the seafloor than over deeper shelves, such as the Canadian Arctic Archipelago and European Arctic shelf. Increased frequency of storms are predicted with a changing Arctic (). This could have major implications for the long-term carbon storage that historically has occurred on the Pacific Arctic shelf. found that significant portions of carbon once thought to be stored in sediments on the Chukchi shelf are mixed up during fall storms. Given these features, we speculate that this system may retain strong coupling between the pelagic and benthic realms, continue to support highly productive pelagic and benthic ecosystems, and act as a strong sink for atmospheric carbon dioxide, possibly mediated by increased frequency of fall and winter storms. If these results prove to be a sustained feature of the rapidly changing Pacific Arctic, the biological carbon pump could represent an important element of resilience for regional ecosystems and biogeochemical cycles.

Conclusion

Measurements from both drifting and moored sediment traps indicate that fluxes of sinking POC on the Pacific Arctic shelf in June 2018 ranged from 0.8 to 2.3 g C m–2 day–1 in BSAW, making them amongst the highest fluxes ever documented in the global oceans. This region was also characterized by high export ratios and low rates of particle associated microbial respiration. These observations indicate that the biological carbon pump on the Pacific Arctic shelf is exceptionally strong and efficient despite a recent multi-year shift to warmer and relatively ice-free conditions (; ; ; ). While the majority of the fluxes we observed during June 2018 were unprecedented relative to the limited number of historical flux measurements from this region, the data are still insufficient to determine whether functioning of the biological carbon pump has changed significantly relative to earlier, colder, and ice-replete years. Nonetheless, our observations do not provide supporting evidence for the common prediction that a weaker biological carbon pump and increased pelagic heterotrophy will prevail on the Pacific Arctic’s continental shelves under future change.

Statements

Data availability statement

The datasets generated for this study section can be found in the DataONE system (10.24431/rw1k46v).

Author contributions

SO, SD, SH, RH, DS, and AM: conceptualization. SD, SH, RH, DS, and AM: funding acquisition and project administration. SO, RH, CL, DS, and AM: methodology. SO, SD, CL, DS, and AM: investigation and formal analysis. SD, SH, and AM: supervision. SO and RH: software. SO: data visualization, writing-original draft. SO, SD, SH, RH, CL, DS, and AM: writing-reviewing and editing.

Funding

Research funded by the NPRB (Grant A91-00a, A91-88 and A91-99a to SD, SH, RH, AM, and DS).

Acknowledgments

We thank Captain and crew of the R/V Sikuliaq as well as all collaborators in the ASGARD project, for making this project possible. Additionally, we thank S. B. Moran for the use of the drifting sediment traps, the Alaska Stable Isotope Facility for technical support and analytical assistance (University of Alaska, Fairbanks), and the two reviewers for their constructive comments. NPRB Publication Number “ArcticIERP-40”.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2020.548931/full#supplementary-material

Abbreviations

  • ACW

    Alaska Coastal Waters

  • ASGARD

    Arctic Shelf Growth Advection, Respiration and Deposition rate experiments project

  • BSAW

    Bering Shelf/Anadyr Waters

  • PN

    particulate nitrogen

  • POC

    particulate organic carbon.

References

  • 1

    AmielD.CochranJ. K. (2008). Terrestrial and marine POC fluxes derived from 234Th distributions and δ13C measurements on the Mackenzie Shelf.J. Geophys. Res.113:C03S06. 10.1029/2007JC004260

  • 2

    AndersenP. (1988). The quantitative importance of the “Microbial loop” in the marine pelagic: a case study for the North Bering/Chukchi Sea.Arch. Hydrobiol. Bieh.31243251.

  • 3

    AndersonL. A.SarmientoJ. L. (1994). Redfield ratios of remineralization determined by nutrient data-analysis.Global Biogeochem. Cycles86580. 10.1029/93GB03318

  • 4

    AndreassenI. J.WassmannP. (1998). Vertical flux of phytoplankton and particulate biogenic matter in the marginal ice zone of the Barents Sea in May 1993.Mar. Ecol. Prog. Ser.170114. 10.3354/meps170001

  • 5

    ArrigoK. R.PerovichD. K.PickartR. S.BrownZ. W.van DijkenG. L.LowryK. E.et al (2014). Phytoplankton blooms beneath the sea ice in the Chukchi sea.Deep Sea Res. II105116. 10.1016/j.dsr2.2014.03.018

  • 6

    AshjianC. J.CampbellR. C.WelchH. E.ButlerM.Van KeurenD. (2003). Annual cycle in abundance, distribution, and size in relation to hydrography of important copepod species in the western Arctic Ocean.Deep Sea Res. I5012351261. 10.1016/S0967-0637(03)00129-8

  • 7

    BatesN. R. (2006). Air-sea CO2 fluxes and the continental shelf pump of carbon in the Chukchi Sea adjacent to the Arctic Ocean.J. Geophys. Res. Ocean111121. 10.1029/2005JC003083

  • 8

    BaumannM. S.MoranS. B.LomasM. W.KellyR. P.BellD. W. (2013). Seasonal decoupling of particulate organic carbon export and net primary production in relation to sea-ice at the shelf break of the eastern Bering Sea: Implications for off-shelf carbon export.J. Geophys. Res. Ocean11855045522. 10.1002/jgrc.20366

  • 9

    BelcherA.IversenM.GieringS.RiouV.HensonS. A.BerlineL.et al (2016a). Depth-resolved particle-associated microbial respiration in the northeast Atlantic.Biogeosciences1349274943. 10.5194/bg-13-4927-2016

  • 10

    BelcherA.IversenM.MannoC.HensonS. A.TarlingG. A.SandersR. (2016b). The role of particle associated microbes in remineralization of fecal pellets in the upper mesopelagic of the Scotia Sea, Antarctica.Limnol. Oceanogr.6110491064. 10.1002/lno.10269

  • 11

    BluhmB. A.GradingerR. R. (2008). Regional variability in food availability.Ecol. Appl.18S77S96. 10.1890/06-0562.1

  • 12

    BuesselerK. O.BaconM. P.LivingstonH. D.CochrantK. (1992). Carbon and nitrogen export during the JGOFS North Atlantic Bloom experiment estimated from 234Th: 238U disequilibria.Deep. Sea Res.391115–1137. 10.1016/0198-0149(92)90060-7

  • 13

    BuesselerK. O.BoydP. W.BlackE. E.SiegelD. A. (2020). Metrics that matter for assessing the ocean biological carbon pump.Proc. Natl. Acad. Sci. U.S.A.11796799687. 10.1071/pnas.1918114117

  • 14

    BuesselerK. O.ChenM.HaradaK.GustafssonO.TrullT.Rutgers van der LoeffM.et al (2007). An assessment of the use of sediment traps for estimating upper ocean particle fluxes.J. Mar. Res.65345416. 10.1357/002224007781567621

  • 15

    ButmanC. A.GrandW. D.StolzenbackK. D. (1986). Predictions of sediment trap biases in turbulent flows, a theoretical analysis based on observations from the literature.J. Mar. Res.44601644. 10.1357/002224086788403024

  • 16

    CampbellR. G.SherrE. B.AshjianC. J.PlourdeS.SherrB. F.HillV.et al (2009). Mesozooplankotn prey preference and grazing impact in the western Arctic Ocean.Deep Sea Res. II5612741289. 10.1016/j.dsr2.2008.10.027

  • 17

    CarrollM. L.CarrollJ. (2003). “The arctic seas,” in Biogeochemistry of Marine Systems, edsBlackK. D.ShimmieldG. B. (Boca Raton, FL: CRC Press), 127156.

  • 18

    ChenC. A.BorgesA. V. (2009). Reconciling opposing views on carbon cycling in the coastal ocean: continental shelves as sinks and near-shore ecosystems as sources of atmospheric CO2.Deep Sea Res. II56578590. 10.1016/j.dsr2.2008.12.009

  • 19

    CollinsJ. R.EdwardsB. R.ThamatrakolnK.OssolinskiJ. E.DitullioG. R.BidleK. D.et al (2015). The multiple fates of sinking particles in the North Atlantic Ocean.Glob. Biogeochem. Cycles2914711494. 10.1002/2014GB005037

  • 20

    CotaG. F.PomeroyL. R.HarrisonW. G.JonesE. P.PetersF.SheldonW. M.et al (1996). Nutrients, primary production and microbial heterotrophy in the southeastern Chukchi Sea: arctic summer nutrient depletion and heterotrophy.Mar. Ecol. Prog. Ser.135247258. 10.3354/meps135247

  • 21

    DanielsonS. L.AhkingaO.AshjianC.BasyukE.CooperL. W.EisnerL.et al (2020). Manifestation and consequences of warming and latered heat fluxes over the Bering and Chukchi Sea continental shelves.Deep Sea Res. II177:104781. 10.1016/j.dsr2.2020/104781

  • 22

    DanielsonS. L.EisnerL.LaddC.MordyC.SousaL.WeingartnerT. J. (2017). A comparison between late summer 2012 and 2013 water masses, macronutrients, and phytoplankton standing crops in the northern Bering and Chukchi Seas.Deep Sea Res. II135726. 10.1016/j.dsr2.2016.05.024

  • 23

    DegerlundM.EilertsenH. C. (2010). Main species characteristics of phytoplankton spring blooms in NE Atlantic and Arctic waters (68-80° N).Estuar. Coasts33242269. 10.1007/s12237-009-9167-7

  • 24

    DugdaleR. C.GoeringJ. J. (1967). Uptake of new and regenerated forms of nitrogen in primary productivity.Limnol. Oceanogr.12196206. 10.4319/lo.1967.12.2.0196

  • 25

    DurkinC. A.EstapaM. L.BuesselerK. O. (2015). Observations of carbon export by small sinking particles in the upper mesopelagic.Mar. Chem.1757281. 10.1016/j.marchem.2015.02.011

  • 26

    EbersbachF.TrullT. W. (2008). Sinking particle properties from polyacrylamide gels during the KErguelen Ocean and Plateau compared Study (KEOPS): zooplankton control of carbon export in an area of persistent natural iron inputs in the Southern Ocean.Limnol. Oceanogr.53212224. 10.4319/lo.2008.53.1.0212

  • 27

    European Centre for Medium-Range Weather Forecasts (2019). Updated Monthly. ERA5 Reanalysis (0.25 Degree Latitude-Longitude Grid). Research Data Archive at the National Center for Atmospheric Research, Computational and Information Systems Laboratory.Reading: European Centre for Medium-Range Weather Forecasts.

  • 28

    FettererF.KnowlesK.MeierW. N.SavoieM.WindnagelA. K. (2017). Updated Daily. Sea Ice Index, Version 3.Boulder, CO: National Snow and Ice Data Center.

  • 29

    FukuchiM.SasakiH.HattoriH.MatsudaO.TanimuraA.HandaN.et al (1993). Temporal variability of particulate flux in the northern Bering Sea.Cont. Shelf Res.13693704. 10.1016/0278-4343(93)90100-C

  • 30

    GarneauM. ÈVincentW. F.TerradoR.LovejoyC. (2009). Importance of particle-associated bacterial heterotrophy in a coastal Arctic ecosystem.J. Mar. Syst.75185197. 10.1016/j.jmarsys.2008.09.002

  • 31

    GradingerR. (1999). Vertical fine structure of the biomass and composition of algal communities in Arctic pack ice.Mar. Biol.133745754. 10.1007/s002270050516

  • 32

    GradingerR. (2009). Sea-ice algae: major contributors to primary production and algal biomass in the Chukchi and Beaufort Seas during May/June 2002.Deep Sea Res. II5612011212. 10.1016/j.dsr2.2008.10.016

  • 33

    GrebmeierJ. M. (2012). Shifting patterns of life in the Pacific Arctic and Sub-Arctic seas.Ann. Rev. Mar. Sci.46378. 10.1146/annurev-marine-120710-100926

  • 34

    GrebmeierJ. M.BarryJ. P. (1991). The influence of oceanographic processes on pelagic-benthic coupling in polar regions?: a benthic perspective.J. Mar. Syst.2495518. 10.1016/0924-7963(91)90049-z

  • 35

    GrebmeierJ. M.McRoyC. P. (1988). Pelagic-benthic coupling on the shelf of the northern Bering and Chukchi Seas. I. Food supply and source and benthic biomass.Mar. Ecol. Prog. Ser.485767. 10.3354/meps053079

  • 36

    GrebmeierJ. M.McRoyC. P. (1989). Pelagic-benthic coupling on the shelf of the northern Bering and Chukchi Seas. III. Benthic food supply and carbon cycling.Mar. Ecol. Prog. Ser.5393100.

  • 37

    GrebmeierJ. M.McRoyC. P.FederH. M. (1988). Pelagic-benthic coupling on the shelf of the northern Bering and Chukchi Seas. I. Food supply and source and benthic biomass.Mar. Ecol. Prog. Ser.485767. 10.3354/meps053079

  • 38

    GustafssonÖLarssonU.AnderssonP.GeltingJ.RoosP. (2013). An assessment of upper ocean carbon and nitrogen export fluxes on the boreal continental shelf: a 3-year study in the open Baltic Sea comparing sediment traps, 234Th proxy, nutrient, and oxygen budgets.Limnol. Oceanogr. Methods11495510. 10.4319/lom.2013.11.495

  • 39

    HamaT.MiyazakiT.OgawaY.IwakumaT.TakahashiM.OtsukiA.et al (1983). Measurement of photosynthetic production of a marine phytoplankton population using a stable 13C isotope.Mar. Biol.733136. 10.1007/BF00396282

  • 40

    HauriC.WinsorP.JuranekL. W.McDonnellA. M. P.TakahashiT.MathisJ. T. (2013). Wind-driven mixing causes a reduction in the strength of the continental shelf carbon pump in the Chukchi Sea.Geophys. Res. Lett.4059325936. 10.1002/2013GL058267

  • 41

    HillV.ArdynaM.LeeS. H.VarelaD. E. (2018). Decadal trends in phytoplankton production in the Pacific Arctic region from 1950 to 2012.Deep Sea Res. II1528294. 10.1016/j.dsr2.2016.12.015

  • 42

    HopcroftR. R.KosobokovaK. N.PinchukA. I. (2010). Zooplankton community patterns in the Chukchi Sea during summer 2004.Deep Sea Res. II572739. 10.1016/j.dsri.2009.08.003

  • 43

    HovelsrudG. K.MckennaM.HuntingtonH. P. (2008). Marine mammal harvests and other interactions with humans.Ecol. Appl.18135147. 10.1890/06-0843.1

  • 44

    HuntingtonH. P.DanielsonS. L.WieseF. K.BakerM.BovengP.CittaJ. J.et al (2020). Evidence suggests potential transformation of the Pacific Arctic ecosystem is underway.Nat. Clim. Chang.10342348. 10.1038/s41558-020-0695-2

  • 45

    IversenM. H.PlougH. (2010). Ballast minerals and the sinking carbon flux in the ocean: carbon-specific respiration rates and sinking velocity of marine snow aggregates.Biogeosciences726132624. 10.5194/bg-7-2613-2010

  • 46

    KitamuraM.AmakasuK.KikuchiT.NishinoS. (2017). Seasonal dynamics of zooplankton in the southern Chukchi Sea revealed from acoustic backscattering strength.Cont. Shelf Res.1334758. 10.1016/j.csr.2016.12.009

  • 47

    LalandeC.GrebmeierJ. M.HopcroftR. R.DanielsonS. L. (2020). Annual cycle of export fluxes of biogenic matter near Hanna Shoal in the northeast Chukchi Sea.Deep Sea Res. II177:104730. 10.1016/j.dsr2.2020.104730

  • 48

    LalandeC.LeporeK.CooperL. W.GrebmeierJ. M.MoranS. B. (2007). Export fluxes of particulate organic carbon in the Chukchi Sea: a comparative study using 234Th/238U disequilibria and drifting sediment traps.Mar. Chem.103185196. 10.1016/j.marchem.2006.07.004

  • 49

    LalandeC.MoranS. B.WassmannP.GrebmeierJ. M.CooperL. W. (2008). 234Th-derived particulate organic carbon fluxes in the northern Barents Sea with comparison to drifting sediment trap fluxes.J. Mar. Syst.73103113. 10.1016/j.marsys.2007.09.004

  • 50

    LalandeC.NöthigE. M.FortierL. (2019). Algal export in the arctic ocean in times of global warming.Geophys. Res. Lett.4659595967. 10.1029/2019GL083167

  • 51

    Le MoigneF. A. C.HensonS. A.SandersR. J.MadsenE. (2013). Global database of surface ocean particulate organic carbon export fluxes diagnosed from the 234Th technique.Earth Syst. Sci. Data5295304. 10.5194/essd-5-295-2013

  • 52

    LeeS. H.JooH. M.LiuZ. L.ChenJ. F.HeJ. F. (2012). Phytoplankton productivity in the newly opened waters of the Western Arctic Ocean.Deep Sea Res II81-841827. 10.1016/j.dsr2.2011.06.005

  • 53

    LeeS. H.WhitledgeT. E.KangS. H. (2007). Recent carbon and nitrogen uptake rates of phytoplankton in Bering Strait and the Chukchi Sea.Cont. Shelf Res.2722312249. 10.1016/j.csr.2007.05.009

  • 54

    LeporeK.MoranS. B.GrebmeierJ. M.CooperL. W.LalandeC.MaslowskiW.et al (2007). Seasonal and interannual changes in particulate organic carbon export and deposition in the Chukchi Sea.J. Geophys. Res. Ocean112114. 10.1029/2006JC003555

  • 55

    LiW. K.McLaughlinF. A.LovejoyC.CarmackE. C. (2009). Smallest algae thrive as the arctic ocean freshens.Science326:539. 10.1126/science.1179798

  • 56

    McDonnellA. M. P.BuesselerK. O. (2010). Variability in the average sinking velocity of marine particles.Limnol. Oceanogr.5520852096. 10.4319/lo.2010.55.5.2085

  • 57

    McDonnellA. M. P.BoydP. W.BuesselerK. O. (2015). Effects of sinking velocities and microbial respiration rates on the attenuation of particulate carbon fluxes through the mesopelagic zone.Glob. Biogeochem. Cycles29175193. 10.1002/2014GB004935

  • 58

    MichelC.NielsenT. G.NozaisC.GosselinM. (2002). Significance of sedimentation and grazing by ice micro- and meiofauna for carbon cycling in annual sea ice (northern Baffin Bay).Aquat. Microb. Ecol.305668. 10.3354/ame030057

  • 59

    MooreS. E.KuletzK. J. (2019). Marine birds and mammals as ecosystem sentinels in and near distributed biological observatory regions: an abbreviated review of published accounts and recommendations for integration to ocean observatories.Deep Sea Res. II162211217. 10.1016/j.dsr2.2018.09.004

  • 60

    MooreS. E.StabenoP. J. (2015). Synthesis of arctic research (SOAR) in marine ecosystems of the Pacific Arctic.Prog. Oceanogr.136111. 10.1016/j.pocean.2015.05.017

  • 61

    MoranS. B.EllisK. M.SmithJ. N. (1997). 234Th/238U disequilibrium in the central Arctic ocean: implications for particulate organic carbon export.Deep Sea Res. II4415931606. 10.1016/S0967-0645(97)00049-0

  • 62

    MoranS. B.KellyR. P.HagstromK.SmithJ. N.GrebmeierJ. M.CooperL. W.et al (2005). Seasonal changes in POC export flux in the Chukchi Sea and implications for water column-benthic coupling in Arctic shelves.Deep Sea Res. II5234273451. 10.1016/j.dsr2.2005.09.011

  • 63

    MoranS. B.KellyR. P.IkenK.MathisJ. T.LomasM. W.GradingerR. (2012). Seasonal succession of net primary productivity, particulate organic carbon export, and autotrophic community composition in the eastern Bering Sea.Deep Sea Res. II65-708497. 10.1016/j.dsr2.2012.02.011

  • 64

    NaiduA. S.CooperL. W.GrebmeierJ. M.WhitledgeT. E.HameediM. J. (2004). “The continental margin of the north Bering-Chukchi Sea: concentrations, sources, fluxes, accumulation and burial rates of organic carbon,” in The Organic Carbon Cycle in the Arctic Ocean, edsSteinR.MacdonaldR. W. (Heidelberg: Springer-Verlag), 193203. 10.1007/978-3-642-18912-8

  • 65

    NeeleyA. R.HarrisL. A.FreyK. E. (2018). Unraveling phytoplankton community dynamics in the Northern Chukchi Sea under sea-ice-covered and sea-ice-free conditions.Geophys. Res. Lett.4576637671. 10.1029/2018GL077684

  • 66

    NodderS. D.CharetteM. A.WaiteA. M.TrullT. W.BoydP. W.ZeldisJ.et al (2001). Particulate transformations and export flux during an in situ iron-stimulated bloom in the Southern Ocean.Geophys. Res. Lett.2824092412. 10.1029/2001GL013008

  • 67

    OlliK.RiserC. W.WassmannP.RatkovaT.ArashkevichE.PasternakA. (2002). Seasonal variation in vertical flux of biogenic matter in the marginal ice zone and the central Barents Sea.J. Mar. Syst.38189204. 10.1016/S0924-7963(02)00177-X

  • 68

    PickartR. S.MooreG. W. K.ChongyuanM.BahrF.NobreC.WeingartnerT. J. (2016). Circulation of winter water on the Chukchi shelf in early summer.Deep Sea Res. II1305675. 10.1016/j.dsr2.2016.05.001

  • 69

    PiepenburgD. (2005). Recent research on Arctic benthos: common notions need to be revised.Polar Biol.28733755. 10.1007/s00300-005-0013-5

  • 70

    PlougH.GrossartH. P. (2000). Bacterial growth and grazing on diatom aggregates: respiratory carbon turnover as a function of aggregate size and sinking velocity.Limnol. Oceanogr.4514671475. 10.4319/lo.2000.45.7.1467

  • 71

    PlougH.JorgensenB. (1999). A net-jet flow system for mass transfer and microsensor studies of sinking aggregates.Mar. Ecol. Prog. Ser.176279290. 10.3354/meps176279

  • 72

    PostD. M. (2002). Using stable isotopes to estimate trophic position: models, methods and assumptions.Ecology83703718. 10.1890/0012-9658(2002)083[0703:usitet]2.0.co;2

  • 73

    Richter-MengeJ.DruckenmillerM. L.JeffriesM.(eds). (2019). Arctic Report Card 2019.https://www.arctic.noaa.gov/Report-Card

  • 74

    SalaM. M.ArrietaJ. M.BorasJ. A.DuarteC. M.VaquéD. (2010). The impact of ice melting on bacterioplankton in the Arctic Ocean.Polar Biol.3316831694. 10.1007/s00300-010-0808-x

  • 75

    SlatsR.OliverC.BahnkeR.BellH.MillerA.PungowiyiD.et al (2019). Voices from the Front Lines of a Changing Bering Sea. Arctic Report Card 2019. Available online at: http://www.arctic.noaa.gov/Report-Card(accessed January 2, 2020).

  • 76

    SpringerA. M.McRoyC. P. (1993). The paradox of pelagic food webs in the northern Bering Sea-III. Patterns of primary production.Cont. Shelf Res.13575599. 10.1016/0278-4343(93)90095-F

  • 77

    StabenoP. J.BellS. W. (2019). Extreme conditions in the bering Sea (2017-2018): record-breaking low sea-ice extent.Geophys. Res. Lett.4689528959. 10.1029/2019gl083816

  • 78

    ThomanR. L.BhattU. S.BieniekP. A.BrettschneiderB. R.BrubakerM.DanielsonS. L.et al (2020). The record low Bering Sea ice extent in 2018: context, impacts, and an assessment of the role of anthropogenic climate change.Bull. Am. Meteorol. Soc.101S53S58. 10.1175/BAMS-D-19-0175.1

  • 79

    VaquéD.LaraE.ArrietaJ. M.HoldingJ.E. L.HendriksI. E.et al (2019). Warming and CO2 enhance arctic heterotrophic microbial activity.Front. Microbiol.10:494. 10.3389/fmicb.2019.00494

  • 80

    VaughanD. G.ComisoJ. C.AllisonI.CarrascoJ.KaserG.KwokR.et al (2013). “Observations: Cryosphere,” in Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment of the Intergovernmental Panel on Climate Change, edsStockerT. F.QinD.PlattnerG.-K.TignorM.AllenS. K.BoschungJ.et al (Cambridge: Cambridge University Press), 317382. 10.1017/CBO9781107415324.012

  • 81

    WalshJ. J.McRoyC. P.CoachmanL. K.GoeringJ. J.NihoulJ. J.WhitledgeT. E.et al (1989). Carbon and nitrogen cycling within the Bering/Chukchi Seas: source regions for organic matter effecting AOU demands of the Arctic Ocean.Prog. Oceanogr.22277359. 10.1016/0079-6611(89)90006-2

  • 82

    WassmannP.ReigstadM. (2011). Future Arctic Ocean seasonal ice zones and implications for pelagic-benthic coupling.Oceanography24220231. 10.5670/oceanog.2011.74

  • 83

    WoodgateR. A.StaffordK. M.PrahlF. G. (2015). A synthesis of year-round interdisciplinary mooring measurements in the Bering Strait (1990-2014) and the RUSALCA years (2004-2011).Oceanography284667. 10.5670/oceanog.2015.57

  • 84

    WoollerM. J.ZazulaG. D.EdwardsM.FroeseD. G.BooneR. D.ParkerC.et al (2007). Stable carbon isotope compositions of eastern beringian grasses and sedges: investigating their potential as paleoenvironmental indicators.Arctic Antarct. Alpine Res.39318331. 10.1657/1523-0430(2007)39[318:scicoe]2.0.co;2

  • 85

    YuW.ChenL.ChengJ.HeJ.YinM.ZengZ. (2010). 234Th-derived particulate organic carbon export flux in the western Arctic Ocean.Chin J. Oceanog. Limnol.2811461151. 10.1007/s00343-010-9933-1

  • 86

    YuW.HeJ.LiY.LinW.ChenL. (2012). Particulate organic carbon export fluxes and validation of steady state model of 234Th export in the Chukchi Sea.Deep Sea Res. II81–846371. 10.1016/j.dsr2.2012.03.003

Summary

Keywords

carbon cycling, particulate organic carbon, Bering and Chukchi Sea Shelves, marine particles, marine snow, Arctic, climate change, biological carbon pump

Citation

O’Daly SH, Danielson SL, Hardy SM, Hopcroft RR, Lalande C, Stockwell DA and McDonnell AMP (2020) Extraordinary Carbon Fluxes on the Shallow Pacific Arctic Shelf During a Remarkably Warm and Low Sea Ice Period. Front. Mar. Sci. 7:548931. doi: 10.3389/fmars.2020.548931

Received

04 April 2020

Accepted

27 October 2020

Published

19 November 2020

Volume

7 - 2020

Edited by

Robyn E. Tuerena, Scottish Association For Marine Science, United Kingdom

Reviewed by

Henry Ruhl, Monterey Bay Aquarium Research Institute (MBARI), United States; Mark Andrew Stevenson, Newcastle University, United Kingdom

Updates

Copyright

*Correspondence: Stephanie H. O’Daly, ;

This article was submitted to Global Change and the Future Ocean, a section of the journal Frontiers in Marine Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics