Preliminary evaluation of MJO simulation in GAMIL3(Grid-point atmospheric model of IAP LASG)

PU Ye,LI Lijuan,LIU Hongbo,DONG Li and WANG Bin

State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics,Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing,China

ABSTRACT This study evaluates the performance of the Grid-point Atmospheric Model of IAP LASG,version 3 (GAMIL3), in simulating the Madden-Julian Oscillation (MJO), based on the CMIP6(phase 6 of the Coupled Model Intercomparison Project) AMIP (Atmospheric Model Intercomparison Project) simulation. Results show that GAMIL3 reasonably captures the main features of the MJO, such as the eastward-propagating signal in the MJO frequency band, the symmetric and asymmetric structures of the MJO, several convectively coupled equatorial waves, and the MJO life cycle. However, GAMIL3 underestimates the MJO amplitude, especially for outgoing longwave radiation, as do most CMIP5 models, and simulates slow eastward propagation.

KEYWORDS MJO;GAMIL3;evaluation

1. Introduction

In the early 1970s,Madden and Julian(1971,1972)found a quasi-periodicoscillationintheatmosphereneartheequator on a time scale of 40-50 days. Subsequent studies have shown 30-60-day quasi-periodic oscillations throughout the tropical atmosphere,and even the global atmosphere,and this is regarded as one of the most important atmospheric circulation systems(Krishinamurti and Subrahmann 1982; Li 1983; Knuston and Weickmann 1987; Liebmann,Hendon, and Glick 1994). Such intraseasonal atmospheric oscillations near the equator are referred to as the Madden-Julian Oscillation(MJO).

Observations and numerical simulations have indicated that the MJO has great influence on tropical weather and climate phenomena, including the Indian Ocean summer monsoon, the Australian winter monsoon, the South China Sea summer monsoon,tropical cyclone activity, and El-Ni?o Southern Oscillation (Yasunari 1979, 1980; Hall, Matthews, and Karoly 2001; Li and Long 2001; Annamalai and Sperber 2005; Lorenz and Hartmann 2006; Wheeler et al. 2009;Wang et al. 2017a). The MJO also influences the weather and climate outside the tropics through atmospheric teleconnection mechanisms (Weickmann,Lussky, and Kutzbach 1985; Ferranti et al. 1990; Lin,Brunet, and Derome 2009; Henderson, Maloney, and Barnes 2016) and modulates the global climate (Wang et al. 2013; Wang, Li, and Chen 2019; Chen, Yu, and Zheng 2016; Chen et al. 2017).

However, current atmospheric models have diffi-culty in providing realistic MJO simulations (Slingo et al. 1996). An evaluation of MJO simulation performance of models participating in phase 5 of the Coupled Model Intercomparison Project (CMIP5; Hung et al. 2013; Ahn et al. 2017) found that most models underestimate the MJO amplitude, and that the simulated MJO propagation is fast compared with observations. Moreover, the symmetry and asymmetry of the convectively coupled Kelvin waves (CCKWs) between the Northern and Southern Hemisphere using the CMIP5 coupled models were evaluated, suggesting that some models overestimate the amplitude of CCKWs in the Southern Hemisphere (Wang and Li,2017). In operational MJO prediction, assimilating the observed surface signals is considered as the minimum requirement (Cui et al. 2020).

The Grid-point Atmospheric Model of IAP LASG(GAMIL) was developed at the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics (LASG), Institute of Atmospheric Physics (IAP), Chinese Academy of Sciences, Beijing, China. The simulation performance for the MJO and other intraseasonal oscillations by early versions of the model (GAMIL1 and GAMIL2)has been evaluated previously (Yang et al. 2009; Mao and Li 2012; Xie et al. 2012). The newest version(GAMIL3) has recently been developed and has completed the Atmospheric Model Intercomparison Project(AMIP) runs of phase 6 of the Coupled Model Intercomparison Project (CMIP6). This paper evaluates the ability of GAMIL3 to simulate the MJO.Descriptions of the model and observational data are presented in section 2. Section 3 evaluates the performance of MJO simulations. Finally, a summary and discussion are presented in section 4.

