1.

## Introduction

The effect of climatic and cryospheric change on carbon cycling is of great research interest due to carbon potential feedback on temperature and ecosystems (Hood et al., 2009; Singer et al., 2012; Lawson et al., 2014; Hood et al., 2015; Musilova et al., 2017). The reservoir of dissolved organic carbon (DOC) in snow and ice is small (4.48 ± 2.79 PgC) in comparison to that stored in gross primary productivity and respiration. However, it is important because DOC in snow and ice is much more bioavailable compared to the net changes in terrestrial carbon storage (Ping et al., 2008; Anesio et al., 2009; Legrand et al., 2013; Fellman et al., 2015; Hood et al., 2015; Liu et al., 2016; Yan et al., 2016), which may impact the carbon cycle of downstream ecosystems (Hood et al., 2009; Spencer et al., 2014). Snow and ice constitute an important link between the terrestrial and aquatic carbon flux due to snow meltwater transportation and is of increasing importance in land-to-ocean fluxes of organic carbon in cryospheric regions (Antony et al., 2014; Lawson et al., 2014; D’Amore et al., 2015; Fellman et al., 2015; Hood et al., 2015; Huntington et al., 2016; Gao et al., 2019, 2020).

Snow cover, a critically important component the Arctic, is greatly affected by climate change (Callaghan et al., 2011; Bokhorst et al., 2016). The duration of the presence of the Arctic snow cover is decreasing rapidly (about 3–5 days per decade) as a result of a delay in the onset of snow cover and the earlier onset of the spring melt (Derksen et al., 2015). This, in turn, affects Arctic sea ice, ecosystems, and societies (Goldenson et al., 2012; Cooper, 2014). In addition to a warmer Arctic climate (Callaghan et al., 2011; Pithan & Mauritsen, 2014; Huang et al., 2017), the presence of light-absorbing organic carbon and black carbon in snow plays a pivotal role in accelerating snow aging and the speed of the snowmelt (Warren & Wiscombe, 1980; Flanner et al., 2007; Doherty et al., 2010; Dumont et al., 2014).

Many studies have shown that DOC concentrations in snow of the polar regions varied considerably across the snowmelt season, both temporally and spatially (Antony et al., 2011; Legrand et al., 2013; Antony et al., 2014; Fellman et al., 2015), with a range of 0.01 to 43.2 mg L−1 (Hood et al., 2015). The highest DOC concentrations were found in basal ice and cryoconite holes, reflecting localized accumulation of organic matter from overrun organic carbon (OC) (Stibal et al., 2012) and in-situ production (Stibal et al., 2010; Telling et al., 2012), respectively. This variability may be the result of the different sources of glacier-derived DOC as mentioned in previous studies (Kim et al., 2005; Hood et al., 2009; Singer et al., 2012; Stubbins et al., 2012; Spencer et al., 2014, 2016), including deposition of aerosols containing soil organic matter and fossil fuel combustion byproducts.

A fraction of DOC known as water-soluble brown carbon is considered to be light-absorber, akin to coloured dissolved organic matter (CDOM) in the atmosphere. Water-soluble CDOM exhibits strong light-absorbing properties in the ultraviolet wavelength (Andreae & Gelencser, 2006; Chen & Bond, 2010; Yan et al., 2016). The radiative forcing (RF) attributed to DOC relative to black carbon (BC) in aerosols has been estimated at 2–10% in a typical polluted area in North China (Kirillova et al., 2014), and approximately 1% in a remote island in the Indian Ocean (Bosch et al., 2014). Using a model simulation, Lin et al. (2014) reported that the RF caused by organic carbon deposited on terrestrial snow and sea ice ranged from +0.0011 to +0.0031 W m−2.

Previous studies have shown that the organic carbon in snow contributes to a decreased albedo and has a high bioavailability (Lin et al., 2014; Fellman et al., 2015). However, data published is still sparse. Here, we focus on documenting the properties of DOC in snow cover at Barrow and across the Alaska, which will enhance our understanding on its characteristics, bioavailability, light-absorbing properties, and potential impact on snowmelt. This study will provide insights into the DOC in the Arctic snow.

2.

## Methodology

2.1.

### Study area and sampling

From April 6, 2017 to May 18, 2017, surface snow samples were collected from snow cover on sea ice in the Chukchi Sea, about 12 km northeast of Barrow, Alaska (Fig. S1), and stored in pre-cleaned polycarbonate bottles. Vertical profile samples throughout the snowpack depth were collected at a typical resolution of 3–5 cm, and seven snowpack profiles were collected on the Chukchi sea ice and two at the Barrow Environmental Observatory (BEO) site, for comparison. To describe the spatial variation in DOC in the Alaskan snow cover, 27 surface snow samples were also collected across central Alaska between April 19 and 23, 2017 (Fig. S1). Detailed information is provided in the Supplementary Table S1 and Text S1.

