Global woodland structure from local interactions: new nearest-neighbour functions for understanding the ontogenesis of global forest structure

A number of hypotheses and theories, such as the Janzen-Connell hypothesis, have been proposed to explain the natural maintenance of biodiversity in tropical and temperate forest ecosystems. However, to date the details of the processes behind this natural maintenance are still unclear. Recently two new nearest-neighbour characteristics were proposed and in this paper we demonstrate how they can contribute to a better understanding of the ontogenesis of global forest structure from localised neighbourhoods. We applied the new species and size segregation functions together with appropriate test procedures to four example woodland data sets from China at Daqingshan, Jiaohe, Jiulongshan and Xiaolongshan forest regions. In addition we quantified the morphology of the new characteristics and modelled a neighbourhood allometric coefficient linking the two functions. The results revealed quite different species segregation patterns with both conspecific and heterospecific attraction. We found these to be generally matched by equivalent size segregation patterns of attraction of similar and different sizes. It was straightforward to model the size segregation function from the knowledge of the species segregation function by estimating a neighbourhood allometric coefficient. The new characteristics have helped to quantify the extent and rate of decline of neighbourhood interactions in terms of spatial species and size diversity. Through the allometric neighbourhood coefficient the analysis highlighted once more how closely related species and size segregation are, thus supporting the mingling-size hypothesis. Using both a traditional and a restricted random-labelling test has provided a valuable tool for understanding the exact nature of species-mingling and size-inequality relationships.


Background
Many theories and hypotheses have been proposed to explain the maintenance of plant and tree diversity in forest ecosystems (Pommerening and Uria-Diez 2017).
The Janzen-Connell hypothesis (Janzen 1970;Connell 1971), for example, proposes that elevated numbers of specialist natural enemies, such as herbivores and pathogens, maintain diversity in plant communities. They reduce the survival rates of conspecific seeds and seedlings located close to reproductive adults or in areas of high conspecific density (Comita et al. 2014) leading to elevated conspecific self-thinning (Yao et al. 2016), i.e. a progressive decline in density in a population of growing individuals of the same species (Begon et al. 2006, p. 157). An important effect of the Janzen-Connell hypothesis is the negative density/distance dependence that occurs when nearby conspecific plants negatively affect performance through mechanisms such as intraspecific competition and pest facilitation (Wills et al. 1997;Wright 2002;Piao et al. 2013;Yao et al. 2016).
Another important ecological hypothesis in this context, the herd immunity hypothesis, focuses on the variation in heterospecific neighbour densities and predicts that species diversity confers protection from natural enemies by making it more difficult for specialist natural enemies to locate host plants (Wills et al. 1997). Thus the spatial spread of an infection can be slowed down or even stopped by mixtures of susceptible and resistant species (Begon et al. 2006). Therefore individual plant fitness should be enhanced in stands involving many species but reduced in stands with low species diversity (Wills et al. 1997;Murphy et al. 2016).
In addition, Ford (1975) and Weiner and Solbrig (1984) demonstrated how the Janzen-Connell and the herd immunity hypotheses can affect natural populations. One observation in many natural plant populations is that self-thinning leads to size inequality (also referred to as local size hierarchies), where heterospecific stands include dominant plants emerging from a first colonisation cohort, which are often surrounded by patches of smaller sized plants of the same cohort. Small plants of these early colonisers are initially often of the same species as the dominant plants and according to the mechanisms of both the Janzen-Connell and herd immunity hypotheses later decrease in numbers due to self-thinning processes. Eventually the small early colonisers are partially or completely replaced by even smaller individuals of other species from subsequent colonisation cohorts. This combined effect of species and size replacement enforces both local size inequality and the mingling of different plant species in a given area or patch and prevents the development of monocultures (Pommerening and Grabarnik 2019).
Synthesising the key elements of the aforementioned ecological hypotheses, Pommerening and Uria-Diez (2017) and Wang et al. (2018) have concluded thatas a consequence of themthere is often a tendency for plants with high species mingling, i.e. those with heterospecific nearest neighbours, to be larger sized plants and referred to this as the mingling-size hypothesis.
Some aspects of these theories could be confirmed in experiments but many processes are still largely unknown. However, the effects described by the Janzen-Connell, herd-immunity, size inequality and minglingsize hypotheses are strictly local effects, i.e. their spatial range hardly exceeds 8-10 m in temperate and boreal forests. It is therefore interesting and important to study how the local woodland structures emerging from these localised effects cooperate in space and time and influence the global structure of a forest stand or woodland. Of particular interest is the question of how localised clusters of tree vegetation evolve in the vicinity of large parent trees, how these localised clusters interact and eventually merge to form overall global stand structure. This question has to be pursued using spatial approaches, since different spatial scales are involved. As typical local clusters of 1-2 mature trees surrounded by smaller offspring essentially constitute neighbourhoods, nearest-neighbour approaches seem most appropriate. These methods are part of point process statistics, one of the main fields of spatial statistics, and good overviews are offered in the textbooks by Illian et al. (2008), Wiegand and Moloney (2014) and Pommerening and Grabarnik (2019).
The objective of this paper is to show how recently proposed nearest-neighbour characteristics can contribute to tracing the ontogenesis of global forest stand structure from localised neighbourhoods and particularly how measures of species and size diversity are related.

