Remote sensing (RS)-derived vegetation indices (VIs) with medium and high spatial resolution have emerged as a promising dataset for fine-scale ecosystem modeling and agricultural monitoring at local or global scales. Before they can be used as reliable inputs for other research, conducting in situ measurements for validation is very critical. However, the spatial heterogeneity due to the diversity of land cover and its spatial organization in the landscape increases the uncertainty of validation, so design of optimal sampling is an important basis for the reliability of the validation. In this paper, we propose an integrative stratified sampling strategy (INTEG-STRAT) based on normalized difference vegetation index (NDVI) data as prior knowledge. The basic idea is to realize a sampling optimization by determining the optimal combination of the spatial sampling method (e.g., simple random sampling (SRS), spatial system sampling (SYS), stratified sampling, generalized random tessellation stratified (GRTS), balanced acceptance sampling (BAS)) and spatial stratification scheme with an objective rule. The objective rule in this paper is to minimize the root mean square error (RMSE) of 10-fold cross validation between estimated values (sample are not included) and the corresponding values on prior knowledge. Relative precision, correlation coefficient, and RMSE are used to compare the effectiveness of the proposed sampling strategy with each sampling method without considering sampling optimization. After comparing, we find that the INTEG-STRAT requires fewer samples to become stable and has higher accuracy. At site 1, when the correlation coefficient between NDVI image and the simulated NDVI surface reached 80%, INTEG-STRAT needed only 70 sampling points while other methods require more sampling points. At the same time, INTEG-STRAT strategy has a smaller RMSE between the estimated values and the corresponding values on prior knowledge image. In general, INTEG-STRAT is an effective method in the selection of representative samples to support the validation of vegetation indices products with medium and high spatial resolution.

The 44 years (1979–2022) of satellite-derived sea ice extent in the Arctic and Antarctic reveals the details and new trends in the process of polar sea ice coverage changes. The speed of Arctic sea ice extent reduction and the interannual difference significantly increased after 2004. Trend analysis suggests that the Arctic Ocean may experience an ice-free period around 2060. The maximum anomaly of Arctic sea ice extent has gradually transitioned from September to October, indicating a trend of prolonged melting period. The center of gravity of sea ice in the Arctic Ocean is biased towards the Pacific side, and the spatial distribution pattern of sea ice is greatly influenced by the Atlantic warm current. The dynamism of the sea ice extent on the Atlantic side is significantly greater than in other regions. Since 2014, the Antarctic sea ice extent has shifted from slow growth to a rapid decreasing trend; the sea ice extent reached a historical minimum in 2022, decreasing by 2.02 × 106 km² compared to 2014. The Antarctic experiences seven months of ice growth each year and five months of ice melting period, the annual change patterns of sea ice extent in the Arctic and Antarctic are slightly different.