Contrasting Trends of Forest Coverage between the Inland and Coastal Urban Groups of China over the Past Decades

: China is building forest urban groups through reforestation and a ﬀ orestation. However, the fast process of urbanization inevitably conﬂicts with multiple vegetated areas around cities. Hence, it is critical to evaluate the changes in regional vegetation cover and its spatial pattern due to complex natural and anthropogenic factors. Nevertheless, systematic studies to quantify and compare the development of forest urban agglomerations were rarely reported. Based on a remote sensing landcover dataset from 1992 to 2015, this study investigated forest cover changes and the impacts on landscape pattern in several urban groups, and tried to explore their di ﬀ erences between the inland and coastal regions of China. The results showed that over the past 24 years, the forest coverage in the coastal urban agglomerations declined (103 km 2 / year) while it increased (26 km 2 / year) in the inland urban agglomerations. There was a certain conﬂict between forest and cropland for the coastal urban agglomerations where the forest area converted to cropland accounted for 61.9% of the total forest loss. The increase in forests coverage in inland urban agglomerations mainly came from grassland which nearly accounted for 66.47% of the total increase. The landscape diversity has also changed in areas where forests have changed signiﬁcantly


Introduction
As the largest developing country in the world, China is in the stage of rapid urban development, which has caused a drastic spatio-temporal impact on land-use/land-cover (LULC) around cities [1,2]. In addition, land cover changes caused by urbanization can directly affect the ecological landscape pattern around cities. Studies have shown that different land cover types have different biomass, water, and soil conditions and, consequently, land cover changes can affect community structure and landscape structure in the environment [3,4]. Many studies, based on landscape ecology, have focused on the internal structural systems of cities and surrounding areas in the context of rapid urban land cover change, and evaluated the ecological environment.
The landscape pattern is closely related to various processes in the landscape. The study of it has laid an important theoretical foundation for the study of regional ecological security pattern [5][6][7]. Landscape pattern is a concept in landscape ecology, which indicates that a series of landscape elements of different sizes and shapes are formed by nature or man-made. It is a complex interaction of various complex physical, biological, and social factors [8][9][10][11][12]. It helps to clarify changes in ecological effects in land-use change. The landscape pattern index is a commonly used method and tool in current research. Many indices, such as the Shannon diversity index, fragmentation index, and aggregation index, are widely used to quantitatively evaluate the relationship between land cover and regional   The Central Plains urban agglomeration (termed CPua) (110 • 9 -118 • 18 E, 32 • 24 -37 • 48 ), the core of the Central Plains economic zone, is located in the central transportation hub of China. It includes 30 prefecture-level cities such as Zhengzhou, Henan Province and Changzhi, Shanxi Province and Handan, Hebei Province and Liaocheng of Shandong Province and Huaibei of Anhui Province. The CPua is typical Eastern monsoon climate with mild climate and distinct seasons. Much of its landform is plain, the terrain is flat, and the main vegetation type is the deciduous broad-leaved forest. The Guanzhong-Tianshui urban agglomeration (termed GTua) (104 • 31 -112 • 48 E, 33 • 09 -37 • 24 N) is located in the Guanzhong Plain area of Shaanxi Province (Figure 1), with Xi'an as the central city and Baoji as the sub-central city, including Tianshui, Weinan, Tongchuan, Longnan, and other secondary cities. It is semi-humid and semi-arid climate. From north to south, loess hills, Guanzhong Plain, and Qinling Mountains go across the whole area. The main vegetation type is the deciduous broad-leaved forest.

