(Table 1) Accumulation and ice discharge of West Antarctic glaciers draining into the Amundsen Sea between 1974 and 2007

The Advanced Land Observation System (ALOS) Phased-Array Synthetic-Aperture Radar (PALSAR) is an L-band frequency (1.27 GHz) radar capable of continental-scale interferometric observations of ice sheet motion. Here, we show that PALSAR data yield excellent measurements of ice motion compared to C-band (5.6 GHz) radar data because of greater temporal coherence over snow and firn. We compare PALSAR velocities from year 2006 in Pine Island Bay, West Antarctica with those spanning years 1974 to 2007. Between 1996 and 2007, Pine Island Glacier sped up 42% and ungrounded over most of its ice plain. Smith Glacier accelerated 83% and ungrounded as well. Their largest speed up are recorded in 2007. Thwaites Glacier is not accelerating but widening with time and its eastern ice shelf doubled its speed. Total ice discharge from these glaciers increased 30% in 12 yr and the net mass loss increased 170% from 39 ± 15 Gt/yr to 105 ± 27 Gt/yr. Longer-term velocity changes suggest only a moderate loss in the 1970s. As the glaciers unground into the deeper, smoother beds inland, the mass loss from this region will grow considerably larger in years to come.

(Fig. 1D, 3E) Sea surface temperature and alkenone C37:4 for sediment core MD01-2444 during the penultimate glaciation

Past glacials can be thought of as natural experiments in which variations in boundary conditions influenced the character of climate change. However, beyond the last glacial, an integrated view of orbital- and millennial-scale changes and their relation to the record of glaciation has been lacking. Here, we present a detailed record of variations in the land-ocean system from the Portuguese margin during the penultimate glacial and place it within the framework of ice-volume changes, with particular reference to European ice-sheet dynamics. The interaction of orbital- and millennial-scale variability divides the glacial into an early part with warmer and wetter overall conditions and prominent climate oscillations, a transitional mid-part, and a late part with more subdued changes as the system entered a maximum glacial state. The most extreme event occurred in the mid-part and was associated with melting of the extensive European ice sheet and maximum discharge from the Fleuve Manche river. This led to disruption of the meridional overturning circulation, but not a major activation of the bipolar seesaw. In addition to stadial duration, magnitude of freshwater forcing, and background climate, the evidence also points to the influence of the location of freshwater discharges on the extent of interhemispheric heat transport.

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.

Seasonally resolved chemical records of six ice cores from the Dronning Maud Land, Antarctica

The ammonium, calcium, and sodium concentrations from three intermediate depth ice cores drilled in the area of Dronning Maud Land, East Antarctica, have been investigated. Since all measurements were performed by a high-resolution Continuous Flow Analysis system, for the first time seasonal signals of chemical trace species could be obtained from the interior of central Antarctica over a period of approximately 2 millennia. Although the elevation as well as the accumulation rate differ between the drilling sites, similar values were obtained by comparing mean concentrations spanning the last 900 years. However, a distinct lack of intersite correlation was found on decadal timescales. Despite a noticeable accumulation change, apparent in one core, no significant concentration change of all three species has occurred. All the measured ions show clear seasonal signals over the whole records. While the sea-salt-related component sodium peaks simultaneously with calcium, the maximum ammonium concentration occurs in the snow with a time lag of 2 months after the sea-salt peak. More than 60% of the calcium concentration can be attributed to an ocean source. Elevated sodium concentrations were found within this millennium compared to mean values of the whole records, but the spatially varying shape of the increase suggests that a possible climatic signal is biased by local deposition effects.

High-resolution ice thickness and bed topography of a land-terminating section of the Greenland Ice Sheet, with links to GeoTIFFs and text files

