Needle bacterial community structure across the species range of limber pine.

Bacteria on and inside leaves can influence forest tree health and resilience. The distribution and limits of a tree species' range can be influenced by various factors, with biological interactions among the most significant. We investigated the processes shaping the bacterial needle community across the species distribution of limber pine, a widespread Western conifer inhabiting a range of extreme habitats. We tested four hypotheses: (i) Needle community structure varies across sites, with site-specific factors more important to microbial assembly than host species selection; (ii) dispersal limitation structures foliar communities across the range of limber pine; (iii) the relative significance of dispersal and selection differs across sites in the tree species range; and (iv) needle age structures bacterial communities. We characterized needle communities from the needle surface and tissue of limber pine and co-occurring conifers across 16 sites in the limber pine distribution. Our findings confirmed that site characteristics shape the assembly of bacterial communities across the host species range and showed that these patterns are not driven by dispersal limitation. Furthermore, the strength of selection by the host varied by site, possibly due to differences in available microbes. Our study, by focusing on trees in their natural setting, reveals real needle bacterial dynamics in forests, which is key to understanding the balance between stochastic and deterministic processes in shaping forest tree-microbe interactions. Such understanding will be necessary to predict or manipulate these interactions to support forest ecosystem productivity or assist plant migration and adaptation in the face of global change.

Allometric relationships and trade-offs in 11 common Mediterranean-climate grasses.

Biomass allocation in plants is the foundation for understanding dynamics in ecosystem carbon balance, species competition, and plant-environment interactions. However, existing work on plant allometry has mainly focused on trees, with fewer studies having developed allometric equations for grasses. Grasses with different life histories can vary in their carbon investment by prioritizing the growth of specific organs to survive, outcompete co-occurring plants, and ensure population persistence. Further, because grasses are important fuels for wildfire, the lack of grass allocation data adds uncertainty to process-based models that relate plant physiology to wildfire dynamics. To fill this gap, we conducted a greenhouse experiment with 11 common California grasses varying in photosynthetic pathway and growth form. We measured plant sizes and harvested above- and belowground biomass throughout the life cycle of annual species, while for the establishment stage of perennial grasses to quantify allometric relationships for leaf, stem, and root biomass, as well as plant height and canopy area. We used basal diameter as a reference measure of plant size. Overall, basal diameter is the best predictor for leaf and stem biomass, height, and canopy area. Including height as another predictor can improve model accuracy in predicting leaf and stem biomass and canopy area. Fine root biomass is a function of leaf biomass alone. Species vary in their allometric relationships, with most variation occurring for plant height, canopy area, and stem biomass. We further explored potential trade-offs in biomass allocation across species between leaf and fine root, leaf and stem, and allocation to reproduction. Consistent with our expectation, we found that fast-growing plants allocated a greater fraction to reproduction. Additionally, plant height and specific leaf area negatively influenced the leaf-to-stem ratio. However, contrary to our hypothesis, there were no differences in root-to-leaf ratio between perennial and annual or C4 and C3 plants. Our study provides species-specific and functional-type-specific allometry equations for both above- and belowground organs of 11 common California grass species, enabling nondestructive biomass assessment in California grasslands. These allometric relationships and trade-offs in carbon allocation across species can improve ecosystem model predictions of grassland species interactions and environmental responses through differences in morphology.

Erratum: Publisher Correction: wMel replacement of dengue-competent mosquitoes is robust to near-term climate change.

[This corrects the article DOI: 10.1038/s41558-023-01746-w.].

Accounting for herbaceous communities in process-based models will advance our understanding of "grassy" ecosystems.

