This data release provides inputs needed to run the LANDIS-II landscape change model, NECN and Base Fire extensions for the Greater Yellowstone Ecosystem (GYE), USA, and simulation results that underlie figures and analysis in the accompanying publication. We ran LANDIS-II simulations for 112 years, from 1988-2100, using interpolated weather station data for 1988-2015 and downscaled output from 5 general circulation models (GCMs) for 2016-2100. We also included a control future scenario with years drawn from interpolated weather station data from 1980-2015. Model inputs include raster maps (250 × 250 m grid cells) of climate regions and tables of monthly temperature and precipitation for each climate region. We provide initial conditions in 1987 as rasters and tables (i.e., species-age cohorts, aboveground biomass, soil carbon and nitrogen in surface litter 3 soil layers, soil percent sand, soil percent clay, soil wilting point, soil field capacity, soil drainage, soil storm flow and base flow fractions, and soil depth), historical fire data for model calibration, climate-inferred lognormal fire size distributions for each simulation year, and LANDIS-II control files including parameters for species and functional groups. Outputs from 10 replicates for each of 5 GCMs and the control scenario are provided as rasters and tables. Tables include spatially-weighted mean annual temperature and precipitation of the GYE for each GCM and the control scenario, summarize annual area burned by scenario, summarize biomass pools, and summarize changes in mean stand age. Rasters include annual simulated fire severity for 2015-2100, simulated total aboveground biomass in 4-year timesteps, aboveground biomass of all species in 4-year timesteps, stand age in 4-year timesteps, maximum and minimum cohort age for three dominant species (Pinus contorta, Picea engelmannii, and Pseudotsuga menziesii) in 4-year timesteps, forest type in 1988 and 2100, net ecosystem exchange in 2040 and 2100, and total ecosystem carbon in 4-year timesteps.

Derived maps of landcover type and crown and stem biomass as model inputs to determine annual evapotranspiration, gross primary production, autotrophic respiration and net primary productivity within the BOREAS SSA-MSA, at a 30 m spatial resolution. Mode

Field measurements, daily meteorological inputs, and previously validated iSnobal simulations were used to run and inform the biogeochemical models Biome-BGC and Biome-BGC MuSo at three aspen stands in the Reynolds Creek Experimental Watershed. iSnobal simulations of snow redistribution were used to modify measured precipitation values to account for the redistribution of snow. Biome-BGC simulations were run under historical conditions (1984-2015) assuming both a uniform and redistributed snow layer. Biome-BGC MuSo simulations were run under historical (1996-2015) and future climate scenarios (2046-2065) and account for the redistribution of snow. Biogeochemical simulation data sets include input files used to run Biome-BGC and Biome-BGC MuSo simulations of aspen at three sites in the Reynolds Creek Experimental Watershed under historical and mid-21st conditions. Input files include .ini files describing site conditions and outputs, .epc files describing ecophysiological parameters, and .met files containing historical and modified climate data used to run simulations under historical and mid-21st century conditions. Variables associated with daily simulation outputs are defined in .ini files. Field measurement datasets include hourly soil moisture measurements, monthly pre-dawn leaf water potential measurements, and leaf area index (LAI) measurements collected at each site between 2012 and 2015. iSnobal simulation data sets include daily iSnobal simulated snow water equivalent (SWE) extracted from a single point in drifts located at each site. Data are displayed based on day of water year. Number of simulation years varies from 2 to 24 years depending on site. Drift factors were calculated based on the ratio of peak SWE extracted from iSnobal simulations and uniform SWE falling across the drift accumulation period. Drift factors were then applied to measured precipitation values based on 0°C day and night time temperature thresholds.

This dataset, derived from the National Forest Carbon Monitoring System (NFCMS), provides estimates of forest carbon stocks and fluxes in the form of aboveground woody biomass (AGB), total live biomass, total ecosystem carbon, aboveground coarse woody debris (CWD), and net ecosystem productivity (NEP) as a function of the number of years since the most recent disturbance (i.e., stand age) for forests of the conterminous U.S. at a 30 m resolution for the benchmark years 1990, 2000, and 2010. The data were derived from an inventory-constrained version of the Carnegie-Ames-Stanford Approach (CASA) carbon cycle process model that accounts for disturbance processes for each combination of forest type, site productivity, and pre-disturbance biomass. Also provided are the core model data inputs including the year of the most recent disturbance according to the North American Forest Dynamics (NAFD) and the Monitoring Trends in Burn Severity (MTBS) data products; the type of disturbance; biomass estimates from the year 2000 according to the National Biomass and Carbon Dataset (NBCD); forest-type group; a site productivity classification; and the number of years since stand-replacing disturbance. The data are useful for a wide range of applications including monitoring and reporting recent dynamics of forest carbon across the conterminous U.S., assessment of recent trends with attribution to disturbance and regrowth drivers, conservation planning, and assessment of climate change mitigation opportunities within the forest sector.

BIOME-BGC is a general ecosystem process model designed to simulate biogeochemical and hydrologic processes across multiple scales. BIOME-BGC is used to estimate daily water and carbon budgets for the BOREAS tower flux sites for 1994.

Various aspects of primary production of a variety of plant species found in natural temperate biomes were compiled from literature and presented for use with process-based ecosystem simulation models or ecosystem studies. Information was selected according to the input parameter needs of the BIOME-BGC process-based simulation model.

This dataset is a synthesis of field plot characterization data, derived above-ground and below-ground combusted carbon, and acquired Fire Weather Index (FWI) System components for burned boreal forest sites across Alaska, USA, the Northwest Territories, and Saskatchewan, Canada from 1983-2016. Unburned plot data are also included. Compiled plot-level characterization data include stand age, disturbance history, tree density, and tree biophysical measurements for calculation of the above-ground (ag) and below-ground (bg) biomass/carbon pools, pre-fire and residual post-fire soil organic layer (SOL) depths and estimates of combustion of tree structural classes. The measured slope and aspect for each site and an assigned moisture class based on topography are also provided. Data from 1019 burned and 152 unburned sites are included. From the estimates of combusted ag and bg carbon pools and SOL losses, the total carbon combusted, the proportion of pre-fire carbon combusted, and the proportion of total carbon combusted were calculated for each plot. FWI System components including moisture and drought codes and indices of fire danger were obtained for each plot from existing data sources based on the plot location, year of burn, and a dynamic start-up date (day of burn, DOB) from the global fire weather database. Data for soil characteristics are included in a separate file.