收藏切换
On the role of wave breaking in ocean dynamics under typhoon Matsa in the Bohai Sea, China
收藏切换
PDF
Menghan Wang1, , Zengan Deng1, 2, *,
Acta Oceanologica Sinica | 2022, 41(9) : 1 - 18
Less
收藏切换
Acta Oceanologica Sinica | 2022, 41(9): 1-18
Physical Oceanography, Marine Meteorology and Marine Physics
On the role of wave breaking in ocean dynamics under typhoon Matsa in the Bohai Sea, China
Full
Menghan Wang1, , Zengan Deng1, 2, *,
Affiliations
  • 1 School of Marine Science and Technology, Tianjin University, Tianjin 300072, China
  • 2 Guangxi Key Laboratory of Marine Environment Change and Disaster in Beibu Gulf, Qinzhou 535011, China
Published: 2022-09-25 doi: 10.1007/s13131-022-1995-3
Outline
收藏切换

The role of wave breaking (WB) in the ocean dynamics in the Bohai Sea, China under typhoon condition is systematically investigated utilizing a coupled wave-current model. The influences of WB on ocean dynamics and processes (mixing coefficient, temperature, mixed layer depth, and current) during the entire typhoon period (including the pre-typhoon, during-typhoon and after-typhoon stages) are comprehensively detected and discussed. Experimental results show that WB greatly enhances the turbulent mixing at about top 10 m depth under typhoon condition, the increase can be up to 10 times that of the normal weather. At the same time, WB generally strengthens the sea surface cooling by ~1.2°C at the during-typhoon stage, about 3 times that in normal weather. The mixed layer depth, is rapidly increased by ~1.6–3.6 m during typhoon due to WB, particularly, the deepening is stronger in the region from 120.5°E to 121.0°E on account of close to the typhoon eye. In addition, WB renders the current speed more uniformly within the entire depth in the Bohai Sea, the change in speed is ~0.2 m/s, whereas the alternation in current vector is generally opposite to the wind direction except for the typhoon eye region, reflecting that WB has an inhibitory effect on the typhoon-forced current change. The effects of WB on vertical mixing coefficient response to the typhoon rapidly, while the impacts of WB on temperature, and mixed layer depth present hysteretic responses to typhoon. Finally, the mechanisms and distribution characteristics of WB-induced mixing and tidal mixing are compared under typhoon condition.

wave breaking  /  tidal mixing  /  turbulent mixing  /  typhoon  /  coupled model
Menghan Wang, Zengan Deng. On the role of wave breaking in ocean dynamics under typhoon Matsa in the Bohai Sea, China[J]. Acta Oceanologica Sinica, 2022 , 41 (9) : 1 -18 . DOI: 10.1007/s13131-022-1995-3
Both winds and tides are the major mechanical energy driving the turbulent mixing, of which winds are mainly acting on the sea surface (Wunsch and Ferrari, 2004), while the mixing near the bottom is predominantly controlled by tides (Meng et al., 2020). In the upper ocean, in addition to turbulence generated by wind-driven currents, wave breaking (WB) can also produce an enormous amount of turbulent kinetic energy (TKE) injecting downward into the ocean. The momentum of the wind is first transferred to the waves and then to the currents mainly through the action of WB (Phillips, 1977). This process significantly enhances the TKE near the surface (Drennan et al., 1992; Yang et al., 2003). About half of WB-generated TKE dissipates within the depth of significant wave height (Hs), and the remaining continues to transmit downward, which will affect the vertical distribution of the upper ocean temperature and circulation (Feddersen, 2012). The overall effect depth can reach several times of the Hs (Agrawal et al., 1992; Scully et al., 2016; Sun et al., 2006; Gerbi et al., 2009). The observations of the wave breaking zones that are mainly controlled by the tides and winds demonstrate that the turbulence generated by WB is critically important to estimate the vertical structure of the water column in shallow water areas (Zippel et al., 2018). On account of the shallow water depth and the complex terrain in the Bohai Sea (BS), both WB and tidal mixing play significant roles in BS dynamic processes (Jones and Monismith., 2008; Zhao et al., 2019). Therefore, the characteristics and distribution patterns of WB induced mixing and tidal mixing are worthy to be explored in the BS.
Typhoon, as an extreme weather event, may cause strong responses in the upper ocean, for example, forcing intense waves with Hs up to 15 m and reducing the sea surface temperature (SST) in a few degrees (Ginis, 2002). BS locates on the west side of the Pacific Ocean, and is often affected by storm surges and occasionally hit by tropical storms. WB, as an important process in the ocean dynamics, has vital influences on upper ocean turbulence, temperature structure, and boundary layer thickness, etc. Therefore, understanding the role of WB under typhoons (tropical storms) is meaningful to accurately simulate BS dynamic processes, and to improve the disaster forecasting models.
The present WB parameterizations are mainly based on the turbulence closure schemes, many are on the basis of the methods proposed by Mellor and Yamada (1974). Craig and Banner (1994) proposed a WB parameterization method based on the Mellor-Yamada 2.5 (M-Y2.5) turbulence closure scheme. In this method, the wave energy factor αCB is introduced into the boundary condition of the turbulent energy equation to represent the TKE injected by WB. Considering the growth of waves, they found that when αCB is 100 the calculated turbulent dissipation is consistent with the observation. Terray et al. (1999) considered that αCB is a function of waveage, and the typical value is ~150. Besides, Kraus and Turner (1967), Denman and Miyake (1973), Gaspar et al. (1988) also put forward different values of αCB according to their respective researches, and some of them have enormous differences.
Taking account of WB, the sea surface roughness will increase markedly, and the sea surface turbulent microscale will be further influenced. At present, the generally accepted determination method of sea surface roughness is referred to the Charnock relationship in the atmosphere (Donelan, 1990; Smith et al., 1992; Toba et al., 2001; Janssen, 2001), and the Charnock number β is introduced into WB parameterization. Indicated by experiments and observations, Bye (1988) suggested that the value of β is ~1 400 in the ocean, while Stacey (1999) proposed β to be O (105) utilizing the Charnock relationship (Charnock. 1955). Mellor and Blumberg (2004) parameterized β to 2 × 105 in numerical simulations. Gemmrich and Farmer (1999) and Drennan et al. (2005) also suggested different values for β. Zhang et al. (2011) utilized M-Y2.5 scheme together with WB parameterization to discuss the sensitivity of SST to β. The study suggested an optimal combination that αCB is about 167 and β is about 4.1 × 105, based on which the daily and monthly changes in SST were more accurately reflected.
Empirical constants are usually used to estimate αCB and β in the present WB parameterization. However, under extremely high wind speeds, empirical αCB and β obtained from the simple wind-wave relationship are not applicable (Zhang et al., 2011). Compared with using empirical constants, applying the wave variables from the coupled wave-current model can reflect the influence of wave-induced turbulence more precisely. Therefore, this study adopts the Simulating Waves Nearshore-Princeton Ocean Model (SWAN-POM) one-way coupled model to simulate the dynamic processes under typhoon event (Matsa), and conducts diagnostic experiments to analyze the mechanism of WB in BS during Matsa. Finally, the results of this study are compared with the research on tidal mixing (Deng and Zhao, 2020) to further discuss the characteristics of WB and tidal mixing under typhoon condition, and explore their influences on BS dynamic processes.
Princeton Ocean Model (POM) is a baroclinic primitive equation model established by Blumberg and Mellor (1987). In this study, the three-dimensional POM with the M-Y2.5 turbulence closure scheme is used. The turbulent energy equation in the turbulence closure scheme of Mellor (2001) (M01),
$ \begin{split}\frac{{\rm{D}}{q}^{2}}{{\rm{D}}t}=&\frac{\partial }{\partial z}\left({K}_{\mathrm{q}}\frac{\partial {q}^{2}}{\partial z}\right)+2{K}_{\mathrm{M}}\left[{\left(\frac{\partial U}{\partial z}\right)}^{2}+{\left(\frac{\partial V}{\partial z}\right)}^{2}\right]+\\&\frac{2g}{{\rho}_{0}}{K}_{\mathrm{H}}\frac{\partial {\rho }}{\partial z}-2\varepsilon +{F}_{{{\rm{q}}}} \end{split}$
is replaced by turbulent energy equation from Mellor and Blumberg (2004):
$ \frac{\partial {q}^{2}}{\partial t}=\frac{\partial }{\partial z}\left({K}_{\mathrm{q}}\frac{\partial {q}^{2}}{\partial z}\right)+2{K}_{\mathrm{M}}{S}^{2}-2{K}_{\mathrm{H}}{N}^{2}-2\frac{{q}^{3}}{{{B}}_{1}l} , $
where $ {q}^{2}/2 $ is the TKE; $K_{\rm{M}} $ is the vertical viscosity coefficient and $K_{\rm{H}} $, $K_{\rm{q}} $ are the vertical diffusion coefficients; (U, V) are the current velocity components in the directions of (x, y); ρ is the density of seawater and $\rho_0 $ is the mean density of seawater; g is the acceleration of gravity; ε indicates the dissipation and $F_{\rm{q}} $ is the horizontal diffusion term of turbulent kinetic energy; z is the vertical coordinate; t is the time; the average velocity shear and the potential density gradient are contained in S2 and N2 respectively, ${S}^{2}\equiv {\left(\dfrac{\partial U}{\partial z}\right)}^{2}+{\left(\dfrac{\partial V}{\partial z}\right)}^{2}$; $ {N}^{2}\equiv -\dfrac{g}{{\rho }_{0}}\dfrac{\partial \widetilde{\rho }}{\partial z} $ ; ${{B}}_{1}=16.6$ is a constant; l is the turbulent macroscale.
The vertical mixing coefficient of momentum, temperature and turbulence follow the scheme of Mellor and Yamada (1974, 1982) ,
$ \left({K}_{\mathrm{M}},{K}_{\mathrm{H}},{K}_{\mathrm{q}}\right)=lq\left({S}_{\mathrm{M}},{S}_{\mathrm{H}},{S}_{\mathrm{q}}\right) , $
where $ {S}_{\mathrm{M}} $ and $ {S}_{\mathrm{H}} $ are functions of $ {\left(lN/q\right)}^{2} $; generally, $ {S}_{\mathrm{q}} $=0.41$ {S}_{\mathrm{H}} $; $ \left({S}_{\mathrm{M}},{S}_{\mathrm{H}},{S}_{\mathrm{q}}\right) $= (0.30, 0.49, 0.20) near the surface.
The study area covers BS (37.083°–41.033°N, 117.520°–122.470°E), including the Liaodong Bay, the Laizhou Bay, the Bohai Bay, the central sea area, and the Bohai Strait (Fig. 1b). The horizontal resolution is 0.05°×0.05°, and the modeling domain is divided into 80×100 orthogonal rectangular grids. The initial fields, including temperature, salinity, water level, and current, are from the reanalysis dataset of Simple Ocean Data Assimilation (SODA, https://climatedataguide.ucar.edu/climate-data/soda-simple-ocean-data-assimilation). The wind forcing from ERA-interim dataset of the European Centre for Medium-Range Weather Forecasts (ECMWF, https://www.ecmwf.int/en/forecasts/dataset/ecmwf- reanalysis-interim) with resolution of 0.125°×0.125° and interval of 6 h, and the heat flux data from the National Center for Environmental Prediction (NCEP, https://psl.noaa.gov/data/gridded/data.ncep.reanalysis2.html) with an interval of 6 h are adopted and interpolated to model grids. The wind field of typhoon Matsa is also from the ERA-interim of ECMWF, with a resolution of 0.125°×0.125° and an interval of 6 h. The bathymetry is derived from the ETOPO2 dataset of the National Oceanic and Atmospheric Administration (NOAA, https://www.ngdc.noaa. gov/mgg/global/etopo2.html). Twenty-one σ-levels are specified in the vertical and all the outputs are interpolated to the 21 Z-layers (0 m, 1.25 m, 2.5 m, 3.75 m, 5 m, 7.5 m, 10 m, 12.5 m, 15 m, 17.5 m, 20 m, 22.5 m, 25 m, 27.5 m, 30 m, 32.5 m, 35 m, 42.5 m, 50 m, 57.5 m, 65 m). The time step of POM is 200 s in internal mode, and 0.33 s in external mode.
The integration time of the spin-up run is a year, from July 1, 2004 to July 31, 2005. The diagnostic experiments started on August 1, 2005. Four main tidal components (semi-diurnal tide M2 and S2, and diurnal tide K1, O1) are added as the boundary conditions at the model opening boundary. The tidal harmonic constants are provided by the three-dimensional barotropic tide model of the Bohai Sea and the East China Sea (Han et al., 2006). The baroclinic water level at the opening boundary is obtained from the reanalysis dataset of SODA (https://climatedataguide.ucar.edu/climate-data/soda-simple-ocean-data-assimilation).
The Simulating WAves Nearshore (SWAN) is a third-generation wave model, which can predict sea surface gravity waves well in coastal areas such as estuaries, lakes, and coasts (Booij et al., 1999; Zhao and Sun, 2013). In this study, a two-dimensional non-constant mode is used in the SWAN model, and the physical processes include shallow water breaking, bottom friction, wave-wave interaction, and wave refraction. The horizontal resolution of SWAN is 0.05°×0.05°, coherent with POM. The wind forcing is the same as POM. The spectral action volume is discretized into 32 logarithmically spaced frequencies, with the smallest frequency of 0.05 Hz, and the spectral action volume direction is equally divided into 36 directions with an interval of 10°. The JONSWAP spectrum is specified as the initial and the boundary conditions. The time step of the SWAN model is 200 s, which is consistent with the internal mode time step of POM, and the time interval of the output wave spectrum is 6 h, which is the same as the input time interval of the driving field of POM.
In general turbulent closure scheme (Mellor and Durbin, 1975), the boundary conditions of the turbulent energy equation on the atmosphere-ocean interface are the same as those on the atmosphere-solid interface (Mellor and Yamada, 1982), such as
$ {q}^{2}={{B}}_{1}^{2/3}{u}_{*}^{2} , $
where $u_*^{ } $ is the friction velocity.
However, the motion and breaking of sea surface wave make more complicated physical processes in the upper ocean, so that the boundary condition (4) cannot resolve the influences of these physical processes on the turbulence structure in the mixed layer. For this reason, Eq. (5) considering the impact of WB is utilized as the boundary condition of the turbulent energy Eq. (2) (Mellor and Blumberg, 2004),
${K}_{\mathrm{q}}\frac{\partial {q}^{2}}{\partial z}=2{\alpha }_{\mathrm{C}\mathrm{B}}{u}_{*\mathrm{w}}^{3} ,\;\; z=0 ,$
where $ {u}_{*\mathrm{w}} $ is the water-side friction velocity, αCB is determined according to the empirical formula (Terray et al., 1996; 1997),
$ {\alpha }_{\mathrm{C}\mathrm{B}}=15{(C}_{\mathrm{p}}/{u}_{*\mathrm{a}})\mathrm{e}\mathrm{x}\mathrm{p}{[-(0.04{C}_{\mathrm{p}}/{u}_{*\mathrm{a}}\left)\right]}^{4} , $
where $ {u}_{*\mathrm{a}} $ is the air-side friction velocity; $ {C}_{\mathrm{p}}/{u}_{*\mathrm{a}} $ is the wave age; $ {C}_{\mathrm{p}} $ is the phase speed of waves at dominant frequent, calculated from the wavelength and the peak wave period.
$ {u}_{*\mathrm{w}} $ and $ {u}_{*\mathrm{a}} $ are transformed by Mellor and Blumberg (2004):
$ {u}_{*\mathrm{a}}^2=\left({\rho }_{\mathrm{w}}/{\rho }_{\mathrm{a}}\right){u}_{*\mathrm{w}}^2 , $
where $ {\rho }_{\mathrm{w}} $ and $ {\rho }_{\mathrm{a}} $ are the density of seawater and air, respectively; here takes $ {\rho }_{\mathrm{a}}/{\rho }_{\mathrm{w}}=1/900 $ (Mellor and Blumberg, 2004).
$ {u}_{*\mathrm{a}} $ is related to wind speed and roughness (Terray et al., 1996),
$ {u}_{*\mathrm{a}}={k}{U}_{10}/\mathrm{ln}\left(10/{z}_{0}\right) , $
where $ {U}_{10} $ is the wind speed at 10 m; $ {z}_{0} $ is the air-side roughness, which is determined based on the empirical regression relationship of the meteorological observation data at the WAVES stations (Donelan, 1990),
$ {z}_{0}/{H}_{\mathrm{s}}=1.38\times {10}^{-4}{\left({U}_{10}/{C}_{\mathrm{p}}\right)}^{2.66} . $
By including WB, the turbulent macroscale l near the sea surface adopts (Terray et al., 1999; Mellor and Blumberg, 2004)
$ l=\mathrm{m}\mathrm{a}\mathrm{x}(\kappa {z}_{\mathrm{{\rm{w}}}},{l}_{\mathrm{z}}) , $
where $ \kappa $ is the von Karman Constant (0.41); $ {H}_{\mathrm{s}} $ is the significant wave height; $ {z}_{\mathrm{w}} $ is the water-side roughness following (Donelan, 1990; Smith et al., 1992; Toba et al., 2001; Janssen, 2001)
$ {z}_{\mathrm{w}}=\beta {u}_{*\mathrm{{\rm{w}}}}^{2}/\mathrm{g} , $
where β is Charnock number, and g=9.8 $ \mathrm{m}/{\mathrm{s}}^{2} $; Mellor and Blumberg (2004) summarized the works of Donelan (1990), Smith et al. (1992), Janssen (2001), and suggested
$ \beta \cong 665{\left({C}_{\mathrm{{\rm{p}}}}/{u}_{*\mathrm{{\rm{a}}}}\right)}^{1.5} . $
In summary, this study introduces the source term of TKE generated by WB to the upper boundary condition and includes sea surface water-side roughness in the turbulent mixing length equation to determine the WB effect.
POM and SWAN are one-way coupled in this study. The specific process is shown in Fig. 2. SWAN simulates the wave variables in BS, and then the output wave parameters from SWAN are used to calculate αCB and β required by WB parameterization in POM. It is noted that in order to clearly detect the role of WB in a typhoon case, the feedback effects from currents to waves are not discussed in this work.
We choose Matsa (the 9th typhoon in 2005) to perform a typhoon case study. Matsa was formed at about 900 km east of the Philippines at 1200 UTC on 31 July (Fig. 1a) (Song, 2012; Shao et al., 2015), and stepped into the BS at 0600 UTC on 8 August, with the maximum wind speed of ~20 m/s in the BS. The typhoon eye moved northeast to the central sea area (Figs 3b, c) on the 8th, and the intensity of Matsa was weakened rapidly on the 9th with the maximum wind speed of ~12 m/s (Gao et al., 2009; Liu et al., 2010).
Matsa mainly affected the BS on 8 August and 9 August. We select 7–10 August (4 days) to represent three different stages of the typhoon action. Specifically, 7 August represents the normal weather in summer (pre-typhoon stage), 8 August and 9 August represent extremely weather (during-typhoon stage), and 10 August represents the weather with the disappearing typhoon (after-typhoon stage).
To ensure the reliability of our simulations, some verifications are carried out. Tide harmonic analysis is performed with the observed from the tide gauge stations (Fig. 1b) and simulated water levels at 10 tide gauge stations (Table 1). The average absolute deviations of amplitudes of M2, S2, K1, O1 between observations and simulations are 6.44 cm, 5.53 cm, 5.29 cm, 5.09 cm, and the average absolute deviations of the phase lag are 8.74°, 7.87°, 8.26° and 8.89°. The simulations fit well with the observations (Fig. 4), and the scatters are distributed on both sides of the line (y=x) randomly, proving the simulations are reliable.
To further testify the simulations, the simulated SSTs are compared with the observations from the Aqua satellite (from NASA OceanColor Web, https://oceancolor.gsfc.nasa.gov) (Fig. 5b). It is manifested that the simulated monthly-mean SSTs are consistent with the observations, especially in the Bohai Bay, the Laizhou Bay, the Liaodong Bay, the Bohai Strait, and some deviations exist in the western and northern central sea area. The mean absolute deviation of SST is 0.83°C in BS.
Wave parameters are very important to WB parameterization, so it is necessary to validate them. Here the simulated Hs, peak wave period (TPS) and wave direction (WD) are compared to the ERA-interim reanalysis product from ECMWF (0.125°×0.125°). At the pre-typhoon stage, the deviations are tiny (Figs 5c, d), and there are small increases during the typhoon action. Although the simulations have certain deviations, especially during the typhoon period, there is no obvious systematic error. The average absolute deviations of hourly-mean Hs, TPS and WD (from 7 to 10 August) are 0.086 m, 0.10 s and 9.81°, respectively. The corrsponding correlation coefficients are 0.952 2, 0.983 5, and 0.970 4. The above verifications all prove that the simulated results from our model can generally reflect the dynamic processes in BS, and are basically accurate and reliable.
Both diagnostic experiment WAVE-WB (with WB) and coordinate experiment WAVE-NON (without WB) are conducted in this study. Hourly model outputs from 7 to 10 August are used to analyze the influences of WB on dynamics and processes (wave variables, mixing coefficient, temperature, mixed layer depth, and current). The 38.333°N section crosses the central sea area and the center of the Bohai Strait, which is the widest section in BS. Therefore, profiles of various parameters along 38.333°N are investigated in detail. The deviations that represent the WB effects are generated by (WAVE-WB)−(WAVE-NON).
According to the statistics of waves from 1950 to 2011 (Ren et al., 2016), the average Hs in BS is 0.7–1.1 m (~0.7 m in 2005) in summer, the smallest in all seasons. Hs simulated by both WAVE-WB (Fig. 6) and WAVE-NON agrees with the value given by Ren et al. (2016). The responses of the Hs to typhoon during 7 and 10 August are basically following the intensity of wind forcing.
αCB is an empirical parameter related to WB-induced TKE, with the magnitude of $ O\left({10}^{2}\right) $ (Craig and Banner, 1994; Terray et al., 1996). For mature waves, αCB is ~57, and is ~146 for younger waves (Terray et al., 1996, 1997; Craig and Banner, 1994; Stacey, 1999; Mellor and Blumberg, 2004). The wave energy factor is roughly negatively correlated with the wind speed (Figs 7a, d, g, j). On 9 August, the typhoon eye is located between the Bohai Strait and the Liaodong Bay, and αCB near the typhoon eye is larger than that in other regions, suggesting the TKE injected into the surface by winds is more intensively inside the cyclone, while αCB at the center of the cyclone is abnormally small, demonstrating that the vertical mixing in this area is extremely weak. Generally, αCB is linearly associated with the typhoon action, and our estimation is thought to be rational.
β is a key parameter for estimating the water-side surface roughness $ {z}_{\mathrm{w}} $, with the magnitude of $ O\left({10}^{5}\right) $ (Stacey, 1999; Mellor and Blumberg, 2004). $ {z}_{\mathrm{w}} $ is also roughly negatively correlated with wind speed as shown in Figs 7b, e, h and k. $ {{z}}_{\mathrm{w}} $ is related to β, and is of the same magnitude as Hs. In this study, $ {z}_{\mathrm{w}} $ is estimated to be in the range of 0.2–1 m (Figs 7c, f, i, l), similar to Hs (Fig. 6). It is indicated that the estimated wave parameters for WB parameterization are reasonable.
Under high wind speed condition, there is a significant increase (Fig. 8) in TPS, which is closely related to WB (Polton et al., 2003). TPS increases as the increasing of wind speed, so that the WB effect depth is deepened. The wind speed reaches the highest value on the 8th but TPS peaks on the 9th, denoting that the response of waves to the typhoon has a lag.
The vertical viscosity coefficient $ {K}_{\mathrm{M}} $, which is obtained from Eq. (3), is an index to evaluate the vertical mixing strength. Without WB, $ {K}_{\mathrm{M}} $ at the sea surface (Figs 9a, c, e, g) is basically positively related to wind speed, since wind can directly affect the vertical mixing. Including WB, $ {K}_{\mathrm{M}} $ is further modulated, consequently resulting in an integral increase of turbulent mixing in BS (Figs 9b, d, f, h).
Under normal weather, the WB-induced $ {K}_{\mathrm{M}} $ change is relatively small, ranging from ~0.001 m2/s to ~0.01 m2/s (Fig. 9b). High wind speeds have a visible effect on aggravating the turbulence mixing generated by WB. At the during-typhoon stage, the surface $ {K}_{\mathrm{M}} $ change rapidly raises to a level of 0.01–0.15 m2/s, approximately 10 times of that at pre-typhoon stage (Fig. 9d). When wind speed is reduced to below 10 m/s, the WB-induced mixing is decreased to an order of 0.01 m2/s (Fig. 9f), and an abnormally small $ {K}_{\mathrm{M}} $ (close to 0 m2/s) appears at the typhoon eye, implying that the vertical mixing in the vicinity of typhoon center is particularly slight. The relatively large $ {K}_{\mathrm{M}} $ is concentrated in the areas with high wind speeds (the Laizhou Bay and the southern Bohai Strait). It is indicated that the influence of WB on the turbulent mixing is highly dependent on the wind speeds. Although the wind speeds in most area of BS are less than 8 m/s at after-typhoon stage (Fig. 3d), WB still notably affects $ {K}_{\mathrm{M}} $, but the trend is toward the normal weather state. At all the three stages, WB significantly enhances the vertical mixing in the region east of 120°E along 38.333°N (Fig. 10).
In addition, the differences in $ {K}_{\mathrm{M}} $ between WAVE-WB and WAVE-NON decrease with the increasing depth (Figs 10c, f, i, l). Forced by Matsa winds, the surface wind energy injection is enhanced, deepening the impact depth of WB from ~5 m to ~10 m. Overall, WB strengthens the turbulent mixing intensity on the upper ocean, and this action responds rapidly to typhoon forcing.
The SSTs simulated by WAVE-WB are smaller than the that by WAVE-NON (Fig. 11). Before the typhoon, the SST cooling induced by WB is relatively slight (~0.4°C), and the impact depth is limited to the range above 20–30 m isobath (Fig. 11a). When the typhoon stepping into the BS, WB is intensified by the increasing winds as well as the deteriorating sea conditions, and the seawater “cold absorption” phenomenon is further enhanced, therefore SST drops rapidly (with maximum of ~1°C) (Fig. 11b). On 9 August, WB causes a wider cooling (Fig. 11c) in the BS. The SST differences in the entire central sea area and the Bohai Strait are added (the maximum exceeds 1.2°C). This phenomenon demonstrates that the continuous effect of WB brings deeper and colder water up to the surface. Although the intensity of typhoon on 9 August is weakened and the $ {K}_{\mathrm{M}} $ reaches its maximum on the 8th, the response of SST gets its maximum on the 9th, evidencing that there is a hysteresis in temperature changes. After the typhoon, the SST differences are still greater than that at pre-typhoon stage. This continuous cooling reveals that WB, under the action of typhoon, not only affects the SST but also impacts the total heat budget in the BS.
In the vertical, a stably vertical temperature stratification exists in the normal weather of summer (Fig. 12a), and WB usually makes SST slightly lower (~0.3°C) and increase the water temperature within ~10 m depth below the sea surface (Fig. 12b). At the beginning of typhoon, since the duration of the typhoon forcing is still short, the vertical temperature structure does not change drastically. The surface layer temperature change is greatest (up to ~1°C) in the region of 120°–122.5°E, and the impact depth is ~15 m (Fig. 12f). One day after the acting of Matsa, the upper seawater is mixed sufficiently, breaking the stratification within 5 m depth, and the impact penetrates deep to ~20 m (Fig. 12j). If the typhoon could last longer, WB would cause a more intensively heat exchange. After the typhoon, the temperature variation caused by WB remains significant, which again reflects the hysteretic effect of temperature in response to the typhoon.
Yablonsky and Ginis (2008) defined mixed layer depth (MLD) as the depth of the water layer at a temperature difference of 0.5°C from the sea surface. Here we calculate the MLD in BS following this definition. In normal summer weather, the vertical temperature stratification in BS is obvious. The temperature differences between the surface and subsurface layers are notable, MLD is about 1–2 m. However, the relatively shallow areas near the shore are sufficiently mixed. As water depth increasing, the MLD in the region east of 119°E along 38.333°N rapidly decreases to about 1 m (Fig. 13a), and is slightly deepened (with a increase of ~0.2 m, but there is an obvious increase of ~1.8 m in 120°–121°E) after adding WB.
By the forcing of Matsa, WB renders a significant deepening effect (~2.4 m) on MLD in the relatively deep water zone especially in 120°–121°E (Fig. 13b). On 9 August, the MLD in east of 121°E region simulated by WAVE-NON is greater than 2 m (Fig. 13c), which denotes that the typhoon can directly deepen the MLD. After adding WB, the MLD rapidly deepens to around 4.5–7.8 m, signaling that the impact depth of WB is about 1.6–3.6 m. Compared with that at the pre-typhoon stage, typhoon markedly intensifies the impact of WB on MLD. Contrasting with the temperature profile, we infer that the water above MLD is sufficiently mixed. When the typhoon is disappeared, the WB impact retreats to the depth of about 0.8–3.5 m, slightly smaller than that at the during-typhoon stage (Fig. 13d).
From the perspective of the horizontal distribution of MLD, WB-induced deepening before the typhoon is limited to the central sea area and is relatively weak (Fig. 14b). On the 8th, the impact is expanded to the central sea area, the Bohai Bay and the Liaodong Bay (Fig. 14d). Since the typhoon eye is close to the Liaodong Bay on the 9th, WB has a noticeable impact on MLD in the Liaodong Bay, the central sea, and the northern Bohai Strait, leading to a significant deepening of the MLD (Fig. 14f). After the typhoon, the isoline values become slightly small (Fig. 14h) and the situation is gradually recovering to the state before typhoon.
Before the typhoon, currents in BS generally flow from the Liaodong Bay, the Bohai Bay and the Bohai Strait to the central sea area and the Laizhou Bay. WB does not change the integral direction of currents, but it alters the local flow patterns (Fig. 15b), particularly in the southern part of Bohai Strait and the central sea area (with a maximum speed change of ~0.07 m/s). In particular in a part of the Bohai Bay and the Bohai Strait, the changes in current vectors (CCVs, obtained by (WAVE-WB)−(WAVE-NON)) are consistent with the wind directions (Fig. 15c).
On the 8th, the current speed in the central sea area is decreased (with maximum change of ~0.2 m/s) by WB. In the most areas of BS, the cosine of the current change and the wind direction is −1 (Fig. 15f), suggesting that the CCVs are opposite to the wind vectors, reflecting an interesting phenomenon that hasn’t been proposed before: WB has an inhibitory effect on the typhoon-forced current change.
One day after Matsa entering BS, WB-induced speed change is about ±0.1 m/s near the cyclone and the Bohai Strait, and it approaches 0 m/s at the center of the typhoon. Because of the influence of winds, the current direction in the Laizhou Bay is gradually parallel to the wind direction (Fig. 15g), implying that during typhoon the wind is the more important factor in affecting currents than tide. In addition, the CCVs are opposite to the wind direction at the periphery of typhoon, and follow the wind direction near the typhoon eye (Fig. 15i). This phenomenon shows that, around the center of the cyclone, the surface current is more intensively affected by winds, and the CCVs follow the wind direction. While the surface current differences at the periphery of typhoon are slight; WB has an inhibitory effect on wind action. After the typhoon, the current change in the cyclone area follows the direction of wind, spreading from the cyclone center to the surrounding area, and the distribution of CCVs (Fig. 15l) return to in line with that at the pre-typhoon stage.
It is worth noting that the areas with obviously current changes are concentrated in the central sea area and the Bohai Strait, basically overlapped with the areas of marked temperature changes, which suggests that WB exerts significant effects on the turbulent mixing in these areas (Figs 15b, e, h, k and 11).
At the pre-typhoon stage, the impact depth of WB is relatively shallow (~7 m) (Fig. 12d). After Matsa entering BS, the current velocities within the depth of –10 m are more evenly distributed under the action of WB, the resulting changes in speed are ±0.15 m/s (Fig. 12h). However, due to the short duration of typhoon action at this stage, TKE is still limited to a few meters below the sea surface, thus the stratification of currents can only be formed in the area with relatively deep water (east of 120°E). Consequently, the velocity differences between the surface and bottom layers are enlarged. But the deviations are decreased obviously after the typhoon lasting for a day due to the injected TKE can be fully transmitted to the depth of ~14 m. As the disappearance of typhoon, the current profile changes from the surface downwards, forming a temporary current stratification (Fig. 12p). WB promotes the synergistic changes of currents in BS. Thus, to a certain extent, WB weakens the current stratification.
Deng and Zhao (2020) performed a case study of the typhoon Matsa to analyze the impacts of tidal mixing on ocean dynamics in BS. To reveal the different influences between WB and tidal mixing on BS dynamics under typhoon condition, here we conduct a comprehensive comparison (Table 2). From a regional point of view, during Matsa, the tidal mixing coefficient is generally increased in the Bohai Bay, the Laizhou Bay, the Liaodong Bay, and the Bohai Strait, and the effect of typhoon on the tidal mixing in the central sea area is relatively small. Whereas, WB-induced mixing during typhoon increases significantly in the central sea area, the Laizhou Bay, and the Bohai Strait. Both WB-induced mixing and tidal mixing can be intensified by approximately 1 order of magnitude by typhoon in BS. The mixing caused by WB under typhoon is generally larger than tidal mixing by 1 order, especially in the central sea area, it is about 3 orders (Fig. 9d). During the typhoon, differences between the two mixings are mainly presented in the central sea area. Tidal mixing is greatly affected by the topography, so that in the coastal area (three bay areas) the enhancement is more obviously; while in the offshore deeper water region (the central sea area and the Bohai Strait), the TKE generated by WB makes the mixed layer thicker.
By the continuous forcing of typhoon, the effects of both tidal mixing and WB on temperature are intensified, and the persistent actions of the two effects facilitate the mixing of surface water with deep colder water, resulting in a significant cooling in the central sea area. Although the tidal mixing in the bay areas is very strong, the tidal mixing-related temperature changes are not that noticeable due to the weak vertical stratification and the overall warm water column. The impact of WB on the temperature is the strongest in the central sea area and the Bohai Strait (Fig. 12). In combination of the profile diagram (Fig. 12f) and temperature differences plots at depths of 0 m, 5 m, 15 m, 25 m (Fig. 16), it can be seen that the WB-induced change at 25 m is less than 0.01°C, a trivial effect, demonstrating that at the during-typhoon stage the impact depth of WB (~15 m) is smaller than that of tidal mixing, and the WB’s influence on temperature (±1°C) is weaker than that of tidal mixing.
The distribution pattern of currents under typhoon is directly related to wind forcing. Tidal mixing has trivial effect on the structure and magnitude of circulation. However, WB exerts an impact on local current speed and direction, the maximum speed change around the central sea area is about −0.2 m/s. In BS, the cosine values of CCVs and wind directions are almost −1 (Fig. 15), revealing that WB induces an inhibitory effect on the typhoon impact, and this phenomenon can not be found in tidal mixing.
Deng and Zhao (2020) investigated the tidal mixing-resulted water level changes at four stations respectively located in the Bohai Bay, the Laizhou Bay, the Liaodong Bay, and the central sea area near the Liaodong Peninsula, and revealed that tidal mixing would partly compensate for wind forcing and cause alternation in tidal phase. We also investigate the effect of WB on water level under typhoon Matsa, and find that the water level change caused by WB on the 8th is O(10−3) (Fig. 17a), and the alternations are within 0.01 m at the four spots from the 8th to 10th (Fig. 17b), which suggests that the impact of WB on the water level is tiny, and there is almost no associated phase change.
In this study, a one-way coupled wave-current model is used to investigate the influences of WB on ocean dynamics under typhoon weather in the BS. Before the diagnostic experiments, some model verifications (including tidal harmonic constant, SST, Hs, TPS, and WD) are carried out, proving that our model is reliable and can generally reproduce the dynamic processes in BS.
The WB-induced turbulent mixing is strongest at the surface layer and gradually weakens as the increasing depth. During the typhoon, $ {K}_{\mathrm{M}} $ is increased apparently, with the largest one of ~0.15 m2/s, reaching ~10 times of that at pre-typhoon stage. Particularly, the WB impact on $ {K}_{\mathrm{M}} $ is close to 0 m2/s at the center of the typhoon. Due to the increase of wind forcing and the deterioration of sea condition, WB brings deeper and colder water up to the surface, consequently increasing the temperature cooling extent (from ~0.4°C to ~1.2°C) and range. In addition, under the continuous acting of the typhoon, WB breaks the vertical stratification at the 5 m water layer, and the impact range deepens from ~10 m to ~20 m.
At the during-typhoon stage, WB deepens the mixed layer for about 1.6–3.6 m, and the increase of MLD is largest near the typhoon eye. The influence range also expands from the central sea area to the central sea area, the Liaodong Bay, and the Bohai Strait. WB doesn’t change the integral current direction in BS, but it impacts the local flow patterns. Under the action of typhoon, the CCVs caused by WB are almost opposite to the wind direction. The impact depth of WB on current is shallow, with a range of ~7 m in normal weather and ~14 m in typhoon condition.
The areas with maximum change in current basically coincide to the regions with maximum temperature change during the entire typhoon period, inferring that WB has the greatest impact on mixing in those areas. Diagnostic experiments evidence that the influences of WB on temperature are cumulative and have an obviously lag.
Comparing with the research of Deng and Zhao (2020) on tidal mixing, we find that tidal mixing has greater influence on temperature structure, with larger impact range and higher intensity, and it will compensate for wind forcing and cause tidal phase alternation. Whereas, WB has a relatively larger impact on currents and generally inhibits the typhoon forcing. Besides, the tidal mixing effect in the three bay areas of BS is relatively strong, and is generally weak in the central sea area and the Bohai Strait, while the situations of WB are opposite.
This study discusses the effects of WB on dynamics in BS under typhoon conditions. Non-breaking waves such as Langmuir turbulence (LT) also have a significant effect on dynamics. In our previous work, we investigated the effects of LT under common weather (Cao et al., 2019). By a briefly comparison of WB and LT effects, we find that LT can impact the entire mixed layer, while WB only acts within the top several meters. The LT-caused change in $ {K}_{\mathrm{M}} $ increases with depth, reaching a maximum (~0.001 m2/s) in the 35 m layer. The effect of WB on $ {K}_{\mathrm{M}} $ at the sea surface (~0.01 m2/s) is greater than that of LT, but this effect decays rapidly with depth, acting within ~5 m below the sea surface. The temperature change caused by LT in BS is greater than 0.2°C even at the 15 m layer, which is greater than that by WB. The impact of LT on current velocity (~0.1 m/s) is also slightly greater than that of WB (~0.07 m/s). In general, both LT and WB make the vertical current velocity distribution more uniformly.
This research is mainly focused on the impacts of WB on BS dynamics in summer typhoon weather. However, the mechanism and distribution characteristics of WB under common weather in spring and autumn are also worth exploring, which will help us further understand the role of WB in BS. The simulation results of this study still have certain deviations from the realistic ocean state, so that the model would be further improved in the future. Limited by experimental conditions, this study only uses a one-way coupled model. In the future, we will detect the impact of WB using an atmosphere-circulation-wave fully coupled model, and further conduct a systematic study on the mechanisms of WB, tidal mixing, and LT in BS.
  • The Grant from Guangxi Key Laboratory of Marine Environment Change and Disaster in Beibu Gulf under contract No. 2021KF03; the National Natural Science Foundation of China under contract Nos 42176020 and 42076007; the Foundation from Key Laboratory of Marine Environmental Information Technology, Ministry of Natural Resources of China under contract No. 2020GKF-0812; the Tianjin Natural Science Foundation under contract No. 18JCYBJC84900.
Agrawal Y C, Terray E A, Donelan M A, et al. 1992. Enhanced dissipation of kinetic energy beneath surface waves. Nature, 359(6392): 219–220, doi: 10.1038/359219a0
Blumberg A F, Mellor G L. 1987. A description of a three-dimensional coastal ocean circulation model. In: Heaps N S, ed. Three-Dimensional Coastal Ocean Models. Washington: American Geophysical Union, 1–16, doi: 10.1029/CO004p0001
Booij N, Ris R C, Holthuijsen L H. 1999. A third-generation wave model for coastal regions: 1. Model description and validation. Journal of Geophysical Research, 104(C4): 7649–7666, doi: 10.1029/98JC02622
Bye J A T. 1988. The coupling of wave drift and wind velocity profiles. Journal of Marine Research, 46(3): 457–472, doi: 10.1357/002224088785113559
Cao Yu, Deng Zengan, Wang Chenxu. 2019. Impacts of surface gravity waves on summer ocean dynamics in Bohai Sea. Estuarine, Coastal and Shelf Science, 230: 106443,
Charnock H. 1955. Wind stress on a water surface. Quarterly Journal of the Royal Meteorological Society, 81(350): 639–640, doi: 10.1002/qj.49708135027
Craig P D, Banner M L. 1994. Modeling wave-enhanced turbulence in the ocean surface layer. Journal of Physical Oceanography, 24(12): 2546–2559, doi: 10.1175/1520-0485(1994)024<2546:MWETIT>2.0.CO;2
Deng Zengan, Zhao Yifei. 2020. Impact of tidal mixing on water mass properties and circulation in the Bohai Sea: a typhoon case. Journal of Marine Systems, 206: 103338, doi: 10.1016/j.jmarsys.2020.103338
Denman K L, Miyake M. 1973. Upper layer modification at ocean station Papa: observations and simulation. Journal of Physical Oceanography, 3(2): 185–196, doi: 10.1175/1520-0485(1973)003<0185:ULMAOS>2.0.CO;2
Donelan M A. 1990. Air-sea interaction. In: Le Mehaute B, Hanes D M, eds. The Sea: Ocean Engineering Science. New York: Wiley-Interscience Press, 239–292
Drennan W M, Kahma K K, Terray E A, et al. 1992. Observations of the enhancement of kinetic energy dissipation beneath breaking wind waves. In: Banner M L, Grimshaw R H J, eds. Breaking Waves: IUTAM Symposium Sydney, Australia 1991. Berlin, Heidelberg: Springer, 95–101
Drennan W M, Taylor P K, Yelland M J. 2005. Parameterizing the sea surface roughness. Journal of Physical Oceanography, 35(5): 835–848, doi: 10.1175/JPO2704.1
Feddersen F. 2012. Observations of the surf-zone turbulent dissipation rate. Journal of Physical Oceanography, 42(3): 386–399, doi: 10.1175/JPO-D-11-082.1
Gao Shuanzhu, Meng Zhiyong, Yang Guiming. 2009. Study on the predictability of the recurvature of Typhoon Matsa (0509) in Bohai. Meteorological Monthly, 35(2): 8–14
Gaspar P. 1988. Modeling the seasonal cycle of the upper ocean. Journal of Physical Oceanography, 18(2): 161–180, doi: 10.1175/1520-0485(1988)018<0161:MTSCOT>2.0.CO;2
Gemmrich J R, Farmer D M. 1999. Near-surface turbulence and thermal structure in a wind-driven sea. Journal of Physical Oceanography, 29(9): 480–499, doi: 10.1175/1520-0485(1999)029<0480:NSTATS>2.0.CO;2
Gerbi G P, Trowbridge J H, Terray E A, et al. 2009. Observations of turbulence in the ocean surface boundary layer: energetics and transport. Journal of Physical Oceanography, 39(5): 1077–1096, doi: 10.1175/2008JPO4044.1
Ginis I. 2002. Tropical cyclone-ocean interactions. In: Perrie W, ed. Atmosphere-Ocean Interactions. Southampton: WIT Press
Han Guijun, Li Wei, He Zhongjie, et al. 2006. Assimilated tidal results of tide gauge and TOPEX/POSEIDON data over the China seas using a variational adjoint approach with an onlinear numerical model. Advances in Atmospheric Sciences, 23(3): 449–460
Janssen P A E M. 2001. Reply. Journal of Physical Oceanography, 31(8): 2537–2544, doi: 10.1175/1520-0485(2001)031<2537:R>2.0.CO;2
Jones N L, Monismith S G. 2008. Modeling the influence of wave-enhanced turbulence in a shallow tide- and wind-driven water column. Journal of Geophysical Research, 113(C3): C03009, doi: 10.1029/2007JC004246
Kraus E B, Turner J S. 1967. A one-dimensional model of the seasonal thermocline II. The general theory and its consequences. Tellus, 19(1): 98–106, doi: 10.3402/tellusa.v19i1.9753
Liu Jiajun, Zhang Lifeng, Guan Jiping, et al. 2010. Application of ensemble prediction products to the forecast of Typhoon Masta. Meteorological Monthly, 36(5): 21–31
Mellor G L. 2001. One-dimensional, ocean surface layer modeling: a problem and a solution. Journal of Physical Oceanography, 31(3): 790–809, doi: 10.1175/1520-0485(2001)031<0790:ODOSLM>2.0.CO;2
Mellor G, Blumberg A. 2004. Wave breaking and ocean surface layer thermal response. Journal of Physical Oceanography, 34(3): 693–698, doi: 10.1175/2517.1
Mellor G L, Durbin P A. 1975. The structure and dynamics of the ocean surface mixed layer. Journal of Physical Oceanography, 5(4): 718–728, doi: 10.1175/1520-0485(1975)005<0718:TSADOT>2.0.CO;2
Mellor G L, Yamada T. 1974. A hierarchy of turbulence closure models for planetary boundary layers. Journal of the Atmospheric Sciences, 31(7): 1791–1806, doi: 10.1175/1520-0469(1974)031<1791:AHOTCM>2.0.CO;2
Mellor G L, Yamada T. 1982. Development of a turbulence closure model for geophysical fluid problems. Reviews of Geophysics, 20(4): 851–875, doi: 10.1029/RG020i004p00851
Meng Qingjun, Li Peiliang, Zhai Fangguo, et al. 2020. The vertical mixing induced by winds and tides over the Yellow Sea in summer: a numerical study in 2012. Ocean Dynamics, 70(7): 847–861, doi: 10.1007/s10236-020-01368-2
Phillips O M. 1977. Dynamtics of the Upper Ocean. Cambridgeshire: Cambridge University Press, 366
Polton J A, Lewis D M, Belcher S E. 2003. The Role of Wave-Induced Coriolis-Stokes Forcing on the Wind-Driven Mixed Layer. Journal of Physical Oceanography, 35(4): 444–457
Ren Huiru, Li Guosheng, Cui Linlin, et al. 2016. Wave climate changes in Bohai Sea related to the East Asian circulation oscillations over the last sixty years. Climatic and Environmental Research, 21(4): 490–502, doi: 10.3878/j.issn.1006-9585.2016.15281
Scully M E, Trowbridge J H, Fisher A W. 2016. Observations of the transfer of energy and momentum to the oceanic surface boundary layer beneath breaking waves. Journal of Physical Oceanography, 46(6): 1823–1837, doi: 10.1175/JPO-D-15-0165.1
Shao Jinchao, Zhao Hui, Shen Chunyan, et al. 2015. Influence of typhoon Matsa on phytoplankton Chlorophyll-a in the Northwest Pacific Ocean offshore and alongshore. Journal of Guangdong Ocean University, 35(4): 67–74
Smith S D, Anderson R J, Oost W A, et al. 1992. Sea surface wind stress and drag coefficients: the hexos results. Boundary-Layer Meteorology, 60(1−2): 109–142, doi: 10.1007/BF00122064
Song Wenpeng. 2009. The analysis of the structure of T-S and the current characteristics in Baohai Sea during winter and summer (in Chinese)[dissertation]. Qingdao: Ocean University of China
Song Yaoling. 2012. Effect of activity characteristics of Typhoon Matsa on aviation meteorological observation and support. Meteorological, Hydrological and Marine Instruments (in Chinese), 29(1): 99–104,
Stacey M W. 1999. Simulation of the wind-forced near-surface circulation in Knight Inlet: a parameterization of the roughness length. Journal of Physical Oceanography, 29(6): 1363–1367, doi: 10.1175/1520-0485(1999)029<1363:SOTWFN>2.0.CO;2
Sun Qun, Guan Changlong, Song Jinbao. 2006. Effect of wave breaking on turbulent energy budgets in ocean surface mixed layer. Oceanologia et Limnologia Sinica, 37(1): 69–74, doi: 10.3321/j.issn:0029-814X.2006.01.011
Terray E A, Donelan M A, Agrawal Y C, et al. 1996. Estimates of kinetic energy dissipation under breaking waves. Journal of Physical Oceanography, 26(5): 792–807, doi: 10.1175/1520-0485(1996)026<0792:EOKEDU>2.0.CO;2
Terray E A, Donelan M A, Agrawal Y C, et al. 1997. Reply. Journal of Physical Oceanography, 27(10): 2308–2309, doi: 10.1175/1520-0485(1997)027<2308:R>2.0.CO;2
Terray E A, Drennan W M, Donelan M A. 1999. The vertical structure of shear and dissipation in the ocean surface layer. In: Proceedings of Symposium on the Wind-driven Air-Sea Interface—Electromagnetic and Acoustic Sensing, Wave Dynamics, and Turbulent Fluxes. Sydney: University of New South Wales, 239–245
Toba Y, Smith S D, Ebuchi N. 2001. Historical drag expressions. In: Jones I S F, Toba Y, eds. Wind Stress over the Ocean. New York: Cambridge University Press
Wunsch C, Ferrari R. 2004. Vertical mixing, energy, and the general circulation of the oceans. Annual Review of Fluid Mechanics, 36: 281–314, doi: 10.1146/annurev.fluid.36.050802.122121
Yablonsky R M, Ginis I. 2008. Improving the ocean initialization of coupled hurricane–ocean models using feature-based data assimilation. Monthly Weather Review, 136(7): 2592–2607, doi: 10.1175/2007MWR2166.1
Yang Yongzeng, Qiao Fangli, Xia Changyong, et al. 2003. Effect of ocean wave momenfum and mixing on upper ocean. Advances in Marine Science, 21(4): 363–368, doi: 10.3969/j.issn.1671-6647.2003.04.001
Zhang Xuefeng, Han Guijun, Wu Xinrong, et al. 2011. Effect of surface wave breaking on upper-ocean structure revealed by assimilating sea temperature data. Journal of Tropical Oceanography, 30(5): 48–54, doi: 10.11978/j.issn.1009-5470.2011.05.048
Zhao Yifei, Deng Zengan, Yu Ting, et al. 2019. Numerical study on tidal mixing in the Bohai Sea. Marine Geodesy, 42(1): 46–63, doi: 10.1080/01490419.2018.1539055
Zhao Xin, Sun Qun. 2013. Influence of reclamation on hydrodynamic environment in Bohai Bay. Advanced Materials Research, 726–731: 3262–3265,
Zippel S F, Thomson J, Farquharson G. 2018. Turbulence from breaking surface waves at a river mouth. Journal of Physical Oceanography, 48(2): 435–453, doi: 10.1175/JPO-D-17-0122.1
Year 2022 volume 41 Issue 9
PDF
75
40
Cite this Article
BibTeX
Article Info
doi: 10.1007/s13131-022-1995-3
  • Receive Date:2021-07-08
  • Online Date:2025-11-21
  • Published:2022-09-25
Article Data
Affiliations
History
  • Received:2021-07-08
  • Accepted:2021-11-29
Funding
The Grant from Guangxi Key Laboratory of Marine Environment Change and Disaster in Beibu Gulf under contract No. 2021KF03; the National Natural Science Foundation of China under contract Nos 42176020 and 42076007; the Foundation from Key Laboratory of Marine Environmental Information Technology, Ministry of Natural Resources of China under contract No. 2020GKF-0812; the Tianjin Natural Science Foundation under contract No. 18JCYBJC84900.
Affiliations
    1 School of Marine Science and Technology, Tianjin University, Tianjin 300072, China
    2 Guangxi Key Laboratory of Marine Environment Change and Disaster in Beibu Gulf, Qinzhou 535011, China

Corresponding:

References
Share
https://castjournals.cast.org.cn/joweb/aos/EN/10.1007/s13131-022-1995-3
Share to
QR

Scan QR to access full text

Cite this article
BibTeX
Citations
表12种不同金属材料的力学参数

Family
属数
Number of
genus
种数
Number of
species
占总种数比例
Percentage of
total species (%)

Genus
种数
Number of
species
占总种数比例
Percentage of total
species (%)
鹅膏菌科Amanitaceae 2 11 5.26 鹅膏菌属 Amanita 10 4.78
小菇科 Mycenaceae 2 12 5.74 丝盖伞属 Inocybe 5 2.39
多孔菌科 Polyporaceae 8 14 6.70 蜡蘑属 Laccaria 5 2.39
红菇科 Russulaceae 3 23 11.00 小皮伞属 Marasmius 6 2.87
小菇属 Mycena 11 5.26
光柄菇属 Pluteus 5 2.39
红菇属 Russula 17 8.13
栓菌属 Trametes 5 2.39
关闭全屏
  • BibTeX
  • EndNote
  • RefWorks
  • TxT