ESA CCI Land Cover Data
This paper used the land cover data of ESA CCI production period of 24 years (1992-2015), with a spatial resolution of 300 m. The data came from http://maps.elie.ucl.ac.be/CCI/viewer/download. php [35]. These images divided the land surface into 37 land cover classes based on the United Nations Land Cover Classification System (UN-LCCS). This unique long-term land cover time series was achieved by combining five different global daily surface albedo observation systems with a pixel accuracy of 71.5%. The highest user accuracy values were found for the classes of rainfed cropland, irrigated cropland, broadleaved evergreen forest, urban areas, bare areas, water bodies, and permanent snow and ice [36]. In addition, this paper mainly focused on the proportion of plant functional type (PFT) in qualitative land cover types, so 37 classes of initial ESA CCI data are aggregated into 0.5 × 0.5 resolution and converted into PFT data. The vegetation types of PFT include forest, shrub, farmland, and grassland. Combining with the actual situation in the study area, we calculated the second-level classification under the first-level classification through the transformation function, including 15 s-level categories, as shown in Table 1.

Land Use Transfer Matrix
In our present work, a transfer matrix method was used to investigate the annual change of land cover pattern during the study period. The land cover transfer matrix uses two land cover images from different periods, taking the pixels whose coordinates in both images are (x, y), and the gray values in the two images are, respectively, G i , G j . All possible gray value combinations (G i , G j )|i, j = 1 . . . n| in the image are mapped to the position of the corresponding element S ij in the matrix, then count the number of occurrences of all gray value combinations and assign it to S ij . The transfer matrix of land use type change and its quantitative information are obtained. The mathematical expression of the transfer matrix is as follows: S ij = S 11 S 12 · · · S 1n S 21 S 22 · · · S 2n · · · · · · · · · · · · S 21 S 22 · · · S nn (1) where S ij is the number of land cover type i converted into type j in the formula, and S ij ≥ 0.

Landscape Pattern Index
The Shannon Diversity Index (SHDI) refers to the diversity of landscape elements or ecosystems in terms of structural function and temporal variation. It reflects the complexity of the landscape and the heterogeneity of the landscape, especially sensitive to the uneven distribution of patch types in the landscape.
In order to explain the results, the paper also used some other landscape pattern indices, including the patch density (PD), the landscape shape index (LSI) and the aggregation index (AI). The PD expresses the number of plaques per unit area, which reflects the degree of fragmentation of the landscape to some extent. The AI is calculated based on the common boundary length between patches of the same type of plaque which indicates the non-randomness or degree of aggregation of different plaque types in the landscape. Additionally, the LSI is calculated based on the degree of deviation between the shape of a patch and a circle or square of the same area, reflecting the complexity of the shape of the patch [37,38]. The formula for the four indices is shown in Table 2. i the landcover class p i the proportion of landscape patch types. n i the area of the patch A the total area of landscape. E the boundary length of patches g ii the number of similar adjacent patches for the corresponding landscape type.

Statistical Analysis Method
(1) According to the frequency of forest occurrence in the study area during 1992-2015, the non-forest area, the area where the forest appeared but changed, and the area was always forest for the 24 years were obtained. (2) In order to assess the trend of forest cover change in the four different urban agglomerations from 1992 to 2015, the entire area was investigated with a 17 × 17 sliding window, and the 24-year trend of forest coverage in the window was analyzed by least squares method. Hence, a spatially continuous distribution of trends of forest changes over the study area was obtained. (3) Using the 17 × 17 sliding window over the whole area, the transfer matrix for each contiguous two years was calculated, and the land cover types were converted to PFTs. Similarly, the SHDI was calculated for each window each year. Finally, the calculated forest transfer amount and SHDI were assigned to the central pixel of the window, then the transfer amount between forest and other vegetation types and the spatial distribution of landscape diversity were obtained for further analysis.

