Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “climate dynamics”

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,063 records · Page 59

Effect of cloud-radiation feedback on the climate of a general circulation model

Experiments are described which show significant changes in the simulated large-scale dynamical circulation of a global model. Fixed clouds acting as zonally asymmetric radiative heat sources increase the generation of eddy available potential energy (EAPE) and the energy's conversion to eddy kinetic energy. Generation of EAPE by net radiative heating increases by 50% (0.11 W/sq m) for the fixed cloud experiment. The increase caused by the stationary component is much larger (approximately 100%), but it is partially compensated by a decrease caused by the transient component. A substantial increase is found in the variances of the planetary-scale stationary waves and the medium-scale waves of 2.7 day period. Although the sea surface temperatures are prescribed identically in both integrations, the changes in evaporation and precipitation are found to be much larger over the oceans than over the land.

Shukla, J.↗

Modeling the small-scale deposition of snow onto structured Arctic sea ice during a MOSAiC storm using snowBedFoam 1.0.

The remoteness and extreme conditions of the Arctic make it a very difficult environment to investigate. In these polar regions covered by sea ice, the wind is relatively strong due to the absence of obstructions and redistributes a large part of the deposited snow mass, which complicates estimates for precipitation hardly distinguishable from blowing or drifting snow. Moreover, the snow mass balance in the sea ice system is still poorly understood, notably due to the complex structure of its surface. Quantitatively assessing the snow distribution on sea ice and its connection to the sea ice surface features is an important step to remove the snow mass balance uncertainties (i.e., snow transport contribution) in the Arctic environment. In this work we introduce snowBedFoam 1.0., a physics-based snow transport model implemented in the open-source fluid dynamics software OpenFOAM. We combine the numerical simulations with terrestrial laser scan observations of surface dynamics to simulate snow deposition in a MOSAiC (Multidisciplinary Drifting Observatory for the Study of Arctic Climate) sea ice domain with a complicated structure typical for pressure ridges. The results demonstrate that a large fraction of snow accumulates in their vicinity, which compares favorably against scanner measurements. However, the approximations imposed by the numerical framework, together with potential measurement errors (precipitation), give rise to quantitative inaccuracies, which should be addressed in future work. The modeling of snow distribution on sea ice should help to better constrain precipitation estimates and more generally assess and predict snow and ice dynamics in the Arctic.

54 ENVIRONMENTAL SCIENCES↗

Discoveries about Tropical Tropospheric Ozone from Satellite and SHADOZ (Southern Hemisphere Additional Ozonesondes) and a Future Perspective on NASA's Ozone Sensors

We have been producing near-real tropical tropospheric ozone ('TTO') data from TOMS since 1997 with Prof. Hudson and students at the University of Maryland. Maps for 1996-2000 for the operational Earth-Probe instrument reside at: . We also have archived 'TTO' data from the Nimbus 7/TOMS satellite (1979-1992). The tropics is a region strongly influenced by natural variability and anthropogenic activity and the satellite data have been used to track biomass burning pollution and to detect interannual variability and climate signals in ozone. We look forward to future ozone sensors from NASA; four will be launched in 2004 as part of the EOS AURA Mission. The satellite view of chemical-dynamical interactions in tropospheric ozone is not adequate to capture vertical variability. Thus, in 1998, NASA's Goddard Space Flight Center, NOAA's Climate Monitoring and Diagnostics Laboratory (CMDL) and a team of international sponsors established the SHADOZ (Southern Hemisphere ADditional OZonesondes) project to address the gap in tropical ozone soundings. SHADOZ augments launches at selected sites and provides a public archive of ozonesonde data from twelve tropical and subtropical stations at http://croc.nsfc.nasa.gov/shadoz. The stations are: Ascension Island; Nairobi, Kenya; Irene, South Africa; R,union Island; Watukosek, Java; Fiji; Tahiti; American Samoa; San Cristobal, Galapagos; Natal, Brazil, Malindi, Kenya; Paramaribo, Surinam. From the first 3-4 years of data (presently greater than 1700 sondes), the following features emerge: (a) highly variable tropospheric ozone; (b) a zonal wave-one pattern in tropospheric column ozone; (c) tropospheric ozone variability over the Indian and Pacific Ocean displays strong convective signatures.

Thompson, Anne↗

Strategic roadmap to assess forest vulnerability under air pollution and climate change

