Retrieval of land surface temperature from the Kalpana-1 VHRR data using a single-channel algorithm and its validation over western India

Indian geostationary satellite Kalpana-1 (K1) offers a potential to capture the diurnal cycle of land surface temperature (LST) through thermal infrared channel (10.5–12.5 μm) observations of the Very High Resolution Radiometer (VHRR) sensor. A study was carried out to retrieve LST by adapting a generalized single-channel (SC) algorithm (Jiménez-Muñoz and Sobrino, 2003) for the VHRR sensor over India. The basis of SC algorithm depends on the concept of Atmospheric Functions (AFs) that are dependent on transmissivity, upwelling and downwelling radiances of the atmosphere. In the present study AFs were computed for the VHRR sensor through the MODTRAN simulations based upon varying atmospheric and surface inputs. The AFs were fitted with the atmospheric columnar water vapour content and a set of coefficients was derived for LST retrieval. The K1-LST derived with the SC algorithm was validated with (a) in situ measurements at two sites located in western parts of India and (b) the MODIS LST products. Comparison of K1-LST with the in situ measurements demonstrated that SC algorithm was successful in capturing the prominent diurnal variations of 283–332 K in the LST at desert and agriculture experimental sites with a rmse of 1.6 K and 2.7 K, respectively. Inter comparison of K1-LST and MODIS LST showed a reasonable agreement between these two retrievals up to LST of 300 K, however a cold bias up to 7.9 K was observed in MODIS LST for higher LST values (310–330 K) over the hot desert region.     

Hyperspectral-based predictive modelling of grapevine water status in the Portuguese Douro wine region

In this study, hyperspectral reflectance (HySR) data derived from a handheld spectroradiometer were used to assess the water status of three grapevine cultivars in two sub-regions of Douro wine region during two consecutive years. A large set of potential predictors derived from the HySR data were considered for modelling/predicting the predawn leaf water potential (Ψ_pd) through different statistical and machine learning techniques. Three HySR vegetation indices were selected as final predictors for the computation of the models and the in-season time trend was removed from data by using a time predictor. The vegetation indices selected were the Normalized Reflectance Index for the wavelengths 554 nm and 561 nm (NRI_554;561), the water index (WI) for the wavelengths 900 nm and 970 nm, and the D1 index which is associated with the rate of reflectance increase in the wavelengths of 706 nm and 730 nm. These vegetation indices covered the green, red edge and the near infrared domains of the electromagnetic spectrum. A large set of state-of-the-art analysis and statistical and machine-learning modelling techniques were tested. Predictive modelling techniques based on generalized boosted model (GBM), bagged multivariate adaptive regression splines (B-MARS), generalized additive model (GAM), and Bayesian regularized neural networks (BRNN) showed the best performance for predicting Ψ_pd, with an average determination coefficient (R²) ranging between 0.78 and 0.80 and RMSE varying between 0.11 and 0.12 MPa. When cultivar Touriga Nacional was used for training the models and the cultivars Touriga Franca and Tinta Barroca for testing (independent validation), the models performance was good, particularly for GBM (R² = 0.85; RMSE = 0.09 MPa). Additionally, the comparison of Ψ_pd observed and predicted showed an equitable dispersion of data from the various cultivars. The results achieved show a good potential of these predictive models based on vegetation indices to support irrigation scheduling in vineyard.     

Development of a 3-D urbanization index using digital terrain models for surface urban heat island effects

