Next Article in Journal
A Method for Extracting Some Key Terrain Features from Shaded Relief of Digital Terrain Models
Next Article in Special Issue
Storm Surge Hazard Assessment of the Levee of a Rapidly Developing City-Based on LiDAR and Numerical Models
Previous Article in Journal
Motion Parameter Identification and Motion Compensation for Shipborne HFSWR by Using the Reference RF Signal Generated at the Shore
Previous Article in Special Issue
Monitoring Suaeda salsa Spectral Response to Salt Conditions in Coastal Wetlands: A Case Study in Dafeng Elk National Nature Reserve, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Enhanced Wave-Induced Mixing on the Ocean Upper Mixed Layer during Typhoon Nepartak in a Regional Model of the Northwest Pacific Ocean

1
Ocean College, Zhejiang University, Zhoushan 316000, China
2
First Institute of Oceanography, Ministry of Natural Resources, Qingdao 266061, China
3
Laboratory for Regional Oceanography and Numerical Modeling, Pilot National Laboratory for Marine Science and Technology, Qingdao 266071, China
4
Key Laboratory of Marine Science and Numerical Modeling (MASNUM), Ministry of Natural Resources, Qingdao 266061, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(17), 2808; https://doi.org/10.3390/rs12172808
Submission received: 30 July 2020 / Revised: 26 August 2020 / Accepted: 27 August 2020 / Published: 30 August 2020
(This article belongs to the Special Issue Coastal Environments and Coastal Hazards)

Abstract

:
To investigate the effect of wave-induced mixing on the upper ocean structure, especially under typhoon conditions, an ocean-wave coupled model is used in this study. Two physical processes, wave-induced turbulence mixing and wave transport flux residue, are introduced. We select tropical cyclone (TC) Nepartak in the Northwest Pacific ocean as a TC example. The results show that during the TC period, the wave-induced turbulence mixing effectively increases the cooling area and cooling amplitude of the sea surface temperature (SST). The wave transport flux residue plays a positive role in reproducing the distribution of the SST cooling area. From the intercomparisons among experiments, it is also found that the wave-induced turbulence mixing has an important effect on the formation of mixed layer depth (MLD). The simulated maximum MLD is increased to 54 m and is only 1 m less than the observed value. The wave transport flux residue shows a dominant role in the mixed layer temperature (MLT) changing. The mean error of the MLT is reduced by 0.19 °C compared with the control experiment without wave mixing effects. The study shows that the effect of wave mixing should be included in the upper ocean structure modeling.

Graphical Abstract

1. Introduction

Vertical mixing of the ocean is known to be one of the most important processes in regulating the SST and mixed layer structure, which are two key ocean factors that control the exchange of heat and momentum between the ocean and atmosphere. For a model to be able to simulate realistic ocean dynamics, it must include an accurate characterization of the vertical mixing process. However, due to the small-scale turbulence processes involved, vertical mixing usually cannot be explicitly expressed in general ocean circulation models and, so far at least, its implementation must be parameterized.
The vertical mixing in many three-dimensional numerical ocean circulation models is often based on two-equation turbulence closure schemes [1,2,3,4,5,6,7,8]. The turbulent kinetic energy (TKE) equation has a firm physical basis, while the equation for length scale, no matter which form is used, is merely patterned after the TKE equation and consequently, its physical basis is tenuous. Unlike the TKE equation, all terms, like production, destruction, and diffusion term in the length scale equation, have to be modeled, which also is a major difficulty in two-equation turbulence models of this type [9]. A common problem of such schemes is an overestimation of SST and underestimation of mixed layer depth (MLD), especially during the summer [1,10,11,12,13,14,15]. For instance, Martin [15] showed that summertime temperatures at Ocean Weather Station Papa exceeded observations by about 2 °C relative to an annual range of 8 °C. Recently, more attention has been paid to the role of the waves which believed to be the primary reason for the insufficient vertical mixing in the upper ocean. By careful tuning, many studies showed some improvement in the final results [9,13,16,17,18]. However, more numerical modeling investigation is needed, especially under different real ocean conditions. At the same time, some studies are concentrated on breaking waves which can be applied to ocean model by using surface boundary conditions. However, further studies have shown that the wave breaking generated turbulence is mostly confined to the surface layer with their energy dissipated locally at a depth of the order of the wave amplitude [19,20]. Yet, all these efforts gave us a strong hint that wave could not be neglected in the study of ocean dynamics.
Typhoons (or tropical cyclones (TCs)) are probably among the most extreme sea–air interactions. These extremes can drastically change the marine environment and seriously affect human activities. Thus, the ability to simulate changes in the marine environment during TCs is vital. The large waves caused by TCs represent a main feature of the sea surface and exert an important influence on the vertical mixing of the upper ocean. Observations and model experiments show that sea surface cooling caused by the TC can exceed 3 °C over the hundreds of kilometers of the TC center, and cooling of 7–11 °C may occur and could last for 1–3 weeks [21,22,23]. Moreover, a strong TC entrains cold water into the mixed layer (ML) from below, resulting in the cooling of the ML and deepening of the MLD [24]. Wang and Qiao [25] introduced the mixing effect of wave-induced turbulence in the ocean current model, and the results show that ocean waves have a significant influence on the decline of SST and deepening of the MLD during a TC. However, many studies have also shown that even with this added mechanism, the strength of thermocline is still weak [26,27,28,29]. There are many other deficiencies in the simulation results under TC conditions, such as large differences in the range, intensity and spatial structure of sea surface cooling, offset of the strong cooling zone, and abnormal structure of the upper mixed layer and so on. The lack of skill of simulation ability may be due to the insufficient consideration of the physical processes. Based on field observations, Toffoli et al. [30] found that variations of a normally shallow mixed layer depth are observed within a relatively short timescale of approximately 10 h after the influence of TC. The quick response of the mixed layer seems to sustain the role of wave-induced turbulence in the upper ocean mixing. Moreover, using coupled hurricane-ocean-wave modeling, Aijaz et al. [31] found that the effect of wave breaking on SST cooling is minor except the center of the TC track. In addition, the model skill is improved by enhancing the turbulence due to nonbreaking waves.
In this study, a 3D wave-current coupled model consisting of the Princeton Ocean Model (POM) and the Marine Science and Numerical Modeling (MASNUM) wave model is used to investigate to what degree the wave-induced mixing can affect the upper ocean superstructure. Except for the wave-induced turbulence mixing mechanism, another wave-induced mixing mechanism called the wave transport flux residue is also introduced into the ocean circulation model. Preliminary research has shown that the latter can affect the structure of the upper ocean as a driving source. The Northwest Pacific ocean is used as the research area, and TC Nepartak is used as a TC example, and the roles of the two mechanisms in ocean mixing are studied through comparative experiments, especially in the case of TC. A series of comparisons are conducted between the experimental results and the available measured data. This paper is structured as follows: The two physical processes are described in detail in Section 2; the descriptions of material and method used in this study are included in Section 3; the model results are presented in Section 4; the relevant discussions are given in Section 5 and the conclusions are summarized in Section 6.

2. Wave-Induced Parameterization Schemes

2.1. Wave-Induced Turbulence Mixing

Surface wave breaking and nonbreaking waves can enhance diffusivity, particularly vertical diffusivity, in the ocean. The turbulence generated by surface wave breaking decays rapidly with distance from the sea surface and is therefore considered as a minor contributor to the overall ocean mixing [20,32]. Studies [33,34] have shown that nonbreaking wave-induced turbulence, which is distributed vertically through the water column at a scale comparable to that of the wave length [35], is a predominant mechanism for turbulence generation and enhancement and can significantly affect the physics of ocean mixing.
Xia et al. [28], Qiao et al. [36], Dai et al. [37] and others have shown that the ability of the model to simulate ocean mixed layers and SST can be significantly improved via a rational consideration of wave-induced turbulence mixing (hereafter WITM) during ocean modeling. Qiao et al. [36], Qiao et al. [38] expressed the vertical mixing, Bv, caused by nonbreaking waves as follow:
B v = α E ( k ) e x p { 2 k z } d k z ( k ω 2 E ( k ) e x p { 2 k z } d k ) 1 2 ,
where the parameter α = ο ( 1 ) and is set to 1 following [36,38]; E ( k ) represents the wave number spectrum, which can be calculated from the numerical wave model; ω is the angular frequency of the wave; k is the wave number, and z is the water depth. The nonbreaking WITM is directly incorporated into the ocean model as follows:
K m = K m c + B v , K h = K h c + B v
where Km and Kh are the vertical viscosity and diffusion coefficients respectively. Kmc and Khc are calculated by the Mellor–Yamada turbulence closure scheme employed in the ocean model.

2.2. Wave Transport Flux Residue

