Sentinel-2数据的冬小麦地上干生物量估算及评价

作物生物量快速精确的监测对于农业资源的合理利用与农田的精准管理具有重要意义。近年来,遥感技术因其独特的优势已被广泛用于作物生物量的估算中。本文主要针对不同宽波段植被指数在冬小麦生物量(文中的生物量均是指地上干生物量)估算方面的表现进行探索。首先利用欧洲空间局最新的Sentinel-2A卫星数据提取出17种常见的植被指数,之后分别构建其与相应时期内采集的冬小麦地上生物量间的最优估算模型,通过分析两者间的相关性与敏感性,获取适宜进行生物量估算的指数。最后,绘制了研究区的生物量空间分布图。结果表明,所选的植被指数均与生物量显著相关。其中,红边叶绿素指数(CI_(re))与生物量的估算精度最高(决定性系数R~2为0.83;均方根误差RMSE为180.29 g·m~(–2))。虽然相关性较高,但部分指数,如归一化差值植被指数(NDVI)等在生物量较高时会出现饱和现象,从而导致生物量的低估。而加入红边波段的指数不仅能够延缓指数的饱和趋势,而且能够提高反演精度。此外,通过敏感性分析发现,归一化差值指数和比值指数分别在作物生长的早期和中后期对生物量的变化保持较高的敏感性。由于红边比值指数(SR_(re))和MERIS叶绿素敏感指数(MTCI)在冬小麦全生长季内一直对生物量的变化保持高灵敏性,二者是生物量估算中最为稳定的指数。     

机载LiDAR数据估算样地和单木尺度森林地上生物量

利用机载激光雷达点云数据,结合大量实测单木结构信息,分别从样地和单木尺度估算了森林地上生物量AGB。首先,利用局部最大值单木提取算法提取了每个样地内的单木结构参数,并针对样地和单木尺度分别计算了一组激光雷达变量。然后,利用激光雷达变量和地上生物量及其两者的对数形式,从样地和单木尺度分别构建了估算模型。最后,针对两种尺度估算过程中存在的不确定性进行了详细讨论。结果表明:(1)样地和单木尺度模型估算的森林地上生物量与地面实测值都具有明显的相关性,且对数模型估算效果要优于非对数模型;(2)样地尺度模型估算效果(R2=0.84,rRMSE=0.23)明显优于单木尺度模型(R2=0.61,rRMSE=0.46);(3)按树木类型分别进行估算可以提高单木地上生物量的估算精度;(4)不论是样地还是单木尺度地上生物量估算都存在一定的不确定性,与样地尺度相比,单木尺度估算过程的不确定性更大,这种不确定性主要来自单木识别过程。     

WorldView-2纹理的森林地上生物量反演

使用高空间分辨率卫星WorldView-2的多光谱遥感影像,构建植被指数和纹理因子等遥感因子与森林地上生物量的关系方程,并计算模型估测精度和均方根误差,探索高分辨率数据的光谱与纹理信息在温带森林地上生物量估测应用中的潜力。以黑龙江省凉水自然保护区温带天然林及天然次生林为研究对象,通过灰度共生矩阵(GLCM)、灰度差分向量(GLDV)及和差直方图(SADH)对高分辨率遥感影像进行纹理信息提取,并利用外业调查的74个样地地上生物量与遥感因子建立参数估计模型。提取的遥感因子包括6种植被指数(比值植被指数RVI、差值植被指数DVI、规一化植被指数NDVI、增强植被指数EVI、土壤调节植被指数SAVI和修正的土壤调节植被指数MSAVI)以及3类纹理因子(GLCM、GLDV和SADH)。为避免特征变量个数较多对估测模型造成过拟合,利用随机森林算法对提取的遥感因子进行特征选择,将最优的特征变量输入模型参与建模估测。采用支持向量回归(SVR)进行生物量建模及验证,结果显示选入模型的和差直方图均值(sadh_mean)、灰度共生矩阵方差(glcm_var)和差值植被指数(DVI)等遥感因子对森林地上生物量有较好的解释效果;植被指数+纹理因子组合的模型获得较精确的AGB估算结果(R2=0.85,RMSE=42.30 t/ha),单独使用植被指数的模型精度则较低(R~2=0.69,RMSE=61.13 t/ha)。     

Estimating aboveground biomass of the mangrove forests on northeast Hainan Island in China using an upscaling method from field plots,UAV-LiDAR data and Sentinel-2 imagery

