Next Article in Journal
Spruce Beetle Biology, Ecology and Management in the Rocky Mountains: An Addendum to Spruce Beetle in the Rockies
Previous Article in Journal
Potentials for Mutually Beneficial Collaboration Between FIA Specialists and IEG-40 Pathologists and Geneticists Working on Fusiform Rust
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Large-Scale Regeneration Patterns of Pinus nigra Subsp. salzmannii: Poor Evidence of Increasing Facilitation Across a Drought Gradient

1
Centro de Capacitación y Experimentación Forestal C/. Vadillo-Castril, s/n., Cazorla 23470, Spain
2
Departamento de Sistemas Físicos, Químicos y Naturales, Universidad Pablo de Olavide, Ctra. de Utrera km 1, Sevilla 41002, Spain
*
Author to whom correspondence should be addressed.
Submission received: 24 October 2013 / Revised: 19 December 2013 / Accepted: 24 December 2013 / Published: 31 December 2013

Abstract

:
Tree recruitment is a key process underlying stand dynamics and sustainability in managed forests. Woody plant cover is known to affect the regeneration success of Pinus nigra, suggesting the existence of facilitative plant-plant interactions. The regeneration patterns of this Mediterranean pine were analyzed across its distribution area, using data from 3226 plots of the Spanish National Forest Inventory. We aimed to test the hypothesis that seedlings establishment occurs under higher values of either canopy or shrub cover in the driest populations, as predicted by the stress-gradient hypothesis. Data were analyzed by means of Generalized Linear Models and multivariate methods. Results revealed that regeneration failure occurs on a regional scale, and that regeneration is facilitated by tree canopy cover of 55%–80%. A non-linear pattern of interaction along an aridity gradient was identified, with competition at the wettest site, high facilitation at the mid-dry sites, and low facilitation at the driest site. Evidence suggests that some shrub species may facilitate recruitment in the harsher areas. Collectively, our results reduce the possibilities of adapting forest management to drying climates by the application of alternative silvicultural prescriptions involving canopy cover.

1. Introduction

Water shortage has been widely reported as the main cause impeding the establishment of juvenile plants in Mediterranean ecosystems [1,2,3]. Seedlings are highly vulnerable to drought stress, and Mediterranean-type ecosystems are characterized by a sharp seasonality of rainfall and temperature, with hot-dry summers [4]. Consequently, seedling mortality is found to be exacerbated both during the summer season and in open areas, where high irradiance increases evaporative demands and further reduces soil-water availability [5,6,7,8]. Studies focused on the recruitment dynamics of Mediterranean tree species have found increasing seedling mortality during the summer, compared with other periods of the year, and in open areas, compared with sites located under the canopy of woody plants [9,10]. Thus, Mediterranean tree species generally benefit from being planted in the proximity of shrubs during reforestation work, as most shrub species may buffer against high temperatures and radiation [8]. Similarly, the risks of excessive evapotranspiration, overheating, and photo-inhibition are ameliorated under the cover of tree canopies [2,10,11]. These positive plant-plant interactions between tree seedlings and woody species are examples of facilitation [12].
According to the stress-gradient hypothesis [13], facilitation should increase across gradients of abiotic stress, i.e., the shade of neighboring plants could be expected to increasingly ameliorate drought stress as conditions become drier. This prediction has been validated for a broad range of environments and species, including species of the Pinus genus [14,15,16]. However, the stress-gradient hypothesis has been challenged by empirical evidence coming from arid and semi-arid areas. In these dry environments, drought stress might increase in the vicinity of other plants, despite the benefits of shading, as neighboring plants reduce the overall water content of the soil by direct water uptake and rainfall interception [17,18]. Furthermore, the net facilitative outcome of shrub-seedling interactions along a climate gradient will depend on the interacting shrub species, as the nurse effect provided by a shrub is modified by its morphological traits, such as size, branch disposition, and spininess [19], and its phylogeny [14,20].
Contrasting with the stress-gradient hypothesis, the trade-off hypothesis predicts that shade will increase the negative effect of drought, as water balance and carbon uptake cannot be maximized simultaneously [21]. Plants living in poor light environments would need to allocate more resources to aerial parts, and less to roots, thereby increasing their potential for light interception, while diminishing their ability to capture water and survive drought. The trade-off hypothesis has been supported by different ecophysiological studies involving Mediterranean species [22,23]. As a way to reconcile the apparent opposing views of the stress-gradient and trade-off hypotheses, it has recently been proposed that drought effects are ameliorated at intermediate irradiance, while they are more severe at higher or lower light levels due to a switch from facilitation to competition at both ends of an aridity gradient. When this situation applies, the expected plant response to drought along the irradiance gradient should be non-linear (a humped-back shape model) [24]. However, it is also possible that biotic interactions become unimportant relative to the effect of the environment when environmental conditions become very severe, generating a collapse of facilitation [25]. Whatever the hypothesis considered, an understanding of facilitative plant-plant interactions is important for forest managers, as tree harvest and/or shrub clearance modify the amount of light reaching the forest floor and, thus, forest regeneration dynamics. This understanding might be of greater importance in a climate change scenario, which implies rising temperatures and altered rainfall patterns in several Mediterranean-type forest ecosystems [26]. In this respect, there is increasing concern about the long-term persistence of some managed Mediterranean mountain pinewoods, which exhibit limited regeneration dynamics [2,27].
Pinus nigra Arn. subsp. salzmannii Dunal (Franco) (P. nigra hereafter) is a Mediterranean mountain pine, commonly harvested [28], and its forests are included in the European Union list of natural habitats requiring specific conservation measures. Several small-scale studies have shown that P. nigra regeneration frequently fails under natural conditions, particularly in open microhabitats, but it is unknown whether such lack of regeneration is a common pattern at the regional scale [6,9,27]. The species is widely distributed along the eastern mountain ranges of the Iberian Peninsula. Here, summer precipitation decreases and temperature increases from north to south, resulting in a latitudinal gradient of aridity [29]. We hypothesize that the establishment of P. nigra seedlings should be more closely related to either canopy or shrub cover in the southern populations than in the northern ones, in accordance with predictions from the stress-gradient hypothesis.
The specific questions addressed in this study were: considering results from small-scale studies, (1) at what scale does regeneration failure of P. nigra occur, and (2) what are the biotic and abiotic factors determining P. nigra regeneration at the regional scale? (3) Is P. nigra regeneration facilitated by tree canopies and/or shrubs, and, if the latter applies, what are the patterns of facilitation across the species distribution area?

2. Materials and Methods

2.1. Species, Data Set, and Study Variables

