www.nature.com/scientificreports open Effects of modes of climate variability on wave power during boreal summer in the western North Pacific 1 2* Sinil Yang & Jai-Ho oh This study investigates extremes of wave climate in the western North Pacific (WNP) as significant responses to modes of climate variability: the El Niño-Southern Oscillation (ENSO) and the Pacific Decadal Oscillation (PDO). Few studies have explicitly investigated significant wave height in this context, and hence, the aim of the present study is to investigate and quantify the responses to a simulated wave data set over the WNP from 1979–2009 by focusing on the combined influences of the ENSO and PDO during the boreal summer. We conducted a composite analysis of sea surface temperatures, sea-level pressure, and extreme anomalies of wave power density (Pw) on different phase combinations of the ENSO and the PDO, and also analyzed the effects of a latitudinal shift of the ITCZ for composite samples with respect to simulated tropical storm (TS) activities. The results demonstrate that the ENSO played a primarily positive role in intensifying anomalous wave climate, while the PDO had the opposite effect. The responses of the peak wave-period were linked to a strengthened anomalous low-pressure and a cooling of sea surface temperature anomalies. The PDO played a significant role in strengthening or weakening of the effects of the ENSO on Pw, thus confirming the findings of previous studies. We found that responses were dependent on whether ENSO and PDO were in or out of phase. These responses can be described by a strengthening of the southeast trade winds that blow across the equator with respect to a latitudinal shift of the Intertropical Convergence Zone (ITCZ). Our findings contribute to the understanding of a relationship between modes of climate variability and TS activities with respect to the status of the ITCZ over the WNP, which can be relevant factors in the lifetime of wave power and related wave parameters in the WNP during the boreal summer. Extreme wave climates are crucial to the development of various natural processes . It is essential to forecast and reproduce regional wave characteristics in the ocean. Additionally, the ocean wave has the highest energy density in coastal marine areas . Thus, considerable research on the regional wave climate has been conducted to evaluate and characterize the wave power energy potential. Specifically, and with reference to climate variability, previous studies have attempted to improve the understanding between significant wave height (hereaer H ft s) and seasonal 3–5 1,6–11 12,13 oscillations on global and regional scales , and for future projections . The effects of long-term wave 14,15 14 climate variability in the western North Pacific (WNP) have also been reported . Sasaki (2012) evaluated the regional wave climate based on 30 years of observational data and identified an significant positive trend that was induced by more frequent swells. Allan and Komar (2000) reported an increasing Hs trend over the eastern North Pacific (ENP) over approximately 25 years previous to their study, which may have been linked to the El Niño-Southern Oscillation (ENSO). Additionally, because the ENSO modulates tropical storm (TS) activities in 16,17 the North Pacific , TS activities in the WNP− in addition to the ENSO− also have a substantial inu fl ence on 1,18,19 the wave climate. Previous studies have demonstrated that the number of TS could be essential for ae ff cting variations in the Hs in the WNP during the boreal summer, which is linked to warm sea surface temperature (SST) anomalies in the Niño 3.4 region [5°S–5°N, 170°–120°W]. However, few studies have explicitly investigated Korea Ocean Satellite Center, Korea Institute of Ocean Science and Technology, Busan, 49111, Republic of Korea. Department of Environmental Atmospheric Sciences, Pukyong National University, Busan, 48513, Republic of Korea. *email: firstname.lastname@example.org Scientific Repo R t S | (2020) 10:5187 | https://doi.org/10.1038/s41598-020-62138-0 1 www.nature.com/scientificreports www.nature.com/scientificreports/ Figure 1. Linear regression spatial patterns of standardized sea surface temperature (SST, left) and sea-level pressure (SLP, right) onto (a,b) the El Niño-Southern Oscillation (ENSO), and (c,d) the Pacific Decadal Oscillation (PDO) indices during June, July, and August (JJA) from 1979–2009. The analysis domain is shown as a black box covering [20°S–70°N, 90°E–150°W]. The hatching represents grid points with significant regression coefficients at the 5% level. wave power density (hereafter Pw), which is defined by a function of the Hs and peak wave period (hereafter 20,21 6,22 Tp) . A longer Tp produces a greater wave run-up, which raises flooding potential . Therefore, Pw is more useful than Hs alone for describing the magnitude of storm effects, which depends on wave energy brought by 23,24 long waves . 9,25 However, previous studies did not address the role of the intertropical convergence zone (ITCZ) and trade winds on extremes of Pw in the WNP during JJA. Consistent trade winds (e.g., northeasterly and southeasterly) are among the tropical climate phenomena. The trade winds regularly blow towards the equator, where extreme precipitation occurs in the ITCZ. However, the ITCZ varies by seasonal latitudinal shifts. Richter et al . explained that the ITCZ moves approximately 8°N from its southernmost latitude during the boreal spring to its northern- most latitude during the boreal summer, which is caused by the strengthening of the southeast trade winds (i.e., 27,28 the equatorial easterlies) blowing across the equator from the southern hemisphere. Previous studies have demonstrated that the latitudinal shift of the ITCZ ae ff cts SST warming to the north of the equator. Furthermore, Vincent et al. suggested that the influence of the ENSO along the ITCZ appears to be associated with intensi - fied TS activities during El-Niño. In particular, the wind variability in the WNP is dominated by subtropical, large-scale atmospheric circulation and ae ff cts the regional wave climate during the boreal summer (i.e., June, July, and August (JJA)). For instance, Tao et al. revealed that the tropical, large-scale atmospheric circulation was accompanied by the development of anomalous cyclones in the WNP. Furthermore, Zhang et al. found a relationship between the central Pacific (CP) cooling and the westerlies (i.e., the Rossby wave) response to the part of the northwestern CP . The aim of the present study is to investigate and quantify the mean distribution and extreme trends of a simulated wave data set over the WNP from 1979–2009, by mainly focusing on the combined influences of the ENSO and PDO during the boreal summer. To investigate influencing factors, we conducted a composite analysis of SST, sea-level pressure (SLP), and extreme anomalies of Pw on different phase combinations of the ENSO and PDO. In particular, we analyze the ee ff cts of a latitudinal shi ft of the ITCZ for each composite sample with respect to TS activities. Results Influence of modes of climate variability. SST and SLP. Figure 1 illustrates the inu fl ence of each mode of climate variability on both SST and SLP as defined by linear regression analysis. The ENSO can be seen as a principal mode among the modes of climate variability for both SST and SLP in the tropical Pacific (Fig. 1a,b). e Th most significant relationship (p < 0.05) between the ENSO and SSTs with a regression slope larger than 0.5 was found over the eastern Tropical Pacific (Fig. 1a). Furthermore, the effects of the ENSO were found to reach far beyond the tropical Pacific, consistent with the findings of Hu et al . . In particular, statistically significant (hereaer a ft t the 5% level) responses of the ENSO to SSTs were recognized from the CP to the eastern Pacific (EP) and the tropical Indian Oceans, where warm conditions prevail, whereas cold conditions prevail in the western tropical Pacific. Additionally, we found a significant relationship (p < 0.05) between the PDO and SSTs with a positive sign in the mid-Pacific Ocean and a negative sign in the WNP (Fig. 1c). Furthermore, the ENSO and the PDO had a substantial impact on the SLP in the anomalous low-pressure region of the northeastern Pacific, and the anomalous high-pressure regions in the CP and Indian Oceans (Fig. 1b,d). Thus, anomalous pressure with a Scientific Repo R t S | (2020) 10:5187 | https://doi.org/10.1038/s41598-020-62138-0 2 www.nature.com/scientificreports www.nature.com/scientificreports/ th Figure 2. Linear regression spatial patterns of standardized 99 percentile of wave parameters: significant wave height (Hs , left), peak wave period (Tp , center), and wave power density (Pw , right) onto (a–c) the El Niño- 99 99 99 Southern Oscillation (ENSO), and (d–f) the Pacific Decadal Oscillation (PDO) indices during June, July, and August (JJA) from 1979–2009. The hatching represents grid points with significant regression coefficients at the 5% level. th sustained wind speed could lead to the 99 percentile of wave parameters: Hs and Pw . The impact of the PDO 99 99 was found to be related to SST warming over the northeastern North Pacific, SST cooling over the northwestern North Pacific, and anomalous low-pressure over the central and northeastern North Pacific (Fig. 1c,d). Extremes of wave climate. Figure 2 presents the spatial regression maps for extremes of the wave climate, which were regressed onto the modes of climate variability. The warm ENSO phase typically exhibited a positive response to wave parameters in the WNP, which was contrary to the influence of the PDO. The influence of the warm ENSO phase on the Pw was statistically significant, and was found to be mainly a positive response in the WNP and a negative response in the CP (Fig. 2c). The strongest statistically significant responses of the Pw were located in the WNP along the 20°N latitude. The positive influences of the ENSO along the region seem to be cor - related with enhanced TS activities during El Niño. This is because increased, sustained wind-speed was a dom- inant factor causing increased Pw , which is consistent with the findings of Yang and Oh (2017) . The response of the PDO to the Pw was approximately opposite to the ENSO (Fig. 2c,f ), but few responses were significant at the 5% level. Finally, for the positive PDO phase, the regions that exhibited negative responses were located in the WNP, which spans the western coast of the WNP to the mid-Pacific and relates to the storm tracks that are developed in the mid-latitudes of the WNP. The positive PDO phase had a significant negative impact on the Pw over the WNP, which covers almost the same region as that influences by the ENSO; hence, there was an inverse relationship (i.e., out of phase) to the ENSO influence in the WNP (Fig. 2a,d). Additionally, the statistically significant areas (p < 0.05) of the ENSO and the PDO influences on the Hs were mainly similar to the response of the Pw because the Hs is a dominant factor when calculating the Pw in Eq. (1) (see Methods). Figure 2a presents a positive response of the Hs over an area that spans from the WNP to the tropical Pacific along 20°N latitudes, and the impact of changes in SST anomalies in the ENP. Contrary to the Hs , the response of the warm ENSO phase on the Tp was higher in the WNP and lower in the central North 99 99 Pacific (Fig. 2b). Most of the positive responses to the warm ENSO phase were found to exist along the western basin of the WNP (e.g., the East China Sea and the Philippine Sea). The notable regions imply a substantial relationship, whereby changes to tropical SST anomalies due to the warm ENSO phase resulted in variations in anomalous long Tp . On the other hand, the response to the positive PDO phase was nearly negative across the WNP and of lesser significance to that of the warm ENSO phase at the 5% level, which was contrary to the ENSO response (Fig. 2b,e). Overall, the northwestern WNP region with a significant response of the Tp to the ENSO was correlated with the strengthening of the anomalous low-pressure (high-pressure) and cooling (warming) of the WNP in in-phase (out-of-phase) between ENSO and PDO. Composite analysis. e a Th forementioned analysis of the effects of the chosen climate indices on the Pw was achieved by applying the influence of individual climate indices to each other. To further investigate the combined influences of the dominant climate mode in the WNP, we conducted a composite analysis of the SST, SLP, and Pw anomalies on different phase combinations of the ENSO and the PDO. In particular, our analyses aimed to determine the effects of TS activities near the ITCZ between the ENSO and the trade wind index (TWI). Scientific Repo R t S | (2020) 10:5187 | https://doi.org/10.1038/s41598-020-62138-0 3 www.nature.com/scientificreports www.nature.com/scientificreports/ Composites Sample years Mean TWI Number of TS Mean origin location +ENSO/+PDO 1982, 1987, and 1997 1.30 34 20. 2°N, 142. 1°E +ENSO/−PDO 1991 and 2002 0.79 23 21. 0°N, 133. 7°E −ENSO/+PDO 1988, 1998, and 2007 −0.99 21 22. 8°N, 136. 3°E −ENSO/−PDO 1999 and 2000 −1.25 10 20. 6°N, 132. 8°E Table 1. Classification and statistics for tropical storm (TS) activities of El Niño (+ENSO) or La Niña (−ENSO) and/or positive or negative PDO (+PDO or −PDO) year combinations > ±1.5 standard deviation occurring in June, July and August (JJA) from 1979–2009. The mean trade wind index (TWI) is the average of the normalized TWI for each case. The number of TS is the sum of the number of TS for all years in each composite. The mean origin location is the averaged location for all TS genesis in each composite. Figure 3. Composite patterns of sea surface temperature (SST, shaded area, in degrees Celsius] and sea-level pressure (SLP, contours, in hPa) anomalies for the four combinations of the El Niño-Southern Oscillation (ENSO) and the Pacific Decadal Oscillation (PDO), which are the (a) El Niño (+ENSO)/+PDO, ( b) +ENSO/− PDO, (c) La Niña (−ENSO)/+PDO, and (d) −ENSO/−PDO years. All anomalies are computed relative to the 31-year period from 1979–2009. Filled and unfilled dots denote significant SLP and SST anomalies at a confidence level of 90%, respectively. Inu fl ence of ENSO and PDO. We selected four composites from El Niño (hereaer ft + ENSO) and La Niña (here- aer ft − ENSO) events based on the Ocean Niño Index (ONI) indices, and positive and negative phases (+PDO and − PDO) of the PDO for comparison to the ±0.5-standard deviation. Table 1 shows that the classified sam- ple year was 3, 2, 3, and 2 for each of the different phase combinations. Figures 3 and 4 show the SST, SLP, and Pw anomalies for the four composites of both the ENSO and the PDO. El Niño and La Niña events have been 34,35 found to present typical SST anomalies in the equatorial tropical Pacific . SST anomalies in the CP and EP are significantly warming during the +ENSO/+PDO (Fig. 3a) and significantly cooling during the −ENSO/−PDO (Fig. 3d). These results indicate strengthening of the ENSO impact when the PDO is in phase with the ENSO. When the ENSO and PDO are out of phase, we found that the spatial patterns of SST were described by the weakening of the ENSO impact, mainly over the CP (Fig. 3b,c). In terms of the influence on SLP, similarly to SST, anomalously low SLPs exist during +ENSO/+PDO and anomalously high SLPs exist during − ENSO/ − PDO over the mid-central North Pacific in comparison to years with other phases (Fig. 3a,d). The results were clearly significant and confirmed that the PDO presents a substantial role in strengthening or weakening the ENSO 9,36 effects on SST and SLP over the North Pacific, as previously reported . Furthermore, Pw anomalies were found to increase during +ENSO/+PDO and decrease during −ENSO/− PDO around the WNP (Fig. 4). In this region, an +ENSO (−ENSO) related to +PDO (−PDO) was observed to cause substantial significant increases Scientific Repo R t S | (2020) 10:5187 | https://doi.org/10.1038/s41598-020-62138-0 4 www.nature.com/scientificreports www.nature.com/scientificreports/ Figure 4. Composite patterns of Pw [shaded, in kW/m] anomalies for the four combinations of the El Niño-Southern Oscillation (ENSO) and the Pacific Decadal Oscillation (PDO), which are the (a) El Niño (+ENSO)/+PDO, (b) +ENSO/−PDO, (c) La Niña (−ENSO)/+PDO, and (d) −ENSO/−PDO years. The blue arrows indicate the significant wind vector anomalies at 850 hPa (vectors, in m/s) at a confidence level of 80%. All anomalies are computed relative to the 31-year period from 1979–2009. The black cross and solid red line represent the TS genesis location and the ITCZ area detected by the ICON model for each case, respectively. The black dot denotes the significant Pw anomalies at a confidence level of 90%. (non-significant decreases) in Pw anomalies. The SST in the central North Pacific was significantly cooling for +ENSO/+PDO events and significantly warming for + ENSO/−PDO events (Fig. 4a,b). These results indi- cate that the −PDO strengthening, when combined with +ENSO events, ae ff cted the sign and peak region of SST anomalies in the CP. Additionally, the significant intensity of the SST and SLP anomalies in the WNP also depends on the phase of the PDO for −ENSO/+PDO and −ENSO/−PDO events (Fig. 3c,d). In the case of −ENSO years, the Pw anomalies in the WNP also represented a negative sign regardless of the PDO phase, thus indicating that the ENSO played a substantial role in strengthening or weakening the Pw in the WNP (Fig. 4c,d). However, the ENSO had the most significant influence on the Pw during +ENSO/+PDO events (Fig. 4a), Pw anomalies during +ENSO/−PDO events were mostly non-significant in the WNP (Fig. 4b). This relationship means that the ENSO did not have a direct impact on the Pw in the WNP, rather an indirect impact. We observed that this indirect impact may be associated with the zonal trade wind and enhanced TS activity, especially during the boreal summer, along with the ITCZ. This relationship is analyzed in the following section. Inu fl ence of ENSO and TS activities. e infl Th uences of the ENSO along the WNP region appear to be correlated with enhanced zonal trade winds because increased zonal wind-speed energy was a dominant factor in inducing changes in the Pw (see Fig. 2a,c). The positive influences of the ENSO, along with the ITCZ, seem to be related to TS activities, thus indicating that increased surface wind-speed energy indeed induced an increased Pw in the WNP. Forecasting properties of TS are important for gaining a better understanding of the modes of climate var- iability. There might be uncertainty because different analyses have used different data assimilation schemes and 37–39 global climate models (GCMs) presented the impacts of ENSO on TS activity using climate models. There might be uncertainty because different analysis uses a different data assimilation scheme and global climate mod- els (GCMs), respectively . The Icosahedral Nonhydrostatic model (ICON) model has been verified by several idealized simulations, and has been used for the operational production of the German Weather Service (DWD) since 2015 . The use of a simulated wave data set of the type used by the Atmospheric Model Intercomparison Project (AMIP) to detect TS presents an opportunity to reduce the uncertainty by supporting prescribed SSTs and sea ice concentrations (SICs) when simulating the atmospheric results of the ICON. In particular, the ITCZ movement varies depending on the phases of ENSO and PDO. Thus, the geographical findings of the ITCZ can be important in interpreting the relationship between the waves and the modes of climate variability. Therefore, we have employed the objective identification of the ITCZ in Berry and Reeder (2014) . To identify the ITCZ Scientific Repo R t S | (2020) 10:5187 | https://doi.org/10.1038/s41598-020-62138-0 5 www.nature.com/scientificreports www.nature.com/scientificreports/ from numerical simulation, we computed the mean divergence between 1000 and 850 layers, the gradient of the divergence, and the wet-bulb potential temperature at 850 hPa, which are obtained from the ICON model. In the present study, we apply the same thresholds of the original method except the Laplacian of the divergence. Through this objective method, we have provided a simulated ITCZ area in each phase between ENSO and PDO. e IT Th CZ area represents solid red lines in each panel in the revised Figure 4. e AMIP Th -typed simulation of the ICON can help to bring more confidence to TS detection used in the pres- ent study. Thus, we detected simulated TS (ranging from tropical storms to typhoons) in the ICON model using the TC tracking method of Murakami and Sugi (2010) . From the composites of Fig. 4, we investigated the rela- tionship between the Pw and SST anomalies with TS activities in JJA (Fig. 4). Over the tropical Pacic fi , the zonal SST gradient between cold water in the eastern part and warm water in the western part, was considered to be the cause of this thermally induced direct circulation, which is called the Walker circulation . The normalized zonal wind index at 850 hPa, TWI, was used to describe the behavior of the near-surface wind over the WNP under the Walker circulation along the equator. Table 1 shows the number and mean origin location of TS for each of the dier ff ent phase combinations, and the TWI. During El Niño years, the westerly zonal trade winds were strong and persistent and the mean TWI was ~ 1.3. The genesis location of the TS shifted eastward across 150°E in the WNP during El Niño years in comparison to La Niña years (Fig. 4a–d). More TS occurred during an El Niño event in the eastern region of the WNP (150°–180°E) than in the western region (120°–150°E). The eastward shift in gene- sis locations was more pronounced during El Niño events, and most of the formation origin locations were found as far east as 110°–180° E along the 20°N latitude (Fig. 4a). Additionally, we found that TCs have a better chance of interacting with synoptic systems in the midlatitudes because of their longer lifespans. The finding is consistent with the results of Chu (2004) , where the eastward and equatorward shift in genesis locations allowed the TS to survive longer during their movement toward the East Asian continent. Most of the origin locations of TS for La Niña events were confined to the west of 15°E during the summer (Table 1). The genesis locations during the sum- mer remained at the midlatitudes (20°–30°N) in the WNP, which was due to the easterly trade winds toward the western Pacific where there are warm pools . In particular, Figure S2 shows the composite means of Tp anomalies and its significance at a confidence level of 90% (see Supplementary Fig. S2). The ENSO and PDO contribute to a magnitude of Tp of tropical waves caused by trade winds blowing from the southern hemisphere. This mechanism is related to the propagation of large swell waves from near the equator to the northern region of the WNP (Fig. S2). This figure helps us understand the wave climate of the North Pacific that includes wind-driven waves and swells. Therefore, substantial variations in the Pw extremes were ae ff cted by TS caused by tropical SST anomalies and anomalous zonal trade winds (Fig. 4). The status of the ITCZ could be relevant to the anomalies of the Pw in the WNP, which is associated with atmospheric circulations in the tropical Pacific. Discussion In summary, several dominant features (e.g., tropical storms, trade winds, monsoons, and modes of climate var- iability) contribute to the wave climate. The contributions include local wind waves ae ff cted by ITCZ movement, swells generated near the equator, and tropical storm activities. Most extreme waves, however, may be related to TS activities, though large swells may also cause some from ITCZ. Several studies have confirmed that the PDO 9,36 plays an important role in the impact of ENSO on wave climate . es Th e findings imply that ENSO may be the sole source for modulating the ee ff cts of low-level wind on the wave climate near ITCZ. This ITCZ movement depends on the different phases of ENSO and PDO. Extremes of the wave climate in the WNP exhibited considerable responses to the ENSO, the PDO, and the northward migration of the ITCZ. The ENSO was found to be the dominant mode of climate variability both in terms of SST and SLP in the tropical Pacific. In particular, the statistically significant areas for the influences of the ENSO and the PDO on the Hs were found to be mainly similar to the response on the Pw , which is because 99 99 the Hs is a dominant factor when calculating the Pw. The notable regions of Tp imply the significant relationship, whereby changes to tropical SST anomalies due to the ENSO result in substantial variations in anomalously long Tp . On the other hand, the response of the Tp to the PDO is almost negative across the WNP and less significant to that of the ENSO. Overall, the statistically significant responses of the Tp were linked to a strengthening of the anomalous low-pressure in the WNP and the cooling of the WNP. In particular, the extremes of the wave param- eters in the WNP could occur as either the warming and cooling of SST or the anomalous high and low-pressure accompanied by a sustained wind. Additionally, we found that extreme wave climates were a response to whether the ENSO and PDO were in phase or out of phase. Thus, our results confirmed that the PDO plays a statisti- cally significant role in strengthening or weakening the ENSO influences on the Pw in the WNP, as previously 9,36 47 reported . Furthermore, Camargo et al. found that SST anomalies in the tropical ocean that were induced by ENSO could alter the atmospheric conditions in the WNP (e.g., atmospheric stability, zonal trade wind, and 16,17 vertical wind shear). Since the ENSO modulates TS activities in the North Pacific , the TS activities in the 1,18,19 WNP, as well as the ENSO, have a substantial impact on the wave climate. Previous studies have shown that the number of TS created could be essential for ae ff cting the variation of the Hs in the WNP during the boreal summer, which was linked to warm SST anomalies in the Niño 3.4 region. In addition to the relevant factors, as mentioned above (e.g., ITCZ movement and southeast trade wind), several studies stated changes in SST, monsoon, and low-level relative vorticity contributed to changes in the fre- 48–50 quency of TSs in the WNP during ENSO events . This change implies that the TS activities can contribute to the wave climate in the WNP. For example, the positive influence of the ENSO appears to be related to enhanced TS activities during El Niño, and indicate that increased wind-speed energy induces increased extremes of Pw, as 29 51 a previous study suggested . In particular, Chan et al. examined the changes in the early summer monsoonal circulation over the South China Sea can be described as the change in the extension of the 500-hPa subtropical high depending on the ENSO coupled with the PDO. The variability of the monsoon trough can be also one of the relevant factors that influence the wave climate in the region of the WNP (e.g., the South China Sea). Further Scientific Repo R t S | (2020) 10:5187 | https://doi.org/10.1038/s41598-020-62138-0 6 www.nature.com/scientificreports www.nature.com/scientificreports/ research may be needed in order to investigate the impact on the wave climate by modulating the monsoonal circulation. Similarly, combined effects of extreme sea levels due to the contribution of waves, tides, and pressure systems lead to extreme events (e.g. coastal inundation) . The main causes of storm surge in the coastal region of the Pacific are from tropical cyclones (TCs), and several studies have investigated how they might change . For example, the results of Wang and Zhou (2016) found the relationship between the important trends in extreme sea levels and the influence of tropical cyclones in Macau and Hong Kong. u Th s, this relationship between ENSO and TCs will help us understand the changes in sea-level height and associated wave parameters for future pro- jection in the WNP. From the composite analysis in this study, we found that the SST and SLP anomalies related to the positive phase of the ENSO and PDO could be relevant to anomalies of extreme wave climate in the WNP. Also, the changes in the wave climate were related to the conditions of equatorial SST and low-level wind fields over the tropical Pacific as well as the TS activities. The response can be explained by a strengthening of the southeast trade winds (i.e., the equatorial easterlies) blowing across the equator in the southern hemisphere with respect to a northward shift of the ITCZ. These Northward shifts of the ITCZ can significantly increase Pw anomalies by increasing Tp anomalies in the WNP, i.e., the ITCZ contributes to the formation of large swells that propagate from the equator to the North Pacific. Thus, the present study contributes to the understanding of the (i) effects of the ENSO and the PDO on extreme wave climate, and (ii) relationship between the ENSO, the PDO, and TS activ- ities with respect to latitudinal shifts of the ITCZ over the WNP. Therefore, the present findings can be important in interpreting the relationship between waves and modes of climate variability. The present study also indicates that the status of the ITCZ in the positive phase of ENSO and PDO can be a relevant factor in the lifetime of extreme wave climates in the WNP during the boreal summer. Data and Methods Data. Global wave data set. Following Yang and Oh (2017) , the present study used the same simulation output in terms of Hs and Tp as a wave parameters. The global wave data set is available from January 1979 to December 2009 on regular latitude-longitude grids at 0.25° × 0.25° resolution at three-hour intervals. The grid - ded wave data set provides commonly used wave parameters applicable to the Hs, mean wave height, Tp, mean wave period, mean and peak wave direction, and wind sea wave period and direction. Yang and Oh (2017) used the full-spectral third-generation numerical wave model WAVEWATCH-III . Their study employed the model to reproduce the global wave characteristics with the AMIP-typed hindcast from an atmospheric general circu- lation model (AGCM); ICON. The ICON model is a joint development project of the Max Planck Institute for Meteorology (MPI-M) and the DWD . A detailed description of the ICON can be found as Supplementary text S1 online. Based on a procedure defined elsewhere , the AMIP-type simulation was run with the prescribed SSTs and SICs observed by the National Center for Atmospheric Research (NCAR) to avoid an unrealistic condition of both SSTs and SICs. e Th Hs among the wave parameters was validated by a 30-year wave hindcast of Climate Forecast System Reanalysis and Reforecast (CFSRR) , in which the area-weighted mean root mean squared error (RMSE) and correlation coefficient for the WNP (100°E–180°E, 0°–60°N) were 0.5 m and 0.69, respectively . This CFSR hindcast provides three hourly global data with a spatial resolution of 0.5°, which is considered a reliable dataset among global hindcasts. Additionally, validation of the simulation results for Hs and Tp was also in agree- ment with the CFSR hindcast (see Supplementary text S2 and Fig. S1). The simulation results, therefore, are useful in the analysis domain and allows us to employ climate research. Ancillary data set. Interannual climate variability in the WNP is partially induced by the behavior and strength of ocean-atmospheric oscillations. We used two climate indices: the PDO and the ENSO, which have 20–30 year 57 58 cycles and 2–7 year periods . Active phases of the PDO and ENSO have brought warmer and drier winters to 59,60 the WNP region, while negative aspects have brought colder and wetter conditions . Therefore, a couple of climate indices, the ONI and the PDO, were used to describe the climate variability of the ENSO and the PDO, respectively. To evaluate the relationship between the simulated wave climate and global climate variability, we obtained the following climate indices during JJA: (i) the ONI was obtained from the United States National Weather Service Climate Prediction Center (https://www.cpc.ncep.noaa.gov/data/indices/), derived as the three-month running average of SST anomalies over the Niño 3.4 region covering [5°S–5°N, 170°–120°W], and (ii) the PDO index was downloaded from the Joint Institute for the Study of the Atmosphere and Ocean (http:// research.jisao.washington.edu/pdo) computed as the principal component of monthly SST anomalies in the spe- cific domain (poleward of 20°N, from 20°–70°N, 110°E–100°W). We also used the global SST data, obtained from the Program for Climate Model Diagnosis and Intercomparison (https://pcmdi.llnl.gov/mips/amip/amip2) to identify the influence on the extremes of the Pw. 20,21 Methodology. By referring to the definition and calculation method described elsewhere , the Pw for an irregular sea-state was derived from the wave spectral parameters using Eq. (1): ρg P =× TH × () wp s 64φ (1) where ρ is the water density, g is the gravitational acceleration with 9.8 m/s, T is the peak wave period, and H is p s the average of the highest third of the waves in a time-series of waves. Eq. (1) estimates the Pw, which increases nonlinearly with H and is proportional to T . The T is sensible to the wavelength and duration of the wave fetch s p p over which strong winds persist. Estimating P helps us to consider the T as well as the H because the effects w p s of storms on coastal processes depend on the wave power energy that is propagated by the waves. Therefore, Scientific Repo R t S | (2020) 10:5187 | https://doi.org/10.1038/s41598-020-62138-0 7 www.nature.com/scientificreports www.nature.com/scientificreports/ contrary to H alone, the Pw provides useful information on wave effects. The summer wave parameters (i.e., Hs, Tp, and Pw) in the WNP covering [20°S–70°N, 90°E–150°W] were then calculated from the wave data set at three-hourly intervals. The analysis domain includes parts of the Southern Hemisphere and ENP. Moreover, the th 99 percentile of wave parameters were defined to be the extremes of the wave parameters (e.g., Hs ). The spatial distribution of climatological averaged wave parameters was computed to identify the regional characteristics during JJA. We examined the responses of SST, SLP, and wave parameters on the ENSO and the PDO from 1979 to 2009 during JJA by using the linear regression method. In the regression analysis, all variables were standardized by the distribution mean and standard deviation for each variable as x′= () x − μ /σ, where x is the original variable, x is the standardized variable, μ is the mean of that variable, and σ is its standard deviation. For SLP, the unit of −1 measurement of the slope of regression is hPa (standard deviation) , i.e., no units. Furthermore, to investigate the combined influences of the dominant climate mode in the WNP during JJA, we conducted a composite anal- ysis of the SST, SLP, and Pw anomalies on different phase combinations of the ENSO and PDO. All anomalies are computed relative to the 31-year period from 1979–2009. To obtain composite samples, we selected several years from +ENSO and −ENSO based on the ONI indices, and +PDO and −PDO for comparison to the stand- ard deviations. We also evaluated the relationship between the Pw anomalies and the TS activities representing the wind-speed energy in the WNP during JJA by considering the averages for each composite. In particular, a detecting method of TS in Murakami and Sugi (2010) was used to detect candidates of TS from 3-hourly raw data in the ICON model. The TWI is derived as latitude-weighted averages of zonal wind anomalies at 850 hPa over the western central Pacific (135°E–180°E, 5°S–5°N), where the positive (negative) values of the TWI indices imply easterly (westerly) anomalies. All atmospheric data sets used in the present study were obtained from the ICON model of Yang and Oh (2017) . Data availability The ONI index is obtained from the Climate Prediction Center (CPC)’s website (https://www.cpc.ncep.noaa. gov/data/indices/) of the National Oceanic and Atmospheric Administration (NOAA). The PDO index is also downloaded from the Joint Institute for the Study of the Atmosphere and Ocean (JISAO)’s website (http:// research.jisao.washington.edu/pdo/PDO.latest). The global SST data for 1979–2009 is obtained from the website (https://pcmdi.llnl.gov/mips/amip/amip2). Received: 1 November 2019; Accepted: 28 February 2020; Published: xx xx xxxx References 1. Yang, S. & Oh, J. H. Long-term changes in the extreme significant wave heights on the western north pacific: Impacts of tropical cyclone activity and enso. Asia-Pacic J fi . Atmospheric Sci. 54, 103–109, https://doi.org/10.1007/s13143-017-0063-y (2018). 2. Saidur, R., Islam, M. R., Rahim, N. A. & Solangi, K. H. A review on global wind energy policy. Renew. Sustain. Energy Rev. 14, 1744–1762, https://doi.org/10.1016/j.rser.2010.03.007 (2010). 3. Espejo, A., Camus, P., Losada, I. J. & Méndez, F. J. Spectral ocean wave climate variability based on atmospheric circulation patterns. J. Phys. Oceanogr. 44, 2139–2152, https://doi.org/10.1175/jpo-d-13-0276.1 (2014). 4. Izaguirre, C., Méndez, F. J., Menéndez, M. & Losada, I. J. Global extreme wave height variability based on satellite data. Geophys. Res. Lett. 38, https://doi.org/10.1029/2011gl047302 (2011). 5. Wang, X. L. & Swail, V. R. Changes of extreme wave heights in northern hemisphere oceans and related atmospheric circulation regimes. J. Clim. 14, 2204–2221, https://doi.org/10.1175/1520-0442 (2001). 6. Bromirski, P. D., Cayan, D. R., Helly, J. & Wittmann, P. Wave power variability and trends across the north pacific. J. Geophys. Res. Ocean. 118, 6329–6348, https://doi.org/10.1002/2013jc009189 (2013). 7. Gleeson, E., Gallagher, S., Clancy, C. & Dias, F. Nao and extreme ocean states in the northeast atlantic ocean. Adv. Sci. Res. 14, 23–33, https://doi.org/10.5194/asr-14-23-2017 (2017). 8. Guo, L., Perrie, W., Long, Z., Toulany, B. & Sheng, J. The impacts of climate change on the autumn north atlantic wave climate. Atmosphere-Ocean 53, 491–509, https://doi.org/10.1080/07055900.2015.1103697 (2015). 9. Kumar, P., Min, S. K., Weller, E., Lee, H. & Wang, X. L. Influence of climate variability on extreme ocean surface wave heights assessed from era-interim and era-20c. J. Clim. 29, 4031–4046, https://doi.org/10.1175/jcli-d-15-0580.1 (2016). 10. Lopez, H. & Kirtman, B. P. Investigating the seasonal predictability of significant wave height in the west pacific and indian oceans. Geophys. Res. Lett. 43, 3451–3458, https://doi.org/10.1002/2016gl068653 (2016). 11. Martínez-Asensio, A. et al. Response of the north atlantic wave climate to atmospheric modes of variability. Int. J. Climatol 36, 1210–1225, https://doi.org/10.1002/joc.4415 (2016). 12. Hemer, M. A., Fan, Y., Mori, N., Semedo, A. & Wang, X. L. Projected changes in wave climate from a multi-model ensemble. Nat. climate change 3, 471–476, https://doi.org/10.1038/NCLIMATE1791 (2013). 13. Shimura, T., Mori, N. & Mase, H. Future projection of ocean wave climate: Analysis of sst impacts on wave climate changes in the western north pacific. J. Clim. 28, 3171–3190, https://doi.org/10.1175/jcli-d-14-00187.1 (2015). 14. Sasaki, W. Changes in wave energy resources around japan. Geophys. Res. Lett. 39, L23702, https://doi.org/10.1029/2012gl053845 (2012). 15. Allan, J. & Komar, P. Are ocean wave heights increasing in the eastern north pacific? Eos, Transactions Am. Geophys. Union 81, 561–567, https://doi.org/10.1029/eo081i047p00561-01 (2000). 16. Jin, F. F., Boucharel, J. & Lin, I. I. Eastern pacific tropical cyclones intensified by el niño delivery of subsurface ocean heat. Nat. 516, 82–85, https://doi.org/10.1038/nature13958 (2014). 17. Toma, V. E. & Webster, P. J. Oscillations of the intertropical convergence zone and the genesis of easterly waves. part i: Diagnostics and theory. Clim. Dyn. 34, 587–604, https://doi.org/10.1007/s00382-009-0584-x (2010). 18. Sasaki, W. Recent increase in summertime extreme wave heights in the western north pacific. Geophys. Res. Lett. 32, L15607, https:// doi.org/10.1029/2005gl023722 (2005). 19. Sasaki, W. & Toshiyuki, H. Interannual variability and predictability of summertime significant wave heights in the western north pacific. J. Oceanogr. 63, 203–213, https://doi.org/10.1007/s10872-007-0022-9 (2007). 20. Cornett, A. M. et al. A global wave energy resource assessment. In The Eighteenth International Offshore and Polar Engineering Conference (International Society of Offshore and Polar Engineers, 2008). Scientific Repo R t S | (2020) 10:5187 | https://doi.org/10.1038/s41598-020-62138-0 8 www.nature.com/scientificreports www.nature.com/scientificreports/ 21. Carballo, R. & Iglesias, G. A methodology to determine the power performance of wave energy converters at a particular coastal location. Energy Convers. Manag. 61, 8–18, https://doi.org/10.1016/j.enconman.2012.03.008 (2012). 22. Stockdon, H. F., Holman, R. A., Howd, P. A. & Sallenger, A. H. Empirical parameterization of setup, swash, and runup. Coast. Eng. 53, 573–588, https://doi.org/10.1016/j.coastaleng.2005.12.005 (2006). 23. Bromirski, P. D. Coastal Flooding-Potential Projections: 2000–2100 (California Energy Commission, 2012). 24. Bromirski, P. D. & Kossin, J. P. Increasing hurricane wave power along the u.s. atlantic and gulf coasts. J. Geophys. Res. Ocean. 113, C07012, https://doi.org/10.1029/2007JC004706 (2008). 25. Bromirski, P. D. & Cayan, D. R. Wave power variability and trends across the north atlantic influenced by decadal climate patterns. J. Geophys. Res. Ocean. 120, 3419–3443, https://doi.org/10.1002/2014jc010440 (2015). 26. Richter, I. et al. Phase locking of equatorial atlantic variability through the seasonal migration of the itcz. Clim. Dyn. 48, 3615–3629, https://doi.org/10.1007/s00382-016-3289-y (2017). 27. Chang, P. et al. Climate uc fl tuations of tropical coupled systems-the role of ocean dynamics. J. Clim. 19, 5122–5174, https://doi. org/10.1175/JCLI3903.1 (2006). 28. Keenlyside, N. S. & Latif, M. Understanding equatorial atlantic interannual variability. J. Clim. 20, 131–142, https://doi.org/10.1175/ JCLI3992.1 (2007). 29. Vincent, E. M. et al. Interannual variability of the south pacific convergence zone and implications for tropical cyclone genesis. J. Clim. 36, 1881–1896, https://doi.org/10.1007/s00382-009-0716-3 (2011). 30. Tao, C., Jiang, H. & Zawislak, J. The relative importance of stratiform and convective rainfall in rapidly intensifying tropical cyclones. Mon. Weather. Rev. 145, 795–809, https://doi.org/10.1175/MWR-D-16-0316.1 (2017). 31. Zhang, Z., Ren, B. & Zheng, J. Leading modes of tropical pacific subsurface ocean temperature and associations with two types of el niño. Sci. reports 7, 42371, https://doi.org/10.1038/srep42371 (2017). 32. Son, H.-Y., Park, J.-Y. & Kug, J.-S. Precipitation variability in september over the korean peninsula during enso developing phase. Clim. dynamics 46, 3419–3430, https://doi.org/10.1007/s00382-015-2776-x (2016). 33. Hu, Z.-Z. et al. An analysis of warm pool and cold tongue el niños: Air-sea coupling processes, global influences, and recent trends. Clim. dynamics 38, 2017–2035, https://doi.org/10.1007/s00382-011-1224-9 (2012). 34. Hoerling, M. P., Kumar, A. & Zhong, M. El niño, la niña, and the nonlinearity of their teleconnections. J. Clim. 10, 1769–1786, https://doi.org/10.1175/1520-0442 (1997). 35. Yeh, S.-W., Kug, J.-S. & An, S.-I. Recent progress on two types of el niño: Observations, dynamics, and future changes. Asia-Pacic J fi . Atmospheric Sci. 50, 69–81, https://doi.org/10.1007/s13143-014-0028-3 (2014). 36. Yu, B. & Zwiers, F. The impact of combined enso and pdo on the pna climate: a 1,000-year climate modeling study. Clim. dynamics 29, 837–851, https://doi.org/10.1007/s00382-007-0267-4 (2007). 37. Zhang, W. et al. Improved simulation of tropical cyclone responses to enso in the western north pacific in the high-resolution gfdl hiflor coupled climate model. J. Clim. 29, 1391–1415, https://doi.org/10.1175/JCLI-D-15-0475.1 (2016). 38. Patricola, C. M., Camargo, S. J., Klotzbach, P. J., Saravanan, R. & Chang, P. The influence of enso flavors on western north pacific tropical cyclone activity. J. Clim. 31, 5395–5416, https://doi.org/10.1175/JCLI-D-17-0678.1 (2018). 39. Krishnamurthy, L. et al. Impact of strong enso on regional tropical cyclone activity in a high-resolution climate model in the north pacific and north atlantic oceans. J. Clim. 29, 2375–2394, https://doi.org/10.1175/JCLI-D-15-0468.1 (2016). 40. Hodges, K., Cobb, A. & Vidale, P. L. How well are tropical cyclones represented in reanalysis datasets? J. Clim. 30, 5243–5264, https:// doi.org/10.1175/JCLI-D-16-0557.1 (2017). 41. Zängl, G., Reinert, D., Rípodas, P. & Baldauf, M. The icon (icosahedral non-hydrostatic) modelling framework of dwd and mpi-m: Description of the non-hydrostatic dynamical core. Q. J. Royal Meteorol. Soc. 141, 563–579, https://doi.org/10.1002/qj.2378 (2015). 42. Berry, G. & Reeder, M. J. Objective identification of the intertropical convergence zone: Climatology and trends from the era- interim. J. Clim. 27, 1894–1909, https://doi.org/10.1175/JCLI-D-13-00339.1 (2014). 43. Murakami, H. & Sugi, M. Effect of model resolution on tropical cyclone climate projections. Sola 6, 73–76, https://doi.org/10.2151/ sola.2010-019 (2010). 44. Bjerknes, J. Atmospheric teleconnections from the equatorial pacific. Mon. weather review 97, 163–172, https://doi. org/10.1175/1520-0493 (1969). 45. Chu, P.-S. Enso and tropical cyclone activity. Hurricanes and typhoons: Past, present, and potential 297–332 (2004). 46. Deser, C. & Wallace, J. M. Large-scale atmospheric circulation features of warm and cold episodes in the tropical pacific. J. Clim. 3, 1254–1281, https://doi.org/10.1175/1520-0442 (1990). 47. Camargo, S. J., Robertson, A. W., Barnston, A. G. & Ghil, M. Clustering of eastern north pacific tropical cyclone tracks: Enso and mjo effects. Geochem. Geophys. Geosystems 9, https://doi.org/10.1029/2007gc001861 (2008). 48. Huang, F. & Xu, S. Super typhoon activity over the western north pacific and its relationship with enso. J. Ocean. Univ. China 9, 123–128, https://doi.org/10.1007/s11802-010-0123-8 (2010). 49. Li, R. C. & Zhou, W. Changes in western pacific tropical cyclones associated with the el niño-southern oscillation cycle. J. Clim. 25, 5864–5878, https://doi.org/10.1175/JCLI-D-11-00430.1 (2012). 50. Li, R. C. & Zhou, W. Interdecadal change in south china sea tropical cyclone frequency in association with zonal sea surface temperature gradient. J. Clim. 27, 5468–5480, https://doi.org/10.1175/JCLI-D-13-00744.1 (2014). 51. Chan, J. C. & Zhou, W. Pdo, enso and the early summer monsoon rainfall over south china. Geophys. Res. Lett. 32, https://doi. org/10.1029/2004GL022015 (2005). 52. Wang, L., Huang, G., Zhou, W. & Chen, W. Historical change and future scenarios of sea level rise in macau and adjacent waters. Adv. Atmospheric Sci. 33, 462–475, https://doi.org/10.1007/s00376-015-5047-1 (2016). 53. Wang, W. & Zhou, W. Statistical modeling and trend detection of extreme sea level records in the pearl river estuary. Adv. Atmospheric Sci. 34, 383–396, https://doi.org/10.1007/s00376-016-6041-y (2017). 54. Tolman, H. L. et al. User manual and system documentation of wavewatch iii tm version 3.14. Tech. note, MMAB Contribution 276, 220 (2009). 55. Hurrell, J. W., Hack, J. J., Shea, D., Caron, J. M. & Rosinski, J. A new sea surface temperature and sea ice boundary dataset for the community atmosphere model. J. Clim. 21, 5145–5153, https://doi.org/10.1175/2008JCLI2292.1 (2008). 56. Saha, S. et al. The ncep climate forecast system reanalysis. Bull. Am. Meteorol. Soc. 91, 1015–1058, https://doi. org/10.1175/2010BAMS3001.1 (2010). 57. Mantua, N. J., Hare, S. R., Zhang, Y., Wallace, J. M. & Francis, R. C. A pacific interdecadal climate oscillation with impacts on salmon production. Bull. american Meteorol. Soc. 78, 1069–1080, https://doi.org/10.1175/1520-0477 (1997). 58. Wolter, K. & Timlin, M. S. El niño/southern oscillation behaviour since 1871 as diagnosed in an extended multivariate enso index (mei. ext). Int. J. Climatol. 31, 1074–1087, https://doi.org/10.1002/joc.2336 (2011). 59. Mote, P. W. et al. Preparing for climatic change: the water, salmon, and forests of the pacific northwest. Clim. change 61, 45–88, https://doi.org/10.1023/A:1026302914358 (2003). 60. Wharton, S. & Falk, M. Climate indices strongly influence old-growth forest carbon exchange. Environ. Res. Lett. 11, 044016, https:// doi.org/10.1088/1748-9326/11/4/044016 (2016). Scientific Repo R t S | (2020) 10:5187 | https://doi.org/10.1038/s41598-020-62138-0 9 www.nature.com/scientificreports www.nature.com/scientificreports/ Acknowledgements This work was funded by the Korea Meteorological Administration Research and Development Program under Grant KMI2018 − 07910. Also, this research was supported by the "Technology development for Practical Applications of Multi-Satellite data to maritime issues" funded by the Ministry of Ocean and Fisheries, Korea. Author contributions Conceptualization, S.Y. and J.-H.O.; methodology, S.Y.; software, S.Y.; validation, S.Y.; writing-original dra ft preparation, S.Y.; writing-review and editing, S.Y. and J.-H.O.; visualization, S.Y.; supervision, J.-H.O.; project administration, J.-H.O.; funding acquisition, J.-H.O. Competing interests e a Th uthor(s) declare no competing interests. Additional information Supplementary information is available for this paper at https://doi.org/10.1038/s41598-020-62138-0. Correspondence and requests for materials should be addressed to J.-H.O. Reprints and permissions information is available at www.nature.com/reprints. Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2020 Scientific Repo R t S | (2020) 10:5187 | https://doi.org/10.1038/s41598-020-62138-0 10
Scientific Reports – Springer Journals
Published: Mar 20, 2020
It’s your single place to instantly
discover and read the research
that matters to you.
Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.
All for just $49/month
Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly
Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.
Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals.
Read from thousands of the leading scholarly journals from SpringerNature, Wiley-Blackwell, Oxford University Press and more.
All the latest content is available, no embargo periods.
“Hi guys, I cannot tell you how much I love this resource. Incredible. I really believe you've hit the nail on the head with this site in regards to solving the research-purchase issue.”Daniel C.
“Whoa! It’s like Spotify but for academic articles.”@Phil_Robichaud
“I must say, @deepdyve is a fabulous solution to the independent researcher's problem of #access to #information.”@deepthiw
“My last article couldn't be possible without the platform @deepdyve that makes journal papers cheaper.”@JoseServera