This study assesses surface urban heat island (SUHI) effects during heat waves in subtropical areas. Two cities in northern Taiwan, Taipei metropolis and its adjacent medium-sized city, Yilan, were selected for this empirical study. Daytime and night time surface temperature and SUHI intensity of both cities in five heat wave cases were obtained from MODIS Land-Surface Temperature (LST) and compared. In order to assess SUHI in finer spatial scale, an innovated three-dimensional Urbanization Index (3DUI) with a 5-m spatial resolution was developed to quantify urbanization from a 3-D perspective using Digital Terrain Models (DTMs). The correlation between 3DUI and surface temperatures were also assessed. The results obtained showed that the highest SUHI intensity in daytime was 10.2 °C in Taipei and 7.5 °C in Yilan. The SUHI intensity was also higher than that in non-heat-wave days (about 5 °C) in Taipei. The difference in SUHI intensity of both cities could be as small as only 1.0 °C, suggesting that SUHI intensity was enhanced in both large and medium-sized cities during heat waves. Moreover, the surface temperatures of rural areas in Taipei and Yilan were elevated in the intense heat wave cases, suggesting that the SUHI may reach a plateau when the heat waves get stronger and last longer. In addition, the correlation coefficient between 3DUI and surface temperature was greater than 0.6. The innovative 3DUI can be employed to assess the spatial variation of temperatures and SUHI intensity in much finer spatial resolutions than measurements obtained from remote sensing and weather stations. In summary, the empirical results demonstrated intensified SUHI in large and medium-sized cities in subtropical areas during heat waves which could result in heat stress risks of residents. The innovative 3DUI can be employed to identify vulnerable areas in fine spatial resolutions for formulation of heat wave adaptation strategies.     

Developing a temporally land cover-based look-up table (TL-LUT) method for estimating land surface temperature based on AMSR-E data over the Chinese landmass

The land surface temperature (LST) is an important parameter when studying the interface between the atmosphere and the Earth's surface. Compared to satellite thermal infrared (TIR) remote sensing, passive microwave (PMW) remote sensing is better able to overcome atmospheric influences and to estimate the LST, especially in cloudy regions. However, methods for estimating PMW LSTs at the country and continental scales are still rare. The necessity of training such methods from a temporally dynamic perspective also needs further investigations. Here, a temporally land cover based look-up table (TL-LUT) method is proposed to estimate the LSTs from AMSR-E data over the Chinese landmass. In this method, the synergies between observations from MODIS (Moderate Resolution Imaging Spectroradiometer) and AMSR-E (Advanced Microwave Scanning Radiometer for EOS), which are onboard the same Aqua satellite, are explored. Validation with the synchronous MODIS LSTs demonstrates that the TL-LUT method has better performances in retrieving LSTs with AMSR-E data than the method that uses a single brightness temperature in 36.5 GHz vertical polarization channel. The accuracy of the TL-LUT method is better than 2.7 K for forest and 3.2 K for cropland. Its accuracy varies according to land cover type, time of day, and season. When compared with the in-situ measured LSTs at four sites without urban warming in the Tibet Plateau, the standard errors of estimation between the estimated AMSR-E LST and in-situ measured LST are from 5.1 K to 6.0 K in the daytime and 3.1 K to 4.5 K in the nighttime. Further comparison with the in-situ measured air temperatures at 24 meteorological stations confirms the good performance of the TL-LUT method. The feasibility of PMW remote sensing in estimating the LST for China can complement the TIR data and can, therefore, aid in the generation of daily LST maps for the entire country. Further study of the penetration of PMW radiation would benefit the LST estimations in barren and other sparsely vegetated environments.     

Agriculture pest and disease risk maps considering MSG satellite data and land surface temperature

Pest risk maps for agricultural use are usually constructed from data obtained from in-situ meteorological weather stations, which are relatively sparsely distributed and are often quite expensive to install and difficult to maintain. This leads to the creation of maps with relatively low spatial resolution, which are very much dependent on interpolation methodologies. Considering that agricultural applications typically require a more detailed scale analysis than has traditionally been available, remote sensing technology can offer better monitoring at increasing spatial and temporal resolutions, thereby, improving pest management results and reducing costs. This article uses ground temperature, or land surface temperature (LST), data distributed by EUMETSAT/LSASAF (with a spatial resolution of 3 × 3 km (nadir resolution) and a revisiting time of 15 min) to generate one of the most commonly used parameters in pest modeling and monitoring: “thermal integral over air temperature (accumulated degree-days)”. The results show a clear association between the accumulated LST values over a threshold and the accumulated values computed from meteorological stations over the same threshold (specific to a particular tomato pest). The results are very promising and enable the production of risk maps for agricultural pests with a degree of spatial and temporal detail that is difficult to achieve using in-situ meteorological stations.     

