Snow Mapping Research Articles

On the Automated Mapping of Snow Cover on Glaciers and Calculation of Snow Line Altitudes from Multi-Temporal Landsat Data

Mapping snow cover (SC) on glaciers at the end of the ablation period provides a possibility to rapidly obtain a proxy for their equilibrium line altitude (ELA) which in turn is a metric for the mass balance. Satellite determination of glacier snow cover, derived over large regions, can reveal its spatial variability and temporal trends. Accordingly, snow mapping on glaciers has been widely applied using several satellite sensors. However, as glacier ice originates from compressed snow, both have very similar spectral properties and standard methods to map snow struggle to distinguish snow on glaciers. Hence, most studies applied manual delineation of snow extent on glaciers. Here we present an automated tool, named ‘ASMAG’ (automated snow mapping on glaciers), to map SC on glaciers and derive the related snow line altitude (SLA) for individual glaciers using multi-temporal Landsat satellite imagery and a digital elevation model (DEM). The method has been developed using the example of the Ötztal Alps, where an evaluation of the method is possible using field-based observations of the annual equilibrium line altitude (ELA) and the accumulation area ratio (AAR) measured for three glaciers for more than 30 years. The tool automatically selects a threshold to map snow on glaciers and robustly calculates the SLA based on the frequency distribution of elevation bins with more than 50% SC. The accuracy of the SC mapping was about 90% and the SLA was determined successfully in 80% of all cases with a mean uncertainty of ±19 m. When cloud-free scenes close to the date of the highest snowline are available, a good to very good agreement of SC ratios (SCR)/SLA with field data of AAR/ELA are obtained, otherwise values are systematically higher/lower as useful images were often acquired too early in the summer season. However, glacier specific differences are still well captured. Snow mapping on glaciers is impeded by clouds and their shadows or when fresh snow is covering the glaciers, so that more frequent image acquisitions (as now provided by Sentinel-2) would improve results.

A Model Setup for Mapping Snow Conditions in High-Mountain Himalaya

Seasonal snow cover is an important source of melt water for irrigation and hydropower production in many regions of the world, but can also be a cause of disasters, such as avalanches and floods. In the remote Himalayan environment there is a great demand for up-to-date information on the snow conditions for the purposes of planned hydropower development and disaster risk reduction initiatives. We describe and evaluate a snow mapping setup for the remote Langtang Valley in the Nepal Himalayas, which can deliver data for snow and water availability mapping all year round. The setup utilizes (1) robust and almost maintenance-free in-situ instrumentation with satellite transmission, (2) a freely available numerical snow model, and (3) estimation of model key parameters from local meteorological and snow observations as well as from freely available climatological data. Novel features in the model include the estimation of melt parameters and solid precipitation from passive gamma-radiation based snow sensor data, as well as improved parameterization and estimation of melt water refreezing (36 % of total melt) within, and sublimation/evaporation (57 mm yr-1) from the snow pack. Evaluation of the model results show a reasonable fit with snow cover data from satellite images. As many of the high-mountain regions in central and eastern Nepal show high correlation (>0.8) with the estimated snow line elevation in the Langtang catchment, the results may provide a first-order approximation of the snow conditions for these areas too.

Dry/wet snow mapping based on the synergistic use of dual polarimetric SAR and multispectral data

We propose a multi-sensor multi-spectral and bi-temporal dual-polarimetric Synthetic Aperture Radar (SAR) data integration scheme for dry/wet snow mapping using Sentinel-2 and Sentinel-1 data which are freely available to the research community. The integration is carried out by incorporating the information retrieved from ratio images of the conventional method for wet snow mapping and the multispectral data in two different frameworks. Firstly, a simple differencing scheme is employed for dry/wet snow mapping, where the snow cover area is derived using the Normalized Differenced Snow Index (NDSI). In the second framework, the ratio images are stacked with the multispectral bands and this stack is used for supervised and unsupervised classification using support vector machines for dry/wet snow mapping. We also investigate the potential of a state of the art backscatter model for the identification of dry/wet snow using Sentinel-1 data. The results are validated using a reference map derived from RADARSAT-2 full polarimetric SAR data. A good agreement was observed between the results and the reference data with an overall accuracy greater than 0.78 for the different blending techniques examined. For all the proposed frameworks, the wet snow was better identified. The coefficient of determination between the snow wetness derived from the backscatter model and the reference based on RADARSAT-2 data was observed to be 0.58 with a significantly higher root mean square error of 1.03% by volume.

