Estimation of Land Surface Albedo from GCOM-C/SGLI Surface Reflectance

XXIV ISPRS Congress “Imaging today, foreseeing tomorrow, ” Commission III2021 edition, 5–9 July 2021This paper examines algorithms for estimating terrestrial albedo from the products of the Global Change Observation Mission – Climate (GCOM-C)/Second-generation Global Imager (SGLI), which was launched in December 2017 by the Japan Aerospace Exploration Agency. We selected two algorithms: one based on a bidirectional reflectance distribution function (BRDF) model and one based on multi-regression models. The former determines kernel-driven BRDF model parameters from multiple sets of reflectance and estimates the land surface albedo from those parameters. The latter estimates the land surface albedo from a single set of reflectance with multi-regression models. The multi-regression models are derived for an arbitrary geometry from datasets of simulated albedo and multi-angular reflectance. In experiments using in situ multi-temporal data for barren land, deciduous broadleaf forests, and paddy fields, the albedos estimated by the BRDF-based and multi-regression-based algorithms achieve reasonable root-mean-square errors. However, the latter algorithm requires information about the land cover of the pixel of interest, and the variance of its estimated albedo is sensitive to the observation geometry. We therefore conclude that the BRDF-based algorithm is more robust and can be applied to SGLI operational albedo products for various applications, including climate-change research

Data-Driven Artificial Intelligence for Calibration of Hyperspectral Big Data

Near-earth hyperspectral big data present both huge opportunities and challenges for spurring developments in agriculture and high-throughput plant phenotyping and breeding. In this article, we present data-driven approaches to address the calibration challenges for utilizing near-earth hyperspectral data for agriculture. A data-driven, fully automated calibration workflow that includes a suite of robust algorithms for radiometric calibration, bidirectional reflectance distribution function (BRDF) correction and reflectance normalization, soil and shadow masking, and image quality assessments was developed. An empirical method that utilizes predetermined models between camera photon counts (digital numbers) and downwelling irradiance measurements for each spectral band was established to perform radiometric calibration. A kernel-driven semiempirical BRDF correction method based on the Ross Thick-Li Sparse (RTLS) model was used to normalize the data for both changes in solar elevation and sensor view angle differences attributed to pixel location within the field of view. Following rigorous radiometric and BRDF corrections, novel rule-based methods were developed to conduct automatic soil removal; and a newly proposed approach was used for image quality assessment; additionally, shadow masking and plot-level feature extraction were carried out. Our results show that the automated calibration, processing, storage, and analysis pipeline developed in this work can effectively handle massive amounts of hyperspectral data and address the urgent challenges related to the production of sustainable bioenergy and food crops, targeting methods to accelerate plant breeding for improving yield and biomass traits

Daily MODIS 500 m Reflectance Anisotropy Direct Broadcast (DB) Products for Monitoring Vegetation Phenology Dynamics

Land surface vegetation phenology is an efficient bio-indicator for monitoring ecosystem variation in response to changes in climatic factors. The primary objective of the current article is to examine the utility of the daily MODIS 500 m reflectance anisotropy direct broadcast (DB) product for monitoring the evolution of vegetation phenological trends over selected crop, orchard, and forest regions. Although numerous model-fitted satellite data have been widely used to assess the spatio-temporal distribution of land surface phenological patterns to understand phenological process and phenomena, current efforts to investigate the details of phenological trends, especially for natural phenological variations that occur on short time scales, are less well served by remote sensing challenges and lack of anisotropy correction in satellite data sources. The daily MODIS 500 m reflectance anisotropy product is employed to retrieve daily vegetation indices (VI) of a 1 year period for an almond orchard in California and for a winter wheat field in northeast China, as well as a 2 year period for a deciduous forest region in New Hampshire, USA. Compared with the ground records from these regions, the VI trajectories derived from the cloud-free and atmospherically corrected MODIS Nadir BRDF (bidirectional reflectance distribution function) adjusted reflectance (NBAR) capture not only the detailed footprint and principal attributes of the phenological events (such as flowering and blooming) but also the substantial inter-annual variability. This study demonstrates the utility of the daily 500 m MODIS reflectance anisotropy DB product to provide daily VI for monitoring and detecting changes of the natural vegetation phenology as exemplified by study regions comprising winter wheat, almond trees, and deciduous forest

