Research Article
Print
Research Article
Analysis of the functional traits of Quercus cerris L. seedlings in the Molise region (southern Italy)
expand article infoLuca Quaranta, Piera Di Marzio, Romeo Di Pietro§, Fabrizio Ferretti|, Umberto Di Salvatore, Paola Fortini
‡ University of Molise, Pesche (IS), Italy
§ Sapienza University of Rome, Rome, Italy
| CREA Research Centre for Forestry and Wood, Arezzo, Italy
¶ CREA Research Centre for Politics and Bioeconomy, Pescara, Italy
Open Access

Abstract

This study deals with the analysis of seedling fitness in three Quercus cerris wood stands, namely Selva di Castiglione (SC), Bosco della Ficora (BF) and Bosco di San Leo (BSL), developed in different lithological and physiographic conditions and subjected to different forestry practices. A phytosociological study was carried out for analysing the coenological features of the forest stands and to classify these latter from a syntaxonomic point of view. The Pignatti-Ellenberg index (PEi) was calculated on the matrix composed of the phytosociological relevés in order to highlight possible ecological differences or gradients among stands. The phenotypic parameters considered were the stem and root length and the leaf area, whereas the plant functional traits (PFTs) were specific leaf area (SLA), leaf dry matter content (LDMC), leaf thickness (Lth) and chlorophyll content (CHL). The results showed that seedlings coming from different sampling sites exhibited similar values in all the phenotypic parameters. Instead, statistically significant differences were observed in the PFTs. The results suggested that the different adaptation strategies implemented by the seedlings are to be related to the physical environment of the sampling sites and to the different forest structures. The Selva di Castiglione forest stand (SC) exhibited better growth conditions for seedlings testified by higher values of SLA and CHL and lower values of LDMC and Lth. These were interpreted as greater investment in carbon production aimed at rapid development and renewal of the seedling rather than carbon storage aimed at ensuring leaf longevity.

Keywords

functional trait, high forest, oak, phenotypic response, phytosociology, southern Apennines

Introduction

Quercus cerris L. (Turkey oak, Fagaceae) acts as a dominant species in a large part of the oak forests of central and southern Italy (Blasi et al. 2004; Di Pietro et al. 2020b; Terzi et al. 2021). In the Molise region, the Q. cerris forests are widely distributed within the hilly, submontane and lower montane belt, ranging in altitude between 300 and 1200 m a.s.l. (Blasi and Paura 1995; Blasi et al. 2004, 2005). In particular, a wide Q. cerris belt separates thermophilous mixed oak woods with Q. cerris, Q. pubescens Willd. and Q. frainetto Ten. of the hilly belt from the beech woods of the montane belt. The Q. cerris belt is linked to Flysch substrates or the occurrence of large karst basins on limestone bordered by mild slopes covered by leached rich in clay soils. However, the abundance of Q. cerris woods in the Molise region is also due to the application of forestry management (e.g., different forest systems and silvicultural systems) that have advantaged this species over other tree species (e.g., Q. petraea (Matt.) Liebl. and Q. frainetto). The high-forest system of the Turkey oak woods, which is typically adopted in many parts of Italy and especially in the Molise region, finds its historical reason in the production of railway sleepers and, only marginally, of firewood (Di Martino 1996). Owing to the importance of Turkey oak for the regional economy, and paying attention to the ongoing climate change, we carried out an in-situ study aimed at identifying possible phenotypic response and plant functional traits (PFTs) variations from Q. cerris seedlings coming from wood stands characterised by different environmental conditions and forestry treatment. PFT variations have proved to be useful for establishing ecological and evolutionary trends for plant species and communities. They have often been used to interpret plant-environment relationships, to make previsions on possible consequences and to quantify a wide range of natural and/or human-driven processes in terms of species or community ecology (Pèrez-Harguindeguy et al. 2013; Adler et al. 2014; Pérez-Ramos et al. 2019). The measurements of the functionality patterns of seedlings of Q. cerris growing under different forest management, could provide useful information on the health state of the woods as well as address forest management policies in the next future. Surprisingly the available literature on such studies (i.e., observations in situ on Q. cerris seedlings) is relatively scarce in Italy. Most of the ecological studies performed on Q. cerris dealt with adult trees (e.g., Gratani and Foti 1998) or analysed, in greenhouses or under controlled conditions, the response of woody species seedlings to various types of induced stress (Otieno et al. 2005; Aref et al. 2013; Karavin et al. 2013; Chiatante et al. 2015; Abdala-Roberts et al. 2018; Torres-Ruis et al. 2019). Other studies have focused on the phenotypic plasticity of a species as a measure of its adaptability and fitness to environmental and climatic changes (Valladares et al. 2000, 2014; Pérez-Ramos et al. 2019). Field studies on the vitality of plants using PFTs, however, are more frequently applied to grassland communities (Bolzan et al. 2007; Catorci et al. 2013; Chelli et al. 2019; Baltieri et al. 2020) than to forest ones.

This study is part of an interdisciplinary and broader line of research that has been undertaken by the Laboratory of Systematic Botany and Floristics of the University of Molise in the last decade, which aims to shed light on the morphological, ecological, and bio-molecular features of Quercus species (Antonecchia et al. 2015; Fortini et al. 2015a, 2015b; Di Pietro et al. 2016; Conte et al. 2019; Di Pietro et al. 2020a, 2020c; 2021). Specifically, the present paper aims to evaluate the response of Quercus cerris seedlings to different ecological and forest management conditions, using PFTs and the analysis of morphological trait variations.

Materials and methods

Study area

The study was carried out in three Quercus cerris woods located in the submontane-lower montane belt of the Molise Region (Italy), namely Selva di Castiglione (SC), Bosco di San Leo (BSL), Bosco della Ficora (BF) (Fig. 1). According to Blasi et al. (2010) and Rivas-Martínez et al. (2011), the macro-bioclimate of the study area can be defined as temperate oceanic/sub-Mediterranean while based on the classifications of the Ecoregions of Italy (Blasi et al. 2018), the study area falls under code 1C3a1. The code identifies the territory according to the following characteristics: 1 = Temperate Division, C = Apennine Province, 3 = Southern Apennine Section, a = Campanian Apennine Subsection. The first wood stand (Selva di Castiglione–SC) is developed within the bottom of a plain on clayey substrates, while the other two stands (Bosco di San Leo–BSL and Bosco della Ficora–BF) are located within slopes with declivity on arenaceous-pelithic substrate. Physiographic information on the sites of collection is summarised in Table 1.

Figure 1. 

A) Location of the three sampled sites (Municipality of Carovilli) in the Province of Isernia in the Molise Region. BF = Bosco della Ficora, BSL = Bosco di San Leo, SC = Selva di Castiglione. B) Location of the Molise Region in Italy.

Table 1.

