Print version ISSN 0103-9016
Sci. agric. (Piracicaba, Braz.) vol.59 no.4 Piracicaba Oct./Dec. 2002
PHYTOMASS MAPPING OF THE "SERIDÓ CAATINGA" VEGETATION BY THE PLANT AREA AND THE NORMALIZED DIFFERENCE VEGETATION INDECES
Thomaz Corrêa e Castro da Costa1*; Luciano José de Oliveira Accioly2; Maria Ap. José de Oliveira3; Nivaldo Burgos2; Flávio Hugo Barreto Batista da Silva2
1Embrapa Solos, R. Jardim Botânico, 1024 - CEP: 22460-001 - Rio de Janeiro, RJ.
2Embrapa Solos Nordeste, R. Antônio Falcão, 402 - CEP: 51020-240 - Recife, PE.
*Corresponding author <email@example.com>
ABSTRACT: Phytomass is a critical information for economic and environmental activities like the establishment of policies for timber resources, forest management, studies of plant nutrient cycling, CO2 sink, among other. The phytomass of a Caatinga area was obtained by an empirical method using normalized difference vegetation index (NDVI) of Landsat images, the plant area index (PAI) and the phytomass inventory. At a first stage, linear, logarithmic and non-linear models were developed and tested. Bush and tree specimens were considered in the study, so that most of the individuals that contribute to the spectral answer detected by satellite images were included. At a second stage, the orbital parameter NDVI was used to map the PAI, which was used to map the phytomass, based on the relationship of this phytomass as a function of PAI. The residues between measurements and estimates based on NDVI varied from 0 to 84%, while the residues of total dry weight of phytomass per ha obtained by mapping and by dendrometrical equations varied from 5 to 104%, with a large trend of 166 and 448% in open Caatinga areas, due to the contribution of the herbaceous stratum to NDVI.
Key words: CO2, LAI, aboveground biomass, desertification
MAPEAMENTO DA FITOMASSA DA CAATINGA DO SERIDÓ PELOS ÍNDICES DE ÁREA DE PLANTA E DE VEGETAÇÃO DA DIFERENÇA NORMALIZADA
RESUMO: A fitomassa, principalmente arbórea, é informação necessária em atividades econômicas e ambientais, como políticas de uso do recurso madeireiro, manejo florestal, estudos de ciclagem de nutrientes, absorção de CO2, entre outros. A finalidade deste estudo foi a verificação de um método empírico para o mapeamento da fitomassa da Caatinga do Seridó, integrando-se um inventário de fitomassa, o índice de área de planta (IAP) e o índice de vegetação da diferença normalizada (NDVI), por meio de imagens Landsat TM. Na primeira etapa foram desenvolvidos e testados modelos lineares, logarítmicos e não lineares. A abordagem de tamanho foi arbustiva e arbórea, incluindo-se a maior parte dos indivíduos que contribuem na resposta espectral mensurada por imagens de satélite. Em uma segunda etapa utilizamos o parâmetro orbital, NDVI, para o mapeamento do IAP, que por sua vez, foi utilizado para mapear a fitomassa. Os desvios entre mensurações de IAP e estimativas a partir do NDVI, variaram de 0 a 84%, enquanto que os desvios entre Peso Seco Total de Fitomassa por ha obtidos pelo mapeamento e por equações dendrométricas, variaram de 5 a 104%, com grandes tendências de 166 e 448% para áreas de caatinga aberta, provocada pela contribuição do estrato herbáceo no NDVI.
Palavras-chave: CO2, LAI, biomassa aérea, desertificação
Phytomass, mainly tree phytomass, is an information required for economic and environmental activities such as the development of policies for the use of timber resources, forest management, nutrient cycling, and CO2 sink, among other (Brown et al., 1989; Souza & Jesus, 1991; Bohrer & Campos, 1993; Silva et al., 1993a; Scolforo et al., 1993; Foody et al., 1996). Since the destructure assessment of phytomass is the most labor-intensive stage of a forest inventory, and considering also the environmental damages, non-destructive measurement methods are required in this type of study.
One of the methods under study (Pereira, 1985; Santos, 1988; Kazmierczak, 1992; Amaral et al., 1996; Jakubauskas & Price, 1997; Santos et al., 1998; Xavier et al., 1998; Linhares et al., 2000; Qi et al., 2000), is based on empirical relationships between phytomass and radiation measured by orbital sensors, allowing for the use of the plant area index (PAI) as an intermediate variable.
PAI is a better term for the commonly used concept of leaf area index (LAI) measured in the field. This new concept was defined by White et al. (2000) and it includes trunks and branches, which in semi-arid regions will probably represent a more significant contribution to PAI, when compared to other biomes. Thus, LAI can be calculated by subtracting the trunk and branch index (IAG) from PAI, measured in the dry season, when all leaves fall.
Santos et al. (1998) established a relationship between backscattering coeficient (dB) of a JERS-1 image and phytomass values for areas of primary forest and Amazon secondary succession forest, with R2=0.429. This was improved by combining bands 3 and 5 of the Landsat TM sensor, to R2=0.77. Using only Landsat TM images, Amaral et al. (1996), for the secondary succession of abandoned areas of the Amazon forest, established linear relationships between the normalized difference vegetation index (NDVI) and LAI, and NDVI and the basal area (cm2 ha-1), with R2 = 0.77. These results were thought to exceed expectations when compared to those of Box et al. (1989), referred by Amaral et al. (1996), with r = 0.234 between LAI and NDVI, and r=0.344 between NDVI and phytomass. For Pinus, Jakubauskas & Price (1997) obtained R2 = 0.62 for LAI as a function of NDVI, and the greenness component of the transformation Tasseled Cap, and R2= 0.59 for the phytomass as a function of NDVI. Kazmierckak (1992) found for Pinus and Araucaria angustifolia that the best variable for the relationship with orbital parameters was the basal area, with Landsat TM band 5 for both populations, however with a poor adjustment (R2=0.3594 for Araucaria angustifolia, and R2=0.6593 for Pinus).
In semi-arid regions, there are fewer possibility of NVDI saturation, because of the low phytomass density, mainly leaf density, a condition that was not found in more advanced stages of the Amazon forest, in which, its correlation with the LAI is limited in more advanced stages of vegetation succession. Sader et al. (1990) and Honzak et al. (1995) cited by Amaral et al. (1996), showed that the NDVI and the LAI are better correlated when they are obtained in areas in initial succession stages or up to 30 years of age, where lower phytomass values occur. The best NDVI-LAI relationships were obtained by Qi et al. (2000), who adjusted a polynomial model to LAI data collected in a semi-arid region of the State of Arizona, obtaining inflexion for LAI > 1.2, when the model takes the exponential form, with constant NDVI results @ 0.6 for LAI > 2.0. The assessed coefficient of determination was 0.94.
To achieve a good correlation of field measured phytomass with PAI and NDVI, in semi-arid regions, it is necessary to adapt the methodology for tree modeling, with the main purpose of estimating the economic portion of the forest resource for energy, timber or pulp uses of planted forests (Paula Neto et al., 1992; Guimarães & Leite, 1992; Leite and Regazzi 1992; Paula Neto & Rezende, 1992; Silva et al., 1993b) and native forests (Scolforo et al., 1993; Silva et al., 1993a; Souza & Jesus, 1991; Brown et al., 1989). This is due to the fact that the radiation reflected to the sensor is related to the whole canopy (tree, shrub and herbaceous strata) and in the Caatinga there is an increased participation of individuals in the lower stratum because of the low leaf density. Thus, size standards and dendrometrical variables were changed, taking into account shrub and tree sizes, and vegetation complexity ( twisting, forking, and deciduous leaves), and this methodology was compared to the method used by Silva (1998) and Ibama (1992a), in phytomass inventories of caatinga biomass.
The objective of this work was to develop a method to map the shrub-tree phytomass at the Seridó caatinga, targeting at the diagnosis of the distribution of phytomass stock as input for desertification zoning, forest management programs and carbon cycling studies.
This study comprises a stage of various investigations related to land use in the caatinga biome, such as desertification, climate change, and non-sustainable use of timber resources, which have been causing losses of the floristic diversity and soil degradation of the Seridó caatinga ecosystem.
MATERIAL AND METHODS
The pilot area for data collection covers a 30 x 25 km (75000 ha) rectangle, with an average altitude of 270 m ranging from 200 to 700 m, inside the Seridó Desertification Nucleus (Rio Grande do Norte/Paraíba, Brazil), (Figure 1).
Phytomass modeling was broken down into three stages (Figure 2). The final result was phytomass mapping assisted by a geographic information system.
To carry out Phytomass Inventories, the selection of areas for data collection was defined in three levels of densities of the Seridó vegetation (Ibama, 1992a). These levels were selected in the field, to mark transects, phyto-sociological inventory and phytomass collection, also using a spectral similarity criterion (SPRING/INPE segmentation module), of Landsat TM bands 3,4 and 5, totaling a 16 area sample (Table 1).
Quadrant point (Brower et. al, 1998) was the chosen sampling method, using georeferenced 500 m transects (Table 1), with accuracy estimated (EPE) between 3 and 5 m. Distance between points was 10 m and the distance between transects 20 m, totaling 3,469 trees measured and 246 trees cut, representing the most abundant species of Seridó.
To build Phytomass Models by Tree, by Species, we chose to model trees considering forking, resulting in the variable "sum of base circumference" (SCBA), which has a better correlation with phytomass than SCBA2 (sum of CBA2), analytically exact, using as tree inclusion the criterion of base diameter (DBA) larger than or equal to 1 cm. Because of the intermediate standards found in the field, we also chose to use models that are not stratified by area.
The correlations between variables were analyzed as a first check of their effectiveness and of the quality of the equations. Trees occurring in the upper right quadrant, at each point, were standardized as trees for cutting. These trees were weighted in a field scale (0-5 kg and 2.5-15 kg), and their volumes obtained with a xilometer (dimensions: D=1 m, H= 0.7 m, 1 mm graduated ruler, Accuracy=1dm3). Leaf dry weight (DWleaf) was measured using a precision scale.
DWbranch, DWleaf and Vol parameters (dependent variables) were estimated on a per tree basis from dendometrical variables DWbranch, SCBA, HT and DC2 (independent variables). For DWleaf, the variable DATA (Dummy) was also considered for intercept, since data were collected on different dates during the rainy season, once in the stage in which plants had 100% leaves (March) and the second time when the dry season was coming to an end (June), with leaves already starting to fall, giving mean data for May.
Linear and logarithmic models were tested by the Stepwise procedure and non linear for the selection of the equations comparing them to two surveys carried out by Silva (1998) and Ibama (1992a). Criteria used to choose the best linear model were the determination coefficient (R2), variation coefficient of (CV), and residues analysis. Considering the heterogeneity of data because of twisting and forking of individuals, and the shrub inclusion level (DBA ³ 1 cm), high variation coefficients were foreseen, and this was decisive for the selection of the best distribution of residues.
One alternative to the use of logarithmic models to reduce the heterogeneity of variances is the transformation suggested by Box & Cox (1964), referred to and modified by Silva et al. (1993a). This model was used to estimate DWbranch as a function of SCBA and HT, resulting in the following function: DWbranch=[b0(b1+b2 *SCBA+b3*HT+e)+1]1/b0, and b0, b1, b2 and b3 are coefficients to be estimated, and e= error.
Another model modified from the non linear Schumacher and Hall model, resulted from an academic exercise, with added effects of intercept and residue , as a variation of the original model. The same variables, DWbranch, SCBA and HT were also tested by this model, resulting in the function: DWbranch =b0+b1*SCBAb2*HT b3+e, and b0, b1, b2 e b3 are coefficients to estimate, and e=error.
The best equations generated by previous surveys carried out at the caatinga by the PNUD/FAO/Ibama project (Ibama, 1992a) and by Silva (1998) were also applied to data for validation, referring to areas 1, 2 and 3. The main differences between these surveys and the survey developed in this study are presented in Table 2.
Phytomass Estimates by Area as a Function of Dendrometrical Data (DWbranch, DWleaf and Vol) were produced from the equations selected from the methods and surveys described in Table 2, for all areas, supplying basic data for the second stage (Figure 2).
Stages 2 and 3
An LAI-2000 was used to Measure the PAI per Area which calculates the foliar area based on the diffuse energy ratio, B (Reading below the canopy)/A (Reading above the canopy) = T (transmitance), which is proportional to the distance from this energy to the detector, and the amount of leaves per crown volume, and leaf orientation. Readings were made with a 45° view cap at sunrise and sunset and under cloudy skies, according to the following sequence: one reading A, and 10 readings B, and each reading B, measured over transect points, 10 m apart.
For Phytomass Estimates as a function of PAI/NDVI (Phytomass Mapping), a NDVI image was generated with LANDSAT 7 TM bands 3 and 4, passage May/00 (rainy season). NDVI is calculated by: NDVI=(TM4-TM3)/(TM4+TM3), in which TM4 and TM3 are radiance values for band 4 (close infrared ) and band 3 (red). The PAI assumes the function of estimator of the phytomass (independent variable) which estimated by the NDVI (dependent variable) enables specializing the phytomass as a function of the NDVI, according to the functions: (PAI) f (NDVI); and (DWtotal, DWleaf, Vol) f (PAI).
RESULTS AND DISCUSSION
Models to Estimate Phytomass by Tree, by Species
The percentage trend between the estimated DWbranch value and DWtotal, and value measured in the field (real) of the models tested for areas 1, 2 and 3 is presented in Figure 3. The inadequacy of equations generated by Silva (1998) is evidenced from the potential model which, in its theoretical form, is inadequate, since its relationship indicates that a reduction in diameter growth do not reduct the acceleration of the increment in phytomass, which occurs for height, and not for the diameter, since growth in height reduces relatively to the diameter until older plants stabilize.
Equations generated by Ibama (1992a) also had strong trends for areas 1 and 2, even with stratification, probably due to the more restricted approach, as compared to our survey, but with better results than the other methods for area 3. This is due to the DWbranch equation for DWbranch for the Pereiro species, which, in area 3 was cut, presenting significant forking. Since Ibama (1992a) modeled forking and not the tree, the accuracy of derived equations was better than ours for this species.
Among equations adjusted with data of the test areas (Stepwise Method, Modified Schumacher & Hall Model, and Silva et al. (1993a)), no sereve, great variation was found between the R2 and CV results (Table 3). Thus, it was assumed that it was not necessary to measure the product of the larger and smaller crown diameters (DC2) from area 8, because this would make data collection more complex, although this variable improves accuracy significantly.
Main equations generated by the Stepwise method, are described in Table 4. Among DWleaf equations, the DATE variable was not significant for Aspidosperma pyrifolium Mart. (Pereiro), more resistant to leaf fall, and Jatropa pohliana (Pinhão Brabo) which still had most leaves in the second collection. In relation to the DC2 variable, the Pinhão Brabo was the only species that was not included as an explanatory variable for DWleaf.
To estimate the phytomass (Dwleaf and DWbranch) and the volume for the 16 areas a combination of the best equations generated by the Stepwise Method, Schumacher & Hall Modified Model and Box & Cox Transformation, Modified by Silva et al. (1993a) were chosen. For other species, less frequent in this region, with standard not similar to the ones that were modeled, the Tavares et al. (1969) 0.6 morphologic coefficient (shape factor) was adopted.
Phytomass Equations as a Function of PAI/NDVI
Correlations for each area between structural variables of vegetation, phytomass, PAI and NDVI are presented in Table 5. DWleaf is the variable with the largest correlation with PAI and NDVI, because of the almost total contribution of leaves to both indeces. Its relationship with the other structural vegetation variables (0.94 for DWbranch and Vol, 0.79 for DT m-2, 0.78 for B ha-1, and 0.75 for HT), explains the correlation of the other variables with NDVI. In addition, the structural and phytomass variables are better correlated to PAI as compared to NDVI, which justifies the use of PAI as an intermediate variable when estimating the phytomass by orbital parameters.
The resulting equation for DWtotal (DWbranch + DWleaf), having as explanatory variable PAI, is:
DWtotal= -980.47+11.851.25*PAI, R2= 76.48, with data distribution expressed in Figure 4.
One drawback of this equation is its sensitivity to PAI variations, which on the other hand, is estimated as a function of NDVI. This means that consistent results require good adjustments of the PAI and NDVI relationships.
The relationship between PAI as a function of NDVI, obtained for bands 3 and 4 of the Landsat 7 ETM+ Sensor, resulted in the equation:
PAI= 0.6401*exp (2.6929*NDVI), with R2=78.02, represented in Figure 5.
Breakdown of DWtotal in DWbranch and DWleaf can also be obtained as according to the equations:
DWleaf = 147.53*exp (0.9518*PAI) (R2=81.62), expressed in Figure 6.
DWbranch = DWtotal - DWleaf
Vol = 0.001117*DWbranch + 0.3537 (R2=99.97)
The errors that contributed most to the low determination coefficients occurred for the intermediate density areas, which correspond to areas with clearings, with a cluster-type distribution of individuals (A9, A10 and A11), and areas with abundant herbaceous plants, mainly "malva" and "capim panasco" (A2, A5 e A6). Areas with clearings may underestimate the NDVI, by means of soil saturation, while the presence of herbaceous type coverage contributes to NDVI, leading to overestimates for the tree-bush portion.
The evaluation of mapping errors by comparing phytomass estimated as a function of PAI/NDVI to the phytomass estimated based on dendrometrical variables is shown in Table 6. The percentage deviation ranged from 0 (A12) to 84% (A5) for PAI, and from 5% (A1 and A12) to 448% (A6) for DWtotal. For DWleaf, errors ranged from 8% (A13 and A16) to 250% (A6). Vol ranged from 0% (A1) to 432% (A6). The significant deviation in area 6 was confirmed for all variables, since this area had a low tree density, but significant herbaceous coverage represented by malva and capim panasco. Low accuracy is then assessed in phytomass mapping. As the final information is generated in a thematic form, in which data are grouped in classes, accuracy is improved by generalization.
The final product of the Categorized DWtotal map is presented in Figure 7 with the following distribution: 27% of the study area corresponds to the following classes: "pasture", "agriculture", "dam", "riverine vegetation" and "urban zone"; while 73% is covered by caatinga.
Relatively to caatinga phytomass, 13% of the area was found to have 0.1 to 5 Mg ha-1 of Total Dry Weight, in the Northwest portion of the area and 30% had 5 to 10 Mg ha-1. These two phytomass classes, distributed in 43% of the area, occurred mainly over non calcic brown soils that in the past were frequently planted with cotton.
Phytomasses between 10 and 20 Mg ha-1 were found in 21% of the area, and values above 20 Mg ha-1, associated mainly to Dense Caatinga, covered 9% of the area, mainly the Queimadas hills on the Southeast. Areas A7, A8 A15 are in this region, with a large diversity of species and the tallest trees, mostly over litholic soils. The fact that most preserved vegetation was found mainly in the hills is evident and this results from restriction to anthropic activity except for mining activities that are very common in the region.
Phytomass as input for Environmental Studies
Phytomass mapping can be a source of data for the zoning of areas susceptible to desertification, when associated to other parameters such as topography, soil, climate, land use, etc. In the Universal soil loss equation (USLE), Troeh et al. (1991), for example, phytomass enables increased accuracy of the C parameter (Land Use/Coverage factor), replacing the rigid structure per weight of vegetation class.
CO2 Absorption by Caatinga
Taking into account the loss of leaves, since semi-arid climate of the caatinga is a canopy typology in which leaves fall completely in the dry season, the average amount of leaves shed (DWleaf) totals 386 kg ha-1 year-1 for this area, ranging from 701 kg ha-1 year-1 for the Dense Caatinga to 276 kg ha-1 year-1 for the open caatinga, on average.
Using a 48% dry carbon factor in the dry matter (Phillips et al., 1998), the gross rate of the carbon sink for this region, considering only the annual fall of leaves, totals 185 kg C ha-1 year-1 or 678 kg CO2 ha-1 year-1, which corresponds to 26% of the carbon entrapped by the tropical forest, which is 710 ± 340 kg C ha-1 year-1 (Phillips et al., 1998), a considerable proportion for caatinga fragments in a desertification nucleus.
The RN Seridó Region was the object of a comprehensive study carried out in the PNUD/FAO/IBAMA/BRA/87/007 project, part of the Brazilian Northeast integrated forest development program (Ibama, 1992a, b e c), which goes from diagnosis to the definition of strategies and simulation of scenarios. A phytomass mapping in a SIG environment could contribute to this study by providing management and processing support, such as the indication of appropriate areas for forest management, the amount of timber available, among other information that can be readily made available.
Focusing on this pilot area, 33% of the area has a volume of timber above 10 m3 ha-1. If we consider a volume above 15 m3 ha-1 this percentage drops to 17% of the area that would have required conditions for sustainable development, and it is also necessary to deduct non usable timber, environmental interest areas and areas in a topographic situation which is not recommended for forestry activities. .
When relating the mean increment rate taken from Ibama (1992a), totaling 530,45 kg ha-1 year-1 in the Seridó region, or 0.5894 m3 ha-1 year-1, using a mean conversion factor, Vol/DWbranch, we have the gross available volume of 7,409 m3 year-1 in the area with stock above 15 m3 ha-1, which is 12,571 ha. Thus, the remaining volume of this pilot area representing 9% of the whole Seridó nucleus, contributed with 18% of the demand projected for 2001 (Ibama, 1992b).
Canopy Coverage Monitoring
Among main experiences of coverage monitoring and land use in the Brazilian territory, the best known is the "Prodes-Projeto de Estimativa do Desflorestamento da Amazônia Legal", (Legal Amazon Deforesting Estimate Project) aiming at the systematic monitoring of the region and assessment of the evolution of the extension and the annual gross deforesting rate. Phytomass mapping associated with a monitoring system focusing on forest exploitation and agriculture-cattle breeding expansion (Costa et al., 2000) adds important information, besides deforested areas and regenerating areas which correspond to suppressed or regenerated phytomass from the use of the land.
Thus, we deem that even if relationships between orbital and biophysical parameters of caatinga vegetation still lack accuracy, the use of remote sensing to estimate the tree-shrub phytomass makes for improved operation at low cost, and represents an alternative to regional applications.
To CNPq for the invested resources and staff. Dra. Maria do Carmo from IPA, for laboratory support, interns Ana Luiza da Silva, Eudmar da S. Alves and Ricardo A. Irmão for the help in data collection and typing, and, Ademar Barros da Silva and Paulo Cardoso de Lima of Embrapa Solos Nordeste; Rinaldo M. Santos and José Antônio Aleixo da Silva, of for the Department of Forest Engineering at UFRPE for the collaboration in the pre-review of this work.
AMARAL, S.; SOARES, J.V.; ALVES, D.S.; MELLO, E.M.K.; ALMEIDA, S.A.S.; SILVA, O.F.; SILVEIRA, A.M. Relações entre índice de área foliar (LAI), área basal e índice de vegetação (NDVI) em relação a diferentes estágios de crescimento secundário na Floresta Amazônica em Rondônia. (Compact disc). In: SIMPÓSIO BRASILEIRO DE SENSORIAMENTO REMOTO, 8., Salvador, 1996. Anais. Salvador: INPE; SELPER, 1996. [ Links ]
BOHRER, C.B.A.; CAMPOS, A.N.S. Estimativa da fitomassa florestal na Amazônia legal a partir de dados provenientes de inventário florestal In: CONGRESSO FLORESTAL PANAMERICANO, 1.; CONGRESSO FLORESTAL BRASILEIRO, 7., Curitiba, 1993. Anais. São Paulo: SBS; SBEF, 1993. v.2, p.511-512. [ Links ]
BROWER, J.E.; ZAR, J.H.; VON ENDE, C.N. Field and laboratory methods for general ecology. New York: McGraw-Hill, 1998. 273p. [ Links ]
BROWN, S.; GILLESPIE, A.J.R.; LUGO, A.E. Biomass estimation methods for tropical forests with applications to forest inventory data. Forest Science, v.35, p.881-902, 1989. [ Links ]
COSTA, T.C.C.; VARGAS FILHO, R.; COELHO, L.M.; MIRANDA, L.H.F.; LOBATO, J.N. Sistema de monitoramento da cobertura vegetal de Minas Gerais Integração com a fiscalização florestal do Instituto Estadual de Florestas IEF/MG. (Compact disc). In: GIS, 6., Salvador, 2000. Anais. Salvador: Fator GIS, 2000. [ Links ]
FOODY, G.M.; PALUBINSKAS, G.; LUCAS, R.M.; CURRAN, P.J.; HONZAK, M. Identifying terrestrial carbon sinks: classification of succession stages in regenerating tropical forest from Landsat TM data. Remote Sensing of Environment, v.55, p.205-216, 1996. [ Links ]
GUIMARÃES, D.P.; LEITE, H.G. Um novo modelo para descrever o perfil do tronco. Revista Árvore, v.16, p.170-180, 1992. [ Links ]
IBAMA. Plano de manejo florestal para a região do Seridó do Rio Grande do Norte. Natal: IBAMA, 1992a. v.1 (Projeto PNUD/FAO/IBAMA) [ Links ]
IBAMA. Plano de manejo florestal para a região do Seridó do Rio Grande do Norte. Natal: IBAMA, 1992b. v.2 (Projeto PNUD/FAO/IBAMA) [ Links ]
IBAMA. Plano de manejo florestal para a região do Seridó do Rio Grande do Norte. Natal: IBAMA, 1992c. v.3 (Projeto PNUD/FAO/IBAMA) [ Links ]
JAKUBAUSKAS, M.E.; PRICE, K.P. Empirical relationships between structural and spectral factors of Yellowstone lodgepole pine forest. Photogrammetric Engineering and Remote Sensing, v.63, p.1375-1381, 1997. [ Links ]
KAZMIERCZAK, M.L. Estimativa de variáveis dendrométricas em florestas homogêneas a partir de dados LANDSAT TM. São José dos Campos, 1992. 186p. Dissertação (Mestrado) - Instituto Nacional de Pesquisas Espaciais. [ Links ]
LEITE, H.G.; REGAZZI, A.J. Métodos estatísticos para avaliar a igualdade de equações volumétricas. Revista Árvore, v.16, p.59-71, 1992. [ Links ]
LI-COR. LAI-2000 Plant Canopy Analyzer: operating manual. Lincoln: LI-COR, 1992. [ Links ]
LINHARES, C.A.; PONZONI, F.J.; SHIMABUKURO, Y.E. Relação entre volume de madeira e índice de área foliar em um povoamento de Pinus spp.: Estudo de Caso. Revista Árvore, v.24, p.47-54, 2000. [ Links ]
PAULA NETO, F.; NUNES, J.R.S.; VITAL, B.R.; SOUZA, A.L. Equações de volume de casca de Eucalyptus de diferentes idades e condições de local, espécie e método de regeneração. Revista Árvore, v.16, p.157-169, 1992. [ Links ]
PAULA NETO, F.; REZENDE, A.V. Equações de fatores de empilhamento e Tables de volumes em metros estere para árvores individuais. Revista Árvore, v.16, p.72-87, 1992. [ Links ]
PEREIRA, M.D.B. Correlação de fitomassa verde de campo cerrado com dados espectrais obtidos pelo sistema MSS-Landsat e por radiometria de campo. São José dos Campos, 1985. Dissertação (Mestrado) - Instituto Nacional de Pesquisas Espaciais. [ Links ]
PHILLIPS, O.L.; MALHI,Y.; HIGUCHI, N.; LAURANCE, W.F.; NÚÑEZ, P.V.; VÁSQUEZ, R.M.; LAURANCE, S.G.; FERREIRA, L.V.; STERN, M.; BROWN, S.; GRACE, J. Changes in the carbon balance of tropical forests: evidence from long-term plots. Science, v.282, p.439-442, 1998. [ Links ]
QI, J.; KERR, Y.H.; MORAN, M.S.; WELTZ, M.; HUETE, A.R.; SOROOSHIAN; S.; BRYANT, R. Leaf area index estimates using remotely sensed data and BRDF model in a semiarid region. Remote Sensing of Environment, v.73, p.18-30, 2000. [ Links ]
SANTOS, J.R. Fitomassa aérea da vegetação de cerrado: estimativa e correlação com dados do sensor "Thematic Mapper" do satélite LANDSAT. Curitiba, 1988. 156p. Tese (Doutorado) - Universidade Federal do Paraná. [ Links ]
SANTOS, J.R.; LACRUZ, M.S.P.; KEIL, M.; KUX, H.J.H.; XAUD, M.R. Análise da imagem JERS-1 para estimativa da fitomassa aérea de florestas tropicais no sudoeste da Amazônia. (Compact disc). In: SIMPÓSIO BRASILEIRO DE SENSORIAMENTO REMOTO, 9., Santos, 1998. Anais. Santos: INPE; SELPER, 1998. [ Links ]
SCOLFORO, J.R.; LIMA, J.T.; SILVA, S.T. Equações de fitomassa e volume para cerrado senso stricto. In: CONGRESSO FLORESTAL PANAMERICANO, 1.; CONGRESSO FLORESTAL BRASILEIRO, 7., Curitiba, 1993. Anais. São Paulo: SBS; SBEF, 1993. v.2, p.508-510. [ Links ]
SILVA, G.C. Relações alométricas de dez espécies vegetais e estimativas de fitomassas aéreas da Caatinga. Recife, 1998. Tese (Doutorado) - Universidade Federal Rural de Pernambuco. [ Links ]
SILVA, J.A.A.; MELO, M.R.C.S.; BORDERS, B.E. A volume equation for mangrove trees in northeast Brazil. Forest Ecology and Management, v.58, p.129-136, 1993a. [ Links ]
SILVA, J.A.A.; MEUNIER, I.M.J.; BORDERS, B.E.; FARIAS, G.G.A.; ASSUNÇÃO, E. P. Equação volumétrica para Eucalyptus camaldulensis, na região de Barbalha, Ceará, usando o volume da primeira tora como variável independente. Revista Árvore, v.17, p.30-37, 1993b. [ Links ]
SOUZA, A.L.; JESUS, R.M. Equações de volume comercial e fator de forma para espécies da Mata Atlântica ocorrentes na reserva florestal da Companhia Vale do Rio Doce, Linhares, ES. Revista Árvore, v.15, p.257-273, 1991. [ Links ]
TAVARES, S.; PAIVA, F.A.F.; TAVARES, E.J.S.; LIMA, J.L.S. Inventário florestal do Ceará. Estudo preliminar das matas remanescentes do Município de Quixadá. Boletim de Recursos Naturais, v.7, p.93-111, 1969. [ Links ]
TROEH, F.R.; HOBBS, J.A.; DONAHUE, R.L. Soil and conservation. New Jersey: Prentice-Hall, 1991. 530p. [ Links ]
WHITE, M.A.; ASNER, G.P.; NEMANI, R.R.; PRIVETTE, J.L.; RUNNING, S.W. Measuring fractional cover and leaf index in arid ecosystems: digital camera, radiation transmittance, and laser altimetry methods. Remote Sensing of Environment, v.74, p.45-57, 2000. [ Links ]
XAVIER, A.C.; SOARES, J.V.; ALMEIDA, A.C.; FREITAS, C.C. Estimativa de LAI de plantações de eucaliptos a partir de dados TM/LANDSAT. (Compact disc). In: SIMPÓSIO BRASILEIRO DE SENSORIAMENTO REMOTO, 9., Santos, 1998. Anais. Santos: INPE; SELPER, 1998. [ Links ]
Received February 16, 2001