The results for the area surveyed, presented hereafter, have been based on the interpretation of the entire recent image in each sampling unit. Sampling unit area in this case does not correspond to that of the change matrices discussed in the previous sections. By referring to the entire recent image the land area observed amounts to 324.05 million hectares, representing an actual sampling intensity of 10.6 percent, versus the 8.2 percent intensity represented by the change matrices.
The resulting forest cover of each SU has been subsequently adjusted to the years 1980 and 1990 by applying the rate of change determined in the area common to both historical and recent images for each specific forest definition. Definitions of forest and their implications are discussed in Section 2.2.1.
In Table (4.3) 1, below, the figures at years 1980 and 1990 represent percentages of land area, while the deforestation rate represents the total area lost in ten years as a percentage of the forest cover in 1980.
These results refer to the land area from which the SUs were selected, which excluded satellite scenes with less than one million hectares of land or estimated, according to the project's pan-tropical vegetation map, to have less than ten percent forest cover (see Section 2.1.1). The area surveyed represents 62 percent of the entire tropical land area (see Annex 4), comprising some 87 percent of total tropical forests, according to most recent FORIS estimates. As a consequence, the forest cover estimates presented here cannot be directly compared to other estimates which refer to total and regional tropical countries' land area1. In this context, for instance, to produce a comparable global estimate, the estimated “(F2) forest cover 1990” mean of 51.73 percent, given below, would correspond to some 32 percent when referred to the total land area, and should be increased to some 36 percent on account of the 13 percent of forest excluded from the area surveyed.
Percentage forest cover at year 1990 | Standard error of 1990 forest cover | Percentage forest cover at year 1980 | Standard error of 1980 forest cover | Rate of deforestation | Standard error of rate of deforestation | |
F-1: Definition of FOREST: Closed Forest | ||||||
AFRICA | 24.65 | ±3.39 | 26.37 | ±3.63 | 6.52 | ±1.25 |
L. AMERICA | 59.10 | ±2.89 | 63.32 | ±2.41 | 6.67 | ±1.48 |
ASIA | 35.08 | ±2.23 | 39.44 | ±2.53 | 11.05 | ±2.24 |
GLOBAL | 40.58 | ±1.84 | 43.83 | ±1.81 | 7.42 | ±1.02 |
F-2: Definition of FOREST: Closed Forest + Open Forest + part Fragmented Forest | ||||||
AFRICA | 43.42 | ±3.18 | 45.82 | ±3.29 | 5.23 | ±0.70 |
L. AMERICA | 65.44 | ±3.14 | 69.96 | ±2.63 | 6.46 | ±1.39 |
ASIA | 40.70 | ±2.12 | 45.19 | ±2.27 | 9.94 | ±1.84 |
GLOBAL | 51.73 | ±1.84 | 55.40 | ±1.75 | 6.62 | ±0.82 |
F-3: Definition of FOREST: Closed + Open Forest + Fragmented Forest + Long Fallow | ||||||
AFRICA | 45.44 | ±3.22 | 47.72 | ±3.32 | 4.77 | ±0.61 |
L.AMERICA | 66.83 | ±3.12 | 71.17 | ±2.64 | 6.10 | ±1.36 |
ASIA | 47.48 | ±3.16 | 51.72 | ±3.19 | 8.20 | ±1.44 |
GLOBAL | 54.45 | ±1.90 | 57.94 | ±1.81 | 6.03 | ±0.75 |
More than the mean value alone, one should consider the range of values within confidence limits, which are determined by the sampling error associated to such mean.
The statistical error associated with the global level estimates of forest cover at year 1990, ranging between ±3.5 and 4.5 percent of the mean, is quite small and very close to the ±3.9 percent value estimated during the design phase of the survey. A higher margin of error is associated with mean deforestation rates, ranging between ±12.3 and ±13.7 percent of the mean rates. This is not too high, however, considering the high variance deriving from the ‘event’ character of deforestation.
Future sample rounds will progressively reduce the sampling error associated to forest cover and deforestation estimates, thus allowing the production of results at more detailed reporting levels.
The graph in Figure (4.3) 1 illustrates the relationship between forest definitions and deforestation rates, highlighting their inverse correlation.
FOREST DEFINITIONS | ||
Fc 19901 | ||
F-1 = Closed Forest | Def + se2 | |
F-2 = Closed + Open + 2/3 Fragmented Forest | Def Rate3 | |
F-3 = Closed + Open + Fragmented Forest + Long Fallow | Def - se4 |
1 Estimated forest cover at year 1990, referred to Y axis on the left.
2 Deforestation rate + standard error, referred to Y axis on the right.
3 Mean deforestation rate, referred to Y axis on the right.
4 Deforestation rate - standard error, referred to Y axis on the right.
There is an apparent decrease in deforestation rate with the broadening of forest definition. This is determined by the fact that the changes within the forest, when a broad forest concept is used, are classified as degradation and fragmentation and not as deforestation. The forest definitions applied in this survey and the definition of forest area changes (deforestation, degradation, fragmentation, amelioration, afforestation, etc.) are discussed in Section 2.1.1.
How reliable are the results of the present survey?
This question cannot be answered by a single statement. There are elements with high reliability and others with low reliability. There are errors that can be quantified and errors that cannot be quantified. It is important, therefore, to discuss all these aspects and to try to understand both the strengths and weaknesses of the survey. This would help not only the user of the data presented here but also those who intend to carry out similar studies in the future. A thorough understanding of the error budget is in fact an essential prerequisite for the development of survey procedures.
The features of the survey that determine the errors, or that can be considered a source of error, can be defined as follows:
The first feature is a source of sampling error while the other two are sources of measurement error.
The sampling errors related to the first feature, viz., the statistical sampling design, are the easiest to estimate and have been reported along with results of forest cover mean estimates and deforestation rates in Section 4.3.
Depending on the definition of forest, the standard errors of forest cover estimates at pan-tropical level range between ±3.5 and ±4.5 percent. It is important to highlight that the result for the intermediate definition of forest, which correspond more closely to the preliminary data set used during the design phase of the survey, has a standard error of ±3.6 percent, which exceeds in precision the originally defined goal of ± 3.9 percent. This shows that both the survey procedure and the sampling intensity were chosen correctly.
The standard error of the estimated rate of deforestation for the same definition of forest is ± 12.4 percent, much higher than that of forest cover. This can be expected, however, as deforestation has a very high coefficient of variation.
An additional source of statistical error relates to the use of dot grids for area measurement (see Section 2.4). The error in area estimates caused by dot grid measurement cannot be quantified exactly but there are empirical formulae which have been developed to obtain a good approximation of such error. Error estimation can be achieved by applying the following formula (Chevrou, 1971), reported in the Manual of Forest Inventory (Lanly, 1973), and modified slightly by Baltaxe in 1991:
e % = 52 (┘k / n3/4)
where: | |
e % | is the percentage error of the measured area at the 95 percent probability level; |
n | is the number of dots counted within the area measured; and |
k | is a factor which depends upon the shape of the area and the regularity of its boundary; its value is between 5 and 7 for regular shapes and increases with the irregularity of the boundary. |
As the land area studied in the present survey is quite large (324.05 million hectares, corresponding to 810 112 dots), and all class measurements have been based on many dots, errors can be considered negligible even when choosing a high k value of 8, as shown in Table (4.4.1) 1:
Closed Forest | Open Forest | Long Fallow | Fragmented Forest | Shrubs | Short Fallow | Other Land Cover | Plantations | Total | |
Dot counted | 262 156 | 60 296 | 12 454 | 36 187 | 32 775 | 30 714 | 188 602 | 4 308 | 627 492 |
Estimated e % | 0.01 | 0.04 | 0.12 | 0.06 | 0.06 | 0.06 | 0.02 | 0.28 | 0.007 |
The measurement error related to the second feature, viz., change assessment, can be of paramount importance. In a comparison of two independent observations, the error associated with the change is higher than the error associated with the two original observations, as the law of propagation of error postulates. In the case of changes of small magnitude the resulting error could be too high to be acceptable.
In consideration of this problem the methodology developed for the assessment of class transitions was based on an interdependent interpretation procedure (see Section 2.3). This procedure ensured the highest level of consistency between the classification of recent and historical images. As a result of this approach the measurement error associated with class transitions can be minimized and remains comparable to that of the recent and historical classifications, as shown in the formulation below:
Errors in the estimation of change
INDEPENDENT versus INTERDEPENDENT interpretation approach
X1 = Interpretation at date 1 | ε1 = noise of image (date 1) |
X2 = Interpretation at date 2 | ε2 = noise of image (date 2) |
Y1 = true class at date 1 | μ1 = error of interpretation (date 1) |
Y2 = true class at date 2 | μ2 = error of interpretation (date 2) |
(ε1 + μ1) and (ε2 + μ2) are the errors ‘e’.
ε and μ may be positively correlated; i.e., much noise in the image makes the interpretation more uncertain.
Therefore: | Var (ε1 + μ1)= σε2 + σμ2 + 2Kov(ε1, μ1) = σe2 |
Assumption: independence between images. |
Independent interpretation | |
X1 = Y1 + ε1 + μ1 | Interpretation at date 1 |
X2 = Y2 + ε2 + μ2 | Interpretation at date 2 |
X1 - X2 = Y1 - Y2 + ε1 - ε2 + μ1 - μ2 | Change estimation |
VAR (state estimation) = σε2 + σμ2 + 2Kov (ε, μ)= σe2 | |
VAR (change estimation) = 2σε2 + 2σμ2 + 4Kov (ε, μ)= 2σe2 |
Interdependent interpretation | |
X1 = Y1 + ε1 + (μ1 + μ2)/2 | Interpretation at date 1 |
X2 = Y2 + ε2 + (μ1 + μ2)/2 | Interpretation at date 2 |
VAR (state estimation where there is no change; i.e., Y1=Y2)
classification according to: X = (X1+X2)/2 = Y1 + (μ1 + μ2)/2 + (ε1+ε2)/2
Var (X)= σε2/2 + σμ2/2 + (2Kov (ε, μ))/2 = σe/2/2
VAR (state estimation where the class has changed)
classification according to: X2 = Y2 + ε2 + (μ1+μ2)/2
Var = σε2 + σμ2/2 + Kov (ε2, μ2) = σe2/2 + σε2/2
VAR (change estimation) = Var (ε2 - ε1) = 2σε2
From the above formulation it would appear that the interdependent approach, where both images are used throughout the interpretation process, is beneficial not only in terms of change assessment but also in state assessment as the consultation of both images provides additional evidence in support of the delineation of the classes in both historical and recent images, as discussed in detail in Section 2.3.
The measurement error related to the interpretation accuracy of all 117 sampling units is very difficult to quantify. In fact, the accuracy of each interpretation could have been estimated exclusively on the basis of “true” reference (i.e. ground truth) collected in each study area at a very high cost; moreover, it is virtually impossible to estimate the interpretation accuracy of historical images.
Nevertheless, it is possible and important to associate the results with a critical evaluation of their reliability based on experience and on “educated guesses”. Following an approach similar to the fuzzy set theory (Zadeh, 1965; Gopal and Woodcock, 1994) it is possible to rank the results at various levels of reliability on the basis of some specific assumptions. Given an adequate knowledge of the study area, a good and coherent perception of the classification used and an acceptable image quality, it can be assumed that:
Following the above assumptions an estimated (indicative) biomass value has been allocated to each land cover class for each of the three broad ecological zones defined in the previous chapter. This value represents aboveground woody vegetation biomass in tonnes per hectare. These values have been defined on the basis of the best available information on average biomass values of closed and open forest classes by ecological zone; the values for the remaining classes in each zone were derived from these two “biomass indicators”, as explained in Section 2.8.
The Tables in Annex 8 give, in the form of transition matrices for each ecological zone, the estimated biomass difference (in absolute value) for all possible class-to-class combinations. This difference is here termed as biomass gradient. The values thus determined can be used in two different contexts:
Table (4.4.3) 1 provides a reliability ranking for the land cover classes in three ecological zones by stratifying the biomass gradients into three groups on the basis of their magnitude. Arbitrarily, the class relations (in space and in time) with a biomass gradient above 100 tonnes have been ranked as 1 (high reliability), those with biomass gradient between 50 and 100 tonnes as 2 (medium reliability) and those with a biomass gradient below 50 tonnes as 3 (low reliability). In view of its distinct spectral signature the class water and all combinations involving it are considered of high reliability.
The estimation of reliability applies equally to the spatial aspect (state assessment) and to the temporal aspect (change assessment); in the former, it represents the reliability in class delineation as a function of the contrast among adjoining classes, while in the latter it represents the reliability of class transition as a function of the contrast between the land cover class at date one and the class at date two.
Table (4.4.3) 1: Reliability ranking | |
1 = High reliability | (Biomass gradient >= 100 tonnes/ha) |
2 = Medium reliability | (Biomass gradient 50 – 100 tonnes/ha) |
3 = Low reliability | (Biomass gradient <50 tonnes/ha) |
Ecological zone: WET and VERY MOIST
Land Cover Categories | Continuous Natural Forest | Fragmented Forest | Other Wooded | Non-Wooded | Plantation | |||||
Land Cover Classes | Closed Forest | Open Forest | Long Fallow | Shrubs | Short Fallow | Other Land Cover | Water | |||
Continuous Natural Forest | Closed Forest | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | |
Open Forest | 1 | 3 | 2 | 2 | 2 | 2 | 1 | 2 | ||
Long Fallow | 1 | 3 | 3 | 2 | 2 | 1 | 1 | 3 | ||
Fragmented Forest | 1 | 2 | 3 | 3 | 2 | 2 | 1 | 2 | ||
Other Wooded | Shrubs | 1 | 2 | 2 | 3 | 3 | 3 | 1 | 1 | |
Short Fallow | 1 | 2 | 2 | 2 | 3 | 3 | 1 | 1 | ||
Non-Wooded | Other Land Cover | 1 | 2 | 1 | 2 | 3 | 3 | 1 | 1 | |
Water | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | ||
Plantation | 1 | 2 | 3 | 2 | 1 | 1 | 1 | 1 |
Ecological zone: MOIST with short and long dry season
Land Cover Categories | Continuous Natural Forest | Fragmented Forest | Other Wooded | Non-Wooded | Plantation | |||||
Land Cover Classes | Closed Forest | Open Forest | Long Fallow | Shrubs | Short Fallow | Other Land Cover | Water | |||
Continuous Natural Forest | Closed Forest | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | |
Open Forest | 1 | 3 | 2 | 2 | 2 | 2 | 1 | 3 | ||
Long Fallow | 1 | 3 | 3 | 3 | 2 | 2 | 1 | 3 | ||
Fragmented Forest | 1 | 2 | 3 | 3 | 3 | 2 | 1 | 3 | ||
Other Wooded | Shrubs | 1 | 2 | 3 | 3 | 3 | 3 | 1 | 2 | |
Short Fallow | 1 | 2 | 2 | 3 | 3 | 3 | 1 | 2 | ||
Non-Wooded | Other Land Cover | 1 | 2 | 2 | 2 | 3 | 3 | 1 | 2 | |
Water | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | ||
Plantation | 1 | 3 | 3 | 3 | 2 | 2 | 2 | 1 |
Ecological zone: SUB-DRY to VERY DRY
Land Cover Categories | Continuous Natural Forest | Fragmented Forest | Other Wooded | Non-Wooded | Plantation | |||||
Land Cover Classes | Closed Forest | Open Forest | Long Fallow | Shrubs | Short Fallow | Other Land Cover | Water | |||
Continuous Natural Forest | Closed Forest | 2 | 2 | 2 | 1 | 1 | 1 | 1 | 3 | |
Open Forest | 2 | 3 | 3 | 3 | 2 | 2 | 1 | 3 | ||
Long Fallow | 2 | 3 | 3 | 3 | 3 | 3 | 1 | 2 | ||
Fragmented Forest | 2 | 3 | 3 | 3 | 3 | 3 | 1 | 2 | ||
Other Wooded | Shrubs | 1 | 3 | 3 | 3 | 3 | 3 | 1 | 2 | |
Short Fallow | 1 | 2 | 3 | 3 | 3 | 3 | 1 | 2 | ||
Non-Wooded | Other Land Cover | 1 | 2 | 3 | 3 | 3 | 3 | 1 | 2 | |
Water | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | ||
Plantation | 3 | 3 | 2 | 2 | 2 | 2 | 2 | 1 |
Frequency of class transition (Percentage of total change in each ecological zone)
No box = less than 1 percent | 1 to 5 percent | 5 to 10 percent | Over 10 percent |
However, the reliability of classification (spatial aspect) is more difficult to relate to the final results discussed here since it depends mainly on the spatial relation existing among classes, which can only be seen on the spatial outputs of each individual SU. On the contrary, this estimation of class transitions reliability relates perfectly to the results since the change matrices report all class-to-class transitions originally observed in each SU.
In order to concentrate attention on the most frequent changes, Table (4.4.3) 1 indicates the class transitions that account for less than one percent of the total change observed in that particular zone, those representing between one and five percent, those representing between five and ten percent and the most important ones, representing more than ten percent of the total change.
Some comments:
Figure (4.4.3) 1 shows the indicative reliability of class transitions, of the pan-tropical change matrix, in relation to the number of changes observed. The pan-tropical reliability ranking has been estimated as the average of the three ecological zones weighted on the transition areas observed.
Figure (4.4.3) 1: Estimated reliability of pan-tropical class transitions |
Although not an estimation of accuracy, this figure can serve as a useful guide when using the results presented or considering specific class transitions.