Next Article in Journal
Spatial and Temporal Variation in Fecundity of Acropora spp. in the Northern Great Barrier Reef
Next Article in Special Issue
Lichen Responses to Disturbance: Clues for Biomonitoring Land-use Effects on Riparian Andean Ecosystems
Previous Article in Journal
No Signs of Genetic Erosion in a 19th Century Genome of the Extinct Paradise Parrot (Psephotellus pulcherrimus)
Previous Article in Special Issue
Testing the Poleotolerance Lichen Response Trait as an Indicator of Anthropic Disturbance in an Urban Environment
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Diversity
Volume 11
Issue 4
10.3390/d11040059
diversity-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Factors Influencing Epiphytic Lichen Species Distribution in a Managed Mediterranean Pinus nigra Arnold Forest

by
Ece Sevgi
1,*,
Osman Yalçın Yılmaz
2,
Gülşah Çobanoğlu Özyiğitoğlu
3,
Hüseyin Barış Tecimen
4 and
Orhan Sevgi
4
1
Department of Pharmaceutical Botany, Faculty of Pharmacy, Bezmialem Vakif University, Istanbul 34093, Turkey
2
Department of Surveying and Cadastre, Faculty of Forestry, Istanbul University-Cerrahpaşa; Istanbul 34473, Turkey
3
Goztepe Campus, Department of Biology, Faculty of Arts and Sciences, Marmara University; Istanbul 34722, Turkey
4
Department of Soil Science and Ecology, Faculty of Forestry, Istanbul University-Cerrahpaşa; Istanbul 34473, Turkey
*
Author to whom correspondence should be addressed.
Diversity 2019, 11(4), 59; https://doi.org/10.3390/d11040059
Submission received: 20 February 2019 / Revised: 5 April 2019 / Accepted: 12 April 2019 / Published: 15 April 2019
(This article belongs to the Special Issue Lichen Diversity and Biomonitoring)
Browse Figures
Versions Notes

Abstract

:
Lichens have important ecological functions in black pine forests, such as nitrogen fixation and nutrient cycling. Understanding lichen diversity could provide a better understanding of black pine ecosystems. The aim of this study was to identify the factors affecting the composition of lichen communities and their specific diversity in Mediterranean black pine forests. Research was conducted in 48 sampling plots. For the analysis, presence–absence and frequency data of lichen species were used. For stand level analysis, four community composition tables were created. We used bioclimate, topography, stand, and parent rock as variables. A total of 33 epiphytic lichen species were identified in the black pine forests from 282 sampled trees. Indicator lichen species were determined according to geographic region and stand age classes. Hypocenomyce scalaris was found to be an indicator species for old forests. Frequency data were more useful for revealing lichen species composition than presence–absence data. Of the topographic variables, elevation was the most prominent and had the highest explanation ratio for the composition of lichen species with a coefficient of correlation (R2) value of 0.49. Significantly positive (p < 0.001) relationships were found between epiphytic lichen richness and tree crown height, tree height, and bark pH. Our results revealed that to retain the trees in the stands rich in lichen species diversity is recommended in the managed forests.

1. Introduction

Mediterranean ecosystems have extremely high species diversity that can be managed within a conservation area [1]. Black pine tree with its five subspecies, is distributed in European and Asian mountainous areas [2,3,4]. The greatest distribution of black pine was recorded in Anatolia (Turkey), accountings for 2.5 million ha of its total worldwide distribution of 3.5 million ha [5]. The most commonly distributed altitude is 300–1800 m above sea level (a.s.l.) in sub-Mediterranean and Mediterranean-montane regions [6]. One of the optimum growing sites for black pine in Turkey is the Alaçam mountains [7] where rarely encountered double-layered stands of these trees can be found [5]. The Alaçam mountains have been found to have substantially high plant biodiversity [8], which is attributed to their elevation. The Alaçam mountains have a more continental climate due to their geographical position compared to Mt. Olympus and Mt. Trapezitsa (Greece). Pinus nigra Arnold is a common tree species in Mediterranean ecosystems as well as in these mountains, and the trees are known for hosting a vast number of lichen species.
Lichen species play an important role in nutrient cycling in forest ecosystems, providing food for wildlife and a nest for microarthropods and insects [9,10,11,12,13]. Lichen physiology is affected by exposure to solar radiation, environmental humidity, and temperature [14]. Therefore, lichens can be used as indicators for monitoring these factors related to forest management regimes, fragmentation, land use, and climate change impacts [15,16,17,18,19,20,21]. The sensitivity and early warning ability of lichen species [22] related to air pollution has led to their use as an indicator species for air quality monitoring programs [23,24,25,26,27,28].
The distribution of lichen species and the factors affecting their distribution must be determined to understand the functions of lichens. Lichen diversity is influenced by both the host tree species (tree height, age, diameter, branch density, pH of tree bark, bark structure, and water holding capacity of bark), by the position of trees in the stand (canopy closure, species mixture, number of trees per ha), and the site quality (aspect, slope, and altitude) [29,30,31,32,33,34,35,36,37,38,39].
Topographical variables influence the water and nutrient budget of a given site [40]. Many topographical variables have been included in models to explain the plant species distribution and plant society [40,41,42,43,44]. To determine the effects on lichen species diversity, some limited topographical variables are used, like altitude, exposure, and slope. The northness, general curvature, catchment slope, and slope length have rarely been incorporated in studies seeking to understand the relationship between lichen species diversity and topographical variables [45].
Burgaz et al. [46] reported that altitude, slope, tree diameter, and tree height affect epiphytic lichen species distribution in Spain. Lichen species diversity changes due to tree species in Mediterranean forests [19,47,48]. Lichen species distribution has been studied on the plane [49], oak [15], and beech [50] species in the Mediterranean. Belinchón et al. [15] found that lichen species vary according to the spatial position of the semi-deciduous forest trees species, i.e., whether they are located at the edge or the center of the stand. Besides the location of the trees, the distance of trees to the streams [51] and management strategies [17] affect the distribution of lichen species in the Mediterranean forests.
Several studies examined lichen species and their distribution in black pine forests [52,53,54]. Pirintsos et al. [52] detected 23 lichen species at Mount Olympus and 63 species were recorded at Mount Trapezitsa, both in Greece [53]. Guvenc et al. [54], recorded 20 species on black pine trees on Uludağ Mountain in Turkey.
In the present study, the distribution of lichen species was investigated at the stand level, and factors affecting the distribution were categorized into groups; climate, topography, stand, and bedrock types. We investigated the relationship of the richness of the lichen species with tree size attributes. Within the current study, we intended to clarify whether (i) bioclimatic, topographic, stand, and parent rock variables influence the lichen species composition; (ii) there is a particular lichen species that acts as an indicator for a geographic region or a life stage of black pine stands; and (iii) black pine tree size and bark properties affect lichen species richness.

2. Materials and Methods

2.1. Study Area

The study area included the Ulus, Eğrigöz, and Alaçam mountains, which together form the Alaçam Mountain range. The Alaçam mountain range is located partly in the cities of Balıkesir and Kütahya. The study site is located 29°15′30–28°15′00 E and 39°38′00–39°07′30 N (Figure 1). Based on the meteorological data gathered from the vicinity of study area, remarkable severe summer droughts occur in the region. The general climate type in the study area is Mediterranean, with the lowest and highest annual precipitation with 458.2 and 860.3 mm in Tavşanlı and Simav towns, respectively [55]. The most abundantly parent rocks, the original rock from which the soil was formed, are granit, tuff-anglomera, mélange, quartzite and dacite.
The single-layer Anatolia black pine (Pinus nigra Arnold subsp. pallasiana (Lamb.) Holmboe) forests in the study area are scattered over an area of 91,744 ha [55]. Some parts of the black pine forests are degraded, and a greater portion of the forests are productive forests managed by the age class method [56]. Black pine forests regenerate through the large area shield method [57]. The silvicultural maintenance period varies between 5 and 10 years, and the removed stem bark remains in the stand.

2.2. Data

2.2.1. Species Data

Epiphytic lichen samples were collected from 95 long-term observation sample plots. All epiphytic lichens were recorded from August to November 2007 within sampled plots. At each plot, five single stems were randomly selected, and the lichen samples were collected from lower trunks (ranging from 0 to 2 m above the ground) of living black pine trees.
The collected epiphytic lichens were treated with chemical spot test and diagnosed with stereomicroscope [58,59,60]. The distribution of species into morphological structures was as follows: 7 foliose, 15 fruticose, 10 crustose, and 1 squamulose. To eliminate the spatial autocorrelation, the sampling plots within 1 km or less from each other were discarded randomly with spThin package [61] in R (R Foundation for Statistical Computing, Vienna, Australia). The autocorrelation threshold was set to 1 km in a study conducted in beech forests harboring a high lichen species diversity in Bayern, Germany [62]. Since the resolution distance of bioclim variables was 1 km, and also regarding the previous study, we applied the autocorrelation threshold value of 1 km. The remaining 48 sampling plots (Dursunbey, n = 26; Sındırgı, n = 16; and Simav, n = 6) were subjected to stand level analyses in our study (Figure 1). Therefore, the species found in the 37 discarded sampling plots, such as Alectoria sarmentosa, Cladonia fimbriata, Cladonia ramulosa, and Parmelia sulcate, were excluded. In the 48 plots, lichen data were collected from a total of 282 trees. The lichen diversity was described using two univariate variables (species presence–absence and species frequency) and one multivariate (species community composition) variable. The presence of epiphytic lichen species in five trees in each sampling plot was converted into the frequency of lichens. A species community composition table was created using both presence–absence (LPA) and frequency (LFQ) data. The most abundant species covered in a portion of 30% were determined using the ordiselect function in the goeveg package [63]. Those selected species included: Bryoria capillaris, Buellia griseovirens, Hypocenomyce scalaris, Hypogymnia farinacea, Hypogymnia physodes, Hypogymnia tubulosa, Ochrolechia turneri, Platismatia glauca, Pseudevernia furfuracea, and Usnea subfloridana. The dataset of the most abundant 30% of species selected from LPA was renamed GLPA and that from LFQ was renamed GLFQ. Thus, four community composition data tables were created and subjected to stand-based analyses. Lichen species richness, tree dimensions, and bark properties were used for tree level analyses in all sampling plots (n = 95). After the missing tree height and canopy thickness data had been discarded, the 373 remaining trees were used in analyses. For tree-level analyses, the number of epiphytic lichen species per tree was used as the response variable.

