High-spatiotemporal reconstruction of biogeochemical dynamics in Australia integrating satellites products and in-situ observations (2000&ndash;2022)

Zhang, Xiaohan; Wang, Lizhe; Yan, Jining; Wang, Sheng

doi:https://doi.org/10.5194/essd-2024-219

Preprints

https://doi.org/10.5194/essd-2024-219

Preprints

02 Jul 2024

| 02 Jul 2024

Status: this preprint was under review for the journal ESSD but the revision was not accepted.

High-spatiotemporal reconstruction of biogeochemical dynamics in Australia integrating satellites products and in-situ observations (2000–2022)

Xiaohan Zhang, Lizhe Wang, Jining Yan, and Sheng Wang

Abstract. The marine biogeochemical time-series products, which include total alkalinity, inorganic carbon, nitrate, phosphate, silicate, and pH, constitute a foundational support mechanism for the ongoing surveillance of oceanic biogeochemical changes. These products play a critical role in facilitating research focused on dynamic monitoring of marine ecosystems and fostering sustainable oceanic development. However, existing monitoring methodologies are hampered by inherent limitations, notably the paucity of observational products that simultaneously offer high spatial and temporal resolutions. Furthermore, the interpolation methods typically employed in these contexts frequently prove low-effective on a large scale, resulting in data with extensive temporal and spatial expanses that are difficulty for applications aimed at monitoring large-scale ocean dynamics. A novel integration of the CANYON-B and Random Forest regression methods was explored to address these challenges in reconstructing key marine biogeochemical parameters. This work reconstructs the concentrations of these marine biogeochemicals at the sea surface within Australia's Exclusive Economic Zone over the period from 2000 to 2022 on a 1-kilometre scale. The approach involves the amalgamation of multi-source in-situ ocean chemistry time-series observations with MODIS Terra ocean reflectance imagery and ocean water colour product distributions. This research highlights the substantial capabilities of machine learning for the large-scale reconstruction of ocean chemistry data, introducing a new, viable method for utilising in-situ measurements and optical imagery in reconstructing marine biogeochemical elements, thereby significantly enhancing our ability to monitor large-scale ocean dynamics. The datasets generated and analysed in this study are available on Science Data Bank (https://doi.org/10.57760/sciencedb.09331) (Zhang et al., 2024)

Received: 04 Jun 2024 – Discussion started: 02 Jul 2024

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this preprint. The responsibility to include appropriate place names lies with the authors.

Download & links

Xiaohan Zhang, Lizhe Wang, Jining Yan, and Sheng Wang

Status: closed

RC1:
'Comment on essd-2024-219', Henry Bittig, 20 Nov 2024

In this work, the authors present a methodology to densify sparse observations of ocean biogeochemistry into a high resolution dataset with the help of machine learning methods. Their starting point is the CANYON-B neural networks, based on GLODAPv2 ocean chemistry sparse sampling, as well as high resolution satellite remote sensing, which are combined by machine learning methods (random forest) to build a high resolution data cube of ocean chemistry in coastal waters around Australia to help ecosystem monitoring and management.

The manuscript is clearly structured, easy to read and well-organized.
The main criticism I have is that the authors want to use open-ocean based parameterizations (CANYON-B) and apply them to highly dynamic coastal systems (Australian EEZ waters), where processes and dynamics do not match together.
Taking the author’s references as indication, the authors have a strong background in remote sensing, less so in sea going biogeochemistry and sample analysis:

The GLODAPv2 data set is a collection of past and recent basin scale, coast-to-coast repeat hydrography cruises, done through different programs such as WOCE or more recently GO-SHIP. The focus is on the deep open ocean. The GO-SHIP target is to cover the world’s ocean with a decadal repeat frequency, i.e., every repeat hydrography section should be covered at least once per 10 years. There are a handful of sections with “higher frequency” (meaning every 1 or 2 years), too, but they are mostly located in the Atlantic basin. In addition, research cruises have a tendency to occur in favourable weather conditions, and in fact there is a seasonal bias of the GLODAPv2 cruises towards summer rather than winter timing. To summarize, the GLODAPv2 dataset is based on multi-decadal observations, with multi-year time resolution between repeats at most. There is no unbiased seasonal sampling yet alone a seasonal resolution. In consequence, with CANYON-B being based entirely on the GLODAPv2 data (and the information captured within these data), CANYON-B’s relations are based on processes acting in the open ocean on multi-annual time scales but neither on seasonal or sub-seasonal processes nor on coastal gradients or complex dynamics.
Like the authors state in lines 40-47:

"Most large-scale ocean chemistry datasets are derived from infrequent ship-based surveys or fixed-point observatories, which are then interpolated to create continuous spatial fields. This interpolation, while necessary, introduces substantial uncertainties, particularly in dynamic regions where biogeochemical properties can vary significantly over short distances and time periods. Traditional interpolation methods [...] may not adequately capture complex gradients or the temporal variability of ocean processes [...]. Such shortcomings can lead to misleading representations of marine biogeochemical environments, potentially skewing our understanding of oceanic processes and their responses to environmental changes."

