This article is a quantitative study of the urban expansion of Beijing in the past two decades and its impact on the city’s heat island effect. The overall idea of this paper is a ‘basic fact description to phenomenon and law reveal-basic cause analysis-mechanism analysis-model simulation’. In this paper, the effects of urbanisation on warming in Beijing are deduced by nonlinear fitting method. Based on Google Earth Engine remote sensing image data within Beijing Inner Sixth Ring Road and using ArcGIS to retrieve Beijing’s surface temperature, the heat island effect intensity index is calculated. This paper quantitatively analyses the temporal and spatial distribution characteristics and development trends of Beijing’s urban heat island effect, combined with land use cover change (LUCC), Land Surface Temperature (LST), Normalised difference Vegetation Index (NDVI), Normalised difference Building Index (NDBI), which can explore the impact of surface vegetation distribution and building density on the urban heat island effect.

#### Keywords

- simulation model
- remote sensing image
- NDVI
- NDBI

#### MSC 2010

- 03E72

Land use cover change (LUCC) has an important impact on the surface-atmosphere system. Urban expansion is one of the most prominent hotspots in China’s current research on LUCC temporal and spatial pattern changes. And the heat island effect caused by urban expansion is also a typical representative of regional climate changes caused by land-use changes (Lu et al., 2020). The current global urbanisation process is becoming more and more intense. While the process of urbanisation brings economic and social benefits to mankind, it also produces a series of ecological and environmental problems. The urban heat island effect is one of the most important problems. With the rapid development of China’s economy, Beijing, the capital of politics, culture, and economy, has seen significant urban expansion and steady population growth, so the thermal environment is bound to be affected (Cui et al., 2016). This paper is a quantitative study of the extent of urban expansion in Beijing in the past two decades and its impact on the city’s heat island effect.

The overall methodology of this paper is following a description of basic facts, the revelation of phenomena and laws, basic cause analysis, mechanism analysis and model simulation. This paper utilises Bigmaper Beijing Sixth Ring Road remote sensing image data as the basic data, with ArcGIS to retrieve the surface temperature of Beijing. The heat island effect intensity index is also calculated. Then it can be used to quantitatively analyse the temporal and spatial distribution characteristics and development trend of the urban heat island effect in Beijing, combining with some factors, including land surface temperature (LST), normalised difference vegetation index (NDVI) and normalised difference building index (NDBI), to explore the impact of surface vegetation distribution and building density on the urban heat island effect. The results show that from 1999 to 2014, the area of urban construction land in Beijing continued to grow, and the growth rate was relatively fast from 2000 to 2003. The rate of urban expansion slowed down significantly from 2004 to 2012. Based on the core functional area, we can say the city gradually expands to the surroundings, showing that the external development is centred on the main urban area. The temperature of urban land gradually started to dominate the temperature level of the whole region, and the difference between the temperature of urban land and the average temperature of the region has decreased significantly. The urban heat island effect of Beijing has gradually increased from 2014 to 2019, and the spread of the heat island is mainly the extension of the urban new development zone and the transportation network. The urban surface temperature has an obvious negative correlation with the normalised vegetation index and an obvious positive correlation with the normalised building index. Further, the correlation between the surface temperature and the normalised vegetation index is particularly prominent. It is recommended that in the process of urban development and construction, increasing urban greening investment, rationally increasing urban vegetation coverage and moderately reducing urban building density can alleviate the urban heat island effect.

In this study, the model that analyse the impact of urban expansion on the urban environment can be determined using literature. From a quantitative point of view, the urban expansion process and basic characteristics of climate change in Beijing for the past 20 years can be analysed through the establishment of multi-period LUCC and meteorological data. Then, the remote sensing data is used to analyse the spatial distribution pattern of urban heat islands and the underlying surface, by mainly focusing on the correlation analysis of LST, NDVI and NDBI. In this paper, we have proposed to focus on analysing the relationship between urban building density distribution and the urban heat island effect. This paper takes Beijing as the research object. Based on high-resolution remote sensing images and through artificial visual interpretation, combined with the surface temperature data obtained from remote sensing inversion, we correlate the internal relationship between the urban building density distribution and the urban heat island effect. This paper discusses the impact of urban building density distribution on the urban heat island effect. It can provide a scientific basis for the planning and layout of urban construction and urban thermal environmental management.

Encroachment of urban lands causes constrained evaporation which arises from decreased thermal inertia and the vegetation index, and thus ultimately leading to reduced heat loss by latent heat flux; it is the flux of heat from the surface to the atmosphere that is associated with evaporation of water and subsequent condensation of water vapour in the troposphere (Goswami & Singh, 2008). UHI is further aggravated by reduction in speed of the environmental winds, congested transportation networks, increased energy demands, higher anthropogenic heat release and other sky-view factors (Chang Cao et al., 2016). Since urban development, occasioned by increased migration into the metropolitan areas, is relatively difficult to control, the only good urban form is the only way which can help sustain the metropolitan area’s environmental condition.

Beijing is one of the leading growth drivers in the country getting some of the country’s huge investments. Besides the manufacturing, processing and manufacturing industries, Beijing’s growth also draws from the real estate segment that forms an integral economic zone, whose growth started from the mid-1990s. The city forms part of the coastal economic zone even though it is geographically located on the North China Plain (Chang Cao et al., 2016). The Ring and Radial Highway System stands as the city’s ideal transportation mechanism model. The system has received a great structural boost since its creation in the 1950s through new legislations like the 1982 and 1993 comprehensive plans for Beijing’s spatial structure. The intercity highway connecting the city’s 14 satellite towns is known as the 6th ring, while the 4th ring road acts as the edge of the city centre and 5th ring road as the linkage between 10 other scattered districts within the metropolitan area and the lower the number of rings, the closer it is to the centre, which also represents major buildings in these areas. Following the adoption of economic reform policies, China’s land use/cover has undergone radical changes that also accelerated economic growth.

