1. Introduction
A halocline within the isothermal layer exists within several basins, including the tropical Atlantic and Pacific Oceans, the Indian Ocean, and the Caribbean Sea, creating a barrier layer (BL) (Sprintall and Tomczak 1992). This layer aids in the suppression of turbulent heat flux from the thermocline to the air–sea interface, affecting air–sea fluxes (Lukas and Lindstrom 1991; McPhaden and Foltz 2013; Mignot et al. 2012; Chi et al. 2014). The role of BLs on air–sea exchanges is important to understand since numerous tropical cyclones (TCs) pass through BL regions and air–sea exchanges that influence TC intensity are based on the upper-ocean dynamics that occur.
Literature has shown that the BL impacts sea surface temperature (SST) response during TC wind forcing via reduced mixing efficiency (Wang et al. 2011; Neetu et al. 2012; Balaguru et al. 2012; Grodsky et al. 2012; Vissa et al. 2013; Reul et al. 2014; Androulidakis et al. 2016; Hernandez et al. 2016; Rudzin et al. 2017; Yan et al. 2017). However, the dynamics behind BL erosion during TC wind forcing has not been studied in detail in the available literature. Hernandez et al. (2016) briefly examined the impact of the BL on SST response in context of differing thermal and haline regimes. Yet, their goal was not to examine the upper-ocean dynamics that contribute to SST response. Yan et al. (2017) investigated mixed layer dynamics that leads to BL erosion but only examined these processes using one type of mixed layer model. Hence, an understanding of different processes that influence BL erosion is missing.
Upper-ocean dynamics are important to understand since modest SST differences during a TC have been shown to dramatically influence air–sea heat and moisture transfer (Cione and Uhlhorn 2003; Jaimes et al. 2015, 2016). Hence, if SST response is misrepresented in a BL environment within coupled numerical forecast models, this could lead to errors in air–sea transfer and potentially forecasted TC intensity. To address the gaps from previous literature, several one-dimensional (1D) mixed layer model experiments are analyzed to investigate how the inclusion of salinity in different thermal regimes impacts the SST response. BL erosion time is estimated to assess the resilience of the upper-ocean thermal and haline structure to each mixing scheme. The use of several mixing schemes highlights how individual upper-ocean processes, such as shear-induced mixing and instantaneous wind erosion, influence SST response. The findings are put in context of a BL baroclinic wave speed to understand how a stratified upper ocean influences ocean coupling during TC passage.
2. Data and methods
Six different temperature profiles and two salinity profiles from Rudzin et al. (2017) are selected to highlight mixing sensitivity between differing ocean regimes commonly found in the Caribbean Sea, a basin with frequent TCs and subsiding BL. These temperature and salinity profiles were measured in September 2014 within the eastern Caribbean Sea under ambient atmospheric conditions. Temperature profiles represent a warm-core eddy (WCE) (Figs. 1a,b), the Caribbean Current (CC) background flow (Figs. 1c,d), and the Amazon–Orinoco River plume (PLUME) (Figs. 1e,f). The two salinity profiles are selected to represent a strong (Figs. 1a,c,e; strong S) and weak (Figs. 1b,d,f; weak S) halocline in the isothermal layer to identify how stratification strength affects SST response.
a. Mixing schemes
The 1D mixing schemes used are the Kraus–Turner (KT) (Kraus and Turner 1967), Price–Weller–Pinkel (PWP) (Price et al. 1986), and Pollard–Rhines–Thompson (PRT) mixing schemes (Pollard et al. 1973). Entrainment occurs in KT via instantaneous wind erosion, whereas PRT mixes through shear-induced currents at the base of the mixed layer. PWP implements buoyancy fluxes and shear-induced mixing. The main difference between KT, PRT, and PWP is that PRT and PWP use dynamic instability criteria to mix based on buoyancy and current shear whereas KT does not. Dynamic instability criteria are based on the bulk Richardson number (mixed layer stability) for PRT and both bulk and gradient Richardson numbers (shear flow stability) for PWP. In PRT and PWP, Earth’s rotation restores flow stability and eventually causes mixing to cease within an inertial period (IP).
KT and PRT schemes are used on two types of experiments: temperature-only (T only) experiments (using only temperature from WCE, CC, and PLUME) and temperature–salinity (T–S) experiments that consider both temperature and salinity profiles. PWP only considers T–S experiments since the scheme originally considers density. The former mixing schemes usually consider only temperature but are modified to also examine both temperature and salinity (section 2b and the appendix).
For all experiments, total heat flux Q is set to zero to isolate the SST response that is due solely to wind forcing. A latitude of 15°N is used to calculate the Coriolis parameter. This latitude is used since the profiles were measured in the Caribbean Sea. Simulations are run out to 0.5 IP, equivalent to 24 h at 15°N. Wind stress is estimated using the bulk aerodynamic formula
b. Definition and calculation of upper-ocean layers
The initial mixing depth h0 for T–S experiments is initialized at h0 = MLD whereas h0 = ILD for T-only experiments such that density is only a function of temperature. For T-only experiments, SST cooling will begin immediately since h0 is at the ILD. For T–S experiments, SST cooling will not occur until mixing has penetrated to the ILD from the MLD. Mixing in these experiments must erode the BL and the erosion time and cooling depend on the BLT.
3. Results and discussion
a. Inclusion of salinity and stratification strength
Less SST cooling occurs with the inclusion of salinity in all mixing schemes (Table 1). The difference in cooling between T-only and T–S experiments depends on both the thermal regime and the mixing scheme. The KT scheme leads to the largest differences in SST cooling between T-only and T–S experiments, with differences ranging from 0.5° to 6.6°C (Table 1). Differences in SST cooling between T-only and T–S experiments for PRT range from 0.1° to 1.3°C. These cooling differences arise because of the extra time for the modeled processes to erode the BL. BL erosion time is defined as the time when the MLD = ILD and SST cooling commences. KT mixing (mechanical wind stirring) takes the longest to erode the BL compared to shear-induced processes such as PRT and PWP; PWP processes result in the least time (Table 1). This indicates that instantaneous wind erosion (from KT) is the most sensitive mixing process to the inclusion of salinity such that increased stratification reduces the efficiency of wind erosion compared to shear-induced currents. Shear-induced currents (as in PRT and PWP) cut through the stratification more efficiently than the former processes. For example, Fig. 2 highlights the evolution of BL erosion for the PLUME strong S regime using PWP. SST cooling begins within 3 h of TC wind forcing, caused by increased mixed layer shear eroding the BL. Since PWP uses both mixed layer stability and shear flow instability criteria, BL erosion initiates the quickest between the schemes compared to PRT, which only uses a mixed layer stability criterion. Although KT processes take longer to erode the BL, they have a more efficient cooling rate than shear-induced processes as in PRT and PWP.
Sea surface temperature cooling (°C) and barrier layer erosion time (h) of 1D mixing experiments using profiles from Fig. 1. Simulations are run to ½IP (~24 h) to estimate ∆SST. The T-only experiments only have one value for each S experiment since T-only experiments do not use salinity. An “N” under barrier layer erosion time indicates that barrier layer did not erode.
The PLUME regime is most sensitive to the inclusion of salinity (largest cooling differences between T-only and T–S), whereas the WCE regime has the smallest differences. Similar results were found in Hernandez et al. (2016) where the authors showed that cyclone-induced SST cooling decreases with increasing BLT and increasing salinity stratification. Additionally, the PLUME thermal regime is also most sensitive to the strength of salinity stratification compared to other thermal regimes. This is attributed to the PLUME thermal regime having a shallower ILD relative to other regimes in this study, and, therefore, entrainment of thermocline waters would take less time compared to deeper ILDs. Interestingly, the shallow isothermal (~30 m) PLUME strong S regime has minimal difference in SST cooling compared to that for the deeper isothermal (~60 m) CC weak S regime (Table 1). This suggests salinity stratification in the PLUME strong S regime is reducing mixing efficiency to the point that the SST response is like that of deeper thermal structure.
b. SST response in context of a BL baroclinic wave speed
For areas that do not have a BL and MLD = ILD,
Though the experiments in this study only consider a stationary wind forcing, the influence of BL salinity stratification on coupling between a passing TC and thermocline can still be speculated using Frs. Including the influence of BL stratification would theoretically decrease the coupling between the sea surface and thermocline response during TC passage compared to that which does not consider the BL. For example, if Uh = c = 2 m s−1, values of Frs for all regimes would be greater than 1 whereas Fr would be equivalent to 1. This indicates a TC would encounter a baroclinic response in BL regions compared to a barotropic response in non-BL regions (Geisler 1970). Furthermore, values of Frs between the CC weak S and PLUME strong S regimes would be comparable for the same Uh and c criteria considering that
4. Summary and conclusions
The results of this study underscore the potential influence of salinity stratification on upper-ocean mixing and SST cooling during TC wind forcing and show the contributions of both vertical thermal and haline structure in these processes. The findings of this study are important since literature has shown that even modest differences in SST have a significant impact on air–sea heat exchange in a TC and TCs frequently pass over BL regimes in the tropical global oceans. In deep isothermal layers, the vertical salinity gradient does not contribute to SST cooling during TC wind forcing because the vertical thermal structure dominates. However, in shallow isothermal layers, SST cooling is sensitive to the vertical salinity gradient. SST cooling in a salinity-stratified environment is also sensitive to which ocean mixed layer (OML) scheme is chosen. KT processes result in the most SST cooling, although the scheme takes the longest to erode the BL; shear instability processes as in PWP erode the BL the fastest. These results are put into context through the creation of a BL baroclinic wave speed
Acknowledgments
The corresponding author acknowledges the generous funding support by NASA (Grant NNX15AG43G) and the insightful suggestions from Dr. Benjamin Jaimes (RSMAS/UM).
APPENDIX
Revisions to KT and PRT Models to Incorporate Salinity
a. KT
b. PRT
PRT and PWP use a two-dimensional wind stress of
REFERENCES
Androulidakis, Y., V. Kourafalou, G. Halliwell, M. Le Hénaff, H. Kang, M. Mehari, and R. Atlas, 2016: Hurricane interaction with the upper ocean in the Amazon-Orinoco plume region. Ocean Dyn., 66, 1559–1588, https://doi.org/10.1007/s10236-016-0997-0.
Balaguru, K., P. Chang, R. Saravanan, L. R. Leung, Z. Xu, M. Li, and J.-S. Hsieh, 2012: Ocean barrier layers’ effect on tropical cyclone intensification. Proc. Natl. Acad. Sci. USA, 109, 14 343–14 347, https://doi.org/10.1073/pnas.1201364109.
Chi, N.-H., R.-C. Lien, E. A. D’Asaro, and B. B. Ma, 2014: The surface mixed layer heat budget from mooring observations in the central Indian Ocean during Madden–Julian oscillation events. J. Geophys. Res. Oceans, 119, 4638–4652, https://doi.org/10.1002/2014JC010192.
Cione, J. J., and E. W. Uhlhorn, 2003: Sea surface temperature variability in hurricanes: Implications with respect to intensity change. Mon. Wea. Rev., 131, 1783–1796, https://doi.org/10.1175//2562.1.
de Boyer Montegut, C., J. Mignot, A. Lazar, and S. Cravatte, 2007: Control of salinity on the mixed layer depth in the world ocean: 1. General description. J. Geophys. Res., 112, C06011, https://doi.org/10.1029/2006JC003953.
Geisler, J. E., 1970: Linear theory on the response of a two layer ocean to a moving hurricane. Geophys. Fluid Dyn., 1, 249–272, https://doi.org/10.1080/03091927009365774.
Grodsky, S. A., and Coauthors, 2012: Haline hurricane wake in the Amazon/Orinoco plume: AQUARIUS/SACD and SMOS observations. Geophys. Res. Lett., 39, L20603, https://doi.org/10.1029/2012GL053335.
Hernandez, O., J. Jouanno, and F. Durand, 2016: Do the Amazon and Orinoco freshwater plumes really matter for hurricane-induced ocean surface cooling? J. Geophys. Res. Oceans, 121, 2119–2141, https://doi.org/10.1002/2015JC011021.
Jacob, S. D., L. K. Shay, A. J. Mariano, and P. G. Black, 2000: The 3D oceanic mixed layer response to Hurricane Gilbert. J. Phys. Oceanogr., 30, 1407–1429, https://doi.org/10.1175/1520-0485(2000)030<1407:TOMLRT>2.0.CO;2.
Jaimes, B., and L. K. Shay, 2009: Mixed layer cooling in mesoscale oceanic eddies during Hurricanes Katrina and Rita. Mon. Wea. Rev., 137, 4188–4207, https://doi.org/10.1175/2009MWR2849.1.
Jaimes, B., L. K. Shay, and E. W. Uhlhorn, 2015: Enthalpy and momentum fluxes during Hurricane Earl relative to underlying ocean features. Mon. Wea. Rev., 143, 111–131, https://doi.org/10.1175/MWR-D-13-00277.1.
Jaimes, B., L. K. Shay, and J. K. Brewster, 2016: Observed air-sea interactions in tropical cyclone Isaac over Loop Current mesoscale eddy features. Dyn. Atmos. Oceans, 76, 306–324, https://doi.org/10.1016/j.dynatmoce.2016.03.001.
Kraus, E. B., and J. S. Turner, 1967: A one-dimensional model of the seasonal thermocline II. The general theory and its consequences. Tellus, 19, 98–106, https://doi.org/10.3402/tellusa.v19i1.9753.
Lukas, R., and E. Lindstrom, 1991: The mixed layer of the western equatorial Pacific Ocean. J. Geophys. Res., 96, 3343–3357, https://doi.org/10.1029/90JC01951.
McPhaden, M. J., and G. R. Foltz, 2013: Intraseasonal variations in the surface layer heat balance of the central equatorial Indian Ocean: The importance of zonal advection and vertical mixing. Geophys. Res. Lett., 40, 2737–2741, https://doi.org/10.1002/grl.50536.
Mignot, J., A. Lazar, and M. Lacarra, 2012: On the formation of barrier layers and associated vertical temperature inversions: A focus on the northwestern tropical Atlantic. J. Geophys. Res., 117, C02010, https://doi.org/10.1029/2011JC007435.
Neetu, S., and Coauthors, 2012: Influence of upper‐ocean stratification on tropical cyclone‐induced surface cooling in the Bay of Bengal. J. Geophys. Res., 117, C12020, https://doi.org/10.1029/2012JC008433.
Pollard, R. T., P. B. Rhines, and R. O. R. Y. Thompson, 1973: The deepening of the wind-mixed layer. Geophys. Fluid Dyn., 3, 381–404, https://doi.org/10.1080/03091927208236105.
Powell, M., P. J. Vickery, and T. A. Reinhold, 2003: Reduced drag coefficient for high wind speeds in tropical cyclones. Nature, 422, 279–283, https://doi.org/10.1038/nature01481.
Price, J. F., R. A. Weller, and R. Pinkel, 1986: Diurnal cycling: Observations and models of the upper ocean response to diurnal heating, cooling, and wind mixing. J. Geophys. Res., 91, 8411–8427, https://doi.org/10.1029/JC091iC07p08411.
Reul, N., Y. Quilfen, B. Chapron, S. Fournier, V. Kurdyavtsev, and R. Sabia, 2014: Multisensor observations of the Amazon-Orinoco River plume interactions with hurricanes. J. Geophys. Res. Oceans, 119, 8271–8295, https://doi.org/10.1002/2014JC010107.
Rudzin, J. E., L. K. Shay, B. Jaimes, and J. K. Brewster, 2017: Upper ocean observations in eastern Caribbean Sea reveal barrier layer within a warm core eddy. J. Geophys. Res. Oceans, 122, 1057–1071, https://doi.org/10.1002/2016JC012339.
Shay, L. K., and E. Uhlhorn, 2008: Loop Current response to Hurricanes Isidore and Lili. Mon. Wea. Rev., 136, 3248–3274, https://doi.org/10.1175/2007MWR2169.1
Shay, L. K., A. J. Mariano, S. D. Jacob, and E. H. Ryan, 1998: Mean and near-inertial ocean current response to Hurricane Gilbert. J. Phys. Oceanogr., 28, 858–889, https://doi.org/10.1175/1520-0485(1998)028<0858:MANIOC>2.0.CO;2.
Sprintall, J., and M. Tomczak, 1992: Evidence of the barrier layer in the surface layer of the tropics. J. Geophys. Res., 97, 7305–7316, https://doi.org/10.1029/92JC00407.
Vissa, N. K., A. N. V. Satyanarayana, and B. P. Kumar, 2013: Response of upper ocean and impact of barrier layer on Sidr cyclone induced sea surface cooling. Ocean Sci. J., 48, 279–288, https://doi.org/10.1007/s12601-013-0026-x.
Wang, X., G. Han, W. Qi, and W. Li, 2011: Impact of barrier layer on typhoon-induced sea surface cooling. Dyn. Atmos. Oceans, 52, 367–385, https://doi.org/10.1016/j.dynatmoce.2011.05.002.
Yan, Y., L. Li, and C. Wang, 2017: The effects of oceanic barrier layer on the upper ocean response to tropical cyclones. J. Geophys. Res. Oceans, 122, 4829–4844, https://doi.org/10.1002/2017JC012694.