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.     

RTM-based dynamic absorption integrals for the retrieval of biochemical vegetation traits

Information about pigment and water contents provides comprehensive insights for evaluating photosynthetic potential and activity of agricultural crops. In this study, we present the concept of using spectral integral ratios (SIR) to retrieve three biochemical traits, namely chlorophyll a and b (C_ab), carotenoids (C_cx), and water (C_w) content, simultaneously from hyperspectral measurements in the wavelength range 460−1100 nm. The SIR concept is based on automatic separation of respective absorption features through local peak and intercept analysis between log-transformed reflectance and convex hulls. The algorithm was tested on two synthetically established databases using a physiologically constrained look-up-table (LUT) generated by (i) the leaf optical properties model PROSPECT and (ii) the canopy radiative transfer model (RTM) PROSAIL. LUT constraints were realized based on natural C_cx-C_ab relations and green peak locations identified in the leaf optical database ANGERS. Linear regression between obtained SIRs and model parameters resulted in coefficients of determination (R²) of 0.66 (i and ii) for C_cx, R² = 0.85 (i) and 0.53 (ii) for C_ab, and R² = 0.97 (i) and 0.67 (ii) for C_w, respectively. Using the model established from the PROSPECT LUT, leaf level validation was carried out based on ANGERS data with reasonable results both in terms of goodness of fit and root mean square error (RMSE) (C_cx: R² = 0.86, RMSE = 2.1 μg cm⁻²; C_ab: R² = 0.67, RMSE = 12.5 μg cm^-2; C_w: R² = 0.89, RMSE = 0.007 cm). The algorithm was applied to airborne spectrometric HyMap data acquired on 12th July 2003 in Barrax, Spain and to AVIRIS-NG data recorded on 2nd July 2018 southwest of Munich, Germany. Mapping of the SIR results as multiband images (3-segment SIR) allows for intuitive visualization of dominant absorptions with respect to the three considered biochemical variables. Barrax in situ validation using linear regression models derived from PROSAIL LUT showed satisfactory results regarding C_ab (R² = 0.84; RMSE = 9.06 μg cm^-2) and canopy water content (CWC, R² = 0.70; RMSE = 0.05 cm). Retrieved C_cx values were reasonable according to C_ab-C_cx-dependence plausibility analysis. Hence, the presented SIR algorithm allows for computationally efficient and RTM supported robust retrievals of the two most important vegetation pigments as well as of water content and is ready to be applied on satellite imaging spectroscopy data available in the near future. The algorithm is publicly available as an interface supported tool within the 'Agricultural Applications' of the EnMAP-Box 3 hyperspectral remote sensing software suite.     

用地基激光雷达提取单木结构参数——以白皮松为例

以白皮松(Pinus bungeana Zucc)为研究对象,针对地基激光雷达TLS扫描的3维点云数据在单株木垂直方向的分布特征,提出了一种基于体元化方法的树干覆盖度变化检测方法,获取单木枝下高;然后根据获取的枝下高引入2维凸包算法获取垂直方向分层树冠轮廓,并计算树冠体积和冠幅;同时获取的单木参数还有胸径与树高。结果表明:单木枝下高的估测精度较高,R2与RMSE分别为0.97 m和0.21 m;胸径估测结果的R2与RMSE分别为0.79 cm和1.07 cm;采用逐步线性回归方法建立单木树冠体积与其他单木参数的相关关系,模型变量包括冠幅、叶子填充树冠长度和胸径,样本数为20,模型的R2与RMSE分别是0.967 m3和2.64 m3。本文方法能较准确地估测枝下高,TLS数据具有对树冠结构3维建模的潜力。     

Spectral behavior of Persian oak under compound stress of water deficit and dust storm

