Stefano Ciavatta

A comparison between the uncertainties in model parameters and in forcing functions: its application to a 3D water-quality model

This paper illustrates the application of both local and global sensitivity analysis techniques to an estimation of the uncertainty in the output of a 3D reaction-diffusion ecological model; the model describes the seasonal dynamics of dissolved Nitrogen and Phosphorous, and those of the phytoplanktonic and zooplanktonic communities in the lagoon of Venice. Two sources of uncertainty were taken into account and compared: (1) uncertainty concerning the parameters of the governing equation; (2) uncertainty concerning the forcing functions. The mean annual concentrations of Dissolved Inorganic Nitrogen (DIN) was regarded as the model output, as it represents the largest fraction of the Total Dissolved Nitrogen, TDN, for which the current Italian legislation sets a quality target in the lagoon of Venice. A local sensitivity analysis was initially used, so as to rank the parameters and provide an initial estimation of the uncertainty, which is a result of an imperfect knowledge of the dynamic of the system. This uncertainty was compared with that induced by an imperfect knowledge of the loads of Nitrogen, which represent the main forcing functions. On the basis of the results of the local analysis, the most important parameters and loads were then taken as the sources of uncertainty, in an attempt to assess their relative contributions. The global uncertainty and sensitivity analyses were carried out by means of a sampling-based Monte Carlo method. The results of the subsequent input-output regression analysis suggest that the variance in the model output could be partitioned among the sources of uncertainty, in accordance with a linear model. Based on this model, 79% of the variance in the mean annual concentration of DIN was accounted for by the uncertainty in the parameters which specify the dynamics of the phytoplankton and zooplankton, and only 5% by the uncertainties in the three main Nitrogen sources.

The Extended Kalman Filter (EKF) as a tool for the assimilation of high frequency water quality data

The Extended Kalman Filter (EKF) was applied to the analysis of high frequency field measurements of dissolved oxygen (DO), water temperature, salinity, collected by multiparametric sensors in the lagoon of Venice. This paper focuses on the practical aspects of the implementation of the EKF as a data assimilation technique and does not deal with the problems associated with the identification of the model. In this regard, the EKF has proved to be a useful tool for the updating of the estimates of the parameters of a simple DO-chlorophyll model, which can be used for linking the high frequency data to meteorological forcings, such as solar radiation and wind, and to other low frequency measurements of water quality parameters, such as the concentrations of Chlorophyll a and nutrients. The model can subsequently be used as a tool for checking the consistency of all this data, and may also be employed for controlling the quality of the data collected by the multiparametric sensors.

Global Uncertainty and Sensitivity Analysis of a Food-web Bioaccumulation Model

A global uncertainty and sensitivity analysis (UA/SA) of a state-of-the-art, food-web bioaccumulation model was carried out. We used an efficient screening analysis technique to identify the subset of the most relevant input factors among the whole set of 227 model parameters. A quantitative UA/SA was then applied to this subset to rank the relevance of the parameters and to partition the variance of the model output among them by means of a nonlinear regression of the outcomes of 1,000 Monte Carlo simulations. The concentrations of four representative persistent organic pollutants (POPs) in two representative species of the coastal marine food web of the Lagoon of Venice (Italy) were taken as model outputs. The screening analysis showed that the ranking was remarkably different in relation to the species and chemical being considered. The subsequent Monte Carlo-based quantitative analysis pointed out that the relationships among some of the parameters and the model outputs were nonlinear. The nonlinear regression showed that the fraction of output variance accounted for by each parameter was strongly dependent on the range of the octanol-water partition coefficient (K(OW)) values being considered. For the less hydrophobic chemicals, the main sources of model uncertainty were the parameters related to the respiratory bioaccumulation, whereas for the more hydrophobic ones, K(OW) and the other parameters related to the dietary uptake explained the largest fractions of the variance of the chemical concentrations in the organisms. The analysis highlighted that efforts are still needed for reducing uncertainty of model parameters to get reliable results from the application of food web bioaccumulation models.

An Objective Framework to Test the Quality of Candidate Indicators of Good Environmental Status

