An operational overview of the EXport Processes in the Ocean An operational overview of the EXport Processes in the Ocean from RemoTe Sensing (EXPORTS) Northeast Pacific field from RemoTe Sensing (EXPORTS) Northeast Pacific field deployment deployment

The goal of the EXport Processes in the Ocean from RemoTe Sensing (EXPORTS) field campaign is to develop a predictive understanding of the export, fate, and carbon cycle impacts of global ocean net primary production. To accomplish this goal, observations of export flux pathways, plankton community composition, food web processes, and optical, physical, and biogeochemical (BGC) properties are needed over a range of ecosystem states. Here we introduce the first EXPORTS field deployment to Ocean Station Papa in the Northeast Pacific Ocean during summer of 2 0 18, providing context for other papers in this special collection. The experiment was conducted with two ships: a Process Ship, focused on ecological rates, BGC fluxes, temporal changes in food web, and BGC and optical properties, that followed an instrumented Lagrangian float; and a Survey Ship that sampled BGC and optical properties in spatial patterns around the Process Ship. An array of autonomous underwater assets provided measurements over a range of spatial and temporal scales, and partnering programs and remote sensing observations provided additional observational context. The oceanographic setting was typical of late-summer conditions at Ocean Station Papa: a shallow mixed layer, strong vertical and weak horizontal gradients in hydrographic properties, sluggish sub-inertial currents, elevated macronutrient concentrations and low phytoplankton abundances. Although nutrient concentrations were consistent with previous observations, mixed layer chlorophyll was lower than typically observed, resulting in a deeper euphotic zone. Analyses of surface layer temperature and salinity found three distinct surface water types, allowing for diagnosis of whether observed changes were spatial or temporal.The 2 0 18 EXPORTS field deployment is among the most comprehensive biological pump studies ever conducted. A second deployment to the North Atlantic Ocean occurred in spring 2 0 21, which will be followed by focused work on data synthesis and modeling using the entire EXPORTS data set.


Introduction
The biological pump exports roughly 10 Pg of organic carbon from the surface ocean to depth each year (nearly equivalent to the global fossil fuel emission rate), sequestering carbon from the atmosphere and affecting the Earth's climate on time scales of months to millennia (e.g., Boyd et al., 2019;DeVries et al., 2019;Friedlingstein et al., 2019).Through these actions, the biological pump also increases the efficiency of the physical pump to sequester atmospheric CO 2 (e.g., Kwon et al., 2009).The biological pump has three components: sinking flux of particulate organic matter, active transport by daily and seasonally migrating animals, and physical transport due to circulation processes on scales from global meridional overturning to centimeter-scale turbulent motions.Thus, the biological pump is regulated by the coupling of ecological, physiological, behavioral, chemical, particle aggregation, and physical processes, making predicting future states of the biological pump one of the most challenging and, due to its role in the global carbon cycle, most important scientific challenges of our time.
The goal of the EXport Processes in the Ocean from RemoTe Sensing (EXPORTS) field campaign is to develop a predictive understanding of the export and fate of global ocean net primary production (NPP) and its implications for the Earth's carbon cycle in present and future climates (EXPORTS Writing Team, 2015;Siegel et al., 2016).Predicting the carbon cycle impacts of the fate of ocean NPP requires knowledge of (1) the flux of organic matter from the well-lit surface ocean, (2) the vertical attenuation of that downward flux within the ocean interior and its respiration back to dissolved CO 2 , and (3) the time scales and paths for these subsurface water masses to ventilate back to the surface.In particular, the vertical transmission of the downward export flux sets the degree to which the biological pump extends into the ocean interior-the deeper this penetration, the longer that organic carbon is sequestered from the atmosphere (e.g., Kwon et al., 2009;Boyd et al., 2019).
The same amount of NPP can lead to very different rates of organic carbon export and vertical flux transmission below the euphotic zone (Figure 1; Buesseler and Boyd, 2009;Buesseler et al., 2020b).A comparison of the fraction of NPP that is exported as vertical carbon flux from the euphotic zone (Ez-ratio) with the fraction of that flux that is transmitted to 100 m below the depth of the euphotic zone (T 100 ) clearly reveals that some regions of the global ocean are efficient at exporting NPP and transporting carbon to great depth (high Ez-ratio and T 100 values), such as the North Atlantic during its spring bloom.Other regions have lower Ez-ratios and T 100 values, including Ocean Station Papa (Station P) in the subarctic North Pacific (Buesseler and Boyd, 2009, using data collected in 1996 from Charette et al., 1999).The 2018 EXPORTS field campaign results (cross-hatched symbol in Figure 1) indicate similar Ezratio values to prior observations at Station P but with a higher T 100 .This difference in T 100 is thought to be due to the differing definitions of the euphotic zone depths used between studies (see Buesseler et al., 2020a, for further details).Assessing the biological carbon pump over a range of conditions will help ensure that the EXPORTS Field Campaign achieves the predictive understanding we seek, with the Northeast Pacific Ocean EXPORTS field deployment acting as a low efficiency endmember.
Three major export pathways link the upper ocean food web with the ocean interior, as depicted in the conceptual "wiring diagram" (Figure 2) from the EXPORTS 2015 Science Plan (EXPORTS Writing Team, 2015).These are (1) the gravitational settling of particulate organic matter, which may be composed of intact algal cells, detrital aggregates, and zooplankton byproducts; (2) the active transport of organic carbon to depth by vertically migrating metazoans; and (3) the advective transport of particulate and dissolved organic matter to depth by ocean circulation processes on a myriad of scales.The simultaneous assessment of the export pathways, along with the food web processes that create and alter this organic matter, is a major aim of the EXPORTS field campaign.
The underlying hypothesis for EXPORTS is that the export and carbon cycle fate of ocean NPP can be quantified knowing the characteristics of surface ocean ecosystems, thereby linking biotic carbon cycling processes to ecosystem properties that can be assessed remotely.To test this hypothesis, the EXPORTS Science Plan (Siegel et al., 2016) proposed three fundamental questions: 1. How do upper ocean ecosystem characteristics determine the vertical transfer of organic matter from the well-lit surface ocean? 2. What controls the efficiency of vertical transfer of organic matter below the well-lit surface ocean? 3. How can the knowledge gained from EX-PORTS be used to reduce uncertainties in contemporary and future estimates of the export and fate of upper ocean net primary production?
The answers to these questions provide a path toward a predictive understanding of the fates of global ocean NPP and its implications for the Earth's carbon cycle in present and future climates.
Here we present the scientific background and programmatic overview of the EXPORTS field campaign and a summary of the oceanographic conditions during the field deployment to Station P to provide context for other papers in this special collection.This introductory paper includes descriptions of campaign logistics, operational timelines, and a comparison of the EXPORTS observations with a climatology for Station P. We also share our plans for the EXPORTS program, as well as ongoing and future analyses of its data.