Abstract Although it is an integral part of global change, most of the research addressing the effects of climate change on forests have overlooked the role of environmental pollution. Similarly, most studies investigating the effects of air pollutants on forests have generally neglected the impacts of climate change. We review the current knowledge on combined air pollution and climate change effects on global forest ecosystems and identify several key research priorities as a roadmap for the future. Specifically, we recommend (1) the establishment of much denser array of monitoring sites, particularly in the South Hemisphere; (2) further integration of ground and satellite monitoring; (3) generation of flux‐based standards and critical levels taking into account the sensitivity of dominant forest tree species; (4) long‐term monitoring of N, S, P cycles and base cations deposition together at global scale; (5) intensification of experimental studies, addressing the combined effects of different abiotic factors on forests by assuring a better representation of taxonomic and functional diversity across the ~73,000 tree species on Earth; (6) more experimental focus on phenomics and genomics; (7) improved knowledge on key processes regulating the dynamics of radionuclides in forest systems; and (8) development of models integrating air pollution and climate change data from long‐term monitoring programs.

54 ENVIRONMENTAL SCIENCES↗

Asymmetries in cloud microphysical properties ascribed to sea ice leads via water vapour transport in the central Arctic

To investigate the influence of sea ice openings like leads on wintertime Arctic clouds, the air mass transport is exploited as a heat and humidity feeding mechanism which can modify Arctic cloud properties. Cloud microphysical properties in the central Arctic are analysed as a function of sea ice conditions during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in 2019–2020. The Cloudnet classification algorithm is used to characterize the clouds based on remote sensing observations and the atmospheric thermodynamic state from the observatory on board the research vessel (RV) Polarstern. To link the sea ice conditions around the observational site with the cloud observations, the water vapour transport (WVT) being conveyed towards RV Polarstern has been utilized as a mechanism to associate upwind sea ice conditions with the measured cloud properties. This novel methodology is used to classify the observed clouds as coupled or decoupled to the WVT based on the location of the maximum vertical gradient of WVT height relative to the cloud-driven mixing layer. Only a conical sub-sector of sea ice concentration (SIC) and the lead fraction (LF) centred on the RV Polarstern location and extending up to 50 km in radius and with an azimuth angle governed by the time-dependent wind direction measured at the maximum WVT is related to the observed clouds. We found significant asymmetries for cases when the clouds are coupled or decoupled to the WVT and selected by LF regimes. Liquid water path of low-level clouds is found to increase as a function of LF, while the ice water path does so only for deep precipitating systems. Clouds coupled to WVT are found to generally have a lower cloud base and larger thickness than decoupled clouds. Thermodynamically, for coupled cases the cloud-top temperature is warmer and accompanied by a temperature inversion at the cloud top, whereas the decoupled cases are found to be closely compliant with the moist adiabatic temperature lapse rate. The ice water fraction within the cloud layer has been found to present a noticeable asymmetry when comparing coupled versus decoupled cases. This novel approach of coupling sea ice to cloud properties via the WVT mechanism unfolds a new tool to study Arctic surface–atmosphere processes. With this formulation, long-term observations can be analysed to enforce the statistical significance of the asymmetries. Furthermore, our results serve as an opportunity to better understand the dynamic linkage between clouds and sea ice and to evaluate its representation in numerical climate models for the Arctic system.

54 ENVIRONMENTAL SCIENCES↗

Decadal Spiciness Variability in the Subtropical-Tropical Pacific in the CESM2 Large Ensemble

Tropical Pacific decadal variations impact weather and climate around the world and are also connected to variations in the global warming trend. The mechanisms driving these long-term modulations, particularly the role of subsurface ocean dynamics, are still debated. Here, we investigate the dynamics of spiciness (density-compensated temperature and salinity) anomalies in the tropical and subtropical Pacific, which are hypothesized as a possible driving mechanism of decadal climate variability. Based on the analysis of 100 realizations from the Community Earth System Model Version 2 Large Ensemble (CESM2-LE), we demonstrate a coupling between the subtropics and the equatorial Pacific by propagating spiciness anomalies at decadal time scales. The CESM2-LE simulates spiciness variability along a subduction path from the subtropics to the equator with frequency spectra that show the highest power at low frequencies and a power decay proportional to a −4 slope for frequencies greater than 0.01 cycles per months, corresponding to periods smaller than ∼8.5 years. Signals that originate in the Southern Hemisphere (SH) dominate and arrive with a larger magnitude at the equator compared to spiciness anomalies from the Northern Hemisphere (NH). Spiciness anomalies from the SH have shorter propagation times and are strengthened along their pathway as stochastic wind stress curl forcing generates anomalous baroclinic ocean pressure gradients. These pressure gradients generate spiciness anomalies via anomalous advection across climatological spiciness gradients in the SH. We conclude that the observed spiciness variance at decadal time scales is consistent with a forcing by stochastic wind variations that are low-pass filtered by ocean dynamics.

