SciELO - Scientific Electronic Library Online

 
vol.23 número2Estimación de la relación de pérdida de suelo para el cultivo de colza (Brassica napus) a partir del modelo USLE/RUSLE en UruguayInfluencia de los precios fundiarios en la reproducción de los ganaderos familiares de Barriga Negra (Lavalleja, Uruguay) índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Journal

Artigo

Links relacionados

Compartilhar


Agrociencia (Uruguay)

versão impressa ISSN 1510-0839versão On-line ISSN 2301-1548

Agrociencia Uruguay vol.23 no.2 Montevideo dez. 2019  Epub 01-Dez-2019

https://doi.org/10.31285/agro.23.85 

Articles

Spatial Distribution and Tree Cover of Hillside and Ravine Forests in Uruguay: the Challenges of Mapping Patchy Ecosystems

Distribución espacial y cobertura arbórea del bosque serrano y de quebrada en Uruguay: los desafíos de mapear ecosistemas parchosos

1Grupo Biodiversidad y Ecología de la Conservación, Instituto de Ecología y Ciencias Ambientales, Facultad de Ciencias, Universidad de la República, Iguá 4225, Piso 8, Montevideo, C.P. 11400.

2Departamento de Producción Forestal y Tecnología de la Madera, Facultad de Agronomía, Universidad de la República, Av. Garzón 780, Montevideo, C.P. 12900.

3Polo de Ecología Fluvial, Departamento del Agua, CENUR Litoral Norte, sede Paysandú, Estación Experimental Dr. Mario A. Cassinoni (EEMAC), Universidad de la República, ruta 3 (km 363), C.P. 60000, Paysandú.

4 Polo de sistemas territoriales complejos, CENUR Noreste, sede Rivera, Universidad de la República, Ituizangó 667, C.P. 40000, Rivera. Email: ctoranza@gmail.com


Summary:

The mapping and monitoring of forest ecosystems on a national scale is key to their management and conservation. Native forests in Uruguay are considered given their importance for biodiversity and ecosystem services. Here we evaluate the spatial distribution of the land cover class ‘hillside and ravine forest’ -a subclass of native forest characterized by patches and transition zones with native grasslands- using Landsat images (30 x 30 m) from 2014 and 2015 and high-resolution images from Google Earth. To evaluate spatial heterogeneity within hillside forests, we then used high-resolution spot images of 1 km2 from 1998-2012 to evaluate differences in the normalized difference vegetation index (ndvi) among canopy coverage categories. The hillside and ravine forest class were characterized as a composite cover class with an average canopy coverage of 69 ± 23%, variability of wich was reflected in ndvi values. The total area of this class in 2015 was estimated as 334,480 ha, somewhat less than an earlier 2008 estimate (384,240 ha). Among the potential errors in delineating hillside forests using Landsat images, there was the classification of «forest» in areas characterized by grassland and a tree canopy cover <25 %. This potential error in delimitation at broader scales led to the overestimation of hillside and ravine forest area in southeastern Uruguay, but an underestimation in northern Uruguay. Our study highlights the large discrepancies in the estimation of the distribution of hillside and ravine forest at different spatial scales, and also indicates the potential of ndvi to evaluate the heterogeneity of this forest within the same cover class.

Keywords: subtropical forests; ndvi; canopy; southeastern South America; Pampean province

Resumen:

El mapeo y monitoreo de la distribución de los ecosistemas a escala nacional es clave para su manejo y conservación. Los bosques nativos de Uruguay son ecosistemas de interés nacional, de gran importancia para la biodiversidad y servicios ecosistémicos. Aquí evaluamos la distribución espacial del «bosque serrano y de quebrada», una formación parcheada y en transición con los pastizales serranos, utilizando imágenes de Landsat (30x30 m) de 2014 y 2015, e imágenes de alta resolución de Google Earth. Adicionalmente, se utilizaron imágenes spot (1 km2), del período 1998-2012, para evaluar diferencias en el Índice de Vegetación de Diferencia Normalizada (ndvi) entre distintas categorías de densidad arbórea. El bosque serrano y de quebrada se caracterizó por presentar una cobertura media de dosel arbóreo de 69 ± 23 %, cuya variabilidad se reflejó en los valores de ndvi de distintas categorías. La distribución espacial del bosque se estimó en 334.480 ha en 2015, siendo una superficie menor respecto a la cobertura institucional de 2008 (384.240 ha). Se detectaron áreas de pastizal serrano con una cobertura arbórea <25 %, clasificadas como bosque con base en las imágenes Landsat. Esto genera un error en la delimitación del bosque a escala gruesa, con una sobrestimación del área boscosa en el sureste y una subestimación en el norte del país. Este trabajo destaca la discrepancia entre las estimaciones del área total de bosque serrano y de quebrada a distintas escalas espaciales y destaca el potencial del ndvi para evaluar su heterogeneidad interna.

Palabras clave: bosques subtropicales; ndvi; dosel arbóreo; sureste de Sudamérica; Provincia Pampeana

Introduction

