Surface velocity fields, digital elevation models and ice front positions derived from multi-mission SAR remote sensing data at Sjögren Inlet glaciers, Antarctic Peninsula

Substantial retreat or disintegration of numerous ice shelves have been observed on the Antarctic Peninsula. The ice shelf in the Prince Gustav Channel retreated gradually since the late 1980's and broke-up in 1995. Tributary glaciers reacted with speed-up, surface lowering and increased ice discharge, consequently contributing to sea level rise. We present a detailed long-term study (1993-2014) on the dynamic response of Sjögren Inlet glaciers to the disintegration of Prince Gustav Ice Shelf. We analyzed various remote sensing datasets to observe the reactions of the glaciers to the loss of the buttressing ice shelf. A strong increase in ice surface velocities was observed with maximum flow speeds reaching 2.82±0.48 m/d in 2007 and 1.50±0.32 m/d in 2004 at Sjögren and Boydell glaciers respectively. Subsequently, the flow velocities decelerated, however in late 2014, we still measured about two times the values of our first measurements in 1996. The tributary glaciers retreated 61.7±3.1 km² behind the former grounding line of the ice shelf. In regions below 1000 m a.s.l., a mean surface lowering of -68±10 m (-3.1 m/a) was observed in the period 1993-2014. The lowering rate decreased to -2.2 m/a in recent years. Based on the surface lowering rates, geodetic mass balances of the glaciers were derived for different time steps. High mass loss rate of -1.21±0.36 Gt/a was found in the earliest period (1993-2001). Due to the dynamic adjustments of the glaciers to the new boundary conditions the ice mass loss reduced to -0.59±0.11 Gt/a in the period 2012-2014, resulting in an average mass loss rate of -0.89±0.16 Gt/a (1993-2014). Including the retreat of the ice front and grounding line, a total mass change of -38.5±7.7 Gt and a contribution to sea level rise of 0.061±0.013 mm were computed. Analysis of the ice flux revealed that available bedrock elevation estimates at Sjögren Inlet are too shallow and are the major uncertainty in ice flux computations. This temporally dense time series analysis of Sjögren Inlet glaciers shows that the adjustments of tributary glaciers to ice shelf disintegration are still going on and provides detailed information of the changes in glacier dynamics.

Snow and ice surfaces on King George Island (Antarctic Peninsula) with high-resolution TerraSAR-X time series, with links to shapfiles and GeoTIFFs

Changes of glaciers and snow cover in polar regions affect a wide range of physical and ecosystem processes on land and in the adjacent marine environment. In this study, we investigate the potential of 11-day repeat high-resolution satellite image time series from the TerraSAR-X mission to derive glaciological and hydrological parameters on King George Island, Antarctica during the period Oct/25/2010 to Apr/19/2011. The spatial pattern and temporal evolution of snow cover extent on ice-free areas can be monitored using multi-temporal coherence images. SAR coherence is used to map glacier extent of land terminating glaciers with an average accuracy of 25 m. Multi-temporal SAR color composites identify the position of the late summer snow line at about 220 m above sea level. Glacier surface velocities are obtained from intensity feature-tracking. Surface velocities near the calving front of Fourcade Glacier were up to 1.8 ± 0.01 m/d. Using an intercept theorem based on fundamental geometric principles together with differential GPS field measurements, the ice discharge of Fourcade Glacier was estimated to 20700 ± 5500 m**3/d (corresponding to ~19 ± 5 kt/d). The rapidly changing surface conditions on King George Island and the lack of high-resolution digital elevation models for the region remain restrictions for the applicability of SAR data and the precision of derived products.

Chemical analysis and argon geochronology from dive SAR_1 on Cornacya seamount of the Sardinia Channel

