Characterizing Land Surface Anisotropic Reflectance over Rugged Terrain: A Review of Concepts and Recent Developments

Vegetation browning: global drivers, impacts, and feedbacks

As global climate conditions continue to change, disturbance regimes and environmental drivers will continue to shift, impacting global vegetation dynamics. Following a period of vegetation greening, there has been a progressive increase in remotely sensed vegetation browning globally. Given the many societal benefits that forests provide, it is critical that we understand vegetation dynamic alterations. Here, we review associative drivers, impacts, and feedbacks, revealing the complexity of browning. Concomitant increases in browning include the weakening of ecosystem services and functions and alterations to vegetation structure and species composition, as well as the development of potential positive climate change feedbacks. Also discussed are the current challenges in browning detection and understanding associated impacts and feedbacks. Finally, we outline recommended strategies.

Full-band, multi-angle, multi-scale, and temporal dynamic field spectral measurements in China

Field-measured spectra are critical for remote sensing physical modelling, retrieval of structural, biophysical, and biochemical parameters, and other practical applications. We present a library of field spectra, which includes (1) portable field spectroradiometer measurements of vegetation, soil, and snow in the full-wave band, (2) multi-angle spectra measurements of desert vegetation, chernozems, and snow with consideration of the anisotropic reflectance of land surface, (3) multi-scale spectra measurements of leaf and canopy of different vegetation cover surfaces, and (4) continuous reflectance spectra time-series data revealing vegetation growth dynamics of maize, rice, wheat, rape, grassland, and so on. To the best of our knowledge, this library is unique in simultaneously providing full-band, multi-angle, multi-scale spectral measurements of the main surface elements of China covering a large spatial extent over a 10-year period. Furthermore, the 101 by 101 satellite pixels of Landsat ETM/OLI and MODIS surface reflectance centered around the field site were extracted, providing a vital linkage between ground measurements and satellite observations. The code language used for this work is Matlab 2016a.

Analysis of the relationship among land surface temperature (LST), land use land cover (LULC), and normalized difference vegetation index (NDVI) with topographic elements in the lower Himalayan region

Land Surface Temperature (LST) affects exchange of energy between earth surface and atmosphere which is important for studying environmental changes. However, research on the relationship between LST, Land Use Land Cover (LULC), and Normalized Difference Vegetation Index (NDVI) with topographic elements in the lower Himalayan region has not been done. Therefore, the present study explored the relationship between LST and NDVI, and LULC types with topographic elements in the lower Himalayan region of Pakistan. The study area was divided into North-South, West-East, North-West to South-East and North-East to South-East directions using ArcMap 3D analysis. The current study used Landsat 8 (OLI/TIRS) data from May 2021 for LULC and LST analysis in the study area. The LST data was obtained from the thermal band of Landsat 8 (TIRS), while the LULC of the study areas was classified using the Maximum Likelihood Classification (MLC) method utilizing Landsat 8 (OLI) data. TIRS collects data for two narrow spectral bands (B10 and B11) with spectral wavelength of 10.6 μm-12.51 μm in the thermal region formerly covered by one wide spectral band (B6) on Landsat 4-7. With 12-bit data products, TIRS data is available in radiometric, geometric, and terrain-corrected file format. The effect of elevation on LST was assessed using LST and elevation data obtained from the USGS website. The LST across LULC types with sunny and shady slopes was analyzed to assess the influence of slope directions. The relationship of LST with elevation and NDVI was examined using correlation analysis. The results indicated that LST decreased from North-South and South-East, while increasing from North-East and South-West directions. The correlation coefficient between LST and elevation was negative, with an R-value of -0.51. The NDVI findings with elevation showed that NDVI increases with an increase in elevation. Zonal analysis of LST for different LULC types showed that built-up and bare soil had the highest mean LST, which was 35.76 °C and 28.08 °C, respectively, followed by agriculture, vegetation, and water bodies. The mean LST difference between sunny and shady slopes was 1.02 °C. The correlation between NDVI and LST was negative for all LULC types except the water body. This study findings can be used to ensure sustainable urban development and minimize urban heat island effects by providing effective guidelines for urban planners, policymakers, and respective authorities in the Lower Himalayan region. The current thermal remote sensing findings can be used to model energy fluxes and surface processes in the study area.