We present ice thickness and bed topography maps with a high spatial resolution (250-500 m) of a land-terminating section of the Greenland Ice Sheet derived from ground-based and airborne radar surveys. The data have a total area of ~12 000 km^2 and cover the whole ablation area of the outlet glaciers of Isunnguata Sermia, Russell, Leverett, Ørkendalen and Isorlersuup up to the long-term mass balance equilibrium line altitude at ~1600 m above sea level. The bed topography shows highly variable subglacial trough systems, and the trough of Isunnguata Sermia Glacier is overdeepened and reaches an elevation of ~500 m below sea level. The ice surface is smooth and only reflects the bedrock topography in a subtle way, resulting in a highly variable ice thickness. The southern part of our study area consists of higher bed elevations compared to the northern part. The compiled data sets of ground-based and airborne radar surveys cover one of the most studied regions of the Greenland Ice Sheet and can be valuable for detailed studies of ice sheet dynamics and hydrology.

EPICA Dronning Maud Land EDML ice core drilling protocol

We report on the EPICA Dronning Maud Land (East Antarctica) deep drilling operation. Starting with the scientific questions that led to the outline of the EPICA project, we introduce the setting of sister drillings at NorthGRIP and EPICA Dome C within the European ice-coring community. The progress of the drilling operation is described within the context of three parallel, deep-drilling operations, the problems that occurred and the solutions we developed. Modified procedures are described, such as the monitoring of penetration rate via cable weight rather than motor torque, and modifications to the system (e.g. closing the openings at the lower end of the outer barrel to reduce the risk of immersing the drill in highly concentrated chip suspension). Parameters of the drilling (e.g. core-break force, cutter pitch, chips balance, liquid level, core production rate and piece number) are discussed. We also review the operational mode, particularly in the context of achieved core length and piece length, which have to be optimized for drilling efficiency and core quality respectively. We conclude with recommendations addressing the design of the chip-collection openings and strictly limiting the cable-load drop with respect to the load at the start of the run.

Digital elevation model (DEM) of the ice sheet in western Dronning Maud Land, Antarctica

In this paper, a new digital elevation model (DEM) is derived for the ice sheet in western Dronning Maud Land, Antarctica. It is based on differential interferometric synthetic aperture radar (SAR) from the European Remote Sensing 1/2 (ERS-1/2) satellites, in combination with ICESat's Geoscience Laser Altimeter System (GLAS). A DEM mosaic is compiled out of 116 scenes from the ERS-1 ice phase in 1994 and the ERS-1/2 tandem mission between 1996 and 1997 with the GLAS data acquired in 2003 that served as ground control. Using three different SAR processors, uncertainties in phase stability and baseline model, resulting in height errors of up to 20 m, are exemplified. Atmospheric influences at the same order of magnitude are demonstrated, and corresponding scenes are excluded. For validation of the DEM mosaic, covering an area of about 130,000 km**2 on a 50-m grid, independent ICESat heights (2004-2007), ground-based kinematic GPS (2005), and airborne laser scanner data (ALS, 2007) are used. Excluding small areas with low phase coherence, the DEM differs in mean and standard deviation by 0.5 +/- 10.1, 1.1 +/- 6.4, and 3.1 +/- 4.0 m from ICESat, GPS, and ALS, respectively. The excluded data points may deviate by more than 50 m. In order to suppress the spatially variable noise below a 5-m threshold, 18% of the DEM area is selectively averaged to a final product at varying horizontal spatial resolution. Apart from mountainous areas, the new DEM outperforms other currently available DEMs and may serve as a benchmark for future elevation models such as from the TanDEM-X mission to spatially monitor ice sheet elevation.

Bedding dip data from core site CRP-3 in Victoria Land, Antarctica