Temporal and Spatial Variations in the Forests
From 1992 to 2015, the change of forest coverage in the urban agglomerations was obvious. The statistics (PFT-based) of changing area show that the forest cover of the YRDua and the PRDua decreased by about 1384.44 km 2 and 883.68 km 2 , respectively; the forest cover of the GTua and the CPua increased by about 311 km 2 and 133.1 km 2 , respectively ( Figure 2). two urban agglomerations. The evergreen deciduous forest in the YRDua decreased by 1493.62 km 2 from 1992 to 2015, and its reduction rate (67.064 km 2 /year) was slightly greater than the regional forest reduction rate (60.706 km 2 /year). The evergreen deciduous forest in the Pearl River Delta urban agglomeration decreased by 1217.11 km 2 , and the reduction rate (58 km 2 /year) was about 1.5 times that of the forest reduction rate (42.355 km 2 /year).
The forest area of the two inland urban agglomerations increased over the past 24 years. The area of forest in the GTua increased from 12,077.8 km 2 (8.26%) in 1992 to 12,210.9 km 2 (8.35%) in 2015. It experienced the fastest growth rate between 1994 and 2006, reaching 14.65 km 2 /year. The growth trend decreased in 2005-2009, and the annual area increased by 5.61 km 2 . After that, the forest coverage tended to be stable. The area of forest in the CPua increased from 4957.7 km 2 (2.38%) in 1992 to 5268.7 km 2 (2.53%) in 2015. It has grown in an "S" shape since 1995. The growth rate was the fastest in 1997-2001, with an average annual growth of 48.82 km 2 . Then it slowed down in 2004-2015 and the area remained stable.  The forest coverage of the two urban agglomerations in the coastal areas declined over the past 24 years. In the YRDua, the area of forests decreased from 11,842.69 km 2 (6.15%) in 1992 to 10,454.25 km 2 (5.44%) in 2015. It had the highest rate of decline in 1994-2000 (147.164 km 2 /year) and in 2001-2010, the rate of decline slowed down and the forest area gradually became stable. The area of forests in the PRDua decreased from 112,786.53 km 2 (12.17%) in 1992 to 11,902.85 km 2 (11.33%) in 2015. It has declined since 1994 with an approximate linear decline trend (43.49 km 2 /year) to 2015. The change in the number of evergreen deciduous forests had a significant impact on the total forests of the above two urban agglomerations. The evergreen deciduous forest in the YRDua decreased by 1493.62 km 2 from 1992 to 2015, and its reduction rate (67.064 km 2 /year) was slightly greater than the regional forest reduction rate (60.706 km 2 /year). The evergreen deciduous forest in the Pearl River Delta urban agglomeration decreased by 1217.11 km 2 , and the reduction rate (58 km 2 /year) was about 1.5 times that of the forest reduction rate (42.355 km 2 /year).
The forest area of the two inland urban agglomerations increased over the past 24 years. The area of forest in the GTua increased from 12,077.8 km 2 (8.26%) in 1992 to 12,210.9 km 2 (8.35%) in 2015. It experienced the fastest growth rate between 1994 and 2006, reaching 14.65 km 2 /year. The growth trend decreased in 2005-2009, and the annual area increased by 5.61 km 2 . After that, the forest coverage tended to be stable. The area of forest in the CPua increased from 4957.7 km 2 (2.38%) in 1992 to 5268.7 km 2 (2.53%) in 2015. It has grown in an "S" shape since 1995. The growth rate was the fastest By the statistics of the 24-year land cover data in the study area, we obtained a picture of forest cover in the four urban agglomerations between 1992 and 2015 ( Figure 3). It can be seen that the forest coverage around both inland urban agglomerations and coastal urban agglomerations in the process of urbanization shows an uneven feature in space.
Zhejiang Province was located in the southern part of the YRDua, in which forests mainly grew. The regions where forests changed during 1992-2015 were mainly distributed in the southwestern part of Hangzhou, Zhejiang Province, the western part of Jinhua, and the eastern part of Taizhou. The forest coverage of the PRDua was characterized by high north and low south. The Nanling Mountains and hilly areas located in the northern part, in which the forests mainly grew. The regions where forests changed were mainly in Shaoguan, Qingyuan, Heyuan, and other cities.
The forest coverage of the CPua showed the characteristics of high in the west and low in the east. Wangwu-Taihang Mountain and Funiu Mountain were in the western region where most forests grew. The regions where forests changed during 1992-2015 were mainly in Changzhi Jincheng Jiaozuo. The forest coverage of GTua was characterized by low north and high south. Forests in the urban agglomerations were mainly grew in the southern Qinling Mountains. The regions where the forests changed were mainly in the southern part of Xi'an in Shaanxi Province and in the eastern part of Yuncheng in Shanxi Province.  changed were mainly in the southern part of Xi'an in Shaanxi Province and in the eastern part of Yuncheng in Shanxi Province. Figure 4 shows the spatially continuous distribution of trends of forest changes in the coastal and the inland urban agglomerations between 1992 and 2015. In the YRDua, the area where the forest area had a slight decreasing trend (−0.01-0 of k) and a slight increasing trend (0-0.01 of k), accounted for 51% and 17%, respectively. In terms of its spatial distribution, the regions with a slight downward trend of forest coverage were mainly concentrated in the hilly areas in Central Zhejiang, including Hangzhou, Shaoxing, Ningbo, and other prefecture-level cities, while the regions with a slightly upward trend of forest coverage were mainly concentrated around the Dabie Mountains, including Anqing, Tongzhou, and Chizhou. Moreover, the regions with a significant downward trend of forest coverage (−0.03-−0.01 of k), accounting for 1.4%, mainly concentrated in the southern part of Huzhou in Zhejiang Province and the coastal areas of Taizhou and Ningbo. In the PRua, the forest change was characterized by spatial agglomeration. The regions with a slight decreasing trend (−0.01-0 of k) and a slight upward trend (0-0.01 of k) of forest coverage accounted for 67% and 28%, respectively. In the past 24 years, the areas where forest coverage slightly reduced were mainly concentrated in Qingyuan and Zhaoqing in the northern part of the urban agglomeration; the areas where forest coverage slightly reduced were mainly concentrated in Huizhou and Shanwei in the southeast coast of Guangdong Province. The area with a significant increase trend of forest coverage (0.01-0.03 of k), accounted for about 3%, concentrated in the northern part of Shanwei and the central part of Guangzhou while the area with significant reduction trend of forest coverage (−0.03-−0.01 of k) accounted for about 1% mainly distributed in Qingyuan and parts of Shaoguan.

