Vegetation Optical Depth Data Research Articles

High-precision sugarcane yield prediction by integrating 10-m Sentinel-1 VOD and Sentinel-2 GRVI indexes

Developing accurate models of sugarcane yield mapping at fine-scale is of paramount importance and will benefit many aspects of managing growth and harvest of sugarcane crops. Here, we combined high-spatial-resolution multi-sensor (optical and microwave) remote sensing data to estimate sugarcane yield using two approaches. First, we retrieved the 10 m resolution vegetation optical depth (VOD) using C-band Sentinel-1 synthetic aperture radar data via the water cloud model, then derived a combined VOD and Sentinel-2 derived green–red vegetation index (GRVI, a traditional vegetation index) time series. Second, we adopted the eXtreme Gradient Boost (XGBoost) machine learning algorithm to estimate total sugarcane production based on the aggregated mean monthly VOD and GRVI time series at the county scale. County-scale sugarcane yield data from 81 counties across Guangxi, China (2018–2020) were used for model training and testing. We then kept the trained best-performance model unchanged to re-predict pixel-scale (10-m resolution) yield using VOD and GRVI images. A SHapley Additive exPlanations (SHAP) approach was employed to test XGBoost knowledge regarding the mechanisms affecting yield. The main outcomes from this study were as follows. (1) The XGBoost model and VOD and GRVI time series unequivocally provided reliable and precise estimates of sugarcane yield, with a low relative error both at the county scale and at the pixel scale. (2) Predictions of sugarcane yield in the growing season using combined VOD and GRVI data yields high accuracy (R2 = 0.815) with specific spatial patterns of sugarcane yields for Guangxi province, and were feasible up to two months prior to harvest. (3) The SHAP model revealed that Sentinel-1 derived VOD data are important for predicting yield. This research portrays the benefits of combining VOD and GRVI data in automatic machine learning models for estimating sugarcane yield at high spatial and temporal resolutions.

Quantifying Uncertainties in Passive Microwave Remote Sensing of Soil Moisture via a Bayesian Probabilistic Inversion Method

Improving the accuracy of remotely sensed soil moisture (SM) is a challenging and popular topic. Quantifying the uncertainty in the SM inversion process and enhancing the confidence of SM retrieval are promising ways to address this challenge but have received little attention. We present a Bayesian probabilistic inversion algorithm that can simultaneously retrieve SM, surface roughness, and vegetation optical depth data and quantify the uncertainty in the inversion. The proposed algorithm is evaluated using airborne polarimetric L-band multibeam radiometer (PLMR) observations. We use three combinations, 3-angular observations at V-polarization (3CV), 3-angular observations at H-polarization (3CH) and 6-channel observations (6CA), to identify the optimal configuration for SM retrieval by taking advantage of the PLMR’s dual-polarization and multiple angles. Uncertainties are quantified by introducing multiple uncertainty quantification metrics into Bayesian posterior distributions of SM retrievals. The estimates are validated against multiscale ground-based measurements, including manual measurements and wireless sensor network (WSN) measurements, and the spatial representativeness of the ground-based reference regarding the validation of pixel-scale SM retrievals is discussed. The 6CA attempt yields the best SM estimates (correlation coefficient (R) ≥ 0.864, root mean square error (RMSE) ≤0.04 m3/m3, and unbiased RMSE (unRMSE) ≤ 0.035 m3/m3), while the 3CH attempt yields the lowest uncertainty. Additionally, dense manual measurements are more representative than sparsely distributed WSN measurements. Overall, combining dual-polarized observations yields the best SM estimates but introduces additional uncertainty. This study highlights uncertainties quantification in SM inversion and thus provides confidence in SM inversion, facilitating improved SM retrieval algorithms.

Leaf surface water, not plant water stress, drives diurnal variation in tropical forest canopy water content.

Variation in canopy water content (CWC) that can be detected from microwave remote sensing of vegetation optical depth (VOD) has been proposed as an important measure of vegetation water stress. However, the contribution of leaf surface water (LWs ), arising from dew formation and rainfall interception, to CWC is largely unknown, particularly in tropical forests and other high-humidity ecosystems. We compared VOD data from the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E) and CWC predicted by a plant hydrodynamics model at four tropical sites in Brazil spanning a rainfall gradient. We assessed how LWs influenced the relationship between VOD and CWC. The analysis indicates that while CWC is strongly correlated with VOD (R2 =0.62 across all sites), LWs accounts for 61-76% of the diurnal variation in CWC despite being < 10% of CWC. Ignoring LWs weakens the near-linear relationship between CWC and VOD and reduces the consistency in diurnal variation. The contribution of LWs to CWC variation, however, decreases at longer, seasonal to inter-annual, time scales. Our results demonstrate that diurnal patterns of dew formation and rainfall interception can be an important driver of diurnal variation in CWC and VOD over tropical ecosystems and therefore should be accounted for when inferring plant diurnal water stress from VOD measurements.

