Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “terrestrial models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,045 records · Page 58

Evaluating the Capacity of Global CO2 Flux and Atmospheric Transport Models to Incorporate New Satellite Observations

As we enter the new era of satellite remote sensing for CO2 and other carbon cyclerelated quantities, advanced modeling and analysis capabilities are required to fully capitalize on the new observations. Model estimates of CO2 surface flux and atmospheric transport are required for initial constraints on inverse analyses, to connect atmospheric observations to the location of surface sources and sinks, and ultimately for future projections of carbon-climate interactions. For application to current, planned, and future remotely sensed CO2 data, it is desirable that these models are accurate and unbiased at time scales from less than daily to multi-annual and at spatial scales from several kilometers or finer to global. Here we focus on simulated CO2 fluxes from terrestrial vegetation and atmospheric transport mutually constrained by analyzed meteorological fields from the Goddard Modeling and Assimilation Office for the period 1998 through 2006. Use of assimilated meteorological data enables direct model comparison to observations across a wide range of scales of variability. The biospheric fluxes are produced by the CASA model at lxi degrees on a monthly mean basis, modulated hourly with analyzed temperature and sunlight. Both physiological and biomass burning fluxes are derived using satellite observations of vegetation, burned area (as in GFED-2), and analyzed meteorology. For the purposes of comparison to CO2 data, fossil fuel and ocean fluxes are also included in the transport simulations. In this presentation we evaluate the model's ability to simulate CO2 flux and mixing ratio variability in comparison to in situ observations at sites in Northern mid latitudes and the continental tropics. The influence of key process representations is inferred. We find that the model can resolve much of the hourly to synoptic variability in the observations, although there are limits imposed by vertical resolution of boundary layer processes. The seasonal cycle and its interannual variations generally respond adequately, but discrepancies in the tropics suggest the need for a refinement of the soil moisture dependence of the respiration flux in CASA. Examples and inferences for interpretation of satellite data will be discussed. In general, the fidelity of these simulations leads us to anticipate incorporation of real-time, highly resolved remote sensing and other observations into quantitative analyses that will reduce uncertainty in the terrestrial CO2 sink and revolutionize our understanding of the key processes controlling atmospheric CO2 and its evolution with time.

Kawa, S. R.↗

Modeling exports of dissolved organic carbon from landscapes: a review of challenges and opportunities

Inland waters receive large quantities of dissolved organic carbon (DOC) from soils and act as conduits for the lateral transport of this terrestrially derived carbon, ultimately storing, mineralizing, or delivering it to oceans. The lateral DOC flux plays a crucial role in the global carbon cycle, and numerous models have been developed to estimate the DOC export from different landscapes. We reviewed 34 published models and compared their characteristics to identify challenges in model applications and opportunities for future model development. We classified these models into three types: indicator-driven, hydrology-forced, and process-based DOC export simulation models. They differ mainly in their environmental inputs, simulation approaches for soil DOC production, leaching from soils to inland waters, and transit through inland waters. It is essential to consider landscape characteristics, climate conditions, available data, and research questions when selecting the most appropriate model. Given the substantial assumptions associated with these models, sufficient measurements are required to benchmark estimates. Accurate accounting of terrestrially derived DOC export to oceans requires incorporating the DOC produced in aquatic ecosystems and deposited with rainwater; otherwise, global export estimates may be overestimated by 40.7%. Additionally, improving the representation of mineralization and burial processes in inland waters allows for more accurate accounting of carbon sequestration through land ecosystems. When all the inland water processes are ignored or assuming DOC leaching is equivalent to DOC export, the loss of soil carbon through this lateral flux could be underestimated by 43.9%.

54 ENVIRONMENTAL SCIENCES↗

Testing Cosmic Inflation

The Cosmic Microwave Background (CMB) has provided a wealth of information about the history and physics of the early Universe. Much progress has been made on uncovering the emerging Standard Model of Cosmology by such experiments as COBE and WMAP, and ESA's Planck Surveyor will likely increase our knowledge even more. Despite the success of this model, mysteries remain. Currently understood physics does not offer a compelling explanation for the homogeneity, flatness, and the origin of structure in the Universe. Cosmic Inflation, a brief epoch of exponential expansion, has been posted to explain these observations. If inflation is a reality, it is expected to produce a background spectrum of gravitational waves that will leave a small polarized imprint on the CMB. Discovery of this signal would give the first direct evidence for inflation and provide a window into physics at scales beyond those accessible to terrestrial particle accelerators. I will briefly review aspects of the Standard Model of Cosmology and discuss our current efforts to design and deploy experiments to measure the polarization of the CMB with the precision required to test inflation.