Sarcya 1 dive explored a previously unknown 12 My old submerged volcano, labelled Cornacya. A well developed fracturation is characterised by the following directions: N 170 to N-S, N 20 to N 40, N 90 to N 120, N 50 to N 70, which corresponds to the fracturation pattern of the Sardinian margin. The sampled lavas exhibit features of shoshonitic suites of intermediate composition and include amphibole-and mica-bearing lamprophyric xenoliths which are geochemically similar to Ti-poor lamproites. Mica compositions reflect chemical exchanges between the lamprophyre and its shoshonitic host rock suggesting their simultaneous emplacement. Nd compositions of the Cornacya K-rich suite indicate that continental crust was largely involved in the genesis of these rocks. The spatial association of the lamprophyre with the shoshonitic rocks is geochemically similar to K-rich and TiO2-poor igneous suites, emplaced in post-collisional settings. Among shoshonitic rocks, sample SAR 1-01 has been dated at 12.6±0.3 My using the 40Ar/39Ar method with a laser microprobe on single grains. The age of the Cornacya shoshonitic suite is similar to that of the Sisco lamprophyre from Corsica, which similarly is located on the western margin of the Tyrrhenian Sea. Thus, the Cornacya shoshonitic rocks and their lamprophyric xenolith and the Sisco lamprophyre could represent post-collisional suites emplaced during the lithospheric extension of the Corsica-Sardinia block, just after its rotation and before the Tyrrhenian sea opening. Drilling on the Sardinia margin (ODP Leg 107) shows that the upper levels of the present day margin (Hole 654) suffered tectonic subsidence before the lower part (Hole 652). The structure of this lower part is interpreted as the result of an eastward migration of the extension during Late Miocene and Early Pliocene times. Data of Cornacya volcano are in good agreement with this model and provide good chronological constraints for the beginning of the phenomenon.

Retreat of top cliff of Kurungnakh Island, Lena Delta, Siberia, Russia, 2010-2014, with links to shapefiles

Thawing-induced cliff top retreat in permafrost landscapes is mainly due to thermo-erosion. Ground-ice-rich permafrost landscapes are specifically vulnerable to thermo-erosion and may show high degradation rates. Within the HGF Alliance Remote Sensing and the FP7 PAGE21 permafrost programs we investigated how SAR and optical remote sensing can contribute to the monitoring of erosion rates of ice-rich cliffs in Arctic Siberia (Lena Delta, Russia). We produced two different vector products: i) Intra-annual cliff top retreat based on TerraSAR-X (TSX) satellite data (2012-2014): High-temporal resolution time series of TSX satellite data allow the inter-annual and intra-annual monitoring of the upper cliff-line retreat also under bad weather conditions and continuous cloud coverage. This published SAR product contains the retreating upper cliff lines of a 1.5 km long part of eroding ice-rich coast of Kurungnakh Island in the central Lena Delta. The upper cliff line was mapped using a thresholding approach for images acquired in the years 2012, 2013 and 2014 for the months June (2013, 2014), July (2013, 2014), August (2012, 2013, 2014) and September (2013, 2014). The cliff top retreat vector product is called 'upper_cliff_TerraSAR-X'. While the 2014 cliff lines show a clear retreat of 2 to 3 m/month, the cliff top lines for 2012 and 2013 are not chronologically ordered. However, lines from the end of the season of a year are always close to the lines from the beginning of the next summer season, indicating low cliff retreat in winter. ii) 4-year cliff top retreat based on optical satellite data (2010-2014): Long-term cliff top retreat could be assessed with two high-spatial resolution optical satellite images (GeoEye-1, 2010-08-05 and Worldview-1, 2014-08-19). The cliff top retreat vector product is called 'upper_cliff_optical'. Results: The long-term cliff top retreat derived from optical satellite data are 35 m cliff retreat within 4 years. The higher-temporal resolution SAR data equivalently show long-term rates of 18 m within 2 years and nearly now degradation activities in winter but maximum erosion rates in summer months.The Intra-seasonal cliff top retreat lines from 2014 show a rate of 2 to 3 m per month.

Chronostratigraphic framework at sites GeoB14403 and GeoB14414 from the Gela Basin, Sicily Channel

