Determination of the Original Coastline on Semarang City and Demak District Using Remote Sensing Approach

E3S Web of Conferences 202, 04001 (2020) https://doi.org/10.1051/e3sconf/202020204001 ICENIS 2020

Determination of the original coastline on Semarang city and Demak district using remote sensing approach

Ignatius Sriyana1*, John Bosco Niyomukiza1,2, Sri Sangkawati1 and Sri Kumala Parahyangsari3 1Department of Civil Engineering, Faculty of Engineering, Diponegoro University, Semarang, Indonesia 2Department of Civil Engineering, Faculty of Engineering, Ndejje University, Kampala, Uganda 3Member of Watershed Coordination Forum in Central Java Province, Central Java, Indonesia

Abstract. Development of the coastal area without environmental consciousness is disrupting coastline equilibrium, thus triggering erosion and accretion. This results in the obscurity of the border area due to coastline change. The objective of this study is to determine the original coastline of Semarang city and Demak district. The method used is remote sensing using Landsat satellite imagery data dating from 1984 to 2018. The results showed that the length of the coastline in Semarang and Demak border coast shows a continuous increase from 48.2, 57.9, 63.7, and 73.5 km in 1984, 1992, 2000, and 2008, respectively to 92.6 km in 2018, and the average annual change rate from 1984 to 1992, 1992 to 2000, 2000 to 2008, and 2008 to 2018 is 1.2 km per year, 0.7 km per year, 1.2 km per year, and 2.1 km per year, respectively. It was concluded that the current coastline condition of Semarang city and Demak district has had an increase of 29.8 km in length and shifted inwards by 0.5-2 km in 2000 – 2018 compared to the original coastline. This shows a deviation of 5-15% from Government Regional Spatial Plan. The study recommends a sustainable coastal management and action plan to be prepared by the government, as it is vital in preserving the coastal regions.

1 Introduction The coastline is one of the vital linear features on the earth’s crust, and due to its dynamic and complex nature, it is sensitive to erosion and deposition [1-3]. It is, therefore, paramount to detect the coastal changes that could occur during specified periods. The disrupted coastal balance makes the coastline undergo frequent changes leading to both erosion and accretion [4]. By definition, a coastline is a boundary between the land and the water surface [1], [5- 7], and it tends to retreat over 70% of the population across the globe that encounter coastal erosion [8]. Changes in coastline sometimes occur for a prolonged period due to geological occurrences. The issues of coastline change are continuously becoming a significant concern

* Corresponding author: [email protected]

© The Authors, published by EDP Sciences. This is an open access article distributed under the terms of the Creative Commons Attribution License 4.0 (http://creativecommons.org/licenses/by/4.0/). E3S Web of Conferences 202, 04001 (2020) https://doi.org/10.1051/e3sconf/202020204001 ICENIS 2020

to the environment, social and economic aspects across the globe [4, 7]. For that matter, the knowledge about coastline changes is essential in coastal management, more especially to countries with many islands like Indonesia. Absence of information about coastline changes lead to many disasters including negative effects to the surrounding environment [9]. Coastline change usually occurs due to different factors ranging from human activities to natural factors [10]. The natural factors include geology, weather condition, slope angle, waves and strength of the wind, and the amount of vegetation, whereas human factors include constructions, industrialization, mining, urbanization, and any other activity that intercepts free passage of sediments along the coast. According to Di et al. [11], it is crucial to detect and map the coastline, as it is important for safe navigation, sustainable coastal development and planning, coastal flood forecasting, resource management, and environmental protection [12-14]. Besides, it is important to determine the changes that occurred during different years for comparison purposes and help acquiring the information to forecast the future coastline changes. Previous studies have discussed about coastline changes, for example the study that was conducted in China showed that coastline change caused the Pearl River Estuary straitened and caused a lot of floods in the area [15]. The study that was done by Hussain et al. [16] showed coastline changes in the eastern part of Meghna Estuary, Bangladesh with a high magnitude and the major cause was sediment discharge from rivers. Alam and Kabir [17] examined morphological variations caused by erosion versus accretion in the coastal areas and offshore island of Bangladesh, and their study revealed that the areas were affected by accretion more than erosion during 1977 – 2010. The cause of this phenomenon is the rising sea level that created accumulation of sediments. Tamassoki et al. [9] monitored shoreline changes in coastal city of Bandar Abbas using Landsat imaginary dated from 1984 to 2009, and noticed a continuous increase in the shoreline length from 272 km to 330 Km in 1984 to 2009, respectively. Their study revealed that the changes were caused by geological and hydrodynamic occurrences. There are different methods of detecting the changes in the coastline, and among them include remote sensing approach, histogram thresholding, ground survey, aerial photography, and other methods [6, 18]. However, remote sensing has proved to be the most effective method to assess the changes in coastal areas [19]. A lot of researches has been conducted using remote sensing to monitor the coastline changes affected by erosion and accretion [1, 7, 20, 21]. Their results yielded success in the areas where they monitored coastlines using remote sensing. According to Hussain et al. [16], remote sensing provides the best solution as far as acquiring and processing data in a quick and faster way compared to other methods like photogrammetry. The objective of this paper is to determine the original coastline of Semarang city and Demak district in order to determine the original border area. This information is essential to the government because knowing the coastline change will enable them to understand the changing of physical environment better, thus developing coastal policies that will promote sustainable management practices. Besides, knowing the location and coastline changes is helpful to engineers while designing coastal protection structures. To achieve the objective of the current study, we used remote sensing approach using Landsat satellite imagery data dating from 1984 to 2018.

2 E3S Web of Conferences 202, 04001 (2020) https://doi.org/10.1051/e3sconf/202020204001 ICENIS 2020 to the environment, social and economic aspects across the globe [4, 7]. For that matter, the 2 Materials and methods knowledge about coastline changes is essential in coastal management, more especially to countries with many islands like Indonesia. Absence of information about coastline changes lead to many disasters including negative effects to the surrounding environment [9]. 2.1 Data and location Coastline change usually occurs due to different factors ranging from human activities to The data needed in this research include; coastline length, area of abrasion and accretion that natural factors [10]. The natural factors include geology, weather condition, slope angle, occurred during different years. The research was conducted along the coast with a length of waves and strength of the wind, and the amount of vegetation, whereas human factors include 48.2 km at Demak district northern coastal area, starting from the coastline in Semarang city constructions, industrialization, mining, urbanization, and any other activity that intercepts located at coordinate 6°55'58.81"S, 110°27'53.20"E and ending at coordinate 6°42'19.98"S, free passage of sediments along the coast. According to Di et al. [11], it is crucial to detect 110°37'2.18"E, as presented in Figure 1. This region marked the northern part of the waters and map the coastline, as it is important for safe navigation, sustainable coastal development of Java Island. and planning, coastal flood forecasting, resource management, and environmental protection [12-14]. Besides, it is important to determine the changes that occurred during different years for comparison purposes and help acquiring the information to forecast the future coastline changes. Previous studies have discussed about coastline changes, for example the study that was conducted in China showed that coastline change caused the Pearl River Estuary straitened and caused a lot of floods in the area [15]. The study that was done by Hussain et al. [16] showed coastline changes in the eastern part of Meghna Estuary, Bangladesh with a high magnitude and the major cause was sediment discharge from rivers. Alam and Kabir [17] examined morphological variations caused by erosion versus accretion in the coastal areas and offshore island of Bangladesh, and their study revealed that the areas were affected by accretion more than erosion during 1977 – 2010. The cause of this phenomenon is the rising sea level that created accumulation of sediments. Tamassoki et al. [9] monitored shoreline changes in coastal city of Bandar Abbas using Landsat imaginary dated from 1984 to 2009, and noticed a continuous increase in the shoreline length from 272 km to 330 Km in 1984 to 2009, respectively. Their study revealed that the changes were caused by geological and hydrodynamic occurrences. There are different methods of detecting the changes in the coastline, and among them include remote sensing approach, histogram thresholding, ground survey, aerial photography, and other methods [6, 18]. However, remote sensing has proved to be the most effective method to assess the changes in coastal areas [19]. A lot of researches has been conducted using remote sensing to monitor the coastline changes affected by erosion and accretion [1, Fig. 1. Location of the study area. 7, 20, 21]. Their results yielded success in the areas where they monitored coastlines using remote sensing. According to Hussain et al. [16], remote sensing provides the best solution 2.2 Sampling method as far as acquiring and processing data in a quick and faster way compared to other methods like photogrammetry. The objective of this paper is to determine the original coastline of The samples were gathered at two distinct spots: the first one at Semarang city at coordinate Semarang city and Demak district in order to determine the original border area. This 6°56'57.24"S - 110°26'34.15"E and the second one at the end of coastline in Demak district information is essential to the government because knowing the coastline change will enable at coordinate 6°42'26.33"S - 110°36'55.84"E. The samples were taken from these coordinates them to understand the changing of physical environment better, thus developing coastal to identify the differences in the coastline from 1984 to 2018. policies that will promote sustainable management practices. Besides, knowing the location and coastline changes is helpful to engineers while designing coastal protection structures. To achieve the objective of the current study, we used remote sensing approach using Landsat 2.3 Analytical technique satellite imagery data dating from 1984 to 2018. Data used in this research was obtained from Landsat satellite imagery recorded on 23 February 1984, 28 June 1992, 23 April 2000, 24 April 2008, and a drone captured on 5 June 2018. The overlay from Landsat images and RGB overhead images were corrected through radiometric and geometric check (see Figure 2) and also integrated with the administrative limit from Morosari district in Demak regency and the border of Semarang City (Tanjungmas district). The processing method of the Landsat images consists of several steps as listed below; • The initial stage of data development is conducted through correlation with the standard procedure issued by the data supplier.

3 E3S Web of Conferences 202, 04001 (2020) https://doi.org/10.1051/e3sconf/202020204001 ICENIS 2020

• Early development began by synchronizing the spatial resolution between Landsat-4 MSS and Landsat series TM.ETM+. This process is done by resampling the Landsat MSS data so that the spatial resolution turns into 30 meters through cubic B-Spline algorithm and then followed by a geometric check on this Landsat data; thus the north and south value of the imaging is identical to the one on map (administration map) and can be overlaid on each other. • The development of RGB (Red Green Blue) composite for each imaging acquisition. • Digitizing the fourth RGB image as an analogy by digitalizing it on screen. • Analysis and calculation are done through the process of digitizing of each year’s image so that it is possible to find out the changes both from accretion and abrasion.

Fig. 2. Landsat Imaging (1984-2018).

3 Results and Discussions

3.1 Landsat Imaging Coastline Analysis The correction results are done through manual works (on-screen digitation). The results of the coastline changes can be seen in Figure 3. According to the coastline length calculation resulted from digitation, the lengths of the border area between Semarang and Demak

4 E3S Web of Conferences 202, 04001 (2020) https://doi.org/10.1051/e3sconf/202020204001 ICENIS 2020

• Early development began by synchronizing the spatial resolution between Landsat-4 (Tanjung Mas, Terboyo Wetan, Terboyo Kulon, and Trimulyo village, and villages along the MSS and Landsat series TM.ETM+. This process is done by resampling the Landsat coast in Demak district) were 48.2 km, 57.9 km, 63.7 km, 73.5 km, and 92.6 km in 1984, MSS data so that the spatial resolution turns into 30 meters through cubic B-Spline 1992, 2000, 2008, and 2018, respectively. As a result, the rate of coastline changes in 1984 algorithm and then followed by a geometric check on this Landsat data; thus the north to 1992, 1992 to 2000, 2000 to 2008, and 2008 to 2018 is 1.2 km per year, 0.7 km per year, and south value of the imaging is identical to the one on map (administration map) and 1.2 km per year, and 2.1 km per year, respectively (see Table 1). It was noticed that a big can be overlaid on each other. incline (28.9 km) occurred in 2000 - 2018. • The development of RGB (Red Green Blue) composite for each imaging acquisition. • Digitizing the fourth RGB image as an analogy by digitalizing it on screen. • Analysis and calculation are done through the process of digitizing of each year’s image so that it is possible to find out the changes both from accretion and abrasion.

Fig. 3. Coastline information, digitation result from Landsat imaging from 1984 to 2018. Table 1. Coastline length and the rate of coastline changes in 1984 to 2018.

Year Coastline length (Km) Rate of coastline changes(Km/year) 1984 48.2 - 1992 57.9 1.2 2000 63.7 0.7 2008 73.5 1.2 2018 92.6 2.1

The occurrence of abrasion and accretion in an area was identified by integrating two results of the digitized coastline imaging from different years. The results of integration and digitization from consecutive years are seen in Figure 3. Fig. 2. Landsat Imaging (1984-2018). Abrasion and accretion occured in Demak district and Semarang city with different sizes in different years, as shown in Table 2. In the period between 1984 to 1992, the largest 3 Results and Discussions abrasion occurred in Sayung sub-district, and a small part of it occurred in Karangtengah, Wedung, Bonang, Genuk, and north Semarang district as seen in Figure 4 (a) within a total area of 846.5 ha. The reason for this big abrasion occurrence in Sayung could be its location 3.1 Landsat Imaging Coastline Analysis that is open without any human or natural protection, thus easily eroded in case the waves advance to the land [3]. Figure 4 (b) depicts the accretion-affected areas, where the most The correction results are done through manual works (on-screen digitation). The results of dominant accretion occurred in Karangtengah, Sayung, Wedung, and north Semarang with a the coastline changes can be seen in Figure 3. According to the coastline length calculation total area of 496.5 ha (see Table 2). According to Xu [10], other than natural factors, man- resulted from digitation, the lengths of the border area between Semarang and Demak made factors also play a major role in coastline changes, for example coastline protection

5 E3S Web of Conferences 202, 04001 (2020) https://doi.org/10.1051/e3sconf/202020204001 ICENIS 2020

structures and beach nourishment projects that are built to reduce coastal erosion may lead to artificial accretion of the coastline. The overall abrasion that occurred in 1984 – 2018 in Semarang city and Demak district is higher than the accretion values. The abrasion affected people highly. The area that was previously used as settlement was covered by water. Table 2. Abrasion and Accretion area at the border of Demak district and Semarang city.

Fig. 4. Satellite imaging result of coast abrasion and accretion information. Figure 5 (a) shows the location where the abrasion process happened in 1992 to 2000. Abrasion took place in Karangtengah, Sayung, Wedung, Bonang, Genuk, and north Semarang with a total area of 349.9 ha. Figure 5 (b) shows the regions where accretion occurred in Karangtengah, Sayung, Bonang, Wedung, and north Semarang with a total area of 241.7 ha.

6 E3S Web of Conferences 202, 04001 (2020) https://doi.org/10.1051/e3sconf/202020204001 ICENIS 2020 structures and beach nourishment projects that are built to reduce coastal erosion may lead to artificial accretion of the coastline. The overall abrasion that occurred in 1984 – 2018 in Semarang city and Demak district is higher than the accretion values. The abrasion affected people highly. The area that was previously used as settlement was covered by water. Table 2. Abrasion and Accretion area at the border of Demak district and Semarang city.

1984-1992 1992-2000 2000-2008 2008-2018 Abr N City/ asio Accr Abra Accr Abra Accr Abra Accr Sub-district o District n etion sion etion sion etion sion etion (Ha) Ha) (Ha) Ha) (Ha) (Ha) (Ha) (Ha) Karangtenga 1 Demak 77.9 11.4 35.5 22.1 8.2 40.1 11.2 8.7 h 258. 128. 390. 460. 2 Sayung Demak 5.9 0.1 9.2 0.4 Fig. 5. Landsat satellite imaging result of coast abrasion and accretion 1992 – 2000. 8 8 2 6 131. 152. 311. 190. Abrasion process that occurred in the period 2000-2008 was relatively the same as the 3 Wedung Demak 1.6 80.1 86.2 82.7 1 4 9 8 previous period (1992-2000) but a bit different, as seen in Figure 6 (a). Abrasion occurred in 175. 103. Karangtengah, Sayung, Wajah, Bonang, Genuk, and north Semarang, with a total area of 4 Bonang Demak 0 18.9 39.7 3.7 69.4 0.1 3 4 682.7 ha. Figure 6 (b) shows the areas where accretion occurs, where the most dominant Semaran 141. 152. occurred in the sub-districts of Wedung, Karangtengah, Sayung, Bonang, Genuk, and north 5 Genuk 80.5 0 45.9 0 17.1 21.7 g 1 9 Semarang with an overall area of 493.2 ha (see Table 2). North Semaran 122. 477. 6 40.7 27.4 53.3 11.5 25.1 13.2 Semarang g 9 6 846. 496. 349. 241. 682. 493. 126. Total 910 5 5 9 7 7 2 8

Fig. 6. Landsat satellite imaging result of coast abrasion and accretion 2000 – 2008.

Figure 7 (a) shows the location where the abrasion process occurred between 2008 and 2018. It was seen that the magnitude of abrasion in this period is relatively larger than the previous periods. Abrasion occurred in Karangtengah, Sayung, Wajah, Bonang, Genuk, and north Semarang district with a total area of 910 ha. Whereas Figure 7 (b) shows the areas where accretion occurred, the most dominant accretion occurred in the sub-districts of Fig. 4. Satellite imaging result of coast abrasion and accretion information. Wedung, Karangtengah, Sayung, Bonang, Genuk, and north Semarang with a total area of Figure 5 (a) shows the location where the abrasion process happened in 1992 to 2000. 126.8 ha (see Table 2). Abrasion took place in Karangtengah, Sayung, Wedung, Bonang, Genuk, and north Semarang with a total area of 349.9 ha. Figure 5 (b) shows the regions where accretion occurred in Karangtengah, Sayung, Bonang, Wedung, and north Semarang with a total area of 241.7 ha.

7 E3S Web of Conferences 202, 04001 (2020) https://doi.org/10.1051/e3sconf/202020204001 ICENIS 2020

Fig. 7. Landsat satellite imaging result of coast abrasion and accretion 2008 – 2018. Figure 8 (a) shows that the coastline accretion process that seems to have been caused by sedimentation processes carried out by river flow, sea, and human activities like fishing, and the permanent settlement around the coast, while the coastline abrasion process as seen in Figure 8 (b) is due to the erosion caused by sea waves, sea level rise, currents, and other factors. By observing Figure 8, it can be deduced that the coastline's accretion is caused by sedimentation and land use changing into dikes, whereas the abrasion process can be assumed to be caused by the erosion due to sea waves.

Fig. 8. Accretion and Abrasion (1984 – 2018). When the results of the digitization of the coastline are fully integrated (see Figure 9), it is seen that the change in the coastline only occurs in the cape. If there is river flow in the cape, then the accretion process occurs in agricultural areas/ponds such as in Sayung sub- district, Bonang, Karangtengah. On the other hand, if there is no river flow, there will be an abrasion process, for example, the one that happened in Wedung sub-district while the accretion process occurred in Tanjung Mas Port, and north Semarang.

8 E3S Web of Conferences 202, 04001 (2020) https://doi.org/10.1051/e3sconf/202020204001 ICENIS 2020

3.2 Comparative analysis of Regional Spatial Plan and Landsat Imaging data

Based on the Spatial Planning of Semarang City and Demak Regency (2011-2031), the length of the coastline in the study location (Semarang city and Demak regency) is 62.8 km. While the length of the coastline based on remote sensing analysis using Landsat satellite image is 92.6 km. This shows that there is a difference of 29.8 km (see Figure 10). Analysis by Landsat satellite image gives a larger length compared to the one provided by Regional Spatial Plan 2011-2031.

9 E3S Web of Conferences 202, 04001 (2020) https://doi.org/10.1051/e3sconf/202020204001 ICENIS 2020

Fig. 10. Coastline change in 2018 according to Regional Spatial Plan and analysis from Landsat Imaging.