Persian oak (Quercus Brantii Lindl.) which is the most widely distributed tree in the Zagros Mountain forests is affected by western dust storms, mostly originating in Iraq, and harsh water stress as well. The objective of this research is to analyze the spectral behavior of Persian oak under water and dust stress scenarios, aiming to pave the way for modeling the stresses of drought and dust storms on oak trees using remote sensing images. Experiments were carried out on 54 two-year old oak tree seedlings, using a portable wind tunnel in greenhouse conditions. Water stress was induced on seedlings by means of changes in irrigation practices, i.e. well-watered (100 % field capacity), medium water deficit condition (40 % field capacity), and severe water deficit condition (20 % field capacity) treatments. Dust stress is also investigated by using three different dust particle concentrations, i.e. 350, 750 and 1500 (μg/m³). The spectrometry experiments were carried out at leaf and canopy levels in dark room by Fieldspec-3-ASD spectrometer. Spectral analysis was conducted using four procedures: (i) narrow-band spectral indices analysis, (ii) geometric indicators extraction from absorption features, (iii) Partial Least Squares Regression (PLSR), and SVM classifier. Results show that water stress could be modeled much better using PLSR statistic (R² = 0.87, RMSE = 0.12), narrow-band indices analysis (R²_cv = 0.75, RMSE_cv = 0.17), and continuum removal (R² = 0.71, RMSE = 0.20), respectively. For dust stress, PLSR (R² = 0.83, RMSE = 0.14) and narrow-band indices (R² _cv = 0.7, RMSE _cv = 0.30) showed the best results, respectively. SVM could successfully separate stressed and not-stressed samples and also the stress types at both leaf and canopy levels, but it could not distinguish the different levels of stresses.     

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.     

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.     

Soil respiration mapped by exclusively use of MODIS data for forest landscapes of Saskatchewan,Canada

Soil respiration (R_s) is of great importance to the global carbon balance. Remote sensing of R_s is challenging because of (1) the lack of long-term R_s data for model development and (2) limited knowledge of using satellite-based products to estimate R_s. Using 8-years (2002–2009) of continuous R_s measurements with nonsteady-state automated chamber systems at a Canadian boreal black spruce stand (SK-OBS), we found that R_s was strongly correlated with the product of the normalized difference vegetation index (NDVI) and the nighttime land surface temperature (LSTn) derived from Moderate Resolution Imaging Spectroradiometer (MODIS) imagery. The coefficients of the linear regression equation of this correlation between R_s and NDVI × LSTn could be further calibrated using the MODIS leaf area index (LAI) product, resulting in an algorithm that is driven solely by remote sensing observations. Modeled R_s closely tracked the seasonal patterns of measured R_s and explained 74–92% of the variance in R_s with a root mean square error (RMSE) less than 1.0 g C/m²/d. Further validation of the model from SK-OBS site at another two independent sites (SK-OA and SK-OJP, old aspen and old jack pine, respectively) showed that the algorithm can produce good estimates of R_s with an overall R² of 0.78 (p < 0.001) for data of these two sites. Consequently, we mapped R_s of forest landscapes of Saskatchewan using entirely MODIS observations for 2003 and spatial and temporal patterns of R_s were well modeled. These results point to a strong relationship between the soil respiratory process and canopy photosynthesis as indicated from the greenness index (i.e., NDVI), thereby implying the potential of remote sensing data for detecting variations in R_s. A combination of both biological and environmental variables estimated from remote sensing in this analysis may be valuable in future investigations of spatial and temporal characteristics of R_s.     

Remote detection of bare soil moisture using a surface-temperature-based soil evaporation transfer coefficient

