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

A Cross Comparison of Spatiotemporally Enhanced Springtime Phenological Measurements From Satellites and Ground in a Northern U.S. Mixed Forest

Cross comparison of satellite-derived land surface phenology (LSP) and ground measurements is useful to ensure the relevance of detected seasonal vegetation change to the underlying biophysical processes. While standard 16-day and 250-m Moderate Resolution Imaging Spectroradiometer (MODIS) vegetation index (VI)-based springtime LSP has been evaluated in previous studies, it remains unclear whether LSP with enhanced temporal and spatial resolutions can capture additional details of ground phenology. In this paper, we compared LSP derived from 500-m daily MODIS and 30-m MODIS-Landsat fused VI data with landscape phenology (LP) in a northern U.S. mixed forest. LP was previously developed from intensively observed deciduous and coniferous tree phenology using an upscaling approach. Results showed that daily MODIS-based LSP consistently estimated greenup onset dates at the study area (625 m × 625 m) level with 4.48 days of mean absolute error (MAE), slightly better than that of using 16-day standard VI (4.63 days MAE). For the observed study areas, the time series with increased number of observations confirmed that post-bud burst deciduous tree phenology contributes the most to vegetation reflectance change. Moreover, fused VI time series demonstrated closer correspondences with LP at the community level (0.1-20 ha) than using MODIS alone at the study area level (390 ha). The fused LSP captured greenup onset dates for respective forest communities of varied sizes and compositions with four days of the overall MAE. This study supports further use of spatiotemporally enhanced LSP for more precise phenological monitoring

Monitoring effects of land cover change on biophysical drivers in rangelands using albedo

This paper explores the relationship between land cover change and albedo, recognized as a regulating ecosystems service. Trends and relationships between land cover change and surface albedo were quantified to characterise catchment water and carbon fluxes, through respectively evapotranspiration (ET) and net primary production (NPP). Moderate resolution imaging spectroradiometer (MODIS) and Landsat satellite data were used to describe trends at catchment and land cover change trajectory level. Peak season albedo was computed to reduce seasonal effects. Different trends were found depending on catchment land management practices, and satellite data used. Although not statistically significant, albedo, NPP, ET and normalised difference vegetation index (NDVI) were all correlated with rainfall. In both catchments, NPP, ET and NDVI showed a weak negative trend, while albedo showed a weak positive trend. Modelled land cover change was used to calculate future carbon storage and water use, with a decrease in catchment carbon storage and water use computed. Grassland, a dominant dormant land cover class, was targeted for land cover change by woody encroachment and afforestation, causing a decrease in albedo, while urbanisation and cultivation caused an increase in albedo. Land cover map error of fragmented transition classes and the mixed pixel effect, affected results, suggesting use of higher-resolution imagery for NPP and ET and albedo as a proxy for land cover

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

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

A General Method to Normalize Landsat Reflectance Data to Nadir BRDF Adjusted Reflectance

The Landsat satellites have been providing spectacular imagery of the Earth\u27s surface for over 40 years. However, they acquire images at view angles ±7.5° from nadir that cause small directional effects in the surface reflectance. There are also variations with solar zenith angle over the year that can cause apparent change in reflectance even if the surface properties remain constant. When Landsat data from adjoining paths, or from long time series are used, a model of the surface anisotropy is required to adjust all Landsat observations to a uniform nadir view (primarily for visual consistency, vegetation monitoring, or detection of subtle surface changes). Here a generalized approach is developed to provide consistent view angle corrections across the Landsat archive. While this approach is not applicable for generation of Landsat surface albedo, which requires a full characterization of the surface bidirectional reflectance distribution function (BRDF), or for correction to a constant solar illumination angle across a wide range of sun angles, it provides Landsat nadir BRDF-adjusted reflectance (NBAR) for a range of terrestrial monitoring applications. The Landsat NBAR is derived as the product of the observed Landsat reflectance and the ratio of the reflectances modeled using MODIS BRDF spectral model parameters for the observed Landsat and for a nadir view and fixed solar zenith geometry. In this study, a total of 567 conterminous United States (CONUS) January and July 2010 Landsat 5 Thematic Mapper (TM) and Landsat 7 Enhanced Thematic Mapper (ETM+) images that have swath edge overlapping paths sensed in alternating backscatter and forward scattering orientations were used. The average difference between Landsat 5 TM and Landsat 7 ETM+ surface reflectance in the forward and backward scatter directions at the overlapping Landsat scan edges was quantified. The CONUS July view zenith BRDF effects were about 0.02 in the Landsat visible bands, and about 0.03, 0.05 and 0.06, in the 2.1 μm, 1.6 μm and nearinfrared bands respectively. Comparisons of Landsat 5 TM and Landsat 7 ETM+ NBAR derived using MODIS BRDF spectral model parameters defined with respect to different spatial and temporal scales, and defined with respect to different land cover types, were undertaken. The results suggest that, because the BRDF shapes of different terrestrial surfaces are sufficiently similar over the narrow 15° Landsat field of view, a fixed set of MODIS BRDF spectral model parameters may be adequate for Landsat NBAR derivation with little sensitivity to the land cover type, condition, or surface disturbance. A fixed set of BRDF spectral model parameters, derived from a global year of highest quality snow-free MODIS BRDF product values, are provided so users may implement the described Landsat NBAR generation method

