Search

Since the year 2000, Greenland ice sheet mass loss has been dominated by a decrease in surface mass balance rather than an increase in solid ice discharge. Southeast Greenland is an important region to understand how high accumulation rates can offset increasing Greenland ice sheet meltwater runoff. To that end, we derive a new 9-year long dataset (2009–17) of accumulation rates in Southeast Greenland using NASA Operation IceBridge snow radar. Our accumulation dataset derived from internal layers focuses on high elevations (1500–3000 m) because at lower elevations meltwater percolation obscured internal layer structure. The uncertainty of the radar-derived accumulation rates is 11% [using Firn Densification Model (FDM) density profiles] and the average accumulation rate ranges from 0.5 to 1.2 m w.e. With our observations spanning almost a decade, we find large inter-annual variability, but no significant trend. Accumulation rates are compared with output from two regional climate models (RCMs), MAR and RACMO2. This comparison shows that the models are underestimating accumulation in Southeast Greenland and the models misrepresent spatial heterogeneity due to an orographically forced bias in snowfall near the coast. Our dataset is useful to fill in temporal and spatial data gaps, and to evaluate RCMs where few in situ measurements are available.

The high Arctic is the fastest warming region on Earth, evidenced by extreme near-surface temperature increase in non-summer seasons, recent rapid sea ice decline and permafrost melting since the early 1990’s. Understanding the impact of climate change on the sensitive Arctic ecosystem to climate change has so far been hampered by the lack of time-constrained, high-resolution records and by implicit climate data analyses. Here, we show evidence of sharp growth in freshwater green algae as well as distinct diatom assemblage changes since ~1995, retrieved from a high-Arctic (80 °N) lake sediment record on Barentsøya (Svalbard). The proxy record approaches an annual to biennial resolution. Combining remote sensing and in-situ climate data, we show that this ecological change is concurrent with, and is likely driven by, the atmospheric warming and a sharp decrease in the length of the sea ice covered period in the region, and throughout the Arctic. Moreover, this research demonstrates the value of palaeoclimate records in pristine environments for supporting and extending instrumental records. Our results reinforce and extend observations from other sites that the high Arctic has already undergone rapid ecological changes in response to on-going climate change, and will continue to do so in the future.

In Earth system models, terrestrial snow is usually modeled by the land surface component. In most cases, these snow models have been developed with an emphasis on seasonal snow. Questions about future sea level rise, however, prompt the need for a realistic representation of perennial snow, as snow processes play a key role in the mass balance of glaciers and ice sheets. Here we enhance realism of modeled polar snow in the Community Land Model (CLM), the land component of the Community Earth System Model (CESM), by implementing (1) new parametrizations for fresh snow density, destructive metamorphism, and compaction by overburden pressure, (2) by allowing for deeper snow packs, and (3) by introducing drifting snow compaction, with a focus on the ice sheet interior. Comparison with Greenlandic and Antarctic snow density observations show that the new physics improve model skill in predicting firn and near-surface density in the absence of melt. Moreover, compensating biases are removed and spurious subsurface melt rates at ice sheets are eliminated. The deeper snow pack enhances refreezing and allows for deeper percolation, raising ice temperatures up to 15 ºC above the skin temperature.

Latest Maastrichtian climate change caused by Deccan volcanism has been invoked as a cause of mass extinction at the Cretaceous-Paleogene (K-Pg) boundary (~66.0 Ma). Yet late Maastrichtian climate and ecological changes are poorly documented, in particular on the Southern Hemisphere. Here we present upper Maastrichtian-lower Danian climate and biotic records from the Bajada del Jagüel (BJ) shelf site (Neuquén Basin, Argentina), employing the TEX86 paleothermometer, marine palynology (dinoflagellate cysts), and micropaleontology (foraminifera). These records are correlated to the astronomically tuned Ocean Drilling Program Site 1262 (Walvis Ridge). Collectively, we use these records to assess climatic and ecological effects of Deccan volcanism in the Southern Atlantic region. Both the TEX86-based sea surface temperature (SST) record at BJ and the bulk carbonate δ18O-based SST record of Site 1262 show a latest Maastrichtian warming of ~2.5-4°C, at 450 to 150 kyr before the K-Pg boundary, coinciding with the a large Deccan outpouring phase. Benthic foraminiferal and dinocyst assemblage changes indicate that this warming resulted in enhanced runoff and stratification of the water column, likely resulting from more humid climate conditions in the Neuquén Basin. These climate conditions could have been caused by an expanding and strengthening thermal low over the South American continent. Biotic changes in response to late Maastrichtian environmental changes are rather limited, when compared to the major turnovers observed at many K-Pg boundary sites worldwide. This suggests that environmental perturbations during the latest Maastrichtian warming event were less severe than those following the K-Pg boundary impact.

Clouds play a pivotal role in the surface energy budget of the polar regions. Here we use two largely independent data sets of cloud and surface downwelling radiation observations derived by satellite remote sensing (2007–2010) to evaluate simulated clouds and radiation over both polar ice sheets and oceans in state-of-the-art atmospheric reanalyses (ERA-Interim and Modern Era Retrospective-Analysis for Research and Applications-2) and the Coupled Model Intercomparison Project Phase 5 (CMIP5) climate model ensemble. First, we show that, compared to Clouds and the Earth's Radiant Energy System-Energy Balanced and Filled, CloudSat-CALIPSO better represents cloud liquid and ice water path over high latitudes, owing to its recent explicit determination of cloud phase that will be part of its new R05 release. The reanalyses and climate models disagree widely on the amount of cloud liquid and ice in the polar regions. Compared to the observations, we find significant but inconsistent biases in the model simulations of cloud liquid and ice water, as well as in the downwelling radiation components. The CMIP5 models display a wide range of cloud characteristics of the polar regions, especially with regard to cloud liquid water, limiting the representativeness of the multimodel mean. A few CMIP5 models (CNRM, GISS, GFDL, and IPSL_CM5b) clearly outperform the others, which enhances credibility in their projected future cloud and radiation changes over high latitudes. Given the rapid changes in polar regions and global feedbacks involved, future climate model developments should target improved representation of polar clouds. To that end, remote sensing observations are crucial, in spite of large remaining observational uncertainties, which is evidenced by the substantial differences between the two data sets.

This study uses output of a high-resolution (5.5 km) regional atmospheric climate model to describe the present-day (1979–2012) climate of Patagonia, with a particular focus on the surface mass balance (SMB) of the Patagonian ice fields. Through a comparison with available in situ observations, it is shown that the model is able to simulate the sharp climate gradients in western Patagonia. The southern Andes are an efficient barrier for the prevalent atmospheric flow, generating strong orographic uplift and precipitation throughout the entire year. The model suggests extreme orographic precipitation west of the Andes divide, with annual precipitation rates of >5 to 34 m w.e. (water equivalent), and a clear rain shadow east of the divide. These modeled precipitation rates are supported qualitatively by available precipitation stations and SMB estimates on the ice fields derived from firn cores. For the period 1979–2012, a slight atmospheric cooling at upper ice field elevations is found, leading to a small but insignificant increase in the ice field SMB.

This paper presents autonomous drifting snow observations performed on the Greenland Ice Sheet in the fall of 2012. High-frequency snow particle counter (SPC) observations at ~ 1 m above the surface provided drifting snow number fluxes and size distributions; these were combined with meteorological observations at six levels. We identify two types of drifting snow events: katabatic events are relatively cold and dry, with prevalent winds from the southeast, whereas synoptic events are short lived, warm and wet. Precipitating snow during synoptic events disturbs the drifting snow measurements. Output of the regional atmospheric climate model RACMO2, which includes the drifting snow routine PIEKTUK-B, agrees well with the observed near-surface climate at the site, as well as with the frequency and timing of drifting snow events. Direct comparisons with the SPC observations at 1 m reveal that the model overestimates the horizontal snow transport at this level, which can be related to an overestimation of saltation and the typical size of drifting snow particles.

In this study the effects of changes in the physics package of the regional atmospheric climate model RACMO2 on the modelled surface energy balance, near-surface temperature and wind speed of Antarctica are presented. The physics package update primarily consists of an improved turbulent and radiative flux scheme and a revised cloud scheme that includes a parameterisation for ice cloud super-saturation. The ice cloud super-saturation has led to more moisture being transported onto the continent, resulting in more and optically thicker clouds and more downward long-wave radiation. Overall, the updated model better represents the surface energy balance, based on a comparison with >750 months of data from nine automatic weather stations located in East Antarctica. Especially the representation of the turbulent sensible heat flux and net long-wave radiative flux has improved with a decrease in biases of up to 40%. As a result, modelled surface temperatures have increased and the bias, when compared to 10 m snow temperatures from 64 ice-core observations, has decreased from −2.3 K to −1.3 K. The weaker surface temperature inversion consequently improves the representation of the sensible heat flux, whereas wind speed biases remain unchanged. However, significant model biases remain, partly because RACMO2 at a resolution of 27 km is unable to resolve steep topography.

Climate change west of the Antarctic Peninsula is the most rapid of anywhere in the Southern Hemisphere, with associated changes in the rates and distributions of freshwater inputs to the ocean. Here, results from the first comprehensive survey of oxygen isotopes in seawater in this region are used to quantify spatial patterns of meteoric water (glacial discharge and precipitation) separately from sea ice melt. High levels of meteoric water are found close to the coast, due to orographic effects on precipitation and strong glacial discharge. Concentrations decrease offshore, driving significant southward geostrophic flows (up to ~30 cm s−1). These produce high meteoric water concentrations at the southern end of the sampling grid, where collapse of the Wilkins Ice Shelf may also have contributed. Sea ice melt concentrations are lower than meteoric water and patchier because of the mobile nature of the sea ice itself. Nonetheless, net sea ice production in the northern part of the sampling grid is inferred; combined with net sea ice melt in the south, this indicates an overall southward ice motion. The survey is contextualized temporally using a decade-long series of isotope data from a coastal Antarctic Peninsula site. This shows a temporal decline in meteoric water in the upper ocean, contrary to expectations based on increasing precipitation and accelerating deglaciation. This is driven by the increasing occurrence of deeper winter mixed layers and has potential implications for concentrations of trace metals supplied to the euphotic zone by glacial discharge. As the regional freshwater system evolves, the continuing isotope monitoring described here will elucidate the ongoing impacts on climate and the ecosystem.

This paper presents the drifting snow climate of the Greenland ice sheet, using output from a high-resolution (∼11 km) regional climate model. Because reliable direct observations of drifting snow do not exist, we evaluate the modeled near-surface climate instead, using automatic weather station (AWS) observations from the K-transect and find that RACMO2 realistically simulates near-surface wind speed and relative humidity, two variables that are important for drifting snow. Integrated over the ice sheet, drifting snow sublimation (SUds) equals 24 ± 3 Gt yr−1, and is significantly larger than surface sublimation (SUs, 16 ± 2 Gt yr−1). SUds strongly varies between seasons, and is only important in winter, when surface sublimation and runoff are small. A rapid transition exists between the winter season, when snowfall and SUds are important, and the summer season, when snowmelt is significant, which increases surface snow density and thereby limits drifting snow processes. Drifting snow erosion (ERds) is only important on a regional scale. In recent decades, following decreasing wind speed and rising near-surface temperatures, SUds exhibits a negative trend (0.1 ± 0.1 Gt yr−1), which is compensated by an increase in SUs of similar magnitude.

The GLACIOCLIM-SAMBA (GS) Antarctic accumulation monitoring network, which extends from the coast of Adelie Land to the Antarctic plateau, has been surveyed annually since 2004. The network includes a 156-km stake-line from the coast inland, along which accumulation shows high spatial and interannual variability with a mean value of 362 mm water equivalent a -1. In this paper, this accumulation is compared with older accumulation reports from between 1971 and 1991. The mean and annual standard deviation and the km-scale spatial pattern of accumulation were seen to be very similar in the older and more recent data. The data did not reveal any significant accumulation trend over the last 40 years. The ECMWF analysis-based forecasts (ERA-40 and ERA-Interim), a stretched-grid global general circulation model (LMDZ4) and three regional circulation models (PMM5, MAR and RACMO2), all with high resolution over Antarctica (27-125 km), were tested against the GS reports. They qualitatively reproduced the meso-scale spatial pattern of the annual-mean accumulation except MAR. MAR significantly underestimated mean accumulation, while LMDZ4 and RACMO2 overestimated it. ERA-40 and the regional models that use ERA-40 as lateral boundary condition qualitatively reproduced the chronology of interannual variability but underestimated the magnitude of interannual variations. Two widely used climatologies for Antarctic accumulation agreed well with the mean GS data. The model-based climatology was also able to reproduce the observed spatial pattern. These data thus provide new stringent constraints on models and other large-scale evaluations of the Antarctic accumulation.

Here we present a long-term (1850–2200) best estimate of Greenland ice sheet (GrIS) freshwater runoff that improves spatial detail of runoff locations and temporal resolution. Ice discharge is taken from observations since 2000 and assumed constant in time. Surface meltwater runoff is retrieved from regional climate model output for the recent past and parameterized for the future based on significant correlations between runoff and midtropospheric (500 hPa) summer temperature changes over the GrIS. The simplicity of this approach enables assimilation of meltwater runoff into coupled climate models, which is demonstrated here in a case study with the medium-resolution (1?) Community Earth System Model. The model results suggest that the decrease in Atlantic Meridional Overturning Circulation (AMOC) is dominated by warming of the surface ocean and enhanced GrIS freshwater forcing leads to a slightly enhanced (?1.2 sverdrup in the 21st century) weakening of the AMOC.

We assess the runoff and surface mass balance (SMB) of the Greenland ice sheet in the Nuuk region (southwest) using output of two regional climate models (RCMs) evaluated by observations. The region encompasses six glaciers that drain into Godthåbsfjord. RCM data (1960–2012) are resampled to a high spatial resolution to include the narrow (relative to the native grid spacing) glacier trunks in the ice mask. Comparing RCM gridded results with automatic weather station (AWS) point measurements reveals that locally models can underestimate ablation and overestimate accumulation by up to tens of per cent. However, comparison with lake discharge indicates that modelled regional runoff totals are more accurate. Model results show that melt and runoff in the Nuuk region have doubled over the past two decades. Regional SMB attained negative values in recent high-melt years. Taking into account frontal ablation of the marine-terminating glaciers, the region lost 10–20 km3 w.e. a–1 in 2010–12. If 2010 melting prevails during the remainder of this century, a low-end estimate of sea-level rise of 5 mm is expected by 2100 from this relatively small section (2.6%) of the ice sheet alone.

This study presents the drifting snow climate of the Earth's ice sheets, Antarctica and Greenland. For that purpose we use a regional atmospheric climate model, RACMO2. We included a routine that is able to calculate the drifting snow fluxes and accounts for the interaction between drifting snow on the one hand and the atmosphere and snow surface on the other. RACMO2 is run at 27 km resolution for Antarctica, and 11 km resolution for Greenland, and forced at its lateral boundaries by ECMWF reanalyses (32 years for Antarctica and 52 years for Greenland). Because direct evaluation for drifting snow is challenging due to sparseness of observational data, we focussed the model evaluation on the ability of RACMO2 to represent near-surface wind climate, temperature, surface mass balance, the extent of ablation areas and remote-sensed drifting snow frequency. We show that RACMO2 is very well able to represent the present-day near-surface climate of Antarctica and Greenland. Drifting snow occurs 20-80% of the time on Antarctica, depending on the local wind climate. Highest frequencies are found in the coastal areas, where drifting snow sublimation (SUds) removes up to 150 mm water equivalent of snow, whereas the high-elevation areas experience little or no SUds. Drifting snow erosion (ER­ds) can be negative (deposition) or positive (erosion), and varies generally between -50 and 50 mm in regions where the wind field convergences and diverges, respectively. Integrated over the ice sheet, SUds removes around 165 Gt of snow, which is equivalent to ~6% of the precipitated snow. The impact of ER­ds on the Antarctic ice sheet SMB is negligible . We found several feedbacks between SUds and the atmosphere. SUds moistens the near-surface atmosphere, limiting its own potential, but also enhancing precipitation in some coastal areas. By removing mass from the snow surface, drifting snow processes increase the top snow layer density, increasing the threshold wind speed for further drifting snow. Since the impacts of drifting snow are included, we are able to quantify all contributors to the Antarctic SMB. The SMB equals 2479 +/- 147 Gt/yr, and is largely determined by snowfall. It varies from

The Canadian Arctic Archipelago (CAA) contains the largest volume of glacier ice on Earth outside of Antarctica and Greenland. In the absence of significant calving, CAA glacier mass balance is governed by the difference between surface snow accumulation and meltwater runoff—surface mass balance. Here we use a coupled atmosphere/snow model to simulate present-day and 21st century CAA glacier surface mass balance. Through comparison with Gravity Recovery and Climate Experiment mass anomalies and in situ observations, we show that the model is capable of representing present-day CAA glacier mass loss, as well as the dynamics of the seasonal snow cover on the CAA tundra. Next, we force this model until 2100 with a moderate climate warming scenario (AR5 RCP4.5). We show that enhanced meltwater runoff from CAA glaciers is not sufficiently compensated by increased snowfall. Extrapolation of these results toward an AR5 multimodel ensemble results in sustained 21st century CAA glacier mass loss in the vast majority (>99%) of the ~7000 temperature realizations.

Enhanced snowfall on the East Antarctic ice sheet is projected to significantly mitigate 21st century global sea level rise. In recent years (2009 and 2011), regionally extreme snowfall anomalies in Dronning Maud Land, in the Atlantic sector of East Antarctica, have been observed. It has been unclear, however, whether these anomalies can be ascribed to natural decadal variability, or whether they could signal the beginning of a long-term increase of snowfall. Here we use output of a regional atmospheric climate model, evaluated with available firn core records and gravimetry observations, and show that such episodes had not been seen previously in the satellite climate data era (1979). Comparisons with historical data that originate from firn cores, one with records extending back to the 18th century, confirm that accumulation anomalies of this scale have not occurred in the past ~60 years, although comparable anomalies are found further back in time. We examined several regional climate model projections, describing various warming scenarios into the 21st century. Anomalies with magnitudes similar to the recently observed ones were not present in the model output for the current climate, but were found increasingly probable toward the end of the 21st century.

Mass loss from the two major ice sheets and their contribution to global sea level rise is accelerating. In Antarctica, mass loss is dominated by increased flow velocities of outlet glaciers, following the thinning or disintegration of coastal ice shelves into which they flow. In contrast, 55% of post-1992 Greenland ice sheet (GrIS) mass loss is accounted for by surface processes, notably increased meltwater runoff. A subtle process in the surface mass balance of the GrIS is the retention and refreezing of meltwater, currently preventing 40% of the meltwater to reach the ocean. Here we force a high-resolution atmosphere/snow model with a mid-range warming scenario (RCP4.5, 1970–2100), to show that rapid loss of firn pore space, by >50% at the end of the 21st century, quickly reduces this refreezing buffer. As a result, GrIS surface mass loss accelerates throughout the 21st century and its contribution to global sea level rise increases to 1.7˙0.5 mm yr–1, more than four times the current value.

We use an airborne-radar method, verified with ice-core accumulation records, to determine the spatiotemporal variations of snow accumulation over Thwaites Glacier, West Antarctica between 1980 and 2009. We also present a regional evaluation of modeled accumulation in Antarctica. Comparisons between radar-derived measurements and model outputs show that three global models capture the interannual variability well (r > 0.9), but a high-resolution regional model (RACMO2) has better absolute accuracy and captures the observed spatial variability (r = 0.86). Neither the measured nor modeled accumulation records over Thwaites Glacier show any trend since 1980. Although an increase in accumulation may potentially accompany the observed warming in the region, the projected trend is too small to detect over the 30 year record.

The Greenland and Antarctic ice sheets have been reported to be losing mass at accelerating rates1, 2. If sustained, this accelerating mass loss will result in a global mean sea-level rise by the year 2100 that is approximately 43 cm greater than if a linear trend is assumed2. However, at present there is no scientific consensus on whether these reported accelerations result from variability inherent to the ice-sheet–climate system, or reflect long-term changes and thus permit extrapolation to the future3. Here we compare mass loss trends and accelerations in satellite data collected between January 2003 and September 2012 from the Gravity Recovery and Climate Experiment to long-term mass balance time series from a regional surface mass balance model forced by re-analysis data. We find that the record length of spaceborne gravity observations is too short at present to meaningfully separate long-term accelerations from short-term ice sheet variability. We also find that the detection threshold of mass loss acceleration depends on record length: to detect an acceleration at an accuracy within ±10 Gt yr−2, a period of 10 years or more of observations is required for Antarctica and about 20 years for Greenland. Therefore, climate variability adds uncertainty to extrapolations of future mass loss and sea-level rise, underscoring the need for continuous long-term satellite monitoring.

To estimate the sea level rise (SLR) originating from changes in surface mass balance (SMB) of the Greenland ice sheet (GrIS), we present 21st century climate projections obtained with the regional climate model MAR (Modèle Atmosphérique Régional), forced by output of three CMIP5 (Coupled Model Intercomparison Project Phase 5) general circulation models (GCMs). Our results indicate that in a warmer climate, mass gain from increased winter snowfall over the GrIS does not compensate mass loss through increased meltwater run-off in summer. Despite the large spread in the projected near-surface warming, all the MAR projections show similar non-linear increase of GrIS surface melt volume because no change is projected in the general atmospheric circulation over Greenland. By coarsely estimating the GrIS SMB changes from GCM output, we show that the uncertainty from the GCM-based forcing represents about half of the projected SMB changes. In 2100, the CMIP5 ensemble mean projects a GrIS SMB decrease equivalent to a mean SLR of +4 ± 2 cm and +9 ± 4 cm for the RCP (Representative Concentration Pathways) 4.5 and RCP 8.5 scenarios respectively. These estimates do not consider the positive melt–elevation feedback, although sensitivity experiments using perturbed ice sheet topographies consistent with the projected SMB changes demonstrate that this is a significant feedback, and highlight the importance of coupling regional climate models to an ice sheet model. Such a coupling will allow the assessment of future response of both surface processes and ice-dynamic changes to rising temperatures, as well as their mutual feedbacks.

This study explores an approach that simultaneously estimates Antarctic mass balance and glacial isostatic adjustment (GIA) through the combination of satellite gravity and altimetry data sets. The results improve upon previous efforts by incorporating reprocessed data sets over a longer period of time, and now include a firn densification model to convert the altimetry volume estimates into mass. When the GIA models created from the combination approach were compared to insitu GPS ground station displacements, the vertical rates estimated showed good agreement after a systematic bias was removed from the computed GIA models. The new models suggest the potential for GIA uplift in the Amundsen Sea Sector, as well as the possible subsidence in large parts of East Antarctica.

