Feedback from plant species change amplifies CO2 enhancement of grassland productivity

Elevated CO2 interacts with nutrient inputs to restructure plant communities in phosphorus-limited grasslands

Globally pervasive increases in atmospheric CO2 and nitrogen (N) deposition could have substantial effects on plant communities, either directly or mediated by their interactions with soil nutrient limitation. While the direct consequences of N enrichment on plant communities are well documented, potential interactions with rising CO2 and globally widespread phosphorus (P) limitation remain poorly understood. We investigated the consequences of simultaneous elevated CO2 (eCO2 ) and N and P additions on grassland biodiversity, community and functional composition in P-limited grasslands. We exposed soil-turf monoliths from limestone and acidic grasslands that have received >25 years of N additions (3.5 and 14 g m-2  year-1 ) and 11 (limestone) or 25 (acidic) years of P additions (3.5 g m-2  year-1 ) to eCO2 (600 ppm) for 3 years. Across both grasslands, eCO2 , N and P additions significantly changed community composition. Limestone communities were more responsive to eCO2 and saw significant functional shifts resulting from eCO2 -nutrient interactions. Here, legume cover tripled in response to combined eCO2 and P additions, and combined eCO2 and N treatments shifted functional dominance from grasses to sedges. We suggest that eCO2 may disproportionately benefit P acquisition by sedges by subsidising the carbon cost of locally intense root exudation at the expense of co-occurring grasses. In contrast, the functional composition of the acidic grassland was insensitive to eCO2 and its interactions with nutrient additions. Greater diversity of P-acquisition strategies in the limestone grassland, combined with a more functionally even and diverse community, may contribute to the stronger responses compared to the acidic grassland. Our work suggests we may see large changes in the composition and biodiversity of P-limited grasslands in response to eCO2 and its interactions with nutrient loading, particularly where these contain a high diversity of P-acquisition strategies or developmentally young soils with sufficient bioavailable mineral P.

Low precipitation due to climate change consistently reduces multifunctionality of urban grasslands in mesocosms

Urban grasslands are crucial for biodiversity and ecosystem services in cities, while little is known about their multifunctionality under climate change. Thus, we investigated the effects of simulated climate change, i.e., increased [CO2] and temperature, and reduced precipitation, on individual functions and overall multifunctionality in mesocosm grasslands sown with forbs and grasses in four different proportions aiming at mimicking road verge grassland patches. Climate change scenarios RCP2.6 (control) and RCP8.5 (worst-case) were simulated in walk-in climate chambers of an ecotron facility, and watering was manipulated for normal vs. reduced precipitation. We measured eight indicator variables of ecosystem functions based on below- and aboveground characteristics. The young grassland communities responded to higher [CO2] and warmer conditions with increased vegetation cover, height, flower production, and soil respiration. Lower precipitation affected carbon cycling in the ecosystem by reducing biomass production and soil respiration. In turn, the water regulation capacity of the grasslands depended on precipitation interacting with climate change scenario, given the enhanced water efficiency resulting from increased [CO2] under RCP8.5. Multifunctionality was negatively affected by reduced precipitation, especially under RCP2.6. Trade-offs arose among single functions that performed best in either grass- or forb-dominated grasslands. Grasslands with an even ratio of plant functional types coped better with climate change and thus are good options for increasing the benefits of urban green infrastructure. Overall, the study provides experimental evidence of the effects of climate change on the functionality of urban ecosystems. Designing the composition of urban grasslands based on ecological theory may increase their resilience to global change.

The stimulatory effect of elevated CO2 on soil respiration is unaffected by N addition

Global atmospheric CO₂ keeps rising and brings about significant effects on ecosystem carbon (C) cycling by altering C processes in soils. Soil C responses to elevated CO₂ are highly uncertain, and how elevated CO₂ interacts with other factors, such as nitrogen (N) availability, to influence soil C flux comprises an important source of this uncertainty, especially for those under-studied ecosystems. By conducting a manipulated CO₂ concentration and N availability experiment on typical alpine grassland (4600 m asl), we combined the five-year in-situ measurement of soil respiration (SR) with an incubation experiment of microbial metabolic efficiency in the lab to explore the response of SR to elevated CO₂ and N availability. The results showed that elevated CO₂ at ambient N conditions and enriched N equally stimulated SR during the experimental period, whereas N supply had no significant effect. Elevated CO₂ enhanced soil dissolved organic C and enzyme activity, while had marginal effects on microbial biomass and C use efficiency (CUE). Strengthened microbial activity dominated SR stimulation under elevated CO₂. Enriched N boosted enzyme activity and microbial CUE. N availability played divergent roles in mediating SR. The negliable regulation of N supply on elevated CO₂ effects on SR was the offset consequences of the negative impacts of enhanced CUE and the positive contribution of heightened enzyme activity. Our findings suggest that rising CO₂ would accelerate soil C cycling of the alpine grassland under various N regimes by stimulating microbial activity instead of lowering microbial metabolic efficiency. Such results are crucial for understanding the role of alpine ecosystems in the global C cycle.