the method used must be fit for the application.
Here, CANYON-B is not fit for seasonal or sub-seasonal application, nor is it fit for coastal waters, i.e., one of the assumptions of this work is invalid, unfortunately. Even if CANYON-B is a machine learning-based method and not a “traditional interpolation method” like the ones listed by the authors, it is a data-driven method that inherits the limits of its parent training dataset. If some information on different orders of magnitude (notably: seasonal dynamics; near-shore coastal processes) is not captured at all by the training data, even a fancy machine learning approach cannot infer that information. In addition, CANYON-B’s predictive skill is high on interior ocean biogeochemistry. It degrades towards the surface and for surface applications, where there is much stronger variability and where the tight biogeochemical coupling between oxygen cycling and CO₂ cycling breaks down in waters in contact with the atmosphere (due to different time scales of air sea gas exchange). (I.e., oxygen becomes a less adequate predictor of the chemical species of interest in surface waters than in the ocean interior.)
I have therefore great concerns towards the validity of the presented data set. It may show some interesting dynamics, as those data interpolation methods can always provide you with an output number, but I wouldn’t assume them at all to be trustworthy or reliable.
The remainder of their approach, to use in-situ Argo or glider data (temperature, salinity, oxygen, chla, turbidity/POC) and to combine them with ocean colour remote sensing by random forest regression to inform about the biogeochemical conditions (of temperature, salinity, oxygen, chla, turbidity/POC) in coastal waters seems valid. Here, the scale of observations (weekly to monthly timescale, some km resolution) matches the scope of the data product aimed at. But as outlined above, the transfer of Argo/glider data to CANYON-B outputs of nutrients or carbon system parameters is an invalid step in the manuscript's setting in the Australian EEZ waters. The monthly trends of chemical variables cannot be considered as robust or reliable.
Only the validation against independent observations (section 4.2.2) would be somewhat suited to hint towards how strongly the CANYON-B derived quantities match with field observations. I would have expected a style like Figure 10 for the Figure 12 comparison, in which is it hard to make out details. In addition, a presentation of data grouped by month would have been helpful to assess the invalidity of the monthly time scale target (which seems to be discernible from Figure 13, where there are repeated (annual?) oscillations in the percentage error). What I can make out from Figure 12/13 and the provided statistics is that the CANYON-B-based products are able to hit the ballpark value for the actually observed carbon parameters and nutrients (as – at some point – the EEZ waters connect to the open ocean), but that the products fail to get the variability and dynamics.
A way to make their approach work would be to use the independent observations of the IMOS National Reference Stations and establish a regionalized/local “CANYON-B”-like parameterization to connect temperature/salinity/oxygen/... with the target parameters of interest (carbon parameters and nutrients) that is applicable to the system and timescales the authors work in.
Further comment: It remains unclear whether this is a surface data product or one that extends into the water column. I believe it captures the surface only, but this should be stated somewhere (e.g., title, abstract).

Citation: https://doi.org/10.5194/essd-2024-219-RC1
- AC1: 'Reply on RC1', xiaohan zhang, 03 Dec 2024
  
  Response to reviewers' comments:
  We are very grateful to the reviewers for their comments on our manuscript. Firstly, it should be clarified that the product produced by our work is indeed a ‘Marine surface product’, as we state in the abstract of the revised manuscript. All added text has been bolded in red, and deleted text has been bolded in blue with strikethrough.
  
  Review comments 1: ‘the CANYON-B output for transferring Argo/glider data to nutrient or carbon system parameters is unreliable’.
  Author's Response:
  We thoroughly investigated this issue before proceeding with our work. Our work is mainly based on the following work:
  The first work is ‘Bittig, Henry C., et al. An alternative to static climatologies: Robust estimation of open ocean CO2 variables and nutrient concentrations from T, S, and O2 data using Bayesian neural networks.’ Frontiers in Marine Science 5 (2018): 328.’ This is the article in which the CANYON-B model is presented, which states that ‘The quantitative analysis of errors between the CANYON-B model simulations and Argo observational values shows that the mean difference for the CANYON-B with the transfer data set is -1 ± 16 µatm.’
  The second work is ‘Mignot, Alexandre, et al. ’Using machine learning and Biogeochemical-Argo (BGC-Argo) floats to assess biogeochemical models and optimise observing system design.’ Biogeosciences 20.7 (2023): 1405-1422.’ This work uses data from Argo floats as CANYON-B inputs, and uses the estimated data from the output CANYON-B outputs and the Argo observations to construct a comprehensive dataset of biogeochemical profiles. The article states that ‘CANYON-B estimates of NO₃and pH were merged with measurements based on the combination of the RMS error of the CANYON-B estimates (NO₃ = 0.7 µ mol kg -1 and pH = 0.013) and the BGC-Argo observation error (NO₃ = 0.5 µ mol kg -1and pH = 0.07) are of the same order of magnitude. ’
  The third work is ‘Asselot, Rémy, et al. “Anthropogenic carbon pathways towards the North Atlantic interior revealed by Argo-O₂, neural networks and back-calculations.” Nature Communications 15.2024: 1630.’ The article uses AT and DIC calculations using CANYON-B to reveal anthropogenic carbon pathways towards the North Atlantic interior.
  The fourth work is ‘Sauzède, Raphaëlle, and Hervé Claustre. ’Ocean Nutrient profiles vertical distribution Product MULTIOBS_GLO_BIO_REP_015_006. (2019).’ The publisher of this document is the European Centre for the Application of Marine and Meteorological Satellites (CMEMS, Copernicus Marine Environment Monitoring Service). This work utilises pressure, temperature, salinity and oxygen concentrations obtained from BGC-Argo data and estimated by the CANYON-B neural network method for nitrate, phosphate and silicate. The accuracy was evaluated as nitrate (NO₃): root mean square error (RMSD) of 0.68 µmol/kg; phosphate (PO₄) RMSD of 0.051 µmol/kg; and silicate (Si(OH)₄): RMSD of 2.3 µmol/kg. The deviations of nitrate ranged from -0.55 to +0.57 µmol/kg with an uncertainty of 0.82 to 1.11 µmol/kg. Phosphate deviations ranged from -0.044 to +0.052 µmol/kg with an uncertainty of 0.059 to 0.080 µmol/kg. Silicate deviations ranged from -1.9 to +2.0 µmol/kg with an uncertainty of 2.4 to 3.9 µmol/kg. The article concludes that ‘BGC-Argo data used as CANYON-B inputs for nutrient concentration estimation at multiple sites around the world show good accuracy and low systematic bias.’
  The four existing studies mentioned above proved that Argo data are feasible to be used as input data to CANYON-B, and the estimation results have good accuracy, low systematic bias results, and high reliability of the estimation results. It is based on this that we determined to input Argo data or glider data into CANYON-B to calculate the concentration of the relevant parameters.
  
  Review comments 2: ‘CANYON-B is not suitable for seasonal or sub-seasonal applications.’
  Author's Response:
  We acknowledge that CANYON-B has this limitation. The aim of our work is to improve the temporal and spatial resolution of biogeochemical parameters by combining remotely sensed imagery with estimates from CANYON-B. MODIS remotely sensed imagery has the advantage of being time-series, having a high temporal resolution, and being able to cover a large area.
  However, inversion of biogeochemical parameter concentrations by remote sensing images requires actual data to establish the relationship between image pixel values and concentrations.
  We did obtain measured values through IMOS, and our initial idea was to use these data to construct a random forest model directly with MODIS imagery to invert concentrations. Unfortunately, these measured data were not spatially distributed uniformly, with very little distribution on the northern side of the study area, which would have had some impact on the concentration estimates. In addition, the IMOS measurements are from 2008. As the reviewer stated ‘the accuracy of data-driven methods depends on the parent training data set’, so it is not reasonable to use an inverse model trained on post-2008 data in conjunction with pre-2008 imagery to infer pre-2008 biogeochemical concentrations.
  That is why we decided to use Argo data and glider data, which are uniformly distributed in both time and space, to estimate biogeochemical concentrations via CANYON-B's and then combine them with imagery to complete this work.
  
  Review comments 3: ‘One way to make their approach work would be to use independent observations from IMOS national reference stations and build a regionalised/localised “CANYON-B” type parameterisation’.
  Author's Response:
  Thank you very much for this comment, and it is very unfortunate that, as answered in 2, problems with the timing of the actual data collection and the non-uniform distribution of the data forced us to abandon the idea of constructing a CANYON-B of parameters specific to the study region for that region.
  Of course, there are three main things we did to ensure the quality of the products we produced.
  The first one is the exact date-specific relational correspondence mentioned in 2 to reduce the effect of time differences on the accuracy of the product.
  The second one is the five-fold cross-validation, which makes full use of the existing observation data for model training and validation by dividing the training set and test set several times to comprehensively use the data to evaluate the model's performance on different data subsets, to effectively measure the model's generalisation ability to reduce the risk of overfitting, and to ensure the accuracy of the constructed random forest inversion model.
  Finally, an observation dataset independent of the training data is used to provide a truer picture of the model performance in real-world applications. The use of independent observation data, such as IMOS national reference stations and Argo buoys, effectively avoids spurious accuracy caused by overfitting of the model to the training data, which directly proves the reliability of the product.
  
  Review comments 4: I would have expected a style like Figure 10 for the Figure 12 comparison, in which is it hard to make out details.
  Author's Response:
          We modified Figure 12 to add the same scatterplot as Figure 10. And we modified the textual description of the manuscript regarding Fig. 12.
  Figure 12 modified（Uncompressed images can be viewed in the attachment）:
  Description of the figure 12 (Line 307-316):
          “……This verification includes six components, labelled (a) to (f)(l), each demonstrating the validation of reconstructed ocean parameters against an ensemble of independent observations collected over time.”
          “In Figure 12 (a) and (b), ……. Figure 12 (b)(c) and (d) …….Proceeding to Figure 12 (c)(e) and (f) , …….
  Figures 12 (d), (e), and (f) (g) and (h) , (i) and (j), (k) and (l) …….”
  
  Further comment: It remains unclear whether this is a surface data product or one that extends into the water column. I believe it captures the surface only, but this should be stated somewhere (e.g., title, abstract).
  Author's Response:
  Abstract modified (Line 9):
  “……. This work reconstructs the concentrations of these marine surface biogeochemicals at the sea surface within Australia's Exclusive Economic Zone over the period from 2000 to 2022 on a 1-kilometre scale. …….”
  
  Citation: https://doi.org/10.5194/essd-2024-219-AC1
