Boreal snow cover variations induced by aerosol emissions in the middle of the 21 st century

Introduction Conclusions References

by 2050, with a decrease ranging from 10 to 100 days from pre-industrial values over large parts of the Northern Hemisphere.This reduction is essentially due to temperature increase, which is quite strong in the RCP8.5 scenario in the absence of climate mitigation policies.Moreover, the projected sea-ice retreat in the next decades will open new routes for shipping in the Arctic.However, a large increase in shipping emissions in the Arctic by the mid 21st century does not lead to significant changes of BC deposition over snow-covered areas in our simulation.Therefore, the MNDWS is clearly not affected through snow darkening effects associated to these Arctic ship emissions.In an experiment without nudging toward atmospheric reanalyses, we simulated however some changes of the MNDWS considering such aerosol ship emissions.These changes are generally not statistically significant in boreal continents, except in the Quebec and in the West Siberian plains, where they range between −5 and −10 days.They are induced both by radiative forcings of the aerosols when they are in the 1 Introduction The boreal regions have been characterized as a region very sensitive to climate change (Lemke et al., IPCC, 2007, chapter 4).One reason for the amplification in Arctic and Subarctic surface warming in response to increased greenhouse gas concentrations is the snow and sea-ice albedo feedback, which decreases surface albedo as snow and sea ice further melt and disappear in response to the warming by greenhouse gases (Serreze et al., 2006;Qu et al., 2007).Both sea-ice and snow-cover extents have been observed to shrink over the last decades in the Northern Hemisphere (Serreze et al., 2007;Shi et al., 2011).Snow-cover extent is expected to decrease further during the 21st century (e.g.Hosaka et al., 2005;Frei and Gong, 2005).However, it is quite difficult to evaluate accurately this decrease using climate models, because of both the complexity of the interactions between the snow and the atmosphere and the uncertainties when predicting future anthropogenic climate forcing (Qu andHall, 2006 and2007;Ghatak et al., 2010).
In contrast with the Antarctic, the Arctic atmosphere is quite polluted.An ensemble of short-lived species emitted in the industrialised mid-latitude regions of the Northern Hemisphere are transported towards the Arctic, where their lifetime increases due to the weak intensity of removal processes, in particular during the winter.The transport of pollutants into the Arctic atmosphere occurs especially in spring, and has been referred to cause the "Arctic Haze" phenomenon (e.g.Shaw, 1995;Stohl et al., 2006).Ozone and aerosols are the main short-lived species transported toward the Arctic that impact significantly the climate of this region, modifying regionally the radiative balance of the atmosphere (Law and Stohl, 2007).Ozone is a strong greenhouse gas, inducing a positive radiative forcing and causing a regional increase of the surface temperature (Shindell et al., 2006).Sulphate, Organic Carbon (OC) and nitrate aerosols are known to scatter solar radiation, inducing a negative radiative forcing at the top of the atmosphere and a cooling of the Earth's surface (Penner et al., 2001;Kanakidou, 2005).Black Carbon (BC) strongly absorbs solar radiation, inducing a positive forcing at the top of the atmosphere and a negative instantaneous forcing at the surface (Reddy et al., 2005).The heating of the atmosphere due to BC induces also an increase in the downward longwave radiation.Over highly reflective surfaces like snow covered areas, this increase in the longwave flux can be higher than the decrease of the shortwave flux induced by atmospheric BC (Quinn et al., 2008).In addition to these direct radiative forcings, aerosols affect clouds microphysics, processes referred to as the aerosol indirect effects.Although uncertain, these effects are thought to induce a negative radiative forcing, both at the top and the bottom of the troposphere (Lohmann et al., 2005).However, it has been suggested that there is also a longwave positive radiative forcing from aerosol-cloud interactions in the Arctic (Garrett and Zhao, 2006;Lubin and Vogelmann, 2006).In addition, once deposited to snow or ice, BC and OC absorb radiation within the snowpack, and cause an earlier snow disappearance or decrease the snow mass, inducing a positive forcing at the surface, through decreased albedo (e.g. Warren and Wiscombe, 1980;Clarke and Noone, 1985;Jacobson, 2004;Hadley and Kirchstetter, 2012).Overall, Shindell and Faluvegi (2009) pointed out that the temperature response to a radiative forcing is not necessarily correlated with the location of this radiative forcing.This is particularly true for the Arctic surface temperature response, which can be of opposite sign to that of the radiative forcing.This points to the necessity to apply Global Circulation Models (GCM) to quantify the surface temperature Introduction

Conclusions References
Tables Figures