Multiple constraints cause positive and negative feedbacks limiting grassland soil CO2 efflux under CO2 enrichment

Terrestrial ecosystems are increasingly enriched with resources such as atmospheric CO₂ that limit ecosystem processes. The consequences for ecosystem carbon cycling depend on the feedbacks from other limiting resources and plant community change, which remain poorly understood for soil CO₂ efflux, J_CO2, a primary carbon flux from the biosphere to the atmosphere. We applied a unique CO₂ enrichment gradient (250 to 500 µL L^-1) for eight years to grassland plant communities on soils from different landscape positions. We identified the trajectory of J_CO2 responses and feedbacks from other resources, plant diversity [effective species richness, exp(H)], and community change (plant species turnover). We found linear increases in J_CO2 on an alluvial sandy loam and a lowland clay soil, and an asymptotic increase on an upland silty clay soil. Structural equation modeling identified CO₂ as the dominant limitation on J_CO2 on the clay soil. In contrast with theory predicting limitation from a single limiting factor, the linear J_CO2 response on the sandy loam was reinforced by positive feedbacks from aboveground net primary productivity and exp(H), while the asymptotic J_CO2 response on the silty clay arose from a net negative feedback among exp(H), species turnover, and soil water potential. These findings support a multiple resource limitation view of the effects of global change drivers on grassland ecosystem carbon cycling and highlight a crucial role for positive or negative feedbacks between limiting resources and plant community structure. Incorporating these feedbacks will improve models of terrestrial carbon sequestration and ecosystem services.

Global change biology: A primer

Because of human action, the Earth has entered an era where profound changes in the global environment are creating novel conditions that will be discernable far into the future. One consequence may be a large reduction of the Earth's biodiversity, potentially representing a sixth mass extinction. With effective stewardship, the global change drivers that threaten the Earth's biota could be alleviated, but this requires clear understanding of the drivers, their interactions, and how they impact ecological communities. This review identifies 10 anthropogenic global change drivers and discusses how six of the drivers (atmospheric CO₂ enrichment, climate change, land transformation, species exploitation, exotic species invasions, eutrophication) impact Earth's biodiversity. Driver impacts on a particular species could be positive or negative. In either case, they initiate secondary responses that cascade along ecological lines of connection and in doing so magnify the initial impact. The unique nature of the threat to the Earth's biodiversity is not simply due to the magnitude of each driver, but due to the speed of change, the novelty of the drivers, and their interactions. Emphasizing one driver, notably climate change, is problematic because the other global change drivers also degrade biodiversity and together threaten the stability of the biosphere. As the main academic journal addressing global change effects on living systems, GCB is well positioned to provide leadership in solving the global change challenge. If humanity cannot meet the challenge, then GCB is positioned to serve as a leading chronicle of the sixth mass extinction to occur on planet Earth.

CO2 enrichment and soil type additively regulate grassland productivity

Atmospheric CO₂ enrichment usually increases the aboveground net primary productivity (ANPP) of grassland vegetation, but the magnitude of the ANPP-CO₂ response differs among ecosystems. Soil properties affect ANPP via multiple mechanisms and vary over topographic to geographic gradients, but have received little attention as potential modifiers of the ANPP-CO₂ response. We assessed the effects of three soil types, sandy loam, silty clay and clay, on the ANPP response of perennial C₃ /C₄ grassland communities to a subambient to elevated CO₂ gradient over 10 yr in Texas, USA. We predicted an interactive, rather than additive, effect of CO₂ and soil type on ANPP. Contrary to prediction, CO₂ and soil additively influenced grassland ANPP. Increasing CO₂ by 250 μl l^-1 increased ANPP by 170 g m^-2 across soil types. Increased clay content from 10% to 50% among soils reduced ANPP by 50 g m^-2 . CO₂ enrichment increased ANPP via a predominant direct effect, accompanied by a smaller indirect effect mediated by a successional shift to increased dominance of the C₄ tallgrass Sorghastrum nutans. Our results indicate a large, positive influence of CO₂ enrichment on grassland productivity that resulted from the direct physiological benefits of CO₂ augmented by species succession, and was expressed similarly across soils of differing physical properties.

