AbstractAirborne radio-echo sounding (RES) surveys are widely used to measure ice-sheet bed topography. Measuring bed topography as accurately and widely as possible is of critical importance to modelling ice dynamics and hence to constraining better future ice response to climate change. Measurement accuracy of RES surveys is influenced both by the geometry of bed topography and the survey design. Here we develop a novel approach for simulating RES surveys over glaciated terrain, to quantify the sensitivity of derived bed elevation to topographic geometry. Furthermore, we investigate how measurement errors influence the quantification of glacial valley geometry. We find a negative bias across RES measurements, where off-nadir return measurement error is typically −1.8 ± 11.6 m. Topographic highlands are under-measured an order of magnitude more than lowlands. Consequently, valley depth and cross-sectional area are largely under-estimated. While overall estimates of ice thickness are likely too high, we find large glacier valley cross-sectional area to be under-estimated by −2.8 ± 18.1%. Therefore, estimates of ice flux through large outlet glaciers are likely too low when this effect is not taken into account. Additionally, bed mismeasurements potentially impact our appreciation of outlet-glacier stability.

AbstractGreenland's bed topography is a primary control on ice flow, grounding line migration, calving dynamics, and subglacial drainage. Moreover, fjord bathymetry regulates the penetration of warm Atlantic water (AW) that rapidly melts and undercuts Greenland's marine‐terminating glaciers. Here we present a new compilation of Greenland bed topography that assimilates seafloor bathymetry and ice thickness data through a mass conservation approach. A new 150 m horizontal resolution bed topography/bathymetric map of Greenland is constructed with seamless transitions at the ice/ocean interface, yielding major improvements over previous data sets, particularly in the marine‐terminating sectors of northwest and southeast Greenland. Our map reveals that the total sea level potential of the Greenland ice sheet is 7.42 ± 0.05 m, which is 7 cm greater than previous estimates. Furthermore, it explains recent calving front response of numerous outlet glaciers and reveals new pathways by which AW can access glaciers with marine‐based basins, thereby highlighting sectors of Greenland that are most vulnerable to future oceanic forcing.

We present a new bed elevation dataset for Greenland derived from a combination of multiple airborne ice thickness surveys undertaken between the 1970s and 2012. Around 420 000 line kilometres of airborne data were used, with roughly 70% of this having been collected since the year 2000, when the last comprehensive compilation was undertaken. The airborne data were combined with satellite-derived elevations for non-glaciated terrain to produce a consistent bed digital elevation model (DEM) over the entire island including across the glaciated-ice free boundary. The DEM was extended to the continental margin with the aid of bathymetric data, primarily from a compilation for the Arctic. Ice thickness was determined where an ice shelf exists from a combination of surface elevation and radar soundings. The across-track spacing between flight lines warranted interpolation at 1 km postings for significant sectors of the ice sheet. Grids of ice surface elevation, error estimates for the DEM, ice thickness and data sampling density were also produced alongside a mask of land/ocean/grounded ice/floating ice. Errors in bed elevation range from a minimum of ±10 m to about ±300 m, as a function of distance from an observation and local topographic variability. A comparison with the compilation published in 2001 highlights the improvement in resolution afforded by the new datasets, particularly along the ice sheet margin, where ice velocity is highest and changes in ice dynamics most marked. We estimate that the volume of ice included in our land-ice mask would raise mean sea level by 7.36 m, excluding any solid earth effects that would take place during ice sheet decay. © 2013 Author(s).

