Generation and evaluation of gross primary productivity using Landsat data through blending with MODIS data

Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) has been used for the blending of Landsat and MODIS data. Specifically, the 30 m Landsat-7 ETM+ (Enhanced Thematic Mapper plus) surface reflectance was predicted for a period of 10 years (2000–2009) as the product of observed ETM+ and MODIS surface reflectance (MOD09A1) on the predicted and observed ETM+ dates. A pixel based analysis for six observed ETM+ dates covering winter and summer crops showed that the prediction method was more accurate for NIR band (mean r² = 0.71, p ≤ 0.01) compared to green band (mean r² = 0.53; p ≤ 0.01). A recently proposed chlorophyll index (CI), which involves NIR and green spectral bands, was used to retrieve gross primary productivity (GPP) as the product of CI and photosynthetic active radiation (PAR). The regression analysis of GPP derived from closet observed and synthetic ETM+ showed a good agreement (r² = 0.85, p ≤ 0.01 and r² = 0.86, p ≤ 0.01) for wheat and sugarcane crops, respectively. The difference between the GPP derived from synthetic and observed ETM+ (prediction residual) was compared with the difference in GPP values from observed ETM+ on the two dates (temporal residual). The prediction residuals (mean value of 1.97 g C/m² in 8 days) was found to be significantly lower than the temporal residuals (mean value of 4.46 g C/m² in 8 days) that correspondence to 12% and 27%, respectively, of GPP values (mean value of 16.53 g C/m² in 8 days) from observed ETM+ data, implying that the prediction method was better than temporal pixel substitution. Investigating the trend in synthetic ETM+ GPP values over a growing season revealed that phenological patterns were well captured for wheat and sugarcane crops. A direct comparison between the GPP values derived from MODIS and synthetic ETM+ data showed a good consistency of the temporal dynamics but a systematic error that can be read as bias (MODIS GPP over estimation). Further, the regression analysis between observed evapotranspiration and synthetic ETM+ GPP showed good agreement (r² = 0.66, p ≤ 0.01).     

Exploring hyperspectral bands and estimation indices for leaf nitrogen accumulation in wheat

Hyperspectral sensing can provide an effective means for fast and non-destructive estimation of leaf nitrogen (N) status in crop plants. The objectives of this study were to design a new method to extract hyperspectral spectrum information, to explore sensitive spectral bands, suitable bandwidth and best vegetation indices based on precise analysis of ground-based hyperspectral information, and to develop regression models for estimating leaf N accumulation per unit soil area (LNA, g N m⁻²) in winter wheat (Triticum aestivum L.). Three field experiments were conducted with different N rates and cultivar types in three consecutive growing seasons, and time-course measurements were taken on canopy hyperspectral reflectance and LNA under the various treatments. Then, normalized difference spectral indices (NDSI) and ratio spectral indices (RSI) based on the original spectrum and the first derivative spectrum were constructed within the range of 350–2500 nm, and their relationships with LNA were quantified. The results showed that both LNA and canopy hyperspectral reflectance in wheat changed with varied N rates, with consistent patterns across different cultivars and seasons. The sensitive spectral bands for LNA existed mainly within visible and near infrared regions. The best spectral indices for estimating LNA in wheat were found to be NDSI (R₈₆₀, R₇₂₀), RSI (R₉₉₀, R₇₂₀), NDSI (FD₇₃₆, FD₅₂₆) and RSI (FD₇₂₅, FD₅₁₆), and the regression models based on the above four spectral indices were formulated as Y = 26.34x^1.887, Y = 5.095x − 6.040, Y = 0.609 e^3.008x and Y = 0.388x^1.260, respectively, with R² greater than 0.81. Furthermore, expanding the bandwidth of NDSI (R₈₆₀, R₇₂₀) and RSI (R₉₉₀, R₇₂₀) from 1 nm to 100 nm at 1 nm interval produced the LNA monitoring models with similar performance within about 33 nm and 23 nm bandwidth, respectively, over which the statistical parameters of the models became less stable. From testing of the derived equations, the model for LNA estimation on NDSI (R₈₆₀, R₇₂₀), RSI (R₉₉₀, R₇₂₀), NDSI (FD₇₃₆, FD₅₂₆) and RSI (FD₇₂₅, FD₅₁₆) gave R² over 0.79 with more satisfactory performance than previously reported models and physical models in wheat. It can be concluded that the present hyperspectral parameters of NDSI (R₈₆₀, R₇₂₀), RSI (R₉₉₀, R₇₂₀), NDSI (FD₇₃₆, FD₅₂₆) and RSI (FD₇₂₅, FD₅₁₆) can be reliably used for estimating LNA in winter wheat.     

