Next Article in Journal
Regeneration Dynamics and Development of Seedlings in Sessile Oak Forests in Relation to the Light Availability and Competing Vegetation
Previous Article in Journal
Deadwood, Soil and Carabid Beetle-Based Interaction Networks—An Initial Case Study from Montane Coniferous Forests in Poland
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reduced-Impact Logging Maintain High Moss Diversity in Temperate Forests

by
Enrique Hernández-Rodríguez
1,2,*,
Luis H. Escalera-Vázquez
3,
Deneb García-Ávila
3,
Miguel Montoro Girona
2,4,5 and
Eduardo Mendoza
1,*
1
Instituto de Investigaciones Sobre los Recursos Naturales, Universidad Michoacana de San Nicolás de Hidalgo, Av. San Juanito Itzícuaro s/n, Morelia 58337, Michoacán, Mexico
2
Institut de Recherche sur les Forêts, Université du Québec en Abitibi Témiscamingue, 445 Boul. de l’Université, Rouyn-Noranda, QC J9X5E4, Canada
3
Facultad de Biología, Universidad Michoacana de San Nicolás de Hidalgo, Av. Gral. Francisco J. Múgica s/n A-1, Morelia 58030, Michoacán, Mexico
4
Groupe de Recherche en Écologie de la MRC-Abitibi (GREMA), Campus of Amos, Université du Québec en Abitibi Témiscamingue, 341 Rue Principal Nord, Amos, QC J9T2L8, Canada
5
Restoration Ecology Group, Department of Wildlife, Fish and Environmental Studies, Swedish University of Agricultural Sciences, SE-901 83 Umeå, Sweden
*
Authors to whom correspondence should be addressed.
Forests 2021, 12(4), 383; https://doi.org/10.3390/f12040383
Submission received: 26 January 2021 / Revised: 14 March 2021 / Accepted: 14 March 2021 / Published: 24 March 2021
(This article belongs to the Section Forest Ecology and Management)

Abstract

:
Forestry harvesting represents an important economic activity around the world. Habitat degradation due to forest harvesting contributes to biodiversity loss; therefore, it is necessary to implement logging management aimed at reducing its impact. Forest management by reduce-impact logging (RIL) involves cutting trees following regulations focused on diminishing the impact on biodiversity by following harvesting plans based on forestry inventories and participation of trained workers. In Mexico, RIL is applied mainly in temperate habitats and its effectiveness has been assessed based on vascular plants. In this study, we analyzed the diversity and community structure of terrestrial and epiphytic mosses in managed (sites number = 3) and conserved (sites number = 3) sites in the temperate forest of Sierra Juárez, Oaxaca, Mexico. Likewise, we evaluated the potential function of mosses as indicators of habitat degradation. Environmental variables were also quantified at local (canopy coverage, altitude, daily temperature, and light) and regional (total annual rainfall, orientation, and slope) scales to evaluate potential relationships with the community and species diversity. We documented 70 mosses species with a diversity (alfa, beta) and community structure similar between managed and conserved sites. For terrestrial mosses, we found marginal differences in their communities, likely related to species coverture variation in managed sites. The diversity and community structure epiphytic mosses were not statistically different in managed and conserved sites. Only the daily variation in light intensity was positively related to the variation of alpha diversity of epiphytic mosses. The species Dicranum sumichrastii Duby and Leptodontium viticulosoides (P. Beauv.) Wijk & Margad. can be considered as ecological indicators for conserved and managed sites, respectively, likely due to their relationship with light and humidity conditions. Our results suggest that that forest management by RIL could be considered as a promising tool to balance timber production and moss diversity.

1. Introduction

Biodiversity is highly threatened by anthropogenic disturbances such as deforestation, land change, soil pollution, and habitat fragmentation around the world [1,2,3,4,5]. Forestry logging is a human activity that affects biodiversity patterns, specially reducing the number of species depending of the disturbance intensity [6,7]. To reduce logging impacts, sustainable forest management represents a tool to develop productive practices and conservation including the protection of species and ecosystems services [8,9,10].
Forest management represents a series of techniques that enable productive practices and depending on the wood harvest intensity and the impact on species richness, six standard silvicultural systems globally can be distinguished, from the least to the most harmful: selection and retention systems, reduced-impact logging, selective logging, clear-cutting, and timber plantations [11].
Reduced-impact logging (RIL) is a silvicultural practice characterized by the selective extraction of trees from previously defined areas to reduce forest damage and protecting of riparian vegetation [12]. Trained workers select and extract individual trees carefully, according to a harvest plan and guidelines reducing the negative impact of tree felling and hauling on the remaining forest [13,14,15,16,17]. RIL has been implemented in the forests of Europe and Asia, showing that silviculture systems reduce impact on local and regional species richness [11,18,19]. In southern Mexico, this RIL management is developed by local communities in tropical and temperate forests [20,21,22,23]. In contrast with forest management strategies such as clear-cutting, where entire tree stands are removed and biodiversity is highly impacted [11,24], the RIL generates small clearings reducing habitat disturbance effects [25,26,27].
Evidence indicates that RIL has a low impact on the species richness of birds, insects, mammals, and vascular plants in tropical regions [28]. However, studies exploring the effects of RIL on native vegetation are focused on vascular plants (e.g., [29,30,31]), whereas effects on non-vascular flora, like mosses have not been assessed. Mosses are associated to specific microenvironmental conditions (e.g., low light, high humidity), structure of the habitat (e.g., host tree species), and type of substrates to establish and thrive (e.g., advanced decaying wood, rocks, soil) [32,33]. These micro conditions are related to richness (alpha diversity) and species turnover (beta diversity) of mosses [34]. Thus, even slight changes in environmental conditions due to RIL could negatively affect the diversity and composition of mosses communities [11,19,35]. Due to this characteristic, mosses can function as indicators of the state of forest where RIL is applied [36]. The conservation of mosses is essential because they play an important role in biochemical cycles and ecosystem dynamics due to the fact they store up to 1400% of their dry weight in water [37], fix up to 400 mg N m−2 year−1 [38], and provide shelter to protozoa, invertebrates, and arthropods [39]. Thus, evaluating the impact of RIL on mosses is a crucial issue to ensure biodiversity preservation in the long-term.
In this study, we evaluated the impact of RIL on moss diversity (terrestrial and epiphytic) in a pine-oak temperate forest in the state of Oaxaca, Mexico. We compared the diversity of mosses from three sites under RIL management with that of three conserved sites. Our hypotheses were: (i) the disturbance caused by RIL will affect the diversity (alpha and beta) of both groups and their community structure due to changes in microhabitat conditions (high light and temperature), and (ii) some moss species will show habitat preferences related to RIL and conserved forest sites.

2. Materials and Methods

2.1. Study Area

Our study area is located in the pine-oak forest of the Sierra Juárez (Sierra Madre) of Oaxaca state, south of Mexico [40] (Figure 1). The climate is temperate sub-humid with rains during the summer; the annual rainfall ranges between 1000 and 2000 mm, and the average annual temperature is 16 °C [41]. Pine-oak vegetation occurs in an altitudinal range from 2000 to 2800 m on humic acrisol soils [42,43]. Forest stands have species of pines (Pinus ayacahuite C. Ehrenb. ex Schltdl., P. leiophylla Schiede ex Schltdl. & Cham., P. patula Schltdl. & Cham., P. pseudostrobus Brongn.), oaks (Quercus crassifolia Bonpl., Q. elliptica Née, Q. laeta Liebm.), and bushes such as Arbutus xalapensis Kunth and Clethra L. spp. The most abundant genus is Pinus followed by Quercus L. [44] and present high associations with mosses like Brachythecium Schimp., Bryum Hedw., Entodon Müll. Hal., and Leptodontium (Müll. Hal.) Hampe ex Lindb. spp. [45,46]. In the forest understory, occur shrubs species including Arctostaphylos pungens Kunth, Gaultheria acuminata Schltdl. & Cham., Litsea glaucescens Kunth and herbaceous species including Alchemilla pectinata Kunth, Stevia lucida Lag., and Dryopteris Adans. spp. Moreover, epiphytes are represented by orchids, bromeliads, and mosses including species of the genus Macromitrium Brid., Meteorium Dozy & Molk., Neckera Hedw. and Orthotrichum Hedw. [44,46]. The area harbors a high bird diversity and threatened mammals species, and constitutes a center of vascular plant endemism in the country [40,47,48,49]. Approximately 74,240 ha of the Sierra Juárez are under management for timber production, almost all pine or pine-oak forest [50].

2.2. Experimental Design

This study was conducted in six sites, three old forest dedicated to conservation (C) in which no logging activities are performed (untreated control), and three managed (M) where the RIL system is applied. The conserved sites differed in size C1 = 7134 ha, C2 = 12,936 ha, C3 = 230 ha. Due to the forest planning, there was variation in the time RIL management was implemented in each site. In M1, logging activities began in 2004 on 107 ha, and about in one year every ten years, diseased and old trees are removed (logging intensity). In M2, RIL started in 2012 on 13.5 ha, and in M3 in 1995 on 10 ha, in both cases, logging is conducted during periods of 10 years every 50 years in which workers prune trees branches and shrubs to reduced plant competition, favoring the increment in diameter and height of pine trees. Time since the last disturbance in each managed site is two years for M1, four for M2, and nine for M3. The tree number removed per hectare in each site during logging period is under five percentage from total trees. Additionally, RIL in the region involves forest inventories to define how the logging procedures are conducted by local trained workers and the preparation of harvest plans to minimizing the impact on vegetation [51]. Total area from these three RIL sites represented around the 0.14% from the pine-oak area in the Sierra Juarez [52].