The mangrove forests of northeast Hainan Island are the most species diverse forests in China and consist of the Dongzhai National Nature Reserve and the Qinglan Provincial Nature Reserve. The former reserve is the first Chinese national nature reserve for mangroves and the latter has the most abundant mangrove species in China. However, to date the aboveground ground biomass (AGB) of this mangrove region has not been quantified due to the high species diversity and the difficulty of extensive field sampling in mangrove habitat. Although three-dimensional point clouds can capture the forest vertical structure, their application to large areas is hindered by the logistics, costs and data volumes involved. To fill the gap and address this issue, this study proposed a novel upscaling method for mangrove AGB estimation using field plots, UAV-LiDAR strip data and Sentinel-2 imagery (named G∼LiDAR∼S2 model) based on a point-line-polygon framework. In this model, the partial-coverage UAV-LiDAR data were used as a linear bridge to link ground measurements to the wall-to-wall coverage Sentinel-2 data. The results showed that northeast Hainan Island has a total mangrove AGB of 312,806.29 Mg with a mean AGB of 119.26 Mg ha⁻¹. The results also indicated that at the regional scale, the proposed UAV-LiDAR linear bridge method (i.e., G∼LiDAR∼S2 model) performed better than the traditional approach, which directly relates field plots to Sentinel-2 data (named the G∼S2 model) (R² = 0.62 > 0.52, RMSE = 50.36 Mg ha⁻¹<56.63 Mg ha⁻¹). Through a trend extrapolation method, this study inferred that the G∼LiDAR∼S2 model could decrease the number of field samples required by approximately 37% in comparison with those required by the G∼S2 model in the study area. Regarding the UAV-LiDAR sampling intensity, compared with the original number of LiDAR plots, 20% of original linear bridges could produce an acceptable accuracy (R² = 0.62, RMSE = 51.03 Mg ha⁻¹). Consequently, this study presents the first investigation of AGB for the mangrove forests on northeast Hainan Island in China and verifies the feasibility of using this mangrove AGB upscaling method for diverse mangrove forests.     

Predicting forest structural parameters using the image texture derived from WorldView-2 multispectral imagery in a dryland forest,Israel

Estimation of forest structural parameters by field-based data collection methods is both expensive and time consuming. Satellite remote sensing is a low-cost alternative in modeling and mapping structural parameters in large forest areas. The current study investigates the potential of using WordView-2 multispectral satellite imagery for predicting forest structural parameters in a dryland plantation forest in Israel. The relationships between image texture features and the several structural parameters such as Number of Trees (NT), Basal Area (BA), Stem Volume (SV), Clark-Evans Index (CEI), Diameter Differentiation Index (DDI), Contagion Index (CI), Gini Coefficient (GC), and Standard Deviation of Diameters at Breast Heights (SDDBH) were examined using correlation analyses. These variables were obtained from 30 m × 30 m square-shaped plots. The Standard Deviation of Gray Levels (SDGL) as a first order texture feature and the second order texture variables based on Gray Level Co-occurrence Matrix (GLCM) were calculated for the pixels that corresponds to field plots. The results of the correlation analysis indicate that the forest structural parameters are significantly correlated with the image texture features. The highest correlation coefficients were calculated for the relationships between the SDDBH and the contrast of red band (r = 0.75, p < 0.01), the BA and the entropy of blue band (r = 0.73, p < 0.01), and the GC and the contrast of blue band (r = 0.71, p < 0.01). Each forest structural parameter was modeled as a function of texture measures derived from the satellite image using stepwise multi linear regression analyses. The determination coefficient (R²) and root mean square error (RMSE) values of the best fitting models, respectively, are 0.38 and 109.56 ha⁻¹ for the NT; 0.54 and 1.79 m² ha⁻¹ for the BA; 0.42 and 27.18 m³ ha⁻¹ for the SV; 0.23 and 0.16 for the CEI; 0.32 and 0.05 for the DDI; 0.25 and 0.06 for the CI; 0.50 and 0.05 for the GC; and 0.67 and 0.70 for the SDDBH. The leave-one-out cross-validation technique was applied for validation of the best-fitted models (R² > 0.50). In conclusion, cross-validated statistics confirmed that the structural parameters including the BA, SDDBH, and GC can be predicted and mapped with a reasonable accuracy using the texture features extracted from the spectral bands of WorldView-2 image.     

Monitoring aboveground forest biomass dynamics over three decades using Landsat time-series and single-date inventory data