High Spatiotemporal Rugged Land Surface Temperature Downscaling over Saihanba Forest Park, China

Satellite-derived rugged land surface temperature (LST) is an important parameter indicating the status of the Earth’s surface energy budget and its seasonal/temporal dynamic change. However, existing LST products from rugged areas are more prone to error when supporting applications in mountainous areas and Earth surface processes that occur at high spatial and temporal resolutions. This research aimed to develop a method for generating rugged LST with a high temporal and spatial resolution by using an improved ensemble LST model combining three regressors, including a random forest, a ridge, and a support vector machine. Different combinations of high-resolution input parameters were also considered in this study. The input datasets included Moderate Resolution Imaging Spectroradiometer (MODIS) LST datasets (MxD11A1) for nighttime, temporal Sentinel-2 Multispectral Instrument (MSI) datasets, and digital elevation model (DEM) datasets. The 30 m rugged LST datasets derived were compared against an in situ LST dataset obtained at Saihanba Forest Park (SFP) sites and an ASTER-derived 90 m LST, respectively. The results with in situ measurements demonstrated significant LST details, with an R2 higher than 0.95 and RMSE around 3.00 K for both Terra/MOD- and Aqua/MYD-based LST datasets, and with slightly better results being obtained from the Aqua/MYD-based LST than that from Terra/MOD. The inter-comparison results with ASTER LST showed that over 80% of the pixels of the difference image for the two datasets were within 2 K. In light of the complex topography and distinct atmospheric conditions, these comparison results are encouraging. The 30 m LST from the method proposed in this study also depicts the seasonality of rugged surfaces.

SGOT: A Simplified Geometric-Optical Model for Crown Scene Components Modeling over Rugged Terrain

Topography affects the fraction of scene components of the canopy and background, resulting in the observed reflectance distortion. Modeling the canopy reflectance over rugged terrain needs to account for topographic effects. For this purpose, the existing models greatly increased the mathematical complexity while improving description of terrain and crown structure, which dramatically decreased the computational efficiency so as to limit their universal application. In this study, we developed a simplified geometric-optical model (SGOT) for simulating the scene components over rugged terrain. The geotropism of tree growth was considered to make SGOT physically sound. The internal structure of crown was simplified to make SGOT mathematically simpler. Scene component observations derived from Persistence of Vision Ray-tracer (POV-Ray) on surfaces with different normal directions and simulations were made using Geometric-Optical and Mutual Shadowing Coupled with Topography Model (GOMST) and Geometric-Optical for Sloping Terrains Model GOST; models were combined to test the SGOT model. In addition, topographic factors and crown density effect on the scene components modeling were analyzed. The results indicated that SGOT has good accuracy (R2 for the areal proportions of sunlit crown (Kc), sunlit background (Kg), shaded crown (Kt), and shaded background (Kz) are 0.853, 0.857, 0.914, and 0.838, respectively) compared with POV-Ray simulation, and performs better than GOMST, especially in scenes with high crown density. Moreover, SGOT outperformed the compared models in computational efficiency (4% faster than GOMST and 29.5% faster than GOST). Finally, the simulations of the scene components distribution in different topographic factors and crown density were further discussed. SGOT and GOST can both capture scene component variations caused by terrain better than GOMST, but comparatively, SGOT provides a more efficient tool to simulate the crown scene components because of its physical soundness and mathematical simplicity, and consequently, it will facilitate the modeling of canopy reflectance over mountainous regions.

Revisiting the Performance of the Kernel-Driven BRDF Model Using Filtered High-Quality POLDER Observations