2.2.

### Laboratory analyses

2.2.1.

#### Measurements of DOC, BC, and major ions

DOC concentrations were determined using a total organic carbon (TOC) analyzer (Shimadzu Corp, Kyoto, Japan), after filtering samples through a 0.45 μm pore-sized polytetrafluoroethylene membrane filter. The detection limit of the TOC analyzer was 15 μg L−1. The TOC analyzer can also measure total nitrogen content (TN) with a standard deviation of less than 3%.

BC concentrations in the snow were measured using a single particle soot photometer (SP2, Droplet Measurement Technology, Boulder, USA), with the detection limit of 0.3 fg/particle. The major cations (Ca2+, Mg2+, Na+, K+, NH4+) and anions (Cl, SO42–, NO3) were measured using the Dionex-600 and the Dionex-2500 Ion Chromatograph Systems (Dionex, USA), respectively. The detection limit was 1 µg L−1 and the precision errors was better than 5%. Detailed information is provided in Supplementary Text S1. Blanks of DOC and TN measurement in this study is listed in Table S2.

2.2.2.

#### Light absorption measurements

The light-absorbing properties of DOC were analyzed using the SpectraMax M5, an ultraviolet/visible absorption spectrophotometer that measures wavelengths between 200–800 nm with a resolution of 5 nm. During sampling, pure water of MilliQ is used as the blank sample. The mass absorption cross section of DOC (MACDOC365) was calculated using the Lambert-Beer law (Bosch et al., 2014; Kirillova et al., 2014; Yan et al., 2016):

((1))
${\mathrm{MAC}}_{\mathit{DOC}}=\frac{-ln|\frac{{I}_{i}}{{I}_{0}}|}{\left(C\right)×L}=\frac{A}{\left(C\right)×L}×\mathrm{ln}\left(10\right)$
where I0 and Ii represent the light intensities of the transmitted light and incident light, respectively; A are the absorbances derived directly from the spectrophotometer at a wavelength of 365 nm; (C) stands for the concentration of DOC (mg L−1); and L is the absorbing path length (0.01 m).

The absorption Ångström exponent (AAE) was calculated using the equation (Eq. 2) (Kirillova et al., 2014):

((2))
$\frac{A\left({\lambda }_{1}\right)}{A\left({\lambda }_{2}\right)}={\left(\frac{{\lambda }_{2}}{{\lambda }_{1}}\right)}^{\mathit{AAE}}$
where λ (nm) represents the wavelength; and A(λ) can be derived directly from the spectrophotometer at λ nm. Measurements of AAE can be obtained at wavelengths (λ) between 330 to 400 nm (Chen & Bond, 2010; Kirillova et al., 2014).

Studies show that BC in snow produces surface warming in the Arctic, further accelerating the snow and ice melting process (Flanner et al., 2007; Dou et al., 2012; Bond et al., 2013; Lin et al., 2014; Qian et al., 2015). Using a simplistic model (Eq. 3), we estimated the RF caused by DOC relative to BC in snow. The fraction (f) can be calculated as:

((3))
where λ (nm) stands for the wavelength; I (λ) is the clear sky Air Mass 1 Global Horizontal (AM1GH) irradiance spectrum on the Earth’s surface, which includes light-absorptive effects of the ozone and water, among others (Levinson et al., 2010); MACDOC365 and MACBC550 represent MAC of DOC at 365 nm and MAC of BC at 550 nm, respectively; hABL is the vertical height of the atmospheric boundary layer; and AAEDOC and AAEBC are the AAE of DOC and BC, respectively. In this model, MACBC550 = 7.5 ± 1.2 m2 g−1 (Bond and Bergstrom, 2006), AAEBC is set to be 1 (Bond and Bergstrom, 2006; Doherty et al., 2010; Kirillova et al., 2014), and hABL is set at 1000 m (Bosch et al., 2014; Kirillova et al., 2014).

Specific ultraviolet absorbance (SUVA254 with unit of L mg C−1 m−1) was defined by dividing the decadal ultraviolet-visible absorption coefficient per meter at the wavelength (λ) of 254 nm by the DOC concentration. SUVA254 can reflect the general absorptivity at λ = 254 nm for DOC, which typically represents an index of DOC aromaticity (Weishaar et al., 2003; O’Donnell et al., 2016).