RC2:
'Comment on essd-2024-219', Anonymous Referee #2, 29 Jan 2025
This work presents a new framework combining the CANYON-B model and the Random Forest Regression approach for reconstructing marine biogeochemical dynamics in Australia's Exclusive Economic Zone (EEZ) at high spatial (1 km) and temporal resolutions. This study fills a gap in traditional interpolation methods for the analysis of dynamic ocean processes at large scales, especially in complex gradient and highly dynamic regions. In addition, the article makes full use of MODIS Terra products with multi-source in situ observations from Argo, gliders and ANMN to achieve accurate modeling across temporal and spatial scales. The innovative approach of combining MODIS satellite remote sensing data with in-situ observations provides an important reference for future ocean monitoring. The manuscript is clearly structured, easy to read and well organized. It is recommended for acceptance with modifications.
Some comments are reported below:
Based on the manuscript's description, it appears that the authors produced a sea surface biogeochemistry data product. It is recommended that the authors explicitly and clearly highlight this in the article to ensure readers fully understand the nature and scope of the data product.

The manuscript describes the use of Argo buoy data as input to the CANYON-B model, followed by validation using the same Argo data. Why not directly use the biogeochemical parameters from Argo for the reconstruction?

In Table 1, the use of the separator “-” in the “Contains Elements” column is potentially confusing, particularly with entries like “Chl-a.” It is suggested to replace the “-” with a comma (”,”) to improve clarity.

Line 145: In Equation 1, it is recommended that the meaning of “Ti(x)” be clearly explained.

Line 385: Some of the references are old (e.g., from 1994 and 2004), and it is suggested to add latest research on random forests.

Why not use IMOS real data directly to train the random forest model?

How to control the date correspondence between point data and MODIS impacts.

