Leaf Traits and Aboveground Biomass Variability of Forest Understory Herbaceous Plant Species

Despite increasing use of trait-based approaches in community ecology, most studies do not account for intraspecific variability of functional traits. Although numerous studies investigated functional traits of species with high economic value, the intraspecific and interspecific (caused by species identity) trait variability of forest understory herbs is still poorly understood. We aimed to assess the variability of specific leaf area (SLA), total leaf area, aboveground biomass and leaf mass fraction among 167 forest understory plant species, and the level of variability explained by species identity and collection site. We hypothesized that the level of intraspecific variability of SLA is underestimated in commonly used trait databases and that the interspecific variability (caused by species identity) is greater than intraspecific variability (site-specific). Our study revealed higher interspecific than intraspecific variability of the traits studied. We also confirmed that level of intraspecific variability available in the LEDA database is underestimated. We confirmed that species identity was the main factor determining the values of all the traits studied, and site-specific random effects explained lower amounts of variation in traits. Use of trait values from databases not acknowledging intraspecific variability is biased by uncertainty about this variability. For that reason, our analysis used mean trait values to reduce uncertainty of the results in the study conducted to assess human impacts on ecosystems. Thus, our study might support the assumption that level of intraspecific variability of functional traits is lower than interspecific variability.

species with high economic value, the intraspecific and interspecific (caused by species identity) trait variability of forest understory herbs is still poorly understood. We aimed to assess the variability of specific leaf area (SLA), total leaf area, aboveground biomass and leaf mass fraction among 167 forest understory plant species, and the level of variability explained by species identity and collection site. We hypothesized that the level of intraspecific variability of SLA is underestimated in commonly used trait databases and that the interspecific variability (caused by species identity) is greater than intraspecific variability (site-specific). Our study revealed higher interspecific than intraspecific variability of the traits studied. We also confirmed that level of intraspecific variability available in the LEDA database is underestimated. We confirmed that species identity was the main factor determining the values of all the traits studied, and site-specific random effects explained lower amounts of variation in traits. Use of trait values from databases not acknowledging intraspecific variability is biased by uncertainty about this variability. For that reason, our analysis used mean trait values to reduce uncertainty of the results in the study conducted to assess human impacts on ecosystems. Thus, our study might support the assumption that level of intraspecific variability of functional traits is lower than interspecific variability.
Key words: biomass allocation; functional traits; herb species; interspecific variation; intraspecific variation; leaf economic spectrum; leaf mass fraction; life strategies.

HIGHLIGHTS
We studied leaf traits and aboveground biomass of 167 herbaceous species. We quantified intraspecific and interspecific trait variability.