2.3.

### In-situ DOC bioavailability experiment

The DOC bioavailability (BDOC) experiment was conducted during the fieldwork between May 7 and 17, 2017, at the sea ice surface near Barrow. In brief, surface snow samples were collected in pre-combusted (550° C for 6 h) aluminum basins. The samples were melted in the field and filtered through pre-combusted glass fiber filters (GF/F 0.7 µm) into 12 pre-cleaned 60 mL polycarbonate bottles that were wrapped in three layers of aluminum foil to avoid solar irradiation. Two samples were refrigerated immediately after filtering to obtain initial DOC concentrations and the others were placed outside at the sampling site A.

2.4.

### Principal component analysis and backward air mass trajectory analysis

Principal component analysis (PCA) is a statistical procedure that uses an orthogonal transformation to convert a set of observations of possibly correlated variables into a set of values of linearly uncorrelated variables called principal components. The resulting vectors are an uncorrelated orthogonal basis set. Zhang et al. (2016) used PCA in discussions of snow chemistry.

Back trajectories were calculated using the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model (Draxler and Rolph, 2010), for the days with recorded precipitation, to determine source regions for air masses arriving in Barrow. Additionally, footprint analyses were also performed using the European Centre for Medium-Range Weather Forecasts (ECMWF) analysis fields, with the Lagrangian analysis tool LAGRANTO (Wernli and Davies, 1997; Sprenger and Wernli, 2015) launched every six hours during the study period.

3.

## Results and discussion

3.1.

### DOC, TN, BC, and concentrations of major ions

3.1.1.

#### Spatial distribution

The spatial distribution of DOC, TN, BC, and major ions is shown in Fig. 1. We found lower concentrations of DOC in Barrow (∼ 0.17 mg L−1) and southern Alaska (near Denali National Park and Preserve) (∼0.24 mg L−1), and higher values occurring in middle Alaska (near Fairbanks) (∼0.33 mg L−1). In terms of TN, higher values were observed in northern Alaska (near Prudhoe Bay) (∼0.15 mg L−1), near the coastal regions; average TN in snow of Barrow is ∼0.13 mg L−1. The concentration of BC in snow was variable across Alaska and exhibited different spatial patterns, showing higher values in Barrow (∼1.55 μg L−1) and central Alaska (∼4.63 μg L−1), and lower values in eastern Alaska (∼0.27 μg L−1). Because of the impact of sea salt aerosols and mineral dust, larger concentrations of major ions appeared in the snow in northern Alaska, especially Cl (∼187 μg L−1), Ca2+ (1753 μg L−1), and Mg2+ (64 μg L−1). The proportion of major ions exhibited differences from the south to the north. In Barrow and southern Alaska, Na+ and Cl- dominated the total ionic concentrations; while in the central to northern Alaska, Ca2+ is the predominant ion (Fig. 1d). These different percentage of major ions indicated the different sources (natural or anthropogenic origins) of these major ions across the study area.

Fig. 1.

Distribution of (a) DOC (mg L−1), (b) TN (mg L−1), (c) BC (μg L−1), and (d) major ions (μg L−1) in Alaskan snow cover in the spring of 2017. (Detailed sampling information can be referred to supplementary Table S1.).

Compared to other studies, DOC concentrations showed a wide range from the values reported in different areas (Table 1). Spatially, the lowest DOC content (about one order of magnitude lower) in snow was reported from the Greenland Ice Sheet (0.06 mg L−1) and Antarctica (Legrabd et al., 2013; Musilova et al., 2017). Our results were comparable to those from glacial surface snow in other Alaskan regions (Fellman et al., 2015), but lower than those reported for snow cover and mountain glaciers in central Asia (Liu et al., 2016; Yan et al., 2016; Li et al., 2018; Zhang et al., 2018, 2019; Gao et al., 2020). This indicated that the DOC in the snow of central Asia was more likely to be influenced by anthropogenic pollutant inputs and desert-sourced mineral dust deposition (Yan et al., 2016). However, the DOC concentration in a glacier in the French Alps was reported to be a little lower than that in central Asia (Legrand et al., 2013), which was comparable to that in our study. Ice core samples indicated that DOC in snow were extremely low Antarctic Ice Sheet (Legrand et al., 2013). However, DOC in dry valley glaciers in the Antarctica, as reported by Hood et al. (2015), was more abundant; this value may have been affected by in situ productivities (Stibal et al., 2012) and may not be representative of the Antarctic Ice Sheet as a whole.

3.1.2.

#### Temporal and vertical variations