Bedding dips in the CRP-3 drillhole were determined in three ways: (1) analysis of a dipmeter log, (2) identification of bed boundaries on borehole televiewer log images, and (3) identification of bed boundaries on digital images of the outer surfaces of oriented cores. All three methods determine both dip magnitude and downdip azimuth of bedding. Dipmeter results document variations in bedding dip throughout the logged interval (20-902 mbsf), whereas core and televiewer results are available at present only for selected depth intervals. Dipmeter data indicate that structural dip is remarkably constant, at 21° dip to azimuth 65°, throughout the Tertiary shelf section, except for the top 100 m where dips appear to be 5-10° shallower. This pattern, in conjunction with the systematically increasing dips throughout CRP-2A, suggests that the growth faulting active during CRP-2A deposition began during the final period of deposition at CRP-3. Normal faults at 260 and 539 mbsf in CRP-3 exhibit neither drag (localized dip steepening) nor significant changes in structural dip across them. Oriented core and televiewer analyses, covering a total of 200 m in the interval 400-900 mbsf, indicate bedding patterns that confirm the dipmeter results. The doleritic breccia at the base of the Tertiary section has steeper dips than overlying structural dips, possibly indicating a sedimentary dip to ENE in these fan sediments. Dip directions in the underlying Devonian Beacon sandstone are surprisingly similar to those in the overlying Tertiary section. Superimposed on the average Beacon dip of 22° to the ENE are localized tilts of up to 20°, probably caused by Tertiary fracturing and brecciation rather than original sedimentary dip variations.

Static GPS measurements in the vicinity of EDML and an elevation model for Dronning Maud Land

The knowledge of ice sheet surface topography and the location of the ice divides are essential for ice dynamic modeling. An improved digital elevation model (DEM) of Dronning Maud Land (DML), Antarctica, is presented in this paper. It is based on ground-based kinematic GPS profiles, airborne radar altimetry, and data of the airborne radio-echo sounding system, as well as spaceborne laser altimetry from NASA's Ice, Cloud and land Elevation Satellite (ICESat). The accuracy of ICESat ice sheet altimetry data in the area of investigation is discussed. The location of the ice divides is derived from aspect calculation of the topography and is verified with several velocity data derived from repeated static GPS measurements.

Per-field land use and land cover data set of the Haean catchment, South Korea

Detailed data on land use and land cover constitute important information for Earth system models, environmental monitoring and ecosystem services research. Global land cover products are evolving rapidly; however, there is still a lack of information particularly for heterogeneous agricultural landscapes. We censused land use and land cover field by field in the agricultural mosaic catchment Haean in South Korea. We recorded the land cover types with additional information on agricultural practice. In this paper we introduce the data, their collection and the post-processing protocol. Furthermore, because it is important to quantitatively evaluate available land use and land cover products, we compared our data with the MODIS Land Cover Type product (MCD12Q1). During the studied period, a large portion of dry fields was converted to perennial crops. Compared to our data, the forested area was underrepresented and the agricultural area overrepresented in MCD12Q1. In addition, linear landscape elements such as waterbodies were missing in the MODIS product due to its coarse spatial resolution. The data presented here can be useful for earth science and ecosystem services research.

Pan-European land cover/use map for 2009

GlobCorine demonstrated an automatic service that can generate in a consistent way land cover / land use maps and land change indicators, based on a CLC-compatible legend. CLC is derived from a visual identification and classification of landscape objects using high resolution images. This methodology provides high thematic accuracy but limits the update rate since it is time-consuming. Therefore, the project evaluated the use of MERIS FR time series, processed automatically to provide a more frequent update of CLC-compatible maps. GlobCorine built upon the experience and resources available through the GlobCover project, to tune the classification chain and adapt it to the EEA needs, covering the pan-European area (including the Mediterranean basin and the European Russia), although the system could be potentially extendable globally. The project delivered two CLC-compatible pan-European land cover maps in less than two years, demonstrating efficient and quick production. The first map is based on Envisat MERIS fine resolution (300m) mode data acquired between end 2004 and mid 2006, while the second used full-year 2009 data. GlobCorine is an initiative of ESA with the partnership of EEA and is implemented by Universite' catholique de Louvain - UCL.