Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.

Gotland Deep (Baltic Sea) independent age model based on Palaeomagnetic Secular Variation (PSV) and lead (Pb) deposition profiles

Dating of sediment cores from the Baltic Sea has proven to be difficult due to uncertainties surrounding the 14C reservoir age and a scarcity of macrofossils suitable for dating. Here we present the results of multiple dating methods carried out on cores in the Gotland Deep area of the Baltic Sea. Particular emphasis is placed on the Littorina stage (8 ka ago to the present) of the Baltic Sea and possible changes in the 14C reservoir age of our dated samples. Three geochronological methods are used. Firstly, palaeomagnetic secular variations (PSV) are reconstructed, whereby ages are transferred to PSV features through comparison with varved lake sediment based PSV records. Secondly, lead (Pb) content and stable isotope analysis are used to identify past peaks in anthropogenic atmospheric Pb pollution. Lastly, 14C determinations were carried out on benthic foraminifera (Elphidium spec.) samples from the brackish Littorina stage of the Baltic Sea. Determinations carried out on smaller samples (as low as 4 µg C) employed an experimental, state-of-the-art method involving the direct measurement of CO2 from samples by a gas ion source without the need for a graphitisation step - the first time this method has been performed on foraminifera in an applied study. The PSV chronology, based on the uppermost Littorina stage sediments, produced ten age constraints between 6.29 and 1.29 cal ka BP, and the Pb depositional analysis produced two age constraints associated with the Medieval pollution peak. Analysis of PSV data shows that adequate directional data can be derived from both the present Littorina saline phase muds and Baltic Ice Lake stage varved glacial sediments. Ferrimagnetic iron sulphides, most likely authigenic greigite (Fe3S4), present in the intermediate Ancylus Lake freshwater stage sediments acquire a gyroremanent magnetisation during static alternating field (AF) demagnetisation, preventing the identification of a primary natural remanent magnetisation for these sediments. An inferred marine reservoir age offset (deltaR) is calculated by comparing the foraminifera 14C determinations to a PSV & Pb age model. This deltaR is found to trend towards younger values upwards in the core, possibly due to a gradual change in hydrographic conditions brought about by a reduction in marine water exchange from the open sea due to continued isostatic rebound.

Raster grids detailing the benthic environment of South Georgia

This dataset contains raster grids in GeoTIFF format describing the benthic environment of South Georgia. The data include topographic layers that are directly calculated from a bathymetry grid (Slope, Aspect, Roughness, Slope, Terrain Ruggedness Index, Topographic Position Index). A benthic classification of the area is included, based on topographic layers. Also included are sea-bed environmental layers that are interpolated from global three dimensional grids (Alkalinity, Apparent Oxygen Utilisation, Omega Aragonite, Omega Calcite, Dissolved Oxygen, Nitrate, pH, Phosphate, Salinity, Silicate, Temperature, and Total CO2). These layers were used to construct a habitat suitability model for Octocorallia. The geographic extent is 43°57'56.65"W - 33°45'38.19"W and 52°47'29.50"S - 56° 9'11.03"S. The spatial resolution is 150m x 150m (except for benthic classification wihch is 450m x 450m). The map projection is EPSG:3762.

Global habitat suitability of deep-water Antipatharia

Antipatharia are a diverse group of corals with many species found in deep water. Many Antipatharia are habitat for associates, have extreme longevity and some species can occur beyond 8500 m depth. As they are major constituents of 'coral gardens', which are Vulnerable Marine Ecosystems (VMEs), knowledge of their distribution and environmental requirements is an important pre-requisite for informed conservation planning particularly where the expense and difficulty of deep-sea sampling prohibits comprehensive surveys. This study uses a global database of Antipatharia distribution data to perform habitat suitability modelling using the Maxent methodology to estimate the global extent of black coral habitat suitability. The model of habitat suitability is driven by temperature but there is notable influence from other variables of topography, surface productivity and oxygen levels. This model can be used to predict areas of suitable habitat, which can be useful for conservation planning. The global distribution of Antipatharia habitat suitability shows a marked contrast with the distribution of specimen observations, indicating that many potentially suitable areas have not been sampled, and that sampling effort has been disproportionate to shallow, accessible areas inside marine protected areas (MPAs). Although 25% of Antipatharia observations are located in MPAs, only 7-8% of predicted suitable habitat is protected, which is short of the Convention on Biological Diversity target to protect 10% of ocean habitats by 2020.