Methods
Nearest-neighbour approaches of species and size diversity Gadow (1993) and Aguirre et al. (2003) extended the species segregation index by Pielou (1977) to general multivariate species patterns involving k neighbours. The mingling index M i (Eq. 1) is defined as the mean heterospecific fraction of plants among the k nearest neighbours of a given plant i. In the analysis, every plant within a given research plot acts once as plant i (Pommerening et al. 2019).
The mingling approach is particularly useful in woodlands with many species, because the only distinction made is between co-and heterospecific pairs of trees.
Expected mingling, EM, (implying independent species marks) as explained in Pommerening et al. (2019) and in analogy to Pielou's segregation index, can now be used to normalise species mingling in such a way that the resulting measure Ψ is independent of the number of species and refers to the null hypothesis of independent species marks according to Pommerening and Uria-Diez (2017) referred to the index defined in Eq. (2) as species segregation index. Here M is mean mingling of a tree population. Consequently, Ψ = 0, if the species are independently or randomly dispersed. If the nearest neighbours and a given tree always share the same species, Ψ = 1 (attraction of the same species leading to species segregation). Ψ = − 1 is approached, when many neighbours are increasingly of a different species than the tree under consideration (attraction of different species).
Since the number k of nearest neighbours is related to different scales of local neighbourhood that can vary considerably within and between forest stands, Pommerening et al. (2019) considered a function Ψ(k) that depends on variable k (Eq. 3).
Ψ ðkÞ ¼ 1− MðkÞ EM ; k ¼ 1; 2; 3; ::: With increasing k, Ψ(k) → 0. Our hypothesis is that the exact way in which Ψ(k) approaches 0 is a better quantitative description of the mingling pattern in a given population than M or Ψ alone. For plant patterns with completely random species dispersal (independent species marks) Ψ(k) = 0 for all k (Pommerening et al. 2019).
To better account for different spatial scales of local neighbourhood, we replaced k by the population mean distance r k ¼ 1 N P N i r ik between any individual i and its k th nearest neighbour. This resulted in the new scale-dependent species segregation function Ψ'(r) (Eq. 4).
The interpretation of the function is straightforward and basically follows the interpretation of Ψ (Eq. 2).
An appropriate test involves the null hypothesis of a priori marking or random superposition, also referred to as population independence (Illian et al. 2008;Pommerening and Grabarnik 2019). This requires simulating n = 2499 spatial patterns with independent marks for estimating global envelopes (Myllymäki and Mrkvička 2019). Since species is the mark of interest, spatial mark independence is simulated by random shifts of species populations (Illian et al. 2008, p. 460 f.;Pommerening and Grabarnik 2019, p. 182 f.). For this purpose we selected as many species as required to shift approximately half of all points in the observation window. In the simulations, all individuals of the selected species are shifted by adding the same random values z x and z y to the x and y coordinates of these individuals. A variant of periodic boundary conditions ensures that all points are inside the observation window.
Size differentiation is defined as the mean of the ratio of smaller and larger plant sizes u of the k nearest neighbours subtracted from one. In the analysis, every plant within a given research plot acts once as plant i (Gadow 1993).
Here u can be any quantifiable plant size measure. The value of T i increases with increasing average size difference between neighbouring trees. T i = 0 implies that neighbouring trees have equal size, whilst larger size differences approach an upper value of 1.
Expected size differentiation across species, ET, is independent of the number of neighbours, k, and can be combined with the arithmetic mean T in an index ϒ expressing the relationship between observed size differentiation and spatially uncorrelated size marks according to Pommerening and Uria-Diez (2017) referred to the index defined in Eq. (6) as size segregation index. Since the index takes the expectation into account, it is normalised for the range of sizes involved in each study population, which the original size differentiation index is not. Consequently, Υ = 0, if plant sizes are independently dispersed without any spatial correlation of marks. If the sizes of the nearest neighbours and a given plant are more frequently of similar size than expected by chance, Υ ≈ 1 (attraction or aggregation of similar sizes leading to a segregation of sizes). Negative values tending towards − 1 are achieved, if more neighbours have sizes quite different from that of a tree under study (attraction of different sizes; aggregation in the classical sense of Pielou's original approach) than expected in a pattern of independently dispersed sizes (Wang et al. 2020).
In analogy to the species segregation function Ψ 0 ðrÞ , we can consider a function Υ(k) and eventually the size segregation function Υ 0 ðrÞ (Wang et al. 2020, Eq. 7).
With increasing r, Υ 0 ðrÞ→ 0, i.e. spatial independence of size marks is approached with distance used in the index calculation. The interpretation of Υ 0 ðrÞ follows that of Υ. For plant patterns with spatially independent plant sizes, Υ 0 ðrÞ ¼ 0 for all r (Wang et al. 2020).
As with any quantitative marks, the null hypothesis relates to a posteriori marking or random labelling. Simulations under this mark independence hypothesis are typically based on fixed point locations and permutated marks. According to the traditional random-labelling method, all size marks are freely permutated without restriction. However, when multivariate patterns involving several species are studied, it is common to restrict random labelling in such a way that plant sizes are only permutated within each species population (Wiegand and Moloney 2014, p. 227 f.; Wang et al. 2020). As a result the non-spatial empirical size distribution of each species are preserved. Also here we applied n = 2499 simulations for estimating global envelopes (Myllymäki and Mrkvička 2019).
Our hypothesis is that the exact way in which Ψ 0 ðrÞ and Υ 0 ðrÞ approach 0, i.e. the morphology of Ψ 0 ðrÞ and Υ 0 ðrÞ, and the two aforementioned random-labelling tests give vital clues about local pattern formation in forests.

Characterising the morphology of species and size segregation functions
For a better quantitative description and interpretation of the function morphology we performed a linear regression on the values of Ψ 0 ðrÞ and Υ 0 ðrÞ at r k for k =1, 2, 3, …, 30. We considered the intercept and the slope of the linear regression as shape parameters, where the intercept is an expression of maximum species or size segregation (Pommerening et al. 2019). The slope parameter expresses the average slope of the (in most cases nonlinear) function, i.e. how rapidly on average the functions tend towards a value of 0 denoting spatial species or size independence. This characteristic is therefore a measure of decay of the species or size effect with distance. In addition we calculated the range r k −r 1 , i.e. the difference between the average distance to the first and the k th nearest neighbour, which we set to k = 30. Short ranges occur in tree patterns with high densities and in patterns that are highly clustered whilst long ranges are common in those with low densities. Finally we also fitted polynomials of degree two to the values of the two functions at r k and calculated the sum of squares of deviation (SSDL, Liu et al. 2009) from linearity (Eq. 8).
The values of f P j ðrÞ and f L j ðrÞ (Eq. 8) were calculated for the j th nearest neighbour from the polynomial (P) and linear (L) models, respectively, where f(r) can refer to both species or size segregation function. Deviation from linearity as expressed by SSDL is a measure of the decreasing effect of species segregation or size inequality with distance. Nonlinear curves of f(r) often indicate a fast decline of neighbourhood effect with distance (Pommerening et al. 2019;Wang et al. 2020).

Allometric relationship between species segregation and size inequality
Spatial species segregation and spatial size inequality act on the same neighbourhood. As work on confirming the size-mingling hypothesis by Pommerening and Uria-Diez (2017) and Wang et al. (2018) has shown, size inequality and species mingling are often closely related. This was an interesting finding highlighting that in mixed-species forests large size inequality can be found where there is large species mingling and vice versa. This association of the two aspects of tree diversity suggests that they can be interpreted as allometric neighbourhood relationships, i.e. it should be possible to estimate one of the two functions from the knowledge of the other one using an allometric coefficient. To estimate the allometric coefficient, first the equivalents of relative growth rates (RGR) are calculated from Ψ 0 ðrÞ and Υ 0 ðrÞ at successive r k thus exchanging time for space and quantifying rates of decline of the respective tree diversity effect with distance: RGRs for Υ 0 ðrÞ, p ðΥ Þ r k , can be calculated in analogy. Finally the allometric coefficient is calculated as yielding the same results as Eq. (10). For all calculations we have written our own code in R (version 3.5.1), C++ and additionally used the GET package (R Development Core Team 2019; Myllymäki and Mrkvička 2019).

Example data
To illustrate the new nearest-neighbour characteristics, the associated statistical tests and other curve characteristics, we selected four example forest sites from China.
Daqingshan forest region (abbreviated as D) forms a part of the Daqingshan Forest Farm of the Experimental Center of Tropical Forestry, Chinese Academy of Forestry. The research area is situated in Pingxiang City, Guangxi Province, which is close to the border between China and Vietnam. The native vegetation type is characterised by southern subtropical monsoon forests and evergreen broad-leaved forests involving diverse broad- Jiaohe forest region (abbreviated as J) is an experimental forest situated in the Dongdapo Nature Reserve (43°51′-44°05′ N and 127°35′-127°51′ E), Jilin Province, north-eastern China. The climate is largely influenced by monsoons from the Pacific Ocean and cold air from the interior Asian continent. The stand, denoted as Jb, is a mixed Fraxinus mandshurica-Juglans mandshurica forest. It mainly includes Fraxinus mandshurica RUPR., Juglans mandshurica MAXIM., Acer mandshurica MAXIM., Carpinus cordata BL., Acer mono MAXIM. and Pinus koraiensis SIEB. ET ZUCC.
Jiulongshan Forest (abbreviated as JS) is located in the western suburbs of Beijing (39°57′ N and 116°05′ E) in the northern branch of Taihang Mountain. The climate in this region is temperate continental and largely influenced by monsoon conditions. The selected stand, JSa, is dominated by planted Platycladus orientalis (L.) FRANCO and is mixed with some naturally regenerated species such as Quercus variabilis BLUME, Broussonetia papyrifera (L.) VENT., Ailanthus altissima (MILL.) Swingle, Prunus davidiana CARR. and Gleditsia sinensis LAM.
Xiaolongshan Forest (abbreviated as XS) is part of the Xiaolongshan Nature Reserve, Gansu Province, northwest China. The forest is situated on the north-facing slopes of the West Qinling Mountain range (33°30′-34°49′ N and 104°22′-106°43′ E) and constitutes a natural mixed pine-oak forest. The research plot labelled as XSc represents a natural deciduous broad-leaved mixed forest and Quercus aliena var. acuteserrata MAXIM., Dendrobenthamia japonica (DC.) FANG. var. chinensis and Acer davidii FRANCH. are the most abundant species. This forest stand was restored since the 1970s after commercial harvesting (Pommerening et al. 2019;Wang et al. 2020).

Results
The four example research plots have tree densities that range from 748 tree per hectare (at Jiaohe, plot b) to 2331 tree per hectare (at Jiulongshan, plot a) and from 20.3 m 2 ·ha − 1 (at Jiulongshan, plot a) to 33.2 m 2 ·ha − 1 (at Daqingshan, plot a) ( Table 1).
The number of species varies from 9 (at Jiulongshan, plot a) to 57 (at Daqingshan, plot a). The stem-diameter coefficient of variation is highest at Jiaohe (0.7, plot b) and lowest at Jiulongshan (0.3, plot a). For the remainder of this text we will mostly drop the plot affiliation for ease of reading. The hat symbol, e.g. inΥ , used throughout the following test denotes estimators.

Mean mingling, b
M , is highest at Daqingshan (0.78) and Jiaohe (0.79) and lowest at Jiulongshan (0.20). In addition the species segregation index,Ψ , highlights that there is weak heterospecific attraction at Daqingshan and mainly conspecific attraction in the other three research plots, particularly at Xiaolongshan. The hints given by the species segregation index are much more clearly elaborated on by the species segregation function Ψ 0 ðrÞ (Fig. 1): There is indeed a mild attraction of heterospecific individuals at Daqingshan Forest up to a distance of r = 5 m, while in all other forest areas there is conspecific attraction only. The different ranges ofΨ 0 ðrÞ reflect the different tree densities at the respective sites in terms of number of trees per hectare. Figure 1 also clearly shows different shapes ofΨ 0 ðrÞ , even where the general attraction pattern is the same. Particularly the decline ofΨ 0 ðrÞ with distance is of im-  portance. At Jiaohe, there is a rather gradual decline of the mingling effect with distance, while the decline is rapid at Xiaolongshan. The random superposition simulations highlight that observedΨ 0 ðrÞ is significant at Daqingshan, Jiaohe, Jiulongshan and Xiaolongshan at least for some distances r. The curves in Fig. 1 also emphasise that even with k = 30 nearest neighbours the correlation range (beyond which species marks are independent) has not yet been reached for all forest patterns with the exception of Daqingshan, since the corresponding curves have not yet closely and asymptotically approached the horizontal line through 0, which would indicate spatial independence.
Like with mean mingling, b M; mean diameter differentiation, b T , is also largest at Jiaohe (0.45) and Daqingshan (0.43) and smallest at Jiulongshan (0.25), thus supporting the mingling-size hypothesis. At Daqingshan Forest, there is a moderate attraction of different tree sizes. HereΨ andΥ are even of the same magnitude. At Jiaohe and Jiulongshan the size segregation index suggests an attraction of similar sizes, which is strongest at Jiulongshan (0.05). In plot c at Xiaolongshan, the observed size marks appear to be largely independent. We can refine the impressions obtained fromΥ by studying the size segregation functionΥ 0 ðrÞ (Fig. 2).
Comparing the curves of Fig. 1 and Fig. 2 we can confirm that size segregation generally operates at a scale smaller than species segregation. In Daqingshan Forest, the attraction of different sizes appears to be much stronger than indexΥ suggested, which supports the use of a function providing more details as opposed to a simple index. It is actually the strongest effect of all four example sites. At Jiaohe and Xiaolongshan there is an attraction of similar sizes with different patterns of decay with increasing distance r: At Jiaohe the differentiation effect decreases much more gradual than at Xiaolongshan. By contrast, at Jiulongshan Forest there is even hardly any change with distance.
Not only for testing the significance of the estimated characteristics but also for understanding the spatial relationship between species and size it is instructive to carry out the aforementioned two different types of random labelling tests: The first test is the "traditional" one (results in the left column of Fig. 2), where the different species occurring in the research plots are ignored and all sizes are permuted regardless of species. An important consequence of this test is that the species-specific mark distributions change from simulation to simulation and can largely deviate from those observed in the original plot. The second variant of the test (results in the right column of Fig. 2) is preferred by many quantitative ecologists and statisticians (Wiegand and Moloney 2014, p. 227 f.). Here the permutations are restricted by species boundaries, i.e. tree sizes are permuted within species populations only. As a consequence species-specific mark distributions are preserved in the simulations, i.e. they are exactly the same in each simulation and correspond with those originally observed. However, for the less abundant species this can lead to situations where within the tiny species populations hardly any permutation of sizes is possible.
When comparing the envelopes in the left with those in right column of Fig. 2 we clearly see that they are centred towards zero for the traditional test but not for the restricted random labelling test. In some cases the envelopes of the restricted test follow the observed curves more closely than those of the traditional test, e.g. for Daqingshan and Xiaolongshan Forests. In other cases, such as for Jiaohe and Jiulongshan, the differences in the envelopes of the two tests appear to be marginal. As a consequence, the strongly significant test outcome indicated by the traditional method when applied to Daqingshan turns out to be merely an artefact and the observed size segregation function is only marginally significant when the correct test procedure, i.e. the restricted random-labelling method, is applied. A similar situation exists at Xiaolongshan. However, the question why the test envelopes differ so markedly between the two test methods for some patterns but not for others remains and will be discussed in the next section, as it is rather important.
The modelling ofΥ 0 ðrÞ using the mean allometric coefficientm is straightforward and leads to good results (Fig. 3). This emphasises that species and size segregation functions are indeed related. RMSE is a fit characteristic and generally shows that the modelling based on only one model parameter, i.e. the mean allometric coefficient, was satisfactory. The curves show the general trend that the values ofΥ 0 ðrÞ are usually much smaller than those ofΨ 0 ðrÞ with the notable exception of Daqingshan Forest where both are of similar magnitude. The mean allometric coefficient is in all four cases close to 1 indicating an isometric decline (increase in Daqingshan Forest) of the size differentiation and mingling effects in tree neighbourhoods. The isometric relationship is particularly strong in Xiaolongshan (m ¼ 1:00 ) and Daqingshan (m ¼ 1:01). Deviations from isometry of 3% occur in Jiaohe and Jiulongshan.
The morphology of the species and size segregation functions is of particular interest for making inference about processes of natural maintenance of tree diversity. Table 2 summarises four characteristics of morphology, intercept and slope from applying a linear regression model, the spatial range of the two segregation functions and SSDL, a measure of deviation from linearity.
The differences in intercept between the two functions reflect the aforementioned differences in the function values. The signs in the slopes are consistent betweenΨ 0 ðrÞ andΥ 0 ðrÞ, only for Jiulongshan the slope of the linear fit is almost zero thus confirming the visual impression we gained from Fig. 2. Interestingly the absolute values of slope for both functions are only close for Daqingshan but not for the other three woodlands. Perhaps the slope of a linear approximation is not the most effective way to estimate average decline of interactions. The range is naturally the same for both Fig. 2 The size segregation functionΥ 0 ðrÞ based on stem diameters as size marks for the four Chinese example research plots Daqingshan (plot a), Jiaohe (plot b), Jiulongshan (plot a) and Xiaolongshan (plot c) including the 95% global envelopes from 2499 traditional random labelling simulations (grey, left) and from 2499 random labelling simulations restricted by species (grey, right) functions and is a reflection of tree density as shown in Table 1. The greatest deviation from linearity for bothΨ 0 ðrÞ andΥ 0 ðrÞ can be confirmed for Xiaolongshan and Daqingshan. This implies that the decline of mingling and size inequality effect with distance is greatest here whilst it is much more gradual for Jiaohe and Jiulongshan.

Discussion and conclusions
Our analyses have clearly shown that the species and the size segregation functions are a very useful addition to existing point process characteristics, since they relate nearest-neighbour characteristics to spatial scales. There are only few existing characteristics in point process statistics that combine neighbourhood and spatial scales (Wiegand and Moloney 2014). Our new segregation functions are nearest-neighbour characteristics that depend on distance and thus give valuable insights on how mingling and size differentiation effects decrease with distance. This information is crucial for understanding how local forest structure evolves and eventually contributes to global forest structure. As part of this, morphology of the segregation functions correlation range and nonlinearity are of particular interest. Applying the new methods to the four example woodlands has demonstrated the extent of neighbourhood effects and that the decline of these effects largely differs between them. Of particular interest is also the correlation range of the two new functions, since it indicates the spatial extent of influence of local neighbourhood. From applying the two functions we have understood that local neighbourhoods in temperate woodlands in China, such as those that we studied, extend to 8-12 m.
Even more information on local neighbourhoods and their structure can be derived from applying species and size segregation functions simultaneously. From this we have learnt that the size segregation effect is often weaker than the species segregation effect. Also, the morphology of both function graphs is usually very similar to such an extent that we can consider this similarity as an allometric neighbourhood relationship. The mean allometric neighbourhood coefficientm clearly con- firmed the close relationship between the decline rates of both characteristics. This finding strongly supports the mingling-size hypothesis suggesting that trees with high spatial species mingling often tend to be larger sized trees with considerable size differentiation in their vicinity.
The two random labelling test methods applied in Fig.  2 revealed another important aspect of size-species relationships: While the restricted random labelling method, which takes the existence of distinct species populations into account, is undoubtedly the better option for testing spatial mark independence, there were two sites, Jiaohe (plot b) and Jiulongshan (plot a) where the choice of test procedure mattered little. This begs the important question why some forest patterns are sensitive to the random labelling test method and others are not. In situations, where the envelopes greatly differ between tests, the non-spatial size-mark distributions greatly  Table 2 Shape characteristics of the estimated species segregation functionΨ 0 ðrÞ and the size segregation functionΥ 0 ðrÞ. Interceptfrom linear regression, expression of maximum size segregation; Slopefrom linear regression, expression of how rapidlyΨ 0 ðrÞ and Υ 0 ðrÞ tend towards a value of 0; Rangedifference between average distance to first and 30th nearest neighbour, SSDLsum of squares of deviation from linearity (Eq. 8) differ between the most abundant species (Fig. 4). This in turn affects spatial size differentiation where local neighbourhoods include several different species, i.e. where attraction of different species is high. Apparently this is the situation in Daqingshan (plot a). In this plot, unrestricted size permutations across species boundaries produce spatial size differentiation among nearest neighbours that does not match observed size differentiation, hence the envelopes of traditional random labelling are markedly different from the observed curves. Thus the strong effect of quite different sizes at close proximity at Daqingshan is mainly a consequence of species with different size-mark distributions mingling at this range. By contrast, the three other forest sites show a conspecific attraction, i.e. here the within-species size variability is more important than the between-species differences in size range, particularly at Xiaolongshan (plot c). However, at Jiaohe (plot b) and Jiulongshan (plot a) also the size distributions of the most abundant species are quite similar (see Fig. 4) which contributes to an outcome where spatial size differentiation in unrestricted random labelling is quite similar to that produced by restricted random labelling. Figure 4 confirms this explanation, except that in the case of our four example plots the size-mark distributions of the two most abundant species appear to be the most influential ones.
Thus the comparison between the two randomlabelling test methods offers vital insights on how species and size segregation are related in mixed-species forests. It tells a particularly interesting story about the interaction between local and global forest structure that otherwise would not have emerged, if only the seemingly "correct" test were applied.