Based on the differences in physical characteristics and time-space scales, Yuan et al. [39] divided the ocean motions into four types: turbulence, wave-like motion, eddy-like motion, and circulation. At the same time, a group of Reynolds averages, including first-fold marked by   S S , second-fold marked by   S M , and third-fold marked by   M M , is defined to decompose the ocean motions. For example, by applying the three Reynolds average operations to the velocity, we can decompose it into four parts as shown in Equations (3)–(6),
U = U ˜ + u S S
U ˜ = U ^ + u S M
U ^ = U ¯ + u M M
U = U ¯ + u M M + u S M + u S S
where the U is the total ocean motion, the U ˜ is the sums of motions except ocean turbulence, the U ^ is the sums of the circulation, eddy-like motion, and the R.H.S. of the Equation (6) represent the circulation, eddy-like motion, wave-like motion, and turbulence, respectively. Moreover, with the first two Reynolds average operations and the ocean motions governing equations, we can get the governing equations for the sum of ocean eddy-like motion and circulation:
U ^ i x i = 0 ,
U ^ i t + U ^ j U ^ i x j 2 ε i j k U ^ j Ω k = 1 ρ 0 p ^ x i g ρ ^ ρ 0 δ i 3 + x j ( ν 0 U ^ i x j ) + x j [ ( u S M j u S M i S M + u S S j u S S i S S S M ) ]
T ^ t + U ^ j T ^ x j = x j ( κ 0 T ^ x j ) + x j [ ( u S M j T S M S M + u S S j T S S S S S M ) ] + Q ^
s ^ t + U ^ j s ^ x j = x j ( D 0 s ^ x j ) + x j [ ( u S M j s S M S M u S S j s S S S S S M ) ]
ρ ^ = ρ ( s ^ , T ^ , p ^ )
and the governing equations for the wave-like ocean motion sub-system:
u S M i x i = 0 ,
u S M i t + U ^ j u S M i x j + u S M j U ^ i x j + x j ( u S M j u S M i u S M j u S M i S M ) 2 ε i j k u S M j Ω k = 1 ρ 0 p S M x i g ρ S M ρ 0 δ i 3 + x j ( ν 0 u S M i x j )   + x j [ ( u S S j u S S i S S u S S j u S S i S S S M ) ]
T S M t + U ^ j T S M x j + u S M j T ^ x j + x j ( u S M j T S M u S M j T S M S M )         = x j ( κ 0 T S M x j ) + x j [ ( u S S j T S S S S u S S j T S S S S S M ) ] + Q S M
s S M t + U ^ j s S M x j + u S M j s ^ x j + x j ( u S M j s S M u S M j s S M S M )         = x j ( D 0 s S M x j ) + x j [ ( u S S j s S S S S u S S j s S S S S S M ) ]
ρ S M = ρ ( T ^ + T S M , s ^ + s S M , p ^ + p S M ) ρ ( T ^ , s ^ , p ^ )
Here x i ( i = 1 , 2 , 3 ) denote the eastward, northward, and vertical of the coordinate, respectively. u S M i ( i = 1 , 2 , 3 ) , T S M , s S M , p S M , ρ S M are the velocity, temperature, salinity, pressure, and density for the wave-like ocean motion, and U ^ i ( i = 1 , 2 , 3 ) , T ^ , s ^ , p ^ , ρ ^ for the sum of ocean eddy-like motion and circulation, which also act as the background field components in the Equations (12)–(16). ν o , κ o , D o are the coefficients of molecular viscosity, thermal conductivity, and the diffusion coefficient of salinity. Q ^ and Q S M are the equivalent temperature source terms. It consists of four components which are net downward shortwave, net downward longwave, sensible, and latent fluxes, and is calculated using the Coupled Ocean–Atmosphere Response Experiment (COARE) bulk formula [40]. The latter three components are assumed to occur at the ocean surface through the surface boundary conditions. The downward shortwave flux is specified as an absorption profile which is interpreted by Paulson and Simpson [41]. Each sub-system has two effects in the ocean dynamic cascade system. The governing equation group of each sub-system gives the expression of these two effects. The first is the advection transport and shear instability generations of larger scale motions through which modulate smaller scale motions (the second and third terms on L.H.S. of Equations (13)–(15)). The large-scale variables in the sub-system are solved in its superior system. The second is the ocean mixing induced by smaller scale motions in the form of transport flux residues through which modulate larger scale motions (the last term on R.H.S. of Equation (8), and the second term on R.H.S. of Equations (9) and (10). The small-scale variables in the sub-system are solved in its subordinate system. The transport flux residue expresses the amount left after applying the Reynold’s average on the flux of small scale motions [39]. The effects of wave-like ocean motions on circulation, we call wave transport flux residue (hereafter WTFR), include:
x j [ u S M i u S M j S M T S M u S M j S M s S M u S M j S M ]
Furthermore, considering the background field shear, we decomposed the wave-like ocean motions into two components, i.e.,
u S M i = u ˜ S M i + u ^ S M i ,   ( i = 1 , 2 , 3 )
T S M = T ˜ S M + T ^ S M
s S M = s ˜ S M + s ^ S M
where u ˜ S M i ( i = 1 , 2 , 3 ) , T ˜ S M , s ˜ S M are the velocity, temperature, and salinity of the dominant gravitational wave components, and u ^ S M i ( i = 1 , 2 , 3 ) , T ^ S M , s ^ S M are the perturbation components generated by background field shear. With Equations (17)–(19), the wave transport flux residue can be expanded as:
( u S M i u S M j ) x j S M = { ( u ˜ S M i u ˜ S M j ) x j S M + ( u ˜ S M j u ^ S M i ) x j S M + ( u ^ S M j u ˜ S M i ) x j S M + ( u ^ S M i u ^ S M j ) x j S M } , i = 1 , 2 , 3
( u S M i T S M ) x i S M = { ( u ˜ S M i T ˜ S M ) x i S M + ( u ˜ S M i T ^ S M ) x i S M + ( u ^ S M i T ˜ S M ) x i S M + ( u ^ S M i T ^ S M ) x i S M } ,
( u S M i s S M ) x i S M = { ( u ˜ S M i s ˜ S M ) x i S M + ( u ˜ S M i s ^ S M ) x i S M + ( u ^ S M i s ˜ S M ) x i S M + ( u ^ S M i s ^ S M ) x i S M } .
On the R.H.S. of Equation (20), the first term is the radiation stress and was widely used. Here, the ocean wave field is assumed to be locally uniform, i.e., T ˜ S M = 0 ,   s ˜ S M = 0 , so that the first and third terms on the R.H.S. of Equations (21) and (22) can be neglected. Moreover, the last terms on the R.H.S. of Equations (20)–(22) are small terms. The role of the third term on the R.H.S. of Equation (20) is not investigated in this study. The second terms on the R.H.S. of Equations (20)–(22) are the main components of WTFR. With the continuity equation for u ˜ S M i ( i = 1 , 2 , 3 ) , they can be expressed as:
( u ˜ S M j u ^ S M i ) x j S M = u ˜ S M j u ^ S M i x j S M , i = 1 , 2 , 3
( u ˜ S M i T ^ S M ) x i S M = u ˜ S M i T ^ S M x i S M
( u ˜ S M i s ^ S M ) x i S M = u ˜ S M i s ^ S M x i S M
Since ocean waves are locally uniform, the horizontal changes of statistic parameters for ocean waves within the mixing length are near zero [36]. Therefore, only the vertical terms in WTFR are relevant, and the final expressions are given in Equations (26)–(28). The more detailed derivation processes are summarized in Yang et al. [42].
u ˜ S M 3 u ^ S M j x 3 S M = U ^ j x 1 k ω k 1 E ( k ) cosh 2 K ( H + x 3 ) sinh 2 K H d k U ^ j x 2 k ω k 2 E ( k ) cosh 2 K ( H + x 3 ) sinh 2 K H d k , j = 1 , 2 , 3
u ˜ S M 3 T ^ S M x 3 S M = T ^ x 1 k ω k 1 E ( k ) cosh 2 K ( H + x 3 ) sinh 2 K H d k T ^ x 2 k ω k 2 E ( k ) cosh 2 K ( H + x 3 ) sinh 2 K H d k
u ˜ S M 3 s ^ S M x 3 S M = s ^ x 1 k ω k 1 E ( k ) cosh 2 K ( H + x 3 ) sinh 2 K H d k s ^ x 2 k ω k 2 E ( k ) cosh 2 K ( H + x 3 ) sinh 2 K H d k
where H is the local bottom topography. The transport flux residue represents the effect of small-scale motion to large scale motion and is dependent linearly on the background field shear, so this value may sometimes become negative. In positive value condition, the WTRF extract momentum, heat, or salinity from the atmosphere to the ocean. The negative value corresponds to upward flux transfer, i.e., from ocean to atmosphere. In other words, the WTFR may uniform or stratify the ocean horizontally or vertically.
Here we define the vertical WTFR terms in momentum, temperature, and salinity equations:
B S M 1 = u ˜ S M 3 u ^ S M 1 x 3 S M , B S M 2 = u ˜ S M 3 u ^ S M 2 x 3 S M , B S M T = u ˜ S M 3 T ^ S M x 3 S M , B S M S = u ˜ S M 3 s ^ S M x 3 S M .
The effects of the WTFR are introduced into the ocean model as follows:
U ^ D t + A d v ( U ^ ) D i f ( U ^ ) f V ^ D + 1 ρ ^ p ^ x = σ ( K M D U ^ σ ) + B S M 1 D
V ^ D t + A d v ( V ^ ) D i f ( V ^ ) + f U ^ D + 1 ρ ^ p ^ y = σ ( K M D V ^ σ ) + B S M 2 D
T ^ D t + A d v ( T ^ ) D i f ( T ^ ) = σ ( K H D T ^ σ ) R σ + B S M T D
S ^ D t + A d v ( S ^ ) D i f ( S ^ ) = σ ( K H D S ^ σ ) + B S M S D
The above four equations represent the momentum, temperature, and salinity equations in the ocean model. D is the total water depth, D = H + η where η is the surface elevation, σ is the vertical direction in the sigma coordinate, f is the Coriolis parameter, p is the pressure, and R is the short wave radiation flux. This scheme is applied in the real ocean situation simulation for the first time, especially under TC condition, we consider only the effect of the WTFR in the temperature equation, which is also beneficial to the analysis of the experimental results.

3. Materials and Methods

3.1. Observation Data

To more objectively analyze the impact of each physical process, we compared the experimental results with the observed data. This subsection lists the observed data we have used and provides a general introduction.

3.1.1. Jason-2

In this study, we used the well-calibrated Geographical Data Record (GDR, https://www.aviso.altimetry.fr) significant wave height (SWH) data measured by the Jason-2 satellite altimeter to test the simulation performance of the wave model. The satellite was launched by NASA and the European Organization for the Exploitation of Meteorological Satellites (EUMETSAT) and is used mainly for measuring sea level and SWH with a high degree of accuracy. These data have been widely used in marine and climate research [43,44]. The prelaunch SWH error budget specification for Jason-2 is 10% or 0.4 m [45]. The later validation studies show that Jason-2 measured SWH has a 2–6% difference compared to buoys data [46,47].

3.1.2. MW_IR

To test the performance of the ocean model, we used the optimally interpolated daily SST data product produced by the US Remote Sensing System (MW_IR, http://www.remss.com). These data combine the through-cloud capabilities of microwave data (MW) with the high spatial resolution of infrared (IR) data and have a resolution of nearly 9 km. The accuracy of the data is guaranteed through quality control and optimal interpolation schemes [48]. The measured SSTs are shown to have a mean bias of −0.07 °C and a standard deviation of 0.57 °C when compared to the buoy data [49].

3.1.3. Argo

The Array for Real-time Geostrophic Oceanography (Argo) floats measure the temperature and salinity of oceans (0–2000 m) globally, even under extreme atmospheric and oceanic conditions, such as a TC passage. Therefore, Argo data can well reflect the upper ocean structure and variability under the impact of a TC [50]. During the passage of TC Nepartak, an Argo float numbered WMO 2902666 was found, and it measured twelve profiles. Furthermore, another 69 profiles near the TC track were also used to study the responses of the upper ocean during TC (Figure 1b). The data were downloaded from the Coriolis Global Data Center in France (ftp://ftp.ifremer.fr).

3.1.4. Buoy

Yang et al. [51] described some very useful observations made a new kind of deep-sea data buoy that is able to survive encounters with TC Nepartak. The observed data are therefore quite valuable and provide a good picture of the rapid changes in the upper ocean temperature field that occur during the direct forcing phase of the TC. In this study, we digitized part of the image data in their paper and used it to compare with the experimental results. The locations of the two buoys, NTU1 and NTU2, are marked in Figure 1b with yellow and blue pentacle stars, respectively.

3.2. Model Descriptions

3.2.1. Ocean Model

The Princeton Ocean Model (POM; [52]) was used to simulate the circulation. POM is a three-dimensional model using Arakawa-C grid staggering, with sigma coordinates in the vertical direction to handle complex topographies and shallow regions. The turbulence closure scheme employed in the POM is the Mellor–Yamada scheme, which is based on the level-2.5 turbulent kinetic energy equations [2].
The model domain covers the geographical area ( 0 50 N ,   99 160 E , see Figure 1a), has a horizontal resolution of ( 1 / 30 ) × ( 1 / 30 ) that corresponds to 3.7 km in latitude and varies between 3.7 and 2.5 km in longitude. A total of 73 vertical sigma layers are included, and the vertical grid spacing had a higher resolution in the surface and bottom boundary layers and avoided artificial diffusion in large water depths. The highest vertical resolution is D × 3.6 × 10 4 meters and the lowest resolution is D × 5.5 × 10 2 meters, where the D is total water depth. The topography of the model is interpolated based on the General Bathymetric Chart of the Oceans (GEBCO_08), and its resolution is ( 1 / 120 ) × ( 1 / 120 ) . The minimum depth was set to 5 m, and the maximum depth was set to 5500 m. The mode splitting technique [53,54] is used in this model for the computer economy. The external model is 2-D with a time step of 3 s, and the internal model is 3-D uses a longer time step of 90 s. The boundary conditions are based on the ( 1 / 12 ) × ( 1 / 12 ) resolution HYCOM analysis data (http://hycom.org/dataserver/). The boundary forcing variables include sea surface level, velocity, salinity, and temperature. The initial condition was obtained from the model result of Zhao et al. [55] with the initial time 1 January 2013.
The model is forced by both wind stress and surface heat fluxes. According to the different types of forcing data, the simulation procedure includes three parts. In the first (from the year 2013 to 2015), Global Forecast System (GFS) wind data provided by the National Centers for Environmental Prediction (NCEP) are used, and wind stresses were computed using a bulk formula. The data set has spatial and temporal resolutions of 0.25 × 0.25 and 3 h, respectively. The drag coefficient formulation which can be used for low-to-moderate winds and high wind speed is adopted [56]. The climatological monthly mean data from the Comprehensive Ocean-Atmosphere Data Set (COADS; [57]) are used as the surface heat fluxes. The model is integrated for 3 years reaching a stable state. For the second part (the year 2016), the climatological monthly mean surface heat fluxes are replaced by real time data. This part is used to assess the performance of both the wave and ocean model, and also to evaluate the impact of the two schemes on the ocean superstructure under general condition. The GFS meteorological variables are used to calculate the surface heat fluxes using the Coupled Ocean–Atmosphere Response Experiment (COARE) bulk formula [40]. For the last part (during TC Nepartak), to better and emphatic analyze the influence of the two schemes during the TC period, a wind field and wind pressure model is used to construct the wind field and pressure field during the TC Nepartak which will be discussed in Section 3.2.3.

3.2.2. Wave Model

The Marine Science and Numerical Modeling (MASNUM) wave model is a third-generation wave model that was originally designed by Yuan et al. [58,59]. Subsequently, Yang et al. [60] updated the model to increase its suitability for the global scale. The model solves the energy balance equation through the wave number spectrum space. The complicated characteristics inlaid scheme is used for wave energy propagation. Source functions including wind input, dissipation, wave–wave interaction, and wave–current interaction are introduced in this wave model. The wind input source function is based on Phillips’ resonance and Miles’ shear flow instability mechanisms [61,62]. The wave spectrum of the model is discretized using 24 directions and 25 frequency bins with a minimum frequency of 0.0412Hz. The frequency increases logarithmically, and the increment factor is 1.21. The in situ measurement based formulation of Wu [63], which is applicable even in hurricanes, is adopted as drag parameterization in MASNUM. The model domain and resolution of the wave model are consistent with that of the POM. The boundary conditions of the wave model are from the JONSWAP spectrum [64]. The MASNUM wave model is integrated for 4 years independently, from the year 2013 to 2016, with the same period as the ocean model. The surface wind stress data used by the wave model are also the same as those used by the ocean model.
In the first two parts, the WITM (Bv, see Equation (1)) and the integral part of WTFR (Bsmt, see Equation (27)) are generated as daily output and introduced into the ocean model subsequently. For the third part, the data output and transferring frequency are increased to 3 h equal to the surface forcing data temporal resolutions.

3.2.3. Wind Field and Wind Pressure Model

During TC Nepartak, the wind and pressure fields are generated by a typhoon wind model. The input conditions for the model are from the Joint Typhoon Warning Center (JTWC, https://metoc.ndbc.noaa.gov/ja/web/guest/jtwc/best_tracks), including the features of the TC, such as warning time, estimation of position, maximum sustained surface wind speed, minimum central pressure, propagation speed and direction, and radius of maximum wind speed. The wind field and wind pressure model gives the pressure P at any radius (r) from the typhoon center as follows [65,66]:
{ P ( r ) P 0 P P 0 = 1 1 1 + 2 ( r R ) 2 , 0 r 2 R P ( r ) P 0 p p 0 = 1 1 1 + r R , 2 R < r <
where P is the ambient pressure; P 0 is the central pressure; P ( r ) is the surface pressure at a distance r from the typhoon center and R is the radius of maximum wind speed.
Using Equation (34), the gradient wind speed can be calculated via gradient wind theory:
V g ( r ) = { f 2 + f 2 4 + 10 3 2 Δ P ρ a R 2 [ 1 + 2 ( r 2 R 2 ) ] 3 2 , 0 r 2 R f 2 + f 2 4 + 10 3 Δ P ρ a ( 1 + r R ) 2 R r , 2 R < r <
where Δ P is the difference between the ambient pressure and the center pressure; f is the Coriolis parameter; and ρ a is the air density (assumed to be constant at 1.2 k g m 3 ).
The wind field derived from the Equation (35) is symmetrical to the typhoon center. Considering the asymmetry of the real wind field, two other factors, the angle of inflow and typhoon forward motion, should be included. In this study, the angle of inflow is set to be constant ( θ = 20 ° ), and the formula proposed by Ueno [67] is employed to calculate the forward wind field:
V f = V x e x p ( π 4 | r R | R ) i + V y e x p ( π 4 | r R | R ) j
where V x and V y are the components of the forward speed of the typhoon speed, with V x being positive eastward and V y being positive northward; and i and j are the unit vectors in the x and y directions, respectively. The model’s wind field comes to:
V M = c 1 V f + c 2 { V g [ ( x x 0 ) sin θ + ( y y 0 ) cos θ ] i + V g [ ( x x 0 ) cos θ + ( y y 0 ) sin θ ] j }
where ( x 0 , y 0 ) is the typhoon center position and ( x , y ) is the calculated position; c 1 and c 2 are the correction coefficients, set to be 1.0 and 0.8, respectively.
The constructed wind field for typhoon is combined with the background wind field as follow [68]:
V C = ( 1 e ) V M + e V Q , e = c 4 1 + c 4 , c = r n R
where V C is the combined wind field and V Q is the background wind field, and n is a constant, set to be 10. The GFS wind field is employed as the background wind field.

3.3. Overview of the Typhoon and Experimental Design

For a better interpretation of the results to be reported soon, it is necessary to provide a short description of the TC that we have considered for the present study. In the next subsection, we describe the experiments we have carried out to test the sensitivity of the model results to the two physical processes described in Section 2.

3.3.1. Overview of the Typhoon

At 1600 UTC 2 July 2016, Nepartak was formed on the surface of the Northwest Pacific Ocean; it then intensified into a strong tropical storm at 1200 UTC 3 July 2016, strengthened at 1800 UTC 4 July 2016 to a typhoon, and then further strengthened at 1800 UTC 4 July 2016 to a super typhoon. At 1200 UTC 6 July 2016, Nepartak increased to its peak intensity of 79.9 m/s and 907 hPa; at 0550 UTC 8 July 2016, it landed at Taitung County, Taiwan Province; at 1345 UTC 9 July 2016, it landed again at Shishi City, Fujian Province; and at 0600 UTC 9 July 2016, it then weakened to a tropical depression (Figure 1b).
TC Nepartak had a strong intensity and reached a super typhoon level, which had a drastic impact on the upper ocean structure. During the TC period and a long time thereafter, the TC transit area was not affected by other TCs or strong weather events, which was more conducive to the analysis of the experimental results.

3.3.2. Experimental Design

To study the effects of WITM and WTFR on ocean vertical mixing, especially during the TC stage, we designed three experiments, denoted as EX1, EX2, and EX3. EX1 is the control experiment, which is the original ocean current model; EX2 introduces WITM into the ocean current model; EX3 introduces both two physical processes (Table 1).

3.4. Measures of Quality for Models Validation

For quantitative comparison, five statistical parameters, including mean error (ME), root mean square error (RMSE), correlation coefficient (r), Nash–Sutcliffe efficiency coefficient (E), and determination coefficient (R2) were used and briefly summarized as below:
M E = 1 N i = 1 N ( M i O i )
R M S E = 1 N i = 1 N ( M i O i ) 2
r = i = 1 N ( M i M ¯ ) ( O i O ¯ ) i = 1 N ( M i M ¯ ) 2 i = 1 N ( O i O ¯ ) 2
E = 1 i = 1 N ( O i M i ) 2 i = 1 N ( O i O ¯ ) 2
R 2 = i = 1 N ( M i M ^ i ) i = 1 N ( M i M ¯ i )
where M represents the model-generated data, O represents the observational data, the bar over M and O denotes their mean values, M ^ denotes the fitted value, and N is the number of collocated data points. The ME is a statistical quantity that reveals any systematic error between the model and buoy data. The RMSE indicates the degree of discretization of data pairs; lower RMSE values indicate a smaller degree of data discretization between the model and observations. The r value indicates how well the model data fit the observations. E is generally used to verify the model simulation results. The value of E is minus infinity to 1, and E is close to 1, indicating the good quality of the model and high reliability of the model. R2 is used to judge the quality of data fitting. Its value range is between 0 and 1, and the larger value shows a higher fitting quality.

4. Results

In this section, first, the comparative results for the entire year of 2016 are used to assess the performance of both the wave and ocean model, the latter also indicates the contribution of the two physical processes to the improvement of the ocean model simulation capability under general condition. Second, a more detailed comparison has been performed using the simulation results of the three experiments and the observation data recorded during the TC Nepartak. The model data are daily averaged during TC period to correspond with the observations which are available as a single data set for each day.

4.1. Evaluation Under General Condition: The year 2016

First, we use the SWH data of the Jason-2 satellite to check the simulation results of the wave model. The scatter distribution statistics of the SWH for the year 2016 in the simulated area are shown in Figure 2a. Most of the data points lie near the diagonal line, indicating that the simulation results are consistent with the satellite observations. The error probability distribution (Figure 2b) shows that as the error increases, the probability distribution decreases rapidly, and most of the errors are within an acceptable range of 0.4 m. The r, E, and R2 value are about 0.95, 0.77, and 0.79, respectively (see Table 2). The statistic comparisons between simulated SWH and the observations show MASNUM’s good simulation capabilities which provide a good basis for our study of the mixing effects of waves.
The box plot of the spatially averaged SST error in the simulated area is given in Figure 3. The observed data used in the comparison are the MW_IR remote sensing observation SST data. The box plot shows the maximum value (upper whisker line), minimum value (lower whisker line), and outliers (external points) of a set of data. The median line inside the box, the bottom of the box and top of the box represents the median, lower quartile, and upper quartile of the data, respectively.
Figure 3 shows that the average error of the SST simulated by the three experiments is within 1 °C, and the median error of each month is within 0.5 °C, which is relatively small. Except for January, the mean SST of each experiment is higher than the observations, with a positive error. However, compared with EX1, the error range of EX2 and EX3 is closer to the zero line, especially in summer (June, July, and August) and early autumn (September). After introducing the mixing effect of the waves, the two experiments (EX2, EX3) show similar simulation results, and the difference between the model simulated SST and the observation is decreased by more than 10% (Table 2). The results of statistical analysis show that the E and R2 value is increased from 0.76 to 0.84 and 0.78 to 0.89, which indicate an improvement of the ability of the model to simulate SST.
The above analyses show that both wave and ocean models perform well and the gradual introduction of two physical processes (in EX2 and EX3) has played a positive role in improving the ocean model simulation capability as a whole. The error of SST has been improved to varying degrees, most notably in summer and early autumn. This is further discussed in Section 5.

4.2. Evaluation Under TC Condition: Nepartak

The simulated wave field is evaluated during the TC period using the orbit data of the Jason-2 satellite (Figure 4). The left panel shows the simulated significant wave height at 0000 UTC 04 July 2016 when Nepartak is at the tropical storm level and the white line shows the TC track. The right panel shows the comparison between the model SWH results and observations during 0424 and 0427 UTC 04 July (tandem track shown in the left panel using gray line). The model results show reasonable agreement with the observations under height wind condition, with a −0.15 m ME (The statistical area contains the latitude range in Figure 4b).
The response of the SST on 11 July 2016 is presented in Figure 5. For this study, similar to Zhang, et al. [69], the response of SST was defined as the SST change relative to its value at pre-typhoon stage (1 July 2016). The model data are interpolated to the observation points before plotting results. All experimental and observation results show a strong SST response to the TC. Comparing with the EX1, the SST cooling area in EX2 and EX3 is significantly increased, which is more consistent with the observations. A noteworthy phenomenon is that the maximum cooling area in observation as well as in EX3 is located to the R.H.S. of the TC path, which is well demonstrated in field observation [70,71] and multiple numerical experiments [25,72,73]. On the right side of TC track, the resonance between wind and current rotations and also between the wind and the TC translation result in larger waves which increase the turbulence and therefore enhance the SST cooling ultimately. However, in EX2, the maximum cooling area shows a left side shift tendency. After introducing the effect of WTFR (Figure 5c), the distribution of the maximum cooling area is considerably improved.
In general, the above comparisons indicate that the ability of the model to simulate SST during the TC gradually increased via the introduction of two different wave effects. The WITM remarkably improves the problem associated with an insufficient SST cooling area and cooling amplitude, although its impact on the distribution of the cooling area is relatively small. The phenomenon, i.e., left shift of the maximum cooling area, can be well modified after the introduction of WTFR.
To clearly show the variations of the upper ocean during TC and evaluate the role of wave mixing in these processes, observations from the two buoys, NTU1 and NTU2, are used [51]. Figure 6a and Figure 7a show the comparison of time series wind and pressure data between observed and model used at NTU1 and NTU2, respectively. The constructed wind and pressure data show great consistency with the observation, although the wind data are flatter to the observed. The upper 25 m-averaged temperature (T25m) obtained from the observations and the simulations are compared (Figure 6b and Figure 7b). The model output results in every 3 h, and the time axis indicates the output time. The comparison shows that with the introduction of wave-induced mixing, there is a marked drop in temperature. Compared with EX1, the results of EX2 have been improved to some extent, but they are still not sufficient; for example, the RSME is decreased from 1.92 °C to 1.43 °C at NUT1 and from 1.67 °C to 0.99 °C at NTU2 (Table 3). EX3’s results are more consistent with the observations, and the RSME becomes 0.68 °C and 0.45 °C at NTU1 and NTU2, respectively. It is worth noting that the EX1 and EX2 give a negative E value at two buoy locations and the model simulation results are unsatisfactory. In EX3, the model performance is greatly improved, and the E value is 0.61 and 0.74, respectively. Figure 6c and Figure 7c are the temperature profiles at the corresponding moment ticked in Figure 6b and Figure 7b which show that the effect of wave-induced mixing is mainly concentrated on the upper 150 m. It can be seen that among the three experiments, the results of the EX3 are more in line with the observations.
Here, we compared the surface heat flux between the model calculated and observed at the buoy station NTU1 and NTU2 (Figure 8). Before and after the effect of TC, the difference of the calculated surface heat flux between the three experiments is relatively small, and the data are in good agreement with the observations. During the TC period, the calculated data are more accurate in the experiments which included the wave effect (EX2 and EX3).
On the other hand, the response of the upper ocean to the TC is also reflected in the changes of the MLD and MLT. To clearly show this structure and its changes as well as the role of wave mixing in these processes, we use Argo float data (Figure 1b) to test the simulation results of the three experiments. The MLD is defined as the depth at which the temperature is 0.5 °C lower than the first layer. The MLT is then defined as the average temperature in this depth range.
The observational results show that at the Argo float WMO 2902666 (the cross in Figure 1b), the MLD rapidly deepens from 25 m due to TC Nepartak and reaches a maximum depth of 55 m (Figure 9a, black line). The MLD in all three experiments shows similar variation trends as the observation. However, the variation in EX1 is too weak and extends from 26 m to a maximum depth of 42 m, showing an insufficient response of the upper ocean layer to the TC. With the effect of the WITM (EX2), the maximum MLD reached 54 m, which is very close to the observed value. The result of EX3 accorded basically to the result by EX2 with only subtle differences in the subsequent oscillation stage. The observational and all three experimental results show that the MLD return to its pre-storm depth approximately 7 days after the passage of the TC (14 July). By comparing three experiments results with the observation, it is found that the WITM has an important effect on the formation of MLD, especially at the deepening stage.
The observed MLT decreased rapidly from 29.9 °C (before the TC impact) to 27.3 °C (Figure 9b), with the maximum cooling exceeding 2 °C. The pre-typhoon MLT of EX1 is 31.1 °C, which is 1.2 °C higher than the observed. Affected by the TC, the minimum MLT declined to 30 °C, with a maximum decrease of 1.1 °C. After considering the effect of WITM (EX2), the simulated MLT cooling amplitude increased to 1.3 °C, decreased from 30.3 °C to 29 °C. With WTFR effect (EX3), the performance of the MLT simulation is further improved with a maximum decrease of 2.1 °C. In EX3, the MLT warmed up to the pre-storm status until 24 July which shows a good agreement with the observation. However, the warming trend of the MLT in EX1 is not obvious and gradually returned to the pre-storm level until the last day of July (purple cross).
Below, three experimental results are compared with another 69 Argo observations during TC Nepartak to further study the effects of the two physical processes on the upper ocean (Figure 10). Shallow observed MLD (10~30 m) and high observed MLT (warmer than 30 °C) which we considered as storm-free status. All three experimental results are relatively consistent with the observation at that status. Affected by the TC, as observed MLD deepened and MLT decreased, three experiments show different performance. The response of the original ocean model (EX1) to the TC is insufficient and this phenomenon is gradually improved with the effects of the wave-induced mixing (EX2, EX3). Table 4 gives the comparisons of the MLD and MLT between Argo float observations and three experimental results during the TC period. In terms of MLD, EX1 presents ME of −8.97 m. In EX2, ME has fallen by more than 68% and EX3 shows a similar result. For the MLT, in the process of introducing the WITM (EX2) and WTFR (EX3), the model results are gradually improved. Especially, in EX3, the ME is decreased by 40% from 0.3 °C to 0.12 °C. WITM can improve the simulated MLD result observably, the E value is increased from 0.32 (EX1) to 0.90 (EX2), while WTFR plays a crucial role in the MLT changing, the E value is increased from 0.49 (EX1) to 0.75 (EX3). This is consistent with the previous results at Argo float WMO 2902666 (Figure 9).
Furthermore, we analyzed the effects of wave-induced mixing on the vertical distribution of temperature. Temperature profiles modeled by the EX2 (Figure 11a) and EX3 (Figure 11b) have been extracted perpendicular to the TC track (shown by a black dashed line in Figure 11b). The difference in temperature between with (EX2, EX3) and without (EX1) wave-induced mixing effects are presented in Figure 11c,d, respectively. The effects of wave-induced mixing show a hamburger-like construction with the upper and lower layer to be a cooling area and the middle layer to be a warming area (Figure 11c,d). The thickness of the upper layer extends to a depth of 50 m which is equivalent to the MLD approximately. These findings are well consistent with the previous study results [31,74]. The effect of wave-induced mixing is more intense on the right side of the track (the black dashed line showed the TC track), while it is weaker on the left side. Moreover, it is shown that the WTFR exerted a more dramatic cooling effect on the right side of track compared with that of the WITM effect.

5. Discussion

In this paper, using the MASNUM wave model and the POM, we studied the impact of two different wave-induced mixing effect processes, wave-induced turbulence mixing (WITM) and wave transport flux residue (WTFR), on the upper ocean structure especially during TCs in the Northwest Pacific Ocean. Unlike the traditional breaking wave studies, these two physical processes are both rely on the nonbreaking wave and are therefore independent of the wind stress.
By comparing the simulation results of the year 2016 with observation, it is shown that the original model tends to overestimate the SST, after including a consideration of the two wave effects, the simulation performance is remarkably improved. This improvement is most pronounced in summer and early autumn (Figure 3). We think this may benefit from the seasonal variation of the MLD. Observation data show that the global ocean has a deeper MLD in winter then in summer, while the spring and fall seasons act as transitional periods [75]. With the same intensity of wave-induced vertical mixing, a shallower MLD will make it easier for the colder water from the thermocline to be entrained into the well-mixed ocean surface and ultimately decrease the SST. The mean SST error shows a negative value in January, this may be caused by the inaccuracy of the surface heat fluxes used during the first part of the simulation (from the year 2013 to 2015). At the start of the year 2016, the SST remains a negative error at that stage, though the real time surface heat fluxes data were used. In addition, the three experiments overestimated SST to varying degrees for most of the year 2016 except January. The sea ice component is ignored during the simulation, and this may be one of the possible reasons causing the SST overestimation. The rainfall effect is unconsidered during surface net heat flux calculating which may also affect the SST simulation results.
During the TC period, the upper ocean always shows a significant response to the air forcing, such as SST cooling, MLD deepening and MLT decreasing [25,31,76,77]. Our study shows that with the WITM effect (EX2) the simulated sea surface cooling area and the cooling amplitude are more consistent with observations. The WTFR resulted in a better match with the observation of cooling area distribution. One of the most important features is that, by its action, the left deviation tendency of the maximum cooling area is suppressed (Figure 5). Here, we give the surface Bsmt distribution (Figure 12) in this area (indicated as a black dashed box in Figure 5b) on 9~11 July. The Figure shows that the left side of the typhoon track is dominated by positive value, indicating a warming effect which inhibits the left deviation tendency of the maximum cooling area. On the contrary, a negative value is observed on the right and means a cooling effect.
For the mixed layer affected by TC, due to the intense and quick effect of the WITM, the MLD rapidly deepened to a similar level of observed which is severely underestimated in the original ocean current model (Figure 9a). WTFR has no great influence on the MLD deepening. However, by introducing the WTFR, the model shows a more realistic MLT cooling and warming period compared with that of the other two experiments (Figure 9b). The results show that the MLD return to its pre-storm depth approximately 7 days after the passage of the TC (14 July), and a longer period of almost twenty days is needed for the MLT recovery. This is similar to the findings of Wu and Chen [78], who also argued that the asymmetric effects of entrainment and detrainment on the ML are responsible for the prolonged time required to restore the MLT compared with that of the MLD. Comparing among three experiments, it is shown that the wave-induced mixing effects on the upper ocean temperature have a three-layer structure (Figure 11), which is cooling, warming, and cooling. This effect is more intense on the right side of the TC track. Especially, WTFR has a more dramatic cooling effect on the right side of track compared with that of WITM.
In the wake of TC, regarded as a “relaxation stage”, the near-inertial waves are excited and the energy is dispersed in different vertical modes [79]. It is a non-negligible fraction of the wind power input and can be radiated down into the interior ocean [80,81]. The inertial response of the ocean to the TC is an important process for subsurface mixing and even for determining large-scale patterns of ocean circulation (like meridional overturning circulation) [82]. Thus, it is worthwhile to perform some analyses of that aspect here. The response of near-inertial currents is more pronounced to the right side of the TC track because of the resonance between the wind stress vector and inertial oscillations [72,83]. In Figure 13, we show the vertical section of the zonal velocity of upper 1000 m at Argo float 2902666 (Figure 1b, cross mask, on the right side of TC track) as a function of time, which is obtained from the present three numerical experiments. Modeled results show a near-inertial oscillation after the passage of TC, and the phase lead of the velocity oscillations clearly increased with depth, indicating downward radiation of energy. It can be seen that with the introduction of the WITM (EX2, Figure 13b) and WTFR (EX2, Figure 13c), the downward radiation velocity obviously increased. Moreover, EX3 shows the strongest oscillation amplitude among the three experiments. The near-inertial oscillation is gradually decayed after five inertial periods, which is consistent with previous findings [83]. Comparing the results of different experiments, EX1 has the shallowest radiation depth, limited in the upper 600 m. With the effect of the wave-induced mixing, the radiation depth is continually increased to almost more than 1000 m in EX3, which is more consistent with many observations [80,84,85]. For instance, Cuypers et al. [80] studied the impact of storm-induced near-inertial internal waves in the upper ocean and found that the near-inertial internal wave packets were identified down to 1000 m, the maximum depth of the measurements. The EX3’s result is in better agreement with their findings, which indirectly reflects the effect of WTFR on the vertical momentum distribution and the mixing processes.
The inertial response of the ocean to the TC is an important process for subsurface mixing. By experimental comparison, it is found that the wave could also exert a considerable effect on this phenomenon. Generally speaking, the wave-induced mixing not only increased the oscillation amplitude and downward radiation depth but also speed the downward radiation velocity. The near-inertial response is significantly enhanced which is more consistent with the previous studies.
Near the surface, the momentum flux (total stress) transferred from the atmosphere to the ocean. Previous studies [86,87] show that the total stress can be divided into wave- and turbulence-generation components. Similarly, the Reynold flux of the temperature is also the combination of these two parts. Because of the experimental design, we will discuss the distribution of the surface temperature Reynold flux here. The total Reynold flux of temperature in this study can be expressed as follows:
F t o t = z ( K h c T z ) + z ( B v T z ) + B S M T = F K h c + F B v + F B s m t
where the first term ( F K h c ) is the turbulence-induced component and the last two terms ( F B v ,   F B s m t ) are the wave-induced components. The EX1, EX2, and EX3 contain the first term, the first two terms, and all three terms, respectively. The total temperature Reynold flux of the three experiments (Figure 14a,d,h) show that the main cooling area is basically located in the right front of the typhoon. Figure 14a also represents the horizontal distribution of the turbulence-induced surface temperature Reynold flux, derived from the EX1. Within the typhoon eye, there exists a near-zero flux region due to the low wind speed [69]. The negative flux region that represents the cooling effect is fanned out toward the typhoon moving direction with the typhoon eye as the center, and the amplitude gradually decreases with the distance from the typhoon eye. In the EX2 and EX3 (Figure 14b,e), the turbulence-induced component is suppressed, and the surface temperature Reynold flux is dominated by the F B v (Figure 14c,f); i.e., the WITM plays a central role in surface cooling that is consistent with our above comparison results (Figure 5b). Comparing to the distribution of F B v , the negative value area of the F B s m t (Figure 14g) is relatively small, and the flux is also declined. However, the distribution of F B s m t is more complex than a simple negative distribution. There are many positive flux areas among the negative and, in particular, an evident positive flux stripe behind the typhoon eye. These special flux distributions can result in a more realistic sea surface cooling area distribution in EX3 (Figure 5c), where the WTFR is considered.
The exchange of momentum and heat between air and sea is a very complicated process, where many different processes are involved especially during the TC period. The physical conditions are always beyond the simplified structure defined in our modeling. In this study, a wind field and wind pressure model is applied to construct the surface forcing data during the typhoon instead of coupling an atmospheric model such as those in the study [14,55]. Although the comparisons showed that the constructed wind and pressure data have a great consistency with the observation, coupled atmosphere-wave-ocean modeling studies might lead to a more comprehensive understanding of the air–sea interaction. In future studies, we will perform further in-depth research to test the effect of the two physical processes on the upper ocean with a fully coupled system.

6. Conclusions

We investigated the effect of two wave-induced mixing processes (WITM and WTFR) on the upper ocean structure under both general and TC conditions. Including these two physical processes generally produces a reduced mean error of SST under general condition, especially in summer (June, July, and August) and early autumn (September) time. After the introduction of the wave effects, the performance of the model to simulate ocean response to TC could also be greatly improved. Specifically, the WITM can effectively improve the simulation results of the SST cooling area and amplitude and has an important effect on the formation of MLD. WTFR plays a positive role in reproducing the distribution of the SST cooling area and is dominated in the MLT changing. The analyses show that the two physical processes could also affect the ocean near-inertial response to the TC. The wave-induced mixing processes should be properly considered during the upper ocean structure modeling.

Author Contributions

Conceptualization, C.Y., Y.Y.; methodology, C.Y., Y.Y. and X.Y.; software, C.Y. and X.Y.; validation, C.Y. and X.Y.; formal analysis, C.Y.; investigation, C.Y.; resources, M.S. and X.Y.; data curation, M.S. and Y.S.; writing—original draft preparation, C.Y.; writing—review and editing, Y.Y. and X.Y.; visualization, C.Y.; supervision, X.Y.; project administration, Y.Y.; funding acquisition, Y.Y. and X.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was jointly funded by the National Key Research Program of China, grant number 2016YFC1402004 and 2017YFC1404201, Public Science and Technology Research Funds Projects of Ocean, grant number 201505013, and National Natural Science Foundation of China, grant number 41106032 and 41876003.

Acknowledgments

We acknowledge the data support from the Joint Typhoon Warning Center (JTWC, https://metoc.ndbc.noaa.gov/ja/web/guest/jtwc/best_tracks), CNES (https://www.aviso.altimetry.fr), Remote Sensing System (RSS, http://www.remss.com) and Coriolis Global Data Center in France (ftp://ftp.ifremer.fr).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kantha, L.H.; Clayson, C.A. An improved mixed layer model for geophysical applications. J. Geophys. Res. Oceans 1994, 99, 25235–25266. [Google Scholar] [CrossRef]
  2. Mellor, G.L.; Yamada, T. Development of a turbulence closure model for geophysical fluid problems. Rev. Geophys. Space Phys. 1982, 20, 851–875. [Google Scholar] [CrossRef] [Green Version]
  3. Rodi, W. Examples of calculation methods for flow and mixing in stratified fluids. J. Geophys. Res. 1987, 92, 5305–5328. [Google Scholar] [CrossRef]
  4. Umlauf, L.; Burchard, H. A generic length-scale equation for geophysical turbulence models. J. Mar. Res. 2003, 61, 235–265. [Google Scholar] [CrossRef]
  5. Wilcox, D.C. Reassessment of the scale-determining equation for advanced turbulence models. Aiaa J. 1988, 26, 1299–1310. [Google Scholar] [CrossRef]
  6. Kärnä, T. Discontinuous Galerkin discretization for two-equation turbulence closure model. Ocean Model. 2020, 150, 101619. [Google Scholar] [CrossRef]
  7. Reffray, G.; Bourdalle-Badie, R.; Calone, C. Modelling turbulent vertical mixing sensitivity using a 1-D version of NEMO. Geosci. Model Dev. 2015, 8, 69–86. [Google Scholar] [CrossRef] [Green Version]
  8. Zhang, Y.; Baptista, A.M. SELFE: A semi-implicit Eulerian–Lagrangian finite-element model for cross-scale ocean circulation. Ocean Model. 2008, 21, 71–96. [Google Scholar] [CrossRef]
  9. Kantha, L.H.; Clayson, C.A. On the effect of surface gravity waves on mixing in the oceanic mixed layer. Ocean Model. 2004, 6, 101–124. [Google Scholar] [CrossRef]
  10. Ezer, T. On the seasonal mixed layer simulated by a basin-scale ocean model and the Mellor-Yamada turbulence scheme. J. Geophys. Res. Oceans 2000, 105, 16843–16855. [Google Scholar] [CrossRef]
  11. Mellor, G.L. One-dimensional, ocean surface layer modeling: A problem and a solution. J. Phys. Oceanogr. 2001, 31, 790–809. [Google Scholar] [CrossRef]
  12. Ardhuin, F.; Rogers, E.; Babanin, A.V.; Filipot, J.F.; Magne, R.; Roland, A.; van der Westhuysen, A.; Queffeulou, P.; Lefevre, J.M.; Aouf, L. Semiempirical Dissipation Source Functions for Ocean Waves. Part I: Definition, Calibration, and Validation. J. Phys. Oceanogr. 2010, 40, 1917–1941. [Google Scholar] [CrossRef] [Green Version]
  13. Uchiyama, Y.; Mcwilliams, J.C.; Restrepo, J.M. Wave-current interaction in nearshore shear instability analyzed with a vortex force formalism. J. Geophys. Res. Oceans 2009, 114. [Google Scholar] [CrossRef] [Green Version]
  14. Kumar, N.; Voulgaris, G.; Warner, J.C.; Olabarrieta, M. Implementation of the vortex force formalism in the Coupled Ocean-Atmosphere-Wave-Sediment Transport (COAWST) modeling system for inner shelf and surf zone applications. Ocean Model. 2012, 47, 65–95. [Google Scholar] [CrossRef]
  15. Martin, P.J. Simulation of the mixed layer at OWS November and Papa with several models. J. Geophys. Res. Oceans 1985, 90, 903–916. [Google Scholar] [CrossRef]
  16. Canuto, V.M.; Howard, A.; Cheng, Y.; Dubovikov, M.S. Ocean Turbulence. Part I: One-Point Closure Model—Momentum and Heat Vertical Diffusivities. J. Phys. Oceanogr. 2001, 31, 1413–1426. [Google Scholar] [CrossRef]
  17. Harcourt, R.R. A Second-Moment Closure Model of Langmuir Turbulence. J. Phys. Oceanogr. 2013, 43, 673–697. [Google Scholar] [CrossRef]
  18. Harcourt, R.R. An Improved Second-Moment Closure Model of Langmuir Turbulence. J. Phys. Oceanogr. 2015, 45, 84–103. [Google Scholar] [CrossRef]
  19. Soloviev, A.; Lukas, R. Observation of wave-enhanced turbulence in the near-surface layer of the ocean during TOGA COARE. Deep Sea Res. Part I Oceanogr. Res. Pap. 2003, 50, 371–395. [Google Scholar] [CrossRef]
  20. Agrawal, Y.C.; Terray, E.A.; Donelan, M.A.; Hwang, P.A.; Williams, A.J.; Drennan, W.M.; Kahma, K.K.; Krtaigorodskii, S.A. Enhanced dissipation of kinetic energy beneath surface waves. Nature 1992, 359, 219–220. [Google Scholar] [CrossRef]
  21. Liu, Z.; Xu, J.; Zhu, B.; Sun, C.; Zhang, L. The upper ocean response to tropical cyclones in the northwestern Pacific analyzed with Argo data. Chin. J. Oceanol. Limnol. 2007, 25, 123–131. [Google Scholar] [CrossRef]
  22. Lai, Q.; Leiming, M.A.; Huang, W.; Liguang, W.U. The ocean response to Typhoon Morakot (2009) in the western North Pacific boundary region. Acta Oceanol. Sin. 2013, 35, 65–77. [Google Scholar]
  23. Liu, Z.; Jianping, X.U.; Sun, C.; Xiaofen, W.U. An upper ocean response to Typhoon Bolaven analyzed with Argo profiling floats. Acta Oceanol. Sin. 2014, 33, 90–101. [Google Scholar] [CrossRef]
  24. Bender, M.A.; Ginis, I.; Kurihara, Y. Numerical simulations of tropical cyclone-ocean interaction with a high-resolution coupled model. J. Geophys. Res. Atmos. 1993, 98, 23245–23263. [Google Scholar] [CrossRef]
  25. Wang, G.; Qiao, F. Ocean temperature responses to Typhoon Mstsa in the East China Sea. Acta Oceanol. Sin. 2008, 27, 26–38. [Google Scholar]
  26. Yang, Y.; Qiao, F.; Xia, C.; Jian, M.; Yuan, Y. Wave-induced mixing in the Yellow Sea. Chin. J. Oceanol. Limnol. 2004, 22, 322–326. [Google Scholar]
  27. Qiao, F.L.; Ma, J.; Xia, C.; Yang, Y.; Yuan, Y.L. Influences of the surface wave-induced mixing and tidal mixing on the vertical temperature structure of the Yellow and East China Seas in summer. Prog. Nat. Sci. Mater. Int. 2006, 16, 739–746. [Google Scholar]
  28. Xia, C.; Qiao, F.; Yang, Y.; Ma, J.; Yuan, Y. Three-dimensional structure of the summertime circulation in the Yellow Sea from a wave-tide-circulation coupled model. J. Geophys. Res. Oceans 2006, 111. [Google Scholar] [CrossRef]
  29. Song, Z.; Qiao, F.; Yang, Y.; Yuan, Y. An improvement of the too cold tongue in the tropical Pacific with the development of an ocean-wave-atmosphere coupled numerical model. Prog. Nat. Sci. Mater. Int. 2007, 17, 576–583. [Google Scholar]
  30. Toffoli, A.; McConochie, J.; Ghantous, M.; Loffredo, L.; Babanin, A.V. The effect of wave-induced turbulence on the ocean mixed layer during tropical cyclones: Field observations on the Australian North-West Shelf. J. Geophys. Res. Oceans 2012, 117. [Google Scholar] [CrossRef] [Green Version]
  31. Aijaz, S.; Ghantous, M.; Babanin, A.V.; Ginis, I.; Thomas, B.; Wake, G. Nonbreaking wave-induced mixing in upper ocean during tropical cyclones using coupled hurricane-ocean-wave modeling. J. Geophys. Res. Oceans 2017, 122, 3939–3963. [Google Scholar] [CrossRef] [Green Version]
  32. Huang, C.J.; Qiao, F.; Song, Z.; Ezer, T. Improving simulations of the upper ocean by inclusion of surface waves in the Mellor-Yamada turbulence scheme. J. Geophys. Res. Oceans 2011, 116. [Google Scholar] [CrossRef] [Green Version]
  33. Babanin, A.V.; Haus, B.K. On the existence of water turbulence induced by nonbreaking surface waves. J. Phys. Oceanogr. 2009, 39, 2675–2679. [Google Scholar] [CrossRef]
  34. Qiao, F.; Yuan, Y.; Jia, D.; Dai, D.; Song, Z. Wave–turbulence interaction-induced vertical mixing and its effects in ocean and climate models. Philos. Trans. R. Soc. A 2016, 374, 20150201. [Google Scholar] [CrossRef] [Green Version]
  35. Babanin, A.V. On a wave-induced turbulence and a wave-mixed upper ocean layer. Geophys. Res. Lett. 2006, 33, 382–385. [Google Scholar] [CrossRef]
  36. Qiao, F.; Yuan, Y.; Ezer, T.; Xia, C.; Yang, Y.; Lü, X.; Song, Z. A three-dimensional surface wave–ocean circulation coupled model and its initial testing. Ocean Dyn. 2010, 60, 1339–1355. [Google Scholar] [CrossRef]
  37. Dai, D.; Qiao, F.; Sulisz, W.; Han, L.; Babanin, A. An experiment on the nonbreaking surface-wave-induced vertical mixing. J. Phys. Oceanogr. 2010, 40, 2180–2188. [Google Scholar] [CrossRef]
  38. Qiao, F.; Yuan, Y.; Yang, Y.; Zheng, Q.; Xia, C.; Ma, J. Wave-induced mixing in the upper ocean: Distribution and application to a global ocean circulation model. Geophys. Res. Lett. 2004, 31, 293–317. [Google Scholar] [CrossRef]
  39. Yuan, Y.L.; Qiao, F.L.; Yin, X.Q.; Han, L.; Lu, M. Establishment of the Ocean Dynamic System with Four Sub-Systems and the Derivation of Their Governing Equation Sets. J. Hydrodyn. 2012, 24, 153–168. [Google Scholar] [CrossRef]
  40. Fairall, C.W.; Bradley, E.F.; Hare, J.E.; Grachev, A.A.; Edson, J.B. Bulk Parameterization of Air Sea Fluxes: Updates and Verification for the COARE Algorithm. J. Clim. 2003, 16, 571–591. [Google Scholar] [CrossRef]
  41. Paulson, C.A.; Simpson, J.J. Irradiance Measurements in the Upper Ocean. J. Phys. Oceanogr. 1977, 7, 952–956. [Google Scholar] [CrossRef]
  42. Yang, Y.; Shi, Y.; Yu, C.; Teng, Y.; Sun, M. Study on surface wave-induced mixing of transport flux residue under typhoon conditions. J. Oceanol. Limnol. 2019, 37, 1837–1845. [Google Scholar] [CrossRef]
  43. Ablain, M.; Cazenave, A.; Larnicol, G.; Balmaseda, M.; Cipollini, P.; Faugère, Y.; Fernandes, M.J.; Henry, O.; Johannessen, J.A.; Knudsen, P. Improved sea level record over the satellite altimetry era (1993–2010) from the Climate Change Initiative Project. Ocean Sci. Discuss. 2015, 11, 2029–2071. [Google Scholar] [CrossRef]
  44. Han, G.; Ma, Z.; Chen, N.; Chen, N.; Yang, J.; Chen, D. Hurricane Isaac storm surges off Florida observed by Jason-1 and Jason-2 satellite altimeters. Remote Sens. Environ. 2017, 198, 244–253. [Google Scholar] [CrossRef]
  45. Lambin, J.; Morrow, R.; Fu, L.-L.; Willis, J.K.; Bonekamp, H.; Lillibridge, J.; Perbos, J.; Zaouche, G.; Vaze, P.; Bannoura, W.; et al. The OSTM/Jason-2 Mission. Mar. Geod. 2010, 33, 4–25. [Google Scholar] [CrossRef]
  46. Abdalla, S.; Janssen, P.A.E.M.; Bidlot, J.R. Jason-2 OGDR Wind and Wave Products: Monitoring, Validation and Assimilation. Mar. Geod. 2010, 33, 239–255. [Google Scholar] [CrossRef]
  47. Ray, R.D.; Beckley, B.D. Calibration of Ocean Wave Measurements by the TOPEX, Jason-1, and Jason-2 Satellites. Mar. Geod. 2012, 35, 238–257. [Google Scholar] [CrossRef]
  48. Reynolds, R.W.; Smith, T.M. Improved global sea surface temperature analyses using optimum interpolation. J. Clim. 1994, 7, 929–948. [Google Scholar] [CrossRef] [Green Version]
  49. Gentemann, C.L.; Wentz, F.J.; Mears, C.A.; Smith, D.K. In situ validation of Tropical Rainfall Measuring Mission microwave sea surface temperatures. J. Geophys. Res. Oceans 2004, 109. [Google Scholar] [CrossRef] [Green Version]
  50. Roemmich, D.H.; Owens, W.B. The ARGO Project: Global Ocean Observations for Understanding and Prediction of Climate Variability. Oceanography 2000, 13, 45–50. [Google Scholar] [CrossRef] [Green Version]
  51. Yang, Y.J.; Chang, M.H.; Hsieh, C.Y.; Chang, H.I.; Jan, S.; Wei, C.L. The role of enhanced velocity shears in rapid ocean cooling during Super Typhoon Nepartak 2016. Nat. Commun. 2019, 10, 1627. [Google Scholar] [CrossRef] [PubMed]
  52. Mellor, G.L. Users Guide for a Three-Dimensional, Primitive Equation, Numerical Ocean Model; Princeton University: Princeton, NJ, USA, 2004. [Google Scholar]
  53. Simos, T.J. Verification of Numerical Models of Lake Ontario Part I. Circulation in Spring and Early Summer. J. Phys. Oceanogr. 1974, 4, 207–523. [Google Scholar]
  54. Madala, R.V.; Piacsek, S.A. A semi-implicit numerical model for baroclinic oceans. J. Comput. Phys. 1977, 23, 167–178. [Google Scholar] [CrossRef]
  55. Zhao, B.; Qiao, F.; Cavaleri, L.; Wang, G.; Bertotti, L.; Liu, L. Sensitivity of typhoon modeling to surface waves and rainfall. J. Geophys. Res. Oceans 2017, 122, 1702–1723. [Google Scholar] [CrossRef]
  56. Oey, L.Y.; Ezer, T.; Wang, D.P.; Fan, S.J.; Yin, X.Q. Loop Current warming by Hurricane Wilma. Geophys. Res. Lett. 2006, 33, 153–172. [Google Scholar] [CrossRef] [Green Version]
  57. Da Silva, A.M.; Young, C.C.; Levitus, S. Atlas of surface marine data 1994, Vol. 1: Algorithms and procedures. Noaa Atlas Nesdis 1994, 6, 20910–23282. [Google Scholar]
  58. Yuan, Y.; Hua, F.; Pan, Z.; Sun, L. LAGFD-WAM numerical wave model-Ⅰ. Basic physical model. Acta Oceanol. Sin. 1991, 10, 483–488. [Google Scholar]
  59. Yuan, Y.; Hua, F.; Pan, Z.; Sun, L. LAGFD-WAM numerical wave model-II. Characteristics inlaid scheme and its application. Acta Oceanol. Sin. 1992, 11, 18–28. [Google Scholar]
  60. Yang, Y.; Qiao, F.L.; Zhao, W.; Teng, Y.; Yuan, Y.L. MASNUM ocean wave numerical model in spherical coordinates and its application. Acta Oceanol. Sin. 2005, 27, 1–7. [Google Scholar]
  61. Phillips, O.M. Spectral and statistical properties of the equilibrium range in wind-generated gravity waves. J. Fluid Mech. 1985, 156, 505–531. [Google Scholar] [CrossRef]
  62. Miles, J.W. On the generation of surface waves by shear flows. J. Fluid Mech. 1959, 7, 185–204. [Google Scholar]
  63. Wu, J. Wind-stress coefficients over sea surface from breeze to hurricane. J. Geophys. Res. Oceans 1982, 87, 9704–9706. [Google Scholar] [CrossRef]
  64. Hasselmann, K.; Barnett, T.P.; Bouws, E.; Carlson, H.; Cartwright, D.E.; Enke, K.; Ewing, J.A.; Gienapp, H.; Hasselmann, D.E.; Kruseman, P.; et al. Measurements of wind-wave growth and swell decay during the Joint North Sea Wave Project (JONSWAP). Ergänzungsheft. 1973, 12, 8–12. [Google Scholar]
  65. Fujita, T. Pressure distribution within typhoon. Geophys. Mag. 1952, 23, 437–451. [Google Scholar]
  66. Takahashi, K. Distribution of pressure and wind in a typhoon. J. Meteor. Soc. Jpn. 1939, 17, 417–421. [Google Scholar]
  67. Ueno, T. Numerical computations of the storm surges in tosa bay. J. Oceanogr. 1981, 37, 61–73. [Google Scholar] [CrossRef]
  68. Zheng, J.H.; Sang, S.; Wang, J.C.; Zhou, C.Y.; Zhao, H.J. Numerical simulation of typhoon-induced storm surge along Jiangsu coast, Part I: Analysis of tropical cyclone. Water Sci. Eng. 2017, 10, 2–7. [Google Scholar] [CrossRef]
  69. Zhang, X.; Chu, P.C.; Li, W.; Liu, C.; Zhang, L.; Shao, C.; Zhang, X.; Chao, G.; Zhao, Y. Impact of Langmuir turbulence on the thermal response of the ocean surface mixed layer to Super Typhoon Haitang (2005). J. Phys. Oceanogr. 2018, 48, 1651–1674. [Google Scholar] [CrossRef]
  70. Sanford, T.B.; Price, J.F.; Girton, J.B. Upper-Ocean Response to Hurricane Frances (2004) Observed by Profiling EM-APEX Floats. J. Phys. Oceanogr. 2011, 41, 1041–1056. [Google Scholar] [CrossRef]
  71. Fan, Y.; Ginis, I.; Hara, T.; Wright, W.C.; Walsh, E.J. Numerical Simulations and Observations of Surface Wave Fields under an Extreme Tropical Cyclone. J. Phys. Oceanogr. 2012, 39, 2097–2116. [Google Scholar] [CrossRef]
  72. Price, J.F. Upper Ocean Response to a Hurricane. J. Phys. Oceanogr. 1981, 11, 153–175. [Google Scholar] [CrossRef] [Green Version]
  73. Moon, I.J.; Ginis, I.; Hara, T.; Tolman, H.L.; Wright, C.W.; Walsh, E.J. Numerical Simulation of Sea Surface Directional Wave Spectra under Hurricane Wind Forcing. J. Phys. Oceanogr. 2003, 33, 1680–1706. [Google Scholar] [CrossRef] [Green Version]
  74. Stoney, L.; Walsh, K.; Babanin, A.V.; Ghantous, M.; Govekar, P.; Young, I. Simulated ocean response to tropical cyclones: The effect of a novel parameterization of mixing from unbroken surface waves. J. Adv. Modeling Earth Syst. 2017, 9, 759–780. [Google Scholar] [CrossRef]
  75. Zhang, Y.; Haiming, X.U.; Qiao, F.; Dong, C. Seasonal variation of the global mixed layer depth: Comparison between Argo data and FIO-ESM. Front. Earth Sci. 2018, 12, 24–36. [Google Scholar] [CrossRef]
  76. Yablonsky, R.M.; Ginis, I.; Thomas, B.; Tallapragada, V.; Sheinin, D.; Bernardet, L. Description and Analysis of the Ocean Component of NOAA’s Operational Hurricane Weather Research and Forecasting Model (HWRF). J. Atmos. Ocean. Technol. 2015, 32, 144–163. [Google Scholar] [CrossRef] [Green Version]
  77. Pan, J.; Sun, Y. Estimate of Ocean Mixed Layer Deepening after a Typhoon Passage over the South China Sea by Using Satellite Data. J. Phys. Oceanogr. 2012, 43, 498–506. [Google Scholar] [CrossRef]
  78. Wu, Q.; Chen, D. Typhoon-induced variability of the oceanic surface mixed layer observed by argo floats in the western north Pacific Ocean. Atmos. Ocean. 2012, 50, 4–14. [Google Scholar] [CrossRef] [Green Version]
  79. Price, J.F.; Sanford, T.B.; Forristall, G.Z. Forced Stage Response to a Moving Hurricane. J. Phys. Oceanogr. 1994, 24, 233–260. [Google Scholar] [CrossRef]
  80. Cuypers, Y.; Vaillant, X.L.; Bouruet-Aubertot, P.; Vialard, J.; Mcphaden, M.J. Tropical storm-induced near-inertial internal waves during the Cirene experiment: Energy fluxes and impact on vertical mixing. J. Geophys. Res. Oceans 2013, 118, 358–380. [Google Scholar] [CrossRef] [Green Version]
  81. Furuichi, N.; Hibiya, T.; Niwa, Y. Model-predicted distribution of wind-induced internal wave energy in the world’s oceans. J. Geophys. Res. 2008, 113. [Google Scholar] [CrossRef]
  82. Huang, R.X. Decadal variability of wind-energy input to the world ocean. Deep Sea Res. Part II Top. Stud. Oceanogr. 2006, 53, 31–41. [Google Scholar] [CrossRef] [Green Version]
  83. Price, J.F. Internal Wave Wake of a Moving Storm. Part I. Scales, Energy Budget and Observations. J. Phys. Oceanogr. 1983, 13, 949–965. [Google Scholar] [CrossRef] [Green Version]
  84. Alford, M.H.; Cronin, M.F.; Klymak, J.M. Annual Cycle and Depth Penetration of Wind-Generated Near-Inertial Internal Waves at Ocean Station Papa in the Northeast Pacific. J. Phys. Oceanogr. 2012, 42, 889–909. [Google Scholar] [CrossRef]
  85. Brooks, D.A. The Wake of Hurricane Allen in the Western Gulf of Mexico. J. Phys. Oceanogr. 1983, 13, 117–129. [Google Scholar] [CrossRef] [Green Version]
  86. Chalikov, D.V.; Makin, V.K. Models of the wave boundary layer. Bound. Layer Meteorol. 1991, 56, 83–99. [Google Scholar] [CrossRef]
  87. Janssen, P.A.E.M. Wave-Induced Stress and the Drag of Air Flow over Sea Waves. J. Phys. Oceanogr. 1989, 19, 745–754. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Model computational domain and bathymetry. The red, green, and purple lines show the 50 m, 200 m, and 3000 m isobaths, respectively. The black dashed box indicates the enlarged region. (b) The best track of TC Nepartak. The colors indicate the stage of TC development. TD: tropical depression; TS: tropical storm; TY: typhoon; ST: super typhoon. Cross represents the Argo float 2902666. Another 69 recording positions are marked with triangles. The yellow and blue pentacle stars show the location of the NTU1 and NTU2 buoys, respectively.
Figure 1. (a) Model computational domain and bathymetry. The red, green, and purple lines show the 50 m, 200 m, and 3000 m isobaths, respectively. The black dashed box indicates the enlarged region. (b) The best track of TC Nepartak. The colors indicate the stage of TC development. TD: tropical depression; TS: tropical storm; TY: typhoon; ST: super typhoon. Cross represents the Argo float 2902666. Another 69 recording positions are marked with triangles. The yellow and blue pentacle stars show the location of the NTU1 and NTU2 buoys, respectively.
Remotesensing 12 02808 g001
Figure 2. Scatter diagrams (a) the color means the number of data in logarithmic form and error probability distribution (b) of the SWH for the year 2016 in the computational domain.
Figure 2. Scatter diagrams (a) the color means the number of data in logarithmic form and error probability distribution (b) of the SWH for the year 2016 in the computational domain.
Remotesensing 12 02808 g002
Figure 3. Box and whisker plots of the SST differences between MW_IR and three experiment results.
Figure 3. Box and whisker plots of the SST differences between MW_IR and three experiment results.
Remotesensing 12 02808 g003
Figure 4. (a) Simulated significant wave height at 0000 UTC 04 July 2016. The white line is the TC track and the gray line is the tandem track of the Jason-2 satellite. (b) significant wave height (SWH) comparison between observation and model during 0424 and 0427 UTC 04 July.
Figure 4. (a) Simulated significant wave height at 0000 UTC 04 July 2016. The white line is the TC track and the gray line is the tandem track of the Jason-2 satellite. (b) significant wave height (SWH) comparison between observation and model during 0424 and 0427 UTC 04 July.
Remotesensing 12 02808 g004
Figure 5. Distribution of SST variations of the three experiments and observations on 11 July ((a) EX1; (b) EX2; (c) EX3; (d) MW_IR). The black line indicates the TC path. The −2 °C contour is marked. The black dashed box indicates the region of Figure 12. The black dashed line shows the profile perpendicular to the track.
Figure 5. Distribution of SST variations of the three experiments and observations on 11 July ((a) EX1; (b) EX2; (c) EX3; (d) MW_IR). The black line indicates the TC path. The −2 °C contour is marked. The black dashed box indicates the region of Figure 12. The black dashed line shows the profile perpendicular to the track.
Remotesensing 12 02808 g005
Figure 6. Comparisons between the model data and the NTU1 data. (a) Wind speed and pressure. (b) Upper 25 m-averaged temperature, the time axis indicates the model output time. (c) Temperature profiles at the corresponding moment ticked in b.
Figure 6. Comparisons between the model data and the NTU1 data. (a) Wind speed and pressure. (b) Upper 25 m-averaged temperature, the time axis indicates the model output time. (c) Temperature profiles at the corresponding moment ticked in b.
Remotesensing 12 02808 g006
Figure 7. Comparisons between the model data and the NTU2 data. (a): Wind speed and pressure. (b): Upper 25 m-averaged temperature, the time axis indicates the model output time. (c): Temperature profiles at the corresponding moment ticked in b.
Figure 7. Comparisons between the model data and the NTU2 data. (a): Wind speed and pressure. (b): Upper 25 m-averaged temperature, the time axis indicates the model output time. (c): Temperature profiles at the corresponding moment ticked in b.
Remotesensing 12 02808 g007
Figure 8. Comparisons between modeled and observed surface heat flux at buoy station NUT1 and NUT2.
Figure 8. Comparisons between modeled and observed surface heat flux at buoy station NUT1 and NUT2.
Remotesensing 12 02808 g008
Figure 9. Plots of MLD (a) and MLT (b) derived from Argo data and numerical experiments at different times, during the processes of TC Nepartak. See Table 1 for the specification of the different experiments. Time is labeled when observation data is available.
Figure 9. Plots of MLD (a) and MLT (b) derived from Argo data and numerical experiments at different times, during the processes of TC Nepartak. See Table 1 for the specification of the different experiments. Time is labeled when observation data is available.
Remotesensing 12 02808 g009
Figure 10. Scatterplots of MLD (a) and MLT (b) between Argo float observations and numerical experiments with the related linear best-fit lines. Observations are from floats within the TC-affected area during TC. The observed locations are shown in Figure 1b and the number of data is 69.
Figure 10. Scatterplots of MLD (a) and MLT (b) between Argo float observations and numerical experiments with the related linear best-fit lines. Observations are from floats within the TC-affected area during TC. The observed locations are shown in Figure 1b and the number of data is 69.
Remotesensing 12 02808 g010
Figure 11. Profiles perpendicular to the track on 8 July 2016. Modeled temperature results for (a) EX2 and (b) EX3. The temperature difference between (c) EX2 and EX1, and between (d) EX3 and EX1. The black dashed line indicates the track position.
Figure 11. Profiles perpendicular to the track on 8 July 2016. Modeled temperature results for (a) EX2 and (b) EX3. The temperature difference between (c) EX2 and EX1, and between (d) EX3 and EX1. The black dashed line indicates the track position.
Remotesensing 12 02808 g011
Figure 12. The distribution of surface BSMT in the dashed black box of Figure 5b. The black line indicates the TC track. (a) 9 July 2016; (b) 10 July 2016; (c) 11 July 2016.
Figure 12. The distribution of surface BSMT in the dashed black box of Figure 5b. The black line indicates the TC track. (a) 9 July 2016; (b) 10 July 2016; (c) 11 July 2016.
Remotesensing 12 02808 g012
Figure 13. Time series of the zonal velocity of upper 1000 m at float 2902666 obtained from the present three numerical experiments. The black dashed line shows the time TC is closest to this point.
Figure 13. Time series of the zonal velocity of upper 1000 m at float 2902666 obtained from the present three numerical experiments. The black dashed line shows the time TC is closest to this point.
Remotesensing 12 02808 g013
Figure 14. Horizontal distributions of the surface temperature Reynold flux of the three experiments at 0600 UTC 6 Jul 2005. (a,d,h) show the total temperature Reynold flux of the EX1, EX2, and EX3. (b,e) show the turbulence-induced temperature Reynold flux of the EX2 and EX3. (c,f) show the WITM-induced temperature Reynold flux of the EX2 and EX3. (g) shows the WTFR-induced temperature Reynold flux of the EX3. The vector shows the model wind data.
Figure 14. Horizontal distributions of the surface temperature Reynold flux of the three experiments at 0600 UTC 6 Jul 2005. (a,d,h) show the total temperature Reynold flux of the EX1, EX2, and EX3. (b,e) show the turbulence-induced temperature Reynold flux of the EX2 and EX3. (c,f) show the WITM-induced temperature Reynold flux of the EX2 and EX3. (g) shows the WTFR-induced temperature Reynold flux of the EX3. The vector shows the model wind data.
Remotesensing 12 02808 g014
Table 1. Experimental design.
Table 1. Experimental design.
WITMWTFR
EX1NoNo
EX2YesNo
EX3YesYes
Table 2. Comparisons of the SWH and SST between the observations and model results.
Table 2. Comparisons of the SWH and SST between the observations and model results.
MERMSErER2
SWH (m)−0.130.470.950.770.79
SST (°C)EX10.440.520.720.760.78
EX20.370.470.820.820.83
EX30.360.460.840.840.89
Table 3. Comparisons of the T25m between the buoy observations and numerical experiments during the TC.
Table 3. Comparisons of the T25m between the buoy observations and numerical experiments during the TC.
ME (°C)RSME (°C)rER2
NTU1EX11.841.920.53−1.080.89
EX21.361.430.63−0.710.95
EX30.570.680.960.610.97
NTU2EX11.551.670.35−1.570.75
EX20.900.990.56−0.260.90
EX3-0.240.451.200.740.91
Table 4. Comparisons of the MLD and MLT between the Argo and numerical experiments during the TC.
Table 4. Comparisons of the MLD and MLT between the Argo and numerical experiments during the TC.
MERSMErER2
MLD (m)EX1−8.9713.240.510.320.67
EX2−2.866.840.720.820.88
EX3-2.605.130.780.900.96
MLT (°C)EX10.310.590.870.490.50
EX20.250.621.090.520.49
EX30.120.571.160.750.56

Share and Cite

MDPI and ACS Style

Yu, C.; Yang, Y.; Yin, X.; Sun, M.; Shi, Y. Impact of Enhanced Wave-Induced Mixing on the Ocean Upper Mixed Layer during Typhoon Nepartak in a Regional Model of the Northwest Pacific Ocean. Remote Sens. 2020, 12, 2808. https://doi.org/10.3390/rs12172808

AMA Style

Yu C, Yang Y, Yin X, Sun M, Shi Y. Impact of Enhanced Wave-Induced Mixing on the Ocean Upper Mixed Layer during Typhoon Nepartak in a Regional Model of the Northwest Pacific Ocean. Remote Sensing. 2020; 12(17):2808. https://doi.org/10.3390/rs12172808

Chicago/Turabian Style

Yu, Chengcheng, Yongzeng Yang, Xunqiang Yin, Meng Sun, and Yongfang Shi. 2020. "Impact of Enhanced Wave-Induced Mixing on the Ocean Upper Mixed Layer during Typhoon Nepartak in a Regional Model of the Northwest Pacific Ocean" Remote Sensing 12, no. 17: 2808. https://doi.org/10.3390/rs12172808

APA Style

Yu, C., Yang, Y., Yin, X., Sun, M., & Shi, Y. (2020). Impact of Enhanced Wave-Induced Mixing on the Ocean Upper Mixed Layer during Typhoon Nepartak in a Regional Model of the Northwest Pacific Ocean. Remote Sensing, 12(17), 2808. https://doi.org/10.3390/rs12172808

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop