Aboveground Biomass of Living Trees Depends on Topographic Conditions and Tree Diversity in Temperate Montane Forests from the Sl ă tioara-Rar ă u Area (Romania)

: The study zone includes one of the largest montane old-growth forests in Europe (Slatioara UNESCO site), and understanding the structure and functioning of sill intact forests in Europe is essential for grounding management strategies for secondary forests. For this reason, we set out to analyze the dependencies between aboveground biomass (AgB), tree species and size diversity and terrain morphology, as well as the relationship between biomass and diversity, since neither of these issues have been sufﬁciently explored. We found that tree species diversity decreases with increased solar radiation and elevation. Tree size heterogeneity reaches its highest mean values at elevations between 1001 and 1100 m, on slopes between 50 and 60 degrees. AgB is differentiated with elevation; the highest mean AgB (293 tonnes per hectare) is recorded at elevations between 801 and 900 m, while it decreases to 79 tonnes per hectare at more than 1500 m a.s.l. It is also inﬂuenced by tree species diversity and tree size heterogeneity, with the highest AgB reached in the most complex forest ecosystems in terms of structural diversity. We showed that intact temperate montane forests develop maximum biomass for optimum species diversity and highest size heterogeneity; all three are modulated mainly by elevation. forest compartment, management type, elevation (m), slope (degrees), aspect (degrees), solar radiation (KWh · m − 2 ) during the growing season, N2000 habitat type, number of tree species, stem density, species diversity, tree size heterogeneity, modiﬁed Lorenz curve asymmetry, aboveground wood volume (m 3 ), and aboveground biomass (tonnes).


Introduction
Intact forests are important repositories of precious irreplaceable and often endangered biodiversity, and they provide a wide range of valuable ecosystem services [1]. The reasons behind the increasing interest in studying intact forests are four-fold: (1) the need to improve the management of second-growth forests [2][3][4][5][6]; (2) the sheer scientific curiosity in studying the most complex natural terrestrial ecosystems still at hand in many geographical regions of the Earth [7,8]; (3) the need to protect fragile habitats, many of them hosted by intact forests [9,10]; and (4) the issue of forest carbon pool [11,12].
Forest biodiversity goes hand in hand with diversified horizontal and vertical structures and the so-called close-to-nature forestry relies to a large extent on the observations and measurements carried out in strict forest reserves [13] or forests affected by disturbances [14]. It is generally accepted that tree diversity enhances high forest productivity and ensures at least moderate levels of ecosystem multifunctionality [15][16][17][18][19][20][21], but the factors influencing this process and the environmental mechanisms that drive it are not fully known and are even controversial [22][23][24][25][26][27].
Global scale studies reveal a consistent positive relationship between tree species diversity and productivity [16], using tree species richness as surrogate for biodiversity because taxonomic diversity is considered to incorporate other traits such as functional, phylogenetic, and genomic diversity [16,28]. One of the complex and multiple implications of forest biodiversity is that globally, forests will produce less biomass and accordingly, their global capacity to absorb and store carbon will decrease [29,30]. Other large-scale studies (in North America and Europe) also indicate a positive relationship between species richness and tree productivity, while addressing the influence of environmental and climatic conditions [19,31]; the role of biodiversity tends to be stronger in harsher climates (dry Mediterranean or boreal) [15,32].
Because both biomass and biodiversity depend on a wide spectrum of factors, some controlled and some uncontrolled through forest management, the scientific hypotheses different scholars have put forward mostly refer to the interplayed roles of forests as carbon sinks [33,34] and shelters for biodiversity [4,35]. An important challenge for forest management and conservation is achieving a proper balance between biodiversity and the aboveground biomass. The adequate balance, wherever it may occur, may help forest managers and policy makers to develop managerial scenarios either for conservation purposes [7,36] or for increasing the resilience of the managed forests [37,38].
A number of studies have highlighted that structurally more complex forests may have higher productivity [39], and hence are at least as profitable as uniform forests [40], they can promote biodiversity of certain animal taxa [41] and are more resilient against disturbances and climate change [20,42]. Consequently, to provide as many ecological services as possible while maintaining their resilience, it is recommended that the composition and structure of managed forests be diversified. Preparing forests for future changing conditions may imply enriching the composition with species complementary in their functional traits [20,43], or controlling stand density, which can be an effective way of managing water availability in expected drier climates [44].
In research performed at smaller scales, in addition to the abovementioned dependencies, forest structure, canopy structural complexity [18,45,46] and the environmental factors modulate the interactions between forest diversity and productivity/biomass [47,48]. For example, it was found that solar radiation impacts on forest primary productivity [49] and influences species richness [50]. Additionally, there is a rich literature documenting: elevation impact on forest diversity [51,52] and forests biomass [53][54][55]; relationship between aspect and biomass [56,57]; or the influence of slope on the forest aboveground biomass [58].
Given the rare opportunity to study these interactions in natural forests located on sites with a wide range of topographical conditions, knowing that they are not fully understood and that Romania has 8% of Europe's intact forests [8], we intend to exploit the existence of a permanent research platform (PRP) within the Slătioara UNESCO site and its surroundings [59] and to expand the research methodology to a wider network of intact forest sites. Through this analysis, we aim to highlight interaction patterns between aboveground biomass, biodiversity and topographical factors, and patterns transcending the micro-site context.
Thus, the following research objectives were addressed: • assessing tree species diversity, tree size heterogeneity and aboveground biomass of living trees at sample plot level; • identifying the effect of some abiotic factors on tree species diversity and aboveground biomass; • identifying how aboveground biomass depends on tree species diversity and size heterogeneity, within the entire study area and at habitat type level.

Study Location
The study is located in Natura 2000 (N2000) site ROSCI0212 Rarău-Giumalău, covering 2205.85 ha, with 1724.57 ha of forests and overlapping Slătioara UNESCO site by 1038.55 ha (of which the core area occupies 609.12 ha) (Figure 1). It is situated in northern Romania, in the North-Eastern Carpathians, on the south-eastern side of the Rarău Mountains, at altitudes between 800 and 1510 m a.s.l. The climate is continental-temperate, with mean annual temperatures between 3.8 • C at high elevations and 5.9 • C at low elevations, and mean annual rainfall between 700 and 810 mm. Duduman et al. [59] provide a detailed description of the study area, in terms of biodiversity and topography. Figure 1. Study location-Slătioara-Rarău area (reserve boundaries from http://mmediu.ro/articol/date-gis/434, (accessed on 15 October 2020). overlaid on the satellite images from ESRI ArcGIS server). Maps throughout this book were created using ArcGIS ® software by Esri. ArcGIS ® and ArcMap™ are the intellectual property of Esri and are used herein under license. Copyright © Esri. All rights reserved. For more information about Esri ® software, please visit www.esri.com (accessed on 15 October 2020).
The landscape is characterized by mixed temperate forests with coniferous trees and European beech at low and middle elevations, and by pure Norway spruce stands at high elevations, where they alternate with pasture lands. The forest included in the UNESCO core area was never subjected to logging operations (609.12 ha). After 2000, when the "Rarău-Pietrele Doamnei" Nature Reserve was established, the strictly protected area increased to 1580.73 ha. The inclusion was feasible due to the fact that the past management was based on very low intensity interventions, with no discernable alterations of natural stand structure. The rest of the forest area (143.84 ha) is still under regular management, under the provisions of the Rarău-Giumalău Natura 2000 site management plan.

Data Collection
Data were collected in 2015 from 193 sample plots, 500 m 2 each (horizontal area), within a permanent research platform (PRP) covering the 1724.57 ha of forests and structured on two levels: the first level of PRP covers the entire area and corresponds to a square grid of 500 × 500 m; the second level of PRP was set up for the core area, according to a square grid of 100 × 100 m (see Figure 3 in [59]). All (living or dead) standing trees with diameter at breast height (dbh) larger than 5 cm were inventoried within every sample plot (SP) but, for this study, only the data regarding the 8296 living standing trees were analyzed. Every tree was tallied at SP level. In order to quantify the aboveground wood biomass, the species was identified, dbh was measured using an appropriate caliper (graded in millimeters) and the tree total height was measured and expressed in decimeters, using a Vertex dendrometer (Haglöf Sweden).
Additionally, in every SP, the forest habitat type was assessed, considering the N2000 habitat classification [60][61][62]. Thus, additional field observations and analysis of forest management plans (FMPs) [63][64][65] were carried out, aiming to determine-at forest compartment level-the habitat type, based on site conditions (elevation, aspect, slope, soil type), management type and management history, vascular species composition of the trees, shrub and herbaceous layers.

Data Processing
The data from the field inventory (FI) were cross-referenced with the forest management planning database (FMPD) and the digital terrain model to extract several site and vegetation data that were subsequently analyzed in relation to aboveground biomass and biodiversity ( Figure 2).

GIS Data Processing
The Digital Elevation Model (DEM) was built using contour lines, elevation points and stream lines extracted from the georeferenced topographic plans (scale 1:5000) within ArcGIS 9.3 software platform. The Topo to Raster module was used to generate the elevation model, which was further used to compute slope and aspect, within the 3D Analyst extension of ArcGIS, for a 10 m cell size.
The position of each SP center was used to extract the corresponding values of elevation, slope and aspect from the three raster files mentioned above within the Zonal statistics as a table function of ArcGIS. The output tables were joined with the initial shapefile containing the SPs. The aspect was expressed in degrees, from 0.1 to 360.0 but, for a proper interpretation of the results, the values larger than 180 degrees were corrected by subtracting them from 360 degrees. Thus, in this case, the values close to 0 correspond to shaded slopes, while those close to 180 degrees indicate sunny slopes. The incident solar radiation per area unit was computed using the Solar Analyst from ArcGIS 9.3, using the georeferenced Digital Elevation Model, average parameters for cloud cover, and the mean length of the growing season considered from 1 April to 15 October.
The forest compartment attributes from the management planning database were assigned correspondingly to SPs using the spatial join function.
For further data processing, discrete classes were distinguished as follows: 100 KWh·m −2 class size in the case of solar radiation, 100 m class size in the case of elevation, 10 degrees class size for slope, 22.5 degrees class size for aspect, 0.4 class size for the species diversity index, and 0.1 class size in the case of the Gini index.

Diversity Indices Computation
The species diversity (D) in every SP was derived from the Shannon [66] index (H), using the exponential of Renyi entropy [67]: (1) where k is the number of species within a given SP; n i represents the number of individuals of the ith species; and N is the total number of individuals within a SP. The D index ranges between 1 (all individuals are form the same species) and k (all individuals are evenly spread between the k species). Tree size heterogeneity was assessed, as recommended in the literature [68][69][70], by using the Gini index, which ranges between 0 and 1. The more homogenous a population is, the smaller the Gini index is [71,72]. The formula used for computing the Gini index [70] is: where b ai (b ai−1 ) is the cumulative fraction of the basal area (%) of the trees from all diameter classes thinner than or equal to the ith (i−1) diameter class (for i = 1, b ai-1 = 0); n i (n i−1 ) is the cumulative fraction of the number of trees (%) from all diameter classes thinner than or equal to the ith (i − 1) diameter class (for i = 1, n i-1 = 0); and k represents the number of 2 cm-diameter classes. For a better interpretation of the results related to size heterogeneity at SP level, the modified Lorenz [73] curve asymmetry (LA) was computed, considering the proportion of basal area accumulated by the trees with dbh above the quadratic mean diameter (M(xQMD)) and the proportion of stem density accumulated by the trees with dbh above the quadratic mean diameter (xQMD) [74]: The quadratic mean diameter (cm) per SP (QMD) is computed with Formula (5): where m = π·10 4 /4, g is the total basal area of a target SP (m 2 ), and n is total number of trees in that SP. As recommended by [70,75], the 0.51 value of Gini index was used as a threshold for separating the uneven-sized balanced forest stands [76] from those with a more simplified size heterogeneity, and also the border corresponding to LA = 0.5 (established by [74]) for differentiating the bimodal plots from reverse J-shaped plots.

Aboveground Biomass Assessment
The total aboveground wood volume per SP (V) was estimated by summing every single tree volume, regardless of species, and considering the dbh threshold of 5 cm: where: k is the number of species within a SP; n is the total number of trees within a SP; vu i;j is the single tree aboveground wood volume, which was computed for the jth single tree from ith species, using the logarithmic equation recommended in the literature [77] with the specific coefficients a 0i , a 1i , a 2i , a 3i , a 4i (Table 1) for the ith species: log 10 vu i;j = a 0i + a 1i · log 10 dbh j + a 2i · log 10 dbh j 2 + a 3i · log 10 h j + a 4i · log 10 h j 2 (7) where dbh j represents the diameter at breast height of the jth tree; h j is the total height of the jth tree. Table 1. Coefficients used for computing the tree aboveground wood volume (according to [77]) and aboveground biomass (according to [78]). W ID is the mean value of wood infra-density. Having computed the aboveground wood volume per tree, the aboveground biomass computed for the jth single tree from ith species (AgB i;j ) is:

Tree Species (i)/Code Coefficients for Equation (7)
where W ID-j is the mean value of wood infra-density (tonnes per cubic meter) established by tree species in [78] (Table 1). Finally, a database was obtained including information for every SP, referring to: SP number, forest district name, forest compartment, management type, elevation (m), slope (degrees), aspect (degrees), solar radiation (KWh·m −2 ) during the growing season, N2000 habitat type, number of tree species, stem density, species diversity, tree size heterogeneity, modified Lorenz curve asymmetry, aboveground wood volume (m 3 ), and aboveground biomass (tonnes).

Data Analysis
In order to compare the distribution of variables among the inventoried SPs, a principal component analysis (PCA) was used as an exploratory data analysis, considering the number of tree species, stem density, species diversity (D), tree size heterogeneity (G), modified Lorenz curve asymmetry, aboveground biomass, elevation, aspect and slope as variables and the 193 plots as observations. The PCA was applied for two different situations: (i) considering all SPs; (ii) considering only the first level of PRP [59]. The FactoMineR package (version 2.3) [79] of R (version 3.6.2) was used.
In relation to the results of PCA, we used only the first level PRP and we selected the next pairs of variables for a more detailed analysis: tree species richness and diversity, tree size heterogeneity and AgB on one hand and elevation, slope and aspect on the other hand. Additionally, we noted the relationship between biomass and tree species and size diversity. The statistical analysis of differences between habitats, data regarding the interaction between tree's diversity and AgB were further detailed separately per Natura 2000 habitat type.
The analysis of the differences among means, in the case of total biomass, number of species, D, and Gini, depending on five factors (habitat type and elevation, slope, aspect, radiation classes), was carried out using the Kruskal-Wallis nonparametric test (because all experimental distributions differ from the normal distribution (Shapiro-Wilk test), even after the data transformation (x' = log (x + 1)). When significant differences were found, the Conover-Iman procedure for multiple pairwise comparations was used to compare means [80]. This analysis was performed with XLSTAT-PRO 2012 (Addinsoft, New York, NY, USA), plugged into EXCEL 2016 (Microsoft Corp., Redmond, Washington, DC, USA). Data processing was performed by employing the ArcGIS 9.3 and R software [81]. Maps throughout this paper were created using ArcGIS ® software (Esri, West Redlands, CA, USA).

Overview of Tree Species, Forest Habitats and Aboveground Biomass
Eleven tree species represented by 8296 living individuals and three N2000 forest habitats were found in 178 SPs, while the remaining 15 SPs were assigned to other types (non-Natura 2000 forest habitats (NN)) ( Table 2). In only one SP (outside core area and classified as NN habitat type), tree composition included non-local species (European larch and Scots pine), artificially introduced 100 years ago [65]. All the 135 SPs from second level PRP, installed in the core area of the Slătioara UNESCO site were assigned to N2000 habitat types, most of them (88) as Dacian Beech forests (Symphyto-Fagion) (91V0), 35 SPs correspond to Luzulo-Fagetum beech forests (9110) and 12 to Acidophilous Picea forests of the montane to alpine levels (Vaccinio-Piceetea) (9410). If considering only the first level of PRP [59], 30 SPs correspond to 9410, 13 SP to 91V0 and 15 SP are NN.
All plots from UNESCO core area (135 SPs from the second level and 12 SPs from the first level PRP) have never been subjected to any logging operations. Another 31 SPs from the first level PRP overlap the buffer area of UNESCO site and the "Rarău-Pietrele Doamnei" Nature Reserve and have been excluded from logging since 2000. Only 15 SPs, classified as NN habitats, are still under regular management.
The Acidophilous Picea forests contain the lowest mean AgB (200.6 tonnes per hectare), while the beech forests reach about 275 tonnes per hectare. The most important three species in descending order of AgB in the habitat type 9110 are Norway spruce (35.8%), European beech (30.8%) and silver fir (29.0%). Within the habitat type 91V0, the European beech scores the highest AgB (36.8%), followed by silver fir (29.0%) and Norway spruce (22.8%), but the other species have an almost triple cumulated AgB compared to the 9110 habitat (11.3% versus 4.4%) and one additional species (wych elm). In the case of Acidophilous Picea forests, the main species is Norway spruce (more than 55% of AgB), followed by silver fir (24.3%), sycamore (10.4%) and European beech (7.6%).

Relationships between the Attributes of the SPs
In the case of the entire PRP (Figure 3a), the variance of principal components is 51.2%, mainly due to the high number of SPs and variables. The variables "species_no", G and D, are positively correlated with the first dimension, while "elevation" and LA are negatively correlated. The variables "stem_density" and "aspect" are positively correlated with the second dimension. Thus, the stem density and aspect are not related to either species diversity or tree size heterogeneity. More than a half of inventory plots are similar, being grouped close to the origin of the graph, while those located on the left part of the first dimension are very different.  If considering only the first level of PRP (Figure 3b), three groups of SPs can be differentiated: the blue dotted line separates the SPs belonging to stands artificially regenerated in the past, with ages between 20 and 120 years; the red dashed line separates the SPs where the Norway spruce is the dominant species in terms of AgB (forest habitat 9410); and the continuous green line encloses the SPs belonging to forest habitats 91V0 and 9110. In both cases, the size heterogeneity is negatively correlated with the modified Lorenz curve asymmetry.

Influence of Ecological Factors on Forest Diversity and AgB
To better capture how the tree diversity and the AgB are influenced by site characteristics, a more detailed analysis was performed to complete the PCA output. The species richness showed a negative trend in relation to the mean solar radiation during the vegetation season (Figure 4a), for plots with one to five species, located in the intact forest area. The trend is relevant for the plots with one to five tree species, and a large number of SPs. The negative trend is maintained also for the diversity index ( Figure 4b) and tree size heterogeneity (Figure 4c), although it is less obvious for the latter. However, the differences between the classes that give these tendencies are insignificant.  The AgB shows important variations between solar radiation classes and also within the same class, but without significant differences between classes. Given such amplitudes within the data collected from the sample plots, we could not establish a simple trend of variation between the two variables.
The species diversity is more obviously related to elevation. Both the species number   Most SPs (87%) are located on slopes below 40 degrees. The maximum number of species per plot is found on slopes between 20 and 30 degrees (Figure 5b1). However, the highest species diversity is found on slopes between 50 and 60 degrees (Figure 5b2). Species diversity is relatively low on plateaus or small slopes, due to the proper site conditions for mixed coniferous-beech forests. In such situations, the dominant species are European beech, silver fir and Norway spruce, while the other tree species have low or very low relative frequencies.
On slopes between 40 and 60 degrees, the site conditions lead to more similar relative frequencies of tree individuals, explaining the higher values of D, while over 60 degrees, one species tends to dominate the relative frequency of individuals.
A comparable trend is observed in the case of size heterogeneity of trees (Figure 5b3) but, in this case, the influence of slope is insignificant (K = 7.1892; p = 0.4094). There is also no significant influence of the slope on AgB (Figure 5b4).
The majority of SPs (84) are located on partially shaded slopes, followed by sunny and shaded slopes. The highest number of species is found on shaded slopes (Figure 5c1). No significant differences were found between aspect class and mean species diversity, although a higher D is found on sunny slopes and partially shaded slopes (Figure 5c2). According to our data (Figure 5c3,4), aspect is not related to tree size heterogeneity (K = 6.0704; p = 0.5316) or AgB (K = 6.1636; p = 0.5208).

Tree Species Richness/Size Heterogeneity and AgB Relations
The number of SPs with a certain number of tree species varies according to a unimodal curve (Figure 6a). The mean AgB significantly differs when analyzed on the number of species within a plot (K = 23.6089; p = 0.0003). It increases from 161 to 287 tonnes per hectare between one and three species' plots, then decreases constantly to about 213 tonnes per hectare in the case of six species' plots, probably as a consequence of asymmetric competition for light. Additionally, the differences between AgB means per species diversity classes are significant (K = 32.1567; p < 0.0001), the highest AgB values being recorded for D between 1.41 and 3.80 (Figure 6b). This is further shown by the spatial representation of species diversity, tree size heterogeneity and AgB of living trees on the digital elevation model (Figure 7): the higher values of D are found in the Slătioara UNESCO site and especially its core area (Figure 7a). In terms of size heterogeneity, the core area is also on the top (Figure 7b), even in SPs with lower species diversity (assigned to 9410 N2000 habitat). A high level of species diversity and size heterogeneity can be found even in SPs with low AgB (Figure 7c). The spatial representation of AgB shows the high heterogeneity of the analyzed forests, which is, at least in the case of the core area of the Slătioara UNESCO site, only a consequence of natural processes.

Influence of Forest Diversity on AgB of Living Trees by Habitat Type
The habitat type significantly influences the tree species richness, diversity, tree size heterogeneity and AgB ( Figure 9); in all situations, the values of analyzed indicators are similar for the 9110 and 91V0 habitats and significantly different from the 9410 habitat ( Figure 9). The highest values of biomass and tree diversity are found in the habitat types with beech (9110 and 91V0).  Tree species diversity influences AgB, with the particularities presented below. For habitat type 91V0, the highest means of AgB (about 260 tonnes per hectare) occur when D is between 1.4 and 3.0 (Figure 10a). If considering habitat type 9110, in all SPs the D index values are higher than 1.4 (Figure 10b). A similarly negative trend is recorded for habitat type 9410 (Figure 10c), even though this is reversed for D classes above 2.21.
The relationship between tree size heterogeneity and AgB, analyzed separately per habitat type, showed relatively similar trends as for all 193 SPs (Figure 6c). Thus, in the case of habitat type 91V0 (Figure 10d), the linear positive trend between tree size heterogeneity on Gini index classes and AgB is maintained.
The reduced number of SPs and Gini index classes recorded in the case of the 9110 habitat type have led to non-conclusive results. For this habitat type, the intact temperate forests are characterized by a high size heterogeneity (0.51 < G < 0.80), while the AgB is more stable and less influenced by the tree size heterogeneity (250 < AgB < 315 tonnes per hectare).
In the case of habitat type 9410, the relationship between tree size heterogeneity and AgB is affected by the larger heterogeneity of tree size achieved in natural conditions (Gini index per SP between 0.21 and 0.90) (Figure 10f).

Tree Species/Size Diversity and AgB Relations
Tree species/size diversity was found to be influenced by topographical parameters. In the case of elevation, there is a significant decrease in diversity towards higher altitudes, the diversity reaching a maximum between 800 and 1200 m, with a sudden decline over 1200 m. The pattern is consistent with other findings in the literature, showing that the biodiversityelevation relationship is modulated by soil parameters and available water capacity [82]. Additionally, [52] found that the tree abundance, richness, and phylogenetic diversity first increased with increasing elevation, then reached maximum values at intermediate elevations, and finally decreased at the highest elevations. The reduced diversity values in the lower part of the study site are partially explained by the more recent silvicultural interventions (last one in 2000, in the buffer area of the UNESCO site), which increased the proportion of coniferous [65]. Similar patterns are shown by the relationship between diversity and slope, which is not linear; the maximum number of species being found in the slope range between 20 and 40 degrees, while tree species diversity is maximum at slopes between 60 and 70 degrees. This is explained by the low declivity areas that are occupied by the dominant species (spruce, fir and beech) and the steeper slopes constrain the distribution of trees of different species between adequate microhabitats. Less favorable conditions for forest vegetation were previously proven to heighten the tree diversity and its contribution to the increase in productivity [83]. This pattern can also be explained by the fewer resources on steep slopes due to shallow soils [84]. In such cases, the main species leave room for less demanding species, which compensate less nutrients with more light, due to a different use of canopy space. The effect of aspect was found to be less important in the variation in biodiversity parameters and AgB, which is characteristic for temperate forests [82], but it was shown to be different in more arid environments [56], where aspect controls diversity and density both for mature stands and regeneration.
The size heterogeneity decreases with elevation, with significant differences between core area of UNESCO site (low-intermediate elevations) and the upper part of the study site, where the spruce forest stands tend to be naturally more homogenous [85,86]. The effects of the aspect and slope on size heterogeneity are less visible, since such microclimate indirect variables might be complementary in terms of nutrients [87] or soil moisture availability [44] and, in such cases, it becomes difficult to accurately distinguish the single variable influence on trees' diversity or AgB.

Aboveground Biomass and Topographical Parameters
Knowing that environmental variation is a key regulator of AgB in forests [88], we investigated how solar radiation and topographic heterogeneity in elevation, slope and aspect indirectly influence the AgB. We could not infer a significant effect of the solar radiation on sample plot features, even though there is a significant gradient of radiation. Such relationships are mentioned in the literature especially in tropical forests, where this parameter is used to explain the decrease in productivity and stock along an altitudinal gradient [50].
The fact that the studied forests stretch over a relatively compact area with a limited elevation range did not allow us to check the way solar radiation influences tree size heterogeneity and AgB. Even though there is a slight decrease in tree species diversity with the increase in solar radiation, radiation use efficiency cannot be considered a primary factor to predict productivity and AgB in the study area as, for example, in studies conducted in Norway spruce stands [89]. Additionally, the reduced effect of slope aspect on AgB could be explained by the species behavior (both silver fir and beech are highly shade-tolerant) and the compensation of radiation differences with a satisfactory water regime [82].
The variation in AgB with topographical parameters follows a complex pattern, related also to the type of habitat. In terms of aspect-AgB relationship, our findings are similar with those obtained by [56], who observed that carbon stocks did not differ between aspects.
There is a significant decrease in AgB towards higher altitudes (above 1300 m), reaching very low values above 1500 m, similar to the observations made by [54], who conclude that tree biomass is highest at intermediate elevation in the montane belt. Such variation is also accompanied by shifts in composition towards stands with one or two species in composition, mainly due to limitations determined by lower temperatures during the growing season. The relationship between altitude and AgB could be used to upscale biomass assessments across wider ranges, as elevation is mentioned in the literature as a top contributor to explaining the biomass variation between topographical parameters [90]. Previous studies also explained that the relationship between elevation and AgB is in some cases monotonic [91] or, in other cases, unimodal, with a peak of biomass within certain altitudinal range, in forests without peak consumers [92].

Tree Diversity Influence on Aboveground Biomass
The relationship between AgB and diversity is complex and non-linear. The maximum biomass is found in plots with three species and intermediate diversity, somewhat in contradiction to the findings of [93], who mention the complex relationship between biodiversity and biomass, but with a general positive correlation. The pattern is explained by the relatively high biodiversity of plots found on steep slopes, where the scarcity of water and nutrients reduces the productivity and wood stock.
There is a positive trend between the AgB and tree size heterogeneity observed in the plots established in the core part of the UNESCO site, where the uneven-aged structure of the forest stands (G > 0.5) is decisive for the high biomass stock. Species diversity seems to be a good predictor for AgB due to functional variation brought about by more species that are able to take advantage of complementarity mechanisms such as niche partitioning [83].

Conclusions
Intact forests represent complex ecosystems to be scientifically exploited, as a source of information for studying bio-ecological processes and a potential model for close-to-nature forestry. We compiled data collected from natural forests in order to figure out the complex relationships between diversity, biomass and terrain conditions across different habitats.
The relationship between the aboveground woody biomass and elevation, slope and tree species diversity in intact forests represents a benchmark for managing secondary forests, especially in the context of climate change, when the effects of the environmental factors could be altered by a defective silviculture.
The biomass of highly diverse and undisturbed temperate forests can be considered a target for cultivated forests, aiming to enhance wood provisioning while optimizing species and size diversity. The analyzed relationships highlighted the optimum levels of tree species diversity and trees size heterogeneity that enable modelling the biomass within a well-defined set of conditions.
We revealed the importance of habitat type in galvanizing the effect of topographical parameters on the aboveground biomass of the living trees, as well as the potential species diversity of each habitat type. These relationships are complex and require the separation of the input data on habitat types and classes established for each parameter.
The applicability of our findings in different locations should be handled with caution, considering the site specificity including management background. Collecting similar data from sites located in different vegetation belts would complete the results, allowing the extension of variation patterns identified based on our dataset.
The outcomes open up further research opportunities that could approach the role of microclimate and soil nutrients in deciphering the complex relationship between topographical parameters, tree diversity and aboveground biomass. The following questions could be addressed to further our results: to what extent an effective policy to monitor biodiversity and AgB can rely on the Gini index? Or, even bolder, how forest ecosystems can be steered in order to produce future stands with optimal structural heterogeneity? The second question was partially addressed, but our results can be used to extend the previous findings, in order to identify the optimal values of the Gini index in terms of both high diversity and high biomass.