Quantifying Crown Morphology of Mixed Pine-Oak Forests Using Terrestrial Laser Scanning

: Mixed forests make up the majority of natural forests, and they are conducive to improving the resilience and resistance of forest ecosystems. Moreover, it is in the crown of the trees where the effect of inter- and intra-speciﬁc interaction between them is evident. However, our knowledge of changes in crown morphology caused by density, competition, and mixture of speciﬁc species is still limited. Here, we provide insight on stand structural complexity based on the study of four response crown variables (Maximum Crown Width Height, MCWH; Crown Base Height, CBH; Crown Volume, CV; and Crown Projection Area, CPA) derived from multiple terrestrial laser scans. Data were obtained from six permanent plots in Northern Spain comprising of two widespread species across Europe; Scots pine ( Pinus sylvestris L.) and sessile oak ( Quercus petraea (Matt.) Liebl.). A total of 193 pines and 256 oaks were extracted from the point cloud. Correlation test were conducted ( ρ ≥ 0.9) and ﬁnally eleven independent variables for each target tree were calculated and categorized into size, density, competition and mixture, which was included as a continuous variable. Linear and non-linear multiple regressions were used to ﬁt models to the four crown variables and the best models were selected according to the lowest AIC Index and biological sense. Our results provide evidence for species plasticity to diverse neighborhoods and show complementarity between pines and oaks in mixtures, where pines have higher MCWH and CBH than oaks but lower CV and CPA, contrary to oaks. The species complementarity in crown variables conﬁrm that mixtures can be used to increase above ground structural diversity. petition, and mixture. We derived four response crown variables (Maximum Crown Width, Crown Base Height, Crown projection Area and Crown Volume) from TLS. Our two hy-potheses were conﬁrmed: (1) pines and oaks differ in size when they grow in mixed conditions compared to pure conditions; where pines are taller and thinner, and oaks are shorter and thicker; and (2) the models suggest that in mixtures with pine, oak crowns are shorter and wider while pine crowns are taller and narrower when growing with oaks. The models we developed provide accurate information on species interactions that can be implemented in forest simulators. This information will help forest managers to design more effective silvicultural prescriptions.


Introduction
Trees determine the living conditions for many other organism groups and, typically, they are the most valuable economical component for forests. Therefore, the tree layer is often the compartment of forest ecosystems that is most strongly modified by humans through their management activity [1]. In this regard, much knowledge exists concerning monocultures, which have been largely studied by the scientific community, even though natural forests are typically mixed forests [1] where species' interaction creates more complex forest structures (which make it harder for their dynamics to be understood). A mixed forest is defined as a forest unit of at least 0.5 ha that is composed of no less than two tree species at any developmental stage, sharing common resources as water, light and, soil nutrients [2].
Trees of the same species grow differently in monocultures and mixtures [3] and empirical evidence suggests that species mixing can improve resource utilization within stands [4] due to the complementarity effect, which could be either due to reduced competition or increased facilitation [5]. Tree species with different growth rates and final heights will likely develop structurally more diverse forests than those composed of only one or few species [3] creating multilayered canopies that absorb more light [6]. An increase in tree species diversity can generate a variety of forest structures and interaction between the component species [7], which may lead to more resistant, resilient, and adaptable forests [2,3], and therefore, an important risk-reduction strategy [8]. For this reason, within the last two decades, to understand and manage species' behavior in mixed conditions, many researchers have focused on the performance of mixed forests species input [2,9,10]. Several studies have already shown that mixed stands may have more biomass production than their respective pure stands [8][9][10][11] compared to monocultures [12,13]. Few have studied the mechanisms underlying the mixing effect on productivity [5,10] and confirmed higher rates of primary production in complex forest structures [6]. However, the understanding of tree species interaction in their structure and functioning is still poor [14]. Considering that mixed forests dynamics vary on a small scale [15] more studies are still needed across a variety of forest types to establish a sound theoretical approach across scales.
To understand forest dynamics, Diameter at Breast Height (DBH) and Total Tree Height (TH) are the two most common and easy variables for measuring, analyzing, and modeling forest stands [16]. They are used separately or together in addition with tree species for estimating other important single-tree attributes such as the stem volume, biomass, or the crown [17]. However, it is important to have an accurate measure and knowledge of the crowns because the crown of a tree accomplishes multiple functions [18] and is very heterogeneous in shape despite similar overall tree dimensions [15]. In a stand, crowns are where the effect of inter and intra-specific interaction between trees have been shown [18][19][20]. The crown might allow the trees within a stand to fully utilize limited resources in different spaces and times compared to a monoculture of any species [2]. The tree crown allometry proved to be essential for understanding, appropriate modelling, and silvicultural regulation of mixed stands. For this reason, continuous allometric analyses of various kinds of tree species mixtures need to be done [21].
Historically, the crown has been measured indirectly through empirical established relationships with conventional approaches (mean diameter, dominant height, or volume per hectare), which are relatively easy to make, repeatable, and transferrable. However, they are laborious and imprecise with conventional methods [15] and lead to some problems. For example, they are difficult to validate, difficult to generalize, and uncertainties are poorly quantified or even unknown [22]. In addition, they ignore the three-dimensional nature of stand structure, its most important characteristic ) basing the results on geometric standard crown shapes [15]. The structural complexity increases in mixed forests due to the different growth rates of each species [23], and other processes resulting in spatial niche-partitioning [24] simplified crown forms are not valid anymore to represent the heterogeneity of this kind of forests [15], and more efficient algorithms need to be developed to calculate tree crown variables to facilitate a complete forest resource survey [18]. To overcome these issues, Terrestrial Laser Scanning (TLS) has shown to hold great potential [25]. TLS enables us to quantify the effect of mixed forest in stand dynamics, e.g., [19,[26][27][28][29][30][31][32]. Species identity modifies the mixture outcomes [2], for that reason all possible species mixtures need to be assessed, thus, an adequate management strategy for each species composition can be defined [2].
Within the last decade understanding and predicting influences of inter-specific competition is becoming increasingly important due to an emphasis on mixed-species man-agement [33]. Most studies on mixed-species dynamics were focused on Pinus sylvestris, the species with the biggest area extension in Europe [34,35], together with other central European species like Picia abies or Fagus sylvatica e.g., [29,31,[36][37][38][39]. And few studies were focused on the interaction pine-oak using conventional methods [4,40,41]. Quercus petraea together with Quercus robur, are the two most widespread oak species in Central Europe [42], and are considered helpful for creating more climate-resilient mixed stands thanks to their broad ecological amplitude [42]. In addition, Quercus petraea has a very high environmental and protector value as montane species which, makes it necessary to consider its conservation [43][44][45].
The purpose of this study was to gain a deeper and more accurate insight into the interaction between these two important and widespread species in Europe, Sessile Oak (Quercus petraea) and Scots Pine, (Pinus sylvestris) [46] by analyzing their crown variables through the use of TLS. The basic parameter to define the crown of a tree, according to Lin et al. [18] are crown height, crown volume, and crown projection. In this study, we have analyzed the crown variables Crown Base Height (CBH) and the Maximum Crown Width Height (MCWH) also known as Height of Maximum Crown Extension but MCWH will be used throughout this paper, and as crown size variables the Crown Projection Area (CPA), and the Crown Volume (CV). Analyzing these crown variables, we want to address and determine how the intra and inter-specific competition affects the crown shape of these two species. For this purpose, six pine-oak mixed stands were scanned, and we determined the following research questions: (1) Have pines and oaks different Diameter Breast Height (DBH) and Total Height of the tree (HT) growing in pure vs. mixed stands? (2) Have the explanatory variables of size, density, competition, and mixture a positive or negative relationship in their crown variables (Response variables)?

Study Area and Experimental Setup
Our study was based on a Pinus silvestrys and Quercus petraea forest located in the region of Castilla y Leon in Northern Spain (42 • 54 48 N, 4 • 14 31 W, which is a neighboring region between continental and Atlantic climate. The area is thus characterized by both continental and Atlantic influences in the climate (mean annual temperature of 9.9 • C and a mean annual precipitation 1044 mm). Pines were planted at the beginning of the '70s, and oaks have naturally regrown and no thinning was done in this area in the last 10 years. In 2017, to gain insight into the forest growth dynamic of these two European species in mixed conditions compared to pure conditions, two triplets were established. Each triplet consists of three plots. All plots are located next to each other and have similar site conditions (Figure 1). The Permanent plots set within each triplet are rectangular, limits are marked with wooden poles 50 cm height in each corner. Plots sizes varied to include at least 40 trees of each species of which at least 20 in total are dominant ( Table 1). In the pure pine stands the proportion of pine stems varies from 73.1% to 90%, and in pure oak stands from 85.3% to 95.3%. Finally, in the mixed stands, pines and oaks vary from 41.1% to 45.7%, and from 52.4% to 58.9% respectively. From here, data were collected, and analyzed as shown in the work diagram ( Figure 2).  Table 1. Inventory data of oak-pine triplets where stand indicates plot condition (pure or mixed) and the letters Ps and Qp stand for Pinus sylvestris and Quercus petraea respectively. n represents the total number of trees within each plot; n pines is the number of Pinus sylvestris, n oaks is the number of Quercus petraea, and n other is the number of other species within the plots different from pines and oaks. Pines% and Oaks% represent the percentage of pines and oaks within each plot.

Conventional Measurements
For each tree with a Diameter at Breast Height (DBH) above 7 cm, species were identified, labeled (tree ID) and its position (Cartesian x and y coordinates) recorded with a Total Station Topcon 220. DBH was measured with a caliper (cm); total height and height to crown base were measured in m with Vertex III (Haglöf, Sweden). Finally, tree crown projection radii (in four directions: N, E, S, W) were measured with tape to the closest cm. Table 2 summarizes the stand characteristics for each species. Table 2. Stand characteristics of the study species taken with conventional methods. n stands for number of total tree species. DBH is the Diameter at Breast Height in cm. TH is the total height of the tree in m. CBH is the Crown Base Height in m. CPA is the Crown Projection Area in cm 2 , and BA is the basal area in m 2 /ha.

Terrestrial Laser Scanning
Terrestrial laser scanning (TLS) data were acquired in February and March 2020. Prior to scanning, we georeferenced the plots with a sub-metric GPS Leica model SR20 single frequency equipment with external antenna reception AT501. We recorded 3 corners of each plot for at least 30 min to minimize errors. Each of the points was identified in the field to be easily recognizable through the scanner images. This way, we could overlap point clouds with the map of trees created by the total station Topcon 220 to identify each tree measured with our point cloud files.
TLS data were captured using a Faro Focus 3D device. The laser scanner was mounted on a tripod at approximately 1.3 m above the ground. The settings characteristics are shown in Table 3. To cover all trees belonging to the plots, a pre-design of multiple-scan approach on each tree-plot map with approximately 12 scanner positions per plot was performed. However, the final amount of scanner positions varied depending on stand density to reduce obstruction and assure that tree stems were captured from all sides. To co-register the scans from all the different perspectives to one single point cloud, 15 white plastic spheres 18 cm in diameter placed on one-meter wooden poles were used. The total time needed for scanning for each plot was about 2 h.

Data Processing
The Faro laser scanner creates .fls files. These files were opened and co-registered into one single point cloud file. xyz extension with the software Faro Scene Version 7.0 (Faro Technologies Inc., Lake Marry, FL, USA). Point clouds were imported to the module IMispect from Polyworks Version 12.1.3 software (InnovMetric Software Inc., Quebec, QC, Canada) [19,47,48] together with a .csv file with the UTM coordinates for each tree taken with the Total Station Topcon 220. Thereby, 3D images could be matched with each tree ID.

Tree Segmentation
Tree segmentation consists of obtaining one point-cloud for each tree belonging to the plot, i.e., we should have as many point clouds as trees in the plot (Table 1). For Triplet 1, this isolation was conducted manually: we edited in IMispect the original point cloud selecting each tree and made a copy of the point cloud for each tree. Triplet 2 trees were first isolated by an algorithm developed by [39] which is based on density spatial clustering algorithm with noise (dbscan) function in the dbscan package of R [49] to detect individual-tree positions. Only the x-and y-axes are used as input data. From here, dbscan is automatically able to find each tree as a cluster, i.e., each recognized cluster (stem) now receives a unique number and can be processed individually. Ground points are recognized by dbscan as noise points because of the horizontal structure in contrast to the vertically grown stems. Then each cluster is individually queried as to, which stem base cluster is closest in distance and whether this distance is close enough (≤0.05 m) to be classified as associated points. After this step, each tree is visually checked for completeness. If necessary, unrecognized tree parts are added manually, and artefacts not belonging to the tree are removed using the software RiSCAN PRO. More detailed information is provided in [39].
In both cases, final refined data for each tree were needed and it was performed manually in IMispect by deleting points that were not part of the trees or by separating trees, that were identified as only one due to their crown proximity. All the trees belonging to the study plots were identified but due to canopy occlusion, especially in oak stands, which made tree separation impossible, the study was finally conducted with 91.2% of the total (Table 4).

Tree Metrics
TLS tree metrics were computed using the software "Mathematica 11" (Wolfram Research Inc., Champaign, IL, USA). In total, for each tree, we obtained eleven variables: Diameter at Breast Height (DBH), Total Height (TH), Lean, sweep and the crown metrics of: Crown base height (CBH), Maximum Crown Width Height (MCWH), Maximum Area, Crown Volume (CV), Crown Surface Area (CSA), Crown Length (CL), and asymmetry of the crown with respect to the stem (asymmetry) (Figure 3). An extensive description of the computing process can be found in [32]
TLS data and conventional data were merged by tree ID. We checked the affinity of TLS data with the field data using the Kolmogorov-Smirnov test to test if conventional and TLS data are not significantly different for each studied variable and to prove equivalence and concordance of data, the Lin Concordance Correlation Coefficient [61] was conducted.
From this point, four models from TLS data were analyzed, setting as response variables Maximum Crown Width Height (MCWH) also known as height of maximum crown extension; Crown Base Height (CBH); Crown Projection Area (CPA); and Crown Volume (CV) the explanatory variables were categorized into four groups: size, density, competition, and mixture. Size category was defined by DBH and TH of target trees. To define density, we used Total Basal Area (BA total ) of each tree, understanding density as competition at stand level. On the other hand, the competition was defined at tree level by the following variables: Largest Basal Area (BAL), both by species separately and total, Asymmetry [32] and the Hegyi Index [62], hereafter referred to as the Competition Index (CI). Finally, to define the mixture, we used the ratio variables of BA, BAL, and the total number of pine trees.
Density, competition, and mixture were defined for every single tree following [29]. The three radii sizes (5, 7.5, and 10 m) were chosen for two reasons: (1) our study plots are maximum 30 × 30 m; and (2) the influence radii for each tree is assumed to be between 1.5 and 2.5 times the target tree crown width (maximum in our sample is 7.41 m). Thanks to this, these three categories were not analyzed as categorical variables distinguishing between pure and mixed stands, but rather as continuous variables.

Explanatory Models
Based on forest modelling literature, different linear and non-linear models were tested to explain response crown variables (Table 5), adapting them to our four category explanatory variables (size, density, competition and mixture). Table 5. Explanatory models defined for this study. s = size, d = density, c = competition and m = mixture, α 0 = intercept, α 1−4 = parameters of each coefficient (s, d, c, and m respectively), ε i = error.
Before fitting models, correlations between explanatory variables were calculated and those that were highly correlated (ρ ≥ 0.9) were removed (oak ratio variables), to avoid multicollinearity problems. Thus, we used pine as a reference for the basal area, the number of trees per plot, and the BAL ratios. In addition, for CPA and CV response variables, the explanatory variable d 2 h (squared diameter at breast height times height) was included as a size variable since this variable represents a proxy for tree volume. Furthermore, for CPA response, the logarithmic transformation of Basal Area was included as a density explanatory variable, as in [63] For each radius of influence considered (5, 7.5, and 10 m) all possible combinations of explanatory models (size, density, competition, and mixture) were tested. We unified all models performed with all influence radii for each response variable, with all resulting models ordered from the lowest to the highest Akaike coefficient (AIC). According to [64] the top five models with the lowest AIC index were selected. Finally, the best model for each response variable and species was selected if it had a coherent biological sense.

Data Analyst
The Kolmogorov Smirnov test was not significant (α = 0.05) and the Lin Concordance Correlation Coefficients for DBH and TH were 0.98 and 0.74, respectively, confirming the distribution of TLS data corresponding to the data taken in the field are consistent and have no significant differences. From then on, the following analysis was performed only with TLS data.
Firstly, descriptive statistics of all trees classified by species and pure and mix plots was performed (Table 6). We observed that oaks show similar tree heights, both in pure and mixed stands, but they are slightly thicker (DBH one cm higher) in the mix with pines than in pure conditions. Unlike pines, DBH is up to almost 3.5 cm thinner in the mix with pines but slightly taller (+0.5 m) compared to pure plots.

Explanatory Models
The best models selected according to the lowest Akaike index (AIC) and biological sense are shown in Table 7 and as equations in Table 8. For all cases, the best model selected was the model with the lowest AIC index, except for the MCWH variable in both species, where models with the second AIC index were selected due to the first ones having no biological sense. Table 7. The explanatory models selected as the best fit according to their lowest Akaike index and biological criteria for each response variable (Variable) and species. r = radius of influence (5, 7.5 and 10 m); s = size; d = density; c = competition; m = mixture; α 0 = intercept; α 1 -4 = the coefficient numbers for each explanatory variable (size, density, competition, and mixture); AIC = Akaike; K-S test = p-value of Kolmogorov-Smirnov test for residuals; R 2 = coefficient of determination of the model.  Table 8. Final equations for each species and crown response variable.

MCWH
Pine We observed that the size of the tree is a significant variable for all cases, with only one exception, MCWH for pines. In all cases, tree size affected the crown shape the bigger the tree, the greater increase in growth of the CBH, MCWH, and their CPA and CV as well. That is to say, the bigger the tree, the higher and wider the crown will be. Competition boosted the MCWH and CPA of pines, but it is we observed for the rest of the variables competition has a negative effect on crown shapes for both species, as crowns were shorter and narrower.
The mixture variable is always significant for oaks. The presence of pines makes the height of the oaks' crowns (MCWH and CBH) smaller, and by contrast, makes their crown projection and volume larger.

Discussion
In this study, we have obtained a more comprehensive understanding of pine-oaks forest dynamics thanks to the use of TLS. We have scanned six plots (two triplets) formed by Quercus petraea and Pinus sylvestris in pure (Two plots for each species) and mixed stands (Two plots). With the scanning four-crown variables: MCWH, CBH, CPA, and CV from a total of 193 pines and 256 oaks and later analyzed them by linear and non-linear regression using tree size, stand density, competition, and mixture proportion as explanatory variables (Table 5). Several models were fitted for each crown response variable, these models were selected by the lowest AIC index and biological sense (one for each species and variable). Results ( Table 8) showed us that tree crowns occupy the gaps between species, and that mixture, competition, density and size play a significant role in the shape of the species' crowns. Oak crowns remain under pines but are wider, while pine crowns are higher and narrower when interacting with oaks.

Mixing Affects Tree Size
Similar to [65] we observed differences in height and diameter of tree species growing in mixed stands vs. monocultures. We have found slightly higher TH for Scots Pines in mixed than in pure stands but no such pattern for their DBH. On the contrary, Sessile Oaks showed larger DBH in mixed stands, but their TH was slightly lower. This is to say that when in interaction, oaks are shorter but thicker and pines are thinner but taller. This explains a complementarity between pines and oaks within the mixed stands which could be related to specific growth strategies [35] when growing shade-tolerant species together with light-demanding species [66]; or to a more efficient use of resources over time [67], as interspecific competition for light decreases in the mixture [68] enabling changes in tree architecture and crown plasticity, and therefore, increase light capture and productivity [6,69]. Shade-tolerant oak trees utilize resources more efficiently in the mixture which could result even in higher radial increment rates [70]. Regarding this aspect, [66] found the higher productivity in Quercus petraea mixed stands when they grow together with light-demanding pioneer species with deeper roots features [71] as it is Pinus sylvestris.

Effect of Size, Density, Competition and Mixture on Species Crown Size
Our results suggest there are both competition and mixture effects on oaks and pines' crown dimensions and shape (Table 8) suggesting that species mixing modifies the crown size and shape, and thereby the canopy space-filling [1]. Some researchers have hypothesized differences in crown shape between deciduous broadleaves and evergreen conifers are at the origin of a positive mixture effect on species productivity in broadleavedconifer mixed stands [72]. Here, the models selected showed us there is a significant effect of size, density, competition, and mixture variables, either negative or positive on crown allometry. [66] highlighted that there exists either a competitive advantage of Q. petraea over a more light-demanding species, like Pinus sylvestris, or a result of the complementary use of resources in the mixture, as it is shown in Figure 4 where it is observed how at a same tree size for both species under two density scenarios (based on BA total ), oak CPA vary less than pines. Regarding mixture, explanatory models revealed that the presence of pines positively affected the crown expansion of oaks (expressed as CV and CPA). At the same time, crown elongation (expressed by CBH and MCWH) was negatively affected (Table 8). On the contrary, pines seem to have higher crowns when growing together with oaks or with bigger pines in the neighborhood (represented by the trees in the three radii considered). Pine crown expansion (CV and CPA) was not significantly affected by the mixture. Our findings suggest that crown oaks, in the presence of pines, remain under the crown pines, occupying the space between pine stems, that may be the reason why crown oaks are shorter but wider, contrary to pines where crowns start at a higher point when growing together with oaks.
Our results suggest a multi-layered canopy (stem diameter and height variety) in oaks-pines stands in Northern Spain. We hypothesize that multi-layered canopies are produced by the mixture complementarity of crown shapes of pines and oaks due to their differential crown architecture. On the one hand, shade-tolerant and shade-intolerant tree species combination [23,73], together with slower growth of oaks lead to differentiation in the canopy. Shade-tolerant species tend to have a crown shape optimized for the capture of light under limiting conditions [74], and in mixtures force non-tolerant species to grow to reach upper canopy level. Additionally, resource use complementarity appears when in winter the lack of leaves in oak allows pine trees to capture light and photosynthesize.
However, the mixture also limits light resources in summer creating interspecific competition, as we have seen in our analysis. This fact may cause changes in tree species crown allometry [66]. Our results suggest a complementarity in canopy space occupation; Oaks produce wider crowns and pines larger stems to be able to capture the light above oaks. This can be explained by differences in foliage persistence during the year [2], and by inter-specific variations in crown architecture and height, when combining species with different shade tolerances or vertically-oriented species with more laterally expanding tree species [5,12].
Finally, the origin of the forest, i.e., plantation or natural regeneration may affect trees and associated vegetation relative to those in monocultures [65], the same way we are seeing in our study combining pine plantation with natural oak resprouting where we have found trees' neighborhood seems to be affecting crown allometry.
This apparent rapport between oaks growing together with pines may agree with relaxing resources competition in mixed forests that leads trees species to temporal diversification and spatial niche partitioning as suggested by [69] and [35]. Wider tree spacing reduced competition and increased resource acquisition capacity per tree, therefore facilitating crown expansion [31] as our models suggest for oaks.
Contrary [23] we have not found the density of shade-intolerant tree species to be the most significant explanatory variable in multiple regression. We found density only explains the CPA and CV of the trees. For these two variables, in line with [75] we saw that trees in stands with a higher BA exhibited greater CV and CPA for both oak and pine, creating a more heterogeneous stand structure.
This study is the most comprehensive work focused on the crown structure on Quercus petraea-Pinus sylvestris mixture contributing to the understanding of this very important mixture. Using TLS, we were able to focus on the crown structural pattern in great detail, in order to gain deeper insight into tree crowns' response to pines and oaks mixture. We were able to quantify geometrical aspects of tree crown: Maximum Crown Width Height, Crown Base Height, Crown Projection Area, and Crown Volume that were conducted on a large number of trees (449 trees, 193 Scots Pines, and 256 Sessile Oak). Therefore, we have obtained rigorous information about crown differences between these two species and how their size, density, competition, and mixture around affect their crown shape, obtaining a plausible answer hidden behind the identified productivity changes between species as [13] stated for pine-pine mixtures. As discussed by [20], the characterization of crown species gives us information on how species tend to occupy the stand canopy and helps to quantify their plasticity. For instance, an increase in oak CPA and CV forces Scots pine trees nearby to raise their crowns to reach upper canopy levels. Unlike other deciduous species, which have a remarkable capacity to adjust their morphology and physiology to a particular set of light conditions [76].
The fitted models of this study suggest mixture increases stand structural complexity index as pointed in [31] due to better use of the available light [26] when growing shadetolerant and intolerant species together. Our study plots had not been thinned in the last 10 years but previous silviculture could have caused a residual effect in our findings due to the memory effect found in the forest stand dynamic [16,77]. The origin of the forest, i.e., plantation or natural regeneration may affect trees and associated vegetation relative to those in monocultures [65]. In our study, combining pine plantation with natural oak resprouting where we have found trees' neighborhood seems to be affecting crown allometry.
Long-term observation of mixed forests and general findings are still rare [8] especially when involving individual tree analysis [78]. Using detailed quantification of the crown can help to gain insight in differences in pure and mixed forests as a result of species' interactions and to define sound management options.

Conclusions
Crown plasticity represents one of the species potential to acquire resources and occupy space [21]. In this study we investigated how the crown shapes of two widespread species in Europe, Pinus sylvestris and Quercus petraea, varied with tree size, density, com-petition, and mixture. We derived four response crown variables (Maximum Crown Width, Crown Base Height, Crown projection Area and Crown Volume) from TLS. Our two hypotheses were confirmed: (1) pines and oaks differ in size when they grow in mixed conditions compared to pure conditions; where pines are taller and thinner, and oaks are shorter and thicker; and (2) the models suggest that in mixtures with pine, oak crowns are shorter and wider while pine crowns are taller and narrower when growing with oaks. The models we developed provide accurate information on species interactions that can be implemented in forest simulators. This information will help forest managers to design more effective silvicultural prescriptions. Funding: This work has financially supported by CARE4C H2020-MSCA-RISE-2017, Grant Agreement no. 778322 and by SMART VA183P20 -Bosques mixtos: Selvicultura, Mitigación, Adaptación, Resiliencia y Trade-offs (project code VA183P20) funded by the Castilla and León Regional Government (Spain) and the the European Regional Development Fund (ERDF).

Informed Consent Statement: Not applicable.
Data Availability Statement: Data sharing not applicable.