Acessibilidade / Reportar erro

Mapping vegetation on ferruginous substrates using ASTER and gamma-spectrometry images in the Iron Quadrangle, Minas Gerais

Mapeamento da vegetação sobre substratos ferruginosos utilizando imagens Aster e Gamaespectrometria no Quadrilátero Ferrífero, Minas Gerais

ABSTRACT

The Iron Quadrangle (IQ) region in Minas Gerais is remarkably geobiodiverse, despite a long history of anthropogenic pressures such as mining and urbanization, but still lacks detailed studies on the distribution of its remaining native vegetation in different substrates. In this study, we utilized Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) images, besides Gamma-spectrometry (Gamma) survey data associated with existing geological mapping (GM) and extensive fieldwork, to discriminate and quantify remnants of vegetation on ferruginous substrates in the IQ. The Maximum Likelihood (ML) algorithm was used to classify the vegetation types, thus named: open Rupestrian Field, shrubby Rupestrian Field, Capão Forest, Cerrado stricto sensu, Cerrado Field, Seasonal Forests, Pastures and Reforestation (the latter three regardless of substrate type) associated with the predominant substrates (ferruginous ironstone, phyllites, and quartzite). The use of ASTER images alone did not allow a reliable separation of ferruginous and non-ferruginous substrates, but the integration of all different data (ASTER-ML + Gamma + GM) allowed the provisional mapping of the vegetation associated with ferruginous substrates, potentially ferruginous and non-ferruginous substrates. The resulting map shows that the vegetation on ferruginous and potentially ferruginous substrates cover 8.7% and 6.9% of the IQ, respectively. The detailed analysis of the distribution and fragmentation of phytophysiognomies on ferruginous substrates is of great importance for developing strategies to conserve the geobiodiversity of the IQ, and need to be further refined by checking and field mapping by novel approaches.

Keywords:
Supervised classification; Remote sensing; Aerogeophysical data

RESUMO

A região do Quadrilátero Ferrífero (QF) em Minas Gerais apresenta notável geobiodiversidade, a despeito de uma longa história de pressões antrópicas (mineração, urbanização), mas ainda carece de estudos detalhados sobre a distribuição das formações vegetais nos diferentes substratos. No presente estudo, utilizou-se imagens Aster, levantamentos gamaespectométricos (Gama), mapeamentos geológicos existentes e trabalhos de campo, para a discriminação e quantificação dos remanescentes da vegetação sobre substratos ferruginosos e outros no QF. Utilizou-se o algoritmo da Máxima Verossimilhança (MaxVer) para a classificação das tipologias de vegetação associadas aos substratos predominantes (canga-ferruginosos, filitos e quartzitos), assim denominados: Campo Rupestre Aberto, Campo Rupestre Arbustivo, Capão Florestal, Cerrado Stricto Sensu, Campo Cerrado, Florestas Estacionais, Pastagens e Reflorestamentos, sendo os três últimos não discretizados quanto ao substrato. A utilização das imagens Aster apenas não permitiu uma separação confiável de substratos ferruginosos dos não ferruginosos, sendo necessária a utilização de produtos de processamento de dados Gama e de geologia. A integração entre as bases de informação (Aster-MaxVer + Gama + geologia) permitiu o mapeamento da vegetação associada a três tipos de substratos: ferruginoso, potencialmente ferruginoso e não-ferruginoso. O mapeamento revelou que a vegetação sobre substrato ferruginoso e potencialmente ferruginoso recobrem, respectivamente, 8,7% e 6,9% no QF. Estes dados preliminares inéditos permitem análises mais detalhadas do padrão de distribuição e fragmentação das fitofisionomias sobre substratos ferruginosos, com grande importância para as estratégias de conservação da geobiodiversidade do QF, e necessitam de checagens de campo e reconhecimento da vegetação, por meio de técnicas inovadoras.

Palavras-Chave:
Classificação supervisionada; Sensoriamento remoto; Dados aerogeofísicos

1.INTRODUCTION

The Iron Quadrangle (IQ) region has a remarkable geobiodiversity (Fernandes, 2016Fernandes GW. Ecology and Conservation of Mountaintop Grasslands in Brazil. Switzerland: Springer; 2016.) but lacks detailed mapping of its vegetation formations. Such mapping could provide technical support for decision-making regarding the conservation of remaining forest and Ferruginous Rupestrian Grassland formations. The region suffers from strong anthropogenic pressures on native ecosystems, mainly as a result of mining and urbanization (Jacobi et al., 2007Jacobi CM, Carmo FF, Vincent RC, Stehmann JR. Plant communities on ironstone outcrops: a diverse and endangered Brazilian ecosystem. Biodivers Conservation. 2007;16:2185-2200.).

