The carbonate system of the NW European shelf: sensitivity and model validation

The ocean plays an important role in regulating the climate, acting as a sink for carbon dioxide, perturbing the carbonate system and resulting in a slow decrease of seawater pH. Understanding the dynamics of the carbonate system in shelf sea regions is necessary to evaluate the impact of Ocean Acidification (OA) in these societally important ecosystems. Complex hydrodynamic and ecosystem coupled models provide a method of capturing the significant heterogeneity of these areas. However rigorous validation is essential to properly assess the reliability of such models. The coupled model POLCOMS–ERSEM has been implemented in the North Western European shelf with a new parameterization for alkalinity explicitly accounting for riverine inputs and the influence of biological processes. The model has been validated in a like with like comparison with North Sea data from the CANOBA dataset. The model shows good to reasonable agreement for the principal variables, physical (temperature and salinity), biogeochemical (nutrients) and carbonate system (dissolved inorganic carbon and total alkalinity), but simulation of the derived variables, pH and pCO2, are not yet fully satisfactory. This high uncertainty is attributed mostly to riverine forcing and primary production. This study suggests that the model is a useful tool to provide information on Ocean Acidification scenarios, but uncertainty on pH and pCO2 needs to be reduced, particularly when impacts of OA on ecosystem functions are included in the model systems.

Simulating The Impact Of Technological-Change On Dissolved Cadmium Distribution In The Severn Estuary

An estuarine model is described which computes the dispersive and advective properties of the Severn Estuary. It was calibrated and validated using 50 measured salinity distributions and then used to predict the magnitude and sitings of the major inputs of dissolved cadmium levels throughout the estuary. The results provided an impetus for implementing tighter controls on effluents and for improving estimates of cadmium discharges from industrial sources. The model has also been used to investigate the sensitivity of the estuarine system to changes in dispersion; by considering large reductions in the dispersion coefficients it is hoped that the results might be indicative of the environmental consequences following the construction of a tidal power generating scheme.

Carbon fluxes and pelagic ecosystem dynamics near two western Antarctic Peninsula Adélie penguin colonies: an inverse model approach

An inverse food-web model for the western Antarctic Peninsula (WAP) pelagic food web was constrained with data from Palmer Long Term Ecological Research (PAL-LTER) project annual austral summer sampling cruises. Model solutions were generated for 2 regions with Adelie penguin Pygoscelis adeliae colonies presenting different population trends (a northern and a southern colony) for a 12 yr period (1995-2006). Counter to the standard paradigm, comparisons of carbon flow through bacteria, microzooplankton, and krill showed that the diatom-krill-top predator food chain is not the dominant pathway for organic carbon exchanges. The food web is more complex, including significant contributions by microzooplankton and the microbial loop. Using both inverse model results and network indices, it appears that in the northern WAP the food web is dominated by the microbial food web, with a temporal trend toward its increasing importance. The dominant pathway for the southern WAP food web varies from year to year, with no detectable temporal trend toward dominance of microzooplankton versus krill. In addition, sensitivity analyses indicated that the northern colony of Adelie penguins, whose population size has been declining over the past 35 yr, appears to have sufficient krill during summer to sustain its basic metabolic needs and rear chicks, suggesting the importance of other processes in regulating the Adelie population decline.

Generalized ocean color inversion model for retrieving marine inherent optical properties

Ocean color measured from satellites provides daily, global estimates of marine inherent optical properties (IOPs). Semi-analytical algorithms (SAAs) provide one mechanism for inverting the color of the water observed by the satellite into IOPs. While numerous SAAs exist, most are similarly constructed and few are appropriately parameterized for all water masses for all seasons. To initiate community-wide discussion of these limitations, NASA organized two workshops that deconstructed SAAs to identify similarities and uniqueness and to progress toward consensus on a unified SAA. This effort resulted in the development of the generalized IOP (GIOP) model software that allows for the construction of different SAAs at runtime by selection from an assortment of model parameterizations. As such, GIOP permits isolation and evaluation of specific modeling assumptions, construction of SAAs, development of regionally tuned SAAs, and execution of ensemble inversion modeling. Working groups associated with the workshops proposed a preliminary default configuration for GIOP (GIOP-DC), with alternative model parameterizations and features defined for subsequent evaluation. In this paper, we: (1) describe the theoretical basis of GIOP; (2) present GIOP-DC and verify its comparable performance to other popular SAAs using both in situ and synthetic data sets; and, (3) quantify the sensitivities of their output to their parameterization. We use the latter to develop a hierarchical sensitivity of SAAs to various model parameterizations, to identify components of SAAs that merit focus in future research, and to provide material for discussion on algorithm uncertainties and future emsemble applications.

Exploiting satellite earth observation to quantify current global oceanic DMS flux and its future climate sensitivity