54 ENVIRONMENTAL SCIENCES↗

Consequences of Altered Root Nutrient Uptake for Soil Carbon Stabilization (Final Report)

Objectives: The primary objectives of this project are to improve our understanding of tropical forest belowground processes, and to increase representation of this biome in ecosystem-scale and global C cycle models. Belowground dynamics present a major source of uncertainty inhibiting our ability to predict C cycle responses to climate change. Representation of tropical forests in these models is particularly lacking, even though these ecosystems hold >25% of terrestrial C stocks. While the majority of humid tropical forests exist on relatively infertile soils (i.e. poor in rock-derived nutrients), enormous gradients in soil nutrient availability exist at landscape-scales due to shifts in geology and soil order. Soil fertility is likely to greatly influence how tropical soil C storage will respond to the declines in precipitation predicted for the tropics. In particular, root dynamics vary across soil fertility gradients, with roots representing the major input of C to subsurface soils. Root characteristics related to soil fertility include biomass, turnover, C exudates, tissue chemistry, and nutrient uptake rates, with each of these likely sensitive to changes in moisture. The proposed project will fully integrate field research with an existing ecosystem model to assess potential effects of drying on belowground tropical C dynamics across a range of soil fertilities. Project Description:. This project will identify key linkages among soil fertility, root nutrient uptake, root dynamics, and soil C storage. The hypothesized framework for understanding these linkages is: Tropical rainforest roots in relatively fertile soils have lesser root biomass, faster turnover, fewer exudates, and improved root tissue quality relative to infertile sites, all resulting from increased root nutrient uptake. Ultimately, these root characteristics in fertile soils promote proportionally greater long-term soil C storage in organo-mineral associations. Specifically, lower root exudates and improved root tissue quality in fertile soils reduce microbial respiration of extant soil C (i.e. “priming) and increase microbial C use efficiency, leading to sorption of protein-rich microbial and plant organic matter on soil mineral surfaces. However, the larger pools of mineral-associated C in fertile tropical soils are more vulnerable to decreased rainfall than smaller C pools in infertile soils, in part because of the lesser root biomass, faster turnover, and softer, more labile root tissues. Also, the shallower lateral rooting structures typical in fertile soils experience greater desiccation and death than deeper roots in infertile sites, with little short-term adaptation to drying.

54 ENVIRONMENTAL SCIENCES↗

Towards an improved representation of carbonaceous aerosols over the Indian monsoon region in a regional climate model: RegCM

Mitigation of carbonaceous aerosol emissions is expected to provide climate and health co-benefits. The accurate representation of carbonaceous aerosols in climate models is critical for reducing uncertainties in their climate feedback. In this regard, emission fluxes and aerosol life cycle processes are the two primary sources of uncertainties. Here, we demonstrate that the incorporation of a dynamic ageing scheme and emission estimates that are updated for the local sources improves the representation of carbonaceous aerosols over the Indian monsoon region in a regional climate model, RegCM, compared with its default configuration. The respective mean black carbon (BC) and organic carbon (OC) surface concentrations in 2010 are estimated to be 4.25 and 10.35 µg m -3 over the Indo-Gangetic Plain (IGP) in the augmented model. The BC column burden over the polluted IGP is found to be 2.47 mg m -2 , 69.95 % higher than in the default model configuration and much closer to available observations. The anthropogenic aerosol optical depth (AOD) increases by more than 19 % over the IGP due to the model enhancement, also leading to a better agreement with observed AOD. The respective top-of-the-atmosphere, surface, and atmospheric anthropogenic aerosol short-wave radiative forcing are estimated at -0.3, -9.3, and 9.0 W m -2 over the IGP and -0.89, -5.33, and 4.44 W m -2 over Peninsular India (PI). Our results suggest that the combined effect of two modifications leads to maximum improvements in the model performance in regions where emissions play a dominant role.

54 ENVIRONMENTAL SCIENCES↗

