Articles | Volume 20, issue 13
Atmos. Chem. Phys., 20, 8139–8156, 2020

Special issue: Arctic mixed-phase clouds as studied during the ACLOUD/PASCAL...

Atmos. Chem. Phys., 20, 8139–8156, 2020

Research article 13 Jul 2020

Research article | 13 Jul 2020

Combining atmospheric and snow radiative transfer models to assess the solar radiative effects of black carbon in the Arctic

Combining atmospheric and snow radiative transfer models to assess the solar radiative effects of black carbon in the Arctic
Tobias Donth1, Evelyn Jäkel1, André Ehrlich1, Bernd Heinold2, Jacob Schacht2, Andreas Herber3, Marco Zanatta3, and Manfred Wendisch1 Tobias Donth et al.
  • 1Leipzig Institute for Meteorology (LIM), University of Leipzig, Leipzig, Germany
  • 2Leibniz Institute for Tropospheric Research, (TROPOS), Leipzig, Germany
  • 3Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research (AWI), Bremerhaven, Germany

Correspondence: Evelyn Jäkel (


The magnitude of solar radiative effects (cooling or warming) of black carbon (BC) particles embedded in the Arctic atmosphere and surface snow layer was explored on the basis of case studies. For this purpose, combined atmospheric and snow radiative transfer simulations were performed for cloudless and cloudy conditions on the basis of BC mass concentrations measured in pristine early summer and more polluted early spring conditions. The area of interest is the remote sea-ice-covered Arctic Ocean in the vicinity of Spitsbergen, northern Greenland, and northern Alaska typically not affected by local pollution. To account for the radiative interactions between the black-carbon-containing snow surface layer and the atmosphere, an atmospheric and snow radiative transfer model were coupled iteratively. For pristine summer conditions (no atmospheric BC, minimum solar zenith angles of 55) and a representative BC particle mass concentration of 5 ng g−1 in the surface snow layer, a positive daily mean solar radiative forcing of +0.2 W m−2 was calculated for the surface radiative budget. A higher load of atmospheric BC representing early springtime conditions results in a slightly negative mean radiative forcing at the surface of about 0.05 W m−2, even when the low BC mass concentration measured in the pristine early summer conditions was embedded in the surface snow layer. The total net surface radiative forcing combining the effects of BC embedded in the atmosphere and in the snow layer strongly depends on the snow optical properties (snow specific surface area and snow density). For the conditions over the Arctic Ocean analyzed in the simulations, it was found that the atmospheric heating rate by water vapor or clouds is 1 to 2 orders of magnitude larger than that by atmospheric BC. Similarly, the daily mean total heating rate (6 K d−1) within a snowpack due to absorption by the ice was more than 1 order of magnitude larger than that of atmospheric BC (0.2 K d−1). Also, it was shown that the cooling by atmospheric BC of the near-surface air and the warming effect by BC embedded in snow are reduced in the presence of clouds.

1 Introduction

Black carbon (BC) aerosol particles, which mostly originate from incomplete combustion of organic material (Bond et al.2013; Petzold et al.2013), absorb and scatter solar radiation in the visible wavelength range and, therefore, influence the atmospheric solar radiative energy budget. The manifold sources of BC particles and their atmospheric transport paths have been studied extensively (Law et al.2014). However, the source strengths of the emissions are hard to quantify, which makes it challenging to quantify the transport of BC particles into the Arctic by simulations (Stohl et al.2013; Arnold et al.2016; Schacht et al.2019). Major sources of BC particles are forest fires, industrial activities, and traffic-related emissions, which are the main factors in lower latitudes: northern parts of Europe, America, and Siberia. The BC particles emitted at the surface of the midlatitudes are lifted and transported into the Arctic, where they can stay for several days and longer (Liu et al.2011). Contrarily, particles produced locally in the Arctic through ship traffic emissions, flaring from the oil industry, or other ground-based activities settle down quickly on the surface and may alter the radiation budget within the snowpack (Bond et al.2013). Today, local sources are only a minor component. In the future, a strong intensification of the ship traffic in the Arctic Ocean and further polluting human activities are expected (Corbett et al.2010). Still, the direct radiative impact by these future additional BC particle emissions is assumed to be of minor importance (Gilgen et al.2018).

The BC magnitude of the atmospheric particle mass concentrations (in units of nanograms per cubic meter) depends on the season and general meteorological conditions. In the case of BC particle plumes reaching the Arctic by long-range transport, atmospheric concentrations of up to 150 ng m−3 were observed (Schulz et al.2019). Sharma et al. (2013) compared atmospheric BC particle mass concentrations measured during different Arctic campaigns. They identified large differences depending on region and season. Measurements in spring 2008 covering Alaska and northern Canada showed values above 200 ng m−3 in higher altitudes, while in spring 2009 more pristine air masses were encountered, showing BC particle mass concentrations of less than 100 ng m−3 integrated over the entire vertical column.

To quantify the amount of BC in a snowpack volume, the BC mass concentration (nanograms of BC in 1 g of snow) is used commonly. Typical values observed in Greenland range between 1 and 10 ng g−1 and in the Canadian Arctic between 5 and 20 ng g−1, and in the northern parts of Russia values may reach 100 ng g−1. Table 1 summarizes observational data of measured BC mass concentrations in snow for different Arctic regions, as reported by Doherty et al. (2010), Forsström et al. (2013), and Pedersen et al. (2015). The numbers given in Table 1 were derived from different measurement methods. More precisely, thermal–optical techniques were applied in Forsström et al. (2013) and Pedersen et al. (2015) and provide the elemental carbon (EC) mass concentration, while filter transmission methods result in BC concentrations (Doherty et al.2010). As a consequence of the different measurement methods, the ratio of the BC to EC concentration in snow can reach values of 1.3 as reported by Dou et al. (2017). A full discussion of the EC  BC terminology can be found in Petzold et al. (2013).

Due to the absorption of solar radiation, BC particles may contribute to the currently ongoing drastic Arctic climate changes (called arctic amplification; e.g., Wendisch et al.2017). The absorption effect can add to the warming of the atmosphere or the snowpack, when the BC particles are suspended either in the air or embedded in the snow. Furthermore, the BC particles may lead to a reduction of the snow surface albedo if the BC sediments on or into the snowpack (Sand et al.2013). Warren (2013) estimated a decrease of 2 % in snow albedo in the visible spectral range for a snowpack with a BC mass concentration of 34 ng g−1, which corresponds to the maximum value observed on the Greenland ice sheet (Doherty et al.2010). More typical BC mass concentrations in Arctic snow range between 5 and 20 ng g−1 (Table 1), which would lead to a reduction of the snow surface albedo of around 1 %. For typical Arctic summer conditions with a downward irradiance of 400 W m−2 at the surface, a snow surface albedo reduction by 1 % would cause an additional absorption of solar radiative energy of 4 W m−2 (Flanner et al.2007). As a further consequence, the absorption by BC particles supports the melting of snow and increases the snow grain size due to an enhanced snow metamorphism, leading to further reduction of the surface albedo. The increase in the snow grain size also feeds back to the absorption by BC particles, which is more efficient for larger snow grain sizes (Warren and Wiscombe1980).

BC particles suspended in the atmosphere influence the absorption and scattering of the incoming solar radiation. If atmospheric BC particles are located in high altitudes, enhanced backscattering and absorption of incoming solar radiation by the BC layer lead to a reduction of the solar radiation reaching the surface. At the same time, the absorbed radiation warms the atmospheric BC layer. In extreme cases, the absorption due to atmospheric BC particles can affect the thermodynamic stability of the BC-containing atmospheric layer (Wendisch et al.2008). The radiative heating of the lofted BC layers and the local cooling of the surface may enhance the thermodynamic stability of the Arctic boundary layer over the snow- and ice-covered areas (Flanner2013).

Several regional and global climate models account for the opposite radiative effects of atmospheric BC particles and snow-embedded BC particles (Samset et al.2014). However, estimates of the total net forcing rely on the accuracy of the distribution of the BC particles assumed in the particular model. Samset et al. (2014) compared 13 aerosol models from the AeroCom Phase II; all of them included BC. They found that modeled atmospheric BC particle mass concentrations often show a spread over more than 1 order of magnitude. In remote regions, dominated by long-range transport, these models tend to overestimate the atmospheric BC concentrations compared to airborne observations. On the other hand, an underestimation of deposition rates induces a lower BC mass concentration in snow (Namazi et al.2015). While this may introduce significant local and temporal uncertainties of the BC concentration and related radiative effects, long-term trends and mean multi-model results are representative for Arctic-wide observations (Sand et al.2017).

Doherty et al. (2010)Doherty et al. (2010)Doherty et al. (2010)Doherty et al. (2010)Doherty et al. (2010)Forsström et al. (2013)Pedersen et al. (2015)Pedersen et al. (2015)Pedersen et al. (2015)Pedersen et al. (2015)Pedersen et al. (2015)

Table 1Values of the black carbon mass concentration in snowpack observed in different regions and seasons in the Arctic. Note that Pedersen et al. (2015) and Forsström et al. (2013) derived the mass concentration of elemental carbon applying a thermal–optical measurement method.

Download Print Version | Download XLSX

Most previous studies quantifying the radiative impact of BC particles focused either on estimates of cooling or warming effects in the atmosphere (e.g., Wendling et al.1985; Samset et al.2013) or on radiative effects of BC in the snow surface layer (Dou and Cun-De2016). In contrast, this paper combines both effects by iteratively coupling radiative transfer simulations in both compartments, the atmosphere and the snowpack. On the basis of measured Arctic BC particle mass concentrations for spring and summer months, the instantaneous radiative forcing of BC particles embedded in the snow surface layer and in the atmosphere was quantified for specific cases. Here, the instantaneous radiative forcing refers to the change of the surface radiation budget caused by the presence of BC particles. With the help of the coupled model, the interaction of radiative effects in the atmosphere and the snowpack was considered. In particular, the role of clouds in the cooling or warming effect caused by BC particles was examined. Due to the fact that clouds enhance the atmospheric multi-scattering between surface and cloud layer, but also enhance the surface albedo (Choudhury and Chang1981), it is expected that clouds also alter the radiative impact by BC particles. To our knowledge, this interaction was not explicitly discussed in previous publications.

The radiative transfer simulations used in this study were based on airborne observations of atmospheric BC concentration in the Arctic, which were taken during three field campaigns in the European and Canadian Arctic. The applied models and observations are introduced in Sect. 2. Section 3 discusses the radiative forcing of BC particles on the surface solar radiative budget. Vertical profiles of heating rates in the atmosphere and in the snowpack are presented for clean and polluted conditions. To estimate the impact of BC particles, effective heating rates are calculated by separating the BC radiative effect from the total heating rates.

2 Configuration of radiative transfer simulations and iterative model coupling

2.1 BC profiles from aircraft campaigns

The input for the radiative transfer simulations was adapted to campaign-specific conditions. The atmospheric BC particle mass concentrations were derived from airborne measurements with the Single Particle Soot Photometer (SP2; Moteki and Kondo2007). Measured profiles of the atmospheric BC were used from three aircraft campaigns representing typical cases with higher BC concentrations (polluted case) in early spring with low Sun and lower BC concentration (pristine conditions) in early summer during the polar day. The Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) spring campaign was performed in April 2008 (Jacob et al.2010; Matsui et al.2011). The aircraft operation of ARCTAS mainly took place in northern Alaska and the Arctic Ocean. Similar SP2 measurements were performed during the Polar Airborne Measurements and Arctic Regional Climate Model Simulation Project (PAMARCMiP) campaigns (Herber et al.2012; Stone et al.2010). In this paper, measurements from the PAMARCMiP 2018 observations conducted from 10 March to 8 April 2018 were analyzed. The research flights, starting from Station Nord, Greenland, were performed above the sea ice in the Arctic Ocean north of Station Nord and the Fram Strait. In contrast to both spring campaigns, the Arctic CLoud Observations Using airborne measurements during polar Day (ACLOUD) campaign was conducted in early summer 2017 characterizing the atmosphere over the Arctic Ocean north and west of Svalbard (Wendisch et al.2019; Ehrlich et al.2019). ACLOUD was coordinated with the Physical Feedbacks of Arctic Boundary Layer, Sea Ice, Cloud and Aerosol (PASCAL) cruise of the research vessel Polarstern which provided a ground-based characterization of snow properties (Wendisch et al.2019).