We used coincident Envisat RA2 and AATSR temperature and wind speed data from 2008/2009 to calculate the global net sea-air flux of dimethyl sulfide (DMS), which we estimate to be 19.6 Tg S a21. Our monthly flux calculations are compared to open ocean eddy correlation measurements of DMS flux from 10 recent cruises, with a root mean square difference of 3.1 lmol m22 day21. In a sensitivity analysis, we varied temperature, salinity, surface wind speed, and aqueous DMS concentration, using fixed global changes as well as CMIP5 model output. The range of DMS flux in future climate scenarios is discussed. The CMIP5 model predicts a reduction in surface wind speed and we estimate that this will decrease the global annual sea-air flux of DMS by 22% over 25 years. Concurrent changes in temperature, salinity, and DMS concentration increase the global flux by much smaller amounts. The net effect of all CMIP5 modelled 25 year predictions was a 19% reduction in global DMS flux. 25 year DMS concentration changes had significant regional effects, some positive (Southern Ocean, North Atlantic, Northwest Pacific) and some negative (isolated regions along the Equator and in the Indian Ocean). Using satellite-detected coverage of coccolithophore blooms, our estimate of their contribution to North Atlantic DMS emissions suggests that the coccolithophores contribute only a small percentage of the North Atlantic annual flux estimate, but may be more important in the summertime and in the northeast Atlantic.

Assessing the sensitivity of subtidal sedimentary habitats to pressures associated with marine activities. Phase 1 Report: Rationale and proposed ecological groupings for Level 5 biotopes against which sensitivity assessments would be best undertaken.

The purpose of this study is to produce a series of Conceptual Ecological Models (CEMs) that represent sublittoral rock habitats in the UK. CEMs are diagrammatic representations of the influences and processes that occur within an ecosystem. They can be used to identify critical aspects of an ecosystem that may be studied further, or serve as the basis for the selection of indicators for environmental monitoring purposes. The models produced by this project are control diagrams, representing the unimpacted state of the environment free from anthropogenic pressures. It is intended that the models produced by this project will be used to guide indicator selection for the monitoring of this habitat in UK waters. CEMs may eventually be produced for a range of habitat types defined under the UK Marine Biodiversity Monitoring R&D Programme (UKMBMP), which, along with stressor models, are designed to show the interactions within impacted habitats, would form the basis of a robust method for indicator selection. This project builds on the work to develop CEMs for shallow sublittoral coarse sediment habitats (Alexander et al 2014). The project scope included those habitats defined as ‘sublittoral rock’. This definition includes those habitats that fall into the EUNIS Level 3 classifications A3.1 Atlantic and Mediterranean high energy infralittoral rock, A3.2 Atlantic and Mediterranean moderate energy infralittoral rock, A3.3 Atlantic and Mediterranean low energy infralittoral rock, A4.1 Atlantic and Mediterranean high energy circalittoral rock, A4.2 Atlantic and Mediterranean moderate energy circalittoral rock, and A4.3 Atlantic and Mediterranean low energy circalittoral rock as well as the constituent Level 4 and 5 biotopes that are relevant to UK waters. A species list of characterising fauna to be included within the scope of the models was identified using an iterative process to refine the full list of species found within the relevant Level 5 biotopes. A literature review was conducted using a pragmatic and iterative approach to gather evidence regarding species traits and information that would be used to inform the models and characterise the interactions that occur within the sublittoral rock habitat. All information gathered during the literature review was entered into a data logging pro-forma spreadsheet that accompanies this report. Wherever possible, attempts were made to collect information from UK-specific peer-reviewed studies, although other sources were used where necessary. All data gathered was subject to a detailed confidence assessment. Expert judgement by the project team was utilised to provide information for aspects of the models for which references could not be sourced within the project timeframe. A multivariate analysis approach was adopted to assess ecologically similar groups (based on ecological and life history traits) of fauna from the identified species to form the basis of the models. A model hierarchy was developed based on these ecological groups. One general control model was produced that indicated the high-level drivers, inputs, biological assemblages, ecosystem processes and outputs that occur in sublittoral rock habitats. In addition to this, seven detailed sub-models were produced, which each focussed on a particular ecological group of fauna within the habitat: ‘macroalgae’, ‘temporarily or permanently attached active filter feeders’, ‘temporarily or permanently attached passive filter feeders’, ‘bivalves, brachiopods and other encrusting filter feeders’, ‘tube building fauna’, ‘scavengers and predatory fauna’, and ‘non-predatory mobile fauna’. Each sub-model is accompanied by an associated confidence model that presents confidence in the links between each model component. The models are split into seven levels and take spatial and temporal scale into account through their design, as well as magnitude and direction of influence. The seven levels include regional to global drivers, water column processes, local inputs/processes at the seabed, habitat and biological assemblage, output processes, local ecosystem functions, and regional to global ecosystem functions. The models indicate that whilst the high level drivers that affect each ecological group are largely similar, the output processes performed by the biota and the resulting ecosystem functions vary both in number and importance between groups. Confidence within the models as a whole is generally high, reflecting the level of information gathered during the literature review. Physical drivers which influence the ecosystem were found to be of high importance for the sublittoral rock habitat, with factors such as wave exposure, water depth and water currents noted to be crucial in defining the biological assemblages. Other important factors such as recruitment/propagule supply, and those which affect primary production, such as suspended sediments, light attenuation and water chemistry and temperature, were also noted to be key and act to influence the food sources consumed by the biological assemblages of the habitat, and the biological assemblages themselves. Output processes performed by the biological assemblages are variable between ecological groups depending on the specific flora and fauna present and the role they perform within the ecosystem. Of particular importance are the outputs performed by the macroalgae group, which are diverse in nature and exert influence over other ecological groups in the habitat. Important output processes from the habitat as a whole include primary and secondary production, bioengineering, biodeposition (in mixed sediment habitats) and the supply of propagules; these in turn influence ecosystem functions at the local scale such as nutrient and biogeochemical cycling, supply of food resources, sediment stability (in mixed sediment habitats), habitat provision and population and algae control. The export of biodiversity and organic matter, biodiversity enhancement and biotope stability are the resulting ecosystem functions that occur at the regional to global scale. Features within the models that are most useful for monitoring habitat status and change due to natural variation have been identified, as have those that may be useful for monitoring to identify anthropogenic causes of change within the ecosystem. Biological, physical and chemical features of the ecosystem have been identified as potential indicators to monitor natural variation, whereas biological factors and those physical /chemical factors most likely to affect primary production have predominantly been identified as most likely to indicate change due to anthropogenic pressures.