Understanding forest biomass dynamics is crucial for carbon and environmental monitoring, especially in the context of climate change. In this study, we propose a robust approach for monitoring aboveground forest biomass (AGB) dynamics by combining Landsat time-series with single-date inventory data. We developed a Random Forest (RF) based kNN model to produce annual maps of AGB from 1988 to 2017 over 7.2 million ha of forests in Victoria, Australia. The model was internally evaluated using a bootstrapping technique. Predictions of AGB and its change were then independently evaluated using multi-temporal Lidar data (2008 and 2016). To understand how natural and anthropogenic processes impact forest AGB, we analysed trends in relation to the history of disturbance and recovery. Specifically, change metrics (e.g., AGB loss and gain, Years to Recovery - Y2R) were calculated at the pixel level to characterise the patterns of AGB change resulting from forest dynamics. The imputation model achieved a RMSE value of 132.9 Mg ha⁻¹ (RMSE% = 46.3%) and R² value of 0.56. Independent assessments of prediction maps in 2008 and 2016 using Lidar-based AGB data achieved relatively high accuracies, with a RMSE of 108.6 Mg ha⁻¹ and 135.9 Mg ha⁻¹ for 2008 and 2016, respectively. Annual validations of AGB maps using un-changed, homogenous Lidar plots suggest that our model is transferable through time (RMSE ranging from 109.65 Mg ha⁻¹ to 112.27 Mg ha⁻¹ and RMSE% ranging from 25.38% to 25.99%). In addition, changes in AGB values associated with forest disturbance and recovery (decrease and increase, respectively) were captured by predicted maps. AGB change metrics indicate that AGB loss and Y2R varied across bioregions and were highly dependent on levels of disturbance severity (i.e., a greater loss and longer recovery time were associated with a higher severity disturbance). On average, high severity fire burnt from 200 Mg ha⁻¹ to 550 Mg ha⁻¹ of AGB and required up to 15 years to recover while clear-fell logging caused a reduction in 250 Mg ha⁻¹ to 600 Mg ha⁻¹ of AGB and required nearly 20 years to recover. In addition, AGB within un-disturbed forests showed statistically significant but monotonic trends, suggesting a mild gradual drop over time across most bioregions. Our methods are designed to support forest managers and researchers in developing forest monitoring systems, especially in developing regions, where only a single date forestry inventory exists.     

Inter-comparison of remote sensing platforms for height estimation of mango and avocado tree crowns

To support the adoption of precision agricultural practices in horticultural tree crops, prior research has investigated the relationship between crop vigour (height, canopy density, health) as measured by remote sensing technologies, to fruit quality, yield and pruning requirements. However, few studies have compared the accuracy of different remote sensing technologies for the estimation of tree height. In this study, we evaluated the accuracy, flexibility, aerial coverage and limitations of five techniques to measure the height of two types of horticultural tree crops, mango and avocado trees. Canopy height estimates from Terrestrial Laser Scanning (TLS) were used as a reference dataset against height estimates from Airborne Laser Scanning (ALS) data, WorldView-3 (WV-3) stereo imagery, Unmanned Aerial Vehicle (UAV) based RGB and multi-spectral imagery, and field measurements. Overall, imagery obtained from the UAV platform were found to provide tree height measurement comparable to that from the TLS (R² = 0.89, RMSE = 0.19 m and rRMSE = 5.37 % for mango trees; R² = 0.81, RMSE = 0.42 m and rRMSE = 4.75 % for avocado trees), although coverage area is limited to 1–10 km² due to battery life and line-of-sight flight regulations. The ALS data also achieved reasonable accuracy for both mango and avocado trees (R² = 0.67, RMSE = 0.24 m and rRMSE = 7.39 % for mango trees; R² = 0.63, RMSE = 0.43 m and rRMSE = 5.04 % for avocado trees), providing both optimal point density and flight altitude, and therefore offers an effective platform for large areas (10 km²–100 km²). However, cost and availability of ALS data is a consideration. WV-3 stereo imagery produced the lowest accuracies for both tree crops (R² = 0.50, RMSE = 0.84 m and rRMSE = 32.64 % for mango trees; R² = 0.45, RMSE = 0.74 m and rRMSE = 8.51 % for avocado trees) when compared to other remote sensing platforms, but may still present a viable option due to cost and commercial availability when large area coverage is required. This research provides industries and growers with valuable information on how to select the most appropriate approach and the optimal parameters for each remote sensing platform to assess canopy height for mango and avocado trees.     

Modeling and mapping aboveground biomass of the restored mangroves using ALOS-2 PALSAR-2 in East Kalimantan,Indonesia