Mean vertical profiles of the measured atmospheric BC particle mass concentrations averaged for each of the three campaigns (ACLOUD, ARCTAS, and PAMARCMiP), are shown in Fig. 1a. The conditions between the individual flights were highly variable (see the standard deviation of each layer in Fig. 1). ACLOUD shows rather low mean BC concentrations, which do not exceed 30 ng m−3. During PAMARCMiP, the background concentrations were similarly low, with the exception of measurements at about 5 km altitude, where more than 100 ng m−3 was recorded. For ARCTAS observations, conducted at lower latitudes, significantly higher BC concentrations of up to 150 ng m−3 were observed. Similar to PAMARCMiP, the maximum concentrations were observed at about 5 km altitude, indicating that the BC particles were linked to long-range transport (Nakoudi et al., 2020). Besides the differences in atmospheric BC concentrations, the range of the daily solar zenith angle (SZA) and, thus, the available incoming solar radiation varied significantly for the three campaign periods. When analyzing the radiative impact of BC on the basis of daily averages, the magnitude of the incident solar radiation and the length of the day play a major role. While the early summer conditions of ACLOUD were characterized by the polar day and SZA between 55 and 78, during ARCTAS the available incoming solar radiation was lower due to lower values of the SZA (minimum at noon of 62.5). PAMARCMiP was conducted in the most northern region and earlier in the year, such that the Sun was about 9.5 h below the horizon and the minimum SZA was 79 at noon. Table 2 summarizes the key characteristics of the three analyzed data sets.

Jacob et al. (2010)Ehrlich et al. (2019)Herber (2019)

Table 2Region, period, solar zenith angle range, and maximum BC particle mass concentration and mean optical depth of BC at the 500 nm wavelength characterizing the three data sets obtained within ARCTAS, ACLOUD, and PAMARCMiP.

Download Print Version | Download XLSX

Therefore, the constructed profiles of PAMARCMiP and ACLOUD are assumed to be representative for Arctic early spring and early summer conditions, respectively.

Figure 1Mean profiles of atmospheric BC particle mass concentration (a) and relative humidity (b) averaged for each of the three campaigns (ACLOUD, ARCTAS, and PAMARCMiP) as used for the radiative transfer simulations. Horizontal bars indicate the standard deviation. The positions of the two implemented cloud layers (blue shaded area) are marked.


2.2 Atmospheric radiative transfer model

To simulate vertical profiles of the spectral upward and downward irradiance, the library for radiative transfer routines and programs (libRadtran, Emde et al.2016; Mayer and Kylling2005) was used (, last access: 10 June 2020). The model also provides the ratio of the direct to global irradiance fdir∕glo, which is required as a boundary condition of the snowpack radiative transfer model. As a solver for the radiative transfer equation, the Discrete Ordinate Radiative Transfer solver (DISORT) 2 (Stamnes et al.2000) routine running with 16 streams was chosen.

For the calculations, a plane-parallel atmosphere was assumed, which is justified for the Arctic conditions during the three campaigns. Using a pseudo-spherical geometry in libRadtran would change the broadband downward irradiance by less than 0.1 % (0.7 %) for a calculation with a SZA of 60 (75). The vertical resolution of the simulated irradiances was adjusted to the measured BC profiles, ranging between 100 m and 1 km. The spectral resolution of the simulations was set to 1 nm covering a wavelength range between 350 and 2400 nm. The extraterrestrial spectrum was taken from Gueymard (2004). The BC optical properties including the refractive index, density, extinction coefficient, single-scattering albedo, and scattering phase function from the OPAC aerosol database were applied (Hess et al.1998). Corresponding to the campaign average BC profiles, the range of the SZA values was set to values representing the campaign conditions (see Table 2).

The meteorological input for the model was based on standard profiles of trace gases, temperature, humidity, and pressure from Anderson et al. (1986). Sub-Arctic summer conditions were chosen for the early summer case (ACLOUD) and subarctic winter conditions for the winter and spring cases (ARCTAS, PAMARCMiP). The standard profiles were adapted to observations from radio soundings near the airborne observations or dropsondes released during the flights and represent the middle of the individual campaign periods. Figure 1b shows the profiles of relative humidity used for the simulations. PAMARCMiP was characterized by rather dry air. Only in the boundary layer was an average humidity up to 60 % often observed linked to boundary layer clouds. ACLOUD and ARCTAS showed a higher relative humidity at a higher altitude of up to 6 km, which indicates the influence of higher-level clouds.

To test the sensitivity of the BC radiative effects with respect to cloud occurrence, two cloud layers were synthetically included in the atmospheric profiles as illustrated in Fig. 1. The cloud layer properties were based on observations by Bierwirth et al. (2013), Leaitch et al. (2016), and Blanchard et al. (2017) to represent typical Arctic cloud conditions. A low-level liquid water cloud was placed between 500 m and 1.4 km representing the humid boundary layer observed during PAMARCMiP. The liquid water content increases from 0.1 g m−3 at cloud base to 0.3 g m−3 at cloud top; the cloud particle effective radius increased from 6 to 12 µm. The second cloud layer represents a thin ice water cloud and was positioned between 5 and 5.5 km representing the higher-level clouds observed during ACLOUD and ARCTAS. This thin cloud was assumed to be homogeneous with an ice water content of 0.006 g m−3 and an effective cloud particle radius of 40 µm, according to airborne measurements reported by Wyser (1998) and Luebke et al. (2013). Optical properties of the liquid cloud droplets were calculated from Mie theory, while the ice crystal optical properties are based on Fu (2007). The assumed cloud properties correspond to a cloud optical thickness of 15 for the water cloud and 0.2 for the thin ice cloud.

2.3 Snowpack radiative transfer model

The Two-streAm Radiative TransfEr in Snow model (TARTES,, last access: 6 February 2020) was used to simulate the radiative transfer through the snowpack (Libois et al.2013, 2014). In TARTES, the snow profile is constructed of a predefined number of horizontally homogeneous snow layers, which allows us to account for the stratification of the snowpack. To consider the single-scattering properties of each layer, the method described by Kokhanovsky and Zege (2004) is applied in TARTES. To solve the radiative transfer equation, the delta-Eddington approximation (Joseph et al.1977) is used. As a result, TARTES computes the spectral surface albedo and the profile of the irradiance within the snowpack. As the boundary condition, the SZA and fdir∕glo have to be predefined. For each of the snow layers, the optical and microphysical properties have to be given, such as the snow density (ρice), the specific surface area (SSA), and the snow grain shape parameters, which represents a mixture of different grains as suggested by Libois et al. (2013). Furthermore, the specific values of the so-called absorption enhancement parameter B=1.6 and the geometric asymmetry factor gG=0.85 were applied. The specific surface area can be translated into the optical snow grain size ropt by

(1) r opt = 3 ρ ice SSA .

TARTES allows us to consider impurities in each snow layer, which are characterized by the impurity type and mass concentration. The impurities are externally mixed and assumed to interact by Rayleigh scattering. To simulate a BC-containing snow layer, the complex refractive index and the density of BC particles given by Bond et al. (2013) are applied.