P. nigra thrives along the eastern calcareous mountains of Spain (Figure 1), under a Mediterranean mountain type climate characterized by cold-wet winters and hot-dry summers [26]. Seed predation, herbivory by ungulates, the presence of vegetation litter and the irregular production of seed crops are factors which negatively affect P. nigra regeneration in natural stands, although summer drought is the major “bottleneck” of P. nigra demography [27,28,30,31].
This study used data from the Third Spanish National Forest Inventory (3SNFI hereafter). The SNFI consists of permanent inventory plots systematically distributed on a one-square kilometer grid. Each plot is divided in four concentric subplots of 5, 10, 15, and 25 m radius. In these subplots, an adult tree is sampled if its diameter at breast height (d.b.h.) is 7.5–12.4, 12.5–22.4, 22.5–42.4, and ≥42.5 cm, respectively. Regeneration is measured separately for seedlings (height < 30 cm) and saplings (height: 30–130 cm) within the 5 m radius subplot, according to a semi-quantitative scale of abundance: absent (0 trees), low (1–4 trees), medium (5–15 trees), and high (>15 trees).
From the 3SNFI provincial databases, we selected inventory plots that (i) were situated within natural stands; (ii) contained at least one P. nigra individual with a d.b.h. ≥20 cm (in order to assure a seed source); (iii) and were located within the provinces of Lleida, Teruel, Cuenca, and Jaén-Granada. This selection of provinces allowed us to cover a latitudinal gradient of 650 km, including data from the mountainous areas with the most extensive forests of P. nigra (study regions hereafter): Prepirineo Catalán (CAT), Alto Maestrazgo-Montes Universales (MAE), Sistema Ibérico Meridional (IBE), and Cordilleras Béticas (BET) (Figure 1). A total of 3226 plots were selected by this filtering process.
Figure 1. Geographical distribution of Pinus nigra ssp. salzmannii in Spain with indication of the four study regions: Prepirineo Catalán (CAT) in yellow, Alto Maestrazago-Montes Universales (MAE) in red, Sistema Ibérico Meridional (IBE) in blue, and Cordilleras Béticas (BET) in green. The inset indicates the distribution area of Pinus nigra sensu lato.
Figure 1. Geographical distribution of Pinus nigra ssp. salzmannii in Spain with indication of the four study regions: Prepirineo Catalán (CAT) in yellow, Alto Maestrazago-Montes Universales (MAE) in red, Sistema Ibérico Meridional (IBE) in blue, and Cordilleras Béticas (BET) in green. The inset indicates the distribution area of Pinus nigra sensu lato.
Forests 05 00001 g001
The region, the abundance of P. nigra regeneration (seedlings and saplings), and the corresponding values of four stand-structural variables and six soil-physiographic variables were noted for the selected plots. The stand-structural variables were: canopy cover, shrub cover, basal area, and depth of the litterfall layer. Canopy cover, i.e., the proportion of the forest floor covered by the vertical projection of the tree crowns, is suitable for describing stand-level microclimate and can easily be estimated, thus, most forest inventories record it routinely [32,33]. Canopy and shrub covers are expected to be correlated in forest stands [34]. The soil and physiographic variables were: texture, content of organic matter, rockiness, slope, altitude, and aspect. The 3SNFI measured cover in percentage, basal area in m2 ha−1, litterfall depth in classes of 0.5 cm, and distinguishes among three classes of soil texture (1-sandy, 2-loam and 3-clay), three classes of organic matter content (1-high, 2-moderate and 3-low), and five classes of rockiness (from 1 to 5: 0%, 1%–25%, 26%–50%, 51%–75%, and >75%, respectively). The type of silvicultural treatment observed within the 25 m radius and the presence of damage caused by wild or domestic ungulates were also noted.
Mean values of temperature, precipitation, and potential evapotranspiration [35] of June, July, and August were used to characterize the climatic conditions at each plot during the summer, i.e., the most critical season for P. nigra seedlings survival [9]. Additionally, we calculated drought length, as the number of months during which the temperature curve (Ti) is above the precipitation curve (Pi) in the climograph (2Ti > Pi), and an Aridity Index as the quotient between precipitation and potential evapotranspiration during the mentioned months [36]. The climatic data were obtained for each inventory plot by an interpolation of the information recorded by the Spanish network of meteorological stations from 1971 to 2000, with 1 km spatial resolution [37]. We also obtained radiation data from [38], by uniformly distributing ten sampling points within each region.

2.3. Statistical Analysis

Both seedling and sapling abundance at each plot were used to address question (1), i.e., what is the scale of P. nigra regeneration failure, as seedlings and saplings will eventually contribute to maintain stocking in managed stands. However, we only considered data of P. nigra seedling abundance in order to address questions (2) and (3) as, working with the smallest and presumably youngest individuals, it could be assumed that the stand-structural variables recorded by the 3SNFI surveyors were equal or rather similar to those experienced by seedlings at the time of establishment (P. nigra surpasses 30 cm height in about six years [39]). Therefore, the stand-structural characteristics observed in plots where seedlings are present should be suitable descriptors of the P. nigra regeneration niche.
Data of seedling abundance from the 3226 selected plots were used to address questions (2) and (3) as follows. In a first step, we analyzed how seedling abundance varies across levels of canopy cover, using histograms and descriptive statistics. Results suggested that seedling establishment tends to fail in open areas, whilst it is more likely to occur in sites with moderate cover, and decrease in places of high cover (Figure 2). Consequently, we built 2 × 3 contingency tables to test whether regeneration success depends on the level of cover across the study regions: CAT, MAE, IBE, and BET. Plots were cross-classified according to the variables seedlings present / total absence of regeneration, and three levels of canopy cover: low (canopy cover < 55%), moderate (canopy cover: 55%–80%), and high (canopy cover > 80%). These levels were quantitatively established after examining the lower and upper quartiles from the distribution of canopy cover in plots where seedlings were present (Figure 2).
Figure 2. Frequency distribution of canopy cover in plots with seedlings present (n = 913). Arrows indicate the position of different percentiles.
Figure 2. Frequency distribution of canopy cover in plots with seedlings present (n = 913). Arrows indicate the position of different percentiles.
Forests 05 00001 g002
In a second step, we generated random numbers to select 170 plots for every combination of region and regeneration level; 170 being the minimum equal number of plots in any of the eight region-regeneration combinations (four regions: BET, IBE, MAE and CAT, and two regeneration levels: seedlings present and total absence of regeneration). This made a randomly selected sample of 1360 plots that was analyzed by Spearman correlation coefficients, Generalized Linear Models (GLM) and multivariate analyses. This selection process was done, as the design needed to be balanced in order to carry out a multivariate analysis of variance (see later). It could be argued that this way of imposing balance dropped part of the available data, but the method is statistically valid [40], and, still, we could analyze data from 136,000 ha of forest.
The relationship between seedling abundance and environmental variables was analyzed by a GLM, with a negative binomial distribution for the error term and log as the link function. A negative binomial distribution was preferred over a Poisson distribution as the response variable showed over-dispersion [41]. The predictor variables included in the full model were the stand-structural, soil, and physiographic variables listed earlier, plus the Aridity Index which integrated the observed variability in meteorological data (see later). Canopy cover and shrub cover were included as a linear and a second-order polynomial term to select the best transformation of these two explicative variables to account for non-linearity. The response variable was seedling abundance transformed into a continuous scale by adopting the lowest value of each class range as a conservative approach, i.e., 0, 1, 5, and 16 seedlings per plot for the 3SNFI classes absent, low, medium, and high, respectively [42]. We did not detect problems of multicollinearity in the fitted model (Variance Inflation Factor ≤ 1.93 for all the predictor variables). A stepwise procedure was run to select the most significant variables using the Akaike Information Criteria (stepAIC library MASS in R [43]).
Additionally, the effect of woody canopies on the presence or absence of P. nigra seedlings across the four study regions was analyzed by a multivariate approach. This option was chosen, because stand-structural variables can act concomitantly. For instance, canopy cover may facilitate seedling establishment, but accumulates litterfall on the ground, which, in turn, might lessen early seedling survival [9,34]. We used a permutational multivariate analysis of variance (PERMANOVA, hereafter). PERMANOVA gives a partitioning of multivariate variation, defined by a distance measure, according to individual factors in any fully balanced multi-way analysis of variance design, with tests done by permutation. No explicit assumptions regarding the distribution of the original variables need to be met [44,45]. In addition to PERMANOVA, a canonical discriminant analysis of principal coordinates (CAP) was conducted to test the hypothesis of no differences among the eight combinations of region and regeneration level, and represent them in a multidimensional space graphic. Specifically, the first two axes of the CAP were used to plot an unconstrained MDS (metric multidimensional scaling) of the data. The axis scores were correlated with the original variables using the Spearman correlation coefficient.
Finally, among-region differences for each study variable in plots with seedlings were tested by one-way analyses of variance (ANOVA), or the Kruskal-Wallis test when variables were in a nominal scale. The assumption of homocedasticity in ANOVA was not met for all the variables, even after data transformation. Nevertheless, ANOVAs were still carried out on raw data because they are robust to departures from the assumptions, if data are balanced (i.e., sizes of samples are all made the same) and samples are large [46].
The analysis were performed with the computer programs PERMANOVA and CAP [47,48], using 9999 permutations and the Gower dissimilarity as distance measure, because the dataset contained a mixture of continuous and categorical variables [49]. GLM and additional analysis were conducted with R 2.11.1 [43]. Throughout the paper, values are means ± standard errors.