Grassland and other herbaceous communities cover significant portions of Earth's terrestrial surface and provide many critical services, such as carbon sequestration, wildlife habitat, and food production. Forecasts of global change impacts on these services will require predictive tools, such as process-based dynamic vegetation models. Yet, model representation of herbaceous communities and ecosystems lags substantially behind that of tree communities and forests. The limited representation of herbaceous communities within models arises from two important knowledge gaps: first, our empirical understanding of the principles governing herbaceous vegetation dynamics is either incomplete or does not provide mechanistic information necessary to drive herbaceous community processes with models; second, current model structure and parameterization of grass and other herbaceous plant functional types limits the ability of models to predict outcomes of competition and growth for herbaceous vegetation. In this review, we provide direction for addressing these gaps by: (1) presenting a brief history of how vegetation dynamics have been developed and incorporated into earth system models, (2) reporting on a model simulation activity to evaluate current model capability to represent herbaceous vegetation dynamics and ecosystem function, and (3) detailing several ecological properties and phenomena that should be a focus for both empiricists and modelers to improve representation of herbaceous vegetation in models. Together, empiricists and modelers can improve representation of herbaceous ecosystem processes within models. In so doing, we will greatly enhance our ability to forecast future states of the earth system, which is of high importance given the rapid rate of environmental change on our planet.

wMel replacement of dengue-competent mosquitoes is robust to near-term change.

Rising temperatures are impacting the range and prevalence of mosquito-borne diseases. A promising biocontrol technology replaces wild mosquitoes with those carrying the virus-blocking Wolbachia bacterium. Because the most widely used strain, wMel, is adversely affected by heat stress, we examined how global warming may influence wMel-based replacement. We simulated interventions in two locations with successful field trials using Coupled Model Intercomparison Project Phase 5 climate projections and historical temperature records, integrating empirical data on wMel's thermal sensitivity into a model of Aedes aegypti population dynamics to evaluate introgression and persistence over one year. We show that in Cairns, Australia, climatic futures necessitate operational adaptations for heatwaves exceeding two weeks. In Nha Trang, Vietnam, projected heatwaves of three weeks and longer eliminate wMel under the most stringent assumptions of that symbiont's thermal limits. We conclude that this technology is generally robust to near-term (2030s) climate change. Accelerated warming may challenge this in the 2050s and beyond.

Linking soil phosphorus with forest litterfall resistance and resilience to cyclone disturbance: A pantropical meta-analysis.

While tropical cyclone regimes are shifting with climate change, the mechanisms underpinning the resistance (ability to withstand disturbance-induced change) and resilience (capacity to return to pre-disturbance reference) of tropical forest litterfall to cyclones remain largely unexplored pantropically. Single-site studies in Australia and Hawaii suggest that litterfall on low-phosphorus (P) soils is more resistant and less resilient to cyclones. We conducted a meta-analysis to investigate the pantropical importance of total soil P in mediating forest litterfall resistance and resilience to 22 tropical cyclones. We evaluated cyclone-induced and post-cyclone litterfall mass (g/m2 /day), and P and nitrogen (N) fluxes (mg/m2 /day) and concentrations (mg/g), all indicators of ecosystem function and essential for nutrient cycling. Across 73 case studies in Australia, Guadeloupe, Hawaii, Mexico, Puerto Rico, and Taiwan, total litterfall mass flux increased from ~2.5 ± 0.3 to 22.5 ± 3 g/m2 /day due to cyclones, with large variation among studies. Litterfall P and N fluxes post-cyclone represented ~5% and 10% of the average annual fluxes, respectively. Post-cyclone leaf litterfall N and P concentrations were 21.6 ± 1.2% and 58.6 ± 2.3% higher than pre-cyclone means. Mixed-effects models determined that soil P negatively moderated the pantropical litterfall resistance to cyclones, with a 100 mg P/kg increase in soil P corresponding to a 32% to 38% decrease in resistance. Based on 33% of the resistance case studies, total litterfall mass flux reached pre-disturbance levels within one-year post-disturbance. A GAMM indicated that soil P, gale wind duration and time post-cyclone jointly moderate the short-term resilience of total litterfall, with the nature of the relationship between resilience and soil P contingent on time and wind duration. Across pantropical forests observed to date, our results indicate that litterfall resistance and resilience in the face of intensifying cyclones will be partially determined by total soil P.

Forest regeneration within Earth system models: current process representations and ways forward.

Earth system models must predict forest responses to global change in order to simulate future global climate, hydrology, and ecosystem dynamics. These models are increasingly adopting vegetation demographic approaches that explicitly represent tree growth, mortality, and recruitment, enabling advances in the projection of forest vulnerability and resilience, as well as evaluation with field data. To date, simulation of regeneration processes has received far less attention than simulation of processes that affect growth and mortality, in spite of their critical role maintaining forest structure, facilitating turnover in forest composition over space and time, enabling recovery from disturbance, and regulating climate-driven range shifts. Our critical review of regeneration process representations within current Earth system vegetation demographic models reveals the need to improve parameter values and algorithms for reproductive allocation, dispersal, seed survival and germination, environmental filtering in the seedling layer, and tree regeneration strategies adapted to wind, fire, and anthropogenic disturbance regimes. These improvements require synthesis of existing data, specific field data-collection protocols, and novel model algorithms compatible with global-scale simulations. Vegetation demographic models offer the opportunity to more fully integrate ecological understanding into Earth system prediction; regeneration processes need to be a critical part of the effort.

Simulating environmentally-sensitive tree recruitment in vegetation demographic models.

Vegetation demographic models (VDMs) endeavor to predict how global forests will respond to climate change. This requires simulating which trees, if any, are able to recruit under changing environmental conditions. We present a new recruitment scheme for VDMs in which functional-type-specific recruitment rates are sensitive to light, soil moisture and the productivity of reproductive trees. We evaluate the scheme by predicting tree recruitment for four tropical tree functional types under varying meteorology and canopy structure at Barro Colorado Island, Panama. We compare predictions to those of a current VDM, quantitative observations and ecological expectations. We find that the scheme improves the magnitude and rank order of recruitment rates among functional types and captures recruitment limitations in response to variable understory light, soil moisture and precipitation regimes. Our results indicate that adopting this framework will improve VDM capacity to predict functional-type-specific tree recruitment in response to climate change, thereby improving predictions of future forest distribution, composition and function.

Warming increased bark beetle-induced tree mortality by 30% during an extreme drought in California.

Quantifying the responses of forest disturbances to climate warming is critical to our understanding of carbon cycles and energy balances of the Earth system. The impact of warming on bark beetle outbreaks is complex as multiple drivers of these events may respond differently to warming. Using a novel model of bark beetle biology and host tree interactions, we assessed how contemporary warming affected western pine beetle (Dendroctonus brevicomis) populations and mortality of its host, ponderosa pine (Pinus ponderosa), during an extreme drought in the Sierra Nevada, California, United States. When compared with the field data, our model captured the western pine beetle flight timing and rates of ponderosa pine mortality observed during the drought. In assessing the influence of temperature on western pine beetles, we found that contemporary warming increased the development rate of the western pine beetle and decreased the overwinter mortality rate of western pine beetle larvae leading to increased population growth during periods of lowered tree defense. We attribute a 29.9% (95% CI: 29.4%-30.2%) increase in ponderosa pine mortality during drought directly to increases in western pine beetle voltinism (i.e., associated with increased development rates of western pine beetle) and, to a much lesser extent, reductions in overwintering mortality. These findings, along with other studies, suggest each degree (°C) increase in temperature may have increased the number of ponderosa pine killed by upwards of 35%-40% °C-1 if the effects of compromised tree defenses (15%-20%) and increased western pine beetle populations (20%) are additive. Due to the warming ability to considerably increase mortality through the mechanism of bark beetle populations, models need to consider climate's influence on both host tree stress and the bark beetle population dynamics when determining future levels of tree mortality.

Hydraulically-vulnerable trees survive on deep-water access during droughts in a tropical forest.

