Process-Based Evaluation of Intraseasonal Oceanic Kelvin Waves in the Pacific Ocean in CMIP6 Models

Abstract Oceanic intraseasonal Kelvin waves (KWs) help modulate upper-ocean thermal characteristics, providing feedbacks to important coupled air–sea phenomena in the tropics. The recent availability of daily thermocline depth fields from several phase 6 of the Coupled Model Intercomparison Project (CMIP6) models makes it possible to evaluate the performance of KWs and identify potential sources of bias. Most models fail to simulate a realistic spatial distribution of KW variability. Models simulate a large variability of KWs in the western or eastern Pacific rather than in the central Pacific as observed. The modeled KWs propagate slowly (about 1.5 m s −1 ) compared to observations (about 2.5 m s −1 ). This slow propagation is also identified in wavenumber–frequency spectra for KWs and meridional KW structures, which is more consistent with a second baroclinic mode structure in models compared to the first baroclinic mode structure in observations. An analysis of the relative contributions of the vertical wavenumber and background ocean stability to KW phase speeds indicates that the high vertical wavenumber bias in models contributes most to the slow propagation, in which the higher-than-observed vertical wavenumbers imply the biased incorporation of higher baroclinic modes in the model KW structure. This finding is further supported by the results of vertical mode decomposition that incorporates background density profiles. These results indicate that a realistic representation of the KW vertical structure is essential to produce realistic KW propagations in models. Significance Statement Oceanic intraseasonal Kelvin waves (KWs) play a significant role in regulating the heat content and temperature of the ocean, which provides feedbacks to coupled ocean–atmosphere phenomena in the tropics such as El Niño–Southern Oscillation (ENSO). Consequently, identifying the biases in KWs and the potential sources of those biases in state-of-the-art models is essential to improve simulations of ENSO and its diversity and advance forecasts of weather extremes around the globe induced by ENSO. We find that KWs in the models propagate slower than observations, mainly due to biases in their vertical structure, with secondary effects due to biases in model ocean stability. These issues may be potential sources for current imperfect ENSO model simulations and predictions.