Fine-scale sampling unveils diazotroph patchiness in the South Pacific Ocean.

Diazotrophs are important contributors to nitrogen availability in the ocean. Oceanographic cruise data accumulated over the past three decades has revealed a heterogeneous distribution of diazotroph species at regional to global scales. However, dynamic fine-scale physical structures likely affect the distribution of diazotrophs at smaller spatiotemporal scales. The interaction between fine-scale ocean dynamics and diazotrophs remains poorly understood due to typically insufficient spatiotemporal sampling resolution and the lack of parallel detailed physical studies. Here we show the distribution of five groups of diazotrophs in the South Pacific at an unprecedented resolution of 7-16 km. We find a patchy distribution of diazotrophs, with each group being differentially affected by parameters describing fine-scale physical structures. The observed variability in species abundance and distribution would be masked by a coarser sampling resolution, highlighting the need to consider fine-scale physics to resolve the distribution of diazotrophs in the ocean.

Optical proxy for particulate organic nitrogen from BGC-Argo floats.

Using biogeochemical-Argo float measurements, we propose, for the first time, an optical proxy for particulate organic nitrogen concentration (PON) in the Western Tropical South Pacific, an area influenced by dinitrogen (N2) fixation. Our results show a significant relationship between the backscattering coefficient at 700 nm (bbp) and PON, especially when the latter is measured using the wet oxidation method (R2=0.87). bbp may be used to estimate PON concentrations (PONopt) between 0.02 and 0.95 µM, allowing for unprecedented monitoring using autonomous profiling floats. The bbpvs PON relationship can be used to study phytoplanktonic biomass dynamics at relevant seasonal temporal scales, with clear evidence of PONopt as a proxy of phytoplanktonic biomass, at least for this specific area. Temporal analyses of PONopt show significant increases (from 0.16 to 0.80 µM) likely related to new production associated to N2 fixation events measured during stratification periods in the Melanesian Archipelago.

Ecological networks: Pursuing the shortest path, however narrow and crooked.

Representing data as networks cuts across all sub-disciplines in ecology and evolutionary biology. Besides providing a compact representation of the interconnections between agents, network analysis allows the identification of especially important nodes, according to various metrics that often rely on the calculation of the shortest paths connecting any two nodes. While the interpretation of a shortest paths is straightforward in binary, unweighted networks, whenever weights are reported, the calculation could yield unexpected results. We analyzed 129 studies of ecological networks published in the last decade that use shortest paths, and discovered a methodological inaccuracy related to the edge weights used to calculate shortest paths (and related centrality measures), particularly in interaction networks. Specifically, 49% of the studies do not report sufficient information on the calculation to allow their replication, and 61% of the studies on weighted networks may contain errors in how shortest paths are calculated. Using toy models and empirical ecological data, we show how to transform the data prior to calculation and illustrate the pitfalls that need to be avoided. We conclude by proposing a five-point check-list to foster best-practices in the calculation and reporting of centrality measures in ecology and evolution studies.

Iron from a submarine source impacts the productive layer of the Western Tropical South Pacific (WTSP).

In the Western Tropical South Pacific, patches of high chlorophyll concentrations linked to the occurrence of N2-fixing organisms are found in the vicinity of volcanic islands. The survival of these organisms relies on a high bioavailable iron supply whose origin and fluxes remain unknown. Here, we measured high dissolved iron (DFe) concentrations (up to 66 nM) in the euphotic layer, extending zonally over 10 degrees longitude (174 E-175 W) at ∼20°S latitude. DFe atmospheric fluxes were at the lower end of reported values of the remote ocean and could not explain the high DFe concentrations measured in the water column in the vicinity of Tonga. We argue that the high DFe concentrations may be sustained by a submarine source, also characterized by freshwater input and recorded as salinity anomalies by Argo float in situ measurements and atlas data. The observed negative salinity anomalies are reproduced by simulations from a general ocean circulation model. Submarine iron sources reaching the euphotic layer may impact nitrogen fixation across the whole region.

On the calculation of betweenness centrality in marine connectivity studies using transfer probabilities.

Betweenness has been used in a number of marine studies to identify portions of sea that sustain the connectivity of whole marine networks. Herein we highlight the need of methodological exactness in the calculation of betweenness when graph theory is applied to marine connectivity studies based on transfer probabilities. We show the inconsistency in calculating betweeness directly from transfer probabilities and propose a new metric for the node-to-node distance that solves it. Our argumentation is illustrated by both simple theoretical examples and the analysis of a literature data set.

A connectivity-based eco-regionalization method of the Mediterranean Sea.

