Patterns of Allocation CSR Plant Functional Types in Northern Europe

Morphological and physiological parameters of 76 vascular plant species typical for Northern Europe were analyzed using Grime’s classification. C species (competitors) have high levels of canopy height, leaf dry weight, and maximal lateral spread. R species (ruderal) have low leaf dry weight, longer flowering period, high rate of photosynthetic capacity and respiration, and high nitrogen content in the leaves. Stress-tolerant (S) species prevailing in habitats with limited resources are small and have low rate of photosynthetic activity and respiration. Principal component analysis (PCA) ordination showed a clear separation of species of different plant functional types according to their morphological and physiological parameters. The first PCA axis showed close relationship with the rate of respiration and photosynthetic activity and allowed us to differentiate S from R species. The second PCA axis correlatedwith morphological parameters associatedwith the size of plants and allowed us to differentiateC species from S and R species. Using PCA ordination, we developed a model that determines plant functional types in Northern Europe and analyzed plant functional types of several species that are not presented in Grime’s classification. The proposed model has higher accuracy (84%) compared to similar models designed for other climatic zones.


Introduction
Prediction of response of individual species, communities, and ecosystems to the global climate changes and anthropogenic impact is one of the main problems of modern ecology [1,2]. In such studies, it is expedient to consider the plant functional groups of assemblages of species with similar responses to changes in environmental conditions. Conception of life strategies developed by Ramenskii and Grime (the theory) is an example of such generalization [3,4]. According to the theory, there are three main types of life strategies: competitors, ruderals, and stress-tolerators.
Species of different plant functional types have different biomass structure, intensity, and duration of growth processes, photosynthetic capacity and respiration rate, use of assimilates, and features of protective, and adaptive mechanisms and reactions [5][6][7][8][9][10].
Grime's system of plant functional types includes many species of vascular plants (more than 1000 species) that are typical for the herbaceous flora of Great Britain [2,11].
The flora of the higher latitudes includes many arctic and alpine species that are not presented in Grime's system. Thus, the classification cannot be used effectively to analyze ecosystems of the northern territories. Moreover, many species have high plasticity and demonstrate different strategic features depending on the environmental conditions (climatic, edaphic, and eco-coenotic), especially on the borders of their distribution areas [12].
Currently, there are several models that identify practical ordination of wild plants within the space. Such models were developed for meadow vegetation of England, Ireland, and France [13]; alpine and subalpine vegetation of Italy [14,15]; and flora of the Great Britain [11,16].
Hodgson et al. 's [11] and Pierce et al. 's [15] models are the most known. These models were used to determine life strategies in different climate zones from alpine meadows in Italy [17] to desert vegetation in Turkey [18]. However, there are no data of these models being used in northern latitudes.
In these models, morphological predictors are used. In our view, morphological parameters are easier to measure 2 International Journal of Ecology than other parameters. However, in some cases, these cannot be used for allocating plant functional types. For example, specific leaf area (SLA) is a key element of the "worldwide leaf economics spectrum" [19]. SLA corresponds to soil fertility and marks stress-tolerant species. However, in Northern Europe where many of the habitats are shaded, the leaf economics spectrum breaks down [20]. We attempt to provide a methodology that gets around this problem using functional (physiological) along with morphological predictors. Functional traits such as photosynthetic capacity, respiration rate, and concentration of carbon and nitrogen in leaves reflect processes of energy transformation and are closely related to plant activity, growth, and productivity [7,9]. Thus, these are closely linked to the plant functional types.
The aim of the present study was to develop mathematical model allocating functional type of vascular plants in Northern Europe based on morphological and functional parameters. . Studies of morphological and physiological parameters were performed during the flowering period (July) in typical communities. Most of the data on physiological parameters (photosynthetic capacity and respiration rate) were obtained previously [10].