Impact of the numerical solution approach of a plant hydrodynamic model (v0.1) on vegetation dynamics

Abstract. Numerous plant hydrodynamic models have started to be implemented in vegetation dynamics models, reflecting the central role of plant hydraulic traits in driving water, energy, and carbon cycles, as well as plant adaptation to climate change. Different numerical approximations of the governing equations of the hydrodynamic models have been documented, but the numerical accuracy of these models and its subsequent effects on the simulated vegetation function and dynamics have rarely been evaluated. Using different numerical solution methods (including implicit and explicit approaches) and vertical discrete grid resolutions, we evaluated the numerical performance of a plant hydrodynamic module in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES-HYDRO version 0.1) based on single-point and global simulations. Our simulation results showed that when near-surface vertical grid spacing is coarsened (grid size >10 cm), the model significantly overestimates aboveground biomass (AGB) in most of the temperate forest locations and underestimates AGB in the boreal forest locations, as compared to a simulation with finer vertical grid spacing. Grid coarsening has a small effect on AGB in the tropical zones of Asia and South America. In particular, coarse surface grid resolution should not be used when there are large and prolonged water content differences among soil layers at depths due to long dry-season duration and/or well-drained soil or when soil evaporation is a dominant fraction of evapotranspiration. Similarly, coarse surface grid resolution should not be used when there is lithologic discontinuity along the soil depth. This information is useful for uncertainty quantification, sensitivity analysis, or the training of surrogate models to design the simulations when computational cost limits the use of ensemble simulations.

54 ENVIRONMENTAL SCIENCES↗

Online real-time learning of dynamical systems from noisy streaming data

Abstract Recent advancements in sensing and communication facilitate obtaining high-frequency real-time data from various physical systems like power networks, climate systems, biological networks, etc. However, since the data are recorded by physical sensors, it is natural that the obtained data is corrupted by measurement noise. In this paper, we present a novel algorithm for online real-time learning of dynamical systems from noisy time-series data, which employs the Robust Koopman operator framework to mitigate the effect of measurement noise. The proposed algorithm has three main advantages: (a) it allows for online real-time monitoring of a dynamical system; (b) it obtains a linear representation of the underlying dynamical system, thus enabling the user to use linear systems theory for analysis and control of the system; (c) it is computationally fast and less intensive than the popular extended dynamic mode decomposition (EDMD) algorithm. We illustrate the efficiency of the proposed algorithm by applying it to identify the Van der Pol oscillator, the chaotic attractor of the Henon map, the IEEE 68 bus system, and a ring network of Van der Pol oscillators.

97 MATHEMATICS AND COMPUTING↗

On the Settling Depth of Meltwater Escaping From Beneath Antarctic Ice Shelves

Antarctic glacial meltwater is thought to play an important role in determining large-scale Southern Ocean climate trends, yet recent modeling efforts have proceeded without a good understanding of how its vertical distribution in the water column is set. To rectify this, here we conduct new large-eddy simulations of the ascent of a buoyant meltwater plume after its escape from beneath an Antarctic ice shelf. We find that the meltwater’s settling depth is primarily a function of the buoyancy forcing per unit width of the source and the ambient stratification, consistent with the classical theory of turbulent buoyant plumes and in contrast to previous work that suggested an important role for centrifugal instability. Our results further highlight the significant role played by localized variability in stratification; this helps explain observed interannual variability in the vertical meltwater distribution near Pine Island Glacier. Because of the vast heterogeneity in mass loss rates and ambient conditions at different Antarctic ice shelves, a dynamic parameterization of meltwater settling depth may be crucial for accurately simulating high-latitude climate in a warming world; we discuss how this may be developed following this work, and where the remaining challenges lie.

Antarctic glacial meltwater↗

Air‐Ice‐Ocean Coupling During a Strong Mid‐Winter Cyclone: Observing Coupled Dynamic Interactions Across Scales