Chuss, David↗

Regolith-structure modeling of lunar facilities

The mechanical properties of lunar regolith are significantly different from terrestrial soils; at low stress levels, the ultimate strength envelope is highly nonlinear, and the material is essentially dilatant. The effects of platen friction, membrane confinement, and material self weight in triaxial tests have been shown to be of great significance in low stress level tests. An inverse identification technique has been devised to account for these effects where the load-displacement response of the entire specimen, including the boundary conditions and the self-weight, was analyzed by nonlinear finite element techniques so as to obtain true material parameters.

Perkins, Steven W.↗

Modeling Phenological Controls on Carbon Dynamics in Dryland Sagebrush Ecosystems

Dryland ecosystems play an important role in determining how precipitation anomalies affect terrestrial carbonfluxes at regional to global scales. Thus, to understand how climate change may affect the global carbon cycle,we must also be able to understand and model its effects on dryland vegetation. Dynamic Global VegetationModels (DGVMs) are an important tool for modeling ecosystem dynamics, but they often struggle to reproduceseasonal patterns of plant productivity. Because the phenological niche of many plant species is linked to bothtotal productivity and competitive interactions with other plants, errors in how process-based models representphenology hinder our ability to predict climate change impacts. This may be particularly problematic in drylandecosystems where many species have developed a complex phenology in response to seasonal variability in bothmoisture and temperature. Here, we examine how uncertainty in key parameters as well as the structure ofexisting phenology routines affect the ability of a DGVM to match seasonal patterns of leaf area index (LAI) andgross primary productivity (GPP) across a temperature and precipitation gradient. First, we optimized modelparameters using a combination of site-level eddy covariance data and remotely-sensed LAI data. Second, wemodified the model to include a semi-deciduous phenology type and added flexibility to the representation ofgrass phenology. While optimizing parameters reduced model bias, the largest gains in model performance wereassociated with the development of our new representation of phenology. This modified model was able to bettercapture seasonal patterns of both leaf area index (R2=0.75) and gross primary productivity (R2=0.84), thoughits ability to estimate total annual GPP depended on using eddy covariance data for optimization. The new modelalso resulted in a more realistic outcome of modeled competition between grass and shrubs. These findingsdemonstrate the importance of improving how DGVMs represent phenology in order to accurately forecastclimate change impacts in dryland ecosystems.

Ecosystem model Phenology Parameter optimization E↗

Validation of Remotely Sensed and Modeled Soil Moisture at Forested and Unforested Sites

Soil moisture is an important driver for forest ecosystems, influencing fire occurrence and extent, insect and pathogen impacts, and tree growth, which creates a need for regular, globally extensive soil moisture information that only satellite-based sensors or models can achieve. However, the reliability of soil moisture measurements in forests is not well understood due to a lack of suitable validation sites (especially relative to unforested ecosystems) and interference caused by high vegetation water content on remotely sensed measurements; although recent studies have started to address this gap [1], [2], [3], [4]. Here we validate the performance of multiyear remotely sensed (SMAP/Sentinel), remotely sensed data assimilation modelled (SMAP-L4), and modelled (NLDAS) surface and root zone (0-1 m) soil moisture datasets with data from in-situ sensors at 39 National Ecological Observatory Network (NEON) sites throughout the contiguous US. Due to differences in spatial resolution, NEON soil moisture (~0.2 km measurement zone) correlations were expected to be stronger with the SMAP/Sentinel product (3 km resolution) than with coarser resolution SMAP-L4 (9 km resolution) or NLDAS products (13 km resolution). However, given the sensitivity of satellite measurements to vegetation water content we expected a deterioration in the correlations based on remotely sensed measurements (SMAP/Sentinel and SMAP-L4) as aboveground biomass increased, whereas the model-based data (NLDAS) was expected to be largely insensitive to vegetation type. We recognize that the SMAP/Sentinel product was developed for unforested regions, therefore our application is outside its primary use case. Soil moisture is measured at up to 8 depths in five soil plots spaced up to 40 m apart at each NEON terrestrial site. Correlation parameters were calculated for the three remotely sensed and modelled data products relative to in-situ measurement following Entekhabi et al. [5]. The datasets comprised 94 (SMAP-L4), 28 (SMAP/Sentinel), and 106 (NLDAS) sites-years for surface soils and 13 (SMAP-L4) and 14 (NLDAS) site-years for the root zone. At unforested sites, the performance of the three remotely sensed and modelled data products was similar for surface soils (Table 1). For example, unbiased RMSD (ubRMSD), which SMAP uses as its primary performance metric [6], ranged from 0.05 to 0.06 m3 m-3 (Table 1), indicating the ability of all three products to track changes in soil moisture over time. The performance of the three products deteriorated at forested sites, however, while the difference in performance was modest for SMAP-L4 and NLDAS, the deterioration in SMAP/Sentinel performance was substantial. For instance, SMAP/Sentinel ubRMSD increased from 0.06 to 0.11 m3 m-3 and absolute mean difference (Abs MD; which includes measurement bias and spatial representativeness errors) increased from 0.06 to 0.16 m3 m-3, indicating both a reduction in ability to track temporal changes and absolute amounts of soil moisture in forest ecosystems. SMAP-L4 and NLDAS had lower unbiased RMSD for root zone (0-1 m) than surface soils at both forested and unforested sites (Tables 1 and 2; SMAP/Sentinel does not produce a root zone measurement). However, in most cases the correlation coefficient (r) was lower for the root zone than surface soils, suggesting the lower unbiased RMSD may be attributed to greater temporal stability of soil moisture in the root zone rather than improved data product performance. Mean difference and absolute mean difference, which encompass measurement bias and spatial representativeness errors, were greater for root zone than surface soils at unforested sites for both data products, but the opposite was generally true at forested sites. As with surface soils, there was relatively little change in the performance of SMAP-L4 and NLDAS between the unforested and forested sites. In summary, all three data products were able to adequately represent soil moisture at unforested sites, at least when aggregating across sites. However, while the performance of all three products deteriorated at forested sites, SMAP-L4 and NLDAS maintained sufficient performance to remain suitable for some use cases (ubRMSD <0.06 m3 m-3 and RMSD <0.13 m3 m-3). In contrast, the relatively poorer performance of the SMAP/Sentinel product at forested sites seems insufficient for most use cases (ubRMSD >0.1 m3 m-3 and RMSD >0.2 m3 m-3). We attribute the large reduction in the performance of the SMAP/Sentinel product in forests to its use of C-band wavelengths, which are particularly sensitive to vegetation interference, and apparently outweighed any gains provided by its higher spatial resolution. A combined SMAP/NISAR soil moisture product may provide improved performance relative to SMAP/Sentinel due to NISAR’s use of L-band wavelengths, which are less sensitive to vegetation (NISAR is scheduled for launch in early 2024).

Edward Ayres↗

Stellar Energetic Particle-driven Production of Biologically Relevant Molecules in Atmospheres of Young Earth-like Exoplanets

The chemistry of N2, CO2, CH4-rich atmospheres of terrestrial-type exoplanets around active G-K stars is a complex problem. The star's energetic inputs in the atmospheres take the form of ionizing radiation in the form of X-ray, Extreme UV, as well as precipitating energetic particles that were accelerated in coronal mass ejection-driven shocks. These inputs lead to the ionization and dissociation of the atmospheric species, which enable complex chemical processes resulting in the creation of molecules of interest for both astrobiology and climate (Airapetian et al., 2016; Hayworth et al. 2022). We developed a set of atmospheric models to simulate the formation of organic molecules that can be considered as atmospheric pre-biosignatures from exoplanets around active stars. They take into account the dissociations and ionizations from all the energetic inputs in the atmosphere as well as their role in the enhanced chemistry in the lower atmosphere. In addition, the evaluation of stellar energetic particle fluxes is done using the latest techniques to accurately highlight the depth at which the chemistry is enabled. These models were used to simulate the chemistry of several young terrestrial-type exoplanets resembling the Hadean Earth to highlight the creation of atmospheric pre-biosignatures and relevant greenhouse gases, in function of different scenarios for its atmospheric evolution. This work shows the importance of space weather in the chemical evolution of atmospheres of astrobiological interest.

Guillaume Gronof↗

Water balance model (WBM) v.1.0.0: a scalable gridded global hydrologic model with water-tracking functionality