The input parameters of the snowpack model are summarized in Table 3. For the bottom layer, a soil albedo of 0.3 was assumed representing the reflection properties below the snowpack. The impact of the soil albedo on the albedo of the snow surface depends on the depth of the overlying snowpack. Sensitivity studies have shown that for snow depths of more than 20 cm the albedo of a snow surface is independent of the choice of the soil albedo below. In this study the snowpack depth was set to 1 m thickness. Reference simulations assuming a pristine homogeneous snow layer were performed. Simulations including BC impurities were based on BC particle mass concentrations summarized in Table 1 (Doherty et al.2010; Forsström et al.2013; Pedersen et al.2015) and observations during PASCAL and PAMARCMiP. For the simulations of a single homogeneous snow layer, typical BC particle mass concentrations of 5 and 20 ng g−1 were chosen. The default values of snow density and SSA were based on measurements during PASCAL and PAMARCMiP and were set to 300 kg m−3 and 20 m2 kg−1, respectively. To analyze the sensitivity of the snow surface albedo with respect to the snow grain size, SSA values of 5 and 60 m2 kg−1 were used. The vertical model resolution was set to 1 cm.

In addition to the simulations of a homogeneously mixed snow layer, a second model setup was used to consider a multilayer snowpack. Pit measurements in Greenland (Doherty et al.2010) identified typical multilayer structures, where BC accumulated in a melting layer approximately 10 cm below the surface. Referring to these measurements, the snowpack of the second model setup consists of three snow layers. The top layer is 5 cm and the BC-containing middle layer is 10 cm thick. The bottom layer below continues to 1 m depth. For this multilayer approach, BC was included in the middle layer, representing an aged melting layer in which impurities had accumulated (SSA = 20 m2 kg−1, snow density of 350 kg m−3, and a BC mass concentration of 15 ng g−1). The top layer was assumed to be of fresh and clean snow with SSA = 40 m2 kg−1, a snow density of 250 kg m−3, and a BC mass concentration of 2 ng g−1, representing measurements from the PASCAL campaign. The aged snow layer at the bottom was characterized by an enhanced snow grain size and density of SSA = 10 m2 kg−1 and ρice=450 kg m−3, respectively, and a BC mass concentration of 2 ng g−1.

Table 3Snowpack model setups for the single-layer and multilayer cases. The default SSA for the single-layer case is 20 m2 kg−1.

Download Print Version | Download XLSX

2.4 Iterative coupling

The surface albedo is an important boundary condition to simulate the radiative transfer in the atmosphere. It depends on the illumination conditions defined by the solar zenith angle, the spectral distribution of downward irradiance, and the ratio of direct-to-global irradiance (e.g., Wiscombe and Warren1980; Gardner and Sharp2010; Stapf et al.2019). The transition from cloudy to cloudless atmospheric conditions increases the direct-to-global ratio (fdir∕glo) and the contribution of short wavelengths to the broadband downward irradiance (Warren1982). Therefore, a cloud cover typically increases the broadband surface albedo. For example, simulations with TARTES assuming cloudless and cloudy conditions changed the broadband snow surface albedo from about 0.8 to 0.9 for a SZA of 60 and a snowpack (no impurities) characterized by SSA =20 m2 kg−1. As clouds absorb solar radiation mostly at wavelengths larger than 1000 nm, the shorter wavelengths, where BC particles strongly absorb solar radiation, become more relevant. Because of the significant surface–cloud interactions, the atmospheric and snowpack radiative transfer models need to be coupled, interactively. Therefore, an iterative method coupling libRadtran and TARTES via their boundary conditions, surface albedo, and direct-to-global ratio (fdir∕glo) of the incident radiation, was applied. Both parameters were transferred between the models as schematically illustrated in Fig. 2.

Figure 2Schematics of the coupling of TARTES (gray box) and libRadtran (blue box) by exchanging the spectral surface albedo and the direct-to-global ratio. The list of varied parameters addresses the variables which were changed between the different realizations. Only the iterated parameters fdir∕glo and αλ were adjusted within an individual iteration cycle.


In the first iteration step, only diffuse radiation was assumed (fdir/glo=0) to calculate the snow surface albedo by TARTES, which subsequently serves as input for the libRadtran simulations. Then a new spectral direct-to-global ratio representing the atmospheric conditions was calculated by libRadtran, which is in turn used to readjust TARTES, starting a revised iteration (n+1) to calculate a new spectral surface albedo αλ(n+1). This procedure was repeated until the deviation of the surface albedo calculated in the previous step (n) and calculated in the revised step (n+1) decreases below 1 %. Exemplarily, Fig. 3 illustrates the change of the spectral surface albedo for a cloudless case without atmospheric BC and a SZA of 60. The BC mass concentration in snow was set to 5 ng g−1. Two iteration steps were necessary in this particular example to match the 1 % termination criterion, which is a typical number for all studied cases. Starting with purely diffuse conditions allows faster calculations in cloudy cases. This quick convergence of the iteration enables the consideration of different cloud properties and atmospheric conditions and facilitates the calculation of the radiative effects of BC particles in the atmosphere and within the snowpack simultaneously. The assumption of a pure diffuse illumination in the initial run caused no significant difference of the calculated visible snow albedo to the first and second iteration step. In contrast, the iterated direct-to-global ratio adjusts the snow albedo in the near infrared because the direct fraction is quickly approaching unity in this spectral range.

Figure 3Change of the spectral snow albedo for cloudless conditions with a SZA of 60 due to the iterative adjustment by the coupled atmosphere and snow radiative transfer models. The initial run assumes a direct-to-global ratio of zero.


2.5 Quantities used to characterize the impact of BC particles

In the following, the surface radiative forcing of BC particles and profiles of heating rates are analyzed. The total radiative forcing at the surface ΔFtot is separated into the forcing of BC particles suspended in the atmosphere ΔFatm and the forcing of BC particles deposited in the snowpack ΔFsnow. ΔFsnow is defined by the difference of the net irradiance (downward minus upward solar irradiance) if BC is considered in the snow layer (Fnet,BC) and a clean reference case without BC in the snow layer (Fnet,clean). Similarly, ΔFatm is defined as the difference between the net irradiances derived for BC in snow and the atmosphere and the atmospheric BC-free reference case:

(2) Δ F i = F net , BC - F net , clean ,

with index “i” standing for “tot”, “atm”, or “snow”. For the separated forcings, Fnet,clean refers to either a clean atmosphere or a clean snow layer, while the other part does consider BC particles. The default case of a clean atmosphere used a BC mass concentration in the snow layer of 5 ng g−1. Vice versa, the default case of a clean snow layer assumed the atmospheric BC profile of the ACLOUD campaign. For ΔFtot, the clean reference assumed both a pristine atmosphere and pristine snow layer.

The calculation of atmospheric and snow heating rate profiles HR(z) (K d−1) was based on the net irradiances at the top (t) and bottom (b) of selected atmospheric or snow layer z, the layer density ρ(z), the specific heat capacity under constant pressure cp, and the layer thickness (ztzb):

(3) HR ( z ) = Δ T Δ t ( z ) = F net z t - F net z b ρ ( z ) c p z t - z b .

For atmospheric profiles, the vertical resolution from the BC profiles was used. Similarly, the heating rate profiles within the snowpack were calculated applying Eq. (3) by accounting for the snow density (set to 300 kg m−3) and the specific heat capacity of ice cp,snow=2060 J kg−1 K−1 at a temperature of 0 C. The layer thickness within TARTES and therefore the resolution of the heating rate profiles is 1 cm.

To separate the contribution of BC particles to the total heating rate, the effective BC heating rate HRBC(z) was calculated as the difference between the total heating rate HRtot(z) and the heating rate of the clean reference case HRclean(z):

(4) HR BC ( z ) = HR tot ( z ) - HR clean ( z ) .

If not indicated differently, radiative effects reported in this study refer to daily means accounting for the change of the SZA and the nighttime. Therefore, simulations were performed for a full diurnal cycle with a temporal resolution of 5 min. The simulated upward and downward irradiance were averaged. Then these daily mean irradiances were applied to calculate mean values of ΔFtot, ΔFatm, ΔFsnow, HRtot(z), and HRBC(z).

3 Results

3.1 Radiative impact of BC at surface level

3.1.1 Effect on surface albedo

The reduction of the snow surface albedo by BC impurities depends on the snow grain size. Here, changes of the snow surface albedo due to the combination of BC impurities and snow grain size variations were evaluated for Arctic conditions. The single-layer snowpack setup, as defined in Sect. 2, was used together with atmospheric properties representing the ACLOUD conditions. The SZA was set to a constant value of 60. Figure 4 shows the spectral snow albedo for variable BC particle mass concentrations (0, 5, and 20 ng g−1) as calculated with TARTES. The selected SSA values represent different snow types, as freshly fallen snow with small snow grains (SSA =60 m2 kg−1), aged snow which has undergone snow metamorphism (SSA =5 m2 kg−1) when surface temperature approaches 0 C, and moderately aged snow without melting (SSA =20 m2 kg−1), which was considered the default case. As expected, the highest values of surface albedo were obtained for the case with clean and fresh snow. Adding BC particles caused a decrease in the spectral surface albedo, in particular in the visible spectral range up to 700 nm, shown in the enlargement of Fig. 4. In contrast, the near-infrared spectral range was dominated by ice absorption, which is affected by the SSA (grain size).

Figure 4Spectral surface albedo of snow for cloudless conditions and a SZA of 60 for different SSA and BC particle mass concentrations. The inlay shows an enlargement of the spectral albedo between 350 and 700 nm.


From the simulations shown in Fig. 4 it becomes apparent that the decrease in surface albedo with increasing BC mass concentration is stronger for aged snow than for fresh snow. Fresh snow with smaller grains leads to an enhanced backscattering of the incident radiation, while larger grains allow for a deeper penetration of the incident radiation into the snowpack. Since the penetration depth for aged snow is deeper, the probability is higher that the radiation gets absorbed by the BC particles, leading to a decrease in the spectral surface albedo.

In the same way, the radiative forcing of BC particles embedded in the snow layer was calculated for overcast cloudy conditions (predefined low-level liquid water cloud case) to assess the relevance of changes of the BC mass concentration compared to variations in SSA and the illumination conditions. To estimate the relevance for the surface energy budget, the solar broadband forcing was analyzed by calculating the broadband albedo αbb. Therefore, the spectral albedo simulated by TARTES and the spectral downward irradiance Fλ(λ) simulated by libRadtran were used:

(5) α bb = α ( λ ) F λ ( λ ) d λ F λ ( λ ) d λ .

The calculated broadband surface albedo values are summarized in Table 4 for the cloudy and cloudless cases, respectively. For both cases, even the most extreme BC mass concentration reduced the surface albedo by less than 1 %. Contrarily, the snow grain size and the presence of clouds cause significant changes of the snow albedo. The difference of the broadband surface albedo between fresh and aged snow ranges up to 0.12 and 0.08 for cloudless and cloudy conditions, respectively, which is of the same order of magnitude as the effect of clouds (0.12 for fresh snow and 0.07 for aged snow). Therefore, for Arctic conditions, the impact of BC impurities on the broadband snow albedo is of minor importance, compared to the impact of modifying the snow grain size. Also, Warren and Wiscombe (1980) and Warren (2013) found only a small reduction of the broadband albedo between 0 % and 1 % for fresh snow and 0 %–3 % for aged snow when adding BC with a mass concentration of 34 ng g−1 to the clean snow.

Table 4Broadband surface albedo (αbb) of fresh (SSA =60 m2 kg−1) and aged snow (SSA =5 and 20 m2 kg−1) depending on the BC particle mass concentration and illumination condition.

Download Print Version | Download XLSX

3.1.2 Surface radiative forcing

The decrease in the snow surface albedo due to an increase in BC particle mass concentration or snow grain size directly alters the surface radiative forcing ΔFsnow. To quantify these radiative effects, ΔFsnow was first calculated for a fixed solar zenith angle of 60. A typical Arctic range of BC particle mass concentrations in snow and SSA values assuming the ACLOUD atmospheric conditions was applied. Figure 5 shows a contour plot of ΔFsnow for combinations of SSA and BC particle mass concentrations. For a BC particle mass concentration of 5 ng g−1 in snow representing clean conditions and a SSA larger than 20 m2kg−1, ΔFsnow ranges between 0.4 and 0.7 W m−2. Higher BC particle mass concentrations increase ΔFsnow depending on the snow grain size or SSA. The strongest increase in the solar radiative warming was calculated for small SSA values, corresponding to larger snow grain sizes. With the larger penetration depth for a smaller SSA, more radiation can be absorbed by the BC particles.

Figure 5Solar surface radiative forcing of BC impurities in snow ΔFsnow calculated for different SSA values and BC particle mass concentrations. The atmospheric conditions correspond to the ACLOUD case with a fixed SZA of 60. Horizontal red lines indicate typical Arctic conditions with rather clean and more polluted snow; the vertical line represents the most typical SSA.


To compare the radiative forcing at the surface of atmospheric BC particle profiles observed during the three aircraft campaigns ACLOUD, PAMARCMiP, and ARCTAS, the daily averaged surface radiative forcing was then analyzed. To limit the degree of freedom, the SSA was set to a default value of SSA = 20 m2 kg−1 representative for snow-covered Arctic sea ice. To estimate the relevance of the atmospheric BC particles, their separated radiative forcing ΔFatm was calculated. Additionally, the total radiative forcing ΔFtot combining the atmospheric and snow BC was analyzed. Figure 6 summarizes the daily averaged ΔFsnow (panel a), ΔFatm (panel b), and ΔFtot (panel c) for different BC particle mass concentrations in snow (0, 5, 20 ng g−1) on cloudless and cloudy conditions.

The BC particles embedded in snow lead to warming effects of up to 0.7 W m−2 for high BC mass concentrations of 20 ng g−1 and ACLOUD conditions. For ARCTAS ΔFsnow was slightly lower and for PAMARCMiP reduced by a factor of about 3. This difference is caused by the lower maximum Sun elevation during PAMARCMiP (location in higher latitude) resulting in a lower amount of available incoming solar irradiance compared to ACLOUD and ARCTAS (see range of SZA in Table 2).

Atmospheric BC particles reduce the incident solar radiation at the surface due to extinction, such that the atmospheric radiative forcing ΔFatm is negative in all scenarios (Figure 6b). This cooling at the surface is strongest with values up to 0.2 W m−2 in cloudless conditions for the ARCTAS case, where the largest atmospheric BC particle concentrations were observed. Despite having a BC optical depth of similar magnitude, the PAMARCMiP case (AODBC=0.006) shows a weaker radiative cooling compared to the ARCTAS case (AODBC=0.008) caused by the higher solar zenith angles in PAMARCMiP. Minor cooling of less than −0.02 W m−2 was observed for the ACLOUD case, where the atmosphere was rather clear with significant reduced atmospheric BC particle concentrations (factor of 10 lower than during ARCTAS). Comparing the simulations with different BC mass concentrations in snow showed only few effects of the surface properties on the radiative forcing of atmospheric BC. A slight decrease in ΔFatm with increasing BC mass concentrations was observed for the ARCTAS case, indicating that a lower surface albedo enhances the radiative forcing of atmospheric BC particles.

The cooling effect of atmospheric BC counteracts the warming effect of BC particles in snow and can lead to a positive and negative total radiative forcing. Figure 6c shows the total radiative forcing ΔFtot for all cases. For BC mass concentration of 20 ng g−1, all cases showed a total warming effect when the warming of BC in the snowpack exceeds the cooling by atmospheric BC. The strongest warming effect of up to 0.7 W m−2 was found for the ACLOUD case which is characterized by the pristine atmospheric conditions in the Arctic summertime. For less polluted snow (5 ng g−1), warming and cooling scenarios can occur depending on the concentration of atmospheric BC (ARCTAS shows a slight cooling) and the solar zenith angle (ACLOUD shows a significant warming effect). ΔFtot calculated for ACLOUD even exceeds the warming effect of PAMARCMiP for the higher BC mass concentration in the snow layer. This clearly demonstrates that the competition between the individual BC radiative forcings ΔFatm and ΔFsnow is strongly driven by solar zenith angle and the available solar radiation and is less affected by the BC concentrations itself.

Figure 6Daily mean of the solar surface radiative forcing simulated for the conditions of the three campaigns ACLOUD, PAMARCMiP, and ARCTAS assuming a fixed SSA of 20 m2 kg−1. The separated forcings by BC embedded in snow (ΔFsnow, panel a), atmospheric BC (ΔFatm, panel b), and the total forcing (ΔFtot, panel c) are shown. The daily mean solar radiative forcing of BC in cloudy conditions is displayed by shaded bars. The difference of ΔFtot between simulations with and without clouds is given in (d).


The available solar irradiance is strongly affected by the presence of clouds. Therefore, the impact of clouds on the BC radiative forcing was analyzed. Two cloud layers as defined in Sect. 2.2 were implemented in the simulations and considered in the calculation of ΔFtot/atm/snow (clean cloudy and polluted cloudy case in Eq. 2) to extract the pure BC radiative forcing. In Fig. 6 the BC radiative forcing of the cloudy scenarios are shown by the shaded bars. The magnitudes of ΔFsnow (panel a) and ΔFatm (panel b) are always reduced by the presence of clouds. ΔFsnow drops by about 15 % in all cases (0.1 W m−2 for ACLOUD and ARCTAS and high BC mass concentration in snow), while ΔFatm increases by more than 50 %, which amounts for ARCTAS to an absolute increase of 0.14 W m−2. Clouds reduce ΔFsnow because less radiation reaches the surface and can be absorbed by BC particles in the snowpack. The shift from a mostly direct illumination of the snow surface by the Sun to a diffuse illumination below the clouds is less significant as demonstrated in Table 4.

These different cloud effects counterbalance in the total radiative forcing ΔFatm (Fig. 6c). To illustrate the total effect by clouds, Fig. 6d shows the difference between cloudy and cloudless simulations. In all scenarios, still slight differences between cloudy and cloudless conditions were observed, but with different directions. For the ACLOUD case, the clouds reduce the warming effect of BC particles mainly due to a reduction of radiation that reaches the surface. As almost no atmospheric BC was present, only ΔFsnow is affected.

For the ARCTAS cases, the clouds always increased ΔFtot. For a BC mass concentration of 5 ng g−1 even the sign shifts from a total cooling to a total warming effect of BC. For ARCTAS, with high atmospheric BC concentrations, the presence of clouds mainly reduce the cooling effect of the atmospheric BC, ΔFatm. As the atmospheric BC layer was located mostly above the cloud, the radiative effect of the clouds, which is typically much stronger than the absorption by the atmospheric BC, reduces the significance of the atmospheric BC forcing. For higher BC mass concentrations in the snow, the increase in ΔFtot by adding a cloud becomes weaker because ΔFsnow simultaneously slightly decreases in cloudy conditions.

The PAMARCMiP case, characterized by the low Sun elevation, in general, showed a reduced effect by clouds. Here, the reduction of the cooling effect of atmospheric BC, ΔFatm, and the increase in the BC snow forcing, ΔFsnow, compete with each other and result in different total cloud radiative effects. Model runs with and without the upper ice cloud layer did not show any significant difference in ΔFtot/atm/snow, which allows the conclusion that mainly the presence of the low liquid water clouds affects the radiative forcing of BC particles.

In summary, the comparison of the radiative forcing by BC particles in snow and the atmosphere with typical concentrations and mass concentrations observed in Arctic spring and summer is rather small compared to other parameters (SZA, grain size) which contribute to solar cooling or heating at the surface level. The highest radiative cooling of BC particles was in the range of 1 W m−2 and is estimated for low SZA, high BC particle mass concentrations, and large grains.

3.2 Vertical radiative impact of BC particles in the atmosphere and snow

3.2.1 Heating rate profiles in the atmosphere