INTRODUCTION
Revealing the effects of plant diversity and interspecific interactions on related ecosystem functions is a challenging task needed to guide conservation efforts and management of ecological resources (Chapin III and others 2000; Díaz and Cabido 2001). However, the immense number of species and their broad ranges of occurrence make collection of the appropriate amount of data difficult due to the time, labor and money needed for fieldwork (Wilson and others 1999). One of the possible solutions is the use of functional traits, as numerous databases are being developed worldwide. Functional traits are morphological, physiological and phenological features influencing the growth, reproduction and survival of species and populations (Violle and others 2007). The most useful functional traits for running analyses and showing strong correlations with other parameters are mainly specific leaf area (SLA), total leaf area, the amount of nitrogen in leaves, seed mass and plant height (Díaz and others 2016). For example, SLA is crucial for quantifying analyses of life strategies of plants (Westoby 1998;Wilson and others 1999) and for assessing the productivity of ecosystems (Smart and others 2017;Ivanova and others 2018). Total leaf area, as well as other traits, for example, seasonality of foliage or canopy structure, may contribute to an explanation of the variability of aboveground biomass in a grassland ecosystem (Schumacher and Roscher 2009). Aboveground biomass can be used as a proxy for plant responses to harsh conditions in high altitudes (Ma and others 2010). Similarly, leaf mass fraction and other traits connected to biomass allocation may reflect important trade-offs and allometric relationships in plant functioning (Kleyer and others 2019;Poorter and others 2015). The recognition and understanding of correlations between those traits provide the ability to assess plant responses to climate change, or to estimate the possible shifts in species ranges, and their potential threats or invasiveness (Ordonez 2014;Canessa and others 2018;Wuest and others 2018). The achievement of this goal is crucial in the context of predicting the impact of global climatic change on the shape, functioning and diversity of the Earth's ecosystems. Moreover, the advantages of trait-based approaches for ecological studies and biomonitoring in comparison with traditional taxonomic approaches are highly significant, as it meets the general need for better inter-taxon and inter-region comparability (Reich and others 1998a;Cornwell and others 2008;Freschet and others 2011). It is particularly important for developing countries, where funding for research is often extremely limited, but where the adverse effects of human activity on the environment are the largest (Schmidhuber and Tubiello 2007;Mü ller and others 2011).
Most studies prove that SLA shows significant diversity among species and is driven by species identity (Westoby 1998;Poorter and De Jong 1999;Dyderski and others 2016;Kwon and others 2016). On the other hand, the impact of other factors (for example, temperature, precipitation and seasonal dynamics of light availability) on the variability of plant traits has also been highlighted (Herná ndez-Calderó n and others 2014; Popović and others 2016; Blonder and others 2018). For that reason, some researchers decide to measure the trait values directly, taking environmental variability into consideration, although this is only possible when examining small-or medium-sized areas (Ma and others 2010;Carlos and Rossatto 2017;McCoy-Sulentic and others 2017;Taylor and others 2018). Yet, in the face of imperfections of currently used methods, it is crucial to examine the characteristics of plant variability and the amount of variance explained by different factors, so that we can use them in further research, while being conscious of their shortcomings and advantages.
This issue has been investigated across different habitats, always indicating the importance of study scale-the more extensive the study area, the less important the intraspecific variability (Albert and others 2011). In our study, we have adopted forest types from the Polish system of forest management divided into four ecological categories of habitats (fresh deciduous, wet deciduous, fresh coniferous, wet coniferous). On the other hand, in places where diversity is low, intraspecific trait variability gains much more importance, even in vast study areas, for example, the tundra biome (Bjorkman and others 2018). Moreover, recently the importance of intraspecific variability has been highlighted as a significant driver of diversity (Violle and others 2012). However, the general pattern shows that the contribution of interspecific variation to species responses is higher than intraspecific, for example in alpine and subalpine species (Kichenin and others 2013). However, intraspecific variability is less well represented in current databases of plants functional traits than interspecific variability.
Poorer representation of intraspecific than interspecific variability results from the aim of numerous surveys of functional traits, which intended to provide a fast and inexpensive method of data collection about numerous taxa. Extensive datasets obtained using those methods have enabled an easy assessment of defined plant features, for example based on allometric models (Enquist and Niklas 2002;Weiner 2004;others 2012, 2015). Recent studies compiled extensive datasets and aimed to explain relationships among traits shaping economic spectra of plants (Díaz and others 2016). Such relationships allow to conclude about numerous species traits, based on measurement of fewer traits (Reich and others 1998b;Westoby 1998;Peeters 2002;Maré chaux and others 2018). Therefore, the complexity and cost of conducting a study of variability of plants has diminished substantially (Czortek and others 2018; Stevens and others 2018).
The use of datasets in scientific research implies, however, an undefined level of uncertainty resulting from variability of traits among different habitat conditions. For example, in the LEDA database (Kleyer and others 2008) mean values provided for species are often based on measurements of few individuals. In contrast, the Bien database (Enquist and others 2016) is composed of individual plant measurements with (usually) data about collection site coordinates. However, this database covers mainly North and South American species, and for Europe it provides data from LEDA (Kleyer and others 2008).
Whereas surveys of woody species of high economic value have been realized worldwide (Poorter and Remkes 1990;Hunt and Cornelissen 1997;Ketterings and others 2001;Segura and Kanninen 2005;Peichl and Arain 2007; Jagodziń ski and others 2018), our knowledge about species of lower economic importance, mainly herbaceous species, remains low, and the topic has received relatively little attention. Forest understory species, although accounting for comparatively low shares of total forest ecosystem biomass, show up to 20-fold higher shares of micro-and macroelement cycling (Gilliam 2007;Muller 2014). Nevertheless, herbaceous forest species have been examined to a lesser degree than woody vegetation and even other herbaceous plants. Although other functional groups and habitats have been examined in that context, collection site has never been used as a random effect. Furthermore, the contribution of both species identity and collection site to variability of functional traits has never been quanti-fied; only the general pattern has been described (Cordlandwehr and others 2013).
Although plant functional traits have been measured worldwide (Castro-Díez and others 2011; Pierce and others 2013; Díaz and others 2016; Jagodziń ski and others 2016), the proportion of their variability explained by collection site is still poorly understood. Improving our knowledge of functional trait variation is a key goal in the context of inferring results for taxa not studied directly, and for using mean values from databases in broader studies concerning functional diversity of vegetation (Rudner and Gross 2012; Czortek and others 2018; Stevens and others 2018). Thus, the aims of this study were to (1) assess the variability of SLA, total leaf area, aboveground biomass and leaf mass fraction among forest understory plant species, (2) to compare SLA obtained from the most widely used LEDA database (Kleyer and others 2008) with field measurements from our study, and (3) to assess the level of variability in all four traits studied, explained by species identity and collection site. We hypothesized that (1) due to our poor recognition of forest understory herbaceous species ecology, the level of intraspecific variability of SLA is underestimated in commonly used trait databases and that (2) the interspecific variability (caused by species identity) is greater than intraspecific variability (site-specific).