Abstract Arctic cyclones are key drivers of sea ice and ocean variability. During the 2019–2020 Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition, joint observations of the coupled air‐ice‐ocean system were collected at multiple spatial scales. Here, we present observations of a strong mid‐winter cyclone that impacted the MOSAiC site as it drifted in the central Arctic pack ice. The sea ice dynamical response showed spatial structure at the scale of the evolving and translating cyclonic wind field. Internal ice stress and ocean stress play significant roles, resulting in timing offsets between the atmospheric forcing and the ice response and post‐cyclone inertial ringing in the ice and ocean. Ice motion in response to the wind field then forces the upper ocean currents through frictional drag. The strongest impacts to the sea ice and ocean from the passing cyclone occur as a result of the surface impacts of a strong atmospheric low‐level jet (LLJ) behind the trailing cold front and changing wind directions between the warm‐sector LLJ and post cold‐frontal LLJ. Impacts of the cyclone are prolonged through the coupled ice‐ocean inertial response. Local impacts of the approximately 120 km wide LLJ occur over a 12 hr period or less and at scales of a kilometer to a few tens of kilometers, meaning that these impacts occur at combined smaller spatial scales and faster time scales than most satellite observations and coupled Earth system models can resolve.

Meteorology & Atmospheric Sciences↗

MITgcm-AD v2: Open source tangent linear and adjoint modeling framework for the oceans and atmosphere enabled by the Automatic Differentiation tool Tapenade

The Massachusetts Institute of Technology General Circulation Model (MITgcm) is widely used by the climate science community to simulate planetary atmosphere and ocean circulations. A defining feature of the MITgcm is that it has been developed to be compatible with an algorithmic differentiation (AD) tool, TAF, enabling the generation of tangent-linear and adjoint models. These provide gradient information which enables dynamics-based sensitivity and attribution studies, state and parameter estimation, and rigorous uncertainty quantification. Importantly, gradient information is essential for computing comprehensive sensitivities and performing efficient large-scale data assimilation, ensuring that observations collected from satellites and in-situ measuring instruments can be effectively used to optimize a large uncertain control space. As a result, the MITgcm forms the dynamical core of a key data assimilation product employed by the physical oceanography research community: Estimating the Circulation and Climate of the Ocean (ECCO) state estimate. Although MITgcm and ECCO are used extensively within the research community, the AD tool TAF is proprietary and hence inaccessible to a large proportion of these users. The new version 2 (MITgcm-AD v2) framework introduced here is based on the source-to-source AD tool Tapenade, which has recently been open-sourced. Another feature of Tapenade is that it stores required variables by default (instead of recomputing them) which simplifies the implementation of efficient, AD-compatible code. The framework has been integrated with the MITgcm model’s main branch and is now freely available.

Adjoints↗

Some Applications of SiFSAP for Wildfires, Weather and Atmospheric Dynamics Studies

The Single Field of View (SFOV) Sounder Atmospheric Product (SiFSAP) derived from Cross-track Infrared Sounder (CrIS) on SNPP and JPSS have a spatial resolution ( ~14.5 km) better than most global weather and climate models. Most recent improvement in the quality of these products enables us to use these high-resolution observation-based sounding data for air quality, weather and atmospheric dynamics studies. A few cases studies through a synergistic use of SiFSAP from S-NPP and J-1 and the reanalysis data from ERA5 and MERRA-2 will be presented, which include a process-oriented analysis of (1) a stratospheric intrusion (SI) event occurred on June 11-13, 2017 in the southwestern US, (2) one extreme cold air break (CAO) episode occurred on 27-31 January 2019 across much of the US Midwest and its impact from stratospheric downward transport, and (3) dynamic transport of carbon monoxide (CO) associated with Australia’s unprecedented fire from the end of 2019 to early 2020. These examples will demonstrate the value of SiFSAP for weather, atmospheric dynamics, air quality, wildfires studies and for model evaluations.

Xiaozhen (Shawn) Xiong↗

Active layer depth and soil properties impact specific leaf area variation and ecosystem productivity in a boreal forest

Specific leaf area (SLA, leaf area per unit dry mass) is a key canopy structural characteristic, a measure of photosynthetic capacity, and an important input into many terrestrial process models. Although many studies have examined SLA variation, relatively few data exist from high latitude, climate-sensitive permafrost regions. We measured SLA and soil and topographic properties across a boreal forest permafrost transition, in which dominant tree species changed as permafrost deepened from 54 to >150 cm over 75 m hillslope transects in Caribou-Poker Creeks Research Watershed, Alaska. We characterized both linear and threshold relationships between topographic and edaphic variables and SLA and developed a conceptual model of these relationships. We found that the depth of the soil active layer above permafrost was significantly and positively correlated with SLA for both coniferous and deciduous boreal tree species. Intraspecific SLA variation was associated with a fivefold increase in net primary production, suggesting that changes in active layer depth due to permafrost thaw could strongly influence ecosystem productivity. While this is an exploratory study to begin understanding SLA variation in a non-contiguous permafrost system, our results indicate the need for more extensive evaluation across larger spatial domains. These empirical relationships and associated uncertainty can be incorporated into ecosystem models that use dynamic traits, improving our ability to predict ecosystem-level carbon cycling responses to ongoing climate change.