AbstractIce discharge from Flade Isblink ice cap (NE Greenland) maintains an ice shelf at the northwestern fringe of the ice cap. The two outlet basins feeding this ice shelf surged during the late 1990s. Ice shelves are rare in Greenland and surges of ice shelf terminating glaciers even rarer. Understanding and explaining the evolution of ice mass changes in the two basins is hampered by a lack of knowledge about processes at their grounding zones. We determined, for the first time, the grounding line locations of these basins and analyzed their variability with time. We further quantified ice discharge and its variability during the period 1988–2020. We found that the grounding lines advanced slightly between 1993 and 1999 during the glacier surges, but showed overall retreats of 2.2 ± 1.3 km in basin 2 and 2.7 ± 0.9 km in basin 3 until 2019 over retrograde sloping beds. The retreats were promoted by increasing buoyancy forces due to increasing water depth, but opposing buttressing forces of the ice shelf induced a differing response of the grounding line in the two basins. Based on the observed patterns of flow and retreat, we characterized the surges as “Svalbard‐type”, modified by buttressing effects of the ice shelf. We calculated total ice discharges over the study period of 1.85 ± 1.59 Gt in basin 2 and 1.38 ± 1.22 Gt in basin 3. We observed reductions in ice discharge of at least 90% after the surges, that persisted for the remainder of the period studied.

Due to Arctic amplification, impacts of surface warming are significantly observed at higher latitudes. Positive feedbacks between periglacial environments and climate have increased thermokarst extent and associated landforms. Subsidence rates depend on active layer sensitivity which controls the geothermal heat balance between the Earth’s surface and interior. Carbon release is important, both gradually as the active layer seasonally thaws and deepens, and more notably as abrupt thaw mobilises deep soil organic carbon (SOC). Hillslope processes responsible for abrupt thaw are less well studied using remote sensing due to detection difficulties. Using high-resolution elevation data, this study increases understandings of relationships between topographic settings and climatic forcings by quantifying rates and magnitudes of geomorphic change for features identified within previous studies. Digital elevation model (DEM) differencing techniques are applied to ice wedges at three sites on Garry Island (Canadian High Arctic), where research has historically been limited to in-situ measurement. The aim is to assess the capabilities of ArcticDEM. Results detected rates of vertical change comparable to field studies. Patterns of positive elevation change were explored through hydrological network analysis and snowblow modelling. The differencing method is applied to mass movements (active layer detachment slides (ALDS) and retrogressive thaw slumps (RTS)) at four sites in the Alaskan Brooks Range. Abrupt thaw events on north facing slopes increased in magnitude, with length exhibiting the most significant rate of change in both features. Semi-automated delineation techniques developed using ArcGIS Model Builder, aided the mapping of five previously undetected features and facilitated carbon release estimates. As the Arctic warms, permafrost will continue to thaw, and hillslopes will become more unstable, impacting spatial extents of periglacial landforms. This thesis uses repeat elevation data at high resolutions in a method which could be implemented across the High Arctic to provide contributions to studies of periglacial environments.

AbstractAirborne radio-echo sounding (RES) surveys are widely used to measure ice-sheet bed topography. Measuring bed topography as accurately and widely as possible is of critical importance to modelling ice dynamics and hence to constraining better future ice response to climate change. Measurement accuracy of RES surveys is influenced both by the geometry of bed topography and the survey design. Here we develop a novel approach for simulating RES surveys over glaciated terrain, to quantify the sensitivity of derived bed elevation to topographic geometry. Furthermore, we investigate how measurement errors influence the quantification of glacial valley geometry. We find a negative bias across RES measurements, where off-nadir return measurement error is typically −1.8 ± 11.6 m. Topographic highlands are under-measured an order of magnitude more than lowlands. Consequently, valley depth and cross-sectional area are largely under-estimated. While overall estimates of ice thickness are likely too high, we find large glacier valley cross-sectional area to be under-estimated by −2.8 ± 18.1%. Therefore, estimates of ice flux through large outlet glaciers are likely too low when this effect is not taken into account. Additionally, bed mismeasurements potentially impact our appreciation of outlet-glacier stability.