Forest monitoring systems at a country national scale are fundamental tools for land management and the mitigation of climate change impacts and land use intensification1. Information on spatial extent and productivity of forests feed global databases and are highly valuable when participating in international initiatives such as REDD+ or the Montreal Protocol2. In addition to their central role in reducing carbon emissions3, forests are one of the most biodiverse terrestrial ecosystems4 and provide multiple ecosystem services, including the control of floods and diseases, they provide timber and non-wood forest products, and offer a myriad of recreational activities5)(6. The quantification of ecosystem services, the evaluation of forest conservation status and the elaboration of management plans depend in part on the quality of data on the spatial distribution of forests, and their availability to the institutions responsible for forest management. Efforts to map forests have increased in recent decades7)(8)(9, mostly focusing on monitoring trends in forest distribution associated with global change10.

Estimating forest area via tree canopy coverage has become one of the main ecological indicators incorporated into forest inventories11)(12. Tree canopy cover is defined as the proportion of the ground covered by the vertical projection of tree crowns13, and it constitutes one of the main structural characteristics of forests, being fundamental for the maintenance of biodiversity and the continuity of ecosystem services14. Tree canopies constitute much of the habitat for flora and fauna associated with forests, they participate in the hydrological cycle through evapotranspiration and rainfall interception and have an important role in climate regulation15. Additionally, in savanna ecosystems, tree canopy cover is also a key player in cycling soil C and N, as well as in fire dynamics16)(17. Likewise, canopy density is strongly linked to carbon reserves in forests and savannas18)(19. Multiple tools have been developed to measure tree canopy cover in the field and, more recently, remote sensing with satellites imagery or LIDAR has become widespread as one of the main techniques for estimating canopy cover at broader scales20)(21.

The need to develop strategies for climate change adaptation and mitigation has promoted the development of monitoring systems to evaluate rates of primary productivity and net carbon storage at country level18. Among the most commonly used indicators for monitoring primary productivity, the Normalized Difference Vegetation Index (ndvi) is derived from red and near infrared bands from satellite imagery22)(23)(24. The utility of the ndvi as an indicator of primary productivity is limited in forests with high canopy density due to saturation25, but it is widely used for the monitoring of forests, savannas and grasslands in southern South America and African savannas26)(27)(28.

In Uruguay, native forests occupy approximately 4.5 % of the land surface area29)(30, a value that varies depending on the estimation method and the scale of analysis used31. According to our knowledge, there is only one study addressing the temporal changes in the spatial extent of native forests countrywide32. An upward trend in the area covered by native forest for the period 1980-2006 was recently reported for Uruguay33. Meanwhile, the General Agricultural Censuses indicate a decrease in forest cover for the period 2000-201134)(35. Other studies at a smaller scale have reported a decrease in native forest area as a result of anthropogenic activities, particularly those pertaining agriculture expansion and tree plantations36)(37. Along these lines, quantifying the recent change in native forest area is one of the objectives of the redd + Project in Uruguay, initiated in 2016.

In order to further contribute to the monitoring of Uruguayan native forests, the objectives of this study are: (1) to generate an updated coverage of the spatial distribution of the hillside and ravine forests in Uruguay (locally referred as monte serrano and monte de quebrada respectively, from now hrf); (2) to quantify the canopy cover of the hrf, and (3) to evaluate the relationship between canopy cover and ndvi in the hrf. To achieve these goals, a series of Geographic Information System (gis) tools and satellite images were used at different spatial resolutions to obtain data of the forests area, ndvi and the degree of coverage in the hrf of Uruguay.

Materials and Methods

Forest cover classes

To generate the hrf land cover class, henceforth referred to as the «hrf-2015 layer», we used three spatial datasets: 1) Landsat 8 images (resolution: 30x30 m), consisting on the Operational Land Imager (oli) and the Thermal Infrared Sensor (tirs), mainly from the austral-summer months (Dec-Feb) of 2014/2015, 2) high resolution images from the summer 201538, and 3) the official hillside forest layer of 2008 from the Uruguayan Land Cover Classification System31.

Landsat 8 images were downloaded from the United States Geological Survey portal39. A set of thirteen images, which together covered the entire land surface of Uruguay, were joined in a mosaic image, which was integrated into a raster and served as the basis for limiting the spatial extent of the analysis to the national boundaries of Uruguay. This national boundary was defined by a vector layer of «terrestrial boundaries» elaborated by the Uruguayan military geographical service. Radiometric correction of Landsat images was deemed unnecessary for several reasons. First, forest patches were strictly identified by visual interpretation of the Landsat images (cross-checking with Google Earth high-resolution images). Second, all Landsat images were obtained from the austral summer months (Dec-Mar) of 2014/2015, which correspond with the growing season for most woody species in the region. Finally, the hrf is dominated by evergreen trees40, that are photosynthetically active throughout most of the year. To reduce the influence of cloud cover, images with a cloud cover of less than 10% were selected for analysis.

Given the patchiness of forest ecosystems and their relatively small size, we used high-resolution images from Google Earth38 from the summer 2015 -viewed at a scale of 1:7000-, as well as Landsat images to delineate polygons based on the presence of forest patches within each polygon.

As a third source of spatial information, we also used an official Land Cover layer from 2008 as a guide for the delineation of hrf patches. This layer of native forest in 2008 was obtained from the land cover classification of Uruguay developed among multiple institutions in 200831, henceforth referred to as «lcc-2008». In the lcc-2008, natural and anthropogenic uses were classified based on Landsat 5 tm images (30x30 m) from 2007 and 2008 (scale 1:100,000). We isolated the hrf class, and superimposed the lcc-2008 layer over the 2014/2015 Landsat 8 satellite images in ArcMap 10.1, working at a 1:30,000 scale.

