IOP Conference Series: Earth and Environmental Science PAPER • OPEN ACCESS Numerical study of tides in the Bay of Sabang To cite this article: H A A Wahyudi et al 2019 IOP Conf. Ser.: Earth Environ. Sci. 348 012096 View the article online for updates and enhancements. This content was downloaded from IP address 116.206.33.45 on 01/12/2019 at 00:55 ICFAES 2019 IOP Publishing IOP Conf. Series: Earth and Environmental Science 348 (2019) 012096 doi:10.1088/1755-1315/348/1/012096 Numerical study of tides in the Bay of Sabang H A A Wahyudi1, Y Haditiar2, M Ikhwan2, R Wafdan1, I Setiawan1,2, M Muhammad1, S Sugianto3 and S Rizal1,2,* 1Department of Marine Sciences, Faculty of Marine and Fisheries, UniversitasSyiah Kuala, Banda Aceh, 23111, Indonesia. 2Graduate School of Mathematics and Applied Science, UniversitasSyiah Kuala, Banda Aceh, 23111, Indonesia. 3Department of Soil Sciences, Faculty of Agriculture, Universitas Syiah Kuala, Banda Aceh, 23111, Indonesia *Email: [email protected] Abstract. This study aims to observe tidal sea behavior in Bay of Sabang by using a two- dimensional hydrodynamic model. The research domain was obtained from SRTM15, which had been interpolated so that it has a spatial resolution Δx = Δy = 0.1 minutes (185 meters). As a generator force, the open boundaries model is given five main tidal components (M2, S2, K1, N2, O1) obtained from TPXO 7.2 data. The model is simulated for 31 days with time step based on CFL condition criteria. Our model verification is quite good when compared with TMD prediction data (r = 0.9996). It shows that our model can be used for tidal hydrodynamics studies in Bay of Sabang. The results show that tides in Bay of Sabang are predominantly mixed prevailing semi-diurnal (F = 0.857). The speed of tidal currents in Bay of Sabang is quite small both for spring and neap tides. The difference in velocity between spring and neap tides reaches 0.6 cm/s. 1. Introduction Bay of Sabang is administratively located on Weh Island, north of Aceh, Indonesia. The hydrodynamics of Sabang Waters are influenced by adjacent waters such as the Andaman Sea, the Malacca Strait, and the Indian Ocean [1]. The Andaman Sea is deep sea waters (an average of 1100 m) located in the northeast of the Indian Ocean, bordering Myanmar, Thailand, Malaysia and Indonesia [2]. Andaman Sea and Indian Ocean are strongly affected by the monsoonal system. It can produce different surface current circulation during NE monsoon and SW monsoon. Likewise, the productivity and temperature of the sea is quite extensive [3, 4]. Meanwhile, the Malacca Strait is very dominantly affected by tides [5]. Tidal phenomena are an event of periodic high and low sea levels that occur throughout the hemisphere due to the attraction of the earth, moon and sun [6]. Tides are the result of the force of gravitational attraction and centrifugal effect (thrust out of the center of rotation) [7].The process of circulating water masses in coastal waters is dominantly generated by tidal currents [8]. Tides in the Indian Ocean, the Bay of Bengal, and the Andaman Sea have an impact on tides in the Bay of Sabang and the Malacca Strait [9, 10]. Tides of the Indian Ocean cross Aceh waters and Andaman Sea. This phenomenon can be seen from the increase in tidal phase changes in shallow waters. Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. Published under licence by IOP Publishing Ltd 1 ICFAES 2019 IOP Publishing IOP Conf. Series: Earth and Environmental Science 348 (2019) 012096 doi:10.1088/1755-1315/348/1/012096 Based on research [1] with a resolution of 18.5 km, it is known that tides in Aceh waters are divided into semi-diurnal and semi-diurnal mixtures. Bathymetry and island shape produce height variations and tidal currents. Thus the resolution of this model is important in getting tidal models in Aceh waters [11, 12, 13]. According to [14], representation of good seawater conditions is needed in coastal management and sea security.Coastal dynamics, such as waves and tidal currents, greatly affects the abiotic content of the sea [15].The development of the marine model is currently developing for the prediction of seawater conditions in coastal areas such as tides and sea temperatures [16]. So far, the tidal information at Sabang Bay is still limited. Therefore, this paper aims toprepare and perform the two-dimensional numeric tide modelsimulation. 2. Materials and Methods Sabang Bay was previously simulated with a 1D model bySetiawan et al., [17]. They showed good results but did not cover most of the Sabang Bay domain. This research uses numerical model that was derived by [18]. This model has consisted of nonlinear equations and linear Navier-Stokes. Models with linear equations have been applied in western of Aceh [19]and nonlinear equations have been successfully applied in the Malacca Strait [12], Gulf of Thailand [20], Aceh Barat Daya Waters [21], and Aceh Besar Waters [22]. In this study, we use the numerical model based on the nonlinear Navier- Stokes equation [18].The resolution of the model is ∆x = ∆y = 0.1 minutes, while the time-step is derived according to the CFL condition. The domain bathymetry is obtained from SRTM 15 at coordinate 95.22 - 95.33 E and 5,835-5.92 N. Figure 1. The domain of model 2 ICFAES 2019 IOP Publishing IOP Conf. Series: Earth and Environmental Science 348 (2019) 012096 doi:10.1088/1755-1315/348/1/012096 The data is interpolated at Δx = Δy = 0.1 minutes so that the grid size is 67x52. In open boundary, it was given main tidal components (M2, S2, K1, N2 and O1) obtained from TPXO 7.2. In the bottom friction, the force was given uniformly with parameter r = 0.0025 and TVD advection scheme. The results of the model are verified with predictions of Tidal Model Driven (TMD) during January 2017. Descriptions of the components of sea tides can be seen in Table 1. Table 1. Tidal component and parameter Amplitude Phase Constant Period Amplitude Phase correction Correction M2 12.4206024 1.01184 126.95103 0.4372 83.65 S2 12.0000000 0.99932 0.13461 0.2060 116.5 K1 23.9344704 0.97829 1.63348 0.1046 114.9 N2 12.6583488 1.01183 47.76974 0.0833 74.8 O1 25.8193416 0.96194 129.8077 0.447 182.58 The tidal motion of the sea water is expressed by the following mathematical equation [17] () = cos( + − ) In the equation, () is the water level at time t, is the average water level, is the astronomical correction factor for the amplitude of the tidal generating element , is the angular velocity, is the correction factor for phase, and is phase. 3. Results and Discussion The tidal type in the Sabang bay waters is a semi-diurnal mixture, where the fromzal number obtained is 0.857. This is in accordance with the research conducted by Rizal [1] where the type of tides in the western waters of Sumatra mostly has semidiurnal or semidiurnal mixed tidal. Figure 2. The sea level of tides that derived from 2D model (blue lines) and TMD (red stars) (5.8767 oN and 95.2733 oE). The 2D model and TMD coincide and show strong correlation (with correlation coefficientr = 0.99). 3 ICFAES 2019 IOP Publishing IOP Conf. Series: Earth and Environmental Science 348 (2019) 012096 doi:10.1088/1755-1315/348/1/012096 The verification results show that the model data equal to TMD predictions (see Figure 2). At spring tide, the range is 1.524 m for the model and 1.5216 m for the TMD, while for neap tides it is 0.4508 m for the model and 0.4637 for TMD. Ebb current is indicated by the seaward flow,while flood current by the inland flow. The highest tidal current occurs during ebb current and flood current in spring tide conditions (Figure 3 and 4). Whereas in neap tide conditions at ebb current and flood current, there is the lowest tidal current (Figure 5 and 6). Spring tide is the highest flood and lowest ebb conditions, while the neap tide is the lowest flood and highest ebb conditions. When the highest flood and the lowest ebb, both conditions are not affected by tidal currents. Figure 3. Ebb current during spring tides (Jan. 13, 2017). The contour lines show the magnitudes of current while the vectors only show directions of current. Figure 4. Flood current during spring tides (Jan. 12, 2017). The contour lines show the magnitudes of current while the vectors only show directions of current. The flood towards the ebb moves out of the coast to the north, but in the northwest, the current turn slightly and move towards the southwest. The highest current velocity is 0.8 cm/s in the east of the Bay of Sabang, and the lowest current velocity is 0.2 cm/s in the west. Current velocity at low tide conditions leads to pairs of Figure 3 not much different even tends to be the same as the ebb tide, where the highest current velocity is 0.8 cm/s and the lowest is 0.2 cm/s. Significant differences occur when flood tide current moves into the shoreline from the north. In other words, the flood tide is 4 ICFAES 2019 IOP Publishing IOP Conf.