Influence of native and exotic tree plantations on biophysical indicators in the Brazilian Savanna1

1 Received: Oct. 02, 2020. Accepted: Feb. 02, 2021. Published: Mar. 22, 2021. DOI: 10.1590/1983-40632021v5165815. 2 Instituto Nacional de Pesquisas Espaciais, Divisão de Observação da Terra e Geoinformática, São José dos Campos, SP, Brasil. E-mail/ORCID: bruno.adorno@inpe.br/0000-0003-0302-7834. 3 Universidade Federal de Goiás, Escola de Agronomia, Setor de Engenharia Florestal, Goiânia, GO, Brasil. E-mail/ORCID: sybelle.barreira@gmail.com/0000-0003-1482-2411. 4 Universidade Federal de Goiás, Instituto de Estudos Socioambientais, Goiânia, GO, Brasil. E-mail/ORCID: manuel@ufg.br/0000-0003-4516-6373. 5 Universidade Federal do Pará, Faculdade de Geografia, Altamira, PA, Brasil. E-mail/ORCID: gabrielveloso.geo@gmail.com/0000-0002-3655-4166. The monitoring of biophysical indicators can show the conservation or recovery status of a landscape. This study aimed to analyze the influence of tree plantations on the dynamics of biophysical indicators (albedo, NDVI, surface temperature and evapotranspiration), in an experimental area of the Cerrado biome (Brazilian Savanna), by applying remote sensing techniques and the SEBAL algorithm. The indicators dynamics were given as a function of changes in the land use, while assessing the response of the environment to the planting of tree species. SEBAL data on areas that underwent changes in land use and cover during this period were analyzed. In the surroundings of the experimental area, albedo and surface temperature decreased in agricultural and exposed soil areas converted to tree plantations, while the NDVI and evapotranspiration increased. The opposite happened in the conversion of native areas destined to agriculture and livestock. In the experimental area, it was confirmed that the plantations contributed not only to the decrease in the surface and albedo temperature, but also to the increase in the NDVI and evapotranspiration. This confirms the positive influence of tree planting in rural properties of the Cerrado as a support to environmental regularization and more sustainable agricultural systems. Moreover, it highlights the potential of the technique applied to assist in monitoring Cerrado ecosystems in areas larger than those commonly monitored in the field.