Contribution to global mean sea level rise by ice sheets, ice caps and glaciers is
accelerating. The total volume of water stored globally in terrestrial ice is
estimated by a multitude of methods but principally by the interpolation of icethickness data. For the ice sheets and large Arctic ice caps, ice thickness is
predominantly measured by airborne radio-echo sounding surveys which use
radio waves to detect the bed of the surveyed ice mass. While such surveys are
now extensive, large portions of ice masses are generally unsurveyed due to their
size. In order to quantify ice thickness and subsequently ice volume over the
entirety of an ice mass, interpolation of the input measurements is used.
Throughout this whole process, uncertainties arise. Initially, from the radio-echo
sounding (RES) survey and subsequently, in the interpolation. Compounding this
is the absence of ground-truthing for measurements and interpolations due to the
inaccessibility of ice mass beds. Hence, there is a requirement to find alternative
means of quantifying uncertainty in ice thickness measurements and
subsequently derived bed topography, and analyses made from these data to
reduce the uncertainty in sea level change projections.
This thesis develops and applies methods which aim to reduce uncertainty in ice
thickness and bed topography datasets. Using high-resolution elevation data, this
study exploits the likely similarity between currently ice-covered topography and
formerly glaciated topography in the Arctic to generate datasets which provide
alternative validation for ice mass bed topography. For the first time topographic
error in RES surveying is quantified and corrections are formulated for treating
future and historic ice thickness and bed topography data. Additionally, the
propagation of these uncertainties through interpolations of bed topography is
quantified and reduced, focussing on the Greenland Ice Sheet. Finally, the full
suite of methods is applied to ice caps in the Canadian Arctic to generate, for the
first time, ice cap wide topography for ice caps in the region that hold
approximately a third of the freshwater outside of the continental ice sheets. By
quantifying and reducing uncertainty in datasets of bed topography and ice
thickness this thesis assesses the perceived stability of the continental ice sheets
and large ice Arctic ice caps. From this, the implications of this for near and far
term global mean sea-level rise are investigated.

Dust storm events occur in arid and semi-arid areas around the world. These result from strong surface winds and blow dust and sand from loose, dry soil surfaces into the atmosphere. Such events can have damaging effects on human health, environment, infrastructure and transport. In the first section of this PhD dissertation, focus on the suitability of the existing of five different MODIS-based methods for detecting airborne dust over the Arabian Peninsula are examined. These are the: (a) Normalized Difference Dust Index (NDDI); (b) Brightness Temperature Difference (BTD) (Band 31–32); (c) BTD (Band 20–31); (d) Middle East Dust Index (MEDI) and (e) Reflective Solar Band (RSB). This work also develops dust detection thresholds for each index by comparing observed values for ‘dust-present’ versus ‘dust-free’ conditions, taking into account various land cover settings and analysing associated temporal trends. The results suggest the most suitable indices for identifying dust storms over different land cover types across the Arabian Peninsula are BTD31–32 and the RSB index. Methods such as NDDI and BTD20 – 31 have limitations in detecting dust over multiple land-cover types. In addition, MEDI was found to be an unsuccessful index for detecting dust storms over all types of land cover in the study area. Furthermore, this thesis explores the spatial and temporal variations of dust storms by using monthly meteorological data from 27 observation stations across Saudi Arabia during the period (2000–2016), considering the associations between dust storm frequency and temperature, precipitation and wind variables. In terms of the frequency of dust in Saudi Arabia, the results show significant spatial, seasonal and inter-annual. In the eastern part of the study area, for example, dust storm events have increased over time, especially in Al-Ahsa. There are evident relationships (p < 0.005) between dust storm occurrence and wind speed, wind direction and precipitation.
This thesis also describes the impact of dust on health, and specifically on respiratory admissions to King Fahad Medical City (KFMC) for the period (February 2015 – January 2016).This study uses dust data from the World Meteorological Or-ganization (WMO) for comparing and analysing the daily weather conditions and hospital admissions. The findings indicate that the total number of emergency respiratory admissions during dust events was higher than background levels by 36% per day on average. Numbers of admissions during ‘widespread dust’ events were 19.62% per day higher than during periods of ‘blowing dust’ activity. The average number of hospital admissions for lower respiratory tract infections (LRTI) was 11.62 per day during widespread dust events and 10.36 per day during blowing dust. The average number of hospital admissions for upper respiratory tract infections (URTI) was 10.25 per day during widespread dust events and 7.87 per day during blowing dust ones. I found clear seasonal variability with a peak in the number of emergency admissions during the months of February to April. Furthermore, qualitative evidence suggests that there is a significant impact on hospital operations due to the increase in patients and pressure on staffing and hospital consumables in this period.
Taken together, these findings suggest the (BTD 31–32) and (RSB) are the most suitable indices of the five different MODIS-based methods for detecting airborne dust over the Arabian Peninsula and over different land cover. There are important spatial and temporal pattern variations, as well as seasonal and inter-annual variability, in the occurrence of dust storms in Saudi Arabia. There is also a seasonal pat-tern to the number of hospital admissions during dust events. This is research in-tended to fill the knowledge gap in the dust detection filed. Here I address the knowledge gap by evaluating the identified dust methods over the whole Arabian Peninsula and by considering different land cover. To my knowledge, this is the first study analysed the temporal trends in indices values considering dust and dust-free conditions.
Previous work has only focused on 13 stations for analysing dust storms over Saudi Arabia. Therefore, this study has analysed the seasonal and inter-annual and spatial variation by using data from 27 observations in Saudi Arabia. This study addresses the relationship between dust storm frequency and the three meteorological factors (i.e. temperature, precipitation and wind variables) which have not yet been clarified in previous studies. In addition, this research fills the gap in the literature by investigating the correlation between different types of dust events such as (wide-spread dust and blowing dust) and their effects on the hospital admissions for upper and lower respiratory tract issues for pediatric. in Riyadh city.

