A probabilistic model estimating oil spill clean-up costs - a case study for the Gulf of Finland

Existing models estimating oil spill costs at sea are based on data from the past, and they usually lack a systematic approach. This make them passive, and limits their ability to forecast the effect of the changes in the oil combating fleet or location of a spill on the oil spill costs. In this paper we make an attempt towards the development of a probabilistic and systematic model estimating the costs of clean-up operations for the Gulf of Finland. For this purpose we utilize expert knowledge along with the available data and information from literature. Then, the obtained information is combined into a framework with the use of a Bayesian Belief Networks. Due to lack of data, we validate the model by comparing its results with existing models, with which we found good agreement. We anticipate that the presented model can contribute to the cost-effective oil-combating fleet optimization for the Gulf of Finland. It can also facilitate the accident consequences estimation in the framework of formal safety assessment (FSA).

Decade Wind and Wave Data for NW Iberia, Cabo Sillerio

Seafloor sediment mobilization on the inner Northwest Iberian continental shelf is caused largely by ocean surface waves. The temporal and spatial variability in the wave height, wave period, and wave direction has a profound effect on local sediment mobilization, leading to distinct sediment mobilization scenarios. Six grain-size specific sediment mobilization scenarios, representing seasonal average and storm conditions, were simulated with a physics-based numerical model. Model inputs included meteorological and oceanographic data in conjunction with seafloor grain-size and the shelf bathymetric data. The results show distinct seasonal variations, most importantly in wave height, leading to sediment mobilization, specifically on the inner shelf shallower than 30 m water depth where up to 49% of the shelf area is mobilized. Medium to severe storm events are modeled to mobilize up to 89% of the shelf area above 150 m water depth. The frequency of each of these seasonal and storm-related sediment mobilization scenarios is addressed using a decade of meteorological and oceanographic data. The temporal and spatial patterns of the modeled sediment mobilization scenarios are discussed in the context of existing geological and environmental processes and conditions to assist scientific, industrial and environmental efforts that are directly affected by sediment mobilization. Examples, where sediment mobilization plays a vital role, include seafloor nutrient advection, recurrent arrival of oil from oil-spill-laden seafloor sediment, and bottom trawling impacts.

Organic compounds and bitumen in surface sediments from the Pechora Bay

Detailed geochemical investigations of bottom sediments in the Pechora Bay were carried out under a monitoring regime. The regional geochemical background is characterized. Organic matter of sediments is found to be clastic, formed under the influence of Pechora River run-off, and determined by genetic and lithofacial factors. An application of geochemical methodology is suggested for estimating hydrocarbon contamination of bottom sediments. As an example consequences of the accidental Usinsky oil spill are studied. They are based on analysis of composition of molecular markers. It is shown that increasing technogenic impact on Pechora Bay sediments is fixed only at molecular level and is not disastrous.

Raster maps for 29 environmental variables in three geographical regions

Aim: Greater understanding of the processes underlying biological invasions is required to determine and predict invasion risk. Two subspecies of olive (Olea europaea subsp. europaea and Olea europaea subsp. cuspidata) have been introduced into Australia from the Mediterranean Basin and southern Africa during the 19th century. Our aim was to determine to what extent the native environmental niches of these two olive subspecies explain the current spatial segregation of the subspecies in their non-native range. We also assessed whether niche shifts had occurred in the non-native range, and examined whether invasion was associated with increased or decreased occupancy of niche space in the non-native range relative to the native range. Location: South-eastern Australia, Mediterranean Basin and southern Africa. Methods: Ecological niche models (ENMs) were used to quantify the similarity of native and non-native realized niches. Niche shifts were characterized by the relative contribution of niche expansion, stability and contraction based on the relative occupancy of environmental space by the native and non-native populations. Results: Native ENMs indicated that the spatial segregation of the two subspecies in their non-native range was partly determined by differences in their native niches. However, we found that environmentally suitable niches were less occupied in the non-native range relative to the native range, indicating that niche shifts had occurred through a contraction of the native niches after invasion, for both subspecies. Main conclusions: The mapping of environmental factors associated with niche expansion, stability or contraction allowed us to identify areas of greater invasion risk. This study provides an example of successful invasions that are associated with niche shifts, illustrating that introduced plant species are sometimes readily able to establish in novel environments. In these situations the assumption of niche stasis during invasion, which is implicitly assumed by ENMs, may be unreasonable.

Environmental variables and multivariate seabed classification maps via k means clustering of Potter Cove, Antarctica, link to input and result files in GeoTIFF format