INTRODUCTION
More and more global targets have been established every year due to the increasing loss of ecosystem services (Sarukhán & Whyte 2005, Wood The monitoring of biophysical indicators can show the conservation or recovery status of a landscape. This study aimed to analyze the influence of tree plantations on the dynamics of biophysical indicators (albedo, NDVI, surface temperature and evapotranspiration), in an experimental area of the Cerrado biome (Brazilian Savanna), by applying remote sensing techniques and the SEBAL algorithm. The indicators dynamics were given as a function of changes in the land use, while assessing the response of the environment to the planting of tree species. SEBAL data on areas that underwent changes in land use and cover during this period were analyzed. In the surroundings of the experimental area, albedo and surface temperature decreased in agricultural and exposed soil areas converted to tree plantations, while the NDVI and evapotranspiration increased. The opposite happened in the conversion of native areas destined to agriculture and livestock. In the experimental area, it was confirmed that the plantations contributed not only to the decrease in the surface and albedo temperature, but also to the increase in the NDVI and evapotranspiration. This confirms the positive influence of tree planting in rural properties of the Cerrado as a support to environmental regularization and more sustainable agricultural systems. Moreover, it highlights the potential of the technique applied to assist in monitoring Cerrado ecosystems in areas larger than those commonly monitored in the field. KEYWORDS: Remote sensing, legal reserve, land use and cover change.
Notwithstanding, many of these areas are under intense conversion, with high deforestation records. This is the case of the Cerrado (Brazilian Savanna) (Brasil 2018). Although considered a biodiversity hotspot (CEPF 2018), this biome has been negatively impacted by the advancement of the agricultural frontier (Alencar et al. 2020). This advancement has initiated in the second half of the twentieth century by government programs such as the Nipo-Brazilian Cooperation Program for the Development of Cerrados (Prodecer) (Kazuhiro 2000). Additionally, water scarcity on hydrographic systems, such as the Araguaia-Tocantins and Meia Ponte, is more severe each year (Mascarenhas et al. 2009, Coe et al. 2011, Latrubesse et al. 2019. To mitigate these adversities, the union of government efforts, civil society and scientific community in decision making and territorial governance is urgent. In the political environment, the Federal Government has promoted management and supporting tools to the regularization of legal reserves and permanent preservation areas, and has encouraged more sustainable production systems, such as agroforestry systems (Brasil 2012). Recently, another widespread action ("Together for Araguaia" or, in Brazilian Portuguese, Juntos pelo Araguaia) has been dealing with forest restoration at the head of the Araguaia River to expand the water availability in areas with higher recharging potential. This action is supported by the Ministry of Regional Development, universities and research institutes (Ferreira Neto 2018).
To the date, the monitoring of ecological indicators for ecosystem restoration has been usually performed with field (in loco) surveys. This process includes high logistical effort and is limited in spatial terms (i.e., coverage of small sample areas). In this context, it becomes interesting and feasible to associate environmental analyses via remote sensing, especially in more extensive areas (Kumar et al. 2015). Sensors on board satellites or aircrafts make it possible to measure biophysical indicators such as albedo, normalized difference vegetation index (NDVI), surface temperature and actual evapotranspiration, which are influenced by the types of land use and cover (Veloso et al. 2017, Veloso et al. 2020. Albedo corresponds to the solar radiation reflection of a given target. Vegetation tends to present low albedo, since it favors the absorption of radiation for photosynthesis (Bala et al. 2007), which also means lower atmospheric heating. The NDVI, corresponding to the normalized ratio between reflectances in the near infrared band and red band, indicates more photosynthetically active vegetation areas with positive values (Ponzoni et al. 2012), showing their level of degradation or conservation. These and other indicators, detailed by Veloso (2014) and Veloso et al. (2020), can measure absorbed radiation. This radiation contributes to measure both the increased surface temperature, in the form of sensible heat, and the phenomena occurring in vegetation areas (such as evapotranspiration), in the form of latent heat (Veloso 2020). That is, these indicators are greatly sensitive to the supply of ecosystem services and to changes in the land use and cover, being, therefore, important in the monitoring of recovery processes in anthropized areas.
Methods that use surface meteorological data and remote sensing techniques to measure these indicators have been developed. Among them, stands out the estimation of radiation and energy balance by the Surface Energy Balance Algorithm for Land -SEBAL (Bastiaanssen 1995), using satellite images as an input parameter in the model to obtain the actual daily evapotranspiration.
The SEBAL method has some advantages over others (e.g., Simple Algorithm for Evapotranspiration Retrieving -SAFER and Mapping Evapotranspiration at High Resolution with Internal Calibration -METRIC; Silva et al. 2019a). Some examples are the greater accuracy of this algorithm (Menezes et al. 2011); greater number of scientific studies using it; better adjustment to the input parameters, with lower dependence on the use of terrestrial data for model calibration; and the possibility of generating information for each image pixel (Bastiaansen 2010 ).
In this perspective, this study monitored the dynamics of vegetation biophysical indicators (albedo, NDVI, surface temperature and evapotranspiration), in an experimental area of the Cerrado, by applying remote sensing techniques and the SEBAL algorithm.

MATERIAL AND METHODS
The study was conducted at the Entre Rios Farm (15º56'50.97"S and 47º28'32.25"W), located in the Federal District of Brazil ( Figure 1). This is an experimental area of the Biomas project, a partnership between the Brazilian Confederation of Agriculture and Livestock (CNA), Brazilian Agricultural Research Corporation (Embrapa) and other interested parties from different Brazilian biomes. The mission of the project is to present to the society, especially to rural producers, tree plantation models for economic purposes and environmental regularization (CNA 2009). Tree species were planted in restricted production areas of permanent preservation and legal reserve, in the last quarter of 2012, in accordance with the Law 12,651 of May 25, 2012 (Brasil 2012) and Law 12,805 of 2013 (Brasil 2013). The aim of this planting was to regularize the rural property and promote the recovery of degraded pastures.
According to the Köppen classification, the climate in the region is Savanna tropical (Aw), with concentrated rainfall from October to April (approximately 1,400 mm, 80 % of the annual total). The relative humidity may range from 20 to 70 %, and the average annual temperature is 22 ºC (Silva et al. 2008).
In August 2015, 12 areas where tree species were planted for various purposes were identified and mapped (Table 1). A Garmin satellite-navigation device (model Etrex) was used to collect control points, which supported the subsequent classification of satellite images. In the absence of a local meteorological station for the two analyzed periods, information from the Brazilian National Institute of Meteorology was used (Brasil 2020), referring to the available point closest to the study area, located at the Águas Emendadas Ecological Station, rural area of the Brazilian Federal District ( Table 2).
The Surface Energy Balance Algorithm for Land (SEBAL), described in Bastiaanssen (1995), was chosen to generate biophysical indicators for each target, allowing a quantitative comparison between them. The methods for its use in Landsat 5 and Landsat 8 images were detailed by Veloso (2014) and Veloso et al. (2020). The algorithm stands out for its greater accuracy in estimating energy balance and evapotranspiration (Menezes et al. 2011).
Still on the data used in SEBAL, Landsat 5 TM and Landsat 8 OLI/TIRS satellite images were obtained from the United States Geological Service (USGS 2020) website, corresponding to path 221, row 71, for November 4, 2008 and September 21, 2015, respectively. The choice of these images is due to the proximity of the date of data collection in the field. Those with cloud cover of up to 10 % were considered. Another criterion was the dry period, with extremes values of biophysical indicators. Additionally, to increase the chance of observing significant changes in the land use, the analysis considered the oldest image with meteorological data available at the Brazilian National Institute of Meteorology (October 02, 2008) (Brasil 2020).
The data processing was divided into two steps: 1) mapping of land use in the farm and its immediate surroundings (limited by the coordinates 15º50'2.11"S and 47º36'13.84"W in the upper left quadrant and 16º03'30.71"S and 47º20'50.22"W in the lower right quadrant) for the two images (2008 and 2015), using the study by Sano et al. (2009) as a reference for the interpretation of Cerrado classes; 2) analysis of the images processed by the SEBAL (limited by the coordinates 15º29'20.45"S and 47º43'49.75"W in the upper left quadrant and 16º09'41.94"S and 47º13'46.33"W in the lower right quadrant).
For mapping, the image pixels were segmented using the feature extraction tool (scale = 90; merge level = 10) of the ENVI software. Polygons were classified by visual interpretation for both periods. Figure 2 shows the aspect of land use and native vegetation observed for visual interpretation within the study site, taking the image of September 2015 as an example.
For comparison purposes, areas outside the experiment (Entre Rios Farm) that underwent changes in the land use and native vegetation were selected and analyzed for the response of biophysical indicators. For that, the average of the indicators values, represented in the pixels delimited by the classes of land use and native vegetation, was compared using the Tukey test (α = 0.05). The same procedures were applied to the changes observed within the farm, with careful sample selection, observing the control points collected in the field. It is important to note that, in this area, the agriculture and livestock classes were grouped as "agro-livestock". This is because the crop-livestock integration started to be practiced as a productive activity in the farm between the analyzed periods.

RESULTS AND DISCUSSION
Based on the mapping of the study area carried out for the two dates (November 2008 and September 2015), the agriculture and livestock classes already predominated in the region for the first date, while the remaining native vegetation was present in only 22 % of the area. The changes that stood out the most in the landscape were selected for further statistical analysis of biophysical indicators by the SEBAL. Figure 3 shows the results of this mapping for the immediate surroundings of the farm, with delimitation of the landscapes that stood out the most, in terms of changes in land use and native vegetation, in the analyzed period. Sano et al. (2009) evaluated the mapping accuracy for different land covers in the Cerrado and achieved a 71 % accuracy, with confusion in the classification of native woody classes (forest and even shrub physiognomies) and classes of anthropic use, such as pasture and agriculture. However, native areas were well differentiated from anthropized areas with 90 % of accuracy.
Knowing the possibility of error to which the classification is subjected, and with few changes observed in the mapping of this study, the comparison of the mappings for the surroundings of the farm between the two satellite images (from 2008 and 2015) does not assure significant changes in the land use and native vegetation, as there was no in loco validation. On the other hand, possible changes on the farm (Figure 4) were confirmed from the control points.
There was a reduction of 2.3 % for the exposed soil class, and 2.4 % for agriculture and livestock, together. On the contrary, there was an increase of 1.5 % for pre-existing forest and savanna on the farm. Native and exotic tree planting areas totaled about 40 hectares (3 %) of the farm, where various planting arrangements are being experimented for economic and environmental purposes. These areas replaced those previously classified as agriculture, as noted in the Entre Rios Farm comparison map ( Figure 4).
Noteworthy, the reduction of exposed soil areas, observed between 2008 and 2015, occurred in   Table 3 presents the average of the results observed in this study for albedo, NDVI, surface temperature and evapotranspiration in the different fragments of land use and cover in the Entre Rios Farm and its surroundings, in the analyzed periods (2008 and 2015).
The results of this study showed a reduction in albedo values for burning areas (10.2 %), followed  by forest areas (planted or native) and savanna areas (13.7 to 15.3 %), and, finally, exposed soil areas, agricultural areas and pastures (17.1 to 20.4 %), respectively. Other authors also found similar rates: 8 to 13 % in burning areas (Lyons et al. 2008), these values being generally observed immediately after a more severe fire (Quintano et al. 2019); 11 to 15 % in forest areas; and 13 to 20 % in pastures (Querino et al. 2006, Giongo et al. 2009, Silva et al. 2015, Veloso et al. 2017; 13 to 26 % in agricultural areas; around 13 % in typical Cerrado vegetation (Giongo et al. 2009, Veloso et al. 2017; and 18 to 45 % in exposed soil areas (Oliveira et al. 2013, Veloso et al. 2017. As for NDVI, which can range from -1 (water) to 1 (maximum photosynthetic activity of vegetation), results are generally greater than 0.600 in dense forests (Bayma & Sano 2015, Martins et al. 2015; between 0.260 and 0.460 in grassland and savanna regions of Cerrado (Trentin et al. 2013, Bayma & Sano 2015; between 0.170 and 0.200 in pastures (Veloso et al. 2017); and between 0.250 and almost 0.800 in agricultural areas during the dry season, but greater than 0.550 when using irrigation (Trentin et al. 2013, Veloso et al. 2017. In this study, the highest NDVI values occurred for tree plantations (average of 0.629). It is noteworthy that, depending on the vegetation density, NDVI may either saturate or have its value decreased in shaded areas such as primary forests (Ponzoni et al. 2012). However, the results were satisfactory for the purpose of characterizing such different land uses and covers. Even so, other vegetation indicators also generated by SEBAL, such as the Soil-Adjusted Vegetation Index (SAVI), can be analyzed in future studies due to their good adjustment for different thematic classes (Silva et al. 2019b).
When studying eucalyptus plantations, Almeida et al. (2015) found NDVI values greater than 0.500 for crops older than one year. Agricultural, exposed soil and livestock areas showed NDVI between 0.220 and 0.266, while Cerrado areas showed average NDVI between 0.266 and 0.353. Evergreen or semideciduous forests, on the other hand, had an average NDVI of 0.413. Burning areas had the lowest NDVI values, what was already expected due to the direct impact on their loss of photosynthetic activity. Souza et al. (2015) calculated NDVI values and showed that the Cerrado typologies may take from 94 to 100 days (with standard deviations of 46 to 50 days, respectively) to recover its vigor after burning.
The highest surface temperature occurred in burning areas (39.6 ºC). Livestock and agricultural areas showed average values between 36.1 and 38.1 ºC, while exposed soil and some Cerrado typologies showed values between 33.6 and 36.2 ºC. In turn, native areas, including some of the denser Cerrado areas, had the lowest values (30.6 to 32 ºC).
Corroborating these values, studies such as that by Eltz & Rovedder (2005) also demonstrated that, between different soil covers, the more forested an area, the lower is the temperature close to the ground in that region. Cerrado typologies, and land uses such as agriculture and livestock, generally have higher values of surface temperature (Gusmão et al. 2013, Martins et al. 2015, Silva et al. 2015, Veloso et al. 2017. Santos et al. (2017) analyzed the Caatinga biome using the SEBAL and showed that the daily evapotranspiration is similar for agricultural areas, grasslands and the characteristic vegetation of that biome, but higher for exposed soil areas and lower for areas with more dense canopies. This pattern tends to be the opposite of that observed for surface temperature. For example, Ning et al. (2017) also adopted the SEBAL in their methodology and showed how vegetation can influence measures of surface temperature and daily evapotranspiration, which are highly correlated (R 2 > 0.95) and inversely proportional.
Evapotranspiration was higher in planted forest (> 4.29 mm day -1 ) and lower in agricultural and exposed soil areas (1.58 to 2.45 mm day -1 ). The other classes showed intermediate values. The results were consistent with the literature values for forest areas, which show evapotranspiration greater than 4 mm day -1 or less than 3 mm day -1 in certain more extreme periods, in semideciduous or deciduous forests (Veloso 2014). However, areas with more sparse native vegetation, pastures or non-irrigated agricultural areas hardly present values greater than 2.5 mm day -1 in the dry period (Veloso et al. 2017).
Finally, the same analysis was performed for the classes of land use and native vegetation mapped on the farm. Therefore, it is important to highlight that planted areas came out of a closer relationship of biophysical indicators with agriculture in 2008, assuming values closer to those of the Cerrado typologies, as observed in 2015, for all indicators. This demonstrates the environmental role of planting trees. Notwithstanding, future analyses for native and exotic trees need to be conducted to understand the behavior of indicators in each case, as they may present different growth rates and interaction with the biophysical environment. CONCLUSIONS 1. Remote sensing, in particular with the use of free satellite images and SEBAL modeling, enables the analysis of biophysical indicators characteristic of different types of land use and cover, and may be an auxiliary tool to monitor the conservation status of Cerrado ecosystems; 2. Planting tree species for the purpose of environmental regularization or sustainable production (with a focus on pasture recovery) increases the NDVI and evapotranspiration, while reducing the albedo and surface temperatures; 3. The analysis method that used SEBAL together with the processing and thematic classification of satellite images proved to be quite efficient, representing a more extensive area than it is normally evaluated in the field.