2.2.2. Environment Data

Environmental variables were grouped into four main topics: bioclimatic, stand structure, topography and parent rock data (Table 1 and Appendix A). The variables for environment were sourced from a 30 arc-second (0.00833°, under 1 km) resolution dataset of 19 bioclimatic variables from WorldClim 1.4 [64] and the ENVIREM [65] dataset, in which an additional 16 climatic and 2 geomorphological variables were included.
Besides the climatic data, we used the microscale topographic variables calculated from field analyses, which also identify the ecological processes. SAGA GIS terrain analysis functions were used, and a total of 40 variables (slope, aspect, and curvature) were obtained from ASTER DEM at a 30-m resolution [66]. GRASS GIS was used to calculate the northness and eastness of the data [67].
The parent rock types, the only categorical variable in this study, were gathered from 1:25,000 scaled geological map provided by Mining Detection and Seeking General Directory (MTA).
In this study, forest structure variables affecting the lichen species richness and composition in Mediterranean forests were used. The ages of trees were determined by counting the year rings on the pencils cored from four to seven trees per stand. The canopy closure values were predicted via visual estimation, ranging between 0.1 and 1.0. The crown height refers to the height from ground level to the lowest level of the crown, and dry branch thickness is the distance between the lowest and highest dry branches on the black pine stem. Both crown height and dry branch thickness were measured using a VERTEX III Laser (Langsele, Sweden) instrument. The tree diameter at breast height (Dbh) were measured using calipers in two perpendicular dimensions.
The pH of black pine bark was measured from a 1:10 diluted deionized (DI) water solution of ground bark [30,68,69]. Total nitrogen content of the barks was determined using the Kjheldahl method [70] with a Buchi Auto Kjheldahl Unit K-370 device.

2.3. Statistical Analyses

Permutational multivariate analysis of variance (PERMANOVA) tests based on the Bray–Curtis ecological distance similarity indices were used to test the statistical significance of the epiphytic lichen community composition variances between study regions, age class (young age, <60 years; mature age, 60–120 years; and old age, >120 years), and parent rock types (granite (g), tuff-aglomera (ta), mélange (m), quartzite (q), and dacite (d)) using the “Adonis” function in the vegan package of R [71].
Indicator species analysis [72,73] was used to explain the relationship between the epiphytic lichen species and study regions and the age class of stands using the multipatt-multi-level pattern analysis function of the indicspecies R package [74]. Indicator species for one region and one age class or combination of region and age classes were calculated. Based on the indicator species analysis, the indicator species for one cluster solely, or combinations of clusters, were calculated. Only significant species at the 0.05 level with IndVal values greater than 0.5 were assessed further.
Ordination methods were used to investigate the relationships between epiphytic lichen species composition and environmental variables. To compare the explanatory power of various environmental variable sets, the variations of the four Y response tables operated with the varpart function in vegan were classified [71]. Partial redundancy analysis (RDA) was used for community matrices, and a partial multiple regression analysis was performed for vector-independent variables.
For the final RDA model, individual explanatory variables were used (Table 1, Appendix A). The significant explanatory variables were selected using the ordistep function in vegan [71]. Before variation partitioning and statistical selection, multicollinearity between the explanatory variables was checked by pairwise correlations [38]. Strongly correlated variables (r > 0.75) were excluded from the selection. The effect of explanatory variables was tested using F-statistics via Monte Carlo simulation with 999 permutations. The accepted significance level was 0.05.
The 8 most frequently found lichen species were detected in 48 investigated sample plots [54]. The variables were determined as being the most frequently encountered in the lichen species composition LPA dataset. Subsequently, species response curves of the most abundant lichen species were drawn using generalized linear models (GLM) [63].
The Pearson pairwise correlation method was applied for the tree-level analyses and significance values are indicated [75]. A linear regression analysis was applied to determine the lichen species richness and single tree properties, like the diameter at breast height, the tree height, crown height, dry branch thickness, the bark pH, and the N content of bark. All statistical analyses were conducted under the framework of R, version 3.5.2 (R Foundation for Statistical Computing, Vienna, Australia, 20 December 2018)-Eggshell Igloo [76].

3. Results

3.1. Epiphytic Lichens in Black Pine Forests

In the context of this study, 33 lichen species were identified from 48 sampled plots. The species Pseudevernia furfuracea, Hypogymnia farinacea, and Hypogymnia tubulosa were the most abundant in the studied sampling plots (Appendix B). The lichen Bryoria capillaris, Platismatia glauca, Hypogymnia physodes, Ochrolechia turneri, Buellia griseovirens, Hypocenomyce scalaris, Usnea subfloridana, Evernia prunastri, Parmelia saxatilis, and Tuckermanopsis chlorophylla were detected as common in over 30% of the sampled plots (Appendix B).
The analysis of the variance of the lichen species composition, including three regions, five main parent rock types, and three age groups, was tested with a permutational multivariate analysis of variance using the Bray–Curtis distance matrix. The analysis showed that the community composition of the lichens was significantly different in terms of compositional similarity (Bray–Curtis, R2ADONIS = 0.15, p = 0.001) in different regions, but no significant difference was detected for either parent rocks (Bray–Curtis, R2ADONIS = 0.11, p = 0.142) nor for age groups (Bray–Curtis, R2ADONIS = 0.06, p = 0.204).
The indicator species analysis of the LPA data only identified two species as significant indicator species in Simav: Bryoria implexa and Letharia vulpina. Platismatia glauca was identified as the only indicator species for the Dursunbey and Sındırgı regions (Table 2). The LFQ data indicator species analysis identified Usnea subfloridana as the only indicator species for the Simav region (Table 2).
Hypocenomyce scalaris was determined as the indicator of old forests. No specific lichen species was assigned as a significant indicator for young or mature forests or for parent rock groups.

3.2. Stand-Level Lichen Species Composition

The climate, topographic, and stand structure variables explained the largest part of the species compositional variance (Figure 2). Frequency data sets explained the lichen species composition more efficiently compared to the presence–absence dataset (Figure 2). The best explanation variation partition score of 0.72 was produced by the GLFQ data; LFQ, LPA, and GLPA resulted in explanation levels of 0.56, 0.34, and 0.26, respectively. Within the GLFQ dataset, topographical variables had an explanation level of 0.49; bioclimatic, parent rock, and stand structure variables had explanation rates of 0.32, 0.30, and 0.15, respectively (Figure 2).
The corrected explanation power values for four datasets were: GLPA R2 = 0.31, GLFQ R2 = 0.25, LPA R2 = 0.25, and LFQ R2 = 0.23 (Table 3). Six topographical variables, five stand variables, and three bioclimatic variables were incorporated into the analyses (Table 3). The altitude and mean diameter at breast height were included in all analyses. Bio15 and mCH were included in three analyses, whereas the remaining variables were included in one or two analyses (Table 3).
After forward selection, 10 explanatory variables for LPA, 7 explanatory variables for GLPA, 7 explanatory variables for LFQ, and 6 explanatory variables for GLPQ were detected as significant in the final RDA model (Table 3). Among the topographical variables, only the elevation (Elev), catchment slope (CtchSlp), downslope curvature (DwnCurv), hill index (HI), tangential curvature (TnCurv), and texture (Txtr) were explanatory for the lichen species composition (Table 3). Stand variables, such as the mean diameter breast height (mDbh), mean crown height (mCH), total basal area of black pine in ha−1 (hecTBA), maximum height (maxH), and standard deviation of crown height (sdCH) revealed significant relationships with the lichen species composition (Table 3). Only three of the bioclimatic variables—precipitation seasonality (Bio15), index of the degree of water deficit below water need (AridxThorn), and count of the number of months with a mean temperature greater than 10 °C (MonthTemp10) (Table 3)—were found to have a significant relationship with the lichen species composition. Parent rock was also found to have a significant relationship with the lichen species composition (Table 3). As a result of the RDA between the LPA data set and environmental variables, a strong and positive relationship with elevation on the first axis of RDA and a negative correlation with AridxThorn, MonthTemp10, and maxH were indicated. H. scalaris was shown to be correlated with AridxThorn and MonthTemp10. T. chlorophylla and H. physodes showed negative correlations with maxH. The second axis of RDA was negatively correlated with Txtr and TnCurv. U. subfloriana and L. vulpina were negatively correlated with Txtr. B. griseovirens was negatively correlated with hecTBA.
According to our results for the RDA between the GLPA data set and environment variables, the first axis of RDA is positively correlated with elevation. The second axis showed a negative correlation with mCH. We observed that H. scalaris and H. physodes species were affected negatively by mDbh together with the climate (Figure 3).
For the RDA performed between LFQ data and environment variables, the first axis of the RDA showed a positive correlation with elevation and a negative correlation with HI (Figure 3).
The first RDA axis indicated a positive correlation with mCH as a result of the RDA performed between GLFQ data and environmental variables. The second axis showed a positive correlation with Bio15. H. physodes was correlated with HI (Figure 3).
A high possibility of the presence of Pseudoverniea furfuracea, Hypogymnia farinacea, and Hypogymnia tubulosa was found regarding their abundance distribution (Figure 4A–F). We estimated that Plastimatia glauca was not present above 1450 m, and Hypogymnia physodes was not present above 1100 m a.s.l. (Figure 4A). The higher the average stand diameter at breast height (mDbh) was, the higher the probability of the presence of Bryoria capillaris and H. physodes; for Platismatia glauca, Ochrolecha turneri, and Buella griseovirens, the probability decreased (Figure 4B). The possibility of the presence of B. capillaris decreased and increased for P. glauca, as long as the increased precipitation seasonality (coefficient of variation) value of Bio15 increased (Figure 4C). Typically, the presence of B. capillaris and B. griseovirens decreased with increasing AridxThorn (Figure 4D). The analyses predicted that B. capillaris and B. griseovirens may disappear above the AridxThorn values of 72 and 68, respectively (Figure 4D). B. griseovirens and O. turneri were present above an altitude of 1100 m a.s.l. The probability of the presence of B. capillaris and H. physodes increased with increasing mCH (Figure 4E). O. turneri, P. glauca, B. capillaris, and B. griseovirens were present with a higher probability with a higher CtchSlp value (Figure 4F).

