1. Introduction
The study of the dynamics of forest systems has recently undergone a resurgence of interest, especially in the Mediterranean region where these systems show a tumultuous history of human action and climate change [1]. The distribution areas of European Mediterranean forests have increased, related to the abandonment of agricultural land in recent decades, which has led to the appearance of secondary forests [2]. By contrast, in North Africa, forest degradation has reached alarming levels and has led to the disappearance of endemic and local species and the regeneration of drought-tolerant species [3].
The Atlas cedar (Cedrus atlantica Manetti) is an endemic species of the highest mountains in North Africa (Algeria and Morocco), where it occupies a very fragmented area [4]. This fragmentation is essentially linked to past climatic variations and anthropogenic action [5]. Palynological studies show that the natural range of this species before and during the early Holocene was more extensive than at present [6]. This species is better suited to humid and sub-humid bioclimates and suffers during long and repeated droughts [7].
The Atlas cedar has an important place in the Mediterranean landscape in Algeria and Morocco [8]. Unfortunately, the status of Atlas cedar stands in Algeria indicates a strong degradation linked primarily to anthropogenic action and the unfavorable climatic conditions marked by persistent drought since the 1980s [9][10]. This situation has led to a loss of vitality and the emergence of tree decline [11][12]. Additionally, a deficiency of regeneration has been reported for several decades, attributed mainly to anthropogenic action (e.g. grazing) and natural factors―such as erratic seed production, fruit consumption and climatic variability [8][13]. Because of the high longevity and large size of C. atlantica adult trees, their populations remain as dominant components of the forest biomass and architecture for long periods, although functioning demographically as remnant populations which lack significant constant regeneration [8][14]. The absence of appropriate forest management practices appears to have exacerbated this situation [15]. In this context, the study of the structure and dynamics of the Atlas cedar stands is indispensable for their future management. Research into tree age coupled with analysis of stand structures, land-use history, climatic data and the ecological attributes of different tree species has proven to be a suitable approach for understanding long-term forest dynamics. In Algeria, most studies on the structure of C. atlantica forests have focused almost exclusively on how spatial distribution has changed in relation to extrinsic factors, specifically climate or disturbances such as grazing or fire [16][17][18]. Additionally, C. atlantica has been widely used in dendrochronology because of its sensitivity to climatic stress and its high longevity [9][10][18][19]. However, no studies on the structural dynamics of C. atlantica in relation to dendrochronological dates, addressing temporal changes, have been performed in this country.
Here, we describe and analyze the forest structure, regeneration and recruitment history of C. atlantica stands in the “Theniet El Had” Atlas cedar forest (northwestern Algeria), in pure and mixed forests located along an altitudinal gradient, as a necessary step prior to the development of forest management alternatives to promote forest conservation at the regional scale. The results will improve our understanding of the structure and dynamics of C. atlantica forests and allow a management plan to be proposed based on a comprehensive understanding of the structures and regeneration processes.
2. Material and Methods
2.1. Study Area
The study was conducted in the National Park of Theniet El Had (hereafter TNP), located in the Eastern extension of the Ouarsenis massif (northwestern Algeria) (coordinates 35˚47'N and 35˚54'N, 01˚54'E and 02˚02'E) (Figure 1). The TNP is an area of 3460 ha. It was the first space to be protected in Algeria by the French colonial government, in 1923. After independence, it was re-created as a National Park by the Algerian government, in 1983. Cedar stands cover an area of 1000 ha, and they are located mainly on the North slope [20]. Since the 19th century, uncontrolled use of the natural forests (fire, grazing and illegal logging) has led to their gradual degradation [15]. The terrain is a NW-SE slope, with the existence of sandstone formations, reflecting extremely rough topography, especially in the Northern exposures where the slope can exceed 50%. The altitude
Figure 1. Distribution map of Cedrus atlantica in the National Park of Theniet El Had (North-West Algeria) and the distribution of the sites studied.
varies between 854 ma.s.l. and 1786 ma.s.l. The geology is characterized by the presence of Oligocene sediments composed of a thick series of sandstones superimposed on thin layers of clay and marl [21]. The dominant soils in the C. atlantica forests are shallow and are sandstone or marl with pH = 7.3 - 8.6, C/N ratio = 4 - 10 and low contents of organic matter ([C]= 13 to 81 kg・ha−1), N (3 - 8.1 kg・ha−1), phosphorus ([P2O5]= 5 - 71 kg・ha−1) and potassium ([K2O]= 0.35 - 1.8 kg・ha−1).
In the TNP, C. atlantica forests occur naturally between 1200 and 1786 ma.s.l. but are rare below 1300 m due to land-use pressure. The stands included in this study are dominated by C. atlantica, often mixed with Quercus ilex L. subsp. ballota (Desf.) Samp and Q. faginea Lam., and, to a lesser extent, Pinus halepensis Mill., Q. suber L. and Juniperus oxycedrus L.
2.2. Climatic Context
The climate is characterized by cool, wet winters and hot, dry summers. Annual precipitation varies between 580 mm (854 ma.s.l.) and 868 mm (1786 ma.s.l.), with only 5% falling between June and August, resulting in late-summer drought. The drier period lasts 3 - 4 months, from May to September, giving the study area a Mediterranean climate [22], with the presence of three variants depending on an altitudinal gradient: semi-arid (<1200 ma.s.l.), sub-humid (between 1200 and 1500 ma.s.l.) and humid (>1500 ma.s.l.). Snow cover usually lasts 20 - 30 days between November and March. A trend towards a drier climate has occurred over the past two decades (1980-2002), with a 25% decrease in annual precipitation [11]. The minimum temperature of the coldest month (January) varies between −2.2˚C and −0.6˚C, and the maximum temperature of the hottest month (July) varies between 27.9˚C and 32.7˚C.
2.3. Experimental Design
Atlas cedar forests in the TNP are spatially heterogeneous and can be classified along an altitudinal gradient. High-altitude stands are characterized by their pure cedar composition and absence of recent human intervention. In medium-altitude forests, the species composition is more complex and human impacts are evident. Although the extent and frequency of logging in these forests are unclear, it is likely that a few moderate, to large-sized trees, have been harvested for local use.
A long this environmental gradient, at each altitude, care was taken to select a plot representative of the site and forest conditions according to the following criteria: 1) at least 10 adult trees were present, 2) the forest patch in which the stand was located was dominated by cedar, 3) reliable information about when the stand was last logged was available, 4) there was no sign of recent disturbance (e.g., fire, windstorm) other than occasional livestock grazing, and 5) it was accessible [14]. Additionally, the selected sites were mapped, and their geographical centers, areas and size were determined (Table 1). In May 2013, six 60 ×
Table 1. Study sites, environmental features and core samples for six Cedrus atlantica populations in the National Park of Theniet El Had (North-West Algeria).
50 m plots were established in C. atlantica stands according to a stratified random design. Forests were classified as a) High altitude and North-facing exposure with pure cedar stands (sites 1, 2 and 3), b) Medium altitude with mixed cedar stands (sites 4, 5 and 6).
In each site, all living and dead trees (diameter at breast height, dbh > 10 cm) were identified, and dbh and height (H) were measured using, respectively, a caliper and a hypsometer (n = 423). Regeneration was classified in four size classes: seedlings (h < 50 cm), short saplings (50 cm ≤ h < 130 cm), tall saplings (h ≥ 130 cm and dbh < 5 cm) and juveniles (5 cm ≤ dbh < 10 cm) [14]. For age determination, two cores at breast height per cedar tree with dbh ≥ 5.0 cm (n = 169) were extracted using a Suunto® increment borer. We added 12 years to the number of years counted at the sampling height to obtain the estimated years of establishment [23]. The age distribution is presented in 10-year intervals.
2.4. Population Structure Analysis
We grouped the C. atlantica trees into 10-cm-diameter size classes, averaged densities across sampling plots and generated diameter frequency distributions for each species at each site. To compare the population size structures quantitatively, the Weibull distribution was fitted to each of the empirical dbh probability distributions obtained from the censuses. We fitted two-parameter Weibull functions [24]using the maximum likelihood method.
Prior to the statistical analysis we examined the dbh and stand density for normality and homoscedasticity using the Kolmogorov-Smirnov test and Levene’s test, respectively. The square-root transformation of dbh was used to stabilize the variances. Once the basic requirements were found to be met, the stand density of adult trees, dbh, height, stand basal area, and regeneration of C. atlantica, Q. ilex, and Q. faginea were compared between pure and mixed stands using analysis of variance (one-way ANOVA) [25]. Null hypotheses were rejected at the P < 0.05 level and all analyses, and the results in tables are shown as means with their standard errors for the untransformed variables. All statistical analyses were performed using R, version 3.4.0 [26].
2.5. Radial Growth Dynamics
From 15 dominant Atlas cedar individuals per plot, two cores at breast height were taken using a Suunto® increment borer, in a direction parallel to the contour lines, to study radial growth dynamics (Table 1). The cores were mounted, sanded and polished; the crossdating was performed under a binocular microscope, following standard dendrochronological methods (Speer 2010). Visual crossdating of samples was conducted before measuring and it was verified using COFECHA after measuring [27]. The ring width was measured to the nearest 0.01 mm using a LINTAB measuring device and the Time Series Analysis Program (TSAP). A 10-year filtering window was chosen to enhance high-frequency signal and to fit a function to the measured chronologies using the 3pBase program [28]. Two chronologies were created and the trend due to the geometrical constraint of adding a volume of wood to a stem of increasing radius was corrected by converting tree-ring widths into basal area increments (BAI) using the formula
where R is the radius of the tree and t is the year of tree-ring formation using the 3pBase program [28].
Cumulative growth curves for each series were established to study the effect of environmental conditions on the productivity of each site [29].
3. Results
3.1. Stand Structure
Among the adult trees, all six plots were dominated by C. atlantica (Table 2). The density of C. atlantica individuals (dbh ≥ 5 cm) ranged from 94 (29% of total tree density) in mixed stands to 203 (76%) trees ha-1in pure stands and was significantly higher in pure cedar stands than in mixed stands (F = 7.79, P = 0.049). The densities of Q. faginea ranged between 81 (25%) and 16 (6%) trees ha−1 and Q. ilex ranged between 21 (6%) and 43 (16%) trees ha−1, without significant differences between forest types (F = 6.63, P = 0.082; F = 1.21, P = 0.332, respectively).
The height of C. atlantica differed significantly (F = 21.81; P < 0.01) between forest types, where the average value in pure stands was 14 m and the average value in mixed stands was 11 m. The tree height of Q. faginea was comparable but for Q. ilex and J. oxycedrus were smaller than that of cedar; however, these differences were not significant between forest types (F = 0.04; P = 0.84; F = 1.63; P = 2.70, respectively).
The mean dbh of C. atlantica ranged from 37.5 cm (mixed stands) to 50.0 cm (pure stands) (Table 2), differing significantly among sites (F = 27.48, P < 0.01). The mean diameters of Q. faginea and Q. ilex ranged from 41.6 cm (mixed
Table 2. Structural characteristics of the sampling sites in Theniet El Had National Park. Values are means ± SE. Statistical differences between forest types based on one-way ANOVA are indicated (ns = non-significant, * = P < 0.05. ** = P < 0.01).
stands) to 36.5 (pure stands) cm and from 28.6 (mixed stands) to 16.7 (pure stands) cm, respectively without significant differences. The basal area of C. atlantica ranged from 14.3 m2・ha−1 (mixed stands) to 45.7 m2・ha−1 (pure stands, F = 9.39, P = 0.03), while those of Q. faginea and Q. ilex ranged from 9.2 m2・ha−1 (mixed stands) to 1.9 m2・ha−1 (pure stands, F = 0.77, P = 0.44) and from 1.6 m2・ha−1 (mixed stands) to 1.2 m2・ha−1 (pure stands, F = 0.22, P = 0.66), respectively. In pure stands, C. atlantica accounted for over 93% of the basal area, but in mixed stands less than 57%. For J. oxycedrus, P. halepensis and Q. suber in mixed stands, the mean basal areas were1, 3.6 and 2.9 m2・ha−1, respectively.
The distribution of the diameter classes of C. atlantica (Figure 2) shows the dominance of trees in the middle classes (40 - 50 cm) in pure stands and in low classes (10 - 20 cm) in mixed stands type. The dbh size-class Weibull distribution was positively-skewed in form in mixed stands (c = 1.3) as result of a peak in abundance for the smallest and intermediate size classes (<30 - 40 cm), and a rapid decline in the number of individuals in larger size classes (Figure 2), which indicates that there were abundant young individuals in the stand and that the regeneration status of C. atlantica was better in this forest. The Weibull shape statistic indicates that the distributions in pure stands (c = 2.0) followed a negative-skewed distribution, meaning few young individuals compared to the individuals of average age and a poor regeneration status. The distribution of the diameter classes of Q. faginea and Q. ilex showed a high abundance of individuals in the smaller size classes and of intermediate-sized trees (20 < dbh < 40 cm) (data not included). For the other species, the dbh distribution shows the dominance of smaller size classes for J. oxycedrus and of intermediate-sized trees for Q. suber and P. halepensis (data not included).
Seedlings, short saplings and tall saplings of C. atlantica were very low for both stand types, ranging from 8.8 trees ha−1 (mixed stands) to 9.9 trees ha−1 (pure stands, F = 0.014, P = 0.931), and they were completely absent in some of the plots (data non-included). The density of juveniles ranged from 3.3 (mixed stands) to 0 trees ha−1 (pure stands, F = 0.60, P = 0.49) for Q. faginea, and from
Figure 2. Diameter-class histograms for Cedrus atlantica stands in the National Park of Theniet El Had (North-West Algeria). The x-axis is uniform and covers 10-cm diameter intervals. Sample size (N), mean diameter and standard deviation (x), and Weibull shape statistics (b, c) are also included.
5.5 (mixed stands) to 9.0 trees ha−1 (pure stands, F = 0.108, P = 0.75) for Q. ilex.
3.2. Cedrus atlantica Colonization and Demography
The age data from all the samples (N = 90) and the recruitment history at the study site suggest that stands of Atlas cedar are old with a clear dominance three age classes: 100 - 120 years (44.7%), 120 - 140 years and 140 - 160 years with 21.1% respectively (Figure 3). The oldest tree in the stand was 172 years old, dating from 1838, but most of the aged trees entered the population during the 20th century; only two of the trees cored were more than 160 years old. The percentage of trees less than 100 years old does not exceed 10.5%. The populations exhibited establishment peaks in the 1860-1880s and 1890-1900s (Figure 4).
3.3. Dynamics of Radial Growth
The curves representing the temporal variations in BAI show the existence of favorable and unfavorable phases for growth (Figure 5). These phases are synchronous for both stands types, suggesting that their origin is common and can be attributed to climatic factors. The stands show a divergence in growth from
Figure 3. Frequency of C. atlantica trees in the National Park of Theniet El Had (North-West Algeria) according to age classes.
Figure 4. Cedrus atlantica individuals established and still surviving over the past 200 years. The tree recruitment peaks correspond to 1860, 1880 and 1900, after the significant mortality phases recorded at the stand.
1910 (Figure 5), the productivity of the pure stand is higher than that of the mixed stand. This situation can be explained by illegal logging in mixed stands that are at low altitude and closer to riparian populations. The land surface curves show a tendency to reduce the land area in both stands during the last three decades; this result can be associated with adverse climatic conditions (drought) and the senescence of stands.
4. Discussion
4.1. Stand Structure
The current structure of C. atlantica stands results from the combined effects of
Figure 5. Cumulative basal area increment (BAI) curves of pure and mixed Cedrus atlantica stands in the National Park of Theniet El Had (North-West Algeria).
site characteristics (altitude, aspect, slope and soil), climate and human pressure. Consequently, we can distinguish two types of Atlas cedar stand in the TNP: pure and mixed. The pure C. atlantica stands occupy the high altitudes (>1500 m a.l.s), where the climatic conditions are more favorable for the species and human activity is moderate. These stands are characterized by high density (≥186.6 trees ha−1) and low values of basal area (≥32.4 m2・ha−1) which indicate competition between trees. In moderate altitude (<1500 m a.s.l), C. atlantica grows in association with Q. ilex and, to a lesser degree, with J. oxycedrus, P. halepensis and Q. suber where human activity has been very intense (grazing, fire, opening of forest areas after illicit logging) and the climatic conditions are less favorable, and in association with Q. faginea in wet and deep soil [4]. The values of the density (≤160 trees ha−1) and basal area (≤25.8 m2・ha−1) were low. These results are like those found in other cedar forests of Algeria and Morocco [14][30]. We find the absence of Taxus baccata and Ilex aquifolium in this cedar forest, while the presence of Acer monspessullanum is sporadic. This bears witness to the more-xeric conditions, in comparison to other northern cedar forests [31].
4.2. Age Structure and Recruitment
The C. atlantica stands are jeopardized by the lack of regeneration and high tree senescence. Juveniles are absent in some sites and their density remains very low in others (3.3 to 26.6 trees ha−1). The reconstruction of tree recruitment in the TNP has shown peaks linked to climate. Indeed, young cedar recruitment has followed phases of adult tree mortality, especially in the 1860s, 1880s and 1890s, due to drought [32]. This mortality likely created gaps in the stands, allowing young cedars to grow after the establishment of favorable climatic conditions. In this context, [33]noted the existence of an enormous amount of dead wood in this forest dating back to the beginning of French colonization (1840) following an anarchic exploitation but from the 1880s the deaths due to drought have amplified the situation.
The role of the climate and the forestry practices during this period are at the origin of the current situation of cedar forest. In 1923, the cedar forest of Theniet El Had was proclaimed a National Park with a total ban on cutting and pasture. This situation would have favored regeneration at first time, but after the closure of the stands and severe drought conditions (during the 1940s) the regeneration process was deteriorated. After the independence of Algeria (1962), the situation did not improve, and the forest continued to undergo all kinds of aggression: illegal cutting, excessive grazing even after its proclamation as a national park in 1983 [20]. Currently, the Atlas cedar forests are the most-important socio-economic factor in the region. Pastoralism, which remains the main activity of the local population, compromises forest regeneration [20]. Additionally, persistent drought since the 1980s [10][11]may have had a contributing role in the lack of regeneration.
Regeneration problems have been mentioned for all stands of cedar throughout its range in Algeria [13][15]and Morocco [12][14][34]. These authors suggested a set of causes to explain this phenomenon: stand senescence, overgrazing, compaction and soil degradation, drought and the lack of appropriate management.
4.3. Dynamics of Radial Growth
The radial growth and the productivity of C. atlantica seem to be better in low-density sites with a Northern exposure, high altitude and little slope. This species is very sensitive to fluctuations in climate, especially rainfall [18][35]. The decrease in radial growth since the 1980s reflects the impact of the drought that has affected the region [19]. This decline was associated with high tree mortality, including the older trees [11][36]. Future climate scenarios indicate an increase in temperature of 3˚C by 2100 and a decrease in rainfall of 25% in North Africa [37], which may compromise the existence of this species in the area.
4.4. Proposal for the Management of Cedar Forests and Suggestions for New Research
The current state of C. atlantica stands requires the development of a management plan for the short- and medium-term, to improve natural regeneration and reduce tree mortality. This plan must consider the particularities of the environmental conditions of cedar forests and should integrate as much as possible the local populations. Algerian forest policy prescribes a regulatory framework for the management of national parks. Indeed, the TNP is sub-divided into several zones depending on the interests of conservation and the socio-economic equilibrium. The cedar forest is mainly found in the integral zone and benefits from a major conservation interest [20]. The application of this policy faces a complex socio-economic situation because local populations use the forest for livestock grazing and to obtain different non-timber products (firewood, resin, etc.). These activities must be moved outside the integral zone to prevent the degradation of the already-fragile cedar forest.
Natural regeneration can be improved through an integrated silvicultural plan. First of all, it will be necessary to apply a thinning program in some stands to reduce density [38]. Once the regeneration is established, seedlings must be protected by reducing grazing [35]. The use of artificial regeneration could be considered in some areas where degradation has reached advanced stages (declining or senescent stands, degraded soils). Direct seeding has already proven successful [8].
5. Conclusion
Site history (fire, forest pests and overgrazing) and lack of management have significantly reduced the extent of old-growth C. atlantica forests in the Theniet El Had National Park during the past two decades. The results of this study show a clear relationship between structures and dynamics of C. atlantica with climate-related mortality events suggesting that structural characteristics of those old-growth woodlands may be threatened under current and projected future conditions. The lack of regeneration has led a reduction in cedar cover and structural complexity. Our analysis suggests that these changes in structure and dynamic together with mortality events may alter the future of C. atlantica forests. However, our results are limited by sampling limits as number of locations, recruitment data, or the inclusion of other perturbation events as fires and forest pests. Therefore, more research is still needed concerning the ecological dynamic and biotic and abiotic impacts on C. atlantica forests to improve sustainable ecosystem management. In particular, those related to pest and climate change impacts on cedar forest in Argelia may be useful in providing data for ecological studies of the relations between pest outbreaks, stand structure, species diversity, site history (fire, logging, grazing, etc.), and climatic events. We suggest silvicultural and conservative management strategies oriented to maintain and preserve not just the overstory canopies, but also the understory complexity and diversity that supports seedling recruitment, so that the woodland may be sustainable in future climates.
Acknowledgements
This work was funded through the PNR project (Algerian Ministry of Higher Education and Scientific Research). The authors thank the Algerian Ministry of Agriculture and Directorate General of Forests for the permission to work in the National Park of Theniet El Had. We are very grateful to Abdelkader Masloub (National Park of Theniet El Had, Algeria) for his valuable assistance during fieldwork. We also acknowledge the institutional support of the University of Cordoba-Campus de Excelencia CEIA3 and ESPECTRAMED (CGL2017-86161-R) project.