Linking glacier annual mass balance and glacier albedo retrieved from MODIS data

Albedo is one of the variables controlling the mass balance of temperate glaciers. Multispectral imagers, such as MODerate Imaging Spectroradiometer (MODIS) on board the TERRA and AQUA satellites, provide a means to monitor glacier surface albedo. In this study, different methods to retrieve broadband glacier surface albedo from MODIS data are compared. The effect of multiple reflections due to the rugged topography and of the anisotropic reflection of snow and ice are particularly investigated. The methods are tested on the Saint Sorlin Glacier (Grandes Rousses area, French Alps). The accuracy of the retrieved albedo is estimated using both field measurements, at two automatic weather stations located on the glacier, and albedo values derived from terrestrial photographs. For summers 2008 and 2009, the root mean square deviation (RMSD) between field measurements and the broadband albedo retrieved from MODIS data at 250 m spatial resolution was found to be 0.052 or about 10% relative error. The RMSD estimated for the MOD10 daily albedo product is about three times higher. One decade (2000–2009) of MODIS data were then processed to create a time series of albedo maps of Saint Sorlin Glacier during the ablation season. The annual mass balance of Saint Sorlin Glacier was compared with the minimum albedo value (average over the whole glacier surface) observed with MODIS during the ablation season. A strong linear correlation exists between the two variables. Furthermore, the date when the average albedo of the whole glacier reaches a minimum closely corresponds to the period when the snow line is located at its highest elevation, thus when the snow line is a good indicator of the glacier equilibrium line. This indicates that this strong correlation results from the fact that the minimal average albedo values of the glacier contains considerable information regarding the relative share of areal surfaces between the ablation zone (i.e. ice with generally low albedo values) and the accumulation zone (i.e. snow with a relatively high albedo). As a consequence, the monitoring of the glacier surface albedo using MODIS data can provide a useful means to evaluate the interannual variability of the glacier mass balance. Finally, the albedo in the ablation area of Saint Sorlin Glacier does not exhibit any decreasing trend over the study period, contrasting with the results obtained on Morteratsch Glacier in the Swiss Alps

Effects of atmospheric, topographic, and BRDF correction on imaging spectroscopy-derived data products

Surface reflectance is an important data product in imaging spectroscopy for obtaining surface information. The complex retrieval of surface reflectance, however, critically relies on accurate knowledge of atmospheric absorption and scattering, and the compensation of these effects. Furthermore, illumination and observation geometry in combination with surface reflectance anisotropy determine dynamics in retrieved surface reflectance not related to surface absorption properties. To the best of authors’ knowledge, no comprehensive assessment of the impact of atmospheric, topographic, and anisotropy effects on derived surface information is available so far.This study systematically evaluates the impact of these effects on reflectance, albedo, and vegetation products. Using three well-established processing schemes (ATCOR F., ATCOR R., and BREFCOR), high-resolution APEX imaging spectroscopy data, covering a large gradient of illumination and observation angles, are brought to several processing states, varyingly affected by mentioned effects. Pixel-wise differences of surface reflectance, albedo, and spectral indices of neighboring flight lines are quantitatively analyzed in their respective overlapping area. We found that compensation of atmospheric effects reveals actual anisotropy-related dynamics in surface reflectance and derived albedo, related to an increase in pixel-wise relative reflectance and albedo differences of more than 40%. Subsequent anisotropy compensation allows us to successfully reduce apparent relative reflectance and albedo differences by up to 20%. In contrast, spectral indices are less affected by atmospheric and anisotropy effects, showing relative differences of 3% to 10% in overlapping regions of flight lines.We recommend to base decisions on the use of appropriate processing schemes on individual use cases considering envisioned data products

LANDSAT-D investigations in snow hydrology

Work undertaken during the contract and its results are described. Many of the results from this investigation are available in journal or conference proceedings literature - published, accepted for publication, or submitted for publication. For these the reference and the abstract are given. Those results that have not yet been submitted separately for publication are described in detail. Accomplishments during the contract period are summarized as follows: (1) analysis of the snow reflectance characteristics of the LANDSAT Thematic Mapper, including spectral suitability, dynamic range, and spectral resolution; (2) development of a variety of atmospheric models for use with LANDSAT Thematic Mapper data. These include a simple but fast two-stream approximation for inhomogeneous atmospheres over irregular surfaces, and a doubling model for calculation of the angular distribution of spectral radiance at any level in an plane-parallel atmosphere; (3) incorporation of digital elevation data into the atmospheric models and into the analysis of the satellite data; and (4) textural analysis of the spatial distribution of snow cover

Active microwave remote sensing of earth/land, chapter 2

Geoscience applications of active microwave remote sensing systems are examined. Major application areas for the system include: (1) exploration of petroleum, mineral, and ground water resources, (2) mapping surface and structural features, (3) terrain analysis, both morphometric and genetic, (4) application in civil works, and (5) application in the areas of earthquake prediction and crustal movements. Although the success of radar surveys has not been widely publicized, they have been used as a prime reconnaissance data base for mineral exploration and land-use evaluation in areas where photography cannot be obtained

Anisotropy parameterization development and evaluation for glacier surface albedo retrieval from satellite observations.

Glacier albedo determines the net shortwave radiation absorbed at the glacier surface and plays a crucial role in glacier energy and mass balance. Remote sensing techniques are efficient means to retrieve glacier surface albedo over large and inaccessible areas and to study its variability. However, corrections of anisotropic reflectance of glacier surface have been established for specific shortwave bands only, such as Landsat 5 Thematic Mapper (L5/TM) band 2 and band 4, which is a major limitation of current retrievals of glacier broadband albedo. In this study, we calibrated and evaluated four anisotropy correction models for glacier snow and ice, applicable to visible, near-infrared and shortwave-infrared wavelengths using airborne datasets of Bidirectional Reflectance Distribution Function (BRDF). We then tested the ability of the best-performing anisotropy correction model, referred to from here on as the ‘updated model’, to retrieve albedo from L5/TM, Landsat 8 Operational Land Imager (L8/OLI) and Moderate Resolution Imaging Spectroradiometer (MODIS) imagery, and evaluated these results with field measurements collected on eight glaciers around the world. Our results show that the updated model: (1) can accurately estimate anisotropic factors of reflectance for snow and ice surfaces; (2) generally performs better than prior approaches for L8/OLI albedo retrieval but is not appropriate for L5/TM; (3) generally retrieves MODIS albedo better than the MODIS standard albedo product (MCD43A3) in both absolute values and glacier albedo temporal evolution, i.e., exhibiting both fewer gaps and better agreement with field observations. As the updated model enables anisotropy correction of a maximum of 10 multispectral bands and is implemented in Google Earth Engine (GEE), it is promising for observing and analyzing glacier albedo at large spatial scales

Recent advances in remote sensing of seasonal snow

Remote sensing offers local, regional and global observations of seasonal snow, providing key information on snowpack processes. This brief review highlights advancements in instrumentation and analysis techniques that have been developed over the past decade. Areas of advancement include improved algorithms for mapping snow-cover extent, snow albedo, snow grain size, snow water equivalent, melt detection and snow depth, as well as new uses of instruments such as multiangular spectroradiometers, scatterometry and lidar. Limitations and synergies of the instruments and techniques are discussed, and remaining challenges such as multisensor mapping, scaling issues, vegetation correction and data assimilation are identified

드론을 활용한 위성 지표반사도 산출물 공간 패턴 분석

학위논문(석사) -- 서울대학교대학원 : 농업생명과학대학 생태조경·지역시스템공학부(생태조경학), 2021.8. 조대솔.High-resolution satellites are assigned to monitor land surface in detail. The reliable surface reflectance (SR) is the fundamental in terrestrial ecosystem modeling so the temporal and spatial validation is essential. Usually based on multiple ground control points (GCPs), field spectroscopy guarantees the temporal continuity. Due to limited sampling, however, it hardly illustrates the spatial pattern. As a map, the pixelwise spatial variability of SR products is not well-documented. In this study, we introduced drone-based hyperspectral image (HSI) as a reference and compared the map with Sentinel 2 and Landsat 8 SR products on a heterogeneous rice paddy landscape. First, HSI was validated by field spectroscopy and swath overlapping, which assured qualitative radiometric accuracy within the viewing geometry. Second, HSI was matched to the satellite SRs. It involves spectral and spatial aggregation, co-registration and nadir bidirectional reflectance distribution function (BRDF)-adjusted reflectance (NBAR) conversion. Then, we 1) quantified the spatial variability of the satellite SRs and the vegetation indices (VIs) including NDVI and NIRv by APU matrix, 2) qualified them pixelwise by theoretical error budget and 3) examined the improvement by BRDF normalization. Sentinel 2 SR exhibits overall good agreement with drone HSI while the two NIRs are biased up to 10%. Despite the bias in NIR, the NDVI shows a good match on vegetated areas and the NIRv only displays the discrepancy on built-in areas. Landsat 8 SR was biased over the VIS bands (-9 ~ -7.6%). BRDF normalization just contributed to a minor improvement. Our results demonstrate the potential of drone HSI to replace in-situ observation and evaluate SR or atmospheric correction algorithms over the flat terrain. Future researches should replicate the results over the complex terrain and canopy structure (i.e. forest).원격탐사에서 지표 반사도(SR)는 지표정보를 비파괴적이고 즉각적인 방법으로 전달해주는 매개체 역할을 한다. 신뢰할 수 있는 SR은 육상 생태계 모델링의 기본이고, 이에 따라 SR의 시공간적 검증이 요구된다. 일반적으로 SR은 여러 지상 기준점(GCP)을 기반으로 하는 현장 분광법을 통해서 시간적 연속성이 보장된다. 그러나 현장 분광법은 제한적인 샘플링으로 공간 패턴을 거의 보여주지 않아, 위성 SR의 픽셀 별 공간 변동성은 잘 분석되지 않았다. 본 연구에서는 드론 기반의 초분광 영상(HSI)을 참고자료로 도입하여, 이를 이질적인 논 경관에서 Sentinel 2 및 Landsat 8 SR과 비교하였다. 우선, 드론 HSI는 현장 분광법 및 경로 중첩을 통해서 관측각도 범위 내에서 정성적인 방사 측정을 보장한다고 검증되었다. 이후, 드론 HSI는 위성 SR의 분광반응특성, 공간해상도 및 좌표계를 기준으로 맞춰졌고, 관측 기하를 통일하기 위해서 드론 HIS와 위성 SR은 각각 양방향반사율분포함수 (BRDF) 정규화 반사도 (NBAR)로 변환되었다. 마지막으로, 1) APU 행렬으로 위성 SR과 NDVI, NIRv를 포함하는 식생지수(VI)의 공간변동성을 정량화 했고, 2) 대기보정의 이론적 오차를 기준으로 SR과 VI를 픽셀별로 평가했고, 3) BRDF 정규화를 통한 개선 사항을 검토했다. Sentinel 2 SR은 드론 HSI와 전반적으로 좋은 일치를 보이나, 두 NIR 채널은 최대 10% 편향되었다. NIR의 편향은 식생지수에서 토지 피복에 따라 다른 영향을 미쳤다. NDVI는 식생에서는 낮은 편향을 보여줬고, NIRv는 도시시설물 영역에서만 높은 편향을 보였다. Landsat 8 SR은 VIS 채널에 대해 편향되었다 (-9 ~ -7.6%). BRDF 정규화는 위성 SR의 품질을 개선했지만, 그 영향은 부수적이었다. 본 연구에서는 평탄한 지형에서 드론 HSI가 현장 관측을 대체할 수 있고, 따라서 위성 SR이나 대기보정 알고리즘을 평가하는데 활용될 수 있다는 것을 보였다. 향후 연구에서는 산림으로 대상지를 확대하여, 지형과 캐노피 구조가 드론 HSI 및 위성 SR에 미치는 영향을 분석할 필요가 있다.Chapter 1. Introduction 1 1.1 Background 1 Chapter 2. Method 3 2.1 Study Site 3 2.2 Drone campaign 4 2.3 Data processing 4 2.3.1 Sensor calibration 5 2.3.2 Bidirectional reflectance factor (BRF) calculation 7 2.3.3 BRDF correction 7 2.3.4 Orthorectification 8 2.3.5 Spatial Aggregation 9 2.3.6 Co-registration 10 2.4 Satellite dataset 10 2.4.2 Landsat 8 12 Chapter 3. Result and Discussion 12 3.1 Drone BRF map quality assessment 12 3.1.1 Radiometric accuracy 12 3.1.2 BRDF effect 15 3.2 Spatial variability in satellite surface reflectance product 16 3.2.1 Sentinel 2B (10m) 17 3.2.2 Sentinel 2B (20m) 22 3.2.3 Landsat 8 26 Chapter 4. Conclusion 28 Supplemental Materials 30 Bibliography 34 Abstract in Korean 43석