Numerical Modeling of the Impact of Sea Level Rise on Tidal Asymmetry in Hangzhou Bay

<jats:p>The tidal asymmetry under the action of sea level rise in Hangzhou Bay. Coastlines dominated by the tide are not only directly affected by the rise of the mean sea level but also by the tidal dynamics. The computational domain of the hydrodynamic model covers the entire Hangzhou Bay and takes into account the feedback between the tidal motion and the erodible bottom. Its main application fields include: the simulation of different sea level rise (SLR) rates, the interaction between tidal duration and skewness and the interaction between tidal range and astronomical tide. The results on tidal asymmetry in Hangzhou Bay is a systematic process consisting of four aspects: Firstly, the tide increase, which is affected by the sea level rise, is between 25% and 50%. Secondly, the value of the sea level rise is about two times the added value of the tidal range, and the tidal range increased to the left side of the tidal wave propagation direction, which accelerated the propagation velocity. Thirdly, the sea level rise amplified the M2 tidal amplitude and delayed the M2 tidal phase in the inner bay, which reduced about 50%. Finally, the change of the tidal range caused by tide level had the same magnitude as the change caused by the mean sea level rise. The purpose of this study was to emphasize the importance of the predictions of the response area affected by tidal asymmetry based on the action of the sea level rise in Hangzhou Bay.</jats:p>

section_content: Tidal asymmetry is a phenomenon of tidal wave deformation that results in an unequal duration of the rise and fall of the tidal amplitude. A tidal wave propagates from the open sea to the coast, and the water depth, runoff and other factors reduce the tidal wave energy. This deformation has an important impact on the sediment movement, topographic variation and material migration. On the role of diurnal tides in contributing to tidal asymmetries in areas of predominantly semi-diurnal tide, the M 2 tidal wave in the open sea and the M 4 and M 6 tidal waves generated near the shore are deformed after superposition, presenting an asymmetric form [1, 2]. The amplitude ratio (M 4 /M 2 ) is usually used to describe the deformation degree of tidal waves, and the relative phase (2M 2 -M 4 ) is used to describe the asymmetric characteristics [3, 4]. Friedrichs [5] used the indicators to quantitatively study the tidal and tidal asymmetry in an estuary. Blanton [6] added the analysis of the M 6 tidal component when studying different types of estuaries. K 1, O 1 and other tidal components replace M 2 and M 4 and other tidal components in the diurnal tide, but the research method in the semidiurnal tide sea area is no longer applicable. Ranasinghe [7] studied the energy transfer between K 1, O 1, M 2 and M 4 tidal components using probability distribution and spectral methods. Elgar [8] and Nidzieko [9] used the "skewness" and "asymmetry" in statistics to study the asymmetry of tides and tidal currents. Song [10] further clarified the calculation method of tidal asymmetry based on Nidzieko's research and studied the tidal asymmetry of multiple stations around the world. Friedriches and Aubrey proposed to use the amplitude ratio of the M 4 component tide and M 2 component tide and quantitatively calculate the direction of the tidal asymmetry. At the same time, the duration of the flood and ebb tide was used to judge which was dominant [11]. Due to the implications in response to human activities in shallow waters, many scholars have applied this method to the study of the Yangtze River Estuary, Bohai Sea, Yellow Sea and so on [12] [13] [14] [15]. The M 2 and M 4 interaction is the main contributor to asymmetry and also the main reason for tidal asymmetry in Hangzhou Bay. 

Sea level rise affects tides on all spatial scales, ranging from the global scale, i. e., oceans [16], to the regional scale, i. e., continental shelves such as the European and Asia continental shelf [17] [18] [19], and the local scale, i. e., coastal waters such as the German Wadden Sea [20] and estuaries [21, 22]. The purpose of our study is bipartite. Firstly, we researched the tidal response to four different SLR scenarios of 0. 44, 0. 56, 0. 68 and 0. 77m in 2100. The global mean sea level rise has accelerated since the late 1960s [23], with an average rate of 2. 3 (1. 6-3. 1) mmyr-1 over the period 1971-2018 and increasing to 3. 7 (3. 2-4. 2) mmyr-1 over the period 2006-2018 [24]. It represents a regional projection of the sea level contributions in 2100 relative to 1995-2014 for SSP5-8. 5 and SSP1-2. 6 [25, 26]. Secondly, we analyzed the hydrodynamic feedbacks due to SLR by 2100. Therefore, we compared the effects of SLR variations in Hangzhou Bay to the effects of tidal asymmetry. The results were mainly analyzed in terms of skewness and the maximum tidal range to get the tidal characteristics. What is more, we discuss the aim of the M 4 amplitude and M 2 amplitude, as well as the ratio between both, as these parameters are the most suitable to demonstrate changes on tidal asymmetry in Hangzhou Bay. The focus on tidal asymmetry of the current velocity is largely motivated by its significance in the residual transport of sediment and the morphologic development of Hangzhou Bay; hence, the formation mechanism of tidal asymmetry is discussed in Hangzhou Bay. 

