Assimilation of remotely-sensed optical properties to improve marine biogeochemistry modelling

In this paper we evaluate whether the assimilation of remotely-sensed optical data into a marine ecosystem model improves the simulation of biogeochemistry in a shelf sea. A localized Ensemble Kalman filter was used to assimilate weekly diffuse light attenuation coefficient data, Kd(443) from SeaWiFs, into an ecosystem model of the western English Channel. The spatial distributions of (unassimilated) surface chlorophyll from satellite, and a multivariate time series of eighteen biogeochemical and optical variables measured in situ at one long-term monitoring site were used to evaluate the system performance for the year 2006. Assimilation reduced the root mean square error and improved the correlation with the assimilated Kd(443) observations, for both the analysis and, to a lesser extent, the forecast estimates, when compared to the reference model simulation. Improvements in the simulation of (unassimilated) ocean colour chlorophyll were less evident, and in some parts of the Channel the simulation of this data deteriorated. The estimation errors for the (unassimilated) in situ data were reduced for most variables with some exceptions, e.g. dissolved nitrogen. Importantly, the assimilation adjusted the balance of ecosystem processes by shifting the simulated food web towards the microbial loop, thus improving the estimation of some properties, e.g. total particulate carbon. Assimilation of Kd(443) outperformed a comparative chlorophyll assimilation experiment, in both the estimation of ocean colour data and in the simulation of independent in situ data. These results are related to relatively low error in Kd(443) data, and because it is a bulk optical property of marine ecosystems. Assimilation of remotely-sensed optical properties is a promising approach to improve the simulation of biogeochemical and optical variables that are relevant for ecosystem functioning and climate change studies.

Southern Annular Mode and westerly-wind-driven changes in Indian-Atlantic exchange mechanisms

The dynamical link between the Indian Ocean and Atlantic Meridional Overturning Circulation (AMOC) remains poorly understood. This partly arises from the complex Agulhas leakage, which occurs via rings, cyclones, and non-eddy flux. Hindcast simulations suggest that leakage has recently increased but have not decomposed this signal into its constituent mechanisms. Here these are isolated in a realistic ocean model. Increases in simulated leakage are attributed to stronger eddy and non-eddy-driven transports, and a strong warming and salinification, especially within Agulhas rings. Variability in both regimes is associated with strengthening Indian Ocean westerly winds, reflecting an increasingly positive Southern Annular Mode. While eddy and non-eddy flux signals are tied through turbulent eddy dissipation, the ratio between the two varies decadally. Consequently, while altimetry suggests a recent increase in retroflection turbulence and implied leakage, non-eddy flux may also play a significant role in modulating the leakage AMOC connection.

Spatio-temporal characteristics of Agulhas leakage: a model inter-comparison study

Investigating the variability of Agulhas leakage, the volume transport of water from the Indian Ocean to the South Atlantic Ocean, is highly relevant due to its potential contribution to the Atlantic Meridional Overturning Circulation as well as the global circulation of heat and salt and hence global climate. Quantifying Agulhas leakage is challenging due to the non-linear nature of this process; current observations are insufficient to estimate its variability and ocean models all have biases in this region, even at high resolution . An Eulerian threshold integration method is developed to examine the mechanisms of Agulhas leakage variability in six ocean model simulations of varying resolution. This intercomparison, based on the circulation and thermo- haline structure at the Good Hope line, a transect to the south west of the southern tip of Africa, is used to identify features that are robust regardless of the model used and takes into account the thermohaline biases of each model. When determined by a passive tracer method, 60 % of the magnitude of Agulhas leakage is captured and more than 80 % of its temporal fluctuations, suggesting that the method is appropriate for investigating the variability of Agulhas leakage. In all simulations but one, the major driver of variability is associated with mesoscale features passing through the section. High resolution (<1/10 deg.) hindcast models agree on the temporal (2–4 cycles per year) and spatial (300–500 km) scales of these features corresponding to observed Agulhas Rings. Coarser resolution models (<1/4 deg.) reproduce similar time scale of variability of Agulhas leakage in spite of their difficulties in representing the Agulhas rings properties. A coarser resolution climate model (2 deg.) does not resolve the spatio-temporal mechanism of variability of Agulhas leakage. Hence it is expected to underestimate the contribution of Agulhas Current System to climate variability.

Spatio-temporal characteristics of Agulhas leakage: a model inter-comparison study

Investigating the variability of Agulhas leakage, the volume transport of water from the Indian Ocean to the South Atlantic Ocean, is highly relevant due to its potential contribution to the Atlantic Meridional Overturning Circulation as well as the global circulation of heat and salt and hence global climate. Quantifying Agulhas leakage is challenging due to the non-linear nature of this process; current observations are insufficient to estimate its variability and ocean models all have biases in this region, even at high resolution . An Eulerian threshold integration method is developed to examine the mechanisms of Agulhas leakage variability in six ocean model simulations of varying resolution. This intercomparison, based on the circulation and thermo- haline structure at the Good Hope line, a transect to the south west of the southern tip of Africa, is used to identify features that are robust regardless of the model used and takes into account the thermohaline biases of each model. When determined by a passive tracer method, 60 % of the magnitude of Agulhas leakage is captured and more than 80 % of its temporal fluctuations, suggesting that the method is appropriate for investigating the variability of Agulhas leakage. In all simulations but one, the major driver of variability is associated with mesoscale features passing through the section. High resolution (<1/10 deg.) hindcast models agree on the temporal (2–4 cycles per year) and spatial (300–500 km) scales of these features corresponding to observed Agulhas Rings. Coarser resolution models (<1/4 deg.) reproduce similar time scale of variability of Agulhas leakage in spite of their difficulties in representing the Agulhas rings properties. A coarser resolution climate model (2 deg.) does not resolve the spatio-temporal mechanism of variability of Agulhas leakage. Hence it is expected to underestimate the contribution of Agulhas Current System to climate variability.

Impacts of the Oceans on Climate Change

The oceans play a key role in climate regulation especially in part buffering (neutralising) the effects of increasing levels of greenhouse gases in the atmosphere and rising global temperatures. This chapter examines how the regulatory processes performed by the oceans alter as a response to climate change and assesses the extent to which positive feedbacks from the ocean may exacerbate climate change. There is clear evidence for rapid change in the oceans. As the main heat store for the world there has been an accelerating change in sea temperatures over the last few decades, which has contributed to rising sea‐level. The oceans are also the main store of carbon dioxide (CO2), and are estimated to have taken up ∼40% of anthropogenic-sourced CO2 from the atmosphere since the beginning of the industrial revolution. A proportion of the carbon uptake is exported via the four ocean ‘carbon pumps’ (Solubility, Biological, Continental Shelf and Carbonate Counter) to the deep ocean reservoir. Increases in sea temperature and changing planktonic systems and ocean currents may lead to a reduction in the uptake of CO2 by the ocean; some evidence suggests a suppression of parts of the marine carbon sink is already underway. While the oceans have buffered climate change through the uptake of CO2 produced by fossil fuel burning this has already had an impact on ocean chemistry through ocean acidification and will continue to do so. Feedbacks to climate change from acidification may result from expected impacts on marine organisms (especially corals and calcareous plankton), ecosystems and biogeochemical cycles. The polar regions of the world are showing the most rapid responses to climate change. As a result of a strong ice–ocean influence, small changes in temperature, salinity and ice cover may trigger large and sudden changes in regional climate with potential downstream feedbacks to the climate of the rest of the world. A warming Arctic Ocean may lead to further releases of the potent greenhouse gas methane from hydrates and permafrost. The Southern Ocean plays a critical role in driving, modifying and regulating global climate change via the carbon cycle and through its impact on adjacent Antarctica. The Antarctic Peninsula has shown some of the most rapid rises in atmospheric and oceanic temperature in the world, with an associated retreat of the majority of glaciers. Parts of the West Antarctic ice sheet are deflating rapidly, very likely due to a change in the flux of oceanic heat to the undersides of the floating ice shelves. The final section on modelling feedbacks from the ocean to climate change identifies limitations and priorities for model development and associated observations. Considering the importance of the oceans to climate change and our limited understanding of climate-related ocean processes, our ability to measure the changes that are taking place are conspicuously inadequate. The chapter highlights the need for a comprehensive, adequately funded and globally extensive ocean observing system to be implemented and sustained as a high priority. Unless feedbacks from the oceans to climate change are adequately included in climate change models, it is possible that the mitigation actions needed to stabilise CO2 and limit temperature rise over the next century will be underestimated.

