Use of landsat ETM+ SLC-off segment-based gap-filled imagery for crop type mapping

Failure of the Scan Line Corrector (SLC) on the Landsat ETM+ sensor has had a major impact on many applications that rely on continuous medium resolution imagery to meet their objectives. The United States Department of Agriculture (USDA) Cropland Data Layer (CDL) program uses Landsat imagery as the primary source of data to produce crop-specific maps for 20 states in the USA. A new method has been developed to fill the image gaps resulting from the SLC failure to support the needs of Landsat users who require coincident spectral data, such as for crop type mapping and monitoring. We tested the new gap-filled method for a CDL crop type mapping project in eastern Nebraska. Scan line gaps were simulated on two Landsat 5 images (spring and late summer 2003) and then gap-filled using landscape boundary models, or segment models, that were derived from 1992 and 2002 Landsat images (used in the gap-fill process). Various date combinations of original and gap-filled images were used to derive crop maps using a supervised classification process. Overall kappa values were slightly higher for crop maps derived from SLC-off gap-filled images compared to crop maps derived from the original imagery (0.3–1.3% higher). Although the age of the segment model used to derive the SLC-off gap-filled product did not negatively impact the overall agreement, differences in individual cover type agreement did increase (?0.8%–1.6% using the 2002 segment model to ?5.0–5.1% using the 1992 segment model). Classification agreement also decreased for most of the classes as the size of the segment used in the gap-fill process increased.     

MODIS植被指数的美国玉米单产遥感估测

针对中国开展的国外农作物产量遥感估测大多依靠中低分辨率耕地信息、省级(州级)或国家级作物产量统计数据的现状,本文以美国玉米为例,探讨利用多年中高分辨率作物分布信息、时序遥感植被指数和县级作物产量统计数据开展国外重点地区作物单产遥感估测技术研究,以期进一步提高中国对国外农作物产量监测精度和精细化水平。首先,利用美国农业部国家农业统计局(NASS/USDA)生产的作物分布数据(CDL)获得多个年份玉米空间分布图,并对相应年份250 m分辨率16天合成的MODIS-NDVI时序数据进行掩膜处理,统计获得每年各县域内玉米主要生育期NDVI均值;其次,以各州为估产区,以多年县级玉米统计单产和县域内玉米主要生育期NDVI均值为基础,建立各州玉米主要生育期NDVI与玉米单产间关系模型;然后,通过主要生育期玉米单产和玉米植被指数间拟合程度,筛选确定各州玉米最佳估产期和最佳估产模型。最终,利用最佳估产模型实现美国各州玉米单产估测和全国玉米单产推算。其中,建模数据覆盖时间为2007年—2010年,验证数据为2011年。结果表明,应用最佳估产模型的2011年美国各州玉米单产估测相对误差在-4.16%—4.92%,均方根误差在148.75—820.93 kg/ha,各州估测结果计算获得全国玉米单产的相对误差仅为2.12%,均方根误差为285.57 kg/ha。可见,本研究的作物单产遥感估测技术方法具有一定可行性,可准确估测全球重点地区作物单产信息。     

Relative performance of landsat MSS,TM and IRS LISS-I data for acreage estimation of Rabi sorghum

Acreage estimation of Rabi sorghum crop in Ahmadnagar, Pune and Solapur districts of central Maharashtra has been attempted using synchronously acquired Landsat MSS and TM data of 1987–88 season and IRS LISS-I data of 1988–89 season; in conjuction with near-synchronous ground truth data. The remote-sensing-based acreage estimations for the districts were compared with the respective estimates by Bureau of Economics and Statistics (BES). As the acreages were underestimated with the classification of standard four-band MSS data, the atmospheric correction of fourband MSS data and normalised differencing (ND) of the atmospheric-corrected MSS data were attempted. The main observations are: (1) the use of Landsat MSS data results in underestimation of sorghum acreage in comparison with BES estimation, (2) the atmospheric correction and ND transformation of MSS data are necessary for bringing acreage estimates in agreement with BES estimates, (3) Mid-IR data in band 1.55 to 1.75 μm are useful in improving the separability of land-use classes, and (4) remote sensing data with radiometric sensitivity comparable to LISS-I or Landsat TM and Signal-to-Noise ratios comparable to LISS-I data are suitable for accurate acreage estimation of sorghum.     