The number of citations in the manuscript appears to be relatively low. It is suggested to add more relevant and recent references to provide broader scientific context and support the study’s methodology and findings.
Citation: https://doi.org/10.5194/essd-2024-219-RC2
- AC2: 'Reply on RC2', xiaohan zhang, 30 Jan 2025
  
  Response to reviewers' comments:
  
  Overall review comments: This work presents a new framework combining the CANYON-B model and the Random Forest Regression approach for reconstructing marine biogeochemical dynamics in Australia's Exclusive Economic Zone (EEZ) at high spatial (1 km) and temporal resolutions. This study fills a gap in traditional interpolation methods for the analysis of dynamic ocean processes at large scales, especially in complex gradient and highly dynamic regions. In addition, the article makes full use of MODIS Terra products with multi-source in situ observations from Argo, gliders and ANMN to achieve accurate modeling across temporal and spatial scales. The innovative approach of combining MODIS satellite remote sensing data with in situ observations provides an important reference for future ocean monitoring. The manuscript is clearly structured, easy to read and well organized. It is recommended for acceptance with modifications.
  Author's Response:
  We sincerely appreciate the reviewer’s positive feedback and their recognition of the significance of our work. We are grateful for the acknowledgment of our novel framework, which integrates the CANYON-B model with Random Forest Regression to reconstruct marine biogeochemical dynamics in Australia’s Exclusive Economic Zone (EEZ) at high spatial and temporal resolutions.
  We particularly appreciate the reviewer’s recognition of our effort to bridge gaps in traditional interpolation methods, especially in highly dynamic and complex gradient regions. The integration of MODIS Terra satellite products with multi-source in situ observations from Argo, gliders, and ANMN indeed enhances the accuracy and robustness of our model across different scales. We are pleased that the reviewer considers our approach to provide an important reference for future ocean monitoring.
  In response to the reviewer’s recommendation for modifications, we have carefully revised the manuscript to further improve clarity, structure, and scientific rigor. The specific modifications are detailed in our response to the comments. We hope that the revised manuscript meets the reviewer’s expectations and look forward to their further insights.
  Once again, we greatly appreciate the reviewer’s constructive comments and their recommendation for acceptance.
  
  Review comments 1: Based on the manuscript's description, it appears that the authors produced a sea surface biogeochemistry data product. It is recommended that the authors explicitly and clearly highlight this in the article to ensure readers fully understand the nature and scope of the data product.
  Author's Response:
  Thanks for your suggestion. The product produced by our work is indeed a “marine surface data product”, as we explain in the abstract of the revised manuscript.
  Line 9:
  “……. This work reconstructs the concentrations of these marine surface biogeochemicals at the sea surface within Australia's Exclusive Economic Zone over the period from 2000 to 2022 on a 1-kilometre scale. …….”
  
  Review comments 2: The manuscript describes using Argo buoy data as input to the CANYON-B model, followed by validation using the same Argo data. Why not directly use the biogeochemical parameters from Argo for the reconstruction?
  Author's Response:
  We made this choice for two primary reasons. First, only a small subset of Argo floats is equipped to monitor biogeochemical parameters, and their spatial and temporal distribution is highly heterogeneous within our study area. Directly incorporating these data into the model would introduce bias in model fitting and reduce the accuracy of the reconstructed biogeochemical products.
  Second, the biogeochemical parameters measured by Argo floats serve as an independent dataset to assess the performance of the CANYON-B model. By using Argo-derived pressure, temperature, salinity, and oxygen concentration as inputs, we can evaluate the model's ability to estimate biogeochemical parameters accurately. This step ensures a rigorous validation process and confirms that Argo data can be reliably used as an input for CANYON-B.
  
  Review comments 3: In Table 1, the use of the separator “-” in the “Contains Elements” column is potentially confusing, particularly with entries like “Chl-a.” It is suggested to replace the “-” with a comma (”,”) to improve clarity.
  Author's Response:
  Thanks for your suggestions, which have been revised in the latest manuscript.
  
  Review comments 4: Line 145: In Equation 1, it is recommended that the meaning of “Ti(x)” be clearly explained.
  Author's Response:
  We provide a clear explanation of the meaning of Ti(x) in the revised manuscript.
  Line 155-156：
  “The predictor variables Ti(x) used in the Random Forest regression include column properties (Pres), water color products (Chl-a, POC, and SST), and MODIS Terra satellite-derived ocean reflectance products from bands 1–8.”
  
  Review comments 5: Line 385, Some of the references are old (e.g., from 1994 and 2004), and it is suggested to add the latest research support on random forests in recent years.
  Author's Response:
  Thanks to your suggestion, the latest manuscript has been added with 2 newest studies having to do with random forest regression.
  Line 146：……decision trees and combining their predictions(Wang et al., 2021).
  Line 154-155：……improves the accuracy and stability of the predictions(Grinsztajn et al., 2022; Guillon et al., 2024).
  Review comments 6: Why not use IMOS real data directly to train the random forest model?.
  Author's Response:
  While the IMOS real data is evenly distributed over the study period, there are two main reasons why it was not used as the sole source of training data. First, IMOS provides biogeochemical concentration data only after 2008, whereas our reconstructed product spans the period from 2000 to 2022. Given that Random Forest is a data-driven approach, training the model exclusively on post-2008 data and extrapolating to pre-2008 conditions would be unreliable.
  Second, the spatial distribution of IMOS observations is uneven, with sparse coverage in the northern part of the study area. Relying solely on IMOS data for training could introduce spatial biases and negatively impact the accuracy of the estimated biogeochemical concentrations.
  
  Review comments 7: How to control the date correspondence between point data and MODIS impacts.
  Author's Response:
  When extracting the pixel values from MODIS images corresponding to in situ biogeochemical concentration measurements, we ensured that the image acquisition date matched the observation date. If no MODIS image was available for the exact date, we used data from the nearest available day, either the preceding or following day, as a substitute.
  For example, temperature (T), salinity (S), and oxygen (O₂) concentrations measured on July 10, 2010, were used as inputs to the CANYON-B model to estimate biogeochemical concentrations. These estimates were then matched with pixel values from MODIS image products acquired on the same day at corresponding locations. If no image was available on July 10, pixel values from July 9 or July 11 were used instead.
  Once a sufficient number of correspondences were established and evenly distributed, Random Forest regression was applied to learn the relationship between biogeochemical concentrations and MODIS pixel values. The trained model was then generalized to other times and locations to reconstruct biogeochemical concentrations across the entire study area.
  
  Review comments 8: The number of citations in the manuscript appears to be relatively low. It is suggested to add more relevant and recent references to provide broader scientific context and support the study’s methodology and findings.
  Author's Response:
  Thanks to your suggestion, citations to three recent studies on marine biochemical parameters have been added to the latest manuscript.
  Line 26-27:
  ……protecting marine life(Fourrier et al., 2022).
  Line 43:
  ……vary significantly over short distances and time periods(Mignot et al., 2023).
  Line 49:
  ……high costs associated with data collection(Asselot et al., 2024).
  
  Citation: https://doi.org/10.5194/essd-2024-219-AC2