3. Results

3.1. Climate and Soil Gradients

Results confirmed that the selected four study regions lay along a latitudinal gradient of summer aridity. Radiation did not differ among regions (F3,36 = 1.25; p = 0.306), but precipitation, mean temperature, potential evapotranspiration and the Aridity Index differed amongst regions for the months of June, July, and August (F3,1356 ≥ 110; p < 0.0001 for the four variables). Mean temperature and potential evapotranspiration increased, north to south, although the proportion of variability in these variables accounted for by the latitudinal gradient was low (R2 ≤ 0.24 for both variables), probably, due to an effect of altitude, which increases towards the south (Table 1). However, the latitudinal gradients for the Aridity Index and for the summer precipitation were clearer (R2 = 0.71 and 0.81, respectively). This north to south gradient of increasing aridity was amplified by the duration of summer drought: 1.3, 2.5, 2.9, and 3.1 months from the northernmost study region to the southernmost one (Figure 1). As a consequence, it is possible to state that P. nigra seedlings are more likely to experience summer drought in localities from the BET region, and less likely in localities from the CAT region. Soil texture was similar in all the study regions, but BET region exhibited the lowest organic matter content (F1,1358 = 16.59; p < 0.0001), and the highest rockiness (F1,1358 = 124.04; p < 0.0001). Indeed, rockiness increased gradually towards the south along the latitudinal gradient (β = −0.025; F1,1358 = 271.6; p < 0.001).

3.2. Scale of Regeneration Failure

P. nigra regeneration, either in the form of seedlings or saplings, was present at 52.9% of all study plots, with saplings being present in a higher proportion of plots in the four study regions (Figure 3). The percentage of plots lacking regeneration increased along the rockiness and aridity gradients (30%, 51%, 46% and 58% for the regions CAT, MAE, IBE, and BET, respectively). Multiple logistic regression showed that summer drought, expressed as the Aridity Index, and rockiness were negatively related to the presence of seedlings (df = 1; Wald ≥ 7.01; p ≤ 0.008). Figure 3 shows both a continuous and gradual decline of frequency from plots lacking regeneration towards those with a high density of seedlings or saplings, and a clear preponderance of plots lacking regeneration or belonging to the low class of regeneration abundance. A weak positive correlation was found between seedling and sapling abundances in the CAT, IBE, and BET regions (Spearman’s rho = 0.21–0.30; p < 0.003), and moderate in the MAE region (Spearman’s rho = 0.47; p < 0.0001), suggesting the existence of a great spatial-variability in the abundance-life stage combinations of regeneration, and that, generally, high densities of seedlings and saplings did not coincide in the same plots. As result, it could be reckoned that regeneration density is at best 509 trees per hectare (four trees/plot) in most forest-stands of the study area.
Table 1. Mean values and standard errors of the abiotic and stand-structural variables measured in 170 plots with presence of seedlings from each study region. Aridity Index = summer precipitation/potential evapotranspiration during summer. Lower values of the index represent drier conditions. The period June-July-August was considered as summer. Different letters between brackets indicate significant differences in a Tukey HDS test (differences were considered significant at p < 0.005 after Bonferroni corrections for 10 variables). *** Significant differences (p < 0.001) in a Kruskal-Wallis test.
Table 1. Mean values and standard errors of the abiotic and stand-structural variables measured in 170 plots with presence of seedlings from each study region. Aridity Index = summer precipitation/potential evapotranspiration during summer. Lower values of the index represent drier conditions. The period June-July-August was considered as summer. Different letters between brackets indicate significant differences in a Tukey HDS test (differences were considered significant at p < 0.005 after Bonferroni corrections for 10 variables). *** Significant differences (p < 0.001) in a Kruskal-Wallis test.
VariablesStudy Region
CAT MAE IBEBET
Summer precipitation (mm)187 ± 1.8(a)134 ± 0.8(b)103 ± 0.8(c)71 ± 0.7(d)
Mean summer temperature (°)19.7 ± 0.1(a)18.2 ± 0.1(b)19.7 ± 0.1(a)20.0 ± 0.1(a)
Aridity Index0.51 ± 0.01(a)0.41 ± 0.00(b)0.29 ± 0.00(c)0.20 ± 0.00(d)
Altitude (m)828 ± 15(d)1269 ± 17(b)1174 ± 12(c)1405 ± 13(a)
Slope (°)9.19 ± 0.02(a)7.45 ± 0.21(b)5.36 ± 0.21(c)8.58 ± 0.20(a)
Rockiness1.73 ± 0.04(d)2.70 ± 0.05(b)2.36 ± 0.05(c)2.96 ± 0.04(a)
Soil texture ***2.26 ± 0.032.14 ± 0.022.26 ± 0.032.21 ± 0.03
Organic Matter content ***2.11 ± 0.022.06 ± 0.022.10 ± 0.031.89 ± 0.03
Canopy cover (%)69.9 ± 0.9(a)60.6 ± 0.9(b)58.2 ± 0.9(b)61.0 ± 1.1(b)
Shrub cover (%)38.6 ± 1.4(b)37.6 ± 1.5(b)51.7 ± 1.7(a)23.6 ± 1.2(c)
Basal Area (m2 ha−1)17.6 ± 0.5(bc)15.8 ± 0.6(cd)13.4 ± 0.6(d)19.3 ± 0.6(ab)
Litter depth (cm)3.09 ± 0.10(a)1.50 ± 0.06(c)2.04 ± 0.10(b)2.01 ± 0.08(b)
Figure 3. Percentage of plots in the four study regions containing 0 individuals, 1–4 individuals, 5–15 individuals, and more than 15 individuals of Pinus nigra subsp. salzmannii in each of two size classes: (a) seedlings (height < 30 cm), and (b) saplings (height: 30–130 cm). n = 3226 plots.
Figure 3. Percentage of plots in the four study regions containing 0 individuals, 1–4 individuals, 5–15 individuals, and more than 15 individuals of Pinus nigra subsp. salzmannii in each of two size classes: (a) seedlings (height < 30 cm), and (b) saplings (height: 30–130 cm). n = 3226 plots.
Forests 05 00001 g003