An approach for estimating soil moisture is presented and tested by using surface-temperature-based soil evaporation transfer coefficient (h_a), a coefficient recently proposed through the equation h_a = (T_s − T_a)/(T_sd − T_a), where T_s, T_sd, and T_a are land surface temperature (LST), reference soil (dry soil without evaporation) surface temperature, and air temperature respectively. Our analysis and controllable experiment indicated that h_a closely related to soil moisture, and therefore, a relationship between field soil moisture and h_a could be developed for soil moisture estimation. Field experiments were carried out to test the relationship between h_a and soil moisture. Time series Aqua-MODIS images were acquired between 11 Sep. 2006 and 1 Nov. 2007. Then, MODIS derived h_a and simultaneous measured soil moisture for different soil depths were used to establish the relations between the two variables. Results showed that there was a logarithmic relationship between soil moisture and h_a (P < 0.01). These logarithmic models were further validated by introducing another ground-truth data gathered from 46 meteorological stations in Hebei Province. Good agreement was observed between the measured and estimated soil moisture with RMSE of 0.0374 cm³/cm³ and 0.0503 cm³/cm³ for surface energy balance method at two soil depths (10 cm and 20 cm), with RMSE of 0.0467 cm³/cm³ and 0.0581 cm³/cm³ for maximum temperature method at two soil depths. For vegetated surfaces, the ratio of h_a and NDVI suggested to be considered. The proposed approach has a great potential for soil moisture and drought evaluation by remote sensing.     

闪电河流域农牧交错带微波遥感土壤水分产品评价

空间网格分辨率为9 km 的 SMAP(Soil Moisture Active and Passive)、0.1D(Degree)的 ASCAT(The Advanced Scatterometer)、25 km 的 FY-3B 以及25 km ESA-CCI(European Space Agency-Climat...     

Different units of measurement of carotenoids estimation in cotton using hyperspectral indices and partial least square regression

This paper examines the use of canopy reflectance for different units of measurements of carotenoids estimation. Field spectral measurements were collected over cotton in different intensive field campaigns organized during the growing seasons of 2010 and 2011. Three units of measurement were evaluated carotenoids expressed as a mass per unit soil surface area (g/m²), a mass per unit leaf area (μg/cm²), and a mass per unit fresh leaf weight (mg/g), respectively. Four methods were compared to retrieve amount of carotenoids: stepwise multiple linear regression (SMLR), published spectral indices, band combination indices, and partial least square regression (PLSR). Results show that maximum sensitivity of reflectance to variation in different units of measurement of carotenoids was found in the green region at 515–550 nm, and at 715 nm and 750 nm regions in the far-red wavelengths. The predictive accuracies of Car (g/m²), Car (μg/cm²) and Car (mg/g) were tested on a validation data set and the results show that the highest R² values between estimations and observations were 0.468 for Car (g/m²), 0.563 for Car (μg/cm²), and 0.456 for Car (mg/g), with relative root mean square error (RMSE%, RMSE/mean) of 48.72%, 22.07% and 21.07%, respectively. Compared to Car (g/m²) and Car (mg/g), the model performance indices for Car (μg/cm²) show a high degree of consistency among the R² values and RMSE% and MAE% values. Further comparison were performed among the estimation accuracies of different unit carotenoids and among the different approaches used in the study by a paired-t-test. The results indicate that although the best estimation results for Car (μg/cm²) and Car (mg/g) were both obtained based on PLSR, they can be estimated by all four adopted methods without significant differences (P > 0.1). Whereas for Car (g/m²), the best estimation results were obtained based on published vegetation indices CIred-edge, which were significantly better than the estimation results based on SMLR (P < 0.000). In summary, the results of this study show that even the carotenoids expressed on concentration (mg/g) or content (μg/cm²) basis at leaf level can be estimated with the same prediction accuracies to the carotenoids expressed as a mass per unit surface area (g/m²) at canopy level using reflectance measurement at canopy level.     

Machine learning approaches to coastal water quality monitoring using GOCI satellite data

Since coastal waters are one of the most vulnerable marine systems to environmental pollution, it is very important to operationally monitor coastal water quality. This study attempts to estimate two major water quality indicators, chlorophyll-a (chl-a) and suspended particulate matter (SPM) concentrations, in coastal environments on the west coast of South Korea using Geostationary Ocean Color Imager (GOCI) satellite data. Three machine learning approaches including random forest, Cubist, and support vector regression (SVR) were evaluated for coastal water quality estimation. In situ measurements (63 samples) collected during four days in 2011 and 2012 were used as reference data. Due to the limited number of samples, leave-one-out cross validation (CV) was used to assess the performance of the water quality estimation models. Results show that SVR outperformed the other two machine learning approaches, yielding calibration R² of 0.91 and CV root-mean-squared-error (RMSE) of 1.74 mg/m³ (40.7%) for chl-a, and calibration R² of 0.98 and CV RMSE of 11.42 g/m³ (63.1%) for SPM when using GOCI-derived radiance data. Relative importance of the predictor variables was examined. When GOCI-derived radiance data were used, the ratio of band 2 to band 4 and bands 6 and 5 were the most influential input variables in predicting chl-a and SPM concentrations, respectively. Hourly available GOCI images were useful to discuss spatiotemporal distributions of the water quality parameters with tidal phases in the west coast of Korea.     

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

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

机载激光雷达及高光谱的森林乔木物种多样性遥感监测

利用机载LiDAR和高光谱数据并结合37个地面调查样本数据,基于结构差异与光谱变异理论,通过相关分析法分别筛选了3个最优林冠结构参数和6个最优光谱指数,在单木尺度上利用自适应C均值模糊聚类算法,在神农架国家自然保护区开展森林乔木物种多样性监测,实现了森林乔木物种多样性的区域成图。研究结果表明,(1)基于结合形态学冠层控制的分水岭算法可以获得较高精度的单木分割结果(R~2=0.88,RMSE=13.17,P0.001);(2)基于LiDAR数据提取的9个结构参数中,95%百分位高度、冠层盖度和植被穿透率为最优结构参数,与Shannon-Wiener指数的相关性达到R~2=0.39—0.42(P0.01);(3)基于机载高光谱数据筛选的16个常用的植被指数中,CRI、OSAVI、Narrow band NDVI、SR、Vogelmann index1、PRI与Shannon-Wiener指数的相关性最高(R~2=0.37—0.45,P0.01);(4)在研究区,利用以30 m×30 m为窗口的自适应模糊C均值聚类算法可预测的最大森林乔木物种数为20,物种丰富度的预测精度为R~2=0.69,RMSE=3.11,Shannon-Wiener指数的预测精度为R~2=0.70,RMSE=0.32。该研究在亚热带森林开展乔木物种多样性监测,是在区域尺度上进行物种多样性成图的重要实践,可有效补充森林生物多样性本底数据的调查手段,有助于实现生物多样性的长期动态监测及科学分析森林物种多样性的现状和变化趋势。     

光学与微波植被指数协同反演农作物叶面积指数的可行性分析

光学遥感是目前反演植被叶面积指数LAI（Leaf Area Index）的主要手段,但是当叶面积指数较大时存在光学遥感信息饱和、反演精度显著降低的问题。叶面积指数和平均叶倾角对光学、微波波段范围内反射和散射特性都有重要影响,主要表现在植被结构参数的变化可以引起冠层孔隙率和消光截面大小的改变。本文以典型农作物玉米为例,通过构建统一的PROSAIL和MIMICS模型输入参数,生成一套玉米全生长期光学二向反射率和全极化微波后向散射系数模拟库和冠层参数库。通过对模拟数据与LAI敏感性和相关性分析得出：（1）光学植被指数MNDVI_{(800 nm,2000 nm)},在LAI为0—3时敏感,基于MNDVI与LAI的回归模型可以估算LAI变化 0.4的情况,RMSE是0.33,R²是0.958。（2）微波植被指数SAR_{SRVI(1.4 GHz HH,9.6 GHz HV)},在LAI为3—6时敏感,基于SAR_SRVI与LAI的回归模型可以估算LAI变化1的情况,RMSE为0.22,R²是0.9839。研究表明,采用分段敏感的植被指数,协同光学和微波遥感反演玉米全生长期叶面积指数是可行的。     

Mapping forest aboveground biomass in the reforested Buffelsdraai landfill site using texture combinations computed from SPOT-6 pan-sharpened imagery

Developing models for estimating aboveground biomass (AGB) in naturally growing forests is critical for climate change modelling. AGB models developed using satellite imagery varies with study area, depending on the complexity of vegetation and landscape structure, which affects the upwelling radiance. We assessed the potential of SPOT-6 imagery in predicting AGB of trees planted at different time periods, using image texture combinations. Image texture variables were computed from the SPOT6 pan-sharpened image data, which is characterised by a 1.5 m spatial resolution. In addition, we incorporated the minimal variance technique to select the optimum window sizes that best captures AGB variation in our study area. The results showed that image texture was able to detect AGB for both mature and young trees, however, models detecting mature trees were more superior, with accuracies of R² = 0.70 and 0.25 for 2009–2011 and 2011–2013 plantation phases, respectively. In addition, our results showed that the three band texture ratios yielded the highest accuracy (R² = 0.88 and RMSE = 54.54 kg m⁻²) compared to two texture (R² = 0.85 and RMSE = 60.65 kg m⁻²) and single texture band combinations (R² = 0.64 and RMSE = 94.13 kg m⁻²). A frequency analysis was also run to determine which bands appeared more frequently in the selected texture band models. The frequency analysis revealed that both the red and green bands appeared more frequently on the selected texture band variables, indicating that they were more sensitive to the variation of AGB in our study area. The results showed high variation in AGB within the Buffelsdraai reforestation site, especially due to varying tree plantation phases as well as topography. In essence, the study demonstrated the possibility of image texture combinations computed from the SPOT-6 image in estimating AGB.     

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.     

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.

     

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.     

Estimating net primary productivity in tropical forest plantations in India using satellite-driven ecosystem model

Net Primary Productivity (NPP) is a significant biophysical vegetation variable to understand the spatio-temporal distribution of carbon and source-sink nature of the ecosystem. This study was carried out in a forest plantation area and aimed to (i) estimate the spatio-temporal patterns of NPP during 2009 and 2010 using Carnegie-Ames-Stanford Approach [CASA] model and (ii) study the effects of climate variables on the NPP using generalized linear modelling (GLM) approach. The total annual NPP varied from 157.21 to 1030.89 gC m^?2 yr^?1 for the year 2009 and from 154.36 to 1124.85 g C m^?2 yr^?1 for the year 2010. The annual NPP was assessed across four major plantation types, where maximum NPP gain (106 and 139 g C m^?2 yr^?1 ) in October was noticed in teak (Tectona grandis) and minimum (77 and 109 g C m^?2 yr^?1 ) in eucalyptus (Eucalyptus hybrid) during 2009 and 2010.The validation, using field-estimated NPP, showed under-estimation of modelled NPP, with maximum MAPE of 34% for eucalyptus and minimum of 13% for teak. The dominant influence of precipitation on the NPP was revealed by GLM explaining more than 20% of variation. CASA model efficiently estimated the annual NPP of plantations. The accuracy could be improved further with inclusion of higher resolution data.     

Mapping grassland leaf area index with airborne hyperspectral imagery: A comparison study of statistical approaches and inversion of radiative transfer models

Statistical and physical models have seldom been compared in studying grasslands. In this paper, both modeling approaches are investigated for mapping leaf area index (LAI) in a Mediterranean grassland (Majella National Park, Italy) using HyMap airborne hyperspectral images. We compared inversion of the PROSAIL radiative transfer model with narrow band vegetation indices (NDVI-like and SAVI2-like) and partial least squares regression (PLS). To assess the performance of the investigated models, the normalized RMSE (nRMSE) and R² between in situ measurements of leaf area index and estimated parameter values are reported. The results of the study demonstrate that LAI can be estimated through PROSAIL inversion with accuracies comparable to those of statistical approaches (R² = 0.89, nRMSE = 0.22). The accuracy of the radiative transfer model inversion was further increased by using only a spectral subset of the data (R² = 0.91, nRMSE = 0.18). For the feature selection wavebands not well simulated by PROSAIL were sequentially discarded until all bands fulfilled the imposed accuracy requirements.