Status: closed

RC1:
'Comment on essd-2024-219', Henry Bittig, 20 Nov 2024

In this work, the authors present a methodology to densify sparse observations of ocean biogeochemistry into a high resolution dataset with the help of machine learning methods. Their starting point is the CANYON-B neural networks, based on GLODAPv2 ocean chemistry sparse sampling, as well as high resolution satellite remote sensing, which are combined by machine learning methods (random forest) to build a high resolution data cube of ocean chemistry in coastal waters around Australia to help ecosystem monitoring and management.

The manuscript is clearly structured, easy to read and well-organized.
The main criticism I have is that the authors want to use open-ocean based parameterizations (CANYON-B) and apply them to highly dynamic coastal systems (Australian EEZ waters), where processes and dynamics do not match together.
Taking the author’s references as indication, the authors have a strong background in remote sensing, less so in sea going biogeochemistry and sample analysis:

The GLODAPv2 data set is a collection of past and recent basin scale, coast-to-coast repeat hydrography cruises, done through different programs such as WOCE or more recently GO-SHIP. The focus is on the deep open ocean. The GO-SHIP target is to cover the world’s ocean with a decadal repeat frequency, i.e., every repeat hydrography section should be covered at least once per 10 years. There are a handful of sections with “higher frequency” (meaning every 1 or 2 years), too, but they are mostly located in the Atlantic basin. In addition, research cruises have a tendency to occur in favourable weather conditions, and in fact there is a seasonal bias of the GLODAPv2 cruises towards summer rather than winter timing. To summarize, the GLODAPv2 dataset is based on multi-decadal observations, with multi-year time resolution between repeats at most. There is no unbiased seasonal sampling yet alone a seasonal resolution. In consequence, with CANYON-B being based entirely on the GLODAPv2 data (and the information captured within these data), CANYON-B’s relations are based on processes acting in the open ocean on multi-annual time scales but neither on seasonal or sub-seasonal processes nor on coastal gradients or complex dynamics.
Like the authors state in lines 40-47:

"Most large-scale ocean chemistry datasets are derived from infrequent ship-based surveys or fixed-point observatories, which are then interpolated to create continuous spatial fields. This interpolation, while necessary, introduces substantial uncertainties, particularly in dynamic regions where biogeochemical properties can vary significantly over short distances and time periods. Traditional interpolation methods [...] may not adequately capture complex gradients or the temporal variability of ocean processes [...]. Such shortcomings can lead to misleading representations of marine biogeochemical environments, potentially skewing our understanding of oceanic processes and their responses to environmental changes."

the method used must be fit for the application.
Here, CANYON-B is not fit for seasonal or sub-seasonal application, nor is it fit for coastal waters, i.e., one of the assumptions of this work is invalid, unfortunately. Even if CANYON-B is a machine learning-based method and not a “traditional interpolation method” like the ones listed by the authors, it is a data-driven method that inherits the limits of its parent training dataset. If some information on different orders of magnitude (notably: seasonal dynamics; near-shore coastal processes) is not captured at all by the training data, even a fancy machine learning approach cannot infer that information. In addition, CANYON-B’s predictive skill is high on interior ocean biogeochemistry. It degrades towards the surface and for surface applications, where there is much stronger variability and where the tight biogeochemical coupling between oxygen cycling and CO₂ cycling breaks down in waters in contact with the atmosphere (due to different time scales of air sea gas exchange). (I.e., oxygen becomes a less adequate predictor of the chemical species of interest in surface waters than in the ocean interior.)
I have therefore great concerns towards the validity of the presented data set. It may show some interesting dynamics, as those data interpolation methods can always provide you with an output number, but I wouldn’t assume them at all to be trustworthy or reliable.
The remainder of their approach, to use in-situ Argo or glider data (temperature, salinity, oxygen, chla, turbidity/POC) and to combine them with ocean colour remote sensing by random forest regression to inform about the biogeochemical conditions (of temperature, salinity, oxygen, chla, turbidity/POC) in coastal waters seems valid. Here, the scale of observations (weekly to monthly timescale, some km resolution) matches the scope of the data product aimed at. But as outlined above, the transfer of Argo/glider data to CANYON-B outputs of nutrients or carbon system parameters is an invalid step in the manuscript's setting in the Australian EEZ waters. The monthly trends of chemical variables cannot be considered as robust or reliable.
Only the validation against independent observations (section 4.2.2) would be somewhat suited to hint towards how strongly the CANYON-B derived quantities match with field observations. I would have expected a style like Figure 10 for the Figure 12 comparison, in which is it hard to make out details. In addition, a presentation of data grouped by month would have been helpful to assess the invalidity of the monthly time scale target (which seems to be discernible from Figure 13, where there are repeated (annual?) oscillations in the percentage error). What I can make out from Figure 12/13 and the provided statistics is that the CANYON-B-based products are able to hit the ballpark value for the actually observed carbon parameters and nutrients (as – at some point – the EEZ waters connect to the open ocean), but that the products fail to get the variability and dynamics.
A way to make their approach work would be to use the independent observations of the IMOS National Reference Stations and establish a regionalized/local “CANYON-B”-like parameterization to connect temperature/salinity/oxygen/... with the target parameters of interest (carbon parameters and nutrients) that is applicable to the system and timescales the authors work in.
Further comment: It remains unclear whether this is a surface data product or one that extends into the water column. I believe it captures the surface only, but this should be stated somewhere (e.g., title, abstract).

