Mass balance of the Greenland ice sheet (2003-2008) from ICESat data:The impact of interpolation, sampling and firn density

ICESat has provided surface elevation measurements of the ice sheets since the launch in January 2003, resulting in a unique dataset for monitoring the changes of the cryosphere. Here, we present a novel method for determining the mass balance of the Greenland ice sheet, derived from ICESat altimetry data. <br><br> Three different methods for deriving elevation changes from the ICESat altimetry dataset are used. This multi-method approach provides a method to assess the complexity of deriving elevation changes from this dataset. <br><br> The altimetry alone can not provide an estimate of the mass balance of the Greenland ice sheet. Firn dynamics and surface densities are important factors that contribute to the mass change derived from remote-sensing altimetry. The volume change derived from ICESat data is corrected for changes in firn compaction over the observation period, vertical bedrock movement and an intercampaign elevation bias in the ICESat data. Subsequently, the corrected volume change is converted into mass change by the application of a simple surface density model, in which some of the ice dynamics are accounted for. The firn compaction and density models are driven by the HIRHAM5 regional climate model, forced by the ERA-Interim re-analysis product, at the lateral boundaries. <br><br> We find annual mass loss estimates of the Greenland ice sheet in the range of 191 ± 23 Gt yr<sup>−1</sup> to 240 ± 28 Gt yr<sup>−1</sup> for the period October 2003 to March 2008. These results are in good agreement with several other studies of the Greenland ice sheet mass balance, based on different remote-sensing techniques

Towards Constraining Glacial Isostatic Adjustment in Greenland Using ICESat and GPS Observations

Constraining glacial isostatic adjustment (GIA) i.e. the Earth’s viscoelastic response to past ice changes, is an important task, because GIA is a significant correction in gravity-based ice sheet mass balance estimates. Here, we investigate how temporal variations in the observed and modeled crustal displacements due to the Earth’s response to ongoing ice mass changes can contribute to the process of constraining GIA. We use mass change grids of the Greenland ice sheet (GrIS) derived from NASA’s high resolution Ice, Cloud and land Elevation Satellite (ICESat) data in three overlapping time spans covering the period 2004–2009 to estimate temporal variations in the elastic response due to present day ice mass loss. The modeled crustal displacements (elastic + GIA) are compared with GPS time series from five permanent sites (KELY, KULU, QAQ1, THU2, and SCOR). We find, that the modeled pattern of elastic crustal displacements shows pronounced variation during the observation period, where an increase in elastic displacement is found at the northwest coast of Greenland, while a decrease is found at the southeast coast. This pattern of temporal changes is supported by the GPS observations. We find, that the temporal behavior of the ICESat-based modeled elastic response agrees well with the GPS observations at the sites KELY, QAQ1, and SCOR. This suggests, that our elastic models are able to resolve the temporal changes in the observed uplift, which indicates that the elastic uplift models are reliable at these sites. Therefore, we conclude that these sites are useful for constraining GIA

Use of expert elicitation to assign weights to climate and hydrological models in climate impact studies

Various methods are available for assessing uncertainties in climate impact studies. Among such methods, model weighting by expert elicitation is a practical way to provide a weighted ensemble of models for specific real-world impacts. The aim is to decrease the influence of improbable models in the results and easing the decision-making process. In this study both climate and hydrological models are analysed, and the result of a research experiment is presented using model weighting with the participation of six climate model experts and six hydrological model experts. For the experiment, seven climate models are a priori selected from a larger EURO-CORDEX (Coordinated Regional Downscaling Experiment - European Domain) ensemble of climate models, and three different hydrological models are chosen for each of the three European river basins. The model weighting is based on qualitative evaluation by the experts for each of the selected models based on a training material that describes the overall model structure and literature about climate models and the performance of hydrological models for the present period. The expert elicitation process follows a three-stage approach, with two individual rounds of elicitation of probabilities and a final group consensus, where the experts are separated into two different community groups: a climate and a hydrological modeller group. The dialogue reveals that under the conditions of the study, most climate modellers prefer the equal weighting of ensemble members, whereas hydrological-impact modellers in general are more open for assigning weights to different models in a multi-model ensemble, based on model performance and model structure. Climate experts are more open to exclude models, if obviously flawed, than to put weights on selected models in a relatively small ensemble. The study shows that expert elicitation can be an efficient way to assign weights to different hydrological models and thereby reduce the uncertainty in climate impact. However, for the climate model ensemble, comprising seven models, the elicitation in the format of this study could only re-establish a uniform weight between climate models

The future potential for wine production in Scotland under high-end climate change

