The Asian Subtropical Westerly Jet Stream in CRA-40, ERA5, and CFSR Reanalysis Data: Comparative Assessment

The Asian subtropical westerly jet (AWJ) exerts crucial influences on Eurasian continent weather and climate. This paper analyzes the advantages and limitations of CRA-40, which is China’s first generation 40-yr (1979–2018) global atmosphere and land reanalysis product, in describing the characteristics of AWJ, compared with the ECMWF Reanalysis version 5 (ERA5) and NCEP Climate Forecast System Reanalysis (CFSR). The results show a close agreement across the three reanalyses on the whole. (1) In terms of climatology, overall differences of 200-hPa zonal wind across the three reanalyses are within ± 0.5 m s−1 (i.e., ± 2%). Large differences with maxima of ± 2 m s−1 (± 5%) appear over the Iranian Plateau and south of the Tibetan Plateau in the mid-upper troposphere in winter. (2) For seasonal cycle, the position and intensity of the AWJ centers in the three reanalyses are highly consistent, with correlation coefficient over 0.98. But there are some discrepancies in the zonal shift of the western AWJ center during the transition season. (3) On the interannual timescale, intensity of all AWJ centers varies consistently among the three reanalyses, while larger differences appear in their meridional displacement, especially in the eastern AWJ center. (4) For long-term variations, the three reanalyses all present a significant northward movement of the westerly jet axis in winter, and a southward displacement over central Asia (40°–80°E) and a northward migration over East Asia (80°–110°E) in summer. Thus, this study has provided confidence that CRA-40 has comparable performance with ERA5 and CFSR in depicting the characteristics of AWJ.


