The social drift of trees. Consequence for growth trend detection, stand dynamics, and silviculture

Recently, many studies worldwide tapped tree ring pattern for detection of growth events and trends caused by weather extremes and climate change. As long-term experiments with permanent survey of all trees are rare, growth trend analyses are mostly based on retrospective measurements of growth via increment coring or stem analyses of the remaining individual trees in older forest stands. However, the growth of the survivor trees in older stands may only unsufficiently represent the course of growth of the dominant trees throughout the stand development. Here, the more than 100 years survey data of the European beech (Fagus sylvatica L.) thinning experiment Fabrikschleichach in South Germany are used to show the long-term changes in social ranking of trees and their consequences for growth trend detection by retrospective tree ring analyses, for stand dynamics and silvicultural management. Firstly, a significant social upwards drift of initially medium-sized trees till 2010 is shown based on the trees' percentiles in the stem diameter distribution in 1904 versus 2010. The social climbing is stronger on the thinned compared to the unthinned plots. Secondly, we show that 40–60% of the 100 tallest trees in 1904 were replaced by social climbers and down-ranked below the 100 tallest trees till 2010. Linear mixed model analyses reveal that the long-term trend of the diameter growth of the 100 dominant survivors in 2010 was on average by 23% steeper than the trend of the 100 tallest starters in 1904. This indicates that the survivors had a steeper and longer lasting growth than the originally dominant trees. Thirdly, the diameter growth trend in the last 20 years, from 1990 to 2010, is analyzed in dependency on the current and past social position. A linear model shows that early subdominance or suppression can significantly steepen the growth trend a century later and vice versa. Finally, we discuss the implications of the social drift for the survivor-based growth trend analyses, for the stand dynamics, and silvicultural management.


Introduction
Trees in seasonal forests record their autobiography in the annual ring width pattern. This autobiography can be harnessed for biomonitoring and bioindication (Markert et al. 2003). The annual or periodical stem growth is an integrative signal of a tree's growth and vitality (Dobbertin 2005). Growth rates primarily represent an unspecific tree response. However, given the impact of silvicultural measures, pathogen attacks, or environmental changes on response patterns, growth trajectories may both indicate and serve to quantify the effect of various external episodic or long-term disturbances of tree growth. Accordingly, a growing body of research worldwide is tapping tree ring chronology as an autobiographical measure to detect growth events and trends, e.g., caused by weather extremes and climate change (Borgaonkar et al. 2011;Camarero et al. 2018;Girardin et al. 2016;Pretzsch 2020a, b).
The detection of growth disturbances improves with continuity of measurement on long-term experimental plots, e.g., using repeated surveys or permanent dendrometer records (Nehrbass-Ahles et al. 2014;Prietzel et al. 2020;Spiecker et al. 2000). Thus, time series data should ideally be obtained over the entire tree lifetime of 100-200 years. However, such data are very rare (Nagel et al. 2012;Pretzsch et al. 2019). Instead, shorter measurement periods or retrospective measurements of growth via increment coring (Camarero et al. 2018;del Río et al. 2020;Pretzsch et al. 2020a andb, Spiecker et al. 2002) or stem analyses (Guttenberg 1915;Pretzsch 2020c;Spiecker et al. 2012;Wichmann 2001) of the remaining individual trees in older forest stands to trace back their growth as a substitute of permanent longterm survey are more common. This is a problem, because the growth of the survivor trees in older stands may not accurately represent the course of growth of the dominant trees throughout the stand development. Indeed, trees remaining in the advanced stand development phase may represent a characteristic latecomer's growth trajectory with low growth in the juvenile but accelerated growth in the adult development phase. Moreover, the growth trend of the survivors may not indicate any climate change induced growth trend but result from the natural social differentiation of the individuals within forest stands and represent the natural growth variation of a subset of latecomer trees (Brienen et al. 2017;Nehrbass-Ahles et al. 2014). Thus, using the surviving trees to make inferences about the long-term growth trend of the mean course of dominant trees, of all trees or the stand in total may be misleading.
The taller the initial tree size compared with neighbors (e.g., quantified by the diameter percentile), the better the primary individual growth potential and the perspective of a tree. Size growth is furthermore intrinsically determined by size (Bertalanffy von 1951;Zeide 1993). However, genetic determination, biotic and abiotic disturbances, and silvicultural treatment may further influence the development of the social position of individual trees, i.e., whether they maintain, improve or lose their relative social ranking that finally determines access to light and other resources. Many studies show that the tree size growth trajectory follows an S-shaped curve depending on age and size; and that quickstarters peak and decelerate earlier, while slow starting trees may accelerate in advanced ages (Guttenberg 1915;Martín-Benito et al. 2008;Smith et al. 2019). However, recent studies show that asymptotic bending of size growth occurs much later than assumed so far (Pretzsch 2020c, Sillett et al. 2015aand b, Stephenson et al. 2014). Nevertheless, initially smaller or suppressed trees in a stand may begin slower but continue to grow longer compared to taller trees of the same age (Assmann 1970;Magin 1959;Wenk et al. 1990).
Social differences between trees and social drift are caused by varying growth trajectories and result in structural heterogeneity. They may affect many ecological, economical and socio-economical forest functions and services such as habitat structure, wood quality, growth stability, recreation and protection functions against windthrow, erosion, or avalanches (Dieler et al. 2017;del Río et al. 2016). Moreover, social heterogeneity and drift over time may generate a continuous structural differentiation where late social climbers may replace the position of social 'losers' and thereby maintain a continuity of structure and function over time. This stands in juxtaposition to clone plantations, where the social ranking and variety of trajectories and timepoints of optimal harvest are kept very similar.
The question to which degree growth trend analyses based on the sampling of survivors are biased and systematically overestimated was raised among others by Cherubini et al. (1998) andNehrbass-Ahles et al. (2014). The authors proposed sampling methods to reduce or avoid the bias. Brienen et al. (2017) further showed that tree demography can generally dominate growth trends inferred from tree ring analyses. Also, Duchesne et al. (2019) reported growth increases in boreal forests as artefacts stemming from sampling biases. Consequently, the study at hand aims to explain and quantify the survivor bias in monospecific and even-aged European beech (Fagus sylvatica L.) stands that are currently in the focus of growth trend analyses (Ammer et al. 2005;Bosela et al. 2016;Geßler et al. 2007;Pretzsch et al. 2020a, b;Rennenberg et al. 2004). Several studies reported a diameter growth increase in dominant trees by 30-40% in the last 100 years, which was mainly ascribed to the extension of the growing season, warming, increase in nitrogen deposition and the atmospheric CO 2 −concentration (Dulamsuren et al. 2017;Hilmers et al. 2019;Pretzsch et al. 2020c). In our study, those revealed growth trends will be contrasted with the magnitude of the potential survivor bias to avoid overestimation of growth trends derived from increment coring of mature dominating European beeches as practiced, i.e. as a substitute for bioindication and trend analyses where longterm experiments are lacking.
To examine this, we used data acquired from the approximately 200-year-old thinning experiment Fabrikschleichach 15 of European beech in South Germany, which contains unthinned (A-grade), and moderately (B-grade) and heavily (C-grade) thinned plots. We will employ these repeated surveys of the social status, growth, and tree morphology from 1904 to 2010 to analyze long-term changes in social ranking and respective growth trends.
First, we will address how social ranking (in terms of stem diameter percentiles) has changed from the first to the present survey (~ 100 years later) and how the collective of the 100 tallest trees changed over this time period. Second, the mean stem diameter growth of the 100 tallest trees in the beginning will be compared to the growth of the 100 tallest survivors at the end of the survey period. Third, we will examine how the growth trend of the survivors over the last 20 years has depended on their past and present social status. The questions were addressed for a 106-year period, for untreated, moderately, and heavily thinned stands, and additional crown measurements were used to characterize the morphological attributes of quickstarters compared with latecomers.
Specifically, we focused on the following hypotheses.

H1
There is a perfect correlation between the social ranking of trees at the beginning and end of the 106-year survey period, i.e., the ranking is stationary.

H2
There are no differences between the course of stem growth of the 100 tallest trees at the beginning and end of the survey period; i.e., there are no social climbers and losers within this dominant and co-dominant collective.

H3
The social ranking at the beginning of the 106-year survey period has no effect on the growth trend of the trees at the end of the survey period, i.e., the growth trend of the survivors is independent from their past social status.
Finally, we will discuss the implications of the social drift of trees in the context of growth trend analyses, forest stand dynamics, and silvicultural management.

Material
The long-term thinning experiment Fabrikschleichach 15 (FAB 15) in European beech was chosen for this study as it represents one of the longest continuously measured forest experimental plots worldwide. It is an ideal example for questioning methodological aspects of the above introduced growth trend analyses in European beech stands. The about 0.4 ha sized unthinned, moderately thinned, and heavily thinned plots (A-, B-, and C-grade plots) of FAB 15 were established in 1870/1871 in a 48-year-old, even-aged European beech stand that originated from natural regeneration by shelterwood cutting (Burschel and Huss 1987, p. 90-92) of a beech stand in 1822. Since their establishment, the plots have been measured 19 times. FAB 15 is located 1 km east of the village Fabrikschleichach and belongs to the Ecoregion 5.2 Steigerwald and to the forest enterprise Ebrach. The location is 10.62°E longitude and 49.95°N latitude, the elevation 470 m a.s.l., the annual precipitation 820 mm (420 mm in the growing season), the mean annual temperature 7.5 °C (14.5 °C in the growing season), and the length of the growing season 150 days. The natural vegetation in this region would be submontane European beech-sessile oak forests. Soils are mainly Loamic Luvisols according to Blume and Schad (2015) with Semionoten sandstone as basic material.
FAB 15 comprises A-, B-, and C-grade thinning plots. Most of the classic long-term experiments comprise socalled A-grade plots (VDFV 1902). "A-grade" is defined by VDFV (1902, § 4) as follows "This is limited to the removal of dying and dead trees, as well as any bowed pole wood. […]. for the purpose of delivering material for comparative growth investigations only". In other words, on A-grade plots nothing more is done than closely monitoring natural mortality and removal of dying or dead trees to prevent possible stand damage coming from dead trees (infestations by fungi or insects). A grade-plots reveal the maximum stand density and self-thinning, and they serve as the reference for quantifying how different levels of stand density regulation effect productivity, carbon sinks and stand structure.
As B-and C-grade thinning moderately and heavily reduce stand density mainly by removing mainly small trees and only tall trees when they are of bad quality, they are called moderate and heavy thinning from below (VDFV 1902, § 4). The B-and C-grade plots were thinned 11 times since their establishment. The stand basal area of the remaining stands ranged between 28.6 and 58.4 m 2 ha −1 on the A-grade plot, 26.6 and 51.0 m 2 ha −1 on the B-grade plot, and 20.2 and 42.8 m 2 ha −1 on the C-grade plot. This decrease in the stand basal area from the A-to the C-grade plot reflects the increasing strength of thinning from below. From the first survey at age 48 years in spring 1871 to the fifth survey at age 77 in autumn 1899, the tree measurements were registered in tally-sheets. The measurements included stem diameter at breast height, tree height, and height to crown base. In addition, the removal trees were recorded. In the course of the sixth survey in 1904 at a stand age of 82 years the trees were numbered individually. Since that time the individual tree development in terms of tree diameter, tree height and height to crown base is known and can be analyzed. This is the reason why the stand characteristics can be reported since 1871 and the individual tree characteristics not before 1904.
The repeated measurements covered in addition to diameter at breast height ( d 1.3 , cm) and tree height ( h , m), the height to crown base ( hcb , m) and at several surveys the crown radii of the trees in eight cardinal directions according to standards described by Pretzsch (2009, pp. 115-118). The eight radii were used to calculate the crown projection area 1 3 in m 2 ( cpa = cr 2 × with cr 2 = √ r 2 1 + r 2 2 + … + r 2 8 ∕8 ). The survey in 1981 included the measurement of the tree coordinates of all trees on the plots and the border trees outside the plots but with crowns reaching into the plots.
The repeated surveys resulted in the data suitable for calculation of all common stand characteristics for each of the 19 successive survey periods since 1870 (Biber 2013). Table 1 gives an overview of the stand characteristics at the last survey (see Pretzsch 2009, pp. 208-210 andPretzsch 2019, p. 309 for overview tables with all relevant stand characteristics). The reported stand level data were derived from the diameter and height records as well as from the removal trees. We used standard evaluation methods according to standard procedures recommended by DVFFA (1988). With the DESER-Norm (Johann 1993), the German Association of Forest Research Institutes, DVFFA, (in German: "Deutscher Verband Forstlicher Forschungsanstalten") has made an important step towards such a standardization. For estimating the merchantable tree volume of European beech, we used the approach by Franz et al (1973) with the equations and coefficients published by Pretzsch (2002, p. 170, Table 7.3).
The mean height at age 100 presented in Table 1 indicates a site index of I.-II. according to the yield of Schober (1967Schober ( , 1975, with other words good-very good site quality and productivity. This is corroborated by the main stand characteristics of the three plots at the last survey (Table 1). The current annual volume growth in Table 1 refers to the last survey period from 2000 to 2010. The mean annual volume growth, the total volume yield, the removed volume, and the intermediate yield in Table 1 refer to the stand development since its establishment in 1870/1871. The stand density, represented by the tree numbers, N, stand basal area, BA and standing stem volume (merchantable wood > 7 cm at the smaller end), V, decreases continuously from the A-to the C-grade due to the thinning interventions. This stand density reduction means more growing space for the remaining trees and had a strong positive effect on the stem diameter, d. The current and mean annual stand volume growth, and the total yield were increased by the thinning on the long term. The last two columns in Table 1 show that the removed volume since stand establishment is about double as high on the C-grade compared with the removed volume on the A-grade plot (trees that were dead and dying due to self-thinning). The percentage of removed stand volume in relation to the total yield at age 188 were increased by the active thinning from 21.2 (A-grade) to 38.2% (C-grade). Table 2 shows that the stem dimensions d 1.3 , h, and v, at the last survey were much higher in moderately and heavily thinned stands than on A-grade plots. Even stronger was the effect of thinning on the crown dimensions, cr, cl, and cpa. The increased size on the B-and C-grade plots is associated with higher stem diameter and stem volume growth in the last survey period. Regarding all variables, the effect size Table 1 Stand characteristics of the A-, B-, and C-grade plots at the last survey in 2010 at the stand age of 188 years N, tree number per hectare; d mean stem diameter; h, mean tree height; BA stand basal area; Current annual volume growth, stand volume growth in terms of merchantable wood (> 7 cm at the smaller end) in the last survey period 2000-2010; Mean annual volume growth, average stand volume growth per year since stand establishment; Total yield, total stand volume growth since stand establishment; Removed volume, volume removed since stand establishment (dead trees on A-grade plot; thinned trees on B-and C-grade plots); Intermediate yield, removed stem volume/total yield ×100 . The Intermediate yield indicates how much of the total yield were already removed by preliminary harvest since stand establishment until age 188 of thinning was much stronger for the transition from A-to B-grade than from B-to C-grade. For more information about the long-term stand development of all three plots see Supplement Stand Description of FAB 15 in the Electronic Supplement. For more details about this long-term experiment see works by Assmann (1950Assmann ( , 1970 and Kennel (1972) regarding thinning reactions, by Franz et al. (1993) and Pretzsch et al. (2019) regarding long-term survey and monitoring and by (Pretzsch 1999, 2014a andb) dealing with growth trend analyses on FAB 15.

Methods and models
H1 There is a perfect correlation between the social ranking of trees at the beginning and end of the 106-year survey period.
For a descriptive view on the dependency between social position and tree growth, all individual tree diameter growth trajectories from 1904 to the last survey were modelled by straight lines ( d = a 0 + a 1 × age ) via linear regressions individually. The relationship d−age was represented by a linear model according to several recent studies on growth of European beech (Pretzsch et al. 2020a-c). Furthermore, based on the stem diameters of all trees of the first individual survey in 1904 the diameter percentiles (from 0 = lowest to 100 = predominant) were calculated. The plot of the slopes of the diameter growth courses over the percentile in 1904 allowed an overview of the nonlinear relationship between social position and course of tree growth for all three thinning grades.
In order to scrutinize any social drift with progressing stand development (H1), the stem diameter percentiles were also calculated for the last survey in 2010. The linear regression was fitted to the stem diameter percentiles of the A-, B-, and C-grade plots separately and also to the data of all three plots together. The resulting slopes a 1 ± 3 t SE indicated whether the slopes deviated significantly (p < 0.001) from a perfect concordance between the diameter ranking at the beginning and end of the survey period. For a strict scrutiny, we chose a 99.9% confidence limits and corresponding tail probabilities of p < 0.001. A stationarity of the slopes would be indicated by slopes of a 1 = 1.0.
H2 There are no differences between the course of stem growth of the 100 tallest trees at the beginning and end of the survey period.
(1) perc 2010,k = a 0 + a 1 × perc 1904,k + k In order to quantify the extent and effect of social drift, the 100 tallest trees per hectare in 1904 were compared with the 100 tallest in 2010. The number of social climbers resulted as the number of trees that were among the 100 tallest in 2010 but not in 1904. This number was divided by 100 and provided the percentage of social 'climbers'. Analogously, the trees available among the 100 tallest in 1904 but no longer in 2010 provided the social 'losers'.
Furthermore, in order to compare the prospective and retrospective view on the collective of dominant trees the course of the diameter growth of the 100 trees that were the tallest at the beginning in 1904 (quickstarters) was compared with the course of the diameter growth of those trees that were the 100 tallest at the last survey in 2010 (latecomers, also called survivors in the following text). Throughout the following text we use the term quickstarters for trees that were among the 100 tallest trees at the beginning of the stand development. We will use the term latecomers or survivors for the trees that belonged to the 100 tallest trees at the end of the 188-year survey period. The growth trajectories of both quickstarters and latecomers were visualized and analyzed in detail.
In addition to their course of diameter growth the mean stem and crown characteristics of the four mentioned collectives, i.e., the quickstarters, survivors, climbers, losers were quantified and the attributes of the starters were tested (via t-test) against the climbers, the survivors against the losers.
In order to reveal any differences in their development, the collectives of the quickstarters and latecomers were analyzed by linear mixed effect models as follows. Models 2 and 3 were applied to analyze how the membership of a tree to the group of quickstarters (group 0) and latecomers (group 1) determined the intercept and slope of the stem diameter-tree age relationship. The fixed effect variables are tree age and group; latter is categorial with group = 0 for the quickstarters and group = 1 for the latecomers. The models in general can be seen as a typical linear relation between tree size in terms of stem diameter and group effects on both, the intercept and the slope. The fixed effects parameters are a 0 − a 3 . The indexes i, j, k in the following equations refer to the levels of experimental plot, tree within plot, and single observation, respectively. The random effects b i and b ij cover the levels of plot and tree.
Model 2 was applied for the analyses at the A-, B-, and C-grade plot levels separately. In this case the random effect variable b i accounts for the temporal autocorrelation (several successive surveys) at the tree level Model 3 was applied for the plot overarching analysis (pooled data of A-, B-, and C-grade). In this model, two random effect variables account for the hierarchical nature of the data. Random variable b i accounts for the spatial correlation (several plots per experiment) between the plots and b ij accounts for the temporal autocorrelation (several successive surveys) at the tree level.
The Models 2 and 3 were applied for stem diameter development of various tree groups depending on tree age, calendar year, and interaction between tree age and calendar year. In order to detect the respective growth behavior on the different plots first and derive a generalized growth pattern of the quickstarters and latecomers in a second step, I refrained from combining models 2 and 3.

H3
The social ranking at the beginning of the 106-year survey period has no effect on the growth trend of the trees at the end of the survey period. The numbers of the Models 1-4 refer to the results in the text and in the captions of the tables.
The statistical software R 3.4.1 was used for all calculations, in particular the function lme from the package nlme (Pinheiro et al. 2018).

Overview of the stand and tree development
The periodical annual stem volume increment, PAI, per hectare and year on the three plots, visualized on Fig. 1, shows the following relevant aspects for this study. Throughout the whole survey period the stand growth oscillates quite strongly but remains on a considerable high level. There is no clear decreasing age trend in the PAI development, slope 1990−2010,k = a 0 + a 1 × perc 1904,k + a 2 × perc 2010,k + k although the stands were already 188 years old at the last survey. The PAI of the A-grade develops at the lowest level, the PAI on the B-, and C-grade plot is in most periods higher than on the A-grade. On the long term the thinning did not reduce the PAI and did not cause a preponed reduction or the PAI development over age. The long continuation of stand growth, the high level beyond the age of 150 years and the high level of growth despite of the strong thinning indicate considerable positive deviations from the development predicted by the common yield tables for European beech by Schober (1967Schober ( , 1975.
The courses of the stem diameter growth (Fig. 2a-c) are an essential basis for the subsequent analysis. The black trajectories represent the remaining trees in the stands; the red trajectories visualize the trees that dropped out by selfthinning on the A-grade plots and were removed by moderate and heavy thinning from below on the B-and C-grade plots. Comparison between the trajectories on the thinned and unthinned plots (Fig. 2b and c versus a) evinces that thinning accelerated the diameter growth. The range of stem diameters at the last survey is about 20-70 cm on the A-grade and 40-90 cm on the C-grade plots. Interestingly, the courses of annual stem diameter growth show hardly any reduction although the trees aged to nearly 200 years till the end of the shown period 1904-2010 (Supplement Fig. 1). The courses are rather linear and intersect in many cases; i.e., the size ranking changed with progressing stand development. The average trend is more or less linearly increasing in case of the diameter development (Fig. 2) and cubic in case of the stem volume development (Fig. 3). Interestingly, the linear trend of the stem diameter development and cubic trend of the stem volume proceeded from middle to advanced tree ages, and the trend is similar for all three thinning grades. The stem volume developments hardly showed any S-shaped trajectories (Fig. 3). This means that the annual individual trees' volume growth increased in most cases continuously until high tree ages (Supplement Fig. 2) and did not yet show a clear age-related reduction as assumed by textbooks by among others Kramer (1988,   A-grade, without thinning and only removal of dead or dying trees; B-grade, moderate thinning from below; C-grade, heavy thinning from below 1 3 p. 66), Bruce andSchumacher (1950, p. 377), Wenk et al. (1990, p. 74), Assmann (1961, p. 80), or Mitscherlich (1970. The red trajectory above the rest of the trees in Figs. 2c and 3c (starting above the rest of the population and ending after the removal of the tree in the 1950ies) represents a hold-over tree from the previous stand generation. This tree was eliminated from the analysis.
Correlations between the social ranking of the trees at the beginning and end of the 106-year survey period (H1) Figure 4 shows the individual trees' slopes ( slope = a 1 ) of their diameter growth from 1904 to 2010 (based on fit by linear regression d = a 0 + a 1 × age ) plotted over the trees' percentile at the first individual survey in 1904. These presented slopes come from individually fitted regression lines; so they do not include random effects. The small, subdominant trees with low percentiles in 1904 have low slopes of diameter growth and preferentially drop out during stand development on the A-grade plots (red dots, Fig. 4a). On the B-and C-grade plots the red dots represent the trees removed according to the guidelines for thinned from below; small trees removed due to their subdominant position and tall trees if of bad quality (VDFV 1902). On all plots the slopes of the diameter growth of the trees that survived until 2010 increased with their social position; i.e., in tendency taller trees have steeper courses of diameter growth.
Interestingly, all point clouds (Fig. 4a-c) widen with increasing percentile in 1904, i.e., trees with percentiles > 30% can both steeply or hardly grow in diameter. Despite their advantageous social positions, some trees show very shallow courses of growth (black points in the lower right edge), other trees show steep courses although they are only co-dominant or subdominant (black points in the upper layer in the middle). In summary this means that very low social position at the beginning mostly caused drop-out, however, with increasing social status the potential of social climbing as well as losing increases continuously.
Thinning generally increased the slope (level of the point clouds), the number of removal trees (red dots), and increased to some extent the probability of small trees to survive or even strongly climb (black points in the left lower edge in Fig. 4b and c). Figure 5 shows the percentiles 2010 plotted over the percentiles 1904 for the three thinning grades. Observations on the bisectoral line would indicate stationarity of the ranking. The black straight line shows the observed non-stationarity of the size ranking. In case of the A-grade the observed slope ( ±SE ) was 0.6807 ± 0.0613, this results in a 99% confidence interval (p < 0.001) of 0.4969-0.8645 that does not include 1.0 and indicates a significant deviation from stationarity ( Table 3). The respective values for the B-grade were 0.3600 ± 0.1018 with an interval of 0.0546-0.6654. The values for the C-grade were 0.5813 ± 0.1088 with an interval of 0.2549-0.9077. The analogous analyses for the pooled data (A-, B-, C-grade combined) yielded 0.5653 ± 0.0487 with an interval of 0.4192-0.7114 (Table 3). This means in all cases a highly significant (p < 0.001) deviation from stationarity, i.e., a high social drift.
All slopes a 1 were significantly < 1.0, and indicate that subdominant trees in 1904 can often socially climb up, whereas dominant ones in 1904 can lose their privileged status. This tendency was stronger on the B-, and C-grade plots than on the A-grade. For more detailed statistical characteristics of the underlying regression and correlation see Table 3.

Differences between the courses of the stem growth of the 100 tallest trees at the beginning and end of the survey period (H2)
It has been frequently shown that small and suppressed trees, which represent low percentiles in the diameter distribution, are frequently outcompeted under self-thinning conditions or removed by moderate or heavy thinning from below as shown in Fig. 1a-c. More interesting is how the dominant trees, e. g., the 100 tallest trees develop, as they dominate the stand structure and they also represent a significant portion of the standing volume and volume growth.
On the A-, B-, and C-grade plots, respectively, 40, 60, and 44% of the 100 tallest trees in 1904 were replaced by others till 2010; with other words, they were replaced by social climbers and down-ranked below the 100 tallest. Figure 6 shows the growth trajectories of the climbers and the social losers. Figure 6a (red broken lines) visualizes for the A-grade (self-thinning) the stem diameter growth of those trees which belonged to the 100 tallest trees in 1904 but lost this status until 2010; most of them were comparatively tall at the beginning but have a shallow course of diameter growth.
The black lines, in contrast, represent the social climbers, i.e., those trees that did not belong to the 100 tallest in 1904 but in 2010. All of them were rather small at the beginning but steep in the diameter growth from 1904 to 2010.
The tendency of the losers and climbers is similar in the stands with moderate and heavy thinning from below (B-and C-grade) (Fig. 6b and c). The losers show even shallower courses of diameter growth and the social climbers have steeper stem diameter curves compared with the A-grade (Fig. 6a). Figure 7 visualizes the mean course of stem diameter growth of the 100 tallest trees at the beginning (black) and at the end of the survey period (red). The differences of the diameter growth trajectories between both social groups (quickstarters versus latecomers) were tested by Model 2 for each thinning grade separately (Supplement Tables 1, 2 and 3) and by Model 3 for all three plots together (Table 4).
On all three plots Model 2 showed that the intercept of the diameter-age relationship of the survivors was significantly lower (p < 0.001) and the slope steeper (p < 0.001) than the diameter-age relationship of the beginners (Supplement Table 1, 2 and 3). The slopes of the survivors were by  Table 3). The bisectoral line would indicate stationarity of the ranking  Table 4 shows the analogous statistical characteristics of the pooled dataset analyzed by Model 3, i.e., the test of survivors versus beginners over all three thinning grades. Compared with the quickstarters (group = 0) the latecomers (group = 1) had a lower intercept ( a 2 = −7.5925 ) and a higher slope ( a 3 = 0.0565 ). This overall analysis underpins the significantly lower intercept (p < 0.001), steeper slope (p < 0.001); the percental increase in the slope was on average 23.37%. The reported increase in steepness in percent resulted from the division of the delta slope of the survivors ( a 1 = 0.0565 ) by the slope of the beginners ( a 0 = 0.2412). In summary, all analyses revealed a significant effect of the social group on both intercept and slope of the d − a ge relationship. Compared with the quickstarters (group = 0) the latecomers (group = 1) had a lower intercept and a steeper slope. With other words, the survivors have a steeper and longer lasting growth then the beginners.
The comparison between the characteristics of the social climbers and the originally dominant trees in Table 5 revealed that the climbers, although subdominant in the past, achieved similar stem and crown characteristics till the last survey in 2010. They had only minor and non-significant shorter crowns and lower crown ratios (Table 5, column 3   Table 4 Results of fitting the linear mixed effect model of tree diameter dependent on tree age of quickstarters (group 0) versus latecomers (group 1) to the combined A-, B-, and C-grade data (Model 3) The underlying model was d ijk = a 0 + a 1 × age ijk + a 2 × group ij + a 3 × age ijk × group ij + +b i + b ij + ijk . AIC comparisons suggested using random effects at the plot and tree level. Number of observations = 3578 from 216 trees on 3 plots.  vs. 4). The social losers, although dominant in 1904, were significantly lower in tree diameter, height, crown diameter and crown projection area compared with the 100 tallest final trees in 2010 (Table 5, column 5 vs. 6). The stem slenderness of the losers was significantly higher (slenderness of 0.78 versus 0.66) due to their decreased social position within the stands.

Effect of the social drift on the slope of the current diameter development (H3)
The effects of the past and present social position within the stand on the current slope of the diameter development according to Model 4 are presented in Table 6. The 30-year slope 1990−2010 of the diameter development was chosen, as growth trend analyses often cover a 20-year period. Table 6 shows that the growth trend in terms of the slope of the diameter development from 1990 to 2010 (mean stem diameter increase in cm yr −1 ) was positively influenced by the current diameter percentile ( perc 2010 ). The higher the current social position of the trees ( perc 2010 ), the steeper their diameter slopes ( a 2 = 0.0026 − 0.0037 ). However, the social tree position in the past ( perc 1904 ) decreased the diameter growth trend at present ( a 1 = −0.0022 to − 0.0012 ). This means that trees that were dominant in the past had a shallower growth trend at present, whereas subdominant trees in the past had steeper diameter growth slopes at present. In all cases the effects were highly significant (p < 0.001) and the R 2 -values ranged from 0.58 to 0.79. Social dominance in 2010 certainly increased the slope of the diameter growth in the period 1990-2010, i. e. taller trees grow better than smaller, subdominant ones. The positive regression coefficients a 2 in Table 6 show that this applies for the trees of all three thinning grades and also for the pooled data of the A-, B-, and C-grade plots together. Figure 6 visualizes this finding for the pooled data (for statistical characteristics see last row in Table 6). However, ceteris paribus, those trees that were subdominant in the past (when first measured individually in 1904) have significantly steeper slopes at present than the former dominant ones. The steepness of the slope of the diameter growth in 1990-2010 increases with decreasing social position in 1904. This means that the present slope is the higher for trees with low ranking in the past and decreases with increasing perc 1904 as follows 5 > 25 > 50 > 75 > 95 (see Fig. 8). Obviously, slow growth in early age is associated with longer continuing growth in advanced tree age. Table 7 summarizes the relative effect of the past social position (percentile in 1904) and the present social position (percentile in 2010) of trees on their current growth trend (slope in 1990-2010). It reflects for instance that currently dominant trees (perc 2010 = 100%) that were co-dominant in the past (perc 1904 = 90%) have a by 8% steeper slope at present. If the same trees were only subdominant or even  Table 6 Table 7 8  9  11  13  16  21  80  16  18  21  26  32  42  70  24  28  32  39  48  64  60  32  37  43  51  64  85  50  40  46  54  64  80  106 suppressed in 1904 (perc 1904 = 80% or 50%), the slope is by 40% steeper than the one of trees that were predominant in 1904 (second column in Table 7). The analogous, but even stronger positive effect of the past social position on the present growth trend is also shown for currently codominant, subdominant, and suppressed trees (perc 2010 = 90-50%, from the third to seventh column).

Growth trend bias by retrospective sampling of the survivors in mature forest stands
This study was mainly motivated by a series of recently published papers which detected growth events and trends methodologically based on increment coring of survivors in mature forest stands (Bosela et al. 2016;del Río et al. 2014del Río et al. , 2017. The revealed growth trends were mainly ascribed to the extension of the growing season, warming, increase in nitrogen deposition and the atmospheric CO 2 -concentration (Dulamsuren et al. 2017;Pretzsch et al. 2020a-c). For instance, Pretzsch et al. (2020a and b) reported for European beech in Central Europe a diameter growth increase of 30-40% in the last 100 years. The lastmentioned studies were based on retrospective tree ring analyses by increment boring of dominant survivor trees in mature and older beech stands and raised the question to what extend the reported growth trends may be a sampling bias. The study in hand found that the growth trend of the survivors in a 188 years old stand was on average by 23% steeper than the growth of the original in 1904 dominant trees. Compared with the on average 23% overestimation resulting from the survivor bias, the observed increase of 30-40% is still considerably higher. However, about a half of the diagnosed growth acceleration ascribed to environmental changes may be simply caused by social drift and results in a survivor bias when stand growth is analyzed. The results are relevant for understanding and improving the sampling design. In stands with target-diameter cutting, which has become a very common method of prolongated final harvest, the tall trees with already shallower diameter growth trajectories are continuously removed. Sampling of the remaining trees in such stands would likely estimate a growth curve and trend that may be representative for the survivors. Table 7 corroborates that the probability of an over-estimation of the current growth trend in the analyzed European beech stand increases when mainly trees are sampled that are currently subdominant. And the over-estimation is the higher the more subdominant the sample trees were in the past. So, the less dominant trees were in the past, the higher is their growth trend at present, and the less representative is their trend for the long-term dominant trees of a stand.
The main cause of the diverging growth trend of survivors from originally dominant trees is that tree growth depends rather on tree size or mass than on tree age (Pretzsch 2020c;Stephenson et al. 2014). The faster trees grow at the beginning the earlier they may slow down later; this has been shown repeatedly for forest (Guttenberg 1915;Pretzsch 2020a-c;Zeide 1993) and urban trees (Smith et al. 2019). Recent studies stress that the growth of many tree species is accelerated and continues longer than expected and reported in textbooks (Pretzsch 2020c ;Sillett et al. 2015a and b;Stephenson et al. 2014). Probably the revealed social drift is reinforced by this growth acceleration, i.e., the survivor bias may even increase with increasing growth trend as the size growth of the dominant trees is accelerated, their slowing down and removal preponed and the social climbing of their latecoming neighbors, preferentially sampled in older stands, may be promoted.

Implication for tree and stand dynamics
The study in hand shows that the effects of social positions and other legacy effects in the past on the present growth are relevant and that their revelation requires long-term survey on experimental spatially explicitly inventored permanent plots. For analyzing and modelling individual tree growth this information is hardly substitutable by temporary plots or retrospective analyses based on tree ring analyses. Best insight into the effects of the trees' history on their behavior at present and for respective model parameterization provides long-term experiments with repeated measurements of all trees, their positions, survival status and sizes (Nagel et al. 2012;Pretzsch et al. 2019).
Our study corroborates that individual tree growth may be followed and described along tree age, but that it is primarily depending on current tree size. Suppressed trees with lagged size development can continue long in growth, quickstarters may slowdown earlier. This confirms the concept of allometrically based growth models that use trees size instead of age as main predictor (Niklas 1994;Pretzsch et al. 2002). The results also suggest that tree history in terms of past social position, social drift, inner tree structure and tree morphology can significantly co-determine current tree growth. That tree history may be relevant for understanding growth was recently shown by Camarero et al. (2018) and Nothdurft and Engel (2020) who introduced lag models for consideration of such effects.
Models that consider the individual trees' growth courses and their determination by present and past characteristics may better predict the asynchronous rhythm of growth trajectories as show for the individual tree level in Fig. 6. The social drift and asynchrony in growth caused by a combination of quickstarters and latecomers may result in a prolonged and stabilized growth as shown for the stand level in Fig. 1. Such a stabilization and persistence of stand growth caused by the interplay of quickstarters and latecomers is hardly predictable by stand level models based on mean tree or sum values (Hanewinkel et al. 2000;Lindner et al. 1997;Pretzsch et al. 2008).

Silvicultural consequences of the social drift
Of considerable relevance for the silviculture of European beech stands were the long persisting courses of individual growth, the high social drift of the trees during the 106-year survey period, and the finding that thinning can considerably break the social order and ranking and can promote the social climbing. Interesting were furthermore the significantly positive effects of the early social status on the growth trend in the advanced age and that the social climbers finally reached similar tree dimensions as the quickstarters. Beech is the dominant tree species in large parts of temperate Europe and a typical example of a late-successional canopy tree, including remarkable shade tolerance at early life stages (Nagel et al. 2014;Leuschner and Ellenberg 2017). European beech has an especially high ability to survive in the shade and high morphological plasticity (Dieler et al. 2013;Barbeito et al. 2017;Bayer et al. 2013). The potential of social drift that we revealed in this study is supported by the shade tolerance and morphological plasticity of European beech.
It is well documented that European beech can reach longevities of 300-500 years Leuschner 2010, p. 104, Di Filippo et al. 2007); however, tree and stand growth are commonly assumed to decrease at advanced ages (Assmann 1961, p. 80, Mitscherlich 1970. Although the analyzed trees in our study were in advanced ages, their growth did not yet show any slowdown in stem diameter growth. This was in line with the results of Pretzsch et al. (2020a and c), who found a linear diameter growth trend for Norway spruce, silver fir and beech trees in mixed mountain forests. This long continuation of individual tree growth in combination with the ability of social upclimbing may considerably increase the silvicultural flexibility. Quickstarters reduce their growth earlier than latecomers. Thus at advanced age any stand growth reductions by quickstarters may be compensated by latecomers that take their social rank or even overtake them in size. Thus, the combination of quickstarters and latecomers, caused by the social drift, may result in a stable and long continuation of stand growth. The ability to alter the social position even in rather late stand development phases results in a reserve of target trees and flexibility in the timing of harvest. The asynchrony of growth rhythm means that the individual trees arrive at their target size for economically optimal harvest asynchronously. In contrast to plantations where all trees are kept dominant and simultaneously in their social rank and course of growth by wide spacing and permanent density reduction, the analyzed naturally regenerated European beech stands represent more asynchronous individual tree growth patterns. The more variable and asynchronous the growth rhythm, the higher the flexibility in thinning from above (Vuokila 1977) or target diameter cutting (Sterba and Zingg 2001;Zell et al. 2004). Non-stationarity of the trees' social status may on the other hand cause structural heterogeneity, that is disadvantageous for wood quality (Rais et al. 2020). Reduction of social rank can mean crown recession and mortality of branches causing wood quality reduction (Pretzsch and Rais 2016). On the other hand, structural heterogeneity is often associated with richness of habitats and species diversity (Dieler et al. 2017).

Conclusion
The social drift of trees during stand development resulted partly from the dependency of the tree growth from tree size. Quickstarters decreased earlier and latecomers started slower but proceeded longer. This may cause a drift of social ranking over time in terms of changing size dominance. It can stabilize stand growth and increase the silvicultural flexibility of stand management regarding final harvest, threshold diameter cutting, natural regeneration, and transformation.
However, the non-stationarity of social status and asynchrony of growth rhythm may cause a considerable sampling bias if the growth trend of dominant trees or of the stand as a whole is deduced from the latecomers and survivors. As their course of growth is still steeper than the main collective of the trees in a stand, sampling of the latecomers and survivors may cause an overestimation of the current growth trend. In this study the sampling bias of growth trend in terms of the slope of the 100 dominant survivors in relation to the slope of the 100 originally dominant trees was + 23% on average. Analyses of growth trends based on long-term experiments may avoid such bias. However, long-term experiments with full tree inventory are rare and science often has to rely on retrospective analyses via increment boring.
In such cases the revealed bias should be kept in mind for cautious interpretation or adjustment of the quantified trends. Sampling that aims at the analyses of growth trends of the dominant trees that are relatively unbiased by competitive effects should therefore select those dominant trees that always belonged to the dominant collective. They are easy to choose in case of permanent plots, in case of sampling by increment coring it may help to select those trees for further analyses that belong to the dominant collective at present and whose reconstructed diameters indicate also a membership to the dominating collective in the past. Trees with ring width pattern that indicate a slow start or even a suppression phase are likely to be latecomers and social climbers that flaw any growth trend analyses. They should be excluded if the study aims at revelation of reliable trends of the long-term dominant trees.
Ongoing studies based on long-term experiments of other species, thinning methods, and site conditions may further quantify the potential survivor sampling bias and other consequences of the social drift of individual trees during stand development.