Plasticity in above- and belowground resource acquisition traits in response to single and multiple environmental factors in three tree species

Functional trait plasticity is a major component of plant adjustment to environmental stresses. Here, we explore how multiple local environmental gradients in resources required by plants (light, water, and nutrients) and soil disturbance together influence the direction and amplitude of intraspecific changes in leaf and fine root traits that facilitate capture of these resources. We measured population-level analogous above- and belowground traits related to resource acquisition, i.e. “specific leaf area”–“specific root length” (SLA–SRL), and leaf and root N, P, and dry matter content (DMC), on three dominant understory tree species with contrasting carbon and nutrient economics across 15 plots in a temperate forest influenced by burrowing seabirds. We observed similar responses of the three species to the same single environmental influences, but partially species-specific responses to combinations of influences. The strength of intraspecific above- and belowground trait responses appeared unrelated to species resource acquisition strategy. Finally, most analogous leaf and root traits (SLA vs. SRL, and leaf versus root P and DMC) were controlled by contrasting environmental influences. The decoupled responses of above- and belowground traits to these multiple environmental factors together with partially species-specific adjustments suggest complex responses of plant communities to environmental changes, and potentially contrasting feedbacks of plant traits with ecosystem properties. We demonstrate that despite the growing evidence for broadly consistent resource-acquisition strategies at the whole plant level among species, plants also show partially decoupled, finely tuned strategies between above- and belowground parts at the intraspecific level in response to their environment. This decoupling within species suggests a need for many species-centred ecological theories on how plants respond to their environments (e.g. competitive/stress-tolerant/ruderal and response-effect trait frameworks) to be adapted to account for distinct plant-environment interactions among distinct individuals of the same species and parts of the same individual.