Around 70-80% of the world’s coastline, and around 60% of the UK’s coastline, can be considered as ‘rocky’. Rocky coasts erode much slower than their softer sedimentary counterparts, but their rates of erosion and their evolutionary history are poorly known. In this dissertation I use a new combination of methods, cosmogenic nuclide exposure dating, structure-from-motion photogrammetry and sea-level modelling, to study a typical stretch of rocky coastline in north Torbay, Devon, southwest England. Torbay’s coast is characterised by the presence of shore platforms and raised beaches above modern sea level, situated on the north headland peninsula, named Hopes Nose. These elevated landforms must relate to a previous interglacial period, with warmer environments and higher sea-levels, and their preservation indicates very slow rates of coastal evolution within the area. I apply exposure dating using 36Cl to determine the degree of geomorphological inheritance from previous high sea-level stands, along north Torbay’s rocky cliffs and across the main body of the raised shore platform at Hopes Nose. I combine this analysis with the measurement of a new digital surface model, collected via drone imagery and structure-from-motion photogrammetry, across the headland to perform a morphometric analysis of the modern and elevated interglacial platform. Lastly, I determine a new estimate of relative sea-level change at the site, considering glacio-isostatic adjustment, using the SELEN sea-level model. Cosmogenic nuclide exposure dating reveals that the rocky coastline around north Torbay has been actively eroding throughout the late Holocene, through a series of stochastic mass movements and some incremental loss. Similarly, the exposure dating of the raised shore platform at Hopes Nose reveals it has been covered by distinctive sediments during the late Pleistocene and hence survived surface erosion. Morphometric analysis of the raised interglacial shore platform and the modern shore platform shows a similar evolutionary history, highlighting the changes in marine and aerial influence over the shore platforms formation. An analysis of the raised platform’s elevation, evaluated relative to modelled relative sea-level change, is most consistent with the platform being formed during the last interglacial period (Marine Isotope Stage 5e). As a result, this research also puts into question the overall height of the sea-level during MIS 5e, or the presence of a double peak within the record. Overall, this research demonstrates that the unique combination of methodologies can quantify coastal erosion and help decipher a rocky coastline’s history under both present and previous sea levels