3.3. Biotic and Abiotic Factors Determining Regeneration

Spearman rank correlation analysis revealed significant univariate relationships between regeneration abundance and several stand-structural, physiographic and climatic variables (Table 2). Table 3 summarizes the frequencies by which seedlings are present or absent across three levels of canopy cover. The observed differences were statistically significant (χ2 > 8.77; p < 0.004 for all the study regions). Thus, seedlings of P. nigra were more frequent than expected by random in the plots of moderate and high canopy cover and less frequent than expected in the plots of low canopy cover, with the exception of the CAT region, where seedlings were more frequent than expected by random at plots of moderate canopy cover only.
Table 2. Spearman rank correlations between biotic and abiotic variables, and seedling density across the four study regions; significant results are noted in bold.
Table 2. Spearman rank correlations between biotic and abiotic variables, and seedling density across the four study regions; significant results are noted in bold.
VariablesRegion
CATMAEIBEBETAll
Canopy cover−0.0500.265 ***0.206 **0.1050.125 ***
Shrub cover0.066−0.067−0.0040.391 ***0.084 *
Stand basal area−0.1300.244 ***0.256 ***0.194 **0.145 ***
Litter depth−0.0420.1380.170 *−0.0470.041
Rockiness−0.0080.201 **−0.1160.178 **0.109 ***
Soil texture0.114−0.001−0.0250.143 *−0.018
Organic Matter0.067−0.0190.1010.0520.050
Aspect0.010−0.0490.0570.0490.017
Slope−0.0190.009−0.0220.006−0.009
Aridity Index−0.114−0.016−0.106−0.006−0.019
* p < 0.05; ** p < 0.01; *** p < 0.001.
Table 3. Contingency table of the seedling frequencies of Pinus nigra subsp. salzmannii in the four study regions. Low (canopy cover < 55%); Moderate (55% ≤ canopy cover ≤ 80%); High (canopy cover > 80%).
Table 3. Contingency table of the seedling frequencies of Pinus nigra subsp. salzmannii in the four study regions. Low (canopy cover < 55%); Moderate (55% ≤ canopy cover ≤ 80%); High (canopy cover > 80%).
RegionSeedlingsClasses of canopy cover
LowModerateHighTotal
CATPresent2114742210
Absent104343114561
MAEPresent3212415171
Absent14419212348
IBEPresent9823231361
Absent40944840897
BETPresent5110119171
Absent19329024507
Total 105218772973226
Canopy cover explained most of the deviance in the GLM (Table 4). Basal area and shrub cover were also significant explanatory factors of P. nigra regeneration, although the deviance explained by canopy cover was four times larger than the deviance explained by shrub cover (canopy cover and shrub cover were both structural variables that might facilitate seedling establishment). Rockiness and Aridity Index were the only two abiotic factors that influenced P. nigra regeneration. In terms of explained deviance, the contribution of rockiness was the second most important. Coefficient estimates indicated a positive effect of canopy cover, shrub cover and basal area on regeneration, while the effects of rockiness and Aridity Index were negative. For canopy cover and shrub cover a model including the second-order polynomial transformations was selected, indicating a non-linear response of P. nigra recruitment to these two structural variables. The parameter of the second-order term was negative for both canopy cover and shrub cover, suggesting that high values of cover exert a negative effect on seedling abundance, and that seedling abundance might be enhanced at intermediate values of cover (Figure 4). This interpretation was further supported by negatively skewed frequency distributions of canopy cover in plots with presence of seedlings (Figure 5).
Table 4. Best model (Generalized Linear Model with a log link and a negative-binomial distribution) of Pinus nigra subsp. salzmannii seedling abundance estimated by a stepwise procedure. The selected transformations poly (2) are polynomial second-order transformations, degree of freedom (df), residual degree of freedom (Resid. df), residual deviance (Resid. Dev.) and probabilities χ2 tests of the effect of the variable are given. Model deviance = 0.064.
Table 4. Best model (Generalized Linear Model with a log link and a negative-binomial distribution) of Pinus nigra subsp. salzmannii seedling abundance estimated by a stepwise procedure. The selected transformations poly (2) are polynomial second-order transformations, degree of freedom (df), residual degree of freedom (Resid. df), residual deviance (Resid. Dev.) and probabilities χ2 tests of the effect of the variable are given. Model deviance = 0.064.
dfDevianceResid. dfResid. Dev.P(χ2 test)
Null 13591282.4
Canopy cover poly (2)242.9113571239.50.00000
Rockiness113.4813531205.80.00024
Shrub cover poly (2)210.3913551229.10.00553
Basal Area19.8513541219.20.00170
Aridity15.1813521200.60.02284
The sampled plots included a list of 85 shrub species. Among them, the highest percentage of shrub cover corresponded to Mediterranean pioneer plants: genus Thymus, Rosmarinus, and Lavandula, followed by legumes: genus Genista, Erinacea, and Dorycnium, and spiny shrubs: genus Rubus, and Rosa. Respectively, these three groups represented 47%, 23%, and 9% of the overall percentage of shrub cover. An additional GLM, carried out to study the relationship between seedling abundance and cover of the mentioned shrub genuses, showed that facilitation effects are species-specific, being negative for Thymus (df = 1; deviance = 8.99; P(χ2) = 0.002), positive for Lavandula (df = 1; deviance = 5.36; P(χ2) = 0.021) and Rosa (df = 1; deviance = 5.82; P(χ2) = 0.015), and not significant for the other genuses (df = 1; deviances ≤ 2.45; P(χ2) ≥ 0.118). Significant interactions were detected between the region of origin and canopy cover (df = 4; deviance = 38.26; P(χ2) < 0.0001) and shrub cover (df = 4; deviance = 17.89; P(χ2) < 0.01). Seedling abundance was positively affected by canopy cover in the four study regions (Z values ≥ 4.95; P(χ2) < 0.0001). However, the positive effect of shrub cover was limited to the BET region (Z value = 4.06; <0.0001), being not significant for the other three regions (Z values ≤ 1.85; P(χ2) ≥ 0.065). A closer examination of data revealed that legumes were as frequent as pioneer plants in the BET area, and that cover of Erinacea was positively correlated to seedling abundance in this region (Spearman’s rho = 0.212; p < 0.001).
The silvicultural treatment applied in 26% of the study plots was selective felling, whilst non-silvicultural treatment was observed in 72% of the plots. Nevertheless, selective felling must be widely implemented in P. nigra forests, as most plots were classified as belonging to multicohort stands by SNFI surveyors (72%, 61%, 79%, and 66% for the regions CAT, MAE, IBE, and BET, respectively). No damage caused by wild or domestic ungulates was observed in the selected plots.

