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.