Forest Transferations
After counting the amount of transfers between the forest and the other types in the urban agglomerations, it can be found that in the YRDua, an average of 6 km 2 of the other vegetation converted into forests each year, while an average of 34 km 2 of forests converted into the other In the GTua, the forest areas with a slight downward trend (−0.01-0 of k) were much more than those with a slight increase trend (0-0.01 of k), at 51% and 36%, respectively. The regions with a slight reduction in forest cover were mainly concentrated in the eastern part of Tianshui and southwest of Baoji. The regions with a slight increase in forest cover were mainly concentrated in the southern part of Xi'an and Yuncheng, Shanxi Province. The regions with a significant increase trend (0.01-0.03 of k) accounted for about 0.3%, mainly distributed in some areas of Yuncheng. In the CPua, the regions with a slight decrease (−0.01-0 of k) and a slight increase of forest coverage (0-0.01 of k) accounted for 27% and 15%, respectively. The regions where forests slightly shrank were mainly distributed in the central part of Changzhi and the western part of Handan, while the majority of those with increasing forest areas were around the Taihang Mountains and Wangwu Mountains. The regions with a significant upward trend (0.01-0.03 of k) accounted for 0.5% of the total, concentrated in the Taihang Mountains, Wangwu Mountains, and other mountainous areas.  Figure 5 shows that between 1992 and 2015, the amount of the vegetation which converted from the forest cover followed the order: cropland > grassland > shrub from large to small, while the cropland is slightly different in the order that the number of the other PFTs transformed into the forest: that in the YRdua followed: farmland > shrub > grassland and that in the PRDua area followed: shrub > grassland > farmland. In the YRDua, the areas where the forest transformed to the other vegetation types were the largest (125.5 km 2 ) between 1994 and 1995, of which the area from the forest to the cropland was 116. accounted for 44.06 km ). The amount of the forest converted to the other vegetation has been small for many years (6.6 km 2 /year), and it peaked at 29.7 km 2 /year in 1994. The area of the other vegetation converted to the forest in the CPua showed the characteristics of first decline and then rise in 1994-2003. It averaged 35 km 2 (grassland accounted for 23.54 km 2 ) per year and peaked in 1998 (79.68 km 2 /year). After 2003, the figure fell to 3.16 km 2 /year. While the area of the forest converted to the other vegetation was less than 10 km 2 in the whole study period (1.76 km 2 /year). If we pay attention to the spatial distribution ( Figure 6) we will find that, in the YRDua, the regions where the forest obviously transformed into the other vegetation (300-1850 m 2 ) were appeared in Qingyuan and Shaoguan, the northern part of Guangdong Province, while other regions where the vegetation transformed into forests (300-1870 m 2 ) were concentrated in Huizhou and Heyuan in the southeast coastal area. In the GTua, the regions with obvious forest transformation were distributed in the southern mountainous area of the Qinling Mountains. In the eastern part of Tianshui, the forest was obviously transformed into the other vegetation (370-780 m 2 ), while in the northern part of Xi'an, the other vegetation transformed obviously into forest (31-1400 m 2 ). In the CPua, the regions where other vegetation converted obviously into the forest (311-1550 m 2 ) were mainly concentrated in the Wangwu-Taihang mountainous area and the Fuxiu mountainous area in Western Henan, while the situation of forest conversion to other vegetation was generally flat (0-260 m 2 ).