EXPORTS and the 2018 Northeast Pacific field deployment A brief history of the EXPORTS field campaign
The history of the EXPORTS field campaign stretches back nearly a decade.The original impetus for EXPORTS came from a scoping project funded by NASA to draft a Science Plan for a large-scale field campaign on the biological pump in 2013.After extensive peer and panel review, the EXPORTS Science Plan was completed (EXPORTS Writing Team, 2015;Siegel et al., 2016).In 2015, NASA formed the EXPORTS Science Definition Team to write an implementation plan that provided guidelines for how the EXPORTS field campaign would be conducted and provided a cost and risk assessment through an ensemble of deployment scenarios (EXPORTS Science Definition Team, 2016).The EXPORTS implementation plan called for a multi-scale sampling approach, combining two ships and autonomous underwater vehicles (AUVs) leveraging ongoing partnership programs, as depicted in Figure 3. NASA structured the EXPORTS field campaign as a multiyear effort with a "Pre-EXPORTS" modeling and datamining activity followed by a first phase with two major field programs and a second synthesis and modeling phase (Figure 4).The "Pre-EXPORTS" projects helped to plan the field campaign (e.g., Rousseaux and Gregg, 2017;Resplandy et al., 2019) and to create data sets mined from the literature (e.g., Bisson et al., 2018;Bisson et al., 2020;Kramer and Siegel, 2019) that can be useful for global syntheses (Table 1).These efforts were followed by Phase 1 of EXPORTS, where the Science Team openly competed to conduct field deployments in the Northeast Pacific in late summer 2018 and in the North Atlantic in spring 2020.Several individual projects also received support from the U.S. National Science Foundation (NSF), with all NSF principal investigators joining the EXPORTS Science Team (Table 1).The objectives of Phase 1 are to collect the necessary data to answer science questions 1 and 2, submit collected data to national repositories (SeaBASS and BCO-DMO) making them publicly available for all interested users, and conduct preliminary syntheses and modeling.
Following conclusion of the field components of Phase 1 and submission of collected data to public repositories, NASA plans to solicit proposals for the Phase 2 synthesis and modeling effort with a focus on further synthesis of the field results from Phase 1 and their application to improve global models of carbon cycle impacts of NPP fates.In March 2020, the COVID-19 pandemic led to the decision to postpone the North Atlantic deployment until the spring of 2021.

