Spectral attenuation coefficients from measurements of light transmission in bare ice on the Greenland Ice Sheet
- 1Department of Geography, University of California, Los Angeles, Los Angeles, California, 90027, USA
- 2Pacific Northwest National Laboratory, Richland, Washington, 99354, USA
- 3Institute at Brown for Environment and Society, Brown University, Providence, Rhode Island, 02912, USA
- 4Department of Earth, Environmental and Planetary Sciences, Brown University, Providence, Rhode Island, 02912, USA
- 5Department of Geography, Rutgers, The State University of New Jersey, New Brunswick, New Jersey, 08854, USA
- 6NASA Goddard Institute for Space Studies, New York, New York, 10025, USA
- 7Lamont-Doherty Earth Observatory, Columbia University, New York, New York, 10964, USA
Correspondence: Matthew G. Cooper (email@example.com)
Light transmission into bare glacial ice affects surface energy balance, biophotochemistry, and light detection and ranging (lidar) laser elevation measurements but has not previously been reported for the Greenland Ice Sheet. We present measurements of spectral transmittance at 350–900 nm in bare glacial ice collected at a field site in the western Greenland ablation zone (67.15∘ N, 50.02∘ W). Empirical irradiance attenuation coefficients at 350–750 nm are ∼ 0.9–8.0 m−1 for ice at 12–124 cm depth. The absorption minimum is at ∼ 390–397 nm, in agreement with snow transmission measurements in Antarctica and optical mapping of deep ice at the South Pole. From 350–530 nm, our empirical attenuation coefficients are nearly 1 order of magnitude larger than theoretical values for optically pure ice. The estimated absorption coefficient at 400 nm suggests the ice volume contained a light-absorbing particle concentration equivalent to ∼ 1–2 parts per billion (ppb) of black carbon, which is similar to pre-industrial values found in remote polar snow. The equivalent mineral dust concentration is ∼ 300–600 ppb, which is similar to values for Northern Hemisphere warm periods with low aeolian activity inferred from ice cores. For a layer of quasi-granular white ice (weathering crust) extending from the surface to ∼ 10 cm depth, attenuation coefficients are 1.5 to 4 times larger than for deeper bubbly ice. Owing to higher attenuation in this layer of near-surface granular ice, optical penetration depth at 532 nm is 14 cm (20 %) lower than asymptotic attenuation lengths for optically pure bubbly ice. In addition to the traditional concept of light scattering on air bubbles, our results imply that the granular near-surface ice microstructure of weathering crust is an important control on radiative transfer in bare ice on the Greenland Ice Sheet ablation zone, and we provide new values of flux attenuation, absorption, and scattering coefficients to support model development and validation.
Understanding the transmission, absorption, and scattering of light in ice is important for snow and ice energy balance modelling (Brandt and Warren, 1993), lidar remote sensing of snow surface elevation and grain size (Deems et al., 2013; Yang et al., 2017), primary productivity beneath sea ice (Frey et al., 2011; Grenfell, 1979), biophotochemistry in ice and snowpack (France et al., 2011), and theoretical predictions of “Snowball Earth” palaeoclimates (Dadic et al., 2013; Warren et al., 2002). Each of these applications requires knowledge of the vertical distribution of light attenuation in ice, which for a medium (such as glacier ice) that both absorbs and scatters light is specified by the spectral attenuation coefficient:
where kabs (m−1) is the absorption coefficient, ksca (m−1) is the scattering coefficient, and all are functions of wavelength λ. This study reports on the irradiance attenuation coefficient katt of bare glacier ice in the Greenland Ice Sheet ablation zone, a critical parameter needed to calculate subsurface absorption and scattering of transmitted radiation that to our knowledge has received no direct field study.
Measurements of katt in snowpack and sea ice indicate three main sources of variation with relevance to geophysical applications. First, the magnitude of katt is primarily controlled by ice microstructure (e.g. the size, shape, orientation, and number of air bubbles, ice grains, and cracks), via its control on ksca (Dadic et al., 2013; Libois et al., 2013; Light et al., 2004, 2008). For the range of air bubble sizes (∼ 10−3–10−4) and ice grain sizes (∼ 10−1–10−3) observed in glacier ice, ksca is effectively independent of wavelength in the visible and near-infrared spectrum (Bohren, 1983; Dadic et al., 2013; Perovich, 1996). Spectrally, katt is low in the near-ultraviolet and blue-green spectrum (∼ 250–600 nm), where kabs is extremely low (< 10−8), and progressively higher for wavelengths > 600 nm, where kabs rapidly increases up to its maximum value (∼ 10−2) at the far end of the solar spectrum (Warren and Brandt, 2008). Vertically, katt is at a maximum at the incident boundary (the snow or ice surface) where a portion of upwelling radiation (i.e. transmitted flux reflected upward) escapes the ice volume before re-reflection downward. Within this near-surface optical boundary layer (Bohren and Barkstrom, 1974), attenuation rates rapidly decrease with depth to an asymptotic value as multiple scattering establishes an isotropic (diffuse) radiation field (Briegleb and Light, 2007; Warren, 1982). For fine-grained dry snow, a few centimetres of depth is typically sufficient to reach the asymptotic regime where monochromatic katt is constant (Brandt and Warren, 1993). For sea ice the depth required is typically larger and can exceed > 20 cm depending on near-surface ice microstructure and the vertical location of the refractive boundary if present (Grenfell, 1991; Grenfell and Maykut, 1977). Attenuation coefficients are also influenced by the horizontal distribution of ice type and surface cover (Frey et al., 2011), but this source of variation is not examined here.
In addition to experimental values obtained from measurements of light transmission in ice or snow, katt is obtained analytically from optical theory (Bohren, 1987; Warren et al., 2006). Light attenuation in pure ice is specified analytically by the complex refractive index:
where mre is the real part of the complex refractive index (∼ 1.31 in the visible), mim is the imaginary part, and
is the absorption coefficient of pure ice (Warren et al., 2006; Warren and Brandt, 2008).
Light attenuation in glacier ice differs from pure ice owing to compositional and structural factors that control scattering and absorption, such as the size, geometry, and vertical distribution of embedded light-absorbing particles (LAPs) and light-scattering air bubbles and ice grains of size larger than wavelength (Askebjer et al., 1997; Picard et al., 2016; Price and Bergström, 1997b; Warren et al., 2006). Analytical methods typically approximate ice and snowpack as homogeneous plane-parallel slabs of spheres having the same volume-to-surface-area ratio (i.e. optically equivalent grain size) as the collection of non-spherical ice grains and air bubbles in realistic ice (Brandt and Warren, 1993; Grenfell and Warren, 1999; Wiscombe and Warren, 1980). Mie theory is used to calculate the single-scattering properties, and two-stream radiative transfer approximations are used to calculate multiple scattering and bulk absorption in the ice volume (Bohren and Barkstrom, 1974; Mullen and Warren, 1988; Wiscombe and Warren, 1980). The single-scattering properties can also be derived from the ratio of surface area to mass (i.e. specific surface area) with or without the assumption of spherical scattering geometry (Kokhanovsky and Zege, 2004; Malinka, 2014), as applied to the highly scattering granular surface layer on sea ice (Malinka et al., 2016). Models of the prior form have been used to calculate subsurface meltwater production caused by penetration of solar radiation in ice in both Greenland (van den Broeke et al., 2008; Kuipers Munneke et al., 2009) and Antarctica (Brandt and Warren, 1993; Hoffman et al., 2014; Liston et al., 1999a, b; Liston and Winther, 2005). However, theoretical values for katt are rarely validated experimentally, and to our knowledge no such experimental values exist for near-surface glacier ice.
In addition to ice surface energy balance, understanding light attenuation in ice is important for interpreting interactions between visible-wavelength light sources and ice surfaces, for example laser altimetry measurements of ice surface elevation (Deems et al., 2013; Gardner et al., 2015; Greeley et al., 2017). The reciprocal of katt is the attenuation length, or the average distance travelled by a photon before attenuation by absorption or scattering (Ackermann et al., 2006). In the context of altimetry, the attenuation length is sometimes referred to as the penetration depth, or the average depth to which the electromagnetic signal penetrates before it is backscattered to the atmosphere (Ridley and Partington, 1988; Rignot et al., 2001; Zebker and Weber Hoen, 2000). The laser altimeter onboard Ice, Cloud, and Land Elevation Satellite (ICESat) transmitted 1064 nm laser pulses to measure the distance (range) between the satellite and ice sheet surfaces (Schutz et al., 2005). Photons with a wavelength of 1064 nm penetrate into snowpack no more than a few centimetres (Brandt and Warren, 1993; Järvinen and Leppäranta, 2013). This length scale is smaller than typical laser altimetry surface elevation errors due to ice and snow surface roughness and geolocation uncertainty (Deems et al., 2013). In contrast, the laser altimeter onboard ICESat-2 transmits 532 nm laser pulses (Markus et al., 2017). Ice is ∼ 10 times more transparent at 532 nm than at 1064 nm (Warren and Brandt, 2008), and photons at 532 nm may penetrate many tens of centimetres into glacier ice. These subsurface scattered photons may introduce a range bias in ICESat-2 surface elevation retrievals over glacier ice, similarly to radar penetration into snow (Brunt et al., 2016; Gardner et al., 2015; Greeley et al., 2017; Smith et al., 2018). To our knowledge no in situ observations of 532 nm optical penetration depth for bare glacier ice exist, precluding field validation of penetration depth obtained from theoretical radiative transfer models.
The purpose of this investigation is to provide experimental values for katt obtained from measurements of solar flux attenuation in bare ice in the Greenland Ice Sheet ablation zone and to compare them with theoretical values for katt obtained from the two-stream analytical solution (see Eq. 26 in Bohren, 1987; Schuster, 1905). We benchmark our field estimates against the two-stream solution because of its wide use in surface energy balance models applied to snow and ice. In Sect. 2 we describe the field measurements and the optical theory used to interpret the solar flux attenuation. In Sect. 3 we report values for katt obtained from our measurements, compare them with values obtained from two-stream theory, and propose a simple empirical model that accounts for enhanced near-surface attenuation. In Sect. 4 we discuss measurement uncertainty, and in Sect. 5 we discuss how our katt values differ from prior experimental values acquired in sea ice, snowpack, and deep South Pole glacial ice and the implication of these differences for modelling radiative transfer in bare glacier ice. To demonstrate the broader implications of our study, we suggest how our findings can be used to improve models for subsurface heating of ablating glacier ice.
2.1 Transmittance measurements
Ice transmittance was measured on 20 and 21 July 2018 in the Kangerlussuaq sector of the western Greenland Ice Sheet. The study site is located ∼ 1 km from the ice sheet margin at 840 m above sea level. (67.15∘ N, 50.02∘ W). Subsurface (in-ice) spectral irradiance was measured at a ∼ 0.35 nm spectral resolution in the wavelength range of 350–900 nm with an Ocean Optics® JAZ spectrometer. Light was guided from the ice interior to the spectrometer with a 3 mm diameter Kevlar-sheathed fibre-optic cable fitted inside a 2 m long insulated white PVC tube (Fig. 1). The fibre was attached at one end to an irradiance sensor consisting of a 90∘ collimating lens adapter and a remote cosine receptor (RCR) with a Spectralon® diffusing element. The RCR lens barrel was wrapped in white PTFE tape and set 2 mm out from the PVC tube exterior to act as a contact horizon between its diffusing element and the ice. The system was operated from a battery-powered computer running the Ocean Optics® OceanView software placed on a tripod platform oriented 180∘ away from the sun and at a 2.5 m horizontal distance from the measurement location.
To access the interior of the ice, holes were drilled horizontally into a 2 m high sidewall of a natural ice feature with a battery-powered hand drill fitted with a 3 cm diameter Kovacs auger bit. To drill these holes, the auger was advanced into the sidewall approximately 20 cm and levelled horizontally with a digital spirit level, and the sequence was repeated to a 2 m horizontal depth. The PVC tube–fibre-optic assembly was then inserted into the hole, RCR facing upward, and a 2 m long ruler was shimmed under the bottom of the PVC tube to ensure the RCR barrel preserved contact with the overlying ice thus minimizing stray light contamination into the RCR field of view. Ice shavings were packed around the drill hole to prevent light reflection into the hole. Spectral irradiance was measured using a 20-scan average with a 0.0228 s integration time per scan, yielding a 0.46 s total integration time per irradiance measurement. Irradiance measurements were recorded at a 1 Hz frequency for 30 s yielding 30 irradiance profiles at each depth, after which the tube was removed; the next hole was drilled, and the sequence was repeated, working from the top hole toward the bottom on 20 July and from the bottom hole toward the top on 21 July. The measurements were completed between 13:45 and 14:35 local time (UTC−3) on 20 July and between 13:09 and 14:00 on 21 July, at solar zenith angles of ∼ 48–51∘. Solar noon at this time and location is ∼ 13:26.
Background upwelling and downwelling spectral irradiances were measured continuously at 2 m height above the ice surface ∼ 3 m away from the in-ice measurements with a dual-channel Ocean Optics JAZ spectrometer. These data were recorded at a 1 min frequency using a 30-scan average with a 0.011 s integration time. Light was guided to the spectrometer via two 3 m fibre-optic cables attached to two RCRs mounted in an upward-looking and downward-looking orientation on a 2 m long horizontally levelled boom attached to a vertical mast frozen into the ice. The horizontal boom became unstable on 21 July, and the upward-looking RCR was moved to the vertical mast; the downward-looking RCR was decommissioned.
The surface-based spectrometer was calibrated for absolute irradiance in a controlled setting prior to the field experiment using an Ocean Optics HL-3P radiometrically calibrated halogen light source. During the field experiment, the in-ice spectrometer was cross-calibrated to the surface spectrometer by holding it level above the ice surface in an upward-looking orientation ∼ 3 m away from the surface spectrometer. Cross-calibration irradiance profiles were collected on 20 and 21 July immediately prior to subsequent in-ice measurements. All in-ice irradiances are cross-calibrated to the surface spectrometer as a pre-processing step prior to further analysis.
Dark-current spectra were recorded prior to each irradiance measurement as input to the OceanView automated dark-current correction module. To measure dark current, the RCR lens barrel was capped with a custom-fit opaque metal cap provided by Ocean Optics. OceanView adjusts these spectra in real time for changes in integration time and for charge leakage if detected, corrects the nonlinear analogue-to-digital response of the linear silicon charge coupled device, and applies a boxcar smoothing over adjacent pixels to further reduce noise. Following these automated corrections, the opaque cap was left in place and residual dark current (noise) was recorded with the reference spectrometer in its experimental setup as described above and with the in-ice spectrometer held level above the ice surface in an upward-looking orientation. These residual dark-current spectra are treated as systematic errors and are subtracted from all irradiance profiles as a pre-processing step prior to analysis (Fig. 2a).
2.2 Weather conditions
The 20 July experiment was conducted under low, thick cloud cover with light rain and no direct sun, ideal conditions for estimating the attenuation of diffuse light in ice. The 21 July experiment was conducted under higher, thinner cloud cover with no rain and very brief periods of intermittent direct sun (see Fig. 1). The effect of intermittent direct sun was easily identified in the in-ice irradiance measurements as a rapid increase in light intensity, which only occurred during the third measurement on 21 July. This was mitigated by averaging over the first 10 in-ice irradiance profiles for that measurement, prior to the rapid increase in light intensity, and discarding the remainder.
2.3 Ice thickness and density
The ice thickness between detector positions was measured to the nearest millimetre with a metre stick and converted to units of solid ice thickness with the relation
where Δh is in situ ice thickness between detector positions, ρ is in situ ice density, ρice is solid ice density (917 kg m−3), and Δz is solid ice thickness between detector positions. Two separate observers made 10 independent measurements of Δh. In addition, one observer made 41 replicate measurements of an ablation stake using the same metre stick, yielding a mean difference and standard error in Δh. The ice density ρ was measured on a 1.2 m ice core extracted at the measurement location with a Kovacs Mark IV corer (https://kovacsicedrillingequipment.com/, last access: 7 April 2021) (Fig. 3). The ice core was split along natural breaks into three segments that were measured to the nearest millimetre with a calliper and weighed to the nearest gram on an Acculab digital scale.
2.4 Experimental asymptotic flux attenuation coefficients and ice surface albedo
Spectral asymptotic flux attenuation coefficients are estimated by fitting a Bouguer-law exponential decay model to the in-ice irradiance profiles (Grenfell and Maykut, 1977):
where katt(λ) is the asymptotic flux attenuation coefficient, Iz(λ) is in-ice spectral irradiance at depth z, I0(λ) is background downwelling spectral irradiance, z0 is the ice surface, and
is spectral transmittance. The optical depth τz(λ) is a dimensionless path length that scales the physical thickness of a layer by its attenuation rate:
Estimates of katt(λ) for each spectral band are obtained by solving a linear equation of the form
where τ0 is a parameter (y intercept), is ice thickness, εΔz is an error term that represents ice thickness measurement uncertainty, and εΔτ is an error term that represents optical path measurement uncertainty. Equation (8) is solved by maximum likelihood estimation (MLE), which gives an unbiased estimate of the slope when measurement errors are present in both the independent and the dependent variables (see Sect. 2.9) (York et al., 2004).
The attenuation length latt(λ) is the inverse of katt(λ) and is analogous to the photon mean free path or transport length (Ackermann et al., 2006). It is equivalent to the path length in ice required to attenuate irradiance to 37 % of its incident intensity, i.e. the path length at which and τ=1:
The ice surface spectral albedo is the ratio of the upwelling spectral irradiance to the downwelling spectral irradiance:
and the broadband albedo is
2.5 Asymptotic flux attenuation coefficients
Theoretical katt(λ) values are calculated using the asymptotic solution to the δ-Eddington two-stream radiative transfer approximation (Joseph et al., 1976; Schuster, 1905):
where Qext(λ) is the extinction efficiency; reff is the effective scattering particle radius (m); g(λ) is the average cosine of the scattering angle, also referred to as the asymmetry parameter; and is the single-scattering albedo,
where σatt(λ) and σsca(λ) are the single-scattering attenuation coefficient (m−1) and scattering coefficient (m−1), respectively. Equation (12) describes light attenuation by multiple scattering and absorption in a homogeneous plane-parallel slab of absorbing spheres far from any boundaries (Bohren, 1987).
To estimate reff, Eq. (12) is inverted and solved by iteration for the value of reff that minimizes the difference between measured and calculated katt at λ=600 nm. This method assumes that all absorption at 600 nm is due to ice (Warren et al., 2006). If absorption was influenced by LAPs, reff would be overestimated. Values for Qext(λ), g(λ), and are obtained from Mie scattering algorithms provided as MATLAB® code (Mätzler, 2002) with input m(λ) from Warren and Brandt (2008). The Mie solutions at each wavelength are integrated over a Gaussian size distribution (N=1000) of scattering radii to eliminate ripples associated with Bessel function solutions to the Mie equations (Gardner and Sharp, 2010). The optimal reff values are ∼ 9.3 and ∼ 10.6 mm with corresponding specific surface areas ∼ 0.35 and ∼ 0.31 m2 kg−1 for 20 and 21 July experimental values, respectively.
2.6 Flux absorption coefficients
Warren et al. (2006) developed a method to estimate kabs for pure ice (i.e. ) from measurements of flux attenuation in snow in Antarctica. The method relies on three assumptions: (1) the value of at the reference wavelength (λ0=600 nm) is known accurately, (2) the value of katt at λ0 is not affected by LAPs in the measured snow or ice, and (3) varies so little as to be effectively independent of wavelength in the spectral range considered (here the near-UV and visible spectra). Warren et al. (2006) verified the validity of these assumptions for the spectral range of 350–600 nm and obtained the following relation (Eq. 15 of that paper) between flux attenuation and flux absorption:
Equation (14) assumes that kabs is not affected by LAPs at the reference wavelength (600 nm), but the relation can be used to estimate kabs at any other wavelength, including those where absorption is affected by LAPs. At those wavelengths, Eq. (14) will predict values for kabs higher than if LAPs are present in the measured snow or ice volume, due to the influence of LAPs on katt.
The inferred kabs values can be related to a mass absorption cross section (MAC) (Doherty et al., 2010):
where β is the spectral MAC (m2 kg−1) and c is the mass mixing ratio of LAPs in the ice volume (g LAPs g−1 ice). We exploit this to interpret differences between our theoretical and experimental values of katt on the basis of differences between (Warren et al., 2006) and the kabs values that we obtain for glacier ice from Eq. (14). To provide context, we use representative values of β for black carbon βBC and insoluble mineral dust (hereafter “dust”) βdust to estimate corresponding equivalent mass mixing ratios ceqBC and ceqdust (Di Mauro et al., 2017; Doherty et al., 2010). The “equivalent” mass mixing ratio is the mass mixing ratio of each LAP species required to explain the difference between and our inferred kabs values at a reference wavelength, and it follows a similar approach to that used to infer LAP absorption in snowpack (Tuzet et al., 2019). For βBC, we use 6 m2 g−1 as a representative MAC at 550 nm and an absorption Ångström exponent range of 0.8–1.9 to scale this value to 400 nm (Doherty et al., 2010). For βdust, we use 0.013 m2 g−1 at 550 nm (Di Mauro et al., 2017) and an absorption Ångström exponent range of 2–3 (Doherty et al., 2010). We note that these descriptive estimates provide context for discussion; actual LAP species concentrations were not measured.
2.7 Near-surface effects
Equations (7) and (12) are applicable at distances far enough from the incident boundary (here the ice surface) that the radiation field is diffuse and katt is constant with depth. Near the ice surface the radiation field is converted via multiple scattering from direct to diffuse flux, and attenuation may be enhanced by direct reflection, enhanced scattering, and/or absorption by the granular near-surface ice microstructure or specular reflection at the ice surface, depending on its roughness (Dadic et al., 2013; Light et al., 2008; Mullen and Warren, 1988). To account for non-diffuse near-surface attenuation, we define a piecewise optical depth:
where k′ is an effective attenuation coefficient for the near-surface non-diffuse layer and z′ is a depth chosen to partition this layer from the interior diffuse region. We estimate k′ with a centred finite-difference form of Eq. (7):
Here, Δz′ is 12 cm and is the 12 cm in-ice irradiance measured on 20 July. Accordingly, the asymptotic attenuation length (Eq. 9) is distinguished from an effective penetration depth dλ to include the effect of near-surface attenuation. The attenuation length is the depth at which τ=1. Setting τ(λ) to 1 in Eq. (16) and solving for z yields
Equation (16) gives estimates of spectral transmittance that account for non-diffuse near-surface attenuation but relies on knowledge of k′, which is sensitive to the spectral composition and directional distribution of I0 and the structure and composition of the near-surface ice (Grenfell and Maykut, 1977; Light et al., 2008). To generalize the magnitude of near-surface attenuation, we calculate the fraction of downwelling spectral irradiance that transmits the non-diffuse layer weighted by the downwelling spectral irradiance:
The χ0 parameter is analogous to the i0 parameter introduced by Grenfell and Maykut (1977) to partition the fraction of solar irradiance absorbed in the upper 10 cm of sea ice, which they termed the “surface scattering layer” (SSL), and the ice interior, in which radiation is exponentially absorbed at a constant rate:
The i0 parameter has been widely adopted in energy balance models of glaciers and sea ice to compute subsurface flux divergence (heating rates) when radiation penetration is considered important (Bintanja and Van Den Broeke, 1995; Hoffman et al., 2014; Holland et al., 2012). For example, the sea ice component of the Community Earth System Model (CESM) uses i0=70 % for the visible (200–700 nm) and i0=0 % for the infrared (700–5000 nm) (Briegleb and Light, 2007). The important distinction is that i0 partitions the absorbed flux whereas χ0 partitions the downward flux (Brandt and Warren, 1993). For both χ0 and i0, we set Δz′ to 10 cm for consistency with prior work (Grenfell and Maykut, 1977; Light et al., 2008; Maykut and Untersteiner, 1971).
2.8 Monte Carlo simulations of detector interference
We developed a Monte Carlo radiative transfer model to estimate the effect of detector interference on measured irradiances and fitted katt values, following methods developed to simulate light propagation in biological tissue, ocean waters, and sea ice (Leathers et al., 2004; Light et al., 2003; Wang et al., 1995). Photon scattering is specified by a Henyey–Greenstein scattering phase function with single-scattering properties Qext(λ), g(λ), and inferred from our optical measurements (Sect. 2.5). A complete technical description is given in the Supplement, where model accuracy is verified by comparison with benchmark solutions to the radiative transfer equation (van de Hulst, 1980).
In the Monte Carlo simulations, photons are launched from an irradiance sensor on a detector rod with dimensions identical to those reported in this study. In the ideal (baseline) simulation, photons originate from an isotropic point source and propagate through ice until they transmit the surface or are terminated by absorption. Detector interference is investigated by repeating the Monte Carlo simulations with an ideal cosine source function describing the angular response to radiance of the RCR and with a non-ideal (empirical) angular response function (Fig. 2), with and without scattering and absorption interference by the PVC detector rod. The detector rod albedo is calculated from the absorption spectra of polyvinyl chloride (Zhang et al., 2020); scattering by the rod is assumed isotropic. The Monte Carlo model is integrated over 10 000 interactions at nine wavelengths in 50 nm increments from 350 to 750 nm, allowing us to fit the wavelength dependence of the estimated systematic uncertainty in simulated katt values.
2.9 Uncertainty propagation
Unless stated otherwise, all statistical uncertainties reported in this paper are standard errors that correspond to 68 % confidence intervals around the mean (Taylor and Kuyatt, 1994). For an individual measurement with standard deviation si and sample size N≥30, the standard error is . For N<30, standard errors are scaled by a critical t value drawn from Student's t distribution. Standard errors for combined quantities are propagated in quadrature and hereafter referred to as combined uncertainty. The combined uncertainty for spectral irradiance I(λ) is
where is the standard deviation of the high-frequency irradiance spectra before dark-noise correction and σD is the standard deviation of the high-frequency dark-noise spectra. An analogous procedure is used to estimate the combined uncertainty for calibrated irradiance. The combined uncertainty for spectral transmittance T(λ) is
where and are the combined uncertainties for calibrated in-ice irradiance and dark-noise-corrected surface downwelling irradiance, respectively. The combined uncertainty for optical depth τλ is
and the combined uncertainty for katt is
Equation (24) gives a first-order description of σk due to statistical propagation of measurement uncertainties, neglecting higher-order interaction terms. A description of the statistical uncertainty in fitted katt(λ) values is given by the MLE of the regression slope of Eq. 8, which can be expressed in terms of an error model as
where and are the true but unobserved (due to measurement error) optical depth and ice thickness and and are normally distributed error terms. Unlike ordinary least squares, MLE gives an unbiased estimate of the slope and standard error of a linear functional relationship between two variables measured with error (York et al., 2004). The method has been used in similar studies to infer optical coefficients (Zieger et al., 2011). The MLE standard errors for katt(λ) are adjusted for N−2 degrees of freedom with a two-sided t statistic (Cantrell, 2008) and combined in quadrature with systematic uncertainty estimated from Monte Carlo simulation to estimate total combined uncertainty for reported katt(λ) values.
3.1 Spectral transmittance
Four in-ice irradiance spectra were collected at 12, 36, 58, and 77 cm depth below the ice surface on 20 July (hereafter referred to as Layer A) (Fig. 4a) and at 53, 67, 82, and 124 cm on 21 July (hereafter referred to as Layer B). At all depths, spectral transmittance T is maximal at 430 nm and maintains relatively stable and high values at wavelengths up to ∼ 500 nm in the visible, beyond which T decreases into the red end of the visible spectrum following the well-known exponential increase in ice absorptivity (Fig. 4c). Maximum T values vary from 78 % at 12 cm to 45 % at 77 cm. For wavelengths > 500 nm, T rapidly decreases both with wavelength and with depth; beyond ∼ 800 nm nearly all incident light is attenuated within 36 cm of the ice surface, although substantial attenuation is apparent in the 12–36 cm depth region (Fig. 4b). The standard deviation of the 1 Hz raw data is < 1 W m−2 nm−1 at all wavelengths, consistent with field observations of thick cloud cover and diffuse light conditions described in Sect. 2.2. Instrumental noise and high-frequency measurement variations propagate as ±1.6 % uncertainty on T for wavelengths between 400–600 nm, ±1 %–8 % for wavelengths between 350–400 nm, where instrumental noise is higher, and ±1 %–12 % uncertainty for wavelengths between 600–750 nm, where noise is higher and light levels are low.
3.2 Experimental flux attenuation coefficients and albedo
Fitted katt values for Layer A range from 0.98 ± 0.17 to 7.86 ± 0.43 m−1 for wavelengths between 350–750 nm (Fig. 5a), with uncertainty bounds that represent combined statistical and systematic uncertainty (see Sect. 4 for a discussion of systematic error). These values correspond to attenuation lengths of 1.02 ± 0.18 to 0.13 ± 0.007 m, respectively. Layer B katt values are ∼ 12 % lower than Layer A values at 350–500 nm and within 1 % at 650 nm (see inset Fig. 6). For Layer A, the minimum in katt is at 390 nm, blue-shifted relative to the maximum in T at 430 nm. For Layer B, the minimum is at 397 nm. The coefficient of determination (r2) ranged from 0.96–1.0 (p<0.01), with a median value of 0.98, suggesting the data are described appropriately by the Bouguer-law exponential decay model up to ∼ 700 nm, beyond which measured values of in-ice irradiance at 58 and 77 cm depth were too low to reliably fit katt values (see Figs. 4b and 5c). For Layer B values, low light levels prevented fits beyond ∼ 650 nm.
Albedo spectra correspond closely to patterns in transmittance and attenuation (Fig. 5c). The near-UV and blue wavelengths that efficiently transmit ice mostly re-emerge as reflected light, owing to the extremely low values of ice absorptivity in the wavelength range 350–500 nm where albedo is maximum (Gardner and Sharp, 2010; He and Flanner, 2020; Warren et al., 2006). The maximum measured albedo value is 0.81 ± 0.004 at 452 nm, further red-shifted from the minimum in katt and the maximum in T. All three quantities have low variability near the minimum; albedo is 0.79 at 390 nm. The broadband albedo α (Eq. 11) for the 350–900 nm wavelength range is 0.70 ± 0.006, which is high but not atypical for melting white ice under overcast skies (Bøggild et al., 2010).
3.3 Theoretical flux attenuation coefficients
Asymptotic katt values predicted by two-stream theory for optically clean bubbly ice are nearly 1 order of magnitude lower than field estimates for wavelengths < 500 nm, where very small concentrations of LAPs in the measured ice volume dominate absorption (compare dotted grey line to solid blue line, Fig. 6) (Warren et al., 2006). In contrast, field estimates and two-stream theory converge at wavelengths > 540 nm where absorption is dominated by grain-size effects (He et al., 2017; Libois et al., 2013). The magnitude of inferred absorption enhancement in the visible due to LAPs (the quantity βcρi in Eq. 15) varies from 0.009–0.015 m−1 at 350–530 nm. The equivalent black carbon concentration ceqBC inferred at 400 nm is 1–2 ng g−1 for both Layer A and Layer B, where the range covers uncertainty in both the absorption spectra and the absorption Ångström exponent (Doherty et al., 2010). The equivalent mineral dust concentration ceqdust is ∼ 344–620 ng g−1 for Layer A and 303–545 ng g−1 for Layer B. Monte Carlo simulations without detector interference replicate both asymptotic theory for clean bubbly ice (i.e. when forced with ) and field estimates when forced with kabs values inferred from our optical measurements (solid blue squares, Fig. 6). Monte Carlo simulations of detector interference are discussed further in Sect. 4.
3.4 Near-surface attenuation and effective penetration depth
Near the ice surface irradiance is not attenuated exponentially and Bouguer's law does not hold, as indicated by the y intercepts of the straight lines in Fig. 5b at values < 100 %. Effective k′ values (Eq. 17) for the quasi-granular 0–12 cm layer are ∼ 1.5 times higher than katt values for interior bubbly ice at 12–77 cm depth for wavelengths > 570 nm and are up to 4 times higher between 400–570 nm (Fig. 7). Owing to higher near-surface attenuation, transmitted irradiance Iz is overestimated by 10 %–60 % if Bouguer's law is applied to the incident downwelling irradiance I0 using asymptotic katt values, with median overestimation of 23 % (Fig. 8a). In contrast, the piecewise optical depth (Eq. 16) predicts Iz to within 12 % of measured values for all wavelengths between 350–700 nm with a median error of 3 %. Integrated over these wavelengths, χ0 is 0.68 and i0 is 0.66, suggesting 66 % of the total incoming irradiance within this spectral region is absorbed at depths below 10 cm. If katt is used rather than k′ to calculate χ0 and i0, the respective values are 0.81 and 0.79.
Stated in terms of penetration depth, deff varies from 12–84 cm between 350–700 nm. These values are 13 %–44 % lower than attenuation lengths latt inferred from empirical asymptotic katt values. Specifically at 532 nm, deff is 52 cm, or 10 cm lower than the 62 cm empirical latt value and 14 cm lower than the 66 cm theoretical latt value for optically pure bubbly ice. These results point to the potential for reduced optical penetration due to enhanced scattering and absorption on or near the ice surface, as well as within the ice volume where small LAP concentrations reduce optical backscattering due to enhanced absorption.
For smooth ice surfaces, attenuation may be enhanced by refraction at the ice–air interface (Mullen and Warren, 1988). If present, a refractive boundary would enhance near-surface attenuation via external specular reflection, and possibly via enhanced near-surface absorption of the internally reflected downward flux. Following Briegleb and Light (2007; Eq. 20–24), we calculate the external diffuse specular reflectivity for a flat ice surface to be 0.063, meaning specular reflection could enhance attenuation by up to 6.3 %. This value is smaller than the 18 %–44 % near-surface attenuation implied by the y intercepts in Fig. 5b, suggesting specular reflection alone cannot explain the discrepancy. Instead, we suggest that enhanced scattering by the granular near-surface ice microstructure, together with absorptive impurities, enhance near-surface light attenuation at our field site, consistent with observations of the granular and porous surface layer on sea ice (Grenfell and Maykut, 1977; Light et al., 2008).
The effect of random and systematic uncertainties on our optical measurements and fitted katt values is evaluated with Monte Carlo simulation and statistical analysis. We considered systematic uncertainties in detector positions, spectrometer sensitivity to dark current, the non-ideal angular response of the irradiance sensor, and attenuation interference by the PVC detector rod.
The detector positions are known to within 0.9 ± 0.4 cm from independent measurements of the vertical ice thickness Δh. The in situ ice density ρ varied from 801–888 kg m−3 between 4–124 cm where irradiances were measured. The variation in ρ was examined by repeating the analysis with Δz values computed with a single depth-weighted average ρ applied to each Δh value and with ρ values estimated for each Δh value from linear and cubic interpolation of the vertical density profile. The maximum Δz difference was 0.9 cm. The katt values differed by < 1 %, and r2 values were nearly identical. We use the depth-weighted average ρ values to calculate Δz, which is 835 and 855 kg m−3 for the measurements collected on 20 and 21 July, respectively.
Detector position uncertainty was further assessed by fitting katt values with an ensemble of 10 000 Δz values perturbed with random errors drawn from a normal distribution . At all wavelengths, the chance of obtaining a fitted katt value > 2 % from the mean value was < 5 %. We take 2 % as a conservative estimate of systematic uncertainty due to ice thickness measurement bias.
As described in Sect. 3, all irradiance spectra are corrected for residual dark noise. The noise may have varied during the experiment, and dark-noise measurements with the in-ice spectrometer were made on the surface, rather than within the ice. To assess possible bias, we fit katt values with and without residual dark-noise correction. The mean difference was % averaged over the 350–700 nm wavelength range. For a few discrete wavelengths between 350–400 nm and 700–750 nm, differences approached 2 %. These wavelengths are those with the highest dark noise in the reference spectrometer (Fig. 2). At wavelengths between 400–700 nm, differences were < 0.5 %.
Monte Carlo simulations indicate a possible +1 %–14 % systematic bias due to detector interference for Layer A values, and +2 %–8 % for Layer B values (Fig. 9; also see purple stars minus solid squares, Fig. 6). The high end of this range applies to the wavelength region of minimum absorption of ∼ 350–450 nm. The simulated bias is within statistical uncertainty at wavelengths > 450 nm for Layer A and at wavelengths > 400 nm for Layer B (Fig. 9). The non-ideal cosine response of the RCR and the presence of the detector rod both tend to increase katt values relative to the ideal case, as expected given the low albedo of the detector rod. However, detector interference is masked somewhat by the presence of LAPs, as indicated by the larger simulated interference for bubbly ice without LAPs (see dotted grey line and associated Monte Carlo values, Fig. 6). Overall, the combined statistical and systematic uncertainty for the 350–450 nm region is < 20 % for Layer A values and < 14 % for Layer B values and as low as ∼ 5 % for wavelengths > 450 nm (Fig. 9).
5.1 Comparison with attenuation spectra for sea ice, snowpack, and deep glacial ice
We report spectral measurements of near-UV and visible light transmission in bare ablating glacier ice. These measurements are used to estimate irradiance attenuation coefficients katt for the spectral range 350–750 nm. Prior studies quantified katt for sea ice and snowpack (e.g. Fisher et al., 2005; Frey et al., 2011; Gerland et al., 2000; Grenfell and Maykut, 1977; Järvinen and Leppäranta, 2013; King and Simpson, 2001; Light et al., 2008; Meirold-Mautner and Lehning, 2004; Pegau and Zaneveld, 2000; Picard et al., 2016; Tuzet et al., 2019; Warren et al., 2006). Scattering and absorption coefficients were quantified for compressed South Pole glacial ice at 800–2350 m depth by the AMANDA (Antarctic Muon and Neutrino Detector Array) experiment (Ackermann et al., 2006; Askebjer et al., 1995, 1997). For South Pole ice at 800–1000 m depth, visible and near-UV light scatters on air bubbles, below which bubbles transition under pressure to non-scattering clathrates and light scatters on dust grains (Price and Bergström, 1997b). In the bubbly-ice regime studied by AMANDA, ksca values at 532 nm are ∼ 1–3 m−1, comparable to the 1.6 m−1 value quantified in this study. Light scattering in the dusty-ice regime (> 1000 m depth) is not comparable to this study; absorption by dust is discussed below.
Figure 10 compares our katt spectra for glacier ice to seven previously published spectra for snowpack and sea ice. In general, glacier ice is the most transparent structure examined with the exception of multi-year and first-year interior sea ice at wavelengths > 540 nm (Grenfell et al., 2006). Light attenuation in sea ice is controlled by its unique vertical composition including brine inclusions, air pockets, solid salts, sea ice algae, dissolved organic matter, water saturation, and radiative interactions between the ice and underlying ocean (Perovich, 1996). The latter factor, together with differences in optically equivalent grain size, may explain the low attenuation at longer wavelengths for sea ice shown here. Relative to snowpack in Greenland and Antarctica (Järvinen and Leppäranta, 2013; Meirold-Mautner and Lehning, 2004; Warren et al., 2006), attenuation by glacial ice has similar spectral structure but is lower at all wavelengths, reflecting the higher specific surface area of fine-grained polar snow. Attenuation within the surface scattering layer (SSL) of sea ice is intermediate, with spectral structure similar to snowpack and glacial ice. Attenuation at 5 cm depth in snow near Summit, Greenland, is the highest of all, possibly due to direct light scattering in the near-surface optical boundary layer. The comparison demonstrates that katt values vary by nearly 2 orders of magnitude at visible wavelengths due to differences in ice structure and composition.
At visible wavelengths between 350–530 nm our field estimates of katt are up to 1 order of magnitude larger than those obtained from two-stream theory for optically pure bubbly ice, consistent with selective absorption by mineral dust, black carbon, and microorganisms found on glaciers and ice sheet surfaces (Bøggild et al., 2010; Ryan et al., 2018; Stibal et al., 2017; Takeuchi, 2002; Yallop et al., 2012). For context, the absorptivity we document at 400 nm for Layer B can be explained by 1.2–1.8 ng g−1 (parts per billion, ppb) equivalent black carbon concentration. Values for Layer A are 1.4–2.0 ppb. Both estimates are relative to pure-ice absorptivity values reported by Warren et al. (2006). These values are within the range of 2±2 ppb reported for clean snow near DYE-2 on the interior Greenland Ice Sheet considered representative of pre-industrial fallout rates (Doherty et al., 2010). The equivalent mineral dust concentration is ∼ 344–620 ppb for Layer A and 303–545 ppb for Layer B.
Relative to South Pole ice, our absorptivity values broadly agree with AMANDA values within two depth regions corresponding to peaks in atmospheric-dust concentration during the Last Glacial Maximum (LGM) and Marine Isotope Stage 4 (MIS-4) glacial periods ∼ 23 000 and ∼ 66 000 years before present (Fig. 11). For these periods in Earth's history, Southern Hemisphere dust concentrations inferred from the Vostok and Dome C ice cores are ∼ 300–1500 ppb (Muhs, 2013; Petit et al., 1999). Hemispherical dust fluxes are generally synchronous at these timescales; similar peaks at the LGM and MIS-4 are observed in Greenland ice cores (Ruth et al., 2003). However, Northern Hemisphere dust concentrations are several times higher (Muhs, 2013; Ruth et al., 2003), meaning correlation with South Pole absorptivity does not map age at our site. Rather, our optical measurements are consistent with the relatively low dust concentrations during Northern Hemisphere warm periods. For the western Greenland ablation zone, alternating bands of visibly dark and bright outcropping ice are associated with periods of higher and lower aeolian activity during both the Early Holocene (post-LGM) and Late Pleistocene, with a characteristic band of older brighter interglacial ice ∼ 0.7–1 km from the margin where our field site is located (Bøggild et al., 2010; Petrenko et al., 2006; Reeh et al., 2002; Wientjes et al., 2012). Taken together, this suggests the optical properties documented here are representative of Pleistocene interglacial ice with relatively low volumetric LAP concentration and smaller crystal diameters than Holocene ice associated with the “dark zone” further inland (Gow et al., 1997; Petrenko et al., 2006; Wientjes et al., 2011).
Regarding pure-ice absorptivity, our kabs values provide support for the lower-bound pure-ice estimate from Warren et al. (2006) (Fig. 12). The steeply sloping high values in the near UV in the laboratory measurements (Grenfell and Perovich, 1981; Perovich and Govoni, 1991) are now understood as signatures of Rayleigh scattering on nanoscale defects in the laboratory-grown ice (Price and Bergström, 1997a). The South Pole values at 1755 m depth and 830 m depth are contaminated by trace dust deposited during the Late Pleistocene and Early Holocene, respectively (Ackermann et al., 2006) (Fig. 11). The lowest values reported by Warren et al. (2006) (hereafter W06) were obtained by applying Eq. (7) to measurements of transmitted radiance in a single snow layer at ∼ 90–135 cm depth near Dome C in Antarctica contaminated by ∼ 0.6 ppb of black carbon (75∘ S, 123∘ E; 3230 m above sea level). Picard et al. (2016) (hereafter P16) repeated the W06 experiment on 56 transmitted radiance profiles collected in snow near Dome C with variable impurity content. The values shown in Fig. 12 are their best estimate of pure-ice absorptivity from radiance profiles collected in snow with low impurity content (see “clean” subset, Fig. 17 of that paper). P16 were unable to reconcile their values with W06, after considering published values for impurity loadings in the vicinity of Dome C, suggesting the W06 values were unreasonably low. Regardless of that discrepancy, our values were undoubtedly influenced by LAPs but are lower than the P16 values. Treating our Monte Carlo estimate of detector interference as a known systematic error would bring our values closer to the AMANDA and W06 values and further from the P16 values (see dotted line, Fig. 12).
Our absorption minimum is at 390 nm for Layer A values and 397 nm for Layer B values, in agreement with W06 and AMANDA. The wavelength shift in the P16 absorption minimum (430 nm) is coincident with our maximum in transmittance. A similar shift is apparent in all attenuation coefficient spectra shown in Fig. 10 that report the surface as a reference horizon but is absent in those that report an interior reference horizon. P16 used an interior reference horizon, excluding radiance measurements within 8 cm of the surface based on Monte Carlo simulations of detector interference and visual inspection of homogeneous attenuation zones, but their shifted minimum may indicate that radiance profiles in the near-UV and blue spectra were impacted disproportionately by detector rod interference and/or other near-surface effects. The same effects may explain the spectral structure in our near-surface (0–12 cm) effective attenuation coefficient profile between 350–400 nm and its broad minimum between 430–490 nm (dotted line, Fig. 10). Similar spectral structure is apparent in diffuse attenuation coefficients obtained in snowpack in the French Alps using the P16 method (cf. Fig. 3b in Tuzet et al., 2019). Differences aside, our inferred absorption spectrum provides new insight into the magnitude of this fundamental but uncertain optical property and supports a conclusion that the minimum is likely < 10−2 m−1 and possibly lower (Ackermann et al., 2006; Picard et al., 2016; Warren et al., 2006).
In addition to the traditional concept of surface melt, visible light transmission provides an energy source for subsurface heating and internal melting of near-surface glacier ice (Cooper et al., 2018; Hoffman et al., 2014; Liston and Winther, 2005; Schuster, 2001). Prior estimates of subsurface meltwater production in bare ice used two-stream theory forced with values of to calculate katt and the absorbed solar flux as a function of depth below the ice surface in both Greenland and Antarctica (van den Broeke et al., 2008; Hoffman et al., 2014; Kuipers Munneke et al., 2009; Liston and Winther, 2005). The influence of LAPs on subsurface meltwater production has not been quantified to our knowledge and is beyond our scope, but our results suggest LAPs enhance subsurface energy absorption in ablating glacier ice, consistent with enhanced surface melt rates caused by LAPs distributed on bare-ice surfaces and within snowpack (Bøggild et al., 1996; Goelles et al., 2015; Goelles and Bøggild, 2017; Tuzet et al., 2019). From a practical perspective, this suggests that kabs values for contaminated ice given here and snowpack elsewhere (Picard et al., 2016) could provide realistic input for radiative transfer models absent explicit knowledge of realistic LAP concentrations. In contrast, simulations that use the canonical values compiled in Warren and Brandt (2008) will likely underestimate light attenuation and misrepresent the distribution of subsurface absorbed flux unless LAP concentrations are otherwise accounted for.
We report the first in situ spectral measurements of near-UV and visible light transmission in near-surface bare glacial ice, collected at a field site in the western Greenland ablation zone on 20–21 July 2018. In general, our empirical irradiance attenuation coefficients are nearly 1 order of magnitude larger in the range of 350–530 nm than predicted by asymptotic two-stream theory using canonical values for the absorption coefficient of pure ice (Warren and Brandt, 2008). The absorption minimum is 0.013–0.014 ± 0.003 m−1 at 390–397 nm, implying absorption length scales of 69–77 m. The volumetric scattering coefficient is 1.6 ± 0.2 m−1 at 532 nm, with an asymptotic attenuation length scale of 0.62 ± 0.08 m. In addition to light scattering on air bubbles, we find that light attenuation is enhanced by a layer of quasi-granular white ice (weathering crust) that extends from the surface to ∼ 10 cm depth at our field site. The effective penetration depth, which accounts for reduced optical transmission through this granular layer relative to deeper bubbly ice, is 0.52 ± 0.07 m at 532 nm. Our co-located measurements of transmittance and albedo suggest that about 34 % of cloudy sky downwelling solar irradiance at 350–700 nm was absorbed within this upper 10 cm surface layer at this time and location, consistent with observations of the semi-granular surface layer on sea ice. The estimated absorption spectrum suggests equivalent black carbon and mineral dust concentrations consistent with pre-industrial and warm interglacial periods with low Northern Hemisphere aeolian activity, and therefore it may provide a reasonable lower bound on volumetric absorption enhancement due to impurities embedded in outcropping glacial ice in the western Greenland ablation zone.
The Monte Carlo model is publicly available at https://doi.org/10.5281/zenodo.4579073 (Cooper, 2021) and can be used to replicate the results presented here.
The data are hosted by the PANGAEA open-access data repository (https://doi.org/10.1594/PANGAEA.930278, Cooper et al., 2021) and are provided in Appendix A.
The supplement related to this article is available online at: https://doi.org/10.5194/tc-15-1931-2021-supplement.
MGC and LCS designed the experiment. MGC, AKR, MT, RM, SZL, and SEM collected the field data. MGC performed the data analysis and wrote the manuscript. LCS, AKR, MT, SZL, and JVF edited the manuscript.
The authors declare that they have no conflict of interest.
The authors thank Steven G. Warren, Richard E. Brandt, and Ghislain Picard for advice on the experimental design and Robert Leathers and Trijntje Downes for guidance on Monte Carlo radiative transfer methods. The authors thank the handling editor and the two anonymous reviewers for thorough and constructive reviews of the manuscript.
This project was funded by the NASA Cryospheric Sciences Program (grant nos. NNX14AH93G and 80NSSC19K0942) managed by Thomas P. Wagner and Thorsten Markus and by a graduate fellowship from the NASA Earth and Space Science Fellowship Program (grant no. 80NSSC17K0374) managed by Lin Chambers. Polar Field Services and Kangerlussuaq International Science Support (KISS) provided field logistical support.
This paper was edited by Mark Flanner and reviewed by two anonymous referees.
Ackermann, M., Ahrens, J., Bai, X., Bartelt, M., Barwick, S. W., Bay, R. C., Becka, T., Becker, J. K., Becker, K.-H., Berghaus, P., Bernardini, E., Bertrand, D., Boersma, D. J., Böser, S., Botner, O., Bouchta, A., Bouhali, O., Burgess, C., Burgess, T., Castermans, T., Chirkin, D., Collin, B., Conrad, J., Cooley, J., Cowen, D. F., Davour, A., Clercq, C. D., Heros, C. P. de los, Desiati, P., DeYoung, T., Ekström, P., Feser, T., Gaisser, T. K., Ganugapati, R., Geenen, H., Gerhardt, L., Goldschmidt, A., Groß, A., Hallgren, A., Halzen, F., Hanson, K., Hardtke, D. H., Harenberg, T., Hauschildt, T., Helbing, K., Hellwig, M., Herquet, P., Hill, G. C., Hodges, J., Hubert, D., Hughey, B., Hulth, P. O., Hultqvist, K., Hundertmark, S., Jacobsen, J., Kampert, K. H., Karle, A., Kestel, M., Kohnen, G., Köpke, L., Kowalski, M., Kuehn, K., Lang, R., Leich, H., Leuthold, M., Liubarsky, I., Lundberg, J., Madsen, J., Marciniewski, P., Matis, H. S., McParland, C. P., Messarius, T., Minaeva, Y., Miočinović, P., Morse, R., Münich, K., Nahnhauer, R., Nam, J. W., Neunhöffer, T., Niessen, P., Nygren, D. R., Olbrechts, P., Pohl, A. C., Porrata, R., Price, P. B., Przybylski, G. T., Rawlins, K., Resconi, E., Rhode, W., Ribordy, M., Richter, S., Martino, J. R., Sander, H.-G., Schlenstedt, S., Schneider, D., Schwarz, R., Silvestri, A., Solarz, M., Spiczak, G. M., Spiering, C., Stamatikos, M., Steele, D., Steffen, P., Stokstad, R.G., Sulanke, K.-H., Taboada, I., Tarasova, O., Thollander, L., Tilav, S., Wagner, W., Walck, C., Walter, M., Wang, Y.-R., Wiebusch, C.H., Wischnewski, R., Wissing, H., and Woschnagg, K.: Optical properties of deep glacial ice at the South Pole, J. Geophys. Res.-Atmos., 111, D13203, https://doi.org/10.1029/2005JD006687, 2006.
Askebjer, P., Barwick, S. W., Bergström, L., Bouchta, A., Carius, S., Coulthard, A., Engel, K., Erlandsson, B., Goobar, A., Gray, L., Hallgren, A., Halzen, F., Hulth, P. O., Jacobsen, J., Johansson, S., Kandhadai, V., Liubarsky, I., Lowder, D., Miller, T., Mock, P. C., Morse, R., Porrata, R., Price, P. B., Richards, A., Rubinstein, H., Schneider, E., Sun, Q., Tilav, S., Walck, C., and Yodh, G.: Optical Properties of the South Pole Ice at Depths Between 0.8 and 1 Kilometer, Science, 267, 1147–1150, https://doi.org/10.1126/science.267.5201.1147, 1995.
Askebjer, P., Barwick, S. W., Bergström, L., Bouchta, A., Carius, S., Dalberg, E., Erlandsson, B., Goobar, A., Gray, L., Hallgren, A., Halzen, F., Heukenkamp, H., Hulth, P. O., Hundertmark, S., Jacobsen, J., Kandhadai, V., Karle, A., Liubarsky, I., Lowder, D., Miller, T., Mock, P., Morse, R., Porrata, R., Price, P. B., Richards, A., Rubinstein, H., Schneider, E., Spiering, C., Streicher, O., Sun, Q., Thon, T., Tilav, S., Wischnewski, R., Walck, C., and Yodh, G.: UV and optical light transmission properties in deep ice at the South Pole, Geophys. Res. Lett., 24, 1355–1358, https://doi.org/10.1029/97GL01246, 1997.
Bintanja, R. and Van Den Broeke, M. R.: The Surface Energy Balance of Antarctic Snow and Blue Ice, J. Appl. Meteor., 34, 902–926, https://doi.org/10.1175/1520-0450(1995)034<0902:TSEBOA>2.0.CO;2, 1995.
Bøggild, C. E., Oerter, H., and Tukiainen, T.: Increased ablation of Wisconsin ice in eastern north Greenland: observations and modelling, Ann. Glaciol., 23, 144–148, https://doi.org/10.3189/S0260305500013367, 1996.
Bøggild, C. E., Brandt, R. E., Brown, K. J. and Warren, S. G.: The ablation zone in northeast Greenland: ice types, albedos and impurities, J. Glaciol., 56, 101–113, https://doi.org/10.3189/002214310791190776, 2010.
Bohren, C. F.: Colors of snow, frozen waterfalls, and icebergs, J. Opt. Soc. Am., J. Opt. Soc. Am., 73, 1646–1652, https://doi.org/10.1364/JOSA.73.001646, 1983.
Bohren, C. F.: Multiple scattering of light and some of its observable consequences, Am. J. Phys., 55, 524–533, https://doi.org/10.1119/1.15109, 1987.
Bohren, C. F. and Barkstrom, B. R.: Theory of the optical properties of snow, J. Geophys. Res., 79, 4527–4535, https://doi.org/10.1029/JC079i030p04527, 1974.
Brandt, R. E. and Warren, S. G.: Solar-heating rates and temperature profiles in Antarctic snow and ice, J. Glaciol., 39, 99–110, https://doi.org/10.3189/S0022143000015756, 1993.
Briegleb, B. P. and Light, B.: A Delta-Eddington Mutiple Scattering Parameterization for Solar Radiation in the Sea Ice Component of the Community Climate System Model, Technical Note, National Center for Atmospheric Research, Boulder, Colorado, https://doi.org/10.5065/D6B27S71, 2007.
Brunt, K. M., Neumann, T. A., Amundson, J. M., Kavanaugh, J. L., Moussavi, M. S., Walsh, K. M., Cook, W. B., and Markus, T.: MABEL photon-counting laser altimetry data in Alaska for ICESat-2 simulations and development, The Cryosphere, 10, 1707–1719, https://doi.org/10.5194/tc-10-1707-2016, 2016.
Cantrell, C. A.: Technical Note: Review of methods for linear least-squares fitting of data and application to atmospheric chemistry problems, Atmos. Chem. Phys., 8, 5477–5487, https://doi.org/10.5194/acp-8-5477-2008, 2008.
Cooper, M. G.: Ice Monte Carlo Radiative Transfer Model v1.0 (Version v1.0), Zenodo, https://doi.org/10.5281/zenodo.4579073, 2021.
Cooper, M. G., Smith, L. C., Rennermalm, A. K., Miège, C., Pitcher, L. H., Ryan, J. C., Yang, K., and Cooley, S. W.: Meltwater storage in low-density near-surface bare ice in the Greenland ice sheet ablation zone, The Cryosphere, 12, 955–970, https://doi.org/10.5194/tc-12-955-2018, 2018.
Cooper, M. G., Smith, L. C., Rennermalm, A. K., Tedesco, M., Muthyala, R., Leidman, S. Z., Moustafa, S. E., and Fayne, J. V.: Optical attenuation coefficients of glacier ice from 350–700 nm and raw irradiance values from 350–900 nm, PANGAEA, https://doi.org/10.1594/PANGAEA.930278, 2021.
Dadic, R., Mullen, P. C., Schneebeli, M., Brandt, R. E., and Warren, S. G.: Effects of bubbles, cracks, and volcanic tephra on the spectral albedo of bare ice near the Transantarctic Mountains: Implications for sea glaciers on Snowball Earth, J. Geophys. Res.-Earth Surf., 118, 1658–1676, https://doi.org/10.1002/jgrf.20098, 2013.
Deems, J. S., Painter, T. H., and Finnegan, D. C.: Lidar measurement of snow depth: a review, J. Glaciol., 59, 467–479, https://doi.org/10.3189/2013JoG12J154, 2013.
Di Mauro, B., Baccolo, G., Garzonio, R., Giardino, C., Massabò, D., Piazzalunga, A., Rossini, M., and Colombo, R.: Impact of impurities and cryoconite on the optical properties of the Morteratsch Glacier (Swiss Alps), The Cryosphere, 11, 2393–2409, https://doi.org/10.5194/tc-11-2393-2017, 2017.
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, https://doi.org/10.5194/acp-10-11647-2010, 2010.
Fisher, F. N., King, M. D., and Lee-Taylor, J.: Extinction of UV-visible radiation in wet midlatitude (maritime) snow: Implications for increased NOx emission, J. Geophys. Res., 110, D21301, https://doi.org/10.1029/2005JD005963, 2005.
France, J. L., King, M. D., Frey, M. M., Erbland, J., Picard, G., Preunkert, S., MacArthur, A., and Savarino, J.: Snow optical properties at Dome C (Concordia), Antarctica; implications for snow emissions and snow chemistry of reactive nitrogen, Atmos. Chem. Phys., 11, 9787–9801, https://doi.org/10.5194/acp-11-9787-2011, 2011.
Frey, K. E., Perovich, D. K., and Light, B.: The spatial distribution of solar radiation under a melting Arctic sea ice cover, Geophys. Res. Lett., 38, L22501, https://doi.org/10.1029/2011GL049421, 2011.
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., 115, F01009, https://doi.org/10.1029/2009JF001444, 2010.
Gardner, A. S., Smith, B. E., Brunt, K. M., Harding, D. J., Neumann, T., and Walsh, K.: ICESat2 subsurface-scattering biases estimated based on the 2015 SIMPL/AVRIS campaign, in AGU Fall Meeting Abstracts, vol. 41, C41C-0710, http://adsabs.harvard.edu/abs/2015AGUFM.C41C0710G (last access: 25 January 2019), 2015.
Gerland, S., Liston, G. E., Winther, J.-G., Ørbæk, J. B., and Ivanov, B. V.: Attenuation of solar radiation in Arctic snow: field observations and modelling, Ann. Glaciol., 31, 364–368, https://doi.org/10.3189/172756400781820444, 2000.
Goelles, T. and Bøggild, C. E.: Albedo reduction of ice caused by dust and black carbon accumulation: a model applied to the K-transect, West Greenland, J. Glaciol., 63, 1063–1076, https://doi.org/10.1017/jog.2017.74, 2017.
Goelles, T., Bøggild, C. E., and Greve, R.: Ice sheet mass loss caused by dust and black carbon accumulation, The Cryosphere, 9, 1845–1856, https://doi.org/10.5194/tc-9-1845-2015, 2015.
Gow, A. J., Meese, D. A., Alley, R. B., Fitzpatrick, J. J., Anandakrishnan, S., Woods, G. A., and Elder, B. C.: Physical and structural properties of the Greenland Ice Sheet Project 2 ice core: A review, J. Geophys. Res., 102, 26559–26575, https://doi.org/10.1029/97JC00165, 1997.
Greeley, A., Kurtz, N. T., Neumann, T., and Markus, T.: Estimating Surface Elevation Bias Due to Subsurface Scattered Photons from Visible Wavelength Laser Altimeters, in AGU Fall Meeting Abstracts, vol. 51, http://adsabs.harvard.edu/abs/2017AGUFM.C51A0961G (last access: 25 January 2019), 2017.
Grenfell, T. C.: The Effects of Ice Thickness on the Exchange of Solar Radiation Over the Polar Oceans, J. Glaciol., 22, 305–320, https://doi.org/10.3189/S0022143000014295, 1979.
Grenfell, T. C.: A radiative transfer model for sea ice with vertical structure variations, J. Geophys. Res.-Oceans, 96, 16991–17001, https://doi.org/10.1029/91JC01595, 1991.
Grenfell, T. C. and Maykut, G. A.: The Optical Properties of Ice and Snow in the Arctic Basin, J. Glaciol., 18, 445–463, https://doi.org/10.3189/S0022143000021122, 1977.
Grenfell, T. C. and Perovich, D. K.: Radiation absorption coefficients of polycrystalline ice from 400–1400 nm, J. Geophys. Res., 86, 7447–7450, https://doi.org/10.1029/JC086iC08p07447, 1981.
Grenfell, T. C. and Warren, S. G.: Representation of a nonspherical ice particle by a collection of independent spheres for scattering and absorption of radiation, J. Geophys. Res., 104, 31697–31709, https://doi.org/10.1029/1999JD900496, 1999.
Grenfell, T. C., Light, B., and Perovich, D. K.: Spectral transmission and implications for the partitioning of shortwave radiation in arctic sea ice, Ann. Glaciol., 44, 1–6, https://doi.org/10.3189/172756406781811763, 2006.
He, C. and Flanner, M.: Snow Albedo and Radiative Transfer: Theory, Modeling, and Parameterization, in: Springer Series in Light Scattering, edited by A. Kokhanovsky, Springer International Publishing, Cham, 67–133, https://doi.org/10.1007/978-3-030-38696-2_3, 2020.
He, C., Takano, Y., Liou, K.-N., Yang, P., Li, Q., and Chen, F.: Impact of Snow Grain Shape and Black Carbon–Snow Internal Mixing on Snow Optical Properties: Parameterizations for Climate Models, J. Climate, 30, 10019–10036, https://doi.org/10.1175/JCLI-D-17-0300.1, 2017.
Hoffman, M. J., Fountain, A. G., and Liston, G. E.: Near-surface internal melting: a substantial mass loss on Antarctic Dry Valley glaciers, J. Glaciol., 60, 361–374, https://doi.org/10.3189/2014JoG13J095, 2014.
Holland, M. M., Bailey, D. A., Briegleb, B. P., Light, B., and Hunke, E.: Improved Sea Ice Shortwave Radiation Physics in CCSM4: The Impact of Melt Ponds and Aerosols on Arctic Sea Ice, J. Climate, 25, 1413–1430, https://doi.org/10.1175/JCLI-D-11-00078.1, 2012.
Järvinen, O. and Leppäranta, M.: Solar radiation transfer in the surface snow layer in Dronning Maud Land, Antarctica, Polar Science, 7, 1–17, https://doi.org/10.1016/j.polar.2013.03.002, 2013.
Joseph, J. H., Wiscombe, W. J., and Weinman, J. A.: The Delta-Eddington Approximation for Radiative Flux Transfer, J. Atmos. Sci., 33, 2452–2459, https://doi.org/10.1175/1520-0469(1976)033<2452:TDEAFR>2.0.CO;2, 1976.
King, M. D. and Simpson, W. R.: Extinction of UV radiation in Arctic snow at Alert, Canada (82∘ N), J. Geophys. Res., 106, 12499–12507, https://doi.org/10.1029/2001JD900006, 2001.
Kokhanovsky, A. A. and Zege, E. P.: Scattering optics of snow, Appl. Opt., 43, 1589, https://doi.org/10.1364/AO.43.001589, 2004.
Kuipers Munneke, P., van den Broeke, M. R., Reijmer, C. H., Helsen, M. M., Boot, W., Schneebeli, M., and Steffen, K.: The role of radiation penetration in the energy budget of the snowpack at Summit, Greenland, The Cryosphere, 3, 155–165, https://doi.org/10.5194/tc-3-155-2009, 2009.
Leathers, R. A., Downes, T. V., Davis, C. O., and Mobley, C. D.: Monte Carlo Radiative Transfer Simulations for Ocean Optics: A Practical Guide, Memorandum, Naval Research Laboratory, Washington, DC, available at: https://www.oceanopticsbook.info/packages/iws_l2h/conversion/files/Leathersetal_NRL2004.pdf (last access: 11 October 2020), 2004.
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, https://doi.org/10.5194/tc-7-1803-2013, 2013.
Light, B., Maykut, G. A., and Grenfell, T. C.: A two-dimensional Monte Carlo model of radiative transfer in sea ice, J. Geophys. Res.-Oceans, 108, 3219, https://doi.org/10.1029/2002JC001513, 2003.
Light, B., Maykut, G. A., and Grenfell, T. C.: A temperature-dependent, structural-optical model of first-year sea ice, J. Geophys. Res., 109, C06013, https://doi.org/10.1029/2003JC002164, 2004.
Light, B., Grenfell, T. C., and Perovich, D. K.: Transmission and absorption of solar radiation by Arctic sea ice during the melt season, J. Geophys. Res., 113, C03023, https://doi.org/10.1029/2006JC003977, 2008.
Liston, G. E. and Winther, J.-G.: Antarctic Surface and Subsurface Snow and Ice Melt Fluxes, J. Climate, 18, 1469–1481, https://doi.org/10.1175/JCLI3344.1, 2005.
Liston, G. E., Bruland, O., Elvehøy, H., and Sand, K.: Below-surface ice melt on the coastal Antarctic ice sheet, J. Glaciology, 45, 273–285, https://doi.org/10.3189/002214399793377130, 1999a.
Liston, G. E., Bruland, O., Winther, J.-G., Elvehøy, H., and Sand, K.: Meltwater production in Antarctic blue-ice areas: sensitivity to changes in atmospheric forcing, Polar Res., 18, 283–290, https://doi.org/10.1111/j.1751-8369.1999.tb00305.x, 1999b.
Malinka, A., Zege, E., Heygster, G., and Istomina, L.: Reflective properties of white sea ice and snow, The Cryosphere, 10, 2541–2557, https://doi.org/10.5194/tc-10-2541-2016, 2016.
Malinka, A. V.: Light scattering in porous materials: Geometrical optics and stereological approach, J. Quant. Spectrosc. Ra., 141, 14–23, https://doi.org/10.1016/j.jqsrt.2014.02.022, 2014.
Markus, T., Neumann, T., Martino, A., Abdalati, W., Brunt, K., Csatho, B., Farrell, S., Fricker, H., Gardner, A., Harding, D., Jasinski, M., Kwok, R., Magruder, L., Lubin, D., Luthcke, S., Morison, J., Nelson, R., Neuenschwander, A., Palm, S., Popescu, S., Shum, C., Schutz, B. E., Smith, B., Yang, Y., and Zwally, J.: The Ice, Cloud, and land Elevation Satellite-2 (ICESat-2): Science requirements, concept, and implementation, Remote Sens. Environ., 190, 260–273, https://doi.org/10.1016/j.rse.2016.12.029, 2017.
Mätzler, C.: MATLAB Functions for Mie Scattering and Absorption, Version 2, Research Report, Institut für Angewandte Physik, Bern, Switzerland, https://doi.org/10.7892/boris.146550, 2002.
Maykut, G. A. and Untersteiner, N.: Some results from a time-dependent thermodynamic model of sea ice, J. Geophys. Res., 76, 1550–1575, https://doi.org/10.1029/JC076i006p01550, 1971.
Meirold-Mautner, I. and Lehning, M.: Measurements and model calculations of the solar shortwave fluxes in snow on Summit, Greenland, Ann. Glaciol., 38, 279–284, https://doi.org/10.3189/172756404781814753, 2004.
Muhs, D. R.: The geologic records of dust in the Quaternary, Aeolian Res., 9, 3–48, https://doi.org/10.1016/j.aeolia.2012.08.001, 2013.
Mullen, P. C. and Warren, S. G.: Theory of the optical properties of lake ice, J. Geophys. Res., 93, 8403–8414, https://doi.org/10.1029/JD093iD07p08403, 1988.
Pegau, W. S. and Zaneveld, J. R. V.: Field measurements of in-ice radiance, Cold Reg. Sci. Technol., 31, 33–46, https://doi.org/10.1016/S0165-232X(00)00004-5, 2000.
Perovich, D. K.: The Optical Properties of Sea Ice., U.S. Army Cold Regions Research and Engineering Laboratory, Hanover, NH, https://apps.dtic.mil/dtic/tr/fulltext/u2/a310586.pdf, 1996.
Perovich, D. K. and Govoni, J. W.: Absorption coefficients of ice from 250 to 400 nm, Geophys. Res. Lett., 18, 1233–1235, https://doi.org/10.1029/91GL01642, 1991.
Petit, J. R., Jouzel, J., Raynaud, D., Barkov, N. I., Barnola, J.-M., Basile, I., Bender, M., Chappellaz, J., Davis, M., Delaygue, G., Delmotte, M., Kotlyakov, V. M., Legrand, M., Lipenkov, V. Y., Lorius, C., PÉpin, L., Ritz, C., Saltzman, E., and Stievenard, M.: Climate and atmospheric history of the past 420,000 years from the Vostok ice core, Antarctica, Nature, 399, 429–436, https://doi.org/10.1038/20859, 1999.
Petrenko, V. V., Severinghaus, J. P., Brook, E. J., Reeh, N., and Schaefer, H.: Gas records from the West Greenland ice margin covering the Last Glacial Termination: a horizontal ice core, Quatern. Sci. Rev., 25, 865–875, https://doi.org/10.1016/j.quascirev.2005.09.005, 2006.
Picard, G., Libois, Q., and Arnaud, L.: Refinement of the ice absorption spectrum in the visible using radiance profile measurements in Antarctic snow, The Cryosphere, 10, 2655–2672, https://doi.org/10.5194/tc-10-2655-2016, 2016.
Price, P. B. and Bergström, L.: Enhanced Rayleigh scattering as a signature of nanoscale defects in highly transparent solids, Philos. Mag. A, 75, 1383–1390, https://doi.org/10.1080/01418619708209861, 1997a.
Price, P. B. and Bergström, L.: Optical properties of deep ice at the South Pole: scattering, Appl. Opt., 36, 4181–4194, https://doi.org/10.1364/AO.36.004181, 1997b.
Reeh, N., Oerter, H., and Thomsen, H. H.: Comparison between Greenland ice-margin and ice-core oxygen-18 records, Ann. Glaciol., 35, 136–144, https://doi.org/10.3189/172756402781817365, 2002.
Ridley, J. K. and Partington, K. C.: A model of satellite radar altimeter return from ice sheets, International J. Remote Sens., 9, 601–624, https://doi.org/10.1080/01431168808954881, 1988.
Rignot, E., Echelmeyer, K., and Krabill, W.: Penetration depth of interferometric synthetic-aperture radar signals in snow and ice, Geophys. Res. Lett., 28, 3501–3504, https://doi.org/10.1029/2000GL012484, 2001.
Ruth, U., Wagenbach, D., Steffensen, J. P., and Bigler, M.: Continuous record of microparticle concentration and size distribution in the central Greenland NGRIP ice core during the last glacial period, J. Geophys. Res.-Atmos., 108, 4098, https://doi.org/10.1029/2002JD002376, 2003.
Ryan, J. C., Hubbard, A. L., Stibal, M., Irvine-Fynn, T. D., Cook, J., Smith, L. C., Cameron, K., and Box, J. E.: Dark zone of the Greenland Ice Sheet controlled by distributed biologically-active impurities, Nat. Commun., 9, 1065, https://doi.org/10.1038/s41467-018-03353-2, 2018.
Schuster, A.: Radiation through a foggy atmosphere, Astrophys. J., 21, 1–22, 1905.
Schuster, C.: Weathering crust processes on melting glacier ice (Alberta, Canada), Theses and Dissertations (Comprehensive), no. 489, available at: http://scholars.wlu.ca/etd/489 (last access: 3 December 2016), 2001.
Schutz, B. E., Zwally, H. J., Shuman, C. A., Hancock, D., and DiMarzio, J. P.: Overview of the ICESat Mission, Geophys. Res. Lett., 32, L21S01, https://doi.org/10.1029/2005GL024009, 2005.
Smith, B. E., Gardner, A., Schneider, A., and Flanner, M.: Modeling biases in laser-altimetry measurements caused by scattering of green light in snow, Remote Sens. Environ., 215, 398–410, https://doi.org/10.1016/j.rse.2018.06.012, 2018.
Stibal, M., Box, J. E., Cameron, K. A., Langen, P. L., Yallop, M. L., Mottram, R. H., Khan, A. L., Molotch, N. P., Chrismas, N. A. M., Quaglia, F. C., Remias, D., Smeets, C. J. P. P., Broeke, M. R. van den, Ryan, J. C., Hubbard, A., Tranter, M., As, D., van and Ahlstrøm, A. P.: Algae Drive Enhanced Darkening of Bare Ice on the Greenland Ice Sheet, Geophys. Res. Lett., 44, 11463–11471, https://doi.org/10.1002/2017GL075958, 2017.
Takeuchi, N.: Optical characteristics of cryoconite (surface dust) on glaciers: the relationship between light absorbency and the property of organic matter contained in the cryoconite, Ann. Glaciol., 34, 409–414, https://doi.org/10.3189/172756402781817743, 2002.
Taylor, B. N. and Kuyatt, C. E.: Guidelines for Evaluating and Expressing the Uncertainty of NIST Measurement Results (NIST Technical Note vol. 1297), National Institue of Standards and Technology, Gaithersburg, MD, http://physics.nist.gov/Pubs/guidelines/TN1297/tn1297s.pdf (last access: 18 January 2021), 1994.
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, https://doi.org/10.5194/tc-13-2169-2019, 2019.
van de Hulst, H. C.: Multiple light scattering: tables, formulas, and applications, Academic Press, New York, 1980.
van den Broeke, M., Smeets, P., Ettema, J., van der Veen, C., van de Wal, R., and Oerlemans, J.: Partitioning of melt energy and meltwater fluxes in the ablation zone of the west Greenland ice sheet, The Cryosphere, 2, 179–189, https://doi.org/10.5194/tc-2-179-2008, 2008.
Wang, L., Jacques, S. L., and Zheng, L.: MCML – Monte Carlo modeling of light transport in multi-layered tissues, Comput. Meth. Prog. Bio., 47, 131–146, https://doi.org/10.1016/0169-2607(95)01640-F, 1995.
Warren, S. G.: Optical properties of snow, Rev. Geophys., 20, 67–89, https://doi.org/10.1029/RG020i001p00067, 1982.
Warren, S. G.: Optical constants of ice from the ultraviolet to the microwave, Appl. Opt., 23, 1206–1225, https://doi.org/10.1364/AO.23.001206, 1984.
Warren, S. G. and Brandt, R. E.: Optical constants of ice from the ultraviolet to the microwave: A revised compilation, J. Geophys. Res., 113, D14220, https://doi.org/10.1029/2007JD009744, 2008.
Warren, S. G., Brandt, R. E., Grenfell, T. C., and McKay, C. P.: Snowball Earth: Ice thickness on the tropical ocean, J. Geophys. Res.-Oceans, 107, 31-1–31-18, https://doi.org/10.1029/2001JC001123, 2002.
Warren, S. G., Brandt, R. E., and Grenfell, T. C.: Visible and near-ultraviolet absorption spectrum of ice from transmission of solar radiation into snow, Appl. Opt., 45, 5320–5334, https://doi.org/10.1364/AO.45.005320, 2006.
Wientjes, I. G. M., Van de Wal, R. S. W., Reichart, G. J., Sluijs, A., and Oerlemans, J.: Dust from the dark region in the western ablation zone of the Greenland ice sheet, The Cryosphere, 5, 589–601, https://doi.org/10.5194/tc-5-589-2011, 2011.
Wientjes, I. G. M., De Van Wal, R. S. W., Schwikowski, M., Zapf, A., Fahrni, S., and Wacker, L.: Carbonaceous particles reveal that Late Holocene dust causes the dark region in the western ablation zone of the Greenland ice sheet, J. Glaciol., 58, 787–794, https://doi.org/10.3189/2012JoG11J165, 2012.
Wiscombe, W. J. and Warren, S. G.: A Model for the Spectral Albedo of Snow. I: Pure Snow, J. Atmos. Sci., 37, 2712–2733, https://doi.org/10.1175/1520-0469(1980)037<2712:AMFTSA>2.0.CO;2, 1980.
Woschnagg, K. and Price, P. B.: Temperature dependence of absorption in ice at 532 nm, Appl. Opt., 40, 2496–2500, https://doi.org/10.1364/AO.40.002496, 2001.
Yallop, M. L., Anesio, A. M., Perkins, R. G., Cook, J., Telling, J., Fagan, D., MacFarlane, J., Stibal, M., Barker, G., Bellas, C., Hodson, A., Tranter, M., Wadham, J., and Roberts, N. W.: Photophysiology and albedo-changing potential of the ice algal community on the surface of the Greenland ice sheet, ISME J., 6, 2302–2313, https://doi.org/10.1038/ismej.2012.107, 2012.
Yang, Y., Marshak, A., Han, M., Palm, S. P., and Harding, D. J.: Snow grain size retrieval over the polar ice sheets with the Ice, Cloud, and land Elevation Satellite (ICESat) observations, J. Quant. Spectrosc. Ra., 188, 159–164, https://doi.org/10.1016/j.jqsrt.2016.03.033, 2017.
York, D., Evensen, N. M., Martınez, M. L., and De Basabe Delgado, J.: Unified equations for the slope, intercept, and standard errors of the best straight line, Am. J. Phys., 72, 367–375, https://doi.org/10.1119/1.1632486, 2004.
Zebker, H. A. and Weber Hoen, E.: Penetration depths inferred from interferometric volume decorrelation observed over the Greenland Ice Sheet, IEEE T. Geosci. Remote Sens., 38, 2571–2583, https://doi.org/10.1109/36.885204, 2000.
Zhang, X., Qiu, J., Li, X., Zhao, J., and Liu, L.: Complex refractive indices measurements of polymers in visible and near-infrared bands, Appl. Opt., 59, 2337, https://doi.org/10.1364/AO.383831, 2020.
Zieger, P., Weingartner, E., Henzing, J., Moerman, M., de Leeuw, G., Mikkilä, J., Ehn, M., Petäjä, T., Clémer, K., van Roozendael, M., Yilmaz, S., Frieß, U., Irie, H., Wagner, T., Shaiganfar, R., Beirle, S., Apituley, A., Wilson, K., and Baltensperger, U.: Comparison of ambient aerosol extinction coefficients obtained from in-situ, MAX-DOAS and LIDAR measurements at Cabauw, Atmos. Chem. Phys., 11, 2603–2624, https://doi.org/10.5194/acp-11-2603-2011, 2011.