Impact of satellite based PAR on estimates of terrestrial net primary productivity

Biospheric productivity regulates the supply of food for mankind and therefore, there is a need to estimate its magnitude. The productivity is controlled by the process of photosynthesis driven by solar radiation, primarily in the visible part of the spectrum (0.4–0.7 μm), known as the photosynthetically active radiation (PAR). Current methods to estimate terrestrial net primary production (NPP) use remotely sensed information on vegetation dynamics. Satellite based estimates of PAR are available at a global scale but have seldom been used for estimating NPP. In this study we show that the use of PAR information from satellites does have an impact on estimates of NPP and that there are detectable differences when compared to similar estimates based on conventional PAR information. Net primary production tends to be higher when compared to estimates based on total shortwave (SW) radiation with PAR assumed to be a constant fraction of SW. We focus on the United States during 2004. Net primary production is generally underestimated in regions with mesic environment while overestimated in areas with xeric environment. The most pronounced underestimated region is the southeast United States. The study demonstrates the usefulness of the satellite-based estimates of PAR for modelling terrestrial primary productivity.     

Adapting a global stratified random sample for regional estimation of forest cover change derived from satellite imagery

A desirable feature of a global sampling design for estimating forest cover change based on satellite imagery is the ability to adapt the design to obtain precise regional estimates, where a region may be a country, state, province, or conservation area. A sampling design stratified by an auxiliary variable correlated with forest cover change has this adaptability. A global stratified random sample can be augmented by additional sample units within a region selected by the same stratified protocol and the resulting sample constitutes a stratified random sample of the region. Stratified sampling allows increasing the sample size in a region by a few to many additional sample units. The additional sample units can be effectively allocated to strata to reduce the standard errors of the regional estimates, even though these strata were not initially constructed for the objective of regional estimation. A complete coverage map of deforestation within the Brazilian Legal Amazon (BLA) is used as a population to evaluate precision of regional estimates obtained by augmenting a global stratified random sample. The standard errors of the regional estimates for the BLA and states within the BLA obtained from the augmented stratified design were generally smaller than those attained by simple random sampling and systematic sampling.     

Surface motion of mountain glaciers derived from satellite optical imagery

A complete and detailed map of the ice-velocity field on mountain glaciers is obtained by cross-correlating SPOT5 optical images. This approach offers an alternative to SAR interferometry, because no present or planned RADAR satellite mission provides data with a temporal separation short enough to derive the displacements of glaciers. The methodology presented in this study does not require ground control points (GCPs). The key step is a precise relative orientation of the two images obtained by adjusting the stereo model of one “slave”' image assuming that the other “master” image is well georeferenced. It is performed with numerous precisely-located homologous points extracted automatically. The strong ablation occurring during summer time on the glaciers requires a correction to obtain unbiased displacements. The accuracy of our measurement is assessed based on a comparison with nearly simultaneous differential GPS surveys performed on two glaciers of the Mont Blanc area (Alps). If the images have similar incidence angles and correlate well, the accuracy is on the order of 0.5 m, or 1/5 of the pixel size. Similar results are also obtained without GCPs. An acceleration event, observed in early August for the Mer de Glace glacier, is interpreted in term of an increase in basal sliding. Our methodology, applied to SPOT5 images, can potentially be used to derive the displacements of the Earth's surface caused by landslides, earthquakes, and volcanoes.     

Classification of Amazonian primary rain forest vegetation using Landsat ETM+ satellite imagery

