Net primary productivity of China's terrestrial ecosystems from a process model driven by remote sensing

The terrestrial carbon cycle is one of the foci in global climate change research. Simulating net primary productivity (NPP) of terrestrial ecosystems is important for carbon cycle research. In this study, China's terrestrial NPP was simulated using the Boreal Ecosystem Productivity Simulator (BEPS), a carbon-water coupled process model based on remote sensing inputs. For these purposes, a national-wide database (including leaf area index, land cover, meteorology, vegetation and soil) at a 1 km resolution and a validation database were established. Using these databases and BEPS, daily maps of NPP for the entire China's landmass in 2001 were produced, and gross primary productivity (GPP) and autotrophic respiration (RA) were estimated. Using the simulated results, we explore temporal-spatial patterns of China's terrestrial NPP and the mechanisms of its responses to various environmental factors. The total NPP and mean NPP of China's landmass were 2.235 GtC and 235.2 gCm(-2)yr(-1), respectively; the total GPP and mean GPP were 4.418 GtC and 465 gCm(-2)yr(-1); and the total RA and mean RA were 2.227 GtC and 234 gCm(-2)yr(-1), respectively. On average, NPP was 50.6% of GPP. In addition, statistical analysis of NPP of different land cover types was conducted, and spatiotemporal patterns of NPP were investigated. The response of NPP to changes in some key factors such as LAI, precipitation, temperature, solar radiation, VPD and AWC are evaluated and discussed.     

Observation and simulation of net primary productivity in Qilian Mountain, western China

We modeled net primary productivity (NPP) at high spatial resolution using an advanced spaceborne thermal emission and reflection radiometer (ASTER) image of a Qilian Mountain study area using the boreal ecosystem productivity simulator (BEPS). Two key driving variables of the model, leaf area index (LAI) and land cover type, were derived from ASTER and moderate resolution imaging spectroradiometer (MODIS) data. Other spatially explicit inputs included daily meteorological data (radiation, precipitation, temperature, humidity), available soil water holding capacity (AWC), and forest biomass. NPP was estimated for coniferous forests and other land cover types in the study area. The result showed that NPP of coniferous forests in the study area was about 4.4 tCha(-1)y(-1). The correlation coefficient between the modeled NPP and ground measurements was 0.84, with a mean relative error of about 13.9%.     

A Comparison of Spatial and Spectral Image Resolution for Mapping Invasive Plants in Coastal California

We explored the potential of detecting three target invasive species: iceplant (Carpobrotus edulis), jubata grass (Cortaderia jubata), and blue gum (Eucalyptus globulus) at Vandenberg Air Force Base, California. We compared the accuracy of mapping six communities (intact coastal scrub, iceplant invaded coastal scrub, iceplant invaded chaparral, jubata grass invaded chaparral, blue gum invaded chaparral, and intact chaparral) using four images with different combinations of spatial and spectral resolution: hyperspectral AVIRIS imagery (174 wavebands, 4 m spatial resolution), spatially degraded AVIRIS (174 bands, 30 m), spectrally degraded AVIRIS (6 bands, 4 m), and both spatially and spectrally degraded AVIRIS (6 bands, 30 m, i.e., simulated Landsat ETM data). Overall success rates for classifying the six classes was 75% (kappa 0.7) using full resolution AVIRIS, 58% (kappa 0.5) for the spatially degraded AVIRIS, 42% (kappa 0.3) for the spectrally degraded AVIRIS, and 37% (kappa 0.3) for the spatially and spectrally degraded AVIRIS. A true Landsat ETM image was also classified to illustrate that the results from the simulated ETM data were representative, which provided an accuracy of 50% (kappa 0.4). Mapping accuracies using different resolution images are evaluated in the context of community heterogeneity (species richness, diversity, and percent species cover). Findings illustrate that higher mapping accuracies are achieved with images possessing high spectral resolution, thus capturing information across the visible and reflected infrared solar spectrum. Understanding the tradeoffs in spectral and spatial resolution can assist land managers in deciding the most appropriate imagery with respect to target invasives and community characteristics.     

Environmental mapping based on spatial variability

