Aerial Biomass Estimation in the Cerrado Biome Using Canopy Height Data

1. Introduction

Savannas cover one fifth of the global land surface, contributing one third of total terrestrial net primary productivity, and are responsible for more than half of global carbon emissions caused by wildfires [1]. The savannas of South America, found in the Cerrado domain, present notable carbon sinks and remain threatened by expansion of agricultural activities and by changing fire regimes [2]. Although the Cerrado biome covers around 23% of the Brazilian territory, only 2.5% of this area is protected [3]. The native Cerrado vegetation has been severely affected by degradation and deforestation and has consequently been fragmented since the 1970s [4]. The conversion of land use to agriculture in this biome is occurring significantly in vast native areas with great potential for terrestrial and soil carbon storage [5]. Despite the rich and complex biodiversity, there are still few studies to estimate the above- and below-ground biomass of Cerrado vegetation [6]. High rates of deforestation and changes in land use place the biome as the second largest Brazilian source of greenhouse gas (GHG) emissions [7]. Due to the heterogeneity of vegetation structure in the Cerrado, the stocks and distribution of aerial biomass are highly variable, making mapping and monitoring efforts non-trivial [8].

Obtaining accurate estimates of the carbon stock in the different phyto-physiognomies that compose the Cerrado is extremely important to support sustainable forest management plans [9]. Modeling growth and biomass production after the implementation of management systems in Cerrado areas is decisive for understanding vegetation succession process, especially regarding changes in biomass accumulation over time [10,11]. Data from multisensory remote sensing platforms combined with data from forest inventories, especially tree height, are efficient technological resources to quantify aerial biomass on a large scale [12]. The use of laser technology combined with remote sensing allows rapid mapping of biomass structure and density, relating forest fire regimes and vegetation dynamics, which is considered beneficial to achieve carbon emissions mitigation strategies in savannah ecosystems [8,13,14].

The predominant vegetation type in the Cerrado biome is Cerrado sensu stricto, which is identified as a savanna forest with 10% to 60% canopy cover and average tree height of 3 to 7 m [15]. In the central regions of Brazil, the distribution of plant biomass in different phytophysiognomies varies according to the structure of the formations, from shrubland, savannas and dense forests [16]. Research carried out in the field to quantify plant biomass in transition zones between biomes is rare, but essential for understanding carbon stocks in the different types of vegetation [17]. It is estimated that the decline in surface litter biomass in tropical forests is being affected by climate change and that the relationship between biomass strata and the dynamics of herbaceous, shrubby and arboreal vegetation remains without strong patterns of comparisons between phytophysiognomies [18]. The resilience of different types of tropical dry forests, commonly called Matas Secas (Dry woods), within the Cerrado biome and in transition zones with the Caatinga, is important for the maintenance and provision of ecosystem services [19]. Therefore, it is crucial to understand the ecological interactions between climate and biomass strata in the different phytophysiognomies of the Cerrado for the development of environmental management strategies in the context of climate change.

Thus, the present study assumes the existence of a linear relationship between tree height and forest biomass stocks in the different phytophysiognomies of the Cerrado biome, as has as its main objective to evaluate methodologies for estimation of total forest biomass. Analyses were performed using two independent sampling data sets, based on past measurements and a recent systematic review of aboveground biomass with sample plot values. The estimate of aboveground biomass was tested using an allometric equation developed from canopy height obtained via digital model and compared with the sample plots. The other biomass components (roots, litter, and necromass) were estimated to add up to the total biomass using aboveground expansion factors obtained from different sources.

2. Materials and Methods

2.1. Methodology flowchart