Landscape Patterns of Forest
The landscape pattern of the YRDua was characterized by low north and high south. A large number of croplands were distributed in the northern Jiangsu Plain and the Jianghuai Plain in the YRDua. Additionally, because of the single land use type, the landscape diversity was low in the northern part of the urban agglomeration. Between 1992 and 2015, the SHDI of the YRDua changed little, and it showed a downward trend (reduced by 0-0.03) in general. The areas with a significant upward trend (increased by 0.03-0.1) were mainly in central Shanghai and eastern Hangzhou. The overall landscape diversity of the PRDua was high, and the SHDI value in most areas was between 1 and 2. From 1992 to 2015, the areas where landscape diversity increased in the PRDua were more than the reduced area, accounting for 73.8% and 20.8%, respectively. In the past 24 years, the landscape diversity in the PRDua changed significantly. The southern part of Qingyuan and the central part of Shaoguan dropped significantly (reduced by 0.2-0.7), while Foshan, Zhongshan, and other regions increased significantly (increased by 0.2-0.8).
The landscape diversity of the GTua was higher in the middle and lower in the north and south. Between 1992 and 2015, in the GTua, the areas where landscape diversity declined were more than

Landscape Patterns of Forest
The landscape pattern of the YRDua was characterized by low north and high south. A large number of croplands were distributed in the northern Jiangsu Plain and the Jianghuai Plain in the YRDua. Additionally, because of the single land use type, the landscape diversity was low in the northern part of the urban agglomeration. Between 1992 and 2015, the SHDI of the YRDua changed little, and it showed a downward trend (reduced by 0-0.03) in general. The areas with a significant upward trend (increased by 0.03-0.1) were mainly in central Shanghai and eastern Hangzhou. The overall landscape diversity of the PRDua was high, and the SHDI value in most areas was between 1 and 2. From 1992 to 2015, the areas where landscape diversity increased in the PRDua were more than the reduced area, accounting for 73.8% and 20.8%, respectively. In the past 24 years, the landscape diversity in the PRDua changed significantly. The southern part of Qingyuan and the central part of Shaoguan dropped significantly (reduced by 0.2-0.7), while Foshan, Zhongshan, and other regions increased significantly (increased by 0.2-0.8).
The landscape diversity of the GTua was higher in the middle and lower in the north and south. Between 1992 and 2015, in the GTua, the areas where landscape diversity declined were more than the increased area, accounting for 62.5% and 14%, respectively. During the 24 years, the landscape diversity of the GTua changed significantly. The SHDI in northern part of Xi'an, Shaanxi Province and Yuncheng, Shanxi Province increased significantly (increased by 0.2-0.5), while the SHDI in central Xi'an was significantly reduced (reduced by 0.2-1). In the CPua, the spatial distribution of landscape diversity was characterized by high northwest and low southeast. The areas where the SHDI was less than one in the urban agglomerations accounted for 61.7% in 2015, while the area where the SHDI was larger than one accounted for 38%.7%. Between 1992 and 2015, the areas where the landscape diversity increased was less than the reduced area in the CPua, accounting for 36.4% and 57%, respectively.
During the past 24 years, the SHDI in Xuchang, Henan Province and parts of Liaocheng, Shanxi Province increased significantly (increased by 0.2-0.5) while the SHDI in Central Zhengzhou and Southern Jincheng decreased significantly (reduced by 0.2-0.8).
The SHDI reflects the richness of land use types. Overall, the abundance of land use in the coastal urban agglomerations was higher than that of the inland urban agglomerations. Additionally, the SHDI could reach a maximum of 2.4 or more, while the SHDI values of the inland urban agglomerations were all around 2.2.