Land surface temperature as potential indicator of burn severity in forest Mediterranean ecosystems

Forest fires are one of the most important causes of environmental alteration in Mediterranean countries. Discrimination of different degrees of burn severity is critical for improving management of fire-affected areas. This paper aims to evaluate the usefulness of land surface temperature (LST) as potential indicator of burn severity. We used a large convention-dominated wildfire, which occurred on 19–21 September, 2012 in Northwestern Spain. From this area, a 1-year series of six LST images were generated from Landsat 7 Enhanced Thematic Mapper (ETM+) data using a single channel algorithm. Further, the Composite Burn Index (CBI) was measured in 111 field plots to identify the burn severity level (low, moderate, and high). Evaluation of the potential relationship between post-fire LST and ground measured CBI was performed by both correlation analysis and regression models. Correlation coefficients were higher in the immediate post-fire LST images, but decreased during the fall of 2012 and increased again with a second maximum value in summer, 2013. A linear regression model between post-fire LST and CBI allowed us to represent spatially predicted CBI (R-squared_adj > 85%). After performing an analysis of variance (ANOVA) between post-fire LST and CBI, a Fisher's least significant difference test determined that two burn severity levels (low-moderate and high) could be statistically distinguished. The identification of such burn severity levels is sufficient and useful to forest managers. We conclude that summer post-fire LST from moderate resolution satellite data may be considered as a valuable indicator of burn severity for large fires in Mediterranean forest ecosytems.     

Defining optimal DEM resolutions and point densities for modelling hydrologically sensitive areas in agricultural catchments dominated by microtopography

Defining critical source areas (CSAs) of diffuse pollution in agricultural catchments depends upon the accurate delineation of hydrologically sensitive areas (HSAs) at highest risk of generating surface runoff pathways. In topographically complex landscapes, this delineation is constrained by digital elevation model (DEM) resolution and the influence of microtopographic features. To address this, optimal DEM resolutions and point densities for spatially modelling HSAs were investigated, for onward use in delineating CSAs. The surface runoff framework was modelled using the Topographic Wetness Index (TWI) and maps were derived from 0.25 m LiDAR DEMs (40 bare-earth points m⁻²), resampled 1 m and 2 m LiDAR DEMs, and a radar generated 5 m DEM. Furthermore, the resampled 1 m and 2 m LiDAR DEMs were regenerated with reduced bare-earth point densities (5, 2, 1, 0.5, 0.25 and 0.125 points m⁻²) to analyse effects on elevation accuracy and important microtopographic features. Results were compared to surface runoff field observations in two 10 km² agricultural catchments for evaluation. Analysis showed that the accuracy of modelled HSAs using different thresholds (5%, 10% and 15% of the catchment area with the highest TWI values) was much higher using LiDAR data compared to the 5 m DEM (70–100% and 10–84%, respectively). This was attributed to the DEM capturing microtopographic features such as hedgerow banks, roads, tramlines and open agricultural drains, which acted as topographic barriers or channels that diverted runoff away from the hillslope scale flow direction. Furthermore, the identification of ‘breakthrough’ and ‘delivery’ points along runoff pathways where runoff and mobilised pollutants could be potentially transported between fields or delivered to the drainage channel network was much higher using LiDAR data compared to the 5 m DEM (75–100% and 0–100%, respectively). Optimal DEM resolutions of 1–2 m were identified for modelling HSAs, which balanced the need for microtopographic detail as well as surface generalisations required to model the natural hillslope scale movement of flow. Little loss of vertical accuracy was observed in 1–2 m LiDAR DEMs with reduced bare-earth point densities of 2–5 points m⁻², even at hedgerows. Further improvements in HSA models could be achieved if soil hydrological properties and the effects of flow sinks (filtered out in TWI models) on hydrological connectivity are also considered.     

A spectral index for land degradation mapping using ASTER data: Application to a semi-arid Mediterranean catchment