Back Close
Full response to different radiative forcings in a particular region.Overall, Shindell (2007) and Shindell and Faluvegi (2009) estimate that both anthropogenic well-mixed greenhouses gases and short-lived species have contributed to the Arctic warming.
The main source of aerosol in the Arctic atmosphere is the transport from polluted regions in North America, Europe and Asia, while local aerosol emissions are very small (Shindell et al., 2008;Browse et al., 2012).Future aerosol concentrations in the Arctic are therefore very dependent on the evolution of the anthropogenic emissions from these regions.According to the Radiative Concentration Pathway (RCP, Moss et al., 2008) emission scenarios, aerosol emissions in Northern America and Europe are estimated to have reached maximum values at different time periods during the 20th century, depending on countries and on the chemical species under consideration (Bond et al., 2007;Smith et al., 2011).These regions now experience a significant decrease in their aerosol emissions.This is not the case of Asian emissions, which are still increasing.Their decrease is projected to take place in the next decades, although the exact timing is quite difficult to estimate, as the projections for energy demand, biofuel consumption and the introduction of new technologies are not set in stone (Ohara et al., 2007).In addition to anthropogenic emissions occurring in densely populated and industrialized regions, it seems that two local sources could affect the Arctic atmosphere in the decades to come: first, ship emissions could increase significantly, as summer sea-ice retreat will open new routes across the Arctic ocean (Corbett et al., 2011); second, biomass burning emissions are expected to become stronger and to occur earlier in the season.This phase advanced is observed in the higher latitudes during warmer and dryer spring periods, in response to climate warming (e.g.Warneke et al., 2009).Flannigan et al. (2009a, b) projected for instance that climate warming will induce an increase of fire activity in temperate and boreal regions, mainly from forest wildfires.
The goal of this study is to estimate snow-cover variations for the middle of the 21st century in the continental pan-Arctic area using simulations with a global coupled atmospheric general circulation and chemistry model prescribed with different aerosol Introduction

Conclusions References
Tables Figures

Back Close
Full local emission scenarios in the boreal region.In addition, using a land surface model enhanced for including the effects of BC on snow albedo, we investigate how the deposition of absorbing aerosols on snow affects snow cover dynamics and feedbacks on regional climate.We evaluate the snow-cover changes for the intensive RCP8.5 scenario (Representative Concentration Pathway 8.5, Moss et al., 2008 and2010;Riahi, 2007) in the 2050 decade, and analyse thereafter the role of possible enhanced aerosols local emissions in the Arctic region.
2 Experimental set-up

Model description
We used the "LMDZ-INCA-ORCHIDEE" atmospheric General Circulation Model to study the interactions between atmosphere, aerosols and snow-covered areas.This model consists of three coupled modules: the LMDZ general circulation model represents the atmospheric component (Hourdin et al., 2006).INCA (Interactions between Chemistry and Aerosols) describes gas-and aqueous-phase chemistry (Hauglustaine et al., 2004;Boucher et al., 2002), as well as aerosol physical properties such as size and hygroscopicity (Balkanski et al., 2010), which control the amount of wet and dry deposition.The coupling between the LMDZ and INCA models allow for an interactive simulation of five aerosol chemical species, namely sulphate, BC, OC, sea-salt and dust.Direct aerosol forcing is taken into account for BC, OC, seasalt and dust, and both direct and indirect effect are taken into account for sulphate, BC and OC aerosol, as described in Deandreis et al. (2012).We used here LMDZ and INCA with a horizontal resolution of 96 × 95 grid points in longitude and latitude, and with a vertical discretisation of 19 layers.Finally the ORCHIDEE land surface model serves as the land surface boundary condition for LMDZ and describes exchanges of energy and water between the atmosphere, the soil and the biosphere (Krinner et al., 2005), including a dynamic snow module.Introduction

Conclusions References
Tables Figures

