Controls on the Basal Water Pressure in Subglacial Channels Near the Margin of the Greenland Ice Sheet

Assuming a channelized drainage system in steady state, we investigate the influence of enhanced surface melting on the water pressure in subglacial channels, compared to that of changes in conduit geometry, ice rheology and catchment variations. The analysis is carried out for a specific part of the western Greenland ice-sheet margin between 66 degrees N and 66 degrees 30' N using new high-resolution digital elevation models of the subglacial topography and the ice-sheet surface, based on an airborne ice-penetrating radar survey in 2003 and satellite repeat-track interferometric synthetic aperture radar analysis of European Remote-sensing Satellite 1 and 2 (ERS-1/-2) imagery, respectively. The water pressure is calculated up-glacier along a likely subglacial channel at distances of 1, 5 and 9 km from the outlet at the ice margin, using a modified version of Rothlisberger's equation. Our results show that for the margin of the western Greenland ice sheet, the water pressure in subglacial channels is not sensitive to realistic variations in catchment size and mean surface water input compared to small changes in conduit geometry and ice rheology.

(Table 1) Basin area, average late-summer speedup, number of sinks and runoff of West Greenland Ice Sheet catchments

We use interferometric synthetic aperture radar observations recorded in a land-terminating sector of western Greenland to characterise the ice sheet surface hydrology and to quantify spatial variations in the seasonality of ice sheet flow. Our data reveal a non-uniform pattern of late-summer ice speedup that, in places, extends over 100 km inland. We show that the degree of late-summer speedup is positively correlated with modelled runoff within the 10 glacier catchments of our survey, and that the pattern of late-summer speedup follows that of water routed at the ice sheet surface. In late-summer, ice within the largest catchment flows on average 48% faster than during winter, whereas changes in smaller catchments are less pronounced. Our observations show that the routing of seasonal runoff at the ice sheet surface plays an important role in shaping the magnitude and extent of seasonal ice sheet speedup.

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

Ice thickness measurements in the Lincoln Sea and adjacent Arctic Ocean

Results of helicopter-borne electromagnetic measurements of total (ice plus snow) sea-ice thickness performed in May 2004 and 2005 in the Lincoln Sea and adjacent Arctic Ocean up to 86° N are presented. Thickness distributions south of 84° N are dominated by multi-year ice with modal thicknesses of 3.9 m in 2004 and 4.2 m in 2005 (mean thicknesses 4.67 and 5.18 m, respectively). Modal and mean snow thickness on multi-year ice amounted to 0.18 and 0.30 m in 2004, and 0.28 and 0.35 m in 2005. There are also considerable amounts of 0.9-2.2 m thick first-year ice (modal thickness), mostly representing ice formed in the recurring, refrozen Lincoln Polynya. Results are in good agreement with ground-based electromagnetic thickness measurements and with ice types demarcated in satellite synthetic aperture radar imagery. Four drifting buoys deployed in 2004 between 86° N and 84.5° N show a similar pattern of a mean southward drift of the ice pack of 83 ± 18 km between May 2004 and April 2005, towards the coast of Ellesmere Island and Nares Strait. The resulting area decrease of 26% between the buoys and the coast is larger than the observed thickness increase south of 84° N. This points to the importance of shear in a narrow band along the coast, and of ice export through Nares Strait in removing ice from the study region.

Mass balance and acceleration of the Greenland ice sheet

Within the last decade, the Greenland ice sheet (GrIS) and its surroundings have experienced record high surface temperatures (Mote, 2007, doi:10.1029/2007GL031976; Box et al., 2010), ice sheet melt extent (Fettweis et al., 2011, doi:10.5194/tc-5-359-2011) and record-low summer sea-ice extent (Nghiem et al., 2007, doi:10.1029/2007GL031138). Using three independent data sets, we derive, for the first time, consistent ice-mass trends and temporal variations within seven major drainage basins from gravity fields from the Gravity Recovery and Climate Experiment (GRACE; Tapley et al., 2004, doi:10.1029/2004GL019920), surface-ice velocities from Inteferometric Synthetic Aperture Radar (InSAR; Rignot and Kanagaratnam, 2006, doi:10.1126/science.1121381) together with output of the regional atmospheric climate modelling (RACMO2/ GR; Ettema et al., 2009, doi:10.1029/2009GL038110), and surface-elevation changes from the Ice, cloud and land elevation satellite (ICESat; Sorensen et al., 2011, doi:10.5194/tc-5-173-2011). We show that changing ice discharge (D), surface melting and subsequent run-off (M/R) and precipitation (P) all contribute, in a complex and regionally variable interplay, to the increasingly negative mass balance of the GrIS observed within the last decade. Interannual variability in P along the northwest and west coasts of the GrIS largely explains the apparent regional mass loss increase during 2002-2010, and obscures increasing M/R and D since the 1990s. In winter 2002/2003 and 2008/2009, accumulation anomalies in the east and southeast temporarily outweighed the losses by M/R and D that prevailed during 2003-2008, and after summer 2010. Overall, for all basins of the GrIS, the decadal variability of anomalies in P, M/R and D between 1958 and 2010 (w.r.t. 1961-1990) was significantly exceeded by the regional trends observed during the GRACE period (2002-2011).