This study’s methodology involved a systematic review of biomass values reported for the Cerrado biome, as well as the comparison of variance in indirect estimates of aboveground biomass using linear allometric equation with canopy height through a digital model, Figure 1. For the estimation of aboveground biomass, the allometric equation proposed by [12] was implemented on the Google Earth Engine platform as a function of canopy height, which was obtained from the results provided by [20]. The pixel values of canopy heights were converted into aboveground biomass values, including necromass in standing trees. Canopy height values defined vertical strata for classifications of shrubland, savanna, and forest formations. To validate the estimated biomass values, two sample sets were tested: the first provided by [13], and the second obtained by composing the mean values of each inventoried fragment found in the systematic review papers, Table 1. The sample sets were evaluated through analysis of variance using the F-test at 5% probability, to identify if there are significant differences between the sample sets and the estimated values. Proportions of biomass stock compartments relative to total aboveground biomass were estimated based on information provided by the 4th Communication Report of the National Greenhouse Gas Inventory [7], grouped into expansion factors for shrubland, savanna, and forest formations. Total biomass for the entire Cerrado biome in Brazil was then estimated by summing the compartments of aboveground biomass (obtained from the canopy heights and regression model), necromass, litter, and belowground biomass of roots (obtained from the expansion factors), allowing the derivation of relationships between above- and belowground parts and between living and dead components.

2.2. Systematic review of forest biomass studies in the Cerrado

The methodology used in this research was based on the systematic review recommendations proposed by [21]. The systematic survey utilized document resources in scientific journals that reported aboveground forest biomass values in the Cerrado, aiming to identify methodologies related to different phytophysiognomies within the biome. The systematic search for scientific article reviews was conducted in June 2023, using the CAPES Journals database platform with the “Advanced Search” feature, entering the terms: “forest AND aboveground AND biomass AND cerrado”. All scientific articles with biomass estimation values from 1998 to 2023 were analyzed, including all publications available in the collections up to the date of June 2023. It is worth noting that the scope of forest biomass research in the Cerrado is broad, with a vast dataset available, and a systematic analysis is recommended to summarize the values found, possibly revealing relevant yet underexplored relationships [22].

The inclusion and exclusion criteria for the results of the articles were analyzed based on the potential use of presented biomass values, involving both live and dead aboveground biomass, with either field sampling or direct biomass estimation in the Cerrado biome. Only estimates from natural formations were considered, excluding plantations or exotic species. When available, the following information was compiled from the articles: location, plot area, sampling date, phytophysiognomy type, measurement method, tree diameter and height values, number of individuals per hectare, employed hypsometric, volumetric and allometric equations, wood basic density of dominant species, basal area, mean and standard deviation values of biomass (aboveground, roots, litter, and total).

The values of aboveground forest biomass were assessed according to the proportions of live and dead forest biomass, both above and below the ground. A total of 67 documents were reviewed based on the search criteria used, and 16 were selected for the analyses of the present study (Table 1).

2.3. Canopy Height values for the Cerrado biome

The model proposed by [20] is the first global canopy height model, estimating height values for the year 2020 as the reference date. This study employed scanning techniques with orbital laser data from the GEDI program to obtain surface models fused on terrain models and vegetation classifications from the Sentinel-2 imaging program. It developed a probabilistic model for estimation of the height of tree top canopies, including associated errors for estimation of canopy height values. The data is available in 10 m spatial resolution pixels in an open-source remote sensing program. This data can be accessed by anyone anywhere with internet access. Thus, it serves as a tool to support biomass and carbon monitoring for all types of forests and was used in this study to support biomass estimation.

Tree height measurements employ optical devices and trigonometric relationships, with models typically adjusted to predict height in relation to diameter, especially in silvicultural settings. However, in native areas, correlations are lower and do not yield as many significant gains in estimates. Ground-based LIDAR equipment can have significant potential to improve height estimates where tree structure allows laser penetration, thereby generating precise estimates of height profiles [23].

2.4. Characterization of the different Cerrado phytophysiognomies.

The Brazilian Cerrado occurs in altitudes ranging from 300 to 1600 meters above sea level, with annual rainfall between 600 and 2000 mm, and well-defined dry and wet seasons. It includes forest formations from Riparian and Gallery Forests, Dry Forest, Evergreen Forest, Semi-Deciduous Forest, Deciduous Forest, to Cerradão. Meanwhile, savanna vegetation formations are categorized into Dense Cerrado, Typical Cerrado and Sparse Cerrado. The shrub and arboreal structures have median strata of approximately 2 meters, and are referred to as “campos sujos” [15]. The classification of phytophysiognomies in the many studies assessed employs distinct definitions, ranging from regional approaches with definitions of open to closed forest structures and canopy cover, to global approaches resembling physiognomic-ecological definitions such as those of dry and wet seasonal forests, in altitude strata, as used in Brazil by [24].