Citation: https://doi.org/10.5194/essd-2024-219-RC1
- AC1: 'Reply on RC1', xiaohan zhang, 03 Dec 2024
  
  Response to reviewers' comments:
  We are very grateful to the reviewers for their comments on our manuscript. Firstly, it should be clarified that the product produced by our work is indeed a ‘Marine surface product’, as we state in the abstract of the revised manuscript. All added text has been bolded in red, and deleted text has been bolded in blue with strikethrough.
  
  Review comments 1: ‘the CANYON-B output for transferring Argo/glider data to nutrient or carbon system parameters is unreliable’.
  Author's Response:
  We thoroughly investigated this issue before proceeding with our work. Our work is mainly based on the following work:
  The first work is ‘Bittig, Henry C., et al. An alternative to static climatologies: Robust estimation of open ocean CO2 variables and nutrient concentrations from T, S, and O2 data using Bayesian neural networks.’ Frontiers in Marine Science 5 (2018): 328.’ This is the article in which the CANYON-B model is presented, which states that ‘The quantitative analysis of errors between the CANYON-B model simulations and Argo observational values shows that the mean difference for the CANYON-B with the transfer data set is -1 ± 16 µatm.’
  The second work is ‘Mignot, Alexandre, et al. ’Using machine learning and Biogeochemical-Argo (BGC-Argo) floats to assess biogeochemical models and optimise observing system design.’ Biogeosciences 20.7 (2023): 1405-1422.’ This work uses data from Argo floats as CANYON-B inputs, and uses the estimated data from the output CANYON-B outputs and the Argo observations to construct a comprehensive dataset of biogeochemical profiles. The article states that ‘CANYON-B estimates of NO₃and pH were merged with measurements based on the combination of the RMS error of the CANYON-B estimates (NO₃ = 0.7 µ mol kg -1 and pH = 0.013) and the BGC-Argo observation error (NO₃ = 0.5 µ mol kg -1and pH = 0.07) are of the same order of magnitude. ’
  The third work is ‘Asselot, Rémy, et al. “Anthropogenic carbon pathways towards the North Atlantic interior revealed by Argo-O₂, neural networks and back-calculations.” Nature Communications 15.2024: 1630.’ The article uses AT and DIC calculations using CANYON-B to reveal anthropogenic carbon pathways towards the North Atlantic interior.
  The fourth work is ‘Sauzède, Raphaëlle, and Hervé Claustre. ’Ocean Nutrient profiles vertical distribution Product MULTIOBS_GLO_BIO_REP_015_006. (2019).’ The publisher of this document is the European Centre for the Application of Marine and Meteorological Satellites (CMEMS, Copernicus Marine Environment Monitoring Service). This work utilises pressure, temperature, salinity and oxygen concentrations obtained from BGC-Argo data and estimated by the CANYON-B neural network method for nitrate, phosphate and silicate. The accuracy was evaluated as nitrate (NO₃): root mean square error (RMSD) of 0.68 µmol/kg; phosphate (PO₄) RMSD of 0.051 µmol/kg; and silicate (Si(OH)₄): RMSD of 2.3 µmol/kg. The deviations of nitrate ranged from -0.55 to +0.57 µmol/kg with an uncertainty of 0.82 to 1.11 µmol/kg. Phosphate deviations ranged from -0.044 to +0.052 µmol/kg with an uncertainty of 0.059 to 0.080 µmol/kg. Silicate deviations ranged from -1.9 to +2.0 µmol/kg with an uncertainty of 2.4 to 3.9 µmol/kg. The article concludes that ‘BGC-Argo data used as CANYON-B inputs for nutrient concentration estimation at multiple sites around the world show good accuracy and low systematic bias.’
  The four existing studies mentioned above proved that Argo data are feasible to be used as input data to CANYON-B, and the estimation results have good accuracy, low systematic bias results, and high reliability of the estimation results. It is based on this that we determined to input Argo data or glider data into CANYON-B to calculate the concentration of the relevant parameters.
  
  Review comments 2: ‘CANYON-B is not suitable for seasonal or sub-seasonal applications.’
  Author's Response:
  We acknowledge that CANYON-B has this limitation. The aim of our work is to improve the temporal and spatial resolution of biogeochemical parameters by combining remotely sensed imagery with estimates from CANYON-B. MODIS remotely sensed imagery has the advantage of being time-series, having a high temporal resolution, and being able to cover a large area.
  However, inversion of biogeochemical parameter concentrations by remote sensing images requires actual data to establish the relationship between image pixel values and concentrations.
  We did obtain measured values through IMOS, and our initial idea was to use these data to construct a random forest model directly with MODIS imagery to invert concentrations. Unfortunately, these measured data were not spatially distributed uniformly, with very little distribution on the northern side of the study area, which would have had some impact on the concentration estimates. In addition, the IMOS measurements are from 2008. As the reviewer stated ‘the accuracy of data-driven methods depends on the parent training data set’, so it is not reasonable to use an inverse model trained on post-2008 data in conjunction with pre-2008 imagery to infer pre-2008 biogeochemical concentrations.
  That is why we decided to use Argo data and glider data, which are uniformly distributed in both time and space, to estimate biogeochemical concentrations via CANYON-B's and then combine them with imagery to complete this work.
  
  Review comments 3: ‘One way to make their approach work would be to use independent observations from IMOS national reference stations and build a regionalised/localised “CANYON-B” type parameterisation’.
  Author's Response:
  Thank you very much for this comment, and it is very unfortunate that, as answered in 2, problems with the timing of the actual data collection and the non-uniform distribution of the data forced us to abandon the idea of constructing a CANYON-B of parameters specific to the study region for that region.
  Of course, there are three main things we did to ensure the quality of the products we produced.
  The first one is the exact date-specific relational correspondence mentioned in 2 to reduce the effect of time differences on the accuracy of the product.
  The second one is the five-fold cross-validation, which makes full use of the existing observation data for model training and validation by dividing the training set and test set several times to comprehensively use the data to evaluate the model's performance on different data subsets, to effectively measure the model's generalisation ability to reduce the risk of overfitting, and to ensure the accuracy of the constructed random forest inversion model.
  Finally, an observation dataset independent of the training data is used to provide a truer picture of the model performance in real-world applications. The use of independent observation data, such as IMOS national reference stations and Argo buoys, effectively avoids spurious accuracy caused by overfitting of the model to the training data, which directly proves the reliability of the product.
  
  Review comments 4: I would have expected a style like Figure 10 for the Figure 12 comparison, in which is it hard to make out details.
  Author's Response:
          We modified Figure 12 to add the same scatterplot as Figure 10. And we modified the textual description of the manuscript regarding Fig. 12.
  Figure 12 modified（Uncompressed images can be viewed in the attachment）:
  Description of the figure 12 (Line 307-316):
          “……This verification includes six components, labelled (a) to (f)(l), each demonstrating the validation of reconstructed ocean parameters against an ensemble of independent observations collected over time.”
          “In Figure 12 (a) and (b), ……. Figure 12 (b)(c) and (d) …….Proceeding to Figure 12 (c)(e) and (f) , …….
  Figures 12 (d), (e), and (f) (g) and (h) , (i) and (j), (k) and (l) …….”
  
  Further comment: It remains unclear whether this is a surface data product or one that extends into the water column. I believe it captures the surface only, but this should be stated somewhere (e.g., title, abstract).
  Author's Response:
  Abstract modified (Line 9):
  “……. This work reconstructs the concentrations of these marine surface biogeochemicals at the sea surface within Australia's Exclusive Economic Zone over the period from 2000 to 2022 on a 1-kilometre scale. …….”
  
  Citation: https://doi.org/10.5194/essd-2024-219-AC1