Circumpolar Scaling Layer for downscaling of Scatterometer derived soil moisture with links to geotiff images and NetCDF files (2005-01 to 2011-12)

This dataset provides scaling information applicable to satellite derived coarse resolution surface soil moisture datasets following the approach by Wagner et al. (2008). It is based on ENVISAT ASAR data and can be utilized to apply the Metop ASCAT dataset (25 km) for local studies as well as to assess the representativeness of in-situ measurement sites and thus their potential for upscaling. The approach based on temporal stability (Wagner et al. 2008) consists of the assessment of the validity of the coarse resolution datasets at medium resolution (1 km, product is the so called 'scaling layer').

Surface velocity and ice discharge of the ice cap on King George Island, Antarctica

Glaciers on King George Island, Antarctica, have shown retreat and surface lowering in recent decades, concurrent with increasing air temperatures. A large portion of the glacier perimeter is ocean-terminating, suggesting possible large mass losses due to calving and submarine melting. Here we estimate the ice discharge into the ocean for the King George Island ice cap. L-band synthetic aperture radar images covering the time-span January 2008 to January 2011 over King George Island are processed using an intensity-tracking algorithm to obtain surface velocity measurements. Pixel offsets from 40 pairs of radar images are analysed and inverted to estimate a weighted average surface velocity field. Ice thicknesses are derived from simple principles of ice flow mechanics using the computed surface velocity fields and in situ thickness data. The maximum ice surface speeds reach mayor que 225 m/yr, and the total ice discharge for the analysed flux gates of King George Island is estimated to be 0.720+/-0.428 Gt/yr, corresponding to a specific mass loss of 0.64+/-0.38 m w.e./yr over the area of the entire ice cap (1127 km2).

Performance of TerraSAR-X for urban subsidence monitoring: Murcia case study

This paper presents an analysis of the performance of TerraSAR-X for subsidence monitoring in urban areas. The city of Murcia has been selected as a test-site due to its high deformation rate and the set of extensometers deployed along the city that provide validation data. The obtained results have been compared with those obtained from ERS/ENVISAT data belonging to the same period and validated with the in-situ measurements.

Influence of Incidence Angle on the Coherent Copolar Polarimetric Response of Rice at X-Band

The coherent nature of the acquisition by TerraSAR-X of both copolar channels (HH and VV) enables the generation of many different polarimetric observables with physical interpretation, as have recently been used for monitoring rice fields. In this letter, the influence of incidence angle upon these polarimetric observables is analyzed by comparing three stacks of images that were acquired simultaneously at different incidence angles (22°, 30°, and 40°) during a whole cultivation campaign. We show that the response of observables related to dominance (entropy, ratios of components) and type of scattering mechanisms (alpha angles) is not greatly influenced by incidence angle at some stages: early and advanced vegetative phases, and maturation. Moreover, the acquisition geometry drives the sensitivity to the presence of the initial stems and tillers, being detected earlier at shallower angles. This analysis is a necessary step before studying potential methodologies for combining different orbits and beams for reducing the time between acquisitions for monitoring purposes.

Remote-Sensing Monitoring of Tide Propagation Through Coastal Wetlands

Tide propagation through coastal wetlands is a complex phenomenon affected by vegetation, channels, and tidal conditions. Generally, tidal flow is studied using stage (water level) observations, which provide good temporal resolution, but they are acquired in limited locations. Here, a remote-sensing technique, wetland InSAR (interferometric synthetic aperture radar), is used to detect tidal flow in vegetated coastal environments over broad spatial scales. The technique is applied to data sets acquired by three radar satellites over the western Everglades in south Florida. Interferometric analysis of the data shows that the greatest water-level changes occur along tidal channels, reflecting a high velocity gradient between fast horizontal flow in the channel and the slow flow propagation through the vegetation. The high-resolution observations indicate that the tidal flushing zone extends 2–3 km on both sides of tidal channels and can extend 3–4 km inland from the end of the channel. The InSAR observations can also serve as quantitative constraints for detailed coastal wetland flow models.

Mass balances of different sectors in Antarctica in 1996, 2000 and 2006

Large uncertainties remain in the current and future contribution to sea level rise from Antarctica. Climate warming may increase snowfall in the continent's interior, but enhance glacier discharge at the coast where warmer air and ocean temperatures erode the buttressing ice shelves. Here, we use satellite interferometric synthetic-aperture radar observations from 1992 to 2006 covering 85% of Antarctica's coastline to estimate the total mass flux into the ocean. We compare the mass fluxes from large drainage basin units with interior snow accumulation calculated from a regional atmospheric climate model for 1980 to 2004. In East Antarctica, small glacier losses in Wilkes Land and glacier gains at the mouths of the Filchner and Ross ice shelves combine to a near-zero loss of 4 ± 61 Gt/yr. In West Antarctica, widespread losses along the Bellingshausen and Amundsen seas increased the ice sheet loss by 59% in 10 years to reach 132 ± 60 Gt/yr in 2006. In the Peninsula, losses increased by 140% to reach 60 ± 46 Gt/yr in 2006. Losses are concentrated along narrow channels occupied by outlet glaciers and are caused by ongoing and past glacier acceleration. Changes in glacier flow therefore have a significant, if not dominant impact on ice sheet mass balance.