The Bidirectional Reflectance Distribution Function (BRDF) is usually used to describe the reflectance anisotropy of a non-Lambertian surface and estimate surface parameters. Among the BRDF models, the kernel-driven models have been extensively used due to their simple form and powerful fitting ability, and their reliability has been validated in some studies. However, existing validation efforts used in situ measurements or limited satellite data, which may be subject to inadequate observational conditions or quality uncertainties. A recently released high-quality BRDF database from Polarization and Directionality of the Earth’s Reflectances (POLDER) provides an opportunity to revisit the performance of the kernel-driven models. Therefore, in order to evaluate the fitting ability of the kernel-driven models under different observational conditions and explore their application direction in the future, we use the filtered high-quality BRDF database to evaluate the fitting ability of the kernel-driven model represented by the RossThick-LiSparseR (RTLSR) kernels in this paper. The results show that the RTLSR model performs well, which shows small fitting residuals under most observational conditions. However, the applicability of the RTLSR model performed differently across land cover types; the RTLSR model exhibited larger fitting residuals, especially over non-vegetated surfaces. Under different sun-sensor geometries, the fitting residuals show a strong positive correlation with the Solar Zenith Angle. The above two factors cause the RTLSR model to exhibit a poorer fitting ability at high latitudes. As an exploration, we designed a model combination strategy that combines the advantages of different models and achieved a better performance at high latitudes. We believe that this study provides a better understanding of the RTLSR model.

Extending the GOSAILT Model to Simulate Sparse Woodland Bi-Directional Reflectance with Soil Reflectance Anisotropy Consideration

Anisotropic canopy reflectance plays a crucial role in estimating vegetation biophysical parameters, whereas soil reflectance anisotropy affects canopy reflectance. However, woodland canopy bidirectional reflectance distribution function (BRDF) models considering soil anisotropy are far from universal, especially for the BRDF models of mountain forest. In this study, a mountain forest canopy model, named geometric-optical and mutual shadowing and scattering from arbitrarily inclined-leaves model coupled with topography (GOSAILT), was extended to consider the soil anisotropic reflectance characteristics by introducing the simple soil directional (SSD) reflectance model. The modified GOSAILT model (named GOSAILT-SSD) was evaluated using unmanned aerial vehicle (UAV) field observations and discrete anisotropic radiative transfer (DART) simulations. Then, the effects of Lambertian soil assumption on simulating the vi-directional reflectance factor (BRF) were evaluated across different fractions of vegetation cover (Cv), view zenith angles (VZA), solar zenith angles (SZA), and spectral bands with the GOSAILT-SSD model. The evaluation results, with the DART simulations, show that the performance of the GOSAILT-SSD model in simulating canopy BRF is significantly improved, with decreasing RMSE, from 0.027 to 0.017 for the red band and 0.051 to 0.037 for the near-infrared (NIR) band. Meanwhile, the GOSAILT-SSD simulations show high consistency with UAV multi-angular observations (R2 = 0.97). Besides, it is also found that the BRF simulation errors caused by Lambertian soil assumption are too large to be neglected, with a maximum relative bias of about 45% for the red band. This inappropriate assumption results in a remarkable BRF underestimation near the hot spot direction and an obvious BRF overestimation for large VZA in the solar principal plane (PP). Meanwhile, this simulation bias decreases with the increase of fraction of vegetation cover. This study provides an effective technique to improve the capability of the mountain forest canopy BRDF model by considering the soil anisotropic characteristics for advancing the modeling of radiative transfer (RT) processes over rugged terrain.

Comparing Three Remotely Sensed Approaches for Simulating Gross Primary Productivity over Mountainous Watersheds: A Case Study in the Wanglang National Nature Reserve, China

Light Use Efficiency (LUE), Vegetation Index (VI)-based, and process-based models are the main approaches for spatially continuous gross primary productivity (GPP) estimation. However, most current GPP models overlook the effects of topography on the vegetation photosynthesis process. Based on the structures of a two-leaf LUE model (TL-LUE), a VI-based model (temperature and greenness, TG), and a process-based model (Boreal Ecosystem Productivity Simulator, BEPS), three models, named mountain TL-LUE (MTL-LUE), mountain TG (MTG), and BEPS-TerrainLab, have been proposed to improve GPP estimation over mountainous areas. The GPP estimates from the three mountain models have been proven to align more closely with tower-based GPP than those from the original models at the site scale, but their abilities to characterize the spatial variation of GPP at the watershed scale are not yet known. In this work, the GPP estimates from three LUE models (i.e., MOD17, TL-LUE, and MTL-LUE), two VI-based models (i.e., TG and MTG), and two process-based models (i.e., BEPS and BEPS-TerrainLab) were compared for a mountainous watershed. At the watershed scale, the annual GPP estimates from MTL-LUE, MTG, and BTL were found to have a higher spatial variation than those from the original models (increasing the spatial coefficient of variation by 6%, 8%, and 22%), highlighting that incorporating topographic information into GPP models might improve understanding of the high spatial heterogeneity of the vegetation photosynthesis process over mountainous areas. Obvious discrepancies were also observed in the GPP estimates from MTL-LUE, MTG, and BTL, with determination coefficients ranging from 0.02–0.29 and root mean square errors ranging from 399–821 gC m−2yr−1. These GPP discrepancies mainly stem from the different (1) structures of original LUE, VI, and process models, (2) assumptions associated with the effects of topography on photosynthesis, (3) input data, and (4) values of sensitive parameters. Our study highlights the importance of considering surface topography when modeling GPP over mountainous areas, and suggests that more attention should be given to the discrepancy of GPP estimates from different models.

Effect of the Solar Zenith Angles at Different Latitudes on Estimated Crop Vegetation Indices

Normalization of anisotropic solar reflectance is an essential factor that needs to be considered for field-based phenotyping applications to ensure reliability, consistency, and interpretability of time-series multispectral data acquired using an unmanned aerial vehicle (UAV). Different models have been developed to characterize the bidirectional reflectance distribution function. However, the substantial variation in crop breeding trials, in terms of vegetation structure configuration, creates challenges to such modeling approaches. This study evaluated the variation in standard vegetation indices and its relationship with ground-reference data (measured crop traits such as seed/grain yield) in multiple crop breeding trials as a function of solar zenith angles (SZA). UAV-based multispectral images were acquired and utilized to extract vegetation indices at SZA across two different latitudes. The pea and chickpea breeding materials were evaluated in a high latitude (46°36′39.92″ N) zone, whereas the rice lines were assessed in a low latitude (3°29′42.43″ N) zone. In general, several of the vegetation index data were affected by SZA (e.g., normalized difference vegetation index, green normalized difference vegetation index, normalized difference red-edge index, etc.) in both latitudes. Nevertheless, the simple ratio index (SR) showed less variability across SZA in both latitude zones amongst these indices. In addition, it was interesting to note that the correlation between vegetation indices and ground-reference data remained stable across SZA in both latitude zones. In summary, SR was found to have a minimum anisotropic reflectance effect in both zones, and the other vegetation indices can be utilized to evaluate relative differences in crop performances, although the absolute data would be affected by SZA.

The Interplay between Canopy Structure and Topography and Its Impacts on Seasonal Variations in Surface Reflectance Patterns in the Boreal Region of Alaska—Implications for Surface Radiation Budget

Forests play an essential role in maintaining the Earth’s overall energy balance. The variability in forest canopy structure, topography, and underneath vegetation background conditions create uncertainty in modeling solar radiation at the Earth’s surface, particularly for boreal regions in high latitude. The purpose of this study is to analyze seasonal variation in visible, near-infrared, and shortwave infrared reflectance with respect to land cover classes, canopy structures, and topography in a boreal region of Alaska. We accomplished this investigation by fusing Landsat 8 images and LiDAR-derived canopy structural data and multivariate statistical analysis. Our study shows that canopy structure and topography interplay and influence reflectance spectra in a complex way, particularly during the snow season. We observed that deciduous trees, also tall with greater rugosity, are more dominant on the southern slope than on the northern slope. Taller trees are typically seen in higher elevations regardless of vegetation types. Surface reflectance in all studied wavelengths shows similar relationships with canopy cover, height, and rugosity, mainly due to close connections between these parameters. Visible and near-infrared reflectance decreases with canopy cover, tree height, and rugosity, especially for the evergreen forest. Deciduous forest shows more considerable variability of surface reflectance in all studied wavelengths, particularly in March, mainly due to the mixing effect of snow and vegetation. The multivariate statistical analysis demonstrates a significant tree shadow effect on surface reflectance for evergreen forests. However, the topographic shadow effect is prominent for deciduous forests during the winter season. These results provide great insight into understanding the role of vegetation structure and topography in surface radiation budget in the boreal region.

USRT: A Solar Radiative Transfer Model Dedicated to Estimating Urban 3D Surface Reflectance