4 Conclusion This study investigated the coastline changes that occurred at the boundary of Semarang city and Demak district. The data from the Regional Spatial plan was compared to the data analysis from Landsat imaging. It was found out that the coastline underwent abrasion (1,592.7 ha) and accretion (620 ha) from 2000 to 2018. The major factors that contributed to the abrasion are sea waves, currents, and storms, while the factors that contributed to accretion are human activities like fishing, and the permanent settlement around the coast. The coastline determination by Government Regional Spatial Plan deviates from the original coastline condition by 5-15%. The study recommends a sustainable coastal management and action plan to be prepared by the government. This will be essential in preserving the coastal regions.

10 E3S Web of Conferences 202, 04001 (2020) https://doi.org/10.1051/e3sconf/202020204001 ICENIS 2020

References 1. A. A. Alesheikh, A. Ghorbanali, and N. Nouri, Int. J. Environ. Sci. Technol., 4(1), 61– 66 (2007) 2. S. Kaliraj, N. Chandrasekar, and K. K. Ramachandran, J. Remote Sens. Sp. Sci., 20(2), 265–282 (2017) 3. I. A, S. Anggoro, T. R. Soeprobowati, M. Helmi, and A. S. Khair, Jurnal Pendidikan IPA Indonesia, 8(1), 1-11 (2019) 4. S. Maiti and A. K. Bhattacharya, Mar. Geol. 257(1), 11–23 (2009) 5. A. Guariglia, A. Buonamassa, A. Losurdo, and R. Saladino, Ann. Geophys., 49(1), 295– 304 (2006) 6. A. Sabuncu, H. Ozener, and B. Turgut, Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci.,41(July), 1169-1174 (2016) 7. U. Natesan, A. Parthasarathy, R. Vishnunath, G. E. J. Kumar, and A. Ferrer, Aquat. Procedia, 4, 325–332 (2015) 8. K. A. Addo, M. Walkden, and J. P. Mills, ISPRS J. Photogramm. Remote Sens., 63, 543–558 (2008) 9. E Tamassoki, H. Amiri, and Z. Soleymani, IOP Conf. Ser. Earth Environ. Sci., 20, 1–9 (2014) 10. N. Xu, Atmosphere (Basel)., 9(3),107 (2018) 11. K. Di, R. Ma, J. Wang, and R. Li, “Coastal Mapping and Change Detection Using High- Resolution IKONOS Satellite Imagery, In Proceedings of the 2003 annual national conference on Digital government research, (2003) 12. S. Kurt, A. Karaburun, and A. Demirci, Sci. Res. Essays, 5(19), 3009-3017 (2010) Fig. 10. Coastline change in 2018 according to Regional Spatial Plan and analysis from Landsat 13. H. Liu and K. C. Jezek, Int. J. Remote Sens., 25(5), 937–958 (2004) Imaging. 14. C. Liu, J. Yang, J. Yin, and W. An, IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens., 9(11), 4908–4920 (2016) 4 Conclusion 15. X. Li and M. C. Damen, J. Mar. Syst., 82, S54–S61 (2010) 16. M. A. Hussain, Y. Tajima, K. Gunasekara, S. Rana, and R. Hasan, IOP Conf. Ser. Earth This study investigated the coastline changes that occurred at the boundary of Semarang city Environ. Sci, 20, 1–10 (2014). and Demak district. The data from the Regional Spatial plan was compared to the data 17. A. M Shamsul and U. Kabir, Am. J. Geogr. Inf. Syst., 2(1), 15–18 (2013) analysis from Landsat imaging. It was found out that the coastline underwent abrasion 18. R. Gens, Int. J. Remote Sens., 31(7), 1819–1836 (2010) (1,592.7 ha) and accretion (620 ha) from 2000 to 2018. The major factors that contributed to 19. E. P. Green, P. J. Mumby, A. J. Edwards, and C. D. Clark, Coast. Manag., 24(1), 1–40 the abrasion are sea waves, currents, and storms, while the factors that contributed to (1996) accretion are human activities like fishing, and the permanent settlement around the coast. 20. P. Chand and P. Acharya, Int. J. Geomatics Geosci., 1(3), 436–455 (2010) The coastline determination by Government Regional Spatial Plan deviates from the original 21. A. Parthasarathy and U. Natesan, Nat. Hazards, 75(2), 1713–1729 (2015) coastline condition by 5-15%. The study recommends a sustainable coastal management and action plan to be prepared by the government. This will be essential in preserving the coastal regions.