The 2015-hrf layer was delineated by hand using tools from the editor toolbar in ArcMap 10.1. In most cases the polygons of the lcc-2008 layer were removed and re-digitized to conform to the shape of forest patches in Landsat 8 images, checking for the presence of forest patches in Google Earth. Given the distinctive nature of the hrf patches (compared to the surrounding vegetation) they were easily identified. Although the reflectance of native forests was similar to those of tree plantations, they can be distinguished by their shape, since tree plantation have geometric spatial patterns. Once this process was completed, a corrected cover of 2015 was obtained, which included the following forest formations recognized in Uruguay: hillside, ravine, escarpment and rocky soils forest, in Spanish: serrano, quebrada, escarpas and mares de piedra41. This new layer was ground-truthed using the locations of plots categorized in the field as hillside and ravine forests (N = 152 plots) in the National Forest Inventory (ifn). The 2015-hrf layer was again revised and corrected to include those IFN plots excluded by the layer, thus reducing omission errors.

Differences in ndvi and canopy cover within the hillside forest cover class

Uruguay ndvi data at a spatial resolution of 1 km2 was obtained from the SPOT 4 & 5 vegetation (vgt-S10) product available on the vito platform42, and previously processed by Ceroni, Achkar, Gazzano & Burgueño24. The index was calculated as: ndvi = (rnir-rred) / (rnir + rred), where rnir was near infrared (nir) reflectance and rred was red reflectance. During the period of April 1998 -March 2012, VGT-S10 products provided 10-day composite ndvi images, applying the maximum-value-composite (mvc) technique to select the «best» ground reflectance values from daily ndvi. Composite data are then post-processed, which includes the incorporation of flags for bad data, clouds, and a land mask to reduce the influence of atmospheric conditions; flag data was obtained from the Status Map layer in the S10 product catalogue42. Further details regarding the processing of ndvi data at the national scale are described in Ceroni, Achkar, Gazzano & Burgueño24.

To select the ndvi data corresponding to hrf, an intersection of the national-scale ndvi pixel layer from spot imagery was made with the hrf-2015 layer. We obtained a total of 439 hrf pixels; those pixels with >90 % of hrf cover were defined as «pure pixels»28. This monthly ndvi time series reflects the behavior of ndvi in areas that are dominated by native forest and have minimal influence from other land cover classes such as tree plantations, bare soil or grasslands. Tree plantations are excluded from this analysis, whereas in previous studies ndvi of native forest and tree plantations are considered as a single class43. Based on this time series data of 168 months, we calculated the annual mean, minimum and maximum ndvi for all hillside pixels.

Finally, to evaluate the relationship between tree canopy cover and ndvi, each of the 439 pixels were assigned to one of the following canopy cover categories: 0-24 %, 25-49 %, 50-74 % and 75-100 % (Figure 1)28. Tree cover class was estimated visually for each pixel, using the tool «create network» to create a grid with cells of 250 x 250 m, which was superimposed with the hillside forest pixels over high resolution esri built-in Basemap Layer «World Imagery» in ArcMap 10.01 at a scale of 1:20,000. As this Basemap Layer includes high-resolution imagery GeoEye ikonos (1m resolution) for most of Uruguay, individual trees were easily visible to visually classify tree cover28. To compare the differences in ndvi between the four tree coverage categories, generalized linear models and the Tukey hsd post-test were used. All analyses were conducted with the software R 3.0.344.

Figure 1: Examples of 1 km2 pixels from SPOT satellite images of hillside and ravine forest (HRF) categorized by tree canopy cover into four categories (%): (a) 1-25, (b) 26-50, (c) 51-75, (d) 76-100. 

Results

The hrf were highly patchy in the Uruguayan landscape, consisting of 4226 polygons with an average area of 45 ha and a median of 24 ha in 2015. We created an updated layer for hrf in Uruguay in 2015, covering a total area of 334,487 ha (1.9 % of the terrestrial surface area of Uruguay; Figure 2). This area was less than the 384,240 ha (2.1 % of Uruguayan territory) as registered in the national land cover class product «lcc-2008»). Using the lcc-2008 layer as a baseline for forest cover, we excluded the area originally categorized as hillside forest, that were subsequently reclassified in 2015 as grasslands and afforestation. The 13 % decline (49,754 ha) in the estimated hillside forest area from 2008 to 2015 was largely due to the elimination of hillside grasslands classified as native forest in the lcc-2008 layer. These potential errors in forest classification based on Landsat 5 images in 2007-2008 increased at regional spatial scales, for example, by eco-region or by department (e.g.: we found that hillside forest cover was overestimated by 27 % in the ecoregion Sierras del Este, while hillside forest cover was underestimated by 34 % in the basaltic soil ecoregion of northern Uruguay; Table 1).

Figure 2: Map of the updated total surface area of hillside and ravine forests in Uruguay (hsr-2015) showing terrestrial ecoregions of Uruguay45

Table 1: Total surface area of hillside and ravine forest in Uruguay by Ecoregion (Fig. 1a) according to the official forest cover classification, lcc-2008, the updated hrf-2015 cover class generated here, and the net difference in area between these two cover classes. Declines in total hrf area are shaded in grey. Ecoregions proposed by Brazeiro and others45

According to the eco-regions of Uruguay45, 75 % of the surface area of hrf were distributed in the eco-region «Sierras del Este», followed by 15 % in the northern eco-region Cuesta basáltica (Figure 2). The departments with the largest forest area were: Maldonado, Lavalleja and Rocha (Table 2).

Table 2: National extent of hillside and ravine forest by administrative divisions (Department, Fig. 1b), according to the official forest cover class lcc-2008, the updated hrf-2015 cover class generated here, and the differences between them. Declines in total hrf area are shaded in grey. 