This study subdivides the Potter Cove, King George Island, Antarctica, into seafloor regions using multivariate statistical methods. These regions are categories used for comparing, contrasting and quantifying biogeochemical processes and biodiversity between ocean regions geographically but also regions under development within the scope of global change. The division obtained is characterized by the dominating components and interpreted in terms of ruling environmental conditions. The analysis includes in total 42 different environmental variables, interpolated based on samples taken during Australian summer seasons 2010/2011 and 2011/2012. The statistical errors of several interpolation methods (e.g. IDW, Indicator, Ordinary and Co-Kriging) with changing settings have been compared and the most reasonable method has been applied. The multivariate mathematical procedures used are regionalized classification via k means cluster analysis, canonical-correlation analysis and multidimensional scaling. Canonical-correlation analysis identifies the influencing factors in the different parts of the cove. Several methods for the identification of the optimum number of clusters have been tested and 4, 7, 10 as well as 12 were identified as reasonable numbers for clustering the Potter Cove. Especially the results of 10 and 12 clusters identify marine-influenced regions which can be clearly separated from those determined by the geological catchment area and the ones dominated by river discharge.

Data compilation on the biological response to ocean acidification: environmental and experimental context of data sets and related literature

The exponential growth of studies on the biological response to ocean acidification over the last few decades has generated a large amount of data. To facilitate data comparison, a data compilation hosted at the data publisher PANGAEA was initiated in 2008 and is updated on a regular basis (doi:10.1594/PANGAEA.149999). By January 2015, a total of 581 data sets (over 4 000 000 data points) from 539 papers had been archived. Here we present the developments of this data compilation five years since its first description by Nisumaa et al. (2010). Most of study sites from which data archived are still in the Northern Hemisphere and the number of archived data from studies from the Southern Hemisphere and polar oceans are still relatively low. Data from 60 studies that investigated the response of a mix of organisms or natural communities were all added after 2010, indicating a welcomed shift from the study of individual organisms to communities and ecosystems. The initial imbalance of considerably more data archived on calcification and primary production than on other processes has improved. There is also a clear tendency towards more data archived from multifactorial studies after 2010. For easier and more effective access to ocean acidification data, the ocean acidification community is strongly encouraged to contribute to the data archiving effort, and help develop standard vocabularies describing the variables and define best practices for archiving ocean acidification data.

Wetland maps including open water extent dynamics based on ENVISAT ASAR WS for Siberia, 2007 and 2008, links to GeoTIFFs

Wetlands store large amounts of carbon, and depending on their status and type, they release specific amounts of methane gas to the atmosphere. The connection between wetland type and methane emission has been investigated in various studies and utilized in climate change monitoring and modelling. For improved estimation of methane emissions, land surface models require information such as the wetland fraction and its dynamics over large areas. Existing datasets of wetland dynamics present the total amount of wetland (fraction) for each model grid cell, but do not discriminate the different wetland types like permanent lakes, periodically inundated areas or peatlands. Wetland types differently influence methane fluxes and thus their contribution to the total wetland fraction should be quantified. Especially wetlands of permafrost regions are expected to have a strong impact on future climate due to soil thawing. In this study ENIVSAT ASAR Wide Swath data was tested for operational monitoring of the distribution of areas with a long-term SW near 1 (hSW) in northern Russia (SW = degree of saturation with water, 1 = saturated), which is a specific characteristic of peatlands. For the whole northern Russia, areas with hSW were delineated and discriminated from dynamic and open water bodies for the years 2007 and 2008. The area identified with this method amounts to approximately 300,000 km**2 in northern Siberia in 2007. It overlaps with zones of high carbon storage. Comparison with a range of related datasets (static and dynamic) showed that hSW represents not only peatlands but also temporary wetlands associated with post-forest fire conditions in permafrost regions. Annual long-term monitoring of change in boreal and tundra environments is possible with the presented approach. Sentinel-1, the successor of ENVISAT ASAR, will provide data that may allow continuous monitoring of these wetland dynamics in the future complementing global observations of wetland fraction.

Secchi depth maps for the Baltic Sea, 2003-2012, with links to GeoTIFs

