Precipitation and snow cover in the Himalaya: from reanalysis to regional climate simulations

We applied a Regional Climate Model (RCM) to simulate precipitation and snow cover over the Himalaya, between March 2000 and December 2002. Due to its higher resolution, our model simulates a more realistic spatial variability of wind and precipitation than those of the reanalysis of the European Centre of Medium range Weather Forecast (ECMWF) used as lateral boundaries. In this region, we found very large discrepancies between the estimations of precipitation provided by reanalysis, rain gauges networks, satellite observations, and our RCM simulation. Our model clearly underestimates precipitation at the foothills of the Himalaya and in its eastern part. However, our simulation provides a first estimation of liquid and solid precipitation in high altitude areas, where satellite and rain gauge networks are not very reliable. During the two years of simulation, our model resembles the snow cover extent and duration quite accurately in these areas. Both snow accumulation and snow cover duration differ widely along the Himalaya: snowfall can occur during the whole year in western Himalaya, due to both summer monsoon and mid-latitude low pressure systems bringing moisture into this region. In Central Himalaya and on the Tibetan Plateau, a much more marked dry season occurs from October to March. Snow cover does not have a pronounced seasonal cycle in these regions, since it depends both on the quite variable duration of the monsoon and on the rare but possible occurrence of snowfall during the extra-monsoon period

Modelling snow accumulation on Greenland in Eemian, glacial inception, and modern climates in a GCM

Changing climate conditions on Greenland influence the snow accumulation rate and surface mass balance (SMB) on the ice sheet and, ultimately, its shape. This can in turn affect local climate via orography and albedo variations and, potentially, remote areas via changes in ocean circulation triggered by melt water or calving from the ice sheet. Examining these interactions in the IPSL global model requires improving the representation of snow at the ice sheet surface. In this paper, we present a new snow scheme implemented in LMDZ, the atmospheric component of the IPSL coupled model. We analyse surface climate and SMB on the Greenland ice sheet under insolation and oceanic boundary conditions for modern, but also for two different past climates, the last glacial inception (115 kyr BP) and the Eemian (126 kyr BP). While being limited by the low resolution of the general circulation model (GCM), present-day SMB is on the same order of magnitude as recent regional model findings. It is affected by a moist bias of the GCM in Western Greenland and a dry bias in the north-east. Under Eemian conditions, the SMB decreases largely, and melting affects areas in which the ice sheet surface is today at high altitude, including recent ice core drilling sites as NEEM. In contrast, glacial inception conditions lead to a higher mass balance overall due to the reduced melting in the colder summer climate. Compared to the widely applied positive degree-day (PDD) parameterization of SMB, our direct modelling results suggest a weaker sensitivity of SMB to changing climatic forcing. For the Eemian climate, our model simulations using interannually varying monthly mean forcings for the ocean surface temperature and sea ice cover lead to significantly higher SMB in southern Greenland compared to simulations forced with climatological monthly means

Evaluation of a high-resolution regional climate simulation over Greenland

A simulation of the 1991 summer has been performed over south Greenland with a coupled atmosphere–snow regional climate model (RCM) forced by the ECMWF re-analysis. The simulation is evaluated with in-situ coastal and ice-sheet atmospheric and glaciological observations. Modelled air temperature, specific humidity, wind speed and radiative fluxes are in good agreement with the available observations, although uncertainties in the radiative transfer scheme need further investigation to improve the model’s performance. In the sub-surface snow-ice model, surface albedo is calculated from the simulated snow grain shape and size, snow depth, meltwater accumulation, cloudiness and ice albedo. The use of snow metamorphism processes allows a realistic modelling of the temporal variations in the surface albedo during both melting periods and accumulation events. Concerning the surface albedo, the main finding is that an accurate albedo simulation during the melting season strongly depends on a proper initialization of the surface conditions which mainly result from winter accumulation processes. Furthermore, in a sensitivity experiment with a constant 0.8 albedo over the whole ice sheet, the average amount of melt decreased by more than 60%, which highlights the importance of a correctly simulated surface albedo. The use of this coupled atmosphere–snow RCM offers new perspectives in the study of the Greenland surface mass balance due to the represented feedback between the surface climate and the surface albedo, which is the most sensitive parameter in energy-balance-based ablation calculations.Peer reviewe

Dynamic recycling of gaseous elemental mercury in the boundary layer of the Antarctic Plateau

International audienceGaseous elemental mercury (Hg0) was investigated in the troposphere and in the interstitial air extracted from the snow at Dome Concordia station (alt. 3320 m) on the Antarctic Plateau during January 2009. Measurements and modeling studies showed evidence of a very dynamic and daily cycling of Hg0 inside the mixing layer with a range of values from 0.2 ng m−3 up to 2.3 ng m−3. During low solar irradiation periods, fast Hg0 oxidation processes in a confined layer were suspected. Unexpectedly high Hg0 concentrations for such a remote place were measured under higher solar irradiation due to snow photochemistry. We suggest that a daily cycling of reemission/oxidation occurs during summer within the mixing layer at Dome Concordia. Hg0 concentrations showed a negative correlation with ozone mixing ratios, which contrasts with atmospheric mer- cury depletion events observed during the Arctic spring. Un- like previous Antarctic studies, we think that atmospheric Hg0 removal may not be the result of advection processes. The daily and dramatic Hg0 losses could be a consequence of surface or snow induced oxidation pathways. It remains however unclear whether halogens are involved. The cycling of other oxidants should be investigated together with Hg species in order to clarify the complex reactivity on the Antarctic plateau

