Missing spatial data is one of the major concerns associated with the application of satellite data. The Data INterpolating Empirical Orthogonal Functions (DINEOF) method has been proven to be an effective tool for filling spatial gaps in various satellite data products. The Ariake Sea, which is a turbid coastal sea, shows the large spatial and temporal variability of chlorophyll-a (Chl-a) and total suspended matter (TSM). However, ocean color satellite data for this region usually have large gaps, which affects the accurate analysis of Chl-a and TSM variability. In this study, we applied the DINEOF method to fill the missing pixels from the regionally tuned Moderate Resolution Imaging Spectroradiometer (MODIS)-Aqua (hereafter, MODIS) Chl-a and MODIS-derived TSM datasets for the period 2002–2017. The validation results showed that the DINEOF reconstructed data were accurate and reliable. Furthermore, the Empirical Orthogonal Functions (EOF) analysis based on the reconstructed data was used to quantitatively analyze the spatial and temporal variability of Chl-a and TSM at both monthly and individual events of spring-neap tidal scales. The first three EOF modes of Chl-a showed seasonal variability mainly caused by precipitation, the sea surface temperature (SST), and river discharge for the first EOF mode and the sea level amplitude for the second. The first three EOF modes of TSM exhibited both seasonal and spring-neap tidal variability. The first and second EOF modes of TSM displayed spring-neap tidal variability caused by the sea level amplitude. The second EOF mode of TSM also showed seasonal variability caused by the sea level amplitude. In this study, we first applied the DINEOF method to reconstruct the satellite data and to capture the major spatial and temporal variability of Chl-a and TSM for the Ariake Sea. Our results demonstrate that the DINEOF method can reconstruct patchy oceanic color datasets and improve spatio-temporal variability analysis.

Shoreline changes are of great importance for evaluating the interaction between humans and ecosystems in coastal areas. They serve as a useful metric for assessing the ecological costs of socioeconomic developmental activities along the coast. In this paper, we present an assessment of shoreline changes along the eastern coast of mainland China from ~1990 to 2019 by applying a novel method recently developed by us. This method which we call the Nearest Distance Method (NDM) is used to make a detailed assessment of shorelines delineated from Landsat Thematic Mapper (TM) and Operational Land Imager (OLI) images. The results indicate a dramatic decline in natural shorelines that correspond to the rapid increase in the construction of artificial shorelines, driven by China’s economic growth. Of the entire coast of mainland China, the biggest change occurred along the Bohai Sea, where artificial shorelines expanded from 42.4% in ~1990 to 81.5% in 2019. Over this period, this study indicates that China lost > 60% of its biological shorelines, a trend that is especially worrisome because these include areas that were once biologically diverse and extremely rich. As anticipated, shoreline losses were greatest where regions of low economic value had been transformed to areas of higher economic value. Overall, this influence of human activities on shorelines in China is unprecedented. The repercussions of these changes on ecosystems, and the susceptibility of new shoreline developments to population growth and sea-level rise, need to be assessed urgently before additional changes are effected.

Airborne hyperspectral data play an important role in remote sensing of coastal waters. However, before their application, atmospheric correction is required to remove or reduce the atmospheric effects caused by molecular and aerosol scattering and absorption. In this study, we first processed airborne hyperspectral CASI-1500 data acquired on 4 May 2019 over the Uljin coast of Korea with Polymer and then compared the performance with the other two widely used atmospheric correction approaches, i.e., 6S and FLAASH, to determine the most appropriate correction technique for CASI-1500 data in coastal waters. Our results show the superiority of Polymer over 6S and FLAASH in deriving the R_rs spectral shape and magnitude. The performance of Polymer was further evaluated by comparing CASI-1500 R_rs data with those obtained from the MODIS-Aqua sensor on 3 May 2019 and processed using Polymer. The spectral shapes of the derived R_rs from CASI-1500 and MODIS-Aqua matched well, but the magnitude of CASI-1500 R_rs was approximately 0.8 times lower than MODIS R_rs. The possible reasons for this difference were time difference (1 day) between CASI-1500 and MODIS data, higher land adjacency effect for MODIS-Aqua than for CASI-1500, and possible errors in MODIS R_rs from Polymer.

The growing mismatch between ecosystem service (ES) supply and demand underscores the importance of thoroughly understanding their spatiotemporal patterns and key drivers to promote ecological civilization and sustainable development at the regional level in China. This study investigates six key ES indicators across mainland China—habitat quality (HQ), carbon sequestration (CS), water yield (WY), sediment delivery ratio (SDR), food production (FP), and nutrient delivery ratio (NDR)—by integrating a suite of analytical approaches. These include a spatiotemporal analysis of trade-offs and synergies in supply, demand, and their ratios; self-organizing maps (SOM) for bundle identification; and interpretable machine learning models. While prior research studies have typically examined ES at a single spatial scale, focusing on supply-side bundles or associated drivers, they have often overlooked demand dynamics and cross-scale interactions. In contrast, this study integrates SOM and SHAP-based machine learning into a dual-scale framework (grid and city levels), enabling more precise identification of scale-dependent drivers and a deeper understanding of the complex interrelationships between ES supply, demand, and their spatial mismatches. The results reveal pronounced spatiotemporal heterogeneity in ES supply and demand at both grid and city scales. Overall, the supply services display a spatial pattern of higher values in the east and south, and lower values in the west and north. High-value areas for multiple demand services are concentrated in the densely populated eastern regions. The grid scale better captures spatial clustering, enhancing the detection of trade-offs and synergies. For instance, the correlation between HQ and NDR supply increased from 0.62 (grid scale) to 0.92 (city scale), while the correlation between HQ and SDR demand decreased from −0.03 to −0.58, indicating that upscaling may highlight broader synergistic or conflicting trends missed at finer resolutions. In the spatiotemporal interaction network of supply–demand ratios, CS, WY, FP, and NDR persistently show low values (below −0.5) in western and northern regions, indicating ongoing mismatches and uneven development. Driver analysis demonstrates scale-dependent effects: at the grid scale, HQ and FP are predominantly influenced by socioeconomic factors, SDR and WY by ecological variables, and CS and NDR by climatic conditions. At the city level, socioeconomic drivers dominate most services. Based on these findings, nine distinct supply–demand bundles were identified at both scales. The largest bundle at the grid scale (B3) occupies 29.1% of the study area, while the largest city-scale bundle (B8) covers 26.5%. This study deepens the understanding of trade-offs, synergies, and driving mechanisms of ecosystem services across multiple spatial scales; reveals scale-sensitive patterns of spatial mismatch; and provides scientific support for tiered ecological compensation, integrated regional planning, and sustainable development strategies.