To quantify the absorption of solar radiation by Arctic atmospheric BC particles and consequent local warming effects, profiles of the heating rates were simulated for the three cases ACLOUD, ARCTAS, and PAMARCMiP. Based on simulations with and without atmospheric BC, the total heating rate HRtot(z) and the effective heating rate of BC particles HRBC(z) were calculated (see Eqs. 3 and 4). Figure 7 shows daily averaged profiles of HRtot(z) and HRBC(z) calculated for the three BC profiles. Solid lines represent the cloudless scenarios while dotted lines show simulations where the two predefined cloud layers were added. The location of the clouds is indicated by the gray shaded area. The highest total heating rates in cloudless conditions were found for the ACLOUD case, with maximum values of more than 1.2 K d−1 at about 2–4 km altitude. This altitude range was characterized by enhanced humidity leading to a stronger absorption of solar radiation by the water vapor. The spring campaigns ARCTAS and PAMARCMiP were characterized by lower water vapor concentrations (factor of 4 and 10 lower than for ACLOUD, respectively) and reduced incident solar radiation due to the time of year and latitude of the observations. This leads to significant lower values of HRtot(z) compared to the ACLOUD case. While ARCTAS showed a similar vertical pattern with maximum HRtot(z) of 0.5 K d−1 in the lower troposphere below 5 km altitude, the conditions during PAMARCMiP lead to a maximum of HRtot(z) of about 0.25 K d−1 located at 5–6 km altitude. This corresponds to the rather dry lower troposphere observed in springtime in the central Arctic. By adding clouds in the simulations, the highest HRtot(z) values were observed within the liquid water cloud layer, where solar radiation is absorbed by the cloud particles. Similar to the cloudless scenarios, the ACLOUD case showed the highest values of HRtot(z) with up to 4.1 K d−1 at cloud top of the lower liquid cloud layer. The absorption in the ice cloud is less pronounced, and the increase in HRtot(z) is significantly lower.

The profiles of the effective BC heating rate HRBC(z) (Fig. 7b) show a completely different pattern compared to HRtot(z). In general, HRBC(z) was about 1 order of magnitude lower than HRtot(z) for all three cases. Significant BC heating rates were observed only for the ARCTAS and PAMARCMiP cases with values up to 0.1 K d−1. The profiles of HRBC(z) are strongly correlated with the vertical distribution of BC particles in the measured profiles. Maximum HRBC(z) values were located in the pollution layers. The pollution layer observed during PAMARCMiP at 5 km and the BC layers of ARCTAS above 5 km altitude showed the largest relative impact of BC particles where nearly one-fifth and one-third, respectively, of the total solar heating is attributed to BC absorption. In lower altitudes of the ARCTAS case, the enhanced absorption by water vapor reduced the relative importance of BC particles. For the summer case of ACLOUD, HRBC(z) was rather small at all altitudes and contributed to the total radiative heating by only 10 %. However, in low altitudes, the absolute values of HRBC(z) were of the same order for both ACLOUD and the PAMARCMiP. This illustrates that the effect of a higher BC particle concentration during PAMARCMiP was compensated for by the dependence of HRBC(z) on the amount of the available incoming solar radiation and the atmospheric water vapor concentration.

Adding clouds in the simulations affected HRBC(z) of the three cases differently. While the clean atmosphere layer of ACLOUD and the PAMARCMiP cases shows almost no differences to cloudless conditions, a minor cloud effect was observed for the ARCTAS case and the polluted layer of PAMARCMiP. In both cases, the ice cloud led to a slight increase in HRBC(z) by about 5 % within and above the cloud layer. This was caused by the enhanced reflection of the incoming radiation which led to additional absorption of the reflected radiation by the atmospheric BC particles. At altitudes between the ice and liquid water clouds no significant effect by the clouds was observed. Within and below the liquid water cloud HRBC(z) was significantly reduced by almost 0.01 K d−1 for the ARCTAS case. This cloud effect was caused by the strong reflection of radiation at the cloud top leading to a reduction of radiation reaching into and below the cloud layer.

Comparing all simulations, it can be concluded that the absolute radiative effects of atmospheric BC particles are potentially strongest in early spring when incoming solar radiation starts to increase and BC particle concentration is still high enough. Furthermore, the surface conditions in spring were dominated by snow and ice coverage, which causes an increase in the amount of upward radiation contributing to the atmospheric heating rate. In late spring and summer, the BC particle concentration decreased rapidly, while the absorption by water vapor became more and more dominant with increasing temperatures.

Figure 7Daily averaged profiles of the total radiative heating rate HRtot(z) (a) and the effective BC heating rate HRBC(z) (b) calculated for the three cases ACLOUD, ARCTAS, and PAMARCMiP. Both cloudless simulations (solid lines) and cloudy scenarios (dotted lines) are shown. The gray shaded areas indicate the location of the cloud layers.


Figure 8(a) Transmissivity profiles of solar radiation within the snowpack for three single layers and one multilayer case assuming ACLOUD conditions. (b) Corresponding daily mean total heating rate profiles within a snowpack HRtot(z). (c) Corresponding effective BC heating rate profiles HRBC(z).


3.2.2 Heating rate profiles in the snowpack

Not only the snow albedo but also the transmission of radiation through the snowpack is affected by BC particles. For the atmospheric boundary conditions of the ACLOUD case, the radiative transfer in the snowpack was simulated and analyzed for different single-layer and multilayer scenarios as introduced in Table 3. The transmissivity was calculated from the ratio of the downward irradiance in the snow layer to the downward irradiance at the top of the snow layer. Figure 8a shows the transmissivity profiles of solar radiation within the snowpack. The homogeneous single-layer reference case without BC particles (SSA =20 m2 kg−1) illustrates the general decrease in transmissivity, which is reduced to 0.3 in 20 cm snow depth. Adding a typical Arctic BC concentration of 5 ng g−1 reduces the transmissivity to almost 0.2. This obviously may have an impact on the radiative processes below the snowpack, in and below the sea ice as discussed by, e.g., Tuzet et al. (2019) and Marks and King (2014). The inhomogeneous multilayer case shows in general lower transmissivities due to the enhanced reflection of the smaller snow grains at the top of the layer (SSA =60 m2 kg−1 down to 5 cm depth) but also indicates a significant dimming effect of the BC particles.

To assess, in which layers of the snowpack the strongest absorption of solar radiation and, therefore, a potential enhancement of the snow metamorphism is located, profiles of the heating rates within the snowpack HRtot(z) were calculated. To quantify how BC particles deposited in snow may change these heating profiles, the effective BC heating rates HRBC(z) were derived in a second step. Figure 8 shows HRtot(z) (panel a) and HRBC(z) (panel b) for all cases in the first 20 cm of the snowpack. For all cases, the total heating rate rapidly decreases by 1 order of magnitude within the first 10 cm of depth. The simulation for the single-layer (solid lines) snowpack shows the maximum values of HRtot(z) which were located in the topmost layers and reach values up to 6.6 K d−1 (note the scale break in Fig. 8a).

Assuming different BC mass concentrations in the single-layer case slightly increases HRtot(z) in the entire column. In the multilayer case, this increase is limited to the upper part of the profile. This contribution of BC particles to the total radiative heating was quantified by HRBC(z) and shown in Fig. 8b. The largest HRBC(z) values were observed for the most polluted single-layer case with a BC mass concentration of 20 ng g−1. For this case, the contribution by BC particles amounts to almost 0.9 K d−1 in the topmost layer, dropping down to a value of less than 0.1 K d−1 at 20 cm snow depth. Compared to the total radiative effect, HRBC(z) contributes about 15 % to the heating rate at the top snow layer and 40 % to the heating in the base layer. For the typical Arctic BC mass concentration of 5 ng g−1, this contribution of BC particles is significantly lower, ranging between 3 % and 20 %.

The multilayer cases are characterized by smaller snow grains in the top layer (SSA =40 m2 kg−1) compared to the single-layer cases (SSA =20 m2 kg−1) and, therefore, shows reduced values of HRtot(z). According to the structure of the snowpack, HRBC(z) is largest in the layer of the highest BC mass concentration. Beneath this layer (z<15 cm) the heating rates for the pristine and polluted cases are almost similar (HRBC(z)≈0 K d−1). In this base layer, the largest snow grains are assumed (lowest SSA), which increases the absorption of radiation by the snow ice water.

Based on these results, it becomes evident that the absorption of solar radiation by the ice water of the snow grains dominates the total heating rate in the snowpack, especially at the top layer, where most radiation is absorbed. Therefore, in Arctic conditions the snow grain size typically plays a larger role than the concentration of BC particles embedded in snow. To estimate if BC particles can accelerate the snow metamorphism, coupled snow physical models need to be applied (e.g., Tuzet et al.2017). However, compared to the results reported by Tuzet et al. (2017), who studied alpine snow with at least a magnitude higher BC mass concentrations, for Arctic conditions it is likely that the self-amplification of the snow metamorphism is dominated by the reduction of the surface albedo.

Simulations in cloudy conditions (not shown here) resulted in a reduced HRtot(z) and HRBC(z) because the clouds reduce the incoming solar radiation. Similarly, a change of the solar zenith angle affects the results by changing the available solar radiation. Therefore, the ACLOUD case used in the simulations presented in this section represents the maximum radiative effects compared to ARCTAS and PAMARCMiP conditions. In general, it can be concluded that the solar heating by BC particles embedded in the snowpack is most effective for low SZA (spring and summer conditions with high amount of available incoming radiation), decreasing SSA (aged snow in conditions near melting temperature), and increasing BC particle mass concentrations (accumulated BC particles caused by melting). Such conditions are mostly linked to late spring and summer, when the Sun is high, snow is close to melting, and BC has accumulated. This suggests that the maximum heating rates due to atmospheric BC and BC embedded in the snowpack typically occur in different periods of the year, early spring and early summer, respectively.

4 Summary and conclusions

This study analyzed the instantaneous solar radiative effect at the surface of Arctic BC particles (suspended in the atmosphere and embedded in the snowpack) over the sea-ice-covered Arctic Ocean. The difference of the BC effects in cloudless and cloudy conditions was compared. For this purpose, an atmospheric and a snow radiative transfer model were iteratively coupled to account for the radiative interactions between both compartments (atmosphere and snow layer). Typical atmospheric BC vertical profiles and BC particle mass concentrations in the snowpack, derived from three field campaigns in the North American and the North Atlantic Arctic, ACLOUD, ARCTAS, and PAMARCMiP, were used in the simulations. These locations typically are not affected by local pollution, but by long-range transport of BC particles. The BC radiative effects were quantified by the surface radiative forcing and profiles of heating rates in the atmosphere and the snow, which were presented on the basis of daily averages. For the surface radiative forcing, the contributions by atmospheric and snow BC particles were separated. For the heating rate profiles, the effective contribution of BC particles to the total heating rates was derived and compared to further atmospheric and snow parameters also leading to a warming or cooling (e.g., water vapor, clouds, snow grain size).