As reported by Qiao et al. (2014), from 1989 to 2010, Beijing’s urban land grew by 775.82 km^{2} at a rate of 184.31%. The mounting intensity and increasing UHI index grew at a high pace from 1989 to 2000 but at a declining pace from 2000 to 2010. The study also reports that ‘the LSTs and urban heat island ratio index (URIs) of urban land in Beijing increased between 1989 and 2010, heat island areas expanded, and the UHI resulting from urban expansion increased’ (Qiao et al., 2014).

China’s land legalities and policies have changed rapidly consistent with the changes in the social and economic conditions. As discussed by Liu, Dunford, Song and Chen (2016), most of these policies are ‘reflected in Hukou reform, urban and rural rights, the Annual Land Use Quota system, public assets, and the local government finance vehicles.’

Committed to Habitat II in 1996 principle, the government or china implemented radical measures associated urban and rural housing. The goal of the policy was to ensure all households have self-contained rooms and to achieve a 9-metre per-capita construction space (Assembly, 2015). The policies have led to an increase in the overall level of urban housing in China. The yearly and persistent large-scale housing construction has reached an unparalleled level in China such that houses constructed increased from 1.22 billion in 1996 to 1.93 billion in 2013. During the same time housing space grew from 395 million to 1.07 billion m^{2} in urban centres (Assembly, 2015).

Internationally recognised l-use regulations such as ‘minimum lot sizes, minimum parking requirements, maximum floor-to-area ratios or floor space index (FSI)’ have had a significant impact on UHI in Chinese cities such as Beijing (Menon et al., 2019). For instance, while FSI limits reduce the ‘density of buildings”, they cannot limit the ‘density of people’ since most individuals opt for smaller and especially informal spaces. Prohibitions on lower height structures inflated prices causing urban sprawl such as the case of Beijing, which caused a 12% expansion of its city boundaries (Deng & Huang, 2004).

China’s land use/cover changes have led to the subdivision of surface structures into three layers based on the specific area’s land use and climate (Siddique et al., 2020). Vegetation cover often acts as a sieve that purifies the air from dust particles and heating thereby improving evapotranspiration which in turn promotes precipitation. In Beijing, after the air passes through the green belt, it is intercepted by the city’s built-up areas and experience solar radiation and absorption, and as a result greenhouse gas emission occurs at the cost of anthropogenic disturbances which makes it warmer, hence urban heat island (Yuan, 2020).

With the intensification of urbanisation, the building area continues to expand, and human activity has increased significantly. These factors may lead to an increase in the intensity of the heat island.

Jing (2019) studied the current situation of the heat island effect in the main urban area of Beijing’s Fifth Ring Road in the summer of 2017, and selected 65 urban green spaces of different sizes and types within Beijing’s Fifth Ring Road through a large amount of data linear fitting and controlled variable method on the influencing factors of the internal temperature of urban green space and the quantitative and influencing factors of the green space cooling range. The study suggests that there is a relatively obvious heat island effect in the main urban area of Beijing. The proportion of high-temperature areas from the Second Ring to the Fifth Ring gradually decreases, and the proportion of low-temperature areas gradually increases. The southwest and southeast quadrants of the study area have the highest proportion of high temperature, and the northeast and the northwest quadrant have the highest proportion of low-temperature areas. Large areas of high-temperature heat island areas are mostly old bungalow residential areas, factories, wholesale markets, train stations, etc. The analysis of the factors affecting the cooling range of green space concludes that the smaller the green area, the smaller the cooling range. When the green area is between 0 ha and 7 ha or the perimeter of the green area is between 0 km and 1.5 km, the cooling range of the green area will follow the increase of green area increases. The percentage of impervious surface in the surrounding land of urban green space also has a significant impact on the cooling range of urban green space. In the case of the same green area, whenever the percentage of the impervious surface area surrounding the green area increases by 1%, the green area cooling range will be decreased by 4.01 m.

Li et al. (2008) used the temperature data from the meteorological observatory in Beijing from 1990 to October 2004 to analyse the characteristics of the urban heat island in autumn in Beijing for the past 15 years. The results showed that the urban heat island in the autumn night in Beijing is stronger than that in the daytime. Besides, a comparative analysis of the characteristics of a strong heat island, a weak heat island and their meteorological influence factors show that there is a strong heat island under certain conditions at night in Beijing in autumn. The formation and maintenance of a strong heat island is the result of the combined effects of multiple factors. During the day and sunny nights, the surface wind field in the suburbs of Beijing is very weak and the atmospheric wind field <47 m in the vertical direction of the urban area continues to be very strong. The rate and amplitude of the atmospheric temperature drop in the suburbs after sunset are much greater than that in the urban area, which promotes the formation and maintenance of a strong heat island at night. The weakening of the urban atmospheric stability and the disappearance of the urban atmospheric inversion are the main reasons for the weakening of strong heat island and its final disappearance at night.

Based on discussing the concept and characteristics of urbanisation, Zhao et al. (2014) analysed the interaction between urbanisation and climate change. Urbanisation has caused significant effects on different climatic elements, but related research is mainly concentrated in terms of local microclimate changes. Climate change is a very complex process.

Assuming the experimental data (_{i},y_{i}

Therefore, the nonlinear fitting method should be used to solve the impact of the urbanisation process in Beijing on the urban warming. The following steps should be followed: First, the scatter diagram should be defined and the function category should be analysed; Second, the initial value of the undetermined parameters is obtained and the best parameters are calculated by using MATLAB software. Third, the determination coefficient is used to compare the effect. The formula of determination coefficient is as follows:
^{2} is to 1, the higher the actual fitting effect is. In MATLAB software to achieve the command of the coefficient of determination:

If the coefficients are polynomial, it is called polynomial regression and the parameters are polynomial coefficients. But if it is of the following forms namely exponential function, logarithmic function, power function, trigonometric function, etc., then it is called nonlinear fitting. Equations with an S-shaped curve include the Rogersti model:
^{−kx}), Richards model: ^{1/δ}, and Weibull model: ^{δ}