Results from early April to late May, 2017, showed fluctuations in DOC, TN, and BC concentrations in snow cover on the sea ice of Barrow, but exhibited no significant trends (Fig. S2a). The parameters appeared to show consistent variations, except on April 27, when DOC showed a higher value, indicating different sources, transport, or deposition processes for DOC, TN, and BC. The lower values (e.g., on April 14 and May 5) were associated with the occurrence of fresh snow. Variations in DOC content differed from the variations in the types of major ions in the snow (Fig. S2b). Especially in late May, when the snow began to melt, DOC concentrations tended to accumulate while the major ions decreased due to elution. Nitrate (NO3) in snow strongly differed compared to the other ions, implying the effects of a post-deposition process. Due to the short sampling term of this study, the seasonality of DOC in the snow at Barrow cannot be discussed further.

The snowpack profiles on the sea ice at Barrow indicated a generally consistent variations (Fig. 2). Overall, the DOC content of the surface snow (0.17 mg L−1) was higher than the subsurface snow (0.14 mg L−1) (Fig. 2). However, the profiles from SP4 (April 17, 2017) were anomalous, with higher values of DOC occurring in the subsurface snow as a result of snowdrift deposition during sampling. For SP1, SP3, and SP10, the highest values of DOC were observed at the bottom of the snowpack (Fig. 2a). Profiles of TN in snow were mostly uniform, with higher values in the surface snow (Fig. 2b). The BC profiles showed different variations across snowpacks (Fig. 2c). Distinct differences in TN and BC were found in SP5 and SP9, showing lower TN and higher BC values in the surface snow. Similar vertical characteristics of BC were also found in Russian Arctic snow (Dou et al., 2012). Ca2+ in snow was usually considered to be an indicator of crustal dust in the mid-latitude regions (Zhang et al., 2015), while in the Polar Regions, Ca2+ may be affected by sea salt from coastal areas (Lambert et al., 2012). Thus, nssCa2+ (non-sea salt Ca2+) were calculated in this study based on method from Lambert et al. (2012), which also showed similar variations with Ca2+ with a coefficient of 0.84 (p < 0.01) (Fig. S3). In Barrow, nssCa2+ were mainly affected by mineral dust. Both Cl and Ca2+ had higher values in the surface snow at Barrow (Fig. 2d,e), indicating that the impact of dust and sea salt aerosol on Cl in snow of the study area. Besides, Cl in snow can be affected by photochemistry within snow (Grannas et al., 2007). The NO3 and NH4+ levels in snow of coastal regions did not exhibit any systematic variations with depth (Fig. 2f), maybe due to the multi sources from primary production in coastal and open ocean waters and anthropogenic input by long range transport. In general, concentrations of DOC and other proxies (BC, TN, and major ions) showed higher values in the surface snow rather than in the subsurface snow (Fig. 3). We assumed that most particulate snow impurities resulted from wet deposition, as suggested through global modeling in studies by Flanner et al. (2007) and Koch et al. (2009). The increasing concentration towards the surface indicated that the mixing ratio of insoluble light-absorbing particulates to snow mass was higher during late spring, when we sampled, compared to earlier in the accumulation season (Doherty et al., 2013). This could be the result of an increase in dry deposition of aerosols in the snow surface at this time.

Fig. 2.

Vertical variations of DOC concentrations and other parameters recorded in the snowpack at Barrow, Alaska.

Fig. 3.

Comparison of DOC (mg L−1), TN (mg L−1), BC (μg L−1), and major ions (μg L−1) in surface and subsurface snow cover at Barrow.

3.2.

### Light-absorption properties of DOC

3.2.1.

#### AAE

The AAE is a critical parameter in determining the impact of absorbing aerosols on the Earth’s radiation balance (Pokhrel et al., 2016), which can indicate the wavelength dependence of DOC’s light-absorbing properties. The AAE cannot be directly parameterized based on fuel type, because fuel type depends strongly on the burn conditions and there are no large datasets available to compare AAE to fuel type. Using Eq. 2, we calculated the AAE of light-absorbing DOC in snow, providing important input data for the RF calculations. The fitted AAE330 − 400 values ranged from 1.16 to 4.22 for snow samples across Alaska, and from 0.23 to 7.21 for the snow at Barrow. The average AAE330 − 440 values were comparable to those of the Arctic snow (Doherty et al., 2010) (Table 2). However, AAE330 − 440 values measured in this study were lower than those reported by Voisin et al. (2012) of HULIS in Arctic snow, and DOC in Tibetan glacial snow (Yan et al., 2016), and much lower than the values for the BrC aerosol (Kirchstetter et al., 2004; Sun et al., 2007; Chen & Bond, 2010; Lambe et al., 2013; Liu et al., 2015; Srinivas et al., 2016). The diversity in the AAE values in snow and aerosols may suggest differences in the chemical composition of DOC as a result of different sources or different degrees of atmospheric processing (Kirillova et al., 2014).

