Seasonal snow covers Arctic lands 6 to 10 months of the
year and is therefore an essential element of the Arctic geosphere and
biosphere. Yet, even the most sophisticated snow physics models are not able
to simulate fundamental physical properties of Arctic snowpacks such as
density, thermal conductivity and specific surface area. The development of
improved snow models is in progress, but testing requires detailed driving
and validation data for high Arctic herb tundra sites, which are presently
not available. We present 6 years of such data for an ice-wedge polygonal
site in the Canadian high Arctic, in Qarlikturvik valley on Bylot Island at
73.15∘ N. The site is on herb tundra with no erect vegetation and
thick permafrost. Detailed soil properties are provided. Driving data are
comprised of air temperature, air relative and specific humidity, wind
speed, shortwave and longwave downwelling radiation, atmospheric pressure,
and precipitation. Validation data include time series of snow depth,
shortwave and longwave upwelling radiation, surface temperature, snow
temperature profiles, soil temperature and water content profiles at five
depths, snow thermal conductivity at three heights, and soil thermal
conductivity at 10 cm depth. Field campaigns in mid-May for 5 of the 6 years
of interest provided spatially averaged snow depths and vertical profiles of
snow density and specific surface area in the polygon of interest and at
other spots in the valley. Data are available at
10.5885/45693CE-02685A5200DD4C38 (Domine et al., 2021). Data
files will be updated as more years of data become available.
Introduction
The seasonal snowpack covers high-latitude regions at low elevation 6 to
10 months of the year (Connolly et al., 2019). Snow physical properties
such as specific surface area (SSA) and density determine albedo (Carmagnola
et al., 2013), an essential component of the surface energy budget. Snow
thermal conductivity determines heat exchanges between the atmosphere and
the ground and therefore impacts the permafrost thermal regime (Zhang,
2005; Barrere et al., 2017). Where permafrost is absent, snow thermal
conductivity determines whether and when the ground freezes, with a very
strong impact on nutrient recycling and the accumulation of organic
compounds in soils (Saccone et al., 2013; Myers-Smith and Hik,
2013; Buckeridge and Grogan, 2008). Snow thermal conductivity also strongly
influences surface air temperature (Domine et al., 2019), and inadequate
simulations of this variable can modify simulated air temperature by up to
4 ∘C, potentially affecting meteorological models and high-latitude weather forecasts. Lastly, snow physical properties affect
wildlife. For example, lemmings, small rodents of the high Arctic, live,
move, feed and reproduce under the snow 9 months of the year in the high
Arctic (Poirier et al., 2019; Bilodeau et al., 2013). Their population cycles
have intrigued scientists for decades (Fauteux et al., 2015), and recent
studies have indicated that snow physical properties, and in particular the
hardness of the snow basal layer, may strongly impact lemming reproductive
success in winter and their summer population dynamics (Domine et al.,
2018b). Likewise, larger Arctic herbivores such as caribou are strongly
affected by snow physical properties, which determine their access to food
under the snow (Langlois et al., 2017).
Adequately simulating snow physical properties is therefore essential for
understanding and/or projecting climate, meteorology, the state of
permafrost, nutrient cycling and carbon storage in soils and therefore
vegetation growth and the carbon budget of Arctic soils and permafrost as well as and
wildlife ecology and population dynamics. Despite these critically high
stakes, there is today no detailed snow physics model capable of simulating
Arctic snowpack physical properties adequately. Domine et al. (2016b) have
shown that both detailed snow models Crocus (Vionnet et al., 2012) and
SNOWPACK (Bartelt and Lehning, 2002) failed to simulate essential
characteristics of the snowpack at the high Arctic site of Bylot Island
(73∘ N). In particular, simulated snow density profiles were
inverted relative to observations. Both models predicted dense basal layers
and light top layers, while most snow observations in the high Arctic have
reported low-density basal layers made of depth hoar and high-density upper
wind slabs (Derksen et al., 2009; Domine et al.,
2002, 2012, 2016b; Gouttevin et al., 2018). The explanation proposed
(Domine et al., 2019, 2016b) is that Crocus and SNOWPACK were
calibrated using data from sites in the Alps, i.e., for mid-latitude warm, thick snowpacks, while the Arctic features cold, thin snowpacks. In the Alps,
an essential driving variable in snowpack vertical profiles of physical
properties is compaction by the snow overburden. In the thin Arctic
snowpack, this process is negligible, and the main process involved in
determining the evolution of the density profile after precipitation and
wind compaction is the upward flux of water vapor, caused mostly by
convection within the snowpack (Trabant and Benson, 1972; Sturm and Benson,
1997; Sturm and Johnson, 1991; Domine et al., 2016b). This flux is driven by
the large vertical temperature gradient, which redistributes mass from lower
to upper layers. This process is so intense and the associated mass loss so
large that it sometimes leads to the collapse of the basal depth hoar layer
(Domine et al., 2016b), even in the low Arctic (Domine et al., 2015). This
process is not simulated by Crocus or SNOWPACK, leading to erroneous
outputs.
Upward water vapor fluxes are also the main determinant of snowpack vertical
profiles of physical properties in many areas of the boreal forest (Sturm
and Benson, 1997). Since Arctic and boreal forest snowpacks together
represent by far the most important seasonal snowpacks on Earth on an areal
basis (Sturm et al., 1995), it is essential that data sets be available
which allow the testing of snow models and their application to high
latitudes. At present, there is not to our knowledge any multi-year high
Arctic data set complete enough to drive and validate in detail a snow physics
model. The northernmost site used in the latest snow model intercomparison
project (SNOWMIP; Krinner et al., 2018) is Sodankylä, Finland,
67∘ N. Although it is classified as “Arctic” in Krinner et al. (2018), Sodankylä is in the boreal forest, while Arctic usually refers to
regions above the tree line. In the boreal forest, the dense wind slabs observed
in the Arctic do not form, and snow properties are markedly different from
those on Arctic tundra (Sturm et al., 1995). Boike et al. (2018) have
provided a 20-year data set of permafrost, active layer and meteorological
data for a site near Ny-Ålesund, Svalbard (78.5∘ N, 11.6∘ E), suitable for
driving land surface and snow models. However, while this data set can be
used for numerous valuable applications, the snow validating data are
limited to snow depth and to snow pit observations in late April or early
May. The snow physical data are comprised of density at several heights and
of the vertical temperature profile when the pit was dug. These data are
useful but are probably not sufficient for thoroughly testing snow physics
model performance under Arctic conditions. Boike et al. (2019) also
presented a 16-year data set of permafrost, active layer and meteorological
conditions for Samoylov Island in Siberia (72.3∘ N, 126.5∘ E). The site is in
the Lena river delta and features ice-wedge polygons, with a very high
fraction of ground ice. The permafrost data, together with a previous paper
(Boike et al., 2013), are extremely detailed so that this data set is
certainly particularly useful for permafrost simulations. Regarding snow
however, data are more limited and comprised of snow depth, time lapse
photographs, and some spring measurements of snow properties such as density
and thermal conductivity. Snow precipitation has not been measured there.
The Samoylov site has been used to test the SNOWPACK model. Gouttevin et
al. (2018) used a 1-year driving data set to simulate snow and used snow
pit data from a field campaign in April as well as ground temperature
monitoring at 5 cm depth as validating data. They modified the SNOWPACK
model to adapt it to Arctic conditions and in particular modified
grain growth rate laws. However, they did not treat upward water vapor fluxes
explicitly. That study constitutes valuable progress towards the elaboration
of an Arctic snow model, but a 1-year test is not sufficient to oversee
the variety of high Arctic conditions. For example, in their study, they
encountered high-density “indurated” depth hoar which is frequent but far
from ubiquitous in the high Arctic. The motivation of the present work is
therefore to provide over an extended period driving and testing data for
snow physics at a high Arctic site so that the ability of snow physics and
land surface schemes can be tested in a variety of meteorological situations
in these high Arctic conditions.
We provide standard meteorological data for driving models: air temperature
and relative humidity, atmospheric pressure, wind speed, shortwave
(SW↓) and longwave (LW↓) incoming radiation, and
precipitation. Detailed soil properties such as density, granulometry,
organic carbon content and thermal conductivity at several depths are also
provided. For validation, we provide continuous monitoring of snow depth,
upwelling shortwave radiation (SW↑), surface albedo, upwelling
longwave radiation (LW↑), surface temperature, snow temperature at
five heights, snow thermal conductivity at three heights, soil temperature
and volume liquid water content at five depths, and soil thermal conductivity
at 10 cm depth. The data cover a 6-year span from 11 July 2013 to 25 June
2019. Data from 2019–2020 could not yet be retrieved due to the COVID-19
pandemic, which prevented access to our site. However, data from future
years will be added to the set as they become available. Furthermore, field
campaigns at this site were possible in May 2014, 2015, 2017, 2018 and 2019,
and we also present snow pit observations of snow stratigraphy and
measurements of vertical profiles of snow density and SSA for those years.
In May 2016, logistical difficulties prevented access to the site.
Instruments used to obtain meteorological, snow and soil data at
Bylot Island.
VariableInstrumentManufacturerInstrument height or depthCommentShortwave radiationdown- and upwellingCNR4 pyranometers, 300–3000 nm with CNF4 ventilator and heaterKipp & Zonen2.3 mCNF4 on for 5 min before hourly measurementLongwave radiationdown- and upwellingCNR4 pyrgeometers, 4.5–42 µm with CNF4 ventilator and heaterKipp & Zonen2.3 mCNF4 on for 5 min before hourly measurementSnow (winter) or soil(summer) surface temperatureIR 120 infrared sensor,8–14 µmCampbell Scientific1.5 mMeasurement everyminute, hourly average recordedAir temperature andrelative humidity(relative to liquid orsupercooled water)HC2-S3-XT sensor, insidewhite PVC tubing, ventilatedRotronic2.3 mVentilator on for 5 minbefore hourly measurementWind speedCup anemometerVector instruments2.3 mMeasurement everyminute, hourly average recordedPrecipitationGeonor 200, complemented with data fromGeonor gauges at PondInlet and Cape Liverpool1.5 mAtmospheric pressureNot measured, data fromPond Inlet and Cape Liverpool usedAverage value of datafrom station at Pond Inlet and Cape LiverpoolSnow depthSR50A acoustic gaugeCampbell scientific2.2 mMeasurement everyminute, hourly average recordedSnow thermalconductivityTP08 heated needle probesHukseflux7, 17, 27 cm in 2013,changed to 2, 12, 22 cm in July 2014Measurement everyother day at 05:00 local summer time (UTC-4)Snow temperaturePt 100 thermistorsHome-assembled sensors2, 7, 17, 27, 37 cm,changed to 0, 5, 15, 25,35 in July 2018Measurement everyminute, hourly average recordedSoil thermalconductivityTP08 heated needle probesHukseflux10 cmMeasurement everyother day at 05:00 local summer timeSoil temperature andvolume water content5TM sensorsDecagon(now Meter)2, 5, 10, 15, 21 cmHourly measurementSceneryTime lapse cameraReconyx1.5 mFour pictures a day
(Left) map of Qarlikturvik valley on Bylot Island in the Canadian
Arctic Archipelago, showing the TUNDRA study site. The Environment and Climate Change Canada (ECCC) meteorological
stations at Cape Liverpool and Pond Inlet as well as the SILA, CAMP and
SALIX stations, all of which were used for data gap filling, are also shown.
View of the polygon where instruments were installed (right). General aerial
view of the polygon area, with arrow pointed at TUNDRA station. Maps from
http://atlas.gc.ca (last access: 15 January 2021).
Study site and instruments
Our study site is in Qarlikturvik valley on Bylot Island, north of Baffin
Island in the Canadian Archipelago (Fig. 1). The nearest community is Pond
Inlet on Baffin Island, 85 km to the southeast, from which our permanent
camp can be accessed by helicopter in summer or snowmobile in spring.
Aircraft landing on skis is occasionally possible in spring or on a nearby
beach in summer with tundra tires. Atmospheric monitoring was initiated in
August 1993 at the Camp Lake site (hereafter CAMP; 73.1567∘ N,
79.9571∘ W), and air temperature and humidity as well as wind
speed (3 m) data are available since that date (CEN, 2020). In July 2004 a
10 m tower (called SILA) was built (at 73.1522∘ N,
79.9886∘ W) 1150 m west-southwest of CAMP and equipped to measure wind
speed and direction (10 m) and air temperature. Data are also available
(CEN, 2020) at the same repository and DOI as the CAMP data. The data
discussed here are from a comprehensive monitoring station established on 7 July 2013 at the TUNDRA site (73.1504∘ N, 80.0046∘ W),
1700 m to the west-southwest of CAMP. The GPS elevation of the site is 20 to 25 m, but
according to the Canada Atlas maps (http://atlas.gc.ca, last access: 15 January 2021), the site is just below the
20 m contour line. Google Earth indicates an elevation of 25 m. The site has
been presented by Domine et al. (2016b). Briefly, the instruments are
within a rather well-drained low-center ice-wedge polygon typical of
permafrost landscape. The polygon is about 11 m in its largest dimension, and
all instruments are within 3 m of its center. Equipment is detailed in Table 1 and includes a tripod supporting meteorological instruments at 2.3 m
height, a vertical polyethylene post supporting three TP08 heated needle
probes for snow thermal conductivity measurements and another post
supporting five thermistors for snow temperature measurements. Below the
surface, 5TM sensors from Decagon (now Meter) measure soil temperature and
volumetric liquid water content at five depths, and a TP08 probe monitors
soil thermal conductivity. The instruments are accessed and maintained once
a year in summer. Instrument failure thus cannot always be fixed
immediately. The CNR4 radiometer was brought back south in 2019 for
recalibration by Kipp & Zonen, and the drift in sensitivities was
considered in data analysis. Some data from 2013 to 2015 have been reported by
Domine et al. (2016b) and have been used by Barrere et al. (2017) to
simulate snow and ground properties, with driving data presented by Barrere
and Domine (2017). Most data were recorded by a Campbell CR1000-XT data
logger, except soil temperature and volume water content, which were recorded
by an Em50 logger from Decagon (now Meter). A Reconyx time lapse camera
taking several pictures a day was installed in summer 2016. It was replaced
and reoriented in July 2018. Similar cameras at other sites within 3 km were
present starting in 2015. In summer 2016, the SALIX meteorological station,
fairly similar to the TUNDRA station described here (except there was no
CNR4), was deployed 9 km up-valley from TUNDRA (at 73.1816∘ N,
79.7454∘ W). Data from that station were occasionally used for
filling data gaps.
There is no small-scale topography within the polygon (Fig. 1) and in
particular no hummock or tussock. The permafrost there has been described by
Fortier and Allard (2005). It is several hundred meters thick with an
active layer 20 to 35 cm deep. The visibly (dark) organic-rich surface soil
layer is 2 to 10 cm thick. Soil samples from a vertical profile taken with a
vertical resolution of 5 cm on 3 July 2017 were analyzed for organic carbon
content using the procedure detailed in Gagnon et al. (2019). We chose a
spot where the thawed layer was deepest, 30 cm, within an observed range of
15–30 cm. There was less moss at this spot than in most other places on the
polygon. The carbon content decreased from 8.3 to 0.3 kg C per kilogram of dry
soil between 0–5 and 20–25 cm depths. A graph is shown in Fig. S1,
along with soil density and water weight fraction. Soil granulometry was
analyzed as mentioned in Domine et al. (2016b) for three depth intervals:
0–5, 10–15 and 20–25 cm (Fig. S2). The 0–5 cm sample was silt loam
with a unimodal asymmetric distribution of grain size centered at 51 µm, the 10–15 cm sample was silt loam with a bimodal distribution at 17.4 and
59.0 µm, and the 20–25 cm sample was sandy loam with a bimodal
distribution at 11.6 and 152.5 µm.
Triplicate measurements of the ground thermal conductivity were performed on
9 July 2013 at 5 cm depth, within 2 m of the TP08 post, with a TP02 heated
needle probe from Hukseflux, yielding values of 0.159±0.004 W m-1 K-1. The ground temperature was 6.3 ∘C, and the
fractional water volume content measured with an EC5 probe from Decagon was
0.151. A vertical profile of the ground thermal conductivity was measured on
29 June 2014, about 5 m away from the TP08 post, showing values increasing
from 0.235 W m-1 K-1 at 3 cm depth to 1.481 W m-1 K-1 at
15 cm depth. Another vertical profile of thermal conductivity was measured
on 3 July 2017, within 2 m of the 29 June 2014 pit, when soil samples for
carbon analysis were taken. Plots of the vertical profiles, along with the
associated profiles of fractional water volume content (obtained with a
Decagon EC5 sensor) and temperature, are shown in Fig. S3. Significant
spatial variability in thermal conductivity is observed since values near
the surface vary by a factor of 3 within a few meters. Over the years,
about eight soil pits were dug in the TUNDRA polygon for physical measurements,
sampling, and instrument installation and maintenance. Variations in soil
color and texture were visible to the eye. In particular the thickness of
the darker surface soil layer, presumably organic-rich, varied between 2 and
10 cm. All soil physical variables mentioned here therefore varied within
the polygon. More extensive pit digging for extra measurements risked
modifying the soil properties in the polygon.
Vegetation in these polygons has been detailed in Gauthier et al. (1995).
It consists mostly of graminoids, sedges and mosses with some prostrate
ligneous species: Salix arctica and S. herbacea. Vegetation height does not exceed 5 to 10 cm as
there is no erect vegetation (Fig. 1). The spectral albedo of the site was
measured on 11 July 2015 around 17:25 UTC. The sky was clear, but the
atmosphere was slightly foggy. The instrument was an SVC spectrometer
equipped with an integrating sphere, one Si photodiode detector for the
visible and near-infrared range, and two InGaAs photodiode detectors for the
shortwave infrared range. Two spectra were recorded over the 346–2513 nm
wavelength range. They are shown in Fig. S4 and are essentially identical.
The broadband albedo derived from the average of these spectra is 0.18.
Time series of air temperature, relative and specific humidity,
wind speed, and atmospheric pressure at the TUNDRA site on Bylot Island.
In May 2014, 2015, 2017, 2018 and 2019, field measurements were performed in
several spots in Qarlikturvik valley. Around the TUNDRA sites, over 100
measurements of snow depth were performed to obtain a more spatially
representative value of snow depth than the one spot measured automatically.
Snow pits were dug to observe the stratigraphy and measure vertical profiles
of density and SSA. Density was measured by weighing snow sampled with a 100 cm3 box cutter. SSA was measured by infrared reflectance at 1310 nm
using the DUFISSS instrument (Gallet et al., 2009). During each campaign, a
pit was dug within 3 m of the thermal conductivity post. Two to seven pits
were dug elsewhere in the valley to assess spatial variability. Logistical
difficulties did not allow a field campaign in May 2016.
Air humidity relative to liquid water. (a) Uncorrected data from
the first RH sensor (red points), used in 2013–2014, and data from the
second RH sensor, used in 2014–2019 (blue points), after multiplication by
1.045 to ensure maximum RH reached 100 % in summer. The ice-saturating RH
values relative to supercooled water, based on the equations of Huang (2018), are also shown. (b) Corrected data for 2013–2014. Both data sets were
corrected with Eq. (1) or (2) to ensure maximum winter values were at
the ice-saturating RH. Only one of the two corrected data sets is shown for
clarity.
Driving data quality check and correction
Several environmental factors and problems with instruments can affect data
quality. For example, since the tripod is on permafrost, ground thawing and
freezing may modify its leveling, which was adjusted in early July every
year. However, further shifting can take place later in the summer. Some
years, this produced a slight offset in snow depth and in the CNR4 leveling.
In winter, frost can build up on the anemometer and block it. Frost can
obscure the CNR4 top windows, producing erroneous measurements. All these
difficulties were thoroughly investigated and corrected for. In 2016, large
surface plates were placed under each tripod leg, and this significantly
reduces tripod movement and tilting. The treatments done to each driving and
validating variable are detailed below.
Air temperature
Air temperature was measured with a ventilated HC2-S3-XT sensor at 2.3 m
height. Data from 2013–2014 were missing because of sensor failure, but we
used CR1000 data logger temperature instead. That sensor was slightly
sensitive to radiation. Based on several years of simultaneous temperature
measurements of the HCS2-S3-XT and CR1000 sensors, we corrected the CR1000
sensor values by adding a linear function of downwelling SW radiation whose
coefficients were optimized to obtain a zero bias. The root mean square deviation (RMSD) was
0.784 ∘C. The sensor was replaced in July 2014, and there was no
other data gap. TUNDRA air temperature data were compared to those from the
SILA, CAMP and SALIX stations. All differences could be readily explained by
topography and basic meteorological concepts, such as katabatic flow at the
bottom of the valley, which led to colder air at TUNDRA in winter. We are
thus confident in the reliability of the air temperature data. The
temperature time series is shown in Fig. 2.
Relative humidity
Relative humidity (RH) was also provided by the HC2-S3-XT sensor. This is a
HUMICAP thin-film capacitive sensor which provides RH relative to liquid or
supercooled water, not ice. It needs to be calibrated, and the calibration
provided by the manufacturer was checked. We found that for the second
sensor, installed in 2014, RH never reached 100 % in summer. We therefore
multiplied the value obtained by 1.045 so that the 100 % value was reached
about as frequently as the first year. Regarding winter data, we observed
that by plotting RH vs. temperature, maximum values deviated from the ice
saturation line (Fig. 3). The first sensor gave lower values for the
2013–2014 period, while the second sensor gave higher values. We corrected
the data so that maximum values for temperature < 0 ∘C
coincided with the ice saturation values. For the ice and supercooled water
saturation vapor pressures, we used the equation of Huang (2018). Huang
does not mention an accuracy for the supercooled water pressure values, but
comparison with measured values available at
https://www.engineeringtoolbox.com/water-supercooled-vapor-pressure-d_1910.html (last access: 5 February 2021) revealed an excellent agreement. At
-40 ∘C, the value of Huang was 1.3 % higher than the measured
value (19.16 vs. 18.91 Pa). At -10 ∘C, the difference is just
0.024 % (286.57 vs. 286.50 Pa). In Fig. 3, we plotted the ice-saturating
RH derived from Huang's equations. To make our data fit the ice values, we
used Eqs. (1) and (2) for the 2014–2019 and 2013–2014 data,
respectively, where T is temperature in Celsius and RH in percent.
1RHcorrected=(RHmeasured×1.045)(0.0031×T+0.77)+212RHcorrected=RHmeasured(0.00065×T+0.75)+24
Specific humidity
Many models use specific humidity rather than RH relative to water as input
variable, and we therefore also provide that variable in grams of water per kilogram of
moist air. To calculate the partial pressure of water vapor, we used
Eq. (17) of Huang (2018). We used PV = nRT for the gas equation of
state, where P is pressure, V the volume considered, n the number of moles in V, R the gas constant and T temperature. Values used in the calculations
are 18.01528 g for the molar mass of water, 8.3145 J K-1 mol-1 for R
and 28.9647 g for the molar mass of dry air. The humidity time series are
shown in Fig. 2.
Raw longwave downwelling radiation measured by the CNR4. The
near-zero values in winter over extended periods indicate the presence of
frost on the pyrgeometer window. The data for those periods are thus
invalid.
Wind speed
Data are from a cup anemometer at 2.3 m height. In winter the anemometer
frosted up during some stable weather periods and was blocked about 2–3 weeks each year. Gaps were filled with data from Young anemometers at SILA,
CAMP or SALIX, adjusted using correlations. Often, two or three of the
anemometers were blocked simultaneously. One gap could not be filled in
November 2018 because all anemometers were blocked. Since these blocking
episodes all happened under low wind speed (< 2 m s-1; usually much
less), that 2-week gap was filled with similar low-value data from another
period. Since the threshold value for cup anemometers is higher than for
Young anemometers, when the cup value was 0 for extended periods, we used
the value provided by the Young from CAMP. Each time, we checked that the
CAMP values were quite low, and in any case < 0.4 m s-1. The wind
speed time series is shown in Fig. 2. A flag indicates gap-filled values
(0: data; 1: gap-filled) in the data repository. A total of 842 out of 52 202 values were gap-filled
(1.6 %).
Atmospheric pressure
We did not measure atmospheric pressure. We relied on measurements performed
by Environment and Climate Change Canada (ECCC, https://climate.weather.gc.ca, last access: 6 January 2021), who operate a station at Pond Inlet
(72.6951∘ N, 77.9600∘ W; 60 m a.s.l.), 84.1 km to the southeast of our site, and
another station at Cape Liverpool (73.6681∘ N, 78.2942∘ W; 2 m a.s.l.), 79.5 km
to the northeast of our site. We present the average of both values, bearing in
mind that the altitude of our site is probably close to 20 m. The
atmospheric pressure time series is shown in Fig. 2.
Longwave downwelling radiation
The pyrgeometer provides a raw signal, which is the difference between
LW↓ and its own LW emission. The raw data were corrected for the
drift in sensor sensitivity, which decreased from 6.57 to 6.055 µV W-1 m2, an 8.5 % decrease, assuming a
constant drift over the 2193 d of use. LW↓ is then obtained
by removing the instrument LW emission, and this is done using the
temperature value provided by the CNR4 temperature sensor. Figure 4 shows
that the raw signal remained close to zero over extended periods. This
indicates similar temperatures for the sensor and the source (upper
atmosphere or clouds), which is only possible if the source is close to the
sensor. This is explained by the presence of frost on the pyrgeometer
window. The 5 min hourly heating by the CNF4 was therefore not sufficient
to prevent frost build-up. Figure 4 however shows some periods in winter
with no frost, e.g., January 2016. The CNR4 temperature sensor data were used
to perform the LW↓ correction and obtain the true
LW↓ data over those periods. These were used to determine the
correlation between ERA5 reanalyses (Hersbach et al., 2020) and our values.
ERA5 values were noticeably lower than CNR4 ones. For winter (i.e., here the
21 October to 14 April periods, which includes all frost-affected events)
non-frost-affected values the correlation is
LW↓CNR4=LW↓ERA5×0.7900+57.814Wm-2.
This is based on 10 850 values and R2=0.50. Invalid data were replaced
with ERA5 values modified with Eq. (3). During the first year the CNR4
temperature sensor was not functioning, preventing the temperature
correction of the raw signal. The correlation between all valid CNR4 values
(including summer) and ERA5 values was
LW↓CNR4=LW↓ERA5×0.9233+41.777Wm-2.
Equation (4), with R2=0.79, was used to obtain a time series of
LW↓ for the first year of this study. It could be argued that
for winter 2013–2014, Eq. (3) should be used. However, for the low
LW↓ values of that period, the difference between Eqs. (3) and (4) is small, between 0 and 10 W m-2. Furthermore, a sudden change in
equation may create an undesirable discontinuity. A flag indicates
gap-filled values (0: CNR4 data; 1: modified ERA5 data) in the data repository. Overall, 21 126
out of 52 202 LW↓ values were gap-filled (40 %).
Time series of downwelling longwave radiation (from our CNR4
pyrgeometer and from ERA5 reanalyses), downwelling shortwave radiation (both
from our CNR4 pyranometer and from ERA5 reanalyses), hourly precipitation and
cumulated seasonal precipitation (snow and snow-free periods).
The CNR4 and ERA5 LW↓ values are shown in the same graph in Fig. 5 for comparison. Over the whole period investigated, LW↓CNR4 values are higher than LW↓ERA5 values by an
average value of 24.5 W m-2. This large difference is discussed and the
validity of our measurements confirmed when we present LW↑ values
in Sect. 4.2.
Shortwave downwelling radiation
The calibration constant for the upward-looking pyranometer of the CNR4
drifted from 15.37 to 15.12 µV W-1 m2, a 1.65 % change. The data were adjusted for this drift. The SW↓ radiation showed
a -2 W m-2 offset, and a value of 2 W m-2 was added to all data. As
for the pyrgeometer, frost built up on the upper pyranometer window.
However, since there is no SW↓ to measure at 73∘ N
most of the winter, data are overall less affected. In late March and early
April and to a lesser extent in late October, some albedo values were
anomalously low, around 0.4. The corresponding SW↓ values were
anomalously high (Fig. S5), much higher than ERA5 values, and this
happened during periods when pyrgeometer data indicated the presence of
frost and when the sky was clear. Our CNR4 is over essentially flat
terrain, so slope effects (Picard et al., 2020) can be ruled out. We propose
that frost sublimated on the south side of the hemispherical pyranometer
window, allowing direct radiation to reach the pyranometer, while it
remained on its north side, scattering extra radiation into the pyranometer.
The pyrgeometer window is flat, so frost there is less likely to be
sublimated by radiation. Detailed data analysis showed that this process was
more intense for solar zenith angles between 71 and 81∘. While we
cannot explain this in full detail, it is consistent with the idea that
there exists an optimal geometry to maximize the scattering effect we
propose. For those SW↓ data, we replaced them with corrected
ERA5 data. For non-frosted conditions, ERA5 and CNR4 values are extremely
well correlated:
SW↓CNR4=SW↓ERA5×0.9956-3.7033Wm-2.
The RMSD is 19.51 W m-2. For frost-affected periods, we therefore used
Eq. (5) to obtain SW↓ values. However, since some ERA5
values were probably underestimated, this resulted in some albedo values
> 1, which is not consistent with a sound radiation budget. Some
data users may decide to modify some of the ERA5-derived SW↓
values presented here to ensure a reasonable albedo value, probably around
0.8. However, we left the ERA5-derived SW↓ values unchanged,
leaving decisions regarding modification to the data users. A flag specifies
when modified ERA5 values were used (0: CNR4 data; 1: modified ERA5 data) in the data repository.
Overall, 2550 out of 52 202 SW↓ values were gap-filled (4.9 %).
Finally, instrument noise yielded non-zero data even during the polar night.
We set SW↓ to 0 when the ERA5 reanalysis value was 0.
A slight tilt due to ground freezing and thawing, always less than
1.5∘ and usually less than 1∘, was observed most years on
the CNR4 during maintenance. Given the generally high solar zenith angle,
this may significantly affect SW↓ radiation measurement under
clear-sky conditions. However, given the excellent correlation with ERA5, we
conclude that these effects were probably not important. We are also
providing ERA5 data for comparison but see no objective reason not to
recommend our CNR4 data over ERA5 data. The SW radiation time series, both
from CNR4 and ERA5, are shown in Fig. 5.
Precipitation
There is a Geonor 200 precipitation gauge with a single alter shield at the
CAMP site, 1.7 km from our TUNDRA site, but most of the time, it did not
function properly. We therefore relied mostly on data from the ECCC Geonor
gauges at Pond Inlet, 84.1 km to the southeast, and Cape Liverpool, 79.5 km to the
northeast (Fig. 1). Geonor gauges measure the weight of a glycol bath into which
precipitation falls. There is noise in the data, for example due to wind,
which produces small positive and negative precipitation events that need to
be corrected and filtered. ECCC does not detail their procedure. The
negative signals from their gauges appear to have been compensated (e.g., by
reducing the positive signal from the subsequent positive event). There are
numerous small (< 0.2 mm in their hourly data) isolated positive
events, and we wondered whether those might just be noise. We examined
isolated events that had no other precipitation 10 h before and after them.
By comparing these events to our time lapse photos and by considering the
observer's remarks at Pond Inlet, we concluded that most of these isolated
events were real, even though about 30 % of them may be noise because
for example, our cameras revealed blue-sky conditions. We estimate that errors due
to such noise amount to less than 4 mm yr-1. Precipitation needs to be
corrected for undercatch under windy conditions, and we used the equations of
Kochendorfer et al. (2017) for rain and snow. The threshold for rain and snow
was set at +0.5 ∘C. Examination of the air temperature when
observations at the Pond Inlet airport indicated that precipitation was snow
(on the ECCC web site) reveals that this threshold is sensible. We therefore
used the gauge data from both ECCC sites. We determined the phase at each site
from the temperature there, also given by ECCC. Lastly, we corrected the amount of
precipitation using the local wind speed given by ECCC. To obtain
precipitation at our site, we averaged both ECCC values and determined the
phase at our site from our temperature measurement.
There were a few data gaps in the ECCC data sets. In that case we just used
data from one of the two sites. There was a gap at Cape Liverpool from 30 November 2017 to 15 April 2018 (115 d) and two gaps at Pond Inlet from 12
to 29 April 2016 (18 d) and 13 February to 11 April 2017 (58 d). There
was also a 24 d period from 8 September to 1 October 2017 when our gauge
at the CAMP site functioned well and measured a greater amount of
precipitation than the ECCC gauges. During that period, we therefore used
the CAMP data. Figure 5 shows hourly precipitation time series, separated as
rain or snow. Overall, considering errors in undercatch correction, the
distance between the instruments and our site, and the instruments' noise,
we estimate the error in the precipitation data provided to be around
20 %.
We also provide cumulated seasonal precipitation data for periods when there
was snow on the ground and periods when the ground was snow-free. Snow onset
is the first day when there is a continuous and permanent snow cover. Often,
the first snowfall melted partially or completely so that there is some
arbitrary character in determining the snow onset date. For example, on 7 September 2017 a significant snowfall resulted in complete snow cover. That
snow had mostly melted when an important snowfall that lasted the whole
season happened on the evening of 17 September so that we retain 17 September as
the snow onset date. A picture on 17 September (Fig. S6) shows what was
left of the 7 September snowfall to illustrate our choice. Meltout date is
when the winter snow cover has almost completely disappeared. Large snow
drifts melt later. A picture in Fig. S6 shows these remaining drifts on 8 June 2019, when we consider the snow to have melted out. Occasional late spring
snowfalls that occur after meltout were added to the summer precipitation.
Snow onset and meltout dates were determined from snow gauges (present at
TUNDRA and CAMP) and, when available, time lapse photographs. For 2013 and
2014, before the deployment of several time lapse cameras in the valley, we
also used satellite images to determine snow dates, as detailed in Domine
et al. (2018b). Table 2 reports the snow onset and meltout dates that we
used. Cumulated seasonal precipitation time series are shown in Fig. 5.
Note that winter 2013–2014 was an exceptionally low-snow year.
Snow onset and meltout dates at the TUNDRA site, used to determine
cumulated seasonal precipitation.
Snow yearSnow onsetMeltout2013–201411 October 20137 June 20142014–201512 September 201413 June 20152015–20161 October 201515 June 20162016–20173 October 201618 June 20172017–201817 September 201715 June 20182018–20198 October 20187 June 2019Validating data quality check and correction
Data for validation consist of monitoring data and snow pit measurements and
observations every year in May except in 2016.
Time series of SW↑ radiation from the CNR4 radiometer,
albedo, LW↑ radiation from the CNR4 radiometer and surface
temperature from the IR 120 infrared sensor. Most
albedo values over snow-covered surfaces that are very different from 0.8
are due to the use of ERA5 SW↓ values.
Shortwave upwelling radiation and albedo
CNR4 values for SW↑ radiation were corrected for sensitivity
drift, similarly to SW↓ radiation. The sensitivity changed from
15.74 to 15.39 µV W-1 m2 between 2013 and 2019, a 2.27 %
change. The downward-looking pyranometer does not frost up. Values were set
to zero when ERA5 SW↓ values were zero. The error due to the
tilt of the CNR4 discussed in the case of the downwelling radiation is
probably negligible here since radiation is diffuse. SW↑ values
are affected by the presence of the tripod and of the solar panel it
carries. Ideally, corrections can be performed, as done, for example, by Wright et
al. (2014). Values given here are uncorrected for the presence of the tripod
and solar panel, but we show in Fig. S7 the geometry of the system so that
the calculations could be performed. However, Wright et al. (2014) had a
geometry less favorable than ours, and their correction was on the order of
1 % so we do not expect this correction to be essential. Values thus
obtained are reported in Fig. 6. The high SW↑ value of 843 W m-2 on 3 June 2018 at noon local summer time (UTC-4) is most likely real since it corresponds to
a high SW↓ value of 1151 W m-2 at the same time (Fig. 5).
Other high values also coincide in the SW↑ and SW↓
data. These high values can be caused by thin clouds over snow which cause
multiple reflections and amplify radiation. Partial cloud cover can also
lead to significant radiation amplification. These effects were predicted long ago (Nack and Green, 1974) and have been evidenced by studies
focusing on UV radiation, with amplification sometimes exceeding a factor of
2 (McKenzie et al., 1998; Weihs et al., 2000; Lee et al., 2015), but the
processes are similar for visible wavelengths and very likely explain these
high values. ERA5 values do not seem to account for these processes, and this
is one reason why we recommend the use of our SW↓ data over
ERA5.
Time series of snow depth monitored by an automatic snow gauge.
The averages of over 100 spot measurements in mid-May around the TUNDRA site
are also shown for 5 of the 6 years.
Albedo obtained from the SW↑/ SW↓ ratio is shown in the
same panel as SW↑ in Fig. 6. Values very different from 0.8
during snow-covered periods are due to the use of modified SW↓
ERA5 values, while SW↑ are all from the CNR4. A flag in the data
file indicates which albedo data used ERA5-derived SW↓ values
(0: CNR4; 1: modified ERA5).
Longwave upwelling radiation and surface temperature
Both the downward-looking pyrgeometer and the IR120 surface temperature
sensor provide information on LW↑ radiation. The CNR4 LW↑ sensor sensitivity changed from 6.38 to 5.95 µV W-1 m2
between 2013 and 2019, a 7.23 % change, and the signals were corrected
accordingly. As for the LW↓, no data are available the first
year. The IR120 provides surface temperature Ts but was not
recalibrated. Ts and LW↑ are linked by the Stefan–Boltzmann
equation so that IR120 data may be used to fill the first year of missing
data from the CNR4 LW↑. Both data sets are quite similar. The RMSD
between the CNR4 LW↑ and the IR120 LW↑ calculated from
Ts using a surface emissivity ε=1 for the 2014–2019
period is 7.51 W m-2. The RMSD can be reduced if different ε values are used during snow-covered and snow-free periods. For the
2014–2015 winter, RMSD = 6.12 W m-2 is obtained for ε=1.027. For the 2015 summer, RMSD = 10.54 W m-2 is obtained for
ε=0.991. Obviously ε cannot be > 1,
and the 1.027 value only indicates a systematic shift between both sensors.
This is not surprising as the wavelength range sensed by the IR120 is
narrower (Table 1), and small errors in the calibrations are inevitable. By
comparing CNR4 drift-corrected data and IR120 uncorrected data over 5 years,
we did not however note any detectable drift in the IR120 sensitivity. For
the first year, the CNR4 LW↑ data gap was filled with the IR120
data, with the optimal emissivities found above for the snow-covered and
snow-free periods. A flag in the data file indicates IR120-filled data. CNR4
LW↑ time series and surface temperature time series from IR120 are
plotted in Fig. 6. IR120 data have not been modified.
We compared our CNR4 LW↑ data with ERA5. ERA5 LW↑ was on
average 16.7 W m-2 lower, showing that ERA5 underestimates the
temperature of the surface. The simplest explanation is that ERA5
LW↓ is underestimated, and this reduces surface warming. This
confirms that our LW↓ data, which are on average 24.5 W m-2
higher than ERA5, are probably correct, and we recommend their use over ERA5
for our site.
Mean values and standard deviation of snow depths measured around
the TUNDRA site in May.
DateMean depthStandard deviation14 May 201416.2 cm13.7 cm12 May 201525.3 cm13.1 cmMay 2016No dataNo data13 May 201741.0 cm10.9 cm14 May 201844.5 cm13.4 cm17 May 201929.5 cm13.8 cm
Time series of snow temperature from Pt100 thermistors and ground
temperature and liquid water volume fraction from 5TM probes. In 2013–2014,
snow temperature data were limited to 7 cm height (low snow height that year)
with a reading from the TP08 probe every other day at 05:00 local summer time (UTC-4). The heights of
the snow temperature sensors were lowered by 2 cm in July 2018. The vertical
black line in July 2018 shows the date of this height change (H. ch. in
legend box). There are negative spikes due to instrumental noise on the soil
temperature at 15 cm depth data until 2016 and on the 21 cm depth data in
summer 2016.
Snow depth
Continuous snow depth data from the TUNDRA snow gauge are shown in Fig. 7.
To facilitate reading, snow-free periods were assigned a snow depth
value of zero. However, snow depth is highly spatially variable because of the
relief at the 10 to 20 m scale in the ice-wedge polygon terrain. Therefore,
additional manual snow depth measurements were taken in May 2014, 2015,
2017, 2018 and 2019 at several hundred random spots around the tundra site.
The means and standard deviations are shown in Table 3. Snow depth
measurements were also done in numerous spots in the whole valley. This
confirmed that spring 2018 was indeed the snowiest year we experienced, and
spring 2014 had by far the lowest snow depth everywhere in the valley. The snow
depth data of Table 3 are therefore representative of the climatology at
least at the 20 km scale.
Snow temperature
Snow temperatures were measured with Pt100 thermistors installed in July 2014. For the 2013–2014 season, we provide temperature given by the TP08
heated needle probe, which produced one value every other day at 05:00 local
summer time, when a thermal conductivity measurement was performed. In
2014–2015, there were only two thermistors, at heights of 2 and 17 cm. In July 2015, thermistors were added at 7, 27 and 37 cm. In July 2018, all five thermistors were lowered by 2 cm to 0, 5, 15, 25 and 35 cm. All data
> 0 ∘C were deleted. Data when no snow was present on
the ground have also been deleted, based on snow height data or time lapse
images. However, the snow gauge is about 6 m away from the thermistor post,
and only the top of the post is in the field of view of the camera. Another
criterion for the presence of snow is the temperature gradient in the set of
sensors. When snow is present, the lowest sensor is expected to be warmer,
at least until spring warm up, when the temperature gradient reverses.
However, all these criteria are not 100 % certain, and we may have
included some data in the absence of snow. Data from upper sensors not
covered by snow have not been deleted. Snow temperature data are shown in Fig. 8.
Time series of snow thermal conductivities at heights of (a) 22 cm, (b) 12 cm and (c) 7 cm for the first winter and 2 cm for subsequent
winters and (d) soil thermal conductivity at 10 cm depth. In 2013–2014
there was insufficient snow to cover the top two TP08 probes, which were at
heights of 17 and 27 cm.
Vertical profiles of snow density and SSA measured within the
TUNDRA polygon in mid-May of 5 years.
Ground temperature and liquid water volume content
These variables were measured using 5TM sensors from Decagon placed within 1 m of the TP08 post. The deepest sensor was placed at the base of the summer
thawed layer. Decagon documentation specifies that “the 5TM determines
volumetric water content (VWC) by measuring the dielectric constant of the
media using capacitance/frequency domain technology. The sensor uses a 70 MHz frequency, which minimizes textural and salinity effects, making the 5TM
accurate in most soils. The 5TM measures temperature using an onboard
thermistor.” Regarding temperature, offsets of up to 0.5 ∘C,
constant over time, were noticed during soil freezing. All temperatures were
corrected so that T=0∘C during the zero-curtain periods.
Regarding VWC, the calibration provided by Decagon was used. For mineral
soils, a 2 % accuracy is claimed by the manufacturer. For other soils,
3 % is claimed. This lower accuracy probably applies to the top two
sensors at 2 and 5 cm depth, where the soil has a significant organics
content. Due to battery failure, there is a data gap between 20 March and 12 July 2016. Before March 2016, the temperature sensor at 15 cm depth showed
very frequent spikes in summer that gave readings lowered by 1.5 to
3 ∘C, which is why the plots appear noisy. The same applies to the
21 cm sensor in August 2016. The causes are unknown. Data are shown in Fig. 8.
Snow and soil thermal conductivity
Measurement methods using the TP08 heated needle probe are detailed in
Domine et al. (2015). Data from the first three winters have already been
reported in Domine et al. (2016b) and Domine et al. (2018a). Figure 9
shows measurements for all 6 years at three heights. In 2013–2014, only the 7 cm
needle was covered. In July 2014 the sensors were lowered to 2, 12 and 22 cm.
Soil thermal conductivity values only show significant variations between
the thawed and frozen state, as frequently observed in soils (Smerdon and
Mendoza, 2010). Thawed and frozen values are around 0.75 and 1.8 W m-1 K-1, respectively, with little variations between years. Values may vary
with water or ice content, but this was not investigated here. In the frozen
state, many heating curves were of insufficient quality because of the
limited heating, and those data were discarded (see Domine et al., 2015, for
details), hence the missing data points.
Snow thermal conductivity is a valuable proxy for snow type. Soft depth hoar
always has a low value, and, for example, the very low thermal conductivity
value at 2 cm height in 2014–2015 (mostly < 0.035 W m-1 K-1) is indicative of the presence of very soft depth hoar, as observed
in May 2015 during the field campaign. In contrast, the high values in
2015–2016 (0.2 to 0.35 W m-1 K-1) indicate that depth hoar was
probably indurated, due either to rain on snow (ROS) that formed a hard
refrozen layer or to high winds during precipitation that formed a hard wind
slab. On 1 October 2015 ROS took place just after snow onset, and on
14–15 October a 36 h storm with wind speeds exceeding 10 m s-1 and
precipitation in excess of 10 mm took place so that either of both options
is possible. We could not get to Bylot Island in spring 2016 for snow
observations. However, snow pit observations near Pond Inlet on 15 May 2016
indeed revealed the presence of a 10 to 15 cm basal layer of indurated depth
hoar.
It has been reported that the heated needle probe method produced a negative
artifact in the measurement of snow thermal conductivity (Riche and
Schneebeli, 2013). A correction
algorithm has just been proposed by Fourteau et al. (2021). Briefly, the amount
of correction decreases with increasing snow density, and a multiplicative
factor of about 1.1 for dense wind slabs and 1.5 for soft depth hoar must be
applied. Data presented are uncorrected. Note that here the depth hoar
thermal conductivity at 2 cm in 2014–2015 has values around 0.02 W m-1 K-1, lower than air, and after correction these values will be
around 0.03, more plausible for very light and incohesive depth hoar.
Field observations of snow
Snow density and SSA profiles cannot be monitored automatically today.
Instead, vertical profiles of these variables were measured at the TUNDRA
site in mid-May 2014, 2015, 2017, 2018 and 2019 during field expeditions.
The snow pits were dug in the actual polygon of the station, within 3 m of
the snow thermal conductivity post. Data are shown in Fig. 10. We stress
here that these profiles are highly variable in space because of wind
erosion and redeposition, which results in heterogeneous and often
discontinuous snow layers. Attempting to reproduce the details of these
profiles using 1-D model simulations is therefore not very meaningful. To
illustrate the spatial variability in these variables we report in Figs. S8 and S9 additional profiles measured in the valley, in the absence of
erect vegetation, i.e., in places where there is no Salix richardsonii as these erect shrubs
significantly affect snow properties (Domine et al., 2016a). The coordinates
and dates of these additional profiles are reported in Table S1.
Data availability
The driving and validating data, including snow pit data, are available on
the Nordicana D repository, 10.5885/45693CE-02685A5200DD4C38
(Domine et al., 2021).
Conclusion
A 6-year time series of driving data for a high Arctic herb tundra site is
presented. A unique set of validation data is provided, which includes time
series of snow and soil thermal conductivity. Vertical profiles of snow
density and specific surface area in mid-May are also provided for all years
except 2016. One important objective of these data is to assist in the
improvement and validation of snow physics models, which today have great
difficulties in simulating high Arctic snowpack properties. We plan to
update the data sets on the Nordicana D repository by adding extra years of
data whenever possible. The COVID-19 pandemic prevented us from accessing
the site in spring and summer 2020 and spring 2021, but we will do our best
to maintain our effort in subsequent years.
The supplement related to this article is available online at: https://doi.org/10.5194/essd-13-4331-2021-supplement.
Author contributions
FD designed the research and obtained funding. DS and FD deployed and maintained
the instruments. FD and GL analyzed the data and prepared the data files. MP, FD and
MBB performed the fieldwork. FD wrote the paper with input from GL and
comments from MBB, DS and MP.
Competing interests
The authors declare that they have no conflict of interest.
Disclaimer
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Acknowledgements
Logistical support was provided by the Polar Continental Shelf Program and by Sirmilik National
Park. We are grateful to Gilles Gauthier and Marie-Christine Cadieux for
their decades-long efforts to build and maintain the research base of the
Centre d'Etudes Nordiques at Bylot Island. Assistance with fieldwork by
Mathieu Barrère, Mikael Gagnon and Marianne Valcourt is gratefully
acknowledged. Christophe Kinnard kindly pointed out an error in the Campbell
program for the CNR4. This allowed us to make the best use of available
radiation data. We acknowledge helpful discussions with Ghislain Picard and
Marie Dumont regarding LW and SW radiation and comments on the manuscript
by Warwick Vincent. We thank Richard Essery and Cécile Ménard (University of
Edinburgh) for expressing interest in this work and encouraging its
completion. Compute Canada assisted with data storage and handling. This
work is a contribution to the IASC project Terrestrial Multidisciplinary
distributed Observatories for the Study of Arctic Connections (T-MOSAiC).
Financial support
This research has been supported by the Natural Sciences and Engineering Research Council of Canada (Discovery Grant), the French Polar Institute (IPEV) (grant no. 1042), the Fondation BNP Paribas (grant no. APT), and the Horizon 2020 program of the European Commission (INTAROS; grant no. 727890).
Review statement
This paper was edited by Kirsten Elger and reviewed by Joshua King and one anonymous referee.
ReferencesBarrere, M. and Domine, F.: Snow, soil and meteorological data at Bylot
Island for simulating the permafrost thermal regime and evaluating output of
the SURFEXv8 land surface scheme, v. 1.0 (1979–2015), Nordicana, D29, 10.5885/45460CE-9B80A99D55F94D95, 2017.Barrere, M., Domine, F., Decharme, B., Morin, S., Vionnet, V., and Lafaysse, M.: Evaluating the performance of coupled snow–soil models in SURFEXv8 to simulate the permafrost thermal regime at a high Arctic site, Geosci. Model Dev., 10, 3461–3479, 10.5194/gmd-10-3461-2017, 2017.
Bartelt, P. and Lehning, M.: A physical SNOWPACK model for the Swiss
avalanche warning Part I: numerical model, Cold Regions Sci. Tech., 35,
123–145, 2002.Bilodeau, F., Gauthier, G., and Berteaux, D.: The effect of snow cover on
lemming population cycles in the Canadian High Arctic, Oecologia, 172,
1007–1016, 10.1007/s00442-012-2549-8, 2013.Boike, J., Kattenstroth, B., Abramova, K., Bornemann, N., Chetverova, A., Fedorova, I., Fröb, K., Grigoriev, M., Grüber, M., Kutzbach, L., Langer, M., Minke, M., Muster, S., Piel, K., Pfeiffer, E.-M., Stoof, G., Westermann, S., Wischnewski, K., Wille, C., and Hubberten, H.-W.: Baseline characteristics of climate, permafrost and land cover from a new permafrost observatory in the Lena River Delta, Siberia (1998–2011), Biogeosciences, 10, 2105–2128, 10.5194/bg-10-2105-2013, 2013.Boike, J., Juszak, I., Lange, S., Chadburn, S., Burke, E., Overduin, P. P., Roth, K., Ippisch, O., Bornemann, N., Stern, L., Gouttevin, I., Hauber, E., and Westermann, S.: A 20-year record (1998–2017) of permafrost, active layer and meteorological conditions at a high Arctic permafrost research site (Bayelva, Spitsbergen), Earth Syst. Sci. Data, 10, 355–390, 10.5194/essd-10-355-2018, 2018.Boike, J., Nitzbon, J., Anders, K., Grigoriev, M., Bolshiyanov, D., Langer, M., Lange, S., Bornemann, N., Morgenstern, A., Schreiber, P., Wille, C., Chadburn, S., Gouttevin, I., Burke, E., and Kutzbach, L.: A 16-year record (2002–2017) of permafrost, active-layer, and meteorological conditions at the Samoylov Island Arctic permafrost research site, Lena River delta, northern Siberia: an opportunity to validate remote-sensing data and land surface, snow, and permafrost models, Earth Syst. Sci. Data, 11, 261–299, 10.5194/essd-11-261-2019, 2019.Buckeridge, K. M. and Grogan, P.: Deepened snow alters soil microbial
nutrient limitations in arctic birch hummock tundra, Appl. Soil Ecol.,
39, 210–222, 10.1016/j.apsoil.2007.12.010, 2008.Carmagnola, C. M., Domine, F., Dumont, M., Wright, P., Strellis, B., Bergin, M., Dibb, J., Picard, G., Libois, Q., Arnaud, L., and Morin, S.: Snow spectral albedo at Summit, Greenland: measurements and numerical simulations based on physical and chemical properties of the snowpack, The Cryosphere, 7, 1139–1160, 10.5194/tc-7-1139-2013, 2013.
CEN: Climate station data from Bylot Island in Nunavut, Canada, v. 1.8
(1992–2019), on: Nordicana D2, 2018, Centre d'Etudes Nordiques, Quebec city,
2020.Connolly, R., Connolly, M., Soon, W., Legates, D. R., Cionco, R. G., and
Herrera, V. M. V.: Northern Hemisphere Snow-Cover Trends (1967–2018): A
Comparison between Climate Models and Observations, Geosciences, 9, 135, 10.3390/geosciences9030135, 2019.Derksen, C., Silis, A., Sturm, M., Holmgren, J., Liston, G. E., Huntington,
H., and Solie, D.: Northwest Territories and Nunavut Snow Characteristics
from a Subarctic Traverse: Implications for Passive Microwave Remote
Sensing, J. Hydrometeorol., 10, 448–463, 10.1175/2008jhm1074.1, 2009.Domine, F., Cabanes, A., and Legagneux, L.: Structure, microphysics, and
surface area of the Arctic snowpack near Alert during the ALERT 2000
campaign, Atmos. Environ., 36, 2753–2765, 10.1016/S1352-2310(02)00108-5,
2002.Domine, F., Gallet, J.-C., Bock, J., and Morin, S.: Structure, specific
surface area and thermal conductivity of the snowpack around Barrow, Alaska,
J. Geophys. Res., 117, D00R14, 10.1029/2011jd016647, 2012.Domine, F., Barrere, M., Sarrazin, D., Morin, S., and Arnaud, L.: Automatic monitoring of the effective thermal conductivity of snow in a low-Arctic shrub tundra, The Cryosphere, 9, 1265–1276, 10.5194/tc-9-1265-2015, 2015.Domine, F., Barrere, M., and Morin, S.: The growth of shrubs on high Arctic tundra at Bylot Island: impact on snow physical properties and permafrost thermal regime, Biogeosciences, 13, 6471–6486, 10.5194/bg-13-6471-2016, 2016a.Domine, F., Barrere, M., and Sarrazin, D.: Seasonal evolution of the effective thermal conductivity of the snow and the soil in high Arctic herb tundra at Bylot Island, Canada, The Cryosphere, 10, 2573–2588, 10.5194/tc-10-2573-2016, 2016b.Domine, F., Belke-Brea, M., Sarrazin, D., Arnaud, L., Barrere, M., and
Poirier, M.: Soil moisture, wind speed and depth hoar formation in the
Arctic snowpack, J. Glaciol., 64, 990–1002, 10.1017/jog.2018.89, 2018a.Domine, F., Gauthier, G., Vionnet, V., Fauteux, D., Dumont, M., and Barrere,
M.: Snow physical properties may be a significant determinant of lemming
population dynamics in the high Arctic, Arct. Sci., 4, 813–826, 10.1139/as-2018-0008, 2018b.Domine, F., Picard, G., Morin, S., Barrere, M., Madore, J.-B., and Langlois,
A.: Major Issues in Simulating Some Arctic Snowpack Properties Using Current
Detailed Snow Physics Models: Consequences for the Thermal Regime and Water
Budget of Permafrost, J. Adv. Model. Earth Sy., 11,
34–44, 10.1029/2018MS001445, 2019.Domine, F., Lackner, G., Sarrazin, D., Poirier, M., and Belke-Brea, M.:
Meteorological, snow and soil data from Bylot Island, Canadian high-Arctic,
for driving and testing snow and land surface models, Nordicana, D86,
10.5885/45693CE-02685A5200DD4C38, 2021.Fauteux, D., Gauthier, G., and Berteaux, D.: Seasonal demography of a cyclic
lemming population in the Canadian Arctic, J. Anim. Ecol., 84,
1412–1422, 10.1111/1365-2656.12385, 2015.Fortier, D. and Allard, M.: Frost-cracking conditions, Bylot Island,
Eastern Canadian Arctic Archipelago, Permafrost Periglac., 16, 145–161, 10.1002/ppp.504, 2005.
Fourteau, K., Hagenmuller, P., Roulle, J., and Domine, F.: On the use of heated needle probes for measuring snow thermal conductivity, J. Glaciol., in press, 2021.Gagnon, M., Domine, F., and Boudreau, S.: The carbon sink due to shrub
growth on Arctic tundra: a case study in a carbon-poor soil in eastern
Canada, Environ. Re. Commun., 1, 091001, 10.1088/2515-7620/ab3cdd, 2019.Gallet, J.-C., Domine, F., Zender, C. S., and Picard, G.: Measurement of the specific surface area of snow using infrared reflectance in an integrating sphere at 1310 and 1550 nm, The Cryosphere, 3, 167–182, 10.5194/tc-3-167-2009, 2009.Gauthier, G., Hughes, R. J., Reed, A., Beaulieu, J., and Rochefort, L.:
Effect of grazing by greater snow geese on the production of graminoids at
an arctic site (Bylot Island, NWT, Canada), J. Ecol., 83, 653–664, 10.2307/2261633, 1995.Gouttevin, I., Langer, M., Löwe, H., Boike, J., Proksch, M., and Schneebeli, M.: Observation and modelling of snow at a polygonal tundra permafrost site: spatial variability and thermal implications, The Cryosphere, 12, 3693–3717, 10.5194/tc-12-3693-2018, 2018.Hersbach, H., Bell, B., Berrisford, P., Hirahara, S., Horányi, A.,
Muñoz-Sabater, J., Nicolas, J., Peubey, C., Radu, R., Schepers, D.,
Simmons, A., Soci, C., Abdalla, S., Abellan, X., Balsamo, G., Bechtold, P.,
Biavati, G., Bidlot, J., Bonavita, M., De Chiara, G., Dahlgren, P., Dee, D.,
Diamantakis, M., Dragani, R., Flemming, J., Forbes, R., Fuentes, M., Geer,
A., Haimberger, L., Healy, S., Hogan, R. J., Hólm, E., Janisková,
M., Keeley, S., Laloyaux, P., Lopez, P., Lupu, C., Radnoti, G., de Rosnay,
P., Rozum, I., Vamborg, F., Villaume, S., and Thépaut, J.-N.: The ERA5
global reanalysis, Q. J. Roy. Meteor. Soc., 146, 1999–2049,
10.1002/qj.3803, 2020.Huang, J.: A Simple Accurate Formula for Calculating Saturation Vapor
Pressure of Water and Ice, J. Appl. Meteorol. Climatol.,
57, 1265–1272, 10.1175/jamc-d-17-0334.1, 2018.Kochendorfer, J., Nitu, R., Wolff, M., Mekis, E., Rasmussen, R., Baker, B., Earle, M. E., Reverdin, A., Wong, K., Smith, C. D., Yang, D., Roulet, Y.-A., Buisan, S., Laine, T., Lee, G., Aceituno, J. L. C., Alastrué, J., Isaksen, K., Meyers, T., Brækkan, R., Landolt, S., Jachcik, A., and Poikonen, A.: Analysis of single-Alter-shielded and unshielded measurements of mixed and solid precipitation from WMO-SPICE, Hydrol. Earth Syst. Sci., 21, 3525–3542, 10.5194/hess-21-3525-2017, 2017.Krinner, G., Derksen, C., Essery, R., Flanner, M., Hagemann, S., Clark, M., Hall, A., Rott, H., Brutel-Vuilmet, C., Kim, H., Ménard, C. B., Mudryk, L., Thackeray, C., Wang, L., Arduini, G., Balsamo, G., Bartlett, P., Boike, J., Boone, A., Chéruy, F., Colin, J., Cuntz, M., Dai, Y., Decharme, B., Derry, J., Ducharne, A., Dutra, E., Fang, X., Fierz, C., Ghattas, J., Gusev, Y., Haverd, V., Kontu, A., Lafaysse, M., Law, R., Lawrence, D., Li, W., Marke, T., Marks, D., Ménégoz, M., Nasonova, O., Nitta, T., Niwano, M., Pomeroy, J., Raleigh, M. S., Schaedler, G., Semenov, V., Smirnova, T. G., Stacke, T., Strasser, U., Svenson, S., Turkov, D., Wang, T., Wever, N., Yuan, H., Zhou, W., and Zhu, D.: ESM-SnowMIP: assessing snow models and quantifying snow-related climate feedbacks, Geosci. Model Dev., 11, 5027–5049, 10.5194/gmd-11-5027-2018, 2018.Langlois, A., Johnson, C. A., Montpetit, B., Royer, A., Blukacz-Richards, E.
A., Neave, E., Dolant, C., Roy, A., Arhonditsis, G., Kim, D. K., Kaluskar,
S., and Brucker, L.: Detection of rain-on-snow (ROS) events and ice layer
formation using passive microwave radiometry: A context for Peary caribou
habitat in the Canadian Arctic, Remote Sens. Environ., 189, 84–95, 10.1016/j.rse.2016.11.006, 2017.Lee, Y. G., Koo, J. H., and Kim, J.: Influence of cloud fraction and snow
cover to the variation of surface UV radiation at King Sejong station,
Antarctica, Atmos. Res., 164, 99–109, 10.1016/j.atmosres.2015.04.020, 2015.McKenzie, R. L., Paulin, K. J., and Madronich, S.: Effects of show cover on
UV irradiance and surface albedo: A case study, J. Geophys. Res., 103,
28785–28792, 10.1029/98jd02704, 1998.Myers-Smith, I. H. and Hik, D. S.: Shrub canopies influence soil
temperatures but not nutrient dynamics: An experimental test of tundra
snow-shrub interactions, Ecol. Evol., 3, 3683–3700, 10.1002/ece3.710, 2013.Nack, M. L. and Green, A. E. S.: Influence of clouds, haze, and smog on
middle ultraviolet reaching ground, Appl. Opt., 13, 2405–2415, 10.1364/ao.13.002405, 1974.Picard, G., Dumont, M., Lamare, M., Tuzet, F., Larue, F., Pirazzini, R., and Arnaud, L.: Spectral albedo measurements over snow-covered slopes: theory and slope effect corrections, The Cryosphere, 14, 1497–1517, 10.5194/tc-14-1497-2020, 2020.Poirier, M., Gauthier, G., and Domine, F.: What guides lemmings movements
through the snowpack?, J. Mammal., 10, 1416–1426, 10.1093/jmammal/gyz129, 2019.Riche, F. and Schneebeli, M.: Thermal conductivity of snow measured by three independent methods and anisotropy considerations, The Cryosphere, 7, 217–227, 10.5194/tc-7-217-2013, 2013.Saccone, P., Morin, S., Baptist, F., Bonneville, J.-M., Colace, M.-P.,
Domine, F., Faure, M., Geremia, R., Lochet, J., Poly, F., Lavorel, S., and
Clément, J.-C.: The effects of snowpack properties and plant strategies
on litter decomposition during winter in subalpine meadows, Plant Soil,
363, 215–229, 10.1007/s11104-012-1307-3, 2013.Smerdon, B. D. and Mendoza, C. A.: Hysteretic freezing characteristics of
riparian peatlands in the Western Boreal Forest of Canada, Hydrol.
Process., 24, 1027–1038, 10.1002/hyp.7544, 2010.Sturm, M. and Johnson, J. B.: Natural-convection in the sub-arctic snow
cover, J. Geophys. Res.-Sol. Ea., 96,
11657–11671, 10.1029/91JB00895, 1991.
Sturm, M., Holmgren, J., and Liston, G. E.: A seasonal snow cover
classification-system for local to global applications, J. Clim., 8,
1261–1283, 1995.Sturm, M. and Benson, C. S.: Vapor transport, grain growth and depth-hoar
development in the subarctic snow, J. Glaciol., 43, 42–59, 10.3189/S0022143000002793, 1997.Trabant, D. and Benson, C. S.: Field experiments on the development of
depth hoar, Geol. Soc. Am. Mem., 135, 309–322,
10.1130/MEM135-p309, 1972.Vionnet, V., Brun, E., Morin, S., Boone, A., Faroux, S., Le Moigne, P., Martin, E., and Willemet, J.-M.: The detailed snowpack scheme Crocus and its implementation in SURFEX v7.2, Geosci. Model Dev., 5, 773–791, 10.5194/gmd-5-773-2012, 2012.
Weihs, P., Webb, A. R., Hutchinson, S. J., and Middleton, G. W.:
Measurements of the diffuse UV sky radiance during broken cloud conditions,
J. Geophys. Res., 105, 4937–4944, 10.1029/1999jd900260, 2000.Wright, P., Bergin, M., Dibb, J., Lefer, B., Domine, F., Carman, T.,
Carmagnola, C., Dumont, M., Courville, Z., Schaaf, C., and Wang, Z.:
Comparing MODIS daily snow albedo to spectral albedo field measurements in
Central Greenland, Remote Sens. Environ., 140, 118–129,
10.1016/j.rse.2013.08.044, 2014.Zhang, T. J.: Influence of the seasonal snow cover on the ground thermal
regime: An overview, Rev. Geophys., 43, RG4002, 10.1029/2004rg000157, 2005.