Previous studies were mostly conducted based on two-dimensional feature space to monitor salinization, while studies on dense long-term salinization monitoring based on three-dimensional feature space have not been reported. Based on Landsat TM/ETM+/OLI images and three-dimensional feature space method, this study introduced six typical salinization surface parameters, including NDVI, salinity index, MSAVI, surface albedo, iron oxide index, wetness index to construct eight different feature space monitoring index. The optimal soil salinization monitoring index model was proposed base on field observed data and then the evolution process of salinization in Yellow River Delta (YRD) were analyzed and revealed during 1984–2022. The salinization monitoring index model of MSAVI-Albedo-I_Fe2O3 feature space had the highest accuracy with R² = 0.93 and RMSE = 0.678g/kg. The spatial distribution of salinization in YRD showed an increasing trend from inland southwest to coastal northeast and the salinization intensity showed an increasing trend during 1984–2022 due to the implements of agricultural measures such as planting salt-tolerant crops, microbial remediation and fertility improvement. The rate of salinization deterioration in the northeast part was greater than others. Zones of salinization improvement were mainly located in cultivated land of the southwest parts.

Accurate information on cropland extent is critical for scientific research and resource management. Several cropland products from remotely sensed datasets are available. Nevertheless, significant inconsistency exists among these products and the cropland areas estimated from these products differ considerably from statistics. In this study, we propose a hierarchical optimization synergy approach (HOSA) to develop a hybrid cropland map of China, circa 2010, by fusing five existing cropland products, i.e., GlobeLand30, Climate Change Initiative Land Cover (CCI-LC), GlobCover 2009, MODIS Collection 5 (MODIS C5), and MODIS Cropland, and sub-national statistics of cropland area. HOSA simplifies the widely used method of score assignment into two steps, including determination of optimal agreement level and identification of the best product combination. The accuracy assessment indicates that the synergy map has higher accuracy of spatial locations and better consistency with statistics than the five existing datasets individually. This suggests that the synergy approach can improve the accuracy of cropland mapping and enhance consistency with statistics.

Cropland fields are the basic spatial units for agricultural management, and information about their distribution is critical for analyzing agricultural investments and management. However, the extraction of cropland fields of smallholder farms is a challenging task because of their irregular shapes and diverse spectrum. In this paper, we proposed a new object-based Genetic Programming (GP) approach to extract cropland fields. The proposed approach used the multiresolution segmentation (MRS) method to acquire objects from a very high resolution (VHR) image, and extracted spectral, shape and texture features as inputs for GP. Then GP was used to automatically evolve the optimal classifier to extract cropland fields. The results show that the proposed approach has obtained high accuracy in two areas with different landscape complexities. Further analysis show that the GP approach significantly outperforms five commonly used classifiers, including K-Nearest Neighbor (KNN), Decision Tree (DT), Naïve Bayes (NB), Support Vector Machine (SVM), and Random Forest (RF). By using different numbers of training samples, GP can maintain high accuracy with any volume of samples compared to other classifiers.

Information on crop spatial distribution is essential for agricultural monitoring and food security. Classification with remote-sensing time series images is an effective way to obtain crop distribution maps across time and space. Optimal features are the precondition for crop classification and are critical to the accuracy of crop maps. Although several approaches are available for extracting spectral, temporal, and phenological features for crop identification, these methods depend heavily on domain knowledge and human experiences, adding uncertainty to the final crop classification. This study proposed a novel Genetic Programming (GP) approach to learning high-level features from time series images for crop classification to address this issue. We developed a new representation of GP to extend the GP tree’s width and depth to dynamically generate either fixed or flexible informative features without requiring domain knowledge. This new GP approach was wrapped with four classifiers, i.e., K-Nearest Neighbor (KNN), Decision Tree (DT), Naive Bayes (NB), and Support Vector Machine (SVM), and was then used for crop classification based on MODIS time series data in Heilongjiang Province, China. The performance of the GP features was compared with the traditional features of vegetation indices (VIs) and the advanced feature learning method Multilayer Perceptron (MLP) to show GP effectiveness. The experiments indicated that high-level features learned by GP improved the classification accuracies, and the accuracies were higher than those using VIs and MLP. GP was more robust and stable for diverse classifiers, different feature numbers, and various training sample sets compared with classification using VI features and the classifier MLP. The proposed GP approach automatically selects valuable features from the original data and uses them to construct high-level features simultaneously. The learned features are explainable, unlike those of a black-box deep learning model. This study demonstrated the outstanding performance of GP for feature learning in crop classification. GP has the potential of becoming a mainstream method to solve complex remote sensing tasks, such as feature transfer learning, image classification, and change detection.

Accurate assessments of the historical and current status of eco-environmental quality (EEQ) are essential for governments to have a comprehensive understanding of regional ecological conditions, formulate scientific policies, and achieve the United Nations Sustainable Development Goals (SDGs). While various approaches to EEQ monitoring exist, they each have limitations and cannot be used universally. Moreover, previous studies lack detailed examinations of EEQ dynamics and its driving factors at national and local levels. Therefore, this study utilized a remote sensing ecological index (RSEI) to assess the EEQ of China from 2001 to 2021. Additionally, an emerging hot-spot analysis was conducted to study the spatial and temporal dynamics of the EEQ of China. The degree of influence of eight major drivers affecting EEQ was evaluated by a GeoDetector model. The results show that from 2001 to 2021, the mean RSEI values in China showed a fluctuating upward trend; the EEQ varied significantly in different regions of China, with a lower EEQ in the north and west and a higher EEQ in the northeast, east, and south in general. The spatio-temporal patterns of hot/cold spots in China were dominated by intensifying hot spots, persistent cold spots, and diminishing cold spots, with an area coverage of over 90%. The hot spots were concentrated to the east of the Hu Huanyong Line, while the cold spots were concentrated to its west. The oscillating hot/cold spots were located in the ecologically fragile agro-pastoral zone, next to the upper part of the Hu Huanyong Line. Natural forces have become the main driving force for changes in China’s EEQ, and precipitation and soil sand content were key variables affecting the EEQ. The interaction between these factors had a greater impact on the EEQ than individual factors.