Empirical orthogonal functions (EOFs) applied to gridded Earth system data enables users to diagnose modes of variability with relative ease. Yet, many challenges to interpretation exist such that they must be used with awareness and intention when applied to gridded climate data, especially with large ensembles. Utilizing data from two different Earth system modelling large ensemble frameworks, the Energy Exoscale Earth System Model and the Community Earth System Model, as well as reanalysis data, common EOF pitfalls are summarized and discussed. Challenges include erroneous mode swapping, sign flipping, and the temporal variability of the centers of action. For modes of variability with similar contribution to variance, mode swapping is not uncommon. Sign flipping can occur with almost any mode where the pattern is correct, but the sign is arbitrary. Although the variability of the center of action is not necessarily problematic, it potentially complicates interpretation over multi-century timescales. A wide variety of alternative methods to EOFs exist, but fitness-for-purpose must be evaluated. Additionally, illustrations of alternative methods and examples of proper use are provided. Alternative methods fit into three categories: EOF variants, linear methods, and multilinear methods.

As of 2023, the use of prescribed fire to manage ecosystems accounts for more than 50% of area burned annually across the United States. Prescribed fire is carried out when meteorological conditions, including temperature, humidity, and wind speed are appropriate for its safe and effective application. However, changes in these meteorological variables associated with future climate change may impact future opportunities to conduct prescribed fire. In this study, we combine climate projections with information on prescribed burning windows for ecoregions across the contiguous United States (CONUS) to compute the number of days when meteorological conditions allow for the safe and effective application of prescribed fire under present-day (2006–2015) and future climate (2051–2060) conditions. The resulting projections, which cover 57% of all vegetated area across the CONUS, indicate fewer days with conditions suitable for prescribed burning across ecoregions of the eastern United States due to rising maximum daily temperatures, but opportunities increase in the northern and northwestern United States, driven primarily by rising minimum temperatures and declining wind speeds.

Future climate presents conflicting implications for forest biomass. We evaluate how plant hydraulic traits, elevated CO2 levels, warming, and changes in precipitation affect forest primary productivity, evapotranspiration, and the risk of hydraulic failure. We used a dynamic vegetation model with plant hydrodynamics (FATES-HYDRO) to simulate the stand-level responses to future climate changes in a wet tropical forest in Barro Colorado Island, Panama. We calibrated the model by selecting plant trait assemblages that performed well against observations. These assemblages were run with temperature and precipitation changes for two greenhouse gas emission scenarios (2086-2100: SSP2-45, SSP5-85) and two CO2 levels (contemporary, anticipated). The risk of hydraulic failure is projected to increase from a contemporary rate of 5.7% to 10.1-11.3% under future climate scenarios, and, crucially, elevated CO2 provided only slight amelioration. By contrast, elevated CO2 mitigated GPP reductions. We attribute a greater variation in hydraulic failure risk to trait assemblages than to either CO2 or climate. Our results project forests with both faster growth (through productivity increases) and higher mortality rates (through increasing rates of hydraulic failure) in the neo-tropics accompanied by certain trait plant assemblages becoming nonviable.

Mesoscale-to-microscale coupling is an important tool to conduct turbulence-resolving multiscale simulations of realistic atmospheric flows, which are crucial for applications ranging from wind energy to wildfire spread studies. Different techniques are used to facilitate the development of realistic turbulence in the large-eddy simulation (LES) domain while minimizing computational cost. Here, we explore the impact of a simple and computationally efficient Stochastic Cell Perturbation method using momentum perturbation (SCPM-M) to accelerate turbulence generation in boundary-coupled LES simulations using the Weather Research and Forecasting (WRF) model. We simulate a convective boundary layer (CBL) to characterize the production and dissipation of turbulent kinetic energy (TKE) and the variation of TKE budget terms. Furthermore, we evaluate the impact of applying momentum perturbations of three magnitudes below, up to, and above the CBL on the TKE budget terms. Momentum perturbations greatly reduce the fetch associated with turbulence generation. When applied to half the vertical extent of the boundary layer, momentum 1 perturbations produce an adequate amount of turbulence. However, when applied above the CBL, additional structures are generated at the top of the CBL, near the inversion layer. The magnitudes of the TKE budgets produced by SCPM-M when applied at varying heights and with different perturbation amplitudes are always higher near the surface and inversion layer than those produced by No-SCPM, as are their contributions to the TKE. This study provides a better understanding of how SCPM-M reduces computational costs and how different budget terms contribute to TKE in a boundary-coupled LES simulation.