Materials and Methods
In addition, we used the data from the TRY DB database (https://www.try-db.org/) [22] in order to eliminate differences in species morphology related to specific habitats. Among all the received data, records from the territory of Central and Northern Europe and Great Britain (more than 17,000 primary records) were selected. The average parameters from all sources (arithmetic mean) were used in the succeeding analysis.

Trait Measurement
2.3.1. Canopy Height. Canopy height was measured from the ground to the top of the generative parts of plants [2]. Measures were conducted with a 30-time replication.

Leaf Dry Weight and Dry Matter Content.
Leaf area and dry weight were measured in mature leaves from the middle of stems. Leaf area was measured using photo images by using ImageJ (https://imagej.nih.gov/ij/). Measures of leaf wet weight were made in the field during 1-2 h after sampling.
Sampled leaves were placed in hermetically sealed plastic bags to reduce the evaporation of water. Leaf dry weight was measured in the laboratory after the leaves were dried for 24 h at a temperature of 80 ∘ C. Dry matter content was calculated as percent ratio of leaf dry weight to leaf wet weight. Measurements of 30-120 samples were replicated according to size and weight of leaves.

Lateral Spread and Phenological
Data. Lateral spread of the species is determined by root system, life form, and distance between the ramets (vegetative individuals) of the plant [11]. We used Raunkiaer's life-form systems to determine the type of lateral spread of each species. Phenology of the model species was studied using http://www.plantarium.ru/ (open atlas of vegetation of Russia and bordering countries).

Functional
Parameters. CO 2 exchange of the leaves was measured using LI-7000 (LI-COR, Inc., Lincoln, NE, USA) and LCPro+ (ADC BioScientific Ltd., Hoddesdon, UK) systems based on the IR analyzer. Small leaves (smaller than the chamber area) were measured in groups. Large leaves (length or width does not fit in the chamber) were cut from the middle without large veins. Gas exchange of excised leaves was measured. Photosynthetic capacity was measured at saturating photosynthetically active radiation (800-1200 mol/m 2 c); respiration rate of the leaves was measured in the dark at a temperature of 18-20 ∘ C. Leaves were sampled from 10 to 15 individual plants. The results for photosynthesis and respiration are presented based on the leaf dry weight.
The content of total nitrogen and carbon was measured in dry ground biomass using elemental CHNS analyzer EA-1110 (Italy) in the ecoanalytical laboratory of the Institute of Biology of Komi SC UB, RAS. Samples were collected in double replication, and each sample was analyzed in double or triple replication.

CSR Classification.
The basis of the classification is three axes that are related to competition, stress tolerance, and ruderality of species (Figure 1(a)). Grime allocated 19 types of primary, secondary, and tertiary plant functional types [2]. It is common to use letter designations for plant functional types (Figure 1(b)). However, in modeling, it is easier to use digital coordinates of species in the three axes (Figure 1(c)).

Model Building Scheme.
In the first step, the status of species (which are presented in the classification) was determined according to Grime's classification [11].

Axis of environmental disturbance
Axis of competition abilities Axis of environmental stress The model was based on multivariate analysis by principal component analysis (PCA) according to their morphological and physiological traits selected in the previous step. The obtained two-dimensional PCA coordinates for all species were then allocated to three axes of space. The axes of space were presented using lines connecting the center of PCA ordination with centroids (mean values of coordinates of species allocated to one group) for , , and groups, respectively.
The projection of PCA coordinates of all species to these three axes was considered as new (adjusted) coordinates in space. In the last step, we identified the strategy for each species as a one of 19 plant functional types that have a minimum distance to the coordinates of species. The distance is calculated by the formula = ( − ) 2 + ( − ) 2 + ( − ) 2 , where , , and are the coordinates of species and , , and are the possible strategies [11].
All analyses were performed for transformed data ( Table 1).
The relationship between parameters was assessed using Spearman's nonparametric rank correlation. Variance analysis, PCA ordination, and correlation analysis were performed using IBM SPSS Statistics 19. Model allocation for plant functional types was realized in Microsoft Excel.

Morphological and Physiological Parameters of Plant
Species of Different Functional Types. Most of the morphological and physiological traits showed significant differences between species of different functional types (Table 2). Mean CH and LDW of species were 2-3 times higher than those of and species. Typical competitors in our studies were perennial grasses (Bromopsis inermis, Calamagrostis  epigejos, and Dactylis glomerata) and nitrophilous plants of rich habitats (Artemisia vulgaris, Urtica dioica, Arctium lappa, and Angelica sylvestris). These species were characterized by a height more than 1 m. SLA did not demonstrate significant differences among plant functional types ( Table 2). The maximal SLA values at 50-60 mm 2 /mg were revealed for species preferring spruce forests (Oxalis acetosella and Trientalis europaea) and ruderal species Stellaria media and Myosotis arvensis (40-50 mm 2 /mg). Minimal SLA values were found for species from the competitor group (15-30 mm 2 /mg). species preferring mire and mountain tundra communities (Rubus chamaemorus, Saussurea alpina, and Hedysarum arcticum) demonstrated low SLA levels.
The functional traits (PN, RD, and LNC) also show a significant difference between functional types. PN and RD of leaves increased in the --species order ( Table 2). The highest level of LNC was revealed in the leaves of species from and plant functional types (3% in dry weight), and this was higher than in leaves of species ( Table 2).
The group of species according to most of the traits (except for FP, FS, and LCC) is in the intermediate position in the system ( Table 2). This fact allows us to use group as the "starting point." Differences in morphological and functional parameters from this "starting point" were interpreted as a "shift" of the species location in Grime's triangle relative to the central point. The next step was to convert the two-dimensional PCA ordination to the relevant , , and axes in Grime's triangle. For this purpose, we projected points to the lines connecting the center of PCA ordination with centroids of , , and groups. The projections of species on these lines were considered as new coordinates in space. For example, the species Dactylis glomerata, with PCA coordinates (0.734 and 1.257), after such projection will receive new coordinates (1.277, −1.360, and −0.022) in the , , and axes, respectively ( Figure 3). These values have a minimum distance to the point (1, −1, and −1), which corresponds to / plant functional type.

Allocating CSR Functional Types in Northern
The developed model is provided in Excel spreadsheet in Supplementary Material Appendix 1.
The coordinates of species in space of our model were closely correlated with their coordinates in Grime's classification. The axis showed a correlation of 0.78 and and axes 0.64 ( < 0.01) ( Table 3).

Analysis of Morphological and Functional Traits and
Their Correlation with CSR Axes. Physical size of plants and Table 3: Correlation between species traits, PCA axes score, and species coordinates in "Grime's triangle."  Table 2. , , and : coordinates of species in the space of Grime's classification. their relative growth rate are considered as key traits in the classification of plant functional types [2]. Competitor species more often establish intact habitats with favorable environmental conditions [2,11]. These plants capture most of the resources and invest them into its own biomass. Leaves of competitor plants are high and/or have large area and weight ( Table 2). Markers of type could be such traits (Table 3) as CH (correlation coefficient with the axis in our model -0.71), LDW (correlation coefficient -0.56), and high LS values of 0.61 In Grime's model, competitors are considered as fastgrowing species with low LDMC values. On the contrary, according to our data, species demonstrate the highest LDMC values. This observation is likely connected with the fact that most of the Northern species cease active growth in July (the month when we take samples) and begin flowering. As growth ceases, cell wall thickening occurs, leading to an increase in the LDMC values.
Stress-tolerators, prevailing in conditions with limited resources, are slow-growing plants with low rate of metabolism [2,5,10,11,23]. Stress-tolerators are small plants with the lowest values of photosynthetic activity and respiration rate that reveal minimal content of nitrogen and maximal content of carbon in the leaves. The highest negative correlations with axis were found for CH, PN, RD, and LNC (Table 3).
It is considered that low SLA values indicate stresstolerant species with relatively low growth rate [24,25]. SLA indicates acquisitive economics of the leaf area development. Low SLA values are typical for plants of type with low metabolism and growth rate [19,26,27]. However, the use of this parameter leads to several issues, called the "Oxalis" effect [16]. The key aspect of this effect is that a group of species includes plants growing under different kinds of stress; therefore, they have different adaptive reactions [28]. In our study, species of tundra and mire communities (Epilobium palustre, Eriophorum vaginatum, Carex dioica, and Potentilla crantzii) growing in an environment with waterlogging and limited nutritional elements were often characterized by thick leathery leaves and demonstrated low SLA values (10-30 mm 2 /mg). However, species preferring spruce forests with limited light (Fragaria vesca, Oxalis acetosella, and Trientalis europaea) had high SLA values (up to 50-60 mm 2 /mg). Therefore, following Gross et al. [29], we considered that SLA should ideally not be used as a marker of plant functional type.
Ruderals prefer rich habitats with weak competition caused by certain disturbance of the vegetation cover. These species demonstrate high growth rate provided by increased metabolic activity: intensive respiration and high nitrogen content in leaves [8,10,23,26,30]. Usually, ruderals are annual and biennial species with seed reproduction, investing most of their resources to reproduction. They occupy disturbed habitats owing to their accelerated ontogenetic development (ephemeroids) and/or extended period of flowering and fruiting [4]. According to our data, low values of lateral spread (LS) and long flowering period (FP) could be markers of plant functional type. Ruderals also demonstrated low LDMC values, high SLA values, and maximal level of PN and RD (Tables 2 and 3). Table 4 shows a comparison of Grime's functional types [2] with plant functional types allocated using Hodgson et al. 's model [11], Pierce et al. 's model [15], and our model. The number of steps (points of tertiary classification, Figure 1(b)) was determined between the species location in our model and that in Grime's system. If functional types matched, there were no steps. Maximal four steps were noted when the coordinates of functional types were in the opposite sides of Grime's triangle.
In general, in Hodgson et al. 's model, the highest deviation (distances 3 and 4) was found for species typical for spruce forests (Oxalis acetosella, Trientalis europaea, Fragaria vesca, and Geum rivale). All these species demonstrated high SLA and low LDMC values and, according to the modeling, were allocated as , / , or / type (Table 4). Obtained difference may be explained by the "Oxalis" effect. Only species with low SLA values (20-30 mm 2 /mg and lower) and high LDMC values (over 30%) may be potentially considered as stress-tolerators based on the analysis of regression equations of Hodgson et al. 's model. In the northern environment, there are species typical for tundra and mire communities, for example, Rubus chamaemorus. It should also be noted that functional types of almost all species allocated using Hodgson et al. 's model were moved to increasing and points and decreasing points (Table 4).
Comparing the model of Pierce et al. [15] with the reference position showed high matching (distances 0 and 1) for 46% of the species (Table 4, Figure 4(b)). Low similarity (distance 2) was revealed for 40% of the species. High deviation (distances 3 and 4) was noted for 14% of the species, including stress-tolerators (Oxalis acetosella and Trientalis europaea) and grass competitors (Bromopsis inermis and Elytrigia repens).
The analysis of the model developed in our research showed the highest matching with the reference species positions given in Grime's classification. Distance between relevant and reference position of species was 0 or 1 for 84% of species. Nine species (16%) demonstrated low matching Table 4: Comparison of reference plant functional type [2] and determined by Hodgson et al. [11], Pierce et al. [15], and our model for * Distance: number of tertiary types between reference strategy and determined by models.
Some species demonstrated high deviation between the reference strategy and our model (distance 2 in Table 4). Epilobium palustre has / type in Grime's classification. However, in our case study, it was allocated in / type. Shifting the species position to the ruderal type was caused by relatively high values of PN (17 mg CO 2 /g dry weight per hour) and RD (3.1 mg CO 2 /g dry weight per hour), which is typical for and species ( Table 2). Change of functional type of Hypericum maculatum from / to / was caused by low values of PN (6.5 mg CO 2 /g dry weight per hour) and RD (1 mg CO 2 /g dry weight per hour), which is typical for plants of functional type.
Some species demonstrating changes in their functional types (Poa annua, Polygonum aviculare, Pimpinella saxifraga, and Vicia sepium) have known not only the position in Grime's classification but also the "area of variability" [2]. For example, P. annua has type in Grime's system. However, its potential "area of variability" (noted by a dotted line in Figure 5(a)) covers the area of --types. In our model, P. annua was allocated in functional type, which is in the area of possible types. Similar results were found for Polygonum aviculare, P. saxifraga, and V. sepium ( Figure 5). Therefore, our model allows to us consider a wide variety of species establishing different habitats with different ecological backgrounds. Averaging of morphological and functional traits for large variety of plants from Russia and Central and Northern Europe allowed us to clarify the impact of specific environment and to analyze the traits determined by genotypic specificity. It should be noted that, by using functional traits closely related to the rate of metabolism, we improved the determination accuracy for stress-tolerators. Comparison of the developed model against other models showed its higher accuracy for the species of Northern Europe. It was revealed that the differences between our model and the reference position were within the limits of potentially possible plant functional types in Grime's triangle. Obtained data on morphological and functional parameters can be used as markers for allocating plant functional types of northern species not presented in Grime's classification. with add-ons. Species with the largest deviation from reference position were selected (distance 2).