Ecological impacts of non-native tree species plantations are broad and heterogeneous : a review of Brazilian research

Non-native tree plantations represent 7% of the world’s forests and 1.24% of the Brazilian vegetation. Planted areas are expected to increase in the near future; thus, it is important to systematize existing knowledge on the ecological effects of plantations to aid forest management and biodiversity conservation. Here, we conducted a systematic review of the ecological literature associated with planted Pinus and Eucalyptus species in Brazil. We compared publication metrics with geographical distribution of species, ecosystems, biomes, studied taxa, and ecological impacts. We found 152 publications from 1992 to 2012. Number of publications positively correlated with area planted, number of plantations with forest certification, number of researchers, and richness of studied kingdoms. Most studies were in terrestrial ecosystems (92.1%), the Atlantic Forest biome (55.3%), and the kingdom Animalia (68.2%). Most impacts of non-native tree plantations were negative (55.9%), followed by positive (27%), and mixed (17.1%). Negative impacts were declines in species richness and abundance, seed bank diversity, and natural regeneration. Positive impacts were increase or mainteinance of seed bank diversity and natural regeneration. Mixed impacts were increases in abundance of native tree plantation pests. Taken together, results suggest forest management can help maintain biodiversity if it considers previous environmental conditions and integrates plantations with surrounding habitats.

Large-scale introductions of non-native species are deemed one of the major agents of global environmental changes (Pereira et al. 2012).On the one hand, recent studies have showed that non-native forest plantations (NNFP) can cause the extinction of species (Pereira et al. 2012) by favoring the occurrence of generalist and invasive species (Martin et al. 2012), by simplifying and homogenizing the structure of the environment, and by undermining ecosystem services (Araújo et al. 2010, Simberlof et al. 2010, Rundel et al. MARCOS O. VALDUGA, RAFAEL D. ZENNI and JEAN R.S. VITULE 2014).Also, these plantations may affect water resources, altering its physical, chemical and biological characteristics (Scott 2005, Stenert et al. 2012).On the other hand, some other studies have showed that NNFP can promote the conservation of biodiversity and ecosystem services (Carnus et al. 2006, Brockerhoff et al. 2013) by maintaining native flora (Abreu et al. 2012) and fauna (Carrara et al. 2007).Based on the positive evidence, NNFP have been recommended for conservation of soil and water, for the production of biofuels, and as mitigation agents of climate change (Evans 2009).
In 2005, non-native forest plantations represented 7% of the Earth's forest cover (Evans 2009), and future estimates suggest they might reach 20% by the end of the century (Brockerhoff et al. 2013).The majority of the increase in planted area has occurred mostly in megadiverse countries, such as Brazil, India and China (GIT 2009).In Brazil, massive investments resulted in an increase in the planted area from 400 thousand to 6.53 million hectares (mha) between 1965 and 2012 and are expected to reach nine million hectares in 2020 (ABRAF 2012).Most plantations use species from the genera Pinus (1.64 mha) and Eucalyptus (4.87 mha).NNFP are present in all Brazilian biomes and States, representing 1.24% of the entire Brazilian vegetation cover, and adding great value for the Brazilian economy (ABRAF 2012).In order to mitigate the negative impacts and promote positive impacts of forest plantations, the majority of Brazilian NNFP have one or more forest certifications, such as FSC (Forest Stewardship Council) (BRACELPA 2012).FSC certification requires monitoring of the impacts of plantations on the native biota, proportionally to the scale, intensity and risk management (FSC 2012).The results found in the monitoring should be incorpored into forest management plans and published (FSC 2012).
Brazil is home to ~20% of the world's biodiversity, has the highest number of endemic species globally (Forzza et al. 2012), and is recognized for the ecological (Evans 2009), social, and economic importance of NNFP (ABRAF 2012).In this context, it is important that all aspects of NNFP are consistently evaluated so that positive impacts may be accentuated and negative impacts be avoided or mitigated.Although largescale tree plantations have important implications for management and conservation of natural resources as well as biodiversity, a systematic review of the subject is lacking.Systematic reviews specify search criteria on a particular issue, using predefined protocols, aiming at reducing bias, allowing rigorous replication, and diagnosing gaps (Lowry et al. 2013).Systematic reviews provide high-quality results, identifying, appraising, and synthesizing all relevant studies on a particular topic (Uman 2011, Lowry et al. 2013, Mercuri et al. 2016).
Through a systematic review, this study sought (i) to identify which journals publish more studies involving Brazilian biodiversity and NNFP; (ii) to verify trends in the number of studies published regarding the increase in planted and FSC certified areas; (iii) to assess the number of studies in relation to Brazilian States, to the cultivated species, and to the number of researchers in each state; (iv) to learn which ecosystems and biomes are the most studied in forestry plantation areas; (v) to check which native organisms are the most studied in forestry plantation areas, and its relation to the estimated number of species in the studied kingdoms; (vi) to identify and quantify impacts these plantations are having on biodiversity; and (vii) to qualify and quantify the likely sources of the impacts.

MatErials and MEtHods
Our review was performed using the ISI -Web of Science (WoS).We searched all WoS databases for studies published between 01/01/1960 and 12/18/2012 with the following string of terms: Topic = (pinus OR pine OR eucal?pt*)AND Topic = ("south america*" OR bra?il OR neotropic*).For the first screening, we used the title and summary of the articles to filter out publications not related to the objectives of this study (e.g., studies reporting only abiotic parameters or growth increment).In a subsequent screening, we selected only studies that measured, mentioned, or assessed the ecological effects of NNFP on the biota.Studies outside of that context were also discarded (e.g., publications that reported measures of community structure only in plantation areas, lacking comparisons with data from natural conditions or other comparative control, and studies focusing exclusively on planted non-native species).For all articles that met all selection criteria, we extracted the following information: year of publication, journal, States where the study was performed, species planted, studied ecosystem (aquatic and/or terrestrial), biomes (following IBGE 2013), taxa (i.e., kingdom, phylum, or division of the studied native species) (Whittaker 1969, ITIS 2013, Stevens 2013), objectives, impacts of plantations, and sources of impacts.In the case of studies involving more than one State, ecosystem, biome, kingdom (e.g., Plantae and Animalia), or phylum, data was entered for all groups.If the study indicated more than one likely source of impacts (e.g., increasing richness with aging planting and use of species planted as a resource for native species), these were quantified individually.Thus, the sum of studies across all categories exceeds the number of articles found in the final search.
To identify which journals publish more studies on the influence of NNFP on Brazilian biodiversity in relationship to the total number of articles published in each journal, we used the methodology proposed by Braga et al. (2012); the method considers the amount of studies found in our search for a particular journal and the total number of articles published in the same journal to calculate the relative weight (w) of each journal in publishing studies related to ecological aspects of NNFP.The calculation was performed as follows w = [n/(p × e × y)] × 1000, where n is the total number of articles found in our search for each journal, p is the sum of the number of articles published in the first issue of each year considered in the search , e is the number of issues per year, and y is the number of years considered in the WoS database search.For the analysis, we used only journals that had four or more publications in the search result for the period between 2004 and 2012.The 2-year impact factor of the journal (IF) was also considered in the analysis, and was obtained from the Journal Citations Reports (2013).
To analyze the influence of the increase in NNFP area and areas with FSC certification (FSC 2013) in publication rates, we grouped articles according to the year of publication.Because of the lack of specific data on planted areas for the period between 1992 and 2000, during which there was a decrease in planted area, a proportional reduction of -1.5% per year on planted area since 1990 was applied, following estimates from ABRAF (2004).For the year 2012, we used the growth rate of plantations as in 2011 and added 0.1% on the area planted in relation to 2011, following data from ABRAF (2012).For FSC certification, we evaluated the increase in number of studies between 1996 and 2010, period for which official data of certified areas was available (BRACELPA 2012).
In order to verify if the distribution of the number of publications follows the spatial distribution of tree plantations, we grouped publications by State and cultivated species (e.g., Pinus spp. or Eucalyptus spp.), and then related the information with the total planted area for each genera across all States with official data on NNFP (FIRJAN 2009, ABRAF 2012).Grouping data by State was required because Pinus spp.plantations are concentrated in the States of the southern region whereas Eucalyptus spp.plantations are concentrated in other States (ABRAF 2012).We MARCOS O. VALDUGA, RAFAEL D. ZENNI and JEAN R.S. VITULE also evaluated the spatial dimension reached by the studies by calculating the ratio of hectares per study for the total planted area, the total area with FSC certification, and the areas planted with Pinus spp.and Eucalyptus spp.
To test the relationship between the number of studies found per State and the number of researchers based on them in 2012, we used data on number of researchers associated with the knowledge areas encompassed by this study, excluding unrelated areas (e.g., medicine and literature).Researcher data were provided by the Brazilian Coordination for the Improvement of Personal Higher Education (CAPES 2013).
To evaluate the proportion of studies on ecosystems (aquatic or terrestrial), biomes (IBGE 2013), and native organisms, we grouped native organisms into taxonomic groups following Whittaker (1969).The greatest taxonomic refine ment was used a posteriori for the two most studied kingdoms, following the Integrated Taxonomic Information System (ITIS 2013), and the Angiosperm Phylogeny Website (Stevens 2013).To estimate the relationship between the numbers of studies found in each kingdom and the number of species in those realms, we used the data from Lewinsohn and Prado (2005).
We also evaluated the proportion of each type of impact that NNFP had on studied ecosystems, biomes, and kingdoms of native species.We considered impact as any change (i.e., fluctuations different than natural variations) in the structural and functional parameters of populations and communities of native taxa and ecosystem functions or services: richness, abundance, reproduction, growth, colonization, density or biomass of organisms, diversity of the seed bank or natural regeneration, biomass or microbial activity in the soil and leaf decomposition (ISO14050 2013).Thus, the impact could have a positive effect (e.g., increase or maintenance of the parameters listed above), a negative effect (e.g., decrease in any of the same parameters listed above) or a mixed effect (positive for one or some of the parameters and negative for one or some of the parameters).To qualify the impacts found and facilitate the understanding of the variations generated, we classified them according to their specificity.Impacts on reproduction, growth, colonization, density, and biomass of organisms were grouped as "population-level impacts."Impacts on the richness, relative abundance, diversity of the seed bank, and natural regeneration were grouped as "community-level impacts."Changes in biomass, mesofauna, microbial activity in soil, decomposition rates in aquatic ecosystems, biological pest control by native organisms, and leaf decomposition were grouped as "ecosystemlevel impacts."Then, based on the conclusion of the studies, we obtained the likely sources for the reported impacts, and quantified the respective proportions of each impact group.

DATA ANALySIS
We conducted graphical analysis of correlation and Pearson correlation tests to exam the relationship between the IF and the number of published studies that matched our search criteria.The same tests were used to investigate the relationship between the number of studies and variation of spatial, temporal and interspecific areas with NNFP, certified by FSC, as well as between the numbers of researchers based on each State.Finally, we checked whether there was a relationship between the number of studies found in the analyzed kingdoms and the estimated number of species for each kingdom.For all correlations performed, data were tested for normality using the Shapiro-Wilk test and log-transformed where appropriate.We performed a Chisquare test to assess differences in the distribution of impacts between ecosystems, biomes and kingdoms.All analyzes were performed in the Origin Pro-8 software (OriginLab Corporation).

rEsults
The initial search resulted in 1,494 articles from which 330 passed the first screening.A total of 152 (10.1%) articles met the inclusion criteria during the second and final screening.These studies were published in 66 journals, 8.5% of them being published in "Forest Ecology and Management," 5.2% published in the "Brazilian Journal of Soil Science," and 4.6% published in "Biodiversity and Conservation" (Fig. 1).The 63 remaining journals published less than 4% of the studies each.The highest relative weights (w) were 1.16 for the "Brazilian Journal of Ornithology," 0.81 for the journal "Mammalia," and 0.69 for the "Journal of Tropical Ecology" (Fig. 1).The 63 remaining journals had w < 0.69.The mean ± standard deviation of the IF was 1.336 ±1.07.The journal with higher IF was the "Journal of Applied Ecology" (IF = 4.7, n = 2), followed by "Conservation Biology" (IF = 4.355, n = 1), and "Perspectives in Plant Ecology, Evolution and Systematics" (IF = 4.158, n = 1).There was no relationship between IF and relative frequency of studies per journal.
All studies matching our criteria were published after 1992, and we detected a substantial increase in the number of publications after 2005 (Fig. 2).The seven-year period concentrated 82.2% of the studies found, correlating with the increase in NNFP area (r = 0.5, p = 0.01).Also, a positive correlation was found between the number of studies published and the increase in FSC certified area (r = 0.7, p < 0.001, Fig. 3a), with a proportion of one study for every 37,000 ha of certified area.
Out of the 17 Brazilian States with NNFP, 14 had studies performed in their territories.Three States (Minas Gerais, São Paulo, and Rio Grande do Sul) held the majority (65.7%) of the studies (Table I); 75% of the studies were conducted in areas planted with Eucalyptus spp.and covered 14 out of the 17 States that plant these species; 66.3% of the studies on Eucalyptus spp.were conducted in Minas Gerais, São Paulo, and Rio Grande do Sul, with the proportion of one study for every 46,800 ha of planted area (Table I).The remaining studies (25%) were conducted in areas planted with Pinus spp.and in eight out of the 12 States growing pine tree plantations, with the majority of studies (75%) being located in Rio Grande do Sul, São Paulo, and Paraná with a proportion of one study for every 28,500 ha of planted area (Table I).The Pearson correlation test showed a positive relationship between number of publications and area planted per State (r = 0.7, p < 0.001), and between number of studies and area planted with Eucalyptus spp.(r = 0.6, p = 0.009) (Fig. 3b) and Pinus spp.(r = 0.8, p = 0.003), repectively (Fig. 3c).We also found a positive relationship between number of studies published and number of researchers based on each State (r = 0.5, p = 0.02) (Table I, Fig. 3d).Minas Gerais, São Paulo, and Rio Grande do Sul were the States with both the greatest production of studies and number of researchers (Fig. 3d).The vast majority of the studies (92.1%) were performed in terrestrial ecosystems, and only 7.9% were carried out in aquatic ecosystems (Table II).The Atlantic Forest was the most studied biome (55.3% of the studies), followed by the Cerrado (28.3%), the Amazon (7.2%), the Pampas (6.6%), and the Caatinga (2.6%).Among the kingdoms, Animalia was the focus of 68.2% of the studies, followed by Plantae (18.2%),Monera (7.8%), and Fungi (5.8%) (Table II).The Pearson correlation test showed a strong relationship between number of studies published and the estimated number of species in each kingdom (r = 0.9, p = 0.01).In the kingdom Animalia, the phylum Arthropoda was the most studied (56.5% of the studies), followed by the phyla Chordata (33%), Mollusca (4.3%), Annelida (3.5%), Platyhelmintes (1.8%) and Cnidaria (0.9%).In the kingdom Plantae, the Angiospermae division was the most studied (57.1% of the studies), followed by the Gimnospermae (20%), Pteridophyta (17.1%), and Bryophyta (5.8%) divisions.