Flagrant soil erosion in Morocco is an alarming sign of soil degradation. Due to the considerable costs of detailed ground surveys of this phenomenon, remote sensing is an appropriate alternative for analyzing and evaluating the risks of the expansion of soil degradation. In this paper, we characterize the state of land degradation in a small Mediterranean watershed using Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data and ground-based spectroradiometric measurements. The two visible, the near-infrared and six shortwave infrared bands of the above sensor were calibrated using ground measurements of the spectral reflectance. Field measurements were carried out in the Saboun experimental basin located in the marl soil region of the Moroccan western Rif. The study leads to the development and evaluation of a new spectral approach to express land degradation. This index called Land degradation index (LDI) is based on the concept of the soil line derived from spectroradiometric ground measurements. In this study, we compare LDI and the spectral angle mapping (SAM) approaches to assess and map land degradation. Results show that LDI provides more accurate results for mapping land degradation (Kappa = 0.79) when compared to the SAM method (Kappa = 0.61). Validation and evaluation of the results are based on the thematic maps derived from the ground data (organic matter, clay, silt and sand) by kriging, DEM, slope gradient and photointerpretation.     

Retrieval of leaf water content spanning the visible to thermal infrared spectra

The objective of this study was to investigate the entire spectra (from visible to the thermal infrared; 0.390–14.0 μm) to retrieve leaf water content in a consistent manner. Narrow-band spectral indices (calculated from all possible two band combinations) and a partial least square regression (PLSR) were used to assess the strength of each spectral region. The coefficient of determination (R²) and root mean square error (RMSE) were used to report the prediction accuracy of spectral indices and PLSR models. In the visible-near infrared and shortwave infrared (VNIR–SWIR), the most accurate spectral index yielded R² of 0.89 and RMSE of 7.60%, whereas in the mid infrared (MIR) the highest R² was 0.93 and RMSE of 5.97%. Leaf water content was poorly predicted using two-band indices developed from the thermal infrared (R² = 0.33). The most accurate PLSR model resulted from MIR reflectance spectra (R² = 0.96, RMSE = 4.74% and RMSE cross validation RMSE_CV = 6.17%) followed by VNIR–SWIR reflectance spectra (R² = 0.91, RMSE = 6.90% and RMSE_CV = 7.32%). Using thermal infrared (TIR) spectra, the PLSR model yielded a moderate retrieval accuracy (R² = 0.67, RMSE = 13.27% and RMSE_CV = 16.39%). This study demonstrated that the mid infrared (MIR) and shortwave infrared (SWIR) domains were the most sensitive spectral region for the retrieval of leaf water content.     

A visible band index for remote sensing leaf chlorophyll content at the canopy scale

Leaf chlorophyll content is an important variable for agricultural remote sensing because of its close relationship to leaf nitrogen content. The triangular greenness index (TGI) was developed based on the area of a triangle surrounding the spectral features of chlorophyll with points at (670 nm, R₆₇₀), (550 nm, R₅₅₀), and (480 nm, R₄₈₀), where R_λ is the spectral reflectance at wavelengths of 670, 550 and 480, respectively. The equation is TGI = −0.5[(670 − 480)(R₆₇₀ − R₅₅₀) − (670 − 550)(R₆₇₀ − R₄₈₀)]. In 1999, investigators funded by NASA's Earth Observations Commercialization and Applications Program collaborated on a nitrogen fertilization experiment with irrigated maize in Nebraska. Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data and Landsat 5 Thematic Mapper (TM) data were acquired along with leaf chlorophyll meter and other data on three dates in July during late vegetative growth and early reproductive growth. TGI was consistently correlated with plot-averaged chlorophyll-meter values at the spectral resolutions of AVIRIS, Landsat TM, and digital cameras. Simulations using the Scattering by Arbitrarily Inclined Leaves (SAIL) canopy model indicate an interaction among TGI, leaf area index (LAI) and soil type at low crop LAI, whereas at high LAI and canopy closure, TGI was only affected by leaf chlorophyll content. Therefore, TGI may be the best spectral index to detect crop nitrogen requirements with low-cost digital cameras mounted on low-altitude airborne platforms.     