Landsat 5 Thematic Mapper Reflectance and NDVI 27-year Time Series Inconsistencies Due to Satellite Orbit Change

The Landsat 5 Thematic Mapper (TM) sensor provided the longest single mission terrestrial remote sensing data record but temporally sparse station keeping maneuvers meant that the Landsat 5 orbit changed over the 27 year mission life. Long-term Landsat 5 TM reflectance inconsistencies may be introduced by orbit change induced solar zenith variations combined with surface reflectance anisotropy, commonly described by the Bi-directional Reflectance Distribution Function (BRDF). This study quantifies the local overpass time and observed solar zenith angle changes for all the Landsat 5 TM images available at two latitudinally separated locations along the same north-south Landsat path (27) in Minnesota (row 26) and Texas (row 42). Over the 27 years the Landsat 5 orbit changed by nearly 1 h and resulted in changes in the Landsat 5 observed solar zenith angle of N10°. The Landsat 5 orbit was relatively stable from 1984 to 1994 and from 2007 to 2011, but changed rapidly from 1995 to 2000, and from 2003 to 2007. Rather than directly examine Landsat 5 TM reflectance time series a modelling approach was used. This was necessary because unambiguous separation of orbit change induced Landsat reflectance variations from other temporal variations is non-trivial. The impact of Landsat 5 orbit induced observed solar zenith angle variations on the red and near-infrared reflectance and derived normalized difference vegetation index (NDVI) values were modelled with respect to different Moderate-Resolution Imaging Spectroradiometer (MODIS) BRDF land cover types. Synthetic nadir BRDF-adjusted reflectance (NBAR) for the Landsat 5 TM observed and a modelled reference year 2011 solar zenith were compared over the 27 years of acquisitions. Ordinary least squares linear regression fits of the NBAR difference values as a function of the acquisition date indicated an increasing trend in red and near-infrared NBAR and a decreasing trend in NDVI NBAR due to orbit changes. The trends are statistically significant but small, no more than 0.0006 NDVI/year. Comparison of certain years of Landsat 5 data may result in significant reflectance and NDVI differences due only to Landsat 5 orbit changes and cause spurious detection of “browning” vegetation events and underestimation of greening trends. The greatest differences will occur when 1995 Landsat 5 TM data are compared with 2007 to 2011 data; NDVI values could be up to 0.11 greater in 1995 than in 2011 for anisotropic land cover types and up to 0.05 greater for average CONUS land cover types. A smaller number of Landsat 5 TM images were also examined and provide support for the modelled based findings. The paper concludes with a discussion of the implications of the research findings for Landsat 5 TM time series analyses

Continuity of Landsat Obersvations: Short Term Considerations

As of writing in mid-2010, both Landsat-5 and -7 continue to function, with sufficient fuel to enable data collection until the launch of the Landsat Data Continuity Mission (LDCM) scheduled for December of 2012. Failure of one or both of Landsat-5 or -7 may result in a lack of Landsat data for a period of time until the 2012 launch. Although the potential risk of a component failure increases the longer the sensor\u27s design life is exceeded, the possible gap in Landsat data acquisition is reduced with each passing day and the risk of Landsat imagery being unavailable diminishes for all except a handful of applications that are particularly data demanding. Advances in Landsat data compositing and fusion are providing opportunities to address issues associated with Landsat-7 SLC-off imagery and to mitigate a potential acquisition gap through the integration of imagery from different sensors. The latter will likely also provide short-term, regional solutions to application-specific needs for the continuity of Landsat-like observations. Our goal in this communication is not to minimize the community\u27s concerns regarding a gap in Landsat observations, but rather to clarify how the current situation has evolved and provide an up-to-date understanding of the circumstances, implications, and mitigation options related to a potential gap in the Landsat data record

Remotely sensed albedo allows the identification of two ecosystem states along aridity gradients in Africa

