With the increasing population and continuation of climate change, an adequate food supply is vital to economic development and social stability. Winter crops are important crop types in China. Changes in winter crops planting areas not only have a direct impact on China’s production and economy, but also potentially affects China’s food security. Therefore, it is necessary to obtain information on the planting of winter crops. In this study, we use the time series data of individual pixels, calculate the temporal statistics of spectral bands and the vegetation indices of optical data based on the phenological characteristics of specific vegetation or crops and record them in the time series data, and apply decision trees and rule-based algorithms to generate annual maps of winter crops. First, we constructed a dataset combining all the available images from Landsat 7/8 and Sentinel-2A/B. Second, we generated an annual map of land cover types to obtain the cropland mask in 2019. Third, we generated a time series of a single cropland pixel, and calculated the phenological indicators for classification by extracting the differences in phenological characteristics of different crops: these phenological indicators include SOS (start of season), SDP (start date of peak), EOS (end of season), GUS (green-up speed) and GSL (growing-season length). Finally, we identified winter crops in 2019 based on their phenological characteristics. The main advantages of the phenology-based algorithm proposed in this study include: (1) Combining multiple sensor data to construct a high spatiotemporal resolution image collection. (2) By analyzing the whole growth season of winter crops, the planting area of winter crops can be extracted more accurately, and (3) the phenological indicators of different periods are extracted, which is conducive to monitoring winter crop planting information and seasonal dynamics. The results show that the algorithm constructed in this study can accurately extract the planting area of winter crops, with user, producer, overall accuracies and Kappa coefficients of 96.61%, 94.13%, 94.56% and 0.89, respectively, indicating that the phenology-based algorithm is reliable for large area crop classification. This research will provide a point of reference for crop area extraction and monitoring.

The spatio-temporal change of the surface water is very important to agricultural, economic, and social development in the Hetao Plain, as well as the structure and function of the ecosystem. To understand the long-term changes of the surface water area in the Hetao Plain, we used all available Landsat images (7534 scenes) and adopted the modified Normalized Difference Water Index (mNDWI), Enhanced Vegetation Index (EVI), and Normalized Difference Vegetation Index (NDVI) to map the open-surface water from 1989 to 2019 in the Google Earth Engine (GEE) cloud platform. We further analyzed precipitation, temperature, and irrigated area, revealing the impact of climate change and human activities on long-term surface water changes. The results show the following. (1) In the last 31 years, the maximum, seasonal, and annual average water body area values in the Hetao Plain have exhibited a downward trend. Meanwhile, the number of maximum, seasonal, and permanent water bodies displayed a significant upward trend. (2) The variation of the surface water area in the Hetao Plain is mainly affected by the maximum water body area, while the variation of the water body number is mainly affected by the number of minimum water bodies. (3) Precipitation has statistically significant positive effects on the water body area and water body number, which has statistically significant negative effects with temperature and irrigation. The findings of this study can be used to help the policy-makers and farmers understand changing water resources and its driving mechanism and provide a reference for water resources management, agricultural irrigation, and ecological protection.

Drought is one of the most complex and least-understood environmental disasters that can trigger environmental, societal, and economic problems. To accurately assess the drought conditions in the Yellow River Basin, this study reconstructed the Land Surface Temperature (LST) using the Annual Temperature Cycle (ATC) model and the Normalized Difference Vegetation Index (NDVI). The Temperature Condition Index (TCI), Vegetation Condition Index (VCI), Vegetation Health Index (VHI), and Temperature-Vegetation Drought Index (TVDI), which are four typical remote sensing drought indices, were calculated. Then, the air temperature, precipitation, and soil moisture data were used to evaluate the applicability of each drought index to different land types. Finally, this study characterized the spatial and temporal patterns of drought in the Yellow River Basin from 2003 to 2019. The results show that: (1) Using the LST reconstructed by the ATC model to calculate the drought index can effectively improve the accuracy of drought monitoring. In most areas, the reconstructed TCI, VHI, and TVDI are more reliable for monitoring drought conditions than the unreconstructed VCI. (2) The four drought indices (TCI, VCI, VH, TVDI) represent the same temporal and spatial patterns throughout the study area. However, in some small areas, the temporal and spatial patterns represented by different drought indices are different. (3) In the Yellow River Basin, the drought level is highest in the northwest and lowest in the southwest and southeast. The dry conditions in the Yellow River Basin were stable from 2003 to 2019. The results in this paper provide a basis for better understanding and evaluating the drought conditions in the Yellow River Basin and can guide water resources management, agricultural production, and ecological protection of this area.

Cropping intensity is a key indicator for evaluating grain production and intensive use of cropland. Timely and accurately monitoring of cropping intensity is of great significance for ensuring national food security and improving the level of national land management. In this study, we used all Sentinel-2 images on the Google Earth Engine cloud platform, and constructed an improved peak point detection method to extract the cropping intensity of a heterogeneous planting area combined with crop phenology. The crop growth cycle profiles were extracted from the multi-temporal normalized difference vegetation index (NDVI) and land surface water index (LSWI) datasets. Results show that by 2020, the area of single cropping, double cropping, and triple cropping in the Henan Province are 52,236.9 km2, 74,334.1 km2, and 1927.1 km2, respectively; the corresponding producer accuracies are 86.12%, 93.72%, and 91.41%, respectively; the corresponding user accuracies are 88.99%, 92.29%, and 71.26%, respectively. The overall accuracy is 90.95%, and the Kappa coefficient is 0.81. Using the sown area in the statistical yearbook data of cities in the Henan Province to verify the extraction results of this paper, the R2 is 0.9717, and the root mean square error is 1715.9 km2. This study shows that using all the Sentinel-2 data, the phenology algorithm, and cloud computing technology has great potential in producing a high spatio-temporal resolution dataset for crop remote sensing monitoring and agricultural policymaking in complex planting areas.

Double cropping is an important cropping system in China, with more than half of China’s cropland adopting the practice. Under the background of global climate change, agricultural policies, and changing farming practices, double-cropping area has changed substantially. However, the spatial-temporal dynamics of double cropping is poorly understood. A better understanding of these dynamics is necessary for the northern limit of double cropping (NLDC) to ensure food security in China and the world and to achieve zero hunger, the second Sustainable Development Goal (SDG). Here, we developed a phenology-based algorithm to identify double-cropping fields by analyzing time-series Moderate Resolution Imaging Spectroradiometer (MODIS) images during the period 2000–2020 using the Google Earth Engine (GEE) platform. We then extracted the NLDC using the kernel density of pixels with double cropping and analyzed the spatial-temporal dynamics of NLDC using the Fishnet method. We found that our algorithm accurately extracted double-cropping fields, with overall, user, and producer accuracies and Kappa coefficients of 95.97%, 96.58%, 92.21%, and 0.91, respectively. Over the past 20 years, the NLDC generally trended southward (the largest movement was 66.60 km) and eastward (the largest movement was 109.52 km). Our findings provide the scientific basis for further development and planning of agricultural production in China.

As the land use issue, caused by urban shrinkage in China, is becoming more and more prominent, research on urban shrinkage and expansion has become particularly challenging and urgent. Based on the points of interest (POI) data, this paper redefines the scope, quantity, and area of natural cities by using threshold methods, which accurately identify the shrinkage and expansion of cities in the Yellow River affected area using night light data in 2013 and 2018. The results show that: (1) there are 3130 natural cities (48,118.75 km2) in the Yellow River affected area, including 604 shrinking cities (8407.50 km2) and 2165 expanding cities (32,972.75 km2). (2) The spatial distributions of shrinking and expanding cities are quite different. The shrinking cities are mainly located in the upper Yellow River affected area, except for the administrative cities of Lanzhou and Yinchuan; the expanding cities are mainly distributed in the middle and lower Yellow River affected area, and the administrative cities of Lanzhou and Yinchuan. (3) Shrinking and expanding cities are typically smaller cities. The research results provide a quick data supported approach for regional urban planning and land use management, for when regional and central governments formulate the outlines of urban development monitoring and regional planning.

As a common natural disaster, drought can significantly affect the agriculture productivity and human life. Compared to Southeast China, Northwest China is short of water year-round and is the most frequent drought disaster area in China. Currently, there are still many controversial issues in drought monitoring of Northwest China in recent decades. To further understand the causes of changes in drought in Northwest China, we chose Shaanxi, Gansu, and Ningxia provinces (SGN) as our study area. We compared the spatiotemporal characteristics of drought intensity and frequency in Northwest China from 2003 to 2020 showed by the Standardized Precipitation Index (SPI), Vegetation Condition Index (VCI), Temperature Condition Index (TCI), Vegetation Health Index (VHI), Normalized Vegetation Supply Water Index (NVSWI), Soil Moisture Condition Index (SMCI), and Soil Moisture Agricultural Drought Index (SMADI). All of these indices showed a wetting trend in the SGN area from 2003 to 2020. The wetting trend of the VCI characterization is the most obvious (R2 = 0.9606, p < 0.05): During the period 2003–2020, the annual average value of the VCI in the SGN region increased from 28.33 to 71.61, with a growth rate of 153.57%. The TCI showed the weakest trend of wetting (R2 = 0.0087), with little change in the annual average value in the SGN region. The results of the Mann–Kendall trend test of the TCI indicated that the SGN region experienced a non-significant (p > 0.05) wetting trend between 2003 and 2020. To explore the effectiveness of different drought indices, we analyzed the Pearson correlation between each drought index and the Palmer Drought Severity Index (PDSI). The PDSI can not only consider the current water supply and demand situation but also consider the impact of the previous dry and wet conditions and their duration on the current drought situation. Using the PDSI as a reference, we can effectively verify the performance of each drought index. SPI-12 showed the best correlation with PDSI, with R values greater than 0.6 in almost all regions and p values less than 0.05 within one-half of the study area. SMADI had the weakest correlation with PDSI, with R values ranging −0.4~−0.2 and p values greater than 0.05 in almost all regions. The results of this study clarified the wetting trend in the SGN region from 2003 to 2020 and effectively analyzed the differences in each drought index. The frequency, duration, and severity of drought are continuously reduced; this helps us to have a more comprehensive understanding of the changes in recent decades and is of significance for the in-depth study of drought disasters in the future.

The rapid development of the web geographic information system (Web GIS) has promoted new vitality in high school geography education, relieved the stress of geography teachers caused by software and technical problems, and made it possible for teachers to devote more energy to geography teaching and research activities. Natural disaster education is not only an important part of the geography curriculum, but also an indispensable aspect of education for sustainable development (ESD) for high school students. The application of Web GIS in the dynamic monitoring, forecast, and early warning of natural disasters is becoming more experienced. Therefore, the application of Web GIS in natural disaster education is quite feasible. How to build a bridge between them is the purpose of this paper. Thus, the paper selects ArcGIS Online, which is not limited by time and space, and analyzes several functions that apply it to geography teaching. These include smart mapping, story maps, 3D web maps, and mobile GIS. Meanwhile, it analyzes the knowledge structure of “natural disasters” in Chinese geography textbooks to guide the subsequent case design. Then, the Web GIS inquiry-based teaching case is formed based on “7.20 Zhengzhou Torrential Rain”. It contains knowledge about natural disasters and designs from many aspects, such as the causes, manifestations, and prevention and control of disasters. The discussion identifies a range of specific educational benefits of applying Web GIS to natural disaster education for teachers and schools. Ultimately, it can provide some reference values for geography teachers and other developers to explore curriculum resources and create quality educational models.