To simulate the impact of drifting snow on the lower atmosphere, surface characteristics and surface mass balance (SMB) of the Antarctic ice sheet regional atmospheric climate model (RACMO2.1/ANT) with horizontal resolution of 27 km is coupled to a drifting snow routine and forced by ERA-Interim fields at its lateral boundaries (1989–2009). This paper evaluates the near-surface and drifting snow climate of RACMO2.1/ANT. Modeled near-surface wind speed (squared correlation coefficient R2 = 0.64) and temperature (R2 = 0.93) agree well with observations. Wind speed is underestimated in topographically complex areas, where observed wind speeds are locally very high (>20 m s!1). Temperature is underestimated in winter in coastal areas due to an underestimation of downward longwave radiation. Near-surface temperature and wind speed are not significantly affected by the inclusion of drifting snow in the model. In contrast, relative humidity with respect to ice increases in regions with strong drifting snow and becomes more consistent with the observations. Drifting snow frequency is the only observable parameter to directly validate drifting snow results; therefore, we derived an empirical relation for fresh snow density, as a function of wind speed and temperature, which determines the threshold wind speed for drifting snow. Modeled drifting snow frequencies agree well with in situ measurements and novel estimates from remote sensing. Finally, we show that including drifting snow is essential to obtaining a realistic extent and spatial distribution of ablation (SMB

Forty years of satellite imagery reveal that meltwater lakes on the margin of the Greenland Ice Sheet have expanded substantially inland to higher elevations with warming. These lakes are important because they provide a mechanism for bringing water to the ice bed, causing sliding. Inland expansion of lakes could accelerate ice flow by bringing water to previously frozen bed, potentially increasing future rates of mass loss. Increasing lake elevations closely follow the rise of the mass balance equilibrium line over much of the ice sheet, suggesting no physical limit on lake expansion. Data are not yet available to detect a corresponding change in ice flow, and the potential effects of lake expansion on ice sheet dynamics are not included in ice sheet models.

Surface snowmelt is widespread in coastal Antarctica. Satellite-based microwave sensors have been observing melt area and duration for over three decades. However, these observations do not reveal the total volume of meltwater produced on the ice sheet. Here we present an Antarctic melt volume climatology for the period 1979–2010, obtained using a regional climate model equipped with realistic snow physics. We find that mean continent-wide meltwater volume (1979–2010) amounts to 89 Gt y�1 with large interannual variability (s = 41 Gt y�1). Of this amount, 57 Gt y�1 (64%) is produced on the floating ice shelves extending from the grounded ice sheet, and 71 Gt y�1 in West-Antarctica, including the Antarctic Peninsula. We find no statistically significant trend in either continentwide or regional meltwater volume for the 31-year period 1979–2010. Citation: Kuipers Munneke, P., G. Picard, M. R. van den Broeke, J. T. M. Lenaerts, and E. van Meijgaard (2012), Insignificant change in Antarctic snowmelt volume since 1979, Geophys. Res. Lett., 39, L01501, doi:10.1029/2011GL050207.

This paper presents a model study of the impact of drifting snow on the lower atmosphere, surface snow characteristics, and surface mass balance of Antarctica. We use the regional atmospheric climate model RACMO2.1/ANT with a high horizontal resolution (27 km), equipped with a drifting snow routine and forced by ERA-Interim (1989–2009) at its lateral and ocean boundaries. Drifting snow sublimation (SUds) is significant in Antarctica, especially in the coastal regions (>150 mm water equivalent yr!1). Integrated over the ice sheet, SUds removes "6% of the annually precipitated snow. Drifting snow interacts with the atmosphere, as it increases the lower atmospheric moisture content and reduces surface sublimation (SUs), and leads to increased snowfall in regions where the atmosphere usually is close to saturation. Drifting snow sublimation (SUds) is smallest in summer, when katabatic wind speeds are lower and melting and surface sublimation consolidate the snow surface. Compared to a simulation without drifting snow, total sublimation (SUds + SUs) doubles on the grounded ice sheet if drifting snow is considered. Drifting snow erosion is locally significant, but can be neglected on a continent-wide scale.

[1] Freshwater (FW) fluxes from river runoff and precipitation minus evaporation for the pan Arctic seas are relatively well documented and prescribed in ocean GCMs. Fluxes from Greenland on the other hand are generally ignored altogether, despite their potential impacts on ocean circulation and marine biology. Here, we present a reconstruction of the spatially distributed FW flux from Greenland for 1958–2010. We find a modest increase into the Arctic Ocean during this period. Fluxes into the Irminger Basin, however, have increased by fifty percent (6.3 ± 0.5 km3 yr−2) in less than twenty years. This greatly exceeds previous estimates. For the ice sheet as a whole the rate of increase since 1992 is 16.9 ± 1.8 km3 yr−2. The cumulative FW anomaly since 1995 is 3200 ± 358 km3, which is about a third of the magnitude of the Great Salinity Anomaly (GSA) of the 1970s. If this trend continues into the future, the anomaly will exceed that of the GSA by about 2025.

Several Antarctic Peninsula (AP) ice shelves have lost significant fractions of their volume over the past decades, coincident with rapid regional climate change. Wilkins Ice Shelf (WIS), on the western side of the AP, is the most recent, experiencing a sequence of large calving events in 2008 and 2009. We analyze the mass balance for WIS for the period 1992 2008 and find that the averaged rate of ice-shelf thinning was 0.8 m a 1, driven by a mean basal melt rate of 〈wb〉 = 1.3 0.4 m a 1. Interannual variability was large, associated with changes in both surface mass accumulation and 〈wb〉. Basal melt rate declined significantly around 2000 from 1.8 0.4 m a 1 for 1992–2000 to 0.75 0.55 m a 1 for 2001–2008; the latter value corresponding to approximately steady-state ice-shelf mass. Observations of ocean temperature T obtained during 2007–2009 by instrumented seals reveal a cold, deep halo of Winter Water (WW; T ≈ 1.6°C) surrounding WIS. The base of the WW in the halo is 170 m, approximately the mean ice draft for WIS. We hypothesize that the transition in 〈wb〉 in 2000 was caused by a small perturbation ( 10–20 m) in the relative depths of the ice base and the bottom of the WW layer in the halo. We conclude that basal melting of thin ice shelves like WIS is very sensitive to upper-ocean and coastal processes that act on shorter time and space scales than those affecting basal melting of thicker West Antarctic ice shelves such as George VI and Pine Island Glacier.

A new, high resolution (27 km) surface mass balance (SMB) map of the Antarctic ice sheet is presented, based on output of a regional atmospheric climate model that includes snowdrift physics and is forced by the most recent reanalysis data from the European Centre for Medium-Range Weather Forecasts (ECMWF), ERA-Interim (1979–2010). The SMB map confirms high accumulation zones in the western Antarctic Peninsula (>1500 mm y!1) and coastal West Antarctica (>1000 mm y!1), and shows low SMB values in large parts of the interior ice sheet (

A new, high resolution (27 km) surface mass balance (SMB) map of the Antarctic ice sheet is presented, based on output of a regional atmospheric climate model that includes snowdrift physics and is forced by the most recent reanalysis data from the European Centre for Medium-Range Weather Forecasts (ECMWF), ERA-Interim (1979–2010). The SMB map confirms high accumulation zones in the western Antarctic Peninsula (>1500 mm y!1) and coastal West Antarctica (>1000 mm y!1), and shows low SMB values in large parts of the interior ice sheet (

We present a sensitivity study of the surface mass balance (SMB) of the Greenland Ice Sheet, as modeled using a regional atmospheric climate model, to various parameter settings in the albedo scheme. The snow albedo scheme uses grain size as a prognostic variable and further depends on cloud cover, solar zenith angle and black carbon concentration. For the control experiment the overestimation of absorbed shortwave radiation (+6%) at the K-transect (west Greenland) for the period 2004–2009 is considerably reduced compared to the previous density-dependent albedo scheme (+22%). To simulate realistic snow albedo values, a small concentration of black carbon is needed, which has strongest impact on melt in the accumulation area. A background ice albedo field derived from MODIS imagery improves the agreement between the modeled and observed SMB gradient along the K-transect. The effect of enhanced meltwater retention and refreezing is a decrease of the albedo due to an increase in snow grain size. As a secondary effect of refreezing the snowpack is heated, enhancing melt and further lowering the albedo. Especially in a warmer climate this process is important, since it reduces the refreezing potential of the firn layer that covers the Greenland Ice Sheet.

This paper presents the impact of model resolution on the simulated wind speed, drifting snow climate and surface mass balance (SMB) of Terre Adélie and its surroundings, East Antarctica. We compare regional climate model simulations at 27 and 5.5 km resolution for the year 2009. The wind speed maxima in Terre Adélie and the narrow glacial valleys of Victoria Land are better represented at 5.5 km resolution, because the topography is better resolved. Drifting snow sublimation is >100 mm a-1 in regions with high wind speeds. Our results indicate a strong feedback between topography, wind gradients and drifting snow erosion. As a result, SMB shows much more local spatial variability at 5.5 km resolution that is controlled by drifting snow erosion, whereas the large-scale SMB gradient is largely determined by precipitation. Drifting snow processes lead to ablation in the narrow glacial valleys of Victoria Land. The integrated SMB equals 86 Gt. Although wind climate, drifting snow processes and SMB variability are better represented at 5.5 km, the area-integrated SMB is not significantly different between the simulations at 27 and 5.5 km. A horizontal resolution of 27 km is sufficient to realistically simulate ice-sheet wide SMB.

[1] A parameterization for broadband snow surface albedo, based on snow grain size evolution, cloud optical thickness, and solar zenith angle, is implemented into a regional climate model for Antarctica and validated against field observations of albedo for the period 1995–2004. Over the Antarctic continent, modeled snow grain size exhibits expected behavior. The agreement between modeled and observed albedo at Neumayer, Dronning Maud Land, is very good, and subtle variability in albedo is well captured by the model. December–February mean differences in modeled and observed net shortwave radiation range from −8.7 to +3.8 W m−2 between 1995 and 2004, with a mean value of −2.7 W m−2. This is a considerable improvement compared to the previous albedo parameterization in the model, which led to overestimates of the net shortwave fluxes by +15.0 to +22.7 W m−2, or 40–55% of the observed net shortwave flux, in the same period.

Until quite recently, the mass balance (MB) of the great ice sheets of Greenland and Antarctica was poorly known and often treated as a residual in the budget of oceanic mass and sea level change. Recent developments in regional climate modelling and remote sensing, especially altimetry, gravimetry and InSAR feature tracking, have enabled us to specifically resolve the ice sheet mass balance components at a near-annual timescale. The results reveal significant mass losses for both ice sheets, caused by the acceleration of marine-terminating glaciers in southeast, west and northwest Greenland and coastal West Antarctica, and increased run-off in Greenland. At the same time, the data show that interannual variability is very significant, masking the underlying trends.

Here, we present a 20‐year record of monthly ice sheet mass balance for Greenland and Antarctica. We examine and reconcile two independent methods for estimating temporal variations in ice sheet mass balance, the mass budget method (MB) and the gravity method, during the last 8 years. The MBM compares the surface mass balance (SMB; i.e., the sum of snowfall minus surface ablation) reconstructed from regional atmospheric models with perimeter loss (D; ice discharge) calculated from a time series of glacier velocity and ice thickness to deduce the rate of mass change, dM/dt [Rignot and Kanagaratnam, 2006; Howat et al., 2007; Rignot et al., 2008a; van den Broeke et al., 2009]. The gravity method employs a monthly time series of time‐variable gravity data from the Gravity Recovery and Climate Experiment (GRACE) to estimate the relative mass as a function of time, M(t) [e.g., Velicogna and Wahr, 2006]. We resolve the differences between the two methods in terms of mass balance, dM(t)/dt, and acceleration of mass loss, d2M/dt2, and conclude by discussing the contribution of the ice sheets to sea level in recent and forthcoming decades.

Acceleration of Greenland's three largest outlet glaciers, Helheim, Kangerdlugssuaq and Jakobshavn Isbræ, accounted for a substantial portion of the ice sheet's mass loss over the past decade. Rapid changes in their discharge, however, make their cumulative mass-change uncertain. We derive monthly mass balance rates and cumulative balance from discharge and surface mass balance (SMB) rates for these glaciers from 2000 through 2010. Despite the dramatic changes observed at Helheim, the glacier gained mass over the period, due primarily to the short-duration of acceleration and a likely longer-term positive balance. In contrast, Jakobshavn Isbræ lost an equivalent of over 11 times the average annual SMB and loss continues to accelerate. Kangerdlugssuaq lost over 7 times its annual average SMB, but loss has returned to the 2000 rate. These differences point to contrasts in the long-term evolution of these glaciers and the danger in basing predictions on extrapolations of recent changes.

In this paper, we estimate the contribution of snowdrift sublimation (SUds) to the surface mass balance at Neumayer, located on the Ekström ice shelf in Eastern Antarctica. A single column version of the RACMO2-ANT model is used as a physical interpolation tool of high-quality radiosonde and surface measurements for a 15-yr period (1993–2007), and combined with a routine to calculate snowdrift sublimation and horizontal snow transport. The site is characterised by a relatively mild, wet and windy climate, so snowdrift is a common phenomenon. The modelled timing and frequency of snowdrift events compares well with observations. This is further illustrated by an additional simulation for Kohnen base, where the timing of snowdrift is realistic, although the modelled horizontal transport is overestimated. Snowdrift sublimation is mainly dependent on wind speed, but also on relative humidity and temperature. During high wind speeds, SUds saturates and cools the air, limiting its own strength. We estimate that SUds removes around 16%±8% of the accumulated snow from the surface. The total sublimation more than triples when snowdrift is considered, although snowdrift sublimation limits sublimation at the surface. SUds shows a strong seasonal cycle, as well as large inter-annual variability. This variability can be related to the variability of the atmospheric conditions in the surface layer.