Simultaneous retrieval of global scale Vegetation Optical Depth, surface roughness, and soil moisture using X-band AMSR-E observations

The radiative transfer scheme implemented for the retrieval of soil moisture from passive microwaves is a function of scattering, polarization mixing and attenuation effects of soil and vegetation. Theses factors are usually represented by Vegetation Optical Depth (VOD), vegetation scattering albedo, and surface roughness parameter, along with soil moisture. The VOD is the degree to which vegetation attenuates the microwave radiation. It has generally the dominant effect from vegetation, whereas scattering is negligible and close to zero. The surface roughness (which varies in space but not much in time) is until recently, often assumed to be a global constant. In this work, we attempted to simultaneously retrieve the VOD, the surface roughness parameter, and the soil moisture at the global scale using the Level 3 daily 0.25° X-band brightness temperatures of the Advanced Microwave Scanning Radiometer – Earth Observing System (AMSR-E) sensor. The methodology, coined as the Simultaneous Parameter Retrieval Algorithm (SPRA), is based on the premise that the vegetation dynamics undergo slower temporal changes than the soil moisture - an assumption, which is successfully used in the past for microwave radiometric retrievals at lower frequencies. Results indicate that the SPRA produces the VOD retrievals with reduced high-frequency noise when compared to the baseline Land Parameter Retrieval Algorithm (LPRM) retrievals. This effect assisted in identifying the influence of precipitation and cropping patterns on the temporal dynamics of the VOD. Good agreement is observed between the mean SPRA VOD and canopy height data (global correlation = 0.75). The spatial patterns of surface roughness parameter agree well with the roughness product (HR map) developed using Soil Moisture Ocean Salinity (SMOS) sensor based data (global correlation = 0.57). Validation of SPRA and LPRM soil moisture products with in-situ observations over the Contiguous United States (CONUS) indicated an improvement in mean ubRMSE with SPRA product (SPRA-0.11 m3/m3 and LPRM-0.18 m3/m3) and comparable mean Pearson correlations between the two products (SPRA-0.36 and LPRM-0.38). Further, a precipitation based consistency evaluation of SPRA and LPRM soil moisture retrievals indicated better skill of the SPRA product over India.

A carbon sink-driven approach to estimate gross primary production from microwave satellite observations

Global estimation of Gross Primary Production (GPP) - the uptake of atmospheric carbon dioxide by plants through photosynthesis - is commonly based on optical satellite remote sensing data. This presents a source-driven approach since it uses the amount of absorbed light, the main driver of photosynthesis, as a proxy for GPP. Vegetation Optical Depth (VOD) estimates obtained from microwave sensors provide an alternative and independent data source to estimate GPP on a global scale, which may complement existing GPP products. Recent studies have shown that VOD is related to aboveground biomass, and that both VOD and temporal changes in VOD relate to GPP. In this study, we build upon this concept and propose a model for estimating GPP from VOD. Since the model is driven by vegetation biomass, as observed through VOD, it presents a carbon sink-driven approach to quantify GPP and, therefore, is conceptually different from common source-driven approaches. The model developed in this study uses single frequencies from active or passive microwave VOD retrievals from C-, X- and Ku-band (Advanced Scatterometer (ASCAT) and Advanced Microwave Scanning Radiometer for Earth Observation (AMSR-E)) to estimate GPP at the global scale. We assessed the ability for temporal and spatial extrapolation of the model using global GPP from FLUXCOM and in situ GPP from FLUXNET. We further performed upscaling of in situ GPP based on different VOD data sets and compared these estimates with the FLUXCOM and MODerate-resolution Imaging Spectroradiometer (MODIS) GPP products. Our results show that the model developed for individual grid cells using VOD and change in VOD as input performs well in predicting temporal patterns in GPP for all VOD data sets. For spatial extrapolation of the model, however, additional input variables are needed to represent the spatial variability of the VOD-GPP relationship due to differences in vegetation type. As additional input variable, we included the grid cell median VOD (as a proxy for vegetation cover), which increased the model performance during cross validation. Mean annual GPP obtained for AMSR-E X-band data tends to overestimate mean annual GPP for FLUXCOM and MODIS but shows comparable latitudinal patterns. Overall, our findings demonstrate the potential of VOD for estimating GPP. The sink-driven approach provides additional information about GPP independent of optical data, which may contribute to our knowledge about the carbon source-sink balance in different ecosystems.

Vegetation optical depth at L-band and above ground biomass in the tropical range: Evaluating their relationships at continental and regional scales