RC2:
'Comment on essd-2024-219', Anonymous Referee #2, 29 Jan 2025
This work presents a new framework combining the CANYON-B model and the Random Forest Regression approach for reconstructing marine biogeochemical dynamics in Australia's Exclusive Economic Zone (EEZ) at high spatial (1 km) and temporal resolutions. This study fills a gap in traditional interpolation methods for the analysis of dynamic ocean processes at large scales, especially in complex gradient and highly dynamic regions. In addition, the article makes full use of MODIS Terra products with multi-source in situ observations from Argo, gliders and ANMN to achieve accurate modeling across temporal and spatial scales. The innovative approach of combining MODIS satellite remote sensing data with in-situ observations provides an important reference for future ocean monitoring. The manuscript is clearly structured, easy to read and well organized. It is recommended for acceptance with modifications.
Some comments are reported below:
Based on the manuscript's description, it appears that the authors produced a sea surface biogeochemistry data product. It is recommended that the authors explicitly and clearly highlight this in the article to ensure readers fully understand the nature and scope of the data product.

The manuscript describes the use of Argo buoy data as input to the CANYON-B model, followed by validation using the same Argo data. Why not directly use the biogeochemical parameters from Argo for the reconstruction?

In Table 1, the use of the separator “-” in the “Contains Elements” column is potentially confusing, particularly with entries like “Chl-a.” It is suggested to replace the “-” with a comma (”,”) to improve clarity.

Line 145: In Equation 1, it is recommended that the meaning of “Ti(x)” be clearly explained.

Line 385: Some of the references are old (e.g., from 1994 and 2004), and it is suggested to add latest research on random forests.

Why not use IMOS real data directly to train the random forest model?

How to control the date correspondence between point data and MODIS impacts.