Introduction
The Asian subtropical westerly jet stream (AWJ), a narrow and strong westerly belt, exists in the upper troposphere and lower stratosphere over the Eurasian continent all year around. As one of the planetary scale systems at middle latitudes, AWJ is often divided into the eastern part (EAWJ) and the western part (WAWJ) according to their jet center locations and affected area. For instance, AWJ in summer has three active centers located in the Northwest Pacific Ocean, north of the Tibetan Plateau, and the Caspian, respectively, along 40°N (Yang and Zhang, 2008a;Du et al., 2016). EAWJ is closely related to Meiyu  and the summer monsoon over East Asia (Yeh et al., 1958), while WAWJ is intimately linked with summer precipitation in Northwest China (Yang and Zhang, 2008b;Zhao et al., 2014a) and central Asia (Schiemann et al., 2009;Zhao et al., 2018). Thus, it is of great importance for the AWJ study given its impacts on regional weather and climate.
Meteorological data with high temporal and spatial resolutions have always been in high demand. Unfortu-nately, earlier studies on AWJ were based on inhomogeneous and raw observations at sounding and pilot balloon (pibal) stations. In the 1950s, owing to the rapid development of the aerological network in China, Chinese researchers were able to carry out a series of synoptic and aerological investigations and made remarkable achievements on the general circulation over East Asia (Staff Members of the Section of Synoptic and Dynamic Meteorology, Institute of Geophysics and Meteorology, Academia Sinica, 1957Sinica, , 1958aZhu et al., 1990). It is proposed that the seasonal evolution of EAWJ, especially its south-north migration, can be a signal for the seasonal transition of the general circulation, the onset and retreat of the summer monsoon, and occurrence and termination of Meiyu over East Asia Yeh et al., 1958;Lu and Schneider, 2017). With the aerological observations increasing, Dao et al. (1965) further illustrated the major features of general circulation over East Asia in the upper troposphere and stratosphere.
With rapid development of numerical weather prediction models, reanalysis products have been prevailing for studying the atmospheric sciences. In the early 21st century, abundant work has been carried out to reveal more refined features of AWJ based on monthly mean variables of modern reanalysis data. For EAWJ, the three-dimensional structures are depicted comprehensively, the seasonal evolution in meridional migration and zonal displacement is reexamined in terms of the intensity and location of its centers and axis, and the thermal mechanism is further investigated by analyzing the meridional difference of air temperature (Kuang et al., 2007). The climate features of AWJ in summer are examined from the perspective of Rossby wave activity (Yang and Zhang, 2008b). In addition, reanalyses data on higher temporal-spatial resolution complement those results based on the monthly mean data. It is noted that the EAWJ core moves westward rapidly during 35th-39th pentads (Zhang et al., 2006) and the median dates of the westerly jet in the Tibetan Plateau point to 28 April and 12 October, by identifying the locations of the jet axis at 6-h intervals (Schiemann et al., 2009).
On interannual and interdecadal timescales, the statistical relationship between AWJ and regional precipitation and associated physical mechanisms have been investigated. In summer, the first two dominant modes of AWJ represent its meridional displacement and southwestnortheast titling (Du et al., 2016). Previous studies demonstrated that the location changes of EAWJ in both meridional and zonal directions have pronounced impacts on the spatial distribution of rainfall anomalies over East Asia through influencing the South Asian high and the western North Pacific subtropical high (Liao et al., 2004;Lin and Lu, 2005;Du et al., 2009;Xie et al., 2015). Moreover, the axis of EAWJ in summer has shifted southward while its intensity in winter enhanced during 1980-2004relative to 1958-1979(Zhang and Huang, 2011. In winter, the affected area goes southward along with the retraction of EAWJ (Li and Sun, 2015;Hunt et al., 2018).
It is also proved that the WAWJ has a significant influence on summer rainfall in Northwest China and central Asia (Yang and Zhang, 2008b;Schiemann et al., 2009;Wei et al., 2017;Ji et al., 2020), and the role of its position is more important than its strength (Zhao et al., 2014a(Zhao et al., , 2018. From the perspective of circulation, the southward displacement of WAWJ is responsible for the rainfall increase in Northwest China, which results in anomalous southerly, water vapor transport, and ascending motion in the lower troposphere (Zhao Y. et al., 2014a(Zhao Y. et al., , b, 2019Huang et al., 2015;Peng and Zhou, 2017). Quite a few reasons are attributed to the interannual variation of AWJ position, which is linked with Arctic Oscillation, sea surface temperature anomalies, the Indian summer monsoon, and the South Asian high (Zhao Y. et al., 2014a(Zhao Y. et al., , 2019Du et al., 2016;Wei et al., 2017). In the recent decade, reanalysis data are more widely used to evaluate the ability of numerical models in reproducing the fundamental features of AWJ and its relationship with regional precipitation (Huang and Liu, 2011;Yan et al., 2014;Du et al., 2017;Zhao et al., 2018;Lin et al., 2019;Yan et al., 2019;Fu et al., 2020).
The structure of the paper is as follows. Section 2 gives a brief description of the data and method used in this study. Subsequently, comparison results are presented in Section 3, focusing on the climatic structure (Section 3.1), seasonal cycle (Section 3.2), interannual variability (Section 3.3), and long-term change (Section 3.4) of AWJ, respectively. To the end, we conclude with a summary and discussion of our results in Section 4. Yu, X. J., L. X. Zhang,T. J. Zhou,et al. 47 2. Data and methods

Data
In this study, we employed three global reanalysis products, CRA-40 from CMA, ERA5 from ECMWF, and CFSR from NCEP. The details of the data are provided as follows.
CRA-40 is the first generation global atmosphere and land reanalysis data produced by CMA. The CMA reanalysis (CRA) project was launched in 2014 (Liu et al., 2017;Zhao B. et al., 2019). The current CRA system is based on the Global System Model (GSM) of the Global Forecast System (GFS) of NCEP and its Gridpoint Statistical Interpolation (GSI) 3DVAR data assimilation system. Compared with other reanalyses, multiple data from conventional observations and satellite instruments, especially over East Asia, have been assimilated into the CRA (Liao et al., 2018;Wang et al., 2018;Yin et al., 2018). CRA-40 has a spacial resolution of 0.312° latitude × 0.312° longitude, with 64 levels in the vertical.
ERA5 is the fifth generation of ECMWF atmospheric reanalysis of the global climate, following FGGE [First GARP (Global Atmospheric Research Program) Global Experiment], ERA-15, ERA-40, and ERA-Interim (Gibson et al., 1997;Uppala et al., 2005;Dee et al., 2011). ERA5 has a horizontal resolution of approximately 31 km and a vertical resolution of 137 levels (Hersbach and Dee, 2016). The ERA data can be obtained at https:// www.ecmwf.int/en/forecasts/datasets/reanalysis-data sets/era5. CFSR is the third generation global reanalysis developed by NCEP, using the Coupled Forecast System (CFS) model (Saha et al., 2010(Saha et al., , 2014. CFSR uses a coupled atmosphere-ocean-sea ice-land model, assimilates conventional observations and satellite radiance, and includes the time evolution of CO 2 concentrations. The CFSR data have the same spatial resolution as CRA-40, which are available at https://climatedataguide.ucar.edu/ climate-data/climate-forecast-system-reanalysis-cfsr.
In the subsequent comparative assessment, monthly mean variables, including zonal wind, meridional wind, and geopotential height, during the period of 1979-2018 from the three reanalysis datasets are employed. For comparison consistency, all datasets are bilinearly interpolated to a horizontal resolution of 0.5° × 0.5° with 17 vertical levels from 1000 to 100 hPa at an interval of 50 hPa.

Methods
During the boreal summer (JJA) and winter (DJF) sea-sons, AWJ is located at its preferred positions, and the boreal spring (MAM) and autumn (SON) seasons are the transition periods. Thus, we mainly concentrate on comparison of results in JJA and DJF seasons.
Various verification statistics including pattern correlation coefficient, Pearson linear correlation coefficient, root mean square difference (RMSD), and relative difference, are applied to assess the performance of the three reanalysis datasets. Standard deviation reveals the variation or dispersion from the average, and the interannual standard deviation (ISD) denotes interannual variations relative to the climatic mean (Huang et al., 2014). In addition, the Empirical Orthogonal Function (EOF) analysis method is used to calculate indices and long-term trends of AWJ (Du et al., 2016;Zhao et al., 2014aZhao et al., , 2018.

Climatic three-dimensional structure of AWJ
We first assess the differences among the three reanalysis datasets in capturing the climatic horizontal distribution of AWJ ( Fig. 1) and the associated circulations In winter, the jet stream is strong and wide, quasi-zonally orientated within 20°-50°N (Fig. 1a), and the westerlies prevail over Eurasia because of the meridional temperature gradient caused by warm tropics and cold Eurasian continent ( Fig. 2a). Two centers appear over western Asia (30°-50°E), and the eastern Asia and western Pacific (110°-160°E) region, respectively. The eastern center (with a maximum wind speed above 70 m s −1 ) is much stronger and larger than the western one (with the maximum wind speed above 50 m s −1 ). The horizontal structures of AWJ in spring and autumn are generally similar to that in winter, except that the jet stream is located slightly northward and the jet centers are weaker (figures omitted). In summer, the horizontal structure of the 200-hPa zonal wind differs remarkably from that in other seasons. The jet stream migrates northward and becomes narrow, occupying over 30°-45°N, and the intensity of the zonal wind is strikingly weakened to 25 m s −1 (Fig.  1b). This is resulted from the summer South Asian high, which is an intense anticyclonic system located over the Iranian and Tibetan plateaus (Fig. 2b). There exhibit three active jet cores located over the Caspian (50°E), north of the Tibetan Plateau (80°-100°E), and the Northwest Pacific Ocean (140°-160°E), respectively, among which the central one is the strongest with a maximum speed above 30 m s −1 .
The climatic horizontal structures of AWJ from ERA5 and CFSR show comparable spatial patterns and intensity to CRA-40 (figures omitted), with their pattern correlation coefficients excelling 0.99 and regional average differences below 0.4 m s −1 . To reveal the detailed differences of the three reanalyses, the differences of the zonal wind and circulations at 200 hPa in winter and summer between each pair are then calculated as shown in Figs. 1c-h; 2c-h. In general, the differences are slightly larger in winter than in summer, although their averages are quite small. In winter, CRA-40 and CFSR exhibit high consistency, with their zonal wind (geopotential height) differences in the range of ± 0.5 m s −1 (± 8 gpm) over almost the whole region (Figs. 1e, 2e), whereas ERA5 shows larger inconsistency compared with them (Figs. 1c, g; 2c, g). One region with zonal wind difference exceeding 2 m s −1 is located over the Iranian Plateau and the south of the Tibetan Plateau (50°-90°E) between the two jet cores, accompanied by positive geopotential height differences and anticyclonic circulation anomalies.  The other region with obvious differences in zonal wind and circulations is situated in western Pacific, southeast of the stronger jet core. In summer (Figs. 1c, e, g; 2c, e, g), differences in zonal wind (geopotential height) between the three reanalyses range within ± 0.5 m s −1 (± 8 gpm) over a majority of regions, along with weaker circulation anomalies.
We then evaluate the performance of the three reanalyses in capturing the vertical structures of climate mean AWJ, as illustrated in Fig. 3. Figures 3a and b show the vertical structures of the zonal wind averaged over 30°-160°E in winter and summer from CRA-40. In winter, the whole atmosphere is almost dominated by the west-erly. The jet center is at 200 hPa around 30°N with the maximum wind speed of about 50 m s −1 . In summer, the jet center migrates northward to around 42°N and the intensity weakens to approximately 20 m s −1 . Meanwhile, the tropical region is occupied by a tilted structure, with relatively strong easterly in the middle and upper troposphere and weak westerly in the lower troposphere.
The vertical structures of AWJ from ERA5 and CFSR resemble CRA-40 with their pattern correlation coefficients exceeding 0.99 (figures omitted). Figures 3c-h further give the vertical structure differences in zonal wind across the three reanalyses. Overall, the results agree with that of horizontal distribution of AWJ (Fig. 1). The  differences in summer change within the range of ± 0.5 m s −1 (Figs. 3c, e, g), weaker than that in winter (Figs. 3d,f,h). In winter, the difference with a maximum > 0.5 m s −1 occupies the south of 40°N, consistent with its horizontal distribution at 200 hPa. It is obvious that the most significant difference between ERA5 and CRA-40 appears in the middle and upper troposphere (above 500 hPa), with the maximum wind speed over 1.5 m s −1 (Fig.  3c). CFSR and CRA-40 differ mainly in the lower troposphere (400-850 hPa), with the maximum wind of approximately 1 m s −1 (Fig. 3e). The inconsistency between ERA5 and CFSR appears at middle and upper levels (above 400 hPa), with the difference values below 1 m s −1 (Fig. 3g).
From the above analysis, we can see that there exhibit some discrepancies among the three reanalyses, especially over the Iranian Plateau and the south of the Tibetan Plateau (around 30°N) in the middle and upper troposphere. Considering that their overall relative differences range within 2% (figures omitted), the three reanalyses also have a high consistency in depicting the threedimensional structures of AWJ.

Seasonal cycle of AWJ
To further compare the performance of CRA-40, ERA5, and CFSR in depicting the seasonal cycle of AWJ, we select latitude, longitude, and intensity (measured by zonal wind speed) of the jet centers as the metrics. Given there are two main AWJ centers in each month (figures omitted), we divide AWJ into two parts: the western section (30°-70°E) and the eastern section (70°-160°E). Figure 4 illustrates the seasonal cycle of AWJ centers in terms of longitude, latitude, and intensity derived from CRA-40, ERA5, and CFSR. AWJ centers identified from the three reanalyses all undergo a pronounced seasonal cycle. The south-north displace- ments of the jet centers are within 10 latitudes, whereas the east-west shifts are more than 50 longitudes.
As for the western center (Figs. 4a,c,e), it stays at about (27°N, 35°E) from January to March, with the maximum zonal wind speed of about 55 m s −1 . During its northward migration from April to July, the center reaches the northernmost position at 44°N in July, moves eastward to near 50°E, and weakens abruptly to about 32 m s −1 . Subsequently (August-November), the center withdraws southward and intensifies gradually, while continues to move eastward and reaches the easternmost position at near 69°E. In December, the center recedes to its southernmost and westernmost position at about (27°N, 35°E), with the zonal wind speed of 48 m s −1 .
The seasonal cycle of the eastern center (Figs. 4b, d, f) exhibits slightly different characteristics, compared with the western one. From January to March, the center is located at about (32°N, 140°E), 5 degrees north compared to the western center, with its intensity decreasing from the maximum of 77 m s −1 in January. During the northward migration from April to August, the eastern center settles in its northernmost location of near 45°N in August, one month later than the western one, while the minimum zonal wind speed of 32 m s −1 appears in July. Meanwhile, the center migrates zonally by a large margin, reaching its westernmost position of 90°E nearby in August. After that (from September to December), the center recedes quickly to its wintertime position (32°N, 140°E), with the zonal wind intensified.
On the whole, the three reanalyses are highly consistent, with most correlation coefficients higher than 0.98 (Table 1), in describing the seasonal cycle characteristics of AWJ. It is worth noting that there are still some recognizable discrepancies between them, especially during the transition period. For the western center, in spring, the AWJ position derived from ERA5 is further north and east than CRA-40 and CFSR, while its intensity is weaker than CRA-40 and CFSR. In autumn, the AWJ center in CRA-40 recedes southward slowly while moves eastward earlier than in the other two reanalyses. The correlation coefficient between ERA5 and CRA-40 (CFSR) is only 0.76 (0.80). As for the eastern center, the differences are much smaller. In spring and winter, the AWJ center in ERA5 sways at a slightly larger range in the latitude and longitude than CRA-40 and CFSR. Therefore, we need to pay attention to the disagreement in the three reanalyses when examining the zonal shift of AWJ, especially the western center.

Interannual variability of AWJ
We further compare the ability of the three reanalyses in depicting the interannual variation of AWJ, which is represented by ISD (figures omitted). In general, ISD is higher, with a wider coverage, in winter than in summer, similar to the distribution of AWJ. Although there are large differences in the amplitude and pattern of ISD between winter and summer, ISD centers reside at the entrance and exit regions along the AWJ axis, which accords with the previous study by Du et al. (2016). In terms of the differences of ISD among the three reanalyses, relatively high values are located along the jet axis, and the differences in winter are still slightly higher than in summer.
As shown in Figs. 1a, b, there are two AWJ centers in winter, denoted as W1 (western) and W2 (eastern), respectively, and three in summer, denoted as S1 (western), S2 (central), and S3 (eastern), respectively. We then examine the interannual variability of the five centers across the three reanalyses, respectively, with respect to latitude, longitude, and intensity (Fig. 5). Every center exhibits obvious interannual variability during the past 40 years, and the consistency across the three reanalyses is generally higher in summer (S1, S2, and S3) than in winter (W1 and W2). For the meridional migration (Figs. 5a-e), the three reanalyses present considerable agreement with each other, with correlation coefficient ≥ 0.88 and RMSD ≤ 0.35 latitudes. As to the zonal displacement (Figs. 5f-j), the differences for four centers (W1, W2, S1, and S2) across the three reanalyses tend to be larger, especially in the recent decade. The correlation coefficients are relatively lower and most RMSDs are approximately 1.0 longitude, some even exceeding 2.0 longitudes. Taking S2 as an example (Fig. 5i), which is located to the north of the Tibetan Plateau, the correlation coefficient between CRA-40 and ERA5 (between ERA5 and CFSR) is only 0.73 (0.79) while the RMSD is 2.02 Yu, X. J., L. X. Zhang, T. J. Zhou, et al.
(1.93) longitudes. Of the three metrics, the intensity for each center from the three reanalyses shows almost complete agreement, with the correlation coefficient approximating 1.0 and RMSD ≤ 1.0 m s −1 (Figs. 5k-o). We further quantitatively assess the performance of the three reanalyses in reproducing the interannual variability of AWJ centers in each month. Here, we still di-vide AWJ into the western and eastern sections in calculation, as in Section 3.2. The statistical results including Pearson correlation coefficient, RMSD, and average difference are summarized in Fig. 6. On the whole, the three reanalyses agree more to each other for the western center than the eastern one, with relatively higher correlation coefficients, lower RMSDs, and lower average dif- Fig. 5. Time series of AWJ centers in terms of (a-e; °) latitude, (f-j; °) longitude, and (k-o; m s −1 ) intensity in winter (W1 and W2) and summer (S1, S2, and S3) during 1979-2018, identified from CRA-40 (green line), ERA5 (red line), and CFSR (blue line). Numbers after C (R) denote the Pearson correlation coefficient (RMSD) between CRA-40 and ERA5, CRA-40 and CFSR, and ERA5 and CFSR in sequence. W1 and W2 indicate the jet centers over the western Asia (30°-50°E) and eastern Asia and western Pacific (110°-160°E) in winter; while S1, S2, and S3 denote the jet centers over the Caspian (50°E), north of the Tibet Plateau (80°-100°E), and the Northwest Pacific Ocean (140°-160°E) in summer, respectively, as shown in Figs. 1a, b.

54
Journal of Meteorological Research Volume 35 ferences. Of the three metrics, the three reanalyses have the highest consistency in the intensity while the largest difference in the zonal displacement of AWJ centers. For intensity in each month, the correlation coefficients all exceed 0.94, the RMSDs are below 2.0 m s −1 , and the average differences are almost within ± 0.5 m s −1 . As to the latitude, the correlation coefficients are within 0.8-0.94, while most RMSDs are below 2.0 latitudes. As regard to the longitude, the correlation coefficients are more scattered (0.7-0.98), and the RMSDs (average differences) exceed 3 (± 0.5) longitudes, especially for the eastern center.
Previous studies have demonstrated that the first two leading EOF modes (EOF1 and EOF2) of 200-hPa zonal wind in summer reflect the changes in position and strength of the subtropical westerly jet over a certain region, thus the corresponding principal components, PC1 and PC2, have often been defined as the position and strength indices, respectively, to investigate the impacts of the subtropical westerly jet to weather and climate (Zhao et al., 2014a(Zhao et al., , 2018Du et al., 2016;Yan et al., 2019). Similarly, by applying the EOF decomposition to    in winter and summer, respectively, the dominant spatial patterns (figures omitted) and PCs (Figs. 7a, b; 9a, b) can be obtained. For convenience, the spatial pattern is represented by the regression of the 200-hPa zonal wind anomalies onto normalized PCs, as illustrated in Figs. 7c-h; 9c-h. Then, the circulations associated with the changes in position and strength of AWJ derived from the three reanalyses are further investigated, as presented in Figs. 8,10. In winter, EOF1 is characterized by a tripole distribution, and explains about 33% of the total variance (Figs.  7a, c, e). The two strongest zonal wind anomaly centers appear over East Asia, and two relatively weaker centers reside over West Asia and to the south and north of the jet axis (30°N). EOF2, accounting for 15%-16% of the total variance, is a dipole pattern with relative weaker magnitude than EOF1 (Figs. 7b, d, f). The zonal wind anomalies are located over the western Pacific and Iranian Plateau. In terms of the circulations related to PC1 (Figs. 8a, c, e), East Asia exhibits a dipole pattern with a positive center located over Mongolia and Northeast China and a negative center over South China. Correspondingly, an anomalous anticyclonic-cyclonic circulation is displayed in the meridional direction. On the other hand, West Asia is occupied by negative values accompanied  In (a, b), the lines indicate the long-term linear trend of the PCs, and p and t are the probability and trend of lines, respectively. In (c, e, g), the red rectangles outline the key areas of change.

56
Journal of Meteorological Research Volume 35 by anomalous cyclonic circulations. When it turns to PC2 (Figs. 8b, d, f), the majority of Eurasia is occupied by negative correlation coefficients and anomalous cyclonic circulations.
In summer, EOF1 is characterized by a meridional dipole structure with zonal wind anomalies located over Northeast Asia and central Asia, accounting for about 24% of the total variance (Figs. 9a, c, e). EOF2 is dominated by a quarter-pole distribution with negative centers appearing over the Iranian and Mongolian plateaus and positive anomalies over central Asia and East China (Figs. 9b, d, f). Then, the 200-hPa circulations associated with PC1 are given in Figs. 10a, c, e. Significant positive correlations are confined to the middle to high latitudes (30°-50°N). The easterlies prevail over 30°-40°N while the westerlies over the high latitude. In terms of PC2, the 200-hPa circulation anomaly exhibits a Rossby wave anomaly, with anticyclonic circulations over the Mongolian Plateau and cyclonic circulations over central Asia. Eurasia is occupied by a tripole pattern in the zonal direction, with a negative center over central Asia and two positive centers over western Siberia and northern China (Figs. 10b, d, f).
Overall, the EOF decomposition results derived from the three reanalyses are highly consistent with each other in both winter and summer. Furthermore, we calculate statistical results for the EOF reanalysis in four seasons, respectively (figures omitted). The three reanalyses products agree extremely with each other. Differences for variances of EOF1 and EOF2 range almost within 1%. The EOF modes and corresponding PCs are also highly similar to each other, with their correlation coefficients excelling 0.97. In addition, the 200-hPa circulations associated with the changes in position and strength of AWJ derived from the three reanalyses are shown with very tiny discrepancies.

Long-term change of AWJ
In addition to interannual variability, the long-term changes of AWJ are also examined based on the PCs of EOF results. In winter (Figs. 7a, b), PC1s from the three reanalyses all show a statistically significant (p ≤ 0.05) decreasing tendency with a trend of −0.027 during 1979-2018, while no long-term trend is observed in PC2s. In summer (Figs. 9a, b), PC1s decrease slightly (t = −0.012 to −0.021) but are statistically insignificant  (p = 0.58-0.80). The PC2s of summertime AWJ from the three reanalyses all show a remarkable decreasing tendency (t = −0.031 to −0.032) statistically significant at the 5% level.
The significant decreasing trends of PC1 in winter and PC2 in summer indicate that the AWJ has changed in the past 40 years. To illustrate the detailed changes of AWJ over the key areas (Figs. 7c,e,g;9d,f,h), we further show the latitude-time cross-sections of 200-hPa zonal wind anomalies and jet axis from the three reanalyses (Fig. 11). Overall, the evolution in 200-hPa zonal wind anomalies from the three reanalyses resemble each other considerably. Interdedacal changes in the zonal wind anomalies over East Asia are obvious, which occur in the late 1990s in winter (Figs. 11a-c) and the mid 1990s in summer . In winter, negative anomalies over East Asia are shown to the north of the AWJ axis before the late 1990s, and turns to positive afterwards. In summer, the positive anomalies over East Asia appear before the mid 1990s, and turn to negative afterwards. For the westerly jet over central Asia in summer (Figs. 11d-f), the interdecadal changes are indistinctive.
To further identify the interdecadal changes of the AWJ axis, the 200-hPa mean zonal wind differences between 2000-2018 (1996-2018) and 1979-1999 (1979-1995) in winter (summer) from the three reanalyses are illustrated in Fig. 12. All of the three reanalyses describe considerably consistent distributions of the 200-hPa mean zonal wind differences, which are characterized by a meridional dipole pattern along the jet axis over

58
Journal of Meteorological Research Volume 35 Fig. 10. As in Fig. 8, but for summer.  Yu, X. J., L. X. Zhang,T. J. Zhou,et al. 59 East Asia. In winter, the zonal winds enhance to the north of the jet axis while weaken to the south of the jet axis, which means that the wintertime jet axis tends to move northward slightly during the past two decades. In summer, the zonal wind change is characterized by a quarter-pole distribution with positive anomalies over the Iranian and Mongolian plateaus and negative anomalies over East Asia and central Asia. As a result, the jet axis migrates southward over central Asia (40°-80°E) while northward over East Asia (80°-110°E) during the past two decades. Across the three reanalyses, the magnitude of zonal wind enhancement from ERA5 (4 m s −1 ) is weaker than CRA-40 and CFSR (5 m s −1 ), especially over the Iranian and Tibetan plateaus regions in winter.

Conclusions and discussion
In this study, we made detailed comparison between CRA-40, ERA5, and CFSR during the period of 1979-2018 with respect to the climatic structure, seasonal cycle, interannual variability, and long-term change of AWJ, aiming at evaluating the performance of CMA reanalysis data in describing AWJ. The results are outlined as follows.
(1) The three reanalyses show comparable spatial patterns and quantity both in the horizontal and vertical structures of AWJ, with their pattern correlation coefficients excelling 0.99 and regional average difference values less than 0.4 m s −1 . Differences (relative differences) of 200-hPa zonal wind range within ± 0.5 m s −1 (± 2%) over a majority of regions, with some distinct differences in winter than in summer. Larger differences with a maximum of ± 2 m s −1 (± 5%) appear over the Iranian and Tibetan plateaus regions (50°-90°E) in the midupper troposphere in winter.
(2) The three datasets are also highly consistent with each other in the seasonal cycle in terms of the latitude, longitude, and intensity of the AWJ centers (with correlation coefficients ≥ 0.98). Some obvious discrepancies appear in MAM and SON, especially for the western center (30°-70°E). The correlation coefficients between zonal shift identified in ERA5 and CRA-40 (CFSR) are less than 0.80.
(4) The dominant spatial patterns and corresponding PCs obtained by the EOF decomposition across the three reanalyses strongly resemble each other with their correlation coefficients exceeding 0.97. By defining the PC1 (PC2) as the AWJ position (strength) index, the related circulations at 200 hPa derived from the three reanalyses are shown with some very tiny discrepancies.
(5) The significant decreasing trends of PC1s in winter and PC2s in summer indicate that the AWJ have evident changes during 1979-2018. In winter, the AWJ moves northward with positive zonal wind anomaly over most parts of China. In summer, the westerly jet undergoes a southward displacement over central Asia (40°-80°E) and a northward migration over East Asia (80°-110°E). It is noticeable that the magnitude of zonal wind enhancement to the north of the AWJ axis is weaker from ERA5 than from CRA-40 and CFSR, especially over the Iranian and Tibetan plateaus in winter.
This study has systematically verified the considerable agreement among the three reanalyses on the above perspectives. Extended discussion is as follows. 1) Due to more observations over East Asia that have been assimilated in the CRA system (Liao et al., 2018;Yin et al., 2018), it is presumable that CRA-40 is more reliable in depicting the atmospheric circulations over East Asia. Nonetheless, some discrepancies in describing the characteristics of EAWJ and its associated atmospheric circulations across the three reanalyses indeed exist. 2) In terms of absolute values, the biggest differences across the three reanalyses appear in winter, which is mainly because AWJ is much stronger in winter than in summer. 3) Regionally, the biggest differences across the three reanalyses lie in the Iranian and Tibetan plateaus. This may be attributable to the complex topography, which has evident dynamic and thermal effects on the troposphere (Yeh, 1982;Xia et al., 2016;Jiang et al., 2020) and depends highly on models (Gao et al., 2017;Li et al., 2018) and data assimilation (Seto et al., 2013;Yu et al., 2018;He et al., 2019). In summary, this study has demonstrated the rationality and reliability of CRA-40 in characterizing the AWJ, and provided support to use of CRA-40 in analyzing the large-scale atmospheric circulations in East Asia in future.