Back Close
Full For this work, we used the detailed representation of snow-cover implemented in ORCHIDEE by Krinner et al. (2006) who studied the interactions between dust aerosol and ice-sheets in Northern Asia during the last glacial maximum (21 000 yr BP).In this scheme, snow albedo and snow cover are described separately for forests and grasslands/deserts, with a subgrid-scale orographic variability to compute accurately the energy balance in mountainous areas (Douville et al., 1995;Roesch et al., 2001).The aerosol content of the snow and its albedo are computed with a two-layer scheme, with a top layer of 8 mm (Snow Water Equivalent, SWE), and a bottom layer containing the remaining snow.A detailed description of the treatment of the snow/aerosol interactions in ORCHIDEE can be found in Krinner et al. (2006).However, only dry-deposited dust aerosol was taken into account in this study.Here, we also take into account BC, as its very absorbing property makes it likely to impact significantly the snowpack energy balance and the snow cover extent (e.g.Jacobson, 2004).Unfortunately, OC deposition on snow is not taken into account in our simulation.This aerosol also absorbs solar radiation, but there remain a lot of uncertainties concerning its radiative properties and its behaviour within the snowpack.We hope to take these processes into account in a further study.BC dry and wet depositions are computed by the INCA atmospheric chemistry module with a six-hourly time step.As in Krinner et al. (2006), dry deposition contributes to increase the aerosol content in the top snow layer.Wet deposition also supplies aerosol to the surface layer, but it should be noted that this process is associated with an entry of fresh snow.If snowfall brings more snow than the maximum height of the snowpack surface layer, then aerosols in this previous surface layer are transferred into the bottom layer (note that we considered a constant snow density of 330 kg m −3 ).Conway et al. (1996) observed that BC could be flushed effectively trough the snow in melting conditions, with velocities strongly dependent on the particle size.snowpack even during melting conditions.Building on this experimental evidence, and in contrast with Krinner et al. (2006), we will consider in this study that both dust and BC do not flush through the snow, and stay at the surface until a new snowfall occurs or until the disappearance of the snow-cover.This assumption could overestimate the magnitude of BC aerosol effects on the snow cover and climate.
Snow albedo is estimated using the parameterisation of Warren and Wiscombe (1980), which is adapted for snow containing aerosols.As in Krinner et al. (2006), the snow albedo of the bottom snowpack layer is computed first for diffuse radiation as a function of the underlying albedo, snow grain size and aerosol content.The spherical albedo of the bottom layer is then used as the underlying albedo for computing the albedo of the surface layer, both for diffuse and direct solar radiation.Snow albedo is averaged separately in the visible and near-infrared parts of the solar spectrum.We adopt the same aerosol physical properties as used in Balkanski et al. (2010) to evaluate their radiative forcings in the atmosphere.Within the snow, we do not know the extent to which aerosols are internally mixed, how they interact with snow grains, and how their hygroscopic and radiative properties evolve in time.Faced with all these uncertainties, we decided to consider simpler physical and radiative properties for aerosols in the snow in comparison with atmospheric aerosols.The size and radiative parameters for dust are the same as used by Krinner et al. (2006), following Guelle et al. (2000) and Balkanski et al. (2007).Black carbon is assumed to follow a log-normal size distribution with a median number diameter of 11.8 nm, characteristic of freshly emitted soot (Dentener et al., 2006;Jacobson et al., 2004).In the real world, this diameter increases quickly, as BC undergoes ageing and coagulation and can be coated by other aerosols in the atmosphere.However, as we do not consider internal mixtures for BC in snow, we consider that BC aerosols regain their initial size when incorporated in the snowpack.We considered a BC density of 1 g cm −3 , and the refractive index for BC is taken to be m = 1.75-0.45i .Refractive indices for ice are taken from the GEISA database (Jacquinet-Husson et al., 1999).Introduction

Conclusions References
Tables Figures

Back Close
Full

Description of simulations
Table 1 describes the eight 11-yr global simulations that we realized to characterize the impact of BC deposition on snow cover both for the present period and for the middle of the 21st century.We exclude from our analysis the first year of simulation, considered as a spin-up period.The two first experiments -designated as S1 and S1B -describe the present-day atmospheric state (1998-2008), using prescribed observed Sea Surface Temperature (SST, see Rayner et al., 2003) with winds nudged toward ERA-40 reanalysis from the European Centre for Medium-range Weather Forecasts (ECMWF).Note that pressure, temperature and humidity are computed with the LMDZ model without nudging in these experiments.The nudging is applied only for horizontal winds as described in Coindreau et al. (2006).These experiments were conducted with the present-day global aerosol emission inventory described in Lamarque et al. (2010), an inventory made for the Coupled Model Inter-comparison Project Phase 5 (CMIP5, CLI-VAR special issue, 2011).In S1B, the BC content in the snow is set to zero, whereas it is computed from aerosol deposition in all the other experiments.The six other experiments were conducted over the period 2050-2060.They are based upon the aerosol and gases intensive emission scenario RCP8.5 (Representative Concentration Pathway 8.5, Moss et al., 2008 and2010;Riahi, 2007), characteristic of a scenario with no climate mitigation policies to limit greenhouse gas emissions.This scenario corresponds to a total anthropogenic forcing in 2100 of approximately 8.5 W m −2 .All six experiments were conducted with prescribed SST for the 2050s decade as produced from a previous coupled ocean-atmosphere simulation using IPSL-CM5A configuration in the context of the CMIP5 exercise (Dufresne et al., 2012).The first one of these six experiments -designated as S2 -has been performed with the aerosol emission inventory corresponding to that defined for the RCP8.5 scenario (Lamarque et al., 2009).
Importantly, none of the RCP emission inventories used in CMIP5 simulations over the 21st century considers variations of "local" emissions in the Arctic, which could be associated to a significant increase in ship traffic in the Arctic or to an intensification of Introduction

Conclusions References
Tables Figures