Introduction
All plants obey the same trade-offs that limit investments of the resources they acquire to different parts of cells, different tissues, and different organs (Herms and Mattson 1992). As a consequence of these and other constraints, there are broadly consistent interspecific pat-terns of resource allocation between plant organs, in terms of total biomass (Enquist and Niklas 2002) and in cell and tissue structure and chemistry (Kerkhoff et al. 2006;Freschet et al. 2010). Nevertheless, environmental conditions have potentially strong impacts on quantitative (Poorter et al. 2012) and qualitative (Liu et al. 2010) patterns of allocation, by driving local selection in plant genotypes and adjustments through plant phenotypic plasticity. At the species level, the relative availability of different resources can trigger adaptive and non-adaptive physiological changes in plants, which can mitigate the constraints imposed by the most limiting resources (Chapin et al. 1987;Van Kleunen and Fischer 2005;Valladares et al. 2007). Plants can therefore adjust biomass allocation to roots or shoots (i.e. quantitative adjustment) according to whether the most limiting resource is above-or belowground (Shipley and Meziane 2002;Poorter et al. 2012). Plants can also adjust the morphology and efficiency of their tissues (i.e. qualitative adjustment) to increase uptake of the most limiting resource (Ryser and Eek 2000;Hill et al. 2006) to the extent that these adjustments might be more important than changes in total mass allocation as a response to abiotic stresses (Poorter et al. 2012).
Qualitative intraspecific adjustments to abiotic stresses are widespread among plants, and are manifest in many plant traits (Chapin 1991). Aboveground, decreasing light availability commonly drives an increase in leaf area per unit mass (specific leaf area, SLA; Rijkers et al. 2000;Evans and Poorter 2001), i.e. there is an increase in the area of light capture for a constant amount of resource invested. Also, decreasing partial pressure in CO 2 leads to higher leaf nitrogen concentrations (LNC), which is indicative of the leaf carboxylation capacity (i.e. concentration of CO 2 -fixing protein Rubisco in leaf) per unit mass invested (Ellsworth et al. 2004). Belowground, increasing nutrient limitation generally drives intraspecific increases in specific root length (SRL) (Hill et al. 2006;Ostonen et al. 2007), which represents the length of root potentially able to explore soil per unit mass invested. Together these examples illustrate how changes in the availability of any one resource, either above-or belowground, trigger changes in functional traits that determine the ability of plants to acquire that resource. However, several other environmental factors can also influence these responses (see Poorter et al. 2009 for an overview of factors influencing SLA). Under natural conditions, plants are exposed to multiple stresses and compete both for above-and belowground resources simultaneously. Furthermore, the acquisition of any one resource (e.g. carbon) generally relies on several others (e.g. light, water, and nutrients). For instance, the effect of light stress on SLA can shift from null to positive along a gradient of nitrogen availability (Meziane and Shipley 1999). Therefore, multiple environmental stresses can drive complex integration of responses at the plant level (Chapin 1991;de Kroon et al. 2005), which are potentially accountable for distinct responses across plant species with contrasting physiologies.
An important component of plant responses to the environment, other than the direction of intraspecific trait adjustments, is the strength of the response (Valladares et al. 2006). This magnitude of response can be quantified as changes in plant biomass, seed production, or any plant trait related to fitness. A high magnitude of response in those traits that determine fitness can potentially help a species maintain fitness in stressful environments and/or increase fitness in favorable environments (Richards et al. 2006). However, strong responses also come with costs and limits to fitness, such as the resources needed to generate responses, intrinsic genetic costs owing to pleiotropy, gene linkages and epistasis, or unstable plant development when environmental signals are unreliable (DeWitt et al. 1998;Van Kleunen and Fischer 2005). Although the adaptive role of intraspecific plant responses has yet to be generally demonstrated, it is likely to play an important role in determining plant fitness (Sultan 2000;Valladares et al. 2006). Several hypotheses have been proposed regarding its connections to various aspects of plant strategies (see Lavorel et al. 2009). In particular, resource-acquisitive species have been hypothesized to express a greater magnitude of intraspecific responses than resource-conservative ones because they generally appear better at exploiting varying levels of resource availability (Crick and Grime 1987;Grassein et al. 2010). However, the extent to which plant intraspecific responses are related to plant economics has received little explicit attention (Grassein et al. 2010) except in the context of plant invasions (Richards et al. 2006;Funk 2008).
In this study, we explore how multiple environmental gradients associated with light (carbon), water, and nutrient resource availability, as well as disturbance, collectively drive intraspecific changes in leaf and root traits related to the capture of these same resources. Firstly, we hypothesize that three tree species competing for the same spatial niche, but with contrasting carbon and nutrient economics, will have broadly similar responses to any single environmental factor, but will have stronger and more specific responses to combinations of several environmental factors. This is because the response of each species results from a complex integration of environmental factors at the whole-plant level (Chapin 1991;de Kroon et al. 2005). Secondly, we hypothesize that the magnitude of intraspecific responses of these species, both above-and belowground, will be related to the extent to which their economies are acquisitive or conservative. Thirdly, we hypothesize that, despite the potentially widespread coordination in leaf and root economics across species (Reich et al. 2008;Liu et al. 2010), there will be uncoordinated shifts in nutrient and carbon economics above-versus belowground within species, because leaf and root traits (associated with the acquisition of aboveand belowground resources respectively) should respond differently to (orthogonally) varying levels of resources above-versus belowground. We test these hypotheses by measuring a set of above-and belowground plant traits representative of plant carbon and nutrient economics on three tree species at the population level within the same ecosystem along strong gradients of light, soil nutrient concentrations, moisture, and disturbance.

Study system
Korapuki is a small, isolated, forested 18 ha island of volcanic origin and relatively steep terrain in the Mercury Islands, off the NE coast of the North Island of New Zealand (36°39′S, 175°50′E, maximum elevation 81 m). The climate is warm temperate with mean monthly temperatures ranging from 14°C in June to 19.5°C in January and annual precipitation of c. 2000 mm. The vegetation cover is generally dense with a tall canopy of Metrosideros excelsa (Myrtaceae) over a sub-canopy of Coprosma macrocarpa (Rubiaceae), Melicytus ramiflorus (Violaceae), and Pittosporum crassifolium (Pittosporaceae)hereafter Coprosma, Melicytus, and Pittosporumwhich is typical of vegetation on warm temperate islands in New Zealand (Atkinson 2004). A variety of burrowing seabirds (primarily Pterodroma macroptera gouldi, but also Pelecanoides urinatrix, Pterodroma pycrofti, Puffinus assimilis haurakiensis, and Puffinus gavia, all Procellariiformes) breed on the island, some of which have become more abundant after the eradication of introduced rabbits (Orcytolagus cuniculus) and Pacific rats (Rattus exulans) in 1986-1987 (Towns 2002). Burrowing by these seabirds causes considerable soil disturbance (Fukami et al. 2006), and as a consequence, on this island ground cover ranges from almost bare soil in heavily burrowed areas to light understorey vegetation and abundant leaf litter in areas with little seabird influence.
Square plots of 10 9 10 m were randomly assigned to the island using ArcGIS 10 software (ESRI, Redlands, CA, USA) from which fifteen plots representing the whole range of seabird burrowing activity (from 0 to 94 burrows per plot) were selected. Our plot selection excluded areas of the island such as rocky shoreline or those that were inaccessible because of topographical constraints. The distance between plots varied between c. 50 and 600 m. The three dominant sub-canopy tree species, Coprosma, Melicytus and Pittosporum were present on all plotsbut only as mature trees in 14, 12, and 15 plots respectivelyand were therefore the focus of this study. The relatively homogeneous cover of these three species across the island creates an ideal system for studying the impact of spatial environmental variations on plant plasticity. Selection of tree species that grow to a similar size has the advantage of limiting potential differences in species responses related to plant size (e.g. the capacity of plants to access heterogeneously distributed patches of light or soil resources; Casper and Jackson 1997). Nevertheless, these three species have contrasting carbon and nutrient economics, as evidenced by contrasting values of leaf and root traits representative of plant economics (see Wardle et al. 2009; this study), and therefore potentially different capacities to overcome resource stresses and disturbances, i.e. potentially species-specific responses. Furthermore, Korapuki is characterized by strong belowground spatial heterogeneity and variation in upper canopy structure over even small spatial scales, making it ideal for testing questions about spatial resource heterogeneity. The presence of seabirds is likely to create strong gradients of disturbances and soil fertility across the island (Bancroft et al. 2005;Mulder et al. 2011) and the occurrence of gaps in the upper canopy leads to considerable spatial variation in the light transmittance to the sub-canopy species. Finally, all sampled tree populations on Korapuki (which is small and isolated) have presumably regrown from the same tree population which has developed since the 1950s after earlier deforestation of the island (Towns 2002). Although new genetic material may have been brought in from distant islands, most passerines involved in the dispersal of fruits from the focal trees in our study disperse most seeds over short distances (tens of meters, Stansbury 2001;Williams 2006). Therefore, differences between tree populations that result from genotypic adaptation to local conditions on the island are presumably very small, and the trait variation reported here is therefore assumed almost entirely to reflect phenotypic plant responses to the environment.

Environmental measures
In July 2010, the number of seabird burrow entrances and eggs (one per breeding pair) were counted per plot. These two variables were strongly positively related across plots (R 2 = 0.87; P < 0.001) and we therefore used only burrow density as a measure of seabird impact. The density of burrows per area is a relative measure both of the nutrient transfer to the soil from the seabirds and disturbance by the seabirds to plant roots resulting from excavation and upkeep of the burrows (Roberts et al. 2007;Wardle et al. 2009). Increasing burrow density is also likely to have some impact on the water content and compaction of soils by favoring water drainage and aeration (Bancroft et al. 2005).
Canopy openness was estimated from two hemispherical photographs per plot (c. 1.3 m height) using Gap ª 2013 The Authors. Published by Blackwell Publishing Ltd.
Light Analyzer 2.0 software (Simon Fraser University, Burnaby, Canada). It therefore represents an index of light availability to leaves of the lower sub-canopy (on average from 2 to 4 m in height), including those of our three focal tree species.
Four 10 9 10 9 10 cm soil pits were sampled per plot in March 2010, pooled, homogenized and used for further soil analysis. The mass and volume of this soil was used to calculate bulk density. Soil moisture content was measured on a sub-sample dried to constant weight at 60°C. Since weather was warm and dry with no precipitation during the entire sampling period, soil moisture content was used here as an index of water availability for plants. The two main form of N available for plants, i.e. NO 3 -N and NH 4 -N, were extracted with KCl solution and measured by auto-analyser (Lachat Instruments Division, Zellweger Analytics Inc., Milwaukee, WI, USA). A measure of available phosphorus for plants, Olsen-P, was measured using sodium bicarbonate digest.

Plant traits
In November 2011, Coprosma, Melicytus and Pittosporum were sampled from each plot for shade leaves and fine superficial roots sampled whenever possible from a minimum of ten distinct plant individuals. Fifty healthy, fully expanded leaves (presumably from the current year) were sampled per species per plot from the inside of the subcanopy (between c. 2-3 m in height). Ten to twenty roots per species per plot were tracked down from the tree stem to the fine root tips and dug out. All leaf and root material was kept cool and wet until cleaning and further processing in the lab. All root branches below 1 mm in diameter, except explorative roots, were clipped off, thoroughly cleaned from residual soil and taken as "fine roots." For each species and plot, all 50 leaves and all fine root material was used to estimate specific leaf area (SLA, mm 2 /mg) and specific root length (SRL, m/mg) as well as leaf and fine root dry matter content (DMC, fresh to dry weight ratio, mg/g). Measurements followed Cornelissen et al. (2003). All leaves and fine roots were then coarsely ground and homogenized manually before taking a subsample for total N and P concentration analysis using acid digest and colorimetric methods of Blakemore et al. (1987).

Data analysis
The single and multiple influences of four environmental parameters (light availability, soil nutrient availability, soil moisture availability, and seabird burrow density) on intraspecific variation in leaf and root functional traits (SLA-SRL, DMC, N and P) were assessed separately for each plant species using simple and multiple linear regressions. In order to maximize the statistical power of our tests and minimize multicollinearity, two highly correlated environmental parameters were transformed into one summarizing index using a principal component analysis (PCA): "soil nutrient availability" was thus estimated as the first axis of a PCA that summarizes 89% of variation in available N (nitrate + ammonium) and available P (Olsen-P). Seabird burrowing activity was included in multiple linear regressions because it represents the impact of disturbances caused by seabird burrowing activity in addition to the effects of seabirds on soil nutrients, so could potentially explain variability that 'soil nutrient availability' cannot. Seabird burrow density was log 10transformed for all single and multiple regressions concerned to comply with normality assumptions. For each species 9 trait combination, the most relevant models were identified using Akaike's Information Criterion (AIC). Models displaying AIC lowest AIC + 1 were considered equivalent. Results of regressions with four parameters were not shown because of the low number of observations per parameter and because their AIC never fell below these with less parameters.
To illustrate the effect of multiple contrasting environmental influences on the intraspecific variation in wholeplant economic strategies of our three tree species, we performed a PCA on all leaf and root traits of all species simultaneously. We used t-tests among pairs of species to determine whether the ranking of these three species along the first axis of the PCA, which represented a "plant economics spectrum" , changed across environmental influences (i.e. across plots).
In order to compare the magnitude of plastic responses among tree species and organs to the same environmental variation, we focused on intraspecific trait responses to multiple environmental influences rather than to single environmental factor or sets of environmental factors (Valladares et al. 2006). This is because intraspecific shifts in plant trait values were generally co-determined by several environmental parameters which often differed across species. Consequently, measuring trait responses of different species to the same set of (subjectively chosen) environmental factors would make little sense in this context. Specifically, for each trait 9 species combination, we transformed trait values at each plot into percentages of deviation from the mean trait value across all plots. Then, for each trait, we used two-sample Fisher's F-tests to test for differences between species in the homogeneity of variance in these percentages of deviation from the mean trait value. For each species, we also used Fisher's F-tests in the same way to assess the homogeneity of intraspecific variance between functionally analogous leaf and root traits (i.e. SLA vs. SRL; root vs. leaf DMC, N and P).
Expressing the magnitude of intraspecific trait responses as a relative deviation from the species mean trait value is equivalent to using coefficients of variation (CV) and allows comparisons of species' variance independently of the mean plant trait value.
Standardized major axes regressions were used to test the coordination between functionally equivalent leaf and root traits across the 15 plots for each species and, when significant, slope comparison procedures were performed between species (SMATR-package; Warton et al. 2006).

Results
Among all possible pairs of four environmental parameters (see Table 1 for an outline of parameters' ranges, means, and variabilities), only soil nutrient availability and seabird burrow density were significantly positively related (R 2 = 0.45; P < 0.01). Light availability was largely unrelated to soil moisture content, soil nutrient availability, and seabird disturbance (R 2 = 0.18, 0.10, 0.20, respectively; P ! 0.10). Soil moisture content was unrelated to seabird disturbance (R 2 = 0.11; P ! 0.10) and only marginally negatively related to soil nutrient availability (R 2 = 0.25; P = 0.06).

Multiple environmental factors codetermine intraspecific trait responses
The PCA primary axis scores for the three species with regard to their traits differed greatly, although there was some overlap between Coprosma and each of the other two species along the plant economics spectrum axis (Fig. 1). The species' ranking along this axis was strictly conserved across the 15 plots, indicating that there were broadly consistent influences of environmental parameters on multiple traits of the three distinct tree species.
More specifically, several models that included either single or multiple environmental factors often explained more than half of the total trait variance across the plots for all three species (Tables 2 and 3; Tables S1 and S2). For most plant traits, although we observed largely consistent trends across all three tree species, the best models often differed at least slightly between species. This was particularly true for SLA, SRL, leaf and root N concentrations, root DMC and root P concentration, but less so for leaf DMC and leaf P concentration.
Among the models with the lowest AIC (i.e. lowest AIC + 1), population level SLA across plots was best fitted by models using seabird disturbance, soil moisture content and light availability for Pittosporum, disturbance and light availability for Coprosma, and disturbance, soil moisture content and soil nutrient availability for Melicy-tus (Table 2). However, when single parameter predictions were considered, disturbance, light and soil moisture on their own consistently explained significant amounts of variation in SLA. While intraspecific variation in SRL was best fitted by models using disturbance, soil nutrient availability and light for both Pittosporum and Melicytus, it was best fitted by models that integrated soil nutrient availability and moisture for Coprosma (Table 2). Also, the sign of disturbance and light parameters of the models appeared reversed between Pittosporum and Melicytus. Thus, in all three cases only soil nutrient availability had a consistent (positive) influence on SRL. When considered alone it explained a large proportion of the total variability in SRL across plots, i.e. 49, 53, and 42% for Pittosporum, Coprosma, and Melicytus, respectively.   For all three species, soil nutrient availability alone explained a substantial proportion of the variability across plots in leaf and root N concentrations (Table 3). However, Pittosporum LNC was most strongly influenced by a combination of soil nutrient availability, light and disturbance, while LNC for Coprosma and Melicytus were most responsive to soil nutrient availability and light only. Intraspecific variability for root N concentration (RNC) was also largely driven by nutrients and light, but while Pittosporum and Coprosma responded most strongly to these in combination with disturbance, Melicytus responded most strongly to them in combination with soil moisture (Table 3).
Intraspecific variability in leaf dry matter content (LDMC) was relatively well fitted by models containing disturbance alone or including both disturbance and light for Coprosma and Melicytus, but no models were satisfactory for Pittosporum (Table S1). In contrast, while largely different models best fitted intraspecific variability in root dry matter content (RDMC) across species, soil nutrient availability emerged as a consistent driver across all three species and explained 64, 63, and 58% of the variability in RDMC of Pittosporum, Coprosma, and Melicytus, respectively (Table S1).
Finally, while a model that included light, soil nutrient availability and moisture explained up to 59% in the intraspecific variability in leaf P concentration (LPC) of Melicytus, none of the environmental parameters could adequately explain the intraspecific variability in LPC for Pittosporum and Coprosma (Table S1). In contrast, the intraspecific variation in root P concentration (RPC) across species appeared mostly driven by disturbance and soil nutrient availability although, depending on the species, either soil moisture or light appeared to slightly improve model fits (higher R 2 but similar AIC).

A broadly similar magnitude of intraspecific trait responses across organs and species
Only leaf and root DMC followed the hypothesized trend of increasing trait variance from Pittosporum to Coprosma to Melicytus, and the only difference in trait variance among species that was statistically significant was for LDMC between Pittosporum and both Coprosma and Melicytus (P = 0.01 in both cases; Fig. 2). Furthermore, variances of six traits, including RDMC, leaf and root N and P concentrations, and SLA, were not significantly different among the three tree species. Also, the intraspecific variance in SRL was significantly higher for Coprosma than for both Pittosporum and Melicytus (P = 0.01 in both cases). At the species level, mostly similar variances were observed between functionally equivalent leaf and root traits, although in three instances root traits showed higher variance than their leaf trait counterparts. Thus, Pittosporum displayed significantly higher variance in root than leaf DMC (P = 0.04) and Coprosma showed significantly higher variance in root than leaf N concentration (P = 0.04) and in SRL than SLA (P < 0.01).
Orthogonal responses of analogous leaf and root traits to multiple environmental influences As described above, for each of the three species, different environmental parameters appeared to drive intraspecific shifts in analogous above-and belowground traits, such as SLA versus SRL (Table 2), LDMC versus RDMC (Table  S1), and LPC versus RPC (Table S2). In contrast, LNC and RNC were driven by a similar combination of parameters, notably nutrient and light availability ( Table 3). As a consequence, while intraspecific variation in LNC and RNC were significantly related (Fig. 3b), all other pairs of leaf and root traits were not (Fig. 3a, c, d). The slope of LNC-RNC relationship for Melicytus was significantly smaller than this for the two other species (P < 0.01 in both cases).

Discussion
Multiple environmental factors often codetermine intraspecific trait responses and drive partially species-specific responses Most traits of the three tree species showed similar responses to the same environmental factors when factors were considered separately (e.g. Wahl et al. 2001;Rozendaal et al. 2006), and, consistent with our first hypothesis, combinations of several factors often explained higher intraspecific trait variation than single factors. For instance, light availability is generally considered as the main driver of intraspecific shifts in SLA (Evans and Poorter 2001;Rozendaal et al. 2006). However, for all three tree species we observed a strong influence of soil disturbance caused by seabird burrowing and soil moisture content (see also Poorter et al. 2009), which together explained more variance than the otherwise significant influence of light availability alone. Also, even though the variance in leaf and root N concentrations explained by light availability was negligible for all three species, the combination of light availability and soil nutrient availability had a much stronger effect on leaf and fine root N concentrations than did soil nutrient availability alone. These results suggest that considering phenotypic responses to single environmental parameters, as done by all but a handful of studies, is inadequate to predict intraspecific plant trait responses to multiple co-occurring environmental factors (see also Meziane and Shipley 1999;Wahl et al. 2001;Valladares et al. 2007;Poorter et al. 2009). Two patterns of trait response to multiple environmental factors differed substantially from those generally observed in studies of single environmental factors. First, in addition to the positive impact of soil nutrient availability on leaf and root N concentrations (Xia and Wan 2008), light availability appeared to co-determine the N concentrations of leaves and fine roots of our three subcanopy species. Increasing leaf N concentrations with decreasing light may correspond to an adaptive investment to improve light capture capacity in situations where carbon is the limiting resource (Evans and Poorter 2001). However, this trend can be more parsimoniously explained by higher carbon assimilation and thereby plant growth in higher light conditions, which potentially creates a N sink in growing plant parts and thereby decreases the 'luxury consumption' of N in both leaf and fine root tissues (Van Wijk et al. 2003). Second, in contrast to most previous studies (e.g. Hill et al. 2006;Ostonen et al. 2007), we observed decreasing SRL and increasing RDMC with decreasing soil nutrient availability (but see Ryser and Lambers 1995;Yano and Kume 2005). This trend may benefit plants by increasing their root lifespan in environments that provide a slower return on investment (Eissenstat et al. 2000;Ryser 2006) as observed across species of low and high resource environments (Ryser and Eek 2000;Freschet et al. 2010). Additionally, decreasing SRL may also represent an increase in fine root mycorrhizal associations (Eissenstat et al. 2000) in conditions of high nutrient stress.
Our results also revealed some species-specific responses to spatial variation in light, soil nutrient, and moisture availabilities, as well as root-damaging disturbance by seabirds. For instance, while RDMC decreased with soil nutrient availability for all three species, nutrient and light availability strongly co-determined RDMC in Pittosporum, and nutrient, light, and seabird burrowing co-determined RDMC in Melicytus. Our results, together with previous studies exploring gradients in soil nutrient availability (Lawrence 2003), soil moisture content (Cornwell and Ackerly 2009), and temperature (Albert et al. 2010), suggest that plant communities could exhibit complex responses to changing environmental conditions (Lavorel et al. 2009), which portends substantial uncertainty in trait-based projections of plant assemblages under scenarios of global change (Thuiller et al. 2008).
Does the strength of intraspecific plant trait responses relate to plant resource acquisition strategy?
We found broadly similar amplitudes of intraspecific plant trait responses across all three tree species to the same multiple environmental influences, irrespective of their nutrient and carbon economics ( Fig. 2; see also Funk 2008). This result did not provide support for our second hypothesis. Only leaf and root DMC displayed the expected trend of increasing strength of intraspecific response (although mostly non-significant) as the species became more resource-acquisitive, which is consistent with results for leaf DMC in alpine grassland communities (Lavorel et al. 2009;Grassein et al. 2010). Our results therefore suggest that the capacity for large intraspecific variation in traits might be equally important for both acquisitive and conservative species. Environmental heterogeneity is widespread within both fertile and unfertile sites and plants characteristic of both environments may have a selective benefit through exhibiting strong intraspecific trait responses. Studies featuring a greater number of species remain necessary to answer whether the magnitude of plant intraspecific responses relates to their resource economics. Also, our study system is unlikely to cover the entire range of environmental conditions that our focal species can occupy. Therefore, the general lack of relationship observed in our study between the strength of plant intraspecific response and plant economics does not exclude the possibility that between-species differences in plasticity could be observed with respect to the physiological limits of plastic responses (rather than the proportional strength of the response).

Multiple environmental influences drive decoupled plant resource acquisition strategies above-and belowground
In support of our third hypothesis, we observed uncoordinated shifts within species for three of four functionally analogous leaf and root traits, which resulted from con-trasting responses of leaf and root traits to the same multiple environmental gradients. Growing evidence suggests that there are broadly coordinated economics of both carbon and nutrients between plant above-and belowground organs at the between-species level (Wright and Westoby 1999;Freschet et al. 2010;Liu et al. 2010), most likely as a consequence of allometric, physiological, and ontogenetic constraints in plant evolutionary selection (Ackerly et al. 2000;Reich et al. 2003b;Kerkhoff et al. 2006). However, these relationships are not always apparent (Reich et al. 2003a;Craine et al. 2005;Tjoelker et al. 2005) and can vary in their slope and intercept across plant communities as a consequence of variation in climatic or edaphic parameters Wright et al. 2005;Liu et al. 2010). The large magnitude and orthogonal direction of intraspecific shifts in some analogous leaf and root trait values (SLA vs. SRL, leaf vs. root DMC and P; Fig. 3) that we observed at the local scale (see also Ryser and Eek 2000) suggests a major role for plant plasticity, as opposed to species turnover, in altering these trait relationships across environmental gradients.
Our findings therefore suggest that global relationships between (species mean) leaf and root traits could potentially occur but effectively disrupted by plant plastic responses to contrasting above-and belowground environmental heterogeneity occurring at multiple spatial scales (Ettema and Wardle 2002). These large partially decoupled plastic adjustments also suggest a need for many species-centred ecological theories on how plants respond to their environments (e.g. r-K strategy theory, competitive/stress-tolerator/ruderal theory, trait-based responseeffect frameworks) to be adapted to account for distinct plant-environment interactions among distinct individuals of the same species and parts of the same individual.

Conclusion
Our results highlight the potentially complex integration of multiple environmental influences by plants, stress the large potential for plastic adjustments both above and belowground and reveal that these adjustments can to some extent differ among coexisting plant species. They also suggest that while there may be broadly consistent strategies at the whole plant level among species (Grime 2001), there are also partially decoupled, finely tuned strategies between above-and belowground plant parts at the intraspecific level. These findings have a number of implications for plant community assembly, plant-soil feedbacks and ecosystem properties. For instance, they indicate that plants are capable of considerable flexibility in passing through multiple environmental filters that may exert contrasting constraints on above-and belowground organs (D ıaz et al. 1998), which may potentially impact on the resilience of ª 2013 The Authors. Published by Blackwell Publishing Ltd.
plant communities to environmental changes. Furthermore, our finding that trait responses to the same underlying environmental influences differed at least partially among coexisting species suggests potentially complex responses of plant communities to shifts in multiple environmental conditions. Also, given the role that plant species have as ecosystem drivers, the decoupled adjustments of leaf and root traits along gradients of environmental stresses may potentially contribute to contrasting feedbacks of litters from above-and belowground plant organs with soil properties (Wardle et al. 2004;Freschet et al. 2012). For all these reasons, our results stress the importance for ecologists to routinely consider the above-and belowground components of plants together and to recognize the major role of functional trait plasticity in driving plant responses to community assembly processes and the effects that plants have on ecosystem properties.

Supporting Information
Additional Supporting Information may be found in the online version of this article: Table S1. Linear models of single and multiple environmental parameters' fit of leaf and fine root dry matter content of three sub-canopy tree species. Table S2. Linear models of single and multiple environmental parameters' fit of leaf and fine root phosphorus content of three sub-canopy tree species.