A lack of spatially and thematically accurate vegetation maps complicates conservation and management planning, as well as ecological research, in tropical rain forests. Remote sensing has considerable potential to provide such maps, but classification accuracy within primary rain forests has generally been inadequate for practical applications. Here we test how accurately floristically defined forest types in lowland tropical rain forests in Peruvian Amazonia can be recognized using remote sensing data (Landsat ETM+ satellite image and STRM elevation model). Floristic data and a vegetation classification with four forest classes were available for eight line transects, each 8 km long, located in an area of ca 800 km². We compared two sampling unit sizes (line transect subunits of 200 and 500 m) and several image feature combinations to analyze their suitability for image classification. Mantel tests were used to quantify how well the patterns in elevation and in the digital numbers of the satellite image correlated with the floristic patterns observed in the field. Most Mantel correlations were positive and highly significant. Linear discriminant analysis was used first to build a function that discriminates between forest classes in the eight field-verified transects on the basis of remotely sensed data, and then to classify those parts of the line transects and the satellite image that had not been visited in the field. Classification accuracy was quantified by 8-fold crossvalidation. Two of the tierra firme (non-inundated) forest types were combined because they were too often misclassified. The remaining three forest types (inundated forest, terrace forest and Pebas formation/intermediate tierra firme forest) could be separated using the 500-m sampling units with an overall classification accuracy of 85% and a Kappa coefficient of 0.62. For the 200-m sampling units, the classification accuracy was clearly lower (71%, Kappa 0.35). The forest classification will be used as habitat data to study wildlife habitat use in the same area. Our results show that remotely sensed data and relatively simple classification methods can be used to produce reasonably accurate forest type classifications, even in structurally homogeneous primary rain forests.     

Estimating census district populations from satellite imagery: Some approaches and limitations

Small-area population densities and counts were estimated for Australian census collection districts (CDs), using Landsat TM imagery. A number of mathematical and statistical refinements to previously reported methods were explored. The robustness of these techniques as a practical methodology for population estimation was investigated and evaluated using a primary image for model development and training, and a second image for validation. Correlations of up to 0.92 in the training set and up to 0.86 in the validation set were obtained between census and remote sensing estimates of CD population density, with median proportional errors of 17.4% and 18.4%, respectively. Total urban populations were estimated with errors of +1% and-3%, respectively. These results indicate a moderate level of accuracy and a substantial degree of robustness. Accuracy was greatest in suburban areas of intermediate population density. There was a general tendency towards attenuation in all models tested, with high densities being under-estimated and low densities being over-estimated. It is concluded that the level of accuracy obtainable with this methodology is limited by heterogeneity within the individual CDs, particularly large rural CDs, and that further improvements are in principle unlikely using the aggregated approach. An alternative statistical approach is foreshadowed.     

Estimating forest stand structure from high resolution imagery using the directional variogram

The height and stocking of forest stands can be estimated with relatively high precision using an empirical model relating parameters extracted from the directional variogram of high resolution images and forest structure parameters. A geometrical-optical model of the forest was first used to generate images of artificial forest stands in order to establish the relation between tree size. tree density and image texture. The resulting equations were then applied on the computer generated images as well as on high resolution MEIS II images to predict the forest structure values. The results show a good concordance between actual and predicted values, even when spatial resolution was degraded from 0·36m to 2·16m.     

Estimating areal means and variances of forest attributes using the k-Nearest Neighbors technique and satellite imagery

The k-Nearest Neighbor (k-NN) technique has become extremely popular for a variety of forest inventory mapping and estimation applications. Much of this popularity may be attributed to the non-parametric, multivariate features of the technique, its intuitiveness, and its ease of use. When used with satellite imagery and forest inventory plot data, the technique has been shown to produce useful estimates of many forest attributes including forest/non-forest, volume, and basal area. However, variance estimators for quantifying the uncertainty of means or sums of k-NN pixel-level predictions for areas of interest (AOI) consisting of multiple pixels have not been reported. The primary objectives of the study were to derive variance estimators for AOI estimates obtained from k-NN predictions and to compare precision estimates resulting from different approaches to k-NN prediction and different interpretations of those predictions. The approaches were illustrated by estimating proportion forest area, tree volume per unit area, tree basal area per unit area, and tree density per unit area for 10-km AOIs. Estimates obtained using k-NN approaches and traditional inventory approaches were compared and found to be similar. Further, variance estimates based on different interpretations of k-NN predictions were similar. The results facilitate small area estimation and simultaneous and consistent mapping and estimation of multiple forest attributes.     