Ecohydrodynamics of cold-water coral reefs: A case study of the Mingulay Reef complex (western Scotland)

Ecohydrodynamics investigates the hydrodynamic constraints on ecosystems across different temporal and spatial scales. Ecohydrodynamics play a pivotal role in the structure and functioning of marine ecosystems, however the lack of integrated complex flow models for deep-water ecosystems beyond the coastal zone prevents further synthesis in these settings. We present a hydrodynamic model for one of Earth's most biologically diverse deep-water ecosystems, cold-water coral reefs. The Mingulay Reef Complex (western Scotland) is an inshore seascape of cold-water coral reefs formed by the scleractinian coral Lophelia pertusa. We applied single-image edge detection and composite front maps using satellite remote sensing, to detect oceanographic fronts and peaks of chlorophyll a values that likely affect food supply to corals and other suspension-feeding fauna. We also present a high resolution 3D ocean model to incorporate salient aspects of the regional and local oceanography. Model validation using in situ current speed, direction and sea elevation data confirmed the model's realistic representation of spatial and temporal aspects of circulation at the reef complex including a tidally driven current regime, eddies, and downwelling phenomena. This novel combination of 3D hydrodynamic modelling and remote sensing in deep-water ecosystems improves our understanding of the temporal and spatial scales of ecological processes occurring in marine systems. The modelled information has been integrated into a 3D GIS, providing a user interface for visualization and interrogation of results that allows wider ecological application of the model and that can provide valuable input for marine biodiversity and conservation applications.

Modelling the Stoichiometric Regulation of C-Rich Toxins in Marine Dinoflagellates

Toxin production in marine microalgae was previously shown to be tightly coupled with cellular stoichiometry. The highest values of cellular toxin are in fact mainly associated with a high carbon to nutrient cellular ratio. In particular, the cellular accumulation of C-rich toxins (i.e., with C:N > 6.6) can be stimulated by both N and P deficiency. Dinoflagellates are the main producers of C-rich toxins and may represent a serious threat for human health and the marine ecosystem. As such, the development of a numerical model able to predict how toxin production is stimulated by nutrient supply/deficiency is of primary utility for both scientific and management purposes. In this work we have developed a mechanistic model describing the stoichiometric regulation of C-rich toxins in marine dinoflagellates. To this purpose, a new formulation describing toxin production and fate was embedded in the European Regional Seas Ecosystem Model (ERSEM), here simplified to describe a monospecific batch culture. Toxin production was assumed to be composed by two distinct additive terms; the first is a constant fraction of algal production and is assumed to take place at any physiological conditions. The second term is assumed to be dependent on algal biomass and to be stimulated by internal nutrient deficiency. By using these assumptions, the model reproduced the concentrations and temporal evolution of toxins observed in cultures of Ostreopsis cf. ovata, a benthic/epiphytic dinoflagellate producing C-rich toxins named ovatoxins. The analysis of simulations and their comparison with experimental data provided a conceptual model linking toxin production and nutritional status in this species. The model was also qualitatively validated by using independent literature data, and the results indicate that our formulation can be also used to simulate toxin dynamics in other dinoflagellates. Our model represents an important step towards the simulation and prediction of marine algal toxicity.