The term “forest” adopted in Brazil considers vegetative cover with an area greater than 0.5 ha, trees of height exceeding 5 m, and canopy cover greater than 10%, or trees capable of reaching these parameters in situ. In the present study, height strata were determined according to the classification of [15], where trees with average canopy height of 1 to 3 m were designated as “campo”; savannas were stratified from 4 to 7 m; and forested areas with trees ≥ 8 m of average canopy height were labeled as arboreal canopy forest. These groupings were used for comparisons between biomass estimates (Table 2).

The definitions of steppe formations or variations with Steppe Savannas are not utilized in the classification by [15], these are associated with seasonal characteristics featuring herbaceous, shrubby, and tree vegetation typical of arid areas. These definitions may be confused with other regional designations such as Carrasco, Capoeira, or Caatinga. Global vegetation cover classifications simplify the grouping into forest and non-forest (shrubland) formations. Such simplification is employed in estimating carbon emissions and removals from plant biomass, as used in the IPCC Good Practice Guidance [25]. The canopy height grouping methodology in this study proposes to identify the forest component occurring in shrubland formations, typically overlooked in global mappings as non-forest.

2.5. Implementation of the allometric model of aerial biomass

The equation used to estimate aerial biomass selected in the systematic review was developed by [12], who related data from forest inventories obtained in the field with data on canopy heights using LIDAR made by high precision aerial mapping. The minimum tree diameter inclusion was DBH ≥ 5 cm, excluding leaf and thin branch biomass. To estimate biomass from the forest inventory data (tree diameter and height), the mixed species biomass equation proposed by [26] for each phytophysiognomy sampled was used. Height data sampled by LIDAR ranged from 1.7 to 38 m in vegetation transitions between savannah and forest, with aerial biomass between 19 and 104 Mg.ha⁻¹. The equation used to estimate the aerial biomass of living and dead trees and shrubs standing above the ground resulted in the linear equation (1), with an adjusted coefficient of determination (R²) of 0.89, root mean square error (RMSE) of 8.12 Mg.ha⁻¹ and bias of 0.43 Mg.ha⁻¹. Where: AGB is the standing live and dead aerial biomass and CHM is the canopy height model value.

AGB = 2.44 + 6.25 *CHM

(1)

In this study, the biomass estimation was obtained from Equation (1). Canopy height input values were obtained using the JavaScript programming language for canopy height mapping, as developed by [20] and accessible through the Google Earth Engine geoprocessing platform, for the ecoregions of the Cerrado biome. Total biomass results were generated by multiplying the aboveground biomass by expansion factors for litter and underground root biomass, in accordance with the groupings of shrubland, savanna, and forest formations stratified by canopy height (Table 2).

2.6. Validation of estimated aboveground biomass values

To assess aboveground biomass estimates, two sample sets were evaluated separately, and analysis of variance was conducted using the F-test at 5% probability to test the hypothesis of variability between field inventory values and estimated values from the digital model based on canopy height. The first sample set tested used forest inventory data from 77 sites, including 893 plots and 95,484 measured trees at different times, as made available by [13]. High-resolution spatial images were used to interpret shrubland, savanna, and forest formations and visually assess the conditions of phytophysiognomies at the 77 sites during the period of September 2020 (date of the canopy height records used in this study).

The second sample set evaluated used the geographic location of statistical parameters related to inventoried populations from 29 forest fragments in different phytophysiognomies, measured at different times, and distributed across the 16 scientific articles found in the Cerrado systematic review (Table 1). Classifications of shrubland, savanna, and forest formations were determined according to the authors’ classifications in the articles. All conditions of the fragments were visually assessed through high-resolution images from June to September 2020.