taBlE i number of studies in each
We found a prevalence of negative effects of NNFP on the communitylevel impacts in different ecosystems, biomes and native organisms studied, with the exception of the kingdom Plantae (Table II).There was an overall high frequency of results showing negative impacts among all the studies (n = 152, χ² = 17.1, df = 2, p < 0.001), as well as per some individual groups: terrestrial ecosystems (n = 140, χ² = 27.1, df = 2, p < 0.001), aquatic ecosystems, (n = 12, χ² = 10.5, df = 2, p = 0.005), Atlantic Rain Forest (n = 84, χ² = 16.9, df = 2, p < 0.001), and Pampa (n = 10, χ² = 8, df = 2, p = 0.02) biomes, and Animalia kingdom (n = 105, χ² = 13.2, df = 2, p < 0.001).Of the studies reporting negative impacts, 65.9% (n = 56) cite some type of suboptimal forest management practices and the resulting reduced structural heterogeneity of the environment as sources of these impacts.Similarly, 24.7% (n = 21) of the studies reported reduced rates of reincorporation of nutrients in the organic matter by mesofauna and microbiota decomposition (nutrient cycling), as well as decreases in the qualities of soil, water, and litter biomass.Another 4.7% (n = 4) of the studies indicated a decrease in species richness of the seed bank along with the development of the rotation cycle.Finally, 4.7% of the studies (n = 4) reported invasion of surrounding areas of native vegetation by the planted species, the increased frequency of fungi diseases in native species, occurrence of hunting, and increased frequency of road kills and fires as generators of negative impacts.
Among studies reporting positive impacts, 53.4% (n = 22) mentioned maintenance of native vegetation in the understory of the plantations and connectivity among native remnants.For 27.6% (n = 11) of the studies, the colonizing behavior of the cultivated species was a source of positive impact, along with maintenance of the seed bank, and 19% (n = 8) of the studies reported the improvement of soil quality and increased fire protection.
In studies showing mixed impacts, 44.9% (n = 12) cited as possible sources the maintenance of the carrying capacity of populations.6.9% (n = 2) pointed to an increase in richness of herbaceous species and a decrease in richness of climax species with aging plantations.3.4% (n = 1) mentioned the use of cultivated areas by native species (i.e., foraging and shelter).Another 3.4% (n = 1) pointed to the role of native species in controling pests that attack NNFP.Finally, 41.4% (n = 10) of the studies did not define the sources of impacts.