(Table 1) Krill (Euphausia superba) biomass estimates at South Georgia Island from 2001-2005

The South Georgia region supports a large biomass of krill that is subject to high interannual variability. The apparent lack of a locally self-maintaining krill population at South Georgia means that understanding the mechanism underlying these observed population characteristics is essential to successful ecosystem-based management of krill fishery in the region. Krill acoustic-density data from surveys conducted in the early, middle and late period of the summers of 2001 to 2005, together with krill population size structure over the same period from predator diet data, were used with a krill population dynamics model to evaluate potential mechanisms behind the observed changes in krill biomass. Krill abundance was highest during the middle of the summer in 3 years and in the late period in 2 years; in the latter there was evidence that krill recruitment was delayed by several months. A model scenario that included empirically derived estimates of both the magnitude and timing of recruitment in each year showed the greatest correlation with the acoustic series. The results are consistent with a krill population with allochthonous recruitment entering a retained adult population; i.e. oceanic transport of adult krill does not appear to be the major factor determining the dynamics of the adult population. The results highlight the importance of the timing of recruitment, especially where this could introduce a mismatch between the peak of krill abundance and the peak demand from predators, which may exacerbate the effects of changes in krill populations arising from commercial harvesting and/or climate change.

Interview forms for integrating education for conservation in marine protected areas management in Colombia

Environmental Education (EE) is a key component in any marine protected area management. However, its visibility and action plans are still poorly developed and structured as a clear element in management procedures. The objective of this study is to contribute with a methodological route that integrates EE to the existing model of management planning and strategies, taking the Colombian National Natural Parks System as a case study. The creation of the route is proposed as a participatory research with different stakeholders in order to respond to the specific conservation needs and goals for the National Parks System. The EE national diagnosis has shown that its integration within the parks management structure is a first priority need, being a converging result on the two case studies on National Parks from the Pacific Coast of Colombia. The diagnosis also demonstrates that communication, participation, training and evaluation have to be reinforced, linking the community and stakeholders involved in the park management to the whole EE process. The proposed methodology route has been agreed upon by the National Parks staff and incorporates advice and recommendations from different stakeholders, in order to better include the park users. This step will help us to advance toward sustainable management in marine and coastal protected areas elsewhere, taking into account not only the biological but also the social-cultural prism. The main challenges in the management and conservation of coastal and marine ecosystems today are discussed.

Occurrence (presence-only) of living Lophelia pertusa reefs in the Irish continental margin

The present data set was used as a training set for a Habitat Suitability Model. It contains occurrence (presence-only) of living Lophelia pertusa reefs in the Irish continental margin, which were assembled from databases, cruise reports and publications. A total of 4423 records were inspected and quality assessed to ensure that they (1) represented confirmed living L. pertusa reefs (so excluding 2900 records of dead and isolated coral colony records); (2) were derived from sampling equipment that allows for accurate (<200 m) geo-referencing (so excluding 620 records derived mainly from trawling and dredging activities); and (3) were not duplicated. A total of 245 occurrences were retained for the analysis. Coral observations are highly clustered in regions targeted by research expeditions, which might lead to falsely inflated model evaluation measures (Veloz, 2009). Therefore, we coarsened the distribution data by deleting all but one record within grid cells of 0.02° resolution (Davies & Guinotte 2011). The remaining 53 points were subject to a spatial cross-validation process: a random presence point was chosen, grouped with its 12 closest neighbour presence points based on Euclidean distance and withheld from model training. This process was repeated for all records, resulting in 53 replicates of spatially non-overlapping sets of test (n=13) and training (n=40) data. The final 53 occurrence records were used for model training.