In the mountainous regions, Lateritic Rupestrian Fields, also termed Ferruginous Rupestrian Fields (or Grasslands), or “canga vegetation,” develop on ferruginous substrates (Vincent, 2004Vincent RC. Florística, fitossociologia e relações entre a vegetação e o solo em áreas de campos ferruginosos no Quadrilátero Ferrífero, Minas Gerais. [Tese de Doutorado]. São Paulo: Instituto de Biociências da Universidade de São Paulo; 2004. ; Viana and Lombardi, 2007Viana PL, Lombardi JA. Florística e caracterização dos Campos Rupestres sobre canga na Serra da Calçada, Minas Gerais, Brasil. Rodriguésia. 2007;58(1):159-177.). These unique ecosystems, colonized by specialist plants adapted to oligotrophic environments, are capable of tolerating a number of severe environmental filters, such as shallow soils, severe water deficits, low fertility, high oxidized iron concentrations, and low water retention, as well as large daily thermal amplitudes, frequent fires, high sun exposure, and constant winds (Vincent, 2004Vincent RC. Florística, fitossociologia e relações entre a vegetação e o solo em áreas de campos ferruginosos no Quadrilátero Ferrífero, Minas Gerais. [Tese de Doutorado]. São Paulo: Instituto de Biociências da Universidade de São Paulo; 2004. ; Jacobi et al., 2007Jacobi CM, Carmo FF, Vincent RC, Stehmann JR. Plant communities on ironstone outcrops: a diverse and endangered Brazilian ecosystem. Biodivers Conservation. 2007;16:2185-2200.; Schaefer et al., 2016Schaefer CE, Candido HG, Corrêa GR, Nunes JA, Arruda DM. Soils Associated with Rupestrian Grasslands. In: Fernandes G(ed.). Ecology and Conservation of Mountaintop grasslands in Brazil. 1ed.: Springer Intern. Publishing; 2016. p.55-69.). The vegetation shows several anatomical, morphological, physiological, and reproductive adaptations that allow it to survive in these environments (Alves and Kolbek, 1994Alves RJV, Kolbek J. Plant species endemism in savanna vegetation on table mountains (campo rupestre) in Brazil. Vegetatio. 1994;113(2):125-139).

The plant communities in the IQ vary as a function of the substrate type, so that the spectral attributes of soils and surface rocks, in some cases, may assist in mapping the vegetation. Remote-sensing techniques, by identifying subtle changes in vegetative cover, can enable the identification of changes in substrate conditions, establishing the rock-soil-vegetation association (Almeida Filho, 1984Almeida Filho R. Multiseasonal and geobotanical approach in remote detection of albitized-greisenized areas in the Serra da Pedra Branca granitic Massif, Goiás State, Brazil. Economic Geology, 1984;79(8):1914-1920). The advent of Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) sensor images has enabled discrimination of geological as well as vegetation targets (Lima et al., 2005Lima BEM, Almeida Filho R, Galvão LS, Meneses PR. Avaliação das imagens do sensor ASTER para discriminação espectral de variações faciológicas no granito Serra Branca, Goiás. In: Anais XII Simpósio Brasileiro de Sensoriamento Remoto; 16-21 de abril 2005;Goiânia, Brasil: INPE; 2005. ; Gil et al., 2014Gil A, Yu Q, Abadi M, Calado H. Using Aster multispectral imagery for mapping woody invasive species in Pico da Vara Natural Reserve (Azores Islands, Portugal). Revista Árvore. 2014;38(3):391-401.). Some studies (Rouskov et al., 2005Rouskov K, Popov K, StoyKov S, Yamaguchi Y. Some applications of the remote sensing in geology by using of ASTER images. Scientific Conference ‘‘Space, Ecology, Safety’’. 2005;167-173.; Rajendran et al., 2011Rajendran S, Thirunavukkarasu A, Balamurugan G, Shankar K. Discrimination of iron ore deposits of granulite terrain of Southern Peninsular India using ASTER data. Journal of Asian Earth Sciences. 2011;41(1):99-106.) used satellite sensors (Landsat Thematic Mapper and ASTER-Terra) for the identification and discrimination of iron-rich deposits through the composition of multispectral indexes. However, this technique becomes limited when such regions are vegetated, which minimizes the effects of the energy reflected by the ferruginous substrates.

Additionally, the gamma-spectrometric (Gamma) data obtained by aerogeophysical surveys enables the elimination of the effects of vegetation cover and the direct discrimination of ferruginous substrates by inferences regarding the geochemical characteristics of the rocks. Also, Gamma spectrometry responds to concentrations of potassium (K) radioisotopes, uranium (U) and thorium (Th) series radioisotopes. In rocks and soils; these concentrations are directly proportional to the intensity of the gamma radiation emitted by their radioactive decay (Wilford et al., 1997Wilford JR, Bierwirth PN, Craig MA. Application of gamma-ray spectrometry in soil/regolith mapping and geomorphology. AGSO Journal of Australian Geology & Geophysics. 1997;17(2):201-216.; Santos et al., 2008Santos ACL, Menezes PTL, Nascimento CTC. Gamaespectrometria aplicada em solo agrícola no noroeste do Estado do Rio de Janeiro. Revista Brasileira de Geofísica. 2008;26(2):181-194 ), and are frequently associated with the geochemical signature of substrates. This allows, for example, the separation of ferruginous and non-ferruginous substrates, with different geochemical compositions. However, the use of Gamma data has some limitations in terms of distinguishing certain substrates based on their similar responses, unresolved radioactive barriers, or differences in soil moisture; therefore, these data should be used with caution, and preferably combined with all available information (Wilford et al., 1997Wilford JR, Bierwirth PN, Craig MA. Application of gamma-ray spectrometry in soil/regolith mapping and geomorphology. AGSO Journal of Australian Geology & Geophysics. 1997;17(2):201-216.) in regolith studies.

In this sense, this work aimed to apply remote-sensing techniques associated with geological data (Lobato et al., 2005Lobato LM, Baltazar OF, Reis LB, Achtschin AB, Baars FJ, Timbó MA, et al Projeto Geologia do Quadrilátero Ferrífero - Integração e Correção Cartográfica em SIG com Nota Explicativa. Belo Horizonte: CODEMIG; 2005. ) with ASTER images and Gamma data to discriminate vegetation remnants on ferruginous substrates and other substrates, in the IQ region.

2.MATERIALS AND METHODS

2.1.STUDY AREA