Deep-water access is arguably the most effective, but under-studied, mechanism that plants employ to survive during drought. Vulnerability to embolism and hydraulic safety margins can predict mortality risk at given levels of dehydration, but deep-water access may delay plant dehydration. Here, we tested the role of deep-water access in enabling survival within a diverse tropical forest community in Panama using a novel data-model approach. We inversely estimated the effective rooting depth (ERD, as the average depth of water extraction), for 29 canopy species by linking diameter growth dynamics (1990-2015) to vapor pressure deficit, water potentials in the whole-soil column, and leaf hydraulic vulnerability curves. We validated ERD estimates against existing isotopic data of potential water-access depths. Across species, deeper ERD was associated with higher maximum stem hydraulic conductivity, greater vulnerability to xylem embolism, narrower safety margins, and lower mortality rates during extreme droughts over 35 years (1981-2015) among evergreen species. Species exposure to water stress declined with deeper ERD indicating that trees compensate for water stress-related mortality risk through deep-water access. The role of deep-water access in mitigating mortality of hydraulically-vulnerable trees has important implications for our predictive understanding of forest dynamics under current and future climates.

Warming acts through earlier snowmelt to advance but not extend alpine community flowering.

Large-scale warming will alter multiple local climate factors in alpine tundra, yet very few experimental studies examine the combined yet distinct influences of earlier snowmelt, higher temperatures and altered soil moisture on alpine ecosystems. This limits our ability to predict responses to climate change by plant species and communities. To address this gap, we used infrared heaters and manual watering in a fully factorial experiment to determine the relative importance of these climate factors on plant flowering phenology, and response differences among plant functional groups. Heating advanced snowmelt and flower initiation, but exposed plants to colder early-spring conditions in the period prior to first flower, indicating that snowmelt timing, not temperature, advances flowering initiation in the alpine community. Flowering duration was largely conserved; heating did not extend average species flowering into the latter part of the growing season but instead flowering was completed earlier in heated plots. Although passive warming experiments have resulted in warming-induced soil drying suggested to advance flower senescence, supplemental water did not counteract the average species advance in flowering senescence caused by heating or extend flowering in unheated plots, and variation in soil moisture had inconsistent effects on flowering periods. Functional groups differed in sensitivity to earlier snowmelt, with flower initiation most advanced for early-season species and flowering duration lengthened only for graminoids and forbs. We conclude that earlier snowmelt, driven by increased radiative heating, is the most important factor altering alpine flowering phenology. Studies that only manipulate summer temperature will err in estimating the sensitivity of alpine flowering phenology to large-scale warming. The wholesale advance in flowering phenology with earlier snowmelt suggests that alpine communities will track warming, but only alpine forbs and graminoids appear able to take advantage of an extended snow-free season.

Pervasive shifts in forest dynamics in a changing world.

Forest dynamics arise from the interplay of environmental drivers and disturbances with the demographic processes of recruitment, growth, and mortality, subsequently driving biomass and species composition. However, forest disturbances and subsequent recovery are shifting with global changes in climate and land use, altering these dynamics. Changes in environmental drivers, land use, and disturbance regimes are forcing forests toward younger, shorter stands. Rising carbon dioxide, acclimation, adaptation, and migration can influence these impacts. Recent developments in Earth system models support increasingly realistic simulations of vegetation dynamics. In parallel, emerging remote sensing datasets promise qualitatively new and more abundant data on the underlying processes and consequences for vegetation structure. When combined, these advances hold promise for improving the scientific understanding of changes in vegetation demographics and disturbances.

Precipitation mediates sap flux sensitivity to evaporative demand in the neotropics.