Large efforts are on-going within the EU to prepare the Marine Strategy Framework Directive’s (MSFD) assessment of the environmental status of the European seas. This assessment will only be as good as the indicators chosen to monitor the eleven descriptors of good environmental status (GEnS). An objective and transparent framework to determine whether chosen indicators actually support the aims of this policy is, however, not yet in place. Such frameworks are needed to ensure that the limited resources available to this assessment optimize the likelihood of achieving GEnS within collaborating states. Here, we developed a hypothesis-based protocol to evaluate whether candidate indicators meet quality criteria explicit to the MSFD, which the assessment community aspires to. Eight quality criteria are distilled from existing initiatives, and a testing and scoring protocol for each of them is presented. We exemplify its application in three worked examples, covering indicators for three GEnS descriptors (1, 5 and 6), various habitat components (seaweeds, seagrasses, benthic macrofauna and plankton), and assessment regions (Danish, Lithuanian and UK waters). We argue that this framework provides a necessary, transparent and standardized structure to support the comparison of candidate indicators, and the decision-making process leading to indicator selection. Its application could help identify potential limitations in currently available candidate metrics and, in such cases, help focus the development of more adequate indicators. Use of such standardized approaches will facilitate the sharing of knowledge gained across the MSFD parties despite context-specificity across assessment regions, and support the evidence-based management of European seas.

Uncertainty in Ocean-Color Estimates of Chlorophyll for Phytoplankton Groups

Over the past decade, techniques have been presented to derive the community structure of phytoplankton at synoptic scales using satellite ocean-colour data. There is a growing demand from the ecosystem modelling community to use these products for model evaluation and data assimilation. Yet, from the perspective of an ecosystem modeller these products are of limited use unless: (i) the phytoplankton products provided by the remote-sensing community match those required by the ecosystem modellers; and (ii) information on per-pixel uncertainty is provided to evaluate data quality. Using a large dataset collected in the North Atlantic, we re-tune a method to estimate the chlorophyll concentration of three phytoplankton groups, partitioned according to size (pico- (20μm)). The method is modified to account for the influence of sea surface temperature, also available from satellite data, on model parameters and on the partitioning of microphytoplankton into diatoms and dinoflagellates, such that the phytoplankton groups provided match those simulated in a state of the art marine ecosystem model (the European Regional Seas Ecosystem Model, ERSEM). The method is validated using another dataset, independent of the data used to parameterise the method, of more than 800 satellite and in situ match-ups. Using fuzzy-logic techniques for deriving per-pixel uncertainty, developed within the ESA Ocean Colour Climate Change Initiative (OC-CCI), the match-up dataset is used to derive the root mean square error and the bias between in situ and satellite estimates of the chlorophyll for each phytoplankton group, for 14 different optical water types (OWT). These values are then used with satellite estimates of OWTs to map uncertainty in chlorophyll on a per pixel basis for each phytoplankton group. It is envisaged these satellite products will be useful for those working on the validation of, and assimilation of data into, marine ecosystem models that simulate different phytoplankton groups.

Decadal reanalysis of biogeochemical indicators and fluxes in the North West European shelf‐sea ecosystem

This paper presents the first decadal reanalysis simulation of the biogeochemistry of the North West European shelf, along with a full evaluation of its skill, confidence, and value. An error-characterized satellite product for chlorophyll was assimilated into a physical-biogeochemical model of the North East Atlantic, applying a localized Ensemble Kalman filter. The results showed that the reanalysis improved the model simulation of assimilated chlorophyll in 60% of the study region. Model validation metrics showed that the reanalysis had skill in matching a large data set of in situ observations for 10 ecosystem variables. Spearman rank correlations were significant and higher than 0.7 for physical-chemical variables (temperature, salinity, and oxygen), ∼0.6 for chlorophyll and nutrients (phosphate, nitrate, and silicate), and significant, though lower in value, for partial pressure of dissolved carbon dioxide (∼0.4). The reanalysis captured the magnitude of pH and ammonia observations, but not their variability. The value of the reanalysis for assessing environmental status and variability has been exemplified in two case studies. The first shows that between 325,000 and 365,000 km2 of shelf bottom waters were vulnerable to oxygen deficiency potentially threatening bottom fishes and benthos. The second application confirmed that the shelf is a net sink of atmospheric carbon dioxide, but the total amount of uptake varies between 36 and 46 Tg C yr−1 at a 90% confidence level. These results indicate that the reanalysis output data set can inform the management of the North West European shelf ecosystem, in relation to eutrophication, fishery, and variability of the carbon cycle.

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.