Estimating aboveground carbon in a catchment of the Siberian forest tundra: Combining satellite imagery and field inventory

This study was part of an interdisciplinary research project on soil carbon and phytomass dynamics of boreal and arctic permafrost landscapes. The 45 ha study area was a catchment located in the forest tundra in northern Siberia, approximately 100 km north of the Arctic Circle.The objective of this study was to estimate aboveground carbon (AGC) and assess and model its spatial variability. We combined multi-spectral high resolution remote sensing imagery and sample based field inventory data by means of the k-nearest neighbor (k-NN) technique and linear regression.Field data was collected by stratified systematic sampling in August 2006 with a total sample size of n = 31 circular nested sample plots of 154 m² for trees and shrubs and 1 m² for ground vegetation. Destructive biomass samples were taken on a sub-sample for fresh weight and moisture content. Species-specific allometric biomass models were constructed to predict dry biomass from diameter at breast height (dbh) for trees and from elliptic projection areas for shrubs.Quickbird data (standard imagery product), acquired shortly before the field campaign and archived ASTER data (Level-1B product) of 2001 were geo-referenced, converted to calibrated radiances at sensor and used as carrier data. Spectral information of the pixels which were located in the inventory plots were extracted and analyzed as reference set. Stepwise multiple linear regression was applied to identify suitable predictors from the set of variables of the original satellite bands, vegetation indices and texture metrics. To produce thematic carbon maps, carbon values were predicted for all pixels of the investigated satellite scenes. For this prediction, we compared the kNN distance-weighted classifier and multiple linear regression with respect to their predictions.The estimated mean value of aboveground carbon from stratified sampling in the field is 15.3 t/ha (standard error SE = 1.50 t/ha, SE% = 9.8%). Zonal prediction from the k-NN method for the Quickbird image as carrier is 14.7 t/ha with a root mean square error RMSE = 6.42 t/ha, RMSE_r = 44%) resulting from leave-one-out cross-validation. The k-NN-approach allows mapping and analysis of the spatial variability of AGC. The results show high spatial variability with AGC predictions ranging from 4.3 t/ha to 28.8 t/ha, reflecting the highly heterogeneous conditions in those permafrost-influenced landscapes. The means and totals of linear regression and k-NN predictions revealed only small differences but some regional distinctions were recognized in the maps.     

Estimating forest variables from top-of-atmosphere radiance satellite measurements using coupled radiative transfer models

Traditionally, it is necessary to pre-process remote sensing data to obtain top of canopy (TOC) reflectances before applying physically-based model inversion techniques to estimate forest variables. Corrections for atmospheric, adjacency, topography, and surface directional effects are applied sequentially and independently, accumulating errors into the TOC reflectance data, which are then further used in the inversion process. This paper presents a proof of concept for demonstrating the direct use of measured top-of-atmosphere (TOA) radiance data to estimate forest biophysical and biochemical variables, by using a coupled canopy-atmosphere radiative transfer model. Advantages of this approach are that no atmospheric correction is needed and that atmospheric, adjacency, topography, and surface directional effects can be directly and more accurately included in the forward modelling.In the case study, we applied both TOC and TOA approaches to three Norway spruce stands in Eastern Czech Republic. We used the SLC soil-leaf-canopy model and the MODTRAN4 atmosphere model. For the TOA approach, the physical coupling between canopy and atmosphere was performed using a generic method based on the 4-stream radiative transfer theory which enables full use of the directional reflectance components provided by SLC. The method uses three runs of the atmosphere model for Lambertian surfaces, and thus avoids running the atmosphere model for each new simulation. We used local sensitivity analysis and singular value decomposition to determine which variables could be estimated, namely: canopy cover, fraction of bark, needle chlorophyll, and dry matter content. TOC and TOA approaches resulted in different sets of estimates, but had comparable performance. The TOC approach, however, was at its best potential because of the flatness and homogeneity of the area. On the contrary, the capacities of the TOA approach would be better exploited in heterogeneous rugged areas. We conclude that, having similar performance, the TOA approach should be preferred in situations where minimizing the pre-processing is important, such as in data assimilation and multi-sensor studies.     