Characteristics of the sampled sites (Cantiani et al. 2010).

Coordinate (DMS) Stand development Age Altitude (m a.s.l.) Aspect Slope (%) Substrate
Bosco di San Leo (BSL) 41°40'33"N
14°15'18"E
Mature high forest 135 810 NW 15 arenaceous-pelithic
Bosco della Ficora (BF) 41°42'55"N
14°18'55"E
Young high forest 55 920 N-NW 15 arenaceous-pelithic
Selva di Castiglione (SC) 41°44'12"N
14°18'50"E
Adult high forest 90 1000 - 0 varicoloured clay

Silvicultural description

The sampled stands have already been the subject of experimentation on the Q. cerris forests by the Council for Research and Analysis of Agricultural Economics (CREA). From a forestry point of view, the three stands have different histories and are currently passing through different silvicultural stages. The BSL stand is an old forest (135 years) and is a result of a natural regeneration stage starting from a senescent high forest subjected to thinning interventions. The average diameter of Turkey oak individuals ranges between 30 and 35 cm. The BF stand is a young high forest (average age 55 years), which having never been subjected to silvicultural interventions, is currently characterised by an excessive density of tree individuals. The SC stand is an adult forest (average age ranging between 60 and 100 years). In 2006 and 2007, the following silvicultural treatments were applied to the three forest stands (Cantiani et al. 2010): BSL, shelterwood cutting (seed-tree); SC, high thinning (grade: heavy) and BF, selective thinning.

Data and methods

To obtain coenological information on the three Q. cerris stands a field-sampling campaign was carried out during the 2019 summer period using the phytosociological approach (Braun-Blanquet 1964). Species and syntaxa nomenclature in the phytosociological table followed Bartolucci et al. (2018) and Biondi et al. (2014), respectively. The following Ellenberg indicator values (L - light availability; T - temperature) and edaphic factors (U - soil moisture; N - soil nutrient) as reported in Pignatti (2005) were also considered to infer possible ecological differences between the three Q. cerris stands investigated. For the indicator values that were missing for some species in Pignatti (2005), reference was made to Sürmen et al. (2014). Subsequently, Pignatti-Ellenberg index was applied to the matrix of phytosociological relevés by calculating the weighted averaging using species abundances to evaluate the average response of each forest community to the abiotic factors associated with Ellenberg indicator values, following what is reported in Gristina and Marcenò (2008) from the unpublished data of Pignatti and Bona (2007).

The in-field collection of Quercus cerris seedlings followed mainly the protocol reported in Gottardini et al. (2016) according to the following steps: 1) identification of a sampling stand, with a random collection of thirty 3 to 4-year-old seedlings, starting from the centre of the site and maintaining a distance of 6–7 m between each seedling collected, choosing seedlings with fully expanded leaves free from damage from herbivores or pathogens (Garnier et al. 2001); 2) eradication of the seedlings by taking an entire clod of soil with a spade up to 30 cm deep and subsequent removing of the soil from the root through water; 3) measurements of the soil pH (classes of pH soil according to National Soil Survey Center 1998) precisely where the individual seedlings were eradicated; 4) measurement of the bulk density of the soil through a metal cylinder (classes of bulk density soil according to National Soil Survey Center 1998); 5) application of a damp cotton cover to each seedling and store it in dark in the cooler and 6) measurements of the chlorophyll leaf content by the at LEAF CHL PLUS instrument (FT Green LLC, Wilmington, US), whose function is comparable to the SPAD 502 (Minolta) (Zhu et al. 2012; Novichonok et al. 2016). The measurements were repeated three times in the central part of the lamina of three leaves directly in the field.

The first stage of field sampling was carried out during three consecutive days (24, 25 and 26 June 2019) under the same climatic conditions (no wind and no clouds). A total of 90 seedlings were collected and further analysed. The chlorophyll contents of the leaves of the seedlings were determined in August 2019 for 30 different seedlings in each stand due to difficulties in obtaining the instrumentation for measuring the chlorophyll content during the June 2019 field campaign. Laboratory measurements followed Pérez-Harguindeguy et al. (2013). Stem length, root length and leaf area of each seedling were measured on fresh materials. Three leaves for each seedling were fixed on a sheet and scanned using the Epson GT15000 scanner, and the leaf area (LA, cm2) was calculated using the ImageJ software (Rasband 1997–2018). Fresh leaf mass saturated (LFMsat.) was calculated (after seven hours from the collection) and dry leaf mass (LDM), was calculated after drying leaves for 72 h at 60 °C until a stable weight was achieved. The leaf measures were used to evaluate the following plant functional traits: specific leaf area (SLA), leaf dry matter content (LDMC), leaf thickness (Lth), leaf water content (LWC) and leaf mass per area (LMA) (Table 2).

It should be noted (Table 2) that some of the PFTs investigated express approximately the same concept providing redundant information (e.g., SLA/LMA or LDMC /LWC). Accordingly, we have decided to make a further selection considering three PFTs only, namely SLA, LDMC and Lth. A spreadsheet (Excel 2016) containing all the field and laboratory data was arranged and processed with XLSTAT 2020 3.1.3 software (Addinsoft 2005–2020). The data underpinning the analysis reported in this paper are deposited at “Mendeley Data” at https://doi.org/10.17632/6m2d4rxc7n.1. All the data were previously subjected to a normality check and then to a non-parametric comparison analysis, according to the Kruskal-Wallis method, with a paired comparison according to the Dunn method (with Bonferroni correction for significance, set at alpha = 0.05).

Table 2.

Description of the functional traits analysed (Pérez-Harguindeguy et al. 2013). §: the dimensionless value is converted into the total chlorophyll value (mg/cm2) using the conversion table https://www.atleaf.com/SPAD.

PFTs Unit Formula
Specific leaf area (SLA): is the area of a fresh leaf, divided by its oven-dry mass cm2/g LA/LDM
Leaf dry matter content (LDMC): is the oven-dry mass of a leaf, divided by its water-saturated fresh mass mg/g LDM/LFMsat.
Leaf thickness (Lth) - 1/SLA*LDMC
Leaf water content (LWC) - 1000-LDMC
Leaf mass per area (LMA) is the oven-dry mass of a leaf, divided by its leaf area g/cm2 1/SLA
Leaf chlorophyll content (CHL) mg/cm2 §

Results

Phytosociological survey

