Dynamic Vegetation Model Dynamic Organic Soil Terrestrial Ecosystem Model (DVM-DOS-TEM) simulations focused on Eight Mile Lake, Alaska and Imnavait Creek, Alaska [2000-2015]

This set of files store model simulations using the biosphere model Dynamic Vegetation Model Dynamic Organic Soil Terrestrial Ecosystem Model (DVM-DOS-TEM), developed to simulate biophysical and biogeochemical interactions between the soil, vegetation and atmosphere. To improve predictions of net carbon releases from thawing permafrost, we tested the sensitivity of a suite of model parameters. We analyzed the responses of ecosystem carbon balances to permafrost thaw by running site-level simulations at two long-term tundra ecological monitoring sites in Alaska: Eight Mile Lake (EML) and Imnavait Creek watershed (IMN). These sites are characterized by similar tussock tundra vegetation but differing soil drainage conditions and climate, IMN consists of well-drained soils, and EML has historically well-drained soils, however permafrost thaw has altered drainage conditions to wetter soils. Simulations were conducted at a 1km resolution, over a 1,000 km2 area (10x10 km square) centered on two long term ecological research sites in Alaska: Eight Mile Lake located in Interior Alaska (63.8900° N, 149.2535° W), and Imnavait creek watershed located on the northern foothills of the Brooks range (68°37′ N, 149°18′ W).

Historical simulations are spanning the 2000 to 2015, and forced using climate simulations from the Climate Research Unit, time series 4.0. We ran 1,000 site level simulations for each model variable. The variables that are produced are gross primary productivity (GPP, in gC.m-2.m-1), net ecosystem exchange (NEE, gC.m-2.m-1), ecosystem respiration (RECO, gC/m2/m-1), active layer thickness (ALT, m), soil temperature (TLAYER,°C) at 5, 10, 40 cm depths, soil moisture (LWCLAYER, m-3/m-3) at 5, 10 cm depths, and snow depth (SNOWDEPTH, m), evapotransipiration(EET, mm/m2/time), potential evapotransipiration (PET, mm/m2/time), leaf area index (LAI, m2/m2), organic layer thickness (OLT, m). The data are stored as compiled csv files, with time as the index, and each model sample output stored in the columns. In addition, there is a postprocessing python script to demonstrate the step and workflow used to generate the individual csv files post processed from the raw model outputs stored as netcdfs.