Typically, urban and industrial fossil fuel emissions have AAE values of 1.0 − 1.5 (Bergstrom et al., 2007), depending on the material being burned and the burn temperature, among other conditions. For example, the AAE value can be larger than 2 for aerosols from inefficient coal burning (Bond, 2001), and tends to be higher (1.5 − 2.5) for biomass-burning aerosols (Kirchstetter et al., 2004; Bergstrom et al., 2007). Mineral dust was common in Arctic snow (Doherty et al., 2010) and the AAE value for dust was poorly constrained and higher (2 − 5) than that for fossil fuel emissions (Bergstrom et al., 2007; Müller et al., 2009). The AAE values that we observed for snow indicated that they were heavily influenced by biomass-burning pollutants or dust.

3.2.2.

#### Macdoc365 and SUVA254

The mass absorption cross section at 365 nm for DOC (MACDOC365) was selected in order to avoid interference from non-organic compounds such as nitrate and to be consistent with previous studies (Kirillova et al., 2014; Cheng et al., 2016; Yan et al., 2016). The average MACDOC365 value was 0.32 ± 0.24 m2 g−1 for the surface snow at Barrow, and 0.37 ± 0.32 m2 g−1 across Alaska (Table 2). Both values were lower than those for water-soluble organic carbon in an outflow from northern China (Kirillova et al., 2014) and for a receptor in the Indian Ocean (Bosch et al., 2014), and comparable to that for precipitation in the Tibetan Plateau (Li et al., 2017) (Table 2). These values were also lower than the MAC reported by Voisin et al. (2012) for Arctic snow, and the DOC in glacial snow of the Tibetan Plateau (Yan et al., 2016), as well as those for typical urban aerosol in Beijing (Cheng et al., 2016) or biomass burning (Kirchstetter et al., 2004).

The MACDOC365 values from different sources, secondary contributions, or photobleaching processes vary widely (Kirillova et al., 2014). Generally, MACDOC365 from biomass smoke can reach 5.0 m2 g−1 (Kirchstetter et al., 2004), while that from secondary organic carbon ranges from 0.001 to 0.088 m2 g−1 (Lambe et al., 2013). This indicates that, at Barrow, DOC in the snow may be affected both by biomass burning and secondary organic carbon, as evidenced by the interval MACDOC365 values. In addition, mineral dust-sourced DOC may contribute to higher MACDOC365 values in snow samples as compared to the MAC values associated with secondary organic aerosols (Yan et al., 2016).

In general, MACDOC365 values observed in the Arctic snow had a negative correlation with DOC concentrations (Fig. 4a), indicating that hydrophobic DOC attached to particles (Meyer et al., 2009). The negative relationship between the MACDOC365 and the AAE330 − 400 values (Fig. 4b) also suggested that the presence of more strongly absorbing DOC might contribute to lower AAE values, as observed in previous studies (Bosch et al., 2014; Kirillova et al., 2014; Yan et al., 2016). Temporal variations in MACDOC365 showed a rough increase during the sampling period, especially in late May when the snow began to melt (Fig. 4), which may indicate that snow melting can cause an increase of light-absorbing DOC. The above phenomenon had also been detected in glaciers in the Tibetan Plateau (Hu et al., 2018) and in snow cover of the Chinese Altai Mountains (Zhang et al., 2019).

Fig. 4.

Relationship between (a) MACDOC365 values and DOC concentrations, (b) MACDOC365 values and AAE330-400, and (c) SUVA254 and DOC concentrations for Alaskan snow samples, and temporal variation in MACDOC365 values (d), and SUVA254 Values (e) in snow cover on the sea ice at Barrow during the study period of April to May in 2017.

SUVA254 was mainly affected by the percent aromaticity of dissolved organic matter (Weishaar et al., 2003), which can be used to obtain information on organic matter in water, in terms of aromaticity (Assaad et al., 2017). SUVA254 for the snow at Barrow ranged from 1.1 to 6.8 L mg C−1 m−1, with an average of 3.5 ± 1.3 L mg C−1 m−1, which showed a negative exponent with DOC concentrations (Fig. 4c). There was no significant trend for SUVA254 during the study period (Fig. 4e), indicating the consistent presence of aromaticity in snow.

