Fire suppression makes wildfires more severe and accentuates impacts of climate change and fuel accumulation.

Fire suppression is the primary management response to wildfires in many areas globally. By removing less-extreme wildfires, this approach ensures that remaining wildfires burn under more extreme conditions. Here, we term this the "suppression bias" and use a simulation model to highlight how this bias fundamentally impacts wildfire activity, independent of fuel accumulation and climate change. We illustrate how attempting to suppress all wildfires necessarily means that fires will burn with more severe and less diverse ecological impacts, with burned area increasing at faster rates than expected from fuel accumulation or climate change. Over a human lifespan, the modeled impacts of the suppression bias exceed those from fuel accumulation or climate change alone, suggesting that suppression may exert a significant and underappreciated influence on patterns of fire globally. Managing wildfires to safely burn under low and moderate conditions is thus a critical tool to address the growing wildfire crisis.

Latitudinal patterns in stabilizing density dependence of forest communities.

Numerous studies have shown reduced performance in plants that are surrounded by neighbours of the same species1,2, a phenomenon known as conspecific negative density dependence (CNDD)3. A long-held ecological hypothesis posits that CNDD is more pronounced in tropical than in temperate forests4,5, which increases community stabilization, species coexistence and the diversity of local tree species6,7. Previous analyses supporting such a latitudinal gradient in CNDD8,9 have suffered from methodological limitations related to the use of static data10-12. Here we present a comprehensive assessment of latitudinal CNDD patterns using dynamic mortality data to estimate species-site-specific CNDD across 23 sites. Averaged across species, we found that stabilizing CNDD was present at all except one site, but that average stabilizing CNDD was not stronger toward the tropics. However, in tropical tree communities, rare and intermediate abundant species experienced stronger stabilizing CNDD than did common species. This pattern was absent in temperate forests, which suggests that CNDD influences species abundances more strongly in tropical forests than it does in temperate ones13. We also found that interspecific variation in CNDD, which might attenuate its stabilizing effect on species diversity14,15, was high but not significantly different across latitudes. Although the consequences of these patterns for latitudinal diversity gradients are difficult to evaluate, we speculate that a more effective regulation of population abundances could translate into greater stabilization of tropical tree communities and thus contribute to the high local diversity of tropical forests.

Mycorrhizal feedbacks influence global forest structure and diversity.

One mechanism proposed to explain high species diversity in tropical systems is strong negative conspecific density dependence (CDD), which reduces recruitment of juveniles in proximity to conspecific adult plants. Although evidence shows that plant-specific soil pathogens can drive negative CDD, trees also form key mutualisms with mycorrhizal fungi, which may counteract these effects. Across 43 large-scale forest plots worldwide, we tested whether ectomycorrhizal tree species exhibit weaker negative CDD than arbuscular mycorrhizal tree species. We further tested for conmycorrhizal density dependence (CMDD) to test for benefit from shared mutualists. We found that the strength of CDD varies systematically with mycorrhizal type, with ectomycorrhizal tree species exhibiting higher sapling densities with increasing adult densities than arbuscular mycorrhizal tree species. Moreover, we found evidence of positive CMDD for tree species of both mycorrhizal types. Collectively, these findings indicate that mycorrhizal interactions likely play a foundational role in global forest diversity patterns and structure.

Reduced fire severity offers near-term buffer to climate-driven declines in conifer resilience across the western United States.

Increasing fire severity and warmer, drier postfire conditions are making forests in the western United States (West) vulnerable to ecological transformation. Yet, the relative importance of and interactions between these drivers of forest change remain unresolved, particularly over upcoming decades. Here, we assess how the interactive impacts of changing climate and wildfire activity influenced conifer regeneration after 334 wildfires, using a dataset of postfire conifer regeneration from 10,230 field plots. Our findings highlight declining regeneration capacity across the West over the past four decades for the eight dominant conifer species studied. Postfire regeneration is sensitive to high-severity fire, which limits seed availability, and postfire climate, which influences seedling establishment. In the near-term, projected differences in recruitment probability between low- and high-severity fire scenarios were larger than projected climate change impacts for most species, suggesting that reductions in fire severity, and resultant impacts on seed availability, could partially offset expected climate-driven declines in postfire regeneration. Across 40 to 42% of the study area, we project postfire conifer regeneration to be likely following low-severity but not high-severity fire under future climate scenarios (2031 to 2050). However, increasingly warm, dry climate conditions are projected to eventually outweigh the influence of fire severity and seed availability. The percent of the study area considered unlikely to experience conifer regeneration, regardless of fire severity, increased from 5% in 1981 to 2000 to 26 to 31% by mid-century, highlighting a limited time window over which management actions that reduce fire severity may effectively support postfire conifer regeneration.

