Analysis of Coastal Zone Data of Northern Yantai Collected by Remote Sensing from 1990 to 2018

Using remote sensing images of different time phases from 1990 to 2018, the surface coverage information of northern Yantai (coastline, 2 km from coastline to land) was extracted by means of average high tide line and visual interpretation. The end point change rate (EPR) and linear regression rate were used to study the coastline change rate, the fractal dimension of the coastline in the study area was analyzed, and the change of the type of coastal surface cover was analyzed by the transition matrix. The results show that: (1) Form 1990 to 2018, a significant trend of a continuous increase in the total length of coastline was observed with an increase of 181.08 km (43.18%). In the study area, the coastline of Laizhou had the greatest change rate with an EPR value of 33.67 m/a, whereas the coastline of Laishan had the smallest change rate with an EPR value of 0.30 m/a. (2) Over the past 30 years, with the rapid economic development of Yantai and the ensuant urbanization, the total surface area of the coastal zone in the study area has increased by 144.94 km2, mainly in the areas covered by structures and forests/grasses, by 112.96 km2 and 96.08 km2, respectively, while the areas of desert/bare land and water have decreased by 92.26 km2 and 12.32 km2, respectively. (3) The changes among different types in the study area were clear, mainly from desert/bare land, cultivated land, and building areas to forests/grasses cover and structures. The change areas were mainly concentrated in Laizhou, Longkou, Zhifu, and Penglai. Frequent human activities are an important factor affecting the continuous expansion of the coastal areas of Jiaodong Peninsula to the sea. Aquaculture, coastal construction, construction of artificial islands, and expansion of port terminals have seriously affected the sustainability of ecological resources in the coastal areas. At the same time, the changes in the ecological environment in the coastal zone will have a greater impact on the health of the coastal zone.


Introduction
A coastal zone is an interface between land and ocean. Coastal zones are the most densely populated type of area and are experiencing rapid economic development [1]. The economy of coastal zones plays a very important role in promoting the overall economic development of a country, in China, it is an important part of the national economy [2]. However, with the continuous increase in population, intensification of human activities, and acceleration of urbanization in coastal areas, coastal zones are facing tremendous pressure and disruption such as shortage of marine resources, deterioration of the ecological environment, rising sea level, and pollution, which seriously restrict the sustainable development of the coastal zone economy [3,4]. Due to its ability to achieve large-scale, fast and real-time dynamic observation, remote sensing technology plays an important role in monitoring the change process of land cover information in the coastal zone and studying the impact of coastal zone changes on human health, especially in densely populated areas. Therefore, it is of great value to study the changes in the resources and environment along coastlines and coastal zones and their driving factors to understand the changes in the coastal ecological environment and even global changes by means of remote sensing technology [5].
In recent years, the evolution of coastlines and coastal zones has attracted extensive attention from global experts and scholars. Most scholars have focused on monitoring the coastline via remote sensing monitoring and its dynamic changes in areas where human activities are frequent and economic strength is significant [6][7][8][9] or areas where coastal changes are significant [10]. Many scholars have used remote sensing technology to study the characteristics of erosion and accretion changes in natural coastlines [11,12]. Some scholars have studied remote sensing extraction methods for coastlines [13]. Esmail et al. [14] compared and analyzed three methods, iso clustering, threshold, and screen digitization and concluded that iso clustering is the best method of coastline extraction, but traditional methods such as visual interpretation and human-computer interactive visual interpretation are still the more common research methods used by experts and scholars to extract coastal information [15]. For remote sensing monitoring of coastal zones, different scholars have defined the scope of a coastal zone in different ways. Ma et al. [16] analyzed the land use and cover change of coastal zones in different buffer zones around the Bohai Sea, and their results help clarify the current status and changes in the surface coverage of different coastal zones. A coastal zone is representative of a fragile ecological environment, which is vulnerable to human activities. The intensity of human activities is closely related to the diversity of the coastal landscape and the change of the coastline [17]. The health of ecosystems and the marine environment of coastal zone areas are important issues affecting human survival, and studying the change of surface coverage in coastal zones can provide helpful information for solving this problem.
Using remote sensing and geographic information system technology to monitor coastline and coastal zone changes has become the topic of much scholarly research in various countries. In many studies, the combination of the two technologies is seldom discussed. One reason for that is due to the difficulty of acquiring data over a long-time span. Other reasons are that many different extraction rules have been used for extracting remote sensing information extraction rules of the coastline and coastal zone, the extraction process is complex, and the relationship between coastline and coastal zone is not always clear. Most scholars only study the coastline or coastal zone and its driving factors. Therefore, based on seven remote sensing images of Yantai from 1990 to 2018, this study uses a visual interpretation method to extract remote sensing information of the coastline as well as the 2 km coastal zone. Quantitative analysis of the Yantai coastal long time series of different spatio-temporal evolution characteristics and response relationship was conducted, and the correlation between coastline and coastal zone was also studied. The research results can provide basic data for coastal protection and coastal zone development and utilization and are of great significance for solving the health ecosystem problems of coastal zones.