3.4. Regional Patterns of Stand-Structural Effects on Regeneration

PERMANOVA detected regional differences in the structural variables considered simultaneously (F = 44.29, p < 0.0001). The structural characteristics of the plots occupied by seedlings differed from those with no regeneration (F = 12.08, p < 0.001), and the interaction region × regeneration density was significant as well (F = 9.26, p < 0.0001). CAP analysis run to test the null hypothesis of no differences among the groups resulting from the region x regeneration interaction further confirmed this result (trace statistic = 0.317, first squared canonical correlation = 0.154; p < 0.0001 for both tests). The first two axes from CAP explained over 87% of the variation observed in the data. The first axis of this ordination showed a strong negative correlation with canopy cover, basal area, and litterfall depth, while the second axis showed a strong negative correlation with shrub cover. These results indicated the existence of gradients along axes 1 and 2, where sites with higher values of canopy and shrub cover were placed at the bottom left side of the graph. In this respect, groups with presence of seedlings tended to dispose themselves around intermediate values of cover and related structural variables, whilst groups representing sites with absence of regeneration exhibited a more disperse pattern (Figure 4). Within regions, seedlings were found to establish in sites with higher values of canopy cover (first CAP axis) in the regions located at mid position along the aridity gradient (MAE and IBE). On the contrary, seedling presence was higher at lower values of canopy cover in the less dry region (CAT), and no differences with regards the effect of canopy cover on the presence or absence of seedlings was found in the driest region (BET).
Figure 4. Canonical discriminant analysis of principal coordinates of the stand-structural variables characterizing the sites where P. nigra subsp. salzmannii seedlings are present (pre) or absent (abs) in the four study regions (CAT, MAE, IBE, and BET). Data are means ± standard errors (n = 1360).
Figure 4. Canonical discriminant analysis of principal coordinates of the stand-structural variables characterizing the sites where P. nigra subsp. salzmannii seedlings are present (pre) or absent (abs) in the four study regions (CAT, MAE, IBE, and BET). Data are means ± standard errors (n = 1360).
Forests 05 00001 g004
Mean values of canopy cover, shrub cover, basal area, and litterfall depth differed among regions (F3,1356 ≥ 18.58, p ≤ 0.0001 for the three variables). Thus, CAT region exhibited the highest mean values of canopy cover and litterfall depth, IBE region the highest value of shrub cover, and CAT and BET regions the highest values of basal area (Table 1).

4. Discussion

4.1. Scale of Regeneration Failure

Our study revealed that regeneration failure of P. nigra might be a common phenomenon across the species distribution area. Nearly half of the plots lacked seedlings or saplings, although the presence of juvenile trees would generally be expected within the observed multi-cohort structure of P. nigra stands. Obviously, the reported patterns of recruitment were influenced by the sample effort applied, just 78 m2 within a larger plot of 1963 m2. Increasing the surveyed area of regeneration in each plot may have provided different recruitment patterns. However, both the area analyzed and the sample size were large enough to accomplish the aims of the study. Thus, results showed a clear decline of regeneration success along the rockiness and summer drought gradients. Seedlings and/or saplings were present in 70% of the plots from the northern region (CAT), but the percentage decreased steadily towards the driest southern edge. Only 42% of the plots from the BET region attained new trees. [42] found significant effects of rockiness on P. nigra recruitment, but, to our knowledge, no other work has previously considered rockiness as an important variable affecting recruitment success in Iberian pine forests. Here, we have found that regeneration is less frequent in plots from the drier and rockier southern edge of P. nigra natural forests, although data did not let us clearly disentangle the contribution of each variable. Other studies have mentioned precipitation as a factor positively correlated with significant increases in recruitment success [42,50]. The regeneration pattern observed at a regional scale in this study is then consistent with results from small-scale studies, which identified drought as the main cause of seedling mortality in rock-free P. nigra stands [6,9,27]. Consequently, summer drought might be the major limiting factor of recruitment throughout the species distribution area. This impression was further sustained by the facts that, first, overgrazing does not seem to be responsible for the lack of P. nigra regeneration, although herbivore pressure is characteristically high in Mediterranean forests [8,50]. Second, although data did not show significant evidences of recent silvicultural treatments (a disturbance that usually triggers regeneration in managed stands), values of stand basal area were within those more suitable for the establishment of new trees [51,52].
Finally, it must be noticed that, where seedlings and saplings are present, current density may be insufficient to restock uneven-aged P. nigra stands. Specifically, the reverse-J-shaped diameter frequency distributions proposed for the species management would need as many as 1600 seedlings per hectare to sustain harvest operations [53]. Likewise, at least 2000 seedlings per hectare are needed to restock even-aged stands [51], but we reckoned that current regeneration density is at best 509 trees per hectare. Concerns about the long-term persistence of P. nigra forests under current climate change seem to be justified [26].

4.2. Biotic and Abiotic Factors Determining Regeneration

While results indicated that summer drought affects the presence or absence of P. nigra juveniles across the species distribution area, summer drought did not seem to produce any effect upon the density of regeneration where it occurs, as indicated by a non-significant Spearman’s correlation coefficient between seedling density and the Aridity Index (Table 2). GLM further confirmed that the density of regeneration is, if at all, only slightly dependent on summer aridity (Table 4). Accordingly, plots with a high density of seedlings and/or saplings were not more frequent at the northern sites than in the southern ones (Figure 3).
Low or non-significant Spearman’s r-values, and low explained variances showed that occurrence of regeneration scarcely varies over topographic, and edaphic site factors in the P. nigra distribution area, with the exception of rockiness (Table 2 and Table 4). The presence of rocks was associated to poor-sandy soils, steep slopes, and open canopies (data not shown). Therefore, and beyond the fact that rocks diminish the surface available for seed germination, the negative correlation found between seedling density and rockiness could be a consequence of water stress due to a combination of both inferior edaphic conditions, and high evaporation under low canopy cover [8]. Canopy cover consequently appeared as the main ecological factor determining seedling establishment in P. nigra stands.

4.3. Patterns of Facilitation across P. nigra Distribution Area