Arctic amplification (AA), characterized by a more rapid surface air temperature (SAT) warming in the Arctic than the global average, is a major feature of global climate warming. Various metrics have been used to quantify AA based on SAT anomalies, trends, or variability, and they can yield quite different conclusions regarding the magnitude and temporal patterns of AA. This study examines and compares various AA metrics for their temporal consistency in the region north of 70 degrees N from the early twentieth to the early 21st century using observational data and reanalysis products. We also quantify contributions of different radiative feedback mechanisms to AA based on short-term climate variability in reanalysis and model data using the Kernel-Gregory approach. Albedo and lapse rate feedbacks are positive and comparable, with albedo feedback being the leading contributor for all AA metrics. The net cloud feedback, which has large uncertainties, depends strongly on the data sets and AA metrics used. By quantifying the influence of internal variability on AA and related feedbacks based on global climate model ensemble simulations, we find that water vapor and cloud feedbacks are most heavily affected by internal variability.

Wildfires are a global crisis, but current fire models fail to capture vegetation response to changing climate. With drought and elevated temperature increasing the importance of vegetation dynamics to fire behavior, and the advent of next generation models capable of capturing increasingly complex physical processes, we provide a renewed focus on representation of woody vegetation in fire models. Currently, the most advanced representations of fire behavior and biophysical fire effects are found in distinct classes of fine-scale models and do not capture variation in live fuel (i.e. living plant) properties. We demonstrate that plant water and carbon dynamics, which influence combustion and heat transfer into the plant and often dictate plant survival, provide the mechanistic linkage between fire behavior and effects. Our conceptual framework linking remotely sensed estimates of plant water and carbon to fine-scale models of fire behavior and effects could be a critical first step toward improving the fidelity of the coarse scale models that are now relied upon for global fire forecasting. This process-based approach will be essential to capturing the influence of physiological responses to drought and warming on live fuel conditions, strengthening the science needed to guide fire managers in an uncertain future.

Outbreaks of several bark beetle species can develop rapidly in response to drought and may result in large transfers of carbon (C) stored in live trees to C stored in dead trees (10s of Tg C yr-1 in the western U.S. alone), which over time will be released back to the atmosphere. The western pine beetle (WPB) outbreak incited by the 2012–2015 mega-drought in the Sierra Nevada, California, U.S., could portend more frequent and/or severe bark beetle outbreaks as the temperature warms and drought frequency and intensity increase in the future. However, changes in the frequency and/or severity (resultant levels of host tree mortality) of beetle outbreaks are difficult to predict as outbreaks are complex with non-linear and eruptive processes primarily driven by interactions among beetle populations, the demography of hosts and other tree species, and climate and weather. Using an insect phenology and tree defense model, we projected the future likelihood of WPB outbreaks in the Sierra Nevada with climate drivers from different Earth System Models. Our goal was to understand how host (ponderosa pine, PIPO) recovery and future warming and drought affect the frequency and severity of WPB outbreaks and their C consequences. Our projections suggested that by 2100 the C stored in live PIPO (mean: 1.98 kg C m-2, 95% CI: 1.74–2.21 kg C m-2) will not return to levels that occurred before the 2012–2015 drought (2012: ∼2.30 kg C m-2) due to future WPB outbreaks. However, differences in climate models indicate a wide range of possible WPB outbreak frequencies and severities. Our results suggest that total plot basal area is the most significant factor in the mortality rate of PIPO by WPB in any given year, followed by drought severity and temperature. High levels of host basal area, higher temperature, and extreme drought all contribute to the frequency and severity of future WPB outbreaks. While PIPO basal area may decline under increased drought and warming, limiting high-stand basal area (>60 m2 ha-1) may reduce the severity of future WPB outbreaks in the Sierra Nevada.