Study Area
The city of Yantai is located in the northeastern part of the Shandong Peninsula, China, between 119 • 34 -121 • 57 and 36 • 16 -38 • 23 . It has a wide sea-front, bordering the Bohai Sea in the north, the Yellow Sea in the east, and the cities of Weifang, Qingdao, and Weihai. There are many small and medium-sized rivers in the city, including the Wulong, Dagujia, and Xin'an Rivers. Located north of the Tropic of Cancer, this region has a typical temperate monsoon climate. Compared with the inland areas of the same latitude, it is characterized by a mild climate, moderate rainfall, and high humidity. The annual average precipitation is approximately 524.9 mm.
The northern region of Yantai was studied, in which eight counties and municipalities are located, Laizhou, Zhaoyuan, Longkou, Penglai, Fushan, Zhifu, Laishan, and Muping. There are four ports along the coast: Yantai Port (or Zhifuwan Port), Penglai Port, Longkou Port, and Laizhou Port. The area is situated at the confluence of the Yellow Sea and the Bohai Sea, and the length of the coastline is approximately 600.40 km. The area that was studied includes the coastline and a 2-km coastal zone. Therefore, the surface area that was monitored was 795.06 km 2 . In order to ensure the integrity of the surface features in the sea area, the whole area of islands, aquaculture areas, and structures within 2 km of the coastal waters and extending outwards continuously is included in the monitoring scope ( Figure 1). By the end of 2017, the regional population reached 43.917 million, accounting for 67.13% of Yantai's population. This area had a GDP of 597.7 billion yuan, accounting for 85.11% of Yantai's GDP. This area is the port area that has the largest population and economic density in Shandong Province. Pervasive human interference from activities such as aquaculture and artificial island construction has had a noticeable impact on the changes of coastline and coastal zone in the past 30 years [18].
Appl. Sci. 2019, 9, x FOR PEER REVIEW 3 of 21 The area is situated at the confluence of the Yellow Sea and the Bohai Sea, and the length of the coastline is approximately 600.40 km. The area that was studied includes the coastline and a 2-km coastal zone. Therefore, the surface area that was monitored was 795.06 km 2 . In order to ensure the integrity of the surface features in the sea area, the whole area of islands, aquaculture areas, and structures within 2 km of the coastal waters and extending outwards continuously is included in the monitoring scope ( Figure 1). By the end of 2017, the regional population reached 43.917 million, accounting for 67.13% of Yantai's population. This area had a GDP of 597.7 billion yuan, accounting for 85.11% of Yantai's GDP. This area is the port area that has the largest population and economic density in Shandong Province. Pervasive human interference from activities such as aquaculture and artificial island construction has had a noticeable impact on the changes of coastline and coastal zone in the past 30 years [18].

