Coastline Fractal Dimension of Mainland, Island, and Estuaries Using Multi-temporal Landsat Remote Sensing Data from 1978 to 2018: A Case Study of the Pearl River Estuary Area

: The Pearl River Estuary Area was selected for this study. For the past 40 years, it has been one of the most complex coasts in China, yet few studies have analyzed the complexity and variations of the area’s di ﬀ erent coastlines. In this investigation, the coastlines of the Pearl River Estuary Area were extracted from multi-temporal Landsat remote sensing data from 1978, 1988, 1997, 2008, and 2018. The coastline of this area was classiﬁed into mainland, island, and estuarine. To obtain more detailed results of the mainland and island, we regarded this area as the main body, rezoned into di ﬀ erent parts. The box-counting dimension was applied to compute the bidimensional (2D) fractal dimension. Coastline length and the fractal dimension of di ﬀ erent types of coastline and di ﬀ erent parts of the main body were calculated and compared. The fractal dimension of the Pearl River Estuary Area was found to have increased signiﬁcantly, from 1.228 to 1.263, and coastline length also increased during the study period. The island and mainland showed the most complex coastlines, while estuaries showed the least complexity during the past forty years. A positive correlation was found between length and 2D-fractal dimension in some parts of the study area. Land reclamation had the strongest inﬂuence on fractal dimension variations.


Introduction
The coastal area is one of the most significant elements of Earth's ecosystems, providing habitats for humans and other creatures [1]. Serving as the boundary between land and sea, the coastline is influenced by natural and anthropogenic factors [2,3]. It is essential to understand the length variations and irregularities of coastlines, as reflected by fractal dimension. Since the implementation of the Reform and Opening-up policy in China in 1978, the coastal area of China has shown dramatic variations due to the significant requirement of land for economic progress [4]. For example, substantial changes appeared in the Pearl River Delta in the southeast of China, which is the pioneer of reformation and development [5][6][7]. With its unique location, frequent interactions between sea and land have occurred in the Pearl River Delta area which have led to significant coastline variations from the natural. In addition to natural influences, anthropogenic activities have taken place in this area over the past 40 years, with the aim of boosting the construction of facilities [8]. The reclamation of land towards the sea has occurred in the coastal area to create more land resources, not only in this particular area, but also in all of China [9]. A detailed understanding of coastlines and their variations in China can provide useful information for coastal management, coastal zoning planning, and coastal preservation.

1.
To classify the coastline of the Pearl River Estuary as estuarine, mainland, and island coastlines, and then calculate the coastline length and FD of each type of coastline.

2.
To analyze the reasons for the length increases and differences between different types of coastlines and discuss comparisons between estuaries and Parts I-IV of the main body. 3.
To analyze the reasons for variations in FD throughout the 40-year time period and how different circumstances influenced the FD. 4.
To discuss how natural and anthropogenic factors impacted on FD and coastline length, and compare the FD of the study site with other coastal areas of China.