The IQ is located in central Minas Gerais, with 7,800 km2 in area. All procedures were conducted considering a 5-km buffer from the IQ limits (Figure 1). According to a geological cartography by Lobato et al. (2005)Lobato LM, Baltazar OF, Reis LB, Achtschin AB, Baars FJ, Timbó MA, et al Projeto Geologia do Quadrilátero Ferrífero - Integração e Correção Cartográfica em SIG com Nota Explicativa. Belo Horizonte: CODEMIG; 2005. , approximately 6.4% of the IQ consists of iron-rich formations (ferruginous substrates) that are associated with generally shallow soils, where a predominant rupestrian vegetation cover ranges from field to cerrado to upper montane forest. A range of characteristic soils can be found in these areas, showing great landscape and geoenvironmental diversity with the occurrence of a ferruginous substrate (Schaefer et al., 2008Schaefer CEGR, Mendonça BAF, Corrêa GR, Ribeiro ASS. Solos desenvolvidos sobre canga ferruginosa no Quadrilátero Ferrífero, Minas Gerais. In: Simpósio Afloramentos Ferruginosos no Quadrilátero Ferrífero, Belo Horizonte: Universidade Federal de Minas Gerais. Instituto de Ciências Biológicas. 2008;1:109-121., Schaefer et al., 2016Schaefer CE, Candido HG, Corrêa GR, Nunes JA, Arruda DM. Soils Associated with Rupestrian Grasslands. In: Fernandes G(ed.). Ecology and Conservation of Mountaintop grasslands in Brazil. 1ed.: Springer Intern. Publishing; 2016. p.55-69.).

Figure 1
Limits of the Iron Quadrangle. ASTER image scenes (with the respective date and numeration: YYMMDD(XX)) and ternary image (CMYK) of the radioisotopes K (C), eTh (M) and eU (Y) of the gamma-ray image.
Figura 1
Limites do Quadrilátero Ferrífero, cenas das imagens Aster (com as respectivas datas e numeração: AAMMDD(XX)) e imagem ternária (CMYK) dos radioisótopos K(C), eTh (M) e eU(Y) da gamaespectometria.

The denomination “ferruginous substrate,” in this study, is broad and more comprehensive than the word “canga”, and comprises a range of substrates (fresh rock, altered rock, sedimentary cover, canga, and soils) associated with the following lithotypes in the geological map (Lobato et al., 2005Lobato LM, Baltazar OF, Reis LB, Achtschin AB, Baars FJ, Timbó MA, et al Projeto Geologia do Quadrilátero Ferrífero - Integração e Correção Cartográfica em SIG com Nota Explicativa. Belo Horizonte: CODEMIG; 2005. ): canga, detrital-lateritic cover, , iron oxide supergene concentrations, hematite bodies, ferruginous dolomites, iron formations, hard hematite, itabirite, dolomitic itabirite, laterite and ferruginous detritus, hard hematite lenses, magnetite, high content iron ore, hematite ore, residual lateritic soil, and associated colluvium and eluvium. However, these lithotypes may occur as secondary component in non-ferruginous geological units, so that the term “ferruginous substrate” is not limited to the delineation of the aforementioned lithotypes.

2.2.PRE-PROCESSING OF ASTER IMAGES

Ten ASTER scenes were used, covering the 5-km IQ buffer (Figure 1). ASTER bands 1 to 9 were those processed comprising the visible and near-infrared (VNIR) and shortwave infrared (SWIR) wavelengths. The SWIR bands, with a 30-m spatial resolution, were resampled using the nearest neighbor method to 15 m for compatibility with the VNIR bands (15 m). The images were georeferenced from terrestrial control points and orthorectified GeoEye images; the final data were projected in zone 23 K SIRGAS2000 UTM. ASTER images were then converted to radiance values based on the maximum and minimum radiance values of each band. These procedures were performed using ArcGIS 10 software. Then, the Fast Line-of-sight Atmospheric Analysis of Hypercubes (FLAASH) method of the ENVI 5.0 software was used for the atmospheric correction of the ASTER images.

2.3.SUPERVISED CLASSIFICATION

In addition to bands 1 to 9 of the ASTER images, multispectral indexes were included in the supervised classification, in which the following band ratios were considered to show certain spectral patterns. 1) To show areas with ferruginous substrates, the red-green-blue (RGB) band composition R = band 2/band 1, G = band 4/band 3, and B = band 4/band 5 (Rouskov et al., 2005Rouskov K, Popov K, StoyKov S, Yamaguchi Y. Some applications of the remote sensing in geology by using of ASTER images. Scientific Conference ‘‘Space, Ecology, Safety’’. 2005;167-173.), was used and 2) to aid in vegetation mapping, the Normalized Difference Vegetation Index (NDVI) was used ((band 3 - band 2) / (band 3 + band 2)) (Rouse et al., 1974Rouse JW, Haas RH, Schell JA, Deering DW. Monitoring vegetation systems in the great plains with ERTS. In: Earth Resources Technology Satellite-1 Symposium, 3, Washington: NASA.1974;1:309-317.). The Maximum Likelihood algorithm (ML) was used and the following vegetation cover classes were separated and subdivided as a function of the substrate type (ferruginous, quarzitic, and phyllitic): Open Rupestrian Field (ORF), Shrubby Rupestrian Field (SRF), Capão Forest (CF), Cerrado Field (CF), and Cerrado Stricto Sensu (CSS). Training and validation samples were collected from all phytophysiognomies based on geological data (Dorr, 1969Dorr JVN. Physiographic, stratigraphic and structural development of Quadrilátero Ferrífero, Minas Gerais, Brazil. Washington: U.S. Government Printing Office; 1969. ; Lobato et al., 2005Lobato LM, Baltazar OF, Reis LB, Achtschin AB, Baars FJ, Timbó MA, et al Projeto Geologia do Quadrilátero Ferrífero - Integração e Correção Cartográfica em SIG com Nota Explicativa. Belo Horizonte: CODEMIG; 2005. ) and associated with the following substrates: ferruginous, quarzitic, or phyllitic (including serpentinites, metabasalts, and shales).