The discrepancy between the spatial extent of the hillside and ravine forest layers of 2008 and 2015 was also reflected in the difference between tree canopy cover for these two years (see Figure 2 as an example). In 2008, hillside and ravine forest 1 km2 pixels had an average canopy cover of 55 ± 25 %, with only 11 % of the pixels in the highest canopy cover class of 76-100 %. In the updated forest layer presented here, the average canopy coverage among pure forest pixels was 69 ± 23 %, with 17 % of pixels displaying 75-100 % tree cover. In this sense, the 2015 layer is more representative of dense forest than the 2008 classification of native hillside forests, due to the inclusion of previously omitted forested areas with high tree cover, as well as the exclusion of areas with scarce tree cover (0-25 %) that corresponded to other land cover classes, mainly hillside grasslands. Nonetheless, hrf remained patchy, with 44% of forest pixels in the 2015 cover class displaying a tree canopy cover less than 50 %.

During the study period, hrf pixels with more than 50% tree cover (N = 194 1 km2 spot pixels) had an overall average ndvi of 0.70 ± 0.03, varying seasonally between an average of 0.74 ± 0.03 in autumn (maximum greenness) and 0.67 ± 0.04 in winter (minimum greenness). For this same subgroup of pixels, mean ndvi in summer and spring were 0.70 ± 0.04 and 0.71 ± 0.03, respectively. The average maximum ndvi -averaged over the 14-year period for each individual pixel- occurred in May for 44% of pixels and in April for 29 % of the pixels. The average minimum ndvi -averaged over the 14-year period for each individual pixel- was largely expressed in August among 52 % of pixels and in September for 23 % of pixels.

Heterogeneity in forest density, measured as canopy cover, was reflected in differences in forest ndvi at the 1 km2 scale. There was an increase in ndvi with an increase in tree canopy cover (Figure 3). Pixels with less than 25 % tree cover had an average ndvi of 0.69 (± 0.03), while those with 100% cover had an average ndvi of 0.72 (± 0.03). The analyses showed significant differences in the average ndvi between the maximum canopy cover category (>75 %) and the two lowest cover categories (<25 % and <50 %; Tukey hsd, p <0.001). Significant differences in minimum ndvi (average of the annual minimum) were also observed among categories, except between the intermediate coverage categories (25-49 %and 50-75 %; Tukey hsd, p <0.01). Alternatively, the maximum values of ndvi were very similar regardless the canopy cover (Tukey hsd, p> 0.05, with the exception of the significant difference in ndvi between 25-49 % and 76-100 %, p <0.05).

Figure 3: Differences in the mean (bars), minimum (triangles) and maximum (dots) ndvi among the four categories of tree canopy cover of 1 km2 spot pixels of hillside and ravine forests (hrf-2015). 

Discussion

National estimates of total forest area are essential for participation in global initiatives for biodiversity conservation and carbon mitigation46. In response to the need to evaluate different methods for forest monitoring via remote sensing, this study highlights that the effective area of native forest depends upon the definition of «forest», and the spatial resolution of the images used for forest cover estimation. Specific concern arose regarding the classification of terrestrial land with less than 25 % tree canopy cover as native forest. Moreover, an adequate method for monitoring forest expansion or deforestation in areas with low tree density is clearly needed. In this regard, this study contributes to a broader discussion about the most efficient and accurate way to estimate the forest surface area in patchy landscapes or savanna biomes, highlighting the difficulties of exclusively using Landsat images (resolution 30x30 m) to quantify national forest cover. We demonstrate that the complementary use of satellite (Landsat) and high-resolution images (Google Earth) with field sampling data (nfi) is a useful strategy to achieve a more accurate delimitation of native forests in Uruguay, particularly in the case of «patchy ecosystems» such as hrf.

Hillside and ravine forest: delimitation, canopy cover and ndvi