Site selection and its oceanographic characteristics
The EXPORTS NE Pacific deployment was conducted at Ocean Station Papa (Station P, nominally 50 N, 145 W).Station P is among the most well-studied openocean sites in the world (e.g., Tabata, 1965;Davis et al., 1981;Miller et al., 1991;D'Asaro et al., 1995;Wong et al. 1995;Boyd et al., 1999;Boyd et al., 2004;Whitney et al., 2005).In the 1940s, Ocean Weather Station Papa was established to provide weather information for forecasting.The first oceanographic samples were collected in 1949 by U.S. vessels, and standard hydrographic sampling commenced in 1956 by Canadian vessels and personnel (Freeland, 2007).Today, Station P is the terminus of Canada's Line P time-series sampling program, constituting one of the longest ongoing oceanographic time series, and is the site of instrumentation from the U.S. Ocean Observing Initiative (OOI), University of Washington Applied Physics Lab (UW-APL) waverider moorings, and the NOAA Pacific Marine Environmental Laboratory (NOAA-PMEL) air-sea interaction mooring (see Figure 5 for a map illustrating the locations of these moored assets, as well as the EXPORTS ship tracks).The Canadian Line P program reoccupied the EXPORTS site roughly 10 days after the EXPORTS ships left the Station P site.
Upper ocean stratification at Station P during the summer can be characterized by three main features: a mixed layer, thermocline, and halocline (Figure 6).The mixed layer typically extends between 20 and 30 m in the summer and deepens to roughly 100 m in winter.During the summer months, the mixed layer overlies a shallow seasonal thermocline that sits at the base of the mixed layer, whereas the wintertime mixed layer depth (MLD) marks the top of a seasonal pycnocline that is dominated by salinity stratification, the halocline.The temperature in the mixed layer undergoes a large seasonal cycle, peaking around 12 C to 14 C in the summer and cooling to 5 C to 6 C in the winter (Figure 7).The surface salinity shows smaller temporal variability spanning 32.4-32.6 in the surface layer with a larger variation of 0.6 in the thermocline layer (Pelland et al., 2017).The low-frequency circulation at Station P is dominated by the North Pacific Current, a slow, eastward-flowing current that separates the North Pacific Subtropical and Subpolar Gyres (e.g., Cummins and Freeland, 2007).Over shorter time scales, the one-dimensional balances described above can be disrupted by other processes that include spatial shifts in the North Pacific Current or the rare passage of mesoscale eddies (Freeland and Cummins, 2005;Jackson et al., 2006;Pelland et al., 2016).Eddy kinetic energy in the northeastern Pacific Ocean is among the lowest in the global ocean (e.g., Chelton et al., 2007;Xu et al., 2014), and coherent mesoscale circulation features are rarely seen, as Station P is noted to be an "eddy desert" (e.g., Chelton et al., 2011).In contrast, instantaneous currents are dominated by near-inertial wave motions due to the relatively strong wind forcing shallow mixed layers and weak mean and mesoscale currents at this site (e.g., D'Asaro et al., 1995).The dominance and relatively large horizontal scales of near-inertial motions makes Station P, in many regards, an ideal location to test one-dimensional models of mixed layer evolution (e.g., Denman and Miyake, 1973;Mellor and Durbin, 1975;Large et al., 1994).
Station P is also located in one of the three major ironlimited, high-nutrient, low-chlorophyll (HNLC) regions of the world ocean.Low iron availability limits phytoplankton production and thereby surface chlorophyll (Chl a) concentrations (e.g., Boyd et al., 2004).Previous to John Martin's demonstration of iron limitation in this area (Martin and Fitzwater, 1988), primary productivity was assumed to be limited by light and phytoplankton biomass to be kept low by zooplankton grazing (Parsons and Lalli, 1988).Martin's finding of iron limitation was at first controversial.However, shipboard carboy experiments at Station P validated that iron controlled primary productivity in late spring and summer (Boyd et al., 1996) and that mesozooplankton grazing could not fully consume the increased phytoplankton biomass when iron was added in carboys (Boyd et al., 1999).The in-situ SERIES iron enrichment experiment conducted in July of 2002 further verified iron limitation through the induction of a large phytoplankton bloom, composed primarily of diatoms, following iron enrichment (Boyd et al., 2004).Climatological surface nitrate concentrations are consistently in excess of 9 mmol L -1 throughout the region surrounding Station P (Figure 5; Whitney and Freeland, 1999), and  dissolved Fe concentrations in these surface waters are typically less than 10 pM (Martin et al., 1989;Schallenberg et al., 2017).Thus, special care is required to ensure that seawater for trace metal composition and biological rate measurements are collected and processed using trace metal clean (TMC) techniques, as was done during the 2018 EXPORTS cruise.
Temperature and salinity profiles from recent late-August cruises at Station P show a shallow well-mixed surface layer, a strong seasonal thermocline from roughly 30-60 m, and a strong stratifying halocline between roughly 100 and 130 m which defines the permanent pycnocline at this site (Figure 6A and B).In August, the mean (+standard deviation) MLD is 23 + 8 m (using   Values of AOU increase to nearly 300 mmol kg -1 at 500 m, illustrating the depletion of dissolved oxygen concentrations at depth due to remineralization (e.g., Bushinsky and Emerson, 2015;Pelland et al., 2018).Nitrate (þ nitrite) concentrations from Line P station P26 are roughly 8 mmol L -1 in the mixed layer and increase to approximately 18 mmol L -1 at 500 m (Figure 6D).Line P silicate concentrations are roughly 15 mmol L -1 in the upper layers and increase to nearly 30 mmol L -1 at a depth of 100 m.Concentrations of both nutrients found during EXPORTS follow these general patterns (Figure 6D and E).Mixed layer Chl a concentrations observed during EXPORTS (0.25 mg L -1 ), however, were lower than the Line P average (0.45 mg L -1 , Figure 6E).The nutrient and relatively low Chl a concentrations observed support the characterization of Station P as an HNLC site.
For August, the average fractional isolume depths for photosynthetically available radiation (PAR) fluxes from  insky and Emerson, 2015;Plant et al., 2016;Yang et al., 2018;Bif et al., 2019;Haskell et al., 2020).Mixed layer temperatures change on average by 6 C to 9 C from winter to summer, with the EXPORTS cruise happening during the warmest time of the year.Mixed layer mean concentrations of Chl a and particulate organic carbon tend to increase from winter through a peak in late summer with a considerable amount of interannual variation about these trends.BGC float-sensed nitrate concentrations typically decrease from March through September by about 7 mmol kg -1 , again with considerable interannual variability.BGC float observations during EXPORTS (Figure 7) occurred in slightly warmer mixed layer temperatures and lower POC, Chl a, and nitrate concentrations.

Operational details for the 2018 EXPORTS Northeast Pacific field deployment
The EXPORTS 2018 field deployment consisted of four major components (depicted in Figures 3 and 8 and listed in Table 2).First, the R/V Roger Revelle functioned as the Process Ship, sampling BGC stocks and fluxes, ecological abundances and rates, and optical properties following a Lagrangian float.Second, the R/V Sally Ride was the Survey Ship and characterized spatial variability about the Process Ship on scales from about 1 to 100 km.Third, a heterogeneous array of AUV platforms was deployed to set the spatial center of the sampling program, to provide horizontal spatial and high-temporal information, and to extend the temporal presence in the area.Last, a longterm sampling presence was created, tying the shipbased observations to climatically relevant time and space scales using BGC floats and partnerships with ongoing research programs.
The mission of the Process Ship, the R/V Revelle, was to constrain the pathways for organic carbon transformation and export, sampling in three distinct 8-day sampling  2).Operations also included the deployment of a multiple opening/closing net (MOCNESS) for assessing depth-discrete day-night differences in zooplankton populations, neutrally buoyant sediment traps (NBSTs), a surface-tethered sediment trap array (STT) that included in situ respiration traps, marine snow catchers for collecting sinking particles and assessing their fluxes, a Wirewalker (WW) for high resolution profiling of physical and bio-optical properties, a TMC towfish, and in situ and underway above-water bio-optics (Table 2).Plankton and microbial community structure were measured using a variety of techniques, including classical microscopy, flow cytometry, highthroughput microscopic imaging systems, metacommunity genomic sequencing, and gel trap-collected sinking particles.Phytoplankton rate measurements were made using TMC sampling methods (e.g., Measures et al., 2008;Mellet and Buck, 2020) with flow-through deckboard incubators screened for in situ light levels, with some incubators chilled to match the temperature below the mixed layer (including microzooplankton grazing, and mesozooplankton fecal pellet production rates).Additional rates were measured in temperature-and light-controlled laboratory incubators (e.g., bacterial production, microbial and zooplankton respiration, and particle sinking velocity).
Operations were conducted in three consecutive sampling cycles or "epochs" of 8-day duration (Figure 9A).Each epoch began with a positioning of the ship proximate to the Lagrangian float (LF) and deployment of the WW and surface-tethered sediment traps on day 1.NBSTs were subsequently deployed, with staggered recovery (3-5 days later) dependent on collection depth (Figure 9A).Shallow (0-150 m) predawn CTD casts (including TMC casts) every other day were used for primary production measurements conducted in deckboard incubators, and deeper (0-1,000 m) casts for multiple measurements occurred mid-day each day, and at night (for diel plankton composition and other comparisons) about every other day.A pair of day/night MOCNESS tows, centered on local noon and midnight, occurred on Days 2 and 7, with other net tows for live zooplankton experiments interspersed.Optics casts centered on the noon hour and occurred on most days, weather permitting.Marine snow catcher deployments and subsequent settling times occurred four days per epoch, which was also the maximum deployment time for the TMC towfish (Figure 9A).The spatial map (Figure 9B) illustrates the ship's track with trajectories of the sediment traps (NBST and STT) and WW deployed in the three epochs superimposed.Typically, NBSTs drifted 5-25 km depending on their deployment depth.Each day the ship embarked on a "poop run" to dump sewage holding tanks outside of a circle of 10-km radius extending from in-water assets (particularly sediment traps).These excursions can be seen in the ship tracks as large excursions from the central trajectory of the ship (Figure 9B) and proved useful to supplement spatial sampling from the ship's underway system.
The second component was the Survey Ship, the R/V Sally Ride.Its mission was to characterize the horizontal and vertical distribution of properties including phytoplankton, particulate and dissolved organic carbon,   2).Vertical profiles were collected from the CTDrosette system on both small (about 8-km spacing) and large (about 20-km spacing) scale grids for each epoch (Figure 10), with more intensive small-scale sampling in Epochs 1 and 3 and more intensive large-scale sampling in Epoch 2. Grids were centered on the location of the LF at the beginning of each epoch and were oriented so that the grid was perpendicular to the mean temperature and sea surface height fields.A complete suite of optical properties, including radiometric and inherent properties, were sampled at dawn and noon.Continuous underway samples for optical, biological, and selected chemical parameters were collected from an approximate 4-m depth with a flow-through system.Underway hyperspectral radiometry was collected during the day from the bow of the ship.Four pumping stations, with large volume pumps collecting particles at six to nine depths between 50 and A second task of the Survey Ship was to cross-calibrate sensors carried by the autonomous platforms and develop proxies for inferring biogeochemical parameters from sensor-based measurements.A total of 16 simultaneous ship CTD and autonomous asset profiles were collected (Figure 10): five for the LF, four for the Seaglider (SG), two for the BCG float, and three for the WW.These profiles were conducted to maximize spatial and temporal collocation of deliberate intercalibration casts by holding the target platform at the surface, moving the Survey Ship, which carried the designated reference sensors, into proximity (typically <300 m), and then simultaneously beginning profiles from both platforms.This approach facilitated alignment of all sensors against the reference units carried aboard the Survey Ship, allowing quantitative calculations that involve the entire array.
Sampling with autonomous platforms complemented the Process and Survey ship observations by resolving physical and biogeochemical variability at shorter temporal and spatial scales, and by expanding the sampling period to capture a broader range of physical and biological conditions (Table 2).A Seaglider (Eriksen et al., 2001), deployed on July 27 from the R/V Ride during the cruise, characterized mesoscale variability at the target site by occupying two repeats of a 25 by 25 km "bow tie" survey (Figure 11), at roughly 5-km horizontal resolution and 1-m vertical resolution, in the 2 weeks prior to the arrival of the EXPORTS ships.An instrumented Lagrangian float (D'Asaro, 2003), deployed from the Survey Ship on 14 August, defined a drifting reference frame used to target measurements through the three epochs and facilitate budget calculations just below the euphotic zone.The float was programmed to drift on an isopycnal near 100 m, with once-a-day surfacing.With the arrival of the ships and start of epoch-based sampling, the SG abandoned the "bowtie" survey to provide profiles at roughly 5-h resolution following the LF.A Wirewalker (Rainville and Pinkel, 2001), deployed near the LF at the start of each of the three epochs, profiled rapidly through the upper 500 m at intervals of 35-40 min.
Long-term sampling was conducted using two BGC profiling floats (Seabird Electronics Navis BGCi þ pH).These were deployed on 15 August near Ocean Station P and on 17 August in proximity of the LF.During the EX-PORTS field campaign BGC floats sampled in bursts, collecting three profiles, one from 2000-m depth and two subsequent profiles from 500-m depth, over 24 h, with a 3-day interval between 24-h bursts.After the EXPORTS campaign, one of the BGC floats was recovered due to malfunction and the other BGC float reverted to the standard 10-day Argo profiling interval.The SG and LF sampled until recovery on December 1, 2018, while the BGC float (BGC-0049; Dory) was still operating as of the submission of this manuscript (March 2021).Additional measurements were provided by two Slocum gliders and moorings operating nearby as part of the OOI Station Papa Array.The BGC floats, LF, SG, and WW provide sampling at temporal resolutions spanning minutes to days and persistence across two annual cycles (Figure 11).The array of autonomous platforms, combined with ship-based observations, also allow characterizing mesoscale and submesoscale spatial variability.
Quantitative interpretation of the sensor-based measurements provided by autonomous platforms relies on careful calibration and intercalibration across the entire sensor array.This task was accomplished by performing bulk, pre-deployment laboratory calibration of all sensors (ship-and autonomous platform-based), deliberate in situ intercalibration (detailed above), and post-deployment laboratory calibration.Unplanned close encounters between platforms offered additional opportunities for intercalibration.
The Canadian Line P program partnered with EXPORTS by adding additional stations and measurements to their regular hydrographic sampling program following the EX-PORTS cruise.The main objectives for this contribution were to: 1) provide an additional time point to the EX-PORTS time series, 2) add measurements unique to the Canadian team, 3) intercalibrate nutrient and highperformance liquid chromatography (HPLC) methods used on the Line P program with those used on EXPORTS, and 4) perform a small spatial survey of mesoscale variability.Sampling was conducted September 20-22, 2018, 12 days after the completion of the main EXPORTS shipboard sampling (Table 3; Figure 5).In addition to sampling at Station P and the NOAA PMEL mooring as part of the regular time series, intensive sampling with multiple casts was conducted at EX-C, the location of the subsurface Lagrangian float (within a ship length for the LF calibration cast), and reduced sampling at an additional four satellite stations, one of which (EX-F1) was the location of BGC float "Dory" (BGC-0049).Sampling included regular hydrographic parameters, productivity rate measurements, 234 Th, large volume pumping, trace metals, organic carbon, carbon system parameters, pigments, dissolved gases, and multi-omics (Table 3).
The large collection of assets as well as the need to intercalibrate across these platforms required both realtime coordination for all platforms and a formal situational awareness capability to inform near-term planning and support adaptive sampling.The location for nearly all elements of the EXPORTS sampling array were tracked in real time and made available to the ship and AUV operators every 5 min.The complexity of the field logistics also required that multiple lines of communication were shared between ship-and shore-based team members.Situational awareness reports were created daily describing weather and ocean conditions as well as the operations performed and planned on the Process and Survey Ships, and the AUV array.The daily reporting allowed progress and status of all platforms to be monitored, facilitated democratic decision-making, increased coordination among measurement efforts and partnering research programs, and supported adaptive sampling and risk mitigation.