Modelling high arctic percent vegetation cover using field digital images and high resolution satellite data

In this study, digital images collected at a study site in the Canadian High Arctic were processed and classified to examine the spatial-temporal patterns of percent vegetation cover (PVC). To obtain the PVC of different plant functional groups (i.e., forbs, graminoids/sedges and mosses), field near infrared-green-blue (NGB) digital images were classified using an object-based image analysis (OBIA) approach. The PVC analyses comparing different vegetation types confirmed: (i) the polar semi-desert exhibited the lowest PVC with a large proportion of bare soil/rock cover; (ii) the mesic tundra cover consisted of approximately 60% mosses; and (iii) the wet sedge consisted almost exclusively of graminoids and sedges. As expected, the PVC and green normalized difference vegetation index (GNDVI; (R_NIR − R_Green)/(R_NIR + R_Green)), derived from field NGB digital images, increased during the summer growing season for each vegetation type: i.e., ∼5% (0.01) for polar semi-desert; ∼10% (0.04) for mesic tundra; and ∼12% (0.03) for wet sedge respectively. PVC derived from field images was found to be strongly correlated with WorldView-2 derived normalized difference spectral indices (NDSI; (R_x − R_y)/(R_x + R_y)), where R_x is the reflectance of the red edge (724.1 nm) or near infrared (832.9 nm and 949.3 nm) bands; R_y is the reflectance of the yellow (607.7 nm) or red (658.8 nm) bands with R²’s ranging from 0.74 to 0.81. NDSIs that incorporated the yellow band (607.7 nm) performed slightly better than the NDSIs without, indicating that this band may be more useful for investigating Arctic vegetation that often includes large proportions of senescent vegetation throughout the growing season.     

Species’ habitat use inferred from environmental variables at multiple scales: How much we gain from high-resolution vegetation data?

Spatial resolution of environmental data may influence the results of habitat selection models. As high-resolution data are usually expensive, an assessment of their contribution to the reliability of habitat models is of interest for both researchers and managers. We evaluated how vegetation cover datasets of different spatial resolutions influence the inferences and predictive power of multi-scale habitat selection models for the endangered brown bear populations in the Cantabrian Range (NW Spain). We quantified the relative performance of three types of datasets: (i) coarse resolution data from Corine Land Cover (minimum mapping unit of 25 ha), (ii) medium resolution data from the Forest Map of Spain (minimum mapping unit of 2.25 ha and information on forest canopy cover and tree species present in each polygon), and (iii) high-resolution Lidar data (about 0.5 points/m²) providing a much finer information on forest canopy cover and height. Despite all the models performed well (AUC > 0.80), the predictive ability of multi-scale models significantly increased with spatial resolution, particularly when other predictors of habitat suitability (e.g. human pressure) were not used to indirectly filter out areas with a more degraded vegetation cover. The addition of fine grain information on forest structure (LiDAR) led to a better understanding of landscape use and a more accurate spatial representation of habitat suitability, even for a species with large spatial requirements as the brown bear, which will result in the development of more effective measures to assist endangered species conservation.     

Detecting Sirex noctilio grey-attacked and lightning-struck pine trees using airborne hyperspectral data,random forest and support vector machines classifiers

The visual progression of sirex (Sirex noctilio) infestation symptoms has been categorized into three distinct infestation phases, namely the green, red and grey stages. The grey stage is the final stage which leads to almost complete defoliation resulting in dead standing trees or snags. Dead standing pine trees however, could also be due to the lightning damage. Hence, the objective of the present study was to distinguish amongst healthy, sirex grey-attacked and lightning-damaged pine trees using AISA Eagle hyperspectral data, random forest (RF) and support vector machines (SVM) classifiers. Our study also presents an opportunity to look at the possibility of separating amongst the previously mentioned pine trees damage classes and other landscape classes on the study area. The results of the present study revealed the robustness of the two machine learning classifiers with an overall accuracy of 74.50% (total disagreement = 26%) for RF and 73.50% (total disagreement = 27%) for SVM using all the remaining AISA Eagle spectral bands after removing the noisy ones. When the most useful spectral bands as measured by RF were exploited, the overall accuracy was considerably improved; 78% (total disagreement = 22%) for RF and 76.50% (total disagreement = 24%) for SVM. There was no significant difference between the performances of the two classifiers as demonstrated by the results of McNemar’s test (chi-squared; χ² = 0.14, and 0.03 when all the remaining ASIA Eagle wavebands, after removing the noisy ones and the most important wavebands were used, respectively). This study concludes that AISA Eagle data classified using RF and SVM algorithms provide relatively accurate information that is important to the forest industry for making informed decision regarding pine plantations health protocols.     

