1. Introduction
Greening land has significant benefits in terms of reducing the land’s thermal properties and regulating ecosystem services. In hot summers, people are often attracted to the cooler climate of coastal cities [
1]. However, in the context of rapid urbanization [
2], reinforced concrete buildings may introduce the urban heat wave agglomeration effect and form a thermal discomfort region for residents in coastal urban areas [
3,
4,
5]. A reasonable layout of urban green spaces can usually alleviate damage from extreme summer temperatures and avoid discomfort, such as heatstroke, for local residents or tourists, thereby improving ecosystem services in this region [
6,
7]. Currently, China is still undergoing a rapid urbanization process, shifting from the old mode of unilateral pursuit of urban expansion to the new mode of rational urban land planning with a harmonious coexistence between people and natural environments [
8,
9]. In the context of the current urban construction of sponge cities, livable cities, innovative cities, and ecological cities in China, greening of land plays an increasingly important role in providing comfortable living environments and regulating ecosystem services [
10]. Especially in the face of the aim of reaching a carbon peak in 2030 and carbon neutrality in 2060 in China, the city’s carbon digestion capacity further requires the improvement of green land distribution patterns [
11,
12,
13]. Studies monitoring the spatiotemporal heterogeneity changes in green spaces and analyzing their impact on residents’ thermal comfort from the perspective of land-surface thermal properties have become hot topics.
Scholars have conducted extensive research on the spatiotemporal heterogeneity of the process of greening land and its change [
14,
15]. Early green land area surveys mainly relied on manpower and manual measuring equipment such as level gauges and theodolites. This type of measurement required many professionals, and the efficiency of the measurement was not high [
16,
17,
18]. With the development of Earth remote sensing observation technology, satellite sensors can be used to obtain Earth’s land surface types, such as green land cover [
19,
20,
21]. The imaging resolution for monitoring land use and cover changes has increased from 1000 m to 300 m, 30 m, 10 m, or even sub-meter. In the context of current big data and cloud platforms, land cover monitoring has become more popular at the local, regional, and global scales [
22]. In this process, many green land mapping products have emerged using remote sensing satellites, such as the global 30 m land mapping products from the National Bureau of Surveying and Mapping in China, the 10 m resolution global land cover mapping, 2017, from Tsinghua University, and the vector dataset of land use/cover changes across China from the Chinese Academy of Sciences [
23]. Based on these green land data monitoring products, global land cover was becoming increasingly green (i.e., greening Earth) [
24], which was especially prominent in China and India. China accounted for 25% of the global net increase in leaf area while accounting for only 6.6% of global vegetated area. The greening in China was mainly from forests (42%) and croplands (32%), but in India, the increased greening of land came mainly from croplands (82%), with a minor contribution from forests (4.4%) [
24]. In urban areas, the cities were also becoming greener, especially in arid and semi-arid regions in Africa [
25]. Meanwhile, among these typical land use products, the vector land product released by the Chinese Academy of Sciences has multiperiod data, including for 1980, 1995, 2000, 2005, 2008, 2010, 2013, 2015, 2018, and 2020, which can be used for a continuous analysis of the land status and its change trend. These data contained 25 land classification types, showing a very extensive land cover mapping ability. In many previous studies, these data achieved good land mapping and land change analysis results in all of China in both regional and local scale investigations [
26,
27,
28]. Considering the relevance of multiple periods and the diversity of land classifications from the land use/cover data from the Chinese Academy of Sciences, these land data were applied in this study.
From the perspective of land surface thermal properties, land use/cover changes alter the surface thermal radiance, causing cooling or warming as well as energy effects, which can be used to assess the thermal comfort of residents. Previously, studies mainly focused on evaluating the effects of a building’s thermal environment on workers [
29]. Subsequently, the Universal Thermal Climate Index (UTCI) was proposed, in order to investigate the relationship between human settlements and factors of land surface convection, thermal radiation, and water evaporation [
30,
31]. Additionally, the impact of greening of land on the comfort of human settlements was continuously explored, considering the impact of vegetation’s horizontal and vertical characteristics on the environment. Specifically, vegetation has a natural physical shielding property, which affects air flow and heat exchange through horizontal blocking effects, resulting in a decrease in ambient temperature, and vegetation further intercepts solar radiation through vertical shading, creating relatively low-temperature areas under vegetation through shading [
31,
32]. Then, scholars conducted research on the comfort of the living environment from a perspective of greening land, using the temperature and humidity index, comfort index, and wind efficiency index in regions where urban spaces were being greened [
33,
34]. The popularization of remote sensing technology has led these studies to move from statistics to spatialized land mapping results, such as spatialized land-surface thermal property maps of land surface temperature, air temperature, net radiation, the Bowen ratio, latent heat flux, available flux, and sensible heat flux [
35,
36,
37,
38], and thus, this has become a popular research approach for investigating the impact of land use changes on surface thermal properties.
The central coastal region of China attracts tourists from all over the country due to the advantages of the natural environment such as a suitable climate and high levels of green vegetation. Currently, the investigation of changes following greening of land in this region is insufficient, which limits the understanding of the impact of green spaces on residents’ comfort in this region. Rizhao city is located in the central coastal region of China. Since 2000, the urbanization process has been accelerating. In 2009, Rizhao won “the United Nations Habitat Award” for its outstanding achievements in “improving housing and infrastructure, and building a green home”, and it was the only region in China to receive this honor that year. Subsequently, Rizhao entered a faster stage of urbanization and obtained more Chinese awards such as that of being a “greening city”. Considering the suitable natural climate and high coverage of green land, Rizhao city is a typical area in terms of green land changes and its thermal property effects on the central coastal region of China.
Currently, the understanding of the spatiotemporal heterogeneity of land changes as a result of greening in the central coastal region of China is still insufficient. In particular, a comparison of green land changes at different urbanization levels and the cooling temperature effect of the different green land densities are lacking in this region, which restricts the exploration of the land surface thermal property effect of green land changes. Therefore, the aims of this study are to effectively monitor spatiotemporal changes in green land and analyze their impact on surface thermal properties. Specifically, the objectives are to (1) provide effective monitoring of the spatiotemporal heterogeneity of green land and capture its new spatiotemporal features from 2000 to 2022; (2) reveal the different levels of green land at different stages of urban development and provide the specific quantitative values; and (3) analyze the impact of green land density and its thermal protective effect on the residential environments and map the spatial distribution characteristics of uncomfortably cold and hot areas for urban residents. Finally, three aspects of are discussed, namely, a very green and livable region on the central coast of China is discussed, the different green land changes in different regions of China and worldwide are compared, and the positive effects of dense greening of land on the comfort of urban residents are investigated.
2. Materials and Methods
2.1. Study Area
Rizhao is located in the mid-latitude region of the Earth’s northern hemisphere, with a total area of 5330.72 km
2 in 2022, a longitude range of 118°25′~119°39′E, and a latitude range of 35°04′~36°04′N (
Figure 1). Regarding the natural environment of the study area, the terrain exhibits a high numerical value in the middle and a flat surrounding area. This region has a temperate monsoon climate with a small temperature difference among the four seasons, meaning that the temperature in summer and winter is moderate, and the area holds abundant water resources. From the perspective of administrative division, a branch of socio-economic environments, Rizhao city contains over 30 administrative sub-regions.
Rizhao is an ecologically coastal, livable, tourism-based, modern port and the portside industrial base in China. It is also known as the “Capital of Water Sports” and “Oriental Sun City”. Rizhao won the United Nations Habitat Award for its beautiful environment. Rizhao was also awarded the titles of National Civilized City, National Forest City, National Environmental Protection Key City, and National Ecological Demonstration Zone Construction Pilot City. The pleasant ecology, livable environment, and high coverage of green land make Rizhao a suitable area for exploring the spatiotemporal features of green land changes and its effect on land thermal properties from the perspective of land surface radiation energy balance on the central coast of China.
2.2. Data Collection and Preprocessing
2.2.1. Collection of Land Use Data
Land use data were very important for this study. Before the start of the study, we reviewed the datasets from Tsinghua University, from Wuhan University, from the Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, from the Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, etc. All these data had good time series with high data consistency to use for studying land use and the monitoring of land use changes, but the data formats that were publicly available on the internet were grids [
39,
40,
41,
42]. Perhaps, the vector format of land use data was more convenient for calculation and statistics. As members of the land dataset group at the Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, we obtained the vector land dataset for the years of 2000 and 2015. To include the most recent land use change monitoring data, we extended the land use time to 2022. Then, the land data for the years of 2000 and 2022 were used in this study.
The land use/cover dataset from the Chinese Academy of Science was developed using the process of human–computer interactive interpretation of remotely sensed information to interpret the Landsat digital images covering China, and to construct the national land-use database [
43,
44]. For areas not covered by Landsat image data or covered with poor-quality data, supplemental data from the CCD multispectral data from the Huanjing-1 satellite (HJ-1) were used. The land changes were labeled with attribute codes using professional geoscience knowledge, and the codes can simultaneously reflect the land use type of the dynamic plots during the different stages. The uniform quality control and integration checking for each dataset was used to ensure high-quality and consistent interpretation. The overall accuracies of datasets were above 91% [
26,
44] and meet the requirement of the user mapping. According to the land use classification system from the Chinese Academy of Sciences data, the first-level land types included 6 categories (i.e., cropland, forest land, grassland, water land, construction land, and unused land), and the second-level types included 25 categories. Rich land classification systems can monitor the land cover and its change more effectively [
42,
45].
2.2.2. Remote Sensing Image Download
The remote sensing images from 2022 were downloaded and preprocessed in this section (
Section 2.2.2) and served as the basic data for producing land use maps for 2022. The United States Geological Survey (USGS) provided the free, universal, and better spectral characteristics of land resource satellites, and this was also the main data source for the land use data production process of the Chinese Academy of Sciences. To match the consistency between this image data source and the Chinese Academy of Sciences’ land data, we used the remote sensing satellite from the USGS to obtain land use maps for 2022. For the satellite data quality inspection, Landsat OLI images with poor quality, bad pixels, and cloud cover were removed. This means that we focused on high-quality images covering the study area. For the time of image acquisition, we selected the summer of 2022, as the different types of land surface covers were relatively easy to identify in summer. Then, images with available observations and of high quality were downloaded for 2022, with the images’ path and row number being 120,035. False color band synthesis was performed on these Landsat images. The data used in this study are displayed below (
Table 1).
2.3. The Technical Process of This Study
Main technical process of this study is shown in
Figure 2. This technical process mainly included four steps. The objectives of steps 1 to 4 were to obtain the updated land use in 2022, the classified land covers within the built-up areas, and the calculated land surface temperature as well as the land surface thermal property indicators, and to perform the comprehensive analysis, respectively. Each step is explained in detail below.
The first step was to obtain the latest land use mapping for the year of 2022, so that we could extract the built-up areas in the years of 2000 and 2022. We obtained land use data for the years of 2000 and 2015 from the Chinese Academy of Sciences. Based on the 2015 land use data, the traditional human–computer interactive interpretation process of remotely sensed information to interpret the Landsat digital images was applied in order to establish new land use maps for 2022 on the ArcGIS software platform. After producing the land use data, the accuracy evaluation was conducted using the method of layered random sampling through creating random points functions in the data management tools. Then, the land use data for the years of 2000 and 2022 were used to extract the built-up areas in the years of 2000 and 2022 for this study.
The second step was to obtain the greening of land within the built-up area, so that the spatiotemporal patterns of green land within the built-up area in different urbanization regions (i.e., the old and new urban regions) could be compared in order to delve deeper into the differences in land changes. For generating land cover within built-up areas, the vegetation–impervious surface area–soil model (V-I-S) from the ENVI software platform was used to obtain the optimal endmembers of the different types of land covers. In this model, the principal component transformation first concentrated the spectra of surface objects in the first three bands. A 2D scatter plot and band choice were used to combine the first three short waves in pairs. In these different combinations of windows, the land type features presented different scatter points. The land type spectrums of vegetation, impervious surface area, and soil can be obtained using the region of interest tool to select the corner points of these scatter points. These endmembers were input in the least squares mixed pixel decomposition model to generate the densities of high albedo objects, low albedo objects, vegetation objects, and bare soil objects. These four densities, as well as other indexes, were used to generate the land cover of impervious surface areas, vegetation, bare soil, and water bodies, using the decision tree classifier method and unsupervised classification method. The vegetation cover was extracted using the select tool in the extract module on the ArcGIS platform. By using spatial overlay analysis tools, we superimposed urban built-up areas with vegetation data, and further analyzed the differences in area and proportion changes of vegetation between original and expanded urban built-up areas.
The third step was to obtain the land surface temperature retrieval and various thermal comfort indicators. The split window algorithm was used to generate the surface temperature pattern. The combination of the surface radiation energy balance model, the pixel component arranging and component algorithm model (PCACA), and temperature vegetation dryness index model (TVDI) was then used to obtain various thermal comfort indicators, such as the sensible heat flux, latent heat flux, and available flux on the ArcGIS and ENVI software platforms.
The fourth step was the comprehensive analysis to provide effective monitoring of the spatiotemporal heterogeneity of green land space and capture its new change characteristics from 2000 to 2022, reveal the different green land levels in different regions of urban development and provide their specific quantitative values, and to reveal the cooling temperature effect within different green land densities as well as the thermal comfort indicators.
2.4. Methodology of Digitizing the Land Use Map, Classifying the Green Land Density and Its Cover, and Retrieving the Land Surface Thermal Properties
2.4.1. Producing the New Land Use Map to Include the New Built-Up Area during 2000–2022
For generating the land use maps for 2022 in this study, we used the traditional process of human–computer interactive interpretation, which was the same method as was used for the data production process of 2000 and 2015 by the Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences. Using the land use map of 2015 as the basic data source, the spatial land map of 2015 was superimposed onto the corresponding images from 2022 on the ArcGIS software platform. In the 2015 land use map, a 2022 attribute table field was added, which was used to obtain the land cover types in the year of 2022. Then, the human–computer interactive interpretation method of remotely sensed information was used to interpret the Landsat digital images using professional knowledge by identifying the color and texture of different land features on the Landsat images. The identified land type results were placed in the attribute table’s 2022 field. After identifying the entire research area, a 2022 land spatial map was generated. To ensure the accuracy of the 2022 data, the land use dynamic patches during 2015–2022 were obtained through the calculation of the attribute table’s 2015 and 2022 fields, and different remote sensing professionals conducted multiple reviews of these dynamic patches. Finally, a spatial land use map for 2022 was generated. Then, the different built-up areas in the years of 2000 and 2022 were extracted and calculated using the spatial statistics analysis tools on the ArcGIS software platform.
After producing the land use data for 2022, the layered random sampling scheme was used to achieve accurate evaluation with a total of 300 sampling points for 2022 (
Figure 3) in the data management tools on the ArcGIS software platform. For the high-resolution Google satellite images, the data were not freely accessible in China. We obtained them on the paid and professional 91 bitmap platform. The images contained the following information: resolution: 0.5 m, date: from 1 July to 30 September 2016 and 2022, sensor: QuickBird, vendor: DigitalGlobe and Google. Then, the ground truth distribution samples, producer’s accuracy, user’s accuracy, total classified pixels, and other indicators of per land class were calculated in the accuracy evaluation matrix (
Table 2). The overall obtained land use accuracy was 92.67% in 2022. This means that the land use data achieved good accuracy.
2.4.2. Green Land Density Generation
Mapping of the density of green land can be used to comprehensively analyze the pattern change in the green land component and also to provide basic land data for investigating the cooling temperature effect from greening of land. To obtain different land type density maps (i.e., the area ratio of a land use type ranged from 0.01% to 100% within each grid pixel), the image band fusion technology was first used to upgrade the resolution of the satellite images from 30 m to 15 m [
47] so that the land terrain pattern could be clearly depicted [
48] on the ENVI software platform. After the band fusion, the minimum noise fraction rotation approach and the principal component analysis technology concentrated the spectral information on ground objects on the first three main bands. Through the V-I-S model, we obtained the endmembers in the land surface types of high and low albedo objects, greening land objects, and soil objects. Then, the endmembers were continuously optimized to obtain different land surface types of pure pixels, since the endmembers of pure pixels can obtain high-density maps of high and low albedo objects, greening land objects, and soil objects. In this process, we also used 0.5 m resolution Google satellite images to obtain different interpretation samples from some sample areas to assist in the density decomposition. After obtaining the endmembers of pure pixels, the least squares mixed pixel decomposition model was applied to obtain the density maps of different surface types (i.e., high and low albedo objects, greening land objects, and soil objects), accompanied by corresponding density values of 0.01–100.00% at each grid pixel after the continuous end element optimization using the ENVI software platform.
In this study, the least squares mixed pixel decomposition model was applied, because it can effectively perform linear decomposition of different land use types with the advantage that the DN values of abundance maps for each end element within each pixel range from 0 to 1. The principle of the model is described as follows:
where
is the albedo in the
i-th pixel from band λ,
is the proportion of the area occupied by
k components in the
i-th pixel,
is the albedo in the
k components from band λ, and
is the residual value.
2.4.3. Green Land Cover Classification
Green land classification was used to analyze the distribution of green land changes within the built-up area, further providing the differences in green land levels within different urbanization regions. To obtain green land classification data, we used a combination of supervised classification (i.e., decision tree) and unsupervised classification to elaborate the urban interior land classification map based on the different surface density data obtained in
Section 2.4.2 using the ENVI software platform, thereby extracting the spatial distribution characteristics of green land. In the process, the improved normalized water body index was first input into the decision tree classification, with a threshold greater than 0 to obtain water coverage. Then, the high and low albedo object densities were used to generate impermeable surface areas and greening land density, normalized difference vegetative index (NDVI) was used to obtain vegetation coverage areas, and soil density was used to obtain soil coverage, resulting in the corresponding thresholds of 0.21, 0.18, 0.17, and 0.30 from impermeable surface areas, green land density, NDVI, and soil density, respectively. There were still some small areas that could not be classified by the remote sensing spectrum. Unsupervised classification divided these small mixed pixel distribution areas into 200 categories. Manual interpretation identified the land cover type that each classification belonged to. Then, the classification of green land within urban areas was composed of two-parts, including the decision tree results and the small patches of unsupervised classification results.
2.4.4. Calculation of Land-Surface Temperature
Considering the differences in the meteorological and environmental conditions during the different years of 2000 and 2022 and further avoiding systematic temperature errors, we used only remote sensing images from the summer of 2022 as an example to calculate the land surface temperature as well as to retrieve the land surface thermal property indicators. Land surface radiation energy is typically used to assess the effect of cold and heat on residents. For the surface energy distribution, the spatial distribution of the surface temperature should be presented first [
36]. The split window algorithm model was used to produce spatiotemporal differences in temperature patterns, mainly based on the thermal infrared band of Landsat images and the related parameters in the python module on the ArcGIS platform. The main principles of the model are as follows [
44]:
where, in Equation (2),
is the temperature of the land surface,
and
are the brightness temperatures obtained from Landsat images, and
,
, and
are the transition parameters that can be calculated from Equations (3)–(7). For Equations (3)–(7),
and
represent the surface emissivity and atmospheric transmittance, which is also obtained from the calculation process of Landsat images.
2.4.5. Retrieval of Land Surface Thermal Property Indicators
On the basis of the calculated spatial surface temperature data, the surface radiation energy balance model was then used to obtain variously spatialized thermal comfort indicators, including the sensible heat flux, latent heat flux, and available flux. In this calculation process, some procedural parameters such as downwelling shortwave/longwave and upwelling longwave radiation were also calculated, using the parameter model to obtain the different flux data. The sensible heat flux was used to identify heat generation regions to explore its impact on residents, latent heat flux was used to explore suitable regions for residents’ leisure activities, and available flux was used to evaluate the total energy that residents obtained on the ArcGIS and ENVI software platforms. The main principles of the model are as follows [
45]:
where, in Equations (8) and (9),
is the land albedo from Landsat band2, band4, band5, band6, and band7; and
,
, and
represent the indicators of shortwave radiation, downwelling longwave radiation and upwelling longwave radiation from the Landsat calculation parameters, respectively. In Equation (10), S, L, and Gsoil are the sensible heat flux, latent heat flux, and soil heat flux. These indicators were obtained from Equations (11)–(15). For Equations (11)–(15), the temperature vegetation dryness index (TVDI) was used as a process parameter for calculating vegetation cover and land surface temperature based on optical and thermal infrared remote sensing channel data in Landsat images. The SAFR and LARF are the proportional parameters that are calculated from the pixel component arranging and component algorithm model (PCACA) using the land surface temperature from Landsat images and NDVI from the near-infrared band and infrared band. In Equation (16), AF is the available flux, which can be calculated using the sensible heat flux and latent heat flux.
4. Discussion
4.1. A Highly Green and Livable Region on the Central Coast of China
To explore the spatiotemporal heterogeneity of green land change on the central coast of China, a typical region (i.e., Rizhao region) was investigated in this study. The first important finding was that the study area was a highly green region in China. We found a high coverage of green land areas, and evident differences in urban greening rates among different urbanization regions were revealed. In 2000, the greening rate was 25.34% in the built-up area, and it increased to 42.97% in 2022, with a percentage increment of 17.63%. This means that the proportion of green spaces in urban development has significantly increased. The green land coverage rate in the study area was amongst the highest in China [
49]. Urban design has paid more attention to the role of green spaces [
49,
50], which are good for ecosystem service function for residents, for example, by reducing land thermal properties, absorbing carbon dioxide, and releasing oxygen [
51,
52].
A high green vegetation rate holds a significant application value. From the monitoring of this study, the urban greening rate reached 46.43% in the new urban area (i.e., the urban expansion area from 2000 to 2022). This means that nearly half of the new urban area was covered by green land/. According to our field investigation, in the newly built-up area, the green land space included centralized and scattered park green spaces, community green spaces, road shade areas, parking lot shaded areas, three-dimensional configuration of grassland landscapes, the green of squares, trellises and flower racks, green boundary regions, greening around sports grounds, and green roofs. A high coverage of green spaces also attracted tourists from other places. The rich and beautiful green space design not only provides local residents with a comfortable feeling, but also attracts large numbers of foreign tourists. Also, it is more advantageous for urban planning departments to carry out greener urban construction.
4.2. Comparison of the Differentiated Green Land Changes in Different Regions of China and Worldwide
Urban green spaces under different climatic backgrounds were first compared in China. In western China, especially in arid/semi-arid regions of northwest China [
53,
54,
55], large-scale greening of areas was usually carried out in urban areas to provide better ecosystem services for residents. However, constrained by the arid climate, the proportion of green land was always less than 20% in urban regions in northwest China. In contrast, in eastern China, a humid and semi-humid region [
56,
57,
58], the proportion of green land areas within cities was relatively high, such as in Beijing, the capital city of China, where the proportion of green spaces was nearly 30% in built-up areas in the year of 2021. Further, in the coastal regions, the proportion of green spaces in built-up areas has further increased, such as in the study area of this paper, and the value rose to approximately 43% in the year of 2022. Under different climate backgrounds, the proportion of urban land that is green land resulted in huge differences.
Then, for different regions of the world, such as in Asia, green spaces in urban areas were then discussed. Relevant research indicated that the proportion of urban green land to the corresponding built-up area boundary decreased by 13.93% and 7.28% in capital cities of East Asia and South Asia from 2000 to 2015 [
59]. By contrast, increases happened in the Middle East and Central Asia, with the incremental values of 12.61% and 15.36%, respectively. As for Europe, the spatial patterns of the cities were basically stable, with relatively small changes in green land. The increase in the proportion of green land was only 1.57%.
4.3. Possible Positive Effects of Densely Green Land Space on Human Welfare
Differences in green land coverage in different urbanized regions were observed, namely, the urban greening rate was 37.78% in the old urban area of 2022, while it was as high as 46.43% in the newly expanded urban area of 2022, showing that more attention should be given to the construction of urban green spaces during urban expansion. This finding shows that the creation of a high green coverage in newly expanded urban areas may have a positive impact on the living environment and residents’ walking and travelling. First, it may increase aerobic exercise and improve the comfort of physical/mental health. Photosynthesis of forest and grassland in spaces with high rates of greening may promote the local production of a high-concentration “oxygen bar” environment [
60], similar to a business place, but for people to breathe oxygen. When residents walk and rest in such a high-oxygen environment, their cerebral cortex and brain activities are strengthened, which may promote the absorption of oxygen, accelerate blood circulation, invigorate the spirit, eliminate fatigue, and improve immunity [
61,
62]. Second, it improves the comfort of leisure and improves mood [
63,
64]. According to our survey, the service radius required for Rizhao residents to reach parks and the centralized green space was less than 500 m. This was similar to the situation in areas of high greening in developed countries in the Americas and Europe, which featured a large population and high coverage of urban green spaces. Service radii shorter than 500 m from the park green space and the centralized green space were convenient for residents to stroll there and relax [
65]. Most of their transportation methods for reaching these green regions are walking or cycling, which are not only green and environmentally friendly, but also have a mood-enhancing effect. However, if the service radius of park green spaces and centralized green spaces to people is large [
66], such as in regions of Africa and other arid regions, people have to use private cars or public transportation to reach green land areas. This not only generates noise pollution and is prone to traffic congestion, which affects residents’ moods, but also makes the exhaust gas from fuel-powered vehicles more likely to pollute the environment and generate greenhouse gas emissions, which is not conducive to the United Nations Sustainable Development Goals (i.e., SDGs [
67]) target.
4.4. Limitations of This Study and Future Prospects
The impact of large-scale regional climate simulations on cities should be thoroughly investigated. This study focused on the central coastal areas of China, taking Rizhao City as an example, to monitor urban green land changes and analyze their effect on land thermal properties from a surface radiation energy balance perspective. Specifically, the spatiotemporal heterogeneity of green land changes in different urbanized regions was first revealed. The urban greening rate was 37.78% in the old urban area of 2022, while it was as high as 46.43% in the newly expanded urban area of 2022, showing that more attention should be given to the construction of urban green spaces during urban expansion. Also, the cooling land surface temperature effect in the middle- and high-green-land-density regions was first studied, showing temperatures that were 1.05 °C and 2.12 °C higher than those in the low-density region. Thermal comfort indicators in the study area were also evaluated in terms of latent heat flux, sensible heat flux, and total available energy. All these results were regional investigations. As is known, Rizhao is a coastal city, and the sea should necessarily be taken into account at the regional climate scale. As a next step, we will elevate the investigation from the regional scale to large regional scale to comprehensively analyze the urban energy effect of green spaces and blue sea spaces using the Weather Research and Forecasting Model (i.e., WRF). Then, the comprehensive impact of green land and blue sea on the regional urban climate will be revealed.
More element indicators and green land types should be added. In the meanwhile, the cooling temperature effect of the density of green land was investigated using green land density data and land surface temperature data. The cold and hot discomfort zones of residents caused by land surface temperature were analyzed. In the next step of research, the proper comfort/discomfort zones based on more energy balance parameters should be investigated to conduct a more in-depth analysis from the comprehensive perspective of energy balance, such as by examining the factors of air temperature and relative humidity. The combination of the air temperature and the relative humidity can usually provide a comprehensive analysis of the impact of climate change on the comfort of urban living environments. Furthermore, in this study, urban green space was taken as a whole to explore the cooling effect of land surface on temperature. Although we obtained the cooling effect of green spaces at different densities (i.e., the low-, medium-, and high- density green land spaces), the different types of urban green such as shrubs, forested areas, and different types of grasslands may have different cooling effects on the land surface temperature. In future investigations, such a scientific issue will be summarized using the combination of high-resolution satellite images and surface radiation energy field observation.
More field observations on land surface temperature and its thermal properties should be investigated to conduct more in-depth research. Further, for the monitoring of green land areas, this study first analyzed the trend characteristics of green land changes throughout the whole Rizhao city using land use data. However, this land use classification system compiled the built-up area into one layer, which means that the monitoring of the green spaces within the built-up area cannot be carried out. The combination of V-I-S model—FCLS model—decision tree classifier—unsupervised classification was established to estimate green land cover and its density data using remote sensing images. Then, the green spaces pattern was monitored, and its cooling temperature effects were revealed within the built-up area. But the cooling temperature effects outside the built-up area were not investigated due to the relatively coarse 1 km resolution of the MODIS land surface temperature product and the inconsistency in the time of different Landsat images. The meteorological stations in the study were also relatively sparse. Therefore, more field observations on land surface temperature and its thermal properties should be conducted for more in-depth research, in order to reveal the cooling temperature effect of different green land types outside the built-up area in our study area. Finally, in this study, when we conducted land use accuracy assessments in 2022, high-resolution Google images were not free, which may have affected the universality of the data for other studies. In subsequent research, we will try our best to find free high-resolution remote sensing satellite data for land use accuracy evaluation.
5. Conclusions
This study investigated the spatiotemporal heterogeneity of green land changes in different urbanized regions and analyzed their effect on land thermal properties using the comprehensive method of human–computer interactive interpretation method—urban interior mixed pixel model—surface radiation energy balance model to compensate for the lack of research on green land and its environmental effects on the central coast of China. For green land changes, remote sensing monitoring indicated that the green land proportion in old urban areas was 25.34% in 2000, and that this proportion increased to 37.78% in 2022. The increase in old urban areas means an improvement in the green land service function for residents. We also found that this proportion was as high as 46.43% in newly expanded urban areas, thus providing better visual and environmental comfort. The land thermal property indicators were also assessed, with values of 0–457.83 W/m2 for latent heat flux, 0–645.09 W/m2 for sensible heat flux, and 254.07–659.42 W/m2 for total available energy in the study area. Furthermore, the cooling temperature effect from green land with different densities was revealed, with the lowered temperature values in the middle- and high-green-land-density regions being 1.05 °C and 2.12 °C higher than those of the low-density region. The minimum temperature was 15.85 °C and was distributed mainly in areas where rivers were concentrated. The maximum temperature was 36.64 °C and concentrated mainly in the polymerization effect area of buildings. Also, the spatial patterns of uncomfortably hot and cold areas were depicted. These results/findings provide a reference for the study of green land patterns and their impact on the land surface thermal properties of human settlements in coastal areas of China. This study also investigated the spatiotemporal heterogeneity of green land changes in different urbanized regions (i.e., the old and new urban regions) and the cooling temperature effect of the different green land densities in the central coast of China.
The contribution of this paper was to first compare the different greening levels in the old and newly expanded urban areas in the period of 2000–2022 and to provide the latest heterogeneity of green land changes on the central coast of China. The novelty was that the cooling land surface temperature effect of different green land densities (i.e., the low-, medium-, and high-density green space areas) was first revealed in the central coastal areas of China, and the relevant indicators of surface radiation energy balance were evaluated.