Measurement and Scaling of Carbon Dioxide (CO<Subscript>2</Subscript>) Exchanges in Wheat Using Flux-Tower and Remote Sensing

The present study investigates the characteristics of CO₂ exchange (photosynthesis and respiration) over agricultural site dominated by wheat crop and their relationship with ecosystem parameters derived from MODIS. Eddy covariance measurement of CO₂ and H₂O exchanges was carried out at 10 Hz interval and fluxes of CO₂ were computed at half-hourly time steps. The net ecosystem exchange (NEE) was partitioned into gross primary productivity (GPP) and ecosystem respiration (R _e) by taking difference between day-time NEE and respiration. Time-series of daily reflectance and surface temperature products at varying resolution (250–1000 m) were used to derive ecosystem variables (EVI, NDVI, LST). Diurnal pattern in Net ecosystem exchange reveals negative NEE during day-time representing CO₂ uptake and positive during night as release of CO₂. The amplitude of the diurnal variation in NEE increased as LAI crop growth advances and reached its peak around the anthesis stage. The mid-day uptake during this stage was around 1.15 mg CO₂ m⁻² s⁻¹ and night-time release was around 0.15 mg CO₂ m⁻² s⁻¹. Linear and non-linear least square regression procedures were employed to develop phenomenological models and empirical fits between flux tower based GPP and NEE with satellite derived variables and environmental parameters. Enhanced vegetation index was found significantly related to both GPP and NEE. However, NDVI showed little less significant relationship with both GPP and NEE. Furthemore, temperature-greenness (TG) model combining scaled EVI and LST was parameterized to estimate daily GPP over dominantly wheat crop site. (R ² = 0.77). Multi-variate analysis shows that inclusion of LST or air temperature with EVI marginally improves variance explained in daily NEE and GPP.     

Estimation of vegetation LAI from hyperspectral reflectance data: Effects of soil type and plant architecture

The retrieval of canopy biophysical variables is known to be affected by confounding factors such as plant type and background reflectance. The effects of soil type and plant architecture on the retrieval of vegetation leaf area index (LAI) from hyperspectral data were assessed in this study. In situ measurements of LAI were related to reflectances in the red and near-infrared and also to five widely used spectral vegetation indices (VIs). The study confirmed that the spectral contrast between leaves and soil background determines the strength of the LAI–reflectance relationship. It was shown that within a given vegetation species, the optimum spectral regions for LAI estimation were similar across the investigated VIs, indicating that the various VIs are basically summarizing the same spectral information for a given vegetation species. Cross-validated results revealed that, narrow-band PVI was less influenced by soil background effects (0.15 ≤ RMSE_cv ≤ 0.56). The results suggest that, when using remote sensing VIs for LAI estimation, not only is the choice of VI of importance but also prior knowledge of plant architecture and soil background. Hence, some kind of landscape stratification is required before using hyperspectral imagery for large-scale mapping of vegetation biophysical variables.     

Terrestrial laser scanning to estimate plot-level forest canopy fuel properties