Environmental maps show the probable environmental states of different types of land use or development of landscape in a geographic context. Remotely sensed data are particularly efficient for environmental mapping in order to outline major environmental types. Multiple schemes of image classification used in environmental mapping are either traditionally statistical or heuristic. While the former methods do not take account of spatial variability in space and aerial data, the latter ones does not lend themselves to optimal solutions we present. Novel probabilistic models of piecewise-homogeneous images are used in environmental mapping to segment real images. The models consider both an image and a land cover map. Such a pair constitutes an example of a Markov random field specified by a joint Gibbs probability distribution of images and maps. Parameters of the model are estimated by using a stochastic approximation technique. Its convergence to the desired values is studied experimentally. Addition of spatial attributes appears to be necessary in most areas where the differences in spatial data between regions in the image occur. Experiments in generating the pairs of images and environmental maps and in segmenting the simulated as well as real images are discussed.     

Spatial scaling between leaf area index maps of different resolutions

We developed algorithms for spatial scaling of leaf area index (LAI) using sub-pixel information. The study area is located near Liping County, Guizhou Province, in China. Methods for LAI spatial scaling were investigated on LAI images with 960 m resolution derived in two ways. LAI from distributed calculation (LAID) was derived using Landsat ETM+ data (30 m), and LAI from lumped calculation (LAIL) was obtained from the coarse (960 m) resolution data derived through resampling the ETM+ data. We found that lumped calculations can be considerably biased compared to the distributed (ETM+) case, suggesting that global and regional LAI maps can be biased if surface heterogeneity within the mapping resolution is ignored. Based on these results, we developed algorithms for removing the biases in lumped LAI maps using sub-pixel land cover-type information, and applied these to correct one coarse resolution LAI product which greatly improved its accuracy.     

Assessing the impact of urbanization on regional net primary productivity in Jiangyin County, China

Urbanization is one of the most important aspects of global change. The process of urbanization has a significant impact on the terrestrial ecosystem carbon cycle. The Yangtze Delta region has one of the highest rates of urbanization in China. In this study, carried out in Jiangyin County as a representative region within the Yangtze Delta, land use and land cover changes were estimated using Landsat TM and ETM+ imagery. With these satellite data and the BEPS process model (Boreal Ecosystem Productivity Simulator), the impacts of urbanization on regional net primary productivity (NPP) and annual net primary production were assessed for 1991 and 2002. Landsat-based land cover maps in 1991 and 2002 showed that urban development encroached large areas of cropland and forest. Expansion of residential areas and reduction of vegetated areas were the major forms of land transformation in Jiangyin County during this period. Mean NPP of the total area decreased from 818 to 699 gCm(-2)yr(-1) during the period of 1991 to 2002. NPP of cropland was only reduced by 2.7% while forest NPP was reduced by 9.3%. Regional annual primary production decreased from 808 GgC in 1991 to 691 GgC in 2002, a reduction of 14.5%. Land cover changes reduced regional NPP directly, and the increasing intensity and frequency of human-induced disturbance in the urbanized areas could be the main reason for the decrease in forest NPP.     

Effects of topography on simulated net primary productivity at landscape scale

Local topography significantly affects spatial variations of climatic variables and soil water movement in complex terrain. Therefore, the distribution and productivity of ecosystems are closely linked to topography. Using a coupled terrestrial carbon and hydrological model (BEPS-TerrainLab model), the topographic effects on the net primary productivity (NPP) are analyzed through four modelling experiments for a 5700 km(2) area in Baohe River basin, Shaanxi Province, northwest of China. The model was able to capture 81% of the variability in NPP estimated from tree rings, with a mean relative error of 3.1%. The average NPP in 2003 for the study area was 741 gCm(-2)yr(-1) from a model run including topographic effects on the distributions of climate variables and lateral flow of ground water. Topography has considerable effect on NPP, which peaks near 1350 m above the sea level. An elevation increase of 100 m above this level reduces the average annual NPP by about 25 gCm(-2). The terrain aspect gives rise to a NPP change of 5% for forests located below 1900 m as a result of its influence on incident solar radiation. For the whole study area, a simulation totally excluding topographic effects on the distributions of climatic variables and ground water movement overestimated the average NPP by 5%.     

