ALBERT

Beschreibung: 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.

Beschreibung: Skilful marine biogeochemical (BGC) models are required to understand a range of coastal and global phenomena such as changes in nitrogen and carbon cycles. The refinement of BGC models through the assimilation of variables calculated from observed in-water inherent optical properties (IOPs), such as phytoplankton absorption, is problematic. Empirically-derived relationships between IOPs and variables such as Chlorophyll-a concentration (Chl-a), Total Suspended Solids (TSS) and Color Dissolved Organic Matter (CDOM) have been shown to have errors that can exceed 100 % of the observed quantity. These errors are greatest in shallow coastal regions, such as the Great Barrier Reef (GBR), due the additional signal from bottom reflectance. Rather than assimilate quantities calculated using error-prone IOP algorithms, this study demonstrates the advantages of assimilating quantities calculated directly from the less error-prone satellite remote-sensing reflectance. The assimilation of a directly-observed quantity, in this case remote-sensing reflectance, is analogous to the assimilation of temperature brightness in Numerical Weather Prediction (NWP), or along-track sea-surface height in hydrodynamic models. To assimilate the observed reflectance, we use an in-water optical model to produce an equivalent simulated remote-sensing reflectance, and calculate the mis-match between the observed and simulated quantities to constrain the BGC model with a Deterministic Ensemble Kalman Filter (DEnKF). Using the assumption that simulated surface Chl-a is equivalent to remotely-sensed OC3M estimate of Chl-a resulted in a forecast error of approximately 75 %. Alternatively, assimilation of remote-sensing reflectance resulted in a forecast error of less than 40 %. Thus, in the coastal waters of the GBR, assimilating remote-sensing reflectance halved the forecast errors. When the analysis and forecast fields from the assimilation system are compared with the non-assimilating model, an independent comparison to in-situ observations of Chl-a, TSS, and dissolved inorganic nutrients (NO3, NH4 and DIP) show that errors are reduced by up to 90 %. In all cases, the assimilation system improves the result compared to the non-assimilating model. This approach allows for the incorporation of vast quantities of remote-sensing observations that have in the past been discarded due to shallow water and/or artefacts introduced by terrestrially-derived TSS and CDOM, or the lack of a calibrated regional IOP algorithm.

Beschreibung: Skillful marine biogeochemical (BGC) models are required to understand a range of coastal and global phenomena such as changes in nitrogen and carbon cycles. The refinement of BGC models through the assimilation of variables calculated from observed in-water inherent optical properties (IOPs), such as phytoplankton absorption, is problematic. Empirically derived relationships between IOPs and variables such as chlorophyll-a concentration (Chl a), total suspended solids (TSS) and coloured dissolved organic matter (CDOM) have been shown to have errors that can exceed 100 % of the observed quantity. These errors are greatest in shallow coastal regions, such as the Great Barrier Reef (GBR), due to the additional signal from bottom reflectance. Rather than assimilate quantities calculated using IOP algorithms, this study demonstrates the advantages of assimilating quantities calculated directly from the less error-prone satellite remote-sensing reflectance (RSR). To assimilate the observed RSR, we use an in-water optical model to produce an equivalent simulated RSR and calculate the mismatch between the observed and simulated quantities to constrain the BGC model with a deterministic ensemble Kalman filter (DEnKF). The traditional assumption that simulated surface Chl a is equivalent to the remotely sensed OC3M estimate of Chl a resulted in a forecast error of approximately 75 %. We show this error can be halved by instead using simulated RSR to constrain the model via the assimilation system. When the analysis and forecast fields from the RSR-based assimilation system are compared with the non-assimilating model, a comparison against independent in situ observations of Chl a, TSS and dissolved inorganic nutrients (NO3, NH4 and DIP) showed that errors are reduced by up to 90 %. In all cases, the assimilation system improves the simulation compared to the non-assimilating model. Our approach allows for the incorporation of vast quantities of remote-sensing observations that have in the past been discarded due to shallow water and/or artefacts introduced by terrestrially derived TSS and CDOM or the lack of a calibrated regional IOP algorithm.

Beschreibung: 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.