시공간 해상도 향상을 통한 식생 변화 모니터링

학위논문(박사) -- 서울대학교대학원 : 환경대학원 협동과정 조경학, 2023. 2. 류영렬.육상 생태계에서 대기권과 생물권의 상호 작용을 이해하기 위해서는 식생 변화의 모니터링이 필요하다. 이 때, 위성영상은 지표면을 관측하여 식생지도를 제공할 수 있지만, 지표변화의 상세한 정보는 구름이나 위성 이미지의 공간 해상도에 의해 제한되었다. 또한 위성영상의 시공간 해상도가 식생지도를 통한 광합성 모니터링에 미치는 영향은 완전히 밝혀지지 않았다. 본 논문에서는 고해상도 식생 지도를 일단위로 생성하기 위성 영상의 시공간 해상도를 향상시키는 것을 목표로 하였다. 고해상도 위성영상을 활용한 식생 변화 모니터링을 시공간적으로 확장하기 위해 1) 정지궤도 위성을 활용한 영상융합을 통해 시간해상도 향상, 2) 적대적생성네트워크를 활용한 공간해상도 향상, 3) 시공간해상도가 높은 위성영상을 토지피복이 균질하지 않은 공간에서 식물 광합성 모니터링을 수행하였다. 이처럼, 위성기반 원격탐지에서 새로운 기술이 등장함에 따라 현재 및 과거의 위성영상은 시공간 해상도 측면에서 향상되어 식생 변화의 모니터링 할 수 있다. 제2장에서는 정지궤도위성영상을 활용하는 시공간 영상융합으로 식물의 광합성을 모니터링 했을 때, 시간해상도가 향상됨을 보였다. 시공간 영상융합 시, 구름탐지, 양방향 반사 함수 조정, 공간 등록, 시공간 융합, 시공간 결측치 보완 등의 과정을 거친다. 이 영상융합 산출물은 경작관리 등으로 식생 지수의 연간 변동이 큰 두 장소(농경지와 낙엽수림)에서 평가하였다. 그 결과, 시공간 영상융합 산출물은 결측치 없이 현장관측을 예측하였다 (R2 = 0.71, 상대 편향 = 5.64% 농경지; R2 = 0.79, 상대 편향 = -13.8%, 활엽수림). 시공간 영상융합은 식생 지도의 시공간 해상도를 점진적으로 개선하여, 식물 생장기동안 위성영상이 현장 관측을 과소 평가를 줄였다. 영상융합은 높은 시공간 해상도로 광합성 지도를 일간격으로 생성하기에 이를 활용하여 위성 영상의 제한된 시공간 해상도로 밝혀지지 않은 식물변화의 과정을 발견하길 기대한다. 식생의 공간분포은 정밀농업과 토지 피복 변화 모니터링을 위해 필수적이다. 고해상도 위성영상으로 지구 표면을 관측하는 것을 용이하게 해졌다. 특히 Planet Fusion은 초소형위성군 데이터를 최대한 활용해 데이터 결측이 없는 3m 공간 해상도의 지표 표면 반사도이다. 그러나 과거 위성 센서(Landsat의 경우 30~60m)의 공간 해상도는 식생의 공간적 변화를 상세 분석하는 것을 제한했다. 제3장에서는 Landsat 데이터의 공간 해상도를 향상하기 위해 Planet Fusion 및 Landsat 8 데이터를 사용하여 이중 적대적 생성 네트워크(the dual RSS-GAN)를 학습시켜, 고해상도 정규화 식생 지수(NDVI)와 식물 근적외선 반사(NIRv)도를 생성하는 한다. 타워기반 현장 식생지수(최대 8년)와 드론기반 초분광지도로 the dual RSS-GAN의 성능을 대한민국 내 두 대상지(농경지와 활엽수림)에서 평가했다. The dual RSS-GAN은 Landsat 8 영상의 공간해상도를 향상시켜 공간 표현을 보완하고 식생 지수의 계절적 변화를 포착했다(R2> 0.96). 그리고 the dual RSS-GAN은 Landsat 8 식생 지수가 현장에 비해 과소 평가되는 것을 완화했다. 현장 관측에 비해 이중 RSS-GAN과 Landsat 8의 상대 편향 값 각각 -0.8% 에서 -1.5%, -10.3% 에서 -4.6% 였다. 이러한 개선은 Planet Fusion의 공간정보를 이중 RSS-GAN로 학습하였기에 가능했다. 헤당 연구 결과는 Landsat 영상의 공간 해상도를 향상시켜 숨겨진 공간 정보를 제공하는 새로운 접근 방식이다. 고해상도에서 식물 광합성 지도는 토지피복이 복잡한 공간에서 탄소 순환 모니터링시 필수적이다. 그러나 Sentinel-2, Landsat 및 MODIS와 같이 태양 동조 궤도에 있는 위성은 공간 해상도가 높거나 시간 해상도 높은 위성영상만 제공할 수 있다. 최근 발사된 초소형위성군은 이러한 해상도 한계을 극복할 수 있다. 특히 Planet Fusion은 초소형위성 자료의 시공간 해상도로 지표면을 관측할 수 있다. 4장에서, Planet Fusion 지표반사도를 이용하여 식생에서 반사된 근적외선 복사(NIRvP)를 3m 해상도 지도를 일간격으로 생성했다. 그런 다음 미국 캘리포니아주 새크라멘토-샌 호아킨 델타의 플럭스 타워 네트워크 데이터와 비교하여 식물 광합성을 추정하기 위한 NIRvP 지도의 성능을 평가하였다. 전체적으로 NIRvP 지도는 습지의 잦은 수위 변화에도 불구하고 개별 대상지의 식물 광합성의 시간적 변화를 포착하였다. 그러나 대상지 전체에 대한 NIRvP 지도와 식물 광합성 사이의 관계는 NIRvP 지도를 플럭스 타워 관측범위와 일치시킬 때만 높은 상관관계를 보였다. 관측범위를 일치시킬 경우, NIRvP 지도는 식물 광합성을 추정하는 데 있어 현장 NIRvP보다 우수한 성능을 보였다. 이러한 성능 차이는 플럭스 타워 관측범위를 일치시킬 때, 연구 대상지 간의 NIRvP-식물 광합성 관계의 기울기가 일관성을 보였기 때문이다. 본 연구 결과는 위성 관측을 플럭스 타워 관측범위와 일치시키는 것의 중요성을 보여주고 높은 시공간 해상도로 식물 광합성을 원격으로 모니터링하는 초소형위성군 자료의 잠재력을 보여준다.Monitoring changes in terrestrial vegetation is essential to understanding interactions between atmosphere and biosphere, especially terrestrial ecosystem. To this end, satellite remote sensing offer maps for examining land surface in different scales. However, the detailed information was hindered under the clouds or limited by the spatial resolution of satellite imagery. Moreover, the impacts of spatial and temporal resolution in photosynthesis monitoring were not fully revealed. In this dissertation, I aimed to enhance the spatial and temporal resolution of satellite imagery towards daily gap-free vegetation maps with high spatial resolution. In order to expand vegetation change monitoring in time and space using high-resolution satellite images, I 1) improved temporal resolution of satellite dataset through image fusion using geostationary satellites, 2) improved spatial resolution of satellite dataset using generative adversarial networks, and 3) showed the use of high spatiotemporal resolution maps for monitoring plant photosynthesis especially over heterogeneous landscapes. With the advent of new techniques in satellite remote sensing, current and past datasets can be fully utilized for monitoring vegetation changes in the respect of spatial and temporal resolution. In Chapter 2, I developed the integrated system that implemented geostationary satellite products in the spatiotemporal image fusion method for monitoring canopy photosynthesis. The integrated system contains the series of process (i.e., cloud masking, nadir bidirectional reflectance function adjustment, spatial registration, spatiotemporal image fusion, spatial gap-filling, temporal-gap-filling). I conducted the evaluation of the integrated system over heterogeneous rice paddy landscape where the drastic land cover changes were caused by cultivation management and deciduous forest where consecutive changes occurred in time. The results showed that the integrated system well predict in situ measurements without data gaps (R2 = 0.71, relative bias = 5.64% at rice paddy site; R2 = 0.79, relative bias = -13.8% at deciduous forest site). The integrated system gradually improved the spatiotemporal resolution of vegetation maps, reducing the underestimation of in situ measurements, especially during peak growing season. Since the integrated system generates daily canopy photosynthesis maps for monitoring dynamics among regions of interest worldwide with high spatial resolution. I anticipate future efforts to reveal the hindered information by the limited spatial and temporal resolution of satellite imagery. Detailed spatial representations of terrestrial vegetation are essential for precision agricultural applications and the monitoring of land cover changes in heterogeneous landscapes. The advent of satellite-based remote sensing has facilitated daily observations of the Earths surface with high spatial resolution. In particular, a data fusion product such as Planet Fusion has realized the delivery of daily, gap-free surface reflectance data with 3-m pixel resolution through full utilization of relatively recent (i.e., 2018-) CubeSat constellation data. However, the spatial resolution of past satellite sensors (i.e., 30–60 m for Landsat) has restricted the detailed spatial analysis of past changes in vegetation. In Chapter 3, to overcome the spatial resolution constraint of Landsat data for long-term vegetation monitoring, we propose a dual remote-sensing super-resolution generative adversarial network (dual RSS-GAN) combining Planet Fusion and Landsat 8 data to simulate spatially enhanced long-term time-series of the normalized difference vegetation index (NDVI) and near-infrared reflectance from vegetation (NIRv). We evaluated the performance of the dual RSS-GAN against in situ tower-based continuous measurements (up to 8 years) and remotely piloted aerial system-based maps of cropland and deciduous forest in the Republic of Korea. The dual RSS-GAN enhanced spatial representations in Landsat 8 images and captured seasonal variation in vegetation indices (R2 > 0.95, for the dual RSS-GAN maps vs. in situ data from all sites). Overall, the dual RSS-GAN reduced Landsat 8 vegetation index underestimations compared with in situ measurements; relative bias values of NDVI ranged from −3.2% to 1.2% and −12.4% to −3.7% for the dual RSS-GAN and Landsat 8, respectively. This improvement was caused by spatial enhancement through the dual RSS-GAN, which captured fine-scale information from Planet Fusion. This study presents a new approach for the restoration of hidden sub-pixel spatial information in Landsat images. Mapping canopy photosynthesis in both high spatial and temporal resolution is essential for carbon cycle monitoring in heterogeneous areas. However, well established satellites in sun-synchronous orbits such as Sentinel-2, Landsat and MODIS can only provide either high spatial or high temporal resolution but not both. Recently established CubeSat satellite constellations have created an opportunity to overcome this resolution trade-off. In particular, Planet Fusion allows full utilization of the CubeSat data resolution and coverage while maintaining high radiometric quality. In Chapter 4, I used the Planet Fusion surface reflectance product to calculate daily, 3-m resolution, gap-free maps of the near-infrared radiation reflected from vegetation (NIRvP). I then evaluated the performance of these NIRvP maps for estimating canopy photosynthesis by comparing with data from a flux tower network in Sacramento-San Joaquin Delta, California, USA. Overall, NIRvP maps captured temporal variations in canopy photosynthesis of individual sites, despite changes in water extent in the wetlands and frequent mowing in the crop fields. When combining data from all sites, however, I found that robust agreement between NIRvP maps and canopy photosynthesis could only be achieved when matching NIRvP maps to the flux tower footprints. In this case of matched footprints, NIRvP maps showed considerably better performance than in situ NIRvP in estimating canopy photosynthesis both for daily sum and data around the time of satellite overpass (R2 = 0.78 vs. 0.60, for maps vs. in situ for the satellite overpass time case). This difference in performance was mostly due to the higher degree of consistency in slopes of NIRvP-canopy photosynthesis relationships across the study sites for flux tower footprint-matched maps. Our results show the importance of matching satellite observations to the flux tower footprint and demonstrate the potential of CubeSat constellation imagery to monitor canopy photosynthesis remotely at high spatio-temporal resolution.Chapter 1. Introduction 2 1. Background 2 1.1 Daily gap-free surface reflectance using geostationary satellite products 2 1.2 Monitoring past vegetation changes with high-spatial-resolution 3 1.3 High spatiotemporal resolution vegetation photosynthesis maps 4 2. Purpose of Research 4 Chapter 2. Generating daily gap-filled BRDF adjusted surface reflectance product at 10 m resolution using geostationary satellite product for monitoring daily canopy photosynthesis 6 1. Introduction 6 2. Methods 11 2.1 Study sites 11 2.2 In situ measurements 13 2.3 Satellite products 14 2.4 Integrated system 17 2.5 Canopy photosynthesis 21 2.6 Evaluation 23 3. Results and discussion 24 3.1 Comparison of STIF NDVI and NIRv with in situ NDVI and NIRv 24 3.2 Comparison of STIF NIRvP with in situ NIRvP 28 4. Conclusion 31 Chapter 3. Super-resolution of historic Landsat imagery using a dual Generative Adversarial Network (GAN) model with CubeSat constellation imagery for monitoring vegetation changes 32 1. Introduction 32 2. Methods 38 2.1 Real-ESRGAN model 38 2.2 Study sites 40 2.3 In situ measurements 42 2.4 Vegetation index 44 2.5 Satellite data 45 2.6 Planet Fusion 48 2.7 Dual RSS-GAN via fine-tuned Real-ESRGAN 49 2.8 Evaluation 54 3. Results 57 3.1 Comparison of NDVI and NIRv maps from Planet Fusion, Sentinel 2 NBAR, and Landsat 8 NBAR data with in situ NDVI and NIRv 57 3.2 Comparison of dual RSS-SRGAN model results with Landsat 8 NDVI and NIRv 60 3.3 Comparison of dual RSS-GAN model results with respect to in situ time-series NDVI and NIRv 63 3.4 Comparison of the dual RSS-GAN model with NDVI and NIRv maps derived from RPAS 66 4. Discussion 70 4.1 Monitoring changes in terrestrial vegetation using the dual RSS-GAN model 70 4.2 CubeSat data in the dual RSS-GAN model 72 4.3 Perspectives and limitations 73 5. Conclusion 78 Appendices 79 Supplementary material 82 Chapter 4. Matching high resolution satellite data and flux tower footprints improves their agreement in photosynthesis estimates 85 1. Introduction 85 2. Methods 89 2.1 Study sites 89 2.2 In situ measurements 92 2.3 Planet Fusion NIRvP 94 2.4 Flux footprint model 98 2.5 Evaluation 98 3. Results 105 3.1 Comparison of Planet Fusion NIRv and NIRvP with in situ NIRv and NIRvP 105 3.2 Comparison of instantaneous Planet Fusion NIRv and NIRvP with against tower GPP estimates 108 3.3 Daily GPP estimation from Planet Fusion -derived NIRvP 114 4. Discussion 118 4.1 Flux tower footprint matching and effects of spatial and temporal resolution on GPP estimation 118 4.2 Roles of radiation component in GPP mapping 123 4.3 Limitations and perspectives 126 5. Conclusion 133 Appendix 135 Supplementary Materials 144 Chapter 5. Conclusion 153 Bibliography 155 Abstract in Korea 199 Acknowledgements 202박

Estimation of land surface directional emissivity in mid-infrared channel around 4.0 mu m from MODIS data

This work addressed the estimate of the directional emissivity in the mid-infrared (MIR) channel around 4.0 mu m from MODIS data. A series of bidirectional reflectances in MODIS channel 22 (3.97 mu m) were retrieved using the method developed by Tang and Li (Int. J. Remote Sens. 29, 4907, 2008) and then were used to estimate the directional emissivity in this channel with the aid of the BRDF model modified by Jiang and Li (Opt. Express 16, 19310, 2008). To validate the estimated directional emissivity, a cross-comparison of MODIS derived emissivities in channel 22 using the proposed method were performed with those provided by the MODIS land surface temperature/emissivity product MYD11B1 data. The results show that the proposed method for estimating the directional emissivity in MIR channel gives results comparable to those of MYD11B1 product with a Mean Error of -0.007 and a Root Mean Square Error of 0.024. (C) 2009 Optical Society of Americ

Daily MODIS 500 m Reflectance Anisotropy Direct Broadcast (DB) Products for Monitoring Vegetation Phenology Dynamics

Land surface vegetation phenology is an efficient bio-indicator for monitoring ecosystem variation in response to changes in climatic factors. The primary objective of the current article is to examine the utility of the daily MODIS 500 m reflectance anisotropy direct broadcast (DB) product for monitoring the evolution of vegetation phenological trends over selected crop, orchard, and forest regions. Although numerous model-fitted satellite data have been widely used to assess the spatio-temporal distribution of land surface phenological patterns to understand phenological process and phenomena, current efforts to investigate the details of phenological trends, especially for natural phenological variations that occur on short time scales, are less well served by remote sensing challenges and lack of anisotropy correction in satellite data sources. The daily MODIS 500 m reflectance anisotropy product is employed to retrieve daily vegetation indices (VI) of a 1 year period for an almond orchard in California and for a winter wheat field in northeast China, as well as a 2 year period for a deciduous forest region in New Hampshire, USA. Compared with the ground records from these regions, the VI trajectories derived from the cloud-free and atmospherically corrected MODIS Nadir BRDF (bidirectional reflectance distribution function) adjusted reflectance (NBAR) capture not only the detailed footprint and principal attributes of the phenological events (such as flowering and blooming) but also the substantial inter-annual variability. This study demonstrates the utility of the daily 500 m MODIS reflectance anisotropy DB product to provide daily VI for monitoring and detecting changes of the natural vegetation phenology as exemplified by study regions comprising winter wheat, almond trees, and deciduous forest