A multi-proxy chronological framework along with sequence-stratigraphic interpretations unveils composite Milankovitch cyclicity in the sedimentary records of the Last GlacialeInterglacial cycle at NE Gela Basin on the Sicilian continental margin. Chronostratigraphic data (including foraminifera-based eco-biostratigraphy and d18O records, tephrochronological markers and 14C AMS radiometric datings) was derived from the shallow-shelf drill sites GeoB14403 (54.6 m recovery) and GeoB14414 (27.5 m), collected with both gravity and drilled MeBo cores in 193 m and 146 m water depth, respectively. The recovered intervals record Marine Isotope Stages and Substages (MIS) from MIS 5 to MIS 1, thus comprising major stratigraphic parts of the progradational deposits that form the last 100-ka depositional sequence. Calibration of shelf sedimentary units with borehole stratigraphies indicates the impact of higher-frequency (20-ka) sea level cycles punctuating this 100-ka cycle. This becomes most evident in the alternation of thick interstadial highstand (HST) wedges and thinner glacial forced-regression (FSST) units mirroring seaward shifts in coastal progradation. Albeit their relatively short-lived depositional phase, these subordinate HST units form the bulk of the 100-ka depositional sequence. Two mechanisms are proposed that likely account for enhanced sediment accumulation ratios (SAR) of up to 200 cm/ka during these intervals: (1) intensified activity of deep and intermediate Levantine Intermediate Water (LIW) associated to the drowning of Mediterranean shelves, and (2) amplified sediment flux along the flooded shelf in response to hyperpycnal plumes that generate through extreme precipitation events during overall arid conditions. Equally, the latter mechanism is thought to be at the origin of undulated features resolved in the acoustic records of MIS 5 Interstadials, which bear a striking resemblance to modern equivalents forming on late-Holocene prodeltas of other Mediterranean shallow-shelf settings.

Physical properties of different sea ice types at ArcticNet, CFL and CASES ice stations, eastern Beaufort Sea

A study of the polarimetric backscattering response of newly formed sea ice types under a large assortment of surface coverage was conducted using a ship-based C-band polarimetric radar system. Polarimetric backscattering results and physical data for 40 stations during the fall freeze-up of 2003, 2006, and 2007 are presented. Analysis of the copolarized correlation coefficient showed its sensitivity to both sea ice thickness and surface coverage and resulted in a statistically significant separation of ice thickness into two regimes: ice less than 6 cm thick and ice greater than 8 cm thick. A case study quantified the backscatter of a layer of snow infiltrated frost flowers on new sea ice, showing that the presence of the old frost flowers can enhance the backscatter by more than 6 dB. Finally, a statistical analysis of a series of temporal-spatial measurements over a visually homogeneous frost-flower-covered ice floe identified temperature as a significant, but not exclusive, factor in the backscattering measurements.

(Table 1) Sea ice and snow characteristics and heat fluxes observed during R/V Aurora Australis cruise to East Antarctica in September/October 2007

The properties of snow on East Antarctic sea ice off Wilkes Land were examined during the Sea Ice Physics and Ecosystem Experiment (SIPEX) in late winter of 2007, focusing on the interaction with sea ice. This observation includes 11 transect lines for the measurement of ice thickness, freeboard, and snow depth, 50 snow pits on 13 ice floes, and diurnal variation of surface heat flux on three ice floes. The detailed profiling of topography along the transects and the d18O, salinity, and density datasets of snow made it possible to examine the snow-sea-ice interaction quantitatively for the first time in this area. In general, the snow displayed significant heterogeneity in types, thickness (mean: 0.14 +- 0.13 m), and density (325 +- 38 kg/m**3), as reported in other East Antarctic regions. High salinity was confined to the lowest 0.1 m. Salinity and d18O data within this layer revealed that saline water originated from the surface brine of sea ice in 20% of the total sites and from seawater in 80%. From the vertical profiles of snow density, bulk thermal conductivity of snow was estimated as 0.15 W/K/m on average, only half of the value used for numerical sea-ice models. Although the upward heat flux within snow estimated with this value was significantly lower than that within ice, it turned out that a higher value of thermal conductivity (0.3 to 0.4 W/K/m) is preferable for estimating ice growth amount in current numerical models. Diurnal measurements showed that upward conductive heat flux within the snow and net long-wave radiation at the surface seem to play important roles in the formation of snow ice from slush. The detailed surface topography allowed us to compare the air-ice drag coefficients of ice and snow surfaces under neutral conditions, and to examine the possibility of the retrieval of ice thickness distribution from satellite remote sensing. It was found that overall snow cover works to enhance the surface roughness of sea ice rather than moderate it, and increases the drag coefficient by about 10%. As for thickness retrieval, mean ice thickness had a higher correlation with ice surface roughness than mean freeboard or surface elevation, which indicates the potential usefulness of satellite L-band SAR in estimating the ice thickness distribution in the seasonal sea-ice zone.

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.