Accurate estimation of forest aboveground biomass (AGB) using remote sensing is a requisite for monitoring, reporting and verification (MRV) system of the United Nations Programme on Reducing Emissions from Deforestation and Forest Degradation. However, attaining high accuracy remains a great challenge in the diverse tropical forests. Among available technologies, l-band Synthetic Aperture Radar (SAR) estimates AGB with reasonably high accuracy in the terrestrial tropical forests. Nevertheless, the accuracy is relatively low in the mangrove forests. In this context, the study was carried out to model and map AGB using backscatter coefficients of Advanced Land Observing Satellite-2 (ALOS-2) Phased Array l-band SAR-2 (PALSAR-2) in part of the restored mangrove forest at Mahakam Delta, Indonesia. PALSAR-2 data was acquired with image scene observation during the peak low tide on 30 July 2018 from Japan Aerospace Exploration Agency. The forest parameters namely tree height and diameter at breast height were measured from 71 field plots in September-October 2018. The parameters were used in mangrove allometry to calculate the field AGB. Finally, HV polarized backscatter coefficients of PALSAR-2 were used to model AGB using linear regression. The model demonstrated a comparatively high performance using three distinct methods viz. independent validation (R² of 0.89 and RMSE of 23.16 tons ha⁻¹), random k-fold cross validation (R² of 0.89 and RMSE of 24.59 tons ha⁻¹) and leave location out cross validation (LLO CV) (R² of 0.88 and RMSE of 24.05 tons ha⁻¹). The high accuracy of the LLO CV indicates no spatial overfitting in the model. Thus, the model based on LLO CV was used to map AGB in the study area. This is the first study that successfully obtains high accuracy in modeling AGB in the mangrove forest. Therefore, it offers a significant contribution to the MRV mechanism for monitoring mangrove forests in the tropics and sub-tropics.     

Evaluation of machine learning models for rice dry biomass estimation and mapping using quad-source optical imagery

ABSTRACT

Several machine learning regression models have been advanced for the estimation of crop biophysical parameters with optical satellite imagery. However, literature on the comparative performances of such models is still limited in range and scope, especially under multiple data sources, despite the potential of multi-source imagery to improving crop monitoring in cloudy areas. To fill in this knowledge gap, this study explored the synergistic use of Landsat-8, Sentinel-2A, China’s environment and disaster monitoring and forecasting satellites (HJ-1 A and B) and Gaofen-1 (GF-1) data to evaluate four machine learning regression models that include Random Forest (RF), Support Vector Machine (SVM), k-Nearest Neighbor (k-NN), and Gradient Boosting Decision Tree (GBDT), for rice dry biomass estimation and mapping. Taking a major rice cultivation area in southeast China as case study during the 2016 and 2017 growing seasons, a cross-calibrated time series of the Enhanced Vegetation Index (EVI) was obtained from the quad-source optical imagery and on which the aforementioned models were applied, respectively. Results indicate that in the before rice heading scenario, the most accurate dry biomass estimates were obtained by the GBDT model (R² of 0.82 and RMSE of 191.8 g/m²) followed by the RF model (R² of 0.79 and RMSE of 197.8 g/m²). After heading, the k-NN model performed best (R² of 0.43 and RMSE of 452.1 g/m²) followed by the RF model (R² of 0.42 and RMSE of 464.7 g/m²). Whist the k-NN model performed least in the before heading scenario, SVM performed least in the after heading scenario. These findings may suggest that machine learning regression models based on an ensemble of decision trees (RF and GBDT) are more suitable for the estimation of rice dry biomass, at least with optical satellite imagery. Studies that would extend the evaluation of these machine learning models, to other parameters like leaf area index, and to microwave imagery, are hereby recommended.     

Estimation of leaf area index using PROSAIL based LUT inversion,MLRA-GPR and empirical models: Case study of tropical deciduous forest plantation,North India

Forests play a vital role in biological cycles and environmental regulation. To understand the key processes of forest canopies (e.g., photosynthesis, respiration and transpiration), reliable and accurate information on spatial variability of Leaf Area Index (LAI), and its seasonal dynamics is essential. In the present study, we assessed the performance of biophysical parameter (LAI) retrieval methods viz. Look-Up Table (LUT)-inversion, MLRA-GPR (Machine Learning Regression Algorithm- Gaussian Processes Regression) and empirical models, for estimating the LAI of tropical deciduous plantation using ARTMO (Automated Radiative Transfer Models Operator) tool and Sentinel-2 satellite images. The study was conducted in Central Tarai Forest Division, Haldwani, located in the Uttarakhand state, India. A total of 49 ESUs (Elementary Sampling Unit) of 30 m × 30 m size were established based on variability in composition and age of plantation stands. In-situ LAI was recorded using plant canopy imager during the leaf growing, peak and senescence seasons. The PROSAIL model was calibrated with site-specific biophysical and biochemical parameters before used to the predicted LAI. The plantation LAI was also predicted by an empirical approach using optimally chosen Sentinel-2 vegetation indices. In addition, Sentinel-2 and MODIS LAI products were evaluated with respect to LAI measurements. MLRA-GPR offered best results for predicting LAI of leaf growing (R² = 0.9, RMSE = 0.14), peak (R² = 0.87, RMSE = 0.21) and senescence (R² = 0.86, RMSE = 0.31) seasons while LUT inverted model outperformed VI’s based parametric regression model. Vegetation indices (VIs) derived from 740 nm, 783 nm and 2190 nm band combinations of Sentinel-2 offered the best prediction of LAI.     