Demographic composition, not demographic diversity, predicts biomass and turnover across temperate and tropical forests.

The growth and survival of individual trees determine the physical structure of a forest with important consequences for forest function. However, given the diversity of tree species and forest biomes, quantifying the multitude of demographic strategies within and across forests and the way that they translate into forest structure and function remains a significant challenge. Here, we quantify the demographic rates of 1961 tree species from temperate and tropical forests and evaluate how demographic diversity (DD) and demographic composition (DC) differ across forests, and how these differences in demography relate to species richness, aboveground biomass (AGB), and carbon residence time. We find wide variation in DD and DC across forest plots, patterns that are not explained by species richness or climate variables alone. There is no evidence that DD has an effect on either AGB or carbon residence time. Rather, the DC of forests, specifically the relative abundance of large statured species, predicted both biomass and carbon residence time. Our results demonstrate the distinct DCs of globally distributed forests, reflecting biogeography, recent history, and current plot conditions. Linking the DC of forests to resilience or vulnerability to climate change, will improve the precision and accuracy of predictions of future forest composition, structure, and function.

Arbuscular mycorrhizal trees influence the latitudinal beta-diversity gradient of tree communities in forests worldwide.

Arbuscular mycorrhizal (AM) and ectomycorrhizal (EcM) associations are critical for host-tree performance. However, how mycorrhizal associations correlate with the latitudinal tree beta-diversity remains untested. Using a global dataset of 45 forest plots representing 2,804,270 trees across 3840 species, we test how AM and EcM trees contribute to total beta-diversity and its components (turnover and nestedness) of all trees. We find AM rather than EcM trees predominantly contribute to decreasing total beta-diversity and turnover and increasing nestedness with increasing latitude, probably because wide distributions of EcM trees do not generate strong compositional differences among localities. Environmental variables, especially temperature and precipitation, are strongly correlated with beta-diversity patterns for both AM trees and all trees rather than EcM trees. Results support our hypotheses that latitudinal beta-diversity patterns and environmental effects on these patterns are highly dependent on mycorrhizal types. Our findings highlight the importance of AM-dominated forests for conserving global forest biodiversity.

Alpine treeline ecotones are potential refugia for a montane pine species threatened by bark beetle outbreaks.

Warming-induced mountain pine beetle (Dendroctonus ponderosae; MPB) outbreaks have caused extensive mortality of whitebark pine (Pinus albicaulis; WBP) throughout the species' range. In the highest mountains where WBP occur, they cross alpine treeline ecotones (ATEs) where growth forms transition from trees to shrub-like krummholz, some of which survived recent MPB outbreaks. This observation motivated the hypothesis that ATEs are refugia for WBP because krummholz growth forms escape MPB attack and have the potential to produce viable seed. To test this hypothesis, we surveyed WBP mortality along transects from the ATE edge (locally highest krummholz WBP) downslope into the forest and, to distinguish if survival mechanisms are unique to ATEs, across other forest ecotones (OFEs) from the edge of WBP occurrence into the forest. We replicated this design at 10 randomly selected sites in the U.S. Northern Rocky Mountains. We also surveyed reproduction in a subset of ATE sites. Mortality was nearly absent in upper ATEs (mean ± SE percent dead across all sites of 0.03% ± 0.03% 0-100 m from the edge and 14.1% ± 1.7% 100-500 m from the edge) but was above 20% along OFEs (21.4 ± 5.2% 0-100 m and 32.4 ± 2.7% 100-500 m from the edge). We observed lower reproduction in upper ATEs (16 ± 9.9 cones/ha and 12.9 ± 5.3 viable seeds/cone 0-100 m from the edge) compared to forests below (317.1 ± 64.4 cones/ha and 32.5 ± 2.5 viable seeds/cone 100-500 m from the edge). Uniquely high WBP survival supports the hypothesis that ATEs serve as refugia because krummholz growth forms escape MPB attack. However, low reproduction suggests ATE refugia function over longer time periods. Beyond our WBP system, we propose that plant populations in marginal environments are candidate refugia if distinct phenotypes result in reduced disturbance impacts.

