SciELO - Scientific Electronic Library Online

vol.26 issue3Meiotic Behavior and Pollen Viability of Spondias mombin L.: Native Fruit Species of the AmazonAnthropogenic Disturbances Affect the Relationship Between Spectral Indices and the Biometric Variables of Brazilian Savannas author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand




Related links


Floresta e Ambiente

Print version ISSN 1415-0980On-line version ISSN 2179-8087

Floresta Ambient. vol.26 no.3 Seropédica  2019  Epub July 15, 2019 

Original Article

Conservation of Nature

Climatic Response of Cedrela fissilis Radial Growth in the Ombrophilous Mixed Forest, Paraná, Brazil

Amanda Koche Marcon1

Tomaz Longhi-Santos2

Franklin Galvão2

Kelly Geronazzo Martins3 

Paulo Cesar Botosso4 

Christopher Thomas Blum2 

1Programa de Pós-Graduação em Engenharia Florestal, Departamento de Ciências Florestais, Universidade Federal do Paraná – UFPR, Curitiba/PR, Brasil

2Departamento de Ciências Florestais, Universidade Federal do Paraná – UFPR, Curitiba/PR, Brasil

3Departamento de Engenharia Ambiental, Universidade Estadual do Centro Oeste – UNICENTRO, Campus de Irati, Irati/PR, Brasil

4Embrapa Florestas, Empresa Brasileira de Pesquisa Agropecuária – EMBRAPA, Colombo/PR, Brasil


We examined the climatic response of Cedrela fissilis radial growth based on precipitation, air relative humidity, temperature and monthly thermal amplitude. It was intended to assess how the secondary growth of C. fissilis is influenced by the climatic variables and which one are the best growth predictors in Ombrophilous Mixed Forest, southern Brazil. Wood cores were processed using classical dendrochronology methodology. Principal Components Analysis, Generalized Linear Models and correlation were used to explore the relationship between radial growth and climate. Our results indicated that the best radial growth predictor is the temperature: the seasonal behavior, under well-watered conditions, seems to have a dominant effect on growth responses. Changes in tree growth corresponding to an increase in temperatures suggest a sensitivity of the species to climate changes. These results are important to help understand how the global warming may influence long-lived pioneer tree growth.

Keywords:  dendrochronology; climate changes; dendroecology; Araucaria forest


Advances in understanding plant autecology across environmental gradients are directly related to the use of environmental variables linked to their growth. The particularities of each species and the environmental conditions may determine occurrence, growth patterns and competitive interactions among species (Tilman, 1985; Pausas & Austin, 2001).

As plants are sessile organisms, their development is profoundly influenced by the environment and its disturbances (Cui et al., 2014). The development patterns are characterized by endogenous rhythms, sometimes driven by environmental variation or others exogenous factors (Walter et al., 2009), where a single factor may limit growth, reproduction or distribution of a species (Billings, 1952). This understanding is also important for the analysis of the possible effects of climate changes (Cusatis et al., 2013).

Advances in the science of climate changes are providing a clearer understanding of the inherent variability of Earth’s climate and its likely response to human and natural influences (Moss et al., 2010). The increase of global mean surface temperature by the end of the 21st century range from 3.7 °C to 4.8 °C, impacting on natural and human systems on all continents and oceans, according to Intergovernmental Panel on Climate Change (IPCC, 2014). From this knowledge, studies have shown that some modification in tree development triggered by drought and/or high temperatures may already be occurring in response to global climate change (Allen et al., 2010). The response of tree growth to a change in temperature may differ from different functional groups and biomes: in high latitude and altitude, for example, tree growth may be temperature-limited and thus benefit from some degree of warming, as opposed to warm-adapted species (Way & Oren, 2010).

Tree-ring chronologies are often used to detect shared signs in a population and reconstruct histories of unmeasured events in the past (Fritts, 1976; Gholami et al., 2015), displaying an important tool to the knowledge on the relationship between climatic conditions and tree growth (Cedro et al., 2013; Gebrekirstos et al., 2014). Thus, tree-ring studies might add important information about past climate variability and may potentially contribute to a better understanding of future climate patterns (Li et al., 2012).

The dendrochronological potential of Cedrela fissilis Vell. (Meliaceae) is evidenced by its distinct and conspicuous annual growth layers, due to its semi-porous rings, visible to the naked eye (López & Villalba, 2016). These rings are delimited by marginal bands of axial parenchyma and earlywood vessels with larger diameters (López & Villalba, 2016; Andreacci et al., 2017). Research this species climatic responses is relevant, considering that it has a wide geographic distribution and allows to compare responses in different climatic types. C. fissilis occurs from Colombia to Brazil, where can be found in the phytogeographic domains of Amazon Rainforest, Cerrado vegetation (Brazilian savanna) and Atlantic Rainforest, occurring in the Mixed and Dense Ombrophilous Forests, Deciduous and Semideciduous Seasonal Forests, as well as in Brazilian savanna formations (Muellner et al., 2010; Flores, 2018).

Cedrela fissilis can be classified as a longer-lived pioneer species, mainly due to its high light demand, height, longevity and denser wood (Corlett, 1995; Lima et al., 2009), although not without controversy (e.g. Gris et al., 2012; Konig Brun et al., 2017). It is characterized by its deciduousness, a morphological adaptation that is tied to the origin of the genus Cedrela at forests with seasonal climate (Muellner et al., 2010).

In this context, the aim of this study was to assess the responses of C. fissilis radial growth to climate in the South Central region of the state of Paraná, Brazil. We aimed to answer the following questions: i) How does climate influence C. fissilis tree growth in this region? ii) Which climatic variables are the best annual growth predictors for this species in this region? In this way, our central hypotheses were that: radial growth of C. fissilis shows a relevant sensitivity to changes in the local climate and temperature is a limiting factor for radial growth in the South Central region of the state of Paraná, Brazil.


2.1. Survey area

The study was carried out in the municipality of Pinhão (PR) in area belonging to COPEL - Companhia Paranaense de Energia (the electric company of the state of Paraná), coordinates 26º 00’ 23” S 51º 40’ 06” W, located in the Terceiro Planalto Paranaense (Paraná`s Third Plateau). The climate is classified as Cfb, with well-distributed rainfall, mild summers and occurrence of frost (Figure 1) (Miranda, 2009; Alvares et al., 2013). The average annual precipitation for the region in the surveyed period is 1,745 mm. The average temperature of the coldest month is below 14 °C and the average temperature in the hottest month does not reach 23 °C (KNMI, 2017; Brasil, 2017).

Figure 1 Climate diagram for study area, in the Ombrophilous Mixed Forest, Pinhão, Paraná, Brazil. Each bar represents the cumulative monthly precipitation (Jan: January; Feb: February etc) from 1945 to 2015 (data source: Brasil, 2017). Solid line represents the average temperature, dashed line represents maximum temperature (Max. temperature) and dotted line represents minimum temperature (Min. temperature), from 1901 to 2014 (data source: KNMI, 2017). 

The local altitude is quite variable, with a gradient of 700 to 1,250 m.a.l.s. It has a hilly relief, with steep slopes (Miranda, 2009). Geologically, the region is part of the São Bento geological group (Grupo São Bento), and the predominant formation is the Serra Geral Formation, which resulted in acid and base rocks, usually basalt intercalated with sedimentary materials (Frank et al., 2009; Chahud & Petri, 2010; Lima, 2012). Local soils are predominantly Latosols, Nitisols and Cambisols (Miranda, 2009).

Vegetation is characterized as Ombrophilous Mixed Forest, although some characteristics of semideciduous seasonal forests occur in lower terrains (Miranda, 2009).

2.2. Data collection

In April 2016, wood samples from 11 (eleven) Cedrela fissilis Vell. (Meliaceae) adult individuals (according to tree size and presence of reproductive structures) were collected through a non-destructive method – 2 to 3 diametrically opposite increment cores per tree using a 5 cm increment borer at 1.3 m above ground level (DBH) (Fritts, 1976; Elling et al., 2009). Trees with a diameter greater than or equal to 40 cm were selected (Figure 2), with similar site conditions in a secondary forest formation with an approximate area of 60 ha.

Figure 2 Adult Cedrela fissilis tree in the Ombrophilous Mixed Forest, in study site (a) and macroscopic image of the annual rings of Cedrela fissilis showing typical tree rings with conspicuous marginal parenchyma bands and larger earlywood vessels (arrows) (b). Scale bar = 1 cm. 

Samples were dried at room temperature and set on wooden holders for the surface polishing process of their cross sections, using successively finer grit sizes (80-1,000 grits/cm2). Tree-ring delimitation was made using a stereo microscope, and samples were subsequently digitized (1200 dpi resolution) using an HP Scanjet G4050 scanner. An image processing software was used to measure the tree-ring width.

To assess the climatic response of radial growth of these trees, historical precipitation data from the Jangada do Sul weather station (Copel) were used, which covered a period from 1945 to 2015. Historical temperature data comes from the KNMI Climate Explorer database and covers the period from 1901 to 2014 (KNMI, 2017). Air relative humidity data were obtained from the National Oceanic & Atmospheric Administration NOAA/ESRL Physical Sciences Division database (Kalnay et al., 1996) for the 1948 to 2015 period.

2.3. Data analysis

Tree ring crossdating was performed using the COFECHA® software (Holmes, 1983) with a 30-year moving window and a 15-year overlap, allowing to check the cross-dating quality, which calculates correlation coefficients between individual tree-ring series as a way to identify absent or false rings (Venegas-González et al., 2018). The standardization of the radial time series, which was done in order to remove biological growth tendencies and to establish the chronologies, was performed by adjusting original tree ring measurements through the application of a cubic smoothing spline function with a 40% frequency cut-off using the ARSTAN® software (Cook & Holmes, 1986). The sample quality was verified based on Expressed Population Signal (EPS) and Rbar. EPS is a measure of the common signal present in a chronology compared to an infinitely replicated chronology and Rbar represents the average correlation coefficient between each series and tree-ring chronology (Wigley et al., 1984).

To assess climatic data, a Principal Component Analysis (PCA) with varimax orthogonal rotation was performed, which included the variables minimum temperature, average temperature, maximum temperature, monthly thermal amplitude (difference between monthly maximum and minimum temperatures), precipitation and air relative humidity. The Kaiser test was used for the determination of the number of significant components (Kaiser, 1958). Additionally, we tested the significance of the variables with PCA main components using the bootstrap randomization method (1000 randomizations).

Generalized Linear Models (GLM) were adjusted using PCA components as predictors and chronology as a response variable, to verify the relationship between the components of the PCA and the chronology. Climatic and tree-ring data of the time period corresponding to the age of the youngest tree examined in this study were used to adjust the GLM. When necessary, we transformed the data. In this way, all climatic variables and trees were evaluated for the same period. The models were compared based on the Akaike Information Criterion (AIC) (Akaike, 1974). In addition, the scores of the PCA components were compared with tree-ring chronology using the Spearman correlation coefficient (rs).

Statistical analyzes were developed in the R software (R Development Core Team, 2017).


The ring-width chronology of C. fissilis covers the period from 1891 to 2014 (Table 1, Figure 2). Of the 11 trees and 24 time series initially sampled, it was possible to crossdate 9 trees and 15 time series. During the tree ring crossdating of C. fissilis, some samples did not result in satisfactory correlations with the master series and were excluded from the data set. Some characteristics that made identification more difficult as: wood anatomical anomalies, inconspicuous rings and false or local missing rings.

Table 1 Descriptive dendrochronological statistics for the Cedrela fissilis tree-ring chronology in the Ombrophilous Mixed Forest, Paraná, Brazil. Numbers in parentheses mean the number of series or trees sampled. Master series represents the period covered by the wood samples. Series intercorrelation measures the association between series (Fritts, 1976) and critical correlation suggest the limit value for a crossdating with 30-year moving window and a 15-year overlap. Mean sensitivity is the average difference between successive ring index values. Mean sensitivity above 0.3 indicate the presence of considerable high-frequency variance (Fritts, 1976). 

Description Chronology
Master series 1891-2014
Number of series dated 15 (24)
Number of trees 9 (11)
Ages (years) 55-124
Average length of time series (years) 86.6
Series intercorrelation 0.313
Critical correlation (p < 0.01) 0.423
Mean sensitivity 0.535

Higher growth was observed in the early years of the chronology (1891-1918, Figure 3). However, only a few series make up the tree-ring chronology in this period, indicating that the higher growth may not be representative.

Figure 3 The mean tree-ring width (dashed black line) and individual growth ring width (solid grey lines) (a); Rbar (dashed black line with circles) and Expressed Population Signal (EPS, solid black line with triangles) statistics (b) and the standard tree-ring chronology (solid black line) with number of series per year (grey area) (c) for Cedrela fissilis in the Ombrophilous Mixed Forest, Paraná, Brazil, for the reference period 1891–2014. The standard tree-ring chronology (tree-ring index) was built using a cubic smoothing spline function with a 40% frequency cut-off. Rbar represents the average correlation coefficient between each series and the chronology and EPS is a measure of the common signal present in a chronology compared to an infinitely replicated chronology. 

Three components were selected based on the analysis of the main components, explained 90.01% of data variability (total of explained variance after sum of eigenvalues ≥90%, Figure 4, Table 2). Correlations (loadings) of the climatic variables with PCA components allowed the separation of data in sets of orthogonal latent variables: PCA1, formed mainly by the minimum, average and maximum temperatures, PCA2, formed basically by precipitation and air relative humidity, and PCA3, formed by monthly thermal amplitude.

Figure 4 Diagram of ordination in the two first axes yielded by Principal Component Analysis (with varimax rotation) of climatic variables, given as vectors: “ThermAmplitude” = monthly thermal amplitude; “AirHumid” = relative air humidity; “Precipitation” = precipitation (rainfall); “MaxTemp” = maximum temperature; “AverageTemp” = average temperature; “MinTemp” = minimum temperature. 

Table 2 Principal Component Analysis (with varimax rotation) statistics. Eigenvalues: quantify the amount of inertia in each component (Jombart et al., 2009); EV (%): explained variance; Variables: show the composition of the latent variables (PCA main components); Loadings: indicate the correlation of the variable with the main component (*95% probability level; **99% probability level, according to bootstrap randomization method). 

Eigenvalues EV
Variables Loadings
PCA1 3.02 50.0 Minimum temperature 0.92*
Maximum temperature 0.91*
Average temperature 0.98*
PCA2 1.44 24.1 Precipitation -0.85**
Air relative humidity -0.75**
PCA3 0.94 15.6 Monthly thermal amplitude 0.92*

The Generalized Linear Models indicated that the best predictor to describe the tree-ring chronology includes, based on the Akaike Information Criterion, the first component of the PCA (PCA1), since it presents the lowest AIC value (Table 3).

Table 3 Models examining effects of latent variables on tree-ring chronology using the Akaike Information Criterium (AIC). Latent variables: PCA1 = component comprising the minimum, maximum and average temperatures; PCA2 = component comprising precipitation and air relative humidity; PCA3 = component comprising the monthly thermal amplitude. Df = degrees of freedom. 

Models Df AIC
PCA1 1 -16.790
PCA2 1 -16.189
PCA3 1 -16.171
PCA1xPCA3 2 -14.975
PCA1xPCA2 2 -14.791
PCA2xPCA3 2 -14.217
PCA1xPCA2xPCA3 3 -13.005

By comparing tree-ring chronology behavior with the PCA latent variables, several points of synchrony among series can be observed, where years with higher scores correspond to higher growth rates (Figure 5). The Spearman correlation between PCA1 and the tree-ring chronology for the period corresponding to the age of the youngest tree collected (1960-2014) was not significant (rs = 0.14; p = 31). However, correlation was significant for the 1960-1979 time period (rs = 0.65; p = 0.005). For the last time period (1980-2014), the correlation with PCA1 was not significant (rs = 0.05; p = 0.76). The other latent variables showed no correlation to the tree-ring chronology.

Figure 5 Tree-ring chronology of Cedrela fissilis and latent variables scores constructed from a Principal Component Analysis for the reference period 1960–2014. Latent variables: PCA1 = component comprising the minimum, maximum and average temperatures; PCA2 = component comprising precipitation and air relative humidity; PCA3 = component comprising the monthly thermal amplitude. Solid black line represents the tree-ring chronology (tree-ring index), solid grey line represents PCA1, dotted grey line represents PCA2 and dashed grey line represents PCA3. 


Our results showed that C. fissilis perform a relevant sensitivity to changes in the local climate, endorsing the possibility of using growth rings as proxies to detect disturbances in the climate scenario (Fan et al., 2009; Venegas-González et al. 2018). The mean sensitivity for the whole set of series is considered high (> 0.3) (Fritts, 1976). Although the time series intercorrelation did not reach critical correlation, this result indicates the high sensitivity of the species to record climate changes (Fritts, 1976).

The EPS, another important parameter to verify the quality of the chronology, reinforces the reliability of the data, especially after 1935, where its value is higher than 0.75. Although it is suggested by Wigley et al. (1984) that the EPS is adequate if higher than 0.85, this parameter can be adjusted depending on the ecological characteristics of the species analyzed (Piraino et al., 2015). Difficulties in crossdating process of C. fissilis, like wood anatomical anomalies, inconspicuous rings and false or local missing rings, were also found by other authors that studied this species (Andreacci et al., 2014; Paredes-Villanueva et al., 2016; López & Villalba, 2016; Cusatis et al., 2013).

PCA main components evidenced the connections between climatic data. The climatic variables formed three main components, which can be seen as latent variables that express most of the information in the data (Wold et al., 1987). This pattern was expected, since it grouped variables correlated with each other. Minimum, average and maximum temperatures are, in essence, collinear, because they originate from the same meteorological variable. Precipitation and air relative humidity are also related to each other, considering that higher volumes of precipitation are related to a more humid atmosphere. Precipitation occurs mainly due to air relative humidity, which can be generated locally, by evapotranspiration or transported from other areas (Quadro et al., 2012). Finally, PCA3 comprises monthly thermal amplitude, the temperature variability - an important variable, since exhibited the role of temperature fluctuations (Folguera et al., 2009).

By graphically observing the behavior of chronology and temperatures, it is possible to distinguish a pattern where higher temperatures contribute to an improved tree growth. Other studies have found temperature as the main growth predictor of several species from Ombrophilous Forests in the South region of Brazil (eg Machado et al., 2010; Kanieski et al., 2012, 2013; Cusatis et al., 2013; Andreacci et al., 2014), where climate variations are related to the seasonality of the temperatures. The South region of Brazil is located almost entirely in the subtropics. Even though it does not have a dry season, the climate presents an important thermal oscillation throughout the year, with a stronger distinction between winter and summer (Rosemback et al., 2010). For this reason, secondary growth occurs mainly in spring and summer, with temperatures responsible for the cambium activate and xylem differentiation (Begum et al., 2013).

Cedrela is a genus of deciduous behavior occurring in Ombrophilous and Seasonal Forests, as well as in Brazilian Cerrado and Caatinga (Muellner et al., 2010; Flores, 2018). In regions where there are droughts, there is a clear seasonality of the growth determined by leaf abscission and cambial dormancy in the driest period (Costa et al., 2013). In regions where temperatures are seasonal, the cold period is decisive (Jeong et al., 2011). In addition, in mid-latitudes seasonal variations in temperature are associated with variations in the photoperiod, with longer days. The light regime is a critical component in determining the biomass production of photoautotrophic beings, which use light as their source of energy (Wahidin et al., 2013). Thus, it is possible that in tropical and subtropical areas, where precipitation is relatively well-distributed throughout the year, growth increase occurs all year long and temperature and photoperiod may promote a stimulus to cambial activity for some species (Santarosa et al., 2007; Shimamoto et al., 2016).

The adjusted model did not elect the component comprised of precipitation and air relative humidity as a significant for the prediction of C. fissilis growth. The fact that there was no well-defined dry season probably contributed to the absence of large variations in air relative humidity and precipitation, not incurring in significant growth responses. Nevertheless, when inferring about growth responses to annual precipitation, there is a tendency that years with stronger rainfall and higher growth match, while fewer rainfall corresponds to a reduction in growth. It can be observed that even in places where precipitation is not seasonal - thus not generating significant responses in the growth prediction model - this climatic variable match with some of the points in the chronology. Increased cambial activity depends on tree water status during periods of high metabolic activity (Zweifel et al., 2006), which could justify the correlation between precipitation and temperature rates and the width of the rings.

Andreacci et al. (2014) evaluated climatic signs in C. fissilis tree rings also in Ombrophilous Forests and suggests that water supply and high temperatures in the growing season have consequences in tree-ring width. Bud break, flowering and cambial activity are events that entail cellular expansion, a process inhibited by the lack of water, even just a moderate one (Borchert, 1994). These events can be even more affected in deciduous species, since one of the basic requirements for metabolic activities is water availability, for the hydrolysis of macromolecules and enzymatic activities (Essiamah & Eschrich, 1986; Borchert et al., 2002). Zanon & Finger (2010), while studying the relationship between temperature and air relative humidity and the growth of Araucaria angustifolia, occurring in the same phytogeographic domain of C. fissilis in this study and having similar light needs in the early stages of life, also found larger increments associated with higher average monthly temperatures and periods with greater precipitation.

In regions where drought and rain periods occur regularly, the influence of precipitation can be seen more clearly. In the Amazon basin, Baker et al. (2017) explored possible variation factors in the periodicity tree-rings from Cedrela species and affirm that growth rates are probably influenced by the seasonality of the precipitations. In the same way, a research developed in a tropical dry forest in eastern Bolivia, showed that C. fissilis was sensitive to precipitation in almost all analyzed time scales, with leaf flushing and radial-increment rate related to a synchronous response to water availability (Mendivelso et al., 2016). In spite of the precipitation influence in tree growth is evident in areas with droughts, the effects may be different in gradients. Venegas-González et al. (2018), indicated that the precipitation affects differentially the growth of Cedrela depending on the sector in the gradient in which they are. In their study, tree growth performance is highly dependent to the dry season rainfall amounts in the most humid sector of the gradient, while sites settled in areas of lower summer temperatures, rainfall is the main determining factor.

The correlation of the PCA1 with the chronology was not significant, considering the entire period evaluated. However, the GLM demonstrated this trend, in which temperatures were the best growth predictors of C. fissilis. When separating in time periods, the correlation was significant until the year 1979, even with a smaller sample, in which the effects of the outliers may be larger (Goodwin & Leech, 2006). In the other periods, there was no significant correlation with this latent variable. The correlation between growth and temperatures in one period, followed by a lack of correspondence in the following years, may indicate some alteration in how tree growth was influenced by the climatic variables studied.

Observing the PCA1 scores, we can see an increase in temperatures from 1982, since this latent variable is composed by the mean, maximum and minimum temperatures. Being in a region with characteristics of seasonal temperatures and no drought, we could expect an increase in the growth of C. fissilis in response to local warming. Although the increase in temperatures did not lead a proportional increase in growth, it may have influenced the climatic response. Whereas no correlation was found from 1980, but there was correlation in the previous period, an indirect impact of temperature increase may have occurred, acting independently or in combination with another factors (Lindner et al., 2010; Bellard et al., 2012). This result could be related to the particularly sensitive of forests to climate change, because the long life-span of trees does not allow for rapid adaptation to environmental changes (Lindner et al., 2010).

Concerning to niche evolution of Cedrela, Koecke et al. (2013) demonstrated that temperature had an increasing impact on ecological diversification of the genus from the Miocene onwards, with recent speciation events may be related to divergence of climatic tolerances. The increase in temperature seems to be a major driver in promoting diversity on Tropical Rain Forests, apparently caused by in situ speciation (Jaramillo & Cárdenas, 2013). But what about longer-lived pioneer species, how do they react? Our results indicated that, from C. fissilis, changes in growth and in correlation with the temperature in response to local warming have occurred. The most negative impacts from a warmer climate are expected for slow-growing trees, and, the least negative for pioneer species, who gain a growth competitive advantage and moving into open habitats (Körner & Basler, 2010; Corlett & Westcott, 2013). Cheesman & Winter (2013) shown that increasing temperatures, especially at night, under well-watered conditions, promote an increase in biomass accumulation on pioneer tree species. On the other hand, C. fissilis is more suitably framed as a long-lived pioneer, because it has long life and dense wood, despite several characteristics that classify it as a pioneer (Corlett, 1995). As a whole, being a long-lived pioneer, C. fissilis did not grow proportionally to the increase in temperature but experienced a change in its growth pattern, suggesting sensitivity to climate changes.

Climatic factors such as temperature and precipitation are directly related to the physiological processes fundamental to plant growth (Reyer et al., 2012; Venegas-González et al., 2016). However, the relationship between these climatic elements and growth is not a direct one, since it is not possible to be sure if the resources will be available (Ribeiro et al., 2011). It is important to understand the factors controlling ring formation at smaller scales, since the variation of these factors between regions may cause trees in different regions to form rings at different times (Baker et al., 2017).

Understand the relative influence of specific climate parameters on tree growth and predict species’ responses to climate changes is a major challenge in ecology (Wolkovich et al., 2012; Williams et al., 2013). Even if studying these relationships is challenging, with obtaining additional tree-ring chronologies from C. fissilis, it should be possible to make climate reconstructions covering the past 150 to 200 years (Pereira et al., 2018). This advance could improve our understanding of tree-rings as proxies for climate reconstructions, identifying past and future trends in climate changes (Fan et al., 2009; Liang et al., 2009).


Temperature is the main climatic variable that affects the radial growth of Cedrela fissilis. The seasonal behavior of this variable seems to be the main reason for its role as the best radial growth predictor of the species in the region. The fact that there is no well-defined dry season contributed to the non-occurrence of significant responses to air relative humidity and precipitation.

Our results indicated a change in the C. fissilis growth pattern after the year of 1980, a corresponding period to an increase in temperatures, suggest a sensitivity of the species to climate changes.


The authors would like to thank Copel (the electric company of the state of Paraná) for providing the study area and for the research support. We thank the Brazilian Research Council (CNPq) for the productivity grant to the third author (process 303849/2009-6) and for the PhD grant to the first author (process 140849/2015-7).

FINANCIAL SUPPORT Conselho Nacional de Desenvolvimento Científico e Tecnológico, (Grant/Award Number: '140849/2015-7','303849/2009-6').


Akaike H. A new look at the statistical model identification. IEEE Transactions on Automatic Control 1974; 19(6): 716-723. [ Links ]

Allen CD, Macalady AK, Chenchouni H, Bachelet D, McDowell N, Vennetier M et al. A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. Forest Ecology and Management 2010; 259(4): 660-684. [ Links ]

Alvares CA, Stape JL, Sentelhas PC, Moraes G, Leonardo J, Sparovek G. Köppen’s climate classification map for Brazil. Meteorologische Zeitschrift 2013; 22(6): 711-728. [ Links ]

Andreacci F, Botosso PC, Galvão F. Fenologia vegetativa e crescimento de Cedrela fissilis na Floresta Atlântica, Paraná, Brasil. Floresta e Ambiente 2017; 24(0): 1-11. [ Links ]

Andreacci F, Botosso PC, Galvão F. Sinais climáticos em anéis de crescimento de Cedrela fissilis em diferentes tipologias de Florestas Ombrófilas do Sul do Brasil. Floresta 2014; 44(2): 323-332. [ Links ]

Baker JCA, Santos GM, Gloor M, Brienen RJW. Does Cedrela always form annual rings? Testing ring periodicity across South America using radiocarbon dating. Trees (Berlin) 2017; 31(6): 1999-2009. [ Links ]

Begum S, Nakaba S, Yamagishi Y, Oribe Y, Funada R. Regulation of cambial activity in relation to environmental conditions: Understanding the role of temperature in wood formation of trees. Physiologia Plantarum 2013; 147(1): 46-54. PMid:22680337. [ Links ]

Bellard C, Bertelsmeier C, Leadley P, Thuiller W, Courchamp F. Impacts of climate change on the future of biodiversity. Ecology Letters 2012; 15(4): 365-377. PMid:22257223. [ Links ]

Billings WD. The environmental complex in relation to plant growth and distribution. The Quarterly Review of Biology 1952; 27(3): 251-265. PMid:13004274. [ Links ]

Borchert R. Water status and development of tropical trees during seasonal drought. Trees 1994; 8(3): 115-125. [ Links ]

Borchert R, Rivera G, Hagnauer W. Modification of Vegetative Phenology in a Tropical Semi-deciduous Forest by Abnormal Drought and Rain. Biotropica 2002; 34(1): 27-39. [ Links ]

Brasil. Agencia Nacional de Águas. Sistema Nacional de Informações sobre Recursos Hídricos. Portal HidroWeb. Brasília; 2017. [cited 2017 June 19]. Available from: ]

Cedro A, Bosiacka B, Myśliwy M. Dendrochronological analysis of three pine species used as pioneer species to stabilize the coastal dunes of the southern Baltic coast. Baltic Forestry 2013; 19(2): 226-235. [ Links ]

Chahud A, Petri S. Anfíbio e Palaeonisciformes da Porção Basal do Membro Taquaral, Formação Irati (Permiano), Estado de São Paulo, Brasil. Geologia USP. Série Científica 2010; 10(1): 29-37. [ Links ]

Cheesman AW, Winter K. Elevated night‐time temperatures increase growth in seedlings of two tropical pioneer tree species. The New Phytologist 2013; 197(4): 1185-1192. PMid:23278464. [ Links ]

Cook ER, Holmes RL. Users manual for program ARSTAN. Tucson: The University of Arizona Press; 1986. [ Links ]

Corlett RT, Westcott DA. Will plant movements keep up with climate change? Trends in Ecology & Evolution 2013; 28(8): 482-488. PMid:23721732. [ Links ]

Corlett RT. Tropical secondary forests. Progress in Physical Geography 1995; 19(2): 159-172. [ Links ]

Costa MS, Vasconcellos TJ, Barros CF, Callado CH. Does growth rhythm of a widespread species change in distinct growth sites? IAWA Journal 2013; 34(4): 498-509. [ Links ]

Cui LG, Shan JX, Shi M, Gao JP, Lin HX. The miR156‐SPL 9‐DFR pathway coordinates the relationship between development and abiotic stress tolerance in plants. The Plant Journal 2014; 80(6): 1108-1117. PMid:25345491. [ Links ]

Cusatis AC, Trazzi PA, Dobner Júnior M, Higa AR. Dendroecologia de Cedrela fissilis na Floresta Ombrófila Mista. Pesquisa Florestal Brasileira 2013; 33(75): 287-297. [ Links ]

Elling W, Dittmar C, Pfaffelmoser K, Rötzer T. Dendroecological assessment of the complex causes of decline and recovery of the growth of silver fir (Abies alba Mill.) in Southern Germany. Forest Ecology and Management 2009; 257(4): 1175-1187. [ Links ]

Essiamah S, Eschrich W. Water uptake in deciduous trees during winter and the role of conducting tissues in spring reactivation. IAWA Journal 1986; 7(1): 31-38. [ Links ]

Fan ZX, Bräuning A, Cao KF, Zhu SD. Growth–climate responses of high-elevation conifers in the central Hengduan Mountains, southwestern China. Forest Ecology and Management 2009; 258(3): 306-313. [ Links ]

Flores TB. Meliaceae in Flora do Brasil 2020 em construção [online]. Rio de Janeiro: Jardim Botânico do Rio de Janeiro; 2018. [cited 2018 out. 01]. Available from: <> [ Links ]

Folguera G, Bastías DA, Bozinovic F. Impact of experimental thermal amplitude on ectotherm performance: Adaptation to climate change variability? Comparative Biochemistry and Physiology. Part A, Molecular & Integrative Physiology 2009; 154(3): 389-393. PMid:19622394. [ Links ]

Frank HT, Gomes MEB, Formoso MLL. Review of the areal extent and the volume of the Serra Geral Formation, Paraná Basin, South America. Pesquisas em Geociências 2009; 36(1): 49-57. [ Links ]

Fritts HC. Tree Rings and Climate. San Diego: Academic Press; 1976. [ Links ]

Gebrekirstos A, Brauning A, Sass-Klassen U, Mbow C. Opportunities and applications of dendrochronology in Africa. Current Opinion in Environmental Sustainability 2014; 6: 48-53. [ Links ]

Gholami V, Chau KW, Fadaee F, Torkaman J, Ghaffari A. Modeling of groundwater level fluctuations using dendrochronology in alluvial aquifers. Journal of Hydrology (Amsterdam) 2015; 529(3): 1060-1069. [ Links ]

Goodwin LD, Leech NL. Understanding correlation: Factors that affect the size of r. Journal of Experiential Education 2006; 74(3): 249-266. [ Links ]

Gris D, Temponi LG, Marcon TR. Native species indicated for degraded area recovery in Western Paraná, Brazil. Revista Árvore 2012; 36(1): 113-125. [ Links ]

Holmes R. Computer-assisted quality control in tree-ring dating and measurement. Tree-Ring Bulletin 1983; 43: 69-78. [ Links ]

Intergovernmental Panel on Climate Change – IPCC. Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Geneva: IPCC; 2014. 151 p. [ Links ]

Jaramillo C, Cárdenas A. Global warming and neotropical rainforests: A historical perspective. Annual Review of Earth and Planetary Sciences 2013; 41(1): 741-766. [ Links ]

Jeong SJ, Ho CH, Gim HJ, Brown ME. Phenology shifts at start vs. end of growing season in temperate vegetation over the Northern Hemisphere for the period 1982-2008. Global Change Biology 2011; 17(7): 2385-2399. [ Links ]

Jombart T, Pontier D, Dufour AB. Genetic markers in the playground of multivariate analysis. Heredity 2009; 102(4): 330-341. PMid:19156164. [ Links ]

Kaiser HF. The varimax criterion for analytic rotation in factor analysis. Psychometrika 1958; 23(3): 187-200. [ Links ]

Kalnay EC, Kanamitsu M, Kistler R, Collins W, Deaven D, Gandin L et al. The NCEP/NCAR Reanalysis 40-year Project. Bulletin of the American Meteorological Society 1996; 77(3): 437-471.<0437:TNYRP>2.0.CO;2. [ Links ]

Kanieski MR, Longhi-Santos T, Milani JEF, Miranda BP, Galvão F, Botosso PC et al. Crescimento diamétrico de Blepharocalyx salicifolius em Remanescente de Floresta Ombrófila Mista Aluvial, Paraná. Floresta e Ambiente 2013; 20(2): 197-206. [ Links ]

Kanieski MR, Santos TL, Neto JG, Souza T, Galvão F, Roderjan CV. Influência da precipitação e da temperatura no incremento diamétrico de espécies florestais aluviais em Araucária-PR. Floresta e Ambiente 2012; 19(1): 17-25. [ Links ]

KNMI. KNMI Climate Explorer [online]. De Bilt: KNMI; 2017 [cited 2017 Sep 3]. Available from: ]

Koecke AV, Muellner‐Riehl AN, Pennington TD, Schorr G, Schnitzler J. Niche evolution through time and across continents: The story of Neotropical Cedrela (Meliaceae). American Journal of Botany 2013; 100(9): 1800-1810. PMid:24018859. [ Links ]

Konig Brun FG, Brun EJ, Longhi SJ, Gorenstein MR, Maria TRBC, Rêgo GMS et al. Vegetação arbórea em remanescentes florestais urbanos: Bosque do Lago da Paz, Dois Vizinhos, PR. Pesquisa Florestal Brasileira 2017; 37(92): 503-512. [ Links ]

Körner C, Basler D. Phenology under global warming. Science 2010; 327(5972): 1461-1462. PMid:20299580. [ Links ]

Li ZS, Zhang QB, Ma K. Tree-ring reconstruction of summer temperature for AD 1475–2003 in the central Hengduan Mountains, Northwestern Yunnan, China. Climatic Change 2012; 110(1-2): 455-467. [ Links ]

Liang EY, Shao XM, Xu Y. Tree-ring evidence of recent abnormal warming on the southeast Tibetan Plateau. Theoretical and Applied Climatology 2009; 98(1-2): 9-18. [ Links ]

Lima AG. Erosão fluvial sobre rochas vulcânicas: algumas inferências a partir de segmentos côncavos de perfis longitudinais. Revista Brasileira de Geociencias 2012; 42(1): 34-41. [ Links ]

Lima JAD, Santana DGD, Nappo ME. Comportamento inicial de espécies na revegetação da mata de galeria na Fazenda Mandaguari, em Indianópolis, MG. Revista Árvore 2009; 33(4): 685-694. [ Links ]

Lindner M, Maroschek M, Netherer S, Kremer A, Barbati A, Garcia-Gonzalo J et al. Climate change impacts, adaptive capacity, and vulnerability of European forest ecosystems. Forest Ecology and Management 2010; 259(4): 698-709. [ Links ]

López L, Villalba R. Reliable estimates of radial growth for eight tropical species based on wood anatomical patterns. Journal of Tropical Forest Science 2016; 28(2): 139-152. [ Links ]

Machado SA, Figura MA, Silva LCR, Nascimento RGM, Quirino SMS, Téo SJ. Dinâmica de crescimento de plantios jovens de Araucaria angustifolia e Pinus taeda. Pesquisa Florestal Brasileira 2010; 30(62): 165-170. [ Links ]

Mendivelso HA, Camarero JJ, Gutiérrez E, Castaño-Naranjo A. Climatic influences on leaf phenology, xylogenesis and radial stem changes at hourly to monthly scales in two tropical dry forests. Agricultural and Forest Meteorology 2016; 216: 20-36. [ Links ]

Miranda TLG. Plano ambiental de conservação e uso do entorno de reservatório artificial UHE Governador Bento Munhoz da Rocha Neto. Curitiba: LACTEC; 2009. [ Links ]

Moss RH, Edmonds JA, Hibbard KA, Manning MR, Rose SK, van Vuuren DP et al. The next generation of scenarios for climate change research and assessment. Nature 2010; 463(7282): 747-756. PMid:20148028. [ Links ]

Muellner AN, Pennington TD, Koecke AV, Renner SS. Biogeography of Cedrela (Meliaceae, Sapindales) in Central and South America. American Journal of Botany 2010; 97(3): 511-518. PMid:21622412. [ Links ]

Paredes-Villanueva K, López L, Navarro Cerrillo RM. Regional chronologies of Cedrela fissilis and Cedrela angustifolia in three forest types and their relation to climate. Trees (Berlin) 2016; 30(5): 1581-1593. [ Links ]

Pausas JG, Austin MP. Patterns of plant species richness in relation to different environments: An appraisal. Journal of Vegetation Science 2001; 12(2): 153-166. [ Links ]

Pereira GDA, Barbosa ACMC, Torbenson MCA, Stahle DW, Granato-Souza D, Santos RMD et al. The climate response of Cedrela fissilis annual ring width in the Rio São Francisco basin, Brazil. Tree-Ring Research 2018; 74(2): 162-171. [ Links ]

Piraino S, Abraham EM, Diblasi A, Roig Juñent FA. Geomorphological-related heterogeneity as reflected in tree growth and its relationships with climate of Monte Desert Prosopis flexuosa DC woodlands. Trees 2015; 29(3): 903-916. [ Links ]

Quadro MFL, Dias MAFS, Herdies DL, Gonçalves LGG. Análise climatológica da precipitação e do transporte de umidade na região da ZCAS através da nova geração de reanálises. Revista Brasileira de Meteorologia 2012; 27(2): 152-162. [ Links ]

R Development Core Team. R: A languageand environment for statistical computing. Vienna: R Core Team; 2017. [ Links ]

Reyer CPO, Leuzinger S, Rammig A, Wolf A, Bartholomeus RP, Bonfante A et al. A plant’s perspective of extremes: terrestrial plant responses to changing climatic variability. Global Change Biology 2012; 19(1): 75-89. PMid:23504722. [ Links ]

Ribeiro CAD, Pezzopane JRM, Pezzopane JEM, Loos RA, Xavier AC, Cecílio RA et al. Delimitação de microrregiões agroclimáticas e suas relações com o potencial produtivo da cultura do eucalipto. Floresta 2011; 41(4): 779-786. [ Links ]

Rosemback R, Ferreira NJ, Shimabukuro YE, Conforte JC. Análise da dinâmica da cobertura vegetal na região Sul do Brasil a partir de dados MODIS/Terra. Revista Brasileira de Cartografia 2010; 62(2): 401-416. [ Links ]

Santarosa E, Oliveira JM, Roig FA, Pillar VD. Crescimento Sazonal em Araucaria angustifolia: Evidências Anatômicas. Revista Brasileira de Biociências 2007; 5(1): 618-620. [ Links ]

Shimamoto CY, Botosso PC, Amano E, Marques MCM. Stem growth rhythms in trees of a tropical rainforest in Southern Brazil. Trees (Berlin) 2016; 30(1): 99-111. [ Links ]

Tilman D. The Resource-Ratio hypothesis of plant succession. American Naturalist 1985; 125(6): 827-852. [ Links ]

Venegas-González A, Chagas MP, Anholetto CR Jr, Alvares CA, Roig FA, Tomazello M Fo. Sensitivity of tree ring growth to local and large-scale climate variability in a region of Southeastern Brazil. Theoretical and Applied Climatology 2016; 123(1-2): 233-245. [ Links ]

Venegas-González A, Roig FA, Lisi CS, Albiero-Junior A, Alvares CA et al. Drought and climate change incidence on hotspot Cedrela forests from the Mata Atlântica biome in southeastern Brazil. Global Ecology and Conservation 2018; 15(e00408): e00408. [ Links ]

Wahidin S, Idris A, Shaleh SRM. The influence of light intensity and photoperiod on the growth and lipid content of microalgae Nannochloropsis sp. Bioresource Technology 2013; 129: 7-11. PMid:23232218. [ Links ]

Walter A, Silk WK, Schurr U. Environmental effects on spatial and temporal patterns of leaf and root growth. Annual Review of Plant Biology 2009; 60(1): 279-304. PMid:19575584. [ Links ]

Way DA, Oren R. Differential responses to changes in growth temperature between trees from different functional groups and biomes: a review and synthesis of data. Tree Physiology 2010; 30(6): 669-688. PMid:20368338. [ Links ]

Wigley TML, Briffa KR, Jones PD. On the average value of correlated time series, with applications in dendroclimatology and hydrometeorology. Journal of Climate and Applied Meteorology 1984; 23(2): 201-213.<0201:OTAVOC>2.0.CO;2. [ Links ]

Williams AP, Allen CD, Macalady AK, Griffin D, Woodhouse CA et al. Temperature as a potent driver of regional forest drought stress and tree mortality. Nature Climate Change 2013; 3(3): 292-297. [ Links ]

Wold S, Esbensen K, Geladi P. Principal component analysis. Chemometrics and Intelligent Laboratory Systems 1987; 2(1): 37-52. [ Links ]

Wolkovich EM, Cook BI, Allen JM, Crimmins TM, Betancourt JL, Travers SE et al. Warming experiments underpredict plant phenological responses to climate change. Nature 2012; 485(7399): 494-497. PMid:22622576. [ Links ]

Zanon MLB, Finger CAG. Relação de variáveis meteorológicas com o crescimento das árvores de Araucaria angustifolia (Bertol.) Kuntze em povoamentos implantados. Ciência Florestal 2010; 20(3): 467-476. [ Links ]

Zweifel R, Zimmermann L, Zeugin F, Newbery DM. Intra-annual radial growth and water relations of trees: implication towards a growth mechanism. Journal of Experimental Botany 2006; 57(6): 1445-1459. PMid:16556628. [ Links ]

Received: August 17, 2018; Accepted: October 30, 2018

Amanda Koche MarconR. São Januário, 223, apto 31, CEP 80210-300, Curitiba, PR, Brasil e-mail:

Creative Commons License This is an Open Access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.