The concept of ecosystem services (ES) supply and demand has attracted increasing attention in science and policy making because it effectively links ecosystem services to human well-being. The imbalance of ES supply and demand in urban areas has become a key issue in regional sustainable development. In this context, we calculated ES supply and demand for Wuhan City, China, using the ES supply and demand ratio (ESDR) and the comprehensive ES supply–demand ratio (CESDR) to express the relationship between ES supply and demand. Ecological zoning was proposed according to the spatial differentiation of the ES supply–demand relationship, and policy recommendations are made. The results show that from the perspective of total ES supply and demand, the water yield supply (S_WY), grain yield supply (S_GY), and recreation services supply (S_RS) are greater than the water yield demand (D_WY), grain yield demand (D_GY), and recreation services demand (D_RS), and that the climate regulation supply (S_CR) is less than the climate regulation demand (D_CR). From a spatial perspective, there are imbalances and mismatches in ES supply and demand, especially in urban central areas. The values of S_WY, S_GY, S_CS, and S_RS per unit area are less than their respective demand values, and the area of mismatch has expanded with the gradual increase of the built-up area. The spatial pattern of ES supply and demand is circular, with the form of “deficit zone–relative equilibrium zone–surplus zone”, which corresponds to “urban central area–near suburbs–distant suburbs and rural areas”.

Current remote sensing-based aerosol optical depth (AOD) products have coarse spatial resolutions, which are useful for studies at continental and global scales, but unsatisfactory for local scale applications, such as urban air pollution monitoring. In this study, we investigated the possibility of using Landsat imagery to develop high-resolution AOD estimations at 30 m based on machine learning algorithms. We assessed the performance of six machine learning algorithms, including Extreme Gradient Boosting, Random Forest, Cascade Random Forest, Gradient Boosted Decision Trees, Extremely Randomized Trees, and Multiple Linear Regression. To obtain accurate AOD estimations, we used prior knowledge from multiple sources as inputs to the machine learning models, including the Global Land Surface Satellite (GLASS) albedo, the 1-km AOD product from MODIS data using the Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm, and meteorological and surface elevation data. A total of 13,624 AOD measurements from Aerosol Robotic Network (AERONET) sites were used for model training and validation. We found that all six algorithms exhibited good performance, with R² values ranging from 0.73 to 0.78 and AOD root-mean-square errors (RMSE) ranging from 0.089 to 0.098. The extremely randomized trees algorithm, however, demonstrated marginally superior performance as compared to the other algorithms; hence, it was used to produce AOD estimates at a 30 m resolution for one Landsat scene coving Beijing in 2013–2019. Through a comparison with overlapping AERONET observations, a high level of accuracy was achieved, with an R² = 0.889 and an RMSE = 0.156. Our method can be potentially used to generate a global high-resolution AOD dataset based on Landsat imagery.

Recent studies have shown that the high-resolution satellite Landsat-8 has the capability to retrieve aerosol optical depth (AOD) over urban areas at a 30 m spatial resolution. However, its long revisiting time and narrow swath limit the coverage and frequency of the high resolution AOD observations. With the increasing number of Earth observation satellites launched in recent years, combining the observations of multiple satellites can provide higher temporal-spatial coverage. In this study, a fusing retrieval algorithm is developed to retrieve high-resolution (30 m) aerosols over urban areas from Landsat-8 and Sentinel-2 A/B satellite measurements. The new fusing algorithm was tested and evaluated over Beijing city and its surrounding area in China. The validation results show that the retrieved AODs show a high level of agreement with the local urban ground-based Aerosol Robotic Network (AERONET) AOD measurements, with an overall high coefficient of determination (R²) of 0.905 and small root mean square error (RMSE) of 0.119. Compared with the operational AOD products processed by the Landsat-8 Surface Reflectance Code (LaSRC-AOD), Sentinel Radiative Transfer Atmospheric Correction code (SEN2COR-AOD), and MODIS Collection 6 AOD (MOD04) products, the AOD retrieved from the new fusing algorithm based on the Landsat-8 and Sentinel-2 A/B observations exhibits an overall higher accuracy and better performance in spatial continuity over the complex urban area. Moreover, the temporal resolution of the high spatial resolution AOD observations was greatly improved (from 16/10/10 days to about two to four days over globe land in theory under cloud-free conditions) and the daily spatial coverage was increased by two to three times compared to the coverage gained using a single sensor.

Ultrahigh resolution fine particulate matter (PM_2.5) mass concentration remote sensing products are crucial for atmospheric environmental monitoring, pollution source verification, health exposure risk assessment, and other fine-scale applications in urban environments. This study developed an ultrahigh resolution retrieval algorithm to estimate 30 m resolution PM_2.5 mass concentrations over urban areas from Landsat-8 and Sentinel-2A/B satellite measurements. The algorithm utilized aerosol optical depth (AOD) products retrieved from the Landsat-8 OLI and Sentinel-2 MSI measurements from 2017 to 2020, combined with multi-source auxiliary data to establish a PM_2.5-AOD relationship model across China. The results showed an overall high coefficient of determination (R²) of 0.82 and 0.76 for the model training accuracy based on samples and stations, respectively. The model prediction accuracy in Beijing and Wuhan reached R² values of 0.86 and 0.85. Applications in both cities demonstrated that ultrahigh resolution PM_2.5 has significant advantages in resolving fine-scale spatial patterns of urban air pollution and pinpointing pollution hotspots. Furthermore, an analysis of point source pollution at a typical heavy pollution emission enterprise confirmed that ultrahigh spatial resolution PM_2.5 can accurately identify the diffusion trend of point source pollution, providing fundamental data support for refined monitoring of urban air pollution and air pollution prevention and control.

The geostationary satellite Himawari-8, possessing the Advanced Himawari Imager (AHI), which features 16 spectral bands from the visible to infrared range, is suitable for aerosol observations. In this study, a new algorithm is introduced to retrieve aerosol optical depth (AOD) over land at a resolution of 2 km from the AHI level 1 data. Considering the anisotropic effects of complex surface structures over land, Moderate Resolution Imaging Spectroradiometer (MODIS) bidirectional reflectance distribution function (BRDF) model parameters product (MCD19A3) is used to calculate the surface reflectance for Himawari-8’s view angle and band. In addition, daily BRDF model parameters are calculated in areas with dense vegetation, considering the rapid variation of surface reflectance caused by vegetation growth. Moreover, aerosol models are constructed based on long duration Aerosol Robotic Network (AERONET) single scattering albedo (SSA) values to stand for aerosol types in the retrieval algorithm. The new algorithm is applied to AHI images over Bohai Rim region from 2018 and is evaluated using the newest AERONET version 3 AOD measurements and the latest MODIS collection 6.1 AOD products. The AOD retrievals from the new algorithm show good agreement with the AERONET AOD measurements, with a correlation coefficient of 0.93 and root mean square error (RMSE) of 0.12. In addition, the new algorithm increases AOD retrievals and retrieval accuracy compared to the Japan Aerospace Exploration Agency (JAXA) aerosol products. The algorithm shows stable performance during different seasons and times, which makes it possible for use in climate or diurnal aerosol variation studies.