The 1979-2005 Greenland ice sheet melt extent from passive microwave data using an improved version of the melt retrieval XPGR algorithm

Analysis of passive microwave satellite observations over the Greenland ice sheet reveals a significant increase in surface melt over the period 1979-2005. Since 1979, the total melt area was found to have increased +1.22 x 10ˆ7 kmˆ2. An improved version of the cross-polarized gradient ratio (XPGR) technique is used to identify the melt from the brightness temperatures. The improvements in the melt retrieval XPGR algorithm as well as the surface melt acceleration are discussed with results from a coupled atmosphere-snow regional climate model. From 1979 to 2005, the ablation period increases everywhere over the melt zone except in the regions where the model simulates an increased summer snowfall. Indeed, more snowfall in summer decreases the liquid water content of the snowpack, raises the albedo and therefore reduces the melt. Finally, this melt acceleration over the Greenland ice sheet is highly correlated with both Greenland and global warming suggesting a continuing surface melt increase in the future.Peer reviewe

Sensitivity of a Greenland ice sheet model to atmospheric forcing fields

International audiencePredicting the climate for the future and how it will impact ice sheet evolution requires coupling ice sheet models with climate models. However, before we attempt to develop a realistic coupled setup, we propose, in this study, to first analyse the impact of a model simulated climate on an ice sheet. We undertake this exercise for a set of regional and global climate models. Modelled near surface air temperature and precipitation are provided as upper boundary conditions to the GRISLI (GRenoble Ice Shelf and Land Ice model) hybrid ice sheet model (ISM) in its Greenland configuration. After 20 kyrs of simulation, the resulting ice sheets highlight the differences between the climate models. While modelled ice sheet sizes are generally comparable to the observed one, there are considerable deviations among the ice sheets on regional scales. These deviations can be explained by biases in temperature and precipitation near the coast. This is especially true in the case of global models. But the deviations between the climate models are also due to the differences in the atmospheric general circulation. To account for these differences in the context of coupling ice sheet models with climate models, we conclude that appropriate downscaling methods will be needed. In some cases, systematic corrections of the climatic variables at the interface may be required to obtain realistic results for the Greenland ice sheet (GIS)

Atmospheric drying as the main driver of dramatic glacier wastage in the southern Indian Ocean

This study was funded by IPEV-1048 GLACIOCLIM-KESAACO and LEFE-INSU KCRuMBLE programs, and the Agence Nationale de la Recherche through contract ANR-14-CE01-0001-01 (ASUMA).The ongoing retreat of glaciers at southern sub-polar latitudes is particularly rapid and widespread. Akin to northern sub-polar latitudes, this retreat is generally assumed to be linked to warming. However, no long-term and well-constrained glacier modeling has ever been performed to confirm this hypothesis. Here, we model the Cook Ice Cap mass balance on the Kerguelen Islands (Southern Indian Ocean, 49°S) since the 1850s. We show that glacier wastage during the 2000s in the Kerguelen was among the most dramatic on Earth. We attribute 77% of the increasingly negative mass balance since the 1960s to atmospheric drying associated with a poleward shift of the mid-latitude storm track. Because precipitation modeling is very challenging for the current generation of climate models over the study area, models incorrectly simulate the climate drivers behind the recent glacier wastage in the Kerguelen. This suggests that future glacier wastage projections should be considered cautiously where changes in atmospheric circulation are expected.Publisher PDFPeer reviewe

"Important role of the mid-tropospheric atmospheric circulation in the recent surface melt increase over the Greenland ice sheet"

Since 2007, there has been a series of surface melt records over the Greenland ice sheet (GrIS), continuing the trend towards increased melt observed since the end of the 1990s. The last two decades are characterized by an increase of negative phases of the North Atlantic Oscillation (NAO) favouring warmer and drier summers than normal over GrIS. In this context, we use a circulation type classification based on daily 500 hPa geopotential height to evaluate the role of atmospheric dynamics in this surface melt acceleration for the last two decades. Due to the lack of direct observations, the interannual melt variability is gauged here by the summer (June-July-August) mean temperature from reanalyses at 700 hPa over Greenland; analogous atmospheric circulations in the past show that similar to 70% of the 1993-2012 warming at 700 hPa over Greenland has been driven by changes in the atmospheric flow frequencies. Indeed, the occurrence of anticyclones centred over the GrIS at the surface and at 500 hPa has doubled since the end of 1990s, which induces more frequent southerly warm air advection along the western Greenland coast and over the neighbouring Canadian Arctic Archipelago (CAA). These changes in the NAO modes explain also why no significant warming has been observed these last summers over Svalbard, where northerly atmospheric flows are twice as frequent as before. Therefore, the recent warmer summers over GrIS and CAA cannot be considered as a long-term climate warming but are more a consequence of NAO variability affecting atmospheric heat transport. Although no global model from the CMIP5 database projects subsequent significant changes in NAO through this century, we cannot exclude the possibility that the observed NAO changes are due to global warming