Methods
The study was conducted in randomly selected forest collection sites in Poland and Czechia. Authors of the study collected plants during their travels connected with study and work, with no assumptions about dataset assembly and species selection. The only requirement was to collect specimens in the forest, i.e., under the canopy of at least a group of trees. The study area represents a climate transitional between maritime and continental. The geographical range of the collection sites studied varied from 49.1338°N to 54.0412°N and from 14.7423°E to 22.7273°E ( Figure 1). Our plots represented mainly lowlands and highlands (elevation range 20-763 m a.s.l., with 75% of observations below 206 m a.s.l.). The collection sites were located in both deciduous and coniferous forests. The mean annual temperature in 1991-2010 was 8.07°C, and mean annual precipitation was 658.28 mm (GUS 2017). We studied 167 species from 112 genera and 40 families of herbaceous plants growing in various environmental condi-Leaf Traits and Aboveground Biomass Variability tions (Table S1). Our study covered species typical not only for forest ecosystems but also for meadows, forest edges and species without preferences to certain ecosystems types. However, all plants were collected from forests.
Field work was conducted in 275 locations throughout both countries and led to collection of 6543 individual plants. Plants were harvested during two growing seasons, from May 27, 2016, to September 5, 2017. We collected material during two growing seasons in order to collect more specimens. We are conscious of that different growing seasons could have an influence on the traits studied, as traits are climate-dependent (Gilliam 2007; Jagodziń ski and others 2016; Blonder and others 2018). However, on our collection sites there were no extreme weather events such as droughts, fires or floods during 2016 and 2017, which could have strongly influenced our results. The material collection took place in different environmental conditions (for example, under closed or low-density canopy, near forest trails or at the forest edge) but exclusively in temperate forest habitats. In each collection site, we aimed to harvest as many species as possible. We collected a minimum of six specimens, from microhabitats as different as possible as recommended by Cornelissen and others (2003). We identified the species in the field, and all identifications were verified in the laboratory. To exclude the influence of disturbance regime on plant functional traits, we decided to collect plant material only in managed forests with no information about previous disturbances or thinning (Dyderski and others 2016). The random design of the study may be both an advantage (independence of collection sites) and a drawback. Specimens were collected on different dates, phenological stages, climate and collection site conditions, which does not allow assessment of the impact of these factors on the traits studied. However, databases of functional traits (for example, LEDA or Bien) have also been compiled from all available observations, so that intraspecific variability in the impact of these factors on traits is probably unevenly characterized among species. For that reason, our study design resembled the process of development of trait databases.
Each plant was cut, divided into components (leaves, stem, and reproductive organs) and then dried in an oven with forced air circulation at 65°C (ULE 600 and UF450, Memmert GmbH + Co. KG, Germany), to a constant mass. Then, leaves were scanned using WinFOLIA 2013 PRO software (Regent Instruments Inc., Quebec, Canada) to measure their area. We measured only a representative subset of mature leaves. Seedlings were excluded from our study. All components were weighed using BP 210 S (Sartorius, Gö ttingen, Germany) and Mettler Toledo PG 1003-S scales, with an accuracy of 0.001 g. We calculated total aboveground biomass (sum of masses of the components; g), specific leaf area (SLA, leaf area/leaf dry biomass; cm 2 g -1 ), total leaf area (SLA 9 total leaf biomass per plant; cm 2 ) and leaf mass fraction (leaf dry biomass/total aboveground biomass; dimensionless).

Data Analysis
To assess the impact of species identity, habitat and collection site on SLA, total aboveground biomass, total leaf area and leaf mass fraction, we performed linear mixed models implemented in the lme4::lmer() function (Bates and others 2015). In each model we used collection site as a random effect, and then we performed four models for each trait studied, changing the fixed effects. As fixed effects we used species identity, habitat (coniferous fresh, coniferous wet, deciduous fresh, deciduous wet) and the interaction between species and habitat. Coniferous habitats refer to less-fertile collection sites, while deciduous-to more fertile, fresh-to sites with low and moderate moisture, and wet-to sites with high moisture and wetlands (Table S2). We also performed a null model (intercept only), to assess how adding the independent variables to the model will increase model quality, expressed by Akaike's information criterion (AIC). To assess the impact of fixed effects, we used Tukey's post hoc test using Bonferroni's correction for multiple comparisons. To evaluate the impact of fixed and random effects we calculated marginal (R m 2 ) and conditional (R c 2 ) coefficients of determination using the MuMIn::r.squaredGLMM() function (Nakagawa and Schielzeth 2013;Bartoń 2018). R m 2 refers to the variability explained by fixed effects, and R c 2 -both by fixed and random effects. Thus, this approach allowed us to exclude and quantify site-specific effects. To compare data obtained in our study with data from the LEDA database, we used the approach of Cordlandwehr and others (2013)-observed versus database SLA plot with a 1:1 line. To assess the differences between mean values provided by the LEDA database (Kleyer and others 2008), we used Student's t-tests using the stats::t.test() function. To decrease the likelihood of incorrectly rejecting a null hypothesis, we applied Bonferroni's corrections to p values, i.e., we assumed statistically significant differences when p = 0.00032, that is, 0.05/number of species compared (156). We decided to use the LEDA database as it is the only open database that covers traits of species studied by us in a geographically similar area. All mean values are followed by standard error (SE). All analyses were conducted in R software (R Core Team 2018).

Comparing SLA Field Measurements and Values from LEDA Database
For 85 of the 156 species studied, the differences of mean species SLA were greater than 20% compared to the LEDA database (Kleyer and others 2008), and for 23 species-greater than 50% (Figure 4, Table S7). We found statistically significant differences (p < 0.00032) between mean SLA from the LEDA database and our study for 89 of 156 species studied. Among these 89 species, in case of Leaf Traits and Aboveground Biomass Variability 22 species relative differences between mean SLA from the LEDA database and our study were lower than 20% (Table S7). The minimum SLA was lower for 79, and the maximum SLA was higher for 60, of the 108 species from LEDA for which SLA ranges were provided. For 78 species, the differences of minimum species SLA were greater than 20%, and for 38 species-greater than 50%. For 68 of 109 species studied, the differences of maximum species SLA were greater than 20% compared to the LEDA database (Kleyer and others 2008), and for 29 species-greater than 50%. For 96 of the species studied, we collected more samples (from 1 up to 221 specimens) for SLA analysis than in the LEDA database.

Variability of Studied Traits Explained by Species Identity and Collection Site
For all traits studied, the best fit models were those including the interaction between species and habitat (Table 1). For the SLA model, 58% of the variability was explained by all the fixed effects together and 82% by both fixed and random effects (Tables S8-S11). For the leaf mass fraction model, 54% of the variability was explained by fixed effects, and 89% by both fixed and random effects (Tables S12-S15). For the total leaf area model, 36% of the variability was explained by fixed effects and 74% by both fixed and random effects (Tables S16-19). For total aboveground biomass, the best model explained 37% of variability by fixed effects, and 75% by both fixed and random effects (Tables S20-S23). However, in all traits studied the level of variability explained by fixed effects between the best models and models including only species identity ranged from only 0.045 (leaf mass fraction) to 0.100 (biomass) (Figure 5). Moreover, for all the traits studied, the models including the interaction of species and habitat explained the highest amount of variability.

DISCUSSION
What are the Sources of Variability?
In our study, we found higher inter-than intraspecific variability of the traits analyzed. Variability of the traits studied may result from Figure 4. Comparison of SLA obtained in our study with SLA from the LEDA database. The diagonal (1:1) line shows where values are the same. The labels present shortened species names, created by joining the first three letters of genus and species names, for example, Imp_par-Impatiens parviflora. To increase the figure's readability only species with a difference in SLA greater than 50% between the two sources compared were labeled. Vertical and horizontal lines show the range (minimum and maximum) of SLA values obtained from the two sources.  Leaf Traits and Aboveground Biomass Variability different factors. As there are numerous possible sources of variability in trait values, we aimed to group the main varying factors into three groups: environmental, biological and geographical factors. Environmental sources of variability in trait values are mainly connected with the characteristics of the collection sites. The features of the collection site may be expressed through stand or species density, stand age, dominant tree species in a stand or position in the forest (that is, core or edge). However, they are mainly connected with one primary factor which is the availability of light (Knight and others 2008;Niinemets 2010;Paź and others 2018). Moreover, light availability can strongly modify SLA, total leaf area and specimen biomass (Jagodziń ski and others 2016). However, studies of functional traits of seedlings of seven tree species revealed that the species identity caused much more variability than light availability (Dyderski and Jagodziń ski 2019). Wide ranges of trait values may result also from collection site differences, differentiated habitats where the plants grow, the phase of the growing season (for example, early, middle, or late) and from year-to-year variation in growing season, and so on ( As one of the most significant biological factors, life strategy influences the variability of parameters studied through functional traits such as height, leaf area, leaf mass per area, leaf nitrogen content per unit mass, seed mass and stem specific density (Díaz and others 2016). For example, at the beginning of the growth season the mean SLA of some species shows an increasing trend (Jagodziń ski and others 2016). On the other hand, with the aging of leaves more carbon-rich compounds occur in their tissues, which makes the tissue thicker and leads to reduction in SLA (Reich and others 1991). Moreover, plants of shaded habitats tend to produce leaves of high SLA, as in their case the growth limiting factor is lack of light (Planchais and Sinoquet 1998;Milla and others 2008;Wyka and others 2012).
The geographical factors, which determine leaf size at the global scale, are amplitudes of leaf and air temperatures during day and night (Wright and others 2017). As total leaf area strongly depends on the geographical latitude, it shows the lowest values for plants growing near the Arctic and Antarctic Circles, and the highest values for equatorial plants (Wright and others 2017). Because the collection of plant specimens in this study took place in the climate transition between maritime and conti-nental, we assumed that the results of leaf size obtained in our study would be determined mainly by the geographical factor. An additional confirmation of that hypothesis was the latitudinal pattern presented by leaves of small area (for example, Vaccinium vitis-idaea, Vaccinium uliginosum) on habitats typical of the boreal climate zone, and leaves of large area (for example, Salvia glutinosa, Scrophularia nodosa) on habitats typical of the temperate zone.
As the impact of each group of factors described is strongly connected with one another, the following observations can be made. Plants of high SLA and low biomass values were mainly connected with deciduous forests of high moisture, as less light availability and higher nitrogen availability allows for thinner leaf blades. On the other hand, the biomass was generally higher in coniferous forests of high moisture, as they relate to high light availability. This may be justified generally by the appearance of the Ericaceae family of high biomass and low SLA values in those habitats. Therefore, as some species occurred on strictly determined habitats (for example, only in coniferous or only in deciduous forests), this had a significant influence on our results. Therefore, the fact that the models including the interaction of species and habitat explained the highest amount of variability may be justified by the ecological preferences of the species studied.

Differences Between SLA Field Measurements and LEDA Database
We found that 57% out of the species studied had significant differences compared to the LEDA database. Moreover, 15% of the species studied differed by more than 50% from the values contained in the LEDA database. This could result from the limited study area (region-scale), as the collection sites were not evenly distributed across the wide ranges of species' occurrences. However, our sampling effort included relatively high habitat diversity-for 62% out of the species studied we collected a higher number of specimens than are included in the LEDA database. Our findings highlight that thorough consideration of database composition must be taken during their use, as the values presented are significantly influenced by the methodology of sampling.

Intra-and Interspecific Variability of Trait Values Differ Among Spatial Scales
We found that SLA and LMF were driven mainly by species identity and that TLA and total aboveground biomass were mainly site-specific. Numerous researchers have shown that the values of traits of different plant species are site-specific and differ among gradients or between different functional groups, producing a functional continuum instead of a clear-cut classification of species (Albert and others 2010a; Jung and others 2010; Mayfield and others 2013;Popović and others 2016;Helsen and others 2018). Significant complex plastic responses of different plant species must be also taken into consideration while analyzing ecological processes at a community level (Kichenin and others 2013). Intraspecific trait variation sources are now the subject of statistical modeling, and assessing the level of intraspecific variability of species might be crucial in some studies (Moran and others 2016). Nevertheless, it is the scale of study that determines the appropriate type of trait data that should be used. For that reason, this decision must always be taken when adapting the methodology to the aim of the survey (Albert and others 2010b). A solution containing a trade-off between sampling effort and conducting research on a wide study area might be the application of recently developed tools integrating trait variability at both intra-and interspecific levels (Violle and others 2012), although some cases require solutions. For example, in the case of testing a physiological hypothesis, it is obvious that direct trait measurements are required, as the intraspecific trait variability may be crucial to properly quantify the correlations studied. Moreover, a thorough understanding of intraspecific trait variability can also contribute to the explanation of ecological processes at the community to global levels (Jung and others 2010; Moran and others 2016; Siefert and others 2015). On the other hand, running ecological surveys of wide spatial scale at the community or greater levels justifies the use of mean trait values from databases, as we showed that the impact of collection site and habitat is much lower than the significance of species identity. It is also worth noting that intraspecific trait variability tends to be lower when comparing organ-level (and not whole-plant) traits, such as leaf-connected traits in this study (Siefert and others 2015), which implies that direct measurements might be dispensable in wide-scale studies. However, all functional traits can differ across a large latitudinal gradient (De Frenne andothers 2011, 2013;Lemke and others 2015). Therefore, using a single trait value and omitting the importance of intraspecific variability can obscure a significant part of the full characteristics of a given species (Albert and others 2010b).

Ecological Implications
Our study revealed higher interspecific than intraspecific variability of the traits studied. This result justifies the use of mean values of SLA, total leaf area, aboveground biomass and leaf mass fraction from databases in studies at the habitat scale and wider spatial units. The analysis of functional traits and functional variability based on calculated mean SLA values can serve for further verification of SLA worthiness as a proxy of correlated traits, such as total leaf area, nitrogen content, seed mass or plant height (Díaz and others 2016). Wide application of SLA as a feature coupled with other parameters by numerous studies (Dyderski and others 2016;Kattenborn and others 2017; Czortek and others 2018) determines the need for thorough recognition and understanding of its variability. We also found that the LEDA database (Kleyer and others 2008), as a commonly used source of mean SLA values, often presents values obtained using only a few specimens. In the context of the given shortcomings, this is a significant limitation. Nevertheless, our study showed that mean trait values depend on site-specific factors to only a small degree.
Our results have increased knowledge of the variability of the functional traits analyzed, which is an important concern for improving the fits of primary net production or carbon sequestration models to the specifics of forest ecosystems (Pompa-García and others 2017; Wright and others 2017; Jagodziń ski and others 2018). Functional traits are also used in nature conservation to assess the level of human impact on the environment, caused, for example, by land-use change, biological invasions or forest management (te Beest and others 2015; Dyderski and others 2016; Czortek and others 2018). Use of trait values from databases not acknowledging intraspecific variability is biased by the uncertainty about this variability. For that reason, our results reduced uncertainty of analysis using mean trait values conducted in order to assess human impacts on ecosystems. Thus, our study might support the assumption that level of intraspecific variability of functional traits is lower than interspecific variability.

AC KNOWLEDGEMENT S
The study was partially financed by the Institute of Dendrology, Polish Academy of Sciences, Kó rnik, Poland. We are grateful to Dr. Lee E. Frelich (Department of Forest Resources, University of Minnesota, USA) for linguistic revision of the manuscript. We are thankful to two anonymous reviewers for valuable and thorough comments on the earlier draft of the manuscript.

OP E N A C C ES S
This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4 .0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.