
Citation: | XIA Changshui, JUNG KyungTae, WANG Guansuo, YIN Xunqiang, GUO Jingsong. Case study on the three-dimensional structure of meso-scale eddy in the South China Sea based on a high-resolution model[J]. Acta Oceanologica Sinica, 2016, 35(2): 29-38. doi: 10.1007/s13131-016-0805-1 |
The Southern Ocean (SO, 30°–85°S, 0°E–180°–0°W) is an important field for the global Meridional Overturning Circulation (MOC) and deep ocean water formation (Kuhlbrodt et al., 2007). The SO is also a key region for the global carbon cycle, air-sea exchange and heat redistribution (Arrigo et al., 2008; Park et al., 2010; Sarmiento et al., 1998). The mixed layer (ML) refers to a near ocean surface water layer with an almost vertical uniformity of temperature, salinity, density and other physical properties, resulting from the wind mixing and effects of ocean internal processes (Carton et al., 2008; Holte and Talley, 2009; Liu et al., 2009). As a bridge between the atmosphere and the ocean, the properties and dynamical processes in the ML have significant effects on the exchange of mass, heat and momentum between air and sea (Kara et al., 2003; de Boyer Montégut et al., 2007). Thus, the mixed layer depth (MLD) in the SO plays a crucial role in global climate regulation. On the other hand, there is a large amount of phytoplankton living in the ML. The variability of MLD has great influences on the availability of light and nutrients, modulating the heat and nutrient exchange between the well-lit surface layer and the dark deep layer, which dominates the growth and distribution of phytoplankton (Mitchell et al., 1991; Sallée et al., 2010). Besides, changes in phytoplankton can further affect primary productivity and even ecosystems through food webs, especially in the SO where the krill is abundant. Dynamic processes of ML in the SO can thus alter the ocean capacity of storing heat and carbon, as well as the availability of light and nutrients for phytoplankton growth. Above all, the SO ML has a pivotal impact on global climate change and ecosystems due to its unique characteristics. Therefore, it is important to study the variability of MLD in the SO, particularly the links with the dominant mode of atmospheric and oceanic variabilities.
Numerous studies in the literature have linked the variability of the SO MLD to sea ice, Amundsen Sea low (ASL), El Niño-Southern Oscillation (ENSO) and the Southern Annular Mode (SAM) through observations and models on various time scales (Gupta and England, 2006; Li et al., 2019; Nardelli et al., 2017; Sallée et al., 2010; Vivier et al., 2010). However, the SO MLD responds most significantly to the SAM and the ENSO, and the influences of the SAM and the ENSO are dominant on the interannual time scale (Li and England, 2020; Sallée et al., 2010). Meanwhile, the response of the SO MLD to the SAM shows a wave-3 pattern in the latitudinal band between 35°S and 55°S, with the maximum anomaly close to New Zealand, south-west of Australia and the east of Chile. The wave-3 pattern is caused by the combined effects of net heat flux and meridional wind (Sallée et al., 2010; Zhang et al., 2018). In addition, the sea surface temperature anomaly (SSTA) in the South Pacific also exhibits a similar tripole-like structure, suggesting that the MLD variability has an impact on SST (Li et al., 2014b).
As for ENSO, although it occurs in the tropics, previous studies have shown that its effects extend to extratropical variations of the southern hemisphere and even high-latitude regions covered by sea ice (Ciasto and England, 2011; Schneider et al., 2012; Turner, 2004). Meanwhile, previous studies have revealed the mechanisms of ENSO-induced anomalies of specific climate variables in the southern hemisphere, such as surface air temperature, sea surface temperature (SST), sea ice concentration, and MLD (Ciasto and England, 2011; Ding et al., 2011; Li and England, 2020). The impact of ENSO on the SO MLD is mainly concentrated in the Pacific sector. Using the model data, Li and England (2020) found that the anomalous ASL intensification associated with the ENSO affects the variation of the MLD in the Pacific sector during the austral winter (June, July, and August; JJA). In turn, the MLD response is characterized by a dipole pattern with a positive MLD anomaly center in the central Pacific (surrounding 120°W) and a negative MLD anomaly center in the southwest Pacific (surrounding New Zealand). This is closely related to the airflow around the anomalous sea level pressure (SLP) center in the Amundsen Sea.
The ENSO can modulate the austral winter MLD in the South Pacific via the teleconnection between it and the ASL (Li and England, 2020). The teleconnection is realized by the Rossby wave train propagation emanating from the equatorial Pacific Ocean into the Antarctic (Hoskins and Karoly, 1981; Turner, 2004), which is known as the Pacific-South American (PSA) teleconnection pattern (Kidson, 1999; Mo and Ghil, 1987; Mo and Higgins, 1998). Positive SSTA in the tropical Pacific related to the positive phase ENSO can be regarded as an effective Rossby wave source, which excites the stationary Rossby wave train in the upper troposphere (Cai et al., 2011; Clem and Fogt, 2015; Fogt et al., 2011; L’Heureux and Thompson, 2006), which propagates poleward, with the alternating occurrence of high pressure anomaly centers and lowpressure anomaly centers, eventually reaching ASL and interacting with the ASL (Li et al., 2021; Meehl et al., 2016; Zhang et al., 2021). The ASL, in turn, is abnormally deepened. As a consequence, colder southerly wind anomaly and warmer northerly wind anomaly appear to the west and east of the anomalous cyclone respectively, thus contributing to the dipole distribution of MLD in the South Pacific (Li and England, 2020).
Moreover, growing number of studies suggest that tropicalpolar teleconnections can also originate from other tropical oceans except for the Pacific Ocean (Cai et al., 2019; Li et al., 2014a, 2015a; Nuncio and Yuan, 2015). They affect the SO on different time scales, either through individual or combined impact (Li and England, 2020; Meehl et al., 2016; Rondanelli et al., 2019; Zhang et al., 2021). In addition to the ENSO, signals from the tropical Indian Ocean, especially the Indian Ocean Dipole (IOD), also play an instrumental role in the generation and propagation of the Rossby wave train (Cai et al., 2014; Luo et al., 2010; Saji et al., 2005). When these two major tropical climate variabilities coexist, such as El Niño accompanied by a positive phase IOD, the Rossby wave train can reach farther to the Bellingshausen Sea. Therefore, they can work together to influence the SO MLD (Cai et al., 2011; Li and England, 2020; Zhang et al., 2021).
The aforementioned studies reveal that the tropical diabatic heating caused by the El Niño can act as the effective Rossby wave source, and trigger a southeastward propagating Rossby wave train in the upper troposphere, thus affecting the SO MLD (Hoskins and Karoly, 1981; Karoly and Hoskins, 1983; Li and England, 2020; Turner, 2004). The different categorization of El Niño was used by previous investigators. For example, Yu and Kao (2007) categorized the El Niño events into two types based on the location of maximum SSTA, namely, the central Pacific (CP) El Niño with maximum SSTA occurring in the central tropical Pacific and the eastern Pacific (EP) El Niño with maximum SSTA in the eastern tropical Pacific (Kao and Yu, 2009; Yu and Kao, 2007). Similarly, Kug et al. (2009) categorized the El Niño events into cold tongue (CT) El Niño and warm pool (WP) El Niño based on the spatial distribution of SSTA. The CT El Niño is typically known as the conventional El Niño, which is distinguished by the presence of considerable SSTA in the Niño 3 region (between 5°S and 5°N and 150°W to 90°W). In contrast, the WP El Niño is characterized by the occurrence of SSTA that are mostly restricted to the Niño 4 region (between 5°S and 5°N and 160°E to 150°W).
This study follows the definitions by Yu and Kao (2007). Using reanalysis data, Wang et al. (2020) examined the effects of EP and CP El Niño events on the SO MLD from 1979 to 1998. They found that the response of the MLD to CP and EP El Niño is relatively strong in austral winter and spring. The MLD response exhibits a seesaw pattern, located on the flanks of South America. The positive MLD anomaly center of the seesaw arises in the Southwest Atlantic and the negative anomaly center occurs in the Southeast Pacific. This feature could be found in two types of El Niño events in spring, and it is mainly controlled by sea surface heat flux anomaly. However, the seesaw pattern disappears during the austral winter for both EP and CP El Niño events. Two types of El Niño events are quite distinct in terms of generation and development mechanisms (Chen et al., 2019), so that have different impacts on the SO MLD. However, previous studies focused mainly on the impacts of classified El Niño events on the SO MLD (Wang et al., 2020). It has not been well understood what are the distinct effects and mechanisms of EP and CP El Niño events on the variation of SO MLD. Compared with the unclassified El Niño event, the impacts of EP and CP El Niño events on the SO MLD are more complex. Thus, this study aims to explore the different responses of the SO MLD to two-type El Niño events during austral winter when the influence of EP and CP El Niño events on the SO MLD is greatest (Li and England, 2020).
The remainder of this paper is organized as follows. Section 2 briefly describes the datasets and methods used in this paper. Next, MLD variabilities induced by EP and CP El Niño events are presented respectively in Section 3. Section 4 investigates the different characteristics of tropical-polar teleconnections during the EP and CP El Niño events. Section 5 further illustrates possible modulation mechanisms of MLD variability during the EP and CP El Niño events. Discussion and conclusions are given in Section 6.
This study uses the monthly data derived from the European Centre for Medium-Range Weather Forecasts (ECMWF) Ocean Reanalysis System version 5 (ORAS5) from January 1980 to December 2021 as a baseline (Zuo et al., 2019). The horizontal resolution is 0.25° × 0.25° with 75 layers in the vertical direction. Variables include potential temperature, salinity, and sea ice concentration (SIC), where potential temperature and salinity are used to estimate the MLD in the SO. Another primary dataset for this study is the fifth-generation ECMWF reanalysis data (ERA5; Bell et al., 2021; Hersbach et al., 2020). The ERA5 has a horizontal resolution of 0.25° × 0.25° and covers the period from January 1980 to December 2021. Geopotential height from 300 hPa, 700 hPa, and 850 hPa, SLP, horizontal wind fields at 300 hPa and 10 m, air temperature at 2 m, surface specific humidity (
In this study, the MLD is defined as the minimum depth of temperature-criteria mixed layer (MLt), salinity-criteria mixed layer (MLs) and density-criteria mixed layer (MLd) (Helber et al., 2012). To unify the criteria for calculating the MLD for each variable so that the calculation results can be compared, 0.2℃ (
$$ \Delta \theta = \theta \left( {{T_{10}} - \Delta T,{S_{10}},{P_0}} \right) - \theta \left( {{T_{10}},{S_{10}},{P_0}} \right), $$ | (1) |
$$ \Delta S = 0.01n,\;n = 1, 2,..., N_{\mathrm{mld}},$$ | (2) |
$$ \Delta {\theta _{\text{s}}} = \theta \left( {{T_{10}},{S_{10}} + \Delta S,{P_0}} \right) - \theta \left( {{T_{10}},{S_{10}},{P_0}} \right), $$ | (3) |
$$ \Delta {\theta _{\text{s}}} = \Delta \theta , $$ | (4) |
where
We use Nep and Ncp to categorize the EP and CP El Niño events according to Ren and Jin (2011). These are presented as
$$ \left\{ {\begin{array}{*{20}{c}} {{N_{{\text{ep}}}} = {N_3} - \alpha {N_4}}, \\ {{N_{{\text{cp}}}} = {N_4} - \alpha {N_3}} , \end{array}} \right. $$ | (5) |
$$\alpha = \left\{ {\begin{array}{*{20}{l}} 2/5,& {N_3}{N_4} > 0 ,\\ 0,& {\text{otherwise}} , \end{array}} \right.$$ | (6) |
and N3 and N4 are the Niño 3 and Niño 4 indices, respectively. The indices are averaged monthly SSTA in the Niño 3 region (5°N–5°S, 150°–90°W) and Niño 4 region (5°N–5°S, 160°E–180°–150°W), respectively.
Before considering the MLD changes caused by two types of El Niño events, it is first important to consider which El Niño years belong to the EP and CP El Niño events. According to previous studies, Niño 3 and Niño 4 indices can be adequately used to represent EP and CP El Niño events (Kug et al., 2009; Yeh et al., 2009). However, Ren and Jin (2011) further proposed two new El Niño indices based on Niño 3 and Niño 4 indices (see Eq. (5)). Figure 1 shows the time evolution of Niño 3 and Niño 4 indices (blue lines) and two new El Niño indices named Nep and Ncp (red lines). Consider that El Niño events typically reach peaks during boreal winter (December, January, and February; DJF) (Zhang et al., 2019), only boreal winter years with indices larger than one standard deviation are selected for this study to avoid interfering with subsequent composite results. One can see that Niño 3 and Nep are highly correlated (Fig. 1a, r = 0.95). The difference between the two is generally minor, especially in years with the strongest EP El Niño, such as 1982/1983 and 1997/1998. It indicates that both Niño 3 and Nep indices can capture EP El Niño events well. Although Niño 4 index is also highly correlated with the Ncp index (Fig. 1b, r = 0.82), the difference between Niño 4 and Ncp indices is more pronounced compared to the difference between Niño 3 and Nep (Fig. 1a). Specifically, during the periods of the strongest El Niño events, Niño 4 and Ncp indices are almost out of phase, such as 1982/1983, 1997/1998, and 2015/2016.
In general, fluctuation amplitudes of Nep and Ncp indices tend to be a little smaller than that of Niño 3 and Niño 4 indices. That is, some years may be classified as the EP El Niño or the CP El Niño events if classified according to Niño 3 and Niño 4 indices, but may not be when classified by Ncp and Nep. In addition, when the CP El Niño events occur, Nep becomes less than one standard deviation in most cases, while the Niño 3 indices are still larger than one standard deviation, such as 1994/1995, 2014/2015, and 2018/2019. This indicates that Ncp and Nep indices are relatively independent, while Niño 3 and Niño 4 indices are interrelated. According to these characteristics, Ncp and Nep indices thus are used as the basis for dividing two types of El Niño events. Ultimately, as listed in Table 1, five EP El Niño and eight CP El Niño events are selected during 1980–2021. It is worth noting that both Ncp and Nep indices exceeded one standard deviation during 1987/1988, so it was classified as mixed type, starting with the EP El Niño and subsequently shifting into the CP El Niño (Wang et al., 2018a). Comparing with the previous finding using the same Nep and Ncp indices, the selection of EP El Niño years is identical. However, the selection of CP El Niño years is different. This study has two more CP Niño years, namely 1990/1991 and 2019/2020. This difference may be attributed to different SST data used to calculate N3 and N4 in Eq. (2), as well as different study periods. The analysis period for the previous study was from 1980 to 2019, while that for this study from 1980 to 2021 (Zhang et al., 2021). In addition, the comparison is made to the conventional definition of two-types of El Niño. The years of EP and CP El Niño are robust over the same research period, except for the mixed type of El Niño event (such as 1987/1988) (Wang et al., 2017; Zhang et al., 2011).
Category | Year |
EP El Niño | 1982/1983, 1991/1992, 1997/1998, 2006/2007, 2015/2016 |
CP El Niño | 1990/1991, 1994/1995, 2002/2003, 2004/2005, 2009/2010, 2014/2015, 2018/2019, 2019/2020 |
Figure 2 shows composite images of seasonal mean MLD anomaly in the SO for five EP and eight CP El Niño events listed in Table 1. Overall, the MLD anomaly in the Pacific sector is evident in both EP and CP El Niño events throughout the year, while the MLD anomaly in the Atlantic and Indian Ocean sectors is weaker. The notable anomalous regions are mainly concentrated in the high latitude region between 40°S and 70°S. In austral spring (September, October, and November; SON) during EP El Niño events, the distribution of MLD anomaly shows a dipole pattern in the Pacific sector, with the positive and negative anomaly centers in the central and eastern Pacific, respectively (Fig. 2a). The other sectors have no obvious distribution features. Although the spatial structure is obvious, the MLD anomaly is not statistically significant at the 95% confidence level. For the austral spring of CP-El Niño, the distributions of MLD anomaly in the SO bear a resemblance to that of EP El Niño events, except that the phases are opposite. Furthermore, only weak and statistically significant negative MLD anomaly occurs in the central Pacific Ocean (Fig. 2e). During the subsequent austral summer (DJF) and autumn (March, April, May; MAM), the distribution characteristics of MLD anomaly in the SO are similar to that during the EP El Niño, with weak and statistically significant negative anomaly centers in the central Pacific. The remaining regions show positive anomaly and are almost statistically insignificant (Figs 2b and c). Similarly, during the CP El Niño, the spatial structure is similar in austral summer and autumn, with the positive anomaly occurring mainly in the central Pacific and southern Australia. The negative anomaly is found in the remaining oceans. Except the Weddell Sea, other areas are statistically insignificant in austral summer and autumn (Figs 2f and g).
In the following austral winter during the years of EP El Niño, the feature of MLD anomaly in the SO exhibits a dipole structure, with positive anomaly (Box 1) over the central and western Pacific (40°–60°S, 160°E–145°W) and negative anomaly (Box 2) over the eastern Pacific (45°–65°S, 140°–70°W). The Indian and Atlantic sectors show a weak and statistically insignificant positive anomaly (Fig. 2d). One can see two negative MLD anomaly centers and one positive anomaly center in austral winter in the CP El Niño cases. Negative anomaly centers occur in the central Pacific and Indian Ocean, while the positive anomaly center occurs in southern Australia. But only the negative anomaly center in the central Pacific (42°–62°S, 165°–100°W) features the strong and statistically significant negative anomaly (Box 3), which is distinct from that of EP El Niño. The above results indicate that both EP and CP El Niño events have dominant impacts on the MLD anomaly in the southern Pacific in austral winter compared with the other three seasons, while the influences are markedly different. We find that the MLD anomaly associated with EP El Niño shows a dipole-like pattern in the Pacific sector during the austral winter (Boxes 1 and 2 in Fig. 2d), while during the CP El Niño, only a separate negative anomaly center (statistically significant) exists in the Pacific sector (Box 3 in Fig. 2h). These response characteristics are different from previous investigations, which indicated that the MLD dipole emerges in spring during EP El Niño. Moreover, they found the MLD dipole located in the Pacific (negative anomaly center) and Atlantic (positive anomaly center) sectors (Wang et al., 2020). Different data sets and study periods, and different methods for defining the EP and CP El Niño and the MLD may all lead to differences in the results. Although the effects of the two types of El Niño on the MLD are significantly different in austral autumn, the impact of CP El Niño on MLD is not statistically significant. In contrast, in austral winter, the effects of two types of El Niño on MLD are statistically significant, which makes comparisons of the different effects valid. On the other hand, austral winter is the deepest season for the MLD, which has a significant impact on both SO primary productivity and global climate change. In addition, although the El Niño during the austral winter is in the decaying phase, the effect on the MLD is strongest at this season. This prompts this study to primarily focus on the different impacts of two types of El Niño events on the MLD of the Pacific sector in austral winter.
Previousnresults have revealed that the two types of El Niño events are accompanied by SSTA at different spots. This section first examines the SSTA triggered by the EP and CP El Niño events in the Pacific sector (Fig. 3), where the austral winter MLD anomaly is most pronounced (Figs 2d and h). During the EP El Niño, the SSTA is most obvious near the equator (Fig. 3a). The positive SSTA occurs in the eastern Pacific south of the equator, extending westward to 180°. The positive SSTA center appears off western Peru, yet the SSTA here is statistically insignificant. The negative SSTA is zonally distributed along the equator with longitudes ranging from 180° to 100°W. It is also not statistically significant. Instead, the statistically significant region appears in the extratropical Pacific between 160°E and 70°W, characterized by a dipole distribution pattern as well. It is noted that the dipole position of SSTA is consistent with that of the MLD dipole (Figs 2d and 3a). The negative SSTA occurs exactly where the ML abnormally deepens (Box 1), while positive SSTA occurs where the ML shoals (Box 2). In addition to the tropical Pacific, a small proportion of significant positive SSTA exists in the northern Indian Ocean near the coast. In comparison, the maximum positive SSTA of CP El Niño emerges in the central Pacific, with 10°–30° of westward shifting compared with the EP El Niño events (Fig. 3b). Although the magnitudes of this positive anomaly are smaller than that of the EP El Niño, they are statistically significant at the 95% confidence level. During the CP El Niño events, the positive SSTA occurs to the northwest of Box 3 and the position of positive SSTA is not completely coincident with Box 3 (Fig. 3b). Additionally, statistically significant positive SSTA is found in almost entire the northern Indian Ocean, which is different from the EP El Niño. The statistically significant SSTA and MLD anomaly are in good agreement, suggesting a possible connection between the two. The mechanisms will be discussed in a later section.
Previous studies have shown that the positive SSTA associated with the El Niño in the tropical Pacific can heat the atmosphere and enhance atmospheric convection (Cai et al., 2011; Clem and Fogt, 2015; Karoly and Hoskins, 1983). It can be regarded as an effective Rossby wave source, which can excite the Rossby wave train propagating poleward in the upper troposphere. Figure 4 shows the austral winter geopotential heights of the upper, middle, and lower troposphere during the EP and CP El Niño events. One can see that in the upper troposphere, the Rossby wave train propagates from the tropical region to the Antarctic region via alternating cyclone anomaly and anticyclone anomaly on the 300 hPa geopotential height anomaly, which is known as the Pacific-South American (PSA) pattern (Figs 4a and b). During the EP El Niño, the 300 hPa geopotential height anomaly extends along an approximate great circle path, originating from the eastern Australia and propagating southeastward into the Southeast Pacific and finally into the Amundsen Sea, characterized by a positive (far east of Australia)-negative (Ross Sea–Amundsen Sea) anomaly centers (Fig. 4a). As a result, low-pressure anomaly center, known as the cyclonic anomaly, dominates the Ross Sea–Amundsen Sea. During the CP El Niño, the Rossby wave train has a similar spatial propagation path to that of EP-Niño, featuring positive (eastern Australia)-negative (southern New Zealand)-positive (Amundsen Sea–Bellingshausen Sea) anomaly centers. In contrast to the EP El Niño, the Amundsen–Bellingshausen seas are controlled by the positive geopotential height anomaly center, namely anticyclonic anomaly (Fig. 4b). Since the interval of geopotential height contours is 5 m, the PSA structure of the EP El Niño is stronger than that in the CP El Niño cases. In the middle and lower troposphere, the phenomenon of Rossby wave train propagation is weaker or even disappears, but the Ross Sea–Amundsen Sea and the Amundsen Sea–Bellingshausen Sea are still dominated by a low-pressure anomaly (EP El Niño) and high-pressure anomaly (CP El Niño) centers, respectively (Figs 4c–f). The above analysis results indicate that tropical Pacific warm SSTA with different magnitudes and placements induce different geopotential height anomalies near the Amundsen Sea in terms of intensity and location.
The Rossby wave train only spreads poleward during three seasons rather than summer, when the subtropical jet is sufficiently strong and properly positioned to act as a waveguide (Li et al., 2015b; Schneider et al., 2012). Figure 5 shows the austral winter 300 hPa zonal wind speed anomaly together with the PV anomaly during the two types of El Niño. During the EP El Niño, the subtropical jet (38°–60°S, 160°E–180°–70°W) at the MLD dipole region is substantially stronger than the climatology, ranging from southeast of Australia to west of Chile around 50°S (Fig. 5a). From Fig. 5c, in the meridional direction, the PV anomaly exhibits distinctive positive-negative-positive distribution patterns in the South Pacific. The positive PV anomaly (25°–50°S, 160°E–180°–70°W) occurs between 25°S and the northern dipole pattern, ranging latitudinally from southeastern Australia to western Chile. Negative PV anomaly (40°–65°S, 160°E–180°–70°W) occurs at the dipole pattern, with the negative anomaly center located on the western border of Box 2. Another positive PV anomaly (65°–80°S, 160°E–110°W) exists to the south of the negative PV anomaly. Moreover, positive (25°–50°S, 160°E–180°–70°W) and negative PV anomalies (40°–65°S, 160°E–70°W) appear to the north and south of the subtropical jet (38°–60°S, 160°E–180°–70°W), respectively (Fig. 5c). This exceptionally strong subtropical jet and the meridional PV gradient in its vicinity facilitate the propagation of the Rossby wave train. However, in the case of the CP El Niño, the subtropical jet is not visible and the meridional PV gradient is substituted by the zonal gradient (Figs 5b and d). This illustrates that the intensity of the PSA pattern may also influence the circulation near Amundsen Sea in addition to the different heat sources in the tropical Pacific. Conditions favoring southeastward propagation of Rossby wave trains may result in a stronger PSA pattern.
The last section investigates poleward propagation of the El Niño-induced Rossby wave trains in the upper troposphere, thus this section further explores responses of near-surface ocean circulation. Figure 6 shows the SLP anomaly and 10 m wind speed anomaly attributed to the EP and CP El Niño. The positive SSTA associated with the EP El Niño excites a Rossby wave train propagating towards the vicinity of the Amundsen Sea, which intensifies the ASL (Fig. 6a). Consequently, a strong anomalous cyclone formed in the Ross Sea–Amundsen Sea–Bellingshausen Sea. Strong southerly and northerly wind speed anomalies occur to the west and east of the anomalous cyclone, respectively, where the MLD dipole occurs (Figs 2d and 6c). In contrast, the CP El Niño events cause anomalous anticyclone only in the Amundsen Sea–Bellingshausen Sea (Fig. 6b). Remarkable northerly wind speed anomaly occurs to the west of this anomalous anticyclone, precisely where the negative MLD anomaly occurs (Figs 2h and 6d). Despite the cold southerly winds to the east of the anomalous anticyclone, it fails to directly affect the MLD in Box 3, as it is located to the east of Box 3. One can see that the strong anomalous cyclone appears in the MLD dipole region (Boxes 1 and 2) during the EP El Niño, while the strong anomalous anticyclone occurs in the region of negative MLD anomaly (Box 3) during the CP El Niño events, which enhances the lateral transport of the corresponding meridional winds of two flanks. It suggests that the meridional winds associated with anomalous cyclones and anticyclones are critical factors to create distinct effects of two types of El Niño events on the MLD anomaly. In addition, the anomalous cyclone is only observed during the EP El Niño events, which is similar to the response of the SLP anomaly field to unclassified ENSO (Li and England, 2020). However, the Amundsen Sea is controlled by the anomalous anticyclone in the CP El Niño, which has not been studied yet.
The different sources of meridional winds exhibit different properties, particularly temperature and humidity. Figure 7 shows temperature anomaly at 2 m and
Furthermore, comparisons of Figs 3a and b to Figs 7a and b show that the distribution of temperature anomaly at 2 m is quite similar to that of the SSTA during the EP and CP El Niño events. However, the fluctuation amplitudes of the SSTA are smaller than that of the 2 m temperature anomaly, implying that the heat exchange exists between the atmosphere and the ocean through sensible heating and then affect the MLD. Therefore, the total heat flux at the air-sea interface is further investigated. Figure 8 shows components of total heat flux anomaly during two types of El Niño. The positive anomaly indicates that the ocean gains the heat from the atmosphere. One can see that during the two types of El Niño, SH and LH anomalies are more pronounced than LW and SW anomalies in the focus regions. During the EP El Niño events, the negative SH anomaly (southeast of New Zealand) and the positive SH anomaly (southwest of Chile) show up in the MLD dipole region (Fig. 8a). In contrast, the negative LH anomaly extends eastward to 100°W compared to the negative SH anomaly (Fig. 8c), which may counteract the warming effect of the positive LH anomaly in Box 2. The LW anomaly is similar to the LH anomaly, but the magnitudes are smaller (Fig. 8e). Comparing to the other three components, the SW anomaly at the MLD dipole region are negligible (Fig. 8g). Accordingly, the spatial distributions of total surface heat flux anomaly are similar to that of the LH anomaly (Figs 8c and i). During the CP El Niño events, except for the SW anomaly, the other three components gain the heat in Box 3 (Figs 8b, d, and f). The SW anomaly here shows a weak and statistically significant negative anomaly with the heat loss (Fig. 8h). The reason for the negative SW anomaly needs to be further explored. Thus, the total surface heat flux anomaly in the CP El Niño cases acquires the heat at the negative MLD anomaly region (Box 3 in Fig. 8j).
These signatures may preliminarily explain the different results of SSTA and MLD anomaly during two types of El Niño events, i.e., the MLD (SSTA) dipole in the EP El Niño cases and separate negative MLD (positive SSTA) anomaly center in the CP El Niño cases. During the EP El Niño events, the ocean loses the heat off the southeastern of New Zealand (Box 1, positive MLD anomaly region), where convection is enhanced and the ML deepens. In contrast, in the southwest of Chile (Box 2, negative MLD anomaly region), the ocean gains the heat from the atmosphere and the ML becomes shallower, even though the heat loss here counteracts the warming effect. As for the CP El Niño events, the ocean acquires the heat from the atmosphere in the central subtropical Pacific (Box 3, negative MLD anomaly region), which enhances stratification and therefore, ML shoals.
Further quantitative analysis of the total surface heat flux anomaly and its components show that during EP El Niño, both the negative surface SH, LH, and SW anomaly contribute to the SST cooling and ML deepening (Box 1), especially the surface SH and LH anomaly (Fig. 9). The regional averaged total heat flux anomaly reaches –4.5 W/m2, accompanied by SH, LH and SW reaching –1.9 W/m2, –1.9 W/m2 and –0.7 W/m2, respectively. In comparison, the contribution of the total heat flux is negligible in the negative MLD anomaly region (Box 2), characterized by values close to 0, due to the counterbalance of heat gain and loss in this region. During the CP El Niño events (Box 3), the LH reaching 3.6 W/m2 is the dominant contributor to SST warming and ML shoaling. In addition, the increased SH and LW of 2.7 W/m2 and 2.6 W/m2 are two other important contributors to SST warming and ML shoaling. The SW cools the SST with –0.8 W/m2. Thus, the total heat flux anomaly is 8 W/m2.
In combination with the surface heat flux anomaly described above, we next examine the austral winter MLD variability driven by surface heat flux anomaly. As shown in Figs 2d and h, the MLD dipole associated with the EP El Niño events occurs in the South Pacific, while only the single negative MLD anomaly center can be found during the CP El Niño events. To comprehensively understand the response of MLD anomaly to El Niño-induced surface heat flux anomaly, the regional means of the three MLD anomaly boxes is shown in Fig. 10. One can see that the MLD increases by 26 m and decreases by 40 m in the MLD dipole region in the EP El Niño cases, while it decreases by 15 m in the CP El Niño cases. It is worth noting that the maximum MLD variation occurs synchronously with the almost zero total heat flux anomaly in EP cases (Box 2 in Figs 9 and 10), suggesting that the surface heat flux is not the primary factor leading to the negative MLD anomaly in Box 2 (southwest of Chile).
The question then arises as to what is responsible for the negative MLD anomaly during the EP El Niño events. The earlier analysis proved that the EP El Niño-related ASL intensifies via the PSA pattern. Subsequently, the enhanced poleward northerly winds to the east of the anomalous cyclone bring warmer air from lower latitudes into the vicinity of the Amundsen Sea. Even though the surface heat flux has no contribution to this negative MLD anomaly in Box 2, the warmer air promotes sea ice melting and facilitates the maintenance of fresh water, which may result in a negative MLD anomaly. To confirm this hypothesis, the SIC anomaly in the EP and CP El Niño cases is shown in Fig. 11. One can see that a strong and statistically significant negative SIC anomaly prevails to the south of Box 2 during the EP El Niño events. Negative SIC anomaly suggests sea ice melting or reduction in sea ice formation, which results in available freshwater on the sea surface. Negative SIC anomaly occurs to the south of Box 2 and does not directly affect MLD change within Box 2. However, this area is dominated by prevailing westerly winds, implying that the northward Ekman transport brings the freshwater from the south into Box 2. As a result, stratification is enhanced and the MLD becomes shallower. In comparison, strong but statistically insignificant positive anomaly and negative SIC anomaly present in the southern Box 1 during the EP El Niño and Box 3 during the CP El Niño, respectively. These suggest that sea ice is a predominant factor in the shallowing of MLD in Box 2, while total surface heat flux anomaly has no effect. Sea ice affects the MLD by varying the strength of stratification. This is inconsistent with previous conclusions, which indicated that the MLD anomaly associated with the EP El Niño is governed by surface heat flux anomaly and wind anomaly (Li and England, 2020; Wang et al., 2020).
In summary, we have shown distinct responses of the austral winter MLD in the SO to the two types of El Niño, i.e., the SO MLD dipole (Box 1 and Box 2) during the EP El Niño events and the isolated negative MLD anomaly center (Box 3) during the CP El Niño events. The corresponding mechanism is summarized in Fig. 12. Positive SSTA in the central and eastern Pacific associated with the EP El Niño events warm the atmosphere and enhance atmospheric convection, which can act as an effective source of the Rossby wave train. The Rossby wave train propagates southeastward along a great circle path in the upper troposphere, ultimately reaching the vicinity of the Amundsen Sea. As a result, the ASL intensifies, enabling the southerly winds to the west of the anomalous cyclone to deliver colder and drier air to the southeast of New Zealand (Box 1). The strong negative air temperature anomaly and weak negative SSTA cause temperature differences at the sea-air interface, resulting in the heat loss on the sea surface, and leading to MLD deepening in Box 1. However, on the other side of the anomalous cyclone, northerly winds transport warmer and wetter air to the southwest of Chile (Box 2). Although the effect of sea surface heat flux is negligible, the warmer air promotes sea ice melting and maintenance of fresh water with subsequent stratification enhancement, thus making MLD shallower in Box 2 (Fig. 12a). During the CP El Niño events, the heat source appears in the central Pacific. The propagation path of the Rossby wave train is similar to that of the EP El Niño, except that the Amundsen Sea is dominated by an anomalous anticyclone. Thus, the poleward winds from the west of the anomalous anticyclone transport warmer and wetter air to the southern Pacific (Box 3). This contributes to warmer air than the ocean. The ocean gains the heat from the atmosphere, thus forming a negative MLD anomaly in Box 3 (Fig. 12b).
From Fig. 10, one can see that the ML variations caused by the EP El Niño events are larger than that of the CP El Niño events. This may be caused by the following three reasons. First, the warm SSTA signal induced by the EP El Niño is stronger than that by the CP El Niño (Fig. 3), which is consistent with the previous results (Kao and Yu, 2009; Ren and Jin, 2011). Therefore, the EP El Niño creates stronger anomalous disturbances in the atmosphere. Second, years of the EP El Niño events coincide with the positive phase of the IOD, and the tropical Indian Ocean also has a Rossby wave train propagating southeastward into the Amundsen Sea (Fig. 4a), resulting in a stronger PSA pattern due to the combined effect of the EP El Niño and the IOD (Zhang et al., 2021). Finally, in the EP El Niño cases, the abnormally strong subtropical jet and the meridional gradient of the PV anomaly are conducive to propagation of the Rossby wave train (Li et al., 2021), which are absent in the CP El Niño cases (Fig. 5).
Based on ORAS5 and ERA5 reanalysis data, we investigate the quite different responses of the MLD to the two types of El Niño. The strongest MLD response to both EP and CP El Niño events occurs in austral winter, accompanied by an EP-induced SO MLD dipole pattern and a separate CP-induced negative MLD anomaly center. The two types of El Niño have different strengths of tropical warm SST anomalous signals. Therefore, the PSA patterns vary in strength, i.e., the southeastward propagation of the Rossby wave train. As a result, the ASL strengthens during the EP El Niño, with an anomalous cyclone dominating the Amundsen Sea, while the ASL weakens during the CP El Niño, with an anomalous anticyclone controlling the Amundsen Sea. During the EP El Niño events, the positive MLD anomaly occurs in the southeast of New Zealand (Box 1, 40°–60°S, 160°E–180°–145°W), while the negative anomaly occurs in the southwest of Chile (Box 2, 45°–65°S, 140°–70°W). The southerly winds to the west of the anomalous cyclone exacerbate the temperature and humidity differences between the air and the sea surface, subsequently causing variations in the surface heat flux. Among the four components of the total surface heat flux, negative SH anomaly and LH anomaly are the main contributors to surface cooling. This makes the MLD deepen. To the east of the anomalous cyclone, the role of sea surface heat fluxes is negligible. However, the warm air from the tropics melts the sea ice and maintains fresh water, thus resulting in the shallower MLD. The separate negative MLD anomaly center in the CP El Niño cases occurs in the far southeast of New Zealand (Box 3, 42°–62°S, 165°–100°W). Warmer and wetter air to the west of the anomalous anticyclone is transported to Box 3, increasing the temperature difference between the ocean and the atmosphere here. As a result, the sea surface gains the heat from the atmosphere through the positive SH, LH, and LW anomalies, eventually the MLD shoals.
Cheng Li, Zhang Zhiwei, Zhao Wei, et al. 2015. Temporal variability of the current in the northeastern South China Sea revealed by 2.5-year-long moored observations. J Oceanogr, 71(4): 361-372
|
Chu P C, Edmons N L, Fan Chenwu. 1999. Dynamical mechanisms for the South China Sea seasonal circulation and thermohaline
|
Da Silva A M, Young C C, Levitus S. 1994. Atlas of Surface Marine Data 1994, vol. 1, Algorithms and Procedures. NOAA Atlas NES-DIS 6, US Dept of Comm., Washington, D C
|
Dale W L. 1956. Wind and drift currents in the South China Sea. The Malaysian Journal of Tropical Geography, 8: 1-31
|
Ducet N, Le Traon P Y, Reverdin G. 2000. Global high resolution mapping of ocean circulation from TOPEX/Poseidon and ERS-1/2. J Geophys Res, 105(C8): 19477-19498
|
Egbert G D, Bennett A F, Foreman M G G. 1994. TOPEX/POSEIDON tides estimated using a global inverse model. J Geophys Res, 99(C12): 24821-24852
|
Ezer T, ArangoH, Shchepetkin A F. 2002 Developments in terrain fol-lowing ocean models: Intercomparisons of numerical aspects, Ocean Modelling, 4: 249-267
|
Fang Guohong, Dwi Suanto, Sugiarta Wirasantosa, et al. 2010. Volume, heat and freshwater transports from the South China Sea to Indonesian Seas in the boreal winter of 2007-2008. Journal of Geophysical Research, 115(C12): doi: 10.1029/2010JC006225
|
Fang Guohong, Fang Wendong, Fang Yue, et al. 1998. A survey of studies on the South China Sea upper ocean circulation. Acta Oceanographica Taiwanica, 37(1): 1-6
|
Fang Wendong, Fang Guohong, Shi Ping, et al. 2002. Seasonal struc-tures of upper layer circulation in the southern South China Sea from in situ observations. Journal of Geophysical Research, 107(C11): doi: 10.1029/2002JC001343
|
He Yinghui, Cai Shuqun, Wang Dongxiao, et al. 2015. A model study of Luzon cold eddies in the northern South China Sea. Deep-Sea Research Part I: Oceanographic Research Papers, 97: 107-123
|
Hu Jianyu, Kawamura H, Hong Huasheng, et al. 2000. A review on the currents in the South China Sea: Seasonal circulation, South China Sea warm current and Kuroshio Intrusion. Journal of Oceanography, 56(6): 607-624
|
Levitus S. 1982. Climatological atlas of the world ocean. NOAA Prof Paper No 13. US Govt Printing Office, 173
|
Li Li, Wu Risheng, Guo Xiaogang. 2000. Seasonal circulation in the South China Sea: A TOPEX/Poseidon satellite altimetry study. Haiyang Xuebao (in Chinese), 22(6): 13-26
|
Lin Xiayan, Dong Changming, Chen Dake, et al. 2015. Three-dimen-sional properties of mesoscale eddies in the South China Sea based on eddy-resolving model output. Deep-Sea Research Part I: Oceanographic Research Papers, 99: 46-64
|
Mellor G L, Oey L Y, Ezer T. 1998. Sigma coordinate pressure gradi-ent errors and the seamount problem. J Atmos Oceanic Techn-ol, 15(5): 1122-1131
|
Mellor G L, Ezer T, Oey L Y. 1994. On the pressure gradient conun-drum of sigma-cordinate ocean models. J Atmos Oceanic Tech-nol, 11(4): 1120-1129
|
Mellor G L, Yamada T. 1982. Development of a turbulence closure model for geophysical fluid problems. Rev Geophys, 20(4): 851-875
|
Qiao Fangli, Chen Shunnan, Li Chengxi, et al. 1999.The study of wind, wave, current extreme parameters and climatic charac-ters of the South China Sea. Journal of Marine Technology, 33(1): 61-68
|
Qiao Fangli, Yang Yongzeng, Xia Changshui, et al. 2008. The role of surface waves in the ocean mixed layer. Acta Oceanologica Sin-ica, 27(3): 30-37
|
Qiao Fangli, Yuan Yeli, Yang Yongzeng, et al. 2004. Wave-induced mixing in the upper ocean: Distribution and application to a global ocean circulation model. Geophys Res Lett, 31(11): doi: 10.1029/2004GL019824
|
Qiu Bo, Chen Shuiming. 2005. Eddy-induced heat transport in the subtropical North Pacific from Argo, TMI, and altimetry meas-urements. Journal of Physical Oceanography, 35(4): 458-473
|
Roemmich D, Gilson J. 2001. Eddy transport of heat and thermocline waters in the North Pacific:.A key to interannual/decadal cli-mate variability?. Journal of Physical Oceanography, 31(3): 675-687
|
Shan Feng, Qiao Fangli, Lv Xingang, et al. 2009a. A numerical study of the wintertime double-warm-tongue structure in the Huanghai (Yellow) Sea. Acta Oceanologica Sinica, 28(4): 8-15
|
Shan Feng, Qiao Fangli, Xia Changshui. 2009b. numerical study of summertime expansion pattern of Changjiang (Yangtze) River diluted water. Acta Oceanologica Sinica, 28(3): 11-16
|
Shaw P T, Chao S Y. 1994. Surface circulation in the South China Sea. Deep-Sea Res I, 41(11-12): 1663-1683
|
Wang Dongxiao, Liu Yun, Qi Yiquan, et al. 2001. Seasonal variability of thermal fronts in the northern South China Sea from satellite data. Geophys Res Lett, 28(30): 3963-3966
|
Wang Lei, Gan Jianping. 2014. Delving into three-dimensional struc-ture of the West Luzon Eddy in a regional ocean model. Deep-Sea Research Part I: Oceanographic Research Papers, 90: 48-61
|
Wang Qiang, Zeng Lili, Zhou Weidong, et al. 2015. Mesoscale eddies cases study at Xisha waters in the South China Sea in 2009/2010. J Geophys Res Oceans, 120(1): 517-532
|
Wyrtki K. 1961. Scientific results of marine investigations of the South China Sea and the Gulf of Thailand 1959-1961. Naga Report, 2: 164-169
|
Xia Changshui, Qiao Fangli, Yang Yongzeng, et al. 2006. Three-di-mensional structure of the summertime circulation in the Yel-low Sea from a wave-tide-circulation coupled model. Journal of Geophysical Research, 111(C11): doi: 10.1029/2005JC003218
|
Yang Haijun, Liu Qinyu, Liu Zhengyu, et al. 2002. A general circula-tion model study of the dynamics of the upper ocean circula-tion of the South China Sea. Journal of Geophysical Research, 107(C7): 22-1-22-14
|
Yang Yongzeng, Qiao Fangli, Zhao Wei, et al. 2005. MASNUM ocean wave model in sphereical coordinate and its application. Acta Oceanologica Sinica, 27(2): 1-7
|
Yuan Yeli, Hua Feng, Pan Zengdi, et al. 1991. LAGDF-WAM numeric-al wave model. Acta Oceanologica Sinica, 10(4): 483-488
|
Yu Weidong, Qiao Fangli, Yuan Yeli, et al. 1997. Numerical modeling of wind and waves for Typhoon Betty (8710). Acta Oceanolo-gica Sinica, 16(4): 459-473
|
Yu Yongqiang, Liu Hailong, Lin Pengfei. 2012. A quasi-global 1/10° eddy-resolving ocean general circulation model and its prelim-inary results. Chinese Science Bulletin 57(30): 3908-3916
|
Zhang Zhiwei, Zhao Wei, Tian Jiwei, et al. 2013. A mesoscale eddy pair southwest of Taiwan and its influence on deep circulation. J Geophys Res Oceans, 118(12): 6479-6494
|
Zhao Chang, Wang Gang, Qiao Fangli, et al. 2015. A numerical invest-igation into the long-term behaviors of Fukushima-derived 137Cs in the ocean. Acta Oceanologica Sinica, 34(12): 37-43
|
Zu Tingting, Gan Jianping, Erofeeva S Y. 2008. Numerical study of the tide and tidal dynamics in the South China Sea. Deep-Sea Re-search Part I: Oceanographic Research Papers, 55(2): 137-154
|
Category | Year |
EP El Niño | 1982/1983, 1991/1992, 1997/1998, 2006/2007, 2015/2016 |
CP El Niño | 1990/1991, 1994/1995, 2002/2003, 2004/2005, 2009/2010, 2014/2015, 2018/2019, 2019/2020 |