Study Area
The PREA is located in the southeast of China ( Figure 1). It is the second largest river system in China, the largest estuary in Southern China, and one of the most complex fluvial networks in the world [5,31]. The annual mean temperature ranges from 14 to 22 • C [32] and the annual mean precipitation is 1470 mm [33]. The dry season is from September to March in the next year, and the wet season is from April to August [34]. During the wet season, the precipitation accounts for 72-88% of the total annual rainfall [32]. The average water depth of the study site is 30 m [9].  In this study, Landsat remote sensing data were used to extract the coastlines of the study site 150 from 1978 to 2018 (Table 1). All datasets were downloaded from the Geospatial Data Cloud Weak tidal conditions and low sediment concentrations are found in this area [34]. The tidal condition of the PREA is mixed semidiurnal [35]. The average tidal range is from 86 cm to 163 cm [9]. Previous observational data suggests that the significant wave height of the PREA is less than 0.5 m in summer [35]. Some scientists considered the tidal condition of the PREA to be weak [34], and low sediment concentrations were found in this area [34]. Recently, the mean annual sediment decreased to 3.04 × 10 7 tons, with nearly 95% of the sediment being transported during the wet season [35]. Coarse sediment was the most dominant near estuaries and channels [35].
The unique geographic characteristics of the PREA include eight estuaries, namely Humen, Jiaomen, Hongqimen, Hengmen, Modaomen, Jitimen, Hutiaomen, and Yamen ( Figure 1). The gross domestic product in this area has soared since 1978 due to the Reform and Opening-up policy [8,36]. More than 8636 reservoirs had been established by the late 1990s [33].
The PREA encompasses an extremely large area, so for a better understanding, the mainland and island were rezoned into smaller parts. The mainland coastline and island coastline of the PREA were regarded as the main body of the entire PREA in this study. The area was divided into four parts: Part I, Part II, Part III, and Part IV (Figure 1), and each part consisted of mainland and island coastline. Due to their unique geological positions, different parts show different coastline variations.
The main body of the PREA was divided into Parts I-IV, based on their unique locations. The mutual characteristic among these parts is human activities that have facilitated coastline change during the study period. Part I is located in the east of the PREA and includes the rapid development of Shenzhen during the past forty years. In this part, land reclamation focused on the mainland area and islands, showing a smooth increase throughout. This part does not contain estuaries. Part II shows significant past changes of island reclamation. This part exhibited the most significant island variations compared with the other three parts. Similar to Part II, the islands of Part III also show significant changes. This part contains the Modaomen and Macau Special Administrative Region. The coastline variations of Part IV differ from the other parts. Land reclamation has filled the area between the mainland and islands since 2008; therefore, the FD of this area shows different results to the other parts.