In order to better realise nonlinear fitting, the function should be defined first: the definition function of inline can be used for numerical calculation of curve fitting. On the one hand, the M file should be built first, and on the other hand,

After combining the factors that influence the urbanisation process of Beijing on the urban temperature increase, the scatter plot is constructed by combining the above nonlinear fitting method, and an appropriate nonlinear regression equation can be selected, as shown below:
_{0} + _{1}_{2}^{2} and _{1} = _{2} = ^{2}, to obtain _{0} + _{1}_{1} + _{2}_{2}. In this case, x_{1} and x_{2} can be regarded as independent variables, and y can be regarded as a linear function of the two variables. On the other hand, it cannot be directly converted into a linear model to calculate parameters, such that _{1} ln(1 + _{2}

In this case,

When calculating the least square method of data fitting, the concept of this theory should be clarified first. Assume that is a set of data given (_{i}_{i}_{i}

And the calculation function can be obtained:

More in line with the following requirements:

The Φ in the

Assuming

Make

And define the inner product:

The system can be called:

Then the normal equations of the function system at discrete points are as follows:

The determinant formed according to the basis function _{r}_{r}_{0}),_{r}_{1}),...,_{r}_{m}

Therefore, the solution of the system is

And

In the case that the basis is ^{j}

When calculating the data fitting method of least square integration, it is necessary to make clear its concrete concept. First, it is clear that n groups of observed data are (_{i}_{i}

Wherein, G (x) refers to the error function of function _{i}_{i}_{i}_{i}_{i}_{j}

Assuming that

Then in the fitting problem of the traditional least square method, the following formula can reach the minimum value if all the parameters contained therein are calculated according to the above formula. This method is called the least square sum method.

By taking the observed data (_{i}_{i}

In general, the fitting function G (x) represents:

Where _{i}_{1},_{2},...,_{i}_{1},_{2},...,_{q}^{dix}

Here, we first analyse the specific situation when g (x) is used as a linear fitting function, and then we can get:

Where, _{i}

Calculate the minimum value, which is

Assuming that

The

In data fitting analysis, two of the most common fitting function is: _{i}_{i}_{i} (x) refers to the weight function of 1 integral inner product of orthogonal polynomials, and as a result of the (x) function hi conform to the requirements of the

The corresponding regression coefficient is:

On the other hand, the basis function is an exponential fitting function. First observe the data (_{i}

In addition, the above coefficient Ai is defined. At the same time, the root _{1},_{2},...,_{m}

According to different characteristic root analysis, clear
_{r}_{+k},_{1} = 1 so we can get:

And _{1},_{2},...,_{m}

Assuming _{i}_{1},_{2},...,_{m}_{i}_{n}_{+1} as the complex conjugate roots of the equation, and we can get:

The corresponding basis function is:

Assuming that there are k real roots and t complex conjugate roots in the characteristic equation of the fitting curve, then the new fitting function will be:

The in the formula needs to be calculated with the least square integration method proposed above.

According to the above steps, we can not only master more influencing factors, but also put forward better solutions based on the numerical changes obtained from calculation and analysis.

The idea of this paper follows the path of problem discovery, subject research and practical reflection, in which both qualitative and quantitative analyses are applied. Figure 1 shows the workflow of this main methodology. One of the main concerns of qualitative researchers is their emphasis on the description of the urbanisation process, and the research analysis model. In this study, literature research can give the key influencing factors of the urban warming effect. Quantitative research is guided by a positivist methodology to explore the different indexes of urbanisation and the heat island effect. At the same time, the model that analyses the impact of urban expansion on the urban environment can be determined from the literature. From a quantitative point of view, the urban expansion process and basic characteristics of climate change in Beijing can be analysed for the past 20 years through the establishment of multi-period LUCC and meteorological data. Then the remote sensing data is used to analyse the spatial distribution pattern of urban heat islands and the underlying surface, mainly focusing on the correlation analysis of LST, NDVI and NDBI. Through this paper, we are hopeful to focus on analysing the relationship between urban building density distribution and the urban heat island effect. This paper takes Beijing as the research object; using high-resolution remote sensing images and through artificial visual interpretation in combination with the surface temperature data obtained from remote sensing inversion, correlation of the internal relationship between the urban building density distribution and the urban heat island effect is obtained. This paper discusses the impact of urban building density distribution on the urban heat island effect. It can provide a scientific basis for the planning and layout of urban construction and urban thermal environmental management.

The measurement system in this paper uses WGS 1984 and the projection uses equal area projection. The landsat5 and landsat8 data are established through the geospatial data cloud platform, including five dates and specific satellites. The cloud coverage is <10%.

Database

999.08.10 | Landsat 5 |

2005.07.25 | Landsat 5 |

2010.08.08 | Landsat 5 |

2014.08.19 | Landsat 8 |

2019.08.17 | Landsat 8 |

Through ENVI, LST and LUCC are established by data processing, sample establishment, classifier selection and post-processing.

Urban heat island remains one of the biggest attention-drawing environmental problems in China (Zheng & Li, 2016). The data of Beijing districts and counties can be downloaded from the National Basic Information Centre. The datasets used for the study of urban heat island, as an impact of urban expansion, can be obtained through periodic data collection at three different times to get a comparative basis for the changes that occur over a given period. Landsat TM images are then chosen and corrected radiantly to retrieve the land cover types followed by visual interpretation to differentiate false-colour composites of TM/ETM.

Beijing city covers approximately 16410 km^{2}, longitudinally and latitudinally lying at about 2”05′ and 1”37′, respectively (Zhi et al., 2014). The city is graced by a sub-humid warm temperate continental monsoon climate that gives rise to four distinct seasons, humid, hot summer, cold and windy winter. Beijing has continued to experience rapid urbanisation, which has consequently led to increased environmental problems like UHI, pollution haze and dust and sandstorms in the urban areas (UA) (Siddiue et al., 2020). The study area is mapped into specific distinctive portions as UA and non-urban areas (NUA). The calculation of urban area land can be performed using various ways depending on the precision and accuracy of the method, using satellite images and change detection. In this study special nets with specific dimensional measurements that match the Moderate-resolution Imaging Spectroradiometer are created for the calculation of the land area percentage.

NDVI is a common parameter in exploring the influencing factors of ground temperature. It can be used to detect vegetation coverage, growth status and seasonal changes. The formula is as follows,

In the formula, NIR is the near-infrared band reflectivity, R is the red-light band reflectivity, which is represented by the 5th and 4th bands respectively in Landsat 8 TIRS. The value of NDVI ranges from −1 to 1. NDVI>0 means all vegetation. The higher the NDVI, the higher the vegetation coverage, which means good growth and the rich types of vegetation in the area.

LST is a good indicator of the energy balance of the earth’s surface and the greenhouse effect. It is a key factor in the physical process of ground strikes on a regional and global scale. It can reflect the energy flow and material exchange of the soil-vegetation-atmosphere system. It is very necessary for many fields such as climate, hydrology, ecology and biogeochemistry. The ground temperature retrieved from satellite data can reflect in more detail the average temperature of the underlying surface of each pixel and can reflect the spatial distribution characteristics of the underlying surface temperature field. Although the error of using remote sensing data to retrieve the surface temperature is inevitable, it is still the most effective and easiest way to obtain the surface temperature of a large area at present. This is also an advantage of LANDSAT images to reflect the difference of the surface temperature distribution in a large area.

Landsat bands can be used in land detection with the wavelength of 433 nm to 2300 nm and the spatial resolution is 30 m. LST product images, with particular desired properties and resolution levels, are used for analytical purposes. For instance, the Landsat product, which is capable of averaging the LST after every 16 days can be used to obtain the daily LST produced by the split-window algorithm, which is the method for LST retrieval from satellite data. On the other hand, the Landsat Vegetation Index product can be used to obtain data for the EVI data by determining the vegetation index for a given period at a specific spatial resolution. The mean EVI, which is commonly used to determine the vegetation abundance, is calculated for the whole period of the study i.e. 1999 to 2019. The NDBI data is generated using a specific Landsat Surface Reflectance product at a particular spatial resolution. This determination provides the estimated surface spectral reflectance of the study bands, which is then corrected for atmospheric conditions. The NDBI is defined by the equation:

NDBI can be used to characterise the density of urban land?where RSWIR is the surface spectral reflectance of the second last band and RNIR is the surface spectral reflectance of the second band used in the study process.

Using the difference in averaged LST value between the urban and rural areas of study, followed by a determination of the urban area’s pixel levels as the difference in the averaged nonurban LST and each pixel’s LST in the urban area, the SUHII levels of the study areas can be determined by calculation. Investigation of the data obtained from temporal trends of variables like SUHII, EVI, LST and NDBI from 1999 to 2019 is performed at every pixel by the use of linear regression models and ordinary least squares (Minghong et al., 2011). Using the LST as the dependent variable and the study year as the independent variable, multiple ordinary least squares models are run at each pixel to examine the regression coefficient that is used in measuring the annual change in the LSTs magnitude and trend. The significance is then tested at a specific level e.g. at 0.05. Thus, any positive trend with p <0.05 indicates a significantly increased LST at that specific pixel from 1999 to 2019.

Spatiotemporal variation of UHI using LST data using remote sensing images formed at different periods of the study is the major focus of this method. Quantification of the UHI is then performed by introducing URI based on bright temperature normalisation (Siddiue et al., 2020). However, the surrounding conditions of the rural areas affect the magnitude of a UHI and the fact that UHI focuses on the LST’s relative spatial intensity. Normalisation of LSTs is done to allow comparison of their spatial distribution to provide a spatiotemporal pattern variation of the HUI because the remote sensing images only change LST values at different study periods as opposed to changing the spatial LST distribution. The normalised LSTs are then classified into five different thermodynamic levels by density segmentation method to characterise their distribution levels and the area of each level calculated. URI is then introduced for the quantification of the rate of distribution of the urban land to UHI using the following equation:

Where: T_{S} is the surface temperature; B (T_{S}) is the blackbody at the same temperature. K1 and K2 can be obtained from the image header file.

Based on Landsat TM data The URI often reflects the degree of UHI development in built-up land; thus, the more extensive the URI, the more severe the UHI effects and vice versa (Jia Wang et al., 2019). Computational analysis of the UHI area to urban land ratio while considering each temperature level’s weighted values provides the URI value.

Classification of heat island intensity

High temperature | T |

Sub high-temperature | u+0.5std < T < u+std |

Mid temperature range | u-0.5std < T < u+0.5std |

Sub low temperature | u-std < T < u-0.5std |

Low temperature | < T < u-std |

u is the average temperature of surface; std is the standard difference of temperature.

According to the Figure 2, the population of Beijing has been rising linearly for the past two decades, from 9,866,000 to 20,035,000, but in terms of growth rate, it can be seen that 2010 is a cut-off point, 4.82%, after 2010. The average increase rate is around 2.1%. The birth rate is 8.12‰, the death rate is 5.49‰ and the natural growth rate is 2.63‰. The permanent population density is 1,312 people per square kilometre.

This study mainly uses six Districts in Beijing as the research area, and the three phases of LandsatTM remote sensing data are utilised in 1999, 2010 and 2019, combined with GIS technology. This aims to use a comprehensive analysis index system of land use dynamic changes, and reveal the land use of the study area in the past 20 years. The magnitude, speed and direction of change are analysed based on the driving force of land-use change. The results show that the urbanisation of Beijing is very rapid. In the past 20 years, due to the expansion of urban population, economic development and the impact of increasing land use demand has led to a sharp decrease in arable land, while construction land, including residential areas, has increased significantly. The built-up land has increased from 38% in 1999 to 72% in 2019. Due to the dry climate, low annual rainfall, and the interference of human activities on the environment, the water area has become slightly smaller. And it is located in UA. The small part of woodland and grassland on the edge has not changed much due to the influence of topography. The woodland is decreasing after 2010. Comparing with the population trend, there is no correlation between population and LUCC.

LST has a significant correlation with NDBI and NDVI. It can be seen from Figures 3–5 that there are two stages of urbanisation in Beijing, while 2010 is the crest.

LST acts as the interface temperature between the atmosphere and the underlying surface of the city as solar radiation. The analysis of spatial pattern characteristics is of great significance for studying the urban thermal environment, microclimate and optimising the spatial configuration of urban ecological land. The LST data is obtained by ETM image inversion of the study area. The spatial analysis results show that the high-temperature patches in the study area are generally distributed in a pattern of more south and less north, and about 80% of the area is located south of the east-west central line, and the high-temperature patches have low landscape connectivity and tend to be regular in shape. Low-temperature patches are concentrated on the boundary of the study area, with relatively high connectivity and their shape is related to the shape of green spaces and water bodies. In terms of the number and area of ??patches, low-temperature patches are dominant in the second ring and outside the fourth ring, and the second ring to the fourth ring road is dominated by hot patches. The water inside the fourth ring road plays a leading role in the generation of low-temperature patches, and outside the fourth ring road is dominated by low-temperature patches with the forest as the underlying surface.

Land Use Classification

^{2} |
||||||
---|---|---|---|---|---|---|

Cropland | 47.44% | 31.13% | 12.97% | 1075.81 | 705.97 | 294.04 |

Woodland | 7.40% | 8.84% | 7.55% | 167.71 | 200.50 | 171.27 |

Grassland | 6.07% | 2.81% | 5.67% | 137.56 | 63.61 | 128.64 |

Water body | 1.07% | 1.12% | 1.67% | 24.20 | 25.40 | 37.87 |

Built-up land | 38.03% | 56.10% | 72.14% | 862.31 | 1272.11 | 1635.79 |

Beijing’s intense heat island centres are mainly distributed in the core areas and expansion areas of the city. The closer to the city centre, the stronger the heat islands. The thermal field variation index of Dongcheng and Xicheng in the centre is more than ten times that of Tongzhou. The heat island effect in the southern part of the entire study area is stronger than that in the northern part. The average surface temperature in the centre of the heat island in different regions is very different from the average surrounding ground temperature. On the whole, the surface temperature values ??in all directions decrease as the distance from the centre increases, and the surface temperature decreases significantly along the south profile line. In the range of 1/3 of the distance from the centre, the LST value changes relatively smoothly, and the latter fluctuates greatly with the increase of the distance, and the closer to the centre, the more frequent LST high values ??appear.

NDBI and NDVI can give a clear explanation and relationship between building intensity and vegetation intensity. NDBI shows that from 1999 to 2014, the building index increased from the range −0.35 to 0, which means the Building intensity was higher. Initially, only the downtown building index was high, while the edge area was low, but in 2014, both the downtown and surrounding building index were high.

The NDVI can accurately reflect the vegetation coverage on the ground. At present, the NDVI time-series data obtained from satellite remote sensing images such as SPOT/VEGETATION and Landsat have been obtained from the research of vegetation dynamic change monitoring, land use/cover change detection, macro-vegetation coverage classification and net primary productivity estimation in various scale areas. Wide range of applications. It can be seen from 1999 to 2010 that the increase in the red area (NDVI from 0.01 to 0.20) means that the vegetation coverage is reduced. After 2010, it can be seen that the area with NDVI >0.50 is getting larger and larger, indicating that the vegetation coverage rate is good.

The results of the study reveal the spatial distribution of the correlation between buildings and temperature in a better manner, thus provide a certain scientific basis for Beijing’s future city and green space planning.

The Beijing urban development developed slowly prior to the reforms and opening up. After the reform and opening up, the growth rate of the Beijing urban development increased by dozens or even a hundred times than prior to the reforms and opening up. From 1999 to 2010, the high-temperature area expands outward from the urban centre, and by 2010, >80% of the area reaches moderate temperature. From 1999 to 2019, the high-temperature area shifted from the middle to the southwest.

Comparing the UHII with the LUCC and LST results, it is observed that the increase in temperature caused by human activities in Beijing has exceeded the increase in temperature caused by natural factors, and thus it can cause the temperature changes in the city.

Since 1999 in Beijing, the increase in urban temperature has been consistent with the development of urbanisation. They have a long-term coordinated and balanced relationship. This helps the planner to formulate urban planning to reduce the urban heat island effect and a reasonable population-scale city development strategy. Based on high-resolution remote sensing images, this paper extracts the distribution of different building density areas in Beijing’s Sixth Ring Road. Through the analysis of the relationship between the distribution of different building density areas and the surface temperature, it is found that Beijing’s urban building density is getting higher year by year, vegetation density distribution, and there is a certain correlation between the urban thermal environment.

Both the building index and the plant index of the study area present a distribution pattern with a high centre and a low edge. However, from 1999 to 2019, the overall trend of centralisation has weakened. The building density of the entire area is getting higher and higher, and the plant index is showing a downward trend of decentralisation. LST has a significant correlation with NDBI and NDVI. In the south and southwest of the city, there are more areas >310 K. This is related to the acceleration of these trends in urbanisation. The higher the building index, the higher the temperature and the building has a positive effect on temperature. Vegetation was affected by the increase in building density, and its density was getting lower and lower, and the value was around 0.20 to 0.40 by 2019.

Affected by the key development of the southern area of ??Beijing, the southern area has an advantageous geographical location. It is connected to the central city of Beijing in the north, the Xiong’an New District in the south and the sub-centre of Beijing city and the Yizhuang Economic Development Zone in the east. Therefore, it can be seen from UHII that the high-temperature area shifted from the central part of 1999 to the south in 2019.

In the past 20 years, Beijing’s urban construction land has gradually expanded to the surrounding area based on the original core functional areas. It is an extensional development cantered on the main urban area. The construction land in the surrounding districts and counties is also based on the original urban construction land. The above shows an extensional expansion to the surroundings. With the progress of urbanisation, the heat island effect has always been also obvious. The temperature of urban land gradually dominates the temperature level of the whole area, and the difference between the temperature of urban land and the regional average temperature has been significantly reduced.

The development of urbanisation is also reflected in many aspects such as changes in the internal industrial structure of cities, changes in population composition and so on (Cui et al., 2016). There are still seasonal changes and day-night differences in the urban thermal environment under the influence. The impact of urbanisation on the regional surface thermal environment needs to be more comprehensive and in-depth research.

The spatial resolution is relatively insufficient. The morphological and structural characteristics of urban heat islands are the basis for the study of urban heat island effects. The spatial resolution of TM/E TM thermal infrared bands of 60 m or 30 m makes it difficult to observe the detailed features of the ground surface on a smaller spatial scale (Zhou et al., 2019). Moreover, the lack of spatial resolution will inevitably cause the problem of mixed pixels, that is, there will be multiple urban feature components in one pixel, in this way, whether it is in the inversion of the true surface temperature or in the morphological structure and formation of the urban heat island. Moreover, there will be certain errors in mechanism of research.

Because the acquisition of remote sensing images is greatly affected by weather conditions, the weather conditions at the time of image acquisition may have a certain impact on the accuracy of the model (Khanal et al., 2020). At the same time, this article only studies the Beijing area. Therefore, whether the relationship between urban heat island and vegetation greenness studied in this paper applies to other areas requires more data to verify.

The goal is not to study the urban heat island effect, but fundamentally to eliminate and weaken the heat island effect. Therefore, in future research, there should be a focus on the use of Landsat data to study the simulation and prediction of the urban heat island process, the ecological effects and regulation of the urban heat island effect. The study aims to seek mitigation and governance measures for the urban heat island on this basis.

#### Database

999.08.10 | Landsat 5 |

2005.07.25 | Landsat 5 |

2010.08.08 | Landsat 5 |

2014.08.19 | Landsat 8 |

2019.08.17 | Landsat 8 |

#### Classification of heat island intensity

High temperature | T |

Sub high-temperature | u+0.5std < T < u+std |

Mid temperature range | u-0.5std < T < u+0.5std |

Sub low temperature | u-std < T < u-0.5std |

Low temperature | < T < u-std |

#### Land Use Classification

^{2} |
||||||
---|---|---|---|---|---|---|

Cropland | 47.44% | 31.13% | 12.97% | 1075.81 | 705.97 | 294.04 |

Woodland | 7.40% | 8.84% | 7.55% | 167.71 | 200.50 | 171.27 |

Grassland | 6.07% | 2.81% | 5.67% | 137.56 | 63.61 | 128.64 |

Water body | 1.07% | 1.12% | 1.67% | 24.20 | 25.40 | 37.87 |

Built-up land | 38.03% | 56.10% | 72.14% | 862.31 | 1272.11 | 1635.79 |

Law of interest rate changes in financial markets based on the differential equation model of liquidity Basalt fibre continuous reinforcement composite pavement reinforcement design based on finite element model Industrial transfer and regional economy coordination based on multiple regression model Satisfactory consistency judgement and inconsistency adjustment of linguistic judgement matrix Spatial–temporal graph neural network based on node attention A contrastive study on the production of double vowels in Mandarin Research of cascade averaging control in hydraulic equilibrium regulation of heating pipe network Mathematical analysis of civil litigation and empirical research of corporate governance Health monitoring of Bridges based on multifractal theory Health status diagnosis of the bridges based on multi-fractal de-trend fluctuation analysis Performance evaluation of college laboratories based on fusion of decision tree and BP neural network Application and risk assessment of the energy performance contracting model in energy conservation of public buildings Sensitivity analysis of design parameters of envelope enclosure performance in the dry-hot and dry-cold areas Analysis of the relationship between industrial agglomeration and regional economic growth based on the multi-objective optimisation model Constraint effect of enterprise productivity based on constrained form variational computing The impact of urban expansion in Beijing and Metropolitan Area urban heat Island from 1999 to 2019 TOPSIS missile target selection method supported by the posterior probability of target recognition Ultrasonic wave promoting ice melt in ice storage tank based on polynomial fitting calculation model The incentive contract of subject librarians in university library under the non-linear task importance Application of Fuzzy Mathematics Calculation in Quantitative Evaluation of Students’ Performance of Basketball Jump Shot Visual error correction of continuous aerobics action images based on graph difference function Application of Higher Order Ordinary Differential Equation Model in Financial Investment Stock Price Forecast Application of Forced Modulation Function Mathematical Model in the Characteristic Research of Reflective Intensity Fibre Sensors Radioactive source search problem and optimisation model based on meta-heuristic algorithm Research on a method of completeness index based on complex model Fake online review recognition algorithm and optimisation research based on deep learning Research on the sustainable development and renewal of Macao inner harbour under the background of digitisation Support design of main retracement passage in fully mechanised coal mining face based on numerical simulation Study on the crushing mechanism and parameters of the two-flow crusher Interaction design of financial insurance products under the Era of AIoT Modeling the pathway of breast cancer in the Middle East Corporate social responsibility fulfilment, product-market competition and debt risk: Evidence from China ARMA analysis of the green innovation technology of core enterprises under the ecosystem – Time series data Reconstruction of multimodal aesthetic critical discourse analysis framework Image design and interaction technology based on Fourier inverse transform What does students’ experience of e-portfolios suggest Research on China interregional industrial transformation slowdown and influencing factors of industrial transformation based on numerical simulation The medical health venture capital network community structure, information dissemination and the cognitive proximity Data mining of Chain convenience stores location The optimal model of employment and entrepreneurship models in colleges and universities based on probability theory and statistics A generative design method of building layout generated by path Parameter Id of Metal Hi-pressure State Equation Analysis of the causes of the influence of the industrial economy on the social economy based on multiple linear regression equation Research of neural network for weld penetration control P-Matrix Reasoning and Information Intelligent Mining Intelligent Recommendation System for English Vocabulary Learning – Based on Crowdsensing Regarding new wave distributions of the non-linear integro-partial Ito differential and fifth-order integrable equations Research on predictive control of students’ performance in PE classes based on the mathematical model of multiple linear regression equation Beam control method for multi-array antennas based on improved genetic algorithm The influence of X fuzzy mathematical method on basketball tactics scoring Application of regression function model based on panel data in bank resource allocation financial risk management Research on aerobics training posture motion capture based on mathematical similarity matching statistical analysis Application of Sobolev-Volterra projection and finite element numerical analysis of integral differential equations in modern art design Influence of displacement ventilation on the distribution of pollutant concentrations in livestock housing Research on motion capture of dance training pose based on statistical analysis of mathematical similarity matching Application of data mining in basketball statistics Application of B-theory for numerical method of functional differential equations in the analysis of fair value in financial accounting Badminton players’ trajectory under numerical calculation method Research on the influence of fuzzy mathematics simulation model in the development of Wushu market Study on audio-visual family restoration of children with mental disorders based on the mathematical model of fuzzy comprehensive evaluation of differential equation Difference-in-differences test for micro effect of technological finance cooperation pilot in China Application of multi-attribute decision-making methods based on normal random variables in supply chain risk management Exploration on the collaborative relationship between government, industry, and university from the perspective of collaborative innovation The impact of financial repression on manufacturing upgrade based on fractional Fourier transform and probability AtanK-A New SVM Kernel for Classification Validity and reliability analysis of the Chinese version of planned happenstance career inventory based on mathematical statistics Visual positioning system for marine industrial robot assembly based on complex variable function Mechanical behaviour of continuous girder bridge with corrugated steel webs constructed by RW Research on the influencing factors of agricultural product purchase willingness in social e-commerce situation Study of a linear-physical-programming-based approach for web service selection under uncertain service quality A mathematical model of plasmid-carried antibiotic resistance transmission in two types of cells Fractional Differential Equations in the Exploration of Geological and Mineral Construction Burnout of front-line city administrative law-enforcing personnel in new urban development areas: An empirical research in China The Law of Large Numbers in Children's Education Data structure simulation for the reform of the teaching process of university computer courses Calculating university education model based on finite element fractional differential equations and macro-control analysis Educational research on mathematics differential equation to simulate the model of children's mental health prevention and control system Analysis of enterprise management technology and innovation based on multilinear regression model Verifying the validity of the whole person model of mental health education activities in colleges based on differential equation RETRACTION NOTE Innovations to Attribute Reduction of Covering Decision System Based on Conditional Information Entropy Research on the mining of ideological and political knowledge elements in college courses based on the combination of LDA model and Apriori algorithm Adoption of deep learning Markov model combined with copula function in portfolio risk measurement Good congruences on weakly U-abundant semigroups Research on the processing method of multi-source heterogeneous data in the intelligent agriculture cloud platform Mathematical simulation analysis of optimal detection of shot-putters’ best path Internal control index and enterprise growth: An empirical study of Chinese listed-companies in the automobile manufacturing industry Determination of the minimum distance between vibration source and fibre under existing optical vibration signals: a study Nonlinear differential equations based on the B-S-M model in the pricing of derivatives in financial markets Nonlinear Differential Equations in the Teaching Model of Educational Informatisation Fed-UserPro: A user profile construction method based on federated learning Smart Communities to Reduce Earthquake Damage: A Case Study in Xinheyuan, China Response Model of Teachers’ Psychological Education in Colleges and Universities Based on Nonlinear Finite Element Equations Institutional investor company social responsibility report and company performance Mathematical analysis of China's birth rate and research on the urgency of deepening the reform of art education Precision Machining Technology of Jewelry on CNC Machine Tool Based on Mathematical Modeling First-principles calculations of magnetic and mechanical properties of Fe-based nanocrystalline alloy Fe _{80}Si_{10}Nb_{6}B_{2}Cu_{2}Computer Vision Communication Technology in Mathematical Modeling The Effect of Children’s Innovative Education Courses Based on Fractional Differential Equations Fractional Differential Equations in the Standard Construction Model of the Educational Application of the Internet of Things Optimization in Mathematics Modeling and Processing of New Type Silicate Glass Ceramics Has the belt and road initiative boosted the resident consumption in cities along the domestic route? – evidence from credit card consumption MCM of Student’s Physical Health Based on Mathematical Cone Attitude control for the rigid spacecraft with the improved extended state observer Sports health quantification method and system implementation based on multiple thermal physiology simulation Research on visual optimization design of machine–machine interface for mechanical industrial equipment based on nonlinear partial equations Research on identifying psychological health problems of college students by logistic regression model based on data mining Abnormal Behavior of Fractional Differential Equations in Processing Computer Big Data Mathematical Modeling Thoughts and Methods Based on Fractional Differential Equations in Teaching A mathematical model of PCNN for image fusion with non-sampled contourlet transform Nonlinear Differential Equations in Computer-Aided Modeling of Big Data Technology The Uniqueness of Solutions of Fractional Differential Equations in University Mathematics Teaching Based on the Principle of Compression Mapping Influence of displacement ventilation on the distribution of pollutant concentrations in livestock housing Cognitive Computational Model Using Machine Learning Algorithm in Artificial Intelligence Environment Application of Higher-Order Ordinary Differential Equation Model in Financial Investment Stock Price Forecast Recognition of Electrical Control System of Flexible Manipulator Based on Transfer Function Estimation Method Automatic Knowledge Integration Method of English Translation Corpus Based on Kmeans Algorithm Real Estate Economic Development Based on Logarithmic Growth Function Model Informatisation of educational reform based on fractional differential equations Financial Crisis Early Warning Model of Listed Companies Based on Fisher Linear Discriminant Analysis Research on the control of quantitative economic management variables under the numerical method based on stochastic ordinary differential equations Network monitoring and processing accuracy of big data acquisition based on mathematical model of fractional differential equation 3D Animation Simulation of Computer Fractal and Fractal Technology Combined with Diamond-Square Algorithm The Summation of Series Based on the Laplace Transformation Method in Mathematics Teaching Optimal Solution of the Fractional Differential Equation to Solve the Bending Performance Test of Corroded Reinforced Concrete Beams under Prestressed Fatigue Load Radial Basis Function Neural Network in Vibration Control of Civil Engineering Structure Optimal Model Combination of Cross-border E-commerce Platform Operation Based on Fractional Differential Equations Research on Stability of Time-delay Force Feedback Teleoperation System Based on Scattering Matrix BIM Building HVAC Energy Saving Technology Based on Fractional Differential Equation Human Resource Management Model of Large Companies Based on Mathematical Statistics Equations Data Forecasting of Air-Conditioning Load in Large Shopping Malls Based on Multiple Nonlinear Regression System dynamics model of output of ball mill Optimisation of Modelling of Finite Element Differential Equations with Modern Art Design Theory Mathematical function data model analysis and synthesis system based on short-term human movement Sensitivity Analysis of the Waterproof Performance of Elastic Rubber Gasket in Shield Tunnel Human gait modelling and tracking based on motion functionalisation Analysis and synthesis of function data of human movement The Control Relationship Between the Enterprise's Electrical Equipment and Mechanical Equipment Based on Graph Theory Financial Accounting Measurement Model Based on Numerical Analysis of Rigid Normal Differential Equation and Rigid Functional Equation Mathematical Modeling and Forecasting of Economic Variables Based on Linear Regression Statistics Design of Morlet wavelet neural network to solve the non-linear influenza disease system Nonlinear Differential Equations in Cross-border E-commerce Controlling Return Rate Differential equation model of financial market stability based on Internet big data 3D Mathematical Modeling Technology in Visualized Aerobics Dance Rehearsal System Children’s cognitive function and mental health based on finite element nonlinear mathematical model Motion about equilibrium points in the Jupiter-Europa system with oblateness Fractional Differential Equations in Electronic Information Models Badminton players’ trajectory under numerical calculation method BIM Engineering Management Oriented to Curve Equation Model Optimal preview repetitive control for impulse-free continuous-time descriptor systems Development of main functional modules for MVB and its application in rail transit Study on the impact of forest fire prevention policy on the health of forest resources Mathematical Method to Construct the Linear Programming of Football Training The Size of Children's Strollers of Different Ages Based on Ergonomic Mathematics Design Stiffness Calculation of Gear Hydraulic System Based on the Modeling of Nonlinear Dynamics Differential Equations in the Progressive Method Relationship Between Enterprise Talent Management and Performance Based on the Structural Equation Model Method Value Creation of Real Estate Company Spin-off Property Service Company Listing Selection by differential mortality rates Digital model creation and image meticulous processing based on variational partial differential equation Dichotomy model based on the finite element differential equation in the educational informatisation teaching reform model Nonlinear Dissipative System Mathematical Equations in the Multi-regression Model of Information-based Teaching The modelling and implementation of the virtual 3D animation scene based on the geometric centre-of-mass algorithm The policy efficiency evaluation of the Beijing–Tianjin–Hebei regional government guidance fund based on the entropy method The transfer of stylised artistic images in eye movement experiments based on fuzzy differential equations Research on behavioural differences in the processing of tenant listing information: An eye-movement experiment A review of the treatment techniques of VOC Some classes of complete permutation polynomials in the form of ( x ^{pm}−x +δ )^{s}+ax ^{pm}+bx overF _{p2m}The consistency method of linguistic information and other four preference information in group decision-making Research on the willingness of Forest Land’s Management Rights transfer under the Beijing Forestry Development A mathematical model of the fractional differential method for structural design dynamics simulation of lower limb force movement step structure based on Sanda movement Fractal structure of magnetic island in tokamak plasma Numerical calculation and study of differential equations of muscle movement velocity based on martial articulation body ligament tension Study on the maximum value of flight distance based on the fractional differential equation for calculating the best path of shot put Sports intensity and energy consumption based on fractional linear regression equation Analysis of the properties of matrix rank and the relationship between matrix rank and matrix operations Study on Establishment and Improvement Strategy of Aviation Equipment Research on Financial Risk Early Warning of Listed Companies Based on Stochastic Effect Mode Characteristics of Mathematical Statistics Model of Student Emotion in College Physical Education Mathematical Calculus Modeling in Improving the Teaching Performance of Shot Put Application of Nonlinear Differential Equation in Electric Automation Control System Nonlinear strategic human resource management based on organisational mathematical model Higher Mathematics Teaching Curriculum Model Based on Lagrangian Mathematical Model Optimization of Color Matching Technology in Cultural Industry by Fractional Differential Equations The Marketing of Cross-border E-commerce Enterprises in Foreign Trade Based on the Statistics of Mathematical Probability Theory The Evolution Model of Regional Tourism Economic Development Difference Based on Spatial Variation Function The Inner Relationship between Students' Psychological Factors and Physical Exercise Based on Structural Equation Model (SEM) Fractional Differential Equations in Sports Training in Universities Higher Education Agglomeration Promoting Innovation and Entrepreneurship Based on Spatial Dubin Model