Transpiration in humid tropical forests modulates the global water cycle and is a key driver of climate regulation. Yet, our understanding of how tropical trees regulate sap flux in response to climate variability remains elusive. With a progressively warming climate, atmospheric evaporative demand [i.e., vapor pressure deficit (VPD)] will be increasingly important for plant functioning, becoming the major control of plant water use in the twenty-first century. Using measurements in 34 tree species at seven sites across a precipitation gradient in the neotropics, we determined how the maximum sap flux velocity (vmax) and the VPD threshold at which vmax is reached (VPDmax) vary with precipitation regime [mean annual precipitation (MAP); seasonal drought intensity (PDRY)] and two functional traits related to foliar and wood economics spectra [leaf mass per area (LMA); wood specific gravity (WSG)]. We show that, even though vmax is highly variable within sites, it follows a negative trend in response to increasing MAP and PDRY across sites. LMA and WSG exerted little effect on vmax and VPDmax, suggesting that these widely used functional traits provide limited explanatory power of dynamic plant responses to environmental variation within hyper-diverse forests. This study demonstrates that long-term precipitation plays an important role in the sap flux response of humid tropical forests to VPD. Our findings suggest that under higher evaporative demand, trees growing in wetter environments in humid tropical regions may be subjected to reduced water exchange with the atmosphere relative to trees growing in drier climates.

A meta-analysis of 1,119 manipulative experiments on terrestrial carbon-cycling responses to global change.

Direct quantification of terrestrial biosphere responses to global change is crucial for projections of future climate change in Earth system models. Here, we synthesized ecosystem carbon-cycling data from 1,119 experiments performed over the past four decades concerning changes in temperature, precipitation, CO2 and nitrogen across major terrestrial vegetation types of the world. Most experiments manipulated single rather than multiple global change drivers in temperate ecosystems of the USA, Europe and China. The magnitudes of warming and elevated CO2 treatments were consistent with the ranges of future projections, whereas those of precipitation changes and nitrogen inputs often exceeded the projected ranges. Increases in global change drivers consistently accelerated, but decreased precipitation slowed down carbon-cycle processes. Nonlinear (including synergistic and antagonistic) effects among global change drivers were rare. Belowground carbon allocation responded negatively to increased precipitation and nitrogen addition and positively to decreased precipitation and elevated CO2. The sensitivities of carbon variables to multiple global change drivers depended on the background climate and ecosystem condition, suggesting that Earth system models should be evaluated using site-specific conditions for best uses of this large dataset. Together, this synthesis underscores an urgent need to explore the interactions among multiple global change drivers in underrepresented regions such as semi-arid ecosystems, forests in the tropics and subtropics, and Arctic tundra when forecasting future terrestrial carbon-climate feedback.

Homoeostatic maintenance of nonstructural carbohydrates during the 2015-2016 El Niño drought across a tropical forest precipitation gradient.

Nonstructural carbohydrates (NSCs) are essential for maintenance of plant metabolism and may be sensitive to short- and long-term climatic variation. NSC variation in moist tropical forests has rarely been studied, so regulation of NSCs in these systems is poorly understood. We measured foliar and branch NSC content in 23 tree species at three sites located across a large precipitation gradient in Panama during the 2015-2016 El Niño to examine how short- and long-term climatic variation impact carbohydrate dynamics. There was no significant difference in total NSCs as the drought progressed (leaf P = 0.32, branch P = 0.30) nor across the rainfall gradient (leaf P = 0.91, branch P = 0.96). Foliar soluble sugars decreased while starch increased over the duration of the dry period, suggesting greater partitioning of NSCs to storage than metabolism or transport as drought progressed. There was a large variation across species at all sites, but total foliar NSCs were positively correlated with leaf mass per area, whereas branch sugars were positively related to leaf temperature and negatively correlated with daily photosynthesis and wood density. The NSC homoeostasis across a wide range of conditions suggests that NSCs are an allocation priority in moist tropical forests.

Snowmelt Timing Regulates Community Composition, Phenology, and Physiological Performance of Alpine Plants.