The number of citations in the manuscript appears to be relatively low. It is suggested to add more relevant and recent references to provide broader scientific context and support the study’s methodology and findings.
Citation: https://doi.org/10.5194/essd-2024-219-RC2
- AC2: 'Reply on RC2', xiaohan zhang, 30 Jan 2025
  
  Response to reviewers' comments:
  
  Overall review comments: This work presents a new framework combining the CANYON-B model and the Random Forest Regression approach for reconstructing marine biogeochemical dynamics in Australia's Exclusive Economic Zone (EEZ) at high spatial (1 km) and temporal resolutions. This study fills a gap in traditional interpolation methods for the analysis of dynamic ocean processes at large scales, especially in complex gradient and highly dynamic regions. In addition, the article makes full use of MODIS Terra products with multi-source in situ observations from Argo, gliders and ANMN to achieve accurate modeling across temporal and spatial scales. The innovative approach of combining MODIS satellite remote sensing data with in situ observations provides an important reference for future ocean monitoring. The manuscript is clearly structured, easy to read and well organized. It is recommended for acceptance with modifications.
  Author's Response:
  We sincerely appreciate the reviewer’s positive feedback and their recognition of the significance of our work. We are grateful for the acknowledgment of our novel framework, which integrates the CANYON-B model with Random Forest Regression to reconstruct marine biogeochemical dynamics in Australia’s Exclusive Economic Zone (EEZ) at high spatial and temporal resolutions.
  We particularly appreciate the reviewer’s recognition of our effort to bridge gaps in traditional interpolation methods, especially in highly dynamic and complex gradient regions. The integration of MODIS Terra satellite products with multi-source in situ observations from Argo, gliders, and ANMN indeed enhances the accuracy and robustness of our model across different scales. We are pleased that the reviewer considers our approach to provide an important reference for future ocean monitoring.
  In response to the reviewer’s recommendation for modifications, we have carefully revised the manuscript to further improve clarity, structure, and scientific rigor. The specific modifications are detailed in our response to the comments. We hope that the revised manuscript meets the reviewer’s expectations and look forward to their further insights.
  Once again, we greatly appreciate the reviewer’s constructive comments and their recommendation for acceptance.
  
  Review comments 1: Based on the manuscript's description, it appears that the authors produced a sea surface biogeochemistry data product. It is recommended that the authors explicitly and clearly highlight this in the article to ensure readers fully understand the nature and scope of the data product.
  Author's Response:
  Thanks for your suggestion. The product produced by our work is indeed a “marine surface data product”, as we explain in the abstract of the revised manuscript.
  Line 9:
  “……. This work reconstructs the concentrations of these marine surface biogeochemicals at the sea surface within Australia's Exclusive Economic Zone over the period from 2000 to 2022 on a 1-kilometre scale. …….”
  
  Review comments 2: The manuscript describes using Argo buoy data as input to the CANYON-B model, followed by validation using the same Argo data. Why not directly use the biogeochemical parameters from Argo for the reconstruction?
  Author's Response:
  We made this choice for two primary reasons. First, only a small subset of Argo floats is equipped to monitor biogeochemical parameters, and their spatial and temporal distribution is highly heterogeneous within our study area. Directly incorporating these data into the model would introduce bias in model fitting and reduce the accuracy of the reconstructed biogeochemical products.
  Second, the biogeochemical parameters measured by Argo floats serve as an independent dataset to assess the performance of the CANYON-B model. By using Argo-derived pressure, temperature, salinity, and oxygen concentration as inputs, we can evaluate the model's ability to estimate biogeochemical parameters accurately. This step ensures a rigorous validation process and confirms that Argo data can be reliably used as an input for CANYON-B.
  
  Review comments 3: In Table 1, the use of the separator “-” in the “Contains Elements” column is potentially confusing, particularly with entries like “Chl-a.” It is suggested to replace the “-” with a comma (”,”) to improve clarity.
  Author's Response:
  Thanks for your suggestions, which have been revised in the latest manuscript.
  
  Review comments 4: Line 145: In Equation 1, it is recommended that the meaning of “Ti(x)” be clearly explained.
  Author's Response:
  We provide a clear explanation of the meaning of Ti(x) in the revised manuscript.
  Line 155-156：
  “The predictor variables Ti(x) used in the Random Forest regression include column properties (Pres), water color products (Chl-a, POC, and SST), and MODIS Terra satellite-derived ocean reflectance products from bands 1–8.”
  
  Review comments 5: Line 385, Some of the references are old (e.g., from 1994 and 2004), and it is suggested to add the latest research support on random forests in recent years.
  Author's Response:
  Thanks to your suggestion, the latest manuscript has been added with 2 newest studies having to do with random forest regression.
  Line 146：……decision trees and combining their predictions(Wang et al., 2021).
  Line 154-155：……improves the accuracy and stability of the predictions(Grinsztajn et al., 2022; Guillon et al., 2024).
  Review comments 6: Why not use IMOS real data directly to train the random forest model?.
  Author's Response:
  While the IMOS real data is evenly distributed over the study period, there are two main reasons why it was not used as the sole source of training data. First, IMOS provides biogeochemical concentration data only after 2008, whereas our reconstructed product spans the period from 2000 to 2022. Given that Random Forest is a data-driven approach, training the model exclusively on post-2008 data and extrapolating to pre-2008 conditions would be unreliable.
  Second, the spatial distribution of IMOS observations is uneven, with sparse coverage in the northern part of the study area. Relying solely on IMOS data for training could introduce spatial biases and negatively impact the accuracy of the estimated biogeochemical concentrations.
  
  Review comments 7: How to control the date correspondence between point data and MODIS impacts.
  Author's Response:
  When extracting the pixel values from MODIS images corresponding to in situ biogeochemical concentration measurements, we ensured that the image acquisition date matched the observation date. If no MODIS image was available for the exact date, we used data from the nearest available day, either the preceding or following day, as a substitute.
  For example, temperature (T), salinity (S), and oxygen (O₂) concentrations measured on July 10, 2010, were used as inputs to the CANYON-B model to estimate biogeochemical concentrations. These estimates were then matched with pixel values from MODIS image products acquired on the same day at corresponding locations. If no image was available on July 10, pixel values from July 9 or July 11 were used instead.
  Once a sufficient number of correspondences were established and evenly distributed, Random Forest regression was applied to learn the relationship between biogeochemical concentrations and MODIS pixel values. The trained model was then generalized to other times and locations to reconstruct biogeochemical concentrations across the entire study area.
  
  Review comments 8: The number of citations in the manuscript appears to be relatively low. It is suggested to add more relevant and recent references to provide broader scientific context and support the study’s methodology and findings.
  Author's Response:
  Thanks to your suggestion, citations to three recent studies on marine biochemical parameters have been added to the latest manuscript.
  Line 26-27:
  ……protecting marine life(Fourrier et al., 2022).
  Line 43:
  ……vary significantly over short distances and time periods(Mignot et al., 2023).
  Line 49:
  ……high costs associated with data collection(Asselot et al., 2024).
  
  Citation: https://doi.org/10.5194/essd-2024-219-AC2

Xiaohan Zhang, Lizhe Wang, Jining Yan, and Sheng Wang

Data sets

Monthly Product of Marine Chemical Data in Australian Waters from 2000 to 2022 Xiaohan Zhang, Lizhe Wang, Jining Yan, and Sheng Wang https://doi.org/10.57760/sciencedb.09331

Xiaohan Zhang, Lizhe Wang, Jining Yan, and Sheng Wang

Viewed

Total article views: 642 (including HTML, PDF, and XML)

HTML	PDF	XML	Total	BibTeX	EndNote
510	103	29	642	29	36

HTML: 510
PDF: 103
XML: 29
Total: 642
BibTeX: 29
EndNote: 36

Views and downloads (calculated since 02 Jul 2024)

Month	HTML	PDF	XML	Total
Jul 2024	135	28	12	175
Aug 2024	66	13	3	82
Sep 2024	41	5	1	47
Oct 2024	33	3	3	39
Nov 2024	56	10	2	68
Dec 2024	45	17	1	63
Jan 2025	56	14	4	74
Feb 2025	43	10	2	55
Mar 2025	35	3	1	39

Cumulative views and downloads (calculated since 02 Jul 2024)

Month	HTML	PDF	XML	Total
Jul 2024	135	28	12	175
Aug 2024	66	13	3	82
Sep 2024	41	5	1	47
Oct 2024	33	3	3	39
Nov 2024	56	10	2	68
Dec 2024	45	17	1	63
Jan 2025	56	14	4	74
Feb 2025	43	10	2	55
Mar 2025	35	3	1	39

Viewed (geographical distribution)

Total article views: 635 (including HTML, PDF, and XML) Thereof 635 with geography defined and 0 with unknown origin.

Country	#	Views	%

Latest update: 01 Apr 2025

Short summary

This study tackles the limitations of marine biogeochemical monitoring by integrating CANYON-B and Random Forest regression methods to reconstruct key biogeochemical concentrations across Australia from 2000 to 2022 at a 1-kilometre resolution. Utilizing multi-source in-situ observations and MODIS products, this approach showcases machine learning's effectiveness in enhancing large-scale ocean chemistry data reconstruction.


Total:	0
HTML:	0
PDF:	0
XML:	0