3.2.3.

#### Fraction of RF by DOC relative to BC

Based on Eq. 3, the fractions of RF by DOC relative to that by BC were calculated. Ratios of RF by DOC to BC ranged from 0.86% to 9.08% (average at 2.28% ± 1.53%) for Arctic snow (Supplementary Table S3). Temporal variations in the ratios showed no obvious trends (Fig. 5b), which was negatively correlated with MACDOC365 in the snow at Barrow (Fig. 5c). This indicated that, during snow melting when high values of MACDOC365 were observed, the RF caused by enriched BC in the snow enhanced much more than that of the DOC. The high RF ratio may be a result of the high ratio of DOC to BC concentration (DOC/BC) (Fig. 5d), which has also been noted in glacial snow in the Tibetan Plateau (Yan et al., 2016).

Fig. 5.

Variation for ratios (a) of RF by DOC to BC, and relationships between ratios of RF by DOC to BC and (b) MACDOC365 and (c) ratios of DOC/BC in the snow at Barrow, Alaska.

Recently, the role of absorbance of light-absorbing OC in snow has been highlighted. Doherty et al. (2010) noted that OC in Arctic snow was typically responsible for ∼40% of the absorption of visible and ultraviolet light. RF by DOC relative to that of BC ranged from 2.1% to 30.4% (average at 9.5% ± 8.4%) in glacial snow in the northern Tibetan Plateau (Yan et al., 2016). Simulations by Lin et al. (2014) indicated that the RF of OC deposited in terrestrial snow and sea ice accounted for 24% of the forcing by BC. In this study, the RF in snow by DOC relative to BC was lower than the observations made by other researchers; this may be due to the low DOC concentrations and the high latitude. Further investigation is needed to understand the post-deposition processes that may reduce the light-absorption of OC, including mixing with other non-absorbing constituents, scavenging efficiency, sublimation, or photochemical reactions.

3.3.

### Bioavailability of DOC

BDOC experiments in the snow found that DOC levels decreased exponentially over time (r2 = 0.81) (Fig. 6), with approximately 60% of DOC degrading within 15 days. Taking into consideration the different experiment conditions, our results showed higher BDOC consumption compared to studies in European Alpine glaciers (59% within 50 days, at 4° C) (Singer et al., 2012) and Tibetan glaciers (26.7% within 15 days at 3.8° C) (Spencer et al., 2014; Yan et al., 2016). A study by Lawson et al. (2014), in a large Greenland Ice Sheet catchment during contrasting melt seasons, indicated a 46 − 64% decline in DOC during glacial incubation (glacier samples only) and exceeded the average DOC decline (17%) in the glacial controls over the same time period. These results suggest the utilization of DOC by microbial heterotrophs in snow and glacier regions (Hood et al., 2009; Yan et al., 2016). For Alpine glaciers, Singer et al. (2012) found that the percentage of BDOC increased with Δ14C DOC, and showed that radiocarbon age explained only 23% of the variation in BDOC across the studied glaciers. This indicated that complementary biogeochemical factors (for example, chemical composition and sources) can influence the bioavailability of snow-linked DOC. Currently, we did not analyze the Δ14C DOC and the biogeochemistry of DOC, thus limiting further discussion.

Fig. 6.

Exponential decrease of DOC concentrations in snow cover on Chukchi sea ice during the biodegradation experiment in May of 2017 at Barrow.

3.4.

### Potential sources of DOC

Backward trajectory analyses indicated that the trajectories arriving at Barrow spread quite a lot, but essentially remain in the northern and southern latitudes (Fig. 7a,b), especially from the Arctic Ocean and Alaskan regions, which may have some distinctly different characteristics. For example, backward air mass trajectories showed that the air masses arriving at Barrow on April 17, 2017 had passed over the Canadian Arctic and Alaska (Fig. 7c), while air masses arriving at Barrow on May 2, 2017 came from the Arctic Ocean (Fig. 7d). This evidenced strong marine and terrestrial sources for organic carbon in Barrow. Distant aerosol transport and the deposition of terrestrial carbon have also been reported in snowpacks and ice cores in other remote locations, such as Greenland (Grannas et al., 2004) and the Russian Arctic (Grannas et al., 2006). Terrestrial DOM in a glacial ecosystem was found in a highly bioavailable form (Lawson et al., 2014; Fellman et al., 2015) and was also a significant driver of bacterial production in aquatic systems (Pollard & Ducklow, 2011).

Fig. 7.