3.3. Tree-Level Lichen Species Richness

A Pearson pairwise correlation analysis was conducted between epiphytic lichen species richness and individual tree properties, including the Dbh, dry branch thickness, crown height, total nitrogen, and bark pH (Figure 5). Correlations between lichen species richness and Dbh and total nitrogen of bark were not significant (Figure 5). We detected positive correlations between the epiphytic lichen richness and crown height and tree height; and negative correlations with bark pH were significant (p < 0.001) (Figure 5). The correlations between epiphytic lichen diversity and dry branch thickness were significant (p < 0.05) (Figure 5).
A linear regression analysis was conducted between lichen species richness and single tree properties, like Dbh, tree height, crown height, and bark pH (Table 4). The models were significant but the R2 values were low; the highest R2 was 0.0577. There was no significant relationship between lichen richness and dry branch thickness and the total nitrogen content of bark (Table 4).

4. Discussion

4.1. Epiphytic Lichen Species in Black Pine Forests

In this study, we investigated the effects of environmental factors on lichen species diversity and species composition. In the black pine forests of the Alaçam mountains, 33 epiphytic lichen species were identified. At Mt. Olympus, which shares the same latitude as the Alaçam mountains, 23 epiphytic lichen species were found [52], and 73 species were detected at Mt. Trapezitsa [53]. The higher number of lichen species found at Mt. Trapezitsa compared to the other two locations could be attributed to the mixture of fir and beech species with black pines in the forests. At Mt. Uludağ, which is one of the closest rangelands to the Alaçam mountains, 20 lichen species were found [54]. The low number of lichen species found in the Alaçam mountains may be due to factors related to the site conditions and the tree morphology of black pine forests. The studies conducted in mixed tree species in pine forests support this hypothesis [38,62,77]. Some tree species host a greater lichen richness than black pines. For example, in a study in the Bolu mountains (Turkey), 72 lichen species were found on only five trees [78]. In a study conducted in Spain, higher lichen species diversity was recorded [48]. The lichen species diversity is reinforced by the exposure of tree stems to sunlight [79,80]. The findings from previous studies positively support our results in the black pine forests with a mono-species composition.
An obvious summer drought occurred in the Alaçam mountains due to their long distance from the oceanic humidity-conveying winds. Black pine forests are generally found in cold and dry sites [81,82]. The presence of typical crustose species in black pine forests is an indicator of these climatic conditions. In the current study, the percentage of crustose species was 27% of the total detected lichen species, which is relatively high in a particularly protected black pine forest. The branch structure affects the humidity level of the stem. Due to the sympodial structure of oak trees, the rainwater leaks through the stem, creating more favorable conditions for lichens. For pine trees, rainwater directly reaches the ground without flowing along tree branches, which in turn, produces dry substrate conditions for lichen species [77].
Hypogymnia physodes and Pseudevernia furfuracea were detected as the most widespread species at Mt. Olympus (Greece) [52]. The same species showed a wide distribution in the northwestern black pine forests at Mt. Uludağ [54]. Our research revealed that the most commonly distributed species in the Alaçam mountains were Pseudevernia furfuracea, Hypogymnia farinacea, Hypogymnia tubulosa, Bryoria capillaris, Platismatia glauca, Hypogymnia physodes, Ochrolechia turneri, Buellia griseovirens, Hypocenomyce scalaris, Usnea subfloridana, Evernia prunastri, Parmelia saxatilis, and Tuckermanopsis chlorophylla. P. furfuracea is a common species in black pine forests, which we also detected in the Alaçam black pine forests. The genus Hypogymnia was represented by three species in the Alaçam black pine forests, which has not been reported in previous studies.
Hypocenomyce scalaris is a species living in acidic tree bark [58]. H. scalaris has been found in the oldest stand life stage (>120 years). H. scalaris has also been found on higher than 75 cm height cedar trees in Mediterranean forests [83]. Sevgi et al. conducted a study in the same region and found that black pine bark pH is inversely correlated with age [84]. Marmor and Randlane concluded that H. scalaris occurs more frequently on tree bark with an air-pollution-induced low pH [85]. In the current study, the presence of H. scalaris in old stands revealed that H. scalaris is an indicator of this tree stand life stage.
Calicioid lichens are an indicator species of the age classes of trees. They are more common in old trees than young trees and are also known as indicators of ecological continuity [34,36]. The calicioid species were assessed to indicate species diversity in pine forests [86]. In the current study, Calicium glaucellum was found in 73-, 75-, and 82-year-old stands, and Chaenotheca chrysocephala was found in 78- and 82-year-old mature age class stands. The average Dbhs of the stands where C. glaucellum was found were 20, 22, and 31 cm, respectively; whereas the average Dbhs of stands where C. chrysocephala was found were 31 and 32 cm, respectively. The presence of these calicioid species indicates that age class diversity was preserved in the stand, demonstrating that old trees are protected by forest management in the black pine forests and this protection has a positive effect on the composition of epiphytic lichen species.

4.2. Stand-Level Lichen Species Composition

