Spatial Pattern of Construction Land Distribution in Bays along the Coast of Vietnam

: Current studies of urban spatial morphology have rarely focused on the particularity of coastal cities, such as sea–land gradient features and bay types. In this paper, we provide a method to analyze the spatial and vertical distribution of construction land in the bay area and discuss the inﬂuence of bay type on the distribution characteristics. Landscape indexes and a clustering algorithm were used to identify the spatial pattern of construction land distribution. Strip division was used to analyze the landscape features of construction land in the sea–land direction. We established eight large bays in Vietnam as the study area. According to the formation and geomorphology of the bay, the eight bays were divided into ﬁve types: delta bay, lagoon bay, structural bay, bedrock erosion-stacking bay and estuary bay. The construction land data were generated from Landsat images captured in 1988 and 2015 or so. The study e ﬀ ectively identiﬁed ﬁve spatial patterns of construction land distribution in bays without prior knowledge. Danang Bay (a bedrock erosion-stacking bay) had a typical high intensity-high concentration pattern, and Hai Loc Bay (a delta bay) had a typical high intensity-low concentration pattern. In the sea–land gradient, the proportions of construction land in the 0–2 km zone were greatest. In Danang Bay, the proportion generally decreased in the sea–land direction; while in Hai Loc Bay, the proportion increased.


Introduction
Urbanization is a spatial process. As a result of urbanization, the change of urban spatial form clearly reflects the process of urbanization. The evolution of urban spatial morphology has always been a hot topic in urban geography. A bay city is a special kind of city type, characterized by three characteristics: highly developed, greatly affected by the ocean, and greatly affected by the type of the bay. First, located at the interface between the sea and the land, bays are a high-value, highly developed and highly fragile ecosystem. Since ancient times, bays have been the gateway to overseas transportation and have acted as industrial centers, playing an extremely important strategic role in economic and social development [1,2]. Many of the world's famous ports, such as the Port of Los Angeles, Shenzhen Port and Bangkok Port, have developed in bays. Second, throughout the history of coastal cities [3], the development of coastal cities often depends on the marine economy. The siting and development of urban centers are influenced by the distance from the ocean and coastal geography. Finally, influenced by the formation of the bay, the types of the bay are various, and the topographical morphology varies greatly. It will affect city orientation and city form. Nowadays, as the population continues to converge to the coastal zone, the demand and scale of construction land in bays is expanding, and the built-up area is expanding rapidly. Quantitative analyses of the development The land-sea gradient and elevation differentiation are two of the most basic spatial pattern characteristics of many geographical elements [27]. As the distance from the sea increased, the physiographic features, biogenic features, and physiochemical features shows a tendency of gradient [28]. Accordingly, the land-use patterns vary in the sea-land direction. Strip division was a common method, and it involves parallel divisions moving the coastline inland with consistent spacing. Hou et al. [27] conducted strip segmentation at an interval of 2 km perpendicular to the coastal zone and chose 15 strips on the landward side as the research areas to analyze the land structure characteristics and differentiation characteristics of different strips. The results showed that the 4-6 km strip had a peak value of the comprehensive land use index, which was followed by the 2-4 km strip and 8-10 km strip. Ding et al. [29,30] developed a model, called sequence-based clustering of coastal land use pattern (SCCLUP), to mine the coastal land use sequence patterns (CLUSPs) along the sea-land direction. The research in the major coastal zone of Bohai Bay and the Yellow River Delta showed that the artificial level of CLUSPs was continuously increasing, and new CLUSPs tended to distribute around port areas. To a certain extent, topography has shaped the distribution pattern of land use [31]. Especially in the coastal areas, the geomorphologic differences in different regions are rather significant. Zhang et al. [32] divided the coastal zone surrounding the South China Sea into 51 zones according to the coastal landform and discussed the spatiotemporal characteristics of construction land expansion based on different geomorphologic backgrounds. The results showed that in areas where there were good conditions for water transport and wide hinterlands, construction land on the delta plain coast and estuary plain coast expanded rapidly.
According to the results of a literature search, the typical urban forms come from empirical exploration. Studies are deficient in the quantification of the spatial pattern of construction land without prior knowledge. Furthermore, current studies of urban spatial morphology rarely focused on the particularity of coastal cities, such as sea-land gradient features and bay types. In this paper, we provided a method to analyze the spatial pattern of construction land in the bay area. It included two parts: spatial analysis and vertical analysis. Landscape indexes and a clustering algorithm were used to identify the spatial pattern of construction land distribution. Strip division was used to analyze the landscape features of construction land in the sea-land direction. The bay type is defined according to the factors underlying the bay formation and is used to analyze the influence on the spatial and vertical distribution of construction land. In order to exclude the influence of different national policies as well as taking into account the richness of bay types, we chose eight bays along the coast of Vietnam as the research areas.

