Next Article in Journal
Composition and Spatial Variation of Germinable Seed Bank in Burned Nothofagus pumilio Forests in Patagonia Argentina
Next Article in Special Issue
Disturbance-Dependent Variation in Functional Redundancy Drives the Species Versus Functional Diversity Relationship across Spatial Scales and Vegetation Layers
Previous Article in Journal
Radiocarbon Investigation of the Historic African Baobabs of Omusati, Namibia
Previous Article in Special Issue
Quality Assessment and Rehabilitation of Mountain Forest in the Chongli Winter Olympic Games Area, China
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Forests
Volume 13
Issue 11
10.3390/f13111900
forests-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Ecological Potential of Mediterranean Habitats for Black Pine (Pinus nigra J.F. Arnold) in Croatia

by
Damir Barčić
1,
Vladimir Hršak
2,
Roman Rosavec
1,* and
Mario Ančić
1
1
Faculty of Forestry and Wood Technology, University of Zagreb, 10000 Zagreb, Croatia
2
Faculty of Science, University of Zagreb, 10000 Zagreb, Croatia
*
Author to whom correspondence should be addressed.
Forests 2022, 13(11), 1900; https://doi.org/10.3390/f13111900
Submission received: 17 October 2022 / Revised: 6 November 2022 / Accepted: 8 November 2022 / Published: 11 November 2022
(This article belongs to the Special Issue Ecological Forestry and Restoration)
Browse Figures
Review Reports Versions Notes

Abstract

:
This study outlines research performed on experimental plots established in forest plantations and in natural black pine forests. The objective of the study was to determine the habitat factors that have the greatest impact on the growth and development of natural black pine forests and black pine forests plantations upon the return of climatozonal vegetation. Plots 625 m2 in size were selected and vegetation inventories made, and the structural elements of black pine plantations examined. Multivariate analysis found that elevation and slope were the two variables that had the highest positive correlation with the floristic composition of the vegetation. Furthermore, one of the goals was to determine the differences with regard to reforestation with black pine. The analysis of the isolated experimental plots showed a clear grouping of plots according to habitat factors. The differences in the synecological factors in the research area resulted in the presence of sub-Mediterranean, epi-Mediterranean and eu-Mediterranean elements within the black pine forest plantations. This also reduces the risk of the occurrence and development of large wildfires.

1. Introduction

In terms of amelioration, the role of black pine throughout the Mediterranean karst area is exceptionally valuable. Black pine is perhaps the most important species for karst reforestation due to the various ecological and social roles it plays, in addition to its commercial importance. It also plays an irreplaceable role in preventing habitat degradation and the prevention erosion from expanding. This has contributed to its widespread distribution, from the western to the eastern Mediterranean. According to Quézel and Médail [1], black pine is represented by six subspecies in the Mediterranean region. Black pine is the most widely distributed pine species in mountain areas of the Mediterranean Basin and commonly used for afforestation [2]. These pine forests are classified as “habitats of European interest” and require specific conservation measures according to the Convention for the Conservation of European Wildlife and Natural Habitats (Resolution 4/1996), due to the lack of successful natural regeneration [3].
Black pine has been used as a reforestation species throughout virtually the entire area the northern Adriatic coast in Croatia. Particular success was achieved during the work and establishment of the Inspectorate (i.e., Inspectorate for the forestation of karst, bare soils and flood management in the northern Adriatic; [4]). During the 20th century, black pine was used as a pioneer species to stop erosion processes, degradation and devastation, not only in the Dinaric karst area, but also in the Apennine Mountains [5,6,7]. Though the ecological and production potential of black pine is prominent in the Dinaric karst, this sub-Mediterranean area is largely threatened with the risk of wildfires [8,9]. Nevertheless, increasing land abandonment in Mediterranean Europe during the last decades has favoured forest recovery [10,11] including Pinus nigra J.F. Arnold. stands [12]. Meanwhile throughout its European distribution range [13,14], black pine has important economic value due to the high quality of the wood products, and also for its strong ecological and protective roles (i.e., black pine forests are included in the EU Habitats Directive under no. 9530*(Sub-Mediterranean pine forests). In sub-Mediterranean habitats, black pine has certain limitations, above all, periods of drought and water deficits. In the western Mediterranean, this is a key factor in the regeneration of forest ecosystems [15]. Projected decreases in precipitation and increases in temperature in future climate scenarios may aggravate their vulnerability and alter their current distribution in the Mediterranean Basin [3]. In addition to the high vulnerability of sub-Mediterranean areas to global climate change [16] and changes in the seasonal distribution of precipitation over the past decade [17], it is evident that these changes have had negative effects on forests and their stability.
Black pine was one of the fundamental reforestation species during the great forest planting project (starting 150 years ago) in Croatian forestry [18]. In Croatia, the distribution, phytocenology of black pine forests and taxonomic research have been examined by many authors [19,20,21]. According to the Flora Europaea [22], it is possible to differentiate two subspecies of black pine in Croatia, ssp. nigra and ssp. dalmatica. Vidaković [23] listed three taxa of black pine in Croatia: Pinus nigra ssp. austriaca, Pinus nigra ssp. dalmatica and Pinus nigra ssp. gočensis var. illyrica. However, according to Liber [21], Vidaković’s report of the taxon ssp. gočensis is questionable, as this is likely just an ecotype characteristic for serpentine substrates.
The habitat conditions of the study area where the black pine stands are found in the western part of the Croatian Dinaric karst are best characterised by their heterogeneity. Erosion is also an important driver of land degradation, in that case afforestation effort in troubled areas is an expensive practice; the success of afforestation is not guaranteed, particulary in degraded soils [24]. This area is encompassed by the Mediterranean phytogeographic area, and is also under the strong influence of the Euro-Siberian–North American phytogeographic area. Moreover, the diversity of environmental conditions, geological and geomorphological factors and edaphic conditions can be associated with a higher degree of vulnerability to various erosion processes. Defining those relationships in karst, and determining the role of black pine stands in those processes, is therefore an important aspect of researching forest amelioration. A complicating circumstance for the management of these forests is the occurrence of forest fires, and subsequent regeneration after fires. Therefore, Santana [9] warns of the key issues to address after fire: when and how to extract the wood mass? The reason for this is the need to reduce all risks of erosion processes. Simultaneously, Mediterranean habitats are under very strong and prolonged drought in the recent period because plants of the Mediterranean region are exposed to drought stress [25,26].
This paper compares and analyses the possibilities of the habitat and potential of black pine stands under varying conditions. Namely, habitat conditions (slope, shallow soil, fire) largely govern the management of black pine in the Dinaric karst in Croatia. It should also be noted that this area also falls within the NATURA 2000 ecological network, where management recommendations require finding a balance between production objectives and conservation needs [27]. Therefore, one of the goals in that study suggested the need for appropriate silviculture works aimed at increasing the quality of black pine stands, and for their future management. In addition, in our study one of the objectives is to find possible differences between black pine plantations (established by reforestation) and natural pine stand. Furthermore, an analysis of the habitat conditions and production potential indicate the possible fulfilment of the ecological (transition towards mixed deciduous stands) and the economic roles (quality of wood products, biomass, etc.).

2. Materials and Methods

2.1. Study Area

Experimental plots were set up in black pine stands in western and southwestern Croatia (northern Adriatic region). Plot selection was based on data from the Forestry Management Plan, and an overview of the geological, pedological and phytocenological maps, taking into account certain factors such as stand age, human impacts, microrelief, elevation, slope and insolation. All plots were set up in stands over 40 years old. Of the 50 experimental plots, the majority (44) were set up in forest plantations, while 6 plots were set up in natural black pine stands (Figure 1 and Figure 2). Study area is represented by the Cfsax” climate according to Köppen. It is a moderately warm rainy climate with hot summers and an average monthly temperature above 22 °C. In the structure of the land cover, skeleton-rich soils of shallow Eutric Cambisol on carbonate rocks, marl and flysch, red soil (terra rossa) and rendzina on carbonate rocks predominate. Surface stoniness and rockiness often contribute to low surface productivity ratings.
The bioclimatically explored area (Table 1) is humid and subhumid with moderate cold climate. The maximum precipitation is in late fall or early winter. North and northeast winds dominate.

2.2. Data Collection

Stand structure analysis was performed on the experimental plots (each 625 m2 in size). For each plot, black pine trees were measured by number of trees, basal area and wood volume. Wood volume was calculated according to the volume tables and yield according to the yield-tables [28]. Vegetation relevés were performed on each experimental plot using the phytosociology methodology [29,30]. The abundance and cover of species was estimated using the expanded scale according to Barkman et al. [31]. Assessments were transformed on an ordinal scale according to van der Maarel [32], prior to analysis. Plant nomenclature was determined according to Nikolić [33,34,35]. Field studies on experimental plots also included the sampling of the humus layer formed from five individual samples within each experimental plot. Samples were taken for the determination of soil reaction, humus and total nitrogen content, and the C/N ratio. Depending on the degree of acidity or pH value of soil, the dominant pedogenic processes in the soil can be determined. To calculate the quantity of forest floor, three samples were taken along the diagonal of a 25 × 25 cm experimental plot. The procedure was performed in accordance with the UNECO guidebook [36].

2.3. Data Analyses

Laboratory analyses were performed on samples of the humus layer:
  • Soil reaction was measured electrometrically with a combined electrode in a suspension of soil in water, i.e., in 0.01 M CaCl2 in the ratio of 1:2.5 for surface mineral and argyl-accumulative layers, and in the ratio of 1:10 for prominently humic surface horizons. A laboratory microprocessor pH meter (MA 5736, Metrel, accuracy ± 0.01 pH) was used for the measurements.
  • Humus content was determined using the bichromatic method according to Tjurin.
  • The total nitrogen content was determined by incineration according to the Kjeldajl procedure and distillation according to Bremner.
  • Samples of the forest floor were ground and homogenised and incinerated in a microwave oven with HNO3 (Milestone Laboratory System). K content was determined by flame-photometry, and all other elements determined using the AAS method (atomic absorption spectrophotometry).

2.3.1. Statistical Analysis

First, cluster analysis was performed using the Soerensen (Bray-Curtis) distance as resemblance measure with the group average as clustering method [37]. The obtained results are arbitrarily divided into 5 clusters, which are further analyzed separately. Indicator Species Analysis (ISA) was used to determine species that are relevant for particular clusters obtained by cluster analysis [38]. Indicative types were those with p < 0.05.

2.3.2. Multivariate Fuzzy Set Ordination

The analysis was done according to [39,40] using PCOrd 7.01 software. As this analysis is an important issue and order of environmental variables, an initial preliminary analysis was made to determine the number and order of the variables involved in the analysis (such as slope, altitude, organic matter, total nitrogen). The analysis was first made for the entire set of data and then each cluster was analyzed individually to determine the importance of individual variables for the vegetation structure. As a resemblance measure in MSFO, the Soerensen (Bray–Curtis) distance was used.

3. Results

According to the cluster analysis, the obtained optimal number of cluster was five. Figure 3 showed the first cluster included 29 plots, the second 7 plots, the third 5 plots, the fourth 5 plots and fifth 4 plots. In the first cluster, favorable habitat conditions positively influenced the presence of sub-Mediterranean vegetation, while the emphasis here is also on the economic potential of certain forest cultures (P22, P27, P28, P29, P31, P32, P33). Wood volume in these plots ranged from 210.56 m3/ha to 411.44 m3/ha. Ecologically, this plot group fell within the sub-Mediterranean and epi-Mediterranean zone, at the border area of the montane vegetation belt that is dominated by common beech.
It is significant that the second cluster included virtually all island plots (island of Rab, Cres, Krk) in the warmer part of the sub-Mediterranean zone (P5, P6, P14, P15, P16, P17, P18). According to the phytocenological relevés, this is a transition area between the eu-Mediterranean and sub-Mediterranean zones. In the economic sense, the state of the forest cultures is quite poor. Wood volume deviates from the prescribed management plans, and ranges from the minimum 53.12 m3/ha to the maximum 282.88 m3/ha. The third cluster covers plots along the seashore at lower altitudes. Plots on the flysch substrate are also grouped here (P39, P40, P41).
The fourth cluster included plots in natural black pine stands (P48, P49, P50) and forest plantations in a protected area (P12 and P13). The joint characteristics of these plots is the greater slope and higher elevation in the natural black pine (Euphorbio triflorae-Pinetum nigrae Trinajstić 1999) stand. In the ecological sense, this is the sub-Mediterranean and epi-Mediterranean vegetation zone. In terms of economic potential, given the habitat quality class (IV), the condition is poor. The fifth cluster included three plots in a natural stand of the association Cotoneastro-Pinetum nigrae Ht. 1938. (P2, P3, P4), and one plot (P19) in the forest culture in the peak zone of the island of Krk. The common ecological characteristic of this cluster is the exceptionally strong impact of the dominant cold northeasterly wind (Bura). This circumstance impacts the poor economic potential of black pine in this area.
The importance of individual gradations is presented with Cum-R2 and Inc.-R2 (Table 2). These values are highest for altitude and slope. The nutrient gradient has a significantly lower value. The R2-Max value shows the association of environmental gradients with the vegetation structure. For altitudes and organic matter, it is smaller than the Cum-R2, which indicates good connectivity. The R2-Max slope gradient is larger than Cum-R2, indicating a weak connection. All this is confirmed by the biplot (Figure 4). The x-axis biplex represents the height gradient and the y-gradient represent nutritient gradient.
Indicator Species Analysis (ISA) shows indicator species for individual clusters (IV = indicator value, p = probability; A, B, C = layers in stand structure). The basis for the ISA was vegetation relevés made for each researched experimental plot.
Indicator species for cluster 1: Acer campestreC (IV = 51.7, p = 0.005), Clematis vitalbaB (IV = 38.8, p = 0.033), Cornus mas (IV = 44.6, p = 0.014).
Indicator species for cluster 2: Asparagus acutifolius (IV = 41.1, p = 0.017), Brachypodium sylvaticum (IV = 38.3, p = 0.04), Dactylis hispanica (IV = 44.3, p = 0.014), Phyllirea latifolia (IV = 28.6, p = 0.032), Quercus ilexC (IV = 42.9, p = 0.012), Rubia peregrina (IV = 28.6, p = 0.023), Rubus ulmifoliusB (IV = 50.0, 0.003), Salvia verbenacea (IV = 28.6, p = 0.034), Sonchus glaucescens (IV = 42.9, p = 0.014), Torilis japonica (IV = 42.9, p = 0.013).
Indicator species for cluster 3: Acer monspessulanumC (IV = 36.5, p = 0.026), Carpinus orientalisB (IV = 72.6, p = 0.002), Centaurea montana (IV = 50.0, p = 0.007), Chrysopogon gryllus (IV = 40.0, p = 0.023), Clematis flammula (IV = 47.9, p = 0.009), Cotinus cogyggria (IV = 49.6, p = 0.008), Festuca ovina (IV = 40.0, p = 0.023), Laurus nobilis (IV = 40.0, p = 0.025).
Ononis spinosa (IV = 40.0, p = 0.023), Quercus pubescensA (IV = 55.2, p = 0.001), Quercus pubescensC (IV = 36.2, 0.014), Rhamnus fallax (IV = 38.9, p = 0.023), Sorbus torminalisB (52.5, p = 0.006), Viburnum lantana (IV = 36.2, p = 0.034).
Indicator species for cluster 4: Acer obtusatumA2 (IV = 37.4, p = 0.028), Amelanchier ovalis (IV = 69.7, p = 0.001), Erica carnea (IV = 80.0, p = 0.001), Hepatica nobilis (IV = 36.8, p = 0.046), Molinia littoralis (IV = 40.0, p = 0.013), Ostrya carpinifoliaA (IV = 45.3, p = 0.015), Ostrya carpinifoliaB (IV = 43.3, p = 0.009), Phyteuma orbiculare (IV = 40.0, p = 0.013), Sesleria caerulea (IV = 40.0, 0.019), Sesleria tenuifolia (IV = 60.0, p = 0.001), Sorbus ariaA2 (IV = 52.3, p = 0.006).
Indicator species for cluster 5: Bromus erectus (IV = 49.1, p = 0.006), Bromus scoparius (IV = 28.6, p = 0.044), Chamaecytisus supinus (IV = 50.0, p = 0.005), Euphorbia seguieriana (IV = 50.0, p = 0.007), Galium sylvaticum (IV = 75.0, p = 0.002), Ligusticum mutellina (IV = 75.0, p = 0.002), Lonicera alpina (IV = 46.8, p = 0.005), Lonicera nigra (IV = 54.4, p = 0.007), Melica ciliata (IV = 59.7, p = 0.003), Moehringia muscosa (IV = 71.4, p = 0.002), Pinus nigraA (IV = 23.6, p = 0.002), Primula columnae (IV = 40.7, p = 0.01), Sorbus aucupariaB (IV = 40.7, 0.029), Thalictrum minus (IV = 61.0, p = 0.003), Valeriana tripteris (IV = 50.0, p = 0.005).
Figure 5 shows that the structure of the pine forest is irregular and has a diameter from 25.5 to 38.2 species per relevé. These data are related to different habitat conditions that determine the arrival of climatogenic vegetation.

4. Discussion

Each habitat factor directly affects all the ecological factors upon which plant life is dependent, while the feedback effect of plants on the habitat is just as important. The influence of plants is most evident in the case of pioneer stands, such as forest plantations. Habitat factors such as climate (air temperature, humidity, precipitation, wind), soil (soil type, mechanical and chemical composition, structure), relief (elevation, slope, exposure) and biotic influences are strongly associated through the interrelationships with ecological factors (light, heat, water, chemical and mechanical factors) [41,42]. The Mediterranean climate has two growing seasons, one in spring and the other in autumn, separated by a dry and hot summer period [43]. In recent decades we have combined effect of warming and drought [44]. The impact of climate change on vulnerability of Mediterranean-type biodiversity hotspots is likely to be highly dependent on the buffering effect of topographical complexity, and ergo, the existence and effectiveness of local refugia and microrefugia [45,46].
In the sub-Mediterranean habitats of Croatia, black pine primarily grows on limestone substrates (smaller ratio on dolomite), though it is also found on flysch substrates. It is most widely distributed in the sub-Mediterranean vegetation zone, where it is the most important species for reforestation of karst. This is also seen in other Mediterranean countries [47,48,49]. The habitat factors for black pine and other coniferous species in the Mediterranean, according to Quezel [50] and M’Hirit [51], were observed through the prism of the Emberger pluviometric coefficient. For black pine, it is significant that it appears in the humid, sub-humid and semi-arid zones.
In that sense, Trinajstić [19] stated that raising black pine plantations has phytogeographic justification only in the border areas between pubescent oak and beech areas, and that the closer these plantations to potential beech areas, the better the results will be. It was confirmed that black pine plantations, with regard to their composition and representation in the climazonal deciduous vegetation, is the best in the areas of pubescent oak forest, towards the borders of beech forests. From the aspect of amelioration, this is where these plantations can have the most direct and positive influence on the habitat conditions. Furthermore, in variable habitat conditions, as only forest plantations and natural stands older than 40 years were examined, there were no records, or only limited and individual findings of sprouting and saplings of black pine in a small number of plots. In the statistical analyses that divided the experimental plots into five main clusters, the largest cluster (29 plots) were situated in such vegetation areas as described above. Taking into account the productivity and economic role, in this cluster there are seven plots (P22, P27, P28, P29, P31, P32, P33) at elevations from 440 m (P22) to 714 m (P27) that stand out in particular. The edaphic conditions in these plots are favourable, both for the regeneration of climazonal vegetation, and for possible raising of new black pine plantations. The assessments depend on the manner and objectives of forest management. For example, silvicultural treatments, normally used for ecological or economic purposes, could have positive effect on biodiversity and biological interactions [48].
The natural stands of black pine (Euphorbio triflorae-Pinetum nigrae Trinajstić 1999) and (Cotoneastro-Pinetum nigrae Ht. 1938) are grouped into the fourth and fifth cluster, and are all at higher elevations (990 to 1082 m). The common ecological characteristics of these plots is the exceptionally strong influence of the dominant, cold, north-easterly Bura wind. In such habitat conditions, black pine has no competition from other species in the tree layer. An analysis of species indicators showed black pine as the indicator species only in the fifth cluster. In the other clusters, most species were characteristic for the sub-Mediterranean zone.
The edaphic factor is certainly decisive in the return of climazonal vegetation, particularly the limited soil horizon, i.e., forest floor. Our results showed (Figure 4 and Table 2) that nutrient gradient has a significantly lower value. We believe that this is to be expected, because it is mainly due to relatively shallow soils, and pines cannot have a positive amelioration effect during one rotation period (80 years). These two gradients of reliefs have the greatest impact on the structure of vegetation in the footage. The feed gradient has a much smaller effect as seen by the length of the vector. This means that the nutrient accumulation, ie the soil melioration effect on the total sample, may have a very limited effect and is generally not significant. It can only be effective in some places that can be determined by analyzing individual clusters.
In the study area, black pine stands are in the landscape with active erosion processes associated with soil degradation and the reduction in valuable land, as in other Mediterranean karst areas [52]. According to Martinović [18], the forest floor is a natural organic fertilizer containing all the nutritional elements required by forest tree species. The forest floor is important for regulating the hydrological conditions of the soil, reducing evaporation and increasing infiltration of precipitation, while it also positively affects temperature relations within the soil. It is closely associated with the general condition of humus in the soil, and its composition and properties depend on numerous factors and processes. One of the elementary processes is humus formation.
The influence of black pine is also considered in the relationship of soil types and productivity [18]. This leads to the production order of soils and specific indicators of soil fertility. There is stronger humus formation in black pine plantation in the northern Adriatic area than under natural forest vegetation, and that the surface layer of soil is richer in total nitrogen and physiologically active potassium. From the research results, and also in relation to the study by Trinajstić [19], in the phytogeographic sense of the productivity of black pine stands, it can be concluded that black pine plantations have their highest productivity on organic mineral black soils (higher sub-Mediterranean belt), followed by brown soil on limestone, and then on rendzina (lower sub-Mediterranean belt). The black pine forests in the Mediterranean have equal commercial and ecological importance, and therefore have been included in the list of Europe’s threatened habitats [27]. However, a problem arises in the lack of successful regeneration of those forests [14], which is also seen in the karst areas of Croatia. As some author mentioned natural regeneration is a main phase of forest dynamics and the only way to guarantee forest persistence in the long term [53]. In the relationship between soil and vegetation on karst, it must be considered how Mediterranean forest ecosystems cannot be considered for their production function alone; however, if these are commercial forests, they must meet the objectives of forest management. In the study area, the advantage is that most of the black pine forest plantation (included the cluster with 29 plots) are from pioneer stands, and near the end of their rotation period, shift in the direction of stand conversion. The process of progressive succession can be seen, with the regeneration of the sub-Mediterranean climazonal deciduous vegetation. Similar results can be seen in Tesei et al. [54], where the structure of the pine forest is irregular and has an average of 29.5 species per relevé (similar to our research, Figure 5). The dominant tree layer is monospecific as Pinus nigra, with an average height of 20 m and mean coverage of 70% to 80% and the shrub layer was well represented. The conversion process of forest cultures and conifer plantations is considered by many authors to be both significant and necessary [55,56,57]. Afforestation efforts to enhance adaptation of these forests in future climate scenarios will then require the selection of more drought-resistant genotypes and species and transformation to mixed, multi-aged forest stands [58]. This process may accelerate silvicultural works and in earlier thinning years to increase the quality of wood in the stand. Unfortunately, these works were not conducted in the study area during the prescribed 80-year harvest cycle (rotation period). In line with the forest management plan, as the document that governs the management of the forests and forest lands, 2000–2500 samplings of black pine per hectare are envisaged for raising a forest plantation. The first thinning cycles are planned between 20 and 30 years, and already at that time the number of trees should be reduced to 1300 to 1500 trees per hectare. Thinning should be continued every 10 years until the end of the rotation period (80 years). In comparing this with the research results, it is evident that there are deviations of wood volume on most of the experimental plots. The main reason for this is the lack of silvicultural works in the stands (cleaning and thinning). These works were only conducted in a small number of forest plantations as weak thinning campaigns and sanitary cutting due to fires in open spaces. Similar works in black pine plantations was also recorded by Zlatanov et al. [59]. Given that these results indicate the underperformance of wood volume, different strategies are need for the management of these forests. Some authors [59] have grouped the silvicultural strategies for the conversion of coniferous plantations. In line with these recommendations, the experimental plots could be grouped according to: (1) the capability for stand conversion after one harvest cycle of black pine, or raising new forest cultures on the same habitat (epi-Mediterranean or higher sections of the sub-Mediterranean zone); (2) stand conversion is possible, though gradually, and with limited and smaller works aimed at obtaining a mixed deciduous and coniferous forest, with later support of succession towards the climazonal vegetation (sub-Mediterranean zone); and (3) conversion after one rotation period is not possible, and at least one more rotation period of black pine is needed, the growth and development of black pine is slowed due to the unfavourable habitat conditions and higher slopes (lower sections of sub-Mediterranean zone towards the eu-Mediterranean zone). Most of the current issues and hindering conditions for conversion are tied to the lack of any silviculture works (cleaning) of forest plantations in the early years and thinning in later years of the pine forest plantation.

5. Conclusions

The sub-Mediterranean habitats in western Croatia are recognisable for two species. Pubescent oak (Quercus pubescens L.) is the climatogenic deciduous species and black pine (Pinus nigra J.F. Arnold) is the main forestation species of this area. The differences in the synecological factors in the research area resulted in the presence of sub-Mediterranean, epi-Mediterranean and eu-Mediterranean elements within the black pine forest plantations. With increasing elevation in zones with higher precipitation, the ecological conditions were better, and the production potential was also more pronounced in those habitats. In the study area, this production potential needs to be improved and better utilised. Forest plantations in which silvicultural works are periodically carried out are better and stronger in the qualitative sense. This also reduces the risk of the occurrence and development of large wildfires. The lack of these works reduces the return of the indigenous deciduous vegetation, while also creating greater quantities of forest fuel (dry lower branches, dry standing trees due to density). These facts should be considered when aiming to achieve higher quality forest plantations of black pine. The research results and previous experience indicate that within the epi-Mediterranean zone of Croatia, new forest cultures of black pine can be raised (after one rotation period), as they can fulfil both the desired ecological and commercial roles. In the sub-Mediterranean zone of Croatia, the conditions are less favourable (edaphic conditions, higher slopes, strong influence of wind), and therefore greater productivity of black pine is lacking these areas.

Author Contributions

Conceptualization, D.B.; methodology, D.B.; software, V.H.; validation, V.H. and D.B.; formal analysis, V.H. and D.B.; investigation, D.B., V.H. and M.A.; resources, D.B. and R.R.; data curation, D.B.; writing—original draft preparation, D.B.; writing—review and editing, D.B. and V.H.; visualization, M.A. and R.R.; supervision, V.H.; project administration, D.B.; funding acquisition, R.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank the employees of the public enterprise „Hrvatske šume“ d.o.o. in the field of research of research (Croatian Forests). This study was supported by human and material resources of public enterprise „Hrvatske šume“ d.o.o.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Quézel, P.; Médail, F. Ecologie et Biogéographie du Bassin Méditerranéen; Elsevier: Amsterdam, The Netherlands, 2003. [Google Scholar]
  2. Barbéro, M.; Oisel, R.; Quézel, P.; Richardson, D.; Romane, F. Pines of the Mediterranean Basin. In Ecology and Biogeography of Pinus; Richardson, D.M., Ed.; Cambridge University Press: Cambridge, UK, 1998. [Google Scholar]
  3. Lucas-Borja, M.E.; Jing, X.; Candel-Perez, D.; Parhizkar, M.; Rocha, F.; Heydari, M.; Muñoz-Rojas, M.; Zema, D.A. Afforestation with Pinus nigra Arn ssp salzmannii along an elevation gradient: Controlling factors and implications for climate change adaptation. Trees 2022, 36, 93–102. [Google Scholar] [CrossRef]
  4. Ivančević, V. Biološko-tehnički radovi na sanaciji Senjske bujice «torrente» i povećanje vodnog kapaciteta. Šum. List. 2005, 129, 91–109. (In Croatian) [Google Scholar]
  5. Piermattei, A.; Renzaglia, F.; Urbinati, C. Recent expansion of Pinus nigra Arn. above the timberline in the central Apennines, Italy. Ann. For. Sci. 2012, 69, 509–517. [Google Scholar] [CrossRef] [Green Version]
  6. Piermattei, A.; Garbarino, M.; Urbinati, C. Structural attributes, tree-ring growth and climate sensitivity of Pinus nigra Arn. At high altitude: Common patterns of a possible treeline shift in the central Apennines (Italy). Dendrochronologia 2014, 32, 210–219. [Google Scholar]
  7. Piermattei, A.; Lingua, E.; Urbinati, C.; Garbarino, M. Pinus nigra anthropogenic treelines in the central Apennines show common pattern of tree recruitment. Eur. J. For. Res. 2016, 135, 1119–1130. [Google Scholar] [CrossRef] [Green Version]
  8. Santana, V.M.; Alday, J.G.; Baeza, M. Effects of fire regime shift in Mediterranean Basin ecosystems: Changes in soil seed bank composition among functional types. Plant Ecol. 2014, 215, 555–566. [Google Scholar] [CrossRef]
  9. Santana, V.M.; González-Pelayo, J.G.; Maia, P.E.E.; Varela, T.M.E.; Valdecantos, A.; Vallejo, V.R.; Keizeri, J.J. Effects of fire recurrence and different salvage logging techniques on carbon storage in Pinus pinaster forests from northern Portugal. Eur. J. For. Res. 2016, 135, 1107–1117. [Google Scholar] [CrossRef]
  10. Muñoz-Rojas, M.; De la Rosa, D.; Zavala, L.; Jordán, A.; Anaya-Romero, M. Changes in land cover and vegetation carbon stocks in Andalusia, Southern Spain (1956–2007). Sci. Total Environ. 2011, 409, 2796–2806. [Google Scholar] [CrossRef]
  11. Morales-Molino, C.; Tinner, W.; García-Antón, M.; Colombaroli, D. The historical demise of Pinus nigra forests in the Northern Iberian Plateau (south-western Europe). J. Ecol. 2016, 105, 634–646. [Google Scholar] [CrossRef] [Green Version]
  12. Sheffer, E. A review of the development of Mediterranean pine-oak ecosystems after land abandonment and afforestation: Are they novel ecosystems? Ann. For. Sci. 2012, 69, 429–443. [Google Scholar] [CrossRef] [Green Version]
  13. Espelta, J.M. La reconstruccio de paisatges forestals afectats per grans incendis. Projecte pilot a l’incendi del Bages-Bergueda del 1994. Silvicultura 1999, 27, 6–9. [Google Scholar]
  14. Espelta, J.M.; Retana, J.; Habrouk, A. An economic and ecological multi-criteria evaluation of reforestation methods to recover burned Pinus nigra forests in NE Spain. For. Ecol. Manag. 2003, 180, 185–198. [Google Scholar] [CrossRef]
  15. Vallejo, V.R.; Bautista, S.; Cortina, J. Restoration for soil protection after disturbances. In Life and Environment in the Mediterranean; Trabaud, L., Ed.; WIT Press: Boston, MA, USA, 2000; pp. 301–343. [Google Scholar]
  16. Allen, H.D. Response of past and present Mediterranean ecosystems to environmental change. Prog. Phys. Geogr. 2003, 27, 359–377. [Google Scholar] [CrossRef]
  17. González-Hidalgo, J.; De Luis, M.; Raventós, J.; Sánchez, J. Spatial distribution of seasonal rainfall trends in a western Mediterranean area. Int. J. Climatol. 2001, 21, 843–860. [Google Scholar] [CrossRef]
  18. Martinović, J. Gospodarenje Šumskim Tlima u Hrvatskoj; Šumarski Institut Jastrebarsko i Hrvatske Šume d.o.o.: Zagreb, Croatia, 2003; p. 521. (In Croatian) [Google Scholar]
  19. Trinajstić, I. Nomenklaturno-sintaksonomska revizija submediteranskih šuma crnoga bora (Pinus nigra Arnold) u Hrvatskoj. Šum. List. 1998, 122, 147–154. (In Croatian) [Google Scholar]
  20. Rauš, Đ.; Trinajstić, I.; Vukelić, J.; Medvedović, J. Biljni svijet hrvatskih šuma. In Šume u Hrvatskoj.; Rauš, Đ., Ed.; Šumarski Fakultet Sveučilišta u Zagrebu i Hrvatske Šume: Zagreb, Croatia, 1992; pp. 33–77. (In Croatian) [Google Scholar]
  21. Liber, Z. Filogenetski i Taksonomski Odnosi Populacija Crnoga Bora (Pinus nigra Arnold) u Hrvatskoj. Ph.D. Thesis, Prirodoslovno-Matematički Fakultet Sveučilišta u Zagrebu (Faculty of Science University of Zagreb), Zagreb, Croatia, 2000. [Google Scholar]
  22. Gaussen, H.; Heywood, V.H.; Chater, A.O.; Pinus, L. Flora Europaea, 2nd ed.; Tutin, I.G., Heywood, V.H., Burges, N.A., Moore, D.M., Valentine, D.H., Walters, S.M., Webb, D.A., Eds.; Cambridge University Press: Cambridge, UK, 1993; p. 42. [Google Scholar]
  23. Vidaković, M. Oblici crnog bora u Jugoslaviji na temelju antomije iglica. Glas. Šum. Pokuse 1957, 13, 111–148. (In Croatian) [Google Scholar]
  24. Yildiz, O.; Altundağ, E.; Çetin, B.; Güner, S.T.; Sarginci, M.; Toprak, B. Experimental arid land afforestation in Central Anatolia, Turkey. Environ. Monit. Assess. 2018, 190, 355. [Google Scholar] [CrossRef]
  25. Bayar, E.; Deligöz, A. Ecophysiological behavior of Mediterranean woody species under summer drought. Bosque 2021, 42, 311–321. [Google Scholar] [CrossRef]
  26. Navarro-Cerrillo, R.M.; Rodriguez-Vallejo, C.; Silveiro, E.; Hortal, A.; Palacios-Rodríguez, G.; Duque-Lazo, J.; Camarero, J.J. Cumulative Drought Stress Leads to a Loss of Growth Resilience and Explains Higher Mortality in Planted than in Naturally Regenerated Pinus pinaster Stands. Forests 2018, 9, 358. [Google Scholar] [CrossRef] [Green Version]
  27. Zaghi, D. Management of Natura 2000 Habitats. 9530 *(Sub)-Mediterranean Pine Forests with Endemic Black Pines; Technical Report 2008 24/24; European Commission: Brussels, Belgium, 2008. [Google Scholar]
  28. Bezak, K. Tablice drvnih masa cera, crnog bora i običnog bora. Radovi 1992, 5, 47–65. (In Croatian) [Google Scholar]
  29. Braun-Blanquet, J. Grundzüge der Vegetationskunde. In Pflanzensoziologie, 3rd ed.; Springer: Wien, Austria, 1964. [Google Scholar]
  30. Dierschke, H. Grundlagen und Methoden. In Pflanzensoziologie; Verlag Eugen Ulmer: Stuttgart, Germany, 1994. [Google Scholar]
  31. Barkman, J.J.; Doing, H.; Segal, S. Kritische Bemerkungen und Vorschläge zur Quantitativen Vegetationsanalysen. Acta Bot. Neerl. 1964, 13, 394–419. [Google Scholar] [CrossRef]
  32. van der Maarel, E. Transformation of cover-abundance values in phytosociology and its effect on community similarity. Vegetatio 1979, 39, 97–114. [Google Scholar]
  33. Nikolić, T. (Ed.) Popis flore Hrvatske, 1. dio. Flora Croatica. Index Florae Croaticae. Pars 1 Nat. Croat. 1994, 3 (Suppl. S2), 1–116. [Google Scholar]
  34. Nikolić, T. (Ed.) Popis flore Hrvatske, 2. dio. Flora Croatica. Index Florae Croaticae. Pars 2 Nat. Croat. 1997, 6 (Suppl. S1), 1–232. [Google Scholar]
  35. Nikolić, T. (Ed.) Popis flore Hrvatske, 3. dio. Flora Croatica. Index Florae Croaticae. Pars 3 Nat. Croat. 2000, 9 (Suppl. S1), 1–324. [Google Scholar]
  36. United Nations Economic Commission for Europe. Manual on Methods and Criteria for Harmonized Sampling, Assessment, Monitoring and Analysis of the Effects of Air Pollution on Forests; Part XI, Sampling and Analysis of Litterfall. 2004, pp. 1–18. Available online: https://www.icp-forests.org/pdf/manual/2000/Chapt11_compl2004.pdf (accessed on 17 September 2022).
  37. McCune, B.; Grace, J.B. Analysis of Ecological Communities; MjM Software Design: Gleneden Beach, OR, USA, 2002. [Google Scholar]
  38. Dufrene, M.; Legendre, P. Species assemblages and indicator species: The need for a flexible asymmetrical approach. Ecol. Monogr. 1997, 67, 354–366. [Google Scholar] [CrossRef]
  39. Roberts, D.W. Statistical analysis of multidimensional fuzzy set ordinations. Ecology 2008, 89, 1246–1260. [Google Scholar] [CrossRef]
  40. Roberts, D.W. Comparison of multidimensional fuzzy set ordination with CCA and DB-RDA. Ecology 2009, 90, 2622–2634. [Google Scholar] [CrossRef]
  41. Ellenberg, H. Wege der Geobotanik zum Verständnis der Pflanzendecke. Naturwissenschaften 1968, 55, 462–470. [Google Scholar] [CrossRef]
  42. Tomašević, A. Meliorativni utjecaj kulture alepskog bora (Pinus halepensis Mill.) i pinije (Pinus pinea Endl.) na degradirano stanište hrasta medunca (Querco pubescentis-Carpinetum orientalis H-ić, Anić 1959) u zadarskom području. Glas. Šum. Pokuse 1994, 30, 223–298. (In Croatian) [Google Scholar]
  43. Enríquez-De-Salamanca, Á. Effects of Climate Change on Forest Regeneration in Central Spain. Atmosphere 2022, 13, 1143. [Google Scholar] [CrossRef]
  44. Allen, C.D.; Macalady, A.K.; Chenchouni, H.; Bachelet, D.; McDowell, N.; Vennetier, M.; Kitzberger, T.; Rigling, A.; Breshears, D.D.; Hogg, E.H.; et al. A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. For. Ecol. Manag. 2010, 259, 660–684. [Google Scholar] [CrossRef] [Green Version]
  45. Trew, B.T.; Maclean, I.M.D. Vulnerability of global biodiversity hotspots to climate change. Glob. Ecol. Biogeogr. 2021, 30, 768–783. [Google Scholar] [CrossRef]
  46. Gil, L.; Prada Aránzazu, M. Los pinos como especies basicas de la restauracion forestal en el medio mediterráneo. Ecologia 1993, 7, 113–125. [Google Scholar]
  47. Serrada Hierro, R. Consideraciones sobre el impacto de la repoblacion forestal en el suelo. Ecologia 1990, 1, 453–462. [Google Scholar]
  48. Grau Corbi, J.M.; Cámara Obregón, A.; Montoto Quinteiro, J.L. Fitoclimatología básica de los Pinus nigra Arn., Pinus sylvestris L. Y Pinus pinaster Ait. Aplicación del modelo de idoneidad. Investig. Agrar. Sist. Recur. For. 1999, 1, 37–51. [Google Scholar]
  49. Quézel, P. Les forêts du pourtour méditerranéen. In Forêts et Maquis Méditerranéens: Écologie, Conservation et Aménagements; Note technique MAB; UNESCO: France, Paris, 1976; Volume 2, pp. 9–33e. [Google Scholar]
  50. M’hirit, O. El bosque mediterráneo: Espacio ecológico, riqueza económica y bien social. Unasylva (Mediterranean Forests, FAO) 1999, 50, 3–15. [Google Scholar]
  51. Cachinero-Vivar, A.M.; Palacios-Rodríguez, G.; Lara-Gómez, M.A.; Navarro-Cerrillo, R.M. Thinning Effect of C Sequestration along an Elevation Gradient of Mediterranean Pinus spp. Plantations. Forests 2021, 12, 1583. [Google Scholar] [CrossRef]
  52. Pardos, M.; Vázquez-Piqué, J.; Benito, L.; Madrigal, G.; Alejano, R.; Fernández, M.; Calama, R. Does the Age of Pinus sylvestris Mother Trees Influence Reproductive Capacity and Offspring Seedling Survival? Forests 2022, 13, 937. [Google Scholar] [CrossRef]
  53. Tesei, G.; Allegrezza, M.; Ballelli, S.; Ciaschetti, G.; Console, C.; Montecchiari, S.; Ottaviani, C.; Pirone, G.; Frattaroli, A.R. The oldest Pinus nigra plantations in the central Apennines: Distribution and floristic-vegetational characterisation. Plant Sociol. 2021, 58, 15–28. [Google Scholar] [CrossRef]
  54. Zerbe, S. Restoration of natural broad-leaved woodland in Central Europe on sites with coniferous forest plantations. For. Ecol. Manag. 2002, 167, 27–42. [Google Scholar] [CrossRef]
  55. Hasenauer, H. Terms and Definitions Relevant for Conservation. Abstracts of the International Conference. In Proceedings of the Question of Conversion of Coniferous Forests, Freiburg, Germany, 27 September–2 October 2003; p. 16. [Google Scholar]
  56. Spiecker, H.; Hansen, J.; Klima, E.; Skovsgaard, J.; Sterba, H.; Teuffel, K. Norway Spruce Conversion-Options and Consequences; European Forest Institute Research Report: Joensuu, Finland, 2004; p. 18. [Google Scholar]
  57. Vizcaíno-Palomar, N.; Revuelta-Eugercios, B.; Zavala, M.A.; Alía, R.; González-Martínez, S.C. The role of population origin and micro environment in seedling emergence and early survival in Mediterranean maritime pine (Pinus pinaster Aiton). PLoS ONE 2014, 9, e109132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Zlatanov, T.; Velichkov, I.; Lexer, M.J.; Dubravac, T. Regeneration dynamics in aging black pine (Pinus nigra Arn.) plantations on the south slopes of the Middle Balkan Range in Bulgaria. New For. 2010, 40, 289–303. [Google Scholar] [CrossRef]
  59. Lüpke, B.; Ammer, C.; Bruciamacchie, M.; Brunner, A.; Ceitel, J.; Collet, C.; Deuleuze, C.; Placido, J.; Huss, J.; Jankovič, J.; et al. Silvicultural Strategies for Conversion. Abstract of the International Conference. In Proceedings of the Question of Conversion of Coniferous Forests, Freiburg, Germany, 27 September–2 October 2003; pp. 25–26. [Google Scholar]
Forests 13 01900 g001 550
Figure 1. Distribution of plots in the study area in Croatia.
Figure 1. Distribution of plots in the study area in Croatia.
Forests 13 01900 g001
Forests 13 01900 g002 550
Figure 2. Black pine stands ((A)—stand on 196 m, (B)—stand on 539 m, (C)—stand on 716 m).
Figure 2. Black pine stands ((A)—stand on 196 m, (B)—stand on 539 m, (C)—stand on 716 m).
Forests 13 01900 g002
Forests 13 01900 g003 550
Figure 3. Cluster analysis of experimental plots obtained results are arbitrarily divided into five clusters.
Figure 3. Cluster analysis of experimental plots obtained results are arbitrarily divided into five clusters.
Forests 13 01900 g003
Forests 13 01900 g004 550
Figure 4. Biplot of multivariate fuzzy set ordination superponed with cluster analysis, the x-axis biplex represents the height gradient and the y-gradient represent nutritient gradient.
Figure 4. Biplot of multivariate fuzzy set ordination superponed with cluster analysis, the x-axis biplex represents the height gradient and the y-gradient represent nutritient gradient.
Forests 13 01900 g004
Forests 13 01900 g005 550
Figure 5. Shows sum of mean cover values for species of sub-Mediterranean plant species.
Figure 5. Shows sum of mean cover values for species of sub-Mediterranean plant species.
Forests 13 01900 g005
Table
Table 1. Vegetation and climate properties of the study area (Source: research and customized data of State hydrometeorological institute).
Table 1. Vegetation and climate properties of the study area (Source: research and customized data of State hydrometeorological institute).
Vegetation BeltVegetation ZonePhytocenosisAverage Annual PrecipitationAverage Annual TemperaturePioneer SpeciesClimatogenic Species
Mediterranean-litoralSub-MediterraneanQuerco pubescenti-Carpinetum orientalis Horvatić 19391200 mm12–14 °CPinus nigraCarpinus orientalis, Quercus pubescens
Mediterranean-montanEpi-MediterraneanAristolochio luteae-Quercetum pubescentis (Horvat 1959) Poldini 20081400 mm10–12 °CPinus nigraOstrya carpinifolia, Quercus pubescens
Table
Table 2. Results of multivariate fuzzy set ordination indicate important values of altitude and slope.
Table 2. Results of multivariate fuzzy set ordination indicate important values of altitude and slope.
Variables Cum-R2 Inc.-R2 Random p GammaR2-Max R2-Min
Altitude 0.19730.19730.00101.00000.19110.0001
Slope 0.32100.10720.06900.89720.34110.0023
Org mat 0.39220.07120.19600.22510.36850.0186
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Barčić, D.; Hršak, V.; Rosavec, R.; Ančić, M. Ecological Potential of Mediterranean Habitats for Black Pine (Pinus nigra J.F. Arnold) in Croatia. Forests 2022, 13, 1900. https://doi.org/10.3390/f13111900

AMA Style

Barčić D, Hršak V, Rosavec R, Ančić M. Ecological Potential of Mediterranean Habitats for Black Pine (Pinus nigra J.F. Arnold) in Croatia. Forests. 2022; 13(11):1900. https://doi.org/10.3390/f13111900

Chicago/Turabian Style

Barčić, Damir, Vladimir Hršak, Roman Rosavec, and Mario Ančić. 2022. "Ecological Potential of Mediterranean Habitats for Black Pine (Pinus nigra J.F. Arnold) in Croatia" Forests 13, no. 11: 1900. https://doi.org/10.3390/f13111900

APA Style

Barčić, D., Hršak, V., Rosavec, R., & Ančić, M. (2022). Ecological Potential of Mediterranean Habitats for Black Pine (Pinus nigra J.F. Arnold) in Croatia. Forests, 13(11), 1900. https://doi.org/10.3390/f13111900

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