discussion
For over a century Brazil has been cultivating NNFP (Shimizu 2006).However, the interest in understanding the influence these plantations have on the biota is recent, even though it has grown substantially in the past decade, as shown by our results.The increased interest in understanding the impacts caused by the large-scale introduction of non-native species on ecosystems and loss of biodiversity (Schlaepfer et al. 2011, Valéry et al. 2013, Vitule et al. 2012) might be the cause for the recent growth in studies in cultivated areas.The diversity and weight of the journals where publications on the studied subject were found show that knowledge of possible implications of NNFP on the biota can be favored by transdisciplinarity work (Crow 2010).The high frequency of publications in journals with low IF may decrease the influence of the studies on other research and on scientific thought (Pyšek et al. 2006, González-Alcaide et al. 2012), as well as potentially limit their dissemination to other sectors of society (Primack et al. 2009), which may impair managerial decision making.
The adoption of FSC certification appears to have had a positive influence in the number of studies per planted area.Our results show an increase in research intensity after the wide adoption of certification by forest companies, probably because of the obligation imposed by FSC for companies to monitor impacts of plantations (FSC 2012).Still, many studies are produced and kept within the reforestation companies, making it difficult for external researchers to access the data ( Van-Kujik et al. 2009, Visseren-Hamakers andPattberg 2013).The availability of nearby research centers, such as universities, can increase the number of studies in NNFP areas.According to our review, the number of researchers per State correlated with the number of studies published.Formal interactions between research centers and reforestation companies are of great importance for the production and dissemination of knowledge, besides the fact that the interaction results in reduced costs for research (Crow 2010, Laursen et al. 2011).We suggest that research produced by the companies to comply with forest certification standards should be published in indexed journals, thus increasing the knowledge on the influence of NNFP and the FSC on biodiversity conservation (Visseren-Hamakers and Pattberg 2013).
The terrestrial ecosystem has been the most studied, as in other studies with non-native species (e.g., Lowry et al. 2013).Factors that may have stimulated the development of studies in terrestrial ecosystems rather than aquatic ones are the large affected area, the direct impact of NNFP on the ecosystem, the greater ease of study (Lowry et al. 2013), the presence of flagship (Heywood 1996) or umbrella species (Wilcox 1984) that attract greater investment for research projects, or even the presence of species of economic interest (e.g., Ramos et al. 2008).Individually or combined, these factors may also have encouraged the highest concentration of studies in the Animalia and Plantae kingdoms, as our results show.Moreover, we expected the number of studies involving the aquatic ecosystem would be greater than what was recorded because NNFP affects hydric resources (Scott 2005) and cover a large number of watersheds (Price 2011).Furthermore, freshwater systems can be more sensitive to the influence of nonnative species compared to marine or terrestrial systems (Cox and Lima 2006).The effects of NNFP on aquatic ecosystems must receive more attention if we are to fully understand the impacts and prevent biodiversity and ecosystem services losses like water supply.
The impacts described by the studies suggest that the effects of NNFP on biodiversity are MARCOS O. VALDUGA, RAFAEL D. ZENNI and JEAN R.S. VITULE contingent on the organism and the conditions within the plantations.Often, the responses of native organisms may depend on various factors that can alter the magnitude and characteristics of the impacts found, such as proximity to areas with native vegetation, protection from hunting, connectivity between patches of native vegetation (Mazzolli 2010), pattern or order of arrival of invasive species in the ecosystem (Simberlof and Vitule 2013, Dickie et al. 2014) or propagule pressure (Rundel et al. 2014).Negative impacts predominate when NNFP replace native vegetation, whereas positive impacts occur when NNFP replace degraded areas; a greater number of mixed impacts are found when NNFP replace agricultural areas (Meyfroidt and Lambin 2011).
The high proportion of negative populationand community-level impacts associated with NNFP is quite worrying.The effects of land use changes can decrease carrying capacity, genetic variability, and fitness of native populations (Reed 2005, Willi et al. 2006); while also increasing the number of species threatened with extinction and causing losses of ecosystem services (e.g., litter fall and leaf decomposition).The environmental simplification in NNFP was appointed as the main of negative impact (e.g.Lyra-Jorge et al. 2010, Martin et al. 2012).Generally, environmental simplification was related to empoverishment of the seed bank (Nobrega et al. 2009) and woody plant diversity (Scolari et al. 2010), resulting in an increase in abundance of generalist organisms (Gonçalves et al. 2008, Gheler-Costa et al. 2012) and native species considered pests in plantations (Zanuncio et al. 1998).Furthermore, the studies reported that biological invasions by the planted species were one of the sources of negative impacts (e.g., Rolon et al. 2011).Some biomes (e.g., Cerrado and Pampa) may be highly susceptible to invasions by some of the most common tree plantations (Abreu and Durigan 2011).The invasive potential can be related to a larger size of planted species in relationship to native vegetation (Pysĕk et al. 2014), to adaptations possessed by some tree species (Zenni et al. 2014), and to ecological and physiological characteristics of planted species (Zenni et al. 2013).For instance, P. eliotti Engelm.dominated and reduced the richness and abundance of native Cerrado vegetation after 22 years of establishment of commercial plantations of this species (Abreu and Durigan 2011).In coastal ponds in the Pampa biome, invasion of riparian vegetation by cultivated species increased the abundance of generalist macroinvertebrates (Stenert et al. 2012).Thus, invasions by forestry trees can jeopardize the stability of ecosystems, affect primary and secondary productivity, and biogeochemical cycles (Simberlof et al. 2010, Rundel et al. 2014, van Wilgen and Richardson 2014).
Despite the prevalence of negative impacts, positive interactions may play an important role in regulating ecosystems (Halpern et al. 2007) and the predominance of positive impacts in the kingdom Plantae should be valued.The positive effects we reviewed can be related to the use of previously degraded areas for the implementation of NNFP (ABRAF 2012).The ability of some non-native tree crops to colonize degraded areas may aid the restoration of these areas, increasing the diversity and abundance of secondary and climax species (Nobrega et al. 2008).Moreover, the legal requirements that ensure the presence of native vegetation between cultivated areas (Federal Brazilian Law 12,651 2012) increase the flow of native propagules, vegetation and the maintenance of seed bank within the plantations (Abreu et al. 2012).

conclusion
In megadiverse countries such as Brazil, the establishment of new NNFP should be limited to degraded rather than more pristine natural areas.There is a need for greater efforts to improve our understanding of the ecological implications of

figure 1 -
figure 1 -Journals with four or more publications with ecological studies in NNFP.The grey bars represent the number of studies and the black squares represent the relative weight (w) of each journal.The numbers on the bars represent the number of studies in each journal.
Brazilian state (n).Where: EPa = area planted with Eucalyptus spp.(ha) (aBraf 2012); nEs = number of studies found for Eucalyptus spp.plantations; sPE = ratio of area planted with Eucalyptus spp.per study; PPa = area planted with Pinus spp.(ha) (aBraf 2012); nPs = number of studies found for Pinus spp.plantations; sPP = ratio of area planted with Pinus spp.per study; rEs = number of researches (caPEs 2013 figure 2 -Relative frequency of studies published and area planted with non-native tree species per year between 1992 and 2012.The grey bars represent the number of studies and the black squares represent the area planted with non-native tree species (million hectares) (ABRAF 2012).The numbers on the bars represent the number of studies in each year.

figure 3 -
figure 3 -Correlations (a) between number of studies and increase in number of FSC certified areas, (b) between number of studies and area planted with Eucalyptus spp. in each Brazilian State, (c) between number of studies and area planted with Pinus spp. each Brazilian State, and (d) between number of studies and number of researchers in each Brazilian State.Dotted lines indicate the lines of best fit.AP = Amapá; BA = Bahia; ES = Espírito Santo; GO = Goiás; MA = Maranhão; MG = Minas Gerais; MS = Mato Grosso do Sul; MT = Mato Grosso; PA = Pará; PE = Pernambuco; PI = Piauí; PR = Paraná; RJ = Rio de Janeiro; RS = Rio Grande do Sul; SC = Santa Catarina; SP = São Paulo; TO = Tocantins.