Heterogeneity in surface fuels produced by overstory trees and understory vegetation is a major driver of fire behavior and ecosystem dynamics. Previous attempts at predicting tree leaf and needle litter accumulation over time have been constrained in scope to probabilistic models that consider a limited number of key factors influencing tree litter dispersal patterns and decomposition processes. We present a mechanistic model for estimating variation in surface fuels called the Distribution of Understory using Elliptical Transport (DUET). DUET uses a pre-generated voxelated canopy array and distributes the leaf litter from the trees in an elliptical shape based on tree species characteristics, wind data, and location-specific features. DUET then calculates grass growth between trees and around the litter patterns and sets moisture levels based on relative humidity. In this study, we use DUET to simulate spatial distributions of leaf litter and illustrate model sensitivity to several parameters, providing inputs for testing fuel spatial heterogeneity impacts on fire behavior using process-based models.

local wind velocity. The second conclusion is that HIGRAD/FIRETEC can be used to develop point-functions for use in HIGRAD/BEHAVE.

Structural diversity is an emerging dimension of biodiversity that accounts for size variations in organs among individuals in a community. Previous studies show significant effects of structural diversity on forest growth, but its effects on forest mortality are not known, particularly at a large scale. To address this knowledge gap, we quantified structural diversity using stem structural diversity (SSD) based on both tree diameter and height. We obtained U.S. Forest Service Forest Inventory and Analysis (FIA) data from over 2400 plots across southcentral U.S. forests that have suffered a recent drought. Using data from multiple sampling times, we calculated SSD and compared the relative importance of SSD, species diversity, functional diversity and other stand attributes in determining tree mortalities caused by fire, insects and diseases. We also used FIRETEC, a physics‐based fire model, to test the effect of SSD on canopy consumption by fire. Our results showed that (1) SSD was positively associated with tree mortalities caused by all three disturbances; (2) species richness was negatively associated with insect‐ and disease‐caused mortalities; (3) functional diversity was negatively associated with fire‐ and disease‐caused mortalities and (4) more phylogenetically related species had more similar mortality rates by insect and disease but not fire. Moreover, the FIRETEC model showed increasing canopy consumption by fire in stands with greater SSD. Together, the different tree mortalities during drought associated with SSD more consistently than the other biodiversity metrics were evaluated. Synthesis. Our results suggest that SSD could be considered in modelling forest dynamics and planning management to sustain forest health under disturbances.

Mesoscale-to-Large Eddy Simulation (LES) grid nesting is an important tool for many atmospheric model applications, ranging from wind energy to wildfire spread studies. Different techniques are used in such applications to accelerate the development of turbulence in the LES domain. Here, we explore the impact of a simple and computationally efficient Stochastic Cell Perturbation method (SCPM) to accelerate the generation of turbulence in the Weather Research and Forecasting (WRF) LES model on the Turbulence Kinetic Energy (TKE) budget. In a convective boundary layer, we study the variation of TKE budget terms under the initial conditions of the Scaled Wind Farm Technology (SWiFT) facility located in West Texas. In this study, WRF LES is used with a horizontal grid resolution of 12 m, and is one-way nested within an idealized mesoscale domain. It is crucial to understand how forced perturbation shifts the balance between the terms of the TKE budget. Here, we quantify the shear production, and buoyant production in an unstable case. Since additional production terms are introduced in the SCPM method, we investigate the dissipation term of TKE. In addition, we also study the generation of turbulent transport. Generally, it integrates over height to null in a planar homogeneous case without subsidence, indicating it is positive over some heights and negative over other heights. Furthermore, we also study the variation of the TKE transport term after extending the random perturbation up to a certain height. The findings of this study will provide a better understanding of the contribution of different budget terms in a forced LES simulation.