The magnitude of the atmospheric BC radiative forcing at the surface derived in this study (up to 0.2 W m−2) agrees quite well with findings from Wendling et al. (1985). They reported a BC-induced solar cooling in the range of 0.0 to 0.5 W m−2 for spring measurements in the Svalbard area. Further, the solar surface radiative forcing due to BC embedded in snow showed solar warming between 0.05 and 0.7 W m−2 depending on the BC mass concentration and incident solar irradiance. For comparison, Dou and Cun-De (2016) deduced an averaged solar warming over Svalbard in spring of 0.54 W m−2 based on a BC mass concentration of 5 ng g−1 in snow.

The simulations suggest that for the specific Arctic cases investigated in our study, the radiative forcing of BC is small compared to the radiative impact of other parameters (water vapor, clouds, snow grain size). The significance of the BC radiative effects shows a strong seasonal dependence. In cloudless conditions, the absorption by atmospheric water vapor shows a much stronger contribution to the atmospheric heating rates than the radiative effect of BC particles. In summer (ACLOUD) and in lower latitudes (ARCTAS), the Arctic shows the most humid conditions, where absorption of water vapor dominates over the BC radiative effects. Similarly, the available incident solar radiation limits the magnitude of the BC radiative effects. Despite the more polluted atmosphere, the low solar zenith angle of the cases of PAMARCMiP (high latitude) and ARCTAS (early spring season) did show lower BC radiative effects than the ACLOUD case. Thus, over the sea-ice-covered Arctic Ocean, the BC radiative effect is about an order of magnitude lower than observed in lower and tropical latitudes, where the pollution level is also typically higher. For example, studies investigating strong pollution conditions in northern India or China reported BC heating rates in the atmosphere larger than 2 K d−1, which may significantly influence the lapse rate and the atmospheric stability (Tripathi et al.2007; Wendisch et al.2008). For the rather pristine Arctic, this study showed significantly lower daily mean BC heating rates of maximum 0.1 K d−1, which do not have the potential to significantly modify the atmospheric stability. However, in other Arctic regions characterized by higher atmospheric BC particle concentrations due to local fires, e.g., northern Siberia, a stronger impact can be expected.

Similarly, the mass concentration of BC particles embedded in the Arctic snowpack is far lower than observed in alpine snow in lower latitudes. Accordingly, the absorption of radiation by the snow water itself dominates the radiative warming in the snowpack. For typical conditions of the central Arctic, the absorption due to BC particles contributes only 3 % to the total heating rate in the uppermost snow layer. These results indicate that the microphysical properties of the snowpack (mainly snow grain sizes) are more important drivers for the degree and strength of the snow metamorphism. It needs to be considered that this picture might change if the accumulation of BC particles is more efficient than it is over the snow-covered Arctic sea ice, where the sea ice and snowpack do not last more than 1 to 3 years. Accumulation of BC on the Greenlandic glaciers will amplify the radiative forcing on a local scale. Furthermore, BC particles are not the only light-absorbing impurities which are transported into the Arctic. The relevance of dust particles and microorganisms is currently the subject of the scientific discussion and may exceed the effect of BC particles (Kylling et al.2018; Skiles et al.2018).

However, the changing relevance of the BC radiative effects suggests that the maximum heating rates due to atmospheric BC and BC embedded in the snowpack typically occur in different periods of the year. While atmospheric BC particles reveal the largest radiative effects in early spring (high concentration of atmospheric BC, medium–high Sun, low water vapor), the BC particles embedded in snow warm more effectively in early summer (accumulation of BC particles in snow, high Sun, large snow grain size). To estimate the role of clouds in the surface warming or cooling by BC particles and the BC heating rates, radiative transfer simulations assuming cloudless and cloudy conditions were compared. Clouds reflect the incident solar radiation and, therefore, reduce the available radiation reaching the surface. This reduces the potential of the warming effect by BC particles embedded in the snow. Similarly, the cooling effect by atmospheric BC on the surface radiative budget is weakened in the presence of clouds. The competition of these two cloud effects depends on the BC concentrations in the snow and atmosphere and is affected by the increased broadband surface albedo and the multiple scattering in the presence of a cloud layer. The profiles of the effective BC heating rates are mainly affected by the ice cloud at higher altitude. Within and above the cloud, the radiation reflected by the cloud enhances the local radiative heating by BC. Contrarily, a low liquid water cloud reduces the available incoming radiation, such that the effective BC radiative effect is lower for the cloudy case compared to the cloudless case. For the same reason, the presence of clouds reduces the radiative heating rates within the snowpack.

For the sea-ice-covered Arctic Ocean, we conclude that (i) the warming effect of BC embedded in the snow overcompensates for the atmospheric BC cooling effect at the surface, (ii) the impact of clouds reduces both the surface cooling by atmospheric BC particles and the warming by BC particles embedded in snow, and (iii) the BC radiative effect is of minor importance compared to other absorbers. However, for the expected increase in BC particle mass concentrations in the future, the relative importance of BC particles might need to be re-evaluated. Additionally, ongoing research, e.g., triggered by the current MOSAiC (Multidisciplinary drifting Observatory for the Study of Arctic Climate) experiment, will enable the quantification of the radiative effects of BC in the eastern and central Arctic also using the methods proposed here.

Data availability

Atmospheric BC mass concentrations for the ARCTAS campaign are available on, last access: 22 January 2020 (Kondo, 2009). PAMARCMiP and ACLOUD data are available on PANGAEA (, Herber et al.2019, and, Zanatta and Herber2019, respectively).

Author contributions

All authors contributed to the editing of the manuscript and to the discussion of the results. MW, AE, and AH designed this study. TD drafted the manuscript, performed the radiative transfer simulations, and prepared the figures. AE, EJ, and BH contributed to the interpretation of the radiative transfer simulations. MZ processed the SP2 data. JS compiled the atmospheric BC profiles.

Competing interests

The authors declare that they have no conflict of interest.

Special issue statement

This article is part of the special issue “Arctic mixed-phase clouds as studied during the ACLOUD/PASCAL campaigns in the framework of (AC)3 (ACP/AMT/ESSD inter-journal SI)”. It is not associated with a conference.


We gratefully acknowledge the funding by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – project ID 268020496 – TRR 172, within the Transregional Collaborative Research Center “ArctiC Amplification: Climate Relevant Atmospheric and SurfaCe Processes, and Feedback Mechanisms (AC)3”. We thank Sho Ohata and Makoto Koike for providing the atmospheric BC profiles measured during the PAMARCMiP campaign 2018.

Financial support

This research has been supported by the Deutsche Forschungsgemeinschaft (grant no. 268020496 – TRR 172).

Review statement

This paper was edited by Jost Heintzenberg and reviewed by Quentin Libois and one anonymous referee.


Anderson, G., Clough, S., Kneizys, F., Chetwynd, J., and Shettle, E.: AFGL Atmospheric Constituent Profiles (0–120 km), Tech. Rep. AFGL-TR-86-0110, AFGL (OPI), Hanscom AFB, MA 01736, 1986. a

Arnold, S., Law, K., Thomas, J., Starckweather, S., von Salzen, K., Stohl, A., Sharma, S., Lund, M., Flanner, M., Petäjä, T., Tanimoto, H., Gamble, J., Dibb, J., Melamad, M., Johnson, N., Fidel, M., Tynkkynen, V.-P., Baklanov, A., Eckhardt, S., Monks, S., Browse, J., and Bozem, H.: Arctic air pollution: Challenges and opportunities for the next decade, ELEMENTA – Science of the Antropocene, 4, 1–17,, 2016. a

Bierwirth, E., Ehrlich, A., Wendisch, M., Gayet, J.-F., Gourbeyre, C., Dupuy, R., Herber, A., Neuber, R., and Lampert, A.: Optical thickness and effective radius of Arctic boundary-layer clouds retrieved from airborne nadir and imaging spectrometry, Atmos. Meas. Tech., 6, 1189–1200,, 2013. a

Blanchard, Y., Royer, A., O'Neill, N. T., Turner, D. D., and Eloranta, E. W.: Thin ice clouds in the Arctic: cloud optical depth and particle size retrieved from ground-based thermal infrared radiometry, Atmos. Meas. Tech., 10, 2129–2147,, 2017. a

Bond, T. C., Doherty, S. J., Fahey, D. W., Forster, P. M., Berntsen, T., DeAngelo, B. J., Flanner, M. G., Ghan, S., Karcher, B., Koch, D., Kinne, S., Kondo, Y., Quinn, P. K., Sarofim, M. C., Schultz, M. G., Schulz, M., Venkataraman, C., Zhang, H., Zhang, S., Bellouin, N., Guttikunda, S. K., Hopke, P. K., Jacobson, M. Z., Kaiser, J. W., Klimont, Z., Lohmann, U., Schwarz, J. P., Shindell, D., Storelvmo, T., Warren, S. G., and Zender, C. S.: Bounding the role of black carbon in the climate system: A scientific assessment, J. Geophys. Res., 118, 5380–5552,, 2013. a, b, c

Choudhury, B. J. and Chang, A. T. C.: The albedo of snow for partially cloudy skies, Bound.-Lay. Meteorol., 20, 371–389,, 1981. a

Corbett, J. J., Lack, D. A., Winebrake, J. J., Harder, S., Silberman, J. A., and Gold, M.: Arctic shipping emissions inventories and future scenarios, Atmos. Chem. Phys., 10, 9689–9704,, 2010. a

Doherty, S. J., Warren, S. G., Grenfell, T. C., Clarke, A. D., and Brandt, R. E.: Light-absorbing impurities in Arctic snow, Atmos. Chem. Phys., 10, 11647–11680,, 2010. a, b, c, d, e, f, g, h, i, j

Dou, T. and Cun-De, X.: An overview of black carbon deposition and its radiative forcing over the Arctic, Adv. Clim. Change Res., 7, 115–122,, 2016. a, b

Dou, T., Xiao, C., Du, Z., Schauer, J., Ren, H., Ge, B., Xie, A., Tan, J., Fu, P., and Zhang, Y.: Sources, evolution and impacts of EC and OC in snow on sea ice: a measurement study in Barrow, Alaska, Sci. Bull., 62, 1547–1554,, 2017. a

Ehrlich, A., Wendisch, M., Lüpkes, C., Buschmann, M., Bozem, H., Chechin, D., Clemen, H.-C., Dupuy, R., Eppers, O., Hartmann, J., Herber, A., Jäkel, E., Järvinen, E., Jourdan, O., Kästner, U., Kliesch, L.-L., Köllner, F., Mech, M., Mertes, S., Neuber, R., Ruiz-Donoso, E., Schnaiter, M., Schneider, J., Stapf, J., and Zanatta, M.: A comprehensive in situ and remote sensing data set from the Arctic CLoud Observations Using airborne measurements during polar Day (ACLOUD) campaign, Earth Syst. Sci. Data, 11, 1853–1881,, 2019. a, b

Emde, C., Buras-Schnell, R., Kylling, A., Mayer, B., Gasteiger, J., Hamann, U., Kylling, J., Richter, B., Pause, C., Dowling, T., and Bugliaro, L.: The libRadtran software package for radiative transfer calculations (version 2.0.1), Geosci. Model Dev., 9, 1647–1672,, 2016. a

Flanner, M. G., Zender, C. S., Randerson, J. T., and Rasch, P. J.: Present-day climate forcing and response from black carbon in snow, J. Geophys. Res, 112, D11202,, 2007. a

Flanner, M. G.: Arctic climate sensitivity to local black carbon, J. Geophys. Res.-Atmos., 118, 1840–1851,, 2013. a

Forsström, S., Isaksson, E., Skeie, R., Ström, J., Pedersen, C., Hudson, S., Berntsen, T., Lihavainen, H., Godtliebsen, F., and Gerland, S.: Elemental carbon measurements in European Arctic snow packs, J. Geophys. Res.-Atmos., 118, 13614–13627,, 2013. a, b, c, d, e

Fu, Q.: A New Parameterization of an Asymmetry Factor of Cirrus Clouds for Climate Models, J. Atmos. Sci., 64, 4140–4150,, 2007. a

Gardner, A. S. and Sharp, M. J.: A review of snow and ice albedo and the development of a new physically based broadband albedo parameterization, J. Geophys. Res.-Earth, 115, F01009,, 2010. a

Gilgen, A., Huang, W. T. K., Ickes, L., Neubauer, D., and Lohmann, U.: How important are future marine and shipping aerosol emissions in a warming Arctic summer and autumn?, Atmos. Chem. Phys., 18, 10521–10555,, 2018. a

Gueymard, C. A.: The sun's total and spectral irradiance for solar energy applications and solar radiation models, Sol. Energ., 76, 423–453,, 2004. a

Herber, A.: Meteorological observations during POLAR 5 campaign PAMARCMIP 2018, PANGAEA,, 2019. a

Herber, A., Bozem, H., Hendricks, S., Holzinger, R., Jäkel, E., Koike, M., Neuber, R., Petzold, A., and Stratmann, F.: Raw data of POLAR 5 campaign PAMARCMIP 2018, PANGAEA, (last access: 28 May 2020), 2019. a

Herber, A. B., Haas, C., Stone, R. S., Bottenheim, J. W., Liu, P., Li, S.‐M., Staebler, R. M., Strapp, J. W., and Dethloff, K.: Regular airborne surveys of Arctic sea ice and atmosphere, EOS Transactions American Geophysical Union, Eos Trans. AGU, 93, 41–48,, 2012. a

Hess, M., Koepke, P., and Schult, I.: Optical properties of aerosols and clouds: the software package OPAC, Bull. Am. Meteorol. Soc., 79, 831–844,<0831:OPOAAC>2.0.CO;2, 1998. a

Jacob, D. J., Crawford, J. H., Maring, H., Clarke, A. D., Dibb, J. E., Emmons, L. K., Ferrare, R. A., Hostetler, C. A., Russell, P. B., Singh, H. B., Thompson, A. M., Shaw, G. E., McCauley, E., Pederson, J. R., and Fisher, J. A.: The Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) mission: design, execution, and first results, Atmos. Chem. Phys., 10, 5191–5212,, 2010. a, b

Joseph, J., Wiscombe, W., and Weinman, J.: The Delta-Eddington Approximation for Radiative Flux Transfer, J. Atmos. Sci., 33, 2452–2459,<2452:TDEAFR>2.0.CO;2, 1977. a

Kokhanovsky, A. A. and Zege, E. P.: Scattering optics of snow, Appl. Opt., 43, 1589–1602,, 2004. a

Kondo, Y.: SP2_DC8 data, available at: (last access 22 January 2020), 2009. 

Kylling, A., Groot Zwaaftink, C. D., and Stohl, A.: Mineral dust instantaneous radiative forcing in the Arctic, Geophys. Res. Lett., 45, 4290–4298,, 2018. a

Law, K. S., Stohl, A., Quinn, P. K., Brock, C. A., Burkhart, J. F., Paris, J.-D., Ancellet, G., Singh, H. B., Roiger, A., Schlager, H., Dibb, J., Jacob, D. J., Arnold, S. R., Pelon, J., and Thomas, J. L.: Arctic Air Pollution: New Insights from POLARCAT-IPY, Bull. Am. Meteorol. Soc., 95, 1873–1895,, 2014. a

Leaitch, W. R., Korolev, A., Aliabadi, A. A., Burkart, J., Willis, M. D., Abbatt, J. P. D., Bozem, H., Hoor, P., Köllner, F., Schneider, J., Herber, A., Konrad, C., and Brauner, R.: Effects of 20–100 nm particles on liquid clouds in the clean summertime Arctic, Atmos. Chem. Phys., 16, 11107–11124,, 2016. a

Libois, Q., Picard, G., France, J. L., Arnaud, L., Dumont, M., Carmagnola, C. M., and King, M. D.: Influence of grain shape on light penetration in snow, The Cryosphere, 7, 1803–1818,, 2013. a, b

Libois, Q., Picard, G., Dumont, M., Arnaud, L., Sergent, C., Pougatch, E., Sudul, M., and Vial, D.: Experimental determination of the absorption enhancement parameter of snow, J. Glaciol., 60, 714–724,, 2014. a

Liu, J., Fan, S., Horowitz, L. W., and Levy II, H.: Evaluation of factors controlling long-range transport of black carbon to the Arctic, J. Geophys. Res.-Atmos., 116, D04307,, 2011. a

Luebke, A. E., Avallone, L. M., Schiller, C., Meyer, J., Rolf, C., and Krämer, M.: Ice water content of Arctic, midlatitude, and tropical cirrus – Part 2: Extension of the database and new statistical analysis, Atmos. Chem. Phys., 13, 6447–6459,, 2013. a

Marks, A. A. and King, M. D.: The effect of snow/sea ice type on the response of albedo and light penetration depth (e-folding depth) to increasing black carbon, The Cryosphere, 8, 1625–1638,, 2014. a

Matsui, H., Kondo, Y., Moteki, N., Takegawa, N., Sahu, L. K., Zhao, Y., Fuelberg, H. E., Sessions, W. R., Diskin, G., Blake, D. R., Wisthaler, A., and Koike, M.: Seasonal variation of the transport of black carbon aerosol from the Asian continent to the Arctic during the ARCTAS aircraft campaign RID F-4366-2011, J. Geophys. Res., 116, D05202,, 2011. a

Mayer, B. and Kylling, A.: Technical note: The libRadtran software package for radiative transfer calculations – description and examples of use, Atmos. Chem. Phys., 5, 1855–1877,, 2005. a

Moteki, N. and Kondo, Y.: Effects of Mixing State on Black Carbon Measurements by Laser-Induced Incandescence, Aerosol Sci. Technol., 41, 398–417,, 2007. a

Nakoudi, K., Ritter, C., Böckmann, C., Kunkel, D., Eppers, O., Rozanov, V., Mei, L., Pefanis, V., Jäkel, E., Herber, A., Maturilli, M., and Neuber, R.: Does the Intra-Arctic Modification of Long-Range Transported Aerosol Affect the Local Radiative Budget? (A Case Study), Remote Sens., 12, 2112,, 2020. 

Namazi, M., von Salzen, K., and Cole, J. N. S.: Simulation of black carbon in snow and its climate impact in the Canadian Global Climate Model, Atmos. Chem. Phys., 15, 10887–10904,, 2015. a

Pedersen, C., Gallet, J.-C., Ström, J., Gerland, S., Hudson, S., Forsström, S., Isaksson, E., and Berntsen, T.: In-situ observations of black carbon in snow and the corresponding spectral surface albedo reduction, J. Geophys. Res.-Atmos., 120, 1476–1489,, 2015. a, b, c, d, e, f, g, h, i

Petzold, A., Ogren, J. A., Fiebig, M., Laj, P., Li, S.-M., Baltensperger, U., Holzer-Popp, T., Kinne, S., Pappalardo, G., Sugimoto, N., Wehrli, C., Wiedensohler, A., and Zhang, X.-Y.: Recommendations for reporting “black carbon” measurements, Atmos. Chem. Phys., 13, 8365–8379,, 2013. a, b

