White paper – On the use of LiDAR data at AmeriFlux sites

Our aim is to inform the AmeriFlux community on existing and upcoming LiDAR technologies (atmospheric Doppler or Raman LiDAR often deployed at flux sites are not considered here), how it is currently used at flux sites, and how we believe it could, in the future, further contribute to the AmeriFlux vision. Heterogeneity in vegetation and ground properties at various spatial scales is omnipresent at flux sites, and 3D mapping of canopy, understory, and ground surface can help move the science forward

Robust estimates of soil moisture and latent heat flux coupling strength obtained from triple collocation

Land surface models (LSMs) are often applied to predict the one-way coupling strength between surface soil moisture (SM) and latent heat (LH) flux. However, the ability of LSMs to accurately represent such coupling has not been adequately established. Likewise, the estimation of SM/LH coupling strength using ground-based observational data is potentially compromised by the impact of independent SM and LH measurements errors. Here we apply a new statistical technique to acquire estimates of one-way SM/LH coupling strength which are nonbiased in the presence of random error using a triple collocation approach based on leveraging the simultaneous availability of independent SM and LH estimates acquired from (1) LSMs, (2) satellite remote sensing, and (3) ground-based observations. Results suggest that LSMs do not generally overestimate the strength of one-way surface SM/LH coupling

Estimation of surface turbulent heat fluxes via variational assimilation of sequences of land surface temperatures from Geostationary Operational Environmental Satellites

Recently, a number of studies have focused on estimating surface turbulent heat fluxes via assimilation of sequences of land surface temperature (LST) observations into variational data assimilation (VDA) schemes. Using the full heat diffusion equation as a constraint, the surface energy balance equation can be solved via assimilation of sequences of LST within a VDA framework. However, the VDA methods have been tested only in limited field sites that span only a few climate and land use types. Hence, in this study, combined-source (CS) and dual-source (DS) VDA schemes are tested extensively over six FluxNet sites with different vegetation covers (grassland, cropland, and forest) and climate conditions. The CS model groups the soil and canopy together as a single source and does not consider their different contributions to the total turbulent heat fluxes, while the DS model considers them to be different sources. LST data retrieved from the Geostationary Operational Environmental Satellites are assimilated into these two VDA schemes. Sensible and latent heat flux estimates from the CS and DS models are compared with the corresponding measurements from flux tower stations. The results indicate that the performance of both models at dry, lightly vegetated sites is better than that at wet, densely vegetated sites. Additionally, the DS model outperforms the CS model at all sites, implying that the DS scheme is more reliable and can characterize the underlying physics of the problem better

A remote sensing and modeling integrated approach for constructing continuous time series of daily actual evapotranspiration

Satellite remote sensing-based surface energy balance (SEB) techniques have emerged as useful tools for quantifying spatialized actual evapotranspiration at various temporal and spatial scales. However, discontinuous data acquisitions and/or gaps in image acquisition due to cloud cover can limit the applicability of satellite remote sensing (RS) in agriculture water management where continuous time series of daily crop actual evapotranspiration (ETc act) are more valued. The aim of the research is to construct continuous time series of daily ETc act starting from temporal estimates of actual evapotranspiration obtained by SEB modelling (ETa eb) on Landsat-TM images. SEBAL model was integrated with the FAO 56 evaporation model, RS-retrieved vegetative biomass dynamics (by NDVI) and on-field measurements of soil moisture and potential evapotranspiration. The procedure was validated by an eddy covariance tower on a vineyard with partial soil coverage in the south of Sardinia Island, Italy. The integrated modeling approach showed a good reproduction of the time series dynamics of observed ETc act (R2 =0.71, MAE=0.54 mm d-1, RMSE=0.73 mm d-1). A daily and a cumulative monthly temporal analysis showed the importance of integrating parameters that capture changes in the soil-plant-atmosphere (SPA) continuum between Landsat acquisitions. The comparison with daily ETc act obtained by the referenced ET fraction (ETrF) method that considers only weather variability (by ETo) confirmed the lead of the proposed procedure in the spring/early summer periods when vegetation biomass changes and soil water evaporation have a significant weight in the ET process. The applied modelling approach was also robust in constructing the missing ETc act data under scenarios of limited cloud-free Landsat acquisitions. The presented integrated approach has a great potential for the near real time monitoring and scheduling of irrigation practices. Further testing of this approach with diverse dataset and the integration with the soil water modeling is to be analyzed in future work

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

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

