Winter Wheat Green-up Date Variation and its Diverse Response on the Hydrothermal Conditions over the North China Plain, Using MODIS Time-Series Data

Vegetation phenology plays a critical role in the dynamic response of terrestrial ecosystems to climate change. However, the relationship between the phenology of winter wheat and hydrothermal factors is inadequate, especially in typical agricultural areas. In this study, the possible effects of preseason climate changes on the green-up date (GUD) of winter wheat over the North China Plain (NCP) was investigated, using the MODIS EVI 8-day time-series data from 2000 to 2015, as well as the concurrent monthly mean temperature (T-m), mean maximum (T-max) and minimum temperature (T-min) and total precipitation (TP) data. Firstly, we quantitatively identified the time lag effects of winter wheat GUD responses to different climatic factors; then, the major driving factors for winter wheat GUD were further explored by applying multiple linear regression models. The results showed that the time lag effects of winter wheat GUD response to climatic factors were site- and climatic parameters-dependent. Negative temperature effects with about a 3-month time lag dominated in most of the NCP, whereas positive temperature effects with a zero-month lag were most common in some of the southern parts. In comparison, total precipitation had a negative zero-month lag effect in the northern region, but two lagged months occurred in the south. Regarding the time lag effects, the explanation power of climatic factors improved relatively by up to 77%, and the explanation area increased by 41.20%. Additionally, change in winter wheat GUD was primarily determined by temperature rather than by TP, with a marked spatial heterogeneity of the T-max and T-min effect. Our results confirmed different time lag effects from different climatic factors on phenological processes in spring, and further suggested that both T-max and T-min should be considered to improve the performance of spring phenology models.