The spatial patterning of alpine plant communities is strongly influenced by the variation in physical factors such as temperature and moisture, which are strongly affected by snow depth and snowmelt patterns. Earlier snowmelt timing and greater soil-moisture limitations may favor wide-ranging species adapted to a broader set of ecohydrological conditions than alpine-restricted species. We asked how plant community composition, phenology, plant water relations, and photosynthetic gas exchange of alpine-restricted and wide-ranging species differ in their responses to a ca. 40-day snowmelt gradient in the Colorado Rocky Mountains (Lewisia pygmaea, Sibbaldia procumbens, and Hymenoxys grandiflora were alpine-restricted and Artemisia scopulorum, Carex rupestris, and Geum rossii were wide-ranging species). As hypothesized, species richness and foliar cover increased with earlier snowmelt, due to a greater abundance of wide-ranging species present in earlier melting plots. Flowering initiation occurred earlier with earlier snowmelt for 12 out of 19 species analyzed, while flowering duration was shortened with later snowmelt for six species (all but one were wide-ranging species). We observed >50% declines in net photosynthesis from July to September as soil moisture and plant water potentials declined. Early-season stomatal conductance was higher in wide-ranging species, indicating a more competitive strategy for water acquisition when soil moisture is high. Even so, there were no associated differences in photosynthesis or transpiration, suggesting no strong differences between these groups in physiology. Our findings reveal that plant species with different ranges (alpine-restricted vs. wide-ranging) could have differential phenological and physiological responses to snowmelt timing and associated soil moisture dry-down, and that alpine-restricted species' performance is more sensitive to snowmelt. As a result, alpine-restricted species may serve as better indicator species than their wide-ranging heterospecifics. Overall, alpine community composition and peak % cover are strongly structured by spatio-temporal patterns in snowmelt timing. Thus, near-term, community-wide changes (or variation) in phenology and physiology in response to shifts in snowmelt timing or rates of soil dry down are likely to be contingent on the legacy of past climate on community structure.

Variation in hydroclimate sustains tropical forest biomass and promotes functional diversity.

The fate of tropical forests under climate change is unclear as a result, in part, of the uncertainty in projected changes in precipitation and in the ability of vegetation models to capture the effects of drought-induced mortality on aboveground biomass (AGB). We evaluated the ability of a terrestrial biosphere model with demography and hydrodynamics (Ecosystem Demography, ED2-hydro) to simulate AGB and mortality of four tropical tree plant functional types (PFTs) that operate along light- and water-use axes. Model predictions were compared with observations of canopy trees at Barro Colorado Island (BCI), Panama. We then assessed the implications of eight hypothetical precipitation scenarios, including increased annual precipitation, reduced inter-annual variation, El Niño-related droughts and drier wet or dry seasons, on AGB and functional diversity of the model forest. When forced with observed meteorology, ED2-hydro predictions capture multiple BCI benchmarks. ED2-hydro predicts that AGB will be sustained under lower rainfall via shifts in the functional composition of the forest, except under the drier dry-season scenario. These results support the hypothesis that inter-annual variation in mean and seasonal precipitation promotes the coexistence of functionally diverse PFTs because of the relative differences in mortality rates. If the hydroclimate becomes chronically drier or wetter, functional evenness related to drought tolerance may decline.

Comparative Transcriptomics Among Four White Pine Species.

Conifers are the dominant plant species throughout the high latitude boreal forests as well as some lower latitude temperate forests of North America, Europe, and Asia. As such, they play an integral economic and ecological role across much of the world. This study focused on the characterization of needle transcriptomes from four ecologically important and understudied North American white pines within the Pinus subgenus Strobus The populations of many Strobus species are challenged by native and introduced pathogens, native insects, and abiotic factors. RNA from the needles of western white pine (Pinus monticola), limber pine (Pinus flexilis), whitebark pine (Pinus albicaulis), and sugar pine (Pinus lambertiana) was sampled, Illumina short read sequenced, and de novo assembled. The assembled transcripts and their subsequent structural and functional annotations were processed through custom pipelines to contend with the challenges of non-model organism transcriptome validation. Orthologous gene family analysis of over 58,000 translated transcripts, implemented through Tribe-MCL, estimated the shared and unique gene space among the four species. This revealed 2025 conserved gene families, of which 408 were aligned to estimate levels of divergence and reveal patterns of selection. Specific candidate genes previously associated with drought tolerance and white pine blister rust resistance in conifers were investigated.