In studies conducted in the same region, annual precipitation was shown to affect the lichen microbiota [29,33,87]. Bioclimatic variables have been used in models to explain the lichen species distribution [88]. In the Mediterranean region, precipitation and temperature are prominent variables that affect the lichen distribution [89]. In our study, bioclimatic variables were responsible for 0.32 of the lichen species distribution variation (Figure 2D). The results of our research revealed that Bio15, AridxThorn, and MonthTemp10 affected the epiphytic lichen species composition. Due to the summer drought in the Alaçam mountains, the humidity regime of the forest was affected by precipitation seasonality (coefficient of variation) and the index of the degree of water deficit below the water need. We inferred that the factors affecting the humidity regime of the forest also affected the epiphytic lichen species composition.
Epiphytic lichen species diversity is affected by microclimatic conditions (air humidity, temperature, and light), and structural factors (canopy closure, vertical structure of the canopy, and shrub layer), both of which affect the microclimate in the stand [80]. The topographic variables, such as Elev, CtchSlp, DwnCurv, HI, TnCurv, and Txtr, were the most important factors that explained the composition of epiphytic lichen species on the stand scale. These variables directly affected the lichen species composition due to their influence on the humidity of the site. Two important macroclimatic factors, temperature and humidity, were clearly identified among the factors influencing the composition of the epiphytic lichen community [90]. Increased light penetration increases the temperature in the stand, altering the epiphytic lichen composition due to the lower humidity [91]. However, Dymytrova et al. reported that the most important factors affecting the diversity of epiphytic lichens in beech forests are aspect and slope [92]. Ardelean et al. noted the importance of northness [45]. In our study, slope, eastness, and northness did not have any significant effect on the lichen species composition.
In our study, altitude (R = 0.26) had the strongest correlation with lichen species composition. Altitude has been the focus of many studies on the effects on lichen species composition. Dymytrova et al. stated that in beech forests (Ukraine), altitude has one of the greatest effects on lichen species density, and at higher altitudes, lichen species density increased [92]. Aragon et al. [17] and Nascimbene et al. [93] reported the same result in oak forests in Spain and spruce alpine forests in Italy, respectively. Our results are in accordance with previous studies that described the effect of altitude on lichen species composition.
The stand variables that affect the climate of the site explain the diversity of epiphytic lichen species. In our study, mDbh and mCH explained more of the lichen species composition, whereas maxH, hecTBA, and sdCH did not. Dymytrova et al. stated that the mean Dbh is a major factor in the determination of lichen species diversity and composition [48,92]. Those stand variables are positively correlated with the size of the crown. The crown structure varies by tree species and significantly affects species composition by supplying various conditions for epiphytic lichens [36,77]. This effect is related to light availability and direct rain protection [36]. As the tree crown radius increases, the shadow effect of trees increases, which increases the epiphytic lichen species composition. Our findings are in accordance with studies that indicated that the width of the crown affects the climate in the stand [36,94]. We found that an increased crown radius (mCH) supports the epiphytic lichen species composition. Our results revealed that the crown size provided shelter from light and direct rain. Thus, the tree crown size positively affected the composition of the lichen species. The shadow of the high trees in the stand diminishes the light penetration and increases the humidity in the stands [95]. The maximum height of the stand had an explanation capacity for the lichen species composition of 0.18.
The species response curve analyses demonstrated that the presence of some lichen species was directly affected by altitude (Figure 4A). The analysis showed that Buellia griseovirens and Ochrolechia turneri may be present above 1100 m a.s.l. in black pine forests. Giordani and Incerti stated that B. griseovirens has a crucial relationship with a cold humid climate in a study focusing on the effects of climate on lichen species at beech forests above 1000 m a.s.l. (Italy) [89]. Platismatia glauca is detected above 1100 m and its presence increases in parallel to precipitation seasonality. In spruce forests, it has been detected above 1400 m a.s.l. [93]. In the current study area, H. physodes was found below 1100 m. The same species has been found above 1400 m in beech and spruce forests [89,93]. The Mediterranean-type climate may be responsible for that difference. The location seems to have different effects on the presence of the same lichen species at varying altitudes.
Although Ochrolechia turneri has been detected in cedar trees with more than 75 cm height in Mediterranean forests [83], in the current study area, the possibility of O. turneri presence decreased in trees with a mDbh greater than 30 cm (Figure 4B). The bark properties of varying trees may differ at a given diameter. The diameter of any tree species may vary regarding to the tree species.

4.3. Tree-Level Lichen Richness

Properties of trees, such as bark pH, structure, and branch density affect the diversity of epiphytic lichens [29,37,38,39]. The lichen species richness was lower on black pine trees than on beech trees, whereas it was higher at breast height than at the base of the stem in beech trees in pine-beech mixed forest [50]. In the Mediterranean region, the number of lichen species changes due the diameter classes of cedar trees [83]. In our study, there was no significant relationship between lichen species diversity and Dbh (Figure 5), whereas in some studies, significant relationships were detected between lichen species diversity and diameter breast height [34,92,96,97,98,99,100] in forests other than black pine forests. Thus, the properties of black pine trees provide different conditions as a substrate for epiphytic lichen richness.
Due to their bark characteristics, coniferous forests generally form an acidic substrate environment for epiphytic lichen species [91] As tree diameter increases in black pine trees, the acid value decreases [79]. Hauck stated that the bark pH of various conifer tree species varies between 3.0 and 4.5 [91]. The bark pH of black pine in the current study area ranged between 3.2 and 4.8 [84]. The average bark pH of black pine in our study was 4.0, varying between 3.3 and 4.9 (Figure 5). The linear regression between the bark pH and lichen richness had an R2 value of 0.05 (p < 0.0001) (Table 4). The linear regression between the bark pH and the lichen richness had an R2 value of 0.05 (p < 0.0001) (Table 4). In the study by Sevgi et. al. (2016), they found a negative correlation between bark pH and tree age. They attributed this relation to the loss of bark leaves over a certain age of black pine trees [84]. Therefore, the decrease in the bark pH as the Dbh increased had a negative effect on the lichen species diversity.
The age and size of the tree affect the epiphytic lichen species diversity. There is greater epiphytic lichen species diversity on larger and older trees than in young trees [34,35,38,92]. There is a significant correlation between tree height and epiphytic lichen species diversity [36]. The results of our study support this relationship.
According to Ódor et al. [38], tree size does not affect epiphytic lichen species diversity. In that, this situation was reported due to the lack of old and large trees in the study area because of forest maintenance applications [38]. However, in our study, the presence of trees of different ages, as a result of forest management practices in stands, enriched the composition of epiphytic lichen species.

5. Conclusions

Both presence–absence and frequency data were used for lichen species composition analyses. We noticed that the frequency of presence data was more powerful for the prediction of lichen species composition than presence–absence data. Topographical variations, such as elevation, catchment slope, downslope curvature, hill index, tangential curvature, and texture, were more explanatory of lichen species composition compared to bioclimatic, stand, and parent rock variables at the stand level. The bioclimatic variables Bio15, MonthTemp10, and AridxThorn and the stand variables mDbh, mCH, hecTBA, maxH, and sdCH were shown to have significant effects on lichen species composition. The lichen Bryoria implexa was found to be an indicator species for the Simav region, and Hypocenomyce scalaris was found to be an indicator species for stands older than 120 years. There was a positive relationship between tree height and lichen species composition, a negative relationship with the pH of tree bark, and no relationship with age. We suggest that the effects of topographical indices on lichen species composition should be determined within various site conditions since they have the potential to explain lichen species composition.

Author Contributions

Conceptualization, O.Y.Y., E.S., and O.S.; methodology, O.Y.Y., E.S., and O.S.; software, O.Y.Y., investigation, E.S., O.S., O.Y.Y., G.Ç.Ö., and H.B.T.; writing—original draft, E.S., O.Y.Y., O.S., G.Ç.Ö., and H.B.T.; writing—review and editing, G.Ç.Ö., H.B.T., E.S., O.S., and O.Y.Y.

Funding

This work was supported by the Research Fund of The University of Istanbul with project number 1527.

Acknowledgments

We highly appreciate the efforts of the MDPI editing service for the language editing of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Used Variables in Varpart after Multicorelation

Variable GroupVariablesDescription
Bioclimatic (n = 6)Bio4 (Temperature seasonality (SD x 100)),
Bio15 (Precipitation seasonality, coefficient of variation),
AridxThorn (index of the degree of water deficit below water need),
MonthTemp10 (count of the number of months with mean temperature greater than 10 °C),
TopoWet (SAGA-GIS topographic wetness index),
Terrough (terrain roughness index),		[64,65]
Topographic (n = 20)Elev (Elevation),
East (Eastness),
Nrth (Northness),
MBI (Mass balance index),
PrI (Protection index),
MRRTF (Multiresolution index of the ridge top flatness),
slp (Slope),
GnCurv (General curvature),
TnCurv (Tangential curvature),
MnCurv (Minimal curvature),
ToCurv (Total curvature),
FLCurv (Flow line curvature),
Txtr (Texture),
Cnvx (Convexity),
TPI (Topographic position index),
DwnCurv (Downslope curvature),
ValDp (Valley depth),
HI (Hill index),
HSI (Hillslope index),
VTR (Vector terrain ruggedness),
ValDp1 (Valley depth 1),
CtchSlp (Catchment slope),
SlpLen (Slope length),		Derived from DEM with SAGA GIS terrain analyses
Stand structuremAge (mean value of increment cores collected from 4–7 trees in stand),
minAge (minimum stand age),
sdAge (Standart Deviation of stand age),
CaCl (Canopy closure degree of trees in the stand)
hecTBA (Total basal area of black pine ha−2),
mDbh (Mean diameter breast height of black pine trees in the stand),
sdDbh (Standard deviation of black pine breast height diameter values),
maxH (Maximum height in stand),
mCH (Mean crown height of black pine trees in the stand),
sdCH (Standard deviation of crown height of black pine trees in the stand),
mDBT (Mean dry branch thickness of black pine trees in the stand),		Calculated from study field data
TreeDbh (Diameter at breast height of tree),
Bark pH (Tree bark pH),
BarkN (Tree bark Nitrogen),
Tree height,
CH (Crown height of black pine tree),
DBT (Dry branch thickness of black pine tree),

Appendix B. Lichen Species According to Region and Type in Black Pine Forests

SpeciesRegionType
Dursunbey (n = 26)Simav (n = 6)Sindirgi (n = 16)
Bryoria capillaris (Ach.) Brodo & D.Hawksw.22611fruticose
Bryoria fuscescens var. fuscescens (Gyeln.) Brodo & D.Hawksw.101fruticose
Bryoria implexa (Hoffm.) Brodo & D.Hawksw.630fruticose
Buellia erubescens Arnold001crustose
Buellia griseovirens (Turner & Borrer ex Sm.) Almb.1368crustose
Calicium glaucellum Ach.300crustose
Chaenotheca chrysocephala (Turner ex Ach.) Th.Fr.110crustose
Cladonia chlorophaea (Flörke ex Sommerf.) Spreng.201crustose
Evernia divaricata (L.) Ach.301fruticose
Evernia prunastri (L.) Ach.816fruticose
Hypocenomyyce scalaris (Ach. ex Lilj.) M.Choisy1413squamulose
Hypogymnia farinacea Zopf26516foliose
Hypogymnia physodes (L.) Nyl.2047foliose
Hypogymnia tubulosa (Schaer.) Hav.24615foliose
Lecanora chlarotera Nyl.100crustose
Lecidella elaeochroma (Ach.) M. Choisy f. elaeochroma001crustose
Letharia vulpina (L.) Hue730fruticose
Ochrolechia turneri (Sm.) Hasselrot1926fruticose
Parmelia saxatilis (L.) Ach.1014foliose
Parmeliopsis ambigua (Wulfen) Nyl.221foliose
Pertusaria albescens (Huds.) M. Choisy & Werner var. albescens400crustose
Pertusaria hemisphaerica (Flörke) Erichsen110crustose
Platismatia glauca (L.) W.L. Culb. & C.F. Culb.23114foliose
Pseudevernia furfuracea var. furfuracea (L.) Zopf26616fruticose
Ramalina farinacea (L.) Ach.200fruticose
Trapeliopsis flexuosa (Fr.) Coppins & P.James100crustose
Tuckermanopsis chlorophylla (Willd.) Hale1203foliose
Usnea florida (L.) Weber ex F.H. Wigg.010fruticose
Usnea glabrescens (Nyl. ex Vain.) Vain.103fruticose
Usnea hirta (L.) Weber ex F.H. Wigg.100fruticose
Usnea scabrata Nyl.713fruticose
Usnea subfloridana Stirt.1043fruticose
Usnea subscabrosa Nyl. ex Motyka212fruticose

References

  1. Numa, C.; Valderrábano, M.; Alcázar, E.; Barrios, V.; Bourlon, S.; Bugalho, M.; Buse, J.; Vecchia, I.D.; Garnier, E.; Kabouya-Loucif, I.; et al. Biodiversity conservation and protected areas in the Mediterranean region. In State of Mediterranean Forests 2018; Food and Agriculture Organization of the United Nations: Rome, Italy; Marseille, France, 2018; pp. 147–162. ISBN 978-92-5-131047-2. [Google Scholar]
  2. Barbero, M.; Loisel, R.; Quezel, P.; Richardson, D.M.; Romane, F. Pines or the Mediterranean Basin. In Ecology and Biogeography of Pinus, 1st ed.; Richardson, D.M., Ed.; Cambridge University Press: Cambridge, UK, 1998; pp. 153–170. ISBN 0521551765. [Google Scholar]
  3. Bussotti, F. Pinus nigra Arnold. In Pines of Silvicultural Importance; CABI Publishing: Wallingford, UK, 2002; pp. 266–286. [Google Scholar]
  4. Akkemik, U.; Yılmaz, H.; Oral, D.; Kaya, A. Pinus. In Türkiye’nin Doğal Gymnospermleri (Açık Tohumlular), 1st ed.; Yaltırık, F., Akkemik, U., Eds.; Turkish Ministry of Enviroment and Forestry Press: Ankara, Turkey, 2011. [Google Scholar]
  5. Isajev, V.; Fady, B.; Semerci, H.; Andonovski, V. EUFORGEN Technical Guidelines for Genetic Conservation and Use for European Black Pine (Pinus nigra); International Plant Genetic Resources Institute: Rome, Italy, 2004; 6p. [Google Scholar]
  6. Mayer, H.; Aksoy, H. Türkiye Ormanları; Orman Bakanlığı Batı Karadeniz Ormancılık Araştırma Enstitüsü Müdürlüğü Yayınları: Bolu, Turkey, 1998; ISBN 975-7829-56-0.
  7. Kalıpsız, A. Türkiye’de Karaçam Meşcerelerinin Tabii Büyümesi ve Verim Kudreti Üzerine Araştırmalar; OGM Publications Sequence Nu: İstanbul, Turkey, 1963. [Google Scholar]
  8. Kavgacı, A.; Sevgi, O.; Tecimen, H.B.; Yılmaz, O.Y.; Serdar Carus, S.; Türker Dündar, T. Classification and Ordination of Pinus nigra Dominated Forests at Alaçam Mountains (Nw Anatolia-Turkey). Eurasian J. For. Sci. 2013, 1, 38–50. [Google Scholar]
  9. Gerson, U.; Seaward, M.R.D. Lichen–Invertebrate Associations. In Lichen Ecology; Seaward, M.R.D., Ed.; Academic Press: London, UK, 1977; pp. 69–120. ISBN 0-12-634350-0. [Google Scholar]
  10. Peck, J.E.; McCune, B. Remmant Trees and Canopy Lichen Communities in Western Oregon: A Retrospective Approach. Ecol. Appl. 1997, 7, 1181–1187. [Google Scholar] [CrossRef]
  11. Huneck, S. The Significance of Lichens and Their Metabolites. Naturwissenschaften 1999, 86, 559–570. [Google Scholar] [CrossRef]
  12. Sillett, S.C.; McCune, B.; Peck, J.E.; Rambo, T.R.; Ruchty, A. Dispersal Limitations of Epiphytic Lichens Result in Species Dependent on Old-growth Forests. Ecol. Appl. 2000, 10, 789–799. [Google Scholar] [CrossRef]
  13. Price, K.; Hochachka, G. Epiphytic Lichen Abundance: Effects of Stand Age and Composition in Cosatal British Columbia. Ecol. Appl. 2001, 11, 904–913. [Google Scholar] [CrossRef]
  14. Nash, T.H. Nitrogen, its metabolism and potential contribution to ecosystems. In Lichen Biology, 1st ed.; Nash, T.H., Ed.; Cambridge University Press: New York, NY, USA, 1996. [Google Scholar]
  15. Belinchón, R.; Martínez, I.; Escudero, A.; Aragón, G.; Valladares, F. Edge Effects on Epiphytic Communities in a Mediterranean Quercus pyrenaica Forest. J. Veg. Sci. 2007, 18, 81–90. [Google Scholar] [CrossRef]
  16. Nöske, N.M.; Hilt, N.; Werner, F.A.; Brehm, G.; Fiedler, K.; Sipman, H.J.M.; Gradstein, S.R. Disturbance effects on diversity of epiphytes and moths in a montane forest in Ecuador. Basic Appl. Ecol. 2008, 9, 4–12. [Google Scholar] [CrossRef]
  17. Aragón, G.; Martínez, I.; Izquierdo, P.; Relinchón, R.; Escudero, A. Effects of forest management on epiphytic lichen diversity in Mediterranean forests. Appl. Veg. Sci. 2010, 13, 183–194. [Google Scholar] [CrossRef] [Green Version]
  18. Pinho, P.; Correia, O.; Lecoq, M.; Munzi, S.; Vasconcelos, S.; Gonçalves, P.; Rebelo, R.; Antunes, C.; Silva, P.; Freitas, C.; et al. Evaluating green infrastructure in urban environments using a multi-taxa and functional diversity approach. Environ. Res. 2016, 147, 601–610. [Google Scholar] [CrossRef]
  19. Nascimbene, J.; Benesperi, R.; Bruniaiti, G.; Catalano, I.; Vedove, M.D.; Grillo, M.; Isocrono, D.; Matteucci, E.; Potenza, G.; Puntillo, D.; et al. Patterns and drivers of ß-diversity and similarity of Lobaria pulmonaria communities in Italian. J. Ecol. 2013, 101, 493–505. [Google Scholar] [CrossRef]
  20. Cardós, J.L.H.; Martínez, I.; Calvo, V.; Aragón, G. Epiphyte communities in Mediterranean fragmented forests: Importance of the fragment size and the surrounding matrix. Landsc. Ecol. 2016, 31, 1975–1995. [Google Scholar] [CrossRef]
  21. Pinho, P.; Dias, T.; Cruz, C.; Sim Tang, Y.; Sutton, M.A.; Martins-Loução, M.A.; Branquinho, C. Using lichen functional diversity to assess the effects of atmospheric ammonia in Mediterranean woodlands. J. Appl. Ecol. 2011, 48, 1107–1116. [Google Scholar] [CrossRef] [Green Version]
  22. Scott, M.G.; Hutchinson, T.C. Experiments and Observations on Epiphytic Lichens as Early Warning Sentinels of Forest Decline. In Biologic Markers of Air Pollution Stress and Damage in Forests; The National Academies Press: Washington, DC, USA, 1989; pp. 205–215. [Google Scholar]
  23. Jackson, L.; Ford, J.; Schwatzman, D. Collection and Chemical Analysis of Lichens for Biomonitoring. In Lichens as Bioindicators of Air Quality; USDA Forest Service General Technical Report RM-224; Huckaby, L.S., Ed.; USDA Forest Service: Fort Collins, CO, USA, 1993; pp. 96–115. [Google Scholar]
  24. Aamlid, D.; Torseth, K.; Venn, K.; Stuanes, A.O.; Solberg, S.; Hylen, G.; Christophersen, N.; Framstad, E. Changes of Forest Health in Norwegian Boreal Forest During 15 Years. For. Ecol. Manag. 2000, 127, 103–118. [Google Scholar] [CrossRef]
  25. McCune, B.; Rogers, P.; Ruchty, A.; Ryan, B. Lichen Communities for Forest Health Monitoring in Colorado, USA; A Report to the USDA Forest Service; Forest Health Monitoring National Office, Southern Research Station, Research Triangle Park: Washington, DC, USA, 1998; 29p.
  26. McCune, B. Lichen Communities as Indicators of Forest Health. Bryologist 2000, 103, 353–356. [Google Scholar] [CrossRef]
  27. Garty, J. Biomonitoring Atmospheric Heavy Metals with Lichens: Theory and Application. CRC Crit. Rev. Plant Sci. 2001, 20, 309–371. [Google Scholar] [CrossRef]
  28. Asta, J.; Erhardt, W.; Ferretti, M.; Fornasier, F.; Kirschbaum, U.; Nimis, P.L.; Purvis, O.W.; Pirintsos, S.; Scheidegger, C.; Haluwyn, C.V.; et al. Mapping Lichen Diversity as in Indicator of Environmental Quality. In Monitoring with Lichens-Monitoring Lichen; Nimis, P.L., Scheidegger, C., Wolseley, P.A., Eds.; NATO Science Series (Series IV: Earth and Environmental Sciences); Springer: Dordrecht, The Netherlands, 2002; Volume 7, pp. 273–279. [Google Scholar]
  29. Barkmann, J.J. Phytosociology and Ecology of Cryptogamic Epiphytes; Van Gorcum: Assen, The Netherlands, 1958. [Google Scholar]
  30. Kuusinen, M. Epiphyte Flora and Diversity on Basal Trunks of Six Old-Growth Forest Tree Species in Southern and Middle Boreal Finland. Lichenologist 1996, 28, 443–463. [Google Scholar] [CrossRef]
  31. Dietrch, M.; Scheidegger, C. Frequency, Diversity and Ecological Strategies of Epiphytic Lichens in the Swiss Central Plateau and Pre-alps. Lichenologist 1997, 29, 237–258. [Google Scholar] [CrossRef]
  32. Neitlich, P.N.; McCune, B. Hostpots of Epiphytic Lichen Diversity in Two Young Managed Forests. Conserv. Biol. 1997, 11, 172–182. [Google Scholar] [CrossRef]
  33. Giordani, P. Variables influencing the distribution of epiphytic lichens in heterogeneous areas: A case study for Liguria, NW Italy. J. Veg. Sci. 2006, 17, 195–206. [Google Scholar] [CrossRef]
  34. Fritz, Ö.; Niklasson, M.; Churski, M. Tree age is a factor for the conservation of epiphytic lichens and bryophytes in beech forests. Appl. Veg. Sci. 2008, 12, 93–106. [Google Scholar] [CrossRef]
  35. Lie, M.H.; Arup, U.; Grytnes, J.A.; Ohlson, M. The importance of host tree age, size and growth rate as determinants of epiphytic lichen diversity in boreal spruce forests. Biodivers. Conserv. 2009, 18, 3579–3596. [Google Scholar] [CrossRef]
  36. Nascimbene, J.; Marini, L.; Motta, R.; Nimis, P.L. Influence of tree age, tree size and crown structure on lichen communities in mature Alpine spruce forests. Biodivers. Conserv. 2009, 18, 1509–1522. [Google Scholar] [CrossRef]
  37. Ellis, C.J. Lichen epiphyte diversity: A species, community and trait-based review. Perspect. Plant Ecol. Evol. Syst. 2012, 14, 131–152. [Google Scholar] [CrossRef]
  38. Ódor, P.; Király, I.; Tinya, F.; Bortignon, F.; Nascimbene, J. Patterns and drivers of species composition of epiphytic bryophytes and lichens in managed temperate forests. For. Ecol. Manag. 2013, 306, 256–265. [Google Scholar] [CrossRef] [Green Version]
  39. Bäcklund, S.; Jönsson, M.; Strengbom, J.; Frisch, A.; Thor, G. A Pine is a Pine and a Spruce is a Spruce—The Effect of Tree Species and Stand Age on Epiphytic Lichen Communities. PLoS ONE 2016, 11, e0147004. [Google Scholar] [CrossRef]
  40. Hörsch, B. Modelling the spatial distribution of montane and subalpie forests in the central Alps using digital elevation models. Ecol. Model. 2003, 168, 267–282. [Google Scholar] [CrossRef]
  41. Hoersch, B.; Braun, G.; Schmidt, U. Relation between landform and vegetation in alpine regions of Wallis, Switzerland. A multiscale remote sensing and GIS approach. Comput. Environ. Urban 2002, 26, 113–139. [Google Scholar] [CrossRef]
  42. Saco, P.M.; Willgoose, G.R.; Hancock, G.R. Eco-geomorphology of banded vegetation patterns in arid and semi-arid regions. Hydrol. Earth Syst. Sci. 2007, 11, 1717–1730. [Google Scholar] [CrossRef] [Green Version]
  43. Engelhardt, B.M.; Weisberg, P.J.; Chambers, J.C. Influences of watershed geomorphology on extent and composition of riparian vegetation. J. Veg. Sci. 2012, 23, 127–139. [Google Scholar] [CrossRef]
  44. Jimenez-Alfaro, B.; Marceno, C.; Bueno, A.; Gavilán, R.; Obeso, J.R. Biogeographic deconstruction of alpine plant communities along altitudinal and topographic gradients. J. Veg. Sci. 2014, 25, 160–171. [Google Scholar] [CrossRef]
  45. Ardelean, I.V.; Keller, C.; Scheidegger, C. Effects of Management on Lichen Species Richness, Ecological Traits and Community Structure in the Rodnei Mountains National Park (Romania). PLoS ONE 2015, 10, e0145808. [Google Scholar] [CrossRef]
  46. Burgaz, A.R.; Fuertes, E.; Escudero, A. Ecology of cryptogamic epiphytes and their communities in deciduous forests in mediterranean Spain. Vegetatio 1994, 112, 73–86. [Google Scholar] [CrossRef]
  47. Fuertes, E.; Burgaz, A.R.; Escudero, A. Pre-climax epiphyte communities of bryophytes and lichens in Mediterranean forests from the Central Plateau (Spain). Vegetatio 1996, 123, 139–151. [Google Scholar] [CrossRef]
  48. Aragón, G.; Belinchón, R.; Martínez, I.; Prieto, M. Estimating epiphytic lichen richness by single families in Mediterranean forests. For. Ecol. Manag. 2013, 310, 187–193. [Google Scholar] [CrossRef]
  49. Christensen, S.N. The epiphytic lichen flora of Platanus orientalis stands in Greece. Willdenowia 2014, 44, 209–227. [Google Scholar] [CrossRef]
  50. Pirintsos, S.A.; Diamantopoulos, J.; Stamou, G.P. Analysis of the Distribution of Epiphytic Lichens within Homogeneous Fagus sylvatica Stands along an Altitudinal Gradient (Mount Olympos, Greece). Vegetatio 1995, 116, 33–40. [Google Scholar]
  51. Belinchón, R.; Martinez, I.; Otalora, M.A.G.; Aragón, G.; Dimas, J.; Escudero, A. Fragment quality and matrix affect epiphytic performance in a mediterranean forest landscape. Am. J. Bot. 2009, 96, 1974–1982. [Google Scholar] [CrossRef] [PubMed]
  52. Pirintsos, S.A.; Diamantopoulos, J.; Stamou, G.P. Analysis of the Vertical Distribution of Epiphytic Lichens on Pinus nigra (MountOlympos, Greece) along an Altitudinal Gradient. Vegetatio 1993, 109, 63–70. [Google Scholar] [CrossRef]
  53. Christensen, S.N. Lichen associated with Pinus nigra on Mt. Trapezitsa, Epirus, NW Greece. Acta Bot. Fenn. 1994, 150, 11–20. [Google Scholar]
  54. Guvenc, S.; Oran, S.; Ozturk, Ş. The Epiphytic Lichens on Anatolian Black Pine [Pinus nigra Arnd. subsp. pallasiana (Lamb.) Holmboe in Mt. Uludag (Bursa–Turkey). JABS 2009, 3, 157–161. [Google Scholar]
  55. Sevgi, O.; Yılmaz, O.; Carus, S.; Dündar, T.; Kavgacı, A.; Tecimen, H.B. Alaçam Dağlarındaki Karaçam Ormanlarının Yükseltiye Göre Beslenme-Büyüme Modelleri ve Odununun Teknolojik Özellikleri; TOVAG 104 O 551 Project Report; TÜBİTAK: Ankara, Turkey, 2010; p. 478.
  56. Pamay, B. Yaş Sınıfları Amenajman Metodunun Türkiye Ormanlarına Uygulanması İmkanları ve Karşılaşılan Güçlükler. İstanbul Üniv. Orman Fakültesi Derg. 1968, 18, 23–41. [Google Scholar]
  57. Atay, İ. Silvikültür II; İstanbul Üniversitesi: İstanbul, Turkey, 1990; 242p, ISBN 975-404-174-1. [Google Scholar]
  58. Purvis, O.W.; Coppins, B.J.; Hawksworth, D.L.; James, P.W.; Moore, D.M. The Lichen Flora of Great Britain and Ireland; Natural History Museum Publications: London, UK, 1992. [Google Scholar]
  59. Wirth, V. Die Flechten Baden-Württem Bergs 2 Vols; Eugen Ulmer: Stuttgart, Germany, 1995. [Google Scholar]
  60. Brodo, I.M.; Sharnoff, S.D.; Sharnoff, S. Lichens of North America; Yale University Press: New Haven/London, UK, 2001. [Google Scholar]
  61. Aiello-Lammens, M.E.; Boria, R.A.; Radosavljevic, A.; Vilela, B.; Anderson, R.P. spThin: Functions for Spatial Thinning of Species Occurrence Records for Use in Ecological Models. R Package Version 0.1.0. 2014. Available online: https://CRAN.R-project.org/package=spThin (accessed on 25 December 2018).
  62. Moning, C.; Werth, S.; Dziock, F.; Bässler, C.; Bradtka, J.; Hothorn, T.; Müller, J. Lichen diversity in temperate montane forests is influenced by forest structure more than climate. For. Ecol. Manag. 2009, 258, 745–751. [Google Scholar] [CrossRef]
  63. Goral, F.; Schellenberg, J. Function for Community Data and Ordinations. R Package Version 0.3.3. 2018. Available online: https://cran.r-project.org/web/packages/goeveg/index.html (accessed on 25 December 2018).
  64. Hijmans, R.J.; Cameron, S.E.; Parra, J.L.; Jones, P.G.; Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 2005, 25, 1965–1978. [Google Scholar] [CrossRef] [Green Version]
  65. Title, P.O.; Bemmels, J.B. ENVIREM: An expanded set of bioclimatic and topographic variables increases flexibility and improves performance of ecological niche modeling. Ecography 2018, 41, 291–307. [Google Scholar] [CrossRef]
  66. Conrad, O.; Bechtel, B.; Bock, M.; Dietrich, H.; Fischer, E.; Gerlitz, L.; Wehberg, J.; Wichmann, V.; Böhner, J. System for automated geoscientific analyses (SAGA) v. 2.1.4. 475. Geosci. Model Dev. 2015, 8, 1991–2007. [Google Scholar] [CrossRef]
  67. Kudrnovsky, H. Geographic Resources Analysis Support System (GRASS) Software. 2018. Available online: https://grass.osgeo.org/grass7/manuals/addons/i.segment.stats.html (accessed on 25 December 2018).
  68. Gauslaa, Y. The Lobarion and Epiphytic Community of Ancient Forests Threatened by Acid Rain. Lichenologist 1995, 27, 59–76. [Google Scholar]
  69. Poikolainen, J. Mosses, Epiphytic Lichens and Tree Barks as Biomonitors for Air Pollutants—Specifically for Heavy Metals in Regional Surveys; The Finnish Forest Research Institute, Oulu University Press: Oulu, Finland, 2004.
  70. Bremmer, J.M.; Mulvaney, C.S. Nitrogen–Total. In Methods of Soil Analysis. Part 2. Chemical and Microbiological Properties, 2nd ed.; Page, A.L., Miller, R.H., Keeney, D.R., Eds.; Soil Science Society of America Inc.: Madison, WI, USA, 1982; pp. 595–624. [Google Scholar]
  71. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. Vegan: Community Ecology Package. R Package Version 2.5-4. 2019. Available online: https://CRAN.R-project.org/package=vegan (accessed on 25 December 2019).
  72. Dufrêne, M.; Legendre, P. Species assemblages and indicator species: The need for a flexible asymmetrical approach. Ecol. Monogr. 1997, 67, 345–366. [Google Scholar] [CrossRef]
  73. De Cáceres, M.; Legendre, P.; Moretti, M. Improving indicator species analysis by combining groups of sites. Oikos 2010, 119, 1674–1684. [Google Scholar] [CrossRef]
  74. De Cáceres, M.; Legendre, P. Associations between species and groups of sites: Indices and statistical inference. Ecology 2009, 119, 1674–1684. [Google Scholar] [CrossRef]
  75. Revelle, W. psych: Procedures for Personality and Psychological Research; Northwestern University: Evanston, IL, USA, 2018. [Google Scholar]
  76. R Development Core Team. A language and Environment for Statistical Computing. In R Foundation for Statistical Computing; R Development Core Team: Vienna, Australia, 2018; Available online: http://www.R-project.org (accessed on 10 June 2018).
  77. Nascimbene, J.; Marini, L.; Nimis, P.L. Influence of tree species on epiphytic macrolichens in temperate mixed forests of northern Italy. Can. J. For. Res. 2009, 39, 785–791. [Google Scholar] [CrossRef]
  78. Çobanoğlu, G.; Sevgi, E.; Sevgi, O. Lichen Mycota along Uludağ Fir (Abies bornmuelleriana MATTF); University of Craiova Seria: Craiova, Romania, 2008; Volume 13, pp. 15–19. [Google Scholar]
  79. Nascimbene, J.; Marini, L.; Ódor, P. Drivers of lichen species richness at multiple spatial scales in temperate forests. Plant Ecol. Div. 2012, 5, 355–363. [Google Scholar] [CrossRef] [Green Version]
  80. Király, I.; Nascimbene, J.; Tinya, F.; Ódor, P. Factors influencing epiphytic bryophyte and lichen species richness at different spatial scales in managed temperate forests. Biodivers. Conserv. 2013, 22, 209–223. [Google Scholar] [CrossRef]
  81. Kreylinga, J.; Wiesenbergb, G.L.B.; Thiela, D.; Wohlfarta, C.; Huberc, G.; Walterd, J.; Jentschd, A.; Konnertc, M.; Beierkuhnleina, C. Cold hardiness of Pinus nigra Arnold as influenced by geographic origin, warming, and extreme summer drought. Environ. Exp. Bot. 2012, 78, 99–108. [Google Scholar] [CrossRef]
  82. Enescu, C.M.; de Rigo, D.; Caudullo, G.; Mauri, A.; Houston Durrant, T. Pinus nigra in Europe: Distribution, habitat, usage and threats. In European Atlas of Forest Tree Species; San-Miguel-Ayanz, J., de Rigo, D., Caudullo, G., Houston Durrant, T., Mauri, A., Eds.; European Commission: Luxembourg, 2016; Volume 6, pp. 126–127. [Google Scholar]
  83. Cobanoglu, G.; Sevgi, O. Analysis of the Distribution of Epiphytic Lichens on Cedrus libani in Elmali Research Forest (Antalya, Turkey). J. Environ. Biol. 2009, 30, 205–212. [Google Scholar]
  84. Sevgi, O.; Tecimen, H.B.; Carus, S.; Akburak, S.; Çaksir, E. Some bark characteristics of Black Pine (Arnold.) and their variation throughout the tree height Pinus nigra. J. Environ. Biol. 2016, 37, 1347–1354. [Google Scholar]
  85. Marmor, L.; Randlane, T.E. Effects of road traffic on bark pH and epiphytic lichens in Tallinn. Folia Cryptog. Est. 2007, 43, 23–37. [Google Scholar]
  86. Selva, S.B. Lichen diversity and stand continuity in the northern hardwoods and spruce-fir forests in northern New England and western New Brunswick. Bryologist 1994, 97, 424–429. [Google Scholar] [CrossRef]
  87. Svoboda, D.; Peksa, O.; Vesela, J. Epiphytic lichen diversity in central European oak forests: Assessment of the effects of natural environmental factors and human influences. Environ. Pollut. 2010, 158, 812–819. [Google Scholar] [CrossRef]
  88. Rubio-Salcedo, M.; Psormas, A.; Prieto, M.; Zimmermann, N.E.; Martinez, I. Case study of the implications of climate change for lichen diversity and distribution. Biodivers. Conserv. 2017, 26, 1121–1141. [Google Scholar] [CrossRef]
  89. Giorani, P.; Incerti, G. The influence of climate on the distribution of lichens: A case study in a borderline area (Liguria, NW Italy). Plant Ecol. 2008, 195, 257–272. [Google Scholar] [CrossRef]
  90. Jovan, S.; McCune, B. Regional Variation in Epiphytic Macrolichen Communities in Northern and Central California Forests. Bryologist 2004, 107, 328–339. [Google Scholar] [CrossRef]
  91. Hauck, M. Site factors controlling epiphytic lichen abundance in northern coniferous forests. Flora 2011, 206, 81–90. [Google Scholar] [CrossRef]
  92. Dymytrova, L.; Nadyeina, O.; Hobi, M.L.; Scheidegger, C. Topographic and forest-stand variables determining epiphytic lichen diversity in the primeval beech forest in the Ukrainian Carpathians. Biodivers. Conserv. 2014, 23, 1367–1394. [Google Scholar] [CrossRef]
  93. Nascimbene, J.; Marini, L. Epiphytic lichen diversity along elevational gradients: Biological traits reveal a complex response to water and energy. J. Biogeogr. 2015, 42, 1222–1232. [Google Scholar] [CrossRef]
  94. Aussenac, G. Interactions between forest stands and microclimate: Ecophysiological aspects and consequences for silviculture. Ann. For. Sci. 2000, 57, 287–301. [Google Scholar] [CrossRef] [Green Version]
  95. Dettki, H.; Klintberg, P.; Esseen, P.A. Are epiphytic lichens in young forests limited by local dispersal? Écoscience 2000, 7, 317–325. [Google Scholar] [CrossRef]
  96. Aude, E.; Poulsen, R.S. Influence of management on the species composition of epiphytic cryptogams in Danish Fagus forests. Appl. Veg. Sci. 2000, 3, 81–88. [Google Scholar] [CrossRef]
  97. Mikhailova, I.; Trubina, M.; Vorobeichik, E.; Scheidegger, C. Influence of environmental factors on the local-scale distribution of cyanobacterial lichens: Case study in the North Urals, Russia. Folia Cryptogam. Est. 2005, 41, 45–54. [Google Scholar]
  98. Löbel, S.; Snäll, T.; Rydin, H. Species richness patterns and metapopulation processes-evidence from epiphyte communities in boreo-nemoral forests. Ecography 2006, 29, 169–182. [Google Scholar] [CrossRef]
  99. Fritz, Ö.; Gustafsson, L.; Larsson, K. Does forest continuity matter in conservation? A study of epiphytic lichens and bryophytes in beech forests of southern Sweden. Biol. Conserv. 2008, 141, 655–668. [Google Scholar] [CrossRef]
  100. Mežaka, A.; Brūmelis, G.; Piterāns, A. Tree and stand-scale factors affecting richness and composition of epiphytic bryophytes and lichens in deciduous woodland key habitats. Biodivers. Conserv. 2012, 21, 3221–3241. [Google Scholar] [CrossRef]
Diversity 11 00059 g001 550
Figure 1. Location of the three study regions (1: Dursunbey, 2: Sındırgı, 3: Simav) of Pinus nigra forest stands in the Alaçam Mountains (top). The richness of each sampling plot is represented by the sizes and colors of circles, which represent different study regions (bottom).
Figure 1. Location of the three study regions (1: Dursunbey, 2: Sındırgı, 3: Simav) of Pinus nigra forest stands in the Alaçam Mountains (top). The richness of each sampling plot is represented by the sizes and colors of circles, which represent different study regions (bottom).
Diversity 11 00059 g001
Diversity 11 00059 g002 550
Figure 2. Variation partitioning of epiphytic lichen species composition among the main variable groups (values <0 not shown): X1 = bioclimatic, X2 = topographical variables, X3 = parent rock, and X4 = stand structure variables: (A) LPA, (B) LFQ, (C) GLPA, and (D) GLFQ.
Figure 2. Variation partitioning of epiphytic lichen species composition among the main variable groups (values <0 not shown): X1 = bioclimatic, X2 = topographical variables, X3 = parent rock, and X4 = stand structure variables: (A) LPA, (B) LFQ, (C) GLPA, and (D) GLFQ.
Diversity 11 00059 g002
Diversity 11 00059 g003 550
Figure 3. Ordination biplot of the stand level redundancy analysis of epiphytic lichen species. Distribution related to environmental variables: (A) LPA, (B) LFQ, (C) GLPA, and (D) GLFQ. Species are indicated by red labels and abbreviated with cepnames. Explanatory variables are indicated using blue arrows.
Figure 3. Ordination biplot of the stand level redundancy analysis of epiphytic lichen species. Distribution related to environmental variables: (A) LPA, (B) LFQ, (C) GLPA, and (D) GLFQ. Species are indicated by red labels and abbreviated with cepnames. Explanatory variables are indicated using blue arrows.
Diversity 11 00059 g003
Diversity 11 00059 g004 550
Figure 4. Species GLM response curves of most of the abundant epiphytic lichen species based on: (A) Elevation, (B) mDbh, (C) Bio15, (D) AridxThorn, (E) mCH, and (F) CtchSlp.
Figure 4. Species GLM response curves of most of the abundant epiphytic lichen species based on: (A) Elevation, (B) mDbh, (C) Bio15, (D) AridxThorn, (E) mCH, and (F) CtchSlp.
Diversity 11 00059 g004
Diversity 11 00059 g005 550
Figure 5. Pearson pairwise correlation analysis between epiphytic lichen species diversity and tree properties.
Figure 5. Pearson pairwise correlation analysis between epiphytic lichen species diversity and tree properties.
Diversity 11 00059 g005
Table
Table 1. Environmental variable groups used in the study area.
Table 1. Environmental variable groups used in the study area.
Explanatory VariableNumber of VariablesSource of the Data
Bioclimatic19 bioclimatic and 18 ENVIREM variables, 37 in total[64,65]
Topographic 40Derived from Aster digital elevation model (DEM) using Saga GIS software [66]
Stand structure10Field survey
Parent rock1 categorical variable1:25,000 parent rock maps
Tree properties6Field survey and laboratory work
Table
Table 2. Results of the indicator epiphytic lichen species analysis for region and age classes.
Table 2. Results of the indicator epiphytic lichen species analysis for region and age classes.
GroupIndicator SpeciesA#B##IndValp-Value
Geographic Region (LPA data)
SimavBryoria implexa0.68420.50000.5850.022 *
Letharia vulpina0.65000.50000.5700.042 *
Dursunbey and SındırgıPlatismatia glauca0.91350.88100.8970.002 **
Geographic Region (LFQ data)
SimavUsnea subfloridana0.73010.66670.6980.008 **
Dursunbey and SındırgıPlatismatia glauca0.96930.88100.9240.001 ***
Age Class (LPA data)
Old-growthHypocenomyce scalaris0.66520.83330.7450.003 **
#A refers to the probability of the suitability of studied site to the target site group, where the concerned species is present. ##B refers to probability of presence of the concerned species in the target site group. Significance codes: *** p < 0.001, ** p < 0.01, and * p < 0.05. LPA: presence–absence data; LFQ: frequency data.
Table
Table 3. The variance of the significant explanatory variables used in the stand-level redundancy analysis (RDA) with LPA, GLPA, LFQ, and GLFQ data sets.
Table 3. The variance of the significant explanatory variables used in the stand-level redundancy analysis (RDA) with LPA, GLPA, LFQ, and GLFQ data sets.
VariableVarianceFp-Value
LPA (R2 = 0.365, adjusted R2 = 0.254)
Elev0.26454.44090.001
maxH0.18663.13170.002
Bio150.14342.40760.003
mDbh0.12062.02520.016
TnCurv0.13042.18930.006
CtchSlp0.12062.02450.009
Txtr0.12712.13420.011
hecTBA0.11831.98650.012
AridxThorn0.09621.61560.054
MonthTemp100.09711.63060.049
GLPA (R2 = 0.389 adjusted R2 = 0.306)
Elev0.15235.94350.001
mDbh0.11054.31460.002
mCH0.08603.35870.001
MonthTemp100.06212.42530.015
Bio150.06042.35770.026
DwnCurv0.05872.29180.032
AridxThorn0.05922.30920.026
LFQ (R2 = 0.397, adjusted R2 = 0.234)
Elev0.09935.98360.001
mCH0.06313.79760.001
mDbh0.04762.86510.002
Parent rock0.10611.59780.008
HI0.03051.83740.050
sdCH0.02981.79670.045
CtchSlp0.02861.72330.068
GLFQ (R2 = 0.342, adjusted R2 = 0.245)
Elev0.05616.93710.001
mDbh0.03404.19710.002
Bio150.02913.59640.006
DwnCurv0.01832.26530.030
HI0.01732.14350.049
mCH0.01722.12300.053
Table
Table 4. The linear regression model parameters between tree properties and lichen richness.
Table 4. The linear regression model parameters between tree properties and lichen richness.
ModelAdjusted R2F-Statisticp-Value
Richness = 4.687085 + 0.010691 × Dbh0.00282.06240.152
Richness = 3.80421 + 0.06814 × Height0.057723.76650.0001 ***
Richness = 5.1316 − 0.5181 × BarkN−0.00210.28650.593
Richness = 12.0355 − 1.7572 × Bark_pH0.051520.87210.0001 ***
Richness = 4.63093 + 0.09175 × DBT0.01165.360.0211*
Notes: *** p < 0.001, ** p < 0.01, and * p < 0.05.

Share and Cite

MDPI and ACS Style

Sevgi, E.; Yılmaz, O.Y.; Çobanoğlu Özyiğitoğlu, G.; Tecimen, H.B.; Sevgi, O. Factors Influencing Epiphytic Lichen Species Distribution in a Managed Mediterranean Pinus nigra Arnold Forest. Diversity 2019, 11, 59. https://doi.org/10.3390/d11040059

AMA Style

Sevgi E, Yılmaz OY, Çobanoğlu Özyiğitoğlu G, Tecimen HB, Sevgi O. Factors Influencing Epiphytic Lichen Species Distribution in a Managed Mediterranean Pinus nigra Arnold Forest. Diversity. 2019; 11(4):59. https://doi.org/10.3390/d11040059

Chicago/Turabian Style

Sevgi, Ece, Osman Yalçın Yılmaz, Gülşah Çobanoğlu Özyiğitoğlu, Hüseyin Barış Tecimen, and Orhan Sevgi. 2019. "Factors Influencing Epiphytic Lichen Species Distribution in a Managed Mediterranean Pinus nigra Arnold Forest" Diversity 11, no. 4: 59. https://doi.org/10.3390/d11040059

APA Style

Sevgi, E., Yılmaz, O. Y., Çobanoğlu Özyiğitoğlu, G., Tecimen, H. B., & Sevgi, O. (2019). Factors Influencing Epiphytic Lichen Species Distribution in a Managed Mediterranean Pinus nigra Arnold Forest. Diversity, 11(4), 59. https://doi.org/10.3390/d11040059

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop