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 geom

Estuaries are increasingly degraded due to coastal urban development and are prone to hypoxia problems. The macro-tidal Gironde Estuary is characterized by a highly concentrated turbidity maximum zone (TMZ). Field observations show that hypoxia occurs in summer in the TMZ at low river flow and a few days after the spring tide peak. In situ data highlight lower dissolved oxygen (DO) concentrations around the city of Bordeaux, located in the upper estuary. Interactions between multiple factors limit the understanding of the processes controlling the dynamics of hypoxia. A 3D biogeochemical model was developed, coupled with hydrodynamics and a sediment transport model, to assess the contribution of the TMZ and the impact of urban effluents through wastewater treatment plants (WWTPs) and sewage overflows (SOs) on hypoxia. Our model describes the transport of solutes and suspended material and the biogeochemical mechanisms impacting oxygen: primary production, degradation of all organic matter (i.e. including phytoplankton respiration, degradation of river and urban watershed matter), nitrification, and gas exchange. The composition and the degradation rates of each variable were characterized by in situ measurements and experimental data from the study area. The DO model was validated against observations in Bordeaux City. The simulated DO concentrations show good agreement with field observations and satisfactorily reproduce the seasonal and neap-spring time scale variations around the city of Bordeaux. Simulations show a spatial and temporal correlation between the formation of summer hypoxia and the location of the TMZ, with minimum DO centered in the vicinity of Bordeaux. To understand the contribution of the urban watershed forcing, different simulations with the presence or absence of urban effluents were compared. Our results show that in summer, a reduction of POC from SO would increase the DO minimum in the vicinity of Bordeaux by 3% of saturation. Omitting dis

The overarching goal of the GISERA marine environmental research program has been to make possible more accurate prediction and understanding of impacts and trends in water quality as well as ecological responses in primary producers (seagrass) and grazers (turtles) that have been assessed as being vulnerable due to expansion of development in the Port of Gladstone. In doing this is has also been the aim of the GISERA marine research program to develop tools that can be used to determine management options that may lead to the reduction of impacts on these key ecological assets in the future, well beyond the Port Curtis and the current phase of development. The GISERA Marine project has made significant progress in integrating environmental and ecological knowledge and towards providing tools, notably a re-locatable seagrass growth model, and a turtle shipping-risk assessment model, that provide for a synoptic picture of conditions within the harbour as well as risks to its key ecological elements. The two major sub-components in the project, 1) Sustaining turtles, dugongs and their habitat - an integrated marine observation system, and 2) Integrated modelling, are presented in five chapters, starting with observations of the biophysical properties of the water column (Chapter I) and seagrasses (Chapter II) that are brought together in a biogeochemical model of seagrass growth (Chapter III), moving up the food web to studies of turtle movements (Chapter IV)and modelling of turtle habitat use and risk from shipping (Chapter V). The optical data detailed in Chapter I has been vitally important for the development and validation of modelling of coastal waters in Gladstone Harbour, as well as more widely. The data provides confidence in both the use of GISERA optical samples for terrestrially-sourced fine sediment parameterisation in the optical model, as well as in the calculations for converting total suspended solids to remotely-sensed reflectance. The improved confi

Integrated multi-trophic aquaculture (IMTA) is being explored on both economic and environmental grounds in many traditional aquaculture regions. To test a variety of suitable macroalgae species and management scenarios, a numerical model is developed to quantify the remediation of dissolved nutrients and production of macroalgae near a nutrient source. Differences in the morphological, physiological, and economic characteristics of different macroalgae species can provide flexibility when considering the cost and benefit of farming macroalgae. Results show that of the three species studied, Macrocystis pyrifera removed 75 % of dissolved inorganic nitrogen (DIN) input from a point source, while Porphyra umbilicalis and Ulva lactuca removed 5 %. Both M. pyrifera and P. umbilicalis have reduced bioremediation capacity at increasing flow rates. U. lactuca showed increased bioremediation potential as flow rate increased from low to moderate flows. Increasing the optical depth increased the bioremediation potential of M. pyrifera for moderate values of the light attenuation coefficient, whereas bioremediation was unaffected by optical depth for both U. lactuca and P. umbilicalis. Harvesting increased bioremediation capacity of all species by up to 25-fold dependent on the establishment phase and harvesting frequency. We conclude that the choice of macroalgae species greatly affects the success of IMTA and that both harvesting and farm arrangements can be used to greatly optimize bioremediation.

PROWQM, a 1-D depth resolving model which couples physical and microbiological processes in the water column with sedimentation/resuspension and benthic mineralisation processes, has been used to simulate seasonal changes of chlorophyll, nutrients and oxygen at the PROVESS north site (59degrees20' N 1 degrees00' E) in the North Sea. PROWQM is derived from the 3-D model COHERENS, and improves COHEREN's benthic and pelagic biology. The physical sub-model of PROWQM implicitly solves turbulence closure equations forced by climatological, or realistic high-frequency, meteorological and tidal data. The pelagic biological sub-model 2MPPD includes a 'diatomy' microplankton (mp1) and a 'flagellatey' (or microbial loop) microplankton (mp2), the cycling of silicon and nitrogen, slow-sinking detritus, and fast-sinking phytodetritus. Phytodetritus is formed by shear-driven aggregation of particulate material, using a simple algorithm. for bulk processes that is derived by considering the interactions of single cells. The microplankton compartments include heterotrophic bacteria and protozoa as well as phytoplankton, and most microplankton rates are specified with the aid of a 'heterotroph fraction' parameter, which was 0.125 for mp1 and 0.6 for mp2. The microbiological system is closed by mesozooplankton grazing pressures imposed as time varying series determined from observed zooplankton abundance. The benthic boundary sub-model includes a superficial fluff layer and a nutrient element reservoir in the consolidated sediment. Particulate material in the fluff layer can be resuspended (in response to bed stress by near-bed flows), mineralised or carried by bioturbation into the underlying, consolidated, sediment, where it is mineralised and its nutrients returned to the water-column at rates mainly dependent on (implicit) macrobenthic pumping. Benthic denitrification can occur when mineralisation rates exceed oxygen supply. Verification of the PROWQM numerical implementation used