The data validity and model setup are briefly described below. What follows in Section 2 is the model validation and analysis methods. The physical processes for applied scenarios (sea level rise) are discussed in Section 3, and the causes of tidal asymmetry are analyzed. Finally, the conclusions are given in Section 4.
section_name: Introduction

section_name: Methods

section_content: Hangzhou Bay, located on the east coast of China, is a typical wide, shallow and funnel-shaped bay, and the strongest tide-dominated coasts often occur to China's east. The tidal changes of Hangzhou Bay are complex and influence the East China Sea [27] with the most significant features of strong tidal currents, large tidal ranges and high suspended sediment concentrations [28]. The hydrodynamic model applied in this study was MIKE21, which is a process-oriented numerical model developed by the DHI [29]. The two-dimensional hydrodynamic model solves the shallow water equations on an unstructured grid consisting of triangles [30]. Its unstructured grid has the advantage of realizing an effective spatial resolution in a complex bathymetric environment, because it can adapt to the morphologic structure and allows highly variable model resolutions. Therefore, it is particularly suitable for simulating the tidal dynamic characteristics of Hangzhou Bay. It is important for the simulation of nonlinear hydrodynamic effects in tidally dominated, shallow coastal waters such as Hangzhou Bay [20]. The Hangzhou Bay model employs the Smagorinsky formulation for horizontal eddy viscosity, and the Smagorinsky coefficient Cs is 0. 28 [31]. 

The study area, as shown in Figure 1a, covers Hangzhou Bay, the Yangtze River estuary and Zhoushan Islands. The data used for validating the numerical modeling results of the tide level and tidal currents (current velocity and current direction) were acquired from the measured data collection [32]. The data collected contained tidal dynamics of approximately 57 tidal cycles with 1 h intervals throughout September 2010. The accuracy of the model can be ensured. Details of the 2D numerical model verification included the tidal component, tide level and tidal current [33]. Jiangyin and Cangqian were limited by the flow in the model established, and the open sea is bounded by the tide level. According to the 2020 water table released on 9 December 2021 (http://www. mwr. gov. cn/sj/tjgb/sqnb/), the annual average river discharges in Datong and Cangqian are 28,850 m 3 /s and 952 m 3 /s, respectively. The open sea boundaries are specified as time-varying tide levels derive from the hydrodynamic model covering the entire East China Sea and its adjacent seas, whose open boundary is driven by the global tidal model [22]. In this paper, the effectiveness of tidal dynamics for setting up conditions along the open sea boundaries of coastal models is verified. The numerical simulation of the impact of the sea level rise on estuarine dynamics usually uses the control of the tide level and the elevation of the tide level boundary to achieve the sea level rise. This method was applied to the Yangtze River estuary [22]. Sea level rise is introduced by adding a constant value to the water levels at the open boundary of the estuary model in Hangzhou Bay. In this way, shallow water effects on the estuary and its potential changes with the sea level rise are included in the forcing of the regional Hangzhou Bay model. The simulated SLR scenario of the two-dimensional hydrodynamic model located in the northwestern East China Sea is 0. 44-0. 7 m. Using these projections as input to construct projections of the GMSL change with IPCC AR5 methods predict that the 95th percentile of GMSL change in 2100 only increases by 3-7 cm [26]. This shows that the model in this paper is consistent with that of the scholars in SLR scenarios. [22]. In this paper, the effectiveness of tidal dynamics for setting up conditions along the open sea boundaries of coastal models is verified. The numerical simulation of the impact of the sea level rise on estuarine dynamics usually uses the control of the tide level and the elevation of the tide level boundary to achieve the sea level rise. This method was applied to the Yangtze River estuary [22]. Sea level rise is introduced by adding a constant value to the water levels at the open boundary of the estuary model in Hangzhou Bay. In this way, shallow water effects on the estuary and its potential changes with the sea level rise are included in the forcing of the regional Hangzhou Bay model. The simulated SLR scenario of the two-dimensional hydrodynamic model located in the northwestern East China Sea is 0. 44-0. 7 m. Using these projections as input to construct projections of the GMSL change with IPCC AR5 methods predict that the 95th percentile of GMSL change in 2100 only increases by 3-7 cm [26]. This shows that the model in this paper is consistent with that of the scholars in SLR scenarios.
section_name: Model Setup and Data Availability

section_content: In the following sections, the validity of the Hangzhou Bay model is analyzed. Figure 2 and Figure 3 show a times series of water levels from model simulations and measurements at twelve different stations within the Hangzhou Bay covering a full spring-neap tidal cycle in September 2010. 

The model well captures the propagation of tidal waves throughout Hangzhou Bay during the whole process, including the increasing tidal range towards the head bay and
section_name: Model Validation

section_content: In the following sections, the validity of the Hangzhou Bay model is analyzed. Figures 2 and 3
section_name: Model Validation

section_content: This section is divided by subheadings. It provides a concise description of the experimental results and their interpretations, as well as the experimental conclusions that can be drawn.
section_name: Methodology

section_content: Following Wu et al. (2018), when the amplitude of the S2 tide is much smaller than that of the M2 tide, as in the case of Hangzhou Bay, is the criterion in Equation (1). According to the changes of the tide level in the sea area of time and space, the tidal properties are divided into semidiurnal and diurnal tidal components. The commonly used tidal type discriminant is as follows:

where H is the amplitude of the harmonic constituent in the subscripts [35]. For regular semidiurnal tides, F = 0-0. 5, for regular semidiurnal tides, F = 0. 5-2. 0, for irregular diurnal tides, F = 2. 0-4. 0 and for regular diurnal tides, F > 3. 0 [36]. 

According to the tidal properties of Hangzhou Bay, the maximum possible tidal dif- The model well captures the propagation of tidal waves throughout Hangzhou Bay during the whole process, including the increasing tidal range towards the head bay and the more stable tidal signal towards the coast. The overall consistency between the simulation and measurement results has been verified. Table 1 shows a good agreement in both the magnitude and phase, and the error of the outer bay is minimal, with all skill values of the tide levels exceeding 0. 94 (Table 2 ). The skill scores (0. 82-0. 99) for several velocities are slightly less than those for the tide level (Figure 3 ); however, the model has captured that the model validation case and the model performance can be regarded as excellent. Therefore, this robust model gives us confidence in the SLR effects forecast. Overall, it can be concluded that the model validity is sufficient to correctly determine regional-scale changes of tidal dynamics in response to changes in the sea level. A high verification accuracy ensures precision of the digital model [34].
section_name: Calculation of Tidal Asymmetry

section_content: This section is divided by subheadings. It provides a concise description of the experimental results and their interpretations, as well as the experimental conclusions that can be drawn. 2018 ), when the amplitude of the S 2 tide is much smaller than that of the M 2 tide, as in the case of Hangzhou Bay, is the criterion in Equation (1). According to the changes of the tide level in the sea area of time and space, the tidal properties are divided into semidiurnal and diurnal tidal components. The commonly used tidal type discriminant is as follows:

where H is the amplitude of the harmonic constituent in the subscripts [35]. For regular semidiurnal tides, F = 0-0. 5, for regular semidiurnal tides, F = 0. 5-2. 0, for irregular diurnal tides, F = 2. 0-4. 0 and for regular diurnal tides, F > 3. 0 [36]. 

According to the tidal properties of Hangzhou Bay, the maximum possible tidal difference is calculated as follows [37] :

The left side of the equation is the maximum possible tidal range. The right side of the equation is where H is the amplitude of the harmonic constituent in the subscripts.
section_name: Methodology

section_content: Tidal skewness is a statistical approach based on the calculations of the probability density function of the water level proposed by Nidzieko and further researched by Song [10]. In this method, the most remarkable feature can be determined by the strength and relative contribution of each of combinations in the total tidal asymmetry based on the amplitudes, frequencies and phases of the tidal constituents. Therefore, the tidal skewness analyzed from the combination of two tidal constituents is given [38] in Equation (3) and the contribution to the total tidal skewness by the triad combination of tidal constituents is obtained [39] in Equation (4) as:

, 2ω 1 = ω 2 (3)

,

where a i, ω i and ϕ i are the amplitude, frequency and phase of the tidal constituent, respectively. A positive β 2 means that the rising tidal is dominant, and a negative β 3 means that the falling tidal is dominant. According to Song's theory, under the condition of secondorder nonlinearity, only when the frequency of two tidal components meets 2ω 1 = ω 2, the frequency of 2 or 3 tidal components meets ω 1 + ω 2 = ω 3, and their interaction can cause tidal asymmetry.
section_name: Tidal Skewness

section_content: In-site measured water elevation and the tidal current were both used to model the validation. To evaluate the model, the parameters were used to quantify the differences between the measured data and simulation results and computed as follows:

where n is the number of variable values, X m represents the time-varying model results and X 0 and X 0 are the time-varying values of the observed results and the time mean values, respectively. The performance of the model depends on the value of S, which is classified as excellent (S > 0. 65), very good (0. 5 < S < 0. 65), good (0. 2 < S < 0. 5) or poor (S < 0. 2) [40].
section_name: Skill Parameters for Quantifying Model Verification

section_content: In this section, the tidal characteristics of the tidal constituents change with the sea level rise, and the influence of skewness on the tidal response in Hangzhou Bay is illustrated. The impacts on the tidal asymmetries by SLR changes at the inner bay and outer bay are elucidated.
section_name: Results and Discussion

section_name: Effects of Tidal Harmonics

section_content: The astronomical tides waves (e. g., M 2 ) that can generate shallow water and compound constituents (e. g., M 4, M 6 and others) spread to the coast shelves, estuaries and bays. Due to the superposition of these components, the tidal current and water level are distorted from their sinusoidal forms, which give rise to a tidal asymmetry [2]. Observations made in the main inner bay indicate that the model is presently ebb-dominant with regards to both the relatively short duration of the ebb tide and also the shorter ebb velocities. The tidal asymmetry plays a major role in shallow estuaries and tidal creeks [41]. There is also an asymmetric pattern for the tide caused by the deformation of the dominant astronomical tidal constituents [42], M 2 and M 4, due to the nonlinear interaction of tide with the irregular estuarine geometry and the erodible bottom. Therefore, the astronomical tidal constituents are considered to discuss the tidal characteristics. 

This section describes the effects of the trumpet-shaped estuary and the rise of the bed elevation applied to the tidal characteristics (Table 3 ). The area in which the tidal characteristics occur consists of an average high to low water level ratio (H max /H min ), and the mean tidal range of Hangzhou Bay quickly rises to the north bank, but the average low water level H min is smaller at the south bank. The mean tidal range irregularities of the bank bottom average 2 to 3 m relief on the western bank (Table 3 ) but can be over 5 m in localized areas at the head of the northern bank. Then, the average flood to tide water level ratio is as high as 1. 4, which shows the flood dominance. Both the northern and western banks of the low water level rise, while the flood tide range decreases. The duration of the flood tide exceeds the duration of the ebb tide, resulting in ebb dominance at Hangzhou Bay. The M 2 and M 4 interaction was the main contribution to tidal asymmetry. The largest constituent was M 2, followed by K 1, O 1 and S 2. The largest shallow water constituent was MK 3, followed by M 4, MS 4, MN 4 and M 6 [43]. We confirmed the tidal constituents in Hangzhou Bay are mainly dominated by semidiurnal tides (e. g., M 2 and S 2 ) and diurnal tides (e. g., K 1 and O 1 ), according to the conclusion by Lafta in the Arabian Gulf. The influence of runoff and the river mouth section shrinks at the head of the Qiantang River Estuary in Hangzhou Bay. It is a much more complicated process of the main tidal constituent within the flood season variations. The amplitude of M 2 slightly decreases from upstream toward downstream. The area in which the main tidal constituents occur consists of the amplitude and phase along the direction from the north bank to the south bank. The amplitude of the four main tidal constituents decreases downstream of the bay along the north bank. The duration of the flood tide varies greatly with the depth and position along the estuary. The phase of the eight tidal gauges decreases along the ebb direction and the duration of the ebb tide exceed the duration of the flood tide. The astronomical tide around Hangzhou Bay is regularly semidiurnal, except near ZH, where it is irregularly semidiurnal. Consequently, the amplitude of the tidal wave transmits from the outer bay to the inner bay, significantly increased by the astronomical tide. The amplitude of the semidiurnal tidal components on the north bank slightly increases and, on the south bank, decreases.
section_name: Tidal Characteristics

section_content: The tidal asymmetry was investigated furtherly by the tidal skewness metrics [43]. Nidzieko proposed that tidal asymmetry be quantified via skewness [38]. In this paper, we extended Nidzieko's analysis to show the skewness distribution. The skewness of H1-H20 stations in Hangzhou Bay (Figure 4 ) and their negative values represent the advantage of the ebb tide. The skewness varies greatly with the coastline shape and estuarine water depth, and tidal gauges have the advantage of the flood tide and a different skewness. Firstly, the skewness gradually increases from the north to south bank in the inner bay. The significant dominance of the flood duration (approximately 7 h per tidal cycle) and tidal asymmetry lead to flood dominance in the inner bay, and the skewness gradually decreases from the north to south bank in the outer bay. Considering the coastline and bathymetrics, Benno [20] proposed that the influence of bathymetric varieties has the same magnitude as the influence of the sea level rise, and the tidal response to SLR and the coastline change rise due to the inclinations of the SZP and ZH areas. Flood dominance has an obvious upward trend from the outer bay to the head bay, and the tide advantage gradually weakens from west to east. It was analyzed that the reason is that the shallower coastline and terrain lead to an increase of the skewness, which causes an increase of the
section_name: Effects of Skewness

section_content: The tidal asymmetry was investigated furtherly by the tidal skewness metrics [43]. Nidzieko proposed that tidal asymmetry be quantified via skewness [38]. In this paper, we extended Nidzieko's analysis to show the skewness distribution. The skewness of H1-H20 stations in Hangzhou Bay (Figure 4 ) and their negative values represent the advantage of the ebb tide. The skewness varies greatly with the coastline shape and estuarine water depth, and tidal gauges have the advantage of the flood tide and a different skewness. Firstly, the skewness gradually increases from the north to south bank in the inner bay. The significant dominance of the flood duration (approximately 7 h per tidal cycle) and tidal asymmetry lead to flood dominance in the inner bay, and the skewness gradually decreases from the north to south bank in the outer bay. Considering the coastline and bathymetrics, Benno [20] proposed that the influence of bathymetric varieties has the same magnitude as the influence of the sea level rise, and the tidal response to SLR and the coastline change rise due to the inclinations of the SZP and ZH areas. Flood dominance has an obvious upward trend from the outer bay to the head bay, and the tide advantage gradually weakens from west to east. It was analyzed that the reason is that the shallower coastline and terrain lead to an increase of the skewness, which causes an increase of the skewness from the open sea to the bay. Therefore, the coastline and water depth have the greatest impact on the SLR scenarios. 

The skewness of the H1-H20 stations is gradually increased from the outer bay to the inner bay (Figure 5 ), and the tide asymmetry increases the advantage of the flood tide along the way. In the four SLR scenarios, the skewness of the GP-CX section significantly increases from the north to south bank and the increase of the skewness is as high as 115% at the ZP section (Figure 5a, b ). The skewness values sharply become larger, though there was a turning point at the H9 station. The skewness of the JSZ-SZP section significantly decreased from the north to south bank in different SLR scenarios, and the turning point was at H14 (Figure 5c ), the skewness increased from the H14 to H15 station in due to the sea level rise (RSLR) and two SLR scenarios (0. 44 m SLR and 0. 56 m SLR). The overall skewness decreased, in short, by 4% for four SLR scenarios (Figure 5d ), and the skewness reached the minimum at the turning point of H18. The tidal calendar time difference of H18 was 0. 1 h with regards to the 0. 77 m SLR scenario. The study showed that the skewness possesses a decreasing trend that nevertheless increases from west to east in The skewness of the H1-H20 stations is gradually increased from the outer bay to the inner bay (Figure 5 ), and the tide asymmetry increases the advantage of the flood tide along the way. In the four SLR scenarios, the skewness of the GP-CX section significantly increases from the north to south bank and the increase of the skewness is as high as 115% at the ZP section (Figure 5a, b ). The skewness values sharply become larger, though there was a turning point at the H9 station. The skewness of the JSZ-SZP section significantly decreased from the north to south bank in different SLR scenarios, and the turning point was at H14 (Figure 5c ), the skewness increased from the H14 to H15 station in due to the sea level rise (RSLR) and two SLR scenarios (0. 44 m SLR and 0. 56 m SLR). The overall skewness decreased, in short, by 4% for four SLR scenarios (Figure 5d ), and the skewness reached the minimum at the turning point of H18. The tidal calendar time difference of H18 was 0. 1 h with regards to the 0. 77 m SLR scenario. The study showed that the skewness possesses a decreasing trend that nevertheless increases from west to east in Hangzhou Bay. The skewness values change from decreasing to increasing along the north bank to the south in the outer bay. The SLR scenarios keep up with the trend of the RSLR scenario increased range of 1~84%. The skewness gradually increases with the rise of the sea level.
section_name: Effects of Skewness