This paper evaluates the potential of a terrestrial laser scanner (TLS) to characterize forest canopy fuel characteristics at plot level. Several canopy properties, namely canopy height, canopy cover, canopy base height and fuel strata gap were estimated. Different approaches were tested to avoid the effect of canopy shadowing on canopy height estimation caused by deployment of the TLS below the canopy. Estimation of canopy height using a grid approach provided a coefficient of determination of R² = 0.81 and an RMSE of 2.47 m. A similar RMSE was obtained using the 99th percentile of the height distribution of the highest points, representing the 1% of the data, although the coefficient of determination was lower (R² = 0.70). Canopy cover (CC) was estimated as a function of the occupied cells of a grid superimposed upon the TLS point clouds. It was found that CC estimates were dependent on the cell size selected, with 3 cm being the optimum resolution for this study. The effect of the zenith view angle on CC estimates was also analyzed. A simple method was developed to estimate canopy base height from the vegetation vertical profiles derived from an occupied/non-occupied voxels approach. Canopy base height was estimated with an RMSE of 3.09 m and an R² = 0.86. Terrestrial laser scanning also provides a unique opportunity to estimate the fuel strata gap (FSG), which has not been previously derived from remotely sensed data. The FSG was also derived from the vegetation vertical profile with an RMSE of 1.53 m and an R² = 0.87.     

Applying different inversion techniques to retrieve stand variables of summer barley with PROSPECT + SAIL

This study describes the retrieval of state variables (LAI, canopy chlorophyll, water and dry matter contents) for summer barley from airborne HyMap data by means of a canopy reflectance model (PROSPECT + SAIL). Three different inversion techniques were applied to explore the impact of the employed method on estimation accuracies: numerical optimization (downhill simplex method), a look-up table (LUT) and an artificial neural network (ANN) approach. By numerical optimization (Num Opt), reliable estimates were obtained for LAI and canopy chlorophyll contents (LAI × C_ab) with r² of 0.85 and 0.94 and RDP values of 1.81 and 2.65, respectively. Accuracies dropped for canopy water (LAI × C_w) and dry matter contents (LAI × C_m). Nevertheless, the range of leaf water contents (C_w) was very narrow in the studied plant material. Prediction accuracies generally decreased in the order Num Opt > LUT > ANN. This decrease in accuracy mainly resulted from an increase in offset in the obtained values, as the retrievals from the different approaches were highly correlated. The same decreasing order in accuracy was found for the difference between the measured spectra and those reconstructed from the retrieved variable values. The parallel application of the different inversion techniques to one collective data set was helpful to identify modelling uncertainties, as shortcomings of the retrieval algorithms themselves could be separated from uncertainties in model structure and parameterisation schemes.     

陆地总初级生产力遥感估算精度分析

准确估算陆地总初级生产力GPP(Gross Primary Productivity)数值对碳循环过程模拟有重要影响。本文介绍了多种基于植被指数以及基于光能利用率的遥感GPP算法,综述了不同算法在其研究区域的估算精度;并分析了MODIS/GPP以及BESS/GPP两种遥感GPP产品在不同植被类型的估算精度。通过对比全球碳通量站网络GPP数据表明,MODIS/GPP产品在全球估算结果具显著相关性(R2=0.59)及中等标准误差(RMSE=2.86 g C/m2/day),估算精度较高的植被类型有落叶阔叶林,草地等;估算精度较低类型包括常绿阔叶林,稀树草原等。本文对GPP产品中存在的不确定性进行分析,通过综述前人研究中发现的遥感估算GPP方法中存在的问题,指出可能的提高卫星遥感GPP产品估算精度的方法及发展趋势。     

Assessment of RapidEye vegetation indices for estimation of leaf area index and biomass in corn and soybean crops

Leaf area index (LAI) and biomass are important indicators of crop development and the availability of this information during the growing season can support farmer decision making processes. This study demonstrates the applicability of RapidEye multi-spectral data for estimation of LAI and biomass of two crop types (corn and soybean) with different canopy structure, leaf structure and photosynthetic pathways. The advantages of Rapid Eye in terms of increased temporal resolution (∼daily), high spatial resolution (∼5 m) and enhanced spectral information (includes red-edge band) are explored as an individual sensor and as part of a multi-sensor constellation. Seven vegetation indices based on combinations of reflectance in green, red, red-edge and near infrared bands were derived from RapidEye imagery between 2011 and 2013. LAI and biomass data were collected during the same period for calibration and validation of the relationships between vegetation indices and LAI and dry above-ground biomass. Most indices showed sensitivity to LAI from emergence to 8 m²/m². The normalized difference vegetation index (NDVI), the red-edge NDVI and the green NDVI were insensitive to crop type and had coefficients of variations (CV) ranging between 19 and 27%; and coefficients of determination ranging between 86 and 88%. The NDVI performed best for the estimation of dry leaf biomass (CV = 27% and r² = 090) and was also insensitive to crop type. The red-edge indices did not show any significant improvement in LAI and biomass estimation over traditional multispectral indices. Cumulative vegetation indices showed strong performance for estimation of total dry above-ground biomass, especially for corn (CV ≤ 20%). This study demonstrated that continuous crop LAI monitoring over time and space at the field level can be achieved using a combination of RapidEye, Landsat and SPOT data and sensor-dependant best-fit functions. This approach eliminates/reduces the need for reflectance resampling, VIs inter-calibration and spatial resampling.     

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.     

New Triangle Vegetation Indices for Estimating Leaf Area Index on Maize

Vegetation index-based methods have been widely used to determine the leaf area index (LAI). Nevertheless, under the high canopy coverage, the estimation ability of current inversion models has been profoundly decreased, due to the “saturation” phenomenon. In this study, the LAI of maize was investigated under various growth conditions. Two new triangular vegetation indices were proposed to improve the inversion ability and estimation accuracy of LAI on maize. The triangle difference vegetation index (TDVI) and triangle ratio vegetation index (TRVI) were constructed, and their accuracies were compared with the present spectral vegetation index models. The result shows that TDVI and TRVI are highly linearly correlated with LAI. The coefficients of determination (R²) and root-mean-square errors are, respectively, 0.92 and 0.94, and 1.42 and 0.92 using the simulated data, while they are, respectively, 0.83 and 0.77, and 0.98 and 1.05 using the measured data. In comparison with other vegetation indices (e.g. MSR, MTVI2, RTVI), TDVI is better able to estimate the LAI of maize. Conversely, TRVI has better inversion ability when the LAI is more than 3. Overall, TDVI is an accurate and robust approach for estimating the LAI of maize. The proposed TDVI and TRVI can be jointly used to retrieve LAI at various canopy coverages.     

Comparative analysis of different retrieval methods for mapping grassland leaf area index using airborne imaging spectroscopy

Fine scale maps of vegetation biophysical variables are useful status indicators for monitoring and managing national parks and endangered habitats. Here, we assess in a comparative way four different retrieval methods for estimating leaf area index (LAI) in grassland: two radiative transfer model (RTM) inversion methods (one based on look-up-tables (LUT) and one based on predictive equations) and two statistical modelling methods (one partly, the other entirely based on in situ data). For prediction, spectral data were used that had been acquired over Majella National Park in Italy by the airborne hyperspectral HyMap instrument. To assess the performance of the four investigated models, the normalized root mean squared error (nRMSE) and coefficient of determination (R²) between estimates and in situ LAI measurements are reported (n = 41). Using a jackknife approach, we also quantified the accuracy and robustness of empirical models as a function of the size of the available calibration data set. The results of the study demonstrate that the LUT-based RTM inversion yields higher accuracies for LAI estimation (R² = 0.91, nRMSE = 0.18) as compared to RTM inversions based on predictive equations (R² = 0.79, nRMSE = 0.38). The two statistical methods yield accuracies similar to the LUT method. However, as expected, the accuracy and robustness of the statistical models decrease when the size of the calibration database is reduced to fewer samples. The results of this study are of interest for the remote sensing community developing improved inversion schemes for spaceborne hyperspectral sensors applicable to different vegetation types. The examples provided in this paper may also serve as illustrations for the drawbacks and advantages of physical and empirical models.     

