Probability inundation map for a scenario combining storm surge and sea-level rise

Coastal managers require reliable spatial data on the extent and timing of potential coastal inundation, particularly in a changing climate. Most sea level rise (SLR) vulnerability assessments are undertaken using the easily implemented bathtub approach, where areas adjacent to the sea and below a given elevation are mapped using a deterministic line dividing potentially inundated from dry areas. This method only requires elevation data usually in the form of a digital elevation model (DEM). However, inherent errors in the DEM and spatial analysis of the bathtub model propagate into the inundation mapping. The aim of this study was to assess the impacts of spatially variable and spatially correlated elevation errors in high-spatial resolution DEMs for mapping coastal inundation. Elevation errors were best modelled using regression-kriging. This geostatistical model takes the spatial correlation in elevation errors into account, which has a significant impact on analyses that include spatial interactions, such as inundation modelling. The spatial variability of elevation errors was partially explained by land cover and terrain variables. Elevation errors were simulated using sequential Gaussian simulation, a Monte Carlo probabilistic approach. 1,000 error simulations were added to the original DEM and reclassified using a hydrologically correct bathtub method. The probability of inundation to a scenario combining a 1 in 100 year storm event over a 1 m SLR was calculated by counting the proportion of times from the 1,000 simulations that a location was inundated. This probabilistic approach can be used in a risk-aversive decision making process by planning for scenarios with different probabilities of occurrence. For example, results showed that when considering a 1% probability exceedance, the inundated area was approximately 11% larger than mapped using the deterministic bathtub approach. The probabilistic approach provides visually intuitive maps that convey uncertainties inherent to spatial data and analysis.

The International Bathymetric Chart of the Southern Ocean (IBCSO) Version 1.0