Modeling soil parameters using hyperspectral image reflectance in subtropical coastal wetlands

Developing spectral models of soil properties is an important frontier in remote sensing and soil science. Several studies have focused on modeling soil properties such as total pools of soil organic matter and carbon in bare soils. We extended this effort to model soil parameters in areas densely covered with coastal vegetation. Moreover, we investigated soil properties indicative of soil functions such as nutrient and organic matter turnover and storage. These properties include the partitioning of mineral and organic soil between particulate (>53 μm) and fine size classes, and the partitioning of soil carbon and nitrogen pools between stable and labile fractions. Soil samples were obtained from Avicennia germinans mangrove forest and Juncus roemerianus salt marsh plots on the west coast of central Florida. Spectra corresponding to field plot locations from Hyperion hyperspectral image were extracted and analyzed. The spectral information was regressed against the soil variables to determine the best single bands and optimal band combinations for the simple ratio (SR) and normalized difference index (NDI) indices. The regression analysis yielded levels of correlation for soil variables with R² values ranging from 0.21 to 0.47 for best individual bands, 0.28 to 0.81 for two-band indices, and 0.53 to 0.96 for partial least-squares (PLS) regressions for the Hyperion image data. Spectral models using Hyperion data adequately (RPD > 1.4) predicted particulate organic matter (POM), silt + clay, labile carbon (C), and labile nitrogen (N) (where RPD = ratio of standard deviation to root mean square error of cross-validation [RMSECV]). The SR (0.53 μm, 2.11 μm) model of labile N with R² = 0.81, RMSECV= 0.28, and RPD = 1.94 produced the best results in this study. Our results provide optimism that remote-sensing spectral models can successfully predict soil properties indicative of ecosystem nutrient and organic matter turnover and storage, and do so in areas with dense canopy cover.     

Monitoring land surface albedo and vegetation dynamics using high spatial and temporal resolution synthetic time series from Landsat and the MODIS BRDF/NBAR/albedo product

Seasonal vegetation phenology can significantly alter surface albedo which in turn affects the global energy balance and the albedo warming/cooling feedbacks that impact climate change. To monitor and quantify the surface dynamics of heterogeneous landscapes, high temporal and spatial resolution synthetic time series of albedo and the enhanced vegetation index (EVI) were generated from the 500 m Moderate Resolution Imaging Spectroradiometer (MODIS) operational Collection V006 daily BRDF/NBAR/albedo products and 30 m Landsat 5 albedo and near-nadir reflectance data through the use of the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). The traditional Landsat Albedo (Shuai et al., 2011) makes use of the MODIS BRDF/Albedo products (MCD43) by assigning appropriate BRDFs from coincident MODIS products to each Landsat image to generate a 30 m Landsat albedo product for that acquisition date. The available cloud free Landsat 5 albedos (due to clouds, generated every 16 days at best) were used in conjunction with the daily MODIS albedos to determine the appropriate 30 m albedos for the intervening daily time steps in this study. These enhanced daily 30 m spatial resolution synthetic time series were then used to track albedo and vegetation phenology dynamics over three Ameriflux tower sites (Harvard Forest in 2007, Santa Rita in 2011 and Walker Branch in 2005). These Ameriflux sites were chosen as they are all quite nearby new towers coming on line for the National Ecological Observatory Network (NEON), and thus represent locations which will be served by spatially paired albedo measures in the near future. The availability of data from the NEON towers will greatly expand the sources of tower albedometer data available for evaluation of satellite products. At these three Ameriflux tower sites the synthetic time series of broadband shortwave albedos were evaluated using the tower albedo measurements with a Root Mean Square Error (RMSE) less than 0.013 and a bias within the range of ±0.006. These synthetic time series provide much greater spatial detail than the 500 m gridded MODIS data, especially over more heterogeneous surfaces, which improves the efforts to characterize and monitor the spatial variation across species and communities. The mean of the difference between maximum and minimum synthetic time series of albedo within the MODIS pixels over a subset of satellite data of Harvard Forest (16 km by 14 km) was as high as 0.2 during the snow-covered period and reduced to around 0.1 during the snow-free period. Similarly, we have used STARFM to also couple MODIS Nadir BRDF Adjusted Reflectances (NBAR) values with Landsat 5 reflectances to generate daily synthetic times series of NBAR and thus Enhanced Vegetation Index (NBAR-EVI) at a 30 m resolution. While normally STARFM is used with directional reflectances, the use of the view angle corrected daily MODIS NBAR values will provide more consistent time series. These synthetic times series of EVI are shown to capture seasonal vegetation dynamics with finer spatial and temporal details, especially over heterogeneous land surfaces.     