The forest fragments inventoried in the articles were overlaid onto the canopy height model to obtain points corresponding to the average canopy height for each type of inventoried phytophysiognomy. Consequently, the mean canopy height values were used as input data in the allometric model to estimate aboveground biomass. The aboveground biomass values were grouped by author, phytophysiognomy, mean value, and standard deviation for comparison with the digital model results. Analysis of variance and F-test were used on the means of aboveground biomass in pairs of sample sets for each type of shrubland, savanna, and forest formation to test hypotheses of variability among the data.

2.7. Proportions of forest biomass in relation to aboveground biomass (AGB)

Table 3 presents dry biomass values in Mg.ha⁻¹ available in the 4th national emissions report with the main biomass compartments by type of phytophysiognomy, using the vegetation classification adopted by [7]. Averages were calculated for each biomass compartment according to each type of plant formation structure

The results are available in georeferenced images, so the aboveground biomass values and relationships with total, underground and litter biomass from this study can be checked with 10 m pixels for any georeferenced field sampling.

3. Results

3.1. Biomass systematic review

Among the 16 selected articles (Table 1), 12 scientific journals were used as sources, with a notable emphasis on the journal Forest Ecology and Management, containing 5 out of the 10 most relevant articles. These journals have a high international impact factor and involve 178 authors in the study theme, with only one review article and no articles with a single author. The most involved countries in the research of the addressed topic are Brazil, Germany, and the United States. When analyzing the publication history of articles from 1998 to 2023, there is a noticeable increase in the number of studies in the last 5 years, with an annual growth rate of 2.81%. The 5 most common keywords were: Aboveground Biomass, Carbon, Savanna, Cerrado, and Vegetation. The articles presented a broad geographical sampling coverage within the biome, contributing to the understanding of results on the variability of forest biomass in different phytophysiognomies. However, there are still regions with less sampling or phytophysiognomies without field sampling, with a lack of actual dendrometric values for the vegetation.

Sampling and data collection from inventories of forest fragments involved demarcating plot areas within the vegetation, systematically or randomly distributed, varying in size (400 to 1000 m²), typically with dimensions of 20 × 50 m. Subdivisions of the plots, often measuring 10 x 10 m, are used to sample herbaceous and shrubby components. Most studies conduct measurements in a single period or specific periods, and the locations of the individual plots are generally not disclosed. Some works included more than one type of sampled phytophysiognomy, such as [8,12,14]. Geographically, 37 forest fragments were identified and grouped into 7 different Cerrado phytophysiognomies, with 17 sampling studies in Typical Cerrado and 11 in Cerradão (Figure 2).

3.2. Estimation of aboveground biomass in the systematic review and the use of the canopy height model

The measurements conducted in the reviewed articles date from 1988 to 2022. Among the evaluated articles, 3 studies used sequentially measured plots [2,11,27], and only one article mentioned vegetation increment values over time [11], reporting an aboveground biomass increase of 0.93 Mg.ha⁻¹.year⁻¹ in savanna. The research presented by [5] included measurements both before and after interventions such as wildfires. The remaining articles provided biomass stock values considering areas without human interventions, such as cutting and burning, during the measurements.

Among the reviewed articles, three studies were identified with direct measurements using destructive tree samples to adjust parameters of volumetric equations, along with wood density determinations in laboratory. In the Federal District of Brazil, 114 trees from 8 species were used in 1990 [13]. In the state of Minas Gerais, 120 trees from the 18 most abundant species were selected in 2011 [4]. In Bahia, 60 trees from the 8 main species were measured in 2019 [9]. The remaining articles employed an indirect methodology for biomass estimation, using existing allometric equations or adapting equations to the sampled phytophysiognomies. The aboveground biomass values reported in the reviewed articles are presented in Table 4.

