The accurate retrieval of chlorophyll-a concentration (Chl-a) from ocean color satellite data is extremely challenging in turbid, optically complex coastal waters. Ariake Bay in Japan is a turbid semi-enclosed bay of great socio-economic significance, but it suffers from serious water quality problems, particularly due to red tide events. Chl-a derived from the MODerate resolution Imaging Spectroradiometer (MODIS) sensor on satellite Aqua in Ariake Bay was investigated, and it was determined that the causes of the errors were from inaccurate atmospheric correction and inappropriate in-water algorithms. To improve the accuracy of MODIS remote sensing reflectance (Rrs) in the blue and green bands, a simple method was adopted using in situ Rrs data. This method assumes that the error in MODIS Rrs(547) is small, and MODIS Rrs(412) can be estimated from MODIS Rrs(547) using a linear relation between in situ Rrs(412) and Rrs(547). We also showed that the standard MODIS Chl-a algorithm, OC3M, underestimated Chl-a, which was mostly due to water column turbidity. A new empirical switching algorithm was generated based on the relationship between in situ Chl-a and the blue-to-green band ratio, max(Rrs(443), Rrs(448)/Rrs(547), which was the same as the OC3M algorithm. The criterion of Rrs(667) of 0.005 sr⁻¹ was used to evaluate the extent of turbidity for the switching algorithm. The results showed that the switching algorithm performed better than OC3M, and the root mean square error (RMSE) of estimated Chl-a decreased from 0.414 to 0.326. The RMSE for MODIS Chl-a using the recalculated Rrs and the switching algorithm was 0.287, which was a significant improvement from the RMSE of 0.610, which was obtained using standard MODIS Chl-a. Finally, the accuracy of our method was tested with an independent dataset collected by the local Fisheries Research Institute, and the results revealed that the switching algorithm with the recalculated Rrs reduced the RMSE of MODIS Chl-a from 0.412 of the standard to 0.335.

A station of AErosol RObotic NETwork Ocean Color (AERONET-OC) has been set on the Ariake Observation Tower of Saga University on April, 2018, for verification of the Second generation Global Imager (SGLI)/Global Change Observation Mission-Climate (GCOM-C). Remote sensing reflectance (Rrs) observed by the AERONET-OC was used for verification of SGLI. SGLI Version 1 data underestimated the shortwave Rrs and Rrs (380) and Rrs(412) were mostly negative, while the estimation was improved by Version 2 with the correction of Rrs(412) to be positive. It was indicated that absorptive aerosol was influenced to SGLI atmospheric correction and caused the underestimation of Rrs. Simple linear correction method to improve shortwave Rrs also worked well for specifically Version 1 data. Water constituents, chlorophyll-a (Chl-a), total suspended matter (TSM) and colored dissolved organic matter (CDOM) of the SGLI were also verified by the ship observation data. All constituents were improved from Version 1 to 2 with the correction of Rrs, although Version 2 underestimated Chl-a and CDOM. Simple regression algorithms were also examined with in situ as well as SGLI Rrs data, and it indicated that more sophisticated algorithms may be required. Time series of water constituents derived from AERONET-OC and SGLI data were compared to river discharge and spring–neap tidal cycle. The results indicated that the behavior, such as the increase of Chl-a after river discharge and interaction of Chl-a and TSM with the spring–neap tidal cycles were observed, although estimations of water constituents were not completely separated by the algorithms.

We investigated the spatio-temporal variability of chlorophyll-a (Chl-a) and total suspended matter (TSM) associated with spring–neap tidal cycles in the Ariake Sea, Japan. Our study relied on significantly improved, regionally-tuned datasets derived from the ocean color sensor Moderate Resolution Imaging Spectroradiometer (MODIS) Aqua over a 16-year period (2002–2017). The results revealed that spring–neap tidal variations in Chl-a and TSM within this macrotidal embayment (the Ariake Sea) are clearly different regionally and seasonally. Generally, the spring–neap tidal variability of Chl-a in the inner part of the Ariake Sea was controlled by TSM for seasons other than summer, whereas it was controlled by river discharge for summer. On the other hand, the contribution of TSM to the variability of Chl-a was not large for two areas in the middle of Ariake Sea where TSM was not abundant. This study demonstrates that ocean color satellite observations of Chl-a and TSM in the macrotidal embayment offer strong advantages for understanding the variations during the spring–neap tidal cycle.

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.