section_content: In this section, SLR enlarges the tidal range and accelerates the wave propagation, resulting in an earlier tidal arrival. The influence of SLR on the tidal range in Hangzhou Bay is discussed. The tidal range changes are included by different SLR scenarios (Figure 6 ). As the sea level rises, the tidal range uniformly increases in the inner and outer bays. The shallow water tidal constituents generated in coastal waters are of great importance to the tidal wave morphology and tidal range. For the SLR scenario by 0. 4m, on the one hand, the tidal wave passes from the outer sea to the Zhoushan Islands and then enters Hangzhou Bay. On the other hand, the tidal range gradually increases from the outer reaches to the upper reaches of the bay. The tidal range will be 0. 16 m~0. 24 m in the LCG-GP section, and the tidal range is as high as 0. 96m near AD. In the 0. 56 m SLR scenario, most of the sea tidal range will be 0. 24m~0. 40m in Hangzhou Bay. The tidal range will increase by 0. 08m at closer to 33. 3%~50%, reaching 0. 12~0. 2 m from the outer bay to the LCG-ZH section (Figure 6a, b ). The tidal range has increased near the GP-CX section by 50%. In the 0. 68m SLR scenario, the overall increase of the tidal range values is 0. 4m~0. 5 and expanded by 25% (Figure 6c ), forming a sharp contrast. In the 0. 77m SLR scenario, the tidal range will attain 0. 4 m~0. 5 m in the inner and outer bays. The maximum tidal range appears near AD and GP (Pan et al. 2009), and the increase of the tidal range will be less than the appreciation of the sea level. In a word, the tidal range will increase with the rise of sea level by 25%~50% in the four SLR scenarios (0. 44, 0. 56, 0. 68 and 0. 77m), and the appreciation of the sea level will be about twice the increase of the tidal range. The tidal wave still spreads from the outer bay to the head bay. Hence, the tidal range will increase towards the left direction of the tidal wave propagation.
section_name: Effects of Tidal Asymmetry

section_content: LeBlond [44] suggested that the shallow estuaries reach a dominant momentum balance between the pressure gradient force and the bottom friction force, so the importance of the two parameters decides the tidal asymmetry in an estuary. We need to explore the ratio of the tidal amplitude to the mean water depth and the ratio of the intertidal storage volume to that of the main river channel [45, 46]. The main effect of the coastline comes from the west to east by the semidiurnal tidal component (e. g., M 2 ) (Table 4 ). Thus, the amplitude and phase of the M 2 tidal component are analyzed (Figure 7 ). In the different SLR scenarios, SLR causes the amplitude of the M 2 tide variation. The amplitude of M 2 in the head bay was concentrated by 0. 1~0. 4. The M 2 amplitude significantly decreased from GP to LCG with an amplitude variation of 1. 2, especially in the RSLR scenario. For the other SLR scenarios, the amplitude variations were about 1. 6~1. 8. In the middle of the inner bay, it reached the maximum value, and the amplitude gradually decreased in the outer sea downstream. Comparisons of their spatial amplitude changes in different SLR scenarios (Figure 7 ) showed the amplitude growth increase with the SLR scale; however, the phase decreased by half with the increase of the SLR scale from the inner bay to the outer bay.
section_name: Tidal Components Behavior