Secchi depth is a measure of water transparency. In the Baltic Sea region, Secchi depth maps are used to assess eutrophication and as input for habitat models. Due to their spatial and temporal coverage, satellite data would be the most suitable data source for such maps. But the Baltic Sea's optical properties are so different from the open ocean that globally calibrated standard models suffer from large errors. Regional predictive models that take the Baltic Sea's special optical properties into account are thus needed. This paper tests how accurately generalized linear models (GLMs) and generalized additive models (GAMs) with MODIS/Aqua and auxiliary data as inputs can predict Secchi depth at a regional scale. It uses cross-validation to test the prediction accuracy of hundreds of GAMs and GLMs with up to 5 input variables. A GAM with 3 input variables (chlorophyll a, remote sensing reflectance at 678 nm, and long-term mean salinity) made the most accurate predictions. Tested against field observations not used for model selection and calibration, the best model's mean absolute error (MAE) for daily predictions was 1.07 m (22%), more than 50% lower than for other publicly available Baltic Sea Secchi depth maps. The MAE for predicting monthly averages was 0.86 m (15%). Thus, the proposed model selection process was able to find a regional model with good prediction accuracy. It could be useful to find predictive models for environmental variables other than Secchi depth, using data from other satellite sensors, and for other regions where non-standard remote sensing models are needed for prediction and mapping. Annual and monthly mean Secchi depth maps for 2003-2012 come with this paper as Supplementary materials.

Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Climate sensitivity across marine domains of life: limits to evolutionary adaptation shape species interactions, supplementary material

Organisms in all domains, Archaea, Bacteria, and Eukarya will respond to climate change with differential vulnerabilities resulting in shifts in species distribution, coexistence, and interactions. The identification of unifying principles of organism functioning across all domains would facilitate a cause and effect understanding of such changes and their implications for ecosystem shifts. For example, the functional specialization of all organisms in limited temperature ranges leads us to ask for unifying functional reasons. Organisms also specialize in either anoxic or various oxygen ranges, with animals and plants depending on high oxygen levels. Here, we identify thermal ranges, heat limits of growth, and critically low (hypoxic) oxygen concentrations as proxies of tolerance in a meta-analysis of data available for marine organisms, with special reference to domain-specific limits. For an explanation of the patterns and differences observed, we define and quantify a proxy for organismic complexity across species from all domains. Rising complexity causes heat (and hypoxia) tolerances to decrease from Archaea to Bacteria to uni- and then multicellular Eukarya. Within and across domains, taxon-specific tolerance limits likely reflect ultimate evolutionary limits of its species to acclimatization and adaptation. We hypothesize that rising taxon-specific complexities in structure and function constrain organisms to narrower environmental ranges. Low complexity as in Archaea and some Bacteria provide life options in extreme environments. In the warmest oceans, temperature maxima reach and will surpass the permanent limits to the existence of multicellular animals, plants and unicellular phytoplankter. Smaller, less complex unicellular Eukarya, Bacteria, and Archaea will thus benefit and predominate even more in a future, warmer, and hypoxic ocean.

Sensitivity to ocean acidification parallels natural pCO2 gradients experienced by Arctic copepods under winter sea ice

The Arctic Ocean is a bellwether for ocean acidification, yet few direct Arctic studies have been carried out and limited observations exist, especially in winter. We present unique under-ice physicochemical data showing the persistence of a mid water column area of high CO2 and low pH through late winter, Zooplankton data demonstrating that the dominant copepod species are distributed across these different physicochemical conditions, and empirical data demonstrating that these copepods show sensitivity to pCO2 that parallels the range of natural pCO2 they experience through their daily vertical migration behavior. Our data, collected as part of the Catlin Arctic Survey, provide unique insight into the link between environmental variability, behavior, and an organism's physiological tolerance to CO2 in key Arctic biota.

Relative sensitivity of soft-bottom intertidal macrofauna to increased CO2 and experimental stress

We used a controlled CO2 perturbation experiment to test hypotheses about changes in diversity, composition and structure of soft-bottom intertidal macrobenthic assemblages, under realistic and locally relevant scenarios of seawater acidification. Patches of undisturbed sediment were collected from 2 types of intertidal sedimentary habitat in the Ria Formosa coastal lagoon (South Portugal) and exposed to 2 levels of seawater acidification (pH reduced by 0.3 and 0.6 units) and 1 unmanipulated (control) level. After 75 d the assemblages differed significantly between the 2 types of sediment and between field controls and the ex situ treatments, but not among the 3 pH levels tested. The naturally high values of total alkalinity buffered seawater from the changes imposed on carbonate chemistry and may have contributed to offsetting acidification at the local scale. Observed differences on biota were strongly related to the organic matter content and grain-size of the sediments, particularly to the fractions of medium and coarse sand. Soft-bottom intertidal macrofauna was significantly affected by the stress of being held in an artificial environment, but not by CO2-induced seawater acidification. Given the previously observed variations in the sensitivities of marine organisms to seawater acidification, direct extrapolations of the present findings to different regions or other types of assemblages do not seem advisable. However, the contribution of ex situ studies to the assessment of ecosystem-level responses to environmental disturbances could generally be improved by incorporating adequate field controls in the experimental design.