Impact of data model and point density on aboveground forest biomass estimation from airborne LiDAR

Background

Accurate estimation of aboveground forest biomass (AGB) and its dynamics is of paramount importance in understanding the role of forest in the carbon cycle and the effective implementation of climate change mitigation policies. LiDAR is currently the most accurate technology for AGB estimation. LiDAR metrics can be derived from the 3D point cloud (echo-based) or from the canopy height model (CHM). Different sensors and survey configurations can affect the metrics derived from the LiDAR data. We evaluate the ability of the metrics derived from the echo-based and CHM data models to estimate AGB in three different biomes, as well as the impact of point density on the metrics derived from them.

Results

Our results show that differences among metrics derived at different point densities were significantly different from zero, with a larger impact on CHM-based than echo-based metrics, particularly when the point density was reduced to 1 point m^?2. Both data models-echo-based and CHM-performed similarly well in estimating AGB at the three study sites. For the temperate forest in the Sierra Nevada Mountains, California, USA, R² ranged from 0.79 to 0.8 and RMSE (relRMSE) from 69.69 (35.59%) to 70.71 (36.12%) Mg ha^?1 for the echo-based model and from 0.76 to 0.78 and 73.84 (37.72%) to 128.20 (65.49%) Mg ha^?1 for the CHM-based model. For the moist tropical forest on Barro Colorado Island, Panama, the models gave R² ranging between 0.70 and 0.71 and RMSE between 30.08 (12.36%) and 30.32 (12.46) Mg ha^?1 [between 0.69–0.70 and 30.42 (12.50%) and 61.30 (25.19%) Mg ha^?1] for the echo-based [CHM-based] models. Finally, for the Atlantic forest in the Sierra do Mar, Brazil, R² was between 0.58–0.69 and RMSE between 37.73 (8.67%) and 39.77 (9.14%) Mg ha^?1 for the echo-based model, whereas for the CHM R² was between 0.37–0.45 and RMSE between 45.43 (10.44%) and 67.23 (15.45%) Mg ha^?1.

Conclusions

Metrics derived from the CHM show a higher dependence on point density than metrics derived from the echo-based data model. Despite the median of the differences between metrics derived at different point densities differing significantly from zero, the mean change was close to zero and smaller than the standard deviation except for very low point densities (1 point m^?2). The application of calibrated models to estimate AGB on metrics derived from thinned datasets resulted in less than 5% error when metrics were derived from the echo-based model. For CHM-based metrics, the same level of error was obtained for point densities higher than 5 points m^?2. The fact that reducing point density does not introduce significant errors in AGB estimates is important for biomass monitoring and for an effective implementation of climate change mitigation policies such as REDD + due to its implications for the costs of data acquisition. Both data models showed similar capability to estimate AGB when point density was greater than or equal to 5 point m^?2.

     

Monitoring early stage invasion of exotic Spartina alterniflora using deep-learning super-resolution techniques based on multisource high-resolution satellite imagery: A case study in the Yellow River Delta,China