section_content: In this section, SLR enlarges the tidal range and accelerates the wave propa resulting in an earlier tidal arrival. The influence of SLR on the tidal range in Han Tidal asymmetry has traditionally been quantified by using the relative phase between constituents (e. g., 2ϕM 2 -ϕM 4 ) to indicate the direction of the tidal asymmetry and using the ratio of the constituent amplitudes (e. g., aM 2 /aM 4 ) to reflect the degree of distortion [46]. Therefore, the amplitude ratio can be briefly described as for the analysis and conclusions. The M 4 constituent plays a crucial role in determining the tidal asymmetry. Davies and Lawrence [47] indicated that a nontidal M 4 signal could generate near the surface between the interaction of steady wind stress and the M 2 tide. It is often said that the M 2 and M 4 constituents are the main tide and the main shallow water tide, respectively. The ratio variation of the tidal is analyzed, and the inner bay of the M 2 amplitude gets larger. Firstly, the variations of the tidal constituent vary 1. 2 m-2. 4 m from the LCG to ZH section, and the enlargement value of the M 2 constituent is 0. 6m near the outer bay. Secondly, for the different SLR scenarios, the M 2 amplitude gradually decreases from the head bay to the outer bay, and the transmission of tidal waves speeds up from downstream to the GP-CX section. Finally, because the M 4 constituent is generated by the propagation to the shallow sea, the increase of the M 4 amplitude positively correlates with the increase of the M 2 amplitude. Tidal asymmetry has traditionally been quantified by using the relative phase between constituents (e. g., 2φM2 -φM4) to indicate the direction of the tidal asymmetry and using the ratio of the constituent amplitudes (e. g., aM2/aM4) to reflect the degree of distortion [46]. Therefore, the amplitude ratio can be briefly described as for the analysis and conclusions. The M4 constituent plays a crucial role in determining the tidal asymmetry. Davies and Lawrence [47] indicated that a nontidal M4 signal could generate near the surface between the interaction of steady wind stress and the M2 tide. It is often said that the and M4 constituents are the main tide and the main shallow water tide, respectively. The ratio variation of the tidal is analyzed, and the inner bay of the M2 amplitude gets larger. Firstly, the variations of the tidal constituent vary 1. 2 m-2. 4 m from the LCG to ZH section, and the enlargement value of the M2 constituent is 0. 6m near the outer bay. Secondly, for the different SLR scenarios, the M2 amplitude gradually decreases from the head bay to the outer bay, and the transmission of tidal waves speeds up from downstream to the GP-CX section. Finally, because the M4 constituent is generated by the propagation to the shallow sea, the increase of the M4 amplitude positively correlates with the increase of the M2 amplitude. As for the original sea level (Figure 8a ), the ratio of the amplitude continuously decreases about 0. 023 from YG to GP and reaches 0. 01~0. 02 at the ZP-JSZ section. However, it increases to 0. 02~0. 04 in the most outer bay. The ratio of the north bank is smaller than that of the south bank. In the 0. 44 m SLR scenario, the ratio decreases by 0. 04 downstream of YG. The ratio of the sea level keeps pace with the original sea level 20km away from GP-CX. The ratio of most seas appears as the minimal value, which occurs in the inner and outer bay at 0. 01~0. 02. Yet, it increases by 0. 01 (Figure 8b ) in the outer bay. In the 0. 56 m SLR scenario, the ratio from YG to GP is basically consistent with Figure 8b. The ratio of the inner and outer bay is at 0. 01~0. 02 and reaches the minimal value. The ratio of the tidal waves passed from the open sea to the LCG-ZH section is varied by 0. 03~0. 04. The tidal waves propagation increases upstream (Figure 8c ). Most of the SLR scenarios basically keep up with the 0. 56 m SLR scenario, but the ratio decreases by about 0. 02 near YG. The tidal wave obviously propagates from the outer bay to the top bay. The propagation interval shifts downstream (Figure 8d ). For the 0. 77 m SLR scenario (Figure 8e ), the tidal wave passes through the head bay to YG. Its ratio integrally increases, and its value is at 0. 01~0. 02. The lower value in SZP appears the minimum at about 0. 005. The results show that the amplitude ratios of the constituents are basically equal. It increases with the increase of the sea level for most of the sea. The propagation range of tidal waves increases, and the change of its value gradually flattens. The tidal waves spread from the outer bay to the head bay in the original sea level, the variation range is much larger than other sea levels, and the ratio near YG is as high as about 1. 0. In addition, an adequate foundation combines the terrain characteristics of Hangzhou Bay. The average water depth is about 11m. Shoreward propagating surface gravity waves evolve substantially in shallow water owing to refraction, shoaling, nonlinear interactions and dissipation [48]. The nonlinear dissipation of the seabed friction and tidal waves is weak, so the tidal energy loss in the M 2 tidal current is small, and M 4 tidal current changes slowly. However, the water depth is small in the central section, and the ratio begins to gradually decrease. On the other hand, the distribution of the ratio is different at both sides of Hangzhou Bay. As the Coriolis force plays an important role in the ratio, the influence of the terrain becomes shallower from north to south. The tidal wave acts to the right, and the incident tidal wave is stronger on the north bank and so is the M 2 tidal component energy priors to diminution. However, the M 2 constituent gets weak due to the influence of reflected tidal waves on the south bank. Finally, the result considering the impact of the portion makes the ratio of the north bank significantly greater than that of the south bank. 0. 01~0. 02. The lower value in SZP appears the minimum at about 0. 005. The results show that the amplitude ratios of the constituents are basically equal. It increases with the increase of the sea level for most of the sea. The propagation range of tidal waves increases, and the change of its value gradually flattens. The tidal waves spread from the outer bay to the head bay in the original sea level, the variation range is much larger than other sea levels, and the ratio near YG is as high as about 1. 0. In addition, an adequate foundation combines the terrain characteristics of Hangzhou Bay. The average water depth is about 11m. Shoreward propagating surface gravity waves evolve substantially in shallow water owing to refraction, shoaling, nonlinear interactions and dissipation [48]. The nonlinear dissipation of the seabed friction and tidal waves is weak, so the tidal energy loss in the M2 tidal current is small, and M4 tidal current changes slowly. However, the water depth is small in the central section, and the ratio begins to gradually decrease. On the other hand, the distribution of the ratio is different at both sides of Hangzhou Bay. As the Coriolis force plays an important role in the ratio, the influence of the terrain becomes shallower from north to south. The tidal wave acts to the right, and the incident tidal wave is stronger on the north bank and so is the M2 tidal component energy priors to diminution. However, the M2 constituent gets weak due to the influence of reflected tidal waves on the south bank. Finally, the result considering the impact of the portion makes the ratio of the north bank significantly greater than that of the south bank.
section_name: Effects of Tidal Asymmetry

