Urban green space takes a dominant role in alleviating the urban heat island (UHI) effect. Most investigations into the effects of cooling factors from urban green spaces on the UHI have evaluated the correlation between each factor and land surface temperature (LST) separately, and the contribution weights of various typical cooling factors in mitigating the thermal environment have rarely been analyzed. For this research, three periods of Landsat 8 data captured between 2014 and 2018 of Xuzhou during the summer and autumn seasons were selected along with corresponding meteorological and flux measurements. The mono-window method was employed to retrieve LST. Based on the characteristics of the vegetation and spatial features of the green space, eight factors related to green space were selected and computed, consisting of three indices that measure vegetation and five metrics that evaluate landscape patterns: vegetation density (VD), evapotranspiration (ET), green space shading degree (GSSD), patch area ratio (PLAND), largest patch index (LPI), patch natural connectivity (COHESION), patch aggregation (AI), and patch mean shape index distribution (SHPAE_MN). Linear regression and bivariate spatial autocorrelation analyses between each green space factor and LST showed that there were significant negative linear and spatial correlations between all factors and LST, which proved that the eight factors were all cooling factors. In addition, LST was strongly correlated with all factors (|r| > 0.5) except for SHPAE_MN, which was moderately correlated (0.3 < |r| < 0.5). Based on this, two principal components were extracted by applying principal component analysis with all standardized green space factors as the original variables. To determine the contribution weight of each green space factor in mitigating the urban heat island (UHI) effect, we multiplied the influence coefficient matrix of the initial variables with the standardized multiple linear regression coefficients between the two principal component variables and LST. The final results indicated that the vegetation indices of green space contribute more to the alleviation of the UHI than its landscape pattern metrics, and the contribution weights are ranked as VD ≥ ET > GSSD > PLAND ≈ LPI > COHESION > AI > SHAPE_MN. Our study suggests that increasing vegetation density is preferred in urban planning to mitigate urban thermal environment, and increasing broadleaf forests with high evapotranspiration and shade levels in urban greening is also an effective way to reduce ambient temperature. For urban green space planning, a priority is to multiply the regional green space proportion or the area of largest patches. Second, improving the connectivity or aggregation among patches of green space can enhance their ability to cool the surrounding environment. Altering the green space spatial shape is likely the least significant factor to consider.

Impervious surfaces (IPS) are the major source of urban heat island effect (UHI), and the relationships between IPS and land surface temperature (LST) have been widely studied. However, the spatial impact of landscape patterns of patches with different IPS density (IPSD) on the thermal environment remains largely unexplored. Based on three Landsat 8 images of the Xuzhou built-up area obtained in May and the corresponding ground observations from 2014 to 2020, the IPSD and LST maps were inversed through a linear spectral mixture analysis and mono-window algorithm, respectively. The landscape patterns of the five IPSD levels were characterized by four landscape-level and five class-level metrics. Finally, the spatial correlation between all landscape metrics and LST were analyzed using bivariate Moran’s I. The results were as follows: (1) The findings revealed that for the landscape-level metrics, LST had significant positive spatial correlations with Shannon’s diversity index (SHDI), Shannon’s evenness index (SHEI), and patch density (PD), while showing a significant negative correlation with contagion index (CONTAG), indicating that increasing the types, even distribution degree, and density of patches, or decreasing the aggregation degree of the five IPSD levels will lead to the enhancement of the thermal environment. (2) Furthermore, the class-level metrics of each IPSD level, percentage of landscape (PLAND), largest patch index (LPI), landscape shape index (LSI), aggregation index (AI), and patch cohesion index (COHESION) showed significant correlations and LST, which signified that the spatial characteristics of patch proportion, predominance degree, shape complexity, aggregation degree, and natural connectivity degree of each IPSD level are important factors affecting UHI. In addition, the spatial correlations between LST and class-level metrics were significantly positive for IPSD levels 4 and 5 with an evidently higher Moran’s I value, indicating that landscape patterns of IPSD levels 4 and 5 were the key factors in UHI enhancement. Furthermore, the impact weights of each class-level metric of IPSD levels 4 and 5 on LST were also analyzed by applying the principal component analysis and the multivariate regression standardization coefficient. These results reveal the importance and impact mechanism of the IPSD spatial patterns on UHI evolution, which may provide a valuable reference for future urban planning and climate management. This study also suggests that regional UHI can be mitigated by reducing the area proportion, natural connectivity, and shape complexity of high-density impervious surfaces.

The landscape patterns of urban green spaces have been proven to be important factors that affect urban thermal environments. However, the spatial effect of the landscape patterns of urban patches with different vegetation densities on land surface temperature (LST) has not been investigated in detail. In this study, the built-up area of Xuzhou City was taken as the study region, and the four phases of Landsat 8 images and their corresponding ground observations from 2014 to 2020 were selected as the basic data. Normalized spectral mixture analysis and an improved mono-window algorithm were used to invert the vegetation component fraction (VF) and LST maps of the study area, respectively, and the surface patches were classified into five levels according to the VF values, from low to high. Four landscape-level and five class-level metrics were then selected to represent the landscape characteristics of each VF-level patch. The tested values of 60 and 780 m were regarded as the best grain size and spatial extent, respectively, in the calculation of all landscape metrics of ALL VF-level patches (VFLM) using the moving-window method. The results of bivariate Moran’s I for VFLM and LST showed the following: (1) for landscape-level metrics, only the Shannon diversity index and patch diversity have substantial negative spatial correlations with LST (with average |Moran’s I| < 0.2), indicating that the types of VF levels and the number of patches exert weak negative effects on the thermal environment for a certain area; (2) for class-level metrics such as percentage of landscape, patch cohesion index, largest patch index, landscape shape index, and aggregation index, only the class-level metrics of sub-high VF (LV4) and extreme-high (LV5) VF levels patches have significant negative spatial correlations with LST (with high Moran’s I value, and high–high and low–high distributions in local indications of spatial association cluster maps), indicating that only the patches of high VF levels can effectively alleviate LST and that patch proportion, natural connectivity degree, predominance degree, shape complexity, and aggregation degree are important landscape factors for regulating the thermal environment. Principal component analysis and multiple linear regression were applied to determine the impact weights of the class-level VFLMs of LV4 and LV5 patches on LST, which revealed the contributions of these landscape metrics to mitigating the urban heat island effect (UHI). These results signify the importance of and differences in the spatial patterns of various VF-level patches for UHI regulation; these patterns can provide new perspectives and references for urban green space planning and climate management.

Water vapor is one of the primary greenhouse gases and significantly impacts the atmosphere. Water vapor is the most active meteorological element and varies rapidly in both the spatial and temporal domains. As a promising means, Global Navigation Satellite Systems (GNSS) tomography has been used to construct the 3D distribution of water vapor in high resolutions. Currently, in the most commonly used node parameterization approaches, the region for the 3D modeling has a preset fixed regular shape for all tomographic epochs. As a result, too many unknown parameters need to be estimated and thus to degrade the performance of the tomographic solution. In this study, an innovative node parameterization approach using a combination of three meshing techniques to dynamically adjust both the boundary of the tomographic region and the position of nodes at each tomographic epoch is proposed. The three meshing techniques were boundary extraction, Delaunay triangulation, and force‐displacement algorithm. The performance of the tomographic model resulting from the new approach was tested using one month GNSS data in May 2015 from the Hong Kong GNSS network and was compared against that of the conventional node parameterization approach. The reference for the validation of the accuracy of the test results were the radiosonde measurements from King's Park Meteorological Station (HKKP) in Hong Kong. Results showed that in terms of root‐mean‐square error the accuracy of the new approach significantly improved in comparison to the traditional approach.

High spatiotemporal resolution atmospheric water vapor can be retrieved using the Global Navigation Satellite System (GNSS) tomography technique, in which the remained ill-posed problem of the tomography system resulting from the acquisition geometry is a vital issue to be addressed. Remote sensing (RS) water vapor data, with high-resolution and global coverage, show great potential for retrieval of slant water vapor (SWV) observations to improve the tomographic geometrical distribution. In this article, we develop a GNSS-RS (GNSS combining RS) tomography model to fully exploit the value of observation signals from GNSS and RS measurements. The two key factors of retrieving the RS SWV are performed by calibrating the original precipitable water vapor (PWV) images and adding the tropospheric horizontal gradients. The results reveal that when introducing the RS SWV observations into the tomography model, the acquisition geometry is significantly improved, with the average rate of voxels crossed by rays from 62% to 95% and the mean number of observation signals from 395 to 508 during the tomographic periods. Independent radiosonde data are used to validate the tomographic water vapor fields. The mean root-mean-square error (RMSE) and bias of the water vapor profiles derived from GNSS-RS solutions are decreased by 28% and 45% with respect to the GNSS-only results, respectively. Such improvements highlight that GNSS-RS troposphere tomography has significant potential to improve the reconstruction of the atmospheric water vapor fields.

Global Navigation Satellite System (GNSS) tomography has developed into an efficient tool for sensing the high spatiotemporal variability of atmospheric water vapor. The integration of GNSS top signals and side rays for tropospheric tomography systems using a novel height factor model (HFM) is proposed and discussed in this paper. Within the HFM, the sectional slant wet delay (SWD) of inside signals (the part of the side signal inside the tomography area), which is considered a key factor for modeling side rays, is separated into isotropic and anisotropic components. Correspondingly, two height factors are defined to calculate the isotropic and anisotropic part of tropospheric delays in the HFM. In addition, the dynamic tomography top boundary is first analyzed and determined based on 30-year radiosonde data to reasonably divide signals into top and side rays. Four special experimental schemes based on different tomography regions of Hong Kong are performed to assess the proposed HFM method, the results of which show increases of 33.42% in the mean utilization of rays, as well as decreases of 0.46 g/m³ in the average root mean square error (RMSE), compared to the traditional approach, revealing the improvement of tomography solutions when the side signals are included in the modeling. Furthermore, compared with the existing correction model for modeling side rays, the water vapor profiles retrieved from the proposed improved model are closer to the radiosonde data, which highlights the advantages of the proposed HFM for optimizing the GNSS tomography model.