The Effect of Mixed Pixels
There is a problem of mixed pixels in low-resolution remote sensing data, which leads to inaccurate results of various landcover classification of statistics. In the LULC data, due to spatial resolution limitations, the image corresponds to a larger area of the actual ground in a pixel. When we focus on the area covered by vegetation, we will find that there are many kinds of vegetation types mixed in our research area and other areas, such as cropland and shrubs interlaced in a certain area, closed forest interlaced with shrubs in a certain area, etc. Mixed pixels are most likely to occur in these regions. Although some special land cover classification, such as mosaic cropland (it is believed that more than 50% of these types of pixels are cropland, and the rest are some mixed natural vegetation, including tree, shrub, and herbaceous cover) are classified in the CCI land cover classification system, these types are not much helpful in calculating the specific content of vegetation in mixed pixels. The LULC data we used in the study has a resolution of 1/360 • , and it is difficult to eliminate the influence of mixed pixels.
In order to deal with this effect, we used the PFT conversion factors provided by ESA for further statistics for more accurate vegetation coverage. The conversion factors refer to the spatial proportion of various plant functional types in each land cover classification, which can reflect the spatial proportion of different vegetation types in mixed pixels to a certain extent. It gives the proportion of different PFTs in some complex land cover classification (such as mosaic cropland) [39,40]. Using this method, we can further calculate the amount of forest resources in the study area. However, the conversion factor is the spatial distribution scale coefficient of large-area scale vegetation. In fact, the proportion of plant function types of mixed pixels may change in different spatial regions.
In order to explore the difference between the results between the results derived from the PFT and classification of land cover, we calculated the area of forest cover for different urban agglomerations in the study area based on these two different approaches. According to the land cover data, the coverage area of the 2010 and 2015 Yangtze River Delta urban agglomerations is 12,400 km 2 and 12,300 km 2 , respectively; the forest coverage of the YRDua obtained from the PFT was 10,700 km 2 and 10,600 km 2 , respectively. In the CPua, the forest coverage based on PFT is 5200 km 2 in the both two years. In short, the forest area obtained by calculating PFT is slightly smaller than that obtained by land cover data. When we focus on the trend of forest area, the difference between the two is almost the same, but the PFT results changes more gradual (Figure 7). This is because a pixel with a land cover classification could contained more than one PFT, which may be not consistent with the plant type of the classification. From the results, the two are different but the difference is small. We could not tell which results are more accurate. However, our research focused on the dynamic process of vegetation in a wide range of areas. Hence, the impact of trend changes was small.