Empirical verification of multiple states in drylands is scarce, impeding the design of indicators to anticipate the onset of desertification. Remote sensing‐derived indicators of ecosystem states are gaining new ground due to the possibilities they bring to be applied inexpensively over large areas. Remotely sensed albedo has been often used to monitor drylands due to its close relationship with ecosystem status and climate. Here, we used a space‐for‐time‐substitution approach to evaluate whether albedo (averaged from 2000 to 2016) can identify multiple ecosystem states in African drylands spanning from the Saharan desert to tropical Africa. By using latent class analysis, we found that albedo showed two states (low and high; the cut‐off level was 0.22 at the shortwave band). Potential analysis revealed that albedo exhibited an abrupt and discontinuous increase with increased aridity (1 − [precipitation/potential evapotranspiration]). The two albedo states co‐occurred along aridity values ranging from 0.72 to 0.78, during which vegetation cover exhibited a rapid, continuous decrease from ~90% to ~50%. At aridity values of 0.75, the low albedo state started to exhibit less attraction than the high albedo state. Low albedo areas beyond this aridity value were considered as vulnerable regions where abrupt shifts in albedo may occur if aridity increases, as forecasted by current climate change models. Our findings indicate that remotely sensed albedo can identify two ecosystem states in African drylands. They support the suitability of albedo indices to inform us about discontinuous responses to aridity experienced by drylands, which can be linked to the onset of land degradation.This work was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant XDA19030500), the National Key Research and Development Program of China (Grant 2016YFC0503302), the European Research Council (BIODESERT project, ERC Grant Agreement 647038), the Joint PhD, Training Program of the University of Chinese Academy of Sciences, and the Research Foundation of Henan University of Technology (Grant 31401178)

Assessing a Multi-Platform Data Fusion Technique in Capturing Spatiotemporal Dynamics of Heterogeneous Dryland Ecosystems in Topographically Complex Terrain

Water-limited ecosystems encompass approximately 40% of terrestrial land mass and play a critical role in modulating Earth’s climate and provisioning ecosystem services to humanity. Spaceborne remote sensing is a critical tool for characterizing ecohydrologic patterns and advancing the understanding of the interactions between atmospheric forcings and ecohydrologic responses. Fine to medium scale spatial and temporal resolutions are needed to capture the spatial heterogeneity and the temporally intermittent response of these ecosystems to environmental forcings. Techniques combining complementary remote sensing datasets have been developed, but the heterogeneous nature of these regions present significant challenges. Here we investigate the capacity of one such approach, the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) algorithm, to map Normalized Difference Vegetation Index (NDVI) at 30 m spatial resolution and at a daily temporal resolution in an experimental watershed in southwest Idaho, USA. The Dry Creek Experimental Watershed captures an ecotone from a sagebrush steppe ecosystem to evergreen needle-leaf forests along an approximately 1000 m elevation gradient. We used STARFM to fuse NDVI retrievals from the MODerate-resolution Imaging Spectroradiometer (MODIS) and Landsat during the course of a growing season (April to September). Specifically we input to STARFM a pair of Landsat NDVI retrievals bracketing a sequence of daily MODIS NDVI retrievals to yield daily estimates of NDVI at resolutions of 30 m. In a suite of data denial experiments we compared these STARFM predictions against corresponding Landsat NDVI retrievals and characterized errors in predicted NDVI. We investigated how errors vary as a function of vegetation functional type and topographic aspect. We find that errors in predicting NDVI were highest during green-up and senescence and lowest during the middle of the growing season. Absolute errors were generally greatest in tree-covered portions of the watershed and lowest in locations characterized by grasses/bare ground. On average, relative errors in predicted average NDVI were greatest in grass/bare ground regions, on south-facing aspects, and at the height of the growing season. We present several ramifications revealed in this study for the use of multi-sensor remote sensing data for the study of spatiotemporal ecohydrologic patterns in dryland ecosystems

Monitoring paddy productivity in North Korea employing geostationary satellite images integrated with GRAMI-rice model

To meet the growing demands of staple crops with a strategy to develop amicable strategic measures that support efficient North Korean relief policies, it is a desirable task to accurately simulate the yield of paddy (Oryza sativa), an important Asian food commodity. We aim to address this with a gridbased crop simulation model integrated with satellite imagery that enables us to monitor the crop productivity of North Korea. Vegetation Indices (VIs), solar insolation, and air temperature data are thus obtained from the Communication Ocean and Meteorological Satellite (COMS), including the reanalysis data of the Korea Local Analysis and Prediction System (KLAPS). Paddy productivities for North Korea are projected based on the bidirectional reflectance distribution function-adjusted VIs and the solar insolation using the grid GRAMI-rice model. The model is calibrated on a 500-m grid paddy field in Cheorwon, and the model simulation performance accuracy is verified for Cheorwon and Paju, located at the borders of North Korea using four years of data from 2011 to 2014. Our results show that the paddy yields are reproduced reasonably accurately within a statistically significant range of accuracy, in comparison with observation data in Cheorwon (p = 0.183), Paju (p = 0.075), and NK (p = 0.101) according to a statistical t-test procedure. We advocate that incorporating a crop model with satellite images for crop yield simulations can be utilised as a reliable estimation technique for the monitoring of crop productivity, particularly in unapproachable, data-sparse regions not only in North Korea, but globally, where estimations of paddy productivity can assist in planning of agricultural activities that support regionally amicable food security strategies