Forcing the SURFEX/Crocus snow model with combined hourly meteorological forecasts and gridded observations in southern Norway
- The Norwegian Meteorological Institute, P.O. Box 43 Blindern, 0313 Oslo, Norway
Correspondence: Hanneke Luijting (email@example.com) and Dagrun Vikhamar-Schuler (firstname.lastname@example.org)
In Norway, 30 % of the annual precipitation falls as snow. Knowledge of the snow reservoir is therefore important for energy production and water resource management. The land surface model SURFEX with the detailed snowpack scheme Crocus (SURFEX/Crocus) has been run with a grid spacing of 1 km over an area in southern Norway for 2 years (1 September 2014–31 August 2016). Experiments were carried out using two different forcing data sets: (1) hourly forecasts from the operational weather forecast model AROME MetCoOp (2.5 km grid spacing) including post-processed temperature (500 m grid spacing) and wind, and (2) gridded hourly observations of temperature and precipitation (1 km grid spacing) combined with meteorological forecasts from AROME MetCoOp for the remaining weather variables required by SURFEX/Crocus. We present an evaluation of the modelled snow depth and snow cover in comparison to 30 point observations of snow depth and MODIS satellite images of the snow-covered area. The evaluation focuses on snow accumulation and snowmelt. Both experiments are capable of simulating the snowpack over the two winter seasons, but there is an overestimation of snow depth when using meteorological forecasts from AROME MetCoOp (bias of 20 cm and RMSE of 56 cm), although the snow-covered area in the melt season is better represented by this experiment. The errors, when using AROME MetCoOp as forcing, accumulate over the snow season. When using gridded observations, the simulation of snow depth is significantly improved (the bias for this experiment is 7 cm and RMSE 28 cm), but the spatial snow cover distribution is not well captured during the melting season. Underestimation of snow depth at high elevations (due to the low elevation bias in the gridded observation data set) likely causes the snow cover to decrease too soon during the melt season, leading to unrealistically little snow by the end of the season. Our results show that forcing data consisting of post-processed NWP data (observations assimilated into the raw NWP weather predictions) are most promising for snow simulations, when larger regions are evaluated. Post-processed NWP data provide a more representative spatial representation for both high mountains and lowlands, compared to interpolated observations. There is, however, an underestimation of snow ablation in both experiments. This is generally due to the absence of wind-induced erosion of snow in the SURFEX/Crocus model, underestimated snowmelt and biases in the forcing data.
Snow is a key element in the hydrological cycle. Seasonal snow covers large areas of the Northern Hemisphere and the Arctic. In these areas the snow cover extent in spring has reduced more rapidly over the past 40 years than over the past 90 years (Brown and Robinson, 2011; Brown et al., 2017). The largest declines in snow cover extent and duration are observed in Arctic coastal areas, e.g. in Scandinavia (Brown et al., 2017; Rasmus et al., 2015). In Norway there is a general trend towards a later start and an earlier end of the snow season, although there are large annual variabilities (Hanssen-Bauer et al., 2015, 2017). Trends in snow depth may vary with elevation, as observed for some Norwegian regions (Dyrrdal et al., 2013; Skaugen et al., 2012). Information about seasonal changes in snow duration and amounts are important for many societal applications and for Arctic ecosystem changes. An overview of changes in snow and impacts due to these changes is provided by Bokhorst et al. (2016).
In Norway 30 % of the annual precipitation falls as snow (Saloranta, 2012). Observations show that changes in the winter climate over the past 50 years, particularly since 2000, produce more winter warming and rainfall events (Kivinen et al., 2017; Vikhamar-Schuler et al., 2016). This, in turn, affects the internal snow structure, producing, for example, more wet snow conditions and ground-ice layering (Johansson et al., 2011; Vikhamar-Schuler et al., 2013). Updated information on the daily local snow properties (e.g. depth, SWE, density profile, crystal structure, etc.) and snowmelt in mountainous and lowland areas is very useful for many applications, notably local flood prediction, hydropower production planning, snow avalanche prediction, tourism and traffic flow management. Typical information needed for these applications include daily forecasts of snow properties and snowmelt (for the next days) but depending on the application, knowledge of snow conditions for the past winter(s), month(s), week(s), 3 days and the previous day is also required.
A wide range of empirically and physically based snow models have been developed and reported in the literature; see e.g. Magnusson et al. (2015). Models differ in several ways, e.g. the parameterization and simplification of snow processes, the spatial and temporal resolution or the need for input data on various weather elements. Input data therefore usually needs to be in a larger quantity and more detailed for physically based models than for the empirically based models. Empirical models often need calibration. Several snow model intercomparison projects have also been carried out, e.g. by Etchevers et al. (2004) and Essery et al. (2013). These studies show that no single model always performs best, and there is no clear link between model complexity and performance. However, physically based models, which includes prognostics of snow density and albedo, tend to perform better (Essery et al., 2013).
In Norway, both the national operational flood forecasting and hydropower companies use the HBV model for hydrological forecasting, which includes an empirical degree-day model for snow simulations (Bergstrøm, 1976; Ruan and Langsholt, 2017; Sælthun, 1996). Snow maps of depth, water equivalent, snowmelt, snow wetness and skiing conditions are also produced operationally on a daily basis and published at http://www.seNorge.no (last access: 18 June 2018) and http://www.xgeo.no (last access: 18 June 2018) (Saloranta, 2016), and these maps are used by the national snow avalanche service (Barfod et al., 2013; Engeset, 2013). Both these applications use gridded near-real-time observations of temperature and precipitation (Lussana et al., 2018a; Mohr, 2008).
Another type of forcing data for snow models are numerical weather forecasts (NWPs) from atmospheric models. NWP data provide all the basic environmental variables required by physically based snowpack models at hourly time steps (e.g. air temperature, relative humidity, wind speed, precipitation rate, incoming short- and long-wave radiation). Using sub-daily (e.g. hourly) data in snow modelling should contribute to improved representation of snowmelt processes (e.g. diurnal freeze–thaw cycles) and precipitation phase. SNOWPACK (Bartelt and Lehning, 2002; Lehning et al., 2002), SURFEX/ISBA/Crocus (Vionnet et al., 2012) and JULES (Best et al., 2011) are examples of models with multi-layer snow schemes of different complexity aiming to simulate the surface energy balance and the internal layering of the snowpack. SNOWPACK and Crocus are used in the operational snow avalanche services in Switzerland and France, respectively (Fierz et al., 2013; Lafaysse et al., 2017). Many studies show that high-resolution NWP data are very valuable for driving these snow models (see e.g. Bellaire et al. (2011, 2013); Horton et al. (2015); Quéno et al. (2016); Vionnet et al. (2016)). NWP data have also been used to drive data in hydrometeorological models (Carrera et al., 2010).
For a point location in the Columbia Mountains, western Canada, Bellaire et al. (2011, 2013) used 15 km resolution weather forecasts from the NWP Global Environmental Multiscale Model (GEM) to force the SNOWPACK model. This study was later extended to a gridded area in the same region by Horton et al. (2015), who forced the SNOWPACK model using 2.5 km resolution NWP data from the limited area model version of GEM (GEM-LAM) model. The use of NWP data as precipitation forcing for snow models was analysed and discussed by Schirmer and Jamieson (2015). They compared two NWP data sets (GEM: 15 km and GEM-LAM: 2.5 km spatial resolution) over complex mountainous terrain during wintertime and found that the highest-resolution data set worked best in terms of precipitation forecasts. Bernier et al. (2011) used a downscaling technique to account for local terrain effects on the surface temperatures that are not resolved by the low-resolution NWP model. With higher spatial resolution NWP models, the orographic precipitation is better reproduced. In the French Alps, high-resolution forecasts (2.5 km) from the AROME NWP model were used to drive Crocus (Vionnet et al., 2016). A similar study using the same AROME NWP model was carried out for the French and the Spanish Pyrenees by Quéno et al. (2016). Both these studies showed that high-resolution NWP data are very useful and promising for producing snow maps with snow models. However, the authors point out some limitations of using only NWP data. Terrain effects are not well enough accounted for on a kilometric scale, whereas future development of sub-kilometric scale NWP data might improve, for example, terrain effects on the incoming solar radiation. Combining NWP data with other data sources (e.g. observations, radar) might improve the forcing data, particularly the precipitation fields. Redistribution of snow due to wind is another difficult issue in mountainous areas. Running snow models with ensemble-based forecasts is a promising method that can be used to account for these uncertainties (Lafaysse et al., 2017; Vernay et al., 2015).
Weather forecasting models are presently evolving fast, and they include more and more detailed parameterization of land surface processes connected to snow and soil. SURFEX (Surface Externaliseé) (Masson et al., 2013) is an example of a land surface model, which can be run both inline as part of an atmospheric weather forecast model, for example AROME MetCoOp (Müller et al., 2017), or offline as a stand-alone model. At the Norwegian Meteorological Institute (MET Norway), the AROME MetCoOp model is run operationally to provide short-term weather forecasts covering large parts of the Nordic region (Müller et al., 2017).
Data assimilation methodologies for snow are presently likewise evolving fast; see e.g. Carrera et al. (2015). Snow analysis incorporating in situ observations, satellite data and estimates from NWP-driven physical snowpack models is an alternative to NWP-driven offline runs. In Norway, observed snow depth at weather stations are assimilated daily into the AROME MetCoOp model in order to improve the predicted surface air temperature. However, the MetCoOp model uses the most simple snow scheme (D95, Douville et al., 1995) of the three snow schemes implemented in the SURFEX model (Boone and Etchevers, 2001; Masson et al., 2013). The D95 scheme models the snowpack as a single layer, with two prognostic variables: SWE and snow density. For many of the above-mentioned applications, information on other snow properties such as internal layering as well as the density, grain size, temperature, wetness etc. of the different layers are of high interest. This is not provided by the D95 snow scheme.
In this study we therefore evaluate the performance of the SURFEX model, using the Crocus snow scheme (Vionnet et al., 2012) for Norwegian snow conditions. Crocus is the most advanced snow scheme implemented in the SURFEX model (Boone and Etchevers, 2001; Masson et al., 2013). SURFEX/Crocus has not previously been run in a gridded stand-alone version for regions in Norway (as a 2-D study). However, the model was tested earlier for single points (1-D study) with observations from weather stations and NWP data (Vikhamar-Schuler et al., 2011). Our study is carried out as part of several research projects within hydropower and flood forecasting. The domain was chosen to cover mountains in southern Norway and to include a cross section from west to east that crosses the watershed in this region. This domain includes catchment areas that are of high interest to hydropower companies.
The aim of our study is to test the effectiveness and benefits of different gridded forcing data sets as input to the SURFEX/Crocus model and validate the simulated snow amounts and snowmelt patterns in the selected domain. The originality of our work is linked to the unique combination of using both raw weather predictions, post-processed weather predictions and gridded observations, which we expect should provide an improved performance of the snow simulations, for example compared to using only raw weather predictions. Combining observations and NWP data for important weather variables (temperature, precipitation and precipitation phase), such as driving data for the snow simulations, should better represent the actual observed weather conditions. Experiments were performed by applying two different data sets from the winters 2014–2015 and 2015–2016 as forcing to the SURFEX/Crocus model: (1) predictions from the AROME MetCoOp model with a grid spacing of 2.5 km (Müller et al., 2017), where both the temperature and the wind data were improved by post-processing algorithms, and (2) gridded observations of precipitation and temperature (GridObs) with a grid spacing of 1 km (Lussana et al., 2018a, b). Both data sets have hourly temporal resolution and are discussed in detail in Sect. 2.3.
Although AROME-SURFEX/Crocus has previously been used over the southern European mountain chains in the French Alps (46∘ N, 9∘ E) and the French and Spanish Pyrenees (42∘ N, 1∘ E) (Quéno et al., 2016; Vionnet et al., 2016), neither of our two data sets described above have been used as forcing for SURFEX/Crocus for Norwegian mountains and lowland regions before. Our study area is located in northern Europe at 61∘ N, 8∘ E, which is at least 15∘ further north. According to the Köppen–Geiger climate classification system (Köppen, 1936), the climate in southern Norway is different from the Pyrenees, while the climate classes are partly the same for southern Norway and the Alps. However, these coarse climate classes generally account for average temperature and precipitation in an area and do not fully account for the differences in probability distribution functions describing the regional climate variability (e.g. precipitation intensities, frequencies and extremes) in these mountainous areas. A west–east transect crossing the mountain chain in southern Norway comprises a climatic transect from maritime to alpine to more continental climate. Snow conditions and stratigraphy vary regionally, e.g. as outlined by Sturm et al. (1995), who defined six snow classes, of which at least two are inside our domain (maritime and alpine). The SURFEX/Crocus snow model may therefore perform differently in individual regions. Our study contributes to a development which can produce new supplementary snow information (including snow stratigraphy) and thereby may contribute to the development of a future system for daily snow mapping. The performance of the SURFEX/Crocus model is also compared with three other snow models, including the seNorge model in a separate study by Skaugen et al. (2018), a study which also shows that there is no “best” snow model.
2.1 The SURFEX/Crocus model
The model used in this study is the detailed snowpack model Crocus (Brun et al., 1992; Vionnet et al., 2012) coupled with the ISBA land surface model within the SURFEX (Surface Externaliseé) interface (Masson et al., 2013). We applied the ISBA-DIF multi-layer soil scheme (Boone et al., 2000; Habets et al., 2003), which uses a diffusive approach for modelling the heat and moisture transport in the soil. The soil was divided into 14 layers, of which the thickness of the individual layers increases with the soil depth. The bottom depth of the lowest layer was 12 m. The HSWD (Harmonized World Soil Database) 1 km resolution database for soil texture (FAO/IIASA/ISRIC/ISS-CAS/JRC, 2012) was used for the soil properties.
The snowpack scheme Crocus models the physical properties of up to 50 dynamic layers within the snowpack, as well as the underlying ground. Once the snowpack reaches a threshold of 1 kg m−2 SWE, the fractional snow cover over a grid point is assumed to be 1. The SURFEX/Crocus model can be run in stand-alone (or offline) mode, or fully coupled to an atmospheric model.
For this study, the SURFEX/Crocus model was used in an offline mode, on a 0.01∘ grid (approximately 1 km), with a 5 min internal time step and output every hour. The orography was taken from GTOPO30 global digital elevation model (DEM) from the US Geological Survey, which has a grid spacing of ∼ 1 km. The transport of snow by wind is not simulated. The SURFEX/Crocus model was run for two winter seasons, from 1 September 2014 to 31 August 2016. These dates were chosen because the hydrological year starts on 1 September, and at that time there is normally no snow in the mountains. In this study, we start a new simulation on 1 September, with no snow present, and with default values for soil properties for both 2014–2015 and 2015–2016. The default soil temperature is 11.9 ∘C for the uppermost surface soil layer for 0 m a.s.l. (sea-level height). The soil temperature is reduced with increasing terrain elevation using a lapse rate of 0.65 ∘C per 100 m, leading to a surface soil temperature of 1.3 ∘C at 1000 m a.s.l. We estimate these surface soil values to be representative of the September climate in our study area. Higher temperatures in the deepest soil layers may, however, represent an uncertain heat contribution for the snow modelling. This effect should be similar for all the experiments, since the initialization is the same.
2.2 The study area
Figure 1 shows the domain over which the SURFEX/Crocus model was run and the elevation over the model domain. The domain covers nearly 20 000 km2 (111 × 175 km) and contains 100 × 330 grid points. As mentioned in the introduction, the study area was chosen to cover the mountains in southern Norway and to include a cross section from west to east that crosses the watershed in this region as well as several catchment areas that are of interest to hydropower companies. The domain covers elevations from 0 m a.s.l. along fjords up to the highest mountain in Norway (2468 m a.s.l.). Therefore, the area includes different vegetation zones, ranging from high mountains above the treeline, sparsely forested and densely forested areas. This makes it a challenging area for snow modelling.
Due to the watershed and the prevailing weather patterns, there is a large gradient in precipitation amount over the domain. The far western parts of the domain receive on average around 1500 mm of precipitation during a winter season, while the eastern parts only receive 100–300 mm (Hanssen-Bauer et al., 2015). The western part of the domain has a maritime climate, while the eastern part has a more inland climate, which means the average temperature during winter is higher at the western part of the domain (around or just below 0 ∘C) compared to the eastern side (around −10 ∘C) (Hanssen-Bauer et al., 2015). This means the gradient in average snowfall amount is not as large as the gradient in precipitation amount, but the western part of the domain still receives significantly more snow than the eastern part (Hanssen-Bauer et al., 2015).
2.3 Forcing data sets
The SURFEX/Crocus model requires atmospheric forcing. For this study, we have used two different sets of forcing data. Table 1 shows an overview of which variables the SURFEX/Crocus model requires and the different sources used in the two experiments: (1) AROME-Crocus and (2) GridObs-Crocus. AROME-Crocus uses both raw and post-processed forecasts from the AROME MetCoOp model (described below in Sect. 2.3.1), while GridObs-Crocus uses a combination of gridded observations of precipitation and temperature, described in Sect. 2.3.2. All forcing data have hourly temporal resolution.
2.3.1 Numerical weather forecasts (AROME-MetCoOp)
AROME MetCoOp is a high-resolution, non-hydrostatic, convective-scale weather prediction model operated by a bilateral cooperative effort (Meteorological Cooperation on Operational Numerical Weather Prediction; MetCoOp) between the Norwegian Meteorological Institute and the Swedish Meteorological and Hydrological Institute (Müller et al., 2017), operational since March 2014. The core of the model is based on the convection-permitting Applications of Research to Operations at Mesoscale (AROME) model developed by Météo-France (Seity et al., 2011). It has been modified and updated to suit advanced high-resolution weather forecasts over the Nordic regions; see Müller et al. (2017) for details. The horizontal grid spacing is 2.5 km and the domain covers the Nordic countries. The atmosphere is divided into 65 vertical levels, with the first level at approximately 12.5 m height. Atmosphere–surface interactions and surface–soil processes are described by SURFEX (Masson et al., 2013). The fluxes computed by SURFEX at the atmosphere–surface interface serve as the lower boundary conditions for the atmosphere within AROME MetCoOp. All surface processes are treated as one-dimensional vertical processes.
AROME MetCoOp operates with a 3 h update cycling, where initial fields of atmospheric and land surface variables are corrected with observations through data assimilation. Observations of air temperature, relative air humidity and snow depth are used in the surface analysis (Müller et al., 2017). At every main cycle (00:00, 06:00, 12:00 and 18:00 UTC) a 66 h forecast is produced. Forcing for our study is taken from the four main cycles, with successive 3–8 h lead time (0–8 h lead time for the 00:00 UTC cycle and 3–5 h lead time for the 18:00 UTC cycle) forecasts combined into a forcing file for each day. These lead times were chosen to avoid the first hours of a cycle when the model might have spin-up issues and to make use of all available cycles with the shortest possible lead time (model error increases with lead time; see for example Homleid and Tveter, 2016).
For temperature and wind speed we used statistically post-processed AROME-MetCoOp forecasts to force SURFEX/Crocus, described by Køltzow (2017). These post-processed weather variables are produced operationally by MET Norway for the weather forecast website YR (https://www.yr.no/, last access: 18 June 2018). The temperature grid has a spatial resolution of 500 m and is produced using a Kalman filter correction at observation stations (Homleid, 1995). Horizontal interpolation is carried out using decreasing weights with increasing distance from the station. The temperature is further corrected for terrain elevation, which also takes into account vertical temperature profiles in inversion situations in wintertime. The AROME-MetCoOp wind speed was statistically post-processed to represent the maximum wind speed at 10 m during the last hour. In addition, correction factors are applied to the wind speed depending on wind direction and region (Køltzow, 2017). The other variables from the raw AROME-MetCoOp 2.5 km forecasts were interpolated to 1 km spatial resolution using bilinear interpolation in order to combine the meteorological forecasts with the gridded observations (with a spatial resolution of 1 km) and to run the SURFEX/Crocus model with 1 km grid spacing. The 500 m post-processed AROME-MetCoOp temperature (Køltzow, 2017) was also interpolated by a bilinear method to 1 km resolution. The spatial interpolation was carried out using the File Interpolation, Manipulation and EXtraction library (http://fimex.met.no, last access: 18 June 2018).
SURFEX/Crocus requires a separate snowfall and rainfall rate. A threshold temperature of +0.5 ∘C was applied to determine snowfall or rainfall. This threshold temperature is commonly used for hydrological purposes in Norway (see, for example, Skaugen, 1998). The post-processed AROME-MetCoOp temperature was used to compute precipitation phase based on the total precipitation predicted from the AROME-MetCoOp model. Correct precipitation phase estimation is crucial for good snow simulations. An additional test was carried out on estimating precipitation phase using the raw AROME-MetCoOp snowfall and rainfall at 2.5 km resolution. AROME's own microphysics should provide good precipitation phase estimates, but these are most representative at 2.5 km spatial resolution on the model's own terrain height. In our study area, the terrain variability is very large, particularly in the western regions where the terrain often rises from 0 m a.s.l. at the fjords to more than 1000 m a.s.l. over very short distances. In these kind of areas, terrain-adjusted precipitation phase determination is necessary. The impact on the snow simulations from using these two different ways of estimating precipitation phase are discussed in Sect. 4.2.
2.3.2 Gridded observations (GridObs)
In an earlier study by Vikhamar-Schuler et al. (2011), it was shown that snow modelling with the SURFEX/Crocus model has the highest sensitivity to the temperature and precipitation input data sets. The best results were obtained when the model was forced with observations of temperature and precipitation, while replacing other input parameters with meteorological forecast data did not notably increase errors.
Hourly gridded observations of temperature and precipitation are available on a 1 km grid over Norway. This data set uses all measurements available in MET Norway's Climate database (Frost, 2018). The station distribution is uneven, with more stations in the southern part of Norway and a sparser network in the north and in the mountains. There is a low elevation bias, because most stations are located at lower elevations (e.g. valley bottoms) and few stations are located above 1000 m a.s.l. (Lussana et al., 2018a). The hourly precipitation values have been obtained by using a two-step procedure. The spatial interpolation method described by Lussana et al. (2018a) has been applied independently to daily and hourly precipitation totals. This method is built on classical methods (such as optimal interpolation and successive-correction schemes) and (spatial) scale-separation using geographical coordinates and elevation as complimentary information in the interpolation. It is based on iterating a statistical interpolation scheme over a decreasing sequence of spatial scales, from synoptic to kilometre scale. Further details can be found in Lussana et al. (2018a). The daily precipitation totals have been disaggregated to an hourly time with a procedure similar to the one described by Vormoor and Skaugen (2013). The two-step procedure has been implemented so that the final hourly product can benefit from the more accurate daily quantitative estimates that are based on a denser network of stations if compared to the hourly ones. The method used to obtain hourly temperature values is described in Lussana et al. (2018b), while the resulting temperature data set is described and evaluated in Lussana et al. (2016).
The resulting gridded temperature data set can be regarded as an unbiased estimate of the true temperature both at grid points and at station locations. Only for the most extreme negative values (temperatures below −30 ∘C) is there a systematic warm bias of about 1 ∘C (Lussana et al., 2016, 2018b). For precipitation, Lussana et al. (2018a) found that the precision of the estimates (at grid points) is about ±20 %, but there is a systematic underestimation of precipitation in data-sparse areas and for intense precipitation.
The first version of these gridded data sets (called seNorge v1.0) included a correction factor for precipitation at individual stations due to undercatch (Lussana et al., 2018a). However, these correction factors were evaluated and found to be unrealistic in the mountains, producing too much precipitation (Saloranta, 2012). In the next version of these gridded data sets (seNorge v.2.0), all undercatch correction factors were removed, and interpolation was based on uncorrected precipitation. For the same reason, the hourly gridded data set used in this study does not contain any correction for precipitation undercatch.
The number of stations included in the gridded data set varies over time. The numbers of stations within the SURFEX/Crocus domain are 20–30 stations for hourly precipitation, 90–100 stations for daily precipitation and 70–100 stations for temperature. Stations just outside the domain are included in this estimate as they are used in the interpolation and are therefore part of the gridded data set used in this study.
Snowfall and rainfall rates were estimated assuming rain/snow separation at +0.5 ∘C (using the gridded observations of temperature available on the same grid), the same threshold as used in the AROME-MetCoOp forcing data set.
2.4 Validation data set
We use two different data sets to validate the results from both experiments: point observations of snow depth and snow cover maps derived from MODIS satellite images.
2.4.1 Snow depth observations
Observations of daily snow depth from 30 stations were selected for verification of the model results (see Fig. 1 for their locations within the domain). Nearly all stations (25 out of 30) are official meteorological stations run by the Norwegian Meteorological Institute, while a few stations are owned by other institutions (municipalities, energy producing companies and Bane Nor, the state-owned company responsible for the Norwegian national railway infrastructure). Data from all stations are freely available from the climate database of the Norwegian Meteorological Institute (Frost, 2018). All stations measure daily snow depth, nearly all (29 out of 30) measure precipitation, and nine stations also measure air temperature. The stations were selected based on the availability of snow depth observations between 1 September 2014 and 31 August 2016. The locations of the stations are reasonably well distributed over the domain (see Fig. 1) and their elevations range between 14 and 1162 metres above sea level. Figure 2 shows the elevation distribution of all stations used in this study. Along with the distribution of elevations of grid points in the SURFEX domain, Fig. 2 shows a typical issue of low elevation bias in the observing network of the Norwegian Meteorological Institute, also illustrated in Lussana et al. (2018a). Stations are located at elevations that seldom exceed 1000 m.
A nearest-neighbour method was used to evaluate the SURFEX/Crocus experiments with the surface snow depth observations. In a domain with deep valleys and high mountains, it is difficult to match the exact elevation of the weather stations with the nearest grid point in the SURFEX/Crocus experiments. As there were only 30 stations with high-quality snow depth observations in the domain, it was decided not to filter out stations based on these elevation differences. The influence of elevation differences in the evaluation is discussed in Sect. 4.1.
Daily snow depth observations taken at 06:00 UTC have been used for direct comparison to snow depth from the SURFEX/Crocus simulations. The observations were also used to calculate the start, length and end of the snow season to compare against model results. The length of the snow season is defined as the number of days with more than 5 cm of snow during a year. The 5 cm threshold was also used by Vionnet et al. (2016), although they used continuous snow on the ground as an additional condition. The start of the snow season is defined as the first day with more than 5 cm of snow and the end of the snow season as the day after the last day with more than 5 cm of snow.
2.4.2 MODIS snow cover images
MODIS (Moderate Resolution Imaging Spectroradiometer; http://modis.gsfc.nasa.gov/, last access: 18 June 2018) snow cover images (Hall and Riggs, 2007; Klein and Stroeve, 2002) with a resolution of 500 m were available and processed for the melt season of the 2014–2015 winter. The same method as described by Lussana et al. (2018a) and Skaugen et al. (2018) was used to obtain estimates of the daily snow cover extent over the domain: the MODIS images were converted to snow-covered area (SCA) on a scale from 0–100 % coverage using a method based on the Norwegian linear reflectance to snow cover algorithm (NLR) (Solberg et al., 2006). The input to the NLR algorithm is the normalized difference snow index signal (NDSI-signal) (Salomonson and Appel, 2004).
The MODIS images were used for visual and quantitative comparison of the snowmelt pattern from satellite images and from both SURFEX/Crocus experiments. For this purpose, four dates with cloud-free conditions were selected from the melt season: 15 March 2015, 20 April 2015, 15 May 2015 and 04 July 2015.
3.1 Snow depth
A density scatter plot of daily observed and simulated snow depth for both experiments and the two winter seasons 2014–2015 and 2015–2016 is shown in Fig. 3. Zero snow depth pairs were excluded. GridObs-Crocus is in reasonably good agreement with the observations (R2= 0.78), although there are cases of over- and underestimation of around 100 cm, while AROME-Crocus shows significantly more variability and overestimation of snow depth (R2= 0.52). To investigate the snow depth at individual stations over a range of station altitudes in more detail, Fig. 4 shows snow depth plots for six locations: two are located below 400 m, two are located between 500 and 900 m, and two are above 900 m (which in our study area means they are located above the treeline). For the location of these six stations within the domain, see Fig. 1 in which they are indicated with blue dots. These six stations show that AROME-Crocus overestimates the snow depth for the highest altitudes (Bakko i Hol and Midtstova, which are both situated above 900 m a.s.l.), while it underestimates the snow depth for the lowest stations (Balestrand and Haukedal). The two lowest stations are located in the western part of the study area where terrain gradients are very steep.
The snow depth from GridObs-Crocus is closer to the observed snow depth, but at times it underestimates the snow depth (most notably for the first winter season at Haukedal, 329 m a.s.l.; and Midtstova, 1162 m a.s.l.). Episodes in which the snow depth decreases during the winter season (apart from snowmelt in spring) are not always well captured by the SURFEX/Crocus experiments, and this issue is partly responsible for the overestimation of snow depth.
The results for Hemsedal II (604 m a.s.l.; see Fig. 4) are of particular interest, as this is the only station measuring snow depth but not precipitation (and therefore not part of the gridded observation data set used as input for GridObs-Crocus). GridObs-Crocus overestimates the snow depth at Hemsedal II but slightly less than AROME-Crocus does. The biases in snow depth at Hemsedal II for the two seasons combined are 25 cm for GridObs-Crocus and 29 cm for AROME-Crocus. When compared to the biases (7 cm for GridObs-Crocus and 20 cm for AROME-Crocus) for all stations for the two seasons combined, it shows that Hemsedal II performs slightly worse than most stations in AROME-Crocus. For GridObs-Crocus, the bias at Hemsedal II is significantly larger than for most other stations. The fact that GridObs-Crocus outperforms AROME-Crocus even at a station that is not part of the gridded observation data set is interesting.
The strong overestimation at Midtstova (see Fig. 4) by AROME-Crocus can be explained by the fact that Midtstova is located in an area with systematic and relatively large overestimation of precipitation in AROME-MetCoOp. In addition, AROME-MetCoOp underestimates the temperature by about 2∘ in this area during winter. This can be seen in verification reports of the AROME MetCoOp model, for example in Homleid and Tveter (2016). In the forcing data for Midtstova we find a bias of −1.5∘ for AROME-Crocus, compared to −0.8∘ for GridObs-Crocus. This bias is larger than the overall bias for all nine stations measuring temperature: −0.5∘ for AROME-Crocus and −0.2∘ for GridObs-Crocus. During the snow accumulation season the temperature at Midtstova is mostly well below freezing level. There are a few episodes each winter with temperatures just above zero, where the underestimated temperature in AROME-MetCoOp means the precipitation during those episodes comes as snow instead of rain, but these do not add up to large amounts. Midtstova is also a high-mountain station, which is very exposed to strong wind. Redistribution of snow due to wind is not captured in the SURFEX/Crocus model. GridObs-Crocus shows much more realistic results for Midtstova, although there is an underestimation during the first part of the 2014–2015 winter. From 27 October 2014 to 26 January 2015, the precipitation sensor at Midtstova was out of order, and the forcing from GridObs-Crocus for Midtstova will therefore be represented by interpolated values from surrounding stations, which might explain the underestimation. An evaluation of the precipitation forcing data for AROME-Crocus and GridObs-Crocus for the 30 weather stations reveals that the AROME-Crocus forcing has about 40 % more rain and 20 % more snow compared to GridObs-Crocus. The differences are largest for the stations above 800 m, which often receive about 50 % more snow. This can clearly be seen for Bakko i Hol and Midtstova in Fig. 4.
Table 2 summarizes the bias over all stations for the two winter seasons (1 September 2014–31 August 2015 and 1 September 2015–31 August 2016). The bias was calculated as the mean of the differences between simulated and observed snow depth, and only for the days on which there is snow present in the observations or at least one of the experiments. GridObs-Crocus shows a significantly smaller bias (4 and 9 cm) compared to AROME-Crocus (23 and 17 cm). The maximum observed snow depth is on average 112 cm for 2014–2015 and 88 cm for 2015–2016. GridObs-Crocus shows a very small bias (0 and 9 cm respectively), while AROME-Crocus overestimates the mean maximum snow depth by 17–18 cm. Table 3 summarizes the RMSE over all stations for the two winter seasons. The RMSE values are larger for AROME-Crocus (compared to GridObs-Crocus) for nearly all variables, except for the date of maximum snow depth for 2015–2016.
As snow depth accumulates over the winter season, a missed (or under/over estimated) snow event can influence the remainder of the season. It can therefore be useful to look at daily snow depth variations instead, as was also done by Quéno et al. (2016) and Schirmer and Jamieson (2015). Figure 5 shows the categorical frequency distribution of the daily change in snow depth for six accumulation categories, five reduction categories and one category centred around zero accumulation, on a logarithmic scale. The first two accumulation categories (up to 10 cm) are overestimated in both GridObs-Crocus and AROME-Crocus. The strongest observed accumulation category (> 40 cm) as well as the strongest reduction category (< −20 cm) are not represented in either of the SURFEX/Crocus experiments.
SURFEX/Crocus in stand-alone mode does not account for wind-induced snow redistribution, which can largely contribute to strong decreases in snow depth. Figure 4 showed that episodes of a reduction in snow depth (not including the snowmelt at the end of the season) were not always captured well by the models, and it could be that blowing snow is the cause of this. Following Quéno et al. (2016), two diagnostics have been applied to look into this issue: blowing snow days and melting snow days. Blowing snow days are defined as days during which the wind speed (at a height of 10 m) (between 06:00 and 06:00 UTC, as this is when snow depth measurements are made) exceeds 8 m s−1 during the past 24 hours, while the snow surface temperature is below 0 ∘C (since only dry snow can be blown away). The temperature of the snow surface is taken from the SURFEX/Crocus output for 12:00 UTC each day. The wind speed is taken from AROME-MetCoOp, which is used as forcing in both SURFEX/Crocus experiments. The modelled wind speed is used because only 6 out of 30 stations used in this study observe wind speed. When comparing the forecasted maximum wind speed from AROME-MetCoOp with the observed maximum wind speed from these 6 stations, we find a slight overestimation by AROME-MetCoOp (a bias of 0.3 m s−1). Blowing snow days and non-blowing snow days are correctly identified in 94 % of all days, with a hit rate of 0.86 and a false alarm rate of 0.04. This shows that the modelled wind speed can be used to determine blowing snow days. The wind threshold of 8 m s−1 for dry snow transport is taken from Li and Pomeroy (1997). Figure 6 shows the cumulated number of the daily changes in snow depth for five categories of decreasing snow depth for blowing snow days and for all days on which the snow depth decreases, for the observations and for GridObs-Crocus, as well as the percentage of snow depth loss due to blowing snow. The cumulated amount of snow reduction is underestimated for nearly all categories. For the strongest decreasing rate (more than 20 cm), the observations indicate that 51 % of the decrease in snow is caused by blowing snow. This category is not represented by GridObs-Crocus. For GridObs-Crocus, blowing snow days only contribute to the smallest reduction categories. In total (over all categories), blowing snow days contribute to 17 % of the cumulated decrease in snow depth in the observations, while this amounts to 10 % in GridObs-Crocus.
Melting snow days are defined as days when the simulated surface temperature of the snow is 0 ∘C. Figure 7 is similar to Fig. 6, but for melting snow days. For GridObs-Crocus, melting snow is the main factor responsible for contributing to a decrease in snow depth. The largest reduction category is not represented by GridObs-Crocus, but for the other categories, melting snow is responsible for 57–100 % of the decrease in snow depth. This is not surprising as SURFEX/Crocus does not represent blowing snow, so a decrease in snow depth is caused by either snowmelt or other processes such as snow compaction. The cumulated daily changes in snow depth for melting snow days as well as all days with a decrease in snow depth are underestimated by GridObs-Crocus for all categories except the smallest one (less than 5 cm loss in snow depth). This shows there is a general underestimation of snow ablation, as well as an underestimation of snowmelt in GridObs-Crocus. The same goes for AROME-Crocus (not shown).
SURFEX/Crocus does have an option to run with sublimation in case of snowdrift. This option has been tested for two stations from Fig. 4: Midtstova and Hemsedal II. In this experiment, SURFEX/Crocus was run twice in 1-D mode for these two locations: one experiment with identical settings to AROME-Crocus, and one nearly identical to the exception of the option for sublimation in case of snowdrift (AROME-Crocus + BS). The results are shown in Fig. 8. For both locations, the snow depth in AROME-Crocus + BS is decreased, as expected. For both locations, this is an improvement. The bias in AROME-Crocus was +13 cm for Hemsedal II, which improved to +5 cm in AROME-Crocus + BS. For Midtstova, AROME-Crocus significantly overestimates the snow depth (bias: +104 cm). This is improved in the AROME-Crocus + BS experiment (+83 cm), although the overestimation is still very large. The length of the snow season is reduced by a few days for both stations and both years, similarly to the results found in Brun et al. (2013). However, for Midtstova, using blowing snow sublimation does not improve the AROME-Crocus experiment to the extent that it performs as well as GridObs-Crocus.
3.2 Characteristics of the snow season
Statistics for the snow season duration are shown for the two snow seasons 2014–2015 and 2015–2016 in Table 2. The length of the snow season is defined as the number of days with more than 5 cm snow during a season. For GridObs-Crocus, the length of the snow season is overestimated by 8–11 days (see Table 2), while AROME-Crocus overestimates the length of the snow season by 17 days in 2014–2015 and underestimates it by only 2 days in 2015–2016. The same positive bias found for AROME-Crocus in the 2014–2015 season was found by Vionnet et al. (2016). One possible explanation of this bias is the fact that the SURFEX/Crocus model assumes a uniform snow cover from the moment snow is present on the ground and therefore shows less variability in snow cover compared to the observations. In the observations, there is often a period in which the snow cover fluctuates - for example thinning to below 5 cm after the first snow has fallen and before a continuous snow cover has been established for the winter season. AROME-Crocus predicts the length of the 2015–2016 season well.
The start of the snow season is defined as the first day with more than 5 cm of snow and the end of the snow season as the day after the last day with more than 5 cm of snow. A negative bias in the start of the snow season means a too-early start, while a positive bias means a too-late start to the snow season. GridObs-Crocus has a bias of only 2 days (negative for the first winter and positive for the second winter) for the start of the snow season, while the snow season starts up to 12 days too-late in AROME-Crocus during the second year (the first year has a bias of −2 days). The season ends on average in late April or early May. In GridObs-Crocus, the season ends 3 days early in the first year and 2 days late in the second year. In AROME-Crocus, this is 7 days late and 3 days early respectively. The observed maximum snow depth occurs on average at the end of January in the first year and in late February in the second year. Both experiments show a later date for the maximum snow depth.
Figure 9 shows the distribution of the bias at the start and end of the snow season as well as the date of maximum snow depth, for all 30 stations and for two winter seasons. Most stations show a bias near zero (between −5 and +5 days) for the start of the snow season. In AROME-Crocus, the snow season sometimes starts much too late. The bias for the end of the snow season shows that GridObs-Crocus often ends the snow season too early, while AROME-Crocus tends to end the season too late. The bias for the date of the maximum snow depth of the season is mostly around zero for most stations and both models, but there are some outliers, especially towards the strong positive bias. This is due to stations such as Midtstova in Fig. 4, where the maximum observed snow depth occurs rather early in the season, while both experiments show a maximum much later in the season.
3.3 Snow cover pattern
Figure 10 shows the spatial pattern of snow cover over the SURFEX/Crocus domain compared to MODIS data over the same area. The snow-covered area is shown on different dates in the snowmelt season: 15 March, 20 April, 15 May and 4 July 2015. On 15 March 2015, nearly the whole area is covered with snow. The only exceptions are areas right besides the fjords (white areas) in the west (well captured by both experiments) and at the bottom of valleys in the east (not captured by the SURFEX/Crocus experiments). On 20 April 2015, the snow has clearly started to melt in the valleys to the east. This is captured well by AROME-Crocus, while GridObs-Crocus shows too little snow around the valleys in the south-east of the domain. By 15 May 2015, a lot of snow had disappeared in the eastern part of the domain, while the western part has not changed much from the previous month. The average date for the end of the snow season for all the 30 weather stations for the 2014–2015 season was 2 May 2015 (see Table 2), but the dates of the end of the snow season for individual stations range from 18 February (Fresvik, 32 m a.s.l.) to 5 July (Midtstova, 1162 m a.s.l.). Again, AROME-Crocus captures the snow cover pattern better than GridObs-Crocus. By 4 July 2015, the snow cover is limited to areas with higher elevation. AROME-Crocus captures the spatial pattern of snow cover very well. In GridObs-Crocus, nearly all snow had melted, and the snow-covered area is underestimated. Earlier it was shown (in Fig. 9) that GridObs-Crocus has a negative bias (too early) for the end of the snow season for the 30 snow depth stations, while AROME-Crocus has a positive bias (too late). As discussed previously, the differences between the snowfall amounts from the two precipitation forcing data sets are largest for the highest parts of the domain, where AROME-Crocus receives about 50 % more snow compared to GridObs-Crocus. This explains why the differences between GridObs-Crocus and AROME-Crocus in Fig. 10 are also largest in this area (especially by the end of the snow season on 4 July).
Table 4 shows the Jaccard indices for the images from Fig. 10. The Jaccard index was also used by Quéno et al. (2016). It is a similarity index applied to the snow cover images which were remapped onto the same grid (which means that the snow cover from the MODIS images used to calculate the Jaccard index has a lower resolution than the one shown in Fig. 10). The Jaccard index is calculated as , where X and Y are the simulated and observed snow cover, respectively. The number of grid points that are snow-covered in both SURFEX/Crocus and in the MODIS image is divided by the total amount of snow-covered grid points (in either SURFEX/Crocus or MODIS). When the Jaccard index equals 1, there is a perfect match between snow-covered grid points, and when the Jaccard index equals 0, there is no match at all. Table 4 shows that AROME-Crocus consistently has higher Jaccard indices compared to GridObs-Crocus. The indices decrease (for both experiments) during the melt season.
Although both experiments are capable of simulating the snowpack over the two winter seasons, the two simulations provide different results regarding the snow depth and the spatial snow-covered area. There is an overestimation of snow depth in the AROME-Crocus experiment, even though the snow-covered area in the melt season is better represented by this experiment. When using gridded observations (GridObs-Crocus), the simulation of snow depth is significantly improved, while the spatial distribution of the snow cover is highly underestimated, particularly late in the snowmelt season. There is an underestimation of snow ablation in both experiments, which is due to a combination of the absence of wind-induced erosion of snow and underestimation of snowmelt in SURFEX/Crocus and biases in the forcing data. Possible causes for these different results are further discussed below, by focusing on the quality of the model validation, the forcing data set and the snowpack model.
4.1 Quality of the model validation
The model validation was carried out using both snow measurements at individual weather stations and MODIS satellite images. Using several data sources to validate simulations is important, as these two sources supplement each other. Stations provide point validations with daily time series, while the satellite images provide images of the snow cover for an entire area for cloud-free days. Even though the GridObs-Crocus simulation provides reasonable results at individual stations, the MODIS images show that the snow cover disappears too fast, particularly late in the snowmelt season. This may indicate that the gridded interpolated observations of temperature and precipitation (the forcing data) are not representing the terrain variability in the study area sufficiently well. The western region is dominated by terrain with steep gradients, which requires a higher density of weather stations representing the full range of terrain elevations compared to smooth landscape areas. As described in Sect. 2.2, there is a low elevation bias in the national observational network with too few stations in areas above 900 m a.s.l. This increases the uncertainty in the precipitation and temperature estimates in the mountainous regions. The quality of the gridded data is obviously highest for locations closest to the stations, providing better results in those areas. An underestimation of the snow depth at high elevations would explain the underestimated snow cover during the melt season.
Nearly all (29) of the 30 stations that measure snow depth also measure precipitation (7 measure hourly precipitation, while 22 measure daily precipitation), which means the observed precipitation from these stations is used in the gridded observation data set used to force GridObs-Crocus. Only 9 out of the 30 stations also measure temperature. Although precipitation is not directly related to snow depth and temperature also plays an important role, it could still be argued that the GridObs-Crocus results are best in the locations of the observations that are included in the gridded data set used to force SURFEX/Crocus. The only station that was not part of the gridded precipitation data set is Hemsedal II. The bias for GridObs-Crocus at Hemsedal II is 25 cm, which is larger than the overall bias for all stations (7 cm), but the RMSE is about the same (27 cm for Hemsedal II and 28 cm for all stations). Although this shows that the performance of a station not included in the gridded precipitation data set is about the same as the performance of stations that are part of this data set, one station is not enough to draw conclusions about the entire domain.
The representativity of a station location is sensitive to the terrain variability. The orography used in the SURFEX/Crocus experiments has a resolution of 1 km, leading to differences between the actual station height and the average height used for the centre of the 1 km grid cell in SURFEX/Crocus. Figure 11 shows the distribution of those differences for all 30 stations. The average bias is 79 m. Most stations are placed at higher elevations in the model than their actual elevation, but it should also be kept in mind that a grid point in SURFEX/Crocus describes a larger area (and range of elevation) compared to the actual observations. Especially in the mountainous region in the west of the domain (see Fig. 1), which has high mountains, steep slopes and deep valleys, there may be large differences in height within a distance range of 1 km. In Sect. 4.2 the sensitivity of terrain effects on precipitation phase computations is discussed.
4.2 Quality of the forcing data sets
Raleigh et al. (2015) showed that snow simulations are more sensitive to biases in forcing data than random errors, and that precipitation bias is the most important factor. There is a negative bias in the gridded precipitation used in GridObs-Crocus (Lussana et al., 2018a), especially for data-sparse areas (e.g. high-mountainous areas) and for intense precipitation. Missing episodes of intense snowfall would explain part of the underestimation of the snow depth in GridObs-Crocus. There are plans to improve the gridded observations of precipitation by adjusting the solid precipitation to account for the wind undercatch and by post-processing the predicted precipitation fields to adjust for bias (Lussana et al., 2018a). Forecasts from the AROME-MetCoOp model are known to overestimate the occurrence of precipitation events of less than 10 mm (Müller et al., 2017), and Fig. 5 showed that AROME-Crocus overestimated daily changes in snow depth up to 15 cm. An evaluation of the accumulated snowfall from the two forcing data sets showed that snowfall amounts in AROME-Crocus are about 20 % higher compared to GridObs-Crocus and even more (∼ 50 %) for stations at altitudes above 800 m a.s.l.
To test the sensitivity of terrain effects on precipitation phase and on the snow depth simulations in regions with steep topography, we compared two methods for determining precipitation phase from the AROME-MetCoOp forecasts. First, precipitation phase was determined as already presented in our AROME-Crocus experiment (by using the post-processed air temperature as a threshold temperature; see Table 1), and second, using the raw AROME-MetCoOp snowfall and rainfall forecasts computed by the atmospheric model's own microphysics. These 2.5 km spatial resolution forecasts were then interpolated to 1 km using standard bilinear interpolation, also described in Sect. 2.3.1. The comparison revealed that high-resolution terrain data improved the results by reducing the amount of snow, particularly in low-altitude areas near steep terrain gradients. The bias in snow depth improved from +42 cm (using raw AROME-MetCoOp precipitation) to +20 cm (using the terrain-adjusted AROME-MetCoOp precipitation), and the RMSE improved from 68 cm (raw) to 56 cm (terrain-adjusted). Figure 12 shows the fraction of snowfall compared to the total precipitation for both methods for the first winter season. It is clear that using the terrain-adjusted AROME precipitation results in a better representation of the terrain in the domain, as many features are lost in the coarser-resolution raw AROME precipitation. Figure 13 illustrates the differences in the computed snowfall from the two methods as a function of elevation and longitude. The west–east transect from steep terrain (including fjords and mountains) in the west to smoother terrain in the east is clearly illustrated, with largest differences in the western part of the domain. In these areas a more realistic precipitation phase (rainfall) was more frequently computed than when using the raw AROME-MetCoOp precipitation phase. This emphasizes the importance of terrain-adjusting the forcing data from NWP models for obtaining more correct precipitation phases.
Precipitation phase in both AROME-Crocus and GridObs-Crocus was determined using a fixed threshold temperature of 0.5 ∘C to distinguish between rainfall and snowfall. This simplification represents an uncertainty which could result in some actual snow events characterized as rainfall and, to a lesser extent, the other way around. Future studies could focus on studying alternative ways to determine the precipitation phase by better exploiting the microphysics of the NWP model. Generally, our experiments show that the highest chance of obtaining good snow simulations for larger regions lies in improving the forcing data, and particularly in improving the raw NWP data by post-processing techniques. The use of gridded interpolated observations alone, as was tested in GridObs-Crocus, displays some limits for snow cover mapping over larger regions. There is a high chance that NWP forcing can be improved by combining different data sources, e.g. assimilation of various observations (weather station data, precipitation radars, etc.).
Sauter and Obleitner (2015) investigated the sensitivity of SURFEX/Crocus snowpack modelling on Svalbard (Arctic Norway) to input parameters and found that, for higher elevations (in the accumulation zone), precipitation and radiation are key factors in the evolution of the snowpack and contribute the most to the model uncertainty. At lower elevations, precipitation was less important but factors such as wind speed or surface roughness increased in importance. Quéno et al. (2017) used satellite products of incoming solar and long-wave radiation to force the SURFEX/Crocus model; however they concluded that improved meteorological forcing does not always lead to more accurate snowpack simulations due to error compensations within the atmospheric forcing and the snowpack model.
4.3 Quality of the snowpack model
The SURFEX/Crocus model assumes a uniform snow cover when SWE reaches the relatively low threshold of 1 kg m−2. The SURFEX/Crocus model was originally developed for use in high alpine regions, where there is not a lot of vegetation. In those areas, the assumption of the uniform snow cover is realistic, as there is no interaction with vegetation, but for areas covered with forest and closer to sea level this could lead to an overestimation of the snow cover. When the snow cover is overestimated, the albedo will be too high and this will slow down the snowmelt at the end of the season. This might explain the underestimated snowmelt in both experiments.
The SURFEX/Crocus model grid is a collection of independent grid points with no transport of snow or other variables between grid points. It is therefore not possible to simulate the redistribution of snow by wind. It can be argued that, with a resolution of approximately 1 km, the drifting snow would in any case be redistributed within the area of a grid point and not transported to neighbouring grid points. Vionnet et al. (2014) showed that, for an explicit simulation of wind-induced snow transport, a spatial resolution of less than 50 m is required. This is currently not a feasible option for snowpack simulations over larger domains. There is an option in the SURFEX/Crocus model to calculate the rate of sublimation in case of snowdrift, which results in a loss of snow. This option was tested for two stations in this study, using the AROME-Crocus forcing data set. As expected, this resulted in a decrease in snow depth and a decrease in season length. This is an improvement, as it reduces the overestimation by AROME-Crocus but only to some degree.
Figure 5 showed that both SURFEX/Crocus experiments underestimate the melting of snow, which is further supplemented by Fig. 7 for the GridObs-Crocus experiment. Underestimated melting was also found by Quéno et al. (2016, 2017) and Vionnet et al. (2016), and complementary studies are needed to investigate the cause of this issue.
Lafaysse et al. (2017) developed an ensemble snowpack model using SURFEX/Crocus called ESCROC (Ensemble System Crocus) to address modelling errors. They found that by using optimal members they were able to explain more than half of the simulation errors, and those ensembles have a significantly better predictive power than the classical deterministic approach. For future work, it would be interesting to use ESCROC and investigate the effect of different physical settings of SURFEX/Crocus. In addition, since November 2016, AROME-MetCoOp is run as an ensemble with 10 members, called MEPS (MetCoOp Ensemble Prediction System). This means that an ensemble of meteorological forcing is another possible direction for future work. Vernay et al. (2015) used the 35 members of the ensemble prediction system based on the French NWP model ARPEGE as forcing to the SURFEX/Crocus model. The results indicated that accounting for the uncertainty in meteorological forecast significantly improves the skill and the usefulness of the model chain.
In this study we have evaluated the performance of the SURFEX/Crocus snow model for a region in southern Norway covering both steep terrain gradients with fjords and high-mountain areas in the western parts as well as smoother terrain in the eastern parts. The experiments tested different types and combinations of forcing data (raw numerical weather predictions, post-processed weather predictions and gridded observations): (1) AROME-Crocus, which used weather forecasts from the AROME-MetCoOp model, including post-processed air temperature and wind speed, and (2) GridObs-Crocus, which used gridded observations of temperature and precipitation combined with meteorological forecasts from AROME-MetCoOp. Snow simulations were carried out for 2 years (1 September 2014–31 August 2016). The main findings are as follows:
GridObs-Crocus provides the best estimates of the snow depth at individual stations with bias of 7 cm and RMSE of 28 cm. AROME-Crocus has a bias of 20 cm and RMSE of 56 cm.
AROME-Crocus provides the best representation of the spatial distribution of snow cover, particularly during the melting season. In GridObs-Crocus the spatial snow cover distribution is captured during winter, but underestimation of snow depth at high elevations (due to the low elevation bias in the gridded observation data set) likely causes the snow cover to decrease too soon during the melt season, leading to unrealistically little snow by the end of the season.
Forcing data consisting of post-processed NWP data (observations assimilated into the raw NWP weather predictions) are most promising for snow simulations, when larger regions are evaluated. Post-processed NWP data (AROME-Crocus) provide a more representative spatial representation for both high mountains and lowlands compared to interpolated observations (GridObs-Crocus).
In regions with steep terrain gradients, terrain adjustment of precipitation phase is highly important for improving the rainfall and snowfall determination when using NWP data.
Blowing snow (which is not simulated by SURFEX/Crocus) contributes to 17 % of all decreases in snow depth and to 50 % of the strongest reductions of more than 20 cm of snow depth loss in a day. Using the option in SURFEX/Crocus of running with sublimation in the case of snowdrift is not enough to address this issue.
To investigate the impact of using gridded observations of temperature and precipitation separately, “leave-one-out” experiments could be carried out (two extra experiments where one uses only gridded observations of temperature, and one uses only gridded observations of precipitation, while all other variables come from AROME-MetCoOp). Using the multi-physical ensemble system ESCROC (Ensemble System Crocus), and/or an ensemble of meteorological forcing would be another interesting topic for future work. Finally, when using AROME-MetCoOp as forcing data for running SURFEX/Crocus at a resolution higher than 2.5 km, terrain adjustment routines should be applied to the generation of forcing data. In this study we accounted for local terrain effects by using post-processed AROME-MetCoOp temperature and wind, but this could be extended to other variables.
The findings in this study have improved our understanding of regional snow modelling in Norway, which is important not only for water resource planning and flood forecasting but also for impact studies related to climate change and winter climate. Running the SURFEX/Crocus model in gridded version for Norwegian conditions using a combination of data sources (raw and post-processed weather predictions and observations) is very promising. The result from this study is very valuable information, which may be used for future development of a system for daily snow mapping in Norway.
Snow depth and meteorological variables from the stations used in this study are freely available through https://frost.met.no/ (Frost, 2018). AROME-MetCoOp forecasts are available through http://thredds.met.no/thredds/metno.html. The gridded data set of temperature and precipitation is available at http://thredds.met.no/thredds/catalog/metusers/senorge2/seNorge2/archive/catalog.html. Hourly temperature and precipitation data are available from 2010 up to the present day. For daily temperature and precipitation data, the archive goes back to 1957 and can be downloaded at http://doi.org/10.5281/zenodo.845733. The data are also shown on the web portals http://www.senorge.no and http://www.xgeo.no (both in Norwegian only). The SURFEX-Crocus simulations for both experiments can be made available for research purposes by contacting the authors.
The authors declare that they have no conflict of interest.
The authors are grateful for the funding of this study by the Research
Council of Norway through the project “Better SNOW models for prediction of
natural hazards and HydropOWer applications” (SNOWHOW), led by Thomas Skaugen
at the Norwegian Water Resources and Energy Directorate (NVE). We would like
to thank our fellow participants in the SNOWHOW project for invaluable help
and discussions: Thomas Skaugen, Tuomo Saloranta (NVE), Karsten Müller (NVE), Kjetil Melvold (NVE) and Sjur Kolberg (SINTEF). In addition, we would
like to thank Tuomo Saloranta for providing the processed MODIS images, and
Cristian Lussana (MET Norway) for valuable help and discussions. We are also
very grateful to two anonymous reviewers for their time, effort and very
helpful suggestions, which resulted in a greatly improved paper.
Edited by: Ross Brown
Reviewed by: two anonymous referees
Barfod, E., Müller, K., Saloranta, T., Andersen, J., Orthe, N., Wartianien, A., Humstad, T., Myrabø, S., and Engeset, R.: The expert tool XGEO and its applications in the Norwegian Avalanche Forecasting Service, in: International Snow Science Workshop Grenoble, 7–11 October 2013, Chamonix Mont-Blanc, France, 2013. a
Bartelt, P. and Lehning, M.: A physical SNOWPACK model for the Swiss avalanche warning: Part I: numerical model, Cold Reg. Sci. Technol., 35, 123–145, http://www.sciencedirect.com/science/article/pii/S0165232X02000745, 2002. a
Bellaire, S., Jamieson, J. B., and Fierz, C.: Corrigendum to “Forcing the snow-cover model SNOWPACK with forecasted weather data” published in The Cryosphere, 5, 1115–1125, 2011, The Cryosphere, 7, 511–513, https://doi.org/10.5194/tc-7-511-2013, 2013. a, b
Bergstrøm, S.: Development and application of a conceptual runoff model for Scandinavian catchments, SMHI report RH07, Swedish Meteorological and Hydrological Institute, Norrköping, Sweden, 1976. a
Bernier, N. B., Bélair, S., Bilodeau, B., and Tong, L.: Near-surface and land surface forecast system of the Vancouver 2010 Winter Olympic and Paralympic Games, J. Hydrometeorol., 12, 508–530, 2011. a
Best, M. J., Pryor, M., Clark, D. B., Rooney, G. G., Essery, R. L. H., Ménard, C. B., Edwards, J. M., Hendry, M. A., Porson, A., Gedney, N., Mercado, L. M., Sitch, S., Blyth, E., Boucher, O., Cox, P. M., Grimmond, C. S. B., and Harding, R. J.: The Joint UK Land Environment Simulator (JULES), model description – Part 1: Energy and water fluxes, Geosci. Model Dev., 4, 677–699, https://doi.org/10.5194/gmd-4-677-2011, 2011. a
Bokhorst, S., Pedersen, S. H., Brucker, L., Anisimov, O., Bjerke, J. W., Brown, R. D., Ehrich, D., Essery, R. L. H., Heilig, A., Ingvander, S., Johansson, C., Johansson, M., Jónsdóttir, I. S., Inga, N., Luojus, K., Macelloni, G., Mariash, H., McLennan, D., Rosqvist, G. N., Sato, A., Savela, H., Schneebeli, M., Sokolov, A., Sokratov, S. A., Terzago, S., Vikhamar-Schuler, D., Williamson, S., Qiu, Y., and Callaghan, T. V.: Changing Arctic snow cover: A review of recent developments and assessment of future needs for observations, modelling, and impacts, Ambio, 45, 516–537, https://doi.org/10.1007/s13280-016-0770-0, 2016. a
Boone, A. and Etchevers, P.: An inter-comparison of three snow schemes of varying complexity coupled to the same land-surface model: Local scale evaluation at an Alpine site, J. Hydrometeorol., 2, 374–394, 2001. a, b
Boone, A., Masson, V., Meyers, T., and Noilhan, J.: The Influence of the Inclusion of Soil Freezing on Simulations by a Soil–Vegetation–Atmosphere Transfer Scheme, J. Appl. Meteorol., 39, 1544–1569, https://doi.org/10.1175/1520-0450(2000)039<1544:TIOTIO>2.0.CO;2, 2000. a
Brown, R., Vikhamar-Schuler, D., Bulygina, O., Derksen, C., Luojus, K., Mudryk, L., Wang, L., and Yang, D.: Arctic terrestrial snow cover. Chapter 3, in: Snow, Water, Ice and Permafrost in the Arctic (SWIPA) 2017, Arctic Monitoring and Assessment Programme (AMAP), Oslo, Norway, 25–64, 2017. a, b
Brown, R. D. and Robinson, D. A.: Northern Hemisphere spring snow cover variability and change over 1922–2010 including an assessment of uncertainty, The Cryosphere, 5, 219–229, https://doi.org/10.5194/tc-5-219-2011, 2011. a
Brun, E., David, P., Sudul, M., and Brunot, G.: A numerical model to simulate snow-cover stratigraphy for operational avalanche forecasting, J. Glaciol., 38, 13–22, https://doi.org/10.3189/S0022143000009552, 1992. a
Brun, E., Vionnet, V., Boone, A., Decharme, B., Peings, Y., Valette, R., Karbou, F., and Morin, S.: Simulation of northern Eurasian local snow depth, mass and density using a detailed snowpack model and meteorological reanalysis, J. Hydrometeorol., 14, 203–219, https://doi.org/10.1175/JHM-D-12-012.1, 2013. a
Carrera, M. L., Bélair, S., Fortin, V., Bilodeau, B., Charpentier, D., and Doré, I.: Evaluation of snowpack simulations over the Canadian Rockies with an experimental hydrometeorological modeling system, J. Hydrometeorol., 11, 1123–1140, 2010. a
Carrera, M. L., Bélair, S., and Bilodeau, B.: The Canadian land data assimilation system (CaLDAS): Description and synthetic evaluation study, J. Hydrometeorol., 16, 1293–1314, 2015. a
Douville, H., Royer, J.-F., and Mahfouf, J.-F.: A new snow parameterization for the Météo-France climate model – Part2: Validation in a 3-D GCM experiment, Clim. Dynam., 12, 37–52, 1995. a
Dyrrdal, A. V., Saloranta, T., Skaugen, T., and Stranden, H. B.: Changes in snow depth in Norway during the period 1961-2010, Hydrol. Res., 44, 169–179, 2013. a
Engeset, R.: National Avalanche Warning Service for Norway, Established 2013, in: International Snow Science Workshop Grenoble, 7–11 October 2013, Chamonix Mont-Blanc, France, 2013. a
Essery, R., Morin, S., Lejeune, Y., and Ménard, C. B.: A comparison of 1701 snow models using observations from an alpine site, Adv. Water Res., 55, 131–148, https://doi.org/10.1016/j.advwatres.2012.07.013, 2013. a, b
Etchevers, P., Martin, E., Brown, R., Fierz, C., Lejeune, Y., Bazile, E., Boone, A., Dai, Y., Essery, R., Fernandez, A., Gusev, Y., Jordan, R., Koren, V., Kowalcyzk, E., Nasonova, N., Pyles, R., Schlosser, A., Shmakin, A., Smirnova, T., Strasser, U., Verseghy, D., Yamazaki, T., and Yang, Z.: Validation of the energy budget of an alpine snowpack simulated by several snow models (SnowMIP project), International Symposium on Snow and Avalanches, Davos, Switzerland, 2–6 June 2003, Ann. Glaciol., 38, 150–158, https://doi.org/10.3189/172756404781814825, 2004. a
FAO/IIASA/ISRIC/ISS-CAS/JRC: Harmonized World Soil Database (version 1.2), Tech. rep., FAO, Rome, Italy and IIASA, Laxenburg, Austria, 2012. a
Fierz, C., Bavay, M., Wever, N., and Lehning, M.: SNOWPACK: where do we stand today?, in: International Snow Science Workshop, EPFL-TALK-197625, 2013. a
Habets, F., Boone, A., and Noilhan, J.: Simulation of a Scandinavian basin using the diffusion transfer version of ISBA, Global Planet. Change, 38, 137–149, https://doi.org/10.1016/S0921-8181(03)00016-X, 2003. a
Hanssen-Bauer, I., Førland, E. J., Haddeland, I., Hisdal, H., Mayer, S., Nesje, A., Nilsen, J., Sandven, S., Sandø, A., Sorteberg, A., and Ådlandsvik, B.: Klima i Norge 2100, Kunnskapsgrunnlag for klimatilpasning oppdatert i 2015, Tech. Rep. 2, Norsk klimaservicesenter, available at: https://cms.met.no/site/2/klimaservicesenteret/klima-i-norge-2100/_attachment/10990 (last access: 18 June 2018), 2015. a, b, c, d
Hanssen-Bauer, I., Førland, E. J., Haddeland, I., Hisdal, H., Lawrence, D., Mayer, S., Nesje, A., Nilsen, J., Sandven, S., Sandø, A., Sorteberg, A., and Ådlandsvik, B.: Climate in Norway 2100 – A knowledge base for climate adaptation, Tech. Rep. 1, Norwegian Climate Service Centre, 2017. a
Homleid, M.: Diurnal corrections of short-term surface temperature forecasts using the Kalman filter, Weather Forecast., 4, 689–707, 1995. a
Homleid, M. and Tveter, F. T.: Verification of Operational Weather Prediction Models December 2015 to February 2016, Met.no report 18, Norwegian Meteorological Institute, Oslo, Norway, available at: https://www.met.no/publikasjoner/met-info, 2016. a, b
Johansson, C., Pohjola, V., Jonasson, C., and Callaghan, T.: Multi-decadal changes in snow characteristics in sub-Arctic Sweden, Ambio, 40, 566–574, 2011. a
Køltzow, M.: MetCoOp Ensemble Prediction System (MEPS), Norwegian Meteorological Institute, Oslo, Norway, available at: https://drive.google.com/file/d/0B-SaEtrDE91WWEJoNkJiUm5TNzg/view (last access: 18 June 2018), 2017. a, b, c
Köppen, W.: Das geographische System der Klimate, Handbuch der Klimatologie, vol. 1, Verlag von Gebrüder Borntraeger, Berlin, Germany, 1936. a
Lafaysse, M., Cluzet, B., Dumont, M., Lejeune, Y., Vionnet, V., and Morin, S.: A multiphysical ensemble system of numerical snow modelling, The Cryosphere, 11, 1173–1198, https://doi.org/10.5194/tc-11-1173-2017, 2017. a, b, c
Lehning, M., Bartelt, P. B., Brown, R. L., Fierz, C., and Satyawali, P.: A physical SNOWPACK model for the Swiss Avalanche Warning Services – Part II: Snow Microstructure, Cold Reg. Sci. Technol., 35, 147–167, 2002. a
Li, L. and Pomeroy, J. W.: Estimates of Threshold Wind Speeds for Snow Transport Using Meteorological Data, J. Appl. Meteorol., 36, 205–213, https://doi.org/10.1175/1520-0450(1997)036<0205:EOTWSF>2.0.CO;2, 1997. a
Lussana, C., Saloranta, T., Skaugen, T., Magnusson, J., Tveito, O. E., and Andersen, J.: seNorge2 daily precipitation, an observational gridded dataset over Norway from 1957 to the present day, Earth Syst. Sci. Data, 10, 235–249, https://doi.org/10.5194/essd-10-235-2018, 2018a. a, b, c, d, e, f, g, h, i, j, k
Magnusson, J., Wever, N., Essery, R., Helbig, N., Winstral, A., and Jonas, T.: Evaluating snow models with varying process representations for hydrological applications, Water Resour. Res., 51, 2707–2723, https://doi.org/10.1002/2014WR016498, 2015. a
Masson, V., Le Moigne, P., Martin, E., Faroux, S., Alias, A., Alkama, R., Belamari, S., Barbu, A., Boone, A., Bouyssel, F., Brousseau, P., Brun, E., Calvet, J.-C., Carrer, D., Decharme, B., Delire, C., Donier, S., Essaouini, K., Gibelin, A.-L., Giordani, H., Habets, F., Jidane, M., Kerdraon, G., Kourzeneva, E., Lafaysse, M., Lafont, S., Lebeaupin Brossier, C., Lemonsu, A., Mahfouf, J.-F., Marguinaud, P., Mokhtari, M., Morin, S., Pigeon, G., Salgado, R., Seity, Y., Taillefer, F., Tanguy, G., Tulet, P., Vincendon, B., Vionnet, V., and Voldoire, A.: The SURFEXv7.2 land and ocean surface platform for coupled or offline simulation of earth surface variables and fluxes, Geosci. Model Dev., 6, 929–960, https://doi.org/10.5194/gmd-6-929-2013, 2013. a, b, c, d
Mohr, M.: New routines for gridding of temperature and precipitation observations for seNorge.no, Met. no Report 8, Norwegian Meteorological Institute, Oslo, Norway, 2008. a
Müller, M., Homleid, M., Ivarsson, K.-I., Køltzow, M. A. Ø., Lindskog, M., Midtbø, K. H., Andrae, U., Aspelien, T., Berggren, L., Bjørge, D., Dahlgren, P., Kristiansen, J., Randriamampianina, R., Ridal, M., and Vignes, O.: AROME-MetCoOp: A Nordic Convective-Scale Operational Weather Prediction Model, Weather Forecast., 32, 609–627, https://doi.org/10.1175/WAF-D-16-0099.1, 2017. a, b, c, d, e, f, g
Quéno, L., Vionnet, V., Dombrowski-Etchevers, I., Lafaysse, M., Dumont, M., and Karbou, F.: Snowpack modelling in the Pyrenees driven by kilometric-resolution meteorological forecasts, The Cryosphere, 10, 1571–1589, https://doi.org/10.5194/tc-10-1571-2016, 2016. a, b, c, d, e, f, g
Quéno, L., Karbou, F., Vionnet, V., and Dombrowski-Etchevers, I.: Satellite products of incoming solar and longwave radiations used for snowpack modelling in mountainous terrain, Hydrol. Earth Syst. Sci. Discuss., https://doi.org/10.5194/hess-2017-563, in review, 2017. a, b
Raleigh, M. S., Lundquist, J. D., and Clark, M. P.: Exploring the impact of forcing error characteristics on physically based snow simulations within a global sensitivity analysis framework, Hydrol. Earth Syst. Sci., 19, 3153–3179, https://doi.org/10.5194/hess-19-3153-2015, 2015. a
Rasmus, S., Boelhouwers, J., Briede, A., Brown, I., Falarz, M., Ingvander, S., Jaagus, J., Kitaev, L., Mercer, A., and Rimkus, E.: Recent change – Terrestrial cryosphere, in: Second Assessment of Climate Change for the Baltic Sea Basin, edited by: The BACC II Author Team, 117–129, https://doi.org/10.1007/978-3-319-16006-1_6, Springer, Cham, 2015. a
Ruan, G. and Langsholt, E.: Rekalibrering av flomvarslingas HBV-modeller med inndata fra seNorge, versjon 2.0, Tech. Rep. 71, NVE Report, Oslo, Norway, 2017. a
Sælthun, N. R.: The Nordic HBV model, NVE Report No. 7, Norwegian Water Resources and Energy Administration, Oslo, Norway, 1996. a
Salomonson, V. V. and Appel, I.: Estimating fractional snow cover from MODIS using the normalized difference snow index, Remote Sens. Environ, 89, 351–360, https://doi.org/10.1016/j.rse.2003.10.016, 2004. a
Sauter, T. and Obleitner, F.: Assessing the uncertainty of glacier mass-balance simulations in the European Arctic based on variance decomposition, Geosci. Model Dev., 8, 3911–3928, https://doi.org/10.5194/gmd-8-3911-2015, 2015. a
Seity, Y., Brousseau, P., Malardel, S., Hello, G., Bénard, P., Bouttier, F., Lac, C., and Masson, V.: The AROME-France convective-scale operational model, Mon. Weather Rev., 139, 976–991, https://doi.org/10.1175/2010MWR3425.1, 2011. a
Skaugen, T.: Studie av skilletemperatur for snø ved hjelp av samlokalisert snøpute, nedbør- og temperaturdata, Norges vassdrags- og energiverk, Oslo, Norway, 1998. a
Skaugen, T., Stranden, H. B., and Saloranta, T.: Trends in snow water equivalent in Norway (1931–2009), Hydrol. Res., 43, 489–499, 2012. a
Skaugen, T., Luijting, H., Saloranta, T., Vikhamar-Schuler, D., and Müller, K.: In search of operational snow model structures for the future – comparing four snow models for 17 catchments in Norway, Hydrol. Res., https://doi.org/10.2166/nh.2018.198, online first, 2018. a, b
Solberg, R., Amlien, J., and Koren, H.: A review of optical snow cover algorithms. Norwegian Computing, Tech. Rep. SAMBA/40/06, Norwegian Computing Center, Oslo, Norway, 2006. a
Sturm, M., Holmgren, J., and Liston, G. E.: A seasonal snow cover classification system for local to global applications, J. Climate, 8, 1261–1283, 1995. a
Vernay, M., Lafaysse, M., Mérindol, L., Giraud, G., and Morin, S.: Ensemble forecasting of snowpack conditions and avalanche hazard, Cold Reg. Sci. Tech., 120, 251–262, https://doi.org/10.1016/j.coldregions.2015.04.010, 2015. a, b
Vikhamar-Schuler, D., Hanssen-Bauer, I., Schuler, T. V., Mathiesen, S. D., and Lehning, M.: Use of a multilayer snow model to assess grazing conditions for reindeer, Ann. Glaciol., 54, 214–226, https://doi.org/10.3189/2013AoG62A306, 2013. a
Vikhamar-Schuler, D., Isaksen, K., Haugen, J. E., Tømmervik, H., Luks, B., Schuler, T. V., and Bjerke, J. W.: Changes in winter warming events in the Nordic Arctic Region, J. Climate, 29, 6223–6244, https://doi.org/10.1175/JCLI-D-15-0763.1, 2016. a
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, https://doi.org/10.5194/gmd-5-773-2012, 2012. a, b, c
Vionnet, V., Martin, E., Masson, V., Guyomarc'h, G., Naaim-Bouvet, F., Prokop, A., Durand, Y., and Lac, C.: Simulation of wind-induced snow transport and sublimation in alpine terrain using a fully coupled snowpack/atmosphere model, The Cryosphere, 8, 395–415, https://doi.org/10.5194/tc-8-395-2014, 2014. a
Vionnet, V., Dombrowski-Etchevers, I., Lafaysse, M., Quéno, L., Seity, Y., and Bazile, E.: Numerical Weather Forecasts at Kilometer Scale in the French Alps: Evaluation and Application for Snowpack Modeling, J. Hydrometeorol., 17, 2591–2614, 2016. a, b, c, d, e, f