Footprint analyses of the backward trajectories launched over the snow sampling site at Barrow, Alaska, during April and May in 2017 for pass by (a) 48 h and (b) 120 h, respectively, and HYSPLIT model seven day back trajectories for air masses arriving in Barrow on (c) April 17, 2017 and (d) May 2, 2017. (Different colors of backward trajectories in c and d represent the air mass arrived at the different date.).

PCA analysis of chemical impurities in the snow at Barrow revealed that the first component (PCA1), which accounted for 34.41% of the variance, was dominated by Cl, Na+, and K+ (Table 3), indicating that they arise mostly from marine sources (i.e., natural sources). PCA2 accounted for 24.04% of the total variability and consisted primarily of anthropogenic impurities (SO42−, NO3, NH4+, and BC). DOC was largely loaded in PCA3 with Ca2+ and TN. When we consider the snow across Alaska, PCA1 (35.48%, mainly loaded with Ca2+, Mg2+, SO42−) and PCA2 (18.92%, mainly loaded with DOC, Na+, K+, and Cl) represented the crustal origins and impacts of marine aerosols/biomass burning, respectively. Beside the impact of sea salt and anthropogenic impact, sulfate in snow may be affected by mineral dust input (Zhang et al., 2016). In this study, nssSO42 was positively correlated with sulfate concentrations (Fig. S4a). There was almost no relationship between DOC and SO42- and nssSO42- (Fig. S4b,c). While, positive correlation was also found between nssSO42- and Ca+, and Mg2+, further indicating the important impact of mineral dust. Previous study of 87Sr/86Sr and εNd(0) data from the sea ice snow samples from the Arctic Ocean and snow samples from Barrow (Alaska) also suggested potential long-range dust transport to the Arctic region (Du et al., 2019).

Fellman et al. (2015) found that DOC in snow collected across Alaska contained high amounts of K+, indicating biomass burning sources, as evidenced by the source-diagnostic dual-isotope composition of the DOC. Dou et al. (2017) observed that polyaromatic hydrocarbons and hopanes were below the detection limits, indicating that fossil fuel combustion (anthropogenic sources) contributes little to the OC in the Barrow snow. Biomass burning is widespread in Alaska (Kelly et al., 2013; Mouteva et al., 2015), and appears to contribute to the OC in the Barrow snow (Dou et al., 2017). Biplot analysis indicated that DOC in the snow in Barrow had a negative relationship with Na+, K+, Cl, Ca2+, Mg2+, and SO42− (Fig. 8a). In terms of snow across Alaska, DOC had a positive correlation with Na+, K+, and Cl-. Since the occurrence of fires in Alaska was quite low (Fig. S5) and many trajectories pass over the Arctic Ocean (Fig. 7a,b), the impact of biomass burning on DOC deposited in the snow at Barrow may not be significant. However, the lack of data related to the chemical composition of the DOM in snow and its carbon isotopic signatures has confined this study to quantifying the different contributions to carbon aerosols.

Fig. 8.

Biplot analyses of DOC with other parameters for (a) snow samples in Barrow, and (b) snow samples across Alaska.

Snowpack cycling of DOC can be important on the local, regional, and global scales (McNeill et al., 2012). OC depositions in snow and sea ice may originate from plant material, microorganisms and their excretions, aerosols, or volatile organic compounds deposited from the atmosphere (Pautler et al., 2011). Fellman et al. (2015) used carbon isotopic signatures (δ13C and Δ14C) combined with a three-source mixing model and found that the DOC deposited in snow across the icefields of southeast Alaska came mainly from fossil fuel combustion products (43 − 73%), and to a smaller extent from marine (21–41%) and terrestrial (1 − 26%) sources. In Antarctica, OC in snow was consistent with materials from microbial sources as well as terrestrial inputs of vascular plant-derived materials (Antony et al., 2014).Carbonaceous isotopes in the snow should be further investigated. Simulations coupled with an emission inventory should also be considered in the future.

3.5.

### Climatic implications of DOC for the Arctic region

The Arctic is of special importance to the world due to its amplification of temperature, consequently changes in atmospheric and oceanic circulation, and as well as reductions in snow cover and sea ice (ACIA, 2004; Derksen et al., 2015). Such changes can affect carbon and nitrogen cycling, as well as methane emissions (AMAP, 2015; Bokhorst et al., 2016), and are thus significant effects for impact assessments and developing adaptation strategies.

Sea ice and snow cover contain BC and other light-absorbing impurities (Flanner et al., 2009; Doherty et al., 2010; Dou et al., 2012; Legrand et al., 2013; Dou et al., 2017; Marks et al., 2017). Light-absorbing impurities in snow can dominate the absorption of solar radiation, reduce the albedo relative to clean snow, contribute to the surface energy budget, and lead to earlier snowmelt (Flanner et al., 2007; Doherty et al., 2010; Derksen et al., 2015; Skiles et al., 2018). In a model simulation by Lin et al. (2014), the OC reduced the snow albedo, resulting in RF of organic aerosols deposited in land snow and sea ice ranging from +0.0011 to +0.0031 W m−2, which could potentially enhance the snowmelt (Flanner et al., 2009; Dumont et al., 2014).

Our findings indicated that light-absorbing DOC has the potential ability to absorb more solar energy and accelerate snow melting. Schmale et al. (2017) reported that snow darkening caused by BC can enhance summer snowmelt rises of up to 6.3% (7 cm w.e. a−1) on the glaciers of central Asia, based on an albedo reduction and snowmelt model. Zhang et al. (2017) also calculated that BC in snow enhanced glacial melt by approximately 15% in the south-eastern Tibetan Plateau. In the Arctic region, a simulation by the Community Earth System Model reported that the deposition of light-absorbing impurities caused 1 − 2° C of surface warming over a large area of the Arctic Ocean and subarctic seas in autumn and winter (Goldenson et al., 2012). Light-absorbing impurities contributed between 1.6%−5.1% to the reduction of albedo with snowmelt during April and May, 2015 (Dou et al., 2017). Changing snow cover and sea ice in the Arctic region has been linked to cold extremes in the mid-latitude and northern continents/subarctic areas (Cohen et al., 2013). The amplified Arctic warming over the past decade has significantly contributed to a continual global warming trend (Huang et al., 2017; IPCC, 2019). Quantifying the impact of light-absorbing impurities on changes in snow cover and sea ice, and understanding the mitigation strategies for these effects are immediate concerns. Improved understanding of processes, additional Arctic observations, and further modeling efforts based on observed data are required to elucidate the connections with the Arctic (Cohen et al., 2014).

4.

## Conclusions

This study investigated the characteristics of DOC in snow cover across the Alaska and on sea ice at Barrow, further confirming that DOC in snow can significantly affect light-absorption in the Arctic regions. Average DOC concentrations in snow cover at Barrow (∼0.17 mg L−1) and in southern Alaska were lower than those measured in central Alaska (0.33 mg L−1). DOC levels in snow cover (0.17 − 0.30 mg L−1) were comparable to those in glacial surface snow in Alaska, but lower than those in mountain glaciers in Central Asia. No significant temporal variation in DOC concentrations was noted at Barrow, due to the observations were recorded over a short time period. TN and BC in snow of Alaska showed spatial distribution with lower concentrations in the eastern Alaska (with 0.05 mg L−1 and 0.24 mg L−1 for TN and BC respectively). While major ions in snow show higher total concentrations in the northern Alaska, with Na+ and Cl- dominated in Barrow and southern Alaska and Ca2+ dominated in the central to northern Alaska due to different source impacts. For snowpack, vertical variations of DOC and other related proxies generally showed higher values in the surface snow than that in the sub-surface snow.

The light-absorption results showed that AAE330 − 400 values of DOC in snow were 2.11 and 2.34 at Barrow and across Alaska, respectively, which were comparable to those in the Arctic snow, but lower than the values of HULIS in Arctic snow, DOC in Tibetan glaciers, and BrC aerosols. MACDOC365 values of DOC were lower than those of aerosol, which was comparable to those for precipitation in the Tibetan Plateau. MACDOC365 was generally negatively correlated with DOC concentrations, indicating hydrophobic DOC attached to particles in the snow. The results show that the ratios of RF by DOC relative to BC ranged from 0.86% to 9.08% for Arctic snow.

DOC concentrations decreased exponentially over time, with approximately 60% of DOC degrading within 15 days. This highlighted the need to further quantify and determine the source of organic matter input to snow meltwater, impacted ecosystems, and heterotrophic activity. Based on the PCA analysis and the modeled back trajectories, we found that marine and terrestrial sources contributed strongly to DOC deposition on the Alaskan snow cover. Our findings also emphasized the importance of the light-absorbing properties of DOC in snow and its climatic implications for the Arctic.

This study focused on the characteristics of DOC in Arctic snow and highlighted its potential roles on RF. Further studies need to be carried out to understand these roles more thoroughly, including expanded sampling across seasons and process studies. Monitoring the snowmelt process at Barrow will help us better understand the redistribution and scavenging of DOC in snow, and a controlled experiment will help verify the radiative-transfer modeling of albedo reduction attributed to DOC.