Abstract
Biomass is the entire foundation of energetic in ecosystems and the ultimate food source for human being, biomass estimation, therefore, provides a sound basis for monitoring the Earth's environment and evaluating human-induced impacts on the environmental change.
This paper presents a promising approach to estimating biomass by integrating satellite data and carbon dynamics model. Anural actual net primary productivity (NPP) is first estimated with monthly composite 1-degre AVHRR Normalized Difference Vegetation Index (NDIV) data using the production efficiency approach in which conversion efficiency of canopy absorbed Phtosynthetically active radiation into primary production is considered. NPP estimates are subsequently incorporated into primary production is considered. NPP estimates are subsequently incorporated into a carbon dynamics model, PHYTOMASS model, to simulate biomass accumulation over succession until equilibrium with climate. Ag global map to terrestrial biomass is finally generated based on the estimation results after their validation with field measurements.
Background
Biomass is the total mass of all living organisms, most of it on the Earth is produced by green plants through photosynthesis. Biomass is of fundamental significance in ecosystems, it provides the entire basis of energy flow and food chain. Biomass is also vital for human being the largest portion of our food supply is from the production of plants.
Due to land overuse and deforestation, however, biomass loss and carbon release, in general, has become very serious; this has led to severe environmental problems such as global warming and it might even threaten human's habitation on the Earth in the log run.
Better understanding of biomass status quo and its change due to human activities is, therefore, important and urgent for suitable land use and natural resource conservation. However, very few researches have been conducted on biomass estimation especially on a global scale. On the other hand, many efforts have made in estimating net primary productivity(NPP), that is, the rate of biomass increment, using regression equations and carbon models,. Remote sensing and satellite data provides new capability of estimating plant production. Recently some researchers relate above-ground biomass to satellite spectral reflectance (Prince et l. 1995); however, this kind of relationship holds weak theoretical basis and is thus far from satisfactory. Furthermore, a global map of biomass is still lacing, though biomass data is often required for global research.
Literature Review and Research Approach
Mainly two kinds of approaches to biomass estimation have been proposed in literature, those are, satellite-based approach and model-based approach.
Researches have shown that spectral vegetation indices derived from spectral signals, such as Normalized Difference Vegetation Index (NDVI), is highly related to several vegetation phenomena, including leaf area index, phonology, and net primary production (Gowrd and Dye, 1987; Box et al., 1989). However, biomass does not have strong relationship with vegetation indices in theory, because nowadays commonly used optical sensor data, like NOAA vegetation indices in theory, because nowadays commonly used optical sensor data, like NOAA AVHRR data, normally reflect the features of plant canopy, not much that of plant woody AVHRR data, normally reflect the features of pant canopy, not much that of plant woody components. Microwave observations(e.g. , SAR data) offer potential for estimating biomass and thus appear promising, since microwave, unlike optical sensors, can penetrate through the
canopy into in to woody parts and provide their spectral information. Nonetheless, the current limited spatial coverage of microwave data and its suffering from topographic effects have so far been obstacles to global application.
On the other hand, biomass is the assimilation and accumulation of carbon into organic mater determined by carbon fixation and carbon release, therefore, biomass can be estimated by modeling carbon dynamics. However, most carbon models, present are primarily driven by climate parameters so that the estimation results give a sense of long-term potentials and can not really reflect the current time biomass status and biomass change induced by human activities.
In this research, an unique approach is proposed to estimating biomass by integrating both satellite -and model -based approaches. This research approach can be divided into two steps: first NPP is estimated with AVHRR NDVI satellite data using production efficiency approach then biomass is simulated with a carbon dynamics model, PHYTOMASS model.
Production Efficiency Approach to NPP Estimation
Monteith (1972) first introduced the concept of "production efficiency"; he suggested that plant production is strongly correlated with the amount of photosynthetically active radiation absorbed or intercepted by plant canopy eventually conversed in to organic matter. Researches have also indicated that vegetation indices, such as NDVI, are highly related to net primary productivity (Goward et al., 1985). Further researches (Gowad and Dye, 1987; Ruimy et al. 1994) have subsequently shown that the fraction of canopy-absorbed radiation (fapar) can be linearly related to NDVI,. Therefore, NPP can be estimated with NDVI in the form:
NPP = Î*å(APAR) =
Î*å( fapar *PAR) =
Î*å ( fapar (NDVI)*PAR)
Where ,
NDVI is normalized difference vegetation index,
PAR is Phtosynthetically active radiation (MJ/m2)
APAR is photosynthetic ally active radiation absorbed by canopy (MJ/m
2)
å (APAAR) is the annual sum of APAR ((MJ/m
2)
NPP is net primary productivity (g/m
2/year)
fapar is the APAR fraction of PAR,
Î is conversing efficiency of APAR into primary production (g/MJ).
In this research, conversion efficiency Î) is treated constant, 1.5g/MJ, based on Monteith's (1977) conclusion and Ruimy et al.'s (1994) experimental result. And the fapar -NDVI relationship equation follows that suggested by Ruimy et al. (1994), which was derived from corrected and filtered NDVI:
fapar = -0/+*NDVI
PHYTHOMASS Model for Biomass Simulation
PHYTOMAS, developed BY Box (unpublished, in progress), is a carbon dynamics model, simulating biomass accumulation and related biometabolism over vegetation succession unit it reaches a climatic steady state. The underlying principle of this model is that the production processes, at broad scales, are controlled mainly by ambient climatic conditions and biomass accumulation is determined primarily by the climatically driven production-respiration balance.
In PHTOMASS, metabolism is based on globally calibrated annual models portioned bioclimatically for monthly potentials which are then modified by stochastic growth models to give running increments and accumulations of standing biomass. The model can be schematically expressed by the following equation.
B-Bseed +
å(DB) = Bseed +
å(NPP-LF-H) = Bseed +
å(GPP-R) - LF -H)
Where
GPP is gross primary productivity (g/m
2/month),
R is respiration (g/m
2/month),
NPP is net primary productivity (g/m
2/month), equal to GPP-R,
LF is litter fall(g/m
2/month),
H is harvest or herbivore (g/m
2/month),
B is biomass (g/m
2),
B
seed is initial biomass((g/m
2),
DB increment of living biomes (g/m
2/month), equal to NPP-LF-H,
å(DB) is total amount of biomes increments over succession until equilibrium (g/m
2)
In the model, initial biomass is set to 10/2 to start simulation. The model runs on a monthly time step; input data required are only lone-term monthly average temprature and precipitation plus elevation data.
The PHYTOMASS model, based on climatic data, was originally written to predict potential biomass, so it is modified to estimate actual biomass by substituting satellite-based NPP estimate for potential NPP estimates for potential NPP outputs produced in running the model and determining the ratio of harvest/herbivory to NPP at each site with respect to its intensity of agricultural cultivation.
Research Framework
Research framework of integration of production efficiency approach and PHYTOMASS model is illustrate in Figure 1.

Figure Research Framework of Integration of Production Efficiency Approach and PHYTOMASS Model
Data
Data required include NDVI and PAR data for NPP estimation, and average monthyly temperature and precipitation, elevation data and global map of cultivation intensity for biomass simulation.
Monthly 1-degree AVHRR NDVI data set for 1989 is used in this research, accessed from the Internet website of NOAA Distributed Active Archive Centere (DAAC). PAR data for 1989 is from the global solar radiation data set developed by Eck and Dye (1991); spatial coverage is within 66o North /South latitude. Leemans and Cramer mean monthly temperature and precipitation data sets are for input in running the PHYTOMAS model.
Since these data sets come from difficult sources as different spatial resolution, all the data sets are thus converted into an uniform resolution, 0.5 degree, in laitude /longitude projection for data integrity and data comparability.