Parameters measured and data availability
The EXPORTS Implementation Plan (EXPORTS Science Definition Team, 2016) provided a table of essential and advantageous measurements needed to answer the EXPORTS science questions (Section 7.2 in EXPORTS Science Definition Team, 2016).These needs were used to estimate the cost of the different scenarios evaluated in the EXPORTS Implementation Plan.Once the EXPORTS science team was selected, parameter working groups were formed to organize the measurement suite and cruise planning.The five working groups were: (1) Stocks, Proxies, and Context; (2) Optics; (3) Particle Characterization; (4) Rates; and (5) Export Pathways.Their goals were to account for all measurements to be made, reconcile any duplication in effort across the science team, assemble sampling and measurement protocols for each measurement, and create a master table of parameters (Table 4).The parameter working group report (EXPORTS Science Team, 2020) provides short methodological descriptions for each parameter and the master table of parameters and data availability can be accessed at https:// sites.google.com/view/oceanexports/home.Data generated during EXPORTS follow the NASA Earth Science Data and Information Policy (http://science.nasa.gov/earth-science/earth-science-data/data-informationpolicy/), which requires data to be available from a designated long-term public data repository within a year of collection.NASA-funded primary data products are archived at SeaWiFS Bio-optical Archive and Storage System (SeaBASS).All EXPORTS data are being archived under one digital object identifier (DOI: http://dx.doi.org/10.5067/SeaBASS/EXPORTS/DATA001) that further expands into the individual data subsets.NSF-funded primary data products are archived and distributed by the Biological and Chemical Oceanographic-Data Management Office (BCO-DMO).The BCO-DMO EXPORTS data-products webpage provides access to all the NSF-funded EXPORTS projects and links to the individual data set (https://www.bcodmo.org/program/757397),where each individual data set has a unique doi.Data sets that cannot be handled Rates of primary productivity (GPP, NPP, NCP) are determined via the incorporation of isotopes (stable and radioactive) over 6-and 24-h incubations to measure uptake of dissolved inorganic carbon, nitrogen and silicon; the production of biogenic gases (e.g., oxygen) to estimate GPP, PP and NCP.Respiration metrics are determined for bacterial and community respiration by O 2 drawdown assays and DOC remineralization assays.Zooplankton respiration rates are assessed by O 2 drawdown, analysis of electron transport system enzyme activity, and converting community composition and biomass measurements to community respiration.Also included are bacterial production, secondary production and grazing, aggregation and sinking rates, and nutrient uptake rates.