Estimating above-ground net primary productivity of the tallgrass prairie ecosystem of the Central Great Plains using AVHRR NDVI

Above-ground net primary productivity (ANPP) is indicative of an ecosystem's ability to capture solar energy and convert it to organic carbon (or biomass), which may be used by consumers or decomposers, or stored in the form of living and nonliving organic matter. Annual and interannual variation in ANPP is often linked to climate dynamics and anthropogenic influences, such as fertilization, irrigation, above-ground biomass harvest, and so on. The Central Great Plains grasslands occupy over 1.5 million km² and are a primary resource for livestock production in North America. The tallgrass prairies are the most productive grasslands in this region, and the Flint Hills of North America represent the largest contiguous area of unploughed tallgrass prairie (1.6 million ha). Measurements of ANPP are of critical importance to the proper management and understanding of climatic and anthropogenic influences on tallgrass prairie. Yet, accurate, detailed, and systematic measurements of ANPP over large geographic regions do not exist for this ecosystem. For these reasons, this study was conducted to investigate the use of the normalized difference vegetation index (NDVI) to model ANPP of the tallgrass prairie. Many studies have established a positive relationship between the NDVI and ANPP, but the strength of this relationship is influenced by vegetation types and can vary significantly from year to year depending on land use and climatic conditions. The goal of this study was to develop a robust model using the Advanced Very High Resolution Radiometer (AVHRR) biweekly NDVI values to predict tallgrass ANPP. This study was conducted using ANPP measurements from a watershed within the Konza Prairie Biological Station (KPBS) as the primary study area, with additional measurements from the Rannells Flint Hills Prairie Preserve (RFHPP) and biennial ANPP measurements by Kansas State University (KSU) students from tallgrass areas near Manhattan, Kansas. Data from the primary study site covered the period of 1989–2005. The optimal period for estimating ANPP using AVHRR NDVI composite data sets was found to be late July. The Tallgrass ANPP Model (TAM) explained 54% (coefficient of determination, R ² = 0.54, p < 0.001) of the year-to-year variation in ANPP. The creation of 1.0 km × 1.0 km resolution ANPP maps for a four-county (～7000 ha) area for years 1989–2007 showed considerable variation in annual and interannual ANPP spatial patterns, suggesting complex interactions among factors influencing ANPP spatially and temporally.     

Technical note. Illustration of the influence of shadowing on high latitude information derived from satellite imagery

Two separate types of contamination by shadowing have been identified following an analysis of two Thematic Mapper scenes from the Arctic. The well-established effect of orographic shadowing is particularly important for cryospheric surfaces. Cirrus clouds, often very difficult to identify by automated techniques, also cast shadows which decrease the radiances detected from snow and ice surfaces. These effects are illustrated here in relation to snow mapping algorithms.     

Topographic effects on estimating net primary productivity of green coniferous forest in complex terrain using Landsat data: a case study of Yoshino Mountain,Japan

In mountainous areas, irregular terrain significantly affects spatial variations of climatic variables and the reflectance of pixels in remote sensing imagery. Consequently, the variations may affect the estimation of net primary productivity (NPP). The light-use efficiency (LUE) model is used to analyse topographic influence on NPP by evaluating topographic effects on primary input data to the model, including both Normalized Difference Vegetation Index (NDVI) and climatic data. A typical green coniferous forest in Yoshino Mountain, Japan, was employed as the study area. The results show that the average NPP is significantly increased after removing topographic influences on NDVI; the average NPP has a relatively minimal change when only topographic effects on climatic data are considered. When both topographic effects on NDVI and climatic data are considered, the average NPP is 1.80 kg m^?2 yr^?1, which is very similar to the ground measurement result of 1.74 kg m^?2 yr^?1.     