The 77 samples provided by [13] with different periods of measurements were distributed across the ecoregions of the Cerrado biome, and had land use and vegetation cover visually classified using satellite images in September 2020, to check if the area was still covered by vegetation or had undergone land-use change. The results excluded 12 locations that underwent changes, which were 2 dam floodings, 9 conversions to pasture, and 1 to agriculture. Aboveground biomass values in the locations studied by [13] ranged from 4.77 to 50.22 Mg.ha⁻¹, showing an average value of 24.54 Mg.ha⁻¹. Vegetation cover varied across almost all phytophysiognomies, from open shrubland to riparian forests. The aboveground biomass values estimated by the canopy height model developed in the study ranged from 14.94 to 121.19 Mg.ha⁻¹, with an average of 62.93 Mg.ha⁻¹, substantially higher for all shrubland, savanna, and forest formations. When comparing the analysis of variance in single factor for mean biomass values, the hypothesis of variability between the biomass values found in the inventory plots by [13] and the mean estimates of the model proposed by [12] can be accepted with significant effect between groups, F (1.110) = 2.53 at p<.05.

The comparison of aboveground biomass estimation by the canopy model in the review articles involved 29 locations of inventoried forest fragments verified through satellite images. Two shrubland fragments, 14 savanna fragments, and 13 forest fragments were sampled. The aboveground biomass values found in the articles ranged from 3.9 to 184.10 Mg.ha⁻¹, with an overall mean of 60.24 Mg.ha⁻¹, including phytophysiognomies distributed across shrubland, savanna, and forest formations. The aboveground biomass values estimated by the canopy height model ranged from 14.94 to 152 Mg.ha⁻¹, with an overall mean of 63.67 Mg.ha⁻¹, and tree heights varying from 2 to 24 m.

Figure 3 presents variations in means and standard deviations across shrubland, savanna, and forest formations between the mean values and respective standard deviations found in scientific articles and the values of aboveground biomass estimated by the canopy height model developed in this study. When comparing the analysis of variance in single factor for mean biomass values, the hypothesis of variability between the review scientific articles and the values estimated by the digital model of biomass values can be rejected without significant effect between groups, F (1.54) = 0.77 at p<.05, considering the tabulated value lower than the critical value. Therefore, this means that the values estimated by the canopy height model are statistically equal to the mean values found in scientific articles.

Figure 4 (a) represents the results of aboveground biomass values in Mg.ha⁻¹ from the 29 samples, comparing the values reported in the scientific article review with the values estimated by the canopy height model developed in this study. The two extreme values (shown in red in Figure 4 (a)) were inspected in the images and found to be in the Pantanal biome, close to the limit with the Cerrado biome; hence, they were considered outliers because they are outside the Cerrado biome and were excluded from the analyses related to the research by [2]. In Figure 4 (b) the distributions of residuals are depicted without these outliers.

3.4. Assessment of aerial biomass in different phytophysiognomies

Figure 5 presents pixel values in pairs of comparisons from aboveground biomass units Mg.ha⁻¹ within the three variations of phytophysiognomies (campo, savanna, and forest formations) in the Cerrado biome related to the experimental farm. The image pairs illustrate vegetation formations, with forests represented by Riparian Forest (a, b), savanna represented by Typical Cerrado (c, d), and Shrublands (e, f). The images on the right side (b, d, f) represent the sampled fragments from the systematic review of scientific articles. Fragments with the results of the estimated biomass values are overlaid on the images on the right side (b, d, f). In each image, arrows with coincident points in the two images can be observed. Each point represents the values of aboveground biomass estimates, with the left side showing averages of the fragments sampled in the articles, and the right side showing the pixel of the average canopy height of the fragment and its respective estimated aboveground biomass value using the canopy height allometric model developed in this study.

Figure 5 (a) illustrates the Mata de Galeria fragment, with data provided by [23], presenting average tree height values of 12.7 m, with a standard deviation of 1.6 m and maximum height of 33.0 m. The fragment’s aerial biomass was estimated as 149.0 Mg.ha⁻¹ and the standard deviation was 74.5 Mg.ha⁻¹. The coordinate points -47.9694 and -15.9519, referring to the average canopy height (b) presents a pixel value of 24 m and a standard deviation of 6 m, with an estimated aerial biomass of 152.44 Mg. ha⁻¹. The illustration of the Typical Cerrado (c, d), also sampled by [23] on the same area, represented the savanna fragment with average tree height values of 2.8 m, standard deviation of 1.4 m and maximum height of 10 m. The aerial biomass value was estimated as 21.7 Mg.ha⁻¹, with a standard deviation of 11.6 Mg.ha⁻¹. The point sampled by the average canopy height to estimate aerial biomass (d) with coordinates -47.8461, -15.8641, resulted in a canopy height of 6 m and a standard deviation of 3 m, with aerial biomass estimated at 39.9 Mg.ha⁻¹. The comparison with the results from the shrublands, carried out on the same farm by [5] resulted in values from the sampled fragments with an average vegetation height of 2.5 m and a deviation of 0.35 m, with a maximum height of 6 m, and an estimated aerial biomass value of 17.6 Mg.ha⁻¹. The estimation data at the coordinate point -47.9166, -15.9500 in this study was 4 m of average canopy height, with 4 m of standard deviation and aerial biomass of 27.4 Mg.ha⁻¹.

3.5. Assessment of total biomass

From the canopy height relationships, it was possible to understand the structure of plant formation in the shrublands, savannah and forests and infer the other proportions of biomass in the Cerrado biome. All the aboveground biomass values reposted and estimated up to this point contain living and dead aerial biomass, thus they already include the percentage of standing necromass. Therefore, to calculate only live aerial biomass, the necromass must be subtracted from the aerial biomass. Thus, to estimate the total biomass of the Cerrado biome, the aerial biomass with dead and living parts were used to expand to the biomass of roots and litter. Table 3 summarizes the groupings of expansion factors for estimates of total biomass using the main compartments of forest biomass in relation to aerial biomass for the underground parts (fine and thick roots), litter, necromass and the possible relationships between the underground part by the total aerial part.

Table 3. Expansion factors (proportion) of biomass compartments in relation to total aerial biomass by type of vegetation.

Structure	AGB (Mg.ha−1)	Roots (%)	litter (%)	Necromass (%)	BGB:Total AGB
Shrubland	4.4 ± 2.2	280% ± 94	93% ± 1.9	1% ± 0.01	2.36
Savanna	22.1 ± 6.2	60% ± 48	12% ± 0.1	13% ± 0.03	0.46
Florest	69.2 ± 25.1	27% ± 12	9% ± 0.1	9% ± 0.04	0.23

Table 3. Expansion factors (proportion) of biomass compartments in relation to total aerial biomass by type of vegetation.

Structure	AGB (Mg.ha−1)	Roots (%)	litter (%)	Necromass (%)	BGB:Total AGB
Shrubland	4.4 ± 2.2	280% ± 94	93% ± 1.9	1% ± 0.01	2.36
Savanna	22.1 ± 6.2	60% ± 48	12% ± 0.1	13% ± 0.03	0.46
Florest	69.2 ± 25.1	27% ± 12	9% ± 0.1	9% ± 0.04	0.23

Using the aboveground biomass estimated by the canopy height model, expansion factors were applied to each type of vegetation formation. Thus, expansion factors for shrublands were applied to canopy heights from 1 to 3 m, savanna factors to canopy heights from 4 to 7 m, and forest factors to canopy heights above 8 m. To calculate the total biomass value in a forest of 100 Mg.ha⁻¹, factors corresponding to 27 Mg.ha⁻¹ for roots and 9 Mg.ha⁻¹ for litter were applied, resulting in a total biomass of 136 Mg.ha⁻¹. It is considered that a part of the total biomass corresponds to 9 Mg.ha⁻¹ of necromass and 91 Mg.ha⁻¹ of live aboveground biomass. Figure 6 presents the map of total biomass for all regions of the Cerrado biome in 10 m spatial resolution pixels and can be accessed through the Google Earth Engine platform through the link in the Supplementary Materials.

The results of the systematic review on underground biomass estimation presented uneven sampling methodologies. Due to the complexity of measuring the compartment of thick roots beneath trees, especially in larger trees found in forests, there is significant variability in averages and a lack of research on roots in forest formations in the Cerrado. The results on sample plots in this study, using root biomass expansion factors multiplied by the aboveground biomass in shrubland formations, resulted in 12.3 Mg.ha⁻¹. In savanna formations, average root values of 13.3 Mg.ha⁻¹ were found, and in forest formations, the average was 18.63 Mg.ha⁻¹.

