Long-term exclusion of invasive ungulates alters tree recruitment and functional traits but not total forest carbon.

Forests are major carbon (C) sinks, but their ability to sequester C and thus mitigate climate change, varies with the environment, disturbance regime, and biotic interactions. Herbivory by invasive, nonnative ungulates can have profound ecosystem effects, yet its consequences for forest C stocks remain poorly understood. We determined the impact of invasive ungulates on C pools, both above- and belowground (to 30 cm), and on forest structure and diversity using 26 paired long-term (>20 years) ungulate exclosures and adjacent unfenced control plots located in native temperate rainforests across New Zealand, spanning 36-41° S. Total ecosystem C was similar between ungulate exclosure (299.93 ± 25.94 Mg C ha-1 ) and unfenced control (324.60 ± 38.39 Mg C ha-1 ) plots. Most (60%) variation in total ecosystem C was explained by the biomass of the largest tree (mean diameter at breast height [dbh]: 88 cm) within each plot. Ungulate exclusion increased the abundance and diversity of saplings and small trees (dbh ≥2.5, <10 cm) compared with unfenced controls, but these accounted for ~5% of total ecosystem C, demonstrating that a few, large trees dominate the total forest ecosystem C but are unaffected by invasive ungulates at a timescale of 20-50 years. However, changes in understory C pools, species composition, and functional diversity did occur following long-term ungulate exclusion. Our findings suggest that, although the removal of invasive herbivores may not affect total forest C at the decadal scale, major shifts in the diversity and composition of regenerating species will have longer term consequences for ecosystem processes and forest C.

Soil biotic and abiotic effects on seedling growth exhibit context-dependent interactions: evidence from a multi-country experiment on Pinus contorta invasion.

The success of invasive plants is influenced by many interacting factors, but evaluating multiple possible mechanisms of invasion success and elucidating the relative importance of abiotic and biotic drivers is challenging, and therefore rarely achieved. We used live, sterile or inoculated soil from different soil origins (native range and introduced range plantation; and invaded plots spanning three different countries) in a fully factorial design to simultaneously examine the influence of soil origin and soil abiotic and biotic factors on the growth of invasive Pinus contorta. Our results displayed significant context dependency in that certain soil abiotic conditions in the introduced ranges (soil nitrogen, phosphorus or carbon content) influenced responses to inoculation treatments. Our findings do not support the enemy release hypothesis or the enhanced mutualism hypothesis, as biota from native and plantation ranges promoted growth similarly. Instead, our results support the missed mutualism hypothesis, as biota from invasive ranges were the least beneficial for seedling growth. Our study provides a novel perspective on how variation in soil abiotic factors can influence plant-soil feedbacks for an invasive tree across broad biogeographical contexts.

Community- and trophic-level responses of soil nematodes to removal of a non-native tree at different stages of invasion.

Success of invasive non-native plant species management is usually measured as changes in the abundance of the invasive plant species or native plant species following invader management, but more complex trophic responses to invader removal are often ignored or assumed. Moreover, the effects of invader removal at different stages of the invasion process is rarely evaluated, despite a growing recognition that invader impacts are density or stage-dependent. Therefore, the effectiveness of invasive species management for restoring community structure and function across trophic levels remains poorly understood. We determined how soil nematode diversity and community composition respond to removal of the globally invasive tree species Pinus contorta at different stages of invasion by reanalysing and expanding an earlier study including uninvaded vegetation (seedlings removed continuously), early invader removal (saplings removed), late removal (trees removed), and no removal (invaded). These treatments allowed us to evaluate the stage-dependent belowground trophic responses to biological invasion and removal. We found that invaded plots had half the nematode taxa richness compared to uninvaded plots, and that tree invasion altered the overall composition of the nematode community. Differences in nematode community composition between uninvaded nematode communities and those under the tree removal strategy tended to dilute higher up the food chain, whereas the composition of uninvaded vs. sapling removal strategies did not differ significantly. Conversely, the composition of invaded compared to uninvaded nematode communities differed across all trophic levels, altering the community structure and function. Specifically, invaded communities were structurally simplified compared to uninvaded communities, and had a higher proportion of short life cycle nematodes, characteristic of disturbed environments. We demonstrate that a shift in management strategies for a globally invasive tree species from removing trees to earlier removal of saplings is needed for maintaining the composition and structure of soil nematode communities to resemble uninvaded conditions.

Microbial and abiotic controls on mineral-associated organic matter in soil profiles along an ecosystem gradient.

Formation of mineral-organic associations is a key process in the global carbon cycle. Recent concepts propose litter quality-controlled microbial assimilation and direct sorption processes as main factors in transferring carbon from plant litter into mineral-organic associations. We explored the pathways of the formation of mineral-associated organic matter (MOM) in soil profiles along a 120-ky ecosystem gradient that developed under humid climate from the retreating Franz Josef Glacier in New Zealand. We determined the stocks of particulate and mineral-associated carbon, the isotope signature and microbial decomposability of organic matter, and plant and microbial biomarkers (lignin phenols, amino sugars and acids) in MOM. Results revealed that litter quality had little effect on the accumulation of mineral-associated carbon and that plant-derived carbon bypassed microbial assimilation at all soil depths. Seemingly, MOM forms by sorption of microbial as well as plant-derived compounds to minerals. The MOM in carbon-saturated topsoil was characterized by the steady exchange of older for recent carbon, while subsoil MOM arises from retention of organic matter transported with percolating water. Overall, MOM formation is not monocausal but involves various mechanisms and processes, with reactive minerals being effective filters capable of erasing chemical differences in organic matter inputs.

Contrasting responses of soil nematode communities to native and non-native woody plant expansion.

Woody plant expansion into grasslands is widespread, driven by both successions to dominance by native woody species or invasion by non-native woody species. These shifts from grass- to woody-dominated systems also have profound effects on both above- and belowground communities and ecosystem processes. Woody-plant expansion should also alter the functional composition of the soil biota, including that of nematodes, which are major drivers of soil food-web structure and belowground processes, but such belowground impacts are poorly understood. We determined whether succession by a widespread native (Kunzea ericoides) and invasion by a non-native woody species (Pinus nigra) into tussock grasslands affect the composition of nematode functional guilds and the structure of nematode-based food webs. Although increasing dominance by woody species in both systems altered the functional guild composition of the nematode community, we found contrasting responses of nematode functional guilds to the different dominant plant species. Specifically, nematode communities reflected conditions of resource enrichment with increasing K. ericoides tree cover, whereas communities became structurally simplified and dominated by stress-tolerant nematode families with increasing P. nigra tree cover. Because nematodes regulate both bacterial- and fungal-dominated food webs in soils, these shifts could in turn alter multiple ecosystem processes belowground such as nutrient cycling. Incorporating species' functional traits into the assessment of habitat-change impacts on communities can greatly improve our understanding of species responses to environmental changes and their consequences in ecosystems.

Global effects of non-native tree species on multiple ecosystem services.

Non-native tree (NNT) species have been transported worldwide to create or enhance services that are fundamental for human well-being, such as timber provision, erosion control or ornamental value; yet NNTs can also produce undesired effects, such as fire proneness or pollen allergenicity. Despite the variety of effects that NNTs have on multiple ecosystem services, a global quantitative assessment of their costs and benefits is still lacking. Such information is critical for decision-making, management and sustainable exploitation of NNTs. We present here a global assessment of NNT effects on the three main categories of ecosystem services, including regulating (RES), provisioning (PES) and cultural services (CES), and on an ecosystem disservice (EDS), i.e. pollen allergenicity. By searching the scientific literature, country forestry reports, and social media, we compiled a global data set of 1683 case studies from over 125 NNT species, covering 44 countries, all continents but Antarctica, and seven biomes. Using different meta-analysis techniques, we found that, while NNTs increase most RES (e.g. climate regulation, soil erosion control, fertility and formation), they decrease PES (e.g. NNTs contribute less than native trees to global timber provision). Also, they have different effects on CES (e.g. increase aesthetic values but decrease scientific interest), and no effect on the EDS considered. NNT effects on each ecosystem (dis)service showed a strong context dependency, varying across NNT types, biomes and socio-economic conditions. For instance, some RES are increased more by NNTs able to fix atmospheric nitrogen, and when the ecosystem is located in low-latitude biomes; some CES are increased more by NNTs in less-wealthy countries or in countries with higher gross domestic products. The effects of NNTs on several ecosystem (dis)services exhibited some synergies (e.g. among soil fertility, soil formation and climate regulation or between aesthetic values and pollen allergenicity), but also trade-offs (e.g. between fire regulation and soil erosion control). Our analyses provide a quantitative understanding of the complex synergies, trade-offs and context dependencies involved for the effects of NNTs that is essential for attaining a sustained provision of ecosystem services.

Symmetric assembly and disassembly processes in an ecological network.

The processes whereby ecological networks emerge, persist and decay throughout ecosystem development are largely unknown. Here we study networks of plant and arbuscular mycorrhizal fungal (AMF) communities along a 120 000 year soil chronosequence, as they undergo assembly (progression) and then disassembly (retrogression). We found that network assembly and disassembly were symmetrical, self-reinforcing processes that together were capable of generating key attributes of network architecture. Plant and AMF species that had short indirect paths to others in the community (i.e. high centrality), rather than many direct interaction partners (i.e. high degree), were best able to attract new interaction partners and, in the case of AMF species, also to retain existing interactions with plants during retrogression. We then show using simulations that these non-random patterns of attachment and detachment promote nestedness of the network. These results have implications for predicting extinction sequences, identifying focal points for invasions and suggesting trajectories for restoration.

Belowground competition drives invasive plant impact on native species regardless of nitrogen availability.

Plant invasions and eutrophication are pervasive drivers of global change that cause biodiversity loss. Yet, how invasive plant impacts on native species, and the mechanisms underpinning these impacts, vary in relation to increasing nitrogen (N) availability remains unclear. Competition is often invoked as a likely mechanism, but the relative importance of the above and belowground components of this is poorly understood, particularly under differing levels of N availability. To help resolve these issues, we quantified the impact of a globally invasive grass species, Agrostis capillaris, on two co-occurring native New Zealand grasses, and vice versa. We explicitly separated above- and belowground interactions amongst these species experimentally and incorporated an N addition treatment. We found that competition with the invader had large negative impacts on native species growth (biomass decreased by half), resource capture (total N content decreased by up to 75%) and even nutrient stoichiometry (native species tissue C:N ratios increased). Surprisingly, these impacts were driven directly and indirectly by belowground competition, regardless of N availability. Higher root biomass likely enhanced the invasive grass's competitive superiority belowground, indicating that root traits may be useful tools for understanding invasive plant impacts. Our study shows that belowground competition can be more important in driving invasive plant impacts than aboveground competition in both low and high fertility ecosystems, including those experiencing N enrichment due to global change. This can help to improve predictions of how two key drivers of global change, plant species invasions and eutrophication, impact native species diversity.

Hierarchical neighbor effects on mycorrhizal community structure and function.

Theory predicts that neighboring communities can shape one another's composition and function, for example, through the exchange of member species. However, empirical tests of the directionality and strength of these effects are rare. We determined the effects of neighboring communities on one another through experimental manipulation of a plant-fungal model system. We first established distinct ectomycorrhizal fungal communities on Douglas-fir seedlings that were initially grown in three soil environments. We then transplanted seedlings and mycorrhizal communities in a fully factorial experiment designed to quantify the direction and strength of neighbor effects by focusing on changes in fungal community species composition and implications for seedling growth (a proxy for community function). We found that neighbor effects on the composition and function of adjacent communities follow a dominance hierarchy. Specifically, mycorrhizal communities established from soils collected in Douglas-fir plantations were both the least sensitive to neighbor effects, and exerted the strongest influence on their neighbors by driving convergence in neighbor community composition and increasing neighbor seedling vigor. These results demonstrate that asymmetric neighbor effects mediated by ecological history can determine both community composition and function.

Soil fertility and disturbance interact to drive contrasting responses of co-occurring native and nonnative species.

Some plant functional groups such as nonnative invasive and nitrogen (N)-fixing plants are widely thought to have consistent, coordinated differences in their functional traits relative to other groups such as native and non -N-fixing plants. Recent evidence suggests that these trait differences between groups can be context dependent, varying with environmental factors such as resource availability and disturbance. However, many previous comparisons among plant groups differing in invasion status have not standardized growth form between groups or have compared species that do not co-occur, which could result in invasion status per se being confounded with other factors. We determined growth and leaf functional trait responses of 20 co-occurring woody species, that is, five species within each of four functional groups (native N-fixers, native non -N-fixers, nonnative [invasive] N-fixers and nonnative [invasive] non-N-fixers), to factorial combinations of soil fertility and defoliation treatments in a mesocosm experiment to test each of two hypotheses. First, we hypothesized that nonnative invasive and N-fixing species will have functional traits associated with rapid resource acquisition whereas natives and non -N-fixing species will have traits linked to resource conservation. Second, we hypothesized that plant growth and leaf traits of nonnative and N-fixing species will be more strongly influenced by environmental factors (i.e., soil fertility and disturbance) than will natives and non-N-fixers. Plant growth, foliar nutrients, and leaf structural traits varied among plant functional groups in a manner consistent with our first hypothesis. Support for our second hypothesis was mixed; origin (native vs. nonnative) and soil fertility rarely interacted to determine plant growth or variation in leaf traits whereas interactions involving N-fixing ability and soil fertility were common. Further, there were no consistent interactive effects between plant groupings and disturbance. Our results demonstrate that variation in growth and functional traits among plant species were driven primarily by the relatively large responses of nonnative N-fixers to soil fertility, rather than by consistent differences between other plant functional groups. These findings highlight the importance of resource availability in determining trait or performance differences among plant functional groups, and provide insights into the assembly of plant functional traits in novel communities of co-occurring native and nonnative species.

Mycorrhizal co-invasion and novel interactions depend on neighborhood context.

Biological invasions are a rapidly increasing driver of global change, yet fundamental gaps remain in our understanding of the factors determining the success or extent of invasions. For example, although most woody plant species depend on belowground mutualists such as mycorrhizal fungi and nitrogen-fixing bacteria, the relative importance of these mutualisms in conferring invasion success is unresolved. Here, we describe how neighborhood context (identity of nearby tree species) affects the formation of belowground ectomycorrhizal partnerships between fungi and seedlings of a widespread invasive tree species, Pseudotsuga menziesii (Douglas-fir), in New Zealand. We found that the formation of mycorrhizal partnerships, the composition of the fungal species involved in these partnerships, and the origin of the fungi (co-invading or native to New Zealand) all depend on neighborhood context. Our data suggest that nearby ectomycorrhizal host trees act as both a reservoir of fungal inoculum and a carbon source for late-successional and native fungi. By facilitating mycorrhization of P. menziesii seedlings, adult trees may alleviate mycorrhizal limitation at the P. menziesii invasion front. These results highlight the importance of studying biological invasions across multiple ecological settings to understand establishment success and invasion speed.

A global meta-analysis of the relative extent of intraspecific trait variation in plant communities.

Recent studies have shown that accounting for intraspecific trait variation (ITV) may better address major questions in community ecology. However, a general picture of the relative extent of ITV compared to interspecific trait variation in plant communities is still missing. Here, we conducted a meta-analysis of the relative extent of ITV within and among plant communities worldwide, using a data set encompassing 629 communities (plots) and 36 functional traits. Overall, ITV accounted for 25% of the total trait variation within communities and 32% of the total trait variation among communities on average. The relative extent of ITV tended to be greater for whole-plant (e.g. plant height) vs. organ-level traits and for leaf chemical (e.g. leaf N and P concentration) vs. leaf morphological (e.g. leaf area and thickness) traits. The relative amount of ITV decreased with increasing species richness and spatial extent, but did not vary with plant growth form or climate. These results highlight global patterns in the relative importance of ITV in plant communities, providing practical guidelines for when researchers should include ITV in trait-based community and ecosystem studies.

Host identity is a dominant driver of mycorrhizal fungal community composition during ecosystem development.

Little is known about the response of arbuscular mycorrhizal fungal communities to ecosystem development. We use a long-term soil chronosequence that includes ecosystem progression and retrogression to quantify the importance of host plant identity as a factor driving fungal community composition during ecosystem development. We identified arbuscular mycorrhizal fungi and plant species from 50 individual roots from each of 10 sites spanning 5-120 000 yr of ecosystem age using terminal restriction fragment length polymorphism (T-RFLP), Sanger sequencing and pyrosequencing. Arbuscular mycorrhizal fungal communities were highly structured by ecosystem age. There was strong niche differentiation, with different groups of operational taxonomic units (OTUs) being characteristic of early succession, ecosystem progression and ecosystem retrogression. Fungal alpha diversity decreased with ecosystem age, whereas beta diversity was high at early stages and lower in subsequent stages. A total of 39% of the variance in fungal communities was explained by host plant and site age, 29% of which was attributed to host and the interaction between host and site (24% and 5%, respectively). The strong response of arbuscular mycorrhizal fungi to ecosystem development appears to be largely driven by plant host identity, supporting the concept that plant and fungal communities are tightly coupled rather than independently responding to habitat.

Belowground legacies of Pinus contorta invasion and removal result in multiple mechanisms of invasional meltdown.

Plant invasions can change soil biota and nutrients in ways that drive subsequent plant communities, particularly when co-invading with belowground mutualists such as ectomycorrhizal fungi. These effects can persist following removal of the invasive plant and, combined with effects of removal per se, influence subsequent plant communities and ecosystem functioning. We used field observations and a soil bioassay with multiple plant species to determine the belowground effects and post-removal legacy caused by invasion of the non-native tree Pinus contorta into a native plant community. Pinus facilitated ectomycorrhizal infection of the co-occurring invasive tree, Pseudotsuga menziesii, but not conspecific Pinus (which always had ectomycorrhizas) nor the native pioneer Kunzea ericoides (which never had ectomycorrhizas). Pinus also caused a major shift in soil nutrient cycling as indicated by increased bacterial dominance, NO3-N (17-fold increase) and available phosphorus (3.2-fold increase) in soils, which in turn promoted increased growth of graminoids. These results parallel field observations, where Pinus removal is associated with invasion by non-native grasses and herbs, and suggest that legacies of Pinus on soil nutrient cycling thus indirectly promote invasion of other non-native plant species. Our findings demonstrate that multi-trophic belowground legacies are an important but hitherto largely unconsidered factor in plant community reassembly following invasive plant removal.

Intraspecific relationships among wood density, leaf structural traits and environment in four co-occurring species of Nothofagus in New Zealand.

Plant functional traits capture important variation in plant strategy and function. Recent literature has revealed that within-species variation in traits is greater than previously supposed. However, we still have a poor understanding of how intraspecific variation is coordinated among different traits, and how it is driven by environment. We quantified intraspecific variation in wood density and five leaf traits underpinning the leaf economics spectrum (leaf dry matter content, leaf mass per unit area, size, thickness and density) within and among four widespread Nothofagus tree species in southern New Zealand. We tested whether intraspecific relationships between wood density and leaf traits followed widely reported interspecific relationships, and whether variation in these traits was coordinated through shared responses to environmental factors. Sample sites varied widely in environmental variables, including soil fertility (25-900 mg kg(-1) total P), precipitation (668-4875 mm yr(-1)), temperature (5.2-12.4 °C mean annual temperature) and latitude (41-46 °S). Leaf traits were strongly correlated with one another within species, but not with wood density. There was some evidence for a positive relationship between wood density and leaf tissue density and dry matter content, but no evidence that leaf mass or leaf size were correlated with wood density; this highlights that leaf mass per unit area cannot be used as a surrogate for component leaf traits such as tissue density. Trait variation was predicted by environmental factors, but not consistently among different traits; e.g., only leaf thickness and leaf density responded to the same environmental cues as wood density. We conclude that although intraspecific variation in wood density and leaf traits is strongly driven by environmental factors, these responses are not strongly coordinated among functional traits even across co-occurring, closely-related plant species.

Declining foliar and litter δ¹5N diverge from soil, epiphyte and input δ¹5N along a 120,000 yr temperate rainforest chronosequence.

Patterns in the natural abundance of nitrogen (N) isotopes (¹5N and ¹4N) can help in the understanding of ecosystem processes along environmental gradients, because some processes fractionate against the heavier isotope. We measured δ¹5N in many components of the Franz Josef soil chronosequence in New Zealand to see how each component varied along the sequence and within sites, and to see what this variation can tell us about how ecosystem processes such as N losses change with soil age. We analyzed δ¹5N in foliage from 18 woody species, abscised leaves from seven woody species, three soil horizons, bryophytes, lichens, bulk deposition, and nodules from the N-fixing tree Coriaria arborea (Coriariaceae). Foliar δ¹5N varied significantly across plant species. Foliage and bulk litter became ¹5N-depleted as soil age increased. Soil N from organic and mineral horizons was significantly more ¹5N-enriched than bulk litter N at each site. Increasing precipitation also decreased foliar and soil δ¹5N. Comparing input and whole ecosystem δ¹5N revealed limited evidence for net fractionation during N losses. These trends are consistent with some combination of increasing fractionation during plant N uptake, mycorrhizal transfer, within-plant processing, and soil decomposition as soils age.

Co-invasion by Pinus and its mycorrhizal fungi.

SUMMARY: *The absence of co-evolved mutualists of plants invading a novel habitat is the logical corollary of the more widely recognized 'enemy escape'. To avoid or overcome the loss of mutualists, plants may co-invade with nonnative mutualists, form novel associations with native mutualists or form associations with native cosmopolitan mutualists, which are native but not novel to the invading plant. *We tested these hypotheses by contrasting the ectomycorrhizal fungal communities associated with invasive Pinus contorta in New Zealand with co-occurring endemic Nothofagus solandri var. cliffortioides. *Fungal communities on Pinus were species poor (14 ectomycorrhizal species) and dominated by nonnative (93%) and cosmopolitan fungi (7%). Nothofagus had a species-rich (98 species) fungal community dominated by native Cortinarius and two cosmopolitan fungi. *These results support co-invasion by mutualists rather than novel associations as an important mechanism by which plants avoid or overcome the loss of mutualists, consistent with invasional meltdown.

Photosynthesis and reflectance indices for rainforest species in ecosystems undergoing progression and retrogression along a soil fertility chronosequence in New Zealand.

Measurements of photosynthesis at saturating irradiance and CO2 partial pressure, Amax, "adjusted" normalised difference vegetation index, RaNDVI, and photochemical reflectance index, RPRI, were made on trees sampled along a soil chronosequence to investigate the relationship between carbon uptake and ecosystem development in relation to nutrient availability. Measurements were made on the three most dominant species at six sites along the sequence in South Westland, New Zealand with soil age ranging from < 6 to 120,000 years resulting from the retreat of the Franz Josef glacier. The decrease in soil phosphorus availability with increasing soil age and high soil nitrogen availability at the two youngest sites, due to the presence of a nitrogen-fixing species, provided marked differences in nutrient availability. Mean Amax was high at the two youngest sites, then decreased markedly with increasing site age. Analysis of the data for individual species within sites revealed separation of groups of species in the response of Amax to Nm and Pm, suggesting complex interactions between the two nutrients. There were strong linear relationships for leaf-level RaNDVI and RPRI with Amax, at high irradiance, showing that measurements of reflectance indices can be used to estimate Amax for foliage with a range in morphology and nutrient concentrations. Notwithstanding the change in species composition from angiosperms to conifers with increasing site age, the presence of nitrogen-fixing species, the variability in foliage morphology from flat leaves to imbricate scales and a wide range in foliar nitrogen and phosphorus concentrations, there were strong positive linear relationships between site average Amax and foliage nitrogen, Nm, and phosphorus, Pm, concentrations on a foliage mass basis. The results provide insights to interpret the regulation of photosynthesis across natural ecosystems with marked gradients in nitrogen and phosphorus availability.