静止卫星图像热带气旋云系自动识别

热带气旋对我国东南沿海地区国民经济和人民生命财产威胁巨大,静止卫星云图是热带气旋实时监测的主要数据源。热带气旋在卫星云图上的纹理特征与其它云系相似度高,为气旋云系的自动准确提取带来困难。本文在矢量矩概念的基础上,提出了旋转系数的概念来表征热带气旋的形态本质特征从而实现热带气旋的自动识别。建立了基于静止卫星图像,运用最大类间方差法确定目标云系分割阈值,结合云系面积和亮温分布特性,利用旋转系数进行热带气旋云系自动识别的方法流程。以1211台风海葵为例,在台风生成发展期、成熟期以及消亡期内,进行了改进前后方法识别率的对比实验,统计发现改进方法的识别率分别为76%、95%、78%,均高于原始方法的59%、90%、63%。实验表明改进方法分割的热带气旋云系更为完整,对各阶段的热带气旋云系识别率均更高。     

Estimating woody browse abundance from aerial imagery

In the deciduous forests of the eastern US, timber harvest programmes are often designed to increase the availability of woody browse for terrestrial wildlife. However, assessing the efficacy of timber harvest at increasingly available browse has traditionally required labour‐intensive field‐based measurements of woody plant growth and abundance. The objective of this study was to use readily available digital aerial imagery to estimate the amount of woody browse in regenerating clearcuts in central West Virginia. Aerial imagery from the National Agriculture Imagery Program and woody browse data collected from 11 regenerating clearcuts in the summer of 2007 were used in this analysis. Red, green and blue visible bands, as well as a simple texture metric, were used to create a multiple linear regression model to predict the amount of woody browse. The final model exhibited large correlation (R ² = 0.94) and was statistically significant (F = 22.48, p = 0.0009), indicating that simple measures of image digital numbers and texture have potential utility in assisting forest and wildlife managers to assess habitat quality in forest regeneration areas.     

Assessing light to moderate grazing effects on grassland production using satellite imagery

Understanding the influences of grazing intensity on grassland production is essential for grassland conservation and management improvement. Grazing at light to moderate intensity theoretically enhances grassland production, thus benefiting grassland ecosystems. However, inconsistent results of the beneficial effects of light to moderate grazing on grassland production were reported due to the lack of accurate and repeatable techniques for discriminating grazing effects from other abiotic factors. Advanced remote-sensing techniques provide a promising tool for filling this gap in grazing effects research due to their high spatial and temporal resolution. In this article, the influences of light to moderate grazing on grassland production in mixed grasslands were investigated for the period 1986–2005, using spectral data derived from satellite images. The effects of precipitation on the detection of grazing-induced production change were also analysed. The results revealed that the normalized canopy index (NCI) showed superior performance in quantifying grassland production in mixed grasslands. Significant differences in grassland production between grazed and ungrazed treatments occurred in the three years with above-average and average growing-season precipitations (April–August), but not in the dry years. Most of the variation in production (75%) was explained by growing-season precipitation for both grazed and ungrazed sites. Our results demonstrate the feasibility of using remote-sensing data to monitor long-term light to moderate grazing effects and the important role of precipitation, especially growing-season precipitation, in modulating production in mixed-grassland ecosystems.     

Estimating energy consumption from night-time DMPS/OLS imagery after correcting for saturation effects

A methodology is presented to accurately estimate electric power consumption from saturated night-time Defense Meteorological Satellite Program (DMSP) Operational Linescan System (OLS) imagery using a stable light correction. An area correction for the stable light image of DMSP/OLS for the year 1999 was performed and the build-up area rate data were used to clarify the intensity distribution characteristics of the stable light. Based on the spatial distribution characteristics of the stable light, the saturation light of the electric power supply area of Japan was corrected using a cubic regression equation. The regression between the correction calculations by the cubic regression equation and the statistical electric power consumption data was applied in Japan and also in China, India and 10 other Asian countries. The correction method was then evaluated. This study confirms that electric power consumption can be estimated with high precision from the stable light.     