No published phytosociological relevés or previous syntaxonomic interpretation were available for the study area. All the three wood stands were assigned to a single association, Roso arvensis-Quercetum cerridis Ubaldi 2003 (syn: Aremonio agrimonioidis-Quercetum cerridis Blasi et al. 2005 ex Terzi et al. 2021). The phytosociological relevés (Table 3) showed some differences among the three sampled sites. Bosco della Ficora (BF) and Selva di Castiglione (SC) woods are characterised by a tree layer composed of Quercus cerris and Carpinus betulus. In central Europe, the abundance of Carpinus betulus is considered as favoured by clumsy cutting forest practices that lead the climax Quercus petraea and Q. robur forest to be fragmented and subsequently invaded by the hornbeam (Oberdorfer, 1992). In the case of our study area, the occurrence of Carpinus betulus in the dominant tree layer is a sign of moist and eutrophic soil conditions and a low degree of disturbance (Blasi et al. 2005). Indeed, the abundance of Carpinus betulus sharply decreases passing from northern Italy to southern Italy and this species normally does not participate in recolonisation stages after oak forest coppicing where it is substituted by Fraxinus ornus, Acer opalus subsp. obtusatum, Corylus avellana and Ostrya carpinifolia (the latter especially on limestone or calcareous marls). The BF and the SC woods show a mesophilous-microthermic aspect due to the presence of several Fagetalia sylvaticae species such as Asperula taurina, Euphorbia amygdaloides, Crataegus laevigata, Cardamine bulbifera, Neottia nidus-avis, Euphorbia dulcis, Viola reichenbachiana (Di Pietro et al. 2004). In one relevé from the BF wood, Abies alba also occurs. The SC stand exhibits abundance of species typical of the Quercus cerris undergrowth such as Fragaria vesca, Festuca heterophylla, Allium pendulinum and Luzula forsteri. Bosco di San Leo wood (BSL) acts as a more thermophilous Q. cerris wood, characterised by a thick secondary tree layer dominated by Carpinus orientalis, Fraxinus ornus and Acer campestre and a shrub layer dominated by Cornus sanguinea, Ligustrum vulgare, Hedera helix and Rubus ulmifolius (Taffetani et al. 2012). In this paper, we have classified the Roso arvensis-Quercetum cerridis in the alliance Erythronio-Carpinion betuli (Fagetalia sylvaticae) due to the occurrence of several Fagetalia species. However, this classification should be considered preliminary. A still unresolved debate is open on the classification of the Quercus cerris woods of the montane belt of central-southern Italy both at the class level (Carpino-Fagetea vs. Quercetea pubescentis) and at the alliance level (Erythronio-Carpinion vs. Physospermo-Quercion cerridis vs. Melittio-Quercion frainetto vs Geranio striati-Fagion). As this is not the paper in which to discuss this issue, we refer to the following bibliography: Biondi et al. 2002; Ubaldi 2003; Blasi et al. 2005; Di Pietro and Fascetti 2005; Rosati et al. 2005; Biondi et al. 2008, 2014; Mucina et al. 2017 and Di Pietro et al. 2020.

Syntaxonomic framework

Querco Roboris-Fagetea sylvaticae Br.-Bl. et Vlieger in Vlieger 1937

Fagetalia sylvaticae Pawlowski in Pawlowski, Sokolowski et Wallish 1928

Erythronio-Carpinion betuli (Horvat, 1958) Marincek in Wallnofer, Mucina et Grass 1993

Roso arvensis-Quercetum cerridis Ubaldi 2003

Pignatti-Ellenberg's indices

Pignatti-Ellenberg's indices calculated for each forest stand were found to be the following: the SC stand: L= 5.179, T= 5.885, U= 5.408 and N= 5.600. For the BF stand, they were L= 5.932, T= 6.590, U= 4.644 and N= 4.962. For the BSL stand, they were L= 5.627, T= 6.742, U= 4.185 and N=3.341. The comparison among stands showed for the SC stand, higher values for soil moisture (U) and nutrients (N) compared to those of the BF and the BSL stands. In contrast, the BF and the BSL stands exhibited higher values for temperature and light (T and L) (Fig. 2).

Bulk density and pH of soil

The bulk density is an indicator of soil compaction. The three sampled stands exhibit similar mean values of bulk density (p-value is not significant), although the BSL and the BF stands showed a wide range of values, whereas the SC stand showed a narrow range of values together with the presence of some outliers (Fig. 3). According to the Natural Resources Conservation Service (USDA-NRCS 2021), the values we have found are considered “ideal conditions for plant growth”. Based on the classes of the pH soil described by the National Soil Survey Center (1998), the SC stand is classified as “extremely acid”, with an average pH value of 4.5, the BF stand is classified “very strongly acid”, with an average value of 4.7 and the BSL stand is classified as “strongly acid”, with an average value of 5.6 (Fig. 3).

Morphological measurements

Seedling root lengths were found to be slightly shorter in the SC stand compared to the other two stands, whereas stem lengths were found to be slightly shorter for the BSL stand. However, these differences in the mean values of the stem and root lengths turned out not to be statistically significant, with the p-value very close to the threshold (alpha=0.05, Fig. 4). Differences in the values of leaf area among the three stands investigated turned out to also be statistically insignificant (Fig. 4).

Plant functional traits

Specific leaf area (SLA) measurements showed the highest mean value in the SC stand (265.551), which differed significantly from the BF stand (215.078) and the BSL stand (205.224) (p=0.000 and p=0.0001 respectively), the values of the latter being very similar (p=0.508) (Fig. 5). Leaf dry matter content (LDMC) measurements showed that the mean value of the SC stand (339.652) differs significantly from those of the BF stand (407.723) and the BSL stand (416.036), with the SC displaying the lowest mean value (p=0.001) (Fig. 5). For Lth, the low values obtained for the SC stand (1.362), differed significantly from those obtained for the BF stand (1.983) and the BSL stand (2.301) (p=0.001) (Fig. 5).

Leaf chlorophyll content

The SC stand exhibited the highest average value for CHL (0.022 mg/cm2). This value differs significantly from those displayed by the BSL (0.016 mg/cm2) and BF stands (0.014 mg/cm2) (Fig. 6).

Table 3.

Roso arvensis-Quercetum cerridis Ubaldi 2003.