Remote Sensing Dataset
In this study, Landsat remote sensing data were used to extract the coastlines of the study site from 1978 to 2018 (Table 1). All datasets were downloaded from the Geospatial Data Cloud (https://www.gsclound.cn/) and the USGS (https://www.usgs.gov/). Data with nearly no cloud cover in the study site were required for this study, as well as dry season time period data. A digital elevation model of Guangdong Province with a spatial resolution of 30 m was downloaded from: http://www.tuxingis.com/. All tidal data were collected from the South China Branch of the State Oceanic Administration (Table 1). slope for the study area, calculated using ArcGIS, 10.3 is 5.39%. The error of the Landsat images would be less than 30 m (spatial resolution) if the tidal difference value is less than 161.7 cm [37], hence, the tidal influence data we used in this study were less than one pixel.
The flow chart of the entire process is shown in Figure 2. The coastline was extracted using automatic and visual interpretation in ArcGIS 10.3 and MATLAB 2014a. The box-counting method was applied via ArcGIS to calculate the 2D-FD of the different coastlines.

Coastline Extraction and Accuracy Assessment
For the coastline extraction, an automatic and visual interpretation method was used with the Landsat dataset (Table 1). Image analysis software (ENVI) was applied for image pre-processing, including radiometric calibration, FLAASH atmospheric correction, and stripe repair for Landsat 7 ETM+. Canny Edge Detection in MATLAB was used to extract the coastline and examine by visual interpretation in ArcGIS 10.3. The coastline length was calculated in ArcGIS 10.3 using the Calculate Geometry tool.
A precision test was applied to evaluate the accuracy of the extracted coastline, and an image was obtained from Google Earth™ Pro. To maintain a consistent time, the Landsat 8 OLI_TIRS data were collected on February 12, 2018, and the Google Earth image was collected on January 18, 2018. Eighty-six ground control points (GCPs) were selected from the Google Earth image ( Figure 3). The distances between each GCP and the corresponding point on Landsat were calculated (Figure 3c) [9]. The average error of all GCPs and the corresponding points on Landsat, the standard deviation, and the root-mean-square error of these points were calculated using the distances obtained between each GCP and corresponding point. Using these data and the theoretical maximum allowable error proposed by Hou et al (2016), the accuracy of the coastline extraction was obtained [9]. The results of the precision test are shown in Table 2.
each GCP and corresponding point. Using these data and the theoretical maximum allowable error 183 proposed by Hou et al (2016), the accuracy of the coastline extraction was obtained [9]. The results of 184 the precision test are shown in Table 2.      The coastline was classified in ArcGIS 10.3 by different locations including the estuary, mainland, and island ( Figure 1) by visual interpretation. The estuary coastline represents approximately eight outlets, and usually extends from the beginning of the junction of the tributaries to the end of the outlet to the sea. The mainland coastline incorporates the boundary between the mainland and sea, while the island coastline is a terrestrial and marine boundary of islands ( Figure 1). It is apparent that the island coastline is the boundary between the island and the sea (Figure 1). In this study, the mainland and island coastline of the PREA were regarded as the main body of the entire PREA. This area was divided into four parts, namely Part I, Part II, Part III, and Part IV ( Figure 1). Each part comprised mainland and island coastline. The FD of the different coastlines were then calculated in ArcGIS 10.3.

Fractal Dimension Calculation
The box-counting method was selected for the following reasons. From previous studies, several researchers commented that the box-counting method exhibited inconsistency whereby the box size and different scales also led to different FD [14,38]. The level of confidence of this method reached 99.9%, more than 99% from previous reports [39,40]. Compared with other methods, the box-counting method could provide the FD of discontinuous coastlines [41], which was the main reason it was selected for this study. In addition, research has found that different FD values would be produced when using different methods [29,41]. For example, the divider method obtained a larger value compared to the box-counting method in the same study area; therefore, the results from different methods cannot compare with each other directly. The result obtained from the box-counting method can compare with other studies calculated using the same method, therefore providing a more comprehensive analysis. Furthermore, Li et al. proposed that the box-counting method is easy to conduct using related automatic software [22].
The box-counting method differs from the divider method using different size of squares to cover the same coastline. When the size of square changed, the number of used squares would change. Equation (1) shows the mathematical process for calculating the fractal dimension with a series of square sizes. N represents the number of squares and r represents the size of the squares. The fractal dimension of coastline obtained is a value between 1 and 2. A larger value represents a more complex coastline.
The coastline length would not change when squares of different sizes are applied. Figure 4 shows a number of different square sizes to cover the same coastline, with the left using 1200 m squares and the right using 2700 m squares. The coastline length of (a) and (b) are the same, being 108.28 km; however, the change of grid size led to a different number of grids required.
Further analysis included a FD comparison among different coastline types, different locations (Part I, Part II, Part III, and Part IV), and different outlets were compared with previous studies. The relationship between coastline length and FD were also analyzed. shows a number of different square sizes to cover the same coastline, with the left using 1200 m 221 squares and the right using 2700 m squares. The coastline length of (a) and (b) are the same, being   When the counting of the number of squares was completed, a linear regression on square quantity log(N) and square size log(r) was conducted in Excel, as well as the R 2 (coefficient of determination). A further analysis of the combined FD and coastline length was also conducted. A linear regression between the coastline length and FD was calculated, and the R 2 (coefficient of correlation) showed the relationship between these two parameters.

Results
In this section, the length and 2D-fractal dimension of different types of coastline and different parts of the main body are described.

Result of Accuracy Assessment
The result of the precision test is shown in Table 2. We defined group data as each GCPs from Google Earth and the corresponding points on the Landsat dataset. Table 2 shows the average error as 16.21 m, which is smaller than the theoretical maximum allowable error of 28.28 m put forward by Hou et al. in 2016 [42]. Therefore, the coastline extracted in this study meets the accuracy standard.

Length of Coastlines
Tables 3 and 4 list the length of different types of coastlines, including the entire study site, mainland, islands, eight estuaries, and four parts. For the four parts, the total length of coastline length represented the sum of the mainland coastline length and the island coastline length of the respective part (Table 4). Tables 3 and 4 demonstrate that the coastline length of the PREA increased from 1133.43 km in 1978 to 1349.73 km in 2018. The mainland, island, and coastlines of the estuaries showed an increasing trend ( Table 2). The mainland coastlines had an increased length, from 371.71 km in 1978 to 510.36 km in 2018. This represents the most rapid growth rate. The length of the island placed in the second sector changed its rate and improved by almost 100 km during the study period. Unlike the significant growth speed of these two types, the estuary coastlines demonstrated a stable trend with a smooth increase of 1-10 km among the eight estuaries (Table 3). Among the eight estuaries, Modaomen showed the most significant variations, with an increase of nearly 10 km during the study period. The Jiaomen, Hongqimen, Hutiaomen, and Yamen estuaries experienced a slight increase of approximately 1 km. The Humen, Hengmen, and Jitimen estuaries increased by nearly 5 km. Compared with the mainland and island coastlines, the estuaries of the PREA remained stable.
For Parts I-IV, different parts showed different trends. The mainland and island of Part I increased over the past four decades, with the mainland increasing by nearly 30 km and the island increasing by nearly 10 km. Part II showed significantly increased island coastlines, advancing nearly 100 km and a slight increase of mainland coastline during the time period. Similar to Part II, the island coastline increased sharply during the time period, increasing approximately 60 km, while the mainland coastline grew gradually from 85.05 km initially to 94.68 km at the end of the forty years. Part IV demonstrated coastline variations that were different to the other three parts. Its total coastline length decreased slightly from 239 km to 234.54 km during the forty-year time period. The island coastline decreased significantly, from 142.81 km in 1978 to 60.58 km in 2018. The mainland coastline increased sharply during the time period. In Part IV, land reclamation filled the gap between the mainland and island, whereby the island coastline became mainland coastline and led to coastline variations of this area.

Fractal Dimension
The Tables 5 and 6 values show that all confidence levels (R 2 ) were greater than 0.99, indicating that coastlines in the PREA can be described using fractal theory [41,43].  Table 5 describes the FD of the PREA, mainland, island, and estuaries. It is clear that the fractal dimension of coastlines, including the entire study area, mainland, island, and Jiaomen shows a growth, with some fluctuations during the time period. Table 6 includes the FD of four parts of the main body shown in Figure 1. Each part included the FD of the entire part, mainland coastline, and island coastline. As FD represented the complexity of the coastline, we can infer that the complexity of the PREA increased during the study period. Furthermore, we found that the island coastline showed the greatest complexity and the mainland showed the second most complexity when comparing different coastline types. All estuary coastlines demonstrated the least complexity.
It can be seen from Tables   When combining Table 5 and Figure 5, the FD of the estuary coastline remained between 1.00 and 1.10. All estuaries exhibited a decline during the study period, with the exception of Jiaomen, which increased from 1.023 in 1978 to 1.027 in 2018. In 1978, this value differed significantly for Jitimen, with the maximum value of 1.097 and Yamen with the minimum value of 1.011. Nevertheless, the difference reduced in 2018. A sharp decline was found during the first decade, from 1978 to 1988. After 1988, the FD of the estuaries exhibited a gradual decrease in the following decades. Based on the complexity reflected by the FD, the complexity of the coastline of the estuaries in 2018 could be found as: Jitimen > Hutiaomen > Hengmen > Hongqimen > Jiaomen > Modaomen > Humen > Yamen. Figures 6 and 7 describe the mainland and island FD of Parts I-IV and the PRE, respectively. In       Figure 8a-d described the FD of Parts I-IV. The island coastline generally obtained a larger FD value than the mainland coastline among these four parts. The FD of the PREA also presented this phenomenon ( Figure 5). Furthermore, the island FD in Part IV decreased during the forty-year time period. Combined with Table 4, we found that the island coastline length of Part IV decreased throughout the time period, but Parts I-III increased. From this, we suggest that there are some relationships between coastline length and FD.
period. Combined with Table 4, we found that the island coastline length of Part IV decreased 341 throughout the time period, but Parts I-III increased. From this, we suggest that there are some 342 relationships between coastline length and FD.      Figure 9 demonstrates the changing trend of coastline length and FD for all types of coastline. When examining the relationship between coastline length and FD, it is found that the variation trend between coastline length and FD of the PREA, mainland, island, and Jiaomen is similar. For other estuaries, these two parameters are inconsistent with Figure 8. The coastline complexity of Humen, Hongqimen, Hengmen, Modaomen, Jitimen, Hutiaomen, and Yamen, therefore, did not correspond with the coastline length change. Figure 10 describes the changing trend of coastline length and FD of the different parts of the main body, from Part I-IV. Nearly all the figures (Figure 10a-l) show that the FD of these parts was almost consistent with coastline length. Therefore, for the main body, the increased coastline length represents a more complex coastal boundary and the decreased coastline length represents a less complex coastal boundary.   coastline variations [44][45][46][47]. In addition to these natural factors, anthropogenic factors also influence 376 the length of coastlines.

377
In this study, we found that the coastline of the PREA increased significantly by 197 km, from

Driving Factors of Coastline Variations
Coastline changes are mainly influenced by natural and anthropogenic factors. Coastal erosion and accretion are natural factors leading to coastlines variations. Other natural factors such as sediment transport, tidal condition, wind, temperature, flow direction, and wave energy also lead to coastline variations [44][45][46][47]. In addition to these natural factors, anthropogenic factors also influence the length of coastlines.
In this study, we found that the coastline of the PREA increased significantly by 197 km, from 1133.43 km to 1349.73 km during the 40-year time period studied. The mainland coastline increased by 139 km, the island coastline increased by 65 km, and the coastline of the estuaries increased by 1-10 km.
In the PREA, the mainland and island coastlines were mainly influenced by anthropogenic factors, particularly land reclamation. Since 1978, land reclamation increased in the PREA [48], producing more land resources and longer artificial coastlines. To some extent, increased artificial coastlines will facilitate an increase in total coastline length. In the entire PREA, land filling was quite common over the forty years, and led to increased land area between the mainland and the island. Parts I-IV showed different coastline change characteristics. The coastlines of Part I mainland, Part II island, Part III island, and Part IV mainland were mainly impacted by land reclamation and increased coastline length to create more land resources. The significantly increased Part IV mainland coastline length is seen in Table 4, with a sharp decrease of Part IV island coastline length. In this part, more land was generated seawards and the gap between the mainland coastline and the island was filled to satisfy the demand for land. Although different parts showed different coastline change trends, land reclamation was the most significant reason leading to this phenomenon.
Unlike the mainland and the island, the coastline of the estuaries was impacted primarily by natural factors with some anthropogenic factors, which have increased gradually during the past four decades. From the coastline extraction, we found that less land reclamation was found around estuaries; therefore, coastline erosion and sediment transport would be the most important factors. Sediment transport via the eight estuaries could reach 3.04 × 10 7 ton per year [35]. A large amount of sediment transport may cause erosion of the PREA. Heise et al. (2010) assumed that a larger amount of sediment was deposited in the PREA than was eroded [49]. Furthermore, dams also impact upon sediment flux [50]. Dai et al. (2008) found that the sediment flux of the Pearl River had decreased since the 1980s [50] and suggested that more sediment was deposited in the reservoirs. As more dams were established in the PREA, less sediment flux would be transported into the sea. Wu et al. also indicated that human activities influence the sediment load carried into the sea [33]. For the coastline of estuaries, anthropogenic factors would therefore influence coastline length indirectly by influencing sediment transport.
The driving factors of coastline length of the PREA depend on different coastline types. Mainland and island coastlines are mainly influenced by land reclamation, but for estuary coastlines, the reason for their variations in length is mainly influenced by natural factors including erosion and sediment transport, and indirectly by human activities.

Driving Factors of FD Variations
The fractal dimension can reflect the complexity of a coastline; therefore, a changed coastline length would lead to a variation in coastline complexity. Figures 8 and 9 demonstrate that coastline length and FD generally had the same trend during the study period. When the coastline length increased, the FD also increased. For estuaries, the trend between the coastline length and the FD was different, with the exception of Jiaomen.
As we analyzed the most significant influences of coastline length, it became apparent that the mainland and the island experienced land reclamation. As more artificial coastline was created around these two types of coastline, the shape of the coastal boundary changed simultaneously. This also led to the variation in coastline complexity, which was reflected by the FD. Increased artificial coastlines facilitate increased total coastline length and also lead to increased FD. A study by Zhang et al. (2015) reached a similar conclusion, whereby the coastline length became longer and the coastline FD of all China increased, with fluctuations as the proportion of artificial coastlines increased over the past forty years [39]. At this time, the original coastline became more complex, although it became smoother after building a wharf, an airport, and other facilities. In addition to the PREA, studies of the Yellow River Delta and all of China also found that as the complexity of coastlines increased, the proportion of smoother artificial coastline simultaneously increased. It could be considered that the total length of coastline grew as the length of artificial coastlines increased, and led to increased coastline complexity. Figure 9 shows detailed information of the mainland and island FD variations and also indicates that the complexity of coastline and coastline length show almost the same changing trend. Therefore, as the length of mainland and island coastline increased, the increased complexity reflects the increasing trend of the FD.
For the four parts of the main body, the eastern coast had a less complex coastline than the western coast. Parts II and III had a larger FD than the other two areas ( Table 3). The western part (Part II and IV) included eight estuaries and described a larger FD than the eastern part. The location of eight estuaries might be one of the natural influences leading to the differing FD among these four parts because sediment deposition would affect the FD. It was observed that Parts II and III had a larger FD than the other two parts, indicating that these two areas displayed a more complex coastline. The coastline length of Parts II and III changed more significantly than Parts I and II (Table 3). Land reclamation and the artificial coastline have increased at the Pearl River Estuary since 1978, and Parts II and III were influenced more by human factors. Anthropogenic activities could, therefore, be a significant reason for the different FD of the different locations.
In contrast to the mainland and island coastlines, the FD of estuaries found a different trend compared with coastline length (Figure 8). The analysis showed that the coastline length of the estuaries remained almost stable, with a slight increase of 1-10 km. Compared with the mainland and island coastlines, little land reclamation would be found among the eight estuaries, and thus the variation of the FD of the estuaries was impacted more by natural factors.
Sediment deposition may decrease the complexity of estuary coastlines because less sediment flux flows past estuaries and the coastline becomes rougher with a more complex condition [46]. In 2008, Su et al. commented that the geographic locations of objects and coastline variations would influence the FD [40]. In light of this, different locations of estuaries may lead to different coastline complexities. Furthermore, the shape of the coastline might be changed by natural factors, and the coastline complexity among the eight estuaries may be reduced. Coastal erosion and sediment transport are the most two significant factors that impact upon the complexity of the coastline.
In general, the FD variations in the PREA were caused by comprehensive natural driving factors, overlaid with anthropogenic factors. Different types of coastlines are influenced by different factors. Here, land reclamation, accompanied by an intense demand for land for economic development, infrastructure construction, and rapid urbanization in this area were the main factors driving the change of coastline complexity of the mainland and the island. Geographic location, coastal erosion, and sediment transport were the main factors in the estuaries. For Parts I-IV, the FD was mainly caused by land reclamation and different areas showed different coastline complexity. The complexity of the eastern coast (Part I) showed less complexity than the western parts (Parts II-IV).

Relationship Between Coastline Length and FD
Figures 9 and 10 demonstrate coastline length and FD variations of the PREA, mainland, island, estuaries, and Parts I-IV. In this section, further analysis of the relationship between coastline length and FD was conducted.
As shown, variations of coastline length and FD have nearly the same trend that increased during the time period. Figure 11 described a linear relationship between these two indicators of the coastlines in the entire PREA, mainland, island, and all estuaries. The closer the R 2 (coefficient of correlation) is to 1, the higher the correlation between these two indicators. This study found a positive correlation between the PREA, mainland, island, and two of the estuaries, including Jiaomen and Hutiaomen (>0.5). The remaining locations obtained a R 2 of <0.5 which is too weak to show a correlation between the coastline length and the FD. Among the areas with a positive relationship between coastline length and FD, the mainland coastline, island coastline, and Jiaomen also showed a consistent trend between these two parameters.
When combined with Figure 12, it is clear that the mainland and the island indicated a better relationship between coastline length and FD because the coastline length and the FD of the mainland and the island show almost the same trend as in Figure 10. The coastline length and complexity were impacted by land reclamation. A more regular artificial coastline was created during the study period, and the mainland and island coastlines were more significantly affected by anthropogenic factors. For these reasons, a positive relationship occurred between the length and FD of the mainland coastline, island coastline, and the majority of the Parts I-IV coastline.

516
The estuary coastline was the least complex coastline of the PREA. For the four parts of the main 517 body, the eastern coast had a less complex coastline than the western coast. Parts II and III had a 518 larger FD than the other two areas. The estuary coastline experiences less influence from artificial factors but more from natural factors. Although a rougher coastline was found around these areas, the coastline length of the estuaries remained mostly stable. Therefore, a negative relationship was found between the coastline length and the coastline complexity of the estuaries.
As the development of the PREA is still in progress, more land reclamation may occur around the mainland and island coastlines in the future; therefore, the complexity of these two types of coastline is likely to increase as the coastline increases. For the estuaries, due to the impacts of erosion, sediment transport, and dams, the complexity of this coastline is likely to decrease and the relationship between coastline length and FD may become weaker.

Conclusions
This study extracted coastlines of the mainland, island, and estuaries of the PREA using automatic and visual interpretation methods from Landsat data for 1978, 1988, 1997, 2008, and 2018. We reclassified the PREA as mainland and island coastline, and reclassified the main body of the PREA as Part I, Part II, Part III, and Part IV. We calculated their lengths and the FD of different coastlines using the box-counting method. The precision test for the coastline extracted in this study met the standard of precision evaluation proposed by Hou et al. (2016).
During the study period, the coastline length and the FD of the entire PREA increased significantly, with nearly the same tendency. A positive relationship between the coastline length and the FD was observed in the entire PREA, mainland, island, and two estuaries, including Jiaomen and Hutiaomen. Both the coastline length and complexity of these particular coastlines have increased over the past four decades, implying that these longer coastlines have led to a larger FD value and greater complexity. The island coastline was the most complex according to the FD result. The estuary coastline was the least complex coastline of the PREA. For the four parts of the main body, the eastern coast had a less complex coastline than the western coast. Parts II and III had a larger FD than the other two areas.
Land reclamation had the greatest influence on FD variations for the mainland and island coastlines. We found that more land filling occurred in the mainland and island, which changed the coastlines of these two areas, but land reclamation was rarely found around the estuaries. Land reclamation may increase the complexity of the coastline from previous studies. The FD value of Parts II and III indicated that more land reclamation was generated in these two areas because they had larger FD values.
For the coastline of the estuaries, coastal erosion, sediment transport, dam establishment, and other geographic objects were the main reasons for the FD variations. The complexity of the estuaries decreased as the sediment flux declined. The FD was mainly influenced by natural factors.
A different relationship between the coastline length and the FD was found among the mainland, island, and estuary coastlines. A positive relationship between length and FD was found for the mainland and the island. Nearly all estuaries showed a negative relationship between length and FD. As more land reclamation occurs around the PREA, government and citizens need to understand the importance of protecting ecosystems and finding a balance between ecological protection and economic development.