A key aspect in determining the spatial distribution of forests in Uruguay is establishing an operational definition of forests. In the delimitation of forests with satellite images, it is useful to establish a minimum canopy cover -at a defined spatial scale-, both to determine forest area as well as to monitor areas of expansion or degradation47)(48. According to the Native Forest Registry of the General Forest Directorate in Uruguay: «Any forested area with a density greater than 200 trees/ha with a canopy cover of 50% is considered as forest»49. Nonetheless, we identified the inclusion of grasslands in the hillside forest cover class in the national land cover classification product lcc-2008 as the main source of error in the estimation of surface area. Through the classification of forest pixels by canopy cover here, we detected that 12 % of pixels had a tree cover of less than 25 % forest in 2015 (e.g.: see Figure 1). These pixels have ndvi values similar to those of pasture and livestock production systems28)(43.

The hrf of Uruguay are naturally patchy within the Uruguayan landscape, and they can be further fragmented by both historical and more recent anthropogenic activities40)(50. These forest patches are located within a transition zone between forest and grassland51, where it can be said that grasslands and forests are alternative states, where the prevalence of one environment or the other depends on complex interactions between climate, fire and land-use activities, especially livestock grazing52)(53)(54.

Our study reaffirms the fragmented and heterogeneous nature of hrf, where about half of the forest pixels analyzed had a tree canopy cover of less than 50 %. This implies that an important area of this ecosystem designated as «forest» strictly corresponds to other land cover classes, such as grasslands and rocks. In fact, a recent study reveals the intimate association between these environments, and provides evidence that the expansion of hillside forests into the surrounding grassland matrix occurs through the nucleation processes around rocks and woody nurseries54. Heterogeneity and the land cover mosaic could influence the function of these forest ecosystems, reflected by the lower ndvi values in areas with lower tree canopy cover. The complexity to quantify wooded areas in heterogeneous landscapes was also observed in African savannas19)(55. Among other tools, the use of high-resolution images and vegetation indices such as ndvi to differentiate woodland types is highlighted19.

Error in the quantification of forest extent

One of the central aspects of forest mapping is obtaining a cover class with a high degree of precision, which allows for the temporal monitoring of deforestation and carbon stocks in native forests18)(56)(57. This is one of the great challenges Uruguay faces to sustain a monitoring program of native forests over time.

Regarding the potential sources of error in estimating cover of native forest in the country, we observed that the 2015 forest cover generated here estimated the hrf area as 19 % and 27 % smaller, relative to this same category in the official layers of Land Cover for 2008 (lcc-2008) and 2011 (lcc-2011) respectively31. These differences are likely due to the images used for analysis (Landsat 5 tm vs. Landsat 8), the methodology used for the classification of forest cover (segmentation vs. visual interpretation and digitization31), and the change in scale of analysis from 1:100,000 to 1:30,000, which allowed for a more precise delimitation of the forests. It should be mentioned that, beyond the apparent decrease in the total area of mapped forests at a national scale, northern and western Uruguay showed a regional increase in hillside forest area associated with a greater inclusion of forest. These adjustments are of vital importance for the monitoring and detection of threats to this type of forest, especially under future scenarios of accelerated changes in land use58.

The overestimation of the hillside forest area by the inclusion of surrounding grasslands and the underestimation of the surface of ravine forests highlights the question of how much error is permissible in the quantification of forest area in Uruguay. If we consider a maximum error of 5 to 10 %59, other indicators of degradation and deforestation such as ndvi60 would have to be considered in order to evaluate the status of native forests.

Given the application of ndvi, among other indices, to monitor vegetation in Uruguay61)(62 it is important to recognize the factors that influence spatial variability within the «native forest» cover class. In addition to the differences identified in the ndvi among native forest types in Uruguay28, we show here that an increase in tree density within the hrf cover class was associated with higher ndvi values. As ndvi correlates with Leaf Area Index, the increase in ndvi may be the product of an increase in leaf area (or photosynthetic tissue). Similarly, in savanna landscapes of Central Africa where tree density varies from 3% to >18%, ndvi has also been identified as a good indicator of tree cover19. This work highlights the heterogeneity that a single forest cover class can have not only in structure, but in function, such as in carbon sequestration and preventing soil erosion.

Recommendations for future research

We propose the following recommendations for the delimitation of native forests in general:

    In this way, efficient national scale monitoring uses methods such as Landsat image analysis to facilitate comparison with global figures of forest extent. A representative fine-scale sub-sampling could be useful to corroborate trends in deforestation and forest expansion, and maintain a system for checking potential errors in national figures.

    Conclusions

    Currently, there is a nationwide discussion regarding the extent and trends of native forests in Uruguay, but with limited efforts to assess the appropriate methods and scales for long-term monitoring of this critical ecosystem. The patched nature of hrf, and their close interaction with grassland and scrubland environments, as well as tree plantations, make the precise delimitation of this land cover class problematic. This difficulty in establishing borders for ecosystems in transition with other land cover classes leads to errors in estimating the total cover of national native forest. Although moderate-scale national land cover classification products are valuable and widely used for both research and decision-making, in the case of spatially restricted ecosystems, a finer scale approach could detect spatiotemporal changes in their cover and function within acceptable ranges of certainty. Such is the case for native forests63, as well as highly threatened palm groves64)(65, biodiverse wetlands and other critical ecosystems that have relatively small areas within Uruguay. In order to capture forest expansion and/or degradation of hrf, we suggest that broad-scale satellite monitoring be complemented by high-resolution images, tree canopy cover indicators and fine-scale monitoring of specific patches representative of each ecoregion. These practices will ultimately contribute towards the monitoring, conservation and management of this ecosystem.

    Acknowledgements:

    We thank Ismael Díaz, Faculty of Sciences, who guided the spatial analysis, and Alejandro Brazeiro, Marcel Achkar and Santiago Baeza for their advice and support; the Biodiversity and Ecology of Conservation Group (bec) of the Faculty of Sciences, the Interdisciplinary Space of the University of the Republic, and cenur Litoral Norte, for logistical support. The authors received financial support from the National Research and Innovation Agency, anii (Post-doc Scholarship 2013-2015 «Prof. Dr. Roberto Caldeyro Barcia Fund.», No. PD_NAC_2012_1_7627), and the Sectoral Commission for Scientific Research (CSIC).

    Bibliography:

    1. World Bank. Forests Sourcebook: Practical Guidance for Sustaining Forests in Development Cooperation. Washington (DC): World Bank; 2008. 369 p. [ Links ]

    2. FAO. National Forest Monitoring Systems: Monitoring and Measurement, Reporting and Verification (M & MRV) in the context of REDD+ Activities (Internet). Roma: FAO; 2013 (cited 2018 Aug 10). 24 p. Available from: Available from: http://www.fao.org/ 3/a-bc395e.pdf . [ Links ]

    3. Malhi Y, Meir P, Brown S. Forests, carbon and global climate. Philos Trans A Math Phys Eng Sci. 2002;360:1567-91. [ Links ]

    4. Joppa LN, Loarie SR, Pimm SL. On the protection of «protected areas». Proc Natl Acad Sci USA. 2008;105:6673-8. [ Links ]

    5. Millennium Ecosystem Assessment. Ecosystem and human well-being: Synthesis. Washington (DC): Millennium Ecosystem Assessment, Island Press; 2005. 137 p. [ Links ]

    6. Mori SA, Lertzman KP, Gustafsson L. Biodiversity and ecosystem services in forest ecosystems: A research agenda for applied forest ecology. J Appl Ecol. 2016;54:12-27. [ Links ]

    7. Lefsky MA, Cohen WB, Spies TA. An evaluation of alternate remote sensing products for forest inventory, monitoring and mapping of Douglas-fir forest in western Oregon. Can J For Res. 2001;31:78-87. [ Links ]

    8. Giri C, Long J. Land-cover characterization and mapping of South America for the year 2010 using Landsat 30 m satellite data. Remote Sens. 2014;6:9494-510. [ Links ]

    9. Achard F, Ebeuchle R, Mayaux P, Stibig H-J, Bodart C, Brink A, Carboni S, Desclée B, Donnay F, Eva HD, Lupi A, Raši R, Seliger R, Simonetti D. Determination of tropical forest deforestation rates and related carbon losses from 1990-2010. Global Change Biol. 2014;20:2540-54. [ Links ]

    10. Hansen MC, Potapov PV, Moore R, Hancher M, Turubanova SA, Tyukavina A, Thau D, Stehman SV, Goetz SJ, Loveland TR, Kommaredy A, Egorov A, Chini L, Justice CO, Townshend JRG. High-Resolution Global Maps of 21st-Century Forest Cover Change. Science. 2013;342:850-3. [ Links ]

    11. Korhonen L, Korhonen KT, Rautainen M, Stenberg P. Estimation of forest canopy cover: a comparison of field measurement techniques. Silva Fenn. 2006;40(4):577-88. [ Links ]

    12. McDowell NG, Coops NC, Beck PSA, Chambers JQ, Gangodagamage C, Hicke JA, Huang C-Y, Kennedy R, Krofcheck DJ, Litvak M, Meddens AJH, Muss J, Negrón-Juarez R, Peng C, Schwantes AM, Swenson JJ, Vernon LJ, Park Williams A, Xu C, Zhao M, Running SW, Allen CD. Global satellite monitoring of climate-induced vegetation disturbances. Trends Plant Sci. 2015;20:114-23. [ Links ]

    13. Jennings SB, Brown ND, Sheil D. Assessing forest canopies and understorey illumination: Canopy closure, canopy cover and other measures. Forestry. 1999;72(1):59-74. [ Links ]

    14. Ozanne CMP, Anhuf D, Boulter SL, Keller M, Kitching RL, Korner C, Meinzer FC, Mitchell AW, Nakashizuka T, Silva Dias PL, Stork NE, Wright SJ, Yoshimura M. Biodiversity meets the atmosphere: A global view of forest canopies. Science. 2003;301:183-6. [ Links ]

    15. Wedeux BMM, Coomes DA. Landscape-scale change in forest canopy structure across a partially logged tropical peat swamp. Biogeosciences. 2015;12:6707-19. [ Links ]

    16. Coetsee C, Bond WJ, February EC. Frequent fire affects soil nitrogen and carbon in an African savanna by changing woody cover. Oecologia. 2010;162:1027-34. [ Links ]

    17. Holdo RM, Mack MC, Arnold SG. Tree canopies explain fire effects on soil nitrogen phosphorus and carbon in a savanna ecosystem. J Veg Sci. 2012;23:352-60. [ Links ]

    18. Asner GP, Knapp DE, Martin RE, Tupayachi R, Anderson CB, Mascaro J, Sinca F, Chadwick KD, Higgins M, Farfan W, Llactayo W, Silman MR. Targeted carbon conservation at national scales with high-resolution monitoring. Proc Natl Acad Sci USA . 2014;10:E5016-22. [ Links ]

    19. Gaughan AE, Holdo RM, Anderson TM. Using short-term MODIS time-series to quantify tree cover in a highly heterogeneous African savanna. Int J Remote Sens . 2013;34:6865-82. [ Links ]

    20. Joshi C, De Leeuw J, Skidmore AK, van Duren IC, van Oosten H. Remotely sensed estimation of forest canopy density: A comparison of the performance of four methods. Int J Appl Earth Obs Geoinf. 2006;8:84-95. [ Links ]

    21. Henry M, Réjou-Méchain M, Cifuentes Jara M, Wayson C, Piotto D, Westfall J, Michel Fuentes, Alice Guier F, Castañeda Lombis H, Castellanos López E, Cuenca Lara R, Cueva Rojas K, Del Águila Pasquel J, Duque Montoya A, Fernández Vega J, Jiménez Galo A, López OR, Gunnar Marklund L, Milla F, Návar Cahidez JJ, Ortiz Malavassi E, Pérez J, Ramírez Zea C, Rangel García L, Rubilar Pons R, Sanquetta C, Scott C, Zapata-Cuartas M, Saint-André L. An overview of existing and promising technologies for national forest monitoring. Ann For Sci. 2015;72:779-88. [ Links ]

    22. Paruelo JM, Garbulsky MF, Guerschman JP, Jobbagy EG. Two decades of Normalized Difference Vegetation Index changes in South America: Identifying the imprint of global change. Int J Remote Sens . 2004;25:2793-806. [ Links ]

    23. Baeza S, Baldassini P, Bagnato C, Pinto P, Paruelo JM. Caracterización del uso/cobertura del suelo en Uruguay a partir de series temporales de imágenes MODIS. Agrociencia Uruguay. 2014;18(2):95-105. [ Links ]

    24. Ceroni M, Achkar M, Gazzano I, Burgueño J. Estudio del NDVI mediante análisis multiescalar y series temporales utilizando imágenes SPOT, durante el período 1998-2012 en el Uruguay. Rev teledetec. 2015;43:31-42. [ Links ]

    25. Huete AR. Vegetation indices, remote sensing and forest monitoring. Geogr Compass. 2012;6:513-32. [ Links ]

    26. Zhu L, Southworth J. Disentangling the relationships between net primary production and precipitation in southern Africa savannas using satellite observations from 1982 to 2010. Remote Sens . 2013;5:3803-25. [ Links ]

    27. Lezama F, Baeza S, Altesor A, Cesa A, Chaneton EJ, Paruelo JM. Variation of grazing-induced vegetation changes across a large-scale productivity gradient. J Veg Sci . 2014;25:8-21. [ Links ]

    28. Lucas C, Ceroni M, Baeza S, Muñoz AA, Brazeiro A. Sensitivity of subtropical forest and savanna productivity to climate variability in South America, Uruguay. J Veg Sci . 2017;28:192-205. [ Links ]

    29. Soust P, Echeverría R, San Román D, Nebel JP. La forestación en Uruguay como complemento del desarrollo. In: Anuario OPYPA 2013(Internet). Montevideo: MGAP; 2013 (cited 2018 Aug 5). p. 623-30 Available from: Available from: http://www.mgap.gub.uy/sites/default/files/diea_anuario_2013.pdf . [ Links ]

    30. Brazeiro A. Los bosques de Uruguay y sus servicios ecosistémicos. In: Caballero N, editor. Memoria de los Foros Técnicos sobre servicios ecosistémicos en Uruguay. Montevideo: IICA; 2014. p. 19-23. [ Links ]

    31. Cal A, Álvarez A, Petraglia C, Dell’ Aqua M, López N, Fernandez VM. Mapa de Cobertura del Suelo de Uruguay. Montevideo: OPP; 2011. 52 p. [ Links ]

    32. Gautreau P. Rethinking the dynamics of woody vegetation in Uruguayan campos, 1800-2000. J Hist Geogr. 2010;36:194-205. [ Links ]

    33. FRA. Evaluación de los Recursos Forestales Mundiales 2015: Evaluación nacional In: Paquete de informe sobre los bosques 2015 (Internet). Montevideo: FAO; 2015 (cited 2017 Nov 26). p. 53-77. Available from: Available from: http://www.mgap.gub.uy/sites/default/files/multimedia/uruguayfra2015.pdf . [ Links ]

    34. Ministerio de Ganadería, Agricultura y Pesca, DIEA (UY). Censo general agropecuario 2000 (Internet). Montevideo: MGAP ; 2000 (cited 2018 Jul 5). Available from: Available from: http://www2.mgap.gub.uy/portal/page.aspx?2,diea,diea-censo-2000-antecedentes,O,es,0,Links ]

    35. Ministerio de Ganadería, Agricultura y Pesca, DIEA (UY). Censo general agropecuario 2011: Resultados definitivos (Internet). Montevideo: MGAP ; 2011 (cited 2018 Jul 5). Available from: Available from: http://www.mgap.gub.uy/sites/default/files/multimedia/censo2011.pdf . [ Links ]

    36. Díaz I, Achkar M. Estimación de superficie de monte nativo en el Litoral Norte de Uruguay mediante la utilización de imágenes satelitales LANDSAT 5TM para los años 2001-2009. Montevideo: Proyecto Monte Nativo; 2010. [ Links ]

    37. Tiscornia G, Achkar M, Brazeiro A. Efectos de la intensificación agrícola sobre la estructura y diversidad del paisaje en la región sojera de Uruguay. Ecol Austral. 2014;24:212-9. [ Links ]

    38. Google LLC. Google Earth (Internet) Version 7.1.2.2041. Mountain View (CA): USA; 2017. (cited 2018 Aug 14). Available from: Available from: https://www.google.com/intl/es-419/earth . [ Links ]

    39. United States Geological Survey. Earth Explorer (Internet) Reston (VA): U.S. Department of the Interior; 2015. (cited 2018 Apr 21). Available from: Available from: https://earthexplorer.usgs.gov . [ Links ]

    40. Brussa C, Grela I. Flora arbórea del Uruguay, con énfasis en especies de Rivera y Tacuarembó. Rivera: COFUSA; 2007. 544 p. [ Links ]

    41. Alonso PE, Bassagoda MA. Aspectos fitogeográficos y diversidad biológica de las formaciones boscosas del Uruguay. Ciência & Ambiente. 2002;24:35-50. [ Links ]

    42. VITO Remote Sens ing. VEGETATION (Internet). VGT-S10. Mol (BE): VITO; 2012 (cited 2017 Nov 3). Available from: Available from: www.spot-vegetation.com . [ Links ]

    43. Texeira M, Oyarzabal M, Pineiro G, Baeza S, Paruelo JM. Land cover and precipitation controls over long term trends in carbon gains in the grassland biome of South America. Ecosph. 2015;6(10)1-21. [ Links ]

    44. R Core Team. R: A Language and Environment for Statistical Computing (Internet). Vienna: R foundation for Statistical Computing; 2012 (cited 2017 Dec 2). Available from: Available from: http:/ /www.R-project.org/ . [ Links ]

    45. Brazeiro A, Panario D, Soutullo A, Gutierrez O, Segura A, Mai P. Identificación y delimitación de eco-regiones de Uruguay. In: Brazeiro A, editor. Eco-Regiones de Uruguay: Biodiversidad, presiones y conservación: Aportes a la Estrategia Nacional de Biodiversidad. Montevideo: Universidad de la República, Facultad de Ciencias; 2015. p. 46-59. [ Links ]

    46. FAO. Map Accuracy Assessment and Area Estimation: A Practical Guide (Internet). Roma: FAO ; 2016 (cited 2018 Mar 10). 60 p. Available from: Available from: http://www.fao.org/3/a-i5601e.pdf . [ Links ]

    47. Sasaki N, Putz FE. Critical need for new definitions of «forest» and «forest degradation» in global climate change agreements. Conserv Lett. 2009;2(5):226-32. [ Links ]

    48. FAO. FRA 2015: Terms and Definitions (Internet). Roma: FAO ; 2015 (cited 2018 Apr 15). 31 p. (Forest resources Assessment Working Paper; 180). Available from: Available from: http://www.fao.org/docrep/017/ap862e/ap862e00.pdf . [ Links ]

    49. Ministerio de Ganadería, Agricultura y Pesca, DGF (UY). Formulario e instructivo para la presentación de la solicitud de registro de bosque nativo (Internet). Montevideo: MGAP ; 2018 (cited 2018 Jun 15). Available from: Available from: http://www.mgap.gub.uy/sites/default/files/multimedia/instructivoregistronativo.pdf . [ Links ]

    50. Gautreau P, Lezama F. Clasificación florística de los bosques y arbustales de las sierras del Uruguay. Ecol Austral . 2009;19:81-92. [ Links ]

    51. Staver AC, Archibald S, Levin S. Tree cover in sub-Saharan Africa: Rainfall and fire constrain forest and savanna as alternative stable states. Ecology. 2011;92:1063-72. [ Links ]

    52. Etchebarne V, Brazeiro A. Effects of livestock exclusion in forests of Uruguay: Soil condition and tree regeneration. Forest Ecol Manag. 2016;362:120-9. [ Links ]

    53. Brussa P. Ecotono bosque-pastizal serrano: Efectos del ganado en la expansión del bosque (grade’s thesis). (Montevideo): Universidad de la República, Facultad de Agronomía; 2018. 45 p. [ Links ]

    54. Brazeiro A, Brussa P, Toranza C. Efectos del ganado en la dinámica del ecotono bosque-pastizal en paisajes serranos de Uruguay. Ecosistemas. 2018;27(3):14-23. [ Links ]

    55. Campo-Bescós MA, Muñoz-Carpena R, Southworth J, Zhu L, Waylen PR, Bunting E. Combined spatial and temporal effects of environmental controls on long-term monthly NDVI in the southern Africa Savanna. Remote Sens . 2013;5(12):6513-38. [ Links ]

    56. Dutrieux LP, Verbesselt J, Kooistra L, Herold M. Monitoring forest cover loss using multiple data streams, a case study of a tropical dry forest in Bolivia. ISPRS J Photogramm Remote Sens . 2015;107:112-25. [ Links ]

    57. Griffiths P, Jakimow B, Hostert P. Reconstructing long term annual deforestation dynamics in Pará and Mato Grosso using the Landsat archive. Remote Sens Environ. 2018;216:497-513. [ Links ]

    58. De Sy V, Herold M, Achard F, Asner GP, Held A, Kellndorfer J, Verbesselt J. Synergies of multiple remote sensing data sources for REDD+ monitoring. Curr Opin Environ Sustain. 2012;4:696-706. [ Links ]

    59. FAO. Voluntary guidelines on National Forest Monitoring (Internet). Roma: FAO ; 2017 (cited 2018 Apr 1). 61 p. Available from: Available from: http://www.fao.org/3/a-i6767e.pdf . [ Links ]

    60. Higginbottom TP, Symeonakis E. Assessing land degradation and desertification using vegetation index data: Current frameworks and future directions. Remote Sens . 2014;6(10):9552-75. [ Links ]

    61. Guido A, Diaz-Varela R, Baldassini P, Paruelo J. Spatial and Temporal Variability in Aboveground Net Primary Production of Uruguayan Grasslands. Rangeland Ecol.& Manag. 2014;67:30-8. [ Links ]

    62. Baeza S, Lezama F, Piñeiro G, Altesor A, Paruelo JM. Spatial variability of above-ground net primary production in Uruguayan grasslands: a remote sensing approach. Appl Veg Sci. 2010;13:72-85. [ Links ]

    63. Bentacourt A, Brazeiro A. Clasificación, mapeo y caracterización general de los bosques de Uruguay. In: Brazeiro A, editor. Seminario Recientes avances en investigación para la gestión y conservación del bosque nativo de Uruguay. Montevideo: Universidad de la República, Facultad de Ciencias ; 2017. p. 55-8. [ Links ]

    64. Rivas M, Barbieri RL. Buenas prácticas para el manejo sostenible del palmar de Butiá. Brasilia (DF): EMBRAPA; 2015. 65 p. [ Links ]

    65. Bortolini SV. Distribución, abundancia y estado de conservación de la palmera Butia yatay en Uruguay: Efectos de las actividades agroforestales y del cambio en el uso del suelo (master’s thesis). (Montevideo): Universidad de la República, Facultad de Agronomía ; 2018. 71 p. [ Links ]

    Author contributions: CT: methodological design, analysis of Landsat 8 images, mapping in geographic information systems (GIS), analysis of results, interpretation of results and writing. CL: analysis of forest cover, statistical processing, participation in drafting the manuscript. MC: analysis and processing of SPOT images, obtaining and processing the NDVI time series.

    Received: April 10, 2019; Accepted: May 15, 2019

    Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License