Date 6/6/2019 6/6/2019 6/6/2019 6/6/2019 7/22/2010 7/22/2010 6/25/2019
Forest Stand acronym SC SC BF BF BSL BSL BSL
Relevé number 1 2 3 4 5 6 7
Aspect - - NNW E N NNE NNE
Elevation (m a.s.l.) 980 980 920 890 810 810 810
Slope (%) - - 10 5 5 10 5
Soil surface aspect, rocky outcrops (%) - - - - - - -
Soil surface aspect, stoniness (%) - - - - 1 - -
Upper tree layer cover (%) 65 65 55 70 65 70 50
Upper tree layer height (m) 20 20 22 17 18 20 28
Lower tree layer cover (%) 80 50 40 30 30 20 70
Lower tree layer height (m) 2.5 2.5 3.4 5.6 4,5 6 4
Shrubs layer cover (%) - - - 10 50 40 90
Shrubs layer height (%) - - - 2 1.5 2 2
Herbaceous layer cover (%) 60 40 95 75 30 40 10
Area (m2) 200 200 200 200 200 200 300
Coverage (%) 100 80 90 90 100 100 100
Species per relevé 26 26 35 26 22 22 24 frq.
Roso arvensis-Quercetum cerridis Ubaldi 2003
Rosa arvensis 1 + 1 1 2 + 1 V
Ligustrum vulgare 1 2 1 2 2 2 2 V
Lonicera caprifolium 2 2 2 2 + + . V
Crataegus laevigata + . . . . . 1 II
Erythronio-Carpinion betuli (Horvat 1958) Marincek in Wallnofer, Mucina & Grass 1993
Carpinus betulus 4 3 2 2 . + . IV
Primula vulgaris + 1 + . 1 + . IV
Cardamine bulbifera 2 2 + + . . . III
Pulmonaria apennina + . + . . + . III
Viola odorata . . . . + . . I
Carex sylvatica . . . . . . 1 I
Geranio versicoloris-Fagion Gentile 1970
Cyclamen hederifolium + + . + . . . III
Allium pendulinum 2 1 + . . . . III
Anemone apennina + + + + . . . III
Ranunculus lanuginosus . . + . 1 . . II
Geranium versicolor . . + . . . . I
Lathyrus venetus . . . + . . . I
Fagetalia sylvaticae Pawłowski in Pawłowski, Sokołowski & Wallisch 1928
Euphorbia amygdaloides . . 2 . . . 1 II
Euphorbia dulcis . . + + . . . II
Fraxinus excelsior +
Asperula taurina . . 1 . . . . I
Neottia nidus-avis . . + . . . . I
Mycelis muralis . + . . . . . I
Galium odoratum . . . . . + . I
Abies alba . . . 1 . . . I
Quercetalia pubescenti-petraeae Klika 1933
Brachypodium rupestre 1 1 2 1 . . + IV
Sorbus torminalis . + . . . 1 + III
Cornus mas . . 1 + . . 1 III
Fraxinus ornus . . . + . 1 1 III
Carpinus orientalis . . . . 1 2 2 III
Veronica chamaedrys + + . . . . . II
Aegonychon purpurocaeruleum . . 2 . . . . I
Helleborus foetidus . . . . . + . I
Scutellaria columnae . . . . . . 1 I
Querco-Fagetea Br.-Bl. & Vlieger in Vlieger 1937
Quercus cerris 4 3 3 4 4 4 4 V
Hedera helix 1 1 1 1 2 3 1 V
Rubus hirtus + 1 2 3 + + . V
Daphne laureola 1 + + . + + . IV
Acer campestre . . 2 2 2 1 3 IV
Ajuga reptans + 1 . . + . 1 III
Festuca heterophylla 2 2 1 + . . . III
Crataegus monogyna . . 1 2 + . + III
Fragaria vesca 2 2 + . . . . III
Cruciata glabra + + + + . . . III
Luzula forsteri 1 1 . . + . . III
Dioscorea communis . + + . . . . II
Digitalis micrantha + . . . . . . I
Lilium bulbiferum . 1 . . . . . I
Melica uniflora . . . + . . . I
Bromopsis ramosa . . . . + . . I
Campanula trachelium . . . . . . 1 I
Epipactis helleborine . . . . . . + I
Rhamno-Prunetea Rivas Goday & Borja ex Tuxen 1962
Euonymus europaeus . . 1 1 + 1 . III
Cornus sanguinea . . . . 3 2 1 III
Pyrus communis . . + + 2 . . III
Clematis vitalba . . 1 . . . 1 II
Prunus spinosa . . . . 1 + . II
Rubus umifolius . . . . . 1 3 II
Ulmus minor . . . . . . 1 I
Other species
Agrimonia eupatoria . . . + + + . III
Carex remota . . . . 1 . 1 II
Hieracium gr. murorum + . . . . . . I
Malva thuringiaca . . + . . . . I
Campanula glomerata . . . . . . + I
Dactylis glomerata . + . . . . . I
Astragalus glyciphyllos . . + . . . . I
Carex sylvatica . . 1 . . . . I
Arum italicum . . . + . . . I
Orobanche hederae . . . . . + . I
Figure 2. 

Ecograms for the three forest stands based on the average values for Pignatti-Ellenberg indices: L – light; T – temperature; K – continentality; U – soil moisture; N – soil nutrients.

Figure 3. 

Bulk density and pH of soil in the three wood stands. Non-parametric comparison analysis (Kruskal-Wallis test) and pair comparison according to the Dunn method with Bonferroni correction for significance. Number of replicates for each stand: 30. Legend: red cross = mean; horizontal bar in the box = median; lower limit of the box = first quartile; upper limit of the box = third quartile; empty dots and stars outer the whiskers' bounds = outliers; filled diamonds = minimum and maximum values. Empty dots: values found in the [Q1 – 3 (Q3 - Q1); Q1 - 1.5 (Q3 - Q1)] interval or in the [Q3 + 1.5 (Q3 - Q1); Q3 + 3 (Q3 - Q1)] interval. See Appendix 1, Table A1 for the summary statistics.

Figure 4. 

Stem length, root length and leaf area values in the three wood stands. Non-parametric comparison analysis (Kruskal-Wallis test) and pair comparison according to the Dunn method with Bonferroni correction for significance. Number of replicates for each stand: 30. Legend: red cross = mean; horizontal bar in the box = median; lower limit of the box = first quartile; upper limit of the box = third quartile; empty dots and stars outer the whiskers' bounds = outliers; filled diamonds = minimum and maximum values. Asterisks: values found outside the [Q1 - 3 (Q3- Q1); Q3 + 3 (Q3 - Q1)] interval. Empty dots: values found in the [Q1 – 3 (Q3 - Q1); Q1 - 1.5 (Q3 - Q1)] interval or in the [Q3 + 1.5 (Q3 - Q1); Q3 + 3 (Q3 - Q1)] interval. See Appendix 1, Table A2 for the summary statistics.

Figure 5. 

SLA, LDMC and Lth values in the three wood stands. Non-parametric comparison analysis (Kruskal-Wallis test) and pair comparison according to the Dunn method with Bonferroni correction for significance. Number of replicates for each stand: 30. Legend: red cross = mean; horizontal bar in the box = median; lower limit of the box = first quartile; upper limit of the box = third quartile; empty dots and stars outer the whiskers' bounds = outliers; filled diamonds = minimum and maximum values. Empty dots: values found in the [Q1 – 3 (Q3 - Q1); Q1 - 1.5 (Q3 - Q1)] interval or in the [Q3 + 1.5 (Q3 - Q1); Q3 + 3 (Q3 - Q1)] interval. See Appendix 1, Table A3 for the summary statistics.

Discussion