Historical data reveal power-law dispersal patterns of invasive aquatic species.

Understanding how invasive species spread is of particular concern in the current era of globalisation and rapid environmental change. The occurrence of super-diffusive movements within the context of Lévy flights has been discussed with respect to particle physics, human movements, microzooplankton, disease spread in global epidemiology and animal foraging behaviour. Super-diffusive movements provide a theoretical explanation for the rapid spread of organisms and disease, but their applicability to empirical data on the historic spread of organisms has rarely been tested. This study focuses on the role of long-distance dispersal in the invasion dynamics of aquatic invasive species across three contrasting areas and spatial scales: open ocean (north-east Atlantic), enclosed sea (Mediterranean) and an island environment (Ireland). Study species included five freshwater plant species, Azolla filiculoides, Elodea canadensis, Lagarosiphon major, Elodea nuttallii and Lemna minuta; and ten species of marine algae, Asparagopsis armata, Antithamnionella elegans, Antithamnionella ternifolia, Codium fragile, Colpomenia peregrina, Caulerpa taxifolia, Dasysiphonia sp., Sargassum muticum, Undaria pinnatifida and Womersleyella setacea. A simulation model is constructed to show the validity of using historical data to reconstruct dispersal kernels. Lévy movement patterns similar to those previously observed in humans and wild animals are evident in the re-constructed dispersal pattern of invasive aquatic species. Such patterns may be widespread among invasive species and could be exacerbated by further development of trade networks, human travel and environmental change. These findings have implications for our ability to predict and manage future invasions, and improve our understanding of the potential for spread of organisms including infectious diseases, plant pests and genetically modified organisms.

The assessment of a global marine ecosystem model on the basis of emergent properties and ecosystem function: a case study with ERSEM

Ecosystem models are often assessed using quantitative metrics of absolute ecosystem state, but these model-data comparisons are disproportionately vulnerable to discrepancies in the location of important circulation features. An alternative method is to demonstrate the models capacity to represent ecosystem function; the emergence of a coherent natural relationship in a simulation indicates that the model may have an appropriate representation of the ecosystem functions that lead to the emergent relationship. Furthermore, as emergent properties are large-scale properties of the system, model validation with emergent properties is possible even when there is very little or no appropriate data for the region under study, or when the hydrodynamic component of the model differs significantly from that observed in nature at the same location and time. A selection of published meta-analyses are used to establish the validity of a complex marine ecosystem model and to demonstrate the power of validation with emergent properties. These relationships include the phytoplankton community structure, the ratio of carbon to chlorophyll in phytoplankton and particulate organic matter, the ratio of particulate organic carbon to particulate organic nitrogen and the stoichiometric balance of the ecosystem. These metrics can also inform aspects of the marine ecosystem model not available from traditional quantitative and qualitative methods. For instance, these emergent properties can be used to validate the design decisions of the model, such as the range of phytoplankton functional types and their behaviour, the stoichiometric flexibility with regards to each nutrient, and the choice of fixed or variable carbon to nitrogen ratios.