59 BASIC BIOLOGICAL SCIENCES↗

Physically-Derived Dynamical Cores in Atmospheric General Circulation Models

The algorithm chosen to represent the advection in atmospheric models is often used as the primary attribute to classify the model. Meteorological models are generally classified as spectral or grid point, with the term grid point implying discretization using finite differences. These traditional approaches have a number of shortcomings that render them non-physical. That is, they provide approximate solutions to the conservation equations that do not obey the fundamental laws of physics. The most commonly discussed shortcomings are overshoots and undershoots which manifest themselves most overtly in the constituent continuity equation. For this reason many climate models have special algorithms to model water vapor advection. This talk focuses on the development of an atmospheric general circulation model which uses a consistent physically-based advection algorithm in all aspects of the model formulation. The shallow-water model is generalized to three dimensions and combined with the physics parameterizations of NCAR's Community Climate Model. The scientific motivation for the development is to increase the integrity of the underlying fluid dynamics so that the physics terms can be more effectively isolated, examined, and improved. The expected benefits of the new model are discussed and results from the initial integrations will be presented.

Rood, Richard B.↗

Physically-Derived Dynamical Cores in Atmospheric General Circulation Models

The algorithm chosen to represent the advection in atmospheric models is often used as the primary attribute to classify the model. Meteorological models are generally classified as spectral or grid point, with the term grid point implying discretization using finite differences. These traditional approaches have a number of shortcomings that render them non-physical. That is, they provide approximate solutions to the conservation equations that do not obey the fundamental laws of physics. The most commonly discussed shortcomings are overshoots and undershoots which manifest themselves most overtly in the constituent continuity equation. For this reason many climate models have special algorithms to model water vapor advection. This talk focuses on the development of an atmospheric general circulation model which uses a consistent physically-based advection algorithm in all aspects of the model formulation. The shallow-water model of Lin and Rood (QJRMS, 1997) is generalized to three dimensions and combined with the physics parameterizations of NCAR's Community Climate Model. The scientific motivation for the development is to increase the integrity of the underlying fluid dynamics so that the physics terms can be more effectively isolated, examined, and improved. The expected benefits of the new model are discussed and results from the initial integrations will be presented.

Rood, Richard B.↗

Importance of Tree- and Species-Level Interactions With Wildfire, Climate, and Soils in Interior Alaska: Implications for Forest Change Under A Warming Climate

The boreal zone of Alaska is dominated by interactions between disturbances, vegetation, and soils. These interactions are likely to change in the future through increasing permafrost thaw, more frequent and intense wildfires, and vegetation change from drought and competition. We utilize an individual tree-based vegetation model, the University of Virginia Forest Model Enhanced (UVAFME), to estimate current and future forest conditions across sites within interior Alaska. We updated UVAFME for application within interior Alaska, including improved simulation of permafrost dynamics, litter decay, nutrient dynamics, fire mortality, and post-fire regrowth. Following these updates, UVAFME output on species-specific biomass and stem density was comparable to inventory measurements at various forest types within interior Alaska. We then simulated forest response to climate change at specific inventory locations and across the Tanana Valley River Basin on a 2 × 2 km2 grid. We derived projected temperature and precipitation from a five-model average taken from the CMIP5 archive under the RCP 4.5 and 8.5 scenarios. Results suggest that climate change and the concomitant impacts on wildfire and permafrost dynamics will result in overall decreases in biomass (particularly for spruce (Picea spp.)) within the interior Tanana Valley, despite increases in quaking aspen (Populus tremuloides) biomass, and a resulting shift towards higher deciduous fraction. Simulation results also predict increases in biomass at cold, wet locations and at high elevations, and decreases in biomass in dry locations, under both moderate (RCP 4.5) and extreme (RCP 8.5) climate change scenarios. These simulations demonstrate that a highly detailed, species interactive model can be used across a large region within Alaska to investigate interactions between vegetation, climate, wildfire, and permafrost. The vegetation changes predicted here have the capacity to feed back to broader scale climate-forest interactions in the North American boreal forest, a region which contributes significantly to the global carbon and energy budgets.

forest modeling↗