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.

Calidad del agua en el área regadía del Río Mendoza (Argentina)

El río Mendoza conforma el oasis norte que es el más importante de la provincia. El crecimiento urbano ha avanzado sobre áreas originalmente agrícolas, rodeando la red de canales y desagües, que también recibe los desagües pluviales urbanos, producto de tormentas convectivas. La actividad antropogénica utiliza el recurso para bebida, saneamiento, riego, recreación, etc., y vuelca sus excedentes a la red, contaminándola. Para conocer la calidad del agua de esta cuenca se seleccionaron, estratégicamente, 15 sitios de muestreo: 3 a lo largo del río y a partir del dique derivador Cipolletti (R_I a R_III), 5 en la red de canales (C_I a C_V) y 7 ubicados en los colectores de drenaje (D_I a D_VII). Se realizaron los siguientes análisis físico-químicos y microbiológicos; en el río y en la red de canales: conductividad eléctrica, temperatura, pH, aniones y cationes (cálculo de RAS), oxígeno disuelto (OD), sólidos sedimentables, demanda química de oxígeno (DQO), bacterias aerobias mesófilas (BAM), coliformes totales y fecales y metales pesados. En la red de drenaje sólo se realizaron los cuatro primeros. Los resultados de los análisis, se incorporaron a una base de datos y se sometieron a un análisis estadístico descriptivo e inferencial. Este último consistió en la aplicación de diversas pruebas en busca de posibles diferencias entre los sitios de muestreo, para cada variable respuesta, a un α = 0.05. Se realizó el análisis de la varianza de efectos fijos y de efectos aleatorios y se probaron los supuestos de homocedasticidad y de normalidad de los errores. En el caso de violación de los supuestos, se utilizó la prueba de Kruskal- Wallis. Se compararon los siguientes sitios de muestreo entre sí: ríos, R_I-canales y drenajes. Se concluyó que hay un aumento significativo de la salinidad y la sodicidad en R_II, que los cambios de calidad ocurridos entre R_II y R_III podrían deberse al aporte de otras aguas. Con respecto a la comparación de los parámetros entre la cabeza del sistema (R_I) y la red de canales se puede decir que los aportes realizados por los escurrimientos urbanos ubicados hacia el oeste del canal Cacique Guaymallén, sumados a los vuelcos de Campo Espejo (detectados en C_II), incrementan significativamente la salinidad (+55 %) y sodicidad del agua (+95 %) respecto del punto R_I, aunque el valor de sodicidad sigue siendo bajo. También se han encontrado incrementos de salinidad (+80 %), de DQO (+1159 %) y BAM (+2873 %) con lógica disminución de OD (-58 %) en el punto C_V (canal Auxiliar Tulumaya) respecto del punto R_I, ocasionados por aportes urbanos (Gran Mendoza) sumados a la carga contaminante del canal Pescara. Los metales pesados no presentan grandes diferencias entre sitios de muestreo.

Rosalia de Castro y la visión poética de Galicia

Este artículo estudia la visión de la poeta española Rosalía de Castro (1837-1885) de su Galicia natal a partir de algunos de sus poemas, en especial del libro En las orillas del Sar (1884), la colección lírica considerada como la más importante de la autora. Desde el análisis de algunas piezas significativas sobre el tema de Galicia, se cubre inicialmente su visión personal de contestación a favor de su tierra natal. La segunda parte del estudio se centra en la temática social en Rosalía de Castro en relación con la defensa del emigrante gallego.

(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.

Accumulation, velocities, fluxes and mass balance of glacial tributaries to the Lambert Glacier-Amery Ice Shelf system, Antarctica

By incorporating recently available remote sensing data, we investigated the mass balance for all individual tributary glacial basins of the Lambert Glacier-Amery Ice Shelf system, East Antarctica. On the basis of the ice flow information derived from SAR interferometry and ICESat laser altimetry, we have determined the spatial configuration of eight tributary drainage basins of the Lambert-Amery glacial system. By combining the coherence information from SAR interferometry and the texture information from SAR and MODIS images, we have interpreted and refined the grounding line position. We calculated ice volume flux of each tributary glacial basin based on the ice velocity field derived from Radarsat three-pass interferometry together with ice thickness data interpolated from Australian and Russian airborne radio echo sounding (RES) surveys and inferred from ICESat laser altimetry data. Our analysis reveals that three tributary basins have a significant net positive imbalance, while five other subbasins are slightly positive or close to zero balance. Overall, in contrast to previous studies, we find that the grounded ice in Lambert Glacier-Amery Ice Shelf system has a positive mass imbalance of 22.9 ± 4.4 Gt/a. The net basal melting for the entire Amery Ice Shelf is estimated to be 27.0 ± 7.0 Gt/a. The melting rate decreases rapidly from the grounding zone to the ice shelf front. Significant basal refreezing is detected in the downstream section of the ice shelf. The mass balance estimates for both the grounded ice sheet and the ice shelf mass differ substantially from other recent estimates.