Application of high spatial resolution satellite imagery for riparian and forest ecosystem classification

Terrestrial Ecosystem Mapping provides critical information to land and resource managers by incorporating information on climate, physiography, surficial material, soil, and vegetation structure. The main objective of this research was to determine the capacity of high spatial resolution satellite image data to discriminate vegetation structural stages in riparian and adjacent forested ecosystems as defined using the British Columbia Terrestrial Ecosystem Mapping (TEM) scheme. A high spatial resolution QuickBird image, captured in June 2005, and coincident field data covering the riparian area of Lost Shoe Creek and adjacent forests on Vancouver Island, British Columbia, was used in this analysis. Semi-variograms were calculated to assess the separability of vegetation structural stages and assess which spatial scales were most appropriate for calculation of grey-level co-occurrence texture measures to maximize structural class separation. The degree of spatial autocorrelation showed that most vegetation structural types in the TEM scheme could be differentiated and that window sizes of 3 × 3 pixels and 11 × 11 pixels were most appropriate for image texture calculations. Using these window sizes, the texture analysis showed that co-occurrence contrast, dissimilarity, and homogeneity texture measures, based on the bands in the visible part of the spectrum, provided the most significant statistical differentiation between vegetation structural classes. Subsequently, an object-oriented classification algorithm was applied to spectral and textural transformations of the QuickBird image data to map the vegetation structural classes. Using both spectral and textural image bands yielded the highest classification accuracy (overall accuracy = 78.95%). The inclusion of image texture increased the classification accuracies of vegetation structure by 2-19%. The results show that information on vegetation structure can be mapped effectively from high spatial resolution satellite image data, providing an additional tool to ongoing aerial photograph interpretation.     

Use of lidar data to simulate forest net primary production

A method has been recently presented to predict the net primary production (NPP) of Mediterranean forests by integrating conventional and remote-sensing data. This method was based on the use of two models, C-Fix and BIOME-BGC, whose outputs are combined with estimates of stem volume and tree age to predict the NPP of the examined ecosystems. This article investigates the possibility of deriving these two forest attributes from airborne high-resolution lidar data. The research was carried out in the San Rossore pine forest, a test site in Central Italy where several investigations have been conducted. First, estimates of stand stem volume and tree age were obtained from lidar data by application of a simplified method based on existing literature and a few ground measurements. The accuracy of these stand attributes was assessed by comparison with the independent ground data derived from a recent forest inventory. Next, the stem volume and tree age estimates were used to drive the NPP modelling strategy, whose outputs were evaluated against the inventory measurements of current annual increment (CAI). The simplified lidar data processing method produces stand stem volume and tree age estimates having moderate accuracy, which are useful to feed the modelling strategy and predict CAI at a stand level. This method's success raises the possibility of integrating ecosystem modelling techniques and lidar data for the simulation of net forest carbon fluxes.     

The consequences of urban land transformation on net primary productivity in the United States

We use data from two satellites and a terrestrial carbon model to quantify the impact of urbanization on the carbon cycle and food production in the US as a result of reduced net primary productivity (NPP). Our results show that urbanization is taking place on the most fertile lands and hence has a disproportionately large overall negative impact on NPP. Urban land transformation in the US has reduced the amount of carbon fixed through photosynthesis by 0.04 pg per year or 1.6% of the pre-urban input. The reduction is enough to offset the 1.8% gain made by the conversion of land to agricultural use, even though urbanization covers an area less than 3% of the land surface in the US and agricultural lands approach 29% of the total land area. At local and regional scales, urbanization increases NPP in resource-limited regions and through localized warming “urban heat” contributes to the extension of the growing season in cold regions. In terms of biologically available energy, the loss of NPP due to urbanization of agricultural lands alone is equivalent to the caloric requirement of 16.5 million people, or about 6% of the US population.