Data
According to the availability of data and the economic development of the research area over the past 30 years, the data used in this study include (1) Three Landsat 5-TM images of Yantai from 1990 to 2010, with a spatial resolution of 30 m, were collected at 10-year intervals. Landsat 5-TM images with good quality covering the study area were obtained through geospatial data cloud (http://www.gscloud.cn/). Then, the tide level data of Yantai city at the corresponding time of each scene image were inquired through the Yantai tide gauge station, from which the relatively high tide level image was selected as the data source for this study. (2) The data from 2015 to 2018 are mainly based on GaoFen-2 (It was successfully launched by China on August 19, 2014) images with an accuracy of 4 m. In the key monitoring area (The key monitoring area refers to the area that can verify the correctness of the extracted information results and the area that is hard to reach by human beings), tilt photography of five-lens unmanned aerial vehicle (DJI M600 UAV) was used as the data source with an accuracy of 0. Seven remote sensing images of Yantai from 1990 to 2018 and UAV images were calibrated and registered in this study. Firstly, the GF-2 images from 2015 to 2018 were geometrically corrected. Then, the image data of other years were corrected according to the GF-2 images in Yantai in 2015. Six typical objects were selected as control points, and the relative geometric accuracy of the images from 1990 to 2010 was corrected and registered respectively to ensure that the error was within one pixel. The key step of UAV image registration was the selection of ground control points. The number of control points depended on the size of monitoring area. In this study, five control points were generally selected for correction. The images taken by UAV were set up in

Data
According to the availability of data and the economic development of the research area over the past 30 years, the data used in this study include (1) Three Landsat 5-TM images of Yantai from 1990 to 2010, with a spatial resolution of 30 m, were collected at 10-year intervals. Landsat 5-TM images with good quality covering the study area were obtained through geospatial data cloud (http://www.gscloud.cn/). Then, the tide level data of Yantai city at the corresponding time of each scene image were inquired through the Yantai tide gauge station, from which the relatively high tide level image was selected as the data source for this study. (2) The data from 2015 to 2018 are mainly based on GaoFen-2 (It was successfully launched by China on August 19, 2014) images with an accuracy of 4 m. In the key monitoring area (The key monitoring area refers to the area that can verify the correctness of the extracted information results and the area that is hard to reach by human beings), tilt photography of five-lens unmanned aerial vehicle (DJI M600 UAV) was used as the data source with an accuracy of 0. Seven remote sensing images of Yantai from 1990 to 2018 and UAV images were calibrated and registered in this study. Firstly, the GF-2 images from 2015 to 2018 were geometrically corrected. Then, the image data of other years were corrected according to the GF-2 images in Yantai in 2015. Six typical objects were selected as control points, and the relative geometric accuracy of the images from 1990 to 2010 was corrected and registered respectively to ensure that the error was within one pixel. The key step of UAV image registration was the selection of ground control points. The number Appl. Sci. 2019, 9,4466 4 of 20 of control points depended on the size of monitoring area. In this study, five control points were generally selected for correction. The images taken by UAV were set up in a unified coordinate system, and then the images were corrected and processed according to the space coordinates of different ground control points.

Coastline Extraction
The coastline is the demarcation line between ocean and land, and it is the link of the average high tide trace [19]. There are many methods to extract coastline data by remote sensing, such as man-machine interactive visual interpretation [20], the normalized water body index [21], and the improved water boundary method [22]. However, considering the large time span, different resolution of image data, and the complex and diverse types of coastline, in this study, we use a visual interpretation method to draw and extract the coastline every 10 years from 1990 to 2018. Due to the decrease of coastal zone variation in recent years, we also extracted the coastline data of the research area from 2015 to 2018, which provides us with high resolution image data for the past four years. The extraction of different types of coastlines is based on defined principles and methods. In order to minimize the human errors in the interpretation of the coastline by different researchers, all of the years of coastline data extraction were completed by the same person. In order to ensure its accuracy, field survey and UAV photography were used to verify the data. First, the coastline types which are difficult to distinguish in remote sensing images are investigated and validated on the spot. Second, a random selection of several areas, using UAV data and remote sensing image data comparative analysis, to verify the accuracy of coastline extraction.
Yantai's coastline can be divided into three categories: natural, artificial, and estuarine. On the basis of the first class, the second-class types were classified. The interpretation signs and basic norms of different types of coastlines mainly refer to the contents and indicators of basic geographic national conditions monitoring in 2018. Detailed classification information is shown in Table 1. Table 1. The coastline classification system of Yantai.

The First Type The Secondary Type Classification Definition
Gravelly boundary line The intertidal sediment of sandy coastline is mainly gravel, which is a relatively straight coastline formed by sediment such as sand and gravel under the long-term action of waves. Natural boundary line Silty boundary line A sea-land boundary consisting of silt or silty mudflat materials. Bedrock boundary line It is the demarcation line between land and sea on the bedrock coast. Breeding boundary line A land-sea boundary constructed artificially for aquaculture.

Salt field boundary
The sea-land dividing line for saline-alkali drying. Port and wharf boundary line The demarcation line between land and sea for the functional use of port terminals.
Artificial boundary line Construction of coastline near the sea The demarcation line between land and sea for the development of coastal industry and other construction purposes.

Road boundary line
The land-sea dividing line of the road constructed artificially.

Protective engineering coastline
It refers to the boundary line formed by the tidal dikes, wave walls and slope protection structures built to prevent the erosion of sea waves and coastal currents and the invasion of natural disasters such as typhoons, cyclones, and cold waves and gales. Other artificial boundary lines Artificial boundary lines that do not belong to the above classification.

Estuary boundary line
The boundary between estuary and ocean, and between river and ocean.
Different types of coastlines are extracted from standard pseudo-color synthetic images taken by satellite remote sensing [23]. Among them, the natural boundary line is the trace line of the average high tide line formed naturally by natural land-sea interaction without human disturbance [24], including the gravel-sandy boundary line, silty boundary line, and bedrock boundary line. The gravel-sandy boundary line is relatively straight, and the boundary line between land and sea is generally located on the top of the ridge to the seaside. Images show a distinct difference between wet and dry due to the influence of tidal water, and the boundary line is located at the dry-wet boundary. The silty boundary line is mainly located on the land-facing side of the broadband where the growth of salt-tolerant plants has changed significantly over the intertidal zone. The acquisition position of the bedrock boundary line is defined at the base of a cliff or at the junction of land and water. The definition of an artificial boundary line is relatively easy, it mainly refers to the sea-land boundary line divided by the periphery of artificial facilities, such as damp-proof dikes, slope protection, wharfs, roads, aquaculture ponds, salt fields, or other water-retaining structures. An estuary boundary line is the demarcation line between an estuary and the ocean. Figure 2 shows the results of various types of coastline extraction methods.
generally located on the top of the ridge to the seaside. Images show a distinct difference between wet and dry due to the influence of tidal water, and the boundary line is located at the dry-wet boundary. The silty boundary line is mainly located on the land-facing side of the broadband where the growth of salt-tolerant plants has changed significantly over the intertidal zone. The acquisition position of the bedrock boundary line is defined at the base of a cliff or at the junction of land and water. The definition of an artificial boundary line is relatively easy, it mainly refers to the sea-land boundary line divided by the periphery of artificial facilities, such as damp-proof dikes, slope protection, wharfs, roads, aquaculture ponds, salt fields, or other water-retaining structures. An estuary boundary line is the demarcation line between an estuary and the ocean. Figure 2 shows the results of various types of coastline extraction methods.

Calculation of Coastline Change
In this study, the end point change rate (EPR) and linear regression rate (LRR) of the digital coastline analysis system module developed by the U.S. Geological Survey are used to analyze the coastline change rate in northern Yantai [25]. The main steps are as follows. First, based on the coastline data of the past years, a baseline is extracted from the land or seaside by establishing a buffer zone, factor transfer line, and adjustment. Then, all of the annual coastline data is integrated into one factor class, that is, all coastlines. Second, equidistant tangents perpendicular to the above baselines and intersecting with all coastlines are made along the baselines. Then, the DSAS system is used to calculate all sections, and the coastline transition rate can be obtained by calculation. The

Calculation of Coastline Change
In this study, the end point change rate (EPR) and linear regression rate (LRR) of the digital coastline analysis system module developed by the U.S. Geological Survey are used to analyze the coastline change rate in northern Yantai [25]. The main steps are as follows. First, based on the coastline data of the past years, a baseline is extracted from the land or seaside by establishing a buffer zone, factor transfer line, and adjustment. Then, all of the annual coastline data is integrated into one factor class, that is, all coastlines. Second, equidistant tangents perpendicular to the above baselines and intersecting with all coastlines are made along the baselines. Then, the DSAS system is used to calculate all sections, and the coastline transition rate can be obtained by calculation. The terminal change rate is used to analyze the change rate of two coastlines and several coastlines in different periods [26]. The formulas are as follows: where E (i,j) is the rate of change of coastline endpoint from the beginning to the end of the study, namely EPR. i is the coastline time at the beginning of the study, j is the coastline time at the end of the study, ∆Y (j,i) is the time interval between the beginning of the study and the end of the study, d j and d i are the distance from the vertical direction of the coastline to the baseline when the time is i and j, respectively. The linear regression rate (LRR) is better than the terminal change rate in studying the changes in multiple coastlines [27]. In this study, LRR analyses the proportional relationship and changes of multiple coastlines at different times, and its mathematical expression are as follows: where y is a dependent variable, which is the spatial position of the coastline, x is the independent variable of the year, a is the constant intercept of the fitting, b is the regression slope, which represents the y change corresponding to each unit x change, namely LRR. In order to clearly see the variation degree of coastline length for different time periods, the following formula is used to express the annual variation information of coastline length, that is, the coastline variation intensity [28]: where LCI ij is the variation intensity of coastline length in the study area from year i to year j, and L i and L j are the coastline lengths in the study area in year i and year j, respectively.

Calculation of Fractal Dimension of the Coastline
In this study, the fractal dimension of the coastline is calculated by the grid method [29]. The basic principle of the grid method is to use square grids of different lengths to cover the measured coastline continuously and with no overlap. When the length of square grids r takes different values, the number of grids N(r) needed to cover the whole coastline will change accordingly [28]. The formula is as follows: By taking the natural logarithm of both sides of equation (6), the following results are obtained: where D is the fractal dimension of the measured coastline, and C is a constant. The fractal dimension D can be calculated by fitting and regression analysis of different R values and corresponding N(r) values.
The larger the value of the fractal dimension, the more complicated and tortuous the coastline. On the contrary, the smaller the value of the fractal dimension, the smaller the tortuosity and complexity of the coastline [30]. Based on the ArcGIS platform and least squares method, the fractal dimension of the coastline is obtained by calculations followed by fitting as well as regression analysis.

Extraction of Coastal Zone Coverage Information
Based on the 2018 coastline, this study defines the coastal zone as the zonal area formed by extending 2 km from coastline to inland and uses the same visual interpretation method as the extraction of the coastline to extract the surface coverage information of the coastal zone. Referring to the content and indicators of China's geographic survey, this study highlights the types and elements of coastal zone development and utilization and forms a classification system of land cover for coastal zone development and utilization change monitoring. The surface coverage of the coastal zone in the study area includes eight first-class classifications: cultivated land, forests/grasses coverage, building areas, roads, structures, manually excavated land, desert/bare land, water areas [31]. Due to the difference in image resolution and subjective error of visual interpretation in different years, this study only considers the change of the first-class coastal zone in the study area and does not analyze the second-class land cover types. In this paper, UAV is used to verify the correctness of information extraction results, and it can also reach places that are difficult for human beings to reach for monitoring ( Figure 3). The classification system of first-class land cover is consistent with the classification standard of Shandong Province's geographic condition monitoring (Table 2). different years, this study only considers the change of the first-class coastal zone in the study area and does not analyze the second-class land cover types. In this paper, UAV is used to verify the correctness of information extraction results, and it can also reach places that are difficult for human beings to reach for monitoring ( Figure 3). The classification system of first-class land cover is consistent with the classification standard of Shandong Province's geographic condition monitoring ( Table 2).

The first type Classification definition Cultivated land
Land reclaimed, cultivated and managed regularly.
Forests/gresses cover It mainly includes mangroves, other woodlands, and grasslands.
Building areas Including housing construction area, independent housing construction and abandoned housing construction.

Roads
It includes railways, highways, urban roads, and rural roads.

Structures
An engineering entity or ancillary building facility built for a purpose of use, in which production and living activities are not carried out directly within it(GB/T 50504-2009). Including hardened surface, industrial facilities, salt ponds, dykes, and other structures.

Manually excavated land
Ground surface covered by human waste or exposed by human excavation for a long time. Desert/bare land Desert and bare surface and low coverage grassland.

The First Type Classification Definition
Cultivated land Land reclaimed, cultivated and managed regularly. Forests/gresses cover It mainly includes mangroves, other woodlands, and grasslands.
Building areas Including housing construction area, independent housing construction and abandoned housing construction.

Roads
It includes railways, highways, urban roads, and rural roads.

Structures
An engineering entity or ancillary building facility built for a purpose of use, in which production and living activities are not carried out directly within it(GB/T 50504-2009). Including hardened surface, industrial facilities, salt ponds, dykes, and other structures.
Manually excavated land Ground surface covered by human waste or exposed by human excavation for a long time. Desert/bare land Desert and bare surface and low coverage grassland.

Waters
It refers to the spatial range of the water body's growth and declines over a long period of time.

Variation of Coastline Length
According to the seven remote sensing images of the study area from 1990 to 2018, the total length of the coastline in each period was obtained by visual interpretation (Figure 4). In order to ensure the accuracy of the extracted coastline data, GF-2 images with a resolution of 4 m were used in this study. The results of coastline extraction were verified by UAV photography, and the accuracy of coastline extraction was calculated. The results showed that the typing accuracy of the coastline extracted in this study reached 97.02% and the length accuracy of the coastline reached 94.78%. In the chart, it is clear that 2015 is a time node for the change of coastline length. Before 2015, the port areas headed by Penglai, Longkou, and Zhifu developed rapidly. The coastline length showed a continuous growth trend, that is, it increased by 174.15 km, an increase of 41.53%. In 2015, the relevant departments began to take effective measures to repair the islands and coastlines in view of the excessive exploitation and utilization of the ocean and coastal zones and the impact of natural disasters, such as typhoons and ocean storm surges [32]. The length of the coastline showed a relatively stable small growth trend with a cumulative increase of 6.93 km, an increase of approximately 1.16% in 2015.
in this study. The results of coastline extraction were verified by UAV photography, and the accuracy of coastline extraction was calculated. The results showed that the typing accuracy of the coastline extracted in this study reached 97.02% and the length accuracy of the coastline reached 94.78%. In the chart, it is clear that 2015 is a time node for the change of coastline length. Before 2015, the port areas headed by Penglai, Longkou, and Zhifu developed rapidly. The coastline length showed a continuous growth trend, that is, it increased by 174.15 km, an increase of 41.53%. In 2015, the relevant departments began to take effective measures to repair the islands and coastlines in view of the excessive exploitation and utilization of the ocean and coastal zones and the impact of natural disasters, such as typhoons and ocean storm surges [32]. The length of the coastline showed a relatively stable small growth trend with a cumulative increase of 6.93 km, an increase of approximately 1.16% in 2015. According to the trend map of coastline length change in the study area from 1990 to 2018, the coastlines of several periods were selected randomly to calculate the intensity of coastline change (Table 3). In the past 30 years, the changing intensity of coastline length in the study area is approximately 1.54%. The changing intensity of coastline length in the study area was the greatest during 2010-2015 with a value of 3.47%. The changing intensity of coastline length in the study area was the smallest during 2015-2018 with a value of 0.39%. According to the trend map of coastline length change in the study area from 1990 to 2018, the coastlines of several periods were selected randomly to calculate the intensity of coastline change (Table 3). In the past 30 years, the changing intensity of coastline length in the study area is approximately 1.54%. The changing intensity of coastline length in the study area was the greatest during 2010-2015 with a value of 3.47%. The changing intensity of coastline length in the study area was the smallest during 2015-2018 with a value of 0.39%.

Spatial and Temporal Variation Characteristics of the Coastline
The end point change rate (EPR) and linear regression rate (LRR) were used to study and analyze the coastline changes of counties and municipalities in the northern region of Yantai from 1990 to 2018. Figures 6 and 7 show that the coastline generally shows a trend of growth extending to the sea, but the trend of coastline changes in different counties and municipalities is different. The coastline changes in Laizhou, Longkou, Penglai, and Zhifu are more noticeable.
From 1990 to 2018, the average annual variation EPR and LRR of the coastal areas in northern Yantai were 6.34 m/a and 7.34 m/a, respectively, and the growth of coastal sediment was dominated by seaward siltation. The annual variation rate of the coastline in each county and municipality ranged from large to small in the following order: Laizhou > Penglai > Longkou > Zhifu > Fushan > Laishan > Mouping > Zhaoyuan. Among them, the annual variation rate of the Laizhou coastline was the largest an EPR of 33.67 m/a and an LRR of 39.37 m/a. The maximum deposition rate of the coastline in the study area also appeared in the coastal section and is located around the northern aquaculture area of Tushan Town, Laizhou with an EPR of 168.88 m/a and an LRR of 200.48 m/a, and the maximum erosion rate of the coastal section was −6.37 m/a (EPR) and −5.83 m/a (LRR). The Zhaoyuan coastline had the smallest annual variation rate with an EPR of −1.91 m/a and an LRR of −1.37 m/a. The maximum deposition rate of the Zhaoyuan coastline was 7.73 m/a (EPR) and 9.56 m/a (LRR), and the maximum erosion rate was −4.93 m/a (EPR) and −4.67 m/a (LRR). The maximum erosion rates of the coastline in the study area were −43.27m/a (EPR) and −47.41 m/a (LRR). In the past 30 years, the coastline of Laizhou, Penglai, Longkou, and Zhifu has become the fastest-changing area of the coastline of the Jiaodong Peninsula due to the continuous influence of human activities, such as reclamation, artificial island construction, dam construction, and port construction.
Human activities are irreversible and destructive to coastline changes, and it is difficult to restore the original natural coastline [33]. From the changes in the coastline of different regions, we can see that human activities have the greatest impact on the changes in the coastline. Especially from 2010 to 2015, the flooding of the Laizhou aquaculture area, construction of the Longkou artificial island, coastal construction of Penglai and Zhifu, and expansion of port terminals have made great changes in the coastline. However, since 2015, due to the public's awareness of a series of problems caused by the destruction of the coastline and reduction of marine resources, various departments have strengthened the control of coastal areas, increasing the control coastline changes and protection of the marine ecological environment.

Characteristics of the Fractal Dimension of Coastline Variation
The shape of the coastline is winding and complex. Therefore, the calculation of the fractal dimension of the coastline is helpful to understand the change of the length and shape of the coastline and can provide an important basis for promoting the protection of the coastline and its manageable development and utilization. Based on the coastline information extracted by visual interpretation from 1990 to 2018, the fractal dimension changes of the coastline in the study area were calculated by ArcMap software and the least squares method. In this study, 10 grid lengths (30,60,90,120,150,180,210,240,270, and 300 m) were used to cover the coastline to be analyzed, and the fractal dimension of the coastline was calculated. Finally, the changes in the fractal dimension of the coastline in the study area in the past 30 years were obtained (Figure 8). Appl. Sci. 2019, 9, x FOR PEER REVIEW 12 of 21   From 1990 to 2018, the average annual variation EPR and LRR of the coastal areas in northern Yantai were 6.34 m/a and 7.34 m/a, respectively, and the growth of coastal sediment was dominated by seaward siltation. The annual variation rate of the coastline in each county and municipality ranged from large to small in the following order: Laizhou > Penglai > Longkou > Zhifu > Fushan > Laishan > Mouping > Zhaoyuan. Among them, the annual variation rate of the Laizhou coastline was the largest an EPR of 33.67 m/a and an LRR of 39.37 m/a. The maximum deposition rate of the coastline in the study area also appeared in the coastal section and is located around the northern aquaculture area of Tushan Town, Laizhou with an EPR of 168.88 m/a and an LRR of 200.48 m/a, and the maximum erosion rate of the coastal section was -6.37 m/a (EPR) and -5.83 m/a (LRR). The Zhaoyuan coastline had the smallest annual variation rate with an EPR of -1.91 m/a and an LRR of -1.37 m/a. The maximum deposition rate of the Zhaoyuan coastline was 7.73 m/a (EPR) and 9.56 m/a The fractal dimension of the study area calculated in this study is close to that of the coastline of Shandong Province studied by Xu Ning [27], and the trend of change is consistent. Generally speaking, the fractal dimension of the coastline in the study area shows an increasing trend with time. However, it can be clearly seen in Figure 7 that before 2015, the coastline fractal dimension increased by 0.0382, but in 2015-2018, the coastline fractal dimension increased by 0.0043. Thus, as people have gradually realized that the destruction of human activities on the oceans has begun to have a serious impact on human survival, the protection of marine resources and careful use of marine resources has become a greater priority [10]. manageable development and utilization. Based on the coastline information extracted by visual interpretation from 1990 to 2018, the fractal dimension changes of the coastline in the study area were calculated by ArcMap software and the least squares method. In this study, 10 grid lengths (30,60,90,120,150,180,210,240,270, and 300 m) were used to cover the coastline to be analyzed, and the fractal dimension of the coastline was calculated. Finally, the changes in the fractal dimension of the coastline in the study area in the past 30 years were obtained. (Figure 8). The fractal dimension of the study area calculated in this study is close to that of the coastline of Shandong Province studied by Xu Ning [27], and the trend of change is consistent. Generally speaking, the fractal dimension of the coastline in the study area shows an increasing trend with time. However, it can be clearly seen in Figure 7 that before 2015, the coastline fractal dimension increased by 0.0382, but in 2015-2018, the coastline fractal dimension increased by 0.0043. Thus, as people have gradually realized that the destruction of human activities on the oceans has begun to have a serious impact on human survival, the protection of marine resources and careful use of marine resources has become a greater priority [10].

Analysis of Surface Cover Change in the Coastal Zone
In 2018, the coastal land cover in the northern Yantai coastal area was dominated by forests/ grasses cover and structures, accounting for 31.40% and 15.45% of the study area, respectively, while other land types accounted for a relatively small proportion. From 1990 to 2018, the area of forests/grasses covers, and structures increased, while the area of desert/bare land and water area decreased. Among them, forests/grasses cover area increased from 153.56 km 2 to 249.65 km 2 , structure area increased from 9.88 km 2 to 122.84 km 2 , desert/bare land area decreased from 129.46 km 2 to 37.20 km 2 , and water area decreased from 93.15 km 2 to 80.83 km 2 (Figure 9) Appl. Sci. 2019, 9, x FOR PEER REVIEW 15 of 21

Analysis of Surface Cover Change in the Coastal Zone
In 2018, the coastal land cover in the northern Yantai coastal area was dominated by forests/ grasses cover and structures, accounting for 31.40% and 15.45% of the study area, respectively, while other land types accounted for a relatively small proportion. From 1990 to 2018, the area of forests/grasses covers, and structures increased, while the area of desert/bare land and water area decreased. Among them, forests/grasses cover area increased from 153.56 km 2 to 249.65 km 2 , structure area increased from 9.88 km 2 to 122.84 km 2 , desert/bare land area decreased from 129.46 km 2 to 37.20 km 2 , and water area decreased from 93.15 km 2 to 80.83 km 2 . (Figure 9) Hou Xiyong's research reported that the coastal land area of Shandong Province is continuously increasing [34]. The eight counties and municipalities involved in the northern coastal area of Yantai studied in this study are all located in the northern part of Shandong Province, and the coastal land area has been increasing over the past 30 years. From 1990 to 2018, the coastal areas of Longkou and Laishan increased by 52.54% and 34.63%, while those of Fushan and Mouping increased by 8.05% Hou Xiyong's research reported that the coastal land area of Shandong Province is continuously increasing [34]. The eight counties and municipalities involved in the northern coastal area of Yantai studied in this study are all located in the northern part of Shandong Province, and the coastal land area has been increasing over the past 30 years. From 1990 to 2018, the coastal areas of Longkou and Laishan increased by 52.54% and 34.63%, while those of Fushan and Mouping increased by 8.05% and 9.04%. Table 6 shows the change of surface coverage of the first-class coastal zones in each county and municipality from 1990 to 2018. The study found that from 1990 to 2018, the area of structures in the study area increased the most, to 112.96 km 2 , of which the area of structures in Laizhou increased by 50.65 km 2 , accounting for 44.84%. The area of desert/bare land decreased the most, by 92.26 km 2 , and Penglai City decreased the most, by 48.46 km 2 , accounting for 52.53%. In contrast, the changes in building areas and water areas are relatively small, i.e., 5.87 km 2 and 12.32 km 2 , respectively. Based on the vector data of the coastal zone extracted by remote sensing, the transfer matrix analysis of various types of coastal zone areas in the study area from 2015 to 2018 was carried out to fully understand the land pattern transformation of the coastal zone in the study area in recent years. As can be seen from Table 7, from 1990 to 2018, forests/grasses cover and structure transfer area were the largest, 139.92 km 2 and 70.58 km 2 , respectively. The increased area of forests/grasses cover was mainly transformed from cultivated land and desert/bare land. The increased area of structure was mainly transformed from water area, desert/bare land, and building area. Desert/bare land has the largest area of transfer, of which 35.72 km 2 is converted into forests/grasses cover, and then into cultivated land, building areas, and structures. In addition, the change of manually excavated land is not clear, and the increase only accounts for 2.70% of the total change. Generally speaking, the region presents a state of transition from desert/bare land, cultivated land, and building areas to forests/grasses cover and structures, but the overall area of the study area experienced little change and is relatively stable.

Relationship between coastline change and coastal surface cover change
The change of land cover types in offshore areas will have a certain impact on the change of coastline types and lengths. Through comparative analysis, it was found that the changes of coastal structures, manually excavated land, and roads in the study area are related to the length of the artificial coastline. From 1990 to 2018, the length of the coastline in the study area increased with the use of structures, manually excavated land and roads as the coastline, and the length of the artificial boundary line in the study area also increased. The correlation is shown in Figure 10.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 17 of 21 use of structures, manually excavated land and roads as the coastline, and the length of the artificial boundary line in the study area also increased. The correlation is shown in Figure 10. Human activities have brought about significant changes in resources and to the environment of the coastal zone, which have not only shaped many types of artificial coastline landforms, such as harbor and wharf, protection engineering, and coastal construction, but they also have had a great impact on the natural evolution process of the coastal zone. The coastline of Yantai is clearly affected by the change of land use in the coastal zone, especially by structures, manually excavated land, and roads. Through comparative analysis of the impact of structures, manually excavated land, roads, and desert/bare land on the coastline from 1990 to 2018, we determined that the artificial boundary line in the northern coastal area of Yantai is directly proportional to the length of the coastline with roads as the coastline, which was the most significant correlation (R 2 = 0.9527, P < 0.01). It can be seen that in the past 30 years, roads have had the greatest impact on the coastline changes in the northern coastal areas of Yantai, and the impact of structures and manually excavated land on the coastline changes is more significant, which is the impact of frequent human activities on the coastline.
Yantai has a superior geographical location and abundant marine resources, and its coastal zone has a long history of development. In recent years, the coastal regional economy of northern Yantai has developed rapidly. The port cities headed by Zhifu, Penglai, and Laizhou have been growing steadily. The construction of the Longkou artificial island group has increased the driving force for the overall economic development of Yantai [18]. However, at the same time, people need to be aware of the series of health problems brought about by the overdevelopment of coastal areas. How to minimize the damage while undergoing developing has become an urgent problem that needs to be solved.

Discussion
The northern coastal area of Yantai is the key development area of the Blue Economic Zone in the Shandong Peninsula. The sensible development and utilization of coastal areas and marine resources is of great significance to the economic development of this area [35]. In the past 30 years, Human activities have brought about significant changes in resources and to the environment of the coastal zone, which have not only shaped many types of artificial coastline landforms, such as harbor and wharf, protection engineering, and coastal construction, but they also have had a great impact on the natural evolution process of the coastal zone. The coastline of Yantai is clearly affected by the change of land use in the coastal zone, especially by structures, manually excavated land, and roads. Through comparative analysis of the impact of structures, manually excavated land, roads, and desert/bare land on the coastline from 1990 to 2018, we determined that the artificial boundary line in the northern coastal area of Yantai is directly proportional to the length of the coastline with roads as the coastline, which was the most significant correlation (R 2 = 0.9527, P < 0.01). It can be seen that in the past 30 years, roads have had the greatest impact on the coastline changes in the northern coastal areas of Yantai, and the impact of structures and manually excavated land on the coastline changes is more significant, which is the impact of frequent human activities on the coastline.
Yantai has a superior geographical location and abundant marine resources, and its coastal zone has a long history of development. In recent years, the coastal regional economy of northern Yantai has developed rapidly. The port cities headed by Zhifu, Penglai, and Laizhou have been growing steadily. The construction of the Longkou artificial island group has increased the driving force for the overall economic development of Yantai [18]. However, at the same time, people need to be aware of the series of health problems brought about by the overdevelopment of coastal areas. How to minimize the damage while undergoing developing has become an urgent problem that needs to be solved.

Discussion
The northern coastal area of Yantai is the key development area of the Blue Economic Zone in the Shandong Peninsula. The sensible development and utilization of coastal areas and marine resources is of great significance to the economic development of this area [35]. In the past 30 years, the economic development of Yantai City has been very fast, especially in the northern coastal areas. Due to the impact of human activities on the coastal zone, the coastal zone of the study area has undergone significant changes, which can be confirmed by the interpretation and analysis of remote sensing images.
Coastal zone information was extracted every 10 years between 1990 and 2018. The analysis shows that the increase range of coastline length changed in the following order: 1990-2000 < 2000-2010 < 2010-2018. Among them, in the period of 2010-2018, the coastline length experienced the greatest change. With the development of the economy, the exploitation and utilization of the coastal zone have been intensified, and the coastline has shown a rapidly increasing trend. However, the analysis results show that the length of the coastline has not been increasing rapidly during this period. During the period from 2015 to 2018, the growth rate of the coastline decreased significantly. Since 2015, in view of a series of problems in coastal areas, relevant departments have begun to take measures to rectify and repair damaged islands and coastal zones, and the over-exploitation of coastal zones has been alleviated to a certain extent [36].
According to Figures 4 and 8, the change of the fractal dimension of the coastline in the study area is consistent with that of its length, first showing a trend of rapid increase and then slowly increasing. There are two main reasons for this trend. One is because of the low resolution of remote sensing images from 1990 to 2010 and the rough coastal information extracted, the coastline length and fractal dimension values are smaller. The resolution of remote sensing images in 2015 and beyond is higher, which can extract coastal information in more detail, and the length and fractal dimension of the extracted coastline are more accurate. The second reason is that with the rapid development of the social economy, especially the economic development of the coastal zone, people's utilization of the ocean has been strengthened. The construction of reclamation, artificial islands and port wharfs have increased the length of the coastline and made it more tortuous, and the fractal dimension value of coastline increases with the length. Since 2015, the coastline has changed slightly, which is because various departments began to attach importance to the development of coastal areas and actively take measures to control the development and utilization of coastal areas and effectively protect marine resources. This awareness of coastlines being adversely impacted by multiple driving forces has accelerated our efforts to assess, monitor and mitigate coastal stressors [37]. The coastal environmental quality has been greatly improved.
In recent years, the artificial exploitation of the coastal zone has been increasing, the carrying capacity of coastal area has become fragile, and the environment on which human beings depend has been challenged. Local government should attach importance to the development of coastal areas and actively take measures to control the development and utilization of coastal areas and effectively protect marine resources. Effective planning and management are the preconditions for sustainable coastal development [38]. With Yantai becoming the backbone city of the Blue Economic Zone in Shandong Peninsula, the land use types in northern coastal areas are continuously changing. It is necessary to develop the economy, protect marine ecological resources and the environment, pay attention to the development status of the coastal zone, and rationally develop and plan coastal areas for the sustainable development of the marine economy in the future. It is also an important precondition to solve the problems of the health ecosystem of coastal zones. Therefore, we should not neglect the value of natural resources for human beings in pursuit of personal interests. Economic development should be based on the balance of ecological resources and the environment. Otherwise, when the economy develops to a certain extent, human beings will face greater survival problems.

Conclusions
In this study, the coastal zone of the Yantai northern region was extracted by remote sensing every 10 years from 1990 to 2018. The coastal zone information of the study area from 2015 to 2018, which can obtain high resolution image data in the past four years, has been extracted more precisely. Then, the changes in the coastline and coastal zones of eight counties and municipalities were analyzed.
(1) In the past 40 years, the spatial and temporal differences of the coastline in the study area are significant: from 1990 to 2018, the length of coastline increased greatly. The increased areas were mainly concentrated in Laizhou, Longkou, Penglai, and Zhifu. The increased coastline types were mainly artificial boundary lines. The length of the natural boundary line decreased, and a large number of natural boundary lines were transformed into artificial boundary lines. Human activities have seriously affected the development of the coastline. Among them, the range of coastline change in the study area decreased significantly from 2015 to 2018. (2) Human activities affect not only the coastline changes but also the types and areas of surface coverage in the coastal zone. In recent years, with the rapid development of the regional economy and urbanization in the northern coastal areas of Yantai, the types of land surface coverage in the coastal zone are mainly transferred from desert/bare land and cultivated land to structure and forests/grasses cover. Influenced by human activities, the monitoring area of various types is increasing, especially the area of structures, manually excavated land, and water area, then health problems in coastal zones increase with changes in coastal zones. (3) The changes of resources and environment in the coastal zone will affect the types and areas of surface coverage, and then affect the types of coastline landforms, so that the types of coastline landforms will change with the changes of the types of surface coverage near the coastline, especially the types of structures and roads constructed to the sea.