Xuewei Li, Dongliang Zhao, Zhongshui Zou. A numerical study of Stokes drift and thermal effects on the oceanic mixed layer[J]. Acta Oceanologica Sinica, 2020, 39(5): 35-45. doi: 10.1007/s13131-019-1448-9
Citation: Xuewei Li, Dongliang Zhao, Zhongshui Zou. A numerical study of Stokes drift and thermal effects on the oceanic mixed layer[J]. Acta Oceanologica Sinica, 2020, 39(5): 35-45. doi: 10.1007/s13131-019-1448-9

A numerical study of Stokes drift and thermal effects on the oceanic mixed layer

doi: 10.1007/s13131-019-1448-9
Funds:  The National Natural Science Foundation of China under contract Nos 41876010 and 41276015; the Public Science and Technology Research Funds Projects of Ocean under contrct No. 201505007; the Joint Project for the National Oceanographic Center by the NSFC and Shandong Government under contract No. U1406402; the National Natural Science Foundation of China under contract No. 41806028.
More Information
  • This study explores the influence of Stokes drift and the thermal effects on the upper ocean bias which occurs in the summer with overestimated sea surface temperature (SST) and shallower mixed layer depth (MLD) using Mellor-Yamada turbulence closure scheme. The upper ocean thermal structures through Princeton ocean model are examined by experiments in the cases of idealized forcing and real observational situation. The results suggest that Stokes drift can generally enhance turbulence kinetic energy and deepen MLD either in summer or in winter. This effect will improve the simulation results in summer, but it will lead to much deeper MLD in winter compared to observational data. It is found that MLD can be correctly simulated by combining Stokes drift and the thermal effects of the cool skin layer and diurnal warm layer on the upper mixing layer. In the case of high shortwave radiation and weak wind speed, which usually occurs in summer, the heat absorbed from sun is blocked in the warm layer and prevented from being transferred downwards. As a result, the thermal effects in summer nearly has no influence on dynamic effect of Stokes drift that leads to deepening MLD. However, when the stratification is weak in winter, the thermal effects will counteract the dynamic effect of Stokes drift through enhancing the strength of stratification and suppress mixing impact. Therefore, the dynamic and thermal effects should be considered simultaneously in order to correctly simulate upper ocean thermal structures in both summer and winter.
  • It is vital for ocean models to correctly reproduce the mixed layer depth (MLD) and sea surface temperature (SST), which are crucial for the air-sea interaction. SST is the lower thermal boundary condition and is usually determined by direct observation at a depth of 2–10 m under sea surface conducted by bucket thermometers of the buoy or the ship’s cooling water intake.

    On the other hand, it is well known that the temperatures at the sea surface are typically a few-tenths degrees Celsius cooler than the temperatures some tens of centimeters below (Fairall et al., 1996a). This layer is referred to as the cool skin layer, which is recognized as an important feature of the ocean viscous layer as a result of new satellite remote sensing methodologies emerging for air-sea flux estimates (Tu and Tsuang, 2005). The warm layer is a region in the upper few meters of the ocean where solar radiation has caused significant warming relative to the deeper mixed layer temperature. Warm layers occur during the day when temperature stratification caused by absorption of the solar radiation is sufficiently strong to suppress shear-induced mixing and can be of the order of 3 K, especially in summer (Halpern and Reed, 1976; Stramma et al., 1986; Fairall et al., 1996a; Soloviev and Lukas, 1997). Fairall et al. (1996a) indicated the peak afternoon warming is about 3.8°C at depth of 0.7 m for a clear day with 10-m wind speeds of 1 m/s. Ward and Donelan (2006) showed strong diurnal warming up to 4.6°C within the upper few meters based on the observation in the Gulf of California. Mattews et al. (2014) also found that the anomaly due to diurnal warm layer can reach to 0.8°C in the afternoon with a daily mean of 0.2°C. In practice, the depth measured SST is well below the cool skin layer and warm layer (Fairall et al., 1996a; Mellor and Blumberg, 2004; Wurl et al., 2019).

    The turbulence closure scheme of Mellor-Yamada (Mellor and Yamada, 1982, hereinafter M-Y) is widely used in atmospheric and oceanic models. However, many studies indicated that MLD simulated by this kind of ocean models in summer is shallower although it is roughly consistent with observations in winter (Martin, 1985; Kantha and Clayson, 1994). Compared with the observations at ocean weather Station Papa, Martin (1985) showed that M-Y scheme underestimates the MLD and overestimates SST in summer. It is shown that insufficient mixing in the upper ocean during summer is a common problem of ocean circulation models, and there are no exceptions for models using M-Y scheme. This problem could not be corrected by adjusting the external air-sea flux parameterizations. For example, a reduction of the surface heat flux will reduce SST in summer, but the simulated MLD will become too high. Many studies have investigated this problem by the inclusion of surface waves, wave breaking, Stokes drift or Langmuir turbulence (Kantha and Clayson, 1994; Huang et al., 2011; Wu et al., 2015; Chen et al., 2018).

    Traditionally, most oceanic mixed layer models ignored the fact that the air-sea interface is a non-rigid, mobile surface with sustaining gravity wave motions. The upper ocean is influenced by all the attendant complex dynamics, including the turbulent kinetic energy injected near the surface by breaking waves and the kinetic energy input from Langmuir circulations (Kantha and Clayson, 1994; Min and Noh, 2004; Huang et al., 2011). The latter is associated with the Stokes drift, which is considered to be an efficient way to enhance both the dissipation rate and the turbulence energy in the upper ocean mixing. Craig and Banner (1994) ascribed the result of insufficient turbulence energy induced by ocean waves, especially the wave breaking, is the warm bias in summertime SSTs in the models. Mellor (2001) introduced a critical Richardson number to directly reduce the turbulence dissipation. As a result, the vertical mixing effect is enhanced and the simulated SSTs decreased in summer. Mellor and Blumberg (2004) further included wave breaking energetics to M-Y scheme to reduce the warm bias in the summer. Recently, many studies suggested that wave breaking-induced turbulence decays rapidly with distance from the surface and hence the overall effects on the mixed layer can be neglected (Kantha et al., 2014; Wu et al., 2015).

    McWilliams et al. (1997) demonstrated the large potential impact of Langmuir turbulence on mixing in the oceanic mixed layer and Stokes-Coriolis force on Ekman turning in the upper ocean. Kantha and Clayson (2004) confirmed the importance of this additional source of energy to turbulent mixing in the oceanic mixed layer, which arises due to the action of Reynolds stresses on the shear of the wave-induced Stokes drift current. Ardhuin and Jenkins (2006) further suggested that the vertical shear of the wave-induced Stokes drift current behaves very similar to the vertical shear of the mean current, as far as the extraction of energy by the Reynolds stress is concerned. While the energy input to turbulence from Langmuir cells elevates the turbulent kinetic energy and mixing throughout the mixed layer, therefore, more effective in deepening the mixed layer (Kantha and Clayson, 2004; Min and Noh, 2004; Huang et al., 2011). It is now widely recognized that a mixed layer model must include Stokes production of turbulent kinetic energy (TKE) for a realistic simulation of properties in the mixed layer (Teixeira and Belcher, 2002; Kumar and Feddersen, 2017). For example, Smith et al. (2013) recognized the importance of including the Stokes drift in their ocean-wave coupled model.

    In order to reduce the warm bias in the summertime, all of previous studies focused on how to enhance the vertical mixing effect through increasing TKE in the models. In this way, it is certain that it can deepen MLD and reduce the warm bias in summer. However, it also would more significantly deepen MLD in winter due to the high wind speed and low insolation. Therefore, it should be considered the other factors that reflect the different characteristics between summer and winter.

    As mentioned above, in addition to the dynamic influence on the oceanic mixed layer, it is believed that the thermal effects on stratification are also important for vertical mixing. It has been known for decades that the upper few micrometers of the ocean are generally cooler by several tenths of degree compared to the well-mixed underlying water masses. This sea surface microlayer is called cool skin layer, in which the heat transfer is primarily conducted by molecular diffusion since turbulent motion is suppressed (Saunders, 1967; Katsaros et al., 1977; Paulson and Simpson, 1981; Wu, 1985; (Fairall et al., 1996b; Ward and Donelan, 2006; Alappattu et al., 2017; Wong and Minnett, 2018; Wurl et al., 2019). The skin SST (hereinafter Tskin) is defined as the radiometric temperature measured across a very small depth of approximately 20 μm. The subskin SST (hereinafter Tsubskin) represents the temperature at the base of the thermal skin layer (~1 000 μm according to the definition in Donlon et al., 2002). According to Wurl et al. (2019), the temperature difference between the skin temperature and the underlying water at 1-m depth, ΔT = TskinT1 m was generally negative with the median = −0.20°C, in which the positive values only occur at times with strong solar radiation (> 600 W/m2) and low wind speeds (< 2 m/s).

    The calculations of bulk fluxes are based on the empirical formula between the turbulent fluxes and the air-sea contrasts of wind, humidity and temperature. Logically, the proper temperature for the determination of sensible and latent turbulent heat fluxes as well as the upwelling longwave emission is Tskin, which is the temperature of the water which is in direct contact with the atmosphere ((Fairall et al., 1996b). However, the radiometric measurements during field campaigns are not always available, thus a model is required to calculate Tskin from oceanographic and meteorological measurements ( Tu and Tsuang, 2005).

    Therefore, it is important to include the skin temperature for accurate calculation of latent, sensible, and net longwave heat fluxes, which needs to decouple Tskin from the conventional in situ bulk SST (hereinafter Tbulk) measured at a depth of a few meters. In general, there is a large difference between Tskin and Tbulk, which will significantly affect the calculation of air-sea fluxes. Accordingly, the cool skin and warm layer correction have been introduced into the standard bulk flux algorithm (COARE 3.0, (Fairall et al., 1996b). In order to correctly simulate the oceanic mixed layer, the vertical resolution in the model must be high enough to capture the thermal effects of the warm layer ((Fairall et al., 1996b; Tseng et al., 2015).

    In this study, it will be investigated the dynamic effect of Stokes drift and thermal effects of the cool-skin and warm layer on the oceanic mixed layer. Section 2 provides a brief description of the model modification due to Stokes drift. The results for idealized forcing experiments and Ocean Weather Station Papa are presented in Section 3 and Section 4, respectively. Finally, Section 5 gives a summary and discussion.

    As shown by Sullivan et al. (2007), the wave influences are shown to be profound on both the mean current profile and turbulent statistics compared to a simulation without these wave influences and forced by an equivalent mean surface stress. The elevation spectrum of wind waves implies a mean Stokes drift profile uS(z), which can be expressed as (Kenyon, 1969)

    $${u_{\rm s}}\left(z \right) = \frac{{16{{\text π} ^2}}}{g}\mathop \int \nolimits_0^{2{\text π} } \mathop \int \nolimits_0^\infty {f\,^3}S\left({f,\theta } \right)\exp \left({\frac{{8{\text π} {f\,^2}}}{g}z} \right)\cos \left({\theta - {\theta _{\rm w} }} \right){\rm{d}}f{\rm{d}}\theta,$$ (1)

    where S(f, θ) is the directional frequency spectrum of wind waves, f is frequency, g is the acceleration of gravity, z is the vertical coordination, and θ, θw are the wave and wind directions, respectively. In general, the development of wind waves depends on wave age or fetch. For a given mean wind speed and sufficiently large fetch, the wave field comes into a fully-developed wind waves, which is independent of wave age or fetch. For simplicity, we just consider the situation of fully-developed wind waves. With wave dispersion relationship in deep water, the Stokes drift is reduced to (Kantha and Clayson, 2004)

    $${u_{\rm s}}\left(z \right) = \frac{{2{{\text π} ^3}H_{\rm s}^2}}{{gT_{\rm s}^3}}{\rm{exp}}\left({\frac{{8{{\text π} ^2}}}{{gT_{\rm s}^2}}z} \right),$$ (2)

    where Hs and Ts are significant wave height and wave period of wind waves, respectively, which correspond to the fully-developed wind waves. They can be estimated by the empirical growth relations proposed by many studies. In this study, we choose the relationships suggested by Wilson (1965),

    $$\begin{split} & {H_{\rm s}} = 0.030\;6U_{10}^2,\\ & {T_{\rm s}} = 0.878{U_{10}}. \end{split}$$ (3)

    Substituted Eq. (3) into Eq. (2), the Stokes drift can be expressed as

    $${u_{\rm s}}\left(z \right) = 0.008\;8{U_{10}}{\rm{exp}}\left({10.45U_{10}^{ - 2}z} \right).$$ (4)

    This result is a little different from that of Harcourt and D’Asaro (2008), which is

    $${u_{\rm s}}\left(z \right) = 0.017\;5{U_{10}}{\rm{exp}}\left({45.44U_{10}^{ - 2}z} \right).$$ (5)

    Equation (5) has a greater magnitude of Stokes drift at sea surface than Eq. (4), but it decreases more quickly with depth z than Eq. (4).

    It has been shown that the influences of Stokes drifts are profound on both the mean current profile and turbulent statistics compared to a simulation without these wave influences. Following Sullivan et al. (2007), we have added the Stokes drifts into the mean current profile. In this way, the dynamic effect of Stokes drift on deepening oceanic mixed layer is included in our simulation. The Princeton ocean model (POM) is used to conduct a series of simulations with different model settings.

    POM is a three-dimensional, primitive equation, numerical ocean model. We set the model domain as 30°S–60°N and 135°–155°W with a regular 0.5°×0.5° longitude-latitude grid. The bottom topography is obtained from ETOPO5, and the maximum water depth is set as 2 000 m. The internal time step used was 6 min.

    It is evident that fine vertical resolution in the model is necessary to reproduce the warm layer (Fairall et al., 1996b; Tseng et al., 2015) . In order to show the influence of various resolutions in POM, we adopt two different vertical resolutions in the simulations. In experiment 1 (hereafter referred as EXP1), the vertical resolution from sea surface to depth of 100 m is set as 2 m, and 13 m from 100 m to 2 000 m. In experiment 2 (hereafter referred as EXP2), except the vertical resolution from sea surface to depth of 10 m is set as 1 m, the others are same as in EXP1. The additional level is placed at 0.5 m to capture the upper ocean warm layer. With the modification of Stokes drift, EXP1 and EXP2 are referred to as “EXP1-S” and “EXP2-S”, respectively (Table 1). We used COARE 3.0 algorithm to calculate surface latent heat, sensible heat flux and ocean skin temperature. The experiments were tested for two types of idealized forcing in order to compare the MLD response to different forcing conditions. The forcing was chosen to represent two basic regimes: the surface heat input which tends to increase stratification with the smaller MLD, and the momentum flux input which tends to destroy the stratification with the greater MLD. The solar radiation fluxes are assumed as a sine function with the amplitudes 200 W/m2 and 800 W/m2, which represent the situations in winter and summer, respectively. The wind speeds U10 vary from 0 to 20 m/s at interval of 1 m/s.

    Table  1.  Summery of the experiments used in this study
    ExperimentVertical resolutionStokes driftsIdealized forcingOWS Papa
    Wind speed
    /m·s–1
    Amplitudes of the solar radiation
    /W·m–2
    Wind speed
    /m·s–1
    Solar radiation
    /W·m–2
    EXP1lowno0–20200/800observedobserved
    EXP2highno0–20200/800observedobserved
    EXP1-Slowyes0–20200/800observedobserved
    EXP2-Shighyes0–20200/800observedobserved
     | Show Table
    DownLoad: CSV

    Figure 1 is shown the comparison between EXP1 and EXP1-S, in which are both with low vertical resolution. It can be seen that SST will decrease when the Stokes drift is considered (Figs 1a, f). The changes in SST are attributed to the enhancement of turbulence kinetic energy. The influences of the Stokes drift are sensitive to the magnitude of wind speed and the property of the stratification. In low wind speeds, the turbulence kinetic energy caused by the Stokes drift is limited to a few meters above the upper level. As the wind speed increases, turbulence kinetic energy extends into much greater depth. The influence of Stokes drift is also dependent on the strength of the mixed layer stratification in the ocean. If the mixed layer is weakly stratified, the Stokes drift can affect deeper. With the same constant wind speed forcing, the turbulence kinetic energy caused by the Stokes drift extends into deeper when the amplitudes of the solar radiation decreases from 800 to 200 W/m2 (Figs 1d, i). The surface heating weakens mixing process mainly because it enhances the stratification of the ocean. In case of low solar radiation, the stratification is weak and the enhanced turbulent energy can propagate effectively downwards, which means that the effect of Stokes drift is deeper in winter than in summer.

    Figure  1.  The difference between experiment EXP1-S and EXP1 in the vertical structure of simulated. a, f. Temperature profile (T); b, c, g, h. horizontal velocities (U and V); d, i. turbulence kinetic energy (TKE, q2); and e, j. vertical kinematic viscosity (KM). The colors represent different wind speed, respectively. The top row shows the result forcing by high shortwave radiation 800 W/m2 (a–e), the other row shows that forcing by low short wave radiation 200 W/m2 (f–j).

    Figure 2 is illustrated the thermal effects of the skin layer and warm layer on the upper ocean during 24 h from the results of EXP2. The temperature of the cool skin layer is usually lower than the sublayer at a depth of 0.5 m. The solar radiation has caused significant warming near the surface and propagates download as it intensifies with increasing solar intensity. Under the light wind conditions, solar heating produces strong warming within a very thin surface layer. The warming peak occurs at noon, and exists until sunset. In the absence of mechanical mixing, the temperature due to solar heating at 2 m-depth can be increased by 2°C (Fig. 2). On diurnal and longer time scales, this difference is not very noticeable in the model results. At midnight, the upper ocean temperature is almost consistent. With solar radiation increasing, the warming leads to a stably stratified surface layer. As shown in Figure 3, the temperature profile for wind speeds from 0 to 10 m/s with internal of 1 m/s. It is clear that the depth of the warm layer becomes deeper and its magnitude decreases with the increasing of wind speed.

    Figure  2.  The evolution of sea temperature profile with 2 m/s wind in high vertical resolution model result (EXP2) are shown. The colors represent temperature profile per 1 h, respectively.
    Figure  3.  The evolution of sea temperature profile per 3 h in high vertical resolution model result (EXP2) are shown. The colors represent the control experiments with wind speed from 0 to 10 m/s, respectively.

    In the diurnal warm layer, the temperature declines rapidly within a thin layer. The vertical structure of sea temperature in the upper ocean has a “double-thermocline structure”, which is a stably stratified surface layer. The double-thermocline structure gradually weakens with the increasing wind speed, and it disappears when wind speed is up to 9 m/s (Fig. 3). It also can be seen that the Tskin is lower than the temperature at depth of 0.5 m (Fig. 2). Therefore, there is a “cold-warm-cold” structure, which can hinder heat propagation downward. The structure is accompanied with warm layer in the upper oceanic layer, and we call it “blocking structure”. This blocking structure can hinder the heat from being transferred between the ocean and the atmosphere.

    Figure 4 shows the comparison of different vertical resolution (EXP1, EXP2) without the dynamic effect of Stokes drift. The differences of the upper ocean temperature are generally confined to the upper 10 m of the ocean (Fig. 4). Clearly, there are no the blocking structure and the double-thermocline structure in EXP1. Compared to EXP1, the temperature profile of the upper ocean in EXP2 is more realistically simulated and the warm layer is correctly detected. Furthermore, the coarser the vertical resolution, the slower the temperature response to the surface heating . This is because the thicker layer needs more heat to alter the temperature. In the heating process, the surface heat is significant stored in the intensive diurnal warm layer. In the cooling process, comparing with EXP1, the blocking structure in EXP2 hinders the heat from transferring upwards. Consequently, the temperature in the upper few meters of the EXP2 is warmer than that in EXP1 (Fig. 4). It suggests that a fine vertical resolution is needed in order to capture temperature variations in the upper ocean.

    Figure  4.  The evolution of sea temperature profile in different control model result. The colors represent temperature profile per 2 h, respectively. Line: EXP1; and line with points: EXP2.

    To investigate the influence of Stokes drift and thermal effects on the depth of mixed layer, 0-20 m/s wind speeds are used in four experiments. The results for wind speeds of 7, 10, 12 and 15 m/s are shown in Figs 5 and 6, which correspond to the solar radiation of 800 W/m2 and 200 W/m2, respectively. The former represents the case in summer, and the latter represents the case in winter. With the increasing wind speed, the effect of Stokes drift can substantially affect the whole mixed layer. The enhanced turbulent energy mixes the surface warm water and transfer the heat downward. Therefore, the SST decreases with the deeped MLD (the black thick lines in Fig. 5). In Fig. 6, the diurnal variation amplitude of SST is relatively small with deeper MLD and smaller SST. When the wind speed is relatively weak (7 m/s), the depth of winter mixing layer is tens of meters deeper than that of summer. The influence of Stokes drift decreases gradually with the depth, which makes the influence of Stokes drift on the mixing layer depth is weaker under the low solar radiation than that under the high solar radiation. But with the increasing of wind speed (15 m/s), the difference of mixing depth caused by solar radiation decreases. At the same depth, Stokes drift has enhanced influence on the mixing depth under the condition of low solar radiation. Therefore, high wind speed and low solar radiation can contribute to enhance and deepen the influence of Stokes drift.

    Figure  5.  The evolution of MLD and SST in two days with idealized forcing: 800 W/m2 solar radiation and wind speeds from weak to strong. Solid black line: EXP1; solid red line: EXP2; black line with points: EXP1-S; and red line with points: EXP2-S.
    Figure  6.  The evolution of MLD and SST in two days with idealized forcing: 200 W/m2 solar radiation and wind speeds from weak to strong. Solid black line: EXP1; solid red line: EXP2; black line with points: EXP1-S; and red line with points: EXP2-S.

    Forced by strong solar radiation and low wind speed, the upper layer of ocean shows a significant diurnal variation. Due to the existence of the cool skin layer and diurnal warm layer, the double-thermocline structure clearly exists. The cool skin layer acts as a hindrance to heat transfer between ocean and atmosphere. The cool skin layer and diurnal warm layer together effectively hinder the transfer of solar radiation downward. The heat is captured by the diurnal warm layer. In other words, the heat of the solar radiation is actually accumulated in the blocking structure. In EXP1, without considering the existence of thermal effects, it is supposed that all solar radiation effectively enters the mixing layer, in which the shallower effect of solar radiation on the mixing layer is overestimated, and leads to the smaller MLD and warmer SST. When the wind speed increases or the solar radiation reduces, the stratification is weak. In this situation, the thermal effects tend to shallow the MLD. Due to the strong mixing, the warm layer cannot accumulate the heat into the surface stably. With the weak stratification, the cool skin layer prevents heat from being transferred upward. It depresses the turbulent kinetic energy of the upper mixed layer, and keep the intensity of the stratification. Therefore, the thermal effects can enhance the strength of stratification.When the wind speed increases to 15 m/s, the double-thermocline layer and blocking structure is collapsed due to the strong shear-induced mixing ascribed to the wind (Fig. 5).

    In winter, the stratification in the upper ocean is weak due to the weak solar radiation and the strong wind speed. It seems that the increased mixing of Stokes drift will cause excessive mixing. However, the thermal effects can compensate for the weak stratification caused by the lack of solar radiation. The cool skin layer and diurnal warm layer increase the stratification, which inhibits the enhancement of turbulent energy induced by Stokes drift to some extent. In this way, it refrains from the problem of over-mixing in winter.

    Since 1956, Ocean Station Papa (OSP) in the Gulf of Alaska at 50°N, 145°W has been in the site of open ocean research. The oceanography and meteorological measurements conducted on board by the Canadian weatherships. It is believed that the motion of the sea water around OSP is dominated by vertical physical processes and the advective effect can be neglected. Thus the data at OSP have been widely applied to evaluate the mixed layer models. We conduct the tests at OSP focusing on the seasonal evolution of the upper ocean thermal structure. The configurations are the same as those in Section 3.

    In this section, the model is forced by observational data from January 1, 2012 to December 31, 2012. The horizontal resolution is 0.25° by 0.25°. The bottom topography is obtained from ETOPO5, and the maximum water depth is set as 2 000 m. The internal (external) time step is 30 min (60 s). The initial velocity is set to zero, and the initial temperature and the salinity profiles are adopted the data on January 1, 2012. Each experiment is run for one year, which is long enough to simulate of the mixed layer. The hourly data of the surface forcing are linearly interpolated into the model grid. The latent, sensible heat flux and sea skin temperature have been computed from the observed 10 m wind speed, the air temperature and the relative humidity using the advanced method of COARE3.0 bulk algorithm (Fairall et al., 1996b).

    The observed and model predicted SST at depth of 1 m and MLD at OSP for the year 2012 are shown in Fig. 7. The MLD is defined as the depth at which the temperature becomes 0.2°C less than the SST. The observed MLD shows a drop of about 40 meters during the middle of January and a sudden increase at the beginning of April (Fig. 7b). Generally, the fluctuations of SST and MLD are controlled by surface forcing that is correctly demonstrated by the model. The EXP1 obtained the summer bias which is similar to that of Martin (1985). Although the results of EXP1 agree well with observations in winter, the simulated SST is warmer and MLD is shallower than the observations from July to August (Fig. 7). In contrast, In EXP1-S, the SST is well consisted with the observations. The improvement of SST in summer can also be seen in EXP2-S from Fig. 7a. However, in terms of MLD from January to April, there is a clear discrepancy between the EXP1-S and the observations. Because of the weak stratification and the heavy wind speed in winter, the over mixing can be expected from the result in EXP1-S (Fig. 7).

    Figure  7.  Observed and simulated SST and MLD at Station Papa for 2012. Black line: observation; blue line: EXP1; green line: EXP1-S; yellow line: EXP2; and red line: EXP2-S.

    Table 2 is listed the bias of SST and MLD of EXP1, EXP2, EXP1-S and EXP2-S compared with the observations. From Table 2, it can be seen that the simulated SST and MLD in EXP2-S are closest to the observations with the summer SST bias of only 0.48°C and MLD bias of 6.17 m. The simulated SST is overestimated about 2°C in EXP1, and the corresponding MLD is about 13 m shallower than the observations (Table 2). Figure 8 shows the isotherm depths at OSP for the simulations and the observations. It is clear that the observed 13℃-isotherm is deeper and sharper than EXP1 (Fig. 8a). The EXP1-S tends to intense the turbulence and the vertical mixing. Strong mixing can transport a large amount of heat from the surface layer to the subsurface, resulting in a decrease in the surface temperature and an increase in subsurface temperature. The summer bias is associated with the cold subsurface temperature bias of 4℃ (Fig. 9). There is less warming accumulates above the MLD during the summer and significant warm bias under the MLD in both summer and winter (Fig. 9b). Although it is believed that the enhanced mixing contributes significantly to mixing in the ocean mixed layer, it is important to note that the effect of Stokes is deeper in winter. Excess heat is transferred to the deep ocean by Stokes in winter.

    Table  2.  Difference between model and observed SST and MLD at Papa for 2012
    BiasEXP1EXP2EXP1-SEXP2-S
    Annual SST/°C1.180.960.570.53
    Summer SST/°C1.981.210.640.48
    Winter MLD/m9.038.3628.94 7.65
    Summer MLD/m13.11 12.11 8.056.17
     | Show Table
    DownLoad: CSV
    Figure  8.  Observed and simulated isotherm depth at Station Papa for 2012. a. EXP1, b. EXP2, c. EXP1-S, d. EXP2-S, and e. observation.
    Figure  9.  Simulated temperature bias profile at Station Papa for 2012. a. EXP1, b. EXP1-S, c. EXP2, and d. EXP2-S.

    The most outstanding thermal effect is the enhanced stratification which prevents the heat from being propagated between ocean and atmosphere. While in summer (winter), the thermal effects decrease (increase) the net heat flux absorbed in the ocean. A reduction of the surface net heat flux improves the SST agreement in summer, but the simulated MLD will be too shallow (Fig. 10). The temperature can be simulated properly to include the thermal effects and the Stokes drift in EXP2-S. The shallowing of the MLD might be induced by inadequate mixing and overestimate of the surface heat flux during summer. In summer, the thermal effects accumulate the heat in the surface layer and result to reduce the net surface heat flux. The dynamic effect of Stokes drift enhances the turbulence kinetic energy and transfer the heat from surface to subsurface. These two factors can effectively solve the problem of insufficient mixing in the classic M–Y scheme and correct the summer simulated bias. In winter, the dynamic effect of Stokes drift due to the strong wind speed and weak stratification could mix much greater depths than in summer. However, the overestimated mixing can be effectively compensated by the thermal effects.

    Figure  10.  The difference of the net heat flux released into the atmosphere compared with EXP1. Yellow line: EXP2; green line: EXP1-S; and red line: EXP2-S.

    In order to improve the insufficient mixing in the ocean model in summer, the impact of Stokes drift, cool skin layer and warm layer are considered. Several simulation experiments are conducted through POM. It is shown that the Stokes drift can deepen the MLD and reduce the SST by enhancing the turbulence kinetic energy. In the case of low solar radiation, the effect of Stokes drift on the mixed layer is more significant compared with the case of high solar radiation at same wind speed. It will be improved the summer bias, and lead to a new bias in winter. We suggest that this contradict can be avoided by the inclusion of thermal effects of the cool skin layer and warm layer in the upper ocean.

    The thermal effects of the cool skin layer and diurnal warm layer on the upper mixing layer is highly dependent on the magnitude of solar radiation and wind speed. In the case of low wind speeds and strong solar radiation, the heat accumulates in the diurnal warm layer. The double-thermocline structure and blocking structure can accumulate heat and prevent the heat from being transferred downwards. With the strong stratification, the thermal effects can enhance the influence of Stokes drift. But when the stratification is weak, the thermal effects can enhance the strength of stratification and counteract Stokes drift by suppressing excessive mixing effect. With strong winds and weak solar radiation, the effect of warm layer can be neglected while the cool skin effect on the net heat flux intensifies. Therefore, with the dynamic effect of Stokes drift and the thermal effects of skin layer and the warm layer together, POM can correctly simulate the SST and MLD either in summer or in winter.

  • [1]
    Alappattu D P, Wang Qing, Yamaguchi R, et al. 2017. Warm layer and cool skin corrections for bulk water temperature measurements for air-sea interaction studies. Journal of Geophysical Research: Oceans, 122(8): 6470–6481. doi: 10.1002/2017JC012688
    [2]
    Ardhuin F, Jenkins A D. 2006. On the interaction of surface waves and upper ocean turbulence. Journal of Physical Oceanography, 36(3): 551–557. doi: 10.1175/JPO2862.1
    [3]
    Chen Siyu, Qiao Fangli, Huang Chuanjiang, et al. 2018. Effects of the non-breaking surface wave-induced vertical mixing on winter mixed layer depth in subtropical regions. Journal of Geophysical Research: Oceans, 123(4): 2934–2944. doi: 10.1002/2017JC013038
    [4]
    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
    [5]
    Donlon C J, Minnett P J, Gentemann C, et al. 2002. Toward improved validation of satellite sea surface skin temperature measurements for climate research. Journal of Climate, 15(4): 353–369. doi: 10.1175/1520-0442(2002)015<0353:TIVOSS>2.0.CO;2
    [6]
    Fairall C W, Bradley E F, Godfrey J S, et al. 1996a. Cool-skin and warm-layer effects on sea surface temperature. Journal of Geophysical Research: Oceans, 101(C1): 1295–1308. doi: 10.1029/95JC03190
    [7]
    Fairall C W, Bradley E F, Rogers D P, et al. 1996b. Bulk parameterization of air-sea fluxes for tropical ocean-global atmosphere coupled-ocean atmosphere response experiment. Journal of Geophysical Research: Oceans, 101(C2): 3747–3764. doi: 10.1029/95JC03205
    [8]
    Halpern D, Reed R K. 1976. Heat Budget of the Upper Ocean Under Light Winds. Journal of Physical Oceanography, : doi: 10.1175/1520-0485(1976)006<0972:hbotuo>2.0.co;2
    [9]
    Harcourt R R, D’Asaro E A. 2008. Large-eddy simulation of Langmuir turbulence in pure wind seas. Journal of Physical Oceanography, 38(7): 1542–1562. doi: 10.1175/2007JPO3842.1
    [10]
    Huang Chuanjiang, Qiao Fangli, Song Zhenya, et al. 2011. Improving simulations of the upper ocean by inclusion of surface waves in the Mellor-Yamada turbulence scheme. Journal of Geophysical Research: Oceans, 116(C1): C01007. doi: 10.1029/2010JC006320
    [11]
    Kantha L H, Clayson C A. 1994. An improved mixed layer model for geophysical applications. Journal of Geophysical Research: Oceans, 99(C12): 25235–25266. doi: 10.1029/94JC02257
    [12]
    Kantha L H, Clayson C A. 2004. On the effect of surface gravity waves on mixing in the oceanic mixed layer. Ocean Modelling, 6(2): 101–124. doi: 10.1016/S1463-5003(02)00062-8
    [13]
    Kantha L, Tamura H, Miyazawa Y. 2014. Comment on “wave-turbulence interaction and its induced mixing in the upper ocean” by Huang and Qiao. Journal of Geophysical Research: Oceans, 119(2): 1510–1515. doi: 10.1002/2013JC009318
    [14]
    Katsaros K B, Liu W T, Businger J A, et al. 1977. Heat thermal structure in the interfacial boundary layer measured in an open tank of water in turbulent free convection. Journal of Fluid Mechanics, 83(2): 311–335. doi: 10.1017/S0022112077001219
    [15]
    Kenyon K E. 1969. Stokes drift for random gravity waves. Journal of Geophysical Research, 74(28): 6991–6994. doi: 10.1029/JC074i028p06991
    [16]
    Kumar N, Feddersen F. 2017. The effect of stokes drift and transient rip currents on the inner shelf. Part II: With stratification. Journal of Physical Oceanography, 47(1): 243–260. doi: 10.1175/JPO-D-16-0077.1
    [17]
    Martin P J. 1985. Simulation of the mixed layer at OWS November and Papa with several models. Journal of Geophysical Research: Oceans, 90(C1): 903–916. doi: 10.1029/JC090iC01p00903
    [18]
    Matthews A J, Baranowski D B, Heywood K J, et al. 2014. The surface diurnal warm layer in the Indian ocean during CINDY/DYNAMO. Journal of Climate, Special issue: 9101–9122. doi: 10.1175/JCLI-D-14-00222.1
    [19]
    McWilliams J C, Sullivan P P, Moeng C H. 1997. Langmuir turbulence in the ocean. Journal of Fluid Mechanics, 334(1): 1–30. doi: 10.1017/S0022112096004375
    [20]
    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
    [21]
    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
    [22]
    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
    [23]
    Min H S, Noh Y. 2004. Influence of the surface heating on Langmuir circulation. Journal of Physical Oceanography, 34(12): 2630–2641. doi: 10.1175/JPOJPO-2654.1
    [24]
    Paulson C A, Simpson J J. 1981. The temperature difference across the cool skin of the ocean. Journal of Geophysical Research: Oceans, 86(C11): 11044–11054. doi: 10.1029/jc086ic11p11044
    [25]
    Saunders P M. 1967. The temperature at the ocean-air interface. Journal of the Atmospheric Sciences, 24(3): 269–273. doi: 10.1175/1520-0469(1967)024<0269:ttatoa>2.0.co;2
    [26]
    Smith T A, Chen Sue, Campbell T, et al. 2013. Ocean-wave coupled modeling in COAMPS-TC: A study of Hurricane Ivan (2004). Ocean Modelling, 69: 181–194. doi: 10.1016/j.ocemod.2013.06.003
    [27]
    Soloviev A, Lukas R. 1997. Observation of large diurnal warming events in the near-surface layer of the western equatorial Pacific warm pool. Deep-Sea Research Part I: Oceanographic Research Papers, 44(6): doi: 10.1016/S0967-0637(96)00124-0
    [28]
    Stramma L, Cornillon P, Weller R A, et al. 1986. Large diurnal sea surface temperature variability: Satellite and In Situ Measurements. Journal of Physical Oceanography, : doi: 10.1175/1520-0485(1986)016<0827:ldsstv>2.0.co;2
    [29]
    Sullivan P P, McWilliams J C, Melville W K. 2007. Surface gravity wave effects in the oceanic boundary layer: Large-eddy simulation with vortex force and stochastic breakers. Journal of Fluid Mechanics, 593: 405–452. doi: 10.1017/S002211200700897X
    [30]
    Teixeira M A C, Belcher S E. 2002. On the distortion of turbulence by a progressive surface wave. Journal of Fluid Mechanics, 458: 229–267. doi: 10.1017/S0022112002007838
    [31]
    Tseng W L, Tsuang B J, Keenlyside N S, et al. 2015. Resolving the upper-ocean warm layer improves the simulation of the Madden-Julian oscillation. Climate Dynamics, 44(5–6): 1487–1503. doi: 10.1007/s00382-014-2315-1
    [32]
    Tu C Y, Tsuang B J. 2005. Cool-skin simulation by a one-column ocean model. Geophysical Research Letters, 32(22): L22602. doi: 10.1029/2005GL024252
    [33]
    Ward B, Donelan M A. 2006. Thermometric measurements of the molecular sublayer at the air-water interface. Geophysical Research Letters, 33(7): L07605. doi: 10.1029/2005GL024769
    [34]
    Wilson B W. 1965. Numerical prediction of ocean waves in the North Atlantic for December, 1959. Deutsche Hydrografische Zeitschrift, 18(3): 114–130. doi: 10.1007/BF02333333
    [35]
    Wong E W, Minnett P J. 2018. The response of the ocean thermal skin layer to variations in incident infrared radiation. Journal of Geophysical Research: Oceans, 123(4): 2475–2493. doi: 10.1002/2017JC013351
    [36]
    Wu Jin. 1985. On the cool skin of the ocean. Boundary-Layer Meteorology, 31(2): 203–207. doi: 10.1007/BF00121179
    [37]
    Wu Lichuan, Rutgersson A, Sahlée E. 2015. Upper-ocean mixing due to surface gravity waves. Journal of Geophysical Research: Oceans, 120(12): 8210–8228. doi: 10.1002/2015JC011329
    [38]
    Wurl O, Landing W M, Mustaffa N I H, et al. 2019. The ocean’s skin layer in the tropics. Journal of Geophysical Research: Oceans, 124(1): 59–74. doi: 10.1029/2018JC014021
  • 加载中

Catalog

    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Figures(10)  / Tables(2)

    Article Metrics

    Article views (276) PDF downloads(3) Cited by()
    Proportional views
    Related

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return