Forest Resources in Different Regions
Urbanization will inevitably lead to conflicts between cities and other land resources, especially the forest resources we were concerned about. The increase in urban area also means that the population growth, which will in turn lead to an increase in demand for food. The demand for food from urban development may lead to an increase in cropland and may, therefore, have a potential impact on the spatial distribution of forests. Moreover, the protection of the ecological environment may also affect the spatial allocation of many other resources. From the above, we can see that with the development of the region, the forest resources will face complex and diverse pressures.
China's LULC data for the past 20 years showed that the urbanization rate varies in different urban agglomerations. In the two coastal urban agglomerations, the urbanization rate has gradually declined after the initial rapid growth of the urban area. However, this feature did not appear in the inland urban agglomerations. In the Yangtze River Delta urban agglomeration, the urban area growth rate was 100 km 2 /year between 2013 and 2015, which was about 2.5 times lower than that in 2005-2012. In contrast, the urban areas of the inland urban agglomerations have maintained a relatively high growth rate during the study period ( Figure 8). China's urbanization rate in 2015 is 56.1% [41,42], which was in the middle and late stages of urbanization (50%-70%) [43]. According to the law of urban development summarized by some scholars, the rate of urbanization in China may gradually slowdown in the future [44]. Statistical data from municipal governments also shows that, in 2015, the urbanization rate of the Yangtze River Delta urban agglomeration was 68.2% and 75.0%, respectively, while the urbanization rate of the Central Plains urban agglomeration and the Guanzhong-Tianshui urban agglomeration was 48.3% and 51.5% respectively. The Yangtze River Delta and Pearl River Delta urban agglomerations may have entered the deceleration stage.  From the perspective of the development of the city, if we assume that the urbanization of the more developed coastal urban agglomerations has entered the later stage of urbanization, when the urban area growth slowed down and the demand for other resources tended to be stable. Additionally, the Yangtze River Delta and the Pearl River Delta in the coastal areas are in the subtropical monsoon climate, and the hydrothermal conditions there are conducive to forest growth. In this way, we have reasons to believe that in the future, the demand of coastal urban agglomerations for various resources may gradually stabilize, and the situation of forest reduction will be alleviated.
On the other hand, some regions, such as the two inland urban agglomerations in the study area, may just enter the middle and late stages of urbanization, and are likely to experience a period of rapid urban growth. More importantly, in these areas, croplands occupied a considerable area, so that changes in the cropland resources may greatly affect forests [45,46]. It is difficult to say if the development of the city will have a greater demand for cropland resources and, thus, change the spatial distribution of cropland. In fact, in the course of urban development, the area of croplands in the inland urban agglomerations even declined after 2000 (Figure 9). This change after 2000 may be mainly due to the nation's "Grain for Green" program [47] and the land set-aside program. In such circumstances, whether the cropland would pose a greater threat to the forest is still worth exploring.

