Land Surface Phenologies of the Northern Great Plains: Possible Futures Arising From Land and Climate Change

Land cover change across the Northern Great Plains of North America over the past three decades has been driven by changes in agricultural management (conservation tillage; irrigation), government incentives (Conservation Reserve Program; subsidies to grain-based ethanol), crop varieties (cold-hardy soybean), and market dynamics (increasing world demand). Climate change across the Northern Great Plains over the past three decades has been evident in trends toward earlier warmth in the spring and a longer frost-free season. Together these land and climate changes induce shifts in local and regional land surface phenologies (LSPs). Any significant shift in LSP may correspond to a significant shift in evapotranspiration, with consequences for regional hydrometeorology. We explored possible future scenarios involving land use and climate change in six steps. First, we defined the nominal draw areas of current and future biorefineries in North Dakota, South Dakota, Nebraska, Minnesota, and Iowa and masked those land cover types within the draw areas that were unlikely to change to agricultural use (open water, settlements, forests, etc.). Second, we estimated the proportion of corn and soybean remaining within the masked draw areas using MODIS-derived crop maps. Third, in each draw area, we modified LSPs to simulate crop changes for a control and two treatment scenarios. In the control, we used LSP profiles identified from MODIS Collection 5 NBAR data. In one treatment, we increased the proportion of tallgrass LSPs in the draw areas to represent widespread cultivation of a perennial cellulosic crop, like switchgrass. In a second treatment, we increased the proportion of corn LSPs in the draw areas to represent increased corn cultivation. Fourth, we characterized the seasonal progression of the thermal regime associated with the LSP profiles using MODIS Land Surface Temperature (LST) products. Fifth, we modeled the LSP profile as a quadratic function of accumulated growing degree-days based on the LST time series. Sixth, we used representative IPCC AR4 mid-century projections to force the quadratic models and produce possible future LSPs. The resulting shifts in potential peak vegetation to earlier dates indicate potential seasonal shifts in evapotranspiration.