2.3. Sampling and Data Compilation

Sampling was conducted during July and August 2016. In each study site, ten sampling quadrats of 1 m2 were randomly established in an area of 5 ha (60 quadrats in total). In each quadrat, terrestrial mosses were those growing in the ground on bare soil, humus, and decaying wood; epiphytic mosses included those on the trunk of the closest tree (diameter at breast height ≥ 30 cm) of the genus Quercus to each of the quadrats. To sample epiphytic mosses, we place a 20 cm2 quadrat on the tree trunk at the height of 1.6 m [45]. The location of the quadrat (N, S, E, or O) was set randomly [45].
To estimate the coverture of the different species of terrestrial and epiphytic, mosses quadrats of 1 and 20 cm2 were divided using a grid of 100 units. We draw the area’s silhouette covered by each species on a plastic sheet placed over the quadrat. The area covered by each species was measured as the percentage of units occupied in the grid. A sample of each moss was collected for taxonomic determination [45,53]. Species were identified using different taxonomical guides [54,55,56,57,58,59,60] and following the nomenclature in the Catalogue of Neotropical Mosses (LATMOSS) [61]. All the collected specimens were deposited in the Herbario de la Escuela de Biología of the Universidad Michoacana de San Nicolás de Hidalgo (EBUM) (Appendix A Table A1).
To assess the relationship between environmental variables with moss diversity, local-scale variables (LSV) and regional-scale variables (RSV) were measured. The LSV included canopy coverage and altitude of each quadrat, as well as daily temperature and light. The coverage of the canopy was measured in the center of each 1 m2 quadrat using a convex spherical densimeter and averaging the percentages of coverage recorded when facing each cardinal point. The temperature and light were recorded in intervals of 15 min during 44 days from 3 March to 17 April 2017 (96 records per day), at each site using two automatic recorders (HOBOS UA-002-64, Onset Computer Corporation, MA, USA). The two recorders were placed in the center of each 5 ha study site separated by 223 m. The altitude was registered with a handheld GPS (GPSMAP 64s, GARMIN, KN, USA). RSV variables included total annual rainfall, orientation, and slope of each site. These were obtained from the Instituto Nacional de Estadística y Geografía (INEGI) databases with a 30 m of resolution [62] and processed using Quantum Geographic Information System (QGIS, Open Source Geospatial Foundation Project) [63]. We obtained the information overlapping the INEGI environmental shapefiles with each site location.

2.4. Data Analyses

To assess sampling completeness in each site, we calculated the Chao 2 species richness index (alfa diversity) [64]. To compare observed and estimated richness, rarefaction/extrapolation (R/E) curves were performed with 95% confidence intervals. We calculated the percentage that the observed species represented in the estimated richness. To assess alpha diversity and perform the R/E curves, the method proposed by Chao [64] based on the Hill number of zero-order (q = 0) was used. Beta diversity was estimated to evaluate how the species replacement (βrepl) and species richness (βrich) affected total beta diversity (βtotal) [65] within and between sites.
We performed a non-metric multidimensional scaling analysis (NMDS), using Bray–Curtis dissimilarities matrices, analysis of similarities (ANOSIM test), and rank abundance curves [66] to compare the moss communities between forest types. For NMDS and ANOSIM, we used data to the quadrant level. We discarded the species occurring in less than 5% of the quadrats to focus only on those most representative and because rare species can have a high influence in community simulation tests [67,68,69]. The percentage of coverture was transformed (by applying the square root and arcsine) prior to analysis to have data with a normal distribution [69]. We conducted the ANOSIM test with 999 permutations. To build rank-abundance curves, we used the coverture data of each species transformed to the logarithm. To evaluate the moss species association with forest types, we applied the indicator value (IndVal index) [70] using the quadrants data per site. This method considers the specificity (A): the degree of association of a species with a site type, and fidelity (B): the degree to which, given that there is a forest conserved or under management, a certain moss species occurs. The statistical significance of each relationship was evaluated with 999 permutations [70].
Finally, to analyze the relationship between abiotic environmental variables and moss diversity (richness values, mean, standard deviation), the averages of each LSV and RSV variable at site level were calculated. Then, their relationship was analyzed through a correlation analysis using the Pearson correlation coefficient [71]. Additionally, the averages values per site were used to evaluate differences in LSV and RSV between management and conservation conditions applying t-student and Wilcoxon test.
All analyses were conducted using R (Version 3.5.0, R Foundation for Statistical Computing, Vienna, Austria) [72]. For richness estimation and rarefaction/extrapolation curves, we used iNEXT 2.0.14 package [73]; for the analysis of beta diversity, the BAT package [74]; for the NMDS and ANOSIM analyzes, we used vegan 2.4-4 package [75]; for the calculation of IndVal index, we used Indicspecies 1.7.1 package [76]; to conduct the correlation analysis, the Performance Analytics package [77].

3. Results

3.1. Alpha and Beta Diversity and Community Structure of Mosses

We obtained 251 samples, which included 70 species of 51 genera and 31 families of mosses (Appendix A Table A1). Fifteen mosses species were terrestrials and epiphytic, 45 were exclusively terrestrial and exclusively 20 epiphytic. Based on the Chao-2 species richness estimator, the sampling’s completeness was, on average of 50% for terrestrial and 48% for epiphytic mosses per site, respectively. Managed sites had, on average, 12 ± 2 species (mean ± SD) of terrestrial mosses, while conserved sites had 18 ± 1 species. In comparison, there was an average of 10 ± 3 species and 10 ± 1 species epiphytic mosses in management and conservation sites (Table 1). Thus, the average of epiphytic moss richness was similar in conserved and management sites.
The R/E curves showed no significant differences (based in the lack of overlap between 95% confidence intervals) between managed and conserved sites for terrestrial and epiphytic mosses. In both groups of mosses, species richness trends were maintained when extrapolating to greater sampling (Figure 2).
Beta diversity values for terrestrial and epiphyte mosses are similar in conserved and managed sites. In conserved sites and for terrestrial mosses, we estimated a total β diversity between 71 and 80%. Between 64 and 79% of β diversity was explained by species replacement (βrepl) and only between 0 to 7% was explained by differences in richness (βrich). In comparison, within managed sites, we estimate a total β between 69 and 81%. From this, between 53 and 67% was explained by βrepl and between 6 to 21% by βrich. Thus, in both forest types, βtotal diversity for terrestrial mosses was mainly associated with βrepl (Table 2). In conserved sites for epiphytic mosses, the βtotal ranged between 62 and 85%. βrepl explained 31 to 80% of this β diversity, and βrich between 5 and 31%. In managed sites, the total β ranged between 67 and 100%, whereas βrepl represented 53 to 90%, and βrich between 10 to 13%. In summary, in both forest types, the βtotal diversity was mostly explained by βrepl, and in conserved sites, βrich was higher than in managed sites (Table 2). Between forest types, the total β diversity of terrestrial and epiphytic mosses varied between 76 to 74%, respectively. For terrestrial mosses, 52% of this β diversity was explained by βrepl, whereas the remaining 24% was explained by βrich. In the case of epiphytic mosses, βtotal was defined entirely by βrepl (Table 2).
For terrestrial and epiphytic mosses, the rank-abundance curves had at each site in both forest types a steep slope. This suggests a low uniformity in the community dominance because a few species have a high abundance. However, the most abundant species, and those with low abundance, differed between sites. For terrestrial mosses, there were 1–5 species with the lowest abundance in the managed forest, while in the conserved forest, there are 6–7 species. There were a few dominant terrestrial moss species, most of them with high abundance, while in epiphytic mosses, there are more dominant species, but few of them had low abundance (Figure 3). Common terrestrial species in both forest types included Braunia squarrulosa (Hampe) Müll. Hal., Bryum billarderii Schwägr., Hypnum amabile (Mitt.) Hampe, Mittenothamnium reptans (Hedw.) Cardot, Sematophyllum swartzii (Schwägr.) W.H. Welch & H.A. Crum, Thuidium delicatulum (Hedw.) Schimp., Trichostomum brachydontium Bruch, and Zygodon ehrenbergiI Müll. Hal. In epiphytic mosses, common species were Holomitrium pulchellum Mitt., Leptodontium viticulosoides (P. Beauv.) Wijk & Margad., Leucodon curvirostris Hampe, Neckera chlorocaulis Müll. Hal., and Zygodon viridissimus (Dicks.) Brid.
The NMDS ordination procedure suggests that the community of terrestrial and epiphytic mosses are similar between managed and conserved sites (Figure 4). However, for the terrestrial moss community, there are higher R values and significant differences (ANOSIM R = 0.14; p = 0.001) in comparison to the epiphytic moss community (ANOSIM R = 0.02; p = 0.199). Differences in terrestrial moss community can be related to higher variation in species coverture in some managed sites. For epiphytic mosses, some quadrants in management sites are distant from the rest due to differences in species composition but without significant differences (Figure 4). The low-stress values in both cases indicate stability in the ordination configuration [78].

3.2. Indicator Species and the Environment Influence on Mosses Diversity

The IndVal index identified the terrestrial mosses Dicranum sumichrasti and Leptodontium viticulosoides as potential indicator species for conservation and management sites, respectively (p < 0.05). D. sumichrasti is an indicator of sites conserved with high specificity (A = 100%) but low fidelity (B = 30%). In contrast, L. viticulosoides is an indicator of sites under management with high specificity (A = 84%) but intermediate fidelity (B = 53%). In the case of the epiphytic mosses, no indicator species were found.
We did not find significant differences in local-scale variables (LSV) and regional-scale variables (RSV) between managed and conserved sites (Table 3). Of the LSV and RSV analyzed by correlation analysis, RSV did not show an effect on mosses diversity. Nevertheless, we found a positive correlation between the variation in alpha diversity of epiphytic mosses and the variation in light availability (LSV) (r = 0.83, p < 0.05) (Appendix A Table A2).

4. Discussion

The evaluation of silvicultural practices to reach sustainable forest management strategies is a major priority in forest ecosystems around the world [79,80,81,82,83]. We are living a critical moment in the adaptation of harvesting methods to maintain the forest biodiversity [84,85,86,87,88,89]. Here, we present the first study evaluating the effect of forest management by RIL on moss diversity in a temperate forest in Mexico. Based on our results, the RIL can keep a great diversity (alpha and beta), but also the community structure, of this group of plants in managed patches.

4.1. Diversity and Community Moss Structure

Our study demonstrated that RIL helps to conserve terrestrial and epiphytic mosses diversity in temperate forest. These results contrast with studies showing a decrease in moss diversity as a consequence of clear-cutting [90,91] and selective logging [92]. Our results also differ from those where other management systems such as shelterwood logging [93], clear-cutting [94], or tree retention [95] are applied and species richness is found to be similar between managed and protected sites but where species composition differs due to changes in the forest structure. These differences are related with disturbance severity of each management system. Compared with those system RIL harvesting areas of Sierra Juárez, they do not show severe impacts on the diversity and moss community structure base on their similarity with conserved sites. However, we found a marginal effect on the terrestrial moss community (Figure 4). In our study, the richness and community structure of terrestrial mosses are similar (Figure 2, Figure 3 and Figure 4; Table 2), but our ANOSIM analysis reveals differences, possibly related to species coverture variation in RIL management sites. This finding highlights the importance of including other bryophyte community attributes, like their coverture, to know disturbance effects that could be overlooked.
The similarity in the richness and community moss structure between conserved and managed forests is likely related to maintaining microhabitat conditions [96,97,98]. These micro conditions are, for example, the humidity and bark substrate provided by Quercus trees keeping and the no severe disturbance on terrestrial substrates in RIL areas [99]. In this way, general characteristics of the structure of moss communities remained similar for both groups in terms of the frequency and composition of dominant and rare species.
Based on our results, RIL does not affect the community level, but it could have effects on some particular species. This possibility is reinforced if we consider that the inventory completeness (alpha diversity) at each site was approximately 50%. In the case of terrestrial mosses in the three sites, rank-abundance curves indicated that RIL could affect some of the less abundant species (or rare local species [100]) because of their specific requirements [101] compared to the most abundant ones (Figure 3a). These requirements are represented, for example, by availability of suitable habitats as deadwood with different decay stages, humidity, and special edaphic conditions [102,103]. Another factor in explaining these results is the time since the last disturbance. For example, site M1 with a recent disturbance (two years) had fewer species observed and expected than other managed sites (Table 1). Some of these terrestrial species do not present in managed sites including Dicranum scoparium Hedw., Heterophyllum affine (Hook. in Kunth) Fleisch., Prionodon luteovirens (Taylor) Mitt., and Zygodon ehrenbergii Müll. Hal. Among epiphytes, the M3 site was different from the rest of the managed and conserved sites which can be related to its stand conditions. This site is the oldest, with less logging area, and presents the highest humidity compared with the rest of the evaluated management sites. An old forest keeps conditions such as humidity and unique subtracts necessary for some species [104,105]. Other studies in boreal forests indicated that differences in the intensity disturbance could change the bryophytes species composition over time [106].

4.2. Indicator Species

Regarding the RIL effects on the species level, we identified D. sumichrastii as an indicator species for conserved sites and L. viticulosoides as an indicator species for sites under management. These mosses have been previously found to be associated with conserved and disturbed sites, respectively [107,108]. They have characteristics as stem size up to 6 cm and growing forms in tufts [60] that can facilitate their determination in the field and, therefore, their use as indicator species by non-specialists. However, since L. viticulosoides also appears as a common epiphytic moss in both forest types, it would be an adequate indicator only when occurring as terrestrial moss. At the same time, this species is common in open sites [58], which are not necessarily linked to management. In the case of D. sumichrastii, it is associated with shaded and partially open conditions [58], characteristics of common occurrence in conserved sites. Given its high percentage of specificity, it is an adequate indicator of this type of forest. For the epiphytic mosses, since the forest management did not significantly affect their communities, no species associated with any forest type were obtained. Thus, we recommend the use of L. viticulosoides and D. sumichrastii as ecological indicators of habitat conditions and forest quality monitoring in temperate forests of the Sierra Juarez, Oaxaca, Mexico.

4.3. Implications for Forest Management and Future Research

Our results supported the documented low effects of RIL on other taxa and determine that these forest practices with low levels of disturbances make it a suitable alternative to maintain high biodiversity [11,109,110]. To guaranty the biodiversity conservation of Sierra Juarez forests, it is crucial to consider the RIL harvesting areas compared to the surrounding conserved forest extensions, keeping the continuity of suitable conditions for mosses [26,111,112,113]. Thus, landscape features should be taken into account to explain the diversity levels in RIL harvesting areas in future studies [50,114]. Because our study is based on three sites under management by RIL, we recommend developing future research adding more study sites in order to increase the completeness of the moss inventory and obtain more fine results. In addition, to obtain a fine understanding of the effects of this silviculture strategy on biodiversity, as well as to study other taxonomical groups as fungi communities. We suggest studying other silvicultural methods in the region which could be an alternative to balance biodiversity conservation and the forest management developed by local communities (e.g., [115,116]).
Including ecologically important taxa such as mosses in the evaluation of RIL impacts will be useful to provide a more comprehensive assessment of the effects of this forest management strategy on the biodiversity which is of special relevance in areas with a high floristic diversity such as our study area.

5. Conclusions

This study supports the use of RIL as a forest management strategy able to maintain the diversity (alpha and beta) and structure of the moss community. This comprises 70 moss species including terrestrial and epiphytic taxa in the temperate forest in the Sierra Juarez, Oaxaca, Mexico. To conserve diversity in managed areas, it will be recommendable to develop strategies to preserve microenvironmental characteristics. For example, to continue with the conservation of Quercus trees in logging sites. Therefore, this research demonstrated the potential of D. sumichrasti and L. viticulosoides as indicator species of conserved and managed sites, respectively. In addition, because RIL is developed by local communities, these results support the positive outcomes in the management and conservation initiatives of their forest.
Finally, although the RIL system seems to help the conservation of a wide diversity of terrestrial and epiphytic mosses, it is recommended that further studies be developed to understand whether this type of management can also benefit, for example, the biomass and coverture of populations, as well as the development of processes such as water storage in which these plants participate. Similarly, the study of other biological groups understudied under RIL management may contribute to a better understanding of the effects of a forest management system that promises to be compatible with the conservation of biodiversity.

Author Contributions

Conceptualization: E.H.-R., E.M. Methodology: E.M., L.H.E.-V., E.H.-R., D.G.-Á. Samples taxonomic determination: E.H.-R., D.G.-Á. Data curation: E.H.-R., L.H.E.-V. Formal analysis: E.H.-R., E.M., L.H.E.-V. Interpretation results: E.H.-R., E.M., L.H.E.-V., M.M.G. Investigation: M.M.G., L.H.E.-V., E.H.-R. Project administration: E.H.-R., E.M. Resources: M.M.G. Supervision: E.M., M.M.G., L.H.E.-V. Validation: E.M., M.M.G., L.H.E.-V. Visualization: M.M.G., L.H.E.-V., E.H.-R. Writing—original draft: E.H.-R., E.M., L.H.E.-V., M.M.G. Writing—review & editing: E.H.-R., E.M., M.M.G., L.H.E.-V., D.G.-Á. Funding acquisition: E.H.-R., E.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Council of Science and Technology (CONACYT) (grant number 595991), The Sigma Xi Grant-in-Aid of Research (grant number G201603152071095), and PRODEP network: “Biodiversity conservation in anthropized environments” (No. Project 180790 of SEP-CONACYT).

Acknowledgments

We are grateful to different municipalities from Sierra Juárez for allow us to develop this work in their forests; to Bruce Allen to provide bibliographic material for the taxonomic determination of samples; to Jesus López Santiago and Violeta Sarai Jiménez Hernández for their support in the fieldwork. Also, we thank Nicole Fenton and Corrige-moi group of the Université du Québec en Abitibi-Témiscamingue for their valuable comments to this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Moss species checklist in conserved and managed areas with reduce-impact logging in the Sierra Juárez, Oaxaca. Species presence in each site is expressed with a T for terrestrial, E for the epiphytic (collected-on Quercus), and TE for those collected in both samplings. The acronyms are used in the rank-abundance curves.
Table A1. Moss species checklist in conserved and managed areas with reduce-impact logging in the Sierra Juárez, Oaxaca. Species presence in each site is expressed with a T for terrestrial, E for the epiphytic (collected-on Quercus), and TE for those collected in both samplings. The acronyms are used in the rank-abundance curves.
AcronymFamilies and SpeciesSites
C1C2C3M1M2M3
Amblystegiaceae
AnacomAnacamptodon compactus (Thér.) W.R. Buck E
CamhisCampylium hispidulum (Brid.) Mitt.TT
HygrobHygrohypnum robinsonii H.A. CrumE
Brachytheciaceae
AercapAerolindgia capillacea (Hornsch.) M. Menzel E
BraoxyBrachythecium oxycladon (Brid.) A. Jaeger T
BracirBrachythecium cirriphylloides K.D. McFarland T
BrarudBrachythecium ruderale (Brid.) W.R. Buck TE
Bryaceae
BrasysBrachymenium systylium (Müll. Hal.) A. Jaeger E
BrybilBryum billarderii Schwägr. TET T
RhohuiRhodobryum huillense (Welw. &Duby) Touw T
Calymperaceae
SyrincSyrrhopodon incompletus Schwägr. T
Daltoniaceae
AdebogAdelothecium bogotense (Hampe) Mitt. E
Dicranaceae
Dic1Dicranaceae 1 T
DicfrigDicranum frigidum Müll. Hal. T
DicscoDicranum scoparium Hedw. T
DicsumDicranum sumichrastii Duby TT
HolarbHolomitrium arboreum Mitt. T
HolpulHolomitrium pulchellum Mitt.T TEE
Entodontaceae
EntbeyEntodon beyrichii (Schwägr.) Müll. Hal.T
ErylonErythrodontium longisetum (Hook.) Paris ET
Fissidentaceae
FiscriFissidens crispus Mont.T
FiseleFissidens elegans Brid.T
Hedwigiaceae
BrasquBraunia squarrulosa (Hampe) Müll. Hal.TETEETETE
Hypnaceae
HercylHerzogiella cylindricarpa (Cardot) Z. Iwats.T TE
HomshaHomomallium sharpii Ando & Higuchi T
HypamaHypnum amabile (Mitt.) HampeT T
HypcupHypnum cupressiforme Hedw. T
MitrepMittenothamnium reptans (Hedw.) CardotTTTTETTE
PylfalPylaisia falcata Schimp.EE E
Lembophyllaceae
PilflePilotrichella flexilis (Hedw.) Ångstr. T
PilmauPilotrichella mauiensis (Sull.) A. Jaeger TE
Leskeaceae
HapangHaplocladium angustifolium (Hampe&Müll. Hal.) Broth. T
Leucodontaceae
LeucurLeucodon curvirostris HampeE EE
Meteoriaceae
MetillMeteorium illecebrum Sull.EE
Mniaceae
PlarosPlagiomnium rostratum (Schrad.) T.J. Kop. T
Neckeraceae
NecchlNeckera chlorocaulis Müll. Hal.ETETEEE
NecehrNeckera ehrenbergii Müll. Hal. E
PorlonPorotrichum longirostre (Hook.) Mitt. E
Orthotrichaceae
MactenMacrocoma tenuis (Hook. &Trev.) Vitt T T T
MacguaMacromitrium guatemalense Müll. Hal. E
MacshaMacromitrium sharpii H.A. Crum ex Vitt EEE
OrthorOrthotrichum hortoniae Vitt T
ZygehrZygodon ehrenbergii Müll. Hal. T E T
ZyglieZygodon liebmannii Schimp. E
ZygobtZygodon obtusifolius Hook. E
ZygreiZygodon reinwardtii (Hornsch.) A. Braun E
ZygvirZygodon viridissimus (Dicks.) Brid.EETEEE
Pilotrichaceae
TrasubTrachyxiphium subfalcatum (Hampe) W.R. Buck T
Polytrichaceae
PogsubPogonatum subflexuosum (Lorentz) Broth. T
PoljunPolytrichum juniperinum Hedw. T
Pottiaceae
HyoinvHyophila involute (Hook.) A. JaegerT
LepvitLeptodontium viticulosoides (P. Beauv.) Wijk & Margad. TTETTET
SynfraSyntrichia fragilis (Taylor) Ochyra E
TimanoTimmiella anomala (Bruch & Schimp.) Limpr. T
TribraTrichostomum brachydontium BruchT T
TritenTrichostomum tenuirostre (Hook. & Taylor) Lindb. E
Prionodontaceae
PridenPrionodon densus (Sw. ex Hedw.) Müll. Hal. E
PrilutPrionodon luteovirens (Taylor) Mitt. E TE
Pylaisiadelphaceae
HetaffHeterophyllium affine (Hook. in Kunth) Fleisch. T
Racopilaceae
RactomRacopilum tomentosum (Hedw.) Brid.T
Rhabdoweisiaceae
SymvagSymblepharis vaginata (Hook.) Wijk & Margad. T
Rhizogoniaceae
PyrspiPyrrhobryum spiniforme (Hedw.) Mitt. T
Rigodiaceae
RigtoxRigodium toxarion (Schwägr.) A. Jaeger E
Sematophyllaceae
AcrlonAcroporium longirostre (Brid.) W.R. Buck T
SemadnSematophyllum adnatum (Michx.) E. BrittonTTTE TE
SemswaSematophyllum swartzii (Schwägr.) W.H. Welch & H.A. CrumT TTTT
Thuidiaceae
CyrminCyrto-hypnum minutulum (Hedw.) W.R. Buck & H.A. CrumTTT
CyrschCyrto-hypnum schistocalyx (Müll. Hal.) W.R. Buck & H.A. Crum E
CyrshaCyrto-hypnum sharpii (H.A. Crum) W.R. Buck & H.A. CrumT
ThudelThuidium delicatulum (Hedw.) Schimp.TTTTTTE
Table A2. Correlation analysis values between alpha diversity and environmental variables. TMR = Terrestrial mosses richness; EMR = Epiphytic mosses richness; M_RTM = Mean richness of terrestrial mosses; M_REM = Mean richness of epiphytic mosses; SD_TMR = Standard deviation of the terrestrial mosses richness; SD_EMR = Standard deviation of epiphytic mosses richness; M_DT = Mean daily temperature; SD_DT = Standard deviation of the daily temperature; M_DL = Mean of the daily light; SD_DL = Standard deviation of the daily light; M_CC = Mean of the canopy coverage; SD_CC = Standard deviation of the canopy coverage; M_Alt = Mean of altitude; M_TAP = Mean of the total annual precipitation; M_Ori = Mean of the orientation; M_Slop = Mean of the slope. Bold numbers indicate correlations > 0.83.
Table A2. Correlation analysis values between alpha diversity and environmental variables. TMR = Terrestrial mosses richness; EMR = Epiphytic mosses richness; M_RTM = Mean richness of terrestrial mosses; M_REM = Mean richness of epiphytic mosses; SD_TMR = Standard deviation of the terrestrial mosses richness; SD_EMR = Standard deviation of epiphytic mosses richness; M_DT = Mean daily temperature; SD_DT = Standard deviation of the daily temperature; M_DL = Mean of the daily light; SD_DL = Standard deviation of the daily light; M_CC = Mean of the canopy coverage; SD_CC = Standard deviation of the canopy coverage; M_Alt = Mean of altitude; M_TAP = Mean of the total annual precipitation; M_Ori = Mean of the orientation; M_Slop = Mean of the slope. Bold numbers indicate correlations > 0.83.
TMREMRM_RTMM_REMSD_TMRSD_EMRM_DTSD_DTM_DLSD_DLM_CCM_AltM_TapM_OriM_SlopSD_CC
TMR1−0.090.80.210.21−0.090.59−0.12−0.38−0.350.57−0.51−0.62−0.020.61−0.63
EMR−0.0910.070.670.120.650.150.270.640.54−0.70.25−0.37−0.20.110.69
M_RTM0.80.0710.630.49−0.280.57−0.53−0.54−0.570.65−0.47−0.35−0.130.72−0.67
M_REM0.210.670.6310.210.080.49−0.210.110−0.05−0.07−0.370.050.450.05
SD_TMR0.210.120.490.211−0.07−0.11−0.81−0.58−0.580.34−0.250.44−0.870.46−0.32
SD_EMR−0.090.65−0.280.08−0.0710.320.460.770.83−0.69−0.2−0.36−0.330.250.71
M_DT0.590.150.570.49−0.110.32100.120.20.26−0.78−0.680.150.82−0.23
SD_DT−0.120.27−0.53−0.21−0.810.46010.80.78−0.670.39−0.590.54−0.470.62
M_DL−0.380.64−0.540.11−0.580.770.120.810.98−0.90.27−0.440.25−0.240.91
SD_DL−0.350.54−0.570−0.580.830.20.780.981−0.840.11−0.420.2−0.140.86
M_CC0.57−0.70.65−0.050.34−0.690.26−0.67−0.9−0.841−0.540.13−0.020.45−0.99
M_Alt−0.510.25−0.47−0.07−0.25−0.2−0.780.390.270.11−0.5410.180.27−0.90.48
M_Tap−0.62−0.37−0.35−0.370.44−0.36−0.68−0.59−0.44−0.420.130.181−0.45−0.29−0.08
M_Ori−0.02−0.2−0.130.05−0.87−0.330.150.540.250.2−0.020.27−0.451−0.39−0.01
M_Slop0.610.110.720.450.460.250.82−0.47−0.24−0.140.45−0.9−0.29−0.391−0.4
SD_CC−0.630.69−0.670.05−0.320.71−0.230.620.910.86−0.990.48−0.08−0.01−0.41

References

  1. Grulke, M.; Valle, P.D.; Calo, I.; Merger, E.; Pawlowski, G.; Wittmann, N. Sustainable Natural Forest Management in the Tropics: Best Practices and Investment Opportunities for Large-Scale Forestry; UNIQUE Forestry and Land Use: Freiburg, Germany, 2016. [Google Scholar]
  2. Cardinale, B.J.; Duffy, J.E.; Gonzalez, A.; Hooper, D.U.; Perrings, C.; Venail, P.; Narwani, A.; Mace, G.M.; Tilman, D.; Wardle, D.A.; et al. Biodiversity loss and its impact on humanity. Nature 2012, 486, 59–67. [Google Scholar] [CrossRef]
  3. Rands, M.R.W.; Adams, W.M.; Bennun, L.; Butchart, S.H.M.; Clements, A.; Coomes, D.; Entwistle, A.; Hodge, I.; Kapos, V.; Scharlemann, J.P.W.; et al. Biodiversity Conservation: Challenges Beyond 2010. Science 2010, 329, 1298–1303. [Google Scholar] [CrossRef] [Green Version]
  4. Foley, J.A.; DeFries, R.; Asner, G.P.; Barford, C.; Bonan, G.; Carpenter, S.R.; Chapin, F.S.; Coe, M.T.; Daily, G.C.; Gibbs, H.K.; et al. Global Consequences of Land Use. Science 2005, 309, 570–574. [Google Scholar] [CrossRef] [Green Version]
  5. Hasan, A.; Montoro Girona, M.; Grosbois, G.; Saha, N.; Halim, M.A. Land Sparing Can Maintain Bird Diversity in Northeastern Bangladesh. Sustainability 2020, 12, 6472. [Google Scholar] [CrossRef]
  6. Kim, S.; Axelsson, E.P.; Montoro Girona, M.; Senior, J.K. Continuous-cover forestry maintains soil fungal communities in Norway spruce dominated boreal forests. For. Ecol. Manag. 2021, 480. [Google Scholar] [CrossRef]
  7. Burivalova, Z.; Sekercioglu, C.H.; Koh, L.P. Thresholds of logging intensity to maintain tropical forest biodiversity. Curr. Biol. CB 2014, 24, 1893–1898. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Imai, N.; Samejima, H.; Langner, A.; Ong, R.C.; Kita, S.; Titin, J.; Chung, A.Y.; Lagan, P.; Lee, Y.F.; Kitayama, K. Co-benefits of sustainable forest management in biodiversity conservation and carbon sequestration. PLoS ONE 2009, 4, e8267. [Google Scholar] [CrossRef]
  9. Bawa, K.S.; Seidler, R. Natural Forest Management and Conservation of Biodiversity in Tropical Forests. Conserv. Biol. 1998, 12, 46–55. [Google Scholar] [CrossRef]
  10. Torres-Rojo, J.M.; Moreno-Sánchez, R.; Mendoza-Briseño, M.A. Sustainable Forest Management in Mexico. Curr. For. Rep. 2016, 2, 93–105. [Google Scholar] [CrossRef] [Green Version]
  11. Chaudhary, A.; Burivalova, Z.; Koh, L.P.; Hellweg, S. Impact of forest management on species richness: Global meta-analysis and economic trade-offs. Sci. Rep. 2016, 6, 23954. [Google Scholar] [CrossRef] [Green Version]
  12. Mongtagnini, F.; Jordan, C.F. Tropical Forest Ecology: The Basis for Conservation and Management; Springer: Berlin, Germany, 2005. [Google Scholar]
  13. Bicknell, J.E.; Struebig, M.J.; Edwards, D.P.; Davies, Z.G. Improved timber harvest techniques maintain biodiversity in tropical forests. Curr. Biol. 2014, 24, R1119–R1120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Günter, S.; Weber, M.; Stimm, B.; Mosandl, R. Silviculture in the Tropics; Springer: Berlin, Germany, 2011. [Google Scholar] [CrossRef]
  15. Hendrison, J. Damage-Controlled Logging in Managed Tropical Rain Forest in Suriname; Wageningen Agricultural University: Wageningen, The Netherlands, 1990. [Google Scholar]
  16. Putz, F.E.; Sist, P.; Fredericksen, T.; Dykstra, D. Reduced-impact logging: Challenges and opportunities. For. Ecol. Manag. 2008, 256, 1427–1433. [Google Scholar] [CrossRef]
  17. Putz, F.E.; Pinard, M.A. Reduced-impact logging as a carbon-offset method. Conserv. Biol. 1993, 7, 755–757. [Google Scholar] [CrossRef]
  18. Asia-Pacific Forestry Commission. Applying Reduced Impact Logging to Advance Sustainable Forest Management; Enters, T., Durst, P.B., Applegate, G.B., Kho, P.C.S., Man, G., Eds.; Food and Agriculture Organization of the United Nations and Regional Office for Asia and the Pacific: Bangkok, Thailand, 2002. [Google Scholar]
  19. Paillet, Y.; Berges, L.; Hjalten, J.; Odor, P.; Avon, C.; Bernhardt-Romermann, M.; Bijlsma, R.J.; De Bruyn, L.; Fuhr, M.; Grandin, U.; et al. Biodiversity differences between managed and unmanaged forests: Meta-analysis of species richness in Europe. Conserv. Biol. J. Soc. Conserv. Biol. 2010, 24, 101–112. [Google Scholar] [CrossRef]
  20. Ellis, E.A.; Montero, S.A.; Hernández Gómez, I.U.; Romero Montero, J.A.; Ellis, P.W.; Rodríguez-Ward, D.; Blanco Reyes, P.; Putz, F.E. Reduced-impact logging practices reduce forest disturbance and carbon emissions in community managed forests on the Yucatán Peninsula, Mexico. For. Ecol. Manag. 2019, 437, 396–410. [Google Scholar] [CrossRef]
  21. Álvarez, S.; Rubio, A. LÍnea Base De Carbono En Bosque Mixto De Pino-Encino De La Sierra JuÁrez (Oaxaca, MÉxico). AplicaciÓn Del Modelo Co2fix V.3.2. Rev. Chapingo Ser. Cienc. For. Ambiente 2013, XIX, 125–137. [Google Scholar] [CrossRef] [Green Version]
  22. Bray, D.B. The struggle for the forest: Conservation and Development in the Sierra Juarez. Grassroots Dev. 1991, 15, 12–25. [Google Scholar]
  23. Bray, D.B.; Pérez, L.M.; Barry, D. Los Bosques Comunitarios de México: Manejo Sustentable de Paisajes Forestales, 1st ed.; Secretaría de Medio Ambiente y Recursos Naturales, Instituto Nacional de Ecología, Instituto de Geografía de la Universidad Nacional Autónoma de México, Consejo Civil Mexicano para la Silvicultura Sostenible A.C., Florida International University: México Distrito Federal, México, 2007. [Google Scholar]
  24. Smith, D.M. The Practice of Silviculture, 8th ed.; Wiley: New York, NY, USA, 1986. [Google Scholar]
  25. Rainforest Alliance. Resumen Público de Certificación de Union de Productores Forestales Zapotecas-Chinantecas de la Sierra de Juárez de R.I. (UZACHI); Rainforest Alliance: New York, NY, USA, 1996; p. 48. [Google Scholar]
  26. Ódor, P.; Király, I.; Tinya, F.; Bortignon, F.; Nascimbene, J. Reprint of: Patterns and drivers of species composition of epiphytic bryophytes and lichens in managed temperate forests. For. Ecol. Manag. 2014, 321, 42–51. [Google Scholar] [CrossRef]
  27. Rainforest Alliance. Resumen Público de Certificación Comunidad Ixtlan de Juarez; Rainforest Alliance: New York, NY, USA, 2006; p. 34. [Google Scholar]
  28. Putz, F.E.; Zuidema, P.A.; Synnott, T.; Peña-Claros, M.; Pinard, M.A.; Sheil, D.; Vanclay, J.K.; Sist, P.; Gourlet-Fleury, S.; Griscom, B.; et al. Sustaining conservation values in selectively logged tropical forests: The attained and the attainable. Conserv. Lett. 2012, 5, 296–303. [Google Scholar] [CrossRef] [Green Version]
  29. Carreño-Rocabado, G.; Peña-Claros, M.; Bongers, F.; Alarcón, A.; Licona, J.-C.; Poorter, L.; Vesk, P. Effects of disturbance intensity on species and functional diversity in a tropical forest. J. Ecol. 2012, 100, 1453–1463. [Google Scholar] [CrossRef]
  30. Imai, N.; Seino, T.; Aiba, S.-I.; Takyu, M.; Titin, J.; Kitayama, K. Effects of selective logging on tree species diversity and composition of Bornean tropical rain forests at different spatial scales. Plant Ecol. 2012, 213, 1413–1424. [Google Scholar] [CrossRef] [Green Version]
  31. Webb, E.L.; Peralta, R. Tree community diversity of lowland swamp forest in Northeast Costa Rica, and changes associated with controlled selective logging. Biodivers. Conserv. 1998, 7, 565–583. [Google Scholar] [CrossRef]
  32. Chen, Y.; Niu, S.; Li, P.; Jia, H.; Wang, H.; Ye, Y.; Yuan, Z. Stand Structure and Substrate Diversity as Two Major Drivers for Bryophyte Distribution in a Temperate Montane Ecosystem. Front. Plant Sci. 2017, 8, 874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Vanderpoorten, A.; Goffinet, B. Introduction to Bryophytes; Cambridge University Press: Cambridge, UK, 2009. [Google Scholar]
  34. Fenton, N.J.; Hylander, K.; Pharo, E.J. Bryophytes in forest ecosystems. In Routledge Handbook of Forest Ecology; Peh, K.S.-H., Corlett, R.T., Bergeron, Y., Eds.; Routledge: London, UK, 2015; pp. 241–251. [Google Scholar]
  35. Söderström, L. Conservation Biology of Bryophytes. Lindbergia 2006, 31, 24–32. [Google Scholar]
  36. Frego, K.A. Bryophytes as potential indicators of forest integrity. For. Ecol. Manag. 2007, 242, 65–75. [Google Scholar] [CrossRef]
  37. Proctor, M.C.F. Physiological ecology. In Bryophyte Biology, 2nd ed.; Goffinet, B., Shaw, A.J., Eds.; Cambridge University Press: Cambridge, UK, 2009; pp. 237–268. [Google Scholar]
  38. Sheridan, R.P. Nitrogenase activity by Hapalosiphon flexuosus associated with Sphagnum erythrocalyx mats in the cloud forest on the volcano La Soufriere, Guadeloupe, French West Indies. Biotropica 1991, 23, 134–140. [Google Scholar] [CrossRef]
  39. Glime, J.M. Bryological Interaction. In Bryophyte Ecology; Michigan Technological University and the International Association of Bryologists: Houghton, MI, USA, 2017; Volume 2. [Google Scholar]
  40. Suárez-Mota, M.E.; Villaseñor, J.L.; Ramírez-Aguirre, M.B. Sitios prioritarios para la conservación de la riqueza florística y el endemismo de la Sierra Norte de Oaxaca, México. Acta Bot. Mex. 2018, 49–74. [Google Scholar] [CrossRef]
  41. Rodríguez-Rivera, V. Distribución, Riqueza y Composición Química de Agallas en dos Bosques Templados: Uso Potencial Como Indicadores Ecológicos; Universidad de la Sierra Juárez: Oaxaca, México, 2014. [Google Scholar]
  42. Dávila, P.; Torres, L.; Torres-Colín, R.; Herrera-McBryde, O. Sierra Juárez, Oaxaca, México. In Centres of Plant Diversity: A Guide and Strategy for their Conservation: The Americas; Davis, S., Heywood, D., Herrera-MacBryde, V.H., Villa-Lobos, O., Hamilton, A.C., Eds.; The World Wide Fund For Nature, International Union for Conservation of Nature, The World Conservation Union: Cambridge, UK, 1997; Volume 3, pp. 135–138. [Google Scholar]
  43. Alfaro-Sánchez, G. Suelos. In Biodiversidad de Oaxaca; García-Mendoza, A.J., Ordóñez-Díaz, M.J., Briones-Salas, M.A., Eds.; Instituto de Biología, Universidad Nacional Autónoma de México, Fondo Oaxaqueño para la Conservación de la Naturaleza, World Wildlife Fund: México Distrito Federal, México, 2004; pp. 55–65. [Google Scholar]
  44. Torres-Colín, R. Tipos de vegetación. In Biodiversidad de Oaxaca; García-Mendoza, A.J., Ordoñez, M.J., Briones-Salas, M., Eds.; Instituto de Biología, Universidad Nacional Autónoma de México, Fondo Oaxaqueño para la Conservación de la Naturaleza, World Wildlife Fund: México Distrito Federal, México, 2004; pp. 105–117. [Google Scholar]
  45. Gradstein, S.R.; Nadkarni, N.M.; Krömer, T.; Holz, I.; Nicole, N. A Protocol for Rapid and Representative Sampling of Vascular and Non-Vascular Epiphyte Diversity of Tropical Rain Forests. Selbyana 2003, 24, 105–111. [Google Scholar]
  46. Delgadillo, M.C. Vegetación de México; Rzedowski, J., Ed.; Limusa: México Distrito Federal, México, 1978. [Google Scholar]
  47. Benitez, H.; Arizmendi, C.; Marquez, L. Base de Datos de las Áreas de Importancia para la Conservación de las Aves. Available online: http://conabioweb.conabio.gob.mx/aicas/doctos/aicas.htm (accessed on 4 August 2018).
  48. Mittermeier, R.A.; Robles, G.P.; Hoffmann, M.; Pilgrim, J.; Brooks, T.; Goettsch Mittermeier, C.; Lamoreux, J.; Da Fonseca, G.A.B. Hotspots Revisited; Cemex, Conservation International and Agrupacion Sierra Madre: Mexico Distrito Federal, Mexico, 2004. [Google Scholar]
  49. Enquist, B.J.; Feng, X.; Boyle, B.; Maitner, B.; Newman, E.A.; Jørgensen, P.M.; Roehrdanz, P.R.; Thiers, B.M.; Burger, J.R.; Corlett, R.T.; et al. The commonness of rarity: Global and future distribution of rarity across land plants. Sci. Adv. 2019, 5, eaaz0414. [Google Scholar] [CrossRef] [Green Version]
  50. Pazos-Almada, B.; Bray, D.B. Community-based land sparing: Territorial land-use zoning and forest management in the Sierra Norte of Oaxaca, Mexico. Land Use Policy 2018, 78, 219–226. [Google Scholar] [CrossRef]
  51. Hernández-Rodríguez, E.; Escalera-Vázquez, L.; Calderón-Patrón, J.M.; Mendoza, E. Mamíferos medianos y grandes en sitios de tala de impacto reducido y de conservación en la sierra Juárez, Oaxaca. Rev. Mex. de Biodivers. 2019, 90. [Google Scholar] [CrossRef]
  52. Navarro Cerrillo, R.M.; Esteves Vieira, D.J.; Ochoa-Gaona, S.; de Jong, B.H.J.; del Mar Delgado Serrano, M. Land cover changes and fragmentation in mountain neotropical ecosystems of Oaxaca, Mexico under community forest management. J. For. Res. 2018, 30, 143–155. [Google Scholar] [CrossRef]
  53. Delgadillo, M.C. Briofitas. In Manual de Herbario. Administración y Manejo de Colecciones, Técnicas de Recolección y Preparación de Ejemplares Botánicos; Lot, A., Chiang, F., Eds.; México: México Distrito Federal, México, 1986; pp. 77–82. [Google Scholar]
  54. Allen, B.H. Moss Flora of Central America Part 3 Anomodontaceae-Symphyodontaceae; Missouri Botanical Garden: St. Louis, MO, USA, 2010; Volume 117. [Google Scholar]
  55. Allen, B.H. Moss Flora of Central America Part 2 Encalyptaceae-Orthotrichaceae; Missouri Botanical Garden: St. Louis, MO, USA, 2002; Volume 90. [Google Scholar]
  56. Allen, B.H. Moss Flora of Central America Part 1 Sphagnaceae-Calymperaceae; Missouri Botanical Garden: St. Louis, MO, USA, 1994; Volume 48. [Google Scholar]
  57. Buck, W.R. Pleurocarpous Mosses of West Indies; New York Botanical Garden: New York, NY, USA, 1998; Volume 82. [Google Scholar]
  58. Gradstein, S.R.; Churchill, S.P.; Salazar-Allen, N. Guide to the Bryophytes of Tropical America; New York Botanical Garden Press: New York, NY, USA, 2001; Volume 86. [Google Scholar]
  59. Ireland, R.R.; Buck, W.R. Some Latin American Genera of Hypnaceae (Musci); Smithsonian Institution Scholarly Press: Washington, DC, USA, 2009; Volume 93. [Google Scholar]
  60. Sharp, A.J.; Crum, H.; Eckel, P.M. The Moss Flora of Mexico; The New York Botanical Garden: New York, NY, USA, 1994; Volume 69. [Google Scholar]
  61. Delgadillo, M.C. LATMOSS 2010. Available online: http://www.ibiologia.unam.mx/briologia/ (accessed on 30 September 2018).
  62. INEGI. Conjunto de Datos Vectoriales de la Carta de usoUso del Suelo y Vegetación, Escala 1:250000, Serie V. Available online: http://www.inegi.org.mx/geo/contenidos/recnat/edafologia/vectorial_serieii.aspx (accessed on 14 March 2017).
  63. QGIS.org. QGIS Geographic Information System. QGIS Association. Available online: https://www.qgis.org/en/site/getinvolved/faq/index.html (accessed on 9 November 2020).
  64. Chao, A.; Gotelli, N.J.; Hsieh, T.C.; Sander, E.L.; Ma, K.H.; Colwell, R.K.; Ellison, A.M. Rarefaction and extrapolation with Hill numbers: A framework for sampling and estimation in species diversity studies. Ecol. Monogr. 2014, 84, 45–67. [Google Scholar] [CrossRef] [Green Version]
  65. Carvalho, J.C.; Cardoso, P.; Gomes, P. Determining the relative roles of species replacement and species richness differences in generating beta-diversity patterns. Glob. Ecol. Biogeogr. 2012, 21, 760–771. [Google Scholar] [CrossRef]
  66. Feinsinger, P. Designing Field Studies for Biodiversity Conservation; Island Press: Washington, DC, USA, 2001. [Google Scholar]
  67. Escalera-Vázquez, L.H.; Zambrano, L. The effect of seasonal variation in abiotic factors on fish community structure in temporary and permanent pools in a tropical wetland. Freshw. Biol. 2010, 55, 2557–2569. [Google Scholar] [CrossRef]
  68. Kobza, R.M.; Trexler, J.C.; Loftus, W.F.; Perry, S.A. Community structure of fishes inhabiting aquatic refuges in a threatened Karst wetland and its implications for ecosystem management. Biol. Conserv. 2004, 116, 153–165. [Google Scholar] [CrossRef]
  69. McCune, B.; Grace, J.B.; Urban, D.L. Analysis of Ecological Communities; MjM Software Design: Gleneden Beach, OR, USA, 2002. [Google Scholar]
  70. Dufrene, M.; Legendre, P. Species assemblages and indicator species: The need for a flexible asymmetrical approach. Ecol. Monogr. 1997, 67, 345–366. [Google Scholar] [CrossRef]
  71. Legendre, P.; Legendre, L. Numerical Ecology, 3rd ed.; Elsevier: Oxford, UK, 2012. [Google Scholar]
  72. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  73. Hsieh, T.C.; Ma, K.H.; Chao, A.; McInerny, G. iNEXT: An R package for rarefaction and extrapolation of species diversity (Hill numbers). Methods Ecol. Evol. 2016, 7, 1451–1456. [Google Scholar] [CrossRef]
  74. Cardoso, P.; Rigal, F.; Carvalho, J.C. BAT: Biodiversity Assessment Tools; R CRAN: Helsinki, Finland, 2017. [Google Scholar]
  75. 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 CRAN: Helsinki, Finland, 2018. [Google Scholar]
  76. Cáceres, M.D.; Legendre, P. Associations between species and groups of sites: Indices and statistical inference. Ecology 2009, 90, 3566–3574. [Google Scholar] [CrossRef]
  77. Peterson, B.G.; Carl, P. PerformanceAnalytics: Econometric Tools for Performance and Risk Analysis; R CRAN: Evanston, IL, USA, 2018. [Google Scholar]
  78. Dexter, E.; Rollwagen-Bollens, G.; Bollens, S.M. The trouble with stress: A flexible method for the evaluation of nonmetric multidimensional scaling. Limnol. Oceanogr. Methods 2018, 16, 434–443. [Google Scholar] [CrossRef]
  79. Montoro Girona, M.; Rossi, S.; Lussier, J.M.; Walsh, D.; Morin, H. Understanding tree growth responses after partial cuttings: A new approach. PLoS ONE 2017, 12, e0172653. [Google Scholar] [CrossRef]
  80. Montoro Girona, M.; Lussier, J.M.; Morin, H.; Thiffault, N. Conifer Regeneration After Experimental Shelterwood and Seed-Tree Treatments in Boreal Forests: Finding Silvicultural Alternatives. Front. Plant Sci. 2018, 9, 1145. [Google Scholar] [CrossRef]
  81. Montoro Girona, M.; Morin, H.; Lussier, J.-M.; Ruel, J.-C. Post-cutting Mortality Following Experimental Silvicultural Treatments in Unmanaged Boreal Forest Stands. Front. For. Glob. Chang. 2019, 2. [Google Scholar] [CrossRef] [Green Version]
  82. D’Amato, A.W.; Bradford, J.B.; Fraver, S.; Palik, B.J. Forest management for mitigation and adaptation to climate change: Insights from long-term silviculture experiments. For. Ecol. Manag. 2011, 262, 803–816. [Google Scholar] [CrossRef]
  83. Bose, A.K.; Harvey, B.D.; Brais, S.; Beaudet, M.; Leduc, A. Constraints to partial cutting in the boreal forest of Canada in the context of natural disturbance-based management: A review. Forestry 2013, 87, 11–28. [Google Scholar] [CrossRef]
  84. Rosenvald, R.; Lõhmus, A. For what, when, and where is green-tree retention better than clear-cutting? A review of the biodiversity aspects. For. Ecol. Manag. 2008, 255, 1–15. [Google Scholar] [CrossRef]
  85. Castro-Arellano, I.; Presley, S.J.; Saldanha, L.N.; Willig, M.R.; Wunderle, J.M. Effects of reduced impact logging on bat biodiversity in terra firme forest of lowland Amazonia. Biol. Conserv. 2007, 138, 269–285. [Google Scholar] [CrossRef]
  86. Johansson, T.; Hjältén, J.; de Jong, J.; von Stedingk, H. Environmental considerations from legislation and certification in managed forest stands: A review of their importance for biodiversity. For. Ecol. Manag. 2013, 303, 98–112. [Google Scholar] [CrossRef]
  87. Fenton, N.J.; Imbeau, L.; Work, T.; Jacobs, J.; Bescond, H.; Drapeau, P.; Bergeron, Y. Lessons learned from 12 years of ecological research on partial cuts in black spruce forests of northwestern Québec. For. Chron. 2013, 89, 350–359. [Google Scholar] [CrossRef] [Green Version]
  88. Elek, Z.; Kovacs, B.; Aszalos, R.; Boros, G.; Samu, F.; Tinya, F.; Odor, P. Taxon-specific responses to different forestry treatments in a temperate forest. Sci. Rep. 2018, 8, 16990. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Yu, D.; Zhou, L.; Zhou, W.; Ding, H.; Wang, Q.; Wang, Y.; Wu, X.; Dai, L. Forest management in Northeast China: History, problems, and challenges. Environ. Manag. 2011, 48, 1122–1135. [Google Scholar] [CrossRef] [PubMed]
  90. Åström, M.; Dynesius, M.; Hylander, K.; Nilsson, C. Effects of slash harvest on bryophytes and vascular plants in southern boreal forest clear-cuts. J. Appl. Ecol. 2005, 42, 1194–1202. [Google Scholar] [CrossRef]
  91. Ódor, P.; Standovár, T. Richness of bryophyte vegetation in near-natural and managed beech stands: The effects of management-induced differences in dead wood. Ecol. Bull. 2001, 49, 219–229. [Google Scholar]
  92. Vellak, K.; Ingerpuu, N. Management Effects on Bryophytes in Estonian Forests. Biodivers. Conserv. 2005, 14, 3255–3263. [Google Scholar] [CrossRef]
  93. Tullus, T.; Rosenvald, R.; Leis, M.; Lõhmus, P. Impacts of shelterwood logging on forest bryoflora: Distinct assemblages with richness comparable to mature forests. For. Ecol. Manag. 2018, 411, 67–74. [Google Scholar] [CrossRef]
  94. Fenton, N.J.; Frego, K.A.; Sims, M.R. Changes in forest floor bryophyte (moss and liverwort) communities 4 years after forest harvest. Can. J. Bot. 2003, 81, 714–731. [Google Scholar] [CrossRef]
  95. Oldén, A.; Ovaskainen, O.; Kotiaho, J.S.; Laaka-Lindberg, S.; Halme, P. Bryophyte species richness on retention aspens recovers in time but community structure does not. PLoS ONE 2014, 9, e93786. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Romero, C. Reduced-impact logging effects on commercial non-vascular pendant epiphyte biomass in a tropical montane forest in Costa Rica. For. Ecol. Manag. 1999, 118, 117–125. [Google Scholar] [CrossRef]
  97. Sporn, S.G.; Bos, M.M.; Hoffstätter-Müncheberg, M.; Kessler, M.; Gradstein, S.R. Microclimate determines community composition but not richness of epiphytic understory bryophytes of rainforest and cacao agroforests in Indonesia. Funct. Plant Biol. 2009, 36, 171–179. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  98. Löbs, N.; Walter, D.; Barbosa, C.G.G.; Brill, S.; Alves, R.P.; Cerqueira, G.R.; de Oliveira Sá, M.; de Araújo, A.C.; de Oliveira, L.R.; Ditas, F.; et al. Microclimatic conditions and water content fluctuations experienced by epiphytic bryophytes in an Amazonian rain forest. Biogeosciences 2020, 17, 5399–5416. [Google Scholar] [CrossRef]
  99. Söderström, L. The occurrence of epixylic bryophyte and lichen species in an old natural and a managed forest stand in Northeast Sweden. Biol. Conserv. 1988, 45, 169–178. [Google Scholar] [CrossRef]
  100. Longton, R.E.; Hedderson, T.A. What are rare species and why conserve them? Lindbergia 2000, 25, 53–61. [Google Scholar]
  101. Gaston, K.J. Rarity; Springer Science and Business Media: Comwall, UK, 1994. [Google Scholar] [CrossRef] [Green Version]
  102. Márialigeti, S.; Németh, B.; Tinya, F.; Ódor, P. The effects of stand structure on ground-floor bryophyte assemblages in temperate mixed forests. Biodivers. Conserv. 2009, 18, 2223–2241. [Google Scholar] [CrossRef]
  103. Goffinet, B.; Shaw, A.J. Bryophyte Biology, 2nd ed.; Cambridge University Press: Cambridge, UK, 2009. [Google Scholar] [CrossRef]
  104. Táborská, M.; Kovács, B.; Németh, C.; Ódor, P.; Wulf, M. The relationship between epixylic bryophyte communities and microclimate. J. Veg. Sci. 2020, 31, 1168–1180. [Google Scholar] [CrossRef]
  105. Boudreault, C.; Paquette, M.; Fenton, N.J.; Pothier, D.; Bergeron, Y. Changes in bryophytes assemblages along a chronosequence in eastern boreal forest of Quebec. Can. J. For. Res. 2018, 48, 821–834. [Google Scholar] [CrossRef]
  106. Fenton, N.J.; Bergeron, Y. Stochastic processes dominate during boreal bryophyte community assembly. Ecology 2013, 94, 1993–2006. [Google Scholar] [CrossRef]
  107. Rovere, A.E.; Calabrese, G.M. Diversidad de musgos en ambientes degradados sujetos a restauración en el Parque Nacional Lago Puelo (Chubut, Argentina). Rev. Chil. Hist. Nat. 2011, 84, 571–580. [Google Scholar] [CrossRef]
  108. Zander, R. Genera of the Pottiaceae: Mosses of Harsh Environments; The Buffalo Society of Natural Sciences: New York, NY, USA, 1993. [Google Scholar]
  109. Melo, F.P.L.; Arroyo-Rodríguez, V.; Fahrig, L.; Martínez-Ramos, M.; Tabarelli, M. On the hope for biodiversity-friendly tropical landscapes. Trends Ecol. Evol. 2013, 28, 462–468. [Google Scholar] [CrossRef]
  110. Arroyo-Rodriguez, V.; Fahrig, L.; Tabarelli, M.; Watling, J.I.; Tischendorf, L.; Benchimol, M.; Cazetta, E.; Faria, D.; Leal, I.R.; Melo, F.P.L.; et al. Designing optimal human-modified landscapes for forest biodiversity conservation. Ecol. Lett. 2020, 23, 1404–1420. [Google Scholar] [CrossRef]
  111. Paltto, H.; Nordén, B.; Götmark, F. Partial cutting as a conservation alternative for oak (Quercus spp.) forest-Response of bryophytes and lichens on dead wood. For. Ecol. Manag. 2008, 256, 536–547. [Google Scholar] [CrossRef]
  112. Vanderpoorten, A.; Engels, P. Patterns of bryophyte diversity and rarity at a regional scale. Biodivers. Conserv. 2003, 12, 545–553. [Google Scholar] [CrossRef]
  113. Benítez, Á.; Prieto, M.; Aragón, G. Large trees and dense canopies: Key factors for maintaining high epiphytic diversity on trunk bases (bryophytes and lichens) in tropical montane forests. Forestry 2015, 88, 521–527. [Google Scholar] [CrossRef] [Green Version]
  114. Robson, J.P.; Berkes, F. Exploring some of the myths of land use change: Can rural to urban migration drive declines in biodiversity? Glob. Environ. Chang. 2011, 21, 844–854. [Google Scholar] [CrossRef]
  115. Bautista, L.J.; Damon, A.; Ochoa-Gaona, S.; Tapia, R.C. Impact of silvicultural methods on vascular epiphytes (ferns, bromeliads and orchids) in a temperate forest in Oaxaca, Mexico. For. Ecol. Manag. 2014, 329, 10–20. [Google Scholar] [CrossRef]
  116. Ramírez Santiago, R.; Ángeles Pérez, G.; Hernández de La Rosa, P.; Cetina Alcalá, V.M.; Plascencia Escalante, O.; Clark-Tapia, R. Efectos del aprovechamiento forestal en la estructura, diversidad y dinámica de rodales mixtos en la Sierra Juárez de Oaxaca, México. Madera y Bosques 2019, 25. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Study sites in Sierra Juárez, Oaxaca (black circles). Conserved sites: Santa María Yavesía (C1), Santa Catarina Lachatao (C2), and San Juan Luvina (C3). Managed sites: Ixtlán de Juárez (M1), San Juan Evangelista Analco (M2), and Capulálpam de Méndez (M3). Stitch line is the road 175 Oaxaca-Tuxtepec. (b) Examples of conservation (C1) and (c) management (M2) sites.
Figure 1. (a) Study sites in Sierra Juárez, Oaxaca (black circles). Conserved sites: Santa María Yavesía (C1), Santa Catarina Lachatao (C2), and San Juan Luvina (C3). Managed sites: Ixtlán de Juárez (M1), San Juan Evangelista Analco (M2), and Capulálpam de Méndez (M3). Stitch line is the road 175 Oaxaca-Tuxtepec. (b) Examples of conservation (C1) and (c) management (M2) sites.
Forests 12 00383 g001
Figure 2. Rarefaction curves (R/E curves) with intrapolation (solid line)/extrapolation (dotted line) of the species richness of (a) terrestrial and (b) epiphytic mosses among the study sites (Guides). Each curve has confidence intervals of 95%. Confidence intervals overlapping indicate no significant differences in the species richness observed and estimated between managed and conservation sites.
Figure 2. Rarefaction curves (R/E curves) with intrapolation (solid line)/extrapolation (dotted line) of the species richness of (a) terrestrial and (b) epiphytic mosses among the study sites (Guides). Each curve has confidence intervals of 95%. Confidence intervals overlapping indicate no significant differences in the species richness observed and estimated between managed and conservation sites.
Forests 12 00383 g002
Figure 3. Rank-abundance curves of terrestrial (a) and epiphytic (b) mosses for each site, in managed and conservation forests. The acronyms correspond to the names of the species in Appendix A Table A1.
Figure 3. Rank-abundance curves of terrestrial (a) and epiphytic (b) mosses for each site, in managed and conservation forests. The acronyms correspond to the names of the species in Appendix A Table A1.
Forests 12 00383 g003
Figure 4. Non-metric multidimensional scaling (NMDS) of the terrestrial (a) and epiphytic (b) mosses in managed sites (hollow black circles and dotted ellipse) and conserved (solid black circles and continuous ellipse).
Figure 4. Non-metric multidimensional scaling (NMDS) of the terrestrial (a) and epiphytic (b) mosses in managed sites (hollow black circles and dotted ellipse) and conserved (solid black circles and continuous ellipse).
Forests 12 00383 g004
Table 1. Observed (Obs.) and estimated (Est.) richness for terrestrial and epiphytic mosses among the study sites (C = conserved; M = managed) by Chao 2 index. CI = Confidence intervals; % spp. record = percentage of species recorded.
Table 1. Observed (Obs.) and estimated (Est.) richness for terrestrial and epiphytic mosses among the study sites (C = conserved; M = managed) by Chao 2 index. CI = Confidence intervals; % spp. record = percentage of species recorded.
SiteTerrestrialEpiphytic
Obs.Est. (CI) (Chao 2)% spp. RecordObs.Est. (CI) (Chao 2)% spp. Record
C11739 (25–78)45718 (10–55)39
C21928 (22–44)681140 (18–123)28
C31749 (28–115)351215 (13–21)80
M11013 (11–20)77913 (10–23)69
M21140 (18–123)281147 (21–150)23
M31428 (18–62)50912 (10–19)75
Table 2. Beta diversity of terrestrial and epiphytic mosses between conserved and managed forests. Species replacement (βrepl), differences in richness (βrich), and total beta diversity (βTot) are expressed in percentage.
Table 2. Beta diversity of terrestrial and epiphytic mosses between conserved and managed forests. Species replacement (βrepl), differences in richness (βrich), and total beta diversity (βTot) are expressed in percentage.
TerrestrialEpiphytic
SitesβreplβrichβTot βrepl βrich βTot
C1-C273780313162
C1-C379079503181
C2-C36477180585
M1-M263669531367
M1-M353217488088
M3-M26714819010100
C-M52247674074
Table 3. Local and regional variables comparison between managed and conserved sites using T-student or Wilcoxon test.
Table 3. Local and regional variables comparison between managed and conserved sites using T-student or Wilcoxon test.
Variablest or Wilcoxon Valuep
Local-scale variables (LSV)Canopy coverage (%)W = 70.35
Altitude (m a.s.l)t = −1.50.22
Daily temperature (°C)t = 2.60.07
Daily light (lx)t = −0.50.66
Regional-scale variables (RSV)Total annual rainfall (mm)t = −1.80.16
Orientation (°)W = 70.40
Slope (°)t = 1.40.23
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hernández-Rodríguez, E.; Escalera-Vázquez, L.H.; García-Ávila, D.; Montoro Girona, M.; Mendoza, E. Reduced-Impact Logging Maintain High Moss Diversity in Temperate Forests. Forests 2021, 12, 383. https://doi.org/10.3390/f12040383

AMA Style

Hernández-Rodríguez E, Escalera-Vázquez LH, García-Ávila D, Montoro Girona M, Mendoza E. Reduced-Impact Logging Maintain High Moss Diversity in Temperate Forests. Forests. 2021; 12(4):383. https://doi.org/10.3390/f12040383

Chicago/Turabian Style

Hernández-Rodríguez, Enrique, Luis H. Escalera-Vázquez, Deneb García-Ávila, Miguel Montoro Girona, and Eduardo Mendoza. 2021. "Reduced-Impact Logging Maintain High Moss Diversity in Temperate Forests" Forests 12, no. 4: 383. https://doi.org/10.3390/f12040383

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