System Engineering Analyses for the Study of Future Multispectral Land Imaging Satellite Sensors for Vegetation Monitoring

Vegetation monitoring is one of the key applications of earth observing systems. Landsat data have spatial resolution of 30 meters, moderate temporal coverage, and reasonable spectral sampling to capture key vegetation features. These characteristics of Landsat make it a good candidate for generating vegetation monitoring products. Recently, the next satellite in the Landsat series has been under consideration and different concepts have been proposed. In this research, we studied the impact on vegetation monitoring of two proposed potential design concepts: a wider field-of-view (FOV) instrument and the addition of red-edge spectral band(s). Three aspects were studied in this thesis: First, inspired by the potential wider FOV design, the impacts of a detector relative spectral response (RSR) central wavelength shift effect at high angles of incidence (AOI) on the radiance signal were studied and quantified. Results indicate: 1) the RSR shift effect is band-dependent and more significant in the green, red and SWIR 2 bands; 2) At high AOI, the impact of the RSR shift effect will exceed sensor noise specifications in all bands except the SWIR 1 band; and 3) The RSR shift will cause SWIR2 band more to be sensitive to atmospheric conditions. Second, also inspired by the potential wider FOV design, the impacts of the potential new wider angular observations on vegetation monitoring scientific products were studied. Both crop classification and biophysical quantity retrieval applications were studied using the simulation code DIRSIG and the canopy radiative transfer model PROSAIL. It should be noted that the RSR shift effect was also considered. Results show that for single view observation based analysis, the higher view angular observations have limited influence on both applications. However, for situations where two different angular observations are available potentially from two platforms, up to 4% improvement for crop classification and 2.9% improvement for leaf chlorophyll content retrieval were found. Third, to quantify the benefits of a potential new design with red-edge band(s), the impact of adding red-edge spectral band(s) in future Landsat instruments on agroecosystem leaf area index (LAI) and canopy chlorophyll content (CCC) retrieval were studied using a real dataset. Three major retrieval approaches were tested, results show that a potential new spectral band located between the Landsat-8 Operational Land Imager (OLI) red and NIR bands slightly improved the retrieval accuracy (LAI: R2 of 0.787 vs. 0.810 for empirical vegetation index regression approach, 0.806 vs. 0.828 for look-up-table inversion approach, and 0.925 vs. 0.933 for machine learning approach; CCC: R2 of 0.853 vs. 0.875 for empirical vegetation index regression approach, 0.500 vs. 0.570 for look-up-table inversion approach, and 0.854 vs. 0.887 for machine learning approach). In general, for the potential wider FOV design, the RSR shift effect was found to cause noticable radiance signal difference that is higher than detector noise in all OLI bands except SWIR1 band, which is not observed in the current OLI design with its 15 degree FOV. Also both the new wider angular observations and potential red-edge band(s) were found to slightly improve the vegetation monitoring product accuracy. In the future, the RSR shift effect in other optical designs should be evaluated since this study assumed the angle reaching the filter array is the same as the angle reaching the sensor. In addition to improve the accuracy of the off angle imaging study, a 3D vegetation geometry model should be explored for vegetation monitoring related studies instead of the 2D PROSAIL model used in this thesis