Ecoregionalization of the ocean is a necessary step for spatial management of marine resources. Previous ecoregionalization efforts were based either on the distribution of species or on the distribution of physical and biogeochemical properties. These approaches ignore the dispersal of species by oceanic circulation that can connect regions and isolates others. This dispersal effect can be quantified through connectivity that is the probability, or time of transport between distinct regions. Here a new regionalization method based on a connectivity approach is described and applied to the Mediterranean Sea. This method is based on an ensemble of Lagrangian particle numerical simulations using ocean model outputs at 1/12° resolution. The domain is divided into square subregions of 50 km size. Then particle trajectories are used to quantify the oceanographic distance between each subregions, here defined as the mean connection time. Finally the oceanographic distance matrix is used as a basis for a hierarchical clustering. 22 regions are retained and discussed together with a quantification of the stability of boundaries between regions. Identified regions are generally consistent with the general circulation with boundaries located along current jets or surrounding gyres patterns. Regions are discussed in the light of existing ecoregionalizations and available knowledge on plankton distributions. This objective method complements static regionalization approaches based on the environmental niche concept and can be applied to any oceanic region at any scale.

Optical Algorithms at Satellite Wavelengths for Total Suspended Matter in Tropical Coastal Waters.

Is it possible to derive accurately Total Suspended Matter concentration or its proxy, turbidity, from remote sensing data in tropical coastal lagoon waters? To investigate this question, hyperspectral remote sensing reflectance, turbidity and chlorophyll pigment concentration were measured in three coral reef lagoons. The three sites enabled us to get data over very diverse environments: oligotrophic and sediment-poor waters in the southwest lagoon of New Caledonia, eutrophic waters in the Cienfuegos Bay (Cuba), and sediment-rich waters in the Laucala Bay (Fiji). In this paper, optical algorithms for turbidity are presented per site based on 113 stations in New Caledonia, 24 stations in Cuba and 56 stations in Fiji. Empirical algorithms are tested at satellite wavebands useful to coastal applications. Global algorithms are also derived for the merged data set (193 stations). The performances of global and local regression algorithms are compared. The best one-band algorithms on all the measurements are obtained at 681 nm using either a polynomial or a power model. The best two-band algorithms are obtained with R412/R620, R443/R670 and R510/R681. Two three-band algorithms based on Rrs620.Rrs681/Rrs412 and Rrs620.Rrs681/Rrs510 also give fair regression statistics. Finally, we propose a global algorithm based on one or three bands: turbidity is first calculated from Rrs681 and then, if < 1 FTU, it is recalculated using an algorithm based on Rrs620.Rrs681/Rrs412. On our data set, this algorithm is suitable for the 0.2-25 FTU turbidity range and for the three sites sampled (mean bias: 3.6 %, rms: 35%, mean quadratic error: 1.4 FTU). This shows that defining global empirical turbidity algorithms in tropical coastal waters is at reach.

Above-water measurements of reflectance and chlorophyll-a algorithms in the Gulf of Lions, NW Mediterranean Sea.

Above-water reflectance and surface chlorophyll a concentrations (Chl a) were measured in the Gulf of Lions, northwestern Mediterranean Sea in 2000 and 2001 in order to test Chl a inversion algorithms. Surface waters were separated in Case 2 waters in the Rhône River plume and proximal Region of Freshwater Influence (ROFI) stations, and Case 1 waters at all the other stations. Case 2 waters were characterized by R443/R555 < R443/R510 < R490/R555 < R490/R510 < 1. In the first part, we compared the concurrent reflectance measurements made with a scanning polarization radiometer (SIMBAD) and a hyperspectral Ocean Optics radiometer. The comparison of the remote-sensing reflectance (Rrs) values at SIMBAD wavelengths shows excellent agreement for Rrs values higher than 0.01 sr-1. Between the two instruments, reflectance ratios, commonly used in Chl a algorithms, show differences smaller than 2% in the Case 2 waters, and smaller than 20% in the Case 1 waters. In the second part, concurrent measurements of Chl a and of hyperspectral reflectance from 6 cruises were used to analyze the statistical performance of global (OC2, OC4) and regional regression algorithms using mainly SeaWiFS bands. The algorithms were tested first over the entire domain, then separately over the Case 1 and Case 2 waters. Chl a algorithms using band ratios such as the one presented in Bricaud et al. (2002) are suitable for the Case 1 waters. However, taking into account the large dispersion of Chl a for very close reflectance ratios in the Case 2 waters, single band ratios are not suitable for deriving Chl a. The use of a 4-wavelength parameter such as Xc, defined by Tassan (1994), leads to better results in the plume and proximal Rhône ROFI.