Over the past decades, Spartina alterniflora, one of the top exotic invasive plants in China, has expanded throughout coastal China. In the Yellow River Delta (YRD), the rapid expansion of S. alterniflora has caused serious negative ecological effects. Current studies have concentrated primarily on mapping the distribution of S. alterniflora with medium-resolution satellite imagery at the regional or landscape scale, which have a limited capability in early detection and monitoring of the invasive process at the patch scale. In this study, we proposed a framework for monitoring the early stage invasion of S. alterniflora patches in the YRD using multiyear multisource high-spatial-resolution satellite imagery with various ground sampling distances (WorldView-2, SPOT-6, GaoFen-1, GaoFen-2, and GaoFen-6 from 2012 to 2019). First, we proposed to use deep-learning-based image super-resolution models to enhance all images to submeter (0.5 m) resolution. Then, we adopted stepwise evolution analysis-based image segmentation and object-based classification rules to detect and delineate S. alterniflora patches from the super-resolved imagery. By investigating Super-Resolution Convolutional Neural Networks (SRCNN) and Fast Super-Resolution Convolutional Neural Networks (FSRCNN) and comparing these methods with the conventional bicubic interpolation method for image resolution enhancement, we concluded that FSRCNN was superior in constructing spectral and structural details from the 1 m/1.5 m/2 m resolution images to 0.5 m resolution. FSRCNN, in particular, was more effective and efficient in discerning and estimating the size of small S. alterniflora patches (<50 m²). Using our method, 76 of 83 field-measured small patches were accurately detected and the delineated S. alterniflora patch perimeters agreed well with the field-measured patch perimeters (root mean square error [RMSE] = 8.29 m, mean absolute percentage error [MAPE] = 23.46 %). The invasion process showed fast expansion from 2012 to 2015 and slow growth from 2016 to 2019. We observed that the landward limits of S. alterniflora patches were influenced by elevation and vicinity to tidal creeks.     

Estimation of higher chlorophylla concentrations using field spectral measurement and HJ-1A hyperspectral satellite data in Dianshan Lake,China

Based on in situ water sampling and field spectral measurements in Dianshan Lake, a semi-analytical three-band algorithm was used to estimate Chlorophylla (Chla) content in case II waters. The three bands selected to estimate Chla for high concentrations included 653, 691 and 748 nm. An equation, based on the difference in reciprocal reflectance between 653 and 691 nm, multiplied by reflectance at 748 nm as [R_rs⁻¹(653) − R_rs⁻¹ (691)] R_rs(748), explained 85.57% of variance in Chla concentration with a root mean square error (RMSE) of <6.56 mg/m³. In order to test the utility of this model with satellite data, HJ-1A Hyperspectral Imager (HSI) data were analyzed using comparable wavelengths selected from the in situ data [B₆₇⁻¹(656) − B₈₀⁻¹(716)] B₈₇(753). This model accounted for 84.3% of Chla variation, estimating Chla concentrations with an RMSE of <4.23 mg/m³. The results illustrate that, based on the determined wavelengths, the spectrum-based model can achieve a high estimation accuracy and can be applied to hyperspectral satellite imagery especially for higher Chla concentration waters.     

Comparing methods for mapping canopy chlorophyll content in a mixed mountain forest using Sentinel-2 data

The Sentinel-2 Multi-Spectral Imager (MSI) has three spectral bands centered at 705, 740, and 783 nm wavelengths that exploit the red-edge information useful for quantifying plant biochemical traits. This sensor configuration is expected to improve the prediction accuracy of vegetation chlorophyll content. In this work, we assessed the performance of several statistical and physical-based methods in retrieving canopy chlorophyll content (CCC) from Sentinel-2 in a heterogeneous mixed mountain forest. Amongst the algorithms presented in the literature, 13 different vegetation indices (VIs), a non-parametric statistical approach, and two radiative transfer models (RTM) were used to assess the CCC prediction accuracy. A field campaign was conducted in July 2017 to collect in situ measurements of CCC in Bavarian forest national park, and the cloud-free Sentinel-2 image was acquired on 13 July 2017. The leave-one-out cross-validation technique was used to compare the VIs and the non-parametric approach. Whereas physical-based methods were calibrated using simulated data and validated using the in situ reference dataset. The statistical-based approaches, such as the modified simple ratio (mSR) vegetation index and the partial least square regression (PLSR) outperformed all other techniques. As such the modified simple ratio (mSR3) (665, 865) gave the lowest cross-validated RMSE of 0.21 g/m² (R² = 0.75). The PLSR resulted in the highest R² of 0.78, and slightly higher RMSE =0.22 g/m² than mSR3. The physical-based approach-INFORM inversion using look-up table resulted in an RMSE =0.31 g/m², and R² = 0.67. Although mapping CCC using these methods revealed similar spatial distribution patterns, over and underestimation of low and high CCC values were observed mainly in the statistical approaches. Further validation using in situ data from different terrestrial ecosystems is imperative for both the statistical and physical-based approaches' effectiveness to quantify CCC before selecting the best operational algorithm to map CCC from Sentinel-2 for long-term terrestrial ecosystems monitoring across the globe.     

Integrating environmental variables and WorldView-2 image data to improve the prediction and mapping of Thaumastocoris peregrinus (bronze bug) damage in plantation forests

This study integrated environmental variables together with high spectral resolution WorldView-2 imagery to detect and map Thaumastocoris peregrinus damage in Eucalypt plantation forests in KwaZulu-Natal, South Africa. The WorldView-2 bands, vegetation indices and environmental variables were entered separately into PLS regression models to predict T. peregrinus damage. The datasets were then integrated to test the collective strength in predicting T. peregrinus damage. Important variables were identified by variable importance (VIP) scores and were re-entered into a PLS regression model. The VIP model was then extrapolated to map the severity of damage and predicted T. peregrinus damage with an R² value of 0.71 and a RMSE of 3.26% on an independent test dataset. The red edge and near-infrared bands of the WorldView-2 sensor together with the temperature dataset were identified as important variables in predicting T. peregrinus damage. The results indicate the potential of integrating WorldView-2 data and environmental variables to improve the mapping and monitoring of insect outbreaks in plantation forests. The result is critical for plantation health monitoring using a new sensor which contains important vegetation wavelengths.     

Evaluation of a three-band model for estimating chlorophyll-a concentration in tidal reaches of the Pearl River Estuary, China

Accurate assessment of phytoplankton chlorophyll-a (Chla) concentration in turbid waters by means of remote sensing was challenging due to the optical complexity of turbid waters. Recently, a conceptual model containing reflectance in three spectral bands in the red and near-infrared range of the spectrum was suggested for retrieving Chla concentrations in turbid productive waters. The objective of this paper was to evaluate the performance of this three-band model to estimate Chla concentration in the Pearl River Estuary (PRE), China. Reflectance spectra of surface water and water samples were collected concurrently. The samples contained variable Chla (4.80-92.60 mg/m³) and total suspended solids (0.4-55.2 mg/L dry wt). Colored dissolved organic matter (CDOM) absorption at 400 nm was 0.40-1.41 m⁻¹; turbidity ranged from 4 to 25 NTU (Nephelometric Turbidity Units). The three-band model was spectrally calibrated by iterative and least-square linear regression methods to select the optimal spectral bands for the most accurate Chla estimation. Strong linear relationships (R²=0.81, RMSE=1.4 mg/m³, N=32) were established between measured Chla and the levels obtained from the calibrated three-band model [R⁻¹(684)-R⁻¹(690)]×R(718), where R(λ) was the reflectance at wavelength λ. The calibrated three-band model was independently validated (R²=0.9521, RMSE=6.44 mg/m³, N=16) and applied to retrieve Chla concentrations from the calibrated EO-1 Hyperion reflectance data in the PRE on December 21, 2006. The EO-1 Hyperion-derived Chla concentrations were further validated using synchronous in situ data collected on the same day (R²=0.64, RMSE=2 mg/m³, N=9). The spatial tendency of Chla distribution mapping by Hyperion showed gradually increased concentrations of Chla farther from the river mouths (although decreasing from east to west), which were disturbed by the combination of river outlets and tidal current in Lingding Bay of the PRE. This observation conformed to previous observations and studies, and could reasonably be explained by geographical changes. Also, results indicated that the slope of the three-band regression line decreased as the Chla concentration increased, resulting in the first sensitive band of the three-band model to move towards short wavelengths. These findings validated the rationale behind the conceptual model and demonstrated the robustness of this algorithm for Chla retrieval from in situ data and the Hyperion satellite sensor in turbid estuarine waters of the PRE, China.     

A radiative transfer model-based method for the estimation of grassland aboveground biomass

This paper presents a novel method to derive grassland aboveground biomass (AGB) based on the PROSAILH (PROSPECT + SAILH) radiative transfer model (RTM). Two variables, leaf area index (LAI, m²m⁻², defined as a one-side leaf area per unit of horizontal ground area) and dry matter content (DMC, gcm⁻², defined as the dry matter per leaf area), were retrieved using PROSAILH and reflectance data from Landsat 8 OLI product. The result of LAI × DMC was regarded as the estimated grassland AGB according to their definitions. The well-known ill-posed inversion problem when inverting PROSAILH was alleviated using ecological criteria to constrain the simulation scenario and therefore the number of simulated spectra. A case study of the presented method was applied to a plateau grassland in China to estimate its AGB. The results were compared to those obtained using an exponential regression, a partial least squares regression (PLSR) and an artificial neural networks (ANN). The RTM-based method offered higher accuracy (R² = 0.64 and RMSE = 42.67 gm⁻²) than the exponential regression (R² = 0.48 and RMSE = 41.65 gm⁻²) and the ANN (R² = 0.43 and RMSE = 46.26 gm⁻²). However, the proposed method offered similar performance than PLSR as presented better determination coefficient than PLSR (R² = 0.55) but higher RMSE (RMSE = 37.79 gm⁻²). Although it is still necessary to test these methodologies in other areas, the RTM-based method offers greater robustness and reproducibility to estimate grassland AGB at large scale without the need to collect field measurements and therefore is considered the most promising methodology.     

Upscaling plot-scale soil respiration in winter wheat and summer maize rotation croplands in Julu County,North China

Soil respiration (R_s) data from 45 plots were used to estimate the spatial patterns of R_s during the peak growing seasons of winter wheat and summer maize in Julu County, North China, by combining satellite remote sensing data, field-measured data, and a support vector regression (SVR) model. The observed R_s values were well reproduced by the model at the plot scale, with a root-mean-square error (RMSE) of 0.31 μmol CO₂ m⁻² s⁻¹ and a coefficient of determination (R²) of 0.73. No significant difference was detected between the prediction accuracy of the SVR model for winter wheat and summer maize. With forcing from satellite remote sensing data and gridded soil property data, we used the SVR model to predict the spatial distributions of R_s during the peak growing seasons of winter wheat and summer maize rotation croplands in Julu County. The SVR model captured the spatial variations of R_s at the county scale. The satellite-derived enhanced vegetation index was found to be the most important input used to predict R_s. Removal of this variable caused an RMSE increase from 0.31 μmol CO₂ m⁻² s⁻¹ to 0.42 μmol CO₂ m⁻² s⁻¹. Soil properties such as soil organic carbon (SOC) content and soil bulk density (SBD) were the second most important factors. Their removal led to an RMSE increase from 0.31 μmol CO₂ m⁻² s⁻¹ to 0.37 μmol CO₂ m⁻² s⁻¹. The SVR model performed better than multiple regression in predicting spatial variations of R_s in winter wheat and summer maize rotation croplands, as shown by the comparison of the R² and RMSE values of the two algorithms. The spatial patterns of R_s are better captured using the SVR model than performing multiple regression, particularly for the relatively high and relatively low R_s values at the center and northeast study areas. Therefore, SVR shows promise for predicting spatial variations of R_s values on the basis of remotely sensed data and gridded soil property data at the county scale.     

Estimation of Seasonal Sun-Induced Fluorescence Dynamics of Indian Tropical Deciduous Forests using SCOPE and Sentinel-2 MSI

This paper presents a new approach to estimate spatial Sun-Induced Fluorescence (SIF) using the empirical relationship between simulated Canopy Chlorophyll Concentration (CCC) and simulated SIF. PROSAIL model [PROpriétésSPECTrales (PROSPECT) and Scattering by Arbitrarily Inclined Leaves (SAIL) models] was used to simulate CCC. CCC maps were generated through an Automated Radiative Transfer Model Operator (ARTMO) using the PROSAIL model and Sentinel-2 Multi-Spectral Imager (MSI) imagery. The Soil Canopy Observation, Photochemistry, and Energy fluxes (SCOPE) model was used to simulate SIF emitted at 740 nm (SIF₇₄₀), at 760 nm (SIF₇₆₀), and top of canopy (SIF_TOC) (640-850 nm). The SCOPE model, configured with the specification of the Sentinel-2 sensor, simulates SIF within the spectrum range of 640-850 nm. A non-linear logarithmic relationship (R²>0.9, p < 0.05) was observed between simulated SIF and simulated CCC. Simulated CCC was linearly related to observed CCC with R² 0.88, 0.92 and 0.89 and RMSE = 0.04, 0.17 and 0.09 gm/m² at p < 0.05 for summer, post-monsoon and early winter respectively. Whereas, the simulated CCC did not capture the full range of CCC variability for the post-monsoon season. Simulated SIF (SIF₇₆₀) was well correlated with SIF from Orbiting Carbon Observatory-2 (OCO-2) satellite with R² 0.68, 0.73 and 0.73 (RMSE = <1 W/m²/sr/μm, p < 0.05) for the month of summer (April), pre-monsoon (May) and early winter season (November) respectively. Temporal SIF_TOC effectively captured the seasonal variability associated with the phenology of deciduous tree species. Among various Sentinel-2 MSI derived VIs, Red Edge NDVI (RENDVI) exhibited maximum sensitivity with SIF (highest monthly average R²> 0.6, p < 0.05). The spatial SIF would serve as an useful link between airborne /satellite derived SIF and in-situ fluorescence measurements to understand multiscale SIF variability of terrestrial vegetation.