Tidal Current Power Resources and Influence of Sea-Level Rise in the Coastal Waters of Kinmen Island, Taiwan

The tidal current power (TCP) resource, the impact of TCP extraction on hydrodynamics and the influence of sea-level rise (SLR) on TCP output in the coastal waters of Kinmen Island (Taiwan) are investigated using a state-of-the-art unstructured-grid depth-integrated numerical model. The model was driven by eight tidal constituents extracted from a global tidal prediction model and verified with time series of measured data for tide level and depth-averaged current. The simulations showed reasonable agreement with the observations; the skill index was in the excellent (0.71–0.93) range with regard to simulating tide level and currents. Model predictions indicated that the channel between Kinmen and Lieyu serves as an appropriate site for deploying the tidal turbines because of its higher tidal current and deeper water depth. The bottom friction approach was utilized to compute the average TCP over a spring-neap cycle (i.e., 15 days). Mean TCP reached its maximum to 45.51 kW for a coverage area of 0.036 km2 when an additional turbine friction coefficient (Ct) increased to 0.08, and a cut-in speed of 0.5 m/s was used. The annual TCP output was estimated to be 1.08 MW. The impact of TCP extraction on the change in current is significant, with a maximum reduction rate of instant current exceeding 60%, and the extent of influence for the average current is 1.26 km in length and 0.30 km in width for the −0.05 m/s contour line. However, the impact of TCP extraction on the change in tide level is insignificant; the maximum change in amplitude is only 0.73 cm for the K2 tide. The influence of SLR on the TCP output in Kinmen waters was also estimated. Modeling assessments showed that due to SLR produces faster tidal current, the annual TCP output increased to 1.52 MW, 2.01 MW, 2.48 MW and 2.97 MW under the same cut-in speed and coverage area conditions when SLR 0.25 m, SLR 0.5 m, SLR 0.75 m and SLR 1.0 m were imposed on the model.


Introduction
Global warming is expected to become more intense during this century due to the increasing emissions of greenhouse gases [1]. The relationship between the burning of fossil fuels and climate change is becoming more established. Therefore, renewable marine power resources, in the form of tides, waves, tidal currents and sea temperature (thermal energy), are potential candidates for possible energy sources [2,3]. The extraction of kinetic energy from tidal currents is a desirable form since it is predictable, regular, and has a high energy density with low greenhouse gas emissions [4].
Recently, many researchers have made great efforts to identify the most appropriate locations for the placement of tidal turbines and to assess how much tidal current power (TCP) can be extracted from a specific site [5][6][7]. Tidal currents are highly dependent on the geographic characteristics. For instance, a channel with varying cross-sections linking with two large water bodies is usually rich in TCP [5][6][7][8][9]. In addition, tidal currents can be accelerated in the area around a headland tip [10][11][12][13][14][15].
Although tidal turbine operation has less negative influences on the marine environment, the hydrodynamics of a tidal system, such as tide level and current speed, might change due to the extraction of TCP in a certain sea or estuarine area [2,16]. Thus, an assessment of the environmental impact should be conducted for a turbine farm before the deployment of tidal turbines. Although global warming may be far worse than thought and sea-level rise (SLR) is expected to accelerate over the next century [17,18], few studies have concentrated on estimating the effect of SLR on TCP output [19,20]. A numerical hydrodynamic model is a powerful tool to assess the distribution of TCP resources; a two-dimensional model (2D model) is required for shallow coastal waters [7,21,22], and a three-dimensional model (3D model) is required for the deep sea [2, 23,24]. In addition, two approaches are commonly incorporated into numerical models to simulate the effect of TCP extraction on hydrodynamics. The first method is the bottom friction approach, and a 2D model is generally employed; the presence of tidal turbines in the water is represented by adding an additional bottom friction in the model [25,26]. The second method is the momentum-sink approach, which usually adopts a 3D model; in this model, the momentum is diminished by adding the TCP extraction effect terms in the momentum equations [2, 23,24]. However, these two approaches produce similar results [24].
Kinmen Island (Kinmen Taiwan), located just off the southeastern coast of mainland China (Figure 1a), covers an area of 153 km 2 and has a population of 135,223 (2017). The Tashan Power Plant (the red solid circle in Figure 1b) on Kinmen Island is a fuel-fired power plant built in 1967. The total power generation capacity of the Tashan Power Plant is approximately 34 MW. It is currently inadequate for meeting the power resource demand for the residents of Kinmen Island. Therefore, the TCP extracted from the coastal waters of Kinmen Island has been considered as an alternative power generation method because of its lower pollution.
for the placement of tidal turbines and to assess how much tidal current power (TCP) can be extracted from a specific site [5][6][7]. Tidal currents are highly dependent on the geographic characteristics. For instance, a channel with varying cross-sections linking with two large water bodies is usually rich in TCP [5][6][7][8][9]. In addition, tidal currents can be accelerated in the area around a headland tip [10][11][12][13][14][15].
Although tidal turbine operation has less negative influences on the marine environment, the hydrodynamics of a tidal system, such as tide level and current speed, might change due to the extraction of TCP in a certain sea or estuarine area [2,16]. Thus, an assessment of the environmental impact should be conducted for a turbine farm before the deployment of tidal turbines. Although global warming may be far worse than thought and sea-level rise (SLR) is expected to accelerate over the next century [17,18], few studies have concentrated on estimating the effect of SLR on TCP output [19,20]. A numerical hydrodynamic model is a powerful tool to assess the distribution of TCP resources; a two-dimensional model (2D model) is required for shallow coastal waters [7,21,22], and a three-dimensional model (3D model) is required for the deep sea [2, 23,24]. In addition, two approaches are commonly incorporated into numerical models to simulate the effect of TCP extraction on hydrodynamics. The first method is the bottom friction approach, and a 2D model is generally employed; the presence of tidal turbines in the water is represented by adding an additional bottom friction in the model [25,26]. The second method is the momentum-sink approach, which usually adopts a 3D model; in this model, the momentum is diminished by adding the TCP extraction effect terms in the momentum equations [2, 23,24]. However, these two approaches produce similar results [24].
Kinmen Island (Kinmen Taiwan), located just off the southeastern coast of mainland China (Figure 1a), covers an area of 153 km 2 and has a population of 135,223 (2017). The Tashan Power Plant (the red solid circle in Figure 1b) on Kinmen Island is a fuel-fired power plant built in 1967. The total power generation capacity of the Tashan Power Plant is approximately 34 MW. It is currently inadequate for meeting the power resource demand for the residents of Kinmen Island. Therefore, the TCP extracted from the coastal waters of Kinmen Island has been considered as an alternative power generation method because of its lower pollution.
The objectives of this study are as follows: (1) investigate the distribution of TCP resources and find an appropriate site for deploying the tidal turbines in the coastal waters of Kinmen Island; (2) calculate maximum mean TCP output over a spring-neap cycle using a high-resolution two-dimensional hydrodynamic model for the turbine site; (3) evaluate the impact of TCP extraction on the tide level and tidal current; and (4) assess the influence of SLR on TCP output.  The objectives of this study are as follows: (1) investigate the distribution of TCP resources and find an appropriate site for deploying the tidal turbines in the coastal waters of Kinmen Island; (2) calculate maximum mean TCP output over a spring-neap cycle using a high-resolution two-dimensional hydrodynamic model for the turbine site; (3) evaluate the impact of TCP extraction on the tide level and tidal current; and (4) assess the influence of SLR on TCP output.

Hydrodynamic Model
A state-of-the-art unstructured-grid hydrodynamic model called Semi-implicit Cross-scale Hydroscience Integrated System Model (SCHISM), has been used to simulate tide levels and tidal currents in the coastal waters of Kinmen Island. SCHISM is a modeling system that is a derivative product of the original SELFE model [27], with many enhancements and upgrades, including a new extension to the large-scale eddying regime and a seamless cross-scale capability from the creek to ocean scale [28].
Similar to the SELFE model, SCHISM is based on an unstructured grid and is designed for seamless simulation of three-dimensional baroclinic circulation across creek-lakeriver-estuary-shelf-ocean scales. A highly efficient, parallel computing, robust and accurate semi-implicit finite-element/finite-volume method with an Eulerian-Lagrangian algorithm is implemented in the SCHISM model to solve Reynolds-averaged Navier-Stokes equations for transport of heat, salt and tracers. The wetting and drying scheme is naturally incorporated into the model for inundation studies. Due to their highly flexible frameworks, the SELFE and SCHISM models have been widely used to solve cross-scale problems: general circulation [29], storm surges [30], tsunami hazards [31], water quality [32], oil spills [33], sediment transport [34,35], biogeochemistry [36,37], and assessment of tidal stream energy resources [9,38]. A two-dimensional version of SCHISM, SCHISM-2D, is employed in this study since Kinmen Island is surrounded by shallow coastal waters (as shown in Figure 2a). The governing equations of the SCHISM-2D model with hydrostatic form and Boussinesq approximation in the Cartesian coordinate system are given as follows: where η(x, y, t) is the free-surface elevation; u(x, y, t) and v(x, y, t) are the depth-averaged velocity in the x, y direction, respectively; f is the Coriolis factor; g is the acceleration due to gravity;ψ is the earth's tidal potential; α is the effective earth elasticity factor; ρ is the density of water; P A (x, y, t) is the atmospheric pressure at the free surface; τ sx and τ sy refer to the wind shear stress in the x, y direction, respectively; τ bx and τ by are the bottom shear stress in the x, y direction, respectively; and H = η + h, where h(x, y) is the bathymetric depth. The formulas of the bottom shear stress are as follows: where C d is the bottom drag coefficient, which is determined through model validation. C t is an additional bottom friction coefficient introduced by tidal turbine. C t is equal to zero when turbine is absent; while is greater than zero when turbine is present.

Calculation of Tidal Current Power
The extractable TCP over a tidal cycle (spring-neap cycle) for a turbine can usually be estimated as follows: where P is the tidal power extraction rate in watts; M is the number of turbines corresponding to the maximum average TCP output; e C is the turbine thrust coefficient; and e A is the swept area of the turbine. The over-bar above velocity term indicates the time average over a spring-neap cycle. Although Equation (5) has widely been used, the effects of additional bottom friction induced by turbine deployment on hydrodynamics cannot be appropriately represented. Garrett and Cummins [5] and Bryden and Couch [39] proposed that an additional drag term on the flow will arise from the deployment of turbines. The tidal current power output ( t P ) can be evaluated considering the turbine-induced additional bottom friction coefficient ( t C ) using the bottom friction approach [9,19,21,25,40]: where t A is the area covered by the tidal turbine field; u and v are obtained from SCHISM-2D.

Model Configuration
The computational domain covers the entire Kinmen Islands group, including Kinmen and Lieyu, which are located just off the southeastern coast of mainland China and are surrounded by the Taiwan Strait (Figure 1a). A global gridded bathymetry dataset with a resolution of 30 arc-second obtained from the General Bathymetric Chart of the Oceans (GEBCO) was used in this study. It was generated by combining quality-controlled ship depth soundings with interpolation between sounding points guided by satellite-derived gravity data. The modeling domain is composed of 197,020 triangle mesh grid cells and 100,714 nodes. Coarse mesh with 1000 m resolution was arranged to open ocean boundaries in the Taiwan Strait, while fine mesh with 50 m resolution was used around Kinmen Island (Figure 2b). Once the unstructured grid was constructed, the gridded bathymetry data were interpolated to each node (Figure 2a). The SCHISM-2D model was driven by harmonic constants of eight major tidal constituents (M2, S2, N2, K2, K1, O1, P1, and Q1), which were extracted from the global inverse tidal model, TPXO [41,42]. A time step of 60 s was set to ensure the stability of the model. The initial conditions for the hydrodynamic model are the null free surface elevation and null velocity. Meteorological and

Calculation of Tidal Current Power
The extractable TCP over a tidal cycle (spring-neap cycle) for a turbine can usually be estimated as follows: where P is the tidal power extraction rate in watts; M is the number of turbines corresponding to the maximum average TCP output; C e is the turbine thrust coefficient; and A e is the swept area of the turbine. The over-bar above velocity term indicates the time average over a spring-neap cycle. Although Equation (5) has widely been used, the effects of additional bottom friction induced by turbine deployment on hydrodynamics cannot be appropriately represented. Garrett and Cummins [5] and Bryden and Couch [39] proposed that an additional drag term on the flow will arise from the deployment of turbines. The tidal current power output (P t ) can be evaluated considering the turbine-induced additional bottom friction coefficient (C t ) using the bottom friction approach [9,19,21,25,40]: where A t is the area covered by the tidal turbine field; u and v are obtained from SCHISM-2D.

Model Configuration
The computational domain covers the entire Kinmen Islands group, including Kinmen and Lieyu, which are located just off the southeastern coast of mainland China and are surrounded by the Taiwan Strait ( Figure 1a). A global gridded bathymetry dataset with a resolution of 30 arc-second obtained from the General Bathymetric Chart of the Oceans (GEBCO) was used in this study. It was generated by combining quality-controlled ship depth soundings with interpolation between sounding points guided by satellite-derived gravity data. The modeling domain is composed of 197,020 triangle mesh grid cells and 100,714 nodes. Coarse mesh with 1000 m resolution was arranged to open ocean boundaries in the Taiwan Strait, while fine mesh with 50 m resolution was used around Kinmen Island (Figure 2b). Once the unstructured grid was constructed, the gridded bathymetry data were interpolated to each node ( Figure 2a). The SCHISM-2D model was driven by harmonic constants of eight major tidal constituents (M 2 , S 2 , N 2 , K 2 , K 1 , O 1 , P 1 , and Q 1 ), which were extracted from the global inverse tidal model, TPXO [41,42]. A time step of 60 s was set to ensure the stability of the model. The initial conditions for the hydrodynamic model are the null free surface elevation and null velocity. Meteorological and hydrological conditions are not considered in the current study. The simulations were only run in barotropic mode, i.e., density gradients are excluded in the model.

Model Performance Indicators
Two criteria were used to evaluate the model performance for the simulation of tide level and the depth-averaged tidal current, which are the mean absolute percentage error (MAPE) and the skill [43]. It should be noted that the MAPE can only be used when there are no zero value among observations. The optimal value of the skill is 1, indicating perfect performance. A skill value ranging between 0.65 and 1.0 indicates excellent performance. Very good performance is in the range of 0.5-0.65; good performance is in the range of 0.2-0.5; and a skill value less than 0.2 expresses poor performance [44]: where N is the total number of data points; O s i is the simulation value; O o i is the measured data; and O o is the mean value of the measured data.

Model Validation for Tide Level and Tidal Currents
The observational data provided by the Taiwan Central Weather Bureau were used to validate the hydrodynamic model for predicting tide level and tidal currents (the gauging station locations are shown in Figure 1b). The date was measured during the winter (from December 2016 to January 2017); this means that the observations are not affected by severe weather conditions (e.g., typhoon). Although the northeast monsoon is sometimes stronger in this area in winter, the effects of winter monsoon winds on tide level and depth-averaged tidal current are not significant. Because the Kinmen buoy is placed at a location with the water depth of 20-25 m. In addition, the unusual values such as −9999.0 were also filtered. The simulation period, a total of 29 days, was from 12 December 2016 to 9 January 2017. The first 14 days are designed to confirm that the model results will not be affected by the initial conditions. Figure 3 shows the model-data comparison for tide level at Suetao (Figure 3a) and Liaolowan (Figure 3b). Excellent agreement between the simulated and measured time series of tide level is observed at these two tidal stations. The model performance in the tidal simulations is illustrated in Table 1. The MAPE and skill values are 17.68% and 0.93 for Suetao and 15.41% and 0.93 for Liaolowan, respectively. The predictions of the depth-averaged tidal currents for the Kinmen buoy were compared with the measurements; the comparison results are presented in Figure 4. The results reveal that the hydrodynamic model faithfully reproduces the peaks and phases of the tidal current. Table 1 indicates that the MAPE is 20.57%, and the skill score of 0.71 indicates an excellent performance of the tidal current monsoon winds on tide level and depth-averaged tidal current are not significant. Because the Kinmen buoy is placed at a location with the water depth of 20-25 m. In addition, the unusual values such as −9999.0 were also filtered. The simulation period, a total of 29 days, was from 12 December 2016 to 9 January 2017. The first 14 days are designed to confirm that the model results will not be affected by the initial conditions. Figure 3 shows the model-data comparison for tide level at Suetao ( Figure 3a) and Liaolowan (Figure 3b). Excellent agreement between the simulated and measured time series of tide level is observed at these two tidal stations. The model performance in the tidal simulations is illustrated in  The predictions of the depth-averaged tidal currents for the Kinmen buoy were compared with the measurements; the comparison results are presented in Figure 4. The results reveal that the hydrodynamic model faithfully reproduces the peaks and phases of the tidal current. Table 1 indicates that the MAPE is 20.57%, and the skill score of 0.71 indicates an excellent performance of the tidal current simulations. The bottom drag coefficient ( d C ) was determined as 0.002 through model validation. The SCHISM-2D model can therefore be practically applied to the TCP calculations.  The distribution of the average tidal current over a spring-neap cycle (i.e., 15 days) around Kinmen Island is shown in Figure 5a. A higher current is observed in a narrow channel between Kinmen and Lieyu, with a maximum current exceeding 2.0 m/s. This suggests that the area is appropriate for the deployment of tidal power turbines. Water depths in the range of 20-50 m are also a requirement, in addition to high tidal velocity [4]. Therefore, Figure 5b shows a desirable site  The predictions of the depth-averaged tidal currents for the Kinmen buoy were compared with the measurements; the comparison results are presented in Figure 4. The results reveal that the hydrodynamic model faithfully reproduces the peaks and phases of the tidal current. Table 1 indicates that the MAPE is 20.57%, and the skill score of 0.71 indicates an excellent performance of the tidal current simulations. The bottom drag coefficient ( d C ) was determined as 0.002 through model validation. The SCHISM-2D model can therefore be practically applied to the TCP calculations.  The distribution of the average tidal current over a spring-neap cycle (i.e., 15 days) around Kinmen Island is shown in Figure 5a. A higher current is observed in a narrow channel between Kinmen and Lieyu, with a maximum current exceeding 2.0 m/s. This suggests that the area is appropriate for the deployment of tidal power turbines. Water depths in the range of 20-50 m are also a requirement, in addition to high tidal velocity [4]. Therefore, Figure 5b   The distribution of the average tidal current over a spring-neap cycle (i.e., 15 days) around Kinmen Island is shown in Figure 5a. A higher current is observed in a narrow channel between Kinmen and Lieyu, with a maximum current exceeding 2.0 m/s. This suggests that the area is appropriate for the deployment of tidal power turbines. Water depths in the range of 20-50 m are also a requirement, in addition to high tidal velocity [4]. Therefore, Figure 5b shows a desirable site (gray box), with an area of 35762.44 m 2 (0.036 km 2 ), depths greater than 20 m and average velocity of 0.43 m/s, for deploying the turbines.

Assessment of Maximum Tidal Current Power
According to a typical power curve, Robins et al. [45] assumed a cut-in velocity of 0.7 m/s and rated velocity of 2.7 m/s in the evaluation of tidal current power. However, this was an idealized representation of the potential performance of a turbine [46]. To assess the maximum potential TCP output in this study, the cut-in speed was therefore set to 0.5 m/s [47] and the rated velocity was not specified because the average tidal currents did not exceed 2.0 m/s at the turbine site.
The depth-averaged tidal currents were calculated by SCHISM-2D. TCP outputs were computed using Equations (6) and (7); the area covered by the tidal turbine field ( t A ) was 35762.44 m 2 (0.036 km 2 ), and the turbine drag term ( t C ) incrementally increased from 0.04 to 0.24. t C was added to d C in Equations (4a) and (4b), this means that the bottom shear stress would be increased with increasing t C , and consequently slow the horizontal velocity in the hydrodynamic model.
Each simulation was carried out for each value of t C . In this way, the average TCP output over a spring-neap cycle increased with increasing t C , resulting in a reduction of the average velocity.
There is a limit to the available TCP output as too many turbines not only create drag but also slow the current and therefore reduce the TCP output [4,17]. The maximum average TCP output is estimated to be 45.51 kW when 0.08 t C = , cut-in speed is 0.5 m/s, turbine coverage area is 35762.44 m 2 and the average velocity is 0.33 m/s (as shown in Figure 6). This indicates that the annual TCP output was estimated to be 1.11 MW ( 45.51 kW 365 /15 /1000 1.11 MW × = ) at the turbine deployment site. The results shown in Figure 6 are similar to those of other studies for the calculations of TCP output using numerical models or theoretical methods [5,9,19,21].

Assessment of Maximum Tidal Current Power
According to a typical power curve, Robins et al. [45] assumed a cut-in velocity of 0.7 m/s and rated velocity of 2.7 m/s in the evaluation of tidal current power. However, this was an idealized representation of the potential performance of a turbine [46]. To assess the maximum potential TCP output in this study, the cut-in speed was therefore set to 0.5 m/s [47] and the rated velocity was not specified because the average tidal currents did not exceed 2.0 m/s at the turbine site.
The depth-averaged tidal currents were calculated by SCHISM-2D. TCP outputs were computed using Equations (6) and (7); the area covered by the tidal turbine field (A t ) was 35762.44 m 2 (0.036 km 2 ), and the turbine drag term (C t ) incrementally increased from 0.04 to 0.24. C t was added to C d in Equations (4a) and (4b), this means that the bottom shear stress would be increased with increasing C t , and consequently slow the horizontal velocity in the hydrodynamic model. Each simulation was carried out for each value of C t . In this way, the average TCP output over a spring-neap cycle increased with increasing C t , resulting in a reduction of the average velocity. There is a limit to the available TCP output as too many turbines not only create drag but also slow the current and therefore reduce the TCP output [4,17]. The maximum average TCP output is estimated to be 45.51 kW when C t = 0.08, cut-in speed is 0.5 m/s, turbine coverage area is 35762.44 m 2 and the average velocity is 0.33 m/s (as shown in Figure 6). This indicates that the annual TCP output was estimated to be 1.11 MW (45.51 kW × 365/15/1000 = 1.11 MW) at the turbine deployment site. The results shown in Figure 6 are similar to those of other studies for the calculations of TCP output using numerical models or theoretical methods [5,9,19,21].

Assessment of Maximum Tidal Current Power
According to a typical power curve, Robins et al. [45] assumed a cut-in velocity of 0.7 m/s and rated velocity of 2.7 m/s in the evaluation of tidal current power. However, this was an idealized representation of the potential performance of a turbine [46]. To assess the maximum potential TCP output in this study, the cut-in speed was therefore set to 0.5 m/s [47] and the rated velocity was not specified because the average tidal currents did not exceed 2.0 m/s at the turbine site.
The depth-averaged tidal currents were calculated by SCHISM-2D. TCP outputs were computed using Equations (6) and (7); the area covered by the tidal turbine field ( t A ) was 35762.44 m 2 (0.036 km 2 ), and the turbine drag term ( t C ) incrementally increased from 0.04 to 0.24. t C was added to d C in Equations (4a) and (4b), this means that the bottom shear stress would be increased with increasing t C , and consequently slow the horizontal velocity in the hydrodynamic model.
Each simulation was carried out for each value of t C . In this way, the average TCP output over a spring-neap cycle increased with increasing t C , resulting in a reduction of the average velocity.
There is a limit to the available TCP output as too many turbines not only create drag but also slow the current and therefore reduce the TCP output [4,17]. The maximum average TCP output is estimated to be 45.51 kW when 0.08 t C = , cut-in speed is 0.5 m/s, turbine coverage area is 35762.44 m 2 and the average velocity is 0.33 m/s (as shown in Figure 6). This indicates that the annual TCP output was estimated to be 1.11 MW ( 45.51 kW 365 /15 /1000 1.11 MW × = ) at the turbine deployment site. The results shown in Figure 6 are similar to those of other studies for the calculations of TCP output using numerical models or theoretical methods [5,9,19,21].

Influence of Tidal Current Power Extraction on Hydrodynamics
TCP extraction may alter the hydrodynamics (e.g., tide level and currents) over an area where the turbines are deployed. To understand the impacts of turbine operation on tidal dynamics, time series and spatial differences of tide level and currents were compared for an area with and without turbines. Due to C t = 0.08 produced the maximum mean TCP, Figure 7 depicts a tidal water level time series for C t = 0 and C t = 0.08 at the turbine site. The turbine-induced tide level change is minor; the difference is difficult to distinguish in Figure 7. The spatial distribution of tide level differences is shown in Figure 8a. The tide level difference at each node was obtained from the average tide level over a spring-neap cycle for C t = 0.08 minus that for C t = 0. The differences in the average tide level between C t = 0.08 and C t = 0 in the spatial distribution are negligible, although an enlargement of the turbine site is shown in Figure 8b.

Influence of Tidal Current Power Extraction on Hydrodynamics
TCP extraction may alter the hydrodynamics (e.g., tide level and currents) over an area where the turbines are deployed. To understand the impacts of turbine operation on tidal dynamics, time series and spatial differences of tide level and currents were compared for an area with and without turbines. Due to 0.08 t C = produced the maximum mean TCP, Figure 7 depicts a tidal water level time series for 0 t C = and 0.08 t C = at the turbine site. The turbine-induced tide level change is minor; the difference is difficult to distinguish in Figure 7. The spatial distribution of tide level differences is shown in Figure 8a. The tide level difference at each node was obtained from the average tide level over a spring-neap cycle for 0.08   TCP extraction produced a very noticeable reduction in the tidal current. Figure 9 demonstrates a time series of the depth-averaged current for The maximum instant tidal current exceeded 1.1 m/s when the turbines are absent in the coastal waters; however, it was reduced to 68% due to the presence of tidal turbines (i.e.,

Influence of Tidal Current Power Extraction on Hydrodynamics
TCP extraction may alter the hydrodynamics (e.g., tide level and currents) over an area where the turbines are deployed. To understand the impacts of turbine operation on tidal dynamics, time series and spatial differences of tide level and currents were compared for an area with and without turbines. Due to 0.08 t C = produced the maximum mean TCP, Figure 7 depicts a tidal water level time series for 0 t C = and 0.08 t C = at the turbine site. The turbine-induced tide level change is minor; the difference is difficult to distinguish in Figure 7. The spatial distribution of tide level differences is shown in Figure 8a. The tide level difference at each node was obtained from the average tide level over a spring-neap cycle for 0.08   TCP extraction produced a very noticeable reduction in the tidal current. Figure 9 demonstrates a time series of the depth-averaged current for The maximum instant tidal current exceeded 1.1 m/s when the turbines are absent in the coastal waters; however, it was reduced to 68% due to the presence of tidal turbines (i.e.,   TCP extraction produced a very noticeable reduction in the tidal current. Figure 9 demonstrates a time series of the depth-averaged current for C t = 0 and C t = 0.08 at the turbine site. The maximum instant tidal current exceeded 1.1 m/s when the turbines are absent in the coastal waters; however, it was reduced to 68% due to the presence of tidal turbines (i.e., C t = 0.08). It is worth noting that the decrease in the time average tidal current was 77% for the comparison between C t = 0 and C t = 0.08 (as shown in Figure 6). Figure 10a exhibits the spatial distribution difference of the average tidal current over a spring-neap cycle for C t = 0.08 minus C t = 0. The extent of the influence is 0.  Figure 10b). The installment of turbines not only significantly retards the tidal current but also affects a certain area adjacent to the TCP extraction site.  Figure 10b). The installment of turbines not only significantly retards the tidal current but also affects a certain area adjacent to the TCP extraction site.

Influence of Sea-Level Rise on Tidal Current Power Output
According to the Intergovernmental Panel on Climate Change (IPCC) AR5 report, the sea levels could rise 0.26 m to 0.54 m for RCP 2.6; 0.32 m to 0.62 m for RCP 4.5; 0.33 m to 0.62 m for RCP6.0; and 0.45 m to 0.81 m for RCP8.5 by 2100 [48]. Therefore the effects of SLR on TCP output were assessed by conducting the scenario simulations for SLR 0.25 m, SLR 0.5 m, SLR 0.75 m and SLR 1.0 m. We assumed that the tidal dynamics outside of the study region do not change. All boundary and bathymetry conditions are kept constant except for increasing the water depth by 0.25 m, 0.5 m, 0.75 m and 1.0 m over the entire computational domain and open boundaries. This straightforward approach has been widely adopted to investigate the effects of SLR on coastal hydrodynamics and tidal energy resources [19,44,49,50]. TCP outputs were also calculated by Equations (6) and (7) with an incremental t C , a cut-in speed of 0.5 m/s and t A of 35762.44 m 2 ; the depth-averaged tidal currents were obtained from the outputs of SCHISM-2D. Figure 11 illustrates   Figure 10b). The installment of turbines not only significantly retards the tidal current but also affects a certain area adjacent to the TCP extraction site.

Influence of Sea-Level Rise on Tidal Current Power Output
According to the Intergovernmental Panel on Climate Change (IPCC) AR5 report, the sea levels could rise 0.26 m to 0.54 m for RCP 2.6; 0.32 m to 0.62 m for RCP 4.5; 0.33 m to 0.62 m for RCP6.0; and 0.45 m to 0.81 m for RCP8.5 by 2100 [48]. Therefore the effects of SLR on TCP output were assessed by conducting the scenario simulations for SLR 0.25 m, SLR 0.5 m, SLR 0.75 m and SLR 1.0 m. We assumed that the tidal dynamics outside of the study region do not change. All boundary and bathymetry conditions are kept constant except for increasing the water depth by 0.25 m, 0.5 m, 0.75 m and 1.0 m over the entire computational domain and open boundaries. This straightforward approach has been widely adopted to investigate the effects of SLR on coastal hydrodynamics and tidal energy resources [19,44,49,50]. TCP outputs were also calculated by Equations (6) and (7) with an incremental t C , a cut-in speed of 0.5 m/s and t A of 35762.44 m 2 ; the depth-averaged tidal currents were obtained from the outputs of SCHISM-2D. Figure 11 illustrates the average TCP outputs corresponding to different t C values. The

Influence of Sea-Level Rise on Tidal Current Power Output
According to the Intergovernmental Panel on Climate Change (IPCC) AR5 report, the sea levels could rise 0.26 m to 0.54 m for RCP 2.6; 0.32 m to 0.62 m for RCP 4.5; 0.33 m to 0.62 m for RCP6.0; and 0.45 m to 0.81 m for RCP8.5 by 2100 [48]. Therefore the effects of SLR on TCP output were assessed by conducting the scenario simulations for SLR 0.25 m, SLR 0.5 m, SLR 0.75 m and SLR 1.0 m. We assumed that the tidal dynamics outside of the study region do not change. All boundary and bathymetry conditions are kept constant except for increasing the water depth by 0.25 m, 0.5 m, 0.75 m and 1.0 m over the entire computational domain and open boundaries. This straightforward approach has been widely adopted to investigate the effects of SLR on coastal hydrodynamics and tidal energy resources [19,44,49,50]. TCP outputs were also calculated by Equations (6) and (7) with an incremental C t , a cut-in speed of 0.5 m/s and A t of 35762.44 m 2 ; the depth-averaged tidal currents were obtained from the outputs of SCHISM-2D. Figure 11 illustrates the average TCP outputs corresponding to different C t values. The maximum average TCP output is evaluated to be 62.

Discussion and Future Work
The effect of TCP extraction on the change in tide level is difficult to determine in the time series and spatial difference figures (Figures 7 and 8). To quantify TCP extraction impact on tide level change, harmonic analysis was conducted to extract the main eight tidal components for (with maximum mean TCP output) at the turbine site. Table 2 lists the amplitudes and phases of eight tidal constituents for . The maximum changes in amplitude and phase caused by turbines are only 0.73 cm and 0.77° for the K2 tide, respectively. The change of tide level for TCP extraction is ignored in both qualitative and quantitative analyses. Yang and Wang [51] compared the tide level changes of a tidal channel with and without turbines. They found that TCP extraction had a minor influence on the tidal elevations. Defne et al. [2] assessed the effects of TCP extraction on water level for an estuarine area, they concluded that both the decrease in the maximum water level and the increase in the minimum are kept below 0.05 m. Chen et al. [38] estimated the influences of TCP extraction on water level for the Taiwan Strait using the momentum sink approach. Their study results found that TCP extraction at a specific farm could significantly affect the tidal current, but the change in tidal amplitude and tidal phase could be neglected. Plew and Stevens [52] also suggested that introducing turbines in the numerical model has only a minor influence on the difference in water level through the Tory Channel. This phenomenon might be due to the water level simulations are mainly determined by the tide elevations specified at open boundaries.  According to the numerical experiments conducted in Section 3.4, the maximum average TCP was increased with SLR because a higher sea level produces faster tidal currents (as shown in Figure 12). The momentum equations in Equations (2) and (3) can be used to explain this phenomenon. The last term on the right-hand side of Equations (2) and (3) [53] estimated the impacts of SLR on tidal dynamics for the Ria de Aveiro lagoon in Portugal, using the ELCIRC 2D model. Chen and Liu [9] assessed the influences of SLR on tidal energy output for the Taiwan Strait. The same findings were also found in their studies, suggesting that SLR produces faster tidal currents and consequently increases TCP output. Uehara et al. [54] used a two-dimensional model to investigate tidal and tide-dependent changes in the northwest European shelf seas over the last twenty thousand years. They adopted two approaches to consider the effect of future ocean-tide variability on the open boundary of the numerical model. One is fixed to the present state; the other is variable according to a global paleotidal model. They found that the M2 amplitude only increases by 1.0 m (0.9-1.9 m) in the west of the Celtic Sea through twenty thousand years. Although more realistic changes in shelf tides during ten thousand years can be obtained in this way, the variation of tidal amplitude might be small 100 years from now. It would be sensible to assume that the change of tidal amplitudes from 2017 to 2100 is ignorable. Therefore, we compared the TCP output between present-day sea level and future SLR using the same boundary conditions in the model. However, modeling assessment of future hydrodynamics might be more faithful if accurate projections of future tide variations for model boundary conditions are available.
Pelling and Green [19] conducted a series of numerical experiments to assess the SLR influence on TCP output for the Gulf of Maine. In their study, SLR was implemented according to two approaches. The first approach was to increase the water depth over the entire model domain and allow low-lying elements on the land to be flooded (referred to as the 'flood run'). The second method was the same as the first except that no area was allowed to be flooded (referred to as the 'no-flood run'). They found that the no-flood run produced higher maximum power output than the flood run for the Minas Passage. This is because TCP became dissipated when new shallow-water areas were created in the flood run. The difference in maximum TCP output is only 8.7% between these two approaches for the Minas Passage. Ward et al. [49] used the flood run

Discussion and Future Work
The effect of TCP extraction on the change in tide level is difficult to determine in the time series and spatial difference figures (Figures 7 and 8). To quantify TCP extraction impact on tide level change, harmonic analysis was conducted to extract the main eight tidal components for C t = 0 (without TCP output) and C t = 0.08 (with maximum mean TCP output) at the turbine site. Table 2 lists the amplitudes and phases of eight tidal constituents for C t = 0 and C t = 0.08. The maximum changes in amplitude and phase caused by turbines are only 0.73 cm and 0.77 • for the K 2 tide, respectively. The change of tide level for TCP extraction is ignored in both qualitative and quantitative analyses. Yang and Wang [51] compared the tide level changes of a tidal channel with and without turbines. They found that TCP extraction had a minor influence on the tidal elevations. Defne et al.
[2] assessed the effects of TCP extraction on water level for an estuarine area, they concluded that both the decrease in the maximum water level and the increase in the minimum are kept below 0.05 m. Chen et al. [38] estimated the influences of TCP extraction on water level for the Taiwan Strait using the momentum sink approach. Their study results found that TCP extraction at a specific farm could significantly affect the tidal current, but the change in tidal amplitude and tidal phase could be neglected. Plew and Stevens [52] also suggested that introducing turbines in the numerical model has only a minor influence on the difference in water level through the Tory Channel. This phenomenon might be due to the water level simulations are mainly determined by the tide elevations specified at open boundaries. According to the numerical experiments conducted in Section 3.4, the maximum average TCP was increased with SLR because a higher sea level produces faster tidal currents (as shown in Figure 12). The momentum equations in Equations (2) and (3) can be used to explain this phenomenon. The last term on the right-hand side of Equations (2) and (3) shows that the horizontal velocities, u(x, y, t) and v(x, y, t), are proportional to the total water depth (H = η + h) when the wind shear stresses (τ sx and τ sy ) are excluded in the model. In other words, deeper water depth (i.e., SLR) leads to faster tidal currents if the bottom drag coefficient remains constant. Lopes and Dias [53] estimated the impacts of SLR on tidal dynamics for the Ria de Aveiro lagoon in Portugal, using the ELCIRC 2D model. Chen and Liu [9] assessed the influences of SLR on tidal energy output for the Taiwan Strait. The same findings were also found in their studies, suggesting that SLR produces faster tidal currents and consequently increases TCP output.
Uehara et al. [54] used a two-dimensional model to investigate tidal and tide-dependent changes in the northwest European shelf seas over the last twenty thousand years. They adopted two approaches to consider the effect of future ocean-tide variability on the open boundary of the numerical model. One is fixed to the present state; the other is variable according to a global paleotidal model. They found that the M 2 amplitude only increases by 1.0 m (0.9-1.9 m) in the west of the Celtic Sea through twenty thousand years. Although more realistic changes in shelf tides during ten thousand years can be obtained in this way, the variation of tidal amplitude might be small 100 years from now. It would be sensible to assume that the change of tidal amplitudes from 2017 to 2100 is ignorable. Therefore, we compared the TCP output between present-day sea level and future SLR using the same boundary conditions in the model. However, modeling assessment of future hydrodynamics might be more faithful if accurate projections of future tide variations for model boundary conditions are available.
Pelling and Green [19] conducted a series of numerical experiments to assess the SLR influence on TCP output for the Gulf of Maine. In their study, SLR was implemented according to two approaches. The first approach was to increase the water depth over the entire model domain and allow low-lying elements on the land to be flooded (referred to as the 'flood run'). The second method was the same as the first except that no area was allowed to be flooded (referred to as the 'no-flood run'). They found that the no-flood run produced higher maximum power output than the flood run for the Minas Passage. This is because TCP became dissipated when new shallow-water areas were created in the flood run. The difference in maximum TCP output is only 8.7% between these two approaches for the Minas Passage. Ward et al. [49] used the flood run method to investigate the impact of SLR on tidal dynamics for the European shelf; they also suggested that different way to implement SLR in the model shown a significant effect on response of tide. Although the flood run approach is much more realistic than the no-flood run approach, we can only adopt the latter due to the scarcity of accurate topography data and land-use data for Kinmen Island and Lieyu Island. In other words, we can suppose that an actual average maximum TCP might be 41.55 kW (45.51 × (1 − 0.087) = 41.55 kW) in the Kinmen coastal waters if the flood run were used.
Due to the amount of TCP output is small and semi-diurnal intermittency of tidal current at the turbine deployment site in the Kinmen coastal waters; it might be possible to conduct storage into the supergrid [15,55]. Besides, Neill et al. [56] proposed that optimization of the TCP and finding the optimal site for parallel development for relative low number of arrays are import for large scale exploitation of TCP and should be conducted in the future work. In addition, Lewis et al. [57] found that the wave-current interaction could potentially alter tidal currents and should be included in hydrodynamic models in future studies.
Using a combination of Equations (5) and (7) to obtain the number of turbines is too simple when maximum mean TCP output was given. Myers and Bahaj [58] and Mestres et al. [59] suggested that the separation between both rows is taken as 10 times the dimeter of the tidal turbine rotor and the inter-lateral turbine separation 2.5 times the diameter. The configuration of tidal turbine array should be carefully considered before the construction of tidal turbine farm; however this requires further examination, which is beyond the scope of this paper.

Conclusions
The hydrodynamics in the coastal waters of Kinmen Island have been simulated using a state-of-the-art unstructured-grid depth-integrated numerical model, SCHISM-2D. The simulations were verified with tide level and depth-averaged tidal current measurements, showing a reasonable agreement with the observations. The validated model was then applied to evaluate TCP resources and appropriate sites for deploying tidal turbines around the coastal waters of Kinmen Island. A narrow channel between Kinmen and Lieyu was found to be a suitable area for turbine deployment because of its higher mean tidal current. The site has an area of 0.036 km 2 , depth greater than 20 m and average velocity of 0.43 m/s, making it suitable as a tidal current turbine farm. The bottom friction approach was employed to calculate the mean TCP over a spring-neap cycle. The maximum average TCP output is estimated as 45.51 kW when the additional turbine friction coefficient C t = 0.08, cut-in speed is 0.5 m/s and average velocity is 0.33 m/s. The impacts of TCP extraction on hydrodynamics were also estimated. The TCP extraction effect can be ignored for tide level but shows a significant reduction in tidal currents. We found that the installation of turbines not only retarded tidal currents but also affected a certain area near the TCP extraction site. Four SLR scenarios, SLR 0.25 m, SLR 0.5 m, SLR 0.75 m and SLR 1.0 m, were conducted to investigate the influence of SLR on TCP output. Compared to the present-day sea-level condition, the annual TCP output increased to 1.52 MW, 2.01 MW, 2.48 MW and 2.97 MW for SLR 0.25 m, SLR 0.5 m, SLR 0.75 m and SLR 1.0 m, respectively, due to the faster tidal current. Acceptable simulations of tidal dynamics for the coastal waters of Kinmen Island were accomplished in this study, however, there are still many uncertainties associated with tidal current simulations, and TCP calculation should be further considered in the future.