Evaluation of sensor, environment and operational factors impacting the use of multiple sensor constellations for long term resource monitoring

Moderate resolution remote sensing data offers the potential to monitor the long and short term trends in the condition of the Earth’s resources at finer spatial scales and over longer time periods. While improved calibration (radiometric and geometric), free access (Landsat, Sentinel, CBERS), and higher level products in reflectance units have made it easier for the science community to derive the biophysical parameters from these remotely sensed data, a number of issues still affect the analysis of multi-temporal datasets. These are primarily due to sources that are inherent in the process of imaging from single or multiple sensors. Some of these undesired or uncompensated sources of variation include variation in the view angles, illumination angles, atmospheric effects, and sensor effects such as Relative Spectral Response (RSR) variation between different sensors. The complex interaction of these sources of variation would make their study extremely difficult if not impossible with real data, and therefore, a simulated analysis approach is used in this study. A synthetic forest canopy is produced using the Digital Imaging and Remote Sensing Image Generation (DIRSIG) model and its measured BRDFs are modeled using the RossLi canopy BRDF model. The simulated BRDF matches the real data to within 2% of the reflectance in the red and the NIR spectral bands studied. The BRDF modeling process is extended to model and characterize the defoliation of a forest, which is used in factor sensitivity studies to estimate the effect of each factor for varying environment and sensor conditions. Finally, a factorial experiment is designed to understand the significance of the sources of variation, and regression based analysis are performed to understand the relative importance of the factors. The design of experiment and the sensitivity analysis conclude that the atmospheric attenuation and variations due to the illumination angles are the dominant sources impacting the at-sensor radiance

Changes in Tall Shrub Abundance on the North Slope of Alaska, 2000-2010

The observed greening of Arctic vegetation and the expansion of shrubs in the last few decades has likely had profound implications for the tundra ecosystem, including feedbacks to climate. Uncertainty surrounding the magnitude, direction, and implications of this vegetation shift calls for monitoring of vegetation structural parameters, such as fractional cover of shrubs. Due to the extent of the North Slope of Alaska and its extreme environments, remote sensing may be the most suitable tool to produce wall-to-wall fractional shrub cover maps for the entire region, however, most regional maps have relied on vegetation indices or needed many years worth of data to cover the whole region. Here, a new mapping approach is presented that uses satellite imagery from the Multi-angle Imaging SpectroRadiometer (MISR) sensor and some landscape variables to predict tall shrub (\u3e 0.5 m) cover with the ultimate goal of evaluating temporal changes in tall shrub fractional cover during the period of 2010-2000. Specifically, we: 1) undertook two field surveys in the North Slope of Alaska to obtain estimates of tall shrub cover, canopy height, crown radius, and total number of shrubs at 26 sites (250 m × 250 m each); 2) evaluated the ability of the semi-automated image interpretation algorithm CANAPI - CANopy Analysis from Panchromatic Imagery, to derive structural data for tall (\u3e 0.5 m) shrubs in the Arctic; 3) constructed a robust reference database with estimates of shrub structural parameters; 4) trained and validated the boosted regression tree model to predict tall shrub fractional cover from moderate resolution imagery; 5) created the 2000 and the 2010 tall shrub fractional cover map for the North Slope of Alaska; and 6) evaluated the changes in shrub abundance during the period 2010-2000 in the North Slope of Alaska. Results from the field surveys suggested that tall shrub fractional cover was less than 5% at 250 m scales. The evaluation of the CANAPI algorithm showed that CANAPI could successfully retrieve fractional cover (R2 = 0.83, P \u3c 0.001), mean crown radius (R2 = 0.81, P \u3c 0.001), and total number of shrubs (R2 = 0.54, P \u3c 0.001) from very-high resolution imagery. As a result, a robust reference database was constructed with estimates of tall shrub fractional cover, canopy radius, and total number of shrubs for 1,039 sites across the domain of the North Slope. After the training and validation of the Boosted Regression Tree (BRT), the best model used 14 predictor variables and explained 52% of the variation in the response variable, fractional cover. The red reflectance, slope, nadir Bidirectional Reflectance Distribution Function (BRDF) adjusted weight of determination, and isotropic scattering kernel were the variables more often used to generate the regression trees, and therefore they contributed the most to the model. The trained BRT model was used to construct the tall shrub fractional cover map for the year 2000 and 2010 using moderate resolution imagery. The maps revealed that cover ranged from 0.00 to 0.21 and about 75% of the sites had a fractional cover less than 0.013. High cover values were predicted along floodplains, creeks, and sloped terrain. The 2000 MISR-derived fractional cover map presented here outperformed the 2000 Landsat-derived tall shrub fractional cover map when compared to the robust validation data set (R2= 0.38, Root Mean Square Error (RMSE) = 0.08). Temporal comparisons of tall shrub abundance in the MISR-derived maps suggested that shrubs expanded during the period 2000-2010. The extent of the area that unequivocally experienced a robust change in tall shrub cover was less than 1 % (1,487 km2) of the total area of the North Slope of Alaska (213,090 km2). It is possible that tall shrubs may have expanded throughout a larger area but there is insufficient precision in the MISR-based estimates to make an unequivocal determination. Nevertheless, it seems that there was a positive trend toward an increase in shrub cover considering that 95% of the locations that had a robust change saw an increase. The tall shrub cover expansion rate varied between 0.006 yr-1 and 0.017 yr-1, being higher along the forest-tundra ecotone, north of the Brooks Range. More research is necessary to determine if the increase in cover corresponded to the advance of the tree line, or to the expansion of the tall shrubs, or both

Influence of forest floor vegetation on the total forest reflectance and its implications for LAI estimation using vegetation indices

Recently a simple analytic canopy bidirectional reflectance factor (BRF) model based on the spectral invariants theory was presented. The model takes into account that the recollision probability in the forest canopy is different for the first scattering than the later ones. Here this model is extended to include the forest floor contribution to the total forest BRF. The effect of the understory vegetation on the total forest BRF as well as on the simple ratio (SR) and the normalized difference (NDVI) vegetation indices is demonstrated for typical cases of boreal forest. The relative contribution of the forest floor to the total BRF was up to 69 % in the red wavelength range and up to 54 % in the NIR wavelength range. Values of SR and NDVI for the forest and the canopy differed within 10 % and 30 % in red and within 1 % and 10 % in the NIR wavelength range. The relative variation of the BRF with the azimuth and view zenith angles was not very sensitive to the forest floor vegetation. Hence, linear correlation of the modelled total BRF and the Ross-thick kernel was strong for dense forests (R2 > 0.9). The agreement between modelled BRF and satellite-based reflectance values was good when measured LAI, clumping index and leaf single scattering albedo values for a boreal forest were used as input to the model.Hiljattain on esitetty yksinkertainen analyyttinen puuston kaksisuuntaisen heijastuskertoimen (BRF) malli, joka perustuu spketristä riippumattomien parametrien teoriaan. Mallissa otetaan huomioon, että fotonin uudelleen siroamisen todennäköisyys metsässä poikkeaa ensimmäisellä kerralla sen myöhemmistä arvoista. Tässä tutkimuksessa mallia on edelleen kehitetty siten, että siinä huomioidaan metsän pohjan osuus metsän BRF:stä. Aluskasvillisuuden vaikutusta BRF:ään ja kasvillisuusindekseihin SR (yksinkertainen suhde) ja NDVI (normalisoitu kasvillisuuden erotusindeksi) havainnollistetaan esimerkeillä tyypillisestä boreaalisesta metsästä. Metsän pohjan suhteellinen osuus BRF:stä ulottui 69 prosenttiin punaisen alueen aallonpituuksilla ja 54 prosenttiin lähiinfrapunan aallonpituusalueella. Metsälle laskettujen SR:n ja NDVI:n arvot poikkesivat pelkälle puustolle lasketuista vastaavista arvoista 10 % ja 30 % punaisella aallonpituusalueella ja 1 % ja 10 % lähi-infrapunan aallonpituusalueella. BRF:n suhteellinen muutos katselukulmien vaihdellessa ei ollut kovin herkkä metsän pohjakasvillisuudelle. Siten mallinnettu metsän BRF oli lineaarisesti verrannollinen Ross:n tiheän metsän kerneliin (selitysaste > 0.9). Mallinnettu BRF ja satellidatasta peräisin olevat reflektanssiarvot vastasivat hyvin toisiaan, kun mallin syöttötietoina käytettiin boreaalisen metsän mitattuja lehtialaindeksin, ryhmittyneisyysindeksin ja lehden albedon arvoja