Spectral mixture analysis to monitor defoliation in mixed-aged Eucalyptus globulus Labill plantations in southern Australia using Landsat5-TM and EO-1Hyperion data

Defoliation is a key parameter of forest health and is associated with reduced productivity and tree mortality. Assessing the health of forests requires regular observations over large areas. Satellite remote sensing provides a cost-effective alternative to traditional ground-based assessment of forest health, but assessing defoliation can be difficult due to mixed pixels where vegetation cover is low or fragmented. In this study we apply a novel spectral unmixing technique, referred to as weighted Multiple Endmember Spectral Mixture Analysis (wMESMA), to Landsat 5-TM and EO-1 Hyperion data acquired over a Eucalyptus globulus (Labill.) plantation in southern Australia. This technique combines an iterative mixture analysis cycle allowing endmembers to vary on a per pixel basis (MESMA) and a weighting algorithm that prioritizes wavebands based on their robustness against endmember variability. Spectral mixture analysis provides an estimate of the physically interpretable canopy cover, which is not necessarily correlated with defoliation in mixed-aged plantations due to natural variation in canopy cover as stands age. There is considerable variability in the degree of defoliation as well as in stand age among sites and in this study we found that results were significantly improved by the inclusion of an age correction algorithm for both the multi-spectral (R²_{no age correction} = 0.55 vs R²_{age correction} = 0.73 for Landsat) and hyperspectral (R²_{no age correction} = 0.12 vs R²_{age correction} = 0.50 for Hyperion) image data. The improved accuracy obtained from Landsat compared to the Hyperion data illustrates the potential of applying SMA techniques for analysis of multi-spectral datasets such as MODIS and SPOT-VEGETATION.     

The potential of the greenness and radiation (GR) model to interpret 8-day gross primary production of vegetation

Remote sensing of vegetation gross primary production (GPP) is an important step to analyze terrestrial carbon (C) cycles in response to changing climate. The availability of global networks of C flux measurements provides a valuable opportunity to develop remote sensing based GPP algorithms and test their performances across diverse regions and plant functional types (PFTs). Using 70 global C flux measurements including 24 non-forest (NF), 17 deciduous forest (DF) and 29 evergreen forest (EF), we present the evaluation of an upscaled remote sensing based greenness and radiation (GR) model for GPP estimation. This model is developed using enhanced vegetation index (EVI) and land surface temperature (LST) from the Moderate Resolution Imaging Spectroradiometer (MODIS) and global course resolution radiation data from the National Center for Environmental Prediction (NCEP). Model calibration was achieved using statistical parameters of both EVI and LST fitted for different PFTs. Our results indicate that compared to the standard MODIS GPP product, the calibrated GR model improved the GPP accuracy by reducing the root mean square errors (RMSE) by 16%, 30% and 11% for the NF, DF and EF sites, respectively. The standard MODIS and GR model intercomparisons at individual sites for GPP estimation also showed that GR model performs better in terms of model accuracy and stability. This evaluation demonstrates the potential use of the GR model in capturing short-term GPP variations in areas lacking ground measurements for most of vegetated ecosystems globally.     

Retrieval of chlorophyll and nitrogen in Norway spruce (Picea abies L. Karst.) using imaging spectroscopy

