Translator Disclaimer
8 December 2017 Estimation of subpixel snow sublimation from multispectral satellite observations
Author Affiliations +
Snow sublimation is an important hydrological process and its spatial and temporal variation remains largely unknown; however, few studies have been conducted to quantify its spatial variability. Our study focuses on the evaluation of two algorithms, Penman–Monteith (P–M) equation and the bulk aerodynamic (BA) parameterization of snow sublimation. The two methods were first evaluated against eddy covariance (EC) measurements of latent heat flux at towers located in the upper reaches of the Heihe River Basin (China). Both methods were in good agreement with the ground observations with high coefficient of determination (<inline-formula> <math display="inline" xmlns:mml=""> <mrow> <msup> <mi>R</mi> <mn>2</mn> </msup> </mrow> </math> </inline-formula>) and low root mean squared error (RMSE). Next, we estimated subpixel snow sublimation using remote sensing data at a <inline-formula> <math display="inline" xmlns:mml=""> <mrow> <mn>1</mn> <mtext>-</mtext> <mi>km</mi> <mo>×</mo> <mn>1</mn> <mtext>-</mtext> <mi>km</mi> </mrow> </math> </inline-formula> spatial resolution. The results based on satellite data were evaluated against ground measurements at the two experimental sites. The P–M equation gave <inline-formula> <math display="inline" xmlns:mml=""> <mrow> <msup> <mi>R</mi> <mn>2</mn> </msup> <mo>=</mo> <mn>0.75</mn> </mrow> </math> </inline-formula>, <inline-formula> <math display="inline" xmlns:mml=""> <mrow> <mi>RMSE</mi> <mo>=</mo> <mn>8.4</mn> <mtext>  </mtext> <mi mathvariant="normal">W</mi> <mtext> </mtext> <msup> <mrow> <mi mathvariant="normal">m</mi> </mrow> <mrow> <mo>−</mo> <mn>2</mn> </mrow> </msup> </mrow> </math> </inline-formula> for Dashalong site and <inline-formula> <math display="inline" xmlns:mml=""> <mrow> <msup> <mi>R</mi> <mn>2</mn> </msup> <mo>=</mo> <mn>0.36</mn> </mrow> </math> </inline-formula>, <inline-formula> <math display="inline" xmlns:mml=""> <mrow> <mi>RMSE</mi> <mo>=</mo> <mn>9.1</mn> <mtext>  </mtext> <mi mathvariant="normal">W</mi> <mtext> </mtext> <msup> <mrow> <mi mathvariant="normal">m</mi> </mrow> <mrow> <mo>−</mo> <mn>2</mn> </mrow> </msup> </mrow> </math> </inline-formula> for the Dadongshu site and performed better than the BA parameterization, which gave <inline-formula> <math display="inline" xmlns:mml=""> <mrow> <msup> <mi>R</mi> <mn>2</mn> </msup> <mo>=</mo> <mn>0.65</mn> </mrow> </math> </inline-formula>, <inline-formula> <math display="inline" xmlns:mml=""> <mrow> <mi>RMSE</mi> <mo>=</mo> <mn>17.5</mn> <mtext>  </mtext> <mi mathvariant="normal">W</mi>



Snow plays an important role in the climate system of high elevations and cold regions because of its high albedo and low heat conductivity and it is also an important water resource for these regions.14 Snow sublimation has been identified as an important hydrological process and results in a significant loss of snowpack water, which involves complex mass and energy exchanges.47 Extreme cases of sublimation have been shown to be very efficient at removing snowpack water with losses of up to 90% of annual snowpack on preferred alpine crests.5 A wide range of snow sublimation rates have been reported in the literature.711 High sublimation rates were observed in the Colorado Front Range with daily values of 2.35 mm and in the alpine region of the Sierra Nevada with daily values of 2.17 mm.8,9 However, sublimation rates were reported with 0.39 and 0.15  mmd1 at wind-exposed and wind-sheltered sites, respectively, at Owyhee Mountains of the USA.7 These various results may be explained by environmental differences among study sites because snow sublimation is affected by multiple factors such as wind speed, net radiation, air temperature, and vapor pressure deficit,5 and the spatial variability of snow sublimation is dependent on the spatial variability of these factors.5,10,11

Sublimation as the moisture fluxes between the snow and the atmosphere is commonly measured by an evaporation pan. The evaporation pan has been widely used with the advantage of simplicity and economy but with more manual work.4,1215 Schmidt et al.13 estimated that sublimation from snowpack was 78-mm water equivalent and represented about 20% of the normal peak water equivalent of the snowpack by a 65-cm diameter pan during a 40-day accumulation period in a Colorado subalpine forest. However, measurements by an evaporation pan were found to be invalid when precipitation or blowing snow occurred.15 The eddy covariance (EC) method is another approach to measure snow sublimation and is widely used.7,1621 EC uses high-frequency (10 to 20 Hz) measurements of wind vector and scalar quantities as the basis for flux computations. A detailed description of the EC method and data processing to determine sublimation can be found in Reba et al.17 EC sensors respond poorly during precipitation events,7,17 however. Additionally, the so-called energy balance closure problem is a widespread occurrence in EC measurements with reported discrepancy of 10% to 30%.22,23 Despite these challenges, EC data can be a very valuable data source for snow research when thoroughly corrected to ensure accuracy.7,17,18,21

In addition to in situ observation, methods to estimate snow sublimation were developed in the past several decades. Initially, empirical formulas were widely adopted, based on the statistical relationship between snow sublimation and meteorological conditions (e.g., wind speed and vapor pressure deficit).24,25 Strasser et al.5 used empirical formulas to estimate sublimation from the snowpack to be about 100-mm snow water equivalent (SWE) during the 2003/2004 winter period in a mountainous region in the southeast of Germany. Empirical methods are simple but have low accuracy and poor physical basis.8,10 The aerodynamic profile method was developed later based on complex turbulent transfer theory and is more reliable than the empirical methods.8,10,14 Hood et al.8 used the aerodynamic profile method to estimate snow sublimation: total net sublimation was 195-mm water equivalent or 15% of annual snow accumulation at Niwot Ridge in the Colorado Front Range during the 1994/1995 snow season. The snowpack is generally assumed to be a saturated surface, which has been proved by many studies,2630 thus the bulk aerodynamic (BA) method was developed to simplify the aerodynamic profile method and remove the need for multilevel observations. The BA model needs only meteorological measurements at a single height above the snowpack and soon became the most widely used method due to a widespread lack of multilevel meteorological data.7,2630 Reba et al.7 found that the measured snow sublimation by EC was very similar to estimated snow sublimation by the BA model with a correlation coefficient of 0.85 and root mean square difference of 0.15  mmd1.

Many of the physical snow sublimation models have been based upon the energy balance approach.21,29,31,32 The relation between the amount of energy used for snow sublimation with net radiation and subsurface heat flux can be described by a conceptual model of the energy balance of a control volume, which applies to a single, thin layer of snow placed at the air–snow interface (Fig. 1).

Fig. 1

Conceptual model of snow surface energy balance.


The snow surface energy balance equation can be written as

Eq. (1)

where the net radiation flux (Rn) is the principal driver of the snow sublimation (Wm2), the sensible heat flux (H) is the energy transfer due to convection in the surface boundary layer (Wm2), the snow subsurface heat flux (Gs) is the energy transfer to the snowpack (Wm2), and the latent heat flux (LE) is the energy used in the phase transition of frozen/liquid water to vapor (Wm2).

The method based on the snow surface energy balance also applies one-level meteorological measurements to estimate snow sublimation with the P–M equation, which has been widely adopted in estimating evaporation.33,34 The P–M method for estimating latent heat flux calculation has been extensively used in hydrological, atmospheric, and environmental modeling, and can be used to estimate snow sublimation as air in the vicinity of the snowpack that is generally assumed to be vapor saturated.29,31,32

Furthermore, many studies on snow sublimation are constrained by the unavailable or very sparse ground observation. At high latitudes, meteorological stations are sparse and located mostly at lower elevations.2 However, snow is more abundant at high latitudes and altitudes. Although many ground-based observations and modeling studies of snow sublimation have provided valuable information (i.e., snow surface temperature and snow depth) at a particular location, it is difficult to extend such information to large scales. The spatial variability of snow sublimation over the vast regions remains largely unknown.2 Satellite observations have long been recognized as effective in providing spatially distributed snow information to estimate snow sublimation.35 Remotely sensed data provide temporally and spatially continuous information over snow surface, i.e., snow surface temperature and fractional snow cover (FSC).36,37 The moderate resolution imaging spectroradiometers (MODIS) onboard NASA’s Terra and Aqua satellites provide unprecedented information regarding snow and surface energies, which can be used for regional- and global-scale snow sublimation estimation in near real-time.38,39 However, to our knowledge, very few studies have been published on the estimation of snow sublimation with satellite observations.

The BA method and P–M equation were employed in this study because both methods assume that the snow surface is vapor saturated and use atmospheric measurements at only one level above the snow surface. This is the foundation of calculating the snow sublimation by remote sensing because we cannot get gradient measurements in the surface boundary layer with satellite observations. We present an approach to estimate subpixel snow sublimation using multiple satellite observations by characterizing the subpixel energy balance of the snow cover fraction. We have evaluated our method using both tower and satellite observations as input and compared alternate parameterizations of snow sublimation. The objectives of our study were: (1) to evaluate the P–M equation and the BA parameterization using tower observations, (2) to analyze the energy balance of the subpixel snow cover fractional area, and (3) to present and validate the regional snow sublimation estimated using satellite data.


Data and Materials


Description of Study Area and Measurements

The upper reaches of the Heihe River Basin (Fig. 2) are located between 37°41′ N to 39°05′ N and 98°34′ E to 101°11′ E, and cover an area of 2.85×104  km2. Elevation ranges from 1631 to 5245 m, the mean annual precipitation is 350 mm and the mean annual air temperature is 2.0°C.40 Glaciers covering 421  km2 are distributed above elevation 4500 m. At elevations above 4000 m, vegetation is very sparse and dominated by cushion plants, while meadows and shrubs occur below 3300 m.41 Temporary snow usually exists under 2700 m, patchy snowpack from 2700 to 3400 m, and permanent snowpack exists above 3400 m.42 Snowmelt provides most water resources in the Heihe River Basin, which supplies agriculture in the middle reaches and the arid ecosystem in the lower reaches.43

Fig. 2

The upper reaches of the Heihe River Basin and location of the research sites.


The data were collected at two observation stations, Dadongshu mountain pass station (4101 m, 100.23°E, 38.01°N) and Dashalong station (3739 m, 98.94°E, 38.84°N), part of the Heihe Watershed Allied Telemetry Experimental Research (HiWATER) project.44,45 These stations were selected because of the locally flat terrain and homogeneous land cover without major obstacles, such as houses and trees. The land cover at these sites is alpine meadow and bare soil. The land surface at these sites is always covered by snow most of the time in the late autumn, winter, and early spring, but snowfall is larger in spring and autumn than in winter.46

The sensible and latent heat flux data during the period between November 1, 2014, and January 31, 2015, at Dadongshu station, and October 17, 2014, and December 31, 2014, at Dashalong station were collected. The fluxes were measured by the EC systems installed at 3 m at the Dadongshu site and 4.5 m at Dashalong site above the ground. The raw data were sampled at 10 Hz, then 30-min mean values were obtained after a series of corrections, such as data filtering, sonic temperature correction, air density correction, and coordinate rotation.44 The auxiliary data, including the net radiation, air temperature, relative humidity, snow depth, air pressure, wind speed, and snow infrared surface temperature, were also observed by automatic weather stations (AWS) near the EC system in the same period. All the data were processed and recorded at 30-min intervals. There were no snow depth measurements at Dashalong station. The latent heat flux at this station was used to validate the estimates based on satellite data.


Satellite and Regional Forcing Data

The satellite and regional forcing data used in this study are listed in Table 1. They are acquired from various sources and comprised of well-validated products.

Table 1

Remote sensing data and regional atmospheric forcing data used for computing snow sublimation at regional scale in the upstream of Heihe River Basin from November 1, 2014, to January 31, 2015.

VariablesSourceTemporal resolutionSpatial resolution
Downward longwave radiationHiWATERHourly5 km
Downward shortwave radiationHiWATERHourly5 km
Air pressureHiWATERHourly5 km
Air temperatureHiWATERHourly5 km
Wind speedHiWATERHourly5 km
Specific humidityHiWATERHourly5 km
PrecipitationHiWATERHourly5 km
Surface temperatureNASA, MOD11A1Daily1 km
FSCNSIDC, MOD10A1Daily500 m
Land cover typeHiWATER1 km

Land surface temperature was retrieved using the radiometric data acquired by the MODIS on Terra and downloaded from National Aeronautics and Space Administration (NASA), which covers the entire upper reach of the Heihe River Basin.

We used the FSC retrieved with the MODIS radiometric data (MOD10A1 product) and downloaded from the National Snow and Ice Data Center (NSIDC). The MOD10A1 FSC is estimated by applying an empirical relationship with the normalized difference snow index (NDSI), established using FSC retrieved at higher spatial resolution with Landsat ETM+ data.37 The original MOD10A1 has a 500-m spatial resolution and was resampled to a 1-km grid.

Land cover type is required to unmix the land surface temperature and estimate the radiometric surface temperature of the snow cover fraction. The land cover map at a spatial resolution of 1 km was obtained from the Cold and Arid Regions Sciences Data Center.47,48 The land cover map was generated by combining multisource land cover/land use classification maps including a 1:1,000,000 scale vegetation map, a 1:100,000 scale land use map for the year 2000, a 1:1,000,000 scale swamp-wetland map, a glacier map, and a MODIS land cover map for China in 2001.48 The dominant land cover types in winter in these regions were forest and bare soil. We aggregated the land cover classes to just two: (a) snow and forest and (b) snow and bare soil.

The regional atmospheric forcing data were obtained from the Cold and Arid Regions Sciences Data Center.4951 Data with a spatial resolution of 5 km were downscaled to 1 km by applying a statistical downscaling approach to hourly downward longwave and shortwave radiations, air temperature, air pressure, and specific humidity.52 Data of wind speed and precipitation were downscaled to 1 km using the bilinear resampling method.



Two methods to estimate snow sublimation were evaluated: the widely used BA method and P–M combination equation adapted to snow cover. The two methods were applied assuming the air near the surface to be vapor saturated and the results were compared with the snow sublimation measured by the EC system in two experimental stations in the upper reaches of the Heihe River Basin in China. The theoretical background of these two methods is described in the following sections.


Penman–Monteith Combination Equation

The Penman53 equation was first derived to estimate evaporation from open water and water-saturated surfaces, where equilibrium at the surface, i.e., under equilibrium with vapor–saturated air. Later, this equation was extended by Monteith33 to nonsaturated surface, i.e., to actual evaporation. The combination equation is obtained by rewriting the surface energy balance equation using the Clausius–Clapeyron equation for (liquid and frozen) water–vapor equilibrium. In Penman and Monteith a liquid water–vapor equilibrium was assumed, while in our case the equilibrium is frozen water (ice)–vapor, i.e., we need to use the appropriate dependence of the saturated water vapor pressure on temperature. The P–M combination equation was applied in the form:33

For full snow cover

Eq. (2)



Eq. (3)

where LE (Wm2) is the latent heat flux with L as latent heat of sublimation (Jkg1) and E (kgs1m2) as sublimation, FSC (–) is the FSC, Δ (PaK1) is the slope of the curve relating saturated vapor pressure in the ice–water vapor equilibrium esat (Pa) to air temperature (Ta), and for ice or snow, Δ=21.87×265.5×esat(265.5+Ta)2, (here Ta is in °C).31 Rn (Wm2) is the net radiation flux, Gs (Wm2) is the snow subsurface heat flux due to the energy transfer in the snowpack, and Gs was estimated using a linear relationship between Gs and Rn obtained from AWS and EC observations. ρ (kgm3) is the air density, Cp (Jkg1K1) is the specific heat capacity of air, esat (Pa) is the saturation vapor pressure at the air temperature (Ta) in the frozen water–air equilibrium, i.e., esat=611×exp(21.87×TaTa+265.5) (here Ta is in °C). e (Pa) is the actual water vapor pressure, ra (sm1) is the aerodynamic resistance, and γ (PaK1) is the psychrometric constant. Equation (2) is a particular form of the combination equation. Menenti derived a general equation and showed that both the Monteith and Penman equations can be obtained as limiting conditions of the Menenti equation. 54 Both the Monteith and the Penman equations applied to a case, where the liquid (frozen) water–vapor phase transition occurs at a surface, where air is vapor saturated at its temperature. More in general, the liquid (frozen) water–vapor phase transition occurs at locations underneath the surface of the evaporating (sublimating) body, such as soil, stomata, or the snowpack particularly when meltwater may be present within the snowpack. Our assumptions, i.e., that the sublimation occurs at the surface of the snowpack and that air in the vicinity of the surface is vapor saturated at the temperature of the snowpack surface, may lead to differences between sublimation estimated with Eq. (2) and total latent heat flux measured by an eddy covariance device. The latter may include contribution due to subsurface sublimation and evaporation of meltwater.

The aerodynamic resistance is calculated as follows:

Eq. (4)

where z (m) is the height of measurement, z0 (m) is the snow surface roughness length, estimated to be z0=0.0002  m according to Dewalle and Rango;30 k is the von Karman’s constant (k=0.4), Uz (ms1) is the wind speed at height z (m), ϕM (−) is the flux–profile relationship for a stable or unstable boundary layer and can be expressed as follows:28,30

Eq. (5)

Stable condition:ϕM=(15  Ri)2,

Eq. (6)

Unstable condition:ϕM=(116  Ri)0.75,
where Ri (−) is the Richardson number given as28

Eq. (7)

where g is the acceleration due to gravity (9.8  ms2), Ta (K) is the air temperature, Tsnow (K) is snow surface temperature, and Tmean=(Ta+Tsnow)/2 in K.

The MOD11A1 land surface radiometric temperature retrieved from MODIS radiometric data is an average applying to an entire pixel. We assumed that in the case of mixed pixels, e.g., composed of snow and other land cover types such as bare soil or vegetation, the fractional cover of snow FSC is the MODIS data product, while the fractional abundance (1FSC) is assigned to the land cover type given by a land cover map for the entire pixel. To retrieve the surface temperature of the snow fractional area, we applied the method proposed by Sun et al.55 The component temperatures of a mixed pixel are assumed to include just two end-members, in our case either snow and vegetation (forest) or snow and bare soil. The end-member radiometric temperatures were estimated using MODIS land surface temperature and the spectral emissivities in the MODIS bands 31 and 32

Eq. (8)


Eq. (9)

where A=ϵ31ϵm32, B=ϵ32ϵm31, C=fmϵm31ϵm32, D=fnϵn31ϵm32, E=fnϵn32ϵm31, F=fnϵn31, and G=fmϵm31; Tm and Tn are the end-member temperatures within a mixed pixel, Ts (K) is the MOD11A1 pixel land surface temperature, fm and fn are the fractional covers of end-members, ϵm31 and ϵn31 are the end-member emissivities in band 31 (0.982319 for snow, 0.9851 for forest, and 0.9832 for soil), ϵm32 and ϵn32 are the end-member emissivities in band 32 (0.962614 for snow, 0.9844 for forest, and 0.9731 for soil), ϵ31 and ϵ32 are the emissivities of the pixel in bands 31 and 32 and were estimated as fractional area weighted average of the end-member emissivities. If FSC=1, then the temperature of the pixel is the snow temperature; if 0<FSC<1, the fractional abundance (1FSC) is assigned to the land cover type given by a land cover map for the entire pixel. Once the fractional abundances fi have been determined, Eqs. (8) and (9) are applied to calculate subpixel snow temperature.

Net radiation flux (Rn) is the principal driver of the snow sublimation and was expressed as

Eq. (10)

where Rswd is downward shortwave radiation (Wm2), ϵ is the spectrally integrated emissivity of the snowpack (–), Rlwd is downward longwave radiation flux (Wm2), and α is the albedo of the snowpack (–) calculated by an empirical formula,56,57

Eq. (11)

For nonmelt condition:αt=αt1τa·t86,400,

Eq. (12)

For melt condition:αt=(αt1αmin)·e0.24×Δt86,400+αmin,
where Δt is the time step after snowfall (s), τa is the empirical coefficient (0.008), and αmin is the minimum snow albedo (0.5). Extensive analyses of snow albedo and snowfall events indicate that 3-mm SWE of snowfall is a threshold for newly fallen snow to cover the old snow and resets the snow surface albedo to a high value.58 In this work, a snowfall exceeding 3-mm SWE is defined as a “new snow” event, and the time between two consecutive “new snow” events is counted as length of elapsed time after the first snowfall. If the precipitation (snowfall here) from the forcing data is >3  mm (SWE), then the snow albedo will be the maximum value of 0.85.


Bulk Aerodynamic Model

According to the BA theory, the latent heat flux can be computed as:28,30

For full snow cover

Eq. (13)



Eq. (14)

where Pa   is the atmospheric pressure (Pa), es (Pa) is the saturation vapor pressure at the snow surface temperature (Tsnow), es=611×exp(21.87×TsnowTsnow+265.5), (here Tsnow is in °C). Ce (ms1) is the bulk transfer coefficient for vapor exchange, which can be expressed as follows:

Eq. (15)



Metrics for Accuracy Assessment

We applied multiple metrics to evaluate the accuracy of our estimated sublimation. The mean relative error (MRE) describes the ratio between the deviation of the estimate from the observation and the observed value. The root mean square error (RMSE) measures the absolute difference between estimation and observation. The coefficient of determination (R2) is a measure of the consistency of estimated with measured values, e.g., in time. These metrics were calculated as follows:

Eq. (16)


Eq. (17)


Eq. (18)

where Pi and Oi are the estimates and observations, respectively, with P and O being the averages and n is the number of estimate–observation pairs.


Results and Discussion


Meteorological Conditions

At Dadongshu station, there were snowfall events from October 10, 2014, as shown in Fig. 3 and the total snowfall was 41.0 mm (SWE). The maximum value was 6.8  mmd1 in October 10, 2014. After several snowfall events in October 2014, the snow depth on 1 November reached a maximum depth of 0.2 m, and then declined until January 31, 2015. On January 31, 2015, soil surface was exposed and LE measurements by EC included soil evaporation. So we choose the period from November 1, 2014, to January 31, 2015, because there was enough snow on the ground, particularly during this period, we could neglect the contribution of soil evaporation to the LE measured by our EC device.

Fig. 3

Time series of daily snowfall at Dadongshu station between October 10, 2014, and January 31, 2015.


Meteorological conditions from November 1, 2014, to January 31, 2015, are shown in Fig. 4. The snow depth near the instrument station was 0.2 m and decreased to 0.12 m during the study period. The mean air temperature was 13.3°C, ranging from highs of 4.4°C to lows near 24.3°C. The mean wind speed for the study period was 4.6  ms1 with a maximum of 13.8  ms1 and a minimum of nearly 0  ms1 at some times. The mean relative humidity at the Dadongshu site was 45.2%, ranging from 9.1% to 87.2%. High relative humidity, which means large amount of water vapor in the air, will inhibit snow sublimation.

Fig. 4

Time series of (a) snow depth, (b) air temperature, (c) wind speed, (d) relative humidity with 30-min interval at Dadongshu station between November 1, 2014, and January 31, 2015.



Snow Subsurface Heat Flux

The snow subsurface heat flux (Gs) was not measured and we estimated as the residual of the energy balance at the surface of the snowpack. This was done only in the Dadongshu site, when it was covered by a homogeneous snowpack during the snow season, for all available complete daytime (11:00 to 17:00) measurements between November 1, 2014, and January 31, 2015. Next, we determined a linear correlation between Gs and Rn, obtaining a regression coefficient of 0.575 and R2=0.63 (Fig. 5)

Eq. (19)


Fig. 5

Linear relationship between snow subsurface heat flux (Gs) and net radiation (Rn) based on the in situ measurements at the Dadongshu site (30–min. data between 11:00 and 17:00 from November 1, 2014, to January 31, 2015, with fully covered snow).


This relationship in Eq. (19) is further applied to estimate regional snow sublimation using satellite data.


Validation of Estimated Snow Sublimation

We evaluated first estimates of sublimation obtained with both Eqs. (2) and (13) using in situ measurements at the Dadongshu site of all the energy balance components, i.e., net radiation, latent heat flux, and sensible heat flux for the period with complete snow cover (FSC=1). The R2 obtained with the P–M equation [Eq. (2)] and the BA parameterization [BA, Eq. (13)] were 0.65 and 0.54, respectively (Fig. 6), while RMSE was 10.48 and 10.21  Wm2, respectively. The MRE for P–M was 28.1%, while it was 1.12% for BA. The relatively small MRE of BA is due to the offset effect of the positive and negative biases. It is clear that the P–M equation overestimates the sublimation of snow more than the BA approach (Fig. 7). Overall, the evaluation against in situ measurements showed that the performance of both P–M and BA is satisfactory, especially taking into account the low RMSE values.

Fig. 6

(a) P–M and (b) BA estimates of snow sublimation against EC measurements of latent heat flux at the Dadongshu site in the Heihe River Basin. Data were collected from 11:00 to 17:00 and averaged over 30-min intervals on all days between November 1, 2014, and January 31, 2015, with fully covered snow.


Fig. 7

The difference of P–M (P–M–EC) and BA (BA–EC) estimates of snow sublimation against EC measurements of latent heat flux at Dadongshu site. Data were collected from 11:00 to 17:00 and averaged over 30-min intervals on all days between November 1, 2014, and January 31, 2015, with fully covered snow.


An overall impression of the accuracy of both P–M and BA estimated snow sublimation is obtained by considering the entire time series of observations (Fig. 8). The mean daytime snow sublimation estimated with P–M was 0.20  mmd1, while it was 0.16  mmd1 with BA, against 0.15  mmd1 with the EC measurements.

Fig. 8

Time series of daytime snow sublimation: EC measurements of latent heat flux; P–M [Eq. (2)] and BA [Eq. (13)]; data were collected during daytime (11:00 to 17:00) from November 1, 2014, till January 31, 2015, at the Dadongshu site.



Sensitivity Analysis

The sensitivity of estimated sublimation to input data and parameters was evaluated to identify the most influential ones and understand how to deal with data of different accuracies. The sensitivity of estimated snow sublimation to net shortwave radiation, snow surface temperature, air temperature, wind speed, and relative humidity was evaluated. The sensitivity analysis (SA) was performed on one variable at a time by assuming a range between ±20% and ±10% of error for each input variable, while keeping the others constant.59 A nondimensional sensitivity coefficient was calculated as

Eq. (20)

where SC is the dimensionless sensitivity coefficient, LE and CV are the reference values of latent heat flux and input variable, respectively, ΔLE is the change of estimated sublimation when ΔCV is changed for ±20%.

The reference value of each variable was the mean value between 11:00 and 17:00 during the period November 1, 2014, till January 31, 2015, at the Dadongshu site. The reference values of both input variables and estimated sublimation are given in Table 2. The sensitivity results are shown in Fig. 9. The SA indicates that the wind speed is the most influential variable when using either P–M or BA. The P–M equation is less sensitive, i.e., SC=0.77, to wind speed, however, than BA, i.e., SC=1.2. This is a clear advantage of the P–M equation for large area estimates because the regional wind field is less accurate than in situ measurements.52 Another influential variable is relative humidity with SC=0.81 and 0.52 for BA and P–M respectively. Relative humidity has a negative impact on LE for both methods. The influence of net shortwave radiation on P–M estimates of sublimation is comparable to relative humidity, i.e., SC=0.52. Both approaches showed the least sensitivity to air temperature and surface temperature. These agreed with former studies on the sensitivity of snow sublimation calculation to the error of input variables, which highlighted that snow sublimation was highly sensitive to error of winds speed and relative humidity, while not sensitive to the error of temperature.31,60,61

Table 2

Reference values of input variables and the SC values for the P–M equation (SC_P) and the BA parameterization (SC_B).

Wind speedRelative humidityAir temperatureSurface temperatureNet shortwave radiation
Reference values4.52 (ms1)44.72 (%)12.29 (°C)14.78 (°C)101.67 (Wm2)

Fig. 9

Sensitivity of estimated sublimation using (a) P–M and (b) BA methods.



Regional Estimates of Snow Sublimation with Remote Sensing Data

The evaluation of the P–M and BA methods with in situ measurements shows that the P–M estimation is more accurate. Since it is derived from the surface energy balance equation, the P–M equation is also most suitable for regional analyses based on the multispectral radiometric data collected by spaceborne imaging spectroradiometers. This notwithstanding, we applied both methods for our regional case study. Regional snow sublimation was estimated for the entire upper reaches of the Heihe River Basin by applying retrievals from satellite data, i.e., subpixel snow surface temperature and FSC, in combination with atmospheric forcing data. The regional analyses were carried out under clear sky conditions within the period from November 1, 2014, to January 31, 2015. The snow component surface temperature was obtained using the method described in Sec. 3.1.

The map of estimated sublimation obtained on a fully clear day (i.e., November 8, 2014) over the entire upper reaches of the Heihe River Basin was chosen to demonstrate our method. The FSC was rather high and the subpixel snow surface temperature was in the range of 260 to 273 K (Fig. 10). The spatial distribution of snow sublimation estimated by the P–M model and BA model is shown in Fig. 10. The estimated sublimation over the study area ranges between 0 and 80  Wm2 with the highest frequency values around 0 to 20  Wm2 when applying the P–M equation [Figs. 11(a) and 11(b)]. The BA estimates were similar, although the highest frequency values were slightly lower, i.e., around 0 to 10  Wm2 [Figs. 11(c) and 11(d)].

Fig. 10

(a) FSC and (b) snow component surface temperature on November 8, 2014.


Fig. 11

Snow sublimation estimated by (a) the P–M equation, (b) its frequency distribution, (c) the BA parameterization, and (d) its frequency distribution, on November 8, 2014.


The regional estimates of snow sublimation under clear sky conditions were evaluated against the in situ EC measurements at the Dadongshu and Dashalong sites. In the period November 1, 2014, till January 31, 2015, there were 18 usable remotely sensed measurements at Dadongshu and 17 measurements at Dashalong. The P–M estimated sublimation at Dashalong was in good agreement with the measurements (Fig. 11) with R2=0.75, RMSE=8.4  Wm2, and MRE=19.8%. The performance of the P–M method was poorer at Dadongshu (Fig. 12) with R2=0.36, RMSE=9.1  Wm2, and MRE=28.7%. The BA estimates of sublimation were poorer at both Dashalong and Dadongshu (Fig. 13), with R2=0.65, RMSE=17.5  Wm2, and MRE=60.0% at Dashalong and R2=0.06, RMSE=21.2  Wm2, and MRE=35.7% at Dadongshu. Overall the P–M equation performed better than BA, while both performed better at Dashalong than at Dadongshu.

Fig. 12

Evaluation of the P–M equation against EC measurements at (a) the Dashalong site and (b) the Dadongshu site for clear sky days between November 1, 2014, and January 31, 2015.


Fig. 13

Evaluation of the BA parameterization against EC measurements at (a) the Dashalong site and (b) the Dadongshu site for clear sky days between November 1, 2014, and January 31, 2015.


It is important to note here that the EC measurements are essentially local-scale measurements representing only very small spatial scales when compared with the 1-km pixel size of the snow sublimation estimated with satellite data. These observations were considered; however, the best option available for evaluating the model performance because we lacked the consistent and long-term snow sublimation data measured by other methods.



The daily mean value of snow sublimation measured by EC from November 1, 2014, till January 31, 2015, at Dadongshu was 0.15  mmd1, compared with 0.41  mmd1 in the Rocky Mountains of Colorado by Molotch et al.62 and 0.16  mmd1 in Mongolia by Zhang et al.10 The magnitude of sublimation has been shown to vary widely across different land surface environments and elevations.5,10,11,62 The different values of snow sublimation may be explained by environmental differences among study sites.5 Our estimates of daily sublimation are well within the bounds of previous research.

Both the P–M and the BA methods performed well at the Dadongshu site. The lower sensitivity to the wind speed gives the P–M method an enormous advantage when applied to regional scales considering the scarce knowledge of the large area wind field at the 1  km×1  km spatial resolution. Despite the ability of the P–M method to capture the spatial and temporal pattern of snow sublimation distribution, there are multiple uncertainties to deal with. First, we did not account for sublimation from blowing snow. Pomeroy and Essery63 found that blowing snow sublimation has been estimated to account for as much as 10% to 50% of snowfall in North American Prairies and Arctic environments. The occurrence of blowing snow in the upper reaches of the Heihe River Basin has most likely significant impacts on the estimates of sublimation during the snow season. The actual amount of sublimation in this catchment could reasonably be higher than what has been reported here.

The second source of error is the accuracy of the spatial input data. These errors can result from both inaccurate forcing data and remote sensing data. Data on hourly 2-m air temperature, 2-m relative humidity, and snow surface temperature were in good agreement with tower measurements, while estimates of 10-m wind speed correlated poorly with in situ measurements (Fig. 14). The latter is a major challenge in regional analyses of sublimation with either method, given the poor accuracy of the large area wind field.64 Errors on relative humidity, air temperature, and snow surface temperature are smaller but not negligible.

Fig. 14

Evaluation of (a) wind speed, (b) relative humidity, (c) air temperature, and (d) snow temperature from remote sensing data and atmospheric forcing data at Dadongshu and Dashalong sites during the satellite overpass time between November 1, 2014, and January 31, 2015.




Snow sublimation plays a significant role in hydrological processes in cold regions and seasons, and it varies with atmospheric and surface conditions. In this study, we developed a method to estimate snow surface sublimation with satellite observations and validated the results with EC measurements of latent heat flux. We applied the P–M equation and the BA parameterization, which use satellite data products to estimate regional snow sublimation at a 1-km spatial resolution. Model estimates have been compared with ground EC measurements. The P–M estimates were in good agreement with measurements at two sites (Dashalong and Dadongshu), with RMSE=8.4 and 9.1  Wm2, respectively. The BA estimates were less accurate. The encouraging performance of the P–M suggested the extension of this approach to large area estimates of snow sublimation based on multiple satellite data. Remote sensing data at moderate spatial resolution enable us to estimate and map snow sublimation over large areas. The applicability and accuracy of the P–M method to estimate snow sublimation relies on the availability and accuracy of a suite of remotely sensed input data products, i.e., snow surface temperature, FSC, and atmospheric forcing data. Anyway, remote sensing has a large advantage to get the snow information with large-scale supervising ability. However, there are still some challenges. For example, clouds will greatly influence the quality of remote sensing data such as snow surface temperature and FSC. It is important to realize that the present study is limited by the unavailability of clear sky remote sensing data. As such, higher quality and accuracy of these input data are vital and expected.


The authors declare no conflict of interest.


This work was supported by the National Key Basic Research Program of China (Grant No. 2015CB953702), the National Natural Science Foundation of China (Grant Nos. 91425303 and 91325203), and the SAFEA Long-Term-Projects of the 1000 Talent Plan for High-Level Foreign Experts (Grant No. WQ20141100224). We gratefully acknowledge the Cold and Arid Regions Sciences Data Center at Lanzhou for providing EC data and atmospheric forcing data. We are grateful to the anonymous reviewer for the useful suggestions.



O. Schulz and C. de Jong, “Snowmelt and sublimation: field experiments and modelling in the High Atlas Mountains of Morocco,” Hydrol. Earth Syst. Sci., 8 (6), 1076 –1089 (2004). Google Scholar


K. Suzuki et al., “Estimation of continental-basin-scale sublimation in the Lena River Basin, Siberia,” Adv. Meteorol., 2015 1 –14 (2015). Google Scholar


G. Zhang et al., “Snow cover dynamics of four lake basins over Tibetan Plateau using time series MODIS data (2001–2010),” Water Resour. Res., 48 (10), W10529 (2012). WRERAQ 0043-1397 Google Scholar


Y. S. Zhang et al., “Sublimation from snow surface in southern mountain taiga of eastern Siberia,” J. Geophys. Res., 109 D21103 (2004). Google Scholar


U. Strasser et al., “Is snow sublimation important in the alpine water balance?,” Cryosphere, 2 (1), 53 –66 (2008). Google Scholar


M. K. MacDonald et al., “On the importance of sublimation to an alpine snow mass balance in the Canadian Rocky Mountains,” Hydrol. Earth Syst., 14 (7), 1401 –1415 (2010). Google Scholar


M. L. Reba et al., “Estimating surface sublimation losses from snowpacks in a mountain catchment using eddy covariance and turbulent transfer calculations,” Hydrol. Processes, 26 (24), 3699 –3711 (2012). HYPRE3 1099-1085 Google Scholar


E. Hood et al., “Sublimation from a seasonal snowpack at a continental, mid-latitude alpine site,” Hydrol. Processes, 13 (12–13), 1781 –1797 (1999). HYPRE3 1099-1085 Google Scholar


D. Marks and J. Dozier, “Climate and energy exchange at the snow surface in the Alpine region of the Sierra Nevada: 2. Snow cover energy balance,” Water Resour. Res., 28 (11), 3043 –3054 (1992). WRERAQ 0043-1397 Google Scholar


Y. S. Zhang et al., “Sublimation from thin snow cover at the edge of the Eurasian cryosphere in Mongolia,” Hydrol. Processes, 22 (18), 3564 –3575 (2008). HYPRE3 1099-1085 Google Scholar


G. A. Sexstone et al., “Comparison of methods for quantifying surface sublimation over seasonally snow-covered terrain,” Hydrol. Processes, 30 (19), 3373 –3389 (2016). HYPRE3 1099-1085 Google Scholar


D. Cline, “Snow surface energy exchanges and snowmelt at a continental, mid-latitude alpine site,” Water Resour. Res., 33 (4), 689 –701 (1997). WRERAQ 0043-1397 Google Scholar


R. A. Schmidt et al., “Sublimation of snowpack in subalpine conifer forests,” Can. J. For. Res., 28 (4), 501 –513 (1998). CJFRAR 0045-5067 Google Scholar


J. E. Box and K. Steffen, “Sublimation on the Greenland ice sheet from automated weather station observations,” J. Geophys. Res., 106 (D24), 33965 –33981 (2001). JGREA2 0148-0227 Google Scholar


H. Y. Li et al., “Influence of blowing snow on snow mass and energy exchange in the Qilian mountainous,” J. Glaciol. Geocryol., 34 (5), 1084 –1090 (2012). Google Scholar


D. Marks et al., “Comparing simulated and measured sensible and latent heat fluxes over snow under a pine canopy,” J. Hydrometeorol., 9 (6), 1506 –1522 (2008). Google Scholar


M. L. Reba et al., “An assessment of corrections for eddy covariance measured turbulent fluxes over snow in mountain environments,” Water Resour. Res., 45 W00D38 (2009). WRERAQ 0043-1397 Google Scholar


D. Bewley et al., “Measurements and modeling of snowmelt and turbulent heat fluxes over shrub tundra,” Hydrol. Earth Syst. Sci., 14 (7), 1331 –1340 (2010). Google Scholar


W. Yang et al., “Summertime surface energy budget and ablation modeling in the ablation zone of a maritime Tibetan glacier,” J. Geophys. Res., 116 D14116 (2011). JGREA2 0148-0227 Google Scholar


N. N. Chen et al., “Influences of snow event on energy balance over temperature meadow in dormant season based on eddy covariance measurements,” J. Hydrol., 399 (1), 100 –107 (2011). JHYDA7 0022-1694 Google Scholar


W. Helgason and J. Pomeroy, “Problems closing the energy balance over a homogeneous snow cover during midwinter,” J. Hydrometeorol., 13 (13), 557 –572 (2012). Google Scholar


T. E. Twine et al., “Correcting eddy-covariance flux underestimates over a grassland,” Agric. For. Meteorol., 103 (3), 279 –300 (2000). 0168-1923 Google Scholar


K. Wilson et al., “Energy balance closure at FLUXNET sites,” Agric. For. Meteorol., 113 (1), 223 –243 (2002). 0168-1923 Google Scholar


L. S. Kuchment and A. N. Gelfan, “The determination of the snowmelt rate and the meltwater outflow from a snowpack for modelling river runoff generation,” J. Hydrol., 179 (1–4), 23 –36 (1996). JHYDA7 0022-1694 Google Scholar


J. C. King et al., “The seasonal cycle of sublimation at Halley, Antarctica,” J. Glaciol., 47 (156), 1 –8 (2001). JOGLAO 0022-1430 Google Scholar


A. N. Gelfan et al., “Modeling forest cover influences on snow accumulation, sublimation, and melt,” J. Hydrometeorol., 5 (5), 785 –803 (2004).<0785:MFCIOS>2.0.CO;2 Google Scholar


R. Moore, “On the use of bulk aerodynamic formulas over melting snow,” Nord. Hydrol., 14 (4), 193 –206 (1983). NOHYBB 0029-1277 Google Scholar


W. P. Kustas et al., “A simple energy budget algorithm for the snowmelt runoff model,” Water Resour. Res., 30 (5), 1515 –1527 (1994). WRERAQ 0043-1397 Google Scholar


L. Mahrt and D. Vickers, “Moisture fluxes over snow with and without protruding vegetation,” Q. J. R. Meteorol. Soc., 131 (607), 1251 –1270 (2005). QJRMAM 0035-9009 Google Scholar


D. R. Dewalle and A. Rango, Principles of Snow Hydrology, 159 –166 Cambridge University Press, Cambridge (2008). Google Scholar


A. Lundberg and S. Halldin, “Evaporation of intercepted snow: analysis of governing factors,” Water Resour. Res., 30 (9), 2587 –2598 (1994). WRERAQ 0043-1397 Google Scholar


A. Lundberg et al., “Evaporation of intercepted snow: measurement and modeling,” J. Hydrol., 206 (3–4), 151 –163 (1998). JHYDA7 0022-1694 Google Scholar


J. L. Monteith, “Evaporation and environment,” in Symp. of the Society for Experimental Biology, 205 –234 (1965). Google Scholar


R. Dhungel et al., “Comparison of latent heat flux using aerodynamic methods and using the Penman–Monteith method with satellite-based surface energy balance,” Remote Sens., 6 (6), 8844 –8877 (2014). Google Scholar


J. Yu et al., “Developing daily cloud-free snow composite products from MODIS Terra-Aqua and IMS for the Tibetan Plateau,” IEEE Trans. Geosci. Remote Sens., 54 (4), 2171 –2180 (2016). IGRSD2 0196-2892 Google Scholar


D. K. Hall et al., “Comparison of satellite-derived and in-situ observations of ice and snow surface temperatures over Greenland,” Remote Sens. Environ., 112 (10), 3739 –3749 (2008). Google Scholar


V. V. Salomonson and I. Appel, “Estimating fractional snow cover from MODIS using the normalized difference snow index,” Remote Sens. Environ., 89 (3), 351 –360 (2004). Google Scholar


C. O. Justice et al., “An overview of MODIS Land data processing and product status,” Remote Sens. Environ., 83 (1), 3 –15 (2002). Google Scholar


Q. Mu et al., “Improvements to a MODIS global terrestrial evapotranspiration algorithm,” Remote Sens. Environ., 115 (8), 1781 –1800 (2011). Google Scholar


Y. H. Wang et al., “Impact of cryosphere hydrological processes on the river runoff in the upper reaches of Heihe River,” J. Hydraul. Eng., 46 (9), 1064 –1072 (2015). JHEND8 0733-9429 Google Scholar


Z. Chen et al., “Environmental isotopic study on the recharge and residence time of groundwater in the Heihe River Basin, northwestern china,” Hydrogeol. J., 14 (8), 1635 –1651 (2006). Google Scholar


J. Wang et al., “Responses of snowmelt runoff to climatic change in an inland river basin, Northwestern China, over the past 50 years,” Hydrol. Earth Syst. Sci., 14 (10), 1979 –1987 (2010). Google Scholar


T. Che et al., “Estimation of snow depth and snow water equivalent distribution using airborne microwave radiometry in the Binggou watershed, the upper reaches of the Heihe River basin,” Int. J. Appl. Earth Obs. Geoinf., 17 (1), 23 –32 (2012). Google Scholar


S. M. Liu et al., “A comparison of eddy-covariance and large aperture scintillometer measurements with respect to the energy balance closure problem,” Hydrol. Earth Syst. Sci., 15 (4), 1291 –1306 (2011). Google Scholar


X. Li et al., “Heihe watershed allied telemetry experimental research (HiWATER) scientific objectives and experimental design,” Bull. Am. Meteorol. Soc., 94 (8), 1145 –1160 (2013). BAMIAT 0003-0007 Google Scholar


H. Y. Li and J. Wang, “Simulation of snow distribution and melt under cloudy conditions in an Alpine watershed,” Hydrol. Earth Syst. Sci., 15 (7), 2195 –2203 (2011). Google Scholar


Y. H. Ran et al., “MICLCover land cover map of the Heihe River basin,” (2011). Google Scholar


Y. H. Ran et al., “Large-scale land cover mapping with the integration of multi-source information based on the Dempster–Shafer theory,” Int. J. Geogr. Inf. Sci., 26 (1), 169 –191 (2012). Google Scholar


X. D. Pan and X. Li, “Validation of WRF model on simulating forcing data for Heihe River Basin,” Sci. Cold Arid Reg., 3 (4), 344 –357 (2011). Google Scholar


X. D. Pan et al., “Dynamic downscaling of near-surface air temperature at the basin scale using WRF—a case study in the Heihe River Basin, China,” Front. Earth Sci., 6 (3), 314 –323 (2012). Google Scholar


X. D. Pan and X. Li, “The atmospheric forcing data from 2000 to 2015 in the Heihe river basin,” (2013). Google Scholar


G. C. Hu and L. Jia, “Monitoring of evapotranspiration in a semi-arid inland river basin by combining microwave and optical remote sensing observations,” Remote Sens., 7 (3), 3056 –3087 (2015). Google Scholar


H. L. Penman, “Natural evaporation from open water, bare soil, and grass,” Proc. R. Soc. London, 193 (1032), 120 –145 (1948). PRSLAZ 0370-1662 Google Scholar


M. Menenti, “Physical aspects and determination of evaporation in deserts applying remote sensing techniques,” 27 –34 Wageningen (1984). Google Scholar


Z. Q. Sun et al., “Methods of retrieving sub-pixel snow surface temperature based on MODIS data,” J. Arid Land Resour. Environ., 28 (5), 166 –171 (2014). Google Scholar


D. Baker et al., “The albedo decay of prairie snows,” J. Appl. Meteorol., 29 (2), 179 –187 (1990).<0179:TADOPS>2.0.CO;2 Google Scholar


H. Douville et al., “A new snow parameterization for the Météo-France Climate Model. 1. Validation in stand-alone experiments,” Clim. Dyn., 12 (1), 21 –35 (1995). CLDYEM 0930-7575 Google Scholar


A. J. Chen et al., “An observational study of snow aging and the seasonal variation of snow albedo by using data from Col de Oorte, France,” Chin. Sci. Bull., 59 (34), 4881 –4889 (2014). CSBUEF 1001-6538 Google Scholar


C. L. Zheng and Q. Wang, “Spatiotemporal variations of reference evapotranspiration in recent five decades in the arid land of Northwestern China,” Hydrol. Processes, 28 (25), 6124 –6134 (2014). HYPRE3 1099-1085 Google Scholar


R. A. Schmidt, “Sublimation of snow intercepted by an artificial conifer,” Agric. For. Meteorol., 54 (1), 1 –27 (1991). 0168-1923 Google Scholar


Y. Nakai, “An observational study on evaporation from intercepted snow on Forest Canopies,” 107 Kyoto (1996). Google Scholar


N. P. Molotch et al., “Estimating sublimation of intercepted and sub-canopy snow using eddy covariance systems,” Hydrol. Processes, 21 (12), 1567 –1575 (2007). HYPRE3 1099-1085 Google Scholar


J. W. Pomeroy and R. L. H. Essery, “Turbulent fluxes during blowing snow: field tests of model sublimation predictions,” Hydrol. Processes, 13 (18), 2963 –2975 (1999). HYPRE3 1099-1085 Google Scholar


R. Dadic et al., “Sensitivity of turbulent fluxes to wind speed over snow surfaces in different climatic setting,” Adv. Water Resour., 55 (3), 178 –189 (2013). AWREDI 0309-1708 Google Scholar


Ning Wang is a PhD candidate at the University of Chinese Academy of Sciences, China. He received his master’s degree in GIS from Guilin University of Technology in 2012. His research interests include remote sensing of snow and image processing.

Li Jia received her PhD in environmental science from Wageningen University of the Netherlands in 2004. She is a professor at the State Key Laboratory of Remote Sensing Science, jointly sponsored by the Institute of Remote Sensing and Digital Earth of Chinese Academy of Sciences and Beijing Normal University. Her research interests are on the study of earth observation and its applications in hydrometeorology, water resources, agriculture, and climate change.

Chaolei Zheng received his PhD in environment energy system at Shizuoka University, Japan, in 2013. Currently, he is an assistant professor at the Institute of Remote Sensing and Digital Earth, Chinese Academy of Science. His research focuses on earth observations of the water cycle, climate change, ecohydrology, and remote sensing.

Massimo Menenti received his PhD from Wageningen Agriculture University of The Netherlands in 1984. He is a full professor at Delft University of Technology of The Netherlands. He is also a professor at the Institute of Remote Sensing and Digital Earth of Chinese Academy. His research interests focus on the use of earth observation to study the hydrology and hydrometeorology of the global land surface.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Ning Wang, Li Jia, Chaolei Zheng, and Massimo Menenti "Estimation of subpixel snow sublimation from multispectral satellite observations," Journal of Applied Remote Sensing 11(4), 046017 (8 December 2017).
Received: 8 May 2017; Accepted: 27 October 2017; Published: 8 December 2017

Back to Top