4. Discussion

The distribution of the samples from the field plots and the inventoried fragments in the scientific articles did not adequately represent the northern portion of the Cerrado biome, specifically in the geographical regions of Alto Parnaíba and to the south in Paraná Guimarães, leaving gaps in these important areas, which could be prioritized for future studies. The standardization of sampling with explicit geographical locations and spatial precision overlaid on pixels in high-resolution images could contribute to increasing the accuracy of digital models for biomass estimates.

The definitions of minimum diameters and the concepts of trees and shrubs showed variability in the application of volumetric and biomass equations, which may impact aerial biomass estimation. [5] considered as shrubs individuals with a basal diameter (BD, measured at 0.30 m above the ground), with minimum values of ≤ 5 cm and a total height of up to 2 meters. Trees, on the other hand, are defined by diameter at breast height (DBH, obtained at 1.3 m (above the ground) with values greater than or equal to 10 cm [5]. In general, the definition of trees and shrubs established DBH ≥ 5 cm and height > 2 meters [2]. Some studies use DBH > 1.6 cm to include shrubs [6,17], while others employ intermediate diameters and separate diameters into classes > 10 cm to improve adjustments and reduce errors associated with biomass estimation equation fittings [8,14].

Plots sampled in earlier dates did not maintain linear correlations over the decades of measurements between the vegetation structures estimated with the digital model in the year 2020. According to [11], the biomass increment in the typical Cerrado can reach 26.11 ± 0.15 Mg.ha⁻¹ after 28 years of clear-cutting. In this context, the fragments inventoried in 1993, grouped by [13], found significant differences between the values of the plots and the values from the digital biomass model. On the other hand, the samples found in the scientific articles showed homogeneity in aerial biomass estimates using the digital canopy height model.

The savanna formation obtained more samples than the forest formation followed by the shrubland formation. It was observed in the results that deviations in heights are more significant in canopies of lower height. For savanna and forest formations, the deviations in canopy heights were less significant. The results found by [8,18] and in aerial biomass mapping in the Cerrado also indicated greater uncertainties with overestimations of values for lower heights in shrubland formations. Field measurements for forest biomass estimation indicate the need for differentiation of allometric models between shrubland and savanna.

Among the biomass compartments evaluated in this study, the underground root compartment is the most variable. The study of [2] presented an average for savannas of 14.2 Mg.ha⁻¹, which is the closest value to the present study (12.3 Mg ha⁻¹). On the other hand, the study of [6] only estimated values for fine roots at 0.7 Mg.ha⁻¹, and in the study of [27], a contrasting value of 3.5 Mg.ha⁻¹ was reported, also for fine roots. For total root estimates, [5] estimated high values of 30.1 Mg.ha⁻¹ for shrubland formations and 52.9 Mg.ha⁻¹ for savannas. Whereas [4] found a value for total roots in savanna formations of 37.5 Mg.ha⁻¹. The results of litter biomass are proportionally less representative of the total biomass and showed significant variation across different regions of the biome and measurements in several seasons. Sampling of necromass, in general, was not addressed in the inventories of review articles. The sampling of standing and fallen dead trees still confounds the inclusion of this biomass portion among the compartments of surface and aerial parts inventoried in review articles.

5. Conclusions

The allometric equation proposed by [12] for the estimation of total aboveground woody biomass, adjusted with the model employing canopy height developed by [20], showed similarity to the results found in 16 scientific articles distributed across the Cerrado biome. The results of these estimates can be accessed freely in georeferenced images with a pixel spatial resolution of 10 m. The mapping of aboveground biomass and associated errors was conducted for all canopy heights in the Cerrado biome. The expansion factors grouped for each vegetation structure of shrubland, savanna, and forest were related to the other compartments for estimation of underground biomass, litter, and total biomass. They proved to be effective and yielded satisfactory results but require improvements in variability control. The next steps will involve evaluations and classifications of native vegetation remnants, distinguishing them from silvicultural activities, and will be assessed through overlays in conservation units or rural properties for more accurate estimates of carbon stocks in the Cerrado biome.