section_content: Previous studies have linked the tidal asymmetry of the sea level rise with the water depth, tidal amplitude, tidal channel and intertidal systems [49, 50]. This paper simulated a sea level rise of 0. 44 m~0. 77 m in 2100 and analyzed the tidal asymmetry response to the SLR. Then, we studied the relationship between the sea level value and tidal range, as well as the formation mechanism and characteristics. It can highly predict the tidal changes and influence in Hangzhou Bay. 

Based on the results of this study, we drew the following conclusions:

(1) The main reason for the unequal duration of the rising and falling tides for the variations of the M 2 constituent is that the flood dominance gradually decreases from west to east, and the skewness gradually increases from the outer bay to the head bay. (2) The M 2 tidal component plays an important role on the maximum tidal range in Hangzhou Bay. The tidal range increases with the rise of the sea level, but the tidal wave propagates from the outer bay to the head bay. Subsequently, the appreciation of the sea level value will be about twice the increase of the tidal range. The tidal range will increase toward the left direction of the tidal wave propagation and accelerate the propagation speed of the tidal waves. (3) The ratio of the tidal amplitude plays a crucial role in determining the asymmetry. 

In the different SLR scenarios, the variations of the M 2 amplitude will significantly increase from the inner bay to the outer bay while the changes of the M 2 phase will be reduced by half from the inner bay to the outer bay. The SLR accelerates the propagation speed of the tidal waves, which will lead to the advance and increase the amplitude of M 2. However, the amplitude of M 4 decreases in the head bay. In the original inner and outer bays, the ratio started to increase from the surroundings to the middle area and gradually increased 0. 01~0. 02, particularly near the SZP area. In the end, the ratio firstly increases, then decreases and, finally, increases from west to east in relation to the topography, seabed friction and nonlinear dissipation of tidal waves. The main influence of the Coriolis force had an impact on the tide wave act toward the right, leading to the ratio of the north bank significantly greater than the south bank. 

Our results showed the tidal asymmetry response to SLR in Hangzhou Bay, and a natural shallow basin provided a reliable numerical simulation. Eventually, a depth understanding the actual variety of the tidal response to the sea level rise in a real system has been proposed and thus provides insight into the impacts of the sea level rise, especially on the tidal asymmetry.
section_name: Conclusions

section_content: This study was financially supported by the National Natural Science Foundation of China, grant number 52101330. Institutional Review Board Statement: Not applicable.

section_name: Data Availability Statement:

section_content: The data are contained within article. 

Acknowledgments: Thank you very much to all the teachers and students who revised the study and put forward their opinions. We are sincerely grateful to the reviewers for their constructive comments. All authors have read and agreed to the published version of the manuscript.
section_name: Data Availability Statement:

section_content: The authors declare no conflict of interest.
section_name: Conflicts of Interest: