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.

Circumpolar dataset of Soil Organic Carbon north of treeline derived from ENVISAT ASAR GM, link to GeoTIFF

A new approach for the estimation of soil organic carbon (SOC) pools north of the tree line has been developed based on synthetic aperture radar (SAR; ENVISAT Advanced SAR Global Monitoring mode) data. SOC values are directly determined from backscatter values instead of upscaling using land cover or soil classes. The multi-mode capability of SAR allows application across scales. It can be shown that measurements in C band under frozen conditions represent vegetation and surface structure properties which relate to soil properties, specifically SOC. It is estimated that at least 29 Pg C is stored in the upper 30 cm of soils north of the tree line. This is approximately 25 % less than stocks derived from the soil-map-based Northern Circumpolar Soil Carbon Database (NCSCD). The total stored carbon is underestimated since the established empirical relationship is not valid for peatlands or strongly cryoturbated soils. The approach does, however, provide the first spatially consistent account of soil organic carbon across the Arctic. Furthermore, it could be shown that values obtained from 1 km resolution SAR correspond to accounts based on a high spatial resolution (2 m) land cover map over a study area of about 7 × 7 km in NE Siberia. The approach can be also potentially transferred to medium-resolution C-band SAR data such as ENVISAT ASAR Wide Swath with ~120 m resolution but it is in general limited to regions without woody vegetation. Global Monitoring-mode-derived SOC increases with unfrozen period length. This indicates the importance of this parameter for modelling of the spatial distribution of soil organic carbon storage.

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.

Ice and snow measurements during Nathaniel B. Palmer cruise NBP0709 (SIMBA experiment)

The Sea Ice Mass Balance in the Antarctic (SIMBA) experiment was conducted from the RVIB N.B. Palmer in September and October 2007 in the Bellingshausen Sea in an area recently experiencing considerable changes in both climate and sea ice cover. Snow and ice properties were observed at 3 short-term stations and a 27-day drift station (Ice Station Belgica, ISB) during the winter-spring transition. Repeat measurements were performed on sea ice and snow cover at 5 ISB sites, each having different physical characteristics, with mean ice (snow) thicknesses varying from 0.6 m (0.1 m) to 2.3 m (0.7 m). Ice cores retrieved every five days from 2 sites and measured for physical, biological, and chemical properties. Three ice mass-balance buoys (IMBs) provided continuous records of snow and ice thickness and temperature. Meteorological conditions changed from warm fronts with high winds and precipitation followed by cold and calm periods through four cycles during ISB. The snow cover regulated temperature flux and controlled the physical regime in which sea ice morphology changed. Level thin ice areas had little snow accumulation and experienced greater thermal fluctuations resulting in brine salinity and volume changes, and winter maximum thermodynamic growth of ~0.6 m in this region. Flooding and snow-ice formation occurred during cold spells in ice and snow of intermediate thickness. In contrast, little snow-ice formed in flooded areas with thicker ice and snow cover, instead nearly isothermal, highly permeable ice persisted. In spring, short-lived cold air episodes did not effectively penetrate the sea ice nor overcome the effect of ocean heat flux, thus favoring net ice thinning from bottom melt over ice thickening from snow-ice growth, in all cases. These warm ice conditions were consistent with regional remote sensing observations of earlier ice breakup and a shorter sea ice season, more recently observed in the Bellingshausen Sea.

Reconstructed glaciation limits in Eurasia (digitized from map Fig.1)

Based on field investigations in northern Russia and interpretation of offshore seismic data, we have made a preliminary reconstruction of the maximum ice-sheet extent in the Barents and Kara Sea region during the Early/Middle Weichselian and the Late Weichselian. Our investigations indicate that the Barents and Kara ice sheets attained their maximum Weichselian positions in northern Russia prior to 50 000 yr BP, whereas the northeastern flank of the Scandinavian Ice Sheet advanced to a maximum position shortly after 17000 calendar years ago. During the Late Weichselian (25 000-10000 yr BP), much of the Russian Arctic remained ice-free. According to our reconstruction, the extent of the ice sheets in the Barents and Kara Sea region during the Late Weichselian glacial maximum was less than half that of the maximum model which, up to now, has been widely used as a boundary condition for testing and refining General Circulation Models (GCMs). Preliminary numerical-modelling experiments predict Late Weichselian ice sheets which are larger than the ice extent implied for the Kara Sea region from dated geological evidence, suggesting very low precipitation.

Fallout and concentration of aerosol-borne radioactivity over the Atlantic (Appendix)

Measurements of atmospheric radioactivity attached to aerosols are described. Fallout was collected in a vessel of large area. Emphasis was on separation of "wet" and "dry" samples. For strontium 90 a ratio of "wet" to "dry" fallout of 5:1 has been found independent of latitude. The total fallout was smaller than comparable values from continents because of very small amounts of rainfall in the equatorial zone. In order to achieve consistency in the global balance a better knowledge not only of radioactivity but also of precipitation over the ocean is required. Fallout of Ra-D clearly shows the ITC as a barrier for the latitudinal movement of near sea-surface air masses. The concentration of short-lived emanation daughters shows large variations according to varying geographic conditions. A variation with time could not be explained. The specific activity of long-lived radioactive substances shows the expected effect of the ITC as well as a seasonal diminuation of average concentration, similar to that measured at Heidelberg.

Reconstruction of the provenance and distribution of terrigenous sediments off NE South America over the last 30ka

We investigate the redistribution of terrigenous materials in the northeastern (NE) South American continental margin during slowdown events of the Atlantic Meridional Overturning Circulation (AMOC). The compilation of stratigraphic data from 108 marine sediment cores collected across the western tropical Atlantic shows an extreme rise in sedimentation rates off the Parnaíba River mouth (about 2°S) during Heinrich Stadial 1 (HS1, 18-15 ka). Sediment core GeoB16206-1, raised offshore the Parnaíba River mouth, documents relatively constant 143Nd/144Nd values (expressed as epsilonNd(0)) throughout the last 30 ka. Whereas the homogeneous epsilonNd(0) data support the input of fluvial sediments by the Parnaíba River from the same source area directly onshore, the increases in Fe/Ca, Al/Si and Rb/Sr during HS1 indicate a marked intensification of fluvial erosion in the Parnaíba River drainage basin. In contrast, the epsilonNd(0) values from sediment core GeoB16224-1 collected off French Guiana (about 7°N) suggest Amazon-sourced materials within the last 30 ka. We attribute the extremely high volume of terrigenous sediments deposited offshore the Parnaíba River mouth during HS1 to (i) an enhanced precipitation in the catchment region and (ii) a reduced North Brazil Current, which are both associated with a weakened AMOC.

Seagrass time-series analysis products, in Moreton Bay, Australia, with links to GeoTIFFs

The spatial and temporal dynamics of seagrasses have been well studied at the leaf to patch scales, however, the link to large spatial extent landscape and population dynamics is still unresolved in seagrass ecology. Traditional remote sensing approaches have lacked the temporal resolution and consistency to appropriately address this issue. This study uses two high temporal resolution time-series of thematic seagrass cover maps to examine the spatial and temporal dynamics of seagrass at both an inter- and intra-annual time scales, one of the first globally to do so at this scale. Previous work by the authors developed an object-based approach to map seagrass cover level distribution from a long term archive of Landsat TM and ETM+ images on the Eastern Banks (~200 km**2), Moreton Bay, Australia. In this work a range of trend and time-series analysis methods are demonstrated for a time-series of 23 annual maps from 1988 to 2010 and a time-series of 16 monthly maps during 2008-2010. Significant new insight was presented regarding the inter- and intra-annual dynamics of seagrass persistence over time, seagrass cover level variability, seagrass cover level trajectory, and change in area of seagrass and cover levels over time. Overall we found that there was no significant decline in total seagrass area on the Eastern Banks, but there was a significant decline in seagrass cover level condition. A case study of two smaller communities within the Eastern Banks that experienced a decline in both overall seagrass area and condition are examined in detail, highlighting possible differences in environmental and process drivers. We demonstrate how trend and time-series analysis enabled seagrass distribution to be appropriately assessed in context of its spatial and temporal history and provides the ability to not only quantify change, but also describe the type of change. We also demonstrate the potential use of time-series analysis products to investigate seagrass growth and decline as well as the processes that drive it. This study demonstrates clear benefits over traditional seagrass mapping and monitoring approaches, and provides a proof of concept for the use of trend and time-series analysis of remotely sensed seagrass products to benefit current endeavours in seagrass ecology.

Calibration and application of B/Ca, Cd/Ca, and d11B in Neogloboquadrina pachyderma (sinistral)

The North Atlantic and Norwegian Sea are prominent sinks of atmospheric CO2 today, but their roles in the past remain poorly constrained. In this study, we attempt to use B/Ca and d11B ratios in the planktonic foraminifera Neogloboquadrina pachyderma (sinistral variety) to reconstruct subsurface water pH and pCO2 changes in the polar North Atlantic during the last deglaciation. Comparison of core-top results with nearby hydrographic data shows that B/Ca in N. pachyderma (s) is mainly controlled by seawater [B(OH)4]?/[HCO3]? with a roughly constant partition coefficient (KD =([B/Ca]of CaCO3)/([B(OH)4]-/[HCO3]-)of seawater) of 1.48 ± 0.15 * 10**-3 (2sigma), and d11B in this species is offset below d11B of the borate in seawater by 3.38 ± 0.71 per mil (2sigma). These values represent our best estimates with the sparse available hydrographic data close to our core-tops. More culturing and sediment trap work is needed to improve our understanding of boron incorporation into N. pachyderma (s). Application of a constant KD of 1.48 * 10**-3 to high resolution N. pachyderma (s) B/Ca records from two adjacent cores off Iceland shows that subsurface pCO2 at the habitat depth of N. pachyderma (s) (~50 m) generally followed the atmospheric CO2 trend but with negative offsets of ~10-50 ppmv during 19-10 ka. These B/Ca-based reconstructions are supported by independent estimates from low-resolution d11B measurements in the same cores. We also calibrate and apply Cd/Ca in N. pachyderma (s) to reconstruct nutrient levels for the same down cores. Like today's North Atlantic, past subsurface pCO2 variability off Iceland was significantly correlated with nutrient changes that might be linked to surface nutrient utilization and mixing within the upper water column. Because surface pCO2 (at 0 m water depth) is always lower than at deeper depths and if the application of a constant KD is valid, our results suggest that the polar North Atlantic has remained a CO2 sink during the calcification seasons of N. pachyderma (s) over the last deglaciation.

Concentration and composition of aerosols over the North Atlantic in August 1995

Data on composition of aerosols are considered. Investigations include electron microscopy, grain size, mineralogical and chemical analyses. Samples of aerosols were collected Cruise 37 of R/V Akademik Mstislav Keldysh along a transatlantic section along 40°-50°N. Variability of concentrations and composition of aerosols with distance from continents is shown: concentrations of aerosols decrease by factor of ten and more. Significant portion of mineral components in aerosols collected over the continent is replaced by organic matter due to mechanical differentiation during transportation. Such anthropogenic components as soot, ash, and combustion spheres were detected in all samples. North African dust was found in one sample in the western part of the section.