The EXport Processes in the Ocean from RemoTe Sensing (EXPORTS) program of National Aeronautics and Space Administration focuses on linking remotely sensed properties from satellites to the mechanisms that control the transfer of carbon from surface waters to depth. Here, the naturally occurring radionuclide thorium-234 was used as a tracer of sinking particle flux. More than 950 234Th measurements were made during August–September 2018 at Ocean Station Papa in the northeast Pacific Ocean. High-resolution vertical sampling enabled observations of the spatial and temporal evolution of particle flux in Lagrangian fashion. Thorium-234 profiles were remarkably consistent, with steady-state (SS) 234Th fluxes reaching 1,450 ± 300 dpm m−2 d−1 at 100 m. Nonetheless, 234Th increased by 6%–10% in the upper 60 m during the cruise, leading to consideration of a non-steady-state (NSS) model and/or horizontal transport, with NSS having the largest impact by decreasing SS 234Th fluxes by 30%. Below 100 m, NSS and SS models overlapped. Particulate organic carbon (POC)/234Th ratios decreased with depth in small (1–5 μm) and mid-sized (5–51 μm) particles, while large particle (>51 μm) ratios remained relatively constant, likely influenced by swimmer contamination. Using an average SS and NSS 234Th flux and the POC/234Th ratio of mid-sized particles, we determined a best estimate of POC flux. Maximum POC flux was 5.5 ± 1.7 mmol C m−2 d−1 at 50 m, decreasing by 70% at the base of the primary production zone (117 m). These results support earlier studies that this site is characterized by a modest biological carbon pump, with an export efficiency of 13% ± 5% (POC flux/net primary production at 120 m) and 39% flux attenuation in the subsequent 100 m (POC flux 220 m/POC flux 120m). This work sets the foundation for understanding controls on the biological carbon pump during this EXPORTS campaign.
1. Introduction
Marine ecosystems are of fundamental importance with regard to their capacity to influence the storage, transformation, and fate of carbon (C) and associated elements in the Earth’s biosphere (Falkowski et al., 2000; Sarmiento and Gruber, 2002; Chavez et al., 2011; Doney et al., 2012). Yet the mechanisms that influence the movement of fixed C from well-lit surface waters to the ocean’s interior remain elusive (Siegel et al., 2016). Even small changes in the magnitude of downward transport of C in the ocean, or export, in global Earth system models can have major impacts on ocean C sequestration and thus on predicted atmospheric CO2 concentrations (Kwon et al., 2009). Subtle differences in export further influence biological food webs with regard to biodiversity and fisheries abundance (e.g., Beaugrand et al., 2010) as well as the fate of other biologically associated and particle-scavenged elements, such as heavy metals and organic contaminants (Shine et al., 1995; Gustafsson et al., 1997).
The primary goal of the EXport Processes in the Ocean from RemoTe Sensing (EXPORTS) program of the National Aeronautics and Space Administration is to develop a predictive understanding of the export and fate of fixed C from global ocean primary production (Siegel et al., 2016). Inherent in these models is the need to reduce uncertainties in C export estimates and to develop constructs that facilitate future C export estimates across a variety of spatial and temporal scales. As such, critical components of the EXPORTS Program are quantitative measurements of sinking particle fluxes and their attenuation versus depth at scales similar to the physical and biological processes that influence fixed C export.
Measuring particle export is confounded by high temporal and spatial variability on vertical scales of <10–100 m, horizontal scales of <3–5 km, and time scales of days to weeks (Benitez-Nelson and McGillicuddy, 2008; Estapa et al., 2015; Stukel et al., 2017). A mechanistic understanding of export processes has been further hampered by the few methods that are able to capture the non-steady-state nature of ocean C flux and its inherent patchiness. These methods are either direct measurement of C fluxes using a variety of sediment trap types (McDonnell et al., 2015), indirect methods such as using radioactive tracers (Buesseler et al., 1998; Verdeny et al., 2009; Haskell et al., 2015) or geochemical mass balance (e.g., Hamme et al., 2019). All these approaches sample over dissimilar spatial and temporal scales and with different operational biases (Burd et al., 2010; Estapa et al., 2015).
Here, we use the radioactive tracer, thorium-234 (half-life, t½ = 24.1 days) to quantify the spatio-temporal variability in particle flux from the well-lit surface layer (euphotic zone, Ez) and its attenuation with depth in the twilight zone below. 234Th is produced by the radioactive decay of uranium-238 (t½ = 4.47 × 109 years) in the water column. As 234Th is highly particle reactive, it is rapidly adsorbed onto sinking particles and released back into the dissolved phase when those particles are remineralized. The disequilibrium from its soluble conservative parent, 238U, provides quantitative information on where particle export (234Th deficiency due to net removal) and remineralization (234Th excess due to solubilization or fragmentation, i.e., conversion of sinking into nonsinking phases) occur in the water column (e.g., Buesseler et al., 1992). As such, 234Th provides rates of downward C fluxes over the smaller vertical and spatial scales required to understand processes that regulate water column export and C attenuation in a dynamic marine system.
This work was part of the EXPORTS field campaign that included satellite remote sensing, autonomous platforms, and observations from survey and process ships to gather key measurements of rates, stocks, and physical and biological processes, including primary production and particle flux (Siegel et al., 2020). We present an extensive set of three-dimensional time-series profiles of total 234Th (n = 61) that were collected in the late summer in the northeast subarctic Pacific near Ocean Station Papa (Station P). This high spatial and temporal resolution allowed an assessment of steady-state and non-steady-state particulate organic C (POC) flux throughout the upper water column using POC/234Th ratios from filtered particles collected using in situ pumps. These are among the earliest EXPORTS results and thus provide a foundation on which to build subsequent assessments of additional C cycle pathways and mechanisms and support EXPORTS goals of obtaining a predictive understanding of the export, fate, and impacts of global ocean net primary production (NPP).
2. Study Site
This first EXPORTS study was conducted near the time-series Station P in the northeast subarctic Pacific (50°N, 145°W) from August 9 to September 13, 2018. Station P was established in 1949 and maintains one of the oldest oceanic time-series programs. This high-nutrient low-chlorophyll region is characterized by a modest spring bloom that typically occurs in March and April and is comprised of a mixture of coccolithophores and diatoms (Boyd and Harrison, 1999; Fassbender et al., 2016). Prior to EXPORTS, particle fluxes were measured for more than two decades using deep-ocean sediment traps (Timothy et al., 2013). In addition, drifting trap arrays were deployed from 50 to 1,000 m during the period of 1987–1998 (Wong et al., 2002). Also important are prior studies at Station P that included water column 234Th measurements for determining sinking C fluxes (Charette et al., 1999; Kawakami et al., 2010; Mackinson et al., 2015).
The water column samples discussed here were collected onboard the R/V Sally Ride, which conducted a series of small-scale (approximately 900 km2) and large-scale (approximately 7,700 km2) surveys, from 49.5 to 50.9°N and 144.5 to 145.8°W, over a period of 28 days (Figure 1). These surveys were centered on a drifting Lagrangian float deployed at approximately 100 m, just below the base of the winter mixed layer and near the 0.1% penetration depth of photosynthetic active radiation (PAR; Siegel et al., 2020). Surveys were conducted over 3 time periods or “epochs,” each lasting 7–8 days. Epochs were set largely by the logistical constraints of sampling and the assumption that surface-derived particles would take 5–10 days to reach 500 m at sinking rates of 50–100 m d−1, which allows surface events to be tied to particle fluxes at depth within each epoch. During an epoch, total 234Th samples were collected from approximately 20 water column profiles using Niskin bottles attached to a rosette with a Sea-Bird 911 system configured with 2 SBE temperature and salinity sensors, 1 SBE 43 oxygen sensor, and 1 Digiquartz pressure sensor. Connected to the Sea-Bird 911 system were a WET Labs C-Star 650-nm beam transmissometer, a Biospherical QSE scalar PAR sensor, and a WET Labs FLBBRTD combination chlorophyll fluorometer and backscatter (700 nm) sensor, which we used to define the Ez. Each epoch also included the collection of size-fractioned particles from large volume in situ pumps (n = 12 casts in total).
3. Methods
3.1. Total thorium-234
Total 234Th sampling and analyses were similar to those described in the study of Buesseler et al. (2001) using a smaller volume of 2-L versus 4-L samples. This method is based on the formation of a manganese dioxide precipitate, which preferentially scavenges 234Th, leaving 238U in the dissolved phase. Unfiltered seawater was collected from 10-L Niskin bottles at 61 stations, each with 13–18 sample depths over the upper 500 m. Samples were immediately acidified to pH 1 with concentrated HNO3 acid and a 1-mL aliquot of 230Th (approximately 50 disintegrations per minute, dpm g−1) added for use as a yield monitor (Pike et al., 2005). Samples were allowed to equilibrate for 6–8 h and then adjusted to a pH of approximately 8.5 using concentrated NH4OH. Reagents KMnO4 and MnCl2 were added to precipitate the Mn oxide. Samples were allowed to stand for approximately 8 h, filtered onto 25-mm diameter quartz microfiber filters (QMA), and dried at 60 °C.
Once dry, the filters were mounted under Mylar and two layers of aluminum foil and counted on 1 of the 6 anticoincidence, 5-sample, low-level Risø β counters. Each β counter was surrounded by 8 cm of lead, resulting in extremely low background counts (average = 0.26 counts per minute, cpm). Samples were typically counted twice at sea with counting errors of <3%, and the average value was used after appropriate decay and ingrowth corrections. Close examination of decay curves and counting errors allowed us to identify outliers in cases where first and second counts disagreed by more than 0.2 dpm L−1. Samples were recounted >6 months after collection for background activity (approximately 0.48 ± 0.24 cpm) stemming from other naturally occurring β-emitting radionuclides that are included in the precipitate (most likely radium daughters; Benitez-Nelson et al., 2001b). Cross calibration and counting efficiencies of β detectors were assessed 2–3 times during and after the cruises using 238U standards that have an activity similar to the measured samples. Standards were rotated and counted on each of the 5 positions for the 6 different Risø detectors. The detector efficiencies varied by <2% in 24 of the 30 detectors and by <5% in the remaining detectors. 234Th values are reported in activity units (dpm L−1) and are provided in the SeaBASS Data Repository (https://seabass.gsfc.nasa.gov/archive/WHOI/buesseler/EXPORTS/EXPORTSNP/archive/).
All total 234Th activities were corrected for the overall efficiency of the Mn oxide precipitation method as determined using 230Th as a yield monitor. Recovery of 230Th was conducted after the final background β count using procedures detailed in the study of Pike et al. (2005), modified to no longer use an ion-exchange column to reduce Mn interferences. Precipitates were dissolved in a nitric acid and hydrogen peroxide solution, and a known amount of 229Th was added to each sample. The 230Th/229Th ratio was then analyzed by inductively coupled plasma mass spectrometry, thus allowing for the determination of 230Th present in each sample and correction for Th loss during processing. Recoveries averaged 88 ± 9%. Samples with recoveries <40% were excluded from the data set (n = 2). For samples with apparent recoveries >100% (n = 18, average recovery = 107.3 ± 8.5%), recoveries were fixed at 100%, and sample errors were increased to 5% to calculate final results. Variations in the amount of 230Th added on different casts required cast-by-cast consideration of final recoveries. The 10 samples collected at depth (2 casts of n = 5 at 3,000 m) confirmed secular equilibrium with 238U (2.41 dpm L−1), with average 234Th activities of 2.40 ± 0.14 and 2.42 ± 0.15 dpm L−1, measured early and late in the cruise, respectively. Errors on individual 234Th activities were propagated from the initial and final counts as well as from the error on the 230Th-derived yield.
3.2. Size-fractionated particles collected using in situ pumps
Size-fractionated particles were collected by Large Volume Water Transfer System battery-powered in situ pumps (McLane Industries) that were modified to accommodate two 142-mm diameter mini-MULVFS (Multiple Unit Large Volume Filtration System) filter holders (Lam et al., 2015). Size-fractionated particles for 234Th and other isotope, element, and compound analyses were collected using a single flow path through a 3-stage filter holder. Pumps were deployed on a total of 12 casts at 50 m and at the depths of the sediment traps: 100, 150, 200, 330, and 500 m ± 5 m (Estapa et al., 2020). A pressure sensor attached to the deepest pump confirmed depths, which can vary depending on the wire angle. Filtration volumes ranged from 700 to 1,500 L for 4–5 hours at an initial flow rate of 8 L min−1. Filter holders were loaded with 51- and 5-µm pore size acid-cleaned Nitex screens followed by a 1-µm pore size precombusted QMA filter. Sampling at these 6 depths at 12 stations generated 72 samples for each of the 3 size classes (1–5, 5–51, and >51 μm). Once per epoch a filter holder containing a full set of filters, but without any active water pumping, was mounted on the deepest in situ pump at 500 m to obtain process (or dipped) blanks (Lam et al., 2015). These blanks were processed in the same manner as the samples. Size-fractionated particles were analyzed for 234Th, particulate carbon (PC) and particulate inorganic carbon (PIC), and other analytes including particulate nitrogen, biogenic silica, phosphorus, pigments, 210Po/210Pb, barium, Ba isotopes, and other organic compounds. Results for these other analytes will be presented elsewhere. During 3 pump deployments, additional samples (n = 5) for the smallest particle size class (1–5 μm) were collected at 20, 85, and 320 m by deploying extra in situ pumps. These samples were collected using a single flow path through a mini-MULVFS filter holder loaded with 51- and 6-µm pore size Nitex screens followed by paired QMA filters and paired 0.3-µm pore size GF75 filters.
Particles were gently rinsed off the Nitex screens onto a 25-mm diameter 1.2-µm pore size silver (Ag) filters using 1-µm filtered seawater collected at depths from 330 to 500 m that was stored refrigerated at 4 °C and collected once per epoch. At these depths, dissolved organic C concentrations were approximately 45 μM (Bif and Hansell, 2019), and appropriate C filter processing blanks were assessed. At 9 stations, the entire screens were rinsed, whereas at 3 stations, the screens were cut in half using a rotary blade to allow for analysis of other analytes according to the GEOTRACES sample-handling protocol (Cutter et al., 2017). QMA filters were subsampled with punches of 21–26 mm diameter for 234Th, PC, and PIC, as particle distribution has been shown to be uniform across the filter (Maiti et al., 2012). All filters were dried in a 60 °C oven, β counted at sea for 234Th activities, and recounted 6 months later onshore for backgrounds. After counting, screen samples that had been rinsed onto Ag filters were split by weight (Lamborg et al., 2008), and one-third of each sample was analyzed for PC and PIC. PC was analyzed using high-temperature combustion on a Thermo Electron FlashEA 1112 C/N analyzer at the WHOI Nutrient Analytical Facility. PIC was determined by coulometric analysis of CO2 after acidification with phosphoric acid (Honjo et al., 1995).
The average of all process blanks was subtracted from total 234Th, PC, and PIC measurements. Process blanks for 234Th and PC were <2% of total 234Th and PC (Ag filters: 0.15 ± 0.11 dpm 234Th, 1/3 Ag filters: 0.17 ± 0.28 μmol C; QMA punches: 0.13 ± 0.06 dpm 234Th, 0.15 ± 0.21 μmol C). In contrast, PIC concentrations were quite low (<3% of total C). Therefore, PIC process blanks were relatively higher, accounting for 82% of the measured PIC on 1/3 Ag filters (0.47 ± 0.03 μmol C) and 42% of PIC on the QMA punches (0.23 ± 0.01 μmol C). Three replicate punches from a 50 and a 500-m sample were analyzed for PC and PIC, obtaining a standard deviation of <5% for both analyses. POC concentrations were obtained from the difference between blank-corrected PC and PIC results. Data for POC have associated uncertainties derived from these blank corrections and weighing errors of the analytical balance used for splitting the Ag filters. Particulate 234Th data have an associated uncertainty from counting and blank corrections. All size-fractionated particle data from the 2018 EXPORTS expedition are found in the SeaBASS Data Repository (https://seabass.gsfc.nasa.gov/archive/WHOI/buesseler/EXPORTS/EXPORTSNP/archive/).
Swimmers (nondetrital zooplankton) observed by the naked eye were picked from the 51-µm screens. However, complete removal of swimmers was not always possible due to the high abundance of tiny copepods (<0.5 mm) that were difficult to remove without eliminating a significant amount of the passive sinking particles. Subsequent photographic analysis confirmed the presence of larger zooplankton on some of the Ag filters, particularly at depths >100 m. Seven (of 72) samples were identified as outliers based on high POC/234Th (hereafter referred to as C/Th) values following the interquartile method by Tukey (1977). Outliers were defined by measurements with values greater than the third quartile by more than 1.5 times the interquartile range (i.e., third quartile–first quartile) and were removed from subsequent >51 µm data analysis. The remaining large particle data were characterized by a normal C/Th distribution, as observed in the other two pump-derived size classes.
4. Results
4.1. Ancillary data
As part of this study, we derived chlorophyll (chl) concentrations from fluorescence on the same casts as those where total 234Th was sampled. These data are shown in Figure 1, with the x-axis shifted to the right for clarity. Chlorophyll was quantified from the raw conductivity, temperature, and depth (CTD) fluorescence sensor data (CTD ChlF, measured in volts, V) using chl (µg L−1) = [raw CTD ChlF (V) – dark value (V)]/[slope of calibration curve (V/(µg L−1)] with a calibration slope = 0.64 ± 0.01 V/(µg L−1) and dark value = 0.0870 ± 0.0009 V, as derived by the EXPORTS calibration team.
Previous work has shown that chl can be used to define the base of the sunlit Ez. Using CTD fluorescence sensor data, Owens et al. (2015) defined this base of the “primary production zone” (PPZ) as the depth at which the signal from in situ fluorescence dropped to 10% of the maximum above that depth (after correction for sensor background signal). On our CTD casts, the PPZ (117 ± 5 m) overlapped with the zone of net 234Th removal (234Th/238U < 1.0; Figure 1) similar to many previous studies (e.g., Buesseler et al., 2020). PAR was also measured as an indicator of the Ez depth and is shown in Figure 1. The depth where 234Th reaches equilibrium fell between the depth of 1% (Ez1.0 = 79 ± 6 m) and 0.1% PAR (Ez0.1 = 118 ± 9 m).
4.2. Thorium-234 activity profiles
The results of all 760 analyses of total 234Th are provided in Figure 1, with the average (± standard deviation) from all three epochs shown versus depth and density. This depiction allows for more careful alignment of individual vertical profiles, given depth variability in physical properties. Average 234Th and 238U activities are also provided in Table 1, with individual data available in SeaBASS and in Table S1 (along with the calculated 234Th and POC fluxes discussed in the following sections). Salinity varied minimally throughout the sampling periods, mainly increasing with greater depth, resulting in 238U activities ranging from 2.22 dpm L−1 at the surface to 2.36 dpm L−1 by 500 m. Overall, 234Th activities averaged 1.37 ± 0.15 dpm L−1 in the upper 10 m and reached equilibrium with 238U around a density of 25.75 kg m−3 (about 80–90 m). Below this depth, average 234Th activities exceeded that of 238U, before returning to equilibrium at a density of 26.5 kg m−3 (about 180 m).
Avg. Depth (m)a . | Range (m)b . | 234Th Activity (dpm L−1) . | 238U Activity (dpm L−1) . | 234Th Flux (dpm m−2 d−1) . | |||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Epoch 1 . | Epoch 2 . | Epoch 3 . | Cruise . | Cruise . | Cruise SSd . | Cruise NSSe . | |||||||||||
Average . | SDc . | n . | Average . | SDc . | n . | Average . | SDc . | n . | Average . | SDc . | Average . | Average . | Uncert . | Average . | Uncert . | ||
5 | 1.35 | 0.17 | 18 | 1.31 | 0.16 | 21 | 1.46 | 0.09 | 19 | 1.37 | 0.15 | 2.22 | 324 | 80 | 221 | 120 | |
20 | 15–24 | 1.37 | 0.12 | 17 | 1.34 | 0.14 | 21 | 1.50 | 0.13 | 16 | 1.40 | 0.14 | 2.22 | 666 | 145 | 445 | 195 |
35 | 25–42 | 1.43 | 0.14 | 14 | 1.41 | 0.16 | 22 | 1.57 | 0.13 | 23 | 1.48 | 0.16 | 2.23 | 964 | 181 | 637 | 254 |
50 | 45–56 | 1.63 | 0.15 | 18 | 1.61 | 0.20 | 21 | 1.73 | 0.13 | 19 | 1.66 | 0.17 | 2.24 | 1,276 | 201 | 789 | 307 |
65 | 59–69 | 1.85 | 0.23 | 12 | 1.91 | 0.22 | 14 | 1.90 | 0.17 | 17 | 1.89 | 0.20 | 2.24 | 1,407 | 254 | 908 | 386 |
80 | 75–83 | 2.13 | 0.21 | 17 | 2.12 | 0.19 | 21 | 2.11 | 0.17 | 21 | 2.12 | 0.19 | 2.25 | 1,444 | 279 | 997 | 469 |
100 | 94–109 | 2.26 | 0.24 | 20 | 2.21 | 0.15 | 23 | 2.24 | 0.13 | 20 | 2.24 | 0.18 | 2.26 | 1,452 | 304 | 1,024 | 545 |
120 | 110–130 | 2.38 | 0.21 | 34 | 2.33 | 0.19 | 40 | 2.34 | 0.18 | 41 | 2.35 | 0.19 | 2.29 | 1,428 | 365 | 1,026 | 672 |
150 | 145–151 | 2.37 | 0.18 | 18 | 2.32 | 0.10 | 22 | 2.38 | 0.15 | 20 | 2.36 | 0.15 | 2.33 | 1,384 | 422 | 941 | 861 |
200 | 194–222 | 2.34 | 0.17 | 22 | 2.38 | 0.19 | 23 | 2.37 | 0.17 | 22 | 2.36 | 0.18 | 2.34 | 1,283 | 585 | 774 | 1,263 |
280 | 249–314 | 2.38 | 0.09 | 14 | 2.36 | 0.29 | 6 | 2.30 | 0.11 | 6 | 2.36 | 0.16 | 2.34 | 1,192 | 447 | 1,100 | 1,382 |
330 | 328–334 | 2.34 | 0.15 | 17 | 2.28 | 0.14 | 21 | 2.32 | 0.15 | 20 | 2.31 | 0.14 | 2.35 | 1,464 | 885 | 1,270 | 1,961 |
500 | 498–503 | 2.35 | 0.17 | 12 | 2.35 | 0.18 | 14 | 2.38 | 0.22 | 13 | 2.36 | 0.19 | 2.36 | 1,527 | 1,119 | 1,116 | 2,437 |
Avg. Depth (m)a . | Range (m)b . | 234Th Activity (dpm L−1) . | 238U Activity (dpm L−1) . | 234Th Flux (dpm m−2 d−1) . | |||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Epoch 1 . | Epoch 2 . | Epoch 3 . | Cruise . | Cruise . | Cruise SSd . | Cruise NSSe . | |||||||||||
Average . | SDc . | n . | Average . | SDc . | n . | Average . | SDc . | n . | Average . | SDc . | Average . | Average . | Uncert . | Average . | Uncert . | ||
5 | 1.35 | 0.17 | 18 | 1.31 | 0.16 | 21 | 1.46 | 0.09 | 19 | 1.37 | 0.15 | 2.22 | 324 | 80 | 221 | 120 | |
20 | 15–24 | 1.37 | 0.12 | 17 | 1.34 | 0.14 | 21 | 1.50 | 0.13 | 16 | 1.40 | 0.14 | 2.22 | 666 | 145 | 445 | 195 |
35 | 25–42 | 1.43 | 0.14 | 14 | 1.41 | 0.16 | 22 | 1.57 | 0.13 | 23 | 1.48 | 0.16 | 2.23 | 964 | 181 | 637 | 254 |
50 | 45–56 | 1.63 | 0.15 | 18 | 1.61 | 0.20 | 21 | 1.73 | 0.13 | 19 | 1.66 | 0.17 | 2.24 | 1,276 | 201 | 789 | 307 |
65 | 59–69 | 1.85 | 0.23 | 12 | 1.91 | 0.22 | 14 | 1.90 | 0.17 | 17 | 1.89 | 0.20 | 2.24 | 1,407 | 254 | 908 | 386 |
80 | 75–83 | 2.13 | 0.21 | 17 | 2.12 | 0.19 | 21 | 2.11 | 0.17 | 21 | 2.12 | 0.19 | 2.25 | 1,444 | 279 | 997 | 469 |
100 | 94–109 | 2.26 | 0.24 | 20 | 2.21 | 0.15 | 23 | 2.24 | 0.13 | 20 | 2.24 | 0.18 | 2.26 | 1,452 | 304 | 1,024 | 545 |
120 | 110–130 | 2.38 | 0.21 | 34 | 2.33 | 0.19 | 40 | 2.34 | 0.18 | 41 | 2.35 | 0.19 | 2.29 | 1,428 | 365 | 1,026 | 672 |
150 | 145–151 | 2.37 | 0.18 | 18 | 2.32 | 0.10 | 22 | 2.38 | 0.15 | 20 | 2.36 | 0.15 | 2.33 | 1,384 | 422 | 941 | 861 |
200 | 194–222 | 2.34 | 0.17 | 22 | 2.38 | 0.19 | 23 | 2.37 | 0.17 | 22 | 2.36 | 0.18 | 2.34 | 1,283 | 585 | 774 | 1,263 |
280 | 249–314 | 2.38 | 0.09 | 14 | 2.36 | 0.29 | 6 | 2.30 | 0.11 | 6 | 2.36 | 0.16 | 2.34 | 1,192 | 447 | 1,100 | 1,382 |
330 | 328–334 | 2.34 | 0.15 | 17 | 2.28 | 0.14 | 21 | 2.32 | 0.15 | 20 | 2.31 | 0.14 | 2.35 | 1,464 | 885 | 1,270 | 1,961 |
500 | 498–503 | 2.35 | 0.17 | 12 | 2.35 | 0.18 | 14 | 2.38 | 0.22 | 13 | 2.36 | 0.19 | 2.36 | 1,527 | 1,119 | 1,116 | 2,437 |
Abbreviation: Uncert = uncertainity.
aAverage sample depth in the bin.
bDepth range of the samples included in the bin.
cStandard deviation.
dCruise-average steady-state (SS) 234Th fluxes were calculated from all the SS 234Th fluxes computed for each station. The uncertainty corresponds to the standard deviation across all of the SS 234Th flux estimates.
eNon-steady-state (NSS) 234Th fluxes were calculated from the depth-averaged 234Th activities of Epochs 1 and 3 following Equation 4. The uncertainty was calculated by propagating the standard deviations of 234Th activity from Epochs 1 and 3 through Equation 4 (see Section 4.5).
The spatio-temporal variability of 234Th activity was fairly consistent per epoch, with a cruise-average standard deviation of 0.17 dpm L−1 per depth (Figure 1), roughly 6%–11% of the mean 234Th activity throughout the water column (Table 1). Mean 234Th activities in Epochs 1 and 2 were similar. In contrast, they significantly increased from Epoch 1 to Epoch 3 by 6%–10% at depths shallower than 60 m or densities less than 25.0 kg m−3 (Figure S1 and Table 1). However, the overall shape of the vertical profiles did not change (Figure 1).
4.3. Thorium-234 flux modeling
The flux of 234Th on sinking particles may be determined using a general thorium activity balance calculated using vertically integrated 234Th activities above a given depth horizon:
where is the change in the vertically integrated 234Th activity (dpm m−2) over a given depth interval (z) over time, t, λ is the radioactive decay constant of 234Th determined from the half-life of 234Th, 0.02876 d−1, V is the sum of physical advection and diffusion processes, and 234ThE is the thorium flux associated with particle export (dpm m−2 d−1) from the integrated depth zone of interest. 238U activities (dpm L−1) at each depth are determined using salinity (S) as measured using the CTD and Equation 2 (Owens et al., 2011):
All fluxes of 234Th were driven by the difference in vertically integrated 234Th versus 238U activities. As such, this difference became smaller with increasing depth, as 234Th approached equilibrium with 238U and 234Th-derived flux errors increased with depth for a given profile (Table 1). Sampling depth intervals provided additional error as depth sampling increased to ≥50 m below 150 m. This wider depth resolution limited definition of the integrated 234Th/238U ratio (Table 1).
4.3.1. Steady-state model
Given the small spatial variability, a 1-dimensional steady-state (SS) model was used to determine 234Th-derived export fluxes (Savoye et al., 2006). Integrating the activity balance from the surface to depth (D) allowed a calculation of the net 234Th export over a specific depth zone in dpm m−2 d−1 (Equation 3).
We computed a SS export flux for each individual 234Th profile (Table S1) and then computed a cruise average for each depth bin (Table 1). The uncertainty for the cruise average was determined using the standard deviation across all 234Th flux SS estimates. No errors were assigned to 238U profiles, as variability in 238U versus salinity has been observed to be low in the vast majority of open ocean data, with the exception of a few locations where rivers, ice, or coastal and low oxygen conditions result in local, nonconservative behavior of 238U (e.g., Owens et al., 2011; Xie et al., 2020).
Fluxes of 234Th increased with a similar slope throughout the upper 50 m at all stations, from an average of 324 ± 80 dpm m−2 d−1 in the upper 10 m to 1,450 ± 300 dpm m−2 d−1 by 80–100 m for the SS model (Figure 2). The highest 234Th-derived SS fluxes within a single station approach 2,500 dpm m−2 d−1 at 125 m, with a few reaching 3,000 dpm m−2 d−1 or higher by 200 m (Figure 2 and Table S1). In contrast, some of the lowest fluxes overlapped with zero at 200 m. At greater depths, there is more variability as SS fluxes continued to increase at some stations, while other SS fluxes decreased more sharply due to a greater 234Th excess below 80–100 m.
While the assumption of SS essentially suggests that no change in 234Th flux occurred over the 28-day sampling period, the SS model can also be used to infer fluxes prior to sampling, depending on the response time of the radionuclide. For 234Th, the SS response time in a given water parcel is shorter than its mean life with respect to decay (1/234λ = 35 days) due to scavenging by sinking particles (1/k, where k is first-order removal rate; Turnewitsch et al., 2008). Here, k was calculated from the removal flux, that is, the rate of scavenging onto sinking particles divided by the total 234Th inventory. Using the upper 50 m 234Th inventory as an example (Table 1), k was calculated to be on the order of approximately 0.02 d−1, and thus the total response time for 234Th was approximately 20 days (residence time of total 234Th = 1/[234λ + k]). We therefore consider the calculated SS fluxes to also represent the average flux over the 3 weeks prior to the start of sampling.
4.3.2. Physical processes versus non-steady-state model
Water column 234Th activities within the upper 60 m increased a small but significant amount between Epochs 1 and 3 (Figure 3 and Figure S1), either due to a temporal weakening in particulate export and/or to the transport of waters with higher 234Th activities into the sampling area from below (upwelling/eddy turbulent mixing) or from horizontal advection and diffusion. We examined each component separately, as the data did not allow for both temporal (non-steady-state, NSS) and physical processes to be resolved using Equation 1.
4.3.3. Non-steady-state model
Non-steady-state model As sampling occurred following a Lagrangian framework, the increase in 234Th activity during the field campaign was most likely due to temporal variability associated with changes in biological processes (see Section 5.1). If we assume that the increase in 234Th activity during the cruise was only due to NSS processes, we can solve Equation 1 and assume that equals the change in 234Th activity observed between epochs:
where the average 234Th activity of Epochs 1 and 3 are 234The1 and 234The3, respectively, and Δt is the time between the midpoints of Epochs 1 and 3 (18 days).
The choice of using Epochs 1 and 3 (and excluding Epoch 2) was motivated by two limitations of the NSS model. First, the NSS model is valid only if the same water mass is sampled over time (Savoye et al., 2006), and if the measurements capture temporal changes in 234Th activity in a Lagrangian framework versus a spatial gradient (Resplandy et al., 2012). Second, the uncertainty on the NSS model is relatively high due to error propagation, such that on timescales < 10 days, the errors are too large to be meaningful (Savoye et al., 2006). Epochs 1 and 3 were sampled 18 days apart and covered a small spatial area (approximately 35 km × 35 km) following a central Lagrangian float. These sampling periods therefore provided the best conditions for application of the NSS model (Resplandy et al., 2012). In contrast, Epoch 2 was a survey designed to capture larger spatial gradients (approximately 88 km × 88 km) that took place within 10 days of both Epochs 1 and 3.
To compute the NSS model, 234Th activity measurements were averaged by depth bin for Epochs 1 and 3. Each depth bin included 12–41 234Th measurements per epoch, with the exception of the 280 m depth bin, which was on average more sparsely sampled (Table 1). NSS integrated export flux at each depth D was then calculated according to Equation 4. The uncertainty on the NSS flux estimate includes the standard deviation of 234Th activity measurements in each depth bin per epoch propagated through Equation 4. The NSS models yield export estimates that are approximately 30% lower than the SS model in the upper 60 m and persist throughout the water column as this lower flux is propagated to deeper depths (Table 1 and Figure 3b).
4.3.4. Consideration of physical transport
Consideration of physical transport If we assume that the increase in 234Th activity during the cruise was only due to physical processes, we can evaluate the influence of the physical transport term, V, as follows:
where u, v, w, and are the zonal, meridional, and vertical components of the velocity and the thorium gradient, and k is the vertical eddy diffusivity coefficient. The first three terms represent horizontal and vertical advection while the last term represents the effect of vertical eddy turbulent mixing.
The vertical transport was estimated using measured vertical 234Th gradients, vertical velocity, and eddy diffusivity estimated at Station P from 2 years of Seaglider data by Pelland et al. (2017). Vertical velocity, w, is largely associated with seasonal Ekman pumping and was equal to 0.25 m d−1 at the time of the cruise. Vertical diffusivity, k, decreased with depth from 10−3 m2 s−1 above the transition layer (<35 m), to 10−4 m2 s−1 within the transition layer (35–80 m), to 10−5 m2 s−1 at depths >80 m. Together vertical advection and vertical eddy mixing supplied less than 20 dpm m−2 d−1 in the upper 25 m and influenced 234Th export by <6% in all epochs (Figure S2).
The horizontal transport was estimated at a scale of 50–100 km, where we had both 234Th activity gradients and horizontal velocities of 2–3 km d−1, mostly to the northeast, derived from acoustic Doppler current profiler measurements in the upper 25 m of the mixed layer (Siegel et al., 2020). Thorium-234 activities varied on scales <10 km during the survey (Figure S3), indicating the presence of small-scale spatio-temporal variations that were not considered in the transport estimates. Here, we used the average difference in 234Th between Epochs 1 and 3 (1.36 dpm L−1 versus 1.48 dpm L−1 in the upper 25 m, the shallowest 2 depth bins in Table 1) and a maximum horizontal velocity of 3 km d−1 (0.035 m s−1) across the distance of the sampling patch (75 km) to calculate a horizontal gradient of 1.6 × 10−6 dpm L−1 m−1 along the direction of the northeastward flow. This yielded a maximum horizontal input term of 0.005 dpm L−1 d−1 or 125 dpm m−2 d−1. This is less than 20% of the 234Th export derived from radioactive production and decay using the SS model (666 dpm m−2 d−1 in upper 25 m, Table 1), decreasing in importance with depth due to a reduction in current velocities.
Combined, the above results suggest that along the Lagrangian trajectory of the float, horizontal and vertical transport added at most 145 dpm m−2 d−1 of 234Th between Epochs 1 and 3 in the upper 25 m, with the magnitude relative to the total 234Th flux decreasing with increasing depth. In other words, transport terms can reduce 234Th fluxes in the surface by about 20% due to the input of higher 234Th activity waters from below or transported horizontally. The SS model with physical transport is qualitatively consistent with 234Th fluxes but would have a smaller impact on them than if calculated by the NSS model. As noted above, we cannot separate NSS from physical transport processes; therefore, we chose to consider 2 end-member models, that is, the SS model without transport (high end-member) and the NSS model (low end-member).
Importantly, the difference between the 2 end-member models (referred to as SS and NSS in the text) becomes indistinguishable within uncertainties at depth (Table 1). Given the error estimates in applying SS versus NSS models and the assumptions inherent in the modeling exercise, we used the average of both models to determine a best estimate of POC flux during the entire cruise, propagating the error of both models on the average 234Th-derived POC flux at a given depth horizon (see Section 4.5).
4.4. Particulate C/234Th ratios
To estimate the flux of POC from 234Th, one must use the C/Th ratio associated with sinking particles. Generally, C/Th ratios are more variable in shallow waters, decrease with depth, and are similar among large size class particles and those found in sediment traps (Buesseler et al., 2006). There are some exceptions, attributed most often to C-rich “swimmers” in the largest particles (Buesseler et al., 2007; Liu et al., 2009). Here we present data from the size-fractionated particles separated via filters and screens into >1–5 µm (small), 5–51 µm (mid), and > 51 µm (large) size classes (Figure 4). No significant temporal changes in C/Th were observed. As such, we focus the remainder of this section on changes in the C/Th ratios across all epochs with increasing depth and with differing size class. Results are summarized in Table 2.
Depth (m) . | POC/234Th (μmol dpm−1)a . | |||||
---|---|---|---|---|---|---|
1–5 μm . | 5–51 μm . | >51 μmb . | ||||
Average . | SD . | Average . | SD . | Average . | SD . | |
50 | 2.95 | 0.47 | 5.30 | 1.90 | 2.16 | 0.36 |
65 | 2.91 | 0.43 | 4.20 | 1.26 | 2.14 | 0.38 |
80 | 2.87 | 0.39 | 3.09 | 0.75 | 2.12 | 0.39 |
100 | 2.82 | 0.34 | 1.62 | 0.26 | 2.09 | 0.41 |
120 | 2.44 | 0.32 | 1.41 | 0.28 | 2.00 | 0.44 |
150 | 1.85 | 0.28 | 1.12 | 0.29 | 1.85 | 0.47 |
200 | 1.30 | 0.21 | 1.03 | 0.06 | 2.15 | 0.38 |
280 | 1.03 | 0.14 | 0.97 | 0.12 | 2.32 | 0.65 |
330 | 0.89 | 0.10 | 0.94 | 0.15 | 2.44 | 0.84 |
500 | 0.82 | 0.10 | 0.88 | 0.16 | 1.83 | 0.51 |
Depth (m) . | POC/234Th (μmol dpm−1)a . | |||||
---|---|---|---|---|---|---|
1–5 μm . | 5–51 μm . | >51 μmb . | ||||
Average . | SD . | Average . | SD . | Average . | SD . | |
50 | 2.95 | 0.47 | 5.30 | 1.90 | 2.16 | 0.36 |
65 | 2.91 | 0.43 | 4.20 | 1.26 | 2.14 | 0.38 |
80 | 2.87 | 0.39 | 3.09 | 0.75 | 2.12 | 0.39 |
100 | 2.82 | 0.34 | 1.62 | 0.26 | 2.09 | 0.41 |
120 | 2.44 | 0.32 | 1.41 | 0.28 | 2.00 | 0.44 |
150 | 1.85 | 0.28 | 1.12 | 0.29 | 1.85 | 0.47 |
200 | 1.30 | 0.21 | 1.03 | 0.06 | 2.15 | 0.38 |
280 | 1.03 | 0.14 | 0.97 | 0.12 | 2.32 | 0.65 |
330 | 0.89 | 0.10 | 0.94 | 0.15 | 2.44 | 0.84 |
500 | 0.82 | 0.10 | 0.88 | 0.16 | 1.83 | 0.51 |
POC = particulate organic carbon.
aPOC/234Th (C/Th) ratios for size-fractionated particles were measured at 50, 100, 150, 200, 330, and 500 m. Average C/Th at these depths were obtained from the arithmetic mean (± standard deviation, SD). Values shown in italics are linearly interpolated (see Section 4.4).
bSeven >51 μm samples (of 72) were identified as outliers due to swimmer contamination (see Section 3.2) and therefore are excluded from the average.
Average C/Th ratios in all size-fractionated particles were determined from the arithmetic mean (± standard deviation) at each sampled depth. C/Th ratios determined from the slopes of linear relationships between POC concentrations and 234Th activities produced comparable results (Table S2). Here, we focus on depth-specific arithmetic means to enable consistency across particle size classes. C/Th ratios at 65, 80, 120, and 280 m were interpolated linearly between sample depths (shown in italics, Table 2).
For small particles, C/Th ratios averaged 2.82–2.95 μmol dpm−1 at 50–100 m and decreased to 0.82–0.89 μmol dpm−1 by 330–500 m (Table 2). Two samples collected within the mixed layer at 20 m (Figure 4a) showed C/Th ratios a factor of 1.3 higher (3.95 ± 0.17 μmol dpm−1) than that measured below the mixed layer at 50–100 m.
For mid-sized particles, the average C/Th ratios at 50 m, 5.3 ± 1.9 μmol dpm−1, were almost double those measured in the small particles at the same depth, yet decreased sharply to 1.12–0.88 μmol dpm−1 by 150–500 m (Table 2 and Figure 5). The highest variability occurred in the upper 50 m (Figure 4b), consistent with enhanced primary production in the Ez.
Large particles showed similar C/Th ratios with depth (Figure 4c), ranging from 1.83 to 2.44 μmol dpm−1 excluding 7 outliers (see Section 3.2). We hypothesize that the lack of a decrease in C/Th ratio with increasing depth is likely due to the presence of zooplankton contamination. While nondetrital zooplankton identified by the naked eye were picked from the 51-µm screens, subsequent photographic analysis of the samples indicated the possible presence of nonpicked swimmers in about 20–30 samples, although only 7 of these samples had statistically higher C/Th ratios (Figure 4c; discussion at end of Section 3.2).
In most prior studies, C/Th ratios decrease with depth due to the preferential heterotrophic utilization of C as a food source (e.g., see reviews by Buesseler et al., 2006; Puigcorbé et al., 2020). The high and invariant C/Th ratio with depth in the >51 µm pump samples suggests pervasive contamination by zooplankton swimmers. During the EXPORTS campaign, individual samples were collected on the process ship using zooplankton nets (D Steinberg and her EXPORTS team). Both large and small pteropods were characterized by C/Th ratios of approximately 24, radiolarians approximately 44, and up to 825 μmol dpm−1 for small amphipods (data not shown). We estimate that the presence of just 1 amphipod or 2 radiolarians in a given sample could have increased the C/Th ratio by a factor of 2. This estimate supports our hypothesis of swimmer contributions resulting in higher C/Th in >51 µm filtered samples. Therefore, we consider the large size class as an unreliable indicator of the C/Th ratio of passively sinking particles.
4.5. POC flux profiles
With a known C/Th ratio versus depth, the POC flux can be calculated at any depth and location where 234Th fluxes are available. This conversion to C flux fails only if the C/Th ratio used differs from the flux-weighted C/Th ratio representative of sinking particles. Given the similarity across all epochs in C/Th ratios for a given size class, especially at depths below 50 m (Figure 4), depth-specific mean C/Th ratios were used to estimate POC flux. As such, station-by-station differences in POC flux were driven by variations in 234Th flux.
We determined a best estimate of POC flux using the average of the POC fluxes calculated from the SS and NSS models and using C/Th from the mid-sized particles (Table 3). Overall, this average cruise estimate showed a POC flux peak of 5.5 ± 1.7 mmol C m−2 d−1 at 50 m. This depth was above the 1% PAR and thus represented a layer where POC was likely still being produced by autotrophs (Figure 6). Note that there is no mid-sized C/Th data above 50 m, and hence, the 50 m ratio can be used to roughly estimate the POC flux at shallower depths where 234Th results were available. However, as a result, POC fluxes at depths shallower than 50 m might be underestimated, and they were certainly less well constrained. There was a sharp decline in the average POC fluxes below the 50 m POC flux maximum to approximately 1.7 mmol C m−2 d−1 at the 0.1% PAR and PPZ depths (120 m). Below these deeper Ez boundaries, there was only a modest decrease in POC fluxes to between 1.1 and 1.3 mmol C m−2 d−1 (Table 3).
. | POC Flux SS (mmol C m−2 d−1)a . | POC Flux NSS (mmol C m−2 d−1)b . | Average POC Flux (mmol C m−2 d−1)c . | |||
---|---|---|---|---|---|---|
. | 5–51 μm . | 5–51 μm . | 5–51 μm . | |||
Depth (m) . | Average . | Uncert . | Average . | Uncert . | Average . | Uncert . |
50 | 6.77 | 2.65 | 4.19 | 2.22 | 5.48 | 1.73 |
65 | 5.91 | 2.07 | 3.81 | 1.98 | 4.86 | 1.43 |
80 | 4.47 | 1.38 | 3.08 | 1.63 | 3.77 | 1.07 |
100 | 2.35 | 0.62 | 1.66 | 0.92 | 2.01 | 0.56 |
120 | 2.02 | 0.66 | 1.45 | 0.99 | 1.74 | 0.60 |
150 | 1.55 | 0.62 | 1.05 | 1.00 | 1.30 | 0.59 |
200 | 1.32 | 0.61 | 0.80 | 1.30 | 1.06 | 0.72 |
280 | 1.15 | 0.45 | 1.06 | 1.34 | 1.11 | 0.71 |
330 | 1.38 | 0.86 | 1.19 | 1.85 | 1.28 | 1.02 |
500 | 1.34 | 1.01 | 0.98 | 2.14 | 1.16 | 1.19 |
. | POC Flux SS (mmol C m−2 d−1)a . | POC Flux NSS (mmol C m−2 d−1)b . | Average POC Flux (mmol C m−2 d−1)c . | |||
---|---|---|---|---|---|---|
. | 5–51 μm . | 5–51 μm . | 5–51 μm . | |||
Depth (m) . | Average . | Uncert . | Average . | Uncert . | Average . | Uncert . |
50 | 6.77 | 2.65 | 4.19 | 2.22 | 5.48 | 1.73 |
65 | 5.91 | 2.07 | 3.81 | 1.98 | 4.86 | 1.43 |
80 | 4.47 | 1.38 | 3.08 | 1.63 | 3.77 | 1.07 |
100 | 2.35 | 0.62 | 1.66 | 0.92 | 2.01 | 0.56 |
120 | 2.02 | 0.66 | 1.45 | 0.99 | 1.74 | 0.60 |
150 | 1.55 | 0.62 | 1.05 | 1.00 | 1.30 | 0.59 |
200 | 1.32 | 0.61 | 0.80 | 1.30 | 1.06 | 0.72 |
280 | 1.15 | 0.45 | 1.06 | 1.34 | 1.11 | 0.71 |
330 | 1.38 | 0.86 | 1.19 | 1.85 | 1.28 | 1.02 |
500 | 1.34 | 1.01 | 0.98 | 2.14 | 1.16 | 1.19 |
Abbreviations: POC = particulate organic carbon; Uncert = uncertainties.
aPOC fluxes were calculated by multiplying the C/Th ratio in 5–51 μm particles (2) by the cruise-average steady-state (SS) 234Th flux (Table 1) at a given depth; uncertainties, by propagating the C/Th ratio and 234Th flux uncertainties.
bPOC fluxes were calculated by multiplying the C/Th ratio in 5–51 μm particles (Table 2) by the non-steady-state (NSS) 234Th flux (Table 1) at a given depth; uncertainties, by propagating the C/Th ratio and 234Th flux uncertainties.
cAverage POC flux was calculated by averaging the SS and NSS POC flux estimates using the C/Th ratio in 5–51 μm particles; uncertainties, by propagating the uncertainties of the SS and NSS POC flux estimates.
We note that the C/Th ratio of mid-sized particles collected using in situ pumps was assumed to be most representative of the population of sinking particles. However, for completeness, we also include POC fluxes derived using the C/Th ratios measured in other size classes and using both SS and NSS 234Th flux models (Table S3). For example, the use of C/Th ratios measured in small-sized particles resulted in a smaller decline in POC flux with depth, with a POC flux maximum at 100 m of 3.5 mmol C m−2 d−1 that decreased to 3.0 mmol C m−2 d−1 by the PPZ (120 m bin). POC fluxes from 200 to 500 m converge to 1.1–1.3 mmol C m−2 d−1, the same as those calculated using the mid-sized particles. Using the large-particle C/Th ratio resulted in calculated POC fluxes from 50 to 500 m that ranged from 2.1 to 3.3 mmol C m−2 d−1 with no discernible depth patterns.
5. Discussion
5.1. POC flux estimates during EXPORTS at Ocean Station Papa
The POC flux derived from our extensive 234Th survey, using C/Th ratios in mid-sized particles and the average of SS and NSS 234Th models, was 5.5 mmol C m−2 d−1 at 50 m decreasing to 1.7 mmol C m−2 d−1 by 120 m (Figure 6 and Table 3), which was near the depth of the PPZ (117 ± 5 m) and Ez0.1 (118 ± 9 m). While individual 234Th depth profiles varied, the calculated differences in the spatial distributions of POC fluxes at 120 m throughout the EXPORTS study area were relatively small (Figure 7). These results were consistent with spatially coherent temperature and salinity measurements that indicate similar water masses were sampled during the cruise (Figure S4). In contrast, the increase in 234Th activity by Epoch 3, consistent with a decline in particle fluxes, aligned with a temporal increase in chlorophyll concentrations (Figure 1) that is consistent with a decrease in C export. Spatially, while 234Th flux profiles showed considerable similarity with increasing depth (Figure 2), closer examination of individual flux profiles showed that relatively distinct differences still existed, with some stations showing a greater decrease with depth (a larger 234Th excess) and others continuing to increase or remain constant below the upper 80–120 m. These results imply that processes below the Ez, or the characteristics of the particles exiting the Ez, such as their lability and sinking rates, were variable and led to differing degrees of POC flux attenuation (see Figure S5 for spatial maps of POC flux at 50, 100, 150, and 200 m).
Sinking particle fluxes were also measured during the EXPORTS campaign using surface-tethered traps and neutrally buoyant sediment traps deployed for 3–5 days per epoch at 5 depths spanning from 100 to 500 m (Estapa et al., 2020). Similar to the largest size fraction, and consistent with many prior sediment trap studies, swimmers increased the measured POC flux in the traps. Large amounts of swimmer C were removed using a 335-µm screen and manual picking under a microscope, though some small swimmers were likely to have passed through the mesh. In the trap samples, this additional swimmer contribution to POC flux was estimated statistically through comparison to mass, 234Th and biogenic silica fluxes, and removed from the trap fluxes reported here (Estapa et al., 2020). The swimmer-corrected sediment trap POC fluxes ranged from 0.45 to 2.8 mmol C m−2 d−1 from 100 to 200 m, lower than 234Th-derived POC fluxes (regardless of C/Th size class) at the same depths (Figure S6 and Table S3). Unlike 234Th-derived estimates, sediment trap POC fluxes were characterized by an increase in POC flux in Epoch 3 relative to Epochs 1 and 2 (from 0.9 to 2.3 mmol C m−2 d−1 at 100 m). This increase in the traps and decrease in the 234Th-derived POC fluxes (i.e., increase in 234Th activities) may be due to the temporal differences in the integration times of the 2 methods (3–5 days for traps versus 3 weeks for 234Th). The higher Epoch 3 trap flux is more similar to our best estimate of the 234Th-derived POC fluxes at 100 m. A more detailed consideration of trap fluxes and a comparison to water column-derived 234Th fluxes will be presented in the study of Estapa et al. (2020).
5.2. Prior studies of POC flux at Ocean Station Papa
A few previous studies have examined 234Th-derived POC fluxes at Station P. Charette et al. (1999) collected 234Th profiles in February, May, and twice in August 1996 along the Canadian “Line P” repeat sampling program that transects from the British Columbian coast to Station P. Their 234Th activities at Station P in August overlap with our data more than 25 years later (Figure 8a); their fluxes of 1,500 dpm m−2 d−1 at approximately 200 m are similar to those measured at the same depths during EXPORTS using a similar SS model (1,300 ± 600 dpm m−2 d−1). This coherence is particularly notable, given that Charette et al. (1999) used a slightly different method. A large volume pumping system was used to scavenge dissolved 234Th onto Mn impregnated cartridges, and particulate 234Th was collected using a 1-µm cartridge filter, though POC was measured on a separate quartz filter, identical to the 1 µm QMA used in this study. Despite these different techniques, their >1 µm C/Th data also overlap with our 5–51 µm size class and showed a clear decrease with depth down to 200–300 m, where C/Th ratios varied from 1.0 to 1.2 μmol dpm−1 (Figure 8b). Charette et al. (1999) reported the maximum POC flux of 8.1 mmol C m−2 d−1 at 40 m for their August sampling at Station P which is a bit higher but similar to our estimates (Table 4).
Data Source . | POC Export Flux (mmol C m−2 d−1) by Depth . | Comments . | |||
---|---|---|---|---|---|
50 m . | 100 m . | 200 m . | 500 m . | ||
234Th-based | |||||
This study | 5.5 ± 1.7 | 2.01 ± 0.56 | 1.06 ± 0.72 | 1.16 ± 1.19 | Aug–Sept, avg. of 61 profiles; >5–51 µm C/Th (Table 3) |
Charette et al. (1999) | 8.1 ± 3.4a | —b | —b | —b | August single profile; >1 µm C/Th |
Kawakami et al. (2010) | —b | 2.3 ± 0.5 | —b | —b | October, single profile; >0.7 μm C/Th |
Mackinson et al. (2015) | 4.1 ± 0.6 | 3.0 ± 0.1 | 2.0 ± 0.8 | —b | June 2011 and 2012; recalculated here from 3 profiles; >53 μm C/Th |
Trap-based | |||||
Wong et al. (2002) | 10.4 ± 3.1 | 3.9 ± 2.5 | 4.0 ± 1.2 | 1.9 ± 0.4 | July–Sept avg.; includes unpublished data (F Whitney, M Robert, pers. comm.) |
Timothy et al. (2013) | —b | —b | 1.3 | —b | Deep moored time-series traps, annual avg. |
Estapa et al. (2020) | —b | 1.38 ± 0.77 | 0.76 ± 0.24 | 0.86 ± 0.43 | Aug–Sept EXPORTS cruise, avg. of 3 casts |
Mackinson et al. (2015) | 11.8 ± 1.7 | 6.2 ± 0.3 | 6.8 ±1.8 | —b | June 2011 and 2012, avg. of 2 deployments, no poison |
Data Source . | POC Export Flux (mmol C m−2 d−1) by Depth . | Comments . | |||
---|---|---|---|---|---|
50 m . | 100 m . | 200 m . | 500 m . | ||
234Th-based | |||||
This study | 5.5 ± 1.7 | 2.01 ± 0.56 | 1.06 ± 0.72 | 1.16 ± 1.19 | Aug–Sept, avg. of 61 profiles; >5–51 µm C/Th (Table 3) |
Charette et al. (1999) | 8.1 ± 3.4a | —b | —b | —b | August single profile; >1 µm C/Th |
Kawakami et al. (2010) | —b | 2.3 ± 0.5 | —b | —b | October, single profile; >0.7 μm C/Th |
Mackinson et al. (2015) | 4.1 ± 0.6 | 3.0 ± 0.1 | 2.0 ± 0.8 | —b | June 2011 and 2012; recalculated here from 3 profiles; >53 μm C/Th |
Trap-based | |||||
Wong et al. (2002) | 10.4 ± 3.1 | 3.9 ± 2.5 | 4.0 ± 1.2 | 1.9 ± 0.4 | July–Sept avg.; includes unpublished data (F Whitney, M Robert, pers. comm.) |
Timothy et al. (2013) | —b | —b | 1.3 | —b | Deep moored time-series traps, annual avg. |
Estapa et al. (2020) | —b | 1.38 ± 0.77 | 0.76 ± 0.24 | 0.86 ± 0.43 | Aug–Sept EXPORTS cruise, avg. of 3 casts |
Mackinson et al. (2015) | 11.8 ± 1.7 | 6.2 ± 0.3 | 6.8 ±1.8 | —b | June 2011 and 2012, avg. of 2 deployments, no poison |
Abbreviations: Avg. = average; EXPORTS = EXport Processes in the Ocean from RemoTe Sensing; POC = particulate organic carbon.
aPOC flux measured at 40 m which was considered to be the base of the Ez0.1 (Charette et al., 1999).
bIndicates data not available.
Kawakami et al. (2010) also measured 234Th at Station P in October 2005. Again, they reported nearly identical 234Th activities (Figure 8a), 234Th fluxes, and the same low and decreasing C/Th ratios at depth (<2 μmol dpm−1) in size-fractionated particulate samples collected with 3- and 53-µm size filters, similar to our EXPORTS results (Figure 8b). Finally, a more recent study at Station P includes 234Th data and C/Th ratios from pumps in June for both 2011 and 2012 (Mackinson et al., 2015). Their data show good agreement as well for both 234Th (234Th fluxes within 25% of EXPORTS in upper 200 m) and similar C/Th ratios on their >10–53 µm and >53 µm filers (Figure 8b). All of these studies reported a decrease in C/Th with depth, which supports our use of the 5–51 µm C/Th ratios in this study. This broad agreement over 25 years of 234Th fluxes and C/Th ratios suggests that similar conditions were sampled at Station P.
In addition to 234Th measurements, sediment trap work has been pursued at Station P for decades, most notably the deep-moored conical sediment trap time series at 3,800 m initiated by Wong in 1982 and ending in 2006 as summarized by Timothy et al. (2013). A second trap at 1,000 m was added in 1983, and a third trap at 200 m was added in 1989. There was a concern that the fluxes at the shallower 200 m depth were too low due to low trap collection efficiencies, as evidenced, for example, in that study by an increase in biogenic silica flux with depth. Overall lower trapping efficiencies are common for moored traps at depths shallower than 1,500 m based upon 230Th and 231Pa budgets (Scholten et al., 2001; Yu et al., 2001; Buesseler et al., 2007). Timothy et al. (2013) reported an annual average POC flux of 1.3 mmol C m−2 d−1 at 200 m (Table 4), very similar to what we derived from 234Th at 200 m, yet their POC fluxes were considered a lower limit due to their concerns that these fluxes were due to low trapping efficiencies, though they did not have 230Th data to quantify this potential bias.
Perhaps more comparable to our 234Th data are a series of shorter, 2-to 5-day deployments between 1987 and 1998 of drifting traps as summarized in the study of Wong et al. (2002) and augmented with additional data through 2004 (F Whitney and M Robert, pers. comm. 10/15/2019). These traps were modeled on the VERTEX cylindrical design (Knauer et al., 1979) and were deployed more than 20 times at up to 10 depths between 50 and 1,000 m, mostly during the spring and summer with a few deployments in autumn and winter. They reported a POC flux as high as 11.3 mmol C m−2 d−1 at 50 m in summer, which decreased rapidly to 6.1 mmol C m−2 d−1 at 100 m and 4.4 mmol C m−2 d−1 by 200 m. Those POC fluxes are higher than the EXPORTS fluxes, though the depth trend is similar (Figure S6). Closer examination of only their July to September deployments shows that average POC fluxes were similar to the annual means, with 10.4 ± 3.1 mmol C m−2 d−1 at the 50 m maximum and 4.0 ± 1.2 mmol C m−2 d−1 by 200 m (Table 4 and Figure S6). Given our concerns with swimmer removal using 335-µm screens (Estapa et al., 2020), we hypothesize that one possible cause for the higher drifting sediment trap POC fluxes is that small swimmers were not removed by either the 500- or 1,000-µm screens that Wong et al. (2002) used to separate swimmers from passively sinking trap material. Additional evidence for swimmer contamination noted by Wong et al. (2002) is the high C/N ratio and a positive intercept of C versus N (positive C flux when N is zero).
While Mackinson et al. (2015) also used similar 3-day drifting traps as in the study of Wong et al. (2002), they did not poison their traps, which is not recommended (Buesseler et al., 2007). Their C flux results were higher than the other trap studies, for example, by a factor of 4–9 below 100 m compared to Estapa et al. (2020; Table 4). Whether the higher fluxes are due to the lack of poison, or perhaps higher flux conditions during the June bloom versus late summer period, is difficult to know. However, their trap fluxes are also higher than their own 234Th-derived C fluxes by a factor of 2–3 (Table 4).
A number of global assessments of POC flux using 234Th have been made (Le Moigne et al., 2013; Henson et al., 2019). Le Moigne et al. (2013) examined 234Th-derived POC fluxes at 100–150 m and noted that their global patterns and magnitudes were similar to those predicted in global models (Laws et al., 2000; Schlitzer, 2004). However, they cautioned that there were considerable data gaps in high export and high productivity settings, particularly in regions dominated by upwelling. Mean POC fluxes averaged 9 mmol C m−2 d−1, albeit across a huge range, from approximately 0 to >100 mmol C m−2 d−1. Within the biogeochemical Longhurst province that includes Station P (East Pacific Subarctic Gyre), POC fluxes averaged 4.1 ± 2.6 mmol C m−2 d−1, about double our results at these depths, although the compilation of Le Moigne et al. (2013) also included locations closer to the coast. Overall, our finding of a POC flux maximum during EXPORTS of 5.5 mmol C m−2 d−1 at 50 m, decreasing to a flux of 1–2 mmol C m−2 d−1 at the base of the Ez and below, places the Station P site as one of the more modest C flux settings.
5.3. Biological carbon pump efficiencies at Ocean Station Papa
Marine ecosystems regulate the global biogeochemical cycling of C through the uptake of inorganic C via photosynthesis and the transformation of this C into a variety of organic forms that are transported to depth via sinking particles, although additional export of organic materials can occur through a variety of other physical and biogeochemical processes as well (Boyd et al., 2019). Here, we tracked the net export of POC using 234Th and quantified the strength of the biological pump and the efficiency of sinking POC transport through the upper twilight zone where most of the remineralization occurs (Buesseler et al., 2020). Combined, these processes determine the overall efficiency of the biological C pump (BCP).
One of the simplest measures of the strength of the BCP is export efficiency or the ratio of export production (EP) to NPP. When EP is determined by sediment traps, it is often referred to as the e-ratio (Murray et al., 1989). When derived using the 234Th approach, as done here, EP is defined as the ThE ratio (Buesseler et al., 1998). Regardless of the method, the depth over which the EP comparison is made is critical, given the rapid changes in the BCP with increasing depth. Buesseler and Boyd (2009) suggested that an Ez-normalized depth is most appropriate for examining the BCP efficiencies because below the Ez, POC losses dominate as the formation of new POC via photosynthesis cannot occur. Defining the depth of the Ez by light or PAR sensors, however, does not work at night and is influenced by time of day and local variations in cloud cover. In contrast, sensors for fluorescence are more common and are not impacted by these issues. As such, using fluorescence as a real-time indicator of Ez is often a more practical way to define the PPZ. As introduced by Owens et al. (2015), the PPZ is defined as the depth where in situ fluorescence declines to 10% of the maximum signal measured in overlying waters. During EXPORTS, the depths of 0.1% PAR (Ez0.1), the PPZ, and where 234Th reaches equilibrium with 238U overlapped at 100–125 m (Figure 1). The depth of the Ez0.1 during EXPORTS can thus be estimated using the PPZ.
During EXPORTS, the average Ez ratio (POC flux/NPP) was 13% ± 5% when using the best estimate of POC fluxes at 120 m (Ez0.1) and average integrated NPP determined from 14C methods using onboard incubations (the average of 11 profiles collected throughout the cruise is 13 ± 2 mmol C m−2 d−1; J Fox, pers. comm. 1/18/2020; 14C data available at https://seabass.gsfc.nasa.gov/archive/OSU/behrenfeld/EXPORTS/EXPORTSNP/archive). We reemphasize the strong dependence of the Ez ratio on the chosen POC flux depth. Here, Ez ratios increased to 29% if assessed at 80 m or at the average depth of Ez1.0. As noted above, the 234Th results and C/Th ratios measured during EXPORTS are similar to the much earlier data presented in Charette et al. (1999). However, Charette et al. also reported much higher rates of NPP, 58 and 85 mmol C m−2 d−1 in May and August, respectively, which results in somewhat lower estimates of the Ez ratio of 3% and 14% (for May and Aug, respectively; Buesseler and Boyd, 2009). Likewise using the 234Th-derived C flux estimates and a higher NPP in June (50–105 mmol C m−2 d−1) resulted in lower Ez ratios of 3%–6% (Ez1.0) in the study by Mackinson et al. (2015) conducted in June (Table 3 in Mackinson et al., 2015).
To assess POC flux attenuation below Ez, we quantified a transfer efficiency, or T100, which is defined as the ratio of the POC flux 100 m below the Ez to the POC flux at the base of the Ez (using PPZ or 0.1% PAR; Buesseler and Boyd, 2009). Using the average POC flux, T100 averages 61% using the 120- and 220-m POC fluxes (220-m flux is approximated from the estimated POC flux at 200 m; Table 3). In other words, about 40% of the POC flux is remineralized in the upper 100 m of the twilight zone. If we used the POC flux peak at 50 m (5.5 mmol C m−2 d−1), T100 would be significantly lower (24%). As stated earlier, the Ez ratio at Station P (14%) reported by Charette et al. (1999) was similar to EXPORTS in August; however, their POC flux attenuation was higher, with a T100 around 30% (70% lost). Although their normalization depth and local Ez0.1 were reported to be at 40 m during their August sampling, in hindsight, the phytoplankton dynamics at Station P (Boyd and Harrison, 1999) suggest that their Ez0.1 was likely much deeper. In effect, their higher POC attenuation reflects a shallower portion of the POC flux profiles and thus an enhanced attenuation bias, as discussed in the study of Buesseler et al. (2020). Mackinson et al. (2015) reported an Ez1.0 depth of 60–85 m in June; extrapolating their average 234Th-derived C fluxes to 75 and 175 m (from Table 4) resulted in a T100 closer to 80% or less POC flux attenuation in June.
During EXPORTS, sediment traps show a similar gradient in POC flux below Ez0.1 as determined using 234Th, although the 234Th-derived gradient was steeper above the shallowest trap flux assessment at 100 m (Figure S6; Estapa et al., 2020). The trap fluxes would result in a cruise-averaged Ez ratio of 10 ± 6% at a nominal 100-m depth (i.e., bracketed by Ez1.0 and Ez0.1), and a T100 of 55 ± 35% between the nominal 100- and 200-m trap depths (Estapa et al., 2020). These results are in good agreement with the 234Th-derived POC fluxes presented here. Earlier Station P moored trap data (Timothy et al., 2013) resulted in an EP/NPP ratio of only 3% at 200 m; however, this ratio is derived from an average NPP of 40 mmol C m−2 d−1; hence, the export ratio is lower due to higher NPP, not lower fluxes. These moored time-series flux data were most useful for observing annual patterns of export. Climatologies of the 200- and 1,000-m Station P sediment trap time-series presented in the study of Timothy et al. (2013) suggest that the August to September time frame is one of the strongly decreasing fluxes, especially at 200 m. The trend toward higher 234Th activities in Epoch 3 during EXPORTS may therefore reflect a seasonal shift. Interestingly, they noted that by the time the sinking particles reached 3,800 m, POC fluxes were 1.6 times higher than the average measured in deep moored traps in the Sargasso Sea. Timothy et al. (2013) attributed the higher deepwater fluxes to higher biogenic silica and carbonate contents of Station P fluxes that facilitated more efficient transport, at least deeper in the water column. These higher BCP efficiencies were not observed in EXPORTS in the upper twilight zone.
Globally, the Ez ratio varies from a few percentage to >50% (Buesseler and Boyd, 2009; Buesseler et al., 2020). Thus, an Ez ratio of 13% (from 234Th) was on the lower end but not as low as some oligotrophic settings. With this Ez ratio and T100 around 61%, the BCP at Station P is relatively inefficient, with only 8% of the organic carbon fixed in the Ez, reaching a depth of 100 m below the Ez. This overall efficiency was slightly higher than that calculated from data of Charette et al. (1999) for conditions at Station P in August (4%) and similar to what has been found in a reassessment of BCP during late bloom conditions in the northwest Pacific during VERTIGO and in summer conditions at the UK Porcupine Abyssal Plain site in the northeast Atlantic (Buesseler et al., 2020). Low overall BCP efficiencies in the summer and late summer in these comparable settings are seasonal in nature. For example, there is much higher export during spring bloom conditions in the northeast Atlantic, where the overall BCP efficiency can approach 50% during the demise of the local diatom spring bloom.
The overall BCP efficiencies here are similar to, or a bit higher than, those measured in the subtropical gyre time-series sites, Station ALOHA (A Long-Term Oligotrophic Habitat Assessment; 5% using 234Th; Buesseler and Boyd, 2009) and Bermuda Atlantic Time-series Study (BATS; 1%–5% using traps; Buesseler et al., 2020). At these sites, both the 234Th fluxes and C/Th ratios were low. Interestingly at Station P, the late summer average 234Th flux (1,230 dpm m−2 d−1 at 120 m) was higher than average 234Th fluxes at Station ALOHA (annual average = 765 dpm m−2 d−1 from 1999 to 2000; Benitez-Nelson et al., 2001a) and BATS (annual average = 600 dpm m−2 d−1 for 1993–1995; Sweeney et al., 2003). Lower C/Th ratios at Station P during EXPORTS, however, resulted in a low POC flux, given the higher relative 234Th flux. Low particulate C/Th ratios are hypothesized to be due to some combination of efficient recycling of organic matter and/or smaller particle sizes (Buesseler et al., 2006). In contrast, the higher C/Th ratios (>10 μmol dpm−1) typical of the Southern Ocean and phytoplankton blooms are often associated with larger diatoms species and high POC export (Buesseler et al., 2001, 2006; Puigcorbé et al., 2020). These larger size classes of plankton were not dominant at Station P during EXPORTS, nor are they common at other sites with low BCP efficiencies (Puigcorbé et al., 2015). The conditions encountered during EXPORTS were similar to what Henson et al. (2019) referred to as a low NPP (<85 mmol C m−2 d−1) and moderate export efficiency regime (EP/NPP 2%–20% at 100 m). These regimes are characterized by ecosystems with small primary producers that are tightly coupled to recycling by zooplankton and bacteria. They are also characteristic of subpolar waters such as Station P and many sites in the oligotrophic and equatorial open ocean of the Atlantic.
6. Conclusions
We have used an extensive set of 234Th measurements to quantify the spatial (horizontal scales of <10 km and vertical scales of tens of meters) and temporal (days to weeks) variability in export fluxes associated with sinking particles and their flux attenuation below the Ez. Our best estimate of the maximum flux of POC, using C/Th ratios measured on 5–51 µm particles and both SS and NSS models, was 5.5 mmol C m−2 d−1 at 50 m, decreasing rapidly with depth within the euphotic zone. This attenuation resulted in < 2 mmol C m−2 d−1 reaching the depth of the Ez0.1, which was near the base of the PPZ (117 m), where the PPZ defined the transition between shallower autotrophic production of particulate organic matter and the twilight zone where heterotrophic and physical processes reduce POC fluxes. Changes with time and space were small, although a statistically significant increase in 234Th activities during the sampling period was observed in the upper 60 m. When compared to NPP, average export efficiencies were relatively low (13%) and characteristic of a modest biological carbon pump. Interestingly, these latest results agreed with prior 234Th studies within the past 25 years at Station P, suggesting similar conditions for BCP efficiencies in the late summer at this site.
These results are the first comprehensive assessment of POC flux from the first EXPORTS campaign and provide a foundation for comparisons with future efforts that focus on the larger suite of EXPORTS measurements that were made during the study, including spatial variability in plankton community structure, particle size and composition, and other biological and physical processes. In addition, the downward flux of other elements and compounds, including particulate nitrogen, biogenic silica, phosphorus, and inorganic C, will be calculated for comparison to POC, using their ratios to 234Th on pump-collected samples and other estimates of net community production. Combined, these results will allow the development of mechanistic and predictive models of POC flux and the inclusion of other C loss terms associated with physical mixing transporting nonsinking POC and DOM to depth, and POC loss associated with the diel and seasonal vertical migrations of zooplankton (e.g., Boyd et al., 2019). These models are of critical importance if we are to better understand and predict the influence of marine ecosystems on the storage, transformation, and fate of carbon and associated elements in the Earth’s biosphere.
Data Accessibility Statement
Buesseler, KO, Benitez-Nelson, CR, Resplandy, L, Clevenger, S, Drysdale, J, Pike, S, Roca-Martí, M, Umhau, B, Wyatt, A (2019) EXPORTS-EXPORTSNP_thorium_234_survey, SeaBASS repository, DOI: 10.5067/SeaBASS/EXPORTS/DATA001, https://seabass.gsfc.nasa.gov/archive/WHOI/buesseler/EXPORTS/EXPORTSNP/archive/
Buesseler, KO, Benitez-Nelson, CR, Resplandy, L, Roca-Martí, M, Pike, S, Umhau, B, Horner, T, Masqué, P (2019) EXPORTS-EXPORTSNP_IN_SITU_PUMPS_SURVEY, SeaBASS repository, DOI: 10.5067/SeaBASS/EXPORTS/DATA001, https://seabass.gsfc.nasa.gov/archive/WHOI/buesseler/EXPORTS/EXPORTSNP/archive/
Estapa, M, Buesseler, KO, Omand, M, Durkin, C (2019) EXPORTS-EXPORTSNP_flux-POC_NBST, EXPORTS-EXPORTSNP_flux-POC_STT, EXPORTS-EXPORTSNP_flux-Th_234_NBST, EXPORTS-EXPORTSNP_flux-Th_234_STT, SeaBASS repository, DOI: 10.5067/SeaBASS/EXPORTS/DATA001, https://seabass.gsfc.nasa.gov/archive/SKIDMORE/estapa/EXPORTS/EXPORTSNP/archive
Supplemental files
The supplemental files for this article can be as follows:
Text S1. Figures S1-S6. Tables S1-S3. PDF.
Acknowledgments
EXPORTS would not have been possible without the conceptual and organizational leadership of Dave Siegel and the many participants in the planning process. On these cruises, we are gratefully appreciative of the captain and crew of the R/V Sally Ride and scientists on board as well as other EXPORTS scientists and support from the R/V Roger Revelle. We also want to thank James Fox, Michael Behrenfeld, and their group for providing the 14C primary production data; Debbie Steinberg and her group for providing samples of plankton from the nets; Colleen Durkin, Melissa Omand, Alyson Santoro, and Pat Kelly for their contributions to the sediment trapping program; and the Canadian scientists who conducted the prior trap studies at Ocean Station Papa, including Frank Whitney and Marie Robert, who provided results for our consideration of the drifting trap fluxes beyond which had been published previously. Finally, we thank the editors of Elementa for their efforts and for finding two very constructive reviewers whose comments helped greatly in improving the revised manuscript.
Author contributions
Contributed to conception and design: KOB, CBN, and LR.
Contributed to acquisition of data: KOB, CBN, MRM, AMW, SJC, JD, SP, MLE, and BPU.
Contributed to analysis and interpretation of data: KOB, CBN, MRM, AMW, LR, and MLE.
Drafted and/or revised the article: KOB, CBN, MRM, AMW, LR, SJC, MLE, and BPU.
Approved the submitted version for publication: KOB, CBN, MRM, AMW, LR, SJC, JD, MLE, SP, and BPU.
Funding
The authors would like to acknowledge support from the National Aeronautics and Space Administration (NASA) as part of the EXport Processes in the Ocean from RemoTe Sensing (EXPORTS) program awards 80NSSC17K0555 and 80NSSC17K0662; the Woods Hole Oceanographic Institution’s Ocean Twilight Zone study for KOB and MRM, and the National Science Foundation Graduate Research Fellowship Program (NSF-GRFP) for funding and support of AW.
References
How to cite this article: Buesseler, KO, Benitez-Nelson, CR, Roca-Martí, M, Wyatt, AM, Resplandy, L, Clevenger, SJ, Drysdale, JA, Estapa, ML, Pike, S, Umhau, BP. 2020. High-resolution spatial and temporal measurements of particulate organic carbon flux using thorium-234 in the northeast Pacific Ocean during the EXport Processes in the Ocean from RemoTe Sensing field campaign. Elem Sci Anth. 8: 1. DOI: https://doi.org/10.1525/elementa.2020.030.
Domain Editor-in-Chief: Jody W. Deming, University of Washington, WA, USA
Guest Editor: Deborah Steinberg, Virginia Institute of Marine Science, William & Mary, Williamsburg, United States
Knowledge Domain: Ocean Science
Part of an Elementa Special Feature: EXPORTS