Evapotranspiration Estimates Derived Using Multi-Platform Remote Sensing in a Semiarid Region

Evapotranspiration (ET) is a key component of the water balance, especially in arid and semiarid regions. The current study takes advantage of spatially-distributed, near real-time information provided by satellite remote sensing to develop a regional scale ET product derived from remotely-sensed observations. ET is calculated by scaling PET estimated from Moderate Resolution Imaging Spectroradiometer (MODIS) products with downscaled soil moisture derived using the Soil Moisture Ocean Salinity (SMOS) satellite and a second order polynomial regression formula. The MODis-Soil Moisture ET (MOD-SMET) estimates are validated using four flux tower sites in southern Arizona USA, a calibrated empirical ET model, and model output from Version 2 of the North American Land Data Assimilation System (NLDAS-2). Validation against daily eddy covariance ET indicates correlations between 0.63 and 0.83 and root mean square errors (RMSE) between 40 and 96 W/m2. MOD-SMET estimates compare well to the calibrated empirical ET model, with a −0.14 difference in correlation between sites, on average. By comparison, NLDAS-2 models underestimate daily ET compared to both flux towers and MOD-SMET estimates. Our analysis shows the MOD-SMET approach to be effective for estimating ET. Because it requires limited ancillary ground-based data and no site-specific calibration, the method is applicable to regions where ground-based measurements are not available

Remote Sensing Monitoring of Land Surface Temperature (LST)

This book is a collection of recent developments, methodologies, calibration and validation techniques, and applications of thermal remote sensing data and derived products from UAV-based, aerial, and satellite remote sensing. A set of 15 papers written by a total of 70 authors was selected for this book. The published papers cover a wide range of topics, which can be classified in five groups: algorithms, calibration and validation techniques, improvements in long-term consistency in satellite LST, downscaling of LST, and LST applications and land surface emissivity research

An Intercomparison of Satellite-Based Daily Evapotranspiration Estimates under Different Eco-Climatic Regions in South Africa

Knowledge of evapotranspiration (ET) is essential for enhancing our understanding of the hydrological cycle, as well as for managing water resources, particularly in semi-arid regions. Remote sensing offers a comprehensive means of monitoring this phenomenon at different spatial and temporal intervals. Currently, several satellite methods exist and are used to assess ET at various spatial and temporal resolutions with various degrees of accuracy and precision. This research investigated the performance of three satellite-based ET algorithms and two global products, namely land surface temperature/vegetation index (TsVI), Penman–Monteith (PM), and the Meteosat Second Generation ET (MET) and the Global Land-surface Evaporation: the Amsterdam Methodology (GLEAM) global products, in two eco-regions of South Africa. Daily ET derived from the eddy covariance system from Skukuza, a sub-tropical, savanna biome, and large aperture boundary layer scintillometer system in Elandsberg, a Mediterranean, fynbos biome, during the dry and wet seasons, were used to evaluate the models. Low coefficients of determination (R2) of between 0 and 0.45 were recorded on both sites, during both seasons. Although PM performed best during periods of high ET at both sites, results show it was outperformed by other models during low ET times. TsVI and MET were similarly accurate in the dry season in Skukuza, as GLEAM was the most accurate in Elandsberg during the wet season. The conclusion is that none of the models performed well, as shown by low R2 and high errors in all the models. In essence, our results conclude that further investigation of the PM model is possible to improve its estimation of low ET measurements

Mapping evapotranspiration with high-resolution aircraft imagery over vineyards using one- and two-source modeling schemes

Thermal and multispectral remote sensing data from low-altitude aircraft can provide high spatial resolution necessary for sub-field ( 10 m) and plant canopy (1 m) scale evapotranspiration (ET) monitoring. In this study, highresolution (sub-meter-scale) thermal infrared and multispectral shortwave data from aircraft are used to map ET over vineyards in central California with the two-source energy balance (TSEB) model and with a simple model having operational immediate capabilities called DATTUTDUT (Deriving Atmosphere Turbulent Transport Useful To Dummies Using Temperature). The latter uses contextual information within the image to scale between radiometric land surface temperature (TR) values representing hydrologic limits of potential ET and a non-evaporative surface. Imagery from 5 days throughout the growing season is used for mapping ET at the sub-field scale. The performance of the two models is evaluated using tower-based measurements of sensible (H) and latent heat (LE) flux or ET. The comparison indicates that TSEB was able to derive reasonable ET estimates under varying conditions, likely due to the physically based treatment of the energy and the surface temperature partitioning between the soil/cover crop inter-row and vine canopy elements. On the other hand, DATTUTDUT performance was somewhat degraded presumably because the simple scaling scheme does not consider differences in the two sources (vine and inter-row) of heat and temperature contributions or the effect of surface roughness on the efficiency of heat exchange. Maps of the evaporative fraction (EFDLE/(H CLE)) from the two models had similar spatial patterns but different magnitudes in some areas within the fields on certain days. Large EF discrepancies between the models were found on 2 of the 5 days (DOY 162 and 219) when there were significant differences with the tower-based ET measurements, particularly using the DATTUTDUT model. These differences in EF between the models translate to significant variations in daily water use estimates for these 2 days for the vineyards. Model sensitivity analysis demonstrated the high degree of sensitivity of the TSEB model to the accuracy of the TR data, while the DATTUTDUT model was insensitive to systematic errors in TR as is the case with contextual-based models. However, it is shown that the study domain and spatial resolution will significantly influence the ET estimation from the DATTUTDUT model. Future work is planned for developing a hybrid approach that leverages the strengths of both modeling schemes and is simple enough to be used operationally with high-resolution imagery

Remote Sensing of Environment: Current status of Landsat program, science, and applications

Formal planning and development of what became the first Landsat satellite commenced over 50 years ago in 1967. Now, having collected earth observation data for well over four decades since the 1972 launch of Landsat- 1, the Landsat program is increasingly complex and vibrant. Critical programmatic elements are ensuring the continuity of high quality measurements for scientific and operational investigations, including ground systems, acquisition planning, data archiving and management, and provision of analysis ready data products. Free and open access to archival and new imagery has resulted in a myriad of innovative applications and novel scientific insights. The planning of future compatible satellites in the Landsat series, which maintain continuity while incorporating technological advancements, has resulted in an increased operational use of Landsat data. Governments and international agencies, among others, can now build an expectation of Landsat data into a given operational data stream. International programs and conventions (e.g., deforestation monitoring, climate change mitigation) are empowered by access to systematically collected and calibrated data with expected future continuity further contributing to the existing multi-decadal record. The increased breadth and depth of Landsat science and applications have accelerated following the launch of Landsat-8, with significant improvements in data quality. Herein, we describe the programmatic developments and institutional context for the Landsat program and the unique ability of Landsat to meet the needs of national and international programs. We then present the key trends in Landsat science that underpin many of the recent scientific and application developments and followup with more detailed thematically organized summaries. The historical context offered by archival imagery combined with new imagery allows for the development of time series algorithms that can produce information on trends and dynamics. Landsat-8 has figured prominently in these recent developments, as has the improved understanding and calibration of historical data. Following the communication of the state of Landsat science, an outlook for future launches and envisioned programmatic developments are presented. Increased linkages between satellite programs are also made possible through an expectation of future mission continuity, such as developing a virtual constellation with Sentinel-2. Successful science and applications developments create a positive feedback loop—justifying and encouraging current and future programmatic support for Landsat