Flowering in grassland predicted by CO2 and resource effects on species aboveground biomass

Continuing enrichment of atmospheric CO₂ may change plant community composition, in part by altering the availability of other limiting resources including soil water, nutrients, or light. The combined effects of CO₂ enrichment and altered resource availability on species flowering remain poorly understood. We quantified flowering culm and ramet production and biomass allocation to flowering culms/ramets for 10 years in C₄ -dominated grassland communities on contrasting soils along a CO₂ concentration gradient spanning pre-industrial to expected mid-21st century levels (250-500 μl/L). CO₂ enrichment explained up to 77% of the variation in flowering culm count across soils for three of the five species, and was correlated with flowering culm count on at least one soil for four of five species. In contrast, allocation to flowering culms was only weakly correlated with CO₂ enrichment for two species. Flowering culm counts were strongly correlated with species aboveground biomass (AGB; R²  = .34-.74), a measure of species abundance. CO₂ enrichment also increased soil moisture and decreased light levels within the canopy but did not affect soil inorganic nitrogen availability. Structural equation models fit across the soils suggested species-specific controls on flowering in two general forms: (1) CO₂ effects on flowering culm count mediated by canopy light level and relative species AGB (species AGB/total AGB) or by soil moisture effects on flowering culm count; (2) effects of canopy light level or soil inorganic nitrogen on flowering and/or relative species AGB, but with no significant CO₂ effect. Understanding the heterogeneity in species responses to CO₂ enrichment in plant communities across soils in edaphically variable landscapes is critical to predict CO₂ effects on flowering and other plant fitness components, and species potential to adapt to future environmental changes.

A portrait of the C4 photosynthetic family on the 50th anniversary of its discovery: species number, evolutionary lineages, and Hall of Fame

Fifty years ago, the C₄ photosynthetic pathway was first characterized. In the subsequent five decades, much has been learned about C₄ plants, such that it is now possible to place nearly all C₄ species into their respective evolutionary lineages. Sixty-one independent lineages of C₄ photosynthesis are identified, with additional, ancillary C₄ origins possible in 12 of these principal lineages. The lineages produced ~8100 C₄ species (5044 grasses, 1322 sedges, and 1777 eudicots). Using midpoints of stem and crown node dates in their respective phylogenies, the oldest and most speciose C₄ lineage is the grass lineage Chloridoideae, estimated to be near 30 million years old. Most C₄ lineages are estimated to be younger than 15 million years. Older C₄ lineages tend to be more speciose, while those younger than 7 million years have <43 species each. To further highlight C₄ photosynthesis for a 50th anniversary snapshot, a Hall of Fame comprised of the 40 most significant C₄ species is presented. Over the next 50 years, preservation of the Earth's C₄ diversity is a concern, largely because of habitat loss due to elevated CO₂ effects, invasive species, and expanded agricultural activities. Ironically, some members of the C₄ Hall of Fame are leading threats to the natural C₄ flora due to their association with human activities on landscapes where most C₄ plants occur.

A portrait of the C4 photosynthetic family on the 50th anniversary of its discovery: species number, evolutionary lineages, and Hall of Fame

Fifty years ago, the C4 photosynthetic pathway was first characterized. In the subsequent five decades, much has been learned about C4 plants, such that it is now possible to place nearly all C4 species into their respective evolutionary lineages. Sixty-one independent lineages of C4 photosynthesis are identified, with additional, ancillary C4 origins possible in 12 of these principal lineages. The lineages produced ~8100 C4 species (5044 grasses, 1322 sedges, and 1777 eudicots). Using midpoints of stem and crown node dates in their respective phylogenies, the oldest and most speciose C4 lineage is the grass lineage Chloridoideae, estimated to be near 30 million years old. Most C4 lineages are estimated to be younger than 15 million years. Older C4 lineages tend to be more speciose, while those younger than 7 million years have <43 species each. To further highlight C4 photosynthesis for a 50th anniversary snapshot, a Hall of Fame comprised of the 40 most significant C4 species is presented. Over the next 50 years, preservation of the Earth's C4 diversity is a concern, largely because of habitat loss due to elevated CO2 effects, invasive species, and expanded agricultural activities. Ironically, some members of the C4 Hall of Fame are leading threats to the natural C4 flora due to their association with human activities on landscapes where most C4 plants occur.

Impacts of warming and elevated CO2 on a semi-arid grassland are non-additive, shift with precipitation, and reverse over time

It is unclear how elevated CO2 (eCO2 ) and the corresponding shifts in temperature and precipitation will interact to impact ecosystems over time. During a 7-year experiment in a semi-arid grassland, the response of plant biomass to eCO2 and warming was largely regulated by interannual precipitation, while the response of plant community composition was more sensitive to experiment duration. The combined effects of eCO2 and warming on aboveground plant biomass were less positive in 'wet' growing seasons, but total plant biomass was consistently stimulated by ~ 25% due to unique, supra-additive responses of roots. Independent of precipitation, the combined effects of eCO2 and warming on C3 graminoids became increasingly positive and supra-additive over time, reversing an initial shift toward C4 grasses. Soil resources also responded dynamically and non-additively to eCO2 and warming, shaping the plant responses. Our results suggest grasslands are poised for drastic changes in function and highlight the need for long-term, factorial experiments.

Soil carbon sequestration is a climate stabilization wedge: comments on Sommer and Bossio (2014)

Sommer and Bossio (2014) model the potential soil organic carbon (SOC) sequestration in agricultural soils (croplands and grasslands) during the next 87 years, concluding that this process cannot be considered as a climate stabilization wedge. We argue, however, that the amounts of SOC potentially sequestered in both scenarios (pessimistic and optimistic) fulfil the requirements for being considered as wedge because in both cases at least 25 GtC would be sequestered during the next 50 years. We consider that it is precisely in the near future, and meanwhile other solutions are developed, when this stabilization effort is most urgent even if after some decades the sequestration rate is significantly reduced. Indirect effects of SOC sequestration on mitigation could reinforce the potential of this solution. We conclude that the sequestration of organic carbon in agricultural soils as a climate change mitigation tool still deserves important attention for scientists, managers and policy makers.

Dominant plant taxa predict plant productivity responses to CO2 enrichment across precipitation and soil gradients

The Earth's atmosphere will continue to be enriched with carbon dioxide (CO2) over the coming century. Carbon dioxide enrichment often reduces leaf transpiration, which in water-limited ecosystems may increase soil water content, change species abundances and increase the productivity of plant communities. The effect of increased soil water on community productivity and community change may be greater in ecosystems with lower precipitation, or on coarser-textured soils, but responses are likely absent in deserts. We tested correlations among yearly increases in soil water content, community change and community plant productivity responses to CO2 enrichment in experiments in a mesic grassland with fine- to coarse-textured soils, a semi-arid grassland and a xeric shrubland. We found no correlation between CO2-caused changes in soil water content and changes in biomass of dominant plant taxa or total community aboveground biomass in either grassland type or on any soil in the mesic grassland (P > 0.60). Instead, increases in dominant taxa biomass explained up to 85 % of the increases in total community biomass under CO2 enrichment. The effect of community change on community productivity was stronger in the semi-arid grassland than in the mesic grassland, where community biomass change on one soil was not correlated with the change in either the soil water content or the dominant taxa. No sustained increases in soil water content or community productivity and no change in dominant plant taxa occurred in the xeric shrubland. Thus, community change was a crucial driver of community productivity responses to CO2 enrichment in the grasslands, but effects of soil water change on productivity were not evident in yearly responses to CO2 enrichment. Future research is necessary to isolate and clarify the mechanisms controlling the temporal and spatial variations in the linkages among soil water, community change and plant productivity responses to CO2 enrichment.

Plant community change mediates the response of foliar δ(15)N to CO 2 enrichment in mesic grasslands

Rising atmospheric CO2 concentration may change the isotopic signature of plant N by altering plant and microbial processes involved in the N cycle. CO2 may increase leaf δ(15)N by increasing plant community productivity, C input to soil, and, ultimately, microbial mineralization of old, (15)N-enriched organic matter. We predicted that CO2 would increase aboveground productivity (ANPP; g biomass m(-2)) and foliar δ(15)N values of two grassland communities in Texas, USA: (1) a pasture dominated by a C4 exotic grass, and (2) assemblages of tallgrass prairie species, the latter grown on clay, sandy loam, and silty clay soils. Grasslands were exposed in separate experiments to a pre-industrial to elevated CO2 gradient for 4 years. CO2 stimulated ANPP of pasture and of prairie assemblages on each of the three soils, but increased leaf δ(15)N only for prairie plants on a silty clay. δ(15)N increased linearly as mineral-associated soil C declined on the silty clay. Mineral-associated C declined as ANPP increased. Structural equation modeling indicted that CO2 increased ANPP partly by favoring a tallgrass (Sorghastrum nutans) over a mid-grass species (Bouteloua curtipendula). CO2 may have increased foliar δ(15)N on the silty clay by reducing fractionation during N uptake and assimilation. However, we interpret the soil-specific, δ(15)N-CO2 response as resulting from increased ANPP that stimulated mineralization from recalcitrant organic matter. By contrast, CO2 favored a forb species (Solanum dimidiatum) with higher δ(15)N than the dominant grass (Bothriochloa ischaemum) in pasture. CO2 enrichment changed grassland δ(15)N by shifting species relative abundances.

Fungal Community Responses to Past and Future Atmospheric CO2 Differ by Soil Type

Soils sequester and release substantial atmospheric carbon, but the contribution of fungal communities to soil carbon balance under rising CO2 is not well understood. Soil properties likely mediate these fungal responses but are rarely explored in CO2 experiments. We studied soil fungal communities in a grassland ecosystem exposed to a preindustrial-to-future CO2 gradient (250 to 500 ppm) in a black clay soil and a sandy loam soil. Sanger sequencing and pyrosequencing of the rRNA gene cluster revealed that fungal community composition and its response to CO2 differed significantly between soils. Fungal species richness and relative abundance of Chytridiomycota (chytrids) increased linearly with CO2 in the black clay (P < 0.04, R(2) > 0.7), whereas the relative abundance of Glomeromycota (arbuscular mycorrhizal fungi) increased linearly with elevated CO2 in the sandy loam (P = 0.02, R(2) = 0.63). Across both soils, decomposition rate was positively correlated with chytrid relative abundance (r = 0.57) and, in the black clay soil, fungal species richness. Decomposition rate was more strongly correlated with microbial biomass (r = 0.88) than with fungal variables. Increased labile carbon availability with elevated CO2 may explain the greater fungal species richness and Chytridiomycota abundance in the black clay soil, whereas increased phosphorus limitation may explain the increase in Glomeromycota at elevated CO2 in the sandy loam. Our results demonstrate that soil type plays a key role in soil fungal responses to rising atmospheric CO2.

Impacts of climate change drivers on C4 grassland productivity: scaling driver effects through the plant community

Climate change drivers affect plant community productivity via three pathways: (i) direct effects of drivers on plants; (ii) the response of species abundances to drivers (community response); and (iii) the feedback effect of community change on productivity (community effect). The contribution of each pathway to driver-productivity relationships depends on functional traits of dominant species. We used data from three experiments in Texas, USA, to assess the role of community dynamics in the aboveground net primary productivity (ANPP) response of C4 grasslands to two climate drivers applied singly: atmospheric CO2 enrichment and augmented summer precipitation. The ANPP-driver response differed among experiments because community responses and effects differed. ANPP increased by 80-120g m(-2) per 100 μl l(-1) rise in CO2 in separate experiments with pasture and tallgrass prairie assemblages. Augmenting ambient precipitation by 128mm during one summer month each year increased ANPP more in native than in exotic communities in a third experiment. The community effect accounted for 21-38% of the ANPP CO2 response in the prairie experiment but little of the response in the pasture experiment. The community response to CO2 was linked to species traits associated with greater soil water from reduced transpiration (e.g. greater height). Community effects on the ANPP CO2 response and the greater ANPP response of native than exotic communities to augmented precipitation depended on species differences in transpiration efficiency. These results indicate that feedbacks from community change influenced ANPP-driver responses. However, the species traits that regulated community effects on ANPP differed from the traits that determined how communities responded to drivers.

Selective grazing modifies previously anticipated responses of plant community composition to elevated CO(2) in a temperate grassland

Our limited understanding of terrestrial ecosystem responses to elevated CO2 is a major constraint on predicting the impacts of climate change. A change in botanical composition has been identified as a key factor in the CO2 response with profound implications for ecosystem services such as plant production and soil carbon storage. In temperate grasslands, there is a strong consensus that elevated CO2 will result in a greater physiological stimulus to growth in legumes and to a lesser extent forbs, compared with C3 grasses, and the presumption this will lead in turn to a greater proportion of these functional groups in the plant community. However, this view is based on data mainly collected in experiments of three or less years in duration and not in experiments where defoliation has been by grazing animals. Grazing is, however, the most common management of grasslands and known in itself to influence botanical composition. In a long-term Free Air Carbon Dioxide Enrichment (FACE) experiment in a temperate grassland managed with grazing animals (sheep), we found the response to elevated CO2 in plant community composition in the first 5 years was consistent with the expectation of increased proportions of legumes and forbs. However, in the longer term, these differences diminished so that the proportions of grasses, legumes and forbs were the same under both ambient and elevated CO2 . Analysis of vegetation before and after each grazing event showed there was a sustained disproportionately greater removal ('apparent selection') of legumes and forbs by the grazing animals. This bias in removal was greater under elevated CO2 than ambient CO2 . This is consistent with sustained faster growth rates of legumes and forbs under elevated CO2 being countered by selective defoliation, and so leading to little difference in community composition.

The response of aboveground net primary productivity of desert vegetation to rainfall pulse in the temperate desert region of northwest China

Rainfall events can be characterized as "pulses", which are discrete and variable episodes that can significantly influence the structure and function of desert ecosystems, including shifts in aboveground net primary productivity (ANPP). To determine the threshold and hierarchical response of rainfall event size on the Normalized Difference Vegetation Index (NDVI, a proxy for ANPP) and the difference across a desert area in northwestern China with two habitats - dune and desert - we selected 17 independent summer rainfall events from 2005 to 2012, and obtained a corresponding NDVI dataset extracted from MODIS images. Based on the threshold-delay model and statistical analysis, the results showed that the response of NDVI to rainfall pulses began at about a 5 mm event size. Furthermore, when the rainfall event size was more than 30 mm, NDVI rapidly increased 3- to 6-fold compared with the response to events of less than 30 mm, suggesting that 30 mm was the threshold for a large NDVI response. These results revealed the importance of the 5 mm and 30 mm rainfall events for plant survival and growth in desert regions. There was an 8- to 16-day lag time between the rainfall event and the NDVI response, and the response duration varied with rainfall event size, reaching a maximum of 32 days. Due to differences in soil physical and mineralogical properties, and to biodiversity structure and the root systems' abilities to exploit moisture, dune and desert areas differed in precipitation responses: dune habitats were characterized by a single, late summer productivity peak; in contrast, deserts showed a multi-peak pattern throughout the growing season.

Net primary productivity and rain-use efficiency as affected by warming, altered precipitation, and clipping in a mixed-grass prairie

Grassland productivity in response to climate change and land use is a global concern. In order to explore the effects of climate change and land use on net primary productivity (NPP), NPP partitioning [fBNPP , defined as the fraction of belowground NPP (BNPP) to NPP], and rain-use efficiency (RUE) of NPP, we conducted a field experiment with warming (+3 °C), altered precipitation (double and half), and annual clipping in a mixed-grass prairie in Oklahoma, USA since July, 2009. Across the years, warming significantly increased BNPP, fBNPP , and RUEBNPP by an average of 11.6%, 2.8%, and 6.6%, respectively. This indicates that BNPP was more sensitive to warming than aboveground NPP (ANPP) since warming did not change ANPP and RUEANPP much. Double precipitation stimulated ANPP, BNPP, and NPP but suppressed RUEANPP , RUEBNPP , and RUENPP while half precipitation decreased ANPP, BNPP, and NPP but increased RUEANPP , RUEBNPP , and RUENPP . Clipping interacted with altered precipitation in impacting RUEANPP , RUEBNPP , and RUENPP , suggesting land use could confound the effects of precipitation changes on ecosystem processes. Soil moisture was found to be a main factor in regulating variation in ANPP, BNPP, and NPP while soil temperature was the dominant factor influencing fBNPP . These findings suggest that BNPP is critical point to future research. Additionally, results from single-factor manipulative experiments should be treated with caution due to the non-additive interactive effects of warming with altered precipitation and land use (clipping).