the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Synthesis Product for Ocean Time Series (SPOTS) – a ship-based biogeochemical pilot
Nico Lange
Björn Fiedler
Marta Álvarez
Alice Benoit-Cattin
Heather Benway
Pier Luigi Buttigieg
Laurent Coppola
Kim Currie
Susana Flecha
Dana S. Gerlach
Makio Honda
I. Emma Huertas
Siv K. Lauvset
Frank Muller-Karger
Arne Körtzinger
Kevin M. O'Brien
Sólveig R. Ólafsdóttir
Fernando C. Pacheco
Digna Rueda-Roa
Ingunn Skjelvan
Masahide Wakita
Angelicque White
Download
- Final revised paper (published on 16 Apr 2024)
- Supplement to the final revised paper
- Preprint (discussion started on 31 Jul 2023)
- Supplement to the preprint
Interactive discussion
Status: closed
-
RC1: 'Comment on essd-2023-238', Anonymous Referee #1, 30 Aug 2023
Review of “Synthesis Product for Ocean Time-Series (SPOTS) – A ship-based biogeochemical pilot” by Nico Lange et al.
This paper provides a synthesis product for ocean biogeochemical time-series data from 12 ship-based time-series programs. Currently, we are still lacking a synthesis time-series dataset. This paper carried out a lot of data processing and the product generated by this work is valuable. More importantly, this product complements existing global ocean carbon data and is anticipated to drive increased recognition and utility of ship-based biogeochemical time-series data. Moreover, the study's recommendations for metadata reporting, quality control measures, and uncertainty estimation aim to promote more standardized and uniform practices within the scientific community, contributing to ongoing advancements in ocean research and data accessibility. I only have a few comments on this work. My major and minor are as below.
Major comments:
There are still some other time series measurements of carbonate chemistry stations such as Hawaii Ocean Time Series (HOT), Bermuda Atlantic Time Series Study (BATS), European Station for Time Series in the Ocean Canary Island, (ESTOC) and South East Asia Time Series study (SEATS). Since this paper provides carbonate chemistry data, why not incorporate the data from these stations into the study? The authors may need to consider adding these time-series dadatasetso this paper.
Table 3. The desired part is not clear to me. Authors may need to consider explaining more details about the method used for each parameter and providing appropriate references for those methods.
Minor comments
Line 135 Table 1: NO3- indicates dissolved inorganic nitrate and NO2 indicates dissolved nitrite. Please be consistent in the text.
Line 225 to 226: The Table 1 shows that K2 and KNOT started in 1999 and 1997, respectively, which is not consistent with the information in your main text Line 225 to 226.
Line 300: Typo in Lauvest and Tanhu, 2015))
Line 355 Table 3. In the variable column, is parameter referring to POC, PON, POP?
Also, carbon and nitrate. Isn’t it carbon and nitrogen? Should clarify the parameters related to particulate matter.
Line 365-385: What about salinity, particulate carbon and particulate nitrogen?
Line 372-373: Please provide more specific information about which carbon variable you are referring to DIC, DOC or POC?
Line 580-581: Why not use the estimated partial pressure of CO2 data from DIC, TA and pH measurements?
Line 808: Typo in de Boyer Montégut et al. (2004))
What does NA mean on table 4 and 5?
Citation: https://doi.org/10.5194/essd-2023-238-RC1 -
AC1: 'Reply on RC1', Nico Lange, 03 Nov 2023
Thank you very much for your feedback and the recommendation for publication after minor corrections. We updated the manuscript according to your specific comments and added extra information that are linked to your major comment regarding Table 3. Furthermore, note that the ODIS Section has been improved.
See below for further detailed replies:
- Missing time-series stations
The HOT data is included in the pilot, however the other station data could not be included due to different reasons. Each of these time-series was contacted and asked for participation in the pilot though.
ESTOC: The ship-based (distinct) data was not ready in time preventing a possible inclusion of these very valuable data in this pilot product. Very recently ESTOC was able to publish parts of their bottle data (carbonate chemistry, with temperature and salinity; https://www.frontiersin.org/articles/10.3389/fmars.2023.1236214/full) covering a time from 1994 until 2023. However, there is still work ongoing regarding their nutrient measurements. Contact with the PIs of ESTOC has been established since the beginning of this project and is still ongoing so that for the planned update of SPOTS, the inclusion of ESTOC`s full bottle datasets are a priority.
BATS: At the early stages of this effort (concept phase), the scientific team behind BATS have informed us that they did not want to be part of this pilot. Despite several discussions we could not convince them to participate. However, we hope that through the publication of this pilot we can convince them of the importance of SPOTS and include BATS in the future.
SEATS: We contacted SEATS together with the other included time-series stations. Unfortunately, we did not receive a reply.
Despite these (and other) time-series data missing in this pilot of SPOTS, we believe that the pilot represents the fundament for a larger, living, and sustained SPOTS in the near future. We further believe that the above listed stations will be included with the upcoming updates.
- Table 3
We added a supplementary table, S9, which links the used references for defining required and desired method recommendations to each variable. Note that in some instances we also relied upon internal discussions and comparisons that are not published. Such instances are indicated accordingly, too. We additionally added an example in Section 3.3.2 (using total alkalinity) to explain the meaning of Table 3, the mentioned methods, and the corresponding flag assignment, in more detail.
Further, the parameters “Particulate Matter” refer to TPN, TPP, and TPC – the manuscript (and product) has been updated accordingly.
- Why not use the estimated partial pressure of CO2 data from DIC, TA and pH measurements?
We have decided to only include measured variables in the product. There are several reasons for this decision, one of them being a clear link between original data and the synthesis product. However, in particular regarding calculating missing inorganic carbon parameters, recent studies indicate that measured and calculated variables can differ significantly, and the estimations also depend strongly on the variables used (e.g., using DIC and TA to calculate pCO2 vs using DIC and pH to calculate pCO2). Accordingly including such calculated variables might harm/bias trend analyses, especially if combining measured data from different sources, in terms of analytical and quality control procedures.
Two review works are being prepared for publication from the OCSIF OCB working group, https://www.us-ocb.org/ocean-carbonate-system-intercomparison-forum/, dealing with the current state and recommendations for CO2 system measurements, the caveats and future work on internal consistency studies and data analysis and reporting. Particular recommendations and warnings about combining measured and calculated CO2 data are given. The main message is that a careful evaluation of the uncertainties in both measured and calculated data should be given to finally commit any trend analysis.
We have chosen to suggest methods on how to increase the data coverage through computations in Section 6.2 (e.g., Lines 861 – 867).
- What about uncertainty estimates for salinity, particulate carbon and particulate nitrogen?
Precision estimates of salinity, TPC, TPN and TPP that are based on replicate measurements are included in the product. Note that for particulate matter neither CVOO nor CARIACO provide such measurements. The latter provides precision values of the instruments which would not provide a fair comparison to precision estimates of replicate samples though. Similarly, all existing accuracy measurements of particulate carbon and nitrogen are either based on instruments (CARIACO) or ratios (TPC/TPN; ALOHA) and were not included in the product as only accuracies based on measurements against reference materials were included. Salinity accuracy estimates were not provided.
Citation: https://doi.org/10.5194/essd-2023-238-AC1
-
AC1: 'Reply on RC1', Nico Lange, 03 Nov 2023
-
RC2: 'Comment on essd-2023-238', Anonymous Referee #2, 23 Sep 2023
Lange et al. present a well-written manuscript synthesizing biogeochemistry data of existing ocean time-series stations on a global scale. The authors have done a great job presenting the data flow and demonstrating the high quality of the compiled dataset. The manuscript is clear and nearly free of typos. The final data product will help detect changes in marine ecosystems globally and greatly benefit the entire oceanography community. I do have some minor comments and suggestions. Overall, I recommend publication with minor revisions.
%--------------
Major interpretation points
First, the authors didn’t include some other easily accessible time-series stations in the current compilation, such as BATS, OSP, and ESTOC. The data product would be more complete if the authors could add those as well. See Figure 2 in Benway et al. (2019) (https://doi.org/10.3389/fmars.2019.00393) for the locations of these stations.
Secondly, I think one aspect of future perspectives that could be emphasized more. Other BGC parameters are also routinely measured at some of the time series stations. For example, parameters such as pigments, dissolved organic nitrogen (DON), and sediment trap particles are measured at Station ALOHA. Likewise, DON, dissolved organic phosphorus, and sediment trap particles are measured at CARIACO.
%----------------
Interpretation points, by line #
Lines 135–136: Which time-series station did you use CTD rather than bottle oxygen data?
Line 272- Figure 2: Why is applying QC tests optional? Based on 3.3.1, QCs from the time-series stations are applied if they are available. To me, it seems like the authors are applying this step to all stations before the BP checks. Therefore, it's not an optional step, and the authors should consider removing "(optional)" here.
Lines 276–277: Are data from Munida and RADCOR directly obtained from the responsible PIs? If so, the authors should either state it here or in Table 1 or S1. In some cases, the data downloaded from data centers do not cover all years mentioned in Table 1. Did the authors get those missing data directly from PIs? Please clarify.
Lines 304–305: If the authors are converting the carbonate system parameters (e.g., pH) that are sensitive to temperature, it makes sense to use reported laboratory temperature for the conversion. However, for nutrients and DOC, it makes more sense to use in-situ T and S measurements for the conversion.
Lines 311–312: The conversion of POP from ug/kg to umol/kg should also be mentioned here.
Lines 313–314: If the authors look at sediment trap data at Station ALOHA which have measurements of both total PC and PIC flux, one can calculate the fraction of PIC in sinking particles. The mean fraction of all data is 10.4±4.5% using all the data. If we make the assumption that suspended particles collected by bottles have the same PIC% as sinking particles, we can easily see that the authors would overestimate the POC concentrations at Station ALOHA by ~11%. This is not a fraction that can be ignored, even at the subtropical Station ALOHA. The authors should consider labeling these data (e.g., PC as POC) with different quality flags to make users understand this caveat or incorporate this in the uncertainty analysis.
Lines 450–455: Do the authors use samples from all depths or just below 1500 db? It's a bit unclear what the authors meant by relaxing the minimum depth requirement.
Line 624: Does this high fraction also represent the natural variabilities of DOC at Station ALOHA? If it's not circulation-driven (since low variability in T and S), is it possible to be caused by biological activities? The same question can be applied to all parameters that have a higher variability than the accuracy of the method.
Lines 634–635: It would be helpful if the authors could list the names of all these stations so that we can determine where to see high natural variabilities.
Line 660: What about 33KB20020923 or 325019971101 or 49HG20010813? They seem to fall within the criteria, at least for the 2˚ spatial box, unless the minimum of 1500 m criterion is applied.
Line 955: I may have missed it. Is this point about O2 and pH discussed at all in the manuscript? If not, please add more detail.
%-----------------
Editorial remarks (by line #):
Line 66: Please double-check this doi number. I cannot find this dataset using this doi.
Table 1- Variables: Why not include temperature as one of the measured properties? Some of these stations have sediment trap measurements, such as Station ALOHA, CARIACO, and BATS.
Table 1- KNOT: Data downloaded from this link is from 1992 to 2008, not 1997 to 2020. Plus, there are no direct measurements of pH reported.
Table 1-K2: Data are only from 1999 to 2008. Neither pH nor DOC values are reported.
Table 1-ALOHA: Data are from 1988 to 2016 rather than 1988 to 2019. Plus, one can download data from Station ALOHA from 1988 to the end of 2022 for the bottle data. Why do the authors only report time ranges between 1988 and 2019?
Table 1-Munida: Are data from the Munida and RADCOR not available?
Table 1-GIFT: There aren't any PO4 and DOC data available associated with the DOI.
Table 1-CVOO: The dataset is currently in review and the DOI link doesn't work now.
Table 1-Irminger Sea: The current time range is from 1983 to 2019, but there are also data available from 2020 to 2022.
Table 1-Iceland Sea: Likewise, there are also measurements from 2020 to 2022. Additionally, the earliest measurements were taken in 1985 rather than 1983.
Line 218: Replace “IS-TS” with “IC-TS”.
Line 290: Replace “fill” with “filling”.
Line 298: The citation should be Jiang et al. (2022) rather than Liqing et al. (2022). Change it throughout the manuscript.
Line 302: In Table S4, particulate matter has the unit of ug/kg in the tab "Product" and umol/kg in the tab "Time-series Stations". Be consistent and change all units to umol/kg.
Line 356- Table 3: Replace “All” with “All except silicate”.
Line 480: Replace "E.g." with "For example".
Line 565: It would be helpful to have figures like Figure 4 or 5 for all parameters in the supplemental, not just for nitrate and TA.
Lines 638 & 653: Use “Hash” rather than “Rhombus”.
Line 653 & 859: Use “Asterisk” rather than “Asterix”.
Line 654: Is it just CTD salinity or CTD oxygen? If it’s just CTD salinity, please be more specific.
Tables 4&5: Add a column of temperature.
Line 654: The authors should consider reporting the mean offset consistently either as an absolute number or a relative fraction, rather than a combination, in Table 5.
Line 654-IcelandSea: Why is there only one number in the parathesis? If it's either 3:0 or 0:3, how can the authors get 1% for the comparison?
Line 705: Suggest separating the dashed line and the text more in Figure 6.
Line 855: Consider using black-outlined empty circles for the monthly means so that observations can be better seen.
Line 860: If the fit is not significant, I don't think the authors need to put a red line fit. It will make it easier for the readers to see which relationships are significant.
Line 891: The data product is not available online yet. The authors need to post the data product before this manuscript is accepted.
Line 906: This link doesn't work. Is it supposed to be "https://data.pmel.noaa.gov/generic/erddap/tabledap/spots_bgc_ts.html"?
Line 915: This link doesn't work. Is it supposed to be "https://data.pmel.noaa.gov/generic/erddap/tabledap/spots_bgc_ts.graph"?
Citation: https://doi.org/10.5194/essd-2023-238-RC2 -
AC2: 'Reply on RC2', Nico Lange, 03 Nov 2023
Thank you very much for your feedback and the recommendation for publication after minor corrections. We updated the manuscript according to your specific comments and added extra information that are linked to your major comments. Furthermore, note that the ODIS Section has been improved, and that the DOI is now working.
See below for further detailed replies:
- "Missing” time-series stations
ESTOC: The ship-based (distinct) data was not ready in time preventing a possible inclusion of these very valuable data in this pilot product. Very recently ESTOC was able to publish parts of their bottle data (carbonate chemistry, with temperature and salinity; https://www.frontiersin.org/articles/10.3389/fmars.2023.1236214/full) covering a time from 1994 until 2023. However, there is still work ongoing regarding their nutrient measurements. Contact with the PIs of ESTOC has been established since the beginning of this project and is still ongoing so that for the planned update of SPOTS, the inclusion of ESTOC`s full bottle datasets are a priority.
BATS: At the early stages of this effort (concept phase), the scientific team behind BATS have informed us that they did not want to be part of this pilot. Despite several discussions we could not convince them to participate. However, we hope that through the publication of this pilot we can convince them of the importance of SPOTS and include BATS in the future.
OSP: OSP is clearly a very important fixed time-series station with a long record including biogeochemical data. However, we chose not to include OSP in the pilot of SPOTS as all DFO Line-P data was recently embedded in GLODAP. Having said that OSP data will be a priority for an update of SPOTS which will become even more important once mooring data are included in SPOTS as well.
Despite these (and other) time-series data missing in this pilot of SPOTS, we believe that the pilot represents the fundament for a larger, living, and sustained SPOTS in the near future. We further believe that the above listed stations will be included with the upcoming updates, as hopefully accompanied by many others.
- Stronger emphasis on other routinely measured BGC variables
We added an extra sentence in the outlook section, now also extending the SPOTS vision to BGC non-EOVs. However, for the pilot we are convinced that focusing on BGC EOVs is appropriate following the focus set by the Framework of Ocean Observations.
- “Optional QC” step in workflow
This step relates to “external” QC, i.e. (additional) QC applied by the SPOTS team (in collaboration with the time-series program). Only Munida and CVOO have chosen to undergo this external check, all other flags have been assigned by the time-series programs themselves, i.e., “internally”. The figure is updated accordingly to highlight that external QC is meant here.
- Data obtained directly from PI
The data sources are provided in S1. The last column (G) shows “PI” if the complete and merged data is obtained directly from the PI (e.g., Munida). The entry “Both” indicates that some data are already archived (i.e., a DOI is present) and that missing/additional data was obtained directly from the PI. Note that for some time-series programs the data was obtained from multiple published data sources, e.g., IrmingerSea, KNOT, RADCOR. For those instances all corresponding DOIs are provided in column F in S1. We added an extra sentence to make this clearer.
- Table 1-ALOHA: Data are from 1988 to 2016 rather than 1988 to 2019. Plus, one can download data from Station ALOHA from 1988 to the end of 2022 for the bottle data. Why do the authors only report time ranges between 1988 and 2019?
The data we used are the most recent one with a DOI (from BCO-DMO, https://www.bco-dmo.org/dataset/3773#data-files, see Suppl. Table 1).
- Table 1-Munida: Are data from the Munida and RADCOR not available?
‘NA’ only denotes that no DOI is available for these datasets.
Particularly for the RADIALES A Coruña (RADCOR) times series, the complete and merged data set including physical, biogeochemical and biological variables is available on request to the web site (https://www.seriestemporales-ieo.net/) or the PI (formerly Antonio Bode, antonio.bode@ieo.csic.es; now Mar Nieto-Cid, mar.nietocid@ieo.csic.es). Particularly, CO2 and ancillary data have been partially discussed (Guallart et al., 2022) and made publicly available (Álvarez et al., 2021). Biological and ancillary data were discussed in Bode et al. (2019 & 2020) with data available in several publications in Pangaea.
- Bode, Antonio; Alvarez, Marta; García, Luz; Louro, M Ángeles; Nieto-Cid, Mar; Ruiz-Villarreal, Manuel; Varela, Marta M (2020): Time series of hydrographic, biogeochemical and plankton variables for a shelf station off A Coruña (NW Spain): 1990-2018. PANGAEA, https://doi.org/10.1594/PANGAEA.919087 -- ------ this one is included in Table 1
- Álvarez, M.; Pérez, F. F.; Fernández Guallart, E.; Castaño, M.; Fajar, N. M.; García-Ibáñez, M. I.; Santiago-Domenech, R.; El Rahman Hassoun (2021). Ten Years of Spanish Ion Carbonate Data in the North Atlantic Ocean and Mediterranean Sea. Digital CSIC. DOI: 10.20350/digitalCSIC/13786 ------ this one is included in Table 1
- Guallart EF, N.M. Fajar, MI. García-Ibáñez, M. Castaño-Carrera, R. Santiago-Domenech, A. El Rahman Hassoun, FF. Pérez, R. Easley, M. Álvarez (2022). Spectrophotometric Measurement of Carbonate Ion in Seawater over a Decade: Dealing with Inconsistencies. Environmental Science & Technology, 56 (12), 7381-7395, https://pubs.acs.org/doi/full/10.1021/acs.est.1c06083.
- Bode, A., M. Álvarez, M. Ruíz-Villarreal and M.M. Varela (2019). Changes in phytoplankton production and upwelling intensity off A Coruña (NW Spain) for the last 28 yr. Ocean Dynamics, 69(7), 861-873, https://doi.org/10.1007/s10236-019-01278-y.
- Bode, A.; Álvarez, M.; García García, L.M.; Louro, M.Á.; Nieto-Cid, M.; Ruíz-Villarreal, M.; Varela, M.M. (2020). Climate and Local Hydrography Underlie Recent Regime Shifts in Plankton Communities off Galicia (NW Spain). Oceans, 1, 181-197. https://doi.org/10.3390/oceans1040014.
- Bode, Antonio; Louro, M Ángeles; Rey, Elena (2020): Monthly data of microzooplankton abundance for a shelf station off A Coruña (NW Spain): 1990-2015. PANGAEA, https://doi.org/10.1594/PANGAEA.919081,
- Bode, Antonio; Louro, M Ángeles; Rey, Elena (2020): Monthly data of mesozooplankton abundance and biomass for a shelf station off A Coruña (NW Spain): 1990-2018. PANGAEA, https://doi.org/10.1594/PANGAEA.919080,
- Bode, Antonio (2019): Monthly series of phytoplankton and zooplankton abundance as well as biomass collected between 1989 and 2016 at A Coruña (NW Spain). PANGAEA, https://doi.org/10.1594/PANGAEA.908815,
- Table 1-Irminger Sea/Iceland Sea: The current time range is from 1983 to 2019, but there are also data available from 2020 to 2022.
The additional data (2019-2022) were archived / provided too late for their inclusion into the pilot (see ‘Lineage’ at metadata landing page).
- Unit conversion using lab-temperatures
We followed the recommendations by Jiang et al, (2022), who explicitly state that for nutrients that the lab-temperature should be used. The reason behind is the following: given the methodology for the nutrient analysis, usually segmented flow analysis, all the solutions (samples and reagents) are assumed to be stable and at the lab temperature that should be constant. TPP was given in nM (CARICAO, see S4) and thus there was no need to use/mention the inverse standard atomic masses of phosphorus.
- Minimum variability and biological activity
Generally, if circulation effects can be excluded through low salinity (and temperature) variabilities, it is reasonable to link relatively high coefficients of variation (in comparison to accuracy) estimates to biological activity. However, to minimize the effects of biological activity we have chosen to estimate the coefficients of variabilities on the most invariable oxygen layers (Lines 415 – 420). Given that for, e.g., ALOHA the coefficient of variation for oxygen is 0.7%, we find it unlikely that the high variabilities of DOC are mainly dominated by biological effects, but rather represent measurement uncertainty not fully captured by the usage of RMs or blank control. The same argument holds for other relatively high coefficients of variations. Also note, that ALOHA accuracy estimates are only provided for four years (2008 – 2011). Nevertheless, even if oxygen variations indicate an established equilibrium, we cannot fully exclude that biological activity can have effects on the calculated variabilities of the other BGC measurements. In the revised version, we have highlighted that users must cautiously interpret these results.
- “Missing” Crossover Cruises
Given the large pool of deep data at ALOHA, we have chosen to set the minimum depth to 2000 dbar (Table 5, column “layer”), but still at least two samples from the same profile are needed below that depths to perform the crossover analysis. For the mentioned and not considered cruises the minimum depths would actually have to be below1000 dbar.
- It would be helpful to have figures like Figure 4 or 5 for all parameters in the supplemental, not just for nitrate and TA.
A set of figures (one for each time-series) has now been attached in the supplemental.
- Table 5: Add a column of temperature
GLODAP has not (yet) implemented temperature offset nor does it provide consistency estimates for temperature. We have accordingly only added temperature to Table 4 and not to Table 5.
- Providing the mean offset consistently either as an absolute number or a relative fraction, rather than a combination
The offsets are given following the example set by GLODAP, which uses a combination of both types depending on the range of expected values.
- Change layout / design of figures 6 and 7
Both figures are generated by existing scripts and tools, accordingly we have not updated their layout.
Citation: https://doi.org/10.5194/essd-2023-238-AC2
-
AC2: 'Reply on RC2', Nico Lange, 03 Nov 2023