Samset, B. H., Myhre, G., Schulz, M., Balkanski, Y., Bauer, S., Berntsen, T. K., Bian, H., Bellouin, N., Diehl, T., Easter, R. C., Ghan, S. J., Iversen, T., Kinne, S., Kirkevåg, A., Lamarque, J.-F., Lin, G., Liu, X., Penner, J. E., Seland, Ø., Skeie, R. B., Stier, P., Takemura, T., Tsigaridis, K., and Zhang, K.: Black carbon vertical profiles strongly affect its radiative forcing uncertainty, Atmos. Chem. Phys., 13, 2423–2434,, 2013. a

Samset, B. H., Myhre, G., Herber, A., Kondo, Y., Li, S.-M., Moteki, N., Koike, M., Oshima, N., Schwarz, J. P., Balkanski, Y., Bauer, S. E., Bellouin, N., Berntsen, T. K., Bian, H., Chin, M., Diehl, T., Easter, R. C., Ghan, S. J., Iversen, T., Kirkevåg, A., Lamarque, J.-F., Lin, G., Liu, X., Penner, J. E., Schulz, M., Seland, Ø., Skeie, R. B., Stier, P., Takemura, T., Tsigaridis, K., and Zhang, K.: Modelled black carbon radiative forcing and atmospheric lifetime in AeroCom Phase II constrained by aircraft observations, Atmos. Chem. Phys., 14, 12465–12477,, 2014. a, b

Sand, M., Berntsen, T., Seland, Ø., and Kristjánsson, J.: Arctic surface temperature change to emissions of black carbon within Arctic or midlatitudes, J. Geophys. Res.-Atmos., 118, 7788–7798,, 2013. a

Sand, M., Samset, B. H., Balkanski, Y., Bauer, S., Bellouin, N., Berntsen, T. K., Bian, H., Chin, M., Diehl, T., Easter, R., Ghan, S. J., Iversen, T., Kirkevåg, A., Lamarque, J.-F., Lin, G., Liu, X., Luo, G., Myhre, G., Noije, T. V., Penner, J. E., Schulz, M., Seland, Ø., Skeie, R. B., Stier, P., Takemura, T., Tsigaridis, K., Yu, F., Zhang, K., and Zhang, H.: Aerosols at the poles: an AeroCom Phase II multi-model evaluation, Atmos. Chem. Phys., 17, 12197–12218,, 2017. a

Schacht, J., Heinold, B., Quaas, J., Backman, J., Cherian, R., Ehrlich, A., Herber, A., Huang, W. T. K., Kondo, Y., Massling, A., Sinha, P. R., Weinzierl, B., Zanatta, M., and Tegen, I.: The importance of the representation of air pollution emissions for the modeled distribution and radiative effects of black carbon in the Arctic, Atmos. Chem. Phys., 19, 11159–11183,, 2019. a

Schulz, H., Zanatta, M., Bozem, H., Leaitch, W. R., Herber, A. B., Burkart, J., Willis, M. D., Kunkel, D., Hoor, P. M., Abbatt, J. P. D., and Gerdes, R.: High Arctic aircraft measurements characterising black carbon vertical variability in spring and summer, Atmos. Chem. Phys., 19, 2361–2384,, 2019. a

Sharma, S., Ishizawa, M., Chan, D., Lavoué, D., Andrews, E., Eleftheriadis, K., and Maksyutov, S.: 16-year simulation of Arctic black carbon: Transport, source contribution, and sensitivity analysis on deposition, J. Geophys. Res.-Atmos., 118, 943–964,, 2013. a

Skiles, S. M., Flanner, M., Cook, J., Dumont, M., and Painter, T. H.: Radiative forcing by light-absorbing particles in snow, Nat. Clim. Change, 8, 964–971,, 2018. a

Stamnes, K., Tsay, S.-C., Wiscombe, W., and Laszlo, I.: DISORT, A General-Purpose Fortran Program for Discrete-Ordinate-Method Radiative Transfer in Scattering and Emitting Layered Media: Documentation of Methodology, Tech. Rep., Dept. of Physics and Engineering Physics, Stevens Institute of Technology, Hoboken, NJ 07030, 2000. a

Stapf, J., Ehrlich, A., Jäkel, E., Lüpkes, C., and Wendisch, M.: Reassessment of the common concept to derive the surface cloud radiative forcing in the Arctic: Consideration of surface albedo – cloud interactions, Atmos. Chem. Phys. Discuss.,, in review, 2019. a

Stohl, A., Klimont, Z., Eckhardt, S., Kupiainen, K., Shevchenko, V. P., Kopeikin, V. M., and Novigatsky, A. N.: Black carbon in the Arctic: the underestimated role of gas flaring and residential combustion emissions, Atmos. Chem. Phys., 13, 8833–8855,, 2013. a

Stone, R., Herber, A., Vitale, V., Mazzola, M., Lupi, A., Schnell, R., Dutton, E., Liu, P., Li, S.-M., Dethloff, K., Lampert, A., Ritter, C., Stock, M., Neuber, R., and Maturilli, M.: A three-dimensional characterization of Arctic aerosol particles from airborne Sun photometer observations: PANARCMIP, April 2009, J. Geophys. Res., 115, D13203,, 2010. a

Tripathi, S., Srivastava, A., Dey, S., Satheesh, S., and Moorthy, K.: The vertical profile of atmospheric heating rate of black carbon aerosols at Kanpur in northern India, Atmos. Environ., 41, 6909–6915,, 2007. a

Tuzet, F., Dumont, M., Lafaysse, M., Picard, G., Arnaud, L., Voisin, D., Lejeune, Y., Charrois, L., Nabat, P., and Morin, S.: A multilayer physically based snowpack model simulating direct and indirect radiative impacts of light-absorbing impurities in snow, The Cryosphere, 11, 2633–2653,, 2017. a, b

Tuzet, F., Dumont, M., Arnaud, L., Voisin, D., Lamare, M., Larue, F., Revuelto, J., and Picard, G.: Influence of light-absorbing particles on snow spectral irradiance profiles, The Cryosphere, 13, 2169–2187,, 2019. a

Warren, S. and Wiscombe, W.: A Model for the Spectral Albedo of Snow, II: Snow containing atmospheric aerosols, J. Atmos. Sci., 37, 2734–2745,<2734:AMFTSA>2.0.CO;2, 1980. a, b

Warren, S. G.: Optical properties of snow, Rev. Geophys., 20, 67–89,, 1982. a

Warren, S. G.: Can black carbon in snow be detected by remote sensing?, J. Geophys. Res., 118, 779–786,, 2013. a, b

Wendisch, M., Hellmuth, O., Ansmann, A., J. Heintzenberg, J., Engelmann, R., Althausen, D., Eichler, H., Müller, D., Hu, M., Zhang, Y., and Mao, J.: Radiative and dynamic effects of absorbing aerosol particles over the Pearl River Delta, China, Atmos. Environ., 42, 6405–6416, 2008. a, b

Wendisch, M., Brückner, M., Burrows, J. P., Crewell, S., Dethloff, K., Ebell, K., Lüpkes, C., Macke, A., Notholt, J., Quaas, J., Rinke, A., and Tegen, I.: Understanding causes and effects of rapid warming in the Arctic, EOS, 98, 22–26,, 2017. a

Wendisch, M., Macke, A., Ehrlich, A., Lüpkes, C., Mech, M., Chechin, D., Dethloff, K., Velasco, C. B., Bozem, H., Brückner, M., Clemen, H.-C., Crewell, S., Donth, T., Dupuy, R., Ebell, K., Egerer, U., Engelmann, R., Engler, C., Eppers, O., Gehrmann, M., Gong, X., Gottschalk, M., Gourbeyre, C., Griesche, H., Hartmann, J., Hartmann, M., Heinold, B., Herber, A., Herrmann, H., Heygster, G., Hoor, P., Jafariserajehlou, S., Jäkel, E., Järvinen, E., Jourdan, O., Kästner, U., Kecorius, S., Knudsen, E. M., Köllner, F., Kretzschmar, J., Lelli, L., Leroy, D., Maturilli, M., Mei, L., Mertes, S., Mioche, G., Neuber, R., Nicolaus, M., Nomokonova, T., Notholt, J., Palm, M., van Pinxteren, M., Quaas, J., Richter, P., Ruiz-Donoso, E., Schäfer, M., Schmieder, K., Schnaiter, M., Schneider, J., Schwarzenböck, A., Seifert, P., Shupe, M. D., Siebert, H., Spreen, G., Stapf, J., Stratmann, F., Vogl, T., Welti, A., Wex, H., Wiedensohler, A., Zanatta, M., and Zeppenfeld, S.: The Arctic Cloud Puzzle: Using ACLOUD/PASCAL Multiplatform Observations to Unravel the Role of Clouds and Aerosol Particles in Arctic Amplification, B. Am. Meteorol. Soc., 100, 841–871,, 2019.  a, b

Wendling, P., Wendling, R., Renger, W., Covert, D., Heintzenberg, J., and Moerl, P.: Calculated radiative effects of Arctic haze during a pollution episode in Spring 1983 based on ground-based and airborne measurements, Atmos. Environ., 19, 2181–2193, 1985. a, b

Wiscombe, W. J. and Warren, S. G.: A Model for the Spectral Albedo of Snow, I: Pure Snow, J. Atmos. Sci., 37, 2712–2733,<2712:AMFTSA>2.0.CO;2, 1980. a

Wyser, K.: The Effective Radius in Ice Clouds, J. Clim., 11, 1793–1802,<1793:TERIIC>2.0.CO;2, 1998. a

Zanatta, M. and Herber, A.: Aircraft measurements of refractory black carbon in the Arctic during the ACLOUD campaign 2017, Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research, Bremerhaven, PANGAEA, (last access: 28 May 2020), 2019. a

Short summary
Solar radiative effects of Arctic black carbon (BC) particles (suspended in the atmosphere and in the surface snowpack) were quantified under cloudless and cloudy conditions. An atmospheric and a snow radiative transfer model were coupled to account for radiative interactions between both compartments. It was found that (i) the warming effect of BC in the snowpack overcompensates for the atmospheric BC cooling effect, and (ii) clouds tend to reduce the atmospheric BC cooling and snow BC warming.
Final-revised paper