All three wood stands investigated were assigned to the same association, which is Roso arvensis-Quercetum cerridis Ubaldi 2003. Nevertheless, owing to their different silvicultural treatment and ecological features, the three forest stands exhibited clear floristic-coenological differences. The BSL stand is dominated by Q. cerris in the upper tree layer and Carpinus orientalis in the lower tree layer, and the coverage of the shrub layer was found to be very high (40–90%). The PEi values were high for T and L. The BF stand is characterised by a dominance of Quercus cerris, with Carpinus betulus having a subordinate role, and by a high coverage by the herb layer (75–90%). The PEi values were found to be high for T and L. The SC stand is characterised by the co-dominance of Quercus cerris and Carpinus betulus in the upper three layer and by an almost totally missing shrub layer. The PEi values were high for U and N and low for T and L.

The stem and root lengths and the leaf areas of the seedlings were found to be rather similar in the three wood stands. Even if the Q. cerris seedlings of the three wood stands investigated are phenotypically very close each other, they showed a different response to environmental conditions if the PFTs values are considered. In the ecological literature, high values of SLA are normally related to environments characterised by being rich in nutrients and/or by shaded conditions. The higher values of SLA showed by the saplings of the SC stand compared to those of the BF and the BSL testify to the propensity of the SC seedlings for an investment in carbon oriented to a rapid growth of the saplings rather than to store material for longevity. High SLA values are related to a lower occurrence of dense tissues, which, as a negative counterpart, makes these seedlings more palatable to grazing herbivores. This negative effect is balanced by a more rapid production of new leaves, linked to a higher availability in nutrients, water, as suggested by the PEi values (see Fig. 3) (Wilson et al. 1999; Shipley et al. 2002; Wright et al. 2004; Pontes et al. 2007).

In both the BSL and the BF stands, the values of LDMC are higher than those measured in the SC, showing, therefore, for the seedlings of the BSL and the BF a greater propensity to longevity (Dijkstra and Lambers 1989; Niemann et al. 1992; Garnier and Laurent 1994; Kazakou et al. 2006). The low SLA and the high LDMC values in the BSL and the BF highlight different strategies compared to the SC stand. We could indirectly address these differences to a lower presence of resources for leaf development. The low Lth values of the SC stand are related to seedlings with thinner leaves. This is another PFT that identifies leaves with a greater propensity to grow. The Lth values are normally negatively correlated with photosynthesis and growth, while being positively correlated with leaf longevity and construction costs (Vile et al. 2005).

Considering that PFT values, such as SLA, LDMC and Lth, are in many cases used to estimate the growth rate of seedlings (Poorter and Garnier 1996; Wright and Westoby 2001), it can be assumed that the SC seedlings grow faster than those of the BF and the BSL. In broad terms, high values in SLA and low values in LDMC correspond to a rapid production of biomass, whereas the low values in SLA and high values in LDMC indicate the conservation of nutrients through long-lived leaves (Reich 1998; Garnier et al. 2001). Assuming that the SC stand seedlings benefit from a greater number of resources that lead these seedlings to exhibit higher SLA values, the next step could be making hypothesis on which kinds of resources are important. In this case, the response derives directly from the PEi, which indicates, for the SC stand, a higher availability of nutrients (N) and soil moisture (U). The flat geomorphological profile of the SC stand and the clayey substrate may play a crucial role in leading to a high degree of soil water retention. While refraining from drawing up a ranking of PEi by importance, we underline the effect that soil moisture can have for the deciduous forest communities occurring in Mediterranean or sub-Mediterranean environments, these latter being usually characterised by more or less prolonged dry periods during late spring and summer.

The pH of the soil, on the other hand, would seem to provide values in contrast, as it goes from 4.5 (SC) to 4.7 (BF) and 5.6 (BSL). These values are more consistent with a typical oligotrophic Quercetalia robori-petraea wood than with eutrophic Quercus cerris woods of Erythronio-Carpinion. Instead, the pH values of the BSL site were perfectly in line with this latter alliance. It would be interesting to understand if these low pH values of the SC and the BF could at least partially explain why the montane-belt Quercus cerris forests often give rise to almost monophytic communities, or in any case, communities lacking other species of oaks (see for example Di Pietro and Tondi 2005). It is possible that in addition to the altitude, the low soil pH values lead Quercus cerris to be less subjected to the competition from other oak species that are common in this area (e.g. Q. pubescens and Q. frainetto), whereas acidophilic oaks (e.g., Q. petraea and Q. robur) that potentially could take advantage from such low pH values are missing due to their rareness in southern Italy.

Another resource that could explain the high SLA values observed for the SC stand could be the light radiation. The PEi values showed a significant difference between the low value (5.19) of the SC stand and the values of the BF and the BSL stands (5.932 and 5.627, respectively).

Observing the headline of the phytosociological table, it emerges that the SC wood stand, while showing a high coverage index at the upper tree layer is characterised by an almost total lack of shrub layer. On the contrary, the BSL wood stand shows the highest cover indices for the shrub layer. It is possible, therefore, that in the SC population, the high coverage of the tree layer was much more effective in reducing the amount of light available for the herb layer than the dense shrub layer occurring in the BSL stand. Accordingly, the high Lth values observed in the BSL stand could be due to the low coverage value of the upper tree layer. In this case, a shrub layer, composed of prostrate species, such as Hedera helix and Rubus ulmifolius occurring in the BSL would not seem to represent an effective impediment to light in reaching the herb layer.

As far as values of CHL are concerned, these were found to be higher in the SC stand than in the other two stands. A greater quantity of chlorophyll favours a greater absorption capacity of solar radiation. Based on published studies (e.g., Curran et al. 1990; Filella et al. 1995; Gitelson et al. 2003) higher chlorophyll content suggests also a higher photosynthetic potential. This is in accordance with other PFT values observed in this study, such as the lower Lth values and with the higher soil moisture, which would turn out to be higher in the SC stand than in the other two stands (see Fig. 2). The soil moisture is one of the factors that is normally considered as influencing the photosynthetic potential (Zhang et al. 2011; Aref et al. 2013). All these results seem to agree that the seedlings of the SC stand, due to the lower light irradiance and higher soil moisture, are more prone to a rapid development than the seedlings of the BF and the BSL stands. The latter (BLS), growing on slightly drier substrates and greater light intensity, exhibits seedlings developing thicker leaves to cope with lower water availability and higher irradiance (Wright et al. 2004; Poorter et al. 2009).

Figure 6. 

CHL values in the three wood stands. Non-parametric comparison analysis (Kruskal-Wallis test) and pair comparison according to the Dunn method with Bonferroni correction for significance. Number of replicates for each stand: 30. Legend: red cross = mean; horizontal bar in the box = median; lower limit of the box = first quartile; upper limit of the box = third quartile; empty dots and stars outer the whiskers' bounds = outliers; filled diamonds = minimum and maximum values. Empty dots: values found in the [Q1 – 3 (Q3 - Q1); Q1 - 1.5 (Q3 - Q1)] interval or in the [Q3 + 1.5 (Q3 - Q1); Q3 + 3 (Q3 - Q1)] interval. See Appendix 1, Table A4 for the summary statistics.

Conclusions

In this study, phenotypic parameters, and some plant functional traits of Quercus cerris seedlings were tested for the first time in natural forest communities to highlight possible adaptation strategies to environmental conditions. The results do not show a significant link between phenotypic plasticity and environmental parameters. It emerged (at least in our case), that the phenotypic expression is not related with the growth strategy of the seedlings. Accordingly, the phenotypic plasticity of the Quercus cerris seedlings should not be considered a target element to detect the fitness to ensure the growth of the tree. On the contrary, PFT values allowed for highlighting significant differences between the three wood stands considered and for hypothesising possible strategies adopted by Quercus cerris in the first years after germination, depending on the environmental characteristics of the sites. Apparently, the PFT values evidenced that the seedlings of Selva di Castiglione (SC) were following a different development strategy than those of the Bosco della Ficora (BF) and Bosco di San Leo (BSL). Thus, the SC seedlings strategy, expressed by higher values of SLA and CHL and lower values of LDMC and Lth leads to the production of thinner leaves and testifies a greater investment in carbon aimed to a rapid growth of the seedling. The presence of more favorable conditions in the SC stand for the development of the seedlings is supported by the Pignatti-Ellenberg index (PEi) values, which confirm for this site a greater availability of nutrients and soil humidity.

Beyond contributing to a relatively still poor national PFTs database with the first data of in situ PFTs from southern Italy regarding a species of enormous territorial value such as Quercus cerris, this study emphasises the importance of investigating on the relationships between young plant development and environmental conditions. The in situ studies on plant functional traits of seedlings of forest species, providing crucial information on their ecology, may assume a great importance to address forest management practices. This information is likely to prove useful in the near future, precisely in the light of recent Italian and European policies linked to the containment of global warming. The guidelines for the “urban and extra-urban forest plan” of the National Recovery and Resilience Plan (Mission 2; Component 4; Investment 3.1) recently published by the Italian Ecological transition Ministry require specific references to natural forest communities and native woody flora to propose models of new urban forests that are consistent with the environmental characteristics of the sites and with the natural potential vegetation types. Yet, the new EU forestry strategy, which plans to plant 3 billion more trees in the European territory by 2030, can become an opportunity to apply the principles of ecology and vegetation science to sustainable land management.

Bibliography

  • Abdala-Roberts L, Galmán A, Petry WK, Covelo F, de la Fuente M, Glauser G, Moreira X (2018) Interspecific variation in leaf functional and defensive traits in oak species and its underlying climatic drivers. PLOS ONE 13: e0202548. https://doi.org/10.1371/journal.pone.0202548
  • Adler PB, Salguero-Gómez R, Compagnoni A, Hsu JS, Ray-Mukherjee J, Mbeau-Ache C, Franco M (2014) Functional traits explain variation in plant life history strategies. Proceedings of the National Academy of Sciences of the United States of America 111: 740–745. https://doi.org/10.1073/pnas.1315179111
  • Antonecchia G, Fortini P, Lepais O, Gerber S, Legér P, Scippa GS, Viscosi V (2015) Genetic structure of a natural oak community in central Italy: evidence of gene flow between three sympatric white oak species (Quercus, Fagaceae). Annals of Forest Research 58: 205–216. https://doi.org/10.15287/afr.2015.415
  • Aref I, El Atta H, El Obeid M, Ahmed A, Khan P, Iqbal M (2013) Effect of water stress on relative water and chlorophyll contents of Juniperus procera Hochst. ex Endlicher in Saudi Arabia. Life Science Journal 10(4): 681–685.
  • Baltieri M, Fantinato E, Vecchio SD, Buffa G (2020) Intraspecific variability of leaf traits and functional strategy of Himantoglossum adriaticum H. Baumann. Plant Sociology 57: 105–112. https://doi.org/10.3897/pls2020572/03
  • Bartolucci F, Peruzzi L, Galasso G, Albano A, Alessandrini A, Ardenghi NMG et al. (2018) An updated checklist of the vascular flora native to Italy. Plant Biosystems 152(2): 79–303. https://doi.org/10.1080/11263504.2017.1419996
  • Biondi E, Casavecchia S, Biscotti N (2008) Forest biodiversity of the Gargano Peninsula and a critical revision of the syntaxonomy of the mesophilous woods of southern Italy. Fitosociologia 45: 93–127.
  • Biondi E, Casavecchia S, Pinzi M, Allegrezza M, Baldoni M (2002) The syntaxonomy of the mesophilous woods of the Central and Northern Apennines (Italy). Fitosociologia 39: 71–93.
  • Blasi C, Capotorti G, Copiz R, Guida D, Mollo B, Smiraglia D, Zavattero L (2018) Terrestrial Ecoregions of Italy. Map and Explanatory notes. Global Map s.r.l., Firenze, Italy.
  • Blasi C, Di R, Filesi L (2004) Syntaxonomical revision of Quercetalia pubescenti-petraeae in the Italian Peninsula. Fitosociologia 41: 87–164.
  • Blasi C, Fortini P, Grossi G, Presti G (2005) Faggete e cerrete mesofile nell'Alto Molise. Fitosociologia 42(2): 67–81.
  • Blasi C, Paura B (1995) Su alcune stazioni a Quercus frainetto Ten. in Campania ed in Molise: analisi fitosociologica e fitogeografica. Annali Botanica (Roma), Studi sul Territorio, 51, Suppl. 10 (2) (1993): 353–366.
  • Bolzan A, Ferroni L, Vecchiettini M, Speranza M (2007) Brachypodium pinnatum, Bromus erectus and Dactylis glomerata leaf functional traits in grassland communities with different environmental conditions. Journal of Mountain Agriculture in the Balkans 10: 8–15.
  • Cantiani P, Ferretti F, Pelleri F, Sansone D, Tagliente G (2010) Le fustaie di cerro del Molise. Analisi del trattamento del passato per le attuali scelte selvicolturali. Annals of Silvicultural Research 36: 25–36.
  • Catorci A, Tardella FM, Cutini M, Luchetti L, Paura B, Vitanzi A (2013) Reproductive traits variation in the herb layer of a submediterranean deciduous forest landscape. Plant Ecology 214: 737–749. https://doi.org/10.1007/s11258-013-0203-0
  • Chelli S, Ottaviani G, Simonetti E, Wellstein C, Canullo R, Carnicelli S, Andreetta A, Puletti N, Bartha S, Cervellini M, Campetella G (2019) Climate is the main driver of clonal and bud bank traits in Italian forest understories. Perspectives in Plant Ecology, Evolution and Systematics 40: 125478. https://doi.org/10.1016/j.ppees.2019.125478
  • Chiatante D, Tognetti R, Scippa GS, Congiu T, Baesso B, Terzaghi M, Montagnoli A (2015) Interspecific variation in functional traits of oak seedlings (Quercus ilex, Quercus trojana, Quercus virgiliana) grown under artificial drought and fire conditions. Journal of Plant Research 128: 595–611. https://doi.org/10.1007/s10265-015-0729-4
  • Conte AL, Di R, Iamonico D, Di Marzio P, Cillis G, Lucia D, Fortini P (2019) Oak decline in the Mediterranean basin: a study case from the southern Apennines (Italy). Plant Sociology 56(2): 69–80. https://doi.org/10.7338/pls2019562/051
  • Di Martino P (1996) Storia del Paesaggio Forestale del Molise (Sec. XIX-XX). Istituto Regionale per gli Studi storici del Molise “V. Cuoco”, Campobasso, 171 pp.
  • Di Pietro R, Fascetti S (2005) A contribution to the knowledge of Abies alba woodlands in the Campania and Basilicata regions (southern Italy). Fitosociologia 42: 71–95.
  • Di Pietro R, Izco J, Blasi C (2004) Contribute to the nomenclatural knowledge of the Fagus sylvatica beech-woodlands communities of southern Italy. Plant Biosystems 138 (1): 27–36-52. https://doi.org/10.1080/11263500410001684099
  • Di Pietro R, Di Marzio P, Antonecchia G, Conte AL, Fortini P (2020a) Preliminary characterization of the Quercus pubescens complex in southern Italy using molecular markers. Acta Botanica Croatica 79(1): 15–25. https://doi.org10.2478/botcro-2019-0026
  • Di Pietro R, Di Marzio P, Misano G, Silletti G, Wagensommer RP, Fortini P (2016) Evidence from multivariate morphometric study of the Quercus pubescens complex in southeast Italy. Botanica Serbica 40(1): 83–100. http://doi.org/10.5281/zenodo.48865
  • Di Pietro R, Fortini P, Ciaschetti G, Rosati L, Viciani D, Terzi M (2020b) A revision of the syntaxonomy of the Apennine-Balkan Quercus cerris and Q. frainetto forests and correct application of the name Melittio-Quercion frainetto. Plant Biosystems 154: 887–909. https://doi.org/10.1080/11263504.2019.1701127
  • Di Pietro R, Conte AL, Di Marzio P, Gianguzzi L, Spampinato G, Caldarella O, Fortini P (2020c) A multivariate morphometric analysis of diagnostic traits in southern Italy and Sicily pubescent oaks. Folia Geobotanica 55: 163–183. https://doi.org/10.1007/s12224-020-09378-0
  • Di Pietro R, Conte AL, Di Marzio P, Fortini P, Farris E, Gianguzzi L, Müller M, Rosati L, Spampinato G, Gailing O (2021) Does the genetic diversity among pubescent white oaks in southern Italy, Sicily and Sardinia islands support the current taxonomic classification? European Journal of Forest Research 140, 355–371. https://doi.org/10.1007/s10342-020-01334-z
  • Fortini P, Antonecchia G, Di Marzio P, Maiuro L, Viscosi V (2015a) Role of micromorphological leaf traits and molecular data in taxonomy of three sympatric white oak species and their hybrids (Quercus L.). Plant Biosystems 149: 546–558. https://doi.org/10.1080/11263504.2013.868374
  • Fortini P, Di Marzio P, Di R (2015b) Differentiation and hybridization of Quercus frainetto, Q. petraea, and Q. pubescens (Fagaceae): insights from macro-morphological leaf traits and molecular data. Plant Systematics and Evolution 301(1): 375–385. https://doi.org/10.1007/s00606-014-1080-2
  • Gitelson AA, Gritz † Y, Merzlyak MN (2003) Relationships between leaf chlorophyll content and spectral reflectance and algorithms for non-destructive chlorophyll assessment in higher plant leaves. Journal of Plant Physiology 160: 271–282. https://doi.org/10.1078/0176-1617-00887
  • Gottardini E, Cristofolini F, Cristofori A, Pollastrini M, Ferretti M (2016) Misura della clorofilla a, contenuto di clorofilla e tratti fogliari: campionamento, raccolta e misurazioni. Guida per studi in campo. Documento del progetto LIFE Future For CoppiceS, Azione B.2, pp 34. https://doi.org/10.13140/RG.2.2.19909.14562
  • Gratani L, Foti I (1998) Estimating forest structure and shade tolerance of the species in a mixed deciduous broad-leaved forest in Abruzzo, Italy. Annales Botanici Fennici 35: 75–83.
  • Gristina AS, Marcenò C. (2008) Gli indici di bioindicazione di Pignatti-Ellenberg nello studio floristico-vegetazionale del promontorio di Capo Zafferano (Sicilia nord-occidentale). Naturalista siciliano, 32(1-2): 61–96.
  • Karavin N (2013) Effects of leaf and plant age on specific leaf area in deciduous tree species Quercus cerris L. var. cerris. Bangladesh Journal of Botany 42: 301–306. https://doi.org/10.3329/bjb.v42i2.18034
  • Kazakou E, Vile D, Shipley B, Gallet C, Garnier E (2006) Co-variations in litter decomposition, leaf traits and plant growth in species from a Mediterranean old-field succession. Functional Ecology 20: 21–30. https://doi.org/10.1111/j.1365-2435.2006.01080.x
  • Mucina L, Bültmann H, Dierßen K, Theurillat JP, Raus T, Čarni A et al. (2016) Vegetation of Europe: hierarchical floristic classification system of vascular plant, bryophyte, lichen, and algal communities. Applied Vegetation Science 19: 3–264. https://doi.org/10.1111/avsc.12257
  • Niemann GJ, Pureveen JBM, Eijkel GB, Poorter H, Boon JJ (1992) Differences in relative growth rate in 11 grasses correlate with differences in chemical composition as determined by pyrolysis mass spectrometry. Oecologia 89: 567–573. https://doi.org/10.1007/BF00317165
  • Novichonok EV, Novichonok AO, Kurbatova JA, Markovskaya EF (2016) Use of the atLEAF+ chlorophyll meter for a non destructive estimate of chlorophyll content. Photosynthetica 54: 130–137. https://doi.org/10.1007/s11099-015-0172-8
  • Oberdofer E (1992) Süddentsche Pflanzengesellschaften. Teil IV: 157–159. G. Fisher. Jena.
  • Otieno DO, Schmidt MWT, Adiku S, Tenhunen J (2005) Physiological and morphological responses to water stress in two Acacia species from contrasting habitats. Tree Physiology 25: 361–371. https://doi.org/10.1093/treephys/25.3.361
  • Pérez-Harguindeguy N, Díaz S, Garnier E, Lavorel S, Poorter H, Jaureguiberry P et al. (2013) New handbook for standardised measurement of plant functional traits worldwide. Australian Journal of Botany 61: 167–234. https://doi.org/10.1071/BT12225
  • Pérez-Ramos IM, Matías L, Gómez-Aparicio L, Godoy Ó (2019) Functional traits and phenotypic plasticity modulate species coexistence across contrasting climatic conditions. Nature Communications 10: 2555. https://doi.org/10.1038/s41467-019-10453-0
  • Pignatti S (2005) Valori di bioindicazione delle piante vascolari della flora d'Italia. Braun-Blanquetia 39: 3–97.
  • Pignatti S, Bona E (2007) Ellenberg 2000. Ecological bioindication engine. EVS, Rome (unpublished)
  • Poorter H, Garnier E (1996) Plant growth analysis: an evaluation of experimental design and computational methods. Journal of Experimental Botany 47: 1343–1351. https://doi.org/10.1093/jxb/47.9.1343
  • Reich PB, Walters MB, Tjoelker MG, Vanderklein D, Buschena C (1998) Photosynthesis and respiration rates depend on leaf and root morphology and nitrogen concentration in nine boreal tree species differing in relative growth rate. Functional Ecology 12: 395–405. https://doi.org/10.1046/j.1365-2435.1998.00209.x
  • Rosati L, Di R, Blasi C (2005) La vegetazione forestale della Regione Temperata del Flysch del Cilento (Italia meridionale). Fitosociologia 42(2): 33–65.
  • Sürmen B, Kutbay HG, Kiliç DD, Hüseyinova R, Kilinç M (2014) Ellenberg's indicator values for soil nitrogen concentration and pH in selected swamp forests in the Central Black Sea region of Turkey. Turkish Journal of Botany 38: 883–895. https://doi.org/10.3906/bot-1311-43
  • Taffetani F, Catorci A, Ciaschetti G, Cutini M, Di Martino L, Frattaroli AR, Paura B, Pirone G, Rismondo M, Zitti S (2012) The Quercus cerris woods of the alliance Carpinion orientalis Horvat 1958 in Italy. Plant Biosystems 146: 918–953. https://doi.org/10.1080/11263504.2012.682613
  • Terzi M, Ciaschetti G, Fortini P, Rosati L, Viciani D, Di R (2021) New syntaxon names for the Italian Quercus cerris woods. Mediterranean Botany 42: e75592. https://doi.org/10.5209/mbot.75592
  • Torres-Ruiz JM, Kremer A, Carins-Murphy MR, Brodribb T, Lamarque LJ, Truffaut L, Bonne F, Ducousso A, Delzon S (2019) Genetic differentiation in functional traits among European sessile oak populations. Tree physiology 39: 1736–1749. https://doi.org/10.1093/treephys/tpz090
  • Ubaldi D (2003) La vegetazione boschiva d'Italia: manuale di Fitosociologia forestale. Clueb, Bologna.
  • Valladares F, Martinez-Ferri E, Balaguer L, Pérez-Corona ME, Manrique E (2000) Low leaf-level response to light and nutrients in Mediterranean evergreen oaks: A conservative resource-use strategy? New Phytologist 148: 79–91. https://doi.org/10.1046/j.1469-8137.2000.00737.x
  • Valladares F, Matesanz S, Guilhaumon F, Araújo MB, Balaguer L, Benito-Garzón M, Cornwell W, Gianoli E, van Kleunen M, Naya DE, Nicotra AB, Poorter H, Zavala MA (2014) The effects of phenotypic plasticity and local adaptation on forecasts of species range shifts under climate change. Ecology Letters 17: 1351–1364. https://doi.org/10.1111/ele.12348
  • Vile D, Garnier É, Shipley B, Laurent G, Navas ML, Roumet C et al. (2005) Specific Leaf Area and Dry Matter Content Estimate Thickness in Laminar Leaves. Annals of Botany 96: 1129–1136. https://doi.org/10.1093/aob/mci264
  • Wright IJ, Reich PB, Westoby M, Ackerly DD, Baruch Z, Bongers F et al. (2004) The worldwide leaf economics spectrum. Nature 428: 821–827. https://doi.org/10.1038/nature02403
  • Wright IJ, Westoby M (2001) Understanding seedling growth relationships through specific leaf area and leaf nitrogen concentration: generalisations across growth forms and growth irradiance. Oecologia 127: 21–29. https://doi.org/10.1007/s004420000554
  • Zhang YJ, Xie ZK, Wang YJ, Su PX, An LP, Gao H (2011) Effect of water stress on leaf photosynthesis, chlorophyll content, and growth of oriental lily. Russian journal of plant physiology 58: 844–850. https://doi.org/10.1134/S1021443711050268
  • Zhu J, Tremblay N, Liang Y (2012) Comparing SPAD and atLEAF values for chlorophyll assessment in crop species. Canadian Journal of Soil Science 92: 645–648. https://doi.org/10.4141/cjss2011-100