Our results indicated that P. nigra juvenile trees can adapt to shady conditions, most likely in order to improve their water balance. Seedlings of this species exhibit no photo-inhibition and a high photosynthetic capacity in full sunlight [54], but recruitment proved to be positively associated with tree cover under natural conditions and across the species distribution area, with the highest frequency of seedling establishment being met at mid canopy cover (between 55% and 80%). Within each study region, the frequency distributions of canopy cover in plots with regeneration exhibited lower coefficients of variation and more negative skew than the plots with no regeneration (Figure 5). Moreover, the environment described by a combination of four stand-structural variables involved in P. nigra recruitment was quite similar for plots with seedling presence across the species distribution area (Figure 4). These findings suggest that cover from the tree canopy is a large-scale habitat condition of the P. nigra regeneration niche, similar to climate or soil conditions. The observation that young P. nigra trees benefit from the special abiotic conditions beneath tree canopies suggests the existence of facilitative interactions between seedlings and canopy trees [6]. Indeed, Spanish foresters have been largely aware of the shade tolerance of P. nigra seedlings, and of the importance of tree cover for P. nigra regeneration [51,53], but this study offers insights into the nature of such facilitative interaction.
Figure 5. Frequency distribution and descriptive statistics of canopy cover in sites where seedlings were present (A) or absent (B) within each study region. SD is the standard deviation; CV is the coefficient of variation.
Figure 5. Frequency distribution and descriptive statistics of canopy cover in sites where seedlings were present (A) or absent (B) within each study region. SD is the standard deviation; CV is the coefficient of variation.
Forests 05 00001 g005
Percentage cover in plots with regeneration did not increase towards the South as we had previously hypothesized based on predictions from the stress-gradient hypothesis. On the contrary, evidence from this study suggests that the interacting effect of light availability and summer drought on P. nigra survival at the seedling stage is nonlinear (humped-back shaped [24]) at both regional and local scales. Overall a humped-back shaped pattern of facilitation was indicated by the inclusion of canopy cover as polynomial of order 2 in the best GLM (Table 4), and by the negative skewness (−0.72 ± 0.07) calculated for the frequency distribution of canopy cover in plots with seedlings present (Figure 2). Within regions, results similarly showed more recruitment at intermediate values of tree cover than at both ends of the canopy gradient (Table 3 and Figure 5).
Most studies about the recruitment dynamics of Mediterranean tree species have identified drought as the major cause of seedling mortality [1,2,3,4], including studies on P. nigra regeneration [6,9]. High irradiance increases evaporative demands and further reduces soil-water availability, therefore, there are obvious reasons to explain the increased seedling mortality observed in open areas [5,6,7,8]. In poor light environments, the eventual death of P. nigra seedlings because of desiccation can be justified by at least two reasons. First, soil water content may be lower than expected under tree canopies because of rainfall interception and, second, a shift from facilitation to competition for water might occur above a given threshold of canopy cover (~80% in this study) [55]. In both cases, seedlings would die in dense stands, because water absorption and light capture cannot be maximized simultaneously, as predicted by the trade-off hypothesis [21] and demonstrated in greenhouse experiments by [54,56] Both reasons could explain why plots with seedling presence from the wetter northern CAT region attained the highest percentages of canopy cover.
However, the trade-off hypothesis [21] and the humped-back shape model [24] proposed a switch from facilitation to competition at both ends of an aridity gradient. This switch was observed in this study in the wettest region but not in the driest one. Specifically, results suggest the existence of competition between seedlings and grown trees in the CAT region: seedling establishment was not more frequent than expected under the more dense tree canopies (Table 3), and regeneration was absent in plots that attained the highest values of canopy cover (Figure 4). On the contrary, sites with and without P. nigra regeneration attained similar values of canopy cover in the BET region, which may suggest a collapse of facilitation in this area with harsher environmental conditions for seedling establishment [25,57]. Indeed, a clear facilitative effect of tree canopy on regeneration was detected in the regions located at mid position on the aridity gradient (MAE and IBE), as shown by the distances between groups in the unconstrained MDS of these two regions, representing sites where seedlings are present or absent (Figure 4).
The facilitative effect of shrubs on P. nigra regeneration was not generally extended or strong at the regional scale. Yet, some studies on pines have reported a nurse effect of shrubs. For instance, [6] observed facilitative effects of Juniperus communis on P. nigra regeneration and two sowing experiments carried out in locations within the P. nigra distribution area revealed facilitative effects from shrubs on the seedling survival of Pinus sylvestris, a pine species that also inhabits Mediterranean mountains [4,10]. Nevertheless, shrub-tree seedling facilitative effects tend to be species-specific [58], and our results are congruent with the poor nurse quality quoted for the most abundant shrub species in the whole study area [59]. Thus, BET was the only region where facilitative effects by shrubs were observed because legume-shrub species were abundant, while the other regions did not benefit from the poor nurse quality of the more frequent Mediterranean pioneer species.

5. Conclusions

This study identifies that P. nigra may have recruited insufficiently during the last decades in order to restock managed stands, this circumstance being more pronounced in the drier southern edge of the species distribution area. The canopy of adult trees facilitates seedling establishment, preferably within an interval of 55%–80% canopy cover. Results suggests that this facilitative interaction is non-linear (humped-back shaped) at the local (within region) and regional (across regions) scales, although this pattern interacts along the aridity gradient showing competition in the wettest site, high facilitation in the mid-dry sites, and low facilitation at the driest site. This pattern resembles predictions from the stress-gradient hypothesis, but with a collapse of facilitation in the most water-stressed environment. We detected these patterns by analysing data along a continuous aridity gradient (GLM), and among regions (contingency tables and multivariate analysis). Evidence suggests that some shrub species may facilitate recruitment in the harsher areas. Collectively, our results reduce the possibilities of adapting forests management to drought by the application of alternative silvicultural prescriptions involving canopy cover. Foresters working with Mediterranean trees at the southern limit of their distribution area will need to know as much about shrub ecology as they already know about the ecology of commercial tree species.

Acknowledgments

