Long time series (1984–2020) of albedo variations on the Greenland ice sheet from harmonized Landsat and Sentinel 2 imagery

Albedo is a key factor in modulating the absorption of solar radiation on ice surfaces. Satellite measurements have shown a general reduction in albedo across the Greenland ice sheet over the past few decades, particularly along the western margin of the ice sheet, a region known as the Dark Zone (albedo < 0.45). Here we chose a combination of Landsat 4–8 and Sentinel 2 imagery to enable us to derive the longest record of albedo variations in the Dark Zone, running from 1984 to 2020. We developed a simple, pragmatic and efficient sensor transformation to provide a long time series of consistent, harmonized satellite imagery. Narrow to broadband conversion algorithms were developed from regression models of harmonized satellite data and in situ albedo from the Program for Monitoring of the Greenland Ice Sheet (PROMICE) automatic weather stations. The albedo derived from the harmonized Landsat and Sentinel 2 data shows that the maximum extent of the Dark Zone expanded rapidly between 2005 and 2007, increasing to ~280% of the average annual maximum extent of 2900 km2 to ~8000 km2 since. The Dark Zone is continuing to darken slowly, with the average annual minimum albedo decreasing at a rate of (p = 0.16, 2001–2020)

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

Cross-comparison of albedo products for glacier surfaces derived from airborne and satellite (sentinel-2 and landsat 8) optical data

Surface albedo partitions the amount of energy received by glacier surfaces from shortwave fluxes and modulates the energy available for melt processes. The ice- albedo feedback, influenced by the contamination of bare-ice surfaces with light- absorbing impurities, plays a major role in the melting of mountain glaciers in a warming climate. However, little is known about the spatial and temporal distribution and variability of bare-ice glacier surface albedo under changing conditions. In this study, we focus on two mountain glaciers located in the western Swiss Alps and perform a cross-comparison of different albedo products. We take advantage of high spectral and spatial resolution (284 bands, 2 m) imaging spectrometer data from the Airborne Prism Experiment (APEX) and investigate the applicability and potential of Sentinel-2 and Landsat 8 data to derive broadband albedo products. The performance of shortwave broadband albedo retrievals is tested and we assess the reliability of published narrow-to-broadband conversion algorithms. The resulting albedo products from the three sensors and different algorithms are further cross-compared. Moreover, the impact of the anisotropy correction is analysed depending on different surface types. While degradation of the spectral resolution impacted glacier-wide mean albedo by about 5%, reducing the spatial resolution resulted in changes of less than 1%. However, in any case, coarser spatial resolution was no longer able to represent small-scale variability of albedo on glacier surfaces. We discuss the implications when using Sentinel-2 and Landsat 8 to map dynamic glaciological processes and to monitor glacier surface albedo on larger spatial and more frequent temporal scales

Measurement and Modeling of Ground-Level Ozone Concentration in Catania, Italy using Biophysical Remote Sensing and GIS

This experimental study examined spatial variation of ground level ozone (O3) in the city of Catania, Italy using thirty passive samplers deployed in a 500-m grid pattern. Significant spatial variation in ground level O3 concentrations (ranging from 12.8 to 41.7 g/m3) was detected across Catania’s urban core and periphery. Biophysical measures derived from satellite imagery and built environment characteristics from GIS were evaluated as correlates of O3 concentrations. A land use regression model based on four variables (land surface temperature, building area, residential street length, and distance to the coast) explained 74% of the variance (adjusted R2) in measured O3. The results of the study suggest that biophysical remote sensing variables are worth further investigation as predictors of ground level O3 (and potentially other air pollutants) because they provide objective measurements that can be tested across multiple locations and over time

Earth Observations for Addressing Global Challenges

"Earth Observations for Addressing Global Challenges" presents the results of cutting-edge research related to innovative techniques and approaches based on satellite remote sensing data, the acquisition of earth observations, and their applications in the contemporary practice of sustainable development. Addressing the urgent tasks of adaptation to climate change is one of the biggest global challenges for humanity. As His Excellency António Guterres, Secretary-General of the United Nations, said, "Climate change is the defining issue of our time—and we are at a defining moment. We face a direct existential threat." For many years, scientists from around the world have been conducting research on earth observations collecting vital data about the state of the earth environment. Evidence of the rapidly changing climate is alarming: according to the World Meteorological Organization, the past two decades included 18 of the warmest years since 1850, when records began. Thus, Group on Earth Observations (GEO) has launched initiatives across multiple societal benefit areas (agriculture, biodiversity, climate, disasters, ecosystems, energy, health, water, and weather), such as the Global Forest Observations Initiative, the GEO Carbon and GHG Initiative, the GEO Biodiversity Observation Network, and the GEO Blue Planet, among others. The results of research that addressed strategic priorities of these important initiatives are presented in the monograph

Chapter Earth Observation for Urban Climate Monitoring: Surface Cover and Land Surface Temperature

The rate at which global climate change is happening is arguably the most pressing environmental challenge of the century, and it affects our cities. Climate change exerts added stress on urban areas through increased numbers of heat waves threatening people’s well-being and, in many cases, human lives. Earth observation (EO) systems and the advances in remote sensing technology increase the opportunities for monitoring the thermal behavior of cities. The Sentinels constitute the first series of operational satellites for Copernicus, a program launched to provide data, information, services, and knowledge in support of Europe’s goals regarding sustainable development and global governance of the environment. This chapter examines the exploitation of EO data for monitoring the urban climate, with particular focus on the urban surface cover and temperature. Two example applications are analyzed: the mapping of the urban surface and its characteristics, using EO data and the estimation of urban temperatures. Approaches, like the ones described in this chapter, can become operational once adapted to Sentinels, since their long-term operation plan guarantees the future supply of satellite observations. Thus, the described methods may support planning activities related to climate change mitigation and adaptation in cities, as well as routine urban planning activities

Comparison of layer-stacking and Dempster-Shafer theory-based methods using Sentinel-1 and Sentinel-2 data fusion in urban land cover mapping

Data fusion has shown potential to improve the accuracy of land cover mapping, and selection of the optimal fusion technique remains a challenge. This study investigated the performance of fusing Sentinel-1 (S-1) and Sentinel-2 (S-2) data, using layer-stacking method at the pixel level and Dempster-Shafer (D-S) theory-based approach at the decision level, for mapping six land cover classes in Thu Dau Mot City, Vietnam. At the pixel level, S-1 and S-2 bands and their extracted textures and indices were stacked into the different single-sensor and multi-sensor datasets (i.e. fused datasets). The datasets were categorized into two groups. One group included the datasets containing only spectral and backscattering bands, and the other group included the datasets consisting of these bands and their extracted features. The random forest (RF) classifier was then applied to the datasets within each group. At the decision level, the RF classification outputs of the single-sensor datasets within each group were fused together based on D-S theory. Finally, the accuracy of the mapping results at both levels within each group was compared. The results showed that fusion at the decision level provided the best mapping accuracy compared to the results from other products within each group. The highest overall accuracy (OA) and Kappa coefficient of the map using D-S theory were 92.67% and 0.91, respectively. The decision-level fusion helped increase the OA of the map by 0.75% to 2.07% compared to that of corresponding S-2 products in the groups. Meanwhile, the data fusion at the pixel level delivered the mapping results, which yielded an OA of 4.88% to 6.58% lower than that of corresponding S-2 products in the groups

Multiscale Imaging of Evapotranspiration

Evapotranspiration (ET; evaporation + transpiration) is central to a wide range of biological, chemical, and physical processes in the Earth system. Accurate remote sensing of ET is challenging due to the interrelated and generally scale dependent nature of the physical factors which contribute to the process. The evaporation of water from porous media like sands and soils is an important subset of the complete ET problem. Chapter 1 presents a laboratory investigation into this question, examining the effects of grain size and composition on the evolution of drying sands. The effects of composition are found to be 2-5x greater than the effects of grain size, indicating that differences in heating caused by differences in reflectance may dominate hydrologic differences caused by grain size variation. In order to relate the results of Chapter 1 to the satellite image archive, however, the question of information loss between hyperspectral (measurements at 100s of wavelength intervals) laboratory measurements and multispectral (≤ 12 wavelength intervals) satellite images must be addressed. Chapter 2 focuses on this question as applied to substrate materials such as sediment, soil, rock, and non-photosynthetic vegetation. The results indicate that the continuum that is resolved by multispectral sensors is sufficient to resolve the gradient between sand-rich and clay-rich soils, and that this gradient is also a dominant feature in hyperspectral mixing spaces where the actual absorptions can be resolved. Multispectral measurements can be converted to biogeophysically relevant quantities using spectral mixture analysis (SMA). However, retrospective multitemporal analysis first requires cross-sensor calibration of the mixture model. Chapter 3 presents this calibration, allowing multispectral image data to be used interchangeably throughout the Landsat 4-8 archive. In addition, a theoretical explanation is advanced for the observed superior scaling properties of SMA-derived fraction images over spectral indices. The physical quantities estimated by the spectral mixture model are then compared to simultaneously imaged surface temperature, as well as to the derived parameters of ET Fraction and Moisture Availability. SMA-derived vegetation abundance is found to produce substantially more informative ET maps, and SMA-derived substrate fraction is found to yield a surprisingly strong linear relationship with surface temperature. These results provide context for agricultural applications. Chapter 5 investigates the question of mapping and monitoring rice agricultural using optical and thermal satellite image time series. Thermal image time series are found to produce more accurate maps of rice presence/absence, but optical image time series are found to produce more accurate maps of rice crop timing. Chapter 6 takes a more global approach, investigating the spatial structure of agricultural networks for a diverse set of landscapes. Surprisingly consistent scaling relations are found. These relations are assessed in the context of a network-based approach to land cover analysis, with potential implications for the scale dependence of ET estimates. In sum, this thesis present a novel approach to improving ET estimation based on a synthesis of complementary laboratory measurements, satellite image analysis, and field observations. Alone, each of these independent sources of information provides novel insights. Viewed together, these insights form the basis of a more accurate and complete geophysical understanding of the ET phenomenon

Bayesian atmospheric correction over land: Sentinel-2/MSI and Landsat 8/OLI

Mitigating the impact of atmospheric effects on optical remote sensing data is critical for monitoring intrinsic land processes and developing Analysis Ready Data (ARD). This work develops an approach to this for the NERC NCEO medium resolution ARD Landsat 8 (L8) and Sentinel 2 (S2) products, called Sensor Invariant Atmospheric Correction (SIAC). The contribution of the work is to phrase and solve that problem within a probabilistic (Bayesian) framework for medium resolution multispectral sensors S2/MSI and L8/OLI and to provide per-pixel uncertainty estimates traceable from assumed top-of-atmosphere (TOA) measurement uncertainty, making progress towards an important aspect of CEOS ARD target requirements. A set of observational and a priori constraints are developed in SIAC to constrain an estimate of coarse resolution (500 m) aerosol optical thickness (AOT) and total column water vapour (TCWV), along with associated uncertainty. This is then used to estimate the medium resolution (10–60 m) surface reflectance and uncertainty, given an assumed uncertainty of 5 % in TOA reflectance. The coarse resolution a priori constraints used are the MODIS MCD43 BRDF/Albedo product, giving a constraint on 500 m surface reflectance, and the Copernicus Atmosphere Monitoring Service (CAMS) operational forecasts of AOT and TCWV, providing estimates of atmospheric state at core 40 km spatial resolution, with an associated 500 m resolution spatial correlation model. The mapping in spatial scale between medium resolution observations and the coarser resolution constraints is achieved using a calibrated effective point spread function for MCD43. Efficient approximations (emulators) to the outputs of the 6S atmospheric radiative transfer code are used to estimate the state parameters in the atmospheric correction stage. SIAC is demonstrated for a set of global S2 and L8 images covering AERONET and RadCalNet sites. AOT retrievals show a very high correlation to AERONET estimates (correlation coefficient around 0.86, RMSE of 0.07 for both sensors), although with a small bias in AOT. TCWV is accurately retrieved from both sensors (correlation coefficient over 0.96, RMSE &lt;0.32 g cm−2). Comparisons with in situ surface reflectance measurements from the RadCalNet network show that SIAC provides accurate estimates of surface reflectance across the entire spectrum, with RMSE mismatches with the reference data between 0.01 and 0.02 in units of reflectance for both S2 and L8. For near-simultaneous S2 and L8 acquisitions, there is a very tight relationship (correlation coefficient over 0.95 for all common bands) between surface reflectance from both sensors, with negligible biases. Uncertainty estimates are assessed through discrepancy analysis and are found to provide viable estimates for AOT and TCWV. For surface reflectance, they give conservative estimates of uncertainty, suggesting that a lower estimate of TOA reflectance uncertainty might be appropriate

A collaborative change detection approach on multi-sensor spatial imagery for desertwetland monitoring after a flash flood in Southern Morocco

© 2019 by the authors. This study aims to present a technique that combines multi-sensor spatial data to monitor wetland areas after a flash-flood event in a Saharan arid region. To extract the most efficient information, seven satellite images (radar and optical) taken before and after the event were used. To achieve the objectives, this study used Sentinel-1 data to discriminate water body and soil roughness, and optical data to monitor the soil moisture after the event. The proposed method combines two approaches: one based on spectral processing, and the other based on categorical processing. The first step was to extract four spectral indices and utilize change vector analysis on multispectral diachronic images from three MSI Sentinel-2 images and two Landsat-8 OLI images acquired before and after the event. The second step was performed using pattern classification techniques, namely, linear classifiers based on support vector machines (SVM) with Gaussian kernels. The results of these two approaches were fused to generate a collaborative wetland change map. The application of co-registration and supervised classification based on textural and intensity information from Radar Sentinel-1 images taken before and after the event completes this work. The results obtained demonstrate the importance of the complementarity of multi-sensor images and a multi-approach methodology to better monitor changes to a wetland area after a flash-flood disaster