Relationships between net primary productivity and stand age for several forest types and their influence on China's carbon balance

Affected by natural and anthropogenic disturbances such as forest fires, insect-induced mortality and harvesting, forest stand age plays an important role in determining the distribution of carbon pools and fluxes in a variety of forest ecosystems. An improved understanding of the relationship between net primary productivity (NPP) and stand age (i.e., age-related increase and decline in forest productivity) is essential for the simulation and prediction of the global carbon cycle at annual, decadal, centurial, or even longer temporal scales. In this paper, we developed functions describing the relationship between national mean NPP and stand age using stand age information derived from forest inventory data and NPP simulated by the BEPS (Boreal Ecosystem Productivity Simulator) model in 2001. Due to differences in ecobiophysical characteristics of different forest types, NPP-age equations were developed for five typical forest ecosystems in China (deciduous needleleaf forest (DNF), evergreen needleleaf forest in tropic and subtropical zones (ENF-S), deciduous broadleaf forest (DBF), evergreen broadleaf forest (EBF), and mixed broadleaf forest (MBF)). For DNF, ENF-S, EBF, and MBF, changes in NPP with age were well fitted with a common non-linear function, with R(2) values equal to 0.90, 0.75, 0.66, and 0.67, respectively. In contrast, a second order polynomial was best suitable for simulating the change of NPP for DBF, with an R(2) value of 0.79. The timing and magnitude of the maximum NPP varied with forest types. DNF, EBF, and MBF reached the peak NPP at the age of 54, 40, and 32 years, respectively, while the NPP of ENF-S maximizes at the age of 13 years. The highest NPP of DBF appeared at 122 years. NPP was generally lower in older stands with the exception of DBF, and this particular finding runs counter to the paradigm of age-related decline in forest growth. Evaluation based on measurements of NPP and stand age at the plot-level demonstrates the reliability and applicability of the fitted NPP-age relationships. These relationships were used to replace the normalized NPP-age relationship used in the original InTEC (Integrated Terrestrial Ecosystem Carbon) model, to improve the accuracy of estimated carbon balance for China's forest ecosystems. With the revised NPP-age relationship, the InTEC model simulated a larger carbon source from 1950-1980 and a larger carbon sink from 1985-2001 for China's forests than the original InTEC model did because of the modification to the age-related carbon dynamics in forests. This finding confirms the importance of considering the dynamics of NPP related to forest age in estimating regional and global terrestrial carbon budgets.     

Modeling Streams and Hydrogeomorphic Attributes in Oregon From Digital and Field Data1

Abstract: Managers, regulators, and researchers of aquatic ecosystems are increasingly pressed to consider large areas. However, accurate stream maps with geo‐referenced attributes are uncommon over relevant spatial extents. Field inventories provide high‐quality data, particularly for habitat characteristics at fine spatial resolutions (e.g., large wood), but are costly and so cover relatively small areas. Recent availability of regional digital data and Geographic Information Systems software has advanced capabilities to delineate stream networks and estimate coarse‐resolution hydrogeomorphic attributes (e.g., gradient). A spatially comprehensive coverage results, but types of modeled outputs may be limited and their accuracy is typically unknown. Capitalizing on strengths in both field and regional digital data, we modeled a synthetic stream network and a variety of hydrogeomorphic attributes for the Oregon Coastal Province. The synthetic network, encompassing 96,000 km of stream, was derived from digital elevation data. We used high‐resolution but spatially restricted data from field inventories and streamflow gauges to evaluate, calibrate, and interpret hydrogeomorphic attributes modeled from digital elevation and precipitation data. The attributes we chose to model (drainage area, mean annual precipitation, mean annual flow, probability of perennial flow, channel gradient, active‐channel width and depth, valley‐floor width, valley‐width index, and valley constraint) have demonstrated value for stream research and management. For most of these attributes, field‐measured, and modeled values were highly correlated, yielding confidence in the modeled outputs. The modeled stream network and attributes have been used for a variety of purposes, including mapping riparian areas, identifying headwater streams likely to transport debris flows, and characterizing the potential of streams to provide high‐quality habitat for salmonids. Our framework and models can be adapted and applied to areas where the necessary field and digital data exist or can be obtained.     

A multi-scale remote sensing approach for monitoring northern peatland hydrology: present possibilities and future challenges

Remote sensing has the potential to provide quantitative spatially explicit hydrological information across northern peatland complexes. This paper details a multi-scale remote sensing approach for assessing the use of Sphagnum mosses as proxy indicators of near-surface hydrology. Several spectral indices developed from the near infra-red (NIR) and shortwave infra-red (SWIR) liquid water absorption bands, as well as a biophysical index can be correlated with measures of near-surface moisture in the laboratory, in the field and from airborne imagery. Data from all platforms revealed similar patterns in the spectral indices in relation to changes in moisture although the strength of correlations was reduced as the spatial scale increased. The rapid collection of temporally and spatially explicit hydrological data means that the technique has potential practical application for environmental managers and peatland scientists at the local scale. The task of up-scaling the technique for use in operational peatland hydrological monitoring to the global scale is challenging but achievable, and requires further investigation into the heterogeneity of near-surface moisture across Sphagnum patches and the application of novel image processing techniques to improve the spatial resolution of currently available satellite imagery.     

Effect of the Spatial Variability of Land Use,Soil Type,and Precipitation on Streamflows in Small Watersheds1

Abstract: The spatial variability of the data used in models includes the spatial discretization of the system into subsystems, the data resolution, and the spatial distribution of hydrologic features and parameters. In this study, we investigate the effect of the spatial distribution of land use, soil type, and precipitation on the simulated flows at the outlet of “small watersheds” (i.e., watersheds with times of concentration shorter than the model computational time step). The Soil and Water Assessment Tool model was used to estimate runoff and hydrographs. Different representations of the spatial data resulted in comparable model performances and even the use of uniform land use and soil type maps, instead of spatially distributed, was not noticeable. It was found that, although spatially distributed data help understand the characteristics of the watershed and provide valuable information to distributed hydrologic models, when the watershed is small, realistic representations of the spatial data do not necessarily improve the model performance. The results obtained from this study provide insights on the relevance of taking into account the spatial distribution of land use, soil type, and precipitation when modeling small watersheds.     

MONITORING SUSPENDED PARTICULATE MATTER IN PUERTO RICO: FIELD MEASUREMENTS AND REMOTE SENSING1

ABSTRACT: The spatial distribution of suspended particulatematter (SPM) was estimated in Mayagüez Bay on the west coast of Puerto Rico by using traditional ship board measurements and remotely sensed data acquired over four days during January 1990. This effort was part of a joint project between NASA and the University of Puerto Rico to develop techniques to monitor the water quality of a Caribbean coastal zone. This paper presents the methods and algorithms developed to map and analyze short term changes in the source and spatial distribution of SPM in Mayagüez Bay by using remotely sensed data acquired by the Calibrated Airborne Multispectral Scanner (CAMS). A PC-based data acquisition system was developed to record continuous ship measurements of select in- water variables. Spectral reflectances derived from CAMS red and near-IR data were corrected for atmospheric effects and then used to generate maps of SPM. These maps displayed SPM plumes associated with the mouths of the bay's three river systems. Significant day-to-day differences in the spatial characteristics were observed, suggesting that changes in river discharge occurred. However, an analysis of estimated river discharge, sediment yield, local wind field, and thermal river plume indicates that observed sediment plumes result primarily from wind-driven resuspension events.     

Selecting Sites for Converting Farmlands to Wetlands in the Sanjiang Plain, Northeast China, Based on Remote Sensing and GIS

Wetlands in the Sanjiang Plain are rich in biodiversity and natural resources in the northeast of China. However, this wetland area has decreased in size and deteriorated in quality owing to expanded agricultural activities since the 1950s. Converting farmlands to wetlands is necessary to improve these conditions. Using Remote Sensing (RS) and Geographic Information Systems (GIS) technologies, we derived farmland productivity data and hydrology data for the Sanjiang Plain. The farmland productivity data were derived from land use and net primary productivity (NPP) data of the MODIS products. We obtained three productivity farmland classes (low, medium, and high) through the NPP anomaly percentage method. We were only concerned with the low-productivity farmland. Hydrology data were modeled with a wetness index, which was derived from Digital Elevation Model (DEM) data. Based on these two data layers, we identified and prioritized sites for the conversion of farmlands to wetlands. The areas with low farmland productivity and medium or high wetness values have potential to support the conversion of farmlands to wetlands. Potential sites were prioritized in terms of patch size and proximity to natural wetlands and water bodies. We obtained three priority classes, among which the high-priority class would be used as the areas for the recent conversion of farmlands to wetlands. The area of this class was 75,888 ha and accounted for 1.3% of the total farmland area.     

Land use changes and net primary production in the Georgia,USA, landscape: 1935–1982

Land use since 1935 was quantified for Georgia, USA, and for a sample of 20 counties from the major physiographic regions within the state. Statistical data on crop production, pasture productivity, and forest growth were used to estimate net primary production. Appropriate harvest indices (ratio of crop yield to total plant production) were used to correct crop yield data for different decades. Net primary production (NPP) of the Georgia landscape increased from 2.5 to 6.4 tonnes/ha from 1935 to 1982, but varied considerably among land uses and physiographic regions. NPP in the piedmont and mountains reached a plateau between 1960 and 1982, but the upper and lower coastal plains showed a continued linear increase in NPP. In all regions, NPP rose most between 1960 and 1982, coinciding with increases in inputs such as fertilizer and irrigation. Natural ecosystem NPP for Georgia is approximately 16–18 tonnes/ha, and the estimated actual NPP is thus considerably less than the potential. Spatial and temporal patterns of NPP may be a useful basis for evaluating the biological performance of a landscape.     

Techniques for accuracy assessment of tree locations extracted from remotely sensed imagery

Remotely sensed imagery is becoming a common source of environmental data. Consequently, there is an increasing need for tools to assess the accuracy and information content of such data. Particularly when the spatial resolution of imagery is fine, the accuracy of image processing is determined by comparisons with field data. However, the nature of error is more difficult to assess. In this paper we describe a set of tools intended for such an assessment when tree objects are extracted and field data are available for comparison. These techniques are demonstrated on individual tree locations extracted from an IKONOS image via local maximum filtering. The locations of the extracted trees are compared with field data to determine the number of found and missed trees. Aspatial and spatial (Voronoi) analysis methods are used to examine the nature of errors by searching for trends in characteristics of found and missed trees. As well, analysis is conducted to assess the information content of found trees.     

Satellite-based indices in the analysis of land cover for municipalities in the province of Siena, Italy

Territorial indicators based on the satellite measured reflected or emitted energy can provide valuable information on the spatial evolution of a territory. Information regarding land use, biomass coverage and radiant temperature were obtained from remotely obtained measurements for each municipality in the Province of Siena. Remotely sensed data were combined into different indices to compare characteristics of land cover between territories and create an information base for continued study. A classification was performed on a municipal level and land use classes were grouped together and a general index of use pressure was created. A vegetation index was used to compare biomass densities. A radiant temperature index was calculated using measure thermal infrared emissions. The results of these analyses allowed researchers to examine the intra and inter-municipal spatial heterogeneity within the Province. A comparison between the results of each index was made and the problems related to the use of such indices were examined.     

OPTICAL DATA PROCESSING AND PROJECTED APPLICATIONS OF THE ERTS-1 IMAGERY COVERING THE 1973 MISSISSIPPI RIVER VALLEY FLOODS1

ABSTRACT. Flooding along the Mississippi River and some of its tributaries was detected by the multispectral scanner (MSS) on the Earth Resources Technology Satellite (ERTS-1) on at least three orbits during the spring of 1973. The ERTS data provided the first opportunity for mapping the regional extent of flooding at the time of the imagery. Special optical data processing techniques were used to produce a variety of multispectral color composites enhancing flood-plain details. One of these, a 2-color composite of near infrared bands 6 and 7, was enlarged and registered to 1:250,000-scale topographic maps and used as the basis for preparation of flood image maps. Two specially filtered 3-color composites of MSS bands 5, 6, and 7 and 4, 5, and 7 were prepared to aid in the interpretation of the data. The extent of the flooding was vividly depicted on a single image by 2-color temporal composites produced on the additive-color viewer using band 7 flood data superimposed on pre-flood band 7 images. On May 24, when the floodwaters at St. Louis receded to bankfull stage, imagery was again obtained by ERTS. Analysis of temporal data composites of the pre-flood and post-flood band 7 images indicate that changes in surface reflectance characteristics caused by the flooding can be delineated, thus making it possible to map the overall area flooded without the necessity of a real-time system to track and image the peak flood waves. Regional planning and disaster relief agencies such as the Corps of Engineers, Office of Emergency Preparedness, Soil Conservation Service, interstate river basin commissions and state agencies, as well as private lending and insurance institutions, have indicated strong potential applications for ERTS image-maps of flood-prone areas.     

Building spectral libraries for wetlands land cover classification and hyperspectral remote sensing

Recent advances in remote sensing provide opportunities to map plant species and vegetation within wetlands at management relevant scales and resolutions. Hyperspectral imagers, currently available on airborne platforms, provide increased spectral resolution over existing space-based sensors that can document detailed information on the distribution of vegetation community types, and sometimes species. Development of spectral libraries of wetland species is a key component needed to facilitate advanced analytical techniques to monitor wetlands. Canopy and leaf spectra at five sites in California, Texas, and Mississippi were sampled to create a common spectral library for mapping wetlands from remotely sensed data. An extensive library of spectra (n=1336) for coastal wetland communities, across a range of bioclimatic, edaphic, and disturbance conditions were measured. The wetland spectral libraries were used to classify and delineate vegetation at a separate location, the Pacheco Creek wetland in the Sacramento Delta, California, using a PROBE-1 airborne hyperspectral data set (5m pixel resolution, 128 bands). This study discusses sampling and collection methodologies for building libraries, and illustrates the potential of advanced sensors to map wetland composition. The importance of developing comprehensive wetland spectral libraries, across diverse ecosystems is highlighted. In tandem with improved analytical tools these libraries provide a physical basis for interpretation that is less subject to conditions of specific data sets. To facilitate a global approach to the application of hyperspectral imagers to mapping wetlands, we suggest that criteria for and compilation of wetland spectral libraries should proceed today in anticipation of the wider availability and eventual space-based deployment of advanced hyperspectral high spatial resolution sensors.     

ASSESSMENT OF REMOTE SENSING INPUT TO HYDROLOGIC MODELS1

Remotely sensed variables such as land cover type and snow-cover extent can currently be used directly and effectively in a few specific hydrologic models. Regression models can also be developed using physiographic and snow-cover data to permit estimation of discharge characteristics over extended periods such as a season or year. Most models, however, are not of an appropriate design to readily accept as input the various types of remote sensing parameters that can be obtained now or in the future. Because this new technology has the potential for producing hydrologic data that has significant information content on an areal basis, both inexpensively and repetitively, effort should be devoted now to either modifying existing models or developing new models that can use these data. Minor modifications would at least allow the remote sensing data to be used in an ancillary way to update the model state variables, whereas major structural modifications or new models would permit direct input of the data through remote sensing compatible algorithms. Although current remote sensing inputs to hydrologic models employ only visible and near infrared data, model modification or development should accommodate microwave and thermal infrared data that will be more widely available in the future.     

The Effects of Urbanization on Net Primary Productivity in Southeastern China

Net primary productivity (NPP) is one of the major ecosystem products on which human societies rely heavily. However, rapid urban sprawl and its associated dense population and economic conditions have generated great pressure on natural resources, food security, and environments. It is valuable to understand how urban expansion and associated demographic and economic conditions affect ecosystem functions. This research conducted a case study in Southeastern China to examine the impacts of urban expansion and demographic and economic conditions on NPP. The data sources used in research include human settlement developed through a combination of MODIS, DMSP-OLS and Landsat ETM+ images, the annual NPP from MODIS, and the population and gross domestic product (GDP) from the 2000 census data. Multiple regression analysis and nonlinear regression analysis were used to examine the relationships of NPP with settlement, population and GDP. This research indicates that settlement, population and GDP have strongly negative correlation with NPP in Southeastern China, but the outcomes were nonlinear when population or GDP reached certain thresholds.