Back Close
Full biomass burning in boreal and temperate regions.For this reason, we performed another simulation -S3 -similar to S2 but replacing the baseline Arctic ship emissions in the RCP8.5 2050 by a scenario that includes important ship traffic over Arctic routes.These higher ship emissions are based on the "high-growth" scenario of Corbett et al. (2010), where they consider a high increase in ship traffic over the Arctic, including current routes conjugated with diversion routes opened during summer following the seasonal retreat of sea-ice expected in the next decades.Finally an S4 simulation was also performed, similar to S2, but with enhanced biomass burning activity.Following Flannigan (2009a, b), we consider an increase of 50 % of BC and other aerosols emitted by fire, together with a 1-month extension of the fire season in the Northern Hemisphere (starting 15 days prior and extending 15 days after the fire season of the present-day).S3 and S4 emission variations are applied to sulphate, BC and OC.S2, S3 and S4 experiments consist of a pair of 11-yr simulations, with initial conditions slightly modified in one of them, to be able to analyze 20 yr of model output, as 10 yr would clearly be insufficient to make comparisons statistically robust.In addition, to evaluate in more details the impact of the future aerosol emissions changes without considering atmospheric feedbacks, we realized three more experiments nudged toward our first 2050-2060 simulation: S2 N, S3 N and S4 N all have winds nudged toward S2, each of them using the same aerosol emissions as respectively S2, S3 and S4.Note that S2 N has been nudged toward itself (S2).This has been done to analyze the difference between simulations induced by the aerosol emissions change and not by the nudging itself.
Current BC emissions are particularly intense over the main industrialized regions of the Northern Hemisphere (Fig. 1a) with 2878 Gg yr −1 of BC emitted north of 30 • N in the S1 simulation.Regarding the difference between S2 and S1 (Fig. 1b), we diagnose that BC emissions are expected to significantly decrease over the major parts of industrialized areas in RCP8.5 (−1588 Gg yr −1 ), except in some regions of Central Introduction

Conclusions References
Tables Figures

Back Close
Full hypothesis that increases in a country's wealth are accompanied with the introduction of new technologies to reduce emissions.This being said, the RCP8.5 projections indicate an increase of emissions over the oceans, associated to an increase in ship and air traffic, which appears inevitable (Eyring et al., 2005;Sovde et al., 2007).Figure 1c shows the increase in BC emissions estimated by Corbett et al. (2010) consequent to the evolution of ship traffic over the Arctic Ocean which could take place in addition to the RCP8.5 emissions for 2050.Note that we consider a diminution of shipping emissions for current routes, as Arctic new routes would partially replace current ones (Corbett et al., 2010).For this reason, the total difference in emissions with the S2 simulation is very small (only +3.9 Gg yr −1 ).Finally we show in Fig. 1d the increase in BC emissions associated to the idealized lengthening (+15 days before and between the fire season) and intensification (+50 %) of biomass burning season applied on top of the RCP8.5 emission scenario (+236 Gg yr −1 north of 30 • N).Note that biomass burning emissions are assumed to be constant during all of the 21st century in the RCP8.5 scenario.

Results
We computed the mean number of days per year with snow at the surface (MNDWS) in all of our simulations as an indicator of the effects of BC changes on snow.Note that dust emissions were constant for all the simulations.In the following, we will not discuss the dust effects on snow.Figure 2a and b represent the MNDWS as observed (NSIDC, 2008) and modelled in our present-day control simulation S1, respectively.The MNDWS ranges from several days at 30 • N to almost a complete year north of 75 • N. The goal of our study is not to analyse in detail the ability of our GCM to describe the snow cover, as we will focus more on the analysis of sensitivity experiments with this GCM.Nevertheless, looking at the Root Mean Square Error (RMSE) between modelled and observed MNDWS (Fig. 2c), we see that our model describes quite well the snow cover duration over flat areas (RMS varying between 5 and 20).This is not

TCD Introduction Conclusions References
Tables Figures