Wetlands in central Siberia from ASAR WS 2003, link to shapefiles

An area in central Siberia (partial coverage of Turukhansky und Yeniseysky districts) was investigated using satellite data. It covers freshwater ecosystems of non-forested peatlands in boreal forests. The satellite data represent the growing seasons of 2003/2004. Microwave data were acquired by the Advanced Synthetic Aperture Radar (ASAR) instrument onboard ENVISAT. The multi-temporal capabilities and resolution (150mx150m in WS mode) of the ASAR wide swath mode enabled the detection of dynamic features >2ha over this vast area. Scatterometer (QuikScat) data could be employed to distinguish hydro-periods. Wetland types have been identified on the basis of seasonal changes in backscatter. Results for peatlands have been compared with Russian forest inventory data which contain information on wetland distribution.

Tile-Based Two-Dimensional Phase Unwrapping for Digital Holography Using a Modular Framework

A variety of physical and biomedical imaging techniques, such as digital holography, interferometric synthetic aperture radar (InSAR), or magnetic resonance imaging (MRI) enable measurement of the phase of a physical quantity additionally to its amplitude. However, the phase can commonly only be measured modulo 2π, as a so called wrapped phase map. Phase unwrapping is the process of obtaining the underlying physical phase map from the wrapped phase. Tile-based phase unwrapping algorithms operate by first tessellating the phase map, then unwrapping individual tiles, and finally merging them to a continuous phase map. They can be implemented computationally efficiently and are robust to noise. However, they are prone to failure in the presence of phase residues or erroneous unwraps of single tiles. We tried to overcome these shortcomings by creating novel tile unwrapping and merging algorithms as well as creating a framework that allows to combine them in modular fashion. To increase the robustness of the tile unwrapping step, we implemented a model-based algorithm that makes efficient use of linear algebra to unwrap individual tiles. Furthermore, we adapted an established pixel-based unwrapping algorithm to create a quality guided tile merger. These original algorithms as well as previously existing ones were implemented in a modular phase unwrapping C++ framework. By examining different combinations of unwrapping and merging algorithms we compared our method to existing approaches. We could show that the appropriate choice of unwrapping and merging algorithms can significantly improve the unwrapped result in the presence of phase residues and noise. Beyond that, our modular framework allows for efficient design and test of new tile-based phase unwrapping algorithms. The software developed in this study is freely available.

Swell dissipation from 10 years of Envisat ASAR in wave mode

Swells are found in all oceans and strongly influence the wave climate and air-sea processes. The poorly known swell dissipation is the largest source of error in wave forecasts and hindcasts. We use synthetic aperture radar data to identify swell sources and trajectories, allowing a statistically significant estimation of swell dissipation. We mined the entire Envisat mission 2003–2012 to find suitable storms with swells (13 < T < 18 s) that are observed several times along their propagation. This database of swell events provides a comprehensive view of swell extending previous efforts. The analysis reveals that swell dissipation weakly correlates with the wave steepness, wind speed, orbital wave velocity, and the relative direction of wind and waves. Although several negative dissipation rates are found, there are uncertainties in the synthetic aperture radar-derived swell heights and dissipation rates. An acceptable range of the swell dissipation rate is −0.1 to 6 × 10−7 m−1 with a median of 1 × 10−7 m−1.

Extrapolating Satellite Winds to Turbine Operating Heights

Ocean wind retrievals from satellite sensors are typically performed for the standard level of 10 m. This restricts their full exploitation for wind energy planning, which requires wind information at much higher levels where wind turbines operate. A new method is presented for the vertical extrapolation of satellite-based wind maps. Winds near the sea surface are obtained from satellite data and used together with an adaptation of the Monin–Obukhov similarity theory to estimate the wind speed at higher levels. The thermal stratification of the atmosphere is taken into account through a long-term stability correction that is based on numerical weather prediction (NWP) model outputs. The effect of the long-term stability correction on the wind profile is significant. The method is applied to Envisat Advanced Synthetic Aperture Radar scenes acquired over the south Baltic Sea. This leads to maps of the long-term stability correction and wind speed at a height of 100 m with a spatial resolution of 0.02°. Calculations of the corresponding wind power density and Weibull parameters are shown. Comparisons with mast observations reveal that NWP model outputs can correct successfully for long-term stability effects and also, to some extent, for the limited number of satellite samples. The satellite-based and NWP-simulated wind profiles are almost equally accurate with respect to those from the mast. However, the satellite-based maps have a higher spatial resolution, which is particularly important in nearshore areas where most offshore wind farms are built.