Abstract. This paper describes the University of New Hampshire Water Balance Model, WBM, a process-based gridded global hydrologic model that simulates the land surface components of the global water cycle and includes water extraction for use in agriculture and domestic sectors. The WBM was first published in 1989; here, we describe the first fully open-source WBM version (v.1.0.0). Earlier descriptions of WBM methods provide the foundation for the most recent model version that is detailed here. We present an overview of the model functionality, utility, and evaluation of simulated global river discharge and irrigation water use. This new version adds a novel suite of water source tracking modules that enable the analysis of flow-path histories on water supply. A key feature of WBM v.1.0.0 is the ability to identify the partitioning of sources for each stock or flux within the model. Three different categories of tracking are available: (1) primary inputs of water to the surface of the terrestrial hydrologic cycle (liquid precipitation, snowmelt, glacier melt, and unsustainable groundwater); (2) water that has been extracted for human use and returned to the terrestrial hydrologic system; and (3) runoff originating from user-defined spatial land units. Such component tracking provides a more fully transparent model in that users can identify the underlying mechanisms generating the simulated behavior. We find that WBM v.1.0.0 simulates global river discharge and irrigation water withdrawals well, even with default parameter settings, and for the first time, we are able to show how the simulation arrives at these fluxes by using the novel tracking functions.

54 ENVIRONMENTAL SCIENCES↗

Accounting for herbaceous communities in process‐based models will advance our understanding of “grassy” ecosystems

Abstract 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.

59 BASIC BIOLOGICAL SCIENCES↗

A carbonate-silicate aqueous geochemical cycle model for Mars

A model for the carbonate-silicate geochemical cycle of an early, wet Mars is under development. The results of this study will be used to constrain models of the geochemical history of Mars and the likely mineralogy of its present surface. Although Mars today is a cold, dry planet, it may once have been much warmer and wetter. Values of total outgassed CO2 from several to about 10 bars are consistent with present knowledge (Pollack et al. 1987), and this amount of CO2 implies an amount of water outgassed at least equal to an equivalent depth of 500-1000 meters (Carr 1986). Pollack et al. (1987), in addition, estimate that a thick CO2 atmosphere may have existed for an extended period of time, perhaps as long as a billion years. The greenhouse effect of such an atmosphere would permit the presence of liquid water on the surface, most likely in the form of a shallow sea in the lowest regions of the planet, such as the northern plains (Schaefer 1990). The treatment of geochemical cycles as complex kinetic chemical reactions has been undertaken for terrestrial systems in recent years with much success (Lasaga 1980, 1981; Berner et al. 1983; Lasaga et al. 1985). Although the Martian system is vastly less well understood, and hence less well-constrained, it is also a much simpler system, due to the lack of biogenic reactions that make the terrestrial system so complex. It should be possible, therefore, to use the same techniques to model the Martian system as have been used for terrestrial systems, and to produce useful results. A diagram of the carbonate-silicate cycle for Mars (simplified from the terrestrial system) is given.

Schaefer, M. W.↗

Modeling Phenological Controls on Carbon Dynamics in Dryland Sagebrush Ecosystems

Dryland ecosystems play an important role in determining how precipitation anomalies affect terrestrial carbon fluxes at regional to global scales. Thus, to understand how climate change may affect the global carbon cycle, we must also be able to understand and model its effects on dryland vegetation. Dynamic Global Vegetation Models (DGVMs) are an important tool for modeling ecosystem dynamics, but they often struggle to reproduce seasonal patterns of plant productivity. Because the phenological niche of many plant species is linked to both total productivity and competitive interactions with other plants, errors in how process-based models represent phenology hinder our ability to predict climate change impacts. This may be particularly problematic in dryland ecosystems where many species have developed a complex phenology in response to seasonal variability in both moisture and temperature. Here, we examine how uncertainty in key parameters as well as the structure of existing phenology routines affect the ability of a DGVM to match seasonal patterns of leaf area index (LAI) and gross primary productivity (GPP) across a temperature and precipitation gradient. First, we optimized model parameters using a combination of site-level eddy covariance data and remotely-sensed LAI data. Second, we modified the model to include a semi-deciduous phenology type and added flexibility to the representation of grass phenology. While optimizing parameters reduced model bias, the largest gains in model performance were associated with the development of our new representation of phenology. This modified model was able to better capture seasonal patterns of both leaf area index (R(exp 2) = 0.75) and gross primary productivity (R(exp 2) = 0.84), though its ability to estimate total annual GPP depended on using eddy covariance data for optimization. The new model also resulted in a more realistic outcome of modeled competition between grass and shrubs. These findings demonstrate the importance of improving how DGVMs represent phenology in order to accurately forecast climate change impacts in dryland ecosystems.

Ecosystem model↗

Comparison of Vertical Surface Deformation Estimates Derived From Space-Based Gravimetry, Ground-Based GPS, and Model-Based Hydrologic Loading Over Snow-Dominated Watersheds in the United States

Spatiotemporal variability in Earth's terrestrial water storage (TWS) causes changes in surface deformation. The potential for using ground-based Global Positioning System (GPS) vertical displacement observations for estimating TWS is explored through a comparison of vertical displacements derived from space-based gravimetric retrievals, ground-based GPS, and model-based hydrologic estimates. The study presented here focuses on two snow-dominated basins in the Western United States for the years 2003–2016. Seasonal variations are observed in the vertical displacements derived from all three data sets, and the variation is coherent with the changes in hydrologic loading. Good consistency is observed between any two of the three data sets with gravimetric retrievals and hydrologic model estimates providing the highest level of agreement (i.e., all examined stations with correlation coefficient R > 0.70). Vertical displacements derived from gravimetric retrievals and ground-based GPS yielded R > 0.70 for more than 89% of the stations. In addition, it is found that both GPS-derived and space-based, gravimetry-derived vertical displacements clearly reflected the impact of climate variation (i.e., heavy precipitation during 2010–2011 winter followed by prolonged drought). Vertical displacements derived from the hydrologic model highlighted the relatively large precipitation convergence phase during late 2010 to early 2011 at some stations but not the prolonged drought that followed. The results indicate that ground-based GPS observations of vertical displacement have the capability to capture variations in TWS changes, which can be systematically merged in conjunction with Gravity Recovery and Climate Experiment (GRACE) into a land surface model to improve TWS estimates in a follow-up study.

Gaohong Yin↗

Modeling Magnetite Reflectance Spectra Using Hapke Theory and Existing Optical Constants

Magnetite is an accessory mineral found in terrestrial environments, some meteorites, and the lunar surface. The reflectance of magnetite powers is relatively low [1], and this property makes it an analog for other dark Fe- or Ti-bearing components, particularly ilmenite on the lunar surface. The real and imaginary indices of refraction (optical constants) for magnetite are available in the literature [2-3], and online [4]. Here we use these values to calculate the reflectance of particulates and compare these model spectra to reflectance measurements of magnetite available on-line [5].

terrestrial environments↗

Properties of the stratospheric aerosol layer studied with a one-dimensional computer model

Aerosol particle effects are often neglected in theoretical studies of stratospheric phenomena. In reality, the particulate matter normally found above the tropopause may influence the terrestrial radiation balance, catalyze heterogeneous chemical interactions, and serve as a tracer of atmospheric motions. The paper proposes a one-dimensional model of the stratospheric aerosol layer, and it is used to compare aerosol theory with observational data. The model considers gaseous sulfur photochemistry and the physical aerosol processes of nucleation, coagulation, sedimentation, and diffusion. Calculations of the effects on the aerosol layer of stratospheric injections of aluminum oxide particles by Space Shuttle engines and of sulfur dioxide molecules by volcanic activity are performed. The relation between measured aerosol variability and changes in stratospheric air temperatures and vertical transport rates are discussed.

Turco, R. P.↗

Chemistry and transport in a multi-dimensional model

The focus of our research program is the achievement of a quantitative understanding of the spatial distribution and temporal variation of chemical species in the terrestrial middle atmosphere, with emphasis on ozone. Although not directed at assessments of anthropogenic impacts, our activities contribute to a refinement of model descriptions of chemical and dynamical processes that are needed for assessment tasks. A unique feature of our research effort is a close interaction with the chemical kinetics group at JPL and the field measurement groups at JPL. The principal tools used in our investigations are l-D and 2-D photochemical models. The advective component of the meridional transport in these models is approximated by the diabatic circulation, which is derived diagnostically from the net radiative heating rates. We are developing accurate radiative transfer algorithms to find the net heating rates and the diabatic circulation for use in the 2-D model.

Yung, Yuk L.↗

The Dynamical Evolution of the Earth-Moon Progenitors: Motivation and Methodology - 1

The giant impact hypothesis was introduced in the mid-1970s after consideration of results from the Apollo missions. This hypothesis best explains the similarity in elemental proportions in lunar and terrestrial rocks, the depletion of lunar volatiles, the lack of lunar Fe, and the large angular momentum in the Earth-Moon system. Comparison between the radiometric ages of inclusions in the most primitive meteorites and in the oldest lunar rocks and the differentiation age of Earth suggests that the Earth-Moon system formed about100 m.y. after the oldest meteorites. In addition, the age of the famous martian meteorite ALH 84001 and an early Martian solidification time obtained by Lee and Halliday suggest that the inner solar system was fairly clear of large bodies about 10 m.y. after the oldest meteorites formed. Thus, the "standard model" suggests that for several tens of millions of years, the terrestrial planet region had few, if any, lunar-sized bodies, and there were five terrestrial planets: Mercury, Venus, the two progenitors of the Earth-Moon system, and Mars. To simulate the dynamics of the solar system before the hypothesized Moon-forming impact, we are integrating the solar system with the Earth-Moon system replaced by two bodies in heliocentric orbits between Venus and Mars. The total (orbital) angular momentum of the Earth-Moon progenitors is that of the present Earth-Moon system, and their total mass is that of the Earth-Moon System. We are looking at ranges in mass ratio and initial values for eccentricity, inclination, and semimajor axis. We are using the SYMBA integrator to integrate these systems until a collision occurs or a time of 200 m.y. elapses. Results are presented in a companion abstract, (also presented at this meeting).

Lissauer, J. J.↗

Files and scripts to support manuscript Shuman et al 2023 FATES-SPITFIRE ecosystem assembly across tropics

The dataset includes the parameter and domain files, relevant output files, and scripts to generate simulations and perform analysis with Jupyter notebooks that support the manuscript Shuman, JK et al 2023 “Dynamic ecosystem assembly and escaping the “fire-trap” in the tropics: Insights from FATES_15.0.0”. We have adapted the fire-behavior and effects module, SPITFIRE, for use with the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), a size-structured vegetation demographic model. We test how climate, fire regime and fire-tolerance plant traits interact to determine the biogeography of tropical forests and grasslands. We assign different fire-tolerance strategies based on crown, leaf and bark characteristics, which are key observed fire-tolerance traits across woody plants. For these simulations, three types of vegetation compete for resources: a fire-vulnerable tree with thin bark, a vulnerable deep crown and fire-intolerant foliage; a fire-tolerant tree with thick bark, a thin crown and fire-tolerant foliage; and a fire-promoting C4 grass. We explore the model sensitivity to a critical parameter governing fuel moisture, and show that drier fuels promote increased burning, an expansion of area for grass and fire-tolerant trees and a reduction of area for fire-vulnerable trees. This conversion to lower biomass or grass areas with increased fuel drying results in increased fire burned area and its effects, which could fee back to local climate variables. Simulated size-based fire mortality for trees less than 20 cm in diameter and those with fire-vulnerable traits is higher than that for larger and/or fire-tolerant trees, in agreement with observations. Fire-disturbed forests demonstrate reasonable productivity and capture observed patterns of aboveground biomass in areas dominated by natural vegetation for the recent historical period, but have a large bias in less disturbed areas. Though the model predicts a greater extent of burned fraction than observed in areas with grass dominance, the resulting biogeography of fire-tolerant, thick-bark trees and fire-vulnerable, thin-bark trees corresponds to observations across the tropics. In areas with more than 2500 mm of precipitation, simulated fire frequency and burned area are low, with fire intensities below 150 kW m-1, consistent with observed understory fire behavior across the Amazon. Areas drier than this demonstrate fire intensities consistent with those measured in savannas and grasslands, with high values up to 4000 kW m-1. The results support a positive grass-fire feedback across the region, and suggest that forests which have existed without frequent burning may be vulnerable at higher fire intensities, which is of greater concern under intensifying climate and land use pressures. The ability of FATES to capture the connection between fire disturbance and plant fire-tolerance strategies in determining biogeography provides a useful tool for assessing the vulnerability and resilience of these critical carbon storage areas under changing conditions across the tropics.

54 ENVIRONMENTAL SCIENCES↗