Study Area
The study area is located along the coast of Vietnam, roughly between the latitudes of 8 • N and 24 • N and the longitudes of 104 • E and 109 • E ( Figure 1). Vietnam is located on the eastern coast of the Indochina Peninsula and bordered by China to the north, Laos and Cambodia to the west, and the South China Sea to the east. Vietnam can be considered a marine nation. The main characteristics are as follows: (1) Hot and humid climate. Vietnam has a tropical monsoon climate with abundant rainfall and perennially high temperatures. The average annual rainfall is approximately 1680 mm on the Red River Delta, with 1650 mm along the central coast and 1980 mm in the Mekong Delta [33]. The annual average temperature is 24 • C. (2) Long coastline. Vietnam presents an S-shaped pattern in a north-south vertically long strip [34].
Vietnam is long from south to north and narrow from east to west. The long coastline is unusually sinuous with many bays, lagoons, estuaries, mangroves and outlying islands [35]. (3) Fertile plains. Topographically, the terrain is high in the west and low in the east. Three-quarters of the region is covered by mountains. There are two great rivers flowing through the region that form two vast low-lying delta plains: the Red River Delta in the north and the Mekong River Delta in the south. In addition, some small and narrow plains are distributed in the coastal zone. The coastal plains were created by deposits of alluvial soils by rivers [33]. The fertile soil and good climate conditions make Vietnam the second-largest agricultural country in Southeast Asia. (4) Rich resources. Vietnam is rich in coal mineral resources, oil and gas, fisheries resources and port resources. Vietnam currently has 3.5 billion tons of coal reserves, and vast coal mines are located in the Red River Delta and northern Quang Ninh Province [34]. The large nearshore oil-gas basins include the Yinggehai basin and Mekong basin [36]. Vietnam has a long tradition in both inshore and deep-sea fisheries, and aquaculture production rose to the third greatest globally (behind only China and India) in 2010 [35,37]. Topographically, the mountains extend directly into the sea, creating a number of protected harbors, including those of the port cities of Da Nang, Qui Nhon, and Nha Trang [33].
ISPRS Int. J. Geo-Inf. 2020, 9, x FOR PEER REVIEW 4 of 21 The coastal plains were created by deposits of alluvial soils by rivers [33]. The fertile soil and good climate conditions make Vietnam the second-largest agricultural country in Southeast Asia. (4) Rich resources. Vietnam is rich in coal mineral resources, oil and gas, fisheries resources and port resources. Vietnam currently has 3.5 billion tons of coal reserves, and vast coal mines are located in the Red River Delta and northern Quang Ninh Province [34]. The large nearshore oilgas basins include the Yinggehai basin and Mekong basin [36]. Vietnam has a long tradition in both inshore and deep-sea fisheries, and aquaculture production rose to the third greatest globally (behind only China and India) in 2010 [35,37]. Topographically, the mountains extend directly into the sea, creating a number of protected harbors, including those of the port cities of Da Nang, Qui Nhon, and Nha Trang [33]. Bays are sea areas with obvious coastal curvature that extends into the mainland or islands. The line connecting the two headlands at the mouth of the bay marks the boundary between the bay and the sea [38]. According to the factors underlying bay formation, bays can be categorized as primary bays or secondary bays ( Table 1). The primary bay is formed in the depression basin, river valley, crater and other depressions dominated by bedrock. Once a primary bay is formed, its morphology changes very slowly. A secondary bay is the sea area enclosed by the accumulation of oceans, rivers, or organisms. The secondary bays are relatively unstable and evolve much faster than the primary bays. Primary bays can be subdivided into structural bays, bedrock erosion-stacking bays, estuary Bays are sea areas with obvious coastal curvature that extends into the mainland or islands. The line connecting the two headlands at the mouth of the bay marks the boundary between the bay and the sea [38]. According to the factors underlying bay formation, bays can be categorized as primary bays or secondary bays ( Table 1). The primary bay is formed in the depression basin, river valley, crater and other depressions dominated by bedrock. Once a primary bay is formed, its morphology changes very slowly. A secondary bay is the sea area enclosed by the accumulation ISPRS Int. J. Geo-Inf. 2020, 9, 707 5 of 21 of oceans, rivers, or organisms. The secondary bays are relatively unstable and evolve much faster than the primary bays. Primary bays can be subdivided into structural bays, bedrock erosion-stacking bays, estuary bays and crater bays. Secondary bays can be subdivided into lagoon bays, tombolo bays, delta bays and atoll bays [39]. There are many bays along the coast of Vietnam, and the bay types vary dramatically. Through the water area of the bay, eight bays with a water area, over 50 square kilometers were established as the study area (Table 2)

