A Synoptic Assessment of the Amazon River-Ocean Continuum during Boreal Autumn: From Physics to Plankton Communities and Carbon Flux

The Amazon generates the world's largest offshore river plume, which covers extensive areas of the tropical Atlantic. The data and samples in this study were obtained during the oceanographic cruise Camadas Finas III in October 2012 along the Amazon River-Ocean Continuum (AROC). The cruise occurred during boreal autumn, when the river plume reaches its maximum eastward extent. In this study, we examine the links between physics, biogeochemistry and plankton community structure along the AROC. Hydrographic results showed very different conditions, ranging from shallow well-mixed coastal waters to offshore areas, where low salinity Amazonian waters mix with open ocean waters. Nutrients, mainly NO3− and SiO2−, were highly depleted in coastal regions, and the magnitude of primary production was greater than that of respiration (negative apparent oxygen utilization). In terms of phytoplankton groups, diatoms dominated the region from the river mouth to the edge of the area affected by the North Brazil Current (NBC) retroflection (with chlorophyll a concentrations ranging from 0.02 to 0.94 mg m−3). The North Equatorial Counter Current (NECC) region, east of retroflection, is fully oligotrophic and the most representative groups are Cyanobacteria and dinoflagellates. Additionally, in this region, blooms of cyanophyte species were associated with diatoms and Mesozooplankton (copepods). A total of 178 zooplankton taxa were observed in this area, with Copepoda being the most diverse and abundant group. Two different zooplankton communities were identified: a low-diversity, high-abundance coastal community and a high-diversity, low-abundance oceanic community offshore. The CO2 fugacity (fCO2sw), calculated from total alkalinity (1,450 < TA < 2,394 μmol kg−1) and dissolved inorganic carbon (1,303 < DIC < 2,062 μmol kg−1) measurements, confirms that the Amazon River plume is a sink of atmospheric CO2 in areas with salinities <35 psu, whereas, in regions with salinities >35 and higher-intensity winds, the CO2 flux is reversed. Lower fCO2sw values were observed in the NECC area. The ΔfCO2 in this region was less than 5 μatm (−0.3 mmol m−2 d−1), while the ΔfCO2 in the coastal region was approximately 50 μatm (+3.7 mmol m−2 d−1). During the cruise, heterotrophic and autotrophic processes were observed and are indicative of the influences of terrestrial material and biological activity, respectively.

The Amazon generates the world's largest offshore river plume, which covers extensive areas of the tropical Atlantic. The data and samples in this study were obtained during the oceanographic cruise Camadas Finas III in October 2012 along the Amazon River-Ocean Continuum (AROC). The cruise occurred during boreal autumn, when the river plume reaches its maximum eastward extent. In this study, we examine the links between physics, biogeochemistry and plankton community structure along the AROC. Hydrographic results showed very different conditions, ranging from shallow well-mixed coastal waters to offshore areas, where low salinity Amazonian waters mix with open ocean waters. Nutrients, mainly NO − 3 and SiO − 2 , were highly depleted in coastal regions, and the magnitude of primary production was greater than that of respiration (negative apparent oxygen utilization). In terms of phytoplankton groups, diatoms dominated the region from the river mouth to the edge of the area affected by the North Brazil Current (NBC) retroflection (with chlorophyll a concentrations ranging from 0.02 to 0.94 mg m −3 ). The North Equatorial Counter Current (NECC) region, east of retroflection, is fully oligotrophic and the most representative groups are Cyanobacteria and dinoflagellates. Additionally, in this region, blooms of cyanophyte species were associated with diatoms and Mesozooplankton (copepods). A total of 178 zooplankton taxa were observed in this area, with Copepoda being the most diverse and abundant group. Two different zooplankton communities were identified: a low-diversity, high-abundance coastal community and a high-diversity, low-abundance oceanic community offshore. The CO 2 fugacity (fCO 2 sw), calculated from total alkalinity (1,450 < TA < 2,394 µmol kg −1 ) and dissolved inorganic carbon (1,303 < DIC < 2,062 µmol kg −1 ) measurements, confirms that the Amazon River plume is a sink of atmospheric CO 2 in areas with salinities <35 psu, whereas, in regions with salinities >35 and higher-intensity winds, the CO 2 flux is reversed. Lower fCO 2 sw values were observed in the NECC area.
The fCO 2 in this region was less than 5 µatm (−0.3 mmol m −2 d −1 ), while the fCO 2 in the coastal region was approximately 50 µatm (+3.7 mmol m −2 d −1 ). During the cruise, heterotrophic and autotrophic processes were observed and are indicative of the influences of terrestrial material and biological activity, respectively.

INTRODUCTION
Each late summer/autumn, the Amazon River plume covers ∼2 × 10 6 km 2 of the western tropical North Atlantic Ocean (WTNA) (DeMaster and Pope, 1996;Smith and Demaster, 1996;Ternon et al., 2000;Körtzinger, 2003;Cooley et al., 2007). The Amazon River has the greatest discharge of any global river and accounts for ∼20% of all of the riverine input to the oceans, more than the next seven largest rivers combined. The mean discharge of the Amazon River is approximately 150,000 m 3 s −1 and is responsible for approximately half of all the freshwater input into the tropical Atlantic (Baumgartner and Reichel, 1975;Yoo and Carton, 1990;Carton, 1991). This rate varies by 50% between a maximum in May-June and a minimum in November-December (Richey et al., 1989;Carton, 1991). The impacts of this plume on the WTNA include nutrients, microorganisms, and fresh water fluxes that contribute to enhanced biological activity and carbon sequestration over a million square kilometers of tropical ocean.
The Amazon River-Ocean Continuum (AROC) is an energetic region subjected to strong geophysical forcing's, including the Amazon River discharge, the North Brazil Current-North Equatorial Counter Current (NBC-NECC) system, macrotides and strong trade winds (Silva et al., 2005(Silva et al., , 2009(Silva et al., , 2010. The dynamics of the WTNA are also affected by the seasonal transposition of the Intertropical Convergence Zone (ITCZ). This region is also known as also an important location of heat exchange through a complicated system of currents and water masses around the equator (Stramma and Schott, 1999).
On the western edge of the WTNA, the northward migration of the ITCZ results in the retroflection of the NBC, which feeds the NECC with eastward-transported waters of the Amazon River plume (Richardson and Reverdin, 1987;Fonseca et al., 2004;Coles et al., 2013). The waters from the Amazon River and the increased rainfall caused by the presence of the ITCZ are the major sources of fresh water along the western edge of the WTNA.
These contributions are sensitive to the quantity and composition of the river discharge itself. The upper Amazon River is a source of CO 2 to the atmosphere (210 ± 60 Tg C year −1 ; Richey et al., 2002), supported by organic matter mineralization and carbon dioxide and organic matter export from flooded wetlands (Abril et al., 2014). When transported through the salt gradient, the high concentrations of nutrients in the river water are diluted via mixing with ocean water, favoring the growth of primary producers and decreasing the amount of associated organic carbon (Chen et al., 2012). The spread of Amazon waters in the tropical Atlantic is also known to support significant N 2 fixation through diatom-diazotroph associations, which represents the main carbon sequestration pathway within the plume (Subramaniam et al., 2008;Yeung et al., 2012).
The different factors that contribute to CO 2 undersaturation in the Amazon River plume and its seasonal variability remain poorly understood. Ternon et al. (2000) estimated that primary production within the plume could be responsible for approximately 30% of the observed CO 2 undersaturation. Additionally, Cooley et al. (2007) suggested that net primary production in the river plume would enhance the observed CO 2 undersaturation by a hundredfold. Elucidation of the physical and biological processes responsible for the modulation of the sea surface carbon dioxide fugacity (fCO 2 ) in the Amazon River plume and the WTNA may help better constrain the role of the tropical Atlantic in the global sea-air CO 2 exchange.
The main objective of this work is to improve our understanding of the processes and organisms responsible for carbon and nutrient cycling along a large-scale tropical riverocean continuum (Amazon River to offshore), focusing on the nearshore and offshore WTNA. This combined river-ocean continuum represents one of the largest environmental gradients on land and in the ocean in the world and stretches across thousands of km from the continental shelf to the middle of the Atlantic.

Hydrography and Currents
The dispersal of Amazon River water forms a brackish water plume that can exceed 10 6 km 2 , reaching latitudes as far from the river mouth as 30 • W (Coles et al., 2013) or even 25 • W when the North Equatorial Countercurrent (NECC) is strong (Lefèvre et al., 1998). Thus, the large areas of fresh sea surface waters (≤33 salinity) observed in the region are primarily due to the Amazon discharge.
The data and samples in this study were obtained during the oceanographic cruise Camadas Finas III (hereafter, CF3) aboard the research vessel NHo. Cruzeiro do Sul -H38 (DHN/Brazilian Navy). This cruise was performed during October 9th-31st, 2012, corresponding to the period when most of the Amazon plume is transported eastward and coinciding with the northernmost annual position of the Intertropical Convergence Zone (ITCZ).
The ship track encompassed the outer estuary portion, the alongshore northwestern NBC region, the NBC retroflection area and the eastward NECC plume transport to 38 • W ( Figure 1A).

Chemical Analysis, Apparent Oxygen Utilization (AOU) and N * -DINxs Indices
Dissolved oxygen (DO) was determined using the modified Winkler method according to Strickland and Parsons (1972) with an accuracy of ±1.3 µmol L −1 . The apparent oxygen utilization (AOU) represents one estimate of the O 2 utilized by biochemical processes relative to a preset value. AOU (mL L −1 ) is calculated as the difference between the O 2 gas solubility (O * 2 ) and the measured O 2 concentration and is expressed as follows: where O * 2 is calculated as a function of in situ temperature and salinity at one atmosphere of total pressure. The O * 2 values were calculated using the equation of Garcia and Gordon (1992) based on the O * 2 values of Benson and Krause (1984). Additionally, O 2 is the measured O 2 concentration (mL L −1 ). Dissolved inorganic nutrients (ammonia (NH 3 +NH + 4 ) + nitrite (NO − 2 ) + nitrate (NO − 3 ), phosphate (PO − 4 ), and reactive silicate (SiO − 2 )) were analyzed according to Grasshoff et al. (1983) Gruber and Sarmiento, 1997;Deutsch et al., 2001;) and DINxs (DINxs = NO − 3 -16PO − 4 ; Bates and Hansell, 2004;Hansell et al., 2004) represent the relative abundances of nitrate and phosphate. These indices measure the departure from classical Redfield ratios of the dissolved inorganic forms of nitrogen and phosphorus (the N * and DINxs indices differ only in the offset of 2.90 µmol L −1 , a value that was intended to fix the global mean N * value to zero). Negative values of DINxs (or N * values < 2.9 µmol L −1 ) indicate a deficit in N relative to P with respect to the requirements for Redfieldian production of organic matter; positive values of DINxs (or N * values >2.9 µmol L −1 ) indicate excess N relative to P. The values obtained for N * are shown in Table 2.

The CO 2 System
Seawater samples were collected for total inorganic carbon (DIC) and total alkalinity (TA) analyses to assess the key parameters of the CO 2 system. DIC and TA were measured via potentiometric titration using a closed cell, following the method of Edmond (1970). Equivalent points were calculated using the code published by (The Department of Energy, 1994). Certified reference material, supplied by Professor A. Dickson (Scripps Institutions of Oceanography, San Diego, USA), was used for calibration. The accuracy was estimated at 3 µmol kg −1 .
The sea surface fCO 2 was calculated from TA, DIC, temperature and salinity using the CO2calc R software (Robbins et al., 2010). The K 1 and K 2 dissociation constants of the carbonic acid used in the calculations were those obtained by Mehrbach et al. (1973) and refitted by Dickson and Millero (1987), and the sulfate dissociation constants were from Dickson (1990a,b).
The monthly averaged atmospheric CO 2 mole fraction (XCO 2 atm, ppm) recorded at the NOAA/Earth System Research Laboratory (ESRL) Global Monitoring Division station closest to the Amazon River plume (Ragged Point, Barbados, 13.17 • N, 59.43 • W; http://www.esrl.noaa.gov/gmd/ccgg/iadv/) was used for the atmospheric fCO 2 (fCO 2 atm) and sea-air CO 2 flux calculations. The fCO 2 atm was calculated as follows: where P is the atmospheric pressure (atm), pH 2 O is the water vapor pressure at 100% humidity (atm) calculated from SST and SSS, and Cf is the fugacity coefficient calculated according to Weiss (1974). The atmospheric pressure in October 2012 was obtained from the National Centers for Environmental Prediction (NCEP)/National Center for Atmospheric Research Reanalysis project, initially at a 2.5 • resolution, and linearly interpolated to the scale of the working grid (1/4 • resolution). The air-sea CO 2 flux (CO 2 fluxes, in mmol m −2 d −1 ) was then calculated as follows: where S o is the solubility of CO 2 (mol kg −1 atm −1 ) as a function of SST and SSS (Weiss, 1974), k is the gas transfer velocity (m d −1 ), and fCO 2 sw is the fCO 2 of the surface ocean waters (calculated from that measured along the CF3 tracks). The term k was calculated according to Sweeney et al. (2007): where Sc is the Schmidt number and U 10 is the wind speed (m d −1 ) at 10 m above sea level. U 10 data for October 2012 were taken from Advanced Scatterometer (ASCAT) observations at 10 m (ftp://ftp.ifremer.fr/ifremer/cersat/products/gridded/ MWF/L3/ASCAT/, 1/4 • resolution). A positive flux value represents net release of CO 2 from the sea surface, and a negative flux value represents absorption of atmospheric CO 2 by the sea.

Plankton Communities and Chl-a
Samples for the study of phytoplankton communities were taken using vertical tows with a plankton net with a 20 µm mesh, starting at 10 m below the recorded deep chlorophyll maximum (DCM) up to the surface. In total, 18 tows were performed. The samples were subsequently fixed with neutral formaldehyde solution to a final concentration of 4%, according to the methods of Newell and Newell (1963). To analyze the phytoplankton composition, the samples were gently mixed to homogeneously suspend the organisms, then 0.5 mL aliquots were qualitatively analyzed using an optical microscope with 100× and 400× magnification. The organisms were identified by consulting specialized literature. The international database Algaebase was used to classify and check the scientific names of the taxa (Guiry and Guiry, 2016). The relative abundances of the taxa were calculated as described by Lobo and Leighton (1986). The specific diversity of the phytoplankton was evaluated using the Shannon index (Shannon, 1948), and the evenness was calculated according to Pielou (1977) using PRIMER 6.0.
The method for determining the Chl-a concentration was the spectrophotometric analysis described in UNESCO (1966).
For the zooplankton sampling, a Bongo frame with four nets (with mesh sizes of 64-, 120-, 300-, and 500 µm) was used. For this study, samples from three nets with the following mesh sizes and diameters were analyzed: 64 µm/30 cm; 120 µm/30 cm, and 300 µm/60 cm. The Bongo net was hauled obliquely at a speed of 2-2.5 knots at depths between 15 m nearshore and from 200 m to the surface at stations beyond the shelf break. A flowmeter (Hydrobios, Kiel) was attached to the opening of each net. Additionally, plankton nets with a 200 µm mesh size were used for vertical hauls. At each station, two vertical hauls were conducted: one shallow vertical haul from the base of the mixed layer to the surface (Tropical Surface Water, TSW) and one deep vertical haul from the mid-thermocline to the surface (TSW and South Atlantic Central Water, TSW+SACW). In areas shallower than 70 m with no thermocline, only a single haul was conducted, from 5 m above the bottom to the surface. Neuston tows were also conducted at all stations. The neuston net used was a David-Hempel nautical aluminum catamaran manufactured by Hydro-Bios (Kiel, Germany). This equipment is composed of two superimposed nets, one at the air/water interface and another 7.5 cm below the interface. The two nets have rectangular mouths, each with a width of 30 cm and a height of 15 cm. The upper net samples 7.5 cm above the airwater interface to 7.5 cm below the interface. The lower net, equipped with a flowmeter (Hydro-Bios, Germany), samples the sub-surface layer (or hyponeuston), from 7.5 cm depth to 22.5 cm depth. The total mouth area is 0.066 m 2 (0.022 m 2 for the upper net and 0.044 m 2 for the lower net). The sampling duration was approximately 20 min. A total of 291 samples were taken and preserved in a 4% buffered formalin-seawater solution, buffered with 0.5 g L −1 sodium tetraborate. The biomass was estimated via the wet-weight method (Omori and Ikeda, 1984).

Statistical Analyses
To test for differences between coastal and oceanic stations and between day and night, we used Mann-Whitney Utests or Student's t-tests, depending on the normality and homoscedasticity of the data (Zar, 1996). Accordingly, Kruskal-Wallis ANOVA, linear regression analysis, and Pearson's correlation (Zar, 1996) were used to analyze the relationships between variables.
Cluster analysis was used to identify spatial divisions within the cruise track. Physical (σ-t), chemical (NO − 3 and DIC) and plankton (phytoplankton and zooplankton biomass) parameters were used for Hierarchical Agglomerative Cluster (CAH) analysis of Pearson's similarity and were agglomerated via the unweighted pair-group average method.
Principal components analysis (PCA) was performed to analyze the structure of the multivariate data using 27 key oceanographic parameters.
All statistical analyses were performed using the XLSTAT R 2010 software.

Temperature, Salinity, and Density
Overall, a mean SST of 28.7 ± 0.9 • C was measured along the cruise track between 50.9 • and 38 • W (Figure 2A). The SST values varied during the cruise, reaching a maximum difference of 3.9 • C and exhibiting a positive linear trend of 0.087 • C (y = 0.087(SST) + 27.62) along the cruise track. However, no significant differences were found between day and night (t-test; p: 0.65; α: 0.05). The large temperature range was caused by two extreme values at stations 8 (26.59 • C) and 15 (30.48 • C) (Figure 2A). The removal of these SST values from the observations does not significantly alter the overall mean value (average without extreme values: 28.7 ± 0.7 • C). In addition, the coefficient of variation (CV) was 0.03%.
The SSS showed typical brackish water values during most of the cruise (>50% of the samples), as reflected in the mean value of only 34.1 ± 3.7 practical salinity units (psu). The day/night cycle did not affect the salinity values (t-test; p: 0.24; α: 0.05). The salinity data showed a higher CV than temperature (0.11%) ( Figure 2B).
The density values (σ-t) (average: 21.5 ± 2.8) were low along the track. For example, the SST average for the track (28.7 • C) and SSS = 35 correspond to an σ-t of 22.16, whereas the highest SST value (30.5 • C) and SSS = 35 correspond to an σ-t of 21.56. The observations during this period showed low density values typical of brackish waters in the initial and final parts of the transect ( Figure 2C). The track showed a minimum σ-t value of 10.86 at station number 1. This station is located near the coast at 0.46 • S and 48.25 • W (Figure 1).
The nitrate concentrations were higher than the ammonia and nitrite concentrations during the entire track. The highest concentrations were in stations 3 and 8 (3.85 and 2.90 µmol L −1 , respectively). The average concentration was 1.02 ± 0.9 µmol L −1 with a CV of 0.94%, and no significant differences were observed between day and night (t-test; p: 0.48; α: 0.05) ( Figure 3C).
The sum of the concentrations of dissolved inorganic nitrogen compounds (ammonia + nitrite + nitrate = DIN) was calculated to obtain a relationship with phosphate concentrations.
The DIN average obtained for this period was 1.25 ± 1.2 µmol L −1 . The contributions of these compounds to DIN showed that the nitrate represented the greatest contribution, at 82.2%, followed by ammonia at 15.7% and nitrite at 2.1%.
The phosphate (PO − 4 ) concentrations featured an average of 0.1 ± 0.06 µmol L −1 for the study period ( Figure 3D). No significant differences were observed between day and night (ttest; p: 0.52; α: 0.05). The highest value (0.3 µmol L −1 ) was recorded at station 8 of the transect. The average ratio between DIN and PO − 4 was 15:1. Large variations in this ratio were observed during the study, mainly due to the low values of some compounds. The calculated CV was 1.6%, and the standard deviation was 40.4. The silicate (SiO − 2 ) showed three considerable peaks along the cruise track ( Figure 3E). The mean concentration obtained was 6.4 ± 6.4 µmol L −1 , and the maximum observed value was 27.0 µmol L −1 at station 16. No significant differences were observed between day and night (t-test; p: 0.68; The mean ratio of SiO

DO, Saturation Rate and AOU
The DO concentrations were always >3.8 mL L −1 , and the mean value was 4.7 ± 0.2 mL L −1 (Figure 4A). The CV was 0.06% for this study. The highest concentration was observed at station 9 (5.5 mL L −1 ). During the day/night cycle, no significant differences were observed (t-test; p: 0.41; α: 0.05). Based on the percent saturation of DO results, 91% of the samples exhibited supersaturation values (>100%), whereas only two samples (stations 1 and 8) exhibited undersaturation values ( Figure 4B). The mean value of the percent saturation was 107.0 ± 6.4%. Similar to the percent saturation, the AOU values at stations 1 and 8 differed from those of other stations. Positive values were observed at stations 1 and 8 (+0.2 and +0.7 mL L −1 , respectively), while the remaining 91% of the samples showed negative values. The mean value was −0.3 ± 0.2 mL L −1 , and the lowest values were observed at stations 9-17 (middle section of the track) ( Figure 4C).
The fCO 2 sw values calculated from the measured DIC and TA values varied between 304.1 and 543.5 µatm during the study period, with a mean value of 407.8 ± 47.8 µatm ( Figure 6A). In October 2012, the average atmospheric fCO 2 was 378.5±0.8 µatm. Regions of low saturation were observed at station 1 and where the SSS values were <35 in the final portion of the track ( Figure 6B).
The fCO 2 sw values did not show a strong correlation with SSS (r 2 = 0.2), especially for SSS values <35. We observed three stations (1, 8, and 9) with values that plot off the trend line. Without these values, the correlation between fCO 2 and SSS rises to r 2 = 0.8.
The calculated CO 2 fluxes varied between −8.6 and +8.4 mmol m −2 d −1 , while the mean was +1.6 ± 3.4 mmol m −2 d −1 , and the CV = 2.1% ( Figure 6C). The variations in the CO 2 fluxes recorded for this period showed that 75% of the samples acted as CO 2 sources to the atmosphere and that 25% acted as CO 2 sinks. Most positive CO 2 fluxes were associated with temperatures in the range of 27.6-28.6 • C.

Plankton Communities and Chl-a
Phytoplankton communities displayed a strong gradient in terms of total individuals. The highest values were observed at stations 9 and 10 (>3,800 individuals), and relatively high values were also observed at stations 8, 15, and 16 (1,000 < individuals < 3,800). The mean value for total phytoplankton was 1,082 ± 2,100 individuals ( Figure 7A).
The zooplankton was composed of 10 phyla: Protozoa, Cnidaria, Mollusca, Annelida, Crustacea, Bryozoa, Brachiopoda, Chaetognatha, Echinodermata, and Chordata. A total of 178 taxa were identified based on the lowest distinguishable taxonomic ranking. Holoplankton dominated in the study area, representing nearly 85% of the zooplanktonic biomass (It means 85% of total identified organisms-biodiversity). The most diverse and abundant group was Copepoda, with 130 species and accounting for more than 60% of the zooplankton. Two zooplankton communities were identified in the area: a low-diversity, high-abundance coastal community present at inshore stations and a high-diversity, low-abundance oceanic community at offshore stations. A maximum abundance zone occurred around the shelf break. Copepoda played a central role in the marine food web of the area. The Copepoda species that dominated offshore areas were Undinula vulgaris, Euchaeta marina, Nannocalanus minor, Clausocalanus furcatus, Scolecithris danae, Calocalanus pavo, Corycaeus (Corycaeus) speciosus, Farranula gracilis, and Oithona plumifera. The Copepod species that were abundant in neritic waters and indicative of the NBC were Subeucalanus pileatus, Rhincalanus cornutus, and Temora stylifera. These species of copepods are known by the literature (Boltovskoy, 1981(Boltovskoy, , 1999 to occur is this water mass. Meroplankton, mainly zoea-stage larvae of brachyuran crabs, were abundant at coastal stations under stronger plume influence. Zooplankton biomass (64, 120, and 300 µm mesh nets) showed a huge range, from 7.6 to 1,605 mg m −3 . Figure 7C shows the spatial distribution of the three compartment types (net size). The lowest values of biomass were associated with the 300 µm mesh (average: 130 ± 174 mg m −3 ; green circles in Figure 7C), while the highest values were associated with the 64 µm mesh (average: 270 ± 265 mg m −3 ; black circles in Figure 7C). The biomass values showed significant differences among the nets (64, 130, and 300 µm, Kruskal-Wallis ANOVA; p: 0.0001), mainly between the 64 and 300 µm nets (Dunn test; p: 0.0001).
Vertical hauls were conducted with plankton nets with a 200 µm mesh size to sample two depth strata (deep hauls and shallow hauls through the upper mixed layer). The vertical tows with the 200 µm mesh net showed high biomass values (>1,000 mg m −3 ) for the upper mixed layer, while the deep hauls yielded much lower values, always less than 1,000 mg m −3 (Figure 7D). The mean values were 1,003 ± 946 and 271 ± 269 mg m −3 for shallow and deep hauls, respectively. The highest deep haul values were observed between stations 10 and 17. Samples from shallow and deep hauls showed significant differences (t-test; p: 0.0006) for the study period. Neuston tows were also conducted at all stations. Two superimposed nets were used, one at the air-water interface (upper neuston) and the other sampling from 7.5 to 22.5 cm below the interface (lower neuston). The highest biomass value was observed at station 13 for the lower net (1,734 ± 371 mg m −3 ; Figure 7E). However, the highest biomass value in the upper net (average: 500 ± 467 mg m −3 ) was associated with a lower net biomass value of 293 ± 371 mg m −3 . The highest combined values for the upper and lower nets (>1,000 mg m −3 ) were observed at stations 8, 12, 13, 14, and 17 ( Figure 7E). The two sets of samples (lower and upper nets) did not show significant differences (t-test; p: 0.09) for the studied period.

Cluster Analysis and PCA
The cluster dendrogram shown in Figure 8A represents the results of the algorithm forming groups of observations and then subgroups of observations. The algorithm successfully grouped all the observations. The dotted line represents the automatic truncation, leading to four groups ( Figure 8A). Groups 1 and 4 (displayed in red and green, respectively) are less homogeneous than the others. Group 2 (displayed in blue) is more homogeneous than groups 1 and 4 (minor variance). Stations 1, 9 and 10 in groups 1 and 4 (red and green colors) were associated with high values of phytoplanktonic and zooplanktonic biomass, low values of density (σ-t), low values of DIC and high values of NO − 3 , whereas group 2 showed the opposite characteristics of groups 1 and 4. Group 2 features the largest number of stations (17). Group 3 contains 4 stations, whereas group 4 features only one station (station 10).
The PCA identified two leading modes that account for 63% of the variability. The first mode (37%) opposes HCO − 3 , DIC, σ-t, SSS, TA, CO 2− 3 , CO 2 fluxes and fCO 2 sw to Mesozooplankton (Mesozoo), Bacillariophyta (Bac), Chla (Chl-a), Dinophyta and Microzooplankton (Microzoo), as shown in the bi-plot of the first two factors (Figure 8B). The principal components (PC) are plotted as a function of the twenty-four stations measured in October 2012 (the green, brown and cyan colors correspond to the third, fourth and fifth components, respectively, in Figures 8B and Table 1). The factor loadings of the five modes are given in Table 1. The first principal component (PC1) highlights the crossshelf variability, with strong differences between stations 1 and 9 and the other stations (blue and red in Figure 8B). The second mode explains 26% of the variance and is dominated by nutrients (NH + 4 , NO − 2 , NO − 3 , DIN, and PO − 4 ), AOU, CO 2 aq and fCO 2 atm (orange in Figure 8B) opposed to SST (+ axis) ( Figure 8B). PC3 (10%) is characterized by SiO − 2 and Cyanobacteria (Cya) (green in Figure 8B). PC4 (7.8%) is characterized by DO, which does not appear to be correlated with any other parameter on this axis (brown in Figure 8B).  Table 1 for PCA analysis.
Macrozooplankton appears in the fifth mode and shows no correlation with any other parameters of this analysis (cyan in Figure 8B). The shading corresponds to factor 1 (positive-blue and negative-red). The following values correspond to factor 2 (positive-purple and negative-yellow), factor 3 (positivegreen), factor 4 (positive-brown) and factor 5 (positive-cyan). The values in bold indicate a significant correlation.

DISCUSSION
The results obtained from this cruise improved our understanding of the processes and organisms responsible for carbon and nutrient cycling along this large-scale river-ocean continuum. The contributions of the different factors explaining the biogeochemical pathways and carbon fluxes in the oceanic area under the influence of the Amazon River plume are extremely complex. Thus, we focused our discussion on systematically linking the variability in the physical and biological processes to the biogeochemistry and carbonate system in this huge area.

Surface Circulation and Distribution of SST and SSS
The hydrographic results revealed very different conditions in the study area, ranging from shallow coastal conditions to offshore areas. The distribution of the sampling stations up to 8 • N-38 • W (station 24) showed that a large number of these stations were located within the complex system of currents in the region, but most were affected by two main currents, the NBC and the NECC. The NECC intensifies during summer (June, July, and August) and fall (September, October, and November), while the NBC retroflects in this period. The NBC retroflection (Wilson et al., 2002), a result of the vorticity balance, is shown in Figure 1B as the clockwise gyre centered at 5 • N, 45 • W. The retroflection is weaker in spring (March, April, and May) and stronger in fall (September, October, and November) (Urbano et al., 2008;Lefevre et al., 2014). We observed the retroflection in fall (October) (Figure 1B). This zonal current system plays an important role in modulating the heat flux in the tropical Atlantic. During this period, SSTs varied considerably along the transect, but this variation did not influence the overall mean value of 28.7 • C. Warmer waters (>28 • C) from the Amazon and the Pará rivers were widely observed along the inner continental shelf, the Amazon and Pará river mouths and offshore areas ( Figure 1C). The values observed here are similar to the historic SSTs in the region 0 • -9 • N, 36 • -52 • W for the period 1958-2011 (28.0 ± 0.3 • C) (Yu et al., 2008) and are in accordance with the findings of other authors, such as Lefevre et al. (2014) and Urbano et al. (2008).
The SST values showed an inverse relationship with nutrients, AOU, and dissolved CO 2 (CO 2 aq) ( Figure 8B) but were not correlated with SSS values, as indicated by Ibánhez et al. (2015) for the plume of the Amazon River. The differences between SST and SSS can be explained by the strong amplitude of SSS (higher than 8 psu) whereas, SSTs exhibited differences of less than 2 • C. Thus, over half of the 24 samples of the CF3 cruise (54%) used here showed a minimum SSS lower than 35 psu ( Figure 2B).
According to Ibánhez et al. (2016), a high correlation exists between rainfall and SSS in the plume region, and this relationship can affect >16% of the Amazon River plume area. The spatial extent of the influence of brackish waters was highly variable. Low SSS values were found at 8 • N, 38 • W, indicating the influence of the local surface water circulation on the spread of Amazon River waters in the WTNA (Figures 1F, 2B). The rainfall in October 2012 in the vicinity of the CF3 transect was approximately 0.5 mm h −1 (Figure 1E).

SSS and Carbonate System Parameters
The SSS was also highly positively correlated with the carbonate system parameters (HCO − 3 , DIC, TA, CO 2− 3 , CO 2 fluxes and fCO 2 sw). According to Koffi et al. (2010); Lefèvre et al. (2010), and Bonou et al. (2016) in the WTNA, SSS has been shown to be highly correlated with TA and DIC.
Both relationships (TA-SSS and DIC-SSS) are quite robust, as shown by the good agreement between the regression and the available observations from previous cruises (TA = 56.46 × SSS + 322.48; r 2 = 0.99 and DIC = 44.74 × SSS + 409.42; r 2 = 0.98). The TA-SSS and DIC-SSS regressions were performed for salinities greater than 19. Our TA-SSS data are in agreement with the relationship of Lefèvre et al. (2010). Ternon et al. (2000) reported TA-SSS and DIC-SSS slopes of 58.85 and 49.48 µmol kg −1 , respectively, using cruise data at the Amazon River mouth. These TA-SSS and DIC-SSS slopes are very close to those obtained in our study (56.4 and 44.7 µmol kg −1 , respectively). Salinity and fCO 2 sw values lower than those in the literature were observed during the CF3 cruise and were due to the influence of Amazon River water and rainfall. The maximum and minimum fCO 2 values were associated with SSS values of <35, both of which were observed in areas influenced by the Amazon plume. The fCO 2 sw values within the Amazon River plume are significantly correlated with SSS Ibánhez et al., 2016). However, Ibánhez et al. (2015) noted that, brackish waters (SSS<35) transported by the NBC and collected south of 8 • N showed the highest SSS-fCO 2 sw discrepancy among consecutive cruises between 2006 and 2013 (the period covering the CF3 cruise). According to these authors, many surface eddies present in this region of the NBC retroflection (a significant path of water mass transport in the area) (Ffield, 2005) may be responsible for the spatial variability found near the coastal zone.
The SSS-HCO − 3 and SSS-CO 2− 3 correlations were positive (r = 0.97 and r = 0.88, respectively) and were linked to the concentrations of TA and DIC. They represented 88% and 11.5% of the DIC, respectively, while the remaining 0.5% corresponded to CO 2 aq, which showed a low correlation with SSS. According to Richey et al. (1990) and Cooley and Yager (2006), the Amazon River mainstream DIC load represents 82% of the HCO − 3 and 18% of the CO 2 aq in the river. CO 2 aq is mainly controlled by pH, SST and SSS and the concentrations decrease along the salt gradient.
The CO 2 fluxes-SSS correlation was positive (r = 0.7) along the ship track. The large variation (17 mmol) was forced by stations 1 and 2, where the magnitude of the CO 2 fluxes was also influenced by the intensity of the winds, with speeds of 7.3 and 9.0 m s −1 , respectively. The majority of positive CO 2 fluxes were associated with temperatures of 27.6-28.6 • C and salinities of 34-37 psu. Our study showed that 75% of the 24 stations exhibited CO 2 oversaturation and that only 25% of the sites exhibited undersaturation. These results may be associated with the route followed by the CF3 cruise. Figure 1F shows that a portion of the stations are located within a region with SSS values of >35 (St. 2-7). Another factor that can affect the direction of the CO 2 fluxes is light. The amount of light affect the growth of phytoplankton in the oceans, and their rate of photosynthesis increases in proportion to the light intensity. During the CF3 cruise, 33.3% of the samples were obtained at night. During this period, productivity is less than respiration, and heterotrophic processes in the water column can increase the CO 2 aq and decrease the DO concentrations. Beyond the retroflection region, we observed a greater fluctuation around zero (oversaturation/undersaturation) in the CO 2 fluxes.
According to Ibánhez et al. (2015), CO 2 undersaturation in the brackish waters transported by the NBC through retroflexion to the north occurs during the second half of the year. In this region and during the boreal autumn, the Amazon discharge and the location of the ITCZ can lead to undersaturation . The migration of the ITCZ is consistent with the salinity distribution with low salinities located north of 2 • N from July to December. The amount of rainfall observed in this study was 10 times lower than that observed by Lefèvre et al. (2010) (Figure 1E).
According to Lefèvre et al. (2010), low salinities (SSS < 35) are encountered from July to December between 2 • and 8 • N. Similar values were recorded by the CF3 cruise during the month of October. Additionally, the low-salinity region is associated with a decrease in fCO 2 . According to Lefèvre et al. (2010), the eastward propagation of Amazon waters is unlikely to explain the CO 2 undersaturation observed throughout the year because Amazon waters might reach 25 • -30 • W only in boreal autumn, depending on the strength of the NECC.
The ∆fCO 2 values varied strongly during the CF3 cruise, passing from highly undersaturated conditions in the coastal region (1 • S) to oversaturated conditions between 0 • and 7 • N. During the final part of the ship track (45 • -38 • W), the fCO 2 varied slightly between positive and negative values (8 • N). Within this region, the CO 2 fluxes varied slightly between source and sink, whereas the fCO 2 sw values were very close to atmospheric fCO 2 (<5 µatm for fCO 2 ). Along the CF3 ship track in October, the calculated CO 2 fluxes (+1.6 ± 3.4 mmol m −2 d −1 ) had the same magnitude as the values reported by Lefèvre et al. (2010) and Ibánhez et al. (2015) for the same period.

Nutrients
Riverine nutrients and in situ organic matter mineralization support primary production in the offshore plume (DeMaster and Pope, 1996;Subramaniam et al., 2008;Yeung et al., 2012). The effects of freshwater inputs on coastal CO 2 and the carbonate system dynamics occur via direct inputs of DIC or through enhanced primary production due to river-borne nutrient inputs (Kitidis et al., 2012).
To further examine the effects of nutrient concentrations on the carbonate system of coastal waters, N * and DINxs indices were applied to determine the relative abundances of DIN and PO − 4 in the study region. The nitrate concentrations were higher than the ammonia and nitrite concentrations throughout the cruise.
Nitrate represented 82.2% of the DIN, ammonia represented 15.7% of the DIN, and nitrite represented only 2.1% of the DIN. Thus, the average DIN:phosphate ratio was 15:1, although this ratio varied during the study due to the low values of some compounds.
The DIN compounds showed relative abundances lower than the limits of the indices DINxs and N * (average: −0.2 and 2.6, respectively) ( Table 2). Negative DINxs values (or N * values < 2.9 µmol L −1 ) indicate a deficit of N relative to P with respect to the requirements for Redfieldian production of organic matter. The values calculated are consistent with those reported by Hansell and Follows (2008) for the tropical Atlantic. Furthermore, the Redfieldian ratios suggest that the Amazon is an important source of "excess" PO 4 − (DIN:PO 4 − < 16) to the WTNA. The compounds are transported thousands of kilometers offshore via the plume. The plume ranges from 5 to 25 m thick (Coles et al., 2013) and supplies allochthonous Si and PO − 4 to offshore regions in the tropical Atlantic. This input of plume-derived Si and PO − 4 into the nitrogen-limited open ocean, with Si:DIN and PO − 4 :DIN ratios in excess of those typically needed by phytoplankton, creates a distinct niche for N 2 fixation by diatom-diazotroph associations (DDAs), leading to enhanced primary production in this region (Subramaniam et al., 2008). A key hypothesis advanced by Subramaniam et al. (2008) is that DDAs represent an effective biological pump in tropical river plumes. In coastal regions, DIN derived from inland waters is rapidly consumed, leaving an extensive area (∼10 6 km 2 ) with lower salinity and excess dissolved SiO − 2 and PO − 4 concentrations relative to Redfield-Brzezinski stoichiometry (i.e., C:Si:N:P = 106:15:16:1) in the tropical Atlantic (Brzezinski, 1985;Shipe et al., 2007;Subramaniam et al., 2008). In the study region, the ratio SiO − 2 :DIN:PO − 4 showed a Redfield-Brzezinsky stoichiometry of 15:80.5:1. The Redfield stoichiometry was lower than 16, indicating that nitrogen was the limiting factor in this region. Plots of nutrients vs. salinity are shown in  Edmond et al. (1981) and Ternon et al. (2000).

Relations between SiO − 2 , DIC, and TA Concentrations
We used the conservative theoretical mixing lines and the mixing line of observed data for the SiO − 2 , TA and DIC concentrations to obtain predictions of the effects of biological activity in the study region.
Based on a comparison between the regression line fitted to observed data and the conservative mixing line, the maximum depression in SiO − 2 ( Figure 1SC; Supplementary Material) due to biological activity was approximately 43 µmol kg −1 (19-62 µmol kg −1 , respectively), which is in agreement with the highest SiO − 2 depletion (greater than 30 µmol kg −1 ) reported by DeMaster and Pope (1996) on the Amazon shelf.
To compare the SiO − 2 and inorganic carbon depletions involved in biological uptake (maximum depressions of 43 and 82 µmol kg −1 for SiO − 2 and DIC, respectively), we need to correct the inorganic carbon depletion for calcification ( Figures 1SC, 2S; Supplementary Material). This correction is always DIC-TA/2, according to Zeebe and Wolf-Gladrow (2001).
The observed TA deviation away from the conservative line (71 µmol kg −1 for SSS = 19.7, Figure 2S, Supplementary Material) is typical of calcium carbonate (CaCO 3 ) production is seawater (Broecker and Peng, 1982). As the change in TA due to carbonate mineral production is twice as large as the change in DIC (Skirrow, 1975) and by neglecting the small changes in TA due to the production and decay of organic matter (Brewer and Goldman, 1976), the organic carbon production corrected for calcification is 82-71/2 = 46 µmol kg −1 . Thus, SiO − 2 and inorganic uptake exhibit a molar ratio of SiO − 2 /DIC = 0.9, which is beyond the range (0.15-0.4) of molar silicon/carbon production ratios reported by  and Ternon et al. (2000) in waters of the Amazon shelf. This high value (0.9), in comparison to the mean value of 0.13 for low-latitude diatoms (Brzezinski, 1985), highlights the dominant role of diatoms (mainly Fragillaria sp. and Pseudo-nitzschia pungens) in the primary production in this zone (stations 8, 9, and 10; Figure  1SC and Figure 7A).

DO and AOU
The decomposition of organic matter changes the concentrations of carbon, nitrogen, phosphorus, oxygen and TA in the ratio 106:16:1:138:-17. We observed areas of DO-supersaturated surface water, indicative of low consumption. DO concentrations were always >3.8 mL L −1 , and the mean value was 4.7 ± 0.2 mL L −1 . AOU concentrations along the CF3 ship track ranged between −0.9 and +0.7 mL L −1 . However, 91% of the samples showed negative values. The mean value was −0.3 ± 0.2 mL L −1 , and the lowest values were observed at stations 9-17 (middle section of the track). Thus, the sampling period was characterized by negative AOU values (production > respiration). The oxygen supersaturation in the river plume (negative AOU) is evidence of high photosynthetic activity. Other authors DeMaster and Aller, 2001;Garcia et al., 2006) also reported negative AOU values.
In addition, biological consumption is one of the processes affecting the variability in the carbon parameters in tropical areas (Cooley et al., 2007;da Cunha and Buitenhuis, 2013;Araujo et al., 2014). We separate the stations with SSS<35 and SSS≥35 along the CF3 ship track. A negative average AOU value of −0.3 mL L −1 was found for both divisions of salinity. According to the PCA (Section Cluster Analysis and PCA above), AOU and CO 2 aq showed a positive correlation in factor 2 (orange color in Figure 8B). Thus, negative AOU values are associated with higher CO 2 aq concentrations. According to Zeebe and Wolf-Gladrow (2001), the release of CO 2 to the atmosphere decreases the DIC concentration, while the TA concentration remains constant. This process leads to a rise (drop) in dissolved CO 2 (CO 2 aq), with the opposite change in pH. The AOU values were characteristic of productive regions and indicated that production was greater than respiration at 91% of the CF3 cruise stations.

Plankton Community and Chl-a
The phylum Miozoa and Bacillariophyta characterized 94% of the floristic diversity in the planktonic flora. Dinoflagellates and diatoms were present at all stations of the ship track, whereas Cyanobacteria were more abundant at stations 15-24 ( Figure 7A). A dominance of diatoms was observed in the region from the river mouth to the beginning of the area affected by the NBC retroflection (Chl-a concentrations ranging from 0.02 to 0.9 mg m −3 ). After retroflection, the NECC region is fully oligotrophic, and the most representative groups are Cyanobacteria and Bacillariophyta.
A comparison of Cyanobacteria vs. dinoflagellates showed a slight correlation (Pearson correlation; ρ: 0.58), whereas diatoms showed significant negative correlations with σ-t (ρ: −0. Increases in phytoplankton accompanied decreases in σ-t and DIC. Thus, the variations in the phytoplankton community were reflected in the concentrations of the parameters of the carbonate system (DIC and TA). In addition, the stations located within the NBC (stations 6-10) system were associated with the Bacillariophyta group, while stations located within the NECC (stations 15-18) featured greater numbers of Cyanobacteria individuals. However, the diversity of groups was always greater after retroflection (Figures 1B, 7B). Furthermore, the highest productivity values (negative AOU) were associated with these phytoplankton groups.
Trichodesmium sp. and Richelia sp. were the main species of Cyanobacteria in the region of the NECC. According to Yeung et al. (2012), phytoplanktonic groups varied along the salinity gradient, and specific groups co-varied. For example, the abundance of Richelia sp. was associated with that of Hemiaulus hauckii. In the NECC region, we observed an association of Richelia sp. and Trichodesmium sp. with Hemiaulus hauckii and Rhizosolenia sp. The planktonic cyanobacteria Trichodesmium sp. is globally distributed in the tropical and subtropical oceans (Luo et al., 2012), where water temperatures are above 20 • C (Detoni et al., 2016). Understanding the global distribution of Trichodesmium is particularly important because of its ability to fix molecular nitrogen (N 2 ) (Yeung et al., 2012). The diatoms Hemiaulus hauckii and Rhizosolenia spp. containing the symbiotic Richelia sp. (DDAs) represented 21% of the total phytoplankton species at the mesohaline (32 ≤ SSS < 35) stations (15-24). According to Subramaniam et al. (2008), the composition of the phytoplankton community changes along the Amazon River plume from the mouth to the open ocean in response to changing nutrient availability. At low-salinity stations, sufficient PO − 4 , SiO − 2 , and N compounds are available at the surface to support coastal diatom species, and very little N 2 fixation occurs in these areas (stations 1-9). As the N compounds are assimilated and the plume is mixed with low-nutrient ocean waters, diazotrophs become significant sources of N. The diatom hosts of Richelia, the dominant diazotroph at the mesohaline stations, require the SiO − 2 and PO − 4 found in the river plume but N is supplied via N 2 fixation. Farther "downstream, " where riverassociated SiO − 2 and PO − 4 are depleted, the species composition transitions to that typical of oligotrophic tropical oceans, and the dominant diazotroph is Trichodesmium.
Thus, we conclude that the non-conservative changes in DIC discussed above were associated with the Cyanobacteria group and consequently with N 2 fixation.
The zooplankton was composed of the phyla Protozoa, Cnidaria, Mollusca, Annelida, Crustacea, Bryozoa, Brachiopoda, Chaetognatha, Echinodermata, and Chordata. In total, 178 taxa were identified, considering the lowest taxonomic unit possible for each phylum. This area was dominated by holoplankton, which represented nearly 85%. Copepoda was the most diverse and abundant group with 130 species, accounting for more than 60% of the zooplankton.
The highest biomass was registered at station 5 at a depth of 46 m followed by station 1 at a depth of 9.2 m, which are strongly affected by the river plume. Both stations exhibited blooms of the diatom Coscinodiscus centralis. Additionally, station 5 was dominated by a high density of medusa followed by Copepoda Undinula vulgaris and Lucicutia flavicornis (adults, copepodite and nauplii), and station 1 featured a high density of characteristic estuarine indicator species (Acartia tonsa, Paracalanus sp., Oithona hebes, Euterpina acutifrons) in addition to numerous Decapoda larvae.
The highest biomass density offshore was registered at station 17 under the plume influence and was dominated by Clausocalanus furcatus, Oithona plumifera, and Oncaea media. The densities of these species were associated with blooms of Trichodesmium sp. in the oceanic area. The neuston biomass densities were higher at the coastal stations 5 and 9 due to medusa blooms and offshore at stations 12, 13, 14, and 17 due to the presence of fish larvae and gelatinous organisms.
Two zooplankton communities were identified in the area: a low-diversity, generally higher-biomass and higher-density coastal community present at inshore stations and a highly diverse, generally low-density oceanic community at offshore stations. A few oceanic stations registered high biomass values due to jellyfish blooms. A maximum biomass/density zone occurs around the shelf break. Throughout the study area, Copepoda play a central role in the marine food web. Meroplankton individuals, mainly Brachyura zoeae, are abundant at coastal stations under plume influence.
According to the multivariate analysis (PCA), zooplankton biomass showed strong correlations with microzooplankton, mesozooplankton, Chl-a and Cyanobacteria, whereas macrozooplankton biomass was not correlated with other parameters (Factor 4).
Recently, Conroy (2016) provided direct evidence demonstrating that two DDAs, Hemiaulus-Richelia and Rhizosolenia-Richelia, are consumed by mesozooplankton. He further showed that calanoid and harpacticoid copepods, as well as some decapod larvae, consume Trichodesmium. Additionally, he showed that unicellular cyanobacteria, particularly nondiazotrophic Synechococcus and Prochlorococcus, as well as diazotrophic (unicellular nitrogen-fixing cyanobacteria, or UCYN-A), are consumed by zooplankton, likely as components of aggregates. Grazing on UCYN-A provides an additional and previously undocumented pathway for diazotrophic nitrogen incorporation into the food web.
Thus, we conclude that the changes in DIC and fCO 2 sw and CO 2 fluxes in the mesohaline stations were also associated with the Cyanobacteria group (phytoplankton) and N 2 fixation (mesozooplankton). Microzooplankton represent 25% of the total zooplanktonic biomass in this region, while mesozooplankton represent 21%. We have to consider that the biomass of station 1 near the mouth of the river is composed of 18% microplankton and 35.5% mesozooplankton.

Clustering Analysis and PCA
Clustering was included to identify spatial divisions within the CF3 ship track. We included the main parameters of each area: biological (phytoplankton and zooplankton biomass), physical (σ-t) and chemical (NO − 3 and DIC). The differences between the main groups were analyzed according to the similarity within the dendrogram. Groups 1 and 4 (red and green colors in Figure 8A, respectively) are associated with low values of σ-t, high values of phytoplankton biomass (diatoms group) and zooplankton biomass, high NO − 3 concentrations, and low DIC concentrations. Group 1 and 4 differ in that station 10 has a higher σ-t value than the other 2 stations. Groups 3 and 4 show a smooth similarity differing mainly because of σ-t. These groups show a mix of stations associated with the NBC and NECC region.
The PCA identifies three leading modes that account for 73% of the variability encountered. The first mode (37%) sets HCO − 3 , DIC, σ-t, SSS, TA, CO 2− 3 , CO 2 fluxes and fCO 2 sw (shown in blue in Figure 8B) in opposition to mesozooplankton, Bacillariophyta, Chl-a, Dinophyta, and microzooplankton (shown in red), as presented in the bi-plot of the first two factors (Figure 8B). In this analysis, we do not consider species.
We found a significant negative correlation between biological parameters and DIC (Table 1). Other parameters of the carbonate system, such as HCO − 3 , CO 2− 3 fCO 2 , and CO 2 fluxes, were also negatively correlated with biological groups. On the other hand, we posit that the composition of the phytoplankton community changes along the Amazon River plume from the mouth to the open ocean in response to changing nutrient availability. Nitrogen compounds and phosphates did not show associations in the first mode (shown in orange in Figure 8B). Stations located near the coast (stations 1 and 8) show a strong association with nutrients and AOU (orange in Figure 8B). Diatom groups that increase productivity (negative AOU) and release CO 2 aq dominate this region (orange in Figure 8B). Based on the cluster analysis, the region exhibits segmentation, as a cyanophyte bloom was observed in the NECC region (stations 15-18). We applied a new PCA using the most abundant phytoplankton and zooplankton groups and species in this region. Additionally, we include NO − 3 , DIC, SSS, and Chl-a. The results of the multivariate analysis showed a strong association among the cyanophytes Trichodesmium sp. and Richelia sp., the diatom group Rhizosolenia sp. and mesozooplankton (Mesozoo) in the first mode (40%) (shown in green in Figure 8C). These species showed a negative correlation with SSS and DIC (shown in blue in Figure 8C). In addition, NO − 3 , Chl-a and microzooplankton (Microzoo) did not show an association with these parameters (Figure 8C). The mesozooplankton is mainly composed of copepods (>60%) in this region.

CONCLUSION
The ship track during the oceanographic cruise Camadas Finas III encompassed the outer Amazon River estuary, the alongshore northwestern NBC region, the NBC retroflection area and the eastern NECC plume transport to 38 • W. The cruise was purposefully planned to take place during boreal autumn (October 2012), when the dispersal of Amazonian waters forms a brackish plume that can reach 25 • W when the NECC is strong.
Hydrographic results showed very different situations, ranging from shallow well-mixed coastal scenarios to offshore areas where low-salinity Amazonian waters induce the formation of barrier layers inhibiting vertical mixing of heat and nutrients. Ship track current measurements noted strong alongshore NBC flow and a meandering NECC, which produces large-scale anticyclonic rings that are transported eastward. Nutrients, mainly NO − 3 and SiO − 2 , were strongly depleted in coastal regions, and the autotrophy was greater than the heterotrophy (negative AOU). In terms of phytoplankton groups, diatoms dominated the region from the river mouth to the edge of the area affected by the NBC retroflection (Chl-a ranging from 0.02 to 0.94 mg m −3 ). Additionally, the NECC region is fully oligotrophic where the most representative groups are cyanobacteria and dinoflagellates (Chl-a ranging from 0.02 to 0.40 mg m −3 ). Copepods were the most diverse and abundant group of the zooplankton, playing a central role in the marine food web: 130 copepod species were identified, and they accounted for more than 60% of the zooplankton abundance. Two different zooplankton communities are represented in the area: a low-diversity, high-density coastal community present at inshore stations and a high-diversity, low-density oceanic community present at offshore stations. Copepods dominated offshore areas, whereas macrozooplankton (mainly Brachyura zoeae) dominated coastal stations under stronger plume influence. Based on the multivariate analysis, phytoplankton and zooplankton showed correlations with carbonate system parameters (DIC, TA, fCO 2 sw, HCO − 3 , CO 2− 3 , and CO 2 fluxes). The fCO 2 sw values reached 543 µatm in the coastal region but oscillated near the value of atmospheric fCO 2 (379 µatm) offshore. Lower fCO 2 sw values were observed in the NECC area. The fCO 2 in this region was less than 5 µatm (−0.3 mmol CO 2 m −2 d −1 ), while in the coastal region before retroflection, the fCO 2 value was approximately 50 µatm (+3.7 mmol CO 2 m −2 d −1 ). The fCO 2 values varied considerably along the CF3 cruise, passing from high undersaturation in the coastal region (0.5 • S) to oversaturation between 0 • and 7 • N. In the final portion of the ship track (45 • -38 • W), fCO 2 varied slightly between positive and negative values (8 • N). Additionally, in the NECC region, blooms of species in the cyanophyte group (Richelia sp. and Trichodesmium sp.) were associated with the diatom group (Rhizosolenia sp.) and mesozooplankton (Copepods).
This study provides foundational data for future processoriented high-resolution numerical modeling experiments, in which physical-biogeochemical mechanisms can be examined together as drivers of the observed geographical, seasonal and interannual variabilities in the AROC. These studies are currently underway.

AUTHOR CONTRIBUTIONS
MA conceived the idea and coordinated the onboard activities during the CF3 cruise. MF and KT performed the chemical analysis. SN, RS, and PM performed the zooplankton analysis.
AO and FF performed the phytoplankton analysis. MA, GH, and JA performed the physical analysis. CN, DV, NL, and LB performed the carbonate parameter analysis, CO 2 analysis, and statistical tests. All authors contributed extensively to the interpretation of the results and to writing the manuscript.