A large part of the area, totaling 894 control points, was covered to collect training and validation samples using a Garmin 60CSx Global Positioning System receiver. In areas of difficult access, samples were collected as field truths based on the interpretation of high-resolution images (GeoEye) of the Native Vegetation Coverage Map of the State of Minas Gerais (Scolforo et al., 2006Scolforo JR, Carvalho LMT. Mapeamento e Inventário da Flora e dos Reflorestamentos de Minas Gerais. Lavras: UFLA/IEF, 2006. ) and through queries to the Google Earth image bank. Additionally, samples were collected in regions where the following vegetation classes had no association with lithology: lakes, pastures, and reforestation (Eucalyptus and/or Pinus). Training and validation samples were collected in a polygon format consisting of 1,152 and 775 polygons, respectively. Classes of urban areas and areas influenced by mining were not included in the supervised classification and were overlapped after final processing via manual mapping over GeoEye images.

To evaluate the classification accuracy, validation samples were used and the Kappa index and Global Accuracy were calculated. The Kappa index is a statistic that measures the agreement between field truths and a classified map, noting the map’s legitimacy (Congalton, 1991Congalton RG. A review of assessing the accuracy of classifications of remotely sensed data. Remote Sensing of Environment. 1991;37(1):35-46). These values were classified as proposed by Monserud and Leemans (1992)Monserud RA, Leemans R. Comparing global vegetation maps with the kappa statistic. Ecological Modelling. 1992;62(4):275-293.. From only the training samples, the Kappa index was evaluated for all band and index combinations to identify the influence of each band or index on the classifier settlement. In addition, following the classification, 40-pixel (~1-ha) groupings were eliminated and replaced by the most representative neighboring pixels. For all these operations, ArcGIS 10 software was used.

2.4. GAMMA-SPECTROMETRY SURVEY

Continuous surfaces were generated through interpolation for data enhancement using the minimal curvature algorithm of the radioisotopes K, eTh, and eU to highlight the lateritic covering and colluvium of ferrous formations (Boyle, 1982Boyle RW. Geochemical prospecting for thorium and uranium deposits. Developments in Economic Geology. 1982;16:71-78; Wildford et al., 1997). Ternary radioisotopic images were also generated using the C(K)M(eTh)Y(eU)K color scale. The ternary Gamma image was cropped to the 5-km IQ buffer (Figure 1) and an unsupervised classification was performed using the IsoCluster method resulting in 10 classes. Then, the classes were superimposed onto the available survey geological units (Lobato et al., 2005Lobato LM, Baltazar OF, Reis LB, Achtschin AB, Baars FJ, Timbó MA, et al Projeto Geologia do Quadrilátero Ferrífero - Integração e Correção Cartográfica em SIG com Nota Explicativa. Belo Horizonte: CODEMIG; 2005. ) and correspondence with the ferruginous and non-ferruginous substrates was visually assessed, allowing the image to be reclassified from 10 to only 2 classes by grouping. Thus, the classes obtained during the Gamma analysis were used to complement and validate the separation of ferruginous and non-ferruginous substrates from the phytophysiognomic classes obtained from the classification of each ASTER scene. These procedures were performed using the ArcGIS 10 software.

2.5. GEOLOGICAL DATA

From the available IQ geological map (Lobato et al., 2005Lobato LM, Baltazar OF, Reis LB, Achtschin AB, Baars FJ, Timbó MA, et al Projeto Geologia do Quadrilátero Ferrífero - Integração e Correção Cartográfica em SIG com Nota Explicativa. Belo Horizonte: CODEMIG; 2005. ), two additional characteristics were processed and used for data refinement and integration as follows: 1) Geological units with subordinate ferruginous formations: these do not correspond to the actual ferruginous substrates, but show inclusions of banded iron formations, represented by the following geological units: the Quebra Osso Group, Nova Lima Group (Catarina Mendes, Córrego da Paina, Córrego do Sítio, Mestre Caetano, Mindé, Morro Vermelho, Ouro Fino and Santa Quitéria Units), Caraça Group (Batatal Formation), Piracicaba Group (Cercadinho Formation and Fecho do Funil Formation), Sabará Group, Itacolomi Group (Santo Antônio Formation and Itacolomi Facies Units), and Cambotas Formation; detritic-lateritic cover unit (Lobato et al., 2005Lobato LM, Baltazar OF, Reis LB, Achtschin AB, Baars FJ, Timbó MA, et al Projeto Geologia do Quadrilátero Ferrífero - Integração e Correção Cartográfica em SIG com Nota Explicativa. Belo Horizonte: CODEMIG; 2005. ). 2) Excluded geological units with no Fe-rich formations: which refer to geological formations without the presence of banded iron formations, used as references to erroneously reclassify the mapped areas as ferruginous substrates. These were the Bação Complex, Belo Horizonte Complex, Bonfim Complex, Caeté Complex, Guanhães Complex, Santa Bárbara Complex, Caraça Group (Moeda Formation, only the core area of Caraça Mountain), Gnaisse Souza Noschese Unit, Granito Borrachudos Unit, Granito Peti Unit and Rocha Intrusiva Unit (Lobato et al., 2005Lobato LM, Baltazar OF, Reis LB, Achtschin AB, Baars FJ, Timbó MA, et al Projeto Geologia do Quadrilátero Ferrífero - Integração e Correção Cartográfica em SIG com Nota Explicativa. Belo Horizonte: CODEMIG; 2005. ). This filter was mainly used for areas where there was no Gamma data (Figure 1), i.e. the north-northwest and north-northeast regions of the study area.