Data Sources
The data used in this study included the Landsat 5 thematic mapper (TM) images and Landsat 8 operational land imager (OLI) images. The images were downloaded from the US Geological Survey Center for Earth Resources Observation and Sciences (USGS, https://earthexplorer.usgs.gov/) at Level 1. TM images were captured in 1988 and 1989. OLI images were captured in 2015-2017. The resolution of the image is 30 m. In total, we collected 9 Landsat 5 TM images and 9 Landsat 8 OLI images. The images are cloudless in most of the bays. In a few bays, the cloud cover is less than 5%. The detailed information of remote sensing images is shown in Table 3. Image clipping and image enhancement were applied before interpretation.

Data Analysis
The research diagram flow was described in Figure 2. First, we acquired the construction land data from Landsat images. Then, we introduced the method of landscape ecology and analyzed the spatial distribution of construction land in two aspects. In the spatial analysis, the K-means clustering algorithm was used to identify the spatial pattern of construction land distribution. In vertical analysis, the method of coastal buffer division was used to analyze the landscape features of construction land in the sea-land direction. Moreover, the influence of bay type on the spatial and vertical distribution of construction land was discussed.

Data Acquisition
The bay area consists of the land area and sea area. The sea area was defined as the waters enclosed by the coastline and the entrance of the bay. For coasts with mangroves, the coastline lies on the landward vegetation line. For coasts of other types, the coastline lies on the edge of the water. The land area was defined as the buffer zone within a certain distance from the coastline. We have adopted a uniform distance for all the bays to ensure a fair assessment. The distribution range of settlement varies in the eight bays. When determining the extent of the land area of the bay, the choice was based on retaining the larger cities along the bay coast. Therefore, the land area was defined as the landward area extending 10 km from the coastline.
Automatic classification and manual visual interpretation were used to generate the land use data. The classification system included eight primary types: cultivated land, forestland, grassland, construction land, industrial land, water, coastal wetland and unutilized land. Overall, fifteen land-use types were identified at the secondary level [41]. Multiscale segmentation and spectral difference segmentation were performed before classification. Six parameters, including the mean value of the red band/thermal infrared band/mid-infrared band, border index, asymmetry and main direction, were used in the maximum-likelihood supervised classification. Based on the classification result, we checked and corrected any apparent classification errors manually.
Together 900 randomly selected reference pixels (200 points for Ganh Rai Bay and 100 points for every other bay) were used to validate the precision. The overall precision and single precision for the construction land are shown in Table 4. The final overall accuracy for each bay was higher than 85%. The construction data were then generated from the land use data.

Landscape Analysis
In our study, on the single class of construction land, we analyzed on class level. The correlation between each landscape index is often very high, and it is necessary to make a reasonable choice according to the actual situation. Cushman et al. [42] used principal component analysis (PCA) to identify independent components of landscape structure and cluster analysis to group the components. Seven class-level landscape structure components were found that were particularly strong (Table 5). They are "aggregation", "large patch dominance", "shape and correlation length of large patches", "patch size variation", "edge/patch density", "mean patch size", and "edge + aggregation". On this basis, we choose six landscape indexes, including aggregation index (AI), largest patch index (LPI), fractal dimension index (FDI), edge density (ED), landscape shape index (LSI), and patch density (PD). The indexes were calculated in FRAGSTATS 4.2 (http://www.umass.edu/landeco/research/ fragstats/fragstats.html). We found that ED was more relevant with PD, and the changes of FDIs in sea-land gradient were irregular. Therefore, we finally used four indexes, including LPI, PD, LSI and AI, to describe the pattern of construction. The explanation of the landscape indexes is shown in Table 6. Table 6. Landscape index and descriptions [43].

LPI
The area (m 2 ) of the largest patch in the landscape is divided by the total landscape area (m 2 ). Size of dominant plaque PD The number of patches of corresponding patch type (class) per unit area. Fragmentation degree of landscape LSI A standardized measure of patch compactness that adjusts for the size of the patch. Shape complexity of the landscape AI The ratio of the observed number of like adjacencies to the maximum possible number of like adjacencies given the proportion of the landscape comprised of each patch type.
Aggregation degree of the homotypic patches

K-means Clustering Algorithm
The K-means clustering algorithm is one of the most widely used partition clustering algorithms and has been proposed for more than 50 years [44]. In the analysis of the spatial pattern of construction land in bays, we classified the bays by K-means clustering. The bays are considered a dataset with four-dimensional data points. It can be described as follows: where n is the number of bays and n equals 8 in our study. The K-means clustering algorithm organizes data objects into K partitions. Each partition represents a class, and each class has a category center. The Euclidean distance was selected as the similarity and distance criterion to calculate the square sum of the distance from each point in the class to the cluster center (J(c k )). The goal of clustering is to minimize the sum of squares of the total distances (J(C)) [45].
where C is the clustering dataset, k is the clustering number, x i is the original data point to be classified, and u k is the category center of class k.

Buffer Division Method
To measure the spatial pattern of construction land in the direction of the coastal zone perpendicular to the coastline, we divided the coastal zone vertically.
The following strategies were used, and the segmentation results are shown in Figure 3.  Table 7 and Figure 4 show the proportions and areas of construction land in each bay. The spatial distributions of construction land in bays for 1988 and 2015 are shown in Figure 5. Across the region, First, land use data were used to generate the coastlines. Second, 2 km, 4 km, 6 km, 8 km and 10 km buffer zones were generated on the landward side, and the boundary lines were obtained.

Regional Trends
Finally, the buffer zone was cut by using the coastline and boundary lines to obtain buffer strips of 0-2 km, 2-4 km, 4-6 km, 6-8 km and 8-10 km on the landward side. Table 7 and Figure 4 show the proportions and areas of construction land in each bay. The spatial distributions of construction land in bays for 1988 and 2015 are shown in Figure 5. Across the region, the total proportion was 4.9% in 1988. The built-up area significantly expanded during the following 30 years. In 2015, the proportion more than doubled, with increases of up to 12.5%. The average annual growth rate was 0.25%.   Among the bays, the proportions of construction land in Hai Loc Bay, Dam Ha Trung Bay and Danang Bay were more than 10% in 1988. These three bays are distributed along the northern coast of Vietnam. In contrast, construction land was distributed sparsely in the five southern bays. The proportions were lower than 5%, especially in Xuan Dai Bay and Rach Giá Bay, where the proportions were lower than 1%. In addition to the above three bays, the proportion of construction land in Ganh Rai Bay climbed to 15.4% in 2015. However, these four bays presented growth that exceeded the average growth. The proportion of construction land was greatest (26.5%) in Danang Bay, followed by Hai Loc Bay (24.5%). The highest growth rate was also in Danang Bay (15.2%), which was twice the average growth (7.6%). It is worth mentioning that although the proportion of construction land in Ganh Rai Bay was not too high (3.5% in 1988 and 15.4% in 2015), the growth rate was drastic (11.9%). The net increase of built-up area in Banh Rai Bay was also the largest, up to 187.4 square kilometers. In addition, the net increase of built-up area in Dam Ha Trung Bay was 110.8 square kilometers, which ranked 2nd by area.   Figure 6 shows the proportions and areas of construction land in each bay. Except for the AIs, the bay-level LPIs, PDs and LSIs in 2015 were greater than those in 1988. This result indicated that the distribution of construction land became more fragmented, and the shape tended to become more complex over 30 years. There was a significant increase in the LPIs of Hai Loc Bay, Dam Ha Trung Bay, Danang Bay and Ganh Rai Bay. In Danang Bay, the LPIs were significantly greater than in the   Figure 6 shows the proportions and areas of construction land in each bay. Except for the AIs, the bay-level LPIs, PDs and LSIs in 2015 were greater than those in 1988. This result indicated that the distribution of construction land became more fragmented, and the shape tended to become more complex over 30 years. There was a significant increase in the LPIs of Hai Loc Bay, Dam Ha Trung Bay, Danang Bay and Ganh Rai Bay. In Danang Bay, the LPIs were significantly greater than in the other bays and increased by more than three times until 2015. In the other four bays, the LPIs were very low with little change. The PDs and LSIs in Hai Loc Bay were much greater than those in other bays in both periods. Especially for the PDs, it increased to 2.0928 in 2015, which was more than 5 times the PDs of the other bays. Moreover, the lowest PDs and LSIs were recorded in Rạch Giá Bay. The changes in the AIs were relatively random. The most obvious change was in Xuan Dai Bay. The AIs increased from 74.6706 to 88.3288, and the ranking increased from 8th to 7th. other bays and increased by more than three times until 2015. In the other four bays, the LPIs were very low with little change. The PDs and LSIs in Hai Loc Bay were much greater than those in other bays in both periods. Especially for the PDs, it increased to 2.0928 in 2015, which was more than 5 times the PDs of the other bays. Moreover, the lowest PDs and LSIs were recorded in Rạch Giá Bay. The changes in the AIs were relatively random. The most obvious change was in Xuan Dai Bay. The AIs increased from 74.6706 to 88.3288, and the ranking increased from 8th to 7th. (3) Medium intensity-medium concentration pattern. Typical features of this pattern were medium LPIs, LSIs and AIs. Dam Ha Trung Bay and Ganh Rai Bay presented this pattern in 2015. The LPIs in the two bays were just lower than those in Danang Bay, and the LSIs were just lower than those in Hai Loc Bay. Among them, Dam Ha Trung Bay was converted from pattern 2 (high intensity-high dispersion pattern) to pattern 3 from 1988 to 2015. The change indicated that Dam Ha Trung Bay experienced moderate expansion, and the expansion tended to be clustered.

Spatial Trends
(4) Low intensity-medium concentration pattern. The PDs, LPIs and LSIs were at a low level, while the AIs were at a medium or high level. Vinh Van Phong Bay, Cam Ranh Bay, Ganh Rai Bay and Rạch Giá Bay presented this pattern in 1988. However, Ganh Rai Bay was converted to pattern 3 (medium intensity-medium concentration pattern) in 2015. This result indicated that Ganh Rai Bay experienced moderate expansion based on the maintenance of a relatively concentrated distribution. In addition, Xuan Dai Bay converted to pattern 4 in 2015.
(5) Low intensity-low concentration pattern. The LPIs, LSIs and AIs were at a low level. Construction land in such bays was sparsely distributed with small patches. Xuan Dai Bay presented this pattern in 1988. The LPI and AI were the lowest among all bays. In 2015, Xuan Dai Bay converted to pattern 4 (low intensity-medium concentration pattern). This result indicated that Xuan Dai Bay experienced moderate expansion and that the expansion tended to be concentrated in certain areas.

Vertical Trends
The overall proportions of construction land in each zone are shown in Table 8. In 1988, the proportion of construction land showed a trend of decrease-increase-decrease with two peaks, one in 0-2 km zone and another in 6-8 km zone. After 30 years of rapid construction sprawl, the proportions in all the zones increased greatly. The rapid pace of the construction sprawl occurred in the 0-2 km zone, and the proportion rose to 17.1%. The growth rate was much greater than that in other zones. In the 2-4 km, 4-6 km and 6-8 km zones, the growth rate was approximately 6.5%. The proportions in the 8-10 km zone were relatively low, and the growth rate was the lowest. The results suggested that the proportions of construction land decreased as the distance from the sea increased. The gap has widened over the following 30 years. The proportions of construction land in each zone and each bay are shown in Figure 7. The results can be discussed in two ways according to the proportion of construction land. In the bays with high construction land density, the distributions of construction land in the five zones were similar in the two periods. The results suggested that the construction land expanded inside the zone. In Hai Loc Bay and Dam Ha Trung Bay, the distribution of construction land in each zone was relatively even. The built-up density in the inland zone was slightly greater than that in the seaward zone. In Danang Bay, the construction land was concentrated in the 0-2 km and 2-4 km zones in both 1988 and 2015. From the 4-6 km zone to the 6-8 km zone to the 8-10 km zone, the proportion of construction land decreased successively. Otherwise, in the bays with low construction land density, the distributions of construction land in the five zones in 2015 were different from the distributions in 1988. The rate of construction sprawl in the 0-2 km zone was obviously greater than that in the other zones. In the other zones, the construction land expanded randomly.        Figure 9 shows the distribution of the PDs in the five buffers of the bays in the two periods. Over the study period, the PDs increased in almost all zones and all bays. Among them, the PDs in Hai Loc Bay were much greater than those in other bays, and the increase of the PDs was also the greatest. From the distribution throughout all zones, the PDs in the 0-2 km zone ranked lowest or second-lowest in Hai Loc Bay, Dam Ha Trung Bay, Danang Bay and Ganh Rai Bay. Accordingly, the four bays had the greatest proportion of construction land. Additionally, the PDs generally showed an increasing trend from the 0-2 km zone to the 8-10 km zone, especially in the first three zones (0-6 km zone). In other bays, the PDs changed little or decreased in the five zones. The results suggested that the fragmentation degree was related to the development degree. In the highly developed bays, the construction land showed an aggregated distribution in the 0-2 km zone. As the distance from the sea increased, the construction land tended to be more fragmented.

Spatial Trends with Distance from Coast
ISPRS Int. J. Geo-Inf. 2020, 9, x FOR PEER REVIEW 14 of 21 Figure 9 shows the distribution of the PDs in the five buffers of the bays in the two periods. Over the study period, the PDs increased in almost all zones and all bays. Among them, the PDs in Hai Loc Bay were much greater than those in other bays, and the increase of the PDs was also the greatest. From the distribution throughout all zones, the PDs in the 0-2 km zone ranked lowest or secondlowest in Hai Loc Bay, Dam Ha Trung Bay, Danang Bay and Ganh Rai Bay. Accordingly, the four bays had the greatest proportion of construction land. Additionally, the PDs generally showed an increasing trend from the 0-2 km zone to the 8-10 km zone, especially in the first three zones (0-6 km zone). In other bays, the PDs changed little or decreased in the five zones. The results suggested that the fragmentation degree was related to the development degree. In the highly developed bays, the construction land showed an aggregated distribution in the 0-2 km zone. As the distance from the sea increased, the construction land tended to be more fragmented.  Figure 10 shows the distribution of the LSIs in the five buffers of the bays in the two periods. Over the period, the LSIs increased in all zones and all bays except Danang Bay. In Danang Bay, the LSIs in the 0-2 km zone and 2-4 km zone decreased in 1988. As the development degree increased, the construction land tended to gather, and the shape tended to be simple. In most bays except Hai Loc Bay and Danang Bay, the LSIs generally showed a decreasing trend from the 0-2 km zone to the 8-10 km zone. The results suggested that as the distance from the sea increased, the proportion of construction land tended to decrease, and the shape tended to be simpler.  Figure 10 shows the distribution of the LSIs in the five buffers of the bays in the two periods. Over the period, the LSIs increased in all zones and all bays except Danang Bay. In Danang Bay, the LSIs in the 0-2 km zone and 2-4 km zone decreased in 1988. As the development degree increased, the construction land tended to gather, and the shape tended to be simple. In most bays except Hai Loc Bay and Danang Bay, the LSIs generally showed a decreasing trend from the 0-2 km zone to the 8-10 km zone. The results suggested that as the distance from the sea increased, the proportion of construction land tended to decrease, and the shape tended to be simpler.  Figure 11 shows the distribution of the AIs in the five buffers of the bays in the two periods. In most bays except Dam Ha Trung Bay, the AIs generally showed a decreasing trend from the 0-2 km zone to the 8-10 km zone. The results suggested that as the distance from the sea increased, the construction land tended to be more discrete.  Figure 11 shows the distribution of the AIs in the five buffers of the bays in the two periods. In most bays except Dam Ha Trung Bay, the AIs generally showed a decreasing trend from the 0-2 km zone to the 8-10 km zone. The results suggested that as the distance from the sea increased, the construction land tended to be more discrete.  Figure 11 shows the distribution of the AIs in the five buffers of the bays in the two periods. In most bays except Dam Ha Trung Bay, the AIs generally showed a decreasing trend from the 0-2 km zone to the 8-10 km zone. The results suggested that as the distance from the sea increased, the construction land tended to be more discrete.

Discussion
Landscapes influence the proportions and spatial forms of construction land in the bays. Table  9 shows the bay types as well as the proportions and spatial patterns of construction land in each bay. The bays with high construction land proportions belonged to delta bay, bedrock erosion-stacking bay and lagoon bay. The influence of these three types of bays on the distribution of construction land is discussed next.

Discussion
Landscapes influence the proportions and spatial forms of construction land in the bays. Table 9 shows the bay types as well as the proportions and spatial patterns of construction land in each bay. The bays with high construction land proportions belonged to delta bay, bedrock erosion-stacking bay and lagoon bay. The influence of these three types of bays on the distribution of construction land is discussed next. The abscissa values are explained as follows: HH: high intensity-high concentration pattern; HL: high intensity-low concentration; MM: medium intensity-medium concentration pattern; LM: low intensity-medium concentration pattern; and LL: low intensity-low concentration pattern.
(1) Delta Bay Deltas are the result of sediment accumulation in estuaries and adjacent areas when rivers enter the sea [46]. Deltas and their estuaries hold both ecological and economic value and are major centers of population and agriculture [47]. A delta bay is a bay formed with the development of a constructive delta [39]. It is easy to convert agricultural infrastructure into urban infrastructure in the urban fringe of deltaic cities [48]. In return, the intensity of development in the delta bay is often high. Hai Loc Bay and Ganh Rai Bay are delta bays located in the two largest deltas in Vietnam-the northern Red River Delta and the southern Mekong River Delta, respectively ( Figure 12). These two deltas are currently agricultural hotspots of Vietnam, contributing 71% of the rice production, 86% of the farmed aquaculture and 65% of the fruit production of the country [49]. Hanoi city and Ho Chi Minh City are the largest cities in the two river deltas, as well as the largest cities in Vietnam. The built-up areas in Hanoi and in Ho Chi Minh have expanded in the past 20 years [50,51]. Correspondingly, the densities of the built-up areas were high in the two bays, and the growth rates of construction land were also very high. The proportions of construction land increased from 15.4% to 24.5% and from 3.5% to 15.4% in Hai Loc Bay and Ganh Rai Bay, respectively. (2) Lagoon Bay A lagoon bay is semi-enclosed water composed of a lagoon. The influx of nutrients from many inland rivers entering the lagoon supports high levels of primary productivity in the lagoon. This productivity, in turn, supports fishery production. Aquaculture and fishery capture are considered the most important sources of income in the region [55]. In addition, the sandbanks outside the lagoon can be used as breakwaters, which are suitable for small and medium-sized ports. This has a radiating effect on the mainland economy. Dam Ha Trung Bay and Cam Ranh Bay are located at Tam Giang-Cau Hai Lagoon and Thuy Trieu Lagoon, respectively ( Figure 13). In terms of spatial form, Dam Ha Trung Bay and Cam Ranh Bay presented the MM pattern. The proportion of construction land in Dam Ha Trung Bay was greater than that in Cam Ranh Bay. The Tam Giang-Cau Hai Lagoon, where Dam Ha Trung Bay is located, lies along the central coastal line of Vietnam. The lagoon runs nearly 70 km along the coast and covers 219.18 km 2 of the water surface, forming the largest lagoon system in Southeast Asia [55,56]. Thuy Trieu Lagoon, where Cam Ranh Bay is located, is a narrow enclosed water body in southeastern Khanh Hoa Province that is approximately 16 km in length, 250 m in width and 0.5 to 6 m in depth [57]. The lagoons experienced the rapid expansion of areas devoted to aquaculture activities (notably, culture ponds) [56,57]. Accordingly, many settlements are scattered along the bay coast. In addition, the developed agriculture in inland areas also gave rise to regional agricultural centers. The distribution of construction land appeared two peaks in 0-2 km zone and 6-8 km zone. From the perspective of the spatial morphology, the construction land in the two bays presented blocked or banded distributions.  The terrain of the bay has a great influence on the spatial forms of construction. The delta area often has a dense water network, flat terrain, fertile soil, developed agriculture. The river network in the Red River Delta is quite dense, with a density of about 0.7 km/km 2 [52]. The Mekong Delta river network is relatively dense and complex and includes natural river systems as well as manmade canals [53]. In this case, it was easy to form dense rural settlements, especially in the narrow border road of the rice field [54]. In terms of spatial form, Hai Loc Bay presented the HL pattern in the two periods. The construction land in Hai Loc Bay had the characteristics of high density and discrete distribution. Meanwhile, the proportion of construction land was even in each zone. At the edge of the bay, there are large uninhabitable depressions, leading to an increase in the proportion of construction land in the sea-land direction. The spatial form of construction land in Ganh Rai Bay was different from that in Hai Loc Bay. Ganh Rai Bay is located in the transitional zone between the southeast platform and the Mekong Delta region. Correspondingly, the construction land in the deltaic area was also relatively discrete, which was similar to that in Hai Loc Bay. In the eastern area of Ganh Rai Bay, the construction land was concentrated in the seaward zone, showing significantly banded distributions.
(2) Lagoon Bay southeastern Khanh Hoa Province that is approximately 16 km in length, 250 m in width and 0.5 to 6 m in depth [57]. The lagoons experienced the rapid expansion of areas devoted to aquaculture activities (notably, culture ponds) [56,57]. Accordingly, many settlements are scattered along the bay coast. In addition, the developed agriculture in inland areas also gave rise to regional agricultural centers. The distribution of construction land appeared two peaks in 0-2 km zone and 6-8 km zone. From the perspective of the spatial morphology, the construction land in the two bays presented blocked or banded distributions. effect on the mainland economy. Dam Ha Trung Bay and Cam Ranh Bay are located at Tam Giang-Cau Hai Lagoon and Thuy Trieu Lagoon, respectively ( Figure 13). In terms of spatial form, Dam Ha Trung Bay and Cam Ranh Bay presented the MM pattern. The proportion of construction land in Dam Ha Trung Bay was greater than that in Cam Ranh Bay. The Tam Giang-Cau Hai Lagoon, where Dam Ha Trung Bay is located, lies along the central coastal line of Vietnam. The lagoon runs nearly 70 km along the coast and covers 219.18 km 2 of the water surface, forming the largest lagoon system in Southeast Asia [55,56]. Thuy Trieu Lagoon, where Cam Ranh Bay is located, is a narrow enclosed water body in southeastern Khanh Hoa Province that is approximately 16 km in length, 250 m in width and 0.5 to 6 m in depth [57]. The lagoons experienced the rapid expansion of areas devoted to aquaculture activities (notably, culture ponds) [56,57]. Accordingly, many settlements are scattered along the bay coast. In addition, the developed agriculture in inland areas also gave rise to regional agricultural centers. The distribution of construction land appeared two peaks in 0-2 km zone and 6-8 km zone. From the perspective of the spatial morphology, the construction land in the two bays presented blocked or banded distributions.  (3) Bedrock erosion-stacking bay A bedrock erosion-stacking bay is mainly a depression or valley formed by long-term erosion of the rock fracture zone and seawater intrusion. The top of the bay often forms a sandy or tidal beach, and the two sides have bedrock capes [40]. Semi-enclosed bedrock bays are natural deep-water harbors, and the tidal channels and scour troughs in these bays are natural channels. The urban layout gives priority to the bay, and the economic structure is mainly based on the bay. Port economy drives the development of the marine industry and the inland economy. Such bays have always been the focus of coastal zone development [58]. Danang Bay is a bedrock erosion bay on the central coast of Vietnam ( Figure 14). Da Nang city, where Danang Bay is located, is the most dynamic city in the middle of central Vietnam [59]. The proportion of construction land was greatest in Danang Bay in 2015, up to 26.5%. In terms of spatial form, Danang Bay presented the HH pattern in the two periods. The bay has lower terrains in the middle and hills on the two sides. The construction land was concentrated in central coastal plains with a high degree of aggregation. In the past 30 years, the construction land has expanded inland from the original core. The proportions of construction land decreased significantly in the sea-land direction. The major seaport in Vietnam-Da Nang port is located at the river mouth [60]. The greatest density of construction land was also located in this region. (3) Bedrock erosion-stacking bay A bedrock erosion-stacking bay is mainly a depression or valley formed by long-term erosion of the rock fracture zone and seawater intrusion. The top of the bay often forms a sandy or tidal beach, and the two sides have bedrock capes [40]. Semi-enclosed bedrock bays are natural deep-water harbors, and the tidal channels and scour troughs in these bays are natural channels. The urban layout gives priority to the bay, and the economic structure is mainly based on the bay. Port economy drives the development of the marine industry and the inland economy. Such bays have always been the focus of coastal zone development [58]. Danang Bay is a bedrock erosion bay on the central coast of Vietnam ( Figure 14). Da Nang city, where Danang Bay is located, is the most dynamic city in the middle of central Vietnam [59]. The proportion of construction land was greatest in Danang Bay in 2015, up to 26.5%. In terms of spatial form, Danang Bay presented the HH pattern in the two periods. The bay has lower terrains in the middle and hills on the two sides. The construction land was concentrated in central coastal plains with a high degree of aggregation. In the past 30 years, the construction land has expanded inland from the original core. The proportions of construction land decreased significantly in the sea-land direction. The major seaport in Vietnam-Da Nang port is located at the river mouth [60]. The greatest density of construction land was also located in this region.

Conclusions
This study used Landscape indexes, clustering method and strip division to analyze the spatial and vertical distribution of construction land in the bay area. Additionally, the influence of bay type

Conclusions
This study used Landscape indexes, clustering method and strip division to analyze the spatial and vertical distribution of construction land in the bay area. Additionally, the influence of bay type on the distribution of construction land was discussed. Our findings were as follows: (1) Across the region, the total proportion of construction land increased from 4.9% to 12.5% from 1988 to 2015. In terms of the spatial distribution, there were significant differences between the northern bays and southern bays. The top three bays with the greatest proportions (Hai Loc Bay, Dam Ha Trung Bay and Danang Bay) were distributed along the northern coast of Vietnam. (2) In the process of the expansion of construction land in the bays, the bay-level LPIs, PDs and LSIs increased during the 30 year study period. This result indicated that the distribution of construction land became more fragmented, and the shape tended to be more complex in the 30 year period. Especially in Hai Loc Bay, construction land was the most densely and fragmented distributed. (3) The method of using the landscape indexes for clustering can identify the spatial pattern of construction land effectively. The eight bays were divided into five patterns (HH, HL, MM, LM and LL) into the two periods. Danang Bay exhibited a typical HH pattern, and Hai Loc Bay exhibited a typical HL pattern. (4) The distance from the sea had a significant influence on the distribution of construction land in bays. There was a trend that the proportions of construction land decreased as the distance from the sea increased in general. The proportions of construction land in the 0-2 km zone were greatest, as well as the growth rate. In contrast, the proportions in the 8-10 km zone were relatively low. Moreover, the construction land tended to be more fragmented and discrete, and the shape tended to be simpler. (5) The bay type had a significant influence on the distribution of construction land in the bays.
The proportions of construction land were high in delta bay, bedrock erosion-stacking bay and lagoon bay. In terms of the spatial distribution, the construction land was densely and fragmented distributed in delta bay, while the construction land was centrally distributed in bedrock erosion-stacking bay.
The study provided a method to analyze the spatial pattern of construction land distribution. In future studies, (i) more bay data should be collected to improve the classification system of the spatial pattern of construction land, and (ii) more landscape indicators should be identified to reveal the pattern of construction land expansion in bays.