This data set contains output from the dynamic vegetation model MC1, as modified to simulate future woody encroachment in the northern Great Plains. Simulations were done for the historical period (1895-2005) and the future period (2006-2100). Separate simulations were done for eastern and western portions of the region, with the eastern simulations using model parameters appropriate for Juniperus virginiana as the major evergreen needle-leaf life form, and the western simulations using model parameters appropriate for Pinus ponderosa as the major evergreen needle-leaf life form. Simulations in each portion were run for two A2 emissions scenario climate projections (CSIRO, representing moderate temperature increases and wetter conditions, and MIROC, representing very hot and dry conditions) crossed with 8 (eastern portion) or 6 (western portion) fire x grazing x tree regeneration capacity (eastern only) scenarios. Output variables provided on a yearly basis are potential evapotranspiration, live aboveground tree carbon and aboveground grass net primary production. Output variables provided as decadal averages are live aboveground tree carbon, tree leaf area index, soil available water for plant survival, surface runoff, potential evapotranspiration, streamflow, and actual evapotranspiration. Child records contain command files for running the model, model parameters, model input, and output from model runs for the equilibrium and spinup stages of model runs (precursors to running historical and future simulations).