Regional Landscapes
Considering the long-term changes in land cover in the urban agglomerations, the impact of human activities on the ecological changes of the urban agglomerations was extremely complicated. Many studies showed that the driving factors of such impact came from social, economic and natural and other aspects. Although quantitative research under the interaction of such complex factors is difficult, the landscape pattern index and its changes can be a good indicator of the results under the influence of multiple factors [48,49]. In order to consider the change of forest landscape in the urban agglomerations, we also used some other landscape indexes to quantitatively describe it.
We selected the patch density (PD), the landscape shape index (LSI), and the aggregation index (AI) at the class level. The results showed that the PD in the two coastal urban agglomerations increased over the years, while that in the inland urban agglomerations showed a declining trend. That is to say, the fragmentation degree of forest landscape in the coastal urban agglomeration area increased, while the inland area decreased in the past 24 years. Moreover, the LSI changed slightly in the different urban agglomeration areas, but there is an upward trend in the coastal urban agglomeration areas and a downward trend in the inland urban agglomeration areas. Additionally, the AI showed a downward trend in the coastal urban agglomeration areas and an upward trend in the inland urban agglomeration areas ( Figure 10). agglomerations, we also used some other landscape indexes to quantitatively describe it.
We selected the patch density (PD), the landscape shape index (LSI), and the aggregation index (AI) at the class level. The results showed that the PD in the two coastal urban agglomerations increased over the years, while that in the inland urban agglomerations showed a declining trend. That is to say, the fragmentation degree of forest landscape in the coastal urban agglomeration area increased, while the inland area decreased in the past 24 years. Moreover, the LSI changed slightly in the different urban agglomeration areas, but there is an upward trend in the coastal urban agglomeration areas and a downward trend in the inland urban agglomeration areas. Additionally, the AI showed a downward trend in the coastal urban agglomeration areas and an upward trend in the inland urban agglomeration areas ( Figure 10). In the coastal urban agglomerations, the forest coverage decreased, but the degree of fragmentation increased and the degree of aggregation decreased. In the inland urban agglomeration In the coastal urban agglomerations, the forest coverage decreased, but the degree of fragmentation increased and the degree of aggregation decreased. In the inland urban agglomeration areas, because of the increasing number of forests, patches may increase or expand outward, resulting in reduced fragmentation and increased aggregation. As for the LSI, although the changes were not so obvious, we could still find that the overall level of the coastal urban agglomerations was higher than that of the inland urban agglomerations, and there was an upward trend. The decrease of forest fragmentation and the increase of forest aggregation indicated that the growth of forests in the inland urban agglomerations was more likely to be a centralized growth forced by human activities, which may be the result of some land policies of local governments such as the Grain for Green Policy. However, the growth of such a model did not make the structure of the forest landscape be more complicated, and the LSI showed a downward trend. However, from another perspective, increasing fragmentation usually means that the forest patches were divided into smaller patches by other land cover classes, which may have a negative impact on the animals living in larger areas of forest. Considering, comprehensively, whether the strategy of development in this region would help to improve regional ecological services and improve biodiversity is still worth exploring.

Conclusions
This research detected the differences of forest coverage between the inland and coastal urban groups of China over the past decades. The forest coverage had changed dramatically over the past 24 years. This study showed that between 1992 and 2015, the forests in the coastal urban agglomerations significantly declined (reduced by 9% compared with 1992), while that in the inland urban agglomerations raised to some extent (increased by 2.6% compared with 1992) between 1992 and 2015. The changing areas of the forest in the coastal urban agglomerations were widely observed in space and varied greatly. In contrast, the forest cover in the inland urban agglomerations changed in few regions with a relatively small magnitude. The regions with obvious changes of forest cover in the Yangtze River Delta and Pearl River Delta urban agglomerations accounted for 1.4% and 4%, respectively. In contrast, that only accounted for 0.3% and 0.5% in the Guanzhong-Tianshui and Central Plains urban agglomerations. In addition, given the spatial distribution of the forest changing trends, the forest coverage in the inland urban agglomeration tended to increase in the northern regions, while that in the coastal urban agglomeration liked to rise in the south. For the transformation of the forest cover, the mutual change between forest and cropland in the two coastal urban agglomerations were very intense, which often dominated the change of regional forest quantity. Nevertheless, grasslands tended to transform into forests in the inland urban agglomerations. The decline of forests in the two coastal urban agglomerations slowed down around 2003, and forest cover remained basically stable after 2003. It may be because that some protection policies for forests have worked. Our finding also showed that the landscape diversity of the coastal urban agglomerations is higher than that of the inland urban agglomerations, but the landscape diversity of inland urban agglomerations had increased over the past 24 years. Given the complexity of the landscape structure, we argue that more careful considerations of the landscape configuration are important for achieving the goal of sustainable development.