Original Research ARTICLE
Evaluation of Primary Production in the Lower Amazon River Based on a Dissolved Oxygen Stable Isotopic Mass Balance
- 1School of Oceanography, University of Washington, Seattle, WA, USA
- 2Marine Sciences Laboratory, Pacific Northwest National Laboratory, Sequim, WA, USA
- 3Whitney Laboratory for Marine Bioscience, University of Florida, St. Augustine, FL, USA
- 4Centro de Energia Nuclear na Agricultura, Universidade de São Paulo, Piracicaba, Brazil
- 5Departamento de Meio Ambiente e Desenvolvimento, Universidade Federal do Amapá, Macapá, Brazil
- 6Instituto Sócio Ambiental e dos Recursos Hídricos, Universidade Federal Rural da Amazonia, Belem, Brazil
- 7Departamento de Sensoriamento Remoto, Instituto Nacional de Pesquisas Espaciais, São José dos Campos, Brazil
The Amazon River outgasses nearly an equivalent amount of CO2 as the rainforest sequesters on an annual basis due to microbial decomposition of terrigenous and aquatic organic matter. Most research performed in the Amazon has been focused on unraveling the mechanisms driving CO2 production since the recognition of a persistent state of CO2 supersaturation. However, although the river system is clearly net heterotrophic, the interplay between primary production and respiration is an essential aspect to understanding the overall metabolism of the ecosystem and potential transfer of energy up trophic levels. For example, an efficient ecosystem is capable of both decomposing high amounts of organic matter at lower trophic levels, driving CO2 emissions, and accumulating energy/biomass in higher trophic levels, stimulating fisheries production. Early studies found minimal evidence for primary production in the Amazon River mainstem and it has since been assumed that photosynthesis is strongly limited by low light penetration attributed to the high sediment load. Here, we test this assumption by measuring the stable isotopic composition of O2 (δ18O-O2) and O2 saturation levels in the lower Amazon River from Óbidos to the river mouth and its major tributaries, the Xingu and Tapajós rivers, during high and low water periods. An oxygen mass balance model was developed to estimate the input of photosynthetic oxygen in the discrete reach from Óbidos to Almeirim, midway to the river mouth. Based on the oxygen mass balance we estimate that primary production occurred at a rate of 0.39 ± 0.24 g O m3 d−1 at high water and 1.02 ± 0.55 g O m3 d−1 at low water. This translates to 41 ± 24% of the rate of O2 drawdown via respiration during high water and 67 ± 33% during low water. These primary production rates are 2–7 times higher than past estimates for the Amazon River mainstem. It is possible that at high water much of this productivity signal is the result of legacy advection from floodplains, whereas limited floodplain connectivity during low water implies that most of this signal is the result of in situ primary production in the Amazon River mainstem.
A considerable body of work over the last several decades has demonstrated that the world's inland waters are generally supersaturated with carbon dioxide (CO2), resulting in a large flux of CO2 from rivers and lakes to the atmosphere (Cole et al., 2007; Battin et al., 2009). This state of supersaturation is driven by the balance between CO2 outputs such as primary production, degassing, and export and CO2 inputs from soils/sediments, the riparian zone, physical and biological decomposition of organic matter (OM) derived from both the terrestrial biosphere (Mayorga et al., 2005; Ward et al., 2013) and aquatic production (Kritzberg et al., 2004; del Giorgio and Pace, 2008), and aquatic plant respiration in floodplains (Abril et al., 2014). Rivers outgas an order of magnitude more CO2 than lakes on a global scale (Raymond et al., 2013). This is commonly attributed to higher rates of primary production relative to respiration in still waters compared to flowing rivers (Tranvik et al., 2009).
Although it is clear that rivers provide a positive flux of CO2 to the atmosphere in most cases, quantifying the relative amount of carbon/energy that is produced vs. consumed in an ecosystem is central to understanding its overall productivity and the potential movement of energy up trophic levels (Odum, 1971). For example, rapid rates of carbon remineralization do not preclude productive fisheries. Lower and higher trophic levels depend on both autochthonous and terrestrial carbon/energy sources to thrive, and the balance in food sources to freshwater ecosystems is tightly linked to hydrology (Hoffman et al., 2008). Likewise, recent evidence has shown that the presence of fresh algal biomass and exudates can stimulate the breakdown of less reactive OM such as terrestrially-derived molecules, implying that systems receiving large inputs of both allochthonous and autochthonous OM are likely highly efficient at recycling carbon (Guenet et al., 2014; Bianchi et al., 2015; Ward et al., 2016).
In large turbid river systems such as the Amazon River it has long been assumed that primary production is minimal relative to heterotrophic respiration due to high sediment loads and low light penetration relative to river depth (Fisher, 1979; Devol et al., 1987; Richey et al., 1990; Hedges et al., 2000). In fact, early measurements of primary production based on 14C uptake indicated rates of primary production that are about an order of magnitude lower than respiration in the Amazon River mainstem, while clearwater tributaries in the basin's lowlands exhibited significantly higher rates of primary production (Wissmar et al., 1981; Benner et al., 1995). However, several studies have more recently found evidence for a higher relative contribution of primary production (and subsequent breakdown of algal OM) to ecosystem metabolism in the Amazon River based on the stable isotopic signature of dissolved oxygen (δ18O-O2) and the presence of algal biomarkers in the mainstem (Quay et al., 1995; Mortillaro et al., 2011; Ellis et al., 2012). Similar observations of high rates of primary production and algal abundance have also been made in the mainstem of the Congo River, a similar large tropical system (Descy et al., 2017).
These emerging results suggest that there is a large gap between the conventional understanding of production in turbid tropical rivers, and what is actually occurring in the environment. Gross primary production and respiration rates are critical components in understanding the production, transformation and export of OM and other biogeochemically important nutrients in aquatic ecosystems. The balance between these two is a key measure of the net metabolic state of the ecosystem. As rivers worldwide are saturated with respect to atmospheric CO2 and outgas this CO2 at a massive scale, understanding the exact dynamics of photosynthesis and respiration remain key to understanding the role of river systems in global biogeochemical cycles and have further implications on ecosystem productivity as previously described (Cole and Caraco, 2001; Butman and Raymond, 2011).
The Amazon River is the largest source of freshwater to the world's oceans and outgases 0.5 Pg C year−1, which is roughly 25% of recent global estimates (Richey et al., 2002; Raymond et al., 2013). The Amazon River ecosystem also supports the world's most diverse fishery with over 2500 unique species and is among the primary sources of protein for the region (Junk et al., 2007). The massive contribution of the Amazon basin to global biogeochemical cycling and regional economics makes it an ideal setting for evaluating large-scale processes. The lower Amazon River, from Óbidos to the river mouth, is an understudied component of the Amazon basin, but is a region where the interplay between primary production and heterotrophic respiration likely rapidly evolves due to a gradient in environmental conditions (Ward et al., 2016). Most historic measurements of the Amazon River have only been made as far downstream as Óbidos, roughly 800 km from the mouth of the river (Hedges and Clark, 1986; Quay et al., 1992; Moreira-Turcq et al., 2013). The lower Amazon watershed represents ~13% of the total drainage basin and is composed of a mixture of wide channels, flooded forests, clear water tributaries and floodplain lakes (Ward et al., 2015). The lower reach of the river is also subject to tidal effects, increasing the residence time of the water and increasing connectivity between flooded forests and the mainstem.
Previous studies have used δ18O-O2 values to estimate gas exchange (Jamieson et al., 2012) and estimate net metabolism (i.e., the ratio of respiration to photosynthesis) in steady-state systems (Quay et al., 1995; Bocaniov et al., 2012) as well as non-steady-state systems (Tobias et al., 2007; Holtgrieve et al., 2010). Past estimates of primary production in the Amazon River are based on a steady-state assumption (i.e., no diel variability in δ18O-O2) (Quay et al., 1995; Ellis et al., 2012); however, if primary production is actually occurring at appreciable rates in the mainstem of the river this steady-state assumption should not hold true.
This study seeks to examine the occurrence and magnitude of primary production within the lower Amazon River relative to respiration, from the historic downstream gaging station, Óbidos, to the river mouth. Measurements of the concentration and stable isotopic composition of dissolved oxygen made along the lower river and throughout a diel cycle were used to construct an oxygen mass balance for a discreet reach of the system in an effort to determine the relative contribution of primary production and respiration to ecosystem metabolism.
Materials and Methods
Measurements of the concentration and stable isotopic composition of dissolved oxygen were made along the lower Amazon River and coupled with measurements of microbial respiration rates (Ward et al., unpublished) and gas transfer velocities (Sawakuchi et al., unpublished) made during the same study period to construct an oxygen mass balance for the lower river. The stable isotopic composition of water (δ18O-H2O) was also measured for calculations associated with δ18O-O2.
Two expeditions were performed along the lower Amazon River during high water (May 2014) and low water (November 2015). The upstream boundary was Óbidos, the historic downstream gaging station, which is roughly 800 km from the river mouth (Figure 1). The downstream boundary was the northern and southern channels near the city of Macapá (NMCP and SMCP, respectively), which is ~650 km downstream of Óbidos. Sampling was also performed in the Amazon River mainstem midway to the mouth near the city of Almeirim and near the mouth of the two primary clearwater tributaries downstream of Óbidos, the Xingu and Tapajós rivers. The Tapajós and Xingu rivers add ~10% to the total discharge of the Amazon River (Sioli, 1985; Molleri et al., 2010). The lower Amazon River is characterized by a mixture of large channels, clear water tributaries, floodplain lakes and flooded forests, representing ~13% of the total Amazon River drainage basin (Ward et al., 2015). River discharge throughout the study period was reported by Ward et al. (2016) and is also shown in Table 1.
Figure 1. Samples were collected in the Amazon River mainstem at Óbidos, Almeirim, and the north and south channels near the city of Macapá. Samples were also collected in the Lagoa Grande de Curuai (high water only, inaccessible at low water), the Tapajós River, and the Xingu River. Adapted from Ward et al. (2015).
Table 1. River discharge, O2 saturation, the stable isotopic composition of dissolved O2, and the stable isotopic composition of H2O in the lower Amazon River.
Water was collected at three equidistant cross-channel sites at each of the five stations on the mainstem (Figure 1) using a Shurflo submersible pump with a 297 μm mesh screen at 50% river depth and from the surface (Ward et al., 2016). Sampling was only performed at one central station in the Tapajós and Xingu rivers at surface and 50% depth. Only one cross channel station was sampled in these tributaries due to their large channel width and relatively low depth, which results in little variability in chemical parameters across the channel. A graduated cylinder was continuously overflowed without bubbles using the submersible pump. Dissolved O2 concentrations were measured using a YSI ProOdo optical probe submerged in the overflowing cylinder and temperature was similarly measured using a Thermo Orion 4-star meter. Continuous measurements of dissolved O2 concentrations and temperature were also made while the ship was underway using a YSI Exo 2 sonde placed in an overflowing bucket with water pumped from the ship's intake.
Samples were collected for analysis of the stable isotopic composition of water (δ18O-H2O) in 15–50 mL centrifuge tubes and do not require preservation before analysis. Samples were collected for analysis of the stable isotopic composition of dissolved O2 in 12 mL exetainers (Barth et al., 2004). Four replicate samples were collected at 50% depth and the surface at all cross channel sites. To assess diel variability in O2 saturation and δ18O-O2, surface water samples were collected on a 3 hr interval for 15 h (evening to morning) at the center of the channel near Óbidos during late high water/early falling water (July, 2015).
Prior to sample collection, exetainers were rinsed with MilliQ water and combusted at 500°C. Fifty microliters of saturated mercuric chloride was added to each exetainer as a preservative. The exetainers were then dried and flushed with helium gas and capped in a pure helium glove box to remove any trace of atmospheric O2. Samples were collected by immersing the exetainers in the overflowing graduated cylinder. The vials were opened, filled completely, and then capped without removing them from the cylinder to prevent bubbles. To ensure the halt of all biological activity, an additional 50 μL of saturated mercuric chloride solution was injected through the septa. Septa were then coated with a thin layer of vacuum grease. Samples were analyzed within 2 months at the University of Washington Oceanography Stable Isotope Lab.
Analysis of δ18O-O2 and δ18O-H2O
For δ18O-O2 samples pure helium was pumped into the exetainer until half of the water was displaced 1 day prior to analysis. The vials were injected with 50 μL of 50% H3PO4 to convert all carbonate species to CO2. Isotopic ratios of headspace gases were determined by simultaneously measuring masses 32, 34, and 40 (i.e., 16O-16O, 18O-16O, and 40Ar) on a Finnegan Delta XL mass spectrometer (Thermo Electron). Dissolved oxygen concentrations were calculated based on the O2:Ar ratio (masses 32 and 40) and the estimated dissolved Ar concentration as a function of water temperature (Weiss, 1970). These values were compared to measurements taken with the YSI ProOdo and YSI Exo 2 sonde described above and were roughly equivalent (i.e., within ~5%). The stable isotopic composition of water (18O-H2O) was measured on 15 mL whole water samples using a Picarro L2130-I Analyzer.
Mass Balance Calculations
A dissolved oxygen mass balance was calculated for the Amazon River mainstem reach between Óbidos and midway to the mouth near the city of Almeirim (Figure 2). Two mass balance models were developed with the same set of inputs/outputs—one model uses only bulk parameters and the other includes the stable isotopic composition of O2 and associated fractionation factors.
Figure 2. Schematic for the mass balance of O2 and 18O-O2 between Óbidos and Almeirim based on measured parameters (green), estimated parameters (blue), and unknown parameters (red).
Inputs of O2 included the measured/calculated flux of O2 through Óbidos, measured/calculated O2 fluxes from the Tapajós River (Figure 1), calculated O2 gas exchange with the atmosphere (concentrations in the river were consistently below atmospheric saturation), estimated floodplain inputs, and primary production, which was the unknown parameter that was solved for. O2 outputs included microbial respiration rates measured by Ward et al. (unpublished) and the measured/calculated flux of O2 through Almeirim. In brief, respiration rates were measured in the dark in rotating incubation chambers interfaced to an YSI Exo 2 sonde equipped with an optical dissolved O2 probe. Respiration rates were calculated based on the linear drawdown of O2 during a roughly 24 hr period in duplicate 2.85 L chambers. Respiration rates include uncertainty associated with differences between replicates and averaged sites (±1 SD) and an additional 3% uncertainty due to the lack of inclusion of benthic respiration rates (more details on this in both the results and discussion sections).
This particular section of the river was chosen because the river maintains one well-constrained channel, for the most part, and river flow is not reversed by tides in this reach (although discharge is dampened semi-diurnally at Almeirim). This mass balance was not performed from Almeirim to Macapá because the river becomes extensively channelized and semi-diurnal tidal variation greatly complicates estimations of floodplain connectivity. The uncertainty associated with each mass balance parameter was propagated to our final calculation of primary production.
Calculation of O2 Gas Exchange
The molecular diffusivity of O2 across the air-water interface is described by the gas transfer velocity parameter (k). Gas transfer velocities can be calculated for other gases (and temperatures) based on the ratio of Schmidt numbers (Sc) (Jähne et al., 1987):
where kSc1 and kSc2 are the gas transfer velocities of a reference gas with a unknown and known value, respectively, and Sc1 and Sc2 are the respective Schmidt numbers for each gas. To determine the gas transfer velocity of O2, we converted calculated kCH4 values that were determined simultaneously (Sawakuchi, unpublished). We calculated Sc for O2 and CH4 at 28°C (measured mean temperature of the water) following Wanninkhof (1992) and used these values and the calculated kCH4 to determine kO2.
The flux of a gas across the air-water interface was calculated as follows:
where F is flux (mol m−2 d−1), k is the gas transfer velocity (m d−1), Cw is the concentration of gas measured in the water (mol m−3) and Ceq is the concentration of a water sample at equilibrium with the atmosphere at in situ temperature. Water temperatures ranged from 27.5 to 28.5°C over the course of our study. As such, we used a constant value of 28°C for all necessary equations. The equilibrium concentration of oxygen at 28°C was calculated using Henry's Law. Fluxes were calculated using varied k values calculated as described above, the calculated equilibration concentration and the average concentration of dissolved O2 for this reach of the river during high and low water periods. Uncertainty the gas transfer was calculated as the combination of analytical error and the standard deviation between observations at Óbidos and Almeirim (± 1 SD).
The water mass balance for this reach can be described by the following equation:
where QOBD is the discharge through Óbidos (m3 s−1), QTAP is discharge from the Tapajós River tributary, Qflood is discharge from floodplains, and QALM is discharge through Almeirim. Discharge was measured across the Amazon River main channel at Óbidos and Almeirim and across the mouth of the Tapajós River using a Sontek River Surveyor M9 Portable nine-beam 3.0 MHz/1.0 MHz/0.5 MHz acoustic Doppler Current Profiler (Ward et al., 2015). Discharge from floodplains was estimated by subtracting the measured discharge at Óbidos and the Tapajós River from the measured discharge at Almeirim. It should be noted that this does not necessarily represent the total amount of water that passed through floodplain complexes, but, rather, additional discharge added along this reach of the river from water stored in floodplains.
O2 Mass Balance
In order to calculate the non-measured/estimated oxygen inputs (i.e., primary production), two mass balance equations were used based on bulk O2 measurements (Equation 4) and both bulk and stable isotopic O2 measurements (Equation 5):
where A is the area of the river, F is the flux of O2 from the atmosphere to the river, V is the volume of the river, R is the average measured respiration rate, and P is the average rate of primary production, which was solved for. O2,ALM, O2,OBD, O2,TAP, and, O2,flood, are the measured concentrations of dissolved O2 at Almeirim, Óbidos, the Tapajós River, and the Lagoa Grande de Curuai floodplain lake, respectively. Uncertainty for the flux terms to/from Óbidos, the Tapajós River, and Almeirim are based on the analytical error of the O2 probe, variability in O2 concentrations with depth and across the channel, and discharge measurements (±1 SD).
Using equation 4 alone can give an estimate of inputs of O2 via primary production between Óbidos and Almeirim, but this can be improved by adding a mass balance for the stable isotopic composition of O2 as follows:
where 18:16O2 is the ratio of 18O to 16O calculated from δ18O-O2 values for each respective term relative to the Vienna Standard Mean Ocean Water (VSMOW) standard, 18:16Og and 18:16Ow are the stable isotopic compositions of atmospheric oxygen and river water, respectively. αR, αP, and αg are the fractionation factors of respiration, photosynthesis, and gas exchange respectively. A value of 0.982 was used for αR, a value of 0.997 was used for αg, and a value of 1.000 was used for αP considering O2 produced via primary production reflects the stable isotopic composition of the water (Quay et al., 1995). The above equations were solved for P, or the rate of primary production within this reach. Uncertainty from each model parameter (±1 SD) was propagated through to the final result. All calculations were performed based on the average temperature of 28°C.
Dissolved O2 Saturation
Dissolved oxygen was under-saturated relative to the atmosphere at all sampling locations at both high and low water with an annual average value of 79.3 ± 18.1% across the study boundaries. The highest and least variable O2 saturation levels were found in the Tapajós and Xingu clearwater rivers, with seasonal average values of 86.7 ± 1.0 and 96.0 ± 1.8% at high and low water, respectively (Table 1). The lowest O2 saturation levels were observed during high water in the Amazon River mainstem, specifically at Óbidos. Across all mainstem sites including the mouth, average O2 saturation was 54.3 ± 6.0% at high water compared to 88.5 ± 6.5% at low water (Table 1).
O2 saturation in the mainstem exhibited a downstream increase from Óbidos to the river mouth with an annual average value of 62.8 ± 23.2% at Óbidos, 71.0 ± 25.4% at Almeirim, and 76.0 ± 19.6%. The downstream increase was most evident during low water, with values increasing from 79.2 ± 1.9% at Óbidos to 93.0 ± 0.3% across the mouth. During the high water period O2 saturation increased from 46.4 ± 4.8% at Óbidos to 59.0 ± 1.3% across the mouth (Table 1).
The partial pressure of CO2 was previously measured along the same study boundaries (Sawakuchi et al., under review; Ward et al., 2016). There was a negative correlation between the molar concentration of dissolved O2 and CO2 across the study boundaries (Figure 3). Spatial variability of O2 and CO2 in the entire study region, including both the Amazon River mainstem and tributary sites, was higher during the high water period. For example, CO2 concentrations ranged from 86.9 to 200.7 μmol L−1 and 36.7 to 61.5 μmol L−1 during high and low water, respectively (Figure 3). O2 concentrations varied from 90.0 to 150.8 μmol L−1 and 160.3 to 197.2 μmol L−1 at high and low water, respectively, across our study boundaries. The slope of the relationship between O2 and CO2 concentrations was 2.2 times higher during low water.
Dissolved O2 Stable Isotopic Composition
The oxygen stable isotopic composition of water (i.e., δ18O-H2O) is influenced by isotopic fractionation during each phase of the hydrologic cycle. Annually-averaged δ18O-H2O was −4.3 ± 1.1‰ across the entire study boundary. δ18O-H2O was the most negative in the Amazon River mainstem at high water with an average value of −5.4 ± 2.1‰ compared to −3.6 ± 0.4‰ at low water (Table 1). δ18O-H2O was more positive in the clearwater tributaries (average = −3.6 ± 1.2‰) compared to the Amazon River mainstem (average = −4.5 ± 1.0‰). There was also a downstream increase in δ18O-H2O from Óbidos to the river mouth. For example, at high water an average δ18O-H2O value of −5.7 ± 0.1‰ was observed at Óbidos compared to −5.2 ± 0.1‰ across the river mouth. The difference between the mainstem and tributaries is likely due to differences in the isotopic composition of rainfall, since much of the water flowing through the mainstem is derived from far upstream.
The annual average stable isotopic composition of dissolved oxygen (i.e., δ18O-O2) was 23.8 ± 1.4‰ across the entire study boundary. The highest δ18O-O2 values were observed in the Amazon River mainstem at high water with an average value of 25.0 ± 1.1‰ compared to 23.9 ± 0.7‰ at low water. The clearwater tributaries had consistently lower and less seasonally variable δ18O-O2 values, varying from 22.9 ± 1.5‰ at high water to 22.2 ± 1.7‰ at low water (Table 1). The Xingu River had the lowest δ18O-O2 values across all sites (annual average = 21.5 ± 0.6‰). δ18O-O2 and oxygen saturation levels were negatively correlated during both sampling seasons (Figure 4). Higher δ18O-O2 values were found at sites with lower oxygen saturation levels and vice versa. The regression slope between δ18O-O2 values and oxygen saturation, while following the same trend, was steeper for samples collected at low water compared to water.
It is typically assumed that d18O-O2 should remain at steady state in the Amazon River mainstem considering the low levels of light penetration and limited primary production (Quay et al., 1995; Ellis et al., 2012). To test this assumption we measured O2 saturation and d18O-O2 throughout a day/night cycle in the center of the channel at Óbidos. While the overall oxygen saturation levels did not vary significantly relative to the signal noise over a 15 hr period, our results indicate that d18O-O2 does follow a diel cycle even in the Amazon River mainstem (Figure 5). The minimum d18O-O2 was observed at 19:00, roughly 1 h after sunset (24.7 ± 0.6‰) and the maximum value was observed at 06:00 just prior to sunset (26.4 ± 0.3‰), followed by a decrease to 25.3 ± 0.7‰ at 08:00 (Figure 5). The difference in d18O-O2 values observed between midnight and 06:00 and the rest of the diel cycle were statistically significant based on an unpaired t-test within a 95% confidence interval (p = 0.02). It should be noted that a complete 24 h cycle was not able to be captured due to logistical constraints of the field campaign.
Figure 5. Diel O2 saturation (red) and δ18O-O2 (blue) measurements made in the center of the channel at Obidos. The shaded portion of the chart represents hours between sunset and sunrise.
Oxygen Mass Balance Model
In the past, δ18O-O2 values have been used to calculate the ratio of photosynthesis to respiration (P:R) under a steady state assumption for the Amazon River (Quay et al., 1995; Ellis et al., 2012). However, considering diel variability of d18O-O2 was observed, a steady state assumption is not accurate. Thus, in order to determine the rates of photosynthetic oxygen input into the mainstem of the lower Amazon River, we used a mass balance model that takes into account the various sources and sinks of O2 within the system (Figure 2). Because the stable isotopic composition of dissolved O2 is also known, we were able to calculate two separate O2 budgets to model the system based on the following: (1) bulk O2 concentrations only (i.e., Equation 4) and (2) the concentration and stable isotopic composition of O2 (i.e., Equation 5). Including d18O-O2 observations and literature-derived fractionation factors for each process expressed in Equation 5 resulted in a slight increase in each model parameter compared to the bulk O2 mass balance, however, the resulting calculated P:R ratios were essentially equivalent (Table 2). For this reason, we report the average value for each parameter from these two different models below.
Table 2. Dissolved O2 fluxes through Óbidos, Almeirim, floodplains, and the Tapajós River and from respiration, gas exchange, and primary production using mass balance models derived from (1) bulk O2 concentrations and (2) O2 stable isotopic compositions.
Although O2 concentrations were considerably lower in the Amazon River mainstem during high water (Table 1), the flux of O2 through Óbidos and Almeirim were greater due to higher river discharge rates. The average flux of O2 through Óbidos at high water, determined using both models, was 901 ± 87 kg O s−1 compared to 739 ± 80 kg O s−1 at low water (Table 2). The flux of O2 through Almeirim was consistently greater than at Óbidos, with values of 1228 ± 112 kg O s−1 at high water and 853 ± 83 kg O s−1. This difference was greatest during high water with a 36% increase in O2 flux from Óbidos to Almeirim compared to a 15% increase during the low water period.
Aside from primary production, gas exchange was the dominant source of O2 to this reach of the river, resulting in an inward flux of 722 ± 169 kg O s−1 at high water and 565 ± 130 kg O s−1 at low water. The Tapajós River provided an additional 81 ± 10 kg O s−1 to the Amazon River mainstem during high water and 25 ± 3 kg O s−1 during low water, which represents about 20–25% of the difference in fluxes between Óbidos and Almeirim during both seasons. We estimate that the flux of O2 from adjacent floodplains to the mainstem was roughly 247 ± 75 kg O s−1 during high water (Table 2). At low water, we considered floodplain discharge to be zero considering that river discharge at Almeirim was lower than the sum discharge from Óbidos and the Tapajós River, indicative of water entering and being stored in the floodplains. Further, the Lagoa Grande de Curuai was inaccessible during this sampling period, illustrating the limited floodplain connectivity to the mainstem during this time.
Aside from the flux through Almeirim, respiration was the primary O2 output. The average water column respiration rate measured across the Amazon River mainstem sampling sites was 0.95 ± 0.10 g O m3 d−1 at high water and 1.52 ± 0.35 g O m3 d−1 at low water (Ward et al., unpublished). Benthic respiration rates were not measured, however, literature data is available from sites further upstream. Average benthic respiration rates of 1.46 g O m2 d−1 have been reported for the central Amazon River (Devol et al., 1987). Considering that the average river depth is ~40 m between Óbidos and Almeirim (Table 1), benthic respiration rates represent only 3% of depth-integrated water column respiration rates. Since benthic respiration was not measured in our study region, we have included this 3% error to the uncertainty in our total respiration estimate rather than adding a separate model parameter similar to calculations made by Devol et al. (1987). Water column respiration rates were multiplied by the volume of water between Óbidos and Almeirim (1.09 × 1011 m3 at high water and 8.17 × 1010 m3 at low water), resulting in a consumption of 1228 ± 112 kg O s−1 at high water and 1447 ± 368 kg O s−1 at low water (Table 2).
With all O2 inputs and outputs other than primary production either measured or estimated we then solved for primary production as the difference between the mass balance outputs and inputs. Primary production was estimated to add 492 ± 291 kg O s−1 to the Amazon River mainstem between Óbidos and Almeirim at high water and 971 ± 407 kg O s−1. This corresponds to a P:R ratio of 0.41 ± 0.24 at high water, 0.67 ± 0.33 at low water, and 0.54 ± 0.41 annually-averaged. When multiplied by the measured respiration rates, this corresponds to a primary production rate of 0.39 ± 0.24 g O m3 d−1 at high water, 1.02 ± 0.55 g O m3 d−1 at low water, and 0.70 ± 0.58 g O m3 d−1 annually-averaged.
Dissolved Oxygen Balance
The stable isotopic composition of dissolved O2 (i.e., δ18O-O2) is primarily regulated by three processes: air-water gas exchange, respiration, and photosynthesis. In a system purely affected by gas-exchange, the δ18O-O2 would be 24.2‰ (relative to VSMOW) due to the dissolution of atmospheric O2 (23.5‰) and the equilibrium fractionation of 0.7‰ that occurs during gas dissolution (Kroopnick, 1975; Benson and Krause, 1984). Photosynthesis produces O2 with a δ18O of the surrounding water because there is little fractionation observed during photosynthesis (Guy et al., 1993). Respiration selectively fractionates oxygen by consuming O2 that is more depleted than the average, which enriches the remaining pool of O2. Thus, photosynthesis is expected to deplete the δ18O-O2 of a system while respiration will enrich the δ18O-O2.
Bulk O2 levels were consistently undersaturated relative to the atmosphere across the entire study boundary. Likewise, δ18O-O2 values were consistently more enriched in the Amazon River mainstem relative to values expected due to gas transfer with the atmosphere alone (i.e., 24.2‰; Table 1), indicating a predominance of respiration relative to photosynthesis. These factors both reflect the net heterotrophic state previously observed in both the Amazon River (Mayorga et al., 2005) and river systems worldwide (Cole et al., 2007), which results in a net flux of CO2 to the atmosphere (Richey et al., 2002; Alin et al., 2011; Raymond et al., 2013).
However, these bulk and stable isotopic O2 observations do not preclude the input of photosynthetically-produced oxygen into the mainstem. In fact, our observations of diel variability in δ18O-O2 values at Óbidos clearly indicate a shift in the balance between primary production and respiration in the mainstem throughout a day/night cycle. Our observations were consistent with daytime production depleting δ18O-O2 and nighttime respiration enriching the signal, implying that O2 does not remain at a steady-state as previously assumed (Quay et al., 1995; Wassenaar, 2012).
The more depleted δ18O-O2 values observed in the Amazon River mainstem during low water indicate a relatively higher amount of primary production relative to respiration during low water. This makes sense considering that suspended sediment concentrations are generally higher during high water in the Amazon River mainstem, limiting primary production (Hedges and Clark, 1986; Moreira-Turcq et al., 2003, 2013).
Respiration rates also vary seasonally, influencing the ratio of production to respiration. The highest respiration rates observed here and in the literature occur during the low water period (Benner et al., 1995; Ellis et al., 2012; Ward et al., 2013), which implies that photosynthesis rates would need to be proportionally higher to achieve a shift in P:R ratios during the low water period. In fact, our mass balance model predicted this difference. We estimate that primary production was 2.0 times higher in the Óbidos to Almeirim reach during low water based on our mass balance models (Table 2). Average measured respiration rates were 1.6 times greater in the mainstem during low water, but when integrated across the entire volume of water in the Óbidos to Almeirim reach, the total respiration flux was only 1.2 times greater at low water compared to high water. These differences result in the observed shift in the ratio of P:R from 0.41 ± 0.24 at high water to 0.67 ± 0.33 at low water.
O2 levels were closer to atmospheric saturation in the Tapajós and Xingu rivers (Table 1), and δ18O-O2 values were below 24.2‰ (Table 2), indicating a predominance of primary production relative to respiration in the less turbid clearwater tributaries. Sediment levels are significantly lower in these tributaries compared to the Amazon River mainstem, and high rates of primary production are reflected in elevated chlorophyll a concentrations relative to the Amazon River (Moreira-Turcq et al., 2013; Ward et al., 2015, 2016). Interestingly, respiration rates have been observed to be higher in clearwater tributaries than in the Amazon River mainstem throughout the Amazon basin (Benner et al., 1995; Ellis et al., 2012), yet these tributaries are overall less net heterotrophic. The elevated respiration rates observed in clearwater tributaries have been attributed to both the breakdown of algal organic matter (Benner et al., 1995; Ellis et al., 2012) and the enhanced breakdown of terrestrially-derived organic matter due to priming effects, i.e., the enhanced breakdown of a recalcitrant substrate due to the presence of a labile substrate (Bianchi et al., 2015; Ward et al., 2016).
Previous studies have used δ18O-O2 values to estimate gas exchange(Jamieson et al., 2012) and estimate net metabolism (e.g., the P:R ratio) in steady-state systems(Quay et al., 1995; Bocaniov et al., 2012) as well as non-steady-state systems (Tobias et al., 2007; Holtgrieve et al., 2010; Hotchkiss and Hall, 2014). For the Amazon River mainstem it has been estimated that the ratio of P:R is roughly 0.25 based on δ18O-O2 values and a steady state assumption (Quay et al., 1995), which is 2.2 times lower than our annual average P:R estimate. Primary production rates have also been estimated based on 14C-NaHCO3 uptake experiments; the average rate of primary production was estimated to be 0.10 ± 0.06 g O m3 d−1 in the Amazon River mainstem and 1.45 ± 0.75 g O m3 d−1 in the Tapajós and Xingu rivers (Wissmar et al., 1981). This is roughly 6.7 times lower than our mass balance estimates and 2.9 lower than estimates based on multiplying our measured rates of respiration by P:R ratios determined by Quay et al. (1995). We argue that this difference is due to the fact that bottle experiments cannot adequately capture the complex ecosystem dynamics occurring along the entire river reach, whereas a mass balance integrates these signals. Our mass balance results are in agreement with observations of a significant abundance of phytoplankton biomarkers in the mainstem of the Amazon (Mortillaro et al., 2011). Likewise, primary production has recently been found to be important in the main channel of the Congo River (Descy et al., 2017).
Limitations and Future Considerations
One factor that remains difficult to constrain is whether this input of isotopically depleted O2 is conclusively from in situ production in the mainstem itself or is the legacy of advection from lateral floodplain lakes. These floodplain lakes have high rates of productivity and there is evidence of their legacy biomass in the Amazon mainstem (de Moraes Novo et al., 2006; Abril et al., 2014). This input of floodplain-derived oxygen is reflected in the ratio of CO2 to O2 measured at high and low water (Figure 3). At high water, the ratio is higher, likely due to the input of root respiration-derived CO2 from the floodplains. At high water, these flooded macrophytes draw down oxygen from the atmosphere to their roots, thereby producing CO2 without the concomitant drawdown of O2 from the water column. Thus, it is possible that the isotopically depleted O2 signals observed at high water is actually the legacy of floodplain production and not mainstem production. However, at low water, there is little evidence of macrophyte root respiration as the river loses connectivity with the floodplains. Thus, it is likely the production signal observed at low water is primarily from in situ production within the mainstem rather than floodplains.
Another factor that we did not consider is the importance of benthic respiration from the riverbed. However, past estimates suggest that benthic respiration in the OM-poor sandy sediments only accounts for ~3% of our measured water column respiration rates (Devol et al., 1987). Thus we did not include a benthic respiration term in our mass balance, but rather added 3% error to our total respiration term. In some estuarine settings benthic respiration has been shown to be the dominant factor controlling dissolved oxygen concentrations and stable isotopic compositions (Lehmann et al., 2009), however this does not appear to be the case in the fast-flowing, deep river setting where reactive organic matter is highly abundant in the water column but relatively depleted in sediments (Richey et al., 1990; Hedges et al., 2000; Moreira-Turcq et al., 2013).
We did not attempt this same mass balance calculation for the river reach between Almeirim and Macapá due to the complex nature of this region. The river evolves into a series of sub-channels in the lower reaches limiting the applicability of a simplistic model, whereas the reach between Óbidos and Almeirim represents a single, well-constrained channel for the most part. Further, semi-diurnal flow reversals due to tides make it even more difficult to quantify connectivity and import/export from floodplain networks near the river mouth. Future efforts will be applied to unraveling the complex tidally-influenced reach of the river. It should also be noted the floodplain discharge used in this mass balance does not represent the total amount of water that passed through floodplain systems, but, rather, excess water entering the river from floodplains between Óbidos and Almeirim. It is difficult to accurately constrain the amount of water that passes through the Lagoa Grande de Curuai, for example, but this will be the focus of future hydrodynamic modeling efforts to further constrain our mass balance estimates.
Early investigations of biogeochemical cycling in the Amazon River found minimal rates of primary production in the turbid mainstem relative to the productive clearwater tributaries (Fisher, 1979; Wissmar et al., 1981). Here, we have shown that there is evidence for primary production within the Amazon River mainstem based on an O2 mass balance, which exceeds past estimates based on a steady-state stable isotopic model by 2.0 times (Quay et al., 1995) and is 6.7 times higher than estimates from radiocarbon uptake experiments (Wissmar et al., 1981). The Amazon River is a major source of CO2 to the atmosphere, but understanding the interplay between photosynthesis and respiration is critical for understanding the fundamental mechanisms driving these fluxes and the overall productivity of the ecosystem. Likewise, in situ primary production has been shown to be an important factor in elevating both bulk respiration rates (Ellis et al., 2012) and the breakdown of terrestrially-derived organic matter (Ward et al., 2016). Although this result does not undermine the idea of the Amazon River as a net heterotrophic system, it does indicate that more work remains to understand the dynamic interplay between production and respiration in large turbid river systems.
WG performed the collection and analysis of δ18O-O2 samples and developed the O2 mass balance. VN, NW, DD, and JD. performed measurements of dissolved O2 concentrations. VN and HS. performed estimations of gas transfer velocities used in the mass balance. NW and WG. performed measurements of respiration used in the O2 mass balance. JR, AD, and DB performed measurements of river discharge. The field sampling scheme was developed and led by JR, AK, RK, NW, and HS. Field and sampling logistics were undertaken by all authors. All authors contributed to the preparation of the manuscript and approved its final submission.
This study was supported by FAPESP Grant # 08/58089-9 and NSF DEB Grant # 1256724. WGM was funded by an NSF IGERT grant DGE-1258485.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank the crew of the B/M Mirage for contributions made during the river cruises, Dr. Gordon Holtgrieve for his advice on oxygen isotopes and the Quay Lab at UW for assistance in isotopic analysis.
Abril, G., Martinez, J. M., Artigas, L. F., Moreira-Turcq, P., Benedetti, M.F., Vidal, L., et al. (2014). Amazon River carbon dioxide outgassing fuelled by wetlands. Nature 505, 395–398. doi: 10.1038/nature12797
Alin, S. R., Rasera, M. D. F. F. L., Salimon, C. I., Richey, J. E., Holtgrieve, G. W., Krusche, A. V., et al. (2011). Physical controls on carbon dioxide transfer velocity and flux in low-gradient river systems and implications for regional carbon budgets. J. Geophys. Res. Biogeosciences 116:G01009. doi: 10.1029/2010jg001398
Benner, R., Opsahl, S., Chin-Leo, G., Richey, J. E., and Forsberg, B. R. (1995). Bacterial carbon metabolism in the Amazon River system. Limnol. Oceanogr. 40, 1262–1270. doi: 10.4319/lo.19184.108.40.2062
Benson, B. B., and Krause, D. (1984). The concentration and isotopic fractionation of oxygen dissolved in freshwater and seawater in equilibrium with the atmosphere. Limnol. Oceanogr. 29, 620–632. doi: 10.4319/lo.1984.29.3.0620
Bianchi, T. S., Thornton, D. C. O., Yvon-lewis, S. A., King, G. M., Eglinton, T. I., Shields, M. R., et al. (2015). Positive priming of terrestrially derived dissolved organic matter in a freshwater microcosm system. Geophys. Res. Lett. 42, 5460–5467. doi: 10.1002/2015GL064765
Bocaniov, S. A., Schiff, S. L., and Smith, R. E. H. (2012). Plankton metabolism and physical forcing in a productive embayment of a large oligotrophic lake : insights from stable oxygen isotopes. Freshw. Biol. 57, 481–496. doi: 10.1111/j.1365-2427.2011.02715.x
Cole, J. J., Prairie, Y. T., Caraco, N. F., McDowell, W. H., Tranvik, L. J., Striegl, R. G., et al. (2007). Plumbing the Global carbon cycle: integrating Inland waters into the terrestrial carbon budget. Ecosystems 10, 172–185. doi: 10.1007/s10021-006-9013-8
del Giorgio, P. A., and Pace, M. L. (2008). Relative independence of organic carbon transport and processing in a large temperate river: the Hudson River as both pipe and reactor. Limnol. Oceanogr. 53, 185–197. doi: 10.4319/lo.2008.53.1.0185
de Moraes Novo, E. M. L., de Farias Barbosa, C. C., de Freitas, R. M., Shimabukuro, Y. E., Melack, J. M., Filho, W. P., et al. (2006). Seasonal changes in chlorophyll distributions in Amazon floodplain lakes derived from MODIS images. Limnol. 7, 153–161. doi: 10.1007/s10201-006-0179-8
Descy, J. P., Darchambeau, F., Lambert, T., Stoyneva-Gaertner, M. P., Bouillon, S., and Borges, A. V. (2017). Phytoplankton dynamics in the Congo River. Freshw. Biol. 62, 87–101. doi: 10.1111/fwb.12851
Devol, A. H., Quay, P. D., and Richey, J. E. (1987). The role of gas exchange in the inorganic carbon, oxygen, and 222Rn budgets of the Amazon River. Limnol. Oceanogr. 32, 235–248. doi: 10.4319/lo.1987.32.1.0235
Ellis, E. E., Richey, J. E., Aufdenkampe, A. K., Krusche, A. V., Quay, P. D., Salimon, C., et al. (2012). Factors controlling water-column respiration in rivers of the central and southwestern Amazon Basin. Limnol. Oceanogr. 57, 527–540. doi: 10.4319/lo.2012.57.2.0527
Guenet, B., Danger, M., Harrault, L., Allard, B., Jauset-Alcala, M., Bardoux, G., et al. (2014). Fast mineralization of land-born C in inland waters: first experimental evidences of aquatic priming effect. Hydrobiologia 721, 35–44. doi: 10.1007/s10750-013-1635-1
Hedges, J. I., Mayorga, E., Tsamakis, E., McClain, M. E., Aufdenkampe, A., Quay, P., et al. (2000). Organic matter in Bolivian tributaries of the Amazon River: a comparison to the lower mainstream. Limnol. Oceanogr. 45, 1449–1466. doi: 10.4319/lo.2000.45.7.1449
Hoffman, J. C., Bronk, D. A., and Olney, J. E. (2008). Organic matter sources supporting lower food web production in the tidal freshwater portion of the York River estuary, Virginia. Estuar. Coast. 31, 898–911. doi: 10.1007/s12237-008-9073-4
Holtgrieve, G. W., Schindler, D. E., Branch, T. A., and Teresa, Z. A. (2010). Simultaneous quantification of aquatic ecosystem metabolism and reaeration using a Bayesian statistical model of oxygen dynamics. Limnol. Oceanogr. 55, 1047–1063. doi: 10.4319/lo.2010.55.3.1047
Hotchkiss, E. R., and Hall R. O. Jr. (2014). High rates of daytime respiration in three streams: Use of δ18OO2 and O2 to model diel ecosystem metabolism. Limnol. Oceanogr. 59, 798–810. doi: 10.4319/lo.2014.59.3.0798
Jähne, B. J., Münnich, K. O. M., Bösinger, R., Dutzi, A., Huber, W., and Libner, P. (1987). On the parameters influencing air-water gas exchange. J. Geophys. Res. 92, 1937–1949. doi: 10.1029/JC092iC02p01937
Jamieson, T. S., Schiff, S. L., and Taylor, W. D. (2012). Using stable isotopes of dissolved oxygen for the determination of gas exchange in the Grand River, Ontario, Canada. Water Res. 47, 781–790. doi: 10.1016/j.watres.2012.11.001
Junk, W. J., Soares, M. G. M., and Bayley, P. B. (2007). Freshwater fishes of the Amazon River basin: their biodiversity, fisheries, and habitats. Aquat. Ecosyst. Health Manag. 10, 153–173. doi: 10.1080/14634980701351023
Kritzberg, E. S., Cole, J. J., Pace, M. L., Grane, W., and Bade, D. L. (2004). Autochthonous versus allochthonous carbon sources of bacteria : results from whole-lake 13C addition experiments. Limnol. Oceanogr. 49, 588–596. doi: 10.4319/lo.2004.49.2.0588
Lehmann, M. F., Barnett, B., Gélinas, Y., Gilbert, D., Maranger, R. J., Mucci, A., et al. (2009). Aerobic respiration and hypoxia in the Lower St. Lawrence Estuary: stable isotope ratios of dissolved oxygen constrain oxygen sink partitioning. Limnol. Oceanogr. 54, 2157–2169. doi: 10.4319/lo.2009.54.6.2157
Mayorga, E., Aufdenkampe, A. K., Masiello, C., A Krusche, A. V., Hedges, J. I., Quay, P. D., et al. (2005). Young organic matter as a source of carbon dioxide outgassing from Amazonian rivers. Nature 436, 538–541. doi: 10.1038/nature03880
Moreira-Turcq, P., Bonnet, M. P., Amorim, M., Bernardes, M., Lagane, C., Maurice, L., et al. (2013). Seasonal variability in concentration, composition, age, and fluxes of particulate organic carbon exchanged between the floodplain and Amazon River. Glob. Biogeochem. Cycles 27, 119–130. doi: 10.1002/gbc.20022
Moreira-Turcq, P., Seyler, P., Guyot, J. L., and Etcheber, H. (2003). Exportation of organic carbon from the Amazon River and its main tributaries. Hydrol. Process. 17, 1329–1344. doi: 10.1002/hyp.1287
Mortillaro, J. M., Abril, G., Moreira-Turcq, P., Sobrinho, R. L., Perez, M., and Meziane, T. (2011). Fatty acid and stable isotope (13C, 15N) signatures of particulate organic matter in the lower Amazon River: seasonal contrasts and connectivity between floodplain lakes and the mainstem. Org. Geochem. 42, 1159–1168. doi: 10.1016/j.orggeochem.2011.08.011
Quay, P. D., Wilbur, D., Richey, J. E., Devol, A. H., Benner, R. and Forsberg, B. R. (1995). The 18O: 16O of dissolved oxygen in rivers and lakes in the Amazon Basin: determining the ratio of respiration to photosynthesis rates in freshwaters. Limnol. Oceanogr. 40, 718–729.
Quay, P. D., Wilbur, D. O., Richey, J. E., Hedges, J. I., Devol, A. H., and Victoria, R. (1992). Carbon cycling in the Amazon River: implications from the 13C compositions of particles and solutes. Limnol. Oceanogr. 37, 857–871. doi: 10.4319/lo.1992.37.4.0857
Richey, J. E., Hedges, J. I., Devol, A. H., Quay, P. D., Victoria, R., Martinelli, L., et al. (1990). Biogeochemistry of carbon in the Amazon River. Limnol. Oceanogr. 35, 352–371. doi: 10.4319/lo.1990.35.2.0352
Richey, J. E., Melack, J. M., Aufdenkampe, A. K., Ballester, V. M., and Hess, L. L. (2002). Outgassing from Amazonian rivers and wetlands as a large tropical source of atmospheric CO2. 6416, 6413–6416. doi: 10.1038/416617a
Tobias, C. R., Böhlke, J. K., and Harvey, J. W. (2007). The oxygen-18 isotope approach for measuring aquatic metabolism in high productivity waters. Limnol. Oceanogr. 52, 1439–1453. doi: 10.4319/lo.2007.52.4.1439
Tranvik, L. J., Downing, J. A., Cotner, J. B., Loiselle, S. A., Striegl, R. G., Ballatore, T. J., et al. (2009). Lakes and reservoirs as regulators of carbon cycling and climate. Limnol. Oceanogr. 54, 2298–2314. doi: 10.4319/lo.2009.54.6_part_2.2298
Ward, N. D., Bianchi, T. S., Sawakuchi, H. O., Gagne-maynard, W., Cunha, A. C., Brito, D. C., et al. (2016). The reactivity of plant-derived organic matter and the potential importance of priming effects along the lower Amazon River. J. Geophys. Res. Biogeosci. 121, 1522–1539. doi: 10.1002/2016jg003342
Ward, N. D., Keil, R. G., Medeiros, P. M., Brito, D. C., Cunha, A. C., Dittmar, T., et al. (2013). Degradation of terrestrially derived macromolecules in the Amazon River. Nat. Geosci. 6, 530–533. doi: 10.1038/ngeo1817
Ward, N. D., Krusche, A. V., Sawakuchi, H. O., Brito, D. C., Cunha, A. C., Moura, J. M. S., et al. (2015). The compositional evolution of dissolved and particulate organic matter along the lower Amazon River-Óbidos to the ocean. Mar. Chem. 177, 244–256. doi: 10.1016/j.marchem.2015.06.013
Keywords: photosynthesis, tropical rivers, respiration, oxygen, mass balance
Citation: Gagne-Maynard WC, Ward ND, Keil RG, Sawakuchi HO, Da Cunha AC, Neu V, Brito DC, Da Silva Less DF, Diniz JEM, De Matos Valerio A, Kampel M, Krusche AV and Richey JE (2017) Evaluation of Primary Production in the Lower Amazon River Based on a Dissolved Oxygen Stable Isotopic Mass Balance. Front. Mar. Sci. 4:26. doi: 10.3389/fmars.2017.00026
Received: 02 November 2016; Accepted: 24 January 2017;
Published: 07 February 2017.
Edited by:Hongbin Liu, The Hong Kong University of Science and Technology, Hong Kong
Reviewed by:Benoit Thibodeau, University of Hong Kong, Hong Kong
Bingzhang Chen, Japan Agency for Marine-Earth Science and Technology, Japan
Copyright © 2017 Gagne-Maynard, Ward, Keil, Sawakuchi, Da Cunha, Neu, Brito, Da Silva Less, Diniz, De Matos Valerio, Kampel, Krusche and Richey. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Nicholas D. Ward, email@example.com