The International Bathymetric Chart of the Southern Ocean (IBCSO) Version 1.0 is a new digital bathymetric model (DBM) portraying the seafloor of the circum-Antarctic waters south of 60° S. IBCSO is a regional mapping project of the General Bathymetric Chart of the Oceans (GEBCO). IBCSO Version 1.0 DBM has been compiled from all available bathymetric data collectively gathered by more than 30 institutions from 15 countries. These data include multibeam and single beam echo soundings, digitized depths from nautical charts, regional bathymetric gridded compilations, and predicted bathymetry. Specific gridding techniques were applied to compile the DBM from the bathymetric data of different origin, spatial distribution, resolution, and quality. The IBCSO Version 1.0 DBM has a resolution of 500 x 500 m, based on a polar stereographic projection, and is publicly available together with a digital chart for printing from the project website (http://www.ibcso.org) and from the two data sets shown at the bottom of this page.

Plant Functional Type (PFT) map over Siberia derived from GlobCover 2005 dataset, link to dataset in NetCDF format

High-latitude ecosystems play an important role in the global carbon cycle and in regulating the climate system and are presently undergoing rapid environmental change. Accurate land cover data sets are required to both document these changes as well as to provide land-surface information for benchmarking and initializing Earth system models. Earth system models also require specific land cover classification systems based on plant functional types (PFTs), rather than species or ecosystems, and so post-processing of existing land cover data is often required. This study compares over Siberia, multiple land cover data sets against one another and with auxiliary data to identify key uncertainties that contribute to variability in PFT classifications that would introduce errors in Earth system modeling. Land cover classification systems from GLC 2000, GlobCover 2005 and 2009, and MODIS collections 5 and 5.1 are first aggregated to a common legend, and then compared to high-resolution land cover classification systems, vegetation continuous fields (MODIS VCFs) and satellite-derived tree heights (to discriminate against sparse, shrub, and forest vegetation). The GlobCover data set, with a lower threshold for tree cover and taller tree heights and a better spatial resolution, tends to have better distributions of tree cover compared to high-resolution data. It has therefore been chosen to build new PFT maps for the ORCHIDEE land surface model at 1 km scale. Compared to the original PFT data set, the new PFT maps based on GlobCover 2005 and an updated cross-walking approach mainly differ in the characterization of forests and degree of tree cover. The partition of grasslands and bare soils now appears more realistic compared with ground truth data. This new vegetation map provides a framework for further development of new PFTs in the ORCHIDEE model like shrubs, lichens and mosses, to represent the water and carbon cycles in northern latitudes better. Updated land cover data sets are critical for improving and maintaining the relevance of Earth system models for assessing climate and human impacts on biogeochemistry and biophysics.

Satellite image map of King George Island, Antarctica

A new topographic database for King George Island, one of the most visited areas in Antarctica, is presented. Data from differential GPS surveys, gained during the summers 1997/98 and 1999/2000, were combined with up to date coastlines from a SPOT satellite image mosaic, and topographic information from maps as well as from the Antarctic Digital Database. A digital terrain model (DTM) was generated using ARC/INFO GIS. From contour lines derived from the DTM and the satellite image mosaic a satellite image map was assembled. Extensive information on data accuracy, the database as well as on the criteria applied to select place names is given in the multilingual map. A lack of accurate topographic information in the eastern part of the island was identified. It was concluded that additional topographic surveying or radar interferometry should be conducted to improve the data quality in this area. In three case studies, the potential applications of the improved topographic database are demonstrated. The first two examples comprise the verification of glacier velocities and the study of glacier retreat from the various input data-sets as well as the use of the DTM for climatological modelling. The last case study focuses on the use of the new digital database as a basic GIS (Geographic Information System) layer for environmental monitoring and management on King George Island.

Seawater carbonate chemistry and biological processes of coral Acropora muricata during experiments and observations in La Saline fringing reef, La Reunion Island, western Indian Ocean, 2011

The effect of decreasing aragonite saturation state (Omega Arag) of seawater (elevated pCO2) on calcification rates of Acropora muricata was studied using nubbins prepared from parent colonies located at two sites of La Saline reef (La Réunion Island, western Indian Ocean): a back-reef site (BR) affected by nutrient-enriched groundwater discharge (mainly nitrate), and a reef flat site (RF) with low terrigenous inputs. Protein and chlorophyll a content of the nubbins, as well as zooxanthellae abundance, were lower at RF than BR. Nubbins were incubated at ~27°C over 2 h under sunlight, in filtered seawater manipulated to get differing initial pCO2 (1,440-340 µatm), Omega Arag (1.4-4.0), and dissolved inorganic carbon (DIC) concentrations (2,100-1,850 µmol/kg). Increasing DIC concentrations at constant total alkalinity (AT) resulted in a decrease in Omega Arag and an increase in pCO2. AT at the beginning of the incubations was kept at a natural level of 2,193 ± 6 µmol/kg (mean ± SD). Net photosynthesis (NP) and calcification were calculated from changes in pH and AT during the incubations. Calcification decrease in response to doubling pCO2 relative to preindustrial level was 22% for RF nubbins. When normalized to surface area of the nubbins, (1) NP and calcification were higher at BR than RF, (2) NP increased in high pCO2 treatments at BR compared to low pCO2 treatments, and (3) calcification was not related to Omega Arag at BR. When normalized to NP, calcification was linearly related to Omega Arag at both sites, and the slopes of the relationships were not significantly different. The increase in NP at BR in the high pCO2 treatments may have increased calcification and thus masked the negative effect of low Omega Arag on calcification. Removing the effect of NP variations at BR showed that calcification declined in a similar manner with decreased Omega Arag (increased pCO2) whatever the nutrient loading.

Seawater carbonate chemistry and protein and chlorophyll a content per nubbin of Acropora muricata during observations of La Saline fringing reef, La Reunion Island, western Indian Ocean, 2011

The effect of decreasing aragonite saturation state (Omega Arag) of seawater (elevated pCO2) on calcification rates of Acropora muricata was studied using nubbins prepared from parent colonies located at two sites of La Saline reef (La Réunion Island, western Indian Ocean): a back-reef site (BR) affected by nutrient-enriched groundwater discharge (mainly nitrate), and a reef flat site (RF) with low terrigenous inputs. Protein and chlorophyll a content of the nubbins, as well as zooxanthellae abundance, were lower at RF than BR. Nubbins were incubated at ~27°C over 2 h under sunlight, in filtered seawater manipulated to get differing initial pCO2 (1,440-340 µatm), Omega Arag (1.4-4.0), and dissolved inorganic carbon (DIC) concentrations (2,100-1,850 µmol/kg). Increasing DIC concentrations at constant total alkalinity (AT) resulted in a decrease in Omega Arag and an increase in pCO2. AT at the beginning of the incubations was kept at a natural level of 2,193 ± 6 µmol/kg (mean ± SD). Net photosynthesis (NP) and calcification were calculated from changes in pH and AT during the incubations. Calcification decrease in response to doubling pCO2 relative to preindustrial level was 22% for RF nubbins. When normalized to surface area of the nubbins, (1) NP and calcification were higher at BR than RF, (2) NP increased in high pCO2 treatments at BR compared to low pCO2 treatments, and (3) calcification was not related to Omega Arag at BR. When normalized to NP, calcification was linearly related to Omega Arag at both sites, and the slopes of the relationships were not significantly different. The increase in NP at BR in the high pCO2 treatments may have increased calcification and thus masked the negative effect of low Omega Arag on calcification. Removing the effect of NP variations at BR showed that calcification declined in a similar manner with decreased Omega Arag (increased pCO2) whatever the nutrient loading.

AWI Bathymetric Chart of the Fram Strait (BCFS) (Scale 1:100,000)

Based on data from R/V Polarstern multibeam sonar surveys between 1984 and 1997 a high resolution bathymetry has been generated for the central Fram Strait. The area ensonified covers approx. 36,500 sqkm between 78°N - 80°N and 0°E - 7.5°E. Basic outcome of the investigation is a Digital Terrain Model (DTM) with 100 m grid spacing which was utilized for contouring and generation of a new series of bathymetric charts at a scale of 1:100,000, the AWI Bathymetric Chart of the Fram Strait (AWI BCFS).

Volume changes of Jacobshavn Isbrae between 1985, 1997 and 2007

Following three decades of relative stability, Jakobshavn Isbrae, West Greenland, underwent dramatic thinning, retreat and speed-up starting in 1998. To assess the amount of ice loss, we analyzed 1985 aerial photos and derived a 40 m grid digital elevation model (DEM). We also obtained a 2007 40 m grid SPOT DEM covering the same region. Comparison of the two DEMs over an area of ~4000 km**2 revealed a total ice loss of 160 ± 4 km**3, with 107 ± 0.2 km**3 in grounded regions (0.27 mm eustatic sea-level rise) and 53 ± 4 km**3 from the disintegration of the floating tongue. Comparison of the DEMs with 1997 NASA Airborne Topographic Mapper data indicates that this ice loss essentially occurred after 1997, with +0.7 ± 5.6 km**3 between 1985 and 1997 and -160 ± 7 km**3 between 1997 and 2007. The latter is equivalent to an average specific mass balance of -3.7 ± 0.2 m/a over the study area. Previously reported thickening of the main glacier during the early 1990s was accompanied by similar-magnitude thinning outside the areas of fast flow, indicating that the land-based ice continued reacting to longer-term climate forcing.

Geological map of Heimefrontfjella, Dronning Maud Land, Antarctica (Scale 1:250,000)

The geological overview map was compiled from 15 geological maps (1 : 25,000) and is based on Jacobs et al. 1996. The topographic basemaps were adapted from unpublished 1:250,000 provisional topographic maps, Institut f. Angewandte Geodäsie, Frankfurt, 1983. Part of the contour lines are from Radarsat (Liu et al. 2001).

Mapping of the Eltanin impact area in the South-East Pacific

Two Polarstern expeditions were conducted in 1995 (ANT-XII/4) and 2001 (ANT-XVIII/5a) to the Bellingshausen Sea and Amundsen Sea and the suspected Eltanin meteorite impact in the SE-Pacific. A survey of the sediment distribution and its acoustic structure along the cruise track was performed. The seafloor topography was sampled using the multibeam sonar system Hydrosweep DS2 which operates on a frequency of 15.5 kHz. The resulting AWI Bathymetric Chart of the Eltanin Meteorite Impact Area is based on a Digital Terrain Model of this area. The mapping was performed using ArcGIS. The Eltanin impact area which covers the 4.100 m high Freden Seamount is visualized by one overview sheet of the scale 1:200,000 and four 1:100,000 subsheets.

AWI Bathymetric Chart of the Gakkel Ridge (AWI BCGR) (Scale 1:150,000)

The 1 : 1,500,000 AWI Bathymetric Chart of the Gakkel Ridge (AWI BCGR) has been developed from multibeam data measured during the Arctic Mid-Ocean Ridge Expedition in 2001 (AMORE 2001, ARK-XVII/2). This expedition was conducted to investigate the Gakkel Ridge in the Arctic Ocean and was carried out by the icebreaking research vessels RV Polarstern and USCGC Healy. Polarstern is equipped with the multibeam sonar system Hydrosweep DS-2, whereas Healy carries Seabeam 2112. During the expedition an area of 8890 km length and 18 - 46 km width, situated between 82°N/8°W and 87°N/75°E, was surveyed simultaneously by both vessels. Water depths ranged from 566 to 5673 meters. Dense sea ice cover derogated the sonar measurements and decreased data quality. Data errors were corrected in an extensive post-processing. The data of two different sonar systems had to be consolidated in order to derive a high resolution bathymetry of the Gakkel Ridge. Final result was a digital terrain model (DTM) with a grid spacing of 100 meters, which was utilized for generating the map series AWI Bathymetric Chart of the Gakkel Ridge, consisting of ten map sheets.

King George Island TanDEM-X DEM, link to GeoTIFF

We have generated a new digital elevation model for entire King George Island, Antarctica, using summer TanDEM-X bistatic SAR satellite data. The data was processed using differential SAR interferometry with an older DEM as reference. 4 TanDEM-X scenes from January 2012 were used as input. The new DEM was referenced to and validated against DGPS measurements. Height values are given in reference to ellipsoid (WGS84).

Seawater carbonate chemistry, photosynthesis and calcification rate during experiments with coral Acropora muricata, 2011

The effect of decreasing aragonite saturation state (Omega Arag) of seawater (elevated pCO2) on calcification rates of Acropora muricata was studied using nubbins prepared from parent colonies located at two sites of La Saline reef (La Réunion Island, western Indian Ocean): a back-reef site (BR) affected by nutrient-enriched groundwater discharge (mainly nitrate), and a reef flat site (RF) with low terrigenous inputs. Protein and chlorophyll a content of the nubbins, as well as zooxanthellae abundance, were lower at RF than BR. Nubbins were incubated at ~27°C over 2 h under sunlight, in filtered seawater manipulated to get differing initial pCO2 (1,440-340 µatm), Omega Arag (1.4-4.0), and dissolved inorganic carbon (DIC) concentrations (2,100-1,850 µmol kg-1). Increasing DIC concentrations at constant total alkalinity (AT) resulted in a decrease in Omega Arag and an increase in pCO2. AT at the beginning of the incubations was kept at a natural level of 2,193 +- 6 µmol kg-1 (mean +- SD). Net photosynthesis (NP) and calcification were calculated from changes in pH and AT during the incubations. Calcification decrease in response to doubling pCO2 relative to preindustrial level was 22% for RF nubbins. When normalized to surface area of the nubbins, (1) NP and calcification were higher at BR than RF, (2) NP increased in high pCO2 treatments at BR compared to low pCO2 treatments, and (3) calcification was not related to Omega Arag at BR. When normalized to NP, calcification was linearly related to Omega Arag at both sites, and the slopes of the relationships were not significantly different. The increase in NP at BR in the high pCO2 treatments may have increased calcification and thus masked the negative effect of low Omega Arag on calcification. Removing the effect of NP variations at BR showed that calcification declined in a similar manner with decreased Omega Arag (increased pCO2) whatever the nutrient loading.

A suite of global, cross-scale topographic variables for environmental and biodiversity modeling, links to files in GeoTIFF format

Topographic variation, the spatial variation in elevation and terrain features, underpins a myriad of patterns and processes in geography and ecology and is key to understanding the variation of life on the planet. The characterization of this variation is scale-dependent, i.e. it varies with the distance over which features are assessed and with the spatial grain (grid cell resolution) of analysis. A fully standardized and global multivariate product of different terrain features has the potential to support many large-scale basic research and analytical applications, however to date, such technique is unavailable. Here we used the digital elevation model products of global 250 m GMTED and near-global 90 m SRTM to derive a suite of topographic variables: elevation, slope, aspect, eastness, northness, roughness, terrain roughness index, topographic position index, vector ruggedness measure, profile and tangential curvature, and 10 geomorphological landform classes. We aggregated each variable to 1, 5, 10, 50 and 100 km spatial grains using several aggregation approaches (median, average, minimum, maximum, standard deviation, percent cover, count, majority, Shannon Index, entropy, uniformity). While a global cross-correlation underlines the high similarity of many variables, a more detailed view in four mountain regions reveals local differences, as well as scale variations in the aggregated variables at different spatial grains. All newly-developed variables are available for download at http://www.earthenv.org and can serve as a basis for standardized hydrological, environmental and biodiversity modeling at a global extent.

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.