3.RESULTS

3.1.SUPERVISED CLASSIFICATION AND VALIDATION

Scene 070826 (02), which covers most of the IQ (Figure 1) and has no cloud cover, was used as a reference for the supervised classification. It was the basis for proposing the best band combinations and multispectral indexes for the Kappa index evaluation. Classification was performed considering all combinations of individual bands and multispectral indexes comprising the best Kappa indexes. Thus, six ASTER bands comprising the best Kappa indices - namely bands 1, 2, 3, 4, 6, and 7 - were selected. ASTER sensor bands 5, 8, and 9 were not included because they did not provide combinations with high Kappa indexes. Band ratios were also selected to show ferruginous substrates according to Rouskov et al. (2005)Rouskov K, Popov K, StoyKov S, Yamaguchi Y. Some applications of the remote sensing in geology by using of ASTER images. Scientific Conference ‘‘Space, Ecology, Safety’’. 2005;167-173. and the NDVI. The band 4/band 5 ratio, suggested by Rouskov et al. (2005)Rouskov K, Popov K, StoyKov S, Yamaguchi Y. Some applications of the remote sensing in geology by using of ASTER images. Scientific Conference ‘‘Space, Ecology, Safety’’. 2005;167-173., did not change the Kappa index value, thereby also showing that it provided no contribution to the classification. Thus, this was the only ratio not adopted for the other scenes.

Thus, by evaluating the confounding matrix of the reference scene, it was observed that, in general, the Ferruginous Rupestrian Fields (Open and Shrubby) were more spectrally separated from the Rupestrian Fields compared to that of the other lithologies, that is, they were better identified. In this regard, a cluster was proposed in which the phytophysiognomies on ferruginous substrates were maintained and the phytophysiognomies on non-ferruginous substrates were grouped (quartzites and phyllites), which improved the Kappa index in all scenes. According to Monserud and Leemans (1992)Monserud RA, Leemans R. Comparing global vegetation maps with the kappa statistic. Ecological Modelling. 1992;62(4):275-293., Kappa values obtained for the final classification of this grouping were considered reasonable, good, or very good except for those scenes with lower spatial expression (Table 1 and Figure 1). For the final analysis of the supervised classification, the Kappa index was calculated through the validation samples for the entire study area considering only the vegetation samples on ferruginous and non-ferruginous substrates. The final Kappa index obtained was 0.61 and the overall accuracy was 0.64, showing good agreement between the field truths and classified map according to Monserud and Leemans (1992)Monserud RA, Leemans R. Comparing global vegetation maps with the kappa statistic. Ecological Modelling. 1992;62(4):275-293..

Table 1
Kappa indices and overall accuracy for each scene considering all classes and groupings.
Tabela 1
Índices Kappa e exatidão global para cada cena considerando todas as classes e o agrupamento.

3.2.INTEGRATION OF DATA

The integration of ML classification data with Gamma and geological data allowed the separation of three substrate types: ferruginous, potentially ferruginous, and non-ferruginous. Figure 2 shows a synthesis of data integration used as a differentiating criterion between ferruginous and potentially ferruginous substrates. After integrating the data, areas in the geological map considered as Fe-rich lithotypes and those indicated as ferruginous substrate in the three information sources were considered areas with vegetation on a ferruginous substrate. In a complementary manner, the ferruginous substrate polygons indicated by the integration of the data extrapolating and intercepting the ferruginous lithotypes presented in the geological map (Lobato et al., 2005Lobato LM, Baltazar OF, Reis LB, Achtschin AB, Baars FJ, Timbó MA, et al Projeto Geologia do Quadrilátero Ferrífero - Integração e Correção Cartográfica em SIG com Nota Explicativa. Belo Horizonte: CODEMIG; 2005. ) had their limits considered as the vegetative cover on the ferruginous substrate.

Figure 2
Flowchart synthesis of the intersection of the information planes. The positive (+) signs indicate the vegetation mapped on ferruginous substrates.
Figura 2
Fluxograma síntese do cruzamento dos planos de informação. Os sinais positivos (+) indicam a vegetação mapeada sobre substratos ferruginosos.

Areas with a potentially ferruginous substrate were associated with a positive indication in two databases, i.e., the ASTER images and the Gamma classifications; by the ASTER image classification and the geological map as a unit with iron formations at the subordinate level; or in the classification of Gamma data and the geological map as a unit with subordinate iron formations. The areas with a non-ferruginous substrates correspond to the remaining areas indicated and also excluded areas of other geological units such as gneisses, granites, and other associations with no known occurrence of Fe-rich formations, particularly for areas without Gamma coverage (Figure 1).

3.3. QUANTIFICATION OF VEGETATION ON FERRUGINOUS SUBSTRATES

The remaining areas of ecosystems developed on a ferruginous substrate cover 676.9 km2 (8.7%) of the IQ (Table 2 and Figure 3). The most representative phytophysiognomies are the Capão Forest /Seasonal Forest and the Shrubby Rupestrian Field followed by the Open Rupestrian Field, Cerrado Stricto Sensu, and Cerrado Field, respectively (Table 2). The areas impacted by anthropogenic activity are the urban areas, areas affected by mining, pastures, and reforestation which together correspond to 24.6% of the total IQ area, proving the great pressure exerted on the natural ecosystems of this region.

Table 2
Total and relative area for each phytophysiognomy associated with ferruginous, potentially ferruginous, and non-ferruginous substrates mapped in the Iron Quadrangle, Minas Gerais.
Tabela 2
Área total e relativa para cada fitofisionomia associada aos substratos ferruginosos, potencialmente ferruginosos e não-ferruginosos mapeadas no Quadrilátero Ferrífero, Minas Gerais.