The research evaluated the information content of spectral reflectance (laboratory and airborne data) for the estimation of needle chlorophyll (C_AB) and nitrogen (C_N) concentration in Norway spruce (Picea abies L. Karst.) needles. To identify reliable predictive models different types of spectral transformations were systematically compared regarding the accuracy of prediction. The results of the cross-validated analysis showed that C_AB can be well estimated from laboratory and canopy reflectance data. The best predictive model to estimate C_AB was achieved from laboratory spectra using continuum-removal transformed data (R²cv = 0.83 and a relative RMSEcv of 8.1%, n = 78) and from hyperspectral HyMap data using band-depth normalised spectra (R²cv = 0.90, relative RMSEcv = 2.8%, n = 13). Concerning the nitrogen concentration, we observed somewhat weaker relations, with however still acceptable accuracies (at canopy level: R²cv = 0.57, relative RMSEcv = 4.6%). The wavebands selected in the regression models to estimate C_AB were typically located in the red edge region and near the green reflectance peak. For C_N, additional wavebands related to a known protein absorption feature at 2350 nm were selected. The portion of selected wavebands attributable to known absorption features strongly depends on the type of spectral transformation applied. A method called “water removal” (WR) produced for canopy spectra the largest percentage of wavebands directly or indirectly related to known absorption features. The derived chlorophyll and nitrogen maps may support the detection and the monitoring of environmental stressors and are also important inputs to many bio-geochemical process models.     

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

光学遥感是目前反演植被叶面积指数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。研究表明,采用分段敏感的植被指数,协同光学和微波遥感反演玉米全生长期叶面积指数是可行的。     

Applications of AnnAGNPS model for soil loss estimation and nutrient loading for Malaysian conditions

The study evaluated the performance and suitability of AnnAGNPS model in assessing runoff, sediment loading and nutrient loading under Malaysian conditions. The watershed of River Kuala Tasik in Malaysia, a combination of two sub-watersheds, was selected as the area of study. The data for the year 2004 was used to calibrate the model and the data for the year 2005 was used for validation purposes. Several input parameters were computed using methods suggested by other researchers and studies carried out in Malaysia. The study shows that runoff was predicted well with an overall R² value of 0.90 and E value of 0.70. Sediment loading was able to produce a moderate result of R² = 0.66 and E = 0.49, nitrogen loading predictions were slightly better with R² = 0.68 and E = 0.53, and phosphorus loading performance was slightly poor with an R² = 0.63 and E = 0.33. The erosion map developed was in agreement with the erosion risk map produced by the Department of Agriculture, Malaysia. Rubber estates and urban areas were found to be the main contributors to soil erosion. The simulation results showed that AnnAGNPS has the potential to be used as a valuable tool for planning and management of watersheds under Malaysian conditions.     

Mapping beech (Fagus sylvatica L.) forest structure with airborne hyperspectral imagery

Estimating forest structural attributes using multispectral remote sensing is challenging because of the saturation of multispectral indices at high canopy cover. The objective of this study was to assess the utility of hyperspectral data in estimating and mapping forest structural parameters including mean diameter-at-breast height (DBH), mean tree height and tree density of a closed canopy beech forest (Fagus sylvatica L.). Airborne HyMap images and data on forest structural attributes were collected from the Majella National Park, Italy in July 2004. The predictive performances of vegetation indices (VI) derived from all possible two-band combinations (VI_(i,j) = (R_i − R_j)/(R_i + R_j), where R_i and R_j = reflectance in any two bands) were evaluated using calibration (n = 33) and test (n = 20) data sets. The potential of partial least squares (PLS) regression, a multivariate technique involving several bands was also assessed. New VIs based on the contrast between reflectance in the red-edge shoulder (756–820 nm) and the water absorption feature centred at 1200 nm (1172–1320 nm) were found to show higher correlations with the forest structural parameters than standard VIs derived from NIR and visible reflectance (i.e. the normalised difference vegetation index, NDVI). PLS regression showed a slight improvement in estimating the beech forest structural attributes (prediction errors of 27.6%, 32.6% and 46.4% for mean DBH, height and tree density, respectively) compared to VIs using linear regression models (prediction errors of 27.8%, 35.8% and 48.3% for mean DBH, height and tree density, respectively). Mean DBH was the best predicted variable among the stand parameters (calibration R² = 0.62 for an exponential model fit and standard error of prediction = 5.12 cm, i.e. 25% of the mean). The predicted map of mean DBH revealed high heterogeneity in the beech forest structure in the study area. The spatial variability of mean DBH occurs at less than 450 m. The DBH map could be useful to forest management in many ways, e.g. thinning of coppice to promote diameter growth, to assess the effects of management on forest structure or to detect changes in the forest structure caused by anthropogenic and natural factors.     

Effects of the partitioning of diffuse and direct solar radiation on satellite-based modeling of crop gross primary production

Modeling crop gross primary production (GPP) is critical to understanding the carbon dynamics of agro-ecosystems. Satellite-based studies have widely used production efficiency models (PEM) to estimate cropland GPP, wherein light use efficiency (LUE) is a key model parameter. One factor that has not been well considered in many PEMs is that canopy LUE could vary with illumination conditions. This study investigates how the partitioning of diffuse and direct solar radiation influences cropland GPP using both flux tower and satellite data. The field-measured hourly LUE under cloudy conditions was 1.50 and 1.70 times higher than that under near clear-sky conditions for irrigated corn and soybean, respectively. We applied a two-leaf model to simulate the canopy radiative transfer process, where modeled photosynthetically active radiation (PAR) absorbed by canopy agreed with tower measurements (R² = 0.959 and 0.914 for corn and soybean, respectively). Derived canopy LUE became similar after accounting for the impact of light saturation on leaf photosynthetic capacity under varied illumination conditions. The impacts of solar radiation partitioning on satellite-based modeling of crop GPP was examined using vegetation indices (VI) derived from MODIS data. Consistent with the field modeling results, the relationship between daily GPP and PAR × VI under varied illumination conditions showed different patterns in terms of regression slope and intercept. We proposed a function to correct the influences of direct and diffuse radiation partitioning and the explained variance of flux tower GPP increased in all experiments. Our results suggest that the non-linear response of leaf photosynthesis to light absorption contributes to higher canopy LUE on cloudy days than on clear days. We conclude that accounting for the impacts of solar radiation partitioning is necessary for modeling crop GPP on a daily or shorter basis.     

A GPP assimilation model for the southeastern Tibetan Plateau based on CO2 eddy covariance flux tower and remote sensing data

The gross primary production (GPP) at individual CO₂ eddy covariance flux tower sites (GPP_Tower) in Dali (DL), Wenjiang (WJ) and Linzhi (LZ) around the southeastern Tibetan Plateau were determined by the net ecosystem exchange of CO₂ (NEE) and ecosystem respiration (R_e). The satellite remote sensing-VPM model estimates of GPP values (GPP_MODIS) used the satellite-derived 8-day surface reflectance product (MOD09A1), including satellite-derived enhanced vegetation index (EVI) and land surface water index (LSWI). In this paper, we assembled a subset of flux tower data at these three sites to calibrate and test satellite-VPM model estimated GPP_MODIS, and introduced the satellite data and site-level environmental factors to develop four new assimilation models. The new assimilation models’ estimates of GPP values were compared with GPP_MODIS and GPP_Tower, and the final optimum model among the four assimilation models was determined and used to calibrate GPP_MODIS. The results showed that GPP_MODIS had similar temporal variations to the GPP_Tower, but GPP_MODlS were commonly higher in absolute magnitude than GPP_Tower with relative error (RE) about 58.85%. While, the assimilation models’ estimates of GPP values (GPP_MODEL) were much more closer to GPP_Tower with RE approximately 6.98%, indicating that the capacity of the simulation in the new assimilation model was greatly improved, the R² and root mean square error (RMSE) of the new assimilation model were 0.57–4.90% higher and 0.74–2.47 g C m⁻² s⁻¹ lower than those of the GPP_MODIS, respectively. The assimilation model was used to predicted GPP dynamics around the Tibetan Plateau and showed a reliable result compared with other researches. This study demonstrated the potential of the new assimilation model for estimating GPP around the Tibetan Plateau and the performances of site-level biophysical parameters in related to satellite-VPM model GPP.