Wine production is climate dependent and highly sensitive to weather variability, which makes the sector a good indicator of ongoing and future climate change impacts. Under high-end climate change (HECC), temperatures in Scotland are projected to increase significantly by the end of the twenty-first century. This raises the possibility of future temperatures becoming sufficiently high to support the growth of wine grapes. In this paper, we explore to what extent Scotland might become suitable for wine production under HECC using a climate analogue approach. Specifically, we address the following questions. What are the projected late twenty-first century temperature changes in Scotland? Where in Europe are current climates (based on summer and annual temperatures) similar to those projected for Scotland by the end of the twenty-first century under HECC? Are any of these locations currently wine grape growing regions? The temperature analogues towards the end of the twenty-first century occurred at more southerly latitudes in Europe, with some variability from west to east arising from the influence of continental climates. Temperature analogues alone match with several current wine grape growing regions of Europe, suggesting that future climates in Scotland could support wine production. However, when precipitation and/or lithology and topography are also taken into account, no matches were found with existing European wine grape growing regions. This study demonstrates how the use of climate analogues in combination with other environmental datasets can be useful in understanding future climate change impacts, especially under HECC

Characterizing, modelling and understanding the climate variability of the deep water formation in the North-Western Mediterranean Sea

Observing, modelling and understanding the climate-scale variability of the deep water formation (DWF) in the North-Western Mediterranean Sea remains today very challenging. In this study, we first characterize the interannual variability of this phenomenon by a thorough reanalysis of observations in order to establish reference time series. These quantitative indicators include 31 observed years for the yearly maximum mixed layer depth over the period 1980–2013 and a detailed multi-indicator description of the period 2007–2013. Then a 1980–2013 hindcast simulation is performed with a fully-coupled regional climate system model including the high-resolution representation of the regional atmosphere, ocean, land-surface and rivers. The simulation reproduces quantitatively well the mean behaviour and the large interannual variability of the DWF phenomenon. The model shows convection deeper than 1000 m in 2/3 of the modelled winters, a mean DWF rate equal to 0.35 Sv with maximum values of 1.7 (resp. 1.6) Sv in 2013 (resp. 2005). Using the model results, the winter-integrated buoyancy loss over the Gulf of Lions is identified as the primary driving factor of the DWF interannual variability and explains, alone, around 50 % of its variance. It is itself explained by the occurrence of few stormy days during winter. At daily scale, the Atlantic ridge weather regime is identified as favourable to strong buoyancy losses and therefore DWF, whereas the positive phase of the North Atlantic oscillation is unfavourable. The driving role of the vertical stratification in autumn, a measure of the water column inhibition to mixing, has also been analyzed. Combining both driving factors allows to explain more than 70 % of the interannual variance of the phenomenon and in particular the occurrence of the five strongest convective years of the model (1981, 1999, 2005, 2009, 2013). The model simulates qualitatively well the trends in the deep waters (warming, saltening, increase in the dense water volume, increase in the bottom water density) despite an underestimation of the salinity and density trends. These deep trends come from a heat and salt accumulation during the 1980s and the 1990s in the surface and intermediate layers of the Gulf of Lions before being transferred stepwise towards the deep layers when very convective years occur in 1999 and later. The salinity increase in the near Atlantic Ocean surface layers seems to be the external forcing that finally leads to these deep trends. In the future, our results may allow to better understand the behaviour of the DWF phenomenon in Mediterranean Sea simulations in hindcast, forecast, reanalysis or future climate change scenario modes. The robustness of the obtained results must be however confirmed in multi-model studies

The future sea-level rise contribution of Greenland’s glaciers and ice caps

We calculate the future sea-level rise contribution from the surface mass balance of all of Greenland's glaciers and ice caps (GICs, ~90 000 km2) using a simplified energy balance model which is driven by three future climate scenarios from the regional climate models HIRHAM5, RACMO2 and MAR. Glacier extent and surface elevation are modified during the mass balance model runs according to a glacier retreat parameterization. Mass balance and glacier surface change are both calculated on a 250 m resolution digital elevation model yielding a high level of detail and ensuring that important feedback mechanisms are considered. The mass loss of all GICs by 2098 is calculated to be 2016 ± 129 Gt (HIRHAM5 forcing), 2584 ± 109 Gt (RACMO2) and 3907 ± 108 Gt (MAR). This corresponds to a total contribution to sea-level rise of 5.8 ± 0.4, 7.4 ± 0.3 and 11.2 ± 0.3 mm, respectively. Sensitivity experiments suggest that mass loss could be higher by 20–30% if a strong lowering of the surface albedo were to take place in the future. It is shown that the sea-level rise contribution from the north-easterly regions of Greenland is reduced by increasing precipitation while mass loss in the southern half of Greenland is dominated by steadily decreasing summer mass balances. In addition we observe glaciers in the north-eastern part of Greenland changing their characteristics towards greater activity and mass turnover