Sensitivity analysis as a tool for the implementation of a water quality regulation based on the Maximum Permissible Loads policy

Abstract This paper shows how local sensitivity analysis, in respect of the parameters which specify the boundary conditions, can be used for relating the total load of non-conservative pollutants to their distributions within a water body. The method is applied to the estimation of the Maximum Permissible Load of inorganic nitrogen in the lagoon of Venice, that is of the maximum load of nitrogen which keeps its average yearly concentration below a prescribed threshold. The use of the spatial distributions of sensitivity coefficients in order to rank the sources of pollution and to forecast the effect of a reduction in the pollution is also discussed.

Assimilation of Ocean‐Color Plankton Functional Types to Improve Marine Ecosystem Simulations

We assimilated plankton functional types (PFTs) derived from ocean colour into a marine ecosystem model, to improve the simulation of biogeochemical indicators and emerging properties in a shelf sea. Error-characterized chlorophyll concentrations of four PFTs (diatoms, dinoflagellates, nanoplankton and picoplankton), as well as total chlorophyll for comparison, were assimilated into a physical-biogeochemical model of the North East Atlantic, applying a localized Ensemble Kalman filter. The reanalysis simulations spanned the years 1998 to 2003. The skill of the reference and reanalysis simulations in estimating ocean colour and in situ biogeochemical data were compared by using robust statistics. The reanalysis outperformed both the reference and the assimilation of total chlorophyll in estimating the ocean-colour PFTs (except nanoplankton), as well as the not-assimilated total chlorophyll, leading the model to simulate better the plankton community structure. Crucially, the reanalysis improved the estimates of not-assimilated in situ data of PFTs, as well as of phosphate and pCO2, impacting the simulation of the air-sea carbon flux. However, the reanalysis increased further the model overestimation of nitrate, in spite of increases in plankton nitrate uptake. The method proposed here is easily adaptable for use with other ecosystem models that simulate PFTs, for, e.g., reanalysis of carbon fluxes in the global ocean and for operational forecasts of biogeochemical indicators in shelf-sea ecosystems.

The Assimilation of Phytoplankton Functional Types for Operational Forecasting in the Northwest European Shelf

This paper proposes the use of assimilation of phytoplankton functional types (PFTs) surface chlorophyll for operational forecasting of biogeochemistry on the North‐West European (NWE) Shelf. We explicitly compare the 5‐day forecasting skill of three runs of a physical‐biogeochemical model: (a) a free reference run, (b) a run with daily data assimilation (DA) of total surface chlorophyll (ChlTot), and (c) a run with daily PFTs DA. We show that small total chlorophyll model bias hides comparatively large biases in PFTs chlorophyll, which ChlTot DA fails to correct. This is because the ChlTot DA splits the assimilated total chlorophyll into PFTs by preserving their simulated ratios, rather than taking account of the observed PFT concentrations. Unlike ChlTot DA, PFTs DA substantially improves model representation of PFTs chlorophyll. During forecasting the DA reanalysis skill in representing PFTs chlorophyll degrades toward the free run skill; however, PFTs DA outperforms free run within the whole 5‐day forecasting period. We validated our results with in situ data, and we demonstrated that (in both DA cases) the DA substantially improves the model representation of CO2 fugacity (PFTs DA more than ChlTot DA). ChlTot DA has a positive impact on the representation of silicate, while the PFTs DA seems to have a negative impact. The impact of DA on nitrate and phosphate is not significant. The implications of using a univariate assimilation method, which preserves the phytoplankton stochiometry, and the impact of model biases on the nonassimilated variables are discussed.