Figure 3
Final vegetation map of ferruginous and potentially ferruginous substrates of the IQ region.
Figura 3
Mapa final da vegetação sobre substrato ferruginoso e potencialmente ferruginoso da região do QF.

4.DISCUSSION

The strong anthropogenic pressure on vegetation remmants associated with the ferruginous substrates in the IQ region (Jacobi et al., 2007Jacobi CM, Carmo FF, Vincent RC, Stehmann JR. Plant communities on ironstone outcrops: a diverse and endangered Brazilian ecosystem. Biodivers Conservation. 2007;16:2185-2200.), which represent extremely important areas for biodiversity conservation in the transition between two large Brazilian hotspots, the cerrado and Atlantic forest (Myers et al., 2000Myers N, Mittermeier RA, Mittermeier CG, Fonseca GAB, Kent J. Biodiversity hotspots for conservation priorities. Nature. 2000;403:853-845.; Fernandes, 2016Fernandes GW. Ecology and Conservation of Mountaintop Grasslands in Brazil. Switzerland: Springer; 2016.), justifies detailed studies of the distribution and conservation degree of these phytophysiognomies. In this sense, the use of remote sensing techniques enables identification of these phytophysiognomies and associated substrates (pedological and geological) (Yamaguchi and Naito, 2003Yamaguchi Y, Naito C. Spectral indices for lithologic discrimination and mapping by using the ASTER SWIR bands. International Journal of Remote Sensing. 2003;24(22):4311-4323; Kalinowski and Oliver, 2004Kalinowski A, Oliver S. ASTER Mineral Index Processing. Remote Sensing Application Geoscience Australia. Australian Government Geoscience Website: http://www.ga.gov.au/image_cache/GA7833.pdf; 2004.
http://www.ga.gov.au/image_cache/GA7833....
; Rouskov et al., 2005Rouskov K, Popov K, StoyKov S, Yamaguchi Y. Some applications of the remote sensing in geology by using of ASTER images. Scientific Conference ‘‘Space, Ecology, Safety’’. 2005;167-173.; Rajendran et al., 2011Rajendran S, Thirunavukkarasu A, Balamurugan G, Shankar K. Discrimination of iron ore deposits of granulite terrain of Southern Peninsular India using ASTER data. Journal of Asian Earth Sciences. 2011;41(1):99-106.). For some remote orbital sensors obtained through reflectance signals from objects on the Earth’s surface (passive sensors), vegetation can interfere or hinder substrate identification, particularly in a closed canopy. However, other sensors, such as gamma spectroscopy, can penetrate vegetation and obtain geochemical data from the upper 30 cm of soils (Minty, 1997Minty BRS. Fundamentals of airborne gamma-ray spectrometry. AGSO Journal of Australian Geology & Geophysics. 1997;17(2):39-50.; Wilford et al., 1997Wilford JR, Bierwirth PN, Craig MA. Application of gamma-ray spectrometry in soil/regolith mapping and geomorphology. AGSO Journal of Australian Geology & Geophysics. 1997;17(2):201-216.).

Regarding the images with spectral data from the ASTER sensor, the combination of bands with the best Kappa indexes-initially calculated with the training samples only-includes six ASTER bands. ASTER sensor bands 5, 8, and 9 were not included because they were not predominant among combinations with higher Kappa indexes. Rowan and Mars (2003)Rowan LC, Mars JC. Lithologic mapping in the Mountain Pass, California area using Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data. Remote Sensing of Environment. 2003;84(3):350-366. suggest a good spectral separation of lithological categories for ASTER sensor bands 5 to 9 except for ferruginous deposits. The VNIR bands (Bands 1, 2, and 3) have information regarding metal absorption, particularly iron (Adams, 1974Adams JB. Visible and near-infrared diffuse reflectance spectra of pyroxenes as applied to remote sensing of solid objects in the solar system. Journal of Geophysical Research. 1974;79(32):4829-4836; Rowan et al., 1995Rowan LC, Bowers TL, Crowley JK, Antón-Pacheco C, Gumiel P, Kingston MJ. Analysis of Airborne Visible-Infrared Imaging Spectrometer (AVIRIS) data of the Iron Hill, Colorado, carbonatite - alkalic igneous rock complex. Economic Geology. 1995;90:1966-1982.), and chlorophyll absorption during vegetation photosynthesis (Knipling, 1970Knipling EB. Physical and physiological basis for the reflectance of visible and near-infrared radiation from vegetation. Remote Sensing of Environment. 1970;1(3):155-159.). These data corroborate the greater importance of these bands for mapping these phytophysiognomies. In this regard, the grouping of phytophysiognomies in classes on ferruginous substrates and differentiated phytophysiognomies on non-ferruginous substrates (quartzites and phyllites) is also important for the improvement in the Kappa index in all scenes (Table 1), as it also favors the aspects that show phytophysiognomies related to the ferruginous substrates.

In addition to the combination of ASTER bands, the inclusion of band ratios to highlight the ferruginous substrates (Rouskov et al., 2005Rouskov K, Popov K, StoyKov S, Yamaguchi Y. Some applications of the remote sensing in geology by using of ASTER images. Scientific Conference ‘‘Space, Ecology, Safety’’. 2005;167-173.) and NDVI were also chosen as those with the highest Kappa indexes. The 4/5 band ratio was also noted by Bierwirth (2002)Bierwirth PN. Evaluation of ASTER satellite data for geological applications. Consultancy Report to Geoscience Australia; 2002. as important in the identification of laterites although it was not influential in this study. NDVI significantly contributes to the identification of different vegetation types as it is influenced by the productivity and photosynthetic dynamics of phytophysiognomies (Rouse et al., 1974Rouse JW, Haas RH, Schell JA, Deering DW. Monitoring vegetation systems in the great plains with ERTS. In: Earth Resources Technology Satellite-1 Symposium, 3, Washington: NASA.1974;1:309-317.; Petorelli et al., 2005Petorelli N, Vik JO, Mysterud A, Gaillard JM, Tucker CJ, Stenseth NC. Using the satellite-derived NDVI to assess ecological responses to environmental change. Trends in Ecology and Evolution. 2005;20(9):503-510.).