Back Close
Full the case in mountainous areas like the Himalayas, the Altay Mountains, the Alps and the Rocky mountains where the RMSE generally exceeds values of 40 and can reach values of 300 days.As a consequence, we have to be very careful when we draw conclusions from the analysis of our simulation in these regions.Such huge errors are clearly due to the coarse resolution of our model, which does not allow a correct representation of the complex topography of these mountain ranges.Note that we did not consider the number of snow days over glaciers, icecaps or sea ice in our study.We discarded as well snow cover variations modelled in grid-cells located just next to icecaps (Greenland) since the representation of these icecaps is also not accurate due to the coarse spatial resolution of our model.
In the following, we discuss the difference of MNDWS between our different simulations.The statistical significance was estimated using a two-sample t-test.This statistical test is applied to validate the hypothesis that the mean of two simulations are different at the 95 % significance level.All areas with statistically significant differences are shaded in grey on Figs. 3 to 7. Regarding present-day conditions, considering the influence of BC deposition on snow albedo induces a decrease of the MNDWS that is statistically significant over a major part of the continents of the Northern Hemisphere (Fig. 3a, difference S1B-S1).This decrease lies within a range of 1 to 10 days over large areas of Eurasia and Northern America.Regarding future conditions, there is a significant decrease of the MNDWS in the S2 simulation for 2050 (Fig. 3b).This reduction is statistically significant, and ranges from 10 to 100 days in most parts of northern continental areas.Due to global warming forced by greenhouse gases, the beginning of the snow-accumulating season (respectively, the beginning of the snow-melting season) is modelled with ORCHIDEE coupled to LMDZ to occur later in autumn (resp.also driven in the model by snowfall variations.As an example, the snow cover duration is less reduced in Eastern Siberia than in Scandinavia, because snowfall is modelled to increase in Eastern Siberia in the middle of the 21st century.We found also a slight increase of the MNDWS compared to present-day over some northern parts of China and over the USA, also induced by a local increase in snowfall for the modelled LMDZ climate in 2050.However, we have to be very careful with this last result, as it concerns mountainous areas, where the GCM coarse resolution cannot provide accurate results as explained above. Considering an increase in aerosol emissions from Arctic ships or from biomass burning in our 2050-2060 nudged experiment induce MNDWS variations quasi equal to zero (see Fig. 3c and d, showing respectively MNDWS differences S3 N-S2 N and S4 N-S2 N).It clearly means that the snow albedo changes associated with this possible increase in aerosol emission is negligible in comparison with the snow albedo changes induced by the present day aerosol emissions in the Northern Hemisphere.However, we have to keep in mind that all of these present-day and future sensitivity experiments were nudged, a process that limits atmospheric feedbacks: these experiments allow to quantify the changes of snow cover duration induced by the aerosol effects on snow albedo, strongly minimizing both the effect of aerosols when they are in the atmosphere and the temperature changes induced by the snow cover variations.The complete effect of aerosols can be evaluated through simulations performed without nudging, as it was done for experiments S3 (with an increase in arctic ship traffic) and S4 (with an increase in biomass burning emissions).Nevertheless, we have to keep in mind that all of these future experiments used the same prescribed SST, which cancel the feedbacks which could be generated through interactions with the ocean.Since our study focuses on the continental response to a continental forcing, the analysis presented here should not be too much affected.Figure 4a and b  simulation performed with the standard RCP aerosol emissions (S2).Note that these variations of MNDWS are not statistically significant according to our two-sample t-test over the major part of the Northern Hemisphere.In other words, it means that the signal induced by the changes of aerosol emissions is too low to affect the highly variable coupled land-atmosphere system.Nevertheless, we obtained a statistically significant decrease of MNDWS in Quebec and in Siberia, both in simulation S3 and S4.These MNDWS local decreases reach 10 days averaged over the decade-long simulation of the 2050s.

Discussion
From the analysis of our nudged and not nudged experiments, we estimate that the possible increase in aerosol emissions from ships or boreal fires will not affect significantly the snow cover directly from snow darkening effects.However, this conclusion may not hold if we had also accounted for the atmospheric effects of aerosols.These effects are however very difficult to quantify: Shindell and Faluvegi (2009) showed that the patterns of temperature response and aerosol radiative forcing do not correspond on a regional basis.The difficulty to answer these complex questions is reinforced by the fact that ships emit different aerosol species (Balkanski et al., 2010), which have differentiated impacts on the climate system: they emit BC, an aerosol which absorbs solar radiation, warming its environment, but they also emit large amount of sulphate, an aerosol which strongly scatter solar radiation, cooling locally the atmosphere via direct end indirect effects (Lohmann, 2005).The sign of the radiative forcing induced by biomass burning, which also emits both BC, OC and sulphate depends also on the height at which the particles are transported (Abel et al., 2005).In front of all these complex questions, we discuss in the following when and how the MNDWS can be affected by increased ship and biomass burning aerosol emissions.
Both the biomass burning lengthened season and increased Arctic ship traffic emissions alternative scenarios produce very low emissions in winter.In summer, the Introduction

Conclusions References
Tables Figures

Back Close
Full Northern Hemisphere experiences a reduced snow cover.During fall, when solar radiation is considerably reduced compared to summer, both atmospheric aerosols and aerosols deposited on snow surface have a weak impact on snow cover (Flanner et al., 2009).Spring is the season when the Arctic atmosphere experiences the most pollution (e.g.Shaw et al., 1995;M én égoz et al., 2012).For all of these reasons, although summer is the period when aerosol concentrations from ship traffic and biomass burning are the largest, it is during the spring that we find the largest significant MNDWS changes associated to aerosol emissions considered in experiments S3 and S4.The significant spring aerosol emissions are concomitant with large residual snow cover over continental regions of the Northern Hemisphere, and thus have the potential to amplify regional warming.This is why we focus the following analysis on the interactions between snow and aerosols during the spring season.

BC deposition on snow
Present-day modelled BC spring deposition reaches 50 mg m −2 month −1 in Europe and Northern America, and exceeds 100 mg m −2 month −1 over South-East Asia (Fig. 5a).
Typical deposition values modelled in the pan-Arctic continental area range between 0.1 and 10 mg m −2 month −1 .In simulation S2 for 2050, a drastic decrease in BC deposition is obtained over the whole Northern Hemisphere (Fig. 5b), with the exception of central Asia and Alaska (where anthropogenic emissions are increasing in the RCP8.5 scenario compared to nowadays level, see Fig. experiment.These changes are certainly responsible for the modelled spatial variations of aerosol deposition during springtime.Note that this very weak signal is not statistically significant, indicating that the increase of ships emissions only generated "noise" in the aerosol spring deposition signal of our sensitivity experiment S3.Such response can be therefore mainly explained by natural variability.By contrast with S3, the earlier fire season considered in simulation S4 causes a significant increase in BC spring deposition over both North America and North Asia (Fig. 5d).The total increase of BC continental deposition in the S4 simulation represents 21 Gg month −1 .Regarding spring aerosol deposition, we can conclude that the MNDWS changes modelled in the S3 experiment is clearly not induced by snow darkening effects by aerosols.They are more due to aerosols when they are in the atmosphere, and to all the possible associated atmospheric feedbacks.Regarding S4 spring aerosol deposition, it is possible that snow darkening effect of BC have impacted the MNDWS via atmospheric feedbacks.

Spring Snow Water Equivalent (SWE)
During the spring, the present-day SWE ranges from 500 to 2000 mm in mountainous areas such as the Rocky Mountains, the Scandinavian mountains, the Ural Mountains or over Kamchatka (Fig. 6a).Elsewhere, over high latitudes continental areas, it takes values on the order of 100 mm.Considering BC deposition on snow in the present-day conditions (S1B-S1) induces only a small SWE decrease over large part of Eurasia an Northern America ranging from 0 to 10 mm (Fig. 6b).However, in a few locations of Western America and Scandinavia, this decrease takes larger values, exceeding 100 mm.The strongest BC induced decrease in present-day SWE appears in regions where the SWE is generally elevated in spring.Overall, spring SWE is modelled to be much lower in the RCP8.5 2050 scenario than under present-day conditions, and the modelled SWE decrease reaches up to 50 % over the major part of snow-covered areas (Fig. 6c).There are very few regions where spring SWE is modelled to increase in S2 compared to S1B, and these exceptions are North Eastern Canadian Islands, the Himalayan region and small parts of Northern Eurasia.An enhancement of ship Introduction

Conclusions References
Tables Figures

Back Close
Full traffic in the Arctic is predicted to induce an extra decrease of the SWE in Alaska, in the Canadian shield, and in large parts of Northern Eurasia, ranging from 10 to 100 mm (Fig. 6d), and in the Baffin Island, reaching 10 mm.In the scenario S4 with an earlier spring biomass burning activity, spring SWE is modelled to decrease in many parts of the continental pan-Arctic areas, by up to 50 mm, except in Baffin Island and in very small regions of Northern Eurasia (Fig. 6e).However, these modelled extra SWE changes in simulations S3 and S4 are not statistically significant according to a twosample t-test, indicating that the signal of the local aerosol emissions taken into account is difficult to be characterized given the large amount of natural climate variability, and the fact that local emissions play a second order role (S3-S2 and S4-S2) compared to the fist order effect of GHG forced future warming effects on SWE (S2-S1B).
The present-day SWE decrease induced by aerosol deposition is quite smaller than the decrease modelled in 2050 under the RCP8.5 scenario (see Fig. 6b and c).The decrease of SWE expected in 2050 is due to the temperature increase associated with the greenhouse gas radiative forcing.This result clearly shows that the drastic reduction of BC deposition in the Northern Hemisphere in 2050 (Fig. 5b) is clearly not sufficient to counteract the decrease of SWE induced by greenhouse gas radiative forcing and its associated temperature increase (Fig. 6c).As explained previously, there are almost no changes in aerosol deposition in the simulation S3 with enhanced ship emissions.The modelled changes in MNDWS and SWE are therefore due to atmospheric aerosol effects, which can experience atmospheric feedbacks.For the S4 simulation with enhanced biomass burning in spring, there is a significant increase of aerosol deposition, which may explain a part of the MNDWS and SWE diminutions in some regions of Northern Eurasia and Northern America.This assumption appears very likely where the MNDWS variations are statistically significant, in North-eastern America as in central and eastern Siberia.However, the SWE variations are generally not statistically significant, and there is no clear correlation between BC deposition and snow cover variations.Therefore, it is likely that part of the SWE changes are also consecutive to surface energy balance changes or to snowfall variations in our simulations.Introduction

Conclusions References
Tables Figures

Back Close
Full

Spring snowfalls
Present-day spring snowfalls are widespread over a large part of the Northern Hemisphere continents (Fig. 7a).In our present-day simulation, snow albedo diminution induced by BC aerosol deposition leads to a slight but statistically significant snowfall diminution (Fig. 7b).A large part of the spring decrease in SWE between 2050 and present-day simulations (Fig. 6c) can be explained by this snowfall feedback (Fig. 7c).
In most part of the spring snow-covered area of the Northern Hemisphere, snowfall decreases by 50 % (see Fig. 7a and c) in S2 compared to S1B.This is mainly due to temperature rise, which transforms snowfall into rainfall.We find only few and small areas, like North Eastern Canadian Islands, parts of the Himalayan region and very small parts of Northern Eurasia where snowfall increases.However, these increases may explain the SWE increases modelled in the same regions.
Based upon the sensitivity experiments S2, S3 and S4, we are able to evaluate the impact of an aerosol emission change in a 2050 scenario.In simulations S3, the spring SWE change exhibits a pattern similar to snowfall change in many continental areas of the Northern Hemisphere, with a general decrease in the pan-arctic area, except in small areas like Baffin Islands and other Northern Canadian islands (Fig. 7d).Therefore, we can assess that the atmospheric perturbations induced by enhanced ship traffic BC emissions in the Arctic induce a small decrease of snowfall over large area of the boreal continents.Even if these variations are not statistically significant according to a two-sample t-test, they partly contribute to the decrease of SWE modelled in the same region.However, it is very difficult to estimate which physical processes link snowfall variations to BC aerosol emissions change, since aerosols from ships contain both absorbing and reflective species which have complex interactions with the atmosphere (Balkanski et al., 2010).Regarding the S4 simulation, we can also assess that the snowfall decreases which take place in the major part of Northern America (Fig. 7e), in North-Eastern Europe and in North-east Asia are responsible for part of the modelled decrease of both MNDWS and SWE in these regions.However, this assumption is not

Conclusions References
Tables Figures

Back Close
Full verified in Northern Central Siberia, where we modelled an increase of snowfall but a decrease of the SWE and the MNDWS.In this region, the SWE decline is certainly induced by an aerosol forcing.It may be due both to a decrease of the snow albedo via aerosol deposition, and to a warming of the atmosphere associated to an increase in the atmospheric concentration of BC.

Conclusions
The snow-cover changes induced by aerosol emissions were evaluated in the boreal continental area both for the present-day and for the middle of the 21st century.The following eight experiments were carried out: two present-day simulations, with one of them not considering the snow albedo variations induced by aerosol deposition, and six 2050-2060 simulations based upon the RCP8.5 gas and aerosol anthropogenic emission inventory.We estimate that current aerosol emissions directly cause a decrease of the MNDWS ranging between 0 and 10 days in large areas of the boreal region.This "snow darkening effect" is essentially due to the BC deposition during the spring, a period of the year when the remaining of snow accumulated during the winter is exposed to both strong solar radiation and large amount of aerosol deposition.This deposition over continents represents 222 Gg month −1 of BC north of 30 • N.
The projected drastic decrease of the anthropogenic aerosol emissions from the RCP scenarios for the middle of the 21st century in the Northern Hemisphere may limit the decrease of snow albedo due to absorbing aerosol deposition.But this response is very much dependent on the quality of the emission scenarios, as no inflexion in BC emissions over Asia has been observed in the past decades.Nonetheless, a major part of snow-cover in the Northern Hemisphere will experience a significant reduction under the GHG forced warming.By comparison with present-day conditions, the MNDWS was found to be reduced by 10 to 100 days over the major part of the continental regions of the Northern Hemisphere by the middle of the 21st century.The main cause

Conclusions References
Tables Figures

Back Close
Full for this decrease is a temperature rise that substitutes snow to rain over several regions and accelerates melting.These conclusions have been reached with a future scenario that considers strong increases in greenhouse gases concentrations.The decrease of the aerosol impact on snow-cover should be relatively more important for a scenario with lower greenhouse gases concentrations.
Considering a significant additional increase in ship traffic in the Arctic by the mid 21st century does not lead to significant changes of the aerosol deposition over snowcovered areas in the most sensitive period for a positive climate feedback, springtime.Therefore, the MNDWS is clearly not affected by snow darkening effects associated to these Arctic ship emissions.This result has been demonstrated using a simulation nudged toward the observed atmosphere, to quantify how aerosol deposition could affect directly the snow cover.In an experiment considering such an increase of ship emissions without nudging toward atmospheric reanalyses, we simulated some changes of the MNDWS.Ships emit absorbing aerosols like BC and to a lesser extent OC, but in comparison a lot more sulphur dioxide, which strongly scatters the incoming solar radiation, thereby cooling the atmosphere.Modifying the atmospheric energy balance by accounting for these aerosols affects the atmospheric circulation and the precipitation pattern.In this experiment, the MNDWS changes are generally not statistically significant in boreal continents, except in the Quebec and in the West Siberian plains, where the MNDWS decrease from 5 to 10 days.
Biomass burning activity proportionally emits more BC and OC aerosol and much less sulphate compared with ship traffic.We modelled a significant increase in BC spring deposition that exceeds 1 mg m −2 month −1 over large parts of America and Eurasia in a 2050-2060 simulation that take into account forest fires that are 50 % stronger and are projected to occur 2 weeks earlier than at present.This increase of BC spring deposition represents 21 Gg month −1 on continents located north of 30 • N.However, with such emissions, we do not simulate a reduction of the MNDWS in an experiment performed with winds nudged toward atmospheric reanalyses.This demonstrates that our biomass burning emission scenario does not induce a significant reduction Introduction

Conclusions References
Tables Figures

Back Close
Full of the snow cover via "snow darkening effects".However, considering all the aerosol forcings and atmospheric feedbacks in an experiment performed without nudging, enhanced fire activity induces a significant decrease of the MNDWS reaching a dozen of days in Quebec and in Eastern Siberia.Due to the snow-albedo feedback, the Arctic is a region very sensitive to climate change.As a consequence of this feedback, Flanner et al. (2009) showed that absorbing aerosol emissions reduced the springtime snow cover as much as anthropogenic greenhouse gases since the pre-industrial period.Consequently, limiting aerosol emissions appears as essential as limiting greenhouse gases emissions to slowdown the snow cover decline observed over the Northern Hemisphere.Foreseeing the possible emissions scenarios in the 21st century, one can envisage for strong aerosol reductions in most industrialized region over the Northern Hemisphere with the introduction of advanced technologies in controlling emissions.However, increases in the emissions and concentrations of greenhouse gases that are projected in most scenarios are expected to significantly reduce the snow cover in the middle of the 21st century.It appears very challenging to estimate accurately the snow cover changes induced by the possible changes in aerosol emissions in the Arctic and the boreal region because of the complex processes linking aerosol forcing, atmosphere response and snow cover dynamics.However, we predict that the likely future aerosol emissions from ships traffic over the Arctic region or an increase in biomass burning will play a minor role in the reduction of continental snow cover area trough snow darkening effects at high Northern latitudes.We have not attempted to predict future changes in sea ice due to these effects but these may be significant.Introduction

Conclusions References
Tables Figures

Back Close
Full   x2) in the period means that the simulation was performed a second time with a slightly modified initial conditions to get 20 yr of simulation as 10 yr would clearly be insufficient to make comparisons statistically robust.Note that all simulations were made with prescribed Sea Surface Temperature (SST, observed for present-day simulations, or simulated from a previous coupled ocean-atmosphere model simulation for future periods).
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | However, theConway et al. (1996) study was based upon experiments with particularly high rates of snow melting since they were performed during summer at altitudes around 2000 m over the Northern United States.More recent observations by Aamaas et al. (2011) in Spitsbergen showed that BC aerosols tend to stay at the surface of the Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | earlier in spring) in most snow-covered northern regions.A negative trend of MNDWS has already been observed during the last decades (e.g.D éry et al., 2007; Roesch et al., 2006; Mote et al., 2005).Moreover, Hosaka al. (2005) and Brutel-Vuilmet et al. (2012) expect an acceleration of this phenomenon into the 21st century.Similar to the results reported by Hosaka et al. (2005), we found that the snow cover changes are Discussion Paper | Discussion Paper | Discussion Paper | show that without nudging the variations in MNDWS with enhanced ship and fire emissions can be positive or negative depending upon the region.They are spatially variable, and reach values ranging from −10 to +10 days per year in comparison with our 2050-2060 Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | 1b).On average over all the continental surfaces of the Northern Hemisphere, this decrease represents half of the present-day spring deposition (decrease of 110 Gg month −1 for a present-day total of 222 Gg month −1 , north of 30 • N).The simulation performed with extra high ship emissions in the Arctic (S3) does not induce significant changes of BC deposition in spring (Fig. 5c) in comparison to the S2 2050 simulation.This is due to the fact that the additional Arctic ship emissions are mainly enhanced in summer, when ships use alternate Arctic routes.Yet, these enhanced ship emissions modify the atmospheric circulation and precipitation via the atmospheric aerosol radiative forcings in our sensitivity Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Table 1 .
Period, aerosol emissions and description of the nudging protocol for our 8 simulations.(