Appendix - Summary statistics

Bulk density (g/cm3) pH
Stand Minimum Maximum Average Standard deviation Minimum Maximum Average Standard deviation
SC 0.627 1 399 1 113 0.181 3.93 4.94 4 469 0.267
BSL 0.820 1 425 1 123 0.180 4.86 6.40 5 599 0.360
BF 0.713 1 472 1 174 0.181 4.26 5.17 4 666 0.233
Stem length (cm) Root length (cm) Leaf area (cm2)
Stand Minimum Maximum Average Standard deviation Minimum Maximum Average Standard deviation Minimum Maximum Average Standard deviation
SC 7 900 24 000 16 697 4 257 5 600 25 100 11 650 3 918 5 114 14 092 10 859 2 448
BSL 9 100 32 100 16 780 4 825 4 200 23 900 13 397 5 102 6 124 22 033 12 799 4 272
BF 11 600 30 500 18 633 3 698 6 800 23 500 13 707 4 512 4 710 23 806 11 495 4 742
SLA (cm2/g) LDMC (mg/g) Lth
Stand Minimum Maximum Average Standard deviation Minimum Maximum Average Standard deviation Minimum Maximum Average Standard deviation
SC 146 135 359 631 265 551 51 847 283 037 416 391 339 652 37 281 0.806 2 854 1 362 0.431
BSL 117 310 334 237 205 224 55 258 240 648 532 341 416 036 74 523 0.720 4 306 2 301 0.939
BF 152 591 299 809 215 078 34 135 298 322 448 317 407 723 31 926 1 022 2 887 1 983 0.419
CHL (mg/cm2)
Stand Minimum Maximum Average Standard deviation
SC 0.017 0.029 0.022 0.003
BSL 0.005 0.027 0.016 0.006
BF 0.008 0.023 0.014 0.004