2. Model and observational data

GAMIL3 uses a hybrid 2D decomposition to improve parallel scalability (Liu et al. 2014). The horizontal resolution is ～2°, with a 180 × 80 longitude-latitude grid, and the vertical resolution is 26 layers from the surface to the top of the model (2.194 hPa) with an σ-coordinate system. To improve water vapor conservation, the weight of the two-step shape-preserving advection scheme of the previous version (GAMIL2; Li et al. 2013) has been modified in GAMIL3. A planetary boundary layer scheme based on turbulence kinetic energy and a stratocumulus cloud-fraction scheme based on estimated inversion strength are adopted in GAMIL3 (Guo and Zhou 2014; Sun, Li, and Wang 2016). In addition, the convective momentum transport and the simple parameterization to describe anthropogenic aerosol effects recommended by CMIP6 have been added (Stevens et al. 2017; Shi,Zhang, and Liu 2019). Following CMIP6 recommendations, the external forcings, such as greenhouse gases, ozone, land use, the solar constant, and aerosols, have also been updated (https://esgf-node.llnl.gov/search/input4mips/).

In this paper, AMIP runs for the years 1980-2005 are analyzed following the CLIVAR (U.S. Climate Variability and Predictability) MJO Working Group’s diagnostics package.We mainly use daily data and apply a 20-100-day Lanczos bandpass filter (Waliser et al. 2009). The analysis will focus on boreal winter.

To verify the model results, the following observation/reanalysis datasets are used: satelliteobserved daily outgoing longwave radiation (OLR)from the Advanced Very High Resolution Radiometer (AVHRR; Liebmann and Smith 1996);and upper-tropospheric (200 hPa) and lower-tropospheric (850 hPa) winds from NCEP-NCAR reanalysis data (Kalnay et al. 1996).

3. Preliminary evaluation of MJO simulation in GAMIL3

The 20-100-day-filtered 850-hPa zonal wind (U850)during boreal winter (November-April) and its percentage contribution to the overall variance are shown in Figure 1. Large percentage values are seen over the tropical Indian Ocean (35%-40%), the Maritime Continent (> 40%), and the tropical and northern Pacific Ocean (35%-40%) in the NCEPNCAR reanalysis (Figure 1(a)). GAMIL3 captures the above spatial distribution of large percentage values well, especially over the central Maritime Continent(Figure 1(c)). The spatial distribution of OLR variance(Figure 1(b,d)) is similar to that ofU850. Two eastwest band-shaped large-variance-percentage belts are located over the northern and southern Indian Ocean in GAMIL3, whereas one region of large variance is located over the tropical Indian Ocean in observations (Figure 1(b,d)). GAMIL3 also has a negative variance percentage bias over the Maritime Continent.

Figure 2 shows the wavenumber-frequency power spectra of 10°S-10°N-averaged OLR andU850, which provides a convenient metric of MJO characteristics.The observed spectral power ofU850 and OLR are concentrated over the domain of eastward wavenumbers 1-3, with periods of 30-80 days (Figure 2(a,c)).The model contains significant eastward-propagating signals in the MJO band for both OLR andU850.However, the simulated energy is smaller than observations of OLR and the dominant peak is outside of the MJO band in U850 (Figure 2(b,d), similar to other atmospheric general circulation model simulations (e.g., ECHAM6; Crueger, Stevens, and Brokopf 2013). The eastward/westward (E/W) ratio, which is obtained by dividing the sum of spectral power over the MJO band by that of its westward propagating counterpart, and the eastward/observation (E/O) ratio, which is formulated by normalizing the sum of spectral power within the MJO band by the observed value,are calculated from frequency-wavenumber spectra,and are metrics used to indicate the robustness of eastward-propagating features of the MJO (Table 1).The observed (simulated) E/W ratios are 3.7 (1.7) and 3.8 (2.9) for OLR andU850, respectively. The E/O ratio is 0.4 for OLR and 0.9 forU850 in GAMIL3. These results indicate that the model generally underestimates the eastward propagation of convection, but the simulation of the eastward propagation of circulation is much better. These results are similar to those from most CMIP5 coupled models(Subramanian et al. 2011; Ahn et al. 2017).

Figure 1.20-100-day band-pass-filtered U850 variance(contours;units:m2 s?2)and its percentage contribution(shaded)to the total variance in(a)NCEP-NCAR reanalysis and(c)GAMIL3.(b,d)As in(a,c)but for OLR(units:W2 m?4)from NOAA AVHRR and GAMIL3,respectively.

Figure 2. November-April (1980-2005) wavenumber-frequency spectra of 10°N-10°S-averaged daily U850 from (a) NCEP and (b)GAMIL3;and daily OLR fields from(c)NOAA satellite OLR and(d)GAMIL3.The dashed lines indicate the MJO band(30 d to 80 d).

Table 1.Eastward/westward(E/W)and eastward/observation(E/O)ratios within the MJO frequency-wavenumber ranges(30-80-day period, wavenumbers 1-3) for the single equatorial fields of OLR and U850.

We can clearly distinguish the main MJO modes and convectively coupled equatorial waves (CCEWs),including Kelvin,equatorial Rossby,mixed Rossby-gravity,and eastward inertia-gravity waves, through the coherencesquared and phase relationship between equatorial OLR andU850 for the symmetric and asymmetric components (Figure 3). In the MJO band, observations show high coherence and an approximate 90° phase lag between convection and 850-hPa winds (Figure 3(a)).The GAMIL3 simulation is almost consistent with observational data and effectively reproduces the MJO and several CCEWs-both the symmetric and asymmetric

Figure 4. Composite MJO OLR (shaded; units W m?2) and 850-hPa wind anomalies (vectors; units: m s?1) in eight phases from November to April from(a)observations and(b)GAMIL3.The reference vector in units of m s?1 is in the upper right,and the phase and number of days used are in the lower right,of each panel.

4. Summary and discussion

This study evaluates MJO simulations of GAMIL3 using CMIP6 AMIP-type runs. Results show that GAMIL3 can simulate the main features of the MJO. For example,the model generally simulates the variance ofU850 and OLR over the tropics. However, the eastward-propagating signal in both OLR andU850 from the power during the MJO band is weaker than in observations.In general, the simulated circulation is better than that of OLR. GAMIL3 is also able to effectively reproduce both symmetric and asymmetric waves. The MJO and several CCEWs are closely consistent with observational data, exhibiting strong coherence in the MJO band. Phase lags between equatorial OLR andU850 are similar to those of observations. The model also simulates the MJO life cycle well. However, it underestimates the MJO amplitude, especially as it relates to OLR, as do most CMIP5 models (Hung et al. 2013;Jiang et al. 2015). This is related to poor climatological mean simulations of Indian Ocean OLR and precipitation in GAMIL3.

Although GAMIL3 can simulate the basic characteristics of the MJO, there is still room for improvement. For example, convective parameterization is considered as the critical process for MJO simulation(Crueger, Stevens, and Brokopf 2013; Wang et al.2017b).However,before concentrating on this aspect,it is necessary to evaluate the contributions of each model improvement(including convective parameterization) in GAMIL3 to MJO simulations relative to the previous model version, applying the new dynamicsoriented diagnostics for the MJO proposed by Wang et al. (2018). In addition, many previous studies have shown air-sea interactions play an important role in the MJO(Crueger,Stevens,and Brokopf 2013;Jiang et al.2015).Further analysis is needed using results of the coupled FGOALS-g3 model.

Disclosure statement

No potential conflict of interest was reported by the authors.

Funding

This study was jointly supported by the National Key Research and Development Program of China [grant number 2017YFA0603903] and the National Natural Science Foundation of China[grant numbers 41622503 and 41775101].