ISSN e:2007-4018 / ISSN print: 2007-3828



Vol. XXII, issue 3 September - December 2016

ISSN: ppub: 2007-3828 epub: 2007-4018

Scientific article

Basal area increment series of dominant trees of Pseudotsuga menziesii (Mirb.) Franco show periodicity according to global climate patterns

Castruita-Esparza, Luis U. 1 ; Correa-Díaz, Arian 1 * ; Gómez-Guerrero, Armando 1 ; Villanueva-Díaz, José 2 ; Ramírez-Guzmán, Martha E. 1 ; Velázquez-Martínez, Alejandro 1 ; Ángeles-Pérez, Gregorio 1

  • 1Colegio de Postgraduados, Campus Montecillo. Carretera México-Texcoco km 36.5. C. P. 56230. Montecillo, Texcoco, Estado de México, México.
  • 2Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias (INIFAP), Centro Nacional de Investigación Disciplinaria en Relación Agua, Suelo, Planta, Atmósfera (CENID RASPA). km 6.5 Margen Derecha Canal Sacramento. C. P. 35140. Gómez Palacio, Durango. México.

Corresponding author.

Received: October 21, 2015; Accepted: July 20, 2016

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


Tree species like Pseudotsuga menziesii (Mirb.) Franco (Douglas-fir) are sensitive to climate cycles and show well defined growth rings. The careful selection of dominant trees with circular Ttrunk section allow the analysis of tree growth trajectories. In this study, we used direct measurements of basal area increment (BAI) to explain biological periodicity and forecast basal area growth of Douglas-fir growing in Western Mexico. To remove the age effect on tree growth we also ran the analysis in terms of cambial age. Results showed significant (P < 0.05) correlation between BAI and precipitation from January to July. We found periodicities in tree growth of 7, 21, 27 and 60 years. However, the 60-year period, was determinant to build an ARIMA model (0,1,1), to forecast BAI for the next decades. Tree growth projections suggest reduced BAI in mature dominant trees for the next decades. Decreased tree-growth is an unexpected result, as BAI in dominant trees remains constant up to the biological age. Our finding is concurrent with a general decrement in tree growth in other forests of the world due to water stress, which suggests that the future climatic variability may worsen health conditions of Douglas-fir forests in North Mexico.

Keywords:Dendrochronology, tree-ring, forest productivity, temperate forests; time series analysis.


Due to the increased variability of the climate in the last three decades (Hansen, Sato, & Ruedy, 2012), more understanding on the relationship between tree-ring width, basal area growth and climatic variables is needed to propose better forest practices. Establishing relationships between tree-ring measurements and basal area increments (BAI) with climatic variables allow the identification of patterns in tree growth. Forest growth depends on tree species, age, stand density, soil quality and climate (Skovsgaard & Vanclay, 2008), and the combined effect of these factors is registered in the ring width of forest species (Esper, Benz, & Pederson, 2012; Matisons, Elferts, & Brūmelis, 2013). Tree-rings of climate sensitive species preserve information about cyclic patterns of tree growth and climate relationships, providing useful information for modeling and forecasting purposes (Matisons et al., 2013; Návar & Lizárraga-Mendiola, 2013).

One of the most important species for dendrochronology purposes has been Douglas-fir (Pseudotsuga menziesii ENT#091;Mirb.ENT#093; Franco), due to its sensitivity in capturing signals of cyclic events of climate (Cerano, Villanueva, Valdéz, Méndez, & Constante, 2011; Griesbauer & Scott, 2010; González-Elizondo et al., 2005; Gugger, González-Rodríguez, Rodríguez-Correa, Sugita, & Cavender-Bares, 2011) but the reliability of that periodicity has not been tested in forecasting models. Time series analysis of tree species like Douglas-fir would also help to confirm extreme past events in atmospheric circulation patterns, water resources shortage, and its socioeconomic impact on the population (Cerano et al., 2011; Návar & Lizárraga-Mendiola, 2013). More recently, the use of time series analysis combined with wood isotopes and dendrochronology have showed to be useful to describe the relationship of tree physiology and climate variability according to increased atmospheric CO2 (Monserud & Marshall, 2001; Roden, Johnstone, & Dawson, 2011; Swetnam & Lynch, 1993).

Basic dendrochronology procedures use standardized the ring width index (RWI) to analyze trends in tree growth; nevertheless, when dominant trees with trunks of circular section are selected, direct measurements of BAI are highly correlated to RWI, allowing the study of tree growth trends overtime. Compared to ring with BAI is more convenient to separate tree growth rates because is a two-dimensional variable (Biondi & Qeadan, 2008). When tree-ring series are transformed into basal area increments (BAI) the practical use is that BAI has units in cm2•year-1 which allow direct interpretation of the tree growth trends; compared to BAI, RWI values are relative units and dimensionless (Rubino & McCarthy, 2000).

Time series analysis tools are frequently used in econometric studies, however, due to recent climate variation and its effects on tree growth, the use of time series analysis should be revaluated as a powerful tool to understand climate change and tree growth relationships.

This paper aimed to analyze trends for BAI of mixed stands of Douglas-fir using time series models to forecast tree growth for the next decades. Our hypothesis is that: BAI series of P. menziesii in the Western Sierra of Mexico reflect regional climate events and are reliable to build a time series model to explain and forecast tree basal area growth.

Material and methods

Study area

The study site is located in the Western Sierra of Mexico in the Chihuahua state (26° 27’ 59.14’’ N and 106° 21’ 14.95’’ W). The altitude ranges from 2,650 to 2,730 m and the climate is temperate sub-humid, with mean annual temperature and precipitation of 12 °C and 780 mm, (Figure 1). Associated tree species of the study site are Picea chihuahuana Martínez, Pinus arizonica Engelm., Pinus durangensis Martínez and Pinus ayacahuite Ehrenb. (García-Arévalo, 2008; González, González, & Márquez, 2007). There is not logging in the study area; it has been under conservation as P. menziesii is an endangered species according to Mexican regulations (NOM-059-SEMARNAT-2010). Total height of the individual of Douglas-fir range in the upper stratum range from 23 to 45 m, with DBH (diameter at breast height) from 60-104 cm. The mixed stand of the study area has an estimated stock of 8.4 m3•ha-1 of bole of the study species.

Wood core sampling

In summer of 2012, 20 healthy and dominant trees carefully selected by their circular section of the bole were sampled with a 12 mm borer (HAGLÖF, BS009, Germany). For each tree, one sample in transversal direction to the slope was collected at 1.3 m above the soil level. In the study site, Douglas-fir is present in mixed and open stand where the canopy cover ranges from 60 to 80 %. Collecting dominant and mature trees we assured lower variation for BAI as the basal area growth of Douglas-fir plateau and remain constant up to the age of 300 years (Poage & Tappeiner, 2002). We avoid old trees close to the biological age of the species because they have low growth rates overtime and poor relationship with climate variation (Black, Colbert, & Pederson, 2008).

Sample measurements

Samples were air dried and polished with sand paper for better contrast of tree-rings. Ring width was measured with a Velmex measuring system (Bloomfield, NY, USA) with a 0.001 mm accuracy (Robinson & Evans, 1980). Increment cores samples were processed following standard dendrochronology techniques (Stokes & Smiley, 1968). Most of the wood core samples included the pith (75 %) and the age of trees without pith was estimated by regression according to the trend of the rings of trees with known age. The total number of tree-ring measurements was 2,480. The average age for the sampled trees was 140 years, ranging from 67 to 203 years. The intercorrelation of the tree-ring measurements series of the wood cores was statistically tested with the COFECHA software (Holmes, 1983). A negative exponential model was used to standardize the ring series for each tree and to develop a tree-ring width chronology using the ARSTAN software (Cook & Holmes, 1984). The “spline” for tree-ring indices was fitted with a ten-year resolution.

Figure 1:. Geographical location of the sampled site in the state of Chihuahua, Mexico.

Basal Area Increment (BAI) and Precipitation data

Basal area increments (BAI) were computed according to the following formula, where R is stem radius (cm) and n is the year of ring formation:

BAI increases in mature stages tend to be stable as long as the tree is not close to the extremes of juvenile stage or the biological senescence (Poage & Tappeiner, 2002; Weiner & Thomas, 2001). Although the age of sampled trees varied from 64 to 204 y, their age is far from the biological senescence for the study species which is at about 300 years (Villanueva-Diaz et al., 2010). In order to construct a representative BAI series for the study site, we averaged the BAI values at a yearly resolution for all samples. The maximum length of the BAI chronology was 203 years (1807 to 2010). The BAI and the standardized tree-ring with series (RWI) were tested for significant correlation (P < 0.05). Furthermore, the RWI and BAI chronology was related with the regional precipitation, for which the program DENDROCLIM 2002 was used (Biondi & Waikul, 2004). We used the information from the four closest meteorological stations as reported by the National Commission for Water in Mexico (Instituto Mexicano de Tecnología del Agua ENT#091;IMTAENT#093;, 2013). Additionally, to remove the effect for different tree ages, we reordered the data from calendar year to cambial age and fitted a five order polynomial equation to estimate the expected BAI according to the mean tree age in our sample (Briffa & Melvin, 2011).

Spectral analysis

A spectral analysis was performed to identify cycles of BAI. The time series was decomposed into the sum of sine and cosine waves with different amplitudes and lengths (Jenkins & Watts, 1968; Priestly, 1981):

The Fisher’s Kappa and Bartlett´s Kolmogorov-Smirnov statistics were computed for the BAI chronology assuring that the series and the spectrum were not “white noise” and were statistically significant. In addition, spectral density and periodogram graphs were developed and tested to determine the presence of significant low-frequency cycles. The spectral density estimate was produced by smoothing the periodogram.

ARMA and ARIMA Box-Jenkins models

ARMA (p, q) Box-Jenkins models rely on second-order stationary; the presence of p autoregressive terms (AR) and q moving average parameters (MA) help to explain the behavior of the response variable (Box & Jenkins, 1970).

Second-order stationary models are mean ENT#091;E(Zt) =µENT#093;, variance ENT#091;Var(Zt) = σ2ENT#093; and covariance ENT#091;Cov(Zt,Zs) = γkENT#093; not time dependent. One extension to the ARMA models is the class of autoregressive integrated moving average, or ARIMA process with dth difference. The d term is the number of times the series must be differentiated to be stationary. The stochastic trend can be removed applying a difference operator ∇d:

where Xt is the d times differenced series to be transformed into stationary series and Zt is the original series. The augmented Dickey-Fuller test of unit root proves that the time series is stationary, testing the hypothesis H0: ρ = 1 and the alternative H1: ρ < 1. If non-stationarity is associated with a variable variance we can use a transformation of power type InZt (Pankrats, 1983). It should be noticed that when it is necessary to remove trend and variable variance, the transformation of actual values into logarithms is required before the differentiation.

ARIMA models can describe cycles for annual data. It is convenient to use a cycle difference operator to achieve cycle stationarity:

where Wt is the D times differenced series (length of the cycle C). Cyclical effects ARIMA models can be proposed under the following general model ARIMA (p, d, q) x (P, D, Q)C :

In this work, we analyzed BAI chronology by autoregressive integrated moving average models (ARIMA Box-Jenkins Models) to find major cycles or periodicity in tree growth.

Estimation of ARIMA model components

The ARIMA modeling involved the comparison of estimated autocorrelation function graphs (ACF) and partial autocorrelation (ACFp) with theoretical ACF and ACFp. The pattern of decay of those functions and the lags with significant values are components to determine the order of parameters for the ARMA model (Box & Jenkins, 1970). We used the ACF to determine the moving average (MA) order (q) and the ACFp was used to determine the order (p) of the autoregressive model. Complementarily, the smallest canonical (SCAN) correlation method was used. This process is helpful to identify the orders of a stationary and non-stationary ARMA model. In general this method consist in identifying a rectangular pattern of a table which has the eigenvalues of a matrix product of the vector Ym,t = (Zt, Zt-1,... Zt-m) where Zt is time series, such that m = Pmin,...,Pmax and j = qmin,...,qmax; where m and j identify the possible order of the AR and MA process, respectively (Tsay & Tiao, 1984). Finally, the estimation of parameters was performed using conditional least squares (CLS) (Statistical Analysis System ENT#091;SASENT#093;, 2011).

Strategy for diagnostic, checking and forecasting

The results of spectral analysis were considered to test the different cyclic component (C) of the ARIMA process and generate an appropriate model for understanding the behavior of tree growth. The statistical reliability of the model was verified through a chi-squared test (Ljung-Box test) by assuring that random shocks and ut were statistically independent (not autocorrelated) and “white noise”. The forecast was calculated with the infinite memory forecast, also called conditional forecast. The term conditional is used because the forecast is computed by assuming that the unknown values of the response series before the start of the data are equal to mean of the series. We performed a forecast of BAI for the next 15-year period and up to the year 2090 at the 95 % confidence level.

Results and discussion

Tree ring series, BAI and precipitation

Results showed that our strategy for selecting dominant and circular trunk section trees was reliable for dendrochronological purposes. The strength of the signal in ARSTAN was high with a value of 0.90, indicating the existence a common response in all the trees sampled (Cook, 1983, 1987). The inter-correlation for the series and the average mean sensitivity computed with COFECHA were adequate (P < 0.01), with values of 0.67 and 0.27, respectively. These values indicate that the ring measurements show a common trend and the suitability for cross-dating of the series (Figure 2a). The correlation between RWI and BAI was statistically significant (r = 0.80; P < 0.01) indicating that in this case direct BAI measurements are reliable to describe tree growth changes and useful for time series analysis (Figure 2b). RWI and BAI were also significantly (P < 0.05) correlated to the precipitation for the period January to July with values of 0.73 and 0.70, respectively. These results also confirms the sensitivity of Douglas-fir to climate variability and the suitability to use their tree ring series to predict tree growth series, as it has been found for P. menziesii trees growing in Durango and Zacatecas, Mexico (González-Elizondo et al., 2005).

Tree growth trends

Mean BAI of P. menziesii was 54 cm2•year-1 (SD ±12). This result is similar to that reported for stands of Douglas-fir growing western Oregon where trees age ranges from 100 to 300 years and plateau at 50 cm2•year-1 (Poage & Tappeiner, 2002). Our findings show a trend of reduced BAI overtime, which is an unexpected results because the sample was composed by healthy and dominant trees with average age of 140 years (range, 67-204 years); this is 150 years less than the natural senescence of the species. A reduction in BAI has also been reported for other forests suggesting that water and nutrient stress are one of the stress factors (Gomez-Guerrero et al., 2013; Park et al., 2013; Peñuelas, Canadell, & Ogaya, 2011; Silva, Gómez-Guerrero, Doane, & Horwath, 2015). A strong link between growth diameter and moisture has also been found for P. menziesii at the Eastern Sierra of Mexico with larger deviations in diameter growth according to extreme wet and dry episodes (Arreola-Ortiz, 2010).

Figure 2:. (a) Douglas-fir standard ring width chronology and sample size (dashed line), covering the period 1850 to 2010. (b) Correlation between tree ring width index and basal area increments.

Analyzing the results at a decadal resolution the reductions in BAI with respect to the previous decade, were 24 and 28 % for the periods of 1890-1899 and 1950-1959, respectively. Favorable periods resulted in increased BAI of about 16 % above the mean value (Figure 3a). In a percent basis and absolute values, the impact of favorable periods is lower than that of dry periods because the benefit of wet periods on BAI is limited to the water storing capacity of the soils (Castruita, 2014). The decadal drought spell of the 1950’s and 1990’s has also been reported in river flow analysis in Mexico (Návar, 2012) and soil moisture content evaluations using hydrological models (Návar, 2015).

The reduction in BAI seen in this study is not related to the age of the trees because according to the cambial age in Figure 3b, the expected average in BAI for trees 140-year old is 50 cm2•year-1, very close to the mean value of 54 cm2•year-1. Thus, the reduced trend in BAI overtime is more likely due to the climate effect because RWI and BAI were statistically correlated and the precipitation from January to July. According to our results, water stress by major droughts might be the reason for important reductions of BAI in 1893, 1953 and 1996 to values of 23, 19 and 25 cm2•year-1, respectively. Prolonged droughts has been reported for these years; for example, a severe dry period around the year 1950 has been documented in North México (Díaz, Therrell, Stahle, & Cleaveland, 2002) affecting the south of United States and that for middle of the 90´s has also been detected in Taxodium mucronatum (Ten.) trees in Central Mexico (Stahle et al., 2009; Stahle et al., 2011).

New revisions for the standardization of tree-ring measurements suggest that the selection of dominant trees could lead to overestimations in tree growth (Nehrbass‐Ahles et al., 2014). We think that this is not the case of our sample of trees because ordering the data according to cambial age we did not see a disproportional difference in BAI between age classes, especially after the cambial age have reached 50 years (Figure 3b). The decreasing trends in tree growth sometimes is due to the inclusion of long-lived trees with slow growth rates (Black et al., 2008). However, the oldest age class of trees in this study (181-204 years) show expected BAI of 50 cm2•year-1 at cambial ages of 150 years, supporting that our sample was composed by mature (no old) trees with average growth rates for this age.

Figure 3: (a). Mean basal area increment (BAI) trend according Chapman-Richards equation Y=(1-e-0.0043*t)1.1 P < 0.05 and (b) BAI and cambial age by age classes.

Spectral Analysis and Periodicity

The Fisher’s Kappa (P < 0.05) and the Bartlett´s Kolmogorov-Smirnov statistic were statistically significant (P < 0.0001), indicating that the series and their spectra are not white noise and the analysis is unbiased (Fuller, 1976). Results indicated periodicities of 7, 21, 27 and 60 years (Figure 4). Because we are analyzing a biological process, these periodicities represent approximated cycles. However, what is very remarkable is that the last cycle is in agreement with droughts recurrence of 60 years reported for northern Mexico (Cerano et al., 2011). Although it is difficult to establish a direct relationship with global circulatory events, the fact that deserves special attention is that the cycles found here have showed to be compatible with ENSO (El Niño Southern Oscillation) and PDO (Pacific Decadal Oscillation) cycles (Hare & Mantua, 2000; Návar & Lizárraga-Mendiola 2013; Trenberth, 1997). One explanation for this results might be the sensitivity of the study species in relation to moisture variations influenced by regional events (Griesbauer & Scott, 2010; Návar, 2015). Our results match with cycles of 7, 21 and 60-y associated to other global circulatory events like ENSO, PDO and AMO (respectively), as suggested by other authors (Návar & Lizárraga-Mendiola, 2013; Návar, 2015).

Figure 4:. (a) Peridiogram and (b) spectral density for basal area increment series of Pseudotsuga menziesii showing the most significant cycles (P < 0.05).

IAB Stationary transformation

The Dickey-Fuller Unit Root Tests and the autocorrelation function (ACF) of the original series indicated nonstationarity (Table 1). The transformation of the BAI into natural logarithm and its first difference (Yt) allowed that the series became grouped around the mean with a homogeneous variance. The logarithmic transformation is convenient to decrease heterocedasticity and the dth difference is suggested when the mean of series changes overtime (Pankrats, 1983). The ACF of the original Douglas-fir BAI chronology decays to no significant levels in three years, suggesting that the growth of the previous three years is correlated to current tree growth and diameter. Transformation and first differencing resulted in different trend of ACF, with significant negative autocorrelation at lag 1 year (Figure 5).

Table 1:. Augmented Dickey-Fuller Unit Root Tests of the original BAI chronology.

Type Lags Rho Pr < Rho Tau Pr < Tau F Pr > F
0 73.40 0.0013 −4.80 <0.0001 26.25 0.001
Single Mean 1 −49.21 0.0013 −5.46 <0.0001 14.97 0.001
2 −53.23 0.0013 −5.49 <0.0002 11.58 0.001
0 −70.75 0.0005 −6.09 <0.0001 32.51 0.001
Trend 1 −73.06 0.0005 −7.24 <0.0001 18.70 0.001
2 −89.03 0.0005 −8.04 <0.0001 15.16 0.001

Figure 5:. Autocorrelation function (ACF) and partial autocorrelation function (ACFp) of the original (top) and transformed basal area increment (BAI) chronology (natural logarithm and first difference) (bottom) of Pseudotsuga menziesii.

Identification and estimation ARIMA models

The pattern of decay for ACF and ACFp of the transformed BAI suggested a MA(1) model. The features of the model were defined by the spike at lag 1, that cuts off to 0 afterwards (Figure 5). The fact that the spike was negative with negative exponential trend (ACFp) indicates that θ > 0 (ACF) (Pankrats, 1983). This result was consistent with the smallest canonical (SCAN) correlation method that suggests an ARIMA (0,1,1) model. After considering the ACF and ACFp the selected model was ARIMA (0,1,1). This model has a moving mean parameter (q) involving random shocks, which is different from the AR(1) or ARMA (1,1) models that are frequently proposed for ring-width series (Meko, 1981). The selected model was evaluated with all the periodicities found (7, 21, 27 and 60) and without them to test their significance. Considering t and P values, all models were significant (C = 7, 21, 27 and 60). The moving average term (θ) for all cases was positive with values from 0.46 to 0.76, while the cyclic moving average term (Θ) varied from 0.45 to 0.83 (Table 2).

Table 2:. Adjustment of models for Pseudotsuga menziesii.

Model Parameter* Standard error T value Pr > |t| Variance estimation Standard error of the estimation AIC** SBC**
( 1 ( 1 (P/(Q
ARIMA (0,1,1) 0.503 -- -- 0.068 7.38 <0.0001 0.043 0.207 -49.23 -46.14
ARIMA (0,1,1) (0,1,1)7 0.458 -- -- 0.071 6.38 <0.0001 0.050 0.225 -20.01 -13.92
-- 0.832 0.047 17.54 <0.0001 -- -- -- --
ARIMA (0,1,1) (0,1,1)21 0.647 -- -- 0.065 9.97 <0.0001 0.054 0.232 -8.81 -2.91
-- 0.645 0.075 8.53 <0.0001 -- -- -- --
ARIMA (0,1,1) (0,1,1)27 0.492 -- -- 0.076 6.49 <0.0001 0.063 0.251 11.75 17.56
-- 0.780 0.074 10.46 <0.0001 -- -- -- --
ARIMA (0,1,1) (0,1,1)60 0.764 -- -- 0.065 11.59 <0.0001 0.060 0.246 5.72 10.97
-- 0.448 0.166 2.70 <0.0001 -- -- -- --
* ( 1 is the noncyclic moving average term, ( 1 is the noncyclic autoregressive operator, (P is the cyclic autoregressive operator of length C and (Q is the cyclic moving average operator of length C. ** Akaike’s information criterion (AIC) and Schwartz’s Bayesian criterion (SBC).

For forecasting purposes the selected model was ARIMA (0,1,1) × (0,1,1)60. Two reasons support this selection: 1) First, its statistical strength is reflected in the significance of its terms, no correlated residuals, the Akaike’s information criterion (AIC) and Schwartz’s Bayesian criterion (SBC) (lowest values), and reliable fitting of the data; 2). Second, the periodicity of 60 years has been related to global circulatory patterns by other authors for the studied species (Cerano et al., 2011) and has a major impact in the tree growth at the study area.

Forecast of BAI

With the best statistical and climatic ARIMA model ENT#091;ARIMA (0,1,1) x (0,1,1)60ENT#093;, we performed a forecast for BAI in the next eight decades using the SAS 9.3 software (SAS, 2011), with the following equations:

Using coefficients from Table 2:

This model ponder BAI based on the function of random shocks, subtracting 76 % of the same component of the previous year, 45 % of the occurred 60 years ago and adding 34 % of the occurred 61 years ago. In the long term, the forecast indicates reductions of 72 % in BAI with an average of 15 cm2•year-1, and a recovery about the year 2025. However the projected growth in 2025 for this species is below the historical mean (54 cm2•year-1) suggesting a net reduction in growth for the next 75 years (Figure 6). Nevertheless, the amplitude for the confidence intervals limits the interpretation for several decades.

Because our results are assuming that the frequency of other stressing agents of the forests, such as pests and fires, remains constant, the health of the studied forest may be at risk with concomitant negative effects from other factors. The prevalence of healthy forests requires a base line study for tree growth as the one performed here.

Figure 6:. Extended forecast of basal area increment of Pseudotsuga menziesii from the ARIMA model. The black line is for real measurements and red line for estimations. The grey area defines the confidence interval at 95 %.

Tree growth trends found here do not show indications of combined and beneficial effects on increased atmospheric CO2 and temperature as found in other studies (Soulé & Knapp, 2015). Our tree sample showed four phases of maximum tree growth which were decreasing overtime (approximately 1980-2000). The major implication of these results is the expected impacts of climate variation in the short-term future. The forecast of the model proposed is not considering implicitly climate variation and the trend for reduced tree growth may be an indication of a generalized response of the forest in other parts of the world (Silva, Anand, & Leithead, 2010).

We think that our data for BAI, RWI and precipitation were highly correlated because we used direct BAI for the year of ring formation. Regardless of the ages of the trees, all rings formed under the same environmental conditions reflect the strong effect of precipitation. Another reason is that tree growth rates derived from a wide span of tree ages included. According to some authors the analysis by cambial age indicated that BAI of Douglas-fir tends to be less sensitive to climate as tree ages (Black et al., 2008), but our results indicated that a sample composed by mature and dominant trees from open stands is reliable to detect a common signal in trees. One of the best methods for selecting a sample of trees for dendrochronological purposes and tree growth estimation is the random selection that include different ages and diameters (Nehrbass‐Ahles et al., 2014). We think that in that respect our selection method for tall and healthy trees in the mixed and open stands of Doulas-fir leads to a sample heterogeneous in age and diameters (40-104 cm) that worked well for the purposes of the study.

The approach used in this study shows potential to be used with other forests species to investigate climate change effects at geographical scales. The trends in tree growth we found are similar to that reported by other authors which indicated that extreme drought events like the ones in 1893 and 1952 were geographically important North Mexico and Southwest USA (Cerano et al., 2011; Stahle et al., 2009).

Although our study was not tested in a regional scale, our approach could have also some potential to correlate series of BAI with plant vigor indexes like the NDVI as proposed by other authors (Berner, Beck, Bunn, & Goetz, 2013). BAI series, should be reconsidered as proxy variable to study climate change effects particularly for sensitive forests species. Indeed, the new standardization process in dendrochronology proposes to standardize in terms of the stable BAI observed in the middle age of forests species (Biondi & Qeadan, 2008).


We found strong correlation among tree ring indices, basal area growth and precipitation for the studied species. Time series analysis of BAI of Douglas-fir of the Western Sierra of Mexico allowed identifying periodicity in tree growth which match with global atmospheric events like the ENSO cycles. The forecast for BAI indicates decreased tree growth within the present decade, although some degree of recovery will be reached by year 2025. However the projections of tree growth in the long-term will not be as high as that seen in the 1960s (60 cm2•year-1). The reduction in BAI seen in this study is not related to the age of the trees because we also analyzed the trends according to the cambial age. The use of time series analysis of BAI from dominant and mature trees can be helpful to predict future changes in tree growth and to propose in advance more suitable management practices.


  • To the FORD foundation, the National Council of Science and Technology (CONACYT) of Mexico and the University of California that provided financial support, through the program UC-MEXUS (Proyecto CN-13-546), for carrying out this study.


Arreola-Ortiz, M. R., & Návar-Cháidez, J. J. (2010). Análisis de sequías y productividad con cronologías de Pseudotsuga menziesii Rob. & Fern., y su asociación con el niño en el nordeste de México. Investigaciones geográficas, 7-20. Retrieved from

Berner, L. T., Beck, P. S. A., Bunn, A. G., & Goetz, S. J. (2013). Plant response to climate change along the forest-tundra ecotone in northeastern Siberia. Global Change Biology, 19(11), 3449-3462. doi: 10.1111/gcb.12304

Biondi, F., & Qeadan, F. (2008). A theory-driven approach to tree-ring standardization: Defining the biological trend from expected basal area increment. Tree-Ring Research, 64(2), 81-96. doi: 10.3959/2008-6.1

Biondi, F., & Waikul, K. (2004). DENDROCLIM2002: A C++ program for statistical calibration of climate signals in tree-ring chronologies. Computers & Geosciences, 30(3), 303-311. doi: 10.1016/j.cageo.2003.11.004

Black, B. A., Colbert, J. J., & Pederson, N. (2008). Relationships between radial growth rates and lifespan within North American tree species. Ecoscience, 15(3), 349-357. doi: 10.2980/15-3-3149

Box, G. E. P., & Jenkins, G. M. (1970). Time series analisys, forecasting and control. San Francisco, USA: Holden-Day.

Briffa, K. R., & Melvin, T. M. (2011). A closer look at regional curve standardization of tree-ring records: Justification of the need, a warning of some pitfalls, and suggested improvements in its application. In K. M., Hughes, W. T., Swetnam & F. H.,Diaz (Eds.), Dendroclimatology: Progress and Prospects (pp. 113-145). Dordrecht, UK: Springer Netherlands. doi: 10.1007/978-1-4020-5725-0_5

Castruita, E. L. U. (2014). Variabilidad climática, eficiencia de uso de agua intrínseca y crecimiento del área basal en Bosques del norte de México. (Tesis doctoral, Colegio de Postgraduados, México. Retrieved from

Cerano, P. J., Villanueva, D. J., Valdéz, C. R., Méndez, G. J., & Constante, G. V. (2011). Sequías reconstruidas en los últimos 600 años para el noreste de México. Revista Mexicana de Ciencias Agrícolas, 2, 235-249. Retrieved from

Cook, E. R. (1983). A time series analysis approach to tree-ring standardization. Ph. D. dissertation. University of Arizona, USA.

Cook, E. R. (1987). The decomposition of tree-ring series for environmental studies. Tree-Ring Bulletin, 47, 37-59. Retrieved from

Cook, E. R., & Holmes, R. L. (1984). Program ARSTAN and users manual. Tucson, AZ, USA: Laboratory of Tree-Ring Research, University of Arizona.

Díaz, S. C., Therrell, M. D., Stahle, D. W., & Cleaveland, M. K. (2002). Chihuahua (Mexico) winter-spring precipitation reconstructed from tree-rings, 1647- 1992. Climate Research, 22(3), 237-244. doi: 10.3354/cr022237

Esper, J., Benz, M., & Pederson, N. (2012). Influence of wood harvest on tree-ring time-series of Picea abies in a temperate forest. Forest Ecology and Management, 284, 86-92. doi: 10.1016/j.foreco.2012.07.047

Fuller, W. A. (1976). Introduction to statical time series. New York, USA: John Wiley & Sons, Inc.

García-Arévalo, A. (2008). Vegetación y flora de un bosque relictual de Picea chihuahuana Martínez del norte de México. Polibotánica, 25, 45-68. Retrieved from

Gomez-Guerrero, A., Silva, L. C. R., Barrera-Reyes, M., Kishchuk, B., Velazquez-Martinez, A., Martinez- Trinidad, T., Horwath, W. R. (2013). Growth decline and divergent tree ring isotopic composition (13C and 18O) contradict predictions of CO2 stimulation in high altitudinal forests. Global Change Biology, 19(6), 1748-1758. doi: 10.1111/gcb.12170

González-Elizondo, M., Jurado, E., Návar, J., González- Elizondo, M. S., Villanueva, J., Aguirre, O., & Jiménez, J. (2005). Tree-rings and climate relationships for douglas-fir chronologies from the sierra madre occidental, mexico: A 1681-2001 rain reconstruction. Forest Ecology and Management, 213, 39- 53. doi: 10.1016/j.foreco.2005.03.012

González, E. M. S., González, M. E., & Márquez, M. A. L. (2007). Vegetación y ecorregiones de Durango. México: Plaza y Valdés.

Griesbauer, H., & Scott, G. D. (2010). Assessing the climatic sensitivity of Douglas-fir at its northern range margins in British Columbia, Canada. Trees, 24(2), 375- 389. doi: 10.1007/s00468-009-0407-z

Gugger, P. F., González-Rodríguez, A., Rodríguez-Correa, H., Sugita, S., & Cavender-Bares, J. (2011). Southward Pleistocene migration of Douglas-fir into Mexico: Phylogeography, ecological niche modeling, and conservation of ‘rear edge’ populations. New Phytologist, 189(4), 1185-1199. doi: 10.1111/j.1469-8137.2010.03559.x

Hansen, J., Sato, M., & Ruedy, R. (2012). Perception of climate change. Proceedings of the National Academy of Sciences, 109(37), E2415-E2423. doi: 10.1073/pnas.1205276109

Hare, S. R., & Mantua, N. J. (2000). Empirical evidence for North Pacific regime shifts in 1977 and 1989. Progress in oceanography, 47(2), 103-145. doi: 10.1016/S0079-6611(00)00033-1

Holmes, R. L. (1983). Computer-assisted quality control in tree-ring dating and measurement. Tree-Ring Bulletin, 43, 69-78.

Instituto Mexicano de Tecnología del Agua (IMTA). (2013). Extractor Rápido de Información Climatológica III (ERIC) version 3.2. Jiutepec, Morelos, México: Author.

Jenkins, G. M., & Watts, D. G. (1968). Spectral analysis and its applications. San Francisco, USA: Holden-Day .

Matisons, R., Elferts, D., & Brūmelis, G. (2013). Pointer years in tree-ring width and earlywood-vessel area time series of Quercus robur-Relation with climate factors near its northern distribution limit. Dendrochronologia, 31(2), 129-139. doi: 10.1016/j.dendro.2012.10.001

Meko, D. M. (1981). Application of Box-Jenkins methods of time series analysis to the reconstruction of drought from tree rings. USA: University of Arizona.

Monserud, R. A., & Marshall, J. D. (2001). Time-series analysis of δ13C from tree rings. Time trends and autocorrelation. Tree Physiology, 21, 1087-1102. doi: 10.1093/treephys/21.15.1087

Návar, J. (2012). Modeling annual discharge of six Mexico’s northern rivers. Revista Ambiente & Água, 7, 36-50. doi: 10.4136/ambi-agua.705

Návar, J. (2015). Hydro-climatic variability and perturbations in mexico’s north-western temperate forests. Ecohydrology, 8, 1065-1072. doi: 10.1002/eco.1564

Návar, J., & Lizárraga-Mendiola, L. (2013). Hydro-climatic variability and forest fires in Mexico’s northern temperate forests. Geofísica internacional, 52, 5-20. doi: 10.1016/S0016-7169(13)71458-2

Nehrbass-Ahles, C., Babst, F., Klesse, S., Nötzli, M., Bouriaud, O., Neukom, R., Frank, D. (2014). The influence of sampling design on tree-ring based quantification of forest growth. Global Change Biology, 20(9), 2867-2885. doi: 10.1111/gcb.12599

Pankrats, A. (1983). Forecasting with univariate Box-Jenkins models. Concepts and cases. USA: John Wiley and Sons.

Park, W. A., Allen, C. D., Macalady, A. K., Griffin, D., Woodhouse, C. A., Meko, D. M., McDowell, N. G. (2013). Temperature as a potent driver of regional forest drought stress and tree mortality. Nature Clim. Change, 3(3), 292-297. doi: 10.1038/nclimate1693

Peñuelas, J., Canadell, J. G., & Ogaya, R. (2011). Increased water-use efficiency during the 20th century did not translate into enhanced tree growth. Global Ecology and Biogeography, 20(4), 597-608. doi: 10.1111/j.1466-8238.2010.00608.x

Poage, N. J., & Tappeiner, J. C. (2002). Long-term patterns of diameter and basal area growth of old-growth Douglas-fir trees in western Oregon. Canadian Journal of Forest Research, 32(7), 1232-1243. doi: 10.1139/x02-045

Priestly, M. B. (1981). Spectral analysis and time series. New York, USA: Academic Press Inc.

Robinson, W. J., & Evans, R. (1980). A microcomputer-based tree-ring measuring system. Tree-Ring Bulletin, 40, 59- 64. Retrieved from

Roden, J. S., Johnstone, J. A., & Dawson, T. E. (2011). Regional and watershed-scale coherence in the stable-oxygen and carbon isotope ratio time series in tree rings of Coast Redwood (Sequoia sempervirens). Tree-Ring Research, 67(2), 71-86. doi: 10.3959/2010-4.1

Rubino, D. L., & McCarthy, B. C. (2000). Dendroclimatological analysis of white oak (Quercus alba L., Fagaceae) from an old-growth forest of Southeastern Ohio, USA. Journal of the Torrey Botanical Society, 127(3), 240-250. doi: 10.2307/3088761

Secretaría de Medio Ambiente y Recursos Naturales (2010). Norma Oficial Mexicana NOM-059-SEMARNAT-2010. Protección ambiental, especies nativas de flora y fauna silvestres de México, categorías de riesgo y especificaciones para su inclusión, exclusión o cambio, y lista de especies en riesgo. México: Diario Oficial de la Federación.

Silva, L. C. R., Anand, M., & Leithead, M. D. (2010). Recent widespread tree growth decline despite increasing atmospheric CO2. PLoS ONE, 5(7), e11543. doi: 10.1371/journal.pone.0011543

Silva, L. C. R., Gómez-Guerrero, A., Doane, T. A., & Horwath, W. R. (2015). Isotopic and nutritional evidence for species- and site-specific responses to N deposition and elevated CO2 in temperate forests. Journal of Geophysical Research: Biogeosciences, 120(6), 1110-1123. doi: 10.1002/2014JG002865

Skovsgaard, J. P., & Vanclay, J. K. (2008). Forest site productivity: A review of the evolution of dendrometric concepts for even-aged stands. Forestry, 81(1), 13-31. doi: 10.1093/forestry/cpm041

Soulé, P. T., & Knapp, P. A. (2015). Analyses of intrinsic water-use efficiency indicate performance differences of ponderosa pine and douglas-fir in response to CO2 enrichment. Journal of Biogeography, 42,144-155. doi: 10.1111/jbi.12408

Stahle, D. W., Cook, E. R., Diaz, J. V., Fye, F. K., Burnette, D. J., Griffin, D., Heim, R. R. (2009). Early 21st- Century drought in Mexico. Eos, Transactions American Geophysical Union, 90(11), 89-90. doi: 10.1029/2009EO110001

Stahle, D. W., Diaz, J. V., Burnette, D. J., Paredes, J., Heim, R., Fye, F. K., Stahle, D. K. (2011). Major Mesoamerican droughts of the past millennium. Geophysical Research Letters, 38(5), 1-4. doi: 10.1029/2010GL046472

Statistical Analysis System (SAS). (2011). SAS/ETS 9.3 User’s Guide. Cary, NC, USA: SAS Institute Inc.

Stokes, M. A., & Smiley, T. L. (1968). An introduction to tree-ring dating, Chicago, USA: University of Arizona Press.

Swetnam, T. W., & Lynch, A. M. (1993). Multicentury, regional-scale patterns of western spruce budworms outbreaks. Ecological Monographis, 63(4), 399-424. doi: 10.2307/2937153

Trenberth, K. E. (1997). The definition of el nino. Bulletin of the American Meteorological Society, 78(12), 2771- 2777. Retrieved from

Tsay, R. S., & Tiao, G. C. (1984). Consistent estimates of autoregressive parameters and extended sample autocorrelation function for stationary and nonstationary ARMA models. Journal of the American Statistical Association, 79(385), 84-96. doi: 10.2307/2288340

Villanueva-Diaz, J., Cerano-Paredes, J., Stahle, D. W., Constante-García, V., Vázquez-Salem, L., Estrada-Ávalos, J., & Benavides-Solorio, J. D. (2010). Árboles longevos de México. Revista Mexicana de Ciencias Forestales, 1(2), 1-23. Retrieved from

Weiner, J., & Thomas, S. C. (2001). The nature of tree growth and the “age-related decline in forest productivity”. Oikos, 94(2), 374-376. doi: 10.1034/j.1600-0706.2001.940219.x


Figure 1:. Geographical location of the sampled site in the state of Chihuahua, Mexico.
Figure 2:. (a) Douglas-fir standard ring width chronology and sample size (dashed line), covering the period 1850 to 2010. (b) Correlation between tree ring width index and basal area increments.
Figure 3: (a). Mean basal area increment (BAI) trend according Chapman-Richards equation Y=(1-e-0.0043*t)1.1 P < 0.05 and (b) BAI and cambial age by age classes.
Figure 4:. (a) Peridiogram and (b) spectral density for basal area increment series of Pseudotsuga menziesii showing the most significant cycles (P < 0.05).
Figure 5:. Autocorrelation function (ACF) and partial autocorrelation function (ACFp) of the original (top) and transformed basal area increment (BAI) chronology (natural logarithm and first difference) (bottom) of Pseudotsuga menziesii.
Figure 6:. Extended forecast of basal area increment of Pseudotsuga menziesii from the ARIMA model. The black line is for real measurements and red line for estimations. The grey area defines the confidence interval at 95 %.


Table 1:. Augmented Dickey-Fuller Unit Root Tests of the original BAI chronology.
Type Lags Rho Pr < Rho Tau Pr < Tau F Pr > F
0 73.40 0.0013 −4.80 <0.0001 26.25 0.001
Single Mean 1 −49.21 0.0013 −5.46 <0.0001 14.97 0.001
2 −53.23 0.0013 −5.49 <0.0002 11.58 0.001
0 −70.75 0.0005 −6.09 <0.0001 32.51 0.001
Trend 1 −73.06 0.0005 −7.24 <0.0001 18.70 0.001
2 −89.03 0.0005 −8.04 <0.0001 15.16 0.001
Table 2:. Adjustment of models for Pseudotsuga menziesii.
Model Parameter* Standard error T value Pr > |t| Variance estimation Standard error of the estimation AIC** SBC**
( 1 ( 1 (P/(Q
ARIMA (0,1,1) 0.503 -- -- 0.068 7.38 <0.0001 0.043 0.207 -49.23 -46.14
ARIMA (0,1,1) (0,1,1)7 0.458 -- -- 0.071 6.38 <0.0001 0.050 0.225 -20.01 -13.92
-- 0.832 0.047 17.54 <0.0001 -- -- -- --
ARIMA (0,1,1) (0,1,1)21 0.647 -- -- 0.065 9.97 <0.0001 0.054 0.232 -8.81 -2.91
-- 0.645 0.075 8.53 <0.0001 -- -- -- --
ARIMA (0,1,1) (0,1,1)27 0.492 -- -- 0.076 6.49 <0.0001 0.063 0.251 11.75 17.56
-- 0.780 0.074 10.46 <0.0001 -- -- -- --
ARIMA (0,1,1) (0,1,1)60 0.764 -- -- 0.065 11.59 <0.0001 0.060 0.246 5.72 10.97
-- 0.448 0.166 2.70 <0.0001 -- -- -- --
* ( 1 is the noncyclic moving average term, ( 1 is the noncyclic autoregressive operator, (P is the cyclic autoregressive operator of length C and (Q is the cyclic moving average operator of length C. ** Akaike’s information criterion (AIC) and Schwartz’s Bayesian criterion (SBC).