Urban 3D surface reflectance is a critical parameter for the modeling of surface biophysical processes. It is of great significance to enhance the accuracy of reflectance in urban areas. Based on the urban solar radiative transfer (USRT) model, this study presents a methodology for estimating urban reflectance using the sky view factor (SVF) derived from airborne LiDAR data. Then, the USRT model was used to retrieve urban 3D surface reflectance from Landsat 8 data over the typical area of Beijing. The reflectance from USRT model was compared with the estimated value obtained from the model without considering the impact of morphological characteristics of the urban underlying surface (flat model). The results showed that the urban sample reflectance estimated by the USRT model was close to the sample reflectance of the suburban underlying surface which was less affected by morphological characteristics. The research summaries are as follows: (1) The definite physical meaning is presented in the USRT model, and can be applied to estimate the physical parameters of the urban underlying surface. (2) The reflectance from the USRT model is slightly larger than the reflectance derived from the flat model, which indicates that the accuracy of urban 3D surface reflectance is improved by the USRT model. (3) The effects of the SVF and building reflectance are different. The SVF presents a strong sensitivity to the estimation of the urban 3D surface reflectance, and the variations of building reflectance setting have little impact on urban reflectance, which is characterized by low sensitivity. Generally, the methodology of estimating urban reflectance proposed in this study can better clarify the impact mechanism of urban geometry on the radiative transfer processes and further promote the application and development of urban quantitative remote sensing.

Evaluation and Normalization of Topographic Effects on Vegetation Indices

The normalization of topographic effects on vegetation indices (VIs) is a prerequisite for their proper use in mountainous areas. We assessed the topographic effects on the normalized difference vegetation index (NDVI), the enhanced vegetation index (EVI), the soil adjusted vegetation index (SAVI), and the near-infrared reflectance of terrestrial vegetation (NIRv) calculated from Sentinel-2. The evaluation was based on two criteria: the correlation with local illumination condition and the dependence on aspect. Results show that topographic effects can be neglected for the NDVI, while they heavily influence the SAVI, EVI, and NIRv: the local illumination condition explains 19.85%, 25.37%, and 26.69% of the variation of the SAVI, EVI, and NIRv, respectively, and the coefficients of variation across different aspects are, respectively, 8.13%, 10.46%, and 14.07%. We demonstrated the applicability of existing correction methods, including statistical-empirical (SE), sun-canopy-sensor with C-correction (SCS + C), and path length correction (PLC), dedicatedly designed for reflectance, to normalize topographic effects on VIs. Our study will benefit vegetation monitoring with VIs over mountainous areas.

Spatiotemporal Variability of Land Surface Albedo over the Tibet Plateau from 2001 to 2019

As an essential climate variable (ECV), land surface albedo plays an important role in the Earth surface radiation budget and regional or global climate change. The Tibetan Plateau (TP) is a sensitive environment to climate change, and understanding its albedo seasonal and inter-annual variations is thus important to help capture the climate change rules. In this paper, we analyzed the large-scale spatial patterns, temporal trends, and seasonal variability of land surface albedo overall the TP, based on the moderate resolution imaging spectroradiometer (MODIS) MCD43 albedo products from 2001 to 2019. Specifically, we assessed the correlations between the albedo anomaly and the anomalies of normalized difference vegetation index (NDVI), the fraction of snow cover (snow cover), and land surface temperature (LST). The results show that there are larger albedo variations distributed in the mountainous terrain of the TP. Approximately 10.06% of the land surface is identified to have been influenced by the significant albedo variation from the year 2001 to 2019. The yearly averaged albedo was decreased significantly at a rate of 0.0007 (Sen’s slope) over the TP. Additionally, the yearly average snow cover was decreased at a rate of 0.0756. However, the yearly average NDVI and LST were increased with slopes of 0.0004 and 0.0253 over the TP, respectively. The relative radiative forcing (RRF) caused by the land cover change (LCC) is larger than that caused by gradual albedo variation in steady land cover types. Overall, the RRF due to gradual albedo variation varied from 0.0005 to 0.0170 W/m2, and the RRF due to LCC variation varied from 0.0037 to 0.0243 W/m2 during the years 2001 to 2019. The positive RRF caused by gradual albedo variation or the LCC can strengthen the warming effects in the TP. The impact of the gradual albedo variations occurring in the steady land cover types was very low between 2001 and 2019 because the time series was short, and it therefore cannot be neglected when examining radiative forcing for a long time series regarding climate change.