Pollen and macrofossils profiles and age determination from three cores in Northeastern Kansas, U.S.A.

Two marshes near Muscotah and Arrington, Atchison County, northeastern Kansas, yielded a pollen sequence covering the last 25,000 yrs of vegetation development. The earliest pollen spectra are comparable with surface pollen spectra from southern Saskatchewan and southeastern Manitoba and might indicate a rather open vegetation but with some pine, spruce, and birch as the most important tree species, with local stands of alder and willow. This type of vegetation changed about 23,000 yrs ago to a spruce forest, which prevailed in the region until at least 15,000 yrs ago. Because of a hiatus, the vegetation changes resulting in the spread of a mixed deciduous forest and prairie, which was present in the region from 11,000 to 9,000 yrs ago, remain unknown. Prairie vegetation, with perhaps a few trees along the valleys, covered the region until about 5,000 yrs ago, when a re-expansion of deciduous trees began in the lowlands.

Pollen and spore counts, influx rates and concentrations at DSDP Hole 64-480

Preliminary analyses of 56 samples from the upper 49 meters of Hole 480 (Cores 1-11) show marked changes in pollen frequencies and concentrations. The largely varved cores (1, 2, 3, 10, and 11) are characterized by low concentrations and pollen types such as Gramineae, Low-spine Compositae and Cheno/Ams. The largely homogeneous section (Cores 3 through 10) contains higher pollen concentrations and is dominated by TCT (probably Juniperus) and Artemisia. Picea pollen is also present in this section. The record as a whole is thought to represent most of the last glacial cycle.

Snow-pit measurements on Atka Bay landfast sea ice during ANT-Land 2012/2013 season

Up to now, snow cover on Antarctic sea ice and its impact on radar backscatter, particularly after the onset of freeze/thaw processes, are not well understood. Here we present a combined analysis of in situ observations of snow properties from the landfast sea ice in Atka Bay, Antarctica, and high-resolution TerraSAR-X backscatter data, for the transition from austral spring (November 2012) to summer (January 2013). The physical changes in the seasonal snow cover during that time are reflected in the evolution of TerraSAR-X backscatter. We are able to explain 76-93% of the spatio-temporal variability of the TerraSAR-X backscatter signal with up to four snowpack parameters with a root-mean-squared error of 0.87-1.62 dB, using a simple multiple linear model. Over the complete study, and especially after the onset of early-melt processes and freeze/thaw cycles, the majority of variability in the backscatter is influenced by changes in snow/ice interface temperature, snow depth and top-layer grain size. This suggests it may be possible to retrieve snow physical properties over Antarctic sea ice from X-band SAR backscatter.

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.

Icebergs in the Amundsen Sea (Antarctica) determined from ENVISAT ASAR WSM data using object-based image analysis approach, link to ESRI GDB files

ENVISAT ASAR WSM images with pixel size 150 × 150 m, acquired in different meteorological, oceanographic and sea ice conditions were used to determined icebergs in the Amundsen Sea (Antarctica). An object-based method for automatic iceberg detection from SAR data has been developed and applied. The object identification is based on spectral and spatial parameters on 5 scale levels, and was verified with manual classification in four polygon areas, chosen to represent varying environmental conditions. The algorithm works comparatively well in freezing temperatures and strong wind conditions, prevailing in the Amundsen Sea during the year. The detection rate was 96% which corresponds to 94% of the area (counting icebergs larger than 0.03 km**2), for all seasons. The presented algorithm tends to generate errors in the form of false alarms, mainly caused by the presence of ice floes, rather than misses. This affects the reliability since false alarms were manually corrected post analysis.