Complications following symptom-limited thoracentesis using suction.

BACKGROUND: Thoracentesis using suction is perceived to have increased risk of complications, including pneumothorax and re-expansion pulmonary oedema (REPO). Current guidelines recommend limiting drainage to 1.5 L to avoid REPO. Our purpose was to examine the incidence of complications with symptom-limited drainage of pleural fluid using suction and identify risk factors for REPO. METHODS: A retrospective cohort study of all adult patients who underwent symptom-limited thoracentesis using suction at our institution between January 1, 2004 and August 31, 2018 was performed, and a total of 10 344 thoracenteses were included. RESULTS: Pleural fluid ≥1.5 L was removed in 19% of the procedures. Thoracentesis was stopped due to chest discomfort (39%), complete drainage of fluid (37%) and persistent cough (13%). Pneumothorax based on chest radiography was detected in 3.98%, but only 0.28% required intervention. The incidence of REPO was 0.08%. The incidence of REPO increased with Eastern Cooperative Oncology Group performance status (ECOG PS) ≥3 compounded with ≥1.5 L (0.04-0.54%; 95% CI 0.13-2.06 L). Thoracentesis in those with ipsilateral mediastinal shift did not increase complications, but less fluid was removed (p<0.01). CONCLUSIONS: Symptom-limited thoracentesis using suction is safe even with large volumes. Pneumothorax requiring intervention and REPO are both rare. There were no increased procedural complications in those with ipsilateral mediastinal shift. REPO increased with poor ECOG PS and drainage ≥1.5 L. Symptom-limited drainage using suction without pleural manometry is safe.

Vegetation dynamics following compound disturbance in a dry pine forest: fuel treatment then bark beetle outbreak.

In the western United States, restoration of forests with historically frequent, low-severity fire regimes often includes fuel reduction that reestablish open, early-seral conditions while reducing fuel continuity and loading. Between 2001 and 2016, fuel reduction (e.g., thinning, prescribed burning, etc.) was implemented on over 26 million hectares of federal lands alone in the United States, reflecting the urgency to mitigate risk from high-severity wildfire. However, between 2001 and 2012, nearly 20 million hectares in the United States were impacted by mountain pine beetle (MPB; Dendroctonus ponderosae), compounding restoration effects in wildfire-hazard-treated stands. Knowledge of the effects of treatments followed by natural disturbance on long-term forest structure and communities is needed, especially considering that fuel treatments are increasingly being implemented and warming climate is predicted to exacerbate disturbance frequency and severity. We tested the interacting effects of treatments designed to reduce high-severity wildfire hazard in stands subsequently challenged by MPB outbreak on vegetation dynamics using a factorial experimental design (control, thin only, burn only, thin + burn) in a ponderosa pine (Pinus ponderosa)-dominated forest. Stands were treated by 2002, then impacted by MPB outbreak from 2005 to 2012. We assessed change in overstory and understory forest community structure, composition, and diversity over time. There were distinct thinning, burning, and year effects. Thinning immediately reduced overstory density; pine density then declined 4.5 times more in unthinned than thinned treatments due to MPB. Burning immediately reduced graminoid, shrub, and total understory cover by as much as 52%, resulting in greater species evenness than unburned treatments, but differences disappeared by 2016 due to growth and MPB outbreak. Similarly, multivariate analyses indicated forest communities were starkly different after treatment but became more similar over time, though key understory and overstory attributes still distinguish control and thin + burn. This study shows the value of long-term silvicultural experiments to evaluate treatment longevity and the compounded effects of treatment and natural disturbance. We demonstrate the homogenizing effects of treatment-induced growth coupled with MPB-caused tree mortality on management strategies that just treat the overstory (thinning) or understory (burning), showing that only combined treatments can provide the unique structural and compositional outcomes expected of restoration.