Cross-Comparison between MODIS and VIIRS Snow Cover Products for the 2016 Hydrological Year

The VIIRS (Visible Infrared Imaging Radiometer Suite) instrument on board the Suomi-NPP (National Polar-Orbiting Partnership) satellite aims to provide long-term continuity of several environmental data series including snow cover initiated with MODIS (Moderate Resolution Imaging Spectroradiometer). Although it is speculated that MODIS and VIIRS snow cover products may differ because of their differing spatial resolutions and spectral coverage, quantitative comparisons between their snow products are currently limited. Therefore, this study intercompares MODIS and VIIRS snow products for the 2016 Hydrological Year over the Midwestern United States and southern Canada. Two hundred and forty-four swath snow products from MODIS/Aqua (MYD10L2) and the VIIRS EDR (Environmental Data Records) (VSCMO/binary) were intercompared using confusion matrices, comparison maps and false color imagery. Thresholding the MODIS NDSI (Normalized Difference Snow Index) Snow Cover product at a snow cover fraction of 30% generated binary snow maps are most comparable to the NOAA VIIRS binary snow product. Overall agreement between MODIS and VIIRS was found to be approximately 98%. This exceeds the VIIRS accuracy requirements of 90% probability of correct typing. The agreement was highest during the winter but lower during late fall and spring. MODIS and VIIRS often mapped snow/no-snow transition zones as a cloud. The assessment of total snow and cloud pixels and comparison snow maps of MODIS and VIIRS indicate that VIIRS is mapping more snow cover and less cloud cover compared to MODIS. This is evidenced by the average area of snow in MYD10L2 and VSCMO being 5.72% and 11.43%, no-snow 26.65% and 28.67% and cloud 65.02% and 59.91%, respectively. While VIIRS and MODIS have a similar capacity to map snow cover, VIIRS has the potential to map snow cover area more accurately, for the successful development of climate data records.

Theia Snow collection: high-resolution operational snow cover maps from Sentinel-2 and Landsat-8 data

Abstract. The Theia Snow collection routinely provides high-resolution maps of the snow-covered area from Sentinel-2 and Landsat-8 observations. The collection covers selected areas worldwide, including the main mountain regions in western Europe (e.g. Alps, Pyrenees) and the High Atlas in Morocco. Each product of the Theia Snow collection contains four classes: snow, no snow, cloud and no data. We present the algorithm to generate the snow products and provide an evaluation of the accuracy of Sentinel-2 snow products using in situ snow depth measurements, higher-resolution snow maps and visual control. The results suggest that the snow is accurately detected in the Theia snow collection and that the snow detection is more accurate than the Sen2Cor outputs (ESA level 2 product). An issue that should be addressed in a future release is the occurrence of false snow detection in some large clouds. The snow maps are currently produced and freely distributed on average 5 d after the image acquisition as raster and vector files via the Theia portal (https://doi.org/10.24400/329360/F7Q52MNK).

Monitoring of Snow Cover Ablation Using Very High Spatial Resolution Remote Sensing Datasets

This study tested the potential of a short time series of very high spatial resolution (cm to dm) remote sensing datasets obtained from unmanned aerial system (UAS)-based photogrammetry and terrestrial laser scanning (TLS) to monitor snow cover ablation in the upper Dischma valley (Davos, Switzerland). Five flight missions (for UAS) and five scans (for TLS) were carried out simultaneously: Four during the snow-covered period (9, 10, 11, and 27 May 2016) and one during the snow-free period (24 June 2016 for UAS and 31 May 2016 for TLS). The changes in both the areal extent of the snow cover and the snow depth (HS) were assessed together in the same case study. The areal extent of the snow cover was estimated from both UAS- and TLS-based orthophotos by classifying pixels as snow-covered and snow-free based on a threshold value applied to the blue band information of the orthophotos. Also, the usage possibility of TLS-based orthophotos for mapping snow cover was investigated in this study. The UAS-based orthophotos provided higher overall classification accuracy (97%) than the TLS-based orthophotos (86%) and allowed for mapping snow cover in larger areas than the ones from TLS scans by preventing the occurrence of gaps in the orthophotos. The UAS-based HS were evaluated and compared to TLS-based HS. Initially, the CANUPO (CAractérisation de NUages de POints) binary classification method, a proposed approach for improving the quality of models to obtain more accurate HS values, was applied to the TLS 3D raw point clouds. In this study, the use of additional artificial ground control points (GCPs) was also proposed to improve the quality of UAS-based digital elevation models (DEMs). The UAS-based HS values were mapped with an error of around 0.1 m during the time series. Most pixels representing change in the HS derived from the UAS data were consistent with the TLS data. The time series used in this study allowed for testing of the significance of the data acquisition interval in the monitoring of snow ablation. Accordingly, this study concluded that both the UAS- and TLS-based high-resolution DSMs were biased in detecting change in HS, particularly for short time spans, such as a few days, where only a few centimeters in HS change occur. On the other hand, UAS proved to be a valuable tool for monitoring snow ablation if longer time intervals are chosen.

Assessment of MODIS-Based Fractional Snow Cover Products Over the Tibetan Plateau

The snow cover distribution is crucial to Earth's climate and hydrological systems. Snow cover maps from Moderate Resolution Imaging Spectroradiometer (MODIS) are widely utilized in climate research and water-source management. Three MODIS-derived fractional snow cover (FSC) products, namely, MOD10A1 version 6, MODSCAG, and MODAGE, are assessed to evaluate the accuracy of FSC data over the Tibetan Plateau from May 2013 to April 2015. A Landsat-8/ Operational Land Imager (OLI) snow map that is based on linear spectral mixture analysis is regarded as the “reference dataset.” A total of 149 OLI images that span the Tibetan Plateau are used. The performance of these three datasets is compared to the “reference dataset” by using binary and fractional metrics over forest, grass, and bare soil, and in the Himalaya region. The binary classification assessment reveals that larger local illumination angles and vegetation coverage at the snow boundaries resulted in omission errors in MOD10A1, and snow patchiness caused this model to miss snow. MODSCAG correctly identified snow pixels, but overestimated snow at the snow boundaries and at local illumination angles that exceeded 30°. MODAGE precisely recognized snow-free areas, but underestimated snow coverage because of snow patchiness. MODSCAG and MODAGE produced smaller root mean square errors (RMSE) than MOD10A1, likely because spectral mixture analysis is superior to the normalized difference snow index-based empirical method. The behaviors of these products differed at various spatial scales. The RMSE significantly declined with decreasing spatial resolution from 500 m to 2 km. However, when aggregated to 5 km, the RMSE increased because fewer pixels were involved in the calculation.

Multi-Source Based Spatio-Temporal Distribution of Snow in a Semi-Arid Headwater Catchment of Northern Mongolia

Knowledge of the duration and distribution of seasonal snow cover is important for understanding the hydrologic regime in mountainous regions within semi-arid climates. In the headwater of the semi-arid Sugnugur catchment (in the Khentii Mountains, northern Mongolia), a spatial analysis of seasonal snow cover duration (SCD) was performed on a 30 m spatial resolution by integrating the spatial resolution of Landsat-7, Landsat-8, and Sentinel-2A images with the daily temporal resolution of Moderate Resolution Imaging Spectroradiometer (MODIS) snow products (2000–2017). Validation was achieved using in situ time series measurements from winter field campaigns and distributed surface temperature loggers. We found a mean increase of SCD with altitude at approximately +6 days/100 m. However, we found no altitude-dependent changes in snow depth during field campaigns. The southern exposed valley slopes are either snow free or covered by intermittent snow throughout the winter months due to high sublimation rates and prevailing wind. The estimated mean SCD ranges from 124 days in the lower parts of the catchment to 226 days on the mountain peaks, with a mean underestimation of 12–13 days. Snow onset and melt dates exhibited large inter-annual variability, but no significant trend in the seasonal SCD was evident. This method can be applied to high-resolution snow mapping in similar mountainous regions.

Integrating camera imagery, crowdsourcing, and deep learning to improve high-frequency automated monitoring of snow at continental-to-global scales.

Snow is important for local to global climate and surface hydrology, but spatial and temporal heterogeneity in the extent of snow cover make accurate, fine-scale mapping and monitoring of snow an enormous challenge. We took 184,453 daily near-surface images acquired by 133 automated cameras and processed them using crowdsourcing and deep learning to determine whether snow was present or absent in each image. We found that the crowdsourced data had an accuracy of 99.1% when compared with expert evaluation of the same imagery. We then used the image classification to train a deep convolutional neural network via transfer learning, with accuracies of 92% to 98%, depending on the image set and training method. The majority of neural network errors were due to snow that was present not being detected. We used the results of the neural networks to validate the presence or absence of snow inferred from the MODIS satellite sensor and obtained similar results to those from other validation studies. This method of using automated sensors, crowdsourcing, and deep learning in combination produced an accurate high temporal dataset of snow presence across a continent. It holds broad potential for real-time large-scale acquisition and processing of ecological and environmental data in support of monitoring, management, and research objectives.

Validating improved-MODIS products from spectral mixture-Landsat snow cover maps in a mountain region in southern Spain

Abstract. Remote sensing is the only feasible data source for distributed modelling of snow in mountain regions on medium to large scales, due to the limited access to these areas together with the lack of dense ground monitoring stations for snow variables. Observations worldwide identify snow cover persistence together with snowfall occurrence as the most affected variables by global warming. In Mediterranean regions, the spatiotemporal evolution of the snow cover can experiment quick changes that result in different accumulation-ablation cycles during the cold season. High frequency sensors are required to adequately monitor such shifts; however, for trend analyses, the Landsat time series constitute the only available source of data, being their frequency low for this regime, especially when cloudy conditions limit the available images. On the other hand, the MODIS daily series provide more than 15 years of continuous snow maps, despite the spatial resolution may pose a constraint in areas with abrupt topography; several approaches have been done to improve their spatial resolution from combining different information. This work presents a methodological approach to validate the improved MODIS daily snow cover maps from Notarnicola et al. (2013a, b), with 250 m spatial resolution, in Sierra Nevada (southern Spain), from a reference data set obtained by spectral mixture analyses of Landsat TM data by Pimentel et al. (2017b). This reference time series of fractional snow maps, with 30 m spatial resolution, were validated from high resolution local time series of snow maps obtained by terrestrial time-lapse cameras. The results show a significantly high correlation between the two snow map products both on a global and basin scales in the Sierra Nevada area. Selected areas and time periods are shown to address the convergence and divergence between both products and assess the development of a fusion algorithm to retrieve daily Landsat-resolution snow maps on a long term basis.

Reflectance variation in boreal landscape during the snow melting period using airborne imaging spectroscopy

We aim a better understanding of the effect of spring-time snow melt on the remotely sensed scene reflectance by using an extensive amount of optical spectral data obtained from an airborne hyperspectral campaign in Northern Finland. We investigate the behaviour of thin snow reflectance for different land cover types, such as open areas, boreal forests and treeless fells. Our results not only confirm the generally known fact that the reflectance of a melting thin snow layer is considerably lower than that of a thick snow layer, but we also present analyses of the reflectance variation over different land covers and in boreal forests as a function of canopy coverage. According to common knowledge, the highly variating reflectance spectra of partially transparent, most likely also contaminated thin snow pack weakens the performance of snow detection algorithms, in particular in the mapping of Fractional Snow Cover (FSC) during the end of the melting period. The obtained results directly support further development of the SCAmod algorithm for FSC retrieval, and can be likewise applied to develop other algorithms for optical satellite data (e.g. spectral unmixing methods), and to perform accuracy assessments for snow detection algorithms.A useful part of this work is the investigation of the competence of Normalized Difference Snow Index (NDSI) in snow detection in late spring, since it is widely used in snow mapping. We conclude, based on the spectral data analysis, that the NDSI -based snow mapping is more accurate in open areas than in forests. However, at the very end of the snow melting period the behavior of the NDSI becomes more unstable and unpredictable in non-forests with shallow snow, increasing the inaccuracy also in non-forested areas. For instance in peatbogs covered by melting snow layer (snow depth < 30 cm) the mean NDSI -0.6 was observed, having coefficient of variation as high as 70%, whereas for deeper snow packs the mean NDSI shows positive values.

SNOW AND CLOUD DISCRIMINATION USING CONVOLUTIONAL NEURAL NETWORKS

Abstract. Snow is an important feature on our planet, and measuring its extent has advantages in climate studies. Snow mapping through satellite remote sensing is often affected by cloud cover. This issue can be resolved by using short wave infrared (SWIR) sensors. In order to obtain an effective cloud mask, our study aims to use SWIR data of a ResourceSat-2 satellite. We employ Convolutional Neural Networks (CNN) to discriminate similar pixels of clouds and snow. The technique is expected to give a high accuracy compared to traditional methods such as thresholding. The cloud mask thus produced can also be used for creating the metadata for Indian Remote Sensing products.

Repeat mapping of snow depth across an alpine catchment with RPAS photogrammetry

Abstract. Being dynamic in time and space, seasonal snow represents a difficult target for ongoing in situ measurement and characterisation. Improved understanding and modelling of the seasonal snowpack requires mapping snow depth at fine spatial resolution. The potential of remotely piloted aircraft system (RPAS) photogrammetry to resolve spatial variability of snow depth is evaluated within an alpine catchment of the Pisa Range, New Zealand. Digital surface models (DSMs) at 0.15 m spatial resolution in autumn (snow-free reference) winter (2 August 2016) and spring (10 September 2016) allowed mapping of snow depth via DSM differencing. The consistency and accuracy of the RPAS-derived surface was assessed by the propagation of check point residuals from the aero-triangulation of constituent DSMs and via comparison of snow-free regions of the spring and autumn DSMs. The accuracy of RPAS-derived snow depth was validated with in situ snow probe measurements. Results for snow-free areas between DSMs acquired in autumn and spring demonstrate repeatability yet also reveal that elevation errors follow a distribution that substantially departs from a normal distribution, symptomatic of the influence of DSM co-registration and terrain characteristics on vertical uncertainty. Error propagation saw snow depth mapped with an accuracy of ±0.08 m (90 % c.l.). This is lower than the characterization of uncertainties on snow-free areas (±0.14 m). Comparisons between RPAS and in situ snow depth measurements confirm this level of performance of RPAS photogrammetry while also highlighting the influence of vegetation on snow depth uncertainty and bias. Semi-variogram analysis revealed that the RPAS outperformed systematic in situ measurements in resolving fine-scale spatial variability. Despite limitations accompanying RPAS photogrammetry, which are relevant to similar applications of surface and volume change analysis, this study demonstrates a repeatable means of accurately mapping snow depth for an entire, yet relatively small, hydrological catchment (∼0.4 km2) at very high resolution. Resolving snowpack features associated with redistribution and preferential accumulation and ablation, snow depth maps provide geostatistically robust insights into seasonal snow processes, with unprecedented detail. Such data will enhance understanding of physical processes controlling spatial distributions of seasonal snow and their relative importance on varying spatial and temporal scales.

A method for monthly mapping of wet and dry snow using Sentinel-1 and MODIS: Application to a Himalayan river basin

Satellite Remote Sensing, with both optical and SAR instruments, can provide distributed observations of snow cover over extended and inaccessible areas. Both instruments are complementary, but there have been limited attempts at combining their measurements. We describe a novel approach to produce monthly maps of dry and wet snow areas through application of data fusion techniques to MODIS fractional snow cover and Sentinel-1 wet snow mask, facilitated by Google Earth Engine. The method is demonstrated in a 55,000 km2 river basin in the Indian Himalayan region over a period of ∼2.5 years, although it can be applied to any areas of the world where Sentinel-1 data are routinely available. The typical underestimation of wet snow area by SAR is corrected using a digital elevation model to estimate the average melting altitude. We also present an empirical model to derive the fractional cover of wet snow from Sentinel-1. Finally, we demonstrate that Sentinel-1 effectively complements MODIS as it highlights a snowmelt phase which occurs with a decrease in snow depth but no/little decrease in snowpack area. Further developments are now needed to incorporate these high resolution observations of snow areas as inputs to hydrological models for better runoff analysis and improved management of water resources and flood risk.

Mapping Snow for Area Where was Detected Ground Deplacement: A Case Studi of Iaşi County

Abstract The main goal of this paper is to detect snow in areas where was detecting and mapping, using Differential Radar Interferometry (DInSAR) technique, ground displacement. DInSAR is a powerful tool to detect and monitor ground deformation. Iaşi county is considered as study area in this research. Study area is geographically situated on latitude 46°48’N to 47°35’N and longitude 26°29’E to 28°07’E. For this paper, to detect and mapping grond displacement, was used Sentinel – 1 images, provided free by The European Space Agency (ESA), for January 2018, with vertical polarization (VV), ascending orbit and Interferometric Wide swath (IW) mode operated. SNAP was used to process the Sentinel – 1 images. Landsat-8 OLI was taken to detect areas cover with snow using Normalized Difference Snow Index (NDSI) - a numerical indicator that shows snow cover over land areas. ArcMap was used to create NDSI map after Landsat-8 data was preprocessed. The presence of snow has been observed both in the areas where it exists vertical displacement positive and negative.

Snow cover mapped daily at 30 meters resolution using a fusion of multi-temporal MODIS NDSI data and Landsat surface reflectance

Snow fall and melt events are complex meteorological phenomena that help chart the effects of climate change and impact many critical environmental processes including hydrologic and biogeographic systems. Daily snow maps, derived from MODIS imagery, provide managers and researchers with vital snow cover information, but only at spatial scales of 500 m or more. Finer resolution time series maps, however, retain large temporal gaps, particularly during recurrent cloud cover. This paper’s authors have developed the novel algorithm MODSAT-NDSI to harness the strengths of both coarse and finer spatial resolution imagery by fusing MODIS and Landsat normalized difference snow index (NDSI) data. Daily 30 m snow cover maps were thus generated for 2000 – 2017 with an overall accuracy of 90%, using 33 validation sites distributed throughout south-central British Columbia. Snow cover trends were analyzed across stratified elevation bands and land cover types, revealing that snow cover persists under lower elevation forests for an average of 23.5 d longer than in adjacent open areas during spring. We conclude that the MODSAT-NDSI approach captures temporal and spatial advantages of freely available snow cover datasets and can be modified to suit a variety of novel investigations relating to snow cover or other spectral indices.

An Algorithm for Subpixel Snow Mapping: Extraction of a Fractional Snow-Covered Area Based on Ten-Day Composited AVHRR\/2 Data of the Qinghai-Tibet Plateau

Advanced very-high-resolution radiometer (AVHRR)/2 sensors lack the e 1.6-nm channel used by Moderate Resolution Imaging Spectroradiometer (MODIS) and AVHRR/3; consequently, the subpixel snow-mapping algorithms for MODIS and AVHRR/3 data cannot be introduced into the subpixel snow mapping for AVHRR/2 data. To formulate subpixel snow mapping for AVHRR/2 data, in this research we determined an algorithm for subpixel snow mapping using ten-day composited AVHRR/2 data from the Qinghai-Tibet Plateau and extending moderate-resolution fractional snow-covered-area data to more than 30 years.

SnowCloudHydro—A New Framework for Forecasting Streamflow in Snowy, Data-Scarce Regions

We tested the efficacy and skill of SnowCloud, a prototype web-based, cloud-computing framework for snow mapping and hydrologic modeling. SnowCloud is the overarching framework that functions within the Google Earth Engine cloud-computing environment. SnowCloudMetrics is a sub-component of SnowCloud that provides users with spatially and temporally composited snow cover information in an easy-to-use format. SnowCloudHydro is a simple spreadsheet-based model that uses Snow Cover Frequency (SCF) output from SnowCloudMetrics as a key model input. In this application, SnowCloudMetrics rapidly converts NASA’s Moderate Resolution Imaging Spectroradiometer (MODIS) daily snow cover product (MOD10A1) into a monthly snow cover frequency for a user-specified watershed area. SnowCloudHydro uses SCF and prior monthly streamflow to forecast streamflow for the subsequent month. We tested the skill of SnowCloudHydro in three snow-dominated headwaters that represent a range of precipitation/snowmelt runoff categories: the Río Elqui in Northern Chile; the John Day River, in the Northwestern United States; and the Río Aragón in Northern Spain. The skill of the SnowCloudHydro model directly corresponded to snowpack contributions to streamflow. Watersheds with proportionately more snowmelt than rain provided better results (R2 values: 0.88, 0.52, and 0.22, respectively). To test the user experience of SnowCloud, we provided the tools and tutorials in English and Spanish to water resource managers in Chile, Spain, and the United States. Participants assessed their user experience, which was generally very positive. While these initial results focus on SnowCloud, they outline methods for developing cloud-based tools that can function effectively across cultures and languages. Our approach also addresses the primary challenges of science-based computing; human resource limitations, infrastructure costs, and expensive proprietary software. These challenges are particularly problematic in countries where scientific and computational resources are underdeveloped.

Quantifying the early snowmelt event of 2015 in the Cascade Mountains, USA by developing and validating MODIS-based snowmelt timing maps

Spring snowmelt serves as the major hydrological contribution to many watersheds of the US West. Since the 1970s the conterminous western USA has seen an earlier arrival of spring snowmelt. The extremely low snowpack and early melt of 2015 in the Cascade Mountains may be a harbinger of winters to come, underscoring the interest in advancements in spring snowmelt monitoring. Target-of-opportunity and point measurements of snowmelt using meteorological stations or stream gauges are common sources of these data, however, there have been few attempts to identify snowmelt timing using remote sensing. In this study, we describe the creation of snowmelt timing maps (STMs) which identify the day of year that each pixel of a remotely sensed image transitions from “snow-covered” to “no snow” during the spring melt season, controlling for cloud coverage and ephemeral spring snow storms. Derived from the 500 m MODerate-resolution Imaging Spectroradiometer (MODIS) standard snow map, MOD10A2, this new dataset provides annual maps of snowmelt timing, with corresponding maps of cloud interference and interannual variability in snow coverage from 2001–2015. We first show that the STMs agree strongly with in-situ snow telemetry (SNOTEL) meteorological station measurements in terms of snowmelt timing. We then use the STMs to investigate the early snowmelt event of 2015 in the Cascade Mountains, USA, highlighting the protected areas of Mt. Rainier, Crater Lake, and Lassen Volcanic National Parks. In 2015 the Cascade Mountains experienced snowmelt 41 days earlier than the 2001–2015 average, with 25% of its land area melting>65 days earlier than average. The upper elevations of the Cascade Mountains experienced the greatest snowmelt anomaly. Our results are relevant to land managers and biologists as they plan adaptation strategies for mitigating the effects of climate change throughout temperate mountains.

Elimination of clouds in AVHRR/2 images of Qinghai-Tibet Plateau

ABSTRACTFor achieving the snow mapping of AVHRR/2 (i.e. the second version of Advanced Very High Resolution Radiometer) data, clouds need to be eliminated in AVHRR/2 images. However, the clouds are difficult to be identified from the AVHRR/2 image with clouds and snow together. Here, clouds were divided into high clouds, medium clouds, low clouds, and thin clouds. Based on the respective thresholds from experiments using AVHRR/2 data of Qinghai–Tibet Plateau, the four kinds of clouds were identified and eliminated from the AVHRR/2 images. According to visual validation, the results of cloud elimination were reliable and accurate.