The exposure of the Great Barrier Reef to ocean acidification

© 2016, Nature Publishing Group. All rights reserved. The Great Barrier Reef (GBR) is founded on reef-building corals. Corals build their exoskeleton with aragonite, but ocean acidification is lowering the aragonite saturation state of seawater (Ωa). The downscaling of ocean acidification projections from global to GBR scales requires the set of regional drivers controlling Ωa to be resolved. Here we use a regional coupled circulation-biogeochemical model and observations to estimate the Ωa experienced by the 3,581 reefs of the GBR, and to apportion the contributions of the hydrological cycle, regional hydrodynamics and metabolism on Ωa variability. We find more detail, and a greater range (1.43), than previously compiled coarse maps of Ωa of the region (0.4), or in observations (1.0). Most of the variability in Ωa is due to processes upstream of the reef in question. As a result, future decline in Ωa is likely to be steeper on the GBR than currently projected by the IPCC assessment report

CSIRO Environmental Modelling Suite (EMS): scientific description of the optical and biogeochemical models (vB3p0)

Abstract. 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. </jats:p

Satellite data assimilation and estimation of a 3D coastal sediment transport model using error-subspace emulators

This paper describes sequential assimilation of data into a three-dimensional coastal ocean model using fast and cheap statistical surrogates of the model (emulators). The model simulates resuspension and deposition of fine sediments in a macro-tidal environment of the Fitzroy Estuary and Keppel Bay, North-East Australian coast. The assimilation algorithm was applied first to synthetic observations produced by a twin model run, and then with real data obtained from satellite observation. The latter are derived from remote sensing algorithms customised to the study region. The main objective of simulations was to test the data assimilation scheme using synthetic observations and identify potential issues and challenges when assimilating real data sets. The assimilation algorithm proved capable of substantially reducing a prior uncertainty of the model for both the scenario with the synthetic observations and the scenario with the satellite data. Significant remaining error in western Keppel Bay after assimilating satellite data is diagnostic of an underlying error in the system conceptualisation in other words, it indicates that the primary source of error is not in the parameter values specified, but in the model structure, in the interpretation of satellite data or in the other input data. The results of our study show the utility of the developed technique for the data assimilation into the three-dimensional sediment transport model of the Fitzroy estuary and Keppel Bay. More research is required to understand the capacity of this technique to generalise to other models and regions. Crown Copyright (C) 2012 Published by Elsevier Ltd. All rights reserved