Consejería de Medio Ambiente (Junta de Andalucía) facilitated the means to carry out the research. Ruiz-Benito P. facilitated access to climatic data. We are grateful to three anonymous referees for comments on an earlier draft of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Herrera, C.M.; Jordano, P.; López-Soria, L.; Amat, J.A. Recruitment of a mast-fruiting, bird-dispersed tree: Bridging frugivore activity and seedling establishment. Ecol. Monogr. 1994, 64, 315–344. [Google Scholar] [CrossRef] [Green Version]
  2. Castro, J.; Zamora, R.; Hódar, J.A.; Gómez, J.M. Seedling establishment of a boreal tree species (Pinus sylvestris) at its southernmost distribution limit: Consequences of being in a marginal Mediterranean habitat. J. Ecol. 2004, 92, 266–277. [Google Scholar] [CrossRef]
  3. Rey, P.J.; Alcántara, J. Recruitment dynamics of a fleshy-fruited plant (Olea europaea): Connecting patterns of seed dispersal to seedling establishment. J. Ecol. 2000, 88, 622–633. [Google Scholar] [CrossRef]
  4. Castro, J.; Zamora, R.; Hódar, J.A.; Gómez, J.M. Alleviation of summer drought boosts establishment success of Pinus sylvestris in a Mediterranean mountain: An experimental approach. Plant Ecol. 2005, 181, 191–202. [Google Scholar] [CrossRef]
  5. Granda, E.; Escudero, A.; de la Cruz, M.; Valladares, F. Juvenile-adult tree associations in a continental Mediterranean ecosystem: No evidence for sustained and general facilitation at increased aridity. J. Veg. Sci. 2012, 23, 164–175. [Google Scholar] [CrossRef] [Green Version]
  6. Tíscar, P.A.; Linares, J.C. Structure and regeneration patterns of Pinus nigra ssp. salzmannii natural forests: A basic knowledge for adaptive management in a changing climate. Forests 2011, 4, 1013–1030. [Google Scholar]
  7. Valladares, F.; Zaragoza-Castells, J.; Sánchez-Gómez, D.; Matesanz, S.; Alonso, B.; Portsmuth, A. Is shade beneficial for Mediterranean shrubs experiencing periods of extreme drought and late-winter frosts? Ann. Bot. 2008, 102, 923–933. [Google Scholar] [CrossRef]
  8. Gómez-Aparicio, L.; Zamora, R.; Castro, J.; Hódar, J.A. Facilitation of tree saplings by nurse plants: Microhabitat amelioration or protection against herbivores? J. Veg. Sci. 2008, 19, 161–172. [Google Scholar] [CrossRef]
  9. Tíscar, P.A. Regeneration dynamics of Pinus nigra subsp. salzmannii in its southern distribution area: Stages, processes and implied factors. For. Syst. 2007, 16, 124–135. [Google Scholar]
  10. Mendoza, I.; Zamora, R.; Castro, J. A seeding experiment for testing tree-community recruitment under variable environments: Implications for forest regeneration and conservation in Mediterranean habitats. Biol. Conserv. 2009, 149, 1491–1499. [Google Scholar] [CrossRef]
  11. Lookingbill, T.; Zavala, M.A. Spatial pattern of Quercus ilex and Quercus pubescens recruitment in Pinus halepensis dominated woodlands. J. Veg. Sci. 2000, 11, 607–612. [Google Scholar] [CrossRef]
  12. Callaway, R.M. Positive interactions among plants. Bot Rev. 1995, 61, 306–349. [Google Scholar] [CrossRef]
  13. Bertness, M.D.; Callaway, R.M. Positive interactions in communities. Trends Ecol. Evol. 1994, 9, 191–193. [Google Scholar] [CrossRef]
  14. Gómez-Aparicio, L.; Zamora, R.; Gómez, J.M.; Hódar, J.A.; Castro, J.; Baraza, E. Applying plant facilitation to forest restoration: A meta-analysis of the use of shrubs as nurse plants. Ecol. Appl. 2004, 14, 1128–1238. [Google Scholar]
  15. Sthultz, C.M.; Gehring, C.A.; Whitham, T.G. Shifts from competition to facilitation between a foundation tree and a pioneer shrub across spatial and temporal scales in a semiarid woodland. New Phytol. 2007, 173, 135–145. [Google Scholar] [CrossRef]
  16. Lucas-Borja, M.L.; Fidalgo, T.; Lousada, J.L.; Silva-Santos, P.; Martínez García, E.; Andrés, M. Natural regeneration of Spanish black pine [Pinus nigra Arn. ssp. salzmannii (Dunal) Franco] at contrasting altitudes in a Mediterranean mountain area. Ecol. Res. 2012, 27, 913–921. [Google Scholar]
  17. Maestre, F.T.; Valladares, F.; Reynolds, J.F. Is the change of plant-plant interactions with abiotic stress predictable? A meta-analysis of field results in arid environments. J. Ecol. 2005, 93, 748–757. [Google Scholar]
  18. Maestre, F.T.; Cortina, J. Are P.halepensis plantations useful as a restoration tool in semiarid Mediterranean areas? For. Ecol. Manag. 2004, 198, 303–317. [Google Scholar] [CrossRef]
  19. Castro, J.; Zamora, R.; Gómez, J.M.; Hódar, J.A.; Gómez-Aparicio, L. Benefits of using shrubs as nurse plants for reforestation in Mediterranean mountains: A 4-year study. Restor. Ecol. 2004, 12, 352–358. [Google Scholar] [CrossRef] [Green Version]
  20. Verdú, M.; Gómez-Aparicio, L.; Valiente-Banuet, A. Phylogenetic relatedness as a tool in restoration ecology: A meta-analysis. Proc. R. Soc. B. 2011. [Google Scholar] [CrossRef] [Green Version]
  21. Smith, T.; Huston, M. A theory of the spatial and temporal dynamics of plant communities. Vegetation 1989, 83, 49–69. [Google Scholar]
  22. Aranda, I.; Pardos, M.; Puértoloas, J.; Jiménez, M.D.; Pardos, J.A. Water-use efficiency in cork oak (Quercus suber) is modified by the interaction of water and light availabilities. Tree Physiol. 2007, 27, 671–677. [Google Scholar] [CrossRef]
  23. Valladares, F.; Pearcy, R.W. Drought can be more critical in the shade than in the sun: A field study of carbon gain and photoinhibition in a Californian shrub during a dry El Niño year. Plant Cell Environ. 2002, 25, 749. [Google Scholar] [CrossRef]
  24. Holmgren, M.; Gómez-Aparicio, L.; Quero, J.L.; Valladares, F. Non-linear effects of drought under shade: Reconciling physiological and ecological models in plant communities. Oecologia 2012, 169, 293–305. [Google Scholar] [CrossRef] [Green Version]
  25. Michalet, R.; Broker, R.W.; Cavieres, L.A.; Kikvidze, Z.; Lortie, C.J.; Pugnaire, F.I.; Valiente-Banuet, A.; Callaway, R.M. Do biotic interactions shape both sides of the humped-back model of species richness in plant communities? Ecol. Lett. 2006, 9, 767–773. [Google Scholar]
  26. Linares, J.C.; Tíscar, P.A. Climate change impacts and vulnerability of the southern populations of Pinus nigra subsp. salzmannii. Tree Physiol. 2010, 30, 795–806. [Google Scholar] [CrossRef]
  27. Lucas-Borja, M.E.; Fonseca, T.; Parresol, B.R.; Silva-Santos, P.; García-Morote, F.A.; Tíscar, P. Modelling Spanish black pine seedling emergence: Establishing management strategies for endangered forest areas. For. Ecol. Manag. 2011, 262, 195–202. [Google Scholar] [CrossRef]
  28. Tíscar, P.A.; Lucas-Borja, M.E.; Candel-Pérez, D. Changes in the structure and composition of two Pinus nigra subsp. salzmannii forests over a century of different silvicultural treatments. For. Syst. 2011, 3, 525–535. [Google Scholar]
  29. Coll, M.; Peñuelas, J.; Ninyerola, M.; Pons, X.; Carnicer, J. Multivariate effect gradients driving forest demographic responses in the Iberian Peninsula. For. Ecol. Manag. 2013, 303, 195–209. [Google Scholar]
  30. Cerro, A.; Lucas-Borja, M.E.; Martínez, E.; López, F.R.; Andrés, M.; García, F.A.; Navarro, R. Influence of stand density and soil treatment on the Spanish Black Pine (Pinus nigra Arn. ssp. salzmannii) regeneration in Spain. Investi. Agro. Sistemas y Recursos Forest. 2009, 18, 167–180. [Google Scholar]
  31. Tíscar, P.A. Conditions and limits to the natural regeneration of a mountain pine forest (Pinus nigra subsp. salzmannii). For. Syst. 2003, 12, 55–64. [Google Scholar]
  32. Jennings, S.B.; Brown, N.D.; Sheil, D. Assessing forest canopies and understory illumination: Canopy closure, canopy cover and other measures. Forestry 1999, 72, 59–73. [Google Scholar]
  33. Marchi, A.; Paletto, A. Relationship between forest canopy and natural regeneration in the subalpine spruce-larch forest (north-east Italy). Folia For. Pol. 2010, 52, 3–12. [Google Scholar]
  34. Castedo-Dorado, F.; Gómez-Vázquez, I.; Fernandes, P.M.; Crecente-Campo, F. Shrub fuel characteristics estimated from overstory variables in NW Spain pine stands. For. Ecol. Manag. 2012, 275, 130–141. [Google Scholar]
  35. Thornthwaite, C.W. An approach toward a rational classification of climate. Geogr. Rev. 1948, 38, 55–94. [Google Scholar] [CrossRef]
  36. United Nations Environment Programme (UNEP). World Atlas of Desertification, 2nd ed.; United Nations Environment Programme: London, UK, 1997. [Google Scholar]
  37. Gonzalo, J. Diagnosis Fitoclimática de la España Peninsular.CD con Bases de Datos Climáticas y Fitoclimaticos; Ministerio de Medio Ambiente, Rural y Marino.: Madrid, Spain, 2010. [Google Scholar]
  38. Ninyerola, M.; Pons, X.; Roure, J.M. Atlas Climático Digital de la Península Ibérica. Metodología y Aplicaciones en Bioclimatología y Geobotánica. Available online: http://opengis.uab.es/wms/iberia/index.htm (accessed on 1 July 2011).
  39. Tíscar, P.A. Estructura, Regeneración y Crecimiento de Pinus Nigra en el Área de Reserva Navahondona-Guadahornillos (Sierra de Cazorla, Jaén). Ph.D. Thesis, Polytechnic University of Madrid, Madrid, Spain, 2004. [Google Scholar]
  40. Shaw, R.G.; Mitchell-Olds, T. Anova for unbalanced data: An overview. Ecology 1993, 74, 1638–1645. [Google Scholar]
  41. Bolker, B. Ecological Models and Data in R; Princenton University Press: Princenton, NJ, USA, 2008. [Google Scholar]
  42. Ruiz-Benito, P.; Gómez-Aparicio, L.; Zavala, M.A. Large-scale assessment of regeneration and diversity in Mediterranean planted pine forests along ecological gradients. Divers. Distribut. 2012, 18, 1092–1106. [Google Scholar]
  43. R Development Core Team. R: A Language and Environment for Statistical Computing; The R Foundation for Statistical Computing: Vienna, Austria, 2013. [Google Scholar]
  44. Anderson, M.J. A new method for non-parametric multivariate analysis of variance. Aust. Ecol. 2001, 26, 32–46. [Google Scholar]
  45. McArdle, B.H.; Anderson, M.J. Fitting multivariate models to community data: A comment on distance-based redundancy analysis. Ecology 2001, 82, 290–297. [Google Scholar] [CrossRef]
  46. Underwood, A.J. Experiments in Ecology. Their Logical Design and Interpretation Using Analysis of Variance; Cambridge University Press: Cambridge, UK, 1997; p. 522. [Google Scholar]
  47. Anderson, M.J.; Robinson, J. Generalised discriminant analysis based on distances. Aust. N. Zeal. J. Stat. 2003, 45, 301–318. [Google Scholar]
  48. Anderson, M.J.; Willis, T.J. Canonical analysis of principal coordinates: A useful method of constrained ordination for ecology. Ecology 2003, 84, 511–552. [Google Scholar] [CrossRef]
  49. Gower, J.C. A general coefficient of similarity and some of its properties. Biometrics 1971, 27, 857–871. [Google Scholar] [CrossRef]
  50. Linares, J.C. Shifting limiting factors for population dynamics and conservation status of the endangered English yew (Taxus baccata L., Taxaceae). For. Ecol. Manag. 2013, 291, 119–127. [Google Scholar] [CrossRef]
  51. Alejano, R.; González, J.M.; Serrada, R. Selvicultura de Pinus nigra Arn. subsp. salzmannii (Dunal) Franco. In Compendio de Selvicultura Aplicada en España; Serrada, R., Montero, G., Reque, J.A., Eds.; INIA: Madrid, Spain, 2008; pp. 313–356. [Google Scholar]
  52. Lucas-Borja, M.E.; Fonseca, T.; Linares, J.C.; García-Morote, F.A.; López-Serrano, F.R. Does the recruitment pattern of Spanish black pine (Pinus nigra Arn ssp. salzmannii) change the regeneration niche over the early life cycle of individuals? For. Ecol. Manag. 2012, 284, 93–99. [Google Scholar] [CrossRef]
  53. González, J.M. Modelos de transformación de masa regular a irregular. For. Syst. 2000, 9, 237–252. [Google Scholar]
  54. Gómez-Aparicio, L.; Valladares, F.; Zamora, R. Differential light responses of Mediterranean tree saplings: Linking ecophysiology with regeneration niche in four co-occurring species. Tree Physiol. 2006, 26, 947–958. [Google Scholar] [CrossRef]
  55. Gómez-Aparicio, L.; Zavala, M.A.; Bonet, F.J.; Zamora., R. Are Pine plantations valid tools for restoring Mediterranean forests? An assessment along abiotic and biotic gradients. Ecol. Appl. 2009, 19, 2124–2141. [Google Scholar] [CrossRef] [Green Version]
  56. Tíscar, P.A.; Lucas-Borja, M.E. Seed mass variation, germination time and seedling performance in a population of Pinus nigra subsp. salzmannii. For. Syst. 2010, 19, 344–353. [Google Scholar]
  57. Michalet, R.; Bagousse-Pinguet, Y.L.; Maalouf, J.P.; Lortie, C.J. Two alternatives to the stress-gradient hypothesis at the edge of life: The collapse of facilitation and the switch from facilitation to competition. J. Veg. Sci. 2013. [Google Scholar] [CrossRef]
  58. Plieninger, T.; Rolo, V.; Moreno, G. Large-scale patterns of Quercus ilex, Quercus Suber, and Quercus pyrenaica regeneration in central-western Spain. Ecosystems 2010, 13, 644–660. [Google Scholar] [CrossRef]
  59. Siles, G.; Rey, P.J.; Alcántara, J.M.; Ramírez, J.M. Assessing the long-term contribution of nurse plants to restoration of Mediterranean forests through Markovian models. J. Appl. Ecol. 2008, 45, 1790–1798. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Tíscar, P.A.; Linares, J.C. Large-Scale Regeneration Patterns of Pinus nigra Subsp. salzmannii: Poor Evidence of Increasing Facilitation Across a Drought Gradient. Forests 2014, 5, 1-20. https://0-doi-org.brum.beds.ac.uk/10.3390/f5010001

AMA Style

Tíscar PA, Linares JC. Large-Scale Regeneration Patterns of Pinus nigra Subsp. salzmannii: Poor Evidence of Increasing Facilitation Across a Drought Gradient. Forests. 2014; 5(1):1-20. https://0-doi-org.brum.beds.ac.uk/10.3390/f5010001

Chicago/Turabian Style

Tíscar, Pedro Antonio, and Juan Carlos Linares. 2014. "Large-Scale Regeneration Patterns of Pinus nigra Subsp. salzmannii: Poor Evidence of Increasing Facilitation Across a Drought Gradient" Forests 5, no. 1: 1-20. https://0-doi-org.brum.beds.ac.uk/10.3390/f5010001

Article Metrics

Back to TopTop