Monthly meteorological data from 27 observation stations provided by the Presidency of Meteorology and Environment (PME) of Saudi Arabia were used to analyze the spatial and temporal distribution of atmospheric dust in Saudi Arabia between 2000 and 2016. These data were used to analyze the effects of environmental forcing on the occurrence of dust storms across Saudi Arabia by considering the relationships between dust storm frequency and temperature, precipitation, and wind variables. We reveal a clear seasonality in the reported incidence of dust storms, with the highest frequency of events during the spring. Our results show significant positive relationships (p < 0.005) between dust storm occurrence and wind speed, wind direction, and precipitation. However, we did not detect a significant relationship with temperature. Our results reveal important spatial patterns, as well as seasonal and inter-annual variations, in the occurrence of dust storms in Saudi Arabia. For instance, the eastern part of the study area experienced an increase in dust storm events over time, especially in the region near Al-Ahsa. Similarly, an increasing trend in dust storms was also observed in the west of the study area near Jeddah. However, the occurrence of dust storm events is decreasing over time in the north, in areas such as Hail and Qaisumah. Overall, the eastern part of Saudi Arabia experiences the highest number of dust storms per year (i.e. 10 to 60 events), followed by the northern region, with the south and the west having fewer dust storm events (i.e. five to 15 events per year). In addition, our results showed that the wind speeds during a dust storm are 15-20 m/s and above, while, on a non-dust day, the wind speeds are approximately 10-15 m/s or lower. Findings of this study provide insight into the relationship between environmental conditions and dust storm occurrence across Saudi Arabia, and a basis for future research into the drivers behind these observed spatio-temporal trends.

Sand and dust storm events (SDEs), which result from strong surface winds in arid and semi-arid areas, exhibiting loose dry soil surfaces are detrimental to human health, agricultural land, infrastructure, and transport. The accurate detection of near-surface dust is crucial for quantifying the spatial and temporal occurrence of SDEs globally. The Arabian Peninsula is an important source region for global dust due to the presence of extensive deserts. This paper evaluates the suitability of five different MODIS-based methods for detecting airborne dust over the Arabian Peninsula: (a) Normalized Difference Dust Index (NDDI); (b) Brightness Temperature Difference (BTD) (31-32); (c) BTD (20-31); (d) Middle East Dust Index (MEDI) and (e) Reflective Solar Band (RSB).We derive detection thresholds for each index by comparing observed values for 'dust-present' versus 'dust-free' conditions, taking into account various land cover settings and analyzing associated temporal trends. Our results suggest that the BTD (31-32) method and the RSB index are the most suitable indices for detecting dust storms over different land-cover types across the Arabian Peninsula. The NDDI and BTD (20-31) methods have limitations in identifying dust over multiple land-cover types. Furthermore, the MEDI has been found to be unsuitable for detecting dust in the study area across all land-cover types.

AbstractGreenland's bed topography is a primary control on ice flow, grounding line migration, calving dynamics, and subglacial drainage. Moreover, fjord bathymetry regulates the penetration of warm Atlantic water (AW) that rapidly melts and undercuts Greenland's marine‐terminating glaciers. Here we present a new compilation of Greenland bed topography that assimilates seafloor bathymetry and ice thickness data through a mass conservation approach. A new 150 m horizontal resolution bed topography/bathymetric map of Greenland is constructed with seamless transitions at the ice/ocean interface, yielding major improvements over previous data sets, particularly in the marine‐terminating sectors of northwest and southeast Greenland. Our map reveals that the total sea level potential of the Greenland ice sheet is 7.42 ± 0.05 m, which is 7 cm greater than previous estimates. Furthermore, it explains recent calving front response of numerous outlet glaciers and reveals new pathways by which AW can access glaciers with marine‐based basins, thereby highlighting sectors of Greenland that are most vulnerable to future oceanic forcing.

Glaciological and hydraulic factors that control the timing and mechanisms of glacier lake outburst floods (GLOFs) remain poorly understood. This study used measurements of lake level at 15min intervals and known lake bathymetry to calculate lake outflow during two GLOF events from the northern margin of Russell Glacier, west Greenland. We used measured ice surface elevation, interpolated subglacial topography and likely conduit geometry to inform a melt enlargement model of the outburst evolution. The model was tuned to best-fit the hydrograph rising limb and timing of peak discharge in both events; it achieved Mean Absolute Errors of