The definition of areas with greater uncertainties, defined as potentially ferruginous substrates, reinforces the difficulty of mapping these vegetation typologies using spectral data, particularly when associated with specific substrates, as in the IQ region where there are great variations in geological units with different degrees of metamorphism (Lobato et al., 2005Lobato LM, Baltazar OF, Reis LB, Achtschin AB, Baars FJ, Timbó MA, et al Projeto Geologia do Quadrilátero Ferrífero - Integração e Correção Cartográfica em SIG com Nota Explicativa. Belo Horizonte: CODEMIG; 2005. ). Thus, Gamma data, in addition to assisting geological and pedological surveys (Vasconcellos et al., 1994Vasconcellos RM, Metelo MJ, Motta AC, Gomes RD. Geofísica em levantamentos geológicos no Brasil. Rio de Janeiro: CPRM, 1994. ; Wilford et al., 1997Wilford JR, Bierwirth PN, Craig MA. Application of gamma-ray spectrometry in soil/regolith mapping and geomorphology. AGSO Journal of Australian Geology & Geophysics. 1997;17(2):201-216.; McBratney et al., 2003Mcbratney AB, Mendonça Santos ML, Minasny B. On digital soil mapping. Geoderma. 2003;117(1-2):3-52.; Santos et al., 2008Santos ACL, Menezes PTL, Nascimento CTC. Gamaespectrometria aplicada em solo agrícola no noroeste do Estado do Rio de Janeiro. Revista Brasileira de Geofísica. 2008;26(2):181-194 ), can be extremely important for mapping phytophysiognomies associated with different substrates.

Therefore, this mapping allows for a more comprehensive analysis of the distribution and fragmentation pattern of the phytophysiognomies on ferruginous substrates which may generate information of great importance for future planning. These data can be combined with more in-depth floristic studies to support future IQ geobiodiversity conservation strategies and provide a more reliable predictability of the impacts caused by anthropogenic actions on these ecosystems. Additional work, such as reducing data using Principal Component Analysis or using more specific remote sensors such as Magnetometry, is required to improve the classification accuracy, particularly for a more precise definition of potentially ferruginous areas. These results should be taken with caution due to confusing identification of potentially ferruginous substrates, and further vegetation mapping combined with field recognition, must be done, employing novel approaches.

5.CONCLUSIONS

Selected ASTER sensor bands, band ratios for ferruginous substrates, and NDVI were combined with Gamma and geological data allowed a consistent and reliable vegetation mapping of ferruginous, non-ferruginous and potentially ferruginous substrates, the latter with lower confidence level. This mapping shows the complexity of the study region and allows detailed analysis of the distribution and vegetation fragmentation pattern of ferruginous substrates, which is of great importance for conservation strategies in the remarkably geobiodiverse IQ region in Minas Gerais, one of the largest iron mining area worldwide.

6.ACKNOWLEDGMENTS

We would like to thank Vale and Bioma Meio Ambiente companies. Financial support from CAPES, CNPq and FAPEMIG is acknowledged.

7.REFERÊNCIAS

  • Adams JB. Visible and near-infrared diffuse reflectance spectra of pyroxenes as applied to remote sensing of solid objects in the solar system. Journal of Geophysical Research. 1974;79(32):4829-4836
  • Almeida Filho R. Multiseasonal and geobotanical approach in remote detection of albitized-greisenized areas in the Serra da Pedra Branca granitic Massif, Goiás State, Brazil. Economic Geology, 1984;79(8):1914-1920
  • Alves RJV, Kolbek J. Plant species endemism in savanna vegetation on table mountains (campo rupestre) in Brazil. Vegetatio. 1994;113(2):125-139
  • Bierwirth PN. Evaluation of ASTER satellite data for geological applications. Consultancy Report to Geoscience Australia; 2002.
  • Boyle RW. Geochemical prospecting for thorium and uranium deposits. Developments in Economic Geology. 1982;16:71-78
  • Congalton RG. A review of assessing the accuracy of classifications of remotely sensed data. Remote Sensing of Environment. 1991;37(1):35-46
  • Dorr JVN. Physiographic, stratigraphic and structural development of Quadrilátero Ferrífero, Minas Gerais, Brazil. Washington: U.S. Government Printing Office; 1969.
  • Fernandes GW. Ecology and Conservation of Mountaintop Grasslands in Brazil. Switzerland: Springer; 2016.
  • Gil A, Yu Q, Abadi M, Calado H. Using Aster multispectral imagery for mapping woody invasive species in Pico da Vara Natural Reserve (Azores Islands, Portugal). Revista Árvore. 2014;38(3):391-401.
  • Jacobi CM, Carmo FF, Vincent RC, Stehmann JR. Plant communities on ironstone outcrops: a diverse and endangered Brazilian ecosystem. Biodivers Conservation. 2007;16:2185-2200.
  • Kalinowski A, Oliver S. ASTER Mineral Index Processing. Remote Sensing Application Geoscience Australia. Australian Government Geoscience Website: http://www.ga.gov.au/image_cache/GA7833.pdf; 2004.
    » http://www.ga.gov.au/image_cache/GA7833.pdf
  • Knipling EB. Physical and physiological basis for the reflectance of visible and near-infrared radiation from vegetation. Remote Sensing of Environment. 1970;1(3):155-159.
  • Lima BEM, Almeida Filho R, Galvão LS, Meneses PR. Avaliação das imagens do sensor ASTER para discriminação espectral de variações faciológicas no granito Serra Branca, Goiás. In: Anais XII Simpósio Brasileiro de Sensoriamento Remoto; 16-21 de abril 2005;Goiânia, Brasil: INPE; 2005.
  • Lobato LM, Baltazar OF, Reis LB, Achtschin AB, Baars FJ, Timbó MA, et al Projeto Geologia do Quadrilátero Ferrífero - Integração e Correção Cartográfica em SIG com Nota Explicativa. Belo Horizonte: CODEMIG; 2005.
  • Mcbratney AB, Mendonça Santos ML, Minasny B. On digital soil mapping. Geoderma. 2003;117(1-2):3-52.
  • Minty BRS. Fundamentals of airborne gamma-ray spectrometry. AGSO Journal of Australian Geology & Geophysics. 1997;17(2):39-50.
  • Monserud RA, Leemans R. Comparing global vegetation maps with the kappa statistic. Ecological Modelling. 1992;62(4):275-293.
  • Myers N, Mittermeier RA, Mittermeier CG, Fonseca GAB, Kent J. Biodiversity hotspots for conservation priorities. Nature. 2000;403:853-845.
  • Petorelli N, Vik JO, Mysterud A, Gaillard JM, Tucker CJ, Stenseth NC. Using the satellite-derived NDVI to assess ecological responses to environmental change. Trends in Ecology and Evolution. 2005;20(9):503-510.
  • Rajendran S, Thirunavukkarasu A, Balamurugan G, Shankar K. Discrimination of iron ore deposits of granulite terrain of Southern Peninsular India using ASTER data. Journal of Asian Earth Sciences. 2011;41(1):99-106.
  • Rouse JW, Haas RH, Schell JA, Deering DW. Monitoring vegetation systems in the great plains with ERTS. In: Earth Resources Technology Satellite-1 Symposium, 3, Washington: NASA.1974;1:309-317.
  • Rouskov K, Popov K, StoyKov S, Yamaguchi Y. Some applications of the remote sensing in geology by using of ASTER images. Scientific Conference ‘‘Space, Ecology, Safety’’. 2005;167-173.
  • Rowan LC, Mars JC. Lithologic mapping in the Mountain Pass, California area using Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data. Remote Sensing of Environment. 2003;84(3):350-366.
  • Rowan LC, Bowers TL, Crowley JK, Antón-Pacheco C, Gumiel P, Kingston MJ. Analysis of Airborne Visible-Infrared Imaging Spectrometer (AVIRIS) data of the Iron Hill, Colorado, carbonatite - alkalic igneous rock complex. Economic Geology. 1995;90:1966-1982.
  • Santos ACL, Menezes PTL, Nascimento CTC. Gamaespectrometria aplicada em solo agrícola no noroeste do Estado do Rio de Janeiro. Revista Brasileira de Geofísica. 2008;26(2):181-194
  • Schaefer CE, Candido HG, Corrêa GR, Nunes JA, Arruda DM. Soils Associated with Rupestrian Grasslands. In: Fernandes G(ed.). Ecology and Conservation of Mountaintop grasslands in Brazil. 1ed.: Springer Intern. Publishing; 2016. p.55-69.
  • Schaefer CEGR, Mendonça BAF, Corrêa GR, Ribeiro ASS. Solos desenvolvidos sobre canga ferruginosa no Quadrilátero Ferrífero, Minas Gerais. In: Simpósio Afloramentos Ferruginosos no Quadrilátero Ferrífero, Belo Horizonte: Universidade Federal de Minas Gerais. Instituto de Ciências Biológicas. 2008;1:109-121.
  • Scolforo JR, Carvalho LMT. Mapeamento e Inventário da Flora e dos Reflorestamentos de Minas Gerais. Lavras: UFLA/IEF, 2006.
  • Vasconcellos RM, Metelo MJ, Motta AC, Gomes RD. Geofísica em levantamentos geológicos no Brasil. Rio de Janeiro: CPRM, 1994.
  • Viana PL, Lombardi JA. Florística e caracterização dos Campos Rupestres sobre canga na Serra da Calçada, Minas Gerais, Brasil. Rodriguésia. 2007;58(1):159-177.
  • Vincent RC. Florística, fitossociologia e relações entre a vegetação e o solo em áreas de campos ferruginosos no Quadrilátero Ferrífero, Minas Gerais. [Tese de Doutorado]. São Paulo: Instituto de Biociências da Universidade de São Paulo; 2004.
  • Wilford JR, Bierwirth PN, Craig MA. Application of gamma-ray spectrometry in soil/regolith mapping and geomorphology. AGSO Journal of Australian Geology & Geophysics. 1997;17(2):201-216.
  • Yamaguchi Y, Naito C. Spectral indices for lithologic discrimination and mapping by using the ASTER SWIR bands. International Journal of Remote Sensing. 2003;24(22):4311-4323

Publication Dates

  • Publication in this collection
    21 Feb 2020
  • Date of issue
    2019

History

  • Received
    20 Sept 2018
  • Accepted
    25 Sept 2019
Sociedade de Investigações Florestais Universidade Federal de Viçosa, CEP: 36570-900 - Viçosa - Minas Gerais - Brazil, Tel: (55 31) 3612-3959 - Viçosa - MG - Brazil
E-mail: rarvore@sif.org.br