Search

Since the mid-1990s, Australia's Commonwealth Science Industry and Research Organisation (CSIRO) has been developing a biogeochemical (BGC) model for coupling with a hydrodynamic and sediment model for application in estuaries, coastal waters and shelf seas. The suite of coupled models is referred to as the CSIRO Environmental Modelling Suite (EMS) and has been applied at tens of locations around the Australian continent. At a mature point in the BGC model's development, this paper presents a full mathematical description, as well as links to the freely available code and user guide. The mathematical description is structured into processes so that the details of new parameterisations can be easily identified, along with their derivation. In EMS, the underwater light field is simulated by a spectrally resolved optical model that calculates vertical light attenuation from the scattering and absorption of 20+ optically active constituents. The BGC model itself cycles carbon, nitrogen, phosphorous and oxygen through multiple phytoplankton, zooplankton, detritus and dissolved organic and inorganic forms in multiple water column and sediment layers. The water column is dynamically coupled to the sediment to resolve deposition, resuspension and benthic–pelagic biogeochemical fluxes. With a focus on shallow waters, the model also includes detailed representations of benthic plants such as seagrass, macroalgae and coral polyps. A second focus has been on, where possible, the use of geometric derivations of physical limits to constrain ecological rates. This geometric approach generally requires population-based rates to be derived from initially considering the size and shape of individuals. For example, zooplankton grazing considers encounter rates of one predator on a prey field based on summing relative motion of the predator with the prey individuals and the search area; chlorophyll synthesis includes a geometrically derived self-shading term; and the bottom coverage of benthic plants is calculated from their biomass using an exponential form derived from geometric arguments. This geometric approach has led to a more algebraically complicated set of equations when compared to empirical biogeochemical model formulations based on populations. But while being algebraically complicated, the model has fewer unconstrained parameters and is therefore simpler to move between applications than it would otherwise be. The version of EMS described here is implemented in the eReefs project that delivers a near-real-time coupled hydrodynamic, sediment and biogeochemical simulation of the Great Barrier Reef, northeast Australia, and its formulation provides an example of the application of geometric reasoning in the formulation of aquatic ecological processes. [ABSTRACT FROM AUTHOR]

Since the mid 1990s, Australia's Commonwealth Science Industry and Research Organisation (CSIRO) has developed a biogeochemical (BGC) model for coupling with a hydrodynamic and sediment model for application in estuaries, coastal waters and shelf seas. The suite of coupled models is referred to as the CSIRO Environmental Modelling Suite (EMS) and has been applied at tens of locations around the Australian continent. At a mature point in the BGC model's development, this paper presents a full mathematical description, as well as links to the freely available code and User Guide. The mathematical description is structured into processes so that the details of new parameterisations can be easily identified, along with their derivation. The EMS BGC model cycles carbon, nitrogen, phosphorous and oxygen through multiple phytoplankton, zooplankton, detritus and dissolved organic and inorganic forms in multiple water column and sediment layers. The underwater light field is simulated by a spectrally-resolved optical model that includes the calculation of water-leaving reflectance for validation with remote sensing. The water column is dynamically coupled to the sediment to resolve deposition, resuspension and benthic-pelagic biogeochemical fluxes. With a focus on shallow waters, the model also includes particularly-detailed representations of benthic plants such as seagrass, macroalgae and coral polyps. A second focus has been on, where possible, the use of geometric derivations of physical limits to constrain ecological rates, which generally requires population-based rates to be derived from initially considering the size and shape of individuals. For example, zooplankton grazing considers encounter rates of one predator on a prey field based on summing relative motion of the predator with the prey individuals and the search area, chlorophyll synthesis includes a geometrically-derived self-shading term, and the bottom coverage of benthic plants is generically-related to their biomass using an exponential form derived from geometric arguments. This geometric approach has led to a more algebraically-complicated set of equations when compared to more empirical biogeochemical model formulations. But while being algebraically-complicated, the model has fewer unconstrained parameters and is therefore simpler to move between applications than it would otherwise be. The version of the biogeochemistry described here is implemented in the eReefs project that is delivering a near real time coupled hydrodynamic, sediment and biogeochemical simulation of the Great Barrier Reef, northeast Australia, and its formulation provides an example of the application of geometric reasoning in the formulation of aquatic ecological processes. [ABSTRACT FROM AUTHOR]

Accurate atmospheric correction (AC) is a prerequisite for quantitative ocean colour remote sensing and remains a challenge in particular over coastal waters. Commonly AC algorithms are validated by establishing a mean retrieval error from match-up analysis, which compares the satellite-derived surface reflectance with concurrent ground radiometric observations. Pixel-based reflectance uncertainties however, are rarely provided by AC algorithms and those for the operational Ocean and Land Colour Instrument (OLCI) marine reflectance product are not yet recommended for use. AC retrieval errors and uncertainties directly determine the quality with which ocean colour products can be estimated from the marine surface reflectance. Increasingly there is also the need for reflectance uncertainty products to be used as data assimilation inputs into biogeochemical models. This paper describes the development of a new coastal AC algorithm for Sentinel-3 OLCI that provides pixel-based estimation of the inherent model inversion uncertainty and sensor noise propagation. The algorithm is a full-spectral model-based inversion of radiative transfer (RT) simulations in a coupled atmosphere–ocean system using an ensemble of artificial neural networks (ANN) that were initialized differently during the training process, but composed of the same network architecture. The algorithm has been validated against in-situ radiometric observations across a wide range of optical water types, and has been compared with the latest EUMETSAT operational Level 2 processor IPF-OL-2 v7.01. In this analysis we found that the ensemble ANN showed improved performance over the operational Level 2 processor with a band-averaged (412–708 nm) mean absolute percentage error (MAPE) of 16% compared to 37% and a four-times lower band-averaged bias of -0.00045 sr-1. In the ensemble inversion process we account for three uncertainty components: (1) the total model variance that describes the variance of the data from the different ANNs, (2) the prediction variance of the mean, which is based on calculations of the RT simulations and (3) the instrument noise variance of the mean by propagating the OLCI spectral signal-to-noise ratios (SNR). To study algorithm performance and to quantify the contribution of the different uncertainty components to the total uncertainty, we applied the algorithm to an optically complex full resolution (FR) test scene covering coastal waters of the Great Barrier Reef, Australia. The uncertainties associated with the instrument noise variance were found to be two orders of magnitude lower than the uncertainty components of the prediction and total model variances. The overall largest uncertainty component in our uncertainty framework is attributed to the total model inversion error from averaging the responses of the slightly different adapted networks in the ensemble. The algorithm is made publicly available as a Python/C plugin for the Sentinel Application Platform (SNAP). • An ANN AC-algorithm that provides pixel-based uncertainty estimation was developed. • The algorithm shows advanced performance compared to the operational processor. • An analysis of the uncertainty components for different water types is provided. • Uncertainty maps provide a useful tool for future algorithm improvement. • The algorithm is publicly available as a plugin for the Sentinel Application Platform. [ABSTRACT FROM AUTHOR]