©2016. American Geophysical Union. All Rights Reserved. Surface meltwater that reaches the base of the Greenland Ice Sheet exerts a fundamental impact on ice flow, but observations of catchment-wide movement and distribution of subglacial water remain limited. Using radar-sounding data from two seasons, we identify the seasonal distribution of subglacial water in western Greenland. Our analysis provides evidence of widespread subglacial water storage beneath Greenland in the wintertime. The winter storage is located primarily on bedrock ridges with higher bed elevations in excess of 200 m. During the melt season water moves to the subglacial troughs. This inverse relationship with topography indicates that the material properties of the glacier bed strongly influence subglacial drainage development. Both the spatial variations in bed properties and the initial state of the subglacial hydrology system at the start of the melt season lead to differing glacier dynamical responses to surface melting across the Greenland Ice Sheet.

Accurate precipitation maps are essential for ecological, environmental, element cycle and hydrological models that have a spatial output component. It is well known that topography has a major influence on the spatial distribution of precipitation and that increasing topographical complexity is associated with increased spatial heterogeneity in precipitation. This means that when mapping precipitation using classical interpolation techniques (e.g. regression, kriging, spline, inverse distance weighting, etc.), a climate measuring network with higher spatial density is needed in mountainous areas in order to obtain the same level of accuracy as compared to flatter regions. In this study, we present a mean total annual precipitation mapping technique that combines topographical information (i.e. elevation and slope orientation) with average total annual rain gauge data in order to overcome this problem. A unique feature of this paper is the identification of the scale at which topography influences the precipitation pattern as well as the direction of the dominant weather circulation. This method was applied for Belgium and surroundings and shows that the identification of the appropriate scale at which topographical obstacles impact precipitation is crucial in order to obtain reliable mean total annual precipitation maps. The dominant weather circulation is determined at 260°. Hence, this approach allows accurate mapping of mean annual precipitation patterns in regions characterized by rather high topographical complexity using a climate data network with a relatively low density and/or when more advanced precipitation measurement techniques, such as radar, aren't available, for example in the case of historical data.

Meltwater delivered to the bed of the Greenland Ice Sheet is a driver of variable ice-motion through changes in effective pressure and enhanced basal lubrication. Ice surface velocities have been shown to respond rapidly both to meltwater production at the surface and to drainage of supraglacial lakes, suggesting efficient transfer of meltwater from the supraglacial to subglacial hydrological systems. Although considerable effort is currently being directed towards improved modelling of the controlling surface and basal processes, modelling the temporal and spatial evolution of the transfer of melt to the bed has received less attention. Here we present the results of spatially distributed modelling for prediction of moulins and lake drainages on the Leverett Glacier in Southwest Greenland. The model is run for the 2009 and 2010 ablation seasons, and for future increased melt scenarios. The temporal pattern of modelled lake drainages are qualitatively comparable with those documented from analyses of repeat satellite imagery. The modelled timings and locations of delivery of meltwater to the bed also match well with observed temporal and spatial patterns of ice surface speed-ups. This is particularly true for the lower catchment (

Abstract. Meltwater delivered to the bed of the Greenland Ice Sheet is a driver of variable ice-motion through changes in effective pressure and enhanced basal lubrication. Ice surface velocities have been shown to respond rapidly both to meltwater production at the surface and to drainage of supraglacial lakes, suggesting efficient transfer of meltwater from the supraglacial to subglacial hydrological systems. Although considerable effort is currently being directed towards improved modelling of the controlling surface and basal processes, modelling the temporal and spatial evolution of the transfer of melt to the bed has received less attention. Here we present the results of spatially-distributed modelling for prediction of moulins and lake drainages on the Leverett Glacier in south-west Greenland. The model is run for the 2009 and 2010 ablation seasons, and for future increased melt scenarios. The temporal and spatial patterns of modelled lake drainages are qualitatively comparable with those seen from analyses of satellite imagery. The modelled timings and locations of delivery of meltwater to the bed match well with observed temporal and spatial patterns of ice surface speed ups. This is particularly true for the lower catchment (&lt; 1000 m a.s.l.) where both the model and observations indicate that the development of moulins is the main mechanism for the transfer of surface meltwater to the bed. At higher elevations (e.g. 1250–1500 m a.s.l.) the development and drainage of supraglacial lakes becomes increasingly important. At these higher elevations, the delay between modelled melt generation and subsequent delivery of melt to the bed matches the observed delay between the peak air temperatures and subsequent velocity speed ups. Although both moulins and lake drainages are predicted to increase in number for future warmer climate scenarios, the lake drainages play an increasingly important role in both expanding the area over which melt accesses the bed and in enabling a greater proportion of surface melt to reach the bed.
.

We present a new bed elevation dataset for Greenland derived from a combination of multiple airborne ice thickness surveys undertaken between the 1970s and 2012. Around 420 000 line kilometres of airborne data were used, with roughly 70% of this having been collected since the year 2000, when the last comprehensive compilation was undertaken. The airborne data were combined with satellite-derived elevations for non-glaciated terrain to produce a consistent bed digital elevation model (DEM) over the entire island including across the glaciated-ice free boundary. The DEM was extended to the continental margin with the aid of bathymetric data, primarily from a compilation for the Arctic. Ice thickness was determined where an ice shelf exists from a combination of surface elevation and radar soundings. The across-track spacing between flight lines warranted interpolation at 1 km postings for significant sectors of the ice sheet. Grids of ice surface elevation, error estimates for the DEM, ice thickness and data sampling density were also produced alongside a mask of land/ocean/grounded ice/floating ice. Errors in bed elevation range from a minimum of ±10 m to about ±300 m, as a function of distance from an observation and local topographic variability. A comparison with the compilation published in 2001 highlights the improvement in resolution afforded by the new datasets, particularly along the ice sheet margin, where ice velocity is highest and changes in ice dynamics most marked. We estimate that the volume of ice included in our land-ice mask would raise mean sea level by 7.36 m, excluding any solid earth effects that would take place during ice sheet decay. © 2013 Author(s).

We present a new method of modelling the growth of supraglacial lakes at the western margin of the Greenland ice sheet, based on routing runoff estimated by a regional climate model across a digital elevation model (DEM) of the ice sheet surface. Using data acquired during the 2003 melt season, we demonstrate that the model is 19 times more likely to correctly predict the presence (or absence) of lakes than it is to make incorrect predictions, within an elevation range of 1100 to 1700 metres above sea level (m a.s.l.), when compared with MODIS satellite imagery. of the 66% of observed lake locations which the model correctly reproduces, the simulated lake onset day is found to be correlated with that observed with a Pearson correlation coefficient of 0.76. Our model accurately simulates maximum cumulative lake area with only a 1.5% overestimate. However, because our model does not simulate processes leading to lake stagnation or decay, such as refreezing or drainage, at present we do not simulate absolute daily lake area. We find that the maximum potential lake-covered ice sheet area is limited by topography to 6.4%. We estimate that this corresponds to a volume of 1.49 km3, 12% of the runoff produced in 2003. This can be taken as an upper bound given uncertainty in the DEM. This study has proved a good first step towards capturing the variability of supraglacial lake evolution with a numerical model. These initial results are promising and suggest that the model is a useful tool for use in analysing the behaviour of supraglacial lakes on the Greenland ice sheet in the present day and potentially beyond. ©2012 Author(s).

We have used Radar Altimeter 2 (RA-2) onboard ESA's EnviSAT and Geosciences Laser Altimeter System (GLAS) onboard NASA's ICESat to map the elevation change of the Flade Isblink Ice Cap (FIIC) in northern Greenland. Based on RA-2 data we show that the mean surface elevation change of the FIIC has been near zero (0.03±0.03 m/a) between fall 2002 and fall 2009. We present the elevation change rate maps and assess the elevation change rates of areas above the late summer snow line (0.09±0.04 m/a) and below it (-0.16±0.05 m/a). The GLAS elevation change rate maps show that some outlet glaciers, previously reported to have been in a surge state, are thickening rapidly. Using the RA-2 measured average elevation change rates for different parts of the ice cap we present a mass change rate estimate of 0.0±0.5 Gt/a for the FIIC. We compare the annual elevation changes with surface mass balance (SMB) estimates from a regional atmospheric climate model RACMO2. We find a strong correlation between the two (R = 0.94 and P < 0.002), suggesting that the surface elevation changes of the FIIC are mainly driven by net SMB. The correlation of modeled net SMB and measured elevation change is strong in the southern areas of the FIIC (R = 0.97 and P < 0.0005), but insignificant in the northern areas (R = 0.38 and P = 0.40). This is likely due to higher variability of glacier flow in the north relative to the south. Copyright 2011 by the American Geophysical Union.

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. © 2010 Elsevier B.V.

We present global positioning system observations that capture the full inland extent of ice motion variations in 2009 along a transect in the west Greenland Ice sheet margin. In situ measurements of air temperature and surface ablation, and satellite monitoring of ice surface albedo and supraglacial lake drainage are used to investigate hydrological controls on ice velocity changes. We find a strong positive correlation between rates of annual ablation and changes in annual ice motion along the transect, with sites nearest the ice sheet margin experiencing greater annual variations in ice motion (15-18%) than those above 1000. m elevation (3-8%). Patterns in the timing and rate of meltwater delivery to the ice-bed interface provide key controls on the magnitude of hydrologically-forced velocity variations at each site. In the lower ablation zone, the overall contribution of variations in ice motion to annual flow rates is limited by evolution in the structure of the subglacial drainage system. At sites in the upper ablation zone, a shorter period of summer melting and delayed establishment of a hydraulic connection between the ice sheet surface and its bed limit the timeframe for velocity variations to occur. Our data suggest that land-terminating sections of the Greenland Ice Sheet will experience increased dynamic mass loss in a warmer climate, as the behaviour that we observe in the lower ablation zone propagates further inland. Findings from this study provide a conceptual framework to understand the impact of hydrologically-forced velocity variations on the future mass balance of land-terminating sections of the Greenland Ice Sheet. © 2011 Elsevier B.V.

We measure hydrological parameters in meltwater draining from an outlet glacier in west Greenland to investigate seasonal changes in the structure and behaviour of the hydrological system of a large catchment in the Greenland ice sheet (GrIS). Our data reveal seasonal upglacier expansion and increase in hydraulic efficiency of the subglacial drainage system, across a catchment >600 km2, to distances >50 km from the ice-sheet margin. This expansion occurs episodically in response to the drainage of surface meltwaters into a hitherto inefficient subglacial drainage system as new input locations become active progressively further upglacier; this system is similar to Alpine glaciers. These observations provide the first synopsis of seasonal hydrological behaviour in the ablation zone of the GrIS. Copyright 2011 by the American Geophysical Union.

The 8500 km2 Flade Isblink ice cap (FIIC) (8115N, 150W) is the largest ice cap in Greenland. We use repeat-pass interferometric synthetic aperture radar (InSAR) techniques to investigate the form and flow of the FIIC. European Remote Sensing satellite (ERS-1 and ERS-2) data acquired in winter 1996 were used to form a 100 m resolution digital elevation model (DEM), which we constrained using Ice Cloud and Elevation satellite (ICESat) laser altimeter elevation measurements from 2007. This InSAR DEM was used to isolate the phase due to motion from seven ERS-tandem (1 day) pairs of SAR scenes acquired between 15 August 1995 and 7 February 1996, to produce one wintertime and two summertime velocity maps. Five of the eight major outlet glaciers draining the FIIC are marine terminating, and two terminate at a lake margin. A maximum ice velocity of 581 m yr-1 was observed in mid-August 1995. Six of the eight major outlet glaciers exhibit seasonal velocity variations between late summer and winter, and flow speeds vary by up to 20% over a 10 day period in August 1995. Our findings show that while marine terminating glaciers flow faster than land terminating glaciers, there is no simple relationship between glacier type and seasonality of ice motion. Copyright 2010 by the American Geophysical Union.