Reconstructions of the 1900–2015 Greenland ice sheet surface mass balance using the regional climate MAR model

With the aim of studying the recent Greenland ice sheet (GrIS) surface mass balance (SMB) decrease relative to the last century, we have forced the regional climate MAR (Modèle Atmosphérique Régional; version 3.5.2) model with the ERA-Interim (ECMWF Interim Re-Analysis; 1979–2015), ERA-40 (1958–2001), NCEP–NCARv1 (National Centers for Environmental Prediction–National Center for Atmospheric Research Reanalysis version 1; 1948–2015), NCEP–NCARv2 (1979–2015), JRA-55 (Japanese 55-year Reanalysis; 1958–2014), 20CRv2(c) (Twentieth Century Reanalysis version 2; 1900–2014) and ERA-20C (1900–2010) reanalyses. While all these forcing products are reanalyses that are assumed to represent the same climate, they produce significant differences in the MAR-simulated SMB over their common period. A temperature adjustment of +1 °C (respectively −1 °C) was, for example, needed at the MAR boundaries with ERA-20C (20CRv2) reanalysis, given that ERA-20C (20CRv2) is ∼ 1 °C colder (warmer) than ERA-Interim over Greenland during the period 1980–2010. Comparisons with daily PROMICE (Programme for Monitoring of the Greenland Ice Sheet) near-surface observations support these adjustments. Comparisons with SMB measurements, ice cores and satellite-derived melt extent reveal the most accurate forcing datasets for the simulation of the GrIS SMB to be ERA-Interim and NCEP–NCARv1. However, some biases remain in MAR, suggesting that some improvements are still needed in its cloudiness and radiative schemes as well as in the representation of the bare ice albedo. Results from all MAR simulations indicate that (i) the period 1961–1990, commonly chosen as a stable reference period for Greenland SMB and ice dynamics, is actually a period of anomalously positive SMB (∼ +40 Gt yr−1) compared to 1900–2010; (ii) SMB has decreased significantly after this reference period due to increasing and unprecedented melt reaching the highest rates in the 120-year common period; (iii) before 1960, both ERA-20C and 20CRv2-forced MAR simulations suggest a significant precipitation increase over 1900–1950, but this increase could be the result of an artefact in the reanalyses that are not well-enough constrained by observations during this period and (iv) since the 1980s, snowfall is quite stable after having reached a maximum in the 1970s. These MAR-based SMB and accumulation reconstructions are, however, quite similar to those from Box (2013) after 1930 and confirm that SMB was quite stable from the 1940s to the 1990s. Finally, only the ERA-20C-forced simulation suggests that SMB during the 1920–1930 warm period over Greenland was comparable to the SMB of the 2000s, due to both higher melt and lower precipitation than normal

Twentieth century temperature and snow cover changes in the French Alps

International audienceChanges in snow cover associated with the warming of the French Alps greatly influence social-ecological systems through their impact on water resources, mountain ecosystems, economic activities, and glacier mass balance. In this study, we investigated trends in snow cover and temperature over the twentieth century using climate model and reanalysis data. The evolution of temperature, precipitation and snow cover in the European Alps has been simulated with the Modèle Atmospherique Régional (MAR) applied with a 7-km horizontal resolution and driven by ERA-20C (1902-2010) and ERA5 (1981–2018) reanalyses data. Snow cover duration and snow water equivalent (SWE) simulated with MAR are compared to the SAFRAN - SURFEX-ISBA-Crocus - MEPRA meteorological and snow cover reanalysis (S2M) data across the French Alps (1958–2018) and in situ glacier mass balance measurements. MAR outputs provide a realistic distribution of SWE and snow cover duration as a function of elevation in the French Alps. Large disagreements are found between the datasets in terms of absolute warming trends over the second part of the twentieth century. MAR and S2M trends are in relatively good agreement for the decrease in snow cover duration, with higher decreases at low elevation ( $\sim$ ∼ 5–10%/decade). Consistent with other studies, the highest warming rates in MAR occur at low elevations ( 2000 m a.s.l) in summer. In spring, warming trends show a maximum at intermediate elevations (1500 to 1800 m). Our results suggest that higher warming at these elevations is mostly linked to the snow-albedo feedback in spring and summer caused by the disappearance of snow cover at higher elevation during these seasons. This work has evidenced that depending on the season and the period considered, enhanced warming at higher elevations may or may not be found. Additional analysis in a physically comprehensive way and more high-quality dataset, especially at high elevations, are still required to better constrain and quantify climate change impacts in the Alps and its relation to elevation