Estimating winter wheat biomass by assimilating leaf area index derived from fusion of Landsat-8 and MODIS data

A sufficient number of satellite acquisitions in a growing season are essential for deriving agronomic indicators, such as green leaf area index (GLAI), to be assimilated into crop models for crop productivity estimation. However, for most high resolution orbital optical satellites, it is often difficult to obtain images frequently due to their long revisit cycles and unfavorable weather conditions. Data fusion algorithms, such as the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) and the Enhanced STARFM (ESTARFM), have been developed to generate synthetic data with high spatial and temporal resolution to address this issue. In this study, we evaluated the approach of assimilating GLAI into the Simple Algorithm for Yield Estimation model (SAFY) for winter wheat biomass estimation. GLAI was estimated using the two-band Enhanced Vegetation Index (EVI2) derived from data acquired by the Operational Land Imager (OLI) onboard the Landsat-8 and a fusion dataset generated by blending the Moderate-Resolution Imaging Spectroradiometer (MODIS) data and the OLI data using the STARFM and ESTARFM models. The fusion dataset had the temporal resolution of the MODIS data and the spatial resolution of the OLI data. Key parameters of the SAFY model were optimised through assimilation of the estimated GLAI into the crop model using the Shuffled Complex Evolution-University of Arizona (SCE-UA) algorithm. A good agreement was achieved between the estimated and field measured biomass by assimilating the GLAI derived from the OLI data (GLAI_L) alone (R² = 0.77 and RMSE = 231 g m⁻²). Assimilation of GLAI derived from the fusion dataset (GLAI_F) resulted in a R² of 0.71 and RMSE of 193 g m⁻² while assimilating the combination of GLAI_L and GLAI_F led to further improvements (R² = 0.76 and RMSE = 176 g m⁻²). Our results demonstrated the potential of using the fusion algorithms to improve crop growth monitoring and crop productivity estimation when the number of high resolution remote sensing data acquisitions is limited.     

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.     

光谱指数趋势面的城市地表温度降尺度转换

针对城市及周边区域建造区和自然地表交织分布的特点,探讨了利用归一化植被指数(NDVI)和归一化建造指数(NDBI)构造趋势面的地表温度(LST)降尺度方法,以北京市市区及周边较平坦区域为例实现了LST自960 m向120 m的降尺度转换。分析了LST空间分布特征及NDVI、NDBI对地物的指示性特征;以北京市四至六环为界分析NDVI、NDBI趋势面对地表温度的拟合程度及各自的适用区域;在120 m、240 m、480 m和960 m 4个尺度上评价了NDVI、NDBI和NDVI+NDBI趋势面对LST的拟合程度和趋势面转换函数的尺度效应;对NDVI、NDBI和NDVI NDBI等3种方法的降尺度结果分覆盖类型、分区域对比评价。实验结果表明结合两种光谱指数的NDVI NDBI方法降尺度转换精度有所改善,改善程度取决于地表覆盖类型组合。     

