Monitoring Land Surface Albedo and Vegetation Dynamics Using High Spatial and Temporal Resolution Synthetic Time Series from Landsat and the MODIS BRDF/NBAR/Albedo Product

Seasonal vegetation phenology can significantly alter surface albedo which in turn affects the global energy balance and the albedo warmingcooling feedbacks that impact climate change. To monitor and quantify the surface dynamics of heterogeneous landscapes, high temporal and spatial resolution synthetic time series of albedo and the enhanced vegetation index (EVI) were generated from the 500-meter Moderate Resolution Imaging Spectroradiometer (MODIS) operational Collection V006 daily BRDF (Bidirectional Reflectance Distribution Function) / NBAR (Nadir BRDF-Adjusted Reflectance) / albedo products and 30-meter Landsat 5 albedo and near-nadir reflectance data through the use of the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). The traditional Landsat Albedo (Shuai et al., 2011) makes use of the MODIS BRDFAlbedo products (MCD43) by assigning appropriate BRDFs from coincident MODIS products to each Landsat image to generate a 30-meter Landsat albedo product for that acquisition date. The available cloud free Landsat 5 albedos (due to clouds, generated every 16 days at best) were used in conjunction with the daily MODIS albedos to determine the appropriate 30-meter albedos for the intervening daily time steps in this study. These enhanced daily 30-meter spatial resolution synthetic time series were then used to track albedo and vegetation phenology dynamics over three Ameriflux tower sites (Harvard Forest in 2007, Santa Rita in 2011 and Walker Branch in 2005). These Ameriflux sites were chosen as they are all quite nearby new towers coming on line for the National Ecological Observatory Network (NEON), and thus represent locations which will be served by spatially paired albedo measures in the near future. The availability of data from the NEON towers will greatly expand the sources of tower albedometer data available for evaluation of satellite products. At these three Ameriflux tower sites the synthetic time series of broadband shortwave albedos were evaluated using the tower albedo measurements with a Root Mean Square Error (RMSE) less than 0.013 and a bias within the range of 0.006. These synthetic time series provide much greater spatial detail than the 500 meter gridded MODIS data, especially over more heterogeneous surfaces, which improves the efforts to characterize and monitor the spatial variation across species and communities. The mean of the difference between maximum and minimum synthetic time series of albedo within the MODIS pixels over a subset of satellite data of Harvard Forest (16 kilometers by 14 kilometers) was as high as 0.2 during the snow-covered period and reduced to around 0.1 during the snow-free period. Similarly, we have used STARFM to also couple MODIS Nadir BRDF-Adjusted Reflectances (NBAR) values with Landsat 5 reflectances to generate daily synthetic times series of NBAR and thus Enhanced Vegetation Index (NBAR-EVI) at a 30-meter resolution. While normally STARFM is used with directional reflectances, the use of the view angle corrected daily MODIS NBAR values will provide more consistent time series. These synthetic times series of EVI are shown to capture seasonal vegetation dynamics with finer spatial and temporal details, especially over heterogeneous land surfaces

Assessing change in the Earth's land surface albedo with moderate resolution satellite imagery

Land surface albedo describes the proportion of incident solar radiant flux that is reflected from the Earth's surface and therefore is a crucial parameter in modeling and monitoring attempts to capture the current climate, hydrological, and biogeochemical cycles and predict future scenarios. Due to the temporal variability and spatial heterogeneity of land surface albedo, remote sensing offers the only realistic method of monitoring albedo on a global scale. While the distribution of bright, highly reflective surfaces (clouds, snow, deserts) govern the vast majority of the fluctuation, variations in the intrinsic surface albedo due to natural and human disturbances such as urban development, fire, pests, harvesting, grazing, flooding, and erosion, as well as the natural seasonal rhythm of vegetation phenology, play a significant role as well. The development of times series of global snow-free and cloud-free albedo from remotely sensed observations over the past decade and a half offers a unique opportunity to monitor and assess the impact of these alterations to the Earth's land surface. By utilizing multiple satellite records from the MODerate-resolution Imaging Spectroradiometer (MODIS), the Multi-angle Imaging Spectroradiometer (MISR) and the Visible Infrared Imaging Radiometer Suite (VIIRS) instruments, and developing innovative spectral conversion coefficients and temporal gap-filling strategies, it has been possible to utilize the strengths of the various sensors to improve the spatial and temporal coverage of global land surface albedo retrievals. The availability of these products is particularly important in tropical regions where cloud cover obscures the forest for significant periods. In the Amazon, field ecologists have noted that some areas of the forest ecosystem respond rapidly with foliage growth at the beginning of the dry season, when sunlight can finally penetrate fully to the surface and have suggested this phenomenon can continue until reductions in water availability (particularly in times of drought) impact the growth cycle. While it has been difficult to capture this variability from individual optical satellite sensors, the temporally gap-filled albedo products developed during this research are used in a case study to monitor the Amazon during the dry season and identify the extent of these regions of foliage growth

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

학위논문(박사) -- 서울대학교대학원 : 환경대학원 협동과정 조경학, 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박

Observing and modeling climate controls and feedbacks on vegetation phenology at local-to-continental scales

Vegetation phenology controls seasonal variation in ecosystem processes and exerts important controls on land-atmosphere exchanges of carbon, water, and energy. However, the ecological processes and interactions between climate and vegetation that control phenology and associated feedbacks to the atmosphere are not fully understood. In this dissertation, I use remote sensing in combination with climate and ecological data to improve understanding of biophysical controls and feedbacks between vegetation phenology and the atmosphere in temperate forest ecosystems of North America. In the first part of this dissertation, I evaluate the agreement and characterize the similarities and differences between land surface phenology products from two remote sensing instruments (MODIS and VIIRS) that are designed to provide long-term continuity of land surface phenology measurements at global scale. Results from this analysis indicate that the VIIRS land surface phenology product provides excellent continuity with the MODIS record despite subtle differences between each instrument and the algorithms used to generate each product. In the second part of this dissertation, a state-space Bayesian modeling framework is applied to seventeen years of MODIS and daily weather data to improve understanding of what controls the timing of springtime phenology in deciduous forests of temperate and boreal North America. Results show that photoperiod is more important in warmer regions than in colder regions, which contradicts a widely held hypothesis that photoperiod provides a key safety mechanism preventing early leaf-out during springtime. In the final part of this dissertation, I use a physically-based attribution method to quantify the relative importance of covarying surface biophysical and atmospheric variables in modifying the surface energy balance during springtime. Results show that the widely observed decrease in the Bowen ratio that occurs with leaf emergence is not solely attributable to changes in surface resistance caused by increasing leaf area during spring. Rather, observed changes in the Bowen ratio reflect the combined effects of changes in surface properties and atmospheric conditions. The results from this dissertation provide an improved foundation for long-term studies focused on observing and modeling springtime vegetation phenology and associated feedbacks to the atmosphere in deciduous forest ecosystems at local-to-continental scales

Vegetation Dynamics Revealed by Remote Sensing and Its Feedback to Regional and Global Climate

This book focuses on some significant progress in vegetation dynamics and their response to climate change revealed by remote sensing data. The development of satellite remote sensing and its derived products offer fantastic opportunities to investigate vegetation changes and their feedback to regional and global climate systems. Special attention is given in the book to vegetation changes and their drivers, the effects of extreme climate events on vegetation, land surface albedo associated with vegetation changes, plant fingerprints, and vegetation dynamics in climate modeling

Evaluating the Spatial Representativeness of the MODerate Resolution Image Spectroradiometer Albedo Product (MCD43) at AmeriFlux Sites

Land surface albedo is a key parameter in regulating surface radiation budgets. The gridded remote sensing albedo product often represents information concerning an area larger than the nominal spatial resolution because of the large viewing angles of the observations. It is essential to quantify the spatial representativeness of remote sensing products to better guide the sampling strategy in field experiments and match products from different sources. This study quantifies the spatial representativeness of the MODerate Resolution Image Spectroradiometer (MODIS) (collection V006) 500 m daily albedo product (MCD43A3) using the high-resolution product as intermediate data for different land cover types. A total of 1820 paired high-resolution Landsat Thematic Mapper (TM) and coarse-resolution (MODIS) albedo data from five land cover types were used. The TM albedo data was used as the spatial-complete high resolution data to evaluate the spatial representativeness of the MODIS albedo product. Semivarioagrams were estimated from 30 m Landsat data at different spatial scales. Surface heterogeneity was evaluated with sill value and relative coefficient of variation. The 30 m Landsat albedo data was aggregated to 450 m–1800 m using two different methods and compared with MODIS albedo product. The spatial representativeness of MODIS albedo product was determined according to the surface heterogeneity and the consistency of MODIS data and the aggregated TM value. Results indicated that for evergreen broadleaf forests, deciduous broadleaf forests, open shrub lands, woody savannas and grasslands, the MODIS 500 m daily albedo product represents a spatial scale of approximately 630 m. For mixed forests and croplands, the representative spatial scale was approximately 690 m. The difference obtained was primarily because of the complexity of the landscape structure. For mixed forests and croplands, the structure of the landscape was relatively complex due to the presence of different forest and plant types in the pixel area, whereas the other landscape structures were considerably simpler