Export
Understand the mechanisms controlling the magnitude and attenuation of the primary export pathway fluxes Sinking particle rates are estimated using arrays of sediment traps, radionuclide mass balances, Marine Snow Catcher deployments and the numerical modeling of suspended and sinking particle organic carbon distributions.Assessments of the zooplankton-mediated migrant fluxes used day-night MOCNESS zooplankton biomass distributions along with assessments of zooplankton respiration and dissolved organic carbon and fecal production rates.Flux composition is assessed using high magnification microscope images of polyacrylamide gel collectors on sediment traps, genetic sequencing, image analyses on collected and in situ samples and contributions to bulk traps of amino acids.Processes regulating the vertical attenuation of sinking particle fluxes were determined using RESPIRE traps, Underwater Vision Profile images and the modeling of aggregate dynamics and coagulation.
increasing to the southeast (Figure 13).The EXPORTS assets were deployed to the northwest of an anticyclonic mesoscale feature creating a north-northeast current that was reflected in the drift of the Lagrangian float (Figure 11C).The SSH and SST patterns were similar for each epoch with a hint of an intensification of the mesoscale feature in time and an increase in SST during Epoch 2 and a cooling in Epoch 3. Minimal satellite ocean color imagery was available during the intensive operations period due to the persistent cloud cover (Figure 14).For August-September multiple satellite composites had better coverage, with most of the clear-sky views occurring in late September.A weak increase in satellite Chl a concentrations (about 0.1 mg m -3 ) appears associated with the location of the anti-cyclonic feature noted in the SSH fields, but even this feature is difficult to distinguish.
Time series of upper layer temperature and salinity from the R/V Revelle's underway flow-through system (Figure 15A) reflects many of the same general patterns observed at the NOAA PMEL mooring (Figure 12); that is, approximately stable SST and sea surface salinity relationships are found in Epoch 1, heating during Epoch 2, and change in surface water masses during Epoch 3. In particular, the influence of rain inputs is clear in this record during Epoch 3 (5 and 6 September).
Acoustic Doppler current profiler observations of mixed layer currents from the R/V Revelle show small sub-inertial currents, typically 2-3 km d -1 (1 km d -1 ¼ 1.157 cm s -1 ; Figure 15B and C).These low values support the notion that Station P has unusually low mesoscale kinetic energy compared with most other open ocean sites.During Epoch 1 and through the first part of Epoch 2, sub-inertial surface currents were comparatively large (about 4 km d -1 ) and flow was consistently to the north (Figure 15B).After that, sub-inertial currents were sluggish ( 2 km d -1 ) until the very end of Epoch 3. Mixed layer currents were dominated by supra-inertial frequency motions, where root mean square of the 15-min averaged current speeds (approximately 14.5 km d -1 ) are more than five times greater than the root mean square of the twoday averaged currents (Figure 15C).
Estimates of MLD were determined using CTD/rosette data from both ships as the first depth where the potential density exceeds the density at 50 m by 0.1 kg m -3 (Figure 15D).Values of MLD varied from 20 to 35 m with a mean depth of 31 + 4 m (n ¼ 227).Changes in time show a temporal pattern roughly similar to SST, where higher SSTs appear to correspond to shallower MLD determinations.
Daily integrated rates of incident PAR fluxes from the R/V Revelle varied nearly four-fold over the experiment (Figure 15E).Integrated daily incident PAR values were derived from the underway PAR data, collected every 15s on each ship, that have been cross-calibrated to the Compact-Optical Profiling System (C-OPS, Biospherical Instruments, Inc) deck sensor data before integration from local sunrise to sunset.Values of incident PAR varied 5fold from 10 to nearly 40 moles photons m -2 d -1 with the highest incident PAR fluxes found in the second epoch and the lowest in the third.
Bio-optical estimates of phytoplankton carbon (C phyto ) and Chl a concentrations were made using the R/V Revelle's underway system (Figure 15F).Particulate backscattering coefficient observations at 470 nm were used to estimate C phyto concentrations following Graff et al. (2015).Estimates of Chl a concentrations were made using particulate absorption line height at 676-nm determinations (Boss et al., 2013) and calibrated with coincident HPLC determinations of total Chl a concentrations collected during the campaign (Chl a ¼ 138.14 Â line_height 1.11 ; r 2 ¼ 0.67; n ¼ 43).In general, near-surface Chl a concentrations were low with mean concentrations of 0.21 mg L -1 (Figure 15F), considerably less than typical August values for Station P (Figure 7F).Variability is observed, as Chl a values ranged from approximately 0.10 to 0.35 mg L -1 , although changes in epoch mean Chl a values were small (approximately 0.01 mg L -1 ).The highest values are found near the middle of the first epoch (around day 232 or August 22), while the lowest values are found in the period between Epochs 1 and 2. Changes in the C phyto concentrations generally follow the behavior of Chl a, but their amplitude of changes are comparatively muted.Large short-duration departures in the underway Chl a concentrations are also observed in the data set at times (for example, see 23 August; Figure 15).These anomalies often coincide with "poop runs" for the R/V Revelle as the ship leaves its Lagrangian sampling frame.
Euphotic zone depths were estimated using depths of 1% incident PAR (Z eu ; Figure 15G).Values of Z eu were  determined from both ships using both PAR sensors mounted on the CTD/rosettes and hand-deployed spectroradiometer casts using C-OPS.These data were available from each ship, although C-OPS profiles were not conducted as often due to adverse weather conditions.Values of the vertical attenuation of PAR were calculated via robust linear regression over the depth range of 10-80 m and converted into estimates of the depth of the 1% incident PAR isolume (also assuming a 4% loss at the air/ water interface).Vertical attenuation coefficients for PAR were averaged if multiple estimates were available for a particular day.Values of Z eu ranged from 70 to 90 m with a good deal of variability with mean 1% PAR depth of 78 + 6 m (n ¼ 149; Figure 15G).Determinations of euphotic zone depth do not show a coherent pattern of change over time.The variation of temperature (T) and salinity (S) properties in the surface mixed layer, as measured from both the ship-based CTD casts and the ships' underway measurements, was small, spanning roughly 1 C (13.5 < T < 14.6 C) and 1 (32.24< S < 32.34) (Figure 16A).The temporal evolution of the measured surface properties in both space and time began with relatively salty surface properties (S ¼ 32.32, T ¼ 14 C) that freshen and cool during most of Epoch 1.As discussed in detail below, the transition to cooler water was sampled during a period of net surface warming, suggesting that these cooler waters had a distinct origin.From the end of Epoch 1 and throughout Epoch 2, surface salinity remained largely constant (32.29), but surface temperature increased by roughly 1 C between year day 232 and 241 (22 and 31 August).This rate of warming agrees well with the predicted surface warming based on surface flux data from the NOAA PMEL mooring (Figure 17).Around year day 244 (1 September) a distinct surface freshening was observed, bringing the salinity to values as low as 32.34, while temperature remained constant around 14.4 C. The magnitude and timing of this freshening signal is Based on these changes and an analysis of the surface flux data, we can identify three distinct groups of water properties in the mixed layer during the experiment (Figure 16A, B, and D).Groups 1 and 3 (blue and green, respectively, in Figure 16) are salty and fresh endmembers that were occupied at the beginning and end of the field program, respectively.The majority of the experiment, roughly year day 228-237 (August 16-25), sampled surface water in Group 2 (red in Figure 16) and covered the entire range of temperatures and a large fraction of the salinity values measured.However, near-surface water property changes during this period were well predicted by fluxes of surface heat, freshwater and momentum, suggesting that measurements made in the Group 2 period were largely in a Lagrangian framework following a fixed patch of surface water.The positions of the two ships with respect to time (year day) and to surface water type are shown in Figure 16C and D.
Isolating the advective changes from meteorological forcing (warming, freshening) in the mixed layer was accomplished by comparison with a one-dimensional Price-Weller-Pinkel model (PWP; Price et al., 1986).The PWP model uses the air-sea heat fluxes and wind stress to predict time-and depth-resolved variations in upper ocean S and T. The model was initialized with a CTD profile collected by the R/V Revelle at the beginning of each of the three water property groups defined above.Modeled T and S (colored lines, as per group, in Figure 17) was compared with T and S from the WW (black lines in Figure 17), both averaged over the upper 10 m, and with underway measurements from the RV Revelle (gray lines in Figure 17).Overall, the PWP model reproduced diel variations and the trend in surface warming, particularly that observed during Group 2. The model also successfully predicted the sharp freshening event that occurred on year day 245 (2 September during Epoch 3) after some prolonged precipitation.
Mean profiles of CTD sensor observations are shown in Figure 18A-D for T, S, potential density (s y ), and dissolved oxygen (O 2 ).The temperature profile (Figure 18A) shows a mixed layer of roughly 30 m and a strong seasonal thermocline, as was seen in the August Line P climatology for this site (Figure 7A).At this scale, few differences are apparent among the three epochs beyond Epoch 1 being slightly cooler than the other two from shows the classic two pycnoclines for this site and a sharp MLD with few inversions.Dissolved oxygen concentrations are elevated in the mixed layer and in the layer directly beneath it (to about 100 m) and decrease rapidly beneath that depth.The particulate beam attenuation coefficient at 700 nm (c p ) and the particulate backscattering coefficient at 700 nm (b bp ; Figure 18F and G) both account for the particle load within the water column, with values of c p representative of larger particles compared to b bp (e.g., Stramski and Kiefer, 1991;Zhang et al., 2020).Values of both c p and b bp are elevated in the upper 40 m and decrease to minimal values beneath that depth.The fractional reduction with depth for c p is much greater than for b bp , suggesting that a background of very small particles supported the b bp signal (Zhang et al., 2020) or that the contributions of organic-dominant particles vs. mineral-dominant particles changed with depth (Cetinic ét al., 2012).Mixed layer particle loads also tended to increase through the experiment.Values of AOU (Figure 18H) show supersaturation of dissolved oxygen concentrations in the upper 75 m of the water column, illustrating net autotrophy in the upper layers, particularly just beneath the mixed layer.Strongly negative AOU values are found at depth, illustrating the consumption of oxygen by respiration.
Nutrient profiles (Figure 19B and C) from water sampling conducted on the R/V Revelle show elevated concentrations (approximately 8 and 15 mmol kg -1 for NO 3 and SiO 4 , respectively) in the upper 40 m increasing with depth (roughly 20 and 25 mmol kg -1 for NO 3 and SiO 4 Profiles of Chl a concentration were determined from night-time fluorometer casts calibrated against total Chl a by HPLC (Figure 18E) or from discrete measurements (Figure 19D).Mixed layer values for the three epochs are roughly 0.2-0.3mg L -1 with the lowest mixed layer Chl a concentration occurring in Epoch 1, consistent with the underway bio-optical proxies from the R/V Revelle (Figure 15E).A subsurface Chl a maximum is found at about 70 m, and Chl a values decrease beneath that to background levels at about 120 m.Changes in mixed layer biogenic silica concentrations largely mirror the changes in Chl a concentrations (Figure 19G).
Both POC and dissolved organic carbon (DOC) were elevated in the mixed layer and decreased below a depth of about 40 m (Figures 19E and F).Values of DOC in the mixed layer are more than four times the POC values.The reduction with depth in DOC is roughly twice as much as the reduction in POC concentrations (roughly 5 vs. 2.5 mmol L -1 for DOC and POC, respectively).DOC concentrations show little change over time, although mixed layer POC concentrations increased slightly in Epoch 3 (similar to Chl a).Values of the absorption coefficient for chromophoric dissolved organic matter (CDOM) at 325 nm (Figure 19H) increase from approximately 0.13 m -1 in the surface layer to 0.2 m -1 at 40 and 60 m.These CDOM values approach the highest seen in the global open ocean (Nelson and Siegel, 2013).Also apparent is the lack of correspondence between the mean vertical profiles for CDOM and DOC, illustrating that CDOM is a poor proxy for DOC concentrations.
Though not all collocated, autonomous platforms provided a longer-term perspective on the evolution of the physical and biogeochemical properties after the ships departed.Because the Seaglider closely followed the Lagrangian float (Figure 11C), changes in Seagliderobserved bio-optical properties (Figure 11E) generally reflect growth and losses in phytoplankton rather than spatial heterogeneity.By mid-October 2018, a phytoplankton bloom had developed, largely in the surface mixed layer.Coinciding with this event, spikes in b bp increased and persisted throughout November and December.The mixed layer reached a maximum of around 100 m in early 2019 (Figure 11B-D), as would be expected from observations in previous years.Through the intensive bottle sampling, proxy development, and particle characterizations conducted during the cruise period, the variation in bio-optical properties observed by the autonomous platforms will be used to understand the seasonal changes in export for the year(s) following the cruise.Additional observations drawn from other sources, including satellite remote sensing, OOI Station Papa assets, the BCG float array, and the Line P program, will be used to expand understanding of export processes over an even broader range of biological and physical variability.

Next steps
We have presented background and goals of the EXPORTS program, highlighted the logistics and operations of the 2018 field campaign in the NE Pacific, and provided a synthesis of oceanographic context during the cruise to set the stage for manuscripts on various core data sets for this special collection in Elementa.Documentation of the accomplishments of the EXPORTS team during the 2018 NE Pacific field campaign, accompanied by the online availability of likely one of the most comprehensive coupled biogeochemical-ecological-oceanographic data sets collected to date (https://sites.google.com/view/oceanexports/home; EXPORTS Science Team, 2020), is critical for NASA's implementation of an EXPORTS Phase 2 program (see timeline in Figure 3).The papers in this special collection of Elementa are a first step towards answering the EXPORTS science questions and delivering on its promise to develop a predictive understanding of the biological pump.
As the first results focused on project level questions are published, the EXPORTS science team has been synthesizing their findings across projects.Over the past 2 years, five Synthesis Working Groups (SWGs) have been meeting monthly.The SWGs are focused on quantifying and understanding the following: (1) export pathways, (2) food web processes, (3) optical proxies, (4) time and space variability, and (5) BGC budgets (Table 5).Some of the early results from the SWGs have been presented here.For example, the mixed-layer water mass analysis and air-sea flux analyses (Figures 15 and 16) were the result of SWG 4. Other working groups have been assembling the flows of carbon through upper ocean and mesopelagic food webs (SWG 2) or assessing biogeochemical carbon budgets (SWG 5).These activities are critical for answering the EXPORTS science questions, but also help to ensure the quality and consistency among the entire measurement suite.The latter is especially critical given that many of the parameters sampled were either measured on different platforms or using different approaches or both.
With the global COVID-19 pandemic continuing at the time of submission of this manuscript, our view into the future was cloudy at best.The North Pacific experiment is only part of the EXPORTS story-an "endmember" in the Ez-ratio and T 100 space illustrated in Figure 1.While the North Atlantic experiment is needed to balance our understanding of how NPP carbon flows through the oceans of today and those of tomorrow, the breadth and complexity of the North Pacific EXPORTS deployment is destined to fuel new discoveries and pave the path for similar programs around the world.As this manuscript reached acceptance (May 2021), we can report that the North Atlantic EXPORTS experiment near the Porcupine Abyssal Plains Sustained Observatory (PAP-SO) site (https://projects.noc.ac.uk/pap/) in the northeast Atlantic Ocean is underway.
In July 2019, a group of researchers gathered in Southampton, UK, to form the Joint Exploration of the Twilight Zone Ocean Network (JETZON; https://jetzon.org).JET-ZON brought together scientists from 13 international projects tackling aspects of the role and function of the mesopelagic region of the oceans, the so-called twilight zone.A major goal of JETZON is to coordinate national studies of the biological carbon pump to optimize the understanding that these individual programs can provide by working together (Martin et al., 2020).To truly develop a predictive understanding of the fate of ocean NPP and its roles in the ocean carbon cycle, data from oceans around the world need to be synthesized.The hope is that by working together we can collect the global data sets needed to quantify the processes driving the ocean's biological pump and thus develop models that predict the present and future impacts of ocean ecological processes on the global carbon cycle.

Figure 1 .
Figure 1.Efficiency of the biological pump measured by export ratio and vertical flux transmission.On the y-axis is the export ratio (Ez-ratio), defined as the fraction of net primary production that is exported as vertical carbon flux from the euphotic zone (Ez); on the x-axis is the fraction of that flux that is transmitted to 100 m below the depth of the euphotic zone (T 100 ).Colors reflect different regional observations of the Ez-ratio and T 100 : orange are from Station P; green, from the North Atlantic bloom site (N Atl bloom); gray, from the Northwest Subarctic Pacific (NW Pacific); blue, from the Southern Ocean (So Ocean); black, from the Bermuda Atlantic Time Series (BATS); red, from the ALOHA site off Hawaii; and brown, from the Equatorial Pacific (EQPAC).Multiple occurrences of the same region illustrate seasonal differences for these sites.The figure is updated from Buesseler et al. (2020b) to include EXPORTS observations from Station P (orange-black hatched symbol); see text for further details and data sources.DOI: https://doi.org/10.1525/elementa.2020.00107.f1

Figure 2 .
Figure 2. Conceptual food web diagram from the EXPORTS Science Plan.Shown here is the conceptual food web diagram or "wiring diagram" from the EXPORTS Science Plan (EXPORTS Writing Team, 2015).The upper portion represents the euphotic zone (EZ) food web while the lower portion represents the mesopelagic zone, sometimes referred to as the twilight zone (TZ).These two regions are connected via the vertical flux pathways due to (A) sinking particles, (B) physical mixing of particles and dissolved organic carbon (DOC), and (C) the vertical migration of macrozooplankton and higher trophic organisms.The EZ food web denotes the input of net primary production by picophytoplankton and microplankton, the grazing by microzooplankton and microzooplankton and the recycling of some of the NPP energy by the microbial loop.Sinking particles can leave the EZ via (1) phytodetritus, (2) aggregates, and (3) fecal materials and zooplankton carcasses.The TZ food web consumes and recycles organic matter transferred to it via the three major pathways.The organic carbon consumption within the TZ food web controls the transmission of the sinking particle flux to depth.DOI: https://doi.org/10.1525/elementa.2020.00107.f2

Figure 3 .
Figure 3. Conceptual diagram of the sampling array for the EXPORTS Northeast Pacific Field Deployment.The sampling array for the EXPORTS Northeast Pacific field deployment is composed of both ship and autonomous sampling platforms.Two ships were deployed: a Process Ship, the R/V Roger Revelle; and a Survey Ship, the R/V Sally Ride.The autonomous sampling array was highly heterogeneous with an instrumented Lagrangian float that the Process Ship followed, a Seaglider, multiple neutrally buoyant sediment traps (NBSTs), a surface-tethered sediment trap array, an instrumented Wirewalker profiler, and two biogeochemical Argo profilers.The EXPORTS Northeast Pacific field deployment was conducted at Station P in the Northeast Pacific (near 45 N 145 W) which is a Global Node of the U.S. NSF's Ocean Observatories Initiative (OOI).DOI: https://doi.org/10.1525/elementa.2020.00107.f3

Figure 4 .
Figure 4. Timeline for the NASA EXPORTS Field Campaign.Shown here is the original NASA time line for the EXPORTS Field Campaign before postponement of the 2020 North Atlantic cruise due to the COVID-19 pandemic.The three boxes arranged vertically illustrate the planning, development and funding period (upper green), the Phase 1 field program (middle tan), and the planned Phase 2 data synthesis and modeling program (bottom blue).The NE Pacific cruise, the subject of this special collection of Elementa, is the first major field deployment for the Phase 1 program.The second field deployment to the North Atlantic was postponed due to the COVID-19 pandemic, but is underway at the time of manuscript acceptance (May 2021).The postponement of the second field deployment also pushes the start date for the Phase 2 program forward by at least a year.DOI: https://doi.org/10.1525/elementa.2020.00107.f4 of South Mississippi), Gray (Naval Research Laboratories), Guidi (LOV), Huot (Universite ´de Sherbrooke) Optically resolving size and composition distributions of particles in the dissolved-particulate continuum from 20 nm to 20 mm to improve the estimate of carbon flux NASA ST a Lead PIs underlined.b ST indicates science team; Pre indicates modeling and data mining activities prior to EXPORTS Phase 1.c Denotes current affiliation.

Figure 5 .
Figure 5. Site map for the EXPORTS Northeast Pacific field deployment.(A) Annual mean surface nitrate concentrations (color bar in mmol kg -1 ) from the 2013 World Ocean Atlas (https://www.nodc.noaa.gov/OC5/woa13/)along with ship tracks of the R/V Revelle (black line) and R/V Ride (gray line) in the boxed area, and Canadian Line P stations from their Fall 2018 cruise (Line P cruise 2018-40; blue circles).(B) Enlargement of inset in (A) showing locations of the NOAA PMEL air-sea interaction mooring (purple star), UW-APL waverider buoy (purple square), and NSF OOI subsurface moorings (purple diamonds) along with locations of proximal Line P stations from their Fall 2018 cruise (Line P cruise 2018-40; https://www.waterproperties.ca/linep/2018-040/;blue circles), as well as ship tracks of the R/V Revelle (black line) and R/V Ride (gray line).DOI: https://doi.org/10.1525/elementa.2020.00107.f5

Figure 7 .
Figure 7. Climatology of mixed layer properties from historical BGC float records near Station P. Mixed layer averages are shown for (A) temperature, (B) Chl a, (C) particulate organic carbon, and (D) nitrate.The thick black line with filled circles is the respective monthly climatology calculated as the average of the monthly means for the years 2008-2020, and the shaded gray envelope represents the associated standard error.Individual yearly records are denoted in the color bar; observations from 2018 and the EXPORTS campaign are shown by the red line and with yellow-filled red circles, respectively.DOI: https://doi.org/10.1525/elementa.2020.00107.f7

Figure 8 .
Figure 8. Experimental timeline for the 2018 EXPORTS Northeast Pacific field deployment.Colored vertical bars indicate the three sampling cycles or "epochs" of the 2018 EXPORTS field campaign, while gray horizontal bars indicate deployment periods for each of the platforms.Note that one of the biogeochemical floats (BGC-0049) was still operational at the time of submission of this manuscript.DOI: https://doi.org/10.1525/elementa.2020.00107.f8

Figure 9 .
Figure 9. Temporal and spatial operations of the R/V Revelle during the EXPORTS Northeast Pacific field deployment.(A) Broad categories of sampling events are marked by color, epoch boundaries are denoted by vertical dashed lines, and autonomous asset deployment periods are shown in gray.The epoch boundaries are given in UTC as Epoch 1, August 14 00:00 to August 23 09:00; Epoch 2, August 23 09:00 to August 31 09:00; and Epoch 3 August 31 09:00 to August 9 18:00.(B) Epochs and their spatial extent are delineated approximately by the autonomous assets deployed and recovered from the R/V Revelle (black dashed, purple, and blue lines) superimposed over the ship tracks (gray), with the first epoch in the south and the last epoch finishing in the most northeast of the EXPORTS experiment region.Locations of specific operations are color-coded to match those in (A).DOI: https://doi.org/10.1525/elementa.2020.00107.f9

Figure 10 .
Figure 10.Sampling by the R/V Ride during the EXPORTS Northeast Pacific field deployment.Panels illustrate the sample locations and types of samples taken from the R/V Ride during the EXPORTS Northeast Pacific field deployment.Columns show panels for each sampling epoch, while top row illustrates on large spatial scales and bottom row shows an expanded view.Symbols indicate the main focus of sampling for each station: black circles, CTD profiles with sensors only (no bottles fired); blue circles, thorium collection; red circles, dawn (only inherent optical properties, IOP) and noon optics suite (radiometric and IOP); yellow-orange, intercalibration casts with the autonomous assets, Wirewalker, and R/V Revelle; and purple, deep CTD casts.Open triangles designate stations for large volume pumping and black x's and plus signs represent locations of radiometric and IOP casts, respectively (when not coincident with optics suite stations).DOI: https://doi.org/10.1525/elementa.2020.00107.f10

Figure 11 .
Figure 11.Autonomous platform operations during the EXPORTS Northeast Pacific field deployment.Time series measurements are shown for (A) nitrate and (B) Chl a concentrations from biogeochemical (BGC) float 0949 profiles (August 16, 2018, to August 15, 2019), along with (C) a map of the positions of autonomous assets operating in the region before, during, and after (up to August 15, 2019) the intensive portion of the EXPORTS field campaign.Also shown are time series profiles of particulate backscattering from (D) BGC float (BGC-0049; August 16, 2018 to August 15, 2019), (E) Seaglider (July 27, 2018 to December 1, 2018), and (F) Wirewalker (August 23, 2018 to August 31, 2018).These panels show preliminary backscattering data (unitless) that have been scaled to allow for qualitative comparison between platforms.The black contour lines in each plot are the potential density surfaces, and the white line is an estimate of the mixed layer depth.Values in parentheses represent the number of profiles contributing to the time series shown for each platform.This number is the same for all BGC (panels A-C).Dots at 190 m in the BGC float panels indicate profile frequency.DOI: https://doi.org/10.1525/elementa.2020.00107.f11

Figure 13 .
Figure 13.Spatial maps of sea surface height and sea surface temperature during the EXPORTS sampling epochs.Spatial maps of (upper) absolute sea surface height from merged daily satellite altimetry and (lower) sea surface temperature (SST) for the three EXPORTS sampling epochs (left to right).Daily, near-real time, and gridded (0.25 resolution) absolute sea surface height were used to create these distributions (resources.marine.copernicus.eu/?option¼com_csw&view¼details&product_id¼SEALEVEL_GLO_PHY_L4_NRT_OBSERVATIONS_008_046).Optimally interpolated passive microwave (cloud penetrating) SST daily products (25-km resolution) were used as the SST products (www.remss.com/measurements/sea-surface-temperature/oisst-description).The black boxes indicate the sampling regions for the EXPORTS field deployment shown in Figure 6B.DOI: https://doi.org/10.1525/elementa.2020.00107.f13

Figure 14 .
Figure 14.Composite satellite ocean color Chl a concentration imagery during the EXPORTS field deployment.Composite log 10 -transformed Chl a concentrations for the ship sampling period at Station P (14 August to September 9, 2018) are shown in the upper and lower left panels (A, C) for the respective months of August and September 2018.Corresponding images for the number of valid days for each composite are shown in the two panels on the right (B, D).The composite Chl a imagery was created from daily merged observations (1/24 at approximate 4km resolution) from all available MODIS and VIIRS images (https://oceandata.sci.gsfc.nasa.gov/opendap/Merged_ATV/L3SMI/2018).The black boxes indicate the sampling regions for the EXPORTS field deployment shown in Figure 6B.DOI: https://doi.org/10.1525/elementa.2020.00107.f14

Figure 15 .
Figure 15.Oceanographic time series from the R/V Revelle and the R/V Ride.(A) Time series of temperature and salinity from the R/V Revelle underway sampling system (approximate 5-m depth intake).(B) Mixed layer currents (averaged between 20 and 35 m) from the R/V Revelle 150 kHz broadband ADCP system.The blue arrows show the low-passed (2-day filtered) mixed layer current magnitude and direction.(C) Root mean square (RMS) current speed from the R/V Revelle ADCP system.The gray line is the RMS of the 15-min averaged current speed, while the blue line is the RMS of the 2-day averaged currents.(D) mixed layer depth (MLD) estimates from the CTD profiles from the R/V Revelle (blue) and the R/V Ride (burnt orange).(E) Incident PAR flux (gray, left axis) and daily integrated PAR fluxes (blue, right axis).(F) Bio-optical estimates of phytoplankton carbon (black, left axis) and Chl a (green, right axis) concentrations from the R/V Revelle's underway optics system.(G) Euphotic zone depths are estimated as the depth of the 1% PAR surface.Estimates from the R/V Revelle (blue) and the R/V Ride (burnt orange) are shown as well as values from the CTDmounted PAR sensors (circles) and C-OPS spectroradiometry casts (stars).DOI: https://doi.org/10.1525/elementa.2020.00107.f15

Figure 16 .
Figure 16.Variations in surface mixed layer water masses.(A) Distribution of temperature ( C) and salinity from underway measurements (both ships, dots) as well as R/V Ride (circles) and R/V Revelle (squares) CTD casts.For all symbols, color indicates the time of the measurement (year day).Three water masses, outlined by blue, red, and green curves, delineate separate water mass properties (see text for details).(B) Time series of the three surfaces water masses encountered by the R/V Ride (circles) and R/V Revelle (squares).The water mass identification is discretized in 0.2-day bins.Maps of the R/V Ride and Revelle positions colored by (C) time (year day) and (D) water mass index are shown in the right-hand panels.DOI: https://doi.org/10.1525/elementa.2020.00107.f16

Figure 17 .
Figure 17.Time series of measured and model predicted temperature and salinity averaged over upper 10 m.Time series of (A) temperature and (B) salinity averaged over the upper 10 m from the PWP model (colored lines; see text for details), the Wirewalker (black lines) and the R/V Revelle underway system (gray lines).The PWP model was initialized at the start of each of the three water property groups.Dashed lines show the model predictions for a few days after each of these transitions, highlighting the departures between the 1D model (meteorological forcing only) and the observations (a combination of lateral advection and meteorological forcing).DOI: https://doi.org/10.1525/elementa.2020.00107.f17

Figure 18 .
Figure 18.Sensor CTD data profiles from both ships.(A) Temperature, (B) salinity, (C) potential density (s y ), (D) dissolved oxygen (O 2 ), (E) chlorophyll fluorescence (Fchl, night time profiles only), (F) particulate beam attenuation coefficient at 700 nm (c p ), (G) particulate backscatter at 700 nm (b bp ), and (H) apparent oxygen utilization (AOU) for the three epochs (blue, orange, and black correspond to Epochs 1, 2, and 3, respectively).Envelopes are 95% confidence intervals for the mean estimates at each depth bin.Data from the R/V Ride are combined with data from the R/V Revelle when both ships were within 25 km of each other.DOI: https://doi.org/10.1525/elementa.2020.00107.f18

Table 2 .
EXPORTS sampling platforms and measurements.DOI: https://doi.org/10.1525/elementa.2020.00107.t2 WirewalkerProfiles every 40 min from surface to about 500 m, instrumented with CTD and O 2 , optical backscatter, chlorophyll fluorescence, CDOM fluorescence, beam attenuation and PAR sensors; deployed at the start and recovered at the end of every epoch