contrasted to predict the fate of whitebark pine over 200 years. Model testing reveals predictions within 10 to 30 percent of observed values.

Coupled fire-atmosphere models are increasingly being used to study low-intensity fires, such as those that are used in prescribed fire applications. Thus, the need arises to evaluate these models for their ability to accurately represent fire spread in marginal burning conditions. In this study, wind and fuel data collected during the Prescribed Fire Combustion and Atmospheric Dynamics Research Experiments (RxCADRE) fire campaign were used to generate initial and boundary conditions for coupled fire-atmosphere simulations. We present a novel method to obtain fuels representation at the model grid scale using a combination of imagery, machine learning, and field sampling. Several methods to generate wind input conditions for the model from eight different anemometer measurements are explored. We find a strong sensitivity of fire outcomes to wind inputs. This result highlights the critical need to include variable wind fields as inputs in modeling marginal fire conditions. This work highlights the complexities of comparing physics-based model results against observations, which are more acute in marginal burning conditions, where stronger sensitivities to local variability in wind and fuels drive fire outcomes.

Atmospheric forcing and interactions between the fire and atmosphere are primary drivers of wildland fire behavior. The atmosphere is known to be a chaotic system that, although deterministic, is very sensitive to small perturbations to initial conditions. We assume that as a result of the tight coupling between fire and atmosphere; wildland fire behavior, in turn, should also be sensitive to perturbations in atmospheric initial conditions. Observations suggest that low intensity prescribed fire, in particular, is susceptible to small perturbations in the wind field, which can significantly alter fire spread. Here, we employ a computational fluid dynamics model of coupled fire‐atmosphere interactions to answer the question: How sensitive is fire behavior to small variations in atmospheric turbulence? We perform ensemble simulations of fires in homogenous grass fuels. The only difference between ensemble members is the state of the turbulent atmosphere provided to the model throughout the simulation. The atmospheric state is a function of the initial conditions applied at the start of the simulation and boundary conditions applied throughout the simulation. We find a wide range of outcomes, with area burned ranging from 2,212 to 11,236 m2 (>400% change), driven primarily by sensitivity to initial conditions, with nonnegligible contributions from boundary condition variability during the initial 30 s of simulation. Our results highlight the need for ensemble simulations, especially when considering fire behavior in marginal burning conditions.

a small parameter. The terms in this series are computed using the symbolic manipulating computer program MACSYMA. 4 refs., 1 fig.

The distribution of fuels is recognised as a key driver of wildland fire behaviour. However, our understanding of how fuel density heterogeneity affects fire behaviour is limited because of the challenges associated with experiments that isolate fuel heterogeneity from other factors. Advances in fire behaviour modelling and computational resources provide a means to explore fire behaviour responses to fuel heterogeneity. Using an ensemble approach to simulate fire behaviour in a coupled fire–atmosphere model, we systematically tested how fuel density fidelity and heterogeneity shape effective wind characteristics that ultimately affect fire behaviour. Results showed that with increased fuel density fidelity and heterogeneity, fire spread and area burned decreased owing to a combination of fuel discontinuities and increased fine-scale turbulent wind structures that blocked forward fire spread. However, at large characteristic length scales of spatial fuel density, the fire spread and area burned increased because local fuel discontinuity decreased, and wind entrainment into the forest canopy maintained near-surface wind speeds that drove forward fire spread. These results demonstrate the importance of incorporating high-resolution fuel fidelity and heterogeneity information to capture effective wind conditions that improve fire behaviour forecasts.

utility of ESMs for predictive understanding and application to societally important problems. This workshop considered (a) dynamic processes that significantly affect terrestrial ecosystems and the coupled Earth system and (b) the data constraints and modeling challenges important for future progress. There were three dominant workshop conclusions: vegetation changes, including disturbance- driven changes, are affecting climate and natural resources; these impacts are expected to increase in the future; and, yet, current models have insufficient data and process representations to adequately predict these changes.