Vegetation phenology in temperate grasslands is highly sensitive to climate change. However, it is still unclear how the timing of vegetation phenology events (especially for autumn phenology) is altered in response to climate change across different grassland types. In this study, we investigated variations of the growing season start (SOS) and end (EOS), derived from Moderate Resolution Imaging Spectroradiometer (MODIS) data (2000–2016), for meadow steppe, typical steppe, and desert steppe in the Inner Mongolian grassland of Northern China. Using gridded climate data (2000–2015), we further analyzed correlations between SOS/EOS and pre-season average air temperature and total precipitation (defined as 90-day period prior to SOS/EOS, i.e., pre-SOS/EOS) in each grid. The results showed that both SOS and EOS occurred later in desert steppe (day of year (doy) 114 and 312) than in meadow steppe (doy 109 and 305) and typical steppe (doy 111 and 307); namely, desert steppe has a relatively late growing season than meadow steppe and typical steppe. For all three grasslands, SOS was mainly controlled by pre-SOS precipitation with the sensitivity being largest in desert steppe. EOS was closely connected with pre-EOS air temperature in meadow steppe and typical steppe, but more closely related to pre-EOS precipitation in desert steppe. During 2000–2015, SOS in typical steppe and desert steppe has significantly advanced by 2.2 days and 10.6 days due to a significant increase of pre-SOS precipitation. In addition, EOS of desert steppe has also significantly advanced by 6.8 days, likely as a result from the combined effects of elevated preseason temperature and precipitation. Our study highlights the diverse responses in the timing of spring and autumn phenology to preceding temperature and precipitation in different grassland types. Results from this study can help to guide grazing systems and to develop policy frameworks for grasslands protection.

Autumn phenology is a crucial indicator for identifying the alpine grassland growing season’s end date on the Qinghai-Tibet Plateau (QTP), which intensely controls biogeochemical cycles in this ecosystem. Although autumn phenology is thought to be mainly influenced by the preseason temperature, precipitation, and insolation in alpine grasslands, the relative contributions of these climatic factors on the QTP remain uncertain. To quantify the impacts of climatic factors on autumn phenology, we built stepwise linear regression models for 91 meteorological stations on the QTP using in situ herb brown-off dates, remotely sensed autumn phenological metrics, and a multi-factor climate dataset during an optimum length period. The results show that autumn precipitation has the most extensive influence on interannual variation in alpine grassland autumn phenology. On average, a 10 mm increase in autumn precipitation during the optimum length period may lead to a delay of 0.2 to 4 days in the middle senescence date (P < 0.05) across the alpine grasslands. The daily minimum air temperature is the second most important controlling factor, namely, a 1 °C increase in the mean autumn minimum temperature during the optimum length period may induce a delay of 1.6 to 9.3 days in the middle senescence date (P < 0.05) across the alpine grasslands. Sunshine duration is the third extensive controlling factor. However, its influence is spatially limited. Moreover, the relative humidity and wind speed also have strong influences at a few stations. Further analysis indicates that the autumn phenology at stations with less autumn precipitation is more sensitive to precipitation variation than at stations with more autumn precipitation. This implies that autumn drought in arid regions would intensely accelerate the leaf senescence of alpine grasslands. This study suggests that precipitation should be considered for improving process-based autumn phenology models in QTP alpine grasslands.

Revealing grassland growing season spatial patterns and their climatic controls is crucial for estimating the spatial heterogeneity of grassland productivity and carbon sequestration. In this study, we first used satellite‐derived normalized difference vegetation index data and a double logistic function to extract the start (SOS), end (EOS), and length of the growing season in midlatitude grasslands of the Northern Hemisphere during 1981–2014. Then, we verified the accuracy of satellite‐derived SOS and EOS using ground‐observed phenological records and gross primary production data at some locations. Moreover, we analyzed the spatial patterns of growing season indicators and their climatic controls. Results show that both SOS and EOS appear first in cool semidesert grasslands (CG), then in temperate grasslands (TG) and alpine grasslands (AG), and finally in warm semidesert grasslands (WG). A delaying tendency of SOS and EOS from north to south was identified in TG of North America. In contrast, an advancing tendency of SOS and EOS from north to south was detected in CG of Central and Western Asia. Further analysis indicates that a spatial opposite effect of spring temperature and precipitation triggers SOS in TG, whereas a spatial synergy effect of spring temperature and precipitation triggers SOS in CG of Asia, WG, and AG. Meanwhile, a spatial synergy effect of autumn temperature and precipitation triggers EOS for TG of North America and AG, whereas a spatial opposite effect of autumn temperature and precipitation determines EOS for CG.

Crop growth in different periods influences the final yield. This study started from the agronomic mechanism of yield formation and aimed to extract useful spectral characteristics in different phenological phases, which could directly describe the final yield and dynamic contributions of different phases to the yield formation. Hyperspectral information of the winter wheat canopy was acquired during three important phases (jointing stage, heading stage, and grain-filling stage). An enhanced 2D correlation spectral analysis method modified by mutual information was proposed to identify the sensitive wavebands. The selected wavebands performed well with good mechanism interpretation and close correlation with important crop growth parameters and main physiological activities related to yield formation. The quantitative contribution proportions of plant growth in three phases to the final yield were estimated by determining the coefficients of partial least square models based on full spectral information. They were then used as single-phase weight factors to merge the selected wavebands. The support vector machine model based on the weighted spectral dataset performed well in yield prediction with satisfactory accuracy and robustness. This result would provide rapid and accurate guidance for agricultural production and would be valuable for the processing of hyperspectral remote sensing data.