The relationship between vegetation optical depth (VOD) retrieved by L-band SMOS radiometer and forest above ground biomass (AGB) was investigated in tropical areas of Africa and South America. VOD was retrieved from the latest version of level 2 SMOS algorithm, while reference AGB was obtained from a pantropical database, encompassing a large number of ground plot data derived from field surveys conducted on both continents. In Africa and South-America, VOD increased with AGB, reaching saturation at about 350 Mg ha−1. The strength of the relation was improved selecting VOD data in appropriate seasons, characterized by a higher dynamic range of values. The capability of VOD data to estimate AGB was further investigated using Random Forest decision trees, adding to VOD selected climate variables from the Climatic Research Unit (temperature, potential evapotranspiration, and precipitation) and water deficit data, and validating regression tests with ground data from the reference AGB database. The results for the five analyzed years indicate that the best estimates of AGB are obtained by the joined use of VOD and potential evapotranspiration input data, but all climate variables brought an improvement in AGB estimates. AGB estimates were relatively stable for the considered period, with limited variations possibly due to changes in biomass and to data quality of VOD and of climate variables. The VOD signal and estimated AGB were also analyzed according to ecological homogeneous units (ecoregions), evidencing data clusters, partially overlapped to each other, in the VOD - AGB plane.

An evaluation of SMOS L-band vegetation optical depth (L-VOD) data sets: high sensitivity of L-VOD to above-ground biomass in Africa

Abstract. The vegetation optical depth (VOD) measured at microwave frequencies is related to the vegetation water content and provides information complementary to visible/infrared vegetation indices. This study is devoted to the characterization of a new VOD data set obtained from SMOS (Soil Moisture and Ocean Salinity) satellite observations at L-band (1.4 GHz). Three different SMOS L-band VOD (L-VOD) data sets (SMOS level 2, level 3 and SMOS-IC) were compared with data sets on tree height, visible/infrared indexes (NDVI, EVI), mean annual precipitation and above-ground biomass (AGB) for the African continent. For all relationships, SMOS-IC showed the lowest dispersion and highest correlation. Overall, we found a strong (R &gt; 0.85) correlation with no clear sign of saturation between L-VOD and four AGB data sets. The relationships between L-VOD and the AGB data sets were linear per land cover class but with a changing slope depending on the class type, which makes it a global non-linear relationship. In contrast, the relationship linking L-VOD to tree height (R = 0.87) was close to linear. For vegetation classes other than evergreen broadleaf forest, the annual mean of L-VOD spans a range from 0 to 0.7 and it is linearly correlated with the average annual precipitation. SMOS L-VOD showed higher sensitivity to AGB compared to NDVI and K/X/C-VOD (VOD measured at 19, 10.7 and 6.9 GHz). The results showed that, although the spatial resolution of L-VOD is coarse ( ∼ 40 km), the high temporal frequency and sensitivity to AGB makes SMOS L-VOD a very promising indicator for large-scale monitoring of the vegetation status, in particular biomass.

Assessing the relationship between microwave vegetation optical depth and gross primary production

At the global scale, the uptake of atmospheric carbon dioxide by terrestrial ecosystems through photosynthesis is commonly estimated through vegetation indices or biophysical properties derived from optical remote sensing data. Microwave observations of vegetated areas are sensitive to different components of the vegetation layer than observations in the optical domain and may therefore provide complementary information on the vegetation state, which may be used in the estimation of Gross Primary Production (GPP). However, the relation between GPP and Vegetation Optical Depth (VOD), a biophysical quantity derived from microwave observations, is not yet known. This study aims to explore the relationship between VOD and GPP. VOD data were taken from different frequencies (L-, C-, and X-band) and from both active and passive microwave sensors, including the Advanced Scatterometer (ASCAT), the Soil Moisture Ocean Salinity (SMOS) mission, the Advanced Microwave Scanning Radiometer for Earth Observation System (AMSR-E) and a merged VOD data set from various passive microwave sensors. VOD data were compared against FLUXCOM GPP and Solar-Induced chlorophyll Fluorescence (SIF) from the Global Ozone Monitoring Experiment-2 (GOME-2). FLUXCOM GPP estimates are based on the upscaling of flux tower GPP observations using optical satellite data, while SIF observations present a measure of photosynthetic activity and are often used as a proxy for GPP. For relating VOD to GPP, three variables were analyzed: original VOD time series, temporal changes in VOD (ΔVOD), and positive changes in VOD (ΔVOD≥0). Results show widespread positive correlations between VOD and GPP with some negative correlations mainly occurring in dry and wet regions for active and passive VOD, respectively. Correlations between VOD and GPP were similar or higher than between VOD and SIF. When comparing the three variables for relating VOD to GPP, correlations with GPP were higher for the original VOD time series than for ΔVOD or ΔVOD≥0 in case of sparsely to moderately vegetated areas and evergreen forests, while the opposite was true for deciduous forests. Results suggest that original VOD time series should be used jointly with changes in VOD for the estimation of GPP across biomes, which may further benefit from combining active and passive VOD data.

Comparing Passive Microwave with Visible-To-Near-Infrared Phenometrics in Croplands of Northern Eurasia

Planting and harvesting times drive cropland phenology. There are few datasets that derive explicit phenological metrics, and these datasets use the visible to near infrared (VNIR) spectrum. Many different methods have been used to derive phenometrics such as Start of Season (SOS) and End of Season (EOS), leading to differing results. This discrepancy is partly due to spatial and temporal compositing of the VNIR satellite data to minimize data gaps resulting from cloud cover, atmospheric aerosols, and solar illumination constraints. Phenometrics derived from the downward Convex Quadratic model (CxQ) include Peak Height (PH) and Thermal Time to Peak (TTP), which are more consistent than SOS and EOS because they are minimally affected by snow and frost and other non-vegetation related issues. Here, we have determined PH using the vegetation optical depth (VOD) in three microwave frequencies (6.925, 10.65 and 18.7 GHz) and accumulated growing degree-days derived from AMSR-E (Advanced Microwave Scanning Radiometer on EOS) data at a spatial resolution of 25 km. We focus on 50 AMSR-E cropland pixels in the major grain production areas of Northern Eurasia (Ukraine, southwestern Russia, and northern Kazakhstan) for 2003–2010. We compared the land surface phenologies of AMSR-E VOD and MODIS NDVI data. VOD time series tracked cropland seasonal dynamics similar to that recorded by the NDVI. The coefficients of determination for the CxQ model fit of the NDVI data were high for all sites (0.78 &lt; R2 &lt; 0.99). The 10.65 GHz VOD (VOD1065GHz) achieved the best linear regression fit (R2 = 0.84) with lowest standard error (SEE = 0.128); it is therefore recommended for microwave VOD studies of cropland land surface phenology. Based on an Analysis of Covariance (ANCOVA) analysis, the slopes from the linear regression fit were not significantly different by microwave frequency, whereas the intercepts were significantly different, given the different magnitudes of the VODs. PHs for NDVI and VOD were highly correlated. Despite their strong correspondence, there was generally a lag of AMSR-E PH VOD10.65GHz by about two weeks compared to MODIS peak greenness. To evaluate the utility of the PH determination based on maximum value, we correlated the CxQ derived and maximum value determined PHs of NDVI and found that they were highly correlated with R2 of 0.87, but with a one-week bias. Considering the one-week bias between the two methods, we find that PH of VOD10.65GHz lags PH of NDVI by three weeks. We conclude, therefore, that maximum-value based PH of VOD can be a complementary phenometric for the CxQ model derived PH NDVI, especially in cloud and aerosol obscured regions of the world.

Analysis of SMOS brightness temperature and vegetation optical depth data with coupled land surface and radiative transfer models in Southern Germany

Abstract. Soil Moisture and Ocean Salinity (SMOS) L1c brightness temperature and L2 optical depth data are analysed with a coupled land surface (PROMET) and radiative transfer model (L-MEB). The coupled models are validated with ground and airborne measurements under contrasting soil moisture, vegetation and land surface temperature conditions during the SMOS Validation Campaign in May and June 2010 in the SMOS test site Upper Danube Catchment in southern Germany. The brightness temperature root-mean-squared errors are between 6 K and 9 K. The L-MEB parameterisation is considered appropriate under local conditions even though it might possibly be further optimised. SMOS L1c brightness temperature data are processed and analysed in the Upper Danube Catchment using the coupled models in 2011 and during the SMOS Validation Campaign 2010 together with airborne L-band brightness temperature data. Only low to fair correlations are found for this comparison (R between 0.1–0.41). SMOS L1c brightness temperature data do not show the expected seasonal behaviour and are positively biased. It is concluded that RFI is responsible for a considerable part of the observed problems in the SMOS data products in the Upper Danube Catchment. This is consistent with the observed dry bias in the SMOS L2 soil moisture products which can also be related to RFI. It is confirmed that the brightness temperature data from the lower SMOS look angles and the horizontal polarisation are less reliable. This information could be used to improve the brightness temperature data filtering before the soil moisture retrieval. SMOS L2 optical depth values have been compared to modelled data and are not considered a reliable source of information about vegetation due to missing seasonal behaviour and a very high mean value. A fairly strong correlation between SMOS L2 soil moisture and optical depth was found (R = 0.65) even though the two variables are considered independent in the study area. The value of coupled models as a tool for the analysis of passive microwave remote-sensing data is demonstrated by extending this SMOS data analysis from a few days during a field campaign to a longer term comparison.