Fixing a snag in carbon emissions estimates from wildfires.

Wildfire is an essential earth-system process, impacting ecosystem processes and the carbon cycle. Forest fires are becoming more frequent and severe, yet gaps exist in the modeling of fire on vegetation and carbon dynamics. Strategies for reducing carbon dioxide (CO2 ) emissions from wildfires include increasing tree harvest, largely based on the public assumption that fires burn live forests to the ground, despite observations indicating that less than 5% of mature tree biomass is actually consumed. This misconception is also reflected though excessive combustion of live trees in models. Here, we show that regional emissions estimates using widely implemented combustion coefficients are 59%-83% higher than emissions based on field observations. Using unique field datasets from before and after wildfires and an improved ecosystem model, we provide strong evidence that these large overestimates can be reduced by using realistic biomass combustion factors and by accurately quantifying biomass in standing dead trees that decompose over decades to centuries after fire ("snags"). Most model development focuses on area burned; our results reveal that accurately representing combustion is also essential for quantifying fire impacts on ecosystems. Using our improvements, we find that western US forest fires have emitted 851 ± 228 Tg CO2 (~half of alternative estimates) over the last 17 years, which is minor compared to 16,200 Tg CO2 from fossil fuels across the region.

Response to Comment on "Plant diversity increases with the strength of negative density dependence at the global scale".

Hülsmann and Hartig suggest that ecological mechanisms other than specialized natural enemies or intraspecific competition contribute to our estimates of conspecific negative density dependence (CNDD). To address their concern, we show that our results are not the result of a methodological artifact and present a null-model analysis that demonstrates that our original findings-(i) stronger CNDD at tropical relative to temperate latitudes and (ii) a latitudinal shift in the relationship between CNDD and species abundance-persist even after controlling for other processes that might influence spatial relationships between adults and recruits.

Response to Comment on "Plant diversity increases with the strength of negative density dependence at the global scale".

Chisholm and Fung claim that our method of estimating conspecific negative density dependence (CNDD) in recruitment is systematically biased, and present an alternative method that shows no latitudinal pattern in CNDD. We demonstrate that their approach produces strongly biased estimates of CNDD, explaining why they do not detect a latitudinal pattern. We also address their methodological concerns using an alternative distance-weighted approach, which supports our original findings of a latitudinal gradient in CNDD and a latitudinal shift in the relationship between CNDD and species abundance.

Plant diversity increases with the strength of negative density dependence at the global scale.

Theory predicts that higher biodiversity in the tropics is maintained by specialized interactions among plants and their natural enemies that result in conspecific negative density dependence (CNDD). By using more than 3000 species and nearly 2.4 million trees across 24 forest plots worldwide, we show that global patterns in tree species diversity reflect not only stronger CNDD at tropical versus temperate latitudes but also a latitudinal shift in the relationship between CNDD and species abundance. CNDD was stronger for rare species at tropical versus temperate latitudes, potentially causing the persistence of greater numbers of rare species in the tropics. Our study reveals fundamental differences in the nature of local-scale biotic interactions that contribute to the maintenance of species diversity across temperate and tropical communities.

Evaluating a new method for reconstructing forest conditions from General Land Office survey records.

Historical forest conditions are often used to inform contemporary management goals because historical forests are considered to be resilient to ecological disturbances. The General Land Office (GLO) surveys of the late 19th and early 20th centuries provide regionally quasi-contiguous data sets of historical forests across much of the Western United States. Multiple methods exist for estimating tree density from point-based sampling such as the GLO surveys, including distance-based and area-based approaches. Area-based approaches have been applied in California mixed-conifer forests but their estimates have not been validated. To assess the accuracy and precision of plotless density estimators with potential for application to GLO data in this region, we imposed a GLO sampling scheme on six mapped forest stands of known densities (159-784 trees/ha) in the Sierra Nevada in California, USA, and Baja California Norte, Mexico. We compared three distance-based plotless density estimators (Cottam, Pollard, and Morisita) as well as two Voronoi area (VA) estimators, the Delincé and mean harmonic Voronoi density (MHVD), to the true densities. We simulated sampling schemes of increasing intensity to assess sampling error. The relative error (RE) of density estimates for the GLO sampling scheme ranged from 0.36 to 4.78. The least biased estimate of tree density in every stand was obtained with the Morisita estimator and the most biased was obtained with the MHVD estimator. The MHVD estimates of tree density were 1.2-3.8 times larger than the true densities and performed best in stands subject to fire exclusion for 100 yr. The Delincé approach obtained accurate estimates of density, implying that the Voronoi approach is theoretically sound but that its application in the MHVD was flawed. The misapplication was attributed to two causes: (1) the use of a crown scaling factor that does not correct for the number of trees sampled and (2) the persistent underestimate of the true VA due to a weak relationship between tree size and VA. The magnitude of differences between true densities and MHVD estimates suggest caution in using results based on the MHVD to inform management and restoration practices in the conifer forests of the American West.

Cutaneous cryptococcosis in a patient taking fingolimod for multiple sclerosis: Here come the opportunistic infections?

BACKGROUND: Fingolimod is an oral disease-modifying therapy for relapsing forms of multiple sclerosis, which acts by sequestering lymphocytes within lymph nodes. OBJECTIVE: To describe a case of extrapulmonary cryptococcosis in a patient taking fingolimod. METHODS: Case report. RESULTS: A 47-year-old man developed a non-healing skin lesion approximately 16 months after starting treatment with fingolimod. Biopsy revealed cryptococcosis. Fingolimod was discontinued and the lesion resolved with antifungal therapy. CONCLUSION: Despite few reported opportunistic infections in the pivotal clinical trials and first few years post-marketing, there has been a recent increase in reported AIDS-defining illnesses in patients taking fingolimod. Neurologists should be alert for opportunistic infections in their patients using this medication.

Spatial aspects of tree mortality strongly differ between young and old-growth forests.

Rates and spatial patterns of tree mortality are predicted to change during forest structural development. In young forests, mortality should be primarily density dependent due to competition for light, leading to an increasingly spatially uniform pattern of surviving trees. In contrast, mortality in old-growth forests should be primarily caused by contagious and spatially autocorrelated agents (e.g., insects, wind), causing spatial aggregation of surviving trees to increase through time. We tested these predictions by contrasting a three-decade record of tree mortality from replicated mapped permanent plots located in young (< 60-year-old) and old-growth (> 300-year-old) Abies amabilis forests. Trees in young forests died at a rate of 4.42% per year, whereas trees in old-growth forests died at 0.60% per year. Tree mortality in young forests was significantly aggregated, strongly density dependent, and caused live tree patterns to become more uniform through time. Mortality in old-growth forests was spatially aggregated, but was density independent and did not change the spatial pattern of surviving trees. These results extend current theory by demonstrating that density-dependent competitive mortality leading to increasingly uniform tree spacing in young forests ultimately transitions late in succession to a more diverse tree mortality regime that maintains spatial heterogeneity through time.

Comparative evolutionary diversity and phylogenetic structure across multiple forest dynamics plots: a mega-phylogeny approach.

Forest dynamics plots, which now span longitudes, latitudes, and habitat types across the globe, offer unparalleled insights into the ecological and evolutionary processes that determine how species are assembled into communities. Understanding phylogenetic relationships among species in a community has become an important component of assessing assembly processes. However, the application of evolutionary information to questions in community ecology has been limited in large part by the lack of accurate estimates of phylogenetic relationships among individual species found within communities, and is particularly limiting in comparisons between communities. Therefore, streamlining and maximizing the information content of these community phylogenies is a priority. To test the viability and advantage of a multi-community phylogeny, we constructed a multi-plot mega-phylogeny of 1347 species of trees across 15 forest dynamics plots in the ForestGEO network using DNA barcode sequence data (rbcL, matK, and psbA-trnH) and compared community phylogenies for each individual plot with respect to support for topology and branch lengths, which affect evolutionary inference of community processes. The levels of taxonomic differentiation across the phylogeny were examined by quantifying the frequency of resolved nodes throughout. In addition, three phylogenetic distance (PD) metrics that are commonly used to infer assembly processes were estimated for each plot [PD, Mean Phylogenetic Distance (MPD), and Mean Nearest Taxon Distance (MNTD)]. Lastly, we examine the partitioning of phylogenetic diversity among community plots through quantification of inter-community MPD and MNTD. Overall, evolutionary relationships were highly resolved across the DNA barcode-based mega-phylogeny, and phylogenetic resolution for each community plot was improved when estimated within the context of the mega-phylogeny. Likewise, when compared with phylogenies for individual plots, estimates of phylogenetic diversity in the mega-phylogeny were more consistent, thereby removing a potential source of bias at the plot-level, and demonstrating the value of assessing phylogenetic relationships simultaneously within a mega-phylogeny. An unexpected result of the comparisons among plots based on the mega-phylogeny was that the communities in the ForestGEO plots in general appear to be assemblages of more closely related species than expected by chance, and that differentiation among communities is very low, suggesting deep floristic connections among communities and new avenues for future analyses in community ecology.

Spatially nonrandom tree mortality and ingrowth maintain equilibrium pattern in an old-growth Pseudotsuga-Tsuga forest.

Mortality processes in old-growth forests are generally assumed to be driven by gap-scale disturbance, with only a limited role ascribed to density-dependent mortality, but these assumptions are rarely tested with data sets incorporating repeated measurements. Using a 12-ha spatially explicit plot censused 13 years apart in an approximately 500-year-old Pseudotsuga-Tsuga forest, we demonstrate significant density-dependent mortality and spatially aggregated tree recruitment. However, the combined effect of these strongly nonrandom demographic processes was to maintain tree patterns in a state of dynamic equilibrium. Density-dependent mortality was most pronounced for the dominant late-successional species, Tsuga heterophylla. The long-lived, early-seral Pseudotsuga menziesii experienced an annual stem mortality rate of 0.84% and no new recruitment. Late-seral species Tsuga and Abies amabilis had nearly balanced demographic rates of ingrowth and mortality. The 2.34% mortality rate for Taxus brevifolia was higher than expected, notably less than ingrowth, and strongly affected by proximity to Tsuga. Large-diameter Tsuga structured both the regenerating conspecific and heterospecific cohorts with recruitment of Tsuga and Abies unlikely in neighborhoods crowded with large-diameter competitors (P < 0.001). Density-dependent competitive interactions strongly shape forest communities even five centuries after stand initiation, underscoring the dynamic nature of even equilibrial old-growth forests.

The importance of large-diameter trees to forest structural heterogeneity.

Large-diameter trees dominate the structure, dynamics and function of many temperate and tropical forests. However, their attendant contributions to forest heterogeneity are rarely addressed. We established the Wind River Forest Dynamics Plot, a 25.6 ha permanent plot within which we tagged and mapped all 30,973 woody stems ≥ 1 cm dbh, all 1,966 snags ≥ 10 cm dbh, and all shrub patches ≥ 2 m(2). Basal area of the 26 woody species was 62.18 m(2)/ha, of which 61.60 m(2)/ha was trees and 0.58 m(2)/ha was tall shrubs. Large-diameter trees (≥ 100 cm dbh) comprised 1.5% of stems, 31.8% of basal area, and 17.6% of the heterogeneity of basal area, with basal area dominated by Tsuga heterophylla and Pseudotsuga menziesii. Small-diameter subpopulations of Pseudotsuga menziesii, Tsuga heterophylla and Thuja plicata, as well as all tree species combined, exhibited significant aggregation relative to the null model of complete spatial randomness (CSR) up to 9 m (P ≤ 0.001). Patterns of large-diameter trees were either not different from CSR (Tsuga heterophylla), or exhibited slight aggregation (Pseudotsuga menziesii and Thuja plicata). Significant spatial repulsion between large-diameter and small-diameter Tsuga heterophylla suggests that large-diameter Tsuga heterophylla function as organizers of tree demography over decadal timescales through competitive interactions. Comparison among two forest dynamics plots suggests that forest structural diversity responds to intermediate-scale environmental heterogeneity and disturbances, similar to hypotheses about patterns of species richness, and richness- ecosystem function. Large mapped plots with detailed within-plot environmental spatial covariates will be required to test these hypotheses.