Remote quantification of phycocyanin in potable water sources through an adaptive model

Cyanobacterial blooms in water supply sources in both central Indiana USA (CIN) and South Australia (SA) are a cause of great concerns for toxin production and water quality deterioration. Remote sensing provides an effective approach for quick assessment of cyanobacteria through quantification of phycocyanin (PC) concentration. In total, 363 samples spanning a large variation of optically active constituents (OACs) in CIN and SA waters were collected during 24 field surveys. Concurrently, remote sensing reflectance spectra (R_rs) were measured. A partial least squares–artificial neural network (PLS–ANN) model, artificial neural network (ANN) and three-band model (TBM) were developed or tuned by relating the R_rs with PC concentration. Our results indicate that the PLS–ANN model outperformed the ANN and TBM with both the original spectra and simulated ESA/Sentinel-3/Ocean and Land Color Instrument (OLCI) and EO-1/Hyperion spectra. The PLS–ANN model resulted in a high coefficient of determination (R²) for CIN dataset (R² = 0.92, R: 0.3–220.7 μg/L) and SA (R² = 0.98, R: 0.2–13.2 μg/L). In comparison, the TBM model yielded an R² = 0.77 and 0.94 for the CIN and SA datasets, respectively; while the ANN obtained an intermediate modeling accuracy (CIN: R² = 0.86; SA: R² = 0.95). Applying the simulated OLCI and Hyperion aggregated datasets, the PLS–ANN model still achieved good performance (OLCI: R² = 0.84; Hyperion: R² = 0.90); the TBM also presented acceptable performance for PC estimations (OLCI: R² = 0.65, Hyperion: R² = 0.70). Based on the results, the PLS–ANN is an effective modeling approach for the quantification of PC in productive water supplies based on its effectiveness in solving the non-linearity of PC with other OACs. Furthermore, our investigation indicates that the ratio of inorganic suspended matter (ISM) to PC concentration has close relationship to modeling relative errors (CIN: R² = 0.81; SA: R² = 0.92), indicating that ISM concentration exert significant impact on PC estimation accuracy.     

Application of 3D triangulations of airborne laser scanning data to estimate boreal forest leaf area index

We propose 3D triangulations of airborne Laser Scanning (ALS) point clouds as a new approach to derive 3D canopy structures and to estimate forest canopy effective LAI (LAI_e). Computational geometry and topological connectivity were employed to filter the triangulations to yield a quasi-optimal relationship with the field measured LAI_e. The optimal filtering parameters were predicted based on ALS height metrics, emulating the production of maps of LAI_e and canopy volume for large areas. The LAI_e from triangulations was validated with field measured LAI_e and compared with a reference LAI_e calculated from ALS data using logarithmic model based on Beer’s law. Canopy transmittance was estimated using All Echo Cover Index (ACI), and the mean projection of unit foliage area (β) was obtained using no-intercept regression with field measured LAI_e. We investigated the influence species and season on the triangulated LAI_e and demonstrated the relationship between triangulated LAI_e and canopy volume. Our data is from 115 forest plots located at the southern boreal forest area in Finland and for each plot three different ALS datasets were available to apply the triangulations. The triangulation approach was found applicable for both leaf-on and leaf-off datasets after initial calibration. Results showed the Root Mean Square Errors (RMSEs) between LAI_e from triangulations and field measured values agreed the most using the highest pulse density data (RMSE = 0.63, the coefficient of determination (R²) = 0.53). Yet, the LAI_e calculated using ACI-index agreed better with the field measured LAI_e (RMSE = 0.53 and R² = 0.70). The best models to predict the optimal alpha value contained the ACI-index, which indicates that within-crown transmittance is accounted by the triangulation approach. The cover indices may be recommended for retrieving LAI_e only, but for applications which require more sophisticated information on canopy shape and volume, such as radiative transfer models, the triangulation approach may be preferred.