Crop classification and acreage estimation in North Korea using phenology features

In North Korea, reliable and timely information on crop acreage and spatial distribution is hard to obtain. In this study, we developed a fast and robust method to estimate crop acreage in North Korea using time-series normalized difference vegetation index (NDVI) derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) data. We proposed a method to identify crop type based on NDVI phenology features using data collected in other areas with similar agri-environmental conditions to mitigate the shortage of ground truth data. Eventually the classification map (MODIScrop) was assessed using the Food and Agriculture Organization (FAO) statistical data and high-resolution crop classification maps derived from one Landsat scene (LScrop). The Pareto boundary method was used to assess the accuracy and crop distribution of the MODIScrop maps. Results showed that acreage derived from the MODIScrop maps was generally consistent with that reported in the FAO data (a relative error <4.1% for rice and <6.1% for maize, and <9.0% for soybean except for in 2004, 2008, and 2009) and the maps derived from the LScrop (a relative error about 5% in 2013, and 7% in 2008 and 2014). The classification accuracy reached 74.4%, 69.8%, and 73.1% of the areas covered by the Landsat images in 2008, 2013, and 2014, respectively. This indicates that features derived from NDVI profiles were able to characterize major crops, and the approaches developed in this study are feasible for crop mapping and acreage estimation in regions with limited ground truth data.     

Measuring land-use and land-cover change using the U.S. department of agriculture’s cropland data layer: Cautions and recommendations

Monitoring agricultural land is important for understanding and managing food production, environmental conservation efforts, and climate change. The United States Department of Agriculture's Cropland Data Layer (CDL), an annual satellite imagery-derived land cover map, has been increasingly used for this application since complete coverage of the conterminous United States became available in 2008. However, the CDL is designed and produced with the intent of mapping annual land cover rather than tracking changes over time, and as a result certain precautions are needed in multi-year change analyses to minimize error and misapplication. We highlight scenarios that require special considerations, suggest solutions to key challenges, and propose a set of recommended good practices and general guidelines for CDL-based land change estimation. We also characterize a problematic issue of crop area underestimation bias within the CDL that needs to be accounted for and corrected when calculating changes to crop and cropland areas. When used appropriately and in conjunction with related information, the CDL is a valuable and effective tool for detecting diverse trends in agriculture. By explicitly discussing the methods and techniques for post-classification measurement of land-cover and land-use change using the CDL, we aim to further stimulate the discourse and continued development of suitable methodologies. Recommendations generated here are intended specifically for the CDL but may be broadly applicable to additional remotely-sensed land cover datasets including the National Land Cover Database (NLCD), Moderate Resolution Imaging Spectroradiometer (MODIS)-based land cover products, and other regional, national, and global land cover classification maps.     

Evaluation of global land cover maps for cropland area estimation in the conterminous United States

Global land cover data could provide continuously updated cropland acreage and distribution information, which is essential to a wide range of applications over large geographical regions. Cropland area estimates were evaluated in the conterminous USA from four recent global land cover products: MODIS land cover (MODISLC) at 500-m resolution in 2010, GlobCover at 300-m resolution in 2009, FROM-GLC and FROM-GLC-agg at 30-m resolution based on Landsat imagery circa 2010 against the US Department of Agriculture survey data. Ratio estimators derived from the 30-m resolution Cropland Data Layer were applied to MODIS and GlobCover land cover products, which greatly improved the estimation accuracy of MODISLC by enhancing the correlation and decreasing mean deviation (MDev) and RMSE, but were less effective on GlobCover product. We found that, in the USA, the CDL adjusted MODISLC was more suitable for applications that concern about the aggregated county cropland acreage, while FROM-GLC-agg gave the least deviation from the survey at the state level. Correlation between land cover map estimates and survey estimates is significant, but stronger at the state level than at the county level. In regions where most mismatches happen at the county level, MODIS tends to underestimate, whereas MERIS and Landsat images incline to overestimate. Those uncertainties should be taken into consideration in relevant applications. Excluding interannual and seasonal effects, R² of the FROM-GLC regression model increased from 0.1 to 0.4, and the slope is much closer to one. Our analysis shows that images acquired in growing season are most suitable for Landsat-based cropland mapping in the conterminous USA.     

Kharif rice acreage estimation in Orissa using NOAA-AVHRR data

Rice-acreage estimation of Orissa state was carried out using single-date NOAA-AVHRR data. Selection of optimum date of data acquisition for this purpose was studied using data of six acquisition dates viz. October, 3, 12, 21, 29, November 7 and 26, 1989. Comparative performance of MXL classification of two NOAA bands (Band-1: 0.58–0.68 μm and Band-2: 0.73–1.10 μm) and Normalised Difference Vegetation Index (NDVI) derived from these two-band data was examined. Acreage thus estimated was compared against Bureau of Economics and Statistics (BES) estimate of the same year. The acreage estimation obtained by two band classification was closer to BES estimate than that based on NDVI. Data acquired in the month of October have given better estimate for state level rice acreage than that acquired in the month of November.     

The use of MODIS data to derive acreage estimations for larger fields: A case study in the south-western Rostov region of Russia

Recent developments in remote sensing technology, in particular improved spatial and temporal resolution, open new possibilities for estimating crop acreage over larger areas. Remotely sensed data allow in some cases the estimation of crop acreage statistics independently of sub-national survey statistics, which are sometimes biased and incomplete. This work focuses on the use of MODIS data acquired in 2001/2002 over the Rostov Oblast in Russia, by the Azov Sea. The region is characterised by large agricultural fields of around 75 ha on average. This paper presents a methodology to estimate crop acreage using the MODIS 16-day composite NDVI product. Particular emphasis is placed on a good quality crop mask and a good quality validation dataset. In order to have a second dataset which can be used for cross-checking the MODIS classification a Landsat ETM time series for four different dates in the season of 2002 was acquired and classified. We attempted to distinguish five different crop types and achieved satisfactory and good results for winter crops. Three hundred and sixty fields were identified to be suitable for the training and validation of the MODIS classification using a maximum likelihood classification. A novel method based on a pure pixel field sampling is introduced. This novel method is compared with the traditional hard classification of mixed pixels and was found to be superior.     

Accuracy assessment in cotton acreage estimation using Indian remote sensing satellite data

The accuracy of cotton crop classification using satellite data has been assessed with respect to a detailed land cover map prepared by field survey. The effect of spatial resolution on classification accuracy was studied using LISS-I (spatial resolution 72.6 m) and LISS-II data (spatial resolution 36.25 m) of the Indian remote sensing satellite IRS-1B. The performances of the maximum likelihood and the minimum distance to mean as classifiers have also been assessed. LISS-II data have been found to give a higher classification accuracy. The estimate of cotton acreage using LISS-II data was closer to that obtained from the base map. The maximum likelihood classifier (MXL) and the minimum distance to mean (MDM) classifier performed equally well.     

全球地图第一版中国土地覆盖产品的生成

介绍了全球地图第一版中国土地覆盖分类产品的遥感数据源、数据处理方法及结果,详细介绍了非监督聚类和决策树分类相结合的区域级土地覆盖分类方法及处理步骤,并以2003年全年的MODIS/TERRA16d合成产品MOD43B4为主要遥感数据源,采用先进行影像聚类,然后以聚类结果为处理单元进行多时相多特征决策树分类,得到了覆盖全国、约1km分辨率的土地覆盖分类结果。采用室内精度评定和外业验证两种方法对分类结果进行了检验,生成的中国土地覆盖分类产品的精度为89.14%。     

Multiple forecasts of kharif rice in orissa state-four year experience of fasal pilot study

Considering the requirement of multiple pre-harvest crop forecasts, the concept of Forecasting Agricultural output using Space, Agrometeorology and Land based observations (FASAL) has been formulated. Development of procedure and demonstration of this technique for four in-season forecasts for kharif rice has been carried out as a pilot study in Orissa State since 1998. As the availability of cloud-free optical remote sensing data during kharif season is very poor for Orissa state, multi-date RADARSAT SCANSAR data were used for acreage estimation of kharif rice. Meteorological models have been developed for early assessment of acreage and prediction of yield at mid and late crop growth season. Four in-season forecasts were made during four kharif seasons (1998-2001); the first forecast of zone level rice acreage at the beginning of kharif crop season using meteorological models, second forecast of district level acreage at mid growth season using two-date RADARSAT SCANSAR data and yield using meteorological models, third forecast at late growth season of district level acreage using three-date RADARSAT SCANSAR data and yield using meteorological models and revised forecast incorporating field observations at maturity. The results of multiple forecasts have shown rice acreage estimation and yield prediction with deviation up to 14 and 11 per cent respectively. This study has demonstrated the potential of FASAL concept to provide inseason multiple forecasts using data of remote sensing, meteorology and land based observations.     

PhenoCrop: An integrated satellite-based framework to estimate physiological growth stages of corn and soybeans

The accurate and timely estimates of crop physiological growth stages are essential for efficient crop management and precise modeling of agricultural systems. Satellite remote sensing has been widely used to retrieve vegetation phenology metrics at local to global scales. However, most of these phenology metrics (e.g., green-up) are different from crop growth stages (e.g., emergence) used in crop management and modeling. As such, an integrated framework referred to as PhenoCrop was developed to: 1) establish a connection between remote sensing-derived phenology metrics and key crop growth stages based on Wang and Engle plant phenology model and 2) use fused MODIS-Landsat 30 m 8-day reflectance data generated using Kalman Filter-based data fusion technique to produce onset dates of key growth stages of corn (Zea mays L.) and soybeans (Glycine max L.) at 30 m spatial resolution. In this paper, we described the PhenoCrop framework, and tested its performance for the State of Nebraska for 2012–2016 by comparison to observations of estimated key growth stages at four experimental sites, and state-level statistical data from Crop Progress Reports (CPRs) published by the United States Department of Agriculture’s (USDA) National Agricultural Statistical Services (NASS). In addition, to evaluate the suitability of using coarse or high spatial resolution satellite imagery, fused MODIS-Landsat-based estimates were compared with those produced using EOS MODIS 250 m (MOD9Q1) reflectance data.The PhenoCrop estimates captured the typical spatial trends of gradual delay in the progression of the growing season from southeast to northwest Nebraska. Also inter-annual differences due to factors such as weather fluctuations and change in management strategies (e.g., early season in 2012) were evident in the estimates. Validation results revealed that average root mean square error (RMSE) of the state-level estimates of corn and soybean growth stages ranged from 1.10 to 4.20 days and from 3.81 to 7.89 days, respectively, while pixel level estimates had a RMSE ranging from 3.72 to 8.51 days for corn and 4.76–9.51 days for soybean growth stages. Although MODIS 250 m based estimates showed similar general spatial patterns observed in the fused MODIS-Landsat based estimates, the accuracy and ability to capture field scale variations was improved with fused MODIS-Landsat data. Overall, results showed the ability of PhenoCrop framework to provide reliable estimates of crop growth stages that can be highly useful in crop modeling and crop management during the growing season.     

Results on potential use of simulated IRS-1C WiFS data for crop monitoring

Some of the basic requirements for cropping system analysis are updated information on crops grown, their phenological behaviour, method and duration of establishment and harvest, inter and intra crop variability, sequential cropping patterns. The next generation Indian Remote Sensing Satellite with high repeat cycle opens new possibility of crop surveys to derive such information. In this study, an attempt has been made to analyse cropping system at district level using simulated IRS-1C Wide Field Sensor (WiFS) data. Data acquired for nineteen dates during 1992–93 season for Bardhaman district, West Bengal has been used. It was feasible to derive accurate information on cropping pattern, crop rotation, crop duration, progress of harvest, crop growth profiles and annual crop acreage using multidate data. It was observed that even a seven to eight day interval of data acquisition during critical growth periods significantly affected classification and identification accuracy.     

极化SAR参数优化与光学波谱相结合的面向对象土地覆盖分类

提出一种基于极化参数优化的面向对象分类方法。该方法结合光学和SAR数据,有效提高了对地物的识别能力。本文方法的关键在于:在■分解中,使用光学影像指导SAR影像选择同质点,使其更精确地估计极化参数并结合光学波谱信息作为输入特征;使用面向对象的分类方法,仅将光学影像作为分割输入,避免SAR噪声引起的分割错误。以美国Bakersfield地区的Sentinel-1/2数据为例,确定7种地物类型,对比分析不同输入与不同分类器对分类结果的影响。研究表明,优化输入参数在纹理丰富区域能够有效提高分类精度;面向对象的分类结果更加稳定并较好地维持地表几何特征;改进分类方法较传统分类方法总体精度提高了近10%,达到92.6%。     

Estimating crop net primary production using national inventory data and MODIS-derived parameters

National estimates of spatially-resolved cropland net primary production (NPP) are needed for diagnostic and prognostic modeling of carbon sources, sinks, and net carbon flux between land and atmosphere. Cropland NPP estimates that correspond with existing cropland cover maps are needed to drive biogeochemical models at the local scale as well as national and continental scales. Existing satellite-based NPP products tend to underestimate NPP on croplands. An Agricultural Inventory-based Light Use Efficiency (AgI-LUE) framework was developed to estimate individual crop biophysical parameters for use in estimating crop-specific NPP over large multi-state regions. The method is documented here and evaluated for corn (Zea mays L.) and soybean (Glycine max L. Merr.) in Iowa and Illinois in 2006 and 2007. The method includes a crop-specific Enhanced Vegetation Index (EVI), shortwave radiation data estimated using the Mountain Climate Simulator (MTCLIM) algorithm, and crop-specific LUE per county. The combined aforementioned variables were used to generate spatially-resolved, crop-specific NPP that corresponds to the Cropland Data Layer (CDL) land cover product. Results from the modeling framework captured the spatial NPP gradient across croplands of Iowa and Illinois, and also represented the difference in NPP between years 2006 and 2007. Average corn and soybean NPP from AgI-LUE was 917 g C m⁻² yr⁻¹ and 409 g C m⁻² yr⁻¹, respectively. This was 2.4 and 1.1 times higher, respectively, for corn and soybean compared to the MOD17A3 NPP product. Site comparisons with flux tower data show AgI-LUE NPP in close agreement with tower-derived NPP, lower than inventory-based NPP, and higher than MOD17A3 NPP. The combination of new inputs and improved datasets enabled the development of spatially explicit and reliable NPP estimates for individual crops over large regional extents.     

国家土地覆盖数据的生产与建库

土地覆盖变化的研究和分析,能够揭示自然与人文过程交叉最密切的一些问题,为此各国均开展了本国的国家级土地覆盖数据的生产。本文在参考国外土地覆盖数据生产相关经验的基础上,结合我国的国情,基于陆地卫星影像,研制了30m土地覆盖分类系统,采用非监督分类、监督分类与人工解译相结合的方法,生产了全国范围的土地覆盖数据;并基于ORACLE数据库;建立了国家土地覆盖数据库。     

Land cover monitoring in Asia by NOAA GVI data

NOAA/AVHRR Global Vegetation Index (GVI) data of Asia in 1983 and 1987 were used to evaluate their usefulness for global land cover monitoring. Color composite images of monthly GVI data and color composite images of principal components from 12 successive monthly GVI data were found to be useful for visual interpretation of seasonal vegetation dynamics. The results of cluster analysis applied to monthly GVI data for a one‐year period, indicate that unsupervised classification method is useful for global or continental land cover classification without ground truth. In order to detect land cover changes, the difference between the 1983 and 1987 12‐month GVI data was calculated. The results show that it is difficult to detect land cover changes due to cloud contamination in monthly GVI data and poor registration of GVI products.     

Crop area monitoring within an advanced agricultural information system

This paper describes a framework for an image processing procedure for operational agricultural crop area estimation. This operational framework has been conceived within the development of an Advanced Agricultural Information System (AAIS) for the “Regione del Veneto “ (RdV ‐ Veneto Region) in northeastern Italy. The objective of this program is to develop the ability to generating timely and accurate area estimates and production information for four major agricultural crops: soybeans, sugar beets, corn, and small grains. AAIS uses state of the art methods in remote sensing and geographic information systems (GIS) technology and integrates a variety of data types including satellite imagery. This paper describes the methodology developed for image and ancillary data processing for the production of crop area statistics. Using a combination of standard unsupervised classification and GIS operations that incorporate knowledge about the agricultural system, a “sequential masking” classification procedure was derived. This sequential masking procedure yielded crop classification accuracies that at the study site level range between 76% and 99% depending on the crop under study. We believe that classification accuracies will improve with full system implementation, along with the incorporation of new and/or improved thematic information and operational experience using AAIS‐based estimation.     

Methods for improving accuracy and extending results beyond periods covered by traditional ground-truth in remote sensing classification of a complex landscape

Successful development of approaches to quantify impacts of diverse landuse and associated agricultural management practices on ecosystem services is frequently limited by lack of historical and contemporary landuse data. We hypothesized that ground truth data from one year could be used to extrapolate previous or future landuse in a complex landscape where cropping systems do not generally change greatly from year to year because the majority of crops are established perennials or the same annual crops grown on the same fields over multiple years. Prior to testing this hypothesis, it was first necessary to classify 57 major landuses in the Willamette Valley of western Oregon from 2005 to 2011 using normal same year ground-truth, elaborating on previously published work and traditional sources such as Cropland Data Layers (CDL) to more fully include minor crops grown in the region. Available remote sensing data included Landsat, MODIS 16-day composites, and National Aerial Imagery Program (NAIP) imagery, all of which were resampled to a common 30 m resolution. The frequent presence of clouds and Landsat7 scan line gaps forced us to conduct of series of separate classifications in each year, which were then merged by choosing whichever classification used the highest number of cloud- and gap-free bands at any given pixel. Procedures adopted to improve accuracy beyond that achieved by maximum likelihood pixel classification included majority-rule reclassification of pixels within 91,442 Common Land Unit (CLU) polygons, smoothing and aggregation of areas outside the CLU polygons, and majority-rule reclassification over time of forest and urban development areas. Final classifications in all seven years separated annually disturbed agriculture, established perennial crops, forest, and urban development from each other at 90 to 95% overall 4-class validation accuracy. In the most successful use of subsequent year ground-truth data to classify prior year landuse, an overall 57-class accuracy of 75% was achieved despite the omission of 10 entire classes, most of which were annually disturbed or perennial crops grown on very few fields. Synthetic ground-truth data for the 2004 harvest year based on the most common landuse classes over the following 7 years classified 49 of 57 categories at an overall accuracy of 96% in a final version that included CLU polygon majority rule, default smoothing and aggregation, and forcing of urban development and forest from multi-year majority-rule.     

Estimation of Snow Cover Area,Snow Physical Properties and Glacier Classification in Parts of Western Himalayas Using C-Band SAR Data

Snow physical properties, snow cover and glacier facies are important parameters which are used to quantify snowpack characteristics, glacier mass balance and seasonal snow and glacier melt. This study has been done using C-band synthetic aperture radar (SAR) data of Indian radar imaging satellite, radar imaging satellite-1 (RISAT)-1, to estimate the seasonal snow cover and retrieve snow physical properties (snow wetness and snow density), and glacier radar zones or facies classification in parts of North West Himalaya (NWH), India. Additional SAR data used are of Radarsat-2 (RS-2) satellite, which was used for glacier facies classification of Smudra Tapu glacier in Himachal Pradesh. RISAT-1 based snow cover area (SCA) mapping, snow wetness and snow density retrieval and glacier facies classification have been done for the first time in NWH region. SAR-based inversion models were used for finding out wet and dry snow dielectric constant, dry and wet SCA, snow wetness and snow density. RISAT-1 medium resolution scan-SAR mode (MRS) in HV polarization was used for first time in NWH for deriving time series of SCA maps in Beas and Bhagirathi river basins for years 2013–2014. The SAR-based inversion models were implemented separately for RISAT-1 quad pol. FRS2, for wet snow and dry snow permittivity retrieval. Masks for layover and shadow were considered in estimating final snow parameters. The overall accuracy in terms of R² value comes out to be 0.74 for snow wetness and 0.72 for snow density based on the limited ground truth data for subset area of Manali sub-basin of Beas River up to Manali for winter of 2014. Accuracy for SCA was estimated to be 95 % when compared with optical remote sensing based SCA maps with error of ±10 %. The time series data of RISAT-1 MRS and hybrid data in RH/RV mode based decompositions were also used for glacier radar zones classification for Gangotri and Samudra Tapu glaciers. The various glaciers radar zones or facies such as debris covered glacier ice, clean or bare glacier ice radar zone, percolation/refreeze radar zone and wet snow, ice wall etc., were identified. The accuracy of classified maps was estimated using ground truth data collected during 2013 and 2014 glacier field work to Samudra Tapu and Gangotri glaciers and overall accuracy was found to be in range of 82–90 %. This information of various glacier radar zones can be utilized in marking firn line of glaciers, which can be helpful for glacier mass balance studies.