A Kernel-Driven BRDF Approach to Correct Airborne Hyperspectral Imagery over Forested Areas with Rugged Topography

Airborne hyper-spectral imaging has been proven to be an efficient means to provide new insights for the retrieval of biophysical variables. However, quantitative estimates of unbiased information derived from airborne hyperspectral measurements primarily require a correction of the anisotropic scattering properties of the land surface depicted by the bidirectional reflectance distribution function (BRDF). Hitherto, angular BRDF correction methods rarely combined viewing-illumination geometry and topographic information to achieve a comprehensive understanding and quantification of the BRDF effects. This is in particular the case for forested areas, frequently underlaid by rugged topography. This paper describes a method to correct the BRDF effects of airborne hyperspectral imagery over forested areas overlying rugged topography, referred in the reminder of the paper as rugged topography-BRDF (RT-BRDF) correction. The local viewing and illumination geometry are calculated for each pixel based on the characteristics of the airborne scanner and the local topography, and these two variables are used to adapt the Ross-Thick-Maignan and Li-Transit-Reciprocal kernels in the case of rugged topography. The new BRDF model is fitted to the anisotropy of multi-line airborne hyperspectral data. The number of pixels is set at 35,000 in this study, based on a stratified random sampling method to ensure a comprehensive coverage of the viewing and illumination angles and to minimize the fitting error of the BRDF model for all bands. Based on multi-line airborne hyperspectral data acquired with the Chinese Academy of Forestry’s LiDAR, CCD, and Hyperspectral system (CAF-LiCHy) in the Pu’er region (China), the results applying the RT-BRDF correction are compared with results from current empirical (C, and sun-canopy-sensor (SCS) adds C (SCS+C)) and semi-physical (SCS) topographic correction methods. Both quantitative assessment and visual inspection indicate that RT-BRDF, C, and SCS+C correction methods all reduce the topographic effects. However, the RT-BRDF method appears more efficient in reducing the variability in reflectance of overlapping areas in multiple flight-lines, with the advantage of reducing the BRDF effects caused by the combination of wide field of view (FOV) airborne scanner, rugged topography, and varying solar illumination angle over long flight time. Specifically, the average decrease in coefficient of variation (CV) is 3% and 3.5% for coniferous forest and broadleaved forest, respectively. This improvement is particularly marked in the near infrared (NIR) region (i.e., >750 nm). This finding opens new possible applications of airborne hyperspectral surveys over large areas.

Assessment of five satellite-derived LAI datasets for GPP estimations through ecosystem models

Ecosystem models have been widely used for obtaining gross primary productivity (GPP) estimations at multiple scales. Leaf area index (LAI) is a critical variable in these models for describing the vegetation canopy structure and predicting vegetation-atmosphere interactions. However, the uncertainties in LAI datasets and the effects of their representation on simulated GPP remain unclear, especially over complex terrain. Here, five most popular datasets, namely the Long-term Global Mapping (GLOBMAP) LAI, Global LAnd Surface Satellite (GLASS) LAI, Geoland2 version 1 (GEOV1) LAI, Global Inventory Monitoring and Modeling System (GIMMS) LAI, and Moderate Resolution Imaging Spectroradiometer (MODIS) LAI, were selected to examine the influences of LAI representation on GPP estimations at 95 eddy covariance (EC) sites. The GPP estimations from the Boreal Ecosystem Productivity Simulator (BEPS) model and the Eddy Covariance Light Use Efficiency (EC-LUE) model were evaluated against EC GPP to assess the performances of LAI datasets. Results showed that MODIS LAI had stronger linear correlations with GLASS and GEOV1 than GIMMS and GLOMAP at the study sites. The GPP estimations from GLASS LAI had a better agreement with EC GPP than those from other four LAI datasets at forest sites, while the GPP estimations from GEOVI LAI matched best with EC GPP at grass sites. Additionally, the GPP estimations from GLASS and GEOVI LAI presented better performances than the other three LAI datasets at crop sites. Besides, the results also showed that complex terrain had larger discrepancies of LAI and GPP estimations, and flat terrain presented better performances of LAI datasets in GPP estimations. Moreover, the simulated GPP from BEPS was more sensitive to LAI than those from EC - LUE, suggesting that LAI datasets can also lead to different uncertainties in GPP estimations from different model structures. Our study highlights that the satellite-derived LAI datasets can cause uncertainties in GPP estimations through ecosystem models.

Evaluating the conservation state of the páramo ecosystem: An object-based image analysis and CART algorithm approach for central Ecuador

Ecuadorian páramo ecosystems (EPEs) function as water sources, contain large soil carbon stores and high levels of biodiversity, and support human populations. The EPEs are mainly herbaceous páramo (HP). To inform policy and management and help drive ecological science toward a better understanding of the HP ecosystem, and the relationships among its multiple ecosystem services, we asked: (1) What is the state of the HP regarding its land use/land cover (LULC)?; and (2) Is the HP being pushed away from its natural state or it is regenerating? To answer these questions, we assessed the LULC in central EPEs using Landsat 8 imagery, Object-Based Image Analysis (OBIA) and a Classification and Regression Trees (CART) algorithm. Results show that two-fifths of the paramo ecosystem remain as native HP (NHP) and two-fifths as anthropogenic HP (AHP). Although the anthropic alteration of the pedogenesis of young paramo soil leads to the establishment of AHP, we found evidence of regeneration and resilience of the NHP. The results of this study will be useful to scientists and decision-makers with interest in páramo ecosystems in central Ecuador. The proposed methodology is simple, fast, and could be implemented in other landscapes to establish comprehensive monitoring systems useful in landscape assessment and planning.

High Spatio-Temporal Resolution CYGNSS Soil Moisture Estimates Using Artificial Neural Networks

This paper presents a learning-based, physics-aware soil moisture (SM) retrieval algorithm for NASA’s Cyclone Global Navigation Satellite System (CYGNSS) mission. The goal of the proposed novel method is to advance CYGNSS-based SM estimations, exploiting the spatio-temporal resolution of the GNSS reflectometry (GNSS-R) signals to its highest potential within a machine learning framework. The methodology employs a fully connected Artificial Neural Network (ANN) regression model to perform SM predictions through learning the nonlinear relations of SM and other land geophysical parameters to the CYGNSS observables. In situ SM measurements from several International SM Network (ISMN) sites are used as reference labels; CYGNSS incidence angles, derived reflectivity and trailing edge slope (TES) values, as well as ancillary data, are exploited as input features for training and validation of the ANN model. In particular, the utilized ancillary data consist of normalized difference vegetation index (NDVI), vegetation water content (VWC), terrain elevation, terrain slope, and h-parameter (surface roughness). Land cover classification and inland water body masks are also used for the intermediate derivations and quality control purposes. The proposed algorithm assumes uniform SM over a 0.0833 ∘ × 0.0833 ∘ (approximately 9 km × 9 km around the equator) lat/lon grid for any CYGNSS observation that falls within this window. The proposed technique is capable of generating sub-daily and high-resolution SM predictions as it does not rely on time-series or spatial averaging of the CYGNSS observations. Once trained on the data from ISMN sites, the model is independent from other SM sources for retrieval. The estimation results obtained over unseen test data are promising: SM predictions with an unbiased root mean squared error of 0.0544 cm 3 /cm 3 and Pearson correlation coefficient of 0.9009 are reported for 2017 and 2018.

Estimating Forest Canopy Height Using MODIS BRDF Data Emphasizing Typical-Angle Reflectances

Forest-canopy height is an important parameter for the estimation of forest biomass and terrestrial carbon flux and climate-change research at regional and global scales. Currently, various methods combining Light Detection and Ranging (LiDAR) data with various auxiliary data, particularly satellite remotely sensed reflectances, have been widely used to produce spatially continuous canopy-height products. However, current methods in use for remote sensing reflectances mainly focus on the nadir view direction, while anisotropic reflectances, which are theoretically more sensitive to the forest canopy height in the multiangle remote sensing field, have rarely been explored. Here, we attempted to examine the potential of using modeled multiangle reflectances at three typical viewing angles (i.e., from the hotspot, darkspot, and nadir directions) to estimate forest-canopy height as auxiliary data sources. First, the sensitivities of the typical angular reflectances as a function of forest canopy height were fully examined using the Extended Fourier Amplitude Sensitivity Test (EFAST) method based on the 4-scale Bidirectional Reflectance Distribution Function (BRDF) model simulations. This indicated that reflectances in the off-nadir viewing directions are generally sensitive to canopy-height variations. Then, the canopy heights were extracted from airborne Laser Vegetation Imaging Sensor (LVIS) data, which were further divided into training and validation data. Moderate Resolution Imaging Spectroradiometer (MODIS) multiangle reflectances at typical viewing angles were calculated from the MODIS BRDF parameter product (MCD43A1, version 6) as partial training-input data, based on a hotspot-adjusted, kernel-driven linear BRDF model. Subsequently, the Random Forest (RF) machine learning model was trained to acquire the relationship between the extracted canopy heights and the corresponding MODIS typical viewing reflectances. The trained model was further applied to estimate the canopy height metrics in the study areas of Howland Forest, Harvard Forest, and Bartlett Forest. Finally, the estimated canopy heights were independently validated by canopy heights extracted from the LVIS data. The results indicate that the canopy heights modeled through this method exhibit generally high accordance with the LVIS-derived canopy heights (R = 0.65−0.67; RMSE = 3.63−5.78). The results suggest that the MODIS multiangle reflectance data at typical observation angles contain important information regarding forest canopy height and can, therefore, be used to estimate forest canopy height for various ecological applications.

A Flight Direction Design Method for Airborne Spectral Imaging Considering the Anisotropy Reflectance of the Target in Rugged Terrain

An excellent mission plan is the prerequisite for the acquisition of high quality airborne hyperspectral images which are useful for environmental research, mining etc. In order to minimize the radiance non-uniformity caused by the anisotropic reflectance of targets, the flight direction is mostly designed on the solar azimuth or 180° from it for whiskbroom and pushbroom imagers. However, the radiance to the observer is determined not only by the reflectance of the target, but also by the terrain, the illumination direction and the observation direction. So, the flight direction which is defined to minimize radiance non-uniformity might change with the terrain. In order to find the best flight direction for rugged terrain, we firstly analyze the causes of the effect of terrain on radiation non-uniformity based on the radiative transfer process. Then, the flight direction design method is proposed for composite sloping terrain. Tested by digital and physical simulation experiments, the radiance non-uniformity is minimized when the aircraft flies in the designated direction. Finally, a workflow for flight direction planning and optimizing is summarized, considering the flight mission planning techniques and the workflow of remote sensing missions.

Recent Progress in Quantitative Land Remote Sensing in China

During the past forty years, since the first book with a title mentioning quantitative and remote sensing was published [1], quantitative land remote sensing has advanced dramatically, and numerous books have been published since then [2–6] although some of them did not use quantitative land remote sensing in their titles. [...]

Diffuse Skylight as a Surrogate for Shadow Detection in High-Resolution Imagery Acquired Under Clear Sky Conditions

An alternative technique for shadow detection and abundance is presented for high spatial resolution imagery acquired under clear sky conditions from airborne/spaceborne sensors. The method, termed Scattering Index (SI), uses Rayleigh scattering principles to create a diffuse skylight vector as a shadow reference. From linear algebra, the proportion of diffuse skylight in each image pixel provides a per pixel measure of shadow extent and abundance. We performed a comparative evaluation against two other methods, first valley detection thresholding (extent) and physics-based unmixing (extent and abundance). Overall accuracy and F-score measures are used to evaluate shadow extent on both Worldview-3 and ADS40 images captured over a common scene. Image subsets are selected to capture objects well documented as shadow detection anomalies, e.g., dark water bodies. Results showed improved accuracies and F-scores for shadow extent and qualitative evaluation of abundance show the method is invariant to scene and sensor characteristics. SI avoids shadow misclassifications by avoiding the use of pixel intensity